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.

The Vienna Atomic Line Database (VALD) has been supplemented with new data and new functionality -- the possibility of taking into account the effect of hyperfine splitting (HFS) of atomic levels in the analysis of line profiles. This has been done through the creation of an ancillary SQL database with the HFS constants for atomic levels of 58 isotopes of 30 neutral and singly-ionized atoms. The completeness of the collected data and new opportunities for studies of stars of various spectral types is analyzed. The database enables analysis of splitting of up to 60\%\ of lines with measurable effects in the ultraviolet ($\lambda\gtrsim1000$~\AA), and up to 100\%\ of such lines in the optical and infrared ranges ($\lambda\lesssim25000$~\AA) for A--M stars. In the spectra of hot O--B stars, it is necessary to use laboratory measurements for atoms in the second and higher stages of ionization.

Context. Studying the spots of late-type stars is crucial for distinguishing between the various proposed dynamo mechanisms believed to be the main cause of starspot activity. For this research it is important to collect observation time series that are long enough to unravel both long- and short-term spot evolution. Doppler imaging is a very efficient method for studying spots of stars that cannot be angularly resolved. Aims. High-resolution spectral observations during 1994-2017 are analysed in order to reveal long- and short-term changes in the spot activity of the FK Comae-type subgiant HD 199178. Methods. Most of the observations were collected with the Nordic Optical Telescope. The Doppler imaging temperature maps were calculated using an inversion technique based on Tikhonov regularisation and utilising multiple spectral lines. Results. We present a unique series of 41 temperature maps spanning more than 23 years. All reliable images show a large cool spot region centred near the visible rotation pole. Some lower latitude cool features are also recovered, although the reliability of these is questionable. There is an expected anti-correlation between the mean surface temperature and the spot coverage. Using the Doppler images, we construct the equivalent of a solar butterfly diagram for HD 199178. Conclusions. HD 199178 clearly has a long-term large and cool spot structure at the rotational pole. This spot structure dominated the spot activity during the years 1994-2017. The size and position of the structure has evolved with time, with a gradual increase during the last years. The lack of lower latitude features prevents the determination of a possible differential rotation.

Context. LQ Hya is one of the most frequently studied young solar analogue stars. Recently, it has been observed to show intriguing behaviour by analysing long-term photometry: during 2003--2009, a coherent spot structure migrating in the rotational frame has been reported by various authors, but since that time the star has entered a chaotic state where coherent structures seem to have disappeared and rapid phase jumps of the photometric minima occur irregularly over time. Aims. LQ Hya is one of the stars included in the SOFIN/FIES long-term monitoring campaign extending over 25 years. Here we publish new temperature maps for the star during 2006--2017, covering the chaotic state of the star. Methods. We use a Doppler imaging technique to derive surface temperature maps from high-resolution spectra. Results. From the mean temperatures of the Doppler maps we see a weak but systematic increase in the surface temperature of the star. This is consistent with the simultaneously increasing photometric magnitude. During nearly all observing seasons we see a high-latitude spot structure which is clearly nonaxisymmetric. The phase behaviour of this structure is very chaotic but agrees reasonably well with the photometry. Equatorial spots are also frequently seen, but many of them we interpret to be artefacts due to the poor to moderate phase coverage. Conclusions. Even during the chaotic phase of the star, the spot topology has remained very similar to the higher activity epochs with more coherent and long-lived spot structures: we see high-latitude and equatorial spot activity, the mid latitude range still being most often void of spots. We interpret the erratic jumps and drifts in phase of the photometric minima to be caused by changes in the high-latitude spot structure rather than the equatorial spots.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

Starspots are important manifestations of stellar magnetic activity. By studying their behaviour in young solar analogues, we can unravel the properties of their magnetic cycles. This gives crucial information of the underlying dynamo process. Comparisons with the solar cycle enable us to infer knowledge about how the solar dynamo has evolved during the Sun's lifetime. Here we study the correlation between photometric brightness variations, spottedness and mean temperature in V889 Her, a young solar analogue. Our data covers 18 years of spectroscopic and 25 years of photometric observations. We use Doppler imaging to derive temperature maps from high-resolution spectra. We use the Continuous Period Search method to retrieve mean V-magnitudes from photometric data. Our Doppler imaging maps show a persistent polar spot structure varying in strength. This structure is centered slightly off the rotational pole. The mean temperature derived from the maps shows an overall decreasing trend, as does the photometric mean brightness, until it reaches its minimum around 2017. The filling factor of cool spots, however, shows only a weak tendency to anti-correlate with the decreasing mean brightness. We interpret V889 Her to have entered into a grand maximum in its activity. The clear relation between the mean temperature of the Doppler imaging surface maps and the mean magnitude supports the reliability of the Doppler images. The lack of correlation between the mean magnitude and the spottedness may indicate that bright features in the Doppler images are real.

Access to molecular and atomic data is critical for many forms of remote sensing analysis across different fields. Many atomic and molecular databases are however highly specialized for their intended application, complicating querying and combination data between sources. The Virtual Atomic and Molecular Data Centre, VAMDC, is an electronic infrastructure that allows each database to register as a "node". Through services such as VAMDC's portal website, users can then access and query all nodes in a homogenized way. Today all major Atomic and Molecular databases are attached to VAMDC. This article describes the software tools we developed to help data providers create and manage a VAMDC node. It gives an overview of the VAMDC infrastructure and of the various standards it uses. The article then discusses the development choices made and how the standards are implemented in practice. It concludes with a full example of implementing a VAMDC node using a real-life case as well as future plans for the node software.

We propose a new data analysis method for obtaining transmission spectra of exoplanet atmospheres and brightness variation across the stellar disk from transit observations. The new method is capable of recovering exoplanet atmosphere absorption spectra and stellar specific intensities without relying on theoretical models of stars and planets. We fit both stellar specific intensity and planetary radius simultaneously directly to transit light curves. This allows stellar models to be removed from the data analysis. Furthermore, we use a data-quality weighted filtering technique to achieve an optimal trade-off between spectral resolution and reconstruction fidelity homogenising the signal to noise ratio across the wavelength range. Such approach is more efficient than conventional data binning onto a low resolution wavelength grid. We demonstrate that our analysis is capable of re-producing results achieved by using explicit quadratic limb darkening equation, and that the filtering technique helps eliminating spurious spectral features in regions with strong telluric absorption. The method is applied to the VLT FORS2 observations of the exoplanets GJ 1214 b and WASP-49 b, and our results are in agreement with previous studies. Comparisons between obtained stellar specific intensity and numerical models indicates that the method is capable of accurately reconstructing the specific intensity. The proposed method enables more robust characterization of exoplanetary atmospheres by separating derivation of planetary transmission and stellar specific intensity spectra (that is model-independent) from chemical and physical interpretation.

The program package SME (Spectroscopy Made Easy), designed to perform an analysis of stellar spectra using spectral fitting techniques, was updated due to adding new functions (isotopic and hyperfine splittins) in VALD and including grids of NLTE calculations for energy levels of few chemical elements. SME allows to derive automatically stellar atmospheric parameters: effective temperature, surface gravity, chemical abundances, radial and rotational velocities, turbulent velocities, taking into account all the effects defining spectral line formation. SME package uses the best grids of stellar atmospheres that allows us to perform spectral analysis with the similar accuracy in wide range of stellar parameters and metallicities - from dwarfs to giants of BAFGK spectral classes.

Today Vienna Atomic Line Database (VALD) is one of main databases of atomic and molecular parameters required for stellar spectra analysis. We present the new features that recently appeared in the VALD3 release, including the effects of isotopic composition and hyperfine splitting. The latest version of VALD contains parameters for several isotopes of Li, Ca, Ti, Cu, Ba, Eu, and hyperfine splitting of 35 isotopes from Li to Eu.

The first generation of E-ELT instruments will include an optical-infrared High Resolution Spectrograph, conventionally indicated as EELT-HIRES, which will be capable of providing unique breakthroughs in the fields of exoplanets, star and planet formation, physics and evolution of stars and galaxies, cosmology and fundamental physics. A 2-year long phase A study for EELT-HIRES has just started and will be performed by a consortium composed of institutes and organisations from Brazil, Chile, Denmark, France, Germany, Italy, Poland, Portugal, Spain, Sweden, Switzerland and United Kingdom. In this paper we describe the science goals and the preliminary technical concept for EELT-HIRES which will be developed during the phase A, as well as its planned development and consortium organisation during the study.

We present a catalog of uniformly determined stellar properties and abundances for 1626 F, G, and K stars using an automated spectral synthesis modeling procedure. All stars were observed using the HIRES spectrograph at Keck Observatory. Our procedure used a single line list to fit model spectra to observations of all stars to determine effective temperature, surface gravity, metallicity, projected rotational velocity, and the abundances of 15 elements (C, N, O, Na, Mg, Al, Si, Ca, Ti, V, Cr, Mn, Fe, Ni, & Y). Sixty percent of the sample had Hipparcos parallaxes and V-band photometry which we combined with the spectroscopic results to obtain mass, radius, and luminosity. Additionally, we used the luminosity, effective temperature, metallicity and alpha-element enhancement to interpolate in the Yonsei-Yale isochrones to derive mass, radius, gravity, and age ranges for those stars. Finally, we determined new relations between effective temperature and macroturbulence for dwarfs and subgiants. Our analysis achieved precisions of 25 K in Teff , 0.01 dex in [M/H], 0.028 dex for log g and 0.5 km/s in v sin ibased on multiple observations of the same stars. The abundance results were similarly precise, between 0.01 and - 0.04 dex, though trends with respect to Teff remained for which we derived empirical corrections. The trends, though small, were much larger than our uncertainties and are shared with published abundances. We show that changing our model atmosphere grid accounts for most of the trend in [M/H] between 5000 K and 5500 K indicating a possible problem with the atmosphere models or opacities.

Context. The Spectroscopy Made Easy (SME) package has become a popular tool for analyzing stellar spectra, often in connection with large surveys or exoplanet research. SME has evolved significantly since it was first described in 1996, but many of the original caveats and potholes still haunt users. The main drivers for this paper are complexity of the modeling task, the large user community, and the massive effort that has gone into SME. Aims. We do not intend to give a comprehensive introduction to stellar atmospheres, but will describe changes to key components of SME: the equation of state, opacities, and radiative transfer. We will describe the analysis and fitting procedure and investigate various error sources that affect inferred parameters. Methods. We review the current status of SME, emphasizing new algorithms and methods. We describe some best practices for using the package, based on lessons learned over two decades of SME usage. We present a new way to assess uncertainties in derived stellar parameters. Results. Improvements made to SME, better line data, and new model atmospheres yield more realistic stellar spectra, but in many cases systematic errors still dominate over measurement uncertainty. Future enhancements are outlined.

We present a description of a new instrument development, HARPS3, planned to be installed on an upgraded and roboticized Isaac Newton Telescope by end-2018. HARPS3 will be a high resolution (R = 115,000) echelle spectrograph with a wavelength range from 380-690 nm. It is being built as part of the Terra Hunting Experiment - a future 10 year radial velocity measurement programme to discover Earth-like exoplanets. The instrument design is based on the successful HARPS spectrograph on the 3.6m ESO telescope and HARPS-N on the TNG telescope. The main changes to the design in HARPS3 will be: a customised fibre adapter at the Cassegrain focus providing a stabilised beam feed and on-sky fibre diameter ~ 1.4 arcsec, the implementation of a new continuous flow cryostat to keep the CCD temperature very stable, detailed characterisation of the HARPS3 CCD to map the effective pixel positions and thus provide an improved accuracy wavelength solution, an optimised integrated polarimeter and the instrument integrated into a robotic operation. The robotic operation will optimise our programme which requires our target stars to be measured on a nightly basis. We present an overview of the entire project, including a description of our anticipated robotic operation.

We performed extensive tests of the accuracy of atmospheric parameter determination for FGK stars based on the spectrum fitting procedure Spectroscopy Made Easy (SME). Our stellar sample consists of 13 objects, including the Sun, in the temperature range 5000--6600~K and metallicity range -1.4 -- +0.4. The analysed stars have the advantage of having parameters derived by interferometry. For each star we use spectra obtained with different spectrographs and different signal-to-noise ratios (S/N). For the fitting we adopted three different sets of constraints and test how the derived parameters depend upon the spectral regions (masks) used in SME. We developed and implemented in SME a new method for estimating uncertainties in the resulting parameters based on fitting residuals, partial derivatives, and data uncertainties. For stars in the 5700--6600 K range the best agreement with the effective temperatures derived by interferometry is achieved when spectrum fitting includes the H$\alpha$ and H$\beta$ lines, while for cooler stars the choice of the mask does not affect the results. The derived atmospheric parameters do not strongly depend on spectral resolution and S/N of the observations, while the uncertainties in temperature and surface gravity increase with increasing effective temperature, with minima at 50~K in Teff and 0.1~dex in log g, for spectra with S/N=150--200. A NLTE analysis of the TiI/TiII and FeI/FeII ionisation equilibria and abundances determined from the atomic CI (NLTE) and molecular CH species supports the parameters we derived with SME by fitting the observed spectra including the hydrogen lines.

We propose a method for observing transiting exoplanets with near-infrared high-resolution spectrometers. We aim to create a robust data analysis method for recovering atmospheric transmission spectra from transiting exoplanets over a wide wavelength range in the near infrared. By using an inverse method approach, combined with stellar models and telluric transmission spectra, the method recovers the transiting exoplanet's atmospheric transmittance at high precision over a wide wavelength range. We describe our method and have tested it by simulating observations. This method is capable of recovering transmission spectra of high enough accuracy to identify absorption features from molecules such as O2, CH4, CO2, and H2O. This accuracy is achievable for Jupiter-size exoplanetsat S/N that can be reached for 8m class telescopes using high-resolution spectrometers (R>20 000) during a single transit, and for Earth-size planets and super-Earths transiting late K or M dwarf stars at S/N reachable during observations of less than 10 transits. We also analyse potential error sources to show the robustness of the method. Detection and characterization of atmospheres of both Jupiter-size planets and smaller rocky planets looks promising using this set-up.

We study the spot distribution on the surface of LQ~Hya during the observing seasons October 1998 -- November 2002. We look for persistent active longitudes, trends in the level of spot activity and compare to photometric data, specifically to the derived time epochs of the lightcurve minima. We apply the Doppler Imaging technique on photospheric spectral lines using an inversion code to retrieve images of the surface temperature. We present new temperature maps using multiple spectral lines for a total of 7 seasons. We calculate spot coverage fractions from each map, and as a result we find a general trend that is interpreted to be an indication of a spot cycle. There is a minimum during the observing season of March 1999. After this the activity increases until November 2000, followed by a general decrease in activity again. We find no evidence for active longitudes persisting over multiple observing seasons. The spot activity appears to be concentrated to two latitude regions. The high latitude spots are particularly strong when the spot coverage is at a maximum. Using the currently accepted rotation period, we find spot structures to show a trend in the phase-time plot, indicative of a need for a longer period. We conclude that the long-term activity of LQ~Hya is more chaotic than that of some magnetically active binary stars analyzed with similar methods, but still with clear indications of an activity cycle.

We demonstrate a new procedure to derive accurate and precise surface gravities from high resolution spectra without the use of external constraints. Our analysis utilizes Spectroscopy Made Easy (SME) with robust spectral line constraints and uses an iterative process to mitigate degeneracies in the fitting process. We adopt an updated radiative transfer code, a new treatment for neutral perturber broadening, a line list with multiple gravity constraints and separate fitting for global stellar properties and abundance determinations. To investigate the sources of temperature dependent trends in determining log g noted in previous studies, we obtained Keck HIRES spectra of 42 Kepler asteroseismic stars. In comparison to asteroseismically determined log g our spectroscopic analysis has a constant offset of 0.01 dex with a root mean square (RMS) scatter of 0.05 dex. We also analyzed 30 spectra which had published surface gravities determined using the $a/R_*$ technique from planetary transits and found a constant offset of 0.06 dex and RMS scatter of 0.07 dex. The two samples covered effective temperatures between 5000K and 6700K with log g between 3.7 and 4.6.

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.

The detection of exoplanets using any method is prone to confusion due to the intrinsic variability of the host star. We investigate the effect of cool starspots on the detectability of the exoplanets around solar-like stars using the radial velocity method. For investigating this activity-caused "jitter" we calculate synthetic spectra using radiative transfer, known stellar atomic and molecular lines, different surface spot configurations, and an added planetary signal. Here, the methods are described in detail, tested and compared to previously published studies. The methods are also applied to investigate the activity jitter in old and young solar-like stars, and over a solar-like activity cycles. We find that the mean full jitter amplitude obtained from the spot surfaces mimicking the solar activity varies during the cycle approximately between 1 m/s and 9 m/s. With a realistic observing frequency a Neptune mass planet on a one year orbit can be reliably recovered. On the other hand, the recovery of an Earth mass planet on a similar orbit is not feasible with high significance. The methods developed in this study have a great potential for doing statistical studies of planet detectability, and also for investigating the effect of stellar activity on recovered planetary parameters.

Exoplanet science is booming. In 20 years our knowledge has expanded considerably, from the first discovery of a Hot Jupiter, to the detection of a large population of Neptunes and super-Earths, to the first steps toward the characterization of exoplanet atmospheres. Between today and 2025, the field will evolve at an even faster pace with the advent of several space-based transit search missions, ground-based spectrographs, high-contrast imaging facilities, and the James Webb Space Telescope. Especially the ESA M-class PLATO mission will be a game changer in the field. From 2024 onwards, PLATO will find transiting terrestrial planets orbiting within the habitable zones of nearby, bright stars. These objects will require the power of Extremely Large Telescopes (ELTs) to be characterized further. The technique of ground-based high-resolution spectroscopy is establishing itself as a crucial pathway to measure chemical composition, atmospheric structure and atmospheric circulation in transiting exoplanets. A high-resolution spectrograph covering the visible and near-IR domains, mounted on the European ELT, will be able to detect molecules such as water vapour, carbon dioxide and oxygen in the atmospheres of habitable planets under favourable circumstances. E-ELT HiRES is the perfect ground-based match to the PLATO space mission and represents a unique opportunity for Europe to lead the world into the era of exploration of exoplanets with habitable conditions. HiRES will also be extremely complementary to other E-ELT planned instruments specialising in different kinds of planets, such as METIS and EPICS.

Intermediate-mass, magnetic chemically peculiar (Ap) stars provide a unique opportunity to study the topology of stellar magnetic fields in detail and to investigate magnetically driven processes of spot formation. Here we aim to derive the surface magnetic field geometry and chemical abundance distributions for the extraordinary Ap star HD 75049. This object hosts a surface field of ~30 kG, one of the strongest known for any non-degenerate star. We used time-series of high-resolution HARPS intensity and circular polarisation observations. These data were interpreted with the help of magnetic Doppler imaging and model atmospheres incorporating effects of a non-solar chemical composition and a strong magnetic field. We refined the rotational period of HD 75049 to Prot=4.048267+/-0.000036 d. We also derived basic stellar parameters, Teff=10250+/-250 K and logg=4.3+/-0.1. Magnetic Doppler imaging revealed that the field topology of HD 75049 is poloidal and dominated by a dipolar contribution with a peak surface field strength of 39 kG. At the same time, deviations from the classical axisymmetric oblique dipolar configuration are significant. Chemical surface maps of Si, Cr, Fe, and Nd show abundance contrasts of 0.5-1.4 dex, which is low compared with many other Ap stars. Nd is found to be enhanced close to the magnetic pole, whereas Si and Cr are concentrated predominantly at the magnetic equator. The iron distribution shows low-contrast features both at the magnetic equator and the pole. The morphology of the magnetic field and the properties of chemical spots in HD 75049 are qualitatively similar to those of Ap stars with weaker fields. Consequently, whatever mechanism forms and sustains global magnetic fields in intermediate-mass main-sequence stars, it operates in the same way over the entire observed range of magnetic field strengths.

Aims: We present a magnetic Doppler imaging study from all Stokes parameters of the cool, chemically peculiar star HD 24712. This is the very first such analysis performed at a resolving power exceeding 10^5. Methods: The analysis is performed on the basis of phase-resolved observations of line profiles in all four Stokes parameters obtained with the HARPSpol instrument attached at the 3.6-m ESO telescope. We use the magnetic Doppler imaging code, INVERS10, which allows us to derive the magnetic field geometry and surface chemical abundance distributions simultaneously. Results: We report magnetic maps of HD 24712 recovered from a selection of FeI, FeII, NdIII, and NaI lines with strong polarization signals in all Stokes parameters. Our magnetic maps successfully reproduce most of the details available from our observation data. We used these magnetic field maps to produce abundance distribution map of Ca. This new analysis shows that the surface magnetic field of HD 24712 has a dominant dipolar component with a weak contribution from higher-order harmonics. The surface abundance distributions of Fe and Ca show enhancements near the magnetic equator with an underabundant patch at the visible (positive) magnetic pole; Nd is highly abundant around the positive magnetic pole. The Na abundance map shows a high overabundance around the negative magnetic pole. Conclusions: Based on our investigation and similar recent four Stokes parameter magnetic mapping studies, we present a tentative evidence for the hypothesis that Ap stars with dipole-like fields are older than stars with magnetic fields that have more small-scale structures. We find that our abundance maps are inconsistent with recent theoretical calculations of atomic diffusion in presence of magnetic fields.

Magnetic Doppler imaging from four Stokes parameter observations has uncovered a new level of complexity of stellar magnetic fields previously not known. This new information is of interest to the evolution and structure of magnetic fields of intermediate and high-mass stars. Recently, we performed magnetic Doppler imaging study of the chemically peculiar star HD 24712 using phase-resolved spectropolarimetric observations of line-profiles in all four Stokes parameters. The observations were obtained with the HARPSpol instrument on the 3.6-m ESO telescope. Our analysis shows that the magnetic field of HD 24712 has very strong dipolar component with weak small-scale contributions. The finding of dipole-like field for HD 24712 gives evidence for the hypothesis that old Ap stars have predominantly dipolar magnetic fields. We plan on performing magnetic Doppler imaging from four Stokes parameter observations for other Ap/Bp stars as well, which will help us understand the nature of magnetic fields in these stars.

CRIRES, the ESO high resolution infrared spectrometer, is a unique instrument which allows astronomers to access a parameter space which up to now was largely uncharted. In its current setup, it consists of a single-order spectrograph providing long-slit, single-order spectroscopy with resolving power up to R=100,000 over a quite narrow spectral range. This has resulted in sub-optimal efficiency and use of telescope time for all the scientific programs requiring broad spectral coverage of compact objects (e.g. chemical abundances of stars and intergalactic medium, search and characterization of extra-solar planets). To overcome these limitations, a consortium was set-up for upgrading CRIRES to a cross-dispersed spectrometer, called CRIRES+. This paper presents the updated optical design of the crossdispersion module for CRIRES+. This new module can be mounted in place of the current pre-disperser unit. The new system yields a factor of >10 increase in simultaneous spectral coverage and maintains a quite long slit (10"), ideal for observations of extended sources and for precise sky-background subtraction.

Based on detailed spectral synthesis we carry out quantitative measurements of the strength and complexity of surface magnetic fields in the four well-known M-dwarfs GJ 388, GJ 729, GJ 285, and GJ 406 populating the mass regime around the boundary between partially and fully convective stars. Very high resolution R=100000, high signal-to-noise (up to 400) near-infrared Stokes I spectra were obtained with CRIRES at ESO's Very Large Telescope covering regions of the FeH Wing-Ford transitions at 1mum. The field distributions in all four stars are characterized by three distinct groups of field components, the data are neither consistent with a smooth distribution of different field strengths, nor with one average field strength covering the full star. We find evidence of a subtle difference in the field distribution of GJ 285 compared to the other three targets. GJ 285 also has the highest average field of 3.5kG and the strongest maximum field component of 7-7.5kG. The maximum local field strengths in our sample seem to be correlated with rotation rate. While the average field strength is saturated, the maximum local field strengths in our sample show no evidence for saturation. We find no difference between the field distributions of partially and fully convective stars. The one star with evidence for a field distribution different to the other three is the most active star (i.e. with largest x-ray luminosity and mean surface magnetic field) rotating relatively fast. A possible explanation is that rotation determines the distribution of surface magnetic fields, and that local field strengths grow with rotation even in stars in which the average field is already saturated.

Current knowledge about stellar magnetic fields relies almost entirely on circular polarization observations. Few objects have been observed in all four Stokes parameters. The magnetic Ap star HD 24712 (DO Eri, HR 1217) was recently observed in the Stokes IQUV parameters with the HARPSpol instrument at the 3.6-m ESO telescope as part of our project at investigating Ap stars in all four Stokes parameters. The resulting spectra have dense phase coverage, resolving power > 10^5, and S/N ratio of 300-600. These are the highest quality full Stokes observations obtained for any star other than the Sun. We present preliminary results from magnetic Doppler imaging of HD 24712. This analysis is the first step towards obtaining detailed 3-D maps of magnetic fields and abundance structures for HD 24712 and other Ap stars that we currently observe with HARPSpol.

Our current knowledge about the stellar magnetic fields relies almost entirely on the circular polarization observations. Very few objects have been observed so far in all four Stokes parameters. The magnetic Ap star HD24712 (DO Eri, HR1217) was recently observed in the Stokes IQUV parameters with the HARPSpol instrument at the 3.6-m ESO telescope. The resulting spectra have dense rotational phase coverage, resolving power > 10^5, and a S/N ratio of 300-600. These are the highest quality full Stokes observations obtained for any star other than the Sun. Here we present preliminary results of the magnetic Doppler imaging of HD24712 using HARPSpol spectra. This analysis is the first step towards obtaining detailed 3-D maps of magnetic fields abundance structures for HD24712 and in other Ap stars that we currently observe with HARPSpol.

High-resolution spectropolarimetric observations provide simultaneous information about stellar magnetic field topologies and three-dimensional distributions of chemical elements. Here we present analysis of a unique full Stokes vector spectropolarimetric data set, acquired for the cool magnetic Ap star HD 24712. The goal of our work is to examine circular and linear polarization signatures inside spectral lines and to study variation of the stellar spectrum and magnetic observables as a function of rotational phase. HD 24712 was observed with the HARPSpol instrument at the 3.6-m ESO telescope over a period of 2010-2011. The resulting spectra have S/N ratio of 300-600 and resolving power exceeding 100000. The multiline technique of least-squares deconvolution (LSD) was applied to combine information from the spectral lines of Fe-peak and rare-earth elements. We used the HARPSPol spectra of HD 24712 to study the morphology of the Stokes profile shapes in individual spectral lines and in LSD Stokes profiles corresponding to different line masks. From the LSD Stokes V profiles we measured the longitudinal component of the magnetic field, <Bz>, with an accuracy of 5-10 G. We also determined the net linear polarization from the LSD Stokes Q and U profiles. We determined an improved rotational period of the star, P_rot = 12.45812 +/- 0.00019d. We measured <Bz> from the cores of Halpha and Hbeta lines. The analysis of <Bz> measurements showed no evidence for a significant radial magnetic field gradient in the atmosphere of HD 24712. We used our <Bz> and net linear polarization measurements to determine parameters of the dipolar magnetic field topology. We found that magnetic observables can be reasonably well reproduced by the dipolar model. We discovered rotational modulation of the Halpha core and related it a non-uniform surface distribution of rare-earth elements.

Several recent spectrophotometric studies failed to detect significant global magnetic fields in late-B HgMn chemically peculiar stars, but some investigations have suggested the presence of strong unstructured or tangled fields in these objects. We used detailed spectrum synthesis analysis to search for evidence of tangled magnetic fields in high-quality observed spectra of 8 slowly rotating HgMn stars and one normal late-B star. We also evaluated recent sporadic detections of weak longitudinal magnetic fields in HgMn stars based on the moment technique. Our analysis of the Zeeman broadening of magnetically sensitive spectral lines reveals no evidence of tangled magnetic fields in any of the studied HgMn or normal stars. We infer upper limits of 200-700 G for the mean magnetic field modulus -- much smaller than the field strengths implied by studies based on differential magnetic line intensification and quadratic field diagnostics. The new HARPSpol longitudinal field measurements for the extreme HgMn star HD 65949 and the normal late-B star 21 Peg are consistent with zero at a precision of 3-6 G. Re-analysis of our Stokes V spectra of the spotted HgMn star HD 11753 shows that the recent moment technique measurements retrieved from the same data are incompatible with the lack of circular polarization signatures in the spectrum of this star. We conclude that there is no evidence for substantial tangled magnetic fields on the surfaces of studied HgMn stars. We cannot independently confirm the presence of very strong quadratic or marginal longitudinal fields for these stars, so results from the moment technique are likely to be spurious.

In atmospheres of magnetic main-sequence stars, the diffusion of chemical elements leads to a number of observed anomalies, such as abundance spots across the stellar surface. The aim of this study was to derive a detailed picture of the surface abundance distribution of the magnetic chemically peculiar star HD 3980. Based on high-resolution, phase-resolved spectroscopic observations of the magnetic A-type star HD 3980, the inhomogeneous surface distribution of 13 chemical elements (Li, O, Si, Ca, Cr, Mn, Fe, La, Ce, Pr, Nd, Eu, and Gd) has been reconstructed. The INVERS12 code was used to invert the rotational variability in line profiles to elemental surface distributions. Assuming a centered, dominantly dipolar magnetic field configuration, we find that Li, O, Mg, Pr, and Nd are mainly concentrated in the area of the magnetic poles and depleted in the regions around the magnetic equator. The high abundance spots of Si, La, Ce, Eu, and Gd are located between the magnetic poles and the magnetic equator. Except for La, which is clearly depleted in the area of the magnetic poles, no obvious correlation with the magnetic field has been found for these elements otherwise. Ca, Cr, and Fe appear enhanced along the rotational equator and the area around the magnetic poles. The intersection between the magnetic and the rotational equator constitutes an exception, especially for Ca and Cr, which are depleted in that region. No obvious correlation between the theoretically predicted abundance patterns and those determined in this study could be found. This can be attributed to a lack of up-to-date theoretical models, especially for rare earth elements.

We present high spectral resolution ($R\approx108,000$) Stokes $V$ polarimetry of the Classical T Tauri stars (CTTSs) GQ Lup and TW Hya obtained with the polarimetric upgrade to the HARPS spectrometer on the ESO 3.6 m telescope. We present data on both photospheric lines and emission lines, concentrating our discussion on the polarization properties of the \ion{He}{1} emission lines at 5876 \AA\ and 6678 \AA. The \ion{He}{1} lines in these CTTSs contain both narrow emission cores, believed to come from near the accretion shock region on these stars, and broad emission components which may come from either a wind or the large scale magnetospheric accretion flow. We detect strong polarization in the narrow component of the two \ion{He}{1} emission lines in both stars. We observe a maximum implied field strength of $6.05 \pm 0.24$ kG in the 5876 \AA\ line of GQ Lup, making it the star with the highest field strength measured in this line for a CTTS. We find field strengths in the two \ion{He}{1} lines that are consistent with each other, in contrast to what has been reported in the literature on at least one star. We do not detect any polarization in the broad component of the \ion{He}{1} lines on these stars, strengthening the conclusion that they form over a substantially different volume relative the formation region of the narrow component of the \ion{He}{1} lines.

Stellar activity signatures such as spots and plage can significantly limit the search for extrasolar planets. Current models of activity-induced radial velocity (RV) signals focused on the impact of temperature contrast in spots predicting the signal to diminish toward longer wavelengths. On the other hand, the relative importance of the Zeeman effect on RV measurements should grow with wavelength, because the Zeeman displacement itself grows with \lambda, and because a magnetic and cool spot contributes more to the total flux at longer wavelengths. We model the impact of active regions on stellar RV measurements including both temperature contrast in spots and Zeeman line broadening. We calculate stellar line profiles using polarized radiative transfer models including atomic and molecular Zeeman splitting from 0.5 to 2.3\mum. Our results show that the amplitude of the RV signal caused by the Zeeman effect alone can be comparable to that caused by temperature contrast. Furthermore, the RV signal caused by cool and magnetic spots increases with wavelength. We also calculate the RV signal due to variations in average magnetic field strength from one observation to the next, but find it unlikely that this can significantly influence the search for extrasolar planets. As an example, we derive the RV amplitude of the active M dwarf AD Leo as a function of wavelength using data from the HARPS spectrograph. The RV signal does not diminish at longer wavelengths but shows evidence for the opposite behavior. We conclude that the RV signal of active stars does not vanish at longer wavelength but sensitively depends on the combination of spot temperature and magnetic field; in active low-mass stars, it is even likely to grow with wavelength.

We present two new accurate and efficient method to compute the formal solution of the polarized radiative transfer equation. In this work, the source function and the absorption matrix are approximated using quadratic and cubic Bezier spline interpolants. These schemes provide 2nd and 3rd order approximation respectively and don't suffer from erratic behavior of the polynomial approximation (overshooting). The accuracy and the convergence of the new method are studied along with other popular solutions of the radiative transfer equation, using stellar atmospheres with strong gradients in the line-of-sight velocity and in the magnetic-field vector.

Mercury-manganese (HgMn) stars have been considered as non-magnetic and non-variable chemically peculiar (CP) stars for a long time. However, recent discoveries of the variability in spectral line profiles have suggested an inhomogeneous surface distribution of chemical elements in some HgMn stars. From the studies of other CP stars it is known that magnetic field plays a key role in the formation of surface spots. All attempts to find magnetic fields in HgMn stars have yielded negative results. In this study, we investigate the possible presence of a magnetic field in phi Phe (HD 11753) and reconstruct surface distribution of chemical elements that show variability in spectral lines.We also test a hypothesis that a magnetic field is concentrated in chemical spots and look into the possibility that some chemical elements are stratified with depth in the stellar atmosphere. Our analysis is based on high-quality spectropolarimetric time-series observations, covering a full rotational period of the star. Spectra were obtained with the HARPSpol at the ESO 3.6-m telescope. Combining information from all suitable spectral lines, we set an upper limit of 4 G on the mean longitudinal magnetic field. For chemical spots, an upper limit on the longitudinal field varies between 8 and 15 G. We confirmed the variability of Y, Sr, and Ti and detected variability in Cr lines. Stratification analysis showed that Y and Ti are not concentrated in the uppermost atmospheric layers. Our spectropolarimetric observations rule out the presence of a strong, globally-organised magnetic field in phi Phe. This implies an alternative mechanism of spot formation, which could be related to a non-equilibrium atomic diffusion. However, the typical time scales of the variation in stratification predicted by the recent time-dependent diffusion models exceed significantly the spot evolution time-scale reported for phi Phe.

Aims. We study the spot activity of II Peg during the years 2004-2010 to determine long- and short-term changes in the magnetic activity. In a previous study, we detected a persistent active longitude, as well as major changes in the spot configuration occurring on a timescale of shorter than a year. The main objective of this study is to determine whether the same phenomena persist in the star during these six years of spectroscopic monitoring. Methods. The observations were collected with the high-resolution SOFIN spectrograph at the Nordic Optical Telescope. The temperature maps were calculated using a Doppler imaging code based on Tikhonov regularization. Results. We present 12 new temperature maps that show spots distributed mainly over high and intermediate latitudes. In each image, 1-3 main active regions can be identified. The activity level of the star is clearly lower than during our previous study for the years 1994-2002. In contrast to the previous observations, we detect no clear drift of the active regions with respect to the rotation of the star. Conclusions. Having shown a systematic longitudinal drift of the spot-generating mechanism during 1994-2002, the star has clearly switched to a low-activity state for 2004-2010, during which the spot locations appear more random over phase space. It could be that the star is near to a minimum of its activity cycle.

The relativley large spread in the derived metallicities ([Fe/H]) of M dwarfs shows that various approaches have not yet converged to consistency. The presence of strong molecular features, and incomplete line lists for the corresponding molecules have made metallicity determinations of M dwarfs difficult. Furthermore, the faint M dwarfs require long exposure times for a signal-to-noise ratio sufficient for a detailed spectroscopic abundance analysis. We present a high-resolution (R~50,000) spectroscopic study of a sample of eight single M dwarfs and three wide-binary systems observed in the infrared J-band. The absence of large molecular contributions allow for a precise continuum placement. We derive metallicities based on the best fit synthetic spectra to the observed spectra. To verify the accuracy of the applied atmospheric models and test our synthetic spectrum approach, three binary systems with a K-dwarf primary and an M-dwarf companion were observed and analysed along with the single M dwarfs. We obtain a good agreement between the metallicities derived for the primaries and secondaries of our test binaries and thereby confirm the reliability of our method of analysing M dwarfs. Our metallicities agree well with certain earlier determinations, and deviate from others. We conclude that spectroscopic abundance analysis in the J band is a reliable method for establishing the metallicity scale for M dwarfs. We recommend its application to a larger sample covering lower as well as higher metallicities. Further prospects of the method include abundance determinations for individual elements.

Chemically peculiar stars of the mercury-manganese (HgMn) type represent a new class of spotted late-B stars, in which evolving surface chemical inhomogeneities are apparently unrelated to the presence of strong magnetic fields but are produced by some hitherto unknown astrophysical mechanism. The goal of this study is to perform a detailed line profile variability analysis and carry out a sensitive magnetic field search for one of the brightest HgMn stars - mu Lep. We acquired a set of very high-quality intensity and polarization spectra of mu Lep with the HARPSpol polarimeter. These data were analyzed with the multiline technique of least-squares deconvolution in order to extract information on the magnetic field and line profile variability. Our spectra show very weak but definite variability in the lines of Sc, all Fe-peak elements represented in the spectrum of mu Lep, as well as Y, Sr, and Hg. Variability might also be present in the lines of Si and Mg. Anomalous profile shapes of Ti II and Y II lines suggest a dominant axisymmetric distribution of these elements. At the same time, we found no evidence of the magnetic field in mu Lep, with the 3 sigma upper limit of only 3 G for the mean longitudinal magnetic field. This is the most stringent upper limit on the possible magnetic field derived for a spotted HgMn star. The very weak variability detected for many elements in the spectrum mu Lep suggests that low-contrast chemical inhomogeneities may be common in HgMn stars and that they have not been recognized until now due to the limited precision of previous spectroscopic observations and a lack of time-series data. The null result of the magnetic field search reinforces the conclusion that formation of chemical spots in HgMn stars is not magnetically driven.

Close M-dwarf binaries and higher multiples allow the investigation of rotational evolution and mean magnetic flux unbiased from scatter in inclination angle and age since the orientation of the spin axis of the components is most likely parallel and the individual systems are coeval. Systems composed of an early (M0.0 -- M4.0) and a late (M4.0 -- M8.0) type component offer the possibility to study differences in rotation and magnetism between partially and fully convective stars. We have selected 10 of the closest dM systems to determine the rotation velocities and the mean magnetic field strengths based on spectroscopic analysis of FeH lines of Wing-Ford transitions at 1 $\mu$m observed with VLT/CRIRES. We also studied the quality of our spectroscopic model regarding atmospheric parameters including metallicity. A modified version of the Molecular Zeeman Library (MZL) was used to compute Land\'e g-factors for FeH lines. Magnetic spectral synthesis was performed with the Synmast code. We confirmed previously reported findings that less massive M-dwarfs are braked less effectively than objects of earlier types. Strong surface magnetic fields were detected in primaries of four systems (GJ 852, GJ 234, LP 717-36, GJ 3322), and in the secondary of the triple system GJ 852. We also confirm strong 2 kG magnetic field in the primary of the triple system GJ 2005. No fields could be accurately determined in rapidly rotating stars with $\vsini>10$ \kms. For slow and moderately rotating stars we find the surface magnetic field strength to increase with the rotational velocity $\vsini$ which is consistent with other results from studying field stars.

Galaxies are believed to be the main providers of Lyman continuum (LyC) photons during the early phases of the cosmic reionization. Little is known however, when it comes to escape fractions and the mechanisms behind the leakage. To learn more one may look at local objects, but so far only one low-z galaxy has shown any sign of emitting LyC radiation. With data from the Far Ultraviolet Spectroscopic Explorer (FUSE), Bergvall et al. (2006) found an absolute escape fraction of ionizing photons (f_esc) of 4-10% for the blue compact galaxy Haro 11. However, using a newer version of the reduction pipeline on the same data set, Grimes et al. (2007) could not confirm this and derived an upper limit of f_esc \leq 2%. Here, using the last version of the pipeline CalFUSE v3.2, we aim at settling the question if Haro 11 is emitting ionizing radiation to a significant level or not. We also investigate the performance of the reduction pipeline for faint targets such as Haro 11. At these faint flux levels both FUSE and CalFUSE are pushed to their limits, and a detailed analysis was undertaken. We show that non-simultaneous background estimates are insuffient when working with low signal-to-noise data, and a new background model was developed based on a direct fit to the detector response. Applying the new background fit, a significant signal in the LyC emerged in both detector segments covering these wavelengths. Thus, the leakage is confirmed with a flux density of f_900=4.0 *10^{-15} erg/s/cm^2/{\AA} (S/N=4.6), measured on the airglow free regions in the LyC for the night only data. This corresponds to an absolute escape fraction of ionizing photons for Haro 11 of f_esc=3.3 \pm 0.7%. The results were confirmed by investigating the 2D data, the count rates and the residual flux in CII(1036{\AA}).

The application of high-resolution spectropolarimetry has led to major progress in understanding the magnetism and activity of late-type stars. During the last decade, magnetic fields have been discovered and mapped for many types of active cool stars using spectropolarimetric data. However, these observations and modeling attempts are fundamentally incomplete since they are based on the interpretation of the circular polarization alone. Taking advantage of the newly built HARPS polarimeter, we have obtained the first systematic observations of several cool active stars in all four Stokes parameters. Here we report the detection of magnetically-induced linear polarization for the primary component of the very active RS CVn binary HR 1099 and for the moderately active K dwarf epsilon Eri. For both stars the amplitude of linear polarization signatures is measured to be ~10^{-4} of the unpolarized continuum, which is approximately a factor of ten lower than for circular polarization. This is the first detection of the linear polarization in line profiles of cool active stars. Our observations of the inactive solar-like star alpha Cen A show neither circular nor linear polarization above the level of 10^{-5}, indicating the absence of a net longitudinal magnetic field stronger than 0.2 G.

According to our current understanding, a subclass of the upper main sequence chemically peculiar stars, called mercury-manganese (HgMn), is non-magnetic. Nevertheless, chemical inhomogeneities were recently discovered on their surfaces. At the same time, no global magnetic fields stronger than 1-100 G are detected by modern studies. The goals of our study are to search for magnetic field in the HgMn binary system 66 Eri and to investigate chemical spots on the stellar surfaces of both components. Our analysis is based on high quality spectropolarimetric time-series observations obtained during 10 consecutive nights with the HARPSpol instrument at the ESO 3.6-m telescope. To increase the sensitivity of the magnetic field search we employed a least-squares deconvolution (LSD). We used spectral disentangling to measure radial velocities and study line profile variability. Chemical spot geometry was reconstructed using multi-line Doppler imaging. We report a non-detection of magnetic field in 66 Eri, with error bars 10-24 G for the longitudinal field. Circular polarization profiles also do not indicate any signatures of complex surface magnetic fields. For a simple dipolar field configuration we estimated an upper limit of the polar field strength to be 60-70 G. For the HgMn component we found variability in spectral lines of Ti, Ba, Y, and Sr with the rotational period equal to the orbital one. The surface maps of these elements reconstructed with the Doppler imaging technique, show relative underabundance on the hemisphere facing the secondary component. The contrast of chemical inhomogeneities ranges from 0.4 for Ti to 0.8 for Ba.

We present a database of well determined orbital parameters of exoplanets. This database comprises spectroscopic orbital elements measured for 427 planets orbiting 363 stars from radial velocity and transit measurements as reported in the literature. We have also compiled fundamental transit parameters, stellar parameters, and the method used for the planets discovery. This Exoplanet Orbit Database includes all planets with robust, well measured orbital parameters reported in peer-reviewed articles. The database is available in a searchable, filterable, and sortable form on the Web at http://exoplanets.org through the Exoplanets Data Explorer Table, and the data can be plotted and explored through the Exoplanets Data Explorer Plotter. We use the Data Explorer to generate publication-ready plots giving three examples of the signatures of exoplanet migration and dynamical evolution: We illustrate the character of the apparent correlation between mass and period in exoplanet orbits, the selection different biases between radial velocity and transit surveys, and that the multiplanet systems show a distinct semi-major axis distribution from apparently singleton systems.

Strong surface magnetic fields are ubiquitously found in M-dwarfs with mean intensities on the order of few thousand Gauss-three orders of magnitude higher than the mean surface magnetic field of the Sun. These fields and their interaction with photospheric convection are the main source of stellar activity, which is of big interest to study links between parent stars and their planets. Moreover, the understanding of stellar magnetism, as well as the role of different dynamo-actions in particular, is impossible without explaining magnetic fields in M-dwarfs. Measuring magnetic field intensities and geometries in such cool objects, however, is strongly limited to our ability to simulate the Zeeman effect in molecular lines. In this work, we present quantitative results of modelling and analysis of the magnetic fields in selected M-dwarfs in FeH Wing-Ford lines and strong atomic lines. Some particular FeH lines are found to be the excellent probes of the magnetic field.

We report the discovery of planets orbiting two bright, nearby early K dwarf stars, HD 97658 and Gl 785. These planets were detected by Keplerian modelling of radial velocities measured with Keck-HIRES for the NASA-UC Eta-Earth Survey. HD 97658 b is a close-in super-Earth with minimum mass Msini = 8.2 +/- 1.2 M_Earth, orbital period P = 9.494 +/- 0.005 d, and an orbit that is consistent with circular. Gl 785 b is a Neptune-mass planet with Msini = 21.6 +/- 2.0 M_Earth, P = 74.39 +/- 0.12 d, and orbital eccentricity 0.30 +/- 0.09. Photometric observations with the T12 0.8 m automatic photometric telescope at Fairborn Observatory show that HD 97658 is photometrically constant at the radial velocity period to 0.09 mmag, supporting the existence of the planet.

We performed a highly sensitive search for magnetic fields on a large set of HgMn stars. With the aid of a new polarimeter attached to the HARPS spectrometer at the ESO 3.6m-telescope, we obtained high-quality circular polarization spectra of 41 single and double HgMn stars. Using a multi-line analysis technique on each star, we co-added information from hundreds of spectral lines resulting in significantly greater sensitivity to the presence of magnetic fields, including very weak fields. For the 47 individual objects studied, including 6 components of SB2 systems, we do not detect any magnetic fields at greater than the 3 sigma level. The lack of detection in the circular polarization profiles indicates that if strong fields are present on these stars, they must have complex surface topologies. For simple global fields, our detection limits imply upper limits to the fields present of 2-10 Gauss in the best cases. We conclude that HgMn stars lack large-scale magnetic fields, typical for spotted magnetic Ap stars, sufficient to form and sustain the chemical spots observed on HgMn stars. Our study confirms that in addition to magnetically altered atomic diffusion, there exists another differentiation mechanism operating in the atmospheres of late-B main sequence stars which can compositional inhomogeneities on their surfaces.

We recently commissioned the polarimetric upgrade of the HARPS spectrograph at ESO's 3.6-m telescope at La Silla, Chile. The HARPS polarimeter is capable of full Stokes spectropolarimetry with large sensitivity and accuracy, taking advantage of the large spectral resolution and stability of HARPS. In this paper we present the instrument design and its polarimetric performance. The first HARPSpol observations show that it can attain a polarimetric sensitivity of ~10^-5 (after addition of many lines) and that no significant instrumental polarization effects are present.

Least squares deconvolution (LSD) is a powerful method of extracting high-precision average line profiles from the stellar intensity and polarization spectra. Despite its common usage, the LSD method is poorly documented and has never been tested using realistic synthetic spectra. In this study we revisit the key assumptions of the LSD technique, clarify its numerical implementation, discuss possible improvements and give recommendations how to make LSD results understandable and reproducible. We also address the problem of interpretation of the moments and shapes of the LSD profiles in terms of physical parameters. We have developed an improved, multiprofile version of LSD and have extended the deconvolution procedure to linear polarization analysis taking into account anomalous Zeeman splitting of spectral lines. This code is applied to the theoretical Stokes parameter spectra. We test various methods of interpreting the mean profiles, investigating how coarse approximations of the multiline technique translate into errors of the derived parameters. We find that, generally, the Stokes parameter LSD profiles do not behave as a real spectral line with respect to the variation of magnetic field and elemental abundance. This problem is especially prominent for the Stokes I variation with abundance and Stokes Q variation with magnetic field. At the same time, the Stokes V LSD spectra closely resemble profile of a properly chosen synthetic line for the magnetic field strength up to 1 kG. We conclude that the usual method of interpreting the LSD profiles by assuming that they are equivalent to a real spectral line gives satisfactory results only in a limited parameter range and thus should be applied with caution. A more trustworthy approach is to abandon the single-line approximation of the average profiles and apply LSD consistently to observations and synthetic spectra.

We present first quantitative results of the surface magnetic field measurements in selected M-dwarfs based on detailed spectra synthesis conducted simultaneously in atomic and molecular lines of the FeH Wing-Ford $F^4\,\Delta-X^4\,\Delta$ transitions. A modified version of the Molecular Zeeman Library (MZL) was used to compute Land\'e g-factors for FeH lines in different Hund's cases. Magnetic spectra synthesis was performed with the Synmast code. We show that the implementation of different Hund's case for FeH states depending on their quantum numbers allows us to achieve a good fit to the majority of lines in a sunspot spectrum in an automatic regime. Strong magnetic fields are confirmed via the modelling of atomic and FeH lines for three M-dwarfs YZ~CMi, EV~Lac, and AD~Leo, but their mean intensities are found to be systematically lower than previously reported. A much weaker field ($1.7-2$~kG against $2.7$~kG) is required to fit FeH lines in the spectra of GJ~1224. Our method allows us to measure average magnetic fields in very low-mass stars from polarized radiative transfer. The obtained results indicate that the fields reported in earlier works were probably overestimated by about $15-30$\%. Higher quality observations are needed for more definite results.

We report the discovery of HD 156668b, an extrasolar planet with a minimum mass of M_P sin i = 4.15 M_Earth. This planet was discovered through Keplerian modeling of precise radial velocities from Keck-HIRES and is the second super-Earth to emerge from the NASA-UC Eta-Earth Survey. The best-fit orbit is consistent with circular and has a period of P = 4.6455 d. The Doppler semi-amplitude of this planet, K = 1.89 m/s, is among the lowest ever detected, on par with the detection of GJ 581e using HARPS. A longer period (P ~ 2.3 yr), low-amplitude signal of unknown origin was also detected in the radial velocities and was filtered out of the data while fitting the short-period planet. Additional data are required to determine if the long-period signal is due to a second planet, stellar activity, or another source. Photometric observations using the Automated Photometric Telescopes at Fairborn Observatory show that HD 156668 (an old, quiet K3 dwarf) is photometrically constant over the radial velocity period to 0.1 mmag, supporting the existence of the planet. No transits were detected down to a photometric limit of ~3 mmag, ruling out transiting planets dominated by extremely bloated atmospheres, but not precluding a transiting solid/liquid planet with a modest atmosphere.

We present precise Doppler measurements of four stars obtained during the past decade at Keck Observatory by the California Planet Survey (CPS). These stars, namely, HD 34445, HD 126614, HD 13931, and Gl 179, all show evidence for a single planet in Keplerian motion. We also present Doppler measurements from the Hobby-Eberly Telescope (HET) for two of the stars, HD 34445 and Gl 179, that confirm the Keck detections and significantly refine the orbital parameters. These planets add to the statistical properties of giant planets orbiting near or beyond the ice line, and merit follow-up by astrometry, imaging, and space-borne spectroscopy. Their orbital parameters span wide ranges of planetary minimum mass (M sin i = 0.38-1.9 M_Jup), orbital period (P = 2.87-11.5 yr), semi-major axis (a = 2.1-5.2 AU), and eccentricity (e = 0.02-0.41). HD 34445b (P = 2.87 yr, M sin i = 0.79 M_Jup, e = 0.27) is a massive planet orbiting an old, G-type star. We announce a planet, HD 126614Ab, and an M dwarf, HD 126614B, orbiting the metal-rich star HD 126614 (which we now refer to as HD 126614A). The planet, HD 126614Ab, has minimum mass M sin i = 0.38 M_Jup and orbits the stellar primary with period P = 3.41 yr and orbital separation a = 2.3 AU. The faint M dwarf companion, HD 126614B, is separated from the stellar primary by 489 mas (33 AU) and was discovered with direct observations using adaptive optics and the PHARO camera at Palomar Observatory. The stellar primary in this new system, HD 126614 A, has the highest measured metallicity (\feh = +0.56) of any known planet-bearing star. HD 13931b (P = 11.5 yr, M sin i = 1.88 M_Jup, e = 0.02) is a Jupiter analog orbiting a near solar twin. Gl 179b (P = 6.3 yr, M sin i = 0.82 M_Jup, e = 0.21) is a massive planet orbiting a faint M dwarf. The high metallicity of Gl 179 is consistent with the planet-metallicity correlation among M dwarfs, as documented recently by Johnson and Apps.

We compare surface maps of the chemically peculiar star HD 50773 produced with a Bayesian technique and based on high quality CoRoT photometry with those derived from rotation phase resolved spectropolarimetry. The goal is to investigate the correlation of surface brightness with surface chemical abundance distribution and the stellar magnetic field. The rotational period of the star was determined from a nearly 60 day long continuous light curve obtained during the initial run of CoRoT. Using a Bayesian approach to star-spot modelling, which in this work is applied for the first time for the photometric mapping of a CP star, we derived longitudes, latitudes and radii of four different spot areas. Additional parameters like stellar inclination and the spot's intensities were also determined. The CoRoT observations triggered an extensive ground-based spectroscopic and spectropolarimetric observing campaign and enabled us to obtain 19 different high resolution spectra in Stokes parameters I and V with NARVAL, ESPaDOnS, and SemelPol spectropolarimeters. Doppler and Magnetic Doppler imaging techniques allowed us to derive the magnetic field geometry of the star and the surface abundance distributions of Mg, Si, Ca, Ti, Cr, Fe, Ni, Y, and Cu. We find a dominant dipolar structure of the surface magnetic field. The CoRoT light curve variations and abundances of most elements mapped are correlated with the aforementioned geometry: Cr, Fe, and Si are enhanced around the magnetic poles and coincide with the bright regions on the surface of HD 50773 as predicted by our light curve synthesis and confirmed by photometric imaging.

The active RS CVn star II Peg has been spectroscopically monitored for almost 18 years with the SOFIN spectrograph at NOT, La Palma, Spain. In this paper we present five new surface temperature maps of the object for the years 1999 (two maps), 2001 (one map) and 2002 (two maps).

We present the first magnetic Doppler images of a rapidly oscillating Ap (roAp) star. We deduce information about magnetic field geometry and abundance distributions of a number of chemical elements on the surface of the hitherto best studied roAp star, HD 24712, using the magnetic Doppler imaging (MDI) code, INVERS10, which allows us to reconstruct simultaneously and consistently the magnetic field geometry and elemental abundance distributions on a stellar surface. For this purpose we analyse time series spectra obtained in Stokes I and V parameters with the SOFIN polarimeter at the Nordic Optical Telescope and recover surface abundance structures of sixteen different chemical elements, respectively ions, including Mg, Ca, Sc, Ti, Cr, Fe, Co, Ni, Y, La, Ce, Pr, Nd, Gd, Tb, and Dy. For the rare earth elements (REE) Pr and Nd separate maps were obtained using lines of the first and the second ionization stage. We find and confirm a clear dipolar structure of the surface magnetic field and an unexpected correlation of elemental abundances with respect to this field: one group of elements accumulates solely where the positive magnetic pole is visible, whereas the other group avoids this region and is enhanced where the magnetic equatorial region dominates the visible stellar surface. We also observe relative shifts of abundance enhancement- or depletion regions between the various elements exhibiting otherwise similar behaviour.

As accretion in protoplanetary disks is enabled by turbulent viscosity, the border between active and inactive (dead) zones constitutes a location where there is an abrupt change in the accretion flow. The gas accumulation that ensues triggers the Rossby wave instability, that in turn saturates into anticyclonic vortices. It was suggested that the trapping of solids within them leads to a burst of planet formation on very short timescales. We perform two-dimensional global simulations of the dynamics of gas and solids in a non-magnetized thin protoplanetary disk with the Pencil Code. We use multiple particle species of radius 1, 10, 30, and 100 cm, solving for the particles' gravitational interaction by a particle-mesh method. The dead zone is modeled as a region of low viscosity. Adiabatic and locally isothermal equations of state are used. We find that the Rossby wave instability is triggered under a variety of conditions, thus making vortex formation a robust process. Inside the vortices, fast accumulation of solids occurs and the particles collapse into objects of planetary mass in timescales as short as five orbits. Because the drag force is size-dependent, aerodynamical sorting ensues within the vortical motion, and the first bound structures formed are composed primarily of similarly-sized particles. In addition to erosion due to ram pressure, we identify gas tides from the massive vortices as a disrupting agent of formed protoplanetary embryos. We also estimate the collisional velocity history of the particles that compose the most massive embryo by the end of the simulation, finding that the vast majority of them never experienced a collision with another particle at speeds faster than 1 m/s.

We present the results on the photometric and spectroscopic monitoring of a luminous Ap star HD103498. The time-series photometric observations were carried out on 17 nights using three-channel fast photometer attached to the 1.04-m optical telescope at ARIES, Nainital. The photometric data of five nights of year 2007 show clear signature of 15-min periodicity. However, the follow-up observations during 2007--2009 could not repeated any such periodicity. To confirm the photometric light variations, the time-series spectroscopic observations were carried out with the 2.56-m Nordic Optical Telescope (NOT) at La Palma on February 2, 2009. Any radial velocity variations were absent in this data set which is in full agreement with the photometric observations taken near the same night. Model atmosphere and abundance analysis of HD103498 show that the star is evolved from the Main Sequence and its atmospheric abundances are similar to two other evolved Ap stars HD133792 and HD204411: large overabundances of Si, Cr, and Fe and moderate overabundances of the rare-earth elements. These chemical properties and a higher effective temperature distinguish HD103498 from any known roAp star.

We report the discovery of the first low-mass planet to emerge from the NASA-UC Eta-Earth Program, a super-Earth orbiting the K0 dwarf HD 7924. Keplerian modeling of precise Doppler radial velocities reveals a planet with minimum mass M_P sin i = 9.26 M_Earth in a P = 5.398 d orbit. Based on Keck-HIRES measurements from 2001 to 2008, the planet is robustly detected with an estimated false alarm probability of less than 0.001. Photometric observations using the Automated Photometric Telescopes at Fairborn Observatory show that HD 7924 is photometrically constant over the radial velocity period to 0.19 mmag, supporting the existence of the planetary companion. No transits were detected down to a photometric limit of ~0.5 mmag, eliminating transiting planets with a variety of compositions. HD 7924b is one of only eight planets known with M_P sin i < 10 M_Earth and as such is a member of an emerging family of low-mass planets that together constrain theories of planet formation.

Centimeter and meter sized solid particles in protoplanetary disks are trapped within long lived high pressure regions, creating opportunities for collapse into planetesimals and planetary embryos. We study the accumulations in the stable Lagrangian points of a giant planet, as well as in the Rossby vortices launched at the edges of the gap it carves. We employ the Pencil Code, tracing the solids with a large number of interacting Lagrangian particles, usually 100,000. For particles of 1 cm to 10 cm radii, gravitational collapse occurs in the Lagrangian points in less than 200 orbits. For 5 cm particles, a 2 Earth mass planet is formed. For 10 cm, the final maximum collapsed mass is around 3 Earth masses. The collapse of the 1 cm particles is indirect, following the timescale of depletion of gas from the tadpole orbits. In the edges of the gap vortices are excited, trapping preferentially particles of 30 cm radii. The rocky planet that is formed is as massive as 17 Earth masses, constituting a Super-Earth. By using multiple particle species, we find that gas drag modifies the streamlines in the tadpole region around the classical L4 and L5 points. As a result, particles of different radii have their stable points shifted to different locations. Collapse therefore takes longer and produces planets of lower mass. Three super-Earths are formed in the vortices, the most massive having 4.4 Earth masses. We conclude that a Jupiter mass planet can induce the formation of other planetary embryos in the outer edge of its gas gap. Trojan Earth mass planets are readily formed, and although not existing in the solar system, might be common in the exoplanetary zoo.

In the borders of the dead zones of protoplanetary disks, the inflow of gas produces a local density maximum that triggers the Rossby wave instability. The vortices that form are efficient in trapping solids. We aim to assess the possibility of gravitational collapse of the solids within the Rossby vortices. We perform global simulations of the dynamics of gas and solids in a low mass non-magnetized self-gravitating thin protoplanetary disk with the Pencil code. We use multiple particle species of radius 1, 10, 30, and 100 cm. The dead zone is modeled as a region of low viscosity. The Rossby vortices excited in the edges of the dead zone are very efficient particle traps. Within 5 orbits after their appearance, the solids achieve critical density and undergo gravitational collapse into Mars sized objects. The velocity dispersions are of the order of 10 m/s for newly formed embryos, later lowering to less than 1 m/s by drag force cooling. After 200 orbits, 38 gravitationally bound embryos were formed inside the vortices, half of them being more massive than Mars. The embryos are composed primarily of same-sized particles. We conclude that the presence of a dead zone naturally gives rise to a population of protoplanetary cores in the mass range of 0.1-0.6 Earth masses, on very short timescales.

We present global 3D MHD simulations of disks of gas and solids, aiming at developing models that can be used to study various scenarios of planet formation and planet-disk interaction in turbulent accretion disks. A second goal is to show that Cartesian codes are comparable to cylindrical and spherical ones in handling the magnetohydrodynamics of the disk simulations, as the disk-in-a-box models presented here develop and sustain MHD turbulence. We investigate the dependence of the magnetorotational instability on disk scale height, finding evidence that the turbulence generated by the magnetorotational instability grows with thermal pressure. The turbulent stresses depend on the thermal pressure obeying a power law of 0.24+/-0.03, compatible with the value of 0.25 found in shearing box calculations. The ratio of stresses decreased with increasing temperature. We also study the dynamics of boulders in the hydromagnetic turbulence. The vertical turbulent diffusion of the embedded boulders is comparable to the turbulent viscosity of the flow. Significant overdensities arise in the solid component as boulders concentrate in high pressure regions.

We present a new implementation of the magnetic Doppler imaging technique, which aims at self-consistent temperature and magnetic mapping of the surface structures in cool active stars. Our magnetic imaging procedure is unique in its capability to model individual spectral features in all four Stokes parameters. We discuss performance and intrinsic limitations of the new magnetic Doppler imaging method. A special emphasis is given to the simultaneous modelling of the magnetically sensitive lines in the optical and infrared regions and to combining information from both atomic and molecular spectral features. These two techniques may, for the first time, give us a tool to study magnetic fields in the starspot interiors.

The exoplanets discovered so far have been mostly around relatively nearby and bright stars. As a result, the host stars are mostly (i) in the Galactic disk, (ii) relatively massive, and (iii) relatively metal rich. The aim of the SWEEPS project is to extend our knowledge to stars which (i) are in a different part of the Galaxy, (ii) have lower masses, and (iii) have a large range of metallicities. To achieve this goal, we used the Hubble Space Telescope to search for transiting planets around F, G, K, and M dwarfs in the Galactic bulge. We photometrically monitored 180,000 stars in a dense bulge field continuously for 7 days. We discovered 16 candidate transiting extrasolar planets with periods of 0.6 to 4.2 days, including a new class of ultra-short period planets (USPPs) with P < 1.2 days. Radial-velocity observations of the two brightest candidates support their planetary nature. These results suggest that planets are as abundant in the Galactic bulge as they are in the solar neighborhood, and they are equally abundant around low-mass stars (within a factor 2). The planet frequency increases with metallicity even for the stars in the Galactic bulge. All the USPP hosts are low-mass stars, suggesting either that close-in planets around higher-mass stars are irradiatively evaporated, or that the planets can migrate to close-in orbits only around such old and low-mass stars.

We present a homogeneous photometric and spectroscopic analysis of 18 stars along the evolutionary sequence of the metal-poor globular cluster NGC 6397 ([Fe/H] = -2), from the main-sequence turnoff point to red giants below the bump. The spectroscopic stellar parameters, in particular stellar-parameter differences between groups of stars, are in good agreement with broad-band and Stroemgren photometry calibrated on the infrared-flux method. The spectroscopic abundance analysis reveals, for the first time, systematic trends of iron abundance with evolutionary stage. Iron is found to be 31% less abundant in the turnoff-point stars than in the red giants. An abundance difference in lithium is seen between the turnoff-point and warm subgiant stars. The impact of potential systematic errors on these abundance trends (stellar parameters, the hydrostatic and LTE approximations) is quantitatively evaluated and found not to alter our conclusions significantly. Trends for various elements (Li, Mg, Ca, Ti and Fe) are compared with stellar-structure models including the effects of atomic diffusion and radiative acceleration. Such models are found to describe the observed element-specific trends well, if extra (turbulent) mixing just below the convection zone is introduced. It is concluded that atomic diffusion and turbulent mixing are largely responsible for the sub-primordial stellar lithium abundances of warm halo stars. Other consequences of atomic diffusion in old metal-poor stars are also discussed.

The formation of long-lasting structures at the surfaces of stars is commonly ascribed to the action of strong magnetic fields. This paradigm is supported by observations of evolving cool spots in the Sun and active late-type stars, and stationary chemical spots in the early-type magnetic stars. However, results of our seven-year monitoring of mercury spots in non-magnetic early-type star alpha Andromedae show that the picture of magnetically-driven structure formation is fundamentally incomplete. Using an indirect stellar surface mapping technique, we construct a series of 2-D images of starspots and discover a secular evolution of the mercury cloud cover in this star. This remarkable structure formation process, observed for the first time in any star, is plausibly attributed to a non-equilibrium, dynamical evolution of the heavy-element clouds created by atomic diffusion and may have the same underlying physics as the weather patterns on terrestrial and giant planets.

Broadening and asymmetry of spectral lines in slowly rotating late A-type stars provide evidence for high-amplitude convective motions. The properties of turbulence observed in the A-star atmospheres are not understood theoretically and contradict results of previous numerical simulations of convection. Here we describe an ongoing effort to understand the puzzling convection signatures of A stars with the help of 3-D hydrodynamic simulations. Our approach combines realistic spectrum synthesis and non-grey hydrodynamic models computed with the CO5BOLD code. We discuss these theoretical predictions and confront them with high-resolution spectra of A stars. Our models have, for the first time, succeeded in reproducing the observed profiles of weak spectral lines without introducing fudge broadening parameters.

A homogeneous spectroscopic analysis of unevolved and evolved stars in the metal-poor globular cluster NGC 6397 with FLAMES-UVES reveals systematic trends of stellar surface abundances that are likely caused by atomic diffusion. This finding helps to understand, among other issues, why the lithium abundances of old halo stars are significantly lower than the abundance found to be produced shortly after the Big Bang.

More than 200 extrasolar planets have been discovered around relatively nearby stars, primarily through the Doppler line shifts owing to the reflex motions of their host stars, and more recently through transits of some planets across the face of the host stars. The detection of planets with the shortest known periods, 1.2 to 2.5 days, has mainly resulted from transit surveys which have generally targeted stars more massive than 0.75 M_sun. Here we report the results from a planetary transit search performed in a rich stellar field towards the Galactic bulge. We discovered 16 candidates with orbital periods between 0.4 and 4.2 days, five of which orbit stars of 0.44 to 0.75 M_sun. In two cases, radial-velocity measurements support the planetary nature of the companions. Five candidates have orbital periods below 1.0 day, constituting a new class of ultra-short-period planets (USPPs), which occur only around stars of less than 0.88 M_sun. This indicates that those orbiting very close to more luminous stars might be evaporatively destroyed, or that jovian planets around lower-mass stars might migrate to smaller radii.

We analyse high-resolution archival UVES data of turnoff and subgiant stars in the nearby globular cluster NGC 6397 ([Fe/H] = -2). Balmer-profile analyses are performed to derive reddening-free effective temperatures. Due to the limited S/N and uncertainties related to blaze removal, we find the data quality insufficient to exclude the existence of gravitational settling. If the newly derived effective temperatures are taken as a basis for an abundance analysis, the photospheric iron (Fe II) abundance in the turnoff stars is 0.11 dex lower than in the (well-mixed) subgiants.

The measurement of the cosmic microwave background has strongly constrained the cosmological parameters of the Universe. When the measured density of baryons (ordinary matter) is combined with standard Big Bang nucleosynthesis calculations, the amounts of hydrogen, helium and lithium produced shortly after the Big Bang can be predicted with unprecedented precision. The predicted primordial lithium abundance is a factor of two to three higher than the value measured in the atmospheres of old stars. With estimated errors of 10 to 25%, this cosmological lithium discrepancy seriously challenges our understanding of stellar physics, Big Bang nucleosynthesis or both. Certain modifications to nucleosynthesis have been proposed, but found experimentally not to be viable. Diffusion theory, however, predicts atmospheric abundances of stars to vary with time, which offers a possible explanation of the discrepancy. Here we report spectroscopic observations of stars in the metalpoor globular cluster NGC 6397 that reveal trends of atmospheric abundance with evolutionary stage for various elements. These element-specific trends are reproduced by stellar-evolution models with diffusion and turbulent mixing. We thus conclude that diffusion is predominantly responsible for the low apparent stellar lithium abundance in the atmospheres of old stars by transporting the lithium deep into the star.

Starspots are usually associated with the action of magnetic fields at the stellar surfaces. However, recently an inhomogeneous chemical distribution of mercury was found for the mercury-manganese (HgMn) star alpha And -- a well-established member of a non-magnetic subclass of the chemically peculiar stars of the upper main sequence. In this study we present first results of the high-resolution survey of the HgII 3984 resonance line in the spectra of rapidly rotating HgMn stars with atmospheric parameters similar to those of alpha And. We use spectrum synthesis modelling and take advantage of the Doppler resolution of the stellar surfaces to probe horizontal structure of mercury distribution. Clear signatures of spots are found in the HgII 3984 line profiles of HR 1185 and HR 8723. Two observations of the latter star separated by two days give evidence for the line profile variability. We conclude that inhomogeneous distribution of Hg is a common phenomenon for the rapidly rotating HgMn stars in the 13000--13800 K effective temperature range independently of the stellar evolutionary stage. These results establish existence of a new class of spectrum variable spotted B-type stars. It is suggested that the observed Hg inhomogeneities arise from dynamical instabilities in the chemical diffusion processes and are unrelated to magnetic phenomena.

We present an analysis of 210 high-resolution time-resolved spectropolarimetric observations of the roAp star gamma Equ obtained over three nights in August and September 2003. Radial velocity variations due to p-mode non-radial pulsations are clearly detected in the lines of rare-earth elements, in particular Pr III, Nd II and Nd III. In contrast, we find absolutely no evidence for the variation of the mean longitudinal magnetic field over the pulsation period in gamma Equ at the level of 110-240 G which was recently reported by Leone & Kurtz (2003). Our investigation of the variability of circularly polarized profiles of 13 Nd III lines demonstrates that, at the 3 sigma confidence level, no magnetic field variation with an amplitude above 40-60 G was present in gamma Equ during our monitoring of this star.

An analysis of H alpha and H beta spectra in a sample of 30 cool dwarf and subgiant stars is presented using MARCS model atmospheres based on the most recent calculations of the line opacities. A detailed quantitative comparison of the solar flux spectra with model spectra shows that Balmer line profile shapes, and therefore the temperature structure in the line formation region, are best represented under the mixing length theory by any combination of a low mixing-length parameter alpha and a low convective structure parameter y. A slightly lower effective temperature is obtained for the sun than the accepted value, which we attribute to errors in models and line opacities. The programme stars span temperatures from 4800 to 7100 K and include a small number of population II stars. Effective temperatures have been derived using a quantitative fitting method with a detailed error analysis. Our temperatures find good agreement with those from the Infrared Flux Method (IRFM) near solar metallicity but show differences at low metallicity where the two available IRFM determinations themselves are in disagreement. Comparison with recent temperature determinations using Balmer lines by Fuhrmann (1998, 2000), who employed a different description of the wing absorption due to self-broadening, does not show the large differences predicted by Barklem et al. (2000). In fact, perhaps fortuitously, reasonable agreement is found near solar metallicity, while we find significantly cooler temperatures for low metallicity stars of around solar temperature.

Details of a theory of self-broadening of hydrogen lines are presented. The main features of the new theory are that the dispersive-inductive components of the interaction (van der Waals forces) have been included, and the resonance components have been computed by perturbation theory without the use of the multipole expansion. The theory is applied to lower Balmer lines and the theoretical and observational impact of the new broadening theory is examined. It is shown that this theory leads to considerable differences in the predicted line profiles in cool stars when compared with previous theories which include only resonance interactions. In particular, the effect is found to be very important in metal poor stars. The theory provides a natural explanation for the behaviour of effective temperatures derived from Balmer lines by others using a theory which includes only resonance broadening. When applied to Balmer lines in the solar spectrum the theory predicts an improved agreement between observed and computed profiles for models which also match limb darkening curves and rules out a model which does not. However significant discrepancies still remain which could be due to inadequacies in our theory or the atmospheric model or both.

In order to clarify the role of the magnetic field in generating abundance inhomogeneities in the atmospheres of Ap stars, we present new abundance Doppler images and an approximate magnetic field geometry for the Ap star iota Cas.

kappa Psc (HD 220825) is a typical Chromium Ap star that happens to have optimal parameters for Doppler imaging (DI). Its short rotational period of less then 2 days, rotational velocity of ~40 km/s, and a moderate inclination of the rotational axis put modest requirements on spectroscopic observations. Anomalies of iron peak elements are clearly present, but small enough to cause significant deviations from model atmospheres with scaled solar abundances. We applied DI to kappa Psc once before, determining the distribution of Cr (Ryabchikova et al. 1996, hereafter Paper I). However, due to strong blending of Fe, the image was based on two short (~2 Angstroems) spectral intervals, dominated by Cr lines. Since the first paper we obtained additional spectra and developed a new code that allows to perform multi-element DI and thus to use larger spectral interval(s). We demonstrate the abilities of the new code and present new maps of Cr and Fe. A much larger time base allowed us to improve the rotational period of kappa Psc as well.

We report the first identification of the Eu III lambda 6666.317 line in optical spectra of CP stars. This line is clearly present in the spectra of HR 4816, 73 Dra, HR 7575, and beta CrB, while it is marginally present or absent in spectra of the roAp stars alpha Cir, gamma Equ, BI Mic, 33 Lib, and HD 24712.

We describe a new method for the computation of opacity distribution functions (ODFs) useful to calculate one-dimensional model atmospheres in local thermal equilibrium (LTE). The new method is fast enough to be applied on current workstations and allows the computation of model atmospheres which deviate significantly from (scaled) solar chemical composition. It has reproduced existing ODFs and model atmospheres for solar abundances. Depending on the type of chemical peculiarity the ``individual'' model atmosphere may have a structure and surface fluxes similar to atmospheres based on (scaled) solar abundances or deviate in a way that cannot be reproduced by any of the conventional models. Examples are given to illustrate this behavior. The availability of models with ``individualized'' abundances is crucial for abundance analyses and Doppler imaging of extreme CP stars.

We present preliminary results from a programme aimed at acquiring linear spectropolarimetry of magnetic A and B stars. Linear polarization in the spectral lines of these objects is due to the Zeeman effect, and should provide detailed new information regarding the structure of their strong magnetic fields. To illustrate the impact of these new data, we compare observed circular and linear polarization line profiles of 53 Cam with the profiles predicted by the magnetic model by Landstreet (1988). Linear polarization in the spectral lines of all stars studied is extremely weak; in most cases, below the threshold of detectability even for very high SNRs. In order to overcome this problem, we employ the Least-Squares Deconvolution (LSD) multi-line analysis technique in order to extract low-noise mean line profiles and polarization signatures from our echelle spectra. Tests show that these mean signatures can be modelled as real spectral lines, and have the potential to lead to high-resolution maps of the magnetic and chemical abundance surface distributions.

The modern Doppler Imaging (DI) technique allows the reconstruction of different stellar surface structures based on accurate calculation of spectra of specific intensity. New applications like the mapping of the magnetic field vector put very stringent requirements on the radiative transfer (RT) solver which should be accurate, fast, and robust against numerical errors. We describe the evaluation of three different algorithms for our new magnetic DI code INVERS10. We also show the first results of numerical experiments made with the new code.