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

A shock tube problem is solved numerically by using one-dimensional full particle-in-cell simulations under the condition that a relatively tenuous and weakly magnetized plasma is continuously pushed by a relatively dense and strongly magnetized plasma having supersonic relative velocity. A forward and a reverse shock and a contact discontinuity are self-consistently reproduced. The spatial width of the contact discontinuity increases as the angle between the discontinuity normal and ambient magnetic field decreases. The inner structure of the discontinuity shows different profiles between magnetic field and plasma density, or pressure, which is caused by a non-MHD effect of the local plasma. The region between the two shocks is turbulent. The fluctuations in the relatively dense plasma are compressible and propagating away from the contact discontinuity, although the fluctuations in the relatively tenuous plasma contain both compressible and incompressible components. The source of the compressible fluctuations in the relatively dense plasma is in the relatively tenuous plasma. Only compressible fast mode fluctuations generated in the relatively tenuous plasma are transmitted through the contact discontinuity and propagate in the relatively dense plasma. These fast mode fluctuations are steepened when passing the contact discontinuity. This wave steepening and probably other effects may cause the broadening of the wave spectrum in the very local interstellar medium plasma. The results are discussed in the context of the heliospheric boundary region or heliopause.

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.

The Parker Solar Probe (PSP) observed an interplanetary coronal mass ejection (ICME) event during its first orbit around the sun, among many other events. This event is analyzed by applying a wavelet analysis technique to obtain the reduced magnetic helicity, cross helicity, and residual energy, the first two of which are magnetohydrodynamics (MHD) invariants. Our results show that the ICME, as a large scale magnetic flux rope, possesses high magnetic helicity, very low cross helicity, and highly negative residual energy, thus pointing to a magnetic fluctuation dominated structure. Using the same technique, we also search for small-scale coherent magnetic flux rope structures during the period from 2018/10/22--2018/11/21, which are intrinsic to quasi-2D MHD turbulence in the solar wind. Multiple structures with duration between 8 and 300 minutes are identified from PSP in-situ spacecraft measurements. The location and scales of these structures are characterized by wavelet spectrograms of the normalized reduced magnetic helicity, normalized cross helicity and normalized residual energy. Transport theory suggests that these small-scale magnetic flux ropes may contribute to the acceleration of charged particles through magnetic reconnection processes, and the dissipation of these structures may be important for understanding the coronal heating processes.

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

The basic physics and recent progresses in theoretical and particle-in-cell (PIC) simulation studies of particle acceleration in multi-island magnetic reconnection are briefly reviewed. Particle acceleration in multi-island magnetic reconnection is considered a plausible mechanism for the acceleration of energetic particles in solar flares and the solar wind. Theoretical studies have demonstrated that such a mechanism can produce the observed power-law energy distribution of energetic particles if the particle motion is sufficiently randomized in the reconnection event. However, PIC simulations seem to suggest that the first-order Fermi acceleration mechanism is unable to produce a power-law particle energy distribution function in mildly relativistic multi-island magnetic reconnections. On the other hand, while simulations of highly relativistic reconnections appear to be able to produce a power-law energy spectrum, the spectral indices obtained are generally harder than the soft power-law spectra with indices $\sim -5$ commonly observed in the solar wind and solar flare events. In addition, the plasma heating due to kinetic instabilities in 3D magnetic reconnection may "thermalize" the power-law particles, making it even more difficult for multi-island reconnections to generate a power-law spectrum. We discuss the possible reasons that may lead to these problems.

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

The current progress in the detection of terrestrial type exoplanets has opened a new avenue in the characterization of exoplanetary atmospheres and in the search for biosignatures of life with the upcoming ground-based and space missions. To specify the conditions favorable for the origin, development and sustainment of life as we know it in other worlds, we need to understand the nature of astrospheric, atmospheric and surface environments of exoplanets in habitable zones around G-K-M dwarfs including our young Sun. Global environment is formed by propagated disturbances from the planet-hosting stars in the form of stellar flares, coronal mass ejections, energetic particles, and winds collectively known as astrospheric space weather. Its characterization will help in understanding how an exoplanetary ecosystem interacts with its host star, as well as in the specification of the physical, chemical and biochemical conditions that can create favorable and/or detrimental conditions for planetary climate and habitability along with evolution of planetary internal dynamics over geological timescales. A key linkage of (astro) physical, chemical, and geological processes can only be understood in the framework of interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary and Earth sciences. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets will significantly expand the current definition of the habitable zone to the biogenic zone and provide new observational strategies for searching for signatures of life. The major goal of this paper is to describe and discuss the current status and recent progress in this interdisciplinary field and to provide a new roadmap for the future development of the emerging field of exoplanetary science and astrobiology.

The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditions favorable for the origin, development and sustainment of life as we know it. This requires the understanding of global (astrospheric) and local (atmospheric, surface and internal) environments of exoplanets in the framework of the physical processes of the interaction between evolving planet-hosting stars along with exoplanetary evolution over geological timescales, and the resulting impact on climate and habitability of exoplanets. Feedbacks between astrophysical, physico-chemical atmospheric and geological processes can only be understood through interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary, Earth sciences, astrobiology, and the origin of life communities. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets and potential exomoons around them may significantly modify the extent and the location of the habitable zone and provide new directions for searching for signatures of life. Thus, characterization of stellar ionizing outputs becomes an important task for further understanding the extent of habitability in the universe. The goal of this white paper is to identify and describe promising key research goals to aid the theoretical characterization and observational detection of ionizing radiation from quiescent and flaring upper atmospheres of planet hosts as well as properties of stellar coronal mass ejections and stellar energetic particle events.

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

It is currently unknown how common life is on exoplanets, or how long planets can remain viable for life. To date, we have a superficial notion of habitability, a necessary first step, but so far lacking an understanding of the detailed interaction between stars and planets over geological timescales, dynamical evolution of planetary systems, and atmospheric evolution on planets in other systems. A planet mass, net insolation, and atmospheric composition alone are insufficient to determine the probability that life on a planet could arise or be detected. The latter set of planetary considerations, among others, underpin the concept of the habitable zone (HZ), defined as the circumstellar region where standing bodies of liquid water could be supported on the surface of a rocky planet. However, stars within the same spectral class are often treated in the same way in HZ studies, without any regard for variations in activity among individual stars. Such formulations ignore differences in how nonthermal emission and magnetic energy of transient events in different stars affect the ability of an exoplanet to retain its atmosphere.In the last few years there has been a growing appreciation that the atmospheric chemistry, and even retention of an atmosphere in many cases, depends critically on the high-energy radiation and particle environments around these stars. Indeed, recent studies have shown stellar activity and the extreme space weather, such as that created by the frequent flares and coronal mass ejections (CMEs) from the active stars and young Sun, may have profoundly affected the chemistry and climate and thus habitability of the early Earth and terrestrial type exoplanets. The goal of this white paper is to identify and describe promising key research goals to aid the field of the exoplanetary habitability for the next 20 years.

The search of life in the Universe is a fundamental problem of astrobiology and a major priority for NASA. A key area of major progress since the NASA Astrobiology Strategy 2015 (NAS15) has been a shift from the exoplanet discovery phase to a phase of characterization and modeling of the physics and chemistry of exoplanetary atmospheres, and the development of observational strategies for the search for life in the Universe by combining expertise from four NASA science disciplines including heliophysics, astrophysics, planetary science and Earth science. The NASA Nexus for Exoplanetary System Science (NExSS) has provided an efficient environment for such interdisciplinary studies. Solar flares, coronal mass ejections and solar energetic particles produce disturbances in interplanetary space collectively referred to as space weather, which interacts with the Earth upper atmosphere and causes dramatic impact on space and ground-based technological systems. Exoplanets within close in habitable zones around M dwarfs and other active stars are exposed to extreme ionizing radiation fluxes, thus making exoplanetary space weather (ESW) effects a crucial factor of habitability. In this paper, we describe the recent developments and provide recommendations in this interdisciplinary effort with the focus on the impacts of ESW on habitability, and the prospects for future progress in searching for signs of life in the Universe as the outcome of the NExSS workshop held in Nov 29 - Dec 2, 2016, New Orleans, LA. This is one of five Life Beyond the Solar System white papers submitted by NExSS to the National Academy of Sciences in support of the Astrobiology Science Strategy for the Search for Life in the Universe.

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.

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.

An essential parameter for models of coronal heating and fast solar wind acceleration that rely on the dissipation of MHD turbulence is the characteristic energy-containing length $\lambda_{\bot}$ of the squared velocity and magnetic field fluctuations ($u^2$ and $b^2$) transverse to the mean magnetic field inside a coronal hole (CH) at the base of the corona. The characteristic length scale defines directly the heating rate. We use a time series analysis of solar granulation and magnetic field measurements inside two CHs obtained with the New Solar Telescope (NST) at Big Bear Solar Observatory. A data set for transverse magnetic fields obtained with the Solar Optical Telescope/Spectro-Polarimeter (SOT/SP) aboard {\it Hinode} spacecraft was utilized to analyze the squared transverse magnetic field fluctuations $b_t^2$. Local correlation tracking (LCT) was applied to derive the squared transverse velocity fluctuations $u^2$. We find that for $u^2$-structures, Batchelor integral scale $\lambda$ varies in a range of 1800 - 2100 km, whereas the correlation length $\varsigma$ and the $e$-folding length $L$ vary between 660 and 1460 km. Structures for $b_t^2$ yield $\lambda \approx 1600$ km, $\varsigma \approx 640$ km, and $L \approx 620$ km. An averaged (over $\lambda, \varsigma$, and $L$) value of the characteristic length of $u^2$-fluctuations is 1260$\pm$500 km, and that of $b_t^2$ is 950$\pm$560 km. The characteristic length scale in the photosphere is approximately 1.5-50 times smaller than that adopted in previous models (3-30$\times10^3$ km). Our results provide a critical input parameter for current models of coronal heating and should yield an improved understanding of fast solar wind acceleration.

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

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.

Because of the dynamic nature of the interstellar medium, the Sun should have encountered a variety of different interstellar environments in its lifetime. As the solar wind interacts with the surrounding interstellar medium to form a heliosphere, different heliosphere shapes, sizes, and particle contents result from the different environments. Some of the large possible interstellar parameter space (density, velocity, temperature) is explored here with the help of global heliosphere models, and the features in the resulting heliospheres are compared and discussed. The heliospheric size, expressed as distance of the nose of the heliopause to the Sun, is set by the solar wind - interstellar pressure balance, even for extreme cases. Other heliospheric boundary locations and neutral particle results correlate with the interstellar parameters as well. If the H0 clouds identified in the Millennium Arecibo survey are typical of clouds encountered by the Sun, then the Sun spends ~99.4% of the time in warm low density ISM, where the typical upwind heliosphere radii are up to two orders of magnitude larger than at present.

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.

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

Measurements of stellar mass loss rates are used to assess how wind strength varies with coronal activity and age for solar-like stars. Mass loss generally increases with activity, but we find evidence that winds suddenly weaken at a certain activity threshold. Very active stars are often observed to have polar starspots, and we speculate that the magnetic field geometry associated with these spots may be inhibiting the winds. Our inferred mass-loss/age relation represents an empirical estimate of the history of the solar wind. This result is important for planetary studies as well as solar/stellar astronomy, since solar wind erosion may have played an important role in the evolution of planetary atmospheres.

We search the Hubble Space Telescope (HST) archive for previously unanalyzed observations of stellar H I Lyman-alpha emission lines, our primary purpose being to look for new detections of Lyman-alpha absorption from the outer heliosphere, and to also search for analogous absorption from the astrospheres surrounding the observed stars. The astrospheric absorption is of particular interest because it can be used to study solar-like stellar winds that are otherwise undetectable. We find and analyze 33 HST Lyman-alpha spectra in the archive. All the spectra were taken with the E140M grating of the Space Telescope Imaging Spectrograph (STIS) instrument on board HST. The HST/STIS spectra yield 4 new detections of heliospheric absorption (70 Oph, Xi Boo, 61 Vir, and HD 165185) and 7 new detections of astrospheric absorption (EV Lac, 70 Oph, Xi Boo, 61 Vir, Delta Eri, HD 128987, and DK UMa), doubling the previous number of heliospheric and astrospheric detections. When combined with previous results, 10 of 17 lines of sight within 10 pc yield detections of astrospheric absorption. This high detection fraction implies that most of the ISM within 10 pc must be at least partially neutral, since the presence of H I within the ISM surrounding the observed star is necessary for an astrospheric detection. In contrast, the detection percentage is only 9.7% (3 out of 31) for stars beyond 10 pc. Our Lyman-alpha analyses provide measurements of ISM H I and D I column densities for all 33 lines of sight, and we discuss some implications of these results. Finally, we measure chromospheric Lyman-alpha fluxes from the observed stars. We use these fluxes to determine how Lyman-alpha flux correlates with coronal X-ray and chromospheric Mg II emission, and we also study how Lyman-alpha emission depends on stellar rotation.

Interstellar material surrounding an extrasolar planetary system interacts with the stellar wind to form the stellar astrosphere, and regulates the properties of the interplanetary medium and cosmic ray fluxes throughout the system. Advanced life and civilization developed on Earth during the time interval when the Sun was immersed in the vacuum of the Local Bubble and the heliosphere was large, and probably devoid of most anomalous and galactic cosmic rays. The Sun entered an outflow of diffuse cloud material from the Sco-Cen Association within the past several thousand years. By analogy with the Sun and solar system, the Galactic environment of an extrasolar planetary system must be a key component in understanding the distribution of systems with stable interplanetary environments, and inner planets which are shielded by stellar winds from interstellar matter (ISM), such as might be expected for stable planetary climates.

Ultraviolet spectra from the GHRS instrument on board the Hubble Space Telescope reveal the presence of a mysterious absorption feature in the Mg II h & k lines of the nearby (d=20.0 pc) K5 III star Alpha Tau. The narrow absorption looks like an interstellar absorption feature but it is in the wrong location based on our knowledge of the local ISM flow vector. Since the absorption is close to the rest frame of the star, it has been interpreted as being from the interaction region between Alpha Tau's massive, cool wind and the interstellar medium, i.e., Alpha Tau's "astrosphere". We compute hydrodynamic models of the Alpha Tau astrosphere in order to see if the models can reproduce the Mg II absorption feature. These models do predict that stellar wind material heated, decelerated, and compressed after passing through a termination shock a few thousand AU from the star should produce a Mg II absorption feature with about the right width at roughly the right velocity. However, our first models underestimate the Mg II column density by an order of magnitude. A much larger parameter search is necessary to see whether the observed Mg II absorption can be reproduced by acceptable changes to the adopted stellar wind and ISM properties.

The heating of the lower solar corona is examined using numerical simulations and theoretical models of magnetohydrodynamic turbulence in open magnetic regions. A turbulent energy cascade to small length scales perpendicular to the mean magnetic field can be sustained by driving with low-frequency Alfven waves reflected from mean density and magnetic field gradients. This mechanism deposits energy efficiently in the lower corona, and we show that the spatial distribution of the heating is determined by the mean density through the Alfven speed profile. This provides a robust heating mechanism that can explain observed high coronal temperatures and accounts for the significant heating (per unit volume) distribution below two solar radius needed in models of the origin of the solar wind. The obtained heating per unit mass on the other hand is much more extended indicating that the heating on a per particle basis persists throughout all the lower coronal region considered here.

Collisions between the winds of solar-like stars and the local ISM result in a population of hot hydrogen gas surrounding these stars. Absorption from this hot H I can be detected in high resolution Lyman-alpha spectra of these stars from the Hubble Space Telescope. The amount of absorption can be used as a diagnostic for the stellar mass loss rate. We present new mass loss rate measurements derived in this fashion for four stars (Epsilon Eri, 61 Cyg A, 36 Oph AB, and 40 Eri A). Combining these measurements with others, we study how mass loss varies with stellar activity. We find that for the solar-like GK dwarfs, the mass loss per unit surface area is correlated with X-ray surface flux. Fitting a power law to this relation yields Mdot ~ Fx^(1.15+/-0.20). The active M dwarf Proxima Cen and the very active RS CVn system Lambda And appear to be inconsistent with this relation. Since activity is known to decrease with age, the above power law relation for solar-like stars suggests that mass loss decreases with time. We infer a power law relation of Mdot ~ t^(-2.00+/-0.52). This suggests that the solar wind may have been as much as 1000 times more massive in the distant past, which may have had important ramifications for the history of planetary atmospheres in our solar system, that of Mars in particular.

We study H I Lyman-alpha absorption observed by the Hubble Space Telescope toward the nearby binary system Alpha Cen (G2 V+K0 V) and its distant companion star Proxima Cen (M5.5 Ve). Absorption from heliospheric H I heated by the solar wind/ISM interaction is observed toward both Alpha Cen and Proxima Cen. Absorption from analogous "astrospheric" material surrounding the stars is detected toward Alpha Cen, but not Proxima Cen. The nondetection of astrospheric absorption toward Proxima Cen suggests that the stellar wind of Proxima Cen must be significantly weaker than that of the Alpha Cen system. We use hydrodynamic models of the astrospheres computed assuming different mass-loss rates to predict astrospheric Lyman-alpha absorption for comparison with the observations. The model that best matches the Alpha Cen data has a mass-loss rate of twice the solar rate, and the models suggest an upper limit of 0.2 solar for Proxima Cen. Finally, we note that the heliospheric absorption observed toward Proxima Cen in 2000 May is identical to the heliospheric absorption observed toward Alpha Cen in 1995 May, implying that the structure of the outer heliosphere does not change significantly during the solar activity cycle.

We use self-consistent kinetic/hydrodynamic models of the heliosphere to predict H I Lyman-alpha absorption profiles for various lines of sight through the heliosphere. These results are compared with Lyman-alpha absorption lines of six nearby stars observed by the Hubble Space Telescope. The directions of these lines of sight range from nearly upwind (36 Oph) to nearly downwind (Epsilon Eri). Only three of the Lyman-alpha spectra (36 Oph, Alpha Cen, and Sirius) actually show evidence for the presence of heliospheric absorption, which is blended with the ubiquitous interstellar absorption, but the other three spectra still provide useful upper limits for the amount of heliospheric absorption for those lines of sight. Most of our models use a Boltzmann particle code for the neutrals, allowing us to estimate neutral velocity distributions throughout the heliosphere, from which we compute model Lyman-alpha absorption profiles. In comparing these models with the data, we find they predict too much absorption in sidewind and downwind directions, especially when higher Mach numbers are assumed for the interstellar wind. Models created assuming different values of the interstellar temperature and proton density fail to improve the agreement. Somewhat surprisingly, a model that uses a multi-fluid treatment of the neutrals rather than the Boltzmann particle code is more consistent with the data, and we speculate as to why this may be the case.

We use high-resolution UV spectra taken by the Space Telescope Imaging Spectrograph instrument on board the Hubble Space Telescope to study the 5.5 pc line of sight to the K0 V star 36 Oph A. The one interstellar component detected for this line of sight has a velocity inconsistent with the local interstellar cloud (LIC) flow vector, but consistent with the flow vector of the so-called G cloud, a very nearby warm cloud in the Galactic Center direction. From our data, we measure the following values for the interstellar temperature, nonthermal velocity, H I column density, and D/H value: T=5900+/-500 K, xi=2.2+/-0.2$ km/s, log N(H)=17.85+/-0.15, and D/H=(1.5+/-0.5)e-5. The temperature of the G cloud is somewhat lower than that of the LIC, and Mg and Fe depletions are also lower, but the D/H value appears to be the same. Based on upper limits for the LIC absorption, we estimate the distance to the edge of the LIC to be d(edge)<0.19 pc, which the Sun will reach in t(edge)<7400 yrs. The H I Lyman-alpha absorption line has properties inconsistent with those of the other absorption lines, indicating the presence of one or more absorption components not seen in the other lines. We present evidence that hot neutral hydrogen local to both the Sun and the star is responsible for the excess Lyman-alpha absorption. This hot H I is created by the interaction between the ISM and the winds of the Sun and 36 Oph A. The observed line of sight lies only 12 degrees from the upwind direction of the LIC flow vector, where hydrodynamic models of the heliosphere suggest that heliospheric H I absorption should be particularly prominent.

The gas-to-dust mass ratios found for interstellar dust within the Solar System, versus values determined astronomically for the cloud around the Solar System, suggest that large and small interstellar grains have separate histories, and that large interstellar grains preferentially detected by spacecraft are not formed exclusively by mass exchange with nearby interstellar gas. Observations by the Ulysses and Galileo satellites of the mass spectrum and flux rate of interstellar dust within the heliosphere are combined with information about the density, composition, and relative flow speed and direction of interstellar gas in the cloud surrounding the solar system to derive an in situ value for the gas-to-dust mass ratio, $R_{g/d} = 94^{+46}_{-38}$. Hubble observations of the cloud surrounding the solar system yield a gas-to-dust mass ratio of Rg/d=551+61-251 when B-star reference abundances are assumed. The exclusion of small dust grains from the heliosheath and heliosphere regions are modeled, increasing the discrepancy between interstellar and in situ observations. The shock destruction of interstellar grains is considered, and comparisons are made with interplanetary and presolar dust grains.

The interaction of the heliosphere with interstellar clouds has attracted interest since the late 1920's, both with a view to explaining apparent quasi-periodic climate "catastrophes" as well as periodic mass extinctions. Until recently, however, models describing the solar wind - local interstellar medium (LISM) interaction self-consistently had not been developed. Here, we describe the results of a two-dimensional (2D) simulation of the interaction between the heliosphere and an interstellar cloud with the same properties as currently, except that the neutral H density is increased from the present value of n(H) ~ 0.2 cm^-3 to 10 cm^-3. The mutual interaction of interstellar neutral hydrogen and plasma is included. The heliospheric cavity is reduced considerably in size (approximately 10 - 14 AU to the termination shock in the upstream direction) and is highly dynamical. The interplanetary environment at the orbit of the Earth changes markedly, with the density of interstellar H increasing to ~2 cm^-3. The termination shock itself experiences periods where it disappears, reforms and disappears again. Considerable mixing of the shocked solar wind and LISM occurs due to Rayleigh-Taylor-like instabilities at the nose, driven by ion-neutral friction. Implications for two anomalously high concentrations of 10Be found in Antarctic ice cores 33 kya and 60 kya, and the absence of prior similar events, are discussed in terms of density enhancements in the surrounding interstellar cloud. The calculation presented here supports past speculation that the galactic environment of the Sun moderates the interplanetary environment at the orbit of the Earth, and possibly also the terrestrial climate.

Redshifted L$\alpha$ absorption toward $\alpha$ Cen has been interpreted by Linsky & Wood (1996) and Frisch et al. (1996) as evidence for decelerated interstellar hydrogen piled up on the upstream side of the heliosphere. We utilize newly developed two-dimensional multi-fluid models of the solar wind interaction with the ISM to corroborate this interpretation by synthesizing the L$\alpha$ absorption profile predicted for this ``hydrogen wall''. Both subsonic and supersonic inflow into the heliosphere are considered, corresponding to one-shock and two-shock global morphologies, respectively. It is found that these two extremes give observably different redward absorption characteristics in the L$\alpha$ profiles, and our preliminary conclusion is that the L$\alpha$ profiles seen toward $\alpha$ Cen favor a barely subsonic model (Mach number 0.9). For such a model to hold, additional interstellar pressure terms, such as cosmic ray or magnetic pressures, must contribute. To make this conclusion more certain, an extended model-parameter survey is required, coupled with L$\alpha$ data along other lines of sight.