We investigate the dynamics and magnetic properties of the plasma, such as line-of-sight velocity (LOS), optical depth, vertical and horizontal magnetic fields, belonging to an erupted solar filament. The filament eruption was observed with the GREGOR Infrared Spectrograph (GRIS) at the 1.5-meter GREGOR telescope on 2016 July 3. Three consecutive full-Stokes slit-spectropolarimetric scans in the He I 10830 \r{A} spectral range were acquired. The Stokes I profiles were classified using the machine learning k-means algorithm and then inverted with different initial conditions using the HAZEL code. The erupting-filament material presents the following physical conditions: (1) ubiquitous upward motions with peak LOS velocities of ~73 km/s; (2) predominant large horizontal components of the magnetic field, on average, in the range of 173-254 G, whereas the vertical components of the fields are much lower, on average between 39-58 G; (3) optical depths in the range of 0.7-1.1. The average azimuth orientation of the field lines between two consecutive raster scans (<2.5 minutes) remained constant. The analyzed filament eruption belonged to the fast rising phase, with total velocities of about 124 km/s. The orientation of the magnetic field lines does not change from one raster scan to the other, indicating that the untwisting phase has not started yet. The untwisting seems to start about 15 min after the beginning of the filament eruption.

Observation from ground based telescopes are affected by the presence of the Earth atmosphere, which severely perturbs them. The use of adaptive optics techniques has allowed us to partly beat this limitation. However, image selection or post-facto image reconstruction methods are routinely needed to reach the diffraction limit of telescopes. Deep learning has been recently used to accelerate these image reconstructions. Currently, these deep neural networks are trained with supervision, so that standard deconvolution algorithms need to be applied a-priori to generate the training sets. Our aim is to propose an unsupervised method which can then be trained simply with observations and check it with data from the FastCam instrument. We use a neural model composed of three neural networks that are trained end-to-end by leveraging the linear image formation theory to construct a physically-motivated loss function. The analysis of the trained neural model shows that multiframe blind deconvolution can be trained self-supervisedly, i.e., using only observations. The output of the network are the corrected images and also estimations of the instantaneous wavefronts. The network model is of the order of 1000 times faster than applying standard deconvolution based on optimization. With some work, the model can bed used on real-time at the telescope.

Aims. We use high-resolution spectropolarimetric observations in the Ca ii 8542 A line obtained with the SST to study the magnetic field in solar spicules. Methods. The equations that result from the application of the Weak Field Approximation (WFA) to the radiative transfer equations are used to infer the LOS component of the magnetic field (BLOS). Two restrictive conditions are imposed on the Stokes I and V profiles at each pixel before they can be used in a Bayesian inversion to compute its BLOS. Results. The LOS magnetic field component has been inferred in six data sets totalling 448 spectral scans in the Ca ii 8542 A line and containing both active region and quiet Sun areas, with values of hundreds of G being abundantly inferred. There seems to be no difference, from the statistical point of view, between the magnetic field strength of spicules in the quiet Sun or near an active region. On the other hand, the BLOS distributions present smaller values on the disk than off-limb, a fact that can be explained by the effect of superposition on the chromosphere of on-disk structures. We find that spicules in the vicinity of a sunspot have a magnetic field polarity (i.e. north or south) equal to that of the sunspot. This paper also contains an analysis of the effect of off-limb overlapping structures on the observed Stokes I and V parameters and the BLOS obtained from the WFA. It is found that this value is equal to or smaller than the largest LOS magnetic field components of the two structures. In addition, using random BLOS, Doppler velocities and line intensities of these two structures leads in ~ 50% of the cases to Stokes I and V parameters unsuitable to be used with the WFA. Conclusions. Our results present a scarcity of LOS magnetic field components smaller than some 50 G, which must not be taken as evidence against the existence of these magnetic field strengths in spicules.

The nonlinear force-free field (NLFFF) modeling has been extensively used to infer the three-dimensional (3D) magnetic field in the solar corona. One of the assumptions in the NLFFF extrapolation is that the plasma beta is low, but this condition is considered to be incorrect in the photosphere. We examine direct measurements of the chromospheric magnetic field in two active regions through spectropolarimetric observations at He I 10830 A, which are compared with the potential fields and NLFFFs extrapolated from the photosphere. The comparisons allow quantitative estimation of the uncertainty in the NLFFF extrapolation from the photosphere. Our analysis shows that observed chromospheric magnetic field may have larger non-potentiality compared to the photospheric magnetic field. Moreover, the large non-potentiality in the chromospheric height may not be reproduced by the NLFFF extrapolation from the photospheric magnetic field. The magnitude of the underestimation of the non-potentiality at chromospheric heights may reach 30-40 degree in shear signed angle in some locations. This deviation may be caused by the non-force-freeness in the photosphere. Our study suggests the importance of the inclusion of measured chromospheric magnetic fields in the NLFFF modeling for the improvement of the coronal extrapolation.

Water is a key volatile that provides insights into the initial stages of planet formation. The low water abundances inferred from water observations toward low-mass protostellar objects may point to a rapid locking of water as ice by large dust grains during star and planet formation. However, little is known about the water vapor abundance in newly formed planet-forming disks. We aim to determine the water abundance in embedded Keplerian disks through spatially-resolved observations of H$_2^{18}$O lines to understand the evolution of water during star and planet formation. We present H$_2^{18}$O line observations with ALMA and NOEMA millimeter interferometers toward five young stellar objects. NOEMA observed the 3$_{1,3}$ - $2_{2,0}$ line (E$_{\rm up}$ = 203.7 K) while ALMA targeted the $4_{1,4}$ - $3_{2,1}$ line (E$_{\rm up}$ = 322.0 K). Water column densities are derived considering optically thin and thermalized emission. Our observations are sensitive to the emission from the known Keplerian disks around three out of the five Class I objects in the sample. No H$_2^{18}$O emission is detected toward any of our five Class I disks. We report upper limits to the integrated line intensities. The inferred water column densities in Class I disks are N < 10$^{15}$ cm$^{-2}$ on 100 au scales which include both disk and envelope. The upper limits imply a disk-averaged water abundance of $\lesssim 10^{-6}$ with respect to H$_2$ for Class I objects. After taking into account the physical structure of the disk, the upper limit to the water abundance averaged over the inner warm disk with $T>$ 100 K is between 10$^{-7}$ up to 10$^{-5}$. Water vapor is not abundant in warm protostellar envelopes around Class I protostars. Upper limits to the water vapor column densities in Class I disks are at least two orders magnitude lower than values found in Class 0 disk-like structures.

The analysis of waves in the visible side of the Sun allows the detection of active regions in the farside through local helioseismology techniques. The knowledge of the magnetism in the whole Sun, including the non-visible hemisphere, is fundamental for several space weather forecasting applications. Seismic identification of farside active regions is challenged by the reduced signal-to-noise, and only large and strong active regions can be reliable detected. Here we develop a new methodology to improve the identification of active region signatures in farside seismic maps. We have constructed a deep neural network that associates the farside seismic maps obtained from helioseismic holography with the probability of presence of active regions in the farside. The network has been trained with pairs of helioseismic phase shift maps and Helioseismic and Magnetic Imager magnetograms acquired half a solar rotation later, which were used as a proxy for the presence of active regions in the farside. The method has been validated using a set of artificial data, and it has also been applied to actual solar observations during the period of minimum activity of the solar cycle 24. Our approach shows a higher sensitivity to the presence of farside active regions than standard methods applied up to date. The neural network can significantly increase the number of detected farside active regions, and will potentially improve the application of farside seismology to space weather forecasting.

We analyze high-resolution spectropolarimetric observations of a flux emerging region (FER) in order to understand its magnetic and kinematic structure. Our spectropolarimetric observations in the He I 1083.0 nm spectral region of a FER are recorded with GRIS at the 1.5 m aperture GREGOR telescope. A Milne-Eddington based inversion code was employed to extract the photospheric information of the Si I spectral line, whereas the He I triplet line was analyzed with the Hazel inversion code, which takes into account the joint action of the Hanle and the Zeeman effect. The spectropolarimetric analysis of Si I line displays a complex magnetic structure near the vicinity of FER. Moreover, we find supersonic downflows of 40 km/sec appears near the footpoints of loops connecting two pores of opposite polarity, whereas a strong upflows of 22 km/sec appears near the apex of the loops. Furthermore, non-force-free field extrapolations were performed separately at two layers in order to understand the magnetic field topology of the FER. We determine, using extrapolations from the photosphere and the observed chromospheric magnetic field, that the average formation height of the He triplet line is 2 Mm from the solar surface. The reconstructed loops using photospheric extrapolations along an arch filament system have a maximum height of 10.5 Mm from the solar surface with a foot-points separation of 19 Mm, whereas the loops reconstructed using chromospheric extrapolations are around 8.4 Mm high from the solar surface with a foot-point separation of 16 Mm at the chromospheric height. The magnetic topology in the FER suggests the presence of small-scale loops beneath the large loops. Under suitable conditions, due to magnetic reconnection, these loops can trigger various heating events in the vicinity of the FER.

The three-dimensional (3-D) shape of a galaxy inevitably is tied to how it has formed and evolved and to its dark matter halo. Local extremely metal-poor galaxies (XMPs; defined as having an average gas-phase metallicity < 0.1 solar) are important objects for understanding galaxy evolution largely because they appear to be caught in the act of accreting gas from the cosmic web, and their 3-D shape may reflect this. Here we report on the 3-D shape of XMPs as inferred from their observed projected minor-to-major axial ratios using a hierarchical Bayesian inference model, which determines the likely shape and orientation of each galaxy while simultaneously inferring the average shape and dispersion. We selected a sample of 149 XMPs and divided it into three sub-samples according to physical size and found that (1) the stellar component of XMPs of all sizes tends to be triaxial, with an intermediate axis approx 0.7 times the longest axis and that (2) smaller XMPs tend to be relatively thicker, with the shortest axis going from approx 0.15 times the longest axis for the large galaxies to approx 0.4 for the small galaxies. We provide the inferred 3-D shape and inclination of the individual XMPs in electronic format. We show that our results for the intermediate axis are not clouded by a selection effect against face-on XMPs. We discuss how an intermediate axis significantly smaller than the longest axis may be produced by several mechanisms, including lopsided gas accretion, non-axisymmetric star formation, or coupling with an elongated dark matter halo. Large relative thickness may reflect slow rotation, stellar feedback, or recent gas accretion.

We analyze the effects of including $\Delta(1232)$ isobars in an equation of state (EoS) for cold, $\beta$-stable neutron star matter, employing relativistic nuclear mean field theory. The selected EoS reproduces the properties of nuclear matter and finite nuclei and, in the astrophysical context, allows for the presence of hyperons in neutron stars having masses larger than 2$M_{\odot}$. We find that the composition and structure of neutron stars is critically influenced by the addition of the $\Delta$ isobars, which allows us to constrain their interaction with the meson fields taking into account astrophysical information. Imposing that the EoS is stable and ensures the existence of 2$M_{\odot}$ neutron stars, as well as requiring agreement with data of $\Delta$ excitation in nuclei, we find that, in the absence of other mechanisms stiffening the EoS at high densities, the interaction of the $\Delta$ isobars with the sigma and omega meson fields must be at least 10\% stronger than that of the nucleons. Moreover, the neutron star moment of inertia turns out to be sensitive to the presence of $\Delta$ isobars, whereas the inclusion of $\Delta$ isobars in the EoS allows for smaller stellar radii and for a lower value of the tidal deformability consistent with the analysis of the GW170817 merger event.

Spectropolarimetric inversions are routinely used in the field of Solar Physics for the extraction of physical information from observations. The application to two-dimensional fields of view often requires the use of supercomputers with parallelized inversion codes. Even in this case, the computing time spent on the process is still very large. Our aim is to develop a new inversion code based on the application of convolutional neural networks that can quickly provide a three-dimensional cube of thermodynamical and magnetic properties from the interpretation of two-dimensional maps of Stokes profiles. We train two different architectures of fully convolutional neural networks. To this end, we use the synthetic Stokes profiles obtained from two snapshots of three-dimensional magneto-hydrodynamic numerical simulations of different structures of the solar atmosphere. We provide an extensive analysis of the new inversion technique, showing that it infers the thermodynamical and magnetic properties with a precision comparable to that of standard inversion techniques. However, it provides several key improvements: our method is around one million times faster, it returns a three-dimensional view of the physical properties of the region of interest in geometrical height, it provides quantities that cannot be obtained otherwise (pressure and Wilson depression) and the inferred properties are decontaminated from the blurring effect of instrumental point spread functions for free. The code is provided for free on a specific repository, with options for training and evaluation.

Our aim is to demonstrate the limitations of using a single component model to study the magnetic field of an active region filament. For that we have analyzed the polarimetric signals of the He I 10830 A multiplet acquired with the infrared spectrograph GRIS of the GREGOR telescope (Tenerife, Spain). After a first analysis of the general properties of the filament using Hazel under the assumption of a single component model atmosphere, in this second part we focus our attention on the observed Stokes profiles and the signatures which cannot be explained with this model. We have found an optically thick filament where the blue and the red components have the same sign in the linear polarization as an indication of radiative transfer effects. Moreover, the circular polarization signals inside the filament show the presence of strong magnetic field gradients. We have also shown that even a filament with such high absorption still shows signatures of the circular polarization generated by the magnetic field below the filament. The reason is that the absorption of the spectral line decays very quickly towards the wings, just where the circular polarization has a larger amplitude. In order to separate both contributions, we explore the possibility of a two component model but the inference becomes impossible to overcome as a high number of solutions is compatible with the observations.

The determination of the magnetic filed vector in solar filaments is possible by interpreting the Hanle and Zeeman effects in suitable chromospheric spectral lines like those of the He I multiplet at 10830 A. We study the vector magnetic field of an active region filament (NOAA 12087). Spectropolarimetric data of this active region was acquired with the GRIS instrument at the GREGOR telescope and studied simultaneously in the chromosphere with the He I 10830 A multiplet and in the photosphere with the Si I 10827 A line. As it is usual from previous studies, only a single component model is used to infer the magnetic properties of the filament. The results are put into a solar context with the help of the Solar Dynamic Observatory images. Some results clearly point out that a more complex inversion had to be done. Firstly, the Stokes $V$ map of He I does not show any clear signature of the presence of the filament. Secondly, the local azimuth map follows the same pattern than Stokes $V$ as if the polarity of Stokes $V$ were conditioning the inference to very different magnetic field even with similar linear polarization signals. This indication suggests that the Stokes $V$ could be dominated by the below magnetic field coming from the active region, and not, from the filament itself. Those and more evidences will be analyzed in depth and a more complex inversion will be attempted in the second part of this series.

We characterize for the first time the torus properties of an ultra-hard X-ray (14-195 keV) volume-limited (DL<40 Mpc) sample of 24 Seyfert (Sy) galaxies (BCS40 sample). The sample was selected from the Swift/BAT nine month catalog. We use high angular resolution nuclear infrared (IR) photometry and N-band spectroscopy, the CLUMPY torus models and a Bayesian tool to characterize the properties of the nuclear dust. In the case of the Sy1s we estimate the accretion disk contribution to the subarcsecond resolution nuclear IR SEDs (~0.4'') which is, on average, 46+-28, 23+-13 and 11+-5% in the J-, H- and K-bands, respectively. This indicates that the accretion disk templates that assume a steep fall for longer wavelengths than 1 micron might underestimate its contribution to the near-IR emission. Using both optical (broad vs narrow lines) and X-ray (unabsorbed vs absorbed) classifications, we compare the global posterior distribution of the torus model parameters. We confirm that Sy2s have larger values of the torus covering factor (CT~0.95) than Sy1s (CT~0.65) in our volume-limited Seyfert sample. These findings are independent of whether we use an optical or X-ray classification. We find that the torus covering factor remains essentially constant within the errors in our luminosity range and there is no clear dependence with the Eddington ratio. Finally, we find tentative evidence that even an ultra hard X-ray selection is missing a significant fraction of highly absorbed type 2 sources with very high covering factor tori.

A classic application of coronal seismology uses transverse oscillations of waveguides to obtain estimates of the magnetic field strength. The procedure requires information on the density of the structures. Often, it ignores the damping of the oscillations. We computed marginal posteriors for parameters such as the waveguide density; the density contrast; the transverse inhomogeneity length-scale; and the magnetic field strength, under the assumption that the oscillations can be modelled as standing magnetohydrodynamic (MHD) kink modes damped by resonant absorption. Our results show that the magnetic field strength can be properly inferred, even if the densities inside and outside the structure are largely unknown. Incorporating observational estimates of plasma density further constrains the obtained posteriors. The amount of information one is willing to include (a priori) for the density and the density contrast influences their corresponding posteriors, but very little the inferred magnetic field strength. The decision to include or leave out the information on the damping and the damping time-scales have a minimal impact on the obtained magnetic field strength. In contrast to the classic method which provides with numerical estimates with error bars or possible ranges of variation for the magnetic field strength, Bayesian methods offer the full distribution of plausibility over the considered range of possible values. The methods are applied to available datasets of observed transverse loop oscillations, can be extended to prominence fine structures or chromospheric spicules and implemented to propagating waves in addition to standing oscillations.

We present two recent parametrizations of the equation of state (FSU2R and FSU2H models) that reproduce the properties of nuclear matter and finite nuclei, fulfill constraints on high-density matter stemming from heavy-ion collisions, produce 2$M_{\odot}$ neutron stars, and generate neutron star radii below 13 km. Making use of these equations of state, cooling simulations for isolated neutron stars are performed. We find that two of the models studied, FSU2R (with nucleons) and, in particular, FSU2H (with nucleons and hyperons), show very good agreement with cooling observations, even without including nucleon pairing. This indicates that cooling observations are compatible with an equation of state that produces a soft nuclear symmetry energy and, thus, generates small neutron star radii. Nevertheless, both schemes produce cold isolated neutron stars with masses above $1.8 M_{\odot}$.

In this study we explore the diagnostic potential of the chromospheric Ca II line at 8542A for studying the magnetic and dynamic properties of solar filaments. We have acquired high spatial resolution spectropolarimetric observations in the Ca II 8542A line using the CRISP instrument at the Swedish 1-m Solar Telescope. We use the NICOLE inversion code to infer physical properties from observations of a solar filament. We discuss the validity of the results due to the assumption of hydrostatic equilibrium. We have used observations from other telescopes such as CHROTEL and SDO, in order to study large scale dynamics and the long term evolution of the filament. We show that the Ca II 8542A line encodes information of the temperature, line-of-sight velocity and magnetic field vector from the region where the filament is located. The current noise level only allow us to estimate an upper limit of 260G for the total magnetic field of the filament. Our study also reveals that if we only consider information from the aforementioned spectral line, the geometric height, the temperature and the density can be degenerated parameters outside the hydrostatic equilibrium approach.

Boyajian's star is an apparently normal main sequence F-type star with a very unusual light curve. The dipping activity of the star, discovered during the Kepler mission, presents deep, asymmetric, and aperiodic events. Here we present high resolution spectroscopic follow-up during some dimming events recorded post-Kepler observations, from ground-based telescopes. We analise data from the HERMES, HARPS-N and FIES spectrographs to characterise the stellar atmosphere and to put some constraints on the hypotheses that have appeared in the literature concerning the occulting elements. The star's magnetism, if existing, is not extreme. The spots on the surface, if present, would occupy 0.02% of the area, at most. The chromosphere, irrespective of the epoch of observation, is hotter than the values expected from radiative equilibrium, meaning that the star has some degree of activity. We find no clear evidence of the interstellar medium nor exocoments being responsible for the dimmings of the light curve. However, we detect at 1-2 sigma level, a decrease of the radial velocity of the star during the first dip recorded after the \emph{\emph{Kepler}} observations. We claim the presence of an optically thick object with likely inclined and high impact parameter orbits that produces the observed Rossiter-McLaughlin effect.

Based on the interstellar CO/H2 ratio, carbon monoxide-based censuses of protoplanetary disks in Lupus, sigma Orionis, and Chamaeleon I found no disks more massive than the minimum-mass solar nebula, which is inconsistent with the existence of exoplanets more massive than Jupiter. Observations and models are converging on the idea that ionization-driven chemistry depletes carbon monoxide in T-Tauri disks. Yet the extent of CO depletion depends on the incident flux of ionizing radiation, and some T-Tauri stars may have winds strong enough to shield their disks from cosmic rays. There is also a range of X-ray luminosities possible for a given stellar mass. Here we use a suite of chemical models, each with a different incident X-ray or cosmic-ray flux, to assess whether CO depletion is a typical outcome for T-Tauri disks. We find that CO dissociation in the outer disk is a robust result for realistic ionization rates, with abundance reductions between 70% and 99.99% over 2 Myr of evolution. Furthermore, after the initial dissociation epoch, the inner disk shows some recovery of the CO abundance from CO2 dissociation. In highly ionized disks, CO recovery in the inner disk combined with depletion in the outer disk creates a centrally peaked CO abundance distribution. The emitting area in rare CO isotopologues may be an indirect ionization indicator: in a cluster of disks with similar ages, those with the most compact CO isotopologue emission see the highest ionization rates.

The quality of images of the Sun obtained from the ground are severely limited by the perturbing effect of the turbulent Earth's atmosphere. The post-facto correction of the images to compensate for the presence of the atmosphere require the combination of high-order adaptive optics techniques, fast measurements to freeze the turbulent atmosphere and very time consuming blind deconvolution algorithms. Under mild seeing conditions, blind deconvolution algorithms can produce images of astonishing quality. They can be very competitive with those obtained from space, with the huge advantage of the flexibility of the instrumentation thanks to the direct access to the telescope. In this contribution we leverage deep learning techniques to significantly accelerate the blind deconvolution process and produce corrected images at a peak rate of ~100 images per second. We present two different architectures that produce excellent image corrections with noise suppression while maintaining the photometric properties of the images. As a consequence, polarimetric signals can be obtained with standard polarimetric modulation without any significant artifact. With the expected improvements in computer hardware and algorithms, we anticipate that on-site real-time correction of solar images will be possible in the near future.

The Chromospheric Lyman-Alpha Spectro-Polarimeter (CLASP) is a suborbital rocket experiment that on 3rd September 2015 measured the linear polarization produced by scattering processes in the hydrogen Ly-$\alpha$ line of the solar disk radiation, whose line-center photons stem from the chromosphere-corona transition region (TR). These unprecedented spectropolarimetric observations revealed an interesting surprise, namely that there is practically no center-to-limb variation (CLV) in the $Q/I$ line-center signals. Using an analytical model, we first show that the geometrical complexity of the corrugated surface that delineates the TR has a crucial impact on the CLV of the $Q/I$ and $U/I$ line-center signals. Secondly, we introduce a statistical description of the solar atmosphere based on a three-dimensional (3D) model derived from a state-of-the-art radiation magneto-hydrodynamic simulation. Each realization of the statistical ensemble is a 3D model characterized by a given degree of magnetization and corrugation of the TR, and for each such realization we solve the full 3D radiative transfer problem taking into account the impact of the CLASP instrument degradation on the calculated polarization signals. Finally, we apply the statistical inference method presented in a previous paper to show that the TR of the 3D model that produces the best agreement with the CLASP observations has a relatively weak magnetic field and a relatively high degree of corrugation. We emphasize that a suitable way to validate or refute numerical models of the upper solar chromosphere is by confronting calculations and observations of the scattering polarization in ultraviolet lines sensitive to the Hanle effect.

On 3rd September 2015, the Chromospheric Lyman-Alpha SpectroPolarimeter (CLASP) successfully measured the linear polarization produced by scattering processes in the hydrogen Lyman-$\alpha$ line of the solar disk radiation, revealing conspicuous spatial variations in the $Q/I$ and $U/I$ signals. Via the Hanle effect the line-center $Q/I$ and $U/I$ amplitudes encode information on the magnetic field of the chromosphere-corona transition region (TR), but they are also sensitive to the three-dimensional structure of this corrugated interface region. With the help of a simple line formation model, here we propose a statistical inference method for interpreting the Lyman-$\alpha$ line-center polarization observed by CLASP.

We perform cooling simulations for isolated neutron stars using recently developed equations of state for their core. The equations of state are obtained from new parametrizations of the FSU2 relativistic mean-field functional that reproduce the properties of nuclear matter and finite nuclei, while fulfilling the restrictions on high-density matter deduced from heavy-ion collisions, measurements of massive 2$M_{\odot}$ neutron stars, and neutron star radii below 13 km. We find that two of the models studied, FSU2R (with nucleons) and in particular FSU2H (with nucleons and hyperons), show very good agreement with cooling observations, even without including extensive nucleon pairing. This suggests that the cooling observations are more compatible with an equation of state that produces a soft nuclear symmetry energy and, hence, generates small neutron star radii. However, both models favor large stellar masses, above $1.8 M_{\odot}$, to explain the colder isolated neutron stars that have been observed, even if nucleon pairing is present.

We present MOLPOP-CEP, a universal line transfer code that allows the exact calculation of multi-level line emission from a slab with variable physical conditions for any arbitrary atom or molecule for which atomic data exist. The code includes error control to achieve any desired level of accuracy, providing full confidence in its results. Publicly available, MOLPOP-CEP employs our recently developed Coupled Escape Probability (CEP) technique, whose performance exceeds other exact methods by orders of magnitude. The program also offers the option of an approximate solution with different variants of the familiar escape probability method. As an illustration of the MOLPOP-CEP capabilities we present an exact calculation of the Spectral Line Energy Distribution (SLED) of the CO molecule and compare it with escape probability results. We find that the popular large-velocity gradient (LVG) approximation is unreliable at large CO column densities. Providing a solution of the multi-level line transfer problem at any prescribed level of accuracy, MOLPOP-CEP is removing any doubts about the validity of its final results.

We present high-precision spectro-polarimetric data with high spatial resolution (0.4$''$) of the very quiet Sun at 1.56$\mu$m obtained with the GREGOR telescope to shed some light on this complex magnetism. Half of our observed quiet-Sun region is better explained by magnetic substructure within the resolution element. However, we cannot distinguish whether this substructure comes from gradients of the physical parameters along the line of sight or from horizontal gradients (across the surface). In these pixels, a model with two magnetic components is preferred, and we find two distinct magnetic field populations. The population with the larger filling factor has very weak ($\sim$150 G) horizontal fields similar to those obtained in previous works. We demonstrate that the field vector of this population is not constrained by the observations, given the spatial resolution and polarimetric accuracy of our data. The topology of the other component with the smaller filling factor is constrained by the observations for field strengths above 250 G: we infer hG fields with inclinations and azimuth values compatible with an isotropic distribution. The filling factors are typically below 30\%. We also find that the flux of the two polarities is not balanced. From the other half of the observed quiet-Sun area $\sim$50\% are two-lobed Stokes $V$ profiles, meaning that 23\% of the field of view can be adequately explained with a single constant magnetic field embedded in a non-magnetic atmosphere. The magnetic field vector and filling factor are reliable inferred in only 50\% based on the regular profiles. Therefore, 12\% of the field of view harbour hG fields with filling factors typically below 30\%. At our present spatial resolution, 70\% of the pixels apparently are non-magnetised.

Chromospheric observations of sunspot umbrae offer an exceptional view of magneto-acoustic shock phenomena and the impact they have on the surrounding magnetically-dominated plasma. We employ simultaneous slit-based spectro-polarimetry and spectral imaging observations of the chromospheric He I 10830{\AA} and Ca II 8542{\AA} lines to examine fluctuations in the umbral magnetic field caused by the steepening of magneto-acoustic waves into umbral flashes. Following the application of modern inversion routines, we find evidence to support the scenario that umbral shock events cause expansion of the embedded magnetic field lines due to the increased adiabatic pressure. The large number statistics employed allow us to calculate the adiabatic index, gamma = 1.12 +/- 0.01, for chromospheric umbral locations. Examination of the vector magnetic field fluctuations perpendicular to the solar normal revealed changes up to ~200 G at the locations of umbral flashes. Such transversal magnetic field fluctuations have not been described before. Through comparisons with non-linear force-free field extrapolations, we find that the perturbations of the transverse field components are orientated in the same direction as the quiescent field geometries. This implies that magnetic field enhancements produced by umbral flashes are directed along the motion path of the developing shock, hence producing relatively small changes, up to a maximum of ~8 degrees, in the inclination and/or azimuthal directions of the magnetic field. Importantly, this work highlights that umbral flashes are able to modify the full vector magnetic field, with the detection of the weaker transverse magnetic field components made possible by high-resolution data combined with modern inversion routines.

The Helioseismic and Magnetic Imager (HMI) provides continuum images and magnetograms with a cadence better than one per minute. It has been continuously observing the Sun 24 hours a day for the past 7 years. The obvious trade-off between full disk observations and spatial resolution makes HMI not enough to analyze the smallest-scale events in the solar atmosphere. Our aim is to develop a new method to enhance HMI data, simultaneously deconvolving and super-resolving images and magnetograms. The resulting images will mimic observations with a diffraction-limited telescope twice the diameter of HMI. Our method, which we call Enhance, is based on two deep fully convolutional neural networks that input patches of HMI observations and output deconvolved and super-resolved data. The neural networks are trained on synthetic data obtained from simulations of the emergence of solar active regions. We have obtained deconvolved and supper-resolved HMI images. To solve this ill-defined problem with infinite solutions we have used a neural network approach to add prior information from the simulations. We test Enhance against Hinode data that has been degraded to a 28 cm diameter telescope showing very good consistency. The code is open source.

We present a photometric detection of the first brightness dips of the unique variable star KIC 8462852 since the end of the Kepler space mission in 2013 May. Our regular photometric surveillance started in October 2015, and a sequence of dipping began in 2017 May continuing on through the end of 2017, when the star was no longer visible from Earth. We distinguish four main 1-2.5% dips, named "Elsie," "Celeste," "Skara Brae," and "Angkor", which persist on timescales from several days to weeks. Our main results so far are: (i) there are no apparent changes of the stellar spectrum or polarization during the dips; (ii) the multiband photometry of the dips shows differential reddening favoring non-grey extinction. Therefore, our data are inconsistent with dip models that invoke optically thick material, but rather they are in-line with predictions for an occulter consisting primarily of ordinary dust, where much of the material must be optically thin with a size scale <<1um, and may also be consistent with models invoking variations intrinsic to the stellar photosphere. Notably, our data do not place constraints on the color of the longer-term "secular" dimming, which may be caused by independent processes, or probe different regimes of a single process.

We reexamine the equation of state for the nucleonic and hyperonic inner core of neutron stars that satisfies the 2$M_{\odot}$ observations as well as the recent determinations of stellar radii below 13 km, while fulfilling the saturation properties of nuclear matter and finite nuclei together with the constraints on the high-density nuclear pressure coming from heavy-ion collisions. The recent nucleonic FSU2R and hyperonic FSU2H models are fine tuned improving the density dependence of pure neutron matter at subsaturation densities. The corresponding nuclear matter properties at saturation, the symmetry energy and its slope turn out to be compatible with recent experimental and theoretical determinations. We obtain the mass, radius and composition of neutron stars for the two updated models and study the impact on these properties of the uncertainties in the hyperon-nucleon couplings estimated from hypernuclear data. We find that the onset of appearance of each hyperon strongly depends on the hyperon-nuclear uncertainties, whereas the maximum masses for neutron stars differ by at most 0.1 $M_{\odot}$, although a larger deviation should be expected tied to the lack of knowledge of the hyperon potentials at the high densities present in the center of $2 M_\odot$ stars. For easier use, we provide tables with the results from the FSU2R and FSU2H models for the equation of state and the neutron star mass-radius relation.

We investigate the properties of a sunspot light-bridge, focusing on the changes produced by the impact of a plasma blob ejected from a C-class flare. We observed a sunspot in active region NOAA 12544 using spectropolarimetric raster maps of the four Fe I lines around 15655 \AA\ with the GREGOR Infrared Spectrograph (GRIS), narrow-band intensity images sampling the Fe I 6173 \AA\ line with the GREGOR Fabry-P\'erot Interferometer (GFPI), and intensity broad band images in G-band and Ca II H band with the High-resolution Fast Imager (HiFI). All these instruments are located at the GREGOR telescope at the Observatorio del Teide, Tenerife, Spain. The data cover the time before, during, and after the flare event. The analysis is complemented with Atmospheric Imaging Assembly (AIA) and Helioseismic and Magnetic Imager (HMI) data from the Solar Dynamics Observatory (SDO). The physical parameters of the atmosphere at differents heights were inferred using spectral-line inversion techniques. We identify photospheric and chromospheric brightenings, heating events, and changes in the Stokes profiles associated to the flare eruption and the subsequent arrival of the plasma blob to the light bridge, after traveling along an active region loop. The measurements suggest that these phenomena are the result of reconnection events driven by the interaction of the plasma blob with the magnetic field topology of the light bridge.

$Context$. The thermal structure of the penumbra below its visible surface (i.e., $\tau_5 \ge 1$) has important implications for our present understanding of sunspots and their penumbrae: their brightness and energy transport, mode conversion of magneto-acoustic waves, sunspot seismology, and so forth. $Aims$. We aim at determining the thermal stratification in the layers immediately beneath the visible surface of the penumbra: $\tau_5 \in [1,3]$ ($\approx 70-80$ km below the visible continuum-forming layer). $Methods$. We analyzed spectropolarimetric data (i.e., Stokes profiles) in three Fe \textsc{i} lines located at 1565 nm observed with the GRIS instrument attached to the 1.5-meter solar telescope GREGOR. The data are corrected for the smearing effects of wide-angle scattered light and then subjected to an inversion code for the radiative transfer equation in order to retrieve, among others, the temperature as a function of optical depth $T(\tau_5)$. $Results$. We find that the temperature gradient below the visible surface of the penumbra is smaller than in the quiet Sun. This implies that in the region $\tau_5 \ge 1$ the penumbral temperature diverges from that of the quiet Sun. The same result is obtained when focusing only on the thermal structure below the surface of bright penumbral filaments. We interpret these results as evidence of a thick penumbra, whereby the magnetopause is not located near its visible surface. In addition, we find that the temperature gradient in bright penumbral filaments is lower than in granules. This can be explained in terms of the limited expansion of a hot upflow inside a penumbral filament relative to a granular upflow, as magnetic pressure and tension forces from the surrounding penumbral magnetic field hinder an expansion like this.

Many phenomena taking place in the solar photosphere are controlled by plasma motions. Although the line-of-sight component of the velocity can be estimated using the Doppler effect, we do not have direct spectroscopic access to the components that are perpendicular to the line-of-sight. These components are typically estimated using methods based on local correlation tracking. We have designed DeepVel, an end-to-end deep neural network that produces an estimation of the velocity at every single pixel and at every time step and at three different heights in the atmosphere from just two consecutive continuum images. We confront DeepVel with local correlation tracking, pointing out that they give very similar results in the time- and spatially-averaged cases. We use the network to study the evolution in height of the horizontal velocity field in fragmenting granules, supporting the buoyancy-braking mechanism for the formation of integranular lanes in these granules. We also show that DeepVel can capture very small vortices, so that we can potentially expand the scaling cascade of vortices to very small sizes and durations.

There is a thin transition region (TR) in the solar atmosphere where the temperature rises from 10,000 K in the chromosphere to millions of degrees in the corona. Little is known about the mechanisms that dominate this enigmatic region other than the magnetic field plays a key role. The magnetism of the TR can only be detected by polarimetric measurements of a few ultraviolet (UV) spectral lines, the Lyman-$\alpha$ line of neutral hydrogen at 121.6 nm (the strongest line of the solar UV spectrum) being of particular interest given its sensitivity to the Hanle effect (the magnetic-field-induced modification of the scattering line polarization). We report the discovery of linear polarization produced by scattering processes in the Lyman-$\alpha$ line, obtained with the Chromospheric Lyman-Alpha Spectro-Polarimeter (CLASP) rocket experiment. The Stokes profiles observed by CLASP in quiet regions of the solar disk show that the $Q/I$ and $U/I$ linear polarization signals are of the order of 0.1 % in the line core and up to a few percent in the nearby wings, and that both have conspicuous spatial variations with scales of $\sim 10$ arcsec. These observations help constrain theoretical models of the chromosphere-corona TR and extrapolations of the magnetic field from photospheric magnetograms. In fact, the observed spatial variation from disk to limb of polarization at the line core and wings already challenge the predictions from three-dimensional magnetohydrodynamical models of the upper solar chromosphere.

We present mid-infrared (MIR, 7.5-13.5 $\mu$m) imaging and spectroscopy observations obtained with the CanariCam (CC) instrument on the 10.4m Gran Telescopio CANARIAS for a sample of 20 nearby, MIR bright and X-ray luminous QSOs. We find that for the majority of QSOs the MIR emission is unresolved at angular scales nearly 0.3 arcsec, corresponding to physical scales $<600$ pc. We find that the higher-spatial resolution CC spectra have similar shapes to those obtained with Spitzer/IRS, and hence we can assume that the spectra are not heavily contaminated by extended emission in the host galaxy. We thus take advantage of the higher signal to noise Spitzer/IRS spectra, as a fair representation of the nuclear emission, to decompose it into a combination of active galactic nuclei (AGN), polycyclic aromatic hydrocarbon (PAH) and stellar components. In most cases the AGN is the dominant component, with a median contribution of 85 per cent of the continuum light at MIR (5-15 $\mu$m) within the IRS slit. This IR AGN emission is well reproduced by clumpy torus models. We find evidence for significant differences in the parameters that describe the dusty tori of QSOs when compared with the same parameters of Seyfert 1 and 2 nuclei. In particular, we find a lower number of clouds ($N_{0}<12$), steeper radial distribution of clouds ($q=1.5-3.0$), and clouds that are less optically thick ($\tau_{V}<100$) than in Seyfert 1, which could be attributed to dusty structures that have been partially evaporated and piled up by the higher radiation field in QSOs. We find that the combination of the angular width $\sigma_{torus}$, viewing angle $i$, and number of clouds along the equatorial line $N_{0}$, produces large escape probabilities ($P_{esc} > 2$ per cent) and low geometrical covering factors ($f_{2}<0.6$), as expected for AGN with broad lines in their optical spectra.

The purpose of this work is to better understand the confidence limits of the photospheric solar oxygen abundance derived from three-dimensional models using the forbidden [OI] line at 6300 \AA , including correlations with other parameters involved. We worked with a three-dimensional empirical model and two solar intensity atlases. We employed Bayesian inference as a tool to determine the most probable value for the solar oxygen abundance given the model chosen. We considered a number of error sources, such as uncertainties in the continuum derivation, in the wavelength calibration and in the abundance/strength of Ni. Our results shows correlations between the effects of several parameters employed in the derivation. The Bayesian analysis provides robust confidence limits taking into account all of these factors in a rigorous manner. We obtain that, given the empirical three-dimensional model and the atlas observations employed here, the most probable value for the solar oxygen abundance is $\log(\epsilon_O) = 8.86\pm0.04$. However, we note that this uncertainty does not consider possible sources of systematic errors due to the model choice.

Solar chromospheric fibrils, as observed in the core of strong chromospheric spectral lines, extend from photospheric field concentrations suggesting that they trace magnetic field lines. These images have been historically used as proxies of magnetic fields for many purposes. We use a Bayesian hierarchical model to analyze several tens of thousands of pixels in spectro-polarimetric chromospheric images of penumbrae and chromospheric fibrils. We compare the alignment between the field azimuth inferred from the linear polarization signals through the transverse Zeeman effect and the direction of the fibrils in the image. We conclude that, in the analyzed fields of view, fibrils are often well aligned with the magnetic field azimuth. Despite this alignment, the analysis also shows that there is a non-negligible dispersion. In penumbral filaments, we find a dispersion with a standard deviation of ~16 degrees, while this dispersion goes up to ~34 degrees in less magnetized regions.

We obtain a new equation of state for the nucleonic and hyperonic inner core of neutron stars that fulfills the 2$M_{\odot}$ observations as well as the recent determinations of stellar radii below 13 km. The nucleonic equation of state is obtained from a new parametrization of the FSU2 relativistic mean-field functional that satisfies these latest astrophysical constraints and, at the same time, reproduces the properties of nuclear matter and finite nuclei while fulfilling the restrictions on high-density matter deduced from heavy-ion collisions. On the one hand, the equation of state of neutron star matter is softened around saturation density, which increases the compactness of canonical neutron stars leading to stellar radii below 13 km. On the other hand, the equation of state is stiff enough at higher densities to fulfill the 2$M_{\odot}$ limit. By a slight modification of the parametrization, we also find that the constraints of 2$M_{\odot}$ neutron stars with radii around 13 km are satisfied when hyperons are considered. The inclusion of the high magnetic fields present in magnetars further stiffens the equation of state. Hyperonic magnetars with magnetic fields in the surface of $ \sim 10^{15}$ G and with values of $\sim 10^{18}$ G in the interior can reach maximum masses of 2$M_{\odot}$ with radii in the 12-13 km range.

Active regions are the most prominent manifestations of solar magnetic fields; their generation and dissipation are fundamental problems in solar physics. Light bridges are commonly present during sunspot decay, but a comprehensive picture of their role in the removal of photospheric magnetic field is still missing. We study the three dimensional configuration of a sunspot and in particular its light bridge during one of the last stages of its decay. We present the magnetic and thermodynamical stratification inferred from full Stokes inversions of the photospheric Si I 10827 \AA\ and Ca I 10839 \AA\ lines obtained with the GREGOR Infrared Spectrograph of the GREGOR telescope at Observatorio del Teide, Tenerife, Spain. The analysis is complemented by a study of continuum images covering the disk passage of the active region, which are provided by the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory. The sunspot shows a light bridge with penumbral continuum intensity that separates the central umbra from a smaller umbra. We find that in this region the magnetic field lines form a canopy with lower magnetic field strength in the inner part. The photospheric light bridge is dominated by gas pressure (high-$\beta$), as opposed to the surrounding umbra where the magnetic pressure is higher. A convective flow is observed in the light bridge. This flow is able to bend the magnetic field lines and to produce field reversals. The field lines close above the light bridge and become as vertical and strong as in the surrounding umbra. We conclude that it develops because of two highly magnetized regions which come closer during the sunspot evolution.

Context. Investigations of the magnetism of the quiet Sun are hindered by extremely weak polarization signals in Fraunhofer spectral lines. Photon noise, straylight, and the systematically different sensitivity of the Zeeman effect to longitudinal and transversal magnetic fields result in controversial results in terms of the strength and angular distribution of the magnetic field vector. Aims. The information content of Stokes measurements close to the diffraction limit of the 1.5 m GREGOR telescope is analyzed. We took the effects of spatial straylight and photon noise into account. Methods. Highly sensitive full Stokes measurements of a quiet-Sun region at disk center in the deep photospheric Fe I lines in the 1.56 {\mu}m region were obtained with the infrared spectropolarimeter GRIS at the GREGOR telescope. Noise statistics and Stokes V asymmetries were analyzed and compared to a similar data set of the Hinode spectropolarimeter (SOT/SP). Simple diagnostics based directly on the shape and strength of the profiles were applied to the GRIS data. We made use of the magnetic line ratio technique, which was tested against MHD simulations. Results. About 80% of the GRIS spectra of a very quiet solar region show polarimetric signals above a 3{\sigma} level. Area and amplitude asymmetries agree well with small-scale surface dynamo MHD simulations. The magnetic line ratio analysis reveals ubiquitous magnetic regions in the ten to hundred Gauss range with some concentrations of kilo-Gauss fields. Conclusions. The GRIS spectropolarimetric data at a spatial resolution of 0.40" are so far unique in the combination of high spatial resolution scans and high magnetic field sensitivity. Nevertheless, the unavoidable effect of spatial straylight and the resulting dilution of the weak Stokes profiles means that inversion techniques still bear a high risk of misinterpretating the data.

Stellar prominences detected in rapidly rotating stars serve as probes of the magnetism in the corona of cool stars. We have synthesized the temporal evolution of the Stokes profiles generated in the He I 10830 and 5876 A triplets during the rotation of a prominence around a star. The synthesis was performed with the HAZEL code using a cloud model in which the prominence is characterized by a slab located at a fixed latitude and height. It accounts for the scattering polarization and Zeeman and Hanle effects. Several cases with different prominence magnetic field strengths and orientations have been analyzed. The results show an emission feature that drifts across the profile while the prominence is out of the stellar disk. When the prominence eclipses the star, the intensity profile shows an absorption. The scattering induced by the prominence generates linear polarization signals in Stokes Q and U profiles, which are modified by the Hanle effect when a magnetic field is present. Due to the Zeeman effect, Stokes V profiles show a signal with very low amplitude when the magnetic field along the line-of-sight is different from zero. The estimated linear polarization signals could potentially be detected with the future spectropolarimeter MIRADAS, to be attached to GTC telescope.

Context. Ellerman bombs (EBs) are short-lived emission features, characterized by extended wing emission in hydrogen Balmer lines. Until now, no distinct signature of EBs has been found in the He I 10830 {\AA} line, and conclusive observations of EBs in He I D 3 have never been reported. Aims. We aim to study the signature of EBs in neutral helium triplet lines. Methods. The observations consist of 10 consecutive SST/TRIPPEL raster scans close to the limb, featuring the H$\beta$, He I D3 and He I 10830 {\AA} spectral regions. We also obtained raster scans with IRIS and make use of the SDO/AIA 1700 {\AA} channel. We use Hazel to invert the neutral helium triplet lines. Results. Three EBs in our data show distinct emission signatures in neutral helium triplet lines, most prominently visible in the He I D3 line. The helium lines have two components: a broad and blue-shifted emission component associated with the EB, and a narrower absorption component formed in the overlying chromosphere. One of the EBs in our data shows evidence of strong velocity gradients in its emission component. The emission component of the other two EBs could be fitted using a constant slab. Our analysis hints towards thermal Doppler motions having a large contribution to the broadening for helium and IRIS lines. We conclude that the EBs must have high temperatures to exhibit emission signals in neutral helium triplet lines. An order of magnitude estimate places our observed EBs in the range of $T\sim 2\cdot 10^4-10^5$ K.

The Mg II h&k lines are powerful diagnostics for studying the solar chromosphere. They have become particularly popular with the launch of the IRIS satellite, and a number of studies that include these lines have lead to great progress in understanding chromospheric heating, in many cases thanks to the support from 3D MHD simulations. In this study we utilize another approach to analyze observations: non-LTE inversions of the Mg II h&k and UV triplet lines including the effects of partial redistribution. Our inversion code attempts to construct a model atmosphere that is compatible with the observed spectra. We have assessed the capabilities and limitations of the inversions using the FALC atmosphere and a snapshot from a 3D radiation-MHD simulation. We find that Mg II h&k allow reconstructing a model atmosphere from the middle photosphere to the transition region. We have also explored the capabilities of a multi-line/multi-atom setup, including the Mg II h&k, the Ca II 854.2 nm and the Fe I 630.25 lines to recover the full stratification of physical parameters, including the magnetic field vector, from the photosphere to the chromosphere. Finally, we present the first inversions of observed IRIS spectra from quiet-Sun, plage and sunspot, with very promising results.

Some models for the topology of the magnetic field in sunspot penumbrae predict the existence of field-free or dynamically weak-field regions in the deep Photosphere. To confirm or rule out the existence of weak-field regions in the deepest photospheric layers of the penumbra. The magnetic field at $\log\tau_5=0$ is investigated by means of inversions of spectropolarimetric data of two different sunspots located very close to disk center with a spatial resolution of approximately 0.4-0.45 arcsec. The data have been recorded using the GRIS instrument attached to the 1.5-meters GREGOR solar telescope at El Teide observatory. It includes three Fe I lines around 1565 nm, whose sensitivity to the magnetic field peaks at half a pressure-scale-height deeper than the sensitivity of the widely used Fe I spectral line pair at 630 nm. Prior to the inversion, the data is corrected for the effects of scattered light using a deconvolution method with several point spread functions. At $\log\tau_5=0$ we find no evidence for the existence of regions with dynamically weak ($B<500$ Gauss) magnetic fields in sunspot penumbrae. This result is much more reliable than previous investigations done with Fe I lines at 630 nm. Moreover, the result is independent of the number of nodes employed in the inversion, and also independent of the point spread function used to deconvolve the data, and does not depend on the amount of straylight (i.e. wide-angle scattered light) considered.

The origin of the unification model for Active Galactic Nuclei (AGN) was the detection of broad hydrogen recombination lines in the optical polarized spectrum of the Seyfert 2 galaxy (Sy2) NGC 1068. Since then, a search for the hidden broad-line region (HBLR) of nearby Sy2s started, but polarized broad lines have only been detected in 30-40% of the nearby Sy2s observed to date. Here we present new VLT/FORS2 optical spectropolarimetry of a sample of 15 Sy2s, including Compton-thin and Compton-thick sources. The sample includes six galaxies without previously published spectropolarimetry, some of them normally treated as non-hidden BLR (NHBLR) objects in the literature, four classified as NHBLR, and five as HBLR based on previous data. We report >=4{\sigma} detections of a HBLR in 11 of these galaxies (73% of the sample) and a tentative detection in NGC 5793, which is Compton-thick according to the analysis of X-ray data performed here. Our results confirm that at least some NHBLRs are misclassified, bringing previous publications reporting differences between HBLR and NHBLR objects into question. We detect broad H{\alpha} and H{\beta} components in polarized light for 10 targets, and just broad H{\alpha} for NGC 5793 and NGC 6300, with line widths ranging between 2100 and 9600 km/s. High bolometric luminosities and low column densities are associated with higher polarization degrees, but not necessarily with the detection of the scattered broad components.

Solar tornadoes are dark vertical filamentary structures observed in the extreme ultraviolet associated with prominence legs and filament barbs. Their true nature and relationship to prominences requires understanding their magnetic structure and dynamic properties. Recently, a controversy has arisen: is the magnetic field organized forming vertical, helical structures or is it dominantly horizontal? And concerning their dynamics, are tornadoes really rotating or is it just a visual illusion? Here, we analyze four consecutive spectropolarimetric scans of a prominence hosting tornadoes on its legs which help us shed some light on their magnetic and dynamical properties. We show that the magnetic field is very smooth in all the prominence, probably an intrinsic property of the coronal field. The prominence legs have vertical helical fields that show slow temporal variation probably related to the motion of the fibrils. Concerning the dynamics, we argue that 1) if rotation exists, it is intermittent, lasting no more than one hour, and 2) the observed velocity pattern is also consistent with an oscillatory velocity pattern (waves).

Polar faculae are bright features that can be detected in solar limb observations and they are related to magnetic field concentrations. Although there is a large number of works studying them, some questions about their nature as their magnetic properties at different heights are still open. Thus, we aim to improve the understanding of solar polar faculae. In that sense, we infer the vertical stratification of the temperature, gas pressure, line of sight velocity and magnetic field vector of polar faculae regions. We performed inversions of the Stokes profiles observed with Hinode/SP after removing the stray light contamination produced by the spatial point spread function of the telescope. Moreover, after solving the azimuth ambiguity, we transform the magnetic field vector to local solar coordinates. The obtained results reveal that the polar faculae are constituted by hot plasma with low line of sight velocities and single polarity magnetic fields in the kilogauss range that are nearly perpendicular to the solar surface. We also found that the spatial location of these magnetic fields is slightly shifted respect to the continuum observations towards the disc centre. We believe that this is due to the hot wall effect that allows detecting photons that come from deeper layers located closer to the solar limb.

Quantitative thermodynamical, dynamical and magnetic properties of the solar and stellar plasmas are obtained by interpreting their emergent non-polarized and polarized spectrum. This inference requires the selection of a set of spectral lines particularly sensitive to the physical conditions in the plasma and a suitable parametric model of the solar/stellar atmosphere. Nonlinear inversion codes are then used to fit the model to the observations. However, the presence of systematic effects like nearby or blended spectral lines, telluric absorption or incorrect correction of the continuum, among others, can strongly affect the results. We present an extension to current inversion codes that can deal with these effects in a transparent way. The resulting algorithm is very simple and can be applied to any existing inversion code with the addition of a few lines of code as an extra step in each iteration.

Recent spectropolarimetric observations of active region filaments have revealed polarization profiles with signatures typical of the strong field Zeeman regime. The conspicuous absence in those observations of scattering polarization and Hanle effect signatures was then pointed out by some authors. This was interpreted either as a signature of mixed "turbulent" field components or as a result of optical thickness. In this article, we present a natural scenario to explain these Zeeman-only spectro-polarimetric observations of active region filaments. We propose a two-component model, one on top of the other. Both components have horizontal fields, the azimuth difference between them being close to 90 degrees. The component that lies lower in the atmosphere is permeated by a strong field of the order of 600 G, while the upper component has much weaker fields, of the order of 10 G. The ensuing scattering polarization signatures of the individual components have opposite signs, so that its combination along the line of sight reduces --and even can cancel out-- the Hanle signatures, giving rise to an apparent only-Zeeman profile. This model is also applicable to other chromospheric structures seen in absorption above active regions.

We present the distributions of geometrical covering factors of active galactic nuclei (AGNs) dusty tori (f2) using an X-ray selected complete sample of 227 AGN drawn from the Bright Ultra-hard XMM-Newton Survey. The AGN have z from 0.05 to 1.7, 2-10 keV luminosities between 10^42 and 10^46 erg/s and Compton-thin X-ray absorption. Employing data from UKIDSS, 2MASS and the Wide-field Infrared Survey Explorer in a previous work we determined the rest-frame 1-20 microns continuum emission from the torus which we model here with the clumpy torus models of Nenkova et al. Optically classified type 1 and type 2 AGN are intrinsically different, with type 2 AGN having on average tori with higher f2 than type 1 AGN. Nevertheless, ~20 per cent of type 1 AGN have tori with large covering factors while ~23-28 per cent of type 2 AGN have tori with small covering factors. Low f2 are preferred at high AGN luminosities, as postulated by simple receding torus models, although for type 2 AGN the effect is certainly small. f2 increases with the X-ray column density, which implies that dust extinction an X-ray absorption takes place in material that shares an overall geometry and most likely belongs to the same structure, the putative torus. Based on our results, the viewing angle, AGN luminosity and also f2 determine the optical appearance of an AGN and control the shape of the rest-frame ~1-20 microns nuclear continuum emission. Thus, the torus geometrical covering factor is a key ingredient of unification schemes.

We show that it is possible to measure polarization with a polarimeter that gets rid of the seeing while still measuring at a frequency well below that of the seeing. We study a standard polarimeter made of two retarders and a beamsplitter. The retarders are modulated at $\sim 500$ Hz, a frequency comparable to that of the variations of the refraction index in the Earth atmosphere, what is usually termed as seeing in astronomical observations. However, we assume that the camera is slow, so that our measurements are time integrations of these modulated signals. In order to recover the time variation of the seeing and obtain the Stokes parameters, we use the theory of compressed sensing to solve the demodulation by impose a sparsity constraint on the Fourier coefficients of the seeing. We demonstrate the feasibility of this sub-Nyquist polarimeter using numerical simulations, both in the case without noise and with noise. We show that a sensible modulation scheme is obtained by randomly changing the fast axis of the modulators or their retardances in specific ways. We finally demonstrate that the value of the Stokes parameters can be recovered with great precision at almost maximum efficiency, although it slightly degrades when the signal-to-noise ratio of the observations increase, a consequence of the multiplexing under the presence of photon noise.

We present an analysis of the nuclear infrared (IR, 1.6 to 18 $\mu$m) emission of the ultraluminous IR galaxy UGC 5101 to derive the properties of its active galactic nucleus (AGN) and its obscuring material. We use new mid-IR high angular resolution ($0.3-0.5$ arcsec) imaging using the Si-2 filter ($\lambda_{C}=8.7\, \mu$m) and $7.5-13$ $\mu$m spectroscopy taken with CanariCam (CC) on the 10.4m Gran Telescopio CANARIAS. We also use archival HST/NICMOS and Subaru/COMICS imaging and Spitzer/IRS spectroscopy. We estimate the near- and mid-IR unresolved nuclear emission by modelling the imaging data with GALFIT. We decompose the Spitzer/IRS and CC spectra using a power-law component, which represents the emission due to dust heated by the AGN, and a starburst component, both affected by foreground extinction. We model the resulting unresolved near- and mid-IR, and the starburst subtracted CC spectrum with the CLUMPY torus models of Nenkova et al. The derived geometrical properties of the torus, including the large covering factor and the high foreground extinction needed to reproduce the deep $9.7\, \mu$m silicate feature, are consistent with the lack of strong AGN signatures in the optical. We derive an AGN bolometric luminosity $L_{bol}\sim1.9\times10^{45}\,$erg s$^{-1}$ that is in good agreement with other estimates in the literature.

Aims. To develop a fully Bayesian least squares deconvolution (LSD) that can be applied to the reliable detection of magnetic signals in noise-limited stellar spectropolarimetric observations using multiline techniques. Methods. We consider LSD under the Bayesian framework and we introduce a flexible Gaussian Process (GP) prior for the LSD profile. This prior allows the result to automatically adapt to the presence of signal. We exploit several linear algebra identities to accelerate the calculations. The final algorithm can deal with thousands of spectral lines in a few seconds. Results. We demonstrate the reliability of the method with synthetic experiments and we apply it to real spectropolarimetric observations of magnetic stars. We are able to recover the magnetic signals using a small number of spectral lines, together with the uncertainty at each velocity bin. This allows the user to consider if the detected signal is reliable. The code to compute the Bayesian LSD profile is freely available.

The spatial variation of physical quantities, such as the mass density, across solar atmospheric waveguides governs the timescales and spatial scales for wave damping and energy dissipation. The direct measurement of the spatial distribution of density, however, is difficult and indirect seismology inversion methods have been suggested as an alternative. We applied Bayesian inference, model comparison, and model-averaging techniques to the inference of the cross-field density structuring in solar magnetic waveguides using information on periods and damping times for resonantly damped magnetohydrodynamic (MHD) transverse kink oscillations. Three commonly employed alternative profiles were used to model the variation of the mass density across the waveguide boundary. Parameter inference enabled us to obtain information on physical quantities such as the Alfv\'en travel time, the density contrast, and the transverse inhomogeneity length scale. The inference results from alternative density models were compared and their differences quantified. Then, the relative plausibility of the considered models was assessed by performing model comparison. Our results indicate that the evidence in favor of any of the three models is minimal, unless the oscillations are strongly damped. In such a circumstance, the application of model-averaging techniques enables the computation of an evidence-weighted inference that takes into account the plausibility of each model in the calculation of a combined inversion for the unknown physical parameters.

Observations of the Sun from the Earth are always limited by the presence of the atmosphere, which strongly disturbs the images. A solution to this problem is to place the telescopes in space satellites, which produce observations without any (or limited) atmospheric aberrations. However, even though the images from space are not affected by atmospheric seeing, the optical properties of the instruments still limit the observations. In the case of diffraction limited observations, the PSF establishes the maximum allowed spatial resolution, defined as the distance between two nearby structures that can be properly distinguished. In addition, the shape of the PSF induce a dispersion of the light from different parts of the image, leading to what is commonly termed as stray light or dispersed light. This effect produces that light observed in a spatial location at the focal plane is a combination of the light emitted in the object at relatively distant spatial locations. We aim to correct the effect produced by the telescope's PSF using a deconvolution method, and we decided to apply the code on Hinode/SP quiet Sun observations. We analyze the validity of the deconvolution process with noisy data and we infer the physical properties of quiet Sun magnetic elements after the deconvolution process.

Proving the magnetic configuration of solar spicules has hitherto been difficult due to the lack of spatial resolution and image stability during off-limb ground-based observations. We report spectropolarimetric observations of spicules taken in the He I 1083 nm spectral region with the Tenerife Infrared Polarimeter II at the German Vacuum Tower Telescope of the Observatorio del Teide (Tenerife; Canary Islands; Spain). The data provide the variation with geometrical height of the Stokes I, Q, U, and V profiles whose encoded information allows the determination of the magnetic field vector by means of the HAZEL inversion code. The inferred results show that the average magnetic field strength at the base of solar spicules is about 80 gauss and then it decreases rapidly with height to about 30 gauss at a height of 3000 km above the visible solar surface. Moreover, the magnetic field vector is close to vertical at the base of the chromosphere and has mid inclinations (about 50 degree) above 2 Mm height.

A properly calibrated longitudinal magnetograph is an instrument that measures circular polarization and gives an estimation of the magnetic flux density in each observed resolution element. This usually constitutes a lower bound of the field strength in the resolution element, given that it can be made arbitrarily large as long as it occupies a proportionally smaller area of the resolution element and/or becomes more transversal to the observer and still produce the same magnetic signal. Yet, we know that arbitrarily stronger fields are less likely --hG fields are more probable than kG fields, with fields above several kG virtually absent-- and we may even have partial information about its angular distribution. Based on a set of sensible considerations, we derive simple formulae based on a Bayesian analysis to give an improved estimation of the magnetic field strength for magnetographs.

Inversion codes are numerical tools used for the inference of physical properties from the observations. Despite their success, the quality of current spectropolarimetric observations and those expected in the near future presents a challenge to current inversion codes. The pixel-by-pixel strategy of inverting spectropolarimetric data that we currently utilize needs to be surpassed and improved. The inverted physical parameters have to take into account the spatial correlation that is present in the data and that contains valuable physical information. We utilize the concept of sparsity or compressibility to develop an new generation of inversion codes for the Stokes parameters. The inversion code uses numerical optimization techniques based on the idea of proximal algorithms to impose sparsity. In so doing, we allow for the first time to exploit the presence of spatial correlation on the maps of physical parameters. Sparsity also regularizes the solution by reducing the number of unknowns. We compare the results of the new inversion code with pixel-by-pixel inversions, demonstrating the increase in robustness of the solution. We also show how the method can easily compensate for the effect of the telescope point spread function, producing solutions with an enhanced contrast.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

We present an exceptional data set acquired with the Vacuum Tower Telescope (Tenerife, Spain) covering the pre-flare, flare, and post-flare stages of an M3.2 flare. The full Stokes spectropolarimetric observations were recorded with the Tenerife Infrared Polarimeter in the He I 1083.0 nm spectral region. The object under study was active region NOAA 11748 on 2013 May 17. During the flare the chomospheric He I 1083.0 nm intensity goes strongly into emission. However, the nearby photospheric Si I 1082.7 nm spectral line profile only gets shallower and stays in absorption. Linear polarization (Stokes Q and U) is detected in all lines of the He I triplet during the flare. Moreover, the circular polarization (Stokes V) is dominant during the flare, being the blue component of the He I triplet much stronger than the red component, and both are stronger than the Si I Stokes V profile. The Si I inversions reveal enormous changes of the photospheric magnetic field during the flare. Before the flare magnetic field concentrations of up to 1500 G are inferred. During the flare the magnetic field strength globally decreases and in some cases it is even absent. After the flare the magnetic field recovers its strength and initial configuration.

We present results from the fitting of infrared (IR) spectral energy distributions of 21 active galactic nuclei (AGN) with clumpy torus models. We compiled high spatial resolution ($\sim 0.3$--$0.7$ arcsec) mid-IR $N$-band spectroscopy, $Q$-band imaging and nuclear near- and mid-IR photometry from the literature. Combining these nuclear near- and mid-IR observations, far-IR photometry and clumpy torus models, enables us to put constraints on the torus properties and geometry. We divide the sample into three types according to the broad line region (BLR) properties; type-1s, type-2s with scattered or hidden broad line region (HBLR) previously observed, and type-2s without any published HBLR signature (NHBLR). Comparing the torus model parameters gives us the first quantitative torus geometrical view for each subgroup. We find that NHBLR AGN have smaller torus opening angles and larger covering factors than those of HBLR AGN. This suggests that the chance to observe scattered (polarized) flux from the BLR in NHBLR could be reduced by the dual effects of (a) less scattering medium due to the reduced scattering volume given the small torus opening angle and (b) the increased torus obscuration between the observer and the scattering region. These effects give a reasonable explanation for the lack of observed HBLR in some type-2 AGN.

Solar prominences are clouds of cool plasma levitating above the solar surface and insulated from the million-degree corona by magnetic fields. They form in regions of complex magnetic topology, characterized by non-potential fields, which can evolve abruptly, disintegrating the prominence and ejecting magnetized material into the heliosphere. However, their physics is not yet fully understood because mapping such complex magnetic configurations and their evolution is extremely challenging, and must often be guessed by proxy from photometric observations.Using state-of-the-art spectro-polarimetric data, we reconstruct the structure of the magnetic field in a prominence. We find that prominence feet harbor helical magnetic fields connecting the prominence to the solar surface below.

This paper shows the first synthetic tomography of the quiet solar chromosphere formed by spatial maps of scattering polarization. It has been calculated for the CaII 8498, 8542 and 3934 A lines by solving the NLTE (non-local thermodynamical equilibrium) RT (radiative transfer) problem of the second kind in a 3D atmosphere model obtained from realistic MHD (magneto-hydrodynamical) simulations. Maps of circular polarization were calculated neglecting atomic polarization. Our investigation focuses on the linear polarization signals induced by kinematics, radiation field anisotropy and Hanle effect in forward-scattering geometry. Thus, instead of considering slit profiles at the limb as normally done in the study of the second solar spectrum, we synthetize and analyze spatial maps of polarization at disk center. It allows us to understand the spatial signatures of dynamics and magnetic field in the linear polarization for discriminating them observationally. Our results suggest new ideas for chromospheric diagnosis that will be developed throughout a serie of papers. In particular, Hanle Polarity Inversion Lines and dynamic Hanle diagrams are two concepts introduced in the present work. We find that chromospheric dynamics and magnetic field topology create spatial fingerprints in the polarization maps that trace the dynamic situation of the plasma and the magnetic field. Based on such spatial features we reconstruct the magnetic field intensity in the middle chromosphere along grooves of null linear polarization. We finally address the problems of diagnosing Hanle saturation and kinematic amplification of scattering signals using Hanle diagrams.

We carried out high-sensitivity spectropolarimetric observations of the central star of the Red Rectangle proto-planetary nebula with the aim of constraining the mechanism that gives its biconical shape. The stellar light of the central binary system is linearly polarised since it is scattered on the dust particles of the nebula. Surprisingly, the linear polarisation in the continuum is aligned with one of the spikes of the biconical outflow. Also, the observed Balmer lines as well as the Ca II K lines are polarised. These observational constraints are used to confirm or reject current theoretical models for the shaping mechanism of the Red Rectangle. We propose that the observed polarisation is very unlikely generated by a uniform biconical stellar wind. Also, the hypothesis of a precessing jet does not completely match the observations since it will require a jet aperture larger than that of the nebula.

Standard statistical analysis of the magnetic properties of the quiet Sun rely on simple histograms of quantities inferred from maximum-likelihood estimations. Because of the inherent degeneracies, either intrinsic or induced by the noise, this approach is not optimal and can lead to highly biased results. We carry out a meta-analysis of the magnetism of the quiet Sun from Hinode observations using a hierarchical probabilistic method. This model allows us to infer the statistical properties of the magnetic field vector over the observed field-of-view consistently taking into account the uncertainties in each pixel due to noise and degeneracies. Our results point out that the magnetic fields are very weak, below 275 G with 95% credibility, with a slight preference for horizontal fields, although the distribution is not far from a quasi-isotropic distribution.

With the advent of a new generation of solar telescopes and instrumentation, the interpretation of chromospheric observations (in particular, spectro-polarimetry) requires new, suitable diagnostic tools. This paper describes a new code, NICOLE, that has been designed for Stokes non-LTE radiative transfer, both for synthesis and inversion of spectral lines and Zeeman-induced polarization profiles, spanning a wide range of atmospheric heights, from the photosphere to the chromosphere. The code fosters a number of unique features and capabilities and has been built from scratch with a powerful parallelization scheme that makes it suitable for application on massive datasets using large supercomputers. The source code is being publicly released, with the idea of facilitating future branching by other groups to augment its capabilities.

Depolarizing collisions are elastic or quasielastic collisions that equalize the populations and destroy the coherence between the magnetic sublevels of atomic levels. In astrophysical plasmas, the main depolarizing collider is neutral hydrogen. We consider depolarizing rates on the lowest levels of neutral and singly ionized alkaly-earths Mg I, Sr I, Ba I, Mg II, Ca II, and Ba II, due to collisions with H. We compute ab initio potential curves of the atom-H system and solve the quantum mechanical dynamics. From the scattering amplitudes we calculate the depolarizing rates for Maxwellian distributions of colliders at temperatures T <10000 K. A comparative analysis of our results and previous calculations in the literature is done. We discuss the effect of these rates on the formation of scattering polarization patterns of resonant lines of alkali-earths in the solar atmosphere, and their effect on Hanle effect diagnostics of solar magnetic fields.

More than about twenty central stars of planetary nebulae (CSPN) have been observed spectropolarimetrically, yet no clear, unambiguous signal of the presence of a magnetic field in these objects has been found. We perform a statistical (Bayesian) analysis of all the available spectropolarimetric observations of CSPN to constrain the magnetic fields on these objects. Assuming that the stellar field is dipolar and that the dipole axis of the objects are oriented randomly (isotropically), we find that the dipole magnetic field strength is smaller than 400 G with 95% probability using all available observations. The analysis introduced allows integration of future observations to further constrain the parameters of the distribution, and it is general, so that it can be easily applied to other classes of magnetic objects. We propose several ways to improve the upper limits found here.

Magnetic field measurements in the upper chromosphere and above, where the gas-to-magnetic pressure ratio $\beta$ is lower than unity, are essential for understanding the thermal structure and dynamical activity of the solar atmosphere. Recent developments in the theory and numerical modeling of polarization in spectral lines have suggested that information on the magnetic field of the chromosphere-corona transition region could be obtained by measuring the linear polarization of the solar disk radiation at the core of the hydrogen Lyman-$\alpha$ line at 121.6~nm, which is produced by scattering processes and the Hanle effect. The Chromospheric Lyman-$\alpha$ Spectropolarimeter (CLASP) sounding rocket experiment aims to measure the intensity (Stokes $I$) and the linear polarization profiles ($Q/I$ and $U/I$) of the hydrogen Lyman-$\alpha$ line. In this paper we clarify the information that the Hanle effect can provide by applying a Stokes inversion technique based on a database search. The database contains all theoretical $Q/I$ and $U/I$ profiles calculated in a one-dimensional semi-empirical model of the solar atmosphere for all possible values of the strength, inclination, and azimuth of the magnetic field vector, though this atmospheric region is highly inhomogeneous and dynamic. We focus on understanding the sensitivity of the inversion results to the noise and spectral resolution of the synthetic observations as well as the ambiguities and limitation inherent to the Hanle effect when only the hydrogen Lyman-$\alpha$ is used. We conclude that spectropolarimetric observations with CLASP can indeed be a suitable diagnostic tool for probing the magnetism of the transition region, especially when complemented with information on the magnetic field azimuth that can be obtained from other instruments.

Time and spatial damping of transverse magnetohydrodynamic (MHD) kink oscillations is a source of information on the cross-field variation of the plasma density in coronal waveguides. We show that a probabilistic approach to the problem of determining the density structuring from the observed damping of transverse oscillations enables us to obtain information on the two parameters that characterise the cross-field density profile. The inference is performed by computing the marginal posterior distributions for density contrast and transverse inhomo- geneity length-scale using Bayesian analysis and damping ratios for transverse oscillations under the assumption that damping is produced by resonant absorption. The obtained distributions show that, for damping times of a few oscillatory periods, low density contrasts and short inho- mogeneity length scales are more plausible in explaining observations. This means that valuable information on the cross-field density profile can be obtained even if the inversion problem, with two unknowns and one observable, is a mathematically ill-posed problem.

Context: The determination of the magnetic field vector in quiescent solar prominences is possible by interpreting the Hanle and Zeeman effects in spectral lines. However, observational measurements are scarce and lack high spatial resolution. Aims: To determine the magnetic field vector configuration along a quiescent solar prominence by interpreting spectropolarimetric measurements in the He I 1083.0 nm triplet obtained with the Tenerife Infrared Polarimeter installed at the German Vacuum Tower Telescope of the Observatorio del Teide. Methods. The He I 1083.0 nm triplet Stokes profiles are analyzed with an inversion code that takes into account the physics responsible of the polarization signals in this triplet. The results are put into a solar context with the help of extreme ultraviolet observations taken with the Solar Dynamic Observatory and the Solar Terrestrial Relations Observatory satellites. Results: For the most probable magnetic field vector configuration, the analysis depicts a mean field strength of 7 gauss. We do not find local variations in the field strength except that the field is, in average, lower in the prominence body than in the prominence feet, where the field strength reaches 25 gauss. The averaged magnetic field inclination with respect to the local vertical is 77 degrees. The acute angle of the magnetic field vector with the prominence main axis is 24 degrees for the sinistral chirality case and 58 degrees for the dextral chirality. These inferences are in rough agreement with previous results obtained from the analysis of data acquired with lower spatial resolutions.

We present high spatial resolution spectropolarimetric observations of a quiescent hedgerow prominence taken in the He I 1083.0 nm triplet. The observation consisted of a time series in sit-and-stare mode of 36 minutes of duration. The spectrograph's slit crossed the prominence body and we recorded the time evolution of individual vertical threads. Eventually, we observed the development of a dark Rayleigh-Taylor plume that propagated upward with a velocity, projected onto the plane of the sky, of 17 km/s. Interestingly, the plume apex collided with the prominence threads pushing them aside. We inferred Doppler shifts, Doppler widths, and magnetic field strength variations by interpreting the He I Stokes profiles with the HAZEL code. The Doppler shifts show that clusters of threads move coherently while individual threads have oscillatory patterns. Regarding the plume we found strong redshifts (9-12 km/s) and large Doppler widths (10 km/s) at the plume apex when it passed through the prominence body and before it disintegrated. We associate the redshifts with perspective effects while the Doppler widths are more likely due to an increase in the local temperature. No local variations of the magnetic field strength associated with the passage of the plume were found; this leads us to conclude that the plumes are no more magnetized than the surroundings. Finally, we found that some of the threads oscillations are locally damped, what allowed us to apply prominence seismology techniques to infer additional prominence physical parameters.

A multi-line Bayesian analysis of the Zeeman broadening in the solar atmosphere is presented. A hierarchical probabilistic model, based on the simple but realistic Milne-Eddington approximation to the solution of the radiative transfer equation, is used to explain the data in the optical and near infrared. Our method makes use of the full line profiles of a more than 500 spectral lines from 4000 $\AA$ to 1.8 $\mu$m. Although the problem suffers from a strong degeneracy between the magnetic broadening and any other remaining broadening mechanism, the hierarchical model allows to isolate the magnetic contribution with reliability. We obtain the cumulative distribution function for the field strength and use it to put reliable upper limits to the unresolved magnetic field strength in the solar atmosphere. The field is below 160-180 G with 90% probability.

We present nuclear spectral energy distributions (SEDs) from 1 to 18 micron of a small sample of nearby, nearly face-on and undisturbed Seyfert galaxies without prominent nuclear dust lanes. These nuclear SEDs probe the central ~35 pc of the galaxies, on average, and include photometric and spectroscopic infrared (IR) data. We use these SEDs, the clumpy torus models of Nenkova et al. and a Bayesian approach to study the sensitivity of different IR wavelengths to the torus parameters. We find that high angular resolution 8-13 micron spectroscopy alone reliably constrains the number of clumps and their optical depth (N0 and tau_V). On the other hand, we need a combination of mid- and near-IR subarcsecond resolution photometry to constrain torus width and inclination, as well as the radial distribution of the clouds (sigma, i and q). For flat radial profiles (q=0,1), it is possible to constrain the extent of the mid-IR-emitting dust within the torus (Y) when N-band spectroscopy is available, in addition to near-IR photometry. Finally, by fitting different combinations of average and individual Seyfert 1 and Seyfert 2 data, we find that, in general, for undisturbed, nearly face-on Seyferts without prominent nuclear dust lanes, the minimum combination of data necessary to reliably constrain all the torus parameters is J+K+M-band photometry + N-band spectroscopy.

While the shape of the extinction curve on the infrared is considered to be set and the extinction ratios between infrared bands are usually taken to be approximately constant, a recent number of studies point either to a spatially variable behavior on the exponent of the power law or to a different extinction law altogether. In this paper, we propose a method to analyze the overall behavior of the interstellar extinction by means of the red-clump population, and we apply it to those areas of the Milky Way where the presence of interstellar matter is heavily felt: areas located in 5 deg<l<30 deg and b=0 deg. We show that the extinction ratios traditionally used for the near infrared could be inappropriate for the inner Galaxy and we analyze the behavior of the extinction law from 1 to 8 microns.

One of the possible mechanisms responsible for the panoply of shapes in planetary nebulae is the presence of magnetic fields that drive the ejection of ionized material during the proto-planetary nebula phase. Therefore, detecting magnetic fields in such objects is of key importance for understanding their dynamics. Still, magnetic fields have not been detected using polarimetry in the central stars of planetary nebulae. Circularly polarized light spectra have been obtained with the Focal Reducer and Low Dispersion Spectrograph at the Very Large Telescope of the European Southern Observatory and the Intermediate dispersion Spectrograph and Imaging System at the William Herschel Telescope. Nineteen planetary nebulae spanning very different morphology and evolutionary stages have been selected. Most of central stars have been observed at different rotation phases to point out evidence of magnetic variability. In this paper, we present the result of two observational campaigns aimed to detect and measure the magnetic field in the central stars of planetary nebulae on the basis of low resolution spectropolarimetry. In the limit of the adopted method, we can state that large scale fields of kG order are not hosted on the central star of planetary nebulae.

We present mid-infrared (MIR) 8-13micron spectroscopy of the nuclear regions of the interacting galaxy Arp299 (IC694+NGC3690) obtained with CanariCam (CC) on the 10.4m Gran Telescopio Canarias (GTC). The high angular resolution (~0.3-0.6arcsec) of the data allows us to probe nuclear physical scales between 60 and 120pc, which is a factor of 10 improvement over previous MIR spectroscopic observations of this system. The GTC/CC spectroscopy displays evidence of deeply embedded Active Galactic Nucleus (AGN) activity in both nuclei. The GTC/CC nuclear spectrum of NGC3690/Arp299-B1 can be explained as emission from AGN-heated dust in a clumpy torus with both a high covering factor and high extinction along the line of sight. The estimated bolometric luminosity of the AGN in NGC3690 is 3.2(+/-0.6)x10^44 erg/s. The nuclear GTC/CC spectrum of IC694/Arp299-A shows 11.3micron polycyclic aromatic hydrocarbon (PAH) emission stemming from a deeply embedded (A_V~24mag) region of less than 120pc in size. There is also a continuum-emitting dust component. If associated with th putative AGN in IC694, we estimate that it would be approximately 5 times less luminous than the AGN in NGC3690. The presence of dual AGN activity makes Arp299 a good example to study such phenomenon in the early coalescence phase of interacting galaxies.

We show preliminary results of an ongoing investigation aimed at determining the configuration of the magnetic field vector in the threads of a quiescent hedgerow solar prominence using high-spatial resolution spectropolarimetric observations taken in the He I 1083.0 nm multiplet. The data consist of a two-dimensional map of a quiescent hedgerow prominence showing vertical threads. The observations were obtained with the Tenerife Infrared Polarimeter attached to the German Vacuum Tower Telescope at the Observatorio del Teide (Spain). The He I 1083.0 nm Stokes signals are interpreted with an inversion code, which takes into account the key physical processes that generate and/or modify circular and linear polarization signals in the He I 1083.0 nm triplet: the Zeeman effect, anisotropic radiation pumping, and the Hanle effect. We present initial results of the inversions, i.e, the strength and orientation of the magnetic field vector along the prominence and in prominence threads.

We describe a simple technique for the acceleration of spectro-polarimetric inversions based on principal component analysis (PCA) of Stokes profiles. This technique involves the indexing of the database models based on the sign of the projections (PCA coefficients) of the first few relevant orders of principal components of the four Stokes parameters. In this way, each model in the database can be attributed a distinctive binary number of $2^{4n}$ bits, where $n$ is the number of PCA orders used for the indexing. Each of these binary numbers (indexes) identifies a group of "compatible" models for the inversion of a given set of observed Stokes profiles sharing the same index. The complete set of the binary numbers so constructed evidently determines a partition of the database. The search of the database for the PCA inversion of spectro-polarimetric data can profit greatly from this indexing. In practical cases it becomes possible to approach the ideal acceleration factor of $2^{4n}$ as compared to the systematic search of a non-indexed database for a traditional PCA inversion. This indexing method relies on the existence of a physical meaning in the sign of the PCA coefficients of a model. For this reason, the presence of model ambiguities and of spectro-polarimetric noise in the observations limits in practice the number $n$ of relevant PCA orders that can be used for the indexing.

We present a Bayesian seismology inversion technique for propagating magnetohydrodynamic (MHD) transverse waves observed in coronal waveguides. The technique uses theoretical predictions for the spatial damping of propagating kink waves in transversely inhomogeneous coronal waveguides. It combines wave amplitude damping length scales along the waveguide with theoretical results for resonantly damped propagating kink waves to infer the plasma density variation across the oscillating structures. Provided the spatial dependence of the velocity amplitude along the propagation direction is measured and the existence of two different damping regimes is identified, the technique would enable us to fully constrain the transverse density structuring, providing estimates for the density contrast and its transverse inhomogeneity length scale.

The derivation of the effective temperature of a star is a critical first step in a detailed spectroscopic analysis. Spectroscopic methods suffer from systematic errors related to model simplifications. Photometric methods may be more robust, but are exposed to the distortions caused by interstellar reddening. Direct methods are difficult to apply, since fundamental data of high accuracy are hard to obtain. We explore a new approach in which the spectrum is used to characterize a star's effective temperature based on a calibration established by a small set of standard stars. We perform principal component analysis on homogeneous libraries of stellar spectra, then calibrate a relationship between the principal components and the effective temperature using a set of stars with reliable effective temperatures. We find that our procedure gives excellent consistency when spectra from a homogeneous set of observations are used. Systematic offsets may appear when combining observations from different sources. Using as reference the spectra of stars with high-quality spectroscopic temperatures in the Elodie library, we define a temperature scale for FG-type disk dwarfs with an internal consistency of about 50 K, in excellent agreement with temperatures from direct determinations, but distinct from widely used scales based on the infrared flux method.

The analysis of quickly damped transverse oscillations of solar coronal loops using magneto-hydrodynamic seismology allow us to infer physical parameters that are difficult to measure otherwise. Under the assumption that such damped oscillations are due to the resonant conversion of global modes into Alfven oscillations of the tube surface, we carry out a global seismological analysis of a large set of coronal loops. A Bayesian hierarchical method is used to obtain distributions for coronal loop physical parameters by means of a global analysis of a large number of observations. The resulting distributions summarise global information and constitute data-favoured information that can be used for the inversion of individual events. The results strongly suggest that internal Alfven travel times along the loop are larger than 100 s and smaller than 540 s with 95% probability. Likewise, the density contrast between the loop interior and the surrounding is larger than 2.3 and below 6.9 with 95% probability.

The detection of multiple mode harmonic kink oscillations in coronal loops enables to obtain information on coronal density stratification and magnetic field expansion using seismology inversion techniques. The inference is based on the measurement of the period ratio between the fundamental mode and the first overtone and theoretical results for the period ratio under the hypotheses of coronal density stratification and magnetic field expansion of the wave guide. We present a Bayesian analysis of multiple mode harmonic oscillations for the inversion of the density scale height and magnetic flux tube expansion, under each of the hypotheses. Then, the two models are compared using a Bayesian model comparison scheme to assess how plausible each one is, given our current state of knowledge.

Using the OSIRIS tunable narrow band imager on the 10.4m GTC (La Palma) we have mapped the SAB(rs)cd galaxy NGC 6946 over a 7.3x7.5 square arcminutes field in the emission lines of the [SII]\lambda\lambda, 6717, 6731 doublet, and in H\alpha. From these maps we have produced catalogs of the H\alpha luminosities and effective radii of 557 HII regions across the disk, and derived the [SII] emission line ratios of 370 of these. The H\alpha observations were used to derive the mean luminosity-weighted electron densities for the regions of the sample, while the [SII] line ratios allowed us to derive values of the in situ electron densities in the denser zones from which the major fraction of the radiation in these lines is emitted, for 58 of the regions. This is by far the largest data set of its kind for a single galaxy. A classical two phase model is used to derive the filling factors of the regions. We find that although the mean electron density decreases with the square root of the radius of the regions, the in situ density is essentially independent of this radius. Thus the filling factor falls systematically, as the radius and the luminosity of the regions increases, with a power law of exponent -2.23 between filling factor and radius. These measurements should enhance the perspectives for more refined physical models of HII regions.

Velocity gradients in a stellar atmospheric plasma have an impact on the anisotropy of the radiation field that illuminates each point within the medium, and this may in principle influence the scattering line polarization that results from the induced atomic level polarization. Here we analyze the emergent linear polarization profiles of the Ca II infrared triplet after solving the radiative transfer problem of scattering polarization in time-dependent hydrodynamical models of the solar chromosphere, taking into account the impact of the plasma macroscopic velocity on the atomic level polarization. We discuss the influence that the velocity and temperature shocks in the considered chromospheric models have on the temporal evolution of the scattering polarization signals of the Ca II infrared lines, as well as on the temporally averaged profiles. Our results indicate that the increase of the linear polarization amplitudes caused by macroscopic velocity gradients may be significant in realistic situations. We also study the effect of the integration time, the microturbulent velocity and the photospheric dynamical conditions, and discuss the feasibility of observing with large-aperture telescopes the temporal variation of the scattering polarization profiles. Finally, we explore the possibility of using the differential Hanle effect in the IR triplet of Ca II with the intention of avoiding the characterization of the zero-field polarization to infer magnetic fields in dynamic situations.

We present observational evidence of apparent plasma rotational motions in the feet of a solar prominence. Our study is based on spectroscopic observations taken in the He I 1083.0 nm multiplet with the Tenerife Infrared Polarimeter attached to the German Vacuum Tower Telescope. We recorded a time sequence of spectra with 34 s cadence placing the slit of the spectrograph almost parallel to the solar limb and crossing two feet of an intermediate size, quiescent hedgerow prominence. The data show opposite Doppler shifts, +/- 6 km/s, at the edges of the prominence feet. We argue that these shifts may be interpreted as prominence plasma rotating counterclockwise around the vertical axis to the solar surface as viewed from above. The evolution of the prominence seen in EUV images taken with the Solar Dynamic Observatory provided us clues to interpret the results as swirling motions. Moreover, time-distance images taken far from the central wavelength show plasma structures moving parallel to the solar limb with velocities of about 10-15 km/s. Finally, the shapes of the observed intensity profiles suggest the presence of, at least, two components at some locations at the edges of the prominence feet. One of them is typically Doppler shifted (up to 20 km/s) with respect to the other, thus suggesting the existence of supersonic counter-streaming flows along the line-of-sight.

We study the presence of reversed polarity magnetic flux in sunspot penumbra. We applied a new regularized method to deconvolve spectropolarimetric data observed with the spectropolarimeter SP onboard Hinode. The new regularization is based on a principal component decomposition of the Stokes profiles. The resulting Stokes profiles were inverted to infer the magnetic field vector using SIR. We find, for the first time, reversed polarity fields at the border of many bright penumbral filaments in the whole penumbra.

We present spectropolarimetric observations in the He I 1083.0 nm multiplet of a quiescent, hedgerow solar prominence. The data were taken with the Tenerife Infrared Polarimeter attached to the German Vacuum Tower Telescope at the Observatorio del Teide (Tenerife; Canary Islands; Spain). The observed He I circular and linear polarization signals are dominated by the Zeeman effect and by atomic level polarization and the Hanle effect, respectively. These observables are sensitive to the strength and orientation of the magnetic field vector at each spatial point of the field of view. We determine the magnetic field vector of the prominence by applying the HAZEL inversion code to the observed Stokes profiles. We briefly discuss the retrieved magnetic field vector configuration.

Extracting magnetic and thermodynamic information from spectropolarimetric observations is a difficult and time consuming task. The amount of science-ready data that will be generated by the new family of large solar telescopes is so large that we will be forced to modify the present approach to inference. In this contribution, I propose several possible ways that might be useful for extracting the thermodynamic and magnetic properties of solar plasmas from such observations quickly.

The analysis of high spectral resolution spectroscopic and spectropolarimetric observations constitute a very powerful way of inferring the dynamical, thermodynamical, and magnetic properties of distant objects. However, these techniques are photon-starving, making it difficult to use them for all purposes. One of the problems commonly found is just detecting the presence of a signal that is buried on the noise at the wavelength of some interesting spectral feature. This is specially relevant for spectropolarimetric observations because typically, only a small fraction of the received light is polarized. We present in this note a Bayesian technique for the detection of spectropolarimetric signals. The technique is based on the application of the non-parametric relevance vector machine to the observations, which allows us to compute the evidence for the presence of the signal and compute the more probable signal. The method would be suited for analyzing data from experimental instruments onboard space missions and rockets aiming at detecting spectropolarimetric signals in unexplored regions of the spectrum such as the Chromospheric Lyman-Alpha Spectro-Polarimeter (CLASP) sounding rocket experiment.

Recent success in explaining several properties of the dusty torus around the central engine of active galactic nuclei has been gathered with the assumption of clumpiness. The properties of such clumpy dusty tori can be inferred by analyzing spectral energy distributions (SEDs), sometimes with scarce sampling given that large aperture telescopes and long integration times are needed to get good spatial resolution and signal. We aim at using the information already present in the data and the assumption of clumpy dusty torus, in particular, the CLUMPY models of Nenkova et al., to evaluate the optimum next observation such that we maximize the constraining power of the new observed photometric point. To this end, we use the existing and barely applied idea of Bayesian adaptive exploration, a mixture of Bayesian inference, prediction and decision theories. The result is that the new photometric filter to use is the one that maximizes the expected utility, which we approximate with the entropy of the predictive distribution. In other words, we have to sample where there is larger variability in the SEDs compatible with the data with what we know of the model parameters. We show that Bayesian adaptive exploration can be used to suggest new observations, and ultimately optimal filter sets, to better constrain the parameters of the clumpy dusty torus models. In general, we find that the region between 10 and 200 um produces the largest increase in the expected utility, although sub-mm data from ALMA also prove to be useful. It is important to note that here we are not considering the angular resolution of the data, which is key when constraining torus parameters. Therefore, the expected utilities derived from this methodology must be weighted with the spatial resolution of the data.

We report Stokes vector observations of solar spicules and a prominence in the He I 1083 nm multiplet carried out with the Tenerife Infrared Polarimeter. The observations show linear polarization profiles that are produced by scattering processes in the presence of a magnetic field. After a careful data reduction, we demonstrate the existence of extremely asymmetric Stokes V profiles in the spicular material that we are able to model with two magnetic components along the line of sight, and under the presence of atomic orientation in the energy levels that give rise to the multiplet. We discuss some possible scenarios that can generate the atomic orientation in spicules. We stress the importance of spectropolarimetric observations across the limb to distinguish such signals from observational artifacts.

We analyze two regions of the quiet Sun (35.6 x 35.6 Mm^2) observed at high spatial resolution (~100 km) in polarized light by the IMaX spectropolarimeter onboard the Sunrise balloon. We identify 497 small-scale (~400 km) magnetic loops, appearing at an effective rate of 0.25 loop h^{-1} arcsec^{-2}; further, we argue that this number and rate are underestimated by ~30%. However, we find that these small dipoles do not appear uniformly on the solar surface: their spatial distribution is rather filamentary and clumpy, creating dead calm areas, characterized by a very low magnetic signal and a lack of organized loop-like structures at the detection level of our instruments, that cannot be explained as just statistical fluctuations of a Poisson spatial process. We argue that this is an intrinsic characteristic of the mechanism that generates the magnetic fields in the very quiet Sun. The spatio-temporal coherences and the clumpy structure of the phenomenon suggest a recurrent, intermittent mechanism for the generation of magnetic fields in the quietest areas of the Sun.

The observation of the polarization emerging from a rotating star at different phases opens up the possibility to map the magnetic field in the stellar surface thanks to the well-known Zeeman Doppler Imaging. When the magnetic field is sufficiently weak, the circular and linear polarization profiles locally in each point of the star are proportional to the first and second derivatives of the unperturbed intensity profile, respectively. We show that the weak-field approximation (for weak lines in the case of linear polarization) can be generalized to the case of a rotating star including the Doppler effect and taking into account the integration on the stellar surface. The Stokes profiles are written as a linear combination of wavelength-dependent terms expressed as series expansions in terms of Hermite polynomials. These terms contain the surface integrated magnetic field and velocity components. The direct numerical evaluation of these quantities is limited to rotation velocities not larger than 8 times the Doppler width of the local absorption profiles. Additionally, we demonstrate that, in a rotating star, the circular polarization flux depends on the derivative of the intensity flux with respect to the wavelength and also on the profile itself. Likewise, the linear polarization depends on the profile and on its first and second derivative with respect to the wavelength. We particularize the general expressions to a rotating dipole.

We report the final results of the Phase II SIMPLE measurements, comprising two run stages of 15 superheated droplet detectors each, the second stage including an improved neutron shielding. The analyses includes a refined signal analysis, and revised nucleation efficiency based on reanalysis of previously-reported monochromatic neutron irradiations. The combined results yield a contour minimum of \sigma_{p} = 4.2 x 10^-3 pb at 35 GeV/c^2 on the spin-dependent sector of WIMP-proton interactions, the most restrictive to date from a direct search experiment and overlapping for the first time results previously obtained only indirectly. In the spin-independent sector, a minimum of 3.6 x 10^-6 pb at 35 GeV/c^2 is achieved, with the exclusion contour challenging the recent CoGeNT region of current interest.

Magnetic field topology, thermal structure and plasma motions are the three main factors affecting the polarization signals used to understand our star. In this theoretical investigation, we focus on the effect that gradients in the macroscopic vertical velocity field have on the non-magnetic scattering polarization signals, establishing the basis for general cases. We demonstrate that the solar plasma velocity gradients have a significant effect on the linear polarization produced by scattering in chromospheric spectral lines. In particular, we show the impact of velocity gradients on the anisotropy of the radiation field and on the ensuing fractional alignment of the CaII levels, and how they can lead to an enhancement of the zero-field linear polarization signals. This investigation remarks the importance of knowing the dynamical state of the solar atmosphere in order to correctly interpret spectropolarimetric measurements, which is important, among other things, for establishing a suitable zero field reference case to infer magnetic fields via the Hanle effect.

Double-peaked line profiles are commonly considered a hallmark of rotating disks, with the distance between the peaks a measure of the rotation velocity. However, double-peaks can arise also from radiative transfer effects in optically thick non-rotating sources. Utilizing exact solutions of the line transfer problem we present a detailed study of line emission from geometrically thin Keplerian disks. We derive the conditions for emergence of kinematic double peaks in optically thin and thick disks, and find that it is generally impossible to disentangle the effects of kinematics and line opacity in observed double-peaked profiles. Unless supplemented by additional information, a double-peaked profile alone is not a reliable indicator of a rotating disk. In certain circumstances, triple and quadruple profiles might be better indicators of rotation in optically thick disks.

Inferring magnetic and thermodynamic information from spectropolarimetric observations relies on the assumption of a parameterized model atmosphere whose parameters are tuned by comparison with observations. Often, the choice of the underlying atmospheric model is based on subjective reasons. In other cases, complex models are chosen based on objective reasons (for instance, the necessity to explain asymmetries in the Stokes profiles) but it is not clear what degree of complexity is needed. The lack of an objective way of comparing models has, sometimes, led to opposing views of the solar magnetism because the inferred physical scenarios are essentially different. We present the first quantitative model comparison based on the computation of the Bayesian evidence ratios for spectropolarimetric observations. Our results show that there is not a single model appropriate for all profiles simultaneously. Data with moderate signal-to-noise ratios favor models without gradients along the line-of-sight. If the observations shows clear circular and linear polarization signals above the noise level, models with gradients along the line are preferred. As a general rule, observations with large signal-to-noise ratios favor more complex models. We demonstrate that the evidence ratios correlate well with simple proxies. Therefore, we propose to calculate these proxies when carrying out standard least-squares inversions to allow for model comparison in the future.

We present near-to-mid-infrared spectral energy distributions (SEDs) for 21 Seyfert galaxies, using subarcsecond resolution imaging data. Our aim is to compare the properties Seyfert 1 (Sy1) and Seyfert 2 (Sy2) tori using clumpy torus models and a Bayesian approach to fit the infrared (IR) nuclear SEDs. These dusty tori have physical sizes smaller than 6 pc radius, as derived from our fits. Active galactic nuclei (AGN) unification schemes account for a variety of observational differences in terms of viewing geometry. However, we find evidence that strong unification may not hold, and that the immediate dusty surroundings of Sy1 and Sy2 nuclei are intrinsically different. The Type 2 tori studied here are broader, have more clumps, and these clumps have lower optical depths than those of Type 1 tori. The larger the covering factor of the torus, the smaller the probability of having direct view of the AGN, and vice-versa. In our sample, Sy2 tori have larger covering factors (C_T=0.95+/-0.02) and smaller escape probabilities than those of Sy1 (C_T=0.5+/-0.1). Thus, on the basis of the results presented here, the classification of a Seyfert galaxy may depend more on the intrinsic properties of the torus rather than on its mere inclination, in contradiction with the simplest unification model.

Full Stokes filter-polarimeters are key instruments for investigating the rapid evolution of magnetic structures on the solar surface. To this end, the image quality is routinely improved using a-posteriori image reconstruction methods. We analyze the robustness of circular polarization asymmetries to phase-diversity image reconstruction techniques. We use snapshots of magneto-hydrodynamical simulations carried out with different initial conditions to synthesize spectra of the magnetically sensitive Fe I line at 5250.2 A. We degrade the synthetic profiles spatially and spectrally to simulate observations with the IMaX full Stokes filter-polarimeter. We also simulate the focused/defocused pairs of images used by the phase-diversity algorithm for reconstruction and the polarimetric modulation scheme. We assume that standard optimization methods are able to infer the projection of the wavefront on the Zernike polynomials with 10% precision. We also consider the less favorable case of 25% precision. We obtain reconstructed monochromatic modulated images that are later demodulated and compared with the original maps. Although asymmetries are often difficult to define in the quiet Sun due to the complexity of the Stokes V profiles, we show how asymmetries are degraded with spatial and spectral smearing. The results indicate that, although image reconstruction techniques reduce the spatial smearing, they can modify the asymmetries of the profiles, mainly caused by the appearance of spatially-correlated noise.

The polarised spectrum of stellar radiation encodes valuable information on the conditions of stellar atmospheres and the magnetic fields that permeate them. In this paper, we give explicit expressions to estimate the magnetic field vector and its associated error from the observed Stokes parameters. We study the solar case where specific intensities are observed and then the stellar case, where we receive the polarised flux. In this second case, we concentrate on the explicit expression for the case of a slow rotator with a dipolar magnetic field geometry. Moreover, we also give explicit formulae to retrieve the magnetic field vector from the LSD profiles without assuming mean values for the LSD artificial spectral line. The formulae have been obtained assuming that the spectral lines can be described in the weak field regime and using a maximum likelihood approach. The errors are recovered by means of the hermitian matrix. The bias of the estimators are analysed in depth.

We perform a Bayesian parameter inference in the context of resonantly damped transverse coronal loop oscillations. The forward problem is solved in terms of parametric results for kink waves in one-dimensional flux tubes in the thin tube and thin boundary approximations. For the inverse problem, we adopt a Bayesian approach to infer the most probable values of the relevant parameters, for given observed periods and damping times, and to extract their confidence levels. The posterior probability distribution functions are obtained by means of Markov Chain Monte Carlo simulations, incorporating observed uncertainties in a consistent manner. We find well localized solutions in the posterior probability distribution functions for two of the three parameters of interest, namely the Alfven travel time and the transverse inhomogeneity length-scale. The obtained estimates for the Alfven travel time are consistent with previous inversion results, but the method enables us to additionally constrain the transverse inhomogeneity length-scale and to estimate real error bars for each parameter. When observational estimates for the density contrast are used, the method enables us to fully constrain the three parameters of interest. These results can serve to improve our current estimates of unknown physical parameters in coronal loops and to test the assumed theoretical model.

We report far-infrared (FIR) imaging of the Seyfert 2 galaxy NGC 3081 in the range 70-500 micron, obtained with an unprecedented angular resolution, using the Herschel Space Observatory instruments PACS and SPIRE. The 11 kpc (~70 arcsec) diameter star-forming ring of the galaxy appears resolved up to 250 micron. We extracted infrared (1.6-500 micron) nuclear fluxes, that is active nucleus-dominated fluxes, and fitted them with clumpy torus models, which successfully reproduce the FIR emission with small torus sizes. Adding the FIR data to the near- and mid-infrared spectral energy distribution (SED) results in a torus radial extent of Ro=4(+2/-1) pc, as well as in a flat radial distribution of the clouds (i.e. the q parameter). At wavelengths beyond 200 micron, cold dust emission at T=28+/-1 K from the circumnuclear star-forming ring of 2.3 kpc (~15 arcsec) in diameter starts making a contribution to the nuclear emission. The dust in the outer parts of the galaxy is heated by the interstellar radiation field (19+/-3 K).

We used the CLUMPY torus models and a Bayesian approach to fit the infrared spectral energy distributions (SEDs) and ground-based high-angular resolution mid-infrared spectroscopy of 13 nearby Seyfert galaxies. This allowed us to put tight constraints on torus model parameters such as the viewing angle, the radial thickness of the torus Y, the angular size of the cloud distribution sigma_torus, and the average number of clouds along radial equatorial rays N_0. The viewing angle is not the only parameter controlling the classification of a galaxy into a type 1 or a type 2. In principle type 2s could be viewed at any viewing angle as long as there is one cloud along the line of sight. A more relevant quantity for clumpy media is the probability for an AGN photon to escape unabsorbed. In our sample, type 1s have relatively high escape probabilities, while in type 2s, as expected, tend to be low. Our fits also confirmed that the tori of Seyfert galaxies are compact with torus model radii in the range 1-6pc. The scaling of the models to the data also provided the AGN bolometric luminosities, which were found to be in good agreement with estimates from the literature. When we combined our sample of Seyfert galaxies with a sample of PG quasars from the literature to span a range of L_bol(AGN)~10^{43}-10^{47}erg/s, we found plausible evidence of the receding torus. That is, there is a tendency for the torus geometrical covering factor to be lower at high AGN luminosities than at low AGN luminosities. This is because at low AGN luminosities the tori appear to have wider angular sizes and more clouds along radial equatorial rays. We cannot, however rule out the possibility that this is due to contamination by extended dust structures not associated with the dusty torus at low AGN luminosities, since most of these in our sample are hosted in highly inclined galaxies. (Abridged)

We present observational evidence for oscillations of magnetic flux density in the quiet areas of the Sun. The majority of magnetic fields on the solar surface have strengths of the order of or lower than the equipartition field (300-500 G). This results in a myriad of magnetic fields whose evolution is largely determined by the turbulent plasma motions. When granules evolve they squash the magnetic field lines together or pull them apart. Here we report on the periodic deformation of the shapes of features in circular polarization observed at high resolution with Sunrise. In particular, we note that the area of patches with constant magnetic flux oscillates with time, which implies that the apparent magnetic field intensity oscillates in antiphase. The periods associated to this oscillatory pattern is compatible with the granular life-time and change abruptly, which suggests that these oscillations might not correspond to characteristic oscillatory modes of magnetic structures, but to the forcing by granular motions. In one particular case, we find three patches around the same granule oscillating in phase, which means that the spatial coherence of these oscillations can reach 1600 km. Interestingly, the same kind of oscillatory phenomenon is found also in the upper photosphere.

Obtaining the magnetic properties of non-resolved structures in the solar photosphere is always challenging and problems arise because the inversion is carried out through the numerical minimization of a merit function that depends on the proposed model. We investigate the reliability of inversions in which the stray-light contamination is obtained from the same observations as a local average. In this case, we show that it is fundamental to include the covariance between the observed Stokes profiles and the stray-light contamination. The ensuing modified merit function of the inversion process penalizes large stray-light contaminations simply because of the presence of positive correlations between the observables and the stray-light, fundamentally produced by spatially variable systematics. We caution that using the wrong merit function, artificially large stray-light contaminations might be inferred. Since this effect disappears if the stray-light contamination is obtained as an average over the full field-of-view, we recommend to take into account stray-light contamination using a global approach.

The weak-field approximation is one of the simplest models that allows us to relate the observed polarization induced by the Zeeman effect with the magnetic field vector present on the plasma of interest. It is usually applied for diagnosing magnetic fields in the solar and stellar atmospheres. A fully Bayesian approach to the inference of magnetic properties in unresolved structures is presented. The analytical expression for the marginal posterior distribution is obtained, from which we can obtain statistically relevant information about the model parameters. The role of a-priori information is discussed and a hierarchical procedure is presented that gives robust results that are almost insensitive to the precise election of the prior. The strength of the formalism is demonstrated through an application to IMaX data. Bayesian methods can optimally exploit data from filter-polarimeters given the scarcity of spectral information as compared with spectro-polarimeters. The effect of noise and how it degrades our ability to extract information from the Stokes profiles is analyzed in detail.

Using a bistatic radar echo sounding (RES) system developed for calibration of the RICE particle astrophysics experiment at the South Pole, we have studied radio frequency (RF) reflections off the bedrock. The total propagation time of ~ns-duration, vertically (z-) broadcast radio signals, as a function of polarization orientation in the horizontal plane, provides a direct probe of the geometry-dependence of the ice permittivity to a depth of 2.8 km. We observe clear birefringent asymmetries along z- in the lowest half of the ice sheet, at a fractional level ~0.3%. This result is in contrast to expectations based on measurements at Dome Fuji, for which birefringence was observed in the upper 1.5 km of the ice sheet. This effect, combined with the increased radio frequency attenuation expected near the bedrock, renders the lower half thickness of South Polar ice less favorable than the upper half of the ice sheet in terms of its ultra-high energy neutrino detection potential.

We present new mid-infrared (MIR) imaging data for three Type-1 Seyfert galaxies obtained with T-ReCS on the Gemini-South Telescope at subarcsecond resolution. Our aim is to enlarge the sample studied in a previous work to compare the properties of Type-1 and Type-2 Seyfert tori using clumpy torus models and a Bayesian approach to fit the infrared nuclear spectral energy distributions (SEDs). Thus, the sample considered here comprises 7 Type-1, 11 Type-2, and 3 intermediate-type Seyferts. The unresolved IR emission of the Seyfert 1 galaxies can be reproduced by a combination of dust heated by the central engine and direct AGN emission, while for the Seyfert 2 nuclei only dust emission is considered. These dusty tori have physical sizes smaller than 6 pc radius, as derived from our fits. Unification schemes of AGN account for a variety of observational differences in terms of viewing geometry. However, we find evidence that strong unification may not hold, and that the immediate dusty surroundings of Type-1 and Type-2 Seyfert nuclei are intrinsically different. The Type-2 tori studied here are broader, have more clumps, and these clumps have lower optical depths than those of Type-1 tori. The larger the covering factor of the torus, the smaller the probability of having direct view of the AGN, and vice-versa. In our sample, Seyfert 2 tori have larger covering factors and smaller escape probabilities than those of Seyfert 1. All the previous differences are significant according to the Kullback-Leibler divergence. Thus, on the basis of the results presented here, the classification of a Seyfert galaxy as a Type-1 or Type-2 depends more on the intrinsic properties of the torus rather than on its mere inclination towards us, in contradiction with the simplest unification model.

We report the clearest detection to date of dusty torus signatures in a Weak-Line Radio Galaxy (WLRG). The deep Spitzer InfraRed Spectrograph (IRS) rest-frame mid-infrared (MIR) spectrum of the WLRG PKS 0043-42 (z=0.116) shows a clear spectral turnover at wavelengths longer than ~20 micron suggestive of warm dust, as well as a 9.7 micron silicate absorption feature. In addition, the hard X-ray results, based on Chandra data, strongly support a picture in which PKS 0043-42 has a torus and accretion disc more typical of Strong-Line Radio Galaxies (SLRGs). The MIR and X-ray spectra are markedly different from those of other WLRGs at similar redshifts, and here we show that the former can be successfully fitted with clumpy torus models with parameters characteristic of Type-2 AGN tori: close to edge-on (i=74 deg) and relatively broad (torus angular width=60 deg), with an outer radius of 2 pc, hydrogen column density ~1.6x10^(23) cm^(-2), and AGN bolometric luminosity ~1.6x10^(44) erg s^(-1). The presence of a compact torus in PKS 0043-42 provides evidence that this WLRG is fuelled by cold, rather than hot, gas accretion. We suggest that WLRGs are a diverse population, and PKS 0043-42 may represent a type of radio galaxy in which the AGN activity has been recently re-triggered as a consequence of intermittent gas supply, or in which the covering factor of the Narrow-Line Region (NLR) clouds is relatively low.

The SIMPLE project uses superheated C2ClF5 liquid detectors to search for particle dark matter candidates. We report the results of the first stage exposure (14.1 kgd) of its latest two-stage, Phase II run, with 15 superheated droplet detectors of total active mass 0.208 kg. In combination with the results of the neutron-spin sensitive XENON10 experiment, these results yield a limit of |a_p| < 0.32, |a_n| < 0.17 for M_W = 50 GeV/c2 on the model-independent, spin-dependent sector of weakly interacting massive particle (WIMP)-nucleus interactions, and together yield a 50% reduction in the previously allowed region of the phase space. The result provides a contour minimum of {\sigma}_p ~ 2.8 x 10-2 pb at M_W = 45 GeV/c2, constituting the most restrictive direct detection limit to date against a spin-dependent WIMP-proton coupling. In the spin-independent sector, the result is seen to offer the prospect of contributing to the question of light mass WIMPs with an improvement in the current understanding of its nucleation efficiency.

Diagnostic techniques for stellar magnetic fields based upon spectropolarimetry. We propose and explore a new technique based upon the linear polarization emitted in Hanle-sensitive lines in disk-integrated stars where a dipolar magnetic field breaks the rotational symmetry of the resonance scattering polarization. A star with a simple dipolar field and a 1-0 spectral line are used to compute polarization amplitudes and angles.Predicted amplitudes are low but within reach of present instruments. A new application of the Hanle effect is proposed and analyzed, a tool that allows measuring of some of the weakest stellar magnetic fields.

We report results of a 14.1 kgd measurement with 15 superheated droplet detectors of total active mass 0.208 kg, comprising the first stage of a 30 kgd Phase II experiment. In combination with the results of the neutron-spin sensitive XENON10 experiment, these results yield a limit of |a_p| < 0.32 for M_W = 50 GeV/c2 on the spin-dependent sector of weakly interacting massive particle-nucleus interactions with a 50% reduction in the previously allowed region of the phase space formerly defined by XENON, KIMS and PICASSO. In the spin-independent sector, a limit of 2.3x10-5 pb at M_W = 45 GeV/c2 is obtained.

Accurate astrophysical polarimetry requires a proper characterization of the polarization properties of the telescope and instrumentation employed to obtain the observations. Determining the telescope and instrument Muller matrix is becoming increasingly difficult with the increase in aperture size of the new and upcoming solar telescopes. We have carried out a detailed multi-wavelength characterization of the Dunn Solar Telescope (DST) at the National Solar Observatory/Sacramento Peak as a case study and explore various possibilites for the determination of its polarimetric properties. We show that the telescope model proposed in this paper is more suitable than that in previous work in that it describes better the wavelength dependence of aluminum-coated mirrors. We explore the adequacy of the degrees of freedom allowed by the model using a novel mathematical formalism. Finally, we investigate the use of polarimeter calibration data taken at different times of the day to characterize the telescope and find that very valuable information on the telescope properties can be obtained in this manner. The results are also consistent with the entrance window polarizer measurements, opening very interesting possibilities for the calibration of future large-aperture solar telescopes such as the ATST or the EST.

We explore whether medium-resolution stellar spectra can be reconstructed from photometric observations, taking advantage of the highly compressible nature of the spectra. We formulate the spectral reconstruction as a least-squares problem with a sparsity constraint. In our test case using data from the Sloan Digital Sky Survey, only three broad-band filters are used as input. We demonstrate that reconstruction using three principal components is feasible with these filters, leading to differences with respect to the original spectrum smaller than 5%. We analyze the effect of uncertainties in the observed magnitudes and find that the available high photometric precision induces very small errors in the reconstruction. This process may facilitate the extraction of purely spectroscopic quantities, such as the overall metallicity, for hundreds of millions of stars for which only photometric information is available, using standard techniques applied to the reconstructed spectra.

The image degradation produced by atmospheric turbulence and optical aberrations is usually alleviated using post-facto image reconstruction techniques, even when observing with adaptive optics systems. These techniques rely on the development of the wavefront using Zernike functions and the non-linear optimization of a certain metric. The resulting optimization procedure is computationally heavy. Our aim is to alleviate this computationally burden. To this aim, we generalize the recently developed extended Zernike-Nijboer theory to carry out the analytical integration of the Fresnel integral and present a natural basis set for the development of the point spread function in case the wavefront is described using Zernike functions. We present a linear expansion of the point spread function in terms of analytic functions which, additionally, takes defocusing into account in a natural way. This expansion is used to develop a very fast phase-diversity reconstruction technique which is demonstrated through some applications. This suggest that the linear expansion of the point spread function can be applied to accelerate other reconstruction techniques in use presently and based on blind deconvolution.

Sunspots are the most spectacular manifestation of solar magnetism, yet, 99% of the solar surface remains 'quiet' at any time of the solar cycle. The quiet sun is not void of magnetic fields, though; they are organized at smaller spatial scales and evolve relatively fast, which makes them difficult to detect. Thus, although extensive quiet Sun magnetism would be a natural driver to a uniform, steady heating of the outer solar atmosphere, it is not clear what the physical processes involved would be due to lack of observational evidence. We report the topology and dynamics of the magnetic field in very quiet regions of the Sun from spectropolarimetric observations of the Hinode satellite, showing a continuous injection of magnetic flux with a well organized topology of Omega-loop from below the solar surface into the upper layers. At first stages, when the loop travels across the photosphere, it has a flattened (staple-like) geometry and a mean velocity ascent of $\sim3$ km/s. When the loop crosses the minimum temperature region, the magnetic fields at the footpoints become almost vertical and the loop topology ressembles a potential field. The mean ascent velocity at chromospheric height is $\sim12$ km/s. The energy input rate of these small-scale loops in the lower boundary of the chromosphere is (at least) of 1.4x10^6-2.2x10^7 erg cm-2 s-1. Our findings provide empirical evidence for solar magnetism as a multi-scale system, in which small-scale low-flux magnetism plays a crucial role, at least as important as active regions, coupling different layers of the solar atmosphere and being an important ingredient for chromospheric and coronal heating models.

The behavior of the observed polarization amplitudes with spatial resolution is a strong constraint on the nature and organization of solar magnetic fields below the resolution limit. We study the polarization of the very quiet Sun at different spatial resolutions using ground- and space-based observations. It is shown that 80% of the observed polarization signals do not change with spatial resolution, suggesting that, observationally, the very quiet Sun magnetism remains the same despite the high spatial resolution of space-based observations. Our analysis also reveals a cascade of spatial scales for the magnetic field within the resolution element. It is manifest that the Zeeman effect is sensitive to the microturbulent field usually associated to Hanle diagnostics. This demonstrates that Zeeman and Hanle studies show complementary perspectives of the same magnetism.

We explore the photospheric and chromospheric magnetic field in an emerging flux region. An image of the equivalent width of the He I 10830 A red blended component shows the presence of filamentary structures that might be interpreted as magnetic loops. We point out that the magnetic field strength in the chromosphere resembles a smoothed version of that found in the photosphere and that it is not correlated at all with the above-mentioned equivalent width map. Lacking other diagnostics, this suggests that one cannot discard the possibility that the chromospheric field we infer from the observations is tracing the lower chromosphere of the active region instead of tracing the magnetic field along loops. If the He I line is formed within magnetic loops, we point out a potential problem that appears when interpreting observations using only one component along the line-of-sight.

Determining the magnetic field related to solar spicules is vital for developing adequate models of these plasma jets, which are thought to play a key role in the thermal, dynamic and magnetic structure of the Chromosphere. Here we report on the magnetic properties of off-limb spicules in a very quiet region of the solar atmosphere, as inferred from new spectropolarimetric observations in the He I 10830 A triplet obtained with the Tenerife Infrared Polarimeter. We have used a novel inversion code for Stokes profiles caused by the joint action of atomic level polarization and the Hanle and Zeeman effects (HAZEL) to interpret the observations. Magnetic fields as strong as ~40G were detected in a very localized area of the slit, which could represent a possible lower value of the field strength of organized network spicules.

We demonstrate through numerical simulations with real data the feasibility of using compressive sensing techniques for the acquisition of spectro-polarimetric data. This allows us to combine the measurement and the compression process into one consistent framework. Signals are recovered thanks to a sparse reconstruction scheme from projections of the signal of interest onto appropriately chosen vectors, typically noise-like vectors. The compressibility properties of spectral lines are analyzed in detail. The results shown in this paper demonstrate that, thanks to the compressibility properties of spectral lines, it is feasible to reconstruct the signals using only a small fraction of the information that is measured nowadays. We investigate in depth the quality of the reconstruction as a function of the amount of data measured and the influence of noise. This change of paradigm also allows us to define new instrumental strategies and to propose modifications to existing instruments in order to take advantage of compressive sensing techniques.

We present subarcsecond resolution mid-infrared (mid-IR) photometry in the wavelength range from 8 to 20 micron of eighteen Seyfert galaxies, reporting high spatial resolution nuclear fluxes for the entire sample. We construct spectral energy distributions (SEDs) that the AGN dominates adding near-IR measurements from the literature at similar angular resolution. The IR SEDs of intermediate-type Seyferts are flatter and present higher 10 to 18 micron ratios than those of Seyfert 2. We fit the individual SEDs with clumpy torus models using the in-house-developed BayesClumpy tool. The models reproduce the high spatial resolution measurements. Regardless of the Seyfert type, even with high spatial resolution data, near- to mid-IR SED fitting poorly constrains the radial extent of the torus. For the Seyfert 2, we find that edge-on geometries are more probable than face-on views, with a number of clouds along equatorial rays of N = 5-15. The 10 micron silicate feature is generally modeled in shallow absorption. For the intermediate-type Seyferts, N and the inclination angle of the torus are lower than those of the Seyfert 2 nuclei, with the silicate feature appearing in weak emission or absent. The columns of material responsible for the X-ray absorption are larger than those inferred from the model fits for most of the galaxies, which is consistent with X-ray absorbing gas being located within the dust sublimation radius whereas the mid-IR flux arises from an area farther from the accretion disc. The fits yield both the bolometric luminosity of the intrinsic AGN and the torus integrated luminosity, from which we derive the reprocessing efficiency of the torus. In the models, the outer radial extent of the torus scales with the AGN luminosity, and we find the tori to be confined to scales less than 5 pc.

Some recent investigations of spectropolarimetric observations of the Zeeman effect in the Fe I lines at 630 nm carried out with the Hinode solar space telescope have concluded that the strength of the magnetic field vector in the internetwork regions of the quiet Sun is in the hG regime and that its inclination is predominantly horizontal. We critically reconsider the analysis of such observations and carry out a complete Bayesian analysis with the aim of extracting as much information as possible from them, including error bars. We apply the recently developed BayesME code that carries out a complete Bayesian inference for Milne-Eddington atmospheres. The sampling of the posterior distribution function is obtained with a Markov Chain Monte Carlo scheme and the marginal distributions are analyzed in detail. The Kullback-Leibler divergence is used to study the extent to which the observations introduce new information in the inference process resulting in sufficiently constrained parameters. Our analysis clearly shows that only upper limits to the magnetic field strength can be inferred with fields in the kG regime completely discarded. Furthermore, the noise level present in the analyzed Hinode observations induces a substantial loss of information for constraining the azimuth of the magnetic field. Concerning the inclination of the field, we demonstrate that some information is available to constrain it for those pixels with the largest polarimetric signal. The results also point out that the field in pixels with small polarimetric signals can be nicely reproduced in terms of a quasi-isotropic distribution.

In order to get a significant Zeeman signature in the polarised spectra of a magnetic star, we usually 'add' the contributions of numerous spectral lines; the ultimate goal is to recover the spectropolarimetric prints of the magnetic field in these line additions. Here we want to clarify the meaning of these techniques of line addition; in particular, we try to interpret the meaning of the 'pseudo-line' formed during this process and to find out why and how its Zeeman signature is still meaningful. We create a synthetic case of line addition and apply well tested standard solar methods routinely used in the research on magnetism in our nearest star. The results are convincing and the Zeeman signatures well detected; Solar methods are found to be quite efficient also for stellar observations. We statistically compare line addition with least-squares deconvolution and demonstrate that they both give very similar results as a consequence of the special statistical properties of the weights. The Zeeman signatures are unequivocally detected in this multiline approach. We may anticipate the outcome that magnetic field detection is reliable well beyond the weak-field approximation. Linear polarisation in the spectra of solar type stars can be detected when the spectral resolution is sufficiently high.

Our aim is to present a fast and general Bayesian inference framework based on the synergy between machine learning techniques and standard sampling methods and apply it to infer the physical properties of clumpy dusty torus using infrared photometric high spatial resolution observations of active galactic nuclei. We make use of the Metropolis-Hastings Markov Chain Monte Carlo algorithm for sampling the posterior distribution function. Such distribution results from combining all a-priori knowledge about the parameters of the model and the information introduced by the observations. The main difficulty resides in the fact that the model used to explain the observations is computationally demanding and the sampling is very time consuming. For this reason, we apply a set of artificial neural networks that are used to approximate and interpolate a database of models. As a consequence, models not present in the original database can be computed ensuring continuity. We focus on the application of this solution scheme to the recently developed public database of clumpy dusty torus models. The machine learning scheme used in this paper allows us to generate any model from the database using only a factor 10^-4 of the original size of the database and a factor 10^-3 in computing time. The posterior distribution obtained for each model parameter allows us to investigate how the observations constrain the parameters and which ones remain partially or completely undetermined, providing statistically relevant confidence intervals. As an example, the application to the nuclear region of Centaurus A shows that the optical depth of the clouds, the total number of clouds and the radial extent of the cloud distribution zone are well constrained using only 6 filters.

Determining the magnetic field of solar spicules is vital for developing adequate models of these plasma jets, which are thought to play a key role in the thermal, dynamic, and magnetic structure of the chromosphere. Here we report on magnetic spicule properties in a very quiet region of the off-limb solar atmosphere, as inferred from new spectropolarimetric observations in the HeI 10830 A triplet. We have used a novel inversion code for Stokes profiles caused by the joint action of atomic level polarization and the Hanle and Zeeman effects (HAZEL) to interpret the observations. Magnetic fields as strong as 40G were unambiguously detected in a very localized area of the slit, which may represent a possible lower value of the field strength of organized network spicules.

The application of Superheated Droplet Detectors (SDDs) to dark matter searches has so far been confined to the light nuclei refrigerants C2ClF5 and C4F10 (SIMPLE and PICASSO, respectively), with a principle sensitivity to spin-dependent interactions. Given the competitive results of these devices, as a result of their intrinsic insensitivity to backgrounds, we have developed a prototype trifluoroiodomethane (CF3I)-loaded SDD with increased sensitivity to spin-independent interactions as well. A low (0.102 kgd) exposure test operation of two high concentration, 1 liter devices is described, and the results compared with leading experiments in both spin-dependent and -independent sectors. Although competitive in both sectors when the difference in exposures is accounted for, a problem with fracturing of the detector gel must be addressed before significantly larger exposures can be envisioned.

Our main objective is to develop a denoising strategy to increase the signal to noise ratio of individual spectral lines of stellar spectropolarimetric observations. We use a multivariate statistics technique called Principal Component Analysis. The cross-product matrix of the observations is diagonalized to obtain the eigenvectors in which the original observations can be developed. This basis is such that the first eigenvectors contain the greatest variance. Assuming that the noise is uncorrelated a denoising is possible by reconstructing the data with a truncated basis. We propose a method to identify the number of eigenvectors for an efficient noise filtering. Numerical simulations are used to demonstrate that an important increase of the signal to noise ratio per spectral line is possible using PCA denoising techniques. It can be also applied for detection of magnetic fields in stellar atmospheres. We analyze the relation between PCA and commonly used well-known techniques like line addition and least-squares deconvolution. Moreover, PCA is very robust and easy to compute.

The polarization analysis of the light is typically carried out using modulation schemes. The light of unknown polarization state is passed through a set of known modulation optics and a detector is used to measure the total intensity passing the system. The modulation optics is modified several times and, with the aid of such several measurements, the unknown polarization state of the light can be inferred. How to find the optimal demodulation process has been investigated in the past. However, since the modulation matrix has to be measured for a given instrument and the optical elements can present problems of repeatability, some uncertainty is present in the elements of the modulation matrix and/or covariances between these elements. We analyze in detail this issue, presenting analytical formulae for calculating the covariance matrix produced by the propagation of such uncertainties on the demodulation matrix, on the inferred Stokes parameters and on the efficiency of the modulation process. We demonstrate that, even if the covariance matrix of the modulation matrix is diagonal, the covariance matrix of the demodulation matrix is, in general, non-diagonal because matrix inversion is a nonlinear operation. This propagates through the demodulation process and induces correlations on the inferred Stokes parameters.

A big challenge in solar and stellar physics in the coming years will be to decipher the magnetism of the solar outer atmosphere (chromosphere and corona) along with its dynamic coupling with the magnetic fields of the underlying photosphere. To this end, it is important to develop rigorous diagnostic tools for the physical interpretation of spectropolarimetric observations in suitably chosen spectral lines. Here we present a computer program for the synthesis and inversion of Stokes profiles caused by the joint action of atomic level polarization and the Hanle and Zeeman effects in some spectral lines of diagnostic interest, such as those of the He I 10830 A and D_3 multiplets. It is based on the quantum theory of spectral line polarization, which takes into account all the relevant physical mechanisms and ingredients (optical pumping, atomic level polarization, Zeeman, Paschen-Back and Hanle effects). The influence of radiative transfer on the emergent spectral line radiation is taken into account through a suitable slab model. The user can either calculate the emergent intensity and polarization for any given magnetic field vector or infer the dynamical and magnetic properties from the observed Stokes profiles via an efficient inversion algorithm based on global optimization methods. The reliability of the forward modeling and inversion code presented here is demonstrated through several applications, which range from the inference of the magnetic field vector in solar active regions to determining whether or not it is canopy-like in quiet chromospheric regions. This user-friendly diagnostic tool called "HAZEL" (from HAnle and ZEeman Light) is offered to the astrophysical community, with the hope that it will facilitate new advances in solar and stellar physics.

The analysis of near infrared spectropolarimetric data at the internetwork at different regions on the solar surface could offer constraints to reject current modeling of these quiet areas. We present spectro-polarimetric observations of very quiet regions for different values of the heliocentric angle for the Fe I lines at 1.56 micron, from disc centre to positions close to the limb. The spatial resolution of the data is 0.7-1". We analyze direct observable properties of the Stokes profiles as the amplitude of circular and linear polarization as well as the total degree of polarization. Also the area and amplitude asymmetries are studied. We do not find any significant variation of the properties of the polarimetric signals with the heliocentric angle. This means that the magnetism of the solar internetwork remains the same regardless of the position on the solar disc. This observational fact discards the possibility of modeling the internetwork as a Network-like scenario. The magnetic elements of internetwork areas seem to be isotropically distributed when observed at our spatial resolution.

The reliability of quiet Sun magnetic field diagnostics based on the \ion{Fe}{1} lines at 6302 \AA has been questioned by recent work. We present here the results of a thorough study of high-resolution multi-line observations taken with the new spectro-polarimeter SPINOR, comprising the 5250 and 6302 \AA spectral domains. The observations were analyzed using several inversion algorithms, including Milne-Eddington, LTE with 1 and 2 components, and MISMA codes. We find that the line-ratio technique applied to the 5250 \AA lines is not sufficiently reliable to provide a direct magnetic diagnostic in the presence of thermal fluctuations and variable line broadening. In general, one needs to resort to inversion algorithms, ideally with realistic magneto-hydrodynamical constrains. When this is done, the 5250 \AA lines do not seem to provide any significant advantage over those at 6302 \AA . In fact, our results point towards a better performance with the latter (in the presence of turbulent line broadening). In any case, for very weak flux concentrations, neither spectral region alone provides sufficient constraints to fully disentangle the intrinsic field strengths. Instead, we advocate for a combined analysis of both spectral ranges, which yields a better determination of the quiet Sun magnetic properties. Finally, we propose the use of two other \ion{Fe}{1} lines (at 4122 and 9000 \AA) with identical line opacities that seem to work much better than the others.

[abridged] Inversion techniques are the most powerful methods to obtain information about the thermodynamical and magnetic properties of solar and stellar atmospheres. In the last years, we have witnessed the development of highly sophisticated inversion codes that are now widely applied to spectro-polarimetric observations. The majority of these inversion codes are based on the optimization of a complicated non-linear merit function. However, no reliable and statistically well-defined confidence intervals can be obtained for the parameters inferred from the inversions. A correct estimation of the confidence intervals for all the parameters that describe the model is mandatory. Additionally, it is fundamental to apply efficient techniques to assess the ability of models to reproduce the observations and to what extent the models have to be refined or can be simplified. Bayesian techniques are applied to analyze the performance of the model to fit a given observed Stokes vector. The posterior distribution, is efficiently sampled using a Markov Chain Monte Carlo method. For simplicity, we focus on the Milne-Eddington approximate solution of the radiative transfer equation and we only take into account the generation of polarization through the Zeeman effect. However, the method is extremely general and other more complex forward models can be applied. We illustrate the ability of the method with the aid of academic and realistic examples. We show that the information provided by the posterior distribution turns out to be fundamental to understand and determine the amount of information available in the Stokes profiles in these particular cases.

We investigate the statistical properties of the extreme events of the solar cycle as measured by the sunspot number. The recent advances in the methodology of the theory of extreme values is applied to the maximal extremes of the time series of sunspots. We focus on the extreme events that exceed a carefully chosen threshold and a generalized Pareto distribution is fitted to the tail of the empirical cumulative distribution. A maximum likelihood method is used to estimate the parameters of the generalized Pareto distribution and confidence levels are also given to the parameters. Due to the lack of an automatic procedure for selecting the threshold, we analyze the sensitivity of the fitted generalized Pareto distribution to the exact value of the threshold. According to the available data, that only spans the previous ~250 years, the cumulative distribution of the time series is bounded, yielding an upper limit of 324 for the sunspot number. We also estimate that the return value for each solar cycle is ~188, while the return value for a century increases to ~228. Finally, the results also indicate that the most probable return time for a large event like the maximum at solar cycle 19 happens once every ~700 years and that the probability of finding such a large event with a frequency smaller than ~50 years is very small. In spite of the essentially extrapolative character of these results, their statistical significance is very large.

%Context: {Previous studies have indicated that the 372.4 GHz ground transition of ortho-H$_2$D$^+$ might be a powerful probe of Proto-Planetary Disks. The line could be especially suited for study of the disk mid-plane, where the bulk of the mass resides and where planet formation takes place.} %Aims: {Provide detailed theoretical predictions for the line intensity, profile and maps expected for representative disk models.} %Methods: {We determine the physical and chemical structure of the disks from the model developed by Ceccarelli & Dominik (2005). The line emission is computed with the new radiative transfer method developed recently by Elitzur & Asensio Ramos (2006).} %Results: {We present intensity maps convolved with the expected ALMA resolution, which delineate the origin of the H$_2$D$^+$ 372.4 GHz line. In the disk inner regions, the line probes the conditions in the mid-plane out to radial distances of a few tens of AU, where Solar-like planetary systems might form. In the disk outermost regions, the line originates from slightly above the mid-plane. When the disk is spatially resolved, the variation of line profile across the image provides important information about the velocity field. Spectral profiles of the entire disk flux show a double peak shape at most inclination angles.} %Conclusions: {Our study confirms that the 372.4 GHz H$_2$D$^+$ line provides powerful diagnostics of the mid-plane of Proto-Planetary Disks. Current submillimeter telescopes are capable of observing this line, though with some difficulties. The future ALMA interferometer will have the sensitivity to observe and even spatially resolve the H$_2$D$^+$ line emission.}

The physical interpretation of spectropolarimetric observations of lines of neutral helium, such as those of the 10830 A multiplet, represents an excellent opportunity for investigating the magnetism of plasma structures in the solar chromosphere. Here we present a powerful forward modeling and inversion code that permits either to calculate the emergent intensity and polarization for any given magnetic field vector or to infer the dynamical and magnetic properties from the observed Stokes profiles. This diagnostic tool is based on the quantum theory of spectral line polarization, which self-consistently accounts for the Hanle and Zeeman effects in the most general case of the incomplete Paschen-Back effect regime. We also take into account radiative transfer effects. An efficient numerical scheme based on global optimization methods has been applied. Our Stokes inversion code permits a fast and reliable determination of the global minimum.

The amount of information available in spectro-polarimetric data is estimated. To this end, the intrinsic dimensionality of the data is inferred with the aid of a recently derived estimator based on nearest-neighbor considerations and obtained applying the principle of maximum likelihood. We show in detail that the estimator correctly captures the intrinsic dimension of artificial datasets with known dimension. The effect of noise in the estimated dimension is analyzed thoroughly and we conclude that it introduces a positive bias that needs to be accounted for. Real simultaneous spectro-polarimetric observations in the visible 630 nm and the near-infrared 1.5 microns spectral regions are also investigated in detail, showing that the near-infrared dataset provides more information of the physical conditions in the solar atmosphere than the visible dataset. Finally, we demonstrate that the amount of information present in an observed dataset is a monotonically increasing function of the number of available spectral lines.

This paper reviews some of the most recent advances in the application of the Hanle effect to solar physics, and how these developments are allowing us to explore the magnetism of the photospheric regions that look ``empty'' in solar magnetograms--that is, the Sun's ``hidden'' magnetism. In particular, we show how a joint analysis of the Hanle effect in atomic and molecular lines indicates that there is a vast amount of hidden magnetic energy and unsigned magnetic flux localized in the (intergranular) downflowing regions of the quiet solar photosphere, carried mainly by tangled fields at sub-resolution scales with strengths between the equipartition field values and 1 kG.

% context heading (optional) {With the development of new polarimeters for large telescopes, the spectro-polarimetric study of astrophysical bodies is becoming feasible and, indeed, more frequent. In particular, this is permitting the observational study of stellar magnetic fields} % aims heading (mandatory) {With the aim to optimize and interpret this kind of observations, we have produced a spectral atlas of circular polarization in a grid of stellar atmospheric models with effective temperatures between 3500 and 10000 K, surface gravities $\log(g)=3.5-5$, metallicities between 10$^{-2}$ and 1, and magnetic field strengths of 100, 1000 and 5000~G} % methods heading (mandatory) {We have computed the emergent Stokes $I$ and $V$ flux spectra in LTE of more than 10$^5$ spectral lines} % results heading (mandatory) {The atlas and several numerical tools are available in electronic format and may be downloaded from http://download.hao.ucar.edu/pub/PSA/. In this paper we review and discuss some of its most relevant features, such as which spectral regions and individual lines harbor the strongest signals, what are interesting lines to observe, how to disentangle field strength from filling factor, etc.}

We report on spectropolarimetric observations of a near-IR line of Mn I located at 15262.702 A whose intensity and polarization profiles are very sensitive to the presence of hyperfine structure. A theoretical investigation of the magnetic sensitivity of this line to the magnetic field uncovers several interesting properties. The most important one is that the presence of strong Paschen-Back perturbations due to the hyperfine structure produces an intensity line profile whose shape changes according to the absolute value of the magnetic field strength. A line ratio technique is developed from the intrinsic variations of the line profile. This line ratio technique is applied to spectropolarimetric observations of the quiet solar photosphere in order to explore the probability distribution function of the magnetic field strength. Particular attention is given to the quietest area of the observed field of view, which was encircled by an enhanced network region. A detailed theoretical investigation shows that the inferred distribution yields information on the average magnetic field strength and the spatial scale at which the magnetic field is organized. A first estimation gives ~250 G for the mean field strength and a tentative value of ~0.45" for the spatial scale at which the observed magnetic field is horizontally organized.

The polarization observed in the spectral lines of the He I 10830 multiplet carries valuable information on the dynamical and magnetic properties of plasma structures in the solar chromosphere and corona, such as spicules, prominences, filaments, emerging magnetic flux regions, etc. Here we investigate the influence of atomic level polarization on the emergent Stokes profiles for a broad range of magnetic field strengths, in both 90 degree and forward scattering geometry. We show that, contrary to a widespread belief, the selective emission and absorption processes caused by the presence of atomic level polarization may have an important influence on the emergent linear polarization, even for magnetic field strengths as large as 1000 G. Consequently, the modeling of the Stokes Q and U profiles should not be done by taking only into account the contribution of the transverse Zeeman effect within the framework of the Paschen-Back effect theory, unless the magnetic field intensity of the observed plasma structure is sensibly larger than 1000 G. We point out also that in low-lying optically thick plasma structures, such as those of active region filaments, the (horizontal) radiation field generated by the structure itself may substantially reduce the positive contribution to the anisotropy factor caused by the (vertical) radiation field coming from the underlying solar photosphere, so that the amount of atomic level polarization may turn out to be negligible. Only under such circumstances may the emergent linear polarization of the He I 10830 multiplet in such regions of the solar atmosphere be dominated by the contribution caused by the transverse Zeeman effect.

It is shown that the two-part Minimum Description Length Principle can be used to discriminate among different models that can explain a given observed dataset. The description length is chosen to be the sum of the lengths of the message needed to encode the model plus the message needed to encode the data when the model is applied to the dataset. It is verified that the proposed principle can efficiently distinguish the model that correctly fits the observations while avoiding over-fitting. The capabilities of this criterion are shown in two simple problems for the analysis of observed spectropolarimetric signals. The first is the de-noising of observations with the aid of the PCA technique. The second is the selection of the optimal number of parameters in LTE inversions. We propose this criterion as a quantitative approach for distinguising the most plausible model among a set of proposed models. This quantity is very easy to implement as an additional output on the existing inversion codes.

We have conducted a study combining H2O lines in two spectral regions. First, Infrared Space Observatory observations of several H2O thermal lines seen in absorption toward Sgr B2(M) at a spectral resolution of 35 kms^-1 have been analyzed. Second, an IRAM-30m telescope map of the para-H2O 3_13-2_20 line at 183.31 GHz, seen in emission, has also been obtained and analyzed. The H2O lines seen in absorption are optically thick and are formed in the outermost gas of the condensations in front of the far-IR continuum sources. They probe a maximum visual extinction of ~5 to 10 mag. Radiative transfer models indicate that these lines are quite insensitive to temperature and gas density, and that IR photons from the dust play a dominant role in the excitation of the involved H2O rotational levels. The water vapor abundance in the region is (1-2)x10^-5. The relatively low H2O/OH abundance ratio in the region, 2-4, is a signature of UV photon dominated surface layers traced by far-IR observations.

The development of fast numerical methods for multilevel radiative transfer (RT) applications often leads to important breakthroughs in astrophysics, because they allow the investigation of problems that could not be properly tackled using the methods previously available. Probably, the most familiar example is the so-called Multilevel Accelerated $\Lambda$-Iteration (MALI) technique of Rybicki & Hummer for the case of a local approximate operator, which is based on Jacobi iteration. However, there are superior operator-splitting methods, based on Gauss-Seidel (GS) and Successive Overrelaxation (SOR) iteration, which provide a dramatic increase in the speed with which non-LTE multilevel transfer problems can be solved in one, two and three-dimensional geometries. Such RT methods, which were introduced by Trujillo Bueno & Fabiani Bendicho ten years ago, are the main subject of the first part of this paper. We show in some detail how they can be applied for solving multilevel RT problems in spherical geometry, for both atomic and molecular line transitions. The second part of the article addresses the issue of the calculation of the molecular number densities when the approximation of instantaneous chemical equilibrium turns out to be inadequate, which happens to be the case whenever the dynamical time scales of the astrophysical plasma under consideration are much shorter than the time needed by the molecules to form.

Analysis of the Hanle effect in solar molecular lines allows us to obtain empirical information on hidden, mixed-polarity magnetic fields at subresolution scales in the (granular) upflowing regions of the `quiet' solar photosphere. Here we report that collisions seem to be very efficient in depolarizing the rotational levels of MgH lines. This has the interesting consequence that in the upflowing regions of the quiet solar photosphere the strength of the hidden magnetic field cannot be sensibly larger than 10 G, assuming the simplest case of a single valued microturbulent field that fills the entire upflowing photospheric volume. Alternatively, an equally good theoretical fit to the observed scattering polarization amplitudes can be achieved by assuming that the rate of depolarizing collisions is an order of magnitude smaller than in the previous collisionally dominated case, but then the required strength of the hidden field in the upflowing regions turns out to be unrealistically high. These constraints reinforce our previously obtained conclusion that there is a vast amount of hidden magnetic energy and unsigned magnetic flux localized in the (intergranular) downflowing regions of the quiet solar photosphere.

We present a new method, the Coupled Escape Probability (CEP), for exact calculation of line emission from multi-level systems, solving only algebraic equations for the level populations. The CEP formulation of the classical two-level problem is a set of {\em linear equations}, and we uncover an exact analytic expression for the emission from two-level optically thick sources that holds as long as they are in the "effectively thin" regime. In comparative study of a number of standard problems, the CEP method outperformed the leading line transfer methods by substantial margins. The algebraic equations employed by our new method are already incorporated in numerous codes based on the escape probability approximation. All that is required for an exact solution with these existing codes is to augment the expression for the escape probability with simple zone-coupling terms. As an application, we find that standard escape probability calculations generally produce the correct cooling emission by the CII 158 mic line but not by the $^3$P lines of OI.

This paper describes a very general approach to the calculation of the Zeeman splitting effect produced by an external magnetic field on the rotational levels of diatomic molecules. The method is valid for arbitrary values of the total electronic spin and of the magnetic field strength -that is, it holds for molecular electronic states of any multiplicity and for both the Zeeman and incomplete Paschen-Back regimes. It is based on an efficient numerical diagonalization of the effective Zeeman Hamiltonian, which can incorporate easily all the contributions one may eventually be interested in, such as the hyperfine interaction of the external magnetic field with the spin motions of the nuclei. The reliability of the method is demonstrated by comparing our results with previous ones obtained via formulae valid only for doublet states. We also present results for molecular transitions arising between non-doublet electronic states, illustrating that their Zeeman patterns show signatures produced by the Paschen-Back effect.

A novel approach is presented for the solution of instantaneous chemical equilibrium problems. The chemical equilibrium can be considered, due to its intrinsically local character, as a mapping of the three-dimensional parameter space spanned by the temperature, hydrogen density and electron density into many one-dimensional spaces representing the number density of each species. We take advantage of the ability of artificial neural networks to approximate non-linear functions and construct neural networks for the fast and efficient solution of the chemical equilibrium problem in typical stellar atmosphere physical conditions. The neural network approach has the advantage of providing an analytic function, which can be rapidly evaluated. The networks are trained with a learning set (that covers the entire parameter space) until a relative error below 1% is reached. It has been verified that the networks are not overtrained by using an additional verification set. The networks are then applied to a snapshot of realistic three-dimensional convection simulations of the solar atmosphere showing good generalization properties.

We report novel spectropolarimetric observations of sunspots carried out with the Tenerife Infrared Polarimeter (TIP) in a near-IR spectral region around 15410 A, which is known to contain two groups of prominent OH lines that show circular polarization signals of opposite polarity. Surrounding these well-known OH lines, we have discovered the presence of CN lines of the $\Delta v=1$ band which show anomalous polarization profiles. Although the Stokes V signals of the OH lines are antisymmetric and with a sizable amplitude, the CN lines show almost negligible circular polarization. On the contrary, the linear polarization signals turn out to be much stronger in the CN lines than in the OH lines. Interestingly, these CN lines present striking antisymmetric linear polarization profiles, which we are able to explain and model via the Paschen-Back effect theory for diatomic molecules. The presence of such peculiar CN lines in the same spectral region of the OH lines may be useful to improve our empirical knowledge of solar magnetic fields via the simultaneous observation and modeling of the transverse and longitudinal Zeeman effects in two different molecular species.

The theory of the generation and transfer of polarized radiation, mainly developed for interpreting solar spectropolarimetric observations, allows to reconsider, in a more rigorous and elegant way, a physical mechanism that has been suggested some years ago to interpret the high degree of polarization often observed in astronomical masers. This mechanism, for which the name of 'dichroic maser' is proposed, can operate when a low density molecular cloud is illuminated by an anisotropic source of radiation (like for instance a nearby star). Here we investigate completely unsaturated masers and show that selective stimulated emission processes are capable of producing highly polarized maser radiation in a non-magnetic environment. The polarization of the maser radiation is linear and is directed tangentially to a ring equidistant to the central star. We show that the Hanle effect due to the presence of a magnetic field can produce a rotation (from the tangential direction) of the polarization by more that 45 degrees for some selected combinations of the strength, inclination and azimuth of the magnetic field vector. However, these very same conditions produce a drastic inhibition of the maser effect. The rotations of about 90 degrees observed in SiO masers in the evolved stars TX Cam by Kemball & Diamond (1997) and IRC+10011 by Desmurs et al (2000) may then be explainedby a local modification of the anisotropy of the radiation field, being transformed from mainly radial to mainly tangential.

Deciphering and understanding the small-scale magnetic activity of the quiet solar photosphere should help to solve many of the key problems of solar and stellar physics, such as the magnetic coupling to the outer atmosphere and the coronal heating. At present, we can see only ${\sim}1%$ of the complex magnetism of the quiet Sun, which highlights the need to develop a reliable way to investigate the remaining 99%. Here we report three-dimensional radiative tranfer modelling of scattering polarization in atomic and molecular lines that indicates the presence of hidden, mixed-polarity fields on subresolution scales. Combining this modelling with recent observational data we find a ubiquitous tangled magnetic field with an average strength of ${\sim}130$ G, which is much stronger in the intergranular regions of solar surface convection than in the granular regions. So the average magnetic energy density in the quiet solar photosphere is at least two orders of magnitude greater than that derived from simplistic one-dimensional investigations, and sufficient to balance radiative energy losses from the solar chromosphere.

Here we report on the first observational investigation of the Zeeman effect in the G band around 4305 A. Our spectropolarimetric observations of sunspots with the Zurich Imaging Polarimeter at the IRSOL facility confirm our previous theoretical prediction that the molecular Zeeman effect produces measurable circular polarization signatures in several CH lines that are not overlapped with atomic transitions. We also find both circular and linear polarization signals produced by atomic lines whose wavelengths lie in the G band spectral region. Together, such molecular and atomic lines are potentially important for empirical investigations on solar and stellar magnetism. For instance, a comparison between observed and calculated Stokes profiles suggests that the thermodynamical and/or magnetic properties of the photospheric regions of sunspot umbrae are horizontally structured with a component that might be associated with umbral dots.

We fit the surface-brightness profiles of 21 elliptical galaxies using both the Sersic function and a new empirical model which combines an inner power law with an outer Sersic function. The profiles are combinations of deconvolved HST profiles from the literature and ellipse fits to the full WFPC2 mosaic images, and thus span a radial range from ~ 0.02 arcsec to ~ twice the half-light radius. We are able to accurately fit the entire profiles using either the Sersic function or our new model. In doing so, we demonstrate that most, if not all, so-called ``power-law'' galaxies are better described as ``Sersic galaxies'' -- they are well modeled by the three-parameter Sersic profile into the limits of HST resolution -- and that ``core'' galaxies are best understood as consisting of an outer Sersic profile with an inner power-law cusp, which is a downward deviation from the inward extrapolation of the Sersic profile. This definition of cores resolves ambiguities that result when the popular ``Nuker law'' is fitted to the profiles of ellipticals and bulges, particularly at lower luminosities. We also find that using the Nuker law to model core-galaxy nuclear profiles systematically overestimates the core radii by factors of 1.5--4.5 and underestimates the inner power-law slope by ~ 20--40% or more.

Here we report the first detection of polarization signals induced by the Zeeman effect in spectral lines of the E^4\Pi-A^4\Pi system of FeH located around 1.6 $\mu$m. Motivated by the tentative detection of this band in the intensity spectrum of late-type dwarfs, we have investigated the full Stokes sunspot spectrum finding circular and linear polarization signatures that we associate with the FeH lines of the E^4\Pi-A^4\Pi band system. We investigate the Zeeman effect in these molecular transitions pointing out that in Hund's case (a) coupling the effective Land\'e factors are never negative. For this reason, the fact that our spectropolarimetric observations indicate that the Land\'e factors of pairs of FeH lines have opposite signs, prompt us to conclude that the E^4\Pi-A^4\Pi system must be in intermediate angular momentum coupling between Hund's cases (a) and (b). We emphasize that theoretical and/or laboratory investigations of this molecular system are urgently needed for exploiting its promising diagnostic capabilities.

We investigate the possibility of measuring magnetic field strength in G-band bright points through the analysis of Zeeman polarization in molecular CH lines. To this end we solve the equations of polarized radiative transfer in the G band through a standard plane-parallel model of the solar atmosphere with an imposed magnetic field, and through a more realistic snapshot from a simulation of solar magneto-convection. This region of the spectrum is crowded with many atomic and molecular lines. Nevertheless, we find several instances of isolated groups of CH lines that are predicted to produce a measurable Stokes V signal in the presence of magnetic fields. In part this is possible because the effective Land\'{e} factors of lines in the stronger main branch of the CH A$^{2}\Delta$--X$^{2}\Pi$ transition tend to zero rather quickly for increasing total angular momentum $J$, resulting in a Stokes $V$ spectrum of the G band that is less crowded than the corresponding Stokes $I$ spectrum. We indicate that, by contrast, the effective Land\'{e} factors of the $R$ and $P$ satellite sub-branches of this transition tend to $\pm 1$ for increasing $J$. However, these lines are in general considerably weaker, and do not contribute significantly to the polarization signal. In one wavelength location near 430.4 nm the overlap of several magnetically sensitive and non-sensitive CH lines is predicted to result in a single-lobed Stokes $V$ profile, raising the possibility of high spatial-resolution narrow-band polarimetric imaging. In the magneto-convection snapshot we find circular polarization signals of the order of 1% prompting us to conclude that measuring magnetic field strength in small-scale elements through the Zeeman effect in CH lines is a realistic prospect.

Employing the most recent parametrization of the baryon-baryon interaction of the Nijmegen group, we investigate, in the framework of the Brueckner--Bethe--Goldstone many-body theory at zero temperature, the influence of neutrino trapping on the composition, equation of state, and structure of neutron stars, relevant to describe the physical conditions of a neutron star immediately after birth (protoneutron star). We find that the presence of neutrinos changes significantly the composition of matter delaying the appearance of hyperons and making the equation of state stiffer. We explore the consequences of neutrino trapping on the early evolution of a neutron star and on the nature of the final compact remnant left by the supernova explosion.

The Nuker law was designed to match the inner few (~3-10) arcseconds of predominantly nearby (< 30 Mpc) early-type galaxy light-profiles; it was never intended to describe an entire profile. The Sersic model, on the other hand, was developed to fit the entire profile; however, due to the presence of partially depleted galaxy cores, the Sersic model cannot always describe the very inner region. We have therefore developed a new empirical model consisting of an inner power-law, a transition region, and an outer Sersic model to connect the inner and outer structure of elliptical galaxies. Moreover, the stability of the Nuker model parameters are investigated. Surprisingly, none are found to be stable quantities; all are shown to vary systematically with a profile's fitted radial extent, and often by more than 100%. Considering elliptical galaxies spanning a range of 7.5 magnitudes, the central stellar density of the underlying host galaxy is observed to increase with galaxy luminosity until the onset of core formation, detected only in the brightest elliptical galaxies. We suggest that the so-called ``power-law'' galaxies may actually be described by the Sersic model over their entire radial range.

Investigating the reliability of the assumption of instantaneous chemical equilibrium (ICE) for calculating the CO number density in the solar atmosphere is of crucial importance for the resolution of the long-standing controversy over the existence of `cool clouds' in the chromosphere, and for determining whether the cool gas owes its existence to CO radiative cooling or to a hydrodynamical process. Here we report the first results of such an investigation in which we have carried out time-dependent gas-phase chemistry calculations in radiation hydrodynamical simulations of solar chromospheric dynamics. We show that while the ICE approximation turns out to be suitable for modeling the observed infrared CO lines at the solar disk center, it may substantially overestimate the `heights of formation' of strong CO lines synthesized close to the edge of the solar disk, especially concerning vigorous dynamic cases resulting from relatively strong photospheric disturbances. This happens because during the cool phases of the hydrodynamical simulations the CO number density in the outer atmospheric regions is smaller than what is stipulated by the ICE approximation, resulting in decreased CO opacity in the solar chromosphere. As a result, the cool CO-bearing gas which produces the observed molecular lines must be located at atmospheric heights not greater than 700 km, approximately. We conclude that taking into account the non-equilibrium chemistry improves the agreement with the available on-disk and off-limb observations, but that the hydrodynamical simulation model has to be even cooler than anticipated by the ICE approximation, and this has to be the case at the `new' (i.e. deeper) formation regions of the rovibrational CO lines.

We have investigated the structural and dynamical properties of triaxial stellar systems whose surface brightness profiles follow the r^{1/n} luminosity law - extending the analysis of Ciotti (1991) who explored the properties of spherical r^{1/n} systems. A new `analytical' expression that accurately reproduces the spatial (i.e. deprojected) luminosity density profiles (error < 0.1%) is presented for detailed modelling of the Sersic family of luminosity profiles. We evaluate both the symmetric and the non--axisymmetric components of the gravitational potential and force, and compute the torques as a function of position. For a given triaxiality, stellar systems with smaller values of n have a greater non-axisymmetric gravitational field component. We also explore the strength of the non-axisymmetric forces produced by bulges with differing n and triaxiality on systems having a range of bulge-to-disc ratios. The increasing disc-to-bulge ratio with increasing galaxy type (decreasing n) is found to heavily reduce the amplitude of the non-axisymmetric terms, and therefore reduce the possibility that triaxial bulges in late-type systems may be the mechanism or perturbation for nonsymmetric structures in the disc. Using seeing-convolved r^{1/n}-bulge plus exponential-disc fits to the K-band data from a sample of 80 nearby disc galaxies, we probe the relations between galaxy type, Sersic index n and the bulge-to-disc luminosity ratio. These relations are shown to be primarily a consequence of the relation between n and the total bulge luminosity. In the K-band, the trend of decreasing bulge-to-disc luminosity ratio along the spiral Hubble sequence is predominantly, although not entirely, a consequence of the change in the total bulge luminosity; the trend between the total disc luminosity and Hubble type is much weaker.

Narrow band imaging in the G-band is commonly used to trace the small magnetic field concentrations of the Sun, although the mechanism that makes them bright has remained unclear. We carry out LTE syntheses of the G-band in an assorted set of semi-empirical model magnetic concentrations. The syntheses include all CH lines as well as the main atomic lines within the band-pass. The model atmospheres produce bright G-band spectra having many properties in common with the observed G-band bright points. In particular, the contrast referred to the quiet Sun is about twice the contrast in continuum wavelengths. The agreement with observations does not depend on the specificities of the model atmosphere, rather it holds from single fluxtubes to MIcro-Structured Magnetic Atmospheres. However, the agreement requires that the real G-band bright points are not spatially resolved, even in the best observations. Since the predicted G-band intensities exceed by far the observed values, we foresee a notable increase of contrast of the G-band images upon improvement of the angular resolution. According to the LTE modeling, the G-band spectrum emerges from the deep photosphere that produces the continuum. Our syntheses also predict solar magnetic concentrations showing up in continuum images but not in the G-band . Finally, we have examined the importance of the CH photo-dissociation in setting the amount of G-band absorption. It turns out to play a minor role.

The highly convergent iterative methods developed by Trujillo Bueno and Fabiani Bendicho (1995) for radiative transfer (RT) applications are generalized to spherical symmetry with velocity fields. These RT methods are based on Jacobi, Gauss-Seidel (GS), and SOR iteration and they form the basis of a new NLTE multilevel transfer code for atomic and molecular lines. The benchmark tests carried out so far are presented and discussed. The main aim is to develop a number of powerful RT tools for the theoretical interpretation of molecular spectra.

We discuss the kaon-nucleon interaction and its consequences for the change of the properties of the kaon in the medium. The onset of kaon condensation in neutron stars under various scenarios as well its effects for neutron star properties are reviewed.