Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum molecular dynamics”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Equation of state for 304L stainless steel

Here, a new equation of state for 304L stainless steel has been constructed and entered in the Los Alamos National Laboratory SESAME database with material ID 4273. The equation of state is based on a multiphase approach. The solid and liquid phases were fitted to available experimental data for thermal expansion, heat capacity, and shock Hugoniot, as well as to data from density functional theory calculations of the cold curve and quantum molecular dynamics of the liquid regime. The electronic thermal contribution over the entire equation of state range is based on Tartarus (Green’s functions) average atom calculations. For this low-carbon stainless steel both the density functional and electronic thermal calculations are comprised of 19 wt% Cr, 10 wt% Ni, and the remainder Fe.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Beam energy dependence of net- Λ fluctuations measured by the STAR experiment at the BNL Relativistic Heavy Ion Collider

The measurements of particle multiplicity distributions have generated considerable interest in understanding the fluctuations of conserved quantum numbers in the quantum chromodynamics (QCD) hadronization regime, in particular near a possible critical point and near the chemical freeze-out. Net-protons and net-kaons have been used as proxies for the net-baryon number and net-strangeness, respectively. We report the measurement of efficiency- and centrality-bin width-corrected cumulant ratios (C 2 /C 1 , C 3 /C 2 ) of net- Λ distributions, in the context of both strangeness and baryon number conservation, as a function of collision energy, centrality, and rapidity. The results are for Au + Au collisions at five beam energies ( √s NN = 19.6 , 27, 39, 62.4, and 200 GeV) recorded with the Solenoidal Tracker at RHIC (STAR). We compare our results to the Poisson and negative binomial (NBD) expectations, as well as to ultrarelativistic quantum molecular dynamics (UrQMD) and hadron resonance gas (HRG) model predictions. Both NBD and Poisson baselines agree with data within the statistical and systematic uncertainties. UrQMD describes the measured net- Λ C 1 and C 3 at 200 GeV reasonably well but deviates from C 2 , and the deviation increases as a function of collision energy. The ratios of the measured cumulants show no features of critical fluctuations. The chemical freeze-out temperatures extracted from a recent HRG calculation, which was successfully used to describe the net-proton, net-kaon, and net-charge data, indicate Λ freeze-out conditions similar to those of kaons. However, large deviations are found when comparing with temperatures obtained from net-proton fluctuations. The net- Λ cumulants show a weak but finite dependence on the rapidity coverage in the acceptance of the detector, which can be attributed to quantum number conservation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Pressure-temperature equation of state of Al 2 ⁢O 3 up to 14 Mbar and 40 kK

Sapphire (Al 2 ⁢O 3 ), known for its remarkable incompressibility at ambient conditions, plays a pivotal role in both static and dynamic compression research. Accurately characterizing its equation of state (EoS) is essential for these applications. Here, we present a complete Hugoniot of Al 2 ⁢ O 3 as locus of experimentally assessed, high-precision, pressure, density and temperature states up to 14 Mbar and 43 kK. The Hugoniot is established with single shock experiments using magnetically launched hyper velocity flyers on the Z Accelerator at Sandia National Laboratories. We explore principal Hugoniot states at very high shock 𝑇 and 𝑝 in the solid phase, tracking the solid-liquid boundary and culminating at 2.4-fold compression, where data provides a direct constraint on the liquid phase. Corresponding shock release data probe thermodynamic states complementary to the Hugoniot and place additional constraints on tabular EoS models. Our findings indicate a significant deviation from existing tabular EoS models for Al 2 ⁢ O 3 dictating a comprehensive overhaul. We develop two advanced EoSs for Al 2 ⁢ O 3 the SESAME 97412 model, featuring an extensive phase diagram that includes three solid phases and the liquid phase, and the updated LEOS 2200m2 model. EoS development is assisted with Quantum Molecular Dynamics simulations. Our experimental data allows for stringent testing of our EoSs. Both models accurately capture the Hugoniot of Al 2 ⁢O 3 up to the highest pressures and temperatures. Rigorous experimental determination of extreme pressures and temperatures, paired with sophisticated models, advances the frontier of EoS development beyond 1 terapascal.

Kalita, Patricia [Sandia National Laboratories (SN↗
DOE OSTI2021

Application of machine learning in the determination of impact parameter in the 132 Sn+ 124 Sn system

Here, 132 Sn + 124 Sn collisions at a beam energy of 270 MeV/nucleon were performed at the Radioactive Isotope Beam Factory (RIBF) in RIKEN to investigate the nuclear equation of state. Reconstructing the impact parameter is one of the important tasks in the experiment as it relates to many observable. In this work, we employ three commonly used algorithms in machine learning, the artificial neural network (ANN), the convolutional neural network (CNN), and the light gradient boosting machine (LightGBM), to determine the impact parameter by analyzing either the charged particle spectra or several features simulated with events from the ultrarelativistic quantum molecular dynamics (UrQMD) model. To closely imitate experimental data and investigate the generalizability of the trained machine learning algorithms, incompressibility of nuclear equation of state and the in-medium nucleon-nucleon cross sections are varied in the UrQMD model to generate the training data. The mean absolute error Δb between the true and the predicted impact parameter is smaller than 0.45 fm if training and testing sets are sampled from the UrQMD model with the same parameter set. However, if training and testing sets are sampled with different parameter sets, Δb would increase to 0.8 fm. The generalizability of the trained machine learning algorithms suggests that these machine learning algorithms can be used reliably to reconstruct the impact parameter in experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Measurement of fifth- and sixth-order fluctuations of (net-)proton number in Au + Au collisions from phase II of the beam energy scan program at RHIC

We report high-statistics measurements of fifth- and sixth-order factorial cumulants and cumulant ratios of (net-)proton multiplicity distributions in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7–27GeV, using data from the STAR experiment collected during the Beam Energy Scan Phase II at RHIC. Protons and antiprotons are identified at midrapidity (|𝑦| < 0.5) with transverse momentum 0.4 < 𝑝 𝑇 < 2.0GeV/𝑐. The proton factorial cumulants 𝜅 4 , 𝜅 5 , and 𝜅 6 increase with order but exhibit no sign alternation within current uncertainties, offering no evidence for a two-component structure in the proton multiplicity distribution, as might be expected near a first-order phase transition. The cumulant ratios 𝐶 5 /𝐶 1 and 𝐶 6 /𝐶 2 fluctuate around zero in collisions at 0–40% centrality. The results are consistent with both the negative predictions from lattice QCD and the positive trends obtained from the ultrarelativistic quantum molecular dynamics (UrQMD) model. At $\sqrt{s_{NN}}$ ≳ 27GeV, the 𝐶 4 /𝐶 2 and 𝐶 5 /𝐶 1 results are compatible with predictions from lattice QCD, functional renormalization group (FRG), and hadron resonance gas (HRG) models, while UrQMD describes the data better at lower energies. Here, these measurements place constraints on baryon number fluctuations and offer valuable insights into the QCD phase structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Differences in Sb 2 Te 3 growth by pulsed laser and sputter deposition

High quality van der Waals chalcogenides are important for phase change data storage, thermoelectrics, and spintronics. Using a combination of statistical design of experiments and density functional theory, we clarify how out-of-equilibrium van der Waals epitaxial deposition methods can improve the crystal quality of Sb 2 Te 3 films. We compare films grown by radio frequency sputtering and pulsed laser deposition (PLD). The growth factors that influence the crystal quality for each method are different. For PLD grown films a thin amorphous Sb 2 Te 3 seed layer most significantly influences the crystal quality. In contrast, the crystalline quality of films grown by sputtering is rather sensitive to the deposition temperature and less affected by the presence of a seed layer. This difference is somewhat surprising as both methods are out-of-thermal-equilibrium plasma-based methods. Non-adiabatic quantum molecular dynamics simulations show that this difference originates from the density of excited atoms in the plasma. The PLD plasma is more intense and with higher energy than that used in sputtering, and this increases the electronic temperature of the deposited atoms, which concomitantly increases the adatom diffusion lengths in PLD. In contrast, the adatom diffusivity is dominated by the thermal temperature for sputter grown films. These results explain the wide range of Sb 2 Te 3 and superlattice crystal qualities observed in the literature. Finally, these results indicate that, contrary to popular belief, plasma-based deposition methods are suitable for growing high quality crystalline chalcogenides.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Evidence of hexagonal germanium grains on annealed monolayer MoS 2

Growing three-dimensional (3D) materials on two-dimensional (2D) van der Waals surface has shown its effectiveness in overcoming materials incompatibility for stacking transferrable membranes toward advanced device manufacturing. Herein, we demonstrate that the nucleation of hexagonal germanium (Ge) grains within a continuous crystalline film, which has been unfeasible through traditional epitaxy techniques, is realized by chemical vapor deposition on top of n-type monolayer molybdenum disulfide (MoS 2 ) substrates. Suggested by quantum molecular dynamics calculation, the hexagonal Ge nucleation is thermodynamically preferable to cubic Ge when growing on monolayer MoS 2 with sulfur vacancies. The strained hexagonal Ge grains have been confirmed by transmission electron microscopy analyses from both real space and reciprocal space. Scanning probe microscopy shows that the hexagonal Ge film possesses higher reflectivity in infrared spectral range, implying a higher carrier concentration resulted from the narrower band gap, as compared to cubic Ge.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Transport model comparison studies of intermediate-energy heavy-ion collisions

Transport models are the main method to obtain physics information on the nuclear equation of state and in-medium properties of particles from low to relativistic-energy heavy-ion collisions. The Transport Model Evaluation Project (TMEP) has been pursued to test the robustness of transport model predictions in reaching consistent conclusions from the same type of physical model. To this end, calculations under controlled conditions of physical input and set-up were performed with various participating codes. These included both calculations of nuclear matter in a box with periodic boundary conditions, which test separately selected ingredients of a transport code, and more realistic calculations of heavy-ion collisions. Over the years, six studies have been performed within this project. In this intermediate review, we summarize and discuss the present status of the project. We also provide condensed descriptions of the 26 participating codes, which contributed to some part of the project. These include the major codes in use today. After a compact description of the underlying transport approaches, we review the main results of the studies completed so far. They show, that in box calculations the differences between the codes can be well understood and a convergence of the results can be reached. These studies also highlight the systematic differences between the two families of transport codes, known under the names of Boltzmann–Uehling–Uhlenbeck (BUU) and Quantum Molecular Dynamics (QMD) type codes. However, when the codes were compared in full heavy-ion collisions using different physical models, as recently for pion production, they still yielded substantially different results. This calls for further comparisons of heavy-ion collisions with controlled models and of box comparisons of important ingredients, like momentum-dependent fields, which are currently underway. Furthermore, our evaluation studies often indicate improved strategies in performing transport simulations and thus can provide guidance to code developers. Results of transport simulations of heavy-ion collisions from a given code will have more significance if the code can be validated against benchmark calculations such as the ones summarized in this review.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Tellurene Photodetector with High Gain and Wide Bandwidth

Two-dimensional (2D) semiconductors have been extensively explored as a new class of materials with great potential. In particular, black phosphorus (BP) has been considered to be a strong candidate for applications such as high-performance infrared photodetectors. However, the scalability of BP thin film is still a challenge, and its poor stability in the air has hampered the progress of the commercialization of BP devices. Herein, we report the use of hydrothermal-synthesized and air-stable 2D tellurene nanoflakes for broadband and ultrasensitive photodetection. The tellurene nanoflakes show high hole mobilities up to 458 cm 2 /V·s at ambient conditions, and the tellurene photodetector presents peak extrinsic responsivity of 383 A/W, 19.2 mA/W, and 18.9 mA/W at 520 nm, 1.55 μm, and 3.39 μm light wavelength, respectively. Because of the photogating effect, high gains up to 1.9 × 10 3 and 3.15 × 10 4 are obtained at 520 nm and 3.39 μm wavelength, respectively. At the communication wavelength of 1.55 μm, the tellurene photodetector exhibits an exceptionally high anisotropic behavior, and a large bandwidth of 37 MHz is obtained. The photodetection performance at different wavelength is further supported by the corresponding quantum molecular dynamics (QMD) simulations. Our approach has demonstrated the air-stable tellurene photodetectors that fully cover the short-wave infrared band with ultrafast photoresponse.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

The structure of liquid carbon elucidated by in situ X-ray diffraction

Carbon has a central role in biology and organic chemistry, and its solid allotropes provide the basis of much of our modern technology. However, the liquid form of carbon remains nearly uncharted, and the structure of liquid carbon and most of its physical properties are essentially unknown. But liquid carbon is relevant for modelling planetary interiors and the atmospheres of white dwarfs, as an intermediate state for the synthesis of advanced carbon materials, inertial confinement fusion implosions, hypervelocity impact events on carbon materials and our general understanding of structured fluids at extreme conditions. Here we present a precise structure measurement of liquid carbon at pressures of around 1 million atmospheres obtained by in situ X-ray diffraction at an X-ray free-electron laser. Our results show a complex fluid with transient bonding and approximately four nearest neighbours on average, in agreement with quantum molecular dynamics simulations. The obtained data substantiate the understanding of the liquid state of one of the most abundant elements in the universe and can test models of the melting line. The demonstrated experimental abilities open the path to performing similar studies of the structure of liquids composed of light elements at extreme conditions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Real-space density kernel method for Kohn–Sham density functional theory calculations at high temperature

Kohn–Sham density functional theory calculations using conventional diagonalization based methods become increasingly expensive as temperature increases due to the need to compute increasing numbers of partially occupied states. In this work, we present a density matrix based method for Kohn–Sham calculations at high temperatures that eliminates the need for diagonalization entirely, thus reducing the cost of such calculations significantly. Specifically, we develop real-space expressions for the electron density, electronic free energy, Hellmann–Feynman forces, and Hellmann–Feynman stress tensor in terms of an orthonormal auxiliary orbital basis and its density kernel transform, the density kernel being the matrix representation of the density operator in the auxiliary basis. Using Chebyshev filtering to generate the auxiliary basis, we next develop an approach akin to Clenshaw–Curtis spectral quadrature to calculate the individual columns of the density kernel based on the Fermi operator expansion in Chebyshev polynomials and employ a similar approach to evaluate band structure and entropic energy components. We implement the proposed formulation in the SPARC electronic structure code, using which we show systematic convergence of the aforementioned quantities to exact diagonalization results, and obtain significant speedups relative to conventional diagonalization based methods. Finally, we employ the new method to compute the self-diffusion coefficient and viscosity of aluminum at 116 045 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Palladium at high pressure and high temperature: A combined experimental and theoretical study

Palladium is one of the most important technological materials, yet its phase diagram remains poorly understood. At ambient conditions, its solid phase is face-centered cubic (fcc). However, another solid phase of Pd, body-centered cubic (bcc), was very recently predicted in two independent theoretical studies to occur at high pressures and temperatures. In this work, we report an experimental study on the room-temperature equation of state (EOS) of Pd to a pressure of 80 GPa, as well as a theoretical study on the phase diagram of Pd including both fcc-Pd and bcc-Pd. Our theoretical approach consists in ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries. We obtain the melting curves of both fcc-Pd and bcc-Pd and an equation for the fcc–bcc solid–solid phase transition boundary as well as the thermal EOS of Pd which is in agreement with experimental data and QMD simulations. We uncover the presence of another solid phase of Pd on its phase diagram, namely, random hexagonal close-packed (rhcp), and estimate the location of the rhcp-bcc solid–solid phase transition boundary and the rhcp–fcc–bcc triple point. We also discuss the topological similarity of the phase diagrams of palladium and silver, the neighbor of Pd in the periodic table. We argue that Pd is a reliable standard for shock-compression studies and present the analytic model of its principal Hugoniot in a wide pressure range.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Alkali hydroxide (LiOH, NaOH, KOH) in water: Structural and vibrational properties, including neutron scattering results

Structural and vibrational properties of aqueous solutions of alkali hydroxides (LiOH, NaOH, and KOH) are computed using quantum molecular dynamics simulations for solute concentrations ranging between 1 and 10M. Element-resolved partial radial distribution functions, neutron and x-ray structure factors, and angular distribution functions are computed for the three hydroxide solutions as a function of concentration. The vibrational spectra and frequency-dependent conductivity are computed from the Fourier transforms of velocity autocorrelation and current autocorrelation functions. Our results for the structure are validated with the available neutron data for 17M concentration of NaOH in water [Semrouni et al., Phys. Chem. Chem. Phys. 21, 6828 (2019)]. We found that the larger ionic radius [rLi+

Chemistry↗
DOE OSTI2024

Atomic cluster expansion potential for large scale simulations of hydrocarbons under shock compression

We present an Atomic Cluster Expansion (ACE) machine learned potential developed for high-fidelity atomistic simulations of hydrocarbons, targeting pressures and temperatures near and above supercritical fluid regimes for molecular fluids. A diverse set of stoichiometries were covered in training, including 1:0 (pure carbon), 1:4 (methane), and 1:1 (benzene), and rich bonding environments sampled at supercritical temperatures, hydrogen rich, reactive mixtures where metastable stoichiometries arise, including 1:2 (ethylene) and 1:3 (ethane). A high-fidelity training database was constructed by performing large-scale quantum molecular dynamic simulations [density functional theory (DFT) MD] of diamond, graphite, methane, and benzene. A novel approach to selecting structures from DFT MD is also presented, which allows for the rapid selection of unique DFT MD frames from complex trajectories. Comparisons to DFT and experimental data demonstrate that the presented ACE potential accurately reproduces isotherms, carbon melting curves, radial distribution functions, and shock Hugoniots for carbon and hydrocarbon systems for pressures up to 100 GPa and temperatures up to 6000 K for hydrocarbon systems and up to 9000 K for pure carbon systems. This work delivers a potential that can be used for accurate, large-scale simulations of shocked hydrocarbons and demonstrates a methodology for fitting and validating machine learning interatomic potentials to complex molecular environments, which can be applied to energetic materials in future works.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Synthesis and investigation into explosive sensitivity for a series of new picramide explosives

Tailoring the molecular properties that govern energetic material sensitivity is essential to improve safety and help develop new energetic materials. Despite this need, understanding the complex chemistry and physics of explosive initiation and propagation is still a challenge. Recent work by our group has reinforced the view that explosive sensitivity under sub-shock conditions is connected to the strength of the weakest covalent bond in the molecule, that is, its “trigger linkage.” These correlations have been observed with different classes of energetic molecules and indicate that “trigger linkage” bond breaking, and heat of explosion are good indicators for the sensitivity trends. Herein we report the synthesis of aliphatic energetic materials with ethane, propane and neopentane backbones. Experimental and computational studies show that the trigger linkage model, based on results from quantum molecular dynamics simulations, correctly predicts trends observed in the impact sensitivity of the molecules. However, while the model predicts the impact sensitivities of the ethane series, the neopentane series has higher impact sensitivities than predicted, which is presumably influenced by crystal packing effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Enhanced dilepton emission from a phase transition in dense matter

It is demonstrated that the presence of a phase transition in heavy-ion collisions, at beam energies that probe dense quantum chromodynamics (QCD) matter, leads to a significant enhancement of the dilepton yield at low invariant mass and requires the embedding of low beam energies per produced pion due to the extended emission time. In addition, the temperature of low-mass dileptons shows a modest decrease due to the mixed phase. The emission of dileptons in the SIS18–SIS100 beam energy range is studied by augmenting the ultra-relativistic quantum molecular dynamics (UrQMD) transport model with a realistic density-dependent equation of state, as well as two different phase transitions. This is achieved by extending the molecular dynamics interaction part of the UrQMD model to a density-dependent interaction potential with a high-density minimum, leading to a phase transition and metastable coexisting high-density states. Together with a high-precision measurement, these simulations will be able to constrain the existence of a phase transition in QCD up to densities of several times nuclear saturation density.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Electric-field-induced crossover of polarization reversal mechanisms in Al 1–x Sc x N ferroelectrics

Scandium-doped aluminum nitride, Al 1–x Sc x N, represents a new class of displacive ferroelectric materials with high polarization and sharp hysteresis along with high-temperature resilience, facile synthesizability and compatibility with standard CMOS fabrication techniques. The fundamental physics behind the transformation of unswitchable piezoelectric AlN into switchable Al–Sc–N ferroelectrics depends upon important atomic properties such as local structure, dopant distributions and the presence of competing mechanism of polarization switching in the presence of an applied electric-field that have not been understood. We computationally synthesize Al 1–x Sc x N to quantify the inhomogeneity of Sc distribution and phase segregation, and characterize its crystal and electronic structure as a function of Sc-doping. Nudged elastic band calculations of the potential energy surface and quantum molecular dynamics simulations of direct electric-field-driven ferroelectric switching reveal a crossover between two polarization reversal mechanisms—inhomogeneous nucleation-and-growth mechanism originating near Sc-rich regions in the limit of low applied fields and nucleation-limited-switching in the high-field regime. Finally, understanding polarization reversal pathways for these two mechanisms as well as the role of local Sc concentration on activation barriers provides design rules to identify other combinations of dopant elements, such as Zr, Mg etc. to synthesize superior AlN-based ferroelectric materials.

77 NANOSCIENCE AND NANOTECHNOLOGY↗