Search NASA⌕ Search

SEARCH · Search NASA

Results for “Electronic Structure Calculations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Direct Evidence of Visible Light-Induced Homolysis in Chlorobis(2,9-dimethyl-1,10-phenanthroline)copper(II)

Developments in the field of photoredox catalysis that leveraged the long-lived excited states of Ir(III) and Ru(II) photosensitizers to enable radical coupling processes paved the way for explorations of synthetic transformations that would otherwise remain unrealized. While first row transition metal photocatalysts have not been as extensively investigated, valuable synthetic transformations covering broad scopes of olefin functionalization have been recently reported featuring photoactivated chlorobis(phenanthroline) Cu(II) complexes. In this study, the photochemical processes underpinning the catalytic activity of [Cu(dmp) 2 Cl]Cl (dmp = 2,9-dimethyl-1,10-phenanthroline) were studied. The combined results from static spectroscopic measurements and conventional photochemistry, ultrafast transient absorption, and electron paramagnetic resonance spin trapping experiments strongly support blue light (λ ex = 427 or 470 nm)-induced Cu–Cl homolytic bond cleavage in [Cu(dmp) 2 Cl] + occurring in <100 fs. On the basis of electronic structure calculations, this bond-breaking photochemistry corresponds to the Cl → Cu(II) ligand-to-metal charge transfer transition, unmasking a Cu(I) species [Cu(dmp) 2 ] + and a Cl atom, thereby serving as a departure point for both Cu(I)- or Cu(II)-based photoredox transformations. No net photochemistry was observed through direct excitation of the ligand-field transitions in the red (λ ex = 785 or 800 nm), and all combined experiments indicated no evidence of Cu–Cl bond cleavage under these conditions. As a result, the underlying visible light-induced homolysis of a metal–ligand bond yielding a one-electron-reduced photosensitizer and a radical species may form the basis for novel transformations initiated by photoinduced homolysis featuring in situ-formed metal–substrate adducts utilizing first row transition metal complexes.

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

Eu 5 Al 3 Sb 6 : Al 4 Tetrahedra Embedded in a Rock-Salt-Like Structure

The new Eu 5 Al 3 Sb 6 phase has been successfully synthesized as a pure phase through Sn flux methods yielding large, high-quality crystals. This structure type features disordered Al clusters that appear in the form of dual tetrahedra. It crystallizes in the monoclinic C2/m space group exhibiting a rock-salt-like Eu–Sb framework with [Al 4 ] tetrahedra replacing some of the cationic Eu atoms (space group: C2/m, a = 8.151(1) Å, b = 14.181(2) Å, c = 8.145(1) Å, β = 109.577(2)°). The structure models the [Al 4 ] as dual tetrahedra with the Al atom sites 37.5% occupied along with Eu present on the central site at 8% occupancy and the remainder of the site being vacant. The presence of the [Al 4 ] cluster is further supported by HRTEM. Electronic structure calculations show that this material is a semimetal with observed band crossings close to the Fermi level. Strong Al–Sb antibonding interactions were found from COHP calculations close to the Fermi level and provide the rationale for the deficiency of the Al cluster. Mössbauer spectroscopy on Eu-151 and Sb-121 provides oxidation states of 2+ and 3– along with the local environment. Magnetic susceptibility measurements can be described well with a Curie–Weiss law where an effective moment of 7.80 μB/mol Eu is obtained, consistent with Eu 2+ , and show canted antiferromagnetic behavior below 10 K. Temperature dependent resistivity shows a Kondo-like low-temperature upturn caused by enhanced scattering of the itinerant electrons with the 4f orbitals of Eu.

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

Photoluminescence of a Uranium(IV) Alkoxide Complex

In this report, we describe the photoluminescence of a homoleptic uranium(IV) alkoxide complex. Excitation of [Li(THF)] 2 [U IV (O t Bu) 6 ] leads to the first example of photoluminescence from a well-defined actinide complex originating from an f–f excitation, supported by second order multiconfigurational electronic structure calculations including spin–orbit coupling. These calculations show strong spin–orbit coupling between the excited triplet and singlet states for the 5f-orbital manifold, which leads to a long-lived excited state lifetime of 0.85 s at low temperature. The photophysical properties of homoleptic uranium(V) and uranium(VI) tertbutoxide complexes are also presented; we find that oxidation of the uranium(IV) alkoxide results in quenching of luminescence in [Li(THF)][U V (O t Bu) 6 ] and [U VI (O t Bu) 6 ]. This is attributed to competing ligand to metal charge transfer absorption processes shifted to lower energy upon oxidation of the actinide center, which mask the relevant f–f transitions in the visible region of the electronic absorption spectrum.

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

Enhanced Menshutkin S N 2 Reactivity in Mesoporous Silica: The Influence of Surface Catalysis and Confinement

A significant enhancement in the Menshutkin S N 2 reaction between 1-methylimidazole (MeIm) and methyl thiocyanate (MeSCN) is observed when the reaction is confined in the nanoscale silica pores of MCM41 and SBA15. The experiments in the silica pores are conducted without surrounding bulk reaction mixture. The influences of temperature, pore radius, and surface chemistry on the kinetics of the confined reaction are analyzed with time-dependent infrared spectroscopy, molecular dynamics simulations, and ab initio calculations. The rate constant of the pseudo-first order reaction increases with decreasing pore size, and the activation energy is found to decrease by 5.6 kJ/mol in the smallest pore studied (2.8 nm) relative to the bulk reaction. The rate constant dependence on pore size is accurately described by a two-state model in which molecules within the 4.6 Å interfacial layer experience a 2.4-fold rate constant increase relative to those reacting at the bulk rate further away from the interface. The removal of polar silanol groups from the silica surface via passivation with trimethylsilyl chloride results in bulk-like kinetics despite a reduction in the pore diameter, demonstrating the role of silanols as catalytic sites. Here, electronic structure calculations of the energy profile on a model silica surface confirm that silanol groups, particularly those of the vicinal type, can reduce the activation energy and reaction endothermicity through the donation of hydrogen bonds to the reactant, transition state, and product complexes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

New layered quaternary Zintl pnictide oxides Ba 2 Zn 2 Pn 2 O ( Pn = Sb, Bi): Discovery, crystal structures, band engineering, and transport properties

Three new heteroanionic oxypnictides, Ba 2 Zn 2 Sb 2 O, Ba 2 Zn 2 Bi 2 O, and the solid solution Ba 2 Zn 2 Sb 2−x Bi x O (x ≈ 1.1–1.6), have been synthesized and structurally characterized. They are isostructural with their Mn-bearing analog, adopting the Ba 2 Mn 2 Sb 2 O-type structure (space group P6 3 /mmc, No. 194), and feature a double-layered 2D $^{2}_{∞}$ [Zn 2 Pn 2 O] 2- substructure (Pn = Sb, Bi, Sb/Bi) composed of corner-sharing, distorted tetrahedral ZnPn 3 O units. Electronic structure calculations reveal a systematic progression from semiconducting Ba 2 Zn 2 Sb 2 O to metallic Ba 2 Zn 2 Bi 2 O as Bi content increases. These trends are corroborated by transport property measurements, with Ba 2 Zn 2 Sb 0.9(1) Bi 1.1 O exhibiting relatively low electrical resistivity, high Hall mobilities of ∼160 cm 2 /V·s, and large Seebeck coefficients from 69 to 132 μV K −1 over the 300–600 K temperature range. Comparison with structurally related Zintl pnictides, such as SrIn 2 As 2 and PrZn 3 As 3 phases, situates Ba 2 Zn 2 Pn 2 O (Pn = Sb, Bi) within a broader family of heteroanionic oxypnictide Zintl compounds, highlighting their structural flexibility and amenability to band engineering. Finally, electronic structure and bonding considerations point to tunable semiconducting behavior and underscore the relevance of these materials for thermoelectric and topological applications.

Band engineering↗
DOE OSTI2020

Electronic structure, magnetic properties, and exchange splitting of gadolinium intermetallics

Here, we report here the electronic structure calculations of the magnetic moment, magneto-crystalline anisotropy energy (MAE), and density of states (DOS) of Gd-metal, GdAl 2 , GdAl 3 , GdCo 2 , GdCo 3 , and GdCo 5 . From the DOS, the exchange splitting energy is derived to map ferrimagnetic to paramagnetic phase transition (Curie temperature, T C ). The employed theoretical method here is based on advanced density functional theory, specifically the linearized augmented plane wave method (LAPW) within the local spin density approximation (LSDA), including spin–orbit coupling (SOC) and Hubbard model parameter (U). DOS analysis of both Gd-Al and Gd-Co systems reveals a strong correlation between exchange splitting energy and T C . The Gd-Co system exhibits an enormous exchange splitting energy among Co-3d electrons, vastly overwhelming the vanishing Gd 5d exchange splitting. Finally, these findings suggest that exchange splitting of Co-3d electrons plays a key role in the ground state magnetism and magnetic transition temperatures in these materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Anisotropic Ferromagnetism in CrAu 3 Sb 6

The crystal structure and properties of CrAu 3 Sb 6 are presented, determined by measurements on single crystal and polycrystalline samples and first-principles calculations. The trigonal structure (space group P$\overline{3}$1m) comprises a CdI 2 -like sublattice of AuSb 2 with Cr occupying octahedral holes in a fully ordered triangular array. It can be viewed as a variation of the interesting and well-known families of partially intercalated transition metal dichalcogenides, but with stronger interactions along the stacking direction evidenced by short Cr–Au distances. The compound is metallic and ferromagnetic with a Curie temperature of 164 K. A strong anomalous contribution to the Hall effect is seen in the ferromagnetic state, and quantum oscillations are observed in magnetization at 2 K. Magnetization measurements show that the ordered moments of 1.5 μ B per Cr are oriented along the c-axis with relatively strong magnetocrystalline anisotropy. Electronic structure calculations confirm this uniaxial anisotropy and the important role of spin–orbit coupling in CrAu 3 Sb 6 and reveal strongly favored ferromagnetic ground state consistent with the measured Curie temperature. Through combined experiment and theory, this work provides a detailed picture of the basic properties and behaviors of this uniquely structured, Cr-based, anisotropic ferromagnet.

Crystal structure↗
DOE OSTI2026

Complex and Unusual Excited-State Relaxation Dynamics of 9,9′-Bifluorenylidene Revealed by Comprehensive Time-Resolved Spectroscopy and MRSF-TDDFT Calculations

Design of novel photochemical molecular motors often requires molecular building blocks that exhibit rather unusual photoactivity, for which conventional analyses of spectroscopic data can lead to conflicting interpretations. We here systematically investigated the excited-state relaxation dynamics of one such molecule, 9,9′-bifluorenylidene (BF), through comprehensive and complementary integration of ultrafast transient absorption (TA) and femtosecond stimulated Raman (FSRS) spectroscopies and first-principles mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT). TA and FSRS identified two sequentially formed transients following photoexcitation. The decay kinetics of the two intermediates differ in response to excitation wavelengths and the viscosity/polarity of solvents. MRSF-TDDFT calculations reveal a direct, barrierless internal-conversion pathway from the bright Franck−Condon state to a dark S1 minimum, where the excited-state population is transiently trapped, accounting for the first transient species observed in spectroscopic experiments. Further tracking down along the PES with MRSF-TDDFT mapped out two nonradiative relaxation pathways via conical intersections that connect the dark S 1 state to three configurations in the ground-state manifolds, within which a ring structure with a C8−C8′ bond and the vibrationally excited ground-state BF were identified from spectroscopic and kinetic data. The complexity of relaxation kinetics was attributed to the flexible torsional and twisting motions about the C9−C9′ bridge bond enabled by the diradical character of the S 1 state. These findings clarify unusual photoactive relaxation dynamics stemming from a novel correlation between structural flexibility and shifting electronic characteristics, and they demonstrate the importance of integrating spectroscopic and advanced electronic structure calculation studies for judicious clarification of complex, competing relaxation pathways of excited states.

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

Effective Hamiltonians derived from equation-of-motion coupled-cluster wave functions: Theory and application to the Hubbard and Heisenberg Hamiltonians

Effective Hamiltonians, which are commonly used for fitting experimental observables, provide a coarse-grained representation of exact many-electron states obtained in quantum chemistry calculations; however, the mapping between the two is not trivial. In this contribution, we apply Bloch’s formalism to equation-of-motion coupled-cluster (EOM-CC) wave functions to rigorously derive effective Hamiltonians in the Bloch’s and des Cloizeaux’ forms. We report the key equations and illustrate the theory by application to systems with two or three unpaired electrons, which give rise to electronic states of covalent and ionic characters. We show that the Hubbard and Heisenberg Hamiltonians can be extracted directly from the so-obtained effective Hamiltonians. By establishing a quantitative connection between many-body states and simple models, the approach facilitates the analysis of the correlated wave functions. We propose a simple diagnostic for assessing the validity of the model space choice based on the overlaps between the target- and model-space states. Finally, artifacts affecting the quality of electronic structure calculations such as spin contamination are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Neutral gas pressure dependence of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations

In this computational study, we describe a self-consistent trajectory simulation approach to capture the effect of neutral gas pressure on ion–ion mutual neutralization (MN) reactions. The electron transfer probability estimated using Landau–Zener (LZ) transition state theory is incorporated into classical trajectory simulations to elicit predictions of MN cross sections in vacuum and rate constants at finite neutral gas pressures. Electronic structure calculations with multireference configuration interaction and large correlation consistent basis sets are used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of background gas on ion transport. For H + – H – and Li + – H(D) – , our approach quantitatively agrees with measured speed-dependent cross sections for up to ~10 5 m/s. For the ion pair Ne + – Cl – , our predictions of the MN rate constant at ~1 Torr are a factor of ~2 to 3 higher than the experimentally measured value. Similarly, for Xe + – F – in the pressure range of ~20 000–80 000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. Here, the paradigm of using trajectory simulations to self-consistently capture the effect of gas pressure on MN reactions advanced here provides avenues for the inclusion of additional nonclassical effects in future work.

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

NbTaV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution ternary alloys niobium-tantalum-vanadium (Nb-Ta-V). The first-principle code that has been used to run the calculations is the closed-source Vienna Ab-Initio Simulation Package (VASP). The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 8. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 10,500 randomized atomic structures over 105 chemical compositions. The calculations have been collected on NERSC-Perlmutter using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for ternary alloys is AXBYCZ, where A, B, and C refer to the constituents, and X, Y, and Z are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. 0.POSCAR: input file that defines the atomic structure of a system 3. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 4. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 6. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 7. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 8. OUTCARL: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 10. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 11. vasprun.xml: contains similar information to OUTCAR, but in an xml format. This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Exciton polaron formation and hot-carrier relaxation in rigid Dion–Jacobson-type two-dimensional perovskites

The efficiency of two-dimensional Dion–Jacobson-type materials relies on the complex interplay between electronic and lattice dynamics; however, questions remain about the functional role of exciton–phonon interactions. Here we establish the robust polaronic nature of the excitons in these materials at room temperature by combining ultrafast spectroscopy and electronic structure calculations. We show that polaronic distortion is associated with low-frequency (30–60 cm –1 ) lead iodide octahedral lattice motions. More importantly, we discover how targeted ligand modification of this two-dimensional perovskite structure manipulates exciton–phonon coupling, exciton polaron population and carrier cooling. At high excitation density, stronger exciton–phonon coupling increases the hot-carrier lifetime, forming a hot-phonon bottleneck. Furthermore, our study provides detailed insight into the exciton–phonon coupling and its role in carrier cooling in two-dimensional perovskites relevant for developing emerging hybrid semiconductor materials with tailored properties.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Colloidal Synthesis Path to 2D Crystalline Quantum Dot Superlattices

By combining colloidal nanocrystal synthesis, self-assembly, and solution phase epitaxial growth techniques, we developed a general method for preparing single dot thick atomically attached quantum dot (QD) superlattices with high-quality translational and crystallographic orientational order along with state-of-the-art uniformity in the attachment thickness. The procedure begins with colloidal synthesis of hexagonal prism shaped core/shell QDs (e.g., CdSe/CdS), followed by liquid subphase self-assembly and immobilization of superlattices on a substrate. Solution phase epitaxial growth of additional semiconductor material fills in the voids between the particles, resulting in a QD-in-matrix structure. The photoluminescence emission spectra of the QD-in-matrix structure retains characteristic 0D electronic confinement. Importantly, annealing of the resulting structures removes inhomogeneities in the QD–QD inorganic bridges, which our atomistic electronic structure calculations demonstrate would otherwise lead to Anderson-type localization. The piecewise nature of this procedure allows one to independently tune the size and material of the QD core, shell, QD–QD distance, and the matrix material. These four choices can be tuned to control many properties (degree of quantum confinement, quantum coupling, band alignments, etc.) depending on the specific applications. Lastly, cation exchange reactions can be performed on the final QD-in-matrix, as demonstrated herein with a CdSe/CdS to HgSe/HgS conversion.

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

Accelerating Embedding Potential Optimization by Reconstructing the Pseudo-Valence Electron Density

Density functional embedding theory (DFET) enables use of electronic structure methods with higher accuracy than density functional theory in a local region, with applications thus far ranging from (photo/electro)catalysis to reactions in solution. DFET partitions a large collection of atoms into smaller groups that interact via a shared embedding (interaction) potential V emb , determined via functional optimization. The optimized effective potential (OEP) process used to optimize V emb is time-consuming and becomes a computational bottleneck due to sharp, oscillating features of V emb near nuclei. Here, similar to pseudopotential theory, by reconstructing electron densities used in the OEP process from smoother pseudo-valence-only (PVO) electron densities as proxies for total densities of the full system and subsystems, we can retain accuracy in the embedded electronic structure calculations while potentially reducing the overhead of V emb construction, within the projector augmented-wave (PAW) formalism. We explore three different chemical reactions as exemplars to test PVO–DFET, namely, H 2 dissociative adsorption on a Cu(111) surface, H 2 O adsorption on a Pt(111) surface, and aqueous [Ca 2+ –SO 4 2– ] ion-pair formation. The PVO approximation works well for all three systems with minimal loss of accuracy (∼10–70 meV error relative to the original exact-derivative (ED) approach) while accelerating V emb generation for the Cu and Pt systems respectively by 20× and 5×. Given proper numerical convergence parameters, the spatial distributions of differences between PVO- and ED-based V emb outside the core regions are small, explaining the exceptional agreement between the two approaches. Finally, we anticipate that this more efficient PVO–DFET approximation will be useful whenever computation of V emb is much more expensive than subsequent embedded high-level electron correlation calculations.

approximation↗
DOE OSTI2023

Accurate calculation of many-body energies in water clusters using a classical geometry-dependent induction model

We incorporate geometry-dependent distributed multipole and polarizability surfaces into an induction model that is used to describe the 3- and 4-body terms of the interaction between water molecules. The expansion is carried out up to hexadecapole with the multipoles distributed on the atom sites. Dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole distributed polarizabilities are used to represent the response of the multipoles to an electric field. We compare the model against two large databases consisting of 43,844 3-body terms and 3,603 4-body terms obtained from high level ab initio calculations previously used to fit the MB-pol and q-AQUA interaction potentials. The classical induction model with no adjustable parameters reproduces the ab-initio 3- and 4-body terms contained in these teo Databases with a Root-Mean-Square-Error (RMSE) of 0.104/0.058 and a Mean-Absolute-Error (MAE) of 0.054/0.026 kcal/mol, respectively, results that are on a par with those obtained 1 by fitting the same data using tens of thousands of Permutationally Invariant Polynomials (PIPs). This demonstrates the accuracy of this physically motivated model in describing the 3- and 4-body terms in the interactions between water molecules with no adjustable parameters. The triple-dipole-dispersion energy was included in the 3-body energy and was found to be small but not quite negligible. The model represents a practical, efficient and transferable approach for obtaining accurate non-additive interactions for multi-component systems without the need of performing tens of thousands of high level electronic structure calculations and fitting them with tens of thousands of PIPs.

molecular analysis water, aqueous↗
DOE OSTI2022

Time-Resolved Optical Pump-Resonant X-ray Probe Spectroscopy of 4-Thiouracil: A Simulation Study

We theoretically monitor the photoinduced ππ* → nπ* internal conversion process in 4-thiouracil (4TU), triggered by an optical pump. The element-sensitive spectroscopic signatures are recorded by a resonant X-ray probe tuned to the sulfur, oxygen, or nitrogen K-edge. We employ high-level electronic structure methods optimized for core-excited electronic structure calculation combined with quantum nuclear wavepacket dynamics computed on two relevant nuclear modes, fully accounting for their quantum nature of nuclear motions. We critically discuss the capabilities and limitations of the resonant technique. For sulfur and nitrogen, we document a pre-edge spectral window free from ground-state background and rich with ππ* and nπ* absorption features. The lowest sulfur K-edge shows strong absorption for both ππ* and nπ*. In the lowest nitrogen K-edge window, we resolve a state-specific fingerprint of the ππ* and an approximate timing of the conical intersection via its depletion. A spectral signature of the nπ* transition, not accessible by UV–vis spectroscopy, is identified. The oxygen K-edge is not sensitive to molecular deformations and gives steady transient absorption features without spectral dynamics. Here, the ππ*/nπ* coherence information is masked by more intense contributions from populations. Altogether, element-specific time-resolved resonant X-ray spectroscopy provides a detailed picture of the electronic excited-state dynamics and therefore a sensitive window into the photophysics of thiobases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Toward a systematic improvement of the fixed-node approximation in diffusion Monte Carlo for solids—A case study in diamond

While Diffusion Monte Carlo (DMC) is in principle an exact stochastic method for ab initio electronic structure calculations, in practice, the fermionic sign problem necessitates the use of the fixed-node approximation and trial wavefunctions with approximate nodes (or zeros). This approximation introduces a variational error in the energy that potentially can be tested and systematically improved. Here, we present a computational method that produces trial wavefunctions with systematically improvable nodes for DMC calculations of periodic solids. These trial wavefunctions are efficiently generated with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. A simple protocol in which both exact and approximate results for finite supercells are used to extrapolate to the thermodynamic limit is introduced. This approach is illustrated in the case of the carbon diamond using Slater–Jastrow trial wavefunctions including up to one million Slater determinants. Fixed-node DMC energies obtained with such large expansions are much improved, and the fixed-node error is found to decrease monotonically and smoothly as a function of the number of determinants in the trial wavefunction, a property opening the way to a better control of this error. The cohesive energy extrapolated to the thermodynamic limit is in close agreement with the estimated experimental value. Interestingly, this is also the case at the single-determinant level, thus, indicating a very good error cancellation in carbon diamond between the bulk and atomic total fixed-node energies when using single-determinant nodes.

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

Atomic structure and Mott nature of the insulating charge density wave phase of 1T-TaS 2

Here, using x-ray pair distribution function (PDF) analysis and computer modeling, we explore structure models for the complex charge density wave (CDW) phases of layered 1T-TaS 2 that both well capture their atomic-level features and are amenable to electronic structure calculations. The models give the most probable position of constituent atoms in terms of 3D repetitive unit cells comprising a minimum number of Ta–S layers. Structure modeling results confirm the emergence of star-of-David (SD) like clusters of Ta atoms in the high-temperature incommensurate (IC) CDW phase and show that, contrary to the suggestions of recent studies, the low-temperature commensurate (C) CDW phase expands upon cooling thus reducing lattice strain. The C-CDW phase is also found to preserve the stacking sequence of Ta–S layers found in the room temperature, nearly commensurate (NC) CDW phase to a large extent. DFT based on the PDF refined model shows that bulk C-CDW 1T-TaS 2 also preserves the insulating state of individual layers of SD clusters, favoring the Mott physics description of the metal-to-insulator (NC-CDW to C-CDW) phase transition in 1T-TaS 2 . Our work highlights the importance of using precise crystal structure models in determining the nature of electronic phases in complex materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗