Search NASA⌕ Search

SEARCH · Search NASA

Results for “MOLECULAR ELECTRONICS”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

TaZr_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 tantalum-zirconium (Ta-Zr). 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 Ta and Zr. 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 80, 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. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10 and between 81 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. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of 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. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system 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. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & 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. & 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. & 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. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2024

Catalyst design for small molecule activation of energy consequence

This project targets the conversion of ubiquitous small molecules (e.g. NO, CO, H 2 O) into viable precursors to synthetic fuels. Current state of the art catalyst design has not directly targeted transition metal complexes capable of mediating the multi-electron redox processes necessary to reduce the overpotential (energy loss) required achieve efficient activation of small molecule substrates. In this vein, a new strategy has been developed for the assembly of polynuclear architectures; allowing for the construction of tunable polymetallic centers that assemble easily within a pre-organized template (conferring stability, selectivity and tunability) that can effect multi-electron redox processes for reactions. Catalyst development has commenced with the following target design elements: (1) catalysts featuring multiple transition metal ions in the same reaction space to greatly expand accessible molecular redox capabilities; (2) catalysts are assembled in a polynucleating ligand framework that permits control over the cluster morphology as well as the local steric and electronic environment of the transition metal ions within the cluster. The high-tunability of the catalyst composition (metal content) and geometric flexibility has permitted a rigorous assessment of electronic-structure-to-function relationship to be developed, further guiding synthetic efforts to realize more potent catalysts. The numerous permutations possible showcase the high degree of generality to this approach with many synthetic handles to tune redox and reaction chemistry. Trinuclear complexes have been synthesized featuring homo- and hetero-trinuclear cores featuring a variety of first row transition metal ions (Cr→Ni). The molecular clusters have been shown to successfully mediate multi-electron redox processes in a cooperative fashion without requiring strong chemical reductants or oxidants. The reactive molecular complexes are being utilized to activate and breakdown the robust bonds within typical waste stream small molecules (e.g., greenhouse gases) and convert them into value-added commodity chemicals. Ultimately, the catalysts developed by this approach will be required to convert energy acquired via renewable resources (e.g., solar or wind) into synthetic fuels as an energy storage mechanism.

10 SYNTHETIC FUELS↗
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 OSTI2024

Emergence of two distinct phase transitions in monolayer CoSe2 on graphene

Abstract Dimensional modifications play a crucial role in various applications, especially in the context of device miniaturization, giving rise to novel quantum phenomena. The many-body dynamics induced by dimensional modifications, including electron-electron, electron-phonon, electron-magnon and electron-plasmon coupling, are known to significantly affect the atomic and electronic properties of the materials. By reducing the dimensionality of orthorhombic CoSe 2 and forming heterostructure with bilayer graphene using molecular beam epitaxy, we unveil the emergence of two types of phase transitions through angle-resolved photoemission spectroscopy and scanning tunneling microscopy measurements. We disclose that the 2 × 1 superstructure is associated with charge density wave induced by Fermi surface nesting, characterized by a transition temperature of 340 K. Additionally, another phase transition at temperature of 160 K based on temperature dependent gap evolution are observed with renormalized electronic structure induced by electron-boson coupling. These discoveries of the electronic and atomic modifications, influenced by electron-electron and electron-boson interactions, underscore that many-body physics play significant roles in understanding low-dimensional properties of non-van der Waals Co-chalcogenides and related heterostructures. Graphical Abstract

36 MATERIALS SCIENCE↗
DOE OSTI2024

Tunable polar distortions and magnetism in Gd 𝑥 ⁢La 1−𝑥 ⁢PtSb epitaxial films

Hexagonal 𝐴⁢𝐵⁢𝐶 intermetallics are predicted to have tunable ferroelectric, topological, and magnetic properties as a function of the polar buckling of 𝐵⁢𝐶 atomic planes. Here, we report the impact of isovalent lanthanide substitution on the buckling, structural phase transitions, and electronic and magnetic properties of Gd 𝑥 ⁢La 1−𝑥 ⁢PtSb films grown by molecular beam epitaxy (MBE) on 𝑐 plane sapphire substrates. The Gd 𝑥 ⁢La 1−𝑥 ⁢PtSb films form a solid solution from 𝑥=0 to 𝑥=1 and retain the polar hexagonal structure (𝑃⁢6 3 ⁢𝑚⁢𝑐) out to 𝑥 ≤ 0.95. With increasing 𝑥, the PtSb buckling increases and the out-of-plane lattice constant 𝑐 decreases due to the lanthanide contraction. While hexagonal LaPtSb is a highly conductive polar metal, the carrier density decreases with 𝑥 until an abrupt phase transition to a zero band overlap semimetal is found for cubic GdPtSb at 𝑥=1. The magnetic susceptibility peaks at small but finite 𝑥, which we attribute to Ruderman-Kittel-Kasuya-Yosida (RKKY) coupling between localized 4⁢𝑓 moments, whose concentration increases with 𝑥, and free carriers that decrease with 𝑥. Samples with 𝑥≥0.3 show antiferromagnetic Curie-Weiss behavior and a Neel temperature that increases with 𝑥. The Gd 𝑥 ⁢La 1−𝑥 ⁢PtSb system provides opportunities to dramatically alter the polar buckling and concentration of local 4⁢𝑓 moments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Spectrochemical Series for Electron Spin Relaxation

Controlling the rate of electron spin relaxation in paramagnetic molecules is essential for contemporary applications in molecular magnetism and quantum information science. However, the physical mechanisms of spin relaxation remain incompletely understood, and new spectroscopic observables play an important role in evaluating spin dynamics mechanisms and structure–property relationships. Here, we use cryogenic magnetic circular dichroism (MCD) spectroscopy and pulse electron paramagnetic resonance (EPR) in tandem to examine the impact of ligand field (d–d) excited states on spin relaxation rates. We employ a broad scope of square-planar Cu(II) compounds with varying ligand field strength, including CuS 4 , CuN 4 , CuN 2 O 2 , and CuO 4 first coordination spheres. An unexpectedly strong correlation exists between spin relaxation rates and the average d–d excitation energy (R 2 = 0.97). The relaxation rate trends as the inverse 11th power of the excited-state energies, whereas simplified theoretical models predict only an inverse second power dependence. These experimental results directly implicate ligand field excited states as playing a critical role in the ground-state spin relaxation mechanism. Furthermore, ligand field strength is revealed to be a particularly powerful design principle for spin dynamics, enabling formation of a spectrochemical series for spin relaxation.

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

Elastic Screening of Pseudogauge Fields in Graphene

Lattice deformations in graphene couple to the low-energy electronic degrees of freedom as effective scalar and gauge fields. Using molecular dynamics simulations, we show that the optical component of the displacement field, i.e., the relative motion of different sublattices, contributes at equal order as the acoustic component and effectively screens the pseudogauge fields. In particular, we consider twisted bilayer graphene and corrugated monolayer graphene. In both cases, optical lattice displacements significantly reduce the overall magnitude of the pseudomagnetic fields. For corrugated graphene, optical contributions also reshape the pseudomagnetic field and significantly modify the electronic bands near charge neutrality. Previous studies based on continuum elasticity, which ignores this effect, have therefore systematically overestimated the strength of the strain-induced pseudomagnetic field. Furthermore, our results have important consequences for the interpretation of experiments and design of straintronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dissociative electron recombination and rotational cooling of the deuterated triatomic hydrogen ions H 2 ⁡D + and D 2 ⁢H +

We have measured the dissociative recombination (DR) of the deuterated triatomic hydrogen ions H 2 ⁡D + and D 2 ⁢H + as a function of storage time at the Cryogenic Storage Ring (CSR). Both molecular ions were stored for up to 1000 s inside the cryogenic vacuum of the CSR prior to the electron recombination measurements, allowing them to cool to their lowest rotational states. We implement a comprehensive model for all relevant processes to predict the internal state evolution of the ions during storage inside the CSR, employing calculated radiative transition strengths and state-selective rate coefficients for electron collisions. Our DR rate coefficient measurements with deuterated triatomic hydrogen ions in defined quantum states allow for meaningful comparisons with state-of-the-art theoretical calculations, paving the way for a better understanding of the complex DR process for polyatomic molecular ions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Ultrathin Boron Growth onto Nanodiamond Surfaces via Electrophilic Boron Precursors

Diamond as a templating substrate is largely unexplored, and the unique properties of diamond, including its large bandgap, thermal conductance, and lack of cytotoxicity, makes it versatile in emergent technologies in medicine and quantum sensing. Surface termination of an inert diamond substrate and its chemical reactivity are key in generating new bonds for nucleation and growth of an overlayer material. Oxidized high-pressure high temperature (HPHT) nanodiamonds (NDs) are largely terminated by alcohols that act as nucleophiles to initiate covalent bond formation when an electrophilic reactant is available. In this work, we demonstrate a templated synthesis of ultrathin boron on ND surfaces using trigonal boron compounds. Boron trichloride (BCl 3 ), boron tribromide (BBr 3 ), and borane (BH 3 ) were found to react with ND substrates at room temperature in inert conditions. BBr 3 and BCl 3 were highly reactive with the diamond surface, and sheet-like structures were produced and verified with electron microscopy. Surface-sensitive spectroscopies were used to probe the molecular and atomic structure of the ND constructs’ surface, and quantification showed the boron shell was less than 1 nm thick after 1–24 h reactions. Observation of the reaction supports a self-terminating mechanism, similar to atomic layer deposition growth, and is likely due to the quenching of alcohols on the diamond surface. X-ray absorption spectroscopy revealed that boron-termination generated midgap electronic states that were originally predicted by density functional theory (DFT) several years ago. DFT also predicted a negative electron surface, which has yet to be confirmed experimentally here. The boron-diamond nanostructures were found to aggregate in dichloromethane and were dispersed in various solvents and characterized with dynamic light scattering for future cell imaging or cancer therapy applications using boron neutron capture therapy (BNCT). The unique templating mechanism based on nucleophilic alcohols and electrophilic trigonal precursors allows for covalent bond formation and will be of interest to researchers using diamond for quantum sensing, additive manufacturing, BNCT, and potentially as an electron emitter.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Transmission spectroscopy of CF 4 molecules in intense x-ray fields

The nonlinear interaction of x rays with matter is at the heart of understanding and controlling ultrafast molecular dynamics from an atom-specific viewpoint, providing new scientific and analytical opportunities to explore the structure and dynamics of small quantum systems. At increasingly high x-ray intensity, the sensitivity of ultrashort x-ray pulses to specific electronic states and emerging short-lived transient intermediates is of particular relevance for our understanding of fundamental multiphoton absorption processes. In this work, intense x-ray free-electron laser (XFEL) pulses at the European XFEL are combined with a gas cell and grating spectrometer for a high-intensity transmission spectroscopy study of multiphoton-induced ultrafast molecular fragmentation dynamics in CF 4 . This approach unlocks the direct intrapulse observation of transient fragments, including neutral atoms, by their characteristic absorption lines in the transmitted broadband x-ray spectrum. The dynamics with and without initially producing fluorine 𝐾-shell holes are studied by tuning the central photon energy. The absorption spectra are measured at different FEL intensities to observe nonlinear effects. Transient isolated fluorine atoms and ions are spectroscopically recorded within the ultrashort pulse duration of a few tens of femtoseconds. An isosbestic point that signifies the correlated transition between intact neutral CF 4 molecules and charged atomic fragments is observed near the fluorine 𝐾 edge. The dissociation dynamics and the multiphoton absorption-induced dynamics encoded in the spectra are theoretically interpreted. Overall, this study demonstrates the potential of high-intensity x-ray transmission spectroscopy to study ultrafast molecular dynamics with sensitivity to specific intermediate species and their electronic structure.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Molecular engineering of ethereal electrolyte for ultrastable Si-based high voltage full cells

The successful application of Si-based high-energy Li-ion batteries (LIBs) depends on our ability to tailor electrolyte properties to achieve long-term stability and reliable performance. In this work, we demonstrate our rationale for the molecular design of ethereal solvents to address low anodic stability issues and produce a highly electrochemically stable electrolyte for Si‖LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811) high-energy full cells. Unlike the trimethylsilyl group, the trifluoromethyl (–CF 3 ) group exerts a very strong electron-withdrawing effect on the glycol ether backbone, reducing the highest occupied molecular orbital (HOMO) energy level of the fluorinated glycol ether (FGE) and significantly enhancing its oxidation potential. The FGE-based electrolyte enables stable cycling of Si‖NMC811 full cells, delivering high specific capacity (900 mA h g −1 ) and coulombic efficiency (>99.78%) over extended (500) cycles. The improved electrochemical performance originates from the terminal fluorination of the diglyme backbone, which strengthens anion coordination in the solvation structure, leading to the preferential reduction of the FSI anion and the formation of robust solid electrolyte interphases (SEIs) on the Si surface. Through molecular engineering of ethereal solvents, we have discovered a promising candidate for a next-generation stable electrolyte, paving the way for the design of practical and commercially viable Si batteries.

Silicon anode↗
DOE OSTI2025

Convergent close-coupling calculations of electron scattering on LiH

Cross sections for electron scattering on the ground state of LiH are calculated using the molecular convergent close-coupling method. The fixed-nuclei approximation is utilized, and calculations are performed at the mean internuclear separation of 3.06 𝑎 0 for projectile energies ranging from 0.1 to 500 eV. Here, a model-potential approach is utilized to treat the LiH molecule as a two-electron system, and good agreement is found with the literature for various properties of the structure model. Cross sections are presented for elastic scattering, excitation to several electronic states, and ionization. Comparison is made with existing results, and qualitative agreement is found with small close-coupling calculations at low projectile energies. Good agreement is found with first-order results at high projectile energies.

74 ATOMIC AND MOLECULAR PHYSICS↗
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 OSTI2024

Nanoscale Control of Intrinsic Magnetic Topological Insulator MnBi 2 Te 4 Using Molecular Beam Epitaxy: Implications for Defect Control

Intrinsic magnetic topological insulators have emerged as a promising platform to study the interplay between the topological surface states and ferromagnetism. This unique interplay can give rise to a variety of exotic quantum phenomena, including the quantum anomalous Hall effect and axion insulating states. Here, in this study, utilizing molecular beam epitaxy (MBE), we present a comprehensive study of the growth of MnBi 2 Te 4 thin films on Si (111), epitaxial graphene, and highly ordered pyrolytic graphite substrates. By combining a suite of in situ characterization techniques, we obtain critical insights into the nanoscale control of MnBi 2 Te 4 epitaxial growth. First, we extract the free energy landscape for the epitaxial relationship as a function of the in-plane angular distribution. Then, by employing an optimized layer-by-layer growth, we determine the chemical potential and Dirac point of the thin film at different thicknesses and how this quantity is manifested by the dopant compensation from different antisite defects. Overall, these results establish a foundation for understanding the growth kinetics of MnBi 2 Te 4 and pave the way for future applications of MBE-grown thin films in emerging topological quantum materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Thorium Monosilicide, ThSi: An Experimental and Theoretical Study

The present theoretical and experimental combination study investigates the ThSi molecule in detail. Computationally, we utilized high-level multireference and coupled-cluster levels of theory conjoined with large correlation consistent basis sets to study a series of electronic and spin–orbit states of ThSi. Here, we report potential energy curves (PECs), electron configurations at equilibrium distances, spectroscopic constants, energetics, and spin–orbit coupling effects for 16 electronic states of ThSi. The studied 16 electronic states are arranged tightly within 0.9 eV, highlighting the complexity of the electronic spectrum of ThSi. The ground electronic state of ThSi is a single-reference 1 1 Σ + state that derives from the 1σ 2 2σ 2 1π 4 electronic configuration. The Ω = 0 + spin–orbit ground state of ThSi is composed of 1 1 Σ + (47%) and 13Π (44%) electronic states. Our measured bond energy (D0) of ThSi, obtained using resonant two-photon ionization (R2PI) spectroscopy is 3.146(4) eV, where the assigned error limit is given in parentheses in units of the last quoted digits. The computed D0 of ThSi (Ω = 0 + ) at the CBS-C-CCSD(T)-δT(Q)-δDK-δSO level (3.181 eV) is in good agreement with the experimental value. Our derived enthalpy of formation for ThSi, Δ f H 0K o (ThSi(g)), is 971.8(6.0) kJ/mol. Finally, we have performed density functional theory (DFT) calculations for ThSi(1 1 Σ + ) using 16 exchange correlation functionals that span multiple rungs of “Jacob’s ladder” of density functional approximation (DFA) to assess the DFT errors on D 0 , r e , and ω e of ThSi with respect to experimental and ab initio coupled-cluster values.

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

MixPI: Mixed-time slicing path integral software for quantized molecular dynamics simulations

We introduce the MixPI software to implement path integral molecular dynamics (PIMD) simulations for the study of condensed phase systems where nuclear quantum effects (NQEs) are important. In contrast to existing PIMD simulation software, MixPI enables the implementation of mixed quantum–classical path integral simulations where only a subset of system degrees of freedom (dofs) are treated quantum mechanically in an extended phase space while the remaining dofs are described classically. We expect this software to be particularly useful for simulations of electron and proton transfer in condensed phase systems, as well as for the study of biological and material systems where only a handful of dofs contribute significantly to the observed NQEs. We demonstrate the use of MixPI in two different systems. The first is a simple water model where we implement a set of mixed quantum–classical simulations to compute average energy and radial distribution functions. We use these simulations to benchmark the effectiveness of MixPI and to demonstrate how it enables systematic investigation into the origin of observed NQEs. We then compute radial distribution functions for a system where MixPI is essential: a solvated metal (M 2+ ) cation described using an explicit quantized electron localized on an M 3+ ion in water.

chemical physics↗
DOE OSTI2024

Ultrafast temporal phase-resolved nonlinear optical spectroscopy in the molecular frame

In an ultrafast nonlinear optical interaction, the electric field of the emitted nonlinear signal provides direct access to the induced nonlinear transient polarization or transient currents and thus carries signatures of ultrafast dynamics in a medium. Measurement of the electric field of such signals offers sensitive observables to track ultrafast electron dynamics in various systems. In this work, we resolve the real-time phase of the electric field of a femtosecond third-order nonlinear optical signal in the molecular frame. The electric field emitted from impulsively pre-aligned gas-phase molecules at room temperature, in a degenerate four-wave mixing scheme, is measured using a spectral interferometry technique. The nonlinear signal is measured around a rotational revival to extract its molecular-frame angle dependence from pump-probe time-delay scans. By comparing these measurements for two linear molecules, carbon dioxide and nitrogen, we show that the measured second-order phase parameter (temporal chirp) of the signal is sensitive to the valence electronic symmetry of the molecules, whereas the amplitude of the signal does not show such sensitivity. We compare measurements to theoretical calculations of the chirp observable in the molecular frame. This work is an important step towards using electric field measurements in nonlinear optical spectroscopy to study ultrafast dynamics of electronically excited molecules in the molecular frame.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Spectral deconvolution without the deconvolution: Extracting temperature from x-ray Thomson scattering spectra without the source-and-instrument function

X-ray Thomson scattering (XRTS) probes the dynamic structure factor of the system, but the measured spectrum is broadened by the combined source-and-instrument function (SIF) of the setup. In order to extract properties such as temperature from an XRTS spectrum, the broadening by the SIF needs to be removed. Recent work [Dornheim et al. Nat. Commun. 13 , 7911 (2022)] has suggested that the SIF may be deconvolved using the two-sided Laplace transform. However, the extracted information can depend strongly on the shape of the input SIF, and the SIF is in practice challenging to measure accurately. Here, we propose an alternative approach: we demonstrate that considering ratios of Laplace-transformed XRTS spectra collected at different scattering angles is equivalent to performing the deconvolution, but without the need for explicit knowledge of the SIF. From these ratios, it is possible to directly extract the temperature from the scattering spectra, when the system is in thermal equilibrium. We find the method to be generally robust to spectral noise and physical differences between the spectrometers, and we explore situations in which the method breaks down. Furthermore, the fact that consistent temperatures can be extracted for systems in thermal equilibrium indicates that non-equilibrium effects could be identified by inconsistent temperatures of a few eV between the ratios of three or more scattering angles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗