Search NASA⌕ Search

SEARCH · Search NASA

Results for “matched field processing”

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 883 records · Page 49

Effects of pulse duration on magnetostimulation thresholds

Purpose: Medical imaging techniques such as magnetic resonance imaging and magnetic particle imaging (MPI) utilize time-varying magnetic fields that are subject to magnetostimulation limits, which often limit the speed of the imaging process. Various human-subject experiments have studied the amplitude and frequency dependence of these thresholds for gradient or homogeneous magnetic fields. Another contributing factor was shown to be number of cycles in a magnetic pulse, where the thresholds decreased with longer pulses. The latter result was demonstrated on two subjects only, at a single frequency of 1.27 kHz. Hence, whether the observed effect was due to the number of cycles or due to the pulse duration was not specified. In addition, a gradient-type field was utilized; hence, whether the same phenomenon applies to homogeneous magnetic fields remained unknown. Here, the authors investigate the pulse duration dependence of magnetostimulation limits for a 20-fold range of frequencies using homogeneous magnetic fields, such as the ones used for the drive field in MPI. Methods: Magnetostimulation thresholds were measured in the arms of six healthy subjects (age: 27 ± 5 yr). Each experiment comprised testing the thresholds at eight different pulse durations between 2 and 125 ms at a single frequency, which took approximately 30–40 min/subject. A total of 34 experiments were performed at three different frequencies: 1.2, 5.7, and 25.5 kHz. A solenoid coil providing homogeneous magnetic field was used to induce stimulation, and the field amplitude was measured in real time. A pre-emphasis based pulse shaping method was employed to accurately control the pulse durations. Subjects reported stimulation via a mouse click whenever they felt a twitching/tingling sensation. A sigmoid function was fitted to the subject responses to find the threshold at a specific frequency and duration, and the whole procedure was repeated at all relevant frequencies and pulse durations. Results: The magnetostimulation limits decreased with increasing pulse duration (T{sub pulse}). For T{sub pulse} < 18 ms, the thresholds were significantly higher than at the longest pulse durations (p < 0.01, paired Wilcoxon signed-rank test). The normalized magnetostimulation threshold (B{sub Norm}) vs duration curve at all three frequencies agreed almost identically, indicating that the observed effect is independent of the operating frequency. At the shortest pulse duration (T{sub pulse} ≈ 2 ms), the thresholds were approximately 24% higher than at the asymptotes. The thresholds decreased to within 4% of their asymptotic values for T{sub pulse} > 20 ms. These trends were well characterized (R{sup 2} = 0.78) by a stretched exponential function given by B{sub Norm}=1+αe{sup −(T{sub p}{sub u}{sub l}{sub s}{sub e}/β){sup γ}}, where the fitted parameters were α = 0.44, β = 4.32, and γ = 0.60. Conclusions: This work shows for the first time that the magnetostimulation thresholds decrease with increasing pulse duration, and that this effect is independent of the operating frequency. Normalized threshold vs duration trends are almost identical for a 20-fold range of frequencies: the thresholds are significantly higher at short pulse durations and settle to within 4% of their asymptotic values for durations longer than 20 ms. These results emphasize the importance of matching the human-subject experiments to the imaging conditions of a particular setup. Knowing the dependence of the safety limits to all contributing factors is critical for increasing the time-efficiency of imaging systems that utilize time-varying magnetic fields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

First upper limits on the radar cross section of cosmic-ray induced extensive air showers

TARA (Telescope Array Radar) is a cosmic ray radar detection experiment colocated with Telescope Array, the conventional surface scintillation detector (SD) and fluorescence telescope detector (FD) near Delta, Utah, U.S.A. Furthermore, the TARA detector combines a 40 kW, 54.1 MHz VHF transmitter and high-gain transmitting antenna which broadcasts the radar carrier over the SD array and within the FD field of view, towards a 250 MS/s DAQ receiver. TARA has been collecting data since 2013 with the primary goal of observing the radar signatures of extensive air showers (EAS). Simulations indicate that echoes are expected to be short in duration (~10 µs) and exhibit rapidly changing frequency, with rates on the order 1 MHz/µs. The EAS radar cross-section (RCS) is currently unknown although it is the subject of over 70 years of speculation. One novel signal search technique is described in which the expected radar echo of a particular air shower is used as a matched filter template and compared to waveforms obtained by triggering the radar DAQ using the Telescope Array fluorescence detector. No evidence for the scattering of radio frequency radiation by EAS is obtained to date. Finally, we report the first quantitative RCS upper limits using EAS that triggered the Telescope Array Fluorescence Detector.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Three Paradigms of Lunar Regolith Evolution

Integration of diverse datasets on the Moon may render some paradigms of lunar science either better-defended or vulnerable. We will consider three paradigms commonly used for understanding the processes of lunar regolith evolution in light of new and accumulated data. Our premise is that all data-sets should converge to a single interpretation if a concept or model is to be accepted as a paradigm. If a convergence is lacking, the paradigm needs fresh scrutiny. SteadyState: Lunar regolith evolution is currently understood in terms of comminution, agglutination, and replenishment as described by McKay and coworkers). Briefly, the model envisages continued micrometeoritic bombardment to comminute exposed soil particles to finer sizes while continued agglutination consumes finer sizes to produce larger constructional particles. Eventually, a balance between these two opposing processes achieves a steady state; soils at steady state maintain their mean grain size (M(sub z)). Episodic higher-energy impacts excavate fresh coarse material from below the soil cover, disturb the steady state, and restart the process to achieve a new steady state. It follows that the thickness of the regolith at any site would control the frequency of replenishment; indeed, the thickness of the regolith at Apollo landing sites was predicted by McKay et al. from the average M(sub z) of local soils. However, replenishment may come also from disintegrating boulders and cobbles at the lunar surface, and rates of comminution and agglutination may depend on the properties of target material. Regression between M(sub z) and I(sub s)/Fe(sup 0) (a measure of maturity or total surface exposure) of Apollo soils at different sites shows the following relations and estimated M(sub z) at a high maturity of I(sub s)/Fe(sup 0)= 100. It is possible that Apollo 12 and 15 sites have the thickest regolith and the Apollo 16 site has the thinnest. It is also possible that Apollo 12 and 15 basalts are comminuted faster than Apollo 16 highland rocks and Apollo 14 and 17 soils are products of mixed parentage. If a soil becomes continually finer as it matures until agglutination catches up, and if comminution is differential-dependent on the physical properties of the constituents, then the composition of the bulk soil has to match the composition of some "fulcrum" grain size fraction, say X Grain size fractions >X and <X will complement each other; their mass balance is the bulk soil. It appears that the 10-20-micron size fraction may be the fulcrum. In general, trace-element chemistry and IR reflectance spectra of this size fraction are closest to that of the bulk soil, regardless of maturity that is surprising. Disaggregated products of regolith breccias may also show similar relationships. If the 10-20 gm is the fulcrum (i.e., X as above) for many soil properties (e.g., major element composition, FMR, solar-wind-implanted elements), then this may be the ultimate mean grain size of lunar soils at steady state. However, different properties of soils may find steady states at different grain size fractions. The steady state of solar-wind-implanted elements, on the other hand, will climb up the grain-size scale as agglutinates transfer surface-correlated components into volume correlated components until a saturation level is reached or the rates of replenishment and implantation become equal. The same will be the case with vapor-deposited reduced metals as they too are incorporated inside constructional particles. Properties that are directly affected by soil-maturation processes will thus have different pathways of achieving steady states. Maturity, i.e., cumulative surface exposure, of lunar soils is best quantified by the amount of nanophase superparamagnetic Fe(sup 0) (np-Fe(sup 0)) normalized to Fe content (=I(sub s)/Fe(sup 0). The majority consensus (paradigm?) for the production of np-Fe(sup 0) is associated with the production of agglutinates. Because large doses of solar-wind H are implanted in all lunar soils upon exposure, any melting (e.g., during agglutinate production) triggers a chemical reduction of Fe-bearing minerals resulting in np-Fe(sup 0) production. The quantity of np-Fe(sup 0) is thus dependent on melting events, (i.e., exposure), and limited by the Fe content of the soil. All freshly produced np-Fe(sup 0) resides in agglutinitic glass, as new TEM images show. Apparently, the correction procedure developed by Lucey et al. to estimate the Fe content of the lunar surface from IR-reflectance spectra depends on accepting the above. However, the process of producing np-Fe(sup 0) may be physical rather than chemical. All np-Fe(sup 0) could be deposits from a vapor produced by micrometeoritic impact on lunar soils. If metal-O bonds in target phases are broken, O being "most volatile" will escape leaving an O-deficient vapor to facilitate the production of np-Fe(sup 0). If so, the quantity of np-Fe(sup 0) is dependent on the vaporizing events, (i.e., exposure), and limited by the efficiency of breaking metal-O bonds and the escape of 0. To the extent that strengths of metal-O bonds are dependent on the local crystal field, production of np-Fe(sup 0) may be limited by the mineral composition of target soils and not by their total Fe content. According to this model, vapor-deposited np-Fe(sup 0) should be found at any retentive sites on lunar soil grains. Indeed, TEM images show np-Fe(sup 0) on plagioclase and ilmenite. Incorporation of such pre-irradiated np-Fe(sup 0)-bearing grains into agglutinates may account for eventual increased emplacement of np-Fe(sup 0) in agglutinates. Such a paradigm shift in understanding the origin of np-Fe(sup 0) will raise questions ranging from the unquestionable use of Is/FeO as the universal maturity parameter of lunar soils to global elemental maps of the Moon from remote-sensing data. Additional information is contained in the original.

Basu, A.↗

Modelling and computational improvements to the simulation of single vector-boson plus jet processes for the ATLAS experiment

This paper presents updated Monte Carlo configurations used to model the production of single electroweak vector bosons (W, Z/γ$^{∗}$) in association with jets in proton-proton collisions for the ATLAS experiment at the Large Hadron Collider. Improvements pertaining to the electroweak input scheme, parton-shower splitting kernels and scale-setting scheme are shown for multi-jet merged configurations accurate to next-to-leading order in the strong and electroweak couplings. The computational resources required for these set-ups are assessed, and approximations are introduced resulting in a factor three reduction of the per-event CPU time without affecting the physics modelling performance. Continuous statistical enhancement techniques are introduced by ATLAS in order to populate low cross-section regions of phase space and are shown to match or exceed the generated effective luminosity. This, together with the lower per-event CPU time, results in a 50% reduction in the required computing resources compared to a legacy set-up previously used by the ATLAS collaboration. The set-ups described in this paper will be used for future ATLAS analyses and lay the foundation for the next generation of Monte Carlo predictions for single vector-boson plus jets production.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Limits to Gauge Coupling in the Dark Sector Set by the Nonobservation of Instanton-Induced Decay of Super-Heavy Dark Matter in the Pierre Auger Observatory Data

Instantons, which are nonperturbative solutions to Yang-Mills equations, provide a signal for the occurrence of quantum tunneling between distinct classes of vacua. They can give rise to decays of particles otherwise forbidden. Using data collected at the Pierre Auger Observatory, we search for signatures of such instanton-induced processes that would be suggestive of super-heavy particles decaying in the Galactic halo. These particles could have been produced during the post-inflationary epoch and match the relic abundance of dark matter inferred today. The nonobservation of the signatures searched for allows us to derive a bound on the reduced coupling constant of gauge interactions in the dark sector: α X ≲ 0.09, for 10 9 ≲ M X /GeV < 10 19 . Conversely, we obtain that, for instance, a reduced coupling constant α X = 0.09 excludes masses M X ≳ 3 x 10 13 GeV. In the context of dark matter production from gravitational interactions alone, we illustrate how these bounds are complementary to those obtained on the Hubble rate at the end of inflation from the nonobservation of tensor modes in the cosmological microwave background.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A temperature profile diagnostic for radiation waves on OMEGA-60

Predicting and matching radiation wave propagation with computational models has proven difficult. Information provided by experiments studying radiation flow has been limited when only radiation breakout is measured. We have developed the COAX (co-axial) diagnostic platform to provide spatial temperature profiles of a radiation wave through low density foams as a more detailed constraint for simulations. COAX uses a standard, laser-driven OMEGA-60 halfraum to drive radiation down a titanium-laden silicon oxide foam. Point-projection X-ray absorption spectroscopy perpendicular to the radiation flow measures the spatial profile of titanium ionization. The spectroscopic measurement utilizes a broadband capsule backlighter. Imaging and streak spectroscopy are used to characterize the size and spectrum of this source. Radiography provides an additional constraint by capturing the developing shock as the radiation flow becomes subsonic. The DANTE diagnostic is used to measure the halfraum temperature. Here, we provide a spectroscopic analysis of COAX data to determine temperature, and we describe experimental sources of uncertainty. The temperature is obtained by comparison to multi-temperature synthetic spectra post-processed from radiation-hydrodynamics simulations. Quantitative comparison between data and synthetic spectra generated from temperature profiles at relevant simulation times enable determination of a peak temperature of 114 ± 8 eV at 265 ± 22.4 μm from the halfraum. This represents an improvement over the temperature uncertainties of previous radiation flow experiments. Further refinements to the spectroscopic analysis could achieve ±4 eV. The combination between space-resolved spectroscopy and radiography enables us to determine the distance from the halfraum of both the radiation front and the shock front at the time of measurement. For the example shown in this paper the radiation front position is 600–630 μm at 3.43 ± 0.16 ns and the shock front position is 633 μm at 3.3 ± 0.24 ns.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Cosmological implications of photon-flux upper limits at ultrahigh energies in scenarios of Planckian-interacting massive particles for dark matter

Here using the data of the Pierre Auger Observatory, we report on a search for signatures that would be suggestive of super-heavy particles decaying in the Galactic halo. From the lack of signal, we present upper limits for different energy thresholds above ≳10 8 GeV on the secondary by-product fluxes expected from the decay of the particles. Assuming that the energy density of these super-heavy particles matches that of dark matter observed today, we translate the upper bounds on the particle fluxes into tight constraints on the couplings governing the decay process as a function of the particle mass. Instantons, which are nonperturbative solutions to Yang-Mills equations, can give rise to decay channels otherwise forbidden and transform stable particles into metastable ones. Assuming such instanton-induced decay processes, we derive a bound on the reduced coupling constant of gauge interactions in the dark sector: α X ≲ 0.09, for 10 9 ≲ M X /GeV < 10 19 . Conversely, we obtain that, for instance, a reduced coupling constant α X = 0.09 excludes masses M X ≳ 3 x 10 13 GeV . In the context of dark matter production from gravitational interactions alone during the reheating epoch, we derive constraints on the parameter space that involves, in addition to M X and α X , the Hubble rate at the end of inflation, the reheating efficiency, and the nonminimal coupling of the Higgs with curvature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

SU-E-T-24: A Simple Correction-Based Method for Independent Monitor Unit (MU) Verification in Monte Carlo (MC) Lung SBRT Plans

Purpose: Lung-SBRT uses hypo-fractionated dose in small non-IMRT fields with tissue-heterogeneity corrected plans. An independent MU verification is mandatory for safe and effective delivery of the treatment plan. This report compares planned MU obtained from iPlan-XVM-Calgorithm against spreadsheet-based hand-calculation using most commonly used simple TMR-based method. Methods: Treatment plans of 15 patients who underwent for MC-based lung-SBRT to 50Gy in 5 fractions for PTV V100%=95% were studied. ITV was delineated on MIP images based on 4D-CT scans. PTVs(ITV+5mm margins) ranged from 10.1- 106.5cc(average=48.6cc). MC-SBRT plans were generated using a combination of non-coplanar conformal arcs/beams using iPlan XVM-Calgorithm (BrainLAB iPlan ver.4.1.2) for Novalis-TX consisting of micro-MLCs and 6MV-SRS (1000MU/min) beam. These plans were re-computed using heterogeneity-corrected Pencil-Beam (PB-hete) algorithm without changing any beam parameters, such as MLCs/MUs. Dose-ratio: PB-hete/MC gave beam-by-beam inhomogeneity-correction-factors (ICFs):Individual Correction. For independent-2nd-check, MC-MUs were verified using TMR-based hand-calculation and obtained an average ICF:Average Correction, whereas TMR-based hand-calculation systematically underestimated MC-MUs by ∼5%. Also, first 10 MC-plans were verified with an ion-chamber measurement using homogenous phantom. Results: For both beams/arcs, mean PB-hete dose was systematically overestimated by 5.5±2.6% and mean hand-calculated MU systematic underestimated by 5.5±2.5% compared to XVMC. With individual correction, mean hand-calculated MUs matched with XVMC by - 0.3±1.4%/0.4±1.4 for beams/arcs, respectively. After average 5% correction, hand-calculated MUs matched with XVMC by 0.5±2.5%/0.6±2.0% for beams/arcs, respectively. Smaller dependence on tumor volume(TV)/field size(FS) was also observed. Ion-chamber measurement was within ±3.0%. Conclusion: PB-hete overestimates dose to lung tumor relative to XVMC. XVMC-algorithm is much more-complex and accurate with tissues-heterogeneities. Measurement at machine is time consuming and need extra resources; also direct measurement of dose for heterogeneous treatment plans is not clinically practiced, yet. This simple correction-based method was very helpful for independent-2nd-check of MC-lung-SBRT plans and routinely used in our clinic. A look-up table can be generated to include TV/FS dependence in ICFs.

60 APPLIED LIFE SCIENCES↗

Noise-resolution tradeoffs in x-ray CT imaging: A comparison of penalized alternating minimization and filtered backprojection algorithms

Purpose: In comparison with conventional filtered backprojection (FBP) algorithms for x-ray computed tomography (CT) image reconstruction, statistical algorithms directly incorporate the random nature of the data and do not assume CT data are linear, noiseless functions of the attenuation line integral. Thus, it has been hypothesized that statistical image reconstruction may support a more favorable tradeoff than FBP between image noise and spatial resolution in dose-limited applications. The purpose of this study is to evaluate the noise-resolution tradeoff for the alternating minimization (AM) algorithm regularized using a nonquadratic penalty function. Methods: Idealized monoenergetic CT projection data with Poisson noise were simulated for two phantoms with inserts of varying contrast (7%-238%) and distance from the field-of-view (FOV) center (2-6.5 cm). Images were reconstructed for the simulated projection data by the FBP algorithm and two penalty function parameter values of the penalized AM algorithm. Each algorithm was run with a range of smoothing strengths to allow quantification of the noise-resolution tradeoff curve. Image noise is quantified as the standard deviation in the water background around each contrast insert. Modulation transfer functions (MTFs) were calculated from six-parameter model fits to oversampled edge-spread functions defined by the circular contrast-insert edges as a metric of local resolution. The integral of the MTF up to 0.5 lp/mm was adopted as a single-parameter measure of local spatial resolution. Results: The penalized AM algorithm noise-resolution tradeoff curve was always more favorable than that of the FBP algorithm. While resolution and noise are found to vary as a function of distance from the FOV center differently for the two algorithms, the ratio of noises when matching the resolution metric is relatively uniform over the image. The ratio of AM-to-FBP image variances, a predictor of dose-reduction potential, was strongly dependent on the shape of the AM's nonquadratic penalty function and was also strongly influenced by the contrast of the insert for which resolution is quantified. Dose-reduction potential, reported here as the fraction (%) of FBP dose necessary for AM to reconstruct an image with comparable noise and resolution, for one penalty parameter value of the AM algorithm was found to vary from 70% to 50% for low-contrast and high-contrast structures, respectively, and from 70% to 10% for the second AM penalty parameter value. However, the second penalty, AM-700, was found to suffer from poor low-contrast resolution when matching the high-contrast resolution metric with FBP. Conclusions: The results of this simulation study imply that penalized AM has the potential to reconstruct images with similar noise and resolution using a fraction (10%-70%) of the FBP dose. However, this dose-reduction potential depends strongly on the AM penalty parameter and the contrast magnitude of the structures of interest. In addition, the authors' results imply that the advantage of AM can be maximized by optimizing the nonquadratic penalty function to the specific imaging task of interest. Future work will extend the methods used here to quantify noise and resolution in images reconstructed from real CT data.

60 APPLIED LIFE SCIENCES↗

Observation of proton-tagged, central (semi)exclusive production of high-mass lepton pairs in pp collisions at 13 TeV with the CMS-TOTEM precision proton spectrometer

The process pp → pℓ$^{+}$ℓ$^{−}$p$^{(*)}$, with ℓ$^{+}$ℓ$^{−}$ a muon or an electron pair produced at midrapidity with mass larger than 110 GeV, has been observed for the first time at the LHC in pp collisions at $ \sqrt{s}=13 $ TeV. One of the two scattered protons is measured in the CMS-TOTEM precision proton spectrometer (CT-PPS), which operated for the first time in 2016. The second proton either remains intact or is excited and then dissociates into a low-mass state p$^{*}$, which is undetected. The measurement is based on an integrated luminosity of 9.4 fb$^{−1}$ collected during standard, high-luminosity LHC operation. A total of 12 μ$^{+}$μ$^{−}$ and 8 e$^{+}$e$^{−}$ pairs with m(ℓ$^{+}$ℓ$^{−}$) > 110 GeV, and matching forward proton kinematics, are observed, with expected backgrounds of 1.49 ± 0.07 (stat) ± 0.53 (syst) and 2.36 ± 0.09 (stat) ± 0.47 (syst), respectively. This corresponds to an excess of more than five standard deviations over the expected background. The present result constitutes the first observation of proton-tagged γγ collisions at the electroweak scale. This measurement also demonstrates that CT-PPS performs according to the design specifications.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Feasibility and robustness of dose painting by numbers in proton therapy with contour-driven plan optimization

Purpose: To prove the ability of protons to reproduce a dose gradient that matches a dose painting by numbers (DPBN) prescription in the presence of setup and range errors, by using contours and structure-based optimization in a commercial treatment planning system. Methods: For two patients with head and neck cancer, voxel-by-voxel prescription to the target volume (GTV{sub PET}) was calculated from {sup 18}FDG-PET images and approximated with several discrete prescription subcontours. Treatments were planned with proton pencil beam scanning. In order to determine the optimal plan parameters to approach the DPBN prescription, the effects of the scanning pattern, number of fields, number of subcontours, and use of range shifter were separately tested on each patient. Different constant scanning grids (i.e., spot spacing = Δx = Δy = 3.5, 4, and 5 mm) and uniform energy layer separation [4 and 5 mm WED (water equivalent distance)] were analyzed versus a dynamic and automatic selection of the spots grid. The number of subcontours was increased from 3 to 11 while the number of beams was set to 3, 5, or 7. Conventional PTV-based and robust clinical target volumes (CTV)-based optimization strategies were considered and their robustness against range and setup errors assessed. Because of the nonuniform prescription, ensuring robustness for coverage of GTV{sub PET} inevitably leads to overdosing, which was compared for both optimization schemes. Results: The optimal number of subcontours ranged from 5 to 7 for both patients. All considered scanning grids achieved accurate dose painting (1% average difference between the prescribed and planned doses). PTV-based plans led to nonrobust target coverage while robust-optimized plans improved it considerably (differences between worst-case CTV dose and the clinical constraint was up to 3 Gy for PTV-based plans and did not exceed 1 Gy for robust CTV-based plans). Also, only 15% of the points in the GTV{sub PET} (worst case) were above 5% of DPBN prescription for robust-optimized plans, while they were more than 50% for PTV plans. Low dose to organs at risk (OARs) could be achieved for both PTV and robust-optimized plans. Conclusions: DPBN in proton therapy is feasible with the use of a sufficient number subcontours, automatically generated scanning patterns, and no more than three beams are needed. Robust optimization ensured the required target coverage and minimal overdosing, while PTV-approach led to nonrobust plans with excessive overdose. Low dose to OARs can be achieved even in the presence of a high-dose escalation as in DPBN.

60 APPLIED LIFE SCIENCES↗

The Distribution of Titanium in Lunar Soils on the Basis of Sensor and In Situ Data Fusion

A variety of remote-sensing measurements have been used to map the distribution of elements on the Moon as a means of providing constraints on the processes from which its crust and major terranes originated. Discussed here is Ti, which is incorporated into refractory minerals such as ilmenite during the latter stages of differentiation, and is thus a most useful element for understanding mare basalt petrogenesis. One of the earliest Ti maps showed Ti variations in nearside maria on the basis of groundbased spectral reflectance measurements. A map of Ti derived from gamma-ray measurements on Apollo 15 and 16 was produced at about the same time, and was improved upon considerably by Davis and coworkers, who effectively removed sources of spurious variation from Fe and Al or REE (e.g., Th) interference, and calibrated Ti on the bases of landing-site soil averages. In recent years, spectral reflectance measurements from Clementine have been used by Lucey and coworkers to produce global Ti distribution maps as well. As we indicated previously, the Lucey and Davis maps agree to first order. Meanwhile, we are using the concept of sensor data fusion to combine measurements from the AGR (Apollo gamma-ray) and CSR (Clementine Spectral Reflectance) techniques with ground truth from lunar soils to utilize the differences between the two maps to understand the distribution of Ti within lunar soil components, as we have done with Fe. This technique should be verified and applied on Lunar Prospector gamma-ray measurements of Ti, as the calibrated data become available within the next couple of years. Lunar Ti is found principally in the mineral ilmenite, and is associated with certain components of lunar soil: crystalline Ilmenite mineral fragments and high Ti-bearing glass. All data indicate that Ti is associated with maria and mafic minerals. In AGR and CSR datasets, Ti is highest on the nearside and in the maria, particularly in southern Serenitatis/northern Tranquillitatis and northwestern Procellarum. Unlike CSR-derived values, AGR Ti values show modest increase (up to 0.7%) on the northern farside. Both techniques show a primarily unimodal distribution with a shoulder in the high Ti direction and a primary mode at approximately 0.2% (0.15% for CSR and 0.25% for AGR data), but the AGR Ti data have more structure in the shoulder, including an apparent minor mode at 1.2% Ti, representing non-highland areas. Recalibrating the two datasets on the basis of matching the peaks and range on the histograms would not account for this additional structure in the AGR data. Low Ti areas, which occur predominantly on the farside, are not represented well in the lunar sample collection, although several meteorites thought to be of lunar highland origin show Ti abundances averaging about 0.2%. (Error bars on AGR Ti values are about 0.5%.) Gamma-ray measurements reflect intrinsic Ti surface composition, regardless of the physical or chemical state of Ti. Lucey and coworkers have attempted to produce an equivalent bulk Ti map by normalizing the spectral feature at 415 nm to remove the effects of physical variations from soil to soil. The normalization is based on laboratory measurements of available samples, primarily nearside maria, and is undoubtedly optimized for the proportions in which Ti- components are found in these soils. In particular, mare soils tend to have a much higher proportion of opaque mineral grains, which appears to be well correlated with Ti abundance. The relationship between Ti abundance and proportion of Ti-bearing glass (weighted for Ti abundance) is not as direct. Also shows the relationship between landing site soil average, AGR, and CSR Ti abundances (with CSR filtered to match resolution and field of view of AGR data). AGR Ti values show the best agreement with Ti soil averages, CSR Ti values being somewhat higher in high-Ti soils, and lower in low-Ti soils. These results would tend to support the argument that the technique for deriving Ti from the CSR measurements is optimized for one of the primary Ti-bearing components, most likely opaque mineral grains, and is thus over-compensating and reporting higher Ti values in mare areas, and underreporting Ti in highland areas with much lower proportions of opaque mineral grains. More information is contained in original.

Clark, P. E.↗

Simultaneous motion estimation and image reconstruction (SMEIR) for 4D cone-beam CT

Purpose: Image reconstruction and motion model estimation in four-dimensional cone-beam CT (4D-CBCT) are conventionally handled as two sequential steps. Due to the limited number of projections at each phase, the image quality of 4D-CBCT is degraded by view aliasing artifacts, and the accuracy of subsequent motion modeling is decreased by the inferior 4D-CBCT. The objective of this work is to enhance both the image quality of 4D-CBCT and the accuracy of motion model estimation with a novel strategy enabling simultaneous motion estimation and image reconstruction (SMEIR).Methods: The proposed SMEIR algorithm consists of two alternating steps: (1) model-based iterative image reconstruction to obtain a motion-compensated primary CBCT (m-pCBCT) and (2) motion model estimation to obtain an optimal set of deformation vector fields (DVFs) between the m-pCBCT and other 4D-CBCT phases. The motion-compensated image reconstruction is based on the simultaneous algebraic reconstruction technique (SART) coupled with total variation minimization. During the forward- and backprojection of SART, measured projections from an entire set of 4D-CBCT are used for reconstruction of the m-pCBCT by utilizing the updated DVF. The DVF is estimated by matching the forward projection of the deformed m-pCBCT and measured projections of other phases of 4D-CBCT. The performance of the SMEIR algorithm is quantitatively evaluated on a 4D NCAT phantom. The quality of reconstructed 4D images and the accuracy of tumor motion trajectory are assessed by comparing with those resulting from conventional sequential 4D-CBCT reconstructions (FDK and total variation minimization) and motion estimation (demons algorithm). The performance of the SMEIR algorithm is further evaluated by reconstructing a lung cancer patient 4D-CBCT.Results: Image quality of 4D-CBCT is greatly improved by the SMEIR algorithm in both phantom and patient studies. When all projections are used to reconstruct a 3D-CBCT by FDK, motion-blurring artifacts are present, leading to a 24.4% relative reconstruction error in the NACT phantom. View aliasing artifacts are present in 4D-CBCT reconstructed by FDK from 20 projections, with a relative error of 32.1%. When total variation minimization is used to reconstruct 4D-CBCT, the relative error is 18.9%. Image quality of 4D-CBCT is substantially improved by using the SMEIR algorithm and relative error is reduced to 7.6%. The maximum error (MaxE) of tumor motion determined from the DVF obtained by demons registration on a FDK-reconstructed 4D-CBCT is 3.0, 2.3, and 7.1 mm along left–right (L-R), anterior–posterior (A-P), and superior–inferior (S-I) directions, respectively. From the DVF obtained by demons registration on 4D-CBCT reconstructed by total variation minimization, the MaxE of tumor motion is reduced to 1.5, 0.5, and 5.5 mm along L-R, A-P, and S-I directions. From the DVF estimated by SMEIR algorithm, the MaxE of tumor motion is further reduced to 0.8, 0.4, and 1.5 mm along L-R, A-P, and S-I directions, respectively.Conclusions: The proposed SMEIR algorithm is able to estimate a motion model and reconstruct motion-compensated 4D-CBCT. The SMEIR algorithm improves image reconstruction accuracy of 4D-CBCT and tumor motion trajectory estimation accuracy as compared to conventional sequential 4D-CBCT reconstruction and motion estimation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Precise measurements of W - and Z -boson transverse momentum spectra with the ATLAS detector using pp collisions at $\sqrt{s} = 5.02$ TeV and 13 TeV

This paper describes measurements of the transverse momentum spectra of W and Z bosons produced in proton–proton collisions at centre-of-mass energies of $\sqrt{s}$ = 5.02 TeV and $\sqrt{s}$ = 13 TeV with the ATLAS experiment at the Large Hadron Collider. Measurements are performed in the electron and muon channels, W → $\ell$$v$ and Z → $\ell$$\ell$ ($\ell$ = e or μ), and for W events further separated by charge. The data were collected in 2017 and 2018, in dedicated runs with reduced instantaneous luminosity, and correspond to 255 and 338 pb -1 at $\sqrt{s}$ = 5.02 TeV and 13 TeV, respectively. These conditions optimise the reconstruction of the W-boson transverse momentum. The distributions observed in the electron and muon channels are unfolded, combined, and compared to QCD calculations based on parton shower Monte Carlo event generators and analytical resummation. The description of the transverse momentum distributions by Monte Carlo event generators is imperfect and shows significant differences largely common to W - , W + and Z production. The agreement is better at $\sqrt{s}$ = 5.02 TeV, especially for predictions that were tuned to Z production data at $\sqrt{s}$ = 7 TeV. Higher-order, resummed predictions based on DYTurbo generally match the data best across the spectra. Distribution ratios are also presented and test the understanding of differences between the production processes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Heterogeneity of Bulk Oxygen Isotopic Compositions in Anhydrous Interplanetary Dust Particles

Introduction: Anhydrous interplanetary dust particles are one of the least altered ancient solar systemmaterials. Their bulk chemical compositions match those of CI chondrites within a factor of 2-3, except for carbonwhich is enriched in both anhydrous and hydrated IDPs by ~4x CI [1]. Hydrous IDPs often show 16O-poor isotopiccompositions, likely due to the interaction with isotopically heavy H2O [2]. Anhydrous IDPs are interpreted tooriginate from comets from the outer solar system [3]. Oxygen isotopic data for these materials is scarce, but resultspoint to a wider range of compositions [4]. This heterogeneity might be a result of mixing between a 16O-rich and a16O-poor reservoir formed via self-shielding and photodissociation of CO [5]. Here, we report bulk oxygen isotopiccompositions of three anhydrous IDPs in context of their mineralogical composition and discuss possible origins ofisotopic heterogeneity. Experimental: Bulk chemical composition and mineralogy of 70 nm thin ultramicrotomed sections of threeanhydrous IDPs (L2099A7, L2099A8, and L2071AB1) were analyzed using a Thermo Scientific Titan Themis G360-300 TEM equipped with a four-quadrant energy-dispersive X-ray detector (Super-X G2) at University ofMünster. Bulk oxygen isotopic compositions were measured using a NanoSIMS 50 at the Max-Planck-Institut fürChemie in Mainz. A focused Cs+ ion beam (~1 pA, ~100 nm) was rastered over a field of view varying from 5 x 5 to6 x 6 μm2, depending on particle size for 15 layers. Negative secondary ions of 16O, 17O, 18O, 12C14N, and 28Si werecollected on electron multipliers. The 16OH contribution to the 17O peak was 1-2 ‰. As a standard, a matrix regionof meteorite CR2 Queen Alexandra Range (QUE) 99177 was used, whose oxygen isotopic composition wasmeasured by [6]. Results and Discussion: Particles A7 and A8 are both fine-grained, consisting mostly of small, 100 nm-sizedequilibrated aggregates (EA). Both particles exhibit discontinuous magnetite rims with thicknesses up to 50 nm inA7 and 100 nm in A8 on the outsides. In Particle AB1, the upper part of the particle consists of a few EAs, while thelower part contains lots of small GEMS grains with 100-200 nm diameter clustered together. Magnetite rims arenearly absent, the rare rims reaching maximum 20 nm thickness, indicating that AB1 is of more primitive naturethan the other two IDPs [1]. Furthermore, it contains several diffuse GEMS-like areas, identified by an amorphoussilicate groundmass with small nano-inclusions of Fe,Ni-metal and Fe-sulfides. All three particles are subsolar forall major element/Si ratios. Mean S/Si of particles A7 and A8 is 0.046 and 0.048, respectively, while it is an order ofmagnitude higher in AB1 (0.184). This reflects the higher degree of thermal alteration in A7 and A8, resulting in theloss of volatile S and oxidation of Fe-sulfides and FeNi-metal to magnetite as present in particle rims. A8 is rather16O-poor with δ17OSMOW = 7.8 ± 3.4 ‰ and δ18OSMOW = 9.5 ± 3.1 ‰ (1σ), followed by A7 with δ17OSMOW = − 0.8 ±4.6 ‰ and δ18OSMOW = −2.2 ± 3.4 ‰. In contrast to that, AB1 has the most 16O-rich bulk composition with δ17OSMOW= −24.2 ± 5.4 ‰ and δ18OSMOW = −25.3 ± 3.6 ‰. All investigated IDPs plot slightly above the CCAM line aspreviously reported for other anhydrous IDPs [4], maybe due to a contribution from circumstellar dust from AGBstars enriched in 17O [7]. The 16O-rich composition of AB1 is similar to the most 16O-rich anhydrous IDP U2015D21measured by [8] which is dominated by GEMS [9] and the GEMS-rich IDP GM4-2 reported by [10]. Nevertheless,it is unlikely that GEMS are the carrier for the 16O-rich composition of AB1 because most GEMS have O isotopiccompositions indistinguishable from terrestrial values, although errors on these analyses are large [11]. The isotopicheterogeneity is best explained by contribution of grains from different regions in the protoplanetary disk, samplingdifferent O isotope reservoirs, created by self-shielding and photodissociation of CO. This process resulted in 16O-rich CO gas and 16O-poor H2O that froze as ice-mantles onto dust grains [12]. Alternatively, the higher degree ofheating in A7 and A8 could have influenced isotopic fractionation, because higher δ17,18O values are reported fromthe regions of the atmosphere where small particles are heated and oxidized, as observed for the thermal alteration ofcosmic spherules [13]. Acknowledgements: We would like to thank NASA Astromaterials Acquisition & Curation Office for providing IDPsamples and DFG for funding this project (VO1816/5-1)

B Schulz↗

Bowtie filters for dedicated breast CT: Theory and computational implementation

Purpose: To design bowtie filters with improved properties for dedicated breast CT to improve image quality and reduce dose to the patient. Methods: The authors present three different bowtie filters designed for a cylindrical 14-cm diameter phantom with a uniform composition of 40/60 breast tissue, which vary in their design objectives and performance improvements. Bowtie design #1 is based on single material spectral matching and produces nearly uniform spectral shape for radiation incident upon the detector. Bowtie design #2 uses the idea of basis material decomposition to produce the same spectral shape and intensity at the detector, using two different materials. Bowtie design #3 eliminates the beam hardening effect in the reconstructed image by adjusting the bowtie filter thickness so that the effective attenuation coefficient for every ray is the same. All three designs are obtained using analytical computational methods and linear attenuation coefficients. Thus, the designs do not take into account the effects of scatter. The authors considered this to be a reasonable approach to the filter design problem since the use of Monte Carlo methods would have been computationally intensive. The filter profiles for a cone-angle of 0° were used for the entire length of each filter because the differences between those profiles and the correct cone-beam profiles for the cone angles in our system are very small, and the constant profiles allowed construction of the filters with the facilities available to us. For evaluation of the filters, we used Monte Carlo simulation techniques and the full cone-beam geometry. Images were generated with and without each bowtie filter to analyze the effect on dose distribution, noise uniformity, and contrast-to-noise ratio (CNR) homogeneity. Line profiles through the reconstructed images generated from the simulated projection images were also used as validation for the filter designs. Results: Examples of the three designs are presented. Initial verification of performance of the designs was done using analytical computations of HVL, intensity, and effective attenuation coefficient behind the phantom as a function of fan-angle with a cone-angle of 0°. The performance of the designs depends only weakly on incident spectrum and tissue composition. For all designs, the dynamic range requirement on the detector was reduced compared to the no-bowtie-filter case. Further verification of the filter designs was achieved through analysis of reconstructed images from simulations. Simulation data also showed that the use of our bowtie filters can reduce peripheral dose to the breast by 61% and provide uniform noise and CNR distributions. The bowtie filter design concepts validated in this work were then used to create a computational realization of a 3D anthropomorphic bowtie filter capable of achieving a constant effective attenuation coefficient behind the entire field-of-view of an anthropomorphic breast phantom. Conclusions: Three different bowtie filter designs that vary in performance improvements were described and evaluated using computational and simulation techniques. Results indicate that the designs are robust against variations in breast diameter, breast composition, and tube voltage, and that the use of these filters can reduce patient dose and improve image quality compared to the no-bowtie-filter case.

60 APPLIED LIFE SCIENCES↗

A technique for estimating 4D-CBCT using prior knowledge and limited-angle projections

Purpose: To develop a technique to estimate onboard 4D-CBCT using prior information and limited-angle projections for potential 4D target verification of lung radiotherapy.Methods: Each phase of onboard 4D-CBCT is considered as a deformation from one selected phase (prior volume) of the planning 4D-CT. The deformation field maps (DFMs) are solved using a motion modeling and free-form deformation (MM-FD) technique. In the MM-FD technique, the DFMs are estimated using a motion model which is extracted from planning 4D-CT based on principal component analysis (PCA). The motion model parameters are optimized by matching the digitally reconstructed radiographs of the deformed volumes to the limited-angle onboard projections (data fidelity constraint). Afterward, the estimated DFMs are fine-tuned using a FD model based on data fidelity constraint and deformation energy minimization. The 4D digital extended-cardiac-torso phantom was used to evaluate the MM-FD technique. A lung patient with a 30 mm diameter lesion was simulated with various anatomical and respirational changes from planning 4D-CT to onboard volume, including changes of respiration amplitude, lesion size and lesion average-position, and phase shift between lesion and body respiratory cycle. The lesions were contoured in both the estimated and “ground-truth” onboard 4D-CBCT for comparison. 3D volume percentage-difference (VPD) and center-of-mass shift (COMS) were calculated to evaluate the estimation accuracy of three techniques: MM-FD, MM-only, and FD-only. Different onboard projection acquisition scenarios and projection noise levels were simulated to investigate their effects on the estimation accuracy.Results: For all simulated patient and projection acquisition scenarios, the mean VPD (±S.D.)/COMS (±S.D.) between lesions in prior images and “ground-truth” onboard images were 136.11% (±42.76%)/15.5 mm (±3.9 mm). Using orthogonal-view 15°-each scan angle, the mean VPD/COMS between the lesion in estimated and “ground-truth” onboard images for MM-only, FD-only, and MM-FD techniques were 60.10% (±27.17%)/4.9 mm (±3.0 mm), 96.07% (±31.48%)/12.1 mm (±3.9 mm) and 11.45% (±9.37%)/1.3 mm (±1.3 mm), respectively. For orthogonal-view 30°-each scan angle, the corresponding results were 59.16% (±26.66%)/4.9 mm (±3.0 mm), 75.98% (±27.21%)/9.9 mm (±4.0 mm), and 5.22% (±2.12%)/0.5 mm (±0.4 mm). For single-view scan angles of 3°, 30°, and 60°, the results for MM-FD technique were 32.77% (±17.87%)/3.2 mm (±2.2 mm), 24.57% (±18.18%)/2.9 mm (±2.0 mm), and 10.48% (±9.50%)/1.1 mm (±1.3 mm), respectively. For projection angular-sampling-intervals of 0.6°, 1.2°, and 2.5° with the orthogonal-view 30°-each scan angle, the MM-FD technique generated similar VPD (maximum deviation 2.91%) and COMS (maximum deviation 0.6 mm), while sparser sampling yielded larger VPD/COMS. With equal number of projections, the estimation results using scattered 360° scan angle were slightly better than those using orthogonal-view 30°-each scan angle. The estimation accuracy of MM-FD technique declined as noise level increased.Conclusions: The MM-FD technique substantially improves the estimation accuracy for onboard 4D-CBCT using prior planning 4D-CT and limited-angle projections, compared to the MM-only and FD-only techniques. It can potentially be used for the inter/intrafractional 4D-localization verification.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Angle-Resolved Photoemission Spectroscopy on Electronic Structure and Electron-Phonon Coupling in Cuprate Superconductors

In addition to the record high superconducting transition temperature (T{sub c}), high temperature cuprate superconductors are characterized by their unusual superconducting properties below T{sub c}, and anomalous normal state properties above T{sub c}. In the superconducting state, although it has long been realized that superconductivity still involves Cooper pairs, as in the traditional BCS theory, the experimentally determined d-wave pairing is different from the usual s-wave pairing found in conventional superconductors. The identification of the pairing mechanism in cuprate superconductors remains an outstanding issue. The normal state properties, particularly in the underdoped region, have been found to be at odd with conventional metals which is usually described by Fermi liquid theory; instead, the normal state at optimal doping fits better with the marginal Fermi liquid phenomenology. Most notable is the observation of the pseudogap state in the underdoped region above T{sub c}. As in other strongly correlated electrons systems, these unusual properties stem from the interplay between electronic, magnetic, lattice and orbital degrees of freedom. Understanding the microscopic process involved in these materials and the interaction of electrons with other entities is essential to understand the mechanism of high temperature superconductivity. Since the discovery of high-T{sub c} superconductivity in cuprates, angle-resolved photoemission spectroscopy (ARPES) has provided key experimental insights in revealing the electronic structure of high temperature superconductors. These include, among others, the earliest identification of dispersion and a large Fermi surface, an anisotropic superconducting gap suggestive of a d-wave order parameter, and an observation of the pseudogap in underdoped samples. In the mean time, this technique itself has experienced a dramatic improvement in its energy and momentum resolutions, leading to a series of new discoveries not thought possible only a decade ago. This revolution of the ARPES technique and its scientific impact result from dramatic advances in four essential components: instrumental resolution and efficiency, sample manipulation, high quality samples and well-matched scientific issues. The purpose of this treatise is to go through the prominent results obtained from ARPES on cuprate superconductors. Because there have been a number of recent reviews on the electronic structures of high-T{sub c} materials, we will mainly present the latest results not covered previously, with a special attention given on the electron-phonon interaction in cuprate superconductors. What has emerged is rich information about the anomalous electron-phonon interaction well beyond the traditional views of the subject. It exhibits strong doping, momentum and phonon symmetry dependence, and shows complex interplay with the strong electron-electron interaction in these materials. ARPES experiments have been instrumental in identifying the electronic structure, observing and detailing the electron-phonon mode coupling behavior, and mapping the doping evolution of the high-T{sub c} cuprates. The spectra evolve from the strongly coupled, polaronic spectra seen in underdoped cuprates to the Migdal-Eliashberg like spectra seen in the optimally and overdoped cuprates. In addition to the marked doping dependence, the cuprates exhibit pronounced anisotropy with direction in the Brillouin zone: sharp quasiparticles along the nodal direction that broaden significantly in the anti-nodal region of the underdoped cuprates, an anisotropic electron-phonon coupling vertex for particular modes identified in the optimal and overdoped compounds, and preferential scattering across the two parallel pieces of Fermi surface in the antinodal region for all doping levels. This also contributes to the pseudogap effect. To the extent that the Migdal-Eliashberg picture applies, the spectra of the cuprates bear resemblance to that seen in established strongly coupled electron-phonon superconductors such as Pb. On the other hand, the cuprates deviate from this conventional picture. In the underdoped regime, the carriers are best understood as small polarons in an antiferromagnetic, highly electron correlated background, while the doped compounds require an anisotropic electron-phonon vertex to detail the prominent mode coupling signatures in the superconducting state. Electronic vertex corrections to the electron-phonon coupling furthermore may enhance, and for certain phonons, determine, the anisotropy of the electron-phonon coupling. A consistent picture emerges of the cuprates, combining strong, anisotropic electron-phonon coupling, particular phonon modes that could give rise to such a coupling, and an electron-electron interaction modifying the el-ph vertex.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗