Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio methods”

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,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

NbV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys niobium-vanadium (Nb-V). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Nb and V. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. 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 4. 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 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit 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. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y 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: FILES contained in each subdirectory with name case-N where N ranges between 11 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 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 5. 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 6. 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 7. 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 8. 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 9. OUTCAR: 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 10. 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 11. 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 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. Subdirectories with name case-N, where N ranges between 1 and 10 (extremes included) contain all the files listed above except 'vaspout.h5'. Subdirectories with name case-N, where N ranges between 41 and 60 (extremes included), contain a duplicate copy of the files listed above except for KPOINTS. The names of the duplicate files end with -bis, and correspond to a second VASP calculation that has converged to a different optimized geometry. 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. REFERENCES (1) Kresse, G. and Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. and Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal-amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. and Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15-50 (1996) (4) Kresse, G. and Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. and Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
NASA NTRS1989

Ionization and thermal equilibrium models for O star winds based on time-independent radiation-driven wind theory

Ab initio ionization and thermal equilibrium models are calculated for the winds of O stars using the results of steady state radiation-driven wind theory to determine the input parameters. Self-consistent methods are used for the roles of H, He, and the most abundant heavy elements in both the statistical and the thermal equilibrium. The model grid was chosen to encompass all O spectral subtypes and the full range of luminosity classes. Results of earlier modeling of O star winds by Klein and Castor (1978) are reproduced and used to motivate improvements in the treatment of the hydrogen equilibrium. The wind temperature profile is revealed to be sensitive to gross changes in the heavy element abundances, but insensitive to other factors considered such as the mass-loss rate and velocity law. The reduced wind temperatures obtained in observing the luminosity dependence of the Si IV lambda 1397 wind absorption profile are shown to eliminate any prospect of explaining the observed O VI lambda 1036 line profiles in terms of time-independent radiation-driven wind theory.

Drew, J. E.↗
DOE OSTI2011

Origin, Evolution And Imaging Of Vortices In Atomic Processes

Calculations of the time-dependent wave function for proton impact on atomic hydrogen using the Lattice-Time-Dependent-Schroedinger-Equation (LTDSE) method follow the wave function from microscopic to macroscopic times. Isolated zeros, now identified as vortices, appear when the target and projectile nuclei are separated by a few atomic units. Such structure has apparently been observed frequently in ab. initio. calculations. Our work shows that such structures persist to macroscopic distances and appear as 'holes' in electron momentum distributions. They can, in principle, be observed in reaction microscope studies. Such observation is formally justified by the 'imaging theorem', which can be derived from first principles. Similarly, two-electron momentum distributions as observed, e.g., in (e,2e) coincidence measurements may show isolated zeros. Calculations using correlated wave functions support the interpretation that a minima in (e,2e) for helium targets corresponds to a vortex. In this case the 'imaging theorem' can be applied to argue that there must be vortices in the two-electron wave function. A general discussion of vortices in quantum mechanics will be illustrated using exact LTDSE calculations, and, for interpretive purposes, the time-dependent first Born approximation. These calculations show that the plane wave B1 amplitudes have no vortices, but the time-dependent B1 amplitudes do. It will be further shown that angular momentum transfer is the key to forming vortices in the time-dependent theory.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Deuteron-induced nucleon transfer reactions within an ab initio framework: First application to p -shell nuclei

Low-energy transfer reactions in which a proton is stripped from a deuteron projectile and dropped into a target play a crucial role in the formation of nuclei in both primordial and stellar nucleosynthesis, as well as in the study of exotic nuclei using radioactive beam facilities and inverse kinematics. Here, ab initio approaches have been successfully applied to describe the 3 H(d,n) 4 He and 3 He(d,p) 4 He fusion processes. An ab initio treatment of transfer reactions would also be desirable for heavier targets. In this work, we extend the ab initio description of (d,p) reactions to processes with light p-shell nuclei. As a first application, we study the elastic scattering of deuterium on 7 Li and the 7 Li(d,p) 8 Li transfer reaction based on a two-body Hamiltonian. We use the no-core shell model to compute the wave functions of the nuclei involved in the reaction, and describe the dynamics between targets and projectiles with the help of microscopic-cluster states in the spirit of the resonating group method. The shapes of the excitation functions for deuterons impinging on 7 Li are qualitatively reproduced up to the deuteron breakup energy. The interplay between d– 7 Li and p– 8 Li particle-decay channels determines some features of the 9 Be spectrum above the d+ 7 Li threshold. Our prediction for the parity of the 17.298 MeV resonance is at odds with the experimental assignment. Deuteron stripping reactions with p-shell targets can now be computed ab initio , but calculations are very demanding. Finally, a quantitative description of the 7 Li(d,p) 8 Li reaction will require further work to include the effect of three-nucleon forces and additional decay channels and to improve the convergence rate of our calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2012

Atomistic simulation and ab initio study of the defect structure of spinel-related Li{sub 0.5−0.5x}Mg{sub x}Fe{sub 2.5−0.5x}O{sub 4}

Graphical abstract: Unit cell of Li0{sub 5−0.5x}Mg{sub x}Fe{sub 2.5−0.5x}O{sub 4}, showing the lowest energy structure obtained using interatomic potential and DFT ab initio calculations. Large white spheres O{sup 2−}; small light grey spheres Mg{sup 2+} (evenly substituting of Li{sup +} and Fe{sup 3+} at octahedral sites); small dark grey spheres Fe{sup 3+}; small black spheres Li{sup +}. Display Omitted Highlights: ► Defect structure of Li{sub 0.5−0.5x}Mg{sub x}Fe{sub 2.5−x}O{sub 4} is studied with atomistic and DFT methods. ► 19 possible defect structure models with ∼60 defect configurations are investigated. ► The most favourable model found is when Mg{sup 2+} ions evenly replace Li{sup +} and octahedral Fe{sup 3+}. ► This defect structure decreases the magnetisation relative to that of Li{sub 0.5}Fe{sub 2.5}O{sub 4}. ► Experimentally-deduced models, at variance with the one obtained here, are discussed. -- Abstract: The position of magnesium ions in Mg{sup 2+}-doped lithium ferrite of the composition Li{sub 0.5−0.5x}Mg{sub x}Fe{sub 2.5−0.5x}O{sub 4}, which has been a matter of uncertainty among some experimentalists, is investigated using interatomic potential and ab initio DFT calculations. Among possible 19 defect structure models, some of which have been reported experimentally to be the most favorable, the lowest energy is found for Mg{sup 2+} ions evenly replacing Li{sup +} and Fe{sup 3+} ion on octahedral sites. This gives a decrease in magnetisation for the Mg{sup 2+}-doped ferrite relative to the un-doped lithium ferrite. The results suggest that some experimental observations of increased magnetisation of spinel lithium ferrite on Mg{sup 2+}-doping could be due to substitution of Mg{sup 2+} or Li{sup +} on tetrahedral sites at the high temperatures used in preparation of the solid and/or the presence of undetected defects in the initial precursors.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Theoretical level energies and transition data for 4p{sup 6}4d{sup 5}, 4p{sup 5}4d{sup 6} and 4p{sup 6}4d{sup 4}4f configurations of W{sup 33+} ion

The ab initio quasirelativistic approach developed specifically for the calculation of spectral parameters of highly charged ions was used to derive transition data for the Nb-like tungsten ion W{sup 33+}. The configuration interaction method was applied to include electron correlation effects. The relativistic effects were taken into account in the Breit–Pauli approximation. The level energies, radiative lifetimes τ, Landé g-factors were determined for the ground configuration 4p{sup 6}4d{sup 5} and two excited configurations 4p{sup 5}4d{sup 6} and 4p{sup 6}4d{sup 4}4f. The radiative transition wavelengths λ and emission transition probabilities A for the electric dipole, electric quadrupole, electric octupole, magnetic dipole, and magnetic quadrupole transitions among the levels of these configurations were produced. The uncertainties of computed spectroscopic parameters were evaluated.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Computational study of the electrostatic potential and charges of multivalent ionic liquid molecules

The electrostatic interactions in ionic liquids (ILs) can be difficult to quantify, but this information is very valuable for improving intermolecular potentials for performing molecular-level simulations. Although polarizable forcefields can provide improved accuracy, the computational efficiency of fixed charge molecular models is often preferred for larger systems. Here, we propose a multi-scale screening approach for analyzing the charge scaling behavior and other electrostatic properties of ILs, which is more attainable than previous analysis methods that rely on the IL crystalline structures. Based on isolated molecules and small clusters, we study the electrostatic potential properties and charges of multivalent IL molecules using density-functional theory calculations and ab initio molecular dynamics simulations. The charges of the IL molecules are calculated and compared using several different charge estimation techniques. Although there are some differences among the estimated partial charges on the atoms, the total anion and cation charges are very similar. In conclusion, the charge scaling factors are substantial, and they are only slightly affected by the IL cluster size model.

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

Electron–Nucleus Hyperfine Coupling Calculated from Restricted Active Space Wavefunctions and an Exact Two-Component Hamiltonian

Exact two-component (X2C) relativistic nuclear hyperfine magnetic field operators were incorporated in X2C ab-initio wavefunction calculations at the multi-reference restricted active space (RAS) level for calculations of nuclear hyperfine magnetic properties. Spin-orbit coupling was treated via RAS state interaction (SO-RASSI). The method was tested by calculations of electron – nucleus hyperfine coupling constants. The approach, implemented in the OpenMolcas program, overcomes a major limitation of a previous SO-RASSI implementation for hyperfine coupling that relied on non-relativistic hyperfine operators [J. Chem. Theor. Comput. 2015, 11, 538–549] and therefore had only limited applicability. Furthermore, results from calculations on systems with light and heavy main group elements, transition metals, lanthanides, and one actinide complex, demonstrate reasonably good agreement with experimental data, where available, as long as the active space can generate sufficient spin polarization.

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

Stimulated X-ray Raman and Absorption Spectroscopy of Iron–Sulfur Dimers

Iron–sulfur complexes play an important role in biological processes such as metabolic electron transport. A detailed understanding of the mechanism of long-range electron transfer requires knowledge of the electronic structure of the complexes, which has traditionally been challenging to obtain, either by theory or by experiment, but the situation has begun to change with advances in quantum chemical methods and intense free electron laser light sources. We compute the spectra for stimulated X-ray Raman spectroscopy (SXRS) and absorption spectroscopy of homovalent and mixed-valence [2Fe–2S] complexes, using the ab initio density matrix renormalization group algorithm. The simulated spectra show clear signatures of the theoretically predicted dense low-lying excited states within the d–d manifold. Finally, the difference in spectral intensity between the absorption-active and Raman-active states provides a potential mechanism to selectively excite states by a proper tuning of the excitation pump, to access the electronic dynamics within this manifold.

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

The Li + CaF → Ca + LiF chemical reaction under cold conditions

The calcium monofluoride (CaF) molecule has emerged as a promising candidate for precision measurements, quantum simulation, and ultracold chemistry experiments. Inelastic and reactive collisions of laser cooled CaF molecules in optical tweezers have recently been reported and collisions of cold Li atoms with CaF are of current experimental interest. In this paper, we report ab initio electronic structure and full-dimensional quantum dynamical calculations of the Li + CaF → LiF + Ca chemical reaction. The electronic structure calculations are performed using the internally contracted multi-reference configuration-interaction method with Davidson correction (MRCI + Q). An analytic fit of the interaction energies is obtained using a many-body expansion method. Furthermore, a coupled-channel quantum reactive scattering approach implemented in hyperspherical coordinates is adopted for the scattering calculations under cold conditions. Results show that the Li + CaF reaction populates several low-lying vibrational levels and many rotational levels of the product LiF molecule and that the reaction is inefficient in the 1–100 mK regime allowing sympathetic cooling of CaF by collisions with cold Li atoms.

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

Atomistic and Ab initio modeling of CaAl{sub 2}O{sub 4} high-pressure polymorphs under Earth’s mantle conditions

Semi-empirical and ab initio theoretical investigation of crystal structure geometry, interatomic distances, phase densities and elastic properties for some CaAl{sub 2}O{sub 4} phases under pressures up to 200 GPa was performed. Two independent simulation methods predicted the appearance of a still unknown super-dense CaAl{sub 2}O{sub 4} modification. In this structure, the Al coordination polyhedron might be described as distorted one with seven vertices. Ca atoms were situated inside polyhedra with ten vertices and Ca–O distances from 1.96 to 2.49 Å. It became the densest modification under pressures of 170 GPa (density functional theory prediction) or 150 GPa (semi-empirical prediction). Both approaches indicated that this super-dense CaAl{sub 2}O{sub 4} modification with a “stuffed α-PbO{sub 2}” type structure could be a probable candidate for mutual accumulation of Ca and Al in the lower mantle. The existence of this phase can be verified experimentally using high pressure techniques.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Defects in SiC for Quantum Computing

Many novel materials are being actively considered for quantum information science and for realizing high-performance qubit operation at room temperature. It is known that deep defects in wide-band gap semiconductors can have spin states and long coherence times suitable for qubit operation. We theoretically investigate from ab-initio density functional theory (DFT) that the defect states in the hexagonal silicon carbide (4H-SiC) are potential qubit materials. The DFT supercell calculations were performed with the local-orbital and pseudopotential methods including hybrid exchange-correlation functionals. Di-vacancies in SiC supercells yielded defect levels in the gap consisting of closely spaced doublet just above the valence band edge, and higher levels in the band gap. The divacancy with a spin state of 1 is charge neutral. Furthermore, the divacancy is characterized by C-dangling bonds and a Si-dangling bonds. Jahn-teller distortions and formation energies as a function of the Fermi level and single photon interactions with these defect levels will be discussed. In contrast, the anti-site defects where C, Si are interchanged have high formation energies of 5.4 eV and have just a single shallow defect level close to the valence band edge, with no spin. We will compare results including the defect levels from both the electronic structure approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Investigation of the pressure dependent thermodynamic and elastic properties of 1,3,5-triamino-2,4,6-trinitrobenzene using dispersion corrected density functional theory

The influence of pressure on the thermodynamic, structural, and elastic properties of the 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) molecular crystal at T = 0 is systematically studied. Calculations are carried out using density functional theory methods in a plane wave basis set with dispersion corrections for the exchange-correlation part of total energy, and ultrasoft pseudopotentials. The equilibrium unit cell parameters, the cold compression curve in the pressure range of 0–50 GPa and the sound speeds are computed. The effect of finite pressure on the molecular structure of TATB is elucidated from the analysis of relative changes in the intra- and intermolecular geometrical parameters. For the first time, the full set of elastic constants of this crystal at zero and non-zero pressures is determined from ab initio calculations. The resulted structural, elastic, and acoustic properties of TATB are shown to be in a good agreement with available experimental and theoretical data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Optimization of an exchange-correlation density functional for water

We describe a method, that we call data projection onto parameter space (DPPS), to optimize an energy functional of the electron density, so that it reproduces a dataset of experimental magnitudes. Our scheme, based on Bayes theorem, constrains the optimized functional not to depart unphysically from existing ab initio functionals. The resulting functional maximizes the probability of being the “correct” parameterization of a given functional form, in the sense of Bayes theory. The application of DPPS to water sheds new light on why density functional theory has performed rather poorly for liquid water, on what improvements are needed, and on the intrinsic limitations of the generalized gradient approximation to electron exchange and correlation. Finally, we present tests of our water-optimized functional, that we call vdW-DF-w, showing that it performs very well for a variety of condensed water systems.

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

High Resolution Infrared Spectroscopy of Highly Reactive Chemical Intermediates: Berkeley Inspiration and a C.B. Moore Retrospective

One of the long-standing paradoxes in chemistry is that the molecular species of greatest interest are often highly reactive chemical intermediates and thus present only at vanishingly small concentrations under typical steady state conditions. Here, this has proven both a key frustration and powerful motivation for physical chemists in probing such species, particularly for high resolution laser spectroscopic methods in the infrared. Over the past 4-5 decades, this fundamental sensitivity challenge has stimulated the development of many novel techniques for efficient generation, cooling, and infrared probing of such transient chemical intermediates, achieving increasingly sophisticated levels of spectroscopic detail, chemical insight, and a testbed for comparison with first principles ab initio quantum calculations. Two of the pioneering contributors to this scientific arena have been George Pimentel and C. Bradley Moore, both responsible for nourishing DJN’s own enduring fascination with spectroscopic IR study of elusive chemical intermediates. This chapter is not intended to provide a comprehensive review of this large and enormously successful field, but it rather more simply attempts to capture a few selective “Moore-centric” snapshots of the Nesbitt group scientific evolution, specifically in i) advances by Pimentel for the original development of rapid scan flash kinetic spectroscopy based on spectrometers with rapidly rotating IR gratings, ii) the modernization/extension of these methods in the Moore group toward the first high resolution IR laser spectroscopy of singlet/triplet methylene 1,3 CH 2 diradical, and finally to iii) the development of slit supersonic discharge expansion methods at JILA for sub-Doppler infrared laser study of multiple radicals, jet cooled molecular ions, and highly reactive chemical intermediates.

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

Simulations of the THz spectrum of liquid water incorporating the effects of intermolecular charge fluxes through hydrogen bonds

The intensity of the band at ∼200 cm{sup −1} (∼6 THz) in the Terahertz spectrum of liquid water mainly arises from the modulations of the extent of intermolecular charge transfer through hydrogen bonds, called intermolecular charge fluxes, occurring upon molecular translations along the O…H hydrogen bonds. To obtain reasonable spectral profiles from simulations, it is necessary to effectively incorporate the effects of those intermolecular charge fluxes, but apparently it is not possible by simple classical molecular dynamics simulations with fixed atomic partial charges even if they are amended by molecular induced dipoles due to intermolecular electrostatic interactions. The present paper shows how we can do reasonably correct spectral simulations, without resorting to fully ab initio molecular dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Vibrational Signatures of Electronic Properties in Renewable-Energy Catalysis

The objective of this research program is to discover and develop new approaches for ab initio computational simulations of molecular motion. Specifically, this program aims to decipher the connection between the underlying electronic structure and the accordant molecular vibrations in reactive ions and radicals for renewable-energy purposes. Recent experimental progress in ion sources and optical spectroscopies has unearthed considerable new inner-sphere detail for these complexes, but the connection between these spectral signatures and mechanistic information often remains elusive, to the continued frustration of experimentalists. For this purpose, new anharmonic vibrational frequency methods, along with a publicly deployed software package, will be developed. Working closely with committed experimental collaborators, this conceptual and computational framework will be used to explain the results of new spectroscopy experiments, focusing specifically on the inner-shell mechanisms of renewable-energy catalysis. The oxidation half of catalytic water-splitting chemistry will be a central focus, along with fundamental studies of the manner in which strong ions and radicals activate solvent as a chemical species. The resulting products of the research program will include openly available software and algorithms for the ab initio simulation of challenging vibrational spectra, as well as critical mechanistic insight into energy-focused catalytic processes that are opaque to other existing analytical techniques.

42 ENGINEERING↗