SEARCH · Engineering Papers

Results for “time dependant density functional theory”

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.

SourceResearch connectionFrom yearThrough year

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

At least 505 records · Page 28

All sources NASA NTRS DOE OSTI

DOE OSTI2022

An Experimental Investigation of Supercritical Methane Injection Characteristics in a CO2 Environment

Abstract Clean energy generation is gaining significant attention from industries, academia, and governments across the globe. The Allam cycle is one such technology that has been under focus due to its efficiency, environmental friendliness, and economics. This is a direct-fired cycle operating at supercritical conditions using carbon dioxide as a working fluid. Fuel or oxidizer jet mixing with CO2 is a vital phenomenon that governs combustion efficiency, and it is not well understood for the Allam cycle conditions. This paper experimentally investigated the jet characteristics of a methane jet injected into a subcritical to supercritical carbon dioxide environment. A wide range of injection pressures and temperatures were targeted between subcritical to supercritical conditions. Unlike previous studies, the current work focused on injecting lower-density (methane) jets into higher-density (carbon dioxide) environments. Schlieren imaging and methane absorption measurements were simultaneously performed with a CMOS high-speed camera and a 3.39 μm infrared laser. Specifically, we looked at the classical injection parameter of jet spreading angle, which was classically established to be mainly a density ratio function. Here, the jet cone angle was obtained from the postprocessed schlieren imaging. The jet cone angle is a critical characteristic parameter that describes the entrainment rate in a jet; thus, it is a crucial parameter in understanding the nature of the jet. The laser measurements were only used as an additional check to confirm the entry time of methane into the chamber filled with carbon dioxide. Notably, this paper makes a detailed comparison between the jet cone angles of jets with a density ratio. The result showed that the classical correlations, such as Abramovich's theory applied to submerged turbulent gas jets developed for low-density ratio jets, were unsuitable for higher-density ratio jets. It was also observed that the divergence angles were dependent not only on density ratio but also on other parameters such as pressure ratios and reduced pressures.

Engineering↗

DOE OSTI2026

Characterizing Defect Dynamics in Silicon Carbide Using Symmetry-Adapted Collective Variables and Machine Learning Interatomic Potentials

Silicon carbide (SiC) divacancies are attractive candidates for spin-defect qubits possessing long coherence times and optical addressability. The high activation barriers associated with SiC defect formation and motion pose challenges for their study by first-principles molecular dynamics. In this work, we develop and deploy machine learning interatomic potentials (MLIPs) to accelerate defect dynamics simulations while retaining ab initio accuracy. We employ an active learning strategy comprising symmetry-adapted collective variable discovery and enhanced sampling to compile configurationally diverse training data, calculation of energies and forces using density functional theory (DFT), and training of an E(3)-equivariant MLIP based on the Allegro model. Here, the trained MLIP reproduces DFT-level accuracy in defect transition activation free energy barriers, enables the efficient and stable simulation of multidefect 216-atom supercells, and permits an analysis of the temperature dependence of defect thermodynamic stability and formation/annihilation kinetics to propose an optimal annealing temperature to maximally stabilize VV divacancies.

Computer simulations↗

DOE OSTI2022

Understanding the Reactivity, Selectivity, and Deactivation of Frustrated Lewis Pairs for Semihydrogenation of Acetylene

A good catalyst for semihydrogenation of alkynes must preclude both over-hydrogenation of alkene to alkane and isomerization to the other alkene isomer. In addition, it should balance the trade-off between selectivity and activity. In 2013, the Repo and Pápai groups reported a frustrated Lewis pair (FLP) (1-NMe 2 -2-B(C 6 F 5 ) 2 -C 6 H 4 ), 1, which is a metal-free catalyst and for the first time shows excellent reactivity for the hydrogenation of internal alkynes. However, it is unreactive for terminal alkynes. In this work, we have designed 13 FLPs, a–m, based on 1 by varying the Lewis base site with N and P and the Lewis acid site with B, Al, Ga, and In and replacing pentafluorophenyl with 1,3,5-trifluorophenyl, phenyl, or trifluoromethyl. We apply density functional theory to study the activity, selectivity, and deactivation of FLP 1-m for acetylene semihydrogenation. The catalytic cycle consists of three steps: (1) alkyne insertion, (2) H 2 heterolysis, and (3) intramolecular protonation. We found the activity does not change much by the modification of bulky ligands, while it decreases with the direct replacement of LA and LB sites. The overall activity depends on steps 1 and 3, which are, respectively, positively and negatively linear correlated with the charge of the Lewis acid site. Most of the FLPs in this work show comparable or better selectivity for semihydrogenation of acetylene than 1. FLP deactivation is due to the strong binding of acetylene and the elimination of electron-withdrawing bulky ligands at the preactivated catalyst rather than at activated catalysts. Furthermore, taking the selectivity and stability of FLPs into account, we predict d and k are potentially active for terminal alkynes.

02 PETROLEUM↗

DOE OSTI2022

Nonequilibrium DMFT + CPA for correlated disordered systems

We present a solution for the nonequilibrium dynamics of an interacting disordered system. The approach adapts the combination of the equilibrium dynamical mean-field theory and the equilibrium coherent potential approximation methods to the nonequilibrium many-body formalism, using the Kadanoff-Baym-Keldysh complex time contour, for the dynamics of interacting disordered systems away from equilibrium. We use our time domain solution to obtain the equilibrium density of states of the disordered interacting system described by the Anderson-Hubbard model, bypassing the necessity for the cumbersome analytical continuation process. We further apply the nonequilibrium solution to the interaction quench problem for an isolated disordered system. Here, the interaction is abruptly changed from zero (noninteracting system) to another constant (finite) value at which it is subsequently kept. We observe via the time dependence of the potential, kinetic, and total energies the effect of disorder on the relaxation of the system as a function of final interaction strength. Finally, the real-time approach has the potential to shed light on the fundamental role of disorder in the nonequilibrium dynamics of interacting quantum systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Scaling and adiabaticity in a rapidly expanding gluon plasma

In this work we aim to gain qualitative insight on the far-from-equilibrium behavior of the gluon plasma produced in the early stages of a heavy-ion collision. It was recently discovered [1] that the distribution functions of quarks and gluons in QCD effective kinetic theory (EKT) exhibit self-similar “scaling” evolution with time-dependent scaling exponents long before those exponents reach their pre-hydrodynamic fixed-point values. In this work we shed light on the origin of this time-dependent scaling phenomenon in the small-angle approximation to the Boltzmann equation. We first solve the Boltzmann equation numerically and find that time-dependent scaling is a feature of this kinetic theory, and that it captures key qualitative features of the scaling of hard gluons in QCD EKT. We then proceed to study scaling analytically and semi-analytically in this equation. We find that an appropriate momentum rescaling allows the scaling distribution to be identified as the instantaneous ground state of the operator describing the evolution of the distribution function, and the approach to the scaling function is described by the decay of the excited states. That is to say, there is a frame in which the system evolves adiabatically. Furthermore, from the conditions for adiabaticity we can derive evolution equations for the time-dependent scaling exponents. In addition to the known free-streaming and BMSS fixed points, we identify a new “dilute” fixed point when the number density becomes small before hydrodynamization. Corrections to the fixed point exponents in the small-angle approximation agree quantitatively with those found previously in QCD EKT and arise from the evolution of the ratio between hard and soft scales.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

A collective elastic fluctuation mechanism for decoupling and stretched relaxation in glassy colloidal and molecular liquids

We propose a microscopic theory for the decoupling of self-diffusion and structural relaxation in glass-forming liquids within the Elastically Collective Nonlinear Langevin Equation (ECNLE) activated dynamics framework. Our central hypothesis is that the heterogeneity relevant to this problem is static fluctuations of local density on the scale of 3–4 particle diameters and how this changes local packing correlations. These fluctuations modify the degree of dynamical cage expansion that mechanistically couples intracage large amplitude hopping and longer range collective elasticity in ECNLE theory. Decoupling only emerges in the deeply supercooled regime where the strongly temperature dependent elastic barrier becomes non-negligible relative to its noncooperative local analog. The theory makes predictions for various aspects of the decoupling phenomenon, including apparent fractional power law Stokes-Einstein behavior, that appear to be consistent with experiments and simulations on hard sphere fluids and molecular liquids. Of central importance is a microscopic connection between the barrier fluctuation variance and most probable barrier height. Sensible results are also obtained for the nonexponential stretching of a generic relaxation time correlation function and its temperature evolution. Nonuniversality can arise from the relative importance of the local and collective barriers (related to fragility) and the precise magnitude of the length scale that defines the transition from local cage to elastic physics. Comparison is made with a traplike model based on a Gaussian distribution of barriers.

Xie, Shi-Jie (ORCID:0000000333174720)↗

DOE OSTI2022

First-Principles Calculations of Thermoelectric Transport Properties of Quaternary and Ternary Bulk Chalcogenide Crystals

Chalcogenide crystals have a wide range of applications, especially as thermoelectric materials for energy conversion. Thermoelectric materials can be used to generate an electric current from a temperature gradient based on the Seebeck effect and based on the Peltier effect, and they can be used in cooling applications. Using first-principles calculations and semiclassical Boltzmann theory, we have computed the Seebeck coefficient, electrical conductivity, electronic thermal conductivity, power factor, and figure of merit of 30 chalcogenide crystals. A Quantum Espresso package is used to calculate the electronic properties and locate the Fermi level. The transport properties are then calculated using the BoltzTraP code. The 30 crystals are divided into two groups. The first group has four crystals with quaternary composition (A 2 BCQ 4 ) (A = Tl; B = Cd, Hg; C = Si, Ge, Sn; Q = S, Se, Te). The second group contains 26 crystals with the ternary composition (A’B’Q 2 ) (A’ = Ag, Cu, Au, Na; B’ = B, Al, Ga, In; Q = S, Se, Te). Among these 30 chalcogenide crystals, the results for 11 crystals: Tl 2 CdGeSe 4 , Tl 2 CdSnSe 4 , Tl 2 HgSiSe 4 , Tl 2 HgSnS 4 , AuBSe 2 , AuBTe 2 , AuAlTe 2 , AuGaTe 2 , AuInTe 2 , AgAlSe 2 , and AgAlTe 2 are revealed for the first time. In addition, temperature-dependent transport properties of pure and doped AgSbSe 2 and AgSbTe 2 crystals with dopant compositions of AgSb 0.94 Cd 0.06 Te 2 and AgSbTe 1.85 Se 0.15 were explored. These results provide an excellent database for bulk chalcogenides crucial for a wide range of potential applications in renewable energy fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Self-Consistent Field Methods for Excited States in Strong Magnetic Fields: A Comparison Between Energy - and Variance-based Approaches

An efficient implementation of geometrical derivatives at the Hartree–Fock (HF) and current-density functional theory (CDFT) levels is presented for the study of molecular structure in strong magnetic fields. The required integral derivatives are constructed using a hybrid McMurchie–Davidson and Rys quadrature approach, which combines the amenability of the former to the evaluation of derivative integrals with the efficiency of the latter for basis sets with high angular momentum. In addition to its application to evaluating derivatives of four-center integrals, this approach is also applied to gradients using the resolution-of-the-identity approximation, enabling efficient optimization of molecular structure for many-electron systems under a strong magnetic field. The CDFT contributions have been implemented for a wide range of density functionals up to and including the meta-GGA level with current-density dependent contributions and (range-separated) hybrids for the first time. Illustrative applications are presented to the OH and benzene molecules, revealing the rich and complex chemistry induced by the presence of an external magnetic field. Challenges for geometry optimization in strong fields are highlighted, along with the requirement for careful analysis of the resulting electronic structure at each stationary point. The importance of correlation effects is examined by comparison of results at the HF and CDFT levels. The present implementation of molecular gradients at the CDFT level provides a cost-effective approach to the study of molecular structure under strong magnetic fields, opening up many new possibilities for the study of chemistry in this regime.

Comparison↗

DOE OSTI2023

A time-dependent momentum-resolved scattering approach to core-level spectroscopies

While new light sources allow for unprecedented resolution in experiments with X-rays, a theoretical understanding of the scattering cross-section is lacking. In the particular case of strongly correlated electron systems, numerical techniques are quite limited, since conventional approaches rely on calculating a response function (Kramers-Heisenberg formula) that is obtained from a perturbative analysis of scattering processes in the frequency domain. This requires a knowledge of a full set of eigenstates in order to account for all intermediate processes away from equilibrium, limiting the applicability to small tractable systems. In this work, we present an alternative paradigm, recasting the problem in the time domain and explicitly solving the time-dependent Schrödinger equation without the limitations of perturbation theory: a faithful simulation of the scattering processes taking place in actual experiments, including photons and core electrons. We show how this approach can yield the full time and momentum resolved Resonant Inelastic X-Ray Scattering (RIXS) spectrum of strongly interacting many-body systems. We demonstrate the formalism with an application to Mott insulating Hubbard chains using the time-dependent density matrix renormalization group method, which does not require a priory knowledge of the eigenstates and can solve very large systems with dozens of orbitals. This approach can readily be applied to systems out of equilibrium without modification and generalized to other spectroscopies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Random Phase Approximation Correlation Energy Using Real-Space Density Functional Perturbation Theory

We present a real-space method for computing the random phase approximation (RPA) correlation energy within Kohn–Sham density functional theory, leveraging the low-rank nature of the frequency-dependent density response operator. In particular, we employ a cubic-scaling formalism based on density functional perturbation theory that circumvents the calculation of the response function matrix, instead relying on the ability to compute its product with a vector through the solution of the associated Sternheimer linear systems. We develop a large-scale parallel implementation of this formalism using the subspace iteration method in conjunction with the spectral quadrature method while employing the Kronecker product-based method for the application of the Coulomb operator and the conjugate orthogonal conjugate gradient method for the solution of the linear systems. We demonstrate convergence with respect to key parameters and verify the method’s accuracy by comparing with plane-wave results. We show that the framework achieves good strong scaling to many thousands of processors, reducing the time to solution for a lithium hydride system with 128 electrons to around 150 s on 4608 processors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Nonharmonic contributions to the high-temperature phonon thermodynamics of Cr

Here, phonon densities of states (DOSs) of body-centered cubic chromium were measured by time-of-flight inelastic neutron scattering at temperatures up to 1493 K. Density functional theory calculations with both quasiharmonic (QH) and anharmonic (AH) methods were performed at temperatures above the Néel temperature. Features in the phonon DOSs decrease in energy (soften) substantially with temperature. A Born–von Kármán analysis using fits to the experimental DOSs reveals a softening of almost 17% of the high-transverse phonon branch between 330 and 1493 K. The low-transverse branch changes by approximately half this amount. The AH calculations capture the observed behavior of the two transverse phonon branches, but the QH calculations give some inverted trends. Vibrational entropies from phonons and electrons are obtained, and their sum is in excellent agreement with the entropy of chromium obtained by calorimetry, indicating that above 330 K, no explicit temperature-dependent magnetic contributions are necessary.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1978

A magnetohydrodynamic theory for the lunar response to time variations in a spatially uniform ambient magnetic field

An analytic theory based on principles of magnetohydrodynamics is derived for interpreting the lunar magnetic response to magnetic field fluctuations in the lobes of the geomagnetic tail. Only the case of magnetic perturbations parallel to a spatially uniform ambient magnetic field is treated. The theoretical frequency-dependent MHD transfer function depends on both the lunar electrical conductivity profile and the Alfven speed of the tail lobe plasma. For circular frequencies much less than the Alfven speed/lunar radius, the MHD transfer function reduces to the vacuum transfer function, while at higher frequencies the finite propagation speed of magnetic disturbances affects the response. These characteristics of the MHD transfer function are in qualitative agreement with observation. A carefully selected set of magnetic field transient events could be used together with this theory to infer lunar electrical conductivity profiles and the tail lobe plasma density (via its dependence on the Alfven speed).

Hood, L. L.↗

DOE OSTI2023

Thermodynamic and electron transport properties of Ca 3 Ru 2 O 7 from first-principles phonon calculations and Boltzmann transport theory

Here, this work demonstrates a first-principles-based approach to obtaining finite temperature thermal and electronic transport properties which can be employed to model and understand mesoscale structural evolution during electronic, magnetic, and structural phase transitions. A computationally tractable model was introduced to estimate electron relaxation time and its temperature dependence. The model is applied to Ca 3 Ru 2 O 7 with a focus on understanding its electrical resistivity across the electronic phase transition at 48 K. A quasiharmonic phonon approach to the lattice vibrations was employed to account for thermal expansion while the Boltzmann transport theory including spin-orbit coupling was used to calculate the electron-transport properties, including the temperature dependence of electrical conductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

VERITAS : A density-functional theory-based multiband kinetic model for understanding x-ray spectroscopy of dense plasmas

X-ray spectroscopy has long been a powerful diagnostic tool for hot, dilute plasmas, providing insights into plasma conditions by measuring line shifts and broadenings of atomic transitions. The technique critically depends on the accuracy of atomic physics models used to interpret spectroscopic measurements for inferring plasma properties such as free-electron density and temperature. Over the past decades, the atomic and plasma physics communities have developed robust atomic physics models to account for various processes in hot, dilute classical plasmas. While these models have been successful in that regime, their applicability becomes uncertain when interpreting x-ray spectroscopy experiments of above-solid-density plasmas. Given that finite-temperature density-functional theory (DFT) offers a more accurate description of dense plasma environments, we present the development of a DFT-based multi-band kinetic model, VERITAS, designed to improve the interpretation of x-ray spectroscopic measurements in high-density plasmas produced by laser-driven spherical implosions. This work details the VERITAS model and its application to both time-integrated and time-resolved x-ray spectra from implosion experiments on OMEGA. The advantages and limitations of the VERITAS model will also be discussed, along with potential directions for advancing x-ray spectroscopy of dense and superdense plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Predicting magnetic properties of single-molecule magnets from self-interaction-free density-functional theory (Final Report)

In this project we investigated electronic structure of an intermediate-sized copper-based molecule and magnetic and hyperfine properties of several small non-magnetic and magnetic molecules including transition-metal elements by applying self-interaction corrections to density functional theory (DFT). In a sufficient number of cases, DFT-calculated electronic structure and magnetic properties of single-molecule magnets qualitatively differ from corresponding experimental data. This is partly due to self-interacting electrons within the DFT formalism. Recently, an efficient method to correct the self-interactions was proposed, i.e., Fermi-Lowdin prbital (FLO) based self-interaction corrected (SIC) methodology, within DFT. Henceforth, this method is referred to as FLO-SIC method which exists in FLOSIC code. We used this FLO-SIC method for our studies of electronic structure and magnetic and hyperfine properties of small magnetic molecules and non-magnetic molecules. Our study will provide insight into predictions of magnetic properties of single-molecule magnets where self-interaction corrections play a critical role. There are two components of this project. In the first work, we studied the electronic structure of a planar mononuclear Cu-based molecule in two oxidation states, using DFT with the FLO-SIC method. We chose this system because it is small enough and it includes a transition metal element. We found that the standard FLO-SIC method takes too much compute time even for the small transition-metal molecule and so we slightly modified the method in order to expedite the process. In the dianionic state, we found that the FLO-SIC spin density agrees quantitatively with accurate quantum chemistry methods, while DFT spin density without self-interaction corrections are severely deviated from the quantum chemistry methods. We also showed that the energy gap between the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) of the dianionic state is larger than that of the monoanionic state. This result is consistent with experimental data. In the second work, we investigated how the interaction between the electron spin and nuclear spin is affected by electron self-interactions within small non-magnetic and magnetic molecules. Such an interaction is called hyperfine interaction. For molecules without significant orbital angular momentum, the hyperfine interaction consists of Fermi contact and dipolar interaction terms. Since the Fermi contact term depends on electron spin density at the nuclear site, it would be highly affected by self-interaction corrections. Therefore, we calculated the hyperfine interaction for the small molecules using DFT with the slightly modified expedited FLO-SIC method which was obtained in the first work, and compared the results to experimental data and DFT calculations without self-interaction corrections. We found significant improvement of the Fermi contact term computed using the FLO-SIC method for small magnetic molecules. Overall, the first and second work provided positive outlook of application of the FLO-SIC method to magnetic molecules and systems including transition-metal elements.

36 MATERIALS SCIENCE↗

NASA NTRS2021

Constraints on Cosmic Strings Using Data from the Third Advanced LIGO-Virgo Observing Run

We search for gravitational-wave signals produced by cosmic strings in the Advanced LIGO and Virgo full O3 dataset. Search results are presented for gravitational waves produced by cosmic string loop features such as cusps, kinks, and, for the first time, kink-kink collisions. A template-based search for short-duration transient signals does not yield a detection. We also use the stochastic gravitational-wave background energy density upper limits derived from the O3 data to constrain the cosmic string tension Gμ as a function of the number of kinks, or the number of cusps, for two cosmic string loop distribution models. Additionally, we develop and test a third model that interpolates between these two models. Our results improve upon the previous LIGO–Virgo constraints on Gμ by 1 to 2 orders of magnitude depending on the model that is tested. In particular, for the one-loop distribution model, we set the most competitive constraints to date: Gμ ≲ 4×10−15. In the case of cosmic strings formed at the end of inflation in the context of grand unified theories, these results challenge simple inflationary models.

J B Camp↗

DOE OSTI2023

Differential charge radii: Proton-neutron interaction effects

The analysis of self-consistency and proton-neutron interaction effects in the buildup of differential charge radii were carried out in covariant density functional theoretical calculations without pairing interaction on the example of selected configurations in the Pb isotopes. The proton-neutron interaction of neutron(s) added to the neutron 𝑁 = 126 core and the protons forming the 𝑍 = 82 proton core is responsible for a major contribution to the buildup of differential charge radii. It depends on the products of proton and neutron wave functions and thus on their nodal structure. This interaction leads to a redistribution of single-particle density of occupied proton states which in turn modifies the charge radii. Furthermore, the microscopic origin of this redistribution and its consequences for differential charge radii were investigated for the first time. Self-consistency effects affecting the shape of proton potential, total proton densities and the energies of the single-particle proton states provide only minor contribution to differential charge radii.

190 ≤ A ≤ 219↗

DOE OSTI2022

Optical and magneto-optical properties of ferromagnetic monolayer CrBr 3 : A first-principles G W and G W plus Bethe-Salpeter equation study

The discovery of atomically thin two-dimensional (2D) magnetic semiconductors has triggered enormous research interest recently. In this paper, we use first-principles many-body perturbation theory to study a prototypical 2D ferromagnetic semiconductor, monolayer chromium tribromide (CrBr 3 ). With broken time-reversal symmetry, spin-orbit coupling, and excitonic effects included through the full-spinor GW and GW plus Bethe-Salpeter equation (GW-BSE) methods, we compute the frequency-dependent layer polarizability tensor and dielectric function tensor that govern the optical and magneto-optical (MO) properties. In addition, we provide a detailed theoretical formalism for simulating magnetic circular dichroism, MO Kerr effect, and Faraday effect, demonstrating the approach with monolayer CrBr 3 . Due to reduced dielectric screening in 2D and the localized nature of the Cr d orbitals, we find strong self-energy effects on the quasiparticle band structure of monolayer CrBr 3 that give a 3.8 eV indirect bandgap. Also, excitonic effects dominate the low-energy optical and MO responses in monolayer CrBr 3 where a large exciton binding energy of 2.3 eV is found for the lowest bright exciton state with excitation energy at 1.5 eV. We further find that the MO signals demonstrate strong dependence on the excitation frequency and substrate refractive index. Furthermore, our theoretical framework for modeling optical and MO effects could serve as a powerful theoretical tool for future study of optoelectronic and spintronics devices consisting of van der Waals 2D magnets.

36 MATERIALS SCIENCE↗