We evaluate an adaptive gaussian quadrature integration scheme that will be suitable for the numerical evaluation of generalized redistribution in frequency functions. The latter are indispensable ingredients for "full non-LTE" radiation transfer computations i.e., assuming potential deviations of the velocity distribution of massive particles from the usual Maxwell-Boltzmann distribution. A first validation is made with computations of the usual Voigt profile.

Empirical correlations between stellar parameters such as rotation or radius and magnetic activity diagnostics require estimates of the effective temperatures and the stellar radii. The aim of this study is to propose simple methods that can be applied to large samples of stars in order to derive estimates of the stellar parameters. Good empirical correlations between Red/Infra-Red colours (e.g. (R-I)$_C$) and effective temperatures have been well established for a long time (e.g. Veeder 1974, Bessell 1979, Leggett 1992). The more recent (R-I)$_C$ colour-$T_{eff}$ correlation using the data of Mann et al. (2015: hereafter M15) and Boyajian et al. (2012: hereafter B12) shows that this colour can be applied as a temperature estimate for large samples of stars. We find that the mean scatter in $T_{eff}$ relative to the (R-I)$_C$-$T_{eff}$ relationship of B12 and M15 data is only $\pm 3\sigma =$44.6~K for K dwarfs and $\pm 3\sigma =$39.4~K for M dwarfs. These figures are small and show that the (R-I)$_C$ colour can be used as a first guess effective temperature estimator for K and M dwarfs. We derive effective temperatures for about 1910 K and M dwarfs using the calibration of (R-I)$_C$ colour-$T_{eff}$ from B12 and M15 data. We also compiled $T_{eff}$ and metallicity measurements available in the literature using the VizieR database. We determine $T_{eff}$ for 441 stars with previously unknown effective temperatures. We also identified 21 new spectroscopic binaries and 1 triple system from our high resolution spectra.

A tutorial associated with the distribution of 2D non-LTE multilevel atom radiative transfer codes.

We present a method for deriving stellar fundamental parameters. It is based on a regularized sliced inverse regression (RSIR). We first tested it on noisy synthetic spectra of A, F, G, and K-type stars, and inverted simultaneously their atmospheric fundamental parameters: Teff, log g, [M/H] and vsini. Different learning databases were calculated using a range of sampling in Teff, log g, vsini, and [M/H]. Combined with a principal component analysis (PCA) nearest neighbors (NN) search, the size of the learning database is reduced. A Tikhonov regularization is applied, given the ill-conditioning of SIR. For all spectral types, decreasing the size of the learning database allowed us to reach internal accuracies better than PCA-based NN-search using larger learning databases. For each analyzed parameter, we have reached internal errors that are smaller than the sampling step of the parameter. We have also applied the technique to a sample of observed FGK and A stars. For a selection of well studied stars, the inverted parameters are in agreement with the ones derived in previous studies. The RSIR inversion technique, complemented with PCA pre-processing proves to be efficient in estimating stellar parameters of A, F, G, and K stars.

This short review aims to clarify upon the origins of the so-called Eddington-Barbier relationships, which relate the emergent specific intensity and the flux to the photospheric source function at specific optical depths. Here we discuss about the assumptions behind the original derivation of Barbier (1943). We also point to the fact that Milne had already formulated these two relations in 1921.

We aim at finding the value of an explanatory variable, through its expression in a large data-vector, without knowing the link function between the explanatory variable and the data-space. Sliced Inverse Regression (SIR) method allows for the projection of a data-vector onto a subspace consistent with the explanatory variable variation. We suggest a method based on the SIR subspace, that gives the most efficient estimation of an unknown explanatory variable.

Our aim is to evaluate fundamental parameters from the analysis of the electromagnetic spectra of stars. We may use $10^3$-$10^5$ spectra; each spectrum being a vector with $10^2$-$10^4$ coordinates. We thus face the so-called "curse of dimensionality". We look for a method to reduce the size of this data-space, keeping only the most relevant information.As a reference method, we use principal component analysis (PCA) to reduce dimensionality. However, PCA is an unsupervised method, therefore its subspace was not consistent with the parameter. We thus tested a supervised method based on Sliced Inverse Regression (SIR), which provides a subspace consistent with the parameter. It also shares analogies with factorial discriminant analysis: the method slices the database along the parameter variation, and builds the subspace which maximizes the inter-slice variance, while standardizing the total projected variance of the data. Nevertheless the performances of SIR were not satisfying in standard usage, because of the non-monotonicity of the unknown function linking the data to the parameter and because of the noise propagation. We show that better performances can be achieved by selecting the most relevant directions for parameter inference. Preliminary tests are performed on synthetic pseudo-line profiles plus noise. Using one direction, we show that compared to PCA, the error associated with SIR is 50$\%$ smaller on a non-linear parameter, and 70$\%$ smaler on a linear parameter. Moreover, using a selected direction, the error is 80$\%$ smaller for a non-linear parameter, and 95$\%$ smaller for a linear parameter.

We study the rotation-activity correlations (RACs) in a sample stars from spectral type dK4 to dM4. We study RACs using chromospheric data and coronal data. We study the Ca\,{\sc ii} line surface fluxes-$P/\sin i$ RACs. We fit the RACs with linear homoscedastic and heteroscedastic regression models. We find that these RACs differ substantially from one spectral sub-type to another. For dM3 and dM4 stars, we find that the RACs cannot be described by a simple model, but instead that there may exist two distinct RAC behaviors for the low activity and the high activity stellar sub-samples respectively. Although these results are preliminary and will need confirmation, the data suggest that these distinct RACs may be associated with different dynamo regimes. We also study R'$_{HK}$ as a function of the Rossby number R$_{0}$. We find that: (i) For dK4 stars, we confirm R'$_{HK}$ as a function of R$_{0}$ agrees well with previous results for F-G-K stars. (ii) In dK6, dM2, dM3 and dM4 stars, we find that, at a given R$_{0}$, the values of R'$_{HK}$ lie a factor of 3, 10, 20 and 90 respectively below the F-G-K RAC. Our results suggest a significant decrease in the efficiency of the dynamo mechanism(s) as regards chromospheric heating before and at dM3, i.e. before and at the TTCC. We also show that the ratio of coronal heating to chromospheric heating $L_{X}/L_{HK}$ increases by a factor of 100 between dK4 and dM4 stars.

Betelgeuse is an M supergiant that harbors spots and giant granules at its surface and presents linear polarization of its continuum. We have previously discovered linear polarization signatures associated with individual lines in the spectra of cool and evolved stars. Here, we investigate whether a similar linearly polarized spectrum exists for Betelgeuse. We used the spectropolarimeter Narval, combining multiple polarimetric sequences to obtain high signal-to-noise ratio spectra of individual lines, as well as the least-squares deconvolution (LSD) approach. We have discovered the existence of a linearly polarized spectrum for Betelgeuse, detecting a rather strong signal (at a few times 10$^{-4}$ of the continuum intensity level), both in individual lines and in the LSD profiles. Studying its properties and the signal observed for the resonant \ion{Na}{i}\,D lines, we conclude that we are mainly observing depolarization of the continuum by the absorption lines. The linear polarization of the Betelgeuse continuum is due to the anisotropy of the radiation field induced by brightness spots at the surface and Rayleigh scattering in the atmosphere. We have developed a geometrical model to interpret the observed polarization, from which we infer the presence of two brightness spots and their positions on the surface of Betelgeuse. We show that applying the model to each velocity bin along the Stokes Q and U profiles allows the derivation of a map of the bright spots. We use the Narval linear polarization observations of Betelgeuse obtained over a period of 1.4 years to study the evolution of the spots. Our study of the linearly polarized spectrum of Betelgeuse provides a novel method for studying the evolution of brightness spots at its surface and complements quasi-simultaneous observations obtained with PIONIER at the VLTI.

Reliable determination of Rotation-Activity Correlation (RAC's) depends on precise measurements of the following stellar parameters: $T_{eff}$, parallax, radius, metallicity, and rotational speed $v\sin i$. In this paper, our goal is to focus on the determination of these parameters for a sample of K and M dwarfs. In a future paper (Paper II), we will combine our rotational data with activity data in order to construct RAC's. Here, we report on a determination of effective temperatures based on the (R-I)$_C$ color from the calibrations of Mann et al. (2015) and Kenyon \& Hartmann (1995) for four samples of late-K, dM2, dM3 and dM4 stars. We also determine stellar parameters ($T_{eff}$, $log(g)$ and [M/H]) using the PCA-based inversion technique for a sample of 105 late-K dwarfs. We compile all effective temperatures from the literature for this sample. We determine empirical radius-[M/H] correlations in our stellar samples. This allows us to propose new effective temperatures, stellar radii, and metallicities for a large sample of 612 late-K and M dwarfs. Our mean radii agree well with those of Boyajian et al. (2012). We analyze HARPS and SOPHIE spectra of 105 late-K dwarfs, and have detected $v\sin i$ in 92 stars. In combination with our previous $v\sin i$ measurements in M and K dwarfs, we now derive $P/sin i$ measures for a sample of 418 K and M dwarfs. We investigate the distributions of $P/sin i$ and we show that they are different from one spectral sub-type to another at a 99.9\% confidence level.

We present an automated procedure that derives simultaneously the effective temperature $T_{eff}$, the surface gravity logg, the metallicity [Fe/H], and the equatorial projected rotational velocity vsini for "normal" A and Am stars. The procedure is based on the principal component analysis inversion method of Paletou et al. (2015a). A sample of 322 high resolution spectra of F0-B9 stars, retrieved from the Polarbase, SOPHIE, and ELODIE databases, were used to test this technique with real data. We have selected the spectral region from 4400-5000\AA\ as it contains many metallic lines and the Balmer H$\beta$ line. Using 3 datasets at resolving powers of R=42000, 65000 and 76000, about 6.6x$10^6$ synthetic spectra were calculated to build a large learning database. The Online Power Iteration algorithm was applied to these learning datasets to estimate the principal components (PC). The projection of spectra onto the few PCs offered an efficient comparison metric in a low dimensional space. The spectra of the well known A0- and A1-type stars, Vega and Sirius A, were used as control spectra in the three databases. Spectra of other well known A-type stars were also employed in order to characterize the accuracy of the inversion technique. All observational spectra were inverted and atmospheric parameters derived. After removal of a few outliers, the PCA-inversion method appears to be very efficient in determining $T_{eff}$, [Fe/H], and vsini for A/Am stars. The derived parameters agree very well with previous determinations. Using a statistical approach, deviations of around 150 K, 0.35 dex, 0.15 dex, and 2 km/s were found for $T_{eff}$, logg, [Fe/H], and vsini with respect to literature values for A-type stars. The PCA-inversion proves to be a very fast, practical, and reliable tool for estimating stellar parameters of FGK and A stars, and deriving effective temperatures of M stars.

This article presents an on-line tool (rttools.irap.omp.eu) and its accompanying software ressources for the numerical solution of basic radiation transfer out of local thermodynamic equilibrium (LTE). State-of-the-art stationary iterative methods such as Accelerated $\Lambda$-Iteration and Gauss-Seidel schemes, using a short characteristics-based formal solver are used. We also comment on typical numerical experiments associated to the basic non-LTE radiation problem. These ressources are intended for the largest use and benefit, in support to more classical radiation transfer lectures usually given at the Master level.

This work addresses a procedure to estimate fundamental stellar parameters such as T eff , logg, [Fe/H], and v sin i using a dimensionality reduction technique called Principal Component Analysis (PCA), applied to a large database of synthetic spectra. This technique shows promising results for inverting stellar parameters of observed targets from Gaia ESO Survey.

Hereafter we describe the activities of the $Grand \, Sud-Ouest$ Data Centre operated for INSU/CNRS by the OMP-IRAP and the Universit\'e Paul Sabatier (Toulouse), in a collaboration with the OASU-LAB (Bordeaux) and OREME-LUPM (Montpellier).

We show how the range of application of the principal component analysis-based inversion method of Paletou et al. (2015) can be extended to late-type stars data. Besides being an extension of its original application domain, for FGK stars, we also used synthetic spectra for our learning database. We discuss our results on effective temperatures against previous evaluations made available from Vizier and Simbad services at CDS.

The theoretically studied impact of rapid rotation on stellar evolution needs to be confronted with the results of high resolution spectroscopy-velocimetry observations. A weak surface magnetic field had recently been detected in the A0 prototype star Vega, potentially leading to a (yet undetected) structured surface. The goal of this article is to present a thorough analysis of the line profile variations and associated estimators in the early-type standard star Vega (A0) in order reveal potential activity tracers, exoplanet companions and stellar oscillations. Vega was monitored in high-resolution spectroscopy with the velocimeter Sophie/OHP. A total of 2588 high S/N spectra was obtained during 5 nights (August 2012) at R = 75000 and covering the visible domain. For each reduced spectrum, Least Square Deconvolved (LSD) equivalent photospheric profiles were calculated with a Teff = 9500 and logg = 4.0 spectral line mask. Several methods were applied to study the dynamic behavior of the profile variations (evolution of radial velocity, bisectors, vspan, 2D profiles, amongst others). We present the discovery of a starspotted stellar surface in an A-type standard star with faint spot amplitudes Delta F/Fc ~5 10^{-4}. A rotational modulation of spectral lines with a period of rotation P = 0.68 d has clearly been exhibited, confirming the results of previous spectropolarimetric studies. Either a very thin convective layer can be responsible for magnetic field generation at small amplitudes, or a new mechanism has to be invoked in order to explain the existence of activity tracing starspots. This first strong evidence that standard A-type stars can show surface structures opens a new field of research and asks the question about a potential link with the recently discovered weak magnetic field discoveries in this category of stars.

The general context of this study is the inversion of stellar fundamental parameters from high-resolution Echelle spectra. We aim indeed at developing a fast and reliable tool for the post-processing of spectra produced by Espadons and Narval spectropolarimeters. Our inversion tool relies on principal component analysis. It allows reduction of dimensionality and the definition of a specific metric for the search of nearest neighbours between an observed spectrum and a set of observed spectra taken from the Elodie stellar library. Effective temperature, surface gravity, total metallicity and projected rotational velocity are derived. Various tests presented in this study, and done from the sole information coming from a spectral band centered around the Mg I b-triplet and with spectra from FGK stars are very promising.

This tutorial is devoted to extending an existing catalogue with data taken elsewhere, either from CDS Vizier or Simbad database. As an example, we used the so-called 'Spectroscopic Survey of Stars in the Solar Neighborhood' (aka. S4N, Allende Prieto et al. 2004) in order to retrieve all objects with available data for the set of fundamental stellar parameters effective temperature, surface gravity and metallicity. Then for each object in this dataset we query Simbad database to retrieve the projected rotational velocity. This combines Vizier and Simbad queries made using Python astroquery module. The tutorial covers remote database access, filtering tables with arbitrary criteria, creating and writing your own tables, and basics of plotting in Python.

The general context of this study is the inversion of stellar fundamental parameters from high-resolution Echelle spectra. We aim at developing a fast and reliable tool for the post-processing of spectra produced, in particular, by the Espadons and Narval spectropolarimeters. Our inversion tool relies on principal component analysis. It allows reduction of dimensionality and the definition of a specific metric for the search of nearest neighbours between an observed spectrum and a set of synthetic spectra. Effective temperature, surface gravity, total metallicity and projected rotation velocity are derived. Our first tests, essentially done from the sole information coming from the spectral band that the RVS spectrometer will soon observe from the GAIA space observatory, and with spectra from mainly FGK-dwarfs are very promising. We also tested our method with a few targets beyond this domain of the H-R diagram.

PolarBase is an evolving data base that contains all stellar data collected with the ESPaDOnS and NARVAL high-resolution spectropolarimeters, in their reduced form, as soon as they become public. As of early 2014, observations of 2,000 stellar objects throughout the Hertzsprung-Russell diagram are available. Intensity spectra are available for all targets, and the majority of the observations also include simultaneous spectra in circular or linear polarization, with the majority of the polarimetric measurements being performed only in circularly polarized light (Stokes V). Observations are associated with a cross-correlation pseudo-line profile in all available Stokes parameters, greatly increasing the detectability of weak polarized signatures. Stokes V signatures are detected for more than 300 stars of all masses and evolutionary stages, and linear polarization is detected in 35 targets. The detection rate in Stokes V is found to be anti-correlated with the stellar effective temperature. This unique set of Zeeman detections offers the first opportunity to run homogeneous magnetometry studies throughout the H-R diagram. The web interface of PolarBase is available at http://polarbase.irap.omp.eu.

The general context of this study concerns the post-processing of multiline spectropolarimetric observations of stars, and in particular these numerical analysis techniques aiming at the detection and the characterization of polarized signatures. Hereafter, using real observational data, we compare and clarify a number of points concerning various methods of analysis. Indeed, simple line addition, least-squares deconvolution and denoising by principal component analysis have been applied, and compared to each other, to polarized stellar spectra available from the TBLegacy database of the Narval spectropolarimeter. Such a comparison between various approaches of distinct sophistication levels allows us to make a safe choice for the next implementation of on-line post-processing of our unique database for the stellar physics community.

Aims. We aim to investigate the long-term temporal evolution of the magnetic field of the solar-type star \xi Bootis A, both from direct magnetic field measurements and from the simultaneous estimate of indirect activity indicators. Methods. We obtained seven epochs of high-resolution, circularly-polarized spectra from the NARVAL spectropolarimeter between 2007 and 2011, for a total of 76 spectra. Using approximately 6,100 photospheric spectral lines covering the visible domain, we employed a cross-correlation procedure to compute a mean polarized line profile from each spectrum. The large-scale photospheric magnetic field of the star was then modelled by means of Zeeman-Doppler Imaging, allowing us to follow the year-to-year evolution of the reconstructed magnetic topology. Simultaneously, we monitored the width of several magnetically sensitive spectral lines, the radial velocity, the line asymmetry of intensity line profiles, and the chromospheric emission in the cores of the Ca II H and Halpha lines. Results. During the highest observed activity states, in 2007 and 2011, the large-scale field of \xi Boo A is almost completely axisymmetric and is dominated by its toroidal component. The magnetic topologies reconstructed for these activity maxima are very similar, suggesting a form of short cyclicity in the large-scale field distribution. Correlated temporal evolution, due to both rotational modulation and seasonal variability, is observed between the Ca II emission, the Halpha emission and the width of magnetically sensitive lines. When measurable, the differential rotation reveals a strong latitudinal shear in excess of 0.2 rad/d.

We aim at a highly sensitive search for weak magnetic fields in main sequence stars of intermediate mass, by scanning classes of stars with no previously reported magnetic members. After detecting a weak magnetic field on the normal, rapidly rotating A-type star Vega, we concentrate here on the bright star Sirius A, taken as a prototypical, chemically peculiar, moderately rotating Am star. We employed the NARVAL and ESPaDOnS high-resolution spectropolarimeters to collect 442 circularly polarized spectra, complemented by 60 linearly polarized spectra. Using a list of about 1,100 photospheric spectral lines, we computed a cross correlation line profile from every spectrum, leading to a signal-to-noise ratio of up to 30,000 in the polarized profile. We report the repeated detection of circularly polarized, highly asymmetric signatures in the line profiles, interpreted as Zeeman signatures of a large-scale photospheric magnetic field, with a line-of-sight component equal to $0.2 \pm 0.1$ G. This is the first polarimetric detection of a surface magnetic field on an Am star. Using rough estimates of the physical properties of the upper layers of Sirius A, we suggest that a dynamo operating in the shallow convective envelope cannot account for the field strength reported here. Together with the magnetic field of Vega, this result confirms that a new class of magnetic objects exists among non Ap/Bp intermediate-mass stars, and it may indicate that a significant fraction of tepid stars are magnetic.

This study concerns the fast and accurate solution of multilevel non-LTE radiation transfer problems. We propose and evaluate an alternative iterative scheme to the classical MALI method. Our study is indeed based on the application of Broyden's method for the solution of nonlinear systems of equations. Comparative tests, in 1D plane-parallel geometry, between the popular MALI method and our alternative method are discussed. The Broyden method is typically 4.5 times faster than MALI. It makes it also fairly competitive with Gauss-Seidel and Successive Over-Relaxation methods developed after MALI.

Magnetic fields of cool stars can be directly investigated through the study of the Zeeman effect on photospheric spectral lines using several approaches. With spectroscopic measurement in unpolarised light, the total magnetic flux averaged over the stellar disc can be derived but very little information on the field geometry is available. Spectropolarimetry provides a complementary information on the large-scale component of the magnetic topology. With Zeeman-Doppler Imaging (ZDI), this information can be retrieved to produce a map of the vector magnetic field at the surface of the star, and in particular to assess the relative importance of the poloidal and toroidal components as well as the degree of axisymmetry of the field distribution. The development of high-performance spectropolarimeters associated with multi-lines techniques and ZDI allows us to explore magnetic topologies throughout the Hertzsprung-Russel diagram, on stars spanning a wide range of mass, age and rotation period. These observations bring novel constraints on magnetic field generation by dynamo effect in cool stars. In particular, the study of solar twins brings new insight on the impact of rotation on the solar dynamo, whereas the detection of strong and stable dipolar magnetic fields on fully convective stars questions the precise role of the tachocline in this process.

This study concerns the fast and accurate solution of the line radiation transfer problem, under non-LTE conditions. We propose and evaluate an alternative iterative scheme to the classical ALI-Jacobi method, and to the more recently proposed Gauss-Seidel and Successive Over-Relaxation (GS/SOR) schemes. Our study is indeed based on the application of a preconditioned bi-conjugate gradient method (BiCG-P). Standard tests, in 1D plane parallel geometry and in the frame of the two-level atom model, with monochromatic scattering, are discussed. Rates of convergence between the previously mentioned iterative schemes are compared, as well as their respective timing properties. The smoothing capability of the BiCG-P method is also demonstrated.

A robust numerical method called the Preconditioned Bi-Conjugate Gradient (Pre-BiCG)method is proposed for the solution of radiative transfer equation in spherical geometry.A variant of this method called Stabilized Preconditioned Bi-Conjugate Gradient (Pre-BiCG-STAB) is also presented. These are iterative methods based on the construction of a set of bi-orthogonal vectors. The application of Pre-BiCG method in some benchmark tests show that the method is quite versatile, and can handle hard problems that may arise in astrophysical radiative transfer theory.

We present test results on the simultaneous inversion of the Stokes profiles of the He I lines at 587.6 nm (D_3) and 1083.0 nm in prominences (90-deg scattering). We created datasets of synthetic Stokes profiles for the case of quiescent prominences (B<200 G), assuming a conservative value of 10^-3 of the peak intensity for the polarimetric sensitivity of the simulated observations. In this work, we focus on the error analysis for the inference of the magnetic field vector, under the usual assumption that the prominence can be assimilated to a slab of finite optical thickness with uniform magnetic and thermodynamic properties. We find that the simultaneous inversion of the two lines significantly reduces the errors on the inference of the magnetic field vector, with respect to the case of single-line inversion. These results provide a solid justification for current and future instrumental efforts with multi-line capabilities for the observations of solar prominences and filaments.

The diagnosis of new high-resolution spectropolarimetric observations of solar prominences made in the visible and near-infrared mainly, requires a radiative modelling taking into account for both multi-dimensional geometry and complex atomic models. Hereafter we contribute to the improvement of the diagnosis based on the observation of He I multiplets, by considering 2D non-LTE unpolarized radiation transfer, and taking also into account the atomic fine structure of helium. It is an improvement and a direct application of the multi-grid Gauss-Seidel/SOR iterative scheme in 2D cartesian geometry developed by us. It allows us to compute realistic emergent intensity profiles for the He I 10830 A and D3 multiplets, which can be directly compared to the simultaneous and high-resolution observations made at THeMIS. A preliminary 2D multi-thread modelling is also discussed.

In his famous monographs, Einar Tandberg-Hanssen writes that "the single, physically most important parameter to study in prominences may be the magnetic field. Shapes, motions, and in fact the very existence of prominences depend on the nature of the magnetic field threading the prominence plasma". Hereafter we sumarize recent contributions and advances in our knowledge about the magnetic field of solar prominences. It mostly relies on high resolution and high sensitivity spectropolarimetry made both in the visible and in the near infrared.

We present here additional results of a spectropolarimetric survey of a small sample of stars ranging from spectral type M0 to M8 aimed at investigating observationally how dynamo processes operate in stars on both sides of the full convection threshold (spectral type M4). The present paper focuses on early M stars (M0--M3), i.e. above the full convection threshold. Applying tomographic imaging techniques to time series of rotationally modulated circularly polarised profiles collected with the NARVAL spectropolarimeter, we determine the rotation period and reconstruct the large-scale magnetic topologies of 6 early M dwarfs. We find that early-M stars preferentially host large-scale fields with dominantly toroidal and non-axisymmetric poloidal configurations, along with significant differential rotation (and long-term variability); only the lowest-mass star of our subsample is found to host an almost fully poloidal, mainly axisymmetric large-scale field ressembling those found in mid-M dwarfs. This abrupt change in the large-scale magnetic topologies of M dwarfs (occuring at spectral type M3) has no related signature on X-ray luminosities (measuring the total amount of magnetic flux); it thus suggests that underlying dynamo processes become more efficient at producing large-scale fields (despite producing the same flux) at spectral types later than M3. We suspect that this change relates to the rapid decrease in the radiative cores of low-mass stars and to the simultaneous sharp increase of the convective turnover times (with decreasing stellar mass) that models predict to occur at M3; it may also be (at least partly) responsible for the reduced magnetic braking reported for fully-convective stars.

We present in this paper the first results of a spectropolarimetric analysis of a small sample (~ 20) of active stars ranging from spectral type M0 to M8, which are either fully-convective or possess a very small radiative core. This study aims at providing new constraints on dynamo processes in fully-convective stars. The present paper focuses on 5 stars of spectral type ~M4, i.e. with masses close to the full convection threshold (~ 0.35 Msun), and with short rotational periods. Tomographic imaging techniques allow us to reconstruct the surface magnetic topologies from the rotationally modulated time-series of circularly polarised profiles. We fnd that all stars host mainly axisymmetric large-scale poloidal fields. Three stars were observed at two different epochs separated by ~1 yr; we find the magnetic topologies to be globally stable on this timescale. We also provide an accurate estimation of the rotational period of all stars, thus allowing us to start studying how rotation impacts the large-scale magnetic field.

We present preliminary results of 2D radiative modelling of He I lines in solar prominences, using a new numerical code developed by us (Leger, Chevallier and Paletou 2007). It treats self-consistently the radiation transfer and the non-LTE statistical equilibrium of H and, in a second stage, the one of He using a detailed atomic model. Preliminary comparisons with new visible plus near-infrared observations made at high spectral resolution with THeMIS are very satisfactory.

From a set of stellar spectropolarimetric observations, we report the detection of surface magnetic fields in a sample of four solar-type stars, namely HD 73350, HD 76151, HD 146233 (18 Sco) and HD 190771. Assuming that the observed variability of polarimetric signal is controlled by stellar rotation, we establish the rotation periods of our targets, with values ranging from 8.8 d (for HD 190771) to 22.7 d (for HD 146233). Apart from rotation, fundamental parameters of the selected objects are very close to the Sun's, making this sample a practical basis to investigate the specific impact of rotation on magnetic properties of Sun-like stars. We reconstruct the large-scale magnetic geometry of the targets as a low-order (l<10) spherical harmonics expansion of the surface magnetic field. From the set of magnetic maps, we draw two main conclusions. (a) The magnetic energy of the large-scale field increases with rotation rate. The increase of chromospheric emission with the mean magnetic field is flatter than observed in the Sun. Since the chromospheric flux is also sensitive to magnetic elements smaller than those contributing to the polarimetric signal, this observation suggests that a larger fraction of the surface magnetic energy is stored in large scales as rotation increases. (b) Whereas the magnetic field is mostly poloidal for low rotation rates, more rapid rotators host a large-scale toroidal component in their surface field. From our observations, we infer that a rotation period lower than ~12 days is necessary for the toroidal magnetic energy to dominate over the poloidal component.

From observations collected with the ESPaDOnS and NARVAL spectropolarimeters, we report the detection of Zeeman signatures on the classical T Tauri star BP Tau. Circular polarisation signatures in photospheric lines and in narrow emission lines tracing magnetospheric accretion are monitored throughout most of the rotation cycle of BP Tau at two different epochs in 2006. We observe that rotational modulation dominates the temporal variations of both unpolarised and circularly polarised spectral proxies tracing the photosphere and the footpoints of accretion funnels. From the complete data sets at each epoch, we reconstruct the large-scale magnetic topology and the location of accretion spots at the surface of BP Tau using tomographic imaging. We find that the field of BP Tau involves a 1.2 kG dipole and 1.6 kG octupole, both slightly tilted with respect to the rotation axis. Accretion spots coincide with the two main magnetic poles at high latitudes and overlap with dark photospheric spots; they cover about 2% of the stellar surface. The strong mainly-axisymmetric poloidal field of BP Tau is very reminiscent of magnetic topologies of fully-convective dwarfs. It suggests that magnetic fields of fully-convective cTTSs such as BP Tau are likely not fossil remants, but rather result from vigorous dynamo action operating within the bulk of their convective zones. Preliminary modelling suggests that the magnetosphere of BP Tau extends to distances of at least 4 R* to ensure that accretion spots are located at high latitudes, and is not blown open close to the surface by a putative stellar wind. It apparently succeeds in coupling to the accretion disc as far out as the corotation radius, and could possibly explain the slow rotation of BP Tau.

A small fraction of intermediate-mass main sequence (A and B type) stars have strong, organised magnetic fields. The large majority of such stars, however, show no evidence for magnetic fields, even when observed with very high precision. In this paper we describe a simple model, motivated by qualitatively new observational results, that provides a natural physical explanation for the small fraction of observed magnetic stars.

BASS 2000 is the French solar database for ground-based instruments. We describe hereafter our organization, our tasks and the products we can deliver to the international community. Our prospects cover data mining into the THeMIS archive, a participation to the EST endeavour and the creation and curation of the ESPaDOnS/NARVAL stellar spectra database.

Most of the solar physicists use very expensive software for data reduction and visualization. We present hereafter a reliable freeware solution based on the Python language. This is made possible by the association of the latter with a small set of additional libraries developed in the scientific community. It provides then a very powerful and economical alternative to other interactive data languages. Although it can also be used for any kind of post-processing of data, we demonstrate the capabities of such a set of freeware tools using THeMIS observations of the second solar spectrum.

We have investigated a sample of 28 well-known spectroscopically-identified magnetic Ap/Bp stars, with weak, poorly-determined or previously undetected magnetic fields, with the aim of exploring the weak part of the magnetic field distribution of Ap/Bp stars. Using the MuSiCoS and NARVAL spectropolarimeters we have obtained 282 LSD Stokes V signatures of our 28 sample stars. All stars were detected, showing clearly that when observed with sufficient precision, all firmly classified Ap/Bp stars show detectable surface magnetic fields. To better characterise the surface magnetic field intensities and geometries of the sample, we have inferred the dipolar field intensity and the magnetic obliquity. The distribution of derived dipole strengths for these stars exhibits a plateau at about 1 kG, falling off to larger and smaller field strengths. Remarkably, in this sample of stars selected for their presumably weak magnetic fields, we find only 2 stars for which the derived dipole strength is weaker than 300 G. We interpret this "magnetic threshold" as a critical value necessary for the stability of large-scale magnetic fields, and develop a simple quantitative model that is able to approximately reproduce the observed threshold characteristics. This scenario leads to a natural explanation of the small fraction of intermediate-mass magnetic stars. It may also explain the near-absence of magnetic fields in more massive B and O-type stars.

From observations collected with the ESPaDOnS spectropolarimeter, we report the discovery of magnetic fields at the surface of the mildly accreting classical T Tauri star V2129 Oph. Zeeman signatures are detected, both in photospheric lines and in the emission lines formed at the base of the accretion funnels linking the disc to the protostar, and monitored over the whole rotation cycle of V2129 Oph. We observe that rotational modulation dominates the temporal variations of both unpolarized and circularly polarized line profiles. We reconstruct the large-scale magnetic topology at the surface of V2129 Oph from both sets of Zeeman signatures simultaneously. We find it to be rather complex, with a dominant octupolar component and a weak dipole of strengths 1.2 and 0.35 kG, respectively, both slightly tilted with respect to the rotation axis. The large-scale field is anchored in a pair of 2-kG unipolar radial field spots located at high latitudes and coinciding with cool dark polar spots at photospheric level. This large-scale field geometry is unusually complex compared to those of non-accreting cool active subgiants with moderate rotation rates. As an illustration, we provide a first attempt at modelling the magnetospheric topology and accretion funnels of V2129 Oph using field extrapolation. We find that the magnetosphere of V2129 Oph must extend to about 7R* to ensure that the footpoints of accretion funnels coincide with the high-latitude accretion spots on the stellar surface. It suggests that the stellar magnetic field succeeds in coupling to the accretion disc as far out as the corotation radius, and could possibly explain the slow rotation of V2129 Oph. The magnetospheric geometry we derive produces X-ray coronal fluxes typical of those observed in cTTSs.

New high-resolution spectropolarimetric observations of solar prominences require improved radiative modelling capabilities in order to take into account both multi-dimensional - at least 2D - geometry and complex atomic models. This makes necessary the use of very fast numerical schemes for the resolution of 2D non-LTE radiative transfer problems considering freestanding and illuminated slabs. The implementation of Gauss-Seidel and successive over-relaxation iterative schemes in 2D, together with a multi-grid algorithm, is thoroughly described in the frame of the short characteristics method for the computation of the formal solution of the radiative transfer equation in cartesian geometry. We propose a new test for multidimensional radiative transfer codes and we also provide original benchmark results for simple 2D multilevel atom cases which should be helpful for the further development of such radiative transfer codes, in general.

Stellar magnetic fields can be investigated using several, very complementary approaches. While conventional spectroscopy is capable of estimating the average magnetic strength of potentially complex field configurations thanks to its low sensitivity to the vector properties of the field, spectropolarimetry can be used to map the medium- and large-scale structure of magnetic topologies. In particular, the latter approach allows one to retrieve information about the poloidal and toroidal components of the large-scale dynamo fields in low-mass stars, and thus to investigate the physical processes that produce them. Similarly, this technique can be used to explore how magnetic fields couple young stars to their massive accretion disc and thus to estimate how much mass and angular momentum are transfered to the newly-born low-mass star. We present here the latest results in this field obtained with spectropolarimetry, with special emphasis on the surprising discoveries obtained on very-low mass fully-convective stars and classical T Tauri stars thanks to the ESPaDOnS spectropolarimeter recently installed on the 3.6m Canada-France-Hawaii Telescope.

The vast majority of recent advances in the field of numerical radiative transfer relies on approximate operator methods better known in astrophysics as Accelerated Lambda-Iteration (ALI). A superior class of iterative schemes, in term of rates of convergence, such as Gauss-Seidel and Successive Overrelaxation methods were therefore quite naturally introduced in the field of radiative transfer by Trujillo Bueno & Fabiani Bendicho (1995); it was thoroughly described for the non-LTE two-level atom case. We describe hereafter in details how such methods can be generalized when dealing with non-LTE unpolarised radiation transfer with multilevel atomic models, in monodimensional geometry.

We report the discovery of a medium-strength (~0.5kG) magnetic field on the young, massive star tauSco (B0.2V), which becomes the third-hottest magnetic star known. Circularly polarized Zeeman signatures are clearly detected in observations collected mostly with the ESPaDOnS spectropolarimeter, recently installed on the 3.6-m Canada-France-Hawaii Telescope; temporal variability is also clearly established in the polarimetry, and can be unambiguously attributed to rotational modulation with a period close to 41d. Archival UV spectra confirm that this modulation repeats over timescales of decades. By reconstructing the large-scale structure of its magnetic topology, we find that the magnetic structure is unusually complex for a hot star. The surface topology is dominated by a potential field, although a moderate toroidal component is probably present. We fail to detect intrinsic temporal variability of the magnetic structure over the 1.5-yr period of our spectropolarimetric observations (in agreement with the stable temporal variations of the UV spectra), and infer that any differential surface rotation must be very small. The topology of the extended magnetic field that we derive from the photospheric magnetic maps is also more complex than a global dipole, and features in particular a significantly warped torus of closed magnetic loops encircling the star (tilted at about 90deg to the rotation axis), with additional, smaller, networks of closed field lines. This topology appears to be consistent with the exceptional Xray properties of tauSco and also provides a natural explanation of the variability observed in wind-formed UV lines. We conclude that its magnetic field is most probably a fossil remnant from the star-formation stage.

We test the accuracy of the ALI method, widely used in stellar atmosphere modelling, by solving exactly a standard radiative transfer problem in plane-parallel geometry. Some recommendations are given for a practical use of this method.

Although published in 1995, the Gauss-Seidel method for solving the non-LTE radiative transfer problem has deserved too little attention in the astrophysical community yet. Further tests of the performances and of the accuracy of the numerical scheme are provided.

We solve the integral equation describing the propagation of light in an isothermal plane-parallel atmosphere of optical thickness $\tau^*$, adopting a uniform thermalization parameter $\epsilon$. The solution given by the ALI method, widely used in the field of stellar atmospheres modelling, is compared to the exact solution. Graphs are given that illustrate the accuracy of the ALI solution as a function of the parameters $\epsilon$, $\tau^*$ and optical depth variable $\tau$.

Models predict that magnetic fields play a crucial role in the physics of astrophysical accretion disks and their associated winds and jets. For example, the rotation of the disk twists around the rotation axis the initially vertical magnetic field, which responds by slowing down the plasma in the disk and by causing it to fall towards the central star. The magnetic energy flux produced in this process points away from the disk, pushing the surface plasma outwards, leading to a wind from the disk and sometimes a collimated jet. But these predictions have hitherto not been supported by observations. Here we report the direct detection of the magnetic field in the core of the protostellar accretion disk FU Orionis. The surface field reaches strengths of about 1 kG close to the centre of the disk, and it includes a significant azimuthal component, in good agreement with recent models. But we find that the field is very filamentary and slows down the disk plasma much more than models predict, which may explain why FU Ori fails to collimate its wind into a jet.

We investigate the magnetic geometry of the active G8 dwarf xi Bootis A, from spectropolarimetric observations obtained in 2003 with the MuSiCoS echelle spectropolarimeter at the Telescope Bernard Lyot (Observatoire du Pic du Midi, France). We repeatedly detect a photospheric magnetic field, with periodic variations consistent with rotational modulation. Circularly polarized (Stokes V) line profiles present a systematic asymmetry, showing up as an excess in amplitude and area of the blue lobe of the profiles. A direct modeling of Stokes V profiles suggests that the global magnetic field is composed of two main components, with an inclined dipole and a large-scale toroidal field. We derive a dipole intensity of about 40 G, with an inclination of 35 degrees of the dipole with respect to the rotation axis. The toroidal field strength is of order of 120 G. A noticeable evolution of the field geometry is observed over the 40 nights of our observing window and results in an increase of the field strength and of the dipole inclination. This study is the first step of a long-term monitoring of xi Bootis A and other active solar-type stars, with the aim to investigate secular fluctuations of stellar magnetic geometries induced by activity cycles.

We present spectropolarimetric observations of the FK Com star HD 199178 obtained between 1998 December and 2003 August at the Telescope Bernard Lyot (Observatoire du Pic du Midi, France). We report the detection of a photospheric magnetic field and reconstruct its distribution by means of Zeeman-Doppler Imaging. We observe large regions where the magnetic field is mainly azimuthal, suggesting that the dynamo processes generating the magnetic activity of HD 199178 may be active very close to the stellar surface. We investigate the rapid evolution of surface brightness and magnetic structures from a continuous monitoring of the star over several weeks in 2002 and 2003. We report that significant changes occur in the distribution of cool spots and magnetic regions on typical timescales of the order of 2 weeks. Our spectropolarimetric observations also suggest that the surface of HD 199178 is sheared by differential rotation, with a difference in rotation rate between equatorial and polar regions of the order of 1.5 times that of the Sun.

In the era of XMM-Newton and Chandra missions, it is crucial to use codes able to compute correctly the line spectrum of X-ray irradiated thick media (Thomson thickness of the order of unity) to build models for the structure and the emission of the central regions of AGN or X-ray binaries. In all photoionized codes except in our code Titan, the line intensities are computed with the "escape probability approximation". In its last version, Titan solves the transfer of a thousand lines and of the continuum with the ``Accelerated Lambda Iteration" method, which is one of the most efficient and most secure for line transfer. We find that for conditions typical of the AGN or X-ray binary emission medium, all escape approximations commonly used lead to an overestimation of the soft X-ray lines which can reach one order of magnitude for intense lines.