A density matrix approach to multiconfiguration calculations
Density matrix components for multiconfiguration wave functions, constructing N electron /spin free/ Hamiltonian configuration interaction matrix
Harriman, J. E.↗
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.
Cu Based Dilute Alloys for Tuning the C 2+ Selectivity of Electrochemical CO 2 Reduction
Electrochemical CO 2 reduction is a promising technology for replacing fossil fuel feedstocks in the chemical industry but further improvements in catalyst selectivity need to be made. So far, only copper-based catalysts have shown efficient conversion of CO 2 into the desired multi-carbon (C 2+ ) products. This work explores Cu-based dilute alloys to systematically tune the energy landscape of CO 2 electrolysis toward C 2+ products. Selection of the dilute alloy components is guided by grand canonical density functional theory simulations using the calculated binding energies of the reaction intermediates CO*, CHO*, and OCCO* dimer as descriptors for the selectivity toward C 2+ products. Here, a physical vapor deposition catalyst testing platform is employed to isolate the effect of alloy composition on the C 2+ /C 1 product branching ratio without interference from catalyst morphology or catalyst integration. Six dilute alloy catalysts are prepared and tested with respect to their C 2+ /C 1 product ratio using different electrolyzer environments including selected tests in a 100-cm 2 electrolyzer. Consistent with theory, CuAl, CuB, CuGa and especially CuSc show increased selectivity toward C 2+ products by making CO dimerization energetically more favorable on the dominant Cu facets, demonstrating the power of using the dilute alloy approach to tune the selectivity of CO 2 electrolysis.
36 MATERIALS SCIENCE↗
Prediction of BiS2-type pnictogen dichalcogenide monolayers for optoelectronics
Abstract In this work, we introduce a 2D materials family with chemical formula MX 2 (M={As, Sb, Bi} and X={S, Se, Te}) having a rectangular 2D lattice. This materials family has been predicted by systematic ab-initio structure search calculations in two dimensions. Using density-functional theory and many-body perturbation theory, we study the structural, vibrational, electronic, optical, and excitonic properties of the predicted MX 2 family. Our calculations reveal that the predicted SbX 2 and BiX 2 monolayers are stable while the AsX 2 layers exhibit an in-plane ferroelectric instability. All materials display strong excitonic effects and good optical absorption within the infrared-to-visible range. Hence, these monolayers can harvest solar energy and serve in optoelectronics applications. Furthermore, our results indicate that exfoliation of the predicted MX 2 monolayers from their bulk counterparts is experimentally viable.
Materials Science↗
Resonance Scattering of Fe XVII X-Ray and EUV Lines
Over the years a number of calculations have been carried out to derive intensities of various X-ray and EUV lines in Fe XVII to compare with observed spectra. The predicted intensities have not agreed with solar observations, particularly for the line at 15.02 Angstroms; resonance scattering has been suggested as the source for much of the disagreement. The atomic data calculated earlier used seven configurations, 2s(sup 2) 2p(sup 6),2s(sup 2) 2p(sup 5)3s, 2s(sup 2) 2p(sup 5)3p, 2s(sup 2)2p(sup 5)3d, 2s2 p(sup 6)3s, 2s2p(sup 6)3p, and 2s2p(sup 6)3d, having n=3 orbitals and the scattering calculations in the distorted wave approximation were carried out only for incident energies above the threshold of the highest fine-structure level. These calculations have now been extended to thirteen configurations by adding six more configurations having n=4 orbitals, namely 2s(sup 2) 2p(sup 5)4s, 2s(sup 2) 2p(sup 5)4p, 2s(sup 2) 2p(sup 5)4d, 2s2p(sup 6)4s, 2s2p(sup 6)4p, and 2s2p(sup 6)4d, giving rise to 73 fine structure levels. The scattering calculations are carried out below as well as above the threshold of the highest fine-structure level. The incident energies below the threshold are 55.8, 65, 70, and 76 Ry and those above are 85, 127.5, 170, 212.5, 255, 340, and 425 Ry. The collision strengths are calculated up to total angular momentum L(sup T) = 33. Level populations and intensity ratios are calculated at various electron temperatures and densities. It is not possible to predict how the ratios would change when configurations with n=5 and 6 orbitals are added but some estimates have been made for a few transitions by Liedhal, who indicates an asymptotic convergence when n=5 and 6 orbitals are added. These improved calculations of Fe XVII change the intensity ratios compared to those obtained earlier, bringing the optically thin F(15.02)/F(16.78) ratio and several other ratios closer to the observed values: However, some disagreement with the solar observations still persists, even though the agreement of the presently calculated optically thin F(15.02)/F(15.26) ratio with the experimental results of Brown et al. and Laming et al. has improved. Some of the remaining discrepancy is still thought to be the effect of opacity, which is consistent with expected physical conditions for solar sources. EUV intensity ratios are also calculated and compared with observations. Level populations and intensity ratios are calculated, as a function of column density of Fe XVII, in the slab and cylindrical geometries. As found previously, the predicted intensities for the resonance lines at 15.02 and 15.26 Angstroms exhibit initial increases in flux relative to the forbidden line at 17.10 Angstroms and the resonance line at 16.78 Angstroms as optical thickness increases. The same behavior is predicted for the lines at 12.262 and 12.122 Angstroms. Predicted intensities for some of the allowed EUV lines are also affected by opacity.
Bhatia, A. K.↗
Collection: TD-DFT and EOM-CCSD Calculations for the GDB-9-Ex Dataset
We present two datasets that contain quantum chemical electronic structure calculations for organic molecules from the GDB-9-Ex dataset. The “GDB-9-Ex_TD-DFT-PBE0” dataset contains calculations performed using the time-dependent density functional theory (TD-DFT) first principles method, and the “GDB-9-Ex_EOMCCSD” dataset contains calculations performed using the equation-of-motion coupled cluster (EOM-CCSD) method. Both types of calculations were performed using the ORCA software and provided ultraviolet-visible spectra with a high level of accuracy.
Mehta, Kshitij [Oak Ridge National Laboratory (ORN↗
Vacancy-mediated transport and segregation tendencies of solutes in fcc nickel under diffusional creep: A density functional theory study
The Nabarro-Herring (NH) diffusional creep theory postulates the vacancy-mediated transport of atoms under a stress gradient as the creep mechanism under low-stress and high-temperature conditions. In multicomponent alloys, we premise that this stress-assisted flow of vacancies to and from grain boundaries will produce elemental segregation. An observation of such segregation, validated with theoretical predictions, can provide the necessary experimental evidence for the occurrence of NH creep. Theoretical calculations of the segregation tendencies via analyzing the dominant solute diffusion mechanisms and the difference in diffusivities of the elements are therefore essential. To this end, this study applies density functional theory calculations of migration barriers and solute-vacancy binding energies as input to the self-consistent mean-field theory to assess the vacancy-mediated diffusion mechanisms, transport coefficients, and segregation tendencies of Co, Cr, Mo, Re, Ta, and W solutes in face-centered-cubic Ni. We find Co, Re, and W to be slow diffusers at high temperatures and Cr, Mo, and Ta to be fast diffusers. Further analysis shows that the slow diffusers tend to always enrich at vacancy sinks over a wide range of temperatures. In contrast, the fast diffusers show a transition from depletion to enrichment as the temperature lowers. Furthermore, our analysis of the segregation tendencies under tensile hydrostatic strains shows that slow diffusers are largely unaffected by the strain and favor enrichment. On the other hand, the fast diffusers exhibit high sensitivity to strain and their segregation tendency can transition from depletion to enrichment at a given temperature. Finally, the transport coefficients calculated in this work are expected to serve as input to mesoscale microstructure models to provide a more rigorous assessment of solute segregation under NH creep conditions.
36 MATERIALS SCIENCE↗
An improved exceedance theory for combined random stresses
An extension is presented of Rice's classic solution for the exceedances of a constant level by a single random process to its counterpart for an n-dimensional vector process. An interaction boundary, analogous to the constant level considered by Rice for the one-dimensional case, is assumed in the form of a hypersurface. The theory for the numbers of boundary exceedances is developed by using a joint statistical approach which fully accounts for all cross-correlation effects. An exact expression is derived for the n-dimensional exceedance density function, which is valid for an arbitrary interaction boundary. For application to biaxial states of combined random stress, the general theory is reduced to the two-dimensional case. An elliptical stress interaction boundary is assumed and the exact expression for the density function is presented. The equations are expressed in a format which facilitates calculating the exceedances by numerically evaluating a line integral. The behavior of the density function for the two-dimensional case is briefly discussed.
Lester, H. C.↗
Iodide Double Perovskites and the Limits of their Structural Stability
AbstractWhile halide double perovskites A2M(I)M(III)X6 have attracted significant attention, examples involving iodides are rare. We examine the limits of the structural stability of iodide double perovskites, presenting the synthesis and single‐crystal structures of Cs2NaScI6 and Cs2NaYI6. Bypassing the common expectation that iodides have small band gaps, these compounds display optical gaps of 3.10 eV [M(III) = Sc] and 3.65 eV [M(III) = Y]. Cs2NaScI6 is the only iodide double perovskite to exhibit a cubic crystal structure at room temperature. Density functional theory‐based electronic structure calculations help understand the role of competing Cs3M(III)2I9 (3 : 2 : 9) phases and provide possible reasons for why iodide double perovskites based around In(III), Sb(III) and Bi(III) cations have proved elusive. We confirm design rules for halide double perovskites based around concepts of the tolerance factor and the radius ratio of the smaller, trivalent ion, but also point to situations such as what is observed for Cs2NaScI6 where a double perovskite can be trapped in a metastable structure.
Mulligan, Anya S.↗
Crystal Growth of Quaternary Rare Earth Selenosilicates by Using the Flux-Assisted Boron Chalcogen Mixture Method: Investigation of Their Magnetic and Optical Properties
A series of rare earth magnesium selenosilicates, RE3Mg 0.5 SiSe 7 (RE = Ce, Pr, Nd, Sm, Gd, Tb, Dy) were obtained as single crystal using the flux assisted boron chalcogen mixture (BCM) method. The structures of the crystals were determined by single-crystal X-ray diffraction. The RE3Mg 0.5 SiSe 7 series crystallizes in the hexagonal crystal system in the space group P6 3 . Polycrystalline powders were synthesized to perform physical property measurements. Magnetic measurements over the 2–300 K temperature range reveal that Ce 3 Mg 0.5 SiSe 7 and Gd 3 Mg 0.5 SiSe 7 exhibit paramagnetic behavior with negative Weiss constants (θ W = −14.50, θ W = −6.13 K, respectively). The optical properties of RE 3 Mg 0.5 SiSe 7 (RE = Ce, Pr, Nd, Sm, Gd) were measured by ultraviolet–visible (UV–vis) diffuse reflectance. Density functional theory (DFT) electronic structure calculations were performed. Furthermore, a second harmonic generation measurement was performed on a polycrystalline powder of Ce 3 Mg 0.5 SiSe 7 and was found to be SHG active with an efficiency of 0.11 times the standard potassium dihydrogen phosphate (KDP).
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Evaluating Material Design Principles for Calcium-Ion Mobility in Intercalation Cathodes
Multivalent-ion batteries offer an alternative to Li-based technologies, with the potential for greater sustainability, improved safety, and higher energy density, primarily due to their rechargeable system featuring a passivating metal anode. Although a system based on the Ca 2+ /Ca couple is particularly attractive given the low electrochemical plating potential of Ca 2+ , the remaining challenge for a viable rechargeable Ca battery is to identify Ca cathodes with fast ion transport. In this work, a high-throughput computational pipeline is adapted to (1) discover novel Ca cathodes in a largely unexplored space of empty intercalation hosts and (2) develop material design rules for Ca-ion mobility. One candidate from the screening, W 2 O 3 (PO 4 ) 2 , is confirmed to have a low Nudged Elastic Band (NEB) barrier of 168 meV within a one-dimensional (1D) ion percolation topology. This candidate is subsequently synthesized and electrochemically tested, achieving reversible Ca cycling with a capacity of 25 mA h/g. To further accelerate the screening for promising Ca intercalation electrodes, machine learning (ML) Random Forest (RF) and Extreme Gradient Boosting (XGB) classification models are created with local environment descriptors based on a large, structurally and chemically diverse dataset of minimum energy pathways, spanning over 5,000 density functional theory (DFT) site energy calculations. Accuracies of 92% are achieved, material design metrics are quantified, ML force-fields are leveraged in an accelerated iteration of the screening, and a total of 27 novel Ca cathode materials are highlighted for further investigation.
25 ENERGY STORAGE↗
How Does Metal Spin State Affect Electronic Communication in Mixed-Valence Dimers? Insights from Ultrafast Near-Infrared and Soft X-ray Transient Absorption Spectroscopy
Recent advancements in photocatalysis, photovoltaics, and quantum information science take advantage of electron spin, and determining how spin multiplicity affects electron transfer is key to understanding these phenomena. Here, in this study, we examine how metal spin state affects electronic communication in an organometallic mixed-valence dimer, ferrocenyl cobaltocenium ([Fe II Cp 2 Co III Cp 2 ] + ). This complex can be photoexcited from its low-spin singlet Fe II ground state to form intermediate-spin triplet Fe II and high-spin quintet Fe II excited states. Using femtosecond optical transient absorption (OTA) spectroscopy with visible (vis), near-infrared (NIR), and short-wave IR (SWIR) probes, supported by time-dependent density functional theory (TD-DFT) calculations, we measure Fe II Co III →Fe III Co II intervalence charge transfer (IVCT) bands in each of the Fe II spin states. Mulliken–Hush analysis of the excited-state IVCT bands was used to compute the electronic coupling between the metal centers in all three spin states, which increased as quintet < triplet < singlet. Meanwhile, the peak energy of the bands, and thus the ΔG of the IVCT transition, increased as triplet < quintet < singlet. Then, to directly probe the electronic structure at each metal center, we employed picosecond soft X-ray transient absorption (XTA) spectroscopy at the Fe and Co L 3 edges. Our results show that the low-spin and high-spin states of [Fe II Cp 2 Co III Cp 2 ] + are valence-localized, while the intermediate-spin state is partially delocalized. The differences in charge delocalization are attributed to differences in orbital occupation and geometry that affect the free energy and superexchange coupling.
Burke, John H. [Univ. of Illinois at Urbana-Champa↗
Electrochemically Mediated Au–C(sp 2 ) Anchors for Molecular Electronics
Terminal anchor groups play a key role in the stability and electronic properties of molecular junctions. Single molecule junctions typically consist of two preinstalled terminal anchors linking organic molecules to metal electrodes. Here, in this work, we show that p -terphenyl derivatives containing only a single terminal anchor show conductance features similar to junctions with two preinstalled terminal anchors. A set of p -terphenyl derivatives with one terminal anchor was prepared using automated chemical synthesis and characterized using single molecule electronics experiments, molecular dynamics (MD) simulations, bulk electrochemistry and spectroscopy, and nonequilibrium Green’s function-density functional theory (NEGF-DFT) calculations. Our results show that 4-amino- p -terphenyl (PPP) and related analogs exhibit a well-defined high conductance state that is diminished or absent in other p-terphenyl derivatives lacking a preinstalled amine terminal anchor or fluorine or methyl substitutions at the terminal para position. However, a low conductance state is observed in all amino- p -terphenyl derivatives with one preinstalled anchor due to molecular junctions formed by noncovalent dimeric π–π stacking interactions. The observed high conductance state diminishes upon the addition of reducing agents and is restored upon the addition of an oxidizing agent. Our results suggest that the high conductance state arises due to Au–C(sp 2 ) bond formation facilitated by a single electron oxidation event at the electrode surface. A series of control experiments with different anchor groups shows that primary amines play a key role in forming Au–C bonds for molecular junctions. Overall, these results suggest that Au–C bond formation gives rise to high conductance pathways in organic molecules containing only one preinstalled terminal anchor. Insights from this work can be leveraged in the design of molecular electronic devices, particularly in understanding the mechanisms of molecular binding and junction formation.
charge transport↗
Characterizing Defects Inside Hexagonal Boron Nitride Using Random Telegraph Signals in van der Waals 2D Transistors
Single-crystal hexagonal boron nitride (hBN) is used extensively in many two-dimensional electronic and quantum devices, where defects significantly impact performance. Therefore, characterizing and engineering hBN defects are crucial for advancing these technologies. Here, we examine the capture and emission dynamics of defects in hBN by utilizing low-frequency noise (LFN) spectroscopy in hBN-encapsulated and graphene-contacted MoS 2 field-effect transistors (FETs). The low disorder of this heterostructure allows the detection of random telegraph signals (RTS) in large device dimensions of 100 μm 2 at cryogenic temperatures. Analysis of gate bias- and temperature-dependent LFN data indicate that RTS originates from a single trap species within hBN. By performing multi-space density functional theory (MS-DFT) calculations on a gated defective hBN/MoS 2 heterostructure model, we assign substitutional carbon atoms in boron sites as the atomistic origin of RTS. This study demonstrates the utility of LFN spectroscopy combined with MS-DFT analysis on a low-disorder all-vdW FET as a powerful means for characterizing the atomistic defects in single-crystal hBN.
77 NANOSCIENCE AND NANOTECHNOLOGY↗
Magnetic Blocking in Fluoflavine Radical-Bridged Dilanthanide Complexes
Magnetic exchange coupling is difficult to foster in polynuclear lanthanide (Ln) complexes and poorly understood. While coupling Ln ions through closed-shell ligands is inherently weak due to the contracted 4f orbitals, placing open-shell ligands instead has proven to promote orders of magnitude stronger coupling, giving rise to single-molecule magnets (SMMs) innate to real magnetic memory effect in the case of the anisotropic Ln ions. Notably, the impact of radical bridges with differing oxidation states on magnetic blocking remains unexplored due to lack of Ln SMMs with radicals in two distinct charge states. Herein, the first dilanthanide complexes (Ln = Gd, Dy) containing fluoflavine (flv) bridges, [(Cp* 2 Ln) 2 (μ-flv z )]X, (where X = [Al(OC{CF 3 } 3 ) 4 ] − (z = 1−•), 1-Ln; X = 0 (z = 2−), 2-Ln; X =[K(crypt-222)] + (z = 3−•), 3-Ln) are reported. 1-Ln and 3-Ln, comprising the flv 1−• and flv 3−• radical bridges, were investigated via single-crystal X-ray diffraction (SCXRD), ultraviolet−visible (UV−vis) spectroscopy, Superconducting Quantum Interference Device (SQUID) magnetometry, high-field electron paramagnetic resonance (HF-EPR) spectroscopy and broken-symmetry density functional theory (BS-DFT) calculations. 1-Dy and 3-Dy constitute the first SMMs innate to radicals in two differing oxidation states. 1-Dy exhibits a spin-reversal barrier U eff of 28.36 cm −1 and open magnetic hysteresis loops below 3 K. By contrast, 3-Dy displays a much higher U eff of 143(2) cm −1 and open hysteresis loops until 9.5 K, representing a record for dilanthanide SMMs containing an organic radical bridge. The boost in SMM properties in 3-Dy is attributed to spin-phonon coupling and improved frontier orbital structure. This study paves the way for advanced design strategies of polynuclear Ln SMMs.
Benner, Florian [Michigan State Univ., East Lansin↗
Performance of Heterostructural TaC/AlGaN Schottky Diodes Based on First Principles Electronic Structure Properties
Advances in ultra-wide bandgap materials, such as high Al-content AlxGa1-xN (AlGaN), are essential for next generation power electronics, but the requirement for lattice matched substrates is currently a significant obstacle. Recently, conductive TaC has emerged as a promising virtual substrate for AlGaN heteroepitaxy, with wurtzite (0001) AlxGa1-xN lattice-matched to rocksalt (111) TaC at x ~ 0.5. Thus, understanding and controlling the electronic properties of the TaC/AlGaN interface is key for developing technological applications based on TaC/AlGaN devices. Using density functional theory and electronic structure calculations, we here investigate TaC/Al0.5Ga0.5N interfaces, where we include explicit alloy models in the slab calculations. We predict the Schottky barrier height and the electric field discontinuity resulting from interface charges. Considering all possible combinations of (Ta, C) substrate termination, (Al/Ga, N) nucleation, and (Al/Ga, N) polarity, we construct a chemical potential phase diagram to identify the stable interfaces that can be accessed through variation of the synthesis conditions. The predicted interface electronic properties are implemented in device performance simulations to demonstrate a practical design for a strain-free, high-efficiency TaC/AlGaN Schottky diode with a low barrier height and without interface charges, underscoring the potential of TaC as a substrate for ultra-wide bandgap devices.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Dynamics of ballistic photocurrents driven by Coulomb scattering in a two-dimensional material
First-principles real-time time-dependent density-functional theory (rt-TDDFT) calculations reveal the existence of ballistic photocurrents generated by Coulomb scattering, which has not previously been considered as a mechanism for the bulk photovoltaic effect. With monolayer GeS as an example, it is predicted that ballistic currents can be comparable to shift currents under experimentally accessible conditions.
2-dimensional systems↗
Neutron skins probed in proton knockout from neutron-rich nuclei
Proton-induced quasifree knockout reactions provide a powerful probe of nuclear single-particle structure and reaction dynamics in both stable and neutron-rich nuclei. Here, in this work, we develop a unified theoretical framework for the calculation of inclusive (𝑝, 2𝑝) and sequential (𝑝, 3𝑝) reaction cross sections and fragment momentum distributions at intermediate and relativistic energies. The approach is based on a probabilistic extension of Glauber multiple-scattering theory combined with microscopic nuclear densities obtained from Hartree-Fock-Bogoliubov calculations using Skyrme energy-density functionals. We focus in particular on the sensitivity of total cross sections and longitudinal momentum dispersions to neutron-skin thickness along isotopic chains. Our results indicate that both (𝑝, 2𝑝) and (𝑝, 3𝑝) reactions exhibit a systematic decrease of cross section and momentum width with increasing neutron excess, reflecting enhanced attenuation and surface bias induced by neutron skins. The effect is significantly stronger for two-proton removal, suggesting that (𝑝, 3𝑝) reactions may offer enhanced sensitivity to isovector nuclear structure. These findings establish proton-induced knockout reactions as complementary hadronic probes of neutron skins and the density dependence of the nuclear symmetry energy.
direct reactions↗
Reversible phase transformations between Pb nanocrystals and a viscous liquid-like phase
Phase transformations have been a prominent topic of study for both fundamental and applied science. Solid-liquid reaction–induced phase transformations can be hard to characterize, and the transformation mechanisms are often not fully understood. Here, we report reversible phase transformations between a metal (Pb) nanocrystal and a viscous liquid-like phase unveiled by in situ liquid cell transmission electron microscopy. The reversible phase transformations are obtained by modulating the electron current density (between 1000 and 3000 electrons Å –2 s –1 ). The metal-organic viscous liquid-like phase exhibits short-range ordering with a preferred Pb-Pb distance of 0.5 nm. Assisted by density functional theory and molecular dynamics calculations, we show that the viscous liquid-like phase results from the reactions of Pb with the CH 3 O fragments from the triethylene glycol solution under electron beam irradiation. Such reversible phase transformations may find broad implementations.
36 MATERIALS SCIENCE↗