Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “correspondence”

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

Anionic Oxygen Redox in the High-Lithium Material Li 8 SnO 6

Lithiated transition metal oxide cathode materials that combine oxygen-anion redox [O 2 – /(O 2 ) n– ] with cation redox offer substantial capacities and higher voltage than cathodes without oxygen redox, but oxygen-anion redox is usually accompanied by the formation of peroxo or superoxo species that may cause oxygen release. To clarify the relationship between the oxygen release and peroxide and superoxide dimers without the complication of transition-metal redox, we performed density functional theory calculations to study the layered Li 8 SnO 6 cathode because it has only oxygen-anion redox during delithiation. Features with O–O distances of ~1.5 Å (corresponding to the formation of a peroxide) were observed at 75% lithium concentration, and further delithiation to 62.5% lithium concentration shortened the O–O bond distance to ~1.3 Å (corresponding to the formation of superoxide). The solely anionic oxygen redox stabilizes the delithiated structure of Li 8–x SnO 6 , and it does not cause structural disorder or release oxygen gas. Ab initio molecular dynamics calculations in the supercell at 300 K indicate that oxygen dimerization may occur at 75% lithium concentration. On the (0001) surface with a high lithium concentration (98%), oxygen dimerization is unfavorable. The calculations by HSE06 show that oxygen redox in Li 8 SnO 6 offers an initial voltage (vs Li/Li + ) greater than 4.2 V. The GGA+U calculations indicate that the Li + -polaron in Li 95 Sn 12 O 72 migrates with a migration barrier of only 0.43 eV. Here, this work also quantifies a solely anionic oxygen redox mechanism in Li 8 SnO 6 , and it provides a deeper understanding of anionic oxygen redox in Li-excess oxide cathode materials.

25 ENERGY STORAGE↗
DOE OSTI2021

Nanoscale Phase Separation and Large Refrigerant Capacity in Magnetocaloric Material LaFe 11.5 Si 1.5

Here, ferromagnetic transitions, the magnetocaloric effect, and the atomic-scale microstructure in an annealed bulk material of LaFe 11.5 Si 1.5 were investigated using magnetic measurements, Mössbauer spectroscopy, high-resolution X-ray diffraction analysis, and high-resolution transmission electron microscopy. The results provide evidence of the coexistence of two ferromagnetic phases with different Curie temperatures. The phase with a Curie temperature of 216 K corresponds to a Si-rich phase of a NaZn 13 -type structure with a small lattice volume, whereas the phase with a Curie temperature of 185 K corresponds to a Si-poor phase of a NaZn 13 -type structure with a large lattice volume. This phase coexistence is observed on a nanometer scale and it can account for inverse thermal hysteresis of the ferromagnetic transition in the Si-rich phase during warming. Furthermore, the phase coexistence helps the bulk material achieve a refrigeration capacity of 170 J kg -1 due to a combination of the magnetocaloric effect in each phase. It is determined that the magnetocaloric effect in the Si-rich and the Si-poor phase brings about a maximum of the total entropy change of -3.55 and -6.80 and J kg -1 K -1 at respective Curie temperatures, respectively. Fundamentally, the phase coexistence indicates coloring of Si atoms in the lattice of a NaZn 13 -type structure and can be attributed to nanoscale phase separation of a disordered precursor during annealing of the bulk material. The nanoscale phase separation may occur in other magnetocaloric materials, thus allowing for improvement of their refrigeration capacity.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Reducing Surface Recombination Velocity of Methylammonium-Free Mixed-Cation Mixed-Halide Perovskites via Surface Passivation

We control surface recombination in the mixed-cation, mixed-halide perovskite, FA0.83Cs0.17Pb(I0.85Br0.15)3, by passivating nonradiative defects with the polymerizable Lewis base (3-aminopropyl)trimethoxysilane (APTMS). We demonstrate average minority carrier lifetimes >4 µs, nearly single exponential monomolecular photoluminescence decays, and high external photoluminescence quantum efficiencies (>20%, corresponding to ~97% of the maximum theoretical quasi-Fermi-level splitting) at low excitation fluence. We confirm both the composition and valence band edge position of the FA0.83Cs0.17Pb(I0.85Br0.15)3 perovskite using multi-institutional, cross-validated, X-ray photoelectron spectroscopy and UV photoelectron spectroscopy measurements. We extend the APTMS surface passivation to higher bandgap double-cation (FA and Cs) compositions (1.7, 1.75, and 1.8 eV) as well as the widely used triple-cation (FA, MA, and Cs) composition. Finally, we demonstrate that the average surface recombination velocity decreases from ~1000 to ~10 cm/s post APTMS passivation for FA0.83Cs0.17Pb(I0.85Br0.15)3. Our results demonstrate that surface-mediated recombination is the primary nonradiative loss pathway in many methylammonium (MA)-free mixed-cation mixed-halide films with a range of different bandgaps, which is a problem observed for a wide range of perovskite active layers and reactive electrical contacts. Our study also provides insights to develop passivating molecules that help reduce surface recombination in MA-free mixed-cation mixed-halide films and indicates that surface passivation and contact engineering will enable near-theoretical device efficiencies with these materials.

14 SOLAR ENERGY↗
DOE OSTI2023

Revealing the Correlation between the Solvation Structures and the Transport Properties of Water-in-Salt Electrolytes

Water-in-salt (WIS) electrolytes containing 21 m lithium bis(trifluoromethane sulfonyl)imide (LiTFSI) have been considered as a safe and environment-friendly alternative to common organic electrolytes used in lithium-ion batteries. However, the relation between the solvation structures and transport properties of these materials remains elusive. Here, for this paper, we performed small-angle X-ray scattering (SAXS), small-angle neutron scattering (SANS), and X-ray pair distribution function (PDF) measurements of LiTFSI aqueous solutions at a wide range of concentrations. Combined with molecular dynamics simulations, the detailed solvation structures from long to short length scale were resolved. We found that the TFSI– solvation structures consist of TFSI– solvated structures and TFSI– networks; the former corresponds to solvent separated ion pairs, while the latter corresponds to contact ion pairs and cation–anion aggregates. In addition, we found that the relaxation time in the q range associated with the anion network structure exhibits the same concentration dependence as the viscosity. By combining the results from the experiments and simulations, this study revealed a correlation between the solvation structures of LiTFSI and the transport properties of the solutions, which is critical to understand the relation between the transport properties and the dynamics of the ions for imide-based lithium-ion salt aqueous electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Spectroscopic Investigation into P(NDI2OD-T2) Charge Localization

P(NDI2OD-T2), commonly referred to as N2200, stands out as a promising electron-transporting (n-type) polymer for low-cost, flexible (photo)electrochemical applications due to its reversible two-electron reduction and high electron mobility. UV–vis spectroelectrochemistry in the tetrabutyl ammonium hexafluorophosphate/acetonitrile electrolyte shows two sets of chemically reversible redox signals in the cyclic voltammetry corresponding to the reduction of the neutral polymer film to polaronic and bipolaronic species. These electrochemical signatures suggest a distinct electronic reorganization upon reduction (polaron/bipolaron formation), highlighting the need for molecular-level insights into how charges are accommodated within the polymer backbone. While it has been previously hypothesized that charge predominantly localizes on the naphthalene diimide (NDI) unit during reductive charging, specific changes in atomic environments that confirm this localization have not been characterized in n-type polymers. Herein, we use near-edge X-ray absorption fine structure (NEXAFS) spectroscopy to probe electronic transitions in an electrochemically charged polymer to deduce charge localization. The O K-edge (1s) spectra exhibit two distinct π* peaks; the intensity of the lower-energy π* a peak that corresponds to an excitation to a largely localized carbonyl state decreases with reductive potentials relative to the higher-energy π* b peak. We corroborate this with Raman spectroscopy at different potentials, which shows a decrease in intensity on the C–C/C═C and C═O stretching bands of NDI as well as a red shift of the carbonyl band due to the formation of a polaron on the NDI. Additionally, new Raman active NDI signals associated with elongated C═O and C═C bonds are observed at lower energy during the formation of charged states. Together with theoretical calculations, these findings show that the injected charge spatially localizes on the NDI units and is dominantly distributed on the carbonyl groups. In conclusion, the combination of NEXAFS, optical and vibrational spectroscopies, and theoretical calculations is generalizable to other pi-conjugated polymers and can identify charge localization for the further development of organic semiconductors.

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

Structural Phase Transitions in the van der Waals Ferromagnets Fe x Pd y Te 2

Two-dimensional van der Waals (vdW) magnets are attracting significant attention, both as platforms for studying fundamental magnetic interactions and for the exciting possibility of utilizing them as building blocks in devices and heterostructures, which may lead to new physical phenomena and functionalities. Here, we provide a detailed study of the crystal structure and physical properties of the recently discovered vdW ferromagnet FePd 2 Te 2 . We find this compound has a relatively wide width of formation, and grow single crystals with compositions Fe x Pd y Te 2 where x ranges from 0.9 to 1.1 and y from 1.8 to 2.5, respectively. Temperature-dependent X-ray diffraction and transport measurements reveal that a first-order structural transition occurs in the range of T = 360–420 K, where the critical temperature, modulation wave vector, and corresponding room-temperature crystal structures all depend on chemical composition. Above the transition, the compounds with Pd fraction y > 2 adopt a disordered derivative of the tetragonal FeTe structure, with the Fe layer showing mixed Fe/Pd occupancy and the extra Pd atoms partially occupying interstitial sites. Below 370 K, the structure is incommensurately modulated, likely associated with the complex ordering of Pd/Fe atoms in the metal layers or the interstitial Pd in the vdW gaps. For y < 2, the composition Fe 1.1 Pd 1.8 Te 2 has monoclinic symmetry at room temperature that is consistent with the reported structure of FePd 2 Te 2 . This phase undergoes a structural transition at 420 K for which the high temperature structure is yet to be determined; however, based on the similarities with the y > 2 compounds, we speculate that its T > 420 K structure is also tetragonal. Importantly, the high temperature, symmetry-breaking structural transition observed here provides a likely explanation for the origin of the structural domains previously observed in FePd 2 Te 2 . All compounds investigated in the Fe x Pd y Te 2 series show metallic behavior, with magnetic characterization indicating that they are easy-plane, hard, ferromagnets with T C spanning 98–180 K. Both the critical temperature for the structural transition and the Curie temperature are moderately suppressed with increasing Pd fraction y and corresponding decreasing Fe fraction x, indicating that synthetic control over x and y paves way for the further exploration of these compounds.

crystal structure↗
DOE OSTI2020

Determination of the N–H Bond Dissociation Free Energy in a Pyridine(diimine)molybdenum Complex Prepared by Proton-Coupled Electron Transfer

The pyridine(diimine) molybdenum bis(imido) complex, ( iPr PDI)Mo(=NTol) 2 (Tol = 4-methylphenyl) was synthesized by addition of two equivalents of 4-methylphenylazide to a corresponding molybdenum benzene derivative, ( iPr PDI)Mo(η 6 -C 6 H 6 ) ( iPr PDI = 2,6-(2,6-iPr 2 C 6 H 3 N=CMe) 2 C 5 H 3 N). Protonation of ( iPr PDI)Mo(=NTol) 2 with 2,6-lutinidum triflate yielded a cationic molybdenum amido complex, [( iPr PDI)Mo(NHTol)(=NTol)][OTf], which was further transformed into the neutral molybdenum amido complex, ( iPr PDI)Mo(NHTol)(=NTol) by reduction with zinc powder. A series of spectroscopic, synthetic and pK a determination studies along with electrochemical measurements by the protonation-reduction pathway were used to establish an N–H bond dissociation free energy (BDFE) between 65-69 kcal/mol for the molybdenum imido-amido compound, ( iPr PDI)Mo(NHTol)(=NTol). Full molecule DFT studies provided a computed value of 61 kcal/mol. By contrast, reduction of ( iPr PDI)Mo(=NTol) 2 with KC 8 afforded the corresponding anionic molybdenum complex, K[( iPr PDI)Mo(=NTol) 2 ] that has a potassium cation intercalated with the pyridine and the tolyl groups. Protonation of K[( iPr PDI)Mo(=NTol) 2 ] with the weak amidinium acids [TBD(H)][BArF 24 ] (TBD = triazabicyclodecene, BArF 24 = B[3,5-(CF 3 ) 2 C 6 H 3 ] 4 ) also produced the neutral molybdenum amido complex, ( iPr PDI)Mo(NHTol)(=NTol). Measurement of the pKa and oxidation potential of K[( iPr PDI)Mo(=NTol) 2 ] provided a range of 69-73 kcal/mol for the N–H BDFE of ( iPr PDI)Mo(NHTol)(=NTol), in good agreement with the protonation-reduction route and completing the square scheme. The similar pK a and redox potentials obtained from each pathway demonstrate that both sequences are energetically feasible for PCET events. This study on the determination of N–H BDFE of the molybdenum amido complex renders fundamental insight into the N 2 reduction cycle by proton-coupled electron transfer.

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

Aqueous TiO 2 Nanoparticles React by Proton-Coupled Electron Transfer

Redox reactions of aqueous colloidal TiO 2 4 nm nanoparticles (NPs) have been examined, including both citratecapped and uncapped NPs (c-TiO 2 and uc-TiO 2 ). Photoreduction gave stable blue colloidal c-TiO 2 R NPs with 10–60 electrons per particle. Equilibration of these reduced NPs with soluble redox reagents such as methylviologen (MV 2+ ) provided measurements of the colloid reduction potential as a function of pH. The potentials of c-TiO 2 from pH 2–9 varied linearly with pH, with a slope of –60 ± 5 mV/pH. Estimates of the potential at pH 12 were consistent with extrapolating that line to high pH. The reduction potentials did not correlate with the zeta potentials (ζ) or the surface charge of the NPs across this pH range. Similar reduction potentials were observed for c- and uc-TiO 2 at low pH even though they have quite different ζ potentials. These results show that the common surface-charging explanation of the pH dependence is not tenable in these systems. Oxidation of reduced c-TiO 2 R with the electron-transfer oxidant potassium triiodide (KI 3 ) occurred with a significant drop in pH, showing that protons were released when the electrons were removed from the NPs. Smaller pH drops were observed for the proton-coupled electron transfer (PCET) reagents O 2 (air) and 4-MeO-TEMPO (4-methoxy-2,2,6,6-tetramethylpiperine-1-oxy radical). The difference in the number of protons released with KI 3 vs O 2 and 4-MeO-TEMPO was roughly one proton per electron removed. Thus, the thermodynamically preferred reactivity of these colloidal TiO 2 NPs is PCET over the pH 2–13 range studied. The measured redox potentials refer to the chemical process TiO 2 + H + + e – → TiO 2 ·e – ,H + ; and therefore they do not correspond with an electronic energy such as a conduction band edge or flat band potential. The 1e – /1H + stoichiometry means that the TiO 2 reduction potentials correspond to a TiO 2 –H bond dissociation free energy (BDFE), determined to be 49 ± 2 kcal mol –1 . The PCET description is consistent with the pH dependence of E(TiO 2 /TiO 2 ·e – ,H + ), the release of protons upon oxidation, the lack of correlation with ζ potentials, the similarity of capped and uncapped NPs, and the small change in the potential and BDFE from the first to the last electron/proton pair (H atom) removed. Furthermore, this behavior is suggested to be the norm for redox-active oxide/water interfaces.

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

Screened Range-Separated Hybrid Functional with Polarizable Continuum Model Overcomes Challenges in Describing Triplet Excitations in the Condensed Phase Using TDDFT

Long range-corrected (LRC) or range-separated hybrid (RSH) functionals where the long-range (LR) limit of electronic interactions is set to the exact exchange have been shown to correct the tendency of traditional density functional theory (DFT) to underestimate the frontier orbital gap. Consequently, the use of such functionals in calculating electronic excited states using linear response based time-dependent DFT (TDDFT) has been successful in correcting the tendency for underestimating the energies of charge transfer states by DFT-based calculations. More recently formulations of functionals that attenuate the LR limit to address condensed-phase effects to polarize the electronic density have been reported. In particular screened RSH (SRSH) combined with polarizable continuum model (PCM) was benchmarked successfully in reproducing the fundamental gap and charge transfer state energies of molecular systems in the condensed phase. Here we use SRSH-PCM to address triplet excited states, and show its success in obtaining correspondence of the low-lying triplet states to the singlet–triplet gap in a similar way that the fundamental orbital gap corresponds to electron removal and addition energies. Importantly, the accuracy of the SRSH-PCM in calculating triplet excitations stands on the polarization consistent framework in addressing the scalar dielectric constant and without affecting the optimal tuning by triplet energies. Finally, the prospect of even further improving the SRSH-PCM accuracy in calculating triplet states can be achieved by optimal tuning on the basis of the spin multiplicity gap.

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

Localized Active Space Pair-Density Functional Theory

Accurate quantum chemical methods for the prediction of spin-state energy gaps for strongly correlated systems are computationally expensive and scale poorly with the size of the system. This makes calculations for many experimentally interesting molecules impractical even with abundant computational resources. Previous work has shown that the localized active space (LAS) self-consistent field (SCF) method can be an efficient way to obtain multiconfiguration SCF wave functions of comparable quality to the corresponding complete active space (CAS) ones. To obtain quantitative results, a post-SCF method is needed to estimate the complete correlation energy. One such method is multiconfiguration pair-density functional theory (PDFT), which calculates the energy based on the density and on-top pair density obtained from a multiconfiguration wave function. In this work, we introduce localized-active-space PDFT, which uses a LAS wave function for subsequent PDFT calculations. The method is tested by computing spin-state energies and gaps in conjugated organic molecules and a bimetallic compound and comparing to the corresponding CAS-PDFT values.

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

DONKEY: A Flexible and Accurate Algorithm for Clustering

We propose an accurate clustering algorithm suitable for the varied and multidimensional data sets that correspond to temporal snapshots from on-the-fly nonadiabatic trajectory-based simulations of photoexcited dynamics. The algorithm approximates the underlying probability density function using variable kernel density estimation, with local maxima corresponding to cluster centers. Each data point is then assigned to one of the maxima by employing a maximization procedure. Finally, clusters artificially separated by minor fluctuations in the probability density are merged. The algorithm does not require parameter tuning, which ensures flexibility and reduces the risk of bias. It is tested on several synthetic data sets, where it consistently outperforms conventional clustering algorithms. As a final example, the algorithm is applied to the excited dynamics of the norbornadiene ⇌ quadricyclane (C 7 H 8 ) molecular photoswitch, demonstrating how distinct reaction pathways can be identified.

algorithms↗
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 OSTI2023

Excited-State Dynamics during Primary C–I Homolysis in Acetyl Iodide Revealed by Ultrafast Core-Level Spectroscopy

In typical carbonyl-containing molecules, bond dissociation events follow initial excitation to $nπ_{C=O}$$^*$ states. However, in acetyl iodide, the iodine atom gives rise to electronic states with mixed $nπ_{C=O}$$^*$ and $nπ_{C–I}$$^*$ character, leading to complex excited-state dynamics, ultimately resulting in dissociation. Using ultrafast extreme ultraviolet (XUV) transient absorption spectroscopy and quantum chemical calculations, we present an investigation of the primary photodissociation dynamics of acetyl iodide via time-resolved spectroscopy of core-to-valence transitions of the I atom after 266 nm excitation. The probed I 4d-to-valence transitions show features that evolve on sub-100-fs time scales, reporting on excited-state wavepacket evolution during dissociation. These features subsequently evolve to yield spectral signatures corresponding to free iodine atoms in their spin–orbit ground and excited states with a branching ratio of 1.1:1 following dissociation of the C–I bond. Calculations of the valence excitation spectrum via equation-of-motion coupled cluster with single and double substitutions (EOM-CCSD) show that initial excited states are of spin-mixed character. From the initially pumped spin-mixed state, we use a combination of time-dependent density functional theory (TDDFT)-driven nonadiabatic ab initio molecular dynamics and EOM-CCSD calculations of the N$_{4,5}$ edge to reveal a sharp inflection point in the transient XUV signal that corresponds to rapid C–I homolysis. Here, by examining the molecular orbitals involved in the core-level excitations at and around this inflection point, we are able to piece together a detailed picture of C–I bond photolysis in which d → σ* transitions give way to d → p excitations as the bond dissociates. We also report theoretical predictions of short-lived, weak 4d → 5d transitions in acetyl iodide, validated by weak bleaching in the experimental transient XUV spectra. This joint experimental–theoretical effort has thus unraveled the detailed electronic structure and dynamics of a strongly spin–orbit coupled system.

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

Cryogenic Photoelectron Spectroscopic and Theoretical Study of the Electronic and Geometric Structures of Undercoordinated Osmium Chloride Anions OsCl n – ( n = 3–5)

A series of anionic transition metal halides OsCl n - (n=3–5) have been investigated using a newly developed, home-constructed, cryogenic anion cluster photoelectron spectroscopy. Here, the target anionic species are generated through collision-induced dissociation in a two-stage ion funnel. The measured vertical detachment energies (VDEs) are 3.48 eV, 4.54 eV, and 4.81 eV for n = 3, 4, and 5, respectively. Density functional theory calculations at the B3LYP-D3(BJ)//aug-cc-pVTZ (-pp) level predict the lowest energy structures of OsCl n - (n = 3–5) to be a quintet triangle, quartet square, and quintet square-based pyramid, respectively. The CCSD(T)-calculated VDEs and corresponding adiabatic detachment energies (ADEs) agree well with our experimental measurements. Analysis of the corresponding frontier molecular orbitals (FMOs) and charge density differences suggests that the d-orbitals of the transition metal Os play a primary role in the single-photon detachment processes, and the detached electrons originating from different molecular orbitals are distinguishable.

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

Concordant Mode Approach (CMA): Vibrational Analysis of New and Upgraded Intermolecular Benchmarks for Noncovalent Bonding

The Concordant Mode Approach (CMA) is a novel method that offers tremendous potential for increasing the system size and the level of theory attainable in quantum chemical computations of molecular vibrational frequencies. To investigate the extension of CMA to intermolecular vibrations, computations with coupled cluster singles and doubles with perturbative triples theory [CCSD(T)] using two augmented correlation-consistent polarized-valence triple-ζ basis sets (aug-cc-pVTZ or h-aug-ccpVTZ) were performed on 17 prototypical loosely bound complexes of hydrogen-bonded, dispersion, and mixed character. These Level A results provide new and upgraded benchmarks for noncovalent bonding and a severe test for CMA vibrational analyses. The Level A target frequencies were recovered remarkably well using second-order Møller−Plesset perturbation theory (MP2) with h-aug-cc-pVTZ for generating the underlying (Level B) normal modes of the CMA scheme. Employing this Level B within the lowest-rung CMA-0A method reproduces the 435 benchmark frequencies with a mean absolute error (MAE) of 0.23 cm −1 and a corresponding standard deviation (σ) of 0.84 cm −1 ; strikingly, the corresponding subset of 106 interfragment frequencies exhibits MAE = 0.34 cm −1 and σ = 0.90 cm −1 . Subsequent application of the higher-rung CMA-2A scheme eliminates all outliers and reduces the overall MAE to a minuscule 0.08 cm −1 with the inclusion of only 3.0% of the off-diagonal couplings not accounted for by CMA-0A. Accordingly, the highly efficient CMA methodology proves to be robust even for vibrations on flat potential energy surfaces.

Aromatic compounds↗
DOE OSTI2023

Molecular Rotations, Multiscale Order, Hyperuniformity, and Signatures of Metastability during the Compression/Decompression Cycles of Amorphous Ices

We model, via large-scale molecular dynamics simulations, the isothermal compression of low-density amorphous ice (LDA) to generate high-density amorphous ice (HDA) and the corresponding decompression extending to negative pressures to recover the low-density amorphous phase (LDA HDA ). Both LDA and HDA are nearly hyperuniform and are characterized by a dynamical HBN, showing that amorphous ices are nonstatic materials and implying that nearly hyperuniformity can be accommodated in dynamical networks. In correspondence with both the LDA-to-HDA and the HDA-to-LDA HDA phase transitions, the (partial) activation of rotational degrees of freedom activates a cascade effect that induces a drastic change in the connectivity and a pervasive reorganization of the HBN topology which, ultimately, break the samples’ hyperuniform character. Key to this effect is the rapid rate at which changes occur, and not their magnitude. The inspection of structural properties from the short- to the long-range shows that signatures of metastability are present at all length-scales, hence providing further solid evidence in support of the liquid–liquid critical point scenario. LDA and LDA HDA differ in terms of HBN and structural properties, implying that they are distinct low-density glasses. Our work unveils the role of molecular rotations in the phase transitions between amorphous ices and shows how the unfreezing of rotational degrees of freedom generates a cascade effect that propagates over multiple length-scales. Our findings greatly improve our basic understanding of water and amorphous ices and can potentially impact the field of molecular network-forming materials at large.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Inelastic Neutron Scattering Study of Magnetic Exchange Pathways in MnS

We report an investigation of the magnetic structure and magnetic exchange pathways in MnS via neutron scattering methods, aided by density functional theory (DFT) modeling. The material has been confirmed to undergo antiferromagnetic (AFM) ordering at 152 K, with the magnetic structure representing AFM stacking of ferromagnetic (FM) (111) planes of Mn magnetic moments. Correspondingly, the magnetic structure is described by a propagation vector k = (1/2, 1/2, 1/2), with the volume of the magnetic unit cell being 8 times larger than the volume of the nuclear unit cell. Analysis of inelastic neutron scattering (INS) data collected on a powder sample of MnS revealed that the next-nearest-neighbor magnetic exchange constant (J 2 ) exceeds the nearest-neighbor exchange constant (J 1 ) by more than 3 times, while in the case of MnO, which exhibits the same nuclear and magnetic structures as MnS, the J 2 /J 1 ratio was reported to be below 1.5. Although for MnO the signs of both J 1 and J 2 indicated AFM exchange interactions, machine-learning INS data analysis in combination with DFT calculations suggests that the INS data collected on MnS are best described with J 1 < 0 and J 2 > 0, corresponding to FM and AFM exchange couplings, respectively. To achieve a satisfactory fit to the experimentally observed data, the Hamiltonian used to model the INS spectra also included the next-next-nearest-neighbor magnetic exchange constant (J 3 ). The best-fit model has been obtained with the values of the exchange constants J 1 = -0.27 meV, J 2 = 1.05 meV, and J 3 = -0.19 meV.

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

Frequency and Time Domain Nuclear-Electronic Orbital Equation-of-Motion Coupled Cluster Methods: Combination Bands and Electronic-Protonic Double Excitations

The accurate description of excited vibronic states is important for modeling a wide range of photoinduced processes. The nuclear-electronic orbital (NEO) approach, which treats specified protons on the same level as the electrons, can describe excited electronic-protonic states. Herein the multicomponent equation-of-motion coupled cluster with singles and doubles (NEO-EOM-CCSD) method and its time-domain counterpart, TD-NEO-EOM-CCSD, are developed and implemented. The application of these methods to the HCN molecule highlights their capabilities. These methods predict qualitatively reasonable energies and intensities for a combination band corresponding to simultaneous excitation of two vibrational modes, as well as an overtone. These methods also describe states with double excitation character, such as excited electronic-protonic states corresponding to the simultaneous excitation of an electron and a proton. The ability of the NEO-EOM-CCSD method and its time-dependent counterpart to describe combination bands, overtones, and double excitations will enable a wide range of photochemical applications.

Pavosevic, Fabijan↗