Approximate input physics for stellar modelling

We present a simple and efficient, yet reasonably accurate, equation of state, which at the moderately low temperatures and high densities found in the interiors of stars less massive than the Sun is substantially more accurate than its predecessor by Eggleton, Faulkner & Flannery. Along with the most recently available values in tabular form of opacities, neutrino loss rates, and nuclear reaction rates for a selection of the most important reactions, this provides a convenient package of input physics for stellar modelling. We briefly discuss a few results obtained with the updated stellar evolution code.Comment: uuencoded compressed postscript. The preprint are also available at http://www.ast.cam.ac.uk/preprint/PrePrint.htm

The nature of the progenitor of the Type II-P supernova 1999em

We present high quality ground-based VRI images of the site of the Type II-P SN1999em (in NGC1637) taken before explosion, which were extracted from the CFHT archive. We determine a precise position of the SN on these images to an accuracy of 0.17''. The host galaxy is close enough (7.5 +/- 0.5 Mpc) that the bright supergiants are resolved as individual objects, however we show that there is no detection of an object at the SN position before explosion that could be interpreted as the progenitor star. By determining the sensitivity limits of the VRI data, we derive bolometric luminosity limits for the progenitor. Comparing these to standard stellar evolutionary tracks which trace evolution up to the point of core carbon ignition, we initially derive an upper mass limit of approximately 12M_sol. However we present evolutionary calculations that follow 7-12M_sol stars throughout their C-burning lifetime and show that we can restrict the mass of the progenitor even further. Our calculations indicate that progenitors initially of 8-10M_sol, undergoing expected mass loss, can also be excluded because a second dredge up sends them to somewhat higher luminosities than a star of initially 12M_sol. These results limit the progenitor's initial main-sequence mass to a very narrow range of 12 +/- 1 M_sol. We discuss the similarities between the Type II-P SNe 1999em and 1999gi and their progenitor mass limits, and suggest that SN Type II-P originate only in intermediate mass stars of 8-12M_sol, which are in the red supergiant region and that higher mass stars produce the other Type II sub-types. (Abridged).Comment: Replaced with accepted version to appear in ApJ, 30 pages, inc. 6 figure

The Zero Age Main Sequence of WIMP burners

We modify a stellar structure code to estimate the effect upon the main sequence of the accretion of weakly interacting dark matter onto stars and its subsequent annihilation. The effect upon the stars depends upon whether the energy generation rate from dark matter annihilation is large enough to shut off the nuclear burning in the star. Main sequence WIMP burners look much like protostars moving on the Hayashi track, although they are in principle completely stable. We make some brief comments about where such stars could be found, how they might be observed and more detailed simulations which are currently in progress. Finally we comment on whether or not it is possible to link the paradoxically young OB stars found at the galactic centre with WIMP burners.Comment: 4 pages, 3 figs. Matches published versio

Cuspons, peakons and regular gap solitons between three dispersion curves

A general wave model with the cubic nonlinearity is introduced to describe a situation when the linear dispersion relation has three branches, which would intersect in the absence of linear couplings between the three waves. Actually, the system contains two waves with a strong linear coupling between them, to which a third wave is then coupled. This model has two gaps in its linear spectrum. Realizations of this model can be made in terms of temporal or spatial evolution of optical fields in, respectively, a planar waveguide or a bulk-layered medium resembling a photonic-crystal fiber. Another physical system described by the same model is a set of three internal wave modes in a density-stratified fluid. A nonlinear analysis is performed for solitons which have zero velocity in the reference frame in which the group velocity of the third wave vanishes. Disregarding the self-phase modulation (SPM) term in the equation for the third wave, we find two coexisting families of solitons: regular ones, which may be regarded as a smooth deformation of the usual gap solitons in a two-wave system, and cuspons with a singularity in the first derivative at their center. Even in the limit when the linear coupling of the third wave to the first two vanishes, the soliton family remains drastically different from that in the linearly uncoupled system; in this limit, regular solitons whose amplitude exceeds a certain critical value are replaced by peakons. While the regular solitons, cuspons, and peakons are found in an exact analytical form, their stability is tested numerically, which shows that they all may be stable. If the SPM terms are retained, we find that there again coexist two different families of generic stable soliton solutions, namely, regular ones and peakons.Comment: a latex file with the text and 10 pdf files with figures. Physical Review E, in pres

Binary coalescence from case A evolution -- mergers and blue stragglers

We constructed some main-sequence mergers from case A binary evolution and studied their characteristics via Eggleton's stellar evolution code. Both total mass and orbital angular momentum are conservative in our binary evolutions. Some mergers might be on the left of the ZAMS as defined by normal surface composition on a CMD because of enhanced surface helium content. The study also shows that central hydrogen content of the mergers is independent of mass. As a consequence, we fit the formula of magnitude and B-V of the mergers when they return back to thermal equilibrium with maximum error 0.29 and 0.037, respectively. Employing the consequences above, we performed Monte Carlo simulations to examine our models in NGC 2682 and NGC 2660. In NGC 2682, binary mergers from our models cover the region with high luminosity, but its importance is much less than that of AML. Our results are well-matched to the observations of NGC2660 if there is about 0.5Mo of mass loss in the merger process.Comment: 14 pages, 12 figures. accepted by MNRA

The Evolutionary Status of SS433

We consider possible evolutionary models for SS 433. We assume that common-envelope evolution is avoided if radiation pressure is able to expel most of a super-Eddington accretion flow from a region smaller than the accretor's Roche lobe. This condition is satisfied, at least initially, for largely radiative donors with masses in the range 4-12 solar masses. For donors more massive than about 5 solar masses, moderate mass ratios q = M_2/M_1 > 1 are indicated, thus tending to favor black-hole accretors. For lower mass donors, evolutionary considerations do not distinguish between a neutron star or black hole accretor. In all cases the mass transfer (and mass loss) rates are much larger than the likely mass-loss rate in the precessing jets. Almost all of the transferred mass is expelled at radii considerably larger than the jet acceleration region, producing the "stationary" H-alpha line, the infrared luminosity, and accounting for the low X-ray luminosity.Comment: 13 pages, Astrophysical Journal Letters, accepte

Alignment Timescale of the Microquasar GRO J1655-40

The microquasar GRO J1655-40 has a black hole with spin angular momentum apparently misaligned to the orbital plane of its companion star. We analytically model the system with a steady state disc warped by Lense-Thirring precession and find the timescale for the alignment of the black hole with the binary orbit. We make detailed stellar evolution models so as to estimate the accretion rate and the lifetime of the system in this state. The secondary can be evolving at the end of the main sequence or across the Hertzsprung gap. The mass-transfer rate is typically fifty times higher in the latter case but we find that, in both cases, the lifetime of the mass transfer state is at most a few times the alignment timescale. The fact that the black hole has not yet aligned with the orbital plane is therefore consistent with either model. We conclude that the system may or may not have been counter-aligned after its supernova kick but that it is most likely to be close to alignment rather than counteralignment now.Comment: Accepted for publication in MNRA

Global dynamics and stability limits for planetary systems around HD 12661, HD 38529, HD 37124 and HD 160691

In order to distinguish between regular and chaotic planetary orbits we apply a new technique called MEGNO in a wide neighbourhood of orbital parameters determined using standard two-body Keplerian fits for HD 12661, HD 38529, HD 37124 and HD 160691 planetary systems. We show that the currently announced orbital parameters place these systems in very different situations from the point of view of dynamical stability. While HD 38529 and HD 37124 are located within large stability zones in the phase space around their determined orbits, the preliminary orbits in HD 160691 are highly unstable. The orbital parameters of the HD 12661 planets are located in a border region between stable and unstable dynamical regimes, so while its currently determined orbital parameters produce stable regular orbits, a minor change within the margin of error of just one parameter may result in a chaotic dynamical system.Comment: 12 pages, 3 figures, accepted ApJ, revised version following the referee's repor

Planet Consumption and Stellar Metallicity Enhancements

The evolution of a giant planet within the stellar envelope of a main-sequence star is investigated as a possible mechanism for enhancing the stellar metallicities of the parent stars of extrasolar planetary systems. Three-dimensional hydrodynamical simulations of a planet subject to impacting stellar matter indicate that the envelope of a Jupiter-like giant planet can be completely stripped in the outer stellar convection zone of a solar-mass star. In contrast, Jupiter-like and less massive Saturn-like giant planets are able to survive through the base of the convection zone of a 1.22 solar-mass star. Although strongly dependent on details of planetary interior models, partial or total dissolution of giant planets can result in significant enhancements in the metallicity of host stars with masses between about 1.0 and 1.3 solar masses. The implications of these results with regard to planetary orbital migration are briefly discussed.Comment: 11 pages, 2 figures, accepted for ApJ Letter

Theory of Nonlinear Dispersive Waves and Selection of the Ground State

A theory of time dependent nonlinear dispersive equations of the Schroedinger / Gross-Pitaevskii and Hartree type is developed. The short, intermediate and large time behavior is found, by deriving nonlinear Master equations (NLME), governing the evolution of the mode powers, and by a novel multi-time scale analysis of these equations. The scattering theory is developed and coherent resonance phenomena and associated lifetimes are derived. Applications include BEC large time dynamics and nonlinear optical systems. The theory reveals a nonlinear transition phenomenon, ``selection of the ground state'', and NLME predicts the decay of excited state, with half its energy transferred to the ground state and half to radiation modes. Our results predict the recent experimental observations of Mandelik et. al. in nonlinear optical waveguides