Search NASA⌕ Search

SEARCH · Search NASA

Results for “Classical 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.

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 OSTI2024

VOTCA: multiscale frameworks for quantum and classical simulations in soft matter

Many physical phenomena in liquids and soft matter are multiscale by nature and can involve processes with quantum and classical degrees of freedom occurring over a vast range of length- and timescales. Examples range from structure formation processes of complex polymers or even polymer blends (Svaneborg & Everaers, 2023) on the classical side to charge and energy transport and conversion processes (Lee et al., 2019) involving explicit electronic and, therefore, quantum information. The Versatile Object-oriented Toolkit for Coarse-graining Applications (VOTCA) provides multiscale frameworks built on a comprehensive set of methods for the development of classical coarse-grained potentials (VOTCA-CSG) as well as state-of-the art excited state electronic structure methods based on density-functional and many-body Green’s function theories, coupled in mixed quantum-classical models and used in kinetic network models (VOTCA-XTP).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Structure and Dynamics of CO 2 at the Air–Water Interface from Classical and Neural Network Potentials

The accurate description of the structure and dynamics of CO 2 at the instantaneous air–water interface, along with the effects of surface fluctuations on the CO 2 -transport processes, is essential for the development of negative emission technologies aimed at minimizing climate change. In this study, we performed molecular dynamics simulations of CO 2 at the air–water interface using neural network potentials (NNPs) trained on ab initio data generated through density-functional-theory-based molecular dynamics simulations. We compared these results with classical force fields to assess their performance in modeling interfacial CO 2 behavior. Our findings revealed that the asymmetric interactions, coupled with thermal surface fluctuations at the air–water interface, significantly influence CO 2 transport into the aqueous phase. The simulations demonstrate that classical force fields underestimate both the free energy of CO 2 transport and the strength of its interactions at the interface compared with the neural network potentials. In conclusion, the free energy and the interfacial dynamics of CO 2 are primarily influenced by the distribution of water within the instantaneous interfacial water layer, responsible for creating an asymmetric intermolecular interaction environment within the interfacial region.

Ab initio molecular dynamics↗
DOE OSTI2023

Extrinsic n-type doping of Cd 3 As 2 thin films

Cd 3 As 2 provides an excellent platform for studying the physics of three-dimensional Dirac semimetals due to its stability as well as its compatibility with thin film growth. Crystals made using both bulk and thin film synthesis are unintentionally doped n-type, and other than introducing Zn to reduce the carrier concentration, no efforts have been reported to alter this intrinsic doping without major changes to the band structure. Here, group VI elements Te and Se are introduced during epitaxy to increase the electron concentration of the films. Starting from an unintentionally doped electron concentration of 1–2 × 10 17 cm -3 , concentrations of up to 3 × 10 18 cm -3 are achieved. Analysis of Shubnikov–de Haas oscillations reveals good agreement in calculated effective mass and Fermi velocity of highly doped films with unintentionally doped single crystals with similar electron concentrations. The density functional theory is also performed to study the effects of group VI substitutions and confirms no strong perturbations in the electronic structure. This work ultimately demonstrates tunability in the carrier concentration using extrinsic dopants without substantial changes in the band structure, allowing for intentional design of Fermi-level position for device applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Absence of High-Pressure Ground-State Reentrant Ferroelectricity in PbTiO 3

We study ferroelectricity in the classic perovskite ferroelectric PbTiO 3 to high pressures with density functional theory (DFT) and experimental diamond-anvil techniques. We use second harmonic generation spectroscopy to detect lack of inversion symmetry. Consistent with early understanding and experiments, we find that ferroelectricity disappears at moderate pressures. However, DFT computations show that the disappearance arises from the overtaking of zone boundary instabilities, and not from the squeezing out of the off-centering ferroelectric displacements with pressure, as previously thought. Moreover, at high pressures the distorted perovskite phases are metastable with respect to a new dense centrosymmetric postperovskite phase with P⁢2 1 /m symmetry and 8-coordinated Ti, which becomes stable at about 70 GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Role of water model on ion dissociation at ambient conditions

In this work, we study ion pair dissociation in water at ambient conditions using a combination of classical and ab initio approaches. The goal of this study is to disentangle the sources of discrepancy observed in computed potentials of mean force. In particular, we aim to understand why some models favor the stability of solvent-separated ion pairs vs contact ion pairs. We found that some observed differences can be explained by non-converged simulation parameters. However, we also unveil that for some models, small changes in the solution density can have significant effects on modifying the equilibrium balance between the two configurations. We conclude that the thermodynamic stability of contact and solvent-separated ion pairs is very sensitive to the dielectric properties of the underlying simulation model. In general, classical models are very robust in providing a similar estimation of the contact ion pair stability, while this is much more variable in density functional theory-based models. The barrier to transition from the solvent-separated to contact ion pair is fundamentally dependent on the balance between electrostatic potential energy and entropy. This reflects the importance of water intra- and inter-molecular polarizability in obtaining an accurate description of the screened ion–ion interactions.

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

Sensitivity of time-dependent density functional theory to initial conditions

Time-dependent density-functional theory is mathematically formulated through nonlinear coupled time-dependent three-dimensional partial differential equations, and it is natural to expect a strong sensitivity of its solutions to variations of the initial conditions, akin to the butterfly effect ubiquitous in classical dynamics. Since the Schrödinger equation for an interacting many-body system is, however, linear and mathematically the exact equations of the density-functional theory reproduce the corresponding one-body properties, it would follow that the Lyapunov exponents are also vanishing within a density-functional theory framework. Whether for realistic implementations of the time-dependent density-functional theory the question of the absence of the butterfly effect and whether the dynamics provided is indeed a predictable theory was never discussed. At the same time, since the time-dependent density-functional theory is a unique tool allowing us to study the nonequilibrium dynamics of strongly interacting many-fermion systems, the question of predictability of this theoretical framework is of paramount importance. Here our analysis, for a number of quantum superfluid many-body systems (unitary Fermi gas, nuclear fission, and heavy-ion collisions) with a classical equivalent number of degrees of freedom O(10 10 ) and larger, suggests that its maximum Lyapunov exponents are negligible for all practical purposes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

The role of surfaces in flexoelectricity

We develop the relationship between the strain derivative of the mean-inner potential and surface contributions to flexoelectricity, identifying the true surface-specific component of the flexoelectric response of finite samples. Density functional theory calculations on a range of experimentally observed, low energy SrTiO 3 , MgO, and Si surfaces demonstrate that the mean-inner potential and its contributions to flexoelectricity are sensitive to small differences in surface structure, chemistry, and adsorbates. We also introduce a method to estimate mean-inner potential contributions to flexoelectricity using electron scattering factors and use this approximation to predict total flexoelectric responses for a variety of insulators. Strategies to experimentally disentangle bulk and surface flexoelectric terms are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Light–matter interaction of a molecule in a dissipative cavity from first principles

Cavity-mediated light–matter coupling can dramatically alter opto-electronic and physico-chemical properties of a molecule. Ab initio theoretical predictions of these systems need to combine non-perturbative, many-body electronic structure theory-based methods with cavity quantum electrodynamics and theories of open-quantum systems. Here, we generalize quantum-electrodynamical density functional theory to account for dissipative dynamics of the cavity and describe coupled cavity–single molecule interactions in the weak-to-strong-coupling regimes. Specifically, to establish this generalized technique, we study excited-state dynamics and spectral responses of benzene and toluene under weak-to-strong light–matter coupling. By tuning the coupling, we achieve cavity-mediated energy transfer between electronically excited states. This generalized ab initio quantum-electrodynamical density functional theory treatment can be naturally extended to describe cavity-mediated interactions in arbitrary electromagnetic environments, accessing correlated light–matter observables and thereby closing the gap between electronic structure theory, quantum optics, and nanophotonics.

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

Deterministic modeling of hybrid nonlinear effects in epsilon-near-zero thin films

In nonlinear optics, significant effort is concentrated on improving the strength and efficiency of interactions; however, experimentally investigating nonlinear materials is a complex, time-consuming, and costly investment. Moreover, it is often challenging to isolate, study, and optimize material parameters in an experiment due to complexities in the growth process. Recently, epsilon-near-zero materials have received a great deal of attention as promising nonlinear optical materials, but like many up-and-coming materials, the ability to explore and optimize their properties has been challenging. Here, we establish a framework to rapidly evaluate the performance of nonlinear epsilon-near-zero materials for both inter- and intraband effects in silico, requiring only an energy-momentum (E-k) diagram, linear optical properties, and experimental conditions. Measured nonlinear reflection and transmission in gallium-doped zinc oxide films are compared to the numerical framework for both intra- and interband excitation to verify accuracy across wavelength and irradiance while two figures of merit (FoMs) are introduced to quickly evaluate the performance of films without a full numerical framework. This capability is used to predict the performance of highly doped gallium nitride, cadmium oxide, zinc oxide, and indium tin oxide films, and efficient intra- and interband operation conditions are identified. Through this numerical framework and the FoMs, the exploration of unstudied epsilon-near-zero materials is enabled without the need for a nonlinear experiment, thereby accelerating the search for more efficient nonlinear materials and excitation conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Ab Initio Transport Calculations: From Normal to Superconducting Current

Applying the Bogoliubov-de Gennes equations with density-functional theory, it is possible to formulate first principles description of current-phase relationships in superconducting/normal (magnetic)/superconducting tri-layers. Such structures are the basis for the superconducting analog of Magnetoresistive random access memory devices (JMRAM). In a recent paper1 we presented results from the first attempt to formulate such a theory, applied to the Nb/Ni/Nb trilayers. In the present work we provide computational details, explaining how to construct key ingredient (scattering matrices SN ) in a framework of linear muffin-tin orbitals (LMTO).

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗
DOE OSTI2021

The role of precursor decomposition in the formation of samarium doped ceria nanoparticles via solid-state microwave synthesis

The impact on the final morphology of ceria (CeO 2 ) nanoparticles made from different precursors (commercial: cerium acetate/nitrate) and in house: cerium tri(methylsilyl)amide (Ce-TMSA)) via a microwave solid state reaction has been determined. In all instances, powder X-ray diffraction indicated that the cubic fluorite CeO 2 phase (PDF# 04–004-9150, with the space group Fm-3 m) had formed. Scanning electron microscopy (SEM) images revealed spherical nanoparticles were produced from the Ce-TMSA precursor. The commercial acetate and nitrate precursors produced particles with irregular morphology. The roles of the precursor decomposition and binding energy in the synthesis of the nanocrystals with various morphologies, as well as a possible growth mechanism, were evaluated based on experimental and computational data. The formation of spherical shaped nanoparticles was determined to be due to the preferential single-step decomposition of the Ce-TMSA as well as the low activation energy to overcome decomposition. Due to the complicated decomposition of the commercial precursors and high activation energy the resulting particles adopted an irregular morphology. Highly uniform samarium doped ceria (Sm x Ce 1-x O 2-δ ) nanospheres were also synthesized from Ce-TMSA and samarium tri(methylsilyl)amide (Sm-TMSA). The effects of reaction time and temperature, on the final morphology were observed through SEM. The rapid single-step decomposition of TMSA-based precursors as observed through thermogravimetric analysis (TGA) and confirmed through the calculation of potential energy surfaces and binding energies from density functional theory (DFT) calculations, indicated that nanoparticle formation follows LaMer’s classical nucleation theory.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Observation of a highly conductive warm dense state of water with ultrafast pump–probe free-electron-laser measurements

The electrical conductivity of water under extreme temperatures and densities plays a central role in modeling planetary magnetic fields. Experimental data are vital to test theories of high-energy-density water and assess the possible development and presence of extraterrestrial life. These states are also important in biology and chemistry studies when specimens in water are confined and excited using ultrafast optical or free-electron lasers (FELs). Here we utilize femtosecond optical lasers to measure the transient reflection and transmission of ultrathin water sheet samples uniformly heated by a 13.6 nm FEL approaching a highly conducting state at electron temperatures exceeding 20 000 K. The experiment probes the trajectory of water through the high-energy-density phase space and provides insights into changes in the index of refraction, charge carrier densities, and AC electrical conductivity at optical frequencies. At excitation energy densities exceeding 10 MJ/kg, the index of refraction falls to n = 0.7, and the thermally excited free-carrier density reaches ne = 5 × 1027 m−3, which is over an order of magnitude higher than that of the electron carriers produced by direct photoionization. Significant specular reflection is observed owing to critical electron density shielding of electromagnetic waves. The measured optical conductivity reaches 2 × 104 S/m, a value that is one to two orders of magnitude lower than those of simple metals in a liquid state. At electron temperatures below 15 000 K, the experimental results agree well with the theoretical calculations using density-functional theory/molecular-dynamics simulations. With increasing temperature, the electron density increases and the system approaches a Fermi distribution. In this regime, the conductivities agree better with predictions from the Ziman theory of liquid metals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Martensitic transformation in superlattices of two non-transforming metals

A number of studies have shown that the coherent integration of a second non-transforming phase in martensitic materials results in unprecedented thermo-mechanical properties. We demonstrate the possibility of martensitic transformation by epitaxially combining two non-transforming materials, Mg (hexagonal closed packed) and MgLi (body centered cubic). Density functional theory calculations predict that MgLi/Mg superlattices undergo a transformation crystallographically equivalent to MgSc alloys and can be expected to result in finite temperature shape memory behavior. We predict the zero-temperature energy of the martensite phase to be 9.3 meV/atom lower than austenite for a 50–50 at. % MgLi/Mg nanolaminate, which is expected to result in a martensitic transformation temperature of approximately 150 K. This energy difference can be tuned between 4.9 and 18.1 meV/atom by varying the fraction of the pure Mg phase in the superlattice, potentially enabling the transformation from cryogenic to room and higher temperatures. While Li diffusion might limit the applications of the systems proposed, this work opens the door to a new approach to create martensitic materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Structure prediction of epitaxial inorganic interfaces by lattice and surface matching with Ogre

We present a new version of the Ogre open source Python package with the capability to perform structure prediction of epitaxial inorganic interfaces by lattice and surface matching. In the lattice matching step, a scan over combinations of substrate and film Miller indices is performed to identify the domain-matched interfaces with the lowest mismatch. Subsequently, surface matching is conducted by Bayesian optimization to find the optimal interfacial distance and in-plane registry between the substrate and the film. For the objective function, a geometric score function is proposed based on the overlap and empty space between atomic spheres at the interface. The score function reproduces the results of density functional theory (DFT) at a fraction of the computational cost. The optimized interfaces are pre-ranked using a score function based on the similarity of the atomic environment at the interface to the bulk environment. Final ranking of the top candidate structures is performed with DFT. Ogre streamlines DFT calculations of interface energies and electronic properties by automating the construction of interface models. The application of Ogre is demonstrated for two interfaces of interest for quantum computing and spintronics, Al/InAs and Fe/InSb.

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

Relativistic density functional theory in the full potential linear muffin tin orbital method

In this write-up, we have outlined the general methodology for treating relativity in the FP-LMTO method. The Schrodinger equation is insufficient for describing the physics of heavy elements because it leaves out any relativistic effects. Therefore DFT based on the Schrodinger equation is not sufficient for studying heavy materials. It is, however, possible to incorporate all of the effects of relativity in DFT if starting with a reformulation of the Hohenberg-Kohn theorem that proceeds from the Dirac equation, rather than the Schrodinger equation. The corresponding theorem in relativistic density functional theory (RDFT) shows that the ground-state density of the system is a functional of not only the charge density as in Schrodinger-based DFT, but also of the current density. This leads to a more elaborate procedure in practice when solving the Dirac-Kohn-Sham equations, the relativistic analogs of the Kohn-Sham equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Pressure-induced structural and dielectric changes in liquid water at room temperature

Understanding the pressure-dependent dielectric properties of water is crucial for a wide range of scientific and practical applications. In this study, we employ a deep neural network trained on density functional theory data to investigate the dielectric properties of liquid water at room temperature across a pressure range of 0.1–1000 MPa. We observe a nonlinear increase in the static dielectric constant ɛ 0 with increasing pressure, a trend that is qualitatively consistent with experimental observations. This increase in ɛ 0 is primarily attributed to the increase in water density under compression, which enhances collective dipole fluctuations within the hydrogen-bonding network as well as the dielectric response. Despite the increase in ɛ 0 , our results reveal a decrease in the Kirkwood correlation factor G K with increasing pressure. Furthermore, this decrease in G K is attributed to pressure-induced structural distortions in the hydrogen-bonding network, which weaken dipolar correlations by disrupting the ideal tetrahedral arrangement of water molecules.

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

Real-space representation of the quasiparticle self-consistent GW self-energy and its application to defect calculations

The quasiparticle self-consistent (QS) GW (G for Green's function, W for screened Coulomb interaction) approach incorporates the corrections of the quasiparticle energies from their Kohn-Sham density functional theory (DFT) eigenvalues by means of an energy-independent and Hermitian self-energy matrix usually given in the basis set of the DFT eigenstates. By expanding these into an atom-centered basis set (specifically here the linearized muffin-tin orbitals) a real space representation of the self-energy corrections becomes possible. In this work, We show that this representation is relatively short-ranged. This offers opportunities to construct the self-energy of a complex system from parts of the system by a cut-and-paste method. Specifically for a point defect, represented in a large supercell, the self-energy can be constructed from those of the host and a smaller defect-containing cell. The self-energy of the periodic host can be constructed simply from a GW calculation for the primitive cell. We show for the case of the As Ga in GaAs that the defect part can already be well represented by a minimal eight-atom cell and allows us to construct the self-energy for a 64-atom cell in good agreement with direct QSGW calculations for the large cell. Using this approach to an even larger 216-atom cell shows the defect band approaches an isolated defect level. The calculations also allow us to identify a second defect band which appears as a resonance near the conduction band minimum. The results on the extracted defect levels agree well with Green's function calculations for an isolated defect and with experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Phase Equilibrium of Water with Hexagonal and Cubic Ice Using the SCAN Functional

Machine learning models are rapidly becoming widely used to simulate complex physicochemical phenomena with ab initio accuracy. Here, we use one such model as well as direct density functional theory (DFT) calculations to investigate the phase equilibrium of water, hexagonal ice (Ih), and cubic ice (Ic), with an eye toward studying ice nucleation. The machine learning model is based on deep neural networks and has been trained on DFT data obtained using the SCAN exchange and correlation functional. We use this model to drive enhanced sampling simulations aimed at calculating a number of complex properties that are out of reach of DFT-driven simulations and then employ an appropriate reweighting procedure to compute the corresponding properties for the SCAN functional. This approach allows us to calculate the melting temperature of both ice polymorphs, the driving force for nucleation, the heat of fusion, the densities at the melting temperature, the relative stability of ices Ih and Ic, and other properties. We find a correct qualitative prediction of all properties of interest. In some cases, quantitative agreement with experiment is better than for state-of-the-art semiempirical potentials for water. Furthermore, our results also show that SCAN correctly predicts that ice Ih is more stable than ice Ic.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗