SEARCH · Engineering Papers
Results for “Density functional 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.
Exploring Anomalous Photoelectron Angular Distributions in the Photoelectron Spectra of Gd 3 O 3 – : Study of Gd 3 O 2 – and Gd 3 O 3 – Using Photoelectron Spectroscopy and Density Functional Theory Calculations
Anion photoelectron (PE) spectra of lanthanide oxide clusters obtained previously have exhibited anomalous photoelectron angular distributions which were attributed to strong PE–valence electron (PEVE) interactions. Here, to further explore this effect, we have obtained the PE spectra of Gd 3 O 2 – and Gd 3 O 3 – , two clusters that have similarly complex electronic structures but contrasting symmetries. The spectra exhibit manifolds of detachment transitions at similar binding energies in a 0.5 eV window of energy. The electron affinity of Gd 3 O 2 is measured to be 1.29 ± 0.05 eV, and that of Gd 3 O 3 is 1.31 ± 0.05 eV. As seen in previous studies on lanthanide oxide cluster anions in lower than conventional oxidation states, transitions in spectra obtained lower photon energies are more congested than those obtained with higher photon energy, a signature of strong PEVE interactions. While the detachment transitions have predominantly parallel photoelectron angular distributions (PAD), the PAD varies across the manifold of transitions in the PE spectrum of Gd 3 O 3 – in a way that suggests four different subgroups of transitions. Results of calculations on Gd 3 O 2 – suggest kite or V-shape structures with antiferromagnetic coupling between one of the 4f 7 subshells with the two others. Calculations on Gd 3 O 3 – more definitively point to ring structures with a nearly isoenergetic ferromagnetically coupled high spin (24-tet) state and a dectet state in which one of the 4f 7 subshells is antiferromagnetically coupled with the other two. Taking these results as qualitative, we propose that strong mixing between the unperturbed states predicted computationally leads to overlapping transitions with different PADs.
Assessment of Adsorbate π-Backbonding in Copper(I) Metal–Organic Frameworks via Multinuclear NMR Spectroscopy and Density Functional Theory Calculations
Not Available
Dynamic and relativistic effects on Pt–Pt indirect spin–spin coupling in aqueous solution studied by ab initio molecular dynamics and two- vs four-component density functional NMR calculations
Treating 195Pt nuclear magnetic resonance parameters in solution remains a considerable challenge from a quantum chemistry point of view, requiring a high level of theory that simultaneously takes into account the relativistic effects, the dynamic treatment of the solvent–solute system, and the dynamic electron correlation. A combination of Car-Parrinello molecular dynamics (CPMD) and relativistic calculations based on two-component zeroth order regular approximation spin–orbit Kohn–Sham (2c-ZKS) and four-component Dirac–Kohn–Sham (4c-DKS) Hamiltonians is performed to address the solvent effect (water) on the conformational changes and JPtPt1 coupling. A series of bridged PtIII dinuclear complexes [L1–Pt2(NH3)4(Am)2–L2]n+ (Am = α–pyrrolidonate and pivalamidate; L = H2O, Cl−, and Br−) are studied. The computed Pt–Pt coupling is strongly dependent on the conformational dynamics of the complexes, which, in turn, is correlated with the trans influence among axial ligands and with the angle N–C–O from the bridging ligands. The J-coupling is decomposed in terms of dynamic contributions. The decomposition reveals that the vibrational and explicit solvation contributions reduce JPtPt1 of diaquo complexes (L1 = L2 = H2O) in comparison to the static gas-phase magnitude, whereas the implicit solvation and bulk contributions correspond to an increase in JPtPt1 in dihalo (L1 = L2 = X−) and aquahalo (L1 = H2O; L2 = X−) complexes. Relativistic treatment combined with CPMD shows that the 2c-ZKS Hamiltonian performs as well as 4c-DKS for the JPtPt1 coupling.
High-pressure phase diagram of beryllium from ab initio free-energy calculations
In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.
Computational Challenges in the development of a surrogate model for Density Functional Theory calculations.
Abstract not provided.
Efficient Calculation of NMR Shielding Constants Using Composite Method Approximations and Locally Dense Basis Sets
This paper presents a systematic study of applying composite method approximations with locally dense basis sets (LDBS) to efficiently calculate NMR shielding constants in small and medium-sized molecules. The pcSseg-n series of basis sets are shown to have similar accuracy to the pcS-n series when n ≥ 1 and can slightly reduce computational costs. We identify two different LDBS partition schemes that perform very effectively for density functional calculations. Here, we select a large subset of the recent NS372 database containing 290 H, C, N, and O shielding values evaluated by reference methods on 106 molecules to carefully assess methods of the high, medium, and low computational costs to make practical recommendations. Our assessment covers conventional electronic structure methods (density functional theory and wave function) with global basis calculations, as well as their use in one of the satisfactory LDBS approaches, and a range of composite approaches, also with and without LDBS. Altogether 99 methods are evaluated. On this basis, we recommend different methods to reach three different levels of accuracy and time requirements across the four nuclei considered.
Covalency of Trivalent Actinide Ions with Different Donor Ligands: Do Density Functional and Multiconfigurational Wavefunction Calculations Corroborate the Observed “Breaks”?
A comprehensive ab initio study of periodic actinide–ligand bonding trends for trivalent actinides is performed. Relativistic density functional theory (DFT) and complete active-space (CAS) self-consistent field wavefunction calculations are used to dissect the chemical bonding in the [AnCl 6 ] 3– , [An(CN) 6 ] 3– , [An(NCS) 6 ] 3– , [An(S 2 PMe 2 ) 3 ], [An(DPA) 3 ] 3– , and [An(HOPO)] – series of actinide (An = U–Es) complexes. Except for some differences for the early actinide complexes with DPA, bond orders and excess 5f-shell populations from donation bonding show qualitatively similar trends in 5f n active-space CAS vs DFT calculations. The influence of spin–orbit coupling on donation bonding is small for the tested systems. Along the actinide series, chemically soft vs chemically harder ligands exhibit clear differences in bonding trends. There are pronounced changes in the 5f populations when moving from Pu to Am or Cm, which correlate with previously noted “breaks” in chemical trends. As a result, bonding involving 5f becomes very weak beyond Cm/Bk. We propose that Cm(III) is a borderline case among the trivalent actinides that can be meaningfully considered to be involved in ground-state 5f covalent bonding.
Iodine Vacancies do not Cause Nonradiative Recombination in Halide Perovskites
The iodine vacancy (V I ) has frequently been discussed as a strong nonradiative recombination center in halide perovskites. This proposition was mainly based on the presence of charge-state transition levels in the band gap, as found in early first-principles calculations. In this work, we perform accurate hybrid-density-functional calculations for V I in CsPbI 3 , CsSnI 3 , and CsGeI 3 and find that V I does not have any transition levels in the band gap in CsPbI 3 , in contrast to the results from calculations based on semilocal functionals. The iodine vacancy V I does introduce levels in the band gap in CsSnI 3 and CsGeI 3 , but our explicitly computed nonradiative capture coefficients demonstrate that V I has a negligible impact on nonradiative recombination. Our study corrects a misunderstanding of the role of V I in the iodide-based perovskites, and shifts the focus toward identifying and mitigating actual recombination centers in order to further improve the optoelectronic performance.
Some Problems in Density Functional Theory
Though calculations based on density functional theory (DFT) are used remarkably widely in chemistry, physics, materials science, and biomolecular research and though the modern form of DFT has been studied for almost 60 years, some mathematical problems remain. From a physical science perspective, it is far from clear whether those problems are of major import. For context, we provide an outline of the basic structure of DFT as it is presented and used conventionally in physical sciences, note some unresolved mathematical difficulties with those conventional demonstrations, then pose several questions regarding both the time-independent and time-dependent forms of DFT that could benefit from attention in applied mathematics. Furthermore, progress on any of these would aid in development of better approximate functionals and in interpretation of DFT.
Giant Magnetic and Optical Anisotropy in Cerium-Substituted M-Type Strontium Hexaferrite Driven by $4f$ Electrons
By performing density-functional calculations, we find a giant magnetocrystalline anisotropy (MCA) constant in abundant element cerium $\mathrm{Ce}$ substituted M-type hexaferrite, in the energetically favorable strontium site, assisted by a quantum confined electron transfer from $\mathrm{Ce}$ to a specific iron $\mathrm{2a}$ site. Remarkably, the calculated electronic structure shows that the electron transfer leads to the formation of $\mathrm{Ce}^{3+}$ and $\mathrm{Fe}^{2+}$ at the $\mathrm{2a}$ site producing an occupied $\mathrm{Ce}$$(4f^1)$ state below the Fermi level that adds a significant contribution to MCA and magnetic moment. A half $\mathrm{Ce}$ substitution forms a metallic state, while a full substitution retains the semiconducting state of the strontium hexaferrite (host). In the latter, the band gap is reduced due to the formation of charge-transferred states in the gap region of the host. The optical absorption coefficient shows an enhanced anisotropy between light polarization in parallel and perpendicular directions. Calculated formation energies, including the analysis of probable competing phases, and elastic constants confirm that both compositions are chemically and mechanically stable. With successful synthesis, the $\mathrm{Ce}$ hexaferrite can be an alternative high-performing critical-element-free permanent magnet material adapted for use in devices such as automotive traction drive motors.
Role of carbon and hydrogen in limiting n -type doping of monoclinic ( Al x Ga 1 – x ) 2 O 3
In this study, we use hybrid density functional calculations to assess n -type doping in monoclinic ( Al x Ga 1 – x ) 2 O 3 alloys. We focus on silicon, the most promising donor dopant, and study the structural properties, formation energies, and charge-state transition levels of its various configurations. We also explore the impact of carbon and hydrogen, which are common impurities in metal-organic chemical vapor deposition (MOCVD). In Ga 2 O 3 , Si Ga is an effective shallow donor, but in Al 2 O 3 Si Al acts as a DX center with a ( + / – ) transition level in the band gap. Interstitial hydrogen acts as a shallow donor in Ga 2 O 3 but behaves as a compensating acceptor in n -type Al 2 O 3 . Interpolation indicates that Si is an effective donor in ( Al x Ga 1 – x ) 2 O 3 up to 70% Al, but it can be compensated by hydrogen already at 1% Al. We also assess the diffusivity of hydrogen and study complex formation. Si cation – H complexes have relatively low binding energies. Substitutional carbon on a cation site acts as a shallow donor in Ga 2 O 3 , but can be stable in a negative charge state in ( Al x Ga 1 – x ) 2 O 3 when x > 5 % . Substitutional carbon on an oxygen site ( C O ) always acts as an acceptor in n -type ( Al x Ga 1 – x ) 2 O 3 , but will incorporate only under relatively oxygen-poor conditions. C O – H complexes can actually incorporate more easily, explaining observations of carbon-related compensation in Ga 2 O 3 grown by MOCVD. We also investigate C cation – H complexes, finding they have high binding energies and act as compensating acceptors when x > 56 % ; otherwise the hydrogen just passivates the unintentional carbon donors. C-H complex formation explains why MOCVD-grown Ga 2 O 3 can exhibit record-low free-carrier concentrations, in spite of the unavoidable incorporation of carbon. Our study highlights that, while Si is in principle a suitable shallow donor in ( Al x Ga 1 – x ) 2 O 3 alloys up to high Al compositions, control of unintentional impurities is essential to avoid compensation.
Spectral-partitioned Kohn-Sham density functional theory
Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.
Interaction energy and isosteric heat of adsorption between hydrogen and magnesium diboride
Hydrogen storage materials form a crucial research topic for future energy utilization employing hydrogen and among those of interest magnesium diboride (MgB 2 ) has shown its prevalence. In this study, a first-principles analytical adsorption model of one hydrogen molecule in the vicinity of various magnesium diboride crystal surfaces was developed in order to obtain surface thermodynamic properties as a function of molecular and lattice properties. Henry's law constant (K H ) and isosteric heat of adsorption (ΔH ads ) indicators of the affinity between a gaseous molecule and a solid surface are thus calculated. The results in this paper not only address questions pertaining to the first stage of hydrogen storage processes but also advance the understanding of physisorption thermodynamics of a neutral molecule (H 2 ) coming in contact with a layered metallic-like surface (MgB 2 ). Although the model is built from a framework of classical calculations, quantum effects are incorporated as the fractional charge of the ions on the free surfaces, which is essential for the calculation of analytic thermodynamic values that approximate calculations from other methods. To benchmark our theoretical models, periodic density functional calculations were performed to determine the interactions between H 2 and different MgB 2 surfaces from first-principles. By considering both the top and sublayers of MgB 2 in calculating interaction energy, we have analytically and computationally calculated the interaction energies of H 2 molecules and MgB 2 's terminated planes, and witnessed the strong dependence of interaction energies on surface charges. We have also observed a dipole flipping phenomenon which explains the discontinuity seen in the interaction energy graph of Mg(0001). Both analytical and computational results showed heat of adsorption at zero coverage varying at a very low range (<7 kJ mol -1 ).
GDB-9-Ex_TD-DFT-PBE0: Dataset containing Time Dependent Density Functional Theory (TDDFT) calculations for organic molecules of the GDB-9-Ex dataset.
This dataset contains data-intensive quantum chemical electronic structure calculations for 96,766 organic molecules of the GDB-9-Ex dataset. Calculations were performed using the Time Dependent Density Functional Theory (TDDFT) first principles method using the ORCA software. It provides UV-vis spectra calculations of molecules with a high level of accuracy. The optical spectra behavior was collected based on the optimized molecular geometries in the DFTB method with 3ob parameters. All calculations utilized the def2-TZVP basis sets with the auxiliary def2/J and def2-TZVP/C basis sets. The time-dependent density-functional theory (TDDFT) approach with the PBE0 exchange-correlation functional and ORCAs default integration grid was employed. For the excitation energy calculations, the lowest 50 excitation states were calculated.
One-bond 13 C– 13 C spin-coupling constants in saccharides: a comparison of experimental and calculated values by density functional theory using solid-state 13 C NMR and X-ray crystallography
Methyl aldohexopyranosides were 13 C-labeled at contiguous carbons, crystallized, and studied by single-crystal X-ray crystallography and solid-state 13 C nuclear magnetic resonance (NMR) spectroscopy to examine the degree to which density functional theory (DFT) can calculate one-bond 13C–13C spin-coupling constants ( 1 J CC ) in saccharides with sufficient accuracy to permit their use in MA'AT analysis, a newly-reported hybrid DFT/NMR method that provides probability distributions of molecular torsion angles in solution (Zhang et al., J. Phys. Chem. B, 2017, 121, 3042–3058; Meredith et al., J. Chem. Inf. Model., 2022, 62, 3135–3141). Experimental 1 J CC values in crystalline samples of the doubly 13 C-labeled compounds were measured by solid-state 13 C NMR and compared to those calculated from five different DFT models: (1) 1 J CC values calculated from single structures identical to those observed in crystalline samples by X-ray crystallography (all atom refinement); (2) 1 J CC values calculated from the single structures in (1) but after Hirshfeld atom refinement (HAR); (3) 1 J CC values calculated from the single structures in (1) after DFT-optimization of hydrogen atoms only; and (4 and 5) 1 J CC values calculated in rotamers of torsion angle θ 2 (C1–C2–O2–O2H) or ω(C4–C5–C6–O6) from which either specific or generalized parameterized equations were obtained and used to calculate 1 J CC values in the specific θ 2 or ω rotamers observed in crystalline samples. Good qualitative agreement was observed between calculated 1 J CC values and those measured by solid-state 13 C NMR regardless of the DFT model, but in no cases were calculated 1 J CC values quantitative, differing (over-estimated) on average by 4–5% from experimental values. These findings, and those reported recently from solution NMR studies (Tetrault et al., J. Phys. Chem. B 2022, 126, 9506–9515), indicate that improvements in DFT calculations are needed before calculated 1 J CC values can be used directly as reliable constraints in MA'AT analyses of saccharides in solution.
Reproducibility of calculations on Li species with correlation-consistent basis sets
Here we have found that two different basis set series, both labeled cc-pVxZ, are widely available and widely used for calculation on Li compounds such that calculations that should be the same are different, causing a reproducibility problem. Here we study the magnitude of the disagreement of these nominally-the-same basis sets for wave function and density functional calculations of Li clusters. We find deviations among the basis sets of as large as 22 millihartrees. In most cases one obtains lower energies with the noncanonical choices. Therefore, it is important to specify which version of the basis set is used for calculations.