Showing 1–50 of 155 results for author: Schleicher, D R G

Supermassive black hole formation via collisions in black hole clusters

Authors: Benjamin Gaete, Dominik R. G. Schleicher, Alessandro Lupi, Bastian Reinoso, Michael Fellhauer, Marcelo C. Vergara

Abstract: More than 300 supermassive black holes have been detected at redshifts larger than 6, and they are abundant in the centers of local galaxies. Their formation mechanisms, however, are still rather unconstrained. A possible origin of these supermassive black holes could be through mergers in dense black hole clusters, forming as a result of mass segregation within nuclear star clusters in the center… ▽ More More than 300 supermassive black holes have been detected at redshifts larger than 6, and they are abundant in the centers of local galaxies. Their formation mechanisms, however, are still rather unconstrained. A possible origin of these supermassive black holes could be through mergers in dense black hole clusters, forming as a result of mass segregation within nuclear star clusters in the center of galaxies. In this study, we present the first systematic investigation of the evolution of such black hole clusters where the effect of an external potential is taken into account. Such a potential could be the result of gas inflows into the central region, for example as a result of galaxy mergers. We show here that the efficiency for the formation of a massive central object is mostly regulated by the ratio of cluster velocity dispersion divided by the speed of light, potentially reaching efficiencies of 0.05-0.08 in realistic systems. Our results show that this scenario is potentially feasible and may provide seeds black hole of at least 10^3 solar masses. We conclude that the formation of seed black holes via this channel should be taken into account in statistical assessments of the black hole population. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: 15 pages, 11 figures, 1 table. Accepted for publication with A&A

Turbulence and the characteristics of circumstellar discs

Authors: R. Riaz, D. R. G. Schleicher, S. Vanaverbeke, R. Klessen, J. Saavedra-Bastidas

Abstract: 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 exte… ▽ More 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. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Efficiency of black hole formation via collisions in stellar systems: An analysis of data from simulations and observations

Authors: M. C. Vergara, D. R. G. Schleicher, A. Escala, B. Reinoso, F. Flammini Dotti, A. W. H. Kamlah, M. Liempi, N. Hoyer, N. Neumayer, R. Spurzem

Abstract: 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 syst… ▽ More 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 ($ε_{\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 ($ε_{\mathrm{BH}}$) is shown to be influenced by the ratio of the initial stellar mass to the critical mass. △ Less

Submitted 11 June, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: 16 pages, 5 figures

Black Hole shadows of $α'$-corrected black holes

Authors: F. Agurto-Sepúlveda, J. Oliva, M. Oyarzo, D. R. G Schleicher

Abstract: 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 $α'$, including also the dilaton. We provide a detailed characterization of the geometry, as well as the innermost stable circular orbit (ISCO) and phot… ▽ More 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 $α'$, 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 $α'$ 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. △ Less

Submitted 24 June, 2024; v1 submitted 19 April, 2024; originally announced April 2024.

Comments: 14 pages, 8 figures. V2: small changes, references added, accepted in PRD

Magnetic field amplification in massive primordial halos: Influence of Lyman-Werner radiation

Authors: V. B. Díaz, D. R. G. Schleicher, M. A. Latif, P. Grete, R. Banerjee

Abstract: 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… ▽ More 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. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: 20 pages, 16 figures, 3 tables. Accepted for publication in A&A

Magnetic braking during direct collapse black hole formation

Authors: Muhammad A. Latif, Dominik R. G. Schleicher

Abstract: Magnetic fields are expected to be efficiently amplified during the formation of the first massive black holes via the small-scale dynamo and in the presence of strong accretion shocks occurring during gravitational collapse. Here, we analyze high-resolution cosmological magneto-hydrodynamical simulations of gravitational collapse in atomic cooling halos, exploring the dynamical role of magnetic f… ▽ More Magnetic fields are expected to be efficiently amplified during the formation of the first massive black holes via the small-scale dynamo and in the presence of strong accretion shocks occurring during gravitational collapse. Here, we analyze high-resolution cosmological magneto-hydrodynamical simulations of gravitational collapse in atomic cooling halos, exploring the dynamical role of magnetic fields, particularly concerning the effect of magnetic braking and angular momentum transport. We find that after the initial amplification, magnetic fields contribute to the transport of angular momentum and reduce it compared to pure hydrodynamical simulations. However, the magnetic and Reynolds torques do not fully compensate for the inward advection of angular momentum, which still accumulates over timescales of $\sim1$~Myr. A Jeans analysis further shows that magnetic pressure strongly contributes to suppressing fragmentation on scales of $0.1-10$~pc. Overall, the presence of magnetic fields thus aids in the transport of angular momentum and favors the formation of massive objects. △ Less

Submitted 2 July, 2023; originally announced July 2023.

Comments: Accepted for Publication in ApJL, comments are welcome

Survival of fossil fields during the pre-main sequence evolution of intermediate-mass stars

Authors: Dominik R. G. Schleicher, Juan Pablo Hidalgo, Daniele Galli

Abstract: Chemically peculiar Ap and Bp stars host strong large-scale magnetic fields in the range of $200$~G up to $30$~kG, which are often considered to be the origin of fossil magnetic fields. We assess the evolution of such fossil fields during the star formation process and the pre-main sequence evolution of intermediate stars, considering fully convective models, models including a transition to a rad… ▽ More Chemically peculiar Ap and Bp stars host strong large-scale magnetic fields in the range of $200$~G up to $30$~kG, which are often considered to be the origin of fossil magnetic fields. We assess the evolution of such fossil fields during the star formation process and the pre-main sequence evolution of intermediate stars, considering fully convective models, models including a transition to a radiative protostar and models with a radiative core. We also examine the implications of the interaction between the fossil field and the core dynamo. We employ analytic and semi-analytic calculations combined with current observational constraints. For fully convective models, we show that magnetic field decay via convection can be expected to be very efficient for realistic parameters of turbulent resistivities. Based on the observed magnetic field strength - density relation, as well as the expected amount of flux loss due to ambipolar diffusion, it appears unlikely that convection could be suppressed via strong enough magnetic fields. On the other hand, a transition from a convective to a radiative core could very naturally explain the survival of a significant amount of flux, along with the presence of a critical mass. We show that in some cases, the interaction of a fossil field with a core dynamo may further lead to changes in the surface magnetic field structure. In the future, it will be important to understand in more detail how the accretion rate evolves as a function of time during the formation of intermediate-mass protostars, including its impact on the protostellar structure. The latter may even allow to derive quantitative predictions concerning the expected population of large scale magnetic fields in radiative stars. △ Less

Submitted 29 June, 2023; originally announced June 2023.

Comments: 10 pages, 6 figures, 3 tables. Accepted for publication with A&A

Origin of supermassive black holes: predictions for the black hole population

Authors: M. Liempi, L. Almonacid, D. R. G. Schleicher, A. Escala

Abstract: 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 t… ▽ More 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. △ Less

Submitted 29 May, 2023; originally announced May 2023.

Comments: 3 pages, 3 figures

Simulations of dynamo action in slowly rotating M dwarfs: Dependence on dimensionless parameters

Authors: C. A. Ortiz-Rodríguez, P. J. Käpylä, F. H. Navarrete, D. R. G Schleicher, R. E. Mennickent, J. P. Hidalgo, B. Toro

Abstract: 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 ma… ▽ More 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. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: 11 pages, 12 figures, submitted to A&A

Effects of the centrifugal force in stellar dynamo simulations

Authors: Felipe H. Navarrete, Petri J. Käpylä, Dominik R. G. Schleicher, Robi Banerjee

Abstract: The centrifugal force is often omitted from simulations of stellar convection either for numerical reasons or because it is assumed to be weak compared to the gravitational force. However, it might be an important factor in rapidly rotating stars, such as solar analogs, due to its $Ω^2$ scaling, where $Ω$ is the rotation rate of the star. We study the effects of the centrifugal force in a set of 2… ▽ More The centrifugal force is often omitted from simulations of stellar convection either for numerical reasons or because it is assumed to be weak compared to the gravitational force. However, it might be an important factor in rapidly rotating stars, such as solar analogs, due to its $Ω^2$ scaling, where $Ω$ is the rotation rate of the star. We study the effects of the centrifugal force in a set of 21 semi-global stellar dynamo simulations with varying rotation rates. Included in the set are three control runs aimed at distinguishing the effects of the centrifugal force from the nonlinear evolution of the solutions. We decomposed the magnetic field into spherical harmonics and studied the migration of azimuthal dynamo waves (ADWs), the energy of different large-scale magnetic modes, and differential rotation. In the regime with the lowest rotation rates, $Ω= 5-10Ω_\odot$, where $Ω_\odot$ is the rotation rate of the Sun, we see no marked changes in either the differential rotation or the magnetic field properties. For intermediate rotation, $Ω= 20-25Ω_\odot$, we identify an increase in the differential rotation as a function of centrifugal force. The axisymmetric magnetic energy tends to decrease with centrifugal force, while the non-axisymmetric one increases. The ADWs are also affected, especially in the propagation direction. In the most rapidly rotating set with $Ω=30Ω_\odot$, these changes are more pronounced, and in one case the propagation direction of the ADW changes from prograde to retrograde. The control runs suggest that the results are a consequence of the centrifugal force and not due to the details of the initial conditions or the history of the run. We find that the differential rotation and the ADWs only change as a function of the centrifugal force when rotation is rapid enough. △ Less

Submitted 14 August, 2023; v1 submitted 2 May, 2023; originally announced May 2023.

Comments: 8 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 678, A9 (2023)

Origin of magnetism in early-type stars

Authors: J. P. Hidalgo, P. J. Käpylä, C. A. Ortiz-Rodríguez, F. H. Navarrete, B. Toro, D. R. G. Schleicher

Abstract: 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 s… ▽ More 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. △ Less

Submitted 19 March, 2023; originally announced March 2023.

Comments: Conference proceeding (Boletín de la Asociación Argentina de Astronomía), 3 pages, 3 figures

The physical origin of super-competitive accretion during the formation of the first supermassive black holes

Authors: Dominik R. G. Schleicher, Bastián Reinoso, Ralf S. Klessen

Abstract: Numerical simulations have shown the occurence of a scenario termed ''super-competitive accretion'', a term that describes a situation where only the central few objects grow supermassive while a larger number of stars compete for the reservoir, with significant accretion flows of $\gtrsim0.1$ M$_\odot$ yr$^{-1}$. This scenario particularly implies that the presence of fragmentation will not neces… ▽ More Numerical simulations have shown the occurence of a scenario termed ''super-competitive accretion'', a term that describes a situation where only the central few objects grow supermassive while a larger number of stars compete for the reservoir, with significant accretion flows of $\gtrsim0.1$ M$_\odot$ yr$^{-1}$. This scenario particularly implies that the presence of fragmentation will not necessarily impeed the formation of a central massive object. We here explore this phenomenon using analytical estimates for growth via collisions and accretion, considering accretion due to self-gravity as well as Bondi-Hoyle accretion. Particularly, we explore under what conditions the accretion onto the central massive object breaks down, and derive a criterion that depends on the mass of the most massive object and the mass in fragments. For compact clusters with sizes about $0.1$~pc, we further find that the mass growth by collisions is comparable to the growth via accretion. Our results are validated through the comparison with numerical simulations, and we overall conclude that super-competitive accretion is a valid mechanism for the formation of very massive objects in the early Universe. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: 9 pages, 5 figures. Accepted for publication in MNRAS

Signatures of winds and jets in the environment of supermassive black holes

Authors: Bidisha Bandyopadhyay, Christian Fendt, Dominik R. G. Schleicher, Javier Lagunas, Javier Pedreros, Neil M. Nagar, Felipe Agurto

Abstract: The Event Horizon Telescope Collaboration (EHTC) has presented first - dynamic-range limited - images of the black hole shadows in M87 and Sgr A*. The next generation Event Horizon Telescope (ngEHT) will provide higher sensitivity and higher dynamic range images (and movies) of these two sources plus image at least a dozen others at $\leq$100 gravitational radii resolution. We here perform an expl… ▽ More The Event Horizon Telescope Collaboration (EHTC) has presented first - dynamic-range limited - images of the black hole shadows in M87 and Sgr A*. The next generation Event Horizon Telescope (ngEHT) will provide higher sensitivity and higher dynamic range images (and movies) of these two sources plus image at least a dozen others at $\leq$100 gravitational radii resolution. We here perform an exploratory study of the appearance of winds and jets in such future observations. To do this we use M87 and Sgr A* as reference systems: we do not aim to exactly reproduce them, but rather to determine how their observed images will depend on specific physical assumptions. Even in the case of similar or the same dynamics, the images depend significantly on global parameters such as the black hole mass and the mass accretion rate. Our results provide guidance in the interpretation of future high-resolution images, particularly if a wind or jet is detected. △ Less

Submitted 12 March, 2023; originally announced March 2023.

Comments: 12 pages, 5 figures, Submitted to MNRAS

V4142 Sgr: a Double Periodic Variable with an accretor surrounded by the accretion-disk's atmosphere

Authors: J. A. Rosales, R. E. Mennickent, G. Djurašević, I. Araya, M. Curé, D. R. G. Schleicher, J. Petrović

Abstract: 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 pr… ▽ More 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$α$ double emission with a persistent $V \leq R$ asymmetry which is considered evidence of a possible wind emergin from the hotspot region... △ Less

Submitted 21 February, 2023; originally announced February 2023.

Journal ref: Astronomy & Astrophysics,670,2023,A94

Non-ideal magnetohydrodynamics of self-gravitating filaments

Authors: Nicol Gutiérrez-Vera, Tommaso Grassi, Stefano Bovino, Alessandro Lupi, Daniele Galli, Dominik R. G. Schleicher

Abstract: Filaments have been studied in detail through observations and simulations. A range of numerical works have separately investigated how chemistry and diffusion effects, as well as magnetic fields and their structure impact the gas dynamics of the filament. However, non-ideal effects have hardly been explored thus far. We investigate how non-ideal magnetohydrodynamic (MHD) effects, combined with a… ▽ More Filaments have been studied in detail through observations and simulations. A range of numerical works have separately investigated how chemistry and diffusion effects, as well as magnetic fields and their structure impact the gas dynamics of the filament. However, non-ideal effects have hardly been explored thus far. We investigate how non-ideal magnetohydrodynamic (MHD) effects, combined with a simplified chemical model affect the evolution and accretion of a star-forming filament. We modeled an accreting self-gravitating turbulent filament using lemongrab, a one-dimensional (1D) non-ideal MHD code that includes chemistry. We explore the influence of non-ideal MHD, the orientation and strength of the magnetic field, and the cosmic ray ionization rate, on the evolution of the filament, with particular focus on the width and accretion rate. We find that the filament width and the accretion rate are determined by the magnetic field properties, including the initial strength, the coupling with the gas controlled by the cosmic ray ionization rate, and the orientation of the magnetic field with respect to the accretion flow direction. Increasing the cosmic-ray ionization rate leads to a behavior closer to that of ideal MHD, reducing the magnetic pressure support and, hence, damping the accretion efficiency with a consequent broadening of the filament width. For the same reason, we obtained a narrower width and a larger accretion rate when we reduced the initial magnetic field strength. Overall, while these factors affect the final results by approximately a factor of~2, removing the non-ideal MHD effects results in a much greater variation (up to a factor of~7). The inclusion of non-ideal MHD effects and the cosmic-ray ionization is crucial for the study of self-gravitating filaments and in determining critical observable quantities, such as the filament width and accretion rate. △ Less

Submitted 28 November, 2022; originally announced November 2022.

Comments: 11 pages, 9 figures, to be published in Astronomy & Astrophysics

Journal ref: A&A 670, A38 (2023)

Formation of metal-free binaries: Impact of H$_{2}$ line cooling and CIE cooling

Authors: Rafeel Riaz, Dominik R. G. Schleicher, Stefano Bovino, Siegfried Vanaverbeke, Ralf S. Klessen

Abstract: During primordial star formation, the main cooling channel is provided by H$_{2}$ and super-molecules, such as H$_{2}$ or H$_{2}$, at sufficiently high densities. When the latter form at $n_{\rm H}$ $\geq$ $10^{14}$~cm$^{-3}$, collision-induced emission (CIE) provides efficient gas cooling. We investigate how CIE cooling affects the formation of metal-free binaries comparing simulations with and w… ▽ More During primordial star formation, the main cooling channel is provided by H$_{2}$ and super-molecules, such as H$_{2}$ or H$_{2}$, at sufficiently high densities. When the latter form at $n_{\rm H}$ $\geq$ $10^{14}$~cm$^{-3}$, collision-induced emission (CIE) provides efficient gas cooling. We investigate how CIE cooling affects the formation of metal-free binaries comparing simulations with and without this process. Irrespective of the cooling mechanism, we find a typical protostellar mass range between 0.01 to 100 M$_{\odot}$. However, models with only H$_{2}$ line cooling produce a greater number of low-mass protostars which exhibit stronger variations in their radial velocities than the high-mass protostars. Similarly, in models with both H$_{2}$ cooling and CIE cooling, significant variations in the radial velocities are found for protostars in the intermediate mass range. The initial number of fragments $N_{\rm max}$ decreases with increasing strength of turbulence. Cooling via super-molecules lets the most massive protobinaries (MMPBs) efficiently accrete mass. The maximum mass accretion rate $\dot M_{\rm max}$ for the MMPBs is more than an order of magnitude higher in the presence of CIE cooling than for pure H$_{2}$ line cooling. As a result, compact binaries with a semi-major axis as small as 3.57 au may form through the H$_{2}$ $-$ H$_{2}$ cooling channel. Our results indicate that in addition to the MMPBs most population III (Pop. III) binaries should be in eccentric i.e. non-circular orbits. This provides an important connection to the eccentric binaries reported in previous studies, which were found to exhibit rich temporal accretion signals during their evolution. △ Less

Submitted 8 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS. Contains 15 pages, 11 figures, 4 tables

Journal ref: Monthly Notices of the Royal Astronomical Society 2022

A model for the infrared-radio correlation of main-sequence galaxies at GHz frequencies and its variation on redshift and stellar mass

Authors: J. Schober, M. T. Sargent, R. S. Klessen, D. R. G. Schleicher

Abstract: 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… ▽ More 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. △ Less

Submitted 19 October, 2023; v1 submitted 14 October, 2022; originally announced October 2022.

Comments: 20 pages, 18 figures

Journal ref: A&A 679, A47 (2023)

Role of magnetic fields in the formation of direct collapse black holes

Authors: Muhammad A. Latif, Dominik R. G. Schleicher, Sadegh Khochfar

Abstract: Direct collapse black holes (DCBHs) are the leading candidates for the origin of the first supermassive black holes. However, the role of magnetic fields during their formation is still unclear as none of the previous studies has been evolved long enough to assess their impact during the accretion phase. Here, we report the results from a suite of 3D cosmological magneto-hydrodynamic (MHD) simulat… ▽ More Direct collapse black holes (DCBHs) are the leading candidates for the origin of the first supermassive black holes. However, the role of magnetic fields during their formation is still unclear as none of the previous studies has been evolved long enough to assess their impact during the accretion phase. Here, we report the results from a suite of 3D cosmological magneto-hydrodynamic (MHD) simulations which are evolved for 1.6 Myrs comparable to the expected lifetime of supermassive stars (SMSs). Our findings suggest that magnetic fields are rapidly amplified by strong accretion shocks irrespective of the initial magnetic field strength and reach the saturation state. They stabilize the accretion disks and significantly reduce fragmentation by enhancing the Jeans mass in comparison with pure hydrodynamical runs. Although the initial clump masses are larger in MHD runs, the rapid coalescence of clumps in non-MHD cases due to the higher degree of fragmentation results in similar masses. Overall, the central clumps have masses of $\rm 10^5~M_{\odot}$ and the mean mass accretion rates of $\rm \sim 0.1 ~M_{\odot}/yr$ are similar in both MHD and non-MHD cases. The multiplicity of SMSs is significantly reduced in MHD simulations. Such strongly amplified magnetic fields are expected to launch Jets and outflows which may be detected with upcoming radio telescopes. △ Less

Submitted 16 February, 2023; v1 submitted 7 October, 2022; originally announced October 2022.

Comments: To be published ApJ, comments are still welcome

Global instability by runaway collisions in nuclear stellar clusters: Numerical tests of a route for massive black hole formation

Authors: M. C. Vergara, A. Escala, D. R. G. Schleicher, B. Reinoso

Abstract: 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 colli… ▽ More 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. △ Less

Submitted 22 May, 2023; v1 submitted 29 September, 2022; originally announced September 2022.

Comments: 15 pages, 9 figures

Simulations of fully convective M dwarfs: dynamo action with varying magnetic Prandtl numbers

Authors: Carolina A. Ortiz-Rodríguez, Dominik R. G. Schleicher, Petri J. Käpylä, Felipe H. Navarrete

Abstract: M dwarfs are low-mass main-sequence stars, the most numerous type of stars in the solar neighbourhood, which are known to have significant magnetic activity. The aim of this work is to explore the dynamo solutions and magnetic fields of fully convective M dwarfs with varying magnetic Prandtl numbers ${\rm Pr_M}$, and a rotation period (${\rm P_{rot}}$) of 43 days. ${\rm Pr_M}$ is known to play an… ▽ More M dwarfs are low-mass main-sequence stars, the most numerous type of stars in the solar neighbourhood, which are known to have significant magnetic activity. The aim of this work is to explore the dynamo solutions and magnetic fields of fully convective M dwarfs with varying magnetic Prandtl numbers ${\rm Pr_M}$, and a rotation period (${\rm P_{rot}}$) of 43 days. ${\rm Pr_M}$ is known to play an important role in the dynamo action; dynamos for low-$\rm Pr_M$ and large-$\rm Pr_M$ have very different properties. We performed three-dimensional magnetohydrodynamical (MHD) numerical simulations with the ``star-in-a-box'' model using stellar parameters for an M5 dwarf with $0.21 M_{\odot}$. We found that the dynamo solutions are sensitive to ${\rm Pr_M}$. The simulations at this rotation period present periodic cycles of the large-scale magnetic field up to ${\rm Pr_M} \leq 2$; for higher values the cycles disappear and irregular solutions to arise. Our results are consistent with previous studies and suggest that the dynamos operating in fully convective stars behave similarly as those in partially convective stars. △ Less

Submitted 28 June, 2022; originally announced June 2022.

Comments: 3 pages, 1 figure, Accepted for publication in Boletín de la Asociación Argentina de Astronomía

Synthetic observations using POLARIS: an application to simulations of massive prestellar cores

Authors: Joaquin Zamponi, Andrea Giannetti, Stefano Bovino, Giovanni Sabatini, Dominik R. G. Schleicher, Bastian Körtgen, Stefan Reissl, Sebastian Wolf

Abstract: Young massive stars are usually found embedded in dense and massive molecular clumps which are known for being highly obscured and distant. During their formation process, the degree of deuteration can be used as a potential indicator of the very early formation stages. This is particularly effective when employing the abundance of H$_2$D$^+$. However, its low abundances and large distances make d… ▽ More Young massive stars are usually found embedded in dense and massive molecular clumps which are known for being highly obscured and distant. During their formation process, the degree of deuteration can be used as a potential indicator of the very early formation stages. This is particularly effective when employing the abundance of H$_2$D$^+$. However, its low abundances and large distances make detections in massive sources hard to achieve. We present an application of the radiative transfer code POLARIS, with the goal to test the observability of the ortho-H$_2$D$^+$ transition $1_{10}$-$1_{11}$ (372.42 GHz) using simulations of high-mass collapsing cores that include deuteration chemistry. We analyzed an early and a late stage of the collapse of a 60 M$_{\odot}$ core, testing different source distances. For all cases, we generated synthetic single-dish and interferometric observations and studied the differences in both techniques. The column densities we derive are comparable to values reported for similar sources. These estimates depend on the extent over which they are averaged, and sources with compact emission they can be highly affected by beam dilution. Combined ALMA-ACA observations improve in signal-to-noise ratio and lead to better column density estimates as compared to ALMA alone. We confirm the feasibility to study ortho-H$_2$D$^+$ emission up to distances of 7 kpc. We provide a proof-of-concept of our framework for synthetic observations and highlight its importance when comparing numerical simulations with real observations. This work also proves how relevant it is to combine single-dish and interferometric measurements to derive appropriate source column densities. △ Less

Submitted 15 June, 2022; originally announced June 2022.

Comments: Accepted for publication in Astrophysics & Space Science. arXiv admin note: substantial text overlap with arXiv:2009.00407

The Fermi paradox: Impact of astrophysical processes and dynamical evolution

Authors: Dominik R. G. Schleicher, Stefano Bovino

Abstract: The Fermi paradox has given rise to various attempts to explain why no evidence of extraterrestrial civilisations was found so far on Earth and in our Solar System. Here, we present a dynamical model for the development of such civilisations, which accounts for self-destruction, colonisation and astrophysical destruction mechanisms of civilisations including gamma-ray bursts, type Ia and type II s… ▽ More The Fermi paradox has given rise to various attempts to explain why no evidence of extraterrestrial civilisations was found so far on Earth and in our Solar System. Here, we present a dynamical model for the development of such civilisations, which accounts for self-destruction, colonisation and astrophysical destruction mechanisms of civilisations including gamma-ray bursts, type Ia and type II supernovae as well as radiation from the supermassive black hole. We adopt conservative estimates regarding the efficiency of such processes and find that astrophysical effects can influence the development of intelligent civilisations and change the number of systems with such civilisations by roughly a factor of 2; potentially more if the feedback is enhanced. Our results show that non-equilibrium evolution allows for solutions in-between extreme cases such as "rare Earth" or extreme colonisation, including scenarios with civilisation fractions between 10^{-2} and 10^{-7}. These would imply still potentially large distances to the next such civilisations, particularly when persistence phenomena are being considered. As previous studies, we confirm that the main uncertainties are due to the lifetime of civilisations as well as the assumed rate of colonisation. For SETI-like studies, we believe that unbiased searches are needed considering both the possibilities that the next civilisations are nearby or potentially very far away. △ Less

Submitted 14 June, 2022; originally announced June 2022.

Comments: 8 pages, 6 figures, 1 table. Accepted for publication in the International Journal of Astrobiology

Origin of eclipsing time variations in Post-Common-Envelope binaries: role of the centrifugal force

Authors: Felipe H. Navarrete, Dominik R. G. Schleicher, Petri J. Käpylä, Carolina Ortiz, Robi Banerjee

Abstract: Eclipsing time variations (ETVs) in post-common-envelope binaries (PCEBs) were proposed to be due to the time-varying component of the stellar gravitational quadrupole moment. This is suggested to be produced by changes in the stellar structure due to an internal redistribution of angular momentum and the effect of the centrifugal force. We examined this hypothesis and present 3D simulations of co… ▽ More Eclipsing time variations (ETVs) in post-common-envelope binaries (PCEBs) were proposed to be due to the time-varying component of the stellar gravitational quadrupole moment. This is suggested to be produced by changes in the stellar structure due to an internal redistribution of angular momentum and the effect of the centrifugal force. We examined this hypothesis and present 3D simulations of compressible magnetohydrodynamics (MHD) performed with the {\sc Pencil Code}. We modeled the stellar dynamo for a solar-mass star with angular velocities of 20 and 30 times solar. We included and varied the strength of the centrifugal force and compared the results with reference simulations without the centrifugal force and with a simulation in which its effect is enhanced. The centrifugal force causes perturbations in the evolution of the numerical model, so that the outcome in the details becomes different as a result of nonlinear evolution. While the average density profile is unaffected by the centrifugal force, a relative change in the density difference between high latitudes and the equator of $\sim10^{-4}$ is found. The power spectrum of the convective velocity is found to be more sensitive to the angular velocity than to the strength of the centrifugal force. The quadrupole moment of the stars includes a fluctuating and a time-independent component, which vary with the rotation rate. As very similar behavior is produced in absence of the centrifugal force, we conclude that it is not the main ingredient for producing the time-averaged and fluctuating quadrupole moment of the star. In a real physical system, we thus expect contributions from both components, that is, from the time-dependent gravitational force from the variation in the quadrupole term and from the spin-orbit coupling that is due to the persistent part of the quadrupole. △ Less

Submitted 30 September, 2022; v1 submitted 6 May, 2022; originally announced May 2022.

Comments: 8 pages, 7 figures. Accepted in A&A

Report number: NORDITA 2022-028

Journal ref: A&A 667, A164 (2022)

Origin of supermassive black holes in massive metal-poor protoclusters

Authors: D. R. G. Schleicher, B. Reinoso, M. Latif, R. S. Klessen, M. Z. C. Vergara, A. Das, P. Alister, V. B. Díaz, P. A. Solar

Abstract: 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_so… ▽ More 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. △ Less

Submitted 31 March, 2022; originally announced April 2022.

Comments: 9 pages, 1 figure, 1 table. Accepted for publication with MNRAS

Fundamental parameters of the eclipsing binary DD CMa and evidence for mass exchange

Authors: J. A. Rosales, R. E. Mennickent, G. Djurašević, J. F. González, I. Araya, M. Cabezas, D. R. G. Schleicher, M. Curé

Abstract: 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$β$ and H$α$ absorptions with weak emission shoulders and we also find color excess in the WISE multiband photometry, interpreted as signatures of circu… ▽ More 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$β$ and H$α$ 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$α$ and H$β$ 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. △ Less

Submitted 1 September, 2021; originally announced September 2021.

Comments: 14 pages with 13 figures

Journal ref: The Astronomical Journal, 161, 2012, 165

Stellar and accretion disk parameters of the close binary HD 50526

Authors: J. A. Rosales, R. E. Mennickent, G. Djurašević, D. R. G. Schleicher, S. Zharikov, I. Araya, L. Celedón, M. Curé

Abstract: 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 algori… ▽ More 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$α$ and H$β$ 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$α$ and H$β$, reveal complex structures of mass fluxes in the system. △ Less

Submitted 1 September, 2021; originally announced September 2021.

Comments: 16 pages with 14 figures

Journal ref: The Astronomical Journal, 162, 2021, 66

Turbulence and its connection to episodic accretion in binary YSOs

Authors: Rafeel Riaz, Dominik R. G. Schleicher, Siegfried Vanaverbeke, Ralf S. Klessen

Abstract: We report signatures of episodic accretion in young stellar objects (YSOs) that emerge in protobinary configurations in a gravoturbulent gas collapse. We find in most of these protobinary systems strong accretion bursts between the two companions with a recurrence time-scale of about 1 kyr. The accretion rate onto the secondary star typically exceeds that onto the primary with a peak value of 2… ▽ More We report signatures of episodic accretion in young stellar objects (YSOs) that emerge in protobinary configurations in a gravoturbulent gas collapse. We find in most of these protobinary systems strong accretion bursts between the two companions with a recurrence time-scale of about 1 kyr. The accretion rate onto the secondary star typically exceeds that onto the primary with a peak value of 2 $\times 10^{-2}$ M$_{\odot}$ yr$^{-1}$ for the former and 6 $\times 10^{-3}$ M$_{\odot}$ yr$^{-1}$ for the latter. We propose that the secondary companion which remains more active in its episodes of accretion bursts, especially for the gas cores with subsonic velocity dispersion, may provide observational opportunities to find traces of episodic accretion in the surrounding gas of the embedded YSOs that are in a binary configuration. Also, protostars evolving as single objects in the same environment show fewer accretion bursts and all together a more steady mass growth history. The prestellar cores with subsonic velocity dispersion exhibit an order of magnitude more intense accretion bursts than in the case of cores with supersonic velocity dispersions. The latter shows the formation of some of the protobinaries in which the primary acts as a more actively accreting companion. This can support these binaries to become systems of extreme mass ratio. Moreover, the YSOs in binary configurations with small semi-major axis $a$ $\approx$ 50 au and high mass ratio $q$ > 0.7 support phases of intense episodic accretion. The eccentricity, however, seems to play no significant role in the occurrence of accretion bursts. △ Less

Submitted 27 August, 2021; originally announced August 2021.

Comments: Accepted in MNRAS, 18 pages, 18 figures, 3 tables

Effect of mass loss due to stellar winds on the formation of supermassive black hole seeds in dense nuclear star clusters

Authors: Arpan Das, Dominik R. G. Schleicher, Shantanu Basu, Tjarda C. N. Boekholt

Abstract: The observations of high redshifts quasars at $z\gtrsim 6$ have revealed that supermassive black holes (SMBHs) of mass $\sim 10^9\,\mathrm{M_{\odot}}$ were already in place within the first $\sim$ Gyr after the Big Bang. Supermassive stars (SMSs) with masses $10^{3-5}\,\mathrm{M_{\odot}}$ are potential seeds for these observed SMBHs. A possible formation channel of these SMSs is the interplay of g… ▽ More The observations of high redshifts quasars at $z\gtrsim 6$ have revealed that supermassive black holes (SMBHs) of mass $\sim 10^9\,\mathrm{M_{\odot}}$ were already in place within the first $\sim$ Gyr after the Big Bang. Supermassive stars (SMSs) with masses $10^{3-5}\,\mathrm{M_{\odot}}$ are potential seeds for these observed SMBHs. A possible formation channel of these SMSs is the interplay of gas accretion and runaway stellar collisions inside dense nuclear star clusters (NSCs). However, mass loss due to stellar winds could be an important limitation for the formation of the SMSs and affect the final mass. In this paper, we study the effect of mass loss driven by stellar winds on the formation and evolution of SMSs in dense NSCs using idealised N-body simulations. Considering different accretion scenarios, we have studied the effect of the mass loss rates over a wide range of metallicities $Z_\ast=[.001-1]\mathrm{Z_{\odot}}$ and Eddington factors $f_{\rm Edd}=L_\ast/L_{\mathrm{Edd}}=0.5,0.7,\,\&\, 0.9$. For a high accretion rate of $10^{-4}\,\mathrm{M_{\odot}yr^{-1}}$, SMSs with masses $\gtrsim 10^3\MSun$ could be formed even in a high metallicity environment. For a lower accretion rate of $10^{-5}\,\mathrm{M_{\odot}yr^{-1}}$, SMSs of masses $\sim 10^{3-4}\,\mathrm{M_{\odot}}$ can be formed for all adopted values of $Z_\ast$ and $f_{\rm Edd}$, except for $Z_\ast=\mathrm{Z_{\odot}}$ and $f_{\rm Edd}=0.7$ or 0.9. For Eddington accretion, SMSs of masses $\sim 10^3\,\mathrm{M_{\odot}}$ can be formed in low metallicity environments with $Z_\ast\lesssim 0.01\mathrm{Z_{\odot}}$. The most massive SMSs of masses $\sim 10^5\,\mathrm{M_{\odot}}$ can be formed for Bondi-Hoyle accretion in environments with $Z_\ast \lesssim 0.5\mathrm{Z_{\odot}}$. △ Less

Submitted 7 May, 2021; originally announced May 2021.

Chemical analysis of prestellar cores in Ophiuchus yields short timescales and rapid collapse

Authors: Stefano Bovino, Alessandro Lupi, Andrea Giannetti, Giovanni Sabatini, Dominik R. G. Schleicher, Friedrich Wyrowski, Karl M. Menten

Abstract: Sun-like stars form from the contraction of cold and dense interstellar clouds. How the collapse proceeds and what are the main physical processes driving it, however, is still under debate and a final consensus on the timescale of the process has not been reached. Does this contraction proceed slowly, sustained by strong magnetic fields and ambipolar diffusion, or is it driven by fast collapse wi… ▽ More Sun-like stars form from the contraction of cold and dense interstellar clouds. How the collapse proceeds and what are the main physical processes driving it, however, is still under debate and a final consensus on the timescale of the process has not been reached. Does this contraction proceed slowly, sustained by strong magnetic fields and ambipolar diffusion, or is it driven by fast collapse with gravity dominating the entire process? One way to answer this question is to measure the age of prestellar cores through statistical methods based on observations or via reliable chemical chronometers, which should better reflect the physical conditions of the cores. Here we report APEX observations of ortho-H$_2$D$^+$ and para-D$_2$H$^+$ for six cores in the Ophiuchus complex and combine them with detailed three-dimensional magneto-hydrodynamical simulations including chemistry, providing a range of ages for the observed cores up to 200 kyr. The outcome of our simulations and subsequent analysis provides a good match with the observational results in terms of physical (core masses and volume densities) and dynamical parameters such as the Mach number and the virial parameter. We show that models of fast collapse successfully reproduce the observed range of chemical abundance ratios as the timescales to reach the observed stages is comparable to the dynamical time of the cores (i.e. the free-fall time) and much shorter than the ambipolar diffusion time, measured from the electron fraction in the simulations. To confirm that this ratio can be used to distinguish between different star-formation scenarios a larger (statistically relevant) sample of star-forming cores should be explored. △ Less

Submitted 17 June, 2021; v1 submitted 5 May, 2021; originally announced May 2021.

Comments: submitted to A&A, 12 pages, 5 figures, 4 tables

Identification of prestellar cores in high-mass star forming clumps via $\rm H_2D^+$ observations with ALMA

Authors: E. Redaelli, S. Bovino, A. Giannetti, G. Sabatini, P. Caselli, F. Wyrowski, D. R. G. Schleicher, D. Colombo

Abstract: 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… ▽ More 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. △ Less

Submitted 13 April, 2021; originally announced April 2021.

Comments: Accepted for publications in A&A

Journal ref: A&A 650, A202 (2021)

Stellar collisions in flattened and rotating Pop. III star clusters

Authors: M. Z. C. Vergara, D. R. G. Schleicher, T. C. N. Boekholt, B. Reinoso, M. Fellhauer, R. S. Klessen, N. W. C. Leigh

Abstract: 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 c… ▽ More 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. △ Less

Submitted 3 April, 2021; originally announced April 2021.

Comments: 10 pages, 7 figures

Origin of eclipsing time variations: Contributions of different modes of the dynamo-generated magnetic field

Authors: Felipe H. Navarrete, Petri J. Käpylä, Dominik R. G. Schleicher, Carolina A. Ortiz, Robi Banerjee

Abstract: The possibility to detect circumbinary planets and to study stellar magnetic fields through eclipsing time variations (ETVs) in binary stars has sparked an increase of interest in this area of research. We revisit the connection between stellar magnetic fields and the gravitational quadrupole moment $Q_{xx}$ and compare different dynamo-generated ETV models with our simulations. We present magneto… ▽ More The possibility to detect circumbinary planets and to study stellar magnetic fields through eclipsing time variations (ETVs) in binary stars has sparked an increase of interest in this area of research. We revisit the connection between stellar magnetic fields and the gravitational quadrupole moment $Q_{xx}$ and compare different dynamo-generated ETV models with our simulations. We present magnetohydrodynamical simulations of solar mass stars with rotation periods of 8.3, 1.2, and 0.8 days and perform a detailed analysis of the magnetic and quadrupole moment using spherical harmonic decomposition. The extrema of $Q_{xx}$ are associated with changes in the magnetic field structure. This is evident in the simulation with a rotation period of 1.2 days. Its magnetic field has a more complex behavior than in the other models, as the large-scale nonaxisymmetric field dominates throughout the simulation and the axisymmetric component is predominantly hemispheric. This triggers variations in the density field that follow the magnetic field asymmetry with respect to the equator, affecting the $zz$ component of the inertia tensor, and thus modulating $Q_{xx}$. The magnetic field of the two other runs are less variable in time and more symmetric with respect to the equator, such that the variations in the density are weaker, and therefore only small variations in $Q_{xx}$ are seen. If interpreted via the classical Applegate mechanism (tidal locking), the quadrupole moment variations obtained in the current simulations are about two orders of magnitude below those deduced from observations of post-common-envelope binaries. However, if no tidal locking is assumed, our results are compatible with the observed ETVs. △ Less

Submitted 27 May, 2022; v1 submitted 22 February, 2021; originally announced February 2021.

Comments: 15 pages, 16 figures. Accepted for publication in A&A

Report number: NORDITA 2021-017

Journal ref: A&A 663, A90 (2022)

General relativistic radiation transport: Implications for VLBI/EHT observations of AGN discs, winds and jets

Authors: Bidisha Bandyopadhyay, Christian Fendt, Dominik R. G. Schleicher, Christos Vourellis

Abstract: In 2019, the Event Horizon Telescope Collaboration (EHTC) has published the first image of a supermassive black hole (SMBH) obtained via the Very Large Baseline Interferometry (VLBI) technique. In the future, it is expected that additional and more sensitive VLBI observations will be pursued for other nearby Active Galactic Nuclei (AGN), and it is therefore important to understand which possible f… ▽ More In 2019, the Event Horizon Telescope Collaboration (EHTC) has published the first image of a supermassive black hole (SMBH) obtained via the Very Large Baseline Interferometry (VLBI) technique. In the future, it is expected that additional and more sensitive VLBI observations will be pursued for other nearby Active Galactic Nuclei (AGN), and it is therefore important to understand which possible features can be expected in such images. In this paper, we post-process General Relativistic Magneto-Hydrodynamical (GR-MHD) simulations which include resistivity, thus providing a self-consistent jet formation model, including resistive mass loading of a wind launched from a disc in Keplerian rotation. The ray-tracing is done using the General Relativistic Ray-Tracing code GRTRANS assuming synchrotron emission. We study the appearance of the black hole environment including the accretion disc, winds and jets under a large range of condition, varying black hole mass, accretion rate, spin, inclination angle, disc parameters and observed frequency. When we adopt M87-like parameters, we show that we can reproduce a ring-like feature (similar as observed by the EHT) for some of our simulations. The latter suggests that such Keplerian disc models thus could be consistent with the observed results. Depending on their masses, accretion rates, spin and the sensitivity of the observation, we note that other SMBHs may show additional features like winds and jets in the observations. △ Less

Submitted 6 September, 2021; v1 submitted 21 February, 2021; originally announced February 2021.

Comments: 22 pages, 19 figures, Version matches with the accepted version for publication to MNRAS

Formation of supermassive black hole seeds in nuclear star clusters via gas accretion and runaway collisions

Authors: Arpan Das, Dominik R. G. Schleicher, Nathan W. C. Leigh, Tjarda C. N. Boekholt

Abstract: More than two hundred supermassive black holes (SMBHs) of masses $\gtrsim 10^9\,\mathrm{M_{\odot}}$ have been discovered at $z \gtrsim 6$. One promising pathway for the formation of SMBHs is through the collapse of supermassive stars (SMSs) with masses $\sim 10^{3-5}\,\mathrm{M_{\odot}}$ into seed black holes which could grow upto few times $10^9\,\mathrm{M_{\odot}}$ SMBHs observed at $z\sim 7$. I… ▽ More More than two hundred supermassive black holes (SMBHs) of masses $\gtrsim 10^9\,\mathrm{M_{\odot}}$ have been discovered at $z \gtrsim 6$. One promising pathway for the formation of SMBHs is through the collapse of supermassive stars (SMSs) with masses $\sim 10^{3-5}\,\mathrm{M_{\odot}}$ into seed black holes which could grow upto few times $10^9\,\mathrm{M_{\odot}}$ SMBHs observed at $z\sim 7$. In this paper, we explore how SMSs with masses $\sim 10^{3-5}\,\mathrm{M_{\odot}}$ could be formed via gas accretion and runaway stellar collisions in high-redshift, metal-poor nuclear star clusters (NSCs) using idealised N-body simulations. We explore physically motivated accretion scenarios, e.g. Bondi-Hoyle-Lyttleton accretion and Eddington accretion, as well as simplified scenarios such as constant accretions. While gas is present, the accretion timescale remains considerably shorter than the timescale for collisions with the most massive object (MMO). However, overall the timescale for collisions between any two stars in the cluster can become comparable or shorter than the accretion timescale, hence collisions still play a crucial role in determining the final mass of the SMSs. We find that the problem is highly sensitive to the initial conditions and our assumed recipe for the accretion, due to the highly chaotic nature of the problem. The key variables that determine the mass growth mechanism are the mass of the MMO and the gas reservoir that is available for the accretion. Depending on different conditions, SMSs of masses $\sim10^{3-5} \,\mathrm{M_{\odot}}$ can form for all three accretion scenarios considered in this work. △ Less

Submitted 14 May, 2021; v1 submitted 2 December, 2020; originally announced December 2020.

Stellar Initial Mass Function over a range of redshifts

Authors: Rafeel Riaz, Dominik R. G. Schleicher, Siegfried Vanaverbeke, Ralf S. Klessen

Abstract: The stellar Initial Mass Function (IMF) seems to be close to universal in the local star-forming regions. However, this quantity of a newborn stellar population responds differently at gas metallicities $Z \sim$ $Z_\odot$ than $Z$ = 0. A view on the cosmic star formation history suggests that the cooling agents in the gas vary both in their types and molecular abundances. For instance, in the prim… ▽ More The stellar Initial Mass Function (IMF) seems to be close to universal in the local star-forming regions. However, this quantity of a newborn stellar population responds differently at gas metallicities $Z \sim$ $Z_\odot$ than $Z$ = 0. A view on the cosmic star formation history suggests that the cooling agents in the gas vary both in their types and molecular abundances. For instance, in the primordial gas environment, the gas temperature can be higher by a factor of 30 as compared to the present day. Stellar radiation feedback and cosmic microwave background (CMB) radiation may even contribute towards increasing the floor temperature of the star-forming gas which subsequently can leave profound impacts on the IMF. We present the contribution of the radiation sources towards the thermodynamical evolution of the Jeans unstable gas cloud which experiences fragmentation and mass accretion. We find evidence which suggests that the latter becomes the dominant process for star formation efficiencies (SFE) above 5 - 7 %, thus increasing the average mass of the stars. We focus on the isolated and binary stellar configurations emerging during the gas collapse. The binary fraction on average remains 0.476 and contributes significantly towards the total SFE of 15 %. △ Less

Submitted 13 October, 2020; originally announced October 2020.

Comments: Accepted for publication in Astronomische Nachrichten. 6 pages, 5 figures, and 1 table. arXiv admin note: text overlap with arXiv:2003.07639

Synthetic observations of deuterated molecules in massive prestellar cores

Authors: Joaquin Zamponi F., Andrea Giannetti, Stefano Bovino, Giovanni Sabatini, Dominik R. G. Schleicher, Bastian Körtgen, Stefan Reissl, Sebastian Wolf

Abstract: Young massive stars are usually found embedded in dense massive molecular clumps and are known for being highly obscured and distant. During their formation process, deuteration is regarded as a potentially good indicator of the very early formation stages. In this work, we test the observability of the ground-state transition of ortho-H$_2$D$^+$ $J_{\rm {K_a, K_c}} = 1_{10}$-$1_{11} $ by performi… ▽ More Young massive stars are usually found embedded in dense massive molecular clumps and are known for being highly obscured and distant. During their formation process, deuteration is regarded as a potentially good indicator of the very early formation stages. In this work, we test the observability of the ground-state transition of ortho-H$_2$D$^+$ $J_{\rm {K_a, K_c}} = 1_{10}$-$1_{11} $ by performing interferometric and single-dish synthetic observations using magneto-hydrodynamic simulations of high-mass collapsing molecular cores, including deuteration chemistry. We studied different evolutionary times and source distances (from 1 to 7 kpc) to estimate the information loss when comparing the column densities inferred from the synthetic observations to the column densities in the model. We mimicked single-dish observations considering an APEX-like beam and interferometric observations using CASA and assuming the most compact configuration for the ALMA antennas. We found that, for centrally concentrated density distributions, the column densities are underestimated by about 51% in the case of high-resolution ALMA observations ($\leqslant$1") and up to 90% for APEX observations (17"). Interferometers retrieve values closer to the real ones, however, their finite spatial sampling results in the loss of contribution from large-scale structures due to the lack of short baselines. We conclude that, the emission of o-H$_2$D$^+$ in distant massive dense cores is faint and would require from $\sim$1 to $\sim$7 hours of observation at distances of 1 and 7 kpc, respectively, to achieve a 14$σ$ detection in the best case scenario. Additionally, the column densities derived from such observations will certainly be affected by beam dilution in the case of single-dishes and spatial filtering in the case of interferometers. △ Less

Submitted 1 September, 2020; originally announced September 2020.

Comments: 12 pages, 8 figures, submitted to MNRAS

Magnetic field amplification in accretion discs around the first stars: implications for the primordial IMF

Authors: Piyush Sharda, Christoph Federrath, Mark R. Krumholz, Dominik R. G. Schleicher

Abstract: Magnetic fields play an important role in the dynamics of present-day molecular clouds. Recent work has shown that magnetic fields are equally important for primordial clouds, which form the first stars in the Universe. While the primordial magnetic field strength on cosmic scales is largely unconstrained, theoretical models strongly suggest that a weak seed field existed in the early Universe. We… ▽ More Magnetic fields play an important role in the dynamics of present-day molecular clouds. Recent work has shown that magnetic fields are equally important for primordial clouds, which form the first stars in the Universe. While the primordial magnetic field strength on cosmic scales is largely unconstrained, theoretical models strongly suggest that a weak seed field existed in the early Universe. We study how the amplification of such a weak field can influence the evolution of accretion discs around first stars, and thus affect the primordial initial mass function (IMF). We perform a suite of 3D ideal magneto-hydrodynamic (MHD) simulations with different initial field strengths and numerical resolutions. We find that, in simulations with sufficient spatial resolution to resolve the Jeans scale during the collapse, even initially weak magnetic fields grow exponentially to become dynamically important due to both the so-called 'small-scale turbulent dynamo' and the 'large-scale mean-field dynamo'. Capturing the small-scale dynamo action depends primarily on how well we resolve the Jeans length, while capturing the large-scale dynamo depends on the Jeans resolution as well as the maximum absolute resolution. Provided enough resolution, we find that fragmentation does not depend strongly on the initial field strength, because even weak fields grow to become strong. However, fragmentation in runs with magnetic fields differs significantly from those without magnetic fields. We conclude that the development of dynamically strong magnetic fields during the formation of the first stars is likely inevitable, and that these fields had a significant impact on the primordial IMF. △ Less

Submitted 22 February, 2021; v1 submitted 6 July, 2020; originally announced July 2020.

Comments: 20 pages, 10 figures, 1 appendix. MNRAS, in press

The effects of a background potential in star cluster evolution: a delay in the relaxation time-scale and runaway collision processes

Authors: B. Reinoso, D. R. G. Schleicher, M. Fellhauer, N. W. C. Leigh, R. S. Klessen

Abstract: 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 explic… ▽ More 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. △ Less

Submitted 15 May, 2020; originally announced May 2020.

Comments: 16 pages. Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 639, A92 (2020)

Do fragmentation and accretion affect the stellar Initial Mass Function?

Authors: R. Riaz, D. R. G. Schleicher, S. Vanaverbeke, Ralf S. Klessen

Abstract: 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… ▽ More 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. △ Less

Submitted 17 March, 2020; originally announced March 2020.

Comments: 12 Pages, 8 Figures, 1 Table; Accepted for publication in MNRAS

A new proxy to estimate the cosmic-ray ionisation rate in dense cores

Authors: Stefano Bovino, Simon Ferrada-Chamorro, Alessandro Lupi, Dominik R. G. Schleicher, Paola Caselli

Abstract: Cosmic rays are a global source of ionisation, and the ionisation fraction represents a fundamental parameter in the interstellar medium. Ions couple to magnetic fields, affect the chemistry, and the dynamics of star-forming regions as well as planetary atmospheres. However, the cosmic-ray ionisation rate represents one of the bottlenecks for astrochemical models, and its determination is one of t… ▽ More Cosmic rays are a global source of ionisation, and the ionisation fraction represents a fundamental parameter in the interstellar medium. Ions couple to magnetic fields, affect the chemistry, and the dynamics of star-forming regions as well as planetary atmospheres. However, the cosmic-ray ionisation rate represents one of the bottlenecks for astrochemical models, and its determination is one of the most puzzling problems in astrophysics. While for diffuse clouds reasonable values have been provided from H$_3^+$ observations, for dense clouds, due to the lack of rotational transitions, this is not possible, and estimates are strongly biased by the employed model. We present here an analytical expression, obtained from first principles, to estimate the cosmic-ray ionisation rate from observational quantities. The theoretical predictions are validated with high-resolution three-dimensional numerical simulations and applied to the well known core L1544; we obtained an estimate of $ζ_2 \sim 2-3 \times 10^{-17}$ s$^{-1}$. Our results and the analytical formulae provided represent the first model-independent, robust tool to probe the cosmic-ray ionisation rate in the densest part of star-forming regions (on spatial scales of $R \leq 0.05$ pc). An error analysis is presented to give statistical relevance to our study. △ Less

Submitted 11 March, 2020; originally announced March 2020.

Comments: Accepted on MNRAS Letters

Formation of SMBH seeds in Pop III star clusters through collisions : the importance of mass loss

Authors: P. J. Alister Seguel, D. R. G. Schleicher, T. C. N. Boekholt, M. Fellhauer, R. S. Klessen

Abstract: 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 cons… ▽ More 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. △ Less

Submitted 13 February, 2020; v1 submitted 3 December, 2019; originally announced December 2019.

Comments: 12 pages, 9 figures, accepted by MNRAS

The 3D structure of CO depletion in high-mass prestellar regions

Authors: S. Bovino, S. Ferrada-Chamorro, A. Lupi, G. Sabatini, A. Giannetti, D. R. G. Schleicher

Abstract: 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 firs… ▽ More 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. △ Less

Submitted 30 October, 2019; originally announced October 2019.

Comments: accepted on ApJ

Resolving accretion flows in nearby active galactic nuclei with the Event Horizon Telescope

Authors: Bidisha Bandyopadhyay, Fu-Guo Xie, Neil M. Nagar, Dominik R. G. Schleicher, Venkatessh Ramakrishnan, Patricia Arévalo, Elena López, Yahelyn Diaz

Abstract: The Event Horizon Telescope (EHT), now with its first ever image of the photon ring around the supermassive black hole of M87, provides a unique opportunity to probe the physics of supermassive black holes through Very Long Baseline Interferometry (VLBI), such as the existence of the event horizon, the accretion processes as well as jet formation in Low Luminosity AGN (LLAGN). We build a theoretic… ▽ More The Event Horizon Telescope (EHT), now with its first ever image of the photon ring around the supermassive black hole of M87, provides a unique opportunity to probe the physics of supermassive black holes through Very Long Baseline Interferometry (VLBI), such as the existence of the event horizon, the accretion processes as well as jet formation in Low Luminosity AGN (LLAGN). We build a theoretical model which includes an Advection Dominated Accretion Flow (ADAF) with emission from thermal and non-thermal electrons in the flow and a simple radio jet outflow. The predicted spectral energy distribution (SED) of this model is compared to sub-arcsec resolution observations to get the best estimates of the model parameters. The model-predicted radial emission profiles at different frequency bands are used to predict whether the inflow can be resolved by the EHT or with telescopes such as the Global 3-mm VLBI array (GMVA). In this work the model is initially tested with high resolution SED data of M87 and then applied to our sample of 5 galaxies (Cen A, M84, NGC 4594, NGC 3998 and NGC 4278). The model then allows us to predict if one can detect and resolve the inflow for any of these galaxies using the EHT or GMVA within an 8 hour integration time. △ Less

Submitted 19 November, 2019; v1 submitted 12 July, 2019; originally announced July 2019.

Comments: 17 pages, 8 figures. Published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 490, Issue 4, p.4606-4621, 2019

Magneto-hydrodynamical origin of eclipsing time variations in post-common-envelope binaries for solar mass secondaries

Authors: Felipe H. Navarrete, Dominik R. G. Schleicher, Petri J. Käpylä, Jennifer Schober, Marcel Völschow, Ronald E. Mennickent

Abstract: Eclipsing time variations have been observed for a wide range of binary systems, including post-common-envelope binaries. A frequently proposed explanation, apart from the possibility of having a third body, is the effect of magnetic activity, which may alter the internal structure of the secondary star, particularly its quadrupole moment, and thereby cause quasi-periodic oscillations. Here we pre… ▽ More Eclipsing time variations have been observed for a wide range of binary systems, including post-common-envelope binaries. A frequently proposed explanation, apart from the possibility of having a third body, is the effect of magnetic activity, which may alter the internal structure of the secondary star, particularly its quadrupole moment, and thereby cause quasi-periodic oscillations. Here we present two compressible non-ideal magneto-hydrodynamical (MHD) simulations of the magnetic dynamo in a solar mass star, one of them with three times the solar rotation rate ("slow rotator"), the other one with twenty times the solar rotation rate ("rapid rotator"), to account for the high rotational velocities in close binary systems. For the slow rotator, we find that both the magnetic field and the stellar quadrupole moment change in a quasi-periodic manner, leading to O-C (observed - corrected times of the eclipse) variations of ~0.025 s. For the rapid rotator, the behavior of the magnetic field as well as the quadrupole moment changes become considerably more complex, due to the less coherent dynamo solution. The resulting O-C variations are of the order 0.13 s. The observed system V471~Tau shows two modes of eclipsing time variations, with amplitudes of 151 s and 20 s, respectively. However, the current simulations may not capture all relevant effects due to the neglect of the centrifugal force and self-gravity. Considering the model limitations and that the rotation of V471 Tau is still a factor of 2.5 faster than our rapid rotator, it may be conceivable to reach the observed magnitudes. △ Less

Submitted 29 October, 2019; v1 submitted 16 June, 2019; originally announced June 2019.

Comments: 16 pages, 24 figures, accepted for publication with MNRAS

Intermittent fragmentation and statistical variations during gas collapse in magnetised atomic cooling haloes

Authors: Philipp Grete, Muhammad A. Latif, Dominik R. G. Schleicher, Wolfram Schmidt

Abstract: Observations reveal the presence of supermassive black holes (SMBH) as early as ~700 million years after the Big Bang. Their formation path is still subject to current debate. We explore the influence of magnetic fields, which are strongly amplified via the turbulent small-scale dynamo, on the formation of SMBH seeds within the direct collapse scenario. In this study, we perform for the first time… ▽ More Observations reveal the presence of supermassive black holes (SMBH) as early as ~700 million years after the Big Bang. Their formation path is still subject to current debate. We explore the influence of magnetic fields, which are strongly amplified via the turbulent small-scale dynamo, on the formation of SMBH seeds within the direct collapse scenario. In this study, we perform for the first time cosmological magnetohydrodynamic large eddy simulations that employ a model for unresolved, compressible MHD turbulence. In total we perform 36 simulations for 9 haloes each with two different initial magnetic field strengths, and with and without employing the unresolved turbulence model. We make use of the adaptive mesh refinement approach to achieve an effective spatial resolution of less than one proper astronomical unit. We consider a regime where cooling is regulated by atomic hydrogen and the molecular hydrogen gets dissociated by a strong radiation field. Our main finding is that the majority of the gas properties in the haloes at the final output are predominantly determined by the run-away gravitational collapse. Turbulence is supersonic and super-Alfvénic in all cases, and magnetic fields are amplified to an approximately dynamically relevant regime. Finally, fragmentation during the collapse is intermittent and mass accretion rates range from 0.2-3 Msun/yr. This suggests that the presence of strongly amplified magnetic fields and turbulence provides additional pressure support on small scales and make the direct collapse a viable scenario for the formation of massive objects under the required ambient conditions. △ Less

Submitted 4 June, 2019; v1 submitted 28 February, 2019; originally announced March 2019.

Comments: 11 pages, 7 figures, accepted for publication in MNRAS

Discovery of a nitrogen-enhanced mildly metal-poor binary system: Possible evidence for pollution from an extinct AGB Star

Authors: José G. Fernández-Trincado, Ronald Mennickent, Mauricio Cabezas, Olga Zamora, Sarah L. Martell, Timothy C. Beers, Vinicius M. Placco, David M. Nataf, Szabolcs Mészáros, Dante Minniti, Dominik R. G. Schleicher, Baitian Tang, Angeles Pérez-Villegas, Annie C. Robin, Céline Reylé, Mario Ortigoza-Urdaneta

Abstract: We report the serendipitous discovery of a nitrogen-rich, mildly metal-poor ([Fe/H]=-1.08) giant star in a single-lined spectroscopic binary system found in the SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE-2) survey, Data Release 14 (DR14). Previous work has assumed that the two percent of halo giants with unusual elemental abundances have been evaporated from globular cl… ▽ More We report the serendipitous discovery of a nitrogen-rich, mildly metal-poor ([Fe/H]=-1.08) giant star in a single-lined spectroscopic binary system found in the SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE-2) survey, Data Release 14 (DR14). Previous work has assumed that the two percent of halo giants with unusual elemental abundances have been evaporated from globular clusters, but other origins for their abundance signatures, including binary mass transfer, must also be explored. We present the results of an abundance re-analysis of the APOGEE-2 high-resolution near-infrared spectrum of 2M12451043+1217401 with the Brussels Automatic Stellar Parameter (BACCHUS) automated spectral analysis code, and re-derive manually the main element families, namely the light elements (C, N), elements (O, Mg, Si), iron-peak element (Fe), \textit{s}-process element (Ce), and the light odd-Z element (Al). Our analysis confirm the N-rich nature of 2M12451043+1217401, which has a [N/Fe] ratio of $+0.69$, and shows that the abundances of C and Al are slightly discrepant from that of a typical mildly metal-poor RGB star, but exhibit Mg, Si, O and \textit{s}-process abundances (Ce) of typical field stars. We also detect a particularly large variability in its radial velocity over the period of the APOGEE-2 observations, and the most likely orbit fit to the radial velocity data has a period of $730.89\pm106.86$ days, a velocity semi-amplitude of $9.92 \pm 0.14$ km s$^{-1}$, and an eccentricity of $\sim 0.1276 \pm0.1174$, which support the hypothesis of a binary companion, and that has probably been polluted by a now-extinct AGB star. △ Less

Submitted 3 September, 2019; v1 submitted 27 February, 2019; originally announced February 2019.

Comments: Accepted for publication in A&A, 13 pages, 8 figures, 7 tables

Journal ref: A&A 631, A97 (2019)

A timeline for massive star-forming regions via combined observation of o-H$_2$D$^+$ and N$_2$D$^+$

Authors: A. Giannetti, S. Bovino, P. Caselli, S. Leurini, D. R. G. Schleicher, B. Koertgen, K. M. Menten, T. Pillai, F. Wyrowski

Abstract: 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,… ▽ More 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. △ Less

Submitted 7 January, 2019; originally announced January 2019.

Accretion models for LLAGNs: Model Parameter Estimation and Prediction of their Detectibility

Authors: Bidisha Bandyopadhyay, Dominik R. G. Schleicher, Neil Nagar, Fu-Guo Xie, Venkatessh Ramakrishnan

Abstract: The Event Horizon Telescope (EHT) provides a unique opportunity to probe the physics of supermassive black holes through Very Large Baseline Interferometry (VLBI), such as the existence of the event horizon, the accretion processes as well as jet formation. We build a theoretical model which includes an Advection Dominated Accretion Flow (ADAF) and a simple radio jet outflow. The predicted spectra… ▽ More The Event Horizon Telescope (EHT) provides a unique opportunity to probe the physics of supermassive black holes through Very Large Baseline Interferometry (VLBI), such as the existence of the event horizon, the accretion processes as well as jet formation. We build a theoretical model which includes an Advection Dominated Accretion Flow (ADAF) and a simple radio jet outflow. The predicted spectral energy distribution (SED) of this model is compared to observations to get the best estimates of the model parameters. Also the model-predicted radial emission profiles at different frequency bands can be used to predict whether the inflow can be resolved by the EHT. △ Less

Submitted 10 December, 2018; originally announced December 2018.

Comments: 3 pages, 3 figures. Submitted to BAAA (proceedings for the Binational AAA-SOCHIAS meeting 2018)

Black hole formation in the context of dissipative dark matter

Authors: M. A. Latif, A. Lupi, D. R. G. Schleicher, G. D'Amico, P. Panci, S. Bovino

Abstract: 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… ▽ More 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. △ Less

Submitted 7 December, 2018; originally announced December 2018.

Formation of massive black holes via collisions and accretion

Authors: D. R. G. Schleicher, M. A. Fellhauer, T. Boekholt, B. Reinoso, R. S. Klessen, M. Z. C. Vergara, P. J. Alister Seguel, S. Bovino, C. Olave, V. B. Díaz, P. Fibla, R. Riaz, B. Bandyopadhyay, R. I. San Martin-Perez, J. Zamponi, L. Haemmerle

Abstract: 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. Sim… ▽ More 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. △ Less

Submitted 5 December, 2018; originally announced December 2018.

Comments: 3 pages, 1 figure, 1 table. Submitted to BAAA (proceedings for the Binational AAA-SOCHIAS meeting 2018)