Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userPartial redistribution in 3D non-LTE radiative transfer in solar atmosphere models
176
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Resonance spectral lines such as H I Ly {alpha}, Mg II h&k, and Ca II H&K that form in the solar chromosphere are influenced by the effects of 3D radiative transfer as well as partial redistribution (PRD). So far no one has modeled these lines including both effects simultaneously owing to the high computing demands of existing algorithms. Such modeling is however indispensable for accurate diagnostics of the chromosphere. We present a computationally tractable method to treat PRD scattering in 3D model atmospheres using a 3D non-LTE radiative transfer code. To make the method memory-friendly, we use the hybrid approximation of Leenaarts et al. (2012) for the redistribution integral. To make it fast, we use linear interpolation on equidistant frequency grids. We verify our algorithm against computations with the RH code and analyze it for stability, convergence, and usefulness of acceleration using model atoms of Mg II with the h&k lines and H I with the Ly {alpha} line treated in PRD. A typical 3D PRD solution can be obtained in a model atmosphere with $252 times 252 times 496$ coordinate points in 50 000--200 000 CPU hours, which is a factor ten slower than computations assuming complete redistribution. We illustrate the importance of the joint action of PRD and 3D effects for the Mg II h&k lines for disk-center intensities as well as the center-to-limb variation. The proposed method allows simulating PRD lines in time series of radiation-MHD models in order to interpret observations of chromospheric lines at high spatial resolution.
rate research
Read More
The reconstruction of the solar spectral irradiance (SSI) on various time scales is essential for the understanding of the Earths climate response to the SSI variability. The driver of the SSI variability is understood to be the intensity contrast of magnetic features present on the Sun with respect to the largely non-magnetic quiet Sun. However, different spectral synthesis codes lead to diverging projections of SSI variability. In this study we compare three different radiative transfer codes and carry out a detailed analysis of their performance. We perform the spectral synthesis at the continuum wavelength of 665 nm with the Code for Solar Irradiance (COSI), and the Rybicki-Hummer (RH), and Max Planck University of Chicago Radiative MHD (MURaM) codes for three 3D MHD simulations snapshots, a non-magnetic case, and MHD simulations with 100 G, and 200 G magnetic field strength. We determine the intensity distributions, the intensity differences and ratios for the spectral synthesis codes. We identify that the largest discrepancies originate in the intergranular lanes where the most field concentration occurs. Overall, the applied radiative transfer codes give consistent intensity distributions. Also, the intensity variation as a function of magnetic field strength for the particular 100 G and 200 G snapshots agree within the 2-3% range.
Nitrogen is an important element in various fields of stellar and Galactic astronomy, and the solar nitrogen abundance is crucial as a yardstick for comparing different objects in the cosmos. In order to obtain a precise and accurate value for this abundance, we carried out N i line formation calculations in a 3D radiative-hydrodynamic STAGGER model solar atmosphere, in full 3D non-local thermodynamic equilibrium (non-LTE), using a model atom that includes physically-motivated descriptions for the inelastic collisions of N i with free electrons and with neutral hydrogen. We selected five N i lines of high excitation energy to study in detail, based on their strengths and on their being relatively free of blends. We found that these lines are slightly strengthened from non-LTE photon losses and from 3D granulation effects, resulting in negative abundance corrections of around $-0.01$ dex and $-0.04$ dex respectively. Our advocated solar nitrogen abundance is $logepsilon_{mathrm{N}} = 7.77$, with the systematic $1sigma$ uncertainty estimated to be $0.05$ dex. This result is consistent with earlier studies after correcting for differences in line selections and equivalent widths.
Context. The pursuit of more realistic spectroscopic modelling and consistent abundances has led us to begin a new series of papers designed to improve current solar and stellar abundances of various atomic species. To achieve this, we have began updating the three-dimensional (3D) non-local thermodynamic equilibrium (non-LTE) radiative transfer code, Multi3D, and the equivalent one-dimensional (1D) non-LTE radiative transfer code, MULTI. Aims. We examine our improvements to these codes by redetermining the solar barium abundance. Barium was chosen for this test as it is an important diagnostic element of the s-process in the context of galactic chemical evolution. New Ba II + H collisional data for excitation and charge exchange reactions computed from first principles had recently become available and were included in the model atom. The atom also includes the effects of isotopic line shifts and hyperfine splitting. Method. A grid of 1D LTE barium lines were constructed with MULTI and fit to the four Ba II lines available to us in the optical region of the solar spectrum. Abundance corrections were then determined in 1D non-LTE, 3D LTE, and 3D non-LTE. A new 3D non-LTE solar barium abundance was computed from these corrections. Results. We present for the first time the full 3D non-LTE barium abundance of $A({rm Ba})=2.27pm0.02pm0.01$, which was derived from four individual fully consistent barium lines. Errors here represent the systematic and random errors, respectively.
Context. Small-scale bright features in the photosphere of the Sun, such as faculae or G-band bright points, appear in connection with small-scale magnetic flux concentrations. Aims. Here we report on a new class of photospheric bright points that are free of magnetic fields. So far, these are visible in numerical simulations only. We explore conditions required for their observational detection. Methods. Numerical radiation (magneto-)hydrodynamic simulations of the near-surface layers of the Sun were carried out. The magnetic field-free simulations show tiny bright points, reminiscent of magnetic bright points, only smaller. A simple toy model for these non-magnetic bright points (nMBPs) was established that serves as a base for the development of an algorithm for their automatic detection. Basic physical properties of 357 detected nMBPs were extracted and statistically evaluated. We produced synthetic intensity maps that mimic observations with various solar telescopes to obtain hints on their detectability. Results. The nMBPs of the simulations show a mean bolometric intensity contrast with respect to their intergranular surroundings of approximately 20%, a size of 60-80 km, and the isosurface of optical depth unity is at their location depressed by 80-100 km. They are caused by swirling downdrafts that provide, by means of the centripetal force, the necessary pressure gradient for the formation of a funnel of reduced mass density that reaches from the subsurface layers into the photosphere. Similar, frequently occurring funnels that do not reach into the photosphere, do not produce bright points. Conclusions. Non-magnetic bright points are the observable manifestation of vertically extending vortices (vortex tubes) in the photosphere. The resolving power of 4-m-class telescopes, such as the DKIST, is needed for an unambiguous detection of them.
The Interface Region Imaging Spectrograph (IRIS) routinely observes the Si IV resonance lines. When analyzing observations of these lines it has typically been assumed they form under optically thin conditions. This is likely valid for the quiescent Sun, but this assumption has also been applied to the more extreme flaring scenario. We used 36 electron beam driven radiation hydrodynamic solar flare simulations, computed using the RADYN code, to probe the validity of this assumption. Using these simulated atmospheres we solved the radiation transfer equations to obtain the non-LTE, non-equilibrium populations, line profiles, and opacities for a model Silicon atom, including charge exchange processes. This was achieved using the `minority species version of RADYN. The inclusion of charge exchange resulted in a substantial fraction of Si IV at cooler temperatures than those predicted by ionisation equilibrium. All simulations with an injected energy flux $F>5times10^{10}$ erg cm$^{-2}$ s$^{-1}$ resulted in optical depth effects on the Si IV emission, with differences in both intensity and line shape compared to the optically thin calculation. Weaker flares (down to $Fapprox5times10^{9}$ erg cm$^{-2}$ s$^{-1}$) also resulted in Si IV emission forming under optically thick conditions, depending on the other beam parameters. When opacity was significant, the atmospheres generally had column masses in excess of $5times10^{-6}$ g cm$^{-2}$ over the temperature range $40$ to $100$ kK, and the Si IV formation temperatures were between $30$ and $60$ kK. We urge caution when analyzing Si IV flare observations, or when computing synthetic emission without performing a full radiation transfer calculation.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
