The Hubble Space Telescope (HST) UV Legacy Survey of Galactic Globular Clusters (GCs) has investigated multiple stellar populations by means of the "chromosome map" (ChM) diagnostic tool that maximises the separation between stars with different chemical composition. One of the most challenging features revealed by ChMs analysis is the apparent inhomogeneity among stars belonging to the first population, a phenomenon largely attributed to He variations. However, this explanation is not supported by the uniformity in p-capture elements of these stars. The HST survey has revealed that the GC NGC 3201 shows an exceptionally wide coverage in the Delta(F275W,F814W) parameter of the ChM. We present a chemical abundance analysis of 24 elements in 18 giants belonging to the first population of this GC, and having a wide range in Delta(F275W,F814W). As far as the p-capture elements are concerned, the chemical abundances are typical of 1G stars, as expected from the location of our targets in the ChM. Based on radial velocities and chemical abundances arguments, we find that the three stars with the lowest Delta(F275W,F814W) values are binary candidates. This suggests that, at least those stars could be explained with binarity. These results are consistent with evidence inferred from multi-band photometry that evolved blue stragglers populate the bluest part of the 1G sequence in the ChM. The remaining 15 spectroscopic targets show a small range in the overall metallicity by ~0.10 dex, with stars at higher Delta(F275W,F814W) values having higher absolute abundances. We suggest that a small variation in metals and binarity govern the color spread of the 1G in the ChM, and that evolved blue stragglers contribute to the bluest tail of the 1G sequence.

The cluster M54 lies at the centre of the Sagittarius dwarf spheroidal galaxy, and therefore may be the closest example of a nuclear star cluster. Either in-situ star formation, inspiralling globular clusters, or a combination have been invoked to explain the wide variety of stellar sub-populations in nuclear star clusters. Globular clusters are known to exhibit light element variations, which can be identified using the photometric construct called a chromosome map. In this letter, we create chromosome maps for three distinct age-metallicity sub-populations in the vicinity of M54. We find that the old, metal-poor population shows the signature of light element variations, while the young and intermediate-age metal rich populations do not. We conclude that the nucleus of Sagittarius formed through a combination of in-situ star formation and globular cluster accretion. This letter demonstrates that properly constructed chromosome maps of iron-complex globular clusters can provide insight into the formation locations of the different stellar populations.

Recent studies show that the majority of blue straggler stars (BSSs) in old open clusters are formed through mass transfer from an evolved star onto a main-sequence companion, resulting in a BSS and white dwarf (WD) in a binary system. We present constraints on the mass transfer histories of two BSS-WD binaries in the open cluster NGC 188, using WD temperatures and surface gravities measured with HST COS far-ultraviolet spectroscopy. Adopting a Gaia-based cluster distance of $1847\pm107$ pc, we determine that one system, WOCS 4540, formed through Case C mass transfer resulting in a CO-core white dwarf with $T_{\text{eff}}=17000^{+140}_{-200}$ K and a $\log g=7.80^{+0.06}_{-0.06}$, corresponding to a mass of $0.53^{+0.03}_{-0.03}$ $M_{\odot}$ and a cooling age of $105^{+6}_{-5}$ Myr. The other system, WOCS 5379, formed through Case B mass transfer resulting in a He-core white dwarf with $T_{\text{eff}}=15500^{+170}_{-150}$ K and a $\log g=7.50^{+0.06}_{-0.05}$, corresponding to a mass of $0.42^{+0.02}_{-0.02}$ $M_{\odot}$ and an age of $250^{+20}_{-20}$ Myr. The WD parameters are consistent across four different cluster distance assumptions. We determine possible progenitor binary systems with a grid of accretion models using MESA, and investigate whether these systems would lead to stable or unstable mass transfer. WOCS 4540 likely resulted from stable mass transfer during periastron passage in an eccentric binary system, while WOCS 5379 challenges our current understanding of the expected outcomes for mass transfer from red giant branch stars. Both systems are examples of the value in using detailed analyses to fine-tune our physical understanding of binary evolutionary processes.

We present a new approach to understanding star-to-star helium abundance variations within globular clusters. We begin with detailed radiation hydrodynamics simulations of cluster formation within giant molecular clouds, and investigate the conditions under which multiple populations could be created. Chemical enrichment occurs dynamically as the cluster is assembled. We test two extreme mechanisms for injection of enriched gas within the clusters, and find that realistic multiple populations can be formed in both mechanisms. The stochastic cluster formation histories are dictated by the inherent randomness of the timing and location of the formation of small clusters, which rapidly merge to build up the larger cluster, in combination with continual accretion of gas from the cloud. These cluster formation histories naturally produce a diversity of abundance patterns across the massive cluster population. We conclude that multiple populations are a natural outcome of the typical mode of star cluster formation.

We investigate the possibility that multiple populations in globular clusters arise as a natural by-product of massive star-cluster formation. We use 3D radiative hydrodynamics simulations for the formation of young massive clusters to track their chemical self-enrichment during their first 5 Myr. These clusters form embedded within filamentary Giant Molecular Clouds by a combination of gas accretion and rapid merging of protoclusters. Chemical enrichment is a dynamic process happening as the young cluster assembles, so that the original (1P) and enriched (2P) subpopulations of stars form almost simultaneously. Here we test two simple and opposite extremes for the injection of enriched material into the intracluster gas: we assume either continuous injection in a way that tracks the star formation rate; or sudden injection by a single instantaneous event. Using helium abundance as a proxy for the enrichment, we find that realistic multiple population features can be reproduced by injecting a total helium mass amounting to a few percent of the cluster's total mass. The differences in individual growth histories can lead to widely differing 1P/2P outcomes. These models suggest that dual or multiple populations can emerge rapidly in massive star clusters undergoing the typical mode of star cluster formation.

In this paper we report on the discovery of a double blue straggler star (BSS) sequence in the core of the core-collapsed cluster M15 (NGC 7078). We performed a detailed photometric analysis of the extremely dense core of the cluster using a set of images secured with the Advanced Camera for Survey in the High Resolution Channel mode on-board the Hubble Space Telescope. The proper combination of the large number of single frames in the near-UV (F220W), and blue (F435W) filters allowed us to perform a superb modeling of the Point Spread Function and an accurate deblending procedure. The Color-Magnitude diagram revealed the presence of two distinct parallel sequences of blue stragglers. In particular, the blue BSS sequence is characterized by the intriguing presence of two different branches. The first branch appears extremely narrow, it extends up to 2.5 magnitudes brighter than the cluster main-sequence turnoff (MS-TO) point, and it is nicely reproduced by a 2 Gyr-old collisional isochrone. The second branch extends up to 1.5 magnitudes from the MS-TO and it is reproduced by a 5.5 Gyr-old collisional isochrone. Our observations suggest that each of these branches is mainly constituted by a population of nearly coeval collisional BSS of different masses generated during two episodes of high collisional activity. We discuss the possibility that the oldest episode corresponds to the core-collapse event (occurred about 5.5 Gyr ago), while the most recent one (occurred about 2 Gyr ago) is associated with a core oscillation in the post-core collapse evolution. The discovery of these features provides further strong evidence in support of the connection between the BSS properties and GC dynamical evolution, and it opens new perspectives on the study of core-collapse and post core-collapse evolution.

Blue straggler stars have been proposeed as powerful indicators to measure the dynamical state of Galactic globular clusters. Here we examine for the first time if this framework of blue straggler stars as dynamical clocks, which was specifically developed for Galactic globular clusters, may also hold for younger globular clusters in the Large Magellanic Cloud. Using Hubble Space Telescope observations, we study seven Large Magellanic Cloud star clusters with ages between ~700 Myr and 7 Gyr. We confirm that our sample clusters are consistent with the empirical correlation previously derived for Galactic globular clusters, within a modest tolerance. We suggest that to further examine if blue straggler stars can measure the dynamical state of their host clusters over an extended range of dynamical ages, more studies of dynamically older Magellanic Cloud clusters are required. We discuss the physical implications of our results in terms of their central, dimensionless King potential, as well as the initial retention fraction of black holes.

Blue stragglers and other mass transfer/collision products are likely born with rapid rotation rates due to angular momentum transfer during mass-transfer, merger or collisional formation. However, less is known about the angular momentum evolution of these stars as they age. Here we compare rotation rates and post-formation ages of mass-transfer products to models of angular momentum evolution for normal main-sequence stars and collisionally formed blue stragglers. In our sample, we include both F- and G-type blue stragglers in the cluster NGC 188 and post-mass-transfer GK main-sequence stars in the field, all binaries with WD companions. We compare ages derived from WD cooling models to photometric rotation periods and/or spectral vsini measurements. We demonstrate that these systems have rapid rotation rates soon after formation. They then spin down as they age much as standard solar-type main-sequence stars do. We discuss the physical implications of this result, which suggests that the spin-down of post-mass transfer stars can be described by standard magnetic-braking prescriptions. This opens up the possibility of using gyrochronology as a method to determine the time since formation of blue straggler stars and other post-mass-transfer binaries.

The Gaia mission has opened a new window into the internal kinematics of young star clusters at the sub-km/s level, with implications for our understanding of how star clusters form and evolve. We use a sample of 28 clusters and associations with ages from 1-5 Myr, where lists of members are available from previous X-ray, optical, and infrared studies. Proper motions from Gaia DR2 reveals that at least 75% of these systems are expanding; however, rotation is only detected in one system. Typical expansion velocities are on the order of ~0.5 km/s, and, in several systems, there is a positive radial gradient in expansion velocity. Systems that are still embedded in molecular clouds are less likely to be expanding than those that are partially or fully revealed. One-dimensional velocity dispersions, which range from 1 to 3 km/s, imply that most of the stellar systems in our sample are supervirial and that some are unbound. In star-forming regions that contain multiple clusters or subclusters, we find no evidence that these groups are coalescing, implying that hierarchical cluster assembly, if it occurs, must happen rapidly during the embedded stage.

For a sample of 38 Galactic globular clusters (GCs), we confront the observed distributions of blue straggler (BS) proper motions and masses (derived from isochrone fitting) from the BS catalog of Simunovic & Puzia with theoretical predictions for each of the two main competing BS formation mechanisms. These are mass transfer from an evolved donor on to a main-sequence (MS) star in a close binary system, and direct collisions involving MS stars during binary encounters. We use the \texttt{FEWBODY} code to perform simulations of single-binary and binary-binary interactions. This provides collisional velocity and mass distributions for comparison to the observed distributions. Most clusters are consistent with BSs derived from a dynamically relaxed population, supportive of the binary mass-transfer scenario. In a few clusters, including all the post-core collapse clusters in our sample, the collisional velocities provide the best fit.

We incorporate a semi-analytic formula for the main sequence lifetime of helium-rich stars in N-body simulations of multiple population globular clusters to investigate how the enriched helium stars impact the dynamics of globular clusters. We show that a globular cluster with a helium-rich concentrated population will be slightly smaller than a globular cluster with a normal-helium second generation, with the largest difference seen in the extended normal-helium population. This effect is shown both for a cluster in isolation and one in a realistic Milky Way tidal field. We show that this effect is a result of mass segregation, and the earlier loss of more massive stars in a helium-rich concentrated population due to their decreased main sequence lifetime. The two populations will therefore become dynamically mixed at a slightly earlier time than if they have the same helium abundance. Furthermore, we find that it is possible for the helium-enriched population to become more extended than the normal-helium population if it forms with a low initial concentration or the cluster is able to evolve for a large number of relaxation times. We conclude that the dynamical effects of helium on stellar masses are modest, and that the initial concentration of the two populations and the strength of the Milky Way tidal field are more important in determining the relative radial distributions of multiple populations.

Line-of-sight kinematic studies indicate that many Galactic globular clusters have a significant degree of internal rotation. However, three-dimensional kinematics from a combination of proper motions and line-of-sight velocities are needed to unveil the role of angular momentum in the formation and evolution of these old stellar systems. Here we present the first quantitative study of internal rotation on the plane-of-the-sky for a large sample of globular clusters using proper motions from Gaia DR2. We detect signatures of rotation in the tangential component of proper motions for 11 out of 51 clusters at a $>$3-sigma confidence level, confirming the detection reported in Gaia collaboration et al. (2018) for 8 clusters, and additionally identify 11 GCs with a 2-sigma rotation detection. For the clusters with a detected global rotation, we construct the two-dimensional rotation maps and proper motion rotation curves, and we assess the relevance of rotation with respect to random motions ($V/\sigma\sim0.08-0.51$). We find evidence of a correlation between the degree of internal rotation and relaxation time, highlighting the importance of long-term dynamical evolution in shaping the clusters current properties. This is a strong indication that angular momentum must have played a fundamental role in the earliest phases of cluster formation. Finally, exploiting the spatial information of the rotation maps and a comparison with line-of-sight data, we provide an estimate of the inclination of the rotation axis for a subset of 8 clusters. Our work demonstrates the potential of Gaia data for internal kinematic studies of globular clusters and provides the first step to reconstruct their intrinsic three-dimensional structure.

We present simulations of the dynamical evolution of young embedded star clusters. Our initial conditions are directly derived from X-ray, infrared, and radio observations of local systems, and our models evolve both gas and stars simultaneously. Our regions begin with both clustered and extended distributions of stars, and a gas distribution which can include a filamentary structure in addition to gas surrounding the stellar subclusters. We find that the regions become spherical, monolithic, and smooth quite quickly, and that the dynamical evolution is dominated by the gravitational interactions between the stars. In the absence of stellar feedback, the gas moves gently out of the centre of our regions but does not have a significant impact on the motions of the stars at the earliest stages of cluster formation. Our models at later times are consistent with observations of similar regions in the local neighbourhood. We conclude that the evolution of young proto-star clusters is relatively insensitive to reasonable choices of initial conditions. Models with more realism, such as an initial population of binary and multiple stars and ongoing star formation, are the next step needed to confirm these findings.

Dynamical evolution drives globular clusters toward core collapse, which strongly shapes their internal properties. Diagnostics of core collapse have so far been based on photometry only, namely on the study of the concentration of the density profiles. Here we present a new method to robustly identify core-collapsed clusters based on the study of their stellar kinematics. We introduce the \textit{kinematic concentration} parameter, $c_k$, the ratio between the global and local degree of energy equipartition reached by a cluster, and show through extensive direct $N$-body simulations that clusters approaching core collapse and in the post-core collapse phase are strictly characterized by $c_k>1$. The kinematic concentration provides a suitable diagnostic to identify core-collapsed clusters, independent from any other previous methods based on photometry. We also explore the effects of incomplete radial and stellar mass coverage on the calculation of $c_k$ and find that our method can be applied to state-of-art kinematic datasets.

The young cluster NGC 6231 (stellar ages ~2-7 Myr) is observed shortly after star-formation activity has ceased. Using the catalog of 2148 probable cluster members obtained from Chandra, VVV, and optical surveys (Paper I), we examine the cluster's spatial structure and dynamical state. The spatial distribution of stars is remarkably well fit by an isothermal sphere with moderate elongation, while other commonly used models like Plummer spheres, multivariate normal distributions, or power-law models are poor fits. The cluster has a core radius of $1.2\pm0.1$ pc and a central density of ~200 stars pc$^{-3}$. The distribution of stars is mildly mass segregated. However, there is no radial stratification of the stars by age. Although most of the stars belong to a single cluster, a small subcluster of stars is found superimposed on the main cluster, and there are clumpy non-isotropic distributions of stars outside ~4 core radii. When the size, mass, and age of NGC 6231 are compared to other young star clusters and subclusters in nearby active star-forming regions, it lies at the high-mass end of the distribution but along the same trend line. This could result from similar formation processes, possibly hierarchical cluster assembly. We argue that NGC 6231 has expanded from its initial size but that it remains gravitationally bound.

Galactic globular clusters (GCs) are believed to have formed in-situ in the Galaxy as well as in dwarf galaxies later accreted onto the Milky Way. However, to date, there is no unambiguous signature to distinguish accreted GCs. Using specifically designed $N$-body simulations of GCs evolving in a variety of time-dependent tidal fields (describing the potential of a dwarf galaxy-Milky Way merger), we analyze the effects imprinted to the internal kinematics of an accreted GC. In particular, we look at the evolution of the velocity anisotropy. Our simulations show that at early phases, the velocity anisotropy is determined by the tidal field of the dwarf galaxy and subsequently the clusters will adapt to the new tidal environment, losing any signature of their original environment in a few relaxation times. At 10 Gyr, GCs exhibit a variety of velocity anisotropy profiles, namely, isotropic velocity distribution in the inner regions and either isotropy or radial/tangential anisotropy in the intermediate and outer regions. Independently of an accreted origin, the velocity anisotropy primarily depends on the strength of the tidal field cumulatively experienced by a cluster. Tangentially anisotropic clusters correspond to systems that have experienced stronger tidal fields and are characterized by higher tidal filling factor, $r_{50}/r_j\gtrsim0.17$, higher mass loss $\gtrsim60\%$ and relaxation times $t_{rel}\lesssim10^9$ Gyr. Interestingly, we demonstrate that the presence of tidal tails can significantly contaminate the measurements of velocity anisotropy when a cluster is observed in projection. Our characterization of the velocity anisotropy profiles in different tidal environments provides a theoretical benchmark for the interpretation of the unprecedented amount of three-dimensional kinematic data progressively available for Galactic GCs.

The age distributions of stellar cluster populations have long been proposed to probe the recent formation history of the host galaxy. However, progress is hampered by the limited understanding of cluster disruption by evaporation and tidal shocks. We study the age distributions of clusters in smoothed particle hydrodynamics simulations of isolated disc galaxies, which include a self-consistent, physical model for the formation and dynamical evolution of the cluster population and account for the variation of cluster disruption in time and space. We show that the downward slope of the cluster age distribution due to disruption cannot be reproduced with a single functional form, because the disruption rate exhibits systematic trends with cluster age (the `cruel cradle effect'). This problem is resolved by using the median cluster age to trace cluster disruption. Across 120 independent galaxy snapshots and simulated cluster populations, we perform two-dimensional power law fits of the median cluster age to various macroscopic physical quantities and find that it scales as $t_{\rm med}\propto \Sigma^{-0.51\pm0.03}\sigma_{\rm 1D}^{-0.85\pm0.10}M_{\rm min}^\gamma$, for the gas surface density $\Sigma$, gas velocity dispersion $\sigma_{\rm 1D}$, and minimum cluster mass $M_{\rm min}$. This scaling accurately describes observed cluster populations and indicates disruption by impulsive tidal shocks from the interstellar medium. The term $M_{\rm min}^\gamma$ provides a model-independent way to measure the mass dependence of the cluster disruption time $\gamma$. Finally, the ensemble-average cluster lifetime depends on the gas density less strongly than the instantaneous disruption time of single clusters. These results reflect the variation of cluster disruption in time and space. We provide quantitative ways of accounting for these physics in cluster population studies.

Sub-subgiants are a new class of stars that are optically redder than normal main-sequence stars and fainter than normal subgiant stars. Sub-subgiants, and the possibly related red stragglers (which fall to the red of the giant branch), occupy a region of the color-magnitude diagram that is predicted to be devoid of stars by standard stellar evolution theory. In previous papers we presented the observed demographics of these sources and defined possible theoretical formation channels through isolated binary evolution, the rapid stripping of a subgiant's envelope, and stellar collisions. Sub-subgiants offer key tests for single- and binary-star evolution and stellar collision models. In this paper, we synthesize these findings to discuss the formation frequencies through each of these channels. The empirical data, our analytic formation rate calculations, and analyses of sub-subgiants in a large grid of Monte Carlo globular cluster models suggest that the binary evolution channels may be the most prevalent, though all channels appear to be viable routes to sub-subgiant creation (especially in higher-mass globular clusters). Multiple formation channels may operate simultaneously to produce the observed sub-subgiant population. Finally, many of these formation pathways can produce stars in both the sub-subgiant and red straggler (and blue straggler) regions of the color-magnitude diagram, in some cases as different stages along the same evolutionary sequence.

The internal dynamics of globular clusters (GCs) is strongly affected by two-body interactions that bring the systems to a state of partial energy equipartition. Using a set of Monte Carlo clusters simulations, we investigate the role of the onset of energy equipartition in shaping the mass-to-light ratio (M/L) in GCs. Our simulations show that the M/L profiles cannot be considered constant and their specific shape strongly depends on the dynamical age of the clusters. Dynamically younger clusters display a central peak up to M/L $\simeq25$ $M_\odot/L_\odot$ caused by the retention of dark remnants; this peak flattens out for dynamically older clusters. Moreover, we find that also the global values of M/L correlate with the dynamical state of a cluster quantified as either the number of relaxation times a system has experienced $n_{rel}$ or the equipartition parameter $m_{eq}$: clusters closer to full equipartition (higher $n_{rel}$ or lower $m_{eq}$) display a lower M/L. We show that the decrease of M/L is primarily driven by the dynamical ejection of dark remnants, rather than by the escape of low-mass stars. The predictions of our models are in good agreement with observations of GCs in the Milky Way and M31, indicating that differences in relaxation state alone can explain variations of M/L up to a factor of $\simeq3$. Our characterization of the M/L as a function of relaxation state is of primary relevance for the application and interpretation of dynamical models.

Sub-subgiants are stars observed to be redder than normal main-sequence stars and fainter than normal subgiant (and giant) stars in an optical color-magnitude diagram. The red straggler stars, which lie redward of the red giant branch, may be related and are often grouped together with the sub-subgiants in the literature. These stars defy our standard theory of single-star evolution, and are important tests for binary evolution and stellar collision models. In total, we identify 65 sub-subgiants and red stragglers in 16 open and globular star clusters from the literature; 50 of these, including 43 sub-subgiants, pass our strict membership selection criteria (though the remaining sources may also be cluster members). In addition to their unique location on the color-magnitude diagram, we find that at least 58% (25/43) of sub-subgiants in this sample are X-ray sources with typical 0.5-2.5 keV luminosities of order 10^30 - 10^31 erg/s. Their X-ray luminosities and optical-to-X-ray flux ratios are similar to those of RS CVn active binaries. At least 65% (28/43) of the sub-subgiants in our sample are variables, 21 of which are known to be radial-velocity binaries. Typical variability periods are <15 days. At least 33% (14/43) of the sub-subgiants are H-alpha emitters. These observational demographics provide strong evidence that binarity is important for sub-subgiant formation. Finally, we find that the number of sub-subgiants per unit mass increases toward lower-mass clusters, such that the open clusters in our sample have the highest specific frequencies of sub-subgiants.

We investigate the shallow increase in globular cluster half-light radii with projected galactocentric distance $R_{gc}$ observed in the giant galaxies M87, NGC 1399, and NGC 5128. To model the trend in each galaxy, we explore the effects of orbital anisotropy and tidally under-filling clusters. While a strong degeneracy exists between the two parameters, we use kinematic studies to help constrain the distance $R_\beta$ beyond which cluster orbits become anisotropic, as well as the distance $R_{f\alpha}$ beyond which clusters are tidally under-filling. For M87 we find $R_\beta > 27$ kpc and $20 < R_{f\alpha} < 40$ kpc and for NGC 1399 $R_\beta > 13$ kpc and $10 < R_{f\alpha} < 30$ kpc. The connection of $R_{f\alpha}$ with each galaxy's mass profile indicates the relationship between size and $R_{gc}$ may be imposed at formation, with only inner clusters being tidally affected. The best fitted models suggest the dynamical histories of brightest cluster galaxies yield similar present-day distributions of cluster properties. For NGC 5128, the central giant in a small galaxy group, we find $R_\beta > 5$ kpc and $R_{f\alpha} > 30$ kpc. While we cannot rule out a dependence on $R_{gc}$, NGC 5128 is well fitted by a tidally filling cluster population with an isotropic distribution of orbits, suggesting it may have formed via an initial fast accretion phase. Perturbations from the surrounding environment may also affect a galaxy's orbital anisotropy profile, as outer clusters in M87 and NGC 1399 have primarily radial orbits while outer NGC 5128 clusters remain isotropic.

Stars form predominantly in groups usually denoted as clusters or associations. The observed stellar groups display a broad spectrum of masses, sizes and other properties, so it is often assumed that there is no underlying structure in this diversity. Here we show that the assumption of an unstructured multitude of cluster or association types might be misleading. Current data compilations of clusters show correlations between cluster mass, size, age, maximum stellar mass etc. In this first paper we take a closer look at the correlation of cluster mass and radius. We use literature data to explore relations in cluster and molecular core properties in the solar neighborhood. We show that for embedded clusters in the solar neighborhood there exists a clear correlation between cluster mass and half-mass radius of the form $M_c = C R_c^{\gamma}$ with gamma = 1.7 +/-0.2. This correlation holds for infra red K band data as well as X-ray sources and for clusters containing a hundred stars up to those consisting of a few tens of thousands of stars. The correlation is difficult to verify for clusters containing <30 stars due to low-number statistics. Dense clumps of gas are the progenitors of the embedded clusters. We find a similar slope for the mass-size relation of dense, massive clumps as for the embedded star clusters. This might point at a direct translation from gas to stellar mass: however, it is difficult to relate size measurements for clusters (stars) to those for gas profiles. Taking into account multiple paths for clump mass into cluster mass, we obtain an average star-formation efficiency of 18%{+9.3}{-5.7} for the embedded clusters in the solar neighborhood. The derived mass-radius relation gives constraints for the theory of clustered star formation. Analytical models and simulations of clustered star formation have to reproduce this relation in order to be realistic (abridged)

We perform $N$-body simulations of star clusters in time-dependant galactic potentials. Since the Milky Way was built-up through mergers with dwarf galaxies, its globular cluster population is made up of clusters formed both during the initial collapse of the Galaxy and in dwarf galaxies that were later accreted. Throughout a dwarf-Milky Way merger, dwarf galaxy clusters are subject to a changing galactic potential. Building on our previous work, we investigate how this changing galactic potential affects the evolution of a cluster's half mass radius. In particular, we simulate clusters on circular orbits around a dwarf galaxy that either falls into the Milky Way or evaporates as it orbits the Milky Way. We find that the dynamical evolution of a star cluster is determined by whichever galaxy has the strongest tidal field at the position of the cluster. Thus, clusters entering the Milky Way undergo changes in size as the Milky Way tidal field becomes stronger and that of the dwarf diminishes. We find that ultimately accreted clusters quickly become the same size as a cluster born in the Milky Way on the same orbit. Assuming their initial sizes are similar, clusters born in the Galaxy and those that are accreted cannot be separated based on their current size alone.

We analyze the resolved stellar populations of the faint stellar system, Crater, based on deep optical imaging taken with the Hubble Space Telescope. The HST/ACS-based color-magnitude diagram (CMD) of Crater extends $\sim$4 magnitudes below the oldest main sequence turnoff, providing excellent leverage on Crater's physical properties. Structurally, Crater has a half-light radius of $\sim$20 pc and shows no evidence for tidal distortions. Crater is well-described by a simple stellar population with an age of $\sim$7.5 Gyr, [M/H]$\sim-1.65$, a M$_{\star}\sim10^4$ M$_{\odot}$, M$_{\rm V}\sim -5.3$, located at a distance of d$\sim$ 145 kpc, with modest uncertainties in these properties due to differences in the underlying stellar evolution models. The sparse sampling of stars above the turnoff and sub-giant branch are likely to be 1.0-1.4 M$_{\odot}$ binary star systems (blue stragglers) and their evolved descendants, as opposed to intermediate age main sequence stars. Confusion of these populations highlights a substantial challenge in accurately characterizing sparsely populated stellar systems. Our analysis shows that Crater is not a dwarf galaxy, but instead is an unusually young cluster given its location in the Milky Way's very outer stellar halo. Crater is similar to SMC cluster Lindsay 38, and its position and velocity are in good agreement with observations and models of the Magellanic stream debris, suggesting it may have accreted from the Magellanic Clouds. However, its age and metallicity are also in agreement with the age-metallicity relationships of lower mass dwarf galaxies such as Leo I or Carina. Despite uncertainty over its progenitor system, Crater appears to have been incorporated into the Galaxy more recently than $z\sim1$ (8 Gyr ago), providing an important new constraint on the accretion history of the Milky Way. [abridged]

We explore the effect of galactic environment on properties of molecular clouds. Using clouds formed in a large-scale galactic disc simulation, we measure the observable properties from synthetic column density maps. We confirm that a significant fraction of unbound clouds forms naturally in a galactic disc environment and that a mixed population of bound and unbound clouds can match observed scaling relations and distributions for extragalactic molecular clouds. By dividing the clouds into inner and outer disc populations, we compare their distributions of properties and test whether there are statistically significant differences between them. We find that clouds in the outer disc have lower masses, sizes, and velocity dispersions as compared to those in the inner disc for reasonable choices of the inner/outer boundary. We attribute the differences to the strong impact of galactic shear on the disc stability at large galactocentric radii. In particular, our Toomre analysis of the disc shows a narrowing envelope of unstable masses as a function of radius, resulting in the formation of smaller, lower mass fragments in the outer disc. We also show that the star formation rate is affected by the environment of the parent cloud, and is particularly influenced by the underlying surface density profile of the gas throughout the disc. Our work highlights the strengths of using galaxy-scale simulations to understand the formation and evolution of cloud properties - and the star formation within them - in the context of their environment.

We present results of a Hubble Space Telescope far-ultraviolet (FUV) survey searching for white dwarf (WD) companions to blue straggler stars (BSSs) in open cluster NGC 188. The majority of NGC 188 BSSs (15 of 21) are single-lined binaries with properties suggestive of mass-transfer formation via Roche lobe overflow, specifically through an asymptotic giant branch star transferring mass to a main sequence secondary, yielding a BSS binary with a WD companion. In NGC 188, a BSS formed by this mechanism within the past 400 Myr will have a WD companion hot and luminous enough to be directly detected as a FUV photometric excess with HST. Comparing expected BSS FUV emission to observed photometry reveals four BSSs with WD companions above 12,000 K (younger than 250 Myr) and three WD companions with temperatures between 11,000-12,000 K. These BSS+WD binaries all formed through recent mass transfer. The location of the young BSSs in an optical color-magnitude diagram (CMD) indicates that distance from the zero-age main sequence does not necessarily correlate with BSS age. There is no clear CMD separation between mass transfer-formed BSSs and those likely formed through other mechanisms, such as collisions. The seven detected WD companions place a lower limit on the mass-transfer formation frequency of 33%. We consider other possible formation mechanisms by comparing properties of the BSS population to theoretical predictions. We conclude that 14 BSS binaries likely formed from mass transfer, resulting in an inferred mass-transfer formation frequency of approximately 67%.

We investigate how the Milky Way tidal field can affect the spatial mixing of multiple stellar populations in the globular cluster NGC 6362. We use $N$-body simulations of multiple population clusters on the orbit of this cluster around the Milky Way. Models of the formation of multiple populations in globular clusters predict that the second population should initially be more centrally concentrated than the first. However, NGC 6362 is comprised of two chemically distinct stellar populations having the same radial distribution. We show that the high mass loss rate experienced on this cluster's orbit significantly accelerates the spatial mixing of the two populations expected from two body relaxation. We also find that for a range of initial second population concentrations, cluster masses, tidal filling factors and fraction of first population stars, a cluster with two populations should be mixed when it has lost 70-80 per cent of its initial mass. These results fully account for the complete spatial mixing of NGC 6362, since, based on its shallow present day mass function, independent studies estimate that the cluster has lost 85 per cent of its initial mass.

We analyze the physical properties of stellar clusters that are detected in massive star-forming regions in the MYStIX project--a comparative, multiwavelength study of young stellar clusters within 3.6 kpc that contain at least one O-type star. Tabulated properties of subclusters in these regions include physical sizes and shapes, intrinsic numbers of stars, absorptions by the molecular clouds, and median subcluster ages. Physical signs of dynamical evolution are present in the relations of these properties, including statistically significant correlations between subcluster size, central density, and age, which are likely the result of cluster expansion after gas removal. We argue that many of the subclusters identified in Paper I are gravitationally bound because their radii are significantly less than what would be expected from freely expanding clumps of stars with a typical initial stellar velocity dispersion of ~3 km/s for star-forming regions. We explore a model for cluster formation in which structurally simpler clusters are built up hierarchically through the mergers of subclusters--subcluster mergers are indicated by an inverse relation between the numbers of stars in a subcluster and their central densities (also seen as a density vs. radius relation that is less steep than would be expected from pure expansion). We discuss implications of these effects for the dynamical relaxation of young stellar clusters.

In this paper, we discuss the origin of the observed correlation between cluster concentration c and present-day mass function (PDMF) slope {\alpha} reported by De Marchi, Paresce & Pulone. This relation can either be reproduced from universal initial conditions combined with some dynamical mechanism(s) that alter(s) the cluster structure and mass function over time, or it must arise early on in the cluster lifetime, such as during the gas-embedded phase of cluster formation. Using a combination of Monte Carlo and N-body models for globular cluster evolution performed with the MOCCA and NBODY6 codes, respectively, we explore a number of dynamical mechanisms that could affect the observed relation. For the range of initial conditions considered here, our results are consistent with an universal initial binary fraction \approx 10% (which does not, however, preclude 100%) and an universal initial stellar mass function resembling the standard Kroupa distribution. Most of the dispersion observed in the c-{\alpha} relation can be attributed to two-body relaxation and Galactic tides. However, dynamical processes alone could not have reproduced the dispersion in concentration, and we require at least some correlation between the initial concentration and the total cluster mass. We argue that the origin of this trend could be connected to the gas-embedded phase of cluster evolution.

We consider the effects that accretion from the interstellar medium onto the particles of an N-body system has on the rate of two-body relaxation. To this end, we derive an accretion-modified relaxation time by adapting Spitzer's two-component model to include the damping effects of accretion. We consider several different mass-dependencies and efficiency factors for the accretion rate, as well as different mass ratios for the two components of the model. The net effect of accretion is to accelerate mass segregation by increasing the average mass $\bar{m}$, since the relaxation time is inversely proportional to $\bar{m}$. Under the assumption that the accretion rate increases with the accretor mass, there are two additional effects that accelerate mass segregation. First, accretion acts to increase the range of any initial mass spectrum, quickly driving the heaviest members to even higher masses. Second, accretion acts to reduce the velocities of the accretors due to conservation of momentum, and it is the heaviest members that are affected the most. Using our two-component model, we quantify these effects as a function of the accretion rate, the total cluster mass, and the component masses. We conclude by discussing the implications of our results for the dynamical evolution of primordial globular clusters, primarily in the context of black holes formed from the most massive stellar progenitors.

(abridged) In this paper, we constrain the properties of primordial binary populations in Galactic globular clusters using the MOCCA Monte Carlo code for cluster evolution. Our results are compared to the observations of Milone et al. (2012) using the photometric binary populations as proxies for the true underlying distributions, in order to test the hypothesis that the data are consistent with an universal initial binary fraction near unity and the binary orbital parameter distributions of Kroupa (1995). With the exception of a few possible outliers, we find that the data are to first-order consistent with the universality hypothesis. Specifically, the present-day binary fractions inside the half-mass radius r$_{\rm h}$ can be reproduced assuming either high initial binary fractions near unity with a dominant soft binary component as in the Kroupa distribution combined with high initial densities (10$^4$-10$^6$ M$_{\odot}$ pc$^{-3}$), or low initial binary fractions ($\sim$ 5-10%) with a dominant hard binary component combined with moderate initial densities near their present-day values (10$^2$-10$^3$ M$_{\odot}$ pc$^{-3}$). This apparent degeneracy can be broken using the binary fractions outside r$_{\rm h}$- only high initial binary fractions with a significant soft component combined with high initial densities can contribute to reproducing the observed anti-correlation between the binary fractions outside r$_{\rm h}$ and the total cluster mass. We further illustrate using the simulated present-day binary orbital parameter distributions and the technique introduced in Leigh et al. (2012) that the relative fractions of hard and soft binaries can be used to further constrain the initial cluster density and mass-density relation. Our results favour an initial mass-density relation of the form r$_{\rm h} \propto$ M$_{\rm clus}^{\alpha}$ with $\alpha <$ 1/3.

We present a multi-passband photometric study of the Blue Straggler Star (BSS) population in the Galactic globular cluster (GC) NGC\,1261, using available space- and ground-based survey data.~The inner BSS population is found to have two distinct sequences in the color-magnitude diagram, similar to double BSS sequences detected in other GCs. These well defined sequences are presumably linked to single short-lived events such as core collapse, which are expected to boost the formation of BSSs.~In agreement with this, we find a BSS sequence in NGC\,1261 which can be well reproduced individually by a theoretical model prediction of a 2 Gyr old population of stellar collision products, which are expected to form in the denser inner regions during short-lived core contraction phases.~Additionally, we report the occurrence of a group of BSSs with unusually blue colours in the CMD, which are consistent with a corresponding model of a 200 Myr old population of stellar collision products.~The properties of the NGC\,1261 BSS populations, including their spatial distributions, suggest an advanced dynamical evolutionary state of the cluster, but the core of this GC does not show the classical signatures of core-collapse.~We argue these apparent contradictions provide evidence for a post-core-collapse bounce state seen in dynamical simulations of old GCs.

We perform N-body simulations of a cluster that forms in a dwarf galaxy and is then accreted by the Milky Way to investigate how a cluster's structure is affected by a galaxy merger. We find that the cluster's half mass radius will respond quickly to this change in potential. When the cluster is placed on an orbit in the Milky Way with a stronger tidal field the cluster experiences a sharp decrease in size in response to increased tidal forces. Conversely, when placed on an orbit with a weaker tidal field the cluster expands since tidal forces decrease and stars moving outwards due to internal effects remain bound at further distances than before. In all cases, we find that the cluster's half mass radius will eventually be indistinguishable from a cluster that has always lived in the Milky Way on that orbit. These adjustments occur within 1-2 half mass relaxation times of the cluster in the dwarf galaxy. We also find this effect to be qualitatively independent of the time that the cluster is taken from the dwarf galaxy. In contrast to the half mass radius, we show the core radius of the cluster is not affected by the potential the cluster lives in. Our work suggests that structural properties of accreted clusters are not distinct from clusters born in the Milky Way. Other cluster properties, such as metallicity and horizontal branch morphology, may be the only way to identify accreted star clusters in the Milky Way.

The structure of molecular clouds can be characterized with the probability distribution function (PDF) of the mass surface density. In particular, the properties of the distribution can reveal the nature of the turbulence and star formation present inside the molecular cloud. In this paper, we explore how these structural characteristics evolve with time and also how they relate to various cloud properties as measured from a sample of synthetic column density maps of molecular clouds. We find that, as a cloud evolves, the peak of its column density PDF will shift to surface densities below the observational threshold for detection, resulting in an underlying lognormal distribution which has been effectively lost at late times. Our results explain why certain observations of actively star-forming, dynamically older clouds, such as the Orion molecular cloud, do not appear to have any evidence of a lognormal distribution in their column density PDFs. We also study the evolution of the slope and deviation point of the power-law tails for our sample of simulated clouds and show that both properties trend towards constant values, thus linking the column density structure of the molecular cloud to the surface density threshold for star formation.

We have performed N-body simulations of tidally filling star clusters with a range of orbits in a Milky Way-like potential to study the effects of orbital inclination and eccentricity on their structure and evolution. At small galactocentric distances Rgc, a non-zero inclination results in increased mass loss rates. Tidal heating and disk shocking, the latter sometimes consisting of two shocking events as the cluster moves towards and away from the disk, help remove stars from the cluster. Clusters with inclined orbits at large Rgc have decreased mass loss rates than the non-inclined case, since the strength the disk potential decreases with Rgc. Clusters with inclined and eccentric orbits experience increased tidal heating due to a constantly changing potential, weaker disk shocks since passages occur at higher Rgc, and an additional tidal shock at perigalacticon. The effects of orbital inclination decrease with orbital eccentricity, as a highly eccentric cluster spends the majority of its lifetime at a large Rgc. The limiting radii of clusters with inclined orbits are best represented by the rt of the cluster when at its maximum height above the disk, where the cluster spends the majority of its lifetime and the rate of change in rt is a minimum. Conversely, the effective radius is independent of inclination in all cases.

We use N-body simulations to explore the influence of orbital eccentricity on the dynamical evolution of star clusters. Specifically we compare the mass loss rate, velocity dispersion, relaxation time, and the mass function of star clusters on circular and eccentric orbits. For a given perigalactic distance, increasing orbital eccentricity slows the dynamical evolution of a cluster due to a weaker mean tidal field. However, we find that perigalactic passes and tidal heating due to an eccentric orbit can partially compensate for the decreased mean tidal field by energizing stars to higher velocities and stripping additional stars from the cluster, accelerating the relaxation process. We find that the corresponding circular orbit which best describes the evolution of a cluster on an eccentric orbit is much less than its semi-major axis or time averaged galactocentric distance. Since clusters spend the majority of their lifetimes near apogalacticon, the properties of clusters which appear very dynamically evolved for a given galactocentric distance can be explained by an eccentric orbit. Additionally we find that the evolution of the slope of the mass function within the core radius is roughly orbit-independent, so it could place additional constraints on the initial mass and initial size of globular clusters with solved orbits. We use our results to demonstrate how the orbit of Milky Way globular clusters can be constrained given standard observable parameters like galactocentric distance and the slope of the mass function. We then place constraints on the unsolved orbits of NGC 1261,NGC 6352, NGC 6496, and NGC 6304 based on their positions and mass functions.

This chapter describes the current state of models of individual blue stragglers. Stellar collisions, binary mergers (or coalescence), and partial or ongoing mass transfer have all been studied in some detail. The products of stellar collisions retain memory of their parent stars and are not fully mixed. Very high initial rotation rates must be reduced by an unknown process to allow the stars to collapse to the main sequence. The more massive collision products have shorter lifetimes than normal stars of the same mass, while products between low mass stars are long-lived and look very much like normal stars of their mass. Mass transfer can result in a merger, or can produce another binary system with a blue straggler and the remnant of the original primary. The products of binary mass transfer cover a larger portion of the colour-magnitude diagram than collision products for two reasons: there are more possible configurations which produce blue stragglers, and there are differing contributions to the blended light of the system. The effects of rotation may be substantial in both collision and merger products, and could result in significant mixing unless angular momentum is lost shortly after the formation event. Surface abundances may provide ways to distinguish between the formation mechanisms, but care must be taking to model the various mixing mechanisms properly before drawing strong conclusions. Avenues for future work are outlined.

The clusters of young stars in massive star-forming regions show a wide range of sizes, morphologies, and numbers of stars. Their highly subclustered structures are revealed by the MYStIX project's sample of 31,754 young stars in nearby sites of star formation (regions at distances <3.6 kpc that contain at least one O-type star.) In 17 of the regions surveyed by MYStIX, we identify subclusters of young stars using finite mixture models -- collections of isothermal ellipsoids that model individual subclusters. Maximum likelihood estimation is used to estimate the model parameters, and the Akaike Information Criterion is used to determine the number of subclusters. This procedure often successfully finds famous subclusters, such as the BN/KL complex behind the Orion Nebula Cluster and the KW-object complex in M17. A catalog of 142 subclusters is presented, with 1 to 20 subclusters per region. The subcluster core radius distribution for this sample is peaked at 0.17 pc with a standard deviation of 0.43 dex, and subcluster core radius is negatively correlated with gas/dust absorption of the stars -- a possible age effect. Based on the morphological arrangements of subclusters, we identify four classes of spatial structure: long chains of subclusters, clumpy structures, isolated clusters with a core-halo structure, and isolated clusters well fit by a single isothermal ellipsoid.

We present a framework that explains the commonly observed variation in light element abundances in globular clusters. If globular clusters form in the centres of dwarf galaxies, they will be pumped onto larger orbits as star formation progresses. The potential well will only retain the moderate velocity AGB ejecta, the expected source of enrichment, but not supernova ejecta. There is no need to increase the initial cluster mass, a requirement of self-enrichment scenarios, as all the stars within the dwarf can contribute. As the clusters move through the dwarf centre they sweep up a mix of AGB ejecta and in-falling pristine gas to form a second generation of stars. The specific mix will vary in time and is thus able to explain the spread in second generation abundances observed in different clusters. The globular clusters will survive to the present day or be stripped as part of the hierarchical merging process of larger galaxies. We illustrate how this process may operate using a high-resolution simulation of a dwarf galaxy at high redshift.

Several possible formation pathways for blue straggler stars have been developed recently, but no one pathway has yet been observationally confirmed for a specific blue straggler. Here we report the first findings from a Hubble Space Telescope ACS/SBC far-UV photometric program to search for white dwarf companions to blue straggler stars. We find three hot and young white dwarf companions to blue straggler stars in the 7-Gyr open cluster NGC 188, indicating that mass transfer in these systems ended less than 300 Myr ago. These companions are direct and secure observational evidence that these blue straggler stars were formed through mass transfer in binary stars. Their existence in a well-studied cluster environment allows for observational constraints of both the current binary system and the progenitor binary system, mapping the entire mass transfer history.

We explore whether observed molecular clouds could include a substantial population of unbound clouds. Using simulations which include only turbulence and gravity, we are able to match observed relations and naturally reproduce the observed scatter in the cloud size-linewidth coefficient, at fixed surface density. We identify the source of this scatter as a spread in the intrinsic virial parameter. Thus these observational trends do not require that clouds exist in a state of dynamical equilibrium. We demonstrate that cloud virial parameters can be accurately determined observationally with an appropriate size estimator. All our simulated clouds eventually form collapsing cores, regardless of whether the cloud is bound overall. This supports the idea that molecular clouds do not have to be bound to form stars or to have observed properties like those of nearby low-mass clouds.

We present new Hubble Space Telescope imaging of the outer regions of M87 in order to study its globular cluster (GC) population out to large galactocentric distances. We discuss particularly the relationship between GC effective radii $r_h$ and projected galactocentric distance $R_{gc}$. The observations suggest a shallow trend $r_h \propto R_{gc}^{0.14}$ out to $R_{gc} \sim 100$ kpc, in agreement with studies of other giant elliptical galaxies. To theoretically reproduce this relationship we simulate GC populations with various distributions of orbits. For an isotropic distribution of cluster orbits we find a steeper trend of $r_h \propto R_{gc}^{0.4}$. Instead we suggest that (a) if the cluster system has an orbital anisotropy profile, where orbits become preferentially radial with increasing galactocentric distance, AND (b) if clusters become more tidally under-filling with galactocentric distance, the observed relationship can be recovered. We also apply this approach to the red and blue GC populations separately and predict that red clusters are preferentially under-filling at large $R_{gc}$ and have a more isotropic distribution of orbits than blue clusters.

We used high-quality images acquired with the WFC3 on board the HST to probe the blue straggler star (BSS) population of the Galactic globular cluster NGC 362. We have found two distinct sequences of BSS: this is the second case, after M 30, where such a feature has been observed. Indeed the BSS location, their extension in magnitude and color and their radial distribution within the cluster nicely resemble those observed in M 30, thus suggesting that the same interpretative scenario can be applied: the red BSS sub-population is generated by mass transfer binaries, the blue one by collisions. The discovery of four new W UMa stars, three of which lying along the red-BSS sequence, further supports this scenario. We also found that the inner portion of the density profile deviates from a King model and is well reproduced by either a mild power-law (\alpha -0.2) or a double King profile. This feature supports the hypothesis that the cluster is currently undergoing the core collapse phase. Moreover, the BSS radial distribution shows a central peak and monotonically decreases outward without any evidence of an external rising branch. This evidence is a further indication of the advanced dynamical age of NGC 362: in fact, together with M 30, NGC 362 belongs to the family of dynamically old clusters (Family III) in the "dynamical clock" classification proposed by Ferraro et al. (2012). The observational evidence presented here strengthens the possible connection between the existence of a double BSS sequence and a quite advanced dynamical status of the parent cluster.

MYStIX (Massive Young Star-Forming Complex Study in Infrared and Xray) seeks to characterize 20 OB-dominated young clusters and their environs at distances d < 4 kpc using imaging detectors on the Chandra X-ray Observatory, Spitzer Space Telescope, and the United Kingdom InfraRed Telescope. The observational goals are to construct catalogs of star-forming complex stellar members with well-defined criteria, and maps of nebular gas (particularly of hot X-ray emitting plasma) and dust. A catalog of MYStIX Probable Complex Members (MPCMs) with several hundred OB stars and > 30, 000 low mass premain sequence is assembled. This sample and related data products will be used to seek new empirical constraints on theoretical models of cluster formation and dynamics, mass segregation, OB star formation, star formation triggering on the periphery of HII regions, the survivability of protoplanetary disks in HII regions. This paper give an introduction and overview of the project, covering the data analysis methodology and application to two star forming regions, NGC 2264 and the Trifid Nebula.

Blue straggler stars (BSS) are abundantly observed in all Galactic globular clusters (GGC) where data exist. However, observations alone cannot reveal the relative importance of various formation channels or the typical formation times for this well studied population of anomalous stars. Using a state-of-the-art H\'enon-type Monte Carlo code that includes all relevant physical processes, we create 128 models with properties typical of the observed GGCs. These models include realistic numbers of single and binary stars, use observationally motivated initial conditions, and span large ranges in central density, concentration, binary fraction, and mass. Their properties can be directly compared with those of observed GGCs. We can easily identify the BSSs in our models and determine their formation channels and birth times. We find that for central densities above ~10^3 Msun/pc^3 the dominant formation channel is stellar collisions while for lower density clusters, mass transfer in binaries provides a significant contribution (up to ~ 60% in our models). The majority of these collisions are binary-mediated, occurring during 3-body and 4-body interactions. As a result a strong correlation between the specific frequency of BSSs and the binary fraction in a cluster can be seen in our models. We find that the number of BSSs in the core shows only a weak correlation with the collision rate estimator \Gamma traditionally used by observers, in agreement with the latest Hubble Space Telescope (ACS) data. Using an idealized "full mixing" prescription for collision products, our models indicate that the BSSs observed today may have formed several Gyrs ago. However, denser clusters tend to have younger (~1 Gyr) BSSs.

We created artificial color-magnitude diagrams of Monte Carlo dynamical models of globular clusters, and then used observational methods to determine the number of blue stragglers in those clusters. We compared these blue stragglers to various cluster properties, mimicking work that has been done for blue stragglers in Milky Way globular clusters to determine the dominant formation mechanism(s) of this unusual stellar population. We find that a mass-based prescription for selecting blue stragglers will choose approximately twice as many blue stragglers than a selection criterion that was developed for observations of real clusters. However, the two numbers of blue stragglers are well-correlated, so either selection criterion can be used to characterize the blue straggler population of a cluster. We confirm previous results that the simplified prescription for the evolution of a collision or merger product in the BSE code overestimates their lifetimes. We show that our model blue stragglers follow similar trends with cluster properties (core mass, binary fraction, total mass, collision rate) as the true Milky Way blue stragglers, as long as we restrict ourselves to model clusters with an initial binary fraction higher than 5%. We also show that, in contrast to earlier work, the number of blue stragglers in the cluster core does have a weak dependence on the collisional parameter Gamma in both our models and in Milky Way globular clusters.

A long-standing issue in the theory of low mass stars is the discrepancy between predicted and observed radii and effective temperatures. In spite of the increasing availability of very precise radius determinations from eclipsing binaries and interferometric measurements of radii of single stars, there is no unanimous consensus on the extent (or even the existence) of the discrepancy and on its connection with other stellar properties (e.g. metallicity, magnetic activity). We investigate the radius discrepancy phenomenon using the best data currently available (accuracy about 5%). We have constructed a grid of stellar models covering the entire range of low mass stars (0.1-1.25 M_sun) and various choices of the metallicity and of the mixing length parameter \alpha. We used an improved version of the Yale Rotational stellar Evolution Code (YREC), implementing surface boundary conditions based on the most up-to-date PHOENIX atmosphere models. Our models are in good agreement with others in the literature and improve and extend the low mass end of the Yale-Yonsei isochrones. Our calculations include rotation-related quantities, such as moments of inertia and convective turnover time scales, useful in studies of magnetic activity and rotational evolution of solar-like stars. Consistently with previous works, we find that both binaries and single stars have radii inflated by about 3% with respect to the theoretical models; among binaries, the components of short orbital period systems are found to be the most deviant. We conclude that both binaries and single stars are comparably affected by the radius discrepancy phenomenon.

We have performed N-body simulations of star clusters orbiting in a spherically symmetric smooth galactic potential. The model clusters cover a range of initial half-mass radii and orbital eccentricities in order to test the historical assumption that the tidal radius of a cluster is imposed at perigalacticon. The traditional assumption for globular clusters is that since the internal relaxation time is larger than its orbital period, the cluster is tidally stripped at perigalacticon. Instead, our simulations show that a cluster with an eccentric orbit does not need to fully relax in order to expand. After a perigalactic pass, a cluster re-captures previously unbound stars, and the tidal shock at perigalacticon has the effect of energizing inner region stars to larger orbits. Therefore, instead of the limiting radius being imposed at perigalacticon, it more nearly traces the instantaneous tidal radius of the cluster at any point in the orbit. We present a numerical correction factor to theoretical tidal radii calculated at perigalacticon which takes into consideration both the orbital eccentricity and current orbital phase of the cluster.

Globular star clusters that formed at the same cosmic time may have evolved rather differently from a dynamical point of view (because that evolution depends on the internal environment) through a variety of processes that tend progressively to segregate stars more massive than the average towards the cluster centre. Therefore clusters with the same chronological age may have reached quite different stages of their dynamical history (that is, they may have different dynamical ages). Blue straggler stars have masses greater than those at the turn-off point on the main sequence and therefore must be the result of either a collision or a mass-transfer event. Because they are among the most massive and luminous objects in old clusters, they can be used as test particles with which to probe dynamical evolution. Here we report that globular clusters can be grouped into a few distinct families on the basis of the radial distribution of blue stragglers. This grouping corresponds well to an effective ranking of the dynamical stage reached by stellar systems, thereby permitting a direct measure of the cluster dynamical age purely from observed properties.

Metal-rich (red) globular clusters in massive galaxies are, on average, smaller than metal-poor (blue) globular clusters. One of the possible explanations for this phenomenon is that the two populations of clusters have different spatial distributions. We test this idea by comparing clusters observed in unusually deep, high signal-to-noise images of M87 with a simulated globular cluster population in which the red and blue clusters have different spatial distributions, matching the observations. We compare the overall distribution of cluster effective radii as well as the relationship between effective radius and galactocentric distance for both the observed and simulated red and blue subpopulations. We find that the different spatial distributions does not produce a significant size difference between the red and blue subpopulations as a whole, or at a given galactocentric distance. These results suggest that the size difference between red and blue globular clusters is likely due to differences during formation or later evolution

(abridged) We use newly available empirical binary fractions for globular clusters to carry out a direct test of the binary evolution hypothesis, and of collisional channels that involve binary stars. More specifically, using the previously reported correlation between blue straggler numbers and core mass as a benchmark, we test for correlations with the number of binary stars, as well as with the rates of single-single, single-binary, and binary-binary encounters. Surprisingly, we find that the simple correlation with core mass remains by far the strongest predictor of blue straggler population size, even in our joint models. This is despite the fact that the binary fractions themselves strongly anti-correlate with core mass, just as expected in the binary evolution model. At first sight, these results do not fit neatly with either binary evolution or collisional models in their simplest forms. Arguably the simplest and most intriguing possibility to explain this unexpected result is that observational errors on the core binary fractions are larger than the true intrinsic dispersion associated with their dependence on core mass. In the context of the binary evolution model, this would explain why the combination of binary fraction and core mass is a poorer predictor of blue straggler numbers than core mass alone. It would also imply that core mass is a remarkably clean predictor of core binary fractions. This would be of considerable importance for the dynamical evolution of globular clusters, and provides an important benchmark for models attempting to understand their present-day properties.

In this paper, we investigate the extent to which observations of molecular clouds can correctly identify and measure star-forming clumps. We produced a synthetic column density map and a synthetic spectral-line data cube from the simulated collapse of a 5000 M$_{\odot}$ molecular cloud. By correlating the clumps found in the simulation to those found in the synthetic observations, clump masses derived from spectral-line data cubes were found to be quite close to the true physical properties of the clumps. We also find that the `observed' clump mass function derived from the column density map is shifted by a factor of ~ 3 higher than the true clump mass function, due to projection of low-density material along the line of sight. Alves et al. (2007) first proposed that a shift of a clump mass function to higher masses by a factor of 3 can be attributed to a star formation efficiency of 30 %. Our results indicate that this finding may instead be due to an overestimate of clump masses determined from column density observations.

We present a new technique to quantify cluster-to-cluster variations in the observed present-day stellar mass functions of a large sample of star clusters. Our method quantifies these differences as a function of both the stellar mass and the total cluster mass, and offers the advantage that it is insensitive to the precise functional form of the mass function. We applied our technique to data taken from the ACS Survey for Globular Clusters, from which we obtained completeness-corrected stellar mass functions in the mass range 0.25-0.75 M$_{\odot}$ for a sample of 27 clusters. The results of our observational analysis were then compared to Monte Carlo simulations for globular cluster evolution spanning a range of initial mass functions, total numbers of stars, concentrations, and virial radii. We show that the present-day mass functions of the clusters in our sample can be reproduced by assuming an universal initial mass function for all clusters, and that the cluster-to-cluster differences are consistent with what is expected from two-body relaxation. A more complete exploration of the initial cluster conditions will be needed in future studies to better constrain the precise functional form of the initial mass function. This study is a first step toward using our technique to constrain the dynamical histories of a large sample of old Galactic star clusters and, by extension, star formation in the early Universe.

Globular clusters have linear sizes (tidal radii) which theory tells us are determined by their masses and by the gravitational potential of their host galaxy. To explore the relationship between observed and expected radii, we utilize the globular cluster population of the Virgo giant M87. Unusually deep, high signal-to-noise images of M87 are used to measure the effective and limiting radii of approximately 2000 globular clusters. To compare with these observations, we simulate a globular cluster population that has the same characteristics as the observed M87 cluster population. Placing these simulated clusters in the well-studied tidal field of M87, the orbit of each cluster is solved and the theoretical tidal radius of each cluster is determined. We compare the predicted relationship between cluster size and projected galactocentric distance to observations. We find that for an isotropic distribution of cluster velocities, theoretical tidal radii are approximately equal to observed limiting radii for Rgc < 10 kpc. However, the isotropic simulation predicts a steep increase in cluster size at larger radii, which is not observed in large galaxies beyond the Milky Way. To minimize the discrepancy between theory and observations, we explore the effects of orbital anisotropy on cluster sizes, and suggest a possible orbital anisotropy profile for M87 which yields a better match between theory and observations. Finally, we suggest future studies which will establish a stronger link between theoretical tidal radii and observed radii.

Observations indicate blue globular clusters have half-light radii systematically larger than those of red globular clusters. In this paper, we test whether the different metallicity-dependent stellar evolution timescales and mass-loss rates within the clusters can impact their early dynamical evolution. By means of N-body simulations including stellar evolution recipes we simulate the early evolution of small centrally concentrated clusters with and without primordial mass segregation. Our simulations include accurate metallicity-dependent mass loss from massive stars. We find blue clusters to be larger than red clusters regardless of whether the clusters have been primordially mass segregated. In addition, the size difference is found to be larger and consistent with observations for initial models with a low central concentration. These results indicate that the systematic size difference found between red and blue clusters can, at least in part, be attributed to the dynamical effects of differing stellar evolution histories, driven by metallicity.

We present an analytic model for blue straggler formation in globular clusters. We assume that blue stragglers are formed only through stellar collisions and binary star evolution, and compare our predictions to observed blue straggler numbers taken from the catalogue of Leigh, Sills & Knigge (2011). We can summarize our key results as follows: (1) Binary star evolution consistently dominates blue straggler production in all our best-fitting models. (2) In order to account for the observed sub-linear dependence of blue straggler numbers on the core masses (Knigge, Leigh & Sills 2009), the core binary fraction must be inversely proportional to the total cluster luminosity and should always exceed at least a few percent. (3) In at least some clusters, blue straggler formation must be enhanced by dynamical encounters (either via direct collisions or by stimulating mass-transfer to occur by altering the distribution of binary orbital parameters) relative to what is expected by assuming a simple population of binaries evolving in isolation. (4) The agreement between the predictions of our model and the observations can be improved by including blue stragglers that form outside the core but later migrate in due to dynamical friction. (5) Longer blue straggler lifetimes are preferred in models that include blue stragglers formed outside the core since this increases the fraction that will have sufficient time to migrate in via dynamical friction.

We present a homogeneous catalogue for blue straggler, red giant branch, horizontal branch and main-sequence turn-off stars in a sample of 35 clusters taken from the ACS Survey for Globular Clusters. As a result of the superior photometry and relatively large field of view offered by the ACS data, this new catalogue is a significant improvement upon the one presented in Leigh, Sills & knigge (2007). Using our catalogue, we study and compare the radial distributions of the different stellar populations. We have confirmed our previous result (Knigge, Leigh & Sills 2009) that there is a clear, but sub-linear, correlation between the number of blue stragglers found in the cluster core and the total stellar mass contained within it. By considering a larger spatial extent than just the core, our results suggest that mass segregation is not the dominant effect contributing to the observed sub-linearity. We also investigate the radial distributions of the different stellar populations in our sample of clusters. Our results are consistent with a linear relationship between the number of stars in these populations and the total mass enclosed within the same radius. Therefore, we conclude that the cluster dynamics does not significantly affect the relative distributions of these populations in our sample.

Close interactions and mass transfer in binary stars can lead to the formation of many different exotic stellar populations, but detailed modeling of mass transfer is a computationally challenging problem. Here, we present an alternate Smoothed Particle Hydrodynamics approach to the modeling of mass transfer in binary systems that allows a better resolution of the flow of matter between main-sequence stars. Our approach consists of modeling only the outermost layers of the stars using appropriate boundary conditions and ghost particles. We arbitrarily set the radius of the boundary and find that our boundary treatment behaves physically and conserves energy well. In particular, when used with our binary relaxation procedure, our treatment of boundary conditions is also shown to evolve circular binaries properly for many orbits. The results of our first simulation of mass transfer are also discussed and used to assess the strengths and limitations of our method. We conclude that it is well suited for the modeling of interacting binary stars. The method presented here represents a convenient alternative to previous hydrodynamical techniques aimed at modeling mass transfer in binary systems since it can be used to model both the donor and the accretor while maintaining the density profiles taken from realistic stellar models.

Despite numerous efforts to better understand binary star evolution, some aspects of it remain poorly constrained. In particular, the evolution of eccentric binaries has remained elusive mainly because the Roche lobe formalism derived for circular binaries does not apply. Here, we report the results of our Smoothed Particle Hydrodynamics simulations of mass transfer in eccentric binaries using an alternate method in which we model only the outermost layers of the stars with appropriate boundary conditions. Using this technique, along with properly relaxed model stars, we characterize the mass transfer episodes of binaries with various orbital parameters. In particular, we show that these episodes can be described by Gaussians with a FWMH of ~0.12 P_orb and that the peak rates occur after periastron, at an orbital phase of ~0.58, independently of the eccentricity and mass of the stars. The accreted material is observed to form a rather sparse envelope around either or both stars. Although the fate of this envelope is not modeled in our simulations, we show that a constant fraction (~5%) of the material transferred is ejected from the systems. We discuss this result in terms of the non-conservative mass transfer scenario. We suggest our results could be incorporated in analytical and binary population synthesis studies to help better understand the evolution of eccentric binaries and the formation of exotic stellar populations.

We present a technique to identify the most probable dynamical formation scenario for an observed binary or triple system containing one or more merger products or, alternatively, to rule out the possibility of a dynamical origin. Our method relies on an analytic prescription for energy conservation during stellar encounters. With this, observations of the multiple star system containing the merger product(s) can be used to work backwards in order to constrain the initial orbital energies of any single, binary or triple systems that went into the encounter. The initial semi-major axes of the orbits provide an estimate for the collisional cross section and therefore the time-scale for the encounter to occur in its host cluster. We have applied our analytic prescription to observed binary and triple systems containing blue stragglers, in particular the triple system S1082 in M67 and the period distribution of the blue straggler binaries in NGC 188. We have shown that both S1082 and most of the blue straggler binaries in NGC 188 could have a dynamical origin, and that encounters involving triples could be a significant contributor to BS populations in old open clusters. In general, our results suggest that encounters involving triples could make up a significant fraction of those dynamical interactions that result in stellar mergers, in particular encounters that produce multiple star systems containing one or more blue stragglers.

We present a technique to identify the most probable dynamical formation scenario for an observed binary or triple system containing one or more merger products or, alternatively, to rule out the possibility of a dynamical origin. Our method relies on an analytic prescription for energy conservation during stellar encounters. With this, observations of the multiple star system containing the merger product(s) can be used to work backwards in order to constrain the initial orbital energies of any single, binary or triple systems that went into the encounter. The initial semi-major axes of the orbits provide an estimate for the collisional cross section and therefore the time-scale for the encounter to occur in its host cluster. We have applied our analytic prescription to observed binary and triple systems containing blue stragglers, in particular the triple system S1082 in M67 and the period distribution of the blue straggler binaries in NGC 188. We have shown that both S1082 and most of the blue straggler binaries in NGC 188 could have a dynamical origin, and that encounters involving triples are a significant contributor to blue straggler populations in old open clusters. In general, our results suggest that encounters involving triples could make up a significant fraction of those dynamical interactions that result in stellar mergers, in particular encounters that produce multiple star systems containing one or more blue stragglers.

Blue stragglers are thought to be formed from the merger or coalescence of two stars, but the details of their formation in clusters has been difficult to disentangle. We discuss the two main formation mechanisms for blue stragglers (stellar collisions or mass transfer in a binary system). We then look at the additional complications caused by the stars living in the dynamically active environment of a star cluster. We review the recent observational and theoretical work which addresses the question "which mechanism dominates?" and conclude that the most likely answer is that both mechanisms are at work, although with different importance in different environments and at different times in the cluster lifetime. We finish with a short discussion of some avenues for future work

We have used high resolution spectra obtained with the spectrograph FLAMES at the ESO Very Large Telescope to determine the kinematical properties and the abundance patterns of 20 blue straggler stars (BSSs) in the globular cluster M4. We found that ~ 40% of the measured BSSs are fast rotators (with rotational velocities > 50 km/s). This is the largest frequency of rapidly rotating BSSs ever detected in a globular cluster. In addition, at odds with what has been found in 47 Tucanae, no evidence of carbon and/or oxygen depletion has been revealed in the sample of 11 BSSs for which we were able to measure the abundances. This could be due either to low statistics, or to a different BSS formation process acting in M4.

We use simulated images of star-forming regions to explore the effects of various image acquisition techniques on the derived clump mass function. In particular, we focus on the effects of finite image angular resolution, the presence of noise, and spatial filtering. We find that, even when the image has been so heavily degraded with added noise and lowered angular resolution that the clumps it contains clearly no longer correspond to pre-stellar cores, still the clump mass function is typically consistent with the stellar initial mass function within their mutual uncertainties. We explain this result by suggesting that noise, source blending, and spatial filtering all randomly perturb the clump masses, biasing the mass function toward a lognormal form whose high-mass end mimics a Salpeter power law. We argue that this is a consequence of the central limit theorem and that it strongly limits our ability to accurately measure the true mass function of the clumps. We support this conclusion by showing that the characteristic mass scale of the clump mass function, represented by the ``break mass'', scales as a simple function of the angular resolution of the image from which the clump mass function is derived. This strongly constrains our ability to use the clump mass function to derive a star formation efficiency. We discuss the potential and limitations of the current and next generation of instruments for measuring the clump mass function.

It has become clear in recent years that globular clusters are not simple stellar populations, but may host chemically distinct sub-populations, typically with an enhanced helium abundance. These helium-rich populations can make up a substantial fraction of all cluster stars. One of the proposed formation channels for blue straggler stars is the physical collision and merger of two stars. In the context of multiple populations, collisions between stars with different helium abundances should occur and contribute to the observed blue straggler population. This will affect the predicted blue straggler colour and luminosity function. We quantify this effect by calculating models of mergers resulting from collisions between stars with different helium abundances and using these models to model a merger population. We then compare these results to four observed clusters, NGC 1851, NGC 2808, NGC 5634 and NGC 6093. As in previous studies our models deviate from the observations, particularly in the colour distributions. However, our results are consistent with observations of multiple populations in these clusters. In NGC 2808, our best fitting models include normal and helium enhanced populations, in agreement with helium enhancement inferred in this cluster. The other three clusters show better agreement with models that do not include helium enhancement. We discuss future prospects to improve the modelling of blue straggler populations and the role that the models we present here can play in such a study.

Globular clusters were thought to be simple stellar populations, but recent photometric and spectroscopic evidence suggests that the clusters' early formation history was more complicated. In particular, clusters show star-to-star abundance variations, and multiple sequences in their colour-magnitude diagrams. These effects seem to be restricted to globular clusters, and are not found in open clusters or the field. In this paper, we combine the two competing models for these multiple populations and include a consideration of the effects of stellar collisions. Collisions are one of the few phenomena which occur solely in dense stellar environments like (proto-)globular clusters. We find that runaway collisions between massive stars can produce material which has abundances comparable to the observed second generations, but that very little total mass is produced by this channel. We then add the contributions of rapidly-rotating massive stars (under the assumption that massive stars are spun up by collisions and interactions), and the contribution of asymptotic giant branch stars. We find that collisions can help produce the extreme abundances which are seen in some clusters. However, the total amount of material produced in these generations is still too small (by at least a factor of 10) to match the observations. We conclude with a discussion of the additional effects which probably need to be considered to solve this particular problem.

Stars in globular clusters are generally believed to have all formed at the same time, early in the Galaxy's history. 'Blue stragglers' are stars massive enough that they should have evolved into white dwarfs long ago. Two possible mechanisms have been proposed for their formation: mass transfer between binary companions and stellar mergers resulting from direct collisions between two stars. Recently, the binary explanation was claimed to be dominant. Here we report that there are two distinct parallel sequences of blue stragglers in M30. This globular cluster is thought to have undergone 'core collapse', during which both the collision rate and the mass transfer activity in binary systems would have been enhanced. We suggest that the two observed sequences arise from the cluster core collapse, with the bluer population arising from direct stellar collisions and the redder one arising from the evolution of close binaries that are probably still experiencing an active phase of mass transfer.

We investigate the relationship between the mass of a globular cluster core and the sizes of its various stellar populations in a sample of 56 globular clusters. The number of core red giant branch stars is found to scale sub-linearly with core mass at the 3-$\sigma$ confidence level, whereas the relation is linear to within one standard deviation for main-sequence turn-off and sub-giant branch stars. We interpret our results as evidence for a surplus of red giant branch stars in the least massive cluster cores which is not seen for main-sequence turn-off and sub-giant branch stars. We explore various possibilities for the source of this discrepancy, discussing our results primarily in terms of the interplay between the cluster dynamics and stellar evolution.

Blue stragglers in globular clusters are abnormally massive stars that should have evolved off the stellar main sequence long ago. There are two known processes that can create these objects: direct stellar collisions and binary evolution. However, the relative importance of these processes has remained unclear. In particular, the total number of blue stragglers found in a given cluster does not seem to correlate with the predicted collision rate, providing indirect support for the binary-evolution model. Yet the radial distributions of blue stragglers in many clusters are bimodal, with a dominant central peak: this has been interpreted as an indication that collisions do dominate blue straggler production, at least in the high-density cluster cores. Here we report that there is a clear, but sublinear, correlation between the number of blue stragglers found in a cluster core and the total stellar mass contained within it. From this we conclude that most blue stragglers, even those found in cluster cores, come from binary systems. The parent binaries, however, may themselves have been affected by dynamical encounters. This may be the key to reconciling all of the seemingly conflicting results found to date.

We study the post-main sequence evolution of products of collisions between main sequence stars (blue stragglers), with particular interest paid to the horizontal branch and asymptotic giant branch phases. We found that the blue straggler progeny populate the colour-magnitude diagram slightly blueward of the red giant branch and between 0.2 and 1 magnitudes brighter than the horizontal branch. We also found that the lifetimes of collision products on the horizontal branch is consistent with the numbers of so-called "evolved blue straggler stars" (E-BSS) identified by various authors in a number of globular clusters, and is almost independent of mass or initial composition profile. The observed ratio of the number of E-BSS to blue stragglers points to a main sequence lifetime for blue stragglers of approximately 1-2 Gyr on average.

The formation of blue stragglers is still not completely understood, particularly the relationship between formation environment and mechanism. We use a large, homogeneous sample of blue stragglers in the cores of 57 globular clusters to investigate the relationships between blue straggler populations and their environments. We use a consistent definition of "blue straggler" based on position in the color-magnitude diagram, and normalize the population relative to the number of red giant branch stars in the core. We find that the previously determined anti-correlation between blue straggler frequency and total cluster mass is present in the purely core population. We find some weak correlations with central velocity dispersion and with half-mass relaxation time. The blue straggler frequency does not show any trend with any other cluster parameter. Even though collisions may be expected to be a dominant blue straggler formation process in globular cluster cores, we find no correlation between the frequency of blue stragglers and the collision rate in the core. We also investigated the blue straggler luminosity function shape, and found no relationship between any cluster parameter and the distribution of blue stragglers in the color-magnitude diagram. Our results are inconsistent with some recent models of blue straggler formation that include collisional formation mechanisms, and may suggest that almost all observed blue stragglers are formed in binary systems.

An Eulerian TVD code and a Lagrangian SPH code are used to simulate the off-axis collision of equal-mass main sequence stars in order to address the question of whether stellar mergers can produce a remnant star where the interior has been replenished with hydrogen due to significant mixing. Each parent main sequence star is chosen to be found near the turnoff, with hydrogen depleted in the core, and is modelled with a M=0.8 solar mass realistic stellar model and as a n=3 polytrope. An ideal fluid description with adiabatic index gamma=5/3 is used for all hydrodynamic calculations. We found good agreement between the simulations for the polytropic case, with the remnant showing strong, non-local mixing throughout. In the interior quarter of the mass, ~35% is mixed in from larger radii and on average the remnant is ~50% fully mixed. For the realistic model, we found less mixing, particularly in the interior and in the SPH simulation. In the inner quarter, ~20% of the contained mass in the TVD case, but only ~3% in the SPH one is mixed in from outside. The simulations give consistent results for the overall profile of the merger remnant and the amount of mass loss, but the differences in mixing suggests that the intrinsic difference between grid and particle based schemes remains a possible artifact. We conclude that both the TVD and SPH schemes can be used equally well for problems that are best suited to their strengths and that care should be taken in interpreting results about fluid mixing.

Various galactic globular clusters display abundance anomalies that affect the morphology of their colour-magnitude diagrams. In this paper we consider the possibility of helium enhancement in the anomalous horizontal branch of NGC 2808. We examine the dynamics of a self-enrichment scenario in which an initial generation of stars with a top-heavy initial mass function enriches the interstellar medium with helium via the low-velocity ejecta of its asymptotic giant branch stars. This enriched medium then produces a second generation of stars which are themselves helium-enriched. We use a direct N-body approach to perform five simulations and conclude that such two-generation clusters are both possible and would not differ significantly from their single-generation counterparts on the basis of dynamics. We find, however, that the stellar populations of such clusters would differ from single-generation clusters with a standard initial mass function and in particular would be enhanced in white dwarf stars. We conclude, at least from the standpoint of dynamics, that two-generation globular clusters are feasible.

We derive detailed theoretical models for 1074 nearby stars from the SPOCS (Spectroscopic Properties of Cool Stars) Catalog. The California and Carnegie Planet Search has obtained high-quality echelle spectra of over 1000 nearby stars taken with the Hamilton spectrograph at Lick Observatory, the HIRES spectrograph at Keck, and UCLES at the Anglo Australian Observatory. A uniform analysis of the high-resolution spectra has yielded precise stellar parameters, enabling systematic error analyses and accurate theoretical stellar modeling. We have created a large database of theoretical stellar evolution tracks using the Yale Stellar Evolution Code (YREC) to match the observed parameters of the SPOCS stars. Our very dense grids of evolutionary tracks eliminate the need for interpolation between stellar evolutionary tracks and allow precise determinations of physical stellar parameters (mass, age, radius, size and mass of the convective zone, etc.). Combining our stellar models with the observed stellar atmospheric parameters and uncertainties, we compute the likelihood for each set of stellar model parameters separated by uniform time steps along the stellar evolutionary tracks. The computed likelihoods are used for a Bayesian analysis to derive posterior probability distribution functions for the physical stellar parameters of interest. We provide a catalog of physical parameters for 1074 stars that are based on a uniform set of high quality spectral observations, a uniform spectral reduction procedure, and a uniform set of stellar evolutionary models. We explore this catalog for various possible correlations between stellar and planetary properties, which may help constrain the formation and dynamical histories of other planetary systems.

We use high-resolution spectra obtained with the ESO Very Large Telescope to measure surface abundance patterns of 43 Blue Stragglers stars (BSS) in 47 Tuc. We discovered that a sub-population of BSS shows a significant depletion of Carbon and Oxygen with respect to the dominant population. This evidence would suggest the presence of CNO burning products on the BSS surface coming from a deeply peeled parent star, as expected in the case of mass-transfer process. This is the first detection of a chemical signature clearly pointing to a specific BSS formation process in a globular cluster.

We investigate the suggestion that there are stellar populations in some globular clusters with enhanced helium (Y from 0.28 to 0.40) compared to the primordial value. We assume that a previous generation of massive Asymptotic Giant Branch (AGB) stars have polluted the cluster. Two independent sets of AGB yields are used to follow the evolution of helium and CNO using a Salpeter initial mass function (IMF) and two top-heavy IMFs. In no case are we able to produce the postulated large Y ~ 0.35 without violating the observational constraint that the CNO content is nearly constant.

We use two stellar populations in the globular cluster 47 Tucanae to trace its dynamical history: blue stragglers and low mass main sequence stars. We assumed that the blue stragglers were formed through stellar collisions in all regions of the cluster. We find that in the core of the cluster, models of collisional blue stragglers agree well with the observations as long as blue stragglers are still continuing to form and the mass function in the cluster is extremely biased towards massive stars (x=-8 where a Salpeter mass function has x=+1.35). We show that such an extreme mass function is supported by direct measurements of the luminosity function of main sequence stars in the centre of the cluster. In the middle region of our dataset (25'' to 130'' from the cluster centre), blue straggler formation seems to have stopped about half a Gyr ago. In the outskirts of the cluster, our models are least successful at reproducing the blue straggler data. Taken at face value, they indicate that blue straggler formation has been insignificant over the past billion years, and that a Salpeter mass function applies. However, it is more likely that the dominant formation mechanism in this part of the cluster is not the collisional one, and that our models are not appropriate for this region of the cluster. We conclude that blue stragglers can be used as tracers of dynamics in globular clusters, despite our incomplete understanding of how and where they were formed.

Using a chemical evolution model we investigate the intriguing suggestion that there are populations of stars in some globular clusters (e.g. NGC 2808, omega Centauri) with enhanced levels of helium (Y from about 0.28 to 0.40) compared to the majority of the population that presumably have a primordial helium abundance. We assume that a previous generation of massive low-metallicity Asymptotic Giant Branch (AGB) stars has polluted the cluster gas via a slow stellar wind. We use two independent sets of AGB yields computed from detailed models to follow the evolution of helium, carbon, nitrogen and oxygen in the cluster gas using a Salpeter initial mass function (IMF) and a number of top-heavy IMFs. In no case were we able to fit the observational constraints, Y > 0.30 and C+N+O approximately constant. Depending on the shape of the IMF and the yields, we either obtained Y approximately greater than 0.30 and large increases in C+N+O or Y < 0.30 and C+N+O approximately constant. These results suggest that either AGB stars alone are not responsible for the large helium enrichment or that any dredge-up from this generation of stars was less than predicted by standard models.

By fitting a flexible stellar anisotropy model to the observed surface brightness and line-of-sight velocity dispersion profiles of Draco we derive a sequence of cosmologically plausible two-component (stars + dark matter) models for this galaxy. The models are consistent with all the available observations and can have either cuspy Navarro-Frenk-White or flat-cored dark matter density profiles. The dark matter halos either formed relatively recently (at z~2...7) and are massive (up to ~5x10^9 M_Sun), or formed before the end of the reionization of the universe (z~7...11) and are less massive (down to ~7x10^7 M_Sun). Our results thus support either of the two popular solutions of the "missing satellites" problem of Lambda cold dark matter cosmology - that dwarf spheroidals are either very massive, or very old. We carry out high-resolution simulations of the tidal evolution of our two-component Draco models in the potential of the Milky Way. The results of our simulations suggest that the observable properties of Draco have not been appreciably affected by the Galactic tides after 10 Gyr of evolution. We rule out Draco being a "tidal dwarf" - a tidally disrupted dwarf galaxy. Almost radial Draco orbits (with the pericentric distance <15 kpc) are also ruled out by our analysis. The case of a harmonic dark matter core can be consistent with observations only for a very limited choice of Draco orbits (with the apocentric-to-pericentric distances ratio of <2.5).

Using a large set of high resolution numerical simulations incorporating non-equilibrium molecular hydrogen chemistry and a constant source of external radiation, we study gas collapse in previously photo-ionized mini-galaxies with virial temperatures less than 10^4 K in the early universe (redshifts z=10-20). We confirm that the mechanism of positive feedback of ionizing radiation on star formation in mini-galaxies proposed by Ricotti, Gnedin, & Shull (2002) can be efficient despite a significant flux of metagalactic photo-dissociating radiation. We derive critical fluxes for the Lyman-Werner background radiation sufficient to prevent the collapse of gas in mini-galaxies as a function of the virial mass of the halo and redshift. In our model, the formation of mini-galaxies in defunct HII regions is most efficient at large redshifts (z>15) and/or for large local gas overdensity delta>10. We show that non-equilibrium chemistry plays an important dynamical role not only during the initial evolutionary phase leading to the gas becoming gravitationally unstable inside the mini-halo, but also at the advanced stages of the core collapse, resulting in efficient gas accretion in the core region. We speculate on a possible connection between our objects and metal-poor globular clusters and dwarf spheroidal galaxies.

We investigate the effects of the addition of pre-main sequence evolution to star cluster simulations. We allowed stars to follow pre-main sequence tracks that begin at the deuterium burning birthline and end at the zero age main sequence. We compared our simulations to ones in which the stars began their lives at the zero age main sequence, and also investigated the effects of particular choices for initial binary orbital parameters. We find that the inclusion of the pre-main sequence phase results in a slightly higher core concentration, lower binary fraction, and fewer hard binary systems. In general, the global properties of star clusters remain almost unchanged, but the properties of the binary star population in the cluster can be dramatically modified by the correct treatment of the pre-main sequence stage.

In this paper, we present the results of smoothed particle hydrodynamic (SPH) simulations of collisions between pre-main sequence stars and a variety of other kinds of stars. Simulations over a range of impact parameters and velocities were performed. We find that pre-main sequence stars tend to ``wrap themselves'' around their impactor. We discuss the probable evolutionary state of products of collisions between pre-main sequence stars and pre-main sequence, main sequence, giant branch, and compact stars. The nature of the collision product does not depend strongly on the impact parameter or the velocity of the collision.

We discuss the high lithium abundances in the secondary stars of X-ray binaries. We show that no lithium production in these stars is necessary, and that the abundances can be explained simply due to the tidally locked rotation of the stars, which lead naturally to slower lithium destruction rates. The differences in abundances of CVs' secondaries from those of LMXBs had previously been put forth as evidence that the compact object was related to the lithium abundance, but this scenario also accounts for the lower lithium abundances in the secondary stars in cataclysmic variable systems (CVs) than in low mass X-ray binaries (LMXBs), since these stars have typically lived much longer before becoming tidally locked short period systems. We point out that if this scenario is correct, then the globular cluster X-ray binaries' donor stars should, as a class, show less lithium enhancement relative to other stars of the same spectral type in the clusters than the field X-ray binaries' donor stars show.

We propose a simple model for the formation of dwarf spheroidal galaxies, in which stars are assumed to have formed from isothermal gas in hydrostatic equilibrium inside extended dark matter halos. After expelling the leftover gas, the stellar system undergoes a dynamical relaxation inside the dark matter halo. These models can adequately describe the observed properties of three (Draco, Sculptor, and Carina) out of four Galactic dwarf spheroidal satellites studied in this paper. We suggest that the fourth galaxy (Fornax), which cannot be fitted well with our model, is observed all the way to its tidal radius. Our best fitting models have virial masses of ~10^9 M_Sun, halo formation redshifts consistent with the age of oldest stars in these dwarfs, and shallow inner dark matter density profiles (with slope gamma -0.5...0). The inferred temperature of gas is ~10^4 K. In our model, the "extratidal" stars observed in the vicinity of some dwarf spheroidal galaxies are gravitationally bound to the galaxies and are a part of the extended stellar halos. The inferred virial masses make Galactic dwarf spheroidals massive enough to alleviate the "missing satellites" problem of LCDM cosmologies.

We investigate the structure and evolution of blue stragglers stars which were formed from direct stellar collisions between main sequence stars in globular clusters. In particular, we look at the rotational evolution of the products of off-axis collisions. As found in previous work, such blue stragglers initially have too high an angular momentum to contract down to the main sequence. We consider angular momentum loss through either disc locking or locking to an outflowing wind and show that both methods allow the merged object to shed sufficient angular momentum to contract down to the main sequence.

In this second paper in our series, we continue to test primordial scenarios of globular cluster formation which predict that globular clusters formed in the early universe in the potential of dark matter minihalos. In this paper we use high-resolution N-body simulations to model tidal stripping experienced by primordial dark-matter dominated globular clusters in the static gravitational potential of the host dwarf galaxy. We test both cuspy Navarro-Frenk-White (NFW) and flat-core Burkert models of dark matter halos. Our primordial globular cluster with an NFW dark matter halo survives severe tidal stripping, and after 10 orbits is still dominated by dark matter in its outskirts. Our cluster with Burkert dark matter halo loses almost all its dark matter to tidal stripping, and starts losing stars at the end of our simulations. The results of this paper reinforce our conclusion in Paper I that current observations of globular clusters are consistent with the primordial picture of globular cluster formation.

In a series of two papers, we test the primordial scenario of globular cluster formation using results of high-resolutions N-body simulations. In this first paper we study the initial relaxation of a stellar core inside a live dark matter minihalo in the early universe. Our dark-matter dominated globular clusters show features which are usually attributed to the action of the tidal field of the host galaxy. Among them are the presence of an apparent cutoff ("tidal radius") or of a "break" in the outer parts of the radial surface brightness profile, and a flat line-of-sight velocity dispersion profile in the outskirts of the cluster. The apparent mass-to-light ratios of our hybrid (stars + dark matter) globular clusters are very close to those of purely stellar clusters. We suggest that additional observational evidence such as the presence of obvious tidal tails is required to rule out the presence of significant amounts of dark matter in present day globular clusters.

We have used high resolution WFPC2-HST and wide field ground-based observations to construct a catalog of blue straggler stars (BSS) which spans the entire radial extent of the globular cluster NGC 6752. The BSS sample is the most extensive ever obtained for this cluster. Though NGC 6752 is a high density cluster with a large binary population, we found that its BSS content is surprisingly low: the specific number of BSS is among the lowest ever measured in a cluster. The BSS distribution is highly peaked in the cluster center, shows a rapid decrease at intermediate radii, and finally rises again at larger distances. This distribution closely resembles those observed in M3 and 47Tuc by Ferraro et al. (1993), Ferraro et al. (2003c). To date, BSS surveys covering the central regions with HST and the outer regions with wide field CCD ground-based observations have been performed for only these three clusters. Despite the different dynamical properties, a bimodal radial distribution has been found in each. A detailed comparison of observed BSS luminosity and temperature distributions with theoretical models reveals a population of luminous, hot BSS which is not easily interpreted.

We investigate a simple model for globular cluster (GC) formation. We simulate the violent relaxation of initially homogeneous isothermal stellar spheres and show that it leads to the formation of clusters with radial density profiles that match the observed profiles of GCs. The best match is achieved for dynamically unevolved clusters. In this model, all the observed correlations between global GC parameters are accurately reproduced if one assumes that all the clusters initially had the same value of the stellar density and the velocity dispersion. This suggests that the gas which formed GCs had the same values of density and temperature throughout the universe.

The recent observation (Ferraro et al. 2003b) of the blue straggler population in 47 Tucanae gives the first detailed characterization of their spatial distribution in the cluster over its entire volume. Relative to the light distribution, blue stragglers appear to be overabundant in the core and at large radii. The observed surface density profile shows a central peak, a zone of avoidance and a rise beyond twenty core radii. In light of these findings we explored the evolution of blue stragglers mimicking their dynamics in a multi-mass King model for 47 Tucanae. We find that the observed spatial distribution can not be explained within a purely collisional scenario in which blue stragglers are generated exclusively in the core through direct mergers. An excellent fit is obtained if we require that a sizable fraction of blue stragglers is generated in the peripheral regions of the cluster inside primordial binaries that evolve in isolation experiencing mass-transfer.

We explore the sensitivity in open cluster ages obtained by the lithium depletion boundary (LDB) technique to the stellar model input physics. The LDB age technique is limited to open clusters with ages ranging from 20 to 200 Myr. Effective 1-sig errors in the LDB technique due to uncertain input physics are roughly 3% at the oldest age increasing to 8% at the youngest age. Bolometric correction uncertainties add an additional 10 to 6% error to the LDB age technique for old and young clusters, respectively. Rotation rates matching the observed fastest rotators in the Pleiades affect LDB ages by less than 2%. The range of rotation rates in an open cluster are expected to ``smear'' the LDB location by only 0.02 mag for a Pleiades age cluster increasing to 0.06 mag for a 20 Myr cluster. Thus, the observational error of locating the LDB (~7-10%) and the bolometric correction uncertainty currently dominate the error in LDB ages. For our base case, we formally derive a LDB age of 148 +- 19 Myr for the Pleiades, where the error includes 8, 3, and 9% contributions from observational, theoretical, and bolometric correction sources, respectively. A maximally plausible 0.3 magnitude shift in the I-band bolometric correction to reconcile main sequence isochrone fits with the observed (V-I) color for the low mass Pleiades members results in an age of 126 +- 11 Myr, where the error includes observational and theoretical errors only. Upper main-sequence-fitting ages that do not include convective core overshoot for the Pleiades (~75 Myr) are ruled out by the LDB age technique.

We have used high resolution WFPC2-HST and wide field ground-based observations to construct a catalog of blue straggler stars (BSS) in the globular cluster 47 Tuc spanning the entire radial extent of the cluster. The BSS distribution is highly peaked in the cluster center, rapidly decreases at intermediate radii, and finally rises again at larger radii. The observed distribution closely resembles that discovered in M3 by Ferraro et al (1993,1997). To date, complete BSS surveys covering the full radial extent (HST in the center and wide field CCD ground based observations of the exterior) have been performed for only these two clusters. Both show a bimodal radial distribution, despite their different dynamical properties. BSS surveys covering the full spatial extent of more globular clusters are clearly required to determine how common bimodality is and what its consequence is for theories of BSS formation and cluster dynamics.

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.

The central regions of six Galactic Globular Clusters (GGCs) (M3, M80, M10, M13, M92 and NGC 288) have been imaged using HST-WFPC2 and the ultraviolet (UV) filters (F255W, F336W). The selected sample covers a large range in both central density and metallicity ([Fe/H]). In this paper, we present a direct cluster-to-cluster comparison of the Blue Stragglers Stars (BSS) population as selected from (m_{255},m_{255}-m_{336}) Color Magnitude Diagrams (CMDs). We have found:(a) BSS in three of the clusters (M3, M80, M92) are much more concentrated toward the center of the cluster than the red giants; because of the smaller BSS samples for the other clusters we can only note that the BSS radial distributions are consistent with central concentration; (b) the specific frequency of BSS varies greatly from cluster to cluster. The most interesting result is that the two clusters with largest BSS specific frequency are at the central density extremes of our sample: NGC 288 (lowest central density) and M80 (highest). This evidence together with the comparison with theoretical collisional models suggests that both stellar interactions in high density cluster cores and at least one other alternate channel operating low density GGCs play an important role in the production of BSS. We also note a possible connection between HB morphology and blue straggler luminosity functions in these six clusters.

We report new measures of radial velocities and rotation rates (v sin i) for 51 F and early-G stars in the open cluster NGC2516, and combine these with previously published data. From high signal-to-noise spectra of two stars, we show that NGC2516 has a relative iron abundance with respect to the Pleiades of delta([Fe/H])= +0.04 +/- 0.07 at the canonical reddening of E(B - V) = 0.12, in contrast to previous photometric studies that placed the cluster 0.2 to 0.4 dex below solar. We construct a color-magnitude diagram based on radial velocity members, and explore the sensitivity of photometric determinations of the metallicity and distance to assumed values of the reddening. For a metal abundance near solar, the Hipparcos distance to NGC2516 is probably underestimated. Finally, we show that the distribution of rotation rates and X-ray emission does not differ greatly from that of the Pleiades, when allowance is made for the somewhat older age of NGC2516.

We performed high-resolution simulations of two stellar collisions relevant for stars in globular clusters. We considered one head-on collision and one off-axis collision between two 0.6 M_sun main sequence stars. We show that a resolution of about 100 000 particles is sufficient for most studies of the structure and evolution of blue stragglers. We demonstrate conclusively that collision products between main-sequence stars in globular clusters do not have surface convection zones larger than 0.004 M_sun after the collision, nor do they develop convection zones during the `pre-main-sequence' thermal relaxation phase of their post-collision evolution. Therefore, any mechanism which requires a surface convection zone (i.e. chemical mixing or angular momentum loss via a magnetic wind) cannot operate in these stars. We show that no disk of material surrounding the collision product is produced in off-axis collisions. The lack of both a convection zone and a disk proves a continuing problem for the angular momentum evolution of blue stragglers in globular clusters.

Collisions of main sequence stars occur frequently in dense star clusters. In open and globular clusters, these collisions produce merger remnants that may be observed as blue stragglers. Detailed theoretical models of this process require lengthy hydrodynamic computations in three dimensions. However, a less computationally expensive approach, which we present here, is to approximate the merger process (including shock heating, hydrodynamic mixing, mass ejection, and angular momentum transfer) with simple algorithms based on conservation laws and a basic qualitative understanding of the hydrodynamics. These algorithms have been fine tuned through comparisons with the results of our previous hydrodynamic simulations. We find that the thermodynamic and chemical composition profiles of our simple models agree very well with those from recent SPH (smoothed particle hydrodynamics) calculations of stellar collisions, and the subsequent stellar evolution of our simple models also matches closely that of the more accurate hydrodynamic models. Our algorithms have been implemented in an easy to use software package, which we are making publicly available (see http://vassun.vassar.edu/~lombardi/mmas/). This software could be used in combination with realistic dynamical simulations of star clusters that must take into account stellar collisions.

We present FUSE H Lyman series spectroscopy of the hot white dwarf companion to the 4th magnitude A1III star Beta Crateris, which shows that is has an unusually low mass, 0.43Msun, and has almost certainly evolved through binary interaction. This system could be a long-sought remnant of Algol-type evolution, although radial velocity measurements appear to show that the pair are not close. Instead, micro-variations in Beta Crt's proper motion as measured by Hipparcos suggests that the period could be as high as 10 years. However, a low mass white dwarf in a system with a period >3 years is difficult to explain by conventional models for binary evolution. We speculate on alternative models for the evolution of this system which involve an eccentric binary or multiple components.

We apply the angular momentum loss rates inferred from open cluster stars to the evolution of cataclysmic variables (CVs). We show that the angular momentum prescriptions used in earlier CV studies are inconsistent with the measured rotation data in open clusters. The timescale for angular momentum loss above the fully convective boundary is ~ 2 orders of magnitude longer than inferred from the older model, and the observed angular momentum loss properties show no evidence for a change in a behavior at the fully convective boundary. This provides evidence against the hypothesis that the period gap is caused by an abrupt change in the angular momentum loss law when the secondary becomes fully convective. It also implies that the timescale for CV evolution is much longer than it was than previously thought, comparable to a Hubble time. For the same reason, it will be more difficult to produce CVs from the products of common envelope evolution and implies a lower space density of CVs. The empirical loss law is consistent with the observed period minimum (1.3 hours) contrary to the minimum predicted by angular momentum loss due to gravitational radiation alone (1.1 hours). We introduce a method to infer the time-averaged mass accretion rate and derive mass-period relation for different evolutionary states of the secondary. The mass-period relationship is more consistent with evolved secondaries than with unevolved secondaries above the period gap. Implications for the CV period gap are discussed, including the possibility that two populations of secondaries could produce the gap.

We continue our exploration of collisionally merged stars in the blue straggler region of the color-magnitude diagram. We report the results of new SPH calculations of parabolic collisions between two main-sequence stars, with the initial structure and composition profiles of the parent stars having been determined from stellar evolution calculations. Parallelization of the SPH code has permitted much higher numerical resolution of the hydrodynamics. We also present evolutionary tracks for the resulting collision products, which emerge as rapidly rotating blue stragglers. The rotating collision products are brighter, bluer and remain on the main sequence longer than their non-rotating counterparts. In addition, they retain their rapid rotation rates throughout their main sequence lifetime. Rotationally-induced mixing strongly affects the evolution of the collision products, although it is not sufficient to mix the entire star. We discuss the implications of these results for studies of blue straggler populations in clusters. This work shows that off-axis collision products cannot become blue stragglers unless they lose a large fraction of their initial angular momentum. The mechanism for this loss is not apparent, although some possibilities are discussed.

We present theoretical models of the angular momentum evolution of very low mass stars (0.1 - 0.5 M_sun) and solar analogues (0.6 - 1.1 M_sun). We investigate the effect of rotation on the effective temperature and luminosity of these stars. We find that the decrease in T_eff and L can be significant at the higher end of our mass range, but becomes negligible below 0.4 M_sun. Formulae for relating T_eff to mass and v_rot are presented. We compare our models to rotational data from young open clusters of different ages to infer the rotational history of low mass stars, and the dependence of initial conditions and rotational evolution on mass. We find that the qualitative conclusions for stars below 0.6 M_sun do not depend on the assumptions about internal angular momentum transport, which makes these low mass stars ideal candidates for the study of the angular momentum loss law and distribution of initial conditions. We find that neither models with solid body nor differential rotation can simultaneously reproduce the observed stellar spin down in the 0.6 to 1.1 M_sun mass range and for stars between 0.1 and 0.6 M_sun. The most likely explanation is that the saturation threshold drops more steeply at low masses than would be predicted with a simple Rossby scaling. In young clusters there is a systematic increase in the mean rotation rate with decreased temperature below 3500 K (0.4 M_sun). This suggests either inefficient angular momentum loss or mass-dependent initial conditions for stars near the fully convective boundary. (abridged)

We have obtained high-resolution spectra of 89 M dwarf members of the Pleiades and Hyades and have derived radial velocities, H-alpha equivalent widths, and spectroscopic rotational velocities for these stars. Typical masses of the newly-observed Pleiades and Hyades stars are ~ 0.4 M_{\sun} and ~ 0.2 M_{\sun}, respectively. We combine our new observations with previously published data to explore the rotational evolution of young stars with M < 0.4 M_\sun. The average rotation rate in the Hyades (age 600 Myr) is about 0.4 that of the Pleiades (110 Myr), and the mean equivalent widths of H-alpha are also lower. As found in previous studies, the correlation between rotation and chromospheric activity is identical in both clusters, implying that the lower activity in the Hyades is a result of the lower rotation rates. We show that a simple scaling of the Pleiades rotational distribution for M \leq 0.4 M_{\sun}, corrected for the effects of structural evolution, matches that of the Hyades if the average angular momentum loss from the Pleiades to the Hyades age is factor of \approx 6. This suggests that the distribution of initial angular momenta and disk-locking lifetimes for the lowest mass stars was similar in both clusters. We argue that this result provides further evidence for a saturation of the angular momentum loss rate at high rotational velocities.

The rotation of horizontal branch stars places important constraints on angular momentum evolution in evolved stars and therefore rotational mixing on the giant branch. Prompted by new observations of rotation rates of horizontal branch stars, we calculate simple models for the angular momentum evolution of a globular cluster giant star from the base of the giant branch to the star's appearance on the horizontal branch. We include mass loss, and infer the accompanied loss of angular momentum for each of four assumptions about the internal angular momentum profile. These models are compared to observations of horizontal branch rotation rates in M13. We find that rapid rotation on the horizontal branch can be reconciled with slow solid body main sequence rotation if giant branch stars have differential rotation in their convective envelopes and a rapidly rotating core, which is then followed by a redistribution of angular momentum on the horizontal branch. We discuss the physical reasons why these very different properties relative to the solar case may exist in giants. Rapid rotation in the core of the main sequence precursors of the rapidly rotating horizontal branch star, or an angular momentum source on the giant branch is required for all cases if the rotational velocity of turnoff stars is less than 4 km s$^{-1}$. We suggest that the observed range in rotation rates on the horizontal branch is caused by internal angular momentum redistribution which occurs on a timescale comparable to the evolution of the stars on the horizontal branch. The apparent lack of rapid horizontal branch rotators hotter than 12 000 K in M13 could be a consequence of gravitational settling, which inhibits internal angular momentum transport. Alternative explanations and observational tests are discussed.

We investigate the effects of changes in the blue straggler formation rate in globular clusters on the blue straggler distribution in the color-magnitude diagram. We find that the blue straggler distribution is highly sensitive to the past formation rate. Comparing our models to new UBV observations of a region close to the core of 47 Tucanae suggests that this cluster may have stopped forming blue straggler formation several Gyr ago. This cessation of formation can be associated with an epoch of primordial binary burning which has been invoked in other clusters to infer the imminence of core collapse.

We present theoretical models of rotating low mass stars (0.1 - 1.0 \Msun) to demonstrate the effect of rotation on the effective temperature and luminosity of stars. The range of rotation rates in our models corresponds to the observed rotation rates in young low mass stars. Rotation lowers the effective temperature and luminosity of the models relative to standard models of the same mass and composition. We find that the decrease in T$_{eff}$ and L can be significant at the higher end of our mass range, but becomes small below 0.4 \Msun. The effects of different assumptions about internal angular momentum transport are discussed. Formulae for relating T$_{eff}$ to mass and v$_{rot}$ are presented. We demonstrate that the kinetic energy of rotation is not a significant contribution to the luminosity of low mass stars.

Using the Yale stellar evolution code, we have calculated theoretical models for nearby stars with planetary-mass companions in short-period nearly circular orbits: 51 Pegasi, Tau Bootis, Upsilon Andromedae, Rho Cancri, and Rho Coronae Borealis. We present tables listing key stellar parameters such as mass, radius, age, and size of the convective envelope as a function of the observable parameters (luminosity, effective temperature, and metallicity), as well as the unknown helium fraction. For each star we construct best models based on recently published spectroscopic data and the present understanding of galactic chemical evolution. We discuss our results in the context of planet formation theory, and, in particular, tidal dissipation effects and stellar metallicity enhancements.

We have performed a series of systematic tests to evaluate the effects of spurious transport in three-dimensional smoothed particle hydrodynamics (SPH) calculations. Our tests investigate (i) particle diffusion, (ii) shock heating, (iii) numerical viscosity, and (iv) angular momentum transport. The results are useful for quantifying the accuracy of the SPH scheme, especially for problems where shear flows or shocks are present, as well as for problems where true hydrodynamic mixing is relevant. We examine the different forms of artificial viscosity (AV) which have been proposed by Monaghan, by Hernquist & Katz, and by Balsara. For each form, our tests suggest a single set of values for the AV parameters $\alpha$ and $\beta$ (coefficients of the linear and quadratic terms) which are appropriate in a large number of situations. We also discuss how these parameters should be adjusted depending on the goals of the particular application. We find that both the Hernquist & Katz and Balsara forms introduce relatively small amounts of numerical viscosity. Furthermore, both Monaghan's and Balsara's AV do well at treating shocks and at limiting the amount of spurious mixing. For these reasons, we endorse the Balsara AV for use in a broad range of applications.

A primary production mechanism for blue stragglers in globular clusters is thought to be collisionally-induced mergers, perhaps mediated by dynamical encounters involving binary stars. We model the formation and evolution of such blue stragglers, and produce theoretical distributions of them in the colour-magnitude diagram. We use a crude representation of cluster dynamics and detailed binary-single star encounter simulations to produce cross sections and rates for a variety of collisions. The results of the collisions are determined based on SPH simulations of realistic star models. The evolution of the collision products are then followed in detail. We use our results to explore the effects of a variety of input assumptions on the number and kind of blue stragglers created by collisions. In particular, we describe the changes in the blue straggler distribution that result from using realistic collision products rather than the ``fully-mixed'' assumption, and from changes in assumptions about the number and orbital period distribution of the primordial binary population. We then apply our models to existing data from the core of M3, where the large blue straggler population is thought to be dominated by collision products. We find that we have difficulty successfully modeling the observed blue stragglers under a single consistent set of assumptions. However, if 3 particularly bright blue stragglers are considered to be part of a different observed population, the remainder can be successfully modeled using realistic encounter products and assuming a 20% binary fraction with plausible period distribution.

We demonstrate the necessity of using realistic stellar models taken from stellar evolution codes, as opposed to polytropes, for starting models in smoothed particle hydrodynamics calculations of collisions between main sequence stars. Evolved stars have mean molecular weight gradients, which affect their entropy profiles and therefore affect how they react during a collision. The structure of stellar collision products of polytrope parent stars is significantly different from that of collision products of realistic parent models. These differences strongly affect the future evolution of the collision products, particularly products of collisions between unequal mass stars which have undergone significant chemical evolution. The use of polytropes as parent star models is likely to result in qualitatively mistaken results for the structure of the collision product.

We explore the evolution of collisionally merged stars in the blue straggler region of the HR diagram. The starting models for our stellar evolution calculations are the results of the smoothed particle hydrodynamics (SPH) simulations of parabolic collisions between main sequence stars performed by Lombardi, Rasio and Shapiro (1996). Since SPH and stellar evolution codes employ different and often contradictory approximations, it is necessary to treat the evolution of these products carefully. The mixture and disparity of the relevant timescales (hydrodynamic, thermal relaxation and nuclear burning) and of the important physical assumptions between the codes makes the combined analysis of the problem challenging, especially during the initial thermal relaxation of the star. In particular, the treatment of convection is important, and semiconvection must be modeled in some detail. The products of seven head-on collisions are evolved through their initial thermal relaxation, and then through the main sequence phase to the base of the giant branch. Their evolutionary tracks are presented. In contrast to the assumptions in previous work, these collision products do not develop substantial convective regions during their thermal relaxation, and therefore are not mixed significantly after the collision.