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 the first stars, and thus affect the primordial initial mass function (IMF). We perform a suite of 3D 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 \textit{small-scale turbulent dynamo} and the \textit{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.

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 perform detailed spectroscopy of the X-ray brightest supernova remnant (SNR) in the Large Magellanic Cloud (LMC), N132D, using Chandra archival observations. By analyzing the spectra of the entire well-defined rim, we determine the mean abundances for O, Ne, Mg, Si, S and Fe for the local LMC environment. We find evidence of enhanced O on the north-western and S on the north-eastern blast wave. By analyzing spectra interior to the remnant, we confirm the presence of a Si-rich relatively hot plasma (> 1.5 kev) that is also responsible for the Fe K emission. Chandra images show that the Fe K emission is distributed throughout the interior of the southern half of the remnant but does not extend out to the blast wave. We estimate the progenitor mass to be $15\pm5\,M_{\odot}$ using abundance ratios in different regions that collectively cover a large fraction of the remnant, as well as from the radius of the forward shock compared with models of an explosion in a cavity created by stellar winds. We fit ionizing and recombining plasma models to the Fe K emission and find that the current data cannot distinguish between the two, hence the origin of the high-temperature plasma remains uncertain. Our analysis is consistent with N132D being the result of a core-collapse supernova in a cavity created by its intermediate mass progenitor.

The adiabatic index of H$_2\,$ ($\gamma_{\mathrm{H_2}}$) is non-constant at temperatures between $100-10^4\,\mathrm{K}$ due to the large energy spacing between its rotational and vibrational modes. For the formation of the first stars at redshifts 20 and above, this variation can be significant because primordial molecular clouds are in this temperature range due to the absence of efficient cooling by dust and metals. We study the possible importance of variations in $\gamma_{\mathrm{H_2}}$ for the primordial initial mass function by carrying out 80 3D gravito-hydrodynamic simulations of collapsing clouds with different random turbulent velocity fields, half using fixed $\gamma_{\rm H_2} = 7/5$ in the limit of classical diatomic gas (used in earlier works) and half using an accurate quantum mechanical treatment of $\gamma_{\mathrm{H_2}}$. We use the adaptive mesh refinement code FLASH with the primordial chemistry network from KROME for this study. The simulation suite produces almost 400 stars, with masses from $0.02 - 50$ M$_\odot$ (mean mass $\sim 10.5\,\mathrm{M_{\odot}}$ and mean multiplicity fraction $\sim 0.4$). While the results of individual simulations do differ when we change our treatment of $\gamma_{\mathrm{H_2}}$, we find no statistically significant differences in the overall mass or multiplicity distributions of the stars formed in the two sets of runs. We conclude that, at least prior to the onset of radiation feedback, approximating H$_2$ as a classical diatomic gas with $\gamma_{\rm H_2} = 7/5$ does not induce significant errors in simulations of the fragmentation of primordial gas. Nonetheless, we recommend using the accurate formulation of the H$_2\,$ adiabatic index in primordial star formation studies since it is not computationally more expensive and provides a better treatment of the thermodynamics.

We probe the star formation properties of the gas in AzTEC-1 in the COSMOS field, one of the best resolved and brightest starburst galaxies at $z \approx 4.3$, forming stars at a rate > 1000 $\mathrm{M_{\odot}}\,\mathrm{yr^{-1}}$. Using recent ALMA observations, we study star formation in the galaxy nucleus and an off-center star-forming clump and measure a median star formation rate (SFR) surface density of $\Sigma^{\mathrm{nucleus}}_{\mathrm{SFR}} = 270\pm54$ and $\Sigma^{\mathrm{sfclump}}_{\mathrm{SFR}} = 170\pm38\,\mathrm{M_{\odot}}\,\mathrm{yr}^{-1}\,\mathrm{kpc}^{-2}$, respectively. Following the analysis by Sharda et al. (2018), we estimate the molecular gas mass, freefall time and turbulent Mach number in these regions to predict $\Sigma_{\mathrm{SFR}}$ from three star formation relations in the literature. The Kennicutt-Schmidt (Kennicutt 1998, KS) relation, which is based on the gas surface density, underestimates the $\Sigma_{\mathrm{SFR}}$ in these regions by a factor 2-3. The $\Sigma_{\mathrm{SFR}}$ we calculate from the single-freefall model of Krumholz et al. 2012 (KDM) is consistent with the measured $\Sigma_{\mathrm{SFR}}$ in the nucleus and the star-forming clump within the uncertainties. The turbulence-regulated star formation relation by Salim et al. 2015 (SFK) agrees slightly better with the observations than the KDM relation. Our analysis reveals that an interplay between turbulence and gravity can help sustain high SFRs in high-redshift starbursts. It can also be extended to other high- and low-redshift galaxies thanks to the high angular resolution and sensitivity of ALMA observations.

Using high-resolution (sub-kiloparsec scale) submillimeter data obtained by ALMA, we analyze the star formation rate (SFR), gas content and kinematics in SDP 81, a gravitationally-lensed star-forming galaxy at redshift 3. We estimate the SFR surface density ($\Sigma_{\mathrm{SFR}}$) in the brightest clump of this galaxy to be $357^{+135}_{-85}\,\mathrm{M_{\odot}\,yr^{-1}\,kpc^{-2}}$, over an area of $0.07\pm0.02\,\mathrm{kpc}^2$. Using the intensity-weighted velocity of CO$\,$(5-4), we measure the turbulent velocity dispersion in the plane-of-the-sky and find $\sigma_{\mathrm{v,turb}} = 37\pm5\,\mathrm{km\,s}^{-1}$ for the star-forming clump, in good agreement with previous estimates along the line of sight. Our measurements of gas surface density, freefall time and turbulent Mach number reveal that the role of turbulence is vital to explaining the observed SFR in this clump. While the Kennicutt Schmidt (KS) relation predicts a SFR surface density of $\Sigma_{\mathrm{SFR,KS}} = 52\pm17\,\mathrm{M_{\odot}\,yr^{-1}\,kpc^{-2}}$, the single-freefall model by Krumholz, Dekel and McKee (KDM) predicts $\Sigma_{\mathrm{SFR,KDM}} = 106\pm37\,\mathrm{M_{\odot}\,yr^{-1}\,kpc^{-2}}$. In contrast, the multi-freefall (turbulence) model by Salim, Federrath and Kewley (SFK) gives $\Sigma_{\mathrm{SFR,SFK}} = 491^{+139}_{-194}\,\mathrm{M_{\odot}\,yr^{-1}\,kpc^{-2}}$. Although the SFK relation overestimates the SFR in this clump (possibly due to the ignorance of magnetic field), it provides the best prediction among the available models. Finally, we compare the star formation and gas properties of this high-redshift galaxy to local star-forming regions and find that the SFK relation provides the best estimates of SFR in both local and high-redshift galaxies.

We compute the concentrations of five transition elements (Cr, Fe, Co, Ni and Zn) in supernova presolar grains (Silicon Carbide Type X) from the time they condense till the end of free expansion phase, via condensation and implantation. We consider relative velocities of these elements with respect to grains as they condense and evolve at temperatures $\le$ 2000 K, use zonal nucleosynthesis yields for three core collapse supernovae models - 15 M\textsubscript{\(\odot\)}, 20 M\textsubscript{\(\odot\)} and 25 M\textsubscript{\(\odot\)} and an ion target simulator SDTrimSP to model their implantation onto the grains. Simulations from SDTrimSP show that maximal implantation in the core of the grain is possible, contrary to previous studies. We find that the 15 M\textsubscript{\(\odot\)} model best explains the measured concentrations of SiC X grains obtained from Murchison meteorite. For grains which contained $\ge$ 300 ppm of Fe and Ni, we find implantation fraction to be $\le$ 0.25 for most probable differential zonal velocities in this phase which implies that condensation is dominant than implantation unless the grain did not spend much time in the supernova environment. We show that radioactive corrections and mixing from the innermost zones is vital to explaining the excess Ni (condensed as well as implanted) obtained in laboratory measurements. This mixing also explains the relative abundances of Co and Ni with respect to Fe simultaneously. The model developed can be used to predict concentrations of all other elements in various presolar grains condensed in supernova ejecta and matched against measured concentrations in grains found in meteorites.