The Gemini NICI planet-finding campaign: the frequency of giant planets around young B and A stars

We have carried out high contrast imaging of 70 young, nearby B and A stars to search for brown dwarf and planetary companions as part of the Gemini NICI Planet-Finding Campaign. Our survey represents the largest, deepest survey for planets around high-mass stars (≈1.5-2.5 M ☉) conducted to date and includes the planet hosts β Pic and Fomalhaut. We obtained follow-up astrometry of all candidate companions within 400 AU projected separation for stars in uncrowded fields and identified new low-mass companions to HD 1160 and HIP 79797. We have found that the previously known young brown dwarf companion to HIP 79797 is itself a tight (3 AU) binary, composed of brown dwarfs with masses 58$^{+21}_{-20}$ M Jup and 55$^{+20}_{-19}$ M Jup, making this system one of the rare substellar binaries in orbit around a star. Considering the contrast limits of our NICI data and the fact that we did not detect any planets, we use high-fidelity Monte Carlo simulations to show that fewer than 20% of 2 M ☉ stars can have giant planets greater than 4 M Jup between 59 and 460 AU at 95% confidence, and fewer than 10% of these stars can have a planet more massive than 10 M Jup between 38 and 650 AU. Overall, we find that large-separation giant planets are not common around B and A stars: fewer than 10% of B and A stars can have an analog to the HR 8799 b (7 M Jup, 68 AU) planet at 95% confidence. We also describe a new Bayesian technique for determining the ages of field B and A stars from photometry and theoretical isochrones. Our method produces more plausible ages for high-mass stars than previous age-dating techniques, which tend to underestimate stellar ages and their uncertainties.

Probing the properties of Be star discs with spectroastrometry and NLTE radiative transfer modelling: beta CMi

While the presence of discs around classical Be stars is well established, their origin is still uncertain. To understand what processes result in the creation of these discs and how angular momentum is transported within them, their physical properties must be constrained. This requires comparing high spatial and spectral resolution data with detailed radiative transfer modelling. We present a high spectral resolution, R similar to 80 000, sub-milliarcsecond precision, spectroastrometric study of the circumstellar disc around the Be star beta CMi. The data are confronted with 3D, non-local thermodynamic equilibrium radiative transfer calculations to directly constrain the properties of the disc. Furthermore, we compare the data to disc models featuring two velocity laws: Keplerian, the prediction of the viscous disc model, and angular momentum conserving rotation. It is shown that the observations of beta CMi can only be reproduced using Keplerian rotation. The agreement between the model and the observed spectral energy distribution, polarization and spectroastrometric signature of beta CMi confirms that the discs around Be stars are well modelled as viscous decretion discs.

Alessi 95 and the short-period Cepheid SU Cassiopeiae

The parameters for the newly discovered open cluster Alessi 95 are established on the basis of available photometric and spectroscopic data, in conjunction with new observations. Colour excesses for spectroscopically observed B- and A-type stars near SU Cas follow a reddening relation described by E(U-B)/E(B-V) = 0.83 + 0.02E(B-V), implying a value of R=AV/E(B-V) ? 2.8 for the associated dust. Alessi 95 has a mean reddening of E(B-V)(B0) = 0.35 +/- 0.02 s.e., an intrinsic distance modulus of V0-MV= 8.16 +/- 0.04 s.e. (+/- 0.21 s.d.), d= 429 +/- 8 pc, and an estimated age of 108.2 yr from zero-age main sequence (ZAMS) fitting of available UBV, CCD BV, NOMAD, and Two Micron All Sky Survey JHKs observations of cluster stars. SU Cas is a likely cluster member, with an inferred space reddening of E(B-V) = 0.33 +/- 0.02 and a luminosity of < MV >=-3.15 +/- 0.07 s.e., consistent with overtone pulsation (PFM= 2.75 d), as also implied by the Cepheids light-curve parameters, rate of period increase and Hipparcos parallaxes for cluster stars. There is excellent agreement of the distance estimates for SU Cas inferred from cluster ZAMS fitting, its pulsation parallax derived from the infrared surface brightness technique and Hipparcos parallaxes, which all agree to within a few per cent.

Polarimetry of the binary PDS 144

Aims. Our goal is to study the circumstellar environment associated with each component of the wide intermediate-mass pre-main sequence binary system PDS 144 using broadband polarimetry. Methods. We present near-infrared (NIR) linear polarimetric observations of PDS 144 gathered with the IAGPOL imaging polarimeter along with the CamIV infrared camera at the Observatorio do Pico dos Dias (OPD). In addition, we re-analyzed OPD archive optical polarization to separate the binary and estimate the interstellar polarization using foreground stars. Results. After discounting the interstellar component, we found that both stars of the binary system are intrinsically polarized. The polarization vectors at optical and NIR bands of both components are aligned with the local magnetic field and the jet axis. These findings indicate an interplay between the interstellar magnetic field and the formation of the binary system. We also found that the PDS 144N is less polarized than its southern companion in the optical. However, in the NIR PDS 144N is more polarized. Our polarization data can only be explained by high inclinations (i greater than or similar to 80 degrees) for the disks of both members. In particular, comparisons of our NIR data with young stellar objects disk models suggest predominantly small grains in the circumstellar environment of PDS 144N. In spite of the different grain types in each component, the infrared spectral indexes indicate a coeval system. We also found evidence of coplanarity between the disks.

He II lambda 4686 IN eta CARINAE: COLLAPSE OF THE WIND-WIND COLLISION REGION DURING PERIASTRON PASSAGE

The periodic spectroscopic events in eta Carinae are now well established and occur near the periastron passage of two massive stars in a very eccentric orbit. Several mechanisms have been proposed to explain the variations of different spectral features, such as an eclipse by the wind-wind collision (WWC) boundary, a shell ejection from the primary star or accretion of its wind onto the secondary. All of them have problems explaining all the observed phenomena. To better understand the nature of the cyclic events, we performed a dense monitoring of eta Carinae with five Southern telescopes during the 2009 low-excitation event, resulting in a set of data of unprecedented quality and sampling. The intrinsic luminosity of the He II lambda 4686 emission line (L similar to 310 L-circle dot) just before periastron reveals the presence of a very luminous transient source of extreme UV radiation emitted in the WWC region. Clumps in the primary's wind probably explain the flare-like behavior of both the X-ray and He II lambda 4686 light curves. After a short-lived minimum, He II lambda 4686 emission rises again to a new maximum, when X-rays are still absent or very weak. We interpret this as a collapse of the WWC onto the "surface" of the secondary star, switching off the hard X-ray source and diminishing the WWC shock cone. The recovery from this state is controlled by the momentum balance between the secondary's wind and the clumps in the primary's wind.

Origin of the heavy elements in HD 140283 Measurement of europium abundance

Context. HD 140283 is a nearby (V = 7.7) subgiant metal-poor star, extensively analysed in the literature. Although many spectra have been obtained for this star, none showed a signal-to-noise (S/N) ratio high enough to enable a very accurate derivation of abundances from weak lines. Aims. The detection of europium proves that the neutron-capture elements in this star originate in the r-process, and not in the s-process, as recently claimed in the literature. Methods. Based on the OSMARCS 1D LTE atmospheric model and with a consistent approach based on the spectrum synthesis code Turbospectrum, we measured the europium lines at 4129 angstrom and 4205 angstrom, taking into account the hyperfine structure of the transitions. The spectrum, obtained with a long exposure time of seven hours at the Canada-France-Hawaii Telescope (CFHT), has a resolving power of 81 000 and a S/N ratio of 800 at 4100 angstrom. Results. We were able to determine the abundance A(Eu) = -2.35 +/- 0.07 dex, compatible with the value predicted for the europium from the r-process. The abundance ratio [Eu/Ba] = +0.58 +/- 0.15 dex agrees with the trend observed in metal-poor stars and is also compatible with a strong r-process contribution to the origin of the neutron-capture elements in HD 140283.

THE FIRST DETERMINATION OF THE VISCOSITY PARAMETER IN THE CIRCUMSTELLAR DISK OF A Be STAR

Be stars possess gaseous circumstellar decretion disks, which are well described using standard alpha-disk theory. The Be star 28 CMa recently underwent a long outburst followed by a long period of quiescence, during which the disk dissipated. Here we present the first time-dependent models of the dissipation of a viscous decretion disk. By modeling the rate of decline of the V-band excess, we determine that the viscosity parameter alpha = 1.0 +/- 0.2, corresponding to a mass injection rate (M) over dot = (3.5 +/- 1.3) x 10(-8) M-circle dot yr(-1). Such a large value of a suggests that the origin of the turbulent viscosity is an instability in the disk whose growth is limited by shock dissipation. The mass injection rate is more than an order of magnitude larger than the wind mass-loss rate inferred from UV observations, implying that the mass injection mechanism most likely is not the stellar wind, but some other mechanism.

Stochastic gravito-inertial modes discovered by CoRoT in the hot Be star HD51452

Context. Be stars are rapidly rotating stars with a circumstellar decretion disk. They usually undergo pressure and/or gravity pulsation modes excited by the kappa-mechanism, i.e. an effect of the opacity of iron-peak elements in the envelope of the star. In the Milky Way, p-modes are observed in stars that are hotter than or equal to the B3 spectral type, while g-modes are observed at the B2 spectral type and cooler. Aims. We observed a B0IVe star, HD51452, with the high-precision, high-cadence photometric CoRoT satellite and high-resolution, ground-based HARPS and SOPHIE spectrographs to study its pulsations in great detail. We also used the lower resolution spectra available in the BeSS database. Methods. We analyzed the CoRoT and spectroscopic data with several methods: CLEAN-NG, FREQFIND, and a sliding window method. We also analyzed spectral quantities, such as the violet over red (V/R) emission variations, to obtain information about the variation in the circumstellar environment. We calculated a stellar structure model with the ESTER code to test the various interpretation of the results. Results. We detect 189 frequencies of variations in the CoRoT light curve in the range between 0 and 4.5 c d(-1). The main frequencies are also recovered in the spectroscopic data. In particular we find that HD51452 undergoes gravito-inertial modes that are not in the domain of those excited by the kappa-mechanism. We propose that these are stochastic modes excited in the convective zones and that at least some of them are a multiplet of r-modes (i.e. subinertial modes mainly driven by the Coriolis acceleration). Stochastically excited gravito-inertial modes had never been observed in any star, and theory predicted that their very low amplitudes would be undetectable even with CoRoT. We suggest that the amplitudes are enhanced in HD51452 because of the very rapid stellar rotation. In addition, we find that the amplitude variations of these modes are related to the occurrence of minor outbursts. Conclusions. Thanks to CoRoT data, we have detected a new kind of pulsations in HD51452, which are stochastically excited gravito-inertial modes, probably due to its very rapid rotation. These modes are probably also present in other rapidly rotating hot Be stars.

The relationship between gamma Cassiopeiae's X-ray emission and its circumstellar environment

gamma Cas is the prototypical classical Be star and is recently best known for its variable hard X-ray emission. To elucidate the reasons for this emission, we mounted a multiwavelength campaign in 2010 centered around four XMM-Newton observations. The observational techniques included long baseline optical interferometry (LBOI) from two instruments at CHARA, photometry carried out by an automated photometric telescope and H alpha observations. Because gamma Cas is also known to be in a binary, we measured radial velocities from the H alpha line and redetermined its period as 203.55 +/- 0.20 days and its eccentricity as near zero. The LBOI observations suggest that the star's decretion disk was axisymmetric in 2010, has an system inclination angle near 45 degrees, and a larger radius than previously reported. In addition, the Be star began an "outburst" at the beginning of our campaign, made visible by a brightening and reddening of the disk during our campaign and beyond. Our analyses of the new high resolution spectra disclosed many attributes also found from spectra obtained in 2001 (Chandra) and 2004 (XMM-Newton). As well as a dominant hot (approximate to 14 keV) thermal component, the familiar attributes included: (i) a fluorescent feature of Fe K even stronger than observed at previous times; (ii) strong lines of N VII and Ne XI lines indicative of overabundances; and (iii) a subsolar Fe abundance from K-shell lines but a solar abundance from L-shell ions. We also found that two absorption columns are required to fit the continuum. While the first one maintained its historical average of 1 x 10(21) cm(-2), the second was very large and doubled to 7.4 x 10(23) cm(-2) during our X-ray observations. Although we found no clear relation between this column density and orbital phase, it correlates well with the disk brightening and reddening both in the 2010 and earlier observations. Thus, the inference from this study is that much (perhaps all?) of the X-ray emission from this source originates behind matter ejected by gamma Cas into our line of sight.

Beyond the diffraction limit of optical/IR interferometers I. Angular diameter and rotation parameters of Achernar from differential phases

Context. Spectrally resolved long-baseline optical/IR interferometry of rotating stars opens perspectives to investigate their fundamental parameters and the physical mechanisms that govern their interior, photosphere, and circumstellar envelope structures. Aims. Based on the signatures of stellar rotation on observed interferometric wavelength-differential phases, we aim to measure angular diameters, rotation velocities, and orientation of stellar rotation axes. Methods. We used the AMBER focal instrument at ESO-VLTI in its high-spectral resolution mode to record interferometric data on the fast rotator Achernar. Differential phases centered on the hydrogen Br gamma line (K band) were obtained during four almost consecutive nights with a continuous Earth-rotation synthesis during similar to 5h/night, corresponding to similar to 60 degrees position angle coverage per baseline. These observations were interpreted with our numerical code dedicated to long-baseline interferometry of rotating stars. Results. By fitting our model to Achernar's differential phases from AMBER, we could measure its equatorial radius R-eq = 11.6 +/- 0.3 R-circle dot, equatorial rotation velocity V-eq = 298 +/- 9 km s(-1), rotation axis inclination angle i = 101.5 +/- 5.2 degrees, and rotation axis position angle (from North to East) PA(rot) = 34.9 +/- 1.6 degrees. From these parameters and the stellar distance, the equatorial angular diameter circle divide(eq) of Achernar is found to be 2.45 +/- 0.09 mas, which is compatible with previous values derived from the commonly used visibility amplitude. In particular, circle divide(eq) and PA(rot) measured in this work with VLTI/AMBER are compatible with the values previously obtained with VLTI/VINCI. Conclusions. The present paper, based on real data, demonstrates the super-resolution potential of differential interferometry for measuring sizes, rotation velocities, and orientation of rotating stars in cases where visibility amplitudes are unavailable and/or when the star is partially or poorly resolved. In particular, we showed that differential phases allow the measurement of sizes up to similar to 4 times smaller than the diffraction-limited angular resolution of the interferometer.

THE GEMINI NICI PLANET-FINDING CAMPAIGN: DISCOVERY OF A MULTIPLE SYSTEM ORBITING THE YOUNG A STAR HD 1160

We report the discovery of two low-mass companions to the young A0V star HD 1160 at projected separations of 81 +/- 5 AU (HD 1160 B) and 533 +/- 25 AU (HD 1160 C) by the Gemini NICI Planet-Finding Campaign. Very Large Telescope images of the system taken over a decade for the purpose of using HD 1160 A as a photometric calibrator confirm that both companions are physically associated. By comparing the system to members of young moving groups and open clusters with well-established ages, we estimate an age of 50(-40)(+50) Myr for HD 1160 ABC. While the UVW motion of the system does not match any known moving group, the small magnitude of the space velocity is consistent with youth. Near-IR spectroscopy shows HD 1160 C to be an M3.5 +/- 0.5 star with an estimated mass of 0.22(-0.04)(+0.03) M-circle dot, while NIR photometry of HD 1160 B suggests a brown dwarf with a mass of 33(-9)(+12) M-Jup. The very small mass ratio (0.014) between the A and B components of the system is rare for A star binaries, and would represent a planetary-mass companion were HD 1160 A to be slightly less massive than the Sun.

Variable rotational line broadening in the Be star Achernar

Aims. The main theoretical problem for the formation of a Keplerian disk around Be stars is how angular momentum is supplied from the star to the disk, even more so since Be stars probably rotate somewhat subcritically. For instance, nonradial pulsation may transport angular momentum to the stellar surface until (part of) this excess supports the disk-formation/replenishment. The nearby Be star Achernar is presently building a new disk and o ers an excellent opportunity to observe this process from relatively close-up. Methods. Spectra from various sources and epochs are scrutinized to identify the salient stellar parameters characterizing the disk life cycle as defined by H emission. The variable strength of the non-radial pulsation is confirmed, but does not a ect the other results. Results. For the first time it is demonstrated that the photospheric line width does vary in a Be star, by as much as v sin i . 35 km However, unlike assumptions in which a photospheric spin-up accumulates during the diskless phase and then is released into the disk as it is fed, the apparent photospheric spin-up is positively correlated with the appearance of H line emission. The photospheric line widths and circumstellar emission increase together, and the apparent stellar rotation declines to the value at quiescence after the H line emission becomes undetectable

Origin of the heavy elements in HD 140283: measurement of europium abundance

Context. HD140283 is a nearby (V = 7:7) subgiant metal-poor star, extensively analysed in the literature. Although many spectra have been obtained for this star, none showed a signal-to-noise (S/N) ratio high enough to enable a very accurate derivation of abundances from weak lines. Aims. The detection of europium proves that the neutron-capture elements in this star originate in the r-process, and not in the s-process, as recently claimed in the literature. Methods. Based on the OSMARCS 1D LTE atmospheric model and with a consistent approach based on the spectrum synthesis code Turbospectrum, we measured the europium lines at 4129 Å and 4205 Å, taking into account the hyperfine structure of the transitions. The spectrum, obtained with a long exposure time of seven hours at the Canada-France-Hawaii Telescope (CFHT), has a resolving power of 81 000 and a S/N ratio of 800 at 4100 Å. Results. We were able to determine the abundance A(Eu) =

Stochastic gravito-inertial modes discovered by CoRoT in the hot Be star HD51452

Context. Be stars are rapidly rotating stars with a circumstellar decretion disk. They usually undergo pressure and/or gravity pulsation modes excited by the κ-mechanism, i.e. an effect of the opacity of iron-peak elements in the envelope of the star. In the Milky Way, p-modes are observed in stars that are hotter than or equal to the B3 spectral type, while g-modes are observed at the B2 spectral type and cooler. Aims. We observed a B0IVe star, HD51452, with the high-precision, high-cadence photometric CoRoT satellite and high-resolution, ground-based HARPS and SOPHIE spectrographs to study its pulsations in great detail. We also used the lower resolution spectra available in the BeSS database. Methods. We analyzed the CoRoT and spectroscopic data with several methods: Clean-NG, FreqFind, and a sliding window method. We also analyzed spectral quantities, such as the violet over red (V/R) emission variations, to obtain information about the variation in the circumstellar environment. We calculated a stellar structure model with the ESTER code to test the various interpretation of the results. Results. We detect 189 frequencies of variations in the CoRoT light curve in the range between 0 and 4.5 c d−1. The main frequencies are also recovered in the spectroscopic data. In particular we find that HD51452 undergoes gravito-inertial modes that are not in the domain of those excited by the κ-mechanism. We propose that these are stochastic modes excited in the convective zones and that at least some of them are a multiplet of r-modes (i.e. subinertial modes mainly driven by the Coriolis acceleration). Stochastically excited gravito-inertial modes had never been observed in any star, and theory predicted that their very low amplitudes would be undetectable even with CoRoT. We suggest that the amplitudes are enhanced in HD51452 because of the very rapid stellar rotation. In addition, we find that the amplitude variations of these modes are related to the occurrence of minor outbursts. Conclusions. Thanks to CoRoT data, we have detected a new kind of pulsations in HD51452, which are stochastically excited gravito-inertial modes, probably due to its very rapid rotation. These modes are probably also present in other rapidly rotating hot Be stars.

OXYGEN ABUNDANCES IN LOW- AND HIGH-alpha FIELD HALO STARS AND THE DISCOVERY OF TWO FIELD STARS BORN IN GLOBULAR CLUSTERS

Oxygen abundances of 67 dwarf stars in the metallicity range -1.6 < [Fe/H] < -0.4 are derived from a non-LTE analysis of the 777 nm O I triplet lines. These stars have precise atmospheric parameters measured by Nissen and Schuster, who find that they separate into three groups based on their kinematics and alpha-element (Mg, Si, Ca, Ti) abundances: thick disk, high-alpha halo, and low-alpha halo. We find the oxygen abundance trends of thick-disk and high-alpha halo stars very similar. The low-alpha stars show a larger star-to-star scatter in [O/Fe] at a given [Fe/H] and have systematically lower oxygen abundances compared to the other two groups. Thus, we find the behavior of oxygen abundances in these groups of stars similar to that of the a elements. We use previously published oxygen abundance data of disk and very metal-poor halo stars to present an overall view (-2.3 < [Fe/H] < +0.3) of oxygen abundance trends of stars in the solar neighborhood. Two field halo dwarf stars stand out in their O and Na abundances. Both G53-41 and G150-40 have very low oxygen and very high sodium abundances, which are key signatures of the abundance anomalies observed in globular cluster (GC) stars. Therefore, they are likely field halo stars born in GCs. If true, we estimate that at least 3% +/- 2% of the local field metal-poor star population was born in GCs.