Using MOST to reveal the secrets of the mischievous Wolf-Rayet binary CV Ser

The WR binary CV Serpentis (= WR113, WC8d + O8-9IV) has been a source of mystery since it was shown that its atmospheric eclipses change with time over decades, in addition to its sporadic dust production. The first high-precision time-dependent photometric observations obtained with the MOST space telescope in 2009 show two consecutive eclipses over the 29d orbit, with varying depths. A subsequent MOST run in 2010 showed a seemingly asymmetric eclipse profile. In order to help make sense of these observations, parallel optical spectroscopy was obtained from the Mont Megantic Observatory (2009, 2010) and from the Dominion Astrophysical Observatory (2009). Assuming these depth variations are entirely due to electron scattering in a beta-law wind, an unprecedented 62% increase in mass-loss rate is observed over one orbital period. Alternatively, no change in mass-loss rate would be required if a relatively small fraction of the carbon ions in the wind globally recombined and coaggulated to form carbon dust grains. However, it remains a mystery as to how this could occur. There also seems to be evidence for the presence of corotating interaction regions (CIR) in the WR wind: a CIR-like signature is found in the light curves, implying a potential rotation period for the WR star of 1.6 d. Finally, a new circular orbit is derived, along with constraints for the wind collision.Comment: 11 pages, 11 figures, 5 table

Pulsation models for the roAp star HD 134214

Precise time-series photometry with the MOST satellite has led to identification of 10 pulsation frequencies in the rapidly oscillating Ap (roAp) star HD 134214. We have fitted the observed frequencies with theoretical frequencies of axisymmetric modes in a grid of stellar models with dipole magnetic fields. We find that, among models with a standard composition of $(X,Z) = (0.70,0.02)$ and with suppressed convection, eigenfrequencies of a $1.65\,{\rm M}_\odot$ model with $\log T_{\rm eff} = 3.858$ and a polar magnetic field strength of 4.1kG agree best with the observed frequencies. We identify the observed pulsation frequency with the largest amplitude as a deformed dipole ($\ell = 1$) mode, and the four next-largest-amplitude frequencies as deformed $\ell = 2$ modes. These modes have a radial quasi-node in the outermost atmospheric layers ($\tau \sim 10^{-3}$). Although the model frequencies agree roughly with observed ones, they are all above the acoustic cut-off frequency for the model atmosphere and hence are predicted to be damped. The excitation mechanism for the pulsations of HD 134214 is not clear, but further investigation of these modes may be a probe of the atmospheric structure in this magnetic chemically peculiar star.Comment: 9 pages, 6 figures; accepted for publication in MNRA

Dynamics and Structure of Three-Dimensional Poloidally Magnetized Supermagnetosonic Jets

A set of 3D MHD simulations of magnetized jets has been performed. The jets contain an equipartition primarily poloidal magnetic field and the effect of jet density on jet dynamics and structure is evaluated. The jet is precessed at the origin to excite Kelvin-Helmholtz unstable helical modes. We extensively compare the structure in these simulations with linear stability theory. The jet that is dense with respect to the external medium develops a high speed core surrounded by a less dense sheath consisting of slower moving jet fluid. These simulations suggest that extended extragalactic jets propagate to such large distances because they are surrounded by a lobe or cocoon whose density is less than the jet density. (Abridged abstract.)Comment: 30 pages, AASTeX, to appear in ApJ, much better versions of Figures 2-5 are available at http://crux.astr.ua.edu/~rosen/hcr/hcr.htm

Massive open star clusters using the VVV survey, I: presentation of the data and description of the approach

Context. The ESO Public Survey “VISTA Variables in the Vía Láctea” (VVV) provides deep multi-epoch infrared observations for unprecedented 562 sq. degrees of the Galactic bulge, and adjacent regions of the disk. Aims. The VVV observations will foster the construction of a sample of Galactic star clusters with reliable and homogeneously derived physical parameters (e.g., age, distance, and mass, etc.). In this first paper in a series, the methodology employed to establish cluster parameters for the envisioned database are elaborated upon by analysing four known young open clusters: Danks 1, Danks 2, RCW 79, and DBS 132. The analysis offers a first glimpse of the information that can be gleaned from the VVV observations for clusters in the final database. Methods. Wide-field, deep JHKs VVV observations, combined with new infrared spectroscopy, are employed to constrain fundamental parameters for a subset of clusters. Results. Results are inferred from VVV near-infrared photometry and numerous low resolution spectra (typically more than 10 per cluster). The high quality of the spectra and the deep wide–field VVV photometry enables us to precisely and independently determine the characteristics of the clusters studied, which we compare to previous determinations. An anomalous reddening law in the direction of the Danks clusters is found, specifically E(J − H)/E(H − Ks) = 2.20 ± 0.06, which exceeds published values for the inner Galaxy. The G305 star forming complex, which includes the Danks clusters, lies beyond the Sagittarius-Carina spiral arm and occupies the Centaurus arm. Finally, the first deep infrared colour-magnitude diagram of RCW 79 is presented, which reveals a sizeable pre-main sequence population. A list of candidate variable stars in G305 region is reported. Conclusions. This study demonstrates the strength of the dataset and methodology employed, and constitutes the first step of a broader study which shall include reliable parameters for a sizeable number of poorly characterised and/or newly discovered clusters.Fil: Chené, A.-N.. Universidad de Valparaíso; ChileFil: Borissova, J.. Universidad de Valparaíso; ChileFil: Clarke, J. R. A.. Universidad de Valparaíso; ChileFil: Bonatto, C.. Universidade Federal do Rio Grande do Sul; BrasilFil: Majaess, D. J.. Saint Marys University; CanadáFil: Moni Bidin, C.. Universidad de Concepción; ChileFil: Sale, S. E.. Universidad de Valparaíso; ChileFil: Mauro, F.. Universidad de Concepción; ChileFil: Kurtev, R.. Universidad de Valparaíso; ChileFil: Baume, Gustavo Luis. Universidad Nacional de La Plata. Facultad de Ciencias Astronómicas y Geofísicas; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto de Astrofísica La Plata. Universidad Nacional de La Plata. Facultad de Ciencias Astronómicas y Geofísicas. Instituto de Astrofísica La Plata; ArgentinaFil: Feinstein Baigorri, Carlos. Universidad Nacional de La Plata. Facultad de Ciencias Astronómicas y Geofísicas; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto de Astrofísica La Plata. Universidad Nacional de La Plata. Facultad de Ciencias Astronómicas y Geofísicas. Instituto de Astrofísica La Plata; ArgentinaFil: Ivanov, V. D.. European Southern Observatory; ChileFil: Geisler, Doug. Universidad de Concepción; ChileFil: Catelan, M.. The Milky Way Millennium Nucleus; ChileFil: Minniti, Dante. The Milky Way Millennium Nucleus; ChileFil: Lucas, P.. University Of Hertfordshire; Reino UnidoFil: de Grijs, R.. Peking University; ChinaFil: Kumar, M. S. N.. Centro de Astrofisica da Universidade do Porto; Portuga

A search for p-modes and other variability in the binary system 85 Pegasi using MOST photometry

Context: Asteroseismology has great potential for the study of metal-poor stars due to its sensitivity to determine stellar ages. Aims: Our goal was to detect p-mode oscillations in the metal-poor sub-dwarf 85 Peg A and to search for other variability on longer timescales. Methods: We have obtained continuous high-precision photometry of the binary system 85 Pegasi with the MOST space telescope in two seasons (2005 & 2007). Furthermore, we redetermined vsini for 85 Peg A using high resolution spectra obtained through the ESO archive, and used photometric spot modeling to interpret long periodic variations. Results: Our frequency analysis yields no convincing evidence for p-modes significantly above a noise level of 4 ppm. Using simulated p-mode patterns we provide upper RMS amplitude limits for 85 Peg A. The light curve shows evidence for variability with a period of about 11 d and this periodicity is also seen in the follow up run in 2007; however, as different methods to remove instrumental trends in the 2005 run yield vastly different results, the exact shape and periodicity of the 2005 variability remain uncertain. Our re-determined vsini value for 85 Peg A is comparable to previous studies and we provide realistic uncertainties for this parameter. Using these values in combination with simple photometric spot models we are able to reconstruct the observed variations. Conclusions: The null-detection of p-modes in 85 Peg A is consistent with theoretical values for pulsation amplitudes in this star. The detected long-periodic variation must await confirmation by further observations with similar or better precision and long-term stability. If the 11 d periodicity is real, rotational modulation of surface features on one of the components is the most likely explanation.Comment: 11 pages, 9 figures, accepted for publication in A&

The Formation of Giant Elliptical Galaxies and Their Globular Cluster Systems

The bimodal globular cluster (GC) metallicity distributions of many giant elliptical galaxies are often cited as evidence for the formation of such galaxies through mergers involving gas-rich spirals. In such models, the metal- rich GCs are assumed to have formed during the merger process. We explore an alternative possibility: that these metal-rich clusters represent the galaxy's intrinsic GC population and that the metal-poor component of the observed GC metallicity distribution arises from the capture of GCs from other galaxies, either through mergers or through tidal stripping. Starting with plausible assumptions for the initial galaxy luminosity function and for the dependence of GC metallicity on parent galaxy luminosity, we show that the growth of a pre-existing seed galaxy through mergers and tidal stripping is accompanied by the capture of metal-poor GCs whose properties are similar to those which are observed to surround giant ellipticals. We describe a method of using the observed number of metal-poor and metal-rich GCs to infer the merger histories of individual elliptical galaxies, and use this technique to derive limits on the number of galaxies and total luminosity accreted to date by M49. We argue that although GC specific frequency is conserved in galaxy mergers, the same may not be true of tidal stripping by the mean field of the host galaxy cluster. Comparisons of model GC metallicity distributions and specific frequencies to those observed for the well-studied galaxies M49 and M87 show that it is possible to explain their bimodal GC metallicity distributions and discordant specific frequencies without resorting to the formation of new GCs in mergers or by invoking multiple bursts of GC formation.Comment: 39 pages AAS Latex and 10 postscript figures. Also available at http://astro.caltech.edu/~pc. Accepted for publication in the Astrophysical Journa

Modeling Convective Core Overshoot and Diffusion in Procyon Constrained by Asteroseismic Data

We compare evolved stellar models, which match Procyons mass and position in the HR diagram, to current ground-based asteroseismic observations. Diffusion of helium and metals along with two conventional core overshoot descriptions and the Kuhfuss nonlocal theory of convection are considered. We establish that one of the two published asteroseismic data reductions for Procyon, which mainly differ in their identification of even versus odd l-values, is a significantly more probable and self-consistent match to our models than the other. The most probable models according to our Bayesian analysis have evolved to just short of turnoff, still retaining a hydrogen convective core. Our most probable models include Y and Z diffusion and have conventional core overshoot between 0.9 and 1.5 pressure scale heights, which increases the outer radius of the convective core by between 22% to 28%, respectively. We discuss the significance of this comparatively higher than expected core overshoot amount in terms of internal mixing during evolution. The parameters of our most probable models are similar regardless of whether adiabatic or nonadiabatic model p-mode frequencies are compared to the observations, although, the Bayesian probabilities are greater when the nonadiabatic model frequencies are used. All the most probable models (with or without core overshoot, adiabatic or nonadiabatic model frequencies, diffusion or no diffusion, including priors for the observed HRD location and mass or not) have masses that are within one sigma of the observed mass 1.497+/-0.037 Msun

The rise of an ionized wind in the Narrow Line Seyfert 1 Galaxy Mrk 335 observed by XMM-Newton and HST

We present the discovery of an outflowing ionized wind in the Seyfert 1 Galaxy Mrk 335. Despite having been extensively observed by most of the largest X-ray observatories in the last decade, this bright source was not known to host warm absorber gas until recent XMM-Newton observations in combination with a long-term Swift monitoring program have shown extreme flux and spectral variability. High resolution spectra obtained by the XMM-Newton RGS detector reveal that the wind consists of three distinct ionization components, all outflowing at a velocity of 5000 km/s. This wind is clearly revealed when the source is observed at an intermediate flux state (2-5e-12 ergs cm^-2 s^-1). The analysis of multi-epoch RGS spectra allowed us to compare the absorber properties at three very different flux states of the source. No correlation between the warm absorber variability and the X-ray flux has been determined. The two higher ionization components of the gas may be consistent with photoionization equilibrium, but we can exclude this for the only ionization component that is consistently present in all flux states (log(xi)~1.8). We have included archival, non-simultaneous UV data from HST (FOS, STIS, COS) with the aim of searching for any signature of absorption in this source that so far was known for being absorption-free in the UV band. In the COS spectra obtained a few months after the X-ray observations we found broad absorption in CIV lines intrinsic to the AGN and blueshifted by a velocity roughly comparable to the X-ray outflow. The global behavior of the gas in both bands can be explained by variation of the covering factor and/or column density, possibly due to transverse motion of absorbing clouds moving out of the line of sight at Broad Line Region scale.Comment: 24 pages, 12 figures, ApJ accepte

The Pulsation Properties of Procyon A

A grid of stellar evolution models for Procyon A has been calculated. These models include the best physics available to us (including the latest opacities and equation of state) and are based on the revised astrometric mass of Girard et al (1996). Models were calculated with helium diffusion and with the combined effects of helium and heavy element diffusion. Oscillation frequencies for l=0,1,2 and 3 p-modes and the characteristic period spacing for the g-modes were calculated for these models. We find that g-modes are sensitive to model parameters which effect the structure of the core, such as convective core overshoot, the heavy element abundance and the evolutionary state (main sequence or shell hydrogen burning) of Procyon A. The p-modes are relatively insensitive to the details of the physics used to model Procyon A, and only depend on the evolutionary state of Procyon A. Hence, observations of p-mode frequencies on Procyon A will serve as a robust test of stellar evolution models.Comment: 4 pages, to appear in ApJ

MaGICC baryon cycle: the enrichment history of simulated disc galaxies

Using cosmological galaxy formation simulations from the MaGICC (Making Galaxies in a Cosmological Context) project, spanning stellar mass from ∼107 to 3 × 1010 M⊙, we trace the baryonic cycle of infalling gas from the virial radius through to its eventual participation in the star formation process. An emphasis is placed upon the temporal history of chemical enrichment during its passage through the corona and circumgalactic medium. We derive the distributions of time between gas crossing the virial radius and being accreted to the star-forming region (which allows for mixing within the corona), as well as the time between gas being accreted to the star-forming region and then ultimately forming stars (which allows for mixing within the disc). Significant numbers of stars are formed from gas that cycles back through the hot halo after first accreting to the star-forming region. Gas entering high-mass galaxies is pre-enriched in low-mass proto-galaxies prior to entering the virial radius of the central progenitor, with only small amounts of primordial gas accreted, even at high redshift (z ∼ 5). After entering the virial radius, significant further enrichment occurs prior to the accretion of the gas to the star-forming region, with gas that is feeding the star-forming region surpassing 0.1 Z⊙ by z = 0. Mixing with halo gas, itself enriched via galactic fountains, is thus crucial in determining the metallicity at which gas is accreted to the disc. The lowest mass simulated galaxy (Mvir ∼ 2 × 1010 M⊙, with M⋆ ∼ 107 M⊙), by contrast, accretes primordial gas through the virial radius and on to the disc, throughout its history. Much like the case for classical analytical solutions to the so-called ‘G-dwarf problem’, overproduction of low-metallicity stars is ameliorated by the interplay between the time of accretion on to the disc and the subsequent involvement in star formation – i.e. due to the inefficiency of star formation. Finally, gas outflow/metal removal rates from star-forming regions as a function of galactic mass are presented