Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “QUASI PARTICLES”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Quasiparticle description of angle-resolved photoemission spectroscopy for SrCuO 2

SrCuO 2 has long been considered a near-archetypal realization of a quasi-one-dimensional (1D) system of interacting electrons with short-range interactions. Within this framework, experimental observations - interpreted through the lens of the 1D Hubbard model - suggest that electron and hole excitations decay into two types of (unphysical) collective bosonic modes: "spinons", which carry the spin degree of freedom, and "holons", which carry the charge degree of freedom. This model, known as spin-charge separation, is most directly evidenced by angle-resolved photoemission spectroscopy (ARPES), where a photo-induced hole decays into a continuum of these excitations. Here we present an alternative perspective grounded in first-principles, self-consistent, and parameter-free many-body perturbation theory. In this revised quasiparticle framework, ARPES can be understood as a one-body effect arising from mild disorder in a long-range antiferromagnetic ground state. the emergence of the so-called spinon branch arises naturally from spin disorder, the anomalous linewidths are accurately captured, and we provide a compelling explanation for the spectral weight observed at the non-magnetic zone boundary. This reinterpretation provides a unified explanation for key experimental signatures previously attributed to spin-charge separation, including features observed in optical conductivity. Additionally, we show that SrCuO 2 exhibits a nontrivial interchain coupling that significantly influences both its one-particle and two-particle spectral functions. By comparing the spectral features of SrCuO 2 with those of La 2 CuO 4 , we argue that SrCuO 2 shares notable similarities with the two-dimensional cuprates-both being rooted in a common CuO 4 plaquette-based molecular orbital framework.

1D cuprates↗
DOE OSTI2023

Chiral photocurrent in a Quasi-1D TiS 3 (001) phototransistor

The presence of in-plane chiral effects, hence spin–orbit coupling, is evident in the changes in the photocurrent produced in a TiS 3 (001) field-effect phototransistor with left versus right circularly polarized light. The direction of the photocurrent is protected by the presence of strong spin–orbit coupling and the anisotropy of the band structure as indicated in NanoARPES measurements. As hsown here, dark electronic transport measurements indicate that TiS 3 is n-type and has an electron mobility in the range of 1–6 cm 2 V –1 s –1 . I–V measurements under laser illumination indicate the photocurrent exhibits a bias directionality dependence, reminiscent of bipolar spin diode behavior. Because the TiS 3 contains no heavy elements, the presence of spin–orbit coupling must be attributed to the observed loss of inversion symmetry at the TiS 3 (001) surface.

Phototransistor Materials↗
NASA NTRS1995

Particle signatures of magnetic topology at the magnetopause: AMPTE/CCE observations

Electron distributions at energies above 50 eV have been found to be a sensitive indicator of magnetic topology for magnetopause crossings of the AMPTE/CCE spacecraft. Progressing from the magnetosheath to the magnetosphere two abrupt transitions occur. First, the magnetosheath electron population directed either parallel or antiparallel to the magnetic field is replaced by a streaming, heated magnetosheath electron population. The other half of the distribution is unchanged. The region with unidirectional, heated magnetosheath electrons is identified as the magnetosheath boundary layer (MSBL). Second, the unheated magnetosheath electron population is replaced by a heated population nearly identical to the population encountered in the MSBL, resulting in a symmetric counterstreaming distribution. The region populated by the bidirectional heated magnetosheath electrons is identified as the low-latitude boundary layer (LLBL). The MSBL and LLBL identified by the electron transitions are the same as the regions identified using ion composition measurements. The magnetosheath-MSBL transition reflects a change in magnetic topology from a solar wind field line to one that threads the magnetopause, and the existence of a magnetosheath-MSBL transition implies that the magnetopause is open. When the current layer is easily identified, the MSBL-LLBL transition coincides with the magnetopause current layer, indicating that the magnetosheath electrons are heated in the current layer. Both magnetosheath-MSBL and MSBL-LLBL transitions are observed for low as well as high magnetic shears. Moreover, the transitions are particularly clear for low shear implying that magnetic topology boundaries are sharp even when abrupt changes in the field and other plasma parameters are absent. Furthermore, for low magnetic shear, solar wind ions with low parallel drift speeds make up the majority of the LLBL population indicating that the magnetosheath plasma has convected directly across the magnetosheath plasma has converted directly across the magnetopause. These observations are consistent with quasi-steady, high-latitude reconnection and indicate that the signatures of this reconnection geometry are commonly present in the subpolar region.

Fuselier, S. A.↗
DOE OSTI2024

Winos from natural SUSY at the high luminosity LHC

In natural supersymmetric models defined by no worse than a part in thirty electroweak fine-tuning, winos and binos are generically expected to be much heavier than Higgsinos. Moreover, the splitting between the Higgsinos is expected to be small, so that the visible decay products of the heavier Higgsinos are soft, rendering the Higgsinos quasi-invisible at the LHC. Within the natural supersymmetry (SUSY) framework, heavy electroweak gauginos decay to W , Z or h bosons plus Higgsinos in the ratio ∼ 2 ∶ 1 ∶ 1 , respectively. This is in sharp contrast to models with a binolike lightest superpartner and very heavy Higgsinos, where the charged (neutral) wino essentially always decays to a W ( h ) boson and an invisible bino. Wino pair production at the LHC, in natural SUSY, thus leads to V V , V h and h h + E T final states ( V = W , Z ) where, for TeV scale winos, the vector bosons and h daughters are considerably boosted. We identify eight different channels arising from the leptonic and hadronic decays of the vector bosons and the decay h → b b ¯ , each of which offers an avenue for wino discovery at the high luminosity LHC (HL-LHC). By combining the signal in all eight channels we find, assuming s = 14 TeV and an integrated luminosity of 3000 fb − 1 , that the discovery reach for winos extends to m ( w i n o ) ∼ 1.1 TeV , while the 95% CL exclusion range extends to a wino mass of almost 1.4 TeV. We also identify “Higgsino specific channels” which could serve to provide 3 σ evidence that winos lighter than 1.2 TeV decay to light Higgsinos rather than to a binolike lightest supersymmetric particle, should a wino signal appear at the HL-LHC. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2021

Energy Dependence in the Neutrino Scattering Data of MINERvA

As we prepare for lower neutrino energy beams where physics is highly energy-dependent, it is important to isolate factors that can contribute largely to these low energy neutrino experiments. The world of neutrino oscillation experiments uses a wide range of energies: from 0.7 GeV beam of T2K and MicroBooNE, 2 GeV beam of NOvA, to the 0.6 to 6.0 GeV beam produced by DUNE in the future. MINERvA’s data are currently the best place to test the upper end of the range for DUNE and can be extrapolated into both DUNE’s and NOvA’s oscillation maximum. The goal of the research is to constrain energy dependence in the neutrino scattering experiment using data from MINERvA and then determine what factors contribute to the observed energy dependence. The research is divided into two parts. First, I will analyze different theoretical models in different channels of interaction, such as the quasi-elastic and delta resonance channels, in terms of neutrino energy dependence. The goal is to study the models at the level of the structure functions, which emphasizes the W2 structure function (the C function for QE) dominates the MINERvA data, but the W3 is increasingly important for NOvA and DUNE. The second half of the research is data-driven and focuses on the different detectors and systematic effects. An experimental effect, the angle acceptance, is larger than the structure function effects. However, it is a detector geometry effect and is well measured and well modeled. Other effects such as the muon energy scale are small and localized. Additionally, the research will determine if the MINERvA GENIE model correctly predicts all of the energy dependence seen in the data, and identify the remaining unmodeled energy dependence between the MINERvA data and its best simulation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Search for a bound di-neutron by comparing 3 He(e,e'p)d and 3 H(e,e'p)X measurements

We report on a search for a bound di-neutron by comparing electron-induced proton-knockout $(e,e'p)$ measurements from Helium-3 ($^3$He) and Tritium ($^3$H). The measurements were performed at Jefferson Lab Hall A with a 4.326 GeV electron beam, and kinematics of large momentum transfer $Q^2 \approx 1.9$ (GeV/$c$)$^2$ and $x_B>1$, to minimize contributions from non quasi-elastic (QE) reaction mechanisms. Analyzing the measured $^3$He missing mass ($M_{miss}$) and missing energy ($E_{miss}$) distributions, we can distinguish the two-body break-up reaction, in which the residual proton-neutron system remains bound as a deuteron. In the $^3$H mirror case, under the exact same kinematic conditions, we do not identify a signature for a bound di-neutron with similar binding energy to that of the deuteron. We calculate exclusion limits as a function of the di-neutron binding energy and find that, for binding equivalent to the deuteron, the two-body break-up cross section on $^3$H is less than 0.9% of that on $^3$He in the measured kinematics at the 95% confidence level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Kinetic simulations comparing quasi-parallel and quasi-perpendicular piston-driven collisionless shock dynamics in magnetized laboratory plasmas

Magnetized collisionless shocks are common in astrophysical systems, and scaled versions can be created in laboratory experiments by utilizing laser-driven piston plasmas to create these shocks in a magnetized background plasma. A key parameter for these experiments is the angle θB between the shock propagation direction and the background magnetic field. We performed quasi-1D piston-driven shock simulations to explore shock formation, evolution, and key observables relevant to laboratory experiments for a range of shock angles between θB=90° to θB=30°. Our results show that the spatial and temporal scales of shock formation for all angles considered are similar when expressed in terms of the perpendicular component of the magnetic field. In a steady state, ion and electron temperatures become more isotropic, and the electron-to-ion temperature ratio is higher for smaller θB. At θB=30°, ion heating parallel to the magnetic field becomes dominant, associated with more ions being reflected at one discontinuity and subsequently trapped by the next discontinuity due to shock reformation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Neutrino-nucleus CC 0 π cross-section tuning in GENIE v3

This article summarizes the state of the art of v μ and $\bar{v}_μ$ CC0π cross-section measurements on carbon and argon and discusses the relevant nuclear models, parametrizations and uncertainties in GENIE v3. The CC0π event topology is common in experiments at a few-GeV energy range. Although its main contribution comes from quasi-elastic interactions, this topology is still not well understood. The GENIE global analysis framework is exploited to analyze CC0π datasets from MiniBooNE, T2K and MINERvA. A partial tune for each experiment is performed, providing a common base for the discussion of tensions between datasets. The results offer an improved description of nuclear CC0π datasets as well as data-driven uncertainties for each experiment. This work is a step towards a GENIE global tune that improves our understanding of neutrino interactions on nuclei. It follows from earlier GENIE work on the analysis of neutrino scattering datasets on hydrogen and deuterium.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2016

First Flight of the Gamma-Ray Imager Polarimeter for Solar Flares (GRIPS) Instrument

The Gamma-Ray Imager/Polarimeter for Solar ares (GRIPS) instrument is a balloon-borne telescope designed to study solar-flare particle acceleration and transport. We describe GRIPS's first Antarctic long-duration flight in January 2016 and report preliminary calibration and science results. Electron and ion dynamics, particle abundances and the ambient plasma conditions in solar flares can be understood by examining hard X-ray (HXR) and gamma-ray emission (20 keV to 10 MeV). Enhanced imaging, spectroscopy and polarimetry of flare emissions in this energy range are needed to study particle acceleration and transport questions. The GRIPS instrument is specifically designed to answer questions including: What causes the spatial separation between energetic electrons producing hard X-rays and energetic ions producing gamma-ray lines? How anisotropic are the relativistic electrons, and why can they dominate in the corona? How do the compositions of accelerated and ambient material vary with space and time, and why? GRIPS's key technological improvements over the current solar state of the art at HXR/gamma-ray energies, the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI), include 3D position-sensitive germanium detectors (3D-GeDs) and a single-grid modulation collimator, the multi-pitch rotating modulator (MPRM). The 3D-GeDs have spectral FWHM resolution of a few hundred keV and spatial resolution less than 1cu mm. For photons that Compton scatter, usually greater or equal to 150 keV, the energy deposition sites can be tracked, providing polarization measurements as well as enhanced background reduction through Compton imaging. Each of GRIPS's detectors has 298 electrode strips read out with ASIC/FPGA electronics. In GRIPS's energy range, indirect imaging methods provide higher resolution than focusing optics or Compton imaging techniques. The MPRM grid-imaging system has a single-grid design which provides twice the throughput of a bi-grid imaging system like RHESSI. The grid is composed of 2.5 cm deep tungsten-copper slats, and quasi-continuous FWHM angular coverage from 12.5-162 arcsecs are achieved by varying the slit pitch between 1-13 mm. This angular resolution is capable of imaging the separate magnetic loop footpoint emissions in a variety of are sizes. In comparison, RHESSI's 35-arcsec resolution at similar energies makes the footpoints resolvable in only the largest ares.

Duncan, Nicole↗
DOE OSTI2025

Toward engineering lattice structures with the material point method (MPM)

This study examines the potential of two variants of the material point method—the generalized interpolation material point (GIMP) and dual domain material point (DDMP) methods—in developing a robust computational framework for engineering lattice structures under different loading conditions. The study begins with assessing the ability of the two methods in predicting elastic buckling phenomena using column geometries with and without initial geometric imperfections. The results indicate that both methods effectively capture buckling phenomena when initial geometric imperfections are introduced. After this verification step, we create several models of tetrahedral lattice structures with varying strut diameter and orientation and subject them to quasi-static loading. We then validate the numerical results using laboratory test results. The results show that, while both methods accurately predict load–displacement curves in the pre-buckling regime, their predictive capabilities diminish in the post-buckling regime. Through visual comparison between the numerical and experimental deformed shapes, it appears that the discrepancies between model and experimental results are attributed to initial geometric imperfections in the lattices that occurred during 3D printing. We then establish a second set of lattice models where different types of initial geometric imperfections are considered. The results from these models show that imperfections have a negligible influence in the pre-buckling regime but affect the behavior considerably in the post-buckling regime. As a final step in this work, we subject the lattice models to impact loading and employ hypothetical soft and stiff materials. These results show that the lattice stiffness, which depends on material stiffness, strut diameter, and orientation, significantly influences the ability of a lattice structure to resist impact. In particular, we find that a stiffer lattice (i.e., one made with a stiff material and thicker struts) is capable of absorbing more energy than a softer one during impact. Although material nonlinearities, inelasticity, and detailed contact formulations are not considered in this study, the findings obtained herein lay the groundwork for engineering lattice structures under extreme loading conditions through a simulation-driven framework based on particle-based methods.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2000

Numerical Model of the Plasma Sheath Generated by the Plasma Source Instrument Aboard the Polar Satellite

The plasma sheath generated by the operation of the Plasma Source Instrument (PSI) aboard the POLAR satellite is studied by using a 3-dimensional Particle-In-Cell (PIC) code. When the satellite passes through the region of low density plasma, the satellite charges to positive potentials as high as 4050Volts due to the photoelectrons emission. In such a case, ambient core ions cannot accurately be measured or detected. The goal of the onboard PSI is to reduce the floating potential of the satellite to a sufficiently low value so that the ions in the polar wind become detectable. When the PSI is operated, an ion-rich Xenon plasma is ejected from the satellite, such that the floating potential of the satellite is reduced and is maintained at about 2Volts. Accordingly, in our 3-dimensional PIC simulation, we considered that the potential of the satellite is 2Volts as a fixed bias. Considering the relatively high density of the Xenon plasma in the sheath (approx. 10 - 10(exp 3)/cc), the ambient plasma of low density (less than 1/cc) is neglected. In the simulations, the electric fields and plasma dynamics are calculated self-consistently. We found that an "Apple" shape positive potential sheath forms surrounding the satellite. In the region near the PSI emission, a high positive potential hill develops. Near the Thermal Ion Detection Experiment (TIDE) detector away from the PSI, the potentials are sufficiently low for the ambient polar wind ions to reach it. In the simulations, it takes about a hundred electron gyroperiods for the sheath to reach a quasi-steady state. This time is approximately the time taken by the heavy Xe(+) ions to expand up to about one average Larmor radius of electrons from the satellite surface. Using the steady state sheath, we performed trajectory calculations to characterize the detector response to a highly supersonic polar wind flow. The detected ions' velocity distribution shows significant deviations from a shifted Maxwellian in the ambient polar wind population. The deviations are caused by the effects of electric fields on the ions' motion as they traverse the sheath.

Leung, Wing C.↗
NASA NTRS2001

Numerical Model of the Plasma Sheath Generated by the Plasma Source Instrument Aboard the Polar Satellite

The plasma sheath generated by the operation of the Plasma Source Instrument (PSI) aboard the Polar satellite is studied by using a three-dimensional particle-in-cell (PIC) code. When the satellite passes through the region of low-density plasma, the satellite charges to positive potentials as high as 40-50 V, owing to the photoelectron emission. In such a case, ambient core ions cannot accurately be measured or detected. The goal of the onboard PSI is to reduce the floating potential of the satellite to a sufficiently low value so that the ions in the polar wind become detectable. When the PSI is operated, ion-rich xenon plasma is ejected from the satellite, such that the floating potential of the satellite is reduced and is maintained at approximately 2 V. Accordingly, in our three-dimensional PIC simulation we considered that the potential of the satellite is 2 V as a fixed bias. Considering the relatively high density of the xenon plasma in the sheath (10-10(exp 3)/cc), the ambient plasma of low density (<1/cc) is neglected. In the simulations the electric fields and plasma dynamics are calculated self-consistently. We found that an 'apple'-shape positive potential sheath forms surrounding the satellite. In the region near the PSI emission a high positive potential hill develops. Near the Thermal Ion Dynamics Experiment detector away from the PSI, the potentials are sufficiently low for the ambient polar wind ions to reach it. In the simulations it takes only about a couple of tens of electron gyroperiods for the sheath to reach a quasi steady state. This time is approximately the time taken by the heavy Xe(+) ions to expand up to about one average Larmor radius of electrons from the satellite surface. After this time the expansion of the sheath in directions transverse to the ambient magnetic field slows down because the electrons are magnetized. Using the quasi steady sheath, we performed trajectory calculations to characterize the detector response to a highly supersonic polar wind flow. The detected ions' velocity distribution shows significant deviations from a shifted Maxwellian in the ambient polar wind population. The deviations are caused by the effects of electric fields on the ions' motion as they traverse the sheath.

Singh, N.↗
NASA NTRS2003

Properties of the Hermean Regolith: New Optical Reflectance Spectra, Comparison with Lunar Anorthosites and Mineralogical Modelling - 5

We present new optical (0.4-0.65 micron) spectra of Mercury and lunar pure anorthosite locations, obtained quasi-simultaneously with the Nordic Optical Telescope (NOT) in 2002. A comparative study is performed with the model of Lucey et al. between iron-poor, mature, pure anorthosite (less than 90% plagioclase feldspar) Clementine spectra from the lunar farside and a combined 0.4-1.0 micron mercurian spectrum, obtained with the NOT, calculated for standard photometric geometry. Mercury is located at more extreme locations in the Lucey ratio-reflectance diagrams than any known lunar soil, specifically with respect to the extremely iron-poor mature anorthosites. Though quantitative prediction of FeO and TiO2 abundances cannot be made without a more generally applicable model, we find qualitatively that the abundances of both these oxides must be near zero for Mercury. We utilize the theory of Hapke, with realistic photometric parameters, to model laboratory spectra of matured mineral powders and lunar soils, and remotely sensed spectra of lunar anorthosites and Mercury. An important difference between fabricated and natural powders is the high value for the internal scattering parameter necessary to interpret the spectra for the former, and the requirement of rough and non-isotropically scattering surfaces in the modelling of the latter. The mature lunar anorthosite spectra were well modelled with binary mixtures of calcic feldspars and olivines, grain sizes of 25-30 micron and a concentration of submicroscopic metallic iron (SMFe) of 0.12-0.15% in grain coatings. The mercurian spectrum is not possible to interpret from terrestrial mineral powder spectra without introducing an average particle scattering function for the bulk soil that increases in backscattering efficiency with wavelength. The observed spectrum is somewhat better predicted with binary mixture models of feldspars and pyroxenes, that single-component regoliths consisting of either albite or diopside. Correct spectral reflectance values were predicted with a concentration of 0.1 wt% SMFe in coatings of 15-30 micron sized grains. Since reasonable cosmogonical formation scenarios for Mercury, or meteoritic infall, predict iron concentrations at least this high, we draw the conclusion that the average grain size of Mercury is about a factor of two smaller than for average returned lunar soil samples. The 0.6-2.5 micron spectrum of McCord and Clark is used to further limit the possible range of mineralogical composition of Mercury. It is found that an intimately mixed and matured 3 : 1 labradorite-to-enstatite regolith composition best matches both the optical and near-infrared spectra, yielding an abundance of approx. 1.2 wt% FeO and -approx. 0 wt% TiO2.

Warell, J.↗
DOE OSTI2025

SULI Intern Final Report: Computationally Investigating Hydrogen Thermo-Diffusion in Yttrium Hydride Using Multiscale Methods

The renaissance of nuclear energy has arrived, heralding an age of abundant inexpensive clean energy, and renewed space exploration. In nuclear-powered spacecraft and microreactors, safety and size are of utmost importance. Yttrium Hydride (YHx) is being researched for its utility as a neutron moderator in nuclear reactors; the hydrogen in YHx slows down neutrons, enabling a continuous nuclear reaction in the reactor. This has the benefit of allowing reactors to be more safe, compact, and efficient. The goal of this effort is to computationally predict the coefficient of temperature-dependent hydrogen diffusion within YHx, the Soret coefficient. This parameter is essential for determining the safe operating modes of YHx moderators. Zirconium Hydride (ZrHx) is used in the Training, Research, Isotopes, General Atomics (TRIGA) reactor, is the reference material for these calculations. In this work, nanoscale atomic modeling in the Vienna Ab initio Simulation Package (VASP) is combined with the mesoscale finite element phase-field module in the Multiphysics Object-Oriented Simulation Environment (MOOSE); this culminates in a new multiscale computational method to simulate Soret diffusion of hydrogen in YHx. This data is useful for predicting experimental outcomes. This workflow involves convergence testing followed by static, Nudged Elastic Band (NEB), Quasi-Harmonic Approximation (QHA), and Molecular Dynamics (MD) calculations - linked with phase field simulation. NEB simulates hydrogen migration, while QHA and MD predict temperature-dependent properties. The static calculations align with literature, and preliminary NEB and QHA calculations yield accurate results. Once the atomic calculations are complete, we will incorporate Electron Backscatter Diffraction (EBSD) images and VASP-generated parameters into the phase field module to simulate intra- and intergranular transport of hydrogen in ZrHx and YHx. Future research will extend our approach to fuel-moderator materials systems such as Uranium-Yttrium Hydride (U-YHx). This work contributes to the development of advanced nuclear energy solutions for space travel.

36 - MATERIALS SCIENCE↗
DOE OSTI2024

A prototype scintillator real‐time beam monitor for ultra‐high dose rate radiotherapy

Background: FLASH Radiotherapy (RT) is an emergent cancer RT modality where an entire therapeutic dose is delivered at more than 1000 times higher dose rate than conventional RT. For clinical trials to be conducted safely, a precise and fast beam monitor that can generate out-of-tolerance beam interrupts is required. This paper describes the overall concept and provides results from a prototype ultra-fast, scintillator-based beam monitor for both proton and electron beam FLASH applications. Purpose: A FLASH Beam Scintillator Monitor (FBSM) is being developed that employs a novel proprietary scintillator material. The FBSM has capabilities that conventional RT detector technologies are unable to simultaneously provide: (1) large area coverage; (2) a low mass profile; (3) a linear response over a broad dynamic range; (4) radiation hardness; (5) real-time analysis to provide an IEC-compliant fast beam-interrupt signal based on true two-dimensional beam imaging, radiation dosimetry and excellent spatial resolution. Methods: The FBSM uses a proprietary low mass, less than 0.5 mm water equivalent, non-hygroscopic, radiation tolerant scintillator material (designated HM: hybrid material) that is viewed by high frame rate CMOS cameras. Folded optics using mirrors enable a thin monitor profile of ∼10 cm. A field programmable gate array (FPGA) data acquisition system generates real-time analysis on a time scale appropriate to the FLASH RT beam modality: 100–1000 Hz for pulsed electrons and 10–20 kHz for quasi-continuous scanning proton pencil beams. An ion beam monitor served as the initial development platform for this work and was tested in low energy heavy-ion beams ( 86 Kr +26 and protons). A prototype FBSM was fabricated and then tested in various radiation beams that included FLASH level dose per pulse electron beams, and a hospital RT clinic with electron beams. Results: Results presented in this report include image quality, response linearity, radiation hardness, spatial resolution, and real-time data processing. Furthermore, the HM scintillator was found to be highly radiation damage resistant. It exhibited a small 0.025%/kGy signal decrease from a 216 kGy cumulative dose resulting from continuous exposure for 15 min at a FLASH compatible dose rate of 237 Gy/s. Measurements of the signal amplitude versus beam fluence demonstrate linear response of the FBSM at FLASH compatible dose rates of >40 Gy/s. Comparison with commercial Gafchromic film indicates that the FBSM produces a high resolution 2D beam image and can reproduce a nearly identical beam profile, including primary beam tails. The spatial resolution was measured at 35–40 µm. Tests of the firmware beta version show successful operation at 20 000 Hz frame rate or 50 µs/frame, where the real-time analysis of the beam parameters is achieved in less than 1 µs. Conclusions: The FBSM is designed to provide real-time beam profile monitoring over a large active area without significantly degrading the beam quality. A prototype device has been staged in particle beams at currents of single particles up to FLASH level dose rates, using both continuous ion beams and pulsed electron beams. Using a novel scintillator, beam profiling has been demonstrated for currents extending from single particles to 10 nA currents. Radiation damage is minimal and even under FLASH conditions would require ≥50 kGy of accumulated exposure in a single spot to result in a 1% decrease in signal output. Beam imaging is comparable to radiochromic films, and provides immediate images without hours of processing. Real-time data processing, taking less than 50 µs (combined data transfer and analysis times), has been implemented in firmware for 20 kHz frame rates for continuous proton beams.

2D beam imaging↗
DOE OSTI2020

Comparing event generator predictions and ab-initio calculations of $\nu$-$^{12}$C neutral current quasi-elastic scattering at 1 GeV

The measurement of neutrino oscillations and exotic physics searches are important parts of the physics program in the near future, with new state-of-the-art experiments planned within the next decade. Future and modern experiments in these fields will make use of nuclear targets. Event Generators (EGs) are software used in the analysis of neutrino oscillation experiments. EGs use to predict kinematic observables for a range of neutrino energies. These simulations may lack physics captured by more rigorous theoretical calculations. This work compares EG performance to nuclear theory calculations by comparing observables generated in the two frameworks. We provide a common set of definitions between theory and experiment and assess the physics contained in EG simulations. Neutral current quasi-elastic (NCQE) scattering events for neutrinos and anti-neutrinos on a $^{12}$C target are simulated with a specific EG, NEUT, used by the T2K experiment for its analysis. The simulated cross sections are compared to analytic calculations from nuclear theory within the factorization scheme. We compare the NEUT implementation of two different models on nuclear spectral functions: the Relativistic Fermi Gas (RFG) and the correlated basis spectral function (CBF) to analytic calculations of the same models in the factorization scheme. For both nuclear physics models, we compare the appearance of features in the distributions relevant to experimental analyses. Qualitatively, the shape of the simulated distribution is similar to the one obtained through theory calculations; however, there are some discrepancies between the theory calculations and the NEUT simulation. While the EG simulations and analytic calculations with the same model of nuclear dynamics show similar overall features, there are still differences between the two.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Flux Ropes, Turbulence, and Collisionless Perpendicular Shock Waves: High Plasma Beta Case

With the onset of solar maximum and the expected increased prevalence of interplanetary shock waves, Parker Solar Probe is likely to observe numerous shocks in the next few years. An outstanding question that has received surprisingly little attention has been how turbulence interacts with collisionless shock waves. Turbulence in the supersonic solar wind is described frequently as a superposition of a majority 2D and a minority slab component. Here we formulate a collisional perpendicular shock-turbulence transmission problem in a way that enables investigation of the interaction and transmission of quasi-perpendicular fluctuations such as magnetic flux ropes/islands and vortices as well as entropy and acoustic modes in the large plasma beta regime. We focus on the transmission of an upstream spectrum of these modes, finding that the downstream spectral amplitude is typically increased significantly (a factor of 10 or more), and that the upstream spectral index of the inertial range, and indeed the general spectral shape, is unchanged for the downstream magnetic variance, kinetic energy, and density variance. A comparison of the theoretically predicted downstream magnetic variance, kinetic energy, and density variance spectra with those observed at 1, 5, and 84 au by Wind, Ulysses, and Voyager 2 shows excellent agreement. The overall theoretically predicted characteristics of the transmission of turbulence across shocks observed in the solar wind appear to be largely consistent with recent observational studies by Pitňa et al. and Borovsky.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2019

Cloud Computing Methods for Near Rectilinear Halo Orbit Trajectory Design

Complicated mission design problems require innovative computational solutions. As spacecraft depart from a proposed Gateway in a Near Rectilinear Halo Orbit (NRHO), recontact analysis is required to avoid risk of collision and ensure safe operations. Escape dynamics from NRHOs are governed by multiple gravitational bodies, yielding a trajectory design space that is exhaustively large. This paper summarizes the recontact analysis for departure from the NRHO and describes how the Deep Space Trajectory Explorer (DSTE) trajectory design software incorporates high performance cloud computing to compute and visualize the orbit design space. Recent focus on exploration missions to cislunar space has kindled accelerated interest in multibody orbit solutions. Trajectory analysis in the presence of multiple gravity fields is complex, and innovative computational tools are needed to simplify complicated design spaces, to generate large quantities of data quickly, and to visualize the output for user accessibility. The Gateway mission is a prime example. The Gateway1 is proposed as a human outpost in deep space. The current baseline orbit for the Gateway is a Near Rectilinear Halo Orbit (NRHO) near the Moon.2 The NRHO exists in a regime that experiences the gravitational effects of the Earth and the Moon simultaneously, complicating orbit analysis. The mission design process benefits greatly from updated computational tools for multibody missions like the Gateway. As an example, consider the problem of assessing the risk of collision in an NRHO. As a staging location to missions to the lunar surface and beyond the Earth-Moon system, the Gateway will experience spacecraft and other objects regularly arriving and departing. Departing objects potentially include spent logistics modules, visiting crew vehicles, debris objects, wastewater particles, and cubesats. Each departure is governed by the dynamics of the Gateway orbit and the surrounding dynamical environment. Over time, any unmaintained object in such an orbit eventually departs due to the small instabilities associated with the NRHOs. A separation maneuver speeds the departure from the NRHO, but the effects of the maneuver on the spacecraft behavior depend on the location, magnitude, and direction of the burn. Escape dynamics from the NRHO with regard to these maneuver options open up an enormous potential trajectory design space where subtle changes in input can produce dramatically large changes in the results. Any departing object must avoid recontacting the Gateway as it leaves the lunar vicinity, and a recontact analysis thus involves a significant number of computations and extensive output data. To explore the dynamics of this extensive design space, the Deep Space Trajectory Explorer3 (DSTE) trajectory design software incorporates new High Performance Computing (HPC) services and novel interactive visualizations. This paper details the HPC and cloud infrastructure techniques that are implemented in the DSTE, applying the new capabilities to analysis of recontact risk with the Gateway in NRHO. NEAR RECTILINEAR HALO ORBITS The Gateway is planned to fly in a lunar NRHO as its baseline orbit. The NRHO families of orbits are subsets of the larger halo families, which originate from planar orbits near the L1 and L2 libration points; the Earth-Moon L2 halo family appears in Figure 1. Each halo orbit is perfectly periodic in the Circular Restricted 3-Body Problem (CR3BP) and becomes a quasi-periodic orbit in a higher fidelity ephemeris force model. The NRHOs are defined as those members of the halo family with bounded stability properties;2 they pass near the Moon at perilune and are nearly polar. Families exist with apolunes located both above the lunar north pole and above the lunar south pole; the Gateway is planned to reside in a southern L2 NRHO in a 9:2 resonance with the lunar synodic period. The 9:2 NRHO is characterized by a period of about 6.5 days, a perilune radius of about 3,500 km, and an apolune radius of about 71,000 km; it is strongly affected by the gravity of both the Earth and the Moon simultaneously. This NRHO offers extended communications with assets on the south pole of the Moon,4 as well as low-cost orbit maintenance and attitude control,5 favorable eclipse avoidance properties,6 and inexpensive transfers from Earth and to other destinations.5,7 The NRHO portion of the southern L2 halo family is highlighted in black in Figure 1, and the 9:2 NRHO appears in blue.

Phillips, Sean M.↗