-
Massive star cluster formation III. Early mass segregation during cluster assembly
Authors:
Brooke Polak,
Mordecai-Mark Mac Low,
Ralf S. Klessen,
Simon Portegies Zwart,
Eric P. Andersson,
Sabrina M. Appel,
Claude Cournoyer Cloutier,
Simon C. O. Glover,
Stephen L. W. McMillan
Abstract:
Mass segregation is seen in many star clusters, but whether massive stars form in the center of a cluster or migrate there dynamically is still debated. N-body simulations have shown that early dynamical mass segregation is possible when sub-clusters merge to form a dense core with a small crossing time. However, the effect of gas dynamics on both the formation and dynamics of the stars could inhi…
▽ More
Mass segregation is seen in many star clusters, but whether massive stars form in the center of a cluster or migrate there dynamically is still debated. N-body simulations have shown that early dynamical mass segregation is possible when sub-clusters merge to form a dense core with a small crossing time. However, the effect of gas dynamics on both the formation and dynamics of the stars could inhibit the formation of the dense core. We aim to study the dynamical mass segregation of star cluster models that include gas dynamics and self-consistently form stars from the dense substructure in the gas. Our models use the Torch framework, which is based on AMUSE and includes stellar and magnetized gas dynamics, as well as stellar evolution and feedback from radiation, stellar winds, and supernovae. Our models consist of three star clusters forming from initial turbulent spherical clouds of mass $10^{4,5,6}\rm~M_\odot$ and radius $11.7\rm~pc$ that have final stellar masses of $3.6\times10^3\rm~M_\odot$, $6.5\times10^4\rm~M_\odot$, and $8.9\times10^5\rm~M_\odot$, respectively. There is no primordial mass segregation in the model by construction. All three clusters become dynamically mass segregated at early times via collapse confirming that this mechanism occurs within sub-clusters forming directly out of the dense substructure in the gas. The dynamics of the embedded gas and stellar feedback do not inhibit the collapse of the cluster. We find that each model cluster becomes mass segregated within $2~$Myr of the onset of star formation, reaching the levels observed in young clusters in the Milky Way. However, we note that the exact values are highly time-variable during these early phases of evolution. Massive stars that segregate to the center during core collapse are likely to be dynamically ejected, a process that can decrease the overall level of mass segregation again.
△ Less
Submitted 26 August, 2024;
originally announced August 2024.
-
On the relation between magnetic field strength and gas density in the interstellar medium: A multiscale analysis
Authors:
David J. Whitworth,
Sundar Srinivasan,
Ralph E. Pudritz,
Mordecai M. Mac Low,
Rowan J. Smith,
Aina Palau,
Kate Pattle,
Gwendoline Eadie,
Hector Robinson,
Rachel Pillsworth,
James Wadsley,
Noe Brucy,
Ugo Lebreuilly,
Patrick Hennebelle,
Philipp Girichidis,
Fred A. Gent,
Jessy Marin,
Lylon Sánchez Valido,
Vianey Camacho,
Ralf S. Klessen,
Enrique Vázquez-Semadeni
Abstract:
The relation between magnetic field strength B and gas density n in the interstellar medium is of fundamental importance to many areas of astrophysics, from protostellar disks to galaxy evolution. We present and compare Bayesian analyses of the B - n relation for a comprehensive observational data set, as well as a large body of numerical MHD simulations. We extend the original Zeeman relation of…
▽ More
The relation between magnetic field strength B and gas density n in the interstellar medium is of fundamental importance to many areas of astrophysics, from protostellar disks to galaxy evolution. We present and compare Bayesian analyses of the B - n relation for a comprehensive observational data set, as well as a large body of numerical MHD simulations. We extend the original Zeeman relation of Crutcher et al. (2010) with a large body of magnetic data that includes 700 observations with the Davis-Chandrasekhar-Fermi method. By using a new multiparameter Bayesian analysis we present a new, more general, time-averaged observational relation: B \propto n^{0.27 \pm 0.017} for n \leq n_0 and B \propto n^{0.54 \pm 0.18} for n \geq n_0, with n_0 = 924^(+145-144) cm^-3. We perform a separate analysis on 19 numerical magnetohydrodynamics simulations that cover a wide range of scales, resolutions, initial conditions, and completed with a variety of codes: arepo, flash, pencil, and ramses. The power law exponents derived from the simulations depend on several physical factors including: dynamo effects, time scales, turbulence, and the initial seed field strength. In particular, early-time simulations where the density, velocity and magnetic fields are unevolved do not match the observational scalings. The simulations that trace the observed density range best, the evolved dwarf galaxy and Milky Way like galaxy simulations, settle into a near consistent exponent of = 0.5 in the dense gas, with variability in the diffuse gas exponent.
△ Less
Submitted 25 July, 2024;
originally announced July 2024.
-
Massive Star Cluster Formation II. Runaway Stars as Fossils of Sub-Cluster Mergers
Authors:
Brooke Polak,
Mordecai-Mark Mac Low,
Ralf S. Klessen,
Simon Portegies Zwart,
Eric P. Andersson,
Sabrina M. Appel,
Claude Cournoyer-Cloutier,
Simon C. O. Glover,
Stephen L. W. McMillan
Abstract:
Two main mechanisms have classically been proposed for the formation of runaway stars. The binary supernova scenario (BSS) suggests that a massive star in a binary explodes as a supernova, ejecting its companion. The dynamical ejection scenario suggests that a star is ejected through a strong dynamical encounter between multiple stars. We propose a third mechanism for the formation of runaway star…
▽ More
Two main mechanisms have classically been proposed for the formation of runaway stars. The binary supernova scenario (BSS) suggests that a massive star in a binary explodes as a supernova, ejecting its companion. The dynamical ejection scenario suggests that a star is ejected through a strong dynamical encounter between multiple stars. We propose a third mechanism for the formation of runaway stars: the sub-cluster ejection scenario (SCES), where an infalling sub-cluster of stars is ejected out of the cluster by a slingshot interaction with the contracting gravitational potential of the assembling cluster. We demonstrate the SCES in a star-by-star simulation of a young massive cluster forming from a $10^6\rm~M_\odot$ gas cloud using the Torch framework. This star cluster forms hierarchically through a sequence of sub-cluster mergers, determined by the initial turbulent, spherical initial conditions of the gas. We find that these mergers drive the formation of runaway stars in our model. Late-forming sub-clusters fall into the central potential, where they are ejected on a slingshot trajectory, forming groups of runaway stars that are distributed highly anisotropically. Runaways formed by the same SCES share similar ages, velocities, and ejection directions. Surveying observations, we identify several SCES candidate groups with anisotropic ejection directions. The SCES is capable of producing runaway binaries: two wide dynamical binaries in infalling sub-clusters were tightened through ejection. This allows for another velocity kick via subsequent BSS, which is a promising avenue for producing hypervelocity stars unbound to the Galaxy. The SCES occurs when sub-cluster formation is resolved. We expect non-spherical initial gas distributions to increase runaway star numbers. The observation of groups of runaway stars formed via SCES thus reveals the assembly history of their natal clusters.
△ Less
Submitted 20 May, 2024;
originally announced May 2024.
-
Computational approaches to modeling dynamos in galaxies
Authors:
Maarit J. Korpi-Lagg,
Mordecai-Mark Mac Low,
Frederick A. Gent
Abstract:
Galaxies are observed to host magnetic fields with a typical total strength of around 15microgauss. A coherent large-scale field constitutes up to a few microgauss of the total, while the rest is built from strong magnetic fluctuations over a wide range of spatial scales. This represents sufficient magnetic energy for it to be dynamically significant. Several questions immediately arise: What is t…
▽ More
Galaxies are observed to host magnetic fields with a typical total strength of around 15microgauss. A coherent large-scale field constitutes up to a few microgauss of the total, while the rest is built from strong magnetic fluctuations over a wide range of spatial scales. This represents sufficient magnetic energy for it to be dynamically significant. Several questions immediately arise: What is the physical mechanism that gives rise to such magnetic fields? How do these magnetic fields affect the formation and evolution of galaxies? In which physical processes do magnetic fields play a role, and how can that role be characterized? Numerical modelling of magnetized flows in galaxies is playing an ever-increasing role in finding those answers. We review major techniques used for these models. Current results strongly support the conclusion that field growth occurs during the formation of the first galaxies on timescales shorter than their accretion timescales due to small-scale turbulent dynamos. The saturated small-scale dynamo maintains field strengths at a few percent of equipartition with turbulence. The subsequent action of large-scale dynamos in differentially rotating discs produces observed modern field strengths in equipartition with the turbulence and having power at large scales. The field structure resulting appears consistent with observations including Faraday rotation and polarisation from synchrotron and dust thermal emission. Major remaining challenges include scaling numerical models toward realistic scale separations and Prandtl and Reynolds numbers.
△ Less
Submitted 8 January, 2024;
originally announced January 2024.
-
Massive Star Cluster Formation I. High Star Formation Efficiency While Resolving Feedback of Individual Stars
Authors:
Brooke Polak,
Mordecai-Mark Mac Low,
Ralf S. Klessen,
Jia Wei Teh,
Claude Cournoyer-Cloutier,
Eric P. Andersson,
Sabrina M. Appel,
Aaron Tran,
Sean C. Lewis,
Maite J. C. Wilhelm,
Simon Portegies Zwart,
Simon C. O. Glover,
Long Wang,
Stephen L. W. McMillan
Abstract:
The mode of star formation that results in the formation of globular clusters and young massive clusters is difficult to constrain through observations. We present models of massive star cluster formation using the Torch framework, which uses AMUSE to couple distinct multi-physics codes that handle star formation, stellar evolution and dynamics, radiative transfer, and magnetohydrodynamics. We upg…
▽ More
The mode of star formation that results in the formation of globular clusters and young massive clusters is difficult to constrain through observations. We present models of massive star cluster formation using the Torch framework, which uses AMUSE to couple distinct multi-physics codes that handle star formation, stellar evolution and dynamics, radiative transfer, and magnetohydrodynamics. We upgrade Torch by implementing the N-body code PeTar, thereby enabling Torch to handle massive clusters forming from $10^6\rm\, M_\odot$ clouds with $\ge10^5$ individual stars. We present results from Torch simulations of star clusters forming from $10^4, 10^5$, and $10^6\rm M_\odot$ turbulent, spherical gas clouds (named M4, M5, M6) of radius $R=11.7$ pc. We find that star formation is highly efficient and becomes more so at higher cloud mass and surface density. For M4, M5, and M6 with initial surface densities $2.325\times 10^{1,2,3}\rm\, M_\odot\, pc^{-2}$, after a free-fall time of $t_{ff}=6.7,2.1,0.67$ Myr, we find that $\sim30\%$, 40%, and 60% of the cloud mass has formed into stars, respectively. The final integrated star formation efficiency is $32\%,\, 65\%$, and 85\% for M4, M5, and M6. Observations of nearby clusters similar to M4 have similar integrated star formation efficiencies of $\leq30\%$. The M5 and M6 models represent a different regime of cluster formation that is more appropriate for the conditions in starburst galaxies and gas-rich galaxies at high redshift, and that leads to a significantly higher efficiency of star formation. We argue that young massive clusters build up through short efficient bursts of star formation in regions that are sufficiently dense ($\ge 10^2 \rm\,M_\odot\,pc^{-2}$) and massive ($\ge10^5\rm\, M_\odot$). In such environments, the dynamical time of the cloud becomes short enough that stellar feedback cannot act quickly enough to slow star formation.
△ Less
Submitted 11 December, 2023;
originally announced December 2023.
-
Pre-supernova feedback sets the star cluster mass function to a power law and reduces the cluster formation efficiency
Authors:
Eric P. Andersson,
Mordecai-Mark Mac Low,
Oscar Agertz,
Florent Renaud,
Hui Li
Abstract:
The star cluster initial mass function is observed to have an inverse power law exponent around 2, yet there is no consensus on what determines this distribution, and why some variation is observed in different galaxies. Furthermore, the cluster formation efficiency covers a range of values, particularly when considering different environments. These clusters are often used to empirically constrai…
▽ More
The star cluster initial mass function is observed to have an inverse power law exponent around 2, yet there is no consensus on what determines this distribution, and why some variation is observed in different galaxies. Furthermore, the cluster formation efficiency covers a range of values, particularly when considering different environments. These clusters are often used to empirically constrain star formation and as fundamental units for stellar feedback models. Detailed galaxy models must therefore accurately capture the basic properties of observed clusters to be considered predictive. We use hydrodynamical simulations of a dwarf galaxy as a laboratory to study star cluster formation. We test different combinations of stellar feedback mechanisms, including stellar winds, ionizing radiation, and supernovae. Each feedback mechanism affects the cluster formation efficiency and cluster mass function. Increasing the feedback budget by combining the different types of feedback decreases the cluster formation efficiency by reducing the number of massive clusters. Ionizing radiation is found to be especially influential. This effect depends on the timing of feedback initiation, as shown by comparing early and late feedback. Early feedback occurs from ionizing radiation and stellar winds with onset immediately after a massive star is formed. Late feedback occurs when energy injection only starts after the main-sequence lifetime of the most massive SN progenitor, a timing that is further influenced by the choice of the most massive SN progenitor. Late feedback alone results in a broad, flat mass function, approaching a log-normal shape in the complete absence of feedback. Early feedback, on the other hand, produces a power-law cluster mass function with lower formation efficiency, albeit with a steeper slope than that usually observed.
△ Less
Submitted 23 August, 2023;
originally announced August 2023.
-
Transition from small-scale to large-scale dynamo in a supernova-driven, multiphase medium
Authors:
Frederick A. Gent,
Mordecai-Mark Mac Low,
Maarit J. Korpi-Lagg
Abstract:
Magnetic fields are widely recognised as critical at many scales to galactic dynamics and structure, including multiphase pressure balance, dust processing, and star formation. Using imposed magnetic fields cannot reliably model the interstellar medium's (ISM) dynamical structure nor phase interactions. Dynamos must be modelled. ISM models exist of turbulent magnetic fields using small-scale dynam…
▽ More
Magnetic fields are widely recognised as critical at many scales to galactic dynamics and structure, including multiphase pressure balance, dust processing, and star formation. Using imposed magnetic fields cannot reliably model the interstellar medium's (ISM) dynamical structure nor phase interactions. Dynamos must be modelled. ISM models exist of turbulent magnetic fields using small-scale dynamo (SSD). Others model the large-scale dynamo (LSD) organising magnetic fields at scale of the disc or spiral arms. Separately, neither can fully describe the galactic magnetic field dynamics nor topology. We model the LSD and SSD together at sufficient resolution to use the low explicit Lagrangian resistivity required. The galactic SSD saturates within 20 Myr. We show that the SSD is quite insensitive to the presence of an LSD and is even stronger in the presence of a large-scale shear flow. The LSD grows more slowly in the presence of SSD, saturating after 5 Gyr vs. 1--2 Gyr in studies where the SSD is weak or absent. The LSD primarily grows in warm gas in the galactic midplane. Saturation of the LSD occurs due to $α$-quenching near the midplane as the growing mean field produces a magnetic $α$ that opposes the kinetic $α$. The magnetic energy in our models of the LSD shows slightly sublinear response to increasing resolution, indicating that we are converging towards the physical solution at 1 pc resolution. Clustering supernovae in OB associations increases the growth rates for both the SSD and the LSD, compared to a horizontally uniform supernova distribution.
△ Less
Submitted 20 November, 2023; v1 submitted 12 June, 2023;
originally announced June 2023.
-
Dust dynamics in current sheets within protoplanetary disks. I. Isothermal models including ambipolar diffusion and Ohmic resistivity
Authors:
U. Lebreuilly,
M. -M. Mac Low,
B. Commerçon,
D. S. Ebel
Abstract:
Context: Chondrules originate from the reprocessing of dust grains. They are key building blocks of telluric planets, yet their formation, which must happen in strongly localized regions of high temperature, remains poorly understood. Aims: We examine the dust spatial distribution near regions of strong local heating produced by current sheets, as a step toward exploring a potential path for chond…
▽ More
Context: Chondrules originate from the reprocessing of dust grains. They are key building blocks of telluric planets, yet their formation, which must happen in strongly localized regions of high temperature, remains poorly understood. Aims: We examine the dust spatial distribution near regions of strong local heating produced by current sheets, as a step toward exploring a potential path for chondrule formation. We further aim to investigate current sheet formation under various conditions in protoplanetary disks in the presence of ambipolar diffusion and Ohmic resistivity and the effect of current sheet morphology on dust dynamics in their vicinity. Methods: We use the RAMSES code including modules for non-ideal magnetohydrodynamics and solution of the dynamics of multiple sizes of dust grains to compute unstratified shearing box simulations of current sheet formation. We investigate, through seven models the effect of the ambipolar diffusion and Ohmic resistivity strength, the initial density, and magnetic field, as well as the resolution and box size. Results: We find that current sheets form in all our models, with typical widths of 0.001-0.01 AU and that strong dust fraction variations occur for millimeter-sized grains. These variations are typically of an order of magnitude and up to two orders of magnitude for the most favorable cases. We also show that the box size and resolution has a strong impact on the current sheet distribution and intensity. Conclusions: The formation of current sheets that can intensely heat their surroundings near strong dynamical dust fraction variations could have important implications for chondrule formation, as it appears likely to happen in regions of large dust fraction.
△ Less
Submitted 26 May, 2023;
originally announced May 2023.
-
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.
-
Radiation shielding of protoplanetary discs in young star-forming regions
Authors:
Martijn J. C. Wilhelm,
Simon Portegies Zwart,
Claude Cournoyer-Cloutier,
Sean C. Lewis,
Brooke Polak,
Aaron Tran,
Mordecai-Mark Mac Low
Abstract:
Protoplanetary discs spend their lives in the dense environment of a star forming region. While there, they can be affected by nearby stars through external photoevaporation and dynamic truncations. We present simulations that use the AMUSE framework to couple the Torch model for star cluster formation from a molecular cloud with a model for the evolution of protoplanetary discs under these two en…
▽ More
Protoplanetary discs spend their lives in the dense environment of a star forming region. While there, they can be affected by nearby stars through external photoevaporation and dynamic truncations. We present simulations that use the AMUSE framework to couple the Torch model for star cluster formation from a molecular cloud with a model for the evolution of protoplanetary discs under these two environmental processes. We compare simulations with and without extinction of photoevaporation-driving radiation. We find that the majority of discs in our simulations are considerably shielded from photoevaporation-driving radiation for at least 0.5 Myr after the formation of the first massive stars. Radiation shielding increases disc lifetimes by an order of magnitude and can let a disc retain more solid material for planet formation. The reduction in external photoevaporation leaves discs larger and more easily dynamically truncated, although external photoevaporation remains the dominant mass loss process. Finally, we find that the correlation between disc mass and projected distance to the most massive nearby star (often interpreted as a sign of external photoevaporation) can be erased by the presence of less massive stars that dominate their local radiation field. Overall, we find that the presence and dynamics of gas in embedded clusters with massive stars is important for the evolution of protoplanetary discs.
△ Less
Submitted 7 February, 2023;
originally announced February 2023.
-
On the distribution of the Cold Neutral Medium in galaxy discs
Authors:
Rowan J. Smith,
Robin Tress,
Juan D. Soler,
Ralf S. Klessen,
Simon C. O. Glover,
Patrick Hennebelle,
Sergio Molinari,
Mordecai-Mark Mac Low,
David Whitworth
Abstract:
The Cold Neutral Medium (CNM) is an important part of the galactic gas cycle and a precondition for the formation of molecular and star forming gas, yet its distribution is still not fully understood. In this work we present extremely high resolution simulations of spiral galaxies with time-dependent chemistry such that we can track the formation of the CNM, its distribution within the galaxy, and…
▽ More
The Cold Neutral Medium (CNM) is an important part of the galactic gas cycle and a precondition for the formation of molecular and star forming gas, yet its distribution is still not fully understood. In this work we present extremely high resolution simulations of spiral galaxies with time-dependent chemistry such that we can track the formation of the CNM, its distribution within the galaxy, and its correlation with star formation. We find no strong radial dependence between the CNM fraction and total HI due to the decreasing interstellar radiation field counterbalancing the decreasing gas column density at larger galactic radii.However, the CNM fraction does increase in spiral arms where the CNM distribution is clumpy, rather than continuous, overlapping more closely with H2. The CNM doesn't extend out radially as far as HI, and the vertical scale height is smaller in the outer galaxy compared to HI with no flaring. The CNM column density scales with total midplane pressure and disappears from the gas phase below values of PT/kB =1000 K/cm3. We find that the star formation rate density follows a similar scaling law with CNM column density to the total gas Kennicutt-Schmidt law. In the outer galaxy we produce realistic vertical velocity dispersions in the HI purely from galactic dynamics but our models do not predict CNM at the extremely large radii observed in HI absorption studies of the Milky Way. We suggest that extended spiral arms might produce isolated clumps of CNM at these radii.
△ Less
Submitted 19 July, 2023; v1 submitted 19 January, 2023;
originally announced January 2023.
-
Fast methods for tracking grain coagulation and ionization. III. Protostellar collapse with non-ideal MHD
Authors:
Pierre Marchand,
Ugo Lebreuilly,
Mordecai-Mark Mac Low,
Vincent Guillet
Abstract:
Dust grains influence many aspects of star formation, including planet formation, opacities for radiative transfer, chemistry, and the magnetic field via Ohmic, Hall, and ambipolar diffusion. The size distribution of the dust grains is the primary characteristic influencing all these aspects. Grain size increases by coagulation throughout the star formation process. We describe here numerical simu…
▽ More
Dust grains influence many aspects of star formation, including planet formation, opacities for radiative transfer, chemistry, and the magnetic field via Ohmic, Hall, and ambipolar diffusion. The size distribution of the dust grains is the primary characteristic influencing all these aspects. Grain size increases by coagulation throughout the star formation process. We describe here numerical simulations of protostellar collapse using methods described in earlier papers of this series. We compute the evolution of the grain size distribution from coagulation and the non-ideal magnetohydrodynamics effects self-consistently and at low numerical cost. We find that the coagulation efficiency is mostly affected by the time spent in high-density regions. Starting from sub-micron radii, grain sizes reach more than 100 μm in an inner protoplanetary disk that is only 1000 years old. We also show that the growth of grains significantly affects the resistivities, and indirectly the dynamics and angular momentum of the disk.
△ Less
Submitted 9 January, 2023; v1 submitted 4 January, 2023;
originally announced January 2023.
-
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.
-
Magnetic fields do not suppress global star formation in low metallicity dwarf galaxies
Authors:
David J. Whitworth,
Rowan J. Smith,
Ralf S. Klessen,
Mordecai-Mark Mac Low,
Simon C. O. Glover,
Robin Tress,
Rudiger Pakmor,
Juan D. Soler
Abstract:
Many studies concluded that magnetic fields suppress star formation in molecular clouds and Milky Way like galaxies. However, most of these studies are based on fully developed fields that have reached the saturation level, with little work on investigating how an initial weak primordial field affects star formation in low metallicity environments. In this paper, we investigate the impact of a wea…
▽ More
Many studies concluded that magnetic fields suppress star formation in molecular clouds and Milky Way like galaxies. However, most of these studies are based on fully developed fields that have reached the saturation level, with little work on investigating how an initial weak primordial field affects star formation in low metallicity environments. In this paper, we investigate the impact of a weak initial field on low metallicity dwarf galaxies. We perform high-resolution AREPO simulations of five isolated dwarf galaxies. Two models are hydrodynamical, two start with a primordial magnetic field of 10$^{-6} μ$G and different sub-solar metallicities, and one starts with a saturated field of 10$^{-2} μ$G. All models include a non-equilibrium, time-dependent chemical network that includes the effects of gas shielding from the ambient ultraviolet field. Sink particles form directly from the gravitational collapse of gas and are treated as star-forming clumps that can accrete gas. We vary the ambient uniform far ultraviolet field, and cosmic ray ionization rate between 1\% and 10\% of solar values. We find that the magnetic field has little impact on the global star formation rate, which is in tension with some previously published results. We further find that the initial field strength has little impact on the global star formation rate. We show that an increase in the mass fractions of both molecular hydrogen and cold gas, along with changes in the perpendicular gas velocity dispersion and the magnetic field acting in the weak-field model, overcome the expected suppression in star formation.
△ Less
Submitted 12 January, 2023; v1 submitted 10 October, 2022;
originally announced October 2022.
-
The small-scale dynamo in a multiphase supernova-driven medium
Authors:
Frederick A. Gent,
Mordecai-Mark Mac Low,
Maarit J. Korpi-Lagg,
Nishant K. Singh
Abstract:
Magnetic fields grow quickly, even at early cosmological times, suggesting the action of a small-scale dynamo (SSD) in the interstellar medium (ISM) of galaxies. Many studies have focused on idealized, isotropic, homogeneous, turbulent driving of the SSD. Here we analyze more realistic simulations of supernova-driven turbulence to understand how it drives an SSD. We find that SSD growth rates are…
▽ More
Magnetic fields grow quickly, even at early cosmological times, suggesting the action of a small-scale dynamo (SSD) in the interstellar medium (ISM) of galaxies. Many studies have focused on idealized, isotropic, homogeneous, turbulent driving of the SSD. Here we analyze more realistic simulations of supernova-driven turbulence to understand how it drives an SSD. We find that SSD growth rates are intermittently variable as a result of the evolving multiphase ISM structure. Rapid growth in the magnetic field typically occurs in hot gas, with the highest overall growth rates occurring when the fractional volume of hot gas is large. SSD growth rates correlate most strongly with vorticity and fluid Reynolds number, which also both correlate strongly with gas temperature. Rotational energy exceeds irrotational energy in all phases, but particularly in the hot phase while SSD growth is most rapid. Supernova (SN) rate does not significantly affect the ISM average kinetic energy density. Rather, higher temperatures associated with high SN rates tend to increase SSD growth rates. SSD saturates with total magnetic energy density around 5% of equipartition to kinetic energy density, increasing slightly with magnetic Prandtl number. While magnetic energy density in the hot gas can exceed that of the other phases when SSD grows most rapidly, it saturates below 5% of equipartition with kinetic energy in the hot gas, while in the cold gas it attains 100%. Fast, intermittent growth of the magnetic field appears to be a characteristic behavior of SN-driven, multiphase turbulence.
△ Less
Submitted 6 December, 2022; v1 submitted 10 October, 2022;
originally announced October 2022.
-
Fast methods to track grain coagulation and ionization. II. Extension to thermal ionization
Authors:
Pierre Marchand,
Vincent Guillet,
Ugo Lebreuilly,
Mordecai-Mark Mac Low
Abstract:
Thermal ionization is a critical process at temperatures T > 10 3 K, particularly during star formation. An increase in ionization leads to a decrease in nonideal magnetohydrodynamics (MHD) resistivities, which has a significant impact on protoplanetary disks and protostar formation. We developed an extension of the fast computational ionization method presented in our recent paper to include ther…
▽ More
Thermal ionization is a critical process at temperatures T > 10 3 K, particularly during star formation. An increase in ionization leads to a decrease in nonideal magnetohydrodynamics (MHD) resistivities, which has a significant impact on protoplanetary disks and protostar formation. We developed an extension of the fast computational ionization method presented in our recent paper to include thermal ionization. The model can be used to inexpensively calculate the density of ions and electrons and the electric charge of each size of grains for an arbitrary size distribution. This tool should be particularly useful for the self-consistent calculation of nonideal MHD resistivities in multidimensional simulations, especially of protostellar collapse and protoplanetary disks.
△ Less
Submitted 30 June, 2022; v1 submitted 23 February, 2022;
originally announced February 2022.
-
Anisotropic Infall and Substructure formation in Embedded Disks
Authors:
Aleksandra Kuznetsova,
Jaehan Bae,
Lee Hartmann,
Mordecai-Mark Mac Low
Abstract:
The filamentary nature of accretion streams found around embedded sources suggest that protostellar disks experience heterogenous infall from the star-forming environment, consistent with the accretion behavior onto star-forming cores in top-down star-cluster formation simulations. This may produce disk substructures in the form of rings, gaps, and spirals continuing to be identified by high-resol…
▽ More
The filamentary nature of accretion streams found around embedded sources suggest that protostellar disks experience heterogenous infall from the star-forming environment, consistent with the accretion behavior onto star-forming cores in top-down star-cluster formation simulations. This may produce disk substructures in the form of rings, gaps, and spirals continuing to be identified by high-resolution imaging surveys in both embedded Class 0/I and later Class II sources. We present a parameter study of anisotropic infall, informed by the properties of accretion flows onto protostellar cores in numerical simulations, and varying the relative specific angular momentum of incoming flows as well as their flow geometry. Our results show that anisotropic infall perturbs the disk and readily launches the Rossby wave instability (RWI). It forms vortices at the inner and outer edge of the infall zone where material is deposited. These vortices drive spiral waves and angular momentum transport, with some models able to drive stresses corresponding to a viscosity parameter on the order of $α\sim 10^{-2}$. The resulting azimuthal shear forms robust pressure bumps that act as barriers to radial drift of dust grains, as demonstrated by post-processing calculations of drift-dominated dust evolution. We discuss how a self-consistent model of anisotropic infall can account for the formation of millimeter rings in the outer disk as well as producing compact dust disks, consistent with observations of embedded sources.
△ Less
Submitted 10 February, 2022;
originally announced February 2022.
-
Symmetry Breaking in Dynamical Encounters in the Disks of Active Galactic Nuclei
Authors:
Yihan Wang,
Barry McKernan,
Saavik Ford,
Rosalba Perna,
Nathan W. C. Leigh,
Mordecai-Mark Mac Low
Abstract:
Active galactic nucleus (AGN) disks may be important sites of binary black hole (BBH) mergers. Here we show via numerical experiments with the high-accuracy, high precision code {\tt SpaceHub} that broken symmetry in dynamical encounters in AGN disks can lead to an asymmetry between prograde and retrograde BBH mergers. The direction of the hardening asymmetry depends on the initial binary semi-maj…
▽ More
Active galactic nucleus (AGN) disks may be important sites of binary black hole (BBH) mergers. Here we show via numerical experiments with the high-accuracy, high precision code {\tt SpaceHub} that broken symmetry in dynamical encounters in AGN disks can lead to an asymmetry between prograde and retrograde BBH mergers. The direction of the hardening asymmetry depends on the initial binary semi-major axis. An asymmetric distribution of mass-weighted projected spin $χ_{\rm eff}$ should therefore be expected in LIGO-Virgo detections of BBH mergers from AGN disks. This channel further predicts that negative $χ_{\rm eff}$ BBH mergers are most likely for massive binaries.
△ Less
Submitted 7 October, 2021;
originally announced October 2021.
-
Gravity Versus Magnetic Fields in Forming Molecular Clouds
Authors:
Juan C. Ibáñez-Mejía,
Mordecai-Mark Mac Low,
Ralf S. Klessen
Abstract:
Magnetic fields are dynamically important in the diffuse interstellar medium. Understanding how gravitationally bound, star-forming clouds form requires modeling of the fields in a self-consistent, supernova-driven, turbulent, magnetized, stratified disk. We employ the FLASH magnetohydrodynamics code to follow the formation and early evolution of clouds with final masses of 3-8…
▽ More
Magnetic fields are dynamically important in the diffuse interstellar medium. Understanding how gravitationally bound, star-forming clouds form requires modeling of the fields in a self-consistent, supernova-driven, turbulent, magnetized, stratified disk. We employ the FLASH magnetohydrodynamics code to follow the formation and early evolution of clouds with final masses of 3-8 $\times 10^3 M_{\odot}$ within such a simulation. We use the code's adaptive mesh refinement capabilities to concentrate numerical resolution in zoom-in regions covering single clouds, allowing us to investigate the detailed dynamics and field structure of individual self-gravitating clouds in a consistent background medium. Our goal is to test the hypothesis that dense clouds are dynamically evolving objects far from magnetohydrostatic equilibrium. We find that the cloud envelopes are magnetically supported with field lines parallel to density gradients and flow velocity, as indicated by the histogram of relative orientations and other statistical measures. In contrast, the dense cores of the clouds are gravitationally dominated, with gravitational energy exceeding internal, kinetic, or magnetic energy and accelerations due to gravity exceeding those due to magnetic or thermal pressure gradients. In these regions field directions vary strongly, with a slight preference towards being perpendicular to density gradients, as shown by three-dimensional histograms of relative orientation.
△ Less
Submitted 17 November, 2021; v1 submitted 10 August, 2021;
originally announced August 2021.
-
Fast methods for tracking grain coagulation and ionization. I. Analytic derivation
Authors:
Pierre Marchand,
Vincent Guillet,
Ugo Lebreuilly,
Mordecai-Mark Mac Low
Abstract:
Dust grains play a major role in many astrophysical contexts. They affect the chemical, magnetic, dynamical, and optical properties of their environment, from galaxies down to the interstellar medium, star-forming regions, and protoplanetary disks. Their coagulation leads to shifts in their size distribution and ultimately to the formation of planets. However, although the coagulation process is r…
▽ More
Dust grains play a major role in many astrophysical contexts. They affect the chemical, magnetic, dynamical, and optical properties of their environment, from galaxies down to the interstellar medium, star-forming regions, and protoplanetary disks. Their coagulation leads to shifts in their size distribution and ultimately to the formation of planets. However, although the coagulation process is reasonably uncomplicated to numerically implement by itself, it is difficult to couple it with multidimensional hydrodynamics numerical simulations because of its high computational cost. We propose here a simple method for tracking the coagulation of grains at far lower cost. Given an initial grain size distribution, the state of the distribution at time t is solely determined by the value of a single variable integrated along the trajectory, independently of the specific path taken by the grains. Although this method cannot account for processes other than coagulation, it is mathematically exact, fast, inexpensive, and can be used to evaluate the effect of grain coagulation in most astrophysical contexts. It is applicable to all coagulation kernels in which local physical conditions and grain properties can be separated. We also describe another method for calculating the average electric charge of grains and the density of ions and electrons in environments that are shielded from radiation fields, given the density and temperature of the gas, the cosmic-ray ionization rate, and the average mass of the ions. The equations we provide are fast to integrate numerically and can be used in multidimensional numerical simulations to self-consistently calculate on the fly the local resistivities that are required to model nonideal magnetohydrodynamics.
△ Less
Submitted 17 May, 2021; v1 submitted 25 February, 2021;
originally announced March 2021.
-
SIGAME v3: Gas Fragmentation in Post-processing of Cosmological Simulations for More Accurate Infrared Line Emission Modeling
Authors:
Karen Pardos Olsen,
Blakesley Burkhart,
Mordecai-Mark Mac Low,
Robin G. Treß,
Thomas R. Greve,
David Vizgan,
Jay Motka,
Josh Borrow,
Gergö Popping,
Romeel Davé,
Rowan J. Smith,
Desika Narayanan
Abstract:
We present an update to the framework called SImulator of GAlaxy Millimeter/submillimeter Emission (SÍGAME). SÍGAME derives line emission in the far-infrared (FIR) for galaxies in particle-based cosmological hydrodynamics simulations by applying radiative transfer and physics recipes via a post-processing step after completion of the simulation. In this version, a new technique is developed to mod…
▽ More
We present an update to the framework called SImulator of GAlaxy Millimeter/submillimeter Emission (SÍGAME). SÍGAME derives line emission in the far-infrared (FIR) for galaxies in particle-based cosmological hydrodynamics simulations by applying radiative transfer and physics recipes via a post-processing step after completion of the simulation. In this version, a new technique is developed to model higher gas densities by parametrizing the gas density probability distribution function (PDF) in higher resolution simulations for use as a look-up table, allowing for more adaptive PDFs than in previous work. SÍGAME v3 is tested on redshift z = 0 galaxies drawn from the SIMBA cosmological simulation for eight FIR emission lines tracing vastly different interstellar medium phases. Including dust radiative transfer with SKIRT and high resolution photo-ionization models with Cloudy, this new method is able to self-consistently reproduce observed relations between line luminosity and star formation rate in all cases, except for [NII]122, [NII]205 and [OI]63, the luminosities of which are overestimated by median factors of 1.6, 1.2 and 1.2 dex, respectively. We attribute the remaining disagreement with observations to the lack of precise attenuation of the interstellar light on subgrid scales (<200 pc).
△ Less
Submitted 3 December, 2021; v1 submitted 4 February, 2021;
originally announced February 2021.
-
Simulations of the star-forming molecular gas in an interacting M51-like galaxy: cloud population statistics
Authors:
Robin G. Tress,
Mattia C. Sormani,
Rowan J. Smith,
Simon C. O. Glover,
Ralf S. Klessen,
Mordecai-Mark Mac Low,
Paul Clark,
Ana Duarte-Cabral
Abstract:
To investigate how molecular clouds react to different environmental conditions at a galactic scale, we present a catalogue of giant molecular clouds resolved down to masses of $\sim 10$~M$_{\odot}$ from a simulation of the entire disc of an interacting M51-like galaxy and a comparable isolated galaxy. Our model includes time-dependent gas chemistry, sink particles for star formation and supernova…
▽ More
To investigate how molecular clouds react to different environmental conditions at a galactic scale, we present a catalogue of giant molecular clouds resolved down to masses of $\sim 10$~M$_{\odot}$ from a simulation of the entire disc of an interacting M51-like galaxy and a comparable isolated galaxy. Our model includes time-dependent gas chemistry, sink particles for star formation and supernova feedback, meaning we are not reliant on star formation recipes based on threshold densities and can follow the physics of the cold molecular phase. We extract giant molecular clouds at a given timestep of the simulations and analyse their properties. In the disc of our simulated galaxies, spiral arms seem to act merely as snowplows, gathering gas and clouds without dramatically affecting their properties. In the centre of the galaxy, on the other hand, environmental conditions lead to larger, more massive clouds. While the galaxy interaction has little effect on cloud masses and sizes, it does promote the formation of counter-rotating clouds. We find that the identified clouds seem to be largely gravitationally unbound at first glance, but a closer analysis of the hierarchical structure of the molecular interstellar medium shows that there is a large range of virial parameters with a smooth transition from unbound to mostly bound for the densest structures. The common observation that clouds appear to be virialised entities may therefore be due to CO bright emission highlighting a specific level in this hierarchical binding sequence. The small fraction of gravitationally bound structures found suggests that low galactic star formation efficiencies may be set by the process of cloud formation and initial collapse.
△ Less
Submitted 28 June, 2021; v1 submitted 10 December, 2020;
originally announced December 2020.
-
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.
-
The Catalogue for Astrophysical Turbulence Simulations (CATS)
Authors:
B. Burkhart,
S. Appel,
S. Bialy,
J. Cho,
A. J. Christensen,
D. Collins,
C. Federrath,
D. Fielding,
D. Finkbeiner,
A. S. Hill,
J. C. Ibanez-Mejia,
M. R. Krumholz,
A. Lazarian,
M. Li,
P. Mocz,
M. -M. Mac Low,
J. Naiman,
S. K. N. Portillo,
B. Shane,
Z. Slepian,
Y. Yuan
Abstract:
Turbulence is a key process in many fields of astrophysics. Advances in numerical simulations of fluids over the last several decades have revolutionized our understanding of turbulence and related processes such as star formation and cosmic ray propagation. However, data from numerical simulations of astrophysical turbulence are often not made public. We introduce a new simulation-oriented databa…
▽ More
Turbulence is a key process in many fields of astrophysics. Advances in numerical simulations of fluids over the last several decades have revolutionized our understanding of turbulence and related processes such as star formation and cosmic ray propagation. However, data from numerical simulations of astrophysical turbulence are often not made public. We introduce a new simulation-oriented database for the astronomical community: The Catalogue for Astrophysical Turbulence Simulations (CATS), located at www.mhdturbulence.com. CATS includes magnetohydrodynamic (MHD) turbulent box simulation data products generated by the public codes athena++, arepo, enzo, and flash. CATS also includes several synthetic observational data sets, such as turbulent HI data cubes. We also include measured power spectra and 3-point correlation functions from some of these data. We discuss the importance of open source statistical and visualization tools for the analysis of turbulence simulations such as those found in CATS.
△ Less
Submitted 21 October, 2020;
originally announced October 2020.
-
Time-Variable Radio Recombination Line Emission in W49A
Authors:
C. G. De Pree,
D. J. Wilner,
L. E. Kristensen,
R. Galván-Madrid,
W. M. Goss,
R. S. Klessen,
M. -M. Mac Low,
T. Peters,
A. Robinson,
S. Sloman,
M. Rao
Abstract:
We present new Jansky Very Large Array (VLA) images of the central region of the W49A star-forming region at 3.6~cm and at 7~mm at resolutions of 0\farcs15 (1650 au) and 0\farcs04 (440 au), respectively. The 3.6~cm data reveal new morphological detail in the ultracompact \ion{H}{2} region population, as well as several previously unknown and unresolved sources. In particular, source A shows elonga…
▽ More
We present new Jansky Very Large Array (VLA) images of the central region of the W49A star-forming region at 3.6~cm and at 7~mm at resolutions of 0\farcs15 (1650 au) and 0\farcs04 (440 au), respectively. The 3.6~cm data reveal new morphological detail in the ultracompact \ion{H}{2} region population, as well as several previously unknown and unresolved sources. In particular, source A shows elongated, edge-brightened, bipolar lobes, indicative of a collimated outflow, and source E is resolved into three spherical components. We also present VLA observations of radio recombination lines at 3.6~cm and 7~mm, and IRAM Northern Extended Millimeter Array (NOEMA) observations at 1.2~mm. Three of the smallest ultracompact \ion{H}{2} regions (sources A, B2 and G2) all show broad kinematic linewidths, with $Δ$V$_{FWHM}\gtrsim$40~km~s$^{-1}$. A multi-line analysis indicates that broad linewidths remain after correcting for pressure broadening effects, suggesting the presence of supersonic flows. Substantial changes in linewidth over the 21 year time baseline at both 3.6 cm and 7 mm are found for source G2. At 3.6 cm, the linewidth of G2 changed from 31.7$\pm$1.8 km s$^{-1}$ to 55.6$\pm$2.7 km s$^{-1}$, an increase of $+$23.9$\pm$3.4 km s$^{-1}$. The G2 source was previously reported to have shown a 3.6~cm continuum flux density decrease of 40\% between 1994 and 2015. This source sits near the center of a very young bipolar outflow whose variability may have produced these changes.
△ Less
Submitted 9 October, 2020;
originally announced October 2020.
-
Small-Scale Dynamo in Supernova-Driven Interstellar Turbulence
Authors:
Frederick A. Gent,
Mordecai-Mark Mac Low,
Maarit J. Kapyla,
Nishant K. Singh
Abstract:
Magnetic fields grow quickly even at early cosmological times, suggesting the action of a small-scale dynamo (SSD) in the interstellar medium of galaxies. Many studies have focused on idealized turbulent driving of the SSD. Here we simulate more realistic supernova-driven turbulence to determine whether it can drive an SSD. Magnetic field growth occurring in our models appears inconsistent with si…
▽ More
Magnetic fields grow quickly even at early cosmological times, suggesting the action of a small-scale dynamo (SSD) in the interstellar medium of galaxies. Many studies have focused on idealized turbulent driving of the SSD. Here we simulate more realistic supernova-driven turbulence to determine whether it can drive an SSD. Magnetic field growth occurring in our models appears inconsistent with simple tangling of magnetic fields, but consistent with SSD action, reproducing and confirming models by Balsara et al. (2004) that did not include physical resistivity $η$. We vary $η$, as well as the numerical resolution and supernova rate, $\dotσ$, to delineate the regime in which an SSD occurs. For a given $\dotσ$ we find convergence for SSD growth rate with resolution of a parsec. For $\dotσ\simeq\dotσ_{\rm sn}$, with $\dotσ_{\rm sn}$ the solar neighbourhood rate, the critical resistivity below which an SSD occurs is $0.005>η_{\rm crit}>0.001\,\rm kpc^{-1}\,\rm km s^{-1}$, and this increases with the supernova rate. Across the modelled range of 0.5--4 pc resolution we find that for $η<η_{\rm crit}$, the SSD saturates at about 5% of kinetic energy equipartition, independent of growth rate. In the range $0.2\dotσ_{\rm sn}\leq \dotσ\leq8\dotσ_{\rm sn}$ growth rate increases with $\dotσ$. SSDs in the supernova-driven interstellar medium commonly exhibit erratic growth.
△ Less
Submitted 12 March, 2021; v1 submitted 5 October, 2020;
originally announced October 2020.
-
Origin of Weak MgII and Higher Ionization Absorption Lines in an Outflow from an Intermediate-Redshift Dwarf Satellite Galaxy
Authors:
Akimi Fujita,
Toru Misawa,
Jane C. Charlton,
Avery Meiksin,
Mordecai-Mark Mac Low
Abstract:
Observations at intermediate redshifts reveal the presence of numerous, compact, weak MgII absorbers with near to super-solar metallicities, often surrounded by more extended regions that produce CIV and/or OVI absorption in the circumgalactic medium at large impact parameters from luminous galaxies. Their origin and nature remains unclear. We hypothesize that undetected, satellite dwarf galaxies…
▽ More
Observations at intermediate redshifts reveal the presence of numerous, compact, weak MgII absorbers with near to super-solar metallicities, often surrounded by more extended regions that produce CIV and/or OVI absorption in the circumgalactic medium at large impact parameters from luminous galaxies. Their origin and nature remains unclear. We hypothesize that undetected, satellite dwarf galaxies are responsible for producing some of these weak MgII absorbers. We test our hypothesis using gas dynamical simulations of galactic outflows from a dwarf satellite galaxy with a halo mass of $5\times10^{9}$ M$_{\odot}$, which could form in a larger $L^{*}$ halo at z=2, to study the gas interaction in the halo. We find that thin, filamentary, weak MgII absorbers are produced in two stages: 1) when shocked core collapse supernova (SNII) enriched gas descending in a galactic fountain gets shock compressed by upward flows driven by subsequent SNIIs and cools (phase 1), and later, 2) during an outflow driven by Type Ia supernovae that shocks and sweeps up pervasive SNII enriched gas, which then cools (phase 2). The width of the filaments and fragments are $\lesssim~100$ pc, and the smallest ones cannot be resolved at 12.8 pc resolution. The MgII absorbers in our simulations are continuously generated for >150 Myr by shocks and cooling, though each cloud survives for only ~60 Myr. Their metallicity is 10-20% solar metallicity and column density is $<10^{12}$ cm$^{-2}$. They are also surrounded by larger (0.5-1 kpc) CIV absorbers that seem to survive longer. In addition, larger-scale (>1 kpc) CIV and OVI clouds are produced in both expanding and shocked SNII enriched gas which is photoionized by the UV metagalactic radiation at intermediate redshift. Our simulation highlights the possibility of dwarf galactic outflows producing highly enriched multiphase gas.
△ Less
Submitted 21 September, 2020;
originally announced September 2020.
-
The Dynamics, Destruction, and Survival of Supernova-Formed Dust Grains
Authors:
Jonathan D. Slavin,
Eli Dwek,
Mordecai-Mark Mac Low,
Alex S. Hill
Abstract:
Observations have demonstrated that supernovae efficiently produce dust. This is consistent with the hypothesis that supernovae and asymptotic giant branch stars are the primary producers of dust in the Universe. However, there has been a longstanding question of how much of the dust detected in the interiors of young supernova remnants can escape into the interstellar medium. We present new hydro…
▽ More
Observations have demonstrated that supernovae efficiently produce dust. This is consistent with the hypothesis that supernovae and asymptotic giant branch stars are the primary producers of dust in the Universe. However, there has been a longstanding question of how much of the dust detected in the interiors of young supernova remnants can escape into the interstellar medium. We present new hydrodynamical calculations of the evolution of dust grains that were formed in dense ejecta clumps within a Cas A-like remnant. We follow the dynamics of the grains as they decouple from the gas after their clump is hit by the reverse shock. They are subsequently subject to destruction by thermal and kinetic sputtering as they traverse the remnant. Grains that are large enough ($\sim 0.25\,μ$m for silicates and $\sim 0.1\,μ$m for carbonaceous grains) escape into the interstellar medium while smaller grains get trapped and destroyed. However, grains that reach the interstellar medium still have high velocities, and are subject to further destruction as they are slowed down. We find that for initial grain size distributions that include large ($\sim 0.25 - 0.5\,μ$m) grains, 10--20\% of silicate grains can survive, while 30--50\% of carbonaceous grains survive even when the initial size distribution cuts off at smaller ($0.25\,μ$m) sizes. For a 19 M$_{\odot}$ star similar to the progenitor of Cas A, up to 0.1 M$_{\odot}$ of dust can survive if the dust grains formed are large. Thus we show that supernovae under the right conditions can be significant sources of interstellar dust.
△ Less
Submitted 3 September, 2020;
originally announced September 2020.
-
The Role of Stellar Feedback in the Chemical Evolution of a Low Mass Dwarf Galaxy
Authors:
Andrew Emerick,
Greg L. Bryan,
Mordecai-Mark Mac Low
Abstract:
We investigate how each aspect of a multi-channel stellar feedback model drives the chemodynamical evolution of a low-mass, isolated dwarf galaxy using a suite of high-resolution simulations. Our model follows individual star particles sampled randomly from an adopted initial mass function, considering independently feedback from: supernovae; stellar radiation causing photoelectric heating of dust…
▽ More
We investigate how each aspect of a multi-channel stellar feedback model drives the chemodynamical evolution of a low-mass, isolated dwarf galaxy using a suite of high-resolution simulations. Our model follows individual star particles sampled randomly from an adopted initial mass function, considering independently feedback from: supernovae; stellar radiation causing photoelectric heating of dust grains, ionization and associated heating, Lyman-Werner (LW) dissociation of H$_2$, and radiation pressure; and winds from massive main sequence (neglecting their energy input) and asymptotic giant branch (AGB) stars. Radiative transfer is done by ray tracing. We consider the effects each of these processes have on regulating the star formation rate, global properties, multi-phase interstellar medium (ISM), and driving of galactic winds. We follow individual metal species from distinct nucleosynthetic enrichment channels (AGB winds, massive star stellar winds, core collapse and Type Ia supernovae) and pay particular attention to how these feedback processes regulate metal mixing in the ISM, the metal content of outflows, and the stellar abundance patterns in our galaxy. We find that---for a low-metallicity, low-mass dwarf galaxy ---stellar radiation, particularly ionizing radiation and LW radiation, are important sources of stellar feedback whose effects dominate over photoelectric heating and HI radiation pressure. However, feedback is coupled non-linearly, and the inclusion or exclusion of each process produces non-negligible effects. We find strong variations with: the star formation history; the ejection fractions of metals, mass, and energy; and the distribution of elements from different nucleosynthetic sources in both the gas and stars.
△ Less
Submitted 7 July, 2020;
originally announced July 2020.
-
Orbital Migration of Interacting Stellar Mass Black Holes in Disks around Supermassive Black Holes II. Spins and Incoming Objects
Authors:
Amy Secunda,
Jillian Bellovary,
Mordecai-Mark Mac Low,
K. E. Saavik Ford,
Barry McKernan,
Nathan W. C. Leigh,
Wladimir Lyra,
Zsolt Sandor,
Jose I. Adorno
Abstract:
The masses, rates, and spins of merging stellar-mass binary black holes (BBHs) detected by aLIGO and Virgo provide challenges to traditional BBH formation and merger scenarios. An active galactic nucleus (AGN) disk provides a promising additional merger channel, because of the powerful influence of the gas that drives orbital evolution, makes encounters dissipative, and leads to migration. Previou…
▽ More
The masses, rates, and spins of merging stellar-mass binary black holes (BBHs) detected by aLIGO and Virgo provide challenges to traditional BBH formation and merger scenarios. An active galactic nucleus (AGN) disk provides a promising additional merger channel, because of the powerful influence of the gas that drives orbital evolution, makes encounters dissipative, and leads to migration. Previous work showed that stellar mass black holes (sBHs) in an AGN disk migrate to regions of the disk, known as migration traps, where positive and negative gas torques cancel out, leading to frequent BBH formation. Here we build on that work by simulating the evolution of additional sBHs that enter the inner disk by either migration or inclination reduction. We also examine whether the BBHs formed in our models have retrograde or prograde orbits around their centers of mass with respect to the disk, determining the orientation, relative to the disk, of the spin of the merged BBHs. Orbiters entering the inner disk form BBHs with sBHs on resonant orbits near the migration trap. When these sBHs reach ~80 Msun, they form BBHs with sBHs in the migration trap, which over 10 Myr reach ~1000 Msun. We find 68% of the BBHs in our simulation orbit in the retrograde direction, which implies BBHs in our merger channel will have small dimensionless aligned spins, χ_eff. Overall, our models produce BBHs that resemble both the majority of BBH mergers detected thus far (0.66 to 120 Gpc^-3 yr^-1 ) and two recent unusual detections, GW190412 (~0.3 Gpc^-3 yr^-1 ) and GW190521 (~0.1 Gpc^-3 yr^-1 ).
△ Less
Submitted 25 September, 2020; v1 submitted 24 April, 2020;
originally announced April 2020.
-
Modelling of the Effects of Stellar Feedback during Star Cluster Formation Using a Hybrid Gas and N-Body Method
Authors:
Joshua E. Wall,
Mordecai-Mark Mac Low,
Stephen L. W. McMillan,
Ralf S. Klessen,
Simon Portegies Zwart,
Andrew Pellegrino
Abstract:
Understanding the formation of stellar clusters requires following the interplay between gas and newly formed stars accurately. We therefore couple the magnetohydrodynamics code FLASH to the N-body code ph4 and the stellar evolution code SeBa using the Astrophysical Multipurpose Software Environment (AMUSE) to model stellar dynamics, evolution, and collisional N-body dynamics and the formation of…
▽ More
Understanding the formation of stellar clusters requires following the interplay between gas and newly formed stars accurately. We therefore couple the magnetohydrodynamics code FLASH to the N-body code ph4 and the stellar evolution code SeBa using the Astrophysical Multipurpose Software Environment (AMUSE) to model stellar dynamics, evolution, and collisional N-body dynamics and the formation of binary and higher-order multiple systems, while implementing stellar feedback in the form of radiation, stellar winds and supernovae in FLASH. We here describe the algorithms used for each of these processes. We denote this integrated package Torch. We then use this novel numerical method to simulate the formation and early evolution of several examples of open clusters of ~1000 stars formed from clouds with a mass range of 10^3-10^5 M_sun. Analyzing the effects of stellar feedback on the gas and stars of the natal clusters, we find that in these examples, the stellar clusters are resilient to disruption, even in the presence of intense feedback. This can even slightly increase the amount of dense, Jeans unstable gas by sweeping up shells; thus, a stellar wind strong enough to trap its own H II region shows modest triggering of star formation. Our clusters are born moderately mass segregated, an effect enhanced by feedback, and retained after the ejection of their natal gas, in agreement with observations.
△ Less
Submitted 10 October, 2020; v1 submitted 19 March, 2020;
originally announced March 2020.
-
Simulations of the star-forming molecular gas in an interacting M51-like galaxy
Authors:
Robin G. Tress,
Rowan J. Smith,
Mattia C. Sormani,
Simon C. O. Glover,
Ralf S. Klessen,
Mordecai-Mark Mac Low,
Paul C. Clark
Abstract:
We present here the first of a series of papers aimed at better understanding the evolution and properties of giant molecular clouds (GMCs) in a galactic context. We perform high resolution, three-dimensional {\sc arepo} simulations of an interacting galaxy inspired by the well-observed M51 galaxy. Our fiducial simulations include a non-equilibrium, time-dependent, chemical network that follows th…
▽ More
We present here the first of a series of papers aimed at better understanding the evolution and properties of giant molecular clouds (GMCs) in a galactic context. We perform high resolution, three-dimensional {\sc arepo} simulations of an interacting galaxy inspired by the well-observed M51 galaxy. Our fiducial simulations include a non-equilibrium, time-dependent, chemical network that follows the evolution of atomic and molecular hydrogen as well as carbon and oxygen self-consistently. Our calculations also treat gas self-gravity and subsequent star formation (described by sink particles), and coupled supernova feedback. In the densest parts of the simulated interstellar medium (ISM) we reach sub-parsec resolution, granting us the ability to resolve individual GMCs and their formation and destruction self-consistently throughout the galaxy. In this initial work we focus on the general properties of the ISM with a particular focus on the cold star-forming gas. We discuss the role of the interaction with the companion galaxy in generating cold molecular gas and controlling stellar birth. We find that while the interaction drives large-scale gas flows and induces spiral arms in the galaxy, it is of secondary importance in determining gas fractions in the different ISM phases and the overall star-formation rate. The behaviour of the gas on small GMC scales instead is mostly controlled by the self-regulating property of the ISM driven by coupled feedback.
△ Less
Submitted 23 September, 2019;
originally announced September 2019.
-
Simulating Metal Mixing of Both Common and Rare Enrichment Sources in a Low Mass Dwarf Galaxy
Authors:
Andrew Emerick,
Greg L. Bryan,
Mordecai-Mark Mac Low
Abstract:
One-zone models constructed to match observed stellar abundance patterns have been used extensively to constrain the sites of nucleosynthesis with sophisticated libraries of stellar evolution and stellar yields. The metal mixing included in these models is usually highly simplified, although it is likely to be a significant driver of abundance evolution. In this work we use high-resolution hydrody…
▽ More
One-zone models constructed to match observed stellar abundance patterns have been used extensively to constrain the sites of nucleosynthesis with sophisticated libraries of stellar evolution and stellar yields. The metal mixing included in these models is usually highly simplified, although it is likely to be a significant driver of abundance evolution. In this work we use high-resolution hydrodynamics simulations to investigate how metals from individual enrichment events with varying source energies $E_{\rm ej}$ mix throughout the multi-phase interstellar medium (ISM) of a low-mass ($M_{\rm gas}=2\times 10^{6}$~M$_{\odot}$), low-metallicity, isolated dwarf galaxy. These events correspond to the characteristic energies of both common and exotic astrophysical sites of nucleosynthesis, including: asymptotic giant branch winds ($E_{\rm ej}\sim$10$^{46}$~erg), neutron star-neutron star mergers ($E_{\rm ej}\sim$10$^{49}$~erg), supernovae ($E_{\rm ej}\sim$10$^{51}$~erg), and hypernovae ($E_{\rm ej}\sim$10$^{52}$~erg). We find the mixing timescales for individual enrichment sources in our dwarf galaxy to be long (100~Myr--1~Gyr), with a clear trend of increasing homogeneity for the more energetic events. Given these timescales, we conclude that the spatial distribution and frequency of events are important drivers of abundance homogeneity on large scales; rare, low $E_{\rm ej}$ events should be characterized by particularly broad abundance distributions. The source energy $E_{\rm ej}$ also correlates with the fraction of metals ejected in galactic winds, ranging anywhere from 60\% at the lowest energy to 95\% for hypernovae. We conclude by examining how the radial position, local ISM density, and global star formation rate influence these results.
△ Less
Submitted 3 January, 2020; v1 submitted 10 September, 2019;
originally announced September 2019.
-
How do velocity structure functions trace gas dynamics in simulated molecular clouds?
Authors:
Roxana-Adela Chira,
Juan~Camilo Ibáñez-Mejía,
Mordecai-Mark Mac Low,
Thomas Henning
Abstract:
Context. Supersonic disordered flows accompany the formation and evolution of MCs. It has been argued that turbulence can support against gravitational collapse and form hierarchical sub-structures. Aims. We study the time evolution of simulated MCs to investigate: What physical process dominates the driving of turbulent flows? How can these flows be characterised? Do the simulated flows agree wit…
▽ More
Context. Supersonic disordered flows accompany the formation and evolution of MCs. It has been argued that turbulence can support against gravitational collapse and form hierarchical sub-structures. Aims. We study the time evolution of simulated MCs to investigate: What physical process dominates the driving of turbulent flows? How can these flows be characterised? Do the simulated flows agree with observations? Methods. We analyse three MCs that have formed self-consistently within kpc-scale numerical simulations of the ISM. The simulated ISM evolves under the influence of physical processes including self-gravity, stratification, magnetic fields, supernova-driven turbulence, and radiative heating and cooling. We characterise the flows using VSFs with/out density weighting or a density cutoff, and computed in one or three dimensions. Results. In regions with sufficient resolution, the density-weighted VSFs initially appear to follow the expectations for uniform turbulence consistent with Larson's size-velocity relationship. SN blast wave impacts produce short-lived coherent motions at large scales, increasing the scaling exponents for a crossing time. Gravitational contraction drives small-scale motions, producing scaling coefficients that drop or even turn negative as small scales become dominant. Conclusions. We conclude that two different effects coincidentally reproduce Larson's size velocity relationship. Initially, uniform turbulence dominates, so the energy cascade produces VSFs consistent with Larson's relationship. Later, contraction dominates, the density-weighted VSFs become much shallower or even inverted, but the relationship of the global average velocity dispersion of the MCs to their radius follows Larson's relationship, reflecting virial equilibrium or free-fall collapse. The injection of energy by shocks is visible in the VSFs, but decays within a crossing time.
△ Less
Submitted 11 August, 2019;
originally announced August 2019.
-
Dynamical Properties of Molecular-forming Gas Clumps in Galaxies at the Epoch of Reionization
Authors:
T. K. Daisy Leung,
Andrea Pallottini,
Andrea Ferrara,
Mordecai-Mark Mac Low
Abstract:
We study the properties of molecular-forming gas clumps (MGCs) at the epoch of reionization using cosmological zoom-in simulations. We identify MGCs in a z=6 prototypical galaxy ("Althaea") using an H2 density-based clump finder. We compare their mass, size, velocity dispersion, gas surface density, and virial parameter (alpha_vir) to observations. In Althaea, the typical MGC mass and size are Mga…
▽ More
We study the properties of molecular-forming gas clumps (MGCs) at the epoch of reionization using cosmological zoom-in simulations. We identify MGCs in a z=6 prototypical galaxy ("Althaea") using an H2 density-based clump finder. We compare their mass, size, velocity dispersion, gas surface density, and virial parameter (alpha_vir) to observations. In Althaea, the typical MGC mass and size are Mgas=10^6.5 Msun and R=45-100 pc, which are comparable to those found in nearby spirals and starburst galaxies. MGCs are highly supersonic and supported by turbulence, with rms velocity dispersions of sigma_gas=20-100 km s^-1 and pressure of P/k_B=10^7.6 K cm^-3 (i.e., >1000x with respect to the Milky Way), similar to those found in nearby and z~2 gas-rich starburst galaxies. In addition, we perform stability analysis to understand the origin and dynamical properties of MGCs. We find that MGCs are globally stable in the main disk of Althaea. Densest regions where star formation is expected to take place in clouds and cores on even smaller scales instead have lower alpha_vir and Toomre-Q values. Detailed studies of the star-forming gas dynamics at the epoch of reionization thus require a spatial resolution of <40 pc (=0.01"), which is within reach with the Atacama Large (sub-)Millimeter Array and the Next Generation Very Large Array.
△ Less
Submitted 22 April, 2020; v1 submitted 5 July, 2019;
originally announced July 2019.
-
The Implications of Local Fluctuations in the Galactic Midplane for Dynamical Analysis in the Gaia Era
Authors:
Angus Beane,
Robyn E. Sanderson,
Melissa K. Ness,
Kathryn V. Johnston,
Douglas Grion Filho,
Mordecai-Mark Mac Low,
Daniel Anglés-Alcázar,
David W. Hogg,
Chervin F. P. Laporte
Abstract:
Orbital properties of stars, computed from their six-dimensional phase space measurements and an assumed Galactic potential, are used to understand the structure and evolution of the Galaxy. Stellar actions, computed from orbits, have the attractive quality of being invariant under certain assumptions and are therefore used as quantitative labels of a star's orbit. We report a subtle but important…
▽ More
Orbital properties of stars, computed from their six-dimensional phase space measurements and an assumed Galactic potential, are used to understand the structure and evolution of the Galaxy. Stellar actions, computed from orbits, have the attractive quality of being invariant under certain assumptions and are therefore used as quantitative labels of a star's orbit. We report a subtle but important systematic error that is induced in the actions as a consequence of local midplane variations expected for the Milky Way. This error is difficult to model because it is non-Gaussian and bimodal, with neither mode peaking on the null value. An offset in the vertical position of the Galactic midplane of $\sim15\,\text{pc}$ for a thin disk-like orbit or $\sim 120\,\text{pc}$ for a thick disk-like orbit induces a $25\%$ systematic error in the vertical action $J_z$. In FIRE simulations of Milky Way-mass galaxies, these variations are on the order of $\sim100\,\text{pc}$ at the solar circle. From observations of the mean vertical velocity variation of $\sim5\text{--}10\,\text{km}\,\text{s}^{-1}$ with radius, we estimate that the Milky Way midplane variations are $\sim60\text{--}170\,\text{pc}$, consistent with three-dimensional dust maps. Action calculations and orbit integrations, which assume the global and local midplanes are identical, are likely to include this induced error, depending on the volume considered. Variation in the local standard of rest or distance to the Galactic center causes similar issues. The variation of the midplane must be taken into account when performing dynamical analysis across the large regions of the disk accessible to Gaia and future missions.
△ Less
Submitted 22 August, 2019; v1 submitted 21 May, 2019;
originally announced May 2019.
-
AGN (and other) astrophysics with Gravitational Wave Events
Authors:
K. E. Saavik Ford,
Imre Bartos,
Barry McKernan,
Zoltan Haiman,
Alessandra Corsi,
Azadeh Keivani,
Szabolcs Marka,
Rosalba Perna,
Matthew Graham,
Nicholas P. Ross,
Daniel Stern,
Jillian Bellovary,
Emanuele Berti,
Matthew O'Dowd,
Wladimir Lyra,
Mordecai-Mark Mac Low,
Zsuzsanna Marka
Abstract:
The stellar mass binary black hole (sBBH) mergers presently detected by LIGO may originate wholly or in part from binary black hole mergers embedded in disks of gas around supermassive black holes. Determining the contribution of these active galactic nucleus (AGN) disks to the sBBH merger rate enables us to uniquely measure important parameters of AGN disks, including their typical density, aspec…
▽ More
The stellar mass binary black hole (sBBH) mergers presently detected by LIGO may originate wholly or in part from binary black hole mergers embedded in disks of gas around supermassive black holes. Determining the contribution of these active galactic nucleus (AGN) disks to the sBBH merger rate enables us to uniquely measure important parameters of AGN disks, including their typical density, aspect ratio, and lifetime, thereby putting unique limits on an important element of galaxy formation. For the first time, gravitational waves will allow us to reveal the properties of the hidden interior of AGN disks, while electromagnetic radiation (EM) probes the disk photosphere. The localization of sBBH merger events from LIGO is generally insufficient for association with a single EM counterpart. However, the contribution to the LIGO event rate from rare source types (such as AGNs) can be determined on a statistical basis. To determine the contribution to the sBBH rate from AGNs in the next decade requires: {\em 1) a complete galaxy catalog for the LIGO search volume, 2) strategic multi-wavelength EM follow-up of LIGO events and 3) significant advances in theoretical understanding of AGN disks and the behavior of objects embedded within them.}
△ Less
Submitted 22 March, 2019;
originally announced March 2019.
-
Planet formation: The case for large efforts on the computational side
Authors:
Wladimir Lyra,
Thomas Haworth,
Bertram Bitsch,
Simon Casassus,
Nicolás Cuello,
Thayne Currie,
Andras Gáspár,
Hannah Jang-Condell,
Hubert Klahr,
Nathan Leigh,
Giuseppe Lodato,
Mordecai-Mark Mac Low,
Sarah Maddison,
George Mamatsashvili,
Colin McNally,
Andrea Isella,
Sebastián Pérez,
Luca Ricci,
Debanjan Sengupta,
Dimitris Stamatellos,
Judit Szulágyi,
Richard Teague,
Neal Turner,
Orkan Umurhan,
Jacob White
, et al. (32 additional authors not shown)
Abstract:
Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been a…
▽ More
Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.
△ Less
Submitted 11 March, 2019;
originally announced March 2019.
-
Collisional N-Body Dynamics Coupled to Self-Gravitating Magnetohydrodynamics Reveals Dynamical Binary Formation
Authors:
Joshua E. Wall,
Stephen L. W. McMillan,
Mordecai-Mark Mac Low,
Ralf S. Klessen,
Simon Portegies Zwart
Abstract:
We describe a star cluster formation model that includes individual star formation from self-gravitating, magnetized gas, coupled to collisional stellar dynamics. The model uses the Astrophysical Multi-purpose Software Environment (AMUSE) to integrate an adaptive-mesh magnetohydrodynamics code (FLASH) with a fourth order Hermite N-body code (ph4), a stellar evolution code (SeBa), and a method for…
▽ More
We describe a star cluster formation model that includes individual star formation from self-gravitating, magnetized gas, coupled to collisional stellar dynamics. The model uses the Astrophysical Multi-purpose Software Environment (AMUSE) to integrate an adaptive-mesh magnetohydrodynamics code (FLASH) with a fourth order Hermite N-body code (ph4), a stellar evolution code (SeBa), and a method for resolving binary evolution (multiples). This combination yields unique star formation simulations that allow us to study binaries formed dynamically from interactions with both other stars and dense, magnetized gas subject to stellar feedback during the birth and early evolution of stellar clusters. We find that for massive stars, our simulations are consistent with the observed dynamical binary fractions and mass ratios. However, our binary fraction drops well below observed values for lower mass stars, presumably due to unincluded binary formation during initial star formation. Further, we observe a build up of binaries near the hard-soft boundary that may be an important mechanism driving early cluster contraction.
△ Less
Submitted 4 January, 2019;
originally announced January 2019.
-
Diffusion and Concentration of Solids in the Dead Zone of a Protoplanetary Disk
Authors:
Chao-Chin Yang,
Mordecai-Mark Mac Low,
Anders Johansen
Abstract:
The streaming instability is a promising mechanism to drive the formation of planetesimals in protoplanetary disks. To trigger this process, it has been argued that sedimentation of solids onto the mid-plane needs to be efficient and therefore that a quiescent gaseous environment is required. It is often suggested that dead-zone or disk-wind structure created by non-ideal magnetohydrodynamical (MH…
▽ More
The streaming instability is a promising mechanism to drive the formation of planetesimals in protoplanetary disks. To trigger this process, it has been argued that sedimentation of solids onto the mid-plane needs to be efficient and therefore that a quiescent gaseous environment is required. It is often suggested that dead-zone or disk-wind structure created by non-ideal magnetohydrodynamical (MHD) effects meets this requirement. However, simulations have shown that the midplane of a dead zone is not completely quiescent. In order to examine the concentration of solids in such an environment, we use the local-shearing-box approximation to simulate a particle-gas system with an Ohmic dead zone including mutual drag force between the gas and the solids. We systematically compare the evolution of the system with ideal or non-ideal MHD, with or without back-reaction drag force from particles on gas, and with varying solid abundances. Similar to previous investigations of dead zone dynamics, we find that particles of dimensionless stopping time $τ_s=0.1$ do not sediment appreciably more than those in ideal magneto-rotational turbulence, resulting in a vertical scale height an order of magnitude larger than in a laminar disk. Contrary to the expectation that this should curb the formation of planetesimals, we nevertheless find that strong clumping of solids still occurs in the dead zone when solid abundances are similar to the critical value for a laminar environment. This can be explained by the weak radial diffusion of particles near the mid-plane. The results imply that the sedimentation of particles to the mid-plane is not a necessary criterion for the formation of planetesimals by the streaming instability.
△ Less
Submitted 11 October, 2018;
originally announced October 2018.
-
Metal Mixing and Ejection in Dwarf Galaxies is Dependent on Nucleosynthetic Source
Authors:
Andrew Emerick,
Greg L. Bryan,
Mordecai-Mark Mac Low,
Benoit Côté,
Kathryn V. Johnston,
Brian W. O'Shea
Abstract:
Using a high resolution simulation of an isolated dwarf galaxy, accounting for multi-channel stellar feedback and chemical evolution on a star-by-star basis, we investigate how each of 15 metal species are distributed within our multi-phase interstellar medium (ISM) and ejected from our galaxy by galactic winds. For the first time, we demonstrate that the mass fraction probability distribution fun…
▽ More
Using a high resolution simulation of an isolated dwarf galaxy, accounting for multi-channel stellar feedback and chemical evolution on a star-by-star basis, we investigate how each of 15 metal species are distributed within our multi-phase interstellar medium (ISM) and ejected from our galaxy by galactic winds. For the first time, we demonstrate that the mass fraction probability distribution functions (PDFs) of individual metal species in the ISM are well described by a piecewise log-normal and power-law distribution. The PDF properties vary within each ISM phase. Hot gas is dominated by recent enrichment, with a significant power-law tail to high metal fractions, while cold gas is predominately log-normal. In addition, elements dominated by asymptotic giant branch (AGB) wind enrichment (e.g. N and Ba) mix less efficiently than elements dominated by supernova enrichment (e.g. $α$ elements and Fe). This result is driven by the differences in source energetics and source locations, particularly the higher chance compared to massive stars for AGB stars to eject material into cold gas. Nearly all of the produced metals are ejected from the galaxy (only 4% are retained), but over 20% of metals dominated by AGB enrichment are retained. In dwarf galaxies, therefore, elements synthesized predominately through AGB winds should be both overabundant and have a larger spread compared to elements synthesized in either core collapse or Type Ia supernovae. We discuss the observational implications of these results, their potential use in developing improved models of galactic chemical evolution, and their generalization to more massive galaxies.
△ Less
Submitted 7 November, 2018; v1 submitted 4 September, 2018;
originally announced September 2018.
-
Stellar Radiation is Critical for Regulating Star Formation and Driving Outflows in Low Mass Dwarf Galaxies
Authors:
Andrew Emerick,
Greg L. Bryan,
Mordecai-Mark Mac Low
Abstract:
Effective stellar feedback is used in models of galaxy formation to drive realistic galaxy evolution. Models typically include energy injection from supernovae as the dominant form of stellar feedback, often in some form of sub-grid recipe. However, it has been recently suggested that pre-SN feedback (stellar winds or radiation) is necessary in high-resolution simulations of galaxy evolution to pr…
▽ More
Effective stellar feedback is used in models of galaxy formation to drive realistic galaxy evolution. Models typically include energy injection from supernovae as the dominant form of stellar feedback, often in some form of sub-grid recipe. However, it has been recently suggested that pre-SN feedback (stellar winds or radiation) is necessary in high-resolution simulations of galaxy evolution to properly regulate star formation and properties of the interstellar medium (ISM). Following these processes is computationally challenging, so many prescriptions model this feedback approximately, accounting for the local destruction of dense gas clouds around newly formed stars in lieu of a full radiative transfer calculation. In this work we examine high resolution simulations (1.8~pc) of an isolated dwarf galaxy with detailed stellar feedback tracked on a star-by-star basis. By following stellar ionizing radiation with an adaptive ray-tracing radiative transfer method, we test its importance in regulating star formation and driving outflows in this galaxy. We find that including ionizing radiation reduces the star formation rate (SFR) by over a factor of 5, and is necessary to produce the ISM conditions needed for supernovae to drive significant outflows. We find that a localized approximation for radiation feedback is sufficient to regulate the SFR on short timescales, but does not allow significant outflows. Short and long range radiation effects are both important in driving the evolution of our low-metallicity, low-mass dwarf galaxy. Generalizing these results to more massive galaxies would be a valuable avenue of future research.
△ Less
Submitted 19 September, 2018; v1 submitted 1 August, 2018;
originally announced August 2018.
-
Flux Density Variations at 3.6 cm in the Massive Star-Forming Region W49A
Authors:
C. G. De Pree,
R. Galván-Madrid,
W. M. Goss,
R. S. Klessen,
M. -M. Mac Low,
T. Peters,
D. Wilner,
J. Bates,
T. Melo,
B. Presler-Marshall,
R. Webb-Forgus
Abstract:
A number of ultracompact H II regions in Galactic star forming environments have been observed to vary significantly in radio flux density on timescales of 10-20 years. Theory predicted that such variations should occur when the accretion flow that feeds a young massive star becomes unstable and clumpy. We have targeted the massive star-forming region W49A with the Karl G. Jansky Very Large Array…
▽ More
A number of ultracompact H II regions in Galactic star forming environments have been observed to vary significantly in radio flux density on timescales of 10-20 years. Theory predicted that such variations should occur when the accretion flow that feeds a young massive star becomes unstable and clumpy. We have targeted the massive star-forming region W49A with the Karl G. Jansky Very Large Array (VLA) for observations at 3.6 cm with the B-configuration at 0.8'' resolution, to compare to nearly identical observations taken almost 21 years earlier (February 2015 and August 1994). Most of the sources in the crowded field of ultracompact and hypercompact H II regions exhibit no significant changes over this time period. However, one source, W49A/G2, decreased by 20% in peak intensity (from 71+/-4 mJy/beam to 57+/-3 mJy/beam), and 40% in integrated flux (from 0.109+/-0.011 Jy to 0.067+/-0.007 Jy), where we cite 5 sigma errors in peak intensity, and 10% errors in integrated flux. We present the radio images of the W49A region at the two epochs, the difference image that indicates the location of the flux density decrease, and discuss explanations for the flux density decrease near the position of W49A/G2.
△ Less
Submitted 27 July, 2018;
originally announced July 2018.
-
Simulating an Isolated Dwarf Galaxy with Multi-Channel Feedback and Chemical Yields from Individual Stars
Authors:
Andrew Emerick,
Greg L. Bryan,
Mordecai-Mark Mac Low
Abstract:
In order to better understand the relationship between feedback and galactic chemical evolution, we have developed a new model for stellar feedback at grid resolutions of only a few parsecs in global disk simulations, using the adaptive mesh refinement hydrodynamics code Enzo. For the first time in galaxy scale simulations, we simulate detailed stellar feedback from individual stars including asym…
▽ More
In order to better understand the relationship between feedback and galactic chemical evolution, we have developed a new model for stellar feedback at grid resolutions of only a few parsecs in global disk simulations, using the adaptive mesh refinement hydrodynamics code Enzo. For the first time in galaxy scale simulations, we simulate detailed stellar feedback from individual stars including asymptotic giant branch winds, photoelectric heating, Lyman-Werner radiation, ionizing radiation tracked through an adaptive ray-tracing radiative transfer method, and core collapse and Type Ia supernovae. We furthermore follow the star-by-star chemical yields using tracer fields for 15 metal species: C, N, O, Na, Mg, Si, S, Ca, Mn, Fe, Ni, As, Sr, Y, and Ba. We include the yields ejected in massive stellar winds, but greatly reduce the winds' velocities due to computational constraints. We describe these methods in detail in this work and present the first results from 500~Myr of evolution of an isolated dwarf galaxy with properties similar to a Local Group, low-mass dwarf galaxy. We demonstrate that our physics and feedback model is capable of producing a dwarf galaxy whose evolution is consistent with observations in both the Kennicutt-Schmidt relationship and extended Schmidt relationship. Effective feedback drives outflows with a greater metallicity than the ISM, leading to low metal retention fractions consistent with observations. Finally, we demonstrate that these simulations yield valuable information on the variation in mixing behavior of individual metal species within the multi-phase interstellar medium.
△ Less
Submitted 1 October, 2018; v1 submitted 18 July, 2018;
originally announced July 2018.
-
Orbital Migration of Interacting Stellar Mass Black Holes in Disks around Supermassive Black Holes
Authors:
Amy Secunda,
Jillian Bellovary,
Mordecai-Mark Mac Low,
K. E. Saavik Ford,
Barry McKernan,
Nathan Leigh,
Wladimir Lyra,
Zsolt Sandor
Abstract:
The merger rate of stellar-mass black hole binaries (sBHBs) inferred by the Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO) suggests the need for an efficient source of sBHB formation. Active galactic nucleus (AGN) disks are a promising location for the formation of these sBHBs, as well as binaries of other compact objects, because of powerful torques exerted by the gas disk. T…
▽ More
The merger rate of stellar-mass black hole binaries (sBHBs) inferred by the Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO) suggests the need for an efficient source of sBHB formation. Active galactic nucleus (AGN) disks are a promising location for the formation of these sBHBs, as well as binaries of other compact objects, because of powerful torques exerted by the gas disk. These gas torques cause orbiting compact objects to migrate towards regions in the disk where inward and outward torques cancel, known as migration traps. We simulate the migration of stellar mass black holes in an example of a model AGN disk, using an augmented N-body code that includes analytic approximations to migration torques, stochastic gravitational forces exerted by turbulent density fluctuations in the disk, and inclination and eccentricity dampening produced by passages through the gas disk, in addition to the standard gravitational forces between objects. We find that sBHBs form rapidly in our model disk as stellar-mass black holes migrate towards the migration trap. These sBHBs are likely to subsequently merge on short time-scales. The process continues, leading to the build-up of a population of over-massive stellar-mass black holes. The formation of sBHBs in AGN disks could contribute significantly to the sBHB merger rate inferred by LIGO.
△ Less
Submitted 12 November, 2019; v1 submitted 8 July, 2018;
originally announced July 2018.
-
Effect of the heating rate on the stability of the three-phase interstellar medium
Authors:
Alex S. Hill,
Mordecai-Mark Mac Low,
Andrea Gatto,
Juan C. Ibáñez-Mejía
Abstract:
We investigate the impact of the far ultraviolet (FUV) heating rate on the stability of the three-phase interstellar medium using three-dimensional simulations of a $1$ kpc$^2$, vertically-extended domain. The FUV heating rate sets the range of thermal pressures across which the cold ($\sim10^2$ K) and warm ($\sim10^4$ K) neutral media (CNM and WNM) can coexist in equilibrium. Even absent a variab…
▽ More
We investigate the impact of the far ultraviolet (FUV) heating rate on the stability of the three-phase interstellar medium using three-dimensional simulations of a $1$ kpc$^2$, vertically-extended domain. The FUV heating rate sets the range of thermal pressures across which the cold ($\sim10^2$ K) and warm ($\sim10^4$ K) neutral media (CNM and WNM) can coexist in equilibrium. Even absent a variable star formation rate regulating the FUV heating rate, the gas physics keeps the pressure in the two-phase regime: because radiative heating and cooling processes happen on shorter timescales than sound wave propagation, turbulent compressions tend to keep the interstellar medium within the CNM-WNM pressure regime over a wide range of heating rates. The thermal pressure is set primarily by the heating rate with little influence from the hydrostatics. The vertical velocity dispersion adjusts as needed to provide hydrostatic support given the thermal pressure: when the turbulent pressure $\langleρ\rangleσ_z^2$ is calculated over scales $\gtrsim500$ pc, the thermal plus turbulent pressure approximately equals the weight of the gas. The warm gas volume filling fraction is $0.2<f_w<0.8$ over a factor of less than three in heating rate, with $f_w$ near unity at higher heating rates and near zero at lower heating rates. We suggest that cosmological simulations that do not resolve the CNM should maintain an interstellar thermal pressure within the two-phase regime.
△ Less
Submitted 5 July, 2018; v1 submitted 14 June, 2018;
originally announced June 2018.
-
Modelling supernova driven turbulence
Authors:
Frederick A. Gent,
Mordecai-Mark Mac Low,
Maarit J. Käpylä,
Graeme R. Sarson,
James F. Hollins
Abstract:
High Mach number shocks are ubiquitous in interstellar turbulence. The Pencil Code is particularly well suited to the study of magnetohydrodynamics in weakly compressible turbulence and the numerical investigation of dynamos because of its high-order advection and time evolution algorithms. However, the high-order algorithms and lack of Riemann solver to follow shocks make it less well suited to h…
▽ More
High Mach number shocks are ubiquitous in interstellar turbulence. The Pencil Code is particularly well suited to the study of magnetohydrodynamics in weakly compressible turbulence and the numerical investigation of dynamos because of its high-order advection and time evolution algorithms. However, the high-order algorithms and lack of Riemann solver to follow shocks make it less well suited to handling high Mach number shocks, such as those produced by supernovae (SNe). Here, we outline methods required to enable the code to efficiently and accurately model SNe, using parameters that allow stable simulation of SN-driven turbulence, in order to construct a physically realistic galactic dynamo model. These include the resolution of shocks with artificial viscosity, thermal conductivity, and mass diffusion; the correction of the mass diffusion terms; and a novel generalization of the Courant condition to include all source terms in the momentum and energy equations. We test our methods with the numerical solution of the one-dimensional (1D) Riemann shock tube (Sod, J. Comput. Phys. 1978, 27), also extended to a 1D adiabatic shock with parameters and Mach number relevant to SN shock evolution, including shocks with radiative losses. We extend our test with the three-dimensional (3D) numerical simulation of individual SN remnant evolution for a range of ambient gas densities typical of the interstellar medium and compare these to the analytical solutions of Sedov-Taylor (adiabatic) and the snowplough and Cioffi, McKee and Bertschinger (Astrophys. J. 1988, 334) results incorporating cooling and heating processes. We show that our new timestep algorithm leads to linear rather than quadratic resolution dependence as the strength of the artificial viscosity varies, because of the corresponding change in the strength of interzone gradients.
△ Less
Submitted 18 June, 2019; v1 submitted 5 June, 2018;
originally announced June 2018.
-
Cosmic Ray Driven Outflows in an Ultraluminous Galaxy
Authors:
Akimi Fujita,
Mordecai-Mark Mac Low
Abstract:
In models of galaxy formation, feedback driven both by supernova (SN) and active galactic nucleus (AGN) is not efficient enough to quench star formation in massive galaxies. Models of smaller galaxies have suggested that cosmic rays (CRs) play a major role in expelling material from the star forming regions by diffusing SN energy to the lower density outskirts. We therefore run gas dynamical simul…
▽ More
In models of galaxy formation, feedback driven both by supernova (SN) and active galactic nucleus (AGN) is not efficient enough to quench star formation in massive galaxies. Models of smaller galaxies have suggested that cosmic rays (CRs) play a major role in expelling material from the star forming regions by diffusing SN energy to the lower density outskirts. We therefore run gas dynamical simulations of galactic outflows from a galaxy contained in a halo with $5\times10^{12}$~M$_{\odot}$ that resembles a local ultraluminous galaxy, including both SN thermal energy and a treatment of CRs using the same diffusion approximation as Salem & Bryan. We find that CR pressure drives a low-density bubble beyond the edge of the shell swept up by thermal pressure, but the main bubble driven by SN thermal pressure overtakes it later, which creates a large-scale biconical outflow. CRs diffusing into the disk are unable to entrain its gas in the outflows, yielding a mass-loading rate of only $\sim0.1~\%$ with varied CR diffusion coefficients. We find no significant difference in mass-loading rates in SN driven outflows with or without CR pressure. Our simulations strongly suggest that it is hard to drive a heavily mass-loaded outflow with CRs from a massive halo potential, although more distributed star formation could lead to a different result.
△ Less
Submitted 26 April, 2018;
originally announced April 2018.
-
Dust concentration and chondrule formation
Authors:
Alexander Hubbard,
Mordecai-Mark Mac Low,
Denton S. Ebel
Abstract:
Meteoritical and astrophysical models of planet formation make contradictory predictions for dust concentration factors in chondrule forming regions of the solar nebula. Meteoritical and cosmochemical models strongly suggest that chondrules, a key component of the meteoritical record, formed in regions with solids-to-gas mass ratios orders of magnitude above background. However, models of dust gra…
▽ More
Meteoritical and astrophysical models of planet formation make contradictory predictions for dust concentration factors in chondrule forming regions of the solar nebula. Meteoritical and cosmochemical models strongly suggest that chondrules, a key component of the meteoritical record, formed in regions with solids-to-gas mass ratios orders of magnitude above background. However, models of dust grain dynamics in protoplanetary disks struggle to surpass factors of a few outside of very brief windows in the lifetime of the dust grains. Worse, those models do not predict significant concentration factors for dust grains the size of chondrule precursors. We briefly develop the difficulty in concentrating dust particles in the context of nebular chondrule formation and show that the disagreement is sufficiently stark that cosmochemists should explore ideas that might revise the concentration factor requirements downwards.
△ Less
Submitted 27 March, 2018;
originally announced March 2018.
-
On the rate of black hole binary mergers in galactic nuclei due to dynamical hardening
Authors:
Nathan W. C. Leigh,
Aaron M. Geller,
B. McKernan,
K. E. S. Ford,
M. -M. Mac Low,
J. Bellovary,
Z. Haiman,
W. Lyra,
J. Samsing,
M. O'Dowd,
B. Kocsis,
S. Endlich
Abstract:
We assess the contribution of dynamical hardening by direct three-body scattering interactions to the rate of stellar-mass black hole binary (BHB) mergers in galactic nuclei. We derive an analytic model for the single-binary encounter rate in a nucleus with spherical and disk components hosting a super-massive black hole (SMBH). We determine the total number of encounters $N_{\rm GW}$ needed to ha…
▽ More
We assess the contribution of dynamical hardening by direct three-body scattering interactions to the rate of stellar-mass black hole binary (BHB) mergers in galactic nuclei. We derive an analytic model for the single-binary encounter rate in a nucleus with spherical and disk components hosting a super-massive black hole (SMBH). We determine the total number of encounters $N_{\rm GW}$ needed to harden a BHB to the point that inspiral due to gravitational wave emission occurs before the next three-body scattering event. This is done independently for both the spherical and disk components. Using a Monte Carlo approach, we refine our calculations for $N_{\rm GW}$ to include gravitational wave emission between scattering events. For astrophysically plausible models we find that typically $N_{\rm GW} \lesssim$ 10.
We find two separate regimes for the efficient dynamical hardening of BHBs: (1) spherical star clusters with high central densities, low velocity dispersions and no significant Keplerian component; and (2) migration traps in disks around SMBHs lacking any significant spherical stellar component in the vicinity of the migration trap, which is expected due to effective orbital inclination reduction of any spherical population by the disk. We also find a weak correlation between the ratio of the second-order velocity moment to velocity dispersion in galactic nuclei and the rate of BHB mergers, where this ratio is a proxy for the ratio between the rotation- and dispersion-supported components. Because disks enforce planar interactions that are efficient in hardening BHBs, particularly in migration traps, they have high merger rates that can contribute significantly to the rate of BHB mergers detected by the advanced Laser Interferometer Gravitational-Wave Observatory.
△ Less
Submitted 28 November, 2017;
originally announced November 2017.