The effects of supernovae on the dynamical evolution of binary stars and star clusters

The effects of supernovae on the dynamical evolution of binary stars and star clusters

Richard J. Parker Richard J. Parker Astrophysics Research Institute,
Liverpool John Moores University,
146 Brownlow Hill,
Liverpool L3 5RF,
United Kingdom.
August 18, 2019

In this chapter I review the effects of supernovae explosions on the dynamical evolution of (1) binary stars and (2) star clusters.

(1) Supernovae in binaries can drastically alter the orbit of the system, sometimes disrupting it entirely, and are thought to be partially responsible for ‘runaway’ massive stars – stars in the Galaxy with large peculiar velocities. The ejection of the lower-mass secondary component of a binary occurs often in the event of the more massive primary star exploding as a supernova. The orbital properties of binaries that contain massive stars mean that the observed velocities of runaway stars (10s – 100s km s) are consistent with this scenario.
(2) Star formation is an inherently inefficient process, and much of the potential in young star clusters remains in the form of gas. Supernovae can in principle expel this gas, which would drastically alter the dynamics of the cluster by unbinding the stars from the potential. However, recent numerical simulations, and observational evidence that gas-free clusters are observed to be bound, suggest that the effects of supernova explosions on the dynamics of star clusters are likely to be minimal.

1 Introduction

Binary stars and star clusters are both ubiquitous outcomes of the star formation process. The collapse of a star-forming core usually transports angular momentum to the outskirts, and the fragmentation results in a binary star system. Furthermore, protostellar discs are prone to fragmentation, which produces a binary system, and binary stars can also form through capture Goodwin07 ; Moeckel10 . On larger scales, the collapse and fragmentation of Giant Molecular Clouds results in groupings of stars, a fraction of which then become bound star clusters Lada03 ; Kruijssen12b .

Massive stars are more likely to be found in binaries than their lower-mass counterparts, whereas most star clusters fully sample the stellar Initial Mass Function Bastian10 and therefore contain massive stars. In this chapter, I will discuss the potential effects of superovae explosions on the dynamical evolution of both binaries and star clusters.

2 Binary stars

The idea that supernovae can directly affect the dynamical evolution of binary star systems – two stars orbiting a common centre of mass – was first advocated by Zwicky in 1957 Zwicky57 and discussed in detail in two papers by Blaauw Blaauw61 and Boersma Boersma61 . These authors were attempting to explain the motion of massive stars (O- or B-type stars with masses 8M) in the Galaxy travelling with very high peculiar velocities, i.e. velocities that differ with respect to the local frame of reference. These fast-moving massive stars were dubbed ‘runaways’ as they could in some instances be traced back to star-forming regions and it was hypothesised that some mechanism had ejected them from these regions.

It is possible that close encounters with other stars can lead to the ejection of massive stars Leonard90 , but this scenario has two immediate caveats:

1) Firstly, ejecting stars of masses 8M requires an interaction with a star at least as massive, if not more so than the ejected star. This comes from the theory of three-body dynamics, where a star can be replaced in a system only if the incoming star has a signifcantly higher energy than the binding energy of the system. If more than three stars are involved in the interaction, this constraint is relaxed somewhat, as the dynamics becomes more chaotic and/or stochastic.

2) Regardless of the number of stars involved, the initial stellar density required to facilitate an interaction that would eject a massive star is likely to be well in excess of 1000 stars pc Allison11 ; Oh15 . Due to their inherently peculiar velocities, runaway OB stars can often be traced back to their natal star-forming regions. These star-forming regions are often so-called OB associations , which have much lower stellar densities than those required to dynamically eject massive stars. It is possible that OB associations were much more dense in the recent past Kroupa01a , although recent observations suggest that they are more likely to have formed with the same spatial and kinematic properties observed today Wright14 .

Given that not all runaway OB stars are likely to have originated in very dense stellar environments, some other disruption mechanism must be responsible for producing them. We will first consider the general effects of stellar evolution (i.e.  mass-loss) on the dynamical evolution of a binary star system, before discussing the effects of rapid, or instantaneous mass-loss due to a supernova.

2.1 Evolution of binaries due to gradual mass loss

If one of the stars in a binary system loses mass, then how will this affect the orbital configuration of the binary? According to Kepler’s third law, the period and semimajor axis are related to the component masses of the system thus:


where is the period of the system, the semimajor axis, the gravitational constant and and the respective masses of the stars. The binding energy of the system is


This binding energy is equivalent to the total energy of the system (potential energy plus kinetic energy). We can write this as the energy per unit reduced mass, , where the reduced mass is


and so the binding energy per unit reduced mass is


In order to assess the effects of mass loss on the orbit of the binary we write the above equation in terms of the semimajor axis:


We can make a differential equation of Eqn. 5 to study the changes in and :


Let us consider the kinetic and potential energies of the binary system:


where and are the kinetic and potential energy per unit reduced mass, respectively. The kinetic energy term is:


where is the velocity of the secondary star, at a given point in the orbit. The potential energy term is:


where is the magnitude of the radius vector. We can now write Eqn. 7 as


In order to compute the value of we need to consider the change in mass with time for the two component stars, using the differential operator . and , hence:


Although it appears that depends on the exact point in the orbit , we can assume that changes in the orbit occur slowly, and that only the time average values are important:


Because we are averaging over all points in the orbit, we have:


and so


Dividing Equations 14 by 4, we obtain


Finally, assuming the secondary star is not losing mass, then we place Eqn 15 into 6 to obtain


Inspection of Eqn. 16 shows that as the primary component loses mass ( is negative), the semimajor axis of the binary orbit increases. A similar analysis can also be performed to show that the period of the binary will also increase Huang56 .

If the mass-loss from one or both components causes the binary’s semimajor axis to increase, here are two implications for the future evolution of the binary system. Firstly, the overall binding energy is reduced, which makes in more susceptible to subsequent dynamical disruption from interactions with passing stars, especially if the binary is still in its (dense) natal star-forming region.

Secondly, the increase in period and semimajor axis leads to a decrease in orbital velocity. As well as the overall lower binding energy, the binary is more prone to encounters where a passing star can exchange into the binary (and eject one of the original stars).

2.2 Evolution of binaries due to rapid/instant mass loss

The above considerations apply to a gradual loss of mass in the binary, which describes the dynamical evolution of the system before any supernova. During the supernova, the more massive component loses mass almost instantaneously, and so the effect on the dynamical evolution of the system is quite different to the slow mass-loss case described above.

Let us re-write Eqn. 7 and drop the subscript ‘red’ (the terms still however, refer to the energy-per-reduced mass) and instead replace them with ‘o’, which refers to the total energy, , kientic energy and potential energy before mass loss:


The potential energy, is


where is the mass of the primary star immediately before the supernova, and is the radius of the orbit (assumed to be circular). The kinetic energy changes slowly as the orbit reacts to the sudden mass-loss, but the potential energy changes immediately.

The mass of the primary (i.e. supernova) decreases instantaneously by a factor , so the new mass of the primary is . The corresponding new potential energy, is thus:


and the change in potential energy is


The new total energy, , is therefore:


Assuming the orbit is circular (), the total energy of the original system before mass-loss () can be written:


If we substitute this into Eqn. 21, we have


which can be re-written as:


This equation describes the new energy of the binary system following the reduction in mass of the component that explodes as a supernova. A closed orbit (i.e. a bound binary system) has a total energy , whereas an open or hyperbolic orbit has a total energy . Therefore, in order to release the secondary component from the system,


or the reduction factor (recall that this is the ratio of post-supernova to pre-supernova mass of the exploding star) must fulfill the following inequality:


It is immediately apparent that must be less than 1/2 and this is only in the limit that ; in reality, the masses of the two components may actually be comparable (see Section 2.3 on binary demographics) and so the reduction factor would need to be even lower.

For an O-type star of initial mass M, the typical pre-supernova mass can be anywhere in the range 10 – 17M Meynet03 . Following the supernova, the remaining compact object (a neutron star) will have a mass of M, implying a reduction factor of . Stars with initial mass M will typically form a black hole with a higher mass – M, but still have a similar pre-supernova mass (10 – 17M) and so the reduction factor could be closer to .

The above considerations – although simplified through various assumptions – suggest that most supernovae in binaries result in the release of the secondary component.

Detailed numerical simulations have shown that the sudden mass-loss resulting from a supernova can often unbind the secondary star from the binary system Boersma61 ; Tauris98 ; Zwart00 . It is also possible, especially in the case of assymetric mass-loss from the primary component, that the velocity kick resulting from the supernova explosion can eject both the secondary star , and the resulting compact object , from the host star-forming region.

2.3 Binary demographics

In addition to the arguments above a pertinent question to ask is what fraction of massive stars may be affected by this mechanism? To answer this, we need to understand the fraction of massive stars that are in binary systems, and the distributions of their orbital semimajor axes and mass ratios.

The binary fraction of stars, is


where , , and are the numbers of single, binary or triple systems, respectively, in either a star-forming region, or a volume–limited region of the Galactic disc. For binary stars where the primary component is of similar mass to the Sun, the binary fraction is roughly 0.5 Raghavan10 . The orbital separation distribution can be approximated by a log-normal with a peak at 50 au and a variance , which means that the binaries can have semimajor axes anywhere between 10au and 10au. The distribution of mass ratios () is flat, i.e. a binary is equally likely to have a secondary component whose mass is 10 per cent of the mass of the primary, as opposed to a secondary component whose mass is equal to the primary mass Reggiani11a .

Unfortunately, the binary statistics for the most massive stars (i.e.  those that may contain a supernova that will affect the dynamics of the binary system) are hampered by the simple fact that very few binaries that contain massive stars exist. Recent research Mason09 ; Sana13 has gone some way to addressing this problem, mainly by targeting binaries containing massive stars in young massive star clusters . These environments are so young that the massive stars have not yet evolved (or gone supernovae), and because they are so massive they contain many massive stars.

Nevertheless, for binary systems that do contain massive stars it is apparent that their properties (overall fraction and orbital parameters) are very different to those of binaries where a solar-mass star is the primary component. A general consensus is that the binary fraction in young massive clusters is between 0.7 and 1.0 Mason09 ; Sana13 , and that their semimajor axis distribution is markedly different from lower-mass binaries, in that there is an excess of very close “spectroscopic” systems. Their semimajor axes follow an Opik Opik24 distribution, i.e. flat in log-space between 0 and 50 au. Furthermore, the mass ratio distribution appears to favour equal-mass systems (), meaning that an O-type star is more likely to be paired with either another O-type star, or a very massive B-type star. .

2.4 Properties of runaway stars

How do the binary properties of massive stars affect the expected velocity distribution of the runaway stars? According to calculations by Boersma61 ; Tauris98 ; Zwart00 , the velocity of the runaway star (the secondary component in the binary) is proportional to the orbital velocity of the binary, :


which for a supernova where the compact object that forms is a neutron star (M),


As before, is the mass of the primary star immediately before the supernova, and the mass of the secondary (runaway) star.

The orbital velocity, , at the point of the supernova is


where and are the pre-supernova masses of the primary and secondary stars, respectively, and is the semimajor axis.

Assuming the observed Opik distribution, then the maximum semimajor axis is 50 au. If we assume both stars have masses of 20M, then the velocity of the runaway after the supernova is 13km s. If the stars are both 40M, then the velocity is slightly higher, at 19km s. Adopting a much smaller semimajor axis (e.g. 1au) results in predicted velocities of order 100km s. However, massive stars that explode as supernovae lose a substantial fraction of their mass before the explosion. Therefore, the pre-supernovae masses and are likely to be in the region of 10 – 17M. Even with these values, a binary with initial semimajor axis 1 au is likely to impart a velocity of 100km s onto the secondary star after the supernova.

2.5 Dependence on star formation environment?

In determining the fraction of binary stars that could be disrupted by supernovae we must implicitly assume that each star forming region has the same binary star properties – a so-called ‘Universal population’ Kroupa95a . However, it is far from established that the binary fraction, and orbital characteristics are independent of environment. For example, in star-forming regions near the Sun the populations of binaries have been argued to be identical at birth, and observed differences are merely due to different amounts of dynamical processing. In this picture, the overall binary fraction and the semimajor axis distributions in two different regions can be altered to different degrees by the regions being born with different stellar densities (a high density region destroys more binaries compared to a lower-density region).

This is important for determining the fraction of binary systems that could be disrupted via supernovae. Dynamical interactions in dense star clusters act over crossing times – the distance travelled by stars in the cluster divided by the velocity; often taken to be the radius and average velocity. In very dense regions the crossing time can be as low as 0.1Myr, i.e. far shorter than the stellar evolution timescales for the most massive stars. Therefore, binaries containing massive stars could be disrupted through dynamical interactions with other massive stars before any supernova explosion.

Alternatively, other authors argue that based on constraints from the spatial distributions of stars in young star clusters, the initial binary population may vary between regions (possibly as a result of a different environment during the star formation process, Horton01 ; King12a ). Young star-forming regions exhibit a large degree of spatial substructureCartwright04 , and this is erased on order of a few crossing times due to dynamical interactions. Therefore, the amount of substructure in a star-forming region places an upper limit on the amount of dynamical interactions that may have taken placeParker14e . If the binary population does vary as a function of environment, then different star-forming regions may produce runaways due to supernovae in different proportions.

3 Star clusters

3.1 What is a star cluster?

Before discussing the effects of supernovae on the dynamical evolution of star clusters, it is worth a digression on the exact definition(s) that exist for classifiying clusters. In undergraduate courses we all learn that Globular Clusters are old (Gyr), massive ( M) dense stellar systems that are believed to have formed early in the evolution of the Universe. Their younger, less massive counterparts are often referred to as Open Clusters, and these have ages between 10 Myr – 1 Gyr, masses between M and fairly low densities (typically Mpc).

However, Type II supernovae occur in stars with progenitor masses 8M, and very massive O-type stars (20M) usually evolve (and explode as supernovae) on timescales less than 10Myr. At ages less than 10Myr, supernovae are therefore not found in Open or Globular clusters but rather in regions of recent and/or ongoing star formation.

As discussed in the previous section, a significant fraction of massive stars are observed in extended complexes (often several to 10s pc across) of star formation referred to as OB associationsBlaauw64 . How these OB associations form is currently the subject of debate. One school of thought is that they form essentially in situ, that is to say the the current observed density and spatial distribution of the stars (and velocity distribution, if measured) is very similar to that set by the initial conditions of star formation in the region in question. The second hypothesis is that OB associations are the expanded remnants of initially very compact (radii less than 1pc), dense groups of stars.

Young stars are typically found still embedded in the giant molecular clouds (GMCs) from which they form. Advances in infra-red detectors in the 1980s and 1990s meant that GMCs were studied in detail for the first time, and significant new groups of stars were discovered (see the review by Lada & Lada,Lada03 ). These were dubbed ‘embedded clusters’ and until recently the prevailing view in the literature was that almost all stars (70 – 90%) formed in these compact, embedded clusters.

However, recent observations suggest the picture may be somewhat more complicated. Alternative theories postulate the stars form in scale-free hierarchically self-similar spatial distributions Efremov95 ; Elmegreen01 ; Kruijssen12b . Some star-forming regions may undergo a dense phase, where dynamical interactions dominate and a dense cluster forms, but it is unclear what fraction of stars experience a clustered phase in their lifetime.

Despite this, we can still think of the estimate of the fraction of stars that form in clusters from Lada & Lada as an upper limit, as it is unclear which type of star-forming region (clustered versus non-clustered) contributes the most stars to the Galactic disc.

A further observation of embedded and young clusters is a sharp drop-off in the numbers of clusters as a function of their age. After 10 Myr, the number of clusters has decreased by 90% when compared to the number of clusters with ages of 1 Myr Lada03 ; Pellerin07 . This has been coined rather distastefully as “infant mortality” in the literature. The implication is that some mechanism(s) is destroying or disrupting the clusters as they age.

3.2 Star formation efficiency and gas expulsion

Even before the discovery of a decrease in the number of embedded clusters with age, it had been postulated that some processes act to destroy star clusters. Star formation is typically an inherently inefficient process. Stars form from the collapse and fragmentation of giant molecular clouds, but usually only around 30% of the gas mass is converted into stars. This means that 70% of the mass remains in the form of molecular gas and the ‘star formation efficiency’ is only 30% .

Given that such a substantial fraction of the system mass is in the form of gas, any process that would act on the gas but not the stars may cause a significant change to the system in question.

Stars appear to form with a mass distribution that is independent of environmental factors, such as the mass, density of gas and velocity motions in the gas. A more contentious issue is whether this mass distribution – the ‘Initial Mass Function’ , or IMF – is also independent of the age of the universe at the point of formation. Nevertheless, in nearby star clusters the IMF is statistically very similar Bastian10 and follows a roughly log-normal distribution for stellar masses ranging from the hydrogen burning limit, 0.08M, to around 1M. At higher masses, the IMF follows a power-law slope of the form


where - the so-called Salpeter slope Salpeter55 .

If we visualise this distribution, we immediately see than the probability of selecting a low-mass star is much higher than the probability of selecting a high-mass star. If we think of stellar mass distributions as an abstract maths problem, we now need to populate a star cluster of a given stellar mass with stars drawn from the IMF until we reach the cluster mass Elmegreen06 .

Therefore, on average a star cluster that contains several massive stars (20 M) will contain hundreds of other, low-mass stars. Massive stars evolve on timescales less than 10 Myr, so we expect they will lose a significant fraction of their mass during the early stages of a star cluster. However, if the IMF has been well-sampled (i.e. the cluster is also populated with many other (low-mass) stars), then the fraction of mass lost by the cluster due to the massive stars evolving and exploding as supernovae is negligible in the first 10 Myr.

However, recall that the vast majority (70% by mass) of gas is not converted into stars. Observations of HII regions – the immediate vicinities of the most massive stars – are often free of gas, which is interpreted as being due to the removal of gas by the intense stellar winds of the most massive stars Matzner02 ; Krumholz06 . The same effect is predicted to occur from the final supernova explosion.

What is the effect of this gas removal on the dynamical evolution of the system? Let us consider a simplistic view of the binding energy, of a star of mass in the system:


where is the distance between the star and the centre of mass of the gas, which has mass . If a significant fraction of the gas mass is removed, then the binding energy of the individual stars will decrease.

If the system is initially in virial equilibrium, then we have


(where is the kinetic energy and is the potential energy) and so


If the cluster experiences sudden mass loss due to instantaneous gas expulsion caused by the supernova, then the total cluster mass (stars + gas) is reduced by a factor of . The potential energy becomes and we can write the total energy as:




In order for the cluster to be unbound, the mean energy per star must be positive, which will happen if (because the potential energy, is negative). (This is analagous to the condition we derived for binary destruction by supernovae in the previous section.) Hence, for a virialised system where the stars and gas are well-mixed, supernovae can unbind the system if a significant fraction of the mass of the cluster (in this case the gas left over from star formation) is expelled.

One can also show that in the case of rapid mass loss, the cluster will expand by a factor


where and are the radii before and after instantaneous mass-loss Hills80 .

3.3 Predictions from simulations

If the binding energy of stars suddenly decreases due to a supernova, their velocities increase to the point that they are no longer bound to the centre of mass of the cluster. This has been shown in numerous theoretical and numerical works where the gas is modelled as a background potential Tutukov78 ; Lada84 ; Goodwin06 ; Baumgardt07 .

In -body simulations of star cluster evolution, it is possible to model the gas as a single potential, which is either static or is time-dependent (i.e. the potential can be made to decrease with time to mimic the depletion of gas). This technique implicitly assumes that the stars and gas are well-mixed but they do not interact apart from in a gravitational sense.

As -body simulations are computationally inexpensive, many calculations that include a background gas potential have been run by various authors. In general, they find that the effect of removing the gas is to unbind a significant fraction of the stellar mass. The action of unbinding this mass is to increase the velocity dispersion of the stars. A system where most of the stellar mass is not bound is said to be ‘supervirial’, and the system will have a higher velocity dispersion than the equivalent system that is in virial equilibrium.

The physical outcome of a supervirial system is for it to expand. When a star cluster expands the stars start to feel less of the cluster’s influence and more of the tidal field of the Galaxy. This so-called ‘Jacobi radius’ is given by


where is the Galactocentric distance, the mass of the Galaxy and is the mass of the cluster.

For a cluster of 1000 M at the position of the Sun in the Milky Way’s potential, the Jacobi radius at which the tidal field of the galaxy dominates over the cluster is at around 7 pc. Any star exceeding this radius is likely to escape the cluster and become a member of the Galactic disc.

Because this Jacobi radius is a function of cluster mass, the more members that are lost, the smaller the radius becomes. For this reason, clusters that are supervirial (expanding) are likely to dissolve more readily than those in virial equilibrium. The chances of a cluster surviving the gas removal phase therefore rest on the star formation efficiency (i.e. how much mass remains as gas) and the efficiency of supernovae to expel the gas and unbind the cluster.

3.4 Observational evidence?

Until recently, observations of star clusters appeared to corroborate this picture. The observed velocity dispersions of many clusters were consistent with the clusters being in a supervirial state when compared with an estimate for the dispersion for the cluster to be in virial equilibrium Bastian06 ; Olczak08 . The virial velocity dispersion is estimated using


where is the graviational constant and is the enclosed mass within a radius . is the structure parameter, and typically is 10 for a dynamically evolved, centrally concentrated cluster. If is significantly higher than the observed velocity dispersion, the region is said to be supervirial, or unbound.

However, estimates of the velocity dispersion in clusters had failed to account for the the orbital motion of binary stars. The fraction of stars that are in binary systems is around 50% in the solar neighbourhood, and is thought to be similar, if not higher in young star clusters. A substantial fraction of binaries have semimajor axes less than 1 au, so-called spectroscopic binary stars. The component stars of these binaries have extremely fast orbital motions, which dominate the velocity of the entire binary system with respect to its motion around the centre of mass of the cluster itself.

If the underlying velocity dispersion for the motion of stars can be approximated as a Gaussian, then the contribution of the orbital motion of binary stars both adds extra signal to the wings of the Gaussian, and widens it. The net effect is to make a measured velocity dispersion appear higher than it really is, and this gives the impression that the star cluster is undergoing supervirial expansion, when it is in fact in virial equilibrium Gieles10 .

3.5 The effects on a realistic potential

When the orbital motion of binary stars is accounted for, most star clusters are observed in virial equilibrium Cottaar12b . Several clusters also appear to be devoid of gas, suggesting that the star formation efficiency is either unusually high (and so the removal of gas has not affected the dynamics of the stars) or that the gas removal has not affected the dynamics.

Clues as to why removing a significant amount of the potential does not influence the subsequent dynamical evolution come from hydrodynamic simulations of the earliest stages of star formation. Typically, these simulations replace overdensities of gas with ‘sink particles’ in order to speed up the calculation. Depending on the resolution limit of the simulation, sink particles are used to model either entire star clusters, or single stars. They are assigned a radius, where any gas crossing this is removed from the simulation and added to the mass of the sink.

In simulations of star formation, several researchers Kruijssen12a ; Dale12b have shown that the locations where the sink particles form are almost entirely devoid of gas. The physical interpretation of this is that the local star formation efficiency is high, so most of the gas is converted to sinks, which are spatially arranged in groups. Therefore, the effects on the dynamics of the stars of instantaneously removing this gas due to supernovae are virtually negligible, because the potential of the remaining gas is so small. If this argument scales up to clusters of 100 – 1000 M, then it is possible that the supernovae explosions of massive stars have no effect on the dynamical evolution of the cluster. Furthermore, any cluster that is gas free at ages of less than 5Myr will never be unbound due to supernovae, because supernovae only start to explode at these ages.

In summary, it appears that the effects of supernovae on the dynamics of star clusters are probably minimal, although this debate is currently ongoing in the literature.

4 Conclusions

I have discussed the effects of rapid mass-loss experienced by a bound binary star system due to a core-collapse supernova and shown that this often leads to the disruption of the system. Furthermore, the secondary, less massive star is often ejected at high velocity – a so-called ‘runaway’ star. Because binary systems which contain stars massive enough to explode as supernovae have a higher likelihood of containing another massive star of comparable mass, the runaway stars produced by supernovae are readily observable. In some cases they can be traced back to their natal star-forming regions, which places constraints on the birth environment of massive star binaries and the initial density of star formation.

The effects of supernovae on the dynamical evolution of star clusters are more modest as the total stellar mass of the star cluster does not significantly decrease during the supernova(e). However, because star formation is so inefficient, a significant fraction of the mass in a star cluster is thought to be in the form of gas left over from the star formation process. Many authors have argued that supernovae could remove this gas and therefore unbind the stars in the cluster, leading to the dissolution of the cluster. Recent observational and theoretical work has shown that the effects of this gas removal on the dynamical evolution of a star cluster are likely to be minimal.

5 Acknowledgements

I am grateful to John Bray for pointing out several typos in equations 24 and 25, which appeared in the original version of this manuscript. I acknowledge support from the Royal Astronomical Society in the form of a research fellowship.


  • (1) S.P. Goodwin, P. Kroupa, A. Goodman, A. Burkert, in Protostars and Planets V, ed. by B. Reipurth, D. Jewitt, K. Keil (2007), pp. 133–147
  • (2) N. Moeckel, M.R. Bate, MNRAS 404, 721 (2010)
  • (3) C.J. Lada, E.A. Lada, ARA&A 41, 57 (2003)
  • (4) J.M.D. Kruijssen, MNRAS 426, 3008 (2012)
  • (5) N. Bastian, K.R. Covey, M.R. Meyer, ARA&A 48, 339 (2010)
  • (6) F. Zwicky, Handbuch der Physik 51, 766 (1957)
  • (7) A. Blaauw, BAIN 15, 265 (1961)
  • (8) J. Boersma, BAIN 15, 291 (1961)
  • (9) P.J.T. Leonard, M.J. Duncan, AJ 99, 608 (1990). DOI 10.1086/115354
  • (10) R.J. Allison, S.P. Goodwin, MNRAS 415, 1967 (2011)
  • (11) S. Oh, P. Kroupa, J. Pflamm-Altenburg, ApJ 805, 92 (2015). DOI 10.1088/0004-637X/805/2/92
  • (12) P. Kroupa, A. Aarseth, J. Hurley, MNRAS 321, 699 (2001)
  • (13) N.J. Wright, R.J. Parker, S.P. Goodwin, J.J. Drake, MNRAS 438, 639 (2014)
  • (14) S.S. Huang, AJ 61, 49 (1956). DOI 10.1086/107290
  • (15) G. Meynet, A. Maeder, A&A 404, 975 (2003). DOI 10.1051/0004-6361:20030512
  • (16) T.M. Tauris, R.J. Takens, A&A 330, 1047 (1998)
  • (17) S.F. Portegies Zwart, ApJ 544, 437 (2000). DOI 10.1086/317190
  • (18) D. Raghavan, H.A. McMaster, T.J. Henry, D.W. Latham, G.W. Marcy, B.D. Mason, D.R. Gies, R.J. White, T.A. ten Brummelaar, ApJSS 190, 1 (2010)
  • (19) M.M. Reggiani, M.R. Meyer, ApJ 738, 60 (2011). DOI 10.1088/0004-637X/738/1/60
  • (20) B.D. Mason, W.I. Hartkopf, D.R. Gies, T.J. Henry, J.W. Helsel, AJ 137, 3358 (2009)
  • (21) H. Sana, A. de Koter, S.E. de Mink, P.R. Dunstall, C.J. Evans, et al., A&A 550, A107 (2013). DOI 10.1051/0004-6361/201219621
  • (22) E. Öpik, Tartu Obs. Publ. 25, 6 (1924)
  • (23) P. Kroupa, MNRAS 277, 1491 (1995)
  • (24) A.J. Horton, M.R. Bate, I.A. Bonnell, MNRAS 321, 585 (2001). DOI 10.1046/j.1365-8711.2001.04071.x
  • (25) R.R. King, R.J. Parker, J. Patience, S.P. Goodwin, MNRAS 421, 2025 (2012)
  • (26) A. Cartwright, A.P. Whitworth, MNRAS 348, 589 (2004)
  • (27) R.J. Parker, MNRAS 445, 4037 (2014)
  • (28) A. Blaauw, ARA&A 2, 213 (1964). DOI 10.1146/annurev.aa.02.090164.001241
  • (29) Y.N. Efremov, AJ 110, 2757 (1995). DOI 10.1086/117728
  • (30) B.G. Elmegreen, D.M. Elmegreen, AJ 121, 1507 (2001)
  • (31) A. Pellerin, M. Meyer, J. Harris, D. Calzetti, ApJL 658, L87 (2007). DOI 10.1086/515437
  • (32) E.E. Salpeter, ApJ 121, 161 (1955)
  • (33) B.G. Elmegreen, ApJ 648, 572 (2006)
  • (34) C.D. Matzner, ApJ 566, 302 (2002). DOI 10.1086/338030
  • (35) M.R. Krumholz, C.D. Matzner, C.F. McKee, ApJ 653, 361 (2006). DOI 10.1086/508679
  • (36) J.G. Hills, ApJ 235, 986 (1980)
  • (37) A.V. Tutukov, A&A 70, 57 (1978)
  • (38) C.J. Lada, M. Margulis, D. Dearborn, ApJ 285, 141 (1984)
  • (39) S.P. Goodwin, N. Bastian, MNRAS 373, 752 (2006)
  • (40) H. Baumgardt, P. Kroupa, MNRAS 380, 1589 (2007). DOI 10.1111/j.1365-2966.2007.12209.x
  • (41) N. Bastian, S.P. Goodwin, MNRAS 369, L9 (2006)
  • (42) C. Olczak, S. Pfalzner, A. Eckart, A&A 488, 191 (2008)
  • (43) M. Gieles, H. Sana, S.F. Portegies Zwart, MNRAS 402, 1750 (2010). DOI 10.1111/j.1365-2966.2009.15993.x
  • (44) M. Cottaar, M.R. Meyer, R.J. Parker, A&A 547, A35 (2012). DOI 10.1051/0004-6361/201219673
  • (45) J.M.D. Kruijssen, T. Maschberger, N. Moeckel, C.J. Clarke, N. Bastian, I.A. Bonnell, MNRAS 419, 841 (2012)
  • (46) J.E. Dale, B. Ercolano, I.A. Bonnell, MNRAS 424, 377 (2012). DOI 10.1111/j.1365-2966.2012.21205.x


Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description