Search NASA⌕ Search

SEARCH · Search NASA

Results for “Chemistry, Physical”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

How Accurate Are Approximate Density Functionals for Noncovalent Interaction of Very Large Molecular Systems?

Noncovalent intermolecular interactions are very important in many research areas. Therefore, it is vital to understand the extent to which approximate density functionals give a proper description of noncovalent interactions. Previous research has demonstrated that some approximate density functionals can predict usefully accurate interaction energies for many noncovalent systems; however, most of that work is limited to small and moderate-sized molecules. Very recently though, accurate benchmarks have become available for some very large molecules. Here, the present work applies 21 approximate density functionals to compute the binding energies of seven large molecular systems that have a number of atoms ranging from 200 to 910. The results are judged by comparison to the recently published CIM-DLPNO-CCSD(T) results, which are assumed to provide a reliable benchmark. The five most accurate methods among those tested are found to be PW6B95-D4, PW6B95-D3(BJ), revM11, M06-L, and MN15.

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

Examination of How Well Long-Range-Corrected Density Functionals Satisfy the Ionization Energy Theorem

For this work, we calculated the vertical ionization energies (VIE) of 99 species in two ways to examine the accuracy of several long-range-corrected (LC) hybrid meta functionals in comparison with a gradient approximation (GA), global hybrids, and doubly hybrids. In the category of LC functionals, we examined both those with meta ingredients (i.e., that depend on the kinetic energy density) and those without them. The LC-hybrid meta functionals examined are M11, revM11, M11plus, and ωB97M-V. The reference data used to assess accuracy consist of 95 molecules and 4 atoms in the GW100 set. The two methods studied are the ΔSCF method (involving the difference of neutral and cation self-consistent field (SCF) energies) and the ionization energy theorem (involving the orbital energy of the highest occupied molecular orbital, HOMO). We calculated linear correlation coefficients (r 2 ) and mean absolute deviations (MADs) between each approach and the reference VIE value from the CCSD(T)/def2-TZVPP level of theory. We compared the new LC-hybrid meta calculations to calculations with the 10 functionals in a previous VIE study by Brémond et al. and to the calculations with LC-BLYP (LC-Becke, Lee–Yang–Parr), CAM-B3LYP (Coulomb-attenuating-method Becke-3-parameter Lee–Yang–Parr), LC-ωHPBE, and ωB97X-D. The results show that Minnesota LC-hybrid meta functionals have the smallest mean absolute deviation of ionization energy theorem VIEs with the reference data; the LC-ωHPBE functional also does quite well in this test. This is very encouraging and indicates that LC-hybrid meta functionals would be the best starting points for the tuning strategy that has been shown to be a very good procedure for improving time-dependent density functional calculations, and it also helps explain the good success of LC-hybrid meta functionals for molecular excitation energies.

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

Theoretical Prediction and Interpretation of 237 Np Mössbauer Isomer Shifts

Different theoretical approaches for the calculation of 237 Np Mössbauer isomer shifts are investigated. The traditional contact density route is compared to a previously proposed alternative approach that uses energy derivatives with respect to the nuclear radius. Both approaches yield similar results as long as suitable basis sets augmented with large exponents and relativistic methods are used. Density functional theory (DFT) calculations do not show a strong dependency of the 237 Np isomer shift on the chosen functional. Wavefunction calculations show that dynamic electron correlation can be important when covalent bonding influences the isomer shift. Effects from spin–orbit coupling are small. The isomer shifts of ionic solids and Np(III) organometallic complexes are largely governed by the oxidation state of Np. Isomer shifts of organometallic Np(IV) complexes are strongly affected by donation bonding. Furthermore, detailed analysis of the wavefunction results with different active spaces demonstrates that correlation among the outer core Np and occupied ligand frontier orbitals contributes significantly to isomer shifts of Np(IV) compounds.

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

Methods to Calculate Electronic Excited-State Dynamics for Molecules on Large Metal Clusters with Many States: Ensuring Fast Overlap Calculations and a Robust Choice of Phase

Here, we present an efficient set of methods for propagating excited-state dynamics involving a large number of configuration interaction singles (CIS) or Tamm-Dancoff approximation (TDA) single-reference excited states. Specifically, (i) following Head-Gordon et al., we implement an exact evaluation of the overlap of singly-excited CIS/TDA electronic states at different nuclear geometries using a biorthogonal basis and (ii) we employ a unified protocol for choosing the correct phase for each adiabat at each geometry. For many-electron systems, the combination of these techniques significantly reduces the computational cost of integrating the electronic Schrodinger equation and imposes minimal overhead on top of the underlying electronic structure calculation. As a demonstration, we calculate the electronic excited-state dynamics for a hydrogen molecule scattering off a silver metal cluster, focusing on high-lying excited states, where many electrons can be excited collectively and crossings are plentiful. Interestingly, we find that the high-lying, plasmon-like collective excitation spectrum changes with nuclear dynamics, highlighting the need to simulate non-adiabatic nuclear dynamics and plasmonic excitations simultaneously. In the future, the combination of methods presented here should help theorists build a mechanistic understanding of plasmon-assisted charge transfer and excitation energy relaxation processes near a nanoparticle or metal surface.

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

Pseudodiagonalization Method for Accelerating Nonlinear Subspace Diagonalization in Density Functional Theory

In density functional theory, each self-consistent field (SCF) nonlinear step updates the discretized Kohn-Sham orbitals by solving a linear eigenvalue problem. The concept of pseudodiagonalization is to solve this linear eigenvalue problem approximately, and specifically utilizing a method involving a small number of Jacobi rotations that takes advantage of the good initial guess to the solution given by the approximation to the orbitals from the previous SCF iteration. The approximate solution to the linear eigenvalue problem can be very rapid, particularly for those steps near SCF convergence. Here, we adapt pseudodiagonalization to finite-temperature and metallic systems, where partially-occupied orbitals must be individually resolved with some accuracy. We apply pseudodiagonalization to the subspace eigenvalue problem that arises in Chebyshev-filtered subspace iteration. In tests on metallic and other systems for a range of temperatures, we show that pseudodiagonalization achieves similar rates of SCF convergence to exact diagonalization.

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

237 Np Mössbauer Isomer Shifts: A Lesson About the Balance of Static and Dynamic Electron Correlation in Heavy Element Complexes

A large set of neptunium compounds with different oxidation states (III to VII) was assembled to study the Mössbauer isomer shift by wave function calculations and better understand covalency in f-elements complexes. The contact density approach was used to calculate the isomer shift using complete active space self-consistent field (CASSCF) multiconfiguration wave functions, as well as matrix product states [from Density Matrix Renormalization Group (DMRG) algorithms] for large active spaces. Dynamic correlation effects for the isomer shifts were treated via CASPT2 energy derivatives with respect to the nuclear radius. The CASSCF calculations appear to produce different orbital overlocalization errors for low and high Np oxidation states. For compounds with low Np oxidation numbers, the errors can be attributed to the overlocalization of the 5f orbitals. For the compounds with high Np oxidation numbers, the main errors arise from an overlocalization of ligand orbitals and concomitant to weak donation bonding. Attempts to mitigate the overlocalization errors with large active spaces using DMRG were only partially successful, showing that explicit treatment of dynamic correlation is necessary for accurate predictions of Mössbauer isomer shifts. The CASPT2 calculations perform very satisfactorily. For a subset of Np compounds, both static and dynamic correlation effects were substantial. Furthermore, a rational active space selection based on orbital entanglement diagrams proved beneficial for determining the optimal reference wave function.

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

Nonadiabatic Dynamics in a Continuous Circularly Polarized Laser Field with Floquet Phase-Space Surface Hopping

Non-adiabatic chemical reactions involving continuous circularly polarized light (cw CPL) have not attracted as much attention as dynamics in unpolarized/linearly polarized light. However, including circularly (in contrast to linearly) polarized light allows one to effectively introduce a complex-valued time-dependent Hamiltonian, which offers a new path for control or exploration through the introduction of Berry forces. Here, we investigate several inexpensive semiclassical approaches for modeling such nonadiabatic dynamics in the presence of a time-dependent complex-valued Hamiltonian, beginning with a straightforward instantaneous adiabatic fewest-switches surface hopping (IA-FSSH) approach (where the electronic states depend on position and time), continuing to a standard Floquet fewest switches surface hopping (F-FSSH) approach (where the electronic states depend on position and frequency), and ending with an exotic Floquet phase-space surface hopping (F-PSSH) approach (wherethe electronic states depend on position, frequency, and momentum). Using a set of model systems with time-dependent complex-valued Hamiltonians, we show that the Floquet phase- space adiabats are the optimal choice of basis as far as accounting for Berry phase effects and delivering accuracy. Furthermore, the F-PSSH algorithm sets the stage for modeling nonadiabatic dynamics under strong externally pumped circular polarization in the future.

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

A Practical Approach to Wave Function Propagation, Hopping Probabilities, and Time Steps in Surface Hopping Calculations

We compare several established approaches for propagating wave functions and calculating hopping probabilities within the fewest switches surface hopping (FSSH) algorithm for difficult cases with many electronic states and many trivial crossings. If only a single time step (Δt c ) is employed, we find that no published approach can accurately capture the dynamics correctly unless Δt c → 0 (which is not computationally feasible). If multiple time steps are employed, for a fixed classical time step (Δt c ), a robust scheme can be found for dynamically choosing quantum time steps (δt q1 and δt q2 ) and calculating hopping probabilities so that one can systematically reduce all errors and achieve maximally efficient accuracy; scattering calculations confirm that one can choose a fairly large classical time step. Furthermore, the robust scheme presented here uses both the “local diabatic” and adiabatic interpolation and thus borrows elements from both the Granucci/Persico and Meek/Levine algorithms. Our findings should be broadly applicable in the future.

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

Nuclear–Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear–electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. Furthermore, these developments and initial examples provide the foundation for future studies of nuclear–electronic quantum dynamics in complex chemical and biological environments.

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

Anion Resonances and Photoelectron Spectroscopy of the Tetracenyl Anion

The photoelectron spectroscopy of the tetracenyl anion using slow electron velocity-map imaging (SEVI) of cryogenically cooled ions is presented here. The total photodetachment yield as a function of photon energy is used to reveal a rich manifold of anion excited states above the detachment threshold. The lowest energy anionic resonance has a sufficiently long lifetime to yield a vibrationally resolved absorption spectrum that can be directly compared with theoretical predictions. Excitation of this state mostly results in electron detachment via thermionic emission. The total photodetachment yield spectrum is used to select photon wavelengths that minimize the indirect detachment signal to allow acquisition of vibrationally-resolved photoelectron spectra that can inform on the neutral tetracenyl radical. Assignment of spectral features corresponding to the ground and first excited state of the neutral 12-tetracenyl isomer is made with the aid of Franck-Condon simulations. This yields adiabatic electron affinity and term energies that differ significantly from the previously reported values. Weak features corresponding to the ground state of the minor 2-teracenyl and 1-tetracenyl isomers are also identified, which allows for the experimental determination of their electron affinities for the first time.

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

Photodissociation Dynamics of Astrophysically Relevant Propyl Derivatives (C 3 H 7 X; X = CN, OH, HCO) at 157 nm Exploiting an Ultracompact Velocity Map Imaging Spectrometer: The (Iso)Propyl Channel

The photodissociation dynamics of astrophysically relevant propyl derivatives (C 3 H 7 X; X = CN, OH, HCO) at 157 nm exploiting an ultracompact velocity map imaging (UVMIS) setup has been reported. The successful operation of UVMIS allowed the exploration of the 157 nm photo dissociation of six (iso)propyl systems - n/i-propyl cyanide (C 3 H 7 CN), n/i-propyl alcohol (C 3 H 7 OH), and (iso)butanal (C 3 H 7 CHO) – to explore the C 3 H 7 loss channel. The distinct center-of-mass translational energy distributions for the i-C 3 H 7 X (X= CN, OH, HCO) could be explained through preferential excitation of the low frequency C-H bending modes of the formyl moiety compared to the higher frequency stretchings of the cyano and hydroxy moieties. Although the ionization energy of the n-C 3 H 7 radical exceeds the energy of a 157 nm photon, C 3 H 7 + was observed in the n-C 3 H 7 X (X= CN, OH, HCO) systems as a result of photoionization of vibrationally "hot" n-C 3 H 7 fragments, photoionization of i-C 3 H 7 after a hydrogen shift in vibrationally "hot" n-C 3 H 7 radicals, and/or two-photon ionization. Our experiments reveal that at least the isopropyl radical (i-C 3 H 7 ) and possibly the normal propyl radical (n-C 3 H 7 ) should be present in the interstellar medium and hence searched for by radio telescopes.

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

Competition between Solvation and Intramolecular Hydrogen-Bonding in Microsolvated Protonated Glycine and β-Alanine

Infrared predissociation (IRPD) spectroscopy is used to reveal and compare the microsolvation motifs of GlyH + (H 2 O) n and β-AlaH + (H 2 O) n . The chemical structure of these amino acids differ only in the length of the carbon chain connecting the amine and carboxyl terminals, which nonetheless leads to a significant difference in the strength of the intramolecular C=O···H-N hydrogen bond in the unsolvated ions. This difference makes them useful in our studies of the competition between solvation and internal hydrogen bonding interactions. Analysis of the IRPD results reveals that the sequential addition of water molecules leads to similar effects on the intramolecular interaction in both GlyH + (H 2 O) n and β-AlaH + (H 2 O) n . Solvation of the -NH 3 + group leads to a weakening of the C=O···H-N hydrogen bond, while solvation of the carboxyl -OH leads to a strengthening of this bond. Additionally, we have found that for β-AlaH + , the addition of a H 2 O to the second solvation shell can still influence the strength of the C=O···H-N hydrogen bonding interaction. Lastly because the C=O···H-N interaction in β-AlaH + is stronger than that in GlyH + , more solvent molecules are needed to sufficiently weaken the intramolecular hydrogen bond such that isomers without this bond begin to be energetically competitive; this occurs at n=5 for β-AlaH + and n=1 for GlyH + .

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

Accelerated Self-Healing and Property Recovery in Brush Particle Solids Featuring Brush Dispersity

Brush particles, hybrid materials consisting of polymer chains tethered to particle surfaces, offer tunable properties that make them promising candidates for advanced functional materials. This study investigated the role of chain dispersity in the viscoelastic self-healing of poly (methyl acrylate) (PMA)-based brush particle solids. Increasing the molecular weight dispersity of grafted chains significantly enhanced both strain-to-fracture and toughness of brush particle solids, while the elastic modulus and glass transition temperature were independent of chain dispersity. Cut-and-adhere testing revealed a significant acceleration of the rate of toughness recovery in high-dispersity systems as compared to low-dispersity analogs for which toughness recovery markedly lagged the recovery of Young’s modulus. The results suggest that structure and property recovery in brush particle solids are sensitive to the dynamical heterogeneity of brush canopies and highlight the role of molecular weight dispersity as a design parameter to enable hybrid materials with advanced self-healing ability.

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

Joint Experimental and Computational Characterization of Sum-Frequency Generation between a Continuous Wave Laser and an Ultrafast Frequency Comb Laser for Tunable Laser Development

Ultrafast optical frequency combs allow for both high spectral and temporal resolution in molecular spectroscopy and have become a powerful tool in many areas of chemistry and physics. Ultrafast lasers and frequency combs generated from ultrafast mode-locked lasers often need to be converted to other wavelengths. Commonly used wavelength conversions are optical parametric oscillators, which require an external optical cavity, and supercontinuum generation combined with optical parametric amplifiers. Whether commercial or home-built, these systems are complex and costly. Here, we investigate an alternative, simple, and easy-to-implement approach to tunable frequency comb ultrafast lasers enabled by new continuous-wave laser technology. Sum-frequency generation between an Nd:YAG continuous-wave laser and a Yb:fiber femtosecond frequency comb in a beta-barium borate (BBO) crystal is explored. The resulting sum-frequency beam is a pulsed frequency comb with the same repetition rate as the Yb:fiber source. SNLO simulation software is used to simulate the results and provide benchmarks for designing future systems to achieve wavelength conversion and tunability in otherwise difficult-to-reach spectral regions.

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

Real-time degradation dynamics of hydrated per- and polyfluoroalkyl substances (PFASs) in the presence of excess electrons

Per- and polyfluoroalkyl substances (PFASs) are synthetic chemicals that are harmful to both the environment and human health. In this work, using self-interaction-corrected Born–Oppenheimer molecular dynamics simulations, we provide the first real-time assessment of PFAS degradation in the presence of excess electrons. In particular, we show that the initial phase of the degradation involves the transformation of an alkane-type C–C bond into an alkene-type C[double bond, length as m-dash]C bond in the PFAS molecule, which is initiated by the trans elimination of fluorine atoms bonded to these adjacent carbon atoms.

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

Formation of resonances and anionic fragments upon electron attachment to benzaldehyde

Benzaldehyde is a simple aromatic aldehyde and has a wide range of applications in the food, pharmaceutical, and chemical industries. The positive electron affinity of this compound suggests that low-energy electrons can be easily trapped by neutral benzaldehyde. In the present study, we investigated the formation of negative ions following electron attachment to benzaldehyde in the gas-phase. Calculations on elastic electron scattering from benzaldehyde indicate a π* valence bound state of the anion at -0.48 eV and three π* shape resonances (0.78, 2.48 and 5.51 eV). The excited state spectrum of the neutral benzaldehyde is also reported to complement our findings. Using mass spectrometry, we observed the formation of the intact anionic benzaldehyde at ~0 eV. We ascribe the detection of the benzaldehyde anion to stabilization of the π* valence bound state upon dissociative electron attachment to a benzaldehyde dimer. In addition, we report the cross sections for nine fragment anions formed through electron attachment to benzaldehyde. Investigations carried out with partially deuterated benzaldehyde show that the hydrogen loss is site-selective with respect to the incident electron energy. In addition, we propose several dissociation pathways, backed up by quantum chemical calculations on their thermodynamic thresholds. The threshold calculations also support that the resonances formed at higher energies lead to fragment anions observable by mass spectrometry, whereas the resonances at low electron energies decay only by electron autodetachment.

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

The water trimer reaction OH + (H 2 O) 3 → (H 2 O) 2 OH + H 2 O

All important stationary points on the potential energy surface (PES) for the reaction OH + (H 2 O) 3 → (H 2 O) 2 OH + H 2 O have been fully optimized using the “gold standard” CCSD(T) method with the large Dunning correlation-consistent cc-pVQZ basis sets. Three types of pathways were found. For the pathway without hydrogen abstraction, the barrier height of the transition state (TS1) is predicted to lie 5.9 kcal mol -1 below the reactants. The two major complexes (H 2 O) 3 ∙∙∙OH (CP1 and CP2a) are found to lie 6.3 and 11.0 kcal mol -1 , respectively, below the reactants [OH + (H 2 O) 3 ]. For one of the H-abstraction pathways the lowest classical barrier height is predicted to be much higher, 6.1 kcal mol -1 (TS2a) above the reactants. For the other H-abstraction pathway the barrier height is even higher, 15.0 (TS3) kcal mol -1 . Vibrational frequencies and the zero-point vibrational energies connected to the PES are also reported. Fianlly, the energy barriers for the H-abstraction pathways are compared with those for the OH + (H 2 O) 2 and OH + H 2 O reactions, and the effects of the third water on the energetics are usually minor (0.2 kcal mol -1 ).

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

High-temperature 2D ferromagnetism in conjugated microporous porphyrin-type polymers

The need for magnetic 2D materials that are stable to the enviroment and have high Curie temperatures is very important for various electronic and spintronic applications. We have found that two-dimensional porphyrin-type aza-conjugated microporous polymer crystals are such a material (Fe-aza-CMPs). Fe-aza-CMPs are stable to CO, CO 2 , and O 2 atmospheres and show unusual adsorption, electronic, and magnetic properties. Indeed, they are semiconductors with small energy band gaps ranging from 0.27 eV to 0.626 eV. CO, CO 2 , and O 2 molecules can be attached in three different ways where single, double, or triple molecules are bound to iron atoms in Fe-aza-CMPs. For different attachment configurations we find that for CO and CO 2 a uniform distribution of the molecules is most energetically favorable while for O 2 molecules aggregation is most energetically preferable. The magnetic moments decrease from 4 to 2 to 0 for singly, doubly, triply occupied configurations for all gasses respectively. The most interesting magnetic properties are found for O 2 molecules attached to the Fe-aza-CMP. For a single attachment configuration we find that an antiferromagnetic state is favorable. When two O 2 molecules are attached, the calculations show the highest exchange integral with a value of J = 1071 μeV. This value has been verified by two independent methods where in the first method J is calculated by the energy difference between ferromagnetic and anitferromagnetic configurations. The second method is based on the frozen magnon approach where the magnon dispersion curve has been fitted by the Ising model. Finally, for the second method J has been estimated at J = 1100 μeV in excellent agreement with the first method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗