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.

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.

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.

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.

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.

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.

Time lags of the kilohertz quasi-periodic oscillations in the low-mass X-ray binaries 4U 1608−52 and 4U 1636−53

We studied the energy and frequency dependence of the Fourier time lags and intrinsic coherence of the kilohertz quasi-periodic oscillations (kHz QPOs) in the neutron-star lowmass X-ray binaries 4U 1608−52 and 4U 1636−53, using a large data set obtained with the Rossi X-ray Timing Explorer. We confirmed that, in both sources, the time lags of the lower kHz QPO are soft and their magnitude increases with energy. We also found that: (i) In 4U 1636−53, the soft lags of the lower kHz QPO remain constant at∼30 μs in the QPO frequency range 500–850 Hz, and decrease to ∼10 μs when the QPO frequency increases further. In 4U 1608−52, the soft lags of the lower kHz QPO remain constant at 40 μs up to 800 Hz, the highest frequency reached by this QPO in our data. (ii) In both sources, the time lags of the upper kHz QPO are hard, independent of energy or frequency and inconsistent with the soft lags of the lower kHz QPO. (iii) In both sources the intrinsic coherence of the lower kHz QPO remains constant at ∼0.6 between 5 and 12 keV, and drops to zero above that energy. The intrinsic coherence of the upper kHz QPO is consistent with being zero across the full energy range. (iv) In 4U 1636−53, the intrinsic coherence of the lower kHz QPO increases from ∼0 at ∼600 Hz to ∼1, and it decreases to ∼0.5 at 920 Hz; in 4U 1608−52, the intrinsic coherence is consistent with the same trend. (v) In both sources the intrinsic coherence of the upper kHz QPO is consistent with zero over the full frequency range of the QPO, except in 4U 1636−53 between 700 and 900 Hz where the intrinsic coherence marginally increases. We discuss our results in the context of scenarios in which the soft lags are either due to reflection off the accretion disc or up-/down-scattering in a hot medium close to the neutron star. We finally explore the connection between, on one hand the time lags and the intrinsic coherence of the kHz QPOs, and on the other the QPOs’ amplitude and quality factor in these two sources.

A cyclic bipolar wind in the interacting binary V 393 Scorpii

V393 Scorpii is a double periodic variable characterized by a relatively stable non-orbital photometric cycle of 253 d. Mennickent et al. argue for the presence of a massive optically thick disc around the more massive B-type component and describe the evolutionary stage of the system. In this paper, we analyse the behaviour of the main spectroscopic optical lines during the long non-orbital photometric cycle. We study the radial velocity of the donor determining its orbital elements and find a small but significant orbital eccentricity (e = 0.04). The donor spectral features are modelled and removed from the spectrum at every observing epoch using the light-curve model given by Mennickent et al. We find that the line emission is larger during eclipses and mostly comes from a bipolar wind. We also find that the long cycle is explained in terms of a modulation of the wind strength; the wind has a larger line and continuum emissivity at the high state. We report the discovery of highly variable chromospheric emission in the donor, as revealed by the Doppler maps of the emission lines Mg II 4481 and C I 6588. We discuss notable and some novel spectroscopic features like discrete absorption components, especially visible at blue depressed O I 7773 absorption wings during the second half-cycle, Balmer double emission with V/R curves showing 'Z-type' and 'S-type' excursions around secondary and main eclipses, respectively, and H beta emission wings extending up to +/- 2000 km s(-1). We also discuss possible causes for these phenomena and for their modulations with the long cycle.

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.

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.

MHD numerical simulations of colliding winds in massive binary systems - I. Thermal versus non-thermal radio emission

In the past few decades detailed observations of radio and X-ray emission from massive binary systems revealed a whole new physics present in such systems. Both thermal and non-thermal components of this emission indicate that most of the radiation at these bands originates in shocks. O and B-type stars and WolfRayet (WR) stars present supersonic and massive winds that, when colliding, emit largely due to the freefree radiation. The non-thermal radio and X-ray emissions are due to synchrotron and inverse Compton processes, respectively. In this case, magnetic fields are expected to play an important role in the emission distribution. In the past few years the modelling of the freefree and synchrotron emissions from massive binary systems have been based on purely hydrodynamical simulations, and ad hoc assumptions regarding the distribution of magnetic energy and the field geometry. In this work we provide the first full magnetohydrodynamic numerical simulations of windwind collision in massive binary systems. We study the freefree emission characterizing its dependence on the stellar and orbital parameters. We also study self-consistently the evolution of the magnetic field at the shock region, obtaining also the synchrotron energy distribution integrated along different lines of sight. We show that the magnetic field in the shocks is larger than that obtained when the proportionality between B and the plasma density is assumed. Also, we show that the role of the synchrotron emission relative to the total radio emission has been underestimated.

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.

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.

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.