Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Magnetic source”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Charm contribution to ultrahigh-energy neutrinos from newborn magnetars

Newborn, strongly magnetized neutron stars (so-called magnetars) surrounded by their stellar or merger ejecta are expected to be sources of ultrahigh-energy neutrinos via decay of mesons produced in hadronic interactions of protons which are accelerated to ultrahigh energies by magnetic dissipation of the spindown energy. In this work, we show that not only pions and kaons but also charm hadrons, which are typically neglected due to their small production cross sections, can represent dominant contributions to neutrino fluence at ultrahigh energies, because of their short lifetimes, while the ultrahigh-energy neutrino fluence from pion and kaon production is suppressed at early times due to their significant cooling before their decay. We show that the next-generation detectors such as Probe Of Extreme Multi-Messenger Astrophysics (POEMMA), Giant Radio Array for Neutrino Detection (GRAND) and IceCube-Gen2 have a good chance of observing neutrinos, primarily originating from charm hadrons, from nearby magnetars. We also show that neutrinos from nearby magnetar-driven merger novae could be observed in the time interval between 10 2 s and 10 3 s, where the charm hadron contribution is dominant for neutrino energies above 10 8 GeV, of relevance to next generation detectors. We also comment on potential impacts of the charm hadron contribution to the diffuse neutrino flux.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Trace impurity transport in multi-species plasmas with large particle fluxes

A burning fusion plasma has a large inwards flux of fuel and outwards flux of ash. Existing impurity transport theories do not account for a steady-state environment of such large fluxes. In this paper, we extend classical transport theory to account for a background of such fluxes. In a mainly two-ion species magnetized plasma, with oppositely directed density gradients maintained by oppositely directed sources of the two species, trace impurity ions are shown to align with the larger mass ions. This unexpected behavior is derived analytically and simulated numerically using the MITNS (Multiple-Ion Transport Numerical Solver) code. Overall, the result suggests a potential advantage of edge-fueled p-B11 fusion in extracting high-Z impurities, and possible applications in plasma-based separation methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulation-trained machine learning models for Lorentz transmission electron microscopy

Understanding the collective behavior of complex spin textures, such as lattices of magnetic skyrmions, is of fundamental importance for exploring and controlling the emergent ordering of these spin textures and inducing phase transitions. It is also critical to understand the skyrmion–skyrmion interactions for applications such as magnetic skyrmion-enabled reservoir or neuromorphic computing. Magnetic skyrmion lattices can be studied using in situ Lorentz transmission electron microscopy (LTEM), but quantitative and statistically robust analysis of the skyrmion lattices from LTEM images can be difficult. In this work, we show that a convolutional neural network, trained on simulated data, can be applied to perform segmentation of spin textures and to extract quantitative data, such as spin texture size and location, from experimental LTEM images, which cannot be obtained manually. This includes quantitative information about skyrmion size, position, and shape, which can, in turn, be used to calculate skyrmion–skyrmion interactions and lattice ordering. We apply this approach to segmenting images of Néel skyrmion lattices so that we can accurately identify skyrmion size and deformation in both dense and sparse lattices. The model is trained using a large set of micromagnetic simulations as well as simulated LTEM images. This entirely open-source training pipeline can be applied to a wide variety of magnetic features and materials, enabling large-scale statistical studies of spin textures using LTEM.

McCray, Arthur R. C. (ORCID:0000000160774698)↗
DOE OSTI2022

Investigation of resistive magnetic field generation by intense proton beams in dense plasmas

Current and future applications of intense proton sources abound, including radiography, cancer therapy, warm dense matter generation, and inertial confinement fusion. With increasingly efficient acceleration and focusing mechanisms, proton current densities may soon approach and exceed 10 10 A/cm 2 , e.g., via intense laser drivers. Simulations have previously shown that in this current density regime, beam-induced field generation plays a significant role in beam transport through dense plasmas. Here, we present a theoretical model for the generation of resistive magnetic fields by intense proton beam transport through solid density plasmas. The theoretical evolution of the magnetic field profile is calculated using an analytic model for aluminum resistivity, heat capacity, and stopping power, applicable from cold matter to hot plasma. The effects of various beam and material parameters on the field are investigated and explained for both monoenergetic and Maxwellian proton beams. For a proton beam with Maxwellian temperature 5 MeV and total energy 10 J, the model calculates resistive magnetic fields up to 150 T in aluminum. In conclusion, the calculated field profiles from several beam cases are compared with 2D hybrid particle-in-cell simulations, with good agreement found in magnitude and time scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Understanding Complex Magnetic Spin Textures with Simulation-Assisted Lorentz Transmission Electron Microscopy

There is an increased interest in topologically nontrivial magnetic spin textures such as skyrmions and chiral domain-wall solitons, both from a point of fundamental physics understanding as well as potential technological interest in low-power memory applications. In order to control their behavior, it is necessary to understand their complex spin texture at the nanoscale. Lorentz transmission electron microscopy (LTEM) is a suitable technique for studying these systems due to its high spatial resolution and capability to simultaneously characterize magnetic texture and microstructure. In this work, we present the application of PyLorentz, an open-source software suite that we have developed, for quantitative image analysis of Neel-type skyrmions in thin-film heterostructures. PyLorentz enhances LTEM capabilities by enabling reconstruction of magnetic induction maps from experimental images, as well as simulating LTEM images using micromagnetic simulation data. We demonstrate this for simulated Neel skyrmions as well as experimental data from [Pt/Co/W] multilayer heterostructures. Finally, we also show how simulation-assisted LTEM analysis is crucial for understanding these complex magnetic spin textures, in which the reconstructed magnetic induction map (seen in the LTEM images) differs significantly from the magnetization configuration.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

ParaStell: parametric modeling and neutronics support for stellarator fusion power plants

The three-dimensional variation inherent to stellarator geometries and fusion sources motivates three-dimensional modeling to obtain accurate results from computational modeling in support of design and analysis of first wall, blanket, and shield (FWBS) systems. Manually constructing stellarator fusion power plant geometries in computer-aided design (CAD) and defining the corresponding fusion source can be cumbersome and challenging. The open-source parametric modeling toolset ParaStell has been developed to automate construction of such geometries in low-fidelity. Low-fidelity modeling is useful during the conceptual phase of engineering design as a means of rapidly exploring the design space of a given device. The modeling capability of ParaStell includes in-vessel components and magnets, for any given stellarator configuration, using a parametric definition and plasma equilibrium data. Furthermore, the toolset automates the generation of detailed, tetrahedral neutron source definitions and DAGMC geometries for use in neutronics modeling. ParaStell assists rapid design iteration, parametric study, and design optimization of stellarator fusion cores. As a demonstration of the design iteration capability, the effect of the three-dimensional parameter space on tritium breeding and magnet shielding is investigated, using the WISTELL-D configuration as a design basis. Blanket and shield thicknesses are varied in three dimensions, using the space available between the plasma edge and magnet coils as a constraint. The corresponding effects on tritium breeding ratio and magnet heating are tallied using the open-source Monte Carlo particle transport code OpenMC. The inclusion of additional and higher-fidelity modeling capabilities is planned for ParaStell’s future, as well as its implementation in machine-driven optimization.

Moreno, Connor↗
DOE OSTI2024

A survey of nuclear quadrupole deformation in order to estimate the nuclear MQM and its relative contribution to the atomic EDM: Searching for the best candidate to focus on an atomic EDM measurement: A nuclear physics perspective

Abstract Measurement of a non-zero permanent electric dipole moment (EDM) in fundamental particles, such as in an electron or a neutron, or in nuclei or atoms, can help us gain a handle on the sources of Charge-Parity (CP) violation, both in the Standard Model (SM) and beyond. The nuclear magnetic quadrupole moment (MQM), the central topic of this work, is also CP, P, and T violating. Nucleons and nuclei have a non-zero MQM from sources within the SM, but the nuclear MQM is dramatically enhanced if the nuclei are structurally quadrupole deformed. Multiple sources contribute to an atomic EDM namely: (i) nuclear EDM through its Schiff moment, which is enhanced by nuclear octupole deformation, (ii) CP violating interactions between the electrons and the nuclei, and (iii) the nuclear MQM that contributes to the atomic EDM in atoms with an unpaired valence electron. Our survey of nuclear quadrupole deformation identified 48 isotopes as ideal systems in which to search for a CP violating EDM via their enhanced nuclear MQM. Of these candidates, $$^{223,225}$$ 223 , 225 Fr, $$^{223}$$ 223 Ra, $$^{223,225,227}$$ 223 , 225 , 227 Ac, $$^{229}$$ 229 Th, and $$^{229}$$ 229 Pa also have maximally enhanced nuclear Schiff moment contribution due to their octupole deformation. Laser cooling of the isotopes of Fr and Ra, among a few others, has already been demonstrated, making $$^{223,225}$$ 223 , 225 Fr and $$^{223}$$ 223 Ra some of the best systems in which to measure an EDM.

MohanMurthy, Prajwal↗
DOE OSTI2020

Effect of magnetic field ripple on parallel electron transport during microwave plasma heating in the Proto-MPEX linear plasma device

The Prototype Material Plasma Exposure eXperiment (Proto-MPEX) is a linear plasma device located at Oak Ridge National Laboratory to develop the plasma source concept for the Materials Plasma Exposure eXperiment (MPEX). Recent experiments have demonstrated the heating of electrons in a high-density deuterium plasma (n e ~ 5 × 10 19 m -3 ) via electron Bernstein waves and source plasma production via helicon waves in this linear configuration. Moreover, experimental observations suggest that the magnetic ripple adversely affects the parallel transport of the heated electrons toward the target where material samples are to be exposed to the plasma. In order to understand the transport process during microwave application, a test-particle Monte-Carlo (MC) code has been developed that incorporates the effects of Coulomb collisions, magnetic mirror adiabatic trapping and electron cyclotron interaction via a quasilinear RF heating operator. MC calculations indicate that the absence of 2nd turning points along the trajectory of cyclotron heated electrons significantly reduces adiabatic trapping and maximizes power transport to the Target surface. Thermalization of cyclotron heating electrons is analyzed with the MC code and it is found that for conditions relevant to ProtoMPEX, a significant fraction of the absorbed heating power is coupled to the Target surface via ballistic fast electrons; however, under certain conditions, up to 60% of the absorbed RF power can be dissipated via thermalization of the fast electrons on the background plasma. Results are compared with experimental measurements. Finally, the implications of the results are discussed in the context of the upcoming MPEX device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Magnetic ac susceptibility of superconducting Ta films for quantum computing

Recently, breakthroughs were reported in the long coherence time (>0.5 ms) transmon qubit using tantalum films (Ta). Identifying the loss mechanism in Ta films will help further improve the performance of Ta-based superconducting qubits. We report a study of superconducting properties in two sets of Ta films grown on c-cut and a-cut sapphires respectively using contact transport and noncontact magnetic ac susceptibility measurements. Although the resistive transition appears to be similar in both films, we found strikingly different responses in their complex magnetic susceptibilities χ' + iχ'$^{'}$ .χ'$^{'}$ in the c-cut films exhibits a sharp peak at superconducting transition Tc and becomes featureless below Tc indicating a strongly coupled superconducting state. In contrast, χ'$^{'}$ \ in the a-cut films exhibits a broad peak near Tc and a second peak appears below Tc, indicating granular superconductivity behavior. This second peak χ'$^{'}$ in is associated with the hysteresis loss at weak-link-type grain boundaries, which is believed to be a leading source of decoherence in the a-cut Ta films. The derived magnetic loss tangent in the a-cut films is about one order of magnitude higher than in the c-cut films. This study demonstrates that ac susceptibility is an excellent probe for characterizing bulk superconducting properties of the constituent superconducting materials used in qubits.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

3D modeling of deep borehole electromagnetic measurements with energized casing source for fracture mapping at the Utah Frontier Observatory for Research in Geothermal Energy

Here, we present a 3D numerical modelling analysis evaluating the deployment of a borehole electromagnetic measurement tool to detect and image a stimulated zone at the Utah Frontier Observatory for Research in Geothermal Energy geothermal site. As the depth to the geothermal reservoir is several kilometres and the size of the stimulated zone is limited to several 100 m, surface-based controlled-source electromagnetic measurements lack the sensitivity for detecting changes in electrical resistivity caused by the stimulation. To overcome the limitation, the study evaluates the feasibility of using a three-component borehole magnetic receiver system at the Frontier Observatory for Research in Geothermal Energy site. To provide sufficient currents inside and around the enhanced geothermal reservoir, we use an injection well as an energized casing source. To efficiently simulate energizing the injection well in a realistic 3D resistivity model, we introduce a novel modelling workflow that leverages the strengths of both 3D cylindrical-mesh-based electromagnetic modelling code and 3D tetrahedral-mesh-based electromagnetic modelling code. The former is particularly well-suited for modelling hollow cylindrical objects like casings, whereas the latter excels at representing more complex 3D geological structures. In this workflow, our initial step involves computing current densities along a vertical steel-cased well using a 3D cylindrical electromagnetic modelling code. Subsequently, we distribute a series of equivalent current sources along the well's trajectory within a complex 3D resistivity model. We then discretize this model using a tetrahedral mesh and simulate the borehole electromagnetic responses excited by the casing source using a 3D finite-element electromagnetic code. This multi-step approach enables us to simulate 3D casing source electromagnetic responses within a complex 3D resistivity model, without the need for explicit discretization of the well using an excessive number of fine cells. We discuss the applicability and limitations of this proposed workflow within an electromagnetic modelling scenario where an energized well is deviated, such as at the Frontier Observatory for Research in Geothermal Energy site. Using the workflow, we demonstrate that the combined use of the energized casing source and the borehole electromagnetic receiver system offer measurable magnetic field amplitudes and sensitivity to the deep localized stimulated zone. The measurements can also distinguish between parallel-fracture anisotropic reservoirs and isotropic cases, providing valuable insights into the fracture system of the stimulated zone. Besides the magnetic field measurements, vertical electric field measurements in the open well sections are also highly sensitive to the stimulated zone and can be used as additional data for detecting and imaging the target. We can also acquire additional multiple-source data by grounding the surface electrode at various locations and repeating borehole electromagnetic measurements. This approach can increase the number of monitoring data by several factors, providing a more comprehensive dataset for analysing the deep-localized stimulated zone. The numerical analysis indicates that it is feasible to use the combination of the energized casing and downhole electromagnetic measurements in monitoring localized stimulated zone at large depths.

58 GEOSCIENCES↗
DOE OSTI2021

Direct evidence of microstructure dependence of magnetic flux trapping in niobium

Abstract Elemental type-II superconducting niobium is the material of choice for superconducting radiofrequency cavities used in modern particle accelerators, light sources, detectors, sensors, and quantum computing architecture. An essential challenge to increasing energy efficiency in rf applications is the power dissipation due to residual magnetic field that is trapped during the cool down process due to incomplete magnetic field expulsion. New SRF cavity processing recipes that use surface doping techniques have significantly increased their cryogenic efficiency. However, the performance of SRF Nb accelerators still shows vulnerability to a trapped magnetic field. In this manuscript, we report the observation of a direct link between flux trapping and incomplete flux expulsion with spatial variations in microstructure within the niobium. Fine-grain recrystallized microstructure with an average grain size of 10–50 µm leads to flux trapping even with a lack of dislocation structures in grain interiors. Larger grain sizes beyond 100–400 µm do not lead to preferential flux trapping, as observed directly by magneto-optical imaging. While local magnetic flux variations imaged by magneto-optics provide clarity on a microstructure level, bulk variations are also indicated by variations in pinning force curves with sequential heat treatment studies. The key results indicate that complete control of the niobium microstructure will help produce higher performance superconducting resonators with reduced rf losses 1 related to the magnetic flux trapping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Study of Magnetic Field and Turbulence in the TeV Halo around the Monogem Pulsar

Magnetic fields are ubiquitous in the interstellar medium, including extended objects such as supernova remnants and diffuse halos around pulsars. Its turbulent characteristics govern the diffusion of cosmic rays and the multiwavelength emission from pulsar wind nebulae (PWNe). However, the geometry and turbulence nature of the magnetic fields in the ambient region of PWN is still unknown. Recent gamma-ray observations from HAWC and synchrotron observations suggest a highly suppressed diffusion coefficient compared to the mean interstellar value. In this study, we present the first direct observational evidence of the orientation of the mean magnetic field and turbulent characteristics by employing a recently developed statistical parameter "Y turb " in the extended halo around the Monogem pulsar. Our study points to two possible scenarios: nearly aligned toward the line of sight (LOS) with compressible modes dominance or high inclination angle toward the LOS and characterized by Alfvénic turbulence. The first scenario appears consistent with other observational signatures. Furthermore, we report that the magnetic field has an observed correlation length of approximately 3 ± 0.6 pc in the Monogem halo. Our study highlights the pivotal role of magnetic field and turbulence in unraveling the physical processes in TeV halos and cosmic-ray transport.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Feasibility of FDEM cross-dipole and TDEM loop sources for monitoring CO 2 at the Kemper CarbonSAFE site

Electromagnetics has been shown to be a viable tool to monitor CO 2 plumes embedded in saline reservoirs. However, the majority of studies focus on measuring the electric field, which requires cumbersome equipment in the field and it is difficult to install permanent measurement stations. Magnetic field receivers offer an opportunity to reduce the form factor of the survey and increase the mobility by utilizing upcoming technologies, such as drones. We explore the use of frequency-domain electric dipole sources, and timedomain loops with a focus on measuring the secondary magnetic field at the surface for a conceptual injection scenario based on the Kemper CarbonSAFE site. We find that electric dipole sources give a response above the sensitivity of current sensor technology and, therefore, be a viable tool for CO 2 monitoring. The time domain loop source does provide fields that are useful for determining the location of the CO 2 plume, however the field magnitude is below the sensitivity of the current generation of instruments. To explore the use of a potential borehole receiver we generate a map of the magnetic field at depth to explore potential borehole placement for monitoring efforts. Finally, we limit the spatial extent of the electric dipole survey to a single parcel of land to help understand how the fields change with survey geometry. We find that the shape of the secondary fields change slightly with the small transmitter, but are still measurable provided that the cultural noise at the site is low. Thus, we conclude that at the Kemper site a frequency-domain cross-dipole source with magnetometer receivers is suitable to monitor the expansion of the CO 2 plume in the saline reservoir, even with a limited transmitter footprint on the surface.

58 GEOSCIENCES↗
DOE OSTI2023

Design of the Lanthanum hexaboride based plasma source for the large plasma device at UCLA

The Large Plasma Device (LAPD) at UCLA (University of California, Los Angeles) produces an 18 m long, magnetized, quiescent, and uniform plasma at a high repetition rate to enable studies of fundamental plasma physics. Here, we report on a major upgrade to the LAPD plasma source that allows for more robust operation and significant expansion of achievable plasma parameters. The original plasma source made use of a heated barium oxide (BaO) coated nickel sheet as an electron emitter. This source had a number of drawbacks, including a limited range of plasma density (≲4.0 × 10 12 cm −3 ), a limited discharge duration (∼10 ms), and susceptibility to poisoning following oxygen exposure. Further, the new plasma source utilizes a 38 cm diameter lanthanum hexaboride (LaB 6 ) cathode, which has a significantly higher emissivity, allowing for a much larger discharge power density, and is robust to exposure to air. Peak plasma density of up to 3.0 × 10 13 cm −33 in helium gas has been achieved. The typical operating pressure is ∼10 −5 Torr, while dynamic pressure can be achieved through the gas-puffing technique. Discharges as long as 70 ms have been produced, enabling a variety of long-time-scale studies of processes, such as turbulent particle transport. The new source has been in continuous operation for 14 months, having survived air leaks, power outages that led to rapid temperature changes on the cathode and heater, and planned machine openings. We describe the design, construction, and initial operation of this novel new large-area LaB6 plasma source.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Annihilation-gamma-based diagnostic techniques for magnetically confined electron–positron pair plasma

Efforts are underway to magnetically confine electron–positron pair plasmas to study their unique behaviour, which is characterized by significant changes in plasma time and length scales, supported waves and unstable modes. However, use of conventional plasma diagnostics presents challenges with these low-density and annihilating matter–antimatter plasmas. To address this problem, we propose to develop techniques based on the distinct emission provided by annihilation. This emission exhibits two spatial correlations: the distance attenuation of isotropic sources and the back-to-back propagation of momentum-preserving 2 $\gamma$ annihilation. We present the results of our analysis of the $\gamma$ emission rate and the spatial profile of the annihilation in a magnetized pair plasma from direct pair collisions, from the formation and decay of positronium as well as from transport processes. In order to demonstrate the effectiveness of annihilation-based techniques, we tested them on annular $\gamma$ emission profiles produced by a $\beta ^+$ radioisotope on a rotating turntable. Direct and positronium-mediated annihilation result in overlapping volumetric $\gamma$ sources, and the 2 $\gamma$ emission from these volumetric sources can be tomographically reconstructed from coincident counts in multiple detectors. Transport processes result in localized annihilation where field lines intersect walls, limiters or internal magnets. These localized sources can be identified by the fractional $\gamma$ counts on spatially distributed detectors.

Physics↗
DOE OSTI2022

The hard x-ray nanotomography microscope at the advanced light source

Beamline 11.3.1 at the Advanced Light Source is a tender/hard (6–17 keV) x-ray bend magnet beamline recently re-purposed with a new full-field, nanoscale transmission x-ray microscope. The microscope is designed to image composite and porous materials possessing a submicrometer structure and compositional heterogeneity that determine materials’ performance and geologic behavior. The theoretical and achieved resolutions are 55 and <100 nm, respectively. The microscope is used in tandem with a <25 nm eccentricity rotation stage for high-resolution volume imaging using nanoscale computed tomography. The system also features a novel bipolar illumination condenser for the illumination of an ~100 μm spot of interest on the sample, followed by a phase-type zone plate magnifying objective of ~52 µm field of view and a phase detection ring. The zone plate serves as the system objective and magnifies the sample with projection onto an indirect x-ray detection system, consisting of a polished single crystal CsI(Tl) scintillator and a range of high-quality Plan Fluorite visible light objectives. The objectives project the final visible light image onto a water-cooled CMOS 2048 × 2048-pixel 2 detector. Here, in this paper, we will discuss the salient features of this instrument and describe early results from imaging the internal three-dimensional microstructure and nanostructure of target materials, including fiber-reinforced composites and geomaterials.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

The Simons Observatory: science goals and forecasts for the enhanced Large Aperture Telescope

We describe updated scientific goals for the wide-field, millimeter-wave survey that will be produced by the Simons Observatory (SO). Significant upgrades to the 6-meter SO Large Aperture Telescope (LAT) are expected to be complete by 2028, and will include a doubled mapping speed with 30,000 new detectors and an automated data reduction pipeline. In addition, a new photovoltaic array will supply most of the observatory's power. The LAT survey will cover about 60% of the sky at a regular observing cadence, with five times the angular resolution and ten times the map depth of the Planck satellite. The science goals are to: (1) determine the physical conditions in the early universe and constrain the existence of new light particles; (2) measure the integrated distribution of mass, electron pressure, and electron momentum in the late-time universe, and, in combination with optical surveys, determine the neutrino mass and the effects of dark energy via tomographic measurements of the growth of structure at redshifts z ≲ 3; (3) measure the distribution of electron density and pressure around galaxy groups and clusters, and calibrate the effects of energy input from galaxy formation on the surrounding environment; (4) produce a sample of more than 30,000 galaxy clusters, and more than 100,000 extragalactic millimeter sources, including regularly sampled AGN light-curves, to study these sources and their emission physics; (5) measure the polarized emission from magnetically aligned dust grains in our Galaxy, to study the properties of dust and the role of magnetic fields in star formation; (6) constrain asteroid regoliths, search for Trans-Neptunian Objects, and either detect or eliminate large portions of the phase space in the search for Planet 9; and (7) provide a powerful new window into the transient universe on time scales of minutes to years, concurrent with observations from the Vera C. Rubin Observatory of overlapping sky.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Towards modelling AR Sco: generalized particle dynamics and strong radiation-reaction regimes

Numerical simulations of relativistic plasmas have become more feasible, popular, and crucial for various astrophysical sources with the availability of computational resources. The necessity for high-accuracy particle dynamics is especially highlighted in pulsar modelling due to the extreme associated electromagnetic fields and particle Lorentz factors. Including the radiation-reaction force in the particle dynamics adds even more complexity to the problem, but is crucial for such extreme astrophysical sources. We have also realized the need for such modelling concerning magnetic mirroring and particle injection models proposed for AR Sco, the first white dwarf pulsar. This paper demonstrates the benefits of using higher-order explicit numerical integrators with adaptive time-step methods to solve the full particle dynamics with radiation-reaction forces included. We show that for standard test scenarios, namely various combinations of uniform E- and B-fields and a static dipole B-field, the schemes we use are equivalent to and in extreme field cases outperform standard symplectic integrators in accuracy. We show that the higher-order schemes have massive computational time improvements due to the adaptive time-steps we implement, especially in non-uniform field scenarios and included radiation reaction where the particle gyro-radius rapidly changes. When balancing accuracy and computational time, we identified the adaptive Dormand–Prince eighth-order scheme to be ideal for our use cases. The schemes we use maintain accuracy and stability in describing the particle dynamics and we indicate how a charged particle enters radiation-reaction equilibrium and conforms to the analytical Aristotelian Electrodynamics expectations.

79 ASTRONOMY AND ASTROPHYSICS↗