Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “eclipse”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

HT Lyn and IR Lyn: Two Semi-detached-type Near-contact Binaries with Stable Orbital Period

Abstract We presented the first photometric and orbital period investigations for two near-contact binaries: HT Lyn and IR Lyn. The light-curves solutions derived from both our ground-based and various surveys’ observations suggested that HT Lyn and IR Lyn are two semi-detached-type near-contact binaries with the secondary and primary components filling the Roche lobe. Combining the eclipse timings derived from several surveys’ data and our observations with those reported in the literature, we modified the orbital periods and revealed that their orbital periods were stable for the last two decades. The absolute physical parameters were well determined by using the empirical relations and the Gaia-distance-based method. Similar to other near-contact binaries, the primary components of HT Lyn and IR Lyn are evolving at the main-sequence stage, while their less-massive secondary components show the over-sized and over-luminosity features and should be more evolved than their primary ones. Combining the stable orbital period with the semi-detached configurations, we infer that the two near-contact binaries may be just located in the critical phase and evolving from the current semi-detached phase to the detached or contact one.

Astronomy & Astrophysics↗
DOE OSTI2022

Assessing the performance of global thermostat adjustment in commercial buildings for load shifting demand response

Abstract Efficiently leveraging new sources of flexibility is critical to mitigating load balancing challenges posed by variable renewable resources. The thermal inertia of commercial buildings allows us to shift their power consumption on minute to hourly timescales to provide demand response to the grid while maintaining occupant comfort. Global thermostat adjustment (GTA) provides a readily available and scalable approach for implementing load shifting demand response using commercial heating, ventilation, and air conditioning (HVAC) systems, since it leverages the inherent sophistication of modern building automation systems. However, there is an incomplete understanding of GTA’s performance for this purpose and its impact on building systems and occupant comfort. In this paper, we explore the performance of GTA by analyzing results from nearly nine hundred experiments on eight university campus buildings in Michigan and North Carolina. Using GTA, we manipulate each building’s thermostat setpoints causing the building to shift its power consumption with respect to its baseline. We quantify the magnitude of HVAC power response, energy use of HVAC subsystems, and impact on occupant comfort. Finally, we connect our experimental results with power system operation using an optimization model that coordinates GTA actions across a large collection of grid-interactive efficient buildings (GEBs) to reduce high ramp rates on the grid and mitigate renewable energy curtailment. Overall, our work finds that the impacts on HVAC subsystems are often complex, and may result in additional energy being consumed by fans and terminal reheat. These effects must be considered when using GTA for load shifting. Additionally, we demonstrate that occupant comfort, as assessed by indoor temperature and humidity, can be maintained during GTA events. From a societal perspective, our modeling work finds that the additional renewable energy that can be integrated through the use of GTA strategies eclipses any additional energy consumed by buildings.

Keskar, Aditya (ORCID:0000000244617980)↗
DOE OSTI2021

High tide: a systematic search for ellipsoidal variables in ASAS-SN

ABSTRACT The majority of non-merging stellar mass black holes are discovered by observing high energy emission from accretion processes. Here, we pursue the large, but still mostly unstudied population of non-interacting black holes and neutron stars by searching for the tidally induced ellipsoidal variability of their stellar companions. We start from a sample of about 200 000 rotational variables, semiregular variables, and eclipsing binary stars from the All-Sky Automated Survey for Supernovae. We use a χ2 ratio test followed by visual inspection to identify 369 candidates for ellipsoidal variability. We also discuss how to combine the amplitude of the variability with mass and radius estimates for observed stars to calculate a minimum companion mass, identifying the most promising candidates for high mass companions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

BEBOP II: sensitivity to sub-Saturn circumbinary planets using radial-velocities

ABSTRACT BEBOP is a radial-velocity survey that monitors a sample of single-lined eclipsing binaries, in search of circumbinary planets by using high-resolution spectrographs. Here, we describe and test the methods we use to identify planetary signals within the BEBOP data and establish how we quantify our sensitivity to circumbinary planets by producing detection limits. This process is made easier and more robust by using a diffusive nested sampler. In the process of testing our methods, we notice that contrary to popular wisdom, assuming circular orbits in calculating detection limits for a radial-velocity survey provides overoptimistic detection limits by up to $40{{\ \rm per\ cent}}$ in semi-amplitude with implications for all radial-velocity surveys. We perform example analyses using three BEBOP targets from our Southern HARPS survey. We demonstrate for the first time a repeated ability to reach a residual root mean squared scatter of $3\, \rm m\, s^{-1}$ (after removing the binary signal), and find that we are sensitive to circumbinary planets with masses down to that of Neptune and Saturn, for orbital periods up to $1000\, \rm d$.

Standing, Matthew R. (ORCID:0000000276088905)↗
DOE OSTI2022

Trumpler 16-26: a new centrifugal magnetosphere star discovered via SDSS/APOGEE H -band spectroscopy

ABSTRACT We report the discovery of a new example of the rare class of highly magnetized, rapidly rotating, helium enhanced, early B stars that produce anomalously wide hydrogen emission due to a centrifugal magnetosphere (CM). The star is Trumpler 16-26, a B1.5 V member of the Trumpler 16 open cluster. A CM was initially suspected based on hydrogen Brackett series emission observed in SDSS/APOGEE H-band spectra. Similar to the other stars of this type, the emission was highly variable and at all times remarkable due to the extreme velocity separations of the double peaks (up to 1300 km s−1.) Another clue lay in the TESS light curve, which shows two irregular eclipses per cycle when phased with the likely 0.971 8115-d rotation period, similar to the behaviour of the well-known CM host star σ Ori E. To confirm a strong magnetic field and rotation-phase-locked variability, we initiated a follow-up campaign consisting of optical spectropolarimetry and spectroscopy. The associated data revealed a longitudinal magnetic field varying between −3.1 and +1.6 kG with the period found from photometry. The optical spectra confirmed rapid rotation (v sin i = 195 km s−1), surface helium enhancement, and wide, variable hydrogen emission. Tr16-26 is thus confirmed as the 20th known, the fourth most rapidly rotating, and the faintest CM host star yet discovered. With a projected dipole magnetic field strength of Bd > 11 kG, Tr16-26 is also among the most magnetic CM stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The TRAPUM L -band survey for pulsars in Fermi -LAT gamma-ray sources

ABSTRACT More than 100 millisecond pulsars (MSPs) have been discovered in radio observations of gamma-ray sources detected by the Fermi Large Area Telescope (LAT), but hundreds of pulsar-like sources remain unidentified. Here, we present the first results from the targeted survey of Fermi-LAT sources being performed by the Transients and Pulsars with MeerKAT (TRAPUM) Large Survey Project. We observed 79 sources identified as possible gamma-ray pulsar candidates by a Random Forest classification of unassociated sources from the 4FGL catalogue. Each source was observed for 10 min on two separate epochs using MeerKAT’s L-band receiver (856–1712 MHz), with typical pulsed flux density sensitivities of $\sim 100\, \mu$Jy. Nine new MSPs were discovered, eight of which are in binary systems, including two eclipsing redbacks and one system, PSR J1526−2744, that appears to have a white dwarf companion in an unusually compact 5 h orbit. We obtained phase-connected timing solutions for two of these MSPs, enabling the detection of gamma-ray pulsations in the Fermi-LAT data. A follow-up search for continuous gravitational waves from PSR J1526−2744 in Advanced LIGO data using the resulting Fermi-LAT timing ephemeris yielded no detection, but sets an upper limit on the neutron star ellipticity of 2.45 × 10−8. We also detected X-ray emission from the redback PSR J1803−6707 in data from the first eROSITA all-sky survey, likely due to emission from an intrabinary shock.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A prescription for the asteroseismic surface correction

In asteroseismology, the surface effect refers to a disparity between the observed and the modelled frequencies in stars with solar-like oscillations. It originates from improper modelling of the surface layers. Correcting the surface effect usually requires using functions with free parameters, which are conventionally fitted to the observed frequencies. On the basis that the correction should vary smoothly across the H–R diagram, we parameterize it as a simple function of surface gravity, effective temperature, and metallicity. We determine this function by fitting a wide range of stars. The absolute amount of the surface correction decreases with luminosity, but the ratio between it and νmax increases, suggesting the surface effect is more important for red giants than dwarfs. Applying the prescription can eliminate unrealistic surface correction, which improves parameter estimations with stellar modelling. Using two open clusters, we found a reduction of scatter in the model-derived ages for each star in the same cluster. As an important application, we provide a new revision for the Δν scaling relation that, for the first time, accounts for the surface correction. The values of the correction factor, f Δν , are up to 2 percent smaller than those determined without the surface effect considered, suggesting decreases of up to 4 percent in radii and up to 8 percent in masses when using the asteroseismic scaling relations. This revision brings the asteroseismic properties into an agreement with those determined from eclipsing binaries. The new correction factor and the stellar models with the corrected frequencies are available at https://www.github.com/parallelpro/surface.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Shedding far-ultraviolet light on the donor star and evolutionary state of the neutron-star LMXB Swift J1858.6−0814

ABSTRACT The evolution of accreting X-ray binary systems is closely coupled to the properties of their donor stars. Consequently, we can constrain the evolutionary track a system is by establishing the nature of its donor. Here, we present far-ultraviolet (far-UV) spectroscopy of the transient neutron-star low-mass X-ray binary J1858 in different accretion states (low-hard, high-hard, and soft). All of these spectra exhibit anomalous N v, C iv, Si iv, and He ii lines, suggesting that its donor star has undergone CNO processing. We also determine the donor’s effective temperature, Td ≃ 5700 K, and radius, Rd ≃ 1.7 R⊙, based on photometric observations obtained during quiescence. Lastly, we leverage the transient nature of the system to set an upper limit of $\dot{M}_{\rm acc} \lesssim 10^{-8.5}~{\rm M}_{\odot }~\mathrm{ yr}^{-1}$ on the present-day mass-transfer rate. Combining these with the orbital period of the system, Porb = 21.3 h, we search for viable evolution paths. The initial donor masses in the allowed solutions span the range 1 M⊙ ≲ Md,i ≲ 3.5 M⊙. All but the lowest masses in this range are consistent with the strong CNO-processing signature in the UV line ratios. The present-day donor mass in the permitted tracks are 0.5 M⊙ ≲ Md,obs ≲ 1.3 M⊙, higher than suggested by eclipse modelling. Since Porb is close to the so-called bifurcation period, both converging and diverging binary tracks are permitted. If the former is confirmed, J1858 will end its life as an ultracompact system with a substellar donor.

Castro Segura, N. (ORCID:0000000258700443)↗
DOE OSTI2024

Short-period post-common envelope binaries with Balmer emission from SDSS and LAMOST based on ZTF photometric data

ABSTRACT We present here 55 short-period post-common envelope binaries (PCEBs) containing a hot white dwarf (WD) and a low-mass main sequence (MS). Based on the photometric data from Zwicky Transient Facility survey data Release 19 (ZTF DR19), the light curves are analysed for about 200 WDMS binaries with emission line(s) identified from the Sloan Digital Sky Survey (SDSS) or the Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST) spectra, in which 55 WDMS binaries are found to exhibit variability in their luminosities with a short period and are thus short-period binaries (i.e. PCEBs). In addition, it is found that the orbital periods of these PCEBs locate in a range from 2.2643 to 81.1526 h. However, only six short-period PCEBs are newly discovered and the orbital periods of 19 PCEBs are improved in this work. Meanwhile, it is found that three objects are newly discovered eclipsing PCEBs, and a object (i.e. SDSS J1541) might be the short-period PCEB with a late M-type star or a brown dwarf companion based on the analysis of its spectral energy distribution. At last, the mechanism(s) being responsible for the emission features in the spectra of these PCEBs are discussed, the emission features arising in their optical spectra might be caused by the stellar activity or an irradiated component owing to a hot WD companion because most of them contain a WD with an effective temperature higher than $\sim$10 000 K.

Li, Lifang↗
DOE OSTI2019

Within-host infectious disease models accommodating cellular coinfection, with an application to influenza

Within-host models are useful tools for understanding the processes regulating viral load dynamics. While existing models have considered a wide range of within-host processes, at their core these models have shown remarkable structural similarity. Specifically, the structure of these models generally consider target cells to be either uninfected or infected, with the possibility of accommodating further resolution (e.g. cells that are in an eclipse phase). Recent findings, however, indicate that cellular coinfection is the norm rather than the exception for many viral infectious diseases, and that cells with high multiplicity of infection are present over at least some duration of an infection. The reality of these cellular coinfection dynamics is not accommodated in current within-host models although it may be critical for understanding within-host dynamics. This is particularly the case if multiplicity of infection impacts infected cell phenotypes such as their death rate and their viral production rates. Here, we present a new class of within-host disease models that allow for cellular coinfection in a scalable manner by retaining the low-dimensionality that is a desirable feature of many current within-host models. The models we propose adopt the general structure of epidemiological ‘macroparasite’ models that allow hosts to be variably infected by parasites such as nematodes and host phenotypes to flexibly depend on parasite burden. Specifically, our within-host models consider target cells as ‘hosts’ and viral particles as ‘macroparasites’, and allow viral output and infected cell lifespans, among other phenotypes, to depend on a cell’s multiplicity of infection. We show with an application to influenza that these models can be statistically fit to viral load and other within-host data, and demonstrate using model selection approaches that they have the ability to outperform traditional within-host viral dynamic models. Important in vivo quantities such as the mean multiplicity of cellular infection and time-evolving reassortant frequencies can also be quantified in a straightforward manner once these macroparasite models have been parameterized. The within-host model structure we develop here provides a mathematical way forward to address questions related to the roles of cellular coinfection, collective viral interactions, and viral complementation in within-host viral dynamics and evolution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Five percent measurement of the gravitational constant in the Large Magellanic Cloud

We perform a novel test of General Relativity by measuring the gravitational constant in the Large Magellanic Cloud (LMC). The LMC contains six well-studied Cepheid variable stars in detached eclipsing binaries. Radial velocity and photometric observations enable a complete orbital solution, and precise measurements of the Cepheids’ periods permit detailed stellar modelling. Both are sensitive to the strength of gravity, the former via Kepler’s third law and the latter through the gravitational free-fall time. We jointly fit the observables for stellar parameters and the gravitational constant. Performing a full Markov Chain Monte Carlo analysis of the parameter space including all relevant nuisance parameters, we constrain the gravitational constant in the Large Magellanic Cloud relative to the Solar System to be G LMC /G SS = $0.93{^{+0.05}_{-0.04}}$. We discuss the implications of this 5% measurement of Newton’s constant in another galaxy for dark energy and modified gravity theories. This result excludes one Cepheid, CEP-1812, which is an outlier and needs further study: it is either a highly unusual system to which our model does not apply, or it prefers G LMC < G SS at 2.6σ. We also obtain new bounds on critical parameters that appear in semi-analytic descriptions of stellar processes. In particular, we measure the mixing length parameter to be $α = 0.90{^{+0.36}_{-0.26}}$ (when assumed to be constant across our sample), and obtain constraints on the parameters describing turbulent dissipation and convective flux.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Probing the Melting of a Two-Dimensional Quantum Wigner Crystal via its Screening Efficiency

One of the most fundamental and yet elusive collective phases of an interacting electron system is the quantum Wigner crystal (WC), an ordered array of electrons expected to form when the electrons’ Coulomb repulsion energy eclipses their kinetic (Fermi) energy. In low-disorder, two-dimensional (2D) electron systems, the quantum WC is known to be favored at very low temperatures ($T$) and small Landau level filling factors ($ν$), near the termination of the fractional quantum Hall states. This WC phase exhibits an insulating behavior, reflecting its pinning by the small but finite disorder potential. An experimental determination of a $T$ vs $ν$ phase diagram for the melting of the WC, however, has proved to be challenging. Here we use capacitance measurements to probe the 2D WC through its effective screening as a function of $T$ and $ν$. We find that, as expected, the screening efficiency of the pinned WC is very poor at very low $T$ and improves at higher $T$ once the WC melts. Surprisingly, however, rather than monotonically changing with increasing $T$, the screening efficiency shows a well-defined maximum at a $T$ that is close to the previously reported melting temperature of the WC. Our experimental results suggest a new method to map out a $T$ vs $ν$ phase diagram of the magnetic-field-induced WC precisely.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Detection of magnetic circular dichroism in amorphous materials utilizing a single-crystalline overlayer

Physicists are fascinated with topological defects in solid-state materials, because by breaking the translational symmetry they offer emerging properties that are not present in their parental phases. For example, edge dislocations—the 2π phase-winding topological defects—in antiferromagnetic NiO crystals can exhibit ferromagnetic behaviors. Herein, we study how these defects could give rise to exotic topological orders when they interact with a high energy electron beam. To probe this interaction, we formed a coherent electron nanobeam in a scanning transmission electron microscope and recorded the far-field transmitted patterns as the beam steps through the edge dislocation core in [001] NiO. Surprisingly, we found the amplitude patterns of the <020> Bragg disks evolve in a similar manner to the evolution of an annular solar eclipse. Using the ptychographic technique, we recovered the missing phase information in the diffraction plane and revealed the topological phase vortices in the diffracted beams. Through atomic topological defects, the wave function of electrons can be converted from plane wave to electron vortex. This approach provides a new perspective for boosting the collection efficiency of magnetic circular dichroism spectra with high spatial resolution and understanding the relationship between symmetry breaking and exotic property of individual topological defect at atomic level.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mapping hidden space-charge distributions across crystalline metal oxide/group IV semiconductor interfaces

The electronic structures of semiconducting heterojunctions are critically dependent on composition including the presence and concentrations of dopants, both intended and unintended. Dopant profiles in the interfacial region can have major effects on band energies which in turn drive transport properties. Here, in this study, we use core-level photoelectron line shapes excited with hard x rays to extract information about electric fields resulting from internal charge transfer in epitaxial La 0.03 Sr 0.97 Zr x Ti 1–x O 3 /Ge(001) (0.1≤x≤0.7) heterostructures. Experiments were carried out for heterojunctions involving both n- and p-type Ge substrates. These heterojunctions were not amenable to electronic characterization of all regions by transport measurements because the doped substrates act as electrical shunts, precluding probing the more resistive films and masking interface conductivity. However, the core-level line shapes were found to be a rich source of information on built-in potentials that exist throughout the heterostructure, and yielded valuable insight into the impact of band bending on band alignment at the buried interfaces. The electronic effects expected for Ge with uniform n- and p-type doping are eclipsed by those of unintended oxygen dopants in the Ge near the interface. This study illustrates the power of hard x-ray photoemission spectroscopy and related modeling to determine electronic structure in material systems for which insight from traditional transport measurements is limited.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Analyzing diffuse scattering with supercomputers. Corrigendum

The study by Michels-Clark et al. (2013 [Michels-Clark, T. M., Lynch, V. E., Hoffmann, C. M., Hauser, J., Weber, T., Harrison, R. & Bürgi, H. B. (2013). J. Appl. Cryst. 46, 1616-1625.]) contains misleading errors which are corrected here. The numerical results reported in that paper and the conclusions given there are not affected and remain unchanged. The transition probabilities in Table 1 (rows 4, 5, 7, 8) and Fig. 2 (rows 1 and 2) of the original paper were different from those used in the numerical calculations. Corrected transition probabilities as used in the computations are given in Table 1 and Fig. 1 of this article. The Δ parameter in the stacking model expresses the preference for the fifth layer in a five-layer stack to be eclipsed with respect to the first layer. This statement corrects the original text on p. 1622, lines 4–7. In the original Fig. 2 the helicity of the layer stacks b L and b R in rows 3 and 4 had been given as opposite to those in rows 1, 2 and 5. Fig. 1 of this article shows rows 3 and 4 corrected to correspond to rows 1, 2 and 5.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Commissioning of intensity modulated neutron radiotherapy (IMNRT)

Purpose: Intensity modulated neutron radiotherapy (IMNRT) has been developed using inhouse treatment planning and delivery systems at the Karmanos Cancer Center/Wayne State University Fast Neutron Therapy facility. The process of commissioning IMNRT for clinical use is presented here. Results of commissioning tests are provided including validation measurements using representative patient plans as well as those from the TG-119 test suite. Methods: IMNRT plans were created using the Varian Eclipse optimization algorithm and an inhouse planning system for calculation of neutron dose distributions. Tissue equivalent ionization chambers and an ionization chamber array were used for point dose and planar dose distribution comparisons with calculated values. Validation plans were delivered to water and virtual water phantoms using TG-119 measurement points and evaluation techniques. Photon and neutron doses were evaluated both inside and outside the target volume for a typical IMNRT plan to determine effects of intensity modulation on the photon dose component. Monitor unit linearity and effects of beam current and gantry angle on output were investigated, and an independent validation of neutron dosimetry was obtained. Results: While IMNRT plan quality is superior to conventional fast neutron therapy plans for clinical sites such as prostate and head and neck, it is inferior to photon IMRT for most TG-119 planning goals, particularly for complex cases. This results significantly from current limitations on the number of segments. Measured and calculated doses for 11 representative plans (six prostate/five head and neck) agreed to within -0.8 {+-} 1.4% and 5.0 {+-} 6.0% within and outside the target, respectively. Nearly all (22/24) ion chamber point measurements in the two phantom arrangements were within the respective confidence intervals for the quantity [(measured-planned)/prescription dose] derived in TG-119. Mean differences for all measurements were 0.5% (max = 7.0%) and 1.4% (max = 4.1%) in water and virtual water, respectively. The mean gamma pass rate for all cases was 92.8% (min = 88.6%). These pass rates are lower than typically achieved with photon IMRT, warranting development of a planar dosimetry system designed specifically for IMNRT and/or the improvement of neutron beam modeling in the penumbral region. The fractional photon dose component did not change significantly in a typical IMNRT plan versus a conventional fast neutron therapy plan, and IMNRT delivery is not expected to significantly alter the RBE. All other commissioning results were considered satisfactory for clinical implementation of IMNRT, including the external neutron dose validation, which agreed with the predicted neutron dose to within 1%. Conclusions: IMNRT has been successfully commissioned for clinical use. While current plan quality is inferior to photon IMRT, it is superior to conventional fast neutron therapy. Ion chamber validation results for IMNRT commissioning are also comparable to those typically achieved with photon IMRT. Gamma pass rates for planar dose distributions are lower than typically observed for photon IMRT but may be improved with improved planar dosimetry equipment and beam modeling techniques. In the meantime, patient-specific quality assurance measurements should rely more heavily on point dose measurements with tissue equivalent ionization chambers. No significant technical impediments are anticipated in the clinical implementation of IMNRT as described here.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

A Monte Carlo approach to validation of FFF VMAT treatment plans for the TrueBeam linac

Purpose: To commission and benchmark a vendor-supplied (Varian Medical Systems) Monte Carlo phase-space data for the 6 MV flattening filter free (FFF) energy mode on a TrueBeam linear accelerator for the purpose of quality assurance of clinical volumetric modulated arc therapy (VMAT) treatment plans. A method for rendering the phase-space data compatible with BEAMnrc/DOSXYZnrc simulation software package is presented. Methods: Monte Carlo (MC) simulations were performed to benchmark the TrueBeam 6 MV FFF phase space data that have been released by the Varian MC Research team. The simulations to benchmark the phase space data were done in three steps. First, the original phase space which was created on a cylindrical surface was converted into a format that was compatible with BEAMnrc. Second, BEAMnrc was used to create field size specific phase spaces located underneath the jaws. Third, doses were calculated with DOSXYZnrc in a water phantom for fields ranging from 1 Multiplication-Sign 1 to 40 Multiplication-Sign 40 cm{sup 2}. Calculated percent depth doses (PDD), transverse profiles, and output factors were compared with measurements for all the fields simulated. After completing the benchmarking study, three stereotactic body radiotherapy (SBRT) VMAT plans created with the Eclipse treatment planning system (TPS) were calculated with Monte Carlo. Ion chamber and film measurements were also performed on these plans. 3D gamma analysis was used to compare Monte Carlo calculation with TPS calculations and with film measurement. Results: For the benchmarking study, MC calculated and measured values agreed within 1% and 1.5% for PDDs and in-field transverse profiles, respectively, for field sizes >1 Multiplication-Sign 1 cm{sup 2}. Agreements in the 80%-20% penumbra widths were better than 2 mm for all the fields that were compared. With the exception of the 1 Multiplication-Sign 1 cm{sup 2} field, the agreement between measured and calculated output factors was within 1%. It is of note that excellent agreement in output factors for all field sizes including highly asymmetric fields was achieved without accounting for backscatter into the beam monitor chamber. For the SBRT VMAT plans, the agreement between Monte Carlo and ion chamber point dose measurements was within 1%. Excellent agreement between Monte Carlo, treatment planning system and Gafchromic film dose distribution was observed with over 99% of the points in the high dose volume passing the 3%, 3 mm gamma test. Conclusions: The authors have presented a method for making the Varian IAEA compliant 6 MV FFF phase space file of the TrueBeam linac compatible with BEAMnrc/DOSXYZnrc. After benchmarking the modified phase space against measurement, they have demonstrated its potential for use in MC based quality assurance of complex delivery techniques.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Development and clinical evaluation of automatic fiducial detection for tumor tracking in cine megavoltage images during volumetric modulated arc therapy

Purpose: Real-time tracking of implanted fiducials in cine megavoltage (MV) imaging during volumetric modulated arc therapy (VMAT) delivery is complicated due to the inherent low contrast of MV images and potential blockage of dynamic leaves configurations. The purpose of this work is to develop a clinically practical autodetection algorithm for motion management during VMAT. Methods: The expected field-specific segments and the planned fiducial position from the Eclipse (Varian Medical Systems, Palo Alto, CA) treatment planning system were projected onto the MV images. The fiducials were enhanced by applying a Laplacian of Gaussian filter in the spatial domain for each image, with a blob-shaped object as the impulse response. The search of implanted fiducials was then performed on a region of interest centered on the projection of the fiducial when it was within an open field including the case when it was close to the field edge or partially occluded by the leaves. A universal template formula was proposed for template matching and normalized cross correlation was employed for its simplicity and computational efficiency. The search region for every image was adaptively updated through a prediction model that employed the 3D position of the fiducial estimated from the localized positions in previous images. This prediction model allowed the actual fiducial position to be tracked dynamically and was used to initialize the search region. The artifacts caused by electronic interference during the acquisition were effectively removed. A score map was computed by combining both morphological information and image intensity. The pixel location with the highest score was selected as the detected fiducial position. The sets of cine MV images taken during treatment were analyzed with in-house developed software written in MATLAB (The Mathworks, Inc., Natick, MA). Five prostate patients were analyzed to assess the algorithm performance by measuring their positioning accuracy during treatment. Results: The algorithm was able to accurately localize the fiducial position on MV images with success rates of more than 90% per case. The percentage of images in which each fiducial was localized in the studied cases varied between 23% and 65%, with at least one fiducial having been localized between 40% and 95% of the images. This depended mainly on the modulation of the plan and fiducial blockage. The prostate movement in the presented cases varied between 0.8 and 3.5 mm (mean values). The maximum displacement detected among all patients was of 5.7 mm. Conclusions: An algorithm for automatic detection of fiducial markers in cine MV images has been developed and tested with five clinical cases. Despite the challenges posed by complex beam aperture shapes, fiducial localization close to the field edge, partial occlusion of fiducials, fast leaf and gantry movement, and inherently low MV image quality, good localization results were achieved in patient images. This work provides a technique for enabling real-time accurate fiducial detection and tumor tracking during VMAT treatments without the use of extra imaging dose.

60 APPLIED LIFE SCIENCES↗