Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “potential energy surface”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Molecular bond stabilization in the strong-field dissociation of O 2 +

We theoretically examine the rotational and vibrational dynamics of O 2 + molecular ions exposed to intense, short laser pulses for conditions realized in contemporary pump-probe experiments. We solve the time-dependent Schrödinger equation within the Born-Oppenheimer approximation for an initial distribution of randomly aligned molecular ions. For fixed peak intensities, our numerical results show that total, angle-integrated O 2 + → O( 3 P) + O + ( 4 S 0 ) dissociation yields do not monotonically increase with increasing infrared-probe pulse duration. We find this pulse-duration-dependent stabilization to be consistent with the transient trapping of nuclear probability density in a light-induced (bond-hardening) potential-energy surface and robust against rotational excitation. We analyze this stabilization effect and its underlying bond-hardening mechanism (i) in the time domain, by following the evolution of partial nuclear probability densities associated with the dipole-coupled O 2 + (a 4 Π u ) and O 2 + (f 4 Π g ) cationic states, and (ii) in the frequency domain, by examining rovibrational quantum-beat spectra for the evolution of the partial nuclear probability densities associated with these states. In this work, our analysis reveals the characteristic timescale for the bond-hardening mechanism in O 2 + and explains the onset of bond stabilization for sufficiently long pulse durations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Deep neural network for the dielectric response of insulators

In this work, we introduce a deep neural network to model in a symmetry preserving way the environmental dependence of the centers of the electronic charge. The model learns from ab initio density functional theory, wherein the electronic centers are uniquely assigned by the maximally localized Wannier functions. When combined with the deep potential model of the atomic potential energy surface, the scheme predicts the dielectric response of insulators for trajectories inaccessible to direct ab initio simulation. The scheme is nonperturbative and can capture the response of a mutating chemical environment. We demonstrate the approach by calculating the infrared spectra of liquid water at standard conditions, and of ice under extreme pressure, when it transforms from a molecular to an ionic crystal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Quantum Monte Carlo determination of the principal Hugoniot of deuterium

Here, we present coupled electron-ion Monte Carlo results for the principal Hugoniot of deuterium together with an accurate study of the initial reference state of shock-wave experiments. We discuss the influence of nuclear quantum effects, thermal electronic excitations, and the convergence of the potential energy surface by wave-function optimization within variational Monte Carlo and projection quantum Monte Carlo methods. Compared to a previous study, our calculations also include low pressure-temperature (P,T) conditions resulting in close agreement with experimental data, while our revised results at higher (P,T) conditions still predict a more compressible Hugoniot than experimentally observed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Heat transport in liquid water from first-principles and deep neural network simulations

In this work, we compute the thermal conductivity of water within linear response theory from equilibrium molecular dynamics simulations, by adopting two different approaches. In one, the potential energy surface (PES) is derived on the fly from the electronic ground state of density functional theory (DFT) and the corresponding analytical expression is used for the energy flux. In the other, the PES is represented by a deep neural network (DNN) trained on DFT data, whereby the PES has an explicit local decomposition and the energy flux takes a particularly simple expression. By virtue of a gauge invariance principle, established by Marcolongo, Umari, and Baroni, the two approaches should be equivalent if the PES were reproduced accurately by the DNN model. We test this hypothesis by calculating the thermal conductivity, at the GGA (PBE) level of theory, using the direct formulation and its DNN proxy, finding that both approaches yield the same conductivity, in excess of the experimental value by approximately 60%. Besides being numerically much more efficient than its direct DFT counterpart, the DNN scheme has the advantage of being easily applicable to more sophisticated DFT approximations, such as meta-GGA and hybrid functionals, for which it would be hard to derive analytically the expression of the energy flux. We find in this way that a DNN model, trained on meta-GGA (SCAN) data, reduces the deviation from experiment of the predicted thermal conductivity by about 50%, leaving the question open as to whether the residual error is due to deficiencies of the functional, to a neglect of nuclear quantum effects in the atomic dynamics, or, likely, to a combination of the two.

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

Magnetic skyrmion diode: Unidirectional skyrmion motion via symmetry breaking of potential energy barriers

In this work, we realize a magnetic skyrmion diode operated by a unidirectional skyrmion transport that flows in only one direction, which is highly significant for information processing in spintronic and nanoelectronic devices. We easily control the skyrmion transport by engineering asymmetric shapes of geometric structures. Here, we present a simple method to describe the underlying mechanism behind the unidirectional skyrmion transport by characterizing the topography of potential energy surfaces from a purely geometric perspective. Our approach enables a deeper physical insight into skyrmion transport manipulation and efficient design of skyrmion-based devices in geometric structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Coherently coupled carrier and phonon dynamics in elemental tellurium probed by extreme ultraviolet transient absorption

The narrow band-gap semiconductor elemental tellurium (Te) has a unique electronic structure due to strong spin-orbit splitting and a lack of inversion symmetry of its helical lattice. Using broadband extreme ultraviolet core-level transient absorption, we simultaneously measure the coherently coupled photoinduced carrier and lattice dynamics at the Te N4,5 edge initiated by a few-cycle NIR pulse. Ultrafast excitation of carriers leads to a coherently excited A1 phonon oscillation and the generation of a hot carrier population distribution that oscillates in temperature, and the phonon excursion and hot carrier temperature are π out of phase with respect to each other. The depths of modulation suggest a significant coupling between the electronic and lattice degrees of freedom in Te. A long-lived shift of the absorption edge suggests an excited state of Te in another equilibrium potential energy surface that lives on the order of the carrier recombination timescale. The observed phonon-induced oscillations of the hot carriers are supportive of a change in the metallicity, whereby Te becomes more metallic with increasing phonon-induced displacement. Additionally, near the Fermi level we observe an energy-dependent phase of the displacive excitation of the A1 phonon mode. The discovery of coherent coupling between the lattice and hot carriers in Te provides the basis to investigate coherent interactions between spin and orbital degrees of freedom. The results spectrally and temporally resolve the correlation between photoexcited hot carriers and coherent lattice excitations, providing insight on the optical manipulation of the Te electronic structure at high carrier densities exceeding 1021cm−3.

Adelman, Jonah R↗
DOE OSTI2020

Identification of new transitions and levels in 163 Gd from $β$-decay studies

Background: Neutron-rich nuclei in the mass region around A = 160 have been and will continue to be of interest for the study of nuclear structure because of the rapid onset of deformation between 88 and 90 neutrons. The observation of detailed changes in nuclear structures within this mass region has provided and will continue to provide insight into the nuclear force. Purpose: Investigations of γ rays emitted following Eu 163 β -decay to Gd 163 have been performed for evaluation of the nuclear structure of Gd 163 . Method: Data were collected at the LeRIBSS station of the Holifield Radioactive Ion Beam Facility at Oak Ridge National Laboratory with an array of four Clover HPGe detectors for γ -rays and two plastic scintillators for β detection. The γ rays were identified as belonging to Gd 163 via mass selection and γ – γ – β , x-ray- γ , or γ – γ coincidences. Results: In total 107 new γ -ray transitions were observed in Gd 163 from 53 newly identified levels. Conclusions: We report the structure of Gd 163 has been identified for the first time. This structure has been evaluated in comparison to projected shell model, and potential energy surface calculations with good agreement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Efficient method for estimation of fission fragment yields of $\textit{r}$-process nuclei

Background: More than half of all the elements heavier than iron are made by the rapid neutron capture process (or $\textit{r}$ process). For very-neutron-rich astrophysical conditions, such at those found in the tidal ejecta of neutron stars, nuclear fission determines the $\textit{r}$-process endpoint, and the fission-fragment yields shape the final abundances of 110 ≤ $\textit{A}$ ≤ 170 nuclei. The knowledge of fission-fragment yields of hundreds of nuclei inhabiting very-neutron-rich regions of the nuclear landscape is thus crucial for the modeling of heavy-element nucleosynthesis. Purpose: In this study, we propose a model for the fast calculation of fission-fragment yields based on the concept of shell-stabilized prefragments defined with help of the nucleonic localization functions. Methods: To generate realistic potential-energy surfaces and nucleonic localizations, we apply Skyrme density-functional theory. In this work, the distribution of the neck nucleons among the two prefragments is obtained by means of a statistical model. Results: We benchmark the method by studying the fission yields of 178 Pt, 240 Pu, 254 Cf, and 254,256,258 Fm and show that it satisfactorily explains the experimental data. We then make predictions for 254 Pu and 290 Fm as two representative cases of fissioning nuclei that are expected to significantly contribute during the $\textit{r}$-process nucleosynthesis occurring in neutron-star mergers. Conclusions: The proposed framework provides an efficient alternative to microscopic approaches based on the evolution of the system in a space of collective coordinates all the way to scission. It can be used to carry out global calculations of fission-fragment distributions across the $\textit{r}$-process region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Isomeric and beta-decay spectroscopy of 173,174 Ho

β-decay spectroscopy of 173,174 Ho (Z = 67, N = 106,107) was conducted at Radioactive Isotope Beam Factory at RIKEN by using in-flight fission of a 345-MeV/u 238 U primary beam. A previously unreported isomeric state at 405 keV with half-life of 3.7(12) μs and a spin and parity of (3/2 + ) is identified in 173 Ho. Moreover, a new state with a spin and parity of 9 - was discovered in 174 Er. The experimental log ft values of 5.84(20) and 5.25(18) suggest an allowed-hindered β decay from the ground state of 174 Ho to the K π = 8 - isomeric state in 174 Er. Configuration-constrained potential energy surface (PES) calculations were performed and the predictions are in reasonable agreement with the experimental results.

150 ≤ A ≤ 189↗
DOE OSTI2020

Merging of transport theory with the time-dependent Hartree-Fock approach: Multinucleon transfer in U + U collisions

Multinucleon transfer mechanism in the collision of 238 U + 238 U system is investigated at E c.m. = 833 MeV in the framework of the quantal diffusion description based on the stochastic mean-field approach. Double cross sections σ (N, Z ) as a function of the neutron and proton numbers, the cross sections σ (Z ) and σ (A) as a function of the atomic numbers and the mass numbers are calculated for production of the primary fragments. The calculation indicates the 238 U + 238 U system may be located at an unstable equilibrium state at the potential energy surface with a slightly negative curvature along the β stability line on the (N, Z ) plane. Here, this behavior may lead to rather large diffusion along the β stability direction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

New transitions and levels for 163 Tb obtained from β-decay studies

Transitions in 163 Tb following β decay of 163 Gd were obtained as part of investigations of γ rays emitted following 163 Euβ decay to 163 Gd. In this work, detailed analysis of the low-energy structure of 163 Tb has been carried out with these data to expand previous β-decay studies and reactions studies of levels in 163 Tb. Data were collected at the LeRIBSS station of the Holifield Radioactive Ion Beam Facility at Oak Ridge National Laboratory with an array of four Clover HPGe detectors for γ rays and two plastic scintillators for β detection. The γ rays were identified as belonging to 163 Gd and 163 Tb via mass selection and γ-γ-β, γ-γ, or γ-x-ray coincidence. In total, 38 new γ-ray transitions were observed in 163 Tb from 15 newly identified levels and 12 previously identified levels. Potential energy surface calculations were performed which support a rigid prolate deformation. Previously identified unplaced transitions in 163 Tb have been placed within the level scheme of 163 Tb and additional states and transitions have been identified.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Proton and neutron exchange as a prelude to fusion at near-barrier energies

Systematic examination of fusion for 39,41,45,47 K + 28 Si and 36,44 Ar + 28 Si provides insight into the impact of neutron and proton exchange on fusion for nuclei at and near the N = 20 and N = 28 shells. Comparison of the reduced excitation functions reveals a marked difference between the behavior of open-shell and closed-shell systems. While coupled channels calculations provide a good description for the closed-shell nuclei, they significantly underpredict the fusion cross section for open-shell nuclei. Furthermore, the observed trends are examined in the context of a potential energy surface, including shell effects, and multinucleon exchange with consideration of Pauli blocking.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Level scheme of Nb 92 and observation of an oblate collective rotational band

Excited states of 92 Nb were investigated using the reactions 82 Se( 14 N, 4$n$) and 89 Y( 6 Li, 1 $p$2$n$). Based on these experimental results, the level scheme was extended abundantly. The structures of both the positive and negative parity states have been compared with shell model calculations using the SNET interaction. Potential energy surface calculation using the configuration-constrained method was performed to study a rotational-like band in 92 Nb, and possible configuration of the band head has been discussed. In conclusion, the calculation indicates that it may be an oblate collective rotational band.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Optimized structure and electronic band gap of monolayer GeSe from quantum Monte Carlo methods

Here, we have used highly accurate quantum Monte Carlo methods to determine the chemical structure and electronic band gaps of monolayer GeSe. Two-dimensional (2D) monolayer GeSe has received a great deal of attention due to its unique thermoelectric, electronic, and optoelectronic properties with a wide range of potential applications. Density functional theory (DFT) methods have usually been applied to obtain optical and structural properties of bulk and 2D GeSe. For the monolayer, DFT typically yields a larger band-gap energy than for bulk GeSe but cannot conclusively determine if the monolayer has a direct or indirect gap. Moreover, the DFT-optimized lattice parameters and atomic coordinates for monolayer GeSe depend strongly on the choice of approximation for the exchange-correlation functional, which makes the ideal structure-and its electronic properties-unclear. In order to obtain accurate lattice parameters and atomic coordinates for the monolayer, we use a surrogate Hessian-based parallel line search within diffusion Monte Carlo to fully optimize the GeSe monolayer structure. The DMC-optimized structure is different from those obtained using DFT, as are calculated band gaps. The potential energy surface has a shallow minimum at the optimal structure. This, combined with the sensitivity of the electronic structure to strain, suggests that the optical properties of monolayer GeSe are highly tunable by strain.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Photoluminescence spectra of point defects in semiconductors: Validation of first-principles calculations

Optically and magnetically active point defects in semiconductors are interesting platforms for the development of solid state quantum technologies. Their optical properties are usually probed by measuring photoluminescence spectra, which provide information on excitation energies and on the interaction of electrons with lattice vibrations. We present a combined computational and experimental study of photoluminescence spectra of defects in diamond and SiC, aimed at assessing the validity of theoretical and numerical approximations used in first-principles calculations, including the use of the Franck-Condon principle and the displaced harmonic oscillator approximation. We focus on prototypical examples of solid state qubits, the divacancy centers in SiC and the nitrogen-vacancy in diamond, and we report computed photoluminescence spectra as a function of temperature that are in very good agreement with the measured ones. As expected we find that the use of hybrid functionals leads to more accurate results than semilocal functionals. Interestingly our calculations show that constrained density functional theory (CDFT) and time-dependent hybrid DFT perform equally well in describing the excited state potential energy surface of triplet states; our findings indicate that CDFT, a relatively cheap computational approach, is sufficiently accurate for the calculations of photoluminescence spectra of the defects studied here. Finally, we find that only by correcting for finite-size effects and extrapolating to the dilute limit can one obtain a good agreement between theory and experiment. Our results provide a detailed validation protocol of first-principles calculations of photoluminescence spectra, necessary both for the interpretation of experiments and for robust predictions of the electronic properties of point defects in semiconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Excitation protocols for nonlinear phononics in bismuth and antimony

Here, we study the optical generation and control of coherent phonons in elemental bismuth (Bi) and antimony (Sb) using a classical equation of motion informed by first-principles calculations of the potential energy surface and the frequency-dependent macroscopic dielectric function along the zone-centered optical phonon coordinates. Using this approach, we demonstrate that phonons with the largest optomechanical couplings also have the strongest degree of anharmonicity among the zone-centered modes, a result of the broken-symmetry structural ground state of Bi and Sb. We show how this anharmonicity, explaining the light-induced phonon softening observed in experiments, prevents the application of standard phonon-amplification and annihilation protocols. We introduce a simple linearization protocol that extends the use of such protocols to the case of anharmonic phonons in broken-symmetry materials, and demonstrate its efficiency at high displacement amplitudes. Our formalism and results provide a path for improving optical control in nonlinear phononics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Ferroelectric phase transition in group-IV monochalcogenides from an equivariant machine learned force field

Group-IV monochalcogenides are a class of layered ferroelectric semiconductors that have demonstrated spontaneous intrinsic polarization above room temperature. Here, in this study, we use the multi-atomic cluster expansion (MACE) machine learning architecture to train and test a force field capable of modeling the structural properties and second-order ferroelectric-to-paraelectric phase transition in a Group-IV monochalcogenide, GeSe. The model captures the double-well potential energy surface associated with the onset of macroscopic polarization in bulk GeSe within 12.5 meV/atom, as well as near-equilibrium properties like the phonon dispersion. The development of this quantitatively accurate force field enables long-time molecular dynamics simulations, which predict the critical temperature of the ferroelectric-to-paraelectric phase transition in bulk GeSe to be T c = 600 K. This study demonstrates the capabilities of equivariant force-fields to accurately describe phenomena associated with structural symmetry breaking.

ferroelectricity↗
DOE OSTI2025

CuXASNet: Rapid and accurate prediction of copper L-edge x-ray absorption spectra using machine learning

In this work, we have developed CuXASNet, a dense neural network that predicts simulated Cu -edge x-ray absorption spectra (XAS) from atomic structures. Featurization of the Cu local environment is performed using a component of M3GNet, a graph neural network developed for predicting the potential energy surface. CuXASNet is trained on simulated spectra from FEFF9 at the multiple scattering level of theory, and can predict the and edges for Cu sites to quantitative accuracy. To validate our approach, we compare 14 experimental spectra extracted from the literature with the predictions of CuXASNet. The agreement of CuXASNet with experiments is shown by an average mean absolute error of 0.125 and an average Spearman's correlation coefficient of 0.891, which is comparable to FEFF9's values of 0.131 and 0.898 for the same metrics. As such, CuXASNet can rapidly predict a large number of -edge XAS spectra at the same accuracy as FEFF9 simulations. This can be used as a drop-in replacement for multiple scattering codes for fast screening of candidate atomic structure models of a measured system. This model establishes a general framework for Cu XAS prediction, and can be extended to more computationally expensive levels of theory and to other transition metal edges.

36 MATERIALS SCIENCE↗