Showing 1–50 of 124 results for author: Sills, A

Massive interacting binaries as an enrichment source for multiple populations in star clusters

Authors: Michelle Nguyen, Alison Sills

Abstract: We present a suite of binary evolution models with massive primaries (10 $\leq$ M$_1$ $\leq$ 40 M$_\odot$) and periods and mass ratios chosen such that the systems undergo non-conservative mass transfer while the primaries have helium cores. We track the total mass and chemical composition of the ejecta from these systems. This material shows the abundance signatures of hot hydrogen burning which… ▽ More We present a suite of binary evolution models with massive primaries (10 $\leq$ M$_1$ $\leq$ 40 M$_\odot$) and periods and mass ratios chosen such that the systems undergo non-conservative mass transfer while the primaries have helium cores. We track the total mass and chemical composition of the ejecta from these systems. This material shows the abundance signatures of hot hydrogen burning which are needed to explain the abundance patterns seen in multiple populations in massive star clusters. We then calculate the total yield of a population of binary stars with masses, mass ratios, and periods consistent with their distribution in a field population. We show that the overall abundance of this material is enriched in helium, nitrogen, sodium, and aluminum, and depleted in carbon, oxygen, and magnesium, by amounts that are consistent with observations. We also show that such a population of binaries will return approximately 25% of its mass in this ejecta (compared to 4% if all the stars were single), over a characteristic timescale of about 12 Myr. We argue that massive binaries must be seriously considered as a contributor to the source of enriched material needed to explain the multiple populations in massive clusters, since essentially all massive stars are formed in binaries or higher order multiples, massive binaries are primarily formed in clusters, and massive binaries naturally produce material of the right composition. △ Less

Submitted 9 May, 2024; originally announced May 2024.

Comments: Accepted by The Astrophysical Journal

Dynamics of Star Cluster Formation: Mergers in Gas Rich Environments

Authors: Jeremy Karam, Alison Sills

Abstract: We perform high resolution simulations of forming star clusters as they merge inside giant molecular clouds (GMCs) using hydrodynamics coupled to N-body dynamics to simultaneously model both the gas and stars. We zoom in to previously run GMC simulations and resolve clusters into their stellar and gas components while including the surrounding GMC environment. We find that GMC gas is important in… ▽ More We perform high resolution simulations of forming star clusters as they merge inside giant molecular clouds (GMCs) using hydrodynamics coupled to N-body dynamics to simultaneously model both the gas and stars. We zoom in to previously run GMC simulations and resolve clusters into their stellar and gas components while including the surrounding GMC environment. We find that GMC gas is important in facilitating the growth of clusters in their embedded phase by promoting cluster mergers. Mergers induce asymmetric expansion of the stellar component of the clusters in our simulations. As well, mergers induce angular momentum in the clusters' stellar and gas components. We find that mergers can lead to an increase in the amount of dense gas present in clusters if a background gas distribution is present. We predict that this can lead to new star formation that can change the overall distribution of cluster stars in velocity space. Our results suggest that subcluster mergers in the presence of background gas can imprint dynamical signatures that can be used to constrain cluster formation histories. △ Less

Submitted 9 April, 2024; originally announced April 2024.

Comments: Accepted for publication in The Astrophysical Journal, 15 pages, 10 figures

Galactic properties that favour star cluster formation: a statistical view

Authors: Samantha C. Berek, Marta Reina-Campos, Gwendolyn Eadie, Alison Sills

Abstract: The presence or absence of star clusters in galaxies, and the properties of star cluster populations compared to their host galaxy properties, are important observables for validating models of cluster formation, galaxy formation, and galaxy assembly. In this work, we apply a Bayesian approach to fit two models to data from surveys of young clusters in star forming galaxies. The first model is a l… ▽ More The presence or absence of star clusters in galaxies, and the properties of star cluster populations compared to their host galaxy properties, are important observables for validating models of cluster formation, galaxy formation, and galaxy assembly. In this work, we apply a Bayesian approach to fit two models to data from surveys of young clusters in star forming galaxies. The first model is a logistic regression, which allows us to include galaxies which do not have any young clusters. The second model is a hurdle model, which includes galaxies with zero clusters and also incorporates information about the total mass in the cluster system. We investigate two predictors (star formation rate and total stellar mass in the galaxy) and look at clusters younger than 10 or 100 Myr. We find that in all cases, star formation rate is the better predictor for both the probability of hosting clusters and the total mass in the cluster system. We compare our results to similar models for old globular clusters, and conclude that star cluster formation was more abundant and more efficient at higher redshifts, likely because of the high gas content of galaxies at that time. △ Less

Submitted 26 July, 2023; originally announced July 2023.

Comments: 10 pages, 5 figures, accepted by MNRAS

Initial sizes of star clusters: implications for cluster dissolution during galaxy evolution

Authors: Marta Reina-Campos, Alison Sills, Godefroy Bichon

Abstract: Massive star clusters are often used as tracers of galaxy formation and assembly. In order to do so, we must understand their properties at formation, and how those properties change with time, galactic environment, and galaxy assembly history. The two most important intrinsic properties that govern star cluster evolution are mass and radius. In this paper, we investigate 10 theoretically and obse… ▽ More Massive star clusters are often used as tracers of galaxy formation and assembly. In order to do so, we must understand their properties at formation, and how those properties change with time, galactic environment, and galaxy assembly history. The two most important intrinsic properties that govern star cluster evolution are mass and radius. In this paper, we investigate 10 theoretically and observationally motivated initial size-mass relations for star clusters, and evolve populations of clusters through galaxy formation models. We compare our results to each other and to observations of cluster populations in M83, M31, and the Milky Way. We find that none of our size-mass relations agree with the observations after 6-10 Gyr of evolution. We can successfully reproduce the cluster mass functions with models that have a small range of initial radii, and which do not allow cluster radii to change with time. However, these models do not agree with our understanding of cluster evolution, which does involve radius evolution, and do not match the observed distributions of radii. We note that there is a region of parameter space where clusters are optimally protected from both tidal shocks and evaporation due to two-body relaxation. Clusters which are allowed to evolve into this parameter space will likely survive. An improved understanding of both mass and radius evolution of star clusters in realistic, time-varying galactic potentials is necessary to appropriately make the connection between present-day cluster properties and their use as tracers of galaxy formation and assembly. △ Less

Submitted 30 June, 2023; originally announced June 2023.

Comments: Accepted by MNRAS

Fast rotating Blue Stragglers prefer loose clusters

Authors: Francesco R. Ferraro, Alessio Mucciarelli, Barbara Lanzoni, Cristina Pallanca, Mario Cadelano, Alex Billi, Alison Sills, Enrico Vesperini, Emanuele Dalessandro, Giacomo Beccari, Lorenzo Monaco, Mario Mateo

Abstract: Blue stragglers are anomalously luminous core hydrogen-burning stars formed through mass-transfer in binary/triple systems and stellar collisions. Their physical and evolutionary properties are largely unknown and unconstrained. Here we analyze 320 high-resolution spectra of blue stragglers collected in eight galactic globular clusters with different structural characteristics and show evidence th… ▽ More Blue stragglers are anomalously luminous core hydrogen-burning stars formed through mass-transfer in binary/triple systems and stellar collisions. Their physical and evolutionary properties are largely unknown and unconstrained. Here we analyze 320 high-resolution spectra of blue stragglers collected in eight galactic globular clusters with different structural characteristics and show evidence that the fraction of fast rotating blue stragglers (with rotational velocities larger than 40 km/s) increases for decreasing central density of the host system. This trend suggests that fast spinning blue stragglers prefer low-density environments and promises to open an unexplored route towards understanding the evolutionary processes of these stars. Since large rotation rates are expected in the early stages of both formation channels, our results provide direct evidence for recent blue straggler formation activity in low-density environments and put strong constraints on the timescale of the collisional blue straggler slow-down processes. △ Less

Submitted 15 May, 2023; originally announced May 2023.

Comments: Published in Nature Communications

Journal ref: Nature Communications 14, 2584 (2023)

Modelling Star Cluster Formation: Gas Accretion

Authors: Jeremy Karam, Alison Sills

Abstract: The formation of star clusters involves the growth of smaller, gas-rich subclusters through accretion of gas from the giant molecular cloud within which the subclusters are embedded. The two main accretion mechanisms responsible for this are accretion of gas from dense filaments, and from the ambient background of the cloud. We perform simulations of both of these accretion processes onto gas-rich… ▽ More The formation of star clusters involves the growth of smaller, gas-rich subclusters through accretion of gas from the giant molecular cloud within which the subclusters are embedded. The two main accretion mechanisms responsible for this are accretion of gas from dense filaments, and from the ambient background of the cloud. We perform simulations of both of these accretion processes onto gas-rich star clusters using coupled smoothed particle hydrodynamics to model the gas, and N-body dynamics to model the stars. We find that, for both accretion processes, the accreting star cluster loses some of its original mass while gaining mass from either the ambient background or the dense filament. The amount of mass lost from both these processes is small compared to the total mass of the cluster. However, in the case of accretion from a background medium, the net effect can be a decrease in the total mass of the cluster if it is travelling fast enough through the ambient medium ($> 4$kms$^{-1}$). We find that the amount of mass lost from the cluster through filamentary accretion is independent of the density, width, or number of filaments funneling gas into the cluster and is always such that the mass of the cluster is constantly increasing with time. We compare our results to idealized prescriptions used to model star cluster formation in larger scale GMC simulations and find that such prescriptions act as an upper limit when describing the mass of the star cluster they represent. △ Less

Submitted 20 March, 2023; originally announced March 2023.

Comments: Accepted for publication in MNRAS, 14 pages, 11 figures

Early Evolution and 3D Structure of Embedded Star Clusters

Authors: Claude Cournoyer-Cloutier, Alison Sills, William E. Harris, Sabrina M. Appel, Sean C. Lewis, Brooke Polak, Aaron Tran, Martijn J. C. Wilhelm, Mordecai-Mark Mac Low, Stephen L. W. McMillan, Simon Portegies Zwart

Abstract: We perform simulations of star cluster formation to investigate the morphological evolution of embedded star clusters in the earliest stages of their evolution. We conduct our simulations with Torch, which uses the AMUSE framework to couple state-of-the-art stellar dynamics to star formation, radiation, stellar winds, and hydrodynamics in FLASH. We simulate a suite of $10^4$ M$_{\odot}$ clouds at… ▽ More We perform simulations of star cluster formation to investigate the morphological evolution of embedded star clusters in the earliest stages of their evolution. We conduct our simulations with Torch, which uses the AMUSE framework to couple state-of-the-art stellar dynamics to star formation, radiation, stellar winds, and hydrodynamics in FLASH. We simulate a suite of $10^4$ M$_{\odot}$ clouds at 0.0683 pc resolution for $\sim$ 2 Myr after the onset of star formation, with virial parameters $α$ = 0.8, 2.0, 4.0 and different random samplings of the stellar initial mass function and prescriptions for primordial binaries. Our simulations result in a population of embedded clusters with realistic morphologies (sizes, densities, and ellipticities) that reproduce the known trend of clouds with higher initial $α$ having lower star formation efficiencies. Our key results are as follows: (1) Cluster mass growth is not monotonic, and clusters can lose up to half of their mass while they are embedded. (2) Cluster morphology is not correlated with cluster mass and changes over $\sim$ 0.01 Myr timescales. (3) The morphology of an embedded cluster is not indicative of its long-term evolution but only of its recent history: radius and ellipticity increase sharply when a cluster accretes stars. (4) The dynamical evolution of very young embedded clusters with masses $\lesssim$ 1000 M$_{\odot}$ is dominated by the overall gravitational potential of the star-forming region rather than by internal dynamical processes such as two- or few-body relaxation. △ Less

Submitted 16 February, 2023; originally announced February 2023.

Comments: 15 pages, 10 figures, to be published in MNRAS

Using Molecular Gas Observations to Guide Initial Conditions for Star Cluster Simulations

Authors: Alison Sills, Steven Rieder, Anne S. M. Buckner, Alvaro Hacar, Simon Portegies Zwart, Paula S. Teixeira

Abstract: The earliest evolution of star clusters involves a phase of co-existence of both newly-formed stars, and the gas from which they are forming. Observations of the gas in such regions provide a wealth of data that can inform the simulations which are needed to follow the evolution of such objects forward in time. We present a method for transforming the observed gas properties into initial condition… ▽ More The earliest evolution of star clusters involves a phase of co-existence of both newly-formed stars, and the gas from which they are forming. Observations of the gas in such regions provide a wealth of data that can inform the simulations which are needed to follow the evolution of such objects forward in time. We present a method for transforming the observed gas properties into initial conditions for simulations that include gas, stars, and ongoing star formation. We demonstrate our technique using the Orion Nebula Cluster. Since the observations cannot provide all the necessary information for our simulations, we make choices for the missing data and assess the impact of those choices. We find that the results are insensitive to the adopted choices of the gas velocity in the plane of the sky. The properties of the surrounding gas cloud (e.g. overall density and size), however, have an effect on the star formation rate and pace of assembly of the resultant star cluster. We also analyze the stellar properties of the cluster and find that the stars become more tightly clustered and in a stronger radial distribution even as new stars form in the filament. △ Less

Submitted 16 December, 2022; originally announced December 2022.

Comments: 11 pages, 12 figures, accepted for publication in MNRAS

Early-Forming Massive Stars Suppress Star Formation and Hierarchical Cluster Assembly

Authors: Sean C. Lewis, Stephen L. W. McMillan, Mordecai-Mark Mac Low, Claude Cournoyer-Cloutier, Brooke Polak, Martijn J. C. Wilhelm, Aaron Tran, Alison Sills, Simon Portegies Zwart, Ralf S. Klessen, Joshua E. Wall

Abstract: Feedback from massive stars plays an important role in the formation of star clusters. Whether a very massive star is born early or late in the cluster formation timeline has profound implications for the star cluster formation and assembly processes. We carry out a controlled experiment to characterize the effects of early-forming massive stars on star cluster formation. We use the star formation… ▽ More Feedback from massive stars plays an important role in the formation of star clusters. Whether a very massive star is born early or late in the cluster formation timeline has profound implications for the star cluster formation and assembly processes. We carry out a controlled experiment to characterize the effects of early-forming massive stars on star cluster formation. We use the star formation software suite \texttt{Torch}, combining self-gravitating magnetohydrodynamics, ray-tracing radiative transfer, $N$-body dynamics, and stellar feedback to model four initially identical $10^4$ M$_\odot$ giant molecular clouds with a Gaussian density profile peaking at $521.5 \mbox{ cm}^{-3}$. Using the \texttt{Torch} software suite through the \texttt{AMUSE} framework we modify three of the models to ensure that the first star that forms is very massive (50, 70, 100 M$_\odot$). Early-forming massive stars disrupt the natal gas structure, resulting in fast evacuation of the gas from the star forming region. The star formation rate is suppressed, reducing the total mass of stars formed. Our fiducial control model without an early massive star has a larger star formation rate and total efficiency by up to a factor of three and a higher average star formation efficiency per free-fall time by up to a factor of seven. Early-forming massive stars promote the buildup of spatially separate and gravitationally unbound subclusters, while the control model forms a single massive cluster. △ Less

Submitted 28 February, 2023; v1 submitted 2 December, 2022; originally announced December 2022.

Comments: 17 pages, 7 figures. Published in ApJ

Journal ref: ApJ 944 211 (2023)

Modelling Star Cluster Formation: Mergers

Authors: Jeremy Karam, Alison Sills

Abstract: Star cluster formation in giant molecular clouds involves the local collapse of the cloud into small gas-rich subclusters, which can then subsequently collide and merge to build up the final star cluster(s). In this paper, we simulate collisions between these subclusters, using coupled smooth particle hydrodynamics for the gas and N-body dynamics for the stars. We are guided by previous radiation… ▽ More Star cluster formation in giant molecular clouds involves the local collapse of the cloud into small gas-rich subclusters, which can then subsequently collide and merge to build up the final star cluster(s). In this paper, we simulate collisions between these subclusters, using coupled smooth particle hydrodynamics for the gas and N-body dynamics for the stars. We are guided by previous radiation hydrodynamics simulations of molecular cloud collapse which provide the global properties of the colliding clusters, such as their stellar and gas masses, and their initial positions and velocities. The subclusters in the original simulation were treated as sink particles which immediately merged into a single entity after the collision. We show that the more detailed treatment provides a more complex picture. At collisional velocities above ~ 10 km/s, the stellar components of the cluster do not form a monolithic cluster within 3 Myr, although the gas may do so. At lower velocities, the clusters do eventually merge but over timescales that may be longer than the time for a subsequent collision. The structure of the resultant cluster is not well-fit by any standard density distribution, and the clusters are not in equilibrium but continue to expand over our simulation time. We conclude that the simple sink particle treatment of subcluster mergers in large-scale giant molecular cloud simulations provides an upper limit on the final cluster properties. △ Less

Submitted 6 May, 2022; originally announced May 2022.

Comments: Accepted for publication in MNRAS, 11 pages, 12 figures

Constraining the shape of dark matter haloes with globular clusters

Authors: Marta Reina-Campos, Sebastian Trujillo-Gomez, Joel L. Pfeffer, Alison Sills, Alis J. Deason, Robert A. Crain, J. M. Diederik Kruijssen

Abstract: We explore how diffuse stellar light and globular clusters (GCs) can be used to trace the matter distribution of their host halo using an observational methodology. For this, we use 117 simulated dark matter (DM) haloes from the $(34.4~\rm cMpc)^3$ periodic volume of the E-MOSAICS project. For each halo, we compare the stellar surface brightness and GC projected number density maps to the surface… ▽ More We explore how diffuse stellar light and globular clusters (GCs) can be used to trace the matter distribution of their host halo using an observational methodology. For this, we use 117 simulated dark matter (DM) haloes from the $(34.4~\rm cMpc)^3$ periodic volume of the E-MOSAICS project. For each halo, we compare the stellar surface brightness and GC projected number density maps to the surface densities of DM and total mass. We find that the dominant structures identified in the stellar light and in the GCs correspond closely with those from the DM and total mass. Our method is unaffected by the presence of satellites and its precision improves with fainter GC samples. We recover tight relations between the profiles of stellar surface brightness and GC number density to those of the DM, suggesting that the profile of DM can be accurately recovered from the stars and GCs ($σ\leq0.5~$dex). We quantify the projected morphology of DM, stars and GCs, and find that the stars and GCs are more flattened than the DM. Additionally, the semi-major axes of the distribution of stars and GCs are typically misaligned by $\sim 10~$degrees from that of DM. We demonstrate that deep imaging of diffuse stellar light and GCs can place constraints on the shape, profile and orientation of their host halo. These results extend down to haloes with central galaxies $M_{\star}\geq10^{10}~M_{\odot}$, and the analysis will be applicable to future data from the Euclid, Roman and the Rubin observatories. △ Less

Submitted 25 April, 2022; originally announced April 2022.

Comments: 17 pages, 15 figures; submitted to MNRAS (25/04/2022) - comments and suggestions welcome!

Introducing EMP-Pathfinder: modelling the simultaneous formation and evolution of stellar clusters in their host galaxies

Authors: Marta Reina-Campos, Benjamin W. Keller, J. M. Diederik Kruijssen, Jindra Gensior, Sebastian Trujillo-Gomez, Sarah M. R. Jeffreson, Joel L. Pfeffer, Alison Sills

Abstract: The formation and evolution of stellar clusters is intimately linked to that of their host galaxies. To study this connection, we present the EMP-Pathfinder suite of cosmological zoom-in Milky Way-mass simulations. These simulations contain a sub-grid description for stellar cluster formation and evolution, allowing us to study the simultaneous formation and evolution of stellar clusters alongside… ▽ More The formation and evolution of stellar clusters is intimately linked to that of their host galaxies. To study this connection, we present the EMP-Pathfinder suite of cosmological zoom-in Milky Way-mass simulations. These simulations contain a sub-grid description for stellar cluster formation and evolution, allowing us to study the simultaneous formation and evolution of stellar clusters alongside their host galaxies across cosmic time. As a key ingredient in these simulations, we include the physics of the multi-phase nature of the interstellar medium (ISM), which enables studies of how the presence of a cold, dense ISM affects cluster formation and evolution. We consider two different star formation prescriptions: a constant star formation efficiency per free-fall time, as well as an environmentally-dependent, turbulence-based prescription. We identify two key results drawn from these simulations. Firstly, we find that tidal shock-driven disruption caused by the graininess of the cold ISM produces old ($τ>10~$Gyr) stellar cluster populations with properties that are in excellent agreement with the observed populations in the Milky Way and M31. Importantly, the addition of the cold ISM addresses the areas of disagreement found in previous simulations that lacked the cold gas phase. Secondly, the formation of stellar clusters is extremely sensitive to the baryonic physics that govern the properties of the cold, dense gas reservoir in the galaxy. This implies that the demographics of stellar cluster populations represent an important diagnostic tool for constraining baryonic physics models in upcoming galaxy formation simulations that also include a description of the cold ISM. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: 30 pages + appendices, 12 figures. MNRAS submitted, comments and suggestions welcome!

The Initial Properties of Young Star Clusters in M83

Authors: Jeremy J. Webb, Alison Sills

Abstract: The initial sizes and masses of massive star clusters provide information about the cluster formation process and also determine how cluster populations are modified and destroyed, which have implications for using clusters as tracers of galaxy assembly. Young massive cluster populations are often assumed to be unchanged since cluster formation, and therefore their distribution of masses and radii… ▽ More The initial sizes and masses of massive star clusters provide information about the cluster formation process and also determine how cluster populations are modified and destroyed, which have implications for using clusters as tracers of galaxy assembly. Young massive cluster populations are often assumed to be unchanged since cluster formation, and therefore their distribution of masses and radii are used as the initial values. However, the first few hundred million years of cluster evolution does change both cluster mass and cluster radius, through both internal and external processes. In this paper, we use a large suite of $N$-body cluster simulations in an appropriate tidal field to determine the best initial mass and initial size distributions of young clusters in the nearby galaxy M83. We find that the initial masses follow a power-law distribution with a slope of -2.7 $\pm$ 0.4, and the half-mass radii follow a log-normal distribution with a mean of 2.57 $\pm$ 0.04 pc and a dispersion of 1.59 $\pm$ 0.01 pc. The corresponding initial projected half-light radius function has a mean of 2.7 $\pm$ 0.3 pc and a dispersion of 1.7 $\pm$ 0.2 pc. The evolution of the initial mass and size distribution functions are consistent with mass loss and expansion due to stellar evolution, independent of the external tidal field and the cluster's initial density profile. Observed cluster sizes and masses should not be used as the initial values, even when clusters are only a few hundred million years old. △ Less

Submitted 14 December, 2020; originally announced December 2020.

Comments: 8 pages, 6 figures, accepted for publication in MNRAS

Implementing Primordial Binaries in Simulations of Star Cluster Formation with a Hybrid MHD and Direct N-Body Method

Authors: Claude Cournoyer-Cloutier, Aaron Tran, Sean Lewis, Joshua E. Wall, William E. Harris, Mordecai-Mark Mac Low, Stephen L. W. McMillan, Simon Portegies Zwart, Alison Sills

Abstract: The fraction of stars in binary systems within star clusters is important for their evolution, but what proportion of binaries form by dynamical processes after initial stellar accretion remains unknown. In previous work, we showed that dynamical interactions alone produced too few low-mass binaries compared to observations. We therefore implement an initial population of binaries in the coupled M… ▽ More The fraction of stars in binary systems within star clusters is important for their evolution, but what proportion of binaries form by dynamical processes after initial stellar accretion remains unknown. In previous work, we showed that dynamical interactions alone produced too few low-mass binaries compared to observations. We therefore implement an initial population of binaries in the coupled MHD and direct N-body star cluster formation code Torch. We compare simulations with, and without, initial binary populations and follow the dynamical evolution of the binary population in both sets of simulations, finding that both dynamical formation and destruction of binaries take place. Even in the first few million years of star formation, we find that an initial population of binaries is needed at all masses to reproduce observed binary fractions for binaries with mass ratios above the $q \geq 0.1$ detection limit. Our simulations also indicate that dynamical interactions in the presence of gas during cluster formation modify the initial distributions towards binaries with smaller primary masses, larger mass ratios, smaller semi-major axes and larger eccentricities. Systems formed dynamically do not have the same properties as the initial systems, and systems formed dynamically in the presence of an initial population of binaries differ from those formed in simulations with single stars only. Dynamical interactions during the earliest stages of star cluster formation are important for determining the properties of binary star systems. △ Less

Submitted 11 November, 2020; originally announced November 2020.

Comments: 15 pages, 14 figures, submitted to MNRAS and edited to address positive referee's report

Chemical abundances along the 1G sequence of the chromosome maps: The Globular Cluster NGC 3201

Authors: A. F. Marino, A. P. Milone, A. Sills, D. Yong, A. Renzini, L. R. Bedin, G. Cordoni, F. D'Antona, H. Jerjen, A. Karakas, E. Lagioia, G. Piotto, M. Tailo

Abstract: 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 popu… ▽ More 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. △ Less

Submitted 7 October, 2019; originally announced October 2019.

Comments: 24 pages, 12 figures, 8 tables. Accepted for publication in ApJ

Light element variations within the different age-metallicity populations in the nucleus of the Sagittarius dwarf

Authors: Alison Sills, Emanuele Dalessandro, Mario Cadelano, Mayte Alfaro-Cuello, J. M. Diederik Kruijssen

Abstract: 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… ▽ More 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. △ Less

Submitted 30 September, 2019; originally announced October 2019.

Comments: Accepted by MNRAS Letters

The Emergence of Multiple Populations in Radiation Hydrodynamics Simulations of Cluster Formation

Authors: Alison Sills

Abstract: 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 mechanism… ▽ More 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. △ Less

Submitted 25 July, 2019; originally announced July 2019.

Comments: 4 pages, 1 figure, to appear in Proceedings IAU Symposium No. 351 Star Clusters: From the Milky Way to the Early Universe, 2019

Constraining Mass Transfer Histories of Blue Straggler Stars with COS Spectroscopy of White Dwarf Companions

Authors: Natalie M. Gosnell, Emily M. Leiner, Robert D. Mathieu, Aaron M. Geller, Christian Knigge, Alison Sills, Nathan W. C. Leigh

Abstract: 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… ▽ More 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. △ Less

Submitted 19 August, 2019; v1 submitted 3 April, 2019; originally announced April 2019.

Comments: 16 pages, 10 figures, 2 tables, submitted to Astrophysical Journal

Discovery of a double Blue Straggler sequence in M15: new insight into the core-collapse process

Authors: Giacomo Beccari, Francesco R. Ferraro, Emanuele Dalessandro, Barbara Lanzoni, Silvia Raso, Livia Origlia, Enrico Vesperini, Jongsuk Hong, Alison Sills, Andrea Dieball, Christian Knigge

Abstract: 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 t… ▽ More 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. △ Less

Submitted 26 March, 2019; originally announced March 2019.

Comments: Accepted for publication on ApJ

Blue straggler stars beyond the Milky Way. IV. Radial distributions and dynamical implications

Authors: Chengyuan Li, Weijia Sun, Jongsuk Hong, Licai Deng, Richard de Grijs, Alison Sills

Abstract: 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 observa… ▽ More 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. △ Less

Submitted 17 December, 2018; originally announced December 2018.

Comments: 16 pages, 15 figures, accepted by the ApJ

Observations of Spin-Down in Post-Mast-Transfer Stars and the Possibility for Blue Straggler Gyrochronology

Authors: Emily Leiner, Robert Mathieu, Natalie Gosnell, Alison Sills

Abstract: 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 fo… ▽ More 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. △ Less

Submitted 5 December, 2018; originally announced December 2018.

Comments: 7 Pages, 1 Figure, 1 Table, accepted to ApJ Letters

Constraints on blue straggler formation mechanisms in Galactic globular clusters from proper motion velocity distributions

Authors: N. W. C. Leigh, T. Panurach, M. Simunovic, A. M. Geller, D. Zurek, M. M. Shara, A. Sills, C. Knigge, N. M. Gosnell, R. Mathieu, T. H. Puzia, J. Ventura, Q. Minor

Abstract: 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… ▽ More 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. △ Less

Submitted 3 October, 2018; originally announced October 2018.

Comments: 10 pages, 8 figures, 1 table; accepted for publication in MNRAS

The Effect of Stellar Helium Abundance on Dynamics of Multiple Populations in Globular Clusters

Authors: Amy Fare, Jeremy Webb, Alison Sills

Abstract: 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 generat… ▽ More 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. △ Less

Submitted 4 September, 2018; originally announced September 2018.

Comments: Accepted by MNRAS

On the Origin of Multiple Populations During Massive Star Cluster Formation

Authors: Corey S. Howard, Ralph E. Pudritz, Alison Sills, William E. Harris

Abstract: 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 accret… ▽ More 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. △ Less

Submitted 5 April, 2019; v1 submitted 21 August, 2018; originally announced August 2018.

Comments: Accepted to MNRAS. 11 pages, 5 figures

Kinematics in Young Star Clusters and Associations with Gaia DR2

Authors: Michael A. Kuhn, Lynne A. Hillenbrand, Alison Sills, Eric D. Feigelson, Konstantin V. Getman

Abstract: 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… ▽ More 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. △ Less

Submitted 11 November, 2018; v1 submitted 5 July, 2018; originally announced July 2018.

Comments: Accepted for publication in the Astrophysical Journal. Supplementary material is available at https://drive.google.com/drive/folders/1Cis_YIbItQ4sbItLUGz_BayOkkCVVkfM?usp=sharing. Changes: A typo in the distance to the Carina Nebula on page 5 was corrected

The internal rotation of globular clusters revealed by Gaia DR2

Authors: P. Bianchini, R. P. van der Marel, A. del Pino, L. L. Watkins, A. Bellini, M. A. Fardal, M. Libralato, A. Sills

Abstract: 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 ro… ▽ More 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/σ\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. △ Less

Submitted 30 August, 2018; v1 submitted 7 June, 2018; originally announced June 2018.

Comments: Accepted for publication in MNRAS

Dynamical evolution of stars and gas of young embedded stellar sub-clusters

Authors: Alison Sills, Steven Rieder, Jennifer Scora, Jessica McCloskey, Sarah Jaffa

Abstract: 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 t… ▽ More 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. △ Less

Submitted 12 March, 2018; originally announced March 2018.

Comments: Accepted for publication in MNRAS. 11 pages, 13 figures, 2 tables

Kinematic fingerprint of core-collapsed globular clusters

Authors: P. Bianchini, J. J. Webb, A. Sills, E. Vesperini

Abstract: 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… ▽ More 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. △ Less

Submitted 23 January, 2018; originally announced January 2018.

Comments: Accepted for publication in MNRAS Letter

The Structure of the Young Star Cluster NGC 6231. II. Structure, Formation, and Fate

Authors: Michael A. Kuhn, Konstantin V. Getman, Eric D. Feigelson, Alison Sills, Mariusz Gromadzki, Nicolás Medina, Jordanka Borissova, Radostin Kurtev

Abstract: 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… ▽ More 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. △ Less

Submitted 4 October, 2017; originally announced October 2017.

Comments: Accepted for publication in The Astronomical Journal; 27 pages, 2 tables, and 14 figures

Characterization of the velocity anisotropy of accreted globular clusters

Authors: Paolo Bianchini, Alison Sills, Meghan Miholics

Abstract: 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),… ▽ More 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. △ Less

Submitted 29 June, 2017; originally announced June 2017.

Comments: Accepted for publication on MNRAS

A tight relation between the age distributions of stellar clusters and the properties of the interstellar medium in the host galaxy

Authors: Meghan Miholics, J. M. Diederik Kruijssen, Alison Sills

Abstract: 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, p… ▽ More 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 Σ^{-0.51\pm0.03}σ_{\rm 1D}^{-0.85\pm0.10}M_{\rm min}^γ$, for the gas surface density $Σ$, gas velocity dispersion $σ_{\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}^γ$ provides a model-independent way to measure the mass dependence of the cluster disruption time $γ$. 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. △ Less

Submitted 30 May, 2017; originally announced May 2017.

Comments: 15 pages, 6 figures, 3 tables; accepted by MNRAS (May 24, 2017)

The relation between the mass-to-light ratio and the relaxation state of globular clusters

Authors: Paolo Bianchini, Alison Sills, Glenn van de Ven, Anna C. Sippel

Abstract: 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 th… ▽ More 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. △ Less

Submitted 5 May, 2017; originally announced May 2017.

Comments: Accepted for publication on MNRAS

On the Origin of Sub-subgiant Stars. III. Formation Frequencies

Authors: Aaron M. Geller, Emily M. Leiner, Sourav Chatterjee, Nathan W. C. Leigh, Robert D. Mathieu, Alison Sills

Abstract: 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 ob… ▽ More 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. △ Less

Submitted 17 May, 2017; v1 submitted 29 March, 2017; originally announced March 2017.

Comments: 14 pages, 6 figures, 1 table, accepted for publication in ApJ

On the Origin of Sub-subgiant Stars. I. Demographics

Authors: Aaron M. Geller, Emily M. Leiner, Andrea Bellini, Robert Gleisinger, Daryl Haggard, Sebastian Kamann, Nathan W. C. Leigh, Robert D. Mathieu, Alison Sills, Laura L. Watkins, David Zurek

Abstract: 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 impor… ▽ More 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. △ Less

Submitted 29 March, 2017; originally announced March 2017.

Comments: 21 pages, 7 figures, 4 tables, accepted for publication in ApJ

Globular Cluster Scale Sizes in Giant Galaxies: Orbital Anisotropy and Tidally Under-filling Clusters in M87, NGC 1399, and NGC 5128

Authors: Jeremy J. Webb, Alison Sills, William E. Harris, Matías Gómez, Maurizio Paolillo, Kristin A. Woodley, Thomas H. Puzia

Abstract: 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… ▽ More 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_β$ beyond which cluster orbits become anisotropic, as well as the distance $R_{fα}$ beyond which clusters are tidally under-filling. For M87 we find $R_β> 27$ kpc and $20 < R_{fα} < 40$ kpc and for NGC 1399 $R_β> 13$ kpc and $10 < R_{fα} < 30$ kpc. The connection of $R_{fα}$ 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_β> 5$ kpc and $R_{fα} > 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. △ Less

Submitted 10 May, 2016; originally announced May 2016.

Comments: 16 pages, 7 figures, 4 tables, Accepted for publication in MNRAS

Observational constraints on star cluster formation theory - I. The mass-radius relation

Authors: S. Pfalzner, H. Kirk, A. Sills, J. S. Urquhart, J. Kauffmann, M. A. Kuhn, A. Bhandare, K. M. Menten

Abstract: 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 clus… ▽ More 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^γ$ 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) △ Less

Submitted 1 December, 2015; originally announced December 2015.

Comments: 13 pages, 7 figures, accepted by A&A

The Dynamical Evolution of Accreted Star Clusters in the Milky Way

Authors: Meghan Miholics, Jeremy J. Webb, Alison Sills

Abstract: 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… ▽ More 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. △ Less

Submitted 15 November, 2015; v1 submitted 12 November, 2015; originally announced November 2015.

Comments: 9 pages, 4 figures, accepted for publication in MNRAS

A Hubble Space Telescope Study of the Enigmatic Milky Way Halo Globular Cluster Crater

Authors: Daniel R. Weisz, Sergey E. Koposov, Andrew E. Dolphin, Vasily Belokurov, Mark Gieles, Mario L. Mateo, Edward W. Olszewski, Alison Sills, Matthew G. Walker

Abstract: 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… ▽ More 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] △ Less

Submitted 3 November, 2015; v1 submitted 28 October, 2015; originally announced October 2015.

Comments: 9 pages, 11 Figures. Submitted to ApJ. Comments welcome

The Properties of Bound and Unbound Molecular Cloud Populations Formed in Galactic Disc Simulations

Authors: Rachel L. Ward, Samantha M. Benincasa, James Wadsley, Alison Sills, H. M. P. Couchman

Abstract: 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 scal… ▽ More 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. △ Less

Submitted 28 October, 2015; originally announced October 2015.

Comments: 12 pages, 7 figures, Accepted for publication by MNRAS

Implications for the Formation of Blue Straggler Stars from HST Ultraviolet Observations of NGC 188

Authors: Natalie M. Gosnell, Robert D. Mathieu, Aaron M. Geller, Alison Sills, Nathan Leigh, Christian Knigge

Abstract: 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… ▽ More 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%. △ Less

Submitted 14 October, 2015; originally announced October 2015.

Comments: 13 pages, 6 figures, accepted to the Astrophysical Journal

The Dynamics of Multiple Populations in the Globular Cluster NGC 6362

Authors: Meghan Miholics, Jeremy J. Webb, Alison Sills

Abstract: 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 concentra… ▽ More 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. △ Less

Submitted 25 September, 2015; v1 submitted 8 September, 2015; originally announced September 2015.

Comments: 8 pages, 5 figures, Accepted for publication in MNRAS

The Spatial Structure of Young Stellar Clusters. III. Physical Properties and Evolutionary States

Authors: Michael A. Kuhn, Eric D. Feigelson, Konstantin V. Getman, Alison Sills, Matthew R. Bate, Jordanka Borissova

Abstract: 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, an… ▽ More 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. △ Less

Submitted 20 July, 2015; originally announced July 2015.

Comments: Accepted for publication in The Astrophysical Journal ; 48 pages, 13 figures, and 6 tables

The state of globular clusters at birth II: primordial binaries

Authors: Nathan W. C. Leigh, Mirek Giersz, Michael Marks, Jeremy J. Webb, Arkadiusz Hypki, Craig O. Heinke, Pavel Kroupa, Alison Sills

Abstract: (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 a… ▽ More (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}^α$ with $α<$ 1/3. △ Less

Submitted 8 October, 2014; originally announced October 2014.

Comments: 17 pages; 5 figures; 1 table; accepted for publication in MNRAS

The Blue Straggler Star Population in NGC 1261: Evidence for a Post-Core-Collapse Bounce State

Authors: Mirko Simunovic, Thomas H. Puzia, Alison Sills

Abstract: 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 si… ▽ More 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. △ Less

Submitted 2 October, 2014; originally announced October 2014.

Comments: 5 pages, 5 figures, accepted for publication in ApJL

The Size of Star Clusters Accreted by the Milky Way

Authors: Meghan Miholics, Jeremy J. Webb, Alison Sills

Abstract: 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 decreas… ▽ More 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. △ Less

Submitted 17 September, 2014; originally announced September 2014.

Comments: 7 pages, 6 figures, Accepted for publication MNRAS

Evolving Molecular Cloud Structure and the Column Density Probability Distribution Function

Authors: Rachel L. Ward, James Wadsley, Alison Sills

Abstract: 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… ▽ More 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. △ Less

Submitted 8 September, 2014; originally announced September 2014.

Comments: 10 pages, 9 figures, Accepted for publication by MNRAS

The Effects of Orbital Inclination on the Scale Size and Evolution of Tidally Filling Star Clusters

Authors: Jeremy J. Webb, Alison Sills, William E. Harris, Jarrod R. Hurley

Abstract: 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 e… ▽ More 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. △ Less

Submitted 2 September, 2014; originally announced September 2014.

Comments: 9 pages, 9 figures, Accepted for publication in MNRAS

Models of Individual Blue Stragglers

Authors: Alison Sills

Abstract: 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 c… ▽ More 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. △ Less

Submitted 13 June, 2014; originally announced June 2014.

Comments: Chapter 12, in Ecology of Blue Straggler Stars, H.M.J. Boffin, G. Carraro & G. Beccari (Eds), Astrophysics and Space Science Library, Springer

The Effect of Orbital Eccentricity on the Dynamical Evolution of Star Clusters

Authors: Jeremy J. Webb, Nathan Leigh, Alison Sills, William E. Harris, Jarrod R. Hurley

Abstract: 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… ▽ More 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. △ Less

Submitted 4 July, 2014; v1 submitted 13 May, 2014; originally announced May 2014.

Comments: 10 pages, 7 figures, Accepted for publication in MNRAS

The Spatial Structure of Young Stellar Clusters. I. Subclusters

Authors: Michael A. Kuhn, Eric D. Feigelson, Konstantin V. Getman, Adrian J. Baddeley, Patrick S. Broos, Alison Sills, Matthew R. Bate, Matthew S. Povich, Kevin L. Luhman, Heather A. Busk, Tim Naylor, Robert R. King

Abstract: 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… ▽ More 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. △ Less

Submitted 17 March, 2014; originally announced March 2014.

Comments: Accepted for publication in The Astrophysical Journal; 67 pages, 9 figures, and 4 tables. Supplemental materials from the online-version of this article, including machine-readable tables and source code in R are available at the public MYStIX website (http://astro.psu.edu/mystix/)