While the stellar Initial Mass Function (IMF) appears to be close to universal within the Milky Way galaxy, it is strongly suspected to be different in the primordial Universe, where molecular hydrogen cooling is less efficient and the gas temperature can be higher by a factor of 30. In between these extreme cases, the gas temperature varies depending on the environment, metallicity and radiation background. In this paper we explore if changes of the gas temperature affect the IMF of the stars considering fragmentation and accretion. The fragmentation behavior depends mostly on the Jeans mass at the turning point in the equation of state where a transition occurs from an approximately isothermal to an adiabatic regime due to dust opacities. The Jeans mass at this transition in the equation of state is always very similar, independent of the initial temperature, and therefore the initial mass of the fragments is very similar. Accretion on the other hand is strongly temperature dependent. We argue that the latter becomes the dominant process for star formation efficiencies above 5 - 7 %, increasing the average mass of the stars.

To explain the observed population of supermassive black holes at z~7, very massive seed black holes or, alternatively, super-Eddington scenarios are needed to reach final masses of the order of 10^9 solar masses. A popular explanation for massive seeds has been the direct collapse model, which predicts the formation of a single massive object due to the direct collapse of a massive gas cloud. Simulations over the last years have however shown that such a scenario is very difficult to achieve. A realistic model of black hole formation should therefore take fragmentation into account, and consider the interaction between stellar-dynamical and gas-dynamical processes. We present here numerical simulations pursued with the AMUSE code, employing an approximate treatment of the gas. Based on these simulations, we show that very massive black holes of 10^4-10^5 solar masses may form depending on the gas supply and the accretion onto the protostars.

The presence of supermassive black holes (SMBHs) of a few billion solar masses at very high redshift has motivated us to study how these massive objects formed during the first billion years after the Big Bang. The most promising model that has been proposed to explain this is the direct collapse of protogalactic gas clouds. In this scenario, very high accretion rates are needed to form massive objects early on and the suppression of $\mathrm H_2$ cooling is important in regulating the fragmentation. Recent studies have shown that if we use a strong radiation background, the hydrogen molecules are destroyed, favoring the high accretion rates and therefore producing objects of very high mass. In this work we study the impact of UV radiation fields in a primordial gas cloud using the recently coupled code GRADSPH-KROME for the modeling of gravitational collapse including primordial chemistry to explore the fragmentation in AU scales and hence the formation of first SMBHs. We found that to suppress the formation of $\mathrm H_2$ a very high value of $J_{21}$ is required, because of that we conclude that the direct collapse black holes (DCBHs) are very unlikely to be an explanation for the formation of the first SMBHs.

We present here a three-dimesional hydrodynamical simulation for star formation. Our aim is to explore the effect of the metal-line cooling on the thermodynamics of the star-formation process. We explore the effect of changing the metallicty of the gas from $Z/Z_{\odot}=10^{-4}$ to $Z/Z_{\odot}=10^{-2}$. Furthermore, we explore the implications of using the observational abundance pattern of a CEMP-no star, which have been considered to be the missing second-generation stars, the so-called Pop. III.2 stars. In order to pursue our aim, we modelled the microphysics by employing the public astrochemistry package KROME, using a chemical network which includes sixteen chemical species (H, H$^{+}$, H$^{-}$, He, He$^{+}$, He$^{++}$, e$^{-}$, H$_{2}$, H$_{2}^{+}$, C, C$^{+}$, O, O$^{+}$, Si, Si$^{+}$, and Si$^{++}$). We couple KROME with the fully three-dimensional Smoothed-particle hydrodynamics (SPH) code GRADSPH. With this framework we investigate the collapse of a metal-enhanced cloud, exploring the fragmentation process and the formation of stars. We found that the metallicity has a clear impact on the thermodynamics of the collapse, allowing the cloud to reach the CMB temperature floor for a metallicity $Z/Z_{\odot}=10^{-2}$, which is in agreement with previous work. Moreover, we found that adopting the abundance pattern given by the star SMSS J031300.36-670839.3 the thermodynamics behavior is very similar to simulations with a metallicity of $Z/Z_{\odot}=10^{-2}$, due to the high carbon abundance. As long as only metal line cooling is considered, our results support the metallicity threshold proposed by previous works, which will very likely regulate the first episode of fragmentation and potentially determine the masses of the resulting star clusters.

The first stars are known to form in primordial gas, either in minihalos with about $10^6$~M$_\odot$ or so-called atomic cooling halos of about $10^8$~M$_\odot$. Simulations have shown that gravitational collapse and disk formation in primordial gas yield dense stellar clusters. In this paper, we focus particularly on the formation of protostellar binary systems, and aim to quantify their properties during the early stage of their evolution. For this purpose, we combine the smoothed particle hydrodynamics code GRADSPH with the astrochemistry package KROME. The GRADSPH-KROME framework is employed to investigate the collapse of primordial clouds in the high-density regime, exploring the fragmentation process and the formation of binary systems. We observe a strong dependence of fragmentation on the strength of the turbulent Mach number $\mathcal{M}$ and the rotational support parameter $\beta{}$. Rotating clouds show significant fragmentation, and have produced several Pop.~III proto-binary systems. We report maximum and minimum mass accretion rates of $2.31 \times 10^{-1}$~M$_{\odot}$ yr$^{-1}$ and $2.18\times 10^{-4}$~M$_{\odot}$ yr$^{-1}$. The mass spectrum of the individual Pop III proto-binary components ranges from $0.88$~M$_{\odot}$ to $31.96$~M$_{\odot}$ and has a sensitive dependence on the Mach number $\mathcal{M}$ as well as on the rotational parameter $\beta{}$. We also report a range from $\sim0.01$ to $\sim1$ for the mass ratio of our proto-binary systems.

The origin of very low-mass stars (VLMS) and brown dwarfs (BDs) is still an unresolved topic of star formation. We here present numerical simulations of the formation of VLMS, BDs, and planet mass objects (planemos) resulting from the gravitational collapse and fragmentation of solar mass molecular cores with varying rotation rates and initial density perturbations. Our simulations yield various types of binary systems including the combinations VLMS-VLMS, BD-BD, planemo-planemo, VLMS-BD, VLMS-planemos, BD-planemo. Our scheme successfully addresses the formation of wide VLMS and BD binaries with semi-major axis up to 441 AU and produces a spectrum of mass ratios closer to the observed mass ratio distribution (q > 0.5). Molecular cores with moderate values of the ratio of kinetic to gravitational potential energy (0.16 <= beta <= 0.21) produce planemos. Solar mass cores with rotational parameters beta outside of this range yield either VLMS/BDs or a combination of both. With regard to the mass ratios we find that for both types of binary systems the mass ratio distribution varies in the range 0.31 <= q <= 0.74. We note that in the presence of radiative feedback, the length scale of fragmentation would increase by approximately two orders of magnitude, implying that the formation of binaries may be efficient for wide orbits, while being suppressed for short-orbit systems.

Observations show a large spread in the luminosities of young protostars, which are frequently explained in the context of episodic accretion. We here test this scenario using numerical simulations following the collapse of a solar mass molecular cloud using the GRADSPH code, varying the strength of the initial perturbations and the temperature of the cores. A specific emphasis of this paper is to investigate the role of binaries and multiple systems in the context of episodic accretion, and to compare their evolution to the evolution in isolated fragments. Our models form a variety of low mass protostellar objects including single, binary and triple systems with binaries more active in exhibiting episodic accretion than isolated protostars. We also find a general decreasing trend for the average mass accretion rate over time, suggesting that the majority of the protostellar mass is accreted within the first 10^5 years. This result can potentially help to explain the surprisingly low average luminosities in the majority of the protostellar population.

We report the results of a numerical study on the initial formation stages of low-mass protostellar binary systems. We determine the separation of protostellar binaries formed as a function of the initial thermal state by varying the initial temperature in a slightly modified version of the Burkert and Bodenheimer collapse test. We find that the outcome is highly sensitive to both the initial temperature of the cloud and the initial amplitude of azimuthal density perturbation A. For A=10 %, variations of only 1 unit Kelvin below 10 K lead to changes of up to 100 AU ( i.e. of order 30 %) in the instantaneous separation, whereas for this small A the initial temperatures above 10 K yield, instead of a binary, a single low-mass fragment that never reaches protostellar densities. Protostellar binaries, however, do emerge when the perturbation amplitude is increased from 10 % to 25 %. We also investigate the impact of the critical density which governs the transition from isothermal to adiabatic thermodynamic behaviour of the collapsing gas. We find that the critical density not only affects the overall structural evolution of the gas envelope, but also the size of the rotating disk structures formed during collapse as well as the number of protostellar fragments resulting from the final fragmentation of the disks. This mechanism can give rise to young protostellar objects constituting bound multiple stellar systems.

An attempt is made here to revisit structure formation in a proto-stellar cloud during the early phase of evolution. Molecular cloud subjected to a set of various initial conditions in terms of initial temperature and amplitude of azimuthal density perturbation is investigated numerically. Special emphasis remained on the analysis of ring and spiral type instabilities that have shown dependence on certain initial conditions chosen for a rotating solar mass cloud of molecular hydrogen. Generally, a star forming hydrogen gas is considered to be initially at 10K. We have found that a possible oscillation around this typical value can affect the fate of a collapsing cloud in terms of its evolving structural properties leading to proto-star formation. We explored the initial temperature range of cloud between 8K to 12K and compared physical properties of each within the first phase of proto-star formation. We suggest that the spiral structures are more likely to form in strongly perturbed molecular cores that initiate their phase of collapse from temperatures below 10K. Whereas, cores with initial temperatures above 10K develop, instead of spiral, a ring type structure which subsequently experiences the fragmentation. A transition from spiral to ring instability can be observed at a typical initial temperature of 10K.