No Arabic abstract
Magnetic fields play an important role in the dynamics of present-day molecular clouds. Recent work has shown that magnetic fields are equally important for primordial clouds, which form the first stars in the Universe. While the primordial magnetic field strength on cosmic scales is largely unconstrained, theoretical models strongly suggest that a weak seed field existed in the early Universe. We study how the amplification of such a weak field can influence the evolution of accretion discs around first stars, and thus affect the primordial initial mass function (IMF). We perform a suite of 3D ideal magneto-hydrodynamic (MHD) simulations with different initial field strengths and numerical resolutions. We find that, in simulations with sufficient spatial resolution to resolve the Jeans scale during the collapse, even initially weak magnetic fields grow exponentially to become dynamically important due to both the so-called small-scale turbulent dynamo and the large-scale mean-field dynamo. Capturing the small-scale dynamo action depends primarily on how well we resolve the Jeans length, while capturing the large-scale dynamo depends on the Jeans resolution as well as the maximum absolute resolution. Provided enough resolution, we find that fragmentation does not depend strongly on the initial field strength, because even weak fields grow to become strong. However, fragmentation in runs with magnetic fields differs significantly from those without magnetic fields. We conclude that the development of dynamically strong magnetic fields during the formation of the first stars is likely inevitable, and that these fields had a significant impact on the primordial IMF.
Disks of gas accreting onto supermassive black holes are thought to power active galactic nuclei (AGN). Stars may form in gravitationally unstable regions of these disks, or may be captured from nuclear star clusters. Because of the dense gas environment, the evolution of such embedded stars can diverge dramatically from those in the interstellar medium. This work extends previous studies of stellar evolution in AGN disks by exploring a variety of ways that accretion onto stars in AGN disks may differ from Bondi accretion. We find that tidal effects from the supermassive black hole significantly alter the evolution of stars in AGN disks, and that our results do not depend critically on assumptions about radiative feedback on the accretion stream. Thus, in addition to depending on $rho/c_s^3$, the fate of stars in AGN disks depends sensitively on the distance to and mass of the supermassive black hole. This affects where in the disk stellar explosions occur, where compact remnants form and potentially merge to produce gravitational waves, and where different types of chemical enrichment take place.
We investigate how a protoplanetary discs susceptibility to gravitational instabilities and fragmentation depends on the mass of its host star. We use 1D disc models in conjunction with 3D SPH simulations to determine the critical disc-to-star mass ratios at which discs become unstable against fragmentation, finding that discs become increasingly prone to the effects of self-gravity as we increase the host star mass. The actual limit for stability is sensitive to the disc temperature, so if the disc is optically thin stellar irradiation can dramatically stabilise discs against gravitational instability. However, even when this is the case we find that discs around $2$M$_{odot}$ stars are prone to fragmentation, which will act to produce wide-orbit giant planets and brown dwarfs. The consequences of this work are two-fold: that low mass stars could in principle support high disc-to-star mass ratios, and that higher mass stars have discs that are more prone to fragmentation, which is qualitatively consistent with observations that favour high-mass wide-orbit planets around higher mass stars. We also find that the initial masses of these planets depends on the temperature in the disc at large radii, which itself depends on the level of stellar irradiation.
Magnetic fields play an important role for the formation of stars in both local and high-redshift galaxies. Recent studies of dynamo amplification in the first dark matter haloes suggest that significant magnetic fields were likely present during the formation of the first stars in the Universe at redshifts of 15 and above. In this work, we study how these magnetic fields potentially impact the initial mass function (IMF) of the first stars. We perform 200 high-resolution, three-dimensional (3D), magneto-hydrodynamic (MHD) simulations of the collapse of primordial clouds with different initial turbulent magnetic field strengths as predicted from turbulent dynamo theory in the early Universe, forming more than 1100 first stars in total. We detect a strong statistical signature of suppressed fragmentation in the presence of strong magnetic fields, leading to a dramatic reduction in the number of first stars with masses low enough that they might be expected to survive to the present day. Additionally, strong fields shift the transition point where stars go from being mostly single to mostly multiple to higher masses. However, irrespective of the field strength, individual simulations are highly chaotic, show different levels of fragmentation and clustering, and the outcome depends on the exact realisation of the turbulence in the primordial clouds. While these are still idealised simulations that do not start from cosmological initial conditions, our work shows that magnetic fields play a key role for the primordial IMF, potentially even more so than for the present-day IMF.
We use a suite of SPH simulations to investigate the susceptibility of protoplanetary discs to the effects of self-gravity as a function of star-disc properties. We also include passive irradiation from the host star using different models for the stellar luminosities. The critical disc-to-star mass ratio for axisymmetry (for which we produce criteria) increases significantly for low-mass stars. This could have important consequences for increasing the potential mass reservoir in a proto Trappist-1 system, since even the efficient Ormel et al. (2017) formation model will be influenced by processes like external photoevaporation, which can rapidly and dramatically deplete the dust reservoir. The aforementioned scaling of the critical $M_d/M_*$ for axisymmetry occurs in part because the Toomre $Q$ parameter has a linear dependence on surface density (which promotes instability) and only an $M_*^{1/2}$ dependence on shear (which reduces instability), but also occurs because, for a given $M_d/M_*$, the thermal evolution depends on the host star mass. The early phase stellar irradiation of the disc (for which the luminosity is much higher than at the zero age main sequence, particularly at low stellar masses) can also play a key role in significantly reducing the role of self-gravity, meaning that even Solar mass stars could support axisymmetric discs a factor two higher in mass than usually considered possible. We apply our criteria to the DSHARP discs with spirals, finding that self-gravity can explain the observed spirals so long as the discs are optically thick to the host star irradiation.
We examine the cosmic evolution of a stellar initial mass function (IMF) in galaxies that varies with the Jeans mass in the interstellar medium, paying particular attention to the K-band stellar mass to light ratio (M/L_K) of present-epoch massive galaxies. We calculate the typical Jeans mass using high-resolution hydrodynamic simulations coupled with a fully radiative model for the ISM, which yields a parameterisation of the IMF characteristic mass as a function of galaxy star formation rate (SFR). We then calculate the star formation histories of galaxies utilising an equilibrium galaxy growth model coupled with constraints on the star formation histories set by abundance matching models. We find that at early times, energetic coupling between dust and gas drive warm conditions in the ISM, yielding bottom-light/top- heavy IMFs associated with large ISM Jeans masses for massive star-forming galaxies. Owing to the remnants of massive stars that formed during the top-heavy phases at early times, the resultant M/L_K(sigma) in massive galaxies at the present epoch is increased relative to the non- varying IMF case. At late times, lower cosmic ray fluxes allow for cooler ISM temperatures in massive galaxies, and hence newly formed clusters will exhibit bottom-heavy IMFs, further increasing M/L_K(sigma). Our central result is hence that a given massive galaxy may go through both top-heavy and bottom-heavy IMF phases during its lifetime, though the bulk of the stars form during a top-heavy phase. Qualitatively, the variations in M/L_K(sigma) with galaxy mass are in agreement with observations, however, our model may not be able to account for bottom-heavy mass functions as indicated by stellar absorption features.