Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex systems theory”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Spin-piston problem for a ferromagnetic thin film: Shock waves and solitons

In this report the unsteady, nonlinear magnetization dynamics induced by spin injection in an easy-plane ferromagnetic channel subject to an external magnetic field are studied analytically. Leveraging a dispersive hydrodynamic description, the Landau-Lifshitz equation is recast in terms of hydrodynamic-like variables for the magnetization's perpendicular component (spin density) and azimuthal phase gradient (fluid velocity). Spin injection acts as a moving piston that generates nonlinear, dynamical spin textures in the ferromagnetic channel with downstream quiescent spin density set by the external field. In contrast to the classical problem of a piston accelerating a compressible gas, here, variable spin injection and field lead to a rich variety of nonlinear wave phenomena from oscillatory spin shocks to solitons and rarefaction (expansion) waves. A full classification of solutions is provided using nonlinear wave modulation theory by identifying two key aspects of the fluid-like dynamics: subsonic/supersonic conditions and convex/nonconvex hydrodynamic flux. Familiar waveforms from the classical piston problem such as rarefaction waves and shocks manifest in their spin-based counterparts as smooth and highly oscillatory transitions, respectively, between two distinct magnetic states. The spin shock is an example of a dispersive shock wave, which arises in many physical systems. New features without a gas dynamics counterpart include composite wave complexes with “contact” spin shocks and rarefactions. Magnetic supersonic conditions lead to two pronounced piston edge behaviors including a stationary soliton and an oscillatory wave train. These coherent wave structures have physical implications for the generation of high frequency spin waves from pulsed injection and persistent, stable stationary and/or propagating solitons in the presence of magnetic damping. The analytical results are favorably compared with numerical simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Block Mott insulating state induced by next-nearest-neighbor hopping in the S = 3 2 zigzag chain BaCoTe 2 O 7

Quasi-one-dimensional correlated electronic multiorbital systems with either ladder or chain geometries continue attracting considerable interest due to their complex electronic phases arising from the interplay of the hopping matrix, the crystal-field splitting, the electronic correlations (Hubbard repulsion U and Hund coupling J H ), and strong quantum fluctuations. Recently, the intriguing cobalt zigzag chain system BaCoTe 2 ⁢ O 7 , with electronic density n=7, was prepared experimentally. Here, we systematically study the electronic and magnetic properties of this quasi-one-dimensional compound from the theoretical perspective. Based on first-principles density functional theory calculations, strongly anisotropic one-dimensional electronic Co 3⁢d bands were found near the Fermi level. By evaluating the relevant hopping amplitudes, we provide the magnitude and origin of the nearest-neighbor (NN) and next-nearest-neighbor (NNN) hopping matrices in BaCoTe 2 ⁢ O 7 . With this information, we constructed a three-orbital electronic Hubbard model for this zigzag chain system, and studied two cases: with only a NN hopping matrix, and with NN plus NNN hopping matrices. Introducing the Hubbard and Hund couplings and studying the model via the density matrix renormalization group method, we constructed the ground-state phase diagram. A robust staggered ↑ - ↓ - ↑ - ↓ antiferromagnetic (AFM) region was found when only the NN hopping matrix in the chain direction was employed. However, for the realistic case where the NNN hopping matrix is also included, the dominant state becomes instead a block AFM ↑ - ↑ - ↓ - ↓ order, in agreement with experiments. The system displays Mott insulator characteristics with three half-filled orbitals, when the block AFM order is stable. In conclusion, our results for BaCoTe 2 ⁢O 7 provide guidance to experimentalists and theorists working on this zigzag one-dimensional chain and related materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Photodriven electron-transfer dynamics in a series of heteroleptic Cu(I)–anthraquinone dyads

Solar fuels catalysis is a promising route to efficiently harvesting, storing, and utilizing abundant solar energy. To achieve this promise, however, molecular systems must be designed with sustainable components that can balance numerous photophysical and chemical processes. To that end, here we report on the structural and photophysical characterization of a series of Cu(I)–anthraquinone-based electron donor–acceptor dyads. The dyads utilized a heteroleptic Cu(I) bis-diimine architecture with a copper(I) bis-phenanthroline chromophore donor and anthraquinone electron acceptor. We characterized the structures of the complexes using x-ray crystallography and density functional theory calculations and the photophysical properties via resonance Raman and optical transient absorption spectroscopy. The calculations and resonance Raman spectroscopy revealed that excitation of the Cu(I) metal-to-ligand charge-transfer (MLCT) transition transfers the electron to a delocalized ligand orbital. The optical transient absorption spectroscopy demonstrated that each dyad formed the oxidized copper–reduced anthraquinone charge-separated state. Unlike most Cu(I) bis-phenanthroline complexes where increasingly bulky substituents on the phenanthroline ligands lead to longer MLCT excited-state lifetimes, here, we observe a decrease in the long-lived charge-separated state lifetime with increasing steric bulk. The charge-separated state lifetimes were best explained in the context of electron-transfer theory rather than with the energy gap law, which is typical for MLCT excited states, despite the complete conjugation between the phenanthroline and anthraquinone moieties.

14 SOLAR ENERGY↗
DOE OSTI2022

A 9.2-GHz clock transition in a Lu(II) molecular spin qubit arising from a 3,467-MHz hyperfine interaction

Spins in molecules are particularly attractive targets for next-generation quantum technologies, enabling chemically programmable qubits and potential for scale-up via self-assembly. Here, we report observation of one of the largest hyperfine interactions for a molecular system, A iso = 3467±50 MHz, along with an associated clock transition of unprecedented magnitude. This is achieved through chemical control of the degree of s-orbital mixing into the spin-bearing d-orbital associated with a series of spin-½ La(II) and Lu(II) complexes. Increased s-orbital character reduces spin-orbit coupling and enhances the electron-nuclear Fermi contact interaction. Both outcomes are advantageous for quantum applications: the former reduces spin-lattice relaxation, while the latter maximizes the hyperfine interaction that, in turn, generates a 9 gigahertz clock transition, leading to an increase in phase memory time from 1.0±0.4 to 12±1 microseconds for one of the Lu(II) complexes. Furthermore, these findings suggest strategies for development of molecular quantum technologies, akin to trapped ion systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

tapir: A tool for topologies, amplitudes, partial fraction decomposition and input for reductions

The demand for precision predictions in the field of high energy physics has dramatically increased over recent years. Experiments conducted at the LHC, as well as precision measurements at the intensity frontier such as Belle II require equally precise theoretical predictions to make full use of the acquired data. To match the experimental precision, second-, third- and, for certain quantities, even higher-order calculations in perturbative quantum field theory are required. To facilitate such calculations, computer software automating as many steps as possible is required. Yet, each calculation poses different challenges and thus, a high level of configurability is required. In this context we present tapir: a tool for identification, manipulation and minimization of Feynman integral families. It is designed to integrate in toolchains based on the computer algebra system FORM, the use of which is common practice in the field. tapir can be used to reduce the complexity of multi-loop problems with cut-filters, topology mapping, partial fraction decomposition and alike. Program Title:tapir CPC Library link to program files:https://doi.org/10.17632/ptc9t46xyn.1 Developer's repository link:https://gitlab.com/tapir-devs/tapir Licensing provisions: GPLv3 Programming language:python 3, C++ Nature of problem: Multi-loop computations require the automatization of a large number of different tasks related to Feynman integral topologies. Among them are the identification and minimization of integral topologies, partial fraction decomposition of topologies in the case of linearly dependent propagators as well as mapping scalar products of loop momenta to scalar functions. Solution method: The minimization of topologies is performed by comparison of their respective Nickel indices [1], even further minimization utilizes Pak's algorithm [2]. To efficiently map scalar products of loop momenta to scalar functions FORM [3] code is generated. Additional comments including restrictions and unusual features: Minimization based on Pak's algorithm slows down for many lines and scales. A coarser minimization using the Nickel indices, however, is still possible. [1]B. Nickel, D. Meiron, G.A.J. Baker, Compilation of 2-pt and 4-pt graphs for continuous spin model, Report, University of Guelph, 1977.[2]A. Pak, J. Phys. Conf. Ser. 368 (2012) 012049, https://doi.org/10.1088/1742-6596/368/1/012049, arXiv:1111.0868.[3]B. Ruijl, T. Ueda, J. Vermaseren, FORM version 4.2, arXiv:1707.06453, 7 2017.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Systematic Evaluation of Counterpoise Correction in Density Functional Theory

A widespread belief persists that the Boys–Bernardi function counterpoise (CP) procedure “overcorrects” supramolecular interaction energies for the effects of basis-set superposition error. To the extent that this is true for correlated wave function methods, it is usually an artifact of low-quality basis sets. The question has not been considered systematically in the context of density functional theory, however, where basis-set convergence is generally less problematic. We present a systematic assessment of the CP procedure for a representative set of functionals and basis sets, considering both benchmark data sets of small dimers and larger supramolecular complexes. The latter include layered composite polymers with ~150 atoms and ligand–protein models with ~300 atoms. Provided that CP correction is used, we find that intermolecular interaction energies of nearly complete-basis quality can be obtained using only double-ζ basis sets. Furthermore, this is less expensive as compared to triple-ζ basis sets without CP correction. CP-corrected interaction energies are less sensitive to the presence of diffuse basis functions as compared to uncorrected energies, which is important because diffuse functions are expensive and often numerically problematic for large systems. Our results upend the conventional wisdom that CP “overcorrects” for basis-set incompleteness. In small basis sets, CP correction is mandatory in order to demonstrate that the results do not rest on error cancellation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Iron Carbonyl Complexes of [2.2.2]Hericene as a Rigid Tris(1,3-diene) Ligand

Hericene is an unusual hexaolefin consisting of three 1,3-diene units located on a rigid bicyclo [2.2.2]octane framework that restricts the geometrical relationships of metal atoms bonded to these olefinic units. In order to explore possible effects of this rigidity limiting metal–metal interaction in polynuclear derivatives possibly stabilizing coordinatively unsaturated species, the structures and energetics of the hericene iron carbonyl complexes (hericene)Fem(CO)n (m = 1, n = 3; m = 2, n = 6, 5; m = 3, n = 9, 8) have been investigated by density functional theory. The lowest-energy (hericene)Fem(CO)3m (m = 1, 2, 3) structures have the cavities of the hericene ligand filled with a single Fe(CO)3 moiety bonded to a 1,3-diolefin unit. Such species have been synthesized by the reaction of Fe2(CO)9 with hericene. For the (hericene)Fe2(CO)5 system, the lowest energy structures are singlet structures with an Fe(CO)3 unit bonded to a 1,3-diene unit in one hericene cavity and an Fe(CO)2 unit in another hericene cavity bonded to three C=C double bonds from two 1,3-diene units. Higher energy (hericene)Fe2(CO)5 structures include a structure in which a single hericene cavity contains a Fe2(CO)4(µ-CO) moiety with each iron atom bonded to a 1,3-diene unit. In addition, both singlet and triplet (hericene)Fe2(CO)5 structures are found in which an Fe(CO)3 moiety and an Fe(CO)2 moiety located in separate hericene cavities are each bonded to a 1,3-diene unit. The lowest-energy (hericene)Fe3(CO)8 structures have two hericene cavities containing Fe(CO)3 moieties fully bonded to 1,3-diene units and a third hericene cavity containing an Fe(CO)2 moiety fully bonded to a 1,3-diene unit.

Luo, Jinfeng↗
DOE OSTI2021

Ruling Out the Onset of Color Transparency up to Q2 = 14.2 GeV2 in Quasielastic 12C(e,e'p) Scattering

Color Transparency (CT) is a prediction of QCD that at high momentum transfer Q2, a system of quarks which would normally interact strongly with nuclear matter could form a small color-neutral object whose compact transverse size would be maintained for some distance, passing through the nucleus undisturbed. A clear signature of CT would be a dramatic rise in nuclear transparency T with increasing Q2. CT emerges as a deviation from Glauber multiple scattering theory, which predicts constant T. While a rise in nuclear transparency would provide an unequivocal validation of QCD factorization theorems, the complex nature of nuclear interactions renders its observation difficult to predict. The E12-06-107 experiment at JLab measured T in quasielastic electron-proton scattering with carbon-12 and liquid hydrogen targets, for Q2 between 8.0 and 14.2 GeV2 , a range over which models of CT predicted that T might differ appreciably from Glauber calculations. Supported in part by US DOE grant DE-FG02-03ER41240.

Matter, John↗
DOE OSTI2024

Importance of Electron Correlation on the Geometry and Electronic Structure of [2Fe–2S] Systems: A Benchmark Study of the [Fe 2 S 2 (SCH 3 ) 4 ] 2–,3–,4– , [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S- p -tol) 4 ] 2– , and [Fe 2 S 2 (S- o -xyl) 4 ] 2– Complexes

Iron–sulfur clusters are crucial for biological electron transport and catalysis. Obtaining accurate geometries, energetics, manifolds of their excited electronic states, and reduction energies is important to understand their role in these processes. Here, using a [2Fe–2S] model complex with Fe II and Fe III oxidation states, which leads to different charges, i.e., [Fe 2 S 2 (SMe) 4 ] 2–,3–,4– , we benchmarked a variety of computational methodologies ranging from density functional theory (DFT) to post-Hartree–Fock methods, including complete active space self-consistent field (CASSCF), multireference configuration interaction, the second-order N-electron valence state perturbation theory (NEVPT2), and the linearized integrand approximation of adiabatic connection (AC0) approaches. Additionally, we studied three experimentally well-characterized complexes, [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S-o-tol) 4 ] 2– , and [Fe 2 S 2 (S-o-xyl) 4 ] 2– , via DFT methods. We conclude that the dynamic electron correlation is important for accurately predicting the geometry of these complexes. Broken symmetry (BS) DFT correctly predicts experimental geometries of low-spin multiplicity, while CASSCF does not. However, BS-DFT significantly overestimates the difference between the low- and high-spin electronic states for a given oxidation state. At the same time, CASSCF underestimates it but provides relative energies closer to the reference NEVPT2 results. Finally, AC0 provides energetics of NEVPT2 quality with the additional advantage of being able to use large CASSCF sizes. NEVPT2 gives the best estimates of the Fe III /Fe III → Fe II /Fe III (4.27 eV) and Fe II /F III → Fe II /F II (7.72 eV) reduction energies. The results provide insight into the electronic structure of these complexes and assist in the understanding of their physical properties.

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

The Cauchy Problem for Boltzmann Bi-linear Systems: The Mixing of Monatomic and Polyatomic Gases

Abstract From a unified vision of vector valued solutions in weighted Banach spaces, this paper establishes the existence and uniqueness for space homogeneous Boltzmann bi-linear systems with conservative collisional forms arising in complex gas dynamical structures. This broader vision is directly applied to dilute multi-component gas mixtures composed of both monatomic and polyatomic gases. Such models can be viewed as extensions of scalar Boltzmann binary elastic flows, as much as monatomic gas mixtures with disparate masses and single polyatomic gases, providing a unified approach for vector valued solutions in weighted Banach spaces. Novel aspects of this work include developing the extension of a general ODE theory in vector valued weighted Banach spaces, precise lower bounds for the collision frequency in terms of the weighted Banach norm, energy identities, angular or compact manifold averaging lemmas which provide coerciveness resulting into global in time stability, a new combinatorics estimate for p -binomial forms producing sharper estimates for the k -moments of bi-linear collisional forms. These techniques enable the Cauchy problem improvement that resolves the model with initial data corresponding to strictly positive and bounded initial vector valued mass and total energy, in addition to only a $$2^+$$ 2 + moment determined by the hard potential rates discrepancy, a result comparable in generality to the classical Cauchy theory of the scalar homogeneous Boltzmann equation.

Physics↗
DOE OSTI2025

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Applications to Materials for Catalysis and Electrochemistry

Modeling complex materials using high-fidelity, ab initio methods at low cost is a fundamental goal for quantum chemical software packages. The GW approximation and random phase approximation (RPA) provide a unified description of both electronic structure and total energies using the same physics in a many-body perturbative approach that can be more accurate than generalized-gradient density functional theory (DFT) methods. However, GW/RPA implementations have historically been limited to either specific materials classes or application toward small chemical systems. Here, the static subspace approximation allows for reduced cost full-frequency GW/RPA calculations and has previously been benchmarked thoroughly for GW calculations. Here, we describe our approach to including partial occupations of electronic orbitals in full-frequency GW and RPA calculations for the study of electrocatalysts. We benchmarked RPA total energy calculations using the subspace approximation across a diverse test suite of materials for a variety of computational parameters. The benchmarking quantifies the impact of different extrapolation procedures for representing the static polarizability at infinite screened cutoff, and shows that using screened cutoffs above 20-25 Ryd result in diminishing accuracy returns for predicting RPA total energies. Additionally, for moderately sized electrocatalytic models, 2-3 times fewer computational resources are used to compute RPA total energies by representing the static polarizability with 20-30% of the static subspace basis, with an error of approximately 0.01 eV or better in RPA adsorption energy calculations. Finally, we show that for these electrochemical models RPA can shift DFT adsorption energy shifts by up to 0.5 eV and that GW can frequently shift DFT eigenvalues of surface and adsorbate states by approximately 0.5-1 eV.

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

Dynamic and relativistic effects on Pt–Pt indirect spin–spin coupling in aqueous solution studied by ab initio molecular dynamics and two- vs four-component density functional NMR calculations

Treating 195Pt nuclear magnetic resonance parameters in solution remains a considerable challenge from a quantum chemistry point of view, requiring a high level of theory that simultaneously takes into account the relativistic effects, the dynamic treatment of the solvent–solute system, and the dynamic electron correlation. A combination of Car-Parrinello molecular dynamics (CPMD) and relativistic calculations based on two-component zeroth order regular approximation spin–orbit Kohn–Sham (2c-ZKS) and four-component Dirac–Kohn–Sham (4c-DKS) Hamiltonians is performed to address the solvent effect (water) on the conformational changes and JPtPt1 coupling. A series of bridged PtIII dinuclear complexes [L1–Pt2(NH3)4(Am)2–L2]n+ (Am = α–pyrrolidonate and pivalamidate; L = H2O, Cl−, and Br−) are studied. The computed Pt–Pt coupling is strongly dependent on the conformational dynamics of the complexes, which, in turn, is correlated with the trans influence among axial ligands and with the angle N–C–O from the bridging ligands. The J-coupling is decomposed in terms of dynamic contributions. The decomposition reveals that the vibrational and explicit solvation contributions reduce JPtPt1 of diaquo complexes (L1 = L2 = H2O) in comparison to the static gas-phase magnitude, whereas the implicit solvation and bulk contributions correspond to an increase in JPtPt1 in dihalo (L1 = L2 = X−) and aquahalo (L1 = H2O; L2 = X−) complexes. Relativistic treatment combined with CPMD shows that the 2c-ZKS Hamiltonian performs as well as 4c-DKS for the JPtPt1 coupling.

Chemistry↗
DOE OSTI2022

Chemistry-mediated Ostwald ripening in carbon-rich C/O systems at extreme conditions

Abstract There is significant interest in establishing a capability for tailored synthesis of next-generation carbon-based nanomaterials due to their broad range of applications and high degree of tunability. High pressure (e.g., shockwave-driven) synthesis holds promise as an effective discovery method, but experimental challenges preclude elucidating the processes governing nanocarbon production from carbon-rich precursors that could otherwise guide efforts through the prohibitively expansive design space. Here we report findings from large scale atomistically-resolved simulations of carbon condensation from C/O mixtures subjected to extreme pressures and temperatures, made possible by machine-learned reactive interatomic potentials. We find that liquid nanocarbon formation follows classical growth kinetics driven by Ostwald ripening (i.e., growth of large clusters at the expense of shrinking small ones) and obeys dynamical scaling in a process mediated by carbon chemistry in the surrounding reactive fluid. The results provide direct insight into carbon condensation in a representative system and pave the way for its exploration in higher complexity organic materials. They also suggest that simulations using machine-learned interatomic potentials could eventually be employed as in-silico design tools for new nanomaterials.

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

Impact of intermolecular interactions on the spectroscopic signals, energetics, and redox behavior of high valent 237 Np (Final Report)

Our studies have led to definitive spectral assignments for neptunyl-neptunyl interactions in aqueous solutions that will be valuable for predicting Np(V) speciation in solution. We have definitively demonstrated how actinyl-hydrogen and actinyl-cation interactions can influence spectral features for uranyl and neptunyl system. In addition, we have developed DFT methodologies to predict enthalpies of formation for uranyl and neptunyl systems and linked stability to chemical descriptors (hydrogen interaction energies or the electrostatic attraction energies). The counterion influenced the rate of Np(V) oxidation to Np(VI) neptunates and [Np(VII)O 4 OH] 3- complexes. Oxidation to Np(VII) during ozonolysis takes place through the presence of hydroxyl radicals.

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

An Ab Initio Multiple Cloning Method for Non-Adiabatic Excited-State Molecular Dynamics in NWChem

The recently developed ab-initio multiple cloning (AIMC) approach based on the multicongurational Ehrenfest (MCE) method provides a powerful and accurate way of describing the excited-state dynamics of molecular systems. The AIMC method is a controlled approximation to non-adiabatic dynamics with a particular strength in the proper description of decoherence effects due to branching of vibrational wavepackets at a level crossing. In this work, we report a new implementation of the AIMC algorithm in the open source NWChem computational chemistry program. The framework combines linear-response time-dependent density functional theory with Ehrenfest mean-eld theory to determine the equations of motion for classical trajectories. The multi-dimensional wave function is decomposed into a superposition of Gaussian coherent states guided by Ehrenfest trajectories (i.e., MCE approach), which can clone with fully quantum mechanical amplitudes and phases. By using an efficient time-derivative based non-adiabatic coupling approach within the AIMC method, all observable are calculated on-the-y in the non-adiabatic molecular dynamics process. As a representative example, we apply our implementation to study the ultrafast photoinduced electronic and vibrational energy transfer in a pyridine molecule. The effects of the cloning procedure on electronic and vibrational coherence, relaxation and unidirectional energy transfer are discussed. This new AIMC implementation provides a high-level non-adiabatic molecular dynamics framework for simulating photoexcited dynamics in complex molecular systems and experimentally relevant ultrafast spectroscopic probes, such as nonlinear coherent optical and X-ray signals.

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

Cities Are Concentrators of Complex, MultiSectoral Interactions Within the Human-Earth System

Cities are concentrators of complex, multi-sectoral interactions. As keystones in the interconnected human-Earth system, cities have an outsized impact on the Earth system. We describe a multi-lens framework for organizing our understanding of the complexity of urban systems and scientific research on urban systems, which may be useful for natural system scientists exploring the ways their work can be made more actionable. We then describe four critical dimensions along which improvements are needed to advance the urban research that addresses urgent climate challenges: (a) solutions-oriented research, (b) equity-centered assessments which rely on fine-scale human and ecological data, (c) co-production of knowledge, and (d) better integration of human and natural systems occurring through theory, observation, and modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Revisiting the reactivity between HCO and CH 3 on interstellar grain surfaces

The formation of interstellar complex organic molecules is currently thought to be dominated by the barrierless coupling between radicals on the interstellar icy grain surfaces. Previous standard density functional theory (DFT) results on the reactivity between CH 3 and HCO on amorphous water surfaces showed that the formation of CH 4 + CO by H transfer from HCO to CH 3 assisted by water molecules of the ice was the dominant channel. In contrast, the adopted description of the electronic structure of the biradical (i.e. CH 3 /HCO) system was inadequate [without the broken-symmetry (BS) approach]. In this work, we revisit the original results by means of BS-DFT both in gas phase and with one water molecule simulating the role of the ice. Results indicate that the adoption of BS-DFT is mandatory to describe properly biradical systems. In the presence of the single water molecule, the water-assisted H transfer exhibits a high energy barrier. In contrast, CH 3 CHO formation is found to be barrierless. Yet, direct H transfer from HCO to CH 3 to give CO and CH 4 presents a very low energy barrier, hence being a potential competitive channel to the radical coupling and indicating, moreover, that the physical insights of the original work remain valid.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Role of Electron Correlation beyond the Active Space in Achieving Quantitative Predictions of Spin-Phonon Relaxation

Single-molecule magnets (SMMs) are promising candidates for molecular-scale data storage and processing due to their strong magnetic anisotropy and long spin relaxation times. However, as temperature rises, interactions between electronic states and lattice vibrations accelerate spin relaxation, significantly limiting their practical applications. Recently, ab initio simulations have made it possible to advance our understanding of phonon-induced magnetic relaxation, but significant deviations from experiments have often been observed. The description of molecules’ electronic structure has been mostly based on complete active space self-consistent field (CASSCF) calculations, and the impact of electron correlation beyond the active space remains largely unexplored. In this study, we provide the first systematic investigation of spin-phonon relaxation in SMMs with post-CASSCF multiconfigurational methods, specifically CAS followed by second-order perturbation theory and multiconfiguration pair-density functional theory. Taking Co(II)- and Dy(III)-based SMMs as case studies, we analyze how electron correlation influences spin-phonon relaxation rates across a range of temperatures, comparing theoretical predictions with experimental observations. Our findings demonstrate that post-CASSCF treatments make it possible to achieve quantitative predictions for Co(II)-based SMMs. For Dy(III)-based systems, however, accurate predictions require consideration of additional effects, underscoring the urgent necessity of further advancing the study of the effects of electronic correlation in these complex systems.

Energy↗