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.
Recent experiments have addressed the nature of the charge order seen in underdoped cuprates. In this paper, I show that x-ray absorption and linear dichroism are excellent probes of such order. Ab initio calculations reveal that a d-wave charge density wave order involving the oxygen ions is a much better description of the data than alternate models.
Zirconium hydride (ZrH x ) is a moderator material for TRIGA reactors and historical space reactor systems, such as SNAP-10A. Thermal neutron scattering laws (TSL) for two phases of this material, δ and ε, have been previously evaluated by Naval Nuclear Laboratory (NNL) and submitted to the National Nuclear Data Center (NNDC) for inclusion in the US national ENDF/B-VIII.1 nuclear data library. In contrast to the current ENDF/B-VIII.0 TSL evaluations, which consider only a single phase, the new evaluations are derived from separate ab initio calculations for both phases and include coherent elastic effects of the zirconium sublattice. To estimate the impact of these changes to the TSL evaluation of this material, comparative critical mass calculations were performed with MC21 for homogenous mixtures of high- enriched uranium (HEU) and ZrH x in bare and water reflected sphere configurations. These calculations yield an impact on the estimated critical mass as a function of 235 U loading density with maximum differences as large as 1% - 5% for over-moderated thermal spectrum systems. Consequently, the NNL TSL evaluations are anticipated to have a small impact on criticality calculations of thermal reactor systems regardless of the material phase. Nevertheless, characteristic differences exist in the predicted thermal spectra as function of energy for the two sets of TSL evaluations, which are attributed to difference in the underlying phonon density of states of hydrogen bound in ZrH x .
The interaction energies for N2-He and N2-H2 are calculated by accurate ab initio methods. The virial coefficient and differential scattering cross section for N2-H2 are calculated; the theoretical results are compared with experimental data. The transport collision integrals for N2-H2 and N2-N2 interactions are calculated and tabulated; the results yield transport coefficients that compare well with measured data. Transport coefficients are found to be determined accurately from the interaction energies for a specific configuration of the molecule formed from the interaction partners. Comparisons with results of measurement and accurate calculations demonstrate that the transport properties of complex molecular interactions can be determined rapidly and fairly accurately from the interaction energies of simpler system using combination rules for the short-range parameters of effective interaction energies and the coefficients for the long-range forces. The coefficients for a two-parameter temperature expansion of diffusion and viscosity are tabulated for a realistic universal potential energy that is based primarily on the results of very accurate calculations of the He-He interaction energy.
Nuclei are used for high-precision tests of the Standard Model and for studies of physics beyond the Standard Model. Without a thorough understanding of nuclei, we will not be able to meaningfully interpret the growing body of experimental data nor will we be able to disentangle new physics signals from underlying nuclear effects. This calls for accurate calculations of nuclear structure and reactions. In this work, we focus on electroweak decays in nuclei with mass number A ≤ 10 and report on ab initio Quantum Monte Carlo calculations of reduced matrix elements entering beta decays and electron captures in nuclei with mass number A ≤ 10. The many-body wave functions are calculated using selected Norfolk two- and three-nucleon potential models and associated one- and two-body axial currents at tree-level obtained from a chiral effective field theory with pions, nucleons, and Δ. The agreement with the experimental data is satisfactory except for transitions in A = 8 nuclei. In this specific case, the theory significantly underpredicts the experimental data, which indicates the need of further improvements in the corresponding nuclear wave functions. In this study, emphasis is placed on the contributions of two-body axial currents that are carefully analyzed using two-body transition densities. This allow us to study the spatial distribution and short-range behavior of two-body dynamics. In particular, the transition densities when scaled to peak at 1.0 exhibit universal short-range behavior across the considered nuclei, while they differ in the long-range tails.
We report a molecular dynamics study of ab initio quality of the ferroelectric phase transition in crystalline PbTiO3. We model anharmonicity accurately in terms of potential energy and polarization surfaces trained on density functional theory data with modern machine learning techniques. Our simulations demonstrate that the transition has a strong order-disorder character, in agreement with diffraction experiments, and provide fresh insight into the approach to equilibrium across the phase transition. We find that the emergence and disappearance of the macroscopic polarization is driven by dipolar switching at the nanometer scale. We also computed the infrared optical absorption spectra in both the ferroelectric and the paraelectric phases, finding good agreement with the experimental Raman frequencies. Often, the almost ideal displacive character of the soft mode detected by Raman scattering in the paraelectric phase has been contrasted with the order-disorder character of the transition suggested by diffraction experiments. We settle this issue by showing that the soft mode coexists with a strong Debye relaxation associated with thermal disordering of the dipoles. The Debye relaxation feature is centered at zero frequency and appears near the transition temperature in both the ferroelectric and the paraelectric phases.
Recent years have seen exciting new developments and progress in nuclear structure theory, reaction theory, and experimental techniques, that allow us to move towards a description of exotic systems and environments, setting the stage for new discoveries. The purpose of the 5-week program was to bring together physicists from the low-energy nuclear structure and reaction communities to identify avenues for achieving reliable and predictive descriptions of reactions involving nuclei across the isotopic chart. The 4-day embedded workshop focused on connecting theory developments to experimental advances and data needs for astrophysics and other applications. Nuclear theory must address phenomena from laboratory experiments to stellar environments, from stable nuclei to weakly-bound and exotic isotopes. Expanding the reach of theory to these regimes requires a comprehensive understanding of the reaction mechanisms involved as well as detailed knowledge of nuclear structure. A recurring theme throughout the program was the desire to produce reliable predictions rooted in either ab initio or microscopic approaches. At the same time it was recognized that some applications involving heavy nuclei away from stability, e.g. those involving fi ssion fragments, may need to rely on simple parameterizations of incomplete data for the foreseeable future. The goal here, however, is to subsequently improve and refine the descriptions, moving to phenomenological, then microscopic approaches. There was overarching consensus that future work should also focus on reliable estimates of errors in theoretical descriptions.
The focus of this work is on carbon clusters, hydrocarbons and other carbonaceous species which are constituents of pyrolysis gases injected into the boundary layer. Some of these molecules are observed to have absorption spectra in the VUV and UV region that match the emission spectra of atomic nitrogen and oxygen. Hence, they can potentially absorb the radiation impinging on the heat shield. Accurate thermodynamic data is not available for many of these species over the relevant temperature range for reentry. The objective of this work is to determine the thermochemical properties of potential radiation absorbing molecules using accurate ab initio quantum chemistry calculations. We have compiled a database of formation enthalpies, anharmonic vibrational frequencies and rotational constants for the ground and low-lying excited electronic states of these carbonaceous species. Using these thermochemical properties, we compute equilibrium mole fractions and observe that the mole fractions of some species of interest are significant, thus they could affect the radiative heat flux. Furthermore, comparisons of the formation enthalpies and mole fractions are made with data from literature.
Light emission under MeV hydrogen and oxygen ions in stoichiometric SrTiO 3 are identified at temperatures of 100 K, 170 K and room-temperature. MeV ions predominately deposit their energies to electrons in SrTiO 3 with energy densities orders of magnitude higher than from UV or x-ray sources but comparable to femtosecond lasers. The ionoluminescence (IL) spectra can be resolved into three main Gaussian bands at 2.0 eV, 2.5 eV and 2.8 eV, whose relative contributions strongly depend on irradiation temperature, electronic energy loss and irradiation fluence. Two main bands, observed at 2.5 eV and 2.8 eV, are intrinsic and associated with electron–hole recombination in the perfect SrTiO 3 lattice. The 2.8 eV band is attributed to recombination of free (conduction) electrons with an in-gap level, possibly related to self-trapped holes. Self-trapped excitons (STEs) are considered suitable candidates for the 2.5 eV emission band, which implies a large energy relaxation in comparison to the intrinsic edge transition. The dynamics of electronic excitation, governs a rapid initial rise of the intensity; whereas, accumulated irradiation damage (competing non-radiative recombination channels) accounts for a subsequent intensity decrease. The previously invoked role of isolated oxygen vacancies for the blue luminescence (2.8 eV) does not appear consistent with the data. An increasing well-resolved band at 2.0 eV dominates at 170 K and below. It has been only previously observed in heavily strained and amorphous SrTiO 3 , and is, here, attributed to transitions from d(t 2g) conduction band levels to d(e g) levels below the gap. In accordance with ab initio theoretical calculations they are associated to trapped electron states in relaxed Ti 3+ centers at an oxygen vacancy within distorted TiO 6 octahedra. The mechanism of defect evolution monitored during real-time IL experiments is presented. In conclusion, the light emission data confirm that IL is a useful tool to investigate lattice disorder in irradiated SrTiO 3 .
We examine the relationship between the asymptotic normalization coefficient (ANC) of 6 Li and other low-energy observables in the α–deuteron system. Our analysis uses a set of calculations carried out within the ab initio no core shell model with continuum (NCSMC) using a variety of inter-nucleon interactions and basis sizes, and yielding 6 Li deuteron separation energies between 1.3 and 1.8 MeV (Hebborn et al 2022 Phys. Rev. Lett. 129 042503). These NCSMC calculations show that the square of the ANC is strongly correlated with the separation energy over this range. In this work, we investigate the origin of this correlation using the phenomenological R-matrix, a single-channel potential and a perturbative approach. We show that this correlation occurs because the depth of the α–deuteron central potential changes by only a small relative amount as the separation energy varies. We then investigate if the ANC can be accurately extracted from α–deuteron phase shifts in an ideal case in which low-energy data are available and there are no experimental errors. We find that both R-matrix and Coulomb-modified effective-range theory (CM-ERE) yield extracted ANCs close to, although not exactly equal to, the NCSMC value, provided the extrapolation is constrained by the known position of the bound-state pole and at least three terms are included in the fit function. The R-matrix approach converges faster than the CM-ERE as the number of parameters increases and is also more robust against the inclusion of low-energy and high-energy phase shift data. Finally, our study also shows that a naive quantification of uncertainties by comparing different truncations used in both theories is not accurate, and suggests the accuracy of ANCs extracted from phase shift data needs further investigation.
Abstract not provided
The generation of reference data for deep learning models is challenging for reactive systems, and more so for combustion reactions due to the extreme conditions that create radical species and alternative spin states during the combustion process. Here, we extend intrinsic reaction coordinate (IRC) calculations with ab initio MD simulations and normal mode displacement calculations to more extensively cover the potential energy surface for 19 reaction channels for hydrogen combustion. A total of ~290,000 potential energies and ~1,270,000 nuclear force vectors are evaluated with a high quality range-separated hybrid density functional, ωB97X-V, to construct the reference data set, including transition state ensembles, for the deep learning models to study hydrogen combustion reaction.
The influence of pressure on the thermodynamic, structural, and elastic properties of the 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) molecular crystal at T = 0 is systematically studied. Calculations are carried out using density functional theory methods in a plane wave basis set with dispersion corrections for the exchange-correlation part of total energy, and ultrasoft pseudopotentials. The equilibrium unit cell parameters, the cold compression curve in the pressure range of 0–50 GPa and the sound speeds are computed. The effect of finite pressure on the molecular structure of TATB is elucidated from the analysis of relative changes in the intra- and intermolecular geometrical parameters. For the first time, the full set of elastic constants of this crystal at zero and non-zero pressures is determined from ab initio calculations. The resulted structural, elastic, and acoustic properties of TATB are shown to be in a good agreement with available experimental and theoretical data.
Ab initio calculations have an essential role in our fundamental understanding of quantum many-body systems across many subfields, from strongly correlated fermions to quantum chemistry and from atomic and molecular systems to nuclear physics. One of the primary challenges is to perform accurate calculations for systems where the interactions may be complicated and difficult for the chosen computational method to handle. Here we address the problem by introducing an approach called wavefunction matching. Wavefunction matching transforms the interaction between particles so that the wavefunctions up to some finite range match that of an easily computable interaction. This allows for calculations of systems that would otherwise be impossible owing to problems such as Monte Carlo sign cancellations. We apply the method to lattice Monte Carlo simulations of light nuclei, medium-mass nuclei, neutron matter and nuclear matter. We use high-fidelity chiral effective field theory interactions and find good agreement with empirical data. These results are accompanied by insights on the nuclear interactions that may help to resolve long-standing challenges in accurately reproducing nuclear binding energies, charge radii and nuclear-matter saturation in ab initio calculations.
Over the last decade, new developments in Similarity Renormalization Group techniques and nuclear many-body methods have dramatically increased the capabilities of ab initio nuclear structure and reaction theory. Ground and excited-state properties can be computed up to the tin region, and from the proton to the presumptive neutron drip lines, providing unprecedented opportunities to confront two- plus three-nucleon interactions from chiral Effective Field Theory with experimental data. In this contribution, I will give a broad survey of the current status of nuclear many-body approaches, and I will use selected results to discuss both achievements and open issues that need to be addressed in the coming decade.
Strong electronic correlations, interfaces, and defects, and disorder each individually challenge the theoretical methods for predictions of materials properties. These challenges are all simultaneously present in complex transition-metal-oxide interfaces and superlattices, which are known to exhibit unique and unusual properties caused by multiple coupled degrees of freedom and strong electronic correlations. Here we show that ab initio quantum Monte Carlo (QMC) solutions of the many-electron problem are now possible for the full complexity of these systems. Within a single nonempirical theoretical approach, we unambiguously establish the site-specific stability of oxygen vacancies in the (LaFeO 3 ) 2 /(SrFeO 3 ) superlattice, accounting for experimental data, and predict their migration pathways. Further, QMC calculations are now capable of playing a major role in the elucidation of many-body phenomena in complex oxides previously out of reach of first-principles theories.
In existing theoretical approaches to core-level excitations of transition-metal ions in solids relaxation and polarization effects due to the inner core hole are often ignored or described phenomenologically. Here, we set up an ab initio computational scheme that explicitly accounts for such physics in the calculation of x-ray absorption and resonant inelastic x-ray scattering spectra. Good agreement is found with experimental transition-metal L-edge data for the strongly correlated d 9 cuprate Li 2 CuO 2 , for which we also determine the absolute scattering intensities. The newly developed methodology opens the way for the investigation of even more complex d n electronic structures of group VI B to VIII B correlated oxide compounds.
We present theoretical ground state energies and their uncertainties for p -shell nuclei obtained from chiral effective field theory internucleon interactions as a function of chiral order, fitted to two- and three-body data only. We apply a Similary Renormalization Group transformation to improve the numerical convergence of the many-body calculations, and discuss both the numerical uncertainties arising from basis truncations and those from omitted induced many-body forces, as well as chiral truncation uncertainties. With complete Next-to-Next-to-Leading (N 2 LO) order two- and three-body interactions, we find significant overbinding for the ground states in the upper p -shell, but using higher-order two-body potentials, in combination with N 2 LO three-body forces, our predictions agree with experiment throughout the p -shell to within our combined estimated uncertainties. The uncertainties due to chiral order truncation are noticeably larger than the numerical uncertainties, but they are expected to become comparable to the numerical uncertainties at complete N 3 LO.
Effective interactions for elastic nucleon-nucleus scattering from first principles require the use of the same nucleon-nucleon interaction in the structure and reaction calculations, as well as a consistent treatment of the relevant operators at each order. Previous work using these interactions has shown good agreement with available data. Here, we study the physical relevance of one of these operators, which involves the spin of the struck nucleon, and examine the interpretation of this quantity in a nuclear structure context. Using the framework of the spectator expansion and the underlying framework of the no-core shell model, we calculate and examine spin-projected, one-body momentum distributions required for effective nucleon-nucleus interactions in $J=0$ nuclear states. The calculated spin-projected, one-body momentum distributions for $^4$He, $^6$He, and $^8$He display characteristic behavior based on the occupation of protons and neutrons in single particle levels, with more nucleons of one type yielding momentum distributions with larger values. Additionally, we find this quantity is strongly correlated to the magnetic moment of the $2^+$ excited state in the ground state rotational band for each nucleus considered. In conclusion, we find that spin-projected, one-body momentum distributions can probe the spin content of a $J=0$ wave function. This feature may allow future ab initio nucleon-nucleus scattering studies to inform spin properties of the underlying nucleon-nucleon interactions. The observed correlation to the magnetic moment of excited states illustrates a previously unknown connection between reaction observables such as the analyzing power and structure observables like the magnetic moment.