Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic physics”

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

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

At least 1,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Pressure dependent magnetic properties on bulk CrBr 3 single crystals

The van der Waals class of materials offer an approach to two-dimensional magnetism as their spin fluctuations can be tuned upon exfoliation of layers. Moreover, it has recently been shown that spin-lattice coupling and long-range magnetic ordering can be modified with pressure in van der Waals materials. In this work, the magnetic properties of quasi two-dimensional CrBr 3 are reported applying hydrostatic pressure. The application of pressure up to 0.844 GPa shows a 1.77% decrease in saturation magnetization with a decrease in the Curie temperature from 33.05 to 30.41 K. Density functional theory calculations with pressure up to 1 GPa show a reduction in volume and interplanar distance as pressure increases. To further understand the magnetic properties with applied pressure, the magnetocrystalline anisotropy energy (MAE) and exchange coupling parameter (J) are calculated. There is small decrease in MAE and the first nearest neighbor interaction (J 1 ) (U = 2.7 eV and J = 0.7 eV) is increasing with respect to increasing pressure. Overall, CrBr 3 displays ferromagnetic interlayer coupling and the calculated exchange coupling and MAE parameters match well with the observations from the experimental work.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Implicit fast sweeping method for hyperbolic systems of conservation laws

Implicit time-accurate methods are often used to integrate stiff problems where explicit schemes impose severe time step restrictions. This paper presents an efficient numerical framework based on the Fast Sweeping Method (FSM) for solving linear and nonlinear hyperbolic systems of conservation laws. The solution at each discrete location is computed by sweeping the numerical domain in several predetermined directions that follow the causality of the characteristic families. The use of a fractional step strategy eliminates the need for a solution selection criterion while one-sided stencils limit the number of sweeps to at most 2 d for d space dimensions. This work focuses on the first-order implicit upwind method since it constitutes the building block for high-order conservative schemes. For problems where the degree of stiffness evolves over time, implicit-explicit hybridization can be accomplished with the same algorithm by simply switching the stencil at each time level. As opposed to traditional implicit solvers, the sweeping method does not require a local time linearization of the fluxes thereby preserving the nonlinear stability properties of the original implicit scheme. It also avoids the large computational and memory requirements associated with solving large block-diagonal systems of equations. Here, a series of one- and two-dimensional test cases are presented for the inviscid Burgers' equation and the reactive Euler equations. The results indicate that the implicit FSM can allow a major reduction in the number of time steps even in the presence of discontinuous solution profiles.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Towards sharp error analysis of extended Lagrangian molecular dynamics

The extended Lagrangian molecular dynamics (XLMD) method provides a useful framework for reducing the computational cost of a class of molecular dynamics simulations with constrained latent variables. The XLMD method relaxes the constraints by introducing a fictitious mass ε for the latent variables and solving a set of singularly perturbed ordinary differential equations. While favorable numerical performance of XLMD has been demonstrated in several different contexts in the past decade, mathematical analysis of the method remains scarce. Here, we propose the first error analysis of the XLMD method in the context of a classical polarizable force field model. While the dynamics with respect to the atomic degrees of freedom are general and nonlinear, the key mathematical simplification of the polarizable force field model is that the constraints on the latent variables are given by a linear system of equations. We prove that when the initial value of the latent variables is compatible in a sense that we define, XLMD converges as the fictitious mass ε is made small with $\mathscr{O}$(ε) error for the atomic degrees of freedom and with $\mathscr{O}$($\sqrt{ε}$) error for the latent variables, when the dimension of the latent variable d' is 1. Furthermore, when the initial value of the latent variables is improved to be optimally compatible in a certain sense, we prove that the convergence rate can be improved to $\mathscr{O}$(ε) for the latent variables as well. Numerical results verify that both estimates are sharp not only for d'=1, but also for arbitrary d'. In the setting of general d', we do obtain convergence, but with the non-sharp rate of $\mathscr{O}$($\sqrt{ε}$) for both the atomic and latent variables.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Multinuclear PFGSTE NMR description of 39 K, 23 Na, 7 Li, and 1 H specific activation energies governing diffusion in alkali nitrite solutions

While pulsed field gradient stimulated echo nuclear magnetic resonance (PFGSTE NMR) spectroscopy has found widespread use in the quantification of self-diffusivity for many NMR-active nuclei, extending this technique to uncommon nuclei with unfavorable NMR properties remains an active area of research. Potassium-39 ( 39 K) is an archetypical NMR nucleus exhibiting an unfavorable gyromagnetic ratio combined with a very low Larmor frequency. Despite these unfavorable properties, this work demonstrates that 39 K PFGSTE NMR experiments are possible in aqueous solutions of concentrated potassium nitrite. Further, analysis of the results indicates that 39 K NMR diffusometry is feasible when the nuclei exhibit spin–lattice and spin–spin relaxation coefficients on the order of 60–100 ms and 50–100 ms, respectively. The diffusivity of 39 K followed Arrhenius behavior, and comparative 23 Na, 7 Li, and 1 H PFGSTE NMR studies of equimolal sodium nitrite and lithium nitrite solutions led to correlations between the enthalpy of hydration with the activation energy governing self-diffusion of the cations and also of water. Realizing the feasibility of 39 K PFGSTE NMR spectroscopy has a widespread impact across energy sciences because potassium is a common alkali element in energy storage materials and other applications.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Multidimensional molecular high-harmonic spectroscopy: A road map for charge migration studies

This work discusses recent experimental and theoretical analyses of high-harmonic spectroscopy in small molecules, with the aim of characterizing charge migration. We discuss the formulation of molecular high-harmonic spectra followed by methods to extract molecular-target-specific information, both in experiments and ab initio simulations. We present measured and simulated high-harmonic spectra from carbon dioxide and carbonyl sulfide to illustrate the necessity for multidimensional analyses for high-harmonic spectroscopies that include both the spectral amplitude and phase. Leveraging these results, we examine how such multidimensional analyses pave the way for the study of charge migration with high-harmonic spectroscopy and illustrate the beneficial role a static molecular feature can play when probing dynamics. Finally, we briefly expand our scope with an outlook on the critical role of integrating theoretical and experimental approaches, beyond just high-harmonic spectroscopy, for the development of versatile harmonic spectroscopic probes of charge migration.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

The 1.66 μm spectrum of the ethynyl radical, CCH

Frequency-modulated diode laser transient absorption spectra of the ethynyl radical have been recorded at wavelengths close to 1.66 μm. The observed spectrum includes strong, regular, line patterns. The two main bands observed originate in the ground $\tilde{χ}$ 2 Σ + state and its first excited bending vibrational level of 2 $Π$ symmetry. The upper states, of 2 Σ + symmetry at 6055.6 cm -1 and 2 $Π$ symmetry at 6413.5 cm -1 , respectively, had not previously been observed and the data were analyzed in terms of an effective Hamiltonian representing their rotational and fine structure levels to derive parameters which can be used to calculate rotational levels up to J = 37/2 for the 2 $Π$ level and J = 29/2 for the 2 Σ one. Additionally, a weaker series of lines have been assigned to absorption from the second excited bending, (0 2 0), level of 2 Σ symmetry, to a previously observed state of 2 $Π$ symmetry near 6819 cm -1 . These strong absorption bands at convenient near-IR laser wavelengths will be useful for monitoring CCH radicals in chemical systems.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Deep glassy state dynamic data challenge glass models: Configurational entropy models

This work was the result of the student's (Dongjie Chen) qualifying exam response. In it the experimental observation from our labs that the relaxation times of ultrastable glasses (a 20 million year old amber and a vapor deposited amorphous teflon) do not follow the expected Vogel-Fulcher or WLF types of divergence as the temperatures fall below the glass transition temperature. In the work we evaluated several theories by considering them in different ways. We used them as generally used and show that they do not capture the non-diverging behavior. Because the models were all entropy-based models, we also considered the possibility that the temperature dependence of the entropy used in the models is incorrect and chose two ways of estimating the entropy that do not follow the classical theoretical models. First we used the actually measured calorimetric response for the two ultrastable glasses, and we used the Milchev model, which does not go to zero above 0 Kelvin. Both of these approaches led to improved agreement between the theories and the experimental data. Finally, we also asked what the temperature dependence of the entropy would have to be so that the models agreed with the observed data. In some cases the entropy so-calculated was reasonable, while in others it was unphysical. In sum, the important point of the work is that the observation of non-diverging time-scales that suggest that there is no ideal glass transition can be reasonably captured by current theories if they take more realistic models of the entropy than is usually the case.

36 MATERIALS SCIENCE↗
DOE OSTI2022

An ab initio molecular dynamics investigation of the thermophysical properties of molten NaCl-MgCl 2

Molten salts have many applications in the nuclear and solar energy industries for thermal storage and heat transfer applications. However, there is a knowledge gap in molten salt thermophysical properties which hinders the technical readiness level of molten salt applications, especially in the nuclear industry. A common method of investigating new materials is through ab initio Molecular Dynamics (AIMD) simulations which is an effective tool to investigate structural and thermophysical properties at realistic temperatures. NaCl-MgCl 2 is an inexpensive salt that is a good candidate for use as a heat transfer medium in solar power applications or in the secondary loop of a nuclear reactor. In this article, the thermophysical properties of NaCl-MgCl 2 are computed via AIMD calculations to supplement the limited experimental studies in the literature. Here a wide range of compositions and temperatures for the pseudo-binary NaCl-MgCl 2 were used to calculate the density, heat capacity, compressibility, enthalpy of mixing, and volumetric thermal expansion coefficient. AIMD is shown to accurately model the densities of molten NaCl-MgCl 2 as there is good agreement with the available literature. This work observed a transition to a monotonic increase of the density with respect to MgCl 2 composition occurring above 1100 K. The heat capacity values increase uniformly with respect to concentration of MgCl 2 at a rate of 2.85 J/mol-K per 10 mol% of MgCl 2 . Select thermophysical properties are fit to a Redlich-Kister expansion for utilization in multiphysics simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Computational determination of a primary diffusion mode in γ U-10Mo under irradiation

Low enriched uranium (< 20 % 235 U)-molybdenum (U-Mo) monolithic fuel is the primary candidate for high performance research and test reactors, and is in the process of being qualified to replace highly enriched uranium (≥ 20 % 23 5U) fuel. As part of the qualification process, it is critical to understand and predict the behavior of fission gas bubbles under irradiation, which affects fuel swelling and fuel failure. Mechanistic fuel models are being developed that can both reproduce the existing experimental data for fuel swelling, and be further applied to irradiation conditions beyond the experimental scope. Diffusion of species under irradiation conditions is an important parameter in the mechanistic fuel models; however, no temperature-relevant experimental diffusion data exists. In the present work, radiation-enhanced diffusion coefficients of U, Mo, and Xe in λU-10wt.%Mo were calculated in the temperature range between 300 K and 1400 K via rate-theory models and molecular dynamics simulations with an embedded-atom method interatomic potential for the U-Mo-Xe system. Accordingly, total diffusion coefficients under relevant irradiation conditions are determined using previously obtained intrinsic thermal diffusion and radiation-driven diffusion coefficients, as well as the newly calculated radiation-enhanced diffusion coefficients presented herein. Radiation-enhanced diffusion of U and Mo was dominant in the intermediate temperature range, whereas radiation-enhanced diffusion of Xe did not significantly contribute to total diffusion of Xe at the relevant fission rate densities. Radiation-enhanced diffusion of Xe became faster than both intrinsic thermal diffusion and radiation-driven diffusion at a fission rate density of 5 x 10 22 fissions/m 3 /s, which is higher than the typical fission rate density range in research reactors. The temperature regime that radiation-enhanced diffusion of each element dominated was dependent on the fission rate density. Finally, the total diffusion coefficients of U, Mo, and Xe, updated in this work, will be utilized as parameters in the mechanistic fuel models to help predict the behavior of fission gas bubbles under irradiation more accurately.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

NMR investigation of proton transport in polybenzimidazole/polyphosphoric acid membranes prepared via novel synthesis route

Here, in this work, we present a molecular-level view of the structural changes in polybenzimidazole (PBI) membranes doped with phosphoric acid (PA) generated by a novel membrane fabrication technique. The modified PA doped membranes displayed unprecedented ionic conductivity at elevated temperatures, comparable to the starting PBI gel membranes prepared by the commonly used PPA process, while also exhibiting enhanced mechanical properties. To elucidate the cause of these effects, we used multi-nuclear ( 1 H, 13 C, 31 P) 1D PFG, MAS and CPMAS, and 2D HETCOR magnetic resonance (NMR) techniques to characterize the structure and dynamics of the modified film and the original gel PBI membranes. 1 H diffusivity measurements show significantly enhanced proton diffusivity, both in magnitude and in lower activation energy, which were consistent with its high ionic conductivity despite the lower PA content compared to the original gel membrane. CPMAS experiments further substantiated that the location of these distinct phosphate environments as being close to the polymer network. Finally, the phosphate groups in the modified membrane were revealed to be strongly bounded to each other and to the PBI polymer backbone in contrast to the combination of weakly and strongly bounded groups of the original gel membrane.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Vacuum-ultraviolet absorption cross-sections of functionalized cyclic hydrocarbons: Five-membered rings

Cross-sections were measured at 50 °C and 800 Torr using differential absorption spectroscopy in the vacuum-ultraviolet from 5.17 – 9.92 eV for a series of five-membered cyclic hydrocarbons: cyclopentane, cyclopentene, methylcyclopentane, tetrahydrofuran, furan, 2-methyltetrahydrofuran, 2,3-dihydro-5-methylfuran, tetrahydrofuran-3-one, 2-methyl-1,3-dioxolane, cyclopentanol, 2-hydroxytetrahydrofuran, cyclopentanone, 2-cyclopenten-1-one, 3-cyclopenten-1-one, 2-methylcyclopentan-1-one, 2-methylcyclopent-2-en-1-one, 1,2-epoxycyclopentane, 2,3-epoxycyclopentanone, and 3,4-epoxycyclopentanone. The majority of the species are functionalized and are directly related to chain-reactions relevant to combustion, such as substituted cyclic ethers and cycloalkenes. Uncertainty calculations were conducted as a function of photon energy to quantify errors associated with experimental repeatability, signal-to-noise ratio, and gas-phase concentration, yielding an upper limit of 5% in spectral regions where absorption is >2%. Here, the primary objective of the present work is to provide quantitative absorption cross-sections of species commonly produced in combustion reactions. Comparing qualitative differences among the five-membered cyclic species, in order to analyze the effect of functional groups on absorption, is a secondary objective. The majority of the spectra are reported for the first time and enable quantitative isomer-resolved speciation measurements in combustion chemistry.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Line coupling and line mixing effects on calculated widths of symmetric-top molecules with the k-degeneracy: A theoretical study of N 2 -, O 2 -, and air-broadened lines of CH 3 I

Calculations of the N 2 -, O 2 -, and air-broadened widths, together with their temperature dependence exponents have been made for transitions of CH 3 I in the ν 5 and ν 6 bands. The calculations are based on a semi-classical line shape formalism developed by the current authors through modifying and refining the Robert-Bonamy formalism. In recent years, we have applied this formalism for linear molecules, symmetric-top molecules with inversion symmetry, and asymmetric-top molecules. For symmetric-top molecules with the k degeneracy such as CH 3 I, the formalism has a new feature. In this case, one should consider each of the CH 3 I transitions labeled by k i or f ≠ 0 as a doublet. Then, one needs to consider the effects of the line mixing process between these two components. Furthermore, comparisons of our theoretical predictions with some data available demonstrate a very reasonable agreement. Finally, we propose new experiments at higher perturber pressures that would enable one to check the theoretically calculated relaxation matrices and to extend the analysis to the inter-doublet mixing effects.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Low-lying states and total internal partition sums of CH

The electronic structure and spin-orbit states of the CH radical have been systematically investigated using multi-reference configuration interaction (MRCI) and single-reference coupled-cluster (CC) methods. These calculations were performed in conjunction with large correlation-consistent basis sets of quadruple-, quintuple-, and sextuple-ζ quality. To achieve high accuracy, electronic energies for all states were extrapolated to the complete basis set (CBS) limit, enabling the detailed construction of potential energy curves and determination of reliable spectroscopic constants. Spin-orbit coupling effects were explicitly incorporated, and vibrational energy levels were computed via Numerov analysis. Furthermore, the resulting values exhibit good to excellent agreement with available experimental data. Dipole moment and transition dipole moment curves were evaluated to assess the opacity characteristics of CH, revealing that transitions such as Χ 2 Π (u′′ = 0) → Α 2 Δ (u′ = 0), Χ 2 Π (u′′ = 0) → Β 2 Σ − (u′ = 0), Χ 2 Π (u′′ = 0) → C 2 Σ + (u′ = 0), and Χ 2 Π (u′′ = 0) → D 2 Σ + (u′ = 3) are particularly probable. Finally, the total internal partition function sum (TIPS) of CH was computed over a broad temperature range (10–30,000 K) based on our high-accuracy ab initio results.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Effect of stoichiometry and hydration level on water domain size and transport in poly(aryl piperidinium) alkaline anion-exchange membranes

Alkaline water electrolysis holds promise in decarbonizing the global economy by enabling renewable hydrogen production with non-precious group metal catalysts. Anion exchange membranes are an important component of alkaline water electrolyzers and would ideally be durable while allowing for high hydroxide conductivity. The poly(aryl piperidinium) (PAP) class of polymers has attracted recent interest due to their good mechanical robustness and high ionic conductivity. Here, in this work, we perform atomistic molecular dynamics (MD) simulations of several PAP polymers at experimentally relevant hydration levels and polymer ion exchange capacities (IECs) to gain nanoscale insight into their properties and to help elucidate the trade-offs that result from tuning the IECs through the polymer stoichiometry. Our MD-predicted macroscopic polymer properties were found to be in good agreement with experimentally available polymer swelling ratios, water-occupied volumes, X-ray scattering, and ionic conductivities. The models show that for hydration levels greater than 8H 2 O per cation a single water cluster will form that percolates through the system. The growth in water cluster size results in large polymer swelling, the creation of larger channels with widths of 7 Å or larger, and nanophase separation between the hydrophilic domains and the polymer with characteristic length scales of approximately 20–30 Å. The experimentally observed lack of a strong X-ray scattering peak at low wavevectors can be explained by a cancellation between the polymer-polymer/water-water and polymer-water correlations and not a loss in nanophase separation. The overlap in coordination environments of the hydroxide oxygen and polymer nitrogen atoms implies that vehicular diffusion between cationic groups could play a role in hydroxide transport. The polymers' hydroxide and water diffusion constants increase by approximately an order of magnitude between hydration levels of 8 and 20H 2 O per cation. However, there are diminishing returns in hydroxide diffusion constant once the IEC exceeds 2.4 meq/g.

08 HYDROGEN↗
DOE OSTI2021

Flow reduction of hydrocarbon liquid in silica nanochannel: Insight from many-body dissipative particle dynamics simulations

A modified many-body dissipative particle dynamics (mDPD) model recently developed for realistic mesoscale multiphase flow simulations is rigorously parameterized, calibrated, and applied for elucidating the flow mechanisms of hydrocarbon liquids (i.e., heptane in this work) in amorphous silica cylindrical nanochannels with inner diameters ranging from 4.5 to 22.5 nm. The simulation results suggest the presence of a strong threshold of pressure gradient under which heptane cannot be driven to flow. The threshold for the 4.5 nm diameter pore is 10 to 100 times as high as for the 9–22.5 nm diameter pore, highlighting a remarkable nanoconfinement effect. Fluid viscosity is found to exhibit a shear-thinning phenomenon with intensity to weaken with increasing channel diameter — a phenomenon not observed in nanochannel flow of liquid water and gas in literature. Most remarkably, the radial profiles of average longitudinal flow velocity fitted by the modified Hagen-Poiseuille equation showed a negative slip length (–2.5% to –0.5% relative to the diameter) and a reduction of apparent permeability by 16% to 23%. Furthermore, this finding suggests silica nanochannels tend to deter hydrocarbon flow, a phenomenon that is opposed to the flow enhancement reported in most of the prior nanochannel flow studies in literature.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Optimizing lossy state preparation for quantum sensing using Hamiltonian engineering

One of the most prominent platforms for demonstrating quantum sensing below the standard quantum limit is the spinor Bose–Einstein condensate. While a quantum advantage using several tens of thousands of atoms has been demonstrated in this platform, it faces an important challenge: atom loss. Atom loss is a Markovian error process modeled by Lindblad jump operators, and a no-go theorem, which we also show here, states that the loss of atoms in all spin components reduces the quantum advantage to a constant factor. Here, we show that this no-go theorem can be circumvented if we constrain atom losses to a single spin component. Moreover, we show that in this case, the maximum quantum Fisher information with N atoms scales as N 3/2 , establishing that a scalable quantum advantage can be achieved despite atom loss. Although Lindblad jump operators are generally non-Hermitian and non-invertible, we use their Moore–Penrose inverse to develop a framework for constructing several states with this scaling of Fisher information in the presence of losses. We use Hamiltonian engineering with realistic Hamiltonians to develop experimental protocols for preparing these states. Finally, we discuss possible experimental techniques to constrain the losses to a single spin mode.

74 ATOMIC AND MOLECULAR PHYSICS↗