Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional calculation”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Integrated Quantum-Classical Protocol for the Realistic Description of Solvated Multinuclear Mixed-Valence Transition-Metal Complexes and Their Solvatochromic Properties

Linear cyanide-bridged polymetallic complexes, which undergo photoinduced metal-to-metal charge transfer, represent prototypical systems for studying long-range electron-transfer reactions and understanding the role played by specific solute–solvent interactions in modulating the excited-state dynamics. Here, to tackle this problem, while achieving a statistically meaningful description of the solvent and of its relaxation, one needs a computational approach capable of handling large polynuclear transition-metal complexes, both in their ground and excited states, as well as the ability to follow their dynamics in several environments up to nanosecond time scales. Here, we present a mixed quantum classical approach, which combines large-scale molecular dynamics (MD) simulations based on an accurate quantum mechanically derived force field (QMD-FF) and self-consistent QMD polarized point charges, with IR and UV–vis spectral calculations to model the solvation dynamics and optical properties of a cyano-bridged trinuclear mixed-valence compound (trans-[(NC) 5 Fe III (μ-CN)Ru II (pyridine) 4 (μ-NC)Fe III (CN) 5 ] 4– ). We demonstrate the reliability of the QMD-FF/MD approach in sampling the solute conformational space and capturing the local solute–solvent interactions by comparing the results with higher-level quantum mechanics/molecular mechanics (QM/MM) MD reference data. The IR spectra calculated along the classical MD trajectories in different solvents correctly predict the red shift of the CN stretching band in the aprotic medium (acetonitrile) and the subtle differences measured in water and methanol, respectively. By explicitly including the solvent molecules around the cyanide ligands and calculating the thermal averaged absorption spectra using time-dependent density functional theory calculations within the Tamm–Dancoff approximation, the experimental solvatochromic shift is quantitatively reproduced going from water to methanol, while it is overestimated for acetonitrile. This discrepancy can likely be traced back to the lack of important dispersion interactions between the solvent cyano groups and the pyridine substituents in our micro solvation model. The proposed protocol is applied to the ground state in water, methanol, and acetonitrile and can be flexibly generalized to study excited-state nonequilibrium solvation dynamics.

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

Direct Deoxygenation of Phenol over Fe-Based Bimetallic Surfaces Using On-the-Fly Surrogate Models

We present an accelerated nudged elastic band (NEB) study of phenol direct deoxygenation (DDO) on Fe-based bimetallic surfaces using a recently developed Gaussian process regression (GPR) calculator. Our test calculations demonstrate that the GPR calculator achieves up to 3 times speedup compared to conventional density functional theory calculations while maintaining high accuracy, with energy barrier errors below 0.015 eV. Using GPR-NEB, we systematically examine the DDO mechanism on pure Fe(110) and surfaces modified with Co and Ni in both top and subsurface layers. Our results show that subsurface Co and Ni substitutions preserve favorable thermodynamics and kinetics for both C–O bond cleavage and C–H bond formation, comparable to those on the pure Fe(110) surface. In contrast, top-layer substitutions generally increase the C–O bond cleavage barrier, render the step endothermic, and result in significantly higher reverse reaction rates, making DDO unfavorable on these surfaces. This work demonstrates the effectiveness of GRR-accelerated transition state searches for complex surface reactions and provides insights into rational design of bimetallic catalysts for selective deoxygenation.

Aromatic compounds↗
DOE OSTI2024

Spectroscopic and Computational Evidence of Uranium Dihydrogen Complexes

Dihydrogen complexation, a phenomenon with robust precedent in the transition metal series, is spectroscopically detected for a uranium(III) complex and thereby extended for the first time to the 5f series. The vacant coordination site and low valence of (C 5 H 4 SiMe 3 ) 3 U prove to be key to the reversible formation of (C 5 H 4 SiMe 3 ) 3 U–H 2 (complex 1), and the paramagnetism of the f 3 center facilitates the detection of complex 1 by NMR spectroscopy. Density functional theory calculations reveal that the delocalization of the 5f electron density from (C 5 H 4 SiMe 3 ) 3 U onto the side-on dihydrogen ligand is crucial to complex formation, an unusual bonding situation for an actinide acid–base complex. Here, the spectroscopic and computational results are compared to those reported for lanthanide metallocenes to yield insight into the nature of–and future possibilities for–f-element dihydrogen complexation.

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

Tuning of the intrinsic magnetic damping parameter in epitaxial CoNi(001) films : Role of the band-filling effect

The band-filling effect, which is due to tuning the electron occupation number near the Fermi level, on the intrinsic magnetic damping parameter is demonstrated by employing high-quality epitaxial Co 1–x Ni x alloy films, in which the Fermi level and the density of states can be continuously tuned by varying the Ni concentration x. The intrinsic magnetic damping parameter, measured by the time-resolved magneto-optical Kerr effect, changes weakly at small x and increases sharply for x > 0.80. The experimental results are well reproduced by the density functional theory calculation. More interestingly, the magnetic damping parameter and the density of states near the Fermi level share similar variation trends, demonstrating their correlation. In this work, the band-filling effect in 3d magnetic transition metal alloys provides a way to tune the magnetic damping parameter as a key element in the controlg of energy loss and speed of spintronic devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Transition to collapsed tetragonal phase in CaFe 2 As 2 single crystals as seen by 57 Fe Mössbauer spectroscopy

Temperature dependent measurements of 57 Fe Mössbauer spectra on CaFe 2 As 2 single crystals in the tetragonal and collapsed tetragonal phases are reported. Clear features in the temperature dependencies of the isomer shift, relative spectra area, and quadrupole splitting are observed at the transition from the tetragonal to the collapsed tetragonal phase. From the temperature dependent isomer shift and spectral area data, an average stiffening of the phonon modes in the collapsed tetragonal phase is inferred. The quadrupole splitting increases by ~25% on cooling from room temperature to ~100 K in the tetragonal phase and is only weakly temperature dependent at low temperatures in the collapsed tetragonal phase, in agreement with the anisotropic thermal expansion in this material. In order to gain microscopic insight about these measurements, we perform ab initio density functional theory calculations of the electric field gradient and the electron density of CaFe 2 As 2 in both phases. By comparing the experimental data with the calculations we are able to fully characterize the crystal structure of the samples in the collapsed-tetragonal phase through determination of the As z coordinate. Furthermore, based on the obtained temperature dependent structural data we are able to propose charge saturation of the Fe-As bond region as the mechanism behind the stabilization of the collapsed-tetragonal phase at ambient pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Electronic structures and formation energies of pentavalent-ion-doped SnO{sub 2}: First-principles hybrid functional calculations

We studied the electronic properties and relative thermodynamic stability of several pentavalent-ion (Ta, Nb, P, Sb, and I) doped SnO{sub 2} systems using first-principles hybrid density functional theory calculations, in order to evaluate their potential as transparent conducting oxides (TCOs). I-doped SnO{sub 2}, though conductive, shows a narrowed optical band gap with respect to the undoped system due to the formation of gap states above the valence band. Nb-doped SnO{sub 2} forms localized impurity states below the conduction band bottom, suggesting that the Nb dopant exists as an Nb{sup 4+}-like cation, which is consistent with the recent experimental finding of the formation of the impurity level below the conduction band bottom [Appl. Phys. Express 5, 061201 (2012)]. Ta- and Sb-doped SnO{sub 2} display n-type conductivity, high charge carrier density, and widened optical band gap. P-doped SnO{sub 2} shows similar n-type electronic properties with that of Sb- and Ta-doped systems, and thus P-doped SnO{sub 2} is proposed as a promising candidate TCO for further experimental validation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Room-temperature ferromagnetism properties of monoclinic VO{sub 2} (M1) nanobelts

Highlights: • The unusual room temperature ferromagnetism of the VO{sub 2}(M1) nanobelts is found in magnetic measurement. • The oxygen vacancies induce the change of Bader charge of V atom from +2.34 to +1.84 (V{sup +4} to V{sup +3}). • The ferromagnetic nature is a result of the oxygen vacancies based on the calculations by DFT. - Abstract: The magnetic properties of monoclinic VO{sub 2}(M1) nanobelts synthesized by the hydrothermal method are investigated. The VO{sub 2}(M1) nanobelts exhibit ferromagnetism at room temperature, which disagrees with previous literatures. The excessive oxygen vacancies induce the ferromagnetic nature of VO{sub 2}(B) nanobelts, which is attributed to the fact that the Bader charge of V atom changes from +2.34 to +1.84 (V{sup +4} to V{sup +3}) based on the calculations by using density functional theory (DFT). The calculations indicate that oxygen vacancy may induce the presence of V{sup 3+} atom, in agreement with experimental observations. The observed ferromagnetism should be a result of oxygen vacancies in VO{sub 2}(M1) nanobelts.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Effect of spin configurations on first-principles-calculated structure and magneto-electric properties of BiFeO{sub 3}

Here we use density functional theory calculations to study the impact of various spin configurations of magnetic ions on structure and properties of BiFeO{sub 3} (BFO). Low-spin and intermediate-spin electronic configurations of Fe in tetragonal BFO impact on the c/a ratio according to the way the t{sub 2g}/e{sub g} orbitals are filled. Such metastable structures could enhance the magnetization in BFO and couple the ferroelectric and magnetic degrees of freedom in novel ways.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Oxidation state and interfacial effects on oxygen vacancies in tantalum pentoxide

First-principles density-functional theory calculations are used to study the atomistic structure, structural energetics, and electron density near the O monovacancy (V{sub O}{sup n}; n = 0,1+,2+) in both bulk, amorphous tantalum pentoxide (a-Ta{sub 2}O{sub 5}), and also at vacuum and metallic Ta interfaces. We calculate multivariate vacancy formation energies to evaluate stability as a function of oxidation state, distance from interface plane, and Fermi energy. V{sub O}{sup n} of all oxidation states preferentially segregates at both Ta and vacuum interfaces, where the metallic interface exhibits global formation energy minima. In a-Ta{sub 2}O{sub 5}, V{sub O}{sup 0} is characterized by structural contraction and electron density localization, while V{sub O}{sup 2+} promotes structural expansion and is depleted of electron density. In contrast, interfacial V{sub O}{sup 0} and V{sub O}{sup 2+} show nearly indistinguishable ionic and electronic signatures indicative of a reduced V{sub O} center. Interfacial V{sub O}{sup 2+} extracts electron density from metallic Ta, indicating that V{sub O}{sup 2+} is spontaneously reduced at the expense of the metal. This oxidation/reduction behavior suggests careful selection and processing of both oxide layer and metal electrodes for engineering memristor device operation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Theoretical oxidation state analysis of Ru-(bpy){sub 3}: Influence of water solvation and Hubbard correction in first-principles calculations

Oxidation state is a powerful concept that is widely used in chemistry and materials physics, although the concept itself is arguably ill-defined quantum mechanically. In this work, we present impartial comparison of four, well-recognized theoretical approaches based on Lowdin atomic orbital projection, Bader decomposition, maximally localized Wannier function, and occupation matrix diagonalization, for assessing how well transition metal oxidation states can be characterized. Here, we study a representative molecular complex, tris(bipyridine)ruthenium. We also consider the influence of water solvation through first-principles molecular dynamics as well as the improved electronic structure description for strongly correlated d-electrons by including Hubbard correction in density functional theory calculations.

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

Catalyst Design in Nitrate Removal

Based on the volcano plot developed by Dr. Goldsmith group (Report linked in submission), we utilized DFT (density functional theory) calculations to search for bimetallic materials in the application of catalysts in aqueous nitrate removal. The calculations are conducted via the high-throughput automated workflow package developed by our group (Github linked in submission) using VASP commercial first-principles calculation software.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2017

Projected Commutator DIIS Method for Accelerating Hybrid Functional Electronic Structure Calculations

The commutator direct inversion of the iterative subspace (commutator DIIS or C-DIIS) method developed by Pulay is an efficient and the most widely used scheme in quantum chemistry to accelerate the convergence of self-consistent field (SCF) iterations in Hartree–Fock theory and Kohn–Sham density functional theory. The C-DIIS method requires the explicit storage of the density matrix, the Fock matrix, and the commutator matrix. Hence, the method can only be used for systems with a relatively small basis set, such as the Gaussian basis set. We develop a new method that enables the C-DIIS method to be efficiently employed in electronic structure calculations with a large basis set such as planewaves for the first time. The key ingredient is the projection of both the density matrix and the commutator matrix to an auxiliary matrix called the gauge-fixing matrix. The resulting projected commutator-DIIS method (PC-DIIS) only operates on matrices of the same dimension as that consists of Kohn–Sham orbitals. The cost of the method is comparable to that of standard charge mixing schemes used in large basis set calculations. The PC-DIIS method is gauge-invariant, which guarantees that its performance is invariant with respect to any unitary transformation of the Kohn–Sham orbitals. We demonstrate that the PC-DIIS method can be viewed as an extension of an iterative eigensolver for nonlinear problems. We use the PC-DIIS method for accelerating Kohn–Sham density functional theory calculations with hybrid exchange–correlation functionals, and demonstrate its superior performance compared to the commonly used nested two-level SCF iteration procedure. Furthermore, we demonstrate that in the context of ab initio molecular dynamics (MD) simulation with hybrid functionals one can extrapolate the gauge-fixing matrix to achieve the goal of extrapolating the entire density matrix implicitly along the MD trajectory. Here, numerical results indicate that the new method significantly reduces the number of SCF iterations per MD step, compared to the commonly used strategy of extrapolating the electron density.

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

Electron-phonon relaxation and excited electron distribution in gallium nitride

We develop a theory of energy relaxation in semiconductors and insulators highly excited by the long-acting external irradiation. We derive the equation for the non-equilibrium distribution function of excited electrons. The solution for this function breaks up into the sum of two contributions. The low-energy contribution is concentrated in a narrow range near the bottom of the conduction band. It has the typical form of a Fermi distribution with an effective temperature and chemical potential. The effective temperature and chemical potential in this low-energy term are determined by the intensity of carriers' generation, the speed of electron-phonon relaxation, rates of inter-band recombination, and electron capture on the defects. In addition, there is a substantial high-energy correction. This high-energy “tail” largely covers the conduction band. The shape of the high-energy “tail” strongly depends on the rate of electron-phonon relaxation but does not depend on the rates of recombination and trapping. We apply the theory to the calculation of a non-equilibrium distribution of electrons in an irradiated GaN. Probabilities of optical excitations from the valence to conduction band and electron-phonon coupling probabilities in GaN were calculated by the density functional perturbation theory. Our calculation of both parts of distribution function in gallium nitride shows that when the speed of the electron-phonon scattering is comparable with the rate of recombination and trapping then the contribution of the non-Fermi “tail” is comparable with that of the low-energy Fermi-like component. So the high-energy contribution can essentially affect the charge transport in the irradiated and highly doped semiconductors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

First-Principles Quantum Transport Modeling of Spin-Transfer and Spin-Orbit Torques in Magnetic Multilayers

A unified approach for computing (i) spin-transfer torque in magnetic trilayers like spin valve and magnetic tunnel junction, where injected charge current flows perpendicularly to interfaces, and (ii) spin-orbit torque in magnetic bilayers of the type ferromagnet/spin-orbit-coupled material, where injected charge current flows parallel to the interface, was reviewed. The experimentally explored and technologically relevant spin-orbit-coupled materials include 5d heavy metals, topological insulators, Weyl semimetals, and transition metal dichalcogenides. This approach requires to construct the torque operator for a given Hamiltonian of the device and the steady-state nonequilibrium density matrix, where the latter is expressed in terms of the nonequilibrium Green’s functions and split into three contributions. Tracing these contributions with the torque operator automatically yields field-like and damping-like components of spin-transfer torque or spin-orbit torque vector, which is particularly advantageous for spin-orbit torque where the direction of these components depends on the unknown-in-advance orientation of the current-driven nonequilibrium spin density in the presence of spin-orbit coupling. Illustrative examples are provided by computing spin-transfer torque in a one-dimensional toy model of a magnetic tunnel junction and realistic Co/Cu/Co spin valve, both of which are described by first-principles Hamiltonians obtained from noncollinear density functional theory calculations, as well as by computing spin-orbit torque in a ferromagnetic layer described by a tight-binding Hamiltonian which includes spin-orbit proximity effect within ferromagnetic monolayers assumed to be generated by the adjacent monolayer transition metal dichalcogenide. In addition, it is shown here how spin-orbit proximity effect, quantified by computing (via first-principles retarded Green’s function) spectral functions and spin textures on monolayers of realistic ferromagnetic material like Co in contact with heavy metal or monolayer transition metal dichalcogenide, can be tailored to enhance the magnitude of spin-orbit torque. Furthermore, errors made in the calculation of spin-transfer torque are quantified when using Hamiltonian from collinear density functional theory, with rigidly rotated magnetic moments to create noncollinear magnetization configurations, instead of proper (but computationally more expensive) self-consistent Hamiltonian obtained from noncollinear density functional theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Distinguishing tautomerism in the crystal structure of (Z)-N-(5-ethyl-2,3-di-hydro-1,3,4-thiadiazol-2-ylidene) -4-methylbenzenesulfonamide using DFT-D calculations and {sup 13}C solid-state NMR

The crystal structure of (Z)-N-(5-ethyl-2,3-di-hydro-1,3,4-thiadiazol-2-ylidene) -4-methylbenzenesulfonamide contains an imine tautomer, rather than the previously reported amine tautomer. The tautomers can be distinguished using dispersion-corrected density functional theory calculations and by comparison of calculated and measured {sup 13}C solid-state NMR spectra. The crystal structure of the title compound, C{sub 11}H{sub 13}N{sub 3}O{sub 2}S{sub 2}, has been determined previously on the basis of refinement against laboratory powder X-ray diffraction (PXRD) data, supported by comparison of measured and calculated {sup 13}C solid-state NMR spectra [Hangan et al. (2010 ▶). Acta Cryst. B66, 615–621]. The mol@@ecule is tautomeric, and was reported as an amine tautomer [systematic name: N-(5-ethyl-1,3,4-thia@@diazol-2-yl)-p-toluene@@sulfonamide], rather than the correct imine tautomer. The protonation site on the mol@@ecule’s 1,3,4-thia@@diazole ring is indicated by the inter@@molecular contacts in the crystal structure: N—H⋯O hydrogen bonds are established at the correct site, while the alternative protonation site does not establish any notable inter molecular inter@@actions. The two tautomers provide essentially identical Rietveld fits to laboratory PXRD data, and therefore they cannot be directly distinguished in this way. However, the correct tautomer can be distinguished from the incorrect one by previously reported qu@@anti@@tative criteria based on the extent of structural distortion on optimization of the crystal structure using dispersion-corrected density functional theory (DFT-D) calculations. Calculation of the {sup 13}C SS-NMR spectrum based on the correct imine tautomer also provides considerably better agreement with the measured {sup 13}C SS-NMR spectrum.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Energy barriers for Dy and H penetrating graphene on 6 H -SiC(0001) and freestanding bilayer graphene from first-principles calculations

Currently, intercalation of foreign guest atoms into two-dimensional (2D) layered van der Waals materials is an active research area motivated in part by the development of next-generation energy-storage technologies and optoelectronic devices. One such extensively studied 2D material is the graphene-on-SiC system. To realize and control the desired intercalated structures, it is fundamentally important to understand the kinetic process of intercalation. For the intercalation of a guest atom into graphene layers on SiC substrate, a critical kinetic parameter is the energy barrier of a guest atom penetrating the perfect graphene top layer into the gallery under it. However, accurate theoretical calculations for such penetration barriers are unavailable in literature. From our first-principles density functional theory calculations, we obtain the global energy barriers of 3.47 and 1.80 eV for single Dy and H atoms penetrating the graphene top layer on a graphene buffer layer supported by a Si-terminated 6H-SiC(0001) substrate, respectively. Furthermore, for comparison as well as for examining the lateral strain effects, we also obtain the global barriers of 5.05 and 1.50 eV for single Dy and H atoms penetrating freestanding bilayer graphene with a tensile strain of about 8.8% to match our model for supported graphene, as well as the global barriers of 7.21 and 4.18 eV for penetrating unstrained freestanding bilayer graphene, respectively. From corresponding minimum energy paths with multiple energy minima and saddle points, we can also obtain various local energy barriers and the global backward barrier from the graphene gallery back to the top surface.

2D materials↗
DOE OSTI2014

Study of electronic transport properties of doped 8AGNR

The electronic and transport properties of 8-armchair graphene nanoribbon (8AGNR) with defect at different sites are investigated by performing first-principles calculations based on density functional theory (DFT). The calculated results show that the 8AGNR are semiconductor. The introduction of 3d transition metals, creates the nondegenerate states in the conduction band, makes 8AGNR metallic. The computed transmission spectrum confirms that AGNR are semiconducting in nature and their band gap remain unchanged and localized states appear when there is vacancy in their structures, and the conductance decreases due to defects compared with the pristine nanoribbon.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Mechanistic Investigation of Dimethylmercury Formation Mediated by a Sulfide Mineral Surface

Mercury (Hg) pollution is a global environmental problem, but many key transformations among 29 various Hg species are not well understood. Although anaerobic microorganisms are the primary 30 producers of monomethylmercury (MMeHg), the abiotic formation of dimethylmercury (DMeHg) 31 from MMeHg have been suggested to account for a large portion of DMeHg in the ocean. It has 32 been shown that abiotic formation of DMeHg from MMeHg can be facilitated by reduced sulfur 33 groups on mineral surfaces. However, a detailed mechanistic understanding of this reaction is 34 lacking. Here we perform density functional theory (DFT) calculations to explore mechanisms of 35 DMeHg formation on the 110 surface of a CdS(s) (hawleyite) nanoparticle. We show that 36 coordination of MMeHg substituents to adjacent reduced sulfur groups protruding from the surface 37 facilitate DMeHg formation and that the reaction proceeds through direct transmethylation from 38 one MMeHg substituent to another. Solvation effects were found to play an important role in the 39 surface reconstruction of the nanoparticle and in decreasing the energetic barrier for DMeHg 40 formation. We also investigate DMeHg formation from MMeHg mediated by organic dithiols and 41 compare the results to those mediated by the mineral surface. In contrast to the surface reaction, 42 the calculations suggest that dithiol-mediated formation of DMeHg instead proceeds through a 43 mechanism in which Hg–Hg interactions activate a methyl group from one MMeHg substituent 44 for intramolecular transfer to the adjacent MMeHg substituent. These findings refine previously 45 proposed mechanisms of abiotic DMeHg formation and fill gaps in our understanding of global 46 Hg cycling in the environment.

Lian, Peng↗