Source author record

M. L. Khodachenko

M. L. Khodachenko 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

24works
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

24 published item(s)

preprint2021arXiv

Simulation of 10830 Å absorption with a 3D hydrodynamic model reveals the solar He abundance in upper atmosphere of WASP-107b

Transmission spectroscopy of WASP-107b revealed 7-8% absorption at the position of metastable HeI triplet at 10830 Å in Doppler velocity range of [-20; 10] km/s, which is stronger than that measured in other exoplanets. With a dedicated 3D self-consistent hydrodynamic multi-fluid model we calculated the expanding upper atmosphere of WASP-107b and reproduced within the observations accuracy the measured HeI absorption profiles, constraining the stellar XUV flux to (6-10) erg cm-2 s-1 at 1 a.u., and the upper atmosphere helium abundance He/H to 0.075-0.15. The radiation pressure acting on the metastable HeI atoms was shown to be an important factor affecting the shape of the absorption profiles. Its effect is counterbalanced by the processes of collisional depopulation of the HeI metastable state. Altogether, the observed HeI absorption in WASP-107b can be interpreted with the expected reasonable parameters of the stellar-planetary system and appropriate account of the electron and atom impact processes.

preprint2021arXiv

The GAPS Programme at TNG. XXXII. The revealing non-detection of metastable HeI in the atmosphere of the hot Jupiter WASP-80b

The hot Jupiter WASP-80b has been identified as a possible excellent target for detecting and measuring HeI absorption in the upper atmosphere. We observed 4 primary transits of WASP-80b in the optical and near-IR using the HARPS-N and GIANO-B high-resolution spectrographs, focusing on the HeI triplet. We further employed a three-dimensional hydrodynamic aeronomy model to understand the observational results. We did not find any signature of planetary absorption at the position of the HeI triplet with an upper limit of 0.7% (i.e. 1.11 planetary radii; 95% confidence level). We re-estimated the stellar high-energy emission that we combined with a stellar photospheric model to generate the input for the hydrodynamic modelling. We obtained that, assuming a solar He to H abundance ratio, HeI absorption should have been detected. Considering a stellar wind 25 times weaker than solar, we could reproduce the non-detection only assuming a He to H abundance ratio about 16 times smaller than solar. Instead, considering a stellar wind 10 times stronger than solar, we could reproduce the non-detection only with a He to H abundance ratio about 10 times smaller than solar. We attempted to understand this result by collecting all past HeI measurements looking for correlations with stellar high-energy emission and planetary gravity, but without finding any. WASP-80b is not the only planet with a sub-solar estimated He to H abundance ratio, suggesting the presence of efficient physical mechanisms (e.g. phase separation, magnetic fields) capable of significantly modifying the He to H content in the upper atmosphere of hot Jupiters. The planetary macroscopic properties and the shape of the stellar spectral energy distribution are not sufficient for predicting the presence or absence of detectable metastable He in a planetary atmosphere, as also the He abundance appears to play a major role.

preprint2020arXiv

Case study on the identification and classification of small-scale flow patterns in flaring active region

We propose a novel methodology to identity flows in the solar atmosphere and classify their velocities as either supersonic, subsonic, or sonic. The proposed methodology consists of three parts. First, an algorithm is applied to the Solar Dynamics Observatory (SDO) image data to locate and track flows, resulting in the trajectory of each flow over time. Thereafter, the differential emission measure inversion method is applied to six AIA channels along the trajectory of each flow in order to estimate its background temperature and sound speed. Finally, we classify each flow as supersonic, subsonic, or sonic by performing simultaneous hypothesis tests on whether the velocity bounds of the flow are larger, smaller, or equal to the background sound speed. The proposed methodology was applied to the SDO image data from the 171 Å spectral line for the date 6 March 2012 from 12:22:00 to 12:35:00 and again for the date 9 March 2012 from 03:00:00 to 03:24:00. Eighteen plasma flows were detected, 11 of which were classified as supersonic, 3 as subsonic, and 3 as sonic at a $70\%$ level of significance. Out of all these cases, 2 flows cannot be strictly ascribed to one of the respective categories as they change from the subsonic state to supersonic and vice versa. We labelled them as a subclass of transonic flows. The proposed methodology provides an automatic and scalable solution to identify small-scale flows and to classify their velocities as either supersonic, subsonic, or sonic. We identified and classified small-scale flow patterns in flaring loops. The results show that the flows can be classified into four classes: sub-, super-, trans-sonic, and sonic. The detected flows from AIA images can be analyzed in combination with the other high-resolution observational data, such as Hi-C 2.1 data, and be used for the development of theories of the formation of flow patterns.

preprint2020arXiv

Evolution of the Earth's Polar Outflow From Mid-Archean to Present

The development of habitable conditions on Earth is tightly connected to the evolution of its atmosphere which is strongly influenced by atmospheric escape. We investigate the evolution of the polar ion outflow from the open field line bundle which is the dominant escape mechanism for the modern Earth. We perform Direct Simulation Monte Carlo (DSMC) simulations and estimate the upper limits on escape rates from the Earth's open field line bundle starting from three gigayears ago (Ga) to present assuming the present-day composition of the atmosphere. We perform two additional simulations with lower mixing ratios of oxygen of 1% and 15% to account for the conditions shortly after the Great Oxydation Event (GOE). We estimate the maximum loss rates due to polar outflow three gigayears ago of $3.3 \times10^{27}$ s$^{-1}$ and $2.4 \times 10^{27}$ s$^{-1}$ for oxygen and nitrogen, respectively. The total integrated mass loss equals to 39% and 10% of the modern atmosphere's mass, for oxygen and nitrogen, respectively. According to our results, the main factors that governed the polar outflow in the considered time period are the evolution of the XUV radiation of the Sun and the atmosphere's composition. The evolution of the Earth's magnetic field plays a less important role. We conclude that although the atmosphere with the present-day composition can survive the escape due to polar outflow, a higher level of CO$_2$ between 3.0 and 2.0~Ga is likely necessary to reduce the escape.

preprint2020arXiv

Global 3D hydrodynamic modeling of absorption in Lyα and He 10830 A lines at transits of GJ3470b

Warm Neptune GJ3470b has been recently observed in 23S-23P transition of metastable helium, yielding absorption of about 1% in Doppler velocity range of [-40; 10] km/s. Along with previous detection of absorption in Lyα with depth of 20-40% in the blue and red wings of the line, it offers a complex target for simulation and testing of the current models. Obtained results suggest that absorption in both these lines comes from interaction of expanding upper planetary atmosphere with stellar plasma wind, allowing to constrain the stellar plasma parameters and the helium abundance in planet atmosphere.

preprint2020arXiv

Three-dimensional hydrodynamic simulations of the upper atmosphere of $π$ Men c: comparison with Ly$α$ transit observations

Aims: We aim at constraining the conditions of the wind and high-energy emission of the host star reproducing the non-detection of Ly$α$ planetary absorption. Methods: We model the escaping planetary atmosphere, the stellar wind, and their interaction employing a multi-fluid, three-dimensional hydrodynamic code. We assume a planetary atmosphere composed of hydrogen and helium. We run models varying the stellar high-energy emission and stellar mass-loss rate, further computing for each case the Ly$α$ synthetic planetary atmospheric absorption and comparing it with the observations. Results: We find that a non-detection of Ly$α$ in absorption employing the stellar high-energy emission estimated from far-ultraviolet and X-ray data requires a stellar wind with a stellar mass-loss rate about six times lower than solar. This result is a consequence of the fact that, for $π$ Men c, detectable Ly$α$ absorption can be caused exclusively by energetic neutral atoms, which become more abundant with increasing the velocity and/or the density of the stellar wind. By considering, instead, that the star has a solar-like wind, the non-detection requires a stellar ionising radiation about four times higher than estimated. This is because, despite the fact that a stronger stellar high-energy emission ionises hydrogen more rapidly, it also increases the upper atmosphere heating and expansion, pushing the interaction region with the stellar wind farther away from the planet, where the planet atmospheric density that remains neutral becomes smaller and the production of energetic neutral atoms less efficient. Conclusions: Comparing the results of our grid of models with what is expected and estimated for the stellar wind and high-energy emission, respectively, we support the idea that the atmosphere of $π$ Men c is likely not hydrogen-dominated.

preprint2019arXiv

Global 3D hydrodynamic modeling of in-transit Lyα absorption of GJ436b

Using a global 3D, fully self-consistent, multi-fluid hydrodynamic model, we simulate the escaping upper atmosphere of the warm Neptune GJ436b, driven by the stellar XUV radiation impact and gravitational forces and interacting with the stellar wind. Under the typical parameters of XUV flux and stellar wind plasma expected for GJ436, we calculate in-transit absorption in Lyα and find that it is produced mostly by Energetic Neutral Atoms outside of the planetary Roche lobe, due to the resonant thermal line broadening. At the same time, the influence of radiation pressure has been shown to be insignificant. The modelled absorption is in good agreement with the observations and reveals such features as strong asymmetry between blue and red wings of the absorbed Lyα line profile, deep transit depth in the high velocity blue part of the line reaching more than 70%, and the timing of early ingress. On the other hand, the model produces significantly deeper and longer egress than in observations, indicating that there might be other processes and factors, still not accounted, that affect the interaction between the planetary escaping material and the stellar wind. At the same time, it is possible that the observational data, collected in different measurement campaigns, are affected by strong variations of the stellar wind parameters between the visits, and therefore, they cannot be reproduced altogether with the single set of model parameters.

preprint2016arXiv

Thermal mass loss of protoplanetary cores with hydrogen-dominated atmospheres: The influences of ionization and orbital distance

We investigate the loss rates of the hydrogen atmospheres of terrestrial planets with a range of masses and orbital distances by assuming a stellar extreme ultraviolet (EUV) luminosity that is 100 times stronger than that of the current Sun. We apply a 1D upper atmosphere radiation absorption and hydrodynamic escape model that takes into account ionization, dissociation and recombination to calculate hydrogen mass loss rates. We study the effects of the ionization, dissociation and recombination on the thermal mass loss rates of hydrogen-dominated super-Earths and compare the results to those obtained by the energy-limited escape formula which is widely used for mass loss evolution studies. Our results indicate that the energy-limited formula can to a great extent over- or underestimate the hydrogen mass loss rates by amounts that depend on the stellar EUV flux and planetary parameters such as mass, size, effective temperature, and EUV absorption radius.

preprint2015arXiv

Atmosphere expansion and mass loss of close-orbit giant exoplanets heated by stellar XUV: I. Modeling of hydrodynamic escape of upper atmospheric material

In the present series of papers we propose a consistent description of the mass loss process. To study the effects of intrinsic magnetic field of a close-orbit giant exoplanet (so-called Hot Jupiter) on the atmospheric material escape and formation of planetary inner magnetosphere in a comprehensive way, we start with a hydrodynamic model of an upper atmosphere expansion presented in this paper. While considering a simple hydrogen atmosphere model, we focus on selfconsistent inclusion of the effects of radiative heating and ionization of the atmospheric gas with its consequent expansion in the outer space. Primary attention is paid to investigation of the role of specific conditions at the inner and outer boundaries of the simulation domain, under which different regimes of material escape (free- and restricted- flow) are formed. Comparative study of different processes, such as XUV heating, material ionization and recombination, H3+ cooling, adiabatic and Lyman-alpha cooling, Lyman-alpha reabsorption is performed. We confirm basic consistence of the outcomes of our modeling with the results of other hydrodynamic models of expanding planetary atmospheres. In particular, we obtain that under the typical conditions of an orbital distance 0.05 AU around a Sun-type star a Hot Jupiter plasma envelope may reach maximum temperatures up to ~9000K with a hydrodynamic escape speed ~9 km/s resulting in the mass loss rates ~(4-7)*10^10 g*s . In the range of considered stellar-planetary parameters and XUV fluxes that is close to mass loss in the energy limited case. The inclusion of planetary intrinsic magnetic fields in the model is a subject of the following up paper (Paper II).

preprint2015arXiv

Investigation of scaling properties of a thin current sheet by means of particle trajectories study

A thin current sheet (TCS), with the width of an order of thermal proton gyroradius, appears a fundamental physical object which plays an important role in structuring of major magnetospheric current systems (magnetotail, magnetodisk, etc.). The TCSs are nowadays under extensive study by means of space missions and theoretical models. We consider a simple model of the TCS separating two half-spaces occupied by a homogenous magnetic field of opposite sign tangential to the TCS; a small normal component of the magnetic field is prescribed. An analytical solution for the electric current and plasma density in the close vicinity of the TCS has been obtained and compared with numerical simulation. The number density and the electric current profiles have two maxima each. The characteristic spatial scale $z_S$ of the maxima location was investigated as a function of initial pitch-angle of an incoming charge particle. The effect of the thermal dispersion of the incoming proton beam have been taken into consideration.

preprint2014arXiv

A search for flares and mass ejections on young late-type stars in the open cluster Blanco-1

We present a search for stellar activity (flares and mass ejections) in a sample of 28 stars in the young open cluster Blanco-1. We use optical spectra obtained with ESO's VIMOS multi-object spectrograph installed on the VLT. From the total observing time of $\sim$ 5 hours, we find four H$α$ flares but no distinct indication of coronal mass ejections (CMEs) on the investigated dK-dM stars. Two flares show "dips" in their light-curves right before their impulsive phases which are similar to previous discoveries in photometric light-curves of active dMe stars. We estimate an upper limit of $<$4 CMEs per day per star and discuss this result with respect to a semi- empirical estimation of the CME rate of main-sequence stars. We find that we should have detected at least one CME per star with a mass of 1-15$\times10^{16}$ g depending on the star's X-ray luminosity, but the estimated H$α$ fluxes associated with these masses are below the detection limit of our observations. We conclude that the parameter which mainly influences the detection of stellar CMEs using the method of Doppler-shifted emission caused by moving plasma is not the spectral resolution or velocity but the flux or mass of the CME.

preprint2014arXiv

Magnetic moment and plasma environment of HD 209458b as determined from Ly$α$ observations

Transit observations of HD 209458b in the stellar Lyman-$α$ (Ly$α$) line revealed strong absorption in both blue and red wings of the line interpreted as hydrogen atoms escaping from the planet's exosphere at high velocities. The following sources for the absorption were suggested: acceleration by the stellar radiation pressure, natural spectral line broadening, charge exchange with stellar wind. We reproduce the observation by means of modelling that includes all aforementioned processes. Our results support a stellar wind with a velocity of $\approx400$ km$\times$s$^{-1}$ at the time of the observation and a planetary magnetic moment of $\approx 1.6 \times 10^{26}$ A$\times$m$^2$.

preprint2014arXiv

Stellar wind interaction and pick-up ion escape of the Kepler-11 "super-Earths"

We study the interactions between stellar wind and the extended hydrogen-dominated upper atmospheres of planets and the resulting escape of planetary pick-up ions from the 5 "super-Earths" in the compact Kepler-11 system and compare the escape rates with the efficiency of the thermal escape of neutral hydrogen atoms. Assuming the stellar wind of Kepler-11 is similar to the solar wind, we use a polytropic 1D hydrodynamic wind model to estimate the wind properties at the planetary orbits. We apply a Direct Simulation Monte Carlo Model to model the hydrogen coronae and the stellar wind plasma interaction around Kepler-11b-f within a realistic expected heating efficiency range of 15-40%. The same model is used to estimate the ion pick-up escape from the XUV heated and hydrodynamically extended upper atmospheres of Kepler-11b-f. From the interaction model we study the influence of possible magnetic moments, calculate the charge exchange and photoionization production rates of planetary ions and estimate the loss rates of pick-up H+ ions for all five planets. We compare the results between the five "super-Earths" and in a more general sense also with the thermal escape rates of the neutral planetary hydrogen atoms. Our results show that for all Kepler-11b-f exoplanets, a huge neutral hydrogen corona is formed around the planet. The non-symmetric form of the corona changes from planet to planet and is defined mostly by radiation pressure and gravitational effects. Non-thermal escape rates of pick-up ionized hydrogen atoms for Kepler-11 "super-Earths" vary between approximately 6.4e30 1/s and 4.1e31 1/s depending on the planet's orbital location and assumed heating efficiency. These values correspond to non-thermal mass loss rates of approximately 1.07e7 g/s and 6.8e7 g/s respectively, which is a few percent of the thermal escape rates.

preprint2014arXiv

The PLATO 2.0 Mission

PLATO 2.0 has recently been selected for ESA's M3 launch opportunity (2022/24). Providing accurate key planet parameters (radius, mass, density and age) in statistical numbers, it addresses fundamental questions such as: How do planetary systems form and evolve? Are there other systems with planets like ours, including potentially habitable planets? The PLATO 2.0 instrument consists of 34 small aperture telescopes (32 with 25 sec readout cadence and 2 with 2.5 sec candence) providing a wide field-of-view (2232 deg2) and a large photometric magnitude range (4-16 mag). It focusses on bright (4-11 mag) stars in wide fields to detect and characterize planets down to Earth-size by photometric transits, whose masses can then be determined by ground-based radial-velocity follow-up measurements. Asteroseismology will be performed for these bright stars to obtain highly accurate stellar parameters, including masses and ages. The combination of bright targets and asteroseismology results in high accuracy for the bulk planet parameters: 2%, 4-10% and 10% for planet radii, masses and ages, respectively. The planned baseline observing strategy includes two long pointings (2-3 years) to detect and bulk characterize planets reaching into the habitable zone (HZ) of solar-like stars and an additional step-and-stare phase to cover in total about 50% of the sky. PLATO 2.0 will observe up to 1,000,000 stars and detect and characterize hundreds of small planets, and thousands of planets in the Neptune to gas giant regime out to the HZ. It will therefore provide the first large-scale catalogue of bulk characterized planets with accurate radii, masses, mean densities and ages. This catalogue will include terrestrial planets at intermediate orbital distances, where surface temperatures are moderate. Coverage of this parameter range with statistical numbers of bulk characterized planets is unique to PLATO 2.0.

preprint2013arXiv

Dynamics of coronal rain and descending plasma blobs in solar prominences: I. Fully ionised case

Observations of active regions and limb prominences often show cold, dense blobs descending with an acceleration smaller than that of free fall. The dynamics of these condensations falling in the solar corona is investigated in this paper using a simple fully ionised plasma model. We find that the presence of a heavy condensation gives rise to a dynamical rearrangement of the coronal pressure that results in the formation of a large pressure gradient that opposes gravity. Eventually this pressure gradient becomes so large that the blob acceleration vanishes or even points upwards. Then, the blob descent is characterised by an initial acceleration phase followed by an essentially constant velocity phase. These two stages can be identified in published time-distance diagrams of coronal rain events. Both the duration of the first stage and the velocity attained by the blob increase for larger values of the ratio of blob to coronal density, for larger blob mass, and for smaller coronal temperature. Dense blobs are characterised by a detectable density growth (up to 60% in our calculations) and by a steepening of the density in their lower part, that could lead to the formation of a shock. They also emit sound waves that could be detected as small intensity changes with periods of the order of 100 s and lasting between a few and about ten periods. Finally, the curvature of the falling path is only relevant when a very dense blob falls along inclined magnetic field lines.

preprint2013arXiv

Torsional Alfvén waves in solar partially ionized plasma: effects of neutral helium and stratification

Ion-neutral collisions may lead to the damping of Alfven waves in chromospheric and prominence plasmas. Neutral helium atoms enhance the damping in certain temperature interval, where the ratio of neutral helium and neutral hydrogen atoms is increased. Therefore, the height-dependence of ionization degrees of hydrogen and helium may influence the damping rate of Alfven waves. We aim to study the effect of neutral helium in the damping of Alfven waves in stratified partially ionized plasma of the solar chromosphere. We consider a magnetic flux tube, which is expanded up to 1000 km height and then becomes vertical due to merging with neighboring tubes, and study the dynamics of linear torsional Alfven waves in the presence of neutral hydrogen and neutral helium atoms. We start with three-fluid description of plasma and consequently derive single-fluid magnetohydrodynamic (MHD) equations for torsional Alfven waves. Thin flux tube approximation allows to obtain the dispersion relation of the waves in the lower part of tubes, while the spatial dependence of steady-state Alfven waves is governed by Bessel type equation in the upper part of tubes. Consecutive derivation of single-fluid MHD equations results in a new Cowling diffusion coefficient in the presence of neutral helium which is different from previously used one. We found that shorter-period (< 5 s) torsional Alfven waves damp quickly in the chromospheric network due to ion-neutral collision. On the other hand, longer-period (> 5 s) waves do not reach the transition region as they become evanescent at lower heights in the network cores. Propagation of torsional Alfven waves through the chromosphere into the solar corona should be considered with caution: low-frequency waves are evanescent due to the stratification, while high-frequency waves are damped due to ion neutral collisions.

preprint2013arXiv

XUV exposed non-hydrostatic hydrogen-rich upper atmospheres of terrestrial planets. Part I: Atmospheric expansion and thermal escape

The recently discovered low-density "super-Earths" Kepler-11b, Kepler-11f, Kepler-11d, Kepler-11e, and planets such as GJ 1214b represent most likely planets which are surrounded by dense H/He envelopes or contain deep H2O oceans also surrounded by dense hydrogen envelopes. Although these "super-Earths" are orbiting relatively close to their host stars, they have not lost their captured nebula-based hydrogen-rich or degassed volatile-rich steam protoatmospheres. Thus it is interesting to estimate the maximum possible amount of atmospheric hydrogen loss from a terrestrial planet orbiting within the habitable zone of late main sequence host stars. For studying the thermosphere structure and escape we apply a 1-D hydrodynamic upper atmosphere model which solves the equations of mass, momentum and energy conservation for a planet with the mass and size of the Earth and for a "super-Earth" with a size of 2 R_Earth and a mass of 10 M_Earth. We calculate volume heating rates by the stellar soft X-ray and EUV radiation and expansion of the upper atmosphere, its temperature, density and velocity structure and related thermal escape rates during planet's life time. Moreover, we investigate under which conditions both planets enter the blow-off escape regime and may therefore experience loss rates which are close to the energy-limited escape. Finally we discuss the results in the context of atmospheric evolution and implications for habitability of terrestrial planets in general.

preprint2013arXiv

XUV exposed, non-hydrostatic hydrogen-rich upper atmospheres of terrestrial planets II: Hydrogen coronae and ion escape

We study the interactions between the stellar wind plasma flow of a typical M star, such as GJ 436, and hydrogen-rich upper atmospheres of an Earth-like planet and a "super-Earth" with the radius of 2 R_Earth and a mass of 10 M_Earth, located within the habitable zone at ~0.24 AU. We investigate the formation of extended atomic hydrogen coronae under the influences of the stellar XUV flux (soft X-rays and EUV), stellar wind density and velocity, shape of a planetary obstacle (e.g., magnetosphere, ionopause), and the loss of planetary pick-up ions on the evolution of hydrogen-dominated upper atmospheres. Stellar XUV fluxes which are 1, 10, 50 and 100 times higher compared to that of the present-day Sun are considered and the formation of high-energy neutral hydrogen clouds around the planets due to the charge-exchange reaction under various stellar conditions have been modeled. Charge-exchange between stellar wind protons with planetary hydrogen atoms, and photoionization, leads to the production of initially cold ions of planetary origin. We found that the ion production rates for the studied planets can vary over a wide range, from ~1.0x10^{25} s^{-1} to ~5.3x10^{30} s^{-1}, depending on the stellar wind conditions and the assumed XUV exposure of the upper atmosphere. Our findings indicate that most likely the majority of these planetary ions are picked up by the stellar wind and lost from the planet. Finally, we estimate the long-time non-thermal ion pick-up escape for the studied planets and compare them with the thermal escape. According to our estimates, non-thermal escape of picked up ionized hydrogen atoms over a planet's lifetime varies between ~0.4 Earth ocean equivalent amounts of hydrogen (EO_H) to <3 EO_H and usually is several times smaller in comparison to the thermal atmospheric escape rates.

preprint2012arXiv

Cut-off wavenumber of Alfven waves in partially ionized plasmas of the solar atmosphere

Alfven wave dynamics in partially ionized plasmas of the solar atmosphere shows that there is indeed a cut-off wavenumber, i.e. the Alfven waves with wavenumbers higher than the cut-off value are evanescent. The cut-off wavenumber appears in single-fluid magnetohydrodynamic (MHD) approximation but it is absent in a multi-fluid approach. Up to now, an explanation for the existence of the cut-off wavenumber is still missing. The aim of this paper is to point out the reason for the appearance of a cut-off wavenumber in single-fluid MHD. Beginning with three-fluid equations (with electrons, protons and neutral hydrogen atoms), we performed consecutive approximations until we obtained the usual single-fluid description is obtained. We solved the dispersion relation of linear Alfven waves at each step and sought the approximation responsible of the cut-off wavenumber appearance. We have found that neglecting inertial terms significantly reduces the real part of the Alfven frequency although it never becomes zero. Therefore, the cut-off wavenumber does not exist at this stage. However, when the inertial terms together with the Hall term in the induction equation are neglected, the real part of the Alfven frequency becomes zero. The appearance of a cut-off wavenumber, when Alfven waves in partially ionized regions of the solar atmosphere are studied, is the result of neglecting inertial and Hall terms, therefore it has no physical origin.

preprint2012arXiv

Probing the Blow-Off Criteria of Hydrogen-Rich "Super-Earths"

The discovery of transiting "super-Earths" with inflated radii and known masses such as Kepler-11b-f, GJ 1214b and 55 Cnc e, indicates that these exoplanets did not lose their nebula-captured hydrogen-rich, degassed or impact-delivered protoatmospheres by atmospheric escape processes. Because hydrodynamic blow-off of atmospheric hydrogen atoms is the most efficient atmospheric escape process we apply a time-dependent numerical algorithm which is able to solve the system of 1-D fluid equations for mass, momentum, and energy conservation to investigate the criteria under which "super-Earths" with hydrogen-dominated upper atmospheres can experience hydrodynamic expansion by heating of the stellar XUV (soft X-rays and extreme ultraviolet) radiation and thermal escape via blow-off. Depending on orbit location, XUV flux, heating efficiency and the planet's mean density our results indicate that the upper atmospheres of all "super-Earths" can expand to large distances, so that besides of Kepler-11c all of them experience atmospheric mass-loss due to Roche lobe overflow. The atmospheric mass-loss of the studied "super-Earths" is one to two orders of magnitude lower compared to that of "hot Jupiters" such as HD 209458b, so that one can expect that these exoplanets cannot lose their hydrogen-envelopes during their remaining lifetimes.

preprint2011arXiv

Damping of Alfven waves in solar partially ionized plasmas: effect of neutral helium in multi-fluid approach

Chromospheric and prominence plasmas contain neutral atoms, which may change the plasma dynamics through collision with ions. Most of the atoms are neutral hydrogen, but a significant amount of neutral helium may also be present in the plasma with a particular temperature. Damping of MHD waves due to ion collision with neutral hydrogen is well studied, but the effects of neutral helium are largely unknown. We aim to study the effect of neutral helium in the damping of Alfven waves in solar partially ionized plasmas. We consider three-fluid magnetohydrodynamic (MHD) approximation, where one component is electron-proton-singly ionized helium and other two components are the neutral hydrogen and neutral helium atoms. We derive the dispersion relation of linear Alfven waves in isothermal and homogeneous plasma. Then we solve the dispersion relation and derive the damping rates of Alfven waves for different plasma parameters. The presence of neutral helium significantly enhances the damping of Alfven waves compared to the damping due to neutral hydrogen at certain values of plasma temperature (10000-40000 K) and ionization. Damping rates have a peak near the ion-neutral collision frequency, but decrease for the higher part of wave spectrum. Collision of ions with neutral helium atoms can be of importance for the damping of Alfven waves in chromospheric spicules and in prominence-corona transition regions.

preprint2011arXiv

Rossby waves and polar spots in rapidly rotating stars: Implications for stellar wind evolution

Rapidly rotating stars show short-period oscillations in magnetic activity and polar appearance of starspots. The aim of this paper is to study large-scale shallow water waves in the tachoclines of rapidly rotating stars and their connection to the periodicity and the formation of starspots at high latitudes. Shallow-water magnetohydrodynamic equations were used to study the dynamics of large-scale waves at the rapidly rotating stellar tachoclines in the presence of toroidal magnetic field. Dispersion relations and latitudinal distribution of wave modes were derived. We found that low-frequency magnetic Rossby waves tend to be located at poles, but high-frequency magnetic Poincare waves are concentrated near the equator in rapidly rotating stars. These results have important implications for the evolution of the stellar wind in young Sun-like stars. Unstable magnetic Rossby waves may lead to the local enhancement of magnetic flux at high latitudes of tachoclines in rapidly rotating stars. The enhanced magnetic flux may rise upwards owing to the magnetic buoyancy in the form of tubes and appear as starspots at polar regions. Magnetic Rossby waves may also cause observed short-term periodicity in the stellar magnetic activity.

preprint2010arXiv

Observation of kink instability during small B5.0 solar flare on 04 June, 2007

Using multi-wavelength observations of SoHO/MDI, SOT-Hinode/blue-continuum (4504 Å), G-band (4305 Å), Ca II H (3968 Å) and TRACE 171 Å, we present the observational signature of highly twisted magnetic loop in AR 10960 during the period 04:43 UT-04:52 UT at 4 June, 2007. SOT-Hinode/blue-continuum (4504 Å) observations show that penumbral filaments of positive polarity sunspot have counter-clock wise twist, which may be caused by the clock-wise rotation of the spot umbrae. The coronal loop, whose one footpoint is anchored in this sunspot, shows strong right-handed twist in chromospheric SOT-Hinode/Ca II H (3968 Å) and coronal TRACE 171 Å\, images. The length and the radius of the loop are $L\sim$80 Mm and $a\sim$4.0 Mm respectively. The distance between neighboring turns of magnetic field lines (i.e. pitch) is estimated as $\approx$ 10 Mm. The total twist angle, $Φ\sim$12$π$ (estimated for the homogeneous distribution of the twist along the loop), is much larger than the Kruskal -Shafranov instability criterion. We detected clear double structure of the loop top during 04:47-04:51 UT on TRACE 171 Å\ images, which is consistent with simulated kink instability in curved coronal loops (T{ö}r{ö}k et al. 2004). We suggest, that the kink instability of this twisted magnetic loop triggered B5.0 class solar flare, which occurred between 04:40 UT and 04:51 UT in this active region.

preprint2010arXiv

Propagation of sausage soliton in the solar lower atmosphere observed by Hinode/SOT

Acoustic waves and pulses propagating from the solar photosphere upwards may quickly develop into shocks due to the rapid decrease of atmospheric density. However, if they propagate along a magnetic flux tube, then the nonlinear steepening may be balanced by tube dispersion effects. This may result in the formation of sausage soliton. The aim of this letter is to report an observational evidence of sausage soliton in the solar chromosphere. Time series of Ca II H line obtained at the solar limb with the Solar Optical Telescope (SOT) on the board of Hinode is analysed. Observations show an intensity blob, which propagates from 500 km to 1700 km above the solar surface with the mean apparent speed of 35 km s$^{-1}$. The speed is much higher than expected local sound speed, therefore the blob can not be a simple pressure pulse. The blob speed, length to width ratio and relative intensity correspond to slow sausage soliton propagating along a magnetic tube. The blob width is increased with height corresponding to the magnetic tube expansion in the stratified atmosphere. Propagation of the intensity blob can be the first observational evidence of slow sausage soliton in the solar atmosphere.