skip to main content

Title: A novel stochastic optimization method for handling misalignments of proton and photon doses in combined treatments

Objective.Combined proton–photon treatments, where most fractions are delivered with photons and only a few are delivered with protons, may represent a practical approach to optimally use limited proton resources. It has been shown that, when organs at risk (OARs) are located within or near the tumor, the optimal multi-modality treatment uses protons to hypofractionate parts of the target volume and photons to achieve near-uniform fractionation in dose-limiting healthy tissues, thus exploiting the fractionation effect. These plans may be sensitive to range and setup errors, especially misalignments between proton and photon doses. Thus, we developed a novel stochastic optimization method to directly incorporate these uncertainties into the biologically effective dose (BED)-based simultaneous optimization of proton and photon plans.Approach.The method considers the expected valueEband standard deviationσbof the cumulative BEDbin every voxel of a structure. For the target, a piecewise quadratic penalty function of the formbminEb2σb+2is minimized, aiming for plans in which the expected BED minus two times the standard deviation exceeds the prescribed BEDbmin.Analogously,Eb+2σbbmax+2is considered for OARs.Main results.Using a spinal metastasis case and a liver cancer patient, it is demonstrated that the novel stochastic optimization method yields robust combined treatment plans. Tumor coverage and a good sparing of the main OARs are maintained despite range and setup errors, and especially misalignments between proton and photon doses. This is achieved without explicitly considering all combinations of proton and photon error scenarios.Significance.Concerns about range and setup errors for safe clinical implementation of optimized proton–photon radiotherapy can be addressed through an appropriate stochastic planning method.

more » « less
Award ID(s):
1847865 1719828
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
IOP Publishing
Date Published:
Journal Name:
Physics in Medicine & Biology
Page Range / eLocation ID:
Article No. 185006
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We present a stellar dynamical mass measurement of a newly detected supermassive black hole (SMBH) at the center of the fast-rotating, massive elliptical galaxy NGC 2693 as part of the MASSIVE survey. We combine high signal-to-noise ratio integral field spectroscopy (IFS) from the Gemini Multi-Object Spectrograph with wide-field data from the Mitchell Spectrograph at McDonald Observatory to extract and model stellar kinematics of NGC 2693 from the central ∼150 pc out to ∼2.5 effective radii. Observations from Hubble Space Telescope WFC3 are used to determine the stellar light distribution. We perform fully triaxial Schwarzschild orbit modeling using the latest TriOS code and a Bayesian search in 6D galaxy model parameter space to determine NGC 2693's SMBH mass (MBH), stellar mass-to-light ratio, dark matter content, and intrinsic shape. We findMBH=1.7±0.4×109Mand a triaxial intrinsic shape with axis ratiosp=b/a= 0.902 ± 0.009 andq=c/a=0.7210.010+0.011, triaxiality parameterT= 0.39 ± 0.04. In comparison, the best-fit orbit model in the axisymmetric limit and (cylindrical) Jeans anisotropic model of NGC 2693 preferMBH=2.4±0.6×109MandMBH=2.9±0.3×109M, respectively. Neither model can account for the non-axisymmetric stellar velocity features present in the IFS data.

    more » « less
  2. Abstract

    We present a detection of 21 cm emission from large-scale structure (LSS) between redshift 0.78 and 1.43 made with the Canadian Hydrogen Intensity Mapping Experiment. Radio observations acquired over 102 nights are used to construct maps that are foreground filtered and stacked on the angular and spectral locations of luminous red galaxies (LRGs), emission-line galaxies (ELGs), and quasars (QSOs) from the eBOSS clustering catalogs. We find decisive evidence for a detection when stacking on all three tracers of LSS, with the logarithm of the Bayes factor equal to 18.9 (LRG), 10.8 (ELG), and 56.3 (QSO). An alternative frequentist interpretation, based on the likelihood ratio test, yields a detection significance of 7.1σ(LRG), 5.7σ(ELG), and 11.1σ(QSO). These are the first 21 cm intensity mapping measurements made with an interferometer. We constrain the effective clustering amplitude of neutral hydrogen (Hi), defined asHI103ΩHIbHI+fμ2, where ΩHiis the cosmic abundance of Hi,bHiis the linear bias of Hi, and 〈fμ2〉 = 0.552 encodes the effect of redshift-space distortions at linear order. We findHI=1.510.97+3.60for LRGs (z= 0.84),HI=6.763.79+9.04for ELGs (z= 0.96), andHI=1.680.67+1.10for QSOs (z= 1.20), with constraints limited by modeling uncertainties at nonlinear scales. We are also sensitive to bias in the spectroscopic redshifts of each tracer, and we find a nonzero bias Δv= − 66 ± 20 km s−1for the QSOs. We split the QSO catalog into three redshift bins and have a decisive detection in each, with the upper bin atz= 1.30 producing the highest-redshift 21 cm intensity mapping measurement thus far.

    more » « less
  3. Abstract

    We conduct a systematic tidal disruption event (TDE) demographics analysis using the largest sample of optically selected TDEs. A flux-limited, spectroscopically complete sample of 33 TDEs is constructed using the Zwicky Transient Facility over 3 yr (from 2018 October to 2021 September). We infer the black hole (BH) mass (MBH) with host galaxy scaling relations, showing that the sampleMBHranges from 105.1Mto 108.2M. We developed a survey efficiency corrected maximum volume method to infer the rates. The rest-frameg-band luminosity function can be well described by a broken power law ofϕ(Lg)Lg/Lbk0.3+Lg/Lbk2.61, withLbk= 1043.1erg s−1. In the BH mass regime of 105.3≲ (MBH/M) ≲ 107.3, the TDE mass function followsϕ(MBH)MBH0.25, which favors a flat local BH mass function (dnBH/dlogMBHconstant). We confirm the significant rate suppression at the high-mass end (MBH≳ 107.5M), which is consistent with theoretical predictions considering direct capture of hydrogen-burning stars by the event horizon. At a host galaxy mass ofMgal∼ 1010M, the average optical TDE rate is ≈3.2 × 10−5galaxy−1yr−1. We constrain the optical TDE rate to be [3.7, 7.4, and 1.6] × 10−5galaxy−1yr−1in galaxies with red, green, and blue colors.

    more » « less
  4. Abstract

    The genericity of Arnold diffusion in the analytic category is an open problem. In this paper, we study this problem in the followinga prioriunstable Hamiltonian system with a time-periodic perturbationHε(p,q,I,φ,t)=h(I)+i=1n±12pi2+Vi(qi)+εH1(p,q,I,φ,t),where(p,q)Rn×Tn,(I,φ)Rd×Tdwithn,d⩾ 1,Viare Morse potentials, andɛis a small non-zero parameter. The unperturbed Hamiltonian is not necessarily convex, and the induced inner dynamics does not need to satisfy a twist condition. Using geometric methods we prove that Arnold diffusion occurs for generic analytic perturbationsH1. Indeed, the set of admissibleH1isCωdense andC3open (a fortiori,Cωopen). Our perturbative technique for the genericity is valid in theCktopology for allk∈ [3, ∞) ∪ {∞,ω}.

    more » « less
  5. Abstract

    We present new measurements of cosmic microwave background (CMB) lensing over 9400 deg2of the sky. These lensing measurements are derived from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) CMB data set, which consists of five seasons of ACT CMB temperature and polarization observations. We determine the amplitude of the CMB lensing power spectrum at 2.3% precision (43σsignificance) using a novel pipeline that minimizes sensitivity to foregrounds and to noise properties. To ensure that our results are robust, we analyze an extensive set of null tests, consistency tests, and systematic error estimates and employ a blinded analysis framework. Our CMB lensing power spectrum measurement provides constraints on the amplitude of cosmic structure that do not depend on Planck or galaxy survey data, thus giving independent information about large-scale structure growth and potential tensions in structure measurements. The baseline spectrum is well fit by a lensing amplitude ofAlens= 1.013 ± 0.023 relative to the Planck 2018 CMB power spectra best-fit ΛCDM model andAlens= 1.005 ± 0.023 relative to the ACT DR4 + WMAP best-fit model. From our lensing power spectrum measurement, we derive constraints on the parameter combinationS8CMBLσ8Ωm/0.30.25ofS8CMBL=0.818±0.022from ACT DR6 CMB lensing alone andS8CMBL=0.813±0.018when combining ACT DR6 and PlanckNPIPECMB lensing power spectra. These results are in excellent agreement with ΛCDM model constraints from Planck or ACT DR4 + WMAP CMB power spectrum measurements. Our lensing measurements from redshiftsz∼ 0.5–5 are thus fully consistent with ΛCDM structure growth predictions based on CMB anisotropies probing primarilyz∼ 1100. We find no evidence for a suppression of the amplitude of cosmic structure at low redshifts.

    more » « less