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

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

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

Using the Sun to estimate Earth-like planets detection capabilities.I. Impact of cold spots

Stellar spots may in some cases produce radial velocity (RV) signatures similar to those of exoplanets. To further investigate the impact of spots, we aim at studying the detectability of Earth mass planets in the habitable zone (HZ) of solar type stars, if covered by spots similar to the sunspots. We have used the Sunspots properties recorded over one solar cycle between 1993 and 2003 to build the RV curve that a solar type star seen edge-on would show, if covered by such spots with Tsun -Tspot = 550K. We also simulate the RV of such a spotted star surrounded by an Earth mass planet located in the HZ. Under present assumptions, the detection of a 1 M Earth planet located between 0.8 and 1.2 AU requires an intensive monitoring (weekly or better), during several years of low activity phasis. The temporal sampling is more crucial than the precision of the data (assuming precisions in the range [1-10] cm/s). Cooler spots may become a problem for such detections. Also, we anticipate that plages, not considered in this paper, could further complicate or even compromise the detections

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

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

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

Multiple spiral patterns in the transitional disk of HD 100546

Protoplanetary disks around young stars harbor many structures related to planetary formation. Of particular interest, spiral patterns were discovered among several of these disks and are expected to be the sign of gravitational instabilities leading to giant planets formation or gravitational perturbations caused by already existing planets. In this context, the star HD100546 presents some specific characteristics with a complex gas and dusty disk including spirals as well as a possible planet in formation. The objective of this study is to analyze high contrast and high angular resolution images of this emblematic system to shed light on critical steps of the planet formation. We retrieved archival images obtained at Gemini in the near IR (Ks band) with the instrument NICI and processed the data using advanced high contrast imaging technique taking advantage of the angular differential imaging. These new images reveal the spiral pattern previously identified with HST with an unprecedented resolution, while the large-scale structure of the disk is mostly erased by the data processing. The single pattern at the southeast in HST images is now resolved into a multi-armed spiral pattern. Using two models of a gravitational perturber orbiting in a gaseous disk we attempted to bring constraints on the characteristics of this perturber assuming each spiral being independent and we derived qualitative conclusions. The non-detection of the northeast spiral pattern observed in HST allows to put a lower limit on the intensity ratio between the two sides of the disk, which if interpreted as forward scattering yields a larger anisotropic scattering than derived in the visible. Also, we found that the spirals are likely spatially resolved with a thickness of about 5-10AU. Finally, we did not detect the candidate forming planet recently discovered in the Lp band, with a mass upper limit of 16-18 MJ.Comment: Accepted for publication in Astronomy and Astrophysics, 10 pages, 8 figure