Search

The heating of the solar corona and solar wind, through suprathermal particles and kinetic Alfv\'en waves within the 0 - 10 $R_{\rm Sun}$ range, has been a subject of great interest for many decades. This study investigates the acceleration and heating of charged particles and the role of KAWs in the solar corona. We investigate how KAWs transport energy and accelerate/heat the charged particles, focusing on the behavior of perturbed EM fields, Poynting flux vectors, net power transfer, resonant particle speed, group speed, and the damping length of KAWs. The study examines how these elements are influenced by suprathermal particles \kappa and the electron-to-ion temperature $T_e/T_i$. We use kinetic plasma theory coupled with the Vlasov-Maxwell model to investigate the dynamics of KAWs and particles. We assume a collisionless, homogeneous, and low-beta electron-ion plasma in which Alfv\'en waves travel in the kinetic limits. The results show the perturbed EM fields are significantly influenced by $\kappa$ and $T_e/T_i$. We evaluate both the parallel and perpendicular Poynting fluxes and find that the parallel Poynting flux dissipates gradually for lower \kappa values. The perpendicular flux dissipates quickly over shorter distances. Power deposition in solar flux tubes is significantly influenced by \kappa and Te/Ti. We find that particles can heat the solar corona over long distances in the parallel direction and short distances in the perpendicular direction. The group velocity of KAWs increases for lower \kappa values, and the damping length is enhanced under lower \kappa, suggesting longer energy transport distances. These findings offer a comprehensive understanding of particle-wave interactions in the solar corona and wind, with potential applications for missions such as the Parker Solar Probe (PSP), and can also apply to other environments., Comment: Accepted for publication in the Astronomy and Astrophysics (A&A) journal (2024)

We present a pioneering model of the interaction between the solar wind and the surrounding interstellar medium that includes the possibility of different pressures in the directions parallel and perpendicular to the magnetic field. The outer heliosheath region is characterized by a low rate of turbulent scattering that would permit a development of pressure anisotropy. The effect is best seen on the interstellar side of the heliopause, where a narrow region is developed with an excessive perpendicular pressure resembling a plasma depletion layer typical of planetary magnetspheres. The magnitude of this effect is relatively small owing to the fact that proton-proton collisions don't permit the anisotropy to approach a stability threshold. We also show that if turbulence in the interstellar medium is weak, a much broader anisotropic boundary layer can exist, with a dominant perpendicular pressure in the southern hemisphere and a dominant parallel pressure in the north., Comment: 13 pages, 3 figures

Parker Solar Probe (PSP) observed sub-Alfvenic solar wind intervals during encounters 8 - 14, and low-frequency magnetohydrodynamic turbulence in these regions may differ from that in super-Alfvenic wind. We apply a new mode-decomposition analysis (Zank et al 2023) to the sub-Alfv\'enic flow observed by PSP on 2021 April 28, identifying and characterizing entropy, magnetic islands, forward and backward Alfv\'en waves, including weakly/non-propagating Alfv\'en vortices, forward and backward fast and slow magnetosonic modes. Density fluctuations are primarily and almost equally entropy and backward propagating slow magnetosonic modes. The mode-decomposition provides phase information (frequency and wavenumber k) for each mode. Entropy-density fluctuations have a wavenumber anisotropy k_{||} >> k_{perp} whereas slow mode density fluctuations have k_{perp} > k_{||}. Magnetic field fluctuations are primarily magnetic island modes (delta B^i) with an O(1) smaller contribution from uni-directionally propagating Alfven waves (delta B^{A+}) giving a variance anisotropy of <{\delta B^i}^2> / = 4.1. Incompressible magnetic fluctuations dominate compressible contributions from fast and slow magnetosonic modes. The magnetic island spectrum is Kolmogorov-like k_{perp}^{-1.6} in perpendicular wavenumber and the uni-directional Alfven wave spectra are k_{||}^{-1.6} and k_{perp}^{-1.5}. Fast magnetosonic modes propagate at essentially the Alfv\'en speed with anti-correlated transverse velocity and magnetic field fluctuations and are almost exclusively magnetic due to beta_p<<1. Transverse velocity fluctuations are the dominant velocity component in fast magnetosonic modes and longitudinal fluctuations dominate in slow modes. Mode-decomposition is an effective tool in identifying the basic building blocks of MHD turbulence and provides detailed phase information about each of the modes.

This paper investigates the self-pulsing of Dielectric Barrier Discharges (DBDs) at low driving frequencies. In particular, (a) the dependence of current on the product pd of gas pressure p and the gas gap length d, (b) the effects of lossy dielectrics (in resistive discharges) and large dielectric permittivity (in ferroelectrics) on current dynamics, (c) the transition from Townsend to a dynamic Capacitively Coupled Plasma (CCP) discharge with changing pd values, and (d) the transition from Townsend to a high-frequency CCP regime with increasing the driving frequency. A one-dimensional fluid model of Argon plasma is coupled to an equivalent RC circuit for lossy dielectrics. Our results show multiple current pulses per AC period in Townsend and CCP discharge modes which are explained by uncoupled electron-ion transport in the absence of quasineutrality and surface charge deposition at dielectric interfaces. The number of current pulses decreases with an increasing applied frequency when the Townsend discharge transforms into the CCP discharge. The resistive barrier discharge with lossy dielectrics exhibits Townsend and glow modes for the same pd value (7.6 Torr cm) for higher and lower resistances, respectively. Finally,we show that ferroelectric materials can amplify discharge current in DBDs. Similarities between current pulsing in DBD, Trichel pulses in corona discharges, and subnormal oscillations in DC discharges are discussed. 1

This Letter reports the first observational estimate of the heating rate in the slowly expanding solar corona. The analysis exploits the simultaneous remote and local observations of the same coronal plasma volume with the Solar Orbiter/Metis and the Parker Solar Probe instruments, respectively, and relies on the basic solar wind magnetohydrodynamic equations. As expected, energy losses are a minor fraction of the solar wind energy flux, since most of the energy dissipation that feeds the heating and acceleration of the coronal flow occurs much closer to the Sun than the heights probed in the present study, which range from 6.3 to 13.3 solar radii. The energy deposited to the supersonic wind is then used to explain the observed slight residual wind acceleration and to maintain the plasma in a non-adiabatic state. As derived in the Wentzel-Kramers-Brillouin limit, the present energy transfer rate estimates provide a lower limit, which can be very useful in refining the turbulence-based modeling of coronal heating and subsequent solar wind acceleration.

HelioSwarm (HS) is a NASA Medium-Class Explorer mission of the Heliophysics Division designed to explore the dynamic three-dimensional mechanisms controlling the physics of plasma turbulence, a ubiquitous process occurring in the heliosphere and in plasmas throughout the universe. This will be accomplished by making simultaneous measurements at nine spacecraft with separations spanning magnetohydrodynamic and sub-ion spatial scales in a variety of near-Earth plasmas. In this paper, we describe the scientific background for the HS investigation, the mission goals and objectives, the observatory reference trajectory and instrumentation implementation before the start of Phase B. Through multipoint, multiscale measurements, HS promises to reveal how energy is transferred across scales and boundaries in plasmas throughout the universe., Comment: Submitted to Space Science Reviews, 60 pages, 15 figures, 2 videos

The transport of energetic particles in response to solar wind turbulence is important for space weather. To understand charged particle transport, it is usually assumed that the phase of the turbulence is randomly distributed (the random phase approximation) in quasi-linear theory and simulations. In this paper, we calculate the coherence index, $C_\phi$, of solar wind turbulence observed by the Helios 2 and Parker Solar Probe spacecraft using the surrogate data technique to check if the assumption is valid. Here, values of $C_\phi=0$ and $1$ indicate that the phase coherence is random and correlated, respectively. We estimate that the coherence index at the resonant scale of energetic ions ($10$ MeV protons) is $0.1$ at $0.87$ and $0.65$ au, $0.18$ at $0.29$ au, and $0.3$ ($0.35$) at $0.09$ au for super (sub)-Alfv\'enic intervals, respectively. Since the random phase approximation corresponds to $C_\phi=0$, this may indicate that the random phase approximation is not valid for the transport of energetic particles in the inner heliosphere, especially very close to the Sun ($\sim0.09$ au)., Comment: 15 pages, 4 figures, published in Frontiers in Astronomy and Space Sciences

Evidence for the presence of ion cyclotron waves, driven by turbulence, at the boundaries of the current sheet is reported in this paper. By exploiting the full potential of the joint observations performed by Parker Solar Probe and the Metis coronagraph on board Solar Orbiter, local measurements of the solar wind can be linked with the large-scale structures of the solar corona. The results suggest that the dynamics of the current sheet layers generates turbulence, which in turn creates a sufficiently strong temperature anisotropy to make the solar-wind plasma unstable to anisotropy-driven instabilities such as the Alfv\'en ion-cyclotron, mirror-mode, and firehose instabilities. The study of the polarization state of high-frequency magnetic fluctuations reveals that ion cyclotron waves are indeed present along the current sheet, thus linking the magnetic topology of the remotely imaged coronal source regions with the wave bursts observed in situ. The present results may allow improvement of state-of-the-art models based on the ion cyclotron mechanism, providing new insights into the processes involved in coronal heating., Comment: 12 pages

We present a model for atmospheric absorption of solar ultraviolet (UV) radiation. The initial motivation for this work is to predict this effect and correct it in Sounding Rocket (SR) experiments. In particular, the Full-sun Ultraviolet Rocket Spectrograph (FURST) is anticipated to launch in mid-2023. FURST has the potential to observe UV absorption while imaging solar spectra between 120-181 nm, at a resolution of R > 2x10$^4$ ($\Delta$ V < $\pm$ 15 km/s), and at altitudes of between 110-255 km. This model uses estimates for density and temperature, as well as laboratory measurements of the absorption cross-section, to predict the UV absorption of solar radiation at high altitudes. Refraction correction is discussed and partially implemented but is negligible for the results presented. Absorption by molecular Oxygen is the primary driver within the UV spectral range of our interest. The model is built with a wide range of applications in mind. The primary result is a method for inversion of the absorption cross-section from images obtained during an instrument flight, even if atmospheric observations were not initially intended. The potential to obtain measurements of atmospheric properties is an exciting prospect, especially since sounding rockets are the only method currently available for probing this altitude in situ. Simulation of noisy spectral images along the FURST flight profile is performed using data from the High-Resolution Telescope and Spectrograph (HRTS) SR and the FISM2 model for comparison. Analysis of these simulated signals allows us to capture the Signal-to-Noise Ratio (SNR) of FURST and the capability to measure atmospheric absorption properties as a function of altitude. Based on the prevalence of distinct spectral features, our calculations demonstrate that atmospheric absorption may be used to perform wavelength calibration from in-flight data., Comment: To be Published in JPCS. Submitted December 2022. Accepted February 2023

The transport of waves and turbulence beyond the photosphere is central to the coronal heating problem. Turbulence in the quiet solar corona has been modeled on the basis of the nearly incompressible magnetohydrodynamic (NI MHD) theory to describe the transport of low-frequency turbulence in open magnetic field regions. It describes the evolution of the coupled majority quasi-2D and minority slab component, driven by the magnetic carpet and advected by a subsonic, sub-Alfvenic flow from the lower corona. In this paper, we couple the NI MHD turbulence transport model with an MHD model of the solar corona to study the heating problem in a coronal loop. In a realistic benchmark coronal loop problem, we find that a loop can be heated to ~1.5 million K by transport and dissipation of MHD turbulence described by the NI MHD model. We also find that the majority 2D component is as important as the minority slab component in the heating of the coronal loop. We compare our coupled MHD/NI MHD model results with a reduced MHD (RMHD) model. An important distinction between these models is that RMHD solves for small-scale velocity and magnetic field fluctuations and obtains the actual viscous/resistive dissipation associated with their evolution whereas NI MHD evolves scalar moments of the fluctuating velocity and magnetic fields and approximates dissipation using an MHD turbulence phenomenology. Despite the basic differences between the models, their simulation results match remarkably well, yielding almost identical heating rates inside the corona., Comment: 23 pages, 6 figures, 1 table, The Astrophysical Journal, in press

Owing to the ever-present solar wind, our vast solar system is full of plasmas. The turbulent solar wind, together with sporadic solar eruptions, introduces various space plasma processes and phenomena in the solar atmosphere all the way to the Earth's ionosphere and atmosphere and outward to interact with the interstellar media to form the heliopause and termination shock. Remarkable progress has been made in space plasma physics in the last 65 years, mainly due to sophisticated in-situ measurements of plasmas, plasma waves, neutral particles, energetic particles, and dust via space-borne satellite instrumentation. Additionally high technology ground-based instrumentation has led to new and greater knowledge of solar and auroral features. As a result, a new branch of space physics, i.e., space weather, has emerged since many of the space physics processes have a direct or indirect influence on humankind. After briefly reviewing the major space physics discoveries before rockets and satellites, we aim to review all our updated understanding on coronal holes, solar flares and coronal mass ejections, which are central to space weather events at Earth, solar wind, storms and substorms, magnetotail and substorms, emphasizing the role of the magnetotail in substorm dynamics, radiation belts/energetic magnetospheric particles, structures and space weather dynamics in the ionosphere, plasma waves, instabilities, and wave-particle interactions, long-period geomagnetic pulsations, auroras, geomagnetically induced currents (GICs), planetary magnetospheres and solar/stellar wind interactions with comets, moons and asteroids, interplanetary discontinuities, shocks and waves, interplanetary dust, space dusty plasmas and solar energetic particles and shocks, including the heliospheric termination shock. This paper is aimed to provide a panoramic view of space physics and space weather., Comment: Accepted for publication in IEEE Transactions on Plasma Science (2022)

On its journey through the Galaxy, the Sun passes through diverse regions of the interstellar medium. High-resolution spectroscopic measurements of interstellar absorption lines in spectra of nearby stars show absorption components from more than a dozen warm partially ionized clouds within 15 pc of the Sun. The two nearest clouds - the Local Interstellar Cloud (LIC) and Galactic (G) cloud - move toward each other. Their bulk heliocentric velocities can be compared with the interstellar neutral helium flow velocity obtained from space-based experiments. We combine recent results from Ulysses, IBEX, and STEREO observations to find a more accurate estimate of the velocity and temperature of the very local interstellar medium. We find that, contrary to the widespread viewpoint that the Sun resides inside the LIC, the locally observed velocity of the interstellar neutral helium is consistent with a linear combination of the velocities of the LIC and G cloud, but not with either of these two velocities. This finding shows that the Sun travels through a mixed-cloud interstellar medium composed of material from both these clouds. Interactions between these clouds explain the substantially higher density of the interstellar hydrogen near the Sun and toward stars located within the interaction region of these two clouds. The observed asymmetry of the interstellar helium distribution function also supports this interaction. The structure and equilibrium in this region require further studies using in situ and telescopic observations., Comment: 13 pages, 5 figures, 1 table, accepted for publication in ApJL

We present a first theoretical modeling of joint Parker Solar Probe (PSP) - Metis/Solar Orbiter (SolO) quadrature observations Telloni et al 2022c. The combined observations describe the evolution of a slow solar wind plasma parcel from the extended solar corona ($3.5-6.3$ R$_\odot$) to the very inner heliosphere (23.2 R$_\odot$). The Metis/SolO instrument remotely measures the solar wind speed finding a range from $96-201$ kms$^{-1}$, and PSP measures the solar wind plasma in situ, observing a radial speed of 219.34 kms$^{-1}$. We find theoretically and observationally that the solar wind speed accelerates rapidly within 3.3 -- 4 R$_\odot$, and then increases more gradually with distance. Similarly, we find that the theoretical solar wind density is consistent with the remotely and in situ observed solar wind density. The normalized cross-helicity and normalized residual energy observed by PSP are 0.96 and -0.07, respectively, indicating that the slow solar wind is very Alfv\'enic. The theoretical NI/slab results are very similar to PSP measurements, which is a consequence of the highly magnetic field-aligned radial flow ensuring that PSP can measure slab fluctuations and not 2D. Finally, we calculate the theoretical 2D and slab turbulence pressure, finding that the theoretical slab pressure is very similar to that observed by PSP., Comment: 12 pages, 5 figures

Switchbacks are sudden, large radial deflections of the solar wind magnetic field, widely revealed in interplanetary space by the Parker Solar Probe. The switchbacks' formation mechanism and sources are still unresolved, although candidate mechanisms include Alfv\'enic turbulence, shear-driven Kelvin-Helmholtz instabilities, interchange reconnection, and geometrical effects related to the Parker spiral. This Letter presents observations from the Metis coronagraph onboard Solar Orbiter of a single large propagating S-shaped vortex, interpreted as first evidence of a switchback in the solar corona. It originated above an active region with the related loop system bounded by open-field regions to the East and West. Observations, modeling, and theory provide strong arguments in favor of the interchange reconnection origin of switchbacks. Metis measurements suggest that the initiation of the switchback may also be an indicator of the origin of slow solar wind.

Zank et al. developed models describing the transport of low frequency incompressible and nearly incompressible turbulence in inhomogeneous flows. The formalism was based on expressing the fluctuating variables in terms of the Els\"assar variables and then taking "moments" subject to various closure hypotheses. The turbulence transport models are different according to whether the plasma beta regime is large or of order 1 or smaller. Here, we show explicitly that the two sets of turbulence transport models admit a conservation representation that resembles the well-known WKB transport equation for Alfv\'en wave energy density after introducing appropriate definitions of the "pressure" associated with the turbulent fluctuations. This includes introducing a distinct turbulent pressure tensor for 3D incompressible turbulence (the large plasma beta limit) and pressure tensors for quasi-2D and slab turbulence (the plasma beta order 1 or small regimes) that generalize the form of the WKB pressure tensor. Various limits of the different turbulent pressure tensors are discussed. However, the analogy between the conservation form of the turbulence transport models and the WKB model is not close for multiple reasons, including that the turbulence models express fully nonlinear physical processes unlike the strictly linear WKB description. The analysis presented here serves both as a check on the validity and correctness of the turbulence transport models and provides greater transparency of the energy dissipation term and the "turbulent pressure" in our models, which is important for many practical applications., Comment: 19 pages, 0 figure

This Letter addresses the first Solar Orbiter (SO) -- Parker Solar Probe (PSP) quadrature, occurring on January 18, 2021, to investigate the evolution of solar wind from the extended corona to the inner heliosphere. Assuming ballistic propagation, the same plasma volume observed remotely in corona at altitudes between 3.5 and 6.3 solar radii above the solar limb with the Metis coronagraph on SO can be tracked to PSP, orbiting at 0.1 au, thus allowing the local properties of the solar wind to be linked to the coronal source region from where it originated. Thanks to the close approach of PSP to the Sun and the simultaneous Metis observation of the solar corona, the flow-aligned magnetic field and the bulk kinetic energy flux density can be empirically inferred along the coronal current sheet with an unprecedented accuracy, allowing in particular estimation of the Alfv\'en radius at 8.7 solar radii during the time of this event. This is thus the very first study of the same solar wind plasma as it expands from the sub-Alfv\'enic solar corona to just above the Alfv\'en surface., Comment: 10 pages, 4 figures

The paper presents a new method for deriving turbulent heating of the solar wind using plasma moments and magnetic field data. We develop the method and then apply it to compute the turbulent heating of the solar wind proton plasma at 1AU. The method employs two physical properties of the expanding solar wind plasma, the wave-particle thermodynamic equilibrium, and the transport of entropic rate. We analyze plasma moments and field datasets taken from Wind S/C, in order to compute (i) the fluctuating magnetic energy, (ii) the corresponding correlation length, and (iii) the turbulent heating rate. We identify their relationships with the solar wind speed, as well as the variation of these relationships relative to solar wind and interplanetary coronal mass ejection plasma., Comment: 16 pages, 12 figures

The physics of the solar chromosphere is complex from both theoretical and modeling perspectives. The plasma temperature from the photosphere to corona increases from ~5,000 K to ~1 million K over a distance of only ~10,000 km from the chromosphere and the transition region. Certain regions of the solar atmosphere have sufficiently low temperature and ionization rates to be considered as weakly-ionized. In particular, this is true at the lower chromosphere. In this paper, we present an overview of our data-driven magnetohydrodynamics model for the weakly-ionized chromosphere and show a benchmark result. It utilizes the Cowling resistivity which is orders of magnitude greater than the Coulomb resistivity. Ohm's law therefore includes anisotropic dissipation. We investigate the effects of the Cowling resistivity on heating and magnetic reconnection in the chromosphere as the flare-producing active region (AR) 11166 evolves. In particular, we analyze a C2.0 flare emerging from AR11166 and find a normalized reconnection rate of 0.12., Comment: 14 pages, 4 figures, Journal of Physics Conference Series - Proceedings of the 19th Annual International Astrophysics Conference, Santa Fe, NM, March 9-13, 2020 (accepted)

The physics of the solar chromosphere is complex from both theoretical and modeling perspectives. The plasma temperature from the photosphere to corona increases from ~5,000 K to ~1 million K over a distance of only ~10,000 km from the chromosphere and the transition region. Certain regions of the solar atmosphere have sufficiently low temperature and ionization rates to be considered as weakly-ionized. In particular, this is true at the lower chromosphere. As a result, the Cowling resistivity is orders of magnitude greater than the Coulomb resistivity. Ohm's law therefore includes anisotropic dissipation. To evaluate the Cowling resistivity, we need to know the external magnetic field strength and to estimate the neutral fraction as a function of the bulk plasma density and temperature. In this study, we determine the magnetic field topology using the non-force-free field (NFFF) extrapolation technique based on SDO/HMI SHARP vector magnetogram data, and the stratified density and temperature profiles from the Maltby-M umbral core model for sunspots. We investigate the variation and effects of Cowling resistivity on heating and magnetic reconnection in the chromosphere as the flare-producing active region (AR) 11166 evolves. In particular, we analyze a C2.0 flare emerging from AR11166 and find a normalized reconnection rate of 0.051., Comment: 15 pages, 5 figures, The Astrophysical Journal Letters (in press)

An observation-based Galactic Cosmic Ray (GCR) spectral model for heavy nuclei is developed. Zhao and Qin (J. Geophys. Res. Space Phys.118, 1837(2013)) proposed an empirical elemental GCR spectra model for nuclear charge 5-28 over the energy range from 30 to 500 MeV/nuc, which is proved to be successful in predicting yearly averaged GCR heavy nuclei spectra.Based on the latest highly statistically precise measurements from ACE/CRIS,a further elemental GCR model with monthly averaged spectra is presented. The model can reproduce the past and predict the futureGCR intensity monthly by correlating model parameters with thecontinuous sunspot number (SSN) record. The effects of solar activity on GCR modulation are considered separately for odd and even solar cycles. Compared with other comprehensive GCR models, our modeling results are satisfyingly consistent with the GCR spectral measurements from ACE/SIS and IMP-8, and have comparable prediction accuracy as the Badhwar & O'Neill 2014 model.A detailed error analysis is also provided.Finally, the GCR carbon and iron nuclei fluxes for the subsequent two solar cycles (SC 25 and 26) are predicted and they show a potential trend in reduced flux amplitude, which is suspected to be relevant to possible weak solar cycles.

Parker Solar Probe (PSP) achieved its first orbit perihelion on November 6, 2018, reaching a heliocentric distance of about 0.165 au (35.55 R$_\odot$). Here, we study the evolution of fully developed turbulence associated with the slow solar wind along the PSP trajectory between 35.55 R$_\odot$ and 131.64 R$_\odot$ in the outbound direction, comparing observations to a theoretical turbulence transport model. Several turbulent quantities, such as the fluctuating kinetic energy and the corresponding correlation length, the variance of density fluctuations, and the solar wind proton temperature are determined from the PSP SWEAP plasma data along its trajectory between 35.55 R$_\odot$ and 131.64 R$_\odot$. The evolution of the PSP derived turbulent quantities are compared to the numerical solutions of the nearly incompressible magnetohydrodynamic (NI MHD) turbulence transport model recently developed by Zank et al. (2017). We find reasonable agreement between the theoretical and observed results. On the basis of these comparisons, we derive other theoretical turbulent quantities, such as the energy in forward and backward propagating modes, the total turbulent energy, the normalized residual energy and cross-helicity, the fluctuating magnetic energy, and the correlation lengths corresponding to forward and backward propagating modes, the residual energy, and the fluctuating magnetic energy.

It is well known that collisionless systems are dissipation free from the perspective of particle collision and thus conserve entropy. On the other hand, processes such as magnetic reconnection and turbulence appear to convert large-scale magnetic energy into heat. In this paper, we investigate the energization and heating of collisionless plasma. The dissipation process is discussed in terms of fluid entropy in both isotropic and gyrotropic forms. Evolution equations for the entropy are derived and they reveal mechanisms that lead to changes in fluid entropy. These equations are verified by a collisionless particle-in-cell simulation of multiple reconnecting current sheets. In addition to previous findings regarding the pressure tensor, we emphasize the role of heat flux in the dissipation process., Comment: 19 pages, 4 figure

Recent detection of superflares on solar-type stars by Kepler mission raised a possibility that they can be associated with energetic coronal mass ejections (CMEs) and energetic particle events (SEPs). These space weather events can impact habitability of exoplanets around these stars. Here we use the improved Particle Acceleration and Transport in the Heliosphere (iPATH) model, to model the time intensity profile and spectrum of SEPs accelerated at CME-driven shocks from stars of different ages traced by their rotation rates. We consider a solar-like (G-type) star with 6 different rotation rates varying from 0.5 to 3.0 times rotation speed of the sun. In all 6 cases, a fast CME is launched with the same speed of 1500 km/sec and the resulting time intensity profiles at 3 locations and and energy spectra at 5 locations at 1 AU are obtained. The maximum particle energy at the shock front as a function of r is also shown. Our results suggest that within 0.8 AU the maximum particle energy at the shock front increases with the rotation rate of the star. However, event integrated spectra for the five selected locations along the CME path show complicated patterns. This is because the Parker magnetic field for rapidly rotating stars is more tightly winded. Our results can be used in estimating the radiation environments of terrestrial-type exoplanets around solar-type stars., Comment: Manuscript has been accepted by The Astrophysical Journal Letters

Magnetic flux ropes are commonly observed throughout the heliosphere, and recent studies suggest that interacting flux ropes are associated with some energetic particle events. In this work, we carry out 2D particle-in-cell (PIC) simulations to study the coalescence of two magnetic flux ropes (or magnetic islands), and the subsequent plasma energization processes. The simulations are initialized with two magnetic islands embedded in a reconnecting current sheet. The two islands collide and eventually merge into a single island. Particles are accelerated during this process as the magnetic energy is released and converted to the plasma energy, including bulk kinetic energy increase by the ideal electric field, and thermal energy increase by the fluid compression and the non-ideal electric field. We find that contributions from these different energization mechanisms are all important and comparable with each other. Fluid shear and a non-gyrotropic pressure tensor also contribute to the energy conversion process. For simulations with different box sizes ranging from $L_x \sim 25$--$100 d_i$ and ion-to-electron mass ratios $m_i / m_e = 25$, 100 and 400, we find that the general evolution is qualitatively the same for all runs, and the energization depends only weakly on either the system size or the mass ratio. The results may help us understand plasma energization in solar and heliospheric environments.

Pickup ions are created when interstellar neutral atoms resonantly exchange charge with the solar wind (SW) ions, especially in the supersonic part of the wind, where they carry most of the plasma pressure. Here we present numerical simulation results of the 3D heliospheric interface treating pickup ions as a separate proton fluid. To satisfy the fundamental conservation laws, we solve the system of equations describing the flow of the mixture of electrons, thermal protons, and pickup ions. To find the density and pressure of pickup ions behind the termination shock, we employ simple boundary conditions that take into account the \emph{Voyager} observations that showed that the decrease in the kinetic energy of the mixture at the termination shock predominantly contributed to the increase in the pressure of pickup ions. We show that this model adequately describes the flow of the plasma mixture and results in a noticeable decrease in the heliosheath width., Comment: Accepted for publication in Journal of Physics, Conference Series

In this paper we discuss conservation laws in ideal magnetohydrodynamics (MHD) and gas dynamics associated with advected invariants. The invariants in some cases, can be related to fluid relabelling symmetries associated with the Lagrangian map. There are different classes of invariants that are advected or Lie dragged with the flow. Simple examples are the advection of the entropy S (a 0-form), and the conservation of magnetic flux (an invariant 2-form advected with the flow). The magnetic flux conservation law is equivalent to Faraday's equation. We discuss the gauge condition required for the magnetic helicity to be advected with the flow. The conditions for the cross helicity to be an invariant are discussed. We discuss the different variants of helicity in fluid dynamics and in MHD, including: fluid helicity, cross helicity and magnetic helicity. The fluid helicity conservation law and the cross helcity conservation law in MHD are derived for the case of a barotropic gas. If the magnetic field lies in the constant entropy surface, then the gas pressure can depend on both the entropy and the density. In these cases the conservation laws are local conservation laws. We obtain nonlocal conservation laws for fluid helicity and cross helicity for non-barotropic fluids using the Clebsch variable formulation of gas dynamics and MHD. Ertel's theorem and potential vorticity, the Hollman invariant, and the Godbillon Vey invariant for special flows for which the magnetic helicity is zero are also discussed., Comment: 40 pages, 0 figures

It has been noted for a long time that the spectra of observed continuum emissions in many solar flares are consistent with double power laws with a hardening at energies $\sim > $ 300 keV. It is now largely believed that at least in electron-dominated events the hardening in photon spectrum reflects an intrinsic hardening in the source electron spectrum. In this paper, we point out that a power law spectrum of electron with a hardening at high energies can be explained by diffusive shock acceleration of electrons at a termination shock with a finite width. Our suggestion is based on an early analytical work by Drury et al., where the steady state transport equation at a shock with a tanh profile was solved for a $p$-independent diffusion coefficient. Numerical simulations with a $p$-dependent diffusion coefficient show hardenings in the accelerated electron spectrum which are comparable with observations. One necessary condition for our proposed scenario to work is that high energy electrons resonate with the inertial range of the MHD turbulence and low energy electrons resonate with the dissipation range of the MHD turbulence at the acceleration site, and the spectrum of the dissipation range $\sim k^{-2.7}$. A $\sim k^{-2.7}$ dissipation range spectrum is consistent with recent solar wind observations., Comment: 20 pages, 2 figures, accepted by ApJ

The solar wind carves a bubble in the surrounding interstellar medium (ISM), known as the heliosphere. Charged interstellar dust grains (ISDG) encountering the heliosphere may be diverted around the heliopause or penetrate it depending on their charge-to-mass ratio. We present new calculations of trajectories of ISDG in the heliosphere, and the dust density distributions that result. We include up-to-date grain charging calculations using a realistic UV radiation field and full 3-D magnetohydrodynamic fluid + kinetic models for the heliosphere. Models with two different (constant) polarities for the solar wind magnetic field (SWMF) are used, with the grain trajectory calculations done separately for each polarity. Small grains a_gr ~ 0.01 micron are completely excluded from the inner heliosphere. Large grains, a_gr ~ 1.0 micron pass into the inner solar system and are concentrated near the Sun by its gravity. Trajectories of intermediate size grains depend strongly on the SWMF polarity. When the field has magnetic north pointing to ecliptic north, the field de-focuses the grains resulting in low densities in the inner heliosphere, while for the opposite polarity the dust is focused near the Sun. The ISDG density outside the heliosphere inferred from applying the model results to in situ dust measurements is inconsistent with local ISM depletion data for both SWMF polarities, but is bracketed by them. This result points to the need to include the time variation in the SWMF polarity during grain propagation. Our results provide valuable insights for interpretation of the in situ dust observations from Ulysses., Comment: 34 pages, 15 figures, Accepted for publication in the Astrophysical Journal

A local turbulence model is developed to study energy cascades in the interstellar medium (ISM) based on self-consistent two-dimensional fluid simulations. The model describes a partially ionized magnetofluid interstellar medium (ISM) that couples a neutral hydrogen fluid with a plasma primarily through charge exchange interactions. Charge exchange interactions are ubiquitous in warm ISM plasma and the strength of the interaction depends largely on the relative speed between the plasma and the neutral fluid. Unlike small length-scale linear collisional dissipation in a single fluid, charge exchange processes introduce channels that can be effective on a variety of length-scales that depend on the neutral and plasma densities, temperature, relative velocities, charge exchange cross section and the characteristic length scales. We find, from scaling arguments and nonlinear coupled fluid simulations, that charge exchange interactions modify spectral transfer associated with large scale energy containing eddies. Consequently, the warm ISM turbulent cascade rate prolongs spectral transfer amongst inertial range turbulent modes. Turbulent spectra associated with the neutral and plasma ISM fluids are therefore steeper than those predicted by Kolmogorov's phenomenology., Comment: Paper is accepted in Astrophysical Journal

New solar wind data from the Voyager 1 and Voyager 2 spacecraft, together with the SOHO SWAN measurements of the direction that neutral hydrogen enters into the inner heliosheath and neutral helium measurements provided by multiple observations are expected to provide more reliable constraints on the ionization ratio of the local interstellar medium (LISM) and the direction and magnitude of the interstellar magnetic field (ISMF). In this paper we use currently the most sophisticated numerical model of the heliospheric interface, which is based on an MHD treatment of the ion flow and kinetic modeling of neutral particles, to analyze an ISMF-induced asymmetry of the heliosphere in the presence of the interplanetary magnetic field and neutral particles. It is shown that secondary hydrogen atoms modify the LISM properties leading to its shock-free deceleration at the heliopause. We determine the deflection of hydrogen atoms from their original trajectory in the unperturbed LISM and show that it occurs not only in the plane defined by the ISMF and LISM velocity vectors, but also, to a lesser extent, perpendicular to this plane. We also consider the possibility of using 2-3 kHz radio emission data to further constrain the ISMF direction., Comment: Paper accepted for publication in the Astrophysical Journal Letters

At present, the heliosphere is embedded in a warm low density interstellar cloud that belongs to a cloud system flowing through the local standard of rest with a velocity near ~18 km/s. The velocity structure of the nearest interstellar material (ISM), combined with theoretical models of the local interstellar cloud (LIC), suggest that the Sun passes through cloudlets on timescales of < 10^3 - 10^4 yr, so the heliosphere has been, and will be, exposed to different interstellar environments over time. By means of a multi-fluid model that treats plasma and neutral hydrogen self-consistently, the interaction of the solar wind with a variety of partially ionized ISM is investigated, with the focus on low density cloudlets such as are currently near the Sun. Under the assumption that the basic solar wind parameters remain/were as they are today, a range of ISM parameters (from cold neutral to hot ionized, with various densities and velocities) is considered. In response to different interstellar boundary conditions, the heliospheric size and structure change, as does the abundance of interstellar and secondary neutrals in the inner heliosphere, and the cosmic ray level in the vicinity of Earth. Some empirical relations between interstellar parameters and heliospheric boundary locations, as well as neutral densities, are extracted from the models., Comment: 24 pages, 9 figures, 2 tables