Braking the Gas in the beta Pictoris Disk

(Abridged) The main sequence star beta Pictoris hosts the best studied circumstellar disk to date. Nonetheless, a long-standing puzzle has been around since the detection of metallic gas in the disk: radiation pressure from the star should blow the gas away, yet the observed motion is consistent with Keplerian rotation. In this work we search for braking mechanisms that can resolve this discrepancy. We find that all species affected by radiation force are heavily ionized and dynamically coupled into a single fluid by Coulomb collisions, reducing the radiation force on species feeling the strongest acceleration. For a gas of solar composition, the resulting total radiation force still exceeds gravity, while a gas of enhanced carbon abundance could be self-braking. We also explore two other braking agents: collisions with dust grains and neutral gas. Grains surrounding beta Pic are photoelectrically charged to a positive electrostatic potential. If a significant fraction of the grains are carbonaceous (10% in the midplane and larger at higher altitudes), ions can be slowed down to satisfy the observed velocity constraints. For neutral gas to brake the coupled ion fluid, we find the minimum required mass to be $\approx$ 0.03 M_\earth, consistent with observed upper limits of the hydrogen column density, and substantially reduced relative to previous estimates. Our results favor a scenario in which metallic gas is generated by grain evaporation in the disk, perhaps during grain-grain collisions. We exclude a primordial origin for the gas, but cannot rule out the possibility of its production by falling evaporating bodies near the star. We discuss the implications of this work for observations of gas in other debris disks.Comment: 19 pages, 12 figures, emulateapj. Accepted for publication in Ap

Comparison of different exoplanet mass detection limit methods using a sample of main-sequence intermediate-type stars

The radial velocity (RV) technique is a powerful tool for detecting extrasolar planets and deriving mass detection limits that are useful for constraining planet pulsations and formation models. Detection limit methods must take into account the temporal distribution of power of various origins in the stellar signal. These methods must also be able to be applied to large samples of stellar RV time series We describe new methods for providing detection limits. We compute the detection limits for a sample of ten main sequence stars, which are of G-F-A type, in general active, and/or with detected planets, and various properties. We use them to compare the performances of these methods with those of two other methods used in the litterature. We obtained detection limits in the 2-1000 day period range for ten stars. Two of the proposed methods, based on the correlation between periodograms and the power in the periodogram of the RV time series in specific period ranges, are robust and represent a significant improvement compared to a method based on the root mean square of the RV signal. We conclude that two of the new methods (correlation-based method and local power analysis, i.e. LPA, method) provide robust detection limits, which are better than those provided by methods that do not take into account the temporal sampling.Comment: 18 pages, 15 figures Accepted in Astronomy & Astrophysic

A new method of correcting radial velocity time series for inhomogeneous convection

Magnetic activity strongly impacts stellar RVs and the search for small planets. We showed previously that in the solar case it induces RV variations with an amplitude over the cycle on the order of 8 m/s, with signals on short and long timescales. The major component is the inhibition of the convective blueshift due to plages. We explore a new approach to correct for this major component of stellar radial velocities in the case of solar-type stars. The convective blueshift depends on line depths; we use this property to develop a method that will characterize the amplitude of this effect and to correct for this RV component. We build realistic RV time series corresponding to RVs computed using different sets of lines, including lines in different depth ranges. We characterize the performance of the method used to reconstruct the signal without the convective component and the detection limits derived from the residuals. We identified a set of lines which, combined with a global set of lines, allows us to reconstruct the convective component with a good precision and to correct for it. For the full temporal sampling, the power in the range 100-500~d significantly decreased, by a factor of 100 for a RV noise below 30 cm/s. We also studied the impact of noise contributions other than the photon noise, which lead to uncertainties on the RV computation, as well as the impact of the temporal sampling. We found that these other sources of noise do not greatly alter the quality of the correction, although they need a better noise level to reach a similar performance level. A very good correction of the convective component can be achieved providing very good RV noise levels combined with a very good instrumental stability and realistic granulation noise. Under the conditions considered in this paper, detection limits at 480~d lower than 1 MEarth could be achieved for RV noise below 15 cm/s.Comment: Accepted in A&A 18 July 201

Interpreting the photometry and spectroscopy of directly imaged planets: a new atmospheric model applied to beta Pictoris b and SPHERE observations

We aim to interpret future photometric and spectral measurements from these instruments, in terms of physical parameters of the planets, with an atmospheric model using a minimal number of assumptions and parameters. We developed Exoplanet Radiative-convective Equilibrium Model (Exo-REM) to analyze the photometric and spectro- scopic data of directly imaged planets. The input parameters are a planet's surface gravity (g), effective temperature (Teff ), and elemental composition. The model predicts the equilibrium temperature profile and mixing ratio profiles of the most important gases. Opacity sources include the H2-He collision-induced absorption and molecular lines from eight compounds (including CH4 updated with the Exomol line list). Absorption by iron and silicate cloud particles is added above the expected condensation levels with a fixed scale height and a given optical depth at some reference wavelength. Scattering was not included at this stage. We applied Exo-REM to photometric and spectral observations of the planet beta Pictoris b obtained in a series of near-IR filters. We derived Teff = 1550 +- 150 K, log(g) = 3.5 +- 1, and radius R = 1.76 +- 0.24 RJup (2-{\sigma} error bars from photometric measurements). These values are comparable to those found in the literature, although with more conservative error bars, consistent with the model accuracy. We were able to reproduce, within error bars, the J- and H-band spectra of beta Pictoris b. We finally investigated the precision to which the above parameterComment: 15 pages, 14 figures, accepted by A&

Gas Absorption Detected from the Edge-on Debris Disk Surrounding HD32297

Near-infrared and optical imaging of HD32297 indicate that it has an edge-on debris disk, similar to beta Pic. I present high resolution optical spectra of the NaI doublet toward HD32297 and stars in close angular proximity. A circumstellar absorption component is clearly observed toward HD32297 at the stellar radial velocity, which is not observed toward any of its neighbors, including the nearest only 0.9 arcmin away. An interstellar component is detected in all stars >90 pc, including HD32297, likely due to the interstellar material at the boundary of the Local Bubble. Radial velocity measurements of the nearest neighbors, BD+07 777s and BD+07 778, indicate that they are unlikely to be physically associated with HD32297. The measured circumstellar column density around HD32997, log N(NaI) ~ 11.4, is the strongest NaI absorption measured toward any nearby main sequence debris disk, even the prototypical edge-on debris disk, beta Pic. Assuming that the morphology and abundances of the gas component around HD32297 are similar to beta Pic, I estimate an upper limit to the gas mass in the circumstellar disk surrounding HD32297 of ~0.3 M_Earth.Comment: 13 pages, 2 figures; Accepted for publication in ApJ Letter

Variability of stellar granulation and convective blueshift with spectral type and magnetic activity. I. K and G main sequence stars

In solar-type stars, the attenuation of convective blueshift by stellar magnetic activity dominates the RV variations over the low amplitude signal induced by low mass planets. Models of stars that differ from the Sun will require a good knowledge of the attenuation of the convective blueshift to estimate its impact on the variations. It is therefore crucial to precisely determine not only the amplitude of the convective blueshift for different types of stars, but also the dependence of this convective blueshift on magnetic activity, as these are key factors in our model producing the RV. We studied a sample of main sequence stars with spectral types from G0 to K2 and focused on their temporally averaged properties: the activity level and a criterion allowing to characterise the amplitude of the convective blueshift. We find the differential velocity shifts of spectral lines due to convection to depend on the spectral type, the wavelength (this dependence is correlated with the Teff and activity level), and on the activity level. This allows us to quantify the dependence of granulation properties on magnetic activity for stars other than the Sun. The attenuation factor of the convective blueshift appears to be constant over the considered range of spectral types. We derive a convective blueshift which decreases towards lower temperatures, with a trend in close agreement with models for Teff lower than 5800 K, but with a significantly larger global amplitude. We finally compare the observed RV variation amplitudes with those that could be derived from our convective blueshift using a simple law and find a general agreement on the amplitude. Our results are consistent with previous results and provide, for the first time, an estimation of the convective blueshift as a function of Teff, magnetic activity, and wavelength, over a large sample of G and K main sequence stars

Reconstructing the solar integrated radial velocity using MDI/SOHO

Searches for exoplanets with radial velocity techniques are increasingly sensitive to stellar activity. It is therefore crucial to characterize how this activity influences radial velocity measurements in their study of the detectability of planets in these conditions. In a previous work we simulated the impact of spots and plages on the radial velocity of the Sun. Our objective is to compare this simulation with the observed radial velocity of the Sun for the same period. We use Dopplergrams and magnetograms obtained by MDI/SOHO over one solar cycle to reconstruct the solar integrated radial velocity in the Ni line 6768 \AA. We also characterize the relation between the velocity and the local magnetic field to interpret our results. We obtain a stronger redshift in places where the local magnetic field is larger (and as a consequence for larger magnetic structures): hence we find a higher attenuation of the convective blueshift in plages than in the network. Our results are compatible with an attenuation of this blueshift by about 50% when averaged over plages and network. We obtain an integrated radial velocity with an amplitude over the solar cycle of about 8 m/s, with small-scale variations similar to the results of the simulation, once they are scaled to the Ni line. The observed solar integrated radial velocity agrees with the result of the simulation made in our previous work within 30%, which validates this simulation. The observed amplitude confirms that the impact of the convective blueshift attenuation in magnetic regions will be critical to detect Earth-mass planets in the habitable zone around solar-like stars.Comment: 17 pages, 11 figures, accepted in Astronomy and Astrophysic

Sparse aperture masking at the VLT II. Detection limits for the eight debris disks stars β\beta Pic, AU Mic, 49 Cet, η\eta Tel, Fomalhaut, g Lup, HD181327 and HR8799

Context. The formation of planetary systems is a common, yet complex mechanism. Numerous stars have been identified to possess a debris disk, a proto-planetary disk or a planetary system. The understanding of such formation process requires the study of debris disks. These targets are substantial and particularly suitable for optical and infrared observations. Sparse Aperture masking (SAM) is a high angular resolution technique strongly contributing to probe the region from 30 to 200 mas around the stars. This area is usually unreachable with classical imaging, and the technique also remains highly competitive compared to vortex coronagraphy. Aims. We aim to study debris disks with aperture masking to probe the close environment of the stars. Our goal is either to find low mass companions, or to set detection limits. Methods. We observed eight stars presenting debris disks ( $\beta$ Pictoris, AU Microscopii, 49 Ceti, $\eta$ Telescopii, Fomalhaut, g Lupi, HD181327 and HR8799) with SAM technique on the NaCo instrument at the VLT. Results. No close companions were detected using closure phase information under 0.5 of separation from the parent stars. We obtained magnitude detection limits that we converted to Jupiter masses detection limits using theoretical isochrones from evolutionary models. Conclusions. We derived upper mass limits on the presence of companions in the area of few times the diffraction limit of the telescope around each target star.Comment: 7 pages, All magnitude detection limits maps are only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5

First NACO observations of the Brown Dwarf LHS 2397aB

Observations of the standard late type M8 star LHS 2397aA were obtained at the ESO-VLT 8m telescope ``Yepun'' using the NAOS/CONICA Adaptive Optics facility. The observations were taken during the NACO commissioning, and the infrared standard star LHS 2397aA was observed in the H, and Ks broad band filters. In both bands the brown dwarf companion LHS2397aB was detected. Using a program recently developed (Bouy et al., 2003) for the detection of stellar binaries we calculated the principal astrometric parameters (angular binary separation and position angle P.A.) and the photometry of LHS 2397aA and LHS 2397aB. Our study largely confirms previous results obtained with the AO-Hokupa'a facility at Gemini-North (Freed et al., 2003); however a few discrepancies are observed.Comment: 5 page