Traditionally, the solar magnetic field has been considered to have a negligible effect in the outer regions of the heliosphere. Recent works have shown that the solar magnetic field may play a crucial role in collimating the plasma in the heliosheath. Interstellar Boundary Explorer (IBEX) observations of the heliotail indicated a latitudinal structure varying with energy in the energetic neutral atom (ENA) fluxes. At energies ~1 keV, the ENA fluxes show an enhancement at low latitudes and a deficit of ENAs near the poles. At energies >2.7 keV, ENA fluxes had a deficit within low latitudes, and lobes of higher ENA flux near the poles. This ENA structure was initially interpreted to be a result of the latitudinal profile of the solar wind during solar minimum. We extend the work of Kornbleuth et al. (2018) by using solar minimum-like conditions and the recently developed SHIELD model. The SHIELD model couples the magnetohydrodynamic (MHD) plasma solution with a kinetic description of neutral hydrogen. We show that while the latitudinal profile of the solar wind during solar minimum contributes to the lobes in ENA maps, the collimation by the solar magnetic field is important in creating and shaping the two high latitude lobes of enhanced ENA flux observed by IBEX. This is the first work to explore the effect of the changing solar magnetic field strength on ENA maps. Our findings suggest that IBEX is providing the first observational evidence of the collimation of the heliosheath plasma by the solar magnetic field.

Neutral hydrogen has been shown to greatly impact the plasma flow in the heliopshere and the location of the heliospheric boundaries. We present the results of the Solar-wind with Hydrogen Ion Exchange and Large-scale Dynamics (SHIELD) model, a new, self-consistent, kinetic-MHD model of the outer heliosphere within the Space Weather Modeling Framework. The charge-exchange mean free path is on order of the size of the heliosphere; therefore, the neutral atoms cannot be described as a fluid. The SHIELD model couples the MHD solution for a single plasma fluid to the kinetic solution from for neutral hydrogen atoms streaming through the system. The kinetic code is based on the Adaptive Mesh Particle Simulator (AMPS), a Monte Carlo method for solving the Boltzmann equation. The SHIELD model accurately predicts the increased filtration of interstellar neutrals into the heliosphere. In order to verify the correct implementation within the model, we compare the results of the SHIELD model to other, well-established kinetic-MHD models. The SHIELD model matches the neutral hydrogen solution of these studies as well as the shift in all heliospheric boundaries closer to the Sun in comparison the the multi-fluid treatment of the neutral hydrogen atoms. Overall the SHIELD model shows excellent agreement to these models and is a significant improvement to the fluid treatment of interstellar hydrogen.

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

This textbook gives a perspective of heliophysics in a way that emphasizes universal processes from a perspective that draws attention to what provides Earth (and similar (exo-)planets) with a relatively stable setting in which life as we know it can thrive. The book is intended for students in physical sciences in later years of their university training and for beginning graduate students in fields of solar, stellar, (exo-)planetary, and planetary-system sciences.

The study of the acceleration of particles is an essential element of research in the heliospheric science. Here, we discuss the predisposition to the particle acceleration around coronal mass ejections (CMEs)-driven shocks with corrugated wave-like features. We adopt these attributes on shocks formed from disturbances due to the bimodal solar wind, CME deflection, irregular CME expansion, and the ubiquitous fluctuations in the solar corona. In order to understand the role of a wavy shock in particle acceleration, we define three initial smooth shock morphologies each one associated with a fast CME. Using polar Gaussian profiles we model these shocks in the low corona. We establish the corrugated appearance on smooth shock by using combinations of wave-like functions that represent the disturbances from medium and CME piston. For both shock types, smooth and corrugated, we calculate the shock normal angles between the shock normal and the radial upstream coronal magnetic field in order to classify the quasi-parallel and quasi-perpendicular regions. We consider that corrugated shocks are predisposed to different process of particle acceleration due to irregular distributions of shock normal angles around of the shock. We suggest that disturbances due to CME irregular expansion may be a decisive factor in origin of particle acceleration. Finally, we regard that accepting these features on shocks may be the start point for investigating some questions in the sheath and shock, like downstream-jets, instabilities, shock thermalization, shock stability, and injection particle process.

The shape of the solar wind bubble within the interstellar medium, the so-called heliosphere, has been explored over six decades. As the Sun moves through the surrounding partially-ionized medium, neutral hydrogen atoms penetrate the heliosphere, and through charge-exchange with the supersonic solar wind, create a population of hot pick-up ions (PUIs). The Termination Shock (TS) crossing by Voyager 2 (V2) data demonstrated that the heliosheath (HS) (the region of shocked solar wind) pressure is dominated by suprathermal particles. Here we use a novel magnetohydrodynamic model that treats the freshly ionized PUIs as a separate fluid from the thermal component of the solar wind. Unlike previous models, the new model reproduces the properties of the PUIs and solar wind ions based on the New Horizon and V2 spacecraft observations. The PUIs charge exchange with the cold neutral H atoms of the ISM in the HS and are quickly depleted. The depletion of PUIs cools the heliosphere downstream of the TS, "deflating" it and leading to a narrower HS and a smaller and rounder shape, in agreement with energetic neutral atom observations by the Cassini spacecraft. The new model, with interstellar magnetic field orientation constrained by the IBEX ribbon, reproduces the magnetic field data outside the HP at Voyager 1(V1). We present the predictions for the magnetic field outside the HP at V2.

A recent study by Opher et al. (2015) suggested the heliosphere has a "croissant" shape, where the heliosheath plasma is confined by the toroidal solar magnetic field. The "croissant" heliosphere is in contrast to the classically accepted view of a comet-like tail. We investigate the effect of the "croissant" heliosphere model on energetic neutral atom (ENA) maps. Regardless of the existence of a split tail, the confinement of the heliosheath plasma should appear in ENA maps. ENA maps from the Interstellar Boundary Explorer (IBEX) have shown two high latitude lobes with excess ENA flux at higher energies in the tail of the heliosphere. These lobes could be a signature of the confinement of the heliosheath plasma, while some have argued they are caused by the fast/slow solar wind profile. Here we present ENA maps of the "croissant" heliosphere, focusing on understanding the effect of the heliosheath plasma collimation by the solar magnetic field while using a uniform solar wind. We incorporate pick-up ions (PUIs) into our model based on Malama et al. (2006) and Zank et al. (2010). We use the neutral solution from our MHD model to determine the angular variation of the PUIs, and include the extinction of PUIs in the heliosheath. In the presence of a uniform solar wind, we find that the collimation in the "croissant" heliosphere does manifest itself into two high latitude lobes of increased ENA flux in the downwind direction.

Based on the difference between the orientation of the interstellar $B_{ISM}$ and the solar magnetic fields, there was an expectation that the magnetic field direction would rotate dramatically across the heliopause (HP). However, the Voyager 1 spacecraft measured very little rotation across the HP. Previously we showed that the $B_{ISM}$ twists as it approaches the HP and acquires a strong T component (East-West). Here we establish that reconnection in the eastern flank of the heliosphere is responsible for the twist. On the eastern flank the solar magnetic field has twisted into the positive N direction and reconnects with the Southward pointing component of the $B_{ISM}$. Reconnection drives a rotational discontinuity (RD) that twists the $B_{ISM}$ into the -T direction and propagates upstream in the interstellar medium towards the nose. The consequence is that the N component of $B_{ISM}$ is reduced in a finite width band upstream of the HP. Voyager 1 currently measures angles ($\delta=sin^{-1}(B_{N}/B)$) close to solar values. We present MHD simulations to support this scenario, suppressing reconnection in the nose region while allowing it in the flanks, consistent with recent ideas about reconnection suppression from diamagnetic drifts. The jump in plasma $\beta$ (the plasma to magnetic pressure) across the nose of HP is much greater than in the flanks because the heliosheath $\beta$ is greater there than in the flanks. Large-scale reconnection is therefore suppressed in the nose but not at the flanks. Simulation data suggest that $B_{ISM}$ will return to its pristine value $10-15~AU$ past the HP.

The misalignment of the solar rotation axis and the magnetic axis of the Sun produces a periodic reversal of the Parker spiral magnetic field and the sectored solar wind. The compression of the sectors is expected to lead to reconnection in the heliosheath (HS). We present particle-in-cell simulations of the sectored HS that reflect the plasma environment along the Voyager 1 and 2 trajectories, specifically including unequal positive and negative azimuthal magnetic flux as seen in the Voyager data \citep{Burlaga03}. Reconnection proceeds on individual current sheets until islands on adjacent current layers merge. At late time bands of the dominant flux survive, separated by bands of deep magnetic field depletion. The ambient plasma pressure supports the strong magnetic pressure variation so that pressure is anti-correlated with magnetic field strength. There is little variation in the magnetic field direction across the boundaries of the magnetic depressions. At irregular intervals within the magnetic depressions are long-lived pairs of magnetic islands where the magnetic field direction reverses so that spacecraft data would reveal sharp magnetic field depressions with only occasional crossings with jumps in magnetic field direction. This is typical of the magnetic field data from the Voyager spacecraft \citep{Burlaga11,Burlaga16}. Voyager 2 data reveals that fluctuations in the density and magnetic field strength are anti-correlated in the sector zone as expected from reconnection but not in unipolar regions. The consequence of the annihilation of subdominant flux is a sharp reduction in the "number of sectors" and a loss in magnetic flux as documented from the Voyager 1 magnetic field and flow data \citep{Richardson13}.

To accurately predict the space weather effects of coronal mass ejection (CME) impacts at Earth one must know if and when a CME will impact Earth, and the CME parameters upon impact. Kay et al. (2015b) presents Forecasting a CME's Altered Trajectory (ForeCAT), a model for CME deflections based on the magnetic forces from the background solar magnetic field. Knowing the deflection and rotation of a CME enables prediction of Earth impacts, and the CME orientation upon impact. We first reconstruct the positions of the 2008 April 10 and the 2012 July 12 CMEs from the observations. The first of these CMEs exhibits significant deflection and rotation (34 degrees deflection and 58 degrees rotation), while the second shows almost no deflection or rotation (<3 degrees each). Using ForeCAT, we explore a range of initial parameters, such as the CME location and size, and find parameters that can successfully reproduce the behavior for each CME. Additionally, since the deflection depends strongly on the behavior of a CME in the low corona (Kay et al. (2015a, 2015b)), we are able to constrain the expansion and propagation of these CMEs in the low corona.

Solar coronal mass ejections (CMEs) produce adverse space weather effects at Earth. Planets in the close habitable zone of magnetically active M dwarfs may experience more extreme space weather than at Earth, including frequent CME impacts leading to atmospheric erosion and leaving the surface exposed to extreme flare activity. Similar erosion may occur for hot Jupiters with close orbits around solar-like stars. We have developed a model, Forecasting a CME's Altered Trajectory (ForeCAT), which predicts a CME's deflection. We adapt ForeCAT to simulate CME deflections for the mid-type M dwarf V374 Peg and hot Jupiters with solar-type hosts. V374 Peg's strong magnetic fields can trap CMEs at the M dwarfs's Astrospheric Current Sheet, the location of the minimum in the background magnetic field. Solar-type CMEs behave similarly, but have much smaller deflections and do not get trapped at the Astrospheric Current Sheet. The probability of planetary impact decreases with increasing inclination of the planetary orbit with respect to the Astrospheric Current Sheet - 0.5 to 5 CME impacts per day for M dwarf exoplanets, 0.05 to 0.5 CME impacts per day for solar-type hot Jupiters. We determine the minimum planetary magnetic field necessary to shield a planet's atmosphere from the CME impacts. M dwarf exoplanets require values between tens and hundreds of Gauss. Hot Jupiters around a solar-type star, however, require a more reasonable <30 G. These values exceed the magnitude required to shield a planet from the stellar wind, suggesting CMEs may be the key driver of atmospheric losses.

Fluctuations in the flow velocity and magnetic fields are ubiquitous in the Solar System. These fluctuations are turbulent, in the sense that they are disordered and span a broad range of scales in both space and time. The study of solar wind turbulence is motivated by a number of factors all keys to the understanding of the Solar Wind origin and thermodynamics. The solar wind spectral properties are far from uniformity and evolve with the increasing distance from the sun. Most of the available spectra of solar wind turbulence were computed at 1 astronomical unit, while accurate spectra on wide frequency ranges at larger distances are still few. In this paper we consider solar wind spectra derived from the data recorded by the Voyager 2 mission during 1979 at about 5 AU from the sun. Voyager 2 data are an incomplete time series with a voids/signal ratio that typically increases as the spacecraft moves away from the sun (45% missing data in 1979), making the analysis challenging. In order to estimate the uncertainty of the spectral slopes, different methods are tested on synthetic turbulence signals with the same gap distribution as V2 data. Spectra of all variables show a power law scaling with exponents between -2.1 and -1.1, depending on frequency subranges. Probability density functions (PDFs) and correlations indicate that the flow has a significant intermittency.

Understanding the trajectory of a coronal mass ejection (CME), including any deflection from a radial path, and the orientation of its magnetic field is essential for space weather predictions. Kay et al. (2015b) developed a model, Forecasting a CME's Altered Trajectory (ForeCAT), of CME deflections and rotation due to magnetic forces, not including the effects of reconnection. ForeCAT is able to reproduce the deflection of observed CMEs (Kay et al. 2015a). The deflecting CMEs tend to show a rapid increase of their angular momentum close to the Sun, followed by little to no increase at farther distances. Here we quantify the distance at which the CME deflection is "determined," which we define as the distance after which the background solar wind has negligible influence on the total deflection. We consider a wide range in CME masses and radial speeds and determine that the deflection and rotation of these CMEs can be well-described by assuming they propagate with constant angular momentum beyond 10 Rs. The assumption of constant angular momentum beyond 10 Rs yields underestimates of the total deflection at 1 AU of only 1% to 5% and underestimates of the rotation of 10%. Since the deflection from magnetic forces is determined by 10 Rs, non-magnetic forces must be responsible for any observed interplanetary deflections or rotations where the CME has increasing angular momentum.

Plasma velocity and magnetic field measurements from the Voyager 2 mission are used to study solar wind turbulence in the slow solar wind at two different heliocentric distances, 5 and 29 astronomical units, sufficiently far apart to provide information on the radial evolution of this turbulence. The magnetic helicity and the cross-helicity, which express the correlation between the plasma velocity and the magnetic field, are used to characterize the turbulence. Wave number spectra are computed by means of the Taylor hypothesis applied to time resolved single point Voyager 2 measurements. The overall picture we get is complex and difficult to interpret. A substantial decrease of the cross-helicity at smaller scales (over 1-3 hours of observation) with increasing heliocentric distance is observed. At 5 AU the only peak in the probability density of the normalized residual energy is negative, near -0.5. At 29 AU the probability density becomes doubly peaked, with a negative peak at -0.5 and a smaller peak at a positive values of about 0.7. A decrease of the cross-helicity for increasing heliocentric distance is observed, together with a reduction of the unbalance toward the magnetic energy of the energy of the fluctuations. For the smaller scales, we found that at 29 AU the normalized polarization is small and positive on average (about 0.1), it is instead zero at 5 AU. For the larger scales, the polarization is low and positive at 5 AU (average around 0.1) while it is negative (around - 0.15) at 29 AU.

The Voyager probes are the furthest, still active, spacecraft ever launched from Earth. During their 38-year trip, they have collected data regarding solar wind properties (such as the plasma velocity and magnetic field intensity). Unfortunately, a complete time evolution of the measured physical quantities is not available. The time series contains many gaps which increase in frequency and duration at larger distances. The aim of this work is to perform a spectral and statistical analysis of the solar wind plasma velocity and magnetic field using Voyager 2 data measured in 1979, when the gaps/signal ratio is of order of unity. This analysis is achieved using four different data reconstruction techniques: averages on linearly interpolated subsets, correlation of linearly interpolated data, compressed sensing spectral estimation, and maximum likelihood data reconstruction. With five frequency decades, the spectra we obtained have the largest frequency range ever computed at 5 astronomical units from the Sun; spectral exponents have been determined for all the components of the velocity and magnetic field fluctuations. Void analysis is also useful in recovering other spectral properties such as integral scales (see for instance Table 4) and, if the confidence level of the measurements is sufficiently high, the decay variation in the small scale range due, for instance, to dissipative effects.

An analytic model of the heliosheath (HS) between the termination shock (TS) and the heliopause (HP) is developed in the limit in which the interstellar flow and magnetic field are neglected. The heliosphere in this limit is axisymmetric and the overall structure of the HS and HP are controlled by the solar magnetic field even in the limit in which the ratio of the plasma to magnetic field pressure, $\beta=8\pi P/B^2$, in the HS is large. The tension of the solar magnetic field produces a drop in the total pressure between the TS and the HP. This same pressure drop accelerates the plasma flow downstream of the TS into the North and South directions to form two collimated jets. The radii of these jets are controlled by the flow through the TS and the acceleration of this flow by the magnetic field -- a stronger solar magnetic field boosts the velocity of the jets and reduces the radii of the jets and the HP. Magnetohydrodynamic (MHD) simulations of the global helioshere embedded in a stationary interstellar medium match well with the analytic model. The results suggest that mechanisms that reduce the HS plasma pressure downstream of the TS can enhance the jet outflow velocity and reduce the HP radius to values more consistent with the Voyager 1 observations than in current global models.

Accurate space weather forecasting requires knowledge of the trajectory of coronal mass ejections (CMEs), including any deflections close to the Sun or through interplanetary space. Kay et al. 2013 introduced ForeCAT, a model of CME deflection resulting from the background solar magnetic field. For a magnetic field solution corresponding to Carrington Rotation (CR) 2029 (declining phase, April-May 2005), the majority of the CMEs deflected to the Heliospheric Current Sheet (HCS), the minimum in magnetic pressure on global scales. Most of the deflection occurred below 4 Rs. Here we extend ForeCAT to include a three dimensional description of the deflecting CME. We attempt to answer the following questions: a) Do all CMEs deflect to the magnetic minimum? and b) Does most deflection occur within the first few solar radii (~4 Rs)? Results for solar minimum and declining phase CMEs show that not every CME deflects to the magnetic minimum and that the deflection is typically determined below 2 Rs. Slow, wide, low mass CMEs in declining phase solar backgrounds with strong magnetic field and magnetic gradients exhibit the largest deflections. Local gradients related to active regions tend to cause the largest deviations from the deflection predicted by global magnetic gradients, but variations can also be seen for CMEs in the quiet sun regions of the declining phase CR. We show the torques due to differential forces along the CME can cause rotation about the CME's toroidal axis.

Decades of observations show that CMEs can deflect from a purely radial trajectory yet no consensus exists as to the cause of these deflections. Many of theories attribute the CME deflection to magnetic forces. We developed ForeCAT (Kay et al. 2013, Kay et al. 2015), a model for CME deflections based solely on magnetic forces, neglecting any reconnection effects. Here we compare ForeCAT predictions to the observed deflection of the 2008 December 12 CME and find that ForeCAT can accurately reproduce the observations. Multiple observations show that this CME deflected nearly 30{\deg} in latitude (Byrne et al. 2010, Gui et al. 2011) and 4.4{\deg} in longitude (Gui et al. 2011). From the observations, we are able to constrain all of the ForeCAT input parameters (initial position, radial propagation speed, and expansion) except the CME mass and the drag coefficient that affects the CME motion. By minimizing the reduced chi-squared, $\chi^2_{\nu}$, between the ForeCAT results and the observations we determine an acceptable mass range between 4.5x10$^{14}$ and 1x10$^{15}$ g and the drag coefficient less than 1.4 with a best fit at 7.5x10$^{14}$ g and 0 for the mass and drag coefficient. ForeCAT is sensitive to the magnetic background and we are also able to constrain the rate at which the quiet sun magnetic field falls to be similar or to or fall slightly slower than the Potential Field Source Surface model.

The classic accepted view of the heliosphere is a quiescent, comet-like shape aligned in the direction of the Sun's travel through the interstellar medium (ISM) extending for 1000's of AUs (AU: astronomical unit). Here we show, based on magnetohydrodynamic (MHD) simulations, that the tension (hoop) force of the twisted magnetic field of the sun confines the solar wind plasma beyond the termination shock and drives jets to the North and South very much like astrophysical jets. These jets are deflected into the tail region by the motion of the Sun through the ISM similar to bent galactic jets moving through the intergalactic medium. The interstellar wind blows the two jets into the tail but is not strong enough to force the lobes into a single comet-like tail, as happens to some astrophysical jets (Morsony et al. 2013). Instead, the interstellar wind flows around the heliosphere and into equatorial region between the two jets. As in some astrophysical jets that are kink unstable (Porth et al. 2014) we show here that the heliospheric jets are turbulent (due to large-scale MHD instabilities and reconnection) and strongly mix the solar wind with the ISM beyond 400 AU. The resulting turbulence has important implications for particle acceleration in the heliosphere. The two-lobe structure is consistent with the energetic neutral atoms (ENAs) images of the heliotail from IBEX (McComas et al. 2013) where two lobes are visible in the North and South and the suggestion from the CASSINI (Krimigis et al. 2009, Dialynas et al. 2013) ENAs that the heliosphere is "tailless".

The acceleration of protons and electrons to high (sometimes GeV/nucleon) energies by solar phenomena is a key component of space weather. These solar energetic particle (SEP) events can damage spacecraft and communications, as well as present radiation hazards to humans. In-depth particle acceleration simulations have been performed for idealized magnetic fields for diffusive acceleration and particle propagation, and at the same time the quality of MHD simulations of coronal mass ejections (CMEs) has improved significantly. However, to date these two pieces of the same puzzle have remained largely decoupled. Such structures may contain not just a shock but also sizable sheath and pileup compression regions behind it, and may vary considerably with longitude and latitude based on the underlying coronal conditions. In this work, we have coupled results from a detailed global three-dimensional MHD time-dependent CME simulation to a global proton acceleration and transport model, in order to study time-dependent effects of SEP acceleration between 1.8 and 8 solar radii in the 2005 May 13 CME. We find that the source population is accelerated to at least 100 MeV, with distributions enhanced up to six orders of magnitude. Acceleration efficiency varies strongly along field lines probing different regions of the dynamically evolving CME, whose dynamics is influenced by the large-scale coronal magnetic field structure. We observe strong acceleration in sheath regions immediately behind the shock.

We perform three-dimensional numerical simulations of stellar winds of early-M dwarf stars. Our simulations incorporate observationally reconstructed large-scale surface magnetic maps, suggesting that the complexity of the magnetic field can play an important role in the angular momentum evolution of the star, possibly explaining the large distribution of periods in field dM stars, as reported in recent works. In spite of the diversity of the magnetic field topologies among the stars in our sample, we find that stellar wind flowing near the (rotational) equatorial plane carries most of the stellar angular momentum, but there is no preferred colatitude contributing to mass loss, as the mass flux is maximum at different colatitudes for different stars. We find that more non-axisymmetric magnetic fields result in more asymmetric mass fluxes and wind total pressures $p_{\rm tot}$ (defined as the sum of thermal, magnetic and ram pressures). Because planetary magnetospheric sizes are set by pressure equilibrium between the planet's magnetic field and $p_{\rm tot}$, variations of up to a factor of $3$ in $p_{\rm tot}$ (as found in the case of a planet orbiting at several stellar radii away from the star) lead to variations in magnetospheric radii of about 20 percent along the planetary orbital path. In analogy to the flux of cosmic rays that impact the Earth, which is inversely modulated with the non-axisymmetric component of the total open solar magnetic flux, we conclude that planets orbiting M dwarf stars like DT~Vir, DS~Leo and GJ~182, which have significant non-axisymmetric field components, should be the more efficiently shielded from galactic cosmic rays, even if the planets lack a protective thick atmosphere/large magnetosphere of their own.

There was an expectation that the magnetic field direction would rotate dramatically across the heliopause (HP). This rotation was used as one of the criteria to determine if Voyager 1 (V1) had crossed the HP. Recently the Voyager team concluded that V1 crossed into interstellar space last year (Gurnett et al. 2013). The question is then why there was no significant rotation in the direction of the magnetic field across the HP (Burlaga et al. 2013). Here we present simulations that reveal that strong rotations in the direction of the magnetic field at the HP at the location of V1 (and Voyager 2) are not expected. The solar magnetic field strongly affects the drapping of the interstellar magnetic field (BISM) around the HP. BISM twists as it approaches the HP and acquires a strong T component (east-west). The strong increase in the T component occurs where the interstellar flow stagnates in front of the HP. At this same location the N component BN is significantly reduced. Above and below the neighboring interstellar magnetic field lines also twist into the T direction. This behavior occurs for a wide range of orientations of BISM. Thus, the BN component at the V1 location and the associated angle delta=a sin(BN/B) are small (around 10-20 degrees), as seen in the observations (Burlaga et al. 2013). Only after some significant distance outside the HP, is the direction of the interstellar field distinguishably different from that of the Parker spiral.

To predict whether a coronal mass ejection (CME) will impact Earth, the effects of the background on the CME's trajectory must be taken into account. We develop a model, ForeCAT (Forecasting a CME's Altered Trajectory), of CME deflection due to magnetic forces. ForeCAT includes CME expansion, a three-part propagation model, and the effects of drag on the CME's deflection. Given the background solar wind conditions, the launch site of the CME, and the properties of the CME (mass, final propagation speed, initial radius, and initial magnetic strength), ForeCAT predicts the deflection of the CME. Two different magnetic backgrounds are considered: a scaled background based on type II radio burst profiles and a Potential Field Source Surface (PFSS) background. For a scaled background where the CME is launched from an active region located between a CH and streamer region the strong magnetic gradients cause a deflection of 8.1 degrees in latitude and 26.4 degrees in longitude for a 1e15 g CME propagating out to 1 AU. Using the PFSS background, which captures the variation of the streamer belt position with height, leads to a deflection of 1.6 degrees in latitude and 4.1 degrees in longitude for the control case. Varying the CME's input parameters within observed ranges leads to the majority of CMEs reaching the streamer belt within the first few solar radii. For these specific backgrounds, the streamer belt acts like a potential well that forces the CME into an equilibrium angular position.

The picture of the heliopause (HP) -- the boundary between the domains of the sun and the local interstellar medium (LISM) -- as a pristine interface with a large rotation in the magnetic field fails to describe recent Voyager 1 (V1) spacecraft data. Magnetohydrodynamic (MHD) simulations of the global heliosphere reveal that the rotation angle of the magnetic field across the HP at V1 is small. Particle-in-cell simulations, based on cuts through the MHD model at the location of V1, suggest that the sectored region of the heliosheath (HS) produces large-scale magnetic islands that reconnect with the interstellar magnetic field and mix LISM and HS plasma. Cuts across the simulation data reveal multiple, anti-correlated jumps in the number densities of LISM and HS particles at the magnetic separatrices of the islands, similar to those observed by V1. A model is presented, based on both the observations and simulation data, of the HP as a porous, multi-layered structure threaded by magnetic fields. This model further suggests that, contrary to the conclusions of recent papers, V1 has already crossed the HP.

After the termination shock (TS) crossing, the Voyager 2 spacecraft has been observing strong variations of the magnetic field and solar wind parameters in the heliosheath. Anomalous cosmic rays, electrons, and galactic cosmic rays present strong intensity fluctuations. Several works suggested that the fluctuations might be attributed to spatial variations within the heliosheath. Additionally, the variability of the solar wind in this region is caused by different temporal events that occur near the Sun and propagate to the outer heliosphere. To understand the spatial and temporal effects in the heliosheath, it is important to study these effects separately. In this work we explore the role of shocks as one type of temporal effects in the dynamics of the heliosheath. Although currently plasma in the heliosheath is dominated by solar minima conditions, with increasing solar cycle shocks associated with transients will play an important role. We used a 3D MHD multi-fluid model of the interaction between the solar wind and the local interstellar medium to study the propagation of a pair of forward-reverse shocks in the supersonic solar wind, interaction with the TS, and propagation to the heliosheath. We found that in the supersonic solar wind the interaction region between the shocks expands, the shocks weaken and decelerate. The fluctuation amplitudes of the plasma parameters vary with heliocentric distance. The interaction of the pair of shocks with the TS creates a variety of new waves and discontinuities in the heliosheath, which produce a highly variable solar wind flow. The collision of the forward shock with the heliopause causes a reflection of fast magnetosonic waves inside the heliosheath.

Data from the Voyager probes and the Interstellar Boundary Explorer have revealed the importance of pick-up ions (PUIs) in understanding the character and behavior of the outer heliosphere, the region of interaction between the solar wind and the interstellar medium. In the outer heliosphere PUIs carry a large fraction of the thermal pressure, which effects the nature of the termination shock, and they are a dominate component of pressure in the heliosheath. This paper describes the development of a new multi-ion, multi-fluid 3-D magnetohydrodynamic model of the outer heliosphere. This model has the added capability of tracking the individual fluid properties of multiple ion populations. For this initial study two ion populations are modeled: the thermal solar wind ions and PUIs produced in the supersonic solar wind. The model also includes 4 neutral fluids that interact through charge-exchange with the ion fluids. The new multi-ion simulation reproduces the significant heating of PUIs at the termination shock, as inferred from Voyager observations, and provides properties of PUIs in the 3-D heliosheath. The thinning of the heliosheath due to the loss of thermal energy in the heliosheath from PUI and neutral interaction is also quantified. In future work the multi-ion, multi-fluid model will be used to simulate energetic neutral atom (ENA) maps for comparison with the Interstellar Boundary Explorer, particularly at PUI energies of less than 1 keV.

Tau Boo is an intriguing planet-host star that is believed to undergo magnetic cycles similar to the Sun, but with a duration that is about one order of magnitude smaller than that of the solar cycle. With the use of observationally derived surface magnetic field maps, we simulate the magnetic stellar wind of Tau Boo by means of three-dimensional MHD numerical simulations. As the properties of the stellar wind depend on the particular characteristics of the stellar magnetic field, we show that the wind varies during the observed epochs of the cycle. Although the mass loss-rates we find (~2.7e-12 Msun/yr) vary less than 3 per cent during the observed epochs of the cycle, our derived angular momentum loss-rates vary from 1.1 to 2.2e32erg. The spin-down times associated to magnetic braking range between 39 and 78Gyr. We also compute the emission measure from the (quiescent) closed corona and show that it remains approximately constant through these epochs at a value of ~10^{50.6} cm^{-3}. This suggests that a magnetic cycle of Tau Boo may not be detected by X-ray observations. We further investigate the interaction between the stellar wind and the planet by estimating radio emission from the hot-Jupiter that orbits at 0.0462 au from Tau Boo. By adopting reasonable hypotheses, we show that, for a planet with a magnetic field similar to Jupiter (~14G at the pole), the radio flux is estimated to be about 0.5-1 mJy, occurring at a frequency of 34MHz. If the planet is less magnetised (field strengths roughly <4G), detection of radio emission from the ground is unfeasible due to the Earth's ionospheric cutoff. According to our estimates, if the planet is more magnetised than that and provided the emission beam crosses the observer line-of-sight, detection of radio emission from Tau Boo b is only possible by ground-based instruments with a noise level of < 1 mJy, operating at low frequencies.

We present a new 3d self-consistent two-component (plasma and neutral hydrogen) model of the solar wind interaction with the local interstellar medium (LISM). This model (K-MHD) combines the MHD treatment of the solar wind and the ionized LISM component, with a kinetic model of neutral interstellar hydrogen (LISH). The local interstellar magnetic field (BLISM) intensity and orientation are chosen based on an early analysis of the heliosheath flows (Opher et al. 2009). The properties of the plasma and neutrals obtained using the (K-MHD) model are compared to previous multi-fluid (Opher et al. 2009) and kinetic models (Izmodenov et al. 2005). The new treatment of LISH revealed important changes in the heliospheric properties not captures by the multi-fluid model. These include a decrease in the heliocentric distance to the termination shock (TS), a thinner heliosheath and a reduced deflection angle ({\theta}) of the heliosheath flows. The asymmetry of the termination shock, however, seems to be unchanged by the kinetic aspect of the LISH.

All the current global models of the heliosphere are based on the assumption that the magnetic field in the heliosheath, in the region close to the heliopause is laminar. We argue that in that region the heliospheric magnetic field is not laminar but instead consists of magnetic bubbles. Recently, we proposed that the annihilation of the "sectored" magnetic field within the heliosheath as it is compressed on its approach to the heliopause produces the anomalous cosmic rays and also energetic electrons. As a product of the annihilation of the sectored magnetic field, densely-packed magnetic islands/bubbles are produced. These magnetic islands/bubbles will be convected with the ambient flows as the sector region is carried to higher latitudes filling the heliosheath. We further argue that the magnetic islands/bubbles will develop upstream within the heliosheath. As a result, the magnetic field in the heliosheath sector region will be disordered well upstream of the heliopause. We present a 3D MHD simulation with very high numerical resolution that captures the north-south boundaries of the sector region. We show that due to the high pressure of the interstellar magnetic field a north-south asymmetry develops such that the disordered sectored region fills a large portion of the northern part of the heliosphere with a smaller extension in the southern hemisphere. We suggest that this scenario is supported by the following changes that occur around 2008 and from 2009.16 onward: a) the sudden decrease in the intensity of low energy electrons detected by Voyager 2; b) a sharp reduction in the intensity of fluctuations of the radial flow; and c) the dramatic differences in intensity trends between GCRs at V1 and 2. We argue that these observations are a consequence of V2 leaving the sector region of disordered field during these periods and crossing into a region of unipolar laminar field.

Recently, surface magnetic field maps had been acquired for a small sample of active M dwarfs, showing that fully convective stars (spectral types ~M4 and later) host intense (~kG), mainly axi-symmetrical poloidal fields. In particular, the rapidly rotating M dwarf V374Peg (M4), believed to lie near the theoretical full convection threshold, presents a stable magnetic topology on a time-scale of 1 yr. The rapid rotation of V374Peg (P=0.44 days) along with its intense magnetic field point toward a magneto-centrifugally acceleration of a coronal wind. In this work, we aim at investigating the structure of the coronal magnetic field in the M dwarf V374Peg by means of three-dimensional magnetohydrodynamical (MHD) numerical simulations of the coronal wind. For the first time, an observationally derived surface magnetic field map is implemented in MHD models of stellar winds for a low-mass star. We self-consistently take into consideration the interaction of the outflowing wind with the magnetic field and vice versa. Hence, from the interplay between magnetic forces and wind forces, we are able to determine the configuration of the magnetic field and the structure of the coronal winds. Our results enable us to evaluate the angular momentum loss of the rapidly rotating M dwarf V374Peg.

The rapid rotation (P=0.44 d) of the M dwarf V374Peg (M4) along with its intense magnetic field point toward magneto-centrifugal acceleration of a coronal wind. In this work, we investigate the structure of the wind of V374Peg by means of 3D magnetohydrodynamical (MHD) numerical simulations. For the first time, an observationally derived surface magnetic field map is implemented in MHD models of stellar winds for a low mass star. We show that the wind of V374Peg deviates greatly from a low-velocity, low-mass-loss rate solar-type wind. We find general scaling relations for the terminal velocities, mass-loss rates, and spin-down times of highly magnetized M dwarfs. In particular, for V374Peg, our models show that terminal velocities across a range of stellar latitudes reach ~(1500-2300) n_{12}^{-1/2} km/s, where n_{12} is the coronal wind base density in units of 10^{12} cm^{-3}, while the mass-loss rates are about 4 x 10^{-10} n_{12}^{1/2} Msun/yr. We also evaluate the angular-momentum loss of V374Peg, which presents a rotational braking timescale ~28 n_{12}^{-1/2} Myr. Compared to observationally derived values from period distributions of stars in open clusters, this suggests that V374Peg may have low coronal base densities (< 10^{11} cm^{-3}). We show that the wind ram pressure of V374Peg is about 5 orders of magnitude larger than for the solar wind. Nevertheless, a small planetary magnetic field intensity (~ 0.1G) is able to shield a planet orbiting at 1 AU against the erosive effects of the stellar wind. However, planets orbiting inside the habitable zone of V374Peg, where the wind ram pressure is higher, might be facing a more significant atmospheric erosion. In that case, higher planetary magnetic fields of, at least, about half the magnetic field intensity of Jupiter, are required to protect the planet's atmosphere.

Based on our previous work (Vidotto et al. 2009a), we investigate the effects on the wind and magnetospheric structures of weak-lined T Tauri stars due to a misalignment between the axis of rotation of the star and its magnetic dipole moment vector. In such configuration, the system loses the axisymmetry presented in the aligned case, requiring a fully 3D approach. We perform 3D numerical MHD simulations of stellar winds and study the effects caused by different model parameters. The system reaches a periodic behavior with the same rotational period of the star. We show that the magnetic field lines present an oscillatory pattern and that by increasing the misalignment angle, the wind velocity increases. Our wind models allow us to study the interaction of a magnetized wind with a magnetized extra-solar planet. Such interaction gives rise to reconnection, generating electrons that propagate along the planet's magnetic field lines and produce electron cyclotron radiation at radio wavelengths. We find that a close-in Jupiter-like planet orbiting at 0.05AU presents a radio power that is ~5 orders of magnitude larger than the one observed in Jupiter, which suggests that the stellar wind from a young star has the potential to generate strong planetary radio emission that could be detected in the near future with LOFAR. This radio power varies according to the phase of rotation of the star. We also analyze whether winds from misaligned stellar magnetospheres could cause a significant effect on planetary migration. Compared to the aligned case, we show that the time-scale tau_w for an appreciable radial motion of the planet is shorter for larger misalignment angles. While for the aligned case tau_w~100Myr, for a stellar magnetosphere tilted by 30deg, tau_w ranges from ~40 to 70Myr for a planet located at a radius of 0.05AU. (Abridged)

We propose that magnetic reconnection at the heliopause only occurs where the interstellar magnetic field points nearly anti-parallel to the heliospheric field. By using large-scale magnetohydrodynamic (MHD) simulations of the heliosphere to provide the initial conditions for kinetic simulations of heliopause (HP) reconnection we show that the energetic pickup ions downstream from the solar wind termination shock induce large diamagnetic drifts in the reconnecting plasma and stabilize non-anti-parallel reconnection. With this constraint the MHD simulations can show where HP reconnection most likely occurs. We also suggest that reconnection triggers the 2-3 kHz radio bursts that emanate from near the HP. Requiring the burst locations to coincide with the loci of anti-parallel reconnection allows us to determine, for the first time, the vector direction of the local interstellar magnetic field. We find it to be oriented towards the southern solar magnetic pole.

The recent observations of the anomalous cosmic ray (ACR) energy spectrum as Voyagers 1 and 2 crossed the heliospheric termination shock have called into question the conventional shock source of these energetic particles. We suggest that the sectored heliospheric magnetic field, which results from the flapping of the heliospheric current sheet, piles up as it approaches the heliopause, narrowing the current sheets that separate the sectors and triggering the onset of collisionless magnetic reconnection. Particle-in-cell simulations reveal that most of the magnetic energy is released and most of this energy goes into energetic ions with significant but smaller amounts of energy going into electrons. The energy gain of the most energetic ions results from their reflection from the ends of contracting magnetic islands, a first order Fermi process. The energy gain of the ions in contracting islands increases their parallel (to the magnetic field ${\bf B}$) pressure $p_\parallel$ until the marginal firehose condition is reached, causing magnetic reconnection and associated particle acceleration to shut down. The model calls into question the strong scattering assumption used to derive the Parker transport equation and therefore the absence of first order Fermi acceleration in incompressible flows. A simple 1-D model for particle energy gain and loss is presented in which the feedback of the energetic particles on the reconnection drive is included. The ACR differential energy spectrum takes the form of a power law with a spectral index slightly above 1.5. The model has the potential to explain several key Voyager observations, including the similarities in the spectra of different ion species.

Here we investigate the contribution of surface Alfven wave damping to the heating of the solar wind in minima conditions. These waves are present in regions of strong inhomogeneities in density or magnetic field (e. g., the border between open and closed magnetic field lines). Using a 3-dimensional Magnetohydrodynamics (MHD) model, we calculate the surface Alfven wave damping contribution between 1-4 solar radii, the region of interest for both acceleration and coronal heating. We consider waves with frequencies lower than those that are damped in the chromosphere and on the order of those dominating the heliosphere. In the region between open and closed field lines, within a few solar radii of the surface, no other major source of damping has been suggested for the low frequency waves we consider here. This work is the first to study surface Alfven waves in a 3D environment without assuming a priori a geometry of field lines or magnetic and density profiles. We determine that waves with frequencies >2.8x10^-4 Hz are damped between 1-4 solar radii. In quiet sun regions, surface Alfven waves are damped at further distances compared to active regions, thus carrying additional wave energy into the corona. We compare the surface Alfven wave contribution to the heating by a variable polytropic index and find that it an order of magnitude larger than needed for quiet sun regions. For active regions the contribution to the heating is twenty percent. As it has been argued that a variable gamma acts as turbulence, our results indicate that surface Alfven wave damping is comparable to turbulence in the lower corona. This damping mechanism should be included self consistently as an energy driver for the wind in global MHD models.

By means of numerical simulations, we investigate magnetized stellar winds of pre-main-sequence stars. In particular we analyze under which circumstances these stars will present elongated magnetic features (e.g., helmet streamers, slingshot prominences, etc). We focus on weak-lined T Tauri stars, as the presence of the tenuous accretion disk is not expected to have strong influence on the structure of the stellar wind. We show that the plasma-beta parameter (the ratio of thermal to magnetic energy densities) is a decisive factor in defining the magnetic configuration of the stellar wind. Using initial parameters within the observed range for these stars, we show that the coronal magnetic field configuration can vary between a dipole-like configuration and a configuration with strong collimated polar lines and closed streamers at the equator (multi-component configuration for the magnetic field). We show that elongated magnetic features will only be present if the plasma-beta parameter at the coronal base is beta<<1. Using our self-consistent 3D MHD model, we estimate for these stellar winds the time-scale of planet migration due to drag forces exerted by the stellar wind on a hot-Jupiter. In contrast to the findings of Lovelace et al. (2008), who estimated such time-scales using the Weber & Davis model, our model suggests that the stellar wind of these multi-component coronae are not expected to have significant influence on hot-Jupiters migration. Further simulations are necessary to investigate this result under more intense surface magnetic field strengths (~2-3 kG) and higher coronal base densities, as well as in a tilted stellar magnetosphere.

We discuss a model of cosmic ray acceleration that accounts for the observations of anomalous cosmic rays by Voyager 1 and 2. The model appeals to fast magnetic reconnection rather than shocks as the driver of acceleration. The ultimate source of energy is associated with magnetic field reversals that occur in the heliosheath. It is expected that the magnetic field reversals will occur throughout the heliosheath, but especially near the heliopause where the flows slows down and diverge in respect to the interstellar wind and also in the boundary sector-in the heliospheric current sheet. While the First Order Fermi acceleration theory within reconnection layers is in its infancy, the available predictions do not contradict the available data on anomalous cosmic ray spectra measured by the spacecrafts. We argue that the Voyager data can be one of the first pieces of evidence favoring the acceleration within regions of fast magnetic reconnection, which we argue is a widely spread astrophysical process.

By means of self-consistent 3D MHD numerical simulations, we analyze magnetized solar-like stellar winds and their dependence on the plasma-beta parameter. We adopt in our simulations a heating parameter described by gamma, which is responsible for the thermal acceleration of the wind. We analyze winds with polar magnetic field intensities ranging from 1 to 20G. We show that the wind structure presents characteristics that are similar to the solar coronal wind. The steady-state magnetic field topology for all cases is similar, presenting a configuration of helmet streamer-type, with zones of closed field lines and open field lines coexisting. Higher magnetic field intensities lead to faster and hotter winds. The increase of the field intensity generates a larger dead zone in the wind, i. e., the closed loops that inhibit matter to escape from latitudes lower than ~45 degrees extend farther away from the star. The Lorentz force leads naturally to a latitude-dependent wind. We show that by increasing the density and maintaining B0=20G, the system recover back to slower and cooler winds. For a fixed gamma, we show that the key parameter in determining the wind velocity profile is the beta-parameter at the coronal base. Therefore, there is a group of magnetized flows that would present the same terminal velocity despite of its thermal and magnetic energy densities, as long as the plasma-beta parameter is the same. This degeneracy, however, can be removed if we compare other physical parameters of the wind, such as the mass-loss rate. We analyze the influence of gamma in our results and we show that it is also important in determining the wind structure. (Abridged)

We investigate magnetized solar-like stellar winds by means of self-consistent three-dimensional (3D) magnetohydrodynamics (MHD) numerical simulations. We analyze winds with different magnetic field intensities and densities as to explore the dependence on the plasma-beta parameter. By solving the fully ideal 3D MHD equations, we show that the plasma-beta parameter is the crucial parameter in the configuration of the steady-state wind. Therefore, there is a group of magnetized flows that would present the same terminal velocity despite of its thermal and magnetic energy densities, as long as the plasma-beta parameter is the same.

Astrophysical fluids, including interstellar and interplanetary medium, are magnetized and turbulent. Their appearance, evolution, and overall properties are determined by the magnetic turbulence that stirs it. We argue that examining magnetic turbulence at a fundamental level is vital to understanding many processes. A point that frequently escapes the attention of researchers is that magnetic turbulence cannot be confidently understood only using "brute force" numerical approaches. In this review we illustrate this point on a number of examples, including intermittent heating of plasma by turbulence, interactions of turbulence with cosmic rays and effects of turbulence on the rate of magnetic reconnection. We show that the properties of magnetic turbulence may vary considerably in various environments, e.g. imbalanced turbulence in solar wind differs from balanced turbulence and both of these differ from turbulence in partially ionized gas. Appealing for the necessity of more observational data on magnetic fields, we discuss a possibility of studying interplanetary turbulence using alignment of Sodium atoms in the tail of comets.

We numerically model the coronal mass ejection (CME) event of October 28, 2003 that erupted from active region 10486 and propagated to Earth in less than 20 hours causing severe geomagnetic storms. The magnetohydrodynamic (MHD) model is formulated by first arriving at a steady state corona and solar wind employing synoptic magnetograms. We initiate two CMEs from the same active region, one approximately a day earlier that preconditions the solar wind for the much faster CME on the 28th. This second CME travels through the corona at a rate of over 2500 km s$^{-1}$ driving a strong forward shock. We clearly identify this shock in an image produced by the Large Angle Spectrometric Coronagraph (LASCO) C3, and reproduce the shock and its appearance in synthetic white light images from the simulation. We find excellent agreement with both the general morphology and the quantitative brightness of the model CME with LASCO observations. These results demonstrate that the CME shape is largely determined by its interaction with the ambient solar wind and may not be sensitive to the initiation process. We then show how the CME would appear as observed by wide-angle coronagraphs onboard the Solar Terrestrial Relations Observatory (STEREO) spacecraft. We find complex time evolution of the white-light images as a result of the way in which the density structures pass through the Thomson sphere. The simulation is performed with the Space Weather Modeling Framework (SWMF).

The orientation of the local interstellar magnetic field introduces asymmetries in the heliosphere that affect the location of heliospheric radio emissions and the streaming direction of ions from the termination shock of the solar wind. We combine observations of radio emissions and energetic particle streaming with extensive 3D MHD computer simulations of magnetic field draping over the heliopause to show that the plane of the local interstellar field is ~ 60-90 degrees from the galactic plane. This suggests that the field orientation in the Local Interstellar Cloud differs from that of a larger scale interstellar magnetic field thought to parallel the galactic plane.

Opher et al. 2006 showed that an interstellar magnetic field parallel to the plane defined by the deflection of interstellar hydrogen atoms can produce a north/south asymmetry in the distortion of the solar wind termination shock. This distortion is consistent with Voyager 1 and Voyager 2 observations of the direction of field-aligned streaming of the termination shock particles upstream the shock. The model also indicates that such a distortion will result in a significant north/south asymmetry in the distance to the shock and the thickness of heliosheath. The two Voyager spacecraft should reveal the nature and degree of the asymmetry in the termination shock and heliosheath.

We show that that an interstellar magnetic field can produce a north/south asymmetry in solar wind termination shock. Using Voyager 1 and 2 measurements, we suggest that the angle $\alpha$ between the interstellar wind velocity and magnetic field is $30^{\circ} < \alpha < 60^{\circ}$. The distortion of the shock is such that termination shock particles could stream outward along the spiral interplanetary magnetic field connecting Voyager 1 to the shock when the spacecraft was within $\sim 2~AU$ of the shock. The shock distortion is larger in the southern hemisphere, and Voyager 2 could be connected to the shock when it is within $\sim 5~AU$ of the shock, but with particles from the shock streaming inward along the field. Tighter constraints on the interstellar magnetic field should be possible when Voyager 2 crosses the shock in the next several years.

There is currently a controversy as to whether Voyager 1 has already crossed the Termination Shock, the first boundary of the Heliosphere. The region between the Termination Shock and the Heliopause, the Helisheath, is one of the most unknown regions theoretically. In the Heliosheath magnetic effects are crucial, as the solar magnetic field is compressed at the Termination Shock by the slowing flow. Recently, our simulations showed that the Heliosheath presents remarkable dynamics, with turbulent flows and the presence of a jet flow at the current sheet that is unstable due to magnetohydrodynamic instabilities \cite{opher,opher1}. In this paper we review these recent results, and present an additional simulation with constant neutral atom background. In this case the jet is still present but with reduced intensity. Further study, e.g., including neutrals and the tilt of the solar rotation from the magnetic axis, is required before we can definitively address how the Heliosheath behaves. Already we can say that this region presents remarkable dynamics, with turbulent flows, indicating that the Heliosheath might be very different from what we previously thought.

To model the interaction between the solar wind and the interstellar wind, magnetic fields must be included. Recently Opher et al. 2003 found that, by including the solar magnetic field in a 3D high resolution simulation using the University of Michigan BATS-R-US code, a jet-sheet structure forms beyond the solar wind Termination Shock. Here we present an even higher resolution three-dimensional case where the jet extends for $150AU$ beyond the Termination Shock. We discuss the formation of the jet due to a de Laval nozzle effect and it's su bsequent large period oscillation due to magnetohydrodynamic instabilities. To verify the source of the instability, we also perform a simplified two dimensional-geometry magnetohydrodynamic calculation of a plane fluid jet embedded in a neutral sheet with the profiles taken from our 3D simulation. We find remarkable agreement with the full three-dimensional evolution. We compare both simulations and the temporal evolution of the jet showing that the sinuous mode is the dominant mode that develops into a velocity-shear-instability with a growth rate of $5 \times 10^{-9} sec^{-1}=0.027 years^{-1}$. As a result, the outer edge of the heliosphere presents remarkable dynamics, such as turbulent flows caused by the motion of the jet. Further study, e.g., including neutrals and the tilt of the solar rotation from the magnetic axis, is required before we can definitively address how this outer boundary behaves. Already, however, we can say that the magnetic field effects are a major player in this region changing our previous notion of how the solar system ends.

The Voyager spacecraft is now approaching the edge of the solar system. Near the boundary between the solar system and the interstellar medium we find that an unstable ``jet-sheet'' forms. The jet-sheet oscillates up and down due to a velocity shear instability. This result is due to a novel application of a state-of-art 3D Magnetohydrodynamic (MHD) code with a highly refined grid. We assume as a first approximation that the solar magnetic and rotation axes are aligned. The effect of a tilt of the magnetic axis with respect to the rotation axis remains to be seen. We include in the model self-consistently magnetic field effects in the interaction between the solar and interstellar winds. Previous studies of this interaction had poorer spatial resolution and did not include the solar magnetic field. This instability can affect the entry of energetic particles into the solar system and the intermixing of solar and interstellar material. The same effect found here is predicted for the interaction of rotating magnetized stars possessing supersonic winds and moving with respect to the interstellar medium, such as O stars.

The effects of the highly damped modes in the energy and reaction rates in a plasma are discussed. These modes, with wavenumbers $k \gg k_{D}$, even being only weakly excited, with less than $k_{B}T$ per mode, make a significant contribution to the energy and screening in a plasma. When the de Broglie wavelength is much less than the distance of closest approach of thermal electrons, a classical analysis of the plasma can

The fluctuation-dissipation theorem is used to evaluate the screening factor of nuclear reactions due to the electromagnetic fluctuations in a plasma. We show that the commonly used Saltpeter factor is obtained if only fluctuations near the plasma eigenfrequency are assumed to be important (\omega \sim \omega_{pe}\ll T (\hbar=k_{B}=1)). By taking into account all the fluctuations, the highly damped ones, with \omega >\omega_{pe}, as well as those with \omega\leq\omega_{pe}, we find that nuclear reaction rates are higher than those obtained using the Saltpeter factor, for many interesting plasmas.

It has recently been argued that there are no dynamic screening corrections to Salpeter's enhancement factor in thermonuclear reactions, in the weak-screening limit. The arguments used were: 1) The Gibbs probability distribution is factorable into two parts, one of which, $exp(-\beta \sum e_{i}e_{j}/r_{ij})$ ($\beta=1/k_{B}T$), is independent of velocity space; and 2) The enhancement factor is $w=1+\beta^{2}e^{2}Z_{1}Z_{2}< \phi^{2} >$ with ${< \phi^{2} >}_{k}={< E^{2} >}_{k}/k^{2}$ and ${< E^{2} >}_{k}/(8\pi)=(T/2)[1-\epsilon^{-1} (0,k)]$. We show that both of these arguments are incorrect.

When \lambda_{T} << d_{T}, where \lambda_{T} is the de Broglie wavelength and d_{T}, the distance of closest approach of thermal electrons, a classical analysis of the energy of a plasma can be made. In all the classical analysis made until now, it was assumed that the frequency of the fluctuations \omega << T (k_{B}=\hbar=1). Using the fluctuation-dissipation theorem, we evaluate the energy of a plasma, allowing the frequency of the fluctuations to be arbitrary. We find that the energy density is appreciably larger than previously thought for many interesting plasmas, such as the plasma of the Universe before the recombination era.

The origin of the magnetic field in galaxies is an open question in astrophysics. Several mechanisms have been proposed related, in general, with the generation of small seed fields amplified by a dynamo mechanism. In general, these mechanisms have difficulty in satisfying both the requirements of a sufficiently high strength for the magnetic field and the necessary large coherent scales. We show that the formation of dense and turbulent shells of matter, in the multiple explosion scenario of Miranda and Opher (1996, 1997) for the formation of the large-scale structures of the Universe, can naturally act as a seed for the generation of a magnetic field. During the collapse and explosion of Population III objects, a temperature gradient not parallel to a density gradient can naturally be established, producing a seed magnetic field through the Biermann battery mechanism. We show that seed magnetic fields $\sim 10^{-12}-10^{-14}G$ can be produced in this multiple explosion scenario on scales of the order of clusters of galaxies (with coherence length $L\sim 1.8Mpc$) and up to $\sim 4.5\times 10^{-10}G$ on scales of galaxies ($L\sim 100 kpc$).

It is assumed, in general, that the electromagnetic spectrum in the Primordial Universe was a blackbody spectrum in vacuum. We derive the electromagnetic spectrum, based on the Fluctuation-Dissipation Theorem that describes the electromagnetic fluctuations in a plasma. Our description includes thermal and collisional effects in a plasma. The electromagnetic spectrum obtained differs from the blackbody spectrum in vacuum at low frequencies. In particular, concentrating on the primordial nucleosynthesis era, it has more energy for frequencies less than 3 to $6\omega_{pe}$, where $\omega_{pe}$ is the electron plasma frequency.

The low frequency magnetic field spectrum in the primordial plasma is of particular interest as a possible origin of magnetic fields in the universe (e.g., Tajima et al. 1992 and Cable and Tajima 1992). We derive the magnetic field spectrum in the primordial plasma, in particular, at the epoch of primordial nucleosynthesis. The pioneering study of Cable and Tajima (1992) of the electromagnetic fluctuations, based on the Fluctuation-Dissipation Theorem, is extended. Our model describes both the thermal and collisional effects in a plasma. It is based on a kinetic description with the BGK collision term. It is shown that the zero-frequency peak found by Cable and Tajima (1992) decreases. At high frequencies, the blackbody spectrum is obtained naturally without the necessity of the link procedure used by them. At low frequencies ($\omega \leq 4\omega_{pe}$, where $\omega_{pe}$ is the electron plasma frequency) it is shown that the magnetic field spectrum has more energy than the blackbody spectrum in vacuum.