Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PP”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Beam focusing and consequences for Doppler backscattering measurements

The phenomenon of focusing of microwave beams in a plasma near a turning-point caustic is discussed by exploiting the analytical solution to the Gaussian beam-tracing equations in the two-dimensional (2-D) linear-layer problem. The location of maximum beam focusing and the beam width at that location are studied in terms of the beam initial conditions. This focusing must be taken into account to interpret Doppler backscattering (DBS) measurements. We find that the filter function that characterises the scattering intensity contribution along the beam path through the plasma is inversely proportional to the beam width, predicting enhanced scattering from the beam focusing region. We show that the DBS signal enhancement for decreasing incident angles between the beam path and the density gradient is due to beam focusing and not due to forward scattering, as was originally proposed by (Gusakov et al., (Plasma Phys. Contr. Fusion, vol. 56, 2014, p. 0250092014, 2017); Plasma Phys. Rep. vol. 43(6), 2017, pp. 605–613). The analytic beam model is used to predict the measurement of the k y density-fluctuation wavenumber power spectrum via DBS, showing that, in an NSTX-inspired example, the spectral exponent of the turbulent, intermediate-to-high k y density-fluctuation spectrum might be quantitatively measurable via DBS, but not the spectral peak corresponding to the driving scale of the turbulent cascade.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Performance and transport in the ARC tokamak

The ARC TM tokamak, a high-field (𝐵 𝑇 = 11.4 T) fusion power plant, under development by Commonwealth Fusion Systems, is studied using a suite of integrated modelling tools to predict its fusion power generation (𝑃𝑓⁡𝑢⁢𝑠), transport and confinement properties. Analysis is based off an ARC operational point scoped first with zero-dimensional (0-D) plasma operational contour (POPCON) modelling to produce 1.13 GW of fusion power. A suite of integrated modelling tools (TRANSP, ASTRA and TORAX) were applied to predict the performance and kinetic profiles of the ARC design point, yielding a range of predicted performance spanning from ∼900 to 1300 MW in rough quantitative agreement with POPCON predictions. The sensitivity of these results to uncertain modelling inputs was probed using scans of pedestal boundary conditions around EPED-predicted values (total pressure and temperature ratios), tungsten concentration and seperatrix density around their nominal assumptions. Pedestal pressure and pedestal top (𝑇 𝑖 /𝑇 𝑒 ) play a large role in 1.5-dimensional performance predictions, able to modify the predicted 𝑃 𝑓⁡𝑢⁢𝑠 by a factor of 2 within reasonable assumptions. High-fidelity core nonlinear gyrokinetic profile predictions, performed using CGYRO (Candy et al. 2016 J. Comput. Phys., vol. 324, pp. 73–93) coupled with the PORTALS (Rodriguez-Fernandez et al. 2024 Nucl. Fusion, vol. 64, 076034; Phys. Plasmas, vol. 31, 2024, 062501) framework, yield substantially lower performance (𝑃 𝑓⁡𝑢⁢𝑠 =677 MW) compared with 0-D and medium-fidelity modelling for nominal assumptions, showing that there is non-negligible uncertainty between models and that future work on SPARC may help resolve discrepancies. Lower overall performance results from significantly reduced volume-averaged densities and temperatures, along with reduced levels of density and temperature peaking. Turbulence and transport are largely dominated by ion temperature gradient across the profile, confirmed by both linear stability and the response of the nonlinear fluxes to changes in gradients, with some impact of kinetic ballooning modes in the deep core. This work represents one of the most complete scoping of potential fusion power plant conditions performed to date. The extensive integrated modelling provides confidence in ARC performance approaching 1 GW, while nonlinear gyrokinetic modelling results in open questions into the physics of density and temperature peaking in fusion-power-plant-relevant operational space. A discussion of results and the role that the SPARC tokamak (Creely et al. 2020 J. Plasma Phys., vol. 86, 865860502) will play in informing ARC design, performance and operation is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Catalytic Upgrading of Pyrolysis Condensables from Postconsumer Polyolefins Using HZSM-5

The conversion of plastic wastes to monomeric olefins is an attractive means for achieving a plastic circular economy. In our study, a fluidized bed reactor converts post-consumer waste high-density polyethylene (HDPE) and polypropylene (PP) to mostly condensed pyrolysis waxes and some oils, preventing carbon loss to gases. The pyrolysis condensables were upgraded to light olefins (C 2 –C 5 ) at carbon yields greater than 76 wt % using the HZSM-5 zeolite catalyst at a post pyrolysis process that employed a micropyrolyzer. These results were comparable to olefin monomer yields from direct ex situ catalytic pyrolysis of the original waste plastics without condensing the vapors, highlighting the potential applicability of this approach in plastic waste recycling. Our results suggest that a centralized catalytic upgrading facility fed by pyrolysis condensables sourced from distributed thermochemical processing plants is a promising pathway to a circular economy. Such an approach enables utilization of available catalytic cracking infrastructure while focusing on setting up distributed thermochemical processing plants close to material recovery facilities. As a result, the energy-dense pyrolysis waxes are more suitable for transportation, contributing to the overall scalability and economic viability of the proposed distributed approach.

10 SYNTHETIC FUELS↗
DOE OSTI2020

Photodissociation and Theory to Investigate Uranium Oxide Cluster Cations

Uranium oxide cluster cations of the form U n O m + are produced by laser vaporization of a depleted uranium rod in a pulsed supersonic expansion. Ions are mass-analyzed and mass-selected with a time-of-flight spectrometer and studied with UV laser multiphoton dissociation. Cations of the stoichiometry UO 2 (UO 3 ) n + were observed as photofragments from all photodissociated cluster cations. (UO 3 ) n + clusters were also observed to result from dissociation of larger (UO 3 ) n + clusters, with UO 3 neutral as a common leaving group. Electronic structure calculations were used to investigate the stability of the prominent uranium oxide cluster cations using density functional theory (DFT) with the hybrid B3LYP exchange-correlation functional and at the CCSD(T) level and cc-pVnZ-PP basis sets (n = D,T), including diffuse orbitals as computational expense and availability permitted. Clustering energies, relative energies and dissociation energies of the cations are reported. The lowest energy neutral (UO 3 ) n clusters up to n = 3 are rings, n = 4 and 5 are chains with very low energy rings, and n = 6 is 3D. Finally, the lowest energy structures for UO 2 (UO 3 ) n + are composed of uranyl-like UO 2 + units bound by bridging oxygens to other UO 2 2+ units for n = 2 and 3, and for n = 4 a more complex 3D structure is predicted.

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

Cryogenic Photoelectron Spectroscopic and Theoretical Study of the Electronic and Geometric Structures of Undercoordinated Osmium Chloride Anions OsCl n – ( n = 3–5)

A series of anionic transition metal halides OsCl n - (n=3–5) have been investigated using a newly developed, home-constructed, cryogenic anion cluster photoelectron spectroscopy. Here, the target anionic species are generated through collision-induced dissociation in a two-stage ion funnel. The measured vertical detachment energies (VDEs) are 3.48 eV, 4.54 eV, and 4.81 eV for n = 3, 4, and 5, respectively. Density functional theory calculations at the B3LYP-D3(BJ)//aug-cc-pVTZ (-pp) level predict the lowest energy structures of OsCl n - (n = 3–5) to be a quintet triangle, quartet square, and quintet square-based pyramid, respectively. The CCSD(T)-calculated VDEs and corresponding adiabatic detachment energies (ADEs) agree well with our experimental measurements. Analysis of the corresponding frontier molecular orbitals (FMOs) and charge density differences suggests that the d-orbitals of the transition metal Os play a primary role in the single-photon detachment processes, and the detached electrons originating from different molecular orbitals are distinguishable.

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

The Remarkable I 2 O 3 Molecule: A New View from Theory

Atmospheric iodine chemistry has garnered increasing attention as a result of increased iodine emissions. A key subset of this chemistry involves iodine oxides (I 2 O 2–5 ), which serve as precursors to particle formation. Among these, I 2 O 3 is the simplest iodine oxide involved in particle formation, but it has remained undetected in the atmosphere. Previous theoretical studies have characterized this peculiar molecule, primarily using energies to refine geometries obtained at low levels of theory. Due to the reemerging interest in I 2 O 3 , this study presents geometries optimized at the CCSD(T)/aug-cc-pwCVTZ-PP level of theory─marking the first instance, to the best of our knowledge, where this system has been studied exclusively with CCSD(T). Harmonic vibrational frequencies were computed at the same level of theory. Final energetics were obtained using the very high level CCSDT(Q) method with basis sets up to quintuple-zeta cardinality (aug-cc-pwCV5Z-PP) and extrapolated to the CBS limit to yield CCSDT(Q)/CBS//CCSD(T)/aug-cc-pwCVTZ-PP energies. These energies include harmonic zero-point vibrational energy corrections and scalar relativistic energy corrections. Additionally, this study discovers new isomers along the I 2 O 3 potential energy surface, a novel contribution to the field. The performance of different computational methods and DFT functionals commonly used in atmospheric chemistry is also assessed relative to high-level theoretical methods.

basis sets↗
DOE OSTI2022

Reactivity of Iridium Complexes of a Triphosphorus-Pincer Ligand Based on a Secondary Phosphine. Catalytic Alkane Dehydrogenation and the Origin of Extremely High Activity

The selective functionalization of alkanes and alkyl groups is a major goal of chemical catalysis. Toward this end, a bulky triphosphine with a central secondary phosphino group, bis(2-di-t-butyl-phosphinophenyl)phosphine ( tBu P H PP), has been synthesized. When complexed to iridium, it adopts a meridional (“pincer”) configuration. The secondary phosphino H atom can undergo migration to iridium to give an anionic phosphido-based–pincer ( tBu PPP) complex. Stoichiometric reactions of the ( tBu PPP)Ir complexes reflect a distribution of steric bulk around the iridium center in which the coordination site trans to the phosphido group is quite crowded; one coordination site cis to the phosphido is even more crowded; and the remaining site is particularly open. The ( tBu PPP)Ir precursors are the most active catalysts reported to date for dehydrogenation of n-alkanes, by about 2 orders of magnitude. The electronic properties of the iridium center are similar to that of well-known analogous ( R PCP)Ir catalysts. Accordingly, DFT calculations predict that ( tBu PPP)Ir and ( tBu PCP)Ir are, intrinsically, comparably active for alkane dehydrogenation. While dehydrogenation by ( R PCP)Ir proceeds through an intermediate trans-(PCP)IrH 2 (alkene), ( tBu PPP)Ir follows a pathway proceeding via cis-(PPP)IrH 2 (alkene), thereby circumventing unfavorable placement of the alkene at the bulky site trans to phosphorus. ( tBu PPP)Ir and ( tBu PCP)Ir, however, have analogous resting states: square planar (pincer)Ir(alkene). Alkene coordination at the crowded trans site is therefore unavoidable in the resting states. Furthermore, the resting state of the ( tBu PPP)Ir catalyst is destabilized by the architecture of the ligand, and this is largely responsible for its unusually high catalytic activity.

02 PETROLEUM↗
DOE OSTI2023

The Stability Transition From Stable to Unstable Frictional Slip With Finite Pore Pressure

Pore fluids are ubiquitous throughout the lithosphere and are commonly invoked as the cause of induced seismicity and slow earthquakes. We perform lab experiments to address these questions for drained fault conditions and low pore pressure. We shear simulated faults at effective normal stress (σ' n ) of 20 MPa and pore pressures Pp from 1 to 4 MPa. We document the full range of lab earthquake behaviors from slow slip to elasto-dynamic rupture and show that slow slip can be explained by the slip rate dependence of the critical rheologic stiffness without dilatancy hardening or other fluid effects. Our fault permeabilities ranges from 10 –18 to 10 –17 m 2 with an initial porosity of 0.1 and estimated fluid diffusion time ≈1 s. Slow slip and quasi-dynamic fault motion may arise from high P p at higher pressures but dilatancy strengthening is not a general requirement.

58 GEOSCIENCES↗
DOE OSTI2022

Muconic acid production from glucose and xylose in Pseudomonas putida via evolution and metabolic engineering

Muconic acid is a bioprivileged molecule that can be converted into direct replacement chemicals for incumbent petrochemicals and performance-advantaged bioproducts. In this study, Pseudomonas putida KT2440 is engineered to convert glucose and xylose, the primary carbohydrates in lignocellulosic hydrolysates, to muconic acid using a model-guided strategy to maximize the theoretical yield. Using adaptive laboratory evolution (ALE) and metabolic engineering in a strain engineered to express the D-xylose isomerase pathway, we demonstrate that mutations in the heterologous D-xylose:H + symporter (XylE), increased expression of a major facilitator superfamily transporter (PP_2569), and overexpression of aroB encoding the native 3-dehydroquinate synthase, enable efficient muconic acid production from glucose and xylose simultaneously. Using the rationally engineered strain, we produce 33.7 g L -1 muconate at 0.18 g L -1 h -1 and a 46% molar yield (92% of the maximum theoretical yield). This engineering strategy is promising for the production of other shikimate pathway-derived compounds from lignocellulosic sugars.

09 BIOMASS FUELS↗
DOE OSTI2025

Mitigating CaCO 3 crystal nucleation and growth through continuous ion displacement via alternating electric fields

Mineral crystal formation poses a challenge on surfaces (e.g., heat exchangers, pipes, membranes, etc.) in contact with super-saturated fluids. Applying alternating currents (AC) to such surfaces can prevent surface crystallization under certain conditions. Here, we demonstrate that ion displacement induced by periodic charging and discharging of the electrical double layer (EDL) inhibits both heterogeneous and homogeneous nucleation (and crystal growth) of CaCO 3 . Titanium sheets (meant to simulate metallic heat exchanger surfaces) are immersed in super-saturated CaCO 3 solutions with a saturation index >11. We show that at relatively high AC frequencies, incomplete EDL formation leads to an alternating electric field that propagates far into the bulk solution, inducing rapid ion migration that overwhelms the Brownian motion of ions. Electrochemical characterization reveals EDL charging/discharging under AC conditions that greatly inhibits precipitation. Operating at 4 V pp , 0.1–10 Hz reduces turbidity by over 96% and reduces CaCO 3 coverage on the metal plates by over 92%. Based on electrokinetic and crystallization models, the ion displacement velocity (exceeding the mean Brownian velocity) and displacement length disrupts ion collision and crystal nucleation. Overall, the technique has potential for preventing mineral crystal formation in heat exchangers and many other industrially relevant systems.

42 ENGINEERING↗
DOE OSTI2025

Addressing genome scale design tradeoffs in Pseudomonas putida for bioconversion of an aromatic carbon source

Genome-scale metabolic models (GSMM) are commonly used to identify gene deletion sets that result in growth coupling and pairing product formation with substrate utilization and can improve strain performance beyond levels typically accessible using traditional strain engineering approaches. However, sustainable feedstocks pose a challenge due to incomplete high-resolution metabolic data for non-canonical carbon sources required to curate GSMM and identify implementable designs. Here we address a four-gene deletion design in the Pseudomonas putida KT2440 strain for the lignin-derived non-sugar carbon source, p-coumarate (p-CA), that proved challenging to implement. We examine the performance of the fully implemented design for p-coumarate to glutamine, a useful biomanufacturing intermediate. In this study glutamine is then converted to indigoidine, an alternative sustainable pigment and a model heterologous product that is commonly used to colorimetrically quantify glutamine concentration. Through proteomics, promoter-variation, and growth characterization of a fully implemented gene deletion design, we provide evidence that aromatic catabolism in the completed design is rate-limited by fumarase hydratase (FUM) enzyme activity in the citrate cycle and requires careful optimization of another fumarate hydratase protein (PP_0897) expression to achieve growth and production. A double sensitivity analysis also confirmed a strict requirement for fumarate hydratase activity in the strain where all genes in the growth coupling design have been implemented. Metabolic cross-feeding experiments were used to examine the impact of complete removal of the fumarase hydratase reaction and revealed an unanticipated nutrient requirement, suggesting additional functions for this enzyme. While a complete implementation of the design was achieved, this study highlights the challenge of completely inactivating metabolic reactions encoded by under-characterized proteins, especially in the context of multi-gene edits.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Observation of deuteron and antideuteron formation from resonance-decay nucleons

High-energy hadronic collisions generate environments characterized by temperatures above 100 MeV (refs.), about 100,000 times hotter than the centre of the Sun. At present, it is therefore unclear how light (anti)nuclei with mass number A of a few units, such as the deuteron, 3 He or 4 He, each bound by only a few MeV, can emerge from these collisions. Here, the ALICE Collaboration reports that deuteron–pion momentum correlations in proton–proton (pp) collisions provide model-independent evidence that about 90% of the observed (anti)deuterons are produced in nuclear reactions following the decay of short-lived resonances, such as the Δ(1232). These findings, obtained by the ALICE Collaboration at the Large Hadron Collider, resolve a gap in our understanding of nucleosynthesis in ultrarelativistic hadronic collisions. Apart from offering insights on how (anti)nuclei are formed in hadronic collisions, the results can be used in the modelling of the production of light and heavy nuclei in cosmic rays and dark-matter decays.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Coulomb-free 1S0 p − p scattering length from the quasi-free p + d → p + p + n reaction and its relation to universality

Abstract The Coulomb-free 1 S 0 proton-proton ( p - p ) scattering length relies heavily on numerous and distinct theoretical techniques to remove the Coulomb contribution. Here, it has been determined from the half-off-the-energy-shell p - p scattering cross section measured at center-of-mass energies below 1 MeV using the quasi-free p + d → p + p + n reaction. A Bayesian data-fitting approach using the expression of the s-wave nucleon-nucleon scattering cross section returned a p - p scattering length $${a}_{pp}=-18.1{7}_{-0.58}^{+0.52}{| }_{stat}\pm 0.0{1}_{syst}$$ a p p = − 18.1 7 − 0.58 + 0.52 ∣ s t a t ± 0.0 1 s y s t fm and effective range r 0 = 2.80 ± 0.05 s t a t ± 0.001 s y s t fm. A model based on universality concepts has been developed to interpret this result. It accounts for the short-range interaction as a whole, nuclear and residual electromagnetic, according to what the s-wave phase-shift δ does in the description of low-energy nucleon-nucleon scattering data. We conclude that our parameters are representative of the short-range physics and propose to assess the charge symmetry breaking of the short-range interaction instead of just the nuclear interaction. This is consistent with the current understanding that the charge dependence of nuclear forces is due to different masses of up-down quarks and their electromagnetic interactions. This achievement suggests that these properties have a lesser than expected impact in the context of the charge symmetry breaking.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

The highly exothermic hydrogen abstraction reaction H 2 Te + OH → H 2 O + TeH: comparison with analogous reactions for H 2 Se and H 2 S

The “gold standard” CCSD(T) method is adopted along with the correlation consistent basis sets up to aug-cc-pV5Z-PP to study the mechanism of the hydrogen abstraction reaction H 2 Te + OH. Here, the predicted geometries and vibrational frequencies for reactants and products are in good agreement with the available experimental results. With the ZPVE corrections, the transition state in the favorable pathway of this reaction energetically lies 1.2 kcal mol -1 below the reactants, which is lower than the analogous relative energies for the H 2 Se + OH reaction (-0.7 kcal mol -1 ), the H 2 S + OH reaction (+0.8 kcal mol -1 ) and the H 2 O + OH reaction (+9.0 kcal mol -1 ). Accordingly, the exothermic reaction energies for these related reactions are predicted to be 47.8 (H 2 Te), 37.7 (H 2 Se), 27.1 (H 2 S), and 0.0 (H 2 O) kcal mol -1 , respectively. Geometrically, the low-lying reactant complexes for H 2 Te + OH and H 2 Se + OH are two-center three-electron hemibonded structures, whereas those for H 2 S + OH and H 2 O + OH are hydrogen-bonded. With ZPVE and spin–orbit coupling corrections, the relative energies for the reactant complex, transition state, product complex, and the products for the H 2 Te + OH reaction are estimated to be -13.1, -1.0, -52.0, and -52.6 kcal mol -1 , respectively. Finally, twenty-eight DFT functionals have been tested systematically to assess their ability in describing the potential energy surface of the H 2 Te + OH reaction. The best of these functionals for the corresponding energetics are -9.9, -1.4, -46.4, and -45.4 kcal mol -1 (MPWB1K), or -13.1, -2.4, -57.1, and -54.6 kcal mol -1 (M06-2X), respectively.

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

Catalytic production of tetrahydropyran (THP): a biomass-derived, economically competitive solvent with demonstrated use in plastic dissolution

Tetrahydropyran (THP) is a five-carbon heterocyclic ether that is non-carcinogenic, non-peroxide forming, biodegradable, and economically competitive with tetrahydrofuran (THF) as a solvent. In this work, THP has been synthesized from renewable biomass at >99.8% selectivity and 98% yield via hydrogenation of furfural-derived 3,4-dihydropyran (DHP) over Ni/SiO 2 in a continuous flow reactor at 150–200 °C. The apparent activation energy of THP formation is 31 kJ mol -1 , and the reaction orders with respect to the partial pressures of H 2 , DHP, and THP are: 2, 1, and -0.3. The kinetic data has been fitted to a Hougen–Watson model where the rate limiting step is the hydrogenation of adsorbed DHP. Ni/SiO 2 is shown to have a low deactivation rate constant of 0.012 h -1 over 100 h time on stream and can be regenerated in situ. As a performance advantage, THP is shown to be resistant to ring opening polymerization under strongly acidic conditions that THF is not, revealing it to be a superior solvent. Further, the minimum selling price of THP is competitive with the market price of THF ($$900 – 1400 per ton ) at a DHP feedstock cost of $1000 per ton. Conductor-like Screening Model for Real Solvents (COSMO-RS) and molecular dynamics (MD) simulations with 1008 solvents and 8 common plastics have demonstrated that THP can serve as an alternative solvent to THF, 2-methyltetrahydrofuran (MeTHF), and cyclopentyl methyl ether (CPME) for plastic dissolution, especially low-density polyethylene (LDPE), polypropylene (PP), polystyrene (PS) and polyvinyl chloride (PVC). This work establishes THP as a green solvent with excellent thermal, chemical and peroxidative stability that can be used for numerous applications, including waste plastic recycling.

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

Radio frequency mixing modules for superconducting qubit room temperature control systems

As the number of qubits in nascent quantum processing units increases, the connectorized RF (radio frequency) analog circuits used in first generation experiments become exceedingly complex. The physical size, cost, and electrical failure rate all become limiting factors in the extensibility of control systems. We have developed a series of compact RF mixing boards to address this challenge by integrating I/Q quadrature mixing, intermediate frequency/LO (local oscillator)/RF power level adjustments, and direct current bias fine tuning on a 40 × 80 mm 2 four-layer printed circuit board with electromagnetic interference shielding. The RF mixing module is designed to work with RF and LO frequencies between 2.5 and 8.5 GHz. The typical image rejection and adjacent channel isolation are measured to be ~27 dBc and ~50 dB. By scanning the drive phase in a loopback test, the module short-term amplitude and phase linearity are typically measured to be 5 ×10 -4 (V pp /V mean ) and 1 ×10 -3 radian (pk-pk). The operation of the RF mixing board was validated by integrating it into the room temperature control system of a superconducting quantum processor and executing randomized benchmarking characterization of single and two qubit gates. We measured a single-qubit process infidelity of 9.3(3) × 10 -4 and a two-qubit process infidelity of 2.7(1) × 10 -2 .

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Supervised learning and the finite-temperature string method for computing committor functions and reaction rates

A central object in the computational studies of rare events is the committor function. Though costly to compute, the committor function encodes complete mechanistic information of the processes involving rare events, including reaction rates and transition-state ensembles. Under the framework of transition path theory, Rotskoff et al. [Proceedings of the 2nd Mathematical and Scientific Machine Learning Conference, Proceedings of Machine Learning Research (PLMR, 2022), Vol. 145, pp. 757–780] proposes an algorithm where a feedback loop couples a neural network that models the committor function with importance sampling, mainly umbrella sampling, which collects data needed for adaptive training. Here, in this work, we show additional modifications are needed to improve the accuracy of the algorithm. The first modification adds elements of supervised learning, which allows the neural network to improve its prediction by fitting to sample-mean estimates of committor values obtained from short molecular dynamics trajectories. The second modification replaces the committor-based umbrella sampling with the finite-temperature string (FTS) method, which enables homogeneous sampling in regions where transition pathways are located. We test our modifications on low-dimensional systems with non-convex potential energy where reference solutions can be found via analytical or finite element methods, and show how combining supervised learning and the FTS method yields accurate computation of committor functions and reaction rates. We also provide an error analysis for algorithms that use the FTS method, using which reaction rates can be accurately estimated during training with a small number of samples. The methods are then applied to a molecular system in which no reference solution is known, where accurate computations of committor functions and reaction rates can still be obtained.

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

Modeling transient edge plasma transport with dynamic recycling

The work presents numerical simulation studies of the role that dynamic plasma recycling on the main wall and divertor target surfaces plays in transient edge plasma transport phenomena, such as edge localized modes (ELMs). The studies are performed by coupling the edge plasma transport code UEDGE [Rognlien et al., J. Nucl. Mater. 196–198, 347 (1992)] and the wall reaction–diffusion transport code FACE [Smirnov et al., Fusion Sci. Technol. 71, 75 (2017)]. The two-dimensional, time-dependent, two-way coupling of the codes, in a realistic tokamak geometry, is accomplished using the Integrated Plasma Simulator framework [Elwasif et al., in 18th Euromicro Conference on Parallel, Distributed and Network-Based Processing (PDP 2010), Pisa, Italy (IEEE, 2010), pp. 419–427] for all modeled material plasma boundaries. The simulations show that dynamic plasma recycling has substantially different characteristics on the main wall and on the divertor plates. It is demonstrated that during an ELM cycle the outer wall can dynamically absorb and release a number of particles comparable to that expelled by the ELM from the core plasma, by far exceeding the dynamic retention capacity of the divertor surfaces. The resulting evolution of the edge and divertor plasma conditions during an ELM cycle is analyzed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗