Star Zeta Research Articles

Time-series photometry of the O4 I(n)fp star ζ Puppis

We report a time-series analysis of the O4 I(n)fp star zeta Pup, based on optical photometry obtained with the SMEI instrument on the Coriolis satellite, 2003--2006. A single astrophysical signal is found, with P = (1.780938 \pm 0.000093) d and a mean semi-amplitude of (6.9 \pm 0.3) mmag. There is no evidence for persistent coherent signals with semi-amplitudes in excess of ca. 2~mmag on any of the timescales previously reported in the literature. In particular, there is no evidence for a signature of the proposed rotation period, ca. 5.1~days; zeta Pup is therefore probably not an oblique magnetic rotator. The 1.8-day signal varies in amplitude by a factor ca. 2 on timescales of 10--100d (and probably by more on longer timescales), and exhibits modest excursions in phase, but there is no evidence for systematic changes in period over the 1000-d span of our observations. Rotational modulation and stellar-wind variability appear to be unlikely candidates for the underlying mechanism; we suggest that the physical origin of the signal may be pulsation associated with low-l oscillatory convection modes.

Amplitude variability in satellite photometry of the non-radially pulsating O9.5 V star ζ Oph

We report a time-series analysis of satellite photometry of the non-radially pulsating Oe star zeta Oph, principally using data from SMEI obtained 2003--2008, but augmented with MOST and WIRE results. Amplitudes of the strongest photometric signals, at 5.18, 2.96, and 2.67/d, each vary independently over the 6-year monitoring period (from ca. 30 to <2 mmag at 5.18/d), on timescales of hundreds of days. Signals at 7.19/d and 5.18/d have persisted (or recurred) for around two decades. Supplementary spectroscopic observations show an H-alpha emission episode in 2006; this coincided with small increases in amplitudes of the three strongest photometric signals.

Dynamics of H ii regions around exiled O stars

At least 25 per cent of massive stars are ejected from their parent cluster, becoming runaways or exiles, travelling with often-supersonic space velocities through the interstellar medium (ISM). Their overpressurised H II regions impart kinetic energy and momentum to the ISM, compress and/or evaporate dense clouds, and can constrain properties of both the star and the ISM. Here we present one-, two-, and (the first) three-dimensional simulations of the H II region around a massive star moving supersonically through a uniform, magnetised ISM, with properties appropriate for the nearby O star Zeta Oph. The H II region leaves an expanding overdense shell behind the star and, inside this, an underdense wake that should be filled with hot gas from the shocked stellar wind. The gas column density in the shell is strongly influenced by the ISM magnetic field strength and orientation. H-alpha emission maps show the H II region remains roughly circular, although the star is displaced somewhat from the centre of emission. For our model parameters, the kinetic energy feedback from the H II region is comparable to the mechanical luminosity of the stellar wind, and the momentum feedback rate is >100X larger than that from the wind and about 10X larger than the total momentum input rate available from radiation pressure. Compared to the star's eventual supernova explosion, the kinetic energy feedback from the H II region over the star's main sequence lifetime is >100X less, but the momentum feedback is up to 4X larger. H II region dynamics are found to have only a small effect on the ISM conditions that a bow shock close to the star would encounter.

GAS DISTRIBUTION, KINEMATICS, AND EXCITATION STRUCTURE IN THE DISKS AROUND THE CLASSICAL Be STARS β CANIS MINORIS AND ζ TAURI

Using CHARA and VLTI near-infrared spectro-interferometry with hectometric baseline lengths (up to 330m) and with high spectral resolution (up to 12000), we studied the gas distribution and kinematics around two classical Be stars. The combination of high spatial and spectral resolution achieved allows us to constrain the gas velocity field on scales of a few stellar radii and to obtain, for the first time in optical interferometry, a dynamical mass estimate using the position-velocity analysis technique known from radio astronomy. For our first target star, Beta Canis Minoris, we model the H+K-band continuum and Br-gamma line geometry with a near-critical rotating stellar photosphere and a geometrically thin equatorial disk. Testing different disk rotation laws, we find that the disk is in Keplerian rotation (v(r)=r^(-0.5+/-0.1)) and derive the disk position angle (140+/-1.7 deg), inclination (38.5+/-1 deg), and the mass of the central star (3.5+/-0.2 M_sun). As a second target star, we observed the prototypical Be star Zeta Tauri and spatially resolved the Br{\gamma} emission as well as nine transitions from the hydrogen Pfund series (Pf14-Pf22). Comparing the spatial origin of the different line transitions, we find that the Brackett (Br-gamma), Pfund (Pf14-17), and Balmer (H-alpha) lines originate from different stellocentric radii (R_cont < R_Pf < R_Br-gamma = R_H-alpha), which we can reproduce with an LTE line radiative transfer computation. Discussing different disk-formation scenarios, we conclude that our constraints are inconsistent with wind compression models predicting a strong outflowing velocity component, but support viscous decretion disk models, where the Keplerian-rotating disk is replenished with material from the near-critical rotating star.

First detection of a magnetic field in the fast rotating runaway Oe star ζ Ophiuchi

AbstractThe star ζ Ophiuchi is one of the brightest massive stars in the northern hemisphere and was intensively studied in various wavelength domains. The currently available observational material suggests that certain observed phenomena are related to the presence of a magnetic field. We acquired spectropolarimetric observations of ζ Oph with FORS 1 mounted on the 8‐m Kueyen telescope of the VLT to investigate if a magnetic field is indeed present in this star. Using all available absorption lines, we detect a mean longitudinal magnetic field 〈Bz〉all = 141 ± 45 G, confirming the magnetic nature of this star. We review the X‐ray properties of ζ Oph with the aim to understand whether the X‐ray emission of ζ Oph is dominated by magnetic or by wind instability processes (© 2011 WILEY‐VCH Verlag GmbH &amp; Co. KGaA, Weinheim)

MULTI-EPOCH NEAR-INFRARED INTERFEROMETRY OF THE SPATIALLY RESOLVED DISK AROUND THE BE STAR ζ TAU

We present interferometric observations of the Be star Zeta Tau obtained using the MIRC beam combiner at the CHARA Array. We resolved the disk during four epochs in 2007-2009. We fit the data with a geometric model to characterize the circumstellar disk as a skewed elliptical Gaussian and the central Be star as a uniform disk. The visibilities reveal a nearly edge-on disk with a FWHM major axis of ~ 1.8 mas in the H-band. The non-zero closure phases indicate an asymmetry within the disk. Interestingly, when combining our results with previously published interferometric observations of Zeta Tau, we find a correlation between the position angle of the disk and the spectroscopic V/R ratio, suggesting that the tilt of the disk is precessing. This work is part of a multi-year monitoring campaign to investigate the development and outward motion of asymmetric structures in the disks of Be stars.

DISCOVERY AND CHARACTERIZATION OF A FAINT STELLAR COMPANION TO THE A3V STAR ζ VIRGINIS

Through the combination of high-order Adaptive Optics and coronagraphy, we report the discovery of a faint stellar companion to the A3V star zeta Virginis. This companion is ~7 magnitudes fainter than its host star in the H-band, and infrared imaging spanning 4.75 years over five epochs indicates this companion has common proper motion with its host star. Using evolutionary models, we estimate its mass to be 0.168+/-.016 solar masses, giving a mass ratio for this system q = 0.082. Assuming the two objects are coeval, this mass suggests a M4V-M7V spectral type for the companion, which is confirmed through integral field spectroscopic measurements. We see clear evidence for orbital motion from this companion and are able to constrain the semi-major axis to be greater than 24.9 AU, the period > 124$ yrs, and eccentricity > 0.16. Multiplicity studies of higher mass stars are relatively rare, and binary companions such as this one at the extreme low end of the mass ratio distribution are useful additions to surveys incomplete at such a low mass ratio. Moreover, the frequency of binary companions can help to discriminate between binary formation scenarios that predict an abundance of low-mass companions forming from the early fragmentation of a massive circumstellar disk. A system such as this may provide insight into the anomalous X-ray emission from A stars, hypothesized to be from unseen late-type stellar companions. Indeed, we calculate that the presence of this M-dwarf companion easily accounts for the X-ray emission from this star detected by ROSAT.

X-ray temperatures and their radial distributions from winds of O-type supergiants: the effects of clumps

Context. Based on the multiphase hydrodynamical flows and the clump model, a two-fluid model of the O stars zeta Ori A and zeta Pup is investigated, in which the flow is described by a set of two components. Aims. I investigate the influence of clumps on the production of X-ray and derive wind parameters from a comparison with the X-ray temperature distributions in the wind of O-type supergiants. Methods. The velocity structure of the clump is assumed to be beta law, thus the transonic velocity structure of the ambient wind can be attained via the shooting method. Results. The mass, filling factor, and velocity structure of the clump are very sensitive to the post-shock temperatures. I find that the filling factor must be about 0.01 for fitting the observed X-ray temperatures. The mass of the clump cannot be smaller than 10(18) g. A clump carries about 10(-4) (M) over dot and 10(-5) (M) over dot for zeta Ori A and zeta Pup, respectively. The model with the exponent of clump velocity beta = 0.5-0.6 for O-type supergiants fits the observations better than one with the exponent beta = 0.8. Conclusions. These results indicate that the forward shock scenario can explain the X-ray emission, except for the hottest X-ray lines for O-type supergiants. Theoretical and observable evidence shows that a low velocity law should be applied to O-type supergiants.

Cyclic variability of the circumstellar disk of the Be starζ Tauri

Emission lines formed in decretion disks of Be stars often undergo long-term cyclic variations, especially in the violet-to-red (V/R) ratio of their primary components. From observations of the bright Be-shell star zeta Tau, the possibly broadest and longest data set illustrating the prototype of this behaviour was compiled from our own and archival observations. It comprises optical and infrared spectra, broad-band polarimetry, and interferometric observations. From 3 V/R cycles between 1997 and 2008, a mean cycle length in H alpha of 1400-1430 days was derived. After each minimum in V/R, the shell absorption weakens and splits into two components, leading to 3 emission peaks. This phase makes the strongest contribution to the variability in cycle length. V/R curves of different lines are shifted in phase. Lines formed on average closer to the central star are ahead of the others. The shell absorption lines fall into 2 categories differing in line width, ionization/excitation potential, and variability of the equivalent width. The interferometry has resolved the continuum and the line emission in Br gamma and He I 2.06. The phasing of the Br gamma emission shows that the photocenter of the line-emitting region lies within the plane of the disk but is offset from the continuum source. The plane of the disk is constant throughout the observed V/R cycles. The self-consistent, one-armed, disk-oscillation model is developed in Paper II.

CO emission and variable CH and CH+absorption towards HD 34078: evidence for a nascent bow shock?

The runaway star HD34078, initially selected to investigate small scale structure in a foreground diffuse cloud has been shown to be surrounded by highly excited H2. We first search for an association between the foreground cloud and HD34078. Second, we extend previous investigations of temporal absorption line variations (CH, CH+, H2) in order to better characterize them. We have mapped the CO(2-1) emission at 12 arcsec resolution around HD34078's position, using the 30 m IRAM antenna. The follow-up of CH and CH+ absorption lines has been extended over 5 more years. In parallel, CH absorption towards the reddened star Zeta Per have been monitored to check the homogeneity of our measurements. Three more FUSE spectra have been obtained to search for N(H2) variations. CO observations show a pronounced maximum near HD34078's position, clearly indicating that the star and diffuse cloud are associated. The optical spectra confirm the reality of strong, rapid and correlated CH and CH+ fluctuations. On the other hand, N(H2, J=0) has varied by less than 5 % over 4 years. We also discard N(CH) variations towards Zeta Per at scales less than 20 AU. Observational constraints from this work and from 24 micron dust emission appear to be consistent with H2 excitation but inconsistent with steady-state bow shock models and rather suggest that the shell of compressed gas surrounding HD34078, is seen at an early stage of the interaction. The CH and CH+ time variations as well as their large abundances are likely due to chemical structure in the shocked gas layer located at the stellar wind/ambient cloud interface. Finally, the lack of variations for both N(H2, J=0) towards HD34078 and N(CH) towards Zeta Per suggests that quiescent molecular gas is not subject to pronounced small-scale structure.

DIRECT SPECTROSCOPIC OBSERVATIONS OF CLUMPING IN O-STAR WINDS

We report the detection and monitoring of transient substructures in the radiation-driven winds of five massive, hot stars in different evolutionary stages. Clumping in the winds of these stars shows up as variable, narrow subpeaks superposed on their wide, wind-broadened (optical) emission lines. Similar patterns of emission-line profile variations are detected in the Of stars zeta Puppis and HD93129A, in the more evolved hydrogen-rich, luminous, Of-like WN stars HD93131 and HD93162, and in the more mass-depleted WC star in gamma2 Velorum. These observations strongly suggest that stochastic wind clumping is a universal phenomenon in the radiation-driven, hot winds from all massive stars, with similar clumping factors in all stages of mass depletion.

Estimating the Porosity of the Interstellar Medium from Three‐dimensional Photoionization Modeling of HiiRegions

We apply our three dimensional photoionization code to model Wisconsin H-alpha Mapper observations of the H II region surrounding the O9.5V star Zeta Oph. Our models investigate the porosity of the interstellar medium around zeta Oph and the effects of 3D densities on the H-alpha surface brightness and variation in the [N II]/H-alpha line ratio. The Zeta Oph H II region has a well characterized ionizing source, so it is an excellent starting point for 3D models of diffuse ionized gas. We investigate various hierarchically clumped density structures, varying the overall smoothness within the clumping algorithm. By simulating the observations, we can estimate the porosity of the medium in the vicinity of Zeta Oph and find that within the context of our hierarchically clumped models, around 50% to 80% of the volume is occupied by clumps surrounded by a low density smooth medium. We also conclude that in order for O stars to ionize the diffuse Warm Ionized Medium, the O star environment must be more porous than that surrounding Zeta Oph, with clumps occupying less than one half of the interstellar volume. Our clumpy models have irregular boundaries, similar to observed H II regions. However, in observed H II regions it is difficult to identify the precise location of the boundary because of the foreground and/or background emission from the widespread Warm Ionized Medium. This complicates the interpretation of the predicted rapid rise of some emission line ratios near the edge of uniform density H II regions and combined with the three dimensional clumpy nature of the interstellar medium may explain the apparent lack of distinctive emission line ratios near H I -- H II interfaces.

Radio and submillimetre observations of wind structure in ζ Puppis

We present radio and submillimetre observations of the O4I(n)f star zeta Pup, and discuss structure in the outer region of its wind (similar to10-100 R-*). The properties of bremsstrahlung, the dominant emission process at these wavelengths, make it sensitive to structure and allow us to study how the amount of structure changes in the wind by comparing the fluxes at different wavelengths. Possible forms of structure at these distances include Corotating Interaction Regions (CIRs), stochastic clumping, a disk or a polar enhancement. As the CIRs are azimuthally asymmetric, they should result in variability at submillimetre or radio wavelengths. To look for this variability, we acquired 3.6 and 6 cm observations with the Australia Telescope Compact Array (ATCA), covering about two rotational periods of the star. We supplemented these with archive observations from the NRAO Very Large Array (VLA), which cover a much longer time scale. We did not find variability at more than the +/-20% level. The long integration time does allow an accurate determination of the fluxes at 3.6 and 6 cm. Converting these fluxes into a mass loss rate, we find. (M)over dot = 3.5 x 10(-6) M-circle dot/yr. This value confirms the significant discrepancy with the mass loss rate derived from the Halpha profile, making zeta Pup an exception to the usually good agreement between the Halpha and radio mass loss rates. To study the run of structure as a function of distance, we supplemented the ATCA data by observing zeta Pup at 850 mum with the James Clerk Maxwell Telescope (JCMT) and at 20 cm with the VLA. A smooth wind model shows that the millimetre fluxes are too high compared to the radio fluxes. While recombination of helium in the outer wind cannot be discounted as an explanation, the wealth of evidence for structure strongly suggests this as the explanation for the discrepancy. Model calculations show that the structure needs to be present in the inner similar to70 R-* of the wind, but that it decays significantly, or maybe even disappears, beyond that radius.

X-ray line emission from a fragmented stellar wind

We discuss X-ray line formation in dense O star winds. A random distribution of wind shocks is assumed to emit X-rays that are partially absorbed by cooler wind gas. The cool gas resides in highly compressed fragments oriented perpendicular to the radial flow direction. For fully opaque fragments, we find that the blueshifted part of X-ray line profiles remains flat-topped even after severe wind attenuation, whereas the red part shows a steep decline. These box-type, blueshifted profiles resemble recent Chandra observations of the O3 star zeta Pup. For partially transparent fragments, the emission lines become similar to those from a homogeneous wind.

The Proper Motion, Parallax, and Origin of the Isolated Neutron Star RX J185635−3754

The isolated neutron star RX J185635-3754 is the closest known neutron star to the Sun. Based on HST WFPC2 obervations over a 3 year baseline, I report its proper motion (332 +/- 1 mas/yr at a position angle of 100.3 +/- 0.1 degrees) and parallax (16.5 +/- 2.3 mas; 61 pc). This proper motion brings the neutron star from the general vicinity of the Sco-Cen OB association. For an assumed neutron star radial velocity between -55 and -60 km/s, the runaway O star zeta Oph, the Upper Sco OB association, and the neutron star come into spatial coincidence between 0.9 and 1.0 million years ago. RX J185635-3754 may be the remnant of the original primary of the zeta Oph system. If so, the space velocity suggests that the neutron star received a kick of about 200 km/s at birth.

Direct detection of pulsations of the Cepheid star zeta Gem and an independent calibration of the period-luminosity relation.

Cepheids are a class of variable (pulsating) stars whose absolute luminosities are related in a simple manner to their pulsational periods. By measuring the period and using the 'period-luminosity' relationship, astronomers can use the observed visual brightness to determine the distance to the star. Because these stars are very luminous, they can be observed in other galaxies, and therefore can be used to help determine the expansion rate of the Universe (the Hubble constant). Calibration of the period-luminosity relation is a necessary first step, but the small number of sufficiently nearby Cepheids has forced the use of a number of indirect means, with associated systematic uncertainties. Here we present a distance to the Cepheid zeta Geminorum, determined using a direct measurement (by an optical interferometer) of its changes in diameter as it pulsates. Within our uncertainty of 15 per cent, our distance is in agreement with previous indirect determinations. Planned improvements to the instrument will allow us to calibrate directly the period-luminosity relation to better than a few per cent.

Classification of O Stars in the Yellow-Green: The Exciting Star VES 735

Acquiring data for spectral classification of heavily reddened stars using traditional criteria in the blue-violet region of the spectrum can be prohibitively time consuming using small to medium sized telescopes. One such star is the Vatican Observatory emission-line star VES 735, which we have found excites the HII region KR140. In order to classify VES 735, we have constructed an atlas of stellar spectra of O stars in the yellow-green (4800 - 5420 A). We calibrate spectral type versus the line ratio HeI 4922:HeII 5411, showing that this ratio should be useful for the classification of heavily reddened O stars associated with HII regions. Application to VES 735 shows that the spectral type is O8.5. The absolute magnitude suggests luminosity class V. Comparison of the rateof emission of ionizing photons and the bolometric luminosity of VES 735, inferred from radio and infrared measurements of the KR 140 region, to recent stellar models gives consistent evidence for a main sequence star of mass 25 M_sun and age less than a few million years with a covering factor 0.4 - 0.5 by the nebular material. Spectra taken in the red (6500 - 6700 A) show that the stellar H-alpha emission is double-peaked about the systemic velocity and slightly variable. H-beta is in absorption, so that the emission-line classification is ``(e)''. However, unlike the case of the more well known O(e) star Zeta Oph, the emission from VES 735 appears to be long-lived rather than episodic.

Imaging binary stars by the cross-correlation technique

We present in this paper a technique for imaging binary stars from speckle data. This technique is based upon the computation of the cross-correlation between the speckle frames and their square. This may be considered as a simple, easy to implement, complementary computation to the autocorrelation function of Labeyrie's technique for a rapid determination of the position angle of binary systems. Angular separation, absolute position angle and relative photometry of binary stars can be derived from this technique. We show an application to the bright double star zeta Sge observed at the 2m Telescope Bernard Lyot.

The [ITAL]IUE[/ITAL] MEGA Campaign: The Rotationally Modulated Wind of ζ Puppis

We discuss 16 days of intensive International Ultraviolet Explorer (IUE) observations of the Si iv doublet (lambda1400) in the spectrum of the O4 I(n) f star zeta Pup. The data show continuous variability throughout the greater part of the blueshifted absorption. Time series analysis of these data reveals significant power at periods of 19.2 hr and 5.2 days, which we identify with the mean recurrence time of Discrete Absorption Components (DACs) and the photospheric rotation period, respectively. These results indicate that the wind has a global longitudinal asymmetry (approaching a factor 2 in optical depth), possibly associated with large-scale magnetic structures, but suggest that the DACs are not directly associated with specific stellar longitudes in this star. There is no significant power in the lines at the 8.5 hr period identified in photospheric absorption-line variability, nor at the 16.7 hr period reported in X-ray observations.

A STUDY OF THE COMETARY GLOBULES IN THE GUM NEBULA

ABSTRACT The first part of the thesis deals with the development of a wide-band mechanically tuned local oscillator using the Gunn diode for use with the 10.4m millimeter-wave radio telescope at the Raman Research Institute. This provides sufficient power to efficiently operate two cryogenic Schottky mixers (dual polarisation) and tunes over the frequency range 75-115 GHz covering most of the 3-mm atmospheric transmission window (W-band). Rotational transitions of many astrophysically important molecules including CO fall in this range. A study of the cometary globulesCGs in the Gum Nubula forms the second part. The CGs are characterised by compact, dusty heads with long faintly luminous tails extending on one side and narrow bright rims on the other side. There exists a significant population of such CGs in the Gum Nebula, distributed over a region ~80 parsec in radius with their tails pointing away from an apparent common center. Some of the heads have embedded young stars. In the region bounded by the CGs there are a few massive hot stars including zeta Puppis believed to be the most luminous star in the southern sky. It has been suspected that the morphological appearance of the CGs may be due to the influence of these stars. In order to understand the kinematics and the origin of the system, a study was undertaken using the first rotational transition of the carbon monoxide molecule. The study consisted of 12CO observations of the heads and the tails of the CGs. In addition, the Globule No.22 was mapped in both 12CO and 13CO. An analysis of this data has led to the following findings: 1. The system of CGs is expanding with respect to a common morphological center at ~12 kms-1. The expansion age is ~6 Myr. 2. Some of the tails observed show systematic velocity gradients. If the tails were formed due to the elongation resulting from these velocity gradients then the estimated stretching age is ~3 Myr. In order to clarify if externally triggered star formation is going on in this region, an analysis of the locations of the embedded young stellar objects (YSOs, identified from the Infra-Red Astronomy Satellite data) in the dark clouds in the Gum-Vela region was undertaken. This study has shown that the YSOs have a statistically significant tendency to fall on the sides of dark clouds facing the morphological center rather than the far sides, supporting external triggering. From the above analysis we come to the following conclusions: 1. The rough agreement between the expansion age and the tail-stretching age suggests a common origin for the expansion and the formation of the tails. The presence of young stars of comparable ages in the heads of some of the globules suggests that the processes responsible for the expansion may have also triggered star formation in them. 2. The radiation pressure from the hot stars in the central region or the stellar winds from them cannot account for the momentum of the expanding globules. It is more likely that the rocket effect arising out of the heating and the consequent anisotropic ablation of the globules supplied the necessary momentum. From the space motion of the star zeta Puppis we suggest that it had a massive companion which exploded as a supernova half a million years ago. The combined effect of the ultra-violet radiation and the stellar wind from this binary as well as from other stars in the neighbourhood significantly affected the parent molecular cloud in which they formed. This resulted in much of the molecular material in the vicinity being blown away except the numerous regions of enhanced density (condensations) in the original molecular cloud. Continued effect of the radiation and stellar winds resulted in these condensations being set in motion, as well as developing cometary tails.