Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TDDFT”

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 55 records · Page 3

Data-Driven Recommendation of Optimal Tuning Scheme for Range-Separated Hybrid Functionals in Solution-Phase UV/Vis Absorption Energy Prediction

Time-dependent density functional theory (TDDFT) combined with range-separated hybrid (RSH) functionals and a tuned range-separation parameter γ offers a computationally economical approach for high-throughput excited- state property predictions. The γ-tuning procedure in the gas phase is well established. However, no agreement on the best γ- tuning procedure has been made when considering the solvent effect with implicit solvent models like the polarizable continuum model (PCM). To answer that question, this study created a diverse dataset with 937 molecules with experimental solutionphase UV/vis absorption spectra. Three γ-tuning methods, the gasphase γ-tuning (GPγT), the partial vertical γ-tuning (PVγT), and the strict vertical γ-tuning (SVγT), were evaluated for the ωPBEh functional over the entire dataset. Additional benchmarks are done for the optimally tuned screened range-separated hybrid combined with the PCM approach (SRSH-PCM) and the solvation-mediated tuning procedure (sol-med-OT). Our findings revealed that the optimal γ-values obtained by the PVγT and the SVγT are significantly smaller than the GPγT. This trend holds consistently across all molecules in our dataset, and we explained the origin of this phenomenon. TDDFT calculations with PVγTand SVγT-tuned γ-values and default global Fock exchange fraction achieve superior performance compared to those using GPγTtuned or default γ and slightly outperform SRSH-PCM and sol-med-OT with similar or lesser computational cost. Furthermore, we found that the smaller γ-values from SVγT captured the expected 1/(εR) asymptotic behavior in the solution phase, resulting in accurate prediction of solution-phase CT excitations, consistent with the screened asymptote behavior encoded in SRSH-PCM. These results show that SVγT is the best scheme for high-throughput UV/vis absorption spectrum calculations using the ωPBEh functional from a data-driven perspective.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electron-Beam-Induced Molecular Plasmon Excitation and Energy Transfer in Silver Molecular Nanowires

In this work, we investigate (1) electron-beam-induced plasmon absorption spectra of Ag molecular nanowire dimers and (2) electron-beam-induced energy transfer between two nanowires. We employ linear-response time-dependent density functional theory (TDDFT) and real-time TDDFT methods to simulate the electron-beam-induced plasmonic excitations, dynamics, and corresponding electron energy loss spectrum for small models of a single molecular nanowire with four Ag atoms and for two Ag nanowires. An array of different relative orientations of nanowires and of different initial excitation conditions resulting from applying an electron beam at different positions with respect to the Ag nanowires is investigated. The results demonstrate (1) an electron beam can induce plasmonic excitations from the molecular Ag nanowire ground state to the excited states that are both optically allowed and forbidden, (2) a tunability for selective excitations that can be controlled by the position of a focused electron beam, and (3) kinetic and dynamic behaviors of time-dependent electron-beam-induced energy transfer between two Ag molecular nanowires depend on the position of the beam source and nanowire separation distance, providing insights into the spatial dependences of plasmonic couplings in nanowire arrays.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Second Linear Response Theory and the Analytic Calculation of Excited-State Properties

We present a method based on second linear response time-dependent density func- tional theory (TDDFT) to calculate permanent and transition multipoles of excited states, which are required to compute excited-state absorption/emission spectra, multi- photon optical processes, among others. In previous work, we examined computations based on second linear response theory in which linear response TDDFT was employed twice. In contrast, the present methodology only requires information from a single linear response calculation to compute the excited-state properties. These are eval- uated analytically through various algebraic operations involving electron repulsion integrals and excitation vectors. The present derivation focuses on full many-body wave functions instead of single orbitals, as in our previous approach. Here, we test the proposed method by applying it to several diatomic and triatomic molecules. This shows that the computed excited-state dipoles are consistent with respect to reference equation-of-motion coupled cluster calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical Investigations on the Plasmon-Mediated Dissociation of Small Molecules in the Presence of Silver Atomic Wires

Plasmonic nanoparticles can promote bond activation in adsorbed molecules under relatively benign conditions via excitation of the nanoparticle’s plasmon resonance. As the plasmon resonance often falls within the visible light region, plasmonic nanomaterials are a promising class of catalysts. However, the exact mechanisms through which plasmonic nanoparticles activate the bonds of nearby molecules are still unclear. Herein, we evaluate Ag 8 –X 2 (X = N, H) model systems via real-time time-dependent density functional theory (RT-TDDFT), linear response time-dependent density functional theory (LR-TDDFT), and Ehrenfest dynamics in order to better understand the bond activation processes of N 2 and H 2 facilitated by the presence of the atomic silver wire under excitation at the plasmon resonance energies. We find that dissociation is possible for both small molecules at high electric field strength. Activation of each adsorbate is symmetry- and electric field-dependent, and H 2 activates at lower electric field strengths than N 2 . Furthermore, this work serves as a step toward understanding the complex time-dependent electron and electron–nuclear dynamics between plasmonic nanowires and adsorbed small molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Excited State Dynamics of a Conformationally Fluxional Copper Coordination Complex

The conversion of solar energy into chemical fuel represents a capstone goal of the 21st century and has the potential to supply terawatts of power in a globally distributed manner. However, the disparate time scales of photodriven charge separation (~fs) and steps in chemical reactions (~μs) represent an inherent bottleneck in solar-to-fuels technology. Here to address this discrepancy, we are developing earth-abundant coordination complexes that undergo light-induced conformational rearrangements such that charge separation (CS) is hastened, while charge recombination (CR) is slowed. To these ends, we report the preparation and characterization of a new series of conformationally fluxional copper coordination complexes that contain a twisted intramolecular charge transfer (TICT) fluorophore as part of their ligand scaffold. Structural and spectroscopic characterization of the Cu(I) and Cu(II) complexes formed with these ligands in their ground states establish oxidation state-dependent conformational dynamicity, while time-resolved emission and transient absorption spectroscopies define the photopysical parameters of photo-induced excited states. Building on initial reports with a related set of molecules, the improved ligand design presented here greatly simplifies the observed photophysics, effectively shutting down unwanted ligand-centered excited states previously observed. Finally, time-dependent density functional theory (TDDFT) analyses reveal an unusual metal-to-TICT electronic transition only reported once before, and though the formation of a CS state is not observed directly through experiments, TDDFT geometry optimizations in the excited states support the formation of transient Cu(II) CS species, lending credence to the potential success of our approach. These studies establish a clear model for the excited state dynamics at play in proof-of-concept systems and clarify key design parameters for future optimizations toward achieving long-lived CS via photoinduced conformational gating.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Near-Quantitative Predictions of the First-Shell Coordination Structure of Hydrated First-Row Transition Metal Ions Using K-Edge X-ray Absorption Near-Edge Spectroscopy

Understanding the solvation structure of transition metal ions is important due to its broad impact on the kinetics, stability, and reactivity in a variety of applications in geochemistry, biochemistry, energy storage, and environmental chemistry. Using water as the common ligand, we study the X-ray absorption pre-edge and near-edge spectra at the K-edge of a near-complete series of hydrated first-row transition metal ions with d-orbital occupancy ranging from d 2 to d 10 . Starting with optimized structures that were derived from an explicit solvation treatment at the density functional theory (DFT) level, we then compute the pre-edge X-ray absorption spectra at the metal ion K-edge at the time dependent density functional theory (TDDFT) and restricted active-space second-order perturbation theory (RASPT2) levels of theory. TDDFT calculations provide accurate results for spectra that are dominated by single excitations, while significant improvements were obtained with RASPT2 calculations that correctly distinguish between singly and doubly excited states, where relevant, with quantitative accuracy compared with experiment. We analyze and assign the dierent pre-edge features for each metal ion in order to reveal the impact of the variations in the d orbital occupancy on the rst-shell coordination environment. For completeness, we also report the lowest energy ligand-field d-d transitions for all the transition metal aqua ions considered in this study using complete active space second order perturbation theory (CASPT2).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Combined Wave Function and Density Functional Approach for K-Edge X-ray Absorption Near-Edge Spectroscopy: A Case Study of Hydrated First-Row Transition Metal Ions

The prediction of X-ray absorption spectra (XAS) of transition metal complexes has important and broad application areas in chemistry and biology. In this letter, we have investigated the predictive ability of multiconfiguration pair-density functional theory (MC-PDFT) for X-ray absorption spectra by calculating the metal K pre-edge features of aquated 3d transition metal ions in common oxidation states. MC-PDFT results were compared with experimentally measured spectra as well as analyzed against results from restricted active-space second-order perturbation theory (RASPT2) and time-dependent density functional theory (TDDFT). As expected, TDDFT performs well for excited states that can be accurately represented by singly excited configurations but fails for excited states where higher order excitations become important. On the other hand, both RASPT2 and MC-PDFT provide quantitatively accurate results for all excited states irrespective of their character. While core-level spectroscopy with RASPT2 is accurate, it is computationally expensive. Finally, our results show that MC-PDFT performs equally well with significantly lower computational cost and is an encouraging alternate approach for X-ray spectroscopies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nuclear–Electronic Orbital Quantum Mechanical/Molecular Mechanical Real-Time Dynamics

Simulating the nuclear–electronic quantum dynamics of large-scale molecular systems in the condensed phase is key for studying biologically and chemically important processes such as proton transfer and proton-coupled electron transfer reactions. Herein, the real-time nuclear–electronic orbital time-dependent density functional theory (RT-NEO-TDDFT) approach is combined with a hybrid quantum mechanical/molecular mechanical (QM/MM) strategy to enable the accurate description of coupled nuclear–electronic quantum dynamics in the presence of heterogeneous environments such as solvent or proteins. The densities of the electrons and quantum protons are propagated in real time, while the other nuclei are propagated classically on the instantaneous electron–proton vibronic surface. This approach is applied to phenol bound to lysozyme, intramolecular proton transfer in malonaldehyde, and nonequilibrium excited-state intramolecular proton transfer in o-hydroxybenzaldehyde. Furthermore, these examples illustrate that the RT-NEO-TDDFT framework, coupled with an atomistic representation of the environment, allows the simulation of condensed-phase systems that exhibit significant nuclear quantum effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Strongly Correlated States of Transition Metal Spin Defects: The Case of an Iron Impurity in Aluminum Nitride

We investigate the electronic properties of an exemplar transition metal impurity in an insulator, with the goal of accurately describing strongly correlated defect states. Here, we consider iron in aluminum nitride, a material of interest for hybrid quantum technologies, and we carry out calculations with quantum embedding methods, density matrix embedding theory (DMET) and quantum defect embedding theory (QDET), and with spin-flip time-dependent density functional theory (TDDFT). We show that both DMET and QDET accurately describe the ground state and low-lying excited states of the defect and that TDDFT yields photoluminescence spectra in agreement with experiments. In addition, we provide a detailed discussion of the convergence of our results as a function of the active space used in the embedding methods, thus defining a protocol to obtain converged data directly comparable with experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Trajectory sampling and finite-size effects in first-principles stopping power calculations

Abstract Real-time time-dependent density functional theory (TDDFT) is presently the most accurate available method for computing electronic stopping powers from first principles. However, obtaining application-relevant results often involves either costly averages over multiple calculations or ad hoc selection of a representative ion trajectory. We consider a broadly applicable, quantitative metric for evaluating and optimizing trajectories in this context. This methodology enables rigorous analysis of the failure modes of various common trajectory choices in crystalline materials. Although randomly selecting trajectories is common practice in stopping power calculations in solids, we show that nearly 30% of random trajectories in an FCC aluminum crystal will not representatively sample the material over the time and length scales feasibly simulated with TDDFT, and unrepresentative choices incur errors of up to 60%. We also show that finite-size effects depend on ion trajectory via “ouroboros” effects beyond the prevailing plasmon-based interpretation, and we propose a cost-reducing scheme to obtain converged results even when expensive core-electron contributions preclude large supercells. This work helps to mitigate poorly controlled approximations in first-principles stopping power calculations, allowing 1–2 order of magnitude cost reductions for obtaining representatively averaged and converged results.

36 MATERIALS SCIENCE↗

Efficient simulations of charge density waves in the transition metal Dichalcogenide TiSe 2

Charge density waves (CDWs) in transition metal dichalcogenides are the subject of growing scientific interest due to their rich interplay with exotic phases of matter and their potential technological applications. Here, using density functional theory with advanced meta-generalized gradient approximations (meta-GGAs) and linear response time-dependent density functional theory (TDDFT) with state-of-the-art exchange-correlation kernels, we investigate the electronic, vibrational, and optical properties in 1T-TiSe 2 with and without CDW. In both bulk and monolayer TiSe 2 , the electronic bands and phonon dispersions in either normal or CDW (semiconducting) phase are described well via meta-GGAs, which separate the valence and conduction bands just as HSE06 does but with significantly more computational feasibility. The experimentally observed humps of electron energy loss spectroscopy are successfully reproduced in TDDFT. Our work opens the door to simulating these complexities in CDW compounds from first principles by revealing meta-GGAs as an accurate low-cost alternative to HSE06.

Electronic properties and materials↗

Analysis and visualization of energy densities. I. Insights from real-time time-dependent density functional theory simulations

In this article, we report a scheme to analyze and visualize the energy density fluctuations during the real-time time-dependent density functional theory (RT-TDDFT) simulations. Using Ag 4 –N 2 complexes as examples, here we show that the grid-based Kohn–Sham energy density can be computed at each time step using a procedure from Nakai and coworkers. Then the instantaneous energy of each molecular fragment (such as Ag 4 and N 2 ) can be obtained by partitioning the Kohn–Sham energy densities using Becke or fragment-based Hirshfeld (FBH) scheme. A strong orientation-dependence is observed for the energy flow between the Ag 4 cluster and a nearby N 2 molecule in the RT-TDDFT simulations. Future applications of such an energy density analysis in electron dynamics simulations are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Valence-to-core X-ray emission spectroscopy of vanadium oxide and lithiated vanadyl phosphate materials

In this work, we report valence-to-core (VTC) X-ray emission spectroscopy (XES) measurements and theoretical calculations of the electrochemical sequence ε-VOPO 4 , ε-LiVOPO 4 , ε-Li 2 VOPO 4 and the reference oxides V 2 O 3 , VO 2 , and V 2 O 5 . In our analysis of these results, we establish a framework for interrogating chemical bonding that is generally applicable to a wide range of systems, including complex, extended inorganic compounds. While this latter regime has garnered less focused application than, e.g. , metalloenzymes in many excellent catalysis studies, we show that the technique provides high utility in materials-focused energy storage research. Here, sensitivities to the local atomic structure and hybridization schemes are discussed in detail. Similarly, the effect of lithiation on oxidation, delocalization, and shifts in ligand valence energy levels are all readily apparent in the analyzed results. Finally, the TDDFT projections clearly reveal the directional dependencies of the valence band at each of the vanadium sites. Our results demonstrate laboratory-based X-ray spectroscopy instrumentation is a viable route for attaining well-resolved VTC-XES features for inorganic compounds of 3d transition metals, even for samples of limited quantity or suffering from sensitivity to the atmosphere. The experimental results are in good agreement with results produced by real-space Green's function and time-dependent density functional theory (TDDFT) methods, respectively. Hence, we propose that VTC-XES, when equipped with appropriate theoretical support, can be a valuable complement to X-ray absorption pre-edge features for more detailed characterization of a compound's electronic structure. We expect similar analyses will find application in a broad range of materials chemistry research and provide both fundamental and applied insights.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Microscopic evidence for scission neutrons

We discuss the neck rupture stage of fission and the emission of particles afterwards, two extremely rapid and highly non-equilibrium processes. Currently, the neck rupture cannot be directly probed by experiment, highlighting the importance of reliable theoretical predictions for this stage of fission. Time-dependent density functional theory (TDDFT) is used to simulate the spontaneous fission of 252 Cf. In conjunction with statistical models, inputs from microscopic calculations can be used to make predictions for prompt neutron and gamma emission spectra, quantities which can be experimentally measured. Here we characterize the scission mechanism within TDDFT and estimate the number of scission neutrons and their kinetic energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Exploring non-adiabaticity to CO reduction reaction through ab initio molecular dynamics simulation

Non-adiabatic chemical reaction refers to the electronic excitation during reactions. This effect cannot be modeled by the ground-state Born-Oppenheimer molecular dynamics (BO-MD), where the electronic structure is at the ground state for every step of ions' movement. Although the non-adiabatic effect has been explored extensively in gas phase reactions, its role in electrochemical reactions, such as water splitting and CO 2 reduction, in electrolyte has been rarely explored. On the other hand, electrochemical reactions usually involve electron transport; thus, a non-adiabatic process can naturally play a significant role. In this work, using one-step CO 2 reduction as an example, we investigated the role of the non-adiabatic effect in the reaction. The reaction barriers were computed by adiabatic BO-MD and non-adiabatic real-time time dependent density functional theory (rt-TDDFT). We found that by including the non-adiabatic effect, rt-TDDFT could increase the reaction barrier up to 6% compared to the BO-MD calculated barrier when the solvent model is used to represent water. Simulations were carried out using explicit water molecules around the reaction site under different overpotentials, and similar non-adiabatic effects were found.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hole–hole Tamm–Dancoff-approximated density functional theory: A highly efficient electronic structure method incorporating dynamic and static correlation

The study of photochemical reaction dynamics requires accurate as well as computationally efficient electronic structure methods for the ground and excited states. While time-dependent density functional theory (TDDFT) is not able to capture static correlation, complete active space self-consistent field methods neglect much of the dynamic correlation. Hence, inexpensive methods that encompass both static and dynamic electron correlation effects are of high interest. Here, we revisit hole–hole Tamm–Dancoff approximated (hh-TDA) density functional theory for this purpose. The hh-TDA method is the hole–hole counterpart to the more established particle–particle TDA (pp-TDA) method, both of which are derived from the particle–particle random phase approximation (pp-RPA). In hh-TDA, the N-electron electronic states are obtained through double annihilations starting from a doubly anionic (N+2 electron) reference state. In this way, hh-TDA treats ground and excited states on equal footing, thus allowing for conical intersections to be correctly described. Furthermore, the treatment of dynamic correlation is introduced through the use of commonly employed density functional approximations to the exchange-correlation potential. Additionally, we show that hh-TDA is a promising candidate to efficiently treat the photochemistry of organic and biochemical systems that involve several low-lying excited states—particularly those with both low-lying ππ* and nπ* states where inclusion of dynamic correlation is essential to describe the relative energetics. In contrast to the existing literature on pp-TDA and pp-RPA, we employ a functional-dependent choice for the response kernel in pp- and hh-TDA, which closely resembles the response kernels occurring in linear response and collinear spin-flip TDDFT.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An optimally tuned range-separated hybrid starting point for ab initio GW plus Bethe–Salpeter equation calculations of molecules

The ab initio GW plus Bethe–Salpeter equation (GW-BSE, where G is the one particle Green's function and W is the screened Coulomb interaction) approach has emerged as a leading method for predicting excitations in both solids and molecules with a predictive power contingent upon several factors. Among these factors are the (1) generalized Kohn–Sham eigensystem used to construct the GW self-energy and to solve the BSE and (2) the efficacy and suitability of the Tamm–Dancoff approximation. Here, we present a detailed benchmark study of low-lying singlet excitations from a generalized Kohn–Sham (gKS) starting point based on an optimally tuned range-separated hybrid (OTRSH) functional. We show that the use of this gKS starting point with one-shot G0W0 and G0W0-BSE leads to the lowest mean absolute errors (MAEs) and mean signed errors (MSEs), with respect to high-accuracy reference values, demonstrated in the literature thus far for the ionization potentials of the GW100 benchmark set and for low-lying neutral excitations of Thiel’s set molecules in the gas phase, without the need for self-consistency. The MSEs and MAEs of one-shot G0W0-BSE@OTRSH excitation energies are comparable to or lower than those obtained with other functional starting points after self-consistency. Additionally, we compare these results with linear-response time-dependent density functional theory (TDDFT) calculations and find GW-BSE to be superior to TDDFT when calculations are based on the same exchange-correlation functional. This work demonstrates tuned range-separated hybrids used in combination with GW and GW-BSE can greatly suppress starting point dependence for molecules, leading to accuracy similar to that for higher-order wavefunction-based theories for molecules without the need for costlier iterations to self-consistency.

McKeon, Caroline A. (ORCID:0000000217373503)↗

Importance of imposing gauge invariance in time-dependent density functional theory calculations with meta-generalized gradient approximations

It has been known for more than a decade that the gauge variance of the kinetic energy density τ leads to additional terms in the magnetic orbital rotation Hessian used in linear-response time-dependent density functional theory (TDDFT), affecting excitation energies obtained with τ-dependent exchange–correlation functionals. While previous investigations found that a correction scheme based on the paramagnetic current density has a small effect on benchmark results, we report more pronounced effects here, in particular, for the popular M06-2X functional and for some other meta-generalized gradient approximations (mGGAs). In the first part of this communication, this is shown by a reassessment of a set of five Ni(II) complexes for which a previous benchmark study that did not impose gauge invariance has found surprisingly large errors for excitation energies obtained with M06-2X. These errors are more than halved by restoring gauge invariance. The variable importance of imposing gauge invariance for different mGGA-based functionals can be rationalized by the derivative of the mGGA exchange energy integrand with respect to τ. In the second part, a large set of valence excitations in small main-group molecules is analyzed. For M06-2X, several selected n → π* and π→π$^{*}_{⊥}$ excitations are heavily gauge-dependent with average changes of –0.17 and –0.28 eV, respectively, while π→π$^{*}_{∥}$ excitations are marginally affected (–0.04 eV). Similar patterns, but of the opposite signs, are found for SCAN0. Here, the results suggest that reevaluation of previous gauge variant TDDFT results based on M06-2X and other mGGA functionals is warranted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗