Stellar systems consisting of multiple stars tend to undergo tidal interactions when the separations between the stars are short. While tidal phenomena have been extensively studied, a certain tidal effect exclusive to hierarchical triples (triples in which one component star has a much wider orbit than the others) has hardly received any attention, mainly due to its complexity and consequent resistance to being modelled. This tidal effect is the tidal perturbation of the tertiary by the inner binary, which in turn depletes orbital energy from the inner binary, causing the inner binary separation to shrink. In this paper, we develop a fully numerical simulation of these "tertiary tides" by modifying established tidal models. We also provide general insight as to how close a hierarchical triple needs to be in order for such an effect to take place, and demonstrate that our simulations can effectively retrieve the orbital evolution for such systems. We conclude that tertiary tides are a significant factor in the evolution of close hierarchical triples, and strongly influence at least ~ 1% of all multiple star systems.

The observed masses, radii and temperatures of 60 medium- to long-period binaries, most of which contain a cool, evolved star and a hotter less-evolved one, are compared with theoretical models which include (a) core convective overshooting, (b)mass loss, possibly driven by dynamo action as in RS CVn binaries, and (c) tidal friction, including its effect on orbital period through magnetic braking. A reasonable fit is found in about 42 cases, but in 11 other cases the primaries appear to have lost either more mass or less mass than the models predict, and in 4 others the orbit is predicted to be either more or less circular than observed. Of the remaining 3 systems, two ($\gamma$ Per and HR 8242) have a markedly `over-evolved' secondary, our explanation being that the primary component is the merged remnant of a former short-period sub-binary in a former triple system. The last system (V695 Cyg) defies any agreement at present. Mention is also made of three other systems (V643 Ori, OW Gem and V453 Cep), which are relevant to our discussion.

We present the Kepler photometric light-variation analysis of the late-type double-lined binary system V568 Lyr that is in the field of the high metallicity old open cluster NGC 6791. The radial velocity and the high-quality short-cadence light curve of the system are analysed simultaneously. The masses, radii and luminosities of the component stars are $M_1 = 1.0886\pm0.0031\, M{\odot}$, $M_2 = 0.8292 \pm 0.0026\, M{\odot}$, $R_1 = 1.4203\pm 0.0058\, R{\odot}$, $R_2 = 0.7997 \pm 0.0015\, R{\odot}$, $L_1 = 1.85\pm 0.15\, L{\odot}$, $L_2 = 0.292 \pm 0.018\, L{\odot}$ and their separation is $a = 31.060 \pm 0.002\, R{\odot}$. The distance to NGC 6791 is determined to be $4.260\pm 0.290\,$kpc by analysis of this binary system. We fit the components of this well-detached binary system with evolution models made with the Cambridge STARS and TWIN codes to test low-mass binary star evolution. We find a good fit with a metallicity of $Z = 0.04$ and an age of $7.704\,$Gyr. The standard tidal dissipation, included in TWIN is insufficient to arrive at the observed circular orbit unless it formed rather circular to begin with.

As part of the BANANA project (Binaries Are Not Always Neatly Aligned), we have found that the eclipsing binary CV Velorum has misaligned rotation axes. Based on our analysis of the Rossiter-McLaughlin effect, we find sky-projected spin-orbit angles of $\beta_{\rm p} = -52\pm6^{\circ}$ and $\beta_{\rm s}= 3\pm7^{\circ}$ for the primary and secondary stars (B2.5V $+$ B2.5V, $P=6.9$ d). We combine this information with several measurements of changing projected stellar rotation speeds ($v \sin i_{\star}$) over the last $30$ years, leading to a model in which the primary star's obliquity is $\approx65^{\circ}$, and its spin axis precesses around the total angular momentum vector with a period of about $140$ years. The geometry of the secondary star is less clear, although a significant obliquity is also implicated by the observed time variations in the $v \sin i_{\star}$. By integrating the secular tidal evolution equations backward in time, we find that the system could have evolved from a state of even stronger misalignment similar to DI Herculis, a younger but otherwise comparable binary.

We have calculated the equilibrium properties of a star in a circular, equatorial orbit about a Super-Massive Black Hole (SMBH), when the star fills and overflows its Roche lobe. The mass transfer time scale is anticipated to be long compared with the dynamical time and short compared with the thermal time of the star, so that the entropy as a function of the interior mass is conserved. We have studied how the stellar entropy, pressure, radius, mean density, and orbital angular momentum vary when the star is evolved adiabatically, for a representative set of stars. We have shown that the stellar orbits change with the stellar mean density. Therefore, sun-like stars, upper main sequence stars and red giants will spiral inward and then outward with respect to the hole in this stable mass transfer process, while lower main sequence stars, brown dwarfs and white dwarfs will always spiral outward.

The multiplicities of stars, and some other properties, were collected recently by Eggleton & Tokovinin, for the set of 4559 stars with Hipparcos magnitude brighter than 6.0 (4558 excluding the Sun). In this paper I give a numerical recipe for constructing, by a Monte Carlo technique, a theoretical ensemble of multiple stars that resembles the observed sample. Only multiplicities up to 8 are allowed; the observed set contains only multiplicities up to 7. In addition, recipes are suggested for dealing with the selection effects and observational uncertainties that attend the determination of multiplicity. These recipes imply, for example, that to achieve the observed average multiplicity of 1.53, it would be necessary to suppose that the real population has an average multiplicity slightly over 2.0. This numerical model may be useful for (a) comparison with the results of star and star cluster formation theory, (b) population synthesis that does not ignore multiplicity above 2, and (c) initial conditions for dynamical cluster simulations.

Luminous X-ray stars are very often found in visual double or multiple stars. Binaries with periods of a few days possess the highest degree of coronal X-ray activity among regular, non-relativistic stars because of their fast, tidally driven rotation. But the orbital periods in visual double stars are too large for any direct interaction between the companions to take place. We suggest that most of the strongest X-ray components in resolved binaries are yet-undiscovered short-period binaries, and that a few are merged remnants of such binaries. The omnipresence of short-period active stars, e.g. of BY-Dra-type binaries, in multiple systems is explained via the dynamical evolution of triple stars with large mutual inclinations. The dynamical perturbation on the inner pair pumps up the eccentricity in a cyclic manner, a phenomenon known as Kozai cycling. At times of close periapsis, tidal friction reduces the angular momentum of the binary, causing it to shrink. When the orbital period of the inner pair drops to a few days, fast surface rotation of the companions is driven by tidal forces, boosting activity by a few orders of magnitude. If the period drops still further, a merger may take place leaving a rapidly-rotating active dwarf with only a distant companion.

We consider the multiplicity of stellar systems with (combined) magnitude brighter than 6.00 in Hipparcos magnitudes. We identify 4559 such bright systems (including the Sun), and the frequencies of multiplicities 1, 2,..., 7 are found to be 2718, 1437, 285, 86, 20, 11 and 2. We discuss the uncertainties, which are substantial. We also consider the distributions of periods of orbits and sub-orbits. We note that for the even more restricted set of 478 systems with V_H <= 4.00 the proportions of higher multiples up to sextuple are progressively larger (213, 179, 54, 19, 8, 5), suggesting substantial incompleteness in even the reasonably well-studied larger sample. This sample can be seen as relatively thoroughly studied for multiplicity, and reasonably representative of stars more massive than the Sun. But the restriction to V_H <= 6 means that our sample contains hardly any systems where all components are low-mass main-sequence stars (K or M). Data on multiplicity is important as a constraint on (a) the star-formation problem, (b) the problem of the evolution of the Galactic stellar population, and (c) the interaction of dynamics and evolution through the effect of Kozai cycles. We discuss these topics briefly.

Three-dimensional stellar modeling has enabled us to identify a deep-mixing mechanism that must operate in all low mass giants. This mixing process is not optional, and is driven by a molecular weight inversion created by the 3He(3He,2p)4He reaction. In this paper we characterize the behavior of this mixing, and study its impact on the envelope abundances. It not only eliminates the problem of 3He overproduction, reconciling stellar and big bang nucleosynthesis with observations, but solves the discrepancy between observed and calculated CNO isotope ratios in low mass giants, a problem of more than 3 decades' standing. This mixing mechanism, which we call `$\delta\mu$-mixing', operates rapidly (relative to the nuclear timescale of overall evolution, ~ 10^8 yrs) once the hydrogen burning shell approaches the material homogenized by the surface convection zone. In agreement with observations, Pop I stars between 0.8 and 2.0$\Msun$ develop 12C/13C ratios of 14.5 +/- 1.5, while Pop II stars process the carbon to ratios of 4.0 +/- 0.5. In stars less than 1.25$\Msun$, this mechanism also destroys 90% to 95% of the 3He produced on the main sequence.

We continue our study of the core helium flash using the three dimensional hydrodynamics code Djehuty. Continuing from earlier calculations, we now take relaxed 3D configurations and add various amounts of rotation. We find that rotation periods consistent with those observed in white dwarfs produce negligible changes in the structure and evolution of the core flash, at least for the very small timescales we have yet been able to investigate. There is no sign of any extra mixing due to the rotation. There is some inconclusive evidence for a slight change in the luminosity, at the 1% level.

Low-mass stars, ~1-2 solar masses, near the Main Sequence are efficient at producing 3He, which they mix into the convective envelope on the giant branch and distribute into the Galaxy by way of envelope loss. This process is so efficient that it is difficult to reconcile the observed cosmic abundance of 3He with the predictions of Big Bang nucleosynthesis. In this paper we find, by modeling a red giant with a fully three-dimensional hydrodynamic code and a full nucleosynthetic network, that mixing arises in the supposedly stable and radiative zone between the hydrogen-burning shell and the base of the convective envelope. This mixing is due to Rayleigh-Taylor instability within a zone just above the hydrogen-burning shell. In this zone the burning of the 3He left behind by the retreating convective envelope is predominantly by the reaction 3He + 3He -> 4He + 2p, a reaction which, untypically for stellar nuclear reactions, {\it lowers} the mean molecular weight, leading to a local minimum. This local minimum leads to Rayleigh-Taylor instability, and turbulent motion is generated which will continue ultimately up into the normal convective envelope. Consequently material from the envelope is dragged down sufficiently close to the burning shell that the 3He in it is progressively destroyed. Thus we are able to remove the threat that 3He production in low-mass stars poses to the Big Bang nucleosynthesis of 3He. Some slow mixing mechanism has long been suspected, that connects the convective envelope of a red giant to the burning shell. It appears to be necessary to account for progressive changes in the 12C/13C and 14N/12C ratios on the First Giant Branch. We suggest that these phenomena are also due to the Rayleigh-Taylor-unstable character of the 3He-burning region.

Low-mass stars, ~1-2 solar masses, near the Main Sequence are efficient at producing He-3, which they mix into the convective envelope on the giant branch and should distribute into the Galaxy by way of envelope loss. This process is so efficient that it is difficult to reconcile the low observed cosmic abundance of He-3 with the predictions of both stellar and Big Bang nucleosynthesis. In this paper we find, by modeling a red giant with a fully three-dimensional hydrodynamic code and a full nucleosynthetic network, that mixing arises in the supposedly stable and radiative zone between the hydrogen-burning shell and the base of the convective envelope. This mixing is due to Rayleigh-Taylor instability within a zone just above the hydrogen-burning shell, where a nuclear reaction lowers the mean molecular weight slightly. Thus we are able to remove the threat that He-3 production in low-mass stars poses to the Big Bang nucleosynthesis of He-3.

We model the core helium flash in a low-mass red giant using Djehuty, a fully three-dimensional (3D) code. The 3D structures were generated from converged models obtained during the 1D evolutionary calculation of a 1$\Msun$ star. Independently of which starting point we adopted, we found that after some transient relaxation the 3D model settled down with a briskly convecting He-burning shell that was not very different from what the 1D model predicted.

We study the dynamical interaction in which the two single runaway stars AE Aurigae and mu Columbae and the binary iota Orionis acquired their unusually high space velocity. The two single runaways move in almost opposite directions with a velocity greater than 100 km/s away from the Trapezium cluster. The star iota Ori is an eccentric (e=0.8) binary moving with a velocity of about 10 km/s at almost right angles with respect to the two single stars. The kinematic properties of the system suggest that a strong dynamical encounter occurred in the Trapezium cluster about 2.5 Myr ago. Curiously enough, the two binary components have similar spectral type but very different masses, indicating that their ages must be quite different. This observation leads to the hypothesis that an exchange interaction occurred in which an older star was swapped into the original iota Orionis binary. We test this hypothesis by a combination of numerical and theoretical techniques, using N-body simulations to constrain the dynamical encounter, binary evolution calculations to constrain the high orbital eccentricity of iota Orionis and stellar evolution calculations to constrain the age discrepancy of the two binary components. We find that an encounter between two low eccentricity (0.4<e<0.6) binaries with comparable binding energy, leading to an exchange and the ionization of the wider binary, provides a reasonable solution to this problem.

Current practice in stellar evolution is to employ one-dimensional calculations that quantitatively apply only to a minority of the observed stars (single non-rotating stars, or well detached binaries). Even in these systems, astrophysicists are dependent on approximations to handle complex three-dimensional processes like convection. Understanding the structure of binary stars, like those that lead to the Type Ia supernovae used to measure the expansion of the universe, are grossly non-spherical and await a 3D treatment. To approach very large problems like multi-dimensional modeling of stars, the Lawrence Livermore National Laboratory has invested in massively parallel computers and invested even more in developing the algorithms to utilize them on complex physics problems. We have leveraged skills from across the lab to develop a 3D stellar evolution code, Djehuty (after the Egyptian god for writing and calculation) that operates efficiently on platforms with thousands of nodes, with the best available physical data (opacities, EOS, etc.). Djehuty has incorporated all basic physics for modeling stars including an accurate equation of state, radiation transport by diffusion, thermonuclear reaction rates, and hydrodynamics, and we have begun testing it in a number of applications.

This is a summary paper of MODEST-2, a workshop held at the Astronomical Institute ``Anton Pannekoek'' in Amsterdam, 16-17 December 2002. MODEST is a loose collaboration of people interested in MOdelling DEnse STellar systems, particularly those interested in modelling these systems using all the available physics (stellar dynamics, stellar evolution, hydrodynamics and the interplay between the three) by defining interfaces between different codes. In this paper, we summarize 1) the main advances in this endeavour since MODEST-1; 2) the main science goals which can be and should be addressed by these types of simulations; and 3) the most pressing theoretical and modelling advances that we identified.

We discuss the evolution of binaries with moderately high masses (about 10 - 30 solar masses), and with periods of about 3 - 300d, corresponding mostly to early Case B. These are usually thought to evolve either by reasonably conservative Roche-lobe overflow, if the initial mass ratio is fairly mild, or else by highly non-conservative common-envelope evolution, with spiral-in to short periods (hours, typically), if the initial mass ratio is rather extreme. We discuss here a handful of binaries from part of this period range (about 50 - 250d), which appear to have followed a different path: we argue that they must have lost a large proportion of initial mass (about 70 - 80%), but without shortening their periods at all. We suggest that their behaviour may be due to the fact that stars of such masses, when evolved also to rather large radii, are not far from the Humphreys-Davidson limit where single stars lose their envelopes spontaneously in P Cygni winds, and so have envelopes which are only lightly bound to the core. These envelopes therefore may be relatively easily dissipated by the perturbing effect of a companion. In addition, some or all of the stars considered here may have been close to the Cepheid instability strip when they filled their Roche lobes. One or other, or both, of high luminosity and Cepheid instability, in combination with an appropriately close binary companion, may be implicated.

We apply a model of dynamo-driven mass loss, magnetic braking and tidal friction to the evolution of stars with cool convective envelopes; in particular we apply it to binary stars where the combination of magnetic braking and tidal friction can cause angular-momentum loss from the {\it orbit}. For the present we consider the simplification that only one component of a binary is subject to these non-conservative effects, but we emphasise the need in some circumstances to permit such effects in {\it both} components. The model is applied to examples of (i) the Sun, (ii) BY Dra binaries, (iii) Am binaries, (iv) RS CVn binaries, (v) Algols, (vi) post-Algols. A number of problems regarding some of these systems appear to find a natural explanation in our model. There are indications from other systems that some coefficients in our model may vary by a factor of 2 or so from system to system; this may be a result of the chaotic nature of dynamo activity.

The DJEHUTY project is an intensive effort at the Lawrence Livermore National Laboratory (LLNL) to produce a general purpose 3-D stellar structure and evolution code to study dynamic processes in whole stars.

We present equations governing the way in which both the orbit and the intrinsic spins of stars in a close binary should evolve subject to a number of perturbing forces, including the effect of a third body in a possibly inclined wider orbit. We illustrate the solutions in some binary-star and triple-star situations: tidal friction in a wide but eccentric orbit of a radio pulsar about a B star, the Darwin and eccentricity instabilities in a more massive but shorter-period massive X-ray binary, and the interaction of tidal friction with Kozai cycles in a triple such as Algol (beta-Per), at an early stage in that star's life when all 3 components were ZAMS stars. We also attempt to model in some detail the interesting triple system SS Lac, which stopped eclipsing in about 1950. We find that our model of SS Lac is quite constrained by the relatively good observational data of this system, and leads to a specific inclination (29 deg) of the outer orbit relative to the inner orbit at epoch zero (1912). Although the intrinsic spins of the stars have little effect on the orbit, the converse is not true: the spin axes can vary their orientation relative to the close binary by up to 120 deg on a timescale of about a century.

Using Eggleton's stellar evolution code, we carry out 150 runs of Pop I binary evolution calculations, with the initial primary mass between 1 and 8 solar masses the initial mass ratio between 1.1 and 4, and the onset of Roche lobe overflow (RLOF) at an early, middle, or late Hertzsprung-gap stage. We assume that RLOF is conservative in the calculations, and find that the remnant mass of the primary may change by more than 40 per cent over the range of initial mass ratio or orbital period, for a given primary mass. This is contrary to the often-held belief that the remnant mass depends only on the progenitor mass if mass transfer begins in the Hertzsprung gap. We fit a formula, with an error less than 3.6 per cent, for the remnant (white dwarf) mass as a function of the initial mass of the primary, the initial mass ratio, and the radius of the primary at the onset of RLOF. We also find that a carbon-oxygen white dwarf with mass as low as 0.33 solar masses may be formed if the primary's initial mass is around 2.5 solar masses.

We undertake a comparison of observed Algol-type binaries with a library of computed Case A binary evolution tracks. The library consists of 5500 binary tracks with various values of initial primary mass M_{10}, mass ratio q_0, and period P_0, designed to sample the phase-space of Case A binaries in the range -0.10 < log M_{10} < 1.7. Each binary is evolved using a standard code with the assumption that both total mass and orbital angular momentum are conserved. This code follows the evolution of both stars until the point where contact or reverse mass transfer occurs. The resulting binary tracks show a rich variety of behavior which we sort into several subclasses of Case A and Case B. We present the results of this classification, the final mass ratio and the fraction of time spent in Roche Lobe overflow for each binary system. The conservative assumption under which we created this library is expected to hold for a broad range of binaries, where both components have spectra in the range G0 to B1 and luminosity class III - V. We gather a list of relatively well-determined observed hot Algol-type binaries meeting this criterion, as well as a list of cooler Algol-type binaries where we expect significant dynamo-driven mass loss and angular momentum loss. We fit each observed binary to our library of tracks using a chi^2-minimizing procedure. We find that the hot Algols display overall acceptable chi^2, confirming the conservative assumption, while the cool Algols show much less acceptable chi^2 suggesting the need for more free parameters, such as mass and angular momentum loss.

We derive from first principles equations governing (a) the quadrupole tensor of a star distorted by both rotation and the presence of a companion in a possibly eccentric orbit, (b) a functional form for the dissipative force of tidal friction, based on the concept that the rate of energy loss from a time-dependent tide should be a positive-definite function of the rate of change of the quadrupole tensor as seen in the frame which rotates with the star, and (c) the equations governing the rates of change of the magnitude and direction of the stellar rotation, and the orbital period and eccentricity, based on the concept of the Laplace-Runge-Lenz vector. Our analysis leads relatively simply to a closed set of equations, valid for arbitrary inclination of the stellar spin to the orbit. The results are equivalent to classical results based on the rather less clear principle that the tidal bulge lags behind the line of centres by some time determined by the rate of dissipation; our analysis gives the effective lag time as a function of the dissipation rate and the quadrupole moment. We discuss briefly some possible applications of the formulation.

In this paper we study the formation, evolution and disruption of hierarchical systems in open clusters. With this purpose, N-body simulations of star clusters containing an initial population of binaries have been carried out using Aarseth's NBODY4 and NBODY5 codes. Stable triples may form from strong interactions of two binaries in which the widest pair is disrupted. The most frequent type of hierarchical systems found in the cluster models are triples in which the outer star is single, but in some cases the outer body is also a binary, giving a hierarchical quadruple. Many triple systems are non-coplanar and the presence of even a very distant and small outer companion may affect the orbital parameters of the inner binary, including a possible mechanism of significant shrinkage if the binary experiences a weak tidal dissipation. The main features of these systems are analyzed in order to derive general properties which can be checked by observations. The formation of exotic objects, such as blue stragglers and white dwarf binaries, inside hierarchical triple systems is particularly interesting.

About 5-15% of stellar systems are at least triple. About 1% of systems with a primary of $\tgs 1 \Mscun$ are triple with a {\it longer} peri od that is less than 30y, and so may in principle be capable of Roche-lobe overflow in both the inner and the outer orbits, at different times. We discuss possible evolutionary paths for these systems, some of which may lead to objects that are difficult to understand in the context of purely binary evolution. An example is OW Gem, a binary containing two supergiants (spectral types F and G) with masses that difffer by a factor of 1.5. There is also a triple-star pathway which could lead rather naturally to low-mass X-ray binaries; whereas binary pathways often appear rather contrived. We also discuss some dynamical pr ocesses involved in the 3-body problem. A number of triple stars are found in clusters. Similar systems can be created by gravitational capture during N-body simulations of Galactic clusters, especially if there is a n assumed primordial binary population. We discuss the properties of these triples , and note that many can be quite long-lived.

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.