BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

T. I. Gombosi

T. I. Gombosi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.SRastro-ph.EPphysics.space-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
22works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

22 published item(s)

preprint2016arXiv

Simulating the Environment Around Planet-Hosting Stars - I. Coronal Structure

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én 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.

0 citations0 reviewsNo code linkedOpen access
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".

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Interaction of Venus-like, M-dwarf Planets with the Stellar Wind of Their Host Star

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Alfvén Wave Turbulence as a Coronal Heating Mechanism: Simultaneously Predicting the Heating Rate and the Wave-Induced Emission Line Broadening

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én 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Magnetospheric Structure and Atmospheric Joule Heating of Habitable Planets Orbiting M-dwarf Stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Global Wave-Driven MHD Solar Model with a Unified Treatment of Open and Closed Magnetic Field Topologies

We describe, analyze and validate the recently developed Alfvén 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én 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én 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Exoplanet Transit Variability: Bow Shocks and Winds Around HD 189733b

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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&#39;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.

0 citations0 reviewsNo code linkedOpen access
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&#39;s atmosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Dynamics of Stellar Coronae Harboring Hot-jupiters I. A Time-dependent MHD Simulation of the Interplanetary Environment in the HD 189733 Planetary System

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Dynamics of Stellar Coronae Harboring Hot-jupiters II. A Space Weather Event on A Hot-jupiter

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Magnetic Structure of Rapidly Rotating FK Comae-Type Coronae

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerical Simulation of the May 12, 1997 CME Event - the Role of Magnetic Reconnection

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Coronal Structure of AB Doradus

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én 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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Impact of Hot Jupiters on the Spin-down of their Host Stars

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én surface. This reduction in angular momentum loss is due to the interaction between the stellar and planetary Alfvén 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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Interactions of the magnetospheres of stars and close-in giant planets

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 &#39;hot spots&#39; 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Study of the April 20, 2007 CME-Comet Interaction Event with an MHD Model

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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)

0 citations0 reviewsNo code linkedOpen access