Source author record

L. S. Anusha

L. S. Anusha 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

16works
3topics
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

16 published item(s)

preprint2022arXiv

Chromospheric Extension of the MURaM Code

Detailed numerical models of chromosphere and corona are required to understand the heating of the solar atmosphere. An accurate treatment of the solar chromosphere is complicated by the effects arising from Non Local Thermodynamic Equilibrium (NLTE) radiative transfer. A small number of strong, highly scattering lines dominate the cooling and heating in the chromosphere. Additionally, the recombination times of ionised hydrogen are longer than the dynamical timescales, requiring a non-equilibrium (NE) treatment of hydrogen ionisation. The MURaM code is extended to include the physical process required for accurate simulation of the solar chromosphere, as implemented in the Bifrost code. This includes a time-dependent treatment of hydrogen ionisation, a scattering multi-group radiation transfer scheme and approximations for NLTE radiative cooling. The inclusion of NE and NLTE physics has a large impact on the structure of the chromosphere; the NE treatment of hydrogen ionisation leads to a higher ionisation fraction and enhanced populations in the first excited state throughout cold inter-shock regions of the chromosphere. Additionally this prevents hydrogen ioniation from buffering energy fluctuations, leading to hotter shocks and cooler inter-shock regions. The hydrogen populations in the ground and first excited state are enhanced by $10^2-10^3$ in the upper chromosphere and up to $10^9$ near the transition region. Including the necessary NLTE physics leads to significant differences in chromospheric structure and dynamics. The thermodynamics and hydrogen populations calculated using the extended version of the MURaM code are consistent with previous non-equilibrium simulations. The electron number and temperature calculated using the non-equilibrium treatment of the chromosphere are required to accurately synthesise chromospheric spectral lines.

preprint2021arXiv

The Dimmest State of the Sun

How the solar electromagnetic energy entering the Earth's atmosphere varied since pre-industrial times is an important consideration in the climate change debate. Detrimental to this debate, estimates of the change in total solar irradiance (TSI) since the Maunder minimum, an extended period of weak solar activity preceding the industrial revolution, differ markedly, ranging from a drop of 0.75 Wm-2 to a rise of 6.3 Wm-2. Consequently, the exact contribution by solar forcing to the rise in global temperatures over the past centuries remains inconclusive. Adopting a novel approach based on state-of-the-art solar imagery and numerical simulations, we establish the TSI level of the Sun when it is in its least-active state to be 2.0 +/- 0.7 Wm-2 below the 2019 level. This means TSI could not have risen since the Maunder minimum by more than this amount, thus restricting the possible role of solar forcing in global warming.

preprint2020arXiv

Importance of Angle-dependent Partial Frequency Redistribution in Hyperfine Structure Transitions Under Incomplete Paschen-Back Effect Regime

Angle-frequency coupling in scattering of polarized light on atoms is represented by the angle-dependent (AD) partial frequency redistribution (PRD) matrices. There are several lines in the linearly polarized solar spectrum, for which PRD combined with quantum interference between hyperfine structure states play a significant role. Here we present the solution of the polarized line transfer equation including the AD-PRD matrix for scattering on a two-level atom with hyperfine structure splitting (HFS) and an unpolarized lower level. We account for the effects of arbitrary magnetic fields (including the incomplete Paschen-Back effect regime) and elastic collisions. For exploratory purposes we consider a self-emitting isothermal planar atmosphere and use atomic parameters that represent an isolated Na\,{\sc i} D$_2$ line. For this case we show that the AD-PRD effects are significant for field strengths below about 30G, but that the computationally much less demanding approximation of angle-averaged (AA) PRD may be used for stronger fields.

preprint2020arXiv

Power spectrum of turbulent convection in the solar photosphere

The solar photosphere provides us with a laboratory for understanding turbulence in a layer where the fundamental processes of transport vary rapidly and a strongly superadiabatic region lies very closely to a subadiabatic layer. Our tools for probing the turbulence are high-resolution spectropolarimetric observations such as have recently been obtained with the two sunrise missions, and numerical simulations. Our aim is to study photospheric turbulence with the help of Fourier power spectra that we compute from observations and simulations. We also attempt to explain some properties of the photospheric overshooting flow with the help of its governing equations and simulations. We find that quiet-Sun observations and smeared simulations exhibit a power-law behavior in the subgranular range of their Doppler velocity power spectra with an index of$~\approx -2$. The unsmeared simulations exhibit a power-law index of$~\approx -2.25$. The smearing considerably reduces the extent of the power-law-like portion of the spectra. Therefore, the limited spatial resolution in some observations might eventually result in larger uncertainties in the estimation of the power-law indices. The simulated vertical velocity power spectra as a function of height show a rapid change in the power-law index from the solar surface to $300$~km above it. A scale-dependent transport of the vertical momentum occurs. At smaller scales, the vertical momentum is more efficiently transported sideways than at larger scales. This results in less vertical velocity power transported upward at small scales than at larger scales and produces a progressively steeper vertical velocity power law below $180$ km. Above this height, the gravity work progressively gains importance at all scales, making the atmosphere progressively more hydrostatic and resulting in a gradually less steep power law.

preprint2020arXiv

Resonance Line Polarization in Spherically Symmetric Moving Media: a Parametric Study

In the present paper we consider the problem of resonance line polarization formed in the spherically symmetric expanding atmospheres. For the solution of the concerned polarized transfer equation we use the comoving frame formulation, and apply the Accelerated Lambda Iteration (ALI) method. We restrict ourselves to the non-relativistic regime of velocities wherein mainly Doppler shift effects are significant. For our studies, we consider the scattering on a two-level atom, including the effects of partial frequency redistribution (PFR). We present the dependence of linearly polarized profiles on different atmospheric and atomic parameters.

preprint2014arXiv

Center to limb observations and modeling of the Ca I 4227 A line

The observed center-to-limb variation (CLV) of the scattering polarization in different lines of the Second Solar Spectrum can be used to constrain the height variation of various atmospheric parameters, in particular the magnetic fields via the Hanle effect. Here we attempt to model non-magnetic CLV observations of the $Q/I$ profiles of the Ca I 4227 A line recorded with the ZIMPOL-3 at IRSOL. For modeling, we use the polarized radiative transfer with partial frequency redistribution with a number of realistic 1-D model atmospheres. We find that all the standard FAL model atmospheres, used by us, fail to simultaneously fit the observed ($I$, $Q/I$) at all the limb distances ($μ$). However, an attempt is made to find a single model which can provide a fit at least to the CLV of the observed $Q/I$ instead of a simultaneous fit to the ($I$, $Q/I$) at all $μ$. To this end we construct a new 1-D model by combining two of the standard models after modifying their temperature structures in the appropriate height ranges. This new combined model closely reproduces the observed $Q/I$ at all the $μ$, but fails to reproduce the observed rest intensity at different $μ$. Hence we find that no single 1-D model atmosphere succeeds in providing a good representation of the real Sun. This failure of 1-D models does not however cause an impediment to the magnetic field diagnostic potential of the Ca I 4227 A line. To demonstrate this we deduce the field strength at various $μ$ positions without invoking the use of radiative transfer.

preprint2014arXiv

Effect of cross-redistribution on the resonance scattering polarization of O {\sc i} line at 1302 Å\,

Oxygen is the most abundant element on the Sun after Hydrogen and Helium. The intensity spectrum of resonance lines of neutral Oxygen namely O {\sc i} (1302, 1305 and 1306 Å\,) has been studied in the literature for chromospheric diagnostics. In this paper we study the resonance scattering polarization in the O {\sc i} line at 1302 Å\, using two-dimensional radiative transfer in a composite atmosphere constructed using a two-dimensional magneto-hydrodynamical snapshot in the photosphere and columns of the one-dimensional FALC atmosphere in the chromosphere. The methods developed by us recently in a series of papers to solve multi-dimensional polarized radiative transfer have been incorporated in our new code POLY2D which we use for our analysis. We find that multi-dimensional radiative transfer including XRD effects is important in reproducing the amplitude and shape of scattering polarization signals of the O {\sc i} line at 1302 Å\,.

preprint2013arXiv

An efficient decomposition technique to solve angle-dependent Hanle scattering problems

Hanle scattering is an important diagnostic tool to study weak solar magnetic fields. Partial frequency redistribution (PRD) is necessary to interpret the linear polarization observed in strong resonance lines. Usually angle-averaged PRD functions are used to analyze linear polarization. However it is established that angle-dependent PRD functions are often necessary to interpret polarization profiles formed in the presence of weak magnetic fields. Our aim is to present an efficient decomposition technique, and the numerical method to solve the concerned angle-dependent line transfer problem. Together with the standard Stokes decomposition technique we employ Fourier expansion over the outgoing azimuth angle to express in a more convenient form, the angle-dependent PRD function for the Hanle effect. It allows the use of angle-dependent frequency domains of Bommier to solve the Hanle transfer problem. Such an approach is self-consistent and accurate compared to a recent approach where angle-averaged frequency domains were used to solve the same problem. We show that it is necessary to incorporate angle-dependent frequency domains instead of angle-averaged frequency domains to solve the Hanle transfer problem accurately, especially for the Stokes U parameter. The importance of using angle-dependent domains has been highlighted by taking the example of Hanle effect in the case of line transfer with vertical magnetic fields in a slab atmosphere. We have also studied the case of polarized line formation when micro-turbulent magnetic fields are present. The difference between angle-averaged and angle-dependent solutions is enhanced by the presence of micro-turbulent fields.

preprint2013arXiv

Multi-dimensional radiative transfer to analyze Hanle effect in Ca {\sc ii} K line at 3933 Å\,

Radiative transfer (RT) studies of the linearly polarized spectrum of the Sun (the second solar spectrum) have generally focused on the line formation, with an aim to understand the vertical structure of the solar atmosphere using one-dimensional (1D) model atmospheres. Modeling spatial structuring in the observations of the linearly polarized line profiles requires the solution of multi-dimensional (multi-D) polarized RT equation and a model solar atmosphere obtained by magneto-hydrodynamical (MHD) simulations of the solar atmosphere. Our aim in this paper is to analyze the chromospheric resonance line Ca {\sc ii} K at 3933 Å using multi-D polarized RT with Hanle effect and partial frequency redistribution in line scattering. We use an atmosphere which is constructed by a two-dimensional snapshot of the three-dimensional MHD simulations of the solar photosphere, combined with columns of an 1D atmosphere in the chromosphere. This paper represents the first application of polarized multi-D RT to explore the chromospheric lines using multi-D MHD atmospheres, with PRD as the line scattering mechanism. We find that the horizontal inhomogeneities caused by MHD in the lower layers of the atmosphere are responsible for strong spatial inhomogeneities in the wings of the linear polarization profiles, while the use of horizontally homogeneous chromosphere (FALC) produces spatially homogeneous linear polarization in the line core.

preprint2013arXiv

Polarized Line Formation in Multi-dimensional Media. IV. A Fourier Decomposition Technique To Formulate The Transfer Equation With Angle-Dependent Partial Frequency Redistribution

To explain the linear polarization observed in spatially resolved structures in the solar atmosphere, the solution of polarized radiative transfer (RT) equation in multi-dimensional (multi-D) geometries is essential. For strong resonance lines partial frequency redistribution (PRD) effects also become important. In a series of papers we have been investigating the nature of Stokes profiles formed in multi-D media including PRD in line scattering. For numerical simplicity so far we restricted our attention to the particular case of PRD functions which are averaged over all the incident and scattered directions. In this paper we formulate the polarized RT equation in multi-D media that takes into account Hanle effect with angle-dependent PRD functions. We generalize here to the multi-D case, the method of Fourier series expansion of angle-dependent PRD functions originally developed for RT in 1D geometry. We show that the Stokes source vector $\bm{S}=(S_I,S_Q,S_U)^T$ and the Stokes vector $\bm{I}=(I,Q,U)^T$ can be expanded in terms of infinite sets of components ${\tilde{\bm{\mathcal{S}}}}^{(k)}$, ${\tilde{\bm{\mathcal{I}}}}^{(k)}$ respectively, $k\in[0,+\infty)$. We show that the components ${\tilde{\bm{\mathcal{S}}}}^{(k)}$ become independent of the azimuthal angle ($φ$) of the scattered ray, whereas the components ${\tilde{\bm{\mathcal{I}}}}^{(k)}$ remain dependent on $φ$ due to the nature of RT in multi-D geometry. We also establish that ${\tilde{\bm{\mathcal{S}}}}^{(k)}$ and ${\tilde{\bm{\mathcal{I}}}}^{(k)}$ satisfy a simple transfer equation, which can be solved by any iterative method like an Approximate Lambda Iteration (ALI) or a Biconjugate-Gradient type projection method provided we truncate the Fourier series to have a finite number of terms.

preprint2013arXiv

Polarized Line Formation in Multi-dimensional Media. V. Effects of Angle-Dependent Partial Frequency Redistribution

The solution of polarized radiative transfer equation with angle-dependent (AD) partial frequency redistribution (PRD) is a challenging problem. Modeling the observed, linearly polarized strong resonance lines in the solar spectrum often requires the solution of the AD line transfer problems in one-dimensional (1D) or multi-dimensional (multi-D) geometries. The purpose of this paper is to develop an understanding of the relative importance of the AD PRD effects and the multi-D transfer effects and particularly their combined influence on the line polarization. This would help in a quantitative analysis of the second solar spectrum (the linearly polarized spectrum of the Sun). We consider both non-magnetic and magnetic media. In this paper we reduce the Stokes vector transfer equation to a simpler form using a Fourier decomposition technique for multi-D media. A fast numerical method is also devised to solve the concerned multi-D transfer problem. The numerical results are presented for a two-dimensional medium with a moderate optical thickness (effectively thin), and are computed for a collisionless frequency redistribution. We show that the AD PRD effects are significant, and can not be ignored in a quantitative fine analysis of the line polarization. These effects are accentuated by the finite dimensionality of the medium (multi-D transfer). The presence of magnetic fields (Hanle effect) modifies the impact of these two effects to a considerable extent.

preprint2012arXiv

J-state interference signatures in the Second Solar Spectrum: Modeling the Cr I triplet at 5204-5208 A

The scattering polarization in the solar spectrum is traditionally modeled with each spectral line treated separately, but this is generally inadequate for multiplets where J-state interference plays a significant role. Through simultaneous observations of all the 3 lines of a Cr I triplet, combined with realistic radiative transfer modeling of the data, we show that it is necessary to include J-state interference consistently when modeling lines with partially interacting fine structure components. Polarized line formation theory that includes J-state interference effects together with partial frequency redistribution for a two-term atom is used to model the observations. Collisional frequency redistribution is also accounted for. We show that the resonance polarization in the Cr I triplet is strongly affected by the partial frequency redistribution effects in the line core and near wing peaks. The Cr I triplet is quite sensitive to the temperature structure of the photospheric layers. Our complete frequency redistribution calculations in semi-empirical models of the solar atmosphere cannot reproduce the observed near wing polarization or the cross-over of the Stokes Q/I line polarization about the continuum polarization level that is due to the J-state interference. When however partial frequency redistribution is included, a good fit to these features can be achieved. Further, to obtain a good fit to the far wings, a small temperature enhancement of the FALF model in the photospheric layers is necessary.

preprint2011arXiv

Observations of the forward scattering Hanle effect in the Ca i 4227 Å line

Chromospheric magnetic fields are notoriously diffcult to measure. The chromospheric lines are broad, while the fields are producing a minuscule Zeeman-effect polarization. A promising diagnostic alternative is provided by the forward-scattering Hanle effect, which can be recorded in chromospheric lines such as the He i 10830 Å and the Ca i 4227 Å lines. We present a set of spectropolarimetric observations of the full Stokes vector obtained near the center of the solar disk in the Ca i 4227 Å line with the ZIMPOL polarimeter at the IRSOL observatory.We detect a number of interesting forward-scattering Hanle effect signatures, which we model successfully using polarized radiative transfer. Here we focus on the observational aspects, while a separate companion paper deals with the theoretical modeling.

preprint2011arXiv

Polarized Line Formation in Multi-Dimensional Media.III. Hanle Effect with Partial Frequency Redistribution

In the previous two papers, namely, \citet{anuknn11} and \citet{anuetal11} we solved the polarized radiative transfer (RT) equation in multi-dimensional (multi-D) geometries, with partial frequency redistribution (PRD) as the scattering mechanism. We assumed Rayleigh scattering as the only source of linear polarization ($Q/I, U/I$) in both these papers. In this paper we extend these previous works to include the effect of weak oriented magnetic fields (Hanle effect) on line scattering. We generalize the technique of Stokes vector decomposition in terms of the irreducible spherical tensors $\mathcal{T}^K_Q$, developed in \citet{anuknn11}, to the case of RT with Hanle effect. A fast iterative method of solution (based on the Stabilized Preconditioned Bi-Conjugate-Gradient technique), developed in \citet{anuetal11}, is now generalized to the case of RT in magnetized three-dimensional media. We use the efficient short-characteristics formal solution method for multi-D media, generalized appropriately to the present context. The main results of this paper are the following: (1) A comparison of emergent $(I, Q/I, U/I)$ profiles formed in one-dimensional (1D) media, with the corresponding emergent, spatially averaged profiles formed in multi-D media, shows that in the spatially resolved structures, the assumption of 1D may lead to large errors in linear polarization, especially in the line wings. (2) The multi-D RT in semi-infinite non-magnetic media causes a strong spatial variation of the emergent $(Q/I, U/I)$ profiles, which is more pronounced in the line wings. (3) The presence of a weak magnetic field modifies the spatial variation of the emergent $(Q/I, U/I)$ profiles in the line core, by producing significant changes in their magnitudes.

preprint2009arXiv

Origin of spatial variations of scattering polarization in the wings of the Ca {\sc i} 4227 Åline

Polarization that is produced by coherent scattering can be modified by magnetic fields via the Hanle effect. According to standard theory the Hanle effect should only be operating in the Doppler core of spectral lines but not in the wings. In contrast, our observations of the scattering polarization in the Ca {\sc i} 4227 Åline reveals the existence of spatial variations of the scattering polarization throughout the far line wings. This raises the question whether the observed spatial variations in wing polarization have a magnetic or non-magnetic origin. A magnetic origin may be possible if elastic collisions are able to cause sufficient frequency redistribution to make the Hanle effect effective in the wings without causing excessive collisional depolarization, as suggested by recent theories for partial frequency redistribution with coherent scattering in magnetic fields. To model the wing polarization we apply an extended version of the technique based on the "last scattering approximation". This model is highly successful in reproducing the observed Stokes $Q/I$ polarization (linear polarization parallel to the nearest solar limb), including the location of the wing polarization maxima and the minima around the Doppler core, but it fails to reproduce the observed spatial variations of the wing polarization in terms of magnetic field effects with frequency redistribution. This null result points in the direction of a non-magnetic origin in terms of local inhomogeneities (varying collisional depolarization, radiation-field anisotropies, and deviations from a plane-parallel atmospheric stratification).

preprint2009arXiv

Preconditioned Bi-Conjugate Gradient Method for Radiative Transfer in Spherical Media

A robust numerical method called the Preconditioned Bi-Conjugate Gradient (Pre-BiCG)method is proposed for the solution of radiative transfer equation in spherical geometry.A variant of this method called Stabilized Preconditioned Bi-Conjugate Gradient (Pre-BiCG-STAB) is also presented. These are iterative methods based on the construction of a set of bi-orthogonal vectors. The application of Pre-BiCG method in some benchmark tests show that the method is quite versatile, and can handle hard problems that may arise in astrophysical radiative transfer theory.