In 2004, the Rosetta spacecraft was sent to comet 67P/Churyumov-Gerasimenko for the first ever long-term investigation of a comet. After its arrival in 2014, the spacecraft spent more than two years in immediate proximity to the comet. During these two years, the ROSINA Double Focusing Mass Spectrometer (DFMS) onboard Rosetta discovered a coma with an unexpectedly complex chemical composition that included many oxygenated molecules. Determining the exact cometary composition is an essential first step to understanding of the organic rich chemistry in star forming regions and protoplanetary disks that are ultimately conserved in cometary ices. In this study a joint approach of laboratory calibration and space data analysis was used to perform a detailed identification and quantification of CHO-compounds in the coma of 67P/Churyumov-Gerasimenko. The goal was to derive the CHO-compound abundances relative to water for masses up to 100 u. For this study, the May 2015 post-equinox period represent the best bulk abundances of comet 67P/Churyumov-Gerasimenko. A wide variety of CHO-compounds were discovered and their bulk abundances were derived. Finally, these results are compared to abundances of CHO-bearing molecules in other comets, obtained mostly from ground-based observations and modelling.

We develop a mixed Long Short Term Memory (LSTM) regression model to predict the maximum solar flare intensity within a 24-hour time window 0$\sim$24, 6$\sim$30, 12$\sim$36 and 24$\sim$48 hours ahead of time using 6, 12, 24 and 48 hours of data (predictors) for each Helioseismic and Magnetic Imager (HMI) Active Region Patch (HARP). The model makes use of (1) the Space-weather HMI Active Region Patch (SHARP) parameters as predictors and (2) the exact flare intensities instead of class labels recorded in the Geostationary Operational Environmental Satellites (GOES) data set, which serves as the source of the response variables. Compared to solar flare classification, the model offers us more detailed information about the exact maximum flux level, i.e. intensity, for each occurrence of a flare. We also consider classification models built on top of the regression model and obtain better results in solar flare classifications. Our results suggest that the most efficient time period for predicting the solar activity is within 24 hours before the prediction time using the SHARP parameters and the LSTM model.

A deep learning network, Long-Short Term Memory (LSTM) network, is used in this work to predict whether the maximum flare class an active region (AR) will produce in the next 24 hours is class $\Gamma$. We considered $\Gamma$ are $\ge M$, $\ge C$ and any flare class. The essence of using LSTM, which is a recurrent neural network, is its capability to capture temporal information of the data samples. The input features are time sequences of 20 magnetic parameters from SHARPs - Space-weather HMI Active Region Patches. We analyzed active regions from June 2010 to Dec 2018, using the Geostationary Operational Environmental Satellite (GOES) X-ray flare catalogs and label the data samples with identified ARs in the GOES X-ray flare catalogs. Our results (i) shows consistent skill scores with recently published results using LSTMs and better than the previous work using single time input (eg. DeFN) (ii) The skill scores from the model show essential differences when different years of data was chosen for training and testing.

Cometary comae are generally depleted in nitrogen. The main carriers for volatile nitrogen in comets are NH3 and HCN. It is known that ammonia readily combines with many acids like e.g. HCN, HNCO, HCOOH, etc. encountered in the interstellar medium as well as in cometary ice to form ammonium salts (NH4+X-) at low temperatures. Ammonium salts, which can play a significant role in prebiotic chemistry, are hard to detect in space as they are unstable in the gas phase and their infrared signature is often hidden by thermal radiation or by e.g. OH in minerals. Here we report the presence of all possible sublimation products of five different ammonium salts at comet 67P/Churyumov-Gerasimenko measured by the ROSINA instrument on Rosetta. The relatively high sublimation temperatures of the salts leads to an apparent lack of volatile nitrogen in the coma. This then also explains the observed trend of higher NH3/H2O ratios with decreasing perihelion distances in comets.

Heliospheric plasmas require multi-scale and multi-physics considerations. On one hand, MHD codes are widely used for global simulations of the solar-terrestrial environments, but do not provide the most elaborate physical description of space plasmas. Hybrid codes, on the other hand, capture important physical processes, such as electric currents and effects of finite Larmor radius, but they can be used locally only, since the limitations in available computational resources do not allow for their use throughout a global computational domain. In the present work, we present a new coupled scheme which allows to switch blocks in the block-adaptive grids from fluid MHD to hybrid simulations, without modifying the self-consistent computation of the electromagnetic fields acting on fluids (in MHD simulation) or charged ion macroparticles (in hybrid simulation). In this way, the hybrid scheme can refine the description in specified regions of interest without compromising the efficiency of the global MHD code.

Simulation of plasmas in the electromagnetic fields requires to solve numerically a kinetic equation, describing the time evolution of the particle distribution function. Here, we propose a finite volume scheme based on the integral relation for the Poisson bracket to solve the most fundamental kinetic equation, namely, the Liouville equation. The proposed scheme conserves the number of particles, maintains the total-variation-diminishing (TVD) property, and provides high-quality numerical results. Some other types of kinetic equations may be also formulated in terms of the Poisson brackets and solved with the proposed method. Among them is the focused transport equation describing the acceleration and propagation of the Solar Energetic Particles (SEPs), which is of practical importance, since the high energy SEPs produce radiation hazards. The newly proposed scheme is demonstrated to be accurate and efficient, which makes it applicable to global simulation systems analysing the space weather. We also discuss a role of focused transport and the accuracy of the diffusive approximation, in application to the SEPs

The Rosetta Orbiter Spectrometer for Ion and Neutral Analysis (ROSINA) suite of instruments operated throughout the over two years of the Rosetta mission operations in the vicinity of comet 67P/Churyumov-Gerasimenko. It measured gas densities and composition throughout the comet's atmosphere, or coma. Here we present two-years' worth of measurements of the relative densities of the four major volatile species in the coma of the comet, H2O. CO2, CO and O2, by one of the ROSINA sub-systems called the Double Focusing Mass Spectrometer (DFMS). The absolute total gas densities were provided by the Comet Pressure Sensor (COPS), another ROSINA sub-system. DFMS is a very high mass resolution and high sensitivity mass spectrometer able to resolve at a tiny fraction of an atomic mass unit. We have analyzed the combined DFMS and COPS measurements using an inversion scheme based on spherical harmonics that solves for the distribution of potential surface activity of each species as the comet rotates, changing solar illumination, over short intervals and as the comet changes distance from the sun and orientation of its spin axis over long time intervals. We also use the surface boundary conditions derived from the inversion scheme to simulate the whole coma with our fully kinetic Direct Simulation Monte Carlo model and calculate the production rates of the four major species throughout the mission. We compare the derived production rates with revised remote sensing observations by the Visible and Infrared Thermal Imaging Spectrometer (VIRTIS) as well as with published observations from the Microwave Instrument for the Rosetta Orbiter (MIRO). Finally we use the variation of the surface production of the major species to calculate the total mass loss over the mission and, for different estimates of the dust/gas ratio, calculate the variation of surface loss over the nucleus.

Comets are considered to be some of the most pristine and unprocessed solar system objects accessible to in-situ exploration. Investigating their molecular and elemental composition takes us on a journey back to the early period of our solar system and possibly even further. In this work, we deduce the bulk abundances of the major volatile species in comet 67P/Churyumov-Gerasimenko, the target of the European Space Agency's Rosetta mission. The basis are measurements obtained with the ROSINA instrument suite on board the Rosetta orbiter during a suitable period of high outgassing near perihelion. The results are combined with both gas and dust composition measurements published in the literature. This provides an integrated inventory of the major elements present in the nucleus of 67P/Churyumov-Gerasimenko. Similar to comet 1P/Halley, which was visited by ESA's Giotto spacecraft in 1986, comet 67P/Churyumov-Gerasimenko also shows near-solar abundances of oxygen and carbon, whereas hydrogen and nitrogen are depleted compared to solar. Still, the degree of devolatilization is lower than that of inner solar system objects, including meteorites and the Earth. This supports the idea that comets are among the most pristine objects in our solar system.

We present several methods towards construction of precursors, which show great promise towards early predictions, of solar flare events in this paper. A data pre-processing pipeline is built to extract useful data from multiple sources, Geostationary Operational Environmental Satellites (GOES) and Solar Dynamics Observatory (SDO)/Helioseismic and Magnetic Imager (HMI), to prepare inputs for machine learning algorithms. Two classification models are presented: classification of flares from quiet times for active regions and classification of strong versus weak flare events. We adopt deep learning algorithms to capture both the spatial and temporal information from HMI magnetogram data. Effective feature extraction and feature selection with raw magnetogram data using deep learning and statistical algorithms enable us to train classification models to achieve almost as good performance as using active region parameters provided in HMI/Space-Weather HMI-Active Region Patch (SHARP) data files. Case studies show a significant increase in the prediction score around 20 hours before strong solar flare events.

The nucleus of the Jupiter-family comet 67P/Churyumov-Gerasimenko was discovered to be bi-lobate in shape when the European Space Agency spacecraft Rosetta first approached it in July 2014. The bi-lobate structure of the cometary nucleus has led to much discussion regarding the possible manner of its formation and on how the composition of each lobe might compare with that of the other. During its two-year-long mission from 2014 to 2016, Rosetta remained in close proximity to 67P/Churyumov-Gerasimenko, studying its coma and nucleus in situ. Based on lobe-specific measurements of HDO and H2O performed with the ROSINA DFMS mass spectrometer on board Rosetta, the Deuterium-to-Hydrogen ratios in water from the two lobes could be compared. No appreciable difference was observed, suggesting that both lobes formed in the same region and are homogeneous in their Deuterium-to-Hydrogen ratios.

The European Space Agency (ESA) spacecraft Rosetta accompanied the Jupiter-family comet (JFC) 67P / Churyumov-Gerasimenko for over two years along its trajectory through the inner solar system. Between 2014 and 2016, it performed almost continuous in-situ measurements of the comet's gaseous atmosphere in close proximity to its nucleus. In this study, the $^{16}$O/$^{18}$O ratio of H$_2$O in the coma of 67P / Churyumov-Gerasimenko, as measured by the ROSINA DFMS mass spectrometer on board Rosetta, was determined from the ratio of H${_2}^{16}$O / H${_2}^{18}$O and $^{16}$OH / $^{18}$OH. The value of 445 $\pm$ 35 represents an $\sim$ 11% enrichment of $^{18}$O compared with the terrestrial ratio of 498.7 $\pm$ 0.1. This cometary value is consistent with the comet containing primordial water, in accordance with leading self-shielding models. These models predict primordial water to be between 5% to 20% enriched in heavier oxygen isotopes compared to terrestrial water. Addendum: The $^{16}$O/$^{17}$O ratio of H$_2$O in the coma of 67P / Churyumov-Gerasimenko.

The ROSINA mass spectrometer DFMS on board ESA's Rosetta spacecraft detected the major isotopes of the noble gases argon, krypton, and xenon in the coma of comet 67P/Churyumov-Gerasimenko. Earlier, it has been shown that xenon exhibits an isotopic composition distinct from anywhere else in the solar system. However, argon isotopes, within error, were shown to be consistent with solar isotope abundances. This discrepancy suggested an additional exotic component of xenon in comet 67P/Churyumov-Gerasimenko. Here we show that also krypton exhibits an isotopic composition close to solar. Furthermore, we found a depletion compared to solar of argon with respect to krypton and of krypton with respect to xenon, which is a necessity to postulate an addition of exotic xenon in the comet.

The history and present state of large-scale magnetohydrodynamic (MHD) modeling of the solar corona and the solar wind with steady or quasi-steady coronal physics is reviewed. We put the evolution of ideas leading to the recognition of the existence of an expanding solar atmosphere into historical context. The development and main features of the first generation of global corona and solar wind models are described in detail. This historical perspective is also applied to the present suite of global corona and solar wind models. We discuss the evolution of new ideas and their implementation into numerical simulation codes. We point out the scientific and computational challenges facing these models and discuss the ways various groups tried to overcome these challenges. Next, we discuss the latest, state-of-the art models and point to the expected next steps in modeling the corona and the interplanetary medium.

Magnetohydrodynamics simulations have been carried out in studying the solar wind and cometary plasma interactions for decades. Various plasma boundaries have been simulated and compared well with observations for comet 1P/Halley. The Rosetta mission, which studies comet 67P/Churyumov-Gerasimenko, challenges our understanding of the solar wind and comet interactions. The Rosetta Plasma Consortium observed regions of very weak magnetic field outside the predicted diamagnetic cavity. In this paper, we simulate the inner coma with the Hall magnetohydrodynamics equations and show that the Hall effect is important in the inner coma environment. The magnetic field topology becomes complex and magnetic reconnection occurs on the dayside when the Hall effect is taken into account. The magnetic reconnection on the dayside can generate weak magnetic filed regions outside the global diamagnetic cavity, which may explain the Rosetta Plasma Consortium observations. We conclude that the substantial change in the inner coma environment is due to the fact that the ion inertial length (or gyro radius) is not much smaller than the size of the diamagnetic cavity.

We present a three-species multi-fluid MHD model (H$^+$, H$_2$O$^+$ and e$^-$), endowed with the requisite atmospheric chemistry, that is capable of accurately quantifying the magnitude of water ion losses from exoplanets. We apply this model to a water world with Earth-like parameters orbiting a Sun-like star for three cases: (i) current normal solar wind conditions, (ii) ancient normal solar wind conditions, and (iii) one extreme "Carrington-type" space weather event. We demonstrate that the ion escape rate for (ii), with a value of 6.0$\times$10$^{26}$ s$^{-1}$, is about an order of magnitude higher than the corresponding value of 6.7$\times$10$^{25}$ s$^{-1}$ for (i). Studies of ion losses induced by space weather events, where the ion escape rates can reach $\sim$ 10$^{28}$ s$^{-1}$, are crucial for understanding how an active, early solar-type star (e.g., with frequent coronal mass ejections) could have accelerated the depletion of the exoplanet's atmosphere. We briefly explore the ramifications arising from the loss of water ions, especially for planets orbiting M-dwarfs where such effects are likely to be significant.

Context. The Rosetta Orbiter Spectrometer for Ion and Neutral Analysis (ROSINA) was designed to measure the composition of the gas in the coma of comet 67P/Churyumov-Gerasimenko, the target of the European Space Agency's Rosetta mission. In addition to the volatiles, ROSINA measured refractories sputtered off the comet by the interaction of solar wind protons with the surface of the comet. Aims. The origin of different solar system materials is still heavily debated. Isotopic ratios can be used to distinguish between different reservoirs and investigate processes occurring during the formation of the solar system. Methods. ROSINA consisted of two mass spectrometers and a pressure sensor. In the ROSINA Double Focusing Mass Spectrometer (DFMS), the neutral gas of cometary origin was ionized and then deflected in an electric and a magnetic field that separated the ions based on their mass-to-charge ratio. The DFMS had a high mass resolution, dynamic range, and sensitivity that allowed detection of rare species and the known major volatiles. Results. We measured the relative abundance of all three stable silicon isotopes with the ROSINA instrument on board the Rosetta spacecraft. Furthermore, we measured $^{13}$C/$^{12}$C in C$_2$H$_4$, C$_2$H$_5$, and CO. The DFMS in situ measurements indicate that the average silicon isotopic composition shows depletion in the heavy isotopes $^{29}$Si and $^{30}$Si with respect to $^{28}$Si and solar abundances, while $^{13}$C to $^{12}$C is analytically indistinguishable from bulk planetary and meteorite compositions. Although the origin of the deficiency of the heavy silicon isotopes cannot be explained unambiguously, we discuss mechanisms that could have contributed to the measured depletion of the isotopes $^{29}$Si and $^{30}$Si.

We present an updated global model of the solar corona, including the transition region. We simulate the realistic tree-dimensional (3D) magnetic field using the data from the photospheric magnetic field measurements and assume the magnetohydrodynamic (MHD) Alfv\'en wave turbulence and its non-linear dissipation to be the only source for heating the coronal plasma and driving the solar wind. In closed field regions the dissipation efficiency in a balanced turbulence is enhanced. In the coronal holes we account for a reflection of the outward propagating waves, which is accompanied by generation of weaker counter-propagating waves. The non-linear cascade rate degrades in strongly imbalanced turbulence, thus resulting in colder coronal holes. The distinctive feature of the presented model is the description of the low corona as almost-steady-state low-beta plasma motion and heat flux transfer along the magnetic field lines. We trace the magnetic field lines through each grid point of the lower boundary of the global corona model, chosen at some heliocentric distance, $R=R_{b}\sim1.1\ R_\odot$ well above the transition region. One can readily solve the plasma parameters along the magnetic field line from 1D equations for the plasma motion and heat transport together with the Alfv\'en wave propagation, which adequately describe physics within the heliocentric distances range, $R_{\odot}<R<R_{b}$, in the low solar corona. By interfacing this threaded-field-lines model with the full MHD global corona model at $r=R_{b}$, we find the global solution and achieve a faster-than-real-time performance of the model on $\sim200$ cores.

Anthropogenic effects on the space environment started in the late 19th century and reached their peak in the 1960s when high-altitude nuclear explosions were carried out by the USA and the Soviet Union. These explosions created artificial radiation belts near Earth that resulted in major damages to several satellites. Another, unexpected impact of the high-altitude nuclear tests was the electromagnetic pulse (EMP) that can have devastating effects over a large geographic area (as large as the continental United States). Other anthropogenic impacts on the space environment include chemical release ex- periments, high-frequency wave heating of the ionosphere and the interaction of VLF waves with the radiation belts. This paper reviews the fundamental physical process behind these phenomena and discusses the observations of their impacts.

We perform and analyze results of a global magnetohydrodyanmic (MHD) simulation of the fast coronal mass ejection (CME) that occurred on 2011 March 7. The simulation is made using the newly developed Alfv\'en Wave Solar Model (AWSoM), which describes the background solar wind starting from the upper chromosphere and extends to 24 R$_{\odot}$. Coupling AWSoM to an inner heliosphere (IH) model with the Space Weather Modeling Framework (SWMF) extends the total domain beyond the orbit of Earth. Physical processes included in the model are multi-species thermodynamics, electron heat conduction (both collisional and collisionless formulations), optically thin radiative cooling, and Alfv\'en-wave turbulence that accelerates and heats the solar wind. The Alfv\'en-wave description is physically self-consistent, including non-Wentzel-Kramers-Brillouin (WKB) reflection and physics-based apportioning of turbulent dissipative heating to both electrons and protons. Within this model, we initiate the CME by using the Gibson-Low (GL) analytical flux rope model and follow its evolution for days, in which time it propagates beyond STEREO A. A detailed comparison study is performed using remote as well as \textit{in situ} observations. Although the flux rope structure is not compared directly due to lack of relevant ejecta observation at 1 AU in this event, our results show that the new model can reproduce many of the observed features near the Sun (e.g., CME-driven extreme ultraviolet (EUV) waves, deflection of the flux rope from the coronal hole, "double-front" in the white light images) and in the heliosphere (e.g., shock propagation direction, shock properties at STEREO A).

We present a first-principles-based coronal mass ejection (CME) model suitable for both scientific and operational purposes by combining a global magnetohydrodynamics (MHD) solar wind model with a flux rope-driven CME model. Realistic CME events are simulated self-consistently with high fidelity and forecasting capability by constraining initial flux rope parameters with observational data from GONG, SOHO/LASCO, and STEREO/COR. We automate this process so that minimum manual intervention is required in specifying the CME initial state. With the newly developed data-driven Eruptive Event Generator Gibson-Low (EEGGL), we present a method to derive Gibson-Low (GL) flux rope parameters through a handful of observational quantities so that the modeled CMEs can propagate with the desired CME speeds near the Sun. A test result with CMEs launched with different Carrington rotation magnetograms are shown. Our study shows a promising result for using the first-principles-based MHD global model as a forecasting tool, which is capable of predicting the CME direction of propagation, arrival time, and ICME magnetic field at 1 AU (see companion paper by Jin et al. 2016b).

We present the results of a comprehensive numerical simulation of the environment around three exoplanet-host stars (HD 1237, HD 22049, and HD 147513). Our simulations consider one of the latest models currently used for space weather studies in the Heliosphere. Large-scale magnetic field maps, recovered with two implementations of the tomographic technique of Zeeman-Doppler imaging, serve to drive steady-state solutions in each system. This paper contains the description of the stellar wind and inner astrosphere, while the coronal structure was previously discussed in Alvarado-G\'omez et al. (2016). The analysis includes the magneto-hydrodynamical properties of the stellar wind, the associated mass and angular momentum loss rates, as well as the topology of the astrospheric current sheet in each system. A systematic comparison among the considered cases is performed, including two reference solar simulations covering activity minimum and maximum. For HD 1237, we investigate the interactions between the structure of the developed stellar wind, and a possible magnetosphere around the Jupiter-mass planet in this system. We find that the process of particle injection into the planetary atmosphere is dominated by the density distribution rather than velocity profile of the stellar wind. In this context, we predict a maximum exoplanetary radio emission of 12 mJy at 40 MHz in this system, assuming the crossing of a high-density streamer during periastron passage. Furthermore, in combination with the analysis performed in Alvarado-G\'omez et al. (2016), we obtain for the first time a fully simulated mass loss-activity relation, which is compared and discussed in the context of the relation based on astrospheric detections proposed by Wood et al. (2005a). Finally, we provide a characterisation of the 3D properties of the stellar wind of these systems, at the inner edges of their habitable zones.

We present the results of a detailed numerical simulation of the circumstellar environment around three exoplanet-hosting stars. A state-of-the-art global magnetohydrodynamic (MHD) model is considered, including Alfv\'en wave dissipation as a self-consistent coronal heating mechanism. This paper contains the description of the numerical set-up, evaluation procedure, and the simulated coronal structure of each system (HD 1237, HD 22049 and HD 147513). The simulations are driven by surface magnetic field maps, recovered with the observational technique of Zeeman Doppler Imaging (ZDI). A detailed comparison of the simulations is performed, where two different implementations of this mapping routine are used to generate the surface field distributions. Quantitative and qualitative descriptions of the coronae of these systems are presented, including synthetic high-energy emission maps in the Extreme Ultra-Violet (EUV) and Soft X-rays (SXR) ranges. Using the simulation results, we are able to recover similar trends as in previous observational studies, including the relation between the magnetic flux and the coronal X-ray emission. Furthermore, for HD 1237 we estimate the rotational modulation of the high-energy emission due to the various coronal features developed in the simulation. We obtain variations, during a single stellar rotation cycle, up to 15\% for the EUV and SXR ranges. The results presented here will be used, in a follow-up paper, to self-consistently simulate the stellar winds and inner astrospheres of these systems.

We study the interaction between the atmospheres of Venus-like, non-magnetized exoplanets orbiting an M-dwarf star, and the stellar wind using a multi-species Magnetohydrodynaic (MHD) model. We focus our investigation on the effect of enhanced stellar wind and enhanced EUV flux as the planetary distance from the star decreases. Our simulations reveal different topologies of the planetary space environment for sub- and super-Alfvenic stellar wind conditions, which could lead to dynamic energy deposition in to the atmosphere during the transition along the planetary orbit. We find that the stellar wind penetration for non-magnetized planets is very deep, up to a few hundreds of kilometers. We estimate a lower limit for the atmospheric mass-loss rate and find that it is insignificant over the lifetime of the planet. However, we predict that when accounting for atmospheric ion acceleration, a significant amount of the planetary atmosphere could be eroded over the course of a billion years.

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 higher charge states found in slow ($<$400km s$^{-1}$) solar wind streams compared to fast streams have supported the hypothesis that the slow wind originates in closed coronal loops, and released intermittently through reconnection. Here we examine whether a highly ionized slow wind can also form along steady and open magnetic field lines. We model the steady-state solar atmosphere using AWSoM, a global magnetohydrodynamic model driven by Alfv{\'e}n waves, and apply an ionization code to calculate the charge state evolution along modeled open field lines. This constitutes the first charge states calculation covering all latitudes in a realistic magnetic field. The ratios $O^{+7}/O^{+6}$ and $C^{+6}/C^{+5}$ are compared to in-situ Ulysses observations, and are found to be higher in the slow wind, as observed; however, they are under-predicted in both wind types. The modeled ion fractions of S, Si, and Fe are used to calculate line-of-sight intensities, which are compared to EIS observations above a coronal hole. The agreement is partial, and suggests that all ionization rates are under-predicted. Assuming the presence of suprathermal electrons improved the agreement with both EIS and Ulysses observations; importantly, the trend of higher ionization in the slow wind was maintained. The results suggest there can be a sub-class of slow wind that is steady and highly ionized. Further analysis shows it originates from coronal hole boundaries (CHB), where the modeled electron density and temperature are higher than inside the hole, leading to faster ionization. This property of CHBs is global, and observationally supported by EUV tomography.

We study the magnetospheric structure and the ionospheric Joule Heating of planets orbiting M-dwarf stars in the habitable zone using a set of magnetohydrodynamic (MHD) models. The stellar wind solution is used to drive a model for the planetary magnetosphere, which is coupled with a model for the planetary ionosphere. Our simulations reveal that the space environment around close-in habitable planets is extreme, and the stellar wind plasma conditions change from sub- to super-Alfvenic along the planetary orbit. As a result, the magnetospheric structure changes dramatically with a bow shock forming in the super-Alfvenic sectors, while no bow shock forms in the sub-Alfvenic sectors. The planets reside most of the time in the sub-Alfvenic sectors with poor atmospheric protection. A significant amount of Joule Heating is provided at the top of the atmosphere as a result of the planetary interaction with the stellar wind. For the steady-state solution, the heating is about 0.1-3\% of the total incoming stellar irradiation, and it is enhanced by 50\% for the time-dependent case. The significant Joule Heating obtained here should be considered in models for the atmospheres of habitable planets in terms of the thickness of the atmosphere, the top-side temperature and density, the boundary conditions for the atmospheric pressure, and particle radiation and transport.

In the present work, we test the predictions of the AWSoM model, a global extended-MHD model capable of calculating the propagation and turbulent dissipation of Alfv\'en waves in any magnetic topology, against high resolution spectra of the quiescent off-disk solar corona. Wave dissipation is the only heating mechanism assumed in this model. Combining 3D model results with the CHIANTI atomic database, we were able to create synthetic line-of-sight spectra which include the effects of emission line broadening due to both thermal and wave-related non-thermal motions. To the best of our knowledge this is the first time a global model is used to obtain synthetic non-thermal line broadening. We obtained a steady-state solution driven by a synoptic magnetogram and compared the synthetic spectra with SUMER observations of a quiescent area above the solar west limb extending between 1.04 and 1.34 solar radii at the equator. Both the predicted line widths and the total line fluxes were consistent with the observations for 5 different ions. Using the 3D solution, we were able to locate the region that contributes the most to the emission used for measuring electron properties; we found that region to be a pseudo-streamer, whose modeled electron temperature and density are consistent with the measured ones. We conclude that the turbulent dissipation assumed in the AWSoM model can simultaneously account for the observed heating rate and the non-dissipated wave energy observed in this region.

We present a new version of the Alfven Wave Solar Model (AWSoM), a global model from the upper chromosphere to the corona and the heliosphere. The coronal heating and solar wind acceleration are addressed with low-frequency Alfven wave turbulence. The injection of Alfven wave energy at the inner boundary is such that the Poynting flux is proportional to the magnetic field strength. The three-dimensional magnetic field topology is simulated using data from photospheric magnetic field measurements. This model does not impose open-closed magnetic field boundaries; those develop self-consistently. The physics includes: (1) The model employs three different temperatures, namely the isotropic electron temperature and the parallel and perpendicular ion temperatures. The firehose, mirror, and ion-cyclotron instabilities due to the developing ion temperature anisotropy are accounted for. (2) The Alfven waves are partially reflected by the Alfven speed gradient and the vorticity along the field lines. The resulting counter-propagating waves are responsible for the nonlinear turbulent cascade. The balanced turbulence due to uncorrelated waves near the apex of the closed field lines and the resulting elevated temperatures are addressed. (3) To apportion the wave dissipation to the three temperatures, we employ the results of the theories of linear wave damping and nonlinear stochastic heating. (4) We have incorporated the collisional and collisionless electron heat conduction. We compare the simulated multi-wavelength EUV images of CR2107 with the observations from STEREO/EUVI and SDO/AIA instruments. We demonstrate that the reflection due to strong magnetic fields in proximity of active regions intensifies the dissipation and observable emission sufficiently.

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.

By analogy with the solar system, it is believed that stellar winds will form bow shocks around exoplanets. For hot Jupiters the bow shock will not form directly between the planet and the star, causing an asymmetric distribution of mass around the exoplanet and hence an asymmetric transit. As the planet orbits thorough varying wind conditions, the strength and geometry of its bow shock will change, thus producing transits of varying shape. We model this process using magnetic maps of HD 189733 taken one year apart, coupled with a 3D stellar wind model, to determine the local stellar wind conditions throughout the orbital path of the planet. We predict the time-varying geometry and density of the bow shock that forms around the magnetosphere of the planet and simulate transit light curves. Depending on the nature of the stellar magnetic field, and hence its wind, we find that both the transit duration and ingress time can vary when compared to optical light curves. We conclude that consecutive near-UV transit light curves may vary significantly and can therefore provide an insight into the structure and evolution of the stellar wind.

We describe, analyze and validate the recently developed Alfv\'en Wave Solar Model (AWSoM), a 3D global model starting from the top of the chromosphere and extending into interplanetary space (up to 1-2 AU). This model solves the extended two temperature magnetohydrodynamics equations coupled to a wave kinetic equation for low frequency Alfv\'en waves. In this picture, heating and acceleration of the plasma are due to wave dissipation and wave pressure gradients, respectively. The dissipation process is described by a fully developed turbulent cascade of counter-propagating waves. We adopt a unified approach for calculating the wave dissipation in both open and closed magnetic field lines, allowing for a self-consistent treatment of any magnetic topology. Wave dissipation is the only heating mechanism assumed in the model, and no geometric heating functions are invoked. Electron heat conduction and radiative cooling are also included. We demonstrate that the large-scale, steady-state (in the co-rotating frame) properties of the solar environment are reproduced, using three adjustable parameters: the Poynting flux of chromospheric Alfv\'en waves, the perpendicular correlation length of the turbulence, and a pseudo-reflection coefficient. We compare model results for Carrington Rotation 2063 (November-December 2007) to remote observations in the EUV and X-ray ranges from STEREO, SOHO and Hinode spacecraft, as well as to in-situ measurements performed by Ulysses. The model results are in good agreement with observations. This is the first global model capable of simultaneously reproducing the multi-wavelength observations of the lower corona and the wind structure beyond Earth's orbit.

The large-scale field of the Sun is well represented by its lowest energy (or potential) state. Recent observations, by comparison, reveal that many solar-type stars show large-scale surface magnetic fields that are highly non-potential - that is, they have been stressed above their lowest-energy state. This non-potential component of the surface field is neglected by current stellar wind models. The aim of this paper is to determine its effect on the coronal structure and wind. We use Zeeman-Doppler surface magnetograms of two stars - one with an almost potential, one with a non-potential surface field - to extrapolate a static model of the coronal structure for each star. We find that the stresses are carried almost exclusively in a band of uni-directional azimuthal field that is confined to mid-latitudes. Using this static solution as an initial state for an MHD wind model, we then find that the final state is determined primarily by the potential component of the surface magnetic field. The band of azimuthal field must be confined close to the stellar surface as it is not compatible with a steady-state wind. By artificially increasing the stellar rotation rate we demonstrate that the observed azimuthal fields can not be produced by the action of the wind but must be due to processes at or below the stellar surface. We conclude that the background winds of solar-like stars are largely unaffected by these highly-stressed surface fields. Nonetheless, the increased flare activity and associated coronal mass ejections that may be expected to accompany such highly-stressed fields may have a significant impact on any surrounding planets.

We present a new global model of the solar corona, including the low corona, the transition region and the top of chromosphere. The realistic 3D magnetic field is simulated using the data from the photospheric magnetic field measurements. The distinctive feature of the new model is incorporating the MHD Alfven wave turbulence. We assume this turbulence and its non-linear dissipation to be the only momentum and energy source for heating the coronal plasma and driving the solar wind. The difference between the turbulence dissipation efficiency in coronal holes and that in closed field regions is because the non-linear cascade rate degrades in strongly anisotropic turbulence in coronal holes (no inward propagating wave), thus resulting in colder coronal holes with the bi-modal solar wind originating from them. The detailed presentation of the theoretical model is illustrated with the synthetic images for multi-wavelength EUV emission compared with the observations from SDO AIA and Stereo EUVI instruments for the Carrington rotation 2107.

Magnetic loops are building blocks of the closed-field corona. While active region loops are readily seen in images taken at EUV and X-ray wavelengths, quiet Sun loops are seldom identifiable and therefore difficult to study on an individual basis. The first analysis of solar minimum (Carrington Rotation 2077) quiet Sun (QS) coronal loops utilizing a novel technique called the Michigan Loop Diagnostic Technique (MLDT) is presented. This technique combines Differential Emission Measure Tomography (DEMT) and a potential field source surface (PFSS) model, and consists of tracing PFSS field lines through the tomographic grid on which the Local Differential Emission Measure (LDEM) is determined. As a result, the electron temperature Te and density Ne at each point along each individual field line can be obtained. Using data from STEREO/EUVI and SOHO/MDI, the MLDT identifies two types of QS loops in the corona: so-called "up" loops in which the temperature increases with height, and so-called "down" loops in which the temperature decreases with height. Up loops are expected, however, down loops are a surprise, and furthermore, they are ubiquitous in the low-latitude corona. Up loops dominate the QS at higher latitudes. The MLDT allows independent determination of the empirical pressure and density scale heights, and the differences between the two remain to be explained. The down loops appear to be a newly discovered property of the solar minimum corona that may shed light on the physics of coronal heating. The results are shown to be robust to the calibration uncertainties of the EUVI instrument.

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 carry out a numerical simulation depicting the effects of a Coronal Mass Ejection (CME) event on a close-in giant planet in an extrasolar system. We drive the CME in a similar manner as in simulations of space weather events on Earth. The simulation includes the planetary orbital motion, which leads to the forming of a comet-like planetary magnetotail which is oriented almost perpendicular to the direction of propagation of the CME. The combination of this feature and the fact that the CME does not expand much by the time it reaches the planet leads to a unique CME-magnetosphere interaction, where the CME itself is highly affected by the presence of the planetary magnetosphere. We find that the planet is well-shielded from CME penetration, even for a relatively weak internal magnetic field. The planetary angular momentum loss associated with such an event is negligible compared to the total planetary angular momentum. We also find that the energy which is deposited in the magnetosphere is much higher than in the case of the Earth, and our simulation suggests there is a large-scale change in the orientation of the magnetosphere-ionosphere current system during the CME event.

We present a numerical investigation of the coronal evolution of a coronal mass ejection (CME) on 2005 August 22 using a 3-D thermodynamics magnetohydrodynamic model, the SWMF. The source region of the eruption was anemone active region (AR) 10798, which emerged inside a coronal hole. We validate our modeled corona by producing synthetic extreme ultraviolet (EUV) images, which we compare to EIT images. We initiate the CME with an out-of-equilibrium flux rope with an orientation and chirality chosen in agreement with observations of a H-alpha filament. During the eruption, one footpoint of the flux rope reconnects with streamer magnetic field lines and with open field lines from the adjacent coronal hole. It yields an eruption which has a mix of closed and open twisted field lines due to interchange reconnection and only one footpoint line-tied to the source region. Even with the large-scale reconnection, we find no evidence of strong rotation of the CME as it propagates. We study the CME deflection and find that the effect of the Lorentz force is a deflection of the CME by about 3 deg/Rsun towards the East during the first 30 minutes of the propagation. We also produce coronagraphic and EUV images of the CME, which we compare with real images, identifying a dimming region associated with the reconnection process. We discuss the implication of our results for the arrival at Earth of CMEs originating from the limb and for models to explain the presence of open field lines in magnetic clouds.

We present a focused parameter study of solar wind - magnetosphere interaction for the young Sun and Earth, $~3.5$ Ga ago, that relies on magnetohydrodynamic (MHD) simulations for both the solar wind and the magnetosphere. By simulating the quiescent young Sun and its wind we are able to propagate the MHD simulations up to Earth's magnetosphere and obtain a physically realistic solar forcing of it. We assess how sensitive the young solar wind is to changes in the coronal base density, sunspot placement and magnetic field strength, dipole magnetic field strength and the Sun's rotation period. From this analysis we obtain a range of plausible solar wind conditions the paleomagnetosphere may have been subject to. Scaling relationships from the literature suggest that a young Sun would have had a mass flux different from the present Sun. We evaluate how the mass flux changes with the aforementioned factors and determine the importance of this and several other key solar and magnetospheric variables with respect to their impact on the paleomagnetosphere. We vary the solar wind speed, density, interplanetary magnetic field strength and orientation as well as Earth's dipole magnetic field strength and tilt in a number of steady-state scenarios that are representative of young Sun-Earth interaction. This study is done as a first step of a more comprehensive effort towards understanding the implications of Sun-Earth interaction for planetary atmospheric evolution.

We carry out the first time-dependent numerical MagnetoHydroDynamic modeling of an extrasolar planetary system to study the interaction of the stellar magnetic field and wind with the planetary magnetosphere and outflow. We base our model on the parameters of the HD 189733 system, which harbors a close-in giant planet. Our simulation reveals a highly structured stellar corona characterized by sectors with different plasma properties. The star-planet interaction varies in magnitude and complexity, depending on the planetary phase, planetary magnetic field strength, and the relative orientation of the stellar and planetary fields. It also reveals a long, comet-like tail which is a result of the wrapping of the planetary magnetospheric tail by its fast orbital motion. A reconnection event occurs at a specific orbital phase, causing mass loss from the planetary magnetosphere that can generate a hot spot on the stellar surface. The simulation also shows that the system has sufficient energy to produce hot-spots observed in Ca II lines in giant planet hosting stars. However, the short duration of the reconnection event suggests that such SPI cannot be observed persistently.

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.

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.

We present a numerical Magnetohydrodynamic (MHD) study of the dependence of stellar mass and angular momentum- loss rates on the orbital distance to close-in giant planets. We find that the mass loss rate drops by a factor of $\approx$1.5-2, while the angular momentum loss rate drops by a factor of $\approx$4 as the distance decreases past the Alfv\'en surface. This reduction in angular momentum loss is due to the interaction between the stellar and planetary Alfv\'en surfaces, which modifies the global structure of the stellar corona and stellar wind on the hemisphere facing the planet, as well as the opposite hemisphere. The simulation also shows that the magnitude of change in angular momentum loss rate depends mostly on the strength of the planetary magnetic field and not on its polarity. The interaction however, begins at greater separation if the overall field topology of the star and the planet are of anti-aligned. Our results are consistent with evidence for excess angular momentum in stars harboring close-in giant planets, and show that the reduction in wind-driven angular momentum loss can compete with, and perhaps dominate, spin-up due to tidal interaction.

We perform a numerical simulation of the corona of the young, rapidly rotating K0 dwarf AB Doradus using a global MHD model. The model is driven by a surface map of the radial magnetic field constructed using Zeeman-Doppler Imaging. We find that the global structure of the stellar corona is dominated by strong azimuthal tangling of the magnetic field due to the rapid rotation. The MHD solution enables us to calculate realistic Alfv\'en surfaces and we can therefore estimate the stellar mass loss rate and angular momentum loss rate without making undue theoretical simplifications. We consider three cases, parametrized by the base density of the corona, that span the range of possible solutions for the system. We find that overall, the mass and angular-momentum loss rates are higher than in the solar case; the mass loss rates are 10 to 500 times higher, and the angular momentum loss rate can be up to $3\times{10}^4$ higher than present day solar values. Our simulations show that this model can be use to constrain the wide parameter space of stellar systems. It also shows that an MHD approach can provide more information about the physical system over the commonly used potential field extrapolation.

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 perform a numerical study of the evolution of a Coronal Mass Ejection (CME) and its interaction with the coronal magnetic field based on the May 12, 1997, CME event using a global MagnetoHydroDynamic (MHD) model for the solar corona. The ambient solar wind steady-state solution is driven by photospheric magnetic field data, while the solar eruption is obtained by superimposing an unstable flux rope onto the steady-state solution. During the initial stage of CME expansion, the core flux rope reconnects with the neighboring field, which facilitates lateral expansion of the CME footprint in the low corona. The flux rope field also reconnects with the oppositely orientated overlying magnetic field in the manner of the breakout model. During this stage of the eruption, the simulated CME rotates counter-clockwise to achieve an orientation that is in agreement with the interplanetary flux rope observed at 1 AU. A significant component of the CME that expands into interplanetary space comprises one of the side lobes created mainly as a result of reconnection with the overlying field. Within 3 hours, reconnection effectively modifies the CME connectivity from the initial condition where both footpoints are rooted in the active region to a situation where one footpoint is displaced into the quiet Sun, at a significant distance ($\approx 1R_\odot$) from the original source region. The expansion and rotation due to interaction with the overlying magnetic field stops when the CME reaches the outer edge of the helmet streamer belt, where the field is organized on a global scale. The simulation thus offers a new view of the role reconnection plays in rotating a CME flux rope and transporting its footpoints while preserving its core structure.

We present a three-dimensional simulation of the corona of an FK Com-type rapidly rotating G giant using a magnetohydrodynamic model that was originally developed for the solar corona in order to capture the more realistic, non-potential coronal structure. We drive the simulation with surface maps for the radial magnetic field obtained from a stellar dynamo model of the FK Com system. This enables us to obtain the coronal structure for different field topologies representing different periods of time. We find that the corona of such an FK Com-like star, including the large scale coronal loops, is dominated by a strong toroidal component of the magnetic field. This is a result of part of the field being dragged by the radial outflow, while the other part remains attached to the rapidly rotating stellar surface. This tangling of the magnetic field,in addition to a reduction in the radial flow component, leads to a flattening of the gas density profile with distance in the inner part of the corona. The three-dimensional simulation provides a global view of the coronal structure. Some aspects of the results, such as the toroidal wrapping of the magnetic field, should also be applicable to coronae on fast rotators in general, which our study shows can be considerably different from the well-studied and well-observed solar corona. Studying the global structure of such coronae should also lead to a better understanding of their related stellar processes, such as flares and coronal mass ejections, and in particular, should lead to an improved understanding of mass and angular momentum loss from such systems.

In this work we describe our implementation of a thermodynamic energy equation into the global corona model of the Space Weather Modeling Framework (SWMF), and its development into the new Lower Corona (LC) model. This work includes the integration of the additional energy transport terms of coronal heating, electron heat conduction, and optically thin radiative cooling into the governing magnetohydrodynamic (MHD) energy equation. We examine two different boundary conditions using this model; one set in the upper transition region (the Radiative Energy Balance model), as well as a uniform chromospheric condition where the transition region can be modeled in its entirety. Via observation synthesis from model results and the subsequent comparison to full sun extreme ultraviolet (EUV) and soft X-Ray observations of Carrington Rotation (CR) 1913 centered on Aug 27, 1996, we demonstrate the need for these additional considerations when using global MHD models to describe the unique conditions in the low corona. Through multiple simulations we examine ability of the LC model to asses and discriminate between coronal heating models, and find that a relative simple empirical heating model is adequate in reproducing structures observed in the low corona. We show that the interplay between coronal heating and electron heat conduction provides significant feedback onto the 3D magnetic topology in the low corona as compared to a potential field extrapolation, and that this feedback is largely dependent on the amount of mechanical energy introduced into the corona.

Since the first discovery of an extrasolar planetary system more than a decade ago, hundreds more have been discovered. Surprisingly, many of these systems harbor Jupiter-class gas giants located close to the central star, at distances of 0.1 AU or less. Observations of chromospheric 'hot spots' that rotate in phase with the planetary orbit, and elevated stellar X-ray luminosities,suggest that these close-in planets significantly affect the structure of the outer atmosphere of the star through interactions between the stellar magnetic field and the planetary magnetosphere. Here we carry out the first detailed three-dimensional MagnetoHydroHynamics (MHD) simulation containing the two magnetic bodies and explore the consequences of such interactions on the steady-state coronal structure. The simulations reproduce the observable features of 1) increase in the total X-ray luminosity, 2) appearance of coronal hot spots, and 3) phase shift of these spots with respect to the direction of the planet. The proximate cause of these is an increase in the density of coronal plasma in the direction of the planet, which prevents the corona from expanding and leaking away this plasma via a stellar wind. The simulations produce significant low temperature heating. By including dynamical effects, such as the planetary orbital motion, the simulation should better reproduce the observed coronal heating.

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 simulate the effect of latitudinal variations in the location of star spots, as well as their magnetic field strength, on stellar angular momentum loss to the stellar wind. We use the Michigan solar corona global MagnetoHydroDynamic model, which incorporates realistic relation between the magnetic field topology and the wind distribution. We find that the spots location significantly affects the stellar wind structure, and as a result, the total mass loss rate and angular momentum loss rate. In particular, we find that the angular momentum loss rate is controlled by the mass flux when spots are located at low latitudes but is controlled by an increased plasma density between the stellar surface and the Alfven surface when spots are located at high latitudes. Our results suggest that there might be a feedback mechanism between the magnetic field distribution, wind distribution, angular momentum loss through the wind, and the motions at the convection zone that generate the magnetic field. This feedback might explain the role of coronal magnetic fields in stellar dynamos.

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)

This study examines the tail disconnection event on April 20, 2007 on comet 2P/Encke, caused by a coronal mass ejection (CME) at a heliocentric distance of 0.34 AU. During their interaction, both the CME and the comet are visible with high temporal and spatial resolution by the STEREO-A spacecraft. Previously, only current sheets or shocks have been accepted as possible reasons for comet tail disconnections, so it is puzzling that the CME caused this event. The MHD simulation presented in this work reproduces the interaction process and demonstrates how the CME triggered a tail disconnection in the April 20 event. It is found that the CME disturbs the comet with a combination of a $180^\circ$ sudden rotation of the interplanetary magnetic field (IMF), followed by a $90^\circ$ gradual rotation. Such an interpretation applies our understanding of solar wind-comet interactions to determine the \textit{in situ} IMF orientation of the CME encountering Encke.

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.

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.

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.