Search

Early-type stars have convective cores due to a steep temperature gradient produced by the CNO cycle. These cores can host dynamos, and the generated magnetic fields can be relevant to explain the magnetism observed in Ap/Bp stars. Our main objective is to characterise the convective core dynamos and differential rotation, and to do the first quantitative analysis of the relation between magnetic activity cycle and rotation period. We use numerical 3D star-in-a-box simulations of a $2.2~M_\odot$ A-type star with a convective core of roughly $20\%$ of the stellar radius surrounded by a radiative envelope. Rotation rates from 8 to 20 days are explored. We use two models of the entire star, and an additional zoom set, where $50\%$ of the radius is retained. The simulations produce hemispheric core dynamos with cycles and typical magnetic field strengths around 60 kG. However, only a very small fraction of the magnetic energy is able to reach the surface. The cores have solar-like differential rotation, and a substantial part of the radiative envelope has quasi-rigid rotation. In the most rapidly rotating cases the magnetic energy in the core is roughly 40\% of the kinetic energy. Finally, we find that the magnetic cycle period $P_\mathrm{cyc}$ increases with decreasing the rotation period $P_\mathrm{rot}$ which is also observed in many simulations of solar-type stars. Our simulations indicate that a strong hemispherical core dynamo arises routinely, but that it is not enough the explain the surface magnetism of Ap/Bp stars. Nevertheless, as the core dynamo produces dynamically relevant magnetic fields it should not be neglected when other mechanisms are explored., Comment: Article accepted for publication in A&A

We investigate the properties of circumstellar discs (CDs) produced in hydrodynamical simulations of gravoturbulent core collapse considering Kolmogorov and Burger-type turbulence. We report that massive discs are more prevalent in the Kolmogorov regime than for Burger-type turbulence. A significant number of discs are formed with a radius of $\sim$ 15 au in both cases. However, the number of extended discs with radii $>$ 15 au is significantly larger in case of Kolmogorov turbulence. The two regimes of turbulence, in general, yield disc radii in the ranges of 7 au $-$ 30 au, and 13 au $-$ 39 au, respectively. The corresponding ranges of the disc masses are 30.37 $M_{\rm Jup}$ $-$ 0.92 M$_{\odot}$, and 2.09 $M_{\rm Jup}$ $-$ 0.13 M$_{\odot}$, respectively. Moreover, the ratio $M_{\rm disc}$/$M_{\rm star}$ is higher in models of Kolmogorov-type turbulence than in models of Burgers-type turbulence. We do not find any correlation between $R_{\rm disc}$ and $M_{\rm disc}$ over the explored range of initial temperatures (8 K $-$ 14 K) and the type of turbulence. Also, for these initial thermal variations, the turbulent circumstellar disc structures do not exhibit signs of turbulent diffusion. Nonetheless, both sub and supersonic velocity dispersions cause variations in the specific angular momentum (AM) of infalling gas, especially for CDs with radii $\sim$ 16 au $-$ 21 au. The radial profiles of CDs do not correlate with the initial conditions.

This paper explores the theoretical relation between star clusters and black holes within, focusing on the potential role of nuclear star clusters (NSCs), globular clusters (GCs), and ultra compact dwarf galaxies (UCDs) as environments that allow for black hole formation via stellar collisions. This study aims to identify the optimal conditions for stellar collisions across a range of stellar systems leading to the formation of very massive stars that subsequently collapse into black holes. We analyze data from numerical simulations and observations of diverse stellar systems, encompassing various initial conditions, initial mass functions, and evolution scenarios. We computed a critical mass, determined by the interplay of collision time, system age, and initial properties of the star cluster. The efficiency of black hole formation ($\epsilon_{\mathrm{BH}}$) is defined as the ratio of initial stellar mass divided by critical mass. We find that stellar systems with a ratio of initial stellar mass over critical mass above 1 exhibit a high efficiencies in terms of black hole formation, ranging from $30-100\%$. While there is some scatter, potentially attributed to complex system histories and the presence of gas, the results highlight the potential for achieving high efficiencies via a purely collisional channel in black hole formation. In conclusion, this theoretical exploration elucidates the connection between star clusters and black hole formation. The study underscores the significance of UCDs, GCs, and NSCs as environments conducive to the black hole formation scenario via stellar collisions. The defined black hole formation efficiency ($\epsilon_{\mathrm{BH}}$) is shown to be influenced by the ratio of the initial stellar mass to the critical mass., Comment: 16 pages, 5 figures

In this paper we study the qualitative features induces by corrections to GR coming from String Theory, on the shadows of rotating black holes. We deal with the slowly rotating black hole solutions up to order $\mathcal{O}(a^3)$, to first order in $\alpha'$, including also the dilaton. We provide a detailed characterization of the geometry, as well as the innermost stable circular orbit (ISCO) and photon ring, and then we proceed to obtain the black hole images within the relativistic thin-disk model. We characterize the images by computing the diameter, displacement and asymmetry. A comparison with the Kerr case, indicates that all these quantities grow due to the $\alpha'$ correction, and that the departure from GR for different observable is enhanced depending on the angle of view, namely for the diameter the maximum departure is obtained when the system is face-on, while for the displacement and asymmetry the departure from GR is maximized for edge-on point of view., Comment: 14 pages, 8 figures. V2: small changes, references added, accepted in PRD

The potential importance of magnetic fields during structure formation and gravitational collapse in the early Universe has been shown in several studies. In particular, magnetic field amplification by the small-scale dynamo plays an important role in addition to the pure amplification expected from gravitational collapse. In this paper, we study the small-scale dynamo for halos of $\gtrsim10^7$ M$_\odot$ collapsing at $z\gtrsim12$, under different ambient conditions due to the strength of the Lyman-Werner background. Additionally, we estimate the approximate saturation level by varying the initial magnetic field strength. We performed cosmological magnetohydrodynamical simulations for three distinct halos of $\sim10^7$ M$_{\odot}$ at $z\geq13$ by varying the Jeans resolution from $32-256$ cells and employed Lyman Werner background flux of strengths $10^2-10^5$ in units of $J_{21}$, where $J_{21}=10^{-21}$ erg$/$cm$^2/$sr$/$s$/$Hz. To follow the chemical and thermal evolution of the gas we made use of the KROME package. In addition to the compression by collapse, we find magnetic field amplification via the dynamo both in the regimes of atomic and molecular hydrogen cooling. Moreover, we find a lower saturation level in the molecular hydrogen cooling regime. This behaviour can be understood due to the generally reduced radial infall velocities and vorticities in this regime, as well as the higher Mach numbers of the gas, which give rise to a smaller saturation ratio. Our results overall suggest that the dynamo operates over a large range of conditions in the collapsing gas., Comment: 20 pages, 16 figures, 3 tables. Accepted for publication in A&A

Black holes (BH) are objects described by General Relativity (GR), particularly by the Kerr solution. However, this solution is based on simplifying assumptions. To achieve a more realistic approach, we investigate the impact of deviations from the Kerr solution on the imaging of supermassive black holes (SMBH). We conduct General Relativistic Radiation Transport simulations using a Kerr-Like metric, which includes four nonlinear free deviation functions and is stationary, axis-symmetric, and asymptotically flat. The RAPTOR I code, an open-source ray-tracing code capable of handling arbitrary spacetimes, is employed to generate the images. We analyze the asymmetry, diameter, and shadow displacement of the BH to compare them with those of a Kerr BH. Our findings confirm that these deviations significantly affect the shape of the shadow. However, we acknowledge the importance of considering the matter distribution and emissivities surrounding the BH, as the results cannot be solely inferred from the shape of the photon ring., Comment: 3 pages, 2 figures,Accepted for publication in Bolet\'in de la Asociaci\'on Argentina de Astronom\'ia

The presence of supermassive black holes at redshift z > 6 raises some questions about their formation and growth in the early universe. Due to the construction of new telescopes like the ELT to observe and detect SMBHs, it will be useful to derive theoretical estimates for the population and to compare observations and model predictions in the future. In consequence our main goal is to estimate the population of SMBHs using a semi-analytic code known as Galacticus which is a code for the formation and evolution of galaxies where we are about to include different scenarios for SMBHs formation indicating the initial mass of the black hole seed, its formation conditions and recipes for the evolution of the components of the galaxies. We found that the principal mechanism of growing SMBHs is is via galaxy mergers and accretion of matter. For the comparison of our results with observations, we calculate the radius of influence of the black hole to estimate which part of the population could be detected, leading to relations similar to the observed ones., Comment: 3 pages, 3 figures

The aim of this study is to explore the magnetic and flow properties of fully convective M dwarfs as a function of rotation period Prot and magnetic Reynolds ReM and Prandlt numbers PrM. We performed three-dimensional simulations of fully convective stars using a star-in-a-box setup. This setup allows global dynamo simulations in a sphere embedded in a Cartesian cube. The equations of non-ideal magnetohydrodynamics were solved with the Pencil Code. We used the stellar parameters of an M5 dwarf with 0.21M_odot at three rotation rates corresponding to rotation periods (Prot): 43, 61 and 90 days, and varied the magnetic Prandtl number in the range from 0.1 to 10. We found systematic differences in the behaviour of the large-scale magnetic field as functions of rotation and PrM. For the simulations with Prot = 43 days and PrM <= 2, we found cyclic large-scale magnetic fields. For PrM > 2 the cycles vanish and field shows irregular reversals. In simulations with Prot = 61 days for PrM <= 2 the cycles are less clear and the reversal are less periodic. In the higher-PrM cases, the axisymmetric mean field shows irregular variations. For the slowest rotation case with Prot = 90 days, the field has an important dipolar component for PrM <= 5. For the highest PrM the large-scale magnetic field is predominantly irregular at mid-latitudes, with quasi-stationary fields near the poles. For the simulations with cycles, the cycle period length slightly increases with increasing ReM., Comment: 11 pages, 12 figures, submitted to A&A

According to our understanding of stellar evolution, early-type stars have radiative envelopes and convective cores due to a steep temperature gradient produced by the CNO cycle. Some of these stars (mainly, the subclasses Ap and Bp) have strong magnetic fields, enough to be directly observed using the Zeeman effect. Here, we present 3D magnetohydrodynamic simulations of an $2 ~M_{\odot}$ A-type star using the star-in-a-box model. Our goal is to explore if the modeled star is able to maintain a magnetic field as strong as the observed ones, via a dynamo driven by its convective core, or via maintaining a stable fossil field configuration coming from its early evolutionary stages, using different rotation rates. We created two models, a partially radiative and a fully radiative one, which are determined by the value of the heat conductivity. Our model is able to explore both scenarios, including convection-driven dynamos., Comment: Conference proceeding (Bolet\'in de la Asociaci\'on Argentina de Astronom\'ia), 3 pages, 3 figures

Context: A detailed study of the close interacting binary V4142\,Sgr based on photometric and spectroscopic analysis is presented.This system belongs to the enigmatic class of Algol-like variables showing a long photometric cycle of unknown nature. Aims: Performing photometric data-mining and spectroscopic observations covering the orbital cycle, we obtain the orbital parameters and the stellar properties of the binary system, along with the physical properties of the accretion disk located around the hot star. Insights on the evolutive path of the system are obtained. Methods: The light curve was modeled through an inverse modeling method using a theoretical light curve of the binary system, considering the light curve contribution of both stars and the accretion disk of the hot star to obtain the fundamental parameters. To constrain the main stellar parameters the mass ratio was fixed, as well as the donor temperature using the obtained values from our spectroscopic analysis including deblending methods to isolate the spectral lines of the stellar components. The system parameters were compared with a grid of binary star evolutive models in order to get insights on the evolutionary history of the system. Results: The orbital period and the long cycle were re-calculated and found to be of $30.633 \pm 0.002 ~\mathrm{days}$ and $1201 \pm 14 ~\mathrm{days}$. The spectral analysis reveals H$\alpha$ double emission with a persistent $V \leq R$ asymmetry which is considered evidence of a possible wind emergin from the hotspot region...

The infrared-radio correlation (IRRC) of star-forming galaxies can be used to estimate their star formation rate (SFR) based on the radio continuum luminosity at MHz-GHz frequencies. For its application in future deep radio surveys, it is crucial to know whether the IRRC persists at high redshift z. Previous works have reported that the 1.4 GHz IRRC correlation of star-forming galaxies is nearly z-invariant up to z=4, but depends strongly on stellar mass M. This should be taken into account for SFR calibrations based on radio luminosity. To understand the physical cause behind the M dependence of the IRRC and its properties at higher z, we constructed a phenomenological model for galactic radio emission. Our model is based on a dynamo-generated magnetic field and a steady-state cosmic ray population. It includes a number of free parameters as well as observed scaling relations. We find that the resulting spread of the infrared-to-radio luminosity ratio, q(z, M), with respect to M is mostly determined by the scaling of the galactic radius with M, while the absolute value of the q(z, M) curves decreases with more efficient conversion of supernova energy to magnetic fields and cosmic rays. Decreasing the slope of the cosmic ray injection spectrum, aCR, results in higher radio luminosity, decreasing the absolute values of the q(z, M) curves. Our model reproduces the observed dependence of the IRRC on M and z when the efficiency of supernova-driven turbulence is 5%, 10% of the kinetic energy is converted into magnetic energy, and aCR = 3. For galaxies with intermediate to high M (10^9.5-10^11 M_sun), our model results in an IRRC that is nearly independent of z. For galaxies with lower M (M=10^8.5 M_sun), we find that the IR-to-radio flux ratio increases with increasing redshift. This matches the observational data in that mass bin which currently, however, only extends to z~1.5., Comment: 20 pages, 18 figures

The centres of galaxies host nuclear stellar clusters, supermassive black holes, or both. The origin of this dichotomy is still a mystery. Nuclear stellar clusters are the densest stellar system in the Universe, so they are ideal places for runaway collisions to occur. Previous studies have proposed the possible existence of a critical mass scale in such clusters, for which the occurrence of collisions becomes very frequent and leads to the formation of a very massive object. While it is difficult to directly probe this scenario with simulations, we here aim for a proof of concept using toy models where the occurrence of such a transition is shown based on simplified compact systems, where the typical evolution time-scales will be faster compared to the real Universe. Indeed our simulations confirm that such a transition takes place and that up to 50 per cent of the cluster mass can go into the formation of a central massive object for clusters that are above the critical mass scale. Our results thus support the proposed new scenario on the basis of idealized simulations. A preliminary analysis of observed nuclear star clusters shows similar trends related to the critical mass as in our simulations. We further discuss the caveats for the application of the proposed scenario in real nuclear star clusters., Comment: 15 pages, 9 figures

While large numbers of supermassive black holes have been detected at z>6, their origin is still essentially unclear. Numerical simulations have shown that the conditions for the classical direct collapse scenario are very restrictive and fragmentation is very difficult to be avoided. We thus consider here a more general case of a dense massive protostar cluster at low metallicity (<~ 10^{-3} Z_solar) embedded in gas. We estimate the mass of the central massive object, formed via collisions and gas accretion, considering the extreme cases of a logarithmically flat and a Salpeter-type initial mass function. Objects with masses of at least 10^4 solar could be formed for inefficient radiative feedback, whereas ~10^3 solar mass objects could be formed when the accretion time is limited via feedback. These masses will vary depending on the environment and could be considerably larger, particularly due to the continuous infall of gas into the cloud. As a result, one may form intermediate mass black holes of ~ 10^4 solar masses or more. Upcoming observations with the James Webb Space Telescope (JWST) and other observatories may help to detect such massive black holes and their environment, thereby shedding additional light on such a formation channel., Comment: 9 pages, 1 figure, 1 table. Accepted for publication with MNRAS

We present a detailed photometric and spectroscopic analysis of DD CMa, based on published survey photometry and new spectroscopic data. We find an improved orbital period of $P_\mathrm{o}= 2.0084530 \pm 0.0000006 ~\mathrm{d}$. Our spectra reveal H$\beta$ and H$\alpha$ absorptions with weak emission shoulders and we also find color excess in the WISE multiband photometry, interpreted as signatures of circumstellar matter. We model the $V$-band orbital light curve derived from the ASAS and ASAS-SN surveys, assuming a semidetached configuration and using the mass ratio and temperature of the hotter star derived from our spectroscopic analysis. Our model indicates that the system consists of a B 2.5 dwarf and a B 9 giant of radii 3.2 and 3.7 $\mathrm{R_{\odot}}$, respectively, orbiting in a circular orbit of radius 6.75 $\mathrm{R_{\odot}}$. We also found $M_{\mathrm{c}} = 1.7 \pm 0.1 ~\mathrm{M_{\odot}}$, $T_{\mathrm{c}} = 11350 \pm 100 ~\mathrm{K}$ and $M_{\mathrm{h}} = 6.4 \pm 0.1 ~\mathrm{M_{\odot}}$, $T_{\mathrm{h}} = 20000 \pm 500 ~\mathrm{K}$, for the cooler and hotter star, respectively. We find broad single emission peaks in H$\alpha$ and H$\beta$ after subtracting the synthetic stellar spectra. Our results are consistent with mass exchange between the stars, and suggest the existence of a stream of gas being accreted onto the early B-type star., Comment: 14 pages with 13 figures

We present a photometric and spectroscopic study of HD 50526, an ellipsoidal binary member of the group Double Periodic Variable stars. Performing data-mining in photometric surveys and conducting new spectroscopic observations with several spectrographs during 2008 to 2015, we obtained orbital and stellar parameters of the system. The radial velocities were analyzed with the genetic PIKAIA algorithm, whereas Doppler tomography maps for the H$\alpha$ and H$\beta$ lines were constructed with the Total Variation Minimization code. An optimized simplex-algorithm was used to solve the inverse-problem adjusting the light curve with the best stellar parameters for the system. We find an orbital period of $6.701 \pm 0.001 ~\mathrm{d}$ and a long photometric cycle of $191 \pm 2 ~\mathrm{d}$. We detected the spectral features of the coldest star, and modeled it with a $\log{g} = 2.79 \pm 0.02 ~\mathrm{dex}$ giant of mass $1.13 \pm 0.02 ~\mathrm{M_{\odot}}$ and effective temperature $10500 \pm 125 ~\mathrm{K}$. In addition, we determine a mass ratio $q= 0.206 \pm 0.033$ and that the hot star is a B-type dwarf of mass $5.48 \pm 0.02 ~\mathrm{M_{\odot}}$. The $V$-band orbital light curve can be modeled including the presence of an accretion disk around the hotter star. This fills the Roche lobe of the hotter star, and has a radius $14.74 \pm 0.02 ~\mathrm{R_{\odot}}$ and temperature at the outer edge $9400 ~\mathrm{K}$. Two bright spots located in the disk account for the global morphology of the light curve. The Doppler tomography maps of H$\alpha$ and H$\beta$, reveal complex structures of mass fluxes in the system., Comment: 16 pages with 14 figures

Context. The different theoretical models concerning the formation of high-mass stars make distinct predictions regarding their progenitors, i.e. the high-mass prestellar cores. However, so far no conclusive observation of such objects has been made. Aims. We aim to study the very early stages of high-mass star formation in two infrared-dark, massive clumps, to identify the core population that they harbour. Methods. We obtained ALMA observations of continuum emission at 0.8mm and of the ortho-$\rm H_2D^+$ transition at 372GHz towards the two clumps. We use the SCIMES algorithm to identify cores in the position-position-velocity space, finding 16 cores. We model their observed spectra in the LTE approximation, deriving the centroid velocity, linewidth, and column density maps. We also study the correlation between the continuum and molecular data, which in general do not present the same structure. Results. We report for the first time the detection of ortho-$\rm H_2D^+$ in high-mass star-forming regions performed with an interferometer. The molecular emission shows narrow and subsonic lines, suggesting that locally the temperature of the gas is less than 10K. From the continuum emission we estimate the cores' total masses, and compare them with the respective virial masses. We also compute the volume density values, which are found to be higher than $10^{6}\, \rm cm^{-3}$. Conclusions. Our data confirm that ortho-$\rm H_2D^+$ is an ideal tracer of cold and dense gas. Interestingly, almost all the $\rm H_2D^+$-identified cores are less massive than 13M_sun , with the exception of one core in AG354. Furthermore, most of them are subvirial and larger than their Jeans masses. These results are difficult to explain in the context of the turbulent accretion models, which predict massive and virialised prestellar cores., Comment: Accepted for publications in A&A

Fragmentation often occurs in disk-like structures, both in the early Universe and in the context of present-day star formation. Supermassive black holes (SMBHs) are astrophysical objects whose origin is not well understood; they weigh millions of solar masses and reside in the centers of galaxies. An important formation scenario for SMBHs is based on collisions and mergers of stars in a massive cluster, in which the most massive star moves to the center of the cluster due to dynamical friction. This increases the rate of collisions and mergers since massive stars have larger collisional cross sections. This can lead to runaway growth of a very massive star which may collapse to become an intermediate-mass black hole. Here we investigate the dynamical evolution of Miyamoto-Nagai models that allow us to describe dense stellar clusters, including flattening and different degrees of rotation. We find that the collisions in these clusters depend mostly on the number of stars and the initial stellar radii for a given radial size of the cluster. By comparison, rotation seems to affect the collision rate by at most $20\%$. For flatness, we compared spherical models with systems that have a scale height of about $10\%$ of their radial extent, in this case finding a change in the collision rate of less than $25\%$. Overall, we conclude that the parameters only have a minor effect on the number of collisions. Our results also suggest that rotation helps to retain more stars in the system, reducing the number of escapers by a factor of $2-3$ depending on the model and the specific realization. After two million years, a typical lifetime of a very massive star, we find that about $630$ collisions occur in typical models with $N=10^4$, $R=100$ $\rm~R_\odot$ and a half-mass radius of $0.1$ $\rm~pc$, leading to a mass of about $6.3\times10^3$ $\rm~M_\odot$ for the most massive object., Comment: 10 pages, 7 figures

Runaway stellar collisions in dense star clusters are invoked to explain the presence of very massive stars or blue stragglers in the center of those systems. This process has also been explored for the first star clusters in the Universe and shown to yield stars that may collapse at some points into an intermediate mass black hole. Although the early evolution of star clusters requires the explicit modeling of the gas out of which the stars form, these calculations would be extremely time-consuming and often the effects of the gas can be accurately treated by including a background potential to account for the extra gravitational force. We apply this approximation to model the early evolution of the first dense star clusters formed in the Universe by performing $N$-body simulations, our goal is to understand how the additional gravitational force affects the growth of a very massive star through stellar mergers in the central parts of the star cluster. Our results show that the background potential increases the velocities of the stars, causing an overall delay in the evolution of the clusters and in the runaway growth of a massive star at the center. The population of binary stars is lower due to the increased kinetic energy of the stars, initially reducing the number of stellar collisions, and we show that relaxation processes are also affected. Despite these effects, the external potential enhances the mass of the merger product by a factor $\sim$2 if the collisions are maintained for long times., Comment: 16 pages. Accepted for publication in Astronomy and Astrophysics

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., Comment: 12 Pages, 8 Figures, 1 Table; Accepted for publication in MNRAS

Runaway collisions in dense clusters may lead to the formation of supermassive black hole (SMBH) seeds, and this process can be further enhanced by accretion, as recent models of SMBH seed formation in Population III star clusters have shown. This may explain the presence of supermassive black holes already at high redshift, $z>6$. However, in this context, mass loss during collisions was not considered and could play an important role for the formation of the SMBH seed. Here, we study the effect of mass loss, due to collisions of protostars, in the formation and evolution of a massive object in a dense primordial cluster. We consider both constant mass loss fractions as well as analytic models based on the stellar structure of the collision components. Our calculations indicate that mass loss can significantly affect the final mass of the possible SMBH seed. Considering a constant mass loss of 5% for every collision, we can lose between 60-80% of the total mass that is obtained if mass loss were not considered. Using instead analytical prescriptions for mass loss, the mass of the final object is reduced by 15-40%, depending on the accretion model for the cluster we study. Altogether, we obtain masses of the order of $10^4M_{\odot}$, which are still massive enough to be SMBH seeds., Comment: 12 pages, 9 figures, accepted by MNRAS

Disentangling the different stages of the star-formation process, in particular in the high-mass regime, is a challenge in astrophysics. Chemical clocks could help alleviating this problem, but their evolution strongly depends on many parameters, leading to degeneracy in the interpretation of the observational data. One of these uncertainties is the degree of CO depletion. We present here the first self-consistent magneto-hydrodynamic simulations of high-mass star-forming regions at different scales, fully coupled with a non-equilibrium chemical network, which includes C-N-O bearing molecules. Depletion and desorption processes are treated time-dependently. The results show that full CO-depletion (i.e. all gas-phase CO frozen-out on the surface of dust grains), can be reached very quickly, in one third or even smaller fractions of the free-fall time, whether the collapse proceeds on slow or fast timescales. This leads to a high level of deuteration in a short time both for typical tracers like N$_2$H$^+$, as well as for the main ion H$_3^+$, the latter being in general larger and more extended. N$_2$ depletion is slightly less efficient, and no direct effects on N-bearing molecules and deuterium fractionation are observed. We show that CO depletion is not the only driver of deuteration, and that there is a strong impact on $D_{frac}$ when changing the grain-size. We finally apply a two-dimensional gaussian Point Spread Function to our results to mimic observations with single-dish and interferometers. Our findings suggest that the low-values observed in high-mass star-forming clumps are in reality masking a full-depletion stage in the inner 0.1 pc region., Comment: accepted on ApJ

Context: In cold and dense gas prior to the formation of young stellar objects, heavy molecular species (including CO) are accreted onto dust grains. Under these conditions H$_3^+$ and its deuterated isotopologues become more abundant, enhancing the deuterium fraction of molecules such as N$_2$H$^+$ that are formed via ion-neutral reactions. Because this process is extremely temperature sensitive, the abundance of these species is likely linked to the evolutionary stage of the source. Aims: We investigate how the abundances of o-H$_2$D$^+$ and N$_2$D$^+$ vary with evolution in high-mass clumps. Methods: We observed with APEX the ground-state transitions of o-H$_2$D$^+$ near 372 GHz, and N$_2$D$^+$(3-2) near 231 GHz for three massive clumps in different evolutionary stages. The sources were selected within the G351.77-0.51 complex to minimise the variation of initial chemical conditions, and to remove distance effects. We modelled their dust continuum emission to estimate their physical properties, and also modelled their spectra under the assumption of local thermodynamic equilibrium to calculate beam-averaged abundances. Results: We find an anticorrelation between the abundance of o-H$_2$D$^+$ and that of N$_2$D$^+$, with the former decreasing and the latter increasing with evolution. With the new observations we are also able to provide a qualitative upper limit to the age of the youngest clump of about 10$^5$ yr, comparable to its current free-fall time. Conclusions: We can explain the evolution of the two tracers with simple considerations on the chemical formation paths, depletion of heavy elements, and evaporation from the grains. We therefore propose that the joint observation and the relative abundance of o-H$_2$D$^+$ and N$_2$D$^+$ can act as an efficient tracer of the evolutionary stages of the star-formation process.

Black holes with masses of $\rm 10^6-10^9~M_{\odot}$ dwell in the centers of most galaxies, but their formation mechanisms are not well known. A subdominant dissipative component of dark matter with similar properties to the ordinary baryons, known as mirror dark matter, may collapse to form massive black holes during the epoch of first galaxies formation. In this study, we explore the possibility of massive black hole formation via this alternative scenario. We perform three-dimensional cosmological simulations for four distinct halos and compare their thermal, chemical and dynamical evolution in both the ordinary and the mirror sectors. We find that the collapse of halos is significantly delayed in the mirror sector due to the lack of $\rm H_2$ cooling and only halos with masses above $ \rm \geq 10^7~ M_{\odot}$ are formed. Overall, the mass inflow rates are $\rm \geq 10^{-2}~M_{\odot}/yr$ and there is less fragmentation. This suggests that the conditions for the formation of massive objects, including black holes, are more favorable in the mirror sector.

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., Comment: 3 pages, 1 figure, 1 table. Submitted to BAAA (proceedings for the Binational AAA-SOCHIAS meeting 2018)

DQ Velorum is a galactic double periodic variable (DPV), this system is a semi-detached binary comprised of a B-type gainer and an A-type donor star plus an extended accretion disc around the gainer. The system also presents an orbital period of $6.08337$~days and a long period of $189$~days whose origin is still under debate. Here we studied the possibility that this period may be driven by a magnetic dynamo investigating the entire evolution of the system. The model matches in a very good way the current state of the system and it can potentially be used to describe the evolution of DQ Velorum. It also predicts an increase of the dynamo number of the donor during epochs of high mass transfer in this system, and a theoretical long/orbital period ratio very close to the observed one at the present system age., Comment: 3 pages, 1 figure, conference proceeding

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., Comment: 3 pages, 1 figure, proceeding of the Second Binational Meeting SOCHIAS-AAA

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., Comment: 3 pages, 2 figures

We present a simple model for the Double Periodic Variable (DPV) V495 Cen, which evolves as a binary system of intermediate mass, where the gainer cannot accrete at high rate, limited by the Eddington accretion rate, leading to the formation of an accretion disc. The theoretical model begins at the zero age main sequence considering the rotation for both stars. For this purpose we used the stellar evolution code MESA, developed to calculate the evolution of stars in a wide range of parameters. We started the model adjusting fundamental parameters published for this system through a chi-square optimization algorithm, and adopting an initial orbital period of 3.9 days and initial masses for the primary component $M_{i,d} = 3.40$ $M_{\odot}$ and $M_{i,g}= 3.18$ $M_{\odot}$ for the gainer, with a metallicity associated to this type of DPV of $Z = 0.02$. The method converged successfully for eight free degrees and 5% of confidence with a chi-square of $\Delta \chi^{2}_{0.95,8}= 0.212$. We describe each evolutionary stage of both components until that the donor reaches 20% core helium depletion as stop criterion. We offer a complementary analysis for understanding the mechanism of the magnetic dynamos inside the donor star using the Tayler-Spruit formalism. Currently, the theoretical model is consistent with the fundamental parameters published for V495 Cen and we discuss how our predictions can help to develop efficients theoretical models for DPV stars.

Since its proposal in 1992, the Applegate mechanism has been discussed as a potential intrinsical mechanism to explain transit timing variations in various kinds of close binary systems. Most analytical arguments presented so far focused on the energetic feasibility of the mechanism, while applying rather crude one- or two-zone prescriptions to describe the exchange of angular momentum within the star. In this paper, we present the most detailed approach to date to describe the physics giving rise to the modulation period from kinetic and magnetic fluctuations. Assuming moderate levels of stellar parameter fluctuations, we find that the resulting binary period variations are one or two orders of magnitude lower than the observed values in RS-CVn like systems, supporting the conclusion of existing theoretical work that the Applegate mechanism may not suffice to produce the observed variations in these systems. The most promising Applegate candidates are low-mass post-common-envelope binaries (PCEBs) with binary separations $\lesssim 1~\mathrm{R}_\odot$ and secondary masses in the range of $0.30~\mathrm{M}_\odot$ and $0.36~\mathrm{M}_\odot$., Comment: 10 pages, 8 figures. Accepted for publication in A&A

Although timing variations in close binary systems have been studied for a long time, their underlying causes are still unclear. A possible explanation is the so-called Applegate mechanism, where a strong, variable magnetic field can periodically change the gravitational quadrupole moment of a stellar component, thus causing observable period changes. One of the systems exhibiting such strong orbital variations is the RS CVn binary HR 1099, whose activity cycle has been studied by various authors via photospheric and chromospheric activity indicators, resulting in contradicting periods. We aim at independently determining the magnetic activity cycle of HR 1099 using archival X-ray data to allow for a comparison to orbital period variations. Archival X-ray data from 80 different observations of HR 1099 acquired with 12 different X-ray facilities and covering almost four decades were used to determine X-ray fluxes in the energy range of 2-10 keV via spectral fitting and flux conversion. Via the Lomb-Scargle periodogram we analyze the resulting long-term X-ray light curve to search for periodicities. We do not detect any statistically significant periodicities within the X-ray data. An analysis of optical data of HR 1099 shows that the derivation of such periods is strongly dependent on the time coverage of available data, since the observed optical variations strongly deviate from a pure sine wave. We argue that this offers an explanation as to why other authors derive such a wide range of activity cycle periods based on optical data. We conclude that our analysis constitutes the longest stellar X-ray activity light curve acquired to date, yet the still rather sparse sampling of the X-ray data, along with stochastic flaring activity, does not allow for the independent determination of an X-ray activity cycle., Comment: 8 pages, 6 figures, 2 tables accepted for publication in A&A

Eclipsing time variations (ETVs) are observed in many close binary systems. In particular, for several post-common-envelope binaries (PCEBs) that consist of a white dwarf and a main sequence star, the O-C diagram suggests that real or apparent orbital period variations are driven by Jupiter-mass planets or as a result of magnetic activity, the so-called Applegate mechanism. The latter explains orbital period variations as a result of changes in the stellar quadrupole moment due to magnetic activity. We explore the feasibility of driving ETVs via the Applegate mechanism for a sample of PCEB systems, including a range of different rotations. Using the MESA code we evolve 12 stars with different masses and rotation rates. We apply a simple dynamo model to their radial profiles to investigate on which scale the predicted activity cycle matches the observed modulation period, and quantify the uncertainty, and further calculate the required energies to drive que Applegate mechanism. We show that the Applegate mechanism is energetically feasible in 5 PCEB systems, and note that these are the systems with the highest rotation rate compared to the critical rotation rate of the main-sequence star. The results suggest that the ratio of physical to critical rotation in the main sequence star is an important indicator for the feasibility of Applegate's mechanism, but exploring larger samples will be necessary to probe this hypothesis., Comment: 9 pages, 5 figures. Accepted for publication in A&A

We present elemental abundance analysis of high-resolution spectra for five giant stars, deriving Fe, Mg, Al, C, N, O, Si and Ce abundances, and spatially located within the innermost regions of the bulge globular cluster NGC 6522, based on H-band spectra taken with the multi-object APOGEE-north spectrograph from the SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. Of the five cluster candidates, two previously unremarked stars are confirmed to have second-generation (SG) abundance patterns, with the basic pattern of depletion in C and Mg simultaneous with enrichment in N and Al as seen in other SG globular cluster populations at similar metallicity. } In agreement with the most recent optical studies, the NGC 6522 stars analyzed exhibit (when available) only mild overabundances of the s-process element Ce, contradicting the idea of the NGC 6522 stars being formed from gas enriched by spinstars and indicating that other stellar sources such as massive AGB stars could be the primary intra-cluster medium polluters. The peculiar abundance signature of SG stars have been observed in our data, confirming the presence of multiple generations of stars in NGC 6522., Comment: 13 pages, 8 figures, 4 tables, matches the accepted version in A&A

Models aiming to explain the formation of massive black hole seeds, and in particular the direct collapse scenario, face substantial difficulties. These are rooted in rather ad hoc and fine-tuned initial conditions, such as the simultaneous requirements of extremely low metallicities and strong radiation backgrounds. Here we explore a modification of such scenarios where a massive primordial star cluster is initially produced. Subsequent stellar collisions give rise to the formation of massive (10^4 - 10^5 solar mass) objects. Our calculations demonstrate that the interplay between stellar dynamics, gas accretion and protostellar evolution is particularly relevant. Gas accretion onto the protostars enhances their radii, resulting in an enhanced collisional cross section. We show that the fraction of collisions can increase from 0.1-1% of the initial population to about 10% when compared to gas-free models or models of protostellar clusters in the local Universe. We conclude that very massive objects can form in spite of initial fragmentation, making the first massive protostellar clusters viable candidate birth places for observed supermassive black holes., Comment: Accepted in MNRAS

Collisions were suggested to potentially play a role in the formation of massive stars in present day clusters, and have likely been relevant during the formation of massive stars and intermediate mass black holes within the first star clusters. In the early Universe, the first stellar clusters were particularly dense, as fragmentation typically only occurred at densities above $10^9$cm$^{-3}$, and the radii of the protostars were enhanced due to the larger accretion rates, suggesting a potentially more relevant role of stellar collisions. We present here a detailed parameter study to assess how the number of collisions as well as the mass growth of the most massive object depends on the properties of the cluster, and we characterize the time evolution with three effective parameters, the time when most collisions occur, the duration of the collisions period, as well as the normalization required to obtain the total number of collisions. We apply our results to typical Population III (Pop.III) clusters of about $1000$M$_\odot$, finding that a moderate enhancement of the mass of the most massive star by a factor of a few can be expected. For more massive Pop.III clusters as expected in the first atomic cooling halos, we expect a more significant enhancement by a factor of $15-32$. We therefore conclude that collisions in massive Pop.III clusters were likely relevant to form the first intermediate mass black holes., Comment: 9 pages, 8 figures, accepted by A&A

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., Comment: 16 pages, 13 figures, 4 tables. Accepted for publication with A&A

We report the peculiar chemical abundance patterns of eleven atypical Milky Way (MW) field red giant stars observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE). These atypical giants exhibit strong Al and N enhancements accompanied by C and Mg depletions, strikingly similar to those observed in the so-called second-generation (SG) stars of globular clusters (GCs). Remarkably, we find low-Mg abundances ([Mg/Fe]$<$0.0) together with strong Al and N overabundances in the majority (5/7) of the metal-rich ([Fe/H]$\gtrsim - 1.0$) sample stars, which is at odds with actual observations of SG stars in Galactic CGs of similar metallicities. This chemical pattern is unique and unprecedented among MW stars, posing urgent questions about its origin. These atypical stars could be former SG stars of dissolved GCs formed with intrinsically lower abundances of Mg and enriched Al (subsequently self-polluted by massive AGB stars) or the result of exotic binary systems. We speculate that the stars Mg-deficiency as well as the orbital properties suggest that they could have an extragalactic origin. This discovery should guide future dedicated spectroscopic searches of atypical stellar chemical patterns in our Galaxy; a fundamental step forward to understand the Galactic formation and evolution., Comment: Accepted for publication in ApJ Letters, 11 pages, 2 Figures, 2 Tables

We have observed 104 gravitationally-lensed quasars at $z\sim1-4$ with Herschel/SPIRE, the largest such sample ever studied. By targeting gravitational lenses, we probe intrinsic far-infrared (FIR) luminosities and star formation rates (SFRs) more typical of the population than the extremely luminous sources that are otherwise accessible. We detect 72 objects with Herschel/SPIRE and find 66 percent (69 sources) of the sample have spectral energy distributions (SEDs) characteristic of dust emission. For 53 objects with sufficiently constrained SEDs, we find a median effective dust temperature of $38^{+12}_{-5}$ K. By applying the radio-infrared correlation, we find no evidence for an FIR excess which is consistent with star-formation-heated dust. We derive a median magnification-corrected FIR luminosity of $3.6^{+4.8}_{-2.4}~\times 10^{11}~{\rm L_{\odot}}$ and median SFR of $120^{+160}_{-80}~{\rm M_{\odot}~yr^{-1}}$ for 94 quasars with redshifts. We find $\sim10$ percent of our sample have FIR properties similar to typical dusty star-forming galaxies at $z\sim2-3$ and a range of SFRs $<20-10000~{\rm M_{\odot}~yr^{-1}}$ for our sample as a whole. These results are in line with current models of quasar evolution and suggests a coexistence of dust-obscured star formation and AGN activity is typical of most quasars. We do not find a statistically-significant difference in the FIR luminosities of quasars in our sample with a radio excess relative to the radio-infrared correlation. Synchrotron emission is found to dominate at FIR wavelengths for $<15$ percent of those sources classified as powerful radio galaxies., Comment: 47 pages, 89 figures, accepted for publication in MNRAS

Even though compressible plasma turbulence is encountered in many astrophysical phenomena, its effect is often not well understood. Furthermore, direct numerical simulations are typically not able to reach the extreme parameters of these processes. For this reason, large-eddy simulations (LES), which only simulate large and intermediate scales directly, are employed. The smallest, unresolved scales and the interactions between small and large scales are introduced by means of a subgrid-scale (SGS) model. We propose and verify a new set of nonlinear SGS closures for future application as an SGS model in LES of compressible magnetohydrodynamics (MHD). We use 15 simulations (without explicit SGS model) of forced, isotropic, homogeneous turbulence with varying sonic Mach number $\mathrm{M_s} = 0.2$ to $20$ as reference data for the most extensive \textit{a priori} tests performed so far in literature. In these tests we explicitly filter the reference data and compare the performance of the new closures against the most widely tested closures. These include eddy-viscosity and scale-similarity type closures with different normalizations. Performance indicators are correlations with the turbulent energy and cross-helicty flux, the average SGS dissipation, the topological structure and the ability to reproduce the correct magnitude and direction of the SGS vectors. We find that only the new nonlinear closures exhibit consistently high correlations (median value \textgreater$0.8$) with the data over the entire parameter space and outperform the other closures in all tests. Moreover, we show that these results are independent of resolution and chosen filter scale. Additionally, the new closures are effectively coefficient-free with a deviation of less than $20\%$., Comment: 17 pages, 10 figures, 2 tables; to be published in Physics of Plasmas

Dust plays a key role in the evolution of the ISM and its correct modelling in numerical simulations is therefore fundamental. We present a new and self-consistent model that treats grain thermal coupling with the gas, radiation balance, and surface chemistry for molecular hydrogen. This method can be applied to any dust distribution with an arbitrary number of grain types without affecting the overall computational cost. In this paper we describe in detail the physics and the algorithm behind our approach, and in order to test the methodology, we present some examples of astrophysical interest, namely (i) a one-zone collapse with complete gas chemistry and thermochemical processes, (ii) a 3D model of a low-metallicity collapse of a minihalo starting from cosmological initial conditions, and (iii) a turbulent molecular cloud with H-C-O chemistry (277 reactions), together with self-consistent cooling and heating solved on the fly. Although these examples employ the publicly available code KROME, our approach can be easily integrated into any computational framework., Comment: 17 pages, 12 figures, MNRAS accepted

Understanding the formation of the extremely metal poor star SDSS-J102915+172927 is of fundamental importance to improve our knowledge on the transition between the first and second generation of stars in the Universe. In this paper, we perform three-dimensional cosmological hydrodynamical simulations of dust-enriched halos during the early stages of the collapse process including a detailed treatment of the dust physics. We employ the astrochemistry package \krome coupled with the hydrodynamical code \textsc{enzo} assuming grain size distributions produced by the explosion of core-collapse supernovae of 20 and 35 M$_\odot$ primordial stars which are suitable to reproduce the chemical pattern of the SDSS-J102915+172927 star. We find that the dust mass yield produced from Population III supernovae explosions is the most important factor which drives the thermal evolution and the dynamical properties of the halos. Hence, for the specific distributions relevant in this context, the composition, the dust optical properties, and the size-range have only minor effects on the results due to similar cooling functions. We also show that the critical dust mass to enable fragmentation provided by semi-analytical models should be revised, as we obtain values one order of magnitude larger. This determines the transition from disk fragmentation to a more filamentary fragmentation mode, and suggests that likely more than one single supernova event or efficient dust growth should be invoked to get such a high dust content., Comment: Accepted on ApJ

The observed eclipsing time variations in post-common-envelope binaries (PCEBs) can be interpreted as potential evidence for massive Jupiter-like planets, or as a result of magnetic activity, leading to quasi-periodic changes in the quadrupole moment of the secondary star. The latter is commonly referred to as the Applegate mechanism. Following Brinkworth et al. (2006), we employ here an improved version of Applegate's model including the angular momentum exchange between a finite shell and the core of the star. The framework is employed to derive the general conditions under which the Applegate mechanism can work, and is subsequently applied to a sample of 16 close binary systems with potential planets, including 11 PCEBs. Further, we present a detailed derivation and study of analytical models which allow for an straightforward extension to other systems. Using our full numerical framework, we show that the Applegate mechanism can clearly explain the observed eclipsing time variations in 4 of the systems, while the required energy to produce the quadrupole moment variations is too high in at least 8 systems. In the remaining 4 systems, the required energy is comparable to the available energy produced by the star, which we consider as borderline cases. Therefore, the Applegate mechanism cannot uniquely explain the observed period time variations for this entire population. Even in systems where the required energy is too high, the Applegate mechanism may provide an additional scatter, which needs to be considered in the derivation and analysis of planetary models., Comment: 12 pages, 7 figures. Accepted for publication in A&A

Observations of quasars at $ z > 6$ suggest the presence of black holes with a few times $\rm 10^9 ~M_{\odot}$. Numerous models have been proposed to explain their existence including the direct collapse which provides massive seeds of $\rm 10^5~M_{\odot}$. The isothermal direct collapse requires a strong Lyman-Werner flux to quench $\rm H_2$ formation in massive primordial halos. In this study, we explore the impact of trace amounts of metals and dust enrichment. We perform three dimensional cosmological simulations for two halos of $\rm > 10^7~M_{\odot}$ with $\rm Z/Z_{\odot}= 10^{-4}-10^{-6}$ illuminated by an intense Lyman Werner flux of $\rm J_{21}=10^5$. Our results show that initially the collapse proceeds isothermally with $\rm T \sim 8000$ K but dust cooling becomes effective at densities of $\rm 10^{8}-10^{12} ~cm^{-3}$ and brings the gas temperature down to a few 100-1000 K for $\rm Z/Z_{\odot} \geq 10^{-6}$. No gravitationally bound clumps are found in $\rm Z/Z_{\odot} \leq 10^{-5}$ cases by the end of our simulations in contrast to the case with $\rm Z/Z_{\odot} = 10^{-4}$. Large inflow rates of $\rm \geq 0.1~M_{\odot}/yr$ are observed for $\rm Z/Z_{\odot} \leq 10^{-5}$ similar to a zero-metallicity case while for $\rm Z/Z_{\odot} = 10^{-4}$ the inflow rate starts to decline earlier due to the dust cooling and fragmentation. For given large inflow rates a central star of $\rm \sim 10^4~M_{\odot}$ may form for $\rm Z/Z_{\odot} \leq 10^{-5}$., Comment: Accepted for publication in ApJ, comments are still welcome

Massive primordial halos exposed to moderate UV backgrounds are the potential birthplaces of supermassive black holes. In such a halo, an initially isothermal collapse will occur, leading to high accretion rates of $\sim0.1$~M$_\odot$~yr$^{-1}$. During the collapse, the gas in the interior will turn into a molecular state, and form an accretion disk due to the conservation of angular momentum. We consider here the structure of such an accretion disk and the role of viscous heating in the presence of high accretion rates for a central star of $10$, $100$ and $10^4$~M$_\odot$. Our results show that the temperature in the disk increases considerably due to viscous heating, leading to a transition from the molecular to the atomic cooling phase. We found that the atomic cooling regime may extend out to several $100$~AU for a $10^4$~M$_\odot$ central star and provides substantial support to stabilize the disk. It therefore favors the formation of a massive central object. The comparison of clump migration and contraction time scales shows that stellar feedback from these clumps may occur during the later stages of the evolution. Overall, viscous heating provides an important pathway to obtain an atomic gas phase within the center of the halo, and helps in the formation of very massive objects. The latter may collapse to form a massive black hole of about $\geq 10^4$~M$_\odot$., Comment: Accepted for publication in Astronomy & Astrophysics, comments are still welcome

Our understanding of population III star formation is still in its infancy. They are formed in dark matter minihalos of $\rm 10^5-10^6 M_{\odot}$ at $z=20-30$. Recent high resolution cosmological simulations show that a protostellar disk forms as a consequence of gravitational collapse and fragments into multiple clumps. However, it is not entirely clear if these clumps will be able to survive to form multiple stars as simulations are unable to follow the disk evolution for longer times. In this study, we employ a simple analytical model to derive the properties of marginally stable steady-state disks. Our results show that the stability of the disk depends on the critical value of the viscous parameter $\alpha$. For $\alpha_{crit} = 1$, the disk is stable for an accretion rate of $\rm \leq 10^{-3} M_{\odot}/yr$ and becomes unstable at radii about $\rm \geq 100 AU$ in the presence of an accretion rate of $\rm 10^{-2} M_{\odot}/yr$. For $0.06 < \alpha_{crit} < 1$, the disk can be unstable for both accretion rates. The comparison of the migration and the Kelvin-Helmholtz time scales shows that clumps are expected to migrate inward before reaching the main sequence. Furthermore, in the presence of a massive central star the clumps within the central 1 AU will be tidally disrupted. We also find that UV feedback from the central star is unable to disrupt the disk, and that photo-evaporation becomes important only once the accretion rate has dropped to $\rm 2 \times 10^{-4} M_{\odot}/yr$. As a result, the central star may reach a mass of 100 $\rm M_{\odot}$ or even higher., Comment: Accepted for publication in MNRAS, comments are still wlecome

Observations of high redshift quasars at $z>6$ indicate that they harbor supermassive black holes (SMBHs) of a billion solar masses. The direct collapse scenario has emerged as the most plausible way to assemble SMBHs. The nurseries for the direct collapse black holes are massive primordial halos illuminated with an intense UV flux emitted by population II (Pop II) stars. In this study, we compute the critical value of such a flux ($J_{21}^{\rm crit}$) for realistic spectra of Pop II stars through three-dimensional cosmological simulations. We derive the dependence of $J_{21}^{\rm crit}$ on the radiation spectra, on variations from halo to halo, and on the impact of X-ray ionization. Our findings show that the value of $J_{21}^{\rm crit}$ is a few times $\rm 10^4$ and only weakly depends on the adopted radiation spectra in the range between $T_{\rm rad}=2 \times 10^4-10^5$ K. For three simulated halos of a few times $\rm 10^{7}$~M$_{\odot}$, $J_{21}^{\rm crit}$ varies from $\rm 2 \times 10^4 - 5 \times 10^4$. The impact of X-ray ionization is almost negligible and within the expected scatter of $J_{21}^{\rm crit}$ for background fluxes of $J_{\rm X,21} \leq 0.1$. The computed estimates of $J_{21}^{\rm crit}$ have profound implications for the quasar abundance at $z=10$ as it lowers the number density of black holes forming through an isothermal direct collapse by a few orders of magnitude below the observed black holes density. However, the sites with moderate amounts of $\rm H_2$ cooling may still form massive objects sufficient to be compatible with observations., Comment: Accepted for publication in MNRAS, comments are welcome

Context. The seeds of the first supermassive black holes may have resulted from the direct collapse of hot primordial gas in $\gtrsim 10^4$ K haloes, forming a supermassive or quasistar as an intermediate stage. Aims. We explore the formation of a protostar resulting from the collapse of primordial gas in the presence of a strong Lyman-Werner radiation background. Particularly, we investigate the impact of turbulence and rotation on the fragmentation behaviour of the gas cloud. We accomplish this goal by varying the initial turbulent and rotational velocities. Methods. We performed 3D adaptive mesh refinement simulations with a resolution of 64 cells per Jeans length using the ENZO code, simulating the formation of a protostar up to unprecedentedly high central densities of $10^{21}$ cm$^{-3}$, and spatial scales of a few solar radii. To achieve this goal, we employed the KROME package to improve modelling of the chemical and thermal processes. Results. We find that the physical properties of the simulated gas clouds become similar on small scales, irrespective of the initial amount of turbulence and rotation. After the highest level of refinement was reached, the simulations have been evolved for an additional ~5 freefall times. A single bound clump with a radius of $2 \times 10^{-2}$ AU and a mass of ~$7 \times 10^{-2}$ M$_{\odot}$ is formed at the end of each simulation, marking the onset of protostar formation. No strong fragmentation is observed by the end of the simulations, regardless of the initial amount of turbulence or rotation, and high accretion rates of a few solar masses per year are found. Conclusions. Given such high accretion rates, a quasistar of $10^5$ M$_{\odot}$ is expected to form within $10^5$ years., Comment: 18 pages, 7 figures, fixed typos, added references and clarified some details; accepted for publication in A&A

We study the Initial Mass Function (IMF) and host halo properties of Intermediate Mass Black Holes (IMBH, 10^{4-6} Msun) formed inside metal-free, UV illuminated atomic cooling haloes (virial temperature T_vir > 10^4 K) either via the direct collapse of the gas or via an intermediate Super Massive Star (SMS) stage. We achieve this goal in three steps: (a) we derive the gas accretion rate for a proto-SMS to undergo General Relativity instability and produce a direct collapse black hole (DCBH) or to enter the ZAMS and later collapse into a IMBH; (b) we use merger-tree simulations to select atomic cooling halos in which either a DCBH or SMS can form and grow, accounting for metal enrichment and major mergers that halt the growth of the proto-SMS by gas fragmentation. We derive the properties of the host halos and the mass distribution of black holes at this stage, and dub it the "Birth Mass Function"; (c) we follow the further growth of the DCBH due to accretion of leftover gas in the parent halo and compute the final IMBH mass.We consider two extreme cases in which minihalos (T_vir < 10^4 K) can (fertile) or cannot (sterile) form stars and pollute their gas leading to a different IMBH IMF. In the (fiducial) fertile case the IMF is bimodal extending over a broad range of masses, M= (0.5-20)x10^5 Msun, and the DCBH accretion phase lasts from 10 to 100 Myr. If minihalos are sterile, the IMF spans the narrower mass range M= (1-2.8)x10^6 Msun, and the DCBH accretion phase is more extended (70-120 Myr). We conclude that a good seeding prescription is to populate halos (a) of mass 7.5 < log (M_h/Msun) < 8, (b) in the redshift range 8 < z < 17, (c) with IMBH in the mass range 4.75 < log (M_BH/Msun) < 6.25., Comment: MNRAS, in press. Comments welcome

Recent discoveries of carbon-enhanced metal-poor stars like SMSS J031300.36-670839.3 provide increasing observational insights into the formation conditions of the first second-generation stars in the Universe, reflecting the chemical conditions after the first supernova explosion. Here, we present the first cosmological simulations with a detailed chemical network including primordial species as well as C, C$^+$, O, O$^+$, Si, Si$^+$, and Si$^{2+}$ following the formation of carbon-enhanced metal poor stars. The presence of background UV flux delays the collapse from $z=21$ to $z=15$ and cool the gas down to the CMB temperature for a metallicity of Z/Z$_\odot$=10$^{-3}$. This can potentially lead to the formation of lower mass stars. Overall, we find that the metals have a stronger effect on the collapse than the radiation, yielding a comparable thermal structure for large variations in the radiative background. We further find that radiative backgrounds are not able to delay the collapse for Z/Z$_\odot$=10$^{-2}$ or a carbon abundance as in SMSS J031300.36-670839.3., Comment: submitted to ApJL