Partial Frequency Redistribution Research Articles

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

Abstract 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 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 i D2 line. For this case we show that the AD-PRD effects are significant for field strengths below about 30 G, but that the computationally much less demanding approximation of angle-averaged PRD may be used for stronger fields.

Polarized Line Formation in Arbitrary Strength Magnetic Fields: The Case of a Two-level Atom with Hyperfine Structure Splitting

Abstract Quantum interference effects, together with partial frequency redistribution (PFR) in line scattering, produce subtle signatures in the so-called Second Solar Spectrum (the linearly polarized spectrum of the Sun). These signatures are modified in the presence of arbitrary strength magnetic fields via the Hanle, Zeeman, and Paschen–Back effects. In the present paper we solve the problem of polarized line formation in a magnetized atmosphere taking into account scattering in a two-level atom with hyperfine structure splitting together with PFR. To this end we incorporate the collisionless PFR matrix derived in Sowmya et al. in the polarized transfer equation. We apply the scattering expansion method to solve this transfer equation. We study the combined effects of PFR and the Paschen–Back effect on polarized line profiles formed in an isothermal one-dimensional planar atmosphere. For this purpose, we consider the cases of D2 lines of Li i and Na i.

Modeling Mg II During Solar Flares. I. Partial Frequency Redistribution, Opacity, and Coronal Irradiation

Abstract The Interface Region Imaging Spectrograph has routinely observed the flaring Mg ii near-ultraviolet (NUV) spectrum, offering excellent diagnostic potential and a window into the location of energy deposition. A number of studies have forward-modeled both the general properties of these lines and specific flare observations. Generally these have forward-modeled radiation via post-processing of snapshots from hydrodynamic flare simulations through radiation transfer codes. There has, however, not been a study of how the physics included in these radiation transport codes affects the solution. A baseline setup for forward-modeling Mg ii in flares is presented and contrasted with approaches that add or remove complexity. It is shown for Mg ii that (1) partial frequency distribution (PRD) is still required during flare simulations despite the increased densities; (2) using full angle-dependent PRD affects the solution but takes significantly longer to process a snapshot; (3) including Mg i in non-LTE (NLTE) results in negligible differences to the Mg ii lines but does affect the NUV quasi-continuum; (4) only hydrogen and Mg ii need to be included in NLTE; (5) ideally the nonequilibrium hydrogen populations, with nonthermal collisional rates, should be used rather than the statistical equilibrium populations; (6) an atom consisting of only the ground state, h and k upper levels, and continuum level is insufficient to model the resonance lines; and (7) irradiation from a hot, dense flaring transition region can affect the formation of Mg ii. We discuss modifications to the RH code allowing straightforward inclusion of the transition region and coronal irradiation in flares.

Magnetic Sensitivity in the Wing Scattering Polarization Signals of the Hydrogen Lyman-α Line of the Solar Disk Radiation

Abstract The linear polarization produced by scattering processes in the hydrogen Lyα line of the solar disk radiation is a key observable for probing the chromosphere–corona transition region (TR) and the underlying chromospheric plasma. While the line-center signal encodes information on the magnetic field and the three-dimensional structure of the TR, the sizable scattering polarization signals that the joint action of partial frequency redistribution and J-state interference produce in the Lyα wings have generally been thought to be sensitive only to the thermal structure of the solar atmosphere. Here we show that the wings of the Q/I and U/I scattering polarization profiles of this line are actually sensitive to the presence of chromospheric magnetic fields, with strengths similar to those that produce the Hanle effect in the line core (i.e., between 5 and 100 G, approximately). In spite of the fact that the Zeeman splitting induced by such weak fields is very small compared to the total width of the line, the magneto-optical effects that couple the transfer equations for Stokes Q and U are actually able to produce sizable changes in the Q/I and U/I wings. We find that magnetic fields with longitudinal components larger than 100 G produce an almost complete depolarization of the wings of the Lyα Q/I profiles within a ±5 Å spectral range around the line center, while stronger fields are required for the U/I wing signals to be depolarized to a similar extent. The theoretical results presented here further expand the diagnostic content of the unprecedented spectropolarimetric observations provided by the Chromospheric Lyman-Alpha Spectro-Polarimeter.

The Response of the Lyα Line in Different Flare Heating Models

Abstract The solar Lyα line is the strongest line in the ultraviolet waveband, and is greatly enhanced during solar flares. Here we present radiative hydrodynamic simulations of solar flares under different heating models, and calculate the response of this line taking into account nonequilibrium ionization of hydrogen and partial frequency redistribution. We find that in nonthermal heating models, the Lyα line can show a red or blue asymmetry corresponding to the chromospheric evaporation or condensation, respectively. The asymmetry may change from red to blue if the electron beam flux is large enough to produce a significant chromospheric condensation region. In the Lyα intensity light curve, a dip appears when the change of asymmetry occurs. In thermal models, the Lyα line intensity peaks quickly and then falls, and the profile has an overall red asymmetry, which is similar to the profiles from heating by a soft electron beam. The Lyα profile shows a single red peak at the end of thermal heating, and the whole line is formed in a very small height range.

Predicting the Extreme Ultraviolet Radiation Environment of Exoplanets around Low-mass Stars: The TRAPPIST-1 System

Abstract The high energy radiation environment around M dwarf stars strongly impacts the characteristics of close-in exoplanet atmospheres, but these wavelengths are difficult to observe due to geocoronal and interstellar contamination. On account of these observational restrictions, a stellar atmosphere model may be used to compute the stellar extreme ultraviolet (EUV; 100–912 Å) spectrum. We construct semiempirical nonlocal thermodynamic equilibrium model spectra of the ultracool M8 star TRAPPIST-1 that span EUV to infrared wavelengths (100 Å−2.5 μm) using the atmosphere code PHOENIX. These upper atmosphere models contain prescriptions for the chromosphere and transition region and include newly added partial frequency redistribution capabilities. In the absence of broadband UV spectral observations, we constrain our models using Hubble Space Telescope Lyman α observations from TRAPPIST-1 and Galaxy Evolution Explorer UV photometric detections from a set of old M8 stars (>1 Gyr). We find that calibrating the models using both data sets separately yield similar far-ultraviolet and NUV fluxes, and EUV fluxes that range from (1.32–17.4) × 10−14 ergs s−1 cm−2. The results from these models demonstrate that the EUV emission is very sensitive to the temperature structure in the transition region. Our lower activity models predict EUV fluxes similar to previously published estimates derived from semiempirical scaling relationships, while the highest activity model predicts EUV fluxes a factor of 10 higher. Results from this study support the idea that the TRAPPIST-1 habitable zone planets likely do not have much liquid water on their surfaces due to the elevated levels of high energy radiation emitted by the host star.

Temperature constraints from inversions of synthetic solar optical, UV, and radio spectra

Context. High-resolution observations of the solar chromosphere at millimeter wavelengths are now possible with the Atacama Large Millimeter Array (ALMA), bringing with them the promise of tackling many open problems in solar physics. Observations from other ground and space-based telescopes will greatly benefit from coordinated endeavors with ALMA, yet the diagnostic potential of combined optical, ultraviolet and mm observations has remained mostly unassessed. Aims. In this paper we investigate whether mm-wavelengths could aid current inversion schemes to retrieve a more accurate representation of the temperature structure of the solar atmosphere. Methods. We performed several non-LTE inversion experiments of the emergent spectra from a snapshot of 3D radiation-MHD simulation. We included common line diagnostics such as Ca II H, K, 8542 Å and Mg II h and k, taking into account partial frequency redistribution effects, along with the continuum around 1.2 mm and 3 mm. Results. We find that including the mm-continuum in inversions allows a more accurate inference of temperature as function of optical depth. The addition of ALMA bands to other diagnostics should improve the accuracy of the inferred chromospheric temperatures between log τ ∼ [−6, −4.5] where the Ca II and Mg II lines are weakly coupled to the local conditions. However, we find that simultaneous multiatom, non-LTE inversions of optical and UV lines present equally strong constraints in the lower chromosphere and thus are not greatly improved by the 1.2 mm band. Nonetheless, the 3 mm band is still needed to better constrain the mid-upper chromosphere.

Doppler Frequency Redistribution upon Coherent Photon Emission by Atoms in an Optically Dense Medium

The problem of transfer of line radiation in sodium vapor under photoexcitation by laser radiation of a resonance line with wavelength λ = 589.6 nm is numerically solved. The influence of Doppler frequency shifts on photon emission coherent in the atomic reference system is taken into account in the formation of an emission line profile using the model of partial frequency redistribution. In the case of a small optical thickness of a medium τ0 ≤ 1, the Doppler frequency redistribution in the laboratory reference system yields higher intensity values in the emission line profile core in comparison with the model of complete frequency redistribution. For an optically dense medium, thermal motion of atoms leads to an increase of reabsorption in the emission line core, and intensity increases in the line wings as compared to the complete redistribution model. The dependence of the Holstein trapping factor on optical thickness is closer to the analytic dependence than the complete redistribution model. This is explained by the fact that the frequencies of emitted photons in the laboratory reference system fall into the spectral profile core, which increases the effects of resonance radiation trapping in a dense medium.

Magneto-optical Effects in the Scattering Polarization Wings of the Ca i 4227 Å Resonance Line

Abstract The linear polarization pattern produced by scattering processes in the Ca i 4227 Å resonance line is a valuable observable for probing the solar atmosphere. Via the Hanle effect, the very significant Q/I and U/I line-center signals are sensitive to the presence of magnetic fields in the lower chromosphere with strengths between 5 and 125 G, approximately. On the other hand, partial frequency redistribution (PRD) produces sizable signals in the wings of the Q/I profile, which have always been thought to be insensitive to the presence of magnetic fields. Interestingly, novel observations of this line revealed a surprising behavior: fully unexpected signals in the wings of the U/I profile and spatial variability in the wings of both Q/I and U/I. We show that the magneto-optical (MO) terms of the Stokes-vector transfer equation produce sizable signals in the wings of U/I and a clear sensitivity of the Q/I and U/I wings to the presence of photospheric magnetic fields with strengths similar to those that produce the Hanle effect in the line core. This radiative transfer investigation on the joint action of scattering processes and the Hanle and Zeeman effects in the Ca i 4227 Å line should facilitate the development of more reliable techniques for exploring the magnetism of stellar atmospheres. To this end, we can now exploit the circular polarization produced by the Zeeman effect, the magnetic sensitivity caused by the above-mentioned MO effects in the Q/I and U/I wings, and the Hanle effect in the line core.

Rayleigh Scattering in Spectral Series with L-term Interference

Abstract We derive a formalism to describe the scattering of polarized radiation over the full spectral range encompassed by atomic transitions belonging to the same spectral series (e.g., the H i Lyman and Balmer series, the UV multiplets of Fe i and Fe ii). This allows us to study the role of radiation-induced coherence among the upper terms of the spectral series, and its contribution to Rayleigh scattering and the polarization of the solar continuum. We rely on previous theoretical results for the emissivity of a three-term atom of the Λ-type, taking into account partially coherent scattering, and generalize its expression in order to describe a “multiple Λ” atomic system underlying the formation of a spectral series. Our study shows that important polarization effects must be expected because of the combined action of partial frequency redistribution and radiation-induced coherence among the terms of the series. In particular, our model predicts the correct asymptotic limit of 100% polarization in the far wings of a complete (i.e., ) group of transitions, which must be expected on the basis of the principle of spectroscopic stability.

Polarized Line Formation in Arbitrary Strength Magnetic Fields Angle-averaged and Angle-dependent Partial Frequency Redistribution

Abstract Magnetic fields in the solar atmosphere leave their fingerprints in the polarized spectrum of the Sun via the Hanle and Zeeman effects. While the Hanle and Zeeman effects dominate, respectively, in the weak and strong field regimes, both these effects jointly operate in the intermediate field strength regime. Therefore, it is necessary to solve the polarized line transfer equation, including the combined influence of Hanle and Zeeman effects. Furthermore, it is required to take into account the effects of partial frequency redistribution (PRD) in scattering when dealing with strong chromospheric lines with broad damping wings. In this paper, we present a numerical method to solve the problem of polarized PRD line formation in magnetic fields of arbitrary strength and orientation. This numerical method is based on the concept of operator perturbation. For our studies, we consider a two-level atom model without hyperfine structure and lower-level polarization. We compare the PRD idealization of angle-averaged Hanle–Zeeman redistribution matrices with the full treatment of angle-dependent PRD, to indicate when the idealized treatment is inadequate and what kind of polarization effects are specific to angle-dependent PRD. Because the angle-dependent treatment is presently computationally prohibitive when applied to realistic model atmospheres, we present the computed emergent Stokes profiles for a range of magnetic fields, with the assumption of an isothermal one-dimensional medium.

Importance of Cross-redistribution in Scattering Polarization of Spectral Lines:The Cases of 3P−3S Triplets of Mg i and Ca i

Abstract Scattering on a multi-level atomic system has dominant contributions from resonance and Raman scattering. While initial and final levels are the same for resonance scattering, they are different for Raman scattering. The frequency redistribution for resonance scattering is described by the usual partial frequency redistribution functions of Hummer, while that for Raman scattering is described by cross-redistribution (XRD) function. In the present paper, we investigate the importance of XRD on linear polarization profiles of 3P−3S triplets of Mg i and Ca i formed in an isothermal one-dimensional atmosphere. We show that XRD produces significant effects on the linear polarization profiles when the wavelength separations between the line components of the multiplet are small, like in the cases of Mg i b and Ca i triplets.

The Transfer of Resonance Line Polarization with Partial Frequency Redistribution in the General Hanle–Zeeman Regime

Abstract The spectral line polarization encodes a wealth of information about the thermal and magnetic properties of the solar atmosphere. Modeling the Stokes profiles of strong resonance lines is, however, a complex problem both from a theoretical and computational point of view, especially when partial frequency redistribution (PRD) effects need to be taken into account. In this work, we consider a two-level atom in the presence of magnetic fields of arbitrary intensity (Hanle–Zeeman regime) and orientation, both deterministic and micro-structured. Working within the framework of a rigorous PRD theoretical approach, we have developed a numerical code that solves the full non-LTE radiative transfer problem for polarized radiation, in one-dimensional models of the solar atmosphere, accounting for the combined action of the Hanle and Zeeman effects, as well as for PRD phenomena. After briefly discussing the relevant equations, we describe the iterative method of solution of the problem and the numerical tools that we have developed and implemented. We finally present some illustrative applications to two resonance lines that form at different heights in the solar atmosphere, and provide a detailed physical interpretation of the calculated Stokes profiles. We find that magneto-optical effects have a strong impact on the linear polarization signals that PRD effects produce in the wings of strong resonance lines. We also show that the weak-field approximation has to be used with caution when PRD effects are considered.

POLARIZED LINE FORMATION IN NON-MONOTONIC VELOCITY FIELDS

ABSTRACT For a correct interpretation of the observed spectro-polarimetric data from astrophysical objects such as the Sun, it is necessary to solve the polarized line transfer problems taking into account a realistic temperature structure, the dynamical state of the atmosphere, a realistic scattering mechanism (namely, the partial frequency redistribution—PRD), and the magnetic fields. In a recent paper, we studied the effects of monotonic vertical velocity fields on linearly polarized line profiles formed in isothermal atmospheres with and without magnetic fields. However, in general the velocity fields that prevail in dynamical atmospheres of astrophysical objects are non-monotonic. Stellar atmospheres with shocks, multi-component supernova atmospheres, and various kinds of wave motions in solar and stellar atmospheres are examples of non-monotonic velocity fields. Here we present studies on the effect of non-relativistic non-monotonic vertical velocity fields on the linearly polarized line profiles formed in semi-empirical atmospheres. We consider a two-level atom model and PRD scattering mechanism. We solve the polarized transfer equation in the comoving frame (CMF) of the fluid using a polarized accelerated lambda iteration method that has been appropriately modified for the problem at hand. We present numerical tests to validate the CMF method and also discuss the accuracy and numerical instabilities associated with it.

THE MAGNETIC SENSITIVITY OF THE Mg ii k LINE TO THE JOINT ACTION OF HANLE, ZEEMAN, AND MAGNETO-OPTICAL EFFECTS

ABSTRACT We highlight the main results of a radiative transfer investigation on the magnetic sensitivity of the solar Mg ii k resonance line at 2795.5 Å, accounting for the joint action of the Hanle and Zeeman effects as well as partial frequency redistribution phenomena. We confirm that at the line center, the linear polarization signals produced by scattering processes are measurable, and that they are sensitive, via the Hanle effect, to magnetic fields with strengths between 5 and 50 G, approximately. We also show that the Zeeman effect produces conspicuous circular polarization signals, especially for longitudinal fields stronger than 50 G, which can be used to estimate the magnetization of the solar chromosphere via the familiar magnetograph formula. The most novel result is that magneto-optical effects produce, in the wings of the line, a decrease of the Q/I scattering polarization pattern and the appearance of U/I signals (i.e., a rotation of the plane of linear polarization). This sensitivity of the Q/I and U/I wing signals to both weak (∼5 G) and stronger magnetic fields expands the scientific interest of the Mg ii k line for probing the chromosphere in quiet and active regions of the Sun.

POLARIZED LINE FORMATION WITH LOWER-LEVEL POLARIZATION AND PARTIAL FREQUENCY REDISTRIBUTION

ABSTRACT In the well-established theories of polarized line formation with partial frequency redistribution (PRD) for a two-level and two-term atom, it is generally assumed that the lower level of the scattering transition is unpolarized. However, the existence of unexplained spectral features in some lines of the Second Solar Spectrum points toward a need to relax this assumption. There exists a density matrix theory that accounts for the polarization of all the atomic levels, but it is based on the flat-spectrum approximation (corresponding to complete frequency redistribution). In the present paper we propose a numerical algorithm to solve the problem of polarized line formation in magnetized media, which includes both the effects of PRD and the lower level polarization (LLP) for a two-level atom. First we derive a collisionless redistribution matrix that includes the combined effects of the PRD and the LLP. We then solve the relevant transfer equation using a two-stage approach. For illustration purposes, we consider two case studies in the non-magnetic regime, namely, the J a = 1, J b = 0 and J a = J b = 1, where J a and J b represent the total angular momentum quantum numbers of the lower and upper states, respectively. Our studies show that the effects of LLP are significant only in the line core. This leads us to propose a simplified numerical approach to solve the concerned radiative transfer problem.

Unified numerical model of collisional depolarization and broadening rates that are due to hydrogen atom collisions

Interpretation of solar polarization spectra accounting for partial or complete frequency redistribution requires data on various collisional processes. Data for depolarization and polarization transfer are needed but often missing, while data for collisional broadening are usually more readily available. Recent work by Sahal-Br\'echot and Bommier concluded that despite underlying similarities in the physics of collisional broadening and depolarization processes, relationships between them are not possible to derive purely analytically. We aim to derive accurate numerical relationships between the collisional broadening rates and the collisional depolarization and polarization transfer rates due to hydrogen atom collisions. Such relationships would enable accurate and efficient estimation of collisional data for solar applications. Using earlier results for broadening and depolarization processes based on general (i.e. not specific to a given atom), semi-classical calculations employing interaction potentials from perturbation theory, genetic programming (GP) has been used to fit the available data and generate analytical functions describing the relationships between them. The predicted relationships from the GP-based model are compared with the original data to estimate the accuracy of the method.

POLARIZED LINE FORMATION IN MOVING ATMOSPHERES WITH PARTIAL FREQUENCY REDISTRIBUTION AND A WEAK MAGNETIC FIELD

The dynamical state of the solar and stellar atmospheres depends on the macroscopic velocity fields prevailing within them. The presence of such velocity fields in the line formation regions strongly affects the polarized radiation field emerging from these atmospheres. Thus it becomes necessary to solve the radiative transfer equation for polarized lines in moving atmospheres. Solutions based on the "observer's frame method" are computationally expensive to obtain, especially when partial frequency redistribution (PRD) in line scattering and large-amplitude velocity fields are taken into account. In this paper we present an efficient alternative method of solution, namely, the comoving frame technique, to solve the polarized PRD line formation problems in the presence of velocity fields. We consider one-dimensional planar isothermal atmospheres with vertical velocity fields. We present a study of the effect of velocity fields on the emergent linear polarization profiles formed in optically thick moving atmospheres. We show that the comoving frame method is far superior when compared to the observer's frame method in terms of the computational speed and memory requirements.

Retaining space and time coherence in radiative transfer models.

A recent model for radiative transfer that accounts for spatial coherence is extended in such a way as to retain temporal coherence. The method employs Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy techniques. Both spatial and temporal coherence are shown to affect the formation of atomic line spectra. Calculations of Lyman α radiation transport in optically thick divertor plasma conditions are reported as an illustration of the model. A possible extension of the formalism to dense media involving correlations between atoms is discussed in an appendix. A link to partial frequency redistribution modeling is also discussed.

RH 1.5D: a massively parallel code for multi-level radiative transfer with partial frequency redistribution and Zeeman polarisation

The emergence of three-dimensional magneto-hydrodynamic (MHD) simulations of stellar atmospheres has sparked a need for efficient radiative transfer codes to calculate detailed synthetic spectra. We present RH 1.5D, a massively parallel code based on the RH code and capable of performing Zeeman polarised multi-level non-local thermodynamical equilibrium (NLTE) calculations with partial frequency redistribution for an arbitrary amount of chemical species. The code calculates spectra from 3D, 2D or 1D atmospheric models on a column-by-column basis (or 1.5D). While the 1.5D approximation breaks down in the cores of very strong lines in an inhomogeneous environment, it is nevertheless suitable for a large range of scenarios and allows for faster convergence with finer control over the iteration of each simulation column. The code scales well to at least tens of thousands of CPU cores, and is publicly available. In the present work we briefly describe its inner workings, strategies for convergence optimisation, its parallelism, and some possible applications.