Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “break”

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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

P not PQ

Parity solutions to the strong CP problem are a compelling alternative to approaches based on Peccei-Quinn symmetry, particularly given the expected violation of global symmetries in a theory of quantum gravity. The most natural of these solutions break parity at a low scale, giving rise to a host of experimentally accessible signals. We assess the status of the simplest parity-based solution in light of LHC data and flavor constraints, highlighting the prospects for near-future tests at colliders, tabletop experiments, and gravitational wave observatories. The origin of parity breaking and associated gravitational effects play crucial roles, providing new avenues for discovery through EDMs and gravity waves. These experimental opportunities underline the promise of generalized parity, rather than Peccei-Quinn symmetry, as a robust and testable solution to the strong CP problem.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Non-invertible global symmetries and completeness of the spectrum

It is widely believed that consistent theories of quantum gravity satisfy two basic kinematic constraints: they are free from any global symmetry, and they contain a complete spectrum of gauge charges. For compact, abelian gauge groups, completeness follows from the absence of a 1-form global symmetry. However, this correspondence breaks down for more general gauge groups, where the breaking of the 1-form symmetry is insufficient to guarantee a complete spectrum. We show that the correspondence may be restored by broadening our notion of symmetry to include non-invertible topological operators, and prove that their absence is sufficient to guarantee a complete spectrum for any compact, possibly disconnected gauge group. In addition, we prove an analogous statement regarding the completeness of twist vortices: codimension-2 objects defined by a discrete holonomy around their worldvolume, such as cosmic strings in four dimensions. We discuss how this correspondence is modified in various, more general contexts, including non-compact gauge groups, Higgsing of gauge theories, and the addition of Chern-Simons terms. Finally, we discuss the implications of our results for the Swampland program, as well as the phenomenological implications of the existence of twist strings.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Studying strangeness and baryon production mechanisms through angular correlations between charged Ξ baryons and identified hadrons in pp collisions at $\sqrt{s}$ = 13 TeV

The angular correlations between charged Ξ baryons and associated identified hadrons (pions, kaons, protons, Λ baryons, and Ξ baryons) are measured in pp collisions at $\sqrt{s}$ = 13 TeV with the ALICE detector to give insight into the particle production mechanisms and balancing of quantum numbers on the microscopic level. In particular, the distribution of strangeness is investigated in the correlations between the doubly-strange Ξ baryon and mesons and baryons that contain a single strange quark, K and Λ. As a reference, the results are compared to Ξπ and Ξp correlations, where the associated mesons and baryons do not contain a strange valence quark. These measurements are expected to be sensitive to whether strangeness is produced through string breaking or in a thermal production scenario. Furthermore, the multiplicity dependence of the correlation functions is measured to look for the turn-on of additional particle production mechanisms with event activity. The results are compared to predictions from the string-breaking model P YTHIA 8, including tunes with baryon junctions and rope hadronisation enabled, the cluster hadronisation model H ERWIG 7, and the core-corona model E POS-LHC . While some aspects of the experimental data are described quantitatively or qualitatively by the Monte Carlo models, no model can match all features of the data. These results provide stringent constraints on the strangeness and baryon number production mechanisms in pp collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

A more natural composite Higgs model

Composite Higgs models provide an attractive solution to the hierarchy problem. However, many realistic models suffer from tuning problems in the Higgs potential. There are often large contributions from the UV dynamics of the composite resonances to the Higgs potential, and tuning between the quadratic term and the quartic term is required to separate the electroweak breaking scale and the compositeness scale. We consider a composite Higgs model based on the SU(6)/Sp(6) coset, where an enhanced symmetry on the fermion resonances can minimize the Higgs quadratic term. Moreover, a Higgs quartic term from the collective symmetry breaking of the little Higgs mechanism can be realized by the partial compositeness couplings between elementary Standard Model fermions and the composite operators, without introducing new elementary fields beyond the Standard Model and the composite sector. The model contains two Higgs doublets, as well as several additional pseudo-Nambu-Goldstone bosons. To avoid tuning, the extra Higgs bosons are expected to be relatively light and may be probed in the future LHC runs. The deviations of the Higgs couplings and the weak gauge boson couplings also provide important tests as they are expected to be close to the current limits in this model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Chern-Weil global symmetries and how quantum gravity avoids them

We draw attention to a class of generalized global symmetries, which we call “Chern-Weil global symmetries,” that arise ubiquitously in gauge theories. The Noether currents of these Chern-Weil global symmetries are given by wedge products of gauge field strengths, such as F 2 $\wedge$ H 3 and tr(F$^2_2$), and their conservation follows from Bianchi identities. As a result, they are not easy to break. However, it is widely believed that exact global symmetries are not allowed in a consistent theory of quantum gravity. As a result, any Chern-Weil global symmetry in a low-energy effective field theory must be either broken or gauged when the theory is coupled to gravity. In this paper, we explore the processes by which Chern-Weil symmetries may be broken or gauged in effective field theory and string theory. We will see that many familiar phenomena in string theory, such as axions, Chern-Simons terms, worldvolume degrees of freedom, and branes ending on or dissolving in other branes, can be interpreted as consequences of the absence of Chern-Weil symmetries in quantum gravity, suggesting that they might be general features of quantum gravity. We further discuss implications of breaking and gauging Chern-Weil symmetries for particle phenomenology and for boundary CFTs of AdS bulk theories. Chern-Weil global symmetries thus offer a unified framework for understanding many familiar aspects of quantum field theory and quantum gravity.

Brane dynamics in gauge theories↗
DOE OSTI2021

Resolving the ($\mathscr{g}$-2)$_μ$ discrepancy with $\mathcal{F}$–SU(5) intersecting D-branes

A discrepancy between the measured anomalous magnetic moment of the muon (g - 2) μ and computed Standard Model value now stands at a combined 4.2σ following experiments at Brookhaven National Lab (BNL) and the Fermi National Accelerator Laboratory (FNAL). A solution to the disagreement is uncovered in flipped SU(5) with additional TeV-Scale vector-like 10 +$\bar{10}$ multiplets and charged singlet derived from local F-Theory, collectively referred to as $\mathcal {F}$-SU(5). Here we engage general No-Scale supersymmetry (SUSY) breaking in $\mathcal {F}$-SU(5) D-brane model building to alleviate the (g-2)μ. tension between the Standard Model and observations. A robust Δa μ (SUSY) is realized via mixing of M 5 and M 1X at the secondary SU(5) x U(1) X unification scale in $\mathcal {F}$-SU(5) emanating from SU(5) breaking and U(1) X flux effects. Calculations unveil Δa μ (SUSY) = 19.0-22.3 x 10 -10 for gluino masses of M($\tilde{g}$ = 2.25-2.56 TeV nd higgsino dark matter, aptly residing within the BNL+FNAL 1σ mean. This (g – 2)μ favorable region of the model space also generates the correct light Higgs boson mass and branching ratios of companion rare decay processes, and is further consistent with all LHC Run 2 constraints. Finally, we also examine the heavy SUSY Higgs boson in light of recent LHC searches for an extended Higgs sector.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

R-parity violating decays of Bino neutralino LSPs at the LHC

The R-parity violating decays of Bino neutralino LSPs are analyzed within the context of the B - L MSSM “heterotic standard model”. These LSPs correspond to statistically determined initial soft supersymmetry breaking parameters which, when evolved using the renormalization group equations, lead to an effective theory satisfying all phenomenological requirements; including the observed electroweak vector boson masses and the Higgs mass. The explicit RPV decay channels of these LSPs into standard model particles, the analytic and numerical decay rates and the associated branching ratios are presented. The analysis of these quantities breaks into two separate calculations; first, for Bino neutralino LSPs with mass larger than MW± and, second, when the Bino neutralino mass is smaller than the electroweak scale. The RPV decay processes in both of these regions is analyzed in detail. The decay lengths of these RPV interactions are discussed. It is shown that for heavy Bino neutralino LSPs the vast majority of these decays are “prompt”, although a small, but calculable, number correspond to “displaced” decays of various lengths. The situation is reversed for light Bino LSPs, only a small number of which can RPV decay promptly. The relation of these results to the neutrino hierarchy–either normal or inverted–is discussed in detail.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A stochastic solver based on the residence time algorithm for crystal plasticity models

Abstract The deformation of crystalline materials by dislocation motion takes place in discrete amounts determined by the Burgers vector. Dislocations may move individually or in bundles, potentially giving rise to intermittent slip. This confers plastic deformation with a certain degree of variability that can be interpreted as being caused by stochastic fluctuations in dislocation behavior. However, crystal plasticity (CP) models are almost always formulated in a continuum sense, assuming that fluctuations average out over large material volumes and/or cancel out due to multi-slip contributions. Nevertheless, plastic fluctuations are known to be important in confined volumes at or below the micron scale, at high temperatures, and under low strain rate/stress deformation conditions. Here, we develop a stochastic solver for CP models based on the residence-time algorithm that naturally captures plastic fluctuations by sampling among the set of active slip systems in the crystal. The method solves the evolution equations of explicit CP formulations, which are recast as stochastic ordinary differential equations and integrated discretely in time. The stochastic CP model is numerically stable by design and naturally breaks the symmetry of plastic slip by sampling among the active plastic shear rates with the correct probability. This can lead to phenomena such as intermittent slip or plastic localization without adding external symmetry-breaking operations to the model. The method is applied to body-centered cubic tungsten single crystals under a variety of temperatures, loading orientations, and imposed strain rates.

42 ENGINEERING↗
DOE OSTI2021

Mechanical Performance and Mesostructure Analysis of Proton-Irradiated Fused Filament Fabrication Acrylonitrile Butadiene Styrene Material

The use of fused filament fabrication (FFF) acrylonitrile butadiene styrene (ABS) and other thermoplastics in radiation environments is beginning to be studied as possible replacements for traditionally manufactured parts and tools. Interlayer adhesion within the mesostructure has been shown in published literature to be an integral component in the strength of FFF ABS. Research of irradiated 3D printed polymers has primarily focused on the influence of gamma irradiation on the mechanical properties of FFF ABS samples, without evaluating its impact on the mesostructure of the samples. The purpose of this paper is to understand the mechanical damage caused by proton radiation on FFF ABS samples through the evaluation of the mesostructure of the ABS samples. To achieve this objective, proton radiation at 40 MeV was applied to FFF ABS samples at radiation doses up to 1.0 MGy at a dose rate of 1 MGy/hr. Following the irradiation, tensile testing was performed on the samples. Here, the fractured surfaces of the tested samples were subsequently observed using a scanning electron microscope. The mechanical test results show an increase in the ultimate tensile strength (UTS) and a decrease in the ductility of the irradiated samples. Statistical analysis on the results shows that there is a statistically significant difference in the UTS and the ductility of unirradiated and 1.0 MGy irradiated samples and between irradiated samples and 1.0 MGy samples. In addition, the percentage elongation at break has statistically significant differences in the means between irradiated samples and 1.0 MGy samples. The UTS has statistically significant differences in the means between 0.1 MGy and 1.0 MGy. The difference in the means between the unirradiated and the 1.0 MGy samples is the most significant for both the UTS and elongation at break. The scanning electron microscopy (SEM) results indicated that interlayer adhesion improved as a function of radiation dose corresponding to the increase in tensile strength. The SEM results also showed that crazing and plastic deformation were reduced; aligning with the loss in ductility observed in the tensile tests results. The proton radiation is causing these mechanical and physical changes through two mechanisms: (1) dose rate effects on ionizing radiation-induced oxidative degradation; and (2) radiation heating effects of high energy (>1 MeV) charged particles.

electron microscopy↗
DOE OSTI2024

Tension-compression asymmetry in superelasticity of SrNi 2 P 2 single crystals and the influence of low temperatures

ThCr 2 Si 2 -type intermetallic compounds are known to exhibit superelasticity associated with structural transitions through lattice collapse and expansion. These transitions occur via the formation and breaking of Si-type bonds, respectively, under uniaxial loading along the [0 0 1] direction. Unlike most ThCr 2 Si 2 -type intermetallic compounds, which have either an uncollapsed tetragonal structure or a collapsed tetragonal structure, SrNi 2 P 2 possesses a third type of collapsed structured: a one-third orthorhombic structure, for which one expects the occurrence of unique structural transitions and superelastic behavior. In this study, uniaxial compression and tension tests were conducted on micron-sized SrNi 2 P 2 single crystalline columns at room temperature, 200 K, and 100 K, to investigate the influence of loading direction and temperature on the superelasticity of SrNi 2 P 2 . Experimental data and density functional theory calculations revealed the presence of tension-compression asymmetry in the structural transitions and superelasticity, as well as an asymmetry in their temperature dependence, due to the opposite superelastic process associated with compression (forming P-P bonds) and tension (breaking P-P bonds). Additionally, following thermodynamics, the observations suggest that this asymmetric superelasticity could lead to an opposite elastocaloric effect between compression and tension, which could be beneficial potentially in obtaining large temperature changes compared to conventional superelastic solids that show the same elastocaloric effect regardless of loading direction. Furthermore, these results provide an important fundamental insight into the structural transitions, superelasticity processes, and potential elastocaloric effects in SrNi 2 P 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2022

Reverse water-gas shift: Na doping of m-ZrO 2 supported Pt for selectivity control

Reverse water-gas shift (RWGS) is a vital step in producing syngas for the chemical conversion of CO 2 to liquid transportation fuels and chemicals. Na-doping of m-ZrO 2 supported Pt catalysts allowed selectivity control by systematically increasing the ratio of relative rates of r CO /r CH4 . This was achieved by facilitating the formation of formate intermediate species, which precedes CO formation, and by suppressing the metallic Pt 0 active sites responsible for CH 4 formation. Here, a 2.5%Na-2%Pt/m-ZrO 2 catalyst was first tested for the forward water-gas shift (FWGS) reaction and found to have 50% higher CO conversion at 285°C compared to the undoped catalyst. Results of DRIFTS spectroscopy of adsorbed CO confirmed a formate ν(CH) band shift to lower wavenumbers (2870–2802 cm —1 ) with the addition of Na and more rapid forward formate decomposition in steam to H 2 and carbonate species, the precursor to CO 2 . This is consistent with C-H bond breaking being the rate limiting step of a FWGS mechanism occurring at the metal-support junction. Consistent with this, DRIFTS of RWGS in 4%CO 2 + 60%H 2 showed more facile formation of formate for the Na-doped catalyst and, once again, the ν(CH) band was shifted to lower wavenumbers (2874–2803 cm —1 ) with Na-doping. In addition, Na doping resulted in a systematic decrease in the Pt-carbonyl band in DRIFTS of adsorbed CO as well as DRIFTS of in-situ RWGS reaction tests, suggesting that Na blocked a fraction of on-top Pt sites, breaking up ensembles of Pt 0 responsible for methanation. The selectivity of the 2.5%Na-doped catalyst, unlike its undoped counterpart, was remarkably resistant to methanation (e.g., selectivity < 0.2% CH 4 with pressures of up to 20 bar).

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

Oxide layer delamination: An energy dissipation mechanism during high-velocity microparticle impacts

The break-up of microparticle surface oxides is needed to generate a clean metal-metal contact for metallurgical bonding during impact events such as seen in cold-spray. However, the mechanism by which this occurs is not clear. Here, using a laser-induced particle impact tester, we conduct site-specific experiment of single Cu microparticle impacts on a polished Cu substrate and characterize the impact sites. We observe that oxide layers on the particles begin to delaminate at velocities high enough to cause jetting of the substrate around the periphery of the impact site. We show that the energy cost of such delamination is ~30% of the energy needed to slow and stop particle as it bonds to the substrate, with the remainder being associated with the jetting process and bonding itself. At higher velocities, the oxide barrier is broken up to permit metallurgical bonding. Closer to the particle south-pole where hoop strain is low, oxide layers break up with few gaps and thus no metallurgical bonding is observed. Away from the south-pole where larger strains develop, oxide breakup permits intermittent bonding by the extrusion of bare metal into the gaps between oxide islands. These observations provide new insights towards understanding impact-induced metallurgical bonding during cold-spray process.

36 MATERIALS SCIENCE↗
DOE OSTI2020

On The Use of Sectional Techniques for the Solution of Depolymerization Population Balances: Results on a Discrete-Continuous Mesh

To study the discrete bond-breaking phenomena of depolymerization, the use of a fully continuous Population Balance Equation (PBE) is inadequate to embody all the inherent characteristics of the process, thus resulting in the need for a discrete-continuous mesh. In this work, the performance of the three most state-of-the-art sectional techniques, i.e. the fixed pivot technique (FPT), cell average technique (CAT) and finite volume scheme (FVS) in approximating discrete depolymerization using discrete-continuous PBEs was extensively compared and evaluated. The solutions from these three methods show different accuracy depending on the breakage mechanisms. For chain-end scission, the FPT and the CAT satisfactorily predict the population densities and moments whereas the FVS fails to predict the population densities but preserves the zeroth and the first moments. In the application of a discrete-continuous model, we identified a previously-not-reported issue of a precipitous drop in the number density at the boundary of discrete and continuous region specifically for chain-end scission. We successfully fixed this problem by employing the alterations proposed in this paper, to the particle allocation functions at the boundary points. For random scission, all three sectional techniques predict the population densities and moments to a high degree of accuracy, even at a very coarse mesh, through the use of our new stoichiometric kernel which is able to closely approximate the inherently discrete bond-breaking depolymerization process. The assessments in this present work intends to provide a clear-cut direction to efficient and economical modelling of depolymerization processes.

Ahamed, Firnaaz↗
DOE OSTI2020

Ring-opening and hydrodenitrogenation of indole under hydrothermal conditions over Ni, Pt, Ru, and Ni-Ru bimetallic catalysts

Here, the activity and selectivity of activated-carbon-supported Ni, Pt, Ru, and Ni-Ru bimetallic catalysts was examined for hydrothermal denitrogenation of indole. The molar yield of pyrrole ring-opening compounds, without an added hydrogen source, are in the order: Ni < Pt < Ni 90 Ru 10 < Ni 75 Ru 25 < Ni 50 Ru 50 ~ Ni 25 Ru 75 ~ Ru. Ru-containing catalysts facilitated production of hydrocarbons (hydrodenitrogenation (HDN) products) when used with added formic acid (hydrogen source). We elucidated catalytic hydrothermal HDN pathways for indole based on the product distributions and the variation of their yields with time. Hydrogenation of indole to indoline is the primary pathway and ring-opening of indoline to form alkyl anilines is faster than forming HDN products (alkyl benzenes). DFT calculations confirmed experimental activity trends, showing Ru is more active than Ni for indole ring opening and for o-toluidine deamination. If no hydrogen source is present, directly breaking the N-C bond in the pyrrole ring is more favorable than breaking the C-N bond with an aromatic carbon. If a H source is provided, the pyrrole ring hydrogenates first, forming indoline, followed by cleavage of the C-N bond.

42 ENGINEERING↗
DOE OSTI2021

Plastic deformation mechanism of calcium-silicate hydrates determined by deviatoric-stress Raman spectroscopy

Highlights: • C-(A-)S-H samples under stress were characterized at the bond and grain level. • Small angle scattering showed that 720 MPa reduced grain thickness by 30%. • In-situ Raman showed that cross-linked C-A-S-H slips along its intralayer. • However, in-situ Raman showed that C-S-H slips along the interlayer. • In-situ Raman also showed a growth in CaCO{sub 3}, indicating fracture in both samples. Creep of the cement matrix affects the structural stability of concrete. In Portland cements, the creep is largely controlled by the binding phase calcium-(aluminum-)silicate-hydrate, or C-(A-)S-H. This phase has a lamellar structure and under deviatoric stress aligns its c-axis with the principal stress. However, the limiting resistance to this reorientation is unknown at the nanocrystalline level. Small-angle X-ray scattering shows that the lamellae thickness decreases under 100's MPa deviatoric stress. Deviatoric stress Raman spectroscopy shows that there are two ways that this break-up can occur. If the material's silicate chains are cross-linked, then strain in SiO bonds does not increase above certain stresses, indicating a relaxation adjacent to the SiO bond. If the chains are not cross-linked, then the silicate chains are broken up by rastering against each other, introducing defects. These results show that the plastic deformation of C-(A-)S-H is relevant for Portland cement creep.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A machine-learning framework for the simulation of nuclear deflection of Planet-Killer-Asteroids

Here, as detection capabilities in astronomy have dramatically improved over the last two decades, concerns over Planet-Killer-Asteroids (PKAs) have become widespread, with nuclear weapons being proposed to destroy or deflect asteroids that are on a short-term projected collision course with Earth. Two main mitigation strategies have been proposed: • Case 1: Break up an incoming asteroid into smaller pieces that will disperse widely, resulting in smaller-scale, less detrimental, Earth-impacts or • Case 2: Deflect an incoming asteroid trajectory to avoid collision altogether. While the two strategies are not mutually exclusive, deflection is a safer strategy, ideally by harnessing all of the released energy from a nuclear device to move the asteroid as a rigid body. However, this case may not be always possible, since the strength of the energy release may break up the asteroid. In this work, the dynamical response of a PKA to a series of ultra-high energy impulses, such as those generated by nuclear devices, is formulated. A rapid iterative Discrete Element Method (DEM) method is developed to describe the deflection and potential breakup of the PKA as a function of a material bonding strength parameter within the asteroid and the magnitude of the applied impulse. The use of DEM allows for fragmentation of the PKA and the ability to compute the trajectories and distribution of the resulting debris field. Finally, a machine-learning algorithm is then developed and combined with the DEM approach to optimize the pulsation strategy for maximum possible safety and success.

42 ENGINEERING↗
DOE OSTI2022

The effect of elastic anisotropy on the symmetry selection of irradiation-induced void superlattices in cubic metals

We report self-organized microstructures and patterns have been widely observed in non-equilibrium physical systems. In particular, irradiation in metals creates far-from-equilibrium environments, in which the competing dynamics of defect production and annihilation can lead to unique self-organized superlattice structures, e.g., void and gas bubble superlattices. From a physical point of view, the superlattice structures are dictated by the intrinsic symmetry breaking in the metals, i.e., anisotropy caused by the breaking of continuous rotational symmetry. In the literature, two distinctive anisotropies, elastic anisotropy and diffusion anisotropy of interstitials, have been proposed to be the origins of superlattice formation. However, it is still unclear which anisotropy dominates the symmetry selection of superlattice structures. In this paper, we study elastic anisotropy and its effect on the symmetry of void superlattices. By using theoretical analyses and phase field simulations, we show that elastic anisotropy in cubic metals can lead to either face-centered cubic or simple cubic superlattices depending on the Zener anisotropy ratio. The superlattices formed under this elastic anisotropy mechanism must form under the influence of spinodal decomposition, as the mechanism requires perturbations in the vacancy concentration field to develop into spatially-static concentration waves. We compare to existing work on symmetry selection in superlattices via diffusion anisotropy and to experimental observations, and we suggest that concentration wave development under the influence of elastic anisotropy is not the mechanism for symmetry selection during the formation of irradiation-induced void superlattices, but that diffusion anisotropy could be the dominant mechanism.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Time-discretization of a plasma-neutral MHD model with a semi-implicit leapfrog algorithm

The semi-implicit leapfrog time-discretization is a workhorse algorithm for initial-value MHD codes to bridge between vastly separated time scales. Inclusion of atomic interactions with neutrals breaks the functional structure of the MHD equations that exploited by the leapfrog. In this work, we address how to best integrate atomic physics into the semi-implicit leapfrog. Following the Crank-Nicolson method, one approach is to time-center the atomic interactions in the linear solver and use a Newton method to include the nonlinear contributions. Alternatively, another family of methods are based on operator-splitting the terms associated with the atomic interactions using a Strang-splitting technique. These methods naturally break equations into constituent ODE and PDE parts and preserve the structure exploited by the semi-implicit leapfrog. We study the accuracy and efficiency of these methods through a battery of 0D and 1D cases and show that a second-order-in-time Douglas-Rachford inspired coupling between the ODE and PDE advances is effective in reducing the time-discretization error to be comparable to that of Crank-Nicolson with Newton iteration of the nonlinear terms. Splitting ODE and PDE parts results in independent matrix solves for each field which reduces the computational cost considerably and provides parallelization over species relative to Crank-Nicolson.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗