Search NASASearch

SEARCH · Search NASA

Results for “quantum phase estimation”

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.

70 records · Page 4

High-performance silicon photonic single-sideband modulators for cold-atom interferometry

The laser system is the most complex component of a light-pulse atom interferometer (LPAI), controlling frequencies and intensities of multiple laser beams to configure quantum gravity and inertial sensors. Its main functions include cold-atom generation, state preparation, state-selective detection, and generating a coherent two-photon process for the light-pulse sequence. To achieve substantial miniaturization and ruggedization, we integrate key laser system functions onto a photonic integrated circuit. Our study focuses on a high-performance silicon photonic suppressed-carrier single-sideband (SC-SSB) modulator at 1560 nanometers, capable of dynamic frequency shifting within the LPAI. By independently controlling radio frequency (RF) channels, we achieve 30-decibel carrier suppression and unprecedented 47.8-decibel sideband suppression at peak conversion efficiency of –6.846 decibels (20.7%). We investigate imbalances in both amplitudes and phases between the RF signals. Using this modulator, we demonstrate cold-atom generation, state-selective detection, and atom interferometer fringes to estimate gravitational acceleration, g ≈ 9.77 ± 0.01 meters per second squared, in a rubidium ( 87 Rb) atom system.

47 OTHER INSTRUMENTATION

Observation of the axion quasiparticle in 2D MnBi 2 Te 4

The axion is a hypothetical fundamental particle that is conjectured to correspond to the coherent oscillation of the θ field in quantum chromodynamics. Its existence would solve multiple fundamental questions, including the strong CP problem of quantum chromodynamics and dark matter, but the axion has never been detected. Electrodynamics of condensed-matter systems can also give rise to a similar θ, so far studied as a static, quantized value to characterize the topology of materials. Coherent oscillation of θ in condensed matter has been proposed to lead to physics directly analogous to the high-energy axion particle—the dynamical axion quasiparticle (DAQ). Here we report the observation of the DAQ in MnBi 2 Te 4 . By combining a two-dimensional electronic device with ultrafast pump–probe optics, we observe a coherent oscillation of θ at about 44 gigahertz, which is uniquely induced by its out-of-phase antiferromagnetic magnon. This represents direct evidence for the presence of the DAQ, which in two-dimensional MnBi 2 Te 4 is found to arise from the magnon-induced coherent modulation of the Berry curvature. The DAQ also has implications in light–matter interaction and coherent antiferromagnetic spintronics, as it might lead to axion polaritons and electric control of ultrafast spin polarization. Finally, the DAQ could be used to detect axion particles. We estimate the detection frequency range and sensitivity in the millielectronvolt regime, which has so far been poorly explored.

topological matter

Signs of nonmonotonic finite-volume corrections to 𝑔 𝐴

We study finite-volume (FV) corrections to determinations of 𝑔 𝐴 via lattice quantum chromodynamics (QCD) using analytic results and numerical analysis. We observe that 𝑆⁢𝑈⁡(2) heavy Baryon chiral perturbation theory does not provide an unambiguous prediction for the sign of the FV correction, which is not surprising when one also considers large-𝑁 𝑐 constraints on the axial couplings. We further show that nonmonotonic FV corrections are naturally allowed when one considers either including explicit Δ-resonance degrees of freedom or one works to higher orders in the chiral expansion. We investigate the potential impact of these FV corrections with a precision study of 𝑔 𝐴 using models of FV corrections that are monotonic and nonmonotonic. Using lattice QCD data that is approximately at the 1% level of precision, we do not see significant evidence of nonmonotonic corrections. Looking forward to the next phase of lattice QCD calculations, we estimate that calculations that are between the 0.1% and 1% level of precision may be sensitive to these FV artifacts. Finally, we present an update of the CalLat prediction of 𝑔 𝐴 in the isospin limit with subpercent precision, 𝑔$^{QCD}_{𝐴}$ = 1.2674⁢(96).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Significance of the structural configuration of B2 disorder in Co and Ti based Heusler alloys

We investigate here structural (at local and global levels) and transport properties for 𝑋 2 ⁢MnAl (𝑋= Co and Ti). Additionally, the magnetic properties were also studied for Ti 2 ⁢MnAl. Our x-ray diffraction results show that both the compounds stabilize in B2 disordered phase with cubic structure of 𝑃⁢𝑚⁢$\overline{3}$⁢𝑚 space group. Further, the structural configuration of the above disordered phase for both the compounds was identified using combined studies of x-ray absorption spectroscopy and multiple scattering calculations at the transition metal 𝐾 edges. Upon such identification, in the case of Co 2 ⁢Mn 1−𝑦⁢ Cr 𝑦 ⁢Al (𝑦= 0, 0.05, 0.1, 0.2) with change in 𝑦, we are able to establish a better connection quantitatively between the inverse of Mn-Co bonds and peak in the temperature-dependent resistivity. This highlights the crucial importance of a detailed understanding of the nature of B2 disorder. In the case of 𝑦=0, in the temperature range of study, the resistivity is driven by the functional form associated with (a) three-dimensional enhanced electron-electron Coulomb interaction scattering mechanism and (b) an unconventional one-magnon process. For Ti 2 ⁢MnAl, the transport shows metallic glasslike behavior at high temperature, while at low temperature it follows both the Cote-Meisel's model and quantum correction model. In this compound, the magnetic studies suggest the formation of superparamagnetic clusters in the paramagnetic matrix at low temperatures. Our density functional theory results are in line with the transport and magnetic properties. In literature, the spin polarization percentage (𝑃) for 𝑋= Co in B2 disordered phase is 76%. However, the present results emphasize the fact that in B2 disordered phase, the value of 𝑃 can range from 90% to 71% depending on the structural configuration introduced by swapping of the atomic positions of Mn and Al. For the compound under study, the value of percentage spin polarization obtained ranges between 82% to 85%. In addition, we also identify the origin of the difference in the shape of the Mn 3⁢𝑑 density of states for both the alloys. In conclusion, our results for 𝑋= Co alloy highlights the importance of identifying the specific structural configuration associated with a particular disorder category especially in the estimation of 𝑇 𝑐 and 𝑃 and for 𝑋= Ti, the physical properties can be tuned by varying the position of 𝐸 𝐹 and thereby its utilization in device applications.

36 MATERIALS SCIENCE

Random insights into the complexity of two-dimensional tensor network calculations

Projected entangled pair states (PEPS) offer memory-efficient representations of some quantum many-body states that obey an entanglement area law and are the basis for classical simulations of ground states in two-dimensional (2d) condensed matter systems. However, rigorous results show that exactly computing observables from a 2d PEPS state is generically a computationally hard problem. Yet approximation schemes for computing properties of 2d PEPS are regularly used, and empirically seen to succeed, for a large subclass of (“not too entangled”) condensed matter ground states. Adopting the philosophy of random matrix theory, in this work, we analyze the complexity of approximately contracting a 2d random PEPS by exploiting an analytic mapping to an effective replicated statistical mechanics model that permits a controlled analysis at a large bond dimension. Through this statistical-mechanics lens, we argue that (i) although approximately sampling wave-function amplitudes of random PEPS faces a computational-complexity phase transition above a critical bond dimension, and (ii) one can generically efficiently estimate the norm and correlation functions for any finite bond dimension. Furthermore, these results are supported numerically for various bond-dimension regimes. It is an important open question whether the above results for random PEPS apply more generally also to PEPS representing physically relevant ground states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Association Kinetics for Perfluorinated n -Alkyl Radicals

Radical-radical reaction channels are important in the pyrolysis and oxidation chemistry of perfluoroalkyl substances (PFAS). In particular, unimolecular dissociation reactions within unbranched n-perfluoroalkyl chains, and their corresponding reverse barrierless association reactions, are expected to be significant contributors to the gas-phase thermal decomposition of families of species such as perfluorinated carboxylic acids and perfluorinated sulfonic acids. Unfortunately, experimental data for these reactions are scarce and uncertain. Furthermore, obtaining reliable theoretical predictions for such reactions is a laborious and computationally intensive task. Here, in this work, the chemical kinetics of the various association/decomposition reactions producing/decomposing the C 2 -C 4 series of unbranched n-perfluoroalkanes (C 2 F 6 , C 3 F 8 , and C 4 F 10 ) are examined using state-of-the-art ab initio transition-state-theory-based master-equation calculations. The variable-reaction-coordinate transition-state theory (VRC-TST) formalism is employed in computing the microcanonical and canonical rates for the association reactions. Reaction thermochemistry is obtained via composite quantum chemistry calculations and the laddering of error-canceling reaction schemes via a connectivity-based hierarchy approach employing ANL1/ANL0-style reference energies. Lennard-Jones collision model parameters for the considered systems were estimated by a direct dynamics approach, and collisional energy transfer parameters were obtained from analogies to systems of similar size and heavy-atom connectivity. A one-dimensional master equation approach was used to convert the microcanonical rate coefficients from the VRC-TST analysis into temperature- and pressure-dependent rate constants for the association reactions and the reverse dissociation reactions. The data are reported in standardized formats for usage in comprehensive chemical kinetic models for PFAS thermal destruction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Nondegenerate Noise-Resilient Superconducting Qubit

We propose a superconducting qubit based on engineering the first and second harmonics of the Josephson energy and phase relation 𝐸 𝐽⁢1 ⁢cos 𝜑 and 𝐸 𝐽⁢2 ⁢cos 2⁢𝜑. By constructing a circuit such that 𝐸 𝐽⁢2 is negative and |𝐸 𝐽⁢1 | ≪ |𝐸 𝐽⁢2 |, we create a periodic potential with two nondegenerate minima. The qubit, which we dub “harmonium,” is formed from the lowest-energy states of each minimum. Bit-flip protection of the qubit arises due to the localization of each qubit state to their respective minima, while phase-flip protection can be understood by considering the circuit within the Born-Oppenheimer approximation. We demonstrate with time-domain simulations that single- and two-qubit gates can be performed in approximately 100 ns. Finally, we compute the qubit coherence times using numerical diagonalization of the complete circuit in conjunction with state-of-the-art noise models. We estimate out-of-manifold heating times on the order of milliseconds, which can be treated as erasure errors using conventional dispersive readout. We estimate pure-dephasing times on the order of many tens of milliseconds, and bit-flip times on the order of seconds.

quantum error correction

Lower bounds on entanglement entropy without twin copy

We discuss the possibility of estimating experimentally the von Neumann entanglement entropy S A v N of a symmetric bipartite quantum system A B by using the basic measurement counts (bitstrings) for a single copy of a prepared state. Using exact diagonalization and analog simulations performed with the publicly available QuEra facilities for chains and ladders of Rydberg atoms, we calculate the Shannon entropy S A B X associated with the bitstrings of adiabatically prepared ground states and the reduced entropies S A X and S B X obtained from the marginal probabilities in A and B . We then calculate the classical mutual information I A B X = S A X + S B X − S A B X , which is a lower bound on S A v N . We show that for a broad range of lattice spacing and detuning, I A B X is typically 20% below S A v N in regions where S A v N is large and a less close bound in regions where S A v N is low. We argue that this use of the easily available bitstrings provides a robust and efficient way to explore empirically the phase diagram of qubit-based quantum simulators and identify critical regions. Published by the American Physical Society 2025

Meurice, Yannick (ORCID:0000000209959694)

Physical-mass calculation of ρ ( 770 ) and K * ( 892 ) resonance parameters via π π and K π scattering amplitudes from lattice QCD

We present our study of the ρ ( 770 ) and K * ( 892 ) resonances from lattice quantum chromodynamics (QCD) employing domain-wall fermions at physical quark masses. We determine the finite-volume energy spectrum in various momentum frames and obtain phase-shift parametrizations via the Lüscher formalism and as a final step the complex resonance poles of the π π and K π elastic scattering amplitudes via an analytical continuation of the models. By sampling a large number of representative sets of underlying energy-level fits, we also assign a systematic uncertainty to our final results. This is a significant extension to data-driven analysis methods that have been used in lattice QCD to date, due to the two-step nature of the formalism. Our final pole positions, M + i Γ / 2 , with all statistical and systematic errors exposed, are M K * = 893 ( 2 ) ( 8 ) ( 54 ) ( 2 ) MeV and Γ K * = 51 ( 2 ) ( 11 ) ( 3 ) ( 0 ) MeV for the K * ( 892 ) resonance and M ρ = 796 ( 5 ) ( 15 ) ( 48 ) ( 2 ) MeV and Γ ρ = 192 ( 10 ) ( 28 ) ( 12 ) ( 0 ) MeV for the ρ ( 770 ) resonance. The four differently grouped sources of uncertainties are, in the order of occurrence: statistical, data-driven systematic, an estimation of systematic effects beyond our computation (dominated by the fact that we employ a single lattice spacing), and the error from the scale-setting uncertainty on our ensemble. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Numerical simulations of liquid jetting with solid inclusions

The dynamics of finite-sized particles in fluids, and their influence on the overall flow, are of great interest across several industrial, environmental, and medical fields. In the context of inkjet printing, the presence of solid inclusions can be either intentional, as in additive manufacturing, or unintentional, as in standard printing processes. These inclusions can strongly impact the jetting process, causing effects such as jet asymmetry, bubble entrapment, and the formation of satellite droplets. Understanding and controlling particle behavior is therefore essential, particularly to predict how and when particles are ejected over multiple jetting cycles. It is therefore critical to develop reliable models that allow for a deeper understanding of the complex interplay between particle and fluid during the whole printing process. To address this, we present a tailored implementation of the Color-Gradient multicomponent Lattice Boltzmann Method for fully resolved three-dimensional (3D) simulations of multicycle liquid jetting with particles. Our method supports realistic parameter settings aligned with industrial inkjet systems, and we provide both qualitative and quantitative validation against experimental data. Additionally, we introduce a simplified model based on the Stokes drag law, in which solid particles are represented as point particles and do not influence the fluid flow. Despite this limitation, the model offers a computationally efficient means to explore the vast parameter space typically encountered in industrial applications, allowing, e.g., identifying critical ejection regions and estimating the number of cycles required for particle release. These qualitative insights are valuable for guiding and complement fully two-way coupled simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC

Manipulating Cu vacancy in Yb triangular lattice Materials: Investigating the quantum disordered compound CuYbSe 2 and spin glass Cu 2.73 ⁢YbP 2

Yb-based triangular lattice materials have attracted significant attention due to their potential for hosting exotic magnetic quantum states. In this paper, we report the structural, magnetic, and thermodynamic properties of two Yb triangular lattice compounds with vacancy disorder at nonmagnetic Cu sites: CuYbSe 2 and Cu 2.73 ⁢YbP 2 . Single-crystal x-ray diffraction determines notable Cu-site vacancies in both materials. Unlike a single Cu site with 0.5 occupancy in CuYbSe 2 , Cu 2.73 ⁢YbP 2 only presents 27% vacancies on one of three Cu sites. Magnetic susceptibility measurements indicate dominant antiferromagnetic interactions between Yb 3+ ions. The effective magnetic moment is consistent with an effective spin-1/2 state of Yb 3+ . CuYbSe 2 exhibits strong magnetic anisotropy between in-plane and out-of-plane directions with the estimated exchange interactions 𝐽 𝑎 /𝑘 B = 6.48 K and 𝐽 𝑐 /𝑘 B = 2.55 K. No magnetic ordering is observed down to 0.4 K, supporting its quantum disordered ground state and candidacy as a quantum spin liquid. In contrast, Cu 2.73 ⁢YbP 2 shows clear evidence of spin freezing at 4.5 K in both magnetization and heat capacity data, indicating a spin glass state. The field-induced magnetic orders are found in CuYbSe 2 . These results highlight the critical role of structural vacancies in perturbing Yb-based triangular lattices. In conclusion, we demonstrate that while the intrinsic Yb triangular lattices are comparable, the presence of Cu vacancies can fundamentally alter the magnetic ground state—tuning the system between quantum disordered and spin glass behavior.

Frustrated magnetism

Photon-odderon interference in exclusive χ c charmonium production at the Electron-Ion Collider

Exclusive C = + 1 scalar, axial-vector, and tensor quarkonium production in high-energy electron-proton scattering requires a C -odd t -channel exchange of a photon or a three-gluon ladder. We derive the expressions for the corresponding amplitudes. The relative phase of the photon vs three-gluon exchange amplitudes is determined by the sign of the light-front matrix element of the eikonal color current operator d a b c J + a J + b J + c at moderate x , and is not affected by small- x QCD evolution. Model calculations predict constructive interference, which is particularly strong for momentum transfer | t | ∼ 1 GeV 2 where the cross section for χ c J production exceeds that for pure photon exchange by up to a factor of 4. Exclusive χ c J electroproduction at the high luminosity Electron-Ion Collider occurs with well measurable rates and measurements of these processes should find evidence for the perturbative odderon exchange. We also compute the total electroproduction cross section as a function of energy and provide first estimates of the number of χ c J events per month at the Electron-Ion Collider design luminosity. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Quantum-enhanced distributed phase sensing with a truncated SU(1,1) interferometer

In recent years, distributed quantum sensing has gained interest for a range of applications requiring networks of sensors, from global-scale clock synchronization to high energy physics. In particular, a network of entangled sensors can improve not only the sensitivity beyond the shot noise limit, but also enable a Heisenberg scaling with the number of sensors. Here, using bright entangled twin beams, we theoretically and experimentally demonstrate the detection of a linear combination of two distributed phases beyond the shot noise limit with a truncated SU(1,1) interferometer. Specifically, we show a quantum noise reduction of 1.7 ± 0.3 dB below what is possible with the corresponding classical configuration. Additionally, we theoretically extend the use of a truncated SU(1,1) interferometer to a multi-phase-distributed sensing scheme that leverages entanglement as a resource to achieve a quantum improvement in the scaling with the number of sensors in the network. Our results pave the way for developing quantum-enhanced sensor networks that can achieve an entanglement-enhanced sensitivity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC

Bridging the time scale in exascale computing of chemical systems (Final Technical Report)

This report summarizes the work carried out with support of the United States Department of Energy under Award DE-SC0019441. The theme of this project was to develop and apply methods that allowed for the acceleration of atomistic calculations, particularly in challenging areas such as multiphase systems, electrified interfaces, uncertainty estimation, and applications requiring chemical accuracy, which tend to be applications where simulation time is severely bottlenecked by the computational time requirements. Much of the focus was on the application of emerging machine-learning methodologies, although a wide range of methodologies were employed. This report has two major sections. The first focuses on the methodological advances themselves. Within this part, we report a number of major advances, a few examples of which are described here. We report the first machine-learning scheme for the acceleration of electronically grand-canonical calculations (that is, those applicable to electrochemistry). We report new methods of performing transfer learning, in which physics-based priors can be used to provide predictions, often with uncertainty estimates, of images well outside of training sets; we also offer ways to fine-tune these transfer-learning models. We provide a new systematic means to generate and apply minimal training data sets to very large (10,000’s of atoms) systems, with only small training sets appropriate for electronic structure. We developed new methodologies to integrate surface vibrations into surface adsorption calculations. We made advances to the applicability of diffusion Monte Carlo methods to allow (learned) force prediction, finite-size error correction, and force-free means of searching for transition states. We integrated machine-learned atomistic predictions into mechanism generation codes. Additionally, we released new software including AmpTorch, a modernized version of our original atomistic machine-learning code Amp. The second part of this report focuses on the scientific applications that accompanied, and were often enabled by, the methodological advances described earlier. A few examples follow, but full details are in the individual chapters of the report. For example, we developed a general theory of phonon-induced friction on molecular adsorbates. We showed fundamentally how solvent influences the adsorption and desorption process and how it differs from the processes typically involved at the solid–gas interface, making aqueous-phase and electrocatalysis different from traditional thermocatalysis. We examined how metal–insulator and magnetic transitions can be probed, and accelerated exciton dynamics via Frenkel Hamiltonian parameters. We showed that the nearsighted force-training approach, developed within this project, can predict both the stability and reactivity of large nanoparticles, and can also lead to insights on catalyst coverage on binding energies and entropies. These applied studies, which generally integrated with our method development, allowed us to push forward the theoretical understanding of several reaction classes.

08 HYDROGEN

Disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y alloys

This dataset contains input and output files from density functional theory (DFT) simulations used to study the disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y (0 ≤ x + y ≤ 1) alloys and their binary end members MnTe, MnSb, and MnBi. The alloys adopt the hexagonal NiAs-type (nickeline) structure and span ternary (MnTexSb1-x, MnTexBi1-x, MnBixSb1-x), and quaternary compositions across the full MnTe–MnSb–MnBi composition triangle. For each alloy composition, the dataset provides DFT calculations in three magnetic configurations: A-type antiferromagnetic (AFM), C-type AFM, and ferromagnetic (FM). Every magnetic configuration folder contains the fully relaxed crystal structure (CONTCAR), VASP input parameters (INCAR), and the main VASP output file (OUTCAR), from which total electronic energies, Mn magnetic moments, lattice parameters, and percent volume changes between magnetic states are extracted. These data are used to construct compositional phase diagrams, evaluate thermodynamic stability (formability), and map magnetoelastic responses across the alloy space. For A-type AFM and FM configurations, additional data are provided as follows: (i) FORCE_CONSTANTS and thermal_properties.yaml files at the top level of A-type_AFM/ and FM/ folders — present only for compositions marked with an asterisk (*) in Table I of the main text. These are derived from Phonopy finite-displacement calculations on full disordered 128-atom supercells and provide vibrational free energy, entropy (Svib)contribution from explicit disorder calculations. (Table I of the associated main manuscript) (ii) A VCA/ subfolder within A-type_AFM/ and FM/, containing FORCE_CONSTANTS and thermal_properties.yaml from Virtual Crystal Approximation phonon calculations (without spin-orbit coupling). VCA data are available for all compositions and are used to estimate vibrational contributions to the Gibbs free energy across the full composition space. (iii) A SOC/ subfolder containing CONTCAR, INCAR, and OUTCAR from spin-orbit coupling calculations, providing relativistic corrections to electronic energies and lattice parameters (Tables S2–S3 of the SM, and Table I of the main manuscript). (iv) A SOC/VCA/ subfolder containing FORCE_CONSTANTS and thermal_properties.yaml from VCA phonon calculations performed within the SOC framework, combining relativistic and vibrational thermodynamic corrections. The computed properties are used to map the AFM–FM magnetic crossover near MnTe0.75Sb0.25, demonstrate disorder- and spin-induced phonon broadening, identify a semiconductor-to-metal crossover, and quantify the pronounced magnetoelastic volume response near the magnetic phase boundary.

36 MATERIALS SCIENCE

Langevin Dynamics modeling of gas-phase ion-ion recombination (Final Technical Report)

A self-consistent trajectory simulation approach to model MN reactions (Fig. 1) which incorporates the probability of electron transfer as a Monte Carlo operator (Fig. 2) was developed and published as Liu et al. J. Chem. Phys. 159, 114111 (2023). The electron transfer probability p ET estimated using the two-state Landau-Zener (LZ) theory was incorporated into classical trajectory simulations to elicit predictions of MN reaction cross-section σ (vacuum) or rate constant β (finite pressure). Electronic structure calculations with multireference configuration interaction (MRCI) and large correlation consistent basis sets were used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the incorporation of electron transfer probability and the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of neutral gas on ions. For H + – H - and Li + – H(D) - pairs, our approach quantitatively agrees with measured speed-dependent cross-sections for up to ~10 5 m/s. For the ion pair Ne + – Cl - , our predictions of the MN rate constant at ~1 torr are a factor of ~2 – 3 higher than the experimentally measured value. Similarly, for Xe + – F - in the pressure range of ~20000 – 80000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. The paradigm of using trajectory simulations to self-consistently model MN reactions is the basis for inclusion of additional non-classical, and static magnetic and electric field effects. Subsequent work, published as Roy et al. focused on modeling recombination rate constant for three ion pairs (rare gas Ar + cation and halide anions): Ar + – Cl - , Ar + – Br - , Ar + – I - , 2) considering spin-orbit couplings in the electronic structure calculations to obtain high-fidelity estimates of the electron transfer probability and incorporated within the classical trajectory simulations to elicit predictions. In addition to calculations of ion-ion recombination rate constants, a classical trajectory simulation technique (published as Roy et al. J. Chem. Phys. 162(9), 094104 (2023)) that uses quaternions to represent orientation of non-spherical particles (ions or aerosol particles) was developed to simulate the recombination of diatomic or more generally, polyatomic molecules. Finally, several other ion pairs such as Ne + – Cl - , Kr + – Cl - , were explored using the developed semi-classical trajectory simulations to understand various challenges in tackling electronic structure calculations. Using empirical approaches to parameterize the electron transfer radius, trajectory simulations were also used to probe the effect of ion number density on MN rate constant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC