Search NASA⌕ Search

SEARCH · Search NASA

Results for “Wave Functions”

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 1,639 records · Page 91

The HANDE-QMC Project: Open-Source Stochastic Quantum Chemistry from the Ground State Up

Building on the success of Quantum Monte Carlo techniques such as diffusion Monte Carlo, alternative stochastic approaches to solve electronic structure problems have emerged over the past decade. The full configuration interaction quantum Monte Carlo (FCIQMC) method allows one to systematically approach the exact solution of such problems, for cases where very high accuracy is desired. The introduction of FCIQMC has subsequently led to the development of coupled cluster Monte Carlo (CCMC) and density matrix quantum Monte Carlo (DMQMC), allowing stochastic sampling of the coupled cluster wave function and the exact thermal density matrix, respectively. In this Article, we describe the HANDE-QMC code, an open-source implementation of FCIQMC, CCMC and DMQMC, including initiator and semistochastic adaptations. We describe our code and demonstrate its use on three example systems; a molecule (nitric oxide), a model solid (the uniform electron gas), and a real solid (diamond). An illustrative tutorial is also included.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Robust, Accurate, and Efficient: Quantum Embedding Using the Huzinaga Level-Shift Projection Operator for Complex Systems

Using wave function (WF) in density functional theory (DFT) embedding methods provides a framework for performing localized, high-accuracy WF calculations on a system, while not incurring the full computational cost of the WF calculation on the full system. To effectively partition a system into localized WF and DFT subsystems, we utilize the Huzinaga level-shift projection operator within an absolutely localized basis. In this work, we study the ability of the absolutely localized Huzinaga level-shift projection operator method to study complex WF and DFT partitions, including partitions between multiple covalent bonds, a double bond, and transition-metal–ligand bonds. We find that our methodology can accurately describe all of these complex partitions. Additionally, we study the robustness of this method with respect to the WF method, specifically where the embedded systems were described using a multiconfigurational WF method. We found that the method is systematically improvable with respect to both the number of atoms in the WF region and the size of the basis set used, with energy errors less than 1 kcal/mol. Additionally, we calculated the adsorption energy of H2 to a model of an iron metal–organic framework (Fe-MOF-74) to within 1 kcal/mol compared to CASPT2 calculations performed on the full model while incurring only a small fraction of the full computational cost. In conclusion, this work demonstrates that the absolutely localized Huzinaga level-shift projection operator method is applicable to very complex systems with difficult electronic structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Improving Excited-State Potential Energy Surfaces via Optimal Orbital Shapes

In this work, we demonstrate that, rather than resorting to high-cost dynamic correlation methods, qualitative failures in excited-state potential energy surface predictions can often be remedied at no additional cost by ensuring that optimal molecular orbitals are used for each individual excited state. This approach also avoids the weighting choices required by state-averaging and dynamic weighting and obviates their need for expensive wave function response calculations when relaxing excited-state geometries. Although multistate approaches are of course preferred near conical intersections, other features of excited-state potential energy surfaces can benefit significantly from our single-state approach. In three different systems, including a double bond dissociation, a biologically relevant amino hydrogen dissociation, and an amino-to-ring intramolecular charge transfer, we show that state-specific orbitals offer qualitative improvements over the state-averaged status quo.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Toward Full Configuration Interaction for Transition-Metal Complexes

In this work, an efficacious approximation to full configuration interaction (FCI) is adapted to calculate singlet-triplet gaps for transition-metal complexes. This strategy, incremental FCI (iFCI), uses a many-body expansion to systematically add correlation to a simple reference wave function and therefore achieves greatly reduced computational costs compared to FCI. iFCI through the 3-body expansion is demonstrated on four model transition-metal complexes involving the metals Zn, V, and Cu. Screening techniques to increase the computational efficiency of iFCI are proposed and tested, showing reduction in the number of 3-body terms by more than 90% with controlled errors. The largest complex treated by iFCI has 142 valence electrons, all of which are correlated among the full set of 444 active orbitals. Computed spin gaps approach experimental results for the four complexes, though room for improvement remains.

74 ATOMIC AND MOLECULAR PHYSICS↗

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↗

Electronic Structure of the Complete Series of Gas-Phase Manganese Acetylacetonates by X-ray Absorption Spectroscopy

Metal centers in transition metal–ligand complexes occur in a variety of oxidation states causing their redox activity and therefore making them relevant for applications in physics and chemistry. The electronic state of these complexes can be studied by X-ray absorption spectroscopy, which is, however, due to the complex spectral signature not always straightforward. Here, we study the electronic structure of gas-phase cationic manganese acetylacetonate complexes Mn(acac) 1–3 + using X-ray absorption spectroscopy at the metal center and ligand constituents. The spectra are well reproduced by multiconfigurational wave function theory, time-dependent density functional theory as well as parameterized crystal field and charge transfer multiplet simulations. This enables us to get detailed insights into the electronic structure of ground-state Mn(acac) 1–3 + and extract empirical parameters such as crystal field strength and exchange coupling from X-ray excitation at both the metal and ligand sites. By comparison to X-ray absorption spectra of neutral, solvated Mn(acac) 2,3 complexes, we also show that the effect of coordination on the L 3 excitation energy, routinely used to identify oxidation states, can contribute about 40–50% to the observed shift, which for the current study is 1.9 eV per oxidation state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Chemist: A Domain-Specific Language by Chemists for Chemists

Managing the complexity of quantum chemistry (QC) software is key to ensuring it remains accessible, maintainable, and reusable. Noticeably missing from the QC ecosystem are modules targeting bottleneck routines. Here we argue that this is likely due to the difficulty in defining interfaces for such modules. To that end, we introduce the open-source, publicly available Chemist library https://github.com/NWChemEx/Chemist. Chemist is a domain-specific language targeting the QC domain. Chemist has been developed focusing on performance and user-friendliness. Using Chemist, QC tasks are defined using familiar domain concepts such as molecules, wave functions, and operators. The domain objects are hierarchical to ensure a systematic encapsulation of information. Key features of Chemist include: extensibility, the ability to alias existing data, and the ability to succinctly define many common QC tasks. The usefulness of Chemist is demonstrated by discussing the interface of NWChemEx’s Fock build module and by showcasing a proof-of-concept self-consistent field algorithm containing uncertainty propagation.

Algorithms↗

Nonunitary Variational Quantum Eigensolver with the Localized Active Space Method and Cost Mitigation

Accurately describing strongly correlated systems with affordable quantum resources remains a central challenge for quantum chemistry applications on near and intermediate term quantum computers. The localized active space self-consistent field (LASSCF) approximates the complete active space self-consistent field (CASSCF) by generating active space-based wave functions within specific fragments while treating interfragment correlation with mean-field approach, hence is computationally less expensive. Hardware-efficient ansatzes (HEA) offer affordable and shallower circuits, yet they often fail to capture the necessary correlation. Previously, Jastrow-factor-inspired nonunitary qubit operators were proposed to use with HEA for variational quantum eigensolver (VQE) calculations (so-called nuVQE), as they do not increase circuit depths and recover correlation beyond the mean-field level for Hartree–Fock initial states. Here, in this study, we explore running nuVQE with LASSCF as the initial state. The method, named LAS-nuVQE, is shown to recover interfragment correlations, reach chemical accuracy with a small number of gates (<70) in both H 4 and square cyclobutadiene (C 4 H 4 ), and produces more accurate energetics than its HEA counterparts at all circuit depths. To further address the inherent symmetry-breaking in HEA, we implemented spin-constrained LAS-nuVQE to extend the capabilities of HEA further and show spin-pure results for square cyclobutadiene. We also mitigate the increased measurement overhead of nuVQE via Pauli grouping and shot-frugal sampling, reducing measurement costs by up to 2 orders of magnitude compared to ungrouped operator, and show that one can achieve better accuracy with a small number of shots (10 3–4 ) per one expectation value calculation compared to noiseless simulations with one or two orders of magnitude more shots. Finally, wall clock time estimates show that, with our measurement mitigation protocols, nuVQE becomes a cheaper and more accurate alternative than vanilla VQE with HEA. Taken together, these developments illustrate a practical pathway toward performing multireference chemical simulations with accuracy and affordable resources on today’s quantum hardware, achieving both accuracy and affordability in challenging correlated systems.

Wang, Qiaohong [Univ. of Chicago, IL (United State↗

Effect of Anisotropic Confinement on Electronic Structure and Dynamics of Band Edge Excitons in Inorganic Perovskite Nanowires

Inorganic lead halide perovskite nanostructures show promise as the active layers in photovoltaics, light emitting diodes, and other optoelectronic devices. They are robust in the presence of oxygen and water, and the electronic structure and dynamics of these nanostructures can be tuned through quantum confinement. Here we create aligned bundles of CsPbBr3 nanowires with widths resulting in quantum confinement of the electronic wave functions and subject them to ultrafast microscopy. We directly image rapid one-dimensional exciton diffusion along the nanowires, and we measure an exciton trap density of roughly one per nanowire. Using transient absorption microscopy, we observe a polarization-dependent splitting of the band edge exciton line, and from the polarized fluorescence of nanowires in solution, we determine that the exciton transition dipole moments are anisotropic in strength. Furthermore, our observations are consistent with a model in which splitting is driven by shape anisotropy in conjunction with long-range exchange.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Photoinduced Intra- and Intermolecular Energy Transfer in Chlorophyll a Dimer

By applying nonadiabatic excited-state molecular dynamics, we investigate excitation energy transfer and exciton localization dynamics in a chlorophyll a (Chl a ) dimer system at the interface of two monomers of light-harvesting complex II trimer. After its optical excitation at the red edge of the Soret (B) band, the Chl a dimer experiences an ultrafast intra- and intermolecular nonradiative relaxation process to the lowest band (Q y ). The energy relaxation is found to run faster in the Chl a dimer than in the Chl a monomer. Once the molecular system reaches the lowest Q y band composed of two lowest excited states S 1 and S 2 , the concluding relaxation step involves the S 2 → S 1 population transfer, resulting in a relatively slower relaxation rate. The strength of thermal fluctuations exceeds intraband electronic coupling between the states belonging to a certain band (B, Q x , and Q y ), producing localized states on individual chromophores. Therefore, time evolution of spatial electronic localization during internal conversion reveals transient trapping on one of the Chl a monomers participating in the events of intermonomeric energy exchange. In the phase space domains where electronic states are strongly coupled, these states become nearly degenerate promoting Frenkel-exciton-like delocalization and interchromophore energy transfer. As energy relaxation occurs, redistribution of the transition density on two Chl a monomers leads to nearly equal distribution of the exciton among the molecules. For a single Chl a , our analysis of excitonic dynamics reveals wave function amplitude transfer from nitrogen and outer carbon atoms to inner carbon atoms during non-radiative relaxation.

59 BASIC BIOLOGICAL SCIENCES↗

Optical Absorption by Charge-Transfer Molecules

An analytical formalism is developed for the calculation of absorption of radiation by charge-transfer molecules with electronic density delocalized between the donor and acceptor parts. The theory consistently incorporates both the vibronic coupling to quantum intramolecular vibrations and electrostatic interactions with a classical polarizable medium. The formulation operates in terms of basis-invariant parameters and can be used for calculations based on both the localized diabatic states and delocalized adiabatic wave function produced by standard quantum-chemistry algorithms. The basis-invariant reorganization energy is particularly important, since it determines observable spectroscopic parameters, in contrast to the adiabatic reorganization energy based on the molecular charge distributions in the adiabatic vacuum states. Analytical formulas are derived for the charge-transfer absorption band and first two spectral moments connecting spectroscopy to properties of charge-transfer molecules. Electronic delocalization requires a new definition of the Huang–Rhys factor traditionally used to construct vibronic band shapes. Theory’s performance is illustrated by application to experimental data for self-exchange charge-transfer optical absorption and for metal-to-ligand emission transitions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Enabling Efficient and Accurate Computational Studies of MOF Reactivity via QM/MM and QM/QM Methods

Electronic structure calculations can provide unique insight into metal–organic framework (MOF) reactivity and defect formation. Such calculations can be broadly categorized as utilizing either periodic or cluster models, each with their respective advantages and disadvantages. In the present study, we demonstrate how multiscale methods can leverage the advantages of both approaches to enable high levels of accuracy and computational efficiency in studies of MOF reactivity. Using defect formation in a zeolitic imidazolate framework (ZIF) as a prototypical example, we benchmark a quantum mechanics/molecular mechanics (QM/MM) scheme that enables highly efficient cluster-based calculations on MOFs. We demonstrate the importance of correctly accounting for the influence of both dative bond cleavage in the QM cluster and long-range mechanical coupling to the bulk to achieving accurate QM/MM studies of MOFs. We subsequently leverage these cluster models in a QM/QM scheme that goes beyond standard DFT to yield “gold-standard” correlated wave function results on MOFs at modest computational cost. Crucially, we find several cases in which the incorporation of these correlated corrections yields qualitatively important corrections over conventional DFT values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Entropy-Driven Charge Separation: A Potential Explanation for the Low Energy Loss Found in OPVs with Non-fullerene Acceptors

The recent development of non-fullerene acceptors (NFA) has led to an abrupt increase in the organic photovoltaic (OPV) efficiency from ~12% to ~20%. NFA OPVs can generate photocurrent efficiently without much energy loss, which contrasts with the large energy loss often needed for generating free charges from excitons in fullerene OPVs. Here, we argue that entropy-driven charge separation, which allows excitons to spontaneously gain heat from the environment, is responsible for the exceptionally low energy loss found in NFA OPVs. The entropic driving force is maximized when the delocalized electron and hole wave functions in the CT exciton intersect each other via pointlike junctions. Furthermore, these pointlike junctions are abundant in the polymer/NFA bulk heterojunction which consists of intertwined, stringlike donor and acceptor domains in which electron delocalizes anisotropically. Such microstructure favors entropy-driven charge separation, which can enable charge separation with a minimal energy loss.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomic Oxygen Recombination at Surface Defects on Reconstructed (0001) α-Quartz Exposed to Atomic and Molecular Oxygen

The surface chemistry of silica is strongly affected by the nature of chemically active sites (or defects) occurring on the surface. Here, we employ quantum mechanical electronic structure calculations to study an uncoordinated silicon defect, a non-bridging oxygen defect, and a peroxyl defect on the reconstructed (0001) surface of α-quartz. We characterized the spin states and energies of the defects, and calculated the reaction profiles for atomic oxygen recombination at the defects. We elucidated the diradical character by analyzing the low-lying excited states using multireference wave function methods. We show that the diradical defects consist of weakly coupled doublet radicals, and the atomic oxygen recombination can take place through a barrierless process at defects. We have delineated the recombination mechanism and computed the formation energy of the peroxyl and non-bridging oxygen defects. We found that key recombination reaction paths are barrierless. In addition, we characterize the electronically excited states that may play a role in the chemical and physical processes that occur during recombination on these surface defect sites.

Environmental Molecular Sciences Laboratory↗

Theoretical Approaches to Describing the Oxygen Reduction Reaction Activity of Single-Atom Catalysts

Single-atom catalysts have recently emerged as promising low-cost alternatives to Pt for the oxygen reduction reaction (ORR). Given the unique properties that distinguish these systems from traditional transition-metal electrocatalysts, it is essential to benchmark and establish appropriate computational approaches to study these novel materials. Herein, we employ multiple levels of theory, including wave function methods, density functional theory (DFT), and classical simulations, to investigate Cu-modified covalent triazine framework catalysts (Cu/CTF). We consider three major aspects of treating this system computationally. First, we present a step-wise approach to predict the ORR mechanism and adsorbate coverages on Cu/CTF. We then benchmark various DFT methods to coupled-cluster theory with the domain-based local pair natural orbital approximation, which indicates that HSE06 and PBE0 hybrid functionals most accurately describe the adsorption energies of ORR adsorbates on Cu/CTF. We finally employ thermodynamic integration and other techniques to consider solvation effects, which play significant roles in predicting the energies of reaction intermediates and the overall ORR pathway. Lastly, our findings indicate that accurate descriptions of both the electronic structure and solvation are necessary to understand the ORR activity of Cu/CTF.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Decoherence in Molecular Electron Spin Qubits: Insights from Quantum Many-Body Simulations

Quantum states are described by wave functions whose phases cannot be directly measured, but which play a vital role in quantum effects such as interference and entanglement. The loss of the relative phase information, termed decoherence, arises from the interactions between a quantum system and its environment. Decoherence is perhaps the biggest obstacle on the path to reliable quantum computing. Here we show that decoherence occurs even in an isolated molecule although not all phase information is lost via a theoretical study of a central electron spin qubit interacting with nearby nuclear spins in prototypical magnetic molecules. The residual coherence, which is molecule-dependent, provides a microscopic rationalization for the nuclear spin diffusion barrier proposed to explain experiments. The contribution of nearby molecules to the decoherence has a non-trivial dependence on separation, peaking at intermediate distances. Molecules that are far away only affect the long-time behavior. Because the residual coherence is simple to calculate and correlates well with the coherence time, it can be used as a descriptor for coherence in magnetic molecules. This work will help establish design principles for enhancing coherence in molecular spin qubits and serve to motivate further theoretical work

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Excited State Intramolecular Proton Transfer with Nuclear-Electronic Orbital Ehrenfest Dynamics

The recent development of the Ehrenfest dynamics approach in the nuclear-electronic orbital (NEO) framework provides a promising way to simulate coupled nuclear-electronic dynamics. Our previous study showed that the NEO-Ehrenfest approach with a semi- classical traveling proton basis method yields accurate predictions of molecular vibrational frequencies. In this work, we provide a more thorough analysis of the semiclassical traveling proton basis method to elucidate its validity and convergence behavior. We also conduct NEO-Ehrenfest dynamics simulations to study an excited state intramolecular proton transfer process. These simulations reveal that nuclear quantum effects influence the predictions of proton transfer reaction rates and kinetic isotope effects due to the intrinsic delocalized nature of the quantum nuclear wave function. This work illustrates the importance of nuclear quantum effects in coupled nuclear-electronic dynamical processes and shows that the NEO-Ehrenfest approach can be a powerful tool to provide insights and predictions for these processes.

Zhao, Luning↗

Internal Conversion between Bright (1 1 B $^+_u$) and Dark (21 A $^-_g$) States in s- trans -Butadiene and s- trans -Hexatriene

Internal conversion (IC) between the two lowest singlet excited states, 1 1 B$_u^+$ and 2 1 A$_g^–$, of s-trans-butadiene and s-trans-hexatriene is investigated using a series of single- and multi- reference wave function and density functional theory (DFT) methodologies. Three independent types of the equation-of-motion coupled-cluster (EOMCC) theory capable of providing an accurate and balanced description of one- as well as two-electron transitions, abbreviated as δ-CR-EOMCC(2,3), DIP-EOMCC(4h2p){N o }, and DEA-EOMCC(4p2h){N u } or DEA-EOMCC(3p1h,4p2h){N u }, consistently predict that the 1 1 B$_u^+$/2 1 A$_g^–$ crossing in both molecules occurs along the bond length alternation coordinate. However, the analogous 1 1 B$_u^+$ and 2 1 A$_g^–$ potentials obtained with some multireference approaches, such as CASSCF and MRCIS(D), as well as with the linear-response formulation of time-dependent DFT (TDDFT), do not cross. Hence, caution needs to be exercised when studying the low-lying singlet excited states of polyenes with conventional multiconfigurational methods and TDDFT. The multistate many-body perturbation theory methods, such as XMCQDPT 2 , do correctly reproduce the curve crossing. Among the simplest and least expensive computational methodologies, the DFT approaches that incorporate the contributions of doubly excited configurations, abbreviated as MRSF (mixed reference spin-flip) TDDFT and SSR( 4,4 ), accurately reproduce our best EOMCC results. This is highly promising for nonadiabatic molecular dynamics simulations in larger systems.

36 MATERIALS SCIENCE↗