Source author record

Ivan Milic

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

ResearcherUnclaimed source record

astro-ph.SR

Catalog footprint

What is connected

5works

1topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Evaluating Non-LTE Spectral Inversions with ALMA and IBIS

We present observations of a solar plage in the millimeter-continuum with the ALMA and in the Ca 8542 and Na 5896 spectral lines with the Interferometric BIdimensional Spectrometer (IBIS). Our goal is to compare the measurement of local gas temperatures provided by ALMA with the temperature diagnostics provided by non-LTE inversions using STIC. In performing these inversions, we find that using column mass as the reference height scale, rather than optical depth, provides more reliable atmospheric profiles above the temperature minimum and that the treatment of non- LTE hydrogen ionization brings the inferred chromospheric temperatures into better agreement with the ALMA measurements. The Band 3 brightness temperatures are higher but well correlated with the inversion-derived temperatures at the height of formation of the Ca 8542 line core. The Band 6 temperatures instead do not show good correlations with the temperatures at any specific layer in the inverted atmospheres. We then performed inversions that included the millimeter continuum intensities as an additional constraint. Incorporating Band 3 generally resulted in atmospheres showing a strong temperature rise in the upper atmosphere, while including Band 6 led to significant regions of anomalously low temperatures at chromospheric heights. This is consistent with the idea that the Band 6 emission can come from a range of heights. The poor constraints on the chromospheric electron density with existing inversion codes introduces difficulties in determining the height(s) of formation of the millimeter continuum as well as uncertainties in the temperatures derived from the spectral lines.

preprint2016arXiv

Inference of magnetic fields in inhomogeneous prominences

Most of the quantitative information about the magnetic field vector in solar prominences comes from the analysis of the Hanle effect acting on lines formed by scattering. As these lines can be of non-negligible optical thickness, it is of interest to study the line formation process further. We investigate the multidimensional effects on the interpretation of spectropolarimetric observations, particularly on the inference of the magnetic field vector. We do this by analyzing the differences between multidimensional models, which involve fully self-consistent radiative transfer computations in the presence of spatial inhomogeneities and velocity fields, and those which rely on simple one-dimensional geometry. We study the formation of a prototype line in ad hoc inhomogeneous, isothermal 2D prominence models. We solve the NLTE polarized line formation problem in the presence of a large-scale oriented magnetic field. The resulting polarized line profiles are then interpreted (i.e. inverted) assuming a simple 1D slab model. We find that differences between input and the inferred magnetic field vector are non-negligible. Namely, we almost universally find that the inferred field is weaker and more horizontal than the input field. Spatial inhomogeneities and radiative transfer have a strong effect on scattering line polarization in the optically thick lines. In real-life situations, ignoring these effects could lead to a serious misinterpretation of spectropolarimetric observations of chromospheric objects such as prominences.

preprint2014arXiv

Accelerating NLTE radiative transfer by means of the Forth-and-Back Implicit Lambda Iteration: A two-level atom line formation in 2D Cartesian coordinates

State-of-the-art methods in multidimensional NLTE radiative transfer are based on the use of local approximate lambda operator within either Jacobi or Gauss-Seidel iterative schemes. Here we propose another approach to the solution of 2D NLTE RT problems, Forth-and-Back Implicit Lambda Iteration (FBILI), developed earlier for 1D geometry. In order to present the method and examine its convergence properties we use the well-known instance of the two-level atom line formation with complete frequency redistribution. In the formal solution of the RT equation we employ short characteristics with two-point algorithm. Using an implicit representation of the source function in the computation of specific intensities, we compute and store the coefficients of the linear relations J = a + bS between the mean intensity J and the corresponding source function S. The use of iteration factors in the 'local' coefficients of these implicit relations in two 'inward' directions, along with the update of the source function in other two, 'outward', directions leads to four times faster solution than the Jacobi's one. Moreover, the update made in all four consecutive sweeps of the grid leads to an acceleration by a factor of 6-7 compared to the Jacobi iterative scheme.

preprint2014arXiv

Scattering line polarization in rotating, optically thick disks

To interpret observations of astrophysical disks it is essential to understand the formation process of the emitted light. If the disk is optically thick, scattering dominated and permeated by a Keplerian velocity field, Non-Local Thermodynamic Equilibrium radiative transfer modeling must be done to compute the emergent spectrum from a given disk model. We investigate Non-local thermodynamic equilibrium polarized line formation in different simple disk models and aim to demonstrate the importance of both radiative transfer effects and scattering as well as the effects of velocity fields. We self-consistently solve the coupled equations of radiative transfer and statistical equilibrium for a two level atom model by means of Jacobi iteration. We compute scattering polarization, that is Q/I and U/I line profiles. The degree of scattering polarization is significantly influenced by the inclination of the disk with respect to observer, but also by the optical thickness of the disk and the presence of rotation. Stokes U shows double-lobed profiles with amplitude which increases with the disk rotation. Our results suggest that the line profiles, especially the polarized ones, emerging from gaseous disks differ significantly from the profiles predicted by simple approximations. The profiles are diverse in shape, but typically symmetric in Stokes Q and antisymmetric in Stokes U. A clear indicator of disk rotation is the presence of Stokes U, which might prove to be a useful diagnostic tool. We also demonstrate that, for moderate rotational velocities, an approximate treatment can be used, where non-local thermodynamic equilibrium radiative transfer is done in the velocity field-free approximation and Doppler shift is applied in the process of spatial integration over the whole emitting surface.

preprint2013arXiv

Boundary conditions for NLTE polarized radiative transfer with incident radiation

Polarized NLTE radiative transfer in the presence of scattering in spectral lines and/or in continua may be cast in a so-called reduced form for six reduced components of the radiation field. In this formalism the six components of the reduced source function are angle-independent quantities. It thus reduces drastically the storage requirement of numerical codes. This approach encounters a fundamental problem when the medium is illuminated by a polarized incident radiation, because there is a priori no way of relating the known (and measurable) Stokes parameters of the incident radiation to boundary conditions for the reduced equations. The origin of this problem is that there is no unique way of deriving the radiation reduced components from its Stokes parameters (only the inverse operation is clearly defined). The method proposed here aims at enabling to work with arbitrary incident radiation field (polarized or unpolarized). In previous works an ad-hoc treatment of the boundary conditions, applying to cases where the incident radiation is unpolarized, has been used. In this note we show that it is possible to account for the incident radiation in a rigorous way, without any assumption on its properties, by expressing the radiation field as the sum of a directly transmitted radiation and of a so-called diffuse radiation. The diffuse radiation field obeys a transfer equation with no incident radiation that may be solved in the reduced form. The first scattering of the incident radiation introduces primary creation terms in the six components of the reduced source function. Once the reduced polarized transfer problem is solved for the diffuse radiation field, its Stokes parameters can be computed. We perform numerical computations of such cases, showing that the emergent line-polarization may be significantly affected by the polarization of the incident radiation.