Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mr”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Enhancing the mass resolving power of FRIB’s proposed high-voltage MR-ToF mass separator and spectrometer: Addressing non-ideal conditions

Multi-reflection time-of-flight mass separators and spectrometers (MR-ToF MSs) are indispensable tools at radioactive ion beam (RIB) facilities. These electrostatic ion beam traps act as highly selective mass separators and high-precision mass spectrometers for rare and exotic nuclei. When well-tuned and designed to minimize higher-order flight-time aberrations, state-of-the-art MR-ToF MSs approach, and slightly exceed, mass resolving powers of m/Δm = 10 6 . Achieving m/Δm ≥ 3 ⋅ 10 6 would provide the ability to resolve > 90 % of all known isomeric states with half-lives above 10 ms. However, the ability to mass separate in all practical setups is limited by non-ideal conditions which place such resolving powers out of reach. Here, to this end, we present a simulated analysis of these conditions in the newly proposed high-voltage MR-ToF MS for the Facility for Rare Isotope Beams (FRIB). It is expected to store ions at 30 keV beam energy and increase ion throughput by two orders of magnitude compared to current devices. Existing efforts to mitigate the effects of non-ideal conditions employed for current MR-ToF devices storing ions at < 3 keV beam energy will already enable mass resolving powers approaching 10 6 for FRIB’s high-voltage MR-ToF device. Simulations of newly proposed mitigation strategies show that even mass resolving powers approaching 10 7 might become feasible.

Electrostatic ion beam traps↗
DOE OSTI2024

ToF-SIMS spectral analysis of Shewanella oneidensis MR-1 biofilms

Analysis of bacterial biofilms is particularly challenging and important with diverse applications from systems biology to biotechnology. Among the variety of techniques that have been applied, time-of-flight secondary ion mass spectrometry (ToF-SIMS) has many powerful features in studying the surface characteristics of biofilms. ToF-SIMS offers high spatial resolution, mass resolution, and mass accuracy, which permit surface sensitive analysis of biofilm components. Thus, ToF-SIMS provides a powerful solution to addressing the challenge of bacterial biofilm analysis. This dataset covers ToF-SIMS analysis of Shewanella oneidensis MR-1 isolated from freshwater lake sediment in New York state. The MR-1 strain is known to have metal and sulfur reducing properties and it can be used for bioremediation and wastewater treatment. There is a current need to identify small molecules and fragments produced from bacterial biofilms, especially those from extracellular polymeric substance (EPS). Static ToF-SIMS spectra of MR-1 were obtained using an IONTOF TOF.SIMS V instrument equipped with a 25 keV Bi$^+_3$ metal ion gun. Identified molecules and molecular fragments are compared against known biological databases and the reported peaks have at least 65 ppm mass accuracy. These molecules range from lipids, fatty acids, flavonoids, and quinolones to other naturally occurring organic compounds. It is anticipated that the mass spectral identification of key peaks will assist detection of metabolites, EPS molecules like polysaccharides, and biologically relevant small organic molecules using ToF-SIMS in future surface and interface research.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

A high-voltage MR-ToF mass spectrometer and separator for the study of exotic isotopes at FRIB

The Facility for Rare Isotope Beams (FRIB) delivers a wide variety of rare isotopes as fast, stopped, or reaccelerated beams to enable forefront research in nuclear structure, astrophysics, and fundamental interactions. To expand the scientific potential of FRIB’s stopped and reaccelerated beam programs, we are designing a Multi-Reflection Time-of-Flight mass spectrometer and separator (MR-ToF MS). It will enable high-precision mass measurements of short-lived isotopes, improve beam diagnostics, and deliver isobarically and isomerically purified beams to downstream experimental stations. It is designed to store ions at a kinetic energy of 30 keV, significantly enhancing ion throughput while maintaining high mass resolving power. In conclusion, we present the scientific motivation, technical design, and simulations demonstrating the expected performance of the system, which has the potential to significantly enhance FRIB’s mass measurement, diagnostic, and mass separation capabilities.

Electrostatic ion beam trap↗
DOE OSTI2025

Radsource Mr: Mixed Reality Planning Tool For Radioactive Recovery

The RadSource MR system leverages Meta Quest 3's advanced mixed reality capabilities to create a comprehensive spatial planning platform for end-of-life sealed radioactive source recovery operations. The application utilizes the Quest 3's high-resolution passthrough cameras and spatial mapping algorithms to generate accurate 3D environmental models. Core technical components include: (1) Real-time spatial measurement algorithms calculating distances, angles, slopes, and surface areas with sub-centimeter accuracy; (2) Virtual object placement system allowing users to position digital representations of recovery equipment (trailers, containment vessels, protective barriers) within the real environment; (3) Voice recording and annotation system for hands-free documentation in protective equipment; (4) 3D mesh capture and storage capabilities for post-operation analysis and regulatory documentation. (5) Procedure documentation is available for viewing in Mixed Reality, providing an innovative and convenient way to access the information during pre-visit and inspection activities. (6) Support for screen capture for the view for real world and virtual objects together to use it later for planning. The system integrates computer vision techniques for environmental understanding, spatial mathematics for precise measurements, and human-computer interaction principles optimized for hazardous environment operations. Data persistence allows teams to save and share planning sessions across multiple stakeholders while maintaining operational security requirements.

Khadka, Rajiv [Idaho National Laboratory (INL), Id↗
DOE OSTI2024

SBIR Phase I Final Technical Report MODULAR NEUTRON DETECTOR ARRAY FOR PROTEIN CRYSTALLOGRAPHY BEAM LINES Award Number DE-SC0024828 Prepared for U.S. Department of Energy, Office of Science, SC-1 For the Period July 22, 2024 – April 21, 2025 Report Date November 22, 2024 Submitted by Vivek V. Nagarkar, Ph.D., Principal Investigator Report Contributors Mr. Nicholas Anastasi, Dr. Charles Sosa, Dr. Jun Wang, Dr. Lakshmi Soundara Pandian and Dr. Vivek V. Nagarkar

Explore the source record for details and available documents.

Nagarkar, Vivek V.↗
DOE OSTI2025

Fermi surface origin of the low-temperature magnetoresistance anomaly

Magnetoresistance (MR) at a fixed field can demonstrate a non-monotonic temperature dependence—an anomaly—in many systems, including low-dimensional chalcogenides, spin- and charge-density-wave metals, and topological semimetals. These systems are often low-carrier-density compensated metals, and the physics are expected to be quasi-classical. Nevertheless, the MR anomaly also exists in the highly conductive metals Cr, Mo, and W for both linear and quadratic field dependence, with their non-saturation attributed to either open orbit or electron-hole compensation. We argue that quantum transport across sharp Fermi surface arcs, but not necessarily the full cyclotron orbit, governs this MR anomaly, thereby accounting for the profound effects of disorder. In Cr, an overlay exists between three temperature dependences: MR at a constant high field, linear MR at a low field, and Shubnikov-de Haas (SdH) oscillations of the smallest orbit. In Mo, the MR anomaly extends beyond the temperature of its SdH oscillations but disappears before Kohler’s scaling reemerges.

Feng, Yejun [Okinawa Institute of Science and Tech↗
DOE OSTI2025

Assessing fuel cycle design options for advanced microreactors

Microreactor (MR) technology have recently gained traction due to their reduced construction cost compared to traditional large reactors. Their small physical geometry results in increased neutron leakage and deteriorated neutron economy which limits the achievable burnup of MR fuel. Refueling strategy currently discussed in MR design community is the battery-like refueling which is involve replacing the entire reactor core at the end of their operational cycle. Here, this paper investigates partial core refueling as means to increase MR fuel burnup and improve fuel cycle economy as a result. Two fuel cycle design strategies have been successful. These are the partial core refueling and partial core refueling combined with assembly rotation. Both strategies have been tested to operate a heatpipe-cooled MR referred to as the eVinci-like core for 72 Effective Full Power Months (EFPMs). The partial core refueling showed to reduce the refueling requirements by 18% compared to the battery-like refueling strategy. This corresponds to 10.5% reduction in the Net Present Value (NPV) of fuel cost over the 72 EFPMs of operation. The partial core refueling combined with rotation case reduces the refueling requirements by 29% (corresponding to 14.5% reduction in NPV of fuel cost) compared to the battery-like refueling strategy.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Feasibility of Recycling Discharged Microreactor Heavy Metal in Light Water and Sodium-Cooled Fast Reactors: A Neutronics Analysis

Nuclear microreactors (MRs) offer unique advantages, such as rapid deployment, potability, low maintenance requirements, and operational flexibility. Their compact size makes them a promising solution for decentralized power generation, particularly in remote areas, military bases, and disaster-stricken regions. However, MRs face challenges, including unutilized fissile material at the end of life, economic inefficiency, increased heavy metal (HM) waste complicating disposal, and the accumulation of plutonium (Pu) with high 239 Pu concentrations raising proliferation risks. Here, this study investigated the neutronics feasibility of a novel three-stage fuel cycle where discharged HM from MRs is recycled and burned in light water reactors and sodium-cooled fast reactors. This approach converts discharged HM into valuable fuel, enhancing the efficiency of MR deployments while improving the safeguardability of their final waste products. Neutronics analysis demonstrated that the safety characteristics of reactor designs in each stage were minimally impacted by the proposed cycle. For two representative MR designs, a fast-spectrum MR with solid pellet fuel and a thermal-spectrum MR with TRISO (TRi-structural-ISOtropic) fuel compacts, the proposed fuel cycle reduced the uranium disposal mass flow rate by ~60%, decreased the 235 U enrichment of the discharge fuel to ~1 wt%, eliminated plutonium disposal, and increased the cumulative fuel burnup to ~580 gigawatt-day per metric ton of initial heavy metal (GWd/t-iHM) or 60% fissions per initial metal atom. Despite the significant differences between the two MR designs, the performance and infrastructure requirements of the developed fuel cycles were remarkably similar, indicating its generalizability to a broader class of MRs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Multireference Equation-of-Motion Driven Similarity Renormalization Group: Theoretical Foundations and Applications to Ionized States

We present a formulation and implementation of an equation-of-motion (EOM) extension of the multireference driven similarity renormalization group (MR-DSRG) formalism for ionization potentials (IP-EOM-DSRG). The IP-EOM-DSRG formalism results in a Hermitian generalized eigenvalue problem, delivering accurate ionization potentials for strongly correlated systems. The EOM step scales as O(N 5 ) with the basis set size N, allowing for efficient calculation of spectroscopic properties, such as transition energies and intensities. The IP-EOM-DSRG formalism is combined with three truncation schemes of the parent MR-DSRG theory: an iterative nonperturbative method with up to two-body excitations [MR-LDSRG(2)] and second- and third-order perturbative approximations [DSRG-MRPT2/3]. We benchmark these variants by computing (1) the vertical valence ionization potentials of a series of small molecules at both equilibrium and stretched geometries; (2) the spectroscopic constants of several low-lying electronic states of the OH, CN, N 2 + , and CO + radicals; and (3) the binding curves of low-lying electronic states of the CN radical. A comparison with experimental data and theoretical results shows that all three IP-EOM-DSRG methods accurately reproduce the vertical ionization potentials and spectroscopic constants of these systems. Notably, the DSRG-MRPT3 and MR-LDSRG(2) versions outperform several state-of-the-art multireference methods of comparable or higher cost.

Hamiltonians↗
DOE OSTI2026

Pressure-Induced Metal-like Transport and Magnetoresistance in a Au 2+ –Au 3+ Halide Perovskite

The Cs 4 Au II Au III 2 Cl 12 perovskite (1), featuring AuCl 4 trimers separated by vacancies, enables the first high-pressure study of Au 2+/3+ mixed-valence. Our computational analysis of the gold frontier orbitals suggests that the Au 2+ →Au 3+ intervalence charge transfer (IVCT) occurs across the vacancies. Computational structures indicate that these vacancies rapidly shrink with pressure and the Au 2+ and Au 3+ coordination spheres become very similar at the phase transition to nearly cubic symmetry at ca. 15 GPa─enabling facile IVCT. Although the activation energy of conductivity of 0.73(4) meV and far-infrared absorption indicate a small but nonzero bandgap, ambient thermal energy drives the IVCT, affording metallic properties: prominent infrared reflectivity and transport values of 10 2 S·cm –1 . This prompted us to perform the first high-pressure studies of magnetoresistance (MR) and Hall effect in halide perovskites. At 16 GPa, the MR increases by 9.3% at 2 K and 9 T; this value is maintained up to 27 GPa, when a local distortion drives electronic localization. By globally fitting the MR and Hall resistance to a two-carrier model we quantify how the carrier densities and mobilities evolve with pressure. Thus, metal-like transport and MR in 1 is driven by a pressure-induced transition from localized to partially delocalized mixed-valence.

Deschene, Christina R. [Stanford University, CA (U↗
DOE OSTI2025

Dual roles of stepped surfaces: Catalytic initiators and stabilizers of oxygen-induced reconstructions

Understanding surface restructuring under reactive conditions is crucial for designing next-generation catalysts with enhanced activity and selectivity. Here, we employ in situ transmission electron microscopy to directly observe the dynamic behavior of Cu(100) and Cu(410) surfaces under both oxidizing and vacuum annealing conditions, revealing a complex interplay among surface crystallography, local oxygen coverage, Cu atom mobility, and step-edge reactivity. The stepped Cu(410) surface acts as an active site for O 2 dissociation, triggering the oscillatory transformation of the c(2 × 2)–O phase into the more stable (2$\sqrt2$ ×$\sqrt2$)R45°–O missing-row (MR) structure on the adjacent flat Cu(100) terrace. Under subsequent vacuum annealing, this same Cu(410) facet exhibits remarkable structural resilience, preserving the MR reconstruction and chemisorbed oxygen. In contrast, the Cu(100) surface undergoes reversible transitions from the MR structure back to the c(2 × 2)–O phase. These results highlight the critical role of surface morphology in directing both the formation and stability of oxygen-induced reconstructions, demonstrating that stepped surfaces serve dual roles as both catalytic initiators and structural stabilizers. Furthermore, this work offers atomic-level insights into the environment-responsive behavior of copper surfaces, establishing a mechanistic basis for designing Cu-based catalysts through facet-specific control of surface reactivity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling Oxidative Dehydrogenation of Propane with Supported Vanadia Catalysts Using Multireference Methods

The oxidative dehydrogenation of propane over supported vanadium oxide catalysts poses significant computational challenges due to complex electronic structure changes along the reaction coordinate, driven primarily by changes in the oxidation states of vanadium. To address these challenges, we systematically test quantum chemical methods, including multireference (MR) approaches, domain-based local pair natural orbital coupled cluster theory (DLPNO-CCSD(T)), and density functional theory (DFT). The initial C–H bond-breaking transition state requires MR treatment due to its multireference character, while subsequent steps permit efficient single-reference calculations. For the rate-limiting C–H activation step mediated by the vanadyl moiety, complete 1 active space second-order perturbation theory (CASPT2) yields an apparent activation barrier (E app 600K ) of 138 kJ/mol, consistent with experimental values (134 ± 4 kJ/mol; Gruene et al. Catal. Today 2010, 157, 137). In contrast, DLPNO-CCSD(T) overestimates this barrier (198 kJ/mol), whereas DFT predictions span 125–150 kJ/mol, depending on the functional. Our multireference investigation of this transition metal oxide-catalyzed process demonstrates that an active space that incorporates the C–H σ and V=O σ/π bonding orbitals, oxygen lone pairs, and their antibonding counterparts adequately captures electronic structure changes along the chemical transformation. Furthermore, these findings provide a general strategy for active space selection in transition metal oxide-catalyzed C/O–H bond activation reactions. The reference dataset from this work, which includes MR calculations with manually selected active spaces for all intermediates and transition states in the propane ODH reaction network, will serve as a benchmark for automating active space selection in similar systems.

Catalysts↗
DOE OSTI2026

Large magnetoresistance and quantum oscillations observed in superconducting β -IrSn 4

Diverse quantum phenomena have been observed in TSn 4 transition metal (T) stannides, including superconductivity, nontrivial topology, and large magnetoresistance (MR) at low temperatures. Here, we report the experimental and theoretical investigation of tetragonal β-IrSn 4 (space group I 4 1 /acd) with the lattice parameters a = 6.362(3) Å and c = 22.723(0) Å. The temperature dependence of the electrical resistivity indicates that β-IrSn 4 is a good metal with Fermi-liquid behavior above the superconducting transition at T c ∼ 0.7 K. Magnetic susceptibility measurements indicate it is non-magnetic. However, positive transverse MR with linear field dependence is observed. At T= 2 K and Μ 0 H = 14 T, the MR reaches 700% without sign of saturation. Shubnikov–de Haas oscillations are observed from proximity detector oscillator measurements up to 60 T. We experimentally extract several frequencies, and through comparison to first-principles calculations, identify corresponding electronic bands. For H // c, two strong oscillations, F 1 ≈ 109 and F 2 ≈ 302 T, allow us to determine effective masses 0.266m 0 and 0.300m 0 (m 0 is the free-electron mass), respectively. The topology of these bands is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Microreactor Assembly Transportation Cask Model Description for Criticality Safety Validation Basis Assessment

Criticality safety analyses are completed on a transportation cask used for microreactor assembly shipment to provide an example of model and analysis to industry for reproducing this type of study on their microreactor fuel shipment. The fuel assembly considered is based on a gas-cooled microreactor (GC-MR), which utilizes HALEU fuel in the form of TRISO particles and utilizes various design options considered in industry designs. Various versions of this GC-MR assembly were studied, with and without YH2 moderator, providing similar conclusions. The shipment cask design is revised based on an existing design ES-3100, developed by Y-12 for the transport of highly enriched uranium (HEU), but is enlarged to hold the GC-MR fuel assembly. Criticality safety analysis for the cask/GC-MR fuel assembly package was performed using the CSAS6 sequence of SCALE6.3.2, utilizing the ENDF/B-VII.1 based continuous energy neutron library, and the analysis strictly follows the guideline from NRC reference reports. Different scenarios, e.g. normal operation, undamaged cask with water flooded, damaged cask with optimal water moderation, have been analyzed and it could be concluded the package would always have a large margin of subcriticality even packed in an infinite array. Sensitivity and similarity analyses are also performed using the TSUNAMI sequence of SCALE6.3.2, and the similarity analysis uses all the experiments from the ICSBEP Handbook with Intermediate and Mixed Enriched Uranium (IEU) and Low Enriched Uranium (LEU) systems together with additional ones that are sponsored by the DNCSH program. These similarity analyses indicate that dry cases have no similar benchmark experiments (ck values greater than 0.8), which may become problematic if more assemblies are shipped together (or a fully loaded core is shipped) and margin to criticality is reduced. However, the damaged cask models with flooded assemblies exhibited similarities to many experiments with ck values greater than 0.8.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Microreactor Assembly Transportation Cask Model Description for Criticality Safety Validation Basis Assessment (Rev. 3)

Criticality safety analyses are completed on a transportation cask used for microreactor assembly shipment to provide an example of model and analysis to industry for reproducing this type of study on their microreactor fuel shipment. The fuel assembly considered is based on a gas-cooled microreactor (GC-MR), which utilizes HALEU fuel in the form of TRISO particles and utilizes various design options considered in industry designs. Various versions of this GC-MR assembly were studied, with and without YH 2 moderator, providing similar conclusions. The shipment cask design is revised based on an existing ES-3100 design, developed by Y-12 for the transport of highly enriched uranium (HEU), but is enlarged to hold the GC-MR fuel assembly. Criticality safety analysis for the cask/GC-MR fuel assembly package was performed using the CSAS6 sequence of SCALE6.3.2, utilizing the ENDF/B-VII.1 based continuous energy neutron library, and the analysis strictly follows the guideline from NRC reference reports. Different scenarios, e.g. normal operation, undamaged cask with water flooded, damaged cask with optimal water moderation, have been analyzed and it could be concluded the package would always have a large margin of subcriticality even packed in an infinite array. Sensitivity and similarity analyses are also performed using the TSUNAMI sequence of SCALE6.3.2, and the similarity analysis uses all the experiments from the ICSBEP Handbook with Highly Enriched Uranium (HEU), Intermediate and Mixed Enriched Uranium (IEU) and Low Enriched Uranium (LEU) systems, together with additional ones that are sponsored by the DNCSH program, and selected IRPhEP experiments using TRISO fuel and graphite moderator. These similarity analyses indicate that the dry nominal design has no similar benchmark experiments (ck values greater than 0.8), which may become problematic if more assemblies are shipped together (or a fully loaded core is shipped) and margin to criticality is reduced. However, the cask models with flooded assemblies exhibited similarities to many experiments with c k values greater than 0.8.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

The Velocity Field of Our Milky Way Outer Stellar Halo Based on DESI DR2

Using 64,000 halo K giants from the Dark Energy Spectroscopic Instrument second Data Release (DR2), we decompose the Milky Way (MW) stellar halo between 3 and 160 kpc into metal-rich (MR) and metal-poor (MP) components via a Gaussian mixture model. The two populations are nearly equal in number but chemically and kinematically distinct: MR stars occupy highly radial orbits with velocity anisotropy of β ≈ 0.94 and metallicity dispersion σ [Fe/H] ≈ 0.17 dex, without obvious dependence on distance, and are mainly contributed by Gaia-Sausage/Enceladus (GSE) debris. The MR component dominates the inner 30 kpc and reemerges beyond 50 kpc, implying GSE debris can extend to ∼70–80 kpc. MP stars exhibit a weaker radial bias of β ≈ 0.46, decreasing to −0.5 beyond 80 kpc, and with a larger metallicity dispersion of σ [Fe/H] ≈ 0.46 dex, showing signatures of multiple minor mergers. Both components exhibit net prograde rotation at ∼10–30 kpc with a stronger azimuthal signal in the MP population. The nonequilibrium motions of the outer halo (>50 kpc) are quantified with a dipole-plus-contraction velocity field. We find that the outer halo is simultaneously contracting (μ compr = −19 km s −1 , distance-independent) and subject to reflex motions (μ dipole increases from −19 to −44 km s −1 with radius), reflecting the perturbation from the Large Magellanic Cloud (LMC). We also confirm a linear dependence of mean polar velocity for the outer stellar halo on μ dipole , a direct consequence of the LMC and MW interaction. Our results provide a quantitative distance-resolved description of the MW’s last major accretion event and its ongoing response to the first infall of the LMC.

Li, Songting [Shanghai Jiao Tong University; Shang↗
DOE OSTI2026

Low-Lying Excited States of Linear All- Trans Polyenes: Insights from Analytic Gradient and Nonadiabatic Coupling Calculations Based on Multireference Configuration Interaction

Polyenes serve as a rigorous test for theoretical models and electronic structure methods, playing a key role in advancing computational and theoretical chemistry. Here, we present a high-level theoretical investigation of linear, all-trans polyenes using energy gradients and nonadiabatic coupling vectors based on an MR-CISD wave function to describe electronic transitions involving the ground state (1 1 A g – ) and three low-lying excited states (2 1 A g – , 1 1 B u + , and 2 1 B u – ) of hexatriene, octatetraene, and decapentaene. This approach enables accurate evaluation of both adiabatic and vertical excitation and emission energies, yielding results in excellent agreement with experiment, as well as locating minima on the crossing seam between adiabatic states. Our results show that vertical excitation energies to the 1 1 B u + state are blue-shifted by 0.2–0.3 eV relative to the experimental absorption maximum, whereas the vertical emission energy from the 2 1 A g – state is red-shifted by ∼0.2 eV relative to the experimental emission maximum. Upon relaxation from the Franck–Condon geometry, the 2 1 A g – state stabilizes by around 1 eV, compared to 0.2–0.3 eV for the 1 1 B u + state. An analysis of the S 1 /S 0 crossing seam in hexatriene shows that its minimum involves asymmetric backbone deformations and provides an efficient channel for ultrafast internal conversion to the ground state, consistent with the absence of detectable fluorescence in this molecule. These results demonstrate the power of analytic gradients and nonadiabatic coupling vectors based on an MR-CISD wave function for accurately characterizing the electronic structure and photophysics of polyenes.

Excited states↗
DOE OSTI2025

COLUMBUS─An Efficient and General Program Package for Ground and Excited State Computations Including Spin–Orbit Couplings and Dynamics

The COLUMBUS program system provides the tools for performing high-level multireference (MR) computations, including the multireference configuration interaction (MRCI) method and its multireference averaged quadratic coupled cluster (MR-AQCC) extension, allowing computations on a wide range of fascinating atomic and molecular systems, including the treatment of open-shells and complicated excited state phenomena. The inclusion of spin−orbit coupling (SOC) directly within the MRCI step enables the description of systems containing heavy elements, such as lanthanides and actinides, whose properties are strongly influenced by SOC. Analytic energy gradients and nonadiabatic couplings at the correlated MRCI level provide the foundation for a variety of dynamics studies, giving insight into ultrafast photochemistry. New and ongoing method developments in COLUMBUS include the computation of spin densities, improved descriptions of ionic states, enhancements to the AQCC method, and the porting of COLUMBUS to graphical processing units (GPUs). New external interfaces enable an enhanced description of electronic resonances and molecules in strong laser fields. This work highlights these new developments while providing a detailed account of the diverse applications of COLUMBUS in recent years.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗