Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kr”

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 361 records · Page 20

Quantifying uncertainty in physics-based predictions of rare-isotope production cross sections via Bayesian-inspired model averaging across nuclear mass tables

Accurate prediction of fragmentation cross sections is essential for rare-isotope beam production, planning new-isotope searches, and designing experiments to study the most exotic regions of the nuclear chart. However, existing reaction models and phenomenological cross-section parametrizations often exhibit significant deviations over broad regions of mass and charge. In this work, a Bayesian-inspired model-averaging framework is developed to combine abrasion-ablation (AA) calculations based on multiple nuclear mass tables into a single statistically weighted estimate. For the calibrated systems, the model weights are assigned empirically according to the relative quality of fit to measured cross sections, thereby reducing systematic model bias while preserving the underlying physics content of the AA description. The weights are constrained using proton-rich fragmentation data for the 78 Kr and 124 Xe projectiles. The resulting parameter trends are then propagated to the 92 Mo and 144 Sm systems through a controlled scaling procedure. In the present implementation, the excitation-energy prescription is fixed, while the averaging is performed across nuclear-mass inputs; the framework provides both weighted cross sections and associated uncertainty estimates. Applied to proton-rich fragmentation, the present approach provides a practical basis for interpolation and limited extrapolation in regions relevant to rare-isotope production. The resulting predictions are used to assess the production of very proton-rich nuclei, and candidate new isotopes are discussed.

Bayesian methods↗

Compliant Intramedullary Stems for Joint Reconstruction

The longevity of current joint replacements is limited by aseptic loosening, which is the primary cause of non-infectious failure for hip, knee, and ankle arthroplasty. Aseptic loosening is typically caused either by osteolysis from particulate wear, or by high shear stresses at the bone-implant interface from over-constraint. Our objective was to demonstrate feasibility of a compliant intramedullary stem that eliminates over-constraint without generating particulate wear. The compliant stem is built around a compliant mechanism that permits rotation about a single axis. We first established several models to understand the relationship between mechanism geometry and implant performance under a given angular displacement and compressive load. We then used a neural network to identify a design space of geometries that would support an expected 100-year fatigue life inside the body. We additively manufactured one representative mechanism for each of three anatomic locations, and evaluated these prototypes on a KR-210 robot. The neural network predicts maximum stress and torsional stiffness with 2.69% and 4.08% error respectively, relative to finite element analysis data. We identified feasible design spaces for all three of the anatomic locations. Simulated peak stresses for the three stem prototypes were below the fatigue limit. Benchtop performance of all three prototypes was within design specifications. Our results demonstrate the feasibility of designing patient- and joint-specific compliant stems that address the root causes of aseptic loosening. Guided by these results, we expect the use of compliant intramedullary stems in joint reconstruction technology to increase implant lifetime.

60 APPLIED LIFE SCIENCES↗

Ta–Zr carbides: Synthesis advances via carbothermal reduction and defect evolution observed through transmission electron microscopy ion irradiation

The thermodynamic stability of six distinct compositions within the Zr—Ta—C ternary system is investigated in this study, marking the first report of their synthesis through carbothermic reduction in vacuum. A prolonged annealing process at 2200°C enabled high densification and phase equilibrium. Detailed phase identification and microstructural characterization through microscopy and X-ray diffraction techniques revealed clear compositional trends and stable phase formations. Two compositions ((Ta 0.2 Zr 0.8 )C 0.6 and (Ta 0.5 Zr 0.5 )C 1 ) were selected for ion irradiation experiments using 200 keV Kr + at 600°C—representing the first-ever irradiation study on the Zr—Ta—C system. The findings indicated defect accumulation and nanoscale cavity formation without any evidence of amorphization, highlighting the system's structural stability under irradiation. Together, the synthesis and irradiation results provide a basis for further investigation of the system and suggest its relevance for applications under extreme environments.

36 MATERIALS SCIENCE↗

The s process in massive stars, a benchmark for neutron capture reaction rates

A clear definition of the contribution from the slow neutron-capture process (s process) to the solar abundances between Fe and the Sr-Zr region is a crucial challenge for nuclear astrophysics. Robust s-process predictions are necessary to disentangle the contribution from other stellar processes producing elements in the same mass region. Nuclear uncertainties are affecting s-process calculations, but most of the needed nuclear input are accessible to present nuclear experiments or they will be in the near future. Neutron-capture rates have a great impact on the s process in massive stars, which is a fundamental source for the solar abundances of the lighter s-process elements heavier than Fe (weak s-process component). In this work we present a new nuclear sensitivity study to explore the impact on the s process in massive stars of 86 neutron-capture rates, including all the reactions between C and Si and between Fe and Zr. We derive the impact of the rates at the end of the He-burning core and at the end of the C-burning shell, where the 22 Ne(α,n) 25 Mg reaction is is the main neutron source. We confirm the relevance of the light isotopes capturing neutrons in competition with the Fe seeds as a crucial feature of the s process in massive stars. For heavy isotopes we study the propagation of the neutron-capture uncertainties, finding a clear difference of the impact of Fe and Co isotope rates with respect to the rates of heavier stable isotopes. The local uncertainty propagation due to the neutron-capture rates at the s-process branching points is also considered, discussing the example of 85 Kr. The complete results of our study for all the 86 neutron-capture rates are available online. Finally, we present the impact on the weak s process of the neutron-capture rates included in the new ASTRAL library (v0.2).

79 ASTRONOMY AND ASTROPHYSICS↗

Demonstration of event position reconstruction based on diffusion in the NEXT-white detector

Noble element time projection chambers are a leading technology for rare event detection in physics, such as for dark matter and neutrinoless double beta decay searches. Time projection chambers typically assign event position in the drift direction using the relative timing of prompt scintillation and delayed charge collection signals, allowing for reconstruction of an absolute position in the drift direction. In this paper, alternate methods for assigning event drift dis tance via quantification of electron diffusion in a pure high pressure xenon gas time projection chamber are explored. Data from the NEXT-White detector demonstrate the ability to achieve good position assignment accuracy for both high and low-energy events. Using point-like energy deposits from 83m Kr calibration electron captures (E ~ 45 keV), the position of origin of low-energy events is determined to 2 cm precision with bias < 1 mm. A convolutional neural network approach is then used to quantify diffusion for longer tracks (E ≥ 1.5 MeV), from radiogenic electrons, yielding a precision of 3 cm on the event barycenter. The precision achieved with these methods indicates the feasibility energy calibrations of better than 1% FWHM at Q ββ in pure xenon, as well as the potential for event fiducialization in large future detectors using an alternate method that does not rely on primary scintillation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Measurement of the electric potential and the magnetic field in the shifted analysing plane of the KATRIN experiment

The projected sensitivity of the effective electron neutrino-mass measurement with the KATRIN experiment is below 0.3 eV (90 % CL) after 5 years of data acquisition. The sensitivity is affected by the increased rate of the background electrons from KATRIN’s main spectrometer. A special shifted-analysing-plane (SAP) configuration was developed to reduce this background by a factor of two. The complex layout of electromagnetic fields in the SAP configuration requires a robust method of estimating these fields. We present in this paper a dedicated calibration measurement of the fields using conversion electrons of gaseous 83m Kr, which enables the neutrino-mass measurements in the SAP configuration.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Cross-polarized common-path temporal interferometry for high-sensitivity strong-field ionization measurements

Absolute density measurements of low-ionization-degree or low-density plasmas ionized by lasers are very important for understanding strong-field physics, atmospheric propagation of intense laser pulses, Lidar etc. A cross-polarized common-path temporal interferometer using balanced detection was developed for measuring plasma density with a sensitivity of ∼0.6 mrad, equivalent to a plasma density-length product of ∼2.6 × 10 13 cm −2 if using an 800 nm probe laser. By using this interferometer, we have investigated strong-field ionization yield versus intensity for various noble gases (Ar, Kr, and Xe) using 800 nm, 55 fs laser pulses with both linear (LP) and circular (CP) polarization. The experimental results were compared to the theoretical models of Ammosov-Delone-Krainov (ADK) and Perelomov-Popov-Terent’ev (PPT). We find that the measured phase change induced by plasma formation can be explained by the ADK theory in the adiabatic tunneling ionization regime, while PPT model can be applied to all different regimes. We have also measured the photoionization and fractional photodissociation of molecular (MO) hydrogen. By comparing our experimental results with PPT and MO-PPT models, we have determined the likely ionization pathways when using three different pump laser wavelengths of 800 nm, 400 nm, and 267 nm.

47 OTHER INSTRUMENTATION↗

Constraining Bedrock Groundwater Residence Times in a Mountain System with Environmental Tracer Observations and Bayesian Uncertainty Quantification: Modeling and Data Package

Groundwater residence times provide fundamental descriptions of hydrologic dynamics and mixing processes in mountainous watersheds. Yet, few observational datasets that can constrain groundwater residence times over broad timescales are available in high elevation mountain systems. Here we present field observations from May 2021 of dissolved noble gases (He, Ne, Ar, Kr, and Xe), Chloroflourcarbons (CFCs), Sulfurhexaflouride (SF6), and tritium (3H) sampled from the Pumphouse Lower Montane study site (wells PLM1, PLM6, and PLM7) within the East River Watershed, Colorado. The presented noble gas (PLM_noblegas_2021.csv) and environmental tracer (PLM_tracers_2021.csv) observation datasets, along with the associated modeling scripts, aide in quantifying groundwater residence times and recharge conditions in a high elevation mountain system. Furthermore, the modeling scripts quantify groundwater residence time and noble gas recharge condition uncertainties using a novel Markov-chain Monte Carlo approach. All data modeling scripts are written in the Python code.

54 ENVIRONMENTAL SCIENCES↗

Fluid – Fluid Interfacial Area and Its Impact on Relative Permeability - A Pore Network Modeling Study

Abstract Relative permeability (kr) is commonly modeled as an empirical function of phase saturation. Although current empirical models can provide a good match of one or two measured relative permeabilities using saturation alone, they are unable to predict relative permeabilities well when there is hysteresis or when physical properties such as wettability change. Further, current models often result in relative permeability discontinuities that can cause convergence and accuracy problems in simulation. To overcome these problems, recent research has modeled relative permeability as a state function of both saturation (S) and phase connectivity (X). Pore network modeling (PNM) data, however, shows small differences in relative permeability for the same S-X value when approached from a different flow direction. This paper examines the impact of one additional Minkowski parameter (Mecke and Arns, 2005), the fluid-fluid interfacial area, on relative permeability to identify if that satisfactorily explains this discrepancy. We calculate the total fluid-fluid interfacial areas (IA) during two-phase (oil/water) flow in porous media using pore network modeling. The area is calculated from PNM simulations using the areas associated with corners and throats in pore elements of different shapes. The pore network is modeled after a Bentheimer sandstone, using square, triangular, and circular pore shapes. Simulations were conducted for numerous primary drainage and imbibition cycles at a constant contact angle of 0° for the wetting phase. Simultaneous measurements of capillary pressure, relative permeability, saturation, and phase connectivity are made for each displacement. Fluid-fluid interfacial area is calculated from the PNM capillary pressure, the fluid location in the pore elements, and the pore element dimensional data. The results show that differences in the relative permeability at the same (S,X) point is explained well by differences in the fluid-fluid interfacial area (IA). That is, for a larger change in IA at these intersection points, the permeability difference is greater. That difference in relative permeability approaches zero as the difference in IA approaches zero. This confirms that relative permeability can be modeled better as a unique function of S, X, and IA. The results also show that an increase in IA restricts flow decreasing the nonwetting (oil) phase permeability. This decrease is caused by an increase in the throat area fraction compared to the corner area as the total area IA increases. The wetting phase relative permeability, however, shows the inverse trend, in that its relative permeability is greater when IA becomes larger owing to a greater fraction of the total area associated with the corners. The area IA, however, impacts the nonwetting phase relative permeability more than the wetting phase relative permeability. Corner flow improves the wetting phase relative permeability because the wetting phase is continuous there. Finally, a sensitivity analysis shows that relative permeability a is more sensitive to change in S than they are for IA for the case studied implying that if only two parameters are used to model relative permeability it is better to choose S and X.

Engineering↗

Fission Gas Monitoring System Efficiency Calibration Summary Report

The Fission Gas Monitoring System was designed to work in conjunction with the Fuel Accident Condition Simulator Furnace and is used to measure fission gas releases during post-irradiation fuel heating tests. A high-purity helium sweep gas is swept past the heated fuel sample in the furnace and routed to a fission gas monitoring system, which cryrogenically traps the krypton (Kr) and xenon (Xe) fission gases. As the radioactive species of the trapped gas decays, high-purity germanium detectors monitor the amount of fission gas in each trap. The primary fission product present in this test is 85Kr, which has a 10.76-year half-life. The presence of 85Kr indicates particle defects and iodine release. There is also the potential for re-irradiating the fuel compacts in the Neutron Radiography Reactor in the Hot Fuel Examination Facility basement. With re-irradiation, it is possible to have xenon fission gas present, specifically 133Xe, which has a 5.243-day half-life. Re-irradiation also produces 131I, which is important for design-basis accidents and regulatory licensing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Effect of Impurity Atoms on Thermonuclear Yield

The presence of impurity atoms in an inertial confinement fusion (ICF) implosion can greatly affect the burning of the thermonuclear fuel. This is a serious concern for attaining ignition as high atomic number atoms, such as carbon, can mix into the hot spot and reach the deuterium–tritium (DT) fuel, absorb energy from the implosion, and radiate away that energy. The temperature of the DT fuel drops and reduces the reactivity of the fuel resulting in lower thermonuclear yield. Fuel impurities can also be a useful diagnostic as when higher-Z atoms, such as Ar, Kr, and Xe, are purposefully mixed into the fuel. From the spectra emitted by these atoms, the temperature and density of the fuel can be deduced. The long history of using such tracer atoms by the national ICF program shows the great advances in experimental technique, diagnostic instruments, and atomic-physics modeling over the past 50 years. A particular set of experiments, the “High-Z” campaign, is reviewed in depth as this was the most systematic study of tracer element composition and quantity ever reported. Difficulties in interpreting the results of those experiments due to the limited code capabilities at that time are discussed. A main concern for the High-Z campaign was the then current state of non-Local Thermodynamic Equilibrium (nLTE) atomic physics models. The present state of nLTE modeling and its implementation in Los Alamos’s multi-physics codes is discussed in this report. To further the understanding of the complex interaction of high-Z atoms in the gas, a survey of relevant existing experimental platforms at the National Ignition Facility (NIF) and the Omega Laser Facility are presented along with the current state of pertinent diagnostics at NIF and Omega. This review concludes with a discussion of the trade-off available between temperature and volume of the imploded ICF capsule with the goal of maximizing the total yield.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Fission Gas Measurements in Spent Fuel: Literature Review

The fission-derived noble gases, xenon(Xe) and krypton(Kr), have long been proposed as a tool for supporting nuclear safeguards and declarations verification. A vast body of literature outlines how these non-reactive fission gases released during reactor operation and shutdowns correlate to burnup level, power history, other reactor specifications and reprocessing activities. However, irradiation conditions can vary appreciably both radially and axially within a fuel pin.

07 ISOTOPE AND RADIATION SOURCES↗

Spatiotemporal Analyses of Groundwater and Shoreline Cr(VI) Concentrations in the 100 Areas at Hanford

Cleanup efforts have been ongoing since the late 1990s to remediate contaminated waste sites and groundwater in the 100 Areas at the U.S. Department of Energy (DOE) Hanford Site. One of the primary contaminants of concern is hexavalent chromium (Cr(VI)), which was used as a corrosion inhibitor in cooling water for nuclear reactors that formerly operated along the shoreline of the Columbia River. Cleanup efforts have included 1) removal, treatment (as needed), and disposal of contaminated sediments; 2) in situ redox manipulation as a permeable reactive barrier; 3) pump-and-treat; 4) soil flushing; and 5) monitored natural attenuation. DOE’s annual groundwater monitoring reports document the significant reductions in Cr(VI) plume areas that have occurred over the past 10 years or more as a result of these cleanup efforts. The Record of Decision for the 100-HR-3 operable unit specified a cleanup level (CUL) for Cr(VI) in groundwater of 48 µg/L to protect human receptors, and a surface water CUL of 10 µg/L to protect aquatic organisms in the Columbia River. The Record of Decision did not specify point-of-compliance locations for the surface water CUL. Data for 2019 from the six groundwater operable units (OUs) in the 100 Areas indicate that the 48 μg/L groundwater CUL has been achieved in 100% of the wells in the 100-BC and 100-NR OUs, and in 89- 97% of the wells in the other OUs (100-KR, 100-HR-D, 100-HR-H, 100-FR). Data for 2019 indicate that 100% of the aquifer tubes monitored for Cr(VI) in the 100 Areas have concentrations below the 48 μg/L groundwater CUL. However, the 10 μg/L standard has not yet been consistently achieved for both inland groundwater monitoring wells and shoreline aquifer tubes. This report describes a series of data analyses performed to identify consistent relationships, if any, between inland well and shoreline Cr(VI) concentrations within the 100 Areas. To this end, select monitoring data for Cr(VI) measured in groundwater and aquifer tubes at the 100 Areas were analyzed for a 10-year period—2010 to 2019. Relationships between inland groundwater plumes and surface-water points of discharge in and along the Columbia River were examined through several analyses that included inland well and aquifer tube concentrations as a function of distance from the shoreline, evaluation of cumulative probability plots, trend analysis, correlation analysis, cluster analysis, and identification of plume trajectories for each of the 100 Areas. The analyses did not identify consistent relationships between inland groundwater Cr(VI) concentrations and shoreline concentrations within the 100 Areas due to several confounding factors influencing groundwater flow directions and Cr(VI) concentrations. The proximity of groundwater Cr(VI) plumes to the river, and the highly dynamic nature of the river, influence the transport behavior of the plumes and create challenges for quantifying attenuation of Cr(VI) between the inland monitoring wells and shoreline concentrations. Other factors contributing to temporal and spatial Cr(VI) concentrations, as supported by some of the data analyses, include the presence of vadose zone sources, variable sorption behavior, and complexities associated with Cr(VI) mass transfer between the upper and lower aquifers and their interactions with the river. Hence, monitoring to assess compliance with target CULs will need to be determined for each area individually since several factors influence Cr(VI) concentrations in the 100 Areas.

54 ENVIRONMENTAL SCIENCES↗

Enabling System Technologies to Improve the Economies and Performance of Existing LWRs and Advanced BWR Plants: Improving Offgas System Performance. Final CRADA Report

The project goal is to study adsorbent materials and adsorption conditions in order to develop recommendations to improve the current BWR off-gas treatment operations with respect to cost savings. ANL tasks included a literature survey on relevant information for the adsorption of radionuclides (e.g., Xe and Kr), including adsorption coefficients on various charcoals, the effect of radionuclide concentration and moisture, and alternative sorption materials; participate in meetings and teleconferences as needed; prepare reports and support efforts INL and other organization Team Members.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Air Pathway Dose Modeling for the E-Area Low-Level Waste Facility

The US Department of Energy (DOE) Order 435.1 performance assessment (PA) process (USDOE 1999) prescribes a performance objective (10 mrem y -1 ) for evaluating atmospheric releases of radionuclides from DOE low-level waste (LLW) disposal facilities. The potential dose to an individual from exposure to radionuclides released into the atmosphere from LLW disposals can be estimated by application of radionuclide-specific dose-release factors (DRFs) (mrem Ci -1 ) to estimated flux rates (Ci y -1 ) at a particular time and location. In Revision 0 of this report, Lee (2006) calculated DRFs for potential atmospheric releases of 15 volatile radionuclides from seven ELLWF disposal units in support of the 2008 ELLWF PA (WSRC 2008). The mainframe version of the US Environmental Protection Agency’s (EPA) dose model CAP88 was used in the Revision 0 assessment. In Revision 1 of this report (Dixon and Minter 2017), the site-specific Savannah River National Laboratory (SRNL) atmospheric dose models MAXDOSE-SR version 2013 (Stone and Jannik 2013a) and MAXINE version 2017 (Bell 2017) were used to calculate DRFs for ten radionuclides. Revision 2 provided updated DRFs using the current version of the US Environmental Protection Agency’s (EPA) dose model CAP88-PC Version 4.1. After Revision 2 was finalized, DRFs were requested for Ar-37, Ar-39, Kr-83m, and Hg-206. Revision 3 provides the additional DRFs calculated using the same methods as used in Revision 2 of this document. The method for estimating new DRFs for the ELLWF disposal units is described in this report.

61 RADIATION PROTECTION AND DOSIMETRY↗

IVEM Contributions to CINR Penn State MT17-12797: In situ Ion Irradiation to Add Irradiation Assisted Grain Growth to the MARMOT Tool

In-situ TEM experiments were performed at the IVEM-Tandem Facility in Argonne National Laboratory to study the irradiation-induced grain growth of nanocrystalline UO 2 and CeO 2 thin films using 1 MeV Kr ions at various temperatures from 50K to 1073K. The grain growth kinetics was studied under both irradiation and thermal annealing conditions. Bright field and dark field TEM images and diffraction patterns were systematically recorded as a function of temperature, irradiation dose and annealing time. A comprehensive grain growth data of UO 2 was collected for further analysis and for the implementation to the MARMOT tool.

36 MATERIALS SCIENCE↗

Chemical Thermodynamic Modeling of Molten Salts to Support Off-Gas Abatement Systems

The reprocessing of used nuclear fuel by any means will liberate gaseous fission products, such as hydrogen ( 3 H), carbon ( 14 C), noble gases ( 85 Kr), and halogens ( 129 I), from the irradiated fuel. These elements will distribute through chemical processing operations and partition into process off-gas streams, and the specific volatile release fractions will be dictated by the chemical and physical properties of the system. A recent assessment found that there were significant knowledge gaps regarding the release quantities of volatile radionuclides from individual unit operations. These knowledge gaps limited the ability to determine what dedicated off-gas treatment technologies could be required for electrochemical-based reprocessing facilities. Unfortunately, experimental efforts to quantify release fractions are limited by the challenges associated with performing experiments using irradiated fuel. This report documents preliminary thermodynamic predictions of iodine and tritium release from chloride-based molten salts as part of an effort to better direct resources toward those experiments (both simulant and irradiated) that will be of the greatest impact. It was predicted that less than 0.5% of tritium was expected to be released and that nearly all of that amount would be released as H 2 . Thermochemical data for hydrogen (H 2 ) are of high fidelity, and no additional validation is recommended. Less than 0.05% of iodine was predicted to be released, with the primary volatile species being Cs 2 I 2 . Unlike the tritium predictions, the data underpinning the iodine release predictions are of low quality. It is recommended that a limited experimental program be dedicated to expanding the current physical property and thermodynamic property data for iodine in electrochemical processing and molten salt conditions, including vapor pressure measurements for the dimerized species predicted to comprise the majority of iodine release. Validating and improving the thermodynamic properties used in predictive modeling can reduce the need for expensive testing with irradiated fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

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↗