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

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 value$Eb$and standard deviation$σb$of the cumulative BED$b$in every voxel of a structure. For the target, a piecewise quadratic penalty function of the form$bmin−Eb−2σb+2$is minimized, aiming for plans in which the expected BED minus two times the standard deviation exceeds the prescribed BED$bmin.$Analogously,is 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.

Authors:
; ; ;
Publication Date:
NSF-PAR ID:
10371467
Journal Name:
Physics in Medicine & Biology
Volume:
67
Issue:
18
Page Range or eLocation-ID:
Article No. 185006
ISSN:
0031-9155
Publisher:
IOP Publishing
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 find$MBH=1.7±0.4×109M⊙$and a triaxial intrinsic shape with axis ratiosp=b/a= 0.902 ± 0.009 and$q=c/a=0.721−0.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 prefer$MBH=2.4±0.6×109M⊙$and$MBH=2.9±0.3×109M⊙$, respectively. Neither model can account for the non-axisymmetric stellar velocity features present inmore »

2. 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 perturbation$Hε(p,q,I,φ,t)=h(I)+∑i=1n±12pi2+Vi(qi)+εH1(p,q,I,φ,t),$where$(p,q)∈Rn×Tn$,$(I,φ)∈Rd×Td$withn,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, ∞) ∪ {∞,ω}.

3. Abstract

We perform particle-in-cell simulations to elucidate the microphysics of relativistic weakly magnetized shocks loaded with electron-positron pairs. Various external magnetizationsσ≲ 10−4and pair-loading factorsZ±≲ 10 are studied, whereZ±is the number of loaded electrons and positrons per ion. We find the following: (1) The shock becomes mediated by the ion Larmor gyration in the mean field whenσexceeds a critical valueσLthat decreases withZ±. AtσσLthe shock is mediated by particle scattering in the self-generated microturbulent fields, the strength and scale of which decrease withZ±, leading to lowerσL. (2) The energy fraction carried by the post-shock pairs is robustly in the range between 20% and 50% of the upstream ion energy. The mean energy per post-shock electron scales as$E¯e∝Z±+1−1$. (3) Pair loading suppresses nonthermal ion acceleration at magnetizations as low asσ≈ 5 × 10−6. The ions then become essentially thermal with mean energy$E¯i$, while electrons form a nonthermal tail, extending from$E∼Z±+1−1E¯i$to$E¯i$. Whenσ= 0, particle acceleration is enhanced by the formation of intense magnetic cavities that populate the precursor during the late stages of shock evolution. Here,more »

4. Abstract

Cosmic reionization was the last major phase transition of hydrogen from neutral to highly ionized in the intergalactic medium (IGM). Current observations show that the IGM is significantly neutral atz> 7 and largely ionized byz∼ 5.5. However, most methods to measure the IGM neutral fraction are highly model dependent and are limited to when the volume-averaged neutral fraction of the IGM is either relatively low ($x¯HI≲10−3$) or close to unity ($x¯HI∼1$). In particular, the neutral fraction evolution of the IGM at the critical redshift range ofz= 6–7 is poorly constrained. We present new constraints on$x¯HI$atz∼ 5.1–6.8 by analyzing deep optical spectra of 53 quasars at 5.73 <z< 7.09. We derive model-independent upper limits on the neutral hydrogen fraction based on the fraction of “dark” pixels identified in the Lyαand Lyβforests, without any assumptions on the IGM model or the intrinsic shape of the quasar continuum. They are the first model-independent constraints on the IGM neutral hydrogen fraction atz∼ 6.2–6.8 using quasar absorption measurements. Our results give upper limits of(1σ),$x¯HI(z=6.5)<0.87±0.03$(1σ), and$x¯HI(z=6.7)<0.94−0.09+0.06$(1σ). The dark pixel fractions atz> 6.1 are consistent with the redshift evolution of the neutral fraction of the IGM derived from Planck 2018.

5. Abstract

The merger of two galaxies, each hosting a supermassive black hole (SMBH) of mass 106Mor more, could yield a bound SMBH binary. For the early-type galaxy NGC 4472, we study how astrometry with a next-generation Very Large Array could be used to monitor the reflex motion of the primary SMBH of massMpri, as it is tugged on by the secondary SMBH of mass$Msec$. Casting the orbit of the putative SMBH binary in terms of its periodP, semimajor axisabin, and mass ratio$q=Msec/Mpri≤1$, we find the following: (1) Orbits with fiducial periods ofP= 4 yr and 40 yr could be spatially resolved and monitored. (2) For a 95% accuracy of 2μas per monitoring epoch, subparsec values ofabincould be accessed over a range of mass ratios notionally encompassing major$q>14$and minor$q<14$galaxy mergers. (3) If no reflex motion is detected forMpriafter 1 (10) yr of monitoring, an SMBH binary with periodP= 4 (40) yr and mass ratioq> 0.01 (0.003) could be excluded. This would suggest no present-day evidence for a past major merger like that recently simulated, where scouring by aq∼ 1 SMBH binary formedmore »