Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional approximations”

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 847 records · Page 47

Carbon X-ray absorption spectra of fluoroethenes and acetone: A study at the coupled cluster, density functional, and static-exchange levels of theory

Near carbon K-edge X-ray absorption fine structure spectra of a series of fluorine-substituted ethenes and acetone have been studied using coupled cluster and density functional theory (DFT) polarization propagator methods, as well as the static-exchange (STEX) approach. With the complex polarization propagator (CPP) implemented in coupled cluster theory, relaxation effects following the excitation of core electrons are accounted for in terms of electron correlation, enabling a systematic convergence of these effects with respect to electron excitations in the cluster operator. Coupled cluster results have been used as benchmarks for the assessment of propagator methods in DFT as well as the state-specific static-exchange approach. Calculations on ethene and 1,1-difluoroethene illustrate the possibility of using nonrelativistic coupled cluster singles and doubles (CCSD) with additional effects of electron correlation and relativity added as scalar shifts in energetics. It has been demonstrated that CPP spectra obtained with coupled cluster singles and approximate doubles (CC2), CCSD, and DFT (with a Coulomb attenuated exchange-correlation functional) yield excellent predictions of chemical shifts for vinylfluoride, 1,1-difluoroethene, trifluoroethene, as well as good spectral features for acetone in the case of CCSD and DFT. Following this, CPP-DFT is considered to be a viable option for the calculation of X-ray absorption spectra of larger {pi}-conjugated systems, and CC2 is deemed applicable for chemical shifts but not for studies of fine structure features. The CCSD method as well as the more approximate CC2 method are shown to yield spectral features relating to {pi}*-resonances in good agreement with experiment, not only for the aforementioned molecules but also for ethene, cis-1,2-difluoroethene, and tetrafluoroethene. The STEX approach is shown to underestimate {pi}*-peak separations due to spectral compressions, a characteristic which is inherent to this method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An ab initio study of multiple phonon scattering resonances in silicon germanium alloys

We have computed here phonon scattering rates and density of states in silicon germanium alloys using Green's function calculations and density functional theory. This method contrasts with the virtual crystal approximation (VCA) used in conjunction with Fermi's golden rule, which cannot capture resonance states occurring through the interaction of substitutional impurities with the host lattice. These resonances are demonstrated by density of states and scattering rate calculations in the dilute limit and show broadening as the concentration increases. Although these deviations become significant from the VCA at high frequencies, the relaxation times obtained for these phonon modes are small in both the full scattering theory and the VCA, resulting in their negligible contribution to thermal transport.

36 MATERIALS SCIENCE↗

Electronic Structure and Quadrupole Interactions in Promising Cathode Materials Na{sub x}M{sub y}(MoO{sub 4}){sub 3}, M = Mn, Fe, Co, and Ni

The electronic structure and the magnetic properties of molybdates Na{sub x}M{sub y}(MoO{sub 4}){sub 3} (M = Mn, Fe, Co, and Ni) which are promising materials for sodium batteries have been studied in the framework of the density functional theory with the GGA and GGA+U approximations for the first time. The calculations show that all the compounds are insulators. An important role of the correlation effects, provided by the on-site Coulomb interactions, was established in the formation of the band gap in these compounds. The quadrupole constants of {sup 23}Na nuclei are calculated in the nonmagnetic and ferromagnetic states within the GGA and GGA+U approaches. It is shown that the quadrupole frequencies for nonequivalent crystallographic positions of sodium are in different frequency ranges, which allows to study the diffusion of sodium in these compounds by the Nuclear Magnetic Resonance method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electronic structure and crystal phase stability of palladium hydrides

The results of electronic structure calculations for a variety of palladium hydrides are presented. The calculations are based on density functional theory and used different local and semilocal approximations. The thermodynamic stability of all structures as well as the electronic and chemical bonding properties are addressed. For the monohydride, taking into account the zero-point energy is important to identify the octahedral Pd-H arrangement with its larger voids and, hence, softer hydrogen vibrational modes as favorable over the tetrahedral arrangement as found in the zincblende and wurtzite structures. Stabilization of the rocksalt structure is due to strong bonding of the 4d and 1s orbitals, which form a characteristic split-off band separated from the main d-band group. Increased filling of the formerly pure d states of the metal causes strong reduction of the density of states at the Fermi energy, which undermines possible long-range ferromagnetic order otherwise favored by strong magnetovolume effects. For the dihydride, octahedral Pd-H arrangement as realized, e.g., in the pyrite structure turns out to be unstable against tetrahedral arrangement as found in the fluorite structure. Yet, from both heat of formation and chemical bonding considerations, the dihydride turns out to be less favorable than the monohydride. Finally, the vacancy ordered defect phase Pd{sub 3}H{sub 4} follows the general trend of favoring the octahedral arrangement of the rocksalt structure for Pd:H ratios less or equal to one.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Supervised learning of a chemistry functional with damped dispersion

Abstract Kohn–Sham density functional theory is widely used in chemistry, but no functional can accurately predict the whole range of chemical properties, although recent progress by some doubly hybrid functionals comes close. Here, we optimized a singly hybrid functional called CF22D with higher across-the-board accuracy for chemistry than most of the existing non-doubly hybrid functionals by using a flexible functional form that combines a global hybrid meta-nonseparable gradient approximation that depends on density and occupied orbitals with a damped dispersion term that depends on geometry. We optimized this energy functional by using a large database and performance-triggered iterative supervised training. We combined several databases to create a very large, combined database whose use demonstrated the good performance of CF22D on barrier heights, isomerization energies, thermochemistry, noncovalent interactions, radical and nonradical chemistry, small and large systems, simple and complex systems and transition-metal chemistry.

Liu, Yiwei (ORCID:0000000288126163)↗

Voltage-Dependent First-Principles Simulation of Insertion of Chloride Ions into Al/Al 2 O 3 Interfaces Using the Quantum Continuum Approximation

Experiments have shown that pitting corrosion can develop in aluminum surfaces at potentials > − 0.5 V relative to the standard hydrogen electrode (SHE). Until recently, the onset of pitting corrosion in aluminum has not been rigorously explored at an atomistic scale because of the difficulty of incorporating a voltage into density functional theory (DFT) calculations. We introduce the Quantum Continuum Approximation (QCA) which self-consistently couples explicit DFT calculations of the metal-insulator and insulator-solution interfaces to continuum Poisson-Boltzmann electrostatic distributions describing the bulk of the insulating region. By decreasing the number of atoms necessary to explicitly simulate with DFT by an order of magnitude, QCA makes the first-principles prediction of the voltage of realistic electrochemical interfaces feasible. After developing this technique, we apply QCA to predict the formation energy of chloride atoms inserting into oxygen vacancies in Al(111)/α-Al 2 O 3 (0001) interfaces as a function of applied voltage. We predict that chloride insertion is only favorable in systems with a grain boundary in the Al 2 O 3 for voltages > − 0.2 V (SHE). Here our results roughly agree with the experimentally demonstrated onset of corrosion, demonstrating QCA's utility in modeling realistic electrochemical systems at reasonable computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The EMC Effect in A=3 Nuclei

The European Muon Collaboration (EMC) discovered an unexpected and puzzling result in 1983 when comparing the deep inelastic scattering(DIS) nuclear structure functions of Deuterium and Iron. The per-nucleon structure functions were found to be different for the two nuclei, rather than a simple average over the proton and neutron structure functions. In subsequent experiments, this phenomenon was confirmed for additional nuclei, with the magnitude approximately scaling with the density. The exact mechanism leading to this in- medium structure function modification has not been identified. I will discuss recent results from JLab exploring the EMC effect by using an electron beam to probe two mirror nuclei, Helium-3 and Tritium. iv

Bane, Jason↗

Distribution of ion current density on a rotating spherical cap substrate during ion-assisted deposition

The uniformity of ion density is critical for applications relying on the ion assisted deposition technique for the fabrication of the high quality thin films. The authors propose and describe here a method allowing one to calculate the ion density distribution on spherical substrate holders under stationary and rotating conditions for different positions of the ion source. The ion beam shape was approximated by a cos{sup n} function, and the ion current density was represented by a function inversely proportional to the distance from the ion source in accordance with our experimental results. As an example, a calculation of the current density distribution on the spherical cap substrate was performed for a broad beam ion source operated with an anode current of 3 A. The authors propose an approach for process optimization with respect to the ion source position and its inclination, in terms of uniformity and absolute value of the ion current density.

42 ENGINEERING↗

Global description of $β$ – decay with the axially deformed Skyrme finite-amplitude method: Extension to odd-mass and odd-odd nuclei

Here we use the finite-amplitude method (FAM), an efficient implementation of the quasiparticle random phase approximation, to compute $β$-decay rates with Skyrme energy-density functionals for 3983 nuclei, essentially all the medium-mass and heavy isotopes on the neutron-rich side of stability. We employ an extension of the FAM that treats odd-mass and odd-odd nuclear ground states in the equal filling approximation. Our rates are in reasonable agreement both with experimental data where available and with rates from other global calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

mBJLDA approach for optical properties of Nowontny-Juza LiMgZ (Z = P, As) compounds

First principles approach has been utilized to investigate electronic and optical properties of LiMgZ (Z = P, As) Newontny-Juza compounds. The exchange and correlation (XC) effects are taken into account by a semi local, orbital independent modified Becke-Johnson (mBJ) potential as coupled with Local Density Approximation (LDA). The local orbital independent mBJ potential could provide better band gap/ HM gap due to its capability to catch the essentials of hybrid functionals. The positive values of real dielectric-function rule out the possibility of showing magnetism by these compounds and stabilize the semiconducting state. The shifting of peaks in optical spectra towards high energy is due to increasing indirect band gap from LiMgP to LiMgAs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Hartree potential dependent exchange functional

We introduce a novel non-local ingredient for the construction of exchange density functionals: the reduced Hartree parameter, which is invariant under the uniform scaling of the density and represents the exact exchange enhancement factor for one- and two-electron systems. The reduced Hartree parameter is used together with the conventional meta-generalized gradient approximation (meta-GGA) semilocal ingredients (i.e., the electron density, its gradient, and the kinetic energy density) to construct a new generation exchange functional, termed u-meta-GGA. This u-meta-GGA functional is exact for the exchange of any one- and two-electron systems, is size-consistent and non-empirical, satisfies the uniform density scaling relation, and recovers the modified gradient expansion derived from the semiclassical atom theory. For atoms, ions, jellium spheres, and molecules, it shows a good accuracy, being often better than meta-GGA exchange functionals. Our construction validates the use of the reduced Hartree ingredient in exchange-correlation functional development, opening the way to an additional rung in the Jacob’s ladder classification of non-empirical density functionals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab-initio calculations of electronic, transport, and structural properties of boron phosphide

We present results from ab-initio, self-consistent density functional theory calculations of electronic and related properties of zinc blende boron phosphide (zb-BP). We employed a local density approximation potential and implemented the linear combination of atomic orbitals formalism. This technique follows the Bagayoko, Zhao, and Williams method, as enhanced by the work of Ekuma and Franklin. The results include electronic energy bands, densities of states, and effective masses. The calculated band gap of 2.02 eV, for the room temperature lattice constant of a=4.5383 Å, is in excellent agreement with the experimental value of 2.02±0.05 eV. Our result for the bulk modulus, 155.7 GPa, agrees with experiment (152–155 GPa). Our predictions for the equilibrium lattice constant and the corresponding band gap, for very low temperatures, are 4.5269 Å and 2.01 eV, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Activation of N 2 on Manganese Nitride-Supported Ni 3 and Fe 3 Clusters and Relevance to Ammonia Formation

Dual-site models were constructed to represent manganese nitride (Mn 4 N)-supported Ni 3 and Fe 3 clusters for NH 3 synthesis. Density functional theory calculations produced an energy barrier of approximately 0.55 eV for N–N bond activation at the interfacial nitrogen vacancy sites (N v ); also, the hydrogenation and removal of interfacial N is promoted by earth-abundant Ni and Fe metals. Steady-state microkinetic modeling revealed that the turnover frequencies of NH 3 production follow an order of Fe 3 @Mn 4 N ≈ Ni 3 @Mn 4 N > Mn 4 N > Fe >> Ni. Moreover, we present clear evidence that, before NH 3 formation, NH migrates from N v onto the metallic sites. Using N binding energy (BE N ) and the transition-state energy of N 2 activation (E TS ) as descriptors, we concluded that the beneficial effects owing to interfacial N v sites are the most pronounced when BE N is either too strong or too weak while E TS is high; otherwise, excessive N v sites may hinder catalyst performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A first principles approach to the electronic properties of liquid and supercritical CO{sub 2}

The electronic absorption spectra of liquid and supercritical CO{sub 2} (scCO{sub 2}) are investigated by coupling a many-body energy decomposition scheme to configurations generated by Born-Oppenheimer molecular dynamics. A Frenkel exciton Hamiltonian formalism was adopted and the excitation energies were calculated with time dependent density functional theory. A red-shift of ∼ 0.2 eV relative to the gas-phase monomer is observed for the first electronic absorption maximum in liquid and scCO{sub 2}. The origin of this shift, which is not very dependent on deviations from the linearity of the CO{sub 2} molecule, is mainly related to polarization effects. However, the geometry changes of the CO{sub 2} monomer induced by thermal effects and intermolecular interactions in condensed phase lead to the appearance of an average monomeric electric dipole moment 〈μ〉 = 0.26 ± 0.04 D that is practically the same at liquid and supercritical conditions. The predicted average quadrupole moment for both liquid and scCO{sub 2} is 〈Θ〉 = − 5.5 D Å, which is increased by ∼ −0.9 D Å relative to its gas-phase value. The importance of investigating the electronic properties for a better understanding of the role played by CO{sub 2} in supercritical solvation is stressed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Non-perturbative calculation of molecular magnetic properties within current-density functional theory

We present a novel implementation of Kohn–Sham density-functional theory utilizing London atomic orbitals as basis functions. External magnetic fields are treated non-perturbatively, which enable the study of both magnetic response properties and the effects of strong fields, using either standard density functionals or current-density functionals—the implementation is the first fully self-consistent implementation of the latter for molecules. Pilot applications are presented for the finite-field calculation of molecular magnetizabilities, hypermagnetizabilities, and nuclear magnetic resonance shielding constants, focusing on the impact of current-density functionals on the accuracy of the results. Existing current-density functionals based on the gauge-invariant vorticity are tested and found to be sensitive to numerical details of their implementation. Furthermore, when appropriately regularized, the resulting magnetic properties show no improvement over standard density-functional results. An advantage of the present implementation is the ability to apply density-functional theory to molecules in very strong magnetic fields, where the perturbative approach breaks down. Comparison with high accuracy full-configuration-interaction results show that the inadequacies of current-density approximations are exacerbated with increasing magnetic field strength. Standard density-functionals remain well behaved but fail to deliver high accuracy. The need for improved current-dependent density-functionals, and how they may be tested using the presented implementation, is discussed in light of our findings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab initio study of the composite phase diagram of Ni–Mn–Ga shape memory alloys

The magnetic and structural properties of a series of nonstoichiometric Ni–Mn–Ga Heusler alloys are theoretically investigated in terms of the density functional theory. Nonstoichiometry is formed in the coherent potential approximation. Concentration dependences of the equilibrium lattice parameter, the bulk modulus, and the total magnetic moment are obtained and projected onto the ternary phase diagram of the alloys. The stable crystalline structures and the magnetic configurations of the austenitic phase are determined.

36 MATERIALS SCIENCE↗

The problem of the universal density functional and the density matrix functional theory

The analysis in this paper shows that the Hohenberg-Kohn theorem is the constellation of two statements: (i) the mathematically rigorous Hohenberg-Kohn lemma, which demonstrates that the same ground-state density cannot correspond to two different potentials of an external field, and (ii) the hypothesis of the existence of the universal density functional. Based on the obtained explicit expression for the nonrel-ativistic particle energy in a local external field, we prove that the energy of the system of more than two non-interacting electrons cannot be a functional of the inhomogeneous density. This result is generalized to the system of interacting electrons. It means that the Hohenberg-Kohn lemma cannot provide justification of the universal density functional for fermions. At the same time, statements of the density functional theory remain valid when considering any number of noninteracting ground-state bosons due to the Bose condensation effect. In the framework of the density matrix functional theory, the hypothesis of the existence of the universal density matrix functional corresponds to the cases of noninteracting particles and to interaction in the Hartree-Fock approximation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nonequilibrium solvent effects in Born-Oppenheimer molecular dynamics for ground and excited electronic states

The effects of solvent on molecular processes such as excited state relaxation and photochemical reaction often occurs in a nonequilibrium regime. Dynamic processes such as these can be simulated using excited state molecular dynamics. In this work, we describe methods of simulating nonequilibrium solvent effects in excited state molecular dynamics using linear-response time-dependent density functional theory and apparent surface charge methods. These developments include a propagation method for solvent degrees of freedom and analytical energy gradients for the calculation of forces. Molecular dynamics of acetaldehyde in water or acetonitrile are demonstrated where the solute-solvent system is out of equilibrium due to photoexcitation and emission.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗