Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ab initio calculations”

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

Semiclassical model for calculating exciton and polaron pair energetics at interfaces

Exciton and polaron pair dissociation is a key functional aspect of photovoltaic devices. To improve upon the current state of interfacial transport models, we augment the existing classical models of dielectric interfaces by incorporating results from ab initio calculations, allowing us to calculate exciton and polaron binding energies more accurately. We demonstrate the predictive capabilities of this new model using two interfaces: (i) the boron subphthalocyanine chloride (SubPc) and C60 interface, which is an archetype for many organic photovoltaic devices; and (ii) pentacene and silicon (1 0 0), which represents a hybrid between organic and inorganic semiconductors. Overall, our calculations predict that the insertion of molecular dipoles at interfaces can be used for improving polaron pair dissociation and that sharp transitions in dielectric permittivity can have a stronger effect on the polaron pair dissociation than even the electron-hole Coulomb interaction.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Short Range Correlations and Few Body Precision Measurements

This talk will discuss the exciting experimental program studying short range correlations and few body dynamics at Jefferson Lab in ?=3 nuclei. Quasi-elastic electron scattering on these light nuclei have readily calculable ground states but can already exhibit complex nuclear behavior including nucleon-nucleon (NN) short range correlations (SRCs). 3H and 3He, being mirror nuclei, exploit the maximum available isospin asymmetry and can be compared directly with ab initio calculations. This talk will summarize the recent Hall A 3H and 3He results. The first measurement of SRC nucleons in 3He using the CLAS detector from Hall B and identifying lead neutrons in ?? pairs at high missing momentum will be shown. These fascinating results will be extended in an approved experiment using the CLAS12 detector to measure high-statistics and absolute cross sections on 2H, 3He, and 3H with a wide range of kinematical coverage. This talk will discuss the recent and proposed measurements and their implications on nuclear models for few body nuclear dynamics and short range structure.

Szumila-Vance, Holly↗
DOE OSTI2020

Ab Initio Nuclear Thermodynamics

Here, we propose a new Monte Carlo method called the pinhole trace algorithm for ab initio calculations of the thermodynamics of nuclear systems. For typical simulations of interest, the computational speedup relative to conventional grand-canonical ensemble calculations can be as large as a factor of one thousand. Using a leading-order effective interaction that reproduces the properties of many atomic nuclei and neutron matter to a few percent accuracy, we determine the location of the critical point and the liquid-vapor coexistence line for symmetric nuclear matter with equal numbers of protons and neutrons. We also present the first ab initio study of the density and temperature dependence of nuclear clustering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Electronic correlations and flattened band in magnetic Weyl semimetal candidate Co 3 Sn 2 S 2

The interplay between electronic correlations and topological protection may offer a rich avenue for discovering emergent quantum phenomena in condensed matter. However, electronic correlations have so far been little investigated in Weyl semimetals (WSMs) by experiments. Here, we report a combined optical spectroscopy and theoretical calculation study on the strength and effect of electronic correlations in a magnet Co 3 Sn 2 S 2 . The electronic kinetic energy estimated from our optical data is about half of that obtained from single-particle ab initio calculations in the ferromagnetic ground state, which indicates intermediate-strength electronic correlations in this system. Furthermore, comparing the energy and side-slope ratios between the interband-transition peaks at high energies in the experimental and single-particle-calculation-derived optical conductivity spectra with the bandwidth-renormalization factors obtained by many-body calculations enables us to estimate the Coulomb-interaction strength (U ~ 4 eV) in Co 3 Sn 2 S 2 . Besides, a sharp experimental optical conductivity peak at low energy, which is absent in the single-particle-calculation-derived spectrum but is consistent with the optical conductivity peaks obtained by many-body calculations with U ~ 4 eV, indicates that an electronic band connecting the two Weyl cones is flattened by electronic correlations and emerges near the Fermi energy in Co 3 Sn 2 S 2 . Our work paves the way for exploring flat-band-generated quantum phenomena in WSMs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Formation and detection of metastable formic acid in a supersonic expansion: High resolution infrared spectroscopy of the jet-cooled cis -HCOOH conformer

High-resolution direct absorption infrared spectra of metastable cis-formic acid (HCOOH) trapped in a cis-well resonance behind a 15 kcal/mol barrier are reported for the first time, with the energetically unstable conformer produced in a supersonic slit plasma expansion of trans-formic acid/H 2 mixtures. We present a detailed high-resolution rovibrational analysis for cis-formic acid species in the OH stretch (ν 1 ) fundamental, providing first precision vibrational band origin, rotational constants, and term values, which in conjunction with ab initio calculations at the couple-cluster with single, double, and perturbative triple [CCSD(T)]/ANOn (n = 0, 1, 2) level support the experimental assignments and establish critical points on the potential energy surface for internal rotor trans-to-cis isomerization. Relative intensities for a- and b-type transitions observed in the spectra permit the transition dipole moment components to be determined in the body fixed frame and prove to be in good agreement with ab initio CCSD(T) theoretical estimates but in poor agreement with simple bond-dipole predictions. As a result, the observed signal dependence on H 2 in the discharge suggests the presence of a novel H atom radical chemical mechanism for strongly endothermic “up-hill” internal rotor isomerization between trans- and cis-formic acid conformers.

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

FeSi binary alloy electronic structure low-Si dataset (1024 atoms)

This dataset contaims the calculated atomic charge density, atomic magnetic moment, and total energy for 1600 configurations of iron-silicon (Fe-Si) binary alloys body-centered cubic (BCC) structures at 3, 6, and 9% Si. These large scale (1024 atom) ab initio calculations were produced with the LSMS code on the OLCF Summit supercomputer. LSMS GitHub repository: https://github.com/mstsuite/lsms

36 MATERIALS SCIENCE↗
DOE OSTI2022

The 1 S 0 Pairing Gap in Neutron Matter

We report ab initio calculations of the S wave pairing gap in neutron matter calculated using realistic nuclear Hamiltonians that include two- and three-body interactions. We use a trial state, properly optimized to capture the essential pairing correlations, from which we extract ground state properties by means of auxiliary field diffusion Monte Carlo simulations. We extrapolate our results to the thermodynamic limit by studying the finite-size effects in the symmetry-restored projected Bardeen-Cooper-Schrieffer (PBCS) theory and compare our results to other ab initio studies done in the past. Our quantum Monte Carlo results for the pairing gap show a modest suppression with respect to the mean-field BCS values. These results can be connected to cold atom experiments, via the unitarity regime where fermionic superfluidity assumes a unified description, and they are important in the prediction of thermal properties and the cooling of neutron stars.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Correlating the three-dimensional atomic defects and electronic properties of two-dimensional transition metal dichalcogenides

The electronic, optical and chemical properties of two-dimensional transition metal dichalcogenides strongly depend on their three-dimensional atomic structure and crystal defects. Using Re-doped MoS 2 as a model system, here we present scanning atomic electron tomography as a method to determine three-dimensional atomic positions as well as positions of crystal defects such as dopants, vacancies and ripples with a precision down to 4 pm. We further measure the three-dimensional bond distortion and local strain tensor induced by single dopants. By directly offering these experimental three-dimensional atomic coordinates to density functional theory, we obtain more accurate electronic band structures than derived from conventional density functional theory calculations that relies on relaxed three-dimensional atomic coordinates. We anticipate that scanning atomic electron tomography not only will be generally applicable to determine the three-dimensional atomic coordinates of two-dimensional materials, but also will enable ab initio calculations to better predict the physical, chemical and electronic properties of these materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Calculation of electron-ion temperature equilibration rates and friction coefficients in plasmas and liquid metals using quantum molecular dynamics

In this work, we discuss a method to calculate with quantum molecular dynamics simulations the rate of energy exchanges between electrons and ions in two-temperature plasmas, liquid metals, and hot solids. Promising results from this method were recently reported for various materials and physical conditions. Like other ab initio calculations, the approach offers a very useful comparison with the experimental measurements and permits an extension into conditions not covered by the experiments. The energy relaxation rate is related to the friction coefficients felt by individual ions due to their nonadiabatic interactions with electrons. Each coefficient satisfies a Kubo relation given by the time integral of the autocorrelation function of the interaction force between an ion and the electrons. These Kubo relations are evaluated using the output of quantum molecular dynamics calculations in which electrons are treated in the framework of finite-temperature density functional theory. The calculation presents difficulties that are unlike those encountered with the Kubo formulas for the electrical and thermal conductivities. In particular, the widely used Kubo-Greenwood approximation is inapplicable here. Indeed, the friction coefficients and the energy relaxation rate diverge in this approximation since it does not properly account for the electronic screening of electron-ion interactions. The inclusion of screening effects considerably complicates the calculations. We discuss the physically motivated approximations we applied to deal with these complications in order to investigate a widest range of materials and physical conditions. Unlike the standard method used for the electronic conductivities, the Kubo formulas are evaluated directly in the time domain and not in the energy domain, which spares one from needing to introduce an extraneous undetermined numerical parameter to account for the discrete character of the numerical density of states. We highlight interesting properties of the energy relaxation rate not shared by other electronic properties, in particular its self-averaging character. We then present a detailed parametric and convergence study with the numerical parameters, including the system size, the number of bands and k-points, and the physical approximations for the dielectric function and the exchange-correlation energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Nanoscale fluctuation of stacking fault energy strengthens multi-principal element alloys

Chemical randomness and the associated energy fluctuation are essential features of multi-principal element alloys (MPEAs). Due to these features, nanoscale stacking fault energy (SFE) fluctuation is a natural and independent contribution to strengthening MPEAs. However, existing models for conventional alloys (i.e., alloys with one principal element) cannot be applied to MPEAs. The extreme values of SFEs required by such models are unknown for MPEAs, which need to calculate the nanoscale volume relevant to the SFE fluctuation. In the present work, we developed an analytic model to evaluate the strengthening effect through the SFE fluctuation, profuse in MPEAs. The model has no adjustable parameters, and all parameters can be determined from experiments and ab initio calculations. Further, this model explains available experimental observations and provides insightful guidance for designing new MPEAs based on the SFE fluctuation. It generally applies to MPEAs in random states and with chemical short-range order.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Refining the thermochemical properties of CF, SiF, and their cations by combining photoelectron spectroscopy, quantum chemical calculations, and the Active Thermochemical Tables approach

Fluorinated species have a pivotal role in semiconductor material chemistry and some of them have been detected beyond the Earth's atmosphere. Achieving good energy accuracy on fluorinated species using quantum chemical calculations has long been a challenge. In addition, obtaining direct experimental thermochemical quantities has also proved difficult. Here, we report the threshold photoelectron and photoion yield spectra of SiF and CF radicals generated with a fluorine reactor. The spectra were analysed with the support of ab initio calculations, resulting in new experimental values for the adiabatic ionisation energies of both CF (9.128 ± 0.006 eV) and SiF (7.379 ± 0.009 eV). Using these values, the underlying thermochemical network of Active Thermochemical Tables was updated, providing further refined enthalpies of formation and dissociation energies of CF, SiF, and their cationic counterparts.

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

One-dimensionality signature in optical conductivity of heavy-fermion CeIr 3 B 2

In low dimensions, the combined effects of interactions and quantum fluctuations can lead to dramatically new physics distinct from that existing in higher dimensions. Here, we investigate the electronic and optical properties of CeIr 3 ⁢B 2 , a quasione-dimensional (1D) Kondo lattice system, using ab initio calculations. The Ce atoms in the hexagonal crystal structure form 1D chains along the c axis, with extremely short Ce-Ce distances. The quasi-1D nature of the crystal structure is well reflected in its electronic structure. Extremely flat bands emerge within the ab plane of the Brillouin zone, yielding sharp optical transitions in the corresponding optical conductivity. Further, our calculations indicate that these prominent peaks in the optical conductivity provide a clear signature of quasi-1D heavy fermion systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Nuclear structure and elastic scattering observables obtained consistently with different N N interactions

Nucleon-nucleon ($NN$) interactions based on chiral effective theories are commonly used in ab initio calculations of light nuclei. Here we present a study based on three different $NN$ interactions (up to next-to-next-to-leading order) for which structure and elastic proton scattering observables are consistently calculated for $^4$He, $^{12}$C, and $^{16}$O. The interactions are compared at the two-body level in terms of Wolfenstein amplitudes, and their predictions for ground state energies, point-proton radii, and charge form factors, as well as proton elastic scattering observables in the leading-order spectator expansion in the energy range between 65 and 160 MeV projectile energy are presented. Finally, to gain further insight into differences visible in elastic scattering observables, we investigate the behavior of the calculated effective nucleon-nucleus interactions for the $^{12}$C nucleus based on the different $NN$ interactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Long-range quadrupole electron-phonon interaction from first principles

Lattice vibrations in materials induce perturbations on the electron dynamics in the form of long-range (dipole and quadrupole) and short-range (octopole and higher) potentials. The dipole Fröhlich term can be included in current first-principles electron-phonon (e-ph) calculations and is present only in polar materials. The quadrupole e-ph interaction is present in both polar and nonpolar materials, but currently it cannot be computed from first principles. In this work, we show an approach to compute the quadrupole e-ph interaction and include it in ab initio calculations of e-ph matrix elements. The accuracy of the approach is demonstrated by comparing with direct density functional perturbation theory calculations. We apply our method to silicon as a case of a nonpolar semiconductor and tetragonal PbTiO 3 as a case of a polar piezoelectric material. In both materials we find that the quadrupole term strongly impacts the e-ph matrix elements. Analysis of e-ph interactions for different phonon modes reveals that the quadrupole term mainly affects optical modes in silicon and acoustic modes in PbTiO 3 , although the quadrupole term is needed for all modes to achieve quantitative accuracy. The effect of the quadrupole e-ph interaction on electron scattering processes and transport is shown to be important. Our approach enables accurate studies of e-ph interactions in broad classes of nonpolar, polar, and piezoelectric materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Is solid copper oxalate a spin chain or a mixture of entangled spin pairs?

Macroscopic assemblies of interacting spins give rise to a broad spectrum of behaviors determined by the spatial arrangement of the magnetic sites and the electronic interactions between them. Compounds of copper (II), in which each copper carries spin 1/2, exhibit a vast variety of physical properties. For antiferromagnetically coupled spin sites, there are two limiting scenarios: spin chains in which the spins can exhibit a long-range order or a mixture of dimers in which the spins within each pair are entangled but do not communicate with the spins from other dimers. In principle, the two types can be distinguished on the basis of experimental observations and modeling using empirically parameterized effective Hamiltonians, but in practice, ambiguity may persist for decades, as is the case for copper oxalate. Here we use high-level ab initio calculations to establish the validity of the nearest-site Heisenberg model and to predict the interaction strength between the magnetic sites. Finally, the computed magnetic susceptibility provides an unambiguous interpretation of magnetic experiments performed throughout half a century, clearly supporting the infinite spin-chain behavior of solid copper oxalate.

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

Non-Boltzmann Effects in Chain Branching and Pathway Branching for Diethyl Ether Oxidation

Low-temperature (LT) engine applications have several potential benefits, including reduced emissions and increased efficiency. Attaining these benefits requires accurate kinetic modeling of LT chain branching, which depends heavily on ketohydroperoxide (KHP) decomposition. For diethyl ether (DEE), a promising biofuel, current estimates of the KHP decomposition rate constant are largely based on empirical fits to data. In this study, we investigate the most important KHP isomer in DEE LT oxidation by applying variable reaction coordinate transition state theory to the main pathway for KHP decomposition: OO bond fission to produce •OH and a keto-alkoxy radical, •OQ'O. We also use ab initio kinetics methods to investigate the decomposition of •OQ'O, where we find dominant branching to acetic acid, with the remaining flux going to CH 3 C(O)OCHO. Additionally, new time-resolved measurements of DEE and acetic acid concentrations during LT (450–600 K) DEE oxidation are obtained in a laser photolysis flow reactor coupled with multiplexed photoionization mass spectrometry. These new experimental data, along with jet-stirred reactor data in the literature, are compared with the predictions of a recent DEE mechanism (Tran et al. Proc. Comb. Inst. 2019, 37, 511-519) that was modified with the newly calculated ab initio rate constants for KHP and •OQ'O decomposition. The predictions of the modified mechanism are quite poor when compared to the experimental data; this is primarily due to the new KHP ⇌ •OQ'O + •OH rate constant, which is 1–2 orders of magnitude slower than empirical values employed in recent mechanisms. To reconcile the new KHP rate constant and the experimental data, we explore and quantify the possible role of non-Boltzmann (nB) reaction sequences. The nB reactions have a substantial effect on both the overall mechanism reactivity and the •OQ'O branching to acetic acid. We also provide guidance on the proper implementation of nB reactions in kinetic mechanisms.

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

First-principles calculation of lattice distortions in four single phase high entropy alloys with experimental validation

Exceptional properties of high entropy alloys (HEAs) are attributed to the disordered random solid solution of multiple alloying elements. Despite numerous studies, the fundamental understanding at the electronic and atomic levels is still missing. We report a comparative study of four fcc HEAs of NiFeCoCr and NiFeCoCrX (X = Mn, Cu, or Pd) based on ab initio calculations using large supercells with 500 atoms in equal composition. After fully optimizing their structures using the VASP package, their electronic structure, interatomic bonding, partial charge distribution, and mechanical properties are calculated and compared, revealing the intricate interdependence among them. A novel parameter based on the quantum mechanical metric for internal cohesion, the total bond order density (TBOD), is used to interpret the calculated properties. The highest TBOD is found in Cantor alloy NiFeCoCrMn but lower in NiFeCoCrPd. The atomic radii vary depending on their local chemical environment. The resulting lattice distortions is validated experimentally in NiFeCoCrMn and NiFeCoCrPd. Moreover, modeling of Cu and Pd clustering in the supercell shows they have lower total energy in agreement with the observation of Cu-enhanced nano-participates in NiFeCoCrCu and Pd-induced concentration wave in NiFeCoCrPd HEAs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Development of neural network force fields for corrosion studies

To fully understand the chemistry and physics of corrosion, novel methods of simulation must be developed. One approach is designing machine learning (ML) algorithms integrated with density functional theory to develop adaptive force fields to gain insight into corrosion behavior namely at the surface of metal oxides. Current methods of modeling corrosion are slow due to the computational cost of resolving both reaction mechanics and mass transport processes. Machine learning methods can be implemented to obtain structure-activity relationships at both the molecular and bulk scale while still retaining the accuracy of density functional theory (DFT) and significantly decreasing the time needed for simulations of complex chemical processes in the various environments of corrosion. Multiscale models are needed for corrosion studies to fully understand its processes not only at the atomic length scale (chemical bonding, energies, and forces), but also at the nano and meso length scales (solid-state physics and material science processes). Current methods of study include DFT, molecular dynamics, and Monte Carlo. The limitation of DFT is that only a small number of atoms or molecules can be simulated at that level of theory. Density functional theory is used to study the electronic structure of atoms and molecules, and calculate the force component of each atom. However, these calculations are limited to about 1000 atoms. Custom periodic boundary conditions (PBC) can be used to describe the various environments and defects that affect the atomic forces to produce a large data set from which a training set can be derived. Machine learning can be utilized to overcome the barrier of modeling macroscopic and multi-scale processes from ab initio calculations through the development of adaptive force fields. Local environments determine the atomic forces of a given system, therefore adaptive force fields must be created to produce reliable quantum mechanical calculations. This can be achieved by developing a learning algorithm that uses the mapped atomic forces or fingerprint as an input to produce energies and magnetic moments as output. A systematic approach was used to begin to build a data set in order to accurately describe the atomic forces in various environments. In Figure 4 below, a simple PBC cell of Fe{sub 2}O{sub 3} was first optimized. A surface optimization was performed next, followed by a hydroxylated surface optimization. Once this calculation has converged, the adsorption of halide species to the hydroxylated surface will be investigated. TensorFlow is an open source platform for machine learning developed by Google. Using a high level application program interface (API) such as Keras allows for building and training ML models easily in a number of different environments and languages. For this project, a neural network was developed within Anaconda in Python. Future Work: Further development of reference data set; Refining neural network and learning algorithm; Fingerprinting atomic environment to enable mapping of atomic force components; Choosing appropriate training set from reference data; Learning from training set and enabling non-linear mapping of training set fingerprints and the atomic forces; Estimation of uncertainty to identify ranges of outside applicability; Testing and analysis of molecular dynamic simulations.

36 MATERIALS SCIENCE↗