Search NASA⌕ Search

SEARCH · Search NASA

Results for “time dependent electronic structure theory”

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 127 records · Page 7

Computational Spectroscopy of Polycyclic Aromatic Hydrocarbons In Support of Laboratory Astrophysics

Polycyclic aromatic hydrocarbons (PAHs) are strong candidates for the molecular carriers of the unidentified infrared bands (UIR) and the diffuse interstellar bands (DIBs). In order to test the PAH hypothesis, we have systematically measured the vibronic spectra of a number of jet-cooled neutral and ionized PAHs in the near ultraviolet (UV) to visible spectral ranges using the cavity ring-down spectroscopy. To support this experimental effort, we have carried out theoretical studies of the spectra obtained in our measurements. Ab initio and (time-dependent) density.functiona1 theory calculations are performed to obtain the geometries, energetics, vibrational frequencies, transition dipole moments, and normal coordinates of these PAH molecules. Franck-Condon (FC) calculations and/or vibronic calculations are then performed using the calculated normal coordinates and vibrational frequencies to simulate the vibronic spectra. It is found that vibronic interactions in these conjugated pi electron systems are often strong enough to cause significant deviations from the Born-Oppenheimer (BO) approximation. For vibronic transitions that are well described by the BO approximation, the vibronic band profiles are simulated by calculating the rotational structure of the vibronic transitions. Vibronic oscillator strength factors are calculated in the frame of the FC approximation from the electronic transition dipole moments and the FC factors. This computational effort together with our experimental measurements provides, for the first time, powerful tools for comparison with space-based data and, hence, a powerful approach to understand the spectroscopy of interstellar PAH analogs and the nature of the UIR and DIBs.

Tan, Xiaofeng↗

QRCODE: Massively parallelized real-time time-dependent density functional theory for periodic systems

We present a new software module, QRCODE (Quantum Research for Calculating Optically Driven Excitations), for massively parallelized real-time time-dependent density functional theory (RT-TDDFT) calculations of periodic systems in the open-source Qbox software package. Our approach utilizes a custom implementation of a fast Fourier transformation scheme that significantly reduces inter-node message passing interface (MPI) communication of the major computational kernel and shows impressive scaling up to 16,344 CPU cores. In addition to improving computational performance, QRCODE contains a suite of various time propagators for accurate RT-TDDFT calculations. As benchmark applications of QRCODE, we calculate the current density and optical absorption spectra of hexagonal boron nitride (h-BN) and photo-driven reaction dynamics of the ozone-oxygen reaction. We also calculate the second and higher harmonic generation of monolayer and multi-layer boron nitride structures as examples of large material systems. Our optimized implementation of RT-TDDFT in QRCODE enables large-scale calculations of real-time electron dynamics of chemical and material systems with enhanced computational performance and impressive scaling across several thousand CPU cores.

97 MATHEMATICS AND COMPUTING↗

Dynamic Molecular Investigation of the Solid-Electrolyte Interphase of an Anode-Free Lithium Metal Battery Using In Situ Liquid SIMS and Cryo-TEM

Solid electrolyte interphase (SEI) has been widely perceived to play a critical role in the stable cycling of rechargeable batteries. However, associated with the fragile and air-sensitive nature of the SEI layer, delineation of the formation process and the nature of SEI remains a big challenge. Here, we use in situ liquid time-of-flight secondary ion mass spectroscopy (TOF-SIMS), cryo- transmission electron microscope (TEM) and density functional theory (DFT) calculation to delineate molecular process on the formation of SEI layer under the dynamic operating condition. We discover that the onset potential for SEI layer formation and the thickness of the SEI show dependence on the solvation shell structure. Using LiCoO 2 as a cathode and Cu film as an anode, the SEI is noticed to start to form at around 2.0 V and reach its final thickness (irreversible part, ~ 40-50 nm) at about 3.0 V in the 1 M LiPF 6 –EC/DMC electrolyte, while for the case of 1 M LiFSI–DME, the SEI starts to form at around 1.5 V and reaches its final thickness (~ 20 nm) at about 2.0 V. The in situ TOF-SIMS clearly indicates the outer SEI layer formation and dissipation upon charging and discharging, implying a continued evolution of electrolyte structure with extended cycling. In conclusion, the present work establishes a direct correlation between the molecular signature of SEI layer with solvation feature of electrolytes in lithium batteries, providing insights for tailoring SEI layer toward improved electrochemical properties of lithium batteries.

25 ENERGY STORAGE↗

Anomalous isotope effect on the optical bandgap in a monolayer transition metal dichalcogenide semiconductor

Isotope effects have received increasing attention in materials science and engineering because altering isotopes directly affects phonons, which can affect both thermal properties and optoelectronic properties of conventional semiconductors. However, how isotopic mass affects the optoelectronic properties in 2D semiconductors remains unclear because of measurement uncertainties resulting from sample heterogeneities. Here, we report an anomalous optical bandgap energy red shift of 13 (±7) milli–electron volts as mass of Mo isotopes is increased in laterally structured 100 MoS 2 - 92 MoS 2 monolayers grown by a two-step chemical vapor deposition that mitigates the effects of heterogeneities. This trend, which is opposite to that observed in conventional semiconductors, is explained by many-body perturbation and time-dependent density functional theories that reveal unusually large exciton binding energy renormalizations exceeding the ground-state renormalization energy due to strong coupling between confined excitons and phonons. The isotope effect on the optical bandgap reported here provides perspective on the important role of exciton-phonon coupling in the physical properties of two-dimensional materials.

36 MATERIALS SCIENCE↗

Photophysical Properties and Electronic Structure of Zinc(II) Porphyrins Bearing 0–4 meso -Phenyl Substituents: Zinc Porphine to Zinc Tetraphenylporphyrin (ZnTPP)

Six zinc(II) porphyrins bearing 0–4 meso-phenyl substituents have been examined spectroscopically and theoretically. Comparisons with previously examined free base analogues afford a deep understanding of the electronic and photophysical effects of systematic addition of phenyl groups in porphyrins containing a central zinc(II) ion versus two hydrogen atoms. Trends in the wavelengths and relative intensities of the absorption bands are generally consistent with predictions from time-dependent density functional theory calculations and simulations from Gouterman’s four-orbital model. These trends derive from a preferential effect of the meso-phenyl groups to raise the energy of the highest occupied molecular orbital. The calculations reveal additional insights, such as a progressive increase in oscillator strength in the violet-red (B-Q) absorption manifold with increasing number of phenyls. Progressive addition of 0–4 phenyl substituents to the zinc porphyrins in O 2 -free toluene engenders a reduction in the measured lifetime of the lowest singlet excited state (2.5–2.1 ns), an increase in the S 1 → S o fluorescence yield (0.022–0.030), a decrease in the yield of S 1 → T 1 intersystem crossing (0.93–0.88), and an increase in the yield of S 1 → S o internal conversion (0.048–0.090). The derived rate constants for S 1 decay reveal significant differences in the photophysical properties of the zinc chelates versus free base forms. Here, the unexpected finding of a larger rate constant for internal conversion for zinc chelates versus free bases is particularly exemplary. Collectively, the findings afford fundamental insights into the photophysical properties and electronic structure of meso-phenylporphyrins, which are widely used as benchmarks for tetrapyrrole-based architectures in solar energy and life sciences research.

14 SOLAR ENERGY↗

Molecule Origins of Dipolar Character in Excited Electronic States

The molecular level origins for symmetry breaking in the excited state of symmetrical quadrupolar molecules, particularly in polar solvents, was investigated using time-dependent density functional theory approaches. Molecules of the form ADA (A/D electron accepting/donating respectively), have been shown to break their symmetry upon excitation, producing an intramolecular charge transfer event and permanent dipole. Current research indicates that polar solvents stabilize the charge transfer event thereby producing asymmetrical solvent dynamics on opposite ends of the molecules. In this work key structural features of the molecule were identified including (1) incorporation of cyano groups, (2) rotation of grafted phenyl rings, and (3) the length of the conjugated R group chain. More specifically, incorporation of cyano groups appears to decrease the magnitude of the dipole in the excited state, thereby indicating that solvent interactions at these groups do not stabilize the charge transfer. While the rotation of the phenyl groups appears to be necessary to break the symmetry of the excited state in the molecule.

36 MATERIALS SCIENCE↗

Time Dependent Density Functional Theory Calculations of Large Compact PAH Cations: Implications for the Diffuse Interstellar Bands

We investigate the electronic absorption spectra of several maximally pericondensed polycyclic aromatic hydrocarbon radical cations with time dependent density functional theory calculations. We find interesting trends in the vertical excitation energies and oscillator strengths for this series containing pyrene through circumcoronene, the largest species containing more than 50 carbon atoms. We discuss the implications of these new results for the size and structure distribution of the diffuse interstellar band carriers.

Weisman, Jennifer L.↗

sp 3 -Functionalization of Single-Walled Carbon Nanotubes Creates Localized Spins

Chemical functionalization-introduced sp 3 quantum defects in single-walled carbon nanotubes (SWCNTs) have shown compelling optical properties for their potential applications in quantum information science and bioimaging. Here, we utilize temperature- and power-dependent electron spin resonance measurements to study the fundamental spin properties of SWCNTs functionalized with well-controlled densities of sp 3 quantum defects. Signatures of isolated spins that are highly localized at the sp 3 defect sites are observed, which we further confirm with density functional theory calculations. Applying temperature-dependent line width analysis and power-saturation measurements, we estimate the spin–lattice relaxation time T 1 and spin dephasing time T 2 to be around 9 μs and 40 ns, respectively. Furthermore, these findings of the localized spin states that are associated with the sp 3 quantum defects not only deepen our understanding of the molecular structures of the quantum defects but could also have strong implications for their applications in quantum information science.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine learning approach for vibronically renormalized electronic band structures

Here, we present a machine learning (ML) method for efficient computation of vibrational thermal expectation values of physical properties from first principles. Our approach is based on the nonperturbative frozen phonon formulation in which stochastic Monte Carlo algorithm is employed to sample configurations of nuclei in a supercell at finite temperatures based on a first-principles phonon model. A deep-learning neural network is trained to accurately predict physical properties associated with sampled phonon configurations, thus bypassing the time-consuming ab initio calculations. To incorporate the point-group symmetry of the electronic system into the ML model, group-theoretical methods are used to develop a symmetry-invariant descriptor for phonon configurations in the supercell. We apply our ML approach to compute the temperature dependent electronic energy gap of silicon based on density functional theory (DFT). We show that, with less than a hundred DFT calculations for training the neural network model, an order of magnitude larger number of sampling can be achieved for the computation of the vibrational thermal expectation values. Our work highlights the promising potential of ML techniques for finite temperature first-principles electronic structure methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Changes in polarization dictate necessary approximations for modeling electronic deexcitation intensity: Application to x-ray emission

Accurate simulation of electronic excitations and deexcitations are critical for complementing complex spectroscopic experiments and can provide validation to theoretical approaches. Here, using a generalized framework, we contrast the accuracy and validity of orbital-constrained and linear-response approaches that build upon Kohn-Sham density functional theory (DFT) to simulate emission spectra of electronic origin and propose an efficient approximation, named many-body x-ray emission spectroscopy or MBXES, for simulating such processes. We show analytically as well as with computed examples that for electronic (de)excitation leading to an appreciable change in polarization (i.e., density rearrangement), the adiabatic approximation in a response-based formalism will be inadequate for the calculation of oscillator strength. Thus, such a change (e.g., in the net electrostatic dipole moment of a finite system) can be used as a metric for evaluating the applicability of the adiabatic response-based approach and can be particularly valuable in x-ray emission spectroscopy. On the other hand, MBXES, the flexible method introduced in this paper, can compute oscillator strengths accurately at a much lower computational expense on the basis of two DFT-based self-consistent field calculations. Using illustrative examples of emission spectra, the efficacy of the MBXES method is demonstrated by comparison with its parent theory, orbital-optimized DFT, and with experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗

The Electron Thermal Conductivity of Pu and Zr Substituted Gamma-Uranium

Uranium alloys are attractive recycled nuclear fuels because of their high thermal conductivity (k) and fissile density; however, the effects of alloying elements on k remain unclear. Here, the electron thermal conductivity (k_e) of U-Pu-Zr compositions are calculated using density functional theory. The electronic structure is evaluated to understand the effects of plutonium (Pu) and zirconium (Zr) substitution on the k_e of ?-U. Alloys of up to 37.5 at. % Pu and 37.5 at. % Zr are examined. Two methods are applied to calculate k_e; we find that the accuracy of each method depends on the electronic and mass similarities between the solute and solvent atoms. Specifically, when the solute atom is similar in electronic structure and mass, the method that applies the electron relaxation time of ?-U is best, while if the elements are dissimilar, a mixed method that mixes several parameters associated with k_e from each element in the alloy is best. The introduction of all alloying elements decreases k_e; however, in binary compounds, Pu and Zr have different effects. Pu generally flattens the electronic bands but compensates for this deleterious effect by increasing electron density near the Fermi level. Zr flattens the electronic bands more severely without adding electron density near the Fermi level. Therefore, Zr decreases the k_e more than Pu in binary compounds. In ternary compounds, the difference between Pu and Zr is minimal due to the phononic change from the large mass change of Zr substitution, even at 12.5 at. %. Thus, we predict that higher loadings of Pu, and potentially other actinides, can be added to U-Pu-Zr compositions for faster recycling of spent fuel with without sacrificing k. We also note that these k_e calculation methods can be applied to non-fuel alloys that require k_e predictions, such as cladding, heat exchanger, and structural materials.

36 MATERIALS SCIENCE↗

In-situ/operando X-ray absorption spectroscopic investigation of the electrode/electrolyte interface on the molecular scale

A method for performing X-ray absorption spectroscopy (XAS) in the soft X-ray region was developed to investigate the surface-sensitive electron signal of electrode/electrolyte interfaces in common multivalent based organic solvents under in situ/operando conditions. Our approach enables us to probe the molecular-scale structure of electrode interfaces by measuring total electron yield and is suitable for redox systems exhibiting low intrinsic electrochemical current. In-situ F K-edge XAS measurements in a 0.5 M Magnesium bis(trifluoromethanesulfonimide)/2-Methyltetrahydrofuran (Mg(TFSI) 2 /2-MeTHF) electrolyte were carried out to determine the evolution of interfacial species during the electrochemical charging/discharging process. Time-dependent density-functional theory (TD-DFT) simulation indicate that the F K-edge evolution is the result of interfacial chemical environment change driven by electrochemical potential. In addition, we performed the “operando XAS” which runs the CV and collects spectra simultaneously. Furthermore, by using this method, the non-equilibrium state of condensed matter interfaces and interfacial dynamical transient process can be revealed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Electron Thermal Conductivity of Pu and Zr Substituted $\mathcal{γ}$-U

Uranium alloys are attractive recycled nuclear fuels because of their high thermal conductivity (𝑘) and fissile density. Limited experimental studies of the 𝑘 of U-Pu-Zr alloys in the range of 15 to 20 wt% Pu and 6 to 15 wt% Zr indicate that increasing the content of either Zr or Pu tends to lower 𝑘. However, which element has the greater effect on 𝑘, and the associated mechanisms, remains unclear. Here, in this study, the electron thermal conductivity (𝑘 𝑒 ) of U-Pu-Zr compositions are calculated using density functional theory. The electronic structure is evaluated to understand the effects of plutonium (Pu) and zirconium (Zr) substitution on the 𝑘 𝑒 of 𝛾-U. Alloys of up to 37.5 at. % Pu and 37.5 at. % Zr are examined. Two methods are applied to calculate 𝑘 𝑒 ; we find that the accuracy of each method depends on the electronic and mass similarities between the solute and solvent atoms. Specifically, when the solute atom is similar in electronic structure and mass, the more accurate method is that which employs the electron relaxation time of 𝛾-U, while if the elements are dissimilar, a mixed method that mixes several parameters associated with JNW_S⁢3033426825100132 from each element in the alloy is best. The introduction of all alloying elements decreases 𝑘 𝑒 ; however, in binary compounds, Pu and Zr have different effects. Pu flattens the electronic bands but compensates for this deleterious effect by increasing electron density near the Fermi level. Zr flattens the electronic bands more severely without adding electron density near the Fermi level. Therefore, Zr decreases 𝑘 𝑒 more than Pu in binary compounds. In ternary compounds, the difference between Pu and Zr is minimal due to the phononic change from the large mass change of Zr substitution, even at 12.5 at. %. Thus, we predict that higher loadings of Pu, and potentially other actinides, can be added to U-Pu-Zr compositions for faster recycling of spent fuel without sacrificing 𝑘. We also note that these 𝑘 𝑒 calculation methods can be applied to non-fuel alloys that require 𝑘 𝑒 predictions, such as cladding, heat exchanger, and structural materials.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

X-ray absorption spectroscopy and theoretical investigations of the effect of extended ligands in potassium organic matter interaction

Potassium (K) is an essential nutrient for plant growth, and despite its abundance in soil, most of the K is structurally bound in minerals, limiting its bioavailability and making this soil K reservoir largely inaccessible to plants. Microbial biochemical weathering has been shown to be a promising pathway to sustainably increase plant available K. However, the mechanisms underpinning microbial K uptake, transformation, storage, and sharing are poorly resolved. Here, to better understand the controls on microbial K transformations, we performed K K-edge x-ray absorption near-edge structure (XANES) spectroscopy on K-organic salts, including acetate, citrate, nitrate, oxalate, and tartrate, which are frequently observed as low molecular weight organic acids secreted by soil microbes, as well as humic acid, which acts as a proxy for higher molecular weight organic acids. The organic salts display feature-rich K XANES spectra, each demonstrating numerous unique features spanning ~13 eV range across the absorption edge. In contrast, the spectra for humic acid have one broad, wide feature across the same energy range. We used a combination of time-dependent density functional theory and the Bethe–Salpeter equation based approach within the OCEAN code to simulate the experimental spectra for K-nitrate (KNO 3 ) and K-citrate [K 3 (C 6 H 5 O 7 )·H 2 O] to identify the electronic transitions that give rise to some of the outlying and unique spectral features in the organic salts. KNO 3 has both the lowest and highest lying energy features, and K 3 (C 6 H 5 O 7 )·H 2 O is produced by several soil microbes and is effective at mineral weathering. Our results analyze the K-organic salt bonding in detail to elucidate why the spectral shapes differ and indicate that the K K-edge XANES spectra are associated with the entire ligand despite similar first-shell bonding environments around the K center. The improved understanding of K bonding environments with organic ligands and their use for interpretation of the K-XANES spectra provides an important toolkit to understand how K is transformed by microbial processes and made bioavailable for plant uptake.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computing molecular excited states on a D-Wave quantum annealer

Abstract The possibility of using quantum computers for electronic structure calculations has opened up a promising avenue for computational chemistry. Towards this direction, numerous algorithmic advances have been made in the last five years. The potential of quantum annealers, which are the prototypes of adiabatic quantum computers, is yet to be fully explored. In this work, we demonstrate the use of a D-Wave quantum annealer for the calculation of excited electronic states of molecular systems. These simulations play an important role in a number of areas, such as photovoltaics, semiconductor technology and nanoscience. The excited states are treated using two methods, time-dependent Hartree–Fock (TDHF) and time-dependent density-functional theory (TDDFT), both within a commonly used Tamm–Dancoff approximation (TDA). The resulting TDA eigenvalue equations are solved on a D-Wave quantum annealer using the Quantum Annealer Eigensolver (QAE), developed previously. The method is shown to reproduce a typical basis set convergence on the example $$\hbox {H}_2$$ H 2 molecule and is also applied to several other molecular species. Characteristic properties such as transition dipole moments and oscillator strengths are computed as well. Three potential energy profiles for excited states are computed for $$\hbox {NH}_3$$ NH 3 as a function of the molecular geometry. Similar to previous studies, the accuracy of the method is dependent on the accuracy of the intermediate meta-heuristic software called qbsolv.

74 ATOMIC AND MOLECULAR PHYSICS↗

Density Matrix Implementation of the Fermi–Löwdin Orbital Self-Interaction Correction Method

The Fermi–Löwdin orbital self-interaction correction (FLOSIC) method effectively provides a transformation from canonical orbitals to localized Fermi–Löwdin orbitals which are used to remove the self-interaction error in the Perdew–Zunger (PZ) framework. This transformation is solely determined by a set of points in space, called Fermi–Löwdin descriptors (FODs), and the occupied canonical orbitals or the density matrix. In this work, we provide a detailed workflow for the implementation of the FLOSIC method for removal of self-interaction error in DFT calculations in an orbital-by-orbital basis that takes advantage of the unitary invariant nature of the FLOSIC method. In this way, it is possible to cast the self-consistent energy minimization at fixed FODs in the same manner than standard Kohn–Sham with one additional term in the Kohn–Sham Hamiltonian that introduces the PZ self-interaction correction. Each energy minimization iteration is divided in two substeps, one for the density matrix and one for the FODs. Expressions for the effective Kohn–Sham matrix and FOD gradients are provided such that its implementation is suitable for most electronic structure codes. Here, we analyze the convergence characteristics of the algorithm and present applications for the evaluation of NMR shielding constants and real-time time-dependent DFT simulations based on the Liouville–von Neumann equation to calculate excitation energies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulating excited-state complex ensembles: Fluorescence and solvatochromism in amine-arene exciplexes

Exciplexes are excited-state complexes formed as a result of partial charge transfer from the donor to the acceptor species when one moiety of the donor–acceptor pair is electronically excited. The arene–amine exciplex formed between oligo-(p-phenylene) (OPP) and triethylamine (TEA) is of interest in the catalytic photoreduction of CO2 because it can compete with complete electron transfer to the OPP catalyst. Therefore, formation of the exciplex can hinder the generation of a radical anion OPP·− necessary for subsequent CO2 reduction. We report an implementation of a workflow automating quantum-chemistry calculations that generate and characterize an ensemble of structures to represent this exciplex state. We use FireWorks, Pymatgen, and Custodian Python packages for high-throughput ensemble generation. The workflow includes time-dependent density functional theory optimization, verification of excited-state minima, and exciplex characterization with natural transition orbitals, exciton analysis, excited-state Mulliken charges, and energy decomposition analysis. Fluorescence spectra computed for these ensembles using Boltzmann-weighted contributions of each structure agree better with experiment than our previous calculations based on a single representative exciplex structure [Kron et al., J. Phys. Chem. A 126, 2319–2329 (2022)]. The ensemble description of the exciplex state also reproduces an experimentally observed red shift of the emission spectrum of [OPP-4–TEA]* relative to [OPP-3–TEA]*. The workflow developed here streamlines otherwise labor-intensive calculations that would require significant user involvement and intervention.

Chemistry↗

Symmetry Breaking in the Lowest-Lying Excited-State of CCl4: Valence Shell Spectroscopy in the 5.0–10.8 eV Photon Energy Range

We report absolute high-resolution vacuum ultraviolet (VUV) photoabsorption cross-sections of carbon tetrachloride (CCl4) in the photon energy range 5.0–10.8 eV (248–115 nm). The molecular spectrum and electronic structure have been comprehensively investigated together with quantum chemical calculations, providing geometries, bond lengths, vertical excitation energies and oscillator strengths. The major electronic excitations have been assigned to valence and Rydberg transitions which are also accompanied by vibrational excitation assigned to degenerate stretching, v3′t2 and degenerate deformation v4′t2 modes. The rather complex nuclear dynamics along the degenerate deformation mode, v4′t2, have been thoroughly investigated by Time-Dependent Density Functional Theory (TD-DFT) method. The relevant Jahn–Teller distortion operative within the lowest-lying electronic excited-state is shown here for the first time in order to yield a weak absorption feature at 6.156 eV. Further calculations on the potential energy curves for the singlet excited-states along the C–Cl stretching coordinate show the relevance of efficient C–Cl bond excision.

Biochemistry & Molecular Biology↗