Source author record

Ansgar Reiners

Ansgar Reiners appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

59works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

59 published item(s)

preprint2026arXiv

Large-scale time-series spectroscopy for stellar ages

To date, Galactic Astronomy has largely concerned itself with astrophysical processes, and with the locations, space motions and compositions of objects. Consider, for example, the elucidation of the components of the Galaxy over the past decades, its mapping as enabled by Gaia and its predecessors, the photometric and spectroscopic characterization of innumerable astrophysical objects in various wavelength ranges, both from the ground and from space, and the expanding discovery and characterization of exoplanets; all focused on the current, static Galaxy. This White Paper proposes a dedicated program to derive stellar ages from time-series spectroscopy to hasten the transformation of this static conception into a dynamical one with age-labeled objects and events.

preprint2025arXiv

Revisiting the atmosphere of HAT-P-70b with CARMENES high-resolution transmission spectroscopy

Owing to hot and inflated envelopes that facilitate atmospheric studies, ultra-hot Jupiters (UHJs) have attracted much attention. Significant progress has been achieved, from enlarging the sample size to broadening the studies to encompass diverse stellar types and ages. Here, we present a transmission spectroscopy study of HAT-P-70b, an UHJ orbiting a young A-type star, through high-resolution observations with CARMENES at the 3.5m Calar Alto telescope. By using the line-by-line technique, we confirm the previous detections of Ha, Na I, and Ca II, report a new tentative detection of K I, and impose an upper limit on the He triplet absorption. Through cross-correlation analysis, we identify the Ca II and Fe I absorptions, both blue-shifted by approximately 5 km/s, indicating a day-to-night side atmospheric wind. Additionally, we find a new tentative detection of K I. We do not see any significant atmospheric molecular signal in the near-infrared data. Putting HAT-P-70b in the context of UHJs from the literature, it turns out that (1) Ha absorption is more common on gas giants orbiting stars younger than 1 Gyr, with a relative detection probability of $P_{\rm Age<1\,Gyr}({\rm Ha})/P_{\rm Age\geq1\,Gyr}({\rm Ha})\sim 3$; (2) any UHJ is likely to exhibit Fe I absorption if it has Ca II.

preprint2022arXiv

A transiting, temperate mini-Neptune orbiting the M dwarf TOI-1759 unveiled by TESS

We report the discovery and characterization of TOI-1759~b, a temperate (400 K) sub-Neptune-sized exoplanet orbiting the M~dwarf TOI-1759 (TIC 408636441). TOI-1759 b was observed by TESS to transit on sectors 16, 17 and 24, with only one transit observed per sector, creating an ambiguity on the orbital period of the planet candidate. Ground-based photometric observations, combined with radial-velocity measurements obtained with the CARMENES spectrograph, confirm an actual period of $18.85019 \pm 0.00014$ d. A joint analysis of all available photometry and radial velocities reveal a radius of $3.17 \pm 0.10\,R_\oplus$ and a mass of $10.8 \pm 1.5\,M_\oplus$. Combining this with the stellar properties derived for TOI-1759 ($R_\star = 0.597 \pm 0.015\,R_\odot$; $M_\star = 0.606 \pm 0.020\,M_\odot$; $T_{\textrm{eff}} = 4065 \pm 51$ K), we compute a transmission spectroscopic metric (TSM) value of over 80 for the planet, making it a good target for transmission spectroscopy studies. TOI-1759 b is among the top five temperate, small exoplanets ($T_\textrm{eq} < 500$ K, $R_p < 4 \,R_\oplus$) with the highest TSM discovered to date. Two additional signals with periods of 80 d and $>$ 200 d seem to be present in our radial velocities. While our data suggest both could arise from stellar activity, the later signal's source and periodicity are hard to pinpoint given the $\sim 200$ d baseline of our radial-velocity campaign with CARMENES. Longer baseline radial-velocity campaigns should be performed in order to unveil the true nature of this long period signal.

preprint2021arXiv

Moderately misaligned orbit of the warm sub-Saturn HD 332231 b

Measurements of exoplanetary orbital obliquity angles for different classes of planets are an essential tool in testing various planet formation theories. Measurements for those transiting planets on relatively large orbital periods ($P$\,$>$\,10\,d) present a rather difficult observational challenge. Here we present the obliquity measurement for the warm sub-Saturn planet HD\,332231\,b, which was discovered through Transiting Exoplanet Survey Satellite (TESS) photometry of sectors 14 and 15, on a relatively large orbital period (18.7\,d). Through a joint analysis of previously obtained spectroscopic data and our newly obtained CARMENES transit observations, we estimated the spin-orbit misalignment angle, \textlambda, to be $-42.0^{+11.3}_{-10.6}$\,$°$, which challenges Laplacian ideals of planet formation. Through the addition of these new radial velocity (RV) data points obtained with CARMENES, we also derived marginal improvements on other orbital and bulk parameters for the planet, as compared to previously published values. We showed the robustness of the obliquity measurement through model comparison with an aligned orbit. Finally, we demonstrated the inability of the obtained data to probe any possible extended atmosphere of the planet, due to a lack of precision, and place the atmosphere in the context of a parameter detection space.

preprint2020arXiv

HIRES, the high-resolution spectrograph for the ELT

HIRES will be the high-resolution spectrograph of the European Extremely Large Telescope at optical and near-infrared wavelengths. It consists of three fibre-fed spectrographs providing a wavelength coverage of 0.4-1.8 mic (goal 0.35-1.8 mic) at a spectral resolution of ~100,000. The fibre-feeding allows HIRES to have several, interchangeable observing modes including a SCAO module and a small diffraction-limited IFU in the NIR. Therefore, it will be able to operate both in seeing and diffraction-limited modes. ELT-HIRES has a wide range of science cases spanning nearly all areas of research in astrophysics and even fundamental physics. Some of the top science cases will be the detection of bio signatures from exoplanet atmospheres, finding the fingerprints of the first generation of stars (PopIII), tests on the stability of Nature's fundamental couplings, and the direct detection of the cosmic acceleration. The HIRES consortium is composed of more than 30 institutes from 14 countries, forming a team of more than 200 scientists and engineers.

preprint2020arXiv

Radial velocity constraints on the long-period transiting planet Kepler-1625 b with CARMENES

The star Kepler-1625 recently attracted considerable attention when an analysis of the stellar photometric time series from the Kepler mission was interpreted as showing evidence of a large exomoon around the transiting Jupiter-sized planet candidate Kepler-1625b. We aim to detect the radial velocity (RV) signal imposed by Kepler-1625b (and its putative moon) on the host star or, as the case may be, determine an upper limit on the mass of the transiting object. We took a total of 22 spectra of Kepler-1625 using CARMENES, 20 of which were useful. Observations were spread over a total of seven nights between October 2017 and October 2018, covering $125\%$ of one full orbit of Kepler-1625b. We used the automatic Spectral Radial Velocity Analyser (SERVAL) pipeline to deduce the stellar RVs and uncertainties. Then we fitted the RV curve model of a single planet on a Keplerian orbit to the observed RVs using a $χ^2$ minimisation procedure. We derive upper limits on the mass of Kepler-1625b under the assumption of a single planet on a circular orbit. In this scenario, the $1\,σ$, $2\,σ$, and $3\,σ$ confidence upper limits for the mass of Kepler-1625b are $2.90\,M_{\rm J}$, $7.15\,M_{\rm J}$, and $11.60\,M_{\rm J}$, respectively. We present strong evidence for the planetary nature of Kepler-1625b, making it the 10th most long-period confirmed planet known today. Our data does not answer the question about a second, possibly more short-period planet that could be responsible for the observed transit timing variation of Kepler-1625b.

preprint2020arXiv

The IAG Solar Flux Atlas: Telluric Correction With a Semi-Empirical Model

Observations of the Sun as a star have been key to guiding models of stellar atmospheres and additionally provide useful insights on the effects of granulation and stellar activity on radial velocity measurements. Most high resolution solar atlases contain telluric lines that span the optical and limit the spectral regions useful for analysis. We present here a telluric-corrected solar atlas covering 0.5-1.0 micron derived from solar spectra taken with a Fourier Transform Spectrograph (FTS) at the Institut für Astrophysik, Göttingen. This atlas is the highest resolution spectrum with a wavelength calibration precise to +/-10 m/s across this 500 nm spectral window. We find that the atlas matches to within 3% of the telluric-corrected Kitt Peak atlas in regions containing telluric absorption weaker than 50% in transmission. The telluric component of the spectral data is fit with a semi-empirical model composed of Lorentz profiles initialized to the HITRAN parameters for each absorption feature. Comparisons between the best-fit telluric parameters describing the Lorentz profile for each absorption feature to the original HITRAN values in general show excellent agreement considering the effects atmospheric pressure and temperature have on our final parameters. However, we identify a small subset of absorption features with larger offsets relative to the catalogued line parameters. We make our final solar atlas available online. We additionally make available the telluric spectra extracted from the data that, given the high resolution of the spectrum, would be useful for studying the time evolution of telluric line shapes and their impact on Doppler measurements.

preprint2020arXiv

The impact of unresolved magnetic spots on high precision radial velocity measurements

The Doppler method of exoplanet detection has been extremely successful, but suffers from contaminating noise from stellar activity. In this work a model of a rotating star with a magnetic field based on the geometry of the K2 star Epsilon Eridani is presented and used to estimate its effect on simulated radial velocity measurements. A number of different distributions of unresolved magnetic spots were simulated on top of the observed large-scale magnetic maps obtained from eight years of spectropolarimetric observations. The radial velocity signals due to the magnetic spots have amplitudes of up to 10 m s$^{-1}$, high enough to prevent the detection of planets under 20 Earth masses in temperate zones of solar type stars. We show that the radial velocity depends heavily on spot distribution. Our results emphasize that understanding stellar magnetic activity and spot distribution is crucial for detection of Earth analogues.

preprint2019arXiv

Radial velocity photon limits for the dwarf stars of spectral classes F--M

The determination of extrasolar planet masses with the radial velocity (RV) technique requires spectroscopic Doppler information from the planet's host star, which varies with stellar brightness and temperature. We analyze Doppler information in spectra of F--M dwarfs utilizing empirical information from HARPS and CARMENES, and from model spectra. We come to the conclusions that an optical setup ($BVR$-bands) is more efficient that a near-infrared one ($YJHK$) in dwarf stars hotter than 3200\,K. We publish a catalogue of 46,480 well-studied F--M dwarfs in the solar neighborhood and compare their distribution to more than one million stars from Gaia DR2. For all stars, we estimate the RV photon noise achievable in typical observations assuming no activity jitter and slow rotation. We find that with an ESPRESSO-like instrument at an 8m-telescope, a photon noise limit of 10\,cm\,s$^{-1}$ or lower can be reached in more than 280 stars in a 5\,min observation. At 4m-telescopes, a photon noise limit of 1\,m\,s$^{-1}$ can be reached in a 10\,min exposure in approx.\ 10,000 predominantly sun-like stars with a HARPS-like (optical) instrument. The same applies to $\sim$3000 stars for a red-optical setup covering the $RIz$-bands, and to $\sim$700 stars for a near-infrared instrument. For the latter two, many of the targets are nearby M dwarfs. Finally, we identify targets in which Earth-mass planets within the liquid water habitable zone can cause RV amplitudes comparable to the RV photon noise. Assuming the same exposure times, we find that an ESPRESSO-like instrument can reach this limit for 1\,M$_\Earth$ planets in more than 1000 stars. The optical, red-optical, and near-infrared configurations reach the limit for 2\,M$_\Earth$ planets in approximately 500, 700, and 200 stars, respectively.

preprint2016arXiv

A CRIRES-search for H3+ emission from the hot Jupiter atmosphere of HD 209458 b

Close-in extrasolar giant planets are expected to cool their thermospheres by producing H3+ emission in the near-infrared (NIR), but simulations predict H3+ emission intensities that differ in the resulting intensity by several orders of magnitude. We want to test the observability of H3+ emission with CRIRES at the Very Large Telescope (VLT), providing adequate spectral resolution for planetary atmospheric lines in NIR spectra. We search for signatures of planetary H3+ emission in the L` band, using spectra of HD 209458 obtained during and after secondary eclipse of its transiting planet HD 209458 b. We searched for H3+ emission signatures in spectra containing the combined light of the star and, possibly, the planet. With the information on the ephemeris of the transiting planet, we derive the radial velocities at the time of observation and search for the emission at the expected line positions and search for planetary signals and use a shift and add technique combining all observed spectra taken after sec. eclipse to calculate an upper emission limit. We do not find signatures of atmospheric H3+ emission in the spectra containing the combined light of HD 209458 and planet b. We calculate the emission limit for the H3+ line at 3953.0 nm (Q(1, 0)) to be 8.32 E18W and a limit of 5.34E18 W for the line at 3985.5 nm (Q(3, 0)). Comparing our emission limits to the theoretical predictions suggests that we lack 1 to 3 magnitudes of sensitivity to measure H3+ emission in our target object. We show that under more favorable weather conditions the data quality can be improved significantly, reaching 5 E16W for star-planet systems that are close to Earth. We estimate that pushing the detection limit down to 1E15W will be possible with ground-based observations with future instrumentation, for example, the E-ELT.

preprint2016arXiv

A terrestrial planet candidate in a temperate orbit around Proxima Centauri

At a distance of 1.295 parsecs, the red-dwarf Proxima Centauri ($α$ Centauri C, GL 551, HIP 70890, or simply Proxima) is the Sun's closest stellar neighbor and one of the best studied low-mass stars. It has an effective temperature of only $\sim$ 3050 K, a luminosity of $\sim$0.1 per cent solar, a measured radius of 0.14 R$_\odot$ and a mass of about 12 per cent the mass of the Sun. Although Proxima is considered a moderately active star, its rotation period is $\sim$ 83 days, and its quiescent activity levels and X-ray luminosity are comparable to the Sun's. New observations reveal the presence of a small planet orbiting Proxima with a minimum mass of 1.3~Earth masses and an orbital period of $\sim$11.2 days. Its orbital semi-major axis is $\sim0.05$ AU, with an equilibrium temperature in the range where water could be liquid on its surface.

preprint2016arXiv

Fundamental M-dwarf parameters from high-resolution spectra using PHOENIX ACES models: I. Parameter accuracy and benchmark stars

M-dwarf stars are the most numerous stars in the Universe; they span a wide range in mass and are in the focus of ongoing and planned exoplanet surveys. To investigate and understand their physical nature, detailed spectral information and accurate stellar models are needed. We use a new synthetic atmosphere model generation and compare model spectra to observations. To test the model accuracy, we compared the models to four benchmark stars with atmospheric parameters for which independent information from interferometric radius measurements is available. We used $χ^2$ -based methods to determine parameters from high-resolution spectroscopic observations. Our synthetic spectra are based on the new PHOENIX grid that uses the ACES description for the equation of state. This is a model generation expected to be especially suitable for the low-temperature atmospheres. We identified suitable spectral tracers of atmospheric parameters and determined the uncertainties in $T_{\rm eff}$, $\log{g}$, and [Fe/H] resulting from degeneracies between parameters and from shortcomings of the model atmospheres. The inherent uncertainties we find are σ$T_{\rm eff}$= 35 K, σ$\log{g}$ = 0.14, and σ[Fe/H] = 0.11. The new model spectra achieve a reliable match to our observed data; our results for $T_{\rm eff}$ and $\log{g}$ are consistent with literature values to within 1σ. However, metallicities reported from earlier photometric and spectroscopic calibrations in some cases disagree with our results by more than 3 σ. A possible explanation are systematic errors in earlier metallicity determinations that were based on insufficient descriptions of the cool atmospheres. At this point, however, we cannot definitely identify the reason for this discrepancy, but our analysis indicates that there is a large uncertainty in the accuracy of M-dwarf parameter estimates.

preprint2016arXiv

Radial velocity observations of the 2015 Mar 20 eclipse - A benchmark Rossiter-McLaughlin curve with zero free parameters

On March 20, 2015, we obtained 159 spectra of the Sun as a star with the solar telescope and the FTS at the Institut für Astrophysik Göttingen, 76 spectra were taken during partial solar eclipse. We obtained RVs using $I_2$ as wavelength reference and determined the RM curve with a peak-to-peak amplitude of almost 1.4 km s$^{-1}$ at typical RV precision better than 1 m s$^{-1}$. We modeled disk-integrated solar RVs using surface velocities, limb darkening, and information about convective blueshift from 3D magneto-hydrodynamic simulations. We confirm that convective blueshift is crucial to understand solar RVs during eclipse. Our best model reproduced the observations to within a relative precision of 10% with residuals less than 30 m s$^{-1}$. We cross-checked parameterizations of velocity fields using a Dopplergram from the Solar Dynamics Observatory and conclude that disk-integration of the Dopplergram does not provide correct information about convective blueshift necessary for m s$^{-1}$ RV work. As main limitation for modeling RVs during eclipses, we identified limited knowledge about convective blueshift and line shape as functions of solar limb angle. We suspect that our model line profiles are too shallow at limb angles larger than $μ= 0.6$ resulting in incorrect weighting of the velocities across the solar disk. Alternative explanations cannot be excluded like suppression of convection in magnetic areas and undiscovered systematics during eclipse observations. Accurate observations of solar line profiles across the solar disk are suggested. We publish our RVs taken during solar eclipse as a benchmark curve for codes calculating the RM effect and for models of solar surface velocities and line profiles.

preprint2016arXiv

State of the Field: Extreme Precision Radial Velocities

The Second Workshop on Extreme Precision Radial Velocities defined circa 2015 the state of the art Doppler precision and identified the critical path challenges for reaching 10 cm/s measurement precision. The presentations and discussion of key issues for instrumentation and data analysis and the workshop recommendations for achieving this precision are summarized here. Beginning with the HARPS spectrograph, technological advances for precision radial velocity measurements have focused on building extremely stable instruments. To reach still higher precision, future spectrometers will need to produce even higher fidelity spectra. This should be possible with improved environmental control, greater stability in the illumination of the spectrometer optics, better detectors, more precise wavelength calibration, and broader bandwidth spectra. Key data analysis challenges for the precision radial velocity community include distinguishing center of mass Keplerian motion from photospheric velocities, and the proper treatment of telluric contamination. Success here is coupled to the instrument design, but also requires the implementation of robust statistical and modeling techniques. Center of mass velocities produce Doppler shifts that affect every line identically, while photospheric velocities produce line profile asymmetries with wavelength and temporal dependencies that are different from Keplerian signals. Exoplanets are an important subfield of astronomy and there has been an impressive rate of discovery over the past two decades. Higher precision radial velocity measurements are required to serve as a discovery technique for potentially habitable worlds and to characterize detections from transit missions. The future of exoplanet science has very different trajectories depending on the precision that can ultimately be achieved with Doppler measurements.

preprint2016arXiv

The Göttingen Solar Radial Velocity Project: Sub-m/s Doppler precision from FTS observations of the Sun as a star

Radial velocity observations of stars are entering the sub-m/s domain revealing fundamental barriers for Doppler precision experiments. Observations of the Sun as a star can easily overcome the m/s photon limit but face other obstacles. We introduce the Göttingen Solar Radial Velocity Project with the goal to obtain high precision (cm/s) radial velocity measurements of the Sun as a star with a Fourier Transform Spectrograph. In this first paper, we present the project and first results. The photon limit of our 2 min observations is at the 2 cm/s level but currently limited by strong instrumental systematics. A drift of a few m/s per h is visible in all observing days probably caused by vignetting of the solar disk in our fiber coupled setup, and imperfections of our guiding system adds further offsets in our data. Binning the data into 30 min groups shows m/s stability after correcting for a daily and linear instrumental trend. Our results show the potential of Sun-as-a-star radial velocity measurements that can possibly be achieved after a substantial upgrade of our spectrograph coupling strategy. Sun-as-a-star observations can provide crucial empirical information about the radial velocity signal of convective motion and stellar activity, and on the wavelength dependence of radial velocity signals caused by stellar line profile variations.

preprint2016arXiv

The habitability of Proxima Centauri b. I. Irradiation, rotation and volatile inventory from formation to the present

Proxima b is a planet with a minimum mass of 1.3 MEarth orbiting within the habitable zone (HZ) of Proxima Centauri, a very low-mass, active star and the Sun's closest neighbor. Here we investigate a number of factors related to the potential habitability of Proxima b and its ability to maintain liquid water on its surface. We set the stage by estimating the current high-energy irradiance of the planet and show that the planet currently receives 30 times more EUV radiation than Earth and 250 times more X-rays. We compute the time evolution of the star's spectrum, which is essential for modeling the flux received over Proxima b's lifetime. We also show that Proxima b's obliquity is likely null and its spin is either synchronous or in a 3:2 spin-orbit resonance, depending on the planet's eccentricity and level of triaxiality. Next we consider the evolution of Proxima b's water inventory. We use our spectral energy distribution to compute the hydrogen loss from the planet with an improved energy-limited escape formalism. Despite the high level of stellar activity we find that Proxima b is likely to have lost less than an Earth ocean's worth of hydrogen before it reached the HZ 100-200 Myr after its formation. The largest uncertainty in our work is the initial water budget, which is not constrained by planet formation models. We conclude that Proxima b is a viable candidate habitable planet.

preprint2015arXiv

Calibrating echelle spectrographs with Fabry-Perot etalons

Over the past decades hollow-cathode lamps have been calibration standards for spectroscopic measurements. Advancing to cm/s radial velocity precisions with the next generation of instruments requires more suitable calibration sources with more lines and less dynamic range problems. Fabry-Perot interferometers provide a regular and dense grid of lines and homogeneous amplitudes making them good candidates for next generation calibrators. We investigate the usefulness of Fabry-Perot etalons in wavelength calibration, present an algorithm to incorporate the etalon spectrum in the wavelength solution and examine potential problems. The quasi periodic pattern of Fabry-Perot lines is used along with a hollow-cathode lamp to anchor the numerous spectral features on an absolute scale. We test our method with the HARPS spectrograph and compare our wavelength solution to the one derived from a laser frequency comb. The combined hollow-cathode lamp/etalon calibration overcomes large distortion (50 m/s) in the wavelength solution of the HARPS data reduction software. Direct comparison to the laser frequency comb bears differences of only maximum 10 m/s. Combining hollow-cathode lamps with Fabry-Perot Interferometers can lead to substantial improvements in the wavelength calibration of echelle spectrographs. Etalons can provide economical alternatives to the laser frequency comb, especially for smaller projects.

preprint2015arXiv

Effect of shear and magnetic field on the heat-transfer efficiency of convection in rotating spherical shells

We study rotating thermal convection in spherical shells. We base our analysis on a set of about 450 direct numerical simulations of the (magneto)hydrodynamic equations under the Boussinesq approximation. The Ekman number ranges from $10^{-3}$ to $10^{-5}$. The supercriticality of the convection reaches about 1000 in some models. Four sets of simulations are considered: non-magnetic simulations and dynamo simulations with either free-slip or no-slip flow boundary conditions. The non-magnetic setup with free-slip boundaries generates the strongest zonal flows. Both non-magnetic simulations with no-slip flow boundary conditions and self-consistent dynamos with free-slip boundaries have drastically reduced zonal-flows. Suppression of shear leads to a substantial gain in heat-transfer efficiency, increasing by a factor of 3 in some cases. Such efficiency enhancement occurs as long as the convection is significantly influenced by rotation. At higher convective driving the heat-transfer efficiency tends towards that of the classical non-rotating Rayleigh-Bénard system. Analysis of the latitudinal distribution of heat flow at the outer boundary reveals that the shear is most effective at suppressing heat-transfer in the equatorial regions. Furthermore, we explore the influence of the magnetic field on the {\em non-zonal} flow components of the convection. For this we compare the heat-transfer efficiency of no-slip non-magnetic cases with that of the no-slip dynamo simulations. We find that at $E=10^{-5}$ magnetic field significantly affects the convection and a maximum gain of about 30\% (as compared to the non-magnetic case) in heat-transfer efficiency is obtained for an Elsasser number of about 3. Our analysis motivates us to speculate that convection in the polar regions in dynamos at $E=10^{-5}$ is probably in a `magnetostrophic' regime.

preprint2015arXiv

Explaining the coexistence of large-scale and small-scale magnetic fields in fully convective stars

Despite the lack of a shear-rich tachocline region low-mass fully convective stars are capable of generating strong magnetic fields, indicating that a dynamo mechanism fundamentally different from the solar dynamo is at work in these objects. We present a self-consistent three dimensional model of magnetic field generation in low-mass fully convective stars. The model utilizes the anelastic magnetohydrodynamic equations to simulate compressible convection in a rotating sphere. A distributed dynamo working in the model spontaneously produces a dipole-dominated surface magnetic field of the observed strength. The interaction of this field with the turbulent convection in outer layers shreds it, producing small-scale fields that carry most of the magnetic flux. The Zeeman-Doppler-Imaging technique applied to synthetic spectropolarimetric data based on our model recovers most of the large-scale field. Our model simultaneously reproduces the morphology and magnitude of the large-scale field as well as the magnitude of the small-scale field observed on low-mass fully convective stars.

preprint2015arXiv

No evidence for activity correlations in the radial velocities of Kapteyn's star

Stellar activity may induce Doppler variability at the level of a few m/s which can then be confused by the Doppler signal of an exoplanet orbiting the star. To first order, linear correlations between radial velocity measurements and activity indices have been proposed to account for any such correlation. The likely presence of two super-Earths orbiting Kapteyn's star was reported in Anglada et al. (2014, MNRAS 443L, 89A), but this claim was recently challenged by Robertson et al. (2015, ApJ 805L, 22R) arguing evidence of a rotation period (143 days) at three times the orbital period of one of the proposed planets (Kapteyn's b, P=48.6 days), and the existence of strong linear correlations between its Doppler signal and activity data. By re-analyzing the data using global optimization methods and model comparison, we show that such claim is incorrect given that; 1) the choice of a rotation period at 143 days is unjustified, and 2) the presence of linear correlations is not supported by the data. We conclude that the radial velocity signals of Kapteyn's star remain more simply explained by the presence of two super-Earth candidates orbiting it. We also advocate for the use of global optimization procedures and objective arguments, instead of claims lacking of a minimal statistical support.

preprint2015arXiv

Precise radial velocities of giant stars VIII. Testing for the presence of planets with CRIRES Infrared Radial Velocities

We have been monitoring 373 very bright (V < 6 mag) G and K giants with high precision optical Doppler spectroscopy for more than a decade at Lick Observatory. Our goal was to discover planetary companions around those stars and to better understand planet formation and evolution around intermediate-mass stars. However, in principle, long-term, g-mode nonradial stellar pulsations or rotating stellar features, such as spots, could effectively mimic a planetary signal in the radial velocity data. Our goal is to compare optical and infrared radial velocities for those stars with periodic radial velocity patterns and to test for consistency of their fitted radial velocity semiamplitudes. Thereby, we distinguish processes intrinsic to the star from orbiting companions as reason for the radial velocity periodicity observed in the optical. Stellar spectra with high spectral resolution have been taken in the H-band with the CRIRES near-infrared spectrograph at ESO's VLT for 20 stars of our Lick survey. Radial velocities are derived using many deep and stable telluric CO2 lines for precise wavelength calibration. We find that the optical and near-infrared radial velocities of the giant stars in our sample are consistent. We present detailed results for eight stars in our sample previously reported to have planets or brown dwarf companions. All eight stars passed the infrared test. We conclude that the planet hypothesis provides the best explanation for the periodic radial velocity patterns observed for these giant stars.

preprint2015arXiv

The IAG solar flux atlas: Accurate wavelengths and absolute convective blueshift in standard solar spectra

We present a new solar flux atlas with the aim to understand wavelength precision and accuracy in solar benchmark data. The atlas covers the wavelength range 405--2300 nm and was observed at the Institut für Astrophysik, Göttingen (IAG) with a Fourier Transform Spectrograph. In contrast to other FTS atlases, the entire visible wavelength range was observed simultaneously using only one spectrograph setting. We compare the wavelength solution of the new atlas to the Kitt Peak solar flux atlases and to the HARPS frequency-comb calibrated solar atlas. Comparison reveals systematics in the two Kitt Peak FTS atlases resulting from their wavelength scale construction, and shows consistency between the IAG and the HARPS atlas. We conclude that the IAG atlas is precise and accurate on the order of $\pm 10$ m s$^{-1}$ in the wavelength range 405--1065 nm while the Kitt Peak atlases show deviations as large as several ten to 100 m s$^{-1}$. We determine absolute convective blueshift across the spectrum from the IAG atlas and report slight differences relative to results from the Kitt Peak atlas that we attribute to the differences between wavelength scales. We conclude that benchmark solar data with accurate wavelength solution are crucial to better understand the effect of convection on stellar RV measurements, which is one of the main limitations of Doppler spectroscopy at m s$^{-1}$ precision.

preprint2015arXiv

Three-dimensional simulations of near-surface convection in main-sequence stars. III. The structure of small-scale magnetic flux concentrations

The convective envelopes of cool main-sequence stars harbour magnetic fields with a complex global and local structure. These fields affect the near-surface convection and the outer stellar atmospheres in many ways and are responsible for the observable magnetic activity of stars. Our aim is to understand the local structure in unipolar regions with moderate average magnetic flux density. These correspond to plage regions covering a substantial fraction of the surface of the Sun (and likely also the surface of other Sun-like stars) during periods of high magnetic activity. We analyse the results of 18 local-box magnetohydrodynamics simulations covering the upper layers of the convection zones and the photospheres of cool main-sequence stars of spectral types F to early M. The average vertical field in these simulations ranges from 20 to 500G. We find a substantial variation of the properties of the surface magnetoconvection between main-sequence stars of different spectral types. As a consequence of a reduced efficiency of the convective collapse of flux tubes, M dwarfs lack bright magnetic structures in unipolar regions of moderate field strength. The spatial correlation between velocity and the magnetic field as well as the lifetime of magnetic structures and their sizes relative to the granules vary significantly along the model sequence of stellar types.

preprint2015arXiv

Three-dimensional simulations of near-surface convection in main-sequence stars. IV. Effect of small-scale magnetic flux concentrations on centre-to-limb variation and spectral lines

Magnetic fields affect the local structure of the photosphere of stars. They can considerably influence the radiative properties near the optical surface, flow velocities, and the temperature and pressure profiles. We aim at understanding qualitatively the influence of small magnetic flux concentrations in unipolar plage regions on the centre-to-limb variation of the intensity and its contrast and on the shape of spectral line profiles in cool main-sequence stars. We analyse the bolometric and continuum intensity and its angular dependence of 24 radiative magnetohydrodynamic simulations of the near-surface layers of main-sequence stars with six different sets of stellar parameters (spectral types F to early M) and four different average magnetic field strengths (including the non-magnetic case). We also calculated disc-integrated profiles of three spectral lines. The small magnetic flux concentrations formed in the magnetic runs of simulations have a considerable impact on the intensity and its centre-to-limb variation. Spectral lines are not only broadened owing to the Zeeman effect, but are also strongly affected by the modified thermodynamical structure and flow patterns. This indirect magnetic impact on the line profiles is often bigger than that of the Zeeman effect. The effects of the magnetic field on the radiation leaving the star can be considerable and is not restricted to spectral line broadening and polarisation by the Zeeman effect. The inhomogeneous structure of the magnetic field on small length scales and its impact on (and spatial correlation with) the local thermodynamical structure and the flow field near the surface influence the measurement of the global field properties and stellar parameters. These effects need to be taken into account in the interpretation of observations.

preprint2014arXiv

Formation of starspots in self-consistent global dynamo models: Polar spots on cool stars

Observations of cool stars reveal dark spot-like features on their surfaces. Compared to sunspots, starspots can be bigger or cover a larger fraction of the stellar surface. While sunspots appear only at low latitudes, starspots are also found in polar regions, in particular on rapidly rotating stars. Sunspots are believed to result from the eruption of magnetic flux-tubes rising from the deep interior of the Sun. The strong magnetic field locally reduces convective heat transport to the solar surface. Such flux-tube models have also been invoked to explain starspot properties. However, these models use several simplifications and so far the generation of either sunspots or starspots has not been demonstrated in a self-consistent simulation of stellar magnetic convection. Here we show that direct numerical simulations of a distributed dynamo operating in a density-stratified rotating spherical shell can spontaneously generate cool spots. Convection in the interior of the model produces a large scale magnetic field which interacts with near surface granular convection leading to strong concentrations of magnetic flux and formation of starspots. Prerequisites for the formation of sizeable high-latitude spots in the model are sufficiently strong density stratification and rapid rotation. Our model presents an alternate mechanism for starspot formation by distributed dynamo action.

preprint2014arXiv

MHD simulations of near-surface convection in cool main-sequence stars

The solar photospheric magnetic field is highly structured owing to its interaction with the convective flows. Its local structure has a strong influence on the profiles of spectral lines not only by virtue of the Zeeman effect, but also through the modification of the thermodynamical structure (e.g. line weakening in hot small-scale magnetic structures). Many stars harbor surface magnetic fields comparable to or larger than the Sun at solar maximum. Therefore, a strong influence of the field on the surface convection and on spectral line profiles can be expected. We carried out 3D local-box MHD simulations of unipolar magnetized regions (average fields of 20, 100, and 500G) with parameters corresponding to six main-sequence stars (spectral types F3V to M2V). The influence of the magnetic field on the convection and the local thermodynamical structure were analyzed in detail. For three spectral lines, we determined the impact of the magnetic field on the disc-integrated Stokes-I profiles. Line weakening has in many cases a stronger impact on the spectral line profiles than the Zeeman effect. Moreover, for some stars, the correlation between the magnetic field and the vertical velocity strongly influences the line shapes. These effects can impair determinations of stellar magnetic fields since currently used methods neglect the local structure of the magnetic field and its interaction with the convective flows. The MHD simulations presented can be used to quantify these effects and thus help to improve magnetic field measurements of cool main-sequence stars.

preprint2014arXiv

Two planets around Kapteyn's star : a cold and a temperate super-Earth orbiting the nearest halo red-dwarf

Exoplanets of a few Earth masses can be now detected around nearby low-mass stars using Doppler spectroscopy. In this paper, we investigate the radial velocity variations of Kapteyn's star, which is both a sub-dwarf M-star and the nearest halo object to the Sun. The observations comprise archival and new HARPS, HIRES and PFS Doppler measurements. Two Doppler signals are detected at periods of 48 and 120 days using likelihood periodograms and a Bayesian analysis of the data. Using the same techniques, the activity indicies and archival ASAS-3 photometry show evidence for low-level activity periodicities of the order of several hundred days. However, there are no significant correlations with the radial velocity variations on the same time-scales. The inclusion of planetary Keplerian signals in the model results in levels of correlated and excess white noise that are remarkably low compared to younger G, K and M dwarfs. We conclude that Kapteyn's star is most probably orbited by two super-Earth mass planets, one of which is orbiting in its circumstellar habitable zone, becoming the oldest potentially habitable planet known to date. The presence and long-term survival of a planetary system seems a remarkable feat given the peculiar origin and kinematic history of Kapteyn's star. The detection of super-Earth mass planets around halo stars provides important insights into planet-formation processes in the early days of the Milky Way.

preprint2013arXiv

A dynamically-packed planetary system around GJ667C with three super-Earths in its habitable zone

Since low-mass stars have low luminosities, orbits at which liquid water can exist on Earth-sized planets are relatively close-in, which produces Doppler signals that are detectable using state-of-the-art Doppler spectroscopy. GJ 667C is already known to be orbited by two super-Earth candidates. We investigate whether the data supports the presence of additional companions. We obtain new Doppler measurements from HARPS extracted spectra and combined them with those obtained from the PFS and HIRES spectrographs. We used Bayesian and periodogram-based methods to re-assess the number of candidates and evaluated the confidence of each detection. Among other tests, we validated the planet candidates by analyzing correlations of each Doppler signal activity indices and investigate quasi-periodicity. Doppler measurements of GJ 667C are described better by six Keplerian-like signals: the two known candidates (b and c); three additional few-Earth mass candidates with periods of 92, 62 and 39 days (d, e and f); a cold super-Earth in a 260-day orbit (g) and tantalizing evidence of a $\sim$ 1 M$_\oplus$ object in a close-in orbit of 17 days (h). We explore whether long-term stable orbits are compatible with the data by integrating 8$\times 10^4$ solutions derived from the Bayesian samplings. The system consisting of six planets is compatible with dynamically stable configurations. As for the solar system, the most stable solutions do not contain mean-motion resonances and are described well by analytic Laplace-Lagrange solutions. The presence of a seventh planet (h) is supported by the fact that it appears squarely centered on the only island of stability left in the six-planet solution. Habitability assessments accounting for the stellar flux, as well as tidal dissipation effects, indicate that three (maybe four) planets are potentially habitable...

preprint2013arXiv

A fast and reliable method to measure stellar differential rotation from photometric data

Co-rotating spots at different latitudes on the stellar surface generate periodic photometric variability and can be useful proxies to detect Differential Rotation (DR). DR is a major ingredient of the solar dynamo but observations of stellar DR are rather sparse. In view of the Kepler space telescope we are interested in the detection of DR using photometric information of the star, and to develop a fast method to determine stellar DR from photometric data. We ran a large Monte-Carlo simulation of differentially rotating spotted stars with very different properties to investigate the detectability of DR. For different noise levels the resulting light curves are prewhitened using Lomb-Scargle periodograms to derive parameters for a global sine fit to detect periodicities. We show under what conditions DR can successfully be detected from photometric data, and in which cases the light curve provides insufficient or even misleading information on the stellar rotation law. In our simulations, the most significant period P1_{out} could be detected in 96.2% of all light curves. Detection of a second period close to P1_{out} is the signature of DR in our model. For the noise-free case, in 64.2% of all stars such a period was found. Calculating the measured latitudinal shear of two distinct spots α_{out}, and comparing it to the known original spot rotation rates shows that the real value is on average 3.2% lower. Comparing the total equator-to-pole shear $α$ to $α_{out}$ we find that $α$ is underestimated by 8.8%, esp. the detection of DR for stars with $α$ < 6% is challenging. Finally, we apply our method to four differentially rotating Kepler stars and find close agreement with results from detailed modeling. Our method is capable of measuring stellar rotation periods and detecting DR with relatively high accuracy and is suitable for large data sets.

preprint2013arXiv

A new extensive library of PHOENIX stellar atmospheres and synthetic spectra

We present a new library of high-resolution synthetic spectra based on the stellar atmosphere code PHOENIX that can be used for a wide range of applications of spectral analysis and stellar parameter synthesis. The spherical mode of PHOENIX was used to create model atmospheres and to derive detailed synthetic stellar spectra from them. We present a new self-consistent way of describing micro-turbulence for our model atmospheres. The synthetic spectra cover the wavelength range from 500AA to 50.000AA with resolutions of R=500.000 in the optical and near IR, R=100.000 in the IR and a step size of 0.1AA in the UV. The parameter space covers 2.300K<=Teff<=12.000K, 0.0<=log(g)<=+6.0, -4.0<=[Fe/H]<=+1.0, and -0.2<=[alpha/Fe]<=+1.2. The library is a work in progress and we expect to extend it up to Teff=25.000 K.

preprint2013arXiv

Magnetic Fields in Low-mass Stars: An Overview of Observational Biases

Stellar magnetic dynamos are driven by rotation, rapidly rotating stars produce stronger magnetic fields than slowly rotating stars do. The Zeeman effect is the most important indicator of magnetic fields, but Zeeman broadening must be disentangled from other broadening mechanisms, mainly rotation. The relations between rotation and magnetic field generation, between Doppler and Zeeman line broadening, and between rotation, stellar radius, and angular momentum evolution introduce several observational biases that affect our picture of stellar magnetism. In this overview, a few of these relations are explicitly shown, and the currently known distribution of field measurements is presented.

preprint2013arXiv

Photometric Variability in Kepler Target Stars. III. Comparison with the Sun on Different Timescales

We utilize Kepler data to study the precision differential photometric variability of solar-type and cooler stars at different timescales, ranging from half an hour to 3 months. We define a diagnostic that characterizes the median differential intensity change between data bins of a given timescale. We apply the same diagnostics to SOHO data that has been rendered comparable to Kepler. The Sun exhibits similar photometric variability on all timescales as comparable solar-type stars in the Kepler field (it is not unusually quiet). The previously-defined photometric "range" serves as our activity proxy (driven by starspot coverage). We revisit the fraction of comparable stars in the Kepler field that are more active than the Sun. The exact active fraction depends on what is meant by "more active than the Sun", and on the magnitude limit of the sample of stars considered. This active fraction is between a quarter and a third (depending on the timescale). We argue that a reliable result requires timescales of half a day or longer and stars brighter than Kepler magnitude of 14, otherwise non-stellar noise distorts it. We also analyze main sequence stars grouped by temperature from 6500-3500K. As one moves to cooler stars, the active fraction of stars becomes steadily larger (greater than 90% for early M dwarfs). The Sun is a good photometric model at all timescales for those cooler stars that have long-term variability within the span of solar variability.

preprint2013arXiv

Radial velocity signatures of Zeeman broadening

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 λ, 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.

preprint2013arXiv

Rotation and differential rotation of active Kepler stars

We present rotation periods for thousands of active stars in the Kepler field derived from Q3 data. In most cases a second period close to the rotation period was detected, which we interpreted as surface differential rotation (DR). Active stars were selected from the whole sample using the range of the variability amplitude. To detect different periods in the light curves we used the Lomb-Scargle periodogram in a pre-whitening approach to achieve parameters for a global sine fit. The most dominant periods from the fit were ascribed to different surface rotation periods, but spot evolution could also play a role. Due to the large number of stars the period errors were estimated in a statistical way. We thus cannot exclude the existence of false positives among our periods. In our sample of 40.661 active stars we found 24.124 rotation periods $P_1$ between 0.5-45 days. The distribution of stars with 0.5 < B-V < 1.0 and ages derived from angular momentum evolution that are younger than 300 Myr is consistent with a constant star-formation rate. A second period $P_2$ within $\pm30$% of the rotation period $P_1$ was found in 18.619 stars (77.2%). Attributing these two periods to DR we found that the relative shear $α=ΔΩ/Ω$ increases with rotation period, and slightly decreases with effective temperature. The absolute shear $ΔΩ$ slightly increases between $T_{eff}=3500-6000$ K. Above 6000 K $ΔΩ$ shows much larger scatter. We found weak dependence of $ΔΩ$ on rotation period. Latitudinal differential rotation measured for the first time in more than 18.000 stars provides a comprehensive picture of stellar surface shear, consistent with major predictions from mean-field theory. To what extent our observations are prone to false positives and selection bias is not fully explored, and needs to be addressed using more Kepler data.

preprint2013arXiv

Three-dimensional simulations of near-surface convection in main-sequence stars - I. Overall structure

The near-surface layers of cool main-sequence stars are structured by convective flows, which are overshooting into the atmosphere. The flows and the associated spatio-temporal variations of density and temperature affect spectral line profiles and thus have an impact on estimates of stellar properties such as effective temperature, gravitational acceleration, and abundances. We aim at identifying distinctive properties of the thermodynamic structure of the atmospheres of different stars and understand their causes. We ran comprehensive 3D radiation hydrodynamics simulations of the near-surface layers of six simulated stars of spectral type F3V to M2V with the MURaM code. We carry out a systematic parameter study of the mean stratifications, flow structures, and the energy flux in these stars.\par Results: We find monotonic trends along the lower main sequence in granule size, flow velocity, and intensity contrast. The convection in the M-star models differs substantially from that of the hotter stars, mainly owing to the more gradual transition from convective to radiative energy transport. While the basic mechanisms driving surface convection in cool stars are the same, the properties of the convection vary along the main sequence. Apart from monotonic trends in rms velocity, intensity contrast, granule size, etc., there is a transition between "naked" and "hidden" granulation around spectral type K5V caused by the (highly non-linear) temperature dependence of the opacity. These variations have to be taken into account when stellar parameters are derived from spectra.

preprint2013arXiv

Three-dimensional simulations of near-surface convection in main-sequence stars - II. Properties of granulation and spectral lines

The atmospheres of cool main-sequence stars are structured by convective flows from the convective envelope that penetrate the optically thin layers and lead to structuring of the stellar atmospheres analogous to solar granulation. The flows have considerable influence on the 3D structure of temperature and pressure and affect the profiles of spectral lines formed in the photosphere. For the set of six 3D radiative (M)HD simulations of cool main-sequence stars described in the first paper of this series, we analyse the near-surface layers. We aim at describing the properties of granulation of different stars and at quantifying the effects on spectral lines of the thermodynamic structure and flows of 3D convective atmospheres. We detected and tracked granules in brightness images from the simulations to analyse their statistical properties, as well as their evolution and lifetime. We calculated spatially resolved spectral line profiles using the line synthesis code SPINOR. To enable a comparison to stellar observations, we implemented a numerical disc-integration, which includes (differential) rotation. Although the stellar parameters change considerably along the model sequence, the properties of the granules are very similar. The impact of the 3D structure of the atmospheres on line profiles is measurable in disc-integrated spectra. Line asymmetries caused by convection are modulated by stellar rotation. The 3D structure of cool stellar atmospheres as shaped by convective flows has to be taken into account when using photospheric lines to determine stellar parameters.

preprint2012arXiv

A catalogue of rotation and activity in early-M stars

We present a catalogue of rotation and chromospheric activity in a sample of 334 M dwarfs of spectral types M0--M4.5 populating the parameter space around the boundary to full convection. We obtained high-resolution optical spectra for 206 targets and determined projected rotational velocity, vsini, and Halpha emission. The data are combined with measurements of vsini in field stars of the same spectral type from the literature. Our sample adds 157 new rotation measurements to the existing literature and almost doubles the sample of available vsini. The final sample provides a statistically meaningful picture of rotation and activity at the transition to full convection in the solar neighborhood. We confirm the steep rise in the fraction of active stars at the transition to full convection known from earlier work. In addition, we see a clear rise in rotational velocity in the same stars. In very few stars, no chromospheric activity but a detection of rotational broadening was reported. We argue that all of them are probably spurious detections; we conclude that in our sample all significantly rotating stars are active, and all active stars are significantly rotating. The rotation-activity relation is valid in partially and in fully convective stars. Thus, we do not observe any evidence for a transition from a rotationally dominated dynamo in partially convective stars to a rotation-independent turbulent dynamo in fully convective stars; turbulent dynamos in fully convective stars of spectral types around M4 are still driven by rotation. Finally, we compare projected rotational velocities of 33 stars to rotational periods derived from photometry in the literature and determine inclinations for a few of them.

preprint2012arXiv

Habitable-zone super-Earth candidate in a six-planet system around the K2.5V star HD 40307

The K2.5 dwarf HD 40307 has been reported to host three super-Earths. The system lacks massive planets and is therefore a potential candidate for having additional low-mass planetary companions. We re-derive Doppler measurements from public HARPS spectra of HD 40307 to confirm the significance of the reported signals using independent data analysis methods. We also investigate these measurements for additional low-amplitude signals. We used Bayesian analysis of our radial velocities to estimate the probability densities of different model parameters. We also estimated the relative probabilities of models with differing numbers of Keplerian signals and verified their significance using periodogram analyses. We investigated the relation of the detected signals with the chromospheric emission of the star. As previously reported for other objects, we found that radial velocity signals correlated with the S-index are strongly wavelength dependent. We identify two additional clear signals with periods of 34 and 51 days, both corresponding to planet candidates with minimum masses a few times that of the Earth. An additional sixth candidate is initially found at a period of 320 days. However, this signal correlates strongly with the chromospheric emission from the star and is also strongly wavelength dependent. When analysing the red half of the spectra only, the five putative planetary signals are recovered together with a very significant periodicity at about 200 days. This signal has a similar amplitude as the other new signals reported in the current work and corresponds to a planet candidate with M sin i = 7 Me (HD 40307 g). ...

preprint2012arXiv

New measurements of rotation and differential rotation in A-F stars: Are there two populations of differentially rotating stars?

We present projected rotational velocities and new measurements of the rotational profile of some 180 nearby stars with spectral types A-F. The overall broadening profile is derived analysing spectral line shape from hundreds of spectral lines by the method of least-squares deconvolution. Rigid and differential rotation can be distinguished in 56 cases. Ten stars with significant differential rotation rates are identified. As of now, 33 differential rotators detected by line profile analysis have been confirmed. The frequency of differential rotators decreases towards high effective temperature and rapid rotation. There is evidence for two populations of differential rotators with a gap in between at spectral type early-F. The gap can only partly be explained by an upper bound found for the horizontal shear of F stars. Apparently, the physical conditions of differential rotation change at early-F spectral types.

preprint2012arXiv

Observations of Cool-Star Magnetic Fields

Cool stars like the Sun harbor convection zones capable of producing substantial surface magnetic fields leading to stellar magnetic activity. The influence of stellar parameters like rotation, radius, and age on cool-star magnetism, and the importance of the shear layer between a radiative core and the convective envelope for the generation of magnetic fields are keys for our understanding of low-mass stellar dynamos, the solar dynamo, and also for other large-scale and planetary dynamos. Our observational picture of cool-star magnetic fields has improved tremendously over the last years. Sophisticated methods were developed to search for the subtle effects of magnetism, which are difficult to detect particularly in cool stars. With an emphasis on the assumptions and capabilities of modern methods used to measure magnetism in cool stars, I review the different techniques available for magnetic field measurements. I collect the analyses on cool-star magnetic fields and try to compare results from different methods, and I review empirical evidence that led to our current picture of magnetic fields and their generation in cool stars and brown dwarfs.

preprint2011arXiv

M-dwarf metallicities - A high-resolution spectroscopic study in the near infrared

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.

preprint2011arXiv

MHD simulations reveal crucial differences between solar and very-cool star magnetic structures

We carried out 3D radiative magnetohydrodynamic simulations of the convective and magnetic structure in the surface layers (uppermost part of the convection zone and photosphere) of main-sequence stars of spectral types F3 to M2. The simulation results were analyzed in terms of sizes and properties of the convection cells (granules) and magnetic flux concentrations as well as velocity, pressure, density, and temperature profiles. Our numerical simulations show for the first time a qualitative difference in the magneto-convection between solar-like stars and M dwarfs. Owing to higher surface gravity, lower opacity (resulting in higher density at optical depth unity), and more stable downflows, small-scale magnetic structures concentrate into pore-like configurations of reduced intensity. This implies that in very cool stars magnetic surface structures like plage regions and starspots significantly differ from the solar example. Such a difference would have major impact on the interpretation of Doppler imaging data and the analysis of M dwarf spectra.

preprint2011arXiv

Radius Dependent Angular Momentum Evolution in Low-Mass Stars. I

Angular momentum evolution in low-mass stars is determined by initial conditions during star formation, stellar structure evolution, and the behaviour of stellar magnetic fields. Here we show that the empirical picture of angular momentum evolution arises naturally if rotation is related to magnetic field strength instead of to magnetic flux, and formulate a corrected braking law based on this. Angular momentum evolution then becomes a strong function of stellar radius, explaining the main trends observed in open clusters and field stars at a few Gyr: the steep transition in rotation at the boundary to full convection arises primarily from the large change in radius across this boundary, and does not require changes in dynamo mode or field topology. Additionally, the data suggest transient core-envelope decoupling among solar-type stars, and field saturation at longer periods in very low-mass stars. For solar-type stars, our model is also in good agreement with the empirical Skumanich law. Finally, in further support of the theory, we show that the predicted age at which low-mass stars spin down from the saturated to unsaturated field regimes in our model corresponds remarkably well to the observed lifetime of magnetic activity in these stars.

preprint2010arXiv

Discovery of a nearby young brown dwarf binary candidate

In near-infrared NaCo observations of the young brown dwarf 2MASS J0041353-562112, we discovered a companion a little less than a magnitude fainter than the primary. The binary candidate has a separation of 143 mas, the spectral types are M6.5 and M9.0 for the two components. Colors and flux ratios are consistent with the components being located at the same distance minimizing the probability of the secondary being a background object. The brown dwarf is known to show Li absorption constraining the age to less than ~200 Myr, and it was suspected to show ongoing accretion, indicating an age as low as ~10 Myr. We estimate distance and orbital parameters of the binary as a function of age. For an age of 10 Myr, the distance to the system is 50 pc, the orbital period is 126 yr, and the masses of the components are ~30 and ~15 MJup. The binary brown dwarf fills a so far unoccupied region in the parameters mass and age; it is a valuable new benchmark object for brown dwarf atmospheric and evolutionary models.

preprint2010arXiv

Investigation of magnetically sensitive FeH lines

M-type stars exhibit strong magnetic fields towards decreasing effective temperatures. The measurement of these fields is complicated due to missing indicators. Molecular FeH lines provide an excellent means to determine magnetic field strengths from the Zeeman broadening of magnetically sensitive lines. Our aim is the investigation of possible dependencies of the amount of sensitivity to magnetic fields from rotational quantum number, branch, and the projection of the total angular momentum onto the internuclear axis (Ω). We also compare results from com-putations with those from observations. We use high resolution CRIRES spectra of the two M dwarfs GJ1002 (M5.5 inactive) and GJ1224 (M4.5 active). Individual lines are fitted by Gaussians and the obtained line depths and widths from the active and inactive star can be compared with each other. In this way, magnetically sensitive lines can be detected. For test purposes, we do the same with computed spectra of FeH. One with zero magnetic field and the other with a 2 kG magnetic field vector used in the disc integration (i.e. pure radial at the disc center). We found, in agreement with theory, that lines with high Ω show strong sensitivity to magnetic fields. No obvious correlation with branch or J was found, which was also expected for lines formed in intermediate Hund's case. The computations agreed in general well with the observations, but in many cases the individual splitting of certain lines can be very different to observations.

preprint2010arXiv

Kepler light curves and stellar rotational periods

The Kepler space telescope monitors over 156.000 stars with an unprecedented photometric precision. We are interested in stellar rotational periods which we find using Lomb-Scargle periodograms. This work focuses on the 306 exoplanet candidate host stars released on June 15, 2010. We present statistics on how many of them show periodic photometric variability, providing preliminary periods and estimates of stellar activity. In the future, our work will focus on spot evolution and differential rotation.

preprint2010arXiv

Magnetic Fields on Cool Stars

Magnetic fields are an important ingredient to cool star physics, and there is great interest in measuring fields and their geometry in order to understand stellar dynamos and their influence on star formation and stellar evolution. During the last few years, a large number of magnetic field measurements became available. Two main approaches are being followed to measure the Zeeman effect in cool stars; 1) the measurement of polarized light, for example to produce magnetic maps, and 2) the measurement of integrated Zeeman broadening to measure the average magnetic field strength on the stellar surface. This article briefly reviews the two methods and compares results between them that are now available for about a dozen M-type stars. It seems that we see a great variety of magnetic geometries and field strengths with typical average fields of a few kG in active M-type stars. The interpretation of geometries, however, has not yet led to a clear picture of magnetic dynamos and field configuration, and work is needed on more observational data but also on the fundamental understanding of our measurements.

preprint2010arXiv

On detectability of Zeeman broadening in optical spectra of F- and G-dwarfs

We investigate the detectability of Zeeman broadening in optical Stokes I spectra of slowly rotating sun-like stars. To this end, we apply the LTE spectral line inversion package SPINOR to very-high quality CES data and explore how fit quality depends on the average magnetic field, Bf . One-component (OC) and two-component (TC) models are adopted. In OC models, the entire surface is assumed to be magnetic. Under this assumption, we determine formal 3σ upper limits on the average magnetic field of 200 G for the Sun, and 150 G for 61 Vir (G6V). Evidence for an average magnetic field of ~ 500 G is found for 59 Vir (G0V), and of ~ 1000 G for HD 68456 (F6V). A distinction between magnetic and non-magnetic regions is made in TC models, while assuming a homogeneous distribution of both components. In our TC inversions of 59 Vir, we investigate three cases: both components have equal temperatures; warm magnetic regions; cool magnetic regions. Our TC model with equal temperatures does not yield significant improvement over OC inversions for 59 Vir. The resulting Bf values are consistent for both. Fit quality is significantly improved, however, by using two components of different temperatures. The inversions for 59 Vir that assume different temperatures for the two components yield results consistent with 0 - 450 G at the formal 3σ confidence level. We thus find a model dependence of our analysis and demonstrate that the influence of an additional temperature component can dominate over the Zeeman broadening signature, at least in optical data. Previous comparable analyses that neglected effects due to multiple temperature components may be prone to the same ambiguities.

preprint2010arXiv

On the kinematic age of brown dwarfs: Radial velocities and space motions of 43 nearby L dwarfs

We present radial velocity measurements of a sample of L0-L8 dwarfs observed with VLT/UVES and Keck/HIRES. We combine these measurements with distance and proper motion from the literature to determine space motions for 43 of our targets. We identify nine candidate members of young moving groups, which have ages of 50-600 Myr according to their space motion. From the total velocity dispersion of the 43 L dwarfs, we calculate a kinematic age of ~5 Gyr for our sample. This age is significantly higher than the ~3 Gyr age known for late M dwarfs in the solar neighbourhood. We find that the distributions of the U and V velocity components of our sample are clearly non-Gaussian, placing the age estimate inferred from the full space motion vector into question. The W-component exhibits a distribution more consistent with a normal distribution, and from W alone we derive an age of ~3 Gyr, which is the same age found for late-M dwarf samples. Our brightness-limited sample is probably contaminated by a number of outliers that predominantly bias the U and V velocity components. The origin of the outliers remain unclear, but we suggest that these brown dwarfs may have gained their high velocities by means of ejection from multiple systems during their formation.

preprint2010arXiv

Rotation and differential rotation in rapidly rotating field stars

We continue our studies on stellar latitudinal differential rotation. The presented work is a sequel of the work of Reiners et al. who studied the spectral line broadening profile of hundreds of stars of spectral types A through G at high rotational speed (vsini > 12 km/s). While most stars were found to be rigid rotators, only a few tens show the signatures of differential rotation. The present work comprises the rotational study of some 180 additional stars. The overall broadening profile is derived according to Reiners et al. from hundreds of spectral lines by least-squares deconvolution, reducing spectral noise to a minimum. Projected rotational velocities vsini are measured for about 120 of the sample stars. Differential rotation produces a cuspy line shape which is best measured in inverse wavelength space by the first two zeros of its Fourier transform. Rigid and differential rotation can be distinguished for more than 50 rapid rotators (vsini > 12 km/s) among the sample stars from the available spectra. Ten stars with significant differential rotation rates of 10-54 % are identified, which add to the few known rapid differential rotators. Differential rotation measurements of 6 % and less for four of our targets are probably spurious and below the detection limit. Including these objects, the line shapes of more than 40 stars are consistent with rigid rotation.

preprint2010arXiv

The CRIRES Search for Planets Around the Lowest-Mass Stars. I. High-Precision Near-Infrared Radial Velocities with an Ammonia Gas Cell

Radial velocities measured from near-infrared spectra are a potentially powerful tool to search for planets around cool stars and sub-stellar objects. However, no technique currently exists that yields near-infrared radial velocity precision comparable to that routinely obtained in the visible. We describe a method for measuring high-precision relative radial velocities of these stars from K-band spectra. The method makes use of a glass cell filled with ammonia gas to calibrate the spectrograph response similar to the "iodine cell" technique that has been used very successfully in the visible. Stellar spectra are obtained through the ammonia cell and modeled as the product of a Doppler-shifted template spectrum of the object and a spectrum of the cell, convolved with a variable instrumental profile model. A complicating factor is that a significant number of telluric absorption lines are present in the spectral regions containing useful stellar and ammonia lines. The telluric lines are modeled simultaneously as well using spectrum synthesis with a time-resolved model of the atmosphere over the observatory. The free parameters in the complete model are the wavelength scale of the spectrum, the instrumental profile, adjustments to the water and methane abundances in the atmospheric model, telluric spectrum Doppler shift, and stellar Doppler shift. Tests of the method based on the analysis of hundreds of spectra obtained for late M dwarfs over six months demonstrate that precisions of ~5 m/s are obtainable over long timescales, and precisions of better than 3 m/s can be obtained over timescales up to a week. The obtained precision is comparable to the predicted photon-limited errors, but primarily limited over long timescales by the imperfect modeling of the telluric lines.

preprint2010arXiv

The CRIRES Search for Planets Around the Lowest-Mass Stars. II. The Proposed Giant Planet Orbiting VB10 Does Not Exist

We present high-precision relative radial velocities of the very low-mass star VB10 that were obtained over a time span of 0.61 yr as part of an ongoing search for planets around stars at the end of the main sequence. The radial velocities were measured from high-resolution near-infrared spectra obtained using the CRIRES instrument on the VLT with an ammonia gas cell. The typical internal precision of the measurements is 10 m/s. These data do not exhibit significant variability and are essentially constant at a level consistent with the measurement uncertainties. Therefore, we do not detect the radial velocity variations of VB10 expected due to the presence of an orbiting giant planet similar to that recently proposed by Pravdo and Shaklan based on apparent astrometric perturbations. In addition, we do not confirm the ~1 km/s radial velocity variability of the star tentatively detected by Zapatero Osorio and colleagues with lower precision measurements. Our measurements rule out planets with M_p > 3 M_Jup and the orbital period and inclination suggested by Pravdo and Shaklan at better than 5 sigma confidence. We conclude that the planet detection claimed by Pravdo and Shaklan is spurious on the basis of this result. Although the outcome of this work is a non-detection, it illustrates the potential of using ammonia cell radial velocities to detect planets around very low-mass stars.

preprint2010arXiv

Zeeman Broadening in Optical Stokes I of Solar-like Stars

We present our analysis aimed at inferring average magnetic fields in slowly-rotating solar-like stars. Using the spectral line inversion code SPINOR, we perform high-accuracy line profile fitting and investigate whether Zeeman broadening can be reliably detected in optical data of unprecedented quality. We argue that our usage of both high- and low-g_{eff} lines does provide a certain sensitivity to magnetic fields that may, indeed, be detected. However, the measurement is subject to a model dependence and prone to ambiguities, e.g. due to spectral blends. Hence, while a field may be successfully recovered, the quantification of this field is subject to large uncertainties, even for the highest-quality optical data.

preprint2009arXiv

A Volume-limited Sample of 63 M7-M9.5 Dwarfs II. Activity, magnetism, and the fade of the rotation-dominated dynamo

In a volume-limited sample of 63 ultracool dwarfs of spectral type M7-M9.5, we have obtained high-resolution spectroscopy with UVES at the Very Large Telescope and HIRES at Keck Observatory. In this second paper, we present projected rotation velocities, average magnetic field strengths, and chromospheric emission from the Halpha line. We confirm earlier results that the mean level of normalized Halpha luminosity decreases with lower temperature, and we find that the scatter among Halpha luminosities is larger at lower temperature. We measure average magnetic fields between 0 and 4kG with no indication for a dependence on temperature between M7 and M9.5. For a given temperature, Halpha luminosity is related to magnetic field strength, consistent with results in earlier stars. A few very slowly rotating stars show very weak magnetic fields and Halpha emission, all stars rotating faster than our detection limit show magnetic fields of at least a few hundred Gauss. In contrast to earlier-type stars, we observe magnetic fields weaker than 1kG in stars rotating faster than ~3km/s, but we find no correlation between rotation and magnetic flux generation among them. We interpret this as a fundamental change in the dynamo mechanism; in ultracool dwarfs, magnetic field generation is predominantly achieved by a turbulent dynamo, while other mechanisms can operate more efficiently at earlier spectral type.

preprint2009arXiv

Activity-induced radial velocity jitter in a flaring M dwarf

We investigate the effect of stellar activity and flares on short-term radial velocity measurements in the mid-M flare star CN Leo. Radial velocity variations are calculated from 181 UVES spectra obtained during three nights. We searched for spectral orders that contain very few atmospheric absorption lines and calibrated them against the telluric A-band from O$_2$ in the Earth's atmosphere. One giant flare occurred during our observations, which has a very strong effect on radial velocity. The apparent radial velocity shift due to the flare is several hundred m s$^{-1}$ and clearly correlated with H$α$ emission. Outside the flare, only spectral orders containing the most prominent emission lines of H, He, and Ca show a correlation to chromospheric activity together with a radial velocity jitter exceeding a few 10 m s$^{-1}$. We identify a number of spectral orders that are free of strong emission lines and show no flaring-related radial velocity jitter, although flares occurred as strong as 0.4 dex in normalized H$α$ luminosity. The mean radial velocity jitter due to moderate flaring is less than 10 m s$^{-1}$. Strong flares are easily recognized directly in the spectra and should be neglected for planet searches.

preprint2009arXiv

Surprisingly Weak Magnetism on Young Accreting Brown Dwarfs

We have measured the surface magnetic flux on four accreting young brown dwarfs and one non-accreting young very low-mass star utilizing high resolution spectra of absorption lines of the FeH molecule. A magnetic field of 1-2 kG had been proposed for one of the brown dwarfs, 2MASS J1207334$-$393254, because of its similarities to higher mass T Tauri stars as manifested in accretion and the presence of a jet. We do not find clear evidence for a kilo-Gauss field in any of our young brown dwarfs but do find a 2 kG field on the young VLM star. Our 3-$σ$ upper limit for the magnetic flux in 2MASS J1207334$-$393254 just reaches 1 kG. We estimate the magnetic field required for accretion in young brown dwarfs given the observed rotations, and find that fields of only a few hundred Gauss are sufficient for magnetospheric accretion. This predicted value is less than our observed upper limit. We conclude that magnetic fields in young brown dwarfs are a factor of five or more lower than in young stars of about one solar mass, and in older stars with spectral types similar to our young brown dwarfs. It is interesting that, during the first few million years, the fields scale down with mass in line with what is needed for magnetospheric accretion, yet no such scaling is observed at later ages within the same effective temperature range. This scaling is opposite to the trend in rotation, with shorter rotation periods for very young accreting brown dwarfs compared with accreting solar-mass objects (and very low Rossby numbers in all cases). We speculate that in young objects a deeper intrinsic connection may exist between magnetospheric accretion and magnetic field strength, or that magnetic field generation in brown dwarfs may be less efficient than in stars. Neither of these currently have an easy physical explanation.

preprint2009arXiv

The Origin of Enhanced Activity in the Suns of M67

We report the results of the analysis of high resolution photospheric line spectra obtained with the UVES instrument on the VLT for a sample of 15 solar-type stars selected from a recent survey of the distribution of H and K chromospheric line strengths in the solar-age open cluster M67. We find upper limits to the projected rotation velocities that are consistent with solar-like rotation (i.e., v sini ~< 2-3 km/s) for objects with Ca II chromospheric activity within the range of the contemporary solar cycle. Two solar-type stars in our sample exhibit chromospheric emission well in excess of even solar maximum values. In one case, Sanders 1452, we measure a minimum rotational velocity of vsini = 4 +/- 0.5 km/s, or over twice the solar equatorial rotational velocity. The other star with enhanced activity, Sanders 747, is a spectroscopic binary. We conclude that high activity in solar-type stars in M67 that exceeds solar levels is likely due to more rapid rotation rather than an excursion in solar-like activity cycles to unusually high levels. We estimate an upper limit of 0.2% for the range of brightness changes occurring as a result of chromospheric activity in solar-type stars and, by inference, in the Sun itself. We discuss possible implications for our understanding of angular momentum evolution in solar-type stars, and we tentatively attribute the rapid rotation in Sanders 1452 to a reduced braking efficiency.

preprint2007arXiv

At the Bottom of the Main Sequence: Activity and magnetic fields beyond the threshold to complete convection

The bottom of the main sequence hosts objects with fundamentally different properties. At masses of about 0.3 M$_{\odot}$, stars become fully convective and at about 0.08 M$_{\odot}$ the hydrogen-burning main sequence ends; less massive objects are brown dwarfs. While stars and brown dwarfs experience very different evolutions, their inner structure has relatively little impact on the atmospheres. The generation of magnetic fields and activity is obviously connected to the threshold between partial and complete convection, because dynamo mechanisms involving a layer of shear like the solar $αΩ$-dynamo must cease. Hence a change in stellar activity can be expected there. Observations of stellar activity do not confirm a rapid break in activity at the convection boundary, but the fraction of active stars and rapid rotators is higher on the fully convective side. I summarize the current picture of stellar activity and magnetic field measurements at the bottom of the main sequence and present recent results on rotational braking beyond.

preprint2005arXiv

Rotation- and temperature-dependence of stellar latitudinal differential rotation

More than 600 high resolution spectra of stars with spectral type F and later were obtained in order to search for signatures of differential rotation in line profiles. In 147 stars, the rotation law could be measured, 28 of them are found to be differentially rotating. Comparison to rotation laws in stars of spectral type A reveals that differential rotation sets in at the convection boundary in the HR-diagram; no star that is significantly hotter than the convection boundary exhibits the signatures of differential rotation. Four late A-/early F-type stars close to the convection boundary and at vsini~100 km/s show extraordinarily strong absolute shear at short rotation periods around one day. It is suggested that this is due to their small convection zone depth and that it is connected to a narrow range in surface velocity. Detection frequencies of differential rotation were analyzed in stars with varying temperature and rotation velocity. Measurable differential rotation is more frequent in late-type stars and slow rotators. The strength of absolute shear and differential rotation are examined as functions of the stellar effective temperature and rotation period. The strongest shear is found at rotation periods between two and three days. In slower rotators, the strongest shear at a given rotation rate is given approximately by DOmega_max ~ P^{-1}. In faster rotators, alpha_max and DOmega_max diminish less rapidly. A comparison with differential rotation measurements in stars of later spectral type shows that F-stars exhibit stronger shear than cooler stars do, the upper boundary in absolute shear DOmega with temperature is consistent with the temperature scaling law found in Doppler Imaging measurements.