SEARCH · Engineering Papers
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.
A formalism for target fragment transport is presented with application to energy loss spectra in thin silicon devices. A nuclear data base is recommended that agrees well with the measurements of McNulty et al. using surface barrier detectors. High-energy events observed by McNulty et al., which are not predicted by intranuclear cascade models, are well represented by the present work.
A systematic analysis of the observed reduced widths obtained in relativistic heavy ion fragmentation reactions is used to develop a phenomenological parametrization of these data. The parametrization is simple, accurate, and completely general in applicability.
An optical potential abrasion-ablation collision model is used to calculate hadronic dissociation cross sections for one, two, and three nucleon removal for the first time for a 14.6A GeV 28Si beam fragmenting in aluminum, tin, and lead targets. These estimates are compared with recent semi-inclusive measurements. Significant differences between some calculated and measured semi-inclusive cross sections exist which cannot be resolved without measurements of the exclusive channel hadronic cross sections. Calculations for each exclusive reaction channel contributing to the semi-inclusive cross sections are presented and discussed.
The characteristics of the hadron-to-quark first-order phase transition differ depending on whether charge neutrality is locally or globally fulfilled. In 𝛽-equilibrated matter, these two possibilities correspond to the Maxwell and Gibbs constructions. Recently, we presented a new framework in which a continuously varying parameter allows one to describe a first-order phase transition in intermediate scenarios to the two extremes of fully local and fully global charge neutrality. In this work, we extend the previous framework to finite temperatures and out-of-𝛽 equilibrium conditions, making it available for simulations of core-collapse supernovae and binary neutron star mergers. We investigate its impact on key thermodynamic quantities across a range of baryon densities, temperatures, and electron fractions. We find that when matter is not in 𝛽 equilibrium, the pressure in the mixed phase is not constant even for the case of fully local charge neutrality. Moreover, we compute the thermal index using three different approaches, demonstrating that the finite-temperature extension of an equation of state using a constant thermal index can be ill defined when applied to the mixed phase.
The particles and energies commonly used for hadron therapy overlap the low end of the charge and energy range of greatest interest for space radiation applications, Z=1-26 and approximately 100-1000 MeV/nucleon. It has been known for some time that the nuclear interactions of the incident ions must be taken into account both in treatment planning and in understanding and addressing the effects of galactic cosmic ray ions on humans in space. Until relatively recently, most of the studies of nuclear fragmentation and transport in matter were driven by the interests of the nuclear physics and later, the hadron therapy communities. However, the experimental and theoretical methods and the accelerator facilities developed for use in heavy ion nuclear physics are directly applicable to radiotherapy and space radiation studies. I will briefly review relevant data taken recently at various accelerators, and discuss the implications of the measurements for radiotherapy, radiobiology and space radiation research.
We investigate the impact of chemical equilibration and the resulting bulk viscosity on nonradial oscillation modes of warm neutron stars at temperatures up to T ≈ 5 MeV, relevant for protoneutron stars and neutron star postmerger remnants. In this regime, the relaxation rate of weak interactions becomes comparable to the characteristic frequencies of composition g-modes in the core, resulting in resonant damping. To capture this effect, we introduce the dynamical sound speed, a complex, frequency-dependent generalization of the adiabatic sound speed that encodes both the restoring force and the dissipative effects of bulk compression. Using realistic weak reaction rates and three representative equations of state, we compute the complex frequencies of composition g-modes with finite-temperature profiles. We find that bulk viscous damping becomes increasingly significant with temperature and can completely suppress composition g-modes. In contrast, the f-mode remains largely unaffected by bulk viscosity due to its nearly divergence-free character. Our results highlight the sensitivity of g-mode behavior to thermal structure, weak reaction rates, and the equation of state, and establish the dynamical sound speed as a valuable descriptor characterizing oscillation properties in dissipative neutron star matter.
This project advanced the understanding of how quarks and gluons produced in high-energy collisions transform into the hadrons observed in particle detectors, a fundamental process known as quantum chromodynamics (QCD) hadronization. By combining theoretical calculations, quantum simulation methods, and modern AI techniques, the research developed new tools to study multi-parton dynamics and nonperturbative effects that are essential for interpreting data from current and future nuclear physics experiments. Key outcomes include new theoretical frameworks for jet and hadron measurements, pioneering quantum simulation algorithms for real-time dynamics in field theories, and the development of advanced machine-learning models, such as diffusion models and explainable classifiers, to simulate and analyze collider events. These results are directly relevant to experiments at Jefferson Lab, Brookhaven National Laboratory, and the future Electron-Ion Collider, and they also have a broader impact in areas such as quantum information science and data-driven modeling of complex systems. The project supported the training of graduate students and postdoctoral fellows and contributed to the broader scientific community through publications, workshops, and collaborative activities. Overall, this work provides new insights into the microscopic mechanisms of hadron formation and establishes a foundation for future studies at the intersection of nuclear physics, artificial intelligence, and quantum computing.
The nuclear equation of state (EOS) is at the center of numerous theoretical and experimental efforts in nuclear physics. With advances in microscopic theories for nuclear interactions, the availability of experiments probing nuclear matter under conditions not reached before, endeavors to develop sophisticated and reliable transport simulations to interpret these experiments, and the advent of multi-messenger astronomy, the next decade will bring new opportunities for determining the nuclear matter EOS, elucidating its dependence on density, temperature, and isospin asymmetry. Among controlled terrestrial experiments, collisions of heavy nuclei at intermediate beam energies (from a few tens of MeV/nucleon to about 25 GeV/nucleon in the fixed-target frame) probe the widest ranges of baryon density and temperature, enabling studies of nuclear matter from a few tenths to about 5 times the nuclear saturation density and for temperatures from a few to well above a hundred MeV, respectively. Collisions of neutron-rich isotopes further bring the opportunity to probe effects due to the isospin asymmetry. However, capitalizing on the enormous scientific effort aimed at uncovering the dense nuclear matter EOS, both at RHIC and at FRIB as well as at other international facilities, depends on the continued development of state-of-the-art hadronic transport simulations. Furthermore, this white paper highlights the essential role that heavy-ion collision experiments and hadronic transport simulations play in understanding strong interactions in dense nuclear matter, with an emphasis on how these efforts can be used together with microscopic approaches and neutron star studies to uncover the nuclear EOS.
Using the Colour Glass Condensate description of electron-nucleus collisions at high energy, we study the diffractive production of a pair of jets with transverse momenta much larger than the nuclear saturation momentum Q s . At leading order in the QCD coupling, the di-jet cross-section exhibits transverse-momentum dependent (TMD) factorisation, with a gluon diffractive TMD distribution (DTMD) which is controlled by gluon saturation and describes the transverse-momentum imbalance between the produced jets. The next-to-leading corrections generate the various quantum evolutions of the diffractive gluon distribution. We focus on the Collins-Soper-Sterman (CSS) evolution which describes the change in the gluon DTMD when increasing the “hard scale” (the typical transverse momentum of the di-jets). We consider two different representations for this equation, one in transverse-momentum space, the other one in transverse-coordinate space. They are not fully equivalent with each other (despite being related by a Fourier transform) because of the respective boundary conditions. These conditions encode the essential physics of gluon saturation together with the effects of two other types of quantum evolution: the BK/JIMWLK evolution over the rapidity gap (“inside the Pomeron”) and the DGLAP evolution outside the rapidity gap (“within the diffractive system”). We demonstrate that, due to gluon saturation, one can compute both the boundary conditions and the CSS solutions mostly from first principles, without the need for a non-perturbative Sudakov. We numerically find a good agreement between the CSS solutions in the two aforementioned representations.
Hard-scattered partons ejected from high-energy proton-proton collisions undergo parton shower and hadronization, resulting in collimated collections of particles that are clustered into jets. A substructure observable that highlights the transition between the perturbative and nonperturbative regimes of jet evolution in terms of the angle between two particles is the two-point energy correlator (EEC). In this Letter, the first measurement of the EEC at RHIC is presented, using data taken from 200 GeV 𝑝 + 𝑝 collisions by the STAR experiment. The EEC is measured both for all the pairs of particles in jets and separately for pairs with like and opposite electric charges. These measurements demonstrate that the transition between perturbative and nonperturbative effects occurs within an angular region that is consistent with expectations of a universal hadronization regime that scales with jet momentum for a given initiator flavor. Additionally, a deviation from Monte Carlo predictions at small angles in the charge-selected sample could result from mechanics of hadronization not fully captured by current models.
Processes that violate baryon number, most notably proton decay and $n\bar{n}$ transitions, are promising probes of physics beyond the Standard Model (BSM) needed to understand the lack of antimatter in the Universe. To interpret current and forthcoming experimental limits, theory input from nuclear matrix elements to UV complete models enters. Thus, an interplay of experiment, effective field theory, lattice QCD, and BSM model building is required to develop strategies to accurately extract information from current and future data and maximize the impact and sensitivity of next-generation experiments. Here, we briefly summarize the main results and discussions from the workshop ‘INT-25-91W: Baryon Number Violation: From Nuclear Matrix Elements to BSM Physics,’ held at the Institute for Nuclear Theory, University of Washington, Seattle, WA, 13–17 January 2025.
The kagome lattice stands as a rich platform for hosting a wide array of correlated quantum phenomena, ranging from charge density waves and superconductivity to electron nematicity and loop current states. Direct detection of loop currents in kagome systems has remained a formidable challenge due to their intricate spatial arrangements and the weak magnetic field signatures they produce, and this has made their identification experimentally subtle. This has left their existence and underlying mechanisms a topic of intense debate. In this work, we uncover signatures compatible with loop currents: spin handedness-selective signals that surpass conventional dichroic, spin, and spin-dichroic responses. We observe this phenomenon in the kagome metal CsTi 3 Bi 5 and we call it the anomalous spin-optical helical effect. This effect arises from the coupling of light's helicity with spin-orbital electron correlations, thereby providing an indirect yet sensitive approach to probe loop-current–related electronic correlations in quantum materials. Our discovery not only enriches the debate surrounding loop currents but also offers new experimental strategies to exploit the electronic phases of quantum materials via light–matter interaction.
Antideuteron-nucleus absorptive cross sections for intermediate to high energies are calculated using an ion-ion optical model. Good agreement with experiment (within 15 percent) is obtained in this same model for mean p-nucleus cross sections at laboratory energies up to 15 GeV. We describe a technique for estimating antinucleus-nucleus cross sections from N mean N data and suggest that further cosmic ray studies to search for antideuterons and other antinuclei be undertaken.
We propose an adaptation of Entanglement Renormalization for quantum field theories that, through the use of discrete wavelet transforms, strongly parallels the tensor network architecture of the Multiscale Entanglement Renormalization Ansatz (a.k.a. MERA). Our approach, called wMERA, has several advantages of over previous attempts to adapt MERA to continuum systems. In particular, (i) wMERA is formulated directly in position space, hence preserving the quasi-locality and sparsity of entanglers; and (ii) it enables a built-in RG flow in the implementation of real-time evolution and in computations of correlation functions, which is key for efficient numerical implementations. As examples, we describe in detail two concrete implementations of our wMERA algorithm for free scalar and fermionic theories in (1+1) spacetime dimensions. Possible avenues for constructing wMERAs for interacting field theories are also discussed.
A comprehensive program of verification and validation has been undertaken to assess the applicability of models to space radiation shielding applications and to track progress as models are developed over time. The models are placed under configuration control, and automated validation tests are used so that comparisons can readily be made as models are improved. Though direct comparisons between theoretical results and experimental data are desired for validation purposes, such comparisons are not always possible due to lack of data. In this work, two uncertainty metrics are introduced that are suitable for validating theoretical models against sparse experimental databases. The nuclear physics models, NUCFRG2 and QMSFRG, are compared to an experimental database consisting of over 3600 experimental cross sections to demonstrate the applicability of the metrics. A cumulative uncertainty metric is applied to the question of overall model accuracy, while a metric based on the median uncertainty is used to analyze the models from the perspective of model development by analyzing subsets of the model parameter space.
The internal structure of protons and neutrons provides insight into both the dynamical behavior of the constitute quarks and gluons, and emergent properties of the nucleons (such as mass, spin, and electromagnetic distributions). Elastic electron-nucleon scattering can probe the elastic electromagnetic form factors of the nucleon. The electric and magnetic form factors, respectively, encode information about the internal charge and magnetization distributions within the nucleon. Precision data for these form factors, over a broad range of the four-momentum transfer squared, Q^2, can benchmark theoretical models describing the strong interaction of nuclear physics. The Super BigBite Spectrometer (SBS) program in Hall A at Jefferson Lab, is a series of high-precision experiments which seek to significantly extend the Q^2 reach of previous data for the nucleon electromagnetic form factors. The first two experiments of this program are known as G_M^n and the neutron Two Photon Exchange (nTPE) and the data were collected from October 2021 to February 2022. Both experiments were conducted with the simultaneous measurement of D(e,e'n) and D(e,e'p) reactions for quasi-elastic electron-deuteron scattering. The scattered electrons were detected in the BigBite Spectrometer, which features multiple large-acceptance Gas Electron Multiplier (GEM) detectors. The Super BigBite Spectrometer provided simultaneous detection of scattered nucleons, and utilized a large acceptance dipole magnet and Hadron Calorimeter (HCal). The G_M^n experiment provides precision measurements of the neutron magnetic form factor, via the ratio method, over a Q^2 range of 3.0 to 13.5 (GeV/c)2. From this data analysis, preliminary values for G_M^n/µ_n G_D are extracted. For Q^2=4.48 (GeV/c)2 we find G_M^n/µ_n G_D=0.9546±0.0132 and for Q^2=4.476 (GeV/c)2 we find G_M^n/µ_n G_D=0.9563±0.0110. These preliminary G_M^n/µ_n G_D values are more precise than existing world data in this Q^2 regime and are consistent with the most recent parameterization of the G_M^n/µ_n G_D world data. The nTPE experiment provides a first measurement of the neutron Rosenbluth Slope and seeks to quantify the two-photon exchange(TPE) contribution to elastic electron-neutron scattering at a fixed Q^2=4.5 (GeV/c)2 with two different beam energies and scattering angle values. For data of the proton form factor ratio, µ_p G_E^p/G_M^p, significant discrepancies exist between values obtained from Rosenbluth Separation and polarization transfer measurement, particularly at large Q^2, and TPE contributions are thought to resolve this discrepancy. The impacts of TPE contributions have not yet been experimentally established for the neutron. From the data analysis presented in this dissertation, a preliminary result for the neutron Rosenbluth Slope is found as S^n=(G_E^n )^2/t_n (G_M^n )^2=0.0916±0.0476 for Q^2=4.48 (GeV/c)2. This value of the neutron Rosenbluth Slope is consistent with the world data extrapolation and the absence of large TPE corrections.
Understanding the equation of state (EOS) of pure neutron matter is necessary for interpreting multimessenger observations of neutron stars. Reliable data analyses of these observations require well-quantified uncertainties for the EOS input, ideally propagating uncertainties from nuclear interactions directly to the EOS. This, however, requires calculations of the EOS for a prohibitively larger number of nuclear Hamiltonians, solving the nuclear many-body problem for each one. Quantum Monte Carlo methods, such as auxiliary-field diffusion Monte Carlo (AFDMC), provide precise and accurate results for the neutron matter EOS, but they are very computationally expensive, making them unsuitable for the fast evaluations necessary for uncertainty propagation. Here, we employ parametric matrix models to develop fast emulators for AFDMC calculations of neutron matter and use them to directly propagate uncertainties of coupling constants in the Hamiltonian to the EOS. As these uncertainties include estimates of the effective field theory truncation uncertainty, this approach provides robust uncertainty estimates for use in astrophysical data analyses. In conclusion, this Letter will enable novel applications such as using astrophysical observations to put constraints on coupling constants for nuclear interactions.
The production of protons in heavy ion collisions through the knockout mechanism (abrasion) is described using the Glauber model. The multiple knockouts from the projectile, including the inelastic collision series with the target, are considered using a closure approximation in treating energy conservation. Calculations for reactions of 12C projectiles with several targets at energies of 1 and 2A GeV arc compared to experiments. For large secondary proton momentum a strong dependence on the target mass is found and attributed to multiple scattering of the projectile knockouts.