Models Of Particle Acceleration Research Articles

The Intrabinary Shock and Companion Star of Redback Pulsar J2215+5135

PSR J2215+5135 (J2215) is a “redback” spider pulsar, where the intrabinary shock (IBS) wraps around the pulsar rather than the stellar-mass companion. Spider orbital light curves are modulated, dominated by their binary companion thermal emission in the optical bands and by IBS synchrotron emission in the X-rays. We report on new XMM-Newton X-ray and U-band observations of J2215. We produce orbital light curves and use them to model the system properties. Our best-fit optical light model gives a neutron star mass M NS = 1.98 ± 0.08 M ⊙, lower than previously reported. However, uncertainty in the stellar atmosphere metallicity, a parameter to which J2215 is unusually sensitive, requires us to consider an acceptable systematic plus statistical range of M NS ∼ 1.85–2.3 M ⊙. From the X-ray analysis, we find that the IBS wraps around the pulsar but with a pulsar-wind-to-companion-wind-momentum ratio unusually close to unity, implying a flatter IBS geometry than seen in other spiders. Estimating the companion wind momentum and speed from the X-ray light curve, we find a companion mass-loss rate of Ṁc≳10−10M⊙ yr−1 so that J2215 may become an isolated millisecond pulsar in ∼1 Gyr. Our X-ray analyses place constraints on the magnetization and particle density of the pulsar wind and support models of magnetic reconnection and particle acceleration in the highly magnetized relativistic IBS.

Nonthermal Signatures of Radiative Supernova Remnants

The end of supernova remnant (SNR) evolution is characterized by a so-called “radiative” stage, in which efficient cooling of the hot bubble inside the forward shock slows expansion, leading to eventual shock breakup. Understanding SNR evolution at this stage is vital for predicting feedback in galaxies, since SNRs are expected to deposit their energy and momentum into the interstellar medium at the ends of their lives. A key prediction of SNR evolutionary models is the formation at the onset of the radiative stage of a cold, dense shell behind the forward shock. However, searches for these shells via their neutral hydrogen emission have had limited success. We instead introduce an independent observational signal of shell formation arising from the interaction between nonthermal particles accelerated by the SNR forward shock (cosmic rays) and the dense shell. Using a semi-analytic model of particle acceleration based on state-of-the-art simulations coupled with a high-resolution hydrodynamic model of SNR evolution, we predict the nonthermal emission that arises from this interaction. We demonstrate that the onset of the radiative stage leads to nonthermal signatures from radio to gamma rays, including radio and gamma-ray brightening by nearly 2 orders of magnitude. Such a signature may be detectable with current instruments, and will be resolvable with the next generation of gamma-ray telescopes (namely, the Cherenkov Telescope Array).

Severe Electron Environment for Surface Charging at Geostationary and Medium Earth Orbits

Electron radiation environment with keV energies on medium earth orbit (MEO) is investigated around the times when worst-case severe environments for surface charging were detected at Los Alamos National Laboratory (LANL) satellites on geostationary earth orbit (GEO). The Inner Magnetosphere Particle Transport and Acceleration Model (IMPTAM) is employed to model 1–100 keV electron fluxes in the inner Earth’s magnetosphere (inside 10RE) during four representative events selected out of 400 LANL worst-case severe environment events. Modeled fluxes are compared to the observed ones on GEO, and the peak electron fluxes are then extracted at L=4.6 which is considered as approximate MEO environment following the LANL events. IMPTAM generally reproduces the 10–50 keV electron fluxes on GEO but underestimates (of factor of 2 to 5) at energies above. It is found that the maximum electron flux on MEO is reached in about 1.5–2 h after the worst-case environment for surface charging was detected at GEO. Whereas the magnetic local time (MLT) location of the worst-case severe environments on GEO varies from 23 to 05, the maximum electron flux on MEO is located on 6–7.8 MLT. The maximum electron flux values on MEO can be 2 to 10 times higher than those observed on GEO and those estimated using European Cooperation for Space Standardization and National Aeronautics and Space Administration guidelines, which can be a significant increase of the risks in terms of surface charging.

Solar Wind With Field Lines and Energetic Particles (SOFIE) Model: Application to Historical Solar Energetic Particle Events

AbstractIn this paper, we demonstrate the applicability of the data‐driven solar energetic particle (SEP) model, SOlar‐wind with FIeld‐lines and Energetic‐particles (SOFIE), to simulate the acceleration and transport processes of SEPs and make forecast of the energetic proton flux at energies ≥10 MeV that will be observed near 1 AU. The SOFIE model is built upon the Space Weather Modeling Framework developed at the University of Michigan. In SOFIE, the background solar wind plasma in the solar corona and interplanetary space is calculated by the Stream‐Aligned Aflvén Wave Solar‐atmosphere Model(‐Realtime) driven by the near‐real‐time hourly updated Global Oscillation Network Group solar magnetograms. In the background solar wind, coronal mass ejections (CMEs) are launched by placing an force‐imbalanced magnetic flux rope on top of the parent active region, using the Eruptive Event Generator using Gibson‐Low model. The acceleration and transport processes are modeled by the Multiple‐Field‐Line Advection Model for Particle Acceleration. In this work, nine SEP events (Solar Heliospheric and INterplanetary Environment challenge/campaign events) are modeled. The three modules in SOFIE are validated and evaluated by comparing with observations, including the steady‐state background solar wind properties, the white‐light image of the CMEs, and the flux of solar energetic protons, at energies of ≥10 MeV.

Testing particle acceleration in blazar jets with continuous high-cadence optical polarization observations

Variability can be the pathway to understanding the physical processes in astrophysical jets. However, the high-cadence observations required to test particle acceleration models are still missing. Here we report on the first attempt to produce continuous, $>24$ hour polarization light curves of blazars using telescopes distributed across the globe, following the rotation of the Earth, to avoid the rising Sun. Our campaign involved 16 telescopes in Asia, Europe, and North America. We observed BL Lacertae and CGRaBS J0211+1051 for a combined 685 telescope hours. We find large variations in the polarization degree and angle for both sources on sub-hour timescales as well as a $ rotation of the polarization angle in CGRaBS J0211+1051 in less than two days. We compared our high-cadence observations to particle-in-cell magnetic reconnection and turbulent plasma simulations. We find that although the state-of-the-art simulation frameworks can produce a large fraction of the polarization properties, they do not account for the entirety of the observed polarization behavior in blazar jets.

Abell 0399–Abell 0401 radio bridge spectral index: First multi-frequency detection

Aims. Recent low-frequency radio observations at 140 MHz discovered a bridge of diffuse emission with a length of 3 Mpc that connects the galaxy clusters Abell 0399 and Abell 0401. We present follow-up observations at 60 MHz to constrain the spectral index of the bridge, which has only been detected at 140 and 144 MHz so far. Methods. We analysed deep (∼18 h) LOw Frequency ARray (LOFAR) Low Band Antenna (LBA) data at 60 MHz to detect the bridge at very low frequencies. We then conducted a multi-frequency study with LOFAR HBA data at 144 MHz and uGMRT data at 400 MHz. Assuming second-order Fermi mechanisms for the re-acceleration of relativistic electrons driven by turbulence in the radio bridge regions, we compared the observed radio spectrum with theoretical synchrotron models. Results. The bridge is detected in the 75″ resolution LOFAR image at 60 MHz, and its emission fully connects the region between the two galaxy clusters. Between 60 MHz and 144 MHz, we found an integrated spectral index value of α60144 = −1.44 ± 0.16 for the bridge emission. For the first time, we produced spectral index and related uncertainties maps for a radio bridge. We produce a radio spectrum that shows a significant steepening between 144 and 400 MHz. Conclusions. This detection at low frequencies provides important information for models of particle acceleration and magnetic field structure on very extended scales. The spectral index gives important clues about the origin of inter-cluster diffuse emission. The steepening of the spectrum above 144 MHz can be explained in a turbulent re-acceleration framework, assuming that the acceleration timescales are longer than ∼200 Myr.

Ground response of rolling dynamic compaction—a finite element modelling approach

Rolling dynamic compaction (RDC) technology utilises a heavy non-circular module (impact roller) to compact the underlying soil dynamically. The stresses imparted to the soil through this technique and the resulting vibrations, have been the subject of investigation in the field. A finite element (FE) model predicting the settlement and densification of a coarse-grained fill material subject to RDC with a BH-1300 4-sided 8 tonne impact roller has been shown to provide good agreement with that observed in the field. This paper presents estimates using the developed FE model for the peak particle velocity and acceleration, and the maximum stresses applied through each impact upon a coarse-grained soil. Distributions of the results and their empirical formulae are presented herein.

The Large Imaging Spectrometer for Solar Accelerated Nuclei (LISSAN): A Next-Generation Solar γ-ray Spectroscopic Imaging Instrument Concept

Models of particle acceleration in solar eruptive events suggest that roughly equal energy may go into accelerating electrons and ions. However, while previous solar X-ray spectroscopic imagers have transformed our understanding of electron acceleration, only one resolved image of γ-ray emission from solar accelerated ions has ever been produced. This paper outlines a new satellite instrument concept—the large imaging spectrometer for solar accelerated nuclei (LISSAN)—with the capability not only to observe hundreds of events over its lifetime, but also to capture multiple images per event, thereby imaging the dynamics of solar accelerated ions for the first time. LISSAN provides spectroscopic imaging at photon energies of 40 keV–100 MeV on timescales of ≲10 s with greater sensitivity and imaging capability than its predecessors. This is achieved by deploying high-resolution scintillator detectors and indirect Fourier imaging techniques. LISSAN is suitable for inclusion in a multi-instrument platform such as an ESA M-class mission or as a smaller standalone mission. Without the observations that LISSAN can provide, our understanding of solar particle acceleration, and hence the space weather events with which it is often associated, cannot be complete.

Acceleration and Spectral Redistribution of Cosmic Rays in Radio-jet Shear Flows

A steady-state, semi-analytical model of energetic particle acceleration in radio-jet shear flows due to cosmic-ray viscosity obtained by Webb et al. is generalized to take into account more general cosmic-ray boundary spectra. This involves solving a mixed Dirichlet–Von Neumann boundary value problem at the edge of the jet. The energetic particle distribution function f 0(r, p) at cylindrical radius r from the jet axis (assumed to lie along the z-axis) is given by convolving the particle momentum spectrum with the Green’s function , which describes the monoenergetic spectrum solution in which as r → ∞ . Previous work by Webb et al. studied only the Green’s function solution for . In this paper, we explore for the first time, solutions for more general and realistic forms for . The flow velocity u = u(r) e z is along the axis of the jet (the z-axis). u is independent of z, and u(r) is a monotonic decreasing function of r. The scattering time in the shear flow region 0 < r < r 2, and , where s > 0 in the region r > r 2 is outside the jet. Other original aspects of the analysis are (i) the use of cosmic ray flow lines in (r, p) space to clarify the particle spatial transport and momentum changes and (ii) the determination of the probability distribution that particles observed at (r, p) originated from r → ∞ with momentum . The acceleration of ultrahigh-energy cosmic rays in active galactic nuclei jet sources is discussed. Leaky box models for electron acceleration are described.

Self-consistent modeling of the energetic storm particle event of November 10, 2012

Context. It is thought that solar energetic ions associated with coronal and interplanetary shock waves are accelerated to high energies by the diffusive shock acceleration mechanism. For this mechanism to be efficient, intense magnetic turbulence is needed in the vicinity of the shock. The enhanced turbulence upstream of the shock can be produced self-consistently by the accelerated particles themselves via streaming instability. Comparisons of quasi-linear-theory-based particle acceleration models that include this process with observations have not been fully successful so far, which has motivated the development of acceleration models of a different nature. Aims. Our aim is to test how well our self-consistent quasi-linear SOLar Particle Acceleration in Coronal Shocks (SOLPACS) simulation code, developed earlier to simulate proton acceleration in coronal shocks, models the particle foreshock region. Methods. We applied SOLPACS to model the energetic storm particle (ESP) event observed by the STEREO A spacecraft on November 10, 2012. Results. All but one main input parameter of SOLPACS are fixed by the in situ plasma measurements from the spacecraft. By comparing a simulated proton energy spectrum at the shock with the observed one, we were able to fix the last simulation input parameter related to the efficiency of particle injection to the acceleration process. A subsequent comparison of simulated proton time-intensity profiles in a number of energy channels with the observed ones shows a very good correspondence throughout the upstream region. Conclusions. Our results strongly support the quasi-linear description of the foreshock region.

The Maximum Energy of Shock-accelerated Cosmic Rays

Identifying the accelerators of Galactic cosmic ray (CR) protons with energies up to a few PeV (1015 eV) remains a theoretical and observational challenge. Supernova remnants (SNRs) represent strong candidates because they provide sufficient energetics to reproduce the CR flux observed at Earth. However, it remains unclear whether they can accelerate particles to PeV energies, particularly after the very early stages of their evolution. This uncertainty has prompted searches for other source classes and necessitates comprehensive theoretical modeling of the maximum proton energy, , accelerated by an arbitrary shock. While analytic estimates of have been put forward in the literature, they do not fully account for the complex interplay between particle acceleration, magnetic field amplification, and shock evolution. This paper uses a multizone, semianalytic model of particle acceleration based on kinetic simulations to place constraints on for a wide range of astrophysical shocks. In particular, we develop relationships between , shock velocity, size, and ambient medium. We find that SNRs can only accelerate PeV particles under a select set of circumstances, namely, if the shock velocity exceeds ∼104 km s−1 and escaping particles drive magnetic field amplification. However, older and slower SNRs may still produce observational signatures of PeV particles due to populations accelerated when the shock was younger. Our results serve as a reference for modelers seeking to quickly produce a self-consistent estimate of the maximum energy accelerated by an arbitrary astrophysical shock. 1 1 Presented as a thesis to the Department of Astronomy and Astrophysics, The University of Chicago, in partial fulfillment of the requirements for a Ph.D. degree.

Proton Synchrotron Origin of the Very-high-energy Emission of GRB 190114C

We consider here a proton-synchrotron model to explain the MAGIC observation of GRB 190114C afterglow in the energy band of 0.2–1 TeV, while the X-ray spectra are explained by electron-synchrotron emission. Given the uncertainty of the particle acceleration process, we consider several variations of the model, and show that they all match the data very well. We find that the values of the uncertain model parameters are reasonable: isotropic explosion energy ∼1054.5 erg, ambient density ∼10–100 cm−3, and the fraction of electrons/protons accelerated to a high-energy power law is of a few percent. All these values are directly derived from the observed teraelectronvolt and X-ray fluxes. The model also requires that protons be accelerated to observed energies as high as a few 1020 eV. Further, assuming that the jet break takes place after 106 s gives the beaming-corrected energy of the burst to be ≈1053 erg, which is one to two orders of magnitude higher than usually inferred, because of the small fraction of electrons accelerated. Our modeling is consistent with both late time data at all bands, from optical to X-rays, and with numerical models of particle acceleration. Our results thus demonstrate the relevance of proton-synchrotron emission to the high-energy observations of gamma-ray bursts during their afterglow phase.

Unexpected energetic particle observations near the Sun by Parker Solar Probe and Solar Orbiter

Solar energetic particles (SEPs) from suprathermal (few keV) up to relativistic (∼few GeV) energies are accelerated at the Sun in association with solar flares and coronal mass ejection-driven shock waves. Although our knowledge of the origin, acceleration, and transport of these particles from close to the Sun through the interplanetary medium has advanced dramatically in the last 40 years, many puzzles have still remained unsolved due to the scarcity of in situ measurements well inside 1 AU. Furthermore, energetic particle intensity enhancements associated with high-speed streams or stream interaction regions (SIRs) have been routinely observed at interplanetary spacecraft near Earth orbit since the 1960s. Since only a small sample of SIR events were observed by the Helios spacecraft inside 1 AU, additional observations well inside 1 AU were also needed to further investigate the energization and transport effects of SIR-associated ions and to compare with expectations from contemporary SIR-associated particle acceleration and transport models and theories. The Solar Orbiter (SolO) and Parker Solar Probe (PSP) pioneering missions have been providing unprecedented measurements of energetic particles in the near-Sun environment. This review presents the unexpected observations of SEP and SIR-related ion events as measured by the PSP/IS⊙IS and SolO/EPD experiments, which revealed surprises that challenge our understanding.

Quiet-time suprathermal ions in the inner heliosphere during the rising phase of solar cycle 25

Context. The Solar Orbiter spacecraft made its first close perihelion passes in 2022, reaching 0.32 au on 26 March and 0.29 au on 12 October. Transient activity was relatively low, making it possible to perform measurements of the quiet-time suprathermal ion pool over multi-day periods. Aims. The inner heliosphere suprathermal ion pool is a source of seed particles accelerated by coronal mass ejection-driven shocks. Determining its constituents and their dependence on solar activity, location, and time is therefore critical to building physical models of particle acceleration from solar events. Methods. By selecting low-activity periods on Solar Orbiter during perihelia passes, and comparing them with a nearly identical monitoring instrument at 1 au, the observed differences in intensities can be related to factors such as distance from the Sun, transient events, and interacting solar wind streams. Results. Below ∼1 MeV/nucleon, the observed quiet-time spectra at 0.32 au for H, 4He, 3He, and Fe rise toward low energies, as observed previously during longer periods at 1 au, and they show a heavy ion composition that has markers of impulsive solar flare material, such as relatively high 3He:4He, and a high Fe/O ratio. The proton and helium abundances are much higher, consistent with a source in corotating interaction regions. Surveying all semi-quiet times during the mission, there is only a modest (∼15%) increase in the fluences in the inner heliosphere compared to 1 au, indicating small gradients in these populations between 1 and 0.3 au.

Transfer learning and multi-fidelity modeling of laser-driven particle acceleration

Computer models of intense, laser-driven ion acceleration require expensive particle-in-cell simulations that may struggle to capture all the multi-scale, multi-dimensional physics involved at reasonable costs. Explored is an approach to ameliorate this deficiency using a multi-fidelity framework that can incorporate physical trends and phenomena at different levels. As the basis for this study, an ensemble of approximately 8000 1D simulations was generated to buttress separate ensembles of hundreds of higher fidelity 1D and 2D simulations. Using transfer learning with deep neural networks, one can reproduce the results of more complex physics at a much lower cost. The networks trained in this fashion can, in turn, act as surrogate models for the simulations themselves, allowing for quick and efficient exploration of the parameter space of interest. Standard figures-of-merit were used as benchmarks such as the hot electron temperature, peak ion energy, conversion efficiency, and so on. We can rapidly identify and explore under what conditions differing fidelities become an important effect and search for outliers in feature space.

Multiple emission components in the Cygnus cocoon detected fromFermi-LAT observations

Context. Star-forming regions may play an important role in the life cycle of Galactic cosmic rays (CRs), notably as home to specific acceleration mechanisms and transport conditions. Gamma-ray observations of Cygnus X have revealed the presence of an excess of hard-spectrum gamma-ray emission, possibly related to a cocoon of freshly accelerated particles.Aims. We seek an improved description of the gamma-ray emission from the cocoon using ~13 yr of observations with theFermi-Large Area Telescope (LAT) and use it to further constrain the processes and objects responsible for the young CR population.Methods. We developed an emission model for a large region of interest, including a description of interstellar emission from the background population of CRs and recent models for other gamma-ray sources in the field. Thus, we performed an improved spectro-morphological characterisation of the residual emission including the cocoon.Results. The best-fit model for the cocoon includes two main emission components: an extended component FCES G78.74+1.56, described by a 2D Gaussian of extensionr68= 4.4° ± 0.1°−0.1°+0.1°and a smooth broken power law spectrum with spectral indices 1.67 ± 0.05−0.01+0.02and 2.12 ± 0.02−0.01+0.00below and above 3.0 ± 0.6−0.2+0.0GeV, respectively; and a central component FCES G80.00+0.50, traced by the distribution of ionised gas within the borders of the photo-dissociation regions and with a power law spectrum of index 2.19 ± 0.03−0.01+0.00that is significantly different from the spectrum of FCES G78.74+1.56. An additional extended emission component FCES G78.83+3.57, located on the edge of the central cavities in Cygnus X and with a spectrum compatible with that of FCES G80.00+0.50, is likely related to the cocoon. For the two brightest components FCES G80.00+0.50 and FCES G78.74+1.56, spectra and radial-azimuthal profiles of the emission can be accounted for in a diffusion-loss framework involving one single population of non-thermal particles with a flat injection spectrum. Particles span the full extent of FCES G78.74+1.56 as a result of diffusion from a central source, and give rise to source FCES G80.00+0.50 by interacting with ionised gas in the innermost region.Conclusions. For this simple diffusion-loss model, viable setups can be very different in terms of energetics, transport conditions, and timescales involved, and both hadronic and leptonic scenarios are possible. The solutions range from long-lasting particle acceleration, possibly in prominent star clusters such as Cyg OB2 and NGC 6910, to a more recent and short-lived release of particles within the last 10–100 kyr, likely from a supernova remnant. The observables extracted from our analysis can be used to perform detailed comparisons with advanced models of particle acceleration and transport in star-forming regions.

Flaring activity from magnetic reconnection in BL Lacertae

ABSTRACT The evolution of the spectral energy distribution during flares constrains models of particle acceleration in blazar jets. The archetypical blazar BL Lacertae provided a unique opportunity to study spectral variations during an extended strong flaring episode from 2020 to 2021. During its brightest γ-ray state, the observed flux (0.1–300 GeV) reached up to $2.15\, \times \, 10^{-5}\, \rm {ph\, cm^{-2}\, s^{-1}}$, with sub-hour-scale variability. The synchrotron hump extended into the X-ray regime showing a minute-scale flare with an associated peak shift of inverse-Compton hump in γ-rays. In shock acceleration models, a high Doppler factor value &amp;gt;100 is required to explain the observed rapid variability, change of state, and γ-ray peak shift. Assuming particle acceleration in minijets produced by magnetic reconnection during flares, on the other hand, alleviates the constraint on required bulk Doppler factor. In such jet-in-jet models, observed spectral shift to higher energies (towards TeV regime) and simultaneous rapid variability arises from the accidental alignment of a magnetic plasmoid with the direction of the line of sight. We infer a magnetic field of ∼0.6 G in a reconnection region located at the edge of broad-line region (∼0.02 pc). The scenario is further supported by lognormal flux distribution arising from merging of plasmoids in reconnection region.

A low-frequency sub-arcsecond view of powerful radio galaxies in rich-cluster environments: 3C 34 and 3C 320

ABSTRACT Models of radio galaxy physics have been primarily based on high frequency (≥1 GHz) observations of their jets, hotspots, and lobes. Without highly resolved low-frequency observations, which provide information on older plasma, our understanding of the dynamics of radio galaxies and their interaction with their environment is limited. Here, we present the first sub-arcsecond (0.3 arcsec) resolution images at 144 MHz of two powerful radio galaxies situated in rich cluster environments, namely 3C 34 and 3C 320, using the International Low Frequency Array Telescope. We detect for the first time at low frequencies a plethora of structures in these objects, including strikingly large filaments across the base of the lobes in both sources, which are spatially associated with dense regions in the ambient medium. For 3C 34, we report a spectral flattening in the region of the central filament, suggesting that the origin of the filaments is related to the presence of large-scale ordered magnetic fields. We also report periodic total intensity and spectral index banding of diffuse emission in the eastern lobe, seen for the first time in radio galaxy lobes. The hotspot complexes are resolved into multiple fragments of varying structure and spectral index; we discuss the implications for particle acceleration and jet termination models. We find at most smooth gradients in the spectral behaviour of the hotspot structure suggesting that particle acceleration, if present, may be occurring throughout the complex, in contrast to simple models, but different jet termination models may apply to both sources.

First-Principles Fermi Acceleration in Magnetized Turbulence.

This Letter provides a concrete implementation of Fermi's model of particle acceleration in magnetohydrodynamic (MHD) turbulence, connecting the rate of energization to the gradients of the velocity of magnetic field lines, which it characterizes within a multifractal picture of turbulence intermittency. It then derives a transport equation in momentum space for the distribution function. This description is shown to be substantiated by a large-scale numerical simulation of strong MHD turbulence. The present general framework can be used to model particle acceleration in a variety of environments.

Operational inner magnetosphere particle transport and acceleration model (IMPTAM) for 1–300 keV electrons

Surface charging, the process of charge deposition on covering insulating surfaces of satellites is directly linked to the space environment at a time scale of a few tens of seconds. Accurate specification of the space environment at different orbits is of key importance. We introduce the model operating online in near-real time for low energy (< 200 keV) electrons in the inner magnetosphere, called Inner Magnetosphere Particle Transport and Acceleration model (IMPTAM). This model in its various versions has been operating online since March 2013 (latest version at imptam.engin.umich.edu) and it is driven by the real time solar wind (solar wind number density, dynamic pressure and velocity) and Interplanetary Magnetic Field (IMF), Y and Z components and total magnitude, parameters and by the real time Dst and Kp indices. The model provides the low energy electron (and proton) flux at all Magnetic Local Times (MLTs) and at L-shells from 3 to 10 and at all satellite orbits inside the modeling region, when necessary. IMPTAM output can serve as an input to assess surface charging levels of spacecraft immersed in severe environments and inside radiation belts models as distributions of electron seed population for further acceleration to MeV energies.