We introduce the stellar sample of the WIYN Open Cluster Study radial-velocity survey for the rich open cluster NGC 7789 (1.6 Gyr, [Fe/H] = +0.02). This sample lies within an 18$'$ circular radius from the cluster center (10 pc in projection, or about 2 core radii), and includes giants, red-clump stars, blue stragglers, red stragglers, sub-subgiants, and main-sequence stars down to one mag below the turnoff. Our survey began in 2005 and comprises more than 9,000 radial-velocity measurements from the Hydra Multi-Object Spectrograph on the WIYN 3.5m telescope. We identify 564 likely cluster members and present the orbital solutions for 83 cluster binary stars with periods between 1.45$\,$d and 4200$\,$d. From the main-sequence binary solutions we fit a circularization period of 8.22$^{+3.51}_{-1.35}\,$d. We calculate an incompleteness-corrected main-sequence binary frequency of 31% $\pm$ 3% for binaries with periods less than 10$^{4}$ days, similar to other WOCS open clusters of all ages. We detect a blue straggler binary frequency of 31% $\pm$ 15%, consistent with the similarly aged cluster NGC 6819.

To fully understand the diverse population of exoplanets, we must study their early lives within open clusters, the birthplace of most stars with masses $>0.5M_\odot$ (including those currently in the field). Indeed, when we observe planets within clustered environments, we notice highly eccentric and odd systems that suggest the importance of dynamical pathways created by interactions with additional bodies (as in the case of HD 285507b). However, it has proven difficult to investigate these effects, as many current numerical solvers for the multi-scale $N$-body problem are simplified and limited in scope. To remedy this, we aim to create a physically complete computational solution to explore the role of stellar close encounters and interplanetary interactions in producing the observed exoplanet populations for both open cluster stars and field stars. We present a new code, TYCHO, which employs a variety of different computational techniques, including multiple $N$-body integration methods, close encounter handling, modified Monte Carlo scattering experiments, and a variety of empirically informed initial conditions. We discuss the methodology in detail, and its implementation within the AMUSE software framework. Approximately 1$\%$ of our systems are promptly disrupted by star-star encounters contributing to the rogue planets occurrence rate. Additionally, we find that close encounters which perturb long-period planets lead to 38.3$\%$ of solar-system-like planetary systems becoming long-term unstable.

The binary fractions of open and globular clusters yield powerful constraints on their dynamical state and evolutionary history. We apply publicly available Bayesian analysis code to a $UBVRIJHK_{S}$ photometric catalog of the old open cluster NGC 188 to detect and characterize photometric binaries along the cluster main sequence. This technique has the advantage that it self-consistently handles photometric errors, missing data in various bandpasses, and star-by-star prior constraints on cluster membership. Simulations are used to verify uncertainties and quantify selection biases in our analysis, illustrating that among binaries with mass ratios >0.5, we recover the binary fraction to better than 7% in the mean, with no significant dependence on binary fraction and a mild dependence on assumed mass ratio distribution. Using our photometric catalog, we recover the majority (65%$\pm$11%) of spectroscopically identified main sequence binaries, including 8 of the 9 with spectroscopically measured mass ratios. Accounting for incompleteness and systematics, we derive a mass ratio distribution that rises toward lower mass ratios (within our $q >$0.5 analysis domain). We observe a raw binary fraction for solar-type main sequence stars with mass ratios $q >$0.5 of 42%$\pm$4%, independent of the assumed mass ratio distribution to within its uncertainties, consistent with literature values for old open clusters but significantly higher than the field solar-type binary fraction. We confirm that the binaries identified by our method are more concentrated than single stars, in agreement with previous studies, and we demonstrate that the binary nature of those candidates which remain unidentified spectroscopically is strongly supported by photometry from Gaia DR2.

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 investigate the old open cluster M67 using ultraviolet photometric data of Ultra-Violet Imaging Telescope in multi-filter far-UV bands. M67, well known for the presence of several blue straggler stars (BSS), has been put to detailed tests to understand their formation pathways. Currently, there are three accepted formation channels: mass transfer due to Roche-lobe overflow in binary systems, stellar mergers either due to dynamical collisions or through coalescence of close binaries. So far, there had not been any confirmed detection of a white dwarf (WD) companion to any of the BSSs in this cluster. Here, we present the detection of WD companions to 5 bright BSSs in M67. The multiwavelength spectral energy distributions covering 0.12 -11.5 $\mu$m range, were found to require binary spectral fits for 5 BSSs, consisting of a cool (BSS) and a hot companion. The parameters (Luminosity, Temperature, Radius and Mass) of the hot companions suggest them to be WDs with mass in the range 0.2 - 0.35 M$_{\odot}$ with T$_{eff}$ $\sim$ 11000 - 24000 K.

The old open cluster M67, populated with blue straggler stars (BSSs), is a well known test bed to study the BSS formation pathways. Here, we report the first direct detection of a white dwarf (WD) companion to a BSS in M67, using far-UV images from the \textit{Ultra Violet Imaging telescope} ({\it UVIT}) on {\it ASTROSAT}. Near-simultaneous observations in three far-UV bands combined with {\it GALEX}, {\it IUE}, ground and space based photometric data covering 0.14 -11.5 $\mu$m range for WOCS1007 were found to require a binary fit to its spectral energy distribution (SED), consisting of a BSS and a hot companion. On the other hand, a single spectral fit was found to be satisfactory for the SEDs of two other BSSs, WOCS1006 and WOCS2011, with the latter showing a deficient far-UV flux. The hot companion of WOCS1007 is found to have a T$_{eff}$ $\sim$ 13250-13750K and a radius of 0.09$\pm$0.01 R$_\odot$. A comparison with WD models suggests it to be a low mass WD ($\sim$ 0.18$M_\odot$), in agreement with the kinematic mass from the literature. As a low mass WD ($<$ 0.4$M_\odot$) necessitates formation through mass transfer (MT) in close binaries, WOCS1007 with a known period of 4.2 days along with its fast rotation, is likely to be formed by a case A or case B binary evolution.

Astrophysics lies at the crossroads of big datasets (such as the Large Synoptic Survey Telescope and Gaia), open source software to visualize and interpret high dimensional datasets (such as Glue, WorldWide Telescope, and OpenSpace), and uniquely skilled software engineers who bridge data science and research fields. At the same time, more than 4,000 planetariums across the globe immerse millions of visitors in scientific data. We have identified the potential for critical synergy across data, software, hardware, locations, and content that -- if prioritized over the next decade -- will drive discovery in astronomical research. Planetariums can and should be used for the advancement of scientific research. Current facilities such as the Hayden Planetarium in New York City, Adler Planetarium in Chicago, Morrison Planetarium in San Francisco, the Iziko Planetarium and Digital Dome Research Consortium in Cape Town, and Visualization Center C in Norrkoping are already developing software which ingests catalogs of astronomical and multi-disciplinary data critical for exploration research primarily for the purpose of creating scientific storylines for the general public. We propose a transformative model whereby scientists become the audience and explorers in planetariums, utilizing software for their own investigative purposes. In this manner, research benefits from the authentic and unique experience of data immersion contained in an environment bathed in context and equipped for collaboration. Consequently, in this white paper we argue that over the next decade the research astronomy community should partner with planetariums to create visualization-based research opportunities for the field. Realizing this vision will require new investments in software and human capital.

At least $70\%$ of massive OBA-type stars reside in binary or higher-order systems. The dynamical evolution of these systems can lend insight into the origins of extreme phenomena such as X-ray binaries and gravitational wave sources. In one such dynamical process, the Eccentric Kozai-Lidov (EKL) Mechanism, a third companion star alters the secular evolution of a binary system. For dynamical stability, these triple systems must have a hierarchical configuration. We explore the effects of a distant third companion's gravitational perturbations on a massive binary's orbital configuration before significant stellar evolution has taken place ($\leq 10$ Myr). We include tidal dissipation and general relativistic precession. With large ($38,000$ total) Monte-Carlo realizations of massive hierarchical triples, we characterize imprints of the birth conditions on the final orbital distributions. Specifically, we find that the final eccentricity distribution over the range $0.1-0.7$ is an excellent indicator of its birth distribution. Furthermore, we find that the period distributions have a similar mapping for wide orbits. Finally, we demonstrate that the observed period distribution for approximately $10$ Myr-old massive stars is consistent with EKL evolution.

We report the discovery of a hot white dwarf (WD) companion to a blue straggler star (BSS) in the globular cluster (GC) NGC 5466, based on observations from the Ultra-Violet Imaging Telescope (UVIT) on board AstroSat. The Spectral Energy Distribution (SED) of the Far-UV detected BSS NH 84 was constructed by combining the flux measurements from 4 filters of UVIT, with GALEX, GAIA and other ground-based observations. The SED of NH 84 reveals the presence of a hot companion to the BSS. The temperature and radius of the BSS (T$_{\mathrm{eff}} = 8000^{+1000}_{-250}$ K, R/R$_\odot = 1.44 \pm 0.05$) derived from Gemini spectra and SED fitting using Kurucz atmospheric models are consistent with each other. The temperature and radius of the hotter companion of NH 84 (T$_{\mathrm{eff}} = 32,000 \pm 2000$ K, R/R$_\odot = 0.021 \pm 0.007$) derived by fitting Koester WD models to the SED suggest that it is likely to be a hot WD. The radial velocity derived from the spectra along with the proper motion from GAIA DR2 confirms NH 84 to be a kinematic member of the cluster. This is the second detection of a BSS-WD candidate in a GC, and the first in the outskirts of a low density GC. The location of this BSS in NGC 5466 along with its dynamical age supports the mass-transfer pathway for BSS formation in low density environments.

About a century ago, Jeans (1919) discovered that if binary stars reach a state approximating energy equipartition, for example through many dynamical encounters that exchange energy, their eccentricity distribution can be described by : dN/de = 2e. This is referred to as the thermal eccentricity distribution, and has been widely used for initial conditions in theoretical investigations of binary stars. However, observations suggest that the eccentricity distributions of most observed binaries, and particularly those with masses < 5 Msun, are flatter than thermal and follow more closely to a uniform distribution. Nonetheless, it is often argued that dynamical interactions in a star cluster would quickly thermalize the binaries, which could justify imposing a thermal eccentricity distribution at birth for all binaries. In this paper we investigate the validity of this assumption. We develop our own rapid semi-analytic model for binary evolution in star clusters, and also compare with detailed N-body and Monte Carlo star cluster models. We show that, for nearly all binaries, dynamical encounters fail to convert an initially uniform eccentricity distribution to thermal within a star cluster's lifetime. Thus, if a thermal eccentricity distribution is observed, it is likely imprinted upon formation rather than through subsequent long-term dynamical processing. Theoretical investigations that initialize all binaries with a thermal distribution will make incorrect predictions for the evolution of the binary population. Such models may overpredict the merger rate for binaries with modest orbital separations by a factor of about two.

NGC 3532 is an extremely rich open cluster embedded in the Galactic disc, hitherto lacking a comprehensive, documented membership list. We provide membership probabilities from new radial velocity observations of solar-type and low-mass stars in NGC 3532, in part as a prelude to a subsequent study of stellar rotation in the cluster. Using extant optical and infra-red photometry we constructed a preliminary photometric membership catalogue, consisting of 2230 dwarf and turn-off stars. We selected 1060 of these for observation with the AAOmega spectrograph at the Anglo-Australian Telescope and 391 stars for observations with the Hydra-South spectrograph at the Victor Blanco Telescope, obtaining spectroscopic observations over a decade for 145 stars. We measured radial velocities for our targets through cross-correlation with model spectra and standard stars, and supplemented them with radial velocities for 433 additional stars from the literature. We also measured log g, Teff, and [Fe/H] from the AAOmega spectra. Together with proper motions from Gaia DR2 we find 660 exclusive members. The members are distributed across the whole cluster sequence, from giant stars to M dwarfs, making NGC 3532 one of the richest Galactic open clusters known to date, on par with the Pleiades. From further spectroscopic analysis of 153 dwarf members we find the metallicity to be marginally sub-solar, with [Fe/H]=-0.07. Exploiting trigonometric parallax measurements from Gaia DR2 we find a distance of $484^{+35}_{-30}$ pc. Based on the membership we provide an empirical cluster sequence in multiple photometric passbands. A comparison of the photometry of the measured cluster members with several recent model isochrones enables us to confirm the 300 Myr cluster age. However, all of the models evince departures from the cluster sequence in particular regions, especially in the lower mass range. (abridged)

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 determine the properties of the binary star V106 in the old open cluster NGC6791. We identify the system to be a blue straggler cluster member by using a combination of ground-based and Kepler photometry and multi-epoch spectroscopy. The properties of the primary component are found to be $M_{\rm p}\sim1.67 \rm M_{\odot}$, more massive than the cluster turn-off, with $R_{\rm p}\sim1.91 \rm R_{\odot}$ and $T_{\rm eff}=7110\pm100$ K. The secondary component is highly oversized and overluminous for its low mass with $M_{\rm s}\sim0.182 \rm M_{\odot}$, $R_{\rm s}\sim0.864 \rm R_{\odot}$ and $T_{\rm eff}=6875\pm200$ K. We identify this secondary star as a bloated (proto) extremely low-mass helium white dwarf. These properties of V106 suggest that it represents a typical Algol-paradox system and that it evolved through a mass-transfer phase which provides insight into its past evolution. We present a detailed binary stellar evolution model for the formation of V106 using the MESA code and find that the mass-transfer phase only ceased about 40 Myr ago. Due to the short orbital period (P=1.4463 d) another mass-transfer phase is unavoidable once the current primary star evolves towards the red giant phase. We argue that V106 will evolve through a common-envelope phase within the next 100 Myr and merge to become a single over-massive giant. The high mass will make it appear young for its true age, which is revealed by the cluster properties. Therefore, V106 is potentially a prototype progenitor of old field giants masquerading as young.

In this paper, the fourth in the series, we continue our study of combinatorics in chaotic Newtonian dynamics. We focus once again on the chaotic four-body problem in Newtonian gravity assuming finite-sized particles, and interactions that produce direct collisions between any two particles. Our long-term goal is to predict the probability of a given collision event occurring over the course of an interaction, as a function of the numbers and properties of the particles. In previous papers, we varied the number of interacting particles, as well as the distributions of particle radii and masses. Here, we refine the methods developed in these preceding studies, and arrive at a final and robust methodology that can be used to study collisional dynamics in a variety of astrophysical contexts, ranging from stars in star clusters, galaxies in galaxy groups and clusters and even the collisional growth of planetesimals in protoplanetary disks. We further present and refine the concept of a Collision Rate Diagram (CRD), the primary tool we use to quantify the relative rates for different collision scenarios to occur. The agreement between our final theoretical model and the results of numerical scattering simulations is excellent.

We assess the contribution of dynamical hardening by direct three-body scattering interactions to the rate of stellar-mass black hole binary (BHB) mergers in galactic nuclei. We derive an analytic model for the single-binary encounter rate in a nucleus with spherical and disk components hosting a super-massive black hole (SMBH). We determine the total number of encounters $N_{\rm GW}$ needed to harden a BHB to the point that inspiral due to gravitational wave emission occurs before the next three-body scattering event. This is done independently for both the spherical and disk components. Using a Monte Carlo approach, we refine our calculations for $N_{\rm GW}$ to include gravitational wave emission between scattering events. For astrophysically plausible models we find that typically $N_{\rm GW} \lesssim$ 10. We find two separate regimes for the efficient dynamical hardening of BHBs: (1) spherical star clusters with high central densities, low velocity dispersions and no significant Keplerian component; and (2) migration traps in disks around SMBHs lacking any significant spherical stellar component in the vicinity of the migration trap, which is expected due to effective orbital inclination reduction of any spherical population by the disk. We also find a weak correlation between the ratio of the second-order velocity moment to velocity dispersion in galactic nuclei and the rate of BHB mergers, where this ratio is a proxy for the ratio between the rotation- and dispersion-supported components. Because disks enforce planar interactions that are efficient in hardening BHBs, particularly in migration traps, they have high merger rates that can contribute significantly to the rate of BHB mergers detected by the advanced Laser Interferometer Gravitational-Wave Observatory.

In this paper, the third in the series, we continue our study of combinatorics in chaotic Newtonian dynamics. We study the chaotic four-body problem in Newtonian gravity assuming finite-sized particles, and we focus on interactions that produce direct collisions between any two stars. Our long-term goal is to construct an equation that gives the probability of a given collision event occurring over the course of the interaction, as a function of the total encounter energy and angular momentum as well as the numbers and properties of the particles. In previous papers, we varied the number of interacting particles and the distribution of particle radii, for all equal mass particles. Here, we focus on the effects of different combinations of particle masses. We develop an analytic formalism for calculating the time-scales for different collision scenarios to occur. Our analytic time-scales reproduce the simulated time-scales when gravitational focusing is included. We present a method for calculating the relative rates for different types of collisions to occur, assuming two different limits for the particle orbits; radial and tangential. These limits yield relative collision probabilities that bracket the probabilities we obtain directly from numerical scattering experiments, and are designed to reveal important information about the (time-averaged) trajectories of the particles as a function of the interaction parameters. Finally, we present a Collision Rate Diagram (CRD), which directly compares the predictions of our analytic rates to the simulations and quantifies the quality of the agreement. The CRD will facilitate refining our analytic collision rates in future work, as we expand in to the remaining parameter space.

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.

Sub-subgiant stars (SSGs) lie to the red of the main-sequence and fainter than the red giant branch in cluster color-magnitude diagrams (CMDs), a region not easily populated by standard stellar evolution pathways. While there has been speculation on what mechanisms may create these unusual stars, no well-developed theory exists to explain their origins. Here we discuss three hypotheses of SSG formation: (1) mass transfer in a binary system, (2) stripping of a subgiant's envelope, perhaps during a dynamical encounter, and (3) reduced luminosity due to magnetic fields that lower convective efficiency and produce large star spots. Using the stellar evolution code MESA, we develop evolutionary tracks for each of these hypotheses, and compare the expected stellar and orbital properties of these models with six known SSGs in the two open clusters M67 and NGC 6791. All three of these mechanisms can create stars or binary systems in the SSG CMD domain. We also calculate the frequency with which each of these mechanisms may create SSG systems, and find that the magnetic field hypothesis is expected to create SSGs with the highest frequency in open clusters. Mass transfer and envelope stripping have lower expected formation frequencies, but may nevertheless create occasional SSGs in open clusters. They may also be important mechanisms to create SSGs in higher mass globular clusters.

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.

In this paper, we study the chaotic four-body problem in Newtonian gravity. Assuming point particles and total encounter energies $\le$ 0, the problem has three possible outcomes. We describe each outcome as a series of discrete transformations in energy space, using the diagrams first presented in Leigh \& Geller (2012; see the Appendix). Furthermore, we develop a formalism for calculating probabilities for these outcomes to occur, expressed using the density of escape configurations per unit energy, and based on the Monaghan description originally developed for the three-body problem. We compare this analytic formalism to results from a series of binary-binary encounters with identical point particles, simulated using the \texttt{FEWBODY} code. Each of our three encounter outcomes produces a unique velocity distribution for the escaping star(s). Thus, these distributions can potentially be used to constrain the origins of dynamically-formed populations, via a direct comparison between the predicted and observed velocity distributions. Finally, we show that, for encounters that form stable triples, the simulated single star escape velocity distributions are the same as for the three-body problem. This is also the case for the other two encounter outcomes, but only at low virial ratios. This suggests that single and binary stars processed via single-binary and binary-binary encounters in dense star clusters should have a unique velocity distribution relative to the underlying Maxwellian distribution (provided the relaxation time is sufficiently long), which can be calculated analytically.

What is the difference between a long-lived unstable (or quasi-stable) multiple star system and a bona fide star cluster? In this paper, we present a possible framework to address this question, by studying the distributions of disruption times for chaotic gravitational encounters as a function of the number of interacting particles. To this end, we perform a series of numerical scattering experiments with the \texttt{FEWBODY} code, to calculate the distributions of disruption times as a function of both the particle number N and the virial coefficient k. The subsequent distributions are fit with a physically-motivated function, consisting of an initial exponential decay followed by a very slowly decreasing tail at long encounter times due to long-lived quasi-stable encounters. We find three primary features characteristic of the calculated distributions of disruption times. These are: (1) the system half-life increases with increasing particle number, (2) the fraction of long-lived quasi-stable encounters increases with increasing particle number and (3) both the system half-life and the fraction of quasi-stable encounters increase with decreasing virial coefficient. We discuss the significance of our results for collisional dynamics, and consider the extrapolation of our results to larger-N systems. We suggest that this could potentially offer a clear and unambiguous distinction between star clusters and (unstable or quasi-stable) multiple star systems. Although we are limited by very small-number statistics, our results tentatively suggest that (for our assumptions) this transition occurs at a critical particle number of order 100.

Stars in star clusters are thought to form in a single burst from a common progenitor cloud of molecular gas. However, massive, old globular clusters -- with ages greater than 10 billion years and masses of several hundred thousand solar masses -- often harbour multiple stellar populations, indicating that more than one star-forming event occurred during their lifetimes. Colliding stellar winds from late-stage, asymptotic-giant-branch stars are often invoked as second-generation star-formation trigger. The initial cluster masses should be at least 10 times more massive than they are today for this to work. However, large populations of clusters with masses greater than a few million solar masses are not found in the local Universe. Here we report on three 1-2 billion-year-old, massive star clusters in the Magellanic Clouds, which show clear evidence of burst-like star formation that occurred a few hundred million years after their initial formation era. We show that such clusters could accrete sufficient gas reservoirs to form new stars if the clusters orbited in their host galaxies' gaseous discs throughout the period between their initial formation and the more recent bursts of star formation. This may eventually give rise to the ubiquitous multiple stellar populations in globular clusters.

As part of our study of the old (~2.5 Gyr) open cluster NGC 6819 in the Kepler field, we present photometric (Kepler and ground-based BVRcIc) and spectroscopic observations of the detached eclipsing binary WOCS 24009 (Auner 665; KIC 5023948) with a short orbital period of 3.6 days. WOCS 24009 is a triple-lined system, and we verify that the brightest star is physically orbiting the eclipsing binary using radial velocities and eclipse timing variations. The eclipsing binary components have masses M_B =1.090+/-0.010 Msun and M_C =1.075+/-0.013 Msun, and radii R_B =1.095+/-0.007 Rsun and R_C =1.057+/-0.008 Rsun. The bright non-eclipsing star resides at the cluster turnoff, and ultimately its mass will directly constrain the turnoff mass: our preliminary determination is M_A =1.251+/-0.057 Msun. A careful examination of the light curves indicates that the fainter star in the eclipsing binary undergoes a very brief period of total eclipse, which enables us to precisely decompose the light of the three stars and place them in the color-magnitude diagram. We also present improved analysis of two previously discussed detached eclipsing stars in NGC 6819 (WOCS 40007 and WOCS 23009) en route to a combined determination of the cluster's distance modulus (m-M)_V=12.38+/-0.04. Because this paper significantly increases the number of measured stars in the cluster, we can better constrain the age of the color-magnitude diagram to be 2.21+/-0.10+/-0.20 Gyr. Additionally, using all measured eclipsing binary star masses and radii, we constrain the age to 2.38+/-0.05+/-0.22 Gyr. The quoted uncertainties are estimates of measurement and systematic uncertainties (due to model physics differences and metal content), respectively.

Binary mass transfer is at the forefront of some of the most exciting puzzles of modern astrophysics, including Type Ia supernovae, gamma-ray bursts, and the formation of most observed exotic stellar populations. Typically, the evolution is assumed to proceed in isolation, even in dense stellar environments such as star clusters. In this paper, we test the validity of this assumption via the analysis of a large grid of binary evolution models simulated with the SeBa code. For every binary, we calculate analytically the mean time until another single or binary star comes within the mean separation of the mass-transferring binary, and compare this time-scale to the mean time for stable mass transfer to occur. We then derive the probability for each respective binary to experience a direct dynamical interruption. The resulting probability distribution can be integrated to give an estimate for the fraction of binaries undergoing mass transfer that are expected to be disrupted as a function of the host cluster properties. We find that for lower-mass clusters ($\lesssim 10^4$ M$_{\odot}$), on the order of a few to a few tens of percent of binaries undergoing mass-transfer are expected to be interrupted by an interloping single, or more often binary, star, over the course of the cluster lifetime, whereas in more massive globular clusters we expect $\ll$ 1% to be interrupted. Furthermore, using numerical scattering experiments performed with the FEWBODY code, we show that the probability of interruption increases if perturbative fly-bys are considered as well, by a factor $\sim 2$.

(Abridged) We present the results of our ongoing radial-velocity (RV) survey of the old (7 Gyr) open cluster NGC 188. Our WIYN 3.5m data set spans a time baseline of 11 years, a magnitude range of 12<=V<=16.5 (1.18-0.94 MSun), and a 1 deg. diameter region on the sky. With the addition of a Dominion Astrophysical Observatory (DAO) data set we extend our bright limit to V = 10.8 and, for some stars, extend our time baseline to 35 years. Our magnitude limits include solar-mass main-sequence stars, subgiants, giants, and blue stragglers (BSs), and our spatial coverage extends radially to 17 pc (~13 core radii). For the WIYN data we find a measurement precision of 0.4 km/s for narrow-lined stars. We have measured RVs for 1046 stars in the direction of NGC 188, finding 473 to be likely cluster members. We detect 124 velocity-variable cluster members, all of which are likely to be dynamically hard-binary stars. Using our single member stars, we find an average cluster RV of -42.36 +/- 0.04 km/s. We use our precise RV and proper-motion membership data to greatly reduce field-star contamination in our cleaned color-magnitude diagram, from which we identify six stars of note that lie far from a standard single-star isochrone. We find the binaries to be centrally concentrated, providing evidence for the presence of mass segregation in NGC 188. We observe the BSs to populate a bimodal spatial distribution that is not centrally concentrated, suggesting that we may be observing two populations of BSs in NGC 188, including a centrally concentrated distribution as well as a halo population. Finally, we find NGC 188 to have a global RV dispersion of 0.64 +/- 0.04 km/s. When corrected for unresolved binaries, the NGC 188 RV dispersion has a nearly isothermal radial distribution. We use this mean-corrected velocity dispersion to derive a virial mass of 2300 +/- 460 MSun.

We present 98 spectroscopic binary orbits resulting from our ongoing radial-velocity survey of the old (7 Gyr) open cluster NGC 188. All but 13 are high-probability cluster members based on both radial-velocity and proper-motion membership analyses. Fifteen of these member binaries are double lined. Our stellar sample spans a magnitude range of 10.8 <=V<= 16.5 (1.14-0.92 Msun) and extends spatially to 17 pc (~13 core radii). All of our binary orbits have periods ranging from a few days to on the order of 10^3 days, and thus are hard binaries that dynamically power the cluster. For each binary, we present the orbital solutions and place constraints on the component masses. Additionally, we discuss a few binaries of note from our sample, identifying a likely blue straggler-blue straggler binary system (7782), a double-lined binary with a secondary star which is under-luminous for its mass (5080), two potential eclipsing binaries (4705 and 5762), and two binaries which are likely members of a quadruple system (5015a and 5015b).

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%.

The formation of Black Hole (BH) Low-Mass X-ray Binaries (LMXB) poses a theoretical challenge, as low-mass companions are not expected to survive the common-envelope scenario with the BH progenitor. Here we propose a formation mechanism that skips the common-envelope scenario and relies on triple-body dynamics. We study the evolution of hierarchical triples, following the secular dynamical evolution up to the octupole-level of approximation, including general relativity, tidal effects and post-main-sequence evolution, such as mass loss, changes to stellar radii and supernovae. During the dynamical evolution of the triple system, the "eccentric Kozai-Lidov" mechanism can cause large eccentricity excitations in the LMXB progenitor, resulting in three main BH-LMXB formation channels. Here we define BH-LMXB candidates as systems where the inner BH companion star crosses its Roche limit. In the "eccentric" channel (~ 81% of the LMXBs in our simulations), the donor star crosses its Roche limit during an extreme eccentricity excitation, while still on a wide orbit. Second, we find a "giant" LMXB channel (~ 11%), where a system undergoes only moderate eccentricity excitations, but the donor star fills its Roche lobe after evolving toward the giant branch. Third, we identify a "classical" channel (~8%), where tidal forces and magnetic braking shrink and circularize the orbit to short periods, triggering mass transfer. Finally, for the giant channel, we predict an eccentric ($\sim 0.3-0.6$), preferably inclined (~40, ~140 degreed) tertiary, typically on a wide enough orbit (~10^4AU), to potentially become unbound later in the triple evolution.

Characterization of the binary fractions in star clusters is of fundamental importance for many fields in astrophysics. Observations indicate that the majority of stars are found in binary systems, while most stars with masses greater than $0.5 M_\odot$ are formed in star clusters. In addition, since binaries are on average more massive than single stars, in resolved star clusters these systems are thought to be good tracers of (dynamical) mass segregation. Over time, dynamical evolution through two-body relaxation will cause the most massive objects to migrate to the cluster center, while the relatively lower-mass objects remain in or migrate to orbits at greater radii. This process will globally dominate a cluster's stellar distribution. However, close encounters involving binary systems may disrupt `soft' binaries. This process will occur more frequently in a cluster's central, dense region than in its periphery, which may mask the effects of mass segregation. Using high resolution Hubble Space Telescope observations, combined with sophisticated $N$-body simulations, we investigate the radial distributions of the main-sequence binary fractions in massive young Large Magellanic Cloud star clusters. We show that binary disruption may play an important role on very short timescales, depending on the environmental conditions in the cluster cores. This may lead to radial binary fractions that initially decline in the cluster centers, which is contrary to the effects expected from dynamical mass segregation.

(Abridged) We present results from 13776 radial-velocity (RV) measurements of 1278 candidate members of the old (4 Gyr) open cluster M67 (NGC 2682). The measurements are the results of a long-term survey that includes data from seven telescopes with observations for some stars spanning over 40 years. For narrow-lined stars, RVs are measured with precisions ranging from about 0.1 to 0.8 km/s. The combined stellar sample reaches from the brightest giants in the cluster down to about 4 magnitudes below the main-sequence turnoff (V = 16.5), covering a mass range of about 1.34 MSun to 0.76 MSun. Spatially, the sample extends to a radius of 30 arcmin (7.4 pc in projection at a distant of 850 pc or 6-7 core radii). We find M67 to have a mean RV of +33.64 km/s (with an internal precision of +/- 0.03 km/s). For stars with >=3 measurements, we derive RV membership probabilities and identify RV variables, finding 562 cluster members, 142 of which show significant RV variability. We use these cluster members to construct a color-magnitude diagram and identify a rich sample of stars that lie far from the standard single star isochrone, including the well-known blue stragglers, sub-subgiants and yellow giants. These exotic stars have a binary frequency of (at least) 80%, more than three times that detected for stars in the remainder of the sample. We confirm that the cluster is mass segregated, finding the binaries to be more centrally concentrated than the single stars in our sample at the 99.8% confidence level. The blue stragglers are centrally concentrated as compared to the solar-type main-sequence single stars in the cluster at the 99.7% confidence level. Accounting for both measurement precision and undetected binaries, we derive a RV dispersion in M67 of 0.59 +0.07 -0.06 km/s, which yields a virial mass for the cluster of 2100 +610 -550 MSun. WIYN Open Cluster Study. LXVII.

Strong encounters between single stars and binaries play a pivotal role in the evolution of star clusters. Such encounters can also dramatically modify the orbital parameters of binaries, exchange partners in and out of binaries, and are a primary contributor to the rate of physical stellar collisions in star clusters. Often, these encounters are studied under the approximation that they happen quickly enough and within a small enough volume to be considered isolated from the rest of the cluster. In this paper, we study the validity of this assumption through the analysis of a large grid of single - binary and binary - binary scattering experiments. For each encounter we evaluate the encounter duration, and compare this with the expected time until another single or binary star will join the encounter. We find that for lower-mass clusters, similar to typical open clusters in our Galaxy, the percent of encounters that will be "interrupted" by an interloping star or binary may be 20-40% (or higher) in the core, though for typical globular clusters we expect <1% of encounters to be interrupted. Thus, the assumption that strong encounters occur in relative isolation breaks down for certain clusters. Instead, many strong encounters develop into more complex "mini-clusters", which must be accounted for in studying, for example, the internal dynamics of star clusters, and the physical stellar collision rate.

We present the Robert C. Byrd Green Bank Telescope discovery of the highly eccentric binary millisecond pulsar PSR J1835$-$3259A in the Fermi Large Area Telescope-detected globular cluster NGC 6652. Timing over one orbit yields the pulse period 3.89 ms, orbital period 9.25 d, eccentricity $\sim 0.95$, and an unusually high companion mass of $0.74\,M_{\odot}$ assuming a $1.4\,M_{\odot}$ pulsar. We caution that the lack of data near periastron prevents a precise measurement of the eccentricity, and that further timing is necessary to constrain this and the other orbital parameters. From tidal considerations, we find that the companion must be a compact object. This system likely formed through an exchange encounter in the dense cluster environment. Our initial timing results predict the measurements of at least two post-Keplerian parameters with long-term phase-connected timing: the rate of periastron advance $\dot{\omega} \sim 0.1^{\circ}\,$yr$^{-1}$, requiring 1 yr of phase connection; and the Einstein delay $\gamma_{\mathrm{GR}} \sim 10\,$ms, requiring 2-3 yr of timing. For an orbital inclination $i > 50^{\circ}$, a measurement of $\sin{i}$ is also likely. PSR J1835$-$3259A thus provides an opportunity to measure the neutron star mass with high precision; to probe the cluster environment; and, depending on the nature of the companion, to investigate the limits of general relativity.

We develop a formalism for calculating probabilities for the outcomes of stellar dynamical interactions, based on results from $N$-body scattering experiments. We focus here on encounters involving up to six particles and calculate probabilities for direct stellar collisions; however our method is in principle valid for larger particle numbers. Our method relies on the binomial theorem, and is applicable to encounters involving any combination of particle radii. We further demonstrate that our base model is valid to within a few percent for any combination of particle masses, provided the minimum mass ratio is within a factor of a few from unity. This method is particularly suitable for models of collisional systems involving large numbers of stars, such as globular clusters, old open clusters and galactic nuclei, where small subsets of stars may regularly have very close encounters, and the direct integration of all such encounters is computationally expensive. Variations of our method may also be used to treat other encounter outcomes, such as ejections and exchanges.

The two Large Magellanic Cloud star clusters NGC 1805 and NGC 1818 are approximately the same chronological age (~30 Myr), but show different radial trends in binary frequency. The F-type stars (1.3 - 2.2 MSun) in NGC 1818 have a binary frequency that decreases towards the core, while the binary frequency for stars of similar mass in NGC 1805 is flat with radius, or perhaps bimodal (with a peak in the core). We show here, through detailed N-body modeling, that both clusters could have formed with the same primordial binary frequency and with binary orbital elements and masses drawn from the same distributions (defined from observations of open clusters and the field of our Galaxy). The observed radial trends in binary frequency for both clusters are best matched with models that have initial substructure. Furthermore, both clusters may be evolving along a very similar dynamical sequence, with the key difference that NGC 1805 is dynamically older than NGC 1818. The F-type binaries in NGC 1818 still show evidence of an initial period of rapid dynamical disruptions (which occur preferentially in the core), while NGC 1805 has already begun to recover a higher core binary frequency, owing to mass segregation (which will eventually produce a distribution in binary frequency that rises only towards the core, as is observed in old Milky Way star clusters). This recovery rate increases for higher-mass binaries, and therefore even at one age in one cluster, we predict a similar dynamical sequence in the radial distribution of the binary frequency as a function of binary primary mass.

This paper provides a detailed comparison of the differences in parameters derived for a star cluster from its color-magnitude diagrams depending on the filters and models used. We examine the consistency and reliability of fitting three widely-used stellar evolution models to fifteen combinations of optical and near-IR photometry for the old open cluster NGC 188. The optical filter response curves match those of the theoretical systems and are thus not the source of fit inconsistencies. NGC 188 is ideally suited to the present study thanks to a wide variety of high-quality photometry and available proper motions and radial velocities which enable us to remove non-cluster members and many binaries. Our Bayesian fitting technique yields inferred values of age, metallicity, distance modulus, and absorption as a function of the photometric band combinations and stellar models. We show that the historically-favored three band combinations of UBV and VRI can be meaningfully inconsistent with each other and with longer baseline datasets such as UBVRIJHKs. Differences among model sets can also be substantial. For instance, fitting Yi et al. (2001) and Dotter et al. (2008) models to UBVRIJHKs photometry for NGC 188 yields the following cluster parameters: age={5.78+ 0.03, 6.45+-0.04} Gyr, [Fe/H]={+0.125+-0.003, -0.077+-0.003} dex, m-M={11.441+-0.007, 11.525+-0.005} mag, and Av={0.162+-0.003, 0.236+-0.003} mag, respectively. Within the formal fitting errors, these two fits are substantially and statistically different. Such differences amongst fits using different filters and models are a cautionary tale regarding our current ability to fit star cluster color-magnitude diagrams. Additional modeling of this kind, with more models and star clusters, and future GAIA parallaxes are critical for isolating and quantifying the most relevant uncertainties in stellar evolutionary models.

We present the current state of the WOCS radial-velocity (RV) survey for the rich open cluster NGC 6819 (2.5 Gyr) including 93 spectroscopic binary orbits with periods ranging from 1.5 to 8,000 days. These results are the product of our ongoing RV survey of NGC 6819 using the Hydra Multi-Object Spectrograph on the WIYN 3.5 m telescope. We also include a detailed analysis of multiple prior sets of optical photometry for NGC 6819. Within a 1-degree field of view, our stellar sample includes the giant branch, the red clump, and blue straggler candidates, and extends to almost 2 mag below the main-sequence (MS) turnoff. For each star observed in our survey we present all RV measurements, the average RV and velocity variability information. Additionally, we discuss notable binaries from our sample, including eclipsing binaries (WOCS 23009, WOCS 24009, and WOCS 40007), stars noted in Kepler asteroseismology studies (WOCS 4008, WOCS 7009, and WOCS 8007), and potential descendants of past blue stragglers (WOCS 1006 and WOCS 6002). We find the incompleteness-corrected binary fraction for all MS binaries with periods less than 10,000 days to be 22% +/- 3% and a tidal circularization period of 6.2 +/- 1.1 days for NGC 6819.

In this paper, we address the question: What is the probability of stellar-mass black hole (BH) binaries co-existing in a globular cluster with an intermediate-mass black hole (IMBH)? Our results suggest that the detection of one or more BH binaries can strongly constrain the presence of an IMBH in most Galactic globular clusters. More specifically, the detection of one or more BH binaries could strongly indicate against the presence of an IMBH more massive than $\gtrsim 10^3$ M$_{\rm \odot}$ in roughly 80\% of the clusters in our sample. To illustrate this, we use a combination of N-body simulations and analytic methods to weigh the rate of formation of BH binaries against their ejection and/or disruption rate via strong gravitational interactions with the central (most) massive BH. The eventual fate of a sub-population of stellar-mass BHs (with or without binary companions) is for all BHs to be ejected from the cluster by the central IMBH, leaving only the most massive stellar-mass BH behind to form a close binary with the IMBH. During each phase of evolution, we discuss the rate of inspiral of the central BH-BH pair as a function of both the properties of the binary and its host cluster.

The old (7 Gyr) open cluster NGC 188 has yielded a wealth of astrophysical insight into its rich blue straggler population. Specifically, the NGC 188 blue stragglers are characterized by: A binary frequency of 80% for orbital periods less than $10^4$ days;Typical orbital periods around 1000 days;Typical secondary star masses of 0.5 M$_{\odot}$; At least some white dwarf companion stars; Modestly rapid rotation; A bimodal radial spatial distribution; Dynamical masses greater than standard stellar evolution masses (based on short-period binaries); Under-luminosity for dynamical masses (short-period binaries). Extensive $N$-body modeling of NGC 188 with empirical initial conditions reproduces the properties of the cluster, and in particular the main-sequence solar-type binary population. The current models also reproduce well the binary orbital properties of the blue stragglers, but fall well short of producing the observed number of blue stragglers. This deficit could be resolved by reducing the frequency of common-envelope evolution during Roche lobe overflow. Both the observations and the $N$-body models strongly indicate that the long-period blue-straggler binaries - which dominate the NGC 188 blue straggler population - are formed by asymptotic-giant (primarily) and red-giant mass transfer onto main sequence stars. The models suggest that the few non-velocity-variable blue stragglers formed from mergers or collisions. Several remarkable short-period double-lined binaries point to the importance of subsequent dynamical exchange encounters, and provide at least one example of a likely collisional origin for a blue straggler.

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.

The massive (13,000-26,000 Msolar), young (15-30 Myr) Large Magellanic Cloud star cluster NGC 1818 reveals an unexpected increasing binary frequency with radius for F-type stars (1.3-2.2 Msolar). This is in contrast to many older star clusters that show a decreasing binary frequency with radius. We study this phenomenon with sophisticated N-body modeling, exploring a range of initial conditions, from smooth virialized density distributions to highly substructured and collapsing configurations. We find that many of these models can reproduce the cluster's observed properties, although with a modest preference for substructured initial conditions. Our models produce the observed radial trend in binary frequency through disruption of soft binaries (with semi-major axes, a > 3000 AU), on approximately a crossing time (~ 5.4 Myr), preferentially in the cluster core. Mass segregation subsequently causes the binaries to sink towards the core. After roughly one initial half-mass relaxation time (t_rh(0) ~ 340 Myr) the radial binary frequency distribution becomes bimodal, the innermost binaries having already segregated towards the core, leaving a minimum in the radial binary frequency distribution that marches outwards with time. After 4-6 t_rh(0), the rising distribution in the halo disappears, leaving a radial distribution that rises only towards the core. Thus, both a radial binary frequency distribution that falls towards the core (as observed for NGC 1818) and one that rises towards the core (as for older star clusters) can arise naturally from the same evolutionary sequence owing to binary disruption and mass segregation in rich star clusters.

We analyze extensive BVR_cI_c time-series photometry and radial-velocity measurements for WOCS 40007 (Auner 259; KIC 5113053), a double-lined detached eclipsing binary and a member of the open cluster NGC 6819. Utilizing photometric observations from the 1-meter telescope at Mount Laguna Observatory and spectra from the WIYN 3.5-meter telescope, we measure precise and accurate masses (~1.6% uncertainty) and radii (~0.5%) for the binary components. In addition, we discover a third star orbiting the binary with a period greater than 3000 days using radial velocities and Kepler eclipse timings. Because the stars in the eclipsing binary are near the cluster turnoff, they are evolving rapidly in size and are sensitive to age. With a metallicity of [Fe/H]=+0.09+/-0.03, we find the age of NGC 6819 to be about 2.4 Gyr from CMD isochrone fitting and 3.1+/-0.4 Gyr by analyzing the mass-radius (M-R) data for this binary. The M-R age is above previous determinations for this cluster, but consistent within 1 sigma uncertainties. When the M-R data for the primary star of the additional cluster binary WOCS 23009 is included, the weighted age estimate drops to 2.5+/-0.2 Gyr, with a systematic uncertainty of at least 0.2 Gyr. The age difference between our CMD and M-R findings may be the result of systematic error in the metallicity or helium abundance used in models, or due to slight radius inflation of one or both stars in the WOCS 40007 binary.

Over the last few decades, observational surveys have revealed that high-order multiple-star systems (e.g. triples, quadruples, etc.), and triples in particular, are common in our Galaxy. In this paper, we consider the dynamical significance of this transformation in our understanding of stellar multiplicity. Using empirically constrained binary and triple fractions in those star clusters for which these values are available in the literature, we compare analytic rates for encounters involving single, binary, and triple stars. Our results show that, \textit{even for relatively low triple fractions, dynamical interactions involving triples occur roughly as often as encounters involving either single or binary stars alone, particularly in low-mass star clusters}. More generally, using empirically-derived multiple star catalogues for the young star-forming association Taurus-Auriga and the Galactic field, we show that the data are consistent with the gravitationally-focused cross section for encounters increasing with increasing multiplicity. Consequently, triple stars, and even higher-order multiples, could be more important than previously realized for a number of astrophysical phenomena, including the formation and destruction of compact binaries and various types of stellar exotica, and the dynamical evolution of star clusters.

We present the discovery of the totally eclipsing long-period (P = 771.8 d) binary system WOCS 23009 in the old open cluster NGC 6819 that contains both an evolved star near central hydrogen exhaustion and a low-mass (0.45 Msun) star. This system was previously known to be a single-lined spectroscopic binary, but the discovery of an eclipse near apastron using data from the Kepler space telescope makes it clear that the system has an inclination that is very close to 90 degrees. Although the secondary star has not been identified in spectra, the mass of the primary star can be constrained using other eclipsing binaries in the cluster. The combination of total eclipses and a mass constraint for the primary star allows us to determine a reliable mass for the secondary star and radii for both stars, and to constrain the cluster age. Unlike well-measured stars of similar mass in field binaries, the low-mass secondary is not significantly inflated in radius compared to model predictions. The primary star characteristics, in combination with cluster photometry and masses from other cluster binaries, indicates a best age of 2.62+/-0.25 Gyr, although stellar model physics may introduce systematic uncertainties at the ~10% level. We find preliminary evidence that the asteroseismic predictions for red giant masses in this cluster are systematically too high by as much as 8%.

(Abridged) Following on from a recently completed radial-velocity survey of the old (7 Gyr) open cluster NGC 188 in which we study in detail the solar-type hard binaries and blue stragglers of the cluster, here we investigate the dynamical evolution of NGC 188 through a sophisticated N-body model. We employ the observed binary properties of the young (150 Myr) open cluster M35, where possible, to guide our choices for parameters of the initial binary population. At 7 Gyr the main-sequence solar-type hard-binary population in the model matches that of NGC 188 in both binary frequency and distributions of orbital parameters. This agreement between the model and observations is in a large part due to the similarities between the NGC 188 and M35 solar-type binaries. Indeed, among the 7 Gyr main-sequence binaries in the model, only those with P>1000 days show potentially observable evidence for modifications by dynamical encounters. This emphasizes the importance of defining accurate initial conditions for star cluster models, which we propose is best accomplished through comparisons with observations of young open clusters like M35. Furthermore, this suggests that observations of the present-day binaries in even old open clusters can provide valuable information on their primordial binary populations. However, despite the model's successes at matching the true cluster, the model underproduces blue stragglers and produces an overabundance of long-period circular main-sequence--white-dwarf binaries as compared to NGC 188. We conclude that improvements in the physics of mass transfer and common envelope may in fact solve both discrepancies with the observations. This project highlights the unique accessibility of open clusters to both comprehensive observational surveys and full-scale N-body simulations, and underscores the importance of open clusters to the study of star cluster dynamics.

In this paper, we study small-N gravitational dynamics involving up to six objects. We perform a large suite of numerical scattering experiments involving single, binary, and triple stars. This is done using the FEWBODY numerical scattering code, which we have upgraded to treat encounters involving triple stars. We focus on outcomes that result in direct physical collisions between stars, within the low angular momentum and high absolute orbital energy regime. The dependence of the collision probability on the number of objects involved in the interaction, N, is found for fixed total energy and angular momentum. Our results are consistent with a collision probability that increases approximately as N^2. Interestingly, this is also what is expected from the mean free path approximation in the limit of very large N. A more thorough exploration of parameter space will be required in future studies to fully explore this potentially intriguing connection. This study is meant as a first step in an on-going effort to extend our understanding of small-N collisional dynamics beyond the three- and four-body problems and into the realm of larger-N.

While hundreds of planets have been discovered around field stars, only a few are known in star clusters. To explain the lack of short-period giant planets in globular clusters (GC), such as 47 Tucane and \omega\ Centauri, it has been suggested that their low metallicities may have prevented planet formation. Alternatively, the high rates of close stellar encounters in these clusters may have influenced the formation and subsequent evolution of planetary systems. How common are planets in clusters around normal main-sequence stars? Here we consider whether this question can be addressed using data from the Kepler mission. The Kepler field of view contains 4 low-density (relative to GCs) open clusters where the metallicities are about solar (or even higher) and stellar encounters are much less frequent than in typical GCs. We provide detailed $N$-body models and show that most planets in Kepler-detectable orbits are not significantly perturbed by stellar encounters in these open clusters. We focus on the most massive cluster, NGC 6791, which has super-solar metallicity, and find that if planets formed in this cluster at the same frequency as observed in the field, Kepler could detect 1 -- 20 transiting planets depending on the planet-size distribution and the duration of data collection. Due to the large distance to NGC 6791 Kepler will have to search relatively faint ($K_p<20$) stars for the full extended mission to achieve such a yield.

(abridged) We present an in-depth study of the hard-binary population of the old (7 Gyr) open cluster NGC 188. The main-sequence solar-type hard binaries in NGC 188 are nearly indistinguishable from similar binaries in the Galactic field. We find a global solar-type main-sequence hard-binary frequency in NGC 188 of 29 +/- 3 % for binaries with periods less than 10^4 days. For main-sequence hard binaries in the cluster we observe a log-period distribution that rises towards our detection limit, a roughly Gaussian eccentricity distribution centered on e = 0.35 (for binaries with periods longer than the circularization period), and a secondary-mass distribution that rises towards lower-mass companions. Importantly, the NGC 188 blue straggler binaries show significantly different characteristics than the solar-type main sequence binaries in NGC 188. We observe a blue straggler hard-binary frequency of 76 +/- 19 %, three times that of the main sequence. The blue straggler binary eccentricity - log period distribution is distinct from that of the main sequence at the 99% confidence level, with the majority of the blue straggler binaries having periods of order 1000 days and lower eccentricities. The secondary-mass distribution for these long-period blue straggler binaries is narrow and peaked with a mean value of about 0.5 Msun. Predictions for mass-transfer products are most closely consistent with the binary properties of these NGC 188 blue stragglers, which comprise two-thirds of the blue straggler population. Additionally we compare the NGC 188 binaries to those evolved within the sophisticated Hurley et al. (2005) N-body open cluster simulation. We find that additional simulations with initial conditions that are better motivated by observations are necessary to properly investigate the dynamical evolution of a rich binary population in open clusters like NGC 188.

We present the first study of the X-ray population of the intermediate-age rich open cluster NGC 6819 using the XMM-Newton Observatory. In the past decade, Chandra X-ray observations have shown a relationship between the X-ray population of globular clusters and their internal dynamics and encounter frequency. We investigate the role dynamics possibly play in the formation of X-ray sources in NGC 6819, and compare our results with known properties of field and globular cluster X-ray populations. We implement a multi wavelength approach to studying the X-ray sources, utilizing X-ray and UV data from XMM observations along with the wealth of photometry and radial-velocity data from the WIYN Open Cluster Study (WOCS) and the CFHT Open Cluster Survey. Within the cluster half-light radius we detect 12 X-ray sources down to a luminosity of 10^30 erg/s for cluster members. The sources include a candidate quiescent low-mass X-ray binary (qLMXB), a candidate cataclysmic variable, and two active binary systems. The presence of a qLMXB in an open cluster is previously unexpected given the known relationships between luminous X-ray sources and encounter frequency in globular clusters, and most likely has a dynamical origin.

In open star clusters, where all members formed at about the same time, blue straggler stars are typically observed to be brighter and bluer than hydrogen-burning main-sequence stars, and therefore should already have evolved into giant stars and stellar remnants. Correlations between blue straggler frequency and cluster binary fraction, core mass, and radial position suggest that mass transfer or mergers in binary stars dominates the production of blue stragglers in open clusters. Analytic models, detailed observations, and sophisticated N-body simulations, however, argue in favor of stellar collisions. Here we report that the blue stragglers in long-period binaries in the old (7 Gyr) open cluster NGC 188 have companions with masses of about half a solar mass, with a surprisingly narrow mass distribution. This conclusively rules out a collisional origin, as the collision hypothesis predicts a companion-mass distribution with significantly higher masses. Mergers in hierarchical triple stars are marginally permitted by the data, but the observations do not favor this hypothesis. The data are closely consistent with a mass transfer origin for the long-period blue straggler binaries in NGC 188, in which the companions would be white dwarfs of about half a solar mass.

We present 5201 radial-velocity measurements of 1144 stars, as part of an ongoing study of the young (150 Myr) open cluster M35 (NGC 2168). We have observed M35 since 1997, using the Hydra Multi-Object Spectrograph on the WIYN 3.5m telescope. Our stellar sample covers main-sequence stars over a magnitude range of 13.0<V<16.5 (1.6 - 0.8 Msun) and extends spatially to a radius of 30 arcminutes (7 pc in projection at a distance of 805 pc or 4 core radii). Due to its youth, M35 provides a sample of late-type stars with a range of rotation periods. Therefore, we analyze the radial-velocity measurement precision as a function of the projected rotational velocity. For narrow-lined stars (v sin i < 10 km/s), the radial velocities have a precision of 0.5 km/s, which degrades to 1.0 km/s for stars with v sin i = 50 km/s. The radial-velocity distribution shows a well-defined cluster peak with a central velocity of -8.16 +/- 0.05 km/s, permitting a clean separation of the cluster and field stars. For stars with >=3 measurements, we derive radial-velocity membership probabilities and identify radial-velocity variables, finding 360 cluster members, 55 of which show significant radial- velocity variability. Using these cluster members, we construct a color-magnitude diagram for our stellar sample cleaned of field star contamination. We also compare the spatial distribution of the single and binary cluster members, finding no evidence for mass segregation in our stellar sample. Accounting for measurement precision, we place an upper limit on the radial-velocity dispersion of the cluster of 0.81 +/- 0.08 km/s. After correcting for undetected binaries, we derive a true radial-velocity dispersion of 0.65 +/- 0.10 km/s.

We present a direct N-body simulation modeling the evolution of the old (7 Gyr) open cluster NGC 188. This is the first N-body open cluster simulation whose initial binary population is directly defined by observations of a specific open cluster: M35 (150 Myr). We compare the simulated color-magnitude diagram at 7 Gyr to that of NGC 188, and discuss the blue stragglers produced in the simulation. We compare the solar-type main sequence binary period and eccentricity distributions of the simulation to detailed observations of similar binaries in NGC 188. These results demonstrate the importance of detailed observations in guiding N-body open cluster simulations. Finally, we discuss the implications of a few discrepancies between the NGC 188 model and observations and suggest a few methods for bringing N-body open cluster simulations into better agreement with observations.

We present the current results from our ongoing radial-velocity survey of the intermediate-age (2.4 Gyr) open cluster NGC 6819. Using both newly observed and other available photometry and astrometry we define a primary target sample of 1454 stars that includes main-sequence, subgiant, giant, and blue straggler stars, spanning a magnitude range of 11<=V<=16.5 and an approximate mass range of 1.1 to 1.6 Msun. Our sample covers a 23 arcminute (13 pc) square field of view centered on the cluster. We have measured 6571 radial velocities for an unbiased sample of 1207 stars in the direction of the open cluster NGC 6819, with a single-measurement precision of 0.4 km/s for most narrow-lined stars. We use our radial-velocity data to calculate membership probabilities for stars with >= 3 measurements, providing the first comprehensive membership study of the cluster core that includes stars from the giant branch through the upper main sequence. We identify 480 cluster members. Additionally, we identify velocity-variable systems, all of which are likely hard binaries that dynamically power the cluster. Using our single cluster members, we find a cluster average radial velocity of 2.34 +/- 0.05 km/s . We use our kinematic cluster members to construct a cleaned color- magnitude diagram from which we identify rich giant, subgiant, and blue straggler populations and a well-defined red clump. The cluster displays a morphology near the cluster turnoff clearly indicative of core convective overshoot. Finally, we discuss a few stars of note, one of which is a short-period red-clump binary that we suggest may be the product of a dynamical encounter.

The mass and chemical composition of a star are the primary determinants of its basic physical properties--radius, temperature, luminosity--and how those properties evolve with time. Thus, two stars born at the same time, from the same natal material, and with the same mass are 'identical twins,' and as such might be expected to possess identical physical attributes. We have discovered in the Orion Nebula a pair of stellar twins in a newborn binary star system. Each star in the binary has a mass of 0.41 +/- 0.01 solar masses, identical to within 2 percent. Here we report that these twin stars have surface temperatures that differ by ~300K (~10%), and luminosities that differ by ~50%, both at high confidence level. Preliminary results indicate that the stars' radii also differ, by 5-10%. These surprising dissimilarities suggest that one of the twins may have been delayed by several hundred thousand years in its formation relative to its sibling. Such a delay could only have been detected in a very young, definitively equal-mass binary system3 such as that reported here. Our findings reveal cosmic limits on the age synchronisation of young binary stars, often used as tests for the age calibrations of star-formation models.