BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

E. Khomenko

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

ResearcherUnclaimed source record
astro-ph.SRphysics.plasm-phastro-phastro-ph.GAphysics.comp-ph

Catalog footprint

What is connected

38works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

38 published item(s)

preprint2022arXiv

Generalized Fluid Models of the Braginskii Type

Several generalizations of the well-known fluid model of Braginskii (Rev. of Plasma Phys., 1965) are considered. We use the Landau collisional operator and the moment method of Grad. We focus on the 21-moment model that is analogous to the Braginskii model, and we also consider a 22-moment model. Both models are formulated for general multi-species plasmas with arbitrary masses and temperatures, where all the fluid moments are described by their evolution equations. The 21-moment model contains two "heat flux vectors" (3rd and 5th-order moments) and two "viscosity-tensors" (2nd and 4th-order moments). The Braginskii model is then obtained as a particular case of a one ion-electron plasma with similar temperatures, with de-coupled heat fluxes and viscosity-tensors expressed in a quasi-static approximation. We provide all the numerical values of the Braginskii model in a fully analytic form (together with the 4th and 5th-order moments). For multi-species plasmas, the model makes calculation of transport coefficients straightforward. Formulation in fluid moments (instead of Hermite moments) is also suitable for implementation into existing numerical codes. It is emphasized that it is the quasi-static approximation which makes some Braginskii coefficients divergent in a weakly-collisional regime. Importantly, we show that the heat fluxes and viscosity-tensors are coupled even in the linear approximation, and that the fully contracted (scalar) perturbations of the 4th-order moment, which are accounted for in the 22-moment model, modify the energy exchange rates. We also provide several Appendices, which can be useful as a guide for deriving the Braginskii model with the moment method of Grad.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Modeling the thermal conduction in the solar atmosphere with the code MANCHA3D

Thermal conductivity is one of the important mechanisms of heat transfer in the solar corona. In the limit of strongly magnetized plasma, it is typically modeled by Spitzer's expression where the heat flux is aligned with the magnetic field. This paper describes the implementation of the heat conduction into the code MANCHA3D with an aim of extending single-fluid MHD simulations from the upper convection zone into the solar corona. Two different schemes to model heat conduction are implemented: (1) a standard scheme where a parabolic term is added to the energy equation, and (2) a scheme where the hyperbolic heat flux equation is solved. The first scheme limits the time step due to the explicit integration of a parabolic term, which makes the simulations computationally expensive. The second scheme solves the limitations on the time step by artificially limiting the heat conduction speed to computationally manageable values. The validation of both schemes is carried out with standard tests in one, two, and three spatial dimensions. Furthermore, we implement the model for heat flux derived by Braginskii (1965) in its most general form, when the expression for the heat flux depends on the ratio of the collisional to cyclotron frequencies of the plasma, and, therefore on the magnetic field strength. Additionally, our implementation takes into account the heat conduction in parallel, perpendicular, and transverse directions, and provides the contributions from ions and electrons separately. The model also transitions smoothly between field-aligned conductivity and isotropic conductivity for regions with a low or null magnetic field. Finally, we present a two-dimensional test for heat conduction using realistic values of the solar atmosphere where we prove the robustness of the two schemes implemented.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Acoustic-gravity wave propagation characteristics in 3D radiation hydrodynamic simulations of the solar atmosphere

There has been tremendous progress in the degree of realism of three-dimensional radiation magneto-hydrodynamic simulations of the solar atmosphere in the past decades. Four of the most frequently used numerical codes are Bifrost, CO5BOLD, MANCHA3D, and MURaM. Here we test and compare the wave propagation characteristics in model runs from these four codes by measuring the dispersion relation of acoustic-gravity waves at various heights. We find considerable differences between the various models. The height dependence of wave power, in particular of high-frequency waves, varies by up to two orders of magnitude between the models, and the phase difference spectra of several models show unexpected features, including $\pm180^\circ$ phase jumps.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Influence of ambipolar and Hall effects on vorticity in 3D simulations of magneto-convection

This paper presents the results of the analysis of 3D simulations of solar magneto-convection that include the joint action of the ambipolar diffusion and the Hall effect. Three simulation-runs are compared: one including both ambipolar diffusion and Hall effect; one including only ambipolar diffusion; and one without any of these two effects. The magnetic field is amplified from initial field to saturation level by the action of turbulent local dynamo. In each of these cases, we study 2 hours of simulated solar time after the local dynamo reaches the saturation regime. We analyze the power spectra of vorticity, of magnetic field fluctuations and of the different components of the magnetic Poynting flux responsible for the transport of vertical or horizontal perturbations. Our preliminary results show that the ambipolar diffusion produces a strong reduction of vorticity in the upper chromospheric layers and that it dissipates the vortical perturbations converting them into thermal energy. The Hall effect acts in the opposite way, strongly enhancing the vorticity. When the Hall effect is included, the magnetic field in the simulations becomes, on average, more vertical and long-lived flux tube-like structures are produced. We trace a single magnetic structure to study its evolution pattern and the magnetic field intensification, and their possible relation to the Hall effect.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Two-fluid simulations of Rayleigh-Taylor instability in a magnetized solar prominence thread. I. Effects of prominence magnetization and mass loading

Solar prominences are formed by partially ionized plasma with inter-particle collision frequencies generally warranting magnetohydrodynamic treatment. In this work, we explore the dynamical impacts and observable signatures of two-fluid effects in the parameter regimes when ion-neutral collisions do not fully couple the neutral and charged fluids. We perform 2.5D two-fluid (charges - neutrals) simulations of the Rayleigh-Taylor instability (RTI) at a smoothly changing interface between a solar prominence thread and the corona. The purpose of this study is to deepen our understanding of the RTI and the effects of the partial ionization on the development of RTI using non-linear two-fluid numerical simulations. Our two-fluid model takes into account viscosity, thermal conductivity, and collisional interaction between neutrals and charges: ionization/recombination, energy and momentum transfer, and frictional heating. In this paper I, the sensitivity of the RTI dynamics to the prominence equilibrium configuration, including the impact of the magnetic field strength and shear supporting the prominence thread, and the amount of prominence mass-loading is explored. We show that, at small scales, a realistically smooth prominence-corona interface leads to qualitatively different linear RTI evolution than that expected for a discontinuous interface, while magnetic field shear has the stabilizing effect of reducing the growth rate or eliminating the instability. In the non-linear phase, we observe that in the presence of field shear the development of the instability leads to formation of coherent and interacting 2.5D magnetic structures, which, in turn, can lead to substantial plasma flow across magnetic field lines and associated decoupling of the fluid velocities of charges and neutrals.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An introductory guide to fluid models with anisotropic temperatures Part 1 -- CGL description and collisionless fluid hierarchy

We present a detailed guide to advanced collisionless fluid models that incorporate kinetic effects into the fluid framework, and that are much closer to the collisionless kinetic description than traditional magnetohydrodynamics. Such fluid models are directly applicable to modeling turbulent evolution of a vast array of astrophysical plasmas, such as the solar corona and the solar wind, the interstellar medium, as well as accretion disks and galaxy clusters. The text can be viewed as a detailed guide to Landau fluid models and it is divided into two parts. Part 1 is dedicated to fluid models that are obtained by closing the fluid hierarchy with simple (non Landau fluid) closures. Part 2 is dedicated to Landau fluid closures. Here in Part 1, we discuss the CGL fluid model in great detail, together with fluid models that contain dispersive effects introduced by the Hall term and by the finite Larmor radius (FLR) corrections to the pressure tensor. We consider dispersive effects introduced by the non-gyrotropic heat flux vectors. We investigate the parallel and oblique firehose instability, and show that the non-gyrotropic heat flux strongly influences the maximum growth rate of these instabilities. Furthermore, we discuss fluid models that contain evolution equations for the gyrotropic heat flux fluctuations and that are closed at the 4th-moment level by prescribing a specific form for the distribution function. For the bi-Maxwellian distribution, such a closure is known as the "normal" closure. We also discuss a fluid closure for the bi-kappa distribution. Finally, by considering one-dimensional Maxwellian fluid closures at higher-order moments, we show that such fluid models are always unstable. The last possible non Landau fluid closure is therefore the "normal" closure, and beyond the 4th-order moment, Landau fluid closures are required.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An introductory guide to fluid models with anisotropic temperatures Part 2 -- Kinetic theory, Padé approximants and Landau fluid closures

In Part 2 of our guide to collisionless fluid models, we concentrate on Landau fluid closures. These closures were pioneered by Hammett and Perkins and allow for the rigorous incorporation of collisionless Landau damping into a fluid framework. It is Landau damping that sharply separates traditional fluid models and collisionless kinetic theory, and is the main reason why the usual fluid models do not converge to the kinetic description, even in the long-wavelength low-frequency limit. We start with a brief introduction to kinetic theory, where we discuss in detail the plasma dispersion function $Z(ζ)$, and the associated plasma response function $R(ζ)=1+ζZ(ζ)=-Z'(ζ)/2$. We then consider a 1D (electrostatic) geometry and make a significant effort to map all possible Landau fluid closures that can be constructed at the 4th-order moment level. These closures for parallel moments have general validity from the largest astrophysical scales down to the Debye length, and we verify their validity by considering examples of the (proton and electron) Landau damping of the ion-acoustic mode, and the electron Landau damping of the Langmuir mode. We proceed by considering 1D closures at higher-order moments than the 4th-order, and as was concluded in Part 1, this is not possible without Landau fluid closures. We show that it is possible to reproduce linear Landau damping in the fluid framework to any desired precision, thus showing the convergence of the fluid and collisionless kinetic descriptions. We then consider a 3D (electromagnetic) geometry in the gyrotropic (long-wavelength low-frequency) limit and map all closures that are available at the 4th-order moment level. In the Appendix A, we provide comprehensive tables with Padé approximants of $R(ζ)$ up to the 8th-pole order, with many given in an analytic form.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Heating of the partially ionized solar chromosphere by waves in magnetic structures

In this paper, we show a "proof of concept" of the heating mechanism of the solar chromosphere due to wave dissipation caused by the effects of partial ionization. Numerical modeling of non-linear wave propagation in a magnetic flux tube, embedded in the solar atmosphere, is performed by solving a system of single-fluid quasi-MHD equations, which take into account the ambipolar term from the generalized Ohm's law. It is shown that perturbations caused by magnetic waves can be effectively dissipated due to ambipolar diffusion. The energy input by this mechanism is continuous and shown to be more efficient than dissipation of static currents, ultimately leading to chromospheric temperature increase in magnetic structures.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the origin of facular brightness

This paper studies the dependence of the CaIIH line core brightness on the strength and inclination of photospheric magnetic field, and on the parameters of convective and wave motions in a facular region at the solar disc center. We use three simultaneous datasets obtained at the German Vacuum Tower Telescope (Observatorio del Teide, Tenerife): (1) spectra of BaII 4554 A line registered with the instrument TESOS to measure the variations of intensity and velocity through the photosphere up to the temperature minimum; (2) spectropolarimetric data in FeI 1.56 $μ$m lines (registered with the instrument TIP II) to measure photospheric magnetic fields; (3) filtergrams in CaIIH that give information about brightness fluctuations in the chromosphere. The results show that the CaIIH brightness in the facula strongly depends on the power of waves with periods in the 5-min range, that propagate upwards, and also on the phase shift between velocity oscillations at the bottom photosphere and around the temperature minimum height, measured from BaII line. The CaIIH brightness is maximum at locations where the phase shift between temperature and velocity oscillations lies within 0-100 degrees. There is an indirect influence of convective motions on the CaIIH brightness. Namely, the higher is the amplitude of convective velocities and the larger is the height where they change their direction of motion, the brighter is the facula. Altogether, our results lead to conclusions that facular regions appear bright not only because of the Wilson depression in magnetic structures, but also due to real heating.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Three-dimensional structure of a sunspot light bridge

Active regions are the most prominent manifestations of solar magnetic fields; their generation and dissipation are fundamental problems in solar physics. Light bridges are commonly present during sunspot decay, but a comprehensive picture of their role in the removal of photospheric magnetic field is still missing. We study the three dimensional configuration of a sunspot and in particular its light bridge during one of the last stages of its decay. We present the magnetic and thermodynamical stratification inferred from full Stokes inversions of the photospheric Si I 10827 Å and Ca I 10839 Å lines obtained with the GREGOR Infrared Spectrograph of the GREGOR telescope at Observatorio del Teide, Tenerife, Spain. The analysis is complemented by a study of continuum images covering the disk passage of the active region, which are provided by the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory. The sunspot shows a light bridge with penumbral continuum intensity that separates the central umbra from a smaller umbra. We find that in this region the magnetic field lines form a canopy with lower magnetic field strength in the inner part. The photospheric light bridge is dominated by gas pressure (high-$β$), as opposed to the surrounding umbra where the magnetic pressure is higher. A convective flow is observed in the light bridge. This flow is able to bend the magnetic field lines and to produce field reversals. The field lines close above the light bridge and become as vertical and strong as in the surrounding umbra. We conclude that it develops because of two highly magnetized regions which come closer during the sunspot evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Evershed flow observed in neutral and singly ionized iron lines

The amplitudes of the Evershed flow are measured using pairs of carefully selected FeI and FeII spectral lines located close in wavelength and registered simultaneously. A sunspot belonging to the NOAA 11582 group was scanned using the spectrograph of the German Vacuum Tower Telescope (Observatorio del Teide, Tenerife). Velocities were extracted from intensity profiles using the lambda-meter technique. The formation heights of the observed spectral lines were calculated using semi-empirical models of a bright and dark penumbral filament taking into account the sunspot location at the limb. Our objective is to compare azimuthally averaged amplitudes of the Evershed flow extracted from neutral and ion lines. We find measurable differences in the radial component of the flow. All five pairs of lines show the same tendency, with a few hundred m/s larger amplitude of the flow measured from FeI lines compared to FeII lines. This tendency is preserved at all photospheric heights and radial distances in the penumbra. We discuss the possible origin of this effect.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the nature of transverse coronal waves revealed by wavefront dislocations

Coronal waves are an important aspect of the dynamics of the plasma in the corona. Wavefront dislocations are topological features of most waves in nature and also of magnetohydrodynamic waves. Are there dislocations in coronal waves? The finding and explanation of dislocations may shed light on the nature and characteristics of the propagating waves, their interaction in the corona and in general on the plasma dynamics. We positively identify dislocations in coronal waves observed by the Coronal Multi-channel Polarimeter (CoMP) as singularities in the Doppler shifts of emission coronal lines. We study the possible singularities that can be expected in coronal waves and try to reproduce the observed dislocations in terms of localization and frequency of appearance. The observed dislocations can only be explained by the interference of a kink and a sausage wave modes propagating with different frequencies along the coronal magnetic field. In the plane transverse to the propagation, the cross-section of the oscillating plasma must be smaller than the spatial resolution, and the two waves result in net longitudinal and transverse velocity components that are mixed through projection onto the line of sight. Alfvén waves can be responsible of the kink mode, but a magnetoacoustic sausage mode is necessary in all cases. Higher (flute) modes are excluded. The kink mode has a pressure amplitude that is smaller than the pressure amplitude of the sausage mode, though its observed velocity is larger. This concentrates dislocations on the top of the loop. To explain dislocations, any model of coronal waves must include the simultaneous propagation and interference of kink and sausage wave modes of comparable but different frequencies, with a sausage wave amplitude much smaller than the kink one.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the robustness of the pendulum model for large-amplitude longitudinal oscillations in prominences

Large-amplitude longitudinal oscillations (LALOs) in prominences are spectacular manifestations of the solar activity. In such events nearby energetic disturbances induce periodic motions on filaments with displacements comparable to the size of the filaments themselves and with velocities larger than 20 km/s. The pendulum model, in which the gravity projected along a rigid magnetic field is the restoring force, was proposed to explain these events. However, it can be objected that in a realistic situation where the magnetic field reacts to the mass motion of the heavy prominence, the simplified pendulum model could be no longer valid. We have performed non-linear time-dependent numerical simulations of LALOs considering a dipped magnetic field line structure. In this work we demonstrate that for even relatively weak magnetic fields the pendulum model works very well. We therefore validate the pendulum model and show its robustness, with important implications for prominence seismology purposes. With this model it is possible to infer the geometry of the dipped field lines that support the prominence.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rayleigh-Taylor instability in partially ionized compressible plasmas: one fluid approach

We study the modification of the classical criterion for the linear onset and growth rate of the Rayleigh-Taylor instability (RTI) in a partially ionized (PI) plasma in the one-fluid description, considering a generalized induction equation. The governing linear equations and appropriate boundary conditions, including gravitational terms, are derived and applied to the case of the RTI in a single interface between two partially ionized plasmas. The boundary conditions lead to an equation for the frequencies in which some of them have positive complex parts, marking the appearance of the RTI. We study the ambipolar term alone first, extending the result to the full induction equation later. We find that the configuration is always unstable because of the presence of a neutral species. In the classical stability regime the growth rate is small, since the collisions prevent the neutral fluid to fully develop the RTI. For parameters in the classical instability regime the growth rate is lowered, but for the considered theoretical values of the collision frequencies and diffusion coefficients for solar prominences the differences with the compressible MHD case are small. We conclude that PI modifies some aspects of the linear RTI instability, since it takes into account that neutrals do not feel the stabilizing effect of the magnetic field. For the set of parameters representative for solar prominences, our model gives the resulting timescale comparable with observed lifetimes of RTI plumes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rayleigh-Taylor instability in prominences from numerical simulations including partial ionization effects

We study the Rayleigh-Taylor instability (RTI) at a prominence-corona transition region in a non-linear regime. Our aim is to understand how the presence of neutral atoms in the prominence plasma influences the instability growth rate, and the evolution of velocity, magnetic field vector and thermodynamic parameters of turbulent drops. We perform 2.5D numerical simulations of the instability initiated by a multi-mode perturbation at the corona-prominence interface using a single-fluid MHD approach including a generalized Ohm's law. The initial equilibrium configuration is purely hydrostatic and contains a homogeneous horizontal magnetic field forming an angle with the direction in which the plasma is perturbed. We analyze simulations with two different orientations of the magnetic field. For each field orientation we compare two simulations, one for the pure MHD case, and one including the ambipolar diffusion in the Ohm's law (AD case). Other than that, both simulations for each field orientation are identical. The numerical results in the initial stage of the instability are compared with the analytical linear calculations. We find that the configuration is always unstable in the AD case. The growth rate of the small-scale modes in the non-linear regime is up to 50% larger in the AD case than in the purely MHD case and the average velocities of flows are a few percent larger. Significant drift momenta are found at the interface between the coronal and the prominence material at all stages of the instability, produced by the faster downward motion of the neutral component with respect to the ionized component. The differences in temperature of the bubbles between the ideal and non-ideal case are also significant, reaching 30%. There is an asymmetry between large rising bubbles and small-scale down flowing fingers, favoring the detection of upward velocities in observations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Synthetic observations of wave propagation in a sunspot umbra

Spectropolarimetric temporal series from Fe I $λ$ 6301.5 Å and Ca II infrared triplet lines are obtained by applying the Stokes synthesis code NICOLE to a numerical simulation of wave propagation in a sunspot umbra from MANCHA code. The analysis of the phase difference between Doppler velocity and intensity core oscillations of the Fe I $λ$ 6301.5 Å line reveals that variations in the intensity are produced by opacity fluctuations rather than intrinsic temperature oscillations, except for frequencies between 5 and 6.5 mHz. On the other hand, the photospheric magnetic field retrieved from the weak field approximation provides the intrinsic magnetic field oscillations associated to wave propagation. Our results suggest that this is due to the low magnetic field gradient of our sunspot model. The Stokes parameters of the chromospheric Ca II infrared triplet lines show striking variations as shock waves travel through the formation height of the lines, including emission self-reversals in the line core and highly abnormal Stokes V profiles. Magnetic field oscillations inferred from the Ca II infrared lines using the weak field approximation appear to be related with the magnetic field strength variation between the photosphere and the chromosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dislocations in magnetohydrodynamic waves in a stellar atmosphere

We describe the presence of wavefront dislocations in magnetohydrodynamic waves in stratified stellar atmospheres. Scalar dislocations such as edges and vortices can appear in Alfvén waves, as well as in general magneto-acoustic waves. We detect those dislocations in observations of magnetohydrodynamic waves in sunspots in the solar chromosphere. Through the measured charge of all the dislocations observed, we can give for the first time estimates of the modal contribution in the waves propagating along magnetic fields in solar sunspots.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Rayleigh-Taylor instability in partially ionized prominence plasma

We study Rayleigh-Taylor instability (RTI) at the coronal-prominence boundary by means of 2.5D numerical simulations in a single-fluid MHD approach including a generalized Ohm's law. The initial configuration includes a homogeneous magnetic field forming an angle with the direction in which the plasma is perturbed. For each field inclination we compare two simulations, one for the pure MHD case, and one including the ambipolar diffusion in the Ohm's law, otherwise identical. We find that the configuration containing neutral atoms is always unstable. The growth rate of the small-scale modes in the non-linear regime is larger than in the purely MHD case.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Properties of convective motions in facular regions

In this paper, we study the properties of solar granulation in a facular region from the photosphere up to the lower chromosphere. Our aim is to investigate the dependence of granular structure on magnetic field strength. We use observations obtained at the German Vacuum Tower Telescope (Observatorio del Teide, Tenerife) using two different instruments: Triple Etalon SOlar Spectrometer (TESOS), in the BaII 4554 A line to measure velocity and intensity variations along the photosphere; and, simultaneously, Tenerife Infrared Polarimeter (TIP-II), in the FeI 1.56 $μ$m lines to the measure Stokes parameters and the magnetic field strength at the lower photosphere. We obtain that the convective velocities of granules in the facular area decrease with magnetic field while the convective velocities of intergranular lanes increase with the field strength. Similar to the quiet areas, there is a contrast and velocity sign reversal taking place in the middle photosphere. The reversal heights depend on the magnetic field strength and are, on average, about 100 km higher than in the quiet regions. The correlation between convective velocity and intensity decreases with magnetic field at the bottom photosphere, but increases in the upper photosphere. The contrast of intergranular lanes observed close to the disc center is almost independent of the magnetic field strength. The strong magnetic field of facular area seems to stabilize the convection and to promote more effective energy transfer in the upper layers of the solar atmosphere, since the convective elements reach larger heights.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Solar Fe abundance and magnetic fields - Towards a consistent reference metallicity

We investigate the impact on Fe abundance determination of including magnetic flux in series of 3D radiation-MHD simulations of solar convection which we used to synthesize spectral intensity profiles corresponding to disc centre. A differential approach is used to quantify the changes in theoretical equivalent width of a set of 28 iron spectral lines spanning a wide range in lambda, excitation potential, oscillator strength, Landé factor, and formation height. The lines were computed in LTE using the spectral synthesis code LILIA. We used input magnetoconvection snapshots covering 50 minutes of solar evolution and belonging to series having an average vertical magnetic flux density of 0, 50, 100 and 200 G. For the relevant calculations we used the Copenhagen Stagger code. The presence of magnetic fields causes both a direct (Zeeman-broadening) effect on spectral lines with non-zero Landé factor and an indirect effect on temperature-sensitive lines via a change in the photospheric T-tau stratification. The corresponding correction in the estimated atomic abundance ranges from a few hundredths of a dex up to |Delta log(Fe)| ~ 0.15 dex, depending on the spectral line and on the amount of average magnetic flux within the range of values we considered. The Zeeman-broadening effect gains relatively more importance in the IR. The largest modification to previous solar abundance determinations based on visible spectral lines is instead due to the indirect effect, i.e., the line-weakening caused by a warmer stratification on an optical depth scale. Our results indicate that the average solar iron abundance obtained when using magnetoconvection models can be 0.03-0.11 dex higher than when using the simpler HD convection approach. We demonstrate that accounting for magnetic flux is important in state-of-the-art solar photospheric abundance determinations based on 3D simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Influence of phase-diversity image reconstruction techniques on circular polarization asymmetries

Full Stokes filter-polarimeters are key instruments for investigating the rapid evolution of magnetic structures on the solar surface. To this end, the image quality is routinely improved using a-posteriori image reconstruction methods. We analyze the robustness of circular polarization asymmetries to phase-diversity image reconstruction techniques. We use snapshots of magneto-hydrodynamical simulations carried out with different initial conditions to synthesize spectra of the magnetically sensitive Fe I line at 5250.2 A. We degrade the synthetic profiles spatially and spectrally to simulate observations with the IMaX full Stokes filter-polarimeter. We also simulate the focused/defocused pairs of images used by the phase-diversity algorithm for reconstruction and the polarimetric modulation scheme. We assume that standard optimization methods are able to infer the projection of the wavefront on the Zernike polynomials with 10% precision. We also consider the less favorable case of 25% precision. We obtain reconstructed monochromatic modulated images that are later demodulated and compared with the original maps. Although asymmetries are often difficult to define in the quiet Sun due to the complexity of the Stokes V profiles, we show how asymmetries are degraded with spatial and spectral smearing. The results indicate that, although image reconstruction techniques reduce the spatial smearing, they can modify the asymmetries of the profiles, mainly caused by the appearance of spatially-correlated noise.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Magneto-acoustic wave energy from numerical simulations of an observed sunspot umbra

We aim at reproducing the height dependence of sunspot wave signatures obtained from spectropolarimetric observations through 3D MHD numerical simulations. A magneto-static sunspot model based on the properties of the observed sunspot is constructed and perturbed at the photosphere introducing the fluctuations measured with the \SiI\ $λ$ 10827 Å line. The results of the simulations are compared with the oscillations observed simultaneously at different heights from the \HeI\ $λ$ 10830 Å line, the \CaIIH\ core and the \FeI\ blends in the wings of the \CaIIH\ line. The simulations show a remarkable agreement with the observations. They reproduce the velocity maps and power spectra at the formation heights of the observed lines, as well as the phase and amplification spectra between several pair of lines. We find that the stronger shocks at the chromosphere are accompanied with a delay between the observed signal and the simulated one at the corresponding height, indicating that shocks shift the formation height of the chromospheric lines to higher layers. Since the simulated wave propagation matches very well the properties of the observed one, we are able to use the numerical calculations to quantify the energy contribution of the magneto-acoustic waves to the chromospheric heating in sunspots. Our findings indicate that the energy supplied by these waves is too low to balance the chromospheric radiative losses. The energy contained at the formation height of the lowermost \SiI\ $λ$ 10827 Å line in the form of slow magneto-acoustic waves is already insufficient to heat the higher layers, and the acoustic energy which reaches the chromosphere is around 3-9 times lower than the required amount of energy. The contribution of the magnetic energy is even lower.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Unnoticed magnetic field oscillations in the very quiet Sun revealed by Sunrise/IMaX

We present observational evidence for oscillations of magnetic flux density in the quiet areas of the Sun. The majority of magnetic fields on the solar surface have strengths of the order of or lower than the equipartition field (300-500 G). This results in a myriad of magnetic fields whose evolution is largely determined by the turbulent plasma motions. When granules evolve they squash the magnetic field lines together or pull them apart. Here we report on the periodic deformation of the shapes of features in circular polarization observed at high resolution with Sunrise. In particular, we note that the area of patches with constant magnetic flux oscillates with time, which implies that the apparent magnetic field intensity oscillates in antiphase. The periods associated to this oscillatory pattern is compatible with the granular life-time and change abruptly, which suggests that these oscillations might not correspond to characteristic oscillatory modes of magnetic structures, but to the forcing by granular motions. In one particular case, we find three patches around the same granule oscillating in phase, which means that the spatial coherence of these oscillations can reach 1600 km. Interestingly, the same kind of oscillatory phenomenon is found also in the upper photosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Magneto-acoustic waves in sunspots from observations and numerical simulations

We study the propagation of waves from the photosphere to the chromosphere of sunspots. From time series of cospatial Ca II H (including its line blends) intensity spectra and polarimetric spectra of Si I 1082.7 nm and He I 1083.0 nm we retrieve the line-of-sight velocity at several heights. The analysis of the phase difference and amplification spectra shows standing waves for frequencies below 4 mHz and propagating waves for higher frequencies, and allows us to infer the temperature and height where the lines are formed. Using these observational data, we have constructed a model of sunspot, and we have introduced the velocity measured with the photospheric Si I 1082.7 nm line as a driver. The numerically propagated wave pattern fits reasonably well with the observed using the lines formed at higher layers, and the simulations reproduce many of the observed features. The observed waves are slow MHD waves propagating longitudinally along field lines.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Magneto-acoustic waves in sunspots: first results from a new 3D nonlinear magnetohydrodynamic code

Waves observed in the photosphere and chromosphere of sunspots show complex dynamics and spatial patterns. The interpretation of high-resolution sunspot wave observations requires modeling of three-dimensional non-linear wave propagation and mode transformation in the sunspot upper layers in realistic spot model atmospheres. Here we present the first results of such modeling. We have developed a 3D non-linear numerical code specially designed to calculate the response of magnetic structures in equilibrium to an arbitrary perturbation. The code solves the 3D nonlinear MHD equations for perturbations; it is stabilized by hyper-diffusivity terms and is fully parallelized. The robustness of the code is demonstrated by a number of standard tests. We analyze several simulations of a sunspot perturbed by pulses of different periods at subphotospheric level, from short periods, introduced for academic purposes, to longer and realistic periods of three and five minutes. We present a detailed description of the three-dimensional mode transformation in a non-trivial sunspot-like magnetic field configuration, including the conversion between fast and slow magneto-acoustic waves and the Alfvén wave, by calculation of the wave energy fluxes. Our main findings are the following: (1) the conversion from acoustic to the Alfvén mode is only observed if the the driving pulse is located out of the sunspot axis, but this conversion is energetically inefficient; (2) as a consequence of the cut-off effects and refraction of the fast magneto-acoustic mode, the energy of the evanescent waves with periods around 5 minutes remains almost completely below the level beta=1; (3) waves with frequencies above the cut-off propagate field-aligned to the chromosphere and their power becomes dominating over that of evanescent 5-minute oscillations, in agreement with observations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Mode transformation and frequency change with height in 3D numerical simulations of magneto-acoustic wave propagation in sunspots

Three-dimensional numerical simulations of magnetoacoustic wave propagation are performed in a sunspot atmosphere with a computational domain covering from the photosphere to the chromosphere. The wave source, with properties resembling the solar spectrum, is located at different distances from the axis of the sunspot for each simulation. These results are compared with the theory of mode transformation and also with observational features. Simulations show that the dominant oscillation frequency in the chromosphere decreases with the radial distance from the sunspot axis. The energy flux of the different wave modes involved, including de Alfvén mode, is evaluated and discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Modeling the Subsurface Structure of Sunspots

While sunspots are easily observed at the solar surface, determining their subsurface structure is not trivial. There are two main hypotheses for the subsurface structure of sunspots: the monolithic model and the cluster model. Local helioseismology is the only means by which we can investigate subphotospheric structure. However, as current linear inversion techniques do not yet allow helioseismology to probe the internal structure with sufficient confidence to distinguish between the monolith and cluster models, the development of physically realistic sunspot models are a priority for helioseismologists. This is because they are not only important indicators of the variety of physical effects that may influence helioseismic inferences in active regions, but they also enable detailed assessments of the validity of helioseismic interpretations through numerical forward modeling. In this paper, we provide a critical review of the existing sunspot models and an overview of numerical methods employed to model wave propagation through model sunspots. We then carry out an helioseismic analysis of the sunspot in Active Region 9787 and address the serious inconsistencies uncovered by \citeauthor{gizonetal2009}~(\citeyear{gizonetal2009,gizonetal2009a}). We find that this sunspot is most probably associated with a shallow, positive wave-speed perturbation (unlike the traditional two-layer model) and that travel-time measurements are consistent with a horizontal outflow in the surrounding moat.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Multi-layer study of wave propagation in sunspots

We analyze the propagation of waves in sunspots from the photosphere to the chromosphere using time series of co-spatial Ca II H intensity spectra (including its line blends) and polarimetric spectra of Si I 10827 and the He I 10830 multiplet. From the Doppler shifts of these lines we retrieve the variation of the velocity along the line-of-sight at several heights. Phase spectra are used to obtain the relation between the oscillatory signals. Our analysis reveals standing waves at frequencies lower than 4 mHz and a continuous propagation of waves at higher frequencies, which steepen into shocks in the chromosphere when approaching the formation height of the Ca II H core. The observed non-linearities are weaker in Ca II H than in He I lines. Our analysis suggests that the Ca II H core forms at a lower height than the He I 10830 line: a time delay of about 20 s is measured between the Doppler signal detected at both wavelengths. We fit a model of linear slow magnetoacoustic wave propagation in a stratified atmosphere with radiative losses according to Newton's cooling law to the phase spectra and derive the difference in the formation height of the spectral lines. We show that the linear model describes well the wave propagation up to the formation height of Ca II H, where non-linearities start to become very important.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerical Simulation of Excitation and Propagation of Helioseismic MHD Waves in Magnetostatic Models of Sunspots

We present comparison of numerical simulations of propagation of MHD waves,excited by subphotospheric perturbations, in two different ("deep" and "shallow") magnetostatic models of the sunspots. The "deep" sunspot model distorts both the shape of the wavefront and its amplitude stronger than the "shallow" model. For both sunspot models, the surface gravity waves (f-mode) are affected by the sunspots stronger than the acoustic p-modes. The wave amplitude inside the sunspot depends on the photospheric strength of the magnetic field and the distance of the source from the sunspot axis. For the source located at 9 Mm from the center of the sunspot, the wave amplitude increases when the wavefront passes through the central part of the sunspot. For the source distance of 12 Mm, the wave amplitude inside the sunspot is always smaller than outside. For the same source distance from the sunspot center but for the models with different strength of the magnetic field, the wave amplitude inside the sunspot increases with the strength of the magnetic field. The simulations show that unlike the case of the uniform inclined background magnetic field, the p- and f-mode waves are not spatially separated inside the sunspot where the magnetic field is strongly non-uniform. These properties have to be taken into account for interpretation of observations of MHD waves traveling through sunspot regions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Solar abundance corrections derived through 3D magnetoconvection simulations

We explore the effect of the magnetic field when using realistic three-dimensional convection experiments to determine solar element abundances. By carrying out magnetoconvection simulations with a radiation-hydro code (the Copenhagen stagger code) and through a-posteriori spectral synthesis of three Fe I lines, we obtain evidence that moderate amounts of mean magnetic flux cause a noticeable change in the derived equivalent widths compared with those for a non-magnetic case. The corresponding Fe abundance correction for a mean flux density of 200 G reaches up to ~0.1 dex in magnitude. These results are based on space- and time-averaged line profiles over a time span of 2.5 solar hours in the statistically stationary regime of the convection. The main factors causing the change in equivalent widths, namely the Zeeman broadening and the modification of the temperature stratification, act in different amounts and, for the iron lines considered here, in opposite directions; yet, the resulting $|Δ\logε_{\odot}(Fe)|$ coincides within a factor two in all of them, even though the sign of the total abundance correction is different for the visible and infrared lines. We conclude that magnetic effects should be taken into account when discussing precise values of the solar and stellar abundances and that an extended study is warranted.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Towards pulsation mode identification in 3-D: theoretical simulations of line profile variations in roAp stars

Time-resolved spectroscopic observations of rapidly oscillating Ap (roAp) stars show a complex picture of propagating magneto-acoustic pulsation waves, with amplitude and phase strongly changing as a function of atmospheric height. We have recently conducted numerical, non-linear MHD simulations to get an insight into the complex atmospheric dynamics of magnetic pulsators. Here we use the resulting time-dependent atmospheric structure and velocity field to predict line profile variations for roAp stars. These calculations use realistic atmospheric structure, account for vertical chemical stratification and treat the line formation in pulsating stellar atmosphere without relying on the simplistic single-layer approximation universally adopted for non-radial pulsators. The new theoretical calculations provide an essential tool for interpreting the puzzling complexity of the spectroscopic pulsations in roAp stars.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Where the granular flows bend

Based on IMaX/Sunrise data, we report on a previously undetected phenomenon in solar granulation. We show that in a very narrow region separating granules and intergranular lanes the spectral line width of the Fe I 5250.2 A line becomes extremely small. We offer an explanation of this observation with the help of magneto-convection simulations. These regions with extremely small line widths correspond to the places where the granular flows bend from mainly upflow in granules to downflow in intergranular lanes. We show that the resolution and image stability achieved by IMaX/Sunrise are important requisites to detect this interesting phenomenon.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Simulations of magneto-acoustic pulsations in atmospheres of rapidly oscillating Ap stars

Rapidly oscillating Ap stars exhibit an astrophysically interesting combination of strong, dipolar-like magnetic fields and high-overtone p-mode pulsations similar to the Sun. Recent time-resolved spectroscopy of these stars unravelled a complex picture of propagating magneto-acoustic pulsation waves, with amplitude and phase strongly changing as a function of atmospheric height. To interpret these observations and gain a new insight into the atmospheric dynamics of roAp stars we have carried out 2-D time-dependent, non-linear magneto-hydrodynamical simulations of waves for a realistic atmospheric stratification of a cool Ap star. We explore a grid of simulations in a wide parameter space, treating oscillations of the velocity, magnetic field and thermodynamic quantities in a self-consistent manner. Our simulations foster a new understanding of the influence of the atmosphere and the magnetic field on the propagation and reflection properties of magneto-acoustic waves, formation of node surfaces, and relative variation of different quantities. Our simulations reproduce all main features of the observed pulsational behavior of roAp stars. We show, for the first time, that the overall depth dependence of the pulsations in roAp atmospheres is strongly influenced by the density inversion at the photospheric base.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Solar BaII 4554 A line as Doppler diagnostics: NLTE analysis in 3D hydrodynamical model

The aim of this paper is to analyse the validity of the Dopplergram and lambda-meter techniques for the Doppler diagnostics of solar photospheric velocities using the BaII 4554 A line. Both techniques are evaluated by means of NLTE radiative transfer calculations of the BaII 4554 A line in a three-dimensional hydrodynamical model of solar convection. We consider the cases of spatially unsmeared profiles and the profiles smeared to the resolution of ground-based observations. We find that: (i) Speckle-reconstructed Dopplergram velocities reproduce the ``true'' velocities well at heights around 300 km, except for intergranular lanes with strong downflows where the velocity can be overestimated. (ii) The lambda-meter velocities give a good representation of the ``true'' velocities through the whole photosphere, both under the original and reduced spatial resolutions. The velocities derived from the inner wing of smeared BaII 4554 A line profiles are more reliable than those for the outer wing. Only under high spatial resolution does the inner wing velocities calculated in intergranular regions give an underestimate (or even a sign reversal) compared with the model velocities. (iii) NLTE effects should be taken into account in modelling the BaII 4554 A line profiles. Such effects are more pronounced in intergranular regions. Our analysis supports the opinion that the Dopplergram technique applied to the BaII 4554 A line is a valuable tool for the Doppler diagnostics of the middle photosphere around 300 km. The λ-meter technique applied to this line gives us a good opportunity to ``trace'' the non-thermal motions along the whole photosphere up to the temperature minimum and lower chromosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Solar granulation from photosphere to low chromosphere observed in BaII 4554 A line

The purpose of this paper is to characterize the statistical properties of solar granulation in the photosphere and low chromosphere up to 650 km. We use velocity and intensity variations obtained at different atmospheric heights from observations in BaII 4554 A. The observations were done during good seeing conditions at the VTT at the Observatorio del Teide on Tenerife. The line core forms rather high in the atmosphere and allows granulation properties to be studied at heights that have been not accessed before in similar studies. In addition, we analyze the synthetic profiles of the BaII 4554 A line by the same method computed taking NLTE effects into account in the 3D hydrodynamical model atmosphere. We suggest a 16-column model of solar granulation depending on the direction of motion and on the intensity contrast measured in the continuum and in the uppermost layer. We calculate the heights of intensity contrast sign reversal and velocity sign reversal. We show that both parameters depend strongly on the granulation velocity and intensity at the bottom photosphere. The larger the two parameters, the higher the reversal takes place in the atmosphere. On average, this happens at about 200-300 km. We suggest that this number also depends on the line depth of the spectral line used in observations. Despite the intensity and velocity reversal, about 40% of the column structure of granulation is preserved up to heights around 650 km.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Sunspot seismic halos generated by fast MHD wave refraction

We suggest an explanation for the high-frequency power excess surrounding active regions known as seismic halos. The idea is based on numerical simulations of magneto-acoustic waves propagation in sunspots. We propose that such an excess can be caused by the additional energy injected by fast mode waves refracted in the higher atmosphere due to the rapid increase of the Alfven speed. Our model qualitatively explains the magnitude of the halo and allows to make some predictions of its behavior that can be checked in future observations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The energy of waves in the photosphere and lower chromosphere: 1. Velocity statistics

Acoustic waves are one of the primary suspects besides magnetic fields for the chromospheric heating process to temperatures above radiative equilibrium (RE). We derived the mechanical wave energy as seen in line-core velocities to obtain a measure of mechanical energy flux with height for a comparison with the energy requirements in a semi-empirical atmosphere model. We analyzed a 1-hour time series and a large-area map of Ca II H spectra on the traces of propagating waves. We analyzed the velocity statistics of several spectral lines in the wing of Ca II H, and the line-core velocity of Ca II H. We converted the velocity amplitudes into volume and mass energy densities. For comparison, we used the increase of internal energy necessary to lift a RE atmosphere to the HSRA temperature stratification. We find that the velocity amplitude grows in agreement with linear wave theory and thus slower with height than predicted from energy conservation. The mechanical energy of the waves above around z~500 km is insufficient to maintain the chromospheric temperature rise in the semi-empirical HSRA model. The intensity variations of the Ca line core (z~1000 km) can be traced back to the velocity variations of the lowermost forming spectral line considered (z~ 250 km). The chromospheric intensity, and hence, (radiation) temperature variations are seen to be induced by passing waves originating in the photosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Theoretical modeling of propagation of magneto-acoustic waves in magnetic regions below sunspots

We use 2D numerical simulations and eikonal approximation, to study properties of MHD waves traveling below the solar surface through the magnetic structure of sunspots. We consider a series of magnetostatic models of sunspots of different magnetic field strengths, from 10 Mm below the photosphere to the low chromosphere. The purpose of these studies is to quantify the effect of the magnetic field on local helioseismology measurements by modeling waves excited by sub-photospheric sources. Time-distance propagation diagrams and wave travel times are calculated for models of various field strength and compared to the non-magnetic case. The results clearly indicate that the observed time-distance helioseismology signals in sunspot regions correspond to fast MHD waves. The slow MHD waves form a distinctly different pattern in the time-distance diagram, which has not been detected in observations. The numerical results are in good agreement with the solution in the short-wavelength (eikonal) approximation, providing its validation. The frequency dependence of the travel times is in a good qualitative agreement with observations.

0 citations0 reviewsNo code linkedOpen access