Source author record

M. Opher

M. Opher appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2022arXiv

First comparison of composite 0.52-55 keV ENA spectra observed by IBEX and Cassini/INCA with simulated ENAs inferred by proton hybrid simulations downstream of the termination shock

We present a first comparison of Energetic Neutral Atom (ENA) heliosheath measurements, remotely sensed by the Interstellar Boundary Explorer (IBEX) mission and the Ion and Neutral Camera (INCA) on the Cassini mission, with modeled ENA inferred from interstellar pickup protons that have been accelerated at the termination shock, using hybrid simulations. The observed ENA intensities are an average value over the time period from 2009 to the end of 2012, along the Voyager 2 trajectory. The hybrid simulations parameters for the solar wind, interstellar pickup ions (PUIs), and magnetic field upstream of the termination shock, where Voyager 2 crossed, are based on observations. We report an energy dependent discrepancy between observed and simulated ENA fluxes, with the observed ENA fluxes, being consistently higher than the simulated ones, and discuss possible causes of this discrepancy.

preprint2021arXiv

The development of a split-tail heliosphere and the role of non-ideal processes: a comparison of the BU and Moscow models

Global models of the heliosphere are critical tools used in the interpretation of heliospheric observations. There are several three-dimensional magnetohydrodynamic (MHD) heliospheric models that rely on different strategies and assumptions. Until now only one paper has compared global heliosphere models, but without magnetic field effects. We compare the results of two different MHD models, the BU and Moscow models. Both models use identical boundary conditions to compare how different numerical approaches and physical assumptions contribute to the heliospheric solution. Based on the different numerical treatments of discontinuities, the BU model allows for the presence of magnetic reconnection, while the Moscow model does not. Both models predict collimation of the solar outflow in the heliosheath by the solar magnetic field and produce a split-tail where the solar magnetic field confines the charged solar particles into distinct north and south columns that become lobes. In the BU model, the ISM flows between the two lobes at large distances due to MHD instabilities and reconnection. Reconnection in the BU model at the port flank affects the draping of the interstellar magnetic field in the immediate vicinity of the heliopause. Different draping in the models cause different ISM pressures, yielding different heliosheath thicknesses and boundary locations, with the largest effects at high latitudes. The BU model heliosheath is 15% thinner and the heliopause is 7% more inwards at the north pole relative to the Moscow model. These differences in the two plasma solutions may manifest themselves in energetic neutral atom measurements of the heliosphere.

preprint2020arXiv

Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena in Solar and Heliospheric Plasmas

Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagency/international partnerships, and close collaborations of the solar, heliospheric, and laboratory plasma communities. These investments will deliver transformative progress in understanding magnetic reconnection and related explosive phenomena including space weather events.

preprint2020arXiv

The Confinement of the Heliosheath Plasma by the Solar Magnetic Field as Revealed by Energetic Neutral Atom Simulations

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

preprint2016arXiv

Probability of CME Impact on Exoplanets Orbiting M Dwarfs and Solar-Like Stars

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

preprint2016arXiv

Turbulence in the solar wind: spectra from Voyager 2 data at 5 AU

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

preprint2016arXiv

Using ForeCAT Deflections and Rotations to Constrain the Early Evolution of CMEs

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

preprint2015arXiv

A model of the heliosphere with jets

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

preprint2015arXiv

Constraining the Mass and the Non-Radial Drag Coefficient of a Solar Coronal Mass Ejection

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

preprint2015arXiv

Cross and magnetic helicity in the outer heliosphere from Voyager 2 observations

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

preprint2015arXiv

Global Trends of CME Deflections Based on CME and Solar Parameters

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

preprint2015arXiv

Magnetized jets driven by the sun: the structure of the heliosphere revisited

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".

preprint2015arXiv

The Heliocentric Distance Where the Deflections and Rotations of Solar Coronal Mass Ejections Occur

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

preprint2013arXiv

A Porous, Layered Heliopause

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

preprint2013arXiv

M-dwarf stellar winds: the effects of realistic magnetic geometry on rotational evolution and planets

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.

preprint2013arXiv

Propagation into the heliosheath of a large-scale solar wind disturbance bounded by a pair of shocks

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

preprint2012arXiv

The stellar wind cycles and planetary radio emission of the Tau Boo system

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.

preprint2011arXiv

Is the magnetic field in the heliosheath laminar or a turbulent bath of bubbles?

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

preprint2011arXiv

Powerful Winds from Low-Mass Stars: V374 Peg

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.

preprint2011arXiv

Understanding the Angular Momentum Loss of Low-Mass Stars: The Case of V374 Peg

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.

preprint2010arXiv

The Vector Direction of the Interstellar Magnetic Field Outside the Heliosphere

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

preprint2009arXiv

Dissipation of the sectored heliospheric magnetic field near the heliopause: a mechanism for the generation of anomalous cosmic rays

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

preprint2009arXiv

Simulations of Winds of Weak-Lined T Tauri Stars: The Magnetic Field Geometry and The Influence of the Wind on Giant Planet Migration

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.

preprint2009arXiv

Surface Alfven Wave Damping in a 3D Simulation of the Solar Wind

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.

preprint2009arXiv

Three-dimensional Numerical Simulations of Magnetized Winds of Solar-Like Stars

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)