Open star clusters play a key role in various aspects of Galactic Archaeology. They are formed by groups of stars born from the same molecular cloud, thus sharing the same age, kinematics, and chemical composition; therefore their properties can be measured much more precisely than those of individual stars.
Moreover, they are among the best test particles to prove the conditions of galaxies both at present and back in time, being located across all the Galactic disc and covering a wide range of ages. In addition, star formation is believed to occur mainly in areas of high density, while stars rarely form in isolation; the majority of clusters then disrupt contributing to the population of field stars. Therefore, open clusters stand at the intersection between the structure of the Galactic disk and the stellar populations it contains.
However, the transition between gravitationally bound stellar populations and field populations is not fully understood, and many open questions still stand: what is the rate of cluster destruction and how does it vary across the Galactic disc? how differently field stars and open clusters trace the Galaxy in space and time? How representative is the current population of open clusters in our Galaxy, especially the oldest ones, of the initial cluster population?
This workshop aims at combining the presentation of new observational results with numerical simulations of cluster disruption and migration. The final goal of the meeting will be to bring an innovative and critical input into the panorama of large spectroscopic surveys and Gaia satellite results, allowing a more detailed understanding of the Galactic stellar populations.
Note:
All talks will last 20 min+ 5 min for questions. At the end of each block there will be a discussion session.
Cluster formation: star formation models, spatial and temporal variation of the rate of clusters formation
Cluster demographics: the Gaia mission, new identification techniques based on artificial intelligence, description of new samples, complete characterisation of the existing ones; presence of halo and tidal tails in clusters.
Cluster disruption and stellar migration: cluster disruption models, role of stellar migration in the redistribution of populations, chemical and kinematic evidence of radial and vertical displacement, differences between field stars and clusters, numerical models of cluster migration and disruption rates.
Clusters in current and next generation instruments/surveys: what are the key open questions that we could tackle with the next generation of instruments and surveys, such as , as e.g., WEAVE@WHT, 4MOST@VISTA, MOONS@VLT
g
Welcome by the Director of INAF-Osservatorio di Arcetri. Brief introduction to the workshop by the organisers.
Galactic studies are currently undergoing a renaissance, thanks to the wealth of data from the Gaia satellite and ground-based surveys. Previously unknown portions of the Galactic disc have been mapped, triggering new interest and questions into the physical mechanisms regulating the evolution of the Galaxy. Recent works based on field stars and open clusters revealed: (i) the segments of the nearest spiral arms out to 5 kpc in heliocentric distance; (ii) substructures in velocity space; (iii) spatially-dependent chemical substructures. In this talk, I will compare the spatial, kinematic and chemical properties of field stars and open clusters, to investigate whether the evolutionary processes that led to their current configuration acted in a similar way.
With the unprecedented increase of known star clusters, quick and modern tools are needed for their analysis. In this work, we develop an artificial neural network trained on synthetic clusters to estimate the age, metallicity, extinction, and distance of Gaia open clusters. We implement a novel technique to extract features from the colour-magnitude diagram of clusters by means of the QuadTree tool and we adopt a multi-band approach. We obtain reliable parameters for ~5400 clusters. We demonstrate the effectiveness of our methodology in accurately determining crucial parameters of Gaia open clusters by performing a comprehensive scientific validation. In particular, with the parameters produced by our neural network, we obtain a Galactic metallicity gradient in agreement with the observed ones thus it demonstrates that our method reliably extracts information on metallicity from color-magnitude diagrams (CMDs) of stellar clusters. For the sample of clusters studied, we find an intriguing systematic older age compared to previous analyses present in the literature. This work introduces a novel approach to feature extraction using a QuadTree algorithm, effectively tracing sequences in CMDs despite photometric errors and outliers. The adoption of ANNs, rather than Convolutional Neural Networks, maintains the full positional information and improves performance, while also demonstrating the potential for deriving clusters' parameters from simultaneous analysis of multiple photometric bands, beneficial for upcoming telescopes like the Vera Rubin Observatory. The implementation of ANN tools with robust isochrone fit techniques could provide further improvements in the quest for open clusters' parameters.
Gaia has revolutionised our ability to identify clusters and to use them as tracers of the Galactic structure. Our view is however severely limited when looking towards the inner regions of the Milky Way, where the clusters we do see are very young and very massive. The differences between the inner and outer disc are in part due to observational biases, but could also indicate that different environments lead to different preferential modes of star formation. Upcoming mid-infrared astrometric missions such as JASMINE and GaiaNIR will allow us to probe the inner disc and shed light on the processes driving cluster formation and dispersal.
The combination of photometric and astrometric data from Gaia with parameters derived by spectroscopic surveys, like Gaia-ESO, improved our knowledge of the properties of young star clusters and our understanding of the processes leading the cluster evolution until its dispersion.
However, current datasets suffer of two main shortcomings: they are mostly based on optical observations, so we are not able to investigate the earliest stages of cluster life, when they are still embedded in molecular clouds, and spectroscopic data coming from multi-object spectrographs are usually limited to the inner and most dense part of the cluster. In this talk, I will present two spectroscopic surveys that will be carried out with the new multi-object spectrographs MOONS at the VLT and 4MOST at VISTA.
The first will gather infrared spectra of about 10 embedded young clusters to investigate their structural and kinematic properties at their earliest stages. While the 4MOST survey will observe ~100,000 candidate young stars within 500 pc selected independently from their position and kinematic to investigate the properties of young unbound stellar populations.
The study of open clusters is a crucial element in the understanding of the processes that led to the formation of the disk in the Milky Way and in general in spiral galaxies. For this reason they have been targeted by all the recent spectroscopic surveys, providing radial velocities and chemical abundances. In the coming future, a number of spectroscopic surveys such as WEAVE and 4MOST will have dedicated open cluster surveys. In this talk I will briefly revise the current status of knowledge on chemical abundances and radial velocities from large spectroscopic surveys. I will focus on the WEAVE and 4MOST open and young cluster survey, providing updates on the status of the projects.
To study stellar clusters and be able to effectively use them to understand how stars and clusters form and evolve and to link this to the Galactic evolution, we need not only large samples, such as those provided by the spectroscopic surveys, but also high resolution and detailed analysis.
As part of the SPA (Stellar Population Astrophysics) Large Programme at the TNG, we observed a large sample of open clusters with GIARPS (GIANO+HARPS-N), obtaining a few to a few tens of stars in each of them. I will present the aims, some highlights of our work, some work in progress and will underline how such project may be very useful also in a large-survey era.
Stars are primarily formed in clustered environments in giant molecular clouds. Stars are also primarily found in binary or higher order multiple systems. Therefore, there is a time during the formation of star clusters when both stellar dynamics and hydrodynamics are needed to describe the dominant populations, and we probably also need to care about radiation, chemistry, dust physics, and magnetic fields at some level too. In this talk, I will describe some recent simulations of star cluster formation and show how a more complete treatment of stellar and gas physics, and their interaction, are necessary to properly model star cluster formation, and can lead to "common sense" solutions to some of the puzzles that remain about the properties of star clusters.
While some young stars are clustered, many other young stars are observed as dispersed,
non-gravitationally bound groups, or in complete isolation. Using N-body simulations,
we study dispersal of young open star clusters of a wide range of initial masses.
We find that observational data are consistent with the assumption that all stars form
in initially embedded star clusters, which loose most of their members at their early
age as the consequence of expulsion of the residual gas. Our model predicts that
the fraction of young stars (age less than 10 Myr) to be observed in gravitationally
bound clusters shows only a weak dependence on the star formation rate per unit area
of the galaxy. This result contrasts to the analytical work due to Kruijssen 2012,
which predicts that the fraction of young stars in gravitationally bound clusters is
a strong function of the star formation rate of the galaxy, and that only the galaxies
with high rates of star formation form a noticeable fraction of stars in clusters.
In the last ~10 yrs, several young stellar clusters have been associated to diffuse gamma-ray emission in the energy range from ~1 GeV up to ~100 TeV. In particular in the region close to the association Cygnus OB2 photons up to the incredible energy od 1.4 PeV have been detected by LHAASO. The orginin of the non-thermal particles responsible for the gamma-ray emission is not completely clear. Here we present a model based on particle acceleration at the wind termination shock of the collective wind generated by massive stars in the cluster. This model requires as an imput the wind kinetic power, which depends on the IMF and age of the clusters. We will critically discuss those aspects.
The initial conditions of star-forming regions will dictate if they will disperse into the field or survive as bound open clusters. Methods have been developed to characterise star-forming regions and infer the initial conditions of them (i.e. initial density, degree of substructure and virial state). These methods have been used to quantify the spatial clustering of stars within star-forming regions and the overall morphologies of such regions. One such method is the Mahalanobis density, which for stars quantifies a 6D (position-velocity) phase space density, which has been used to infer the initial formation conditions of exoplanet host stars. I have applied the Mahalanobis density to simulations with different initial conditions and find that it is unable to reliably infer the initial conditions of these regions when used by itself. I will show that using multiple methods together enables the initial conditions of star-forming regions to be better constrained.
In our quest to gain a deeper understanding of our Solar Neighborhood, Gaia has proven instrumental in elucidating the 3D spatial structure of the local interstellar medium (ISM) and the distribution of young stellar clusters. However, to unravel the origin and evolution of nearby young structures, a crucial dimension is the measurement of their 3D space motions. These motions allow us to trace back the orbits of stellar clusters and molecular clouds within our Galaxy, and investigate the relative space motions within single complexes.
To address these challenges, we have recently developed a machine learning-based clustering tool (SigMA) tailored for the selection of stellar clusters in the 5D phase space as provided by Gaia. This innovative approach was initially applied to the Scorpio-Centaurus complex, leading to the identification of previously unnoticed substructures and richer stellar populations than earlier established. Combined with updated cluster ages, this framework enables a more comprehensive investigation of the history of this nearby region, including the evolution of velocity dispersion within a single complex and the propagation of star formation influenced by feedback from massive stars.
To achieve this we analyze the 6D phase-space of individual entities such as clusters and clouds, utilizing Gaia DR3 data in combination with ancillary radial velocity (RV) data. This approach allows us to reconstruct the formation history of star-forming complexes, also shedding light on past feedback processes. Consequently, we can study the interaction of stars with the ISM and the formation and evolution of feedback-driven bubbles, (e.g., Orion-Eridanus Superbubble, Local Bubble). Hence, our research provides a quantifiable assessment of the impact of feedback from massive stars on nearby regions, contributing to a more comprehensive understanding of star formation in the Milky Way.
We perform simulations of cluster formation in regions taken from galaxy scale simulations, including photoionization and supernovae feedback. We simulate regions with different densities, and from different galactic environments. In all our simulations, clusters undergo mergers and splits during their formation. More massive clusters form in regions of spiral arms with stronger converging flows, and in bars and inner spiral arm regions. In inter-arm or outer galaxy regions, we tend to see looser groups more characteristic of associations. In all simulations, the most massive clusters are formed by mergers. Feedback has a greater impact on cluster formation at lower densities, where star formation occurs over a longer timescale. In our lowest density case, the resulting clusters (in terms of number, mass) are very different with and without feedback, and one of our simulations forms an association which has a similar size, mass and morphology to the Orion OB1 association. We also show that photoionizing feedback is necessary to produce clusters with the observed cluster mass relation. Initial results with magnetic fields show that magnetic fields can have a strong impact on cluster formation, depending on field strength. Lastly, we place our simulated clusters into the GAIA catalogue and try to re-identify them. We find that the spatial properties of clusters are surprisingly robust even at relatively large distances and when relatively few members of the original cluster are observed. However dynamical measurements tend to be much less reliable.
Star formation is a fundamental process that impacts many fields of astrophysics, from the formation and evolution of planets to galaxies. The interaction between the natal cloud and the newborn stars is one of the least understood star formation processes and has an important impact on the final star-formation efficiency and cluster dynamics. I will present an innovative methodology to measure the timescale of the gas-embedded phase by comparing stellar ages derived with two independent methods: dynamical tracebacks and isochrone fitting. In this new framework, the dynamical-traceback “clock” initiates when a stellar cluster or association begins to expand after expelling most of the gas, while the isochronal “clock” initiates earlier when most stars form. Measuring this difference accurately and understanding its variations across different environments provides new information on the impact of local conditions and stellar feedback on the formation and dispersal of stellar clusters.
The formation and evolution of young star clusters and OB associations is fundamental to our understanding of the star formation process, the conditions faced by young binary and planetary systems, and the formation of long-lived open and globular clusters. The Gaia-ESO Survey (GES) has spectroscopically observed 18 young star forming regions, star clusters and OB associations, providing spectroscopic indicators of youth and radial velocities to complement Gaia DR3 proper motions. We have conducted a 3D dynamical study of ~2500 young stars across the 18 young star clusters observed by GES, combining Gaia proper motions with GES radial velocities and youth indicators. This is the first large-scale and comprehensive 3D dynamical study of multiple star clusters. We measure 3D velocity dispersions, quantifying virial states, and search for evidence of (and quantify) expansion, rotation, anisotropy and energy equipartition. We compare the dynamical properties of these systems with their mass, density and age, revealing trends that hint at the formation of these systems, their dynamical evolution, and their eventual dispersal into the field.
Young clusters are a formidable tool to study star formation, stellar evolution, binary stars, as well as the formation and evolution of clusters themselves. The availability of data from Gaia, coupled with additional ground-based data and clustering analysis methods, allow us for the first time to discover and study a wide range of clusters with unprecedented details. In the region around Gamma Vel, we have identified a 260pc long filamentary structure of coeval stars that connects several 30 Myr old associations and known clusters, possible relict of the formation of stars along a more that 100 pc-long filament. Similarly, we identified a long (90pc), thin (about 10pc width), co-eval (10Myr) stellar structure in the Orion star-forming region which is likely a relic of star formation along a molecular cloud filament, naming it the Orion relic filament. In both cases, we exclude that such populations originates by the same mechanism responsible to create tidal streams around older clusters. We are instead observing a constellations of young and coeval star clusters and associations few million years after their formation, bridged on a more than 100pc scale by a filament of coeval stars. Such structures are likely the outcome of a mechanism of filamentary star formation in a Giant Molecular Clouds. Clearly, thanks to Gaia, the study of young stellar clusters on a galactic scale is undergoing a complete revolution.
The EWOCS (Extended Westerlund 1 and 2 Open Clusters Survey) project has the objective of studying star and planet formation, and early stellar evolution, in super massive star clusters (SSCs). With a mass in excess of 10^4 solar masses, the very few SSCs known in the Milky Way represent the most accessible examples of starburst regions, which are very rare in our Galaxy today, but common in galaxies experiencing epochs of intense star formation. These regions are characterized by very high stellar density, and they are dominated by a rich and compact ensemble of massive stars that produce an environment dominated by energetic radiation and particles. With a distance of 3.87 kpc and 4.5 kpc, respectively, the Westerlund 1 and 2 clusters are the closest SSCs to the Sun, and thus the best targets to study how stars and planets form in the most energetic star forming environment known. In this talk, I will present the motivations, status and the preliminary results of the EWOCS project, which is mainly based on a 1Msec Chandra/ACIS-I Large Project and a cycle 1 JWST observation of Westerlund 1, a cycle 2 JWST observation of Westerlund 2, and other data at high spatial resolution of the two clusters.”
In just five years, Gaia has revolutionised the census of star clusters in the Milky Way - including allowing for the detection of thousands of new clusters. However, many of these clusters are much smaller and more dispersed than canonical open clusters like the Pleiades, and it is not clear how many of these clusters are actually gravitationally bound. This presents critical challenges when trying to interpret the results of cluster detection algorithms.
In practice, measuring the boundness of a star cluster observationally is challenging, and dynamical studies of star clusters have so far been limited to analysing no more than a few dozen reliable objects at a time. In this talk, I will present the results of a study quantifying the boundness of around 7000 clusters in the recent catalogue of Hunt & Reffert 2023. I will overview the methods I used, including how masses and dispersions were measured for these objects. I will also discuss which physical theories (such as the virial theorem) were the most useful in measuring the boundness of star clusters, as well as discussing their caveats and areas that require more work in the future. Finally, I will present the overall results of this work, which show that at least 21% of the clusters in Hunt & Reffert 2023 are not gravitationally bound clusters, with most of these unbound clusters being within a few hundred parsecs of the Sun. The methods developed in this work should have wide applications to other studies of the boundness of star clusters in the Gaia era.
Accurate mass determinations for open clusters (OCs) are essential for
studying star formation as well as OC dissolution processes. Using the Gaia DR3 OC catalogue of Hunt & Reffert (2023), we construct an extensive sample of OC mass estimates. We analyze more than 2000 OCs within 2 kpc, estimating their tidal masses and fitting their mass functions where possible (including the effect of unresolved binaries). Our analysis is in general agreement with previous Gaia studies of smaller samples, and indicates that the mass function break shifts towards higher masses as the cluster ages. We also show that the complex selection function of the Gaia DR3 OC census is the main missing ingredient for tighter measurements of the cluster formation and destruction rates in the solar vicinity.
Deducing the masses and mass functions of stellar clusters (SCs) is an important step to understanding their formation and evolution and to comprehending the formation of the Galaxy. This contribution explores, how unresolved binaries affect not only the deduction of these parameters for SCs, but also for their tails. If the binaries in SCs and their tails cannot be resolved, their masses are expected to be underestimated by up to 25%, while the power law index of the deduced mass function is underestimated by about 0.2. It is also found that since the Galactic field stars were born in SCs, populations of stars with large velocity dispersions would be expected to host larger binary fractions. However, the large scatter in the velocity dispersions makes it hard to observe this effect.
Fluorine remains an enigmatic element in galactic archaeology, with a scarcity that camouflages its importance in understanding chemical and stellar evolution. Various astrophysical sites, including massive stars, AGB stars, Wolf-Rayet stars, and Novae, have been proposed as fluorine producers, yet the element's origins remain a subject of debate. To shed light on this, we have initiated a specialized project that leverages the capabilities of the GIANO-B instrument at the TNG telescope to study fluorine's chemical evolution in the Galactic disk using open clusters as key tracers. Our research will uniquely provide fluorine abundance data as a function of both galactocentric distance and cluster age, a feat not achievable with current or forthcoming stellar surveys. By doing so, we aim to offer critical observational constraints to models of stellar and galactic chemical evolution. This project is poised to become the most exhaustive database on fluorine abundances in open clusters to date.
The recent detection of large diffuse $\gamma$-ray halos in coincidence with a dozen young massive stellar clusters (YMSCs) has strengthened the hypothesis of stellar clusters as cosmic ray factories and $\gamma$-ray sources. Noticeably, the observed $\gamma$-ray emission is, on average, remarkably extended (~1-3°) and of the same size as the dimension of the superbubble developed by YMSCs. Such a large diffuse emission, characterized by a low surface brightness, suggests that non-resolved YMSCs could significantly contribute to the diffuse $\gamma$-ray emission observed along the Galactic plane. To estimate this contribution, a robust modeling of both the galactic population of YMSCs and the feedback from stellar winds in terms of wind power and mass loss rate is required.
In this work, we present a novel method for simulating a synthetic population of galactic YMSCs based on observations of local stellar clusters. We additionally estimate the power of stellar wind and mass loss rate by building a mock stellar population for every cluster and empirically modeling the stellar wind physics.
The stars in a galaxy form in compact embedded star clusters which expand after removal of their residual gas. The subsequent virialisation leads to a fraction of the stars condensing into an open star cluster. The open star cluster dissolves as it orbits in the potential of the hosting galaxy through the energy equipartition process. In the Newtonian gravitational theory, the ensuing tidal tails contain about as many stars in the leading tail as in the trailing tail. In Milgromian gravitation the leading tail is predicted to contain more stars as a result of the non-linearity of this gravitational theory. The observed tidal tails of COIN-Gaia13, Coma Berenices, the Praesepe and the Hyades each contain more stars in the leading tail, implying the validity of Milgromian rather than Newtonian gravitation. I will also touch upon two new methods to assess the ages of open stars clusters using tidal tails from the gas expulsion process and the kissing instability in the interiors of M dwarf stars.
Compared to the globular clusters and their tidal tails, open cluster tails are ~1000 smaller in terms of star counts and must be found in an environment ~1000 times denser than the halo. Hence finding stars in tidal tails of open clusters is a difficult task even in the era of Gaia. I will present a probabilistic method where we start by simulating the dissolution of a cluster, which defines a likelihood of where the tidal tails should appear in some parameter space. We then select most likely stars belonging to the tidal tails from the Gaia catalogue as well as from a simulated population of stars. Comparing the two selections of stars allows us to define properly normalized membership probabilities for Gaia stars to be part of the tidal tails of a cluster. Normalized membership probabilities have not been calculated before. They allow quantitative measure for matching models of cluster formation and dissolution to observations and open several possibilities for more complete studies of stellar populations in clusters. We used Gaia's positions, distances, and proper motions to search for tidal tails of 500 open clusters and found tails up to thousands of parsecs long. I will also show the diversity of the shapes and structures of tidal tails and discuss the numerous dynamical processes that can be inferred from it.
The tidal tails of stellar clusters are an important tool for studying the clusters’ birth conditions, their evolution, coupling, and interaction with the Galactic potential, and to understand how field stars populate the Milky Way. Thanks to Gaia, much progress has been accomplished in finding tails of open clusters. I will show here that such tidal tails are much longer than previously observed, and that their identification requires not only a sophisticated analysis of the Gaia catalogue, using the convergent point method and clustering algorithms, but ideally, the use of N-body simulations and the new compact convergent point method. I will highlight recent results about the tails of several open clusters, which extend over several hundreds of parsecs.
We study the kinematics of the young stars around the spiral arms in the Milky Way, using Gaia DR3. We first demonstrate that we can measure the circular velocity at the position of the Sun in the Milky Way using the young OB stars, but there is a clear sign of the impact of the Local arm. From the kinematics of these stars and also of old stars, we find that the Local arm is not a minor arm, but a significant major spiral arm. In addition, the stellar kinematics indicates that the Local arm is a transient dynamic arm, and it is likely to be a growing phase. We also discuss the kinematics of the young stars around the Perseus and Outer arms. The kinematics of the young star clusters are likely to be also affected by the spiral arms in similar way to these young stars, and they will be a great chemo-dynamical tracer to further study of the nature of the spiral arm. Finally, we also briefly discuss the potential impact of the Local arm on the nearby star clusters.
I will discuss the results from numerical simulations of star cluster formation and evolution. I will review what has been learnt from numerical simulations regarding how the properties of low-mass stellar clusters and their stars depend on initial conditions, such as metallicity, molecular cloud density, magnetic fields, and turbulence. I will also discuss how these clusters may be expected to disperse into the galactic population of field stars, and possible kinematic signatures of dispersing stellar populations.
Understanding the factors that influence the survival of open clusters within the Milky Way is a complex endeavor. This study delves into the role of key parameters such as mass, density, size, and the galactic environment, including the Galactic bar, spiral structures, and molecular clouds. Particular attention is given to how open clusters and field stars respond differently to perturbations, with a special focus on the longevity of the oldest clusters. Our analysis compares kinematic properties between two Gaia DR3 datasets: ~40 open clusters and a vast sample of field stars. These datasets share quality criteria, thin disc membership, a similar metallicity range and locations. Employing a statistical approach, we explore the behavior of their kinematic properties. Our findings indicate that younger clusters (< 2-3 Gyr) exhibit greater resilience to perturbations, following quasi-circular orbits. In contrast, the surviving older clusters (> 3 Gyr) display more eccentric and inclined orbits.
The wealth and homogeneity of Gaia data have allowed the discovery of several open clusters with signs of disruption.
We do a morphological, kinematic and chemical analysis of the disrupting cluster UBC 274 (2.5 Gyr, d = 1778 pc), with the objective of studying its global properties.
A new membership study up to 50 pc from its center and un to magnitude G=19 using GaiaEDR3 data, shows that the cluster has a highly eccentric (0.93) component, tilted 10 deg with respect to the plane of the Galaxy, which is morphologically compatible with the result of a test-particle simulation of a disrupting cluster.
We find a significant sign of mass segregation where the most massive stars appear 1.5 times more concentrated than other stars.
We obtained high resolution and high signal-to-noise spectra of 6 giants and subgiants. Our abundance analysis shows that the cluster has a slightly subsolar metallicity of [Fe=H] = -0.08\pm0.02. Its chemical pattern is compatible with that of Ruprecht 147, of similar age but located closer to the Sun, with the remarkable exception of neutron-capture elements which present an overabundance.
M67 is a cornerstone open cluster for stellar astrophysicists. It is the same age and metallicity as the Sun - indeed its stars are so similar to the Sun that some have speculated that the Sun may have been born there. Furthermore its location - relatively close to use but well above the plane of the Milky Way - makes it easy to distinguish from the background field population. However this poses a question - how did M67 reach its present-day Galactic orbit? I will describe our investigations into different dynamical processes that affect cluster orbits in the Milky Way, how they could have acted to take M67 onto its present orbit, and the present-day orbits of the stars lost along the way. Our conclusion is that dynamics alone does not prevent a solar origin in M67, but makes it rather unlikely. Finally, I will briefly mention an interesting side-effect of M67's age and metallicity - the low mass of its surface convection zone - which makes it uniquely valuable as a place to study the ingestion of planets into stellar photospheres.
Numerical and analytical evidences that collisionless and collisional relaxation processes work differently in Newtonian and modified Newtonian dynamics (MOND) theories of gravitation imply that processes such as the disruption of satellites falling onto the parent galaxy should, in principle, behave differently in the two scenarios. In particular, it is not clear whether a stronger (barionic) dynamical friction in MOND systems could compensate or not the absence of the contribution of dark matter with respect to the parent (equivalent) Newtonian system.
Aiming at shedding some light on this issue, here we present some preliminary simulations with a hybrid stochastic-N-body method where the the internal dynamics of the stellar cluster is solved with a standard direct N-body scheme, while the effect of host galaxy is modelled, in both Newtonian and MOND gravities, with a smooth potential plus noise and dynamical friction (accounting for the discreteness effects) implemented via a Langevin-like method.
The sun is located at a Galactocentric distance where the transition from Newtonian to MoNDian dynamics is expected to occur. Due to their proximity solar neighbourhood open star clusters are ideal kinematical laboratories to discriminate between the validity of Newtonian or MoNDian dynamics. Direct MoNDian N-body integrations of open star clusters embedded in an external field are performed by the application of Milgrom’s law. The evolution of open star clusters in Newtonian and MoNDian dynamics show three main differences: i) The tidal tails are populated asymmetrically in MoND. The leading tail hosts more members than the trailing arm, whereas in Newton the tidal arms are populated nearly symmetrically. ii) Open star clusters dissolve approximately 30% faster in MoNDian dynamics than in Newtonian dynamics. iii) Open star cluster appear slightly super-virial when their internal kinematics is interpreted in Newtonian dynamics.