The correct interpretation of the large amount of complex data from next-generation (in particular, space-based) observational facilities requires a very strong theoretical underpinning. One can thus predict that in the near future the use of atmospheric models obtained through the use of three-dimensional (3-D) radiation magnetohydrodynamics (RMHD) codes, coupled with advanced radiative transfer treatment including non-local thermodynamic equilibrium (non-LTE) effects and polarisation, will become the norm. In particular, stellar brightness variability in cool (i. e., spectral type F, G, K, and M) stars can be induced by several different effects, besides pulsation. We briefly discuss some literature results and mention some of our recent progress. Finally, we attempt to peek into the future of understanding this important aspect of the life of stars.

Theory surrounding the origin of the dust-laden winds from evolved stars remains mired in controversy. Characterizing the formation loci and the dust distribution within approximately the first stellar radius above the surface is crucial for understanding the physics that underlie the mass-loss phenomenon. By exploiting interferometric polarimetry, we derive the fundamental parameters that govern the dust structure at the wind base of a red supergiant. We present near-infrared aperture-masking observations of Betelgeuse in polarimetric mode obtained with the NACO/SAMPol instrument. We used both parametric models and radiative transfer simulations to predict polarimetric differential visibility data and compared them to SPHERE/ZIMPOL measurements. Using a thin dust shell model, we report the discovery of a dust halo that is located at only 0.5 R$_{\star}$ above the photosphere (i.e. an inner radius of the dust halo of 1.5 R$_{\star}$). By fitting the data under the assumption of Mie scattering, we estimate the grain size and density for various dust species. By extrapolating to the visible wavelengths using radiative transfer simulations, we compare our model with SPHERE/ZIMPOL data and find that models based on dust mixtures that are dominated by forsterite are most favored. Such a close dusty atmosphere has profound implications for the dust formation mechanisms around red supergiants.

High spatial resolution is the key for the understanding various astrophysical phenomena. But even with the future E-ELT, single dish instruments are limited to a spatial resolution of about 4 mas in the visible. For the closest objects within our Galaxy most of the stellar photosphere remains smaller than 1 mas. With the success of long baseline interferometry these limitations were soom overcome. Today low and high resolution interferometric instruments on the VLTI and CHARA offer an immense range of astrophysical studies. Combining more telescopes and moving to visible wavelengths broadens the science cases even more. With the idea of developing strong science cases for a future visible interferometer, we organized a science group around the following topics: pre-main sequence and main sequence stars, fundamental parameters, asteroseismology and classical pulsating stars, evolved stars, massive stars, active galactic nuclei (AGNs) and imaging techniques. A meeting was organized on the 15th and 16th of January, 2015 in Nice with the support of the Action Specific in Haute Resolution Angulaire (ASHRA), the Programme National en Physique Stellaire (PNPS), the Lagrange Laboratory and the Observatoire de la Cote d'Azur, in order to present these cases and to discuss them further for future visible interferometers. This White Paper presents the outcome of the exchanges. This book is dedicated to the memory of our colleague Olivier Chesneau who passed away at the age of 41.

One of the necessary parameters needed for the interpretation of the light curves of transiting exoplanets or eclipsing binaries, as well as interferometric measurements of a star or microlensing events is how the intensity and polarization of a light change from the center to the limb. Scattering and absorption processes in stellar atmosphere affect both the center-to limb variation of intensity (CLVI) and polarization (CLVP). In this paper, we present a study of the CLVI and CLVP in continuum spectra considering different contributions of scattering and absorption opacity for different spectral type stars with spherical atmospheres. We solve the polarized radiative transfer equation in the presence of continuum scattering, considering spherical stellar model atmospheres. We developed two independent codes based on Feautrier and short characteristics methods to cross-check our results. We calculate the CLVI and CLVP in continuum for the Phoenix grid of spherical stellar model atmospheres for a range of $T_{eff} = 4000 - 7000 \rm K$, $\log g = 1.0 - 5.5$ and $\lambda = 4000 - 7000 \rm \AA$, which are tabulated and available at the CDS. For sub-giant and dwarf stars ($\log g = 3.0 - 4.5$), lower $\log g$ and lower $T_{eff}$ of a star lead to higher limb polarization of the star. For giant and supergiant stars ($\log g = 1.0 - 2.5$), the highest effective temperature yields the largest polarization. By decreasing of the $T_{eff}$ of a star down to $4500 - 5500 \rm K$ (depending on $\log g$) the limb polarization decreases and reaches a local minimum. It increases again down to $T_{eff}$ of $4000 \rm K$. For the most compact dwarf stars ($\log g = 5.0 - 5.5$) the limb polarization degree shows a maximum for models with $T_{eff}$ in the range $4200 - 4600 \rm K$ (depending on $\log g$) and decreases toward higher and lower temperatures.

We present results of numerical simulations of flux and linear polarization variations in transiting exoplanetary systems, caused by the host star disk symmetry breaking. We consider different configurations of planetary transits depending on orbital parameters. Starspot contribution to the polarized signal is also estimated. Applying the method to known systems and simulating observational conditions, a number of targets is selected where transit polarization effects could be detected. We investigate several principal benefits of the transit polarimetry, particularly, for determining orbital spatial orientation and distinguishing between grazing and near-grazing planets. Simulations show that polarization parameters are also sensitive to starspots, and they can be used to determine spot positions and sizes.

In radial velocity observations, a pair of extrasolar planets near a 2:1 orbital resonance can be misinterpreted as a single eccentric planet, if data are sparse and measurement precision insufficient to distinguish between these models. We determine the fraction of alleged single-planet RV detected systems for which a 2:1 resonant pair of planets is also a viable model and address the question of how the models can be disentangled. By simulation we quantified the mismatch arising from applying the wrong model. Model alternatives are illustrated using the supposed single-planet system HD 27894 for which we also study the dynamical stability of near-2:1 resonant solutions. From the data scatter around the fitted single-planet Keplerians, we find that for $74\% $ of the $254$ putative single-planet systems, a 2:1 resonant pair cannot be excluded as a viable model, since the error due to the wrong model is smaller than the scatter. For $187$ stars $\chi ^2$-probabilities can be used to reject the Keplerian models with a confidence of $95\% $ for $54\% $ of the stars and with $99.9\% $ for $39\% $ of the stars. For HD 27894 a considerable fit improvement is obtained when adding a low-mass planet near half the orbital period of the known Jovian planet. Dynamical analysis demonstrates that this system is stable when both planets are initially placed on circular orbits. For fully Keplerian orbits a stable system is only obtained if the eccentricity of the inner planet is constrained to $<0.3$. A large part of the allegedly RV detected single-planet systems should be scrutinized in order to determine the fraction of systems containing near-2:1 resonant pairs of planets. Knowing the abundance of such systems will allow us to revise the eccentricity distribution for extrasolar planets and provide direct constraints for planetary system formation.

Context. Scattering and absorption processes in stellar atmosphere affect the center-to-limb variations of the intensity (CLVI) and the linear polarization (CLVP) of stellar radiation. Aims. There are several theoretical and observational studies of CLVI using different stellar models, however, most studies of CLVP have concentrated on the solar atmosphere and have not considered the CLVP in cooler non-gray stellar atmospheres at all. In this paper, we present a theoretical study of the CLV of the intensity and the linear polarization in continuum spectra of different spectral type stars. Methods. We solve the radiative transfer equations for polarized light iteratively assuming no magnetic field and considering a plane-parallel model atmospheres and various opacities. Results. We calculate the CLVI and the CLVP for Phoenix stellar model atmospheres for the range of effective temperatures (4500K - 6900K), gravities (log g = 3.0 - 5.0), and wavelengths (4000 - 7000 {\AA}), which are tabulated and available at the Strasbourg astronomical Data Center (CDS). In addition, we present several tests of our code and compare our results with measurements and calculations of CLVI and the CLVP for the Sun. The resulting CLVI are fitted with polynomials and their coefficients are presented in this paper. Conclusions. For the stellar model atmospheres with lower gravity and effective temperature the CLVP is larger.

We present and discuss a polarimetric effect caused by a planet transiting the stellar disk thus breaking the symmetry of the light distribution and resulting in linear polarization of the partially eclipsed star. Estimates of this effect for transiting planets have been made only recently. In particular, we demonstrate that the maximum polarization during transits depends strongly on the centre-to-limb variation of the linear polarization of the host star. However, observational and theoretical studies of the limb polarization have largely concentrated on the Sun. Here we solve the radiative transfer problem for polarized light and calculate the centre-to-limb polarization for one of the brightest transiting planet host HD189733 taking into account various opacities. Using that we simulate the transit effect and estimate the variations of the flux and the linear polarization for HD189733 during the event. As the spots on the stellar disk also break the limb polarization symmetry we simulate the flux and polarization variation due to the spots on the stellar disk.

In this paper, we develop a spectral differential technique with which the dynamical mass of low-mass companions can be found. This method aims at discovering close companions to late-type stars by removing the stellar spectrum through a subtraction of spectra obtained at different orbital phases and discovering the companion spectrum in the difference spectrum in which the companion lines appear twice (positive and negative signal). The resulting radial velocity difference of these two signals provides the true mass of the companion, if the orbital solution for the radial velocities of the primary is known. We select the CO line region in the K-band for our study, because it provides a favourable star-to-companion brightness ratio for our test case GJ1046, an M2V dwarf with a low-mass companion that most likely is a brown dwarf. Furthermore, these lines remain largely unblended in the difference spectrum so that the radial velocity amplitude of the companion can be measured directly. Only if the companion rotates rapidly and has a small radial velocity due to a high mass, does blending occur for all lines so that our approach fails. We also consider activity of the host star, and show that the companion difference flux can be expected to have larger amplitude than the residual signal from the active star so that stellar activity does not inhibit the determination of the companion mass. In addition to determining the companion mass, we restore the single companion spectrum from the difference spectrum using singular value decomposition.

The method of the optical parameter estimations of the nonisothermal giant planet atmospheres by using intensity data of Raman scattering features was used. We applied this method to the observational data of Uranus' geometric albedo spectra from 1981, 1993 and 1995 and obtained the spectral values of the optical depth, namely, $\tau_a/\tau_R$ and $\tau_\kappa/\tau_S$ (where $\tau_a$, $\tau_R$ are aerosol and gas components, $\tau_S=\tau_a+\tau_R$ and $\tau_\kappa$ is absorption component of the effective optical depth of the intensity of diffuse - reflected irradiation forming). We showed that these ratios are different for the three years. The conclusion is that this effect can be due to the horizontal inhomogeneity of aerosol component of optical depth over the Uranus' disk.

A method of the calculation of optical parameters of the nonisothermal giant planet atmospheres was developed using detailed intensity data of Raman scattering. We have used the model of Morozhenko (A.V. Morozhenko, 1997) as a baseline. In such a way, using observational data of Uranus and Neptune (E.Karkoschka, 1994), the spectral values of ratio of optical depth components: aerosol and gas components \tau a/ \tau R, absorbing and scattering components \tau a/ \tau R, and also single scattering albedo of aerosol component corrected for Raman scattering \omega were obtained (where \tau a, \tau R are aerosol and gas components, and \tau ? is absorbing components of effective optical depths of the formation of diffusely reflected irradiation). The averaged value of ratio \tau a/ \tau R is 0.96 but it slowly decreases in the spectral range of 350-450nm for Uranus and \tau a/ \tau R is 1.35 for Neptune.