NLTE modeling of Stokes vector center-to-limb variations in the CN violet system

The solar surface magnetic field is connected with and even controls most of the solar activity phenomena. Zeeman effect diagnostics allow for measuring only a small fraction of the fractal-like structured magnetic field. The remaining hidden magnetic fields can only be accessed with the Hanle effect. Molecular lines are very convenient for applying the Hanle effect diagnostics thanks to the broad range of magnetic sensitivities in a narrow spectral region. With the UV version of the Zurich Imaging Polarimeter ZIMPOL II installed at the 45 cm telescope of the Istituto Ricerche Solari Locarno (IRSOL), we simultaneously observed intensity and linear polarization center-to-limb variations in two spectral regions containing the (0,0) and (1,1) bandheads of the CN B 2 {\Sigma} - X 2 {\Sigma} system. Here we present an analysis of these observations. We have implemented coherent scattering in molecular lines into a NLTE radiative transfer code. A two-step approach was used. First, we separately solved the statistical equilibrium equations and compute opacities and intensity while neglecting polariza- tion. Then we used these quantities as input for calculating scattering polarization and the Hanle effect. We have found that it is impossible to fit the intensity and polarization simultaneously at different limb angles in the frame- work of standard 1D modeling. The atmosphere models that provide correct intensity center-to-limb variations fail to fit linear polar- ization center-to-limb variations due to lacking radiation field anisotropy. We had to increase the anisotropy by means of a specially introduced free parameter. This allows us to successfully interpret our observations. We discuss possible reasons for underestimating the anisotropy in the 1D modeling.Comment: 15 pages, 10 figures, accepted for publication in Astronomy&Astrophysic

NLTE solar irradiance modeling with the COSI code

Context. The solar irradiance is known to change on time scales of minutes to decades, and it is suspected that its substantial fluctua- tions are partially responsible for climate variations. Aims. We are developing a solar atmosphere code that allows the physical modeling of the entire solar spectrum composed of quiet Sun and active regions. This code is a tool for modeling the variability of the solar irradiance and understanding its influence on Earth. Methods. We exploit further development of the radiative transfer code COSI that now incorporates the calculation of molecular lines. We validated COSI under the conditions of local thermodynamic equilibrium (LTE) against the synthetic spectra calculated with the ATLAS code. The synthetic solar spectra were also calculated in non-local thermodynamic equilibrium (NLTE) and compared to the available measured spectra. In doing so we have defined the main problems of the modeling, e.g., the lack of opacity in the UV part of the spectrum and the inconsistency in the calculations of the visible continuum level, and we describe a solution to these problems. Results. The improved version of COSI allows us to reach good agreement between the calculated and observed solar spectra as measured by SOLSTICE and SIM onboard the SORCE satellite and ATLAS 3 mission operated from the Space Shuttle. We find that NLTE effects are very important for the modeling of the solar spectrum even in the visual part of the spectrum and for its variability over the entire solar spectrum. In addition to the strong effect on the UV part of the spectrum, NLTE effects influence the concentration of the negative ion of hydrogen, which results in a significant change of the visible continuum level and the irradiance variability.Comment: 14 pages, 14 figures, accepted for publication in Astronomy&Astrophysic

Eclipses observed by LYRA - a sensitive tool to test the models for the solar irradiance

We analyze the light curves of the recent solar eclipses measured by the Herzberg channel (200-220 nm) of the Large Yield RAdiometer (LYRA) onboard PROBA-2. The measurements allow us to accurately retrieve the center- to-limb variations (CLV) of the solar brightness. The formation height of the radiation depends on the observing angle so the examination of the CLV provide information about a broad range of heights in the solar atmosphere. We employ the 1D NLTE radiative transfer COde for Solar Irradiance (COSI) to model the measured light curves and corresponding CLV dependencies. The modeling is used to test and constrain the existing 1D models of the solar atmosphere, e.g. the temperature structure of the photosphere and the treatment of the pseudo- continuum opacities in the Herzberg continuum range. We show that COSI can accurately reproduce not only the irradiance from the entire solar disk, but also the measured CLV. It hence can be used as a reliable tool for modeling the variability of the spectral solar irradiance.Comment: 19 pages, 9 figures, Solar Physic

Preparing for low surface brightness science with the Vera C. Rubin Observatory:Characterization of tidal features from mock images

Tidal features in the outskirts of galaxies yield unique information about their past interactions and are a key prediction of the hierarchical structure formation paradigm. The Vera C. Rubin Observatory is poised to deliver deep observations for potentially millions of objects with visible tidal features, but the inference of galaxy interaction histories from such features is not straightforward. Utilizing automated techniques and human visual classification in conjunction with realistic mock images produced using the NewHorizon cosmological simulation, we investigate the nature, frequency, and visibility of tidal features and debris across a range of environments and stellar masses. In our simulated sample, around 80 per cent of the flux in the tidal features around Milky Way or greater mass galaxies is detected at the 10-yr depth of the Legacy Survey of Space and Time (30-31 mag arcsec-2), falling to 60 per cent assuming a shallower final depth of 29.5 mag arcsec-2. The fraction of total flux found in tidal features increases towards higher masses, rising to 10 per cent for the most massive objects in our sample (M* ∼1011.5 M⊙). When observed at sufficient depth, such objects frequently exhibit many distinct tidal features with complex shapes. The interpretation and characterization of such features varies significantly with image depth and object orientation, introducing significant biases in their classification. Assuming the data reduction pipeline is properly optimized, we expect the Rubin Observatory to be capable of recovering much of the flux found in the outskirts of Milky Way mass galaxies, even at intermediate redshifts (z < 0.2)

Preparing for low surface brightness science with the Vera C. Rubin Observatory: characterisation of tidal features from mock images

Tidal features in the outskirts of galaxies yield unique information about their past interactions and are a key prediction of the hierarchical structure formation paradigm. The Vera C. Rubin Observatory is poised to deliver deep observations for potentially of millions of objects with visible tidal features, but the inference of galaxy interaction histories from such features is not straightforward. Utilising automated techniques and human visual classification in conjunction with realistic mock images produced using the NEWHORIZON cosmological simulation, we investigate the nature, frequency and visibility of tidal features and debris across a range of environments and stellar masses. In our simulated sample, around 80 per cent of the flux in the tidal features around Milky Way or greater mass galaxies is detected at the 10-year depth of the Legacy Survey of Space and Time (30-31 mag / sq. arcsec), falling to 60 per cent assuming a shallower final depth of 29.5 mag / sq. arcsec. The fraction of total flux found in tidal features increases towards higher masses, rising to 10 per cent for the most massive objects in our sample (M*~10^{11.5} Msun). When observed at sufficient depth, such objects frequently exhibit many distinct tidal features with complex shapes. The interpretation and characterisation of such features varies significantly with image depth and object orientation, introducing significant biases in their classification. Assuming the data reduction pipeline is properly optimised, we expect the Rubin Observatory to be capable of recovering much of the flux found in the outskirts of Milky Way mass galaxies, even at intermediate redshifts (z<0.2)

A method to estimate the effect of line blanketing in NLTE radiative transfer calculations

We present a method to estimate the contribution of line opacity to the total opacity as a function of wavelength. The estimated line-opacity function can then be used to simulate line-blanketing in NLTE radiative transfer calculations. Given a reference flux distribution (either observed or theoretical), our method allows to obtain a good estimate of the spectrum without the need for considering in detail all the millions of lines contributing to line blanketing. We applied the method to the spectra computed from a sample of photospheric models with effective temperatures $T_\mathrm{eff}$ $= 4200$, 5200 and 6200 K, $\log g$ $= 4.0$, 4.5, 5.0 and [A/H] $= 0.0$, -1.0, -2.0, taken from the NextGen database (Allard & Hauschildt [CITE]). The computed flux distributions agree quite well with the corresponding LTE line-blanketed NextGenfluxes when we introduce the estimated line-opacity contribution as a multiplicative factor of the continuum opacity in the radiative transfer calculations. In particular we discuss the importance of a correct estimate of the continuum flux, mainly in the UV, in the NLTE formation of the Cai

ARCO: a program for Automatic Reduction of CCD Observations

The huge amount of CCD photometric data collected by robotic telescopes (as for instance the ICE-T telescope to be installed at Dome-C) requires a fully automated approach to the reduction and analysis procedures. To this end, we are developing a pipeline, making use of IRAF, DAOPHOTII and tasks build up by us, which will enable to automatically extract differential magnitude time series. Both aperture and PSF photometry methods are used to build a proper sample of comparison stars, with the best method between the two being automatically selected. Finally, the search of rotational periods is performed. To test our package, we present results obtained using a 4-year dataset of a field in the Orion Nebula Cluster (ONC) containing about 350 stars (150 of which turned out to be periodic variables)