Observational Tests and Predictive Stellar Evolution II: Non-standard Models

We examine contributions of second order physical processes to results of stellar evolution calculations amenable to direct observational testing. In the first paper in the series (Young et al. 2001) we established baseline results using only physics which are common to modern stellar evolution codes. In the current paper we establish how much of the discrepancy between observations and baseline models is due to particular elements of new physics. We then consider the impact of the observational uncertainties on the maximum predictive accuracy achievable by a stellar evolution code. The sun is an optimal case because of the precise and abundant observations and the relative simplicity of the underlying stellar physics. The Standard Model is capable of matching the structure of the sun as determined by helioseismology and gross surface observables to better than a percent. Given an initial mass and surface composition within the observational errors, and no additional constraints for which the models can be optimized, it is not possible to predict the sun's current state to better than ~7%. Convectively induced mixing in radiative regions, seen in multidimensional hydrodynamic simulations, dramatically improves the predictions for radii, luminosity, and apsidal motions of eclipsing binaries while simultaneously maintaining consistency with observed light element depletion and turnoff ages in young clusters (Young et al. 2003). Systematic errors in core size for models of massive binaries disappear with more complete mixing physics, and acceptable fits are achieved for all of the binaries without calibration of free parameters. The lack of accurate abundance determinations for binaries is now the main obstacle to improving stellar models using this type of test.Comment: 33 pages, 8 figures, accepted for publication in the Astrophysical Journa

The Long Period, Massive Binaries HD 37366 and HD 54662: Potential Targets for Long Baseline Optical Interferometry

We present the results from an optical spectroscopic analysis of the massive stars HD 37366 and HD 54662. We find that HD 37366 is a double-lined spectroscopic binary with a period of 31.8187 +/- 0.0004 days, and HD 54662 is also a double lined binary with a much longer period of 557.8 +/- 0.3 days. The primary of HD 37366 is classified as O9.5 V, and it contributes approximately two-thirds of the optical flux. The less luminous secondary is a broad-lined, early B-type main-sequence star. Tomographic reconstruction of the individual spectra of HD 37366 reveals absorption lines present in each component, enabling us to constrain the nature of the secondary and physical characteristics of both stars. Tomographic reconstruction was not possible for HD 54662; however, we do present mean spectra from our observations that show that the secondary component is approximately half as bright as the primary. The observed spectral energy distributions (SEDs) were fit with model SEDs and galactic reddening curves to determine the angular sizes of the stars. By assuming radii appropriate for their classifications, we determine distance ranges of 1.4 - 1.9 and 1.2 - 1.5 kpc for HD 37366 and HD 54662, respectively.Comment: 27 pages, 8 figures, Accepted for publication in Ap

Radial Velocities of Six OB Stars

We present new results from a radial velocity study of six bright OB stars with little or no prior measurements. One of these, HD 45314, may be a long-period binary, but the velocity variations of this Be star may be related to changes in its circumstellar disk. Significant velocity variations were also found for HD 60848 (possibly related to nonradial pulsations) and HD 61827 (related to wind variations). The other three targets, HD 46150, HD 54879, and HD 206183, are constant velocity objects, but we note that HD 54879 has H$\alpha$ emission that may originate from a binary companion. We illustrate the average red spectrum of each target.Comment: Accepted for publication in PASP July 2007 issu

A Spectroscopic Study of Field and Runaway OB Stars

Identifying binaries among runaway O- and B-type stars offers valuable insight into the evolution of open clusters and close binary stars. Here we present a spectroscopic investigation of 12 known or suspected binaries among field and runaway OB stars. We find new orbital solutions for five single-lined spectroscopic binaries (HD 1976, HD 14633, HD 15137, HD 37737, and HD 52533), and we classify two stars thought to be binaries (HD 30614 and HD 188001) as single stars. In addition, we reinvestigate their runaway status using our new radial velocity data with the UCAC2 proper motion catalogs. Seven stars in our study appear to have been ejected from their birthplaces, and at least three of these runaways are spectroscopic binaries and are of great interest for future study.Comment: 21 pages, 1 figure, 7 tables; Accepted to Ap

Measurement of the CMS Magnetic Field

The measurement of the magnetic field in the tracking volume inside the superconducting coil of the Compact Muon Solenoid (CMS) detector under construction at CERN is done with a fieldmapper designed and produced at Fermilab. The fieldmapper uses 10 3-D B-sensors (Hall probes) developed at NIKHEF and calibrated at CERN to precision 0.05% for a nominal 4 T field. The precise fieldmapper measurements are done in 33840 points inside a cylinder of 1.724 m radius and 7 m long at central fields of 2, 3, 3.5, 3.8, and 4 T. Three components of the magnetic flux density at the CMS coil maximum excitation and the remanent fields on the steel-air interface after discharge of the coil are measured in check-points with 95 3-D B-sensors located near the magnetic flux return yoke elements. Voltages induced in 22 flux-loops made of 405-turn installed on selected segments of the yoke are sampled online during the entire fast discharge (190 s time-constant) of the CMS coil and integrated offline to provide a measurement of the initial magnetic flux density in steel at the maximum field to an accuracy of a few percent. The results of the measurements made at 4 T are reported and compared with a three-dimensional model of the CMS magnet system calculated with TOSCA.Comment: 4 pages, 5 figures, 15 reference

Physical Orbit for Lambda Virginis and a Test of Stellar Evolution Models

Lambda Virginis (LamVir) is a well-known double-lined spectroscopic Am binary with the interesting property that both stars are very similar in abundance but one is sharp-lined and the other is broad-lined. We present combined interferometric and spectroscopic studies of LamVir. The small scale of the LamVir orbit (~20 mas) is well resolved by the Infrared Optical Telescope Array (IOTA), allowing us to determine its elements as well as the physical properties of the components to high accuracy. The masses of the two stars are determined to be 1.897 Msun and 1.721 Msun, with 0.7% and 1.5% errors respectively, and the two stars are found to have the same temperature of 8280 +/- 200 K. The accurately determined properties of LamVir allow comparisons between observations and current stellar evolution models, and reasonable matches are found. The best-fit stellar model gives LamVir a subsolar metallicity of Z=0.0097, and an age of 935 Myr. The orbital and physical parameters of LamVir also allow us to study its tidal evolution time scales and status. Although currently atomic diffusion is considered to be the most plausible cause of the Am phenomenon, the issue is still being actively debated in the literature. With the present study of the properties and evolutionary status of LamVir, this system is an ideal candidate for further detailed abundance analyses that might shed more light on the source of the chemical anomalies in these A stars.Comment: 43 Pages, 13 figures. Accepted for publication in Ap

Modeling of the Super-Eddington Phase for Classical Novae: Five IUE Novae

We present a light curve model for the super-Eddington luminosity phase of five classical novae observed with IUE. Optical and UV light curves are calculated based on the optically thick wind theory with a reduced effective opacity for a porous atmosphere. Fitting a model light curve with the UV 1455 \AA light curve, we determine the white dwarf mass and distance to be (1.3 M_sun, 4.4 kpc) for V693 CrA, (1.05 M_sun, 1.8 kpc) for V1974 Cyg, (0.95 M_sun, 4.1 kpc) for V1668 Cyg, (1.0 M_sun, 2.1 kpc) for V351 Pup, and (1.0 M_sun, 4.3 kpc) for OS And.Comment: 9 pages including 8 figures, to appear in the Astrophysical Journa

The mass of the neutron star in Vela X-1 and tidally induced non-radial oscillations in GP Vel

We report new radial velocity observations of GP Vel/HD77581, the optical companion to the eclipsing X-ray pulsar Vela X-1. Using data spanning more than two complete orbits of the system, we detect evidence for tidally induced non-radial oscillations on the surface of GP Vel, apparent as peaks in the power spectrum of the residuals to the radial velocity curve fit. By removing the effect of these oscillations (to first order) and binning the radial velocities, we have determined the semi-amplitude of the radial velocity curve of GP Vel to be K_o=22.6+/-1.5 km/s. Given the accurately measured semi-amplitude of the pulsar's orbit, the mass ratio of the system is 0.081+/-0.005. We are able to set upper and lower limits on the masses of the component stars as follows. Assuming GP Vel fills its Roche lobe then the inclination angle of the system, i=70.1+/-2.6 deg. In this case we obtain the masses of the two stars as M_x=2.27 +/-0.17 M_sun for the neutron star and M_o=27.9+/-1.3 M_sun for GP Vel. Conversely, assuming the inclination angle is i=90 deg, the ratio of the radius of GP Vel to the radius of its Roche lobe is beta=0.89+/-0.03 and the masses of the two stars are M_x=1.88+/-0.13 M_sun and M_o=23.1+/-0.2 M_sun. A range of solutions between these two sets of limits is also possible, corresponding to other combinations of i and beta. In addition, we note that if the zero phase of the radial velocity curve is allowed as a free parameter, rather than constrained by the X-ray ephemeris, a significantly improved fit is obtained with an amplitude of 21.2+/-0.7 km/s and a phase shift of 0.033+/-0.007 in true anomaly. The apparent shift in the zero phase of the radial velocity curve may indicate the presence of an additional radial velocity component at the orbital period.Comment: Accepted for publication in Astronomy & Astrophysic

Modelling CO emission from Mira's wind

We have modelled the circumstellar envelope of {\it o} Ceti (Mira) using new observational constraints. These are obtained from photospheric light scattered in near-IR vibrational-rotational lines of circumstellar CO molecules at 4.6 micron: absolute fluxes, the radial dependence of the scattered intensity, and two line ratios. Further observational constraints are provided by ISO observations of far-IR emission lines from highly excited rotational states of the ground vibrational state of CO, and radio observations of lines from rotational levels of low excitation of CO. A code based on the Monte-Carlo technique is used to model the circumstellar line emission. We find that it is possible to model the radio and ISO fluxes, as well as the highly asymmetric radio-line profiles, reasonably well with a spherically symmetric and smooth stellar wind model. However, it is not possible to reproduce the observed NIR line fluxes consistently with a `standard model' of the stellar wind. This is probably due to incorrectly specified conditions of the inner regions of the wind model, since the stellar flux needs to be larger than what is obtained from the standard model at the point of scattering, i.e., the intermediate regions at approximately 100-400 stellar radii (2"-7") away from the star. Thus, the optical depth in the vibrational-rotational lines from the star to the point of scattering has to be decreased. This can be accomplished in several ways. For instance, the gas close to the star (within approximately 2") could be in such a form that light is able to pass through, either due to the medium being clumpy or by the matter being in radial structures (which, further out, developes into more smooth or shell-like structures).Comment: 18 pages, 3 figures, accepted for publication in Ap

Predictors of training-related improvement in visuomotor performance in patients with multiple sclerosis: A behavioural and MRI study

Background: The development of tailored recovery-oriented strategies in multiple sclerosis requires early identification of an individual’s potential for functional recovery. Objective: To identify predictors of visuomotor performance improvements, a proxy of functional recovery, using a predictive statistical model that combines demographic, clinical and magnetic resonance imaging (MRI) data. Methods: Right-handed multiple sclerosis patients underwent baseline disability assessment and MRI of the brain structure, function and vascular health. They subsequently undertook 4 weeks of right upper limb visuomotor practice. Changes in performance with practice were our outcome measure. We identified predictors of improvement in a training set of patients using lasso regression; we calculated the best performing model in a validation set and applied this model to a test set. Results: Patients improved their visuomotor performance with practice. Younger age, better visuomotor abilities, less severe disease burden and concurrent use of preventive treatments predicted improvements. Neuroimaging localised outcome-relevant sensory motor regions, the microstructure and activity of which correlated with performance improvements. Conclusion: Initial characteristics, including age, disease duration, visuo-spatial abilities, hand dexterity, self-evaluated disease impact and the presence of disease-modifying treatments, can predict functional recovery in individual patients, potentially improving their clinical management and stratification in clinical trials. MRI is a correlate of outcome, potentially supporting individual prognosis