Coronal Mass Ejections (CMEs) are one of the primary drivers of extreme space weather. They are large eruptions of mass and magnetic field from the solar corona and can travel the distance between Sun and Earth in half a day to a few days. Predictions of CMEs at 1 Astronomical Unit (AU), in terms of both its arrival time and magnetic field configuration, are very important for predicting space weather. Magnetohydrodynamic (MHD) modeling of CMEs, using flux-rope-based models is a promising tool for achieving this goal. In this study, we present one such model for CME simulations, based on spheromak magnetic field configuration. We have modified the spheromak solution to allow for independent input of poloidal and toroidal fluxes. The motivation for this is a possibility to estimate these fluxes from solar magnetograms and extreme ultraviolet (EUV) data from a number of different approaches. We estimate the poloidal flux of CME using post eruption arcades (PEAs) and toroidal flux from the coronal dimming. In this modified spheromak, we also have an option to control the helicity sign of flux ropes, which can be derived from the solar disk magnetograms using the magnetic tongue approach. We demonstate the applicability of this model by simulating the 12 July 2012 CME in the solar corona.

The magnetic fields of interplanetary coronal mass ejections (ICMEs), which originate close to the Sun in the form of a flux rope, determine their geoeffectiveness. Therefore, robust flux rope-based models of CMEs are required to perform magnetohydrodynamic (MHD) simulations aimed at space weather predictions. We propose a modified spheromak model and demonstrate its applicability to CME simulations. In this model, such properties of a simulated CME as the poloidal and toroidal magnetic fluxes, and the helicity sign can be controlled with a set of input parameters. We propose a robust technique for introducing CMEs with an appropriate speed into a background, MHD solution describing the solar wind in the inner heliosphere. Through a parametric study, we find that the speed of a CME is much more dependent on its poloidal flux than on the toroidal flux. We also show that the CME speed increases with its total energy, giving us control over its initial speed. We further demonstrate the applicability of this model to simulations of CME-CME collisions. Finally, we use this model to simulate the 12 July 2012 CME and compare the plasma properties at 1 AU with observations. The predicted CME properties agree reasonably with observational data.

Since the launch on 2018/08/12, Parker Solar Probe (PSP) has completed its first and second orbits around the Sun, having reached down to 35.7 solar radii at each perihelion. In anticipation of the exciting new data at such unprecedented distances, we have simulated the global 3D heliosphere using an MHD model coupled with a semi-empirical coronal model using the best available photospheric magnetograms as input. We compare our heliospheric MHD simulation results with in situ measurements along the PSP trajectory from its launch to the completion of the second orbit, with particular emphasis on the solar wind structure around the first two solar encounters. Furthermore, we show our model prediction for the third perihelion, which occurred on 2019/09/01. Comparison of the MHD results with PSP observations provides a new insight on the solar wind acceleration. Moreover, PSP observations reveal how accurately the ADAPT-WSA predictions work throughout the inner heliosphere.

Energetic neutral atoms (ENA) are an important tool for investigating the structure of the heliosphere. Recently, it was observed that fluxes of ENAs (with energy $\le$ 55 keV) coming from the upwind and downwind regions of the heliosphere are similar in strength. This led the authors of these observations to hypothesize that the heliosphere is bubble-like rather than comet-like, meaning that it has no extended tail. We investigate the directional distribution of the ENA flux for a wide energy range (3--88 keV) including the observations from IBEX (Interstellar Boundary Explorer), INCA (Ion and Neutral Camera, on board Cassini), and HSTOF (High energy Suprathermal Time Of Flight sensor, on board SOHO, Solar and Heliospheric Observatory). An essential element is the model of pickup ion acceleration at the termination shock (TS) proposed by Zank. We use state of the art models of the global heliosphere, interstellar neutral gas density, and pickup ion distributions. The results, based on the "comet-like" model of the heliosphere, are close in flux magnitude to ENA observations by IBEX, HSTOF and partly by INCA (except for the 5.2-13.5 keV energy channel). We find that the ENA flux from the tail dominates at high energy (in agreement with HSTOF, but not INCA). At low energy, our comet-like model produces the similar strengths of the ENA fluxes from the upwind and downwind directions, which, therefore, removes this as a compelling argument for a bubble-like heliosphere.

It is currently believed that the turbulent fluctuations pervade the outermost heliosphere. Turbulence, magnetic reconnection, and their link may be responsible for magnetic energy conversion in these regions. The governing mechanisms of such anisotropic and compressible magnetic turbulence in the inner heliosheath (IHS) and in the local interstellar medium (LISM) still lack a thorough description. The present literature mainly concerns large scales which are not representative of the inertial-cascade dynamics of turbulence. Moreover, lack of broadband spectral analysis makes the IHS dynamics critically understudied. Our recent study shows that 48 s magnetic-field data from the Voyager mission are appropriate for a spectral analysis over a frequency range of six decades, from 5 x 10-8 Hz to 10-2 Hz. Here, focusing on the Voyager 2 observation interval from 2013.824 to 2016.0, we describe the structure of turbulence in a sector zone of the IHS. A spectral break around 7 x 10-7 Hz (magnetic structures with size l~1.3 Astronomical Units) separates the energy-injection regime from the inertial-cascade regime of turbulence. A second scale is observed around 6 x 10-5 Hz (l~ 0.017 AU) and corresponds to a peak of compressibility and intermittency of fluctuations.

The understanding of inertial-scale dynamics in the heliosheath is not yet thorough. Magnetic field fluctuations across the inner heliosheath (IHS) and the local interstellar medium (LISM) are here considered to provide accurate and highly resolved statistics over different plasma conditions between 88 and 136 au. By using the unique in situ 48 s measurements from the Voyager Interstellar Mission, we investigate different fluctuation regimes at the magnetohydrodynamic (MHD) scales, down to the MHD-to-kinetic transition. We focus on a range of scales exceeding five frequency decades (5 x 10-8 < f < 10-2 Hz), which is unprecedented in literature analysis. A set of magnetic field data for eight intervals in the IHS, in both unipolar and sector regions, and four intervals in the LISM is used for the analysis. Results are set forth in terms of the power spectral density, spectral compressibility, structure functions, and intermittency of magnetic field increments. In the heliosheath, we identify the energy-injection regime displaying a ~1/f energy decay, and the inertial-cascade regime. Here, the power spectrum is anisotropic and dominated by compressive modes, with intermittency that can reach kurtosis values of up to 10. In the interstellar medium the structure of turbulence is anisotropic as well, with transverse fluctuations clearly prevailing after 2015 May. Here, we show that intermittent features occur only at scales smaller than 10-6 Hz.

We demonstrate how the parameters of a Gibson-Low flux-rope-based coronal mass ejection (CME) can be constrained using remote observations. Our Multi Scale Fluid-Kinetic Simulation Suite (MS-FLUKSS) has been used to simulate the propagation of a CME in a data driven solar corona background computed using the photospheric magnetogram data. We constrain the CME model parameters using the observations of such key CME properties as its speed, orientation, and poloidal flux. The speed and orientation are estimated using multi-viewpoint white-light coronagraph images. The reconnected magnetic flux in the area covered by the post eruption arcade is used to estimate the poloidal flux in the CME flux rope. We simulate the partial halo CME on 7 March 2011 to demonstrate the efficiency of our approach. This CME erupted with the speed of 812 km/s and its poloidal flux, as estimated from source active region data, was 4.9e21 Mx. Using our approach, we were able to simulate this CME with the speed 840 km/s and the poloidal flux of 5.1e21 Mx, in remarkable agreement with the observations.

Coronal Mass Ejections (CMEs) are major drivers of extreme space weather conditions, this being a matter of serious concern for our modern technologically-dependent society. Development of numerical approaches that would simulate CME generation and propagation through the interplanetary space is an important step towards our capability to predict CME arrival times at Earth and their geo-effectiveness. In this paper, we utilize a data-constrained Gibson--Low (GL) flux rope model to generate CMEs. We derive the geometry of the initial GL flux rope using the Graduated Cylindrical Shell (GCS) method. This method uses multiple viewpoints from STEREO A & B Cor1/Cor2, and SOHO/LASCO C2/C3 coronagraphs to determine the size and orientation of a CME flux rope as it starts to erupt from the Sun. A flux rope generated in this way is inserted into a quasi-steady global magnetohydrodynamics (MHD) background solar wind flow driven by SDO/HMI line-of-sight magnetogram data, and erupts immediately. Numerical results obtained with the Multi-Scale Fluid-Kinetic Simulation Suite (MS-FLUKSS) code are compared with STEREO and SOHO/LASCO coronagraph observations in particular in terms of the CME speed, acceleration, and magnetic field structure.

The heliosphere is formed due to interaction between the solar wind (SW) and local interstellar medium (LISM). The shape and position of the heliospheric boundary, the heliopause, in space depend on the parameters of interacting plasma flows. The interplay between the asymmetrizing effect of the interstellar magnetic field and charge exchange between ions and neutral atoms plays an important role in the SW-LISM interaction. By performing three-dimensional, MHD plasma / kinetic neutral atom simulations, we determine the width of the outer heliosheath - the LISM plasma region affected by the presence of the heliosphere - and analyze quantitatively the distributions in front of the heliopause. It is shown that charge exchange modifies the LISM plasma to such extent that the contribution of a shock transition to the total variation of plasma parameters becomes small even if the LISM velocity exceeds the fast magnetosonic speed in the unperturbed medium. By performing adaptive mesh refinement simulations, we show that a distinct boundary layer of decreased plasma density and enhanced magnetic field should be observed on the interstellar side of the heliopause. We show that this behavior is in agreement with the plasma oscillations of increasing frequency observed by the plasma wave instrument onboard Voyager 1. We also demonstrate that Voyager observations in the inner heliosheath between the heliospheric termination shock and the heliopause are consistent with dissipation of the heliospheric magnetic field. The choice of LISM parameters in this analysis is based on the simulations that fit observations of energetic neutral atoms performed by IBEX.

The magnetometer (MAG) on Voyager 1 (V1) has been sampling the interstellar magnetic field (ISMF) since August 2012. The V1 MAG observations have shown draped ISMF in the very local interstellar medium disturbed occasionally by significant enhancements in magnetic field strength. Using a three-dimensional, data driven, multi-fluid model, we investigated these magnetic field enhancements beyond the heliopause that are supposedly associated with solar transients. To introduce time-dependent effects at the inner boundary at 1 astronomical unit, we used daily averages of the solar wind parameters from the OMNI data set. The model ISMF strength, direction, and proton number density are compared with V1 data beyond the heliopause. The model reproduced the large-scale fluctuations between 2012.652 and 2016.652, including major events around 2012.9 and 2014.6. The model also predicts shocks arriving at V1 around 2017.395 and 2019.502. Another model driven by OMNI data with interplanetary coronal mass ejections (ICMEs) removed at the inner boundary suggests that ICMEs may play a significant role in the propagation of shocks into the interstellar medium.

We have developed a data-driven magnetohydrodynamic (MHD) model of the global solar corona which uses characteristically-consistent boundary conditions (BCs) at the inner boundary. Our global solar corona model can be driven by different observational data including Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI) synoptic vector magnetograms together with the horizontal velocity data in the photosphere obtained by the time-distance helioseismology method, and the line-of-sight (LOS) magnetogram data obtained by HMI, Solar and Heliospheric Observatory/Michelson Doppler Imager (SOHO/MDI), National Solar Observatory/Global Oscillation Network Group (NSO/GONG) and Wilcox Solar Observatory (WSO). We implemented our model in the Multi-Scale Fluid-Kinetic Simulation Suite (MS-FLUKSS) - a suite of adaptive mesh refinement (AMR) codes built upon the Chombo AMR framework developed at the Lawrence Berkeley National Laboratory. We present an overview of our model, characteristic BCs, and two results we obtained using our model: A benchmark test of relaxation of a dipole field using characteristic BCs, and relaxation of an initial PFSS field driven by HMI LOS magnetogram data, and horizontal velocity data obtained by the time-distance helioseismology method using a set of non-characteristic BCs.

The aurorae of Uranus were recently detected in the far ultraviolet with the Hubble Space Telescope (HST) providing a new, so far unique, means to remotely study the asymmetric Uranian magnetosphere from Earth. We analyze here two new HST Uranus campaigns executed in Sept. 2012 and Nov. 2014 with different temporal coverage and under variable solar wind conditions numerically predicted by three different MHD codes. Overall, HST images taken with the Space Telescope Imaging Spectrograph reveal auroral emissions in three pairs of successive images (one pair acquired in 2012 and two in 2014), hence six additional auroral detections in total, including the most intense Uranian aurorae ever seen with HST. The detected emissions occur close the expected arrival of interplanetary shocks. They appear as extended spots at southern latitudes, rotating with the planet. They radiate 5-24 kR and 1.3-8.8 GW of ultraviolet emission from H$_2$, last for tens of minutes and vary on timescales down to a few seconds. Fitting the 2014 observations with model auroral ovals constrains the longitude of the southern (northern) magnetic pole to 104+\-26{\deg} (284-\-26{\deg}) in the Uranian Longitude System. We suggest that the Uranian near-equinoctial aurorae are pulsed cusp emissions possibly triggered by large-scale magnetospheric compressions.

This paper summarizes the results obtained by the team "Heliosheath Processes and the Structure of the Heliopause: Modeling Energetic Particles, Cosmic Rays, and Magnetic Fields" supported by the International Space Science Institute in Bern, Switzerland. We focus on the physical processes occurring in the outer heliosphere, especially at its boundary called the heliopause (HP), and in the LISM. The importance of magnetic field, charge exchange between atoms and ions, and solar cycle on the heliopause topology and observed heliocentric distances to different heliospheric discontinuities are discussed. It is shown that time-dependent boundary conditions are necessary to describe the heliospheric asymmetries detected by the Voyager spacecraft. We also discuss the structure of the HP, especially due to its instability and magnetic reconnection. It is demonstrated that the Rayleigh-Taylor instability of the nose of the HP creates consecutive layers of the interstellar and heliospheric plasma which are magnetically connected to different sources. This may be a possible explanation of abrupt changes in the galactic and anomalous cosmic ray fluxes observed by Voyager 1 when it was crossing the HP structure for a period of about one month in the summer of 2012. This paper also discusses the plausibility of fitting simulation results to a number of observational data sets obtained by in situ and remote measurements. The distribution of magnetic field in the vicinity of the HP is discussed in the context of Voyager measurements. We discuss the transport of energetic particles in the inner and outer heliosheath, concentrating on the anisotropic spatial diffusion diffusion tensor and the pitch-angle dependence of perpendicular diffusion and demonstrate that the latter can explain the observed pitch-angle anisotropies of both the anomalous and galactic cosmic rays in the outer heliosheath.

We performed numerical calculations to test the suggestion by Desiati & Lazarian (2013) that the anisotropies of TeV cosmic rays may arise from their interactions with the heliosphere. For this purpose, we used a magnetic field model of the heliosphere and performed direct numerical calculations of particle trajectories. Unlike earlier papers testing the idea, we did not employ time-reversible techniques that are based on Liouville's theorem. We showed numerically that for scattering by the heliosphere the conditions of Liouville's theorem are not satisfied and the adiabatic approximation and time-reversibility of the particle trajectories are not valid. Our results indicate sensitivity to the magnetic structure of the heliospheric magnetic field, and we expect that this will be useful for probing this structure in future research.

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.

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.

The Interstellar Boundary Explorer (IBEX) spacecraft is providing the first all-sky maps of the energetic neutral atoms (ENAs) produced by charge-exchange between interstellar neutral \HI\ atoms and heliospheric solar wind and pickup ions in the heliosphere boundary regions. The 'edge' of the interstellar cloud presently surrounding the heliosphere extends less than 0.1 pc in the upwind direction, terminating at an unknown distance, indicating that the outer boundary conditions of the heliosphere could change during the lifetime of the IBEX satellite. Using reasonable values for future outer heliosphere boundary conditions, ENA fluxes are predicted for one possible source of ENAs coming from outside of the heliopause. The ENA production simulations use three-dimensional MHD plasma models of the heliosphere that include a kinetic description of neutrals and a Lorentzian distribution for ions. Based on this ENA production model, it is then shown that the sensitivities of the IBEX 1.1 keV skymaps are sufficient to detect the variations in ENA fluxes that are expected to accompany the solar transition into the next upwind cloud. Approximately 20% of the IBEX 1.1 keV pixels appear capable of detecting the predicted model differences at the $ 3 \sigma$ level, with these pixels concentrated in the Ribbon region. Regardless of the detailed ENA production model, the success of the modeled \BdotR\ directions in reproducing the Ribbon locus, together with our results, indicate that the Ribbon phenomenon traces the variations in the heliosphere distortion caused by the relative pressures of the interstellar magnetic and gaseous components.

The distribution of interstellar dust grains (ISDG) observed in the Solar System depends on the nature of the interstellar medium-solar wind interaction. The charge of the grains couples them to the interstellar magnetic field (ISMF) resulting in some fraction of grains being excluded from the heliosphere while grains on the larger end of the size distribution, with gyroradii comparable to the size of the heliosphere, penetrate the termination shock. This results in a skewing the size distribution detected in the Solar System. We present new calculations of grain trajectories and the resultant grain density distribution for small ISDGs propagating through the heliosphere. We make use of detailed heliosphere model results, using three-dimensional (3-D) magnetohydrodynamic/kinetic models designed to match data on the shape of the termination shock and the relative deflection of interstellar neutral H and He flowing into the heliosphere. We find that the necessary inclination of the ISMF relative to the inflow direction results in an asymmetry in the distribution of the larger grains (0.1 micron) that penetrate the heliopause. Smaller grains (0.01 micron) are completely excluded from the Solar System at the heliopause.

We investigate heliosheath energetic neutral atom (ENA) fluxes at keV energies, by assuming that the heliosheath proton distribution can be approximated by a kappa-distribution. The choice of the kappa parameter derives from observational data of the solar wind (SW). This has direct applications to the upcoming IBEX mission. We will look at all-sky ENA maps within the IBEX energy range (10 eV to 6 keV), as well as ENA energy spectra in several directions. We find that the use of kappa, as opposed to a Maxwellian, gives rise to greatly increased ENA fluxes above 1 keV, while medium energy fluxes are somewhat reduced. We show how IBEX data can be used to estimate the spectral slope in the heliosheath, and that the use of kappa reduces the differences between ENA maps at different energies. We also investigate the effect introducing a kappa-distribution has on the global interaction between the SW and the local interstellar medium (LISM), and find that there is generally an increase in energy transport from the heliosphere into the LISM, due to the modified profile of ENA's energies. This results in a termination shock that moves out by 4 AU, a heliopause that moves in by 9 AU and a bow shock 25 AU farther out, in the nose direction.

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.

We study axisymmetric accretion of slowly rotating matter onto a gravitation center. The process of the flow restructuring is presented as the angular velocity of accreting gas approaches the Keplerian velocity. For spherically-symmetric accretion we found conditions of existence of steady-state supersonic solutions for various external boundary conditions. These conditions are used for choosing proper initial and boundary conditions for accretion of the rotating matter. The study is performed on the basis of the numerical solution of the Euler gas dynamic equations using the second order extension of the Lax-Friedrichs scheme with the implicit approximation of the source term.