Source author record

Ralf Klessen

Ralf Klessen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2022arXiv

Primordial magnetic fields in Population III star formation: a magnetised resolution study

Population III stars form in groups due to the fragmentation of primordial gas. While uniform magnetic fields have been shown to support against fragmentation in present day star formation, it is unclear whether realistic k^3/2 primordial fields can have the same effect. We bypass the issues associated with simulating the turbulent dynamo by introducing a saturated magnetic field at equipartition with the velocity field when the central densities reaches 10-13 g cm-3. We test a range of sink particle creation densities from 10-10-10-8 g cm-3. Within the range tested, the fields did not suppress fragmentation of the gas and hence could not prevent the degree of fragmentation from increasing with increased resolution. The number of sink particles formed and total mass in sink particles was unaffected by the magnetic field across all seed fields and resolutions. The magnetic pressure remained sub-dominant to the gas pressure except in the highest density regions of the simulation box, where it became equal to but never exceeded gas pressure. Our results suggest that the inclusion of magnetic fields in numerical simulations of Pop III star formation is largely unimportant.

preprint2020arXiv

The molecular cloud lifecycle

Giant molecular clouds (GMCs) and their stellar offspring are the building blocks of galaxies. The physical characteristics of GMCs and their evolution are tightly connected to galaxy evolution. The macroscopic properties of the interstellar medium propagate into the properties of GMCs condensing out of it, with correlations between e.g. the galactic and GMC scale gas pressures, surface densities and volume densities. That way, the galactic environment sets the initial conditions for star formation within GMCs. After the onset of massive star formation, stellar feedback from e.g. photoionisation, stellar winds, and supernovae eventually contributes to dispersing the parent cloud, depositing energy, momentum and metals into the surrounding medium, thereby changing the properties of galaxies. This cycling of matter between gas and stars, governed by star formation and feedback, is therefore a major driver of galaxy evolution. Much of the recent debate has focused on the durations of the various evolutionary phases that constitute this cycle in galaxies, and what these can teach us about the physical mechanisms driving the cycle. We review results from observational, theoretical, and numerical work to build a dynamical picture of the evolutionary lifecycle of GMC evolution, star formation, and feedback in galaxies.

preprint2019arXiv

When H II Regions are Complicated: Considering Perturbations from Winds, Radiation Pressure, and Other Effects

We explore to what extent simple algebraic models can be used to describe H II regions when winds, radiation pressure, gravity and photon breakout are included. We a) develop algebraic models to describe the expansion of photoionised H II regions under the influence of gravity and accretion in power-law density fields with $ρ\propto r^{-w}$, b) determine when terms describing winds, radiation pressure, gravity and photon breakout become significant enough to affect the dynamics of the H II region where $w=2$, and c) solve these expressions for a set of physically-motivated conditions. We find that photoionisation feedback from massive stars is the principal mode of feedback on molecular cloud scales, driving accelerating outflows from molecular clouds in cases where the peaked density structure around young massive stars is considered at radii between $\sim$0.1 and 10-100 pc. Under a large range of conditions the effect of winds and radiation on the dynamics of H II regions is around 10% of the contribution from photoionisation. The effect of winds and radiation pressure are most important at high densities, either close to the star or in very dense clouds such as those in the Central Molecular Zone of the Milky Way. Out to $\sim$0.1 pc they are the principal drivers of the H II region. Lower metallicities make the relative effect of photoionisation even stronger as the ionised gas temperature is higher.

preprint2016arXiv

Distribution and mass of diffuse and dense CO gas in the Milky Way

Emission from carbon monoxide (CO) is ubiquitously used as a tracer of dense star forming molecular clouds. There is, however, growing evidence that a significant fraction of CO emission originates from diffuse molecular gas. Quantifying the contribution of diffuse CO-emitting gas is vital for understanding the relation between molecular gas and star formation. We examine the Galactic distribution of two CO-emitting gas components, a high column density component detected in 13CO and 12CO, and a low column density component detected in 12CO, but not in 13CO. The "diffuse" and "dense" components are identified using a combination of smoothing, masking, and erosion/dilation procedures, making use of three large-scale 12CO and 13CO surveys of the Inner and Outer Milky Way. The diffuse component, which globally represents 25 (1.5x1e8 Mo) of the total molecular gas mass (6.5x1e8 Mo), is more extended perpendicular to the Galactic plane. The fraction of diffuse gas increases from 15% at a galactocentric radius of 3 kpc to 50% at 15 kpc, and increases with decreasing surface density. In the Inner Galaxy, a yet denser component traced by CS emission represents 14% of the total molecular gas mass traced by 12CO emission. Only 14% of the molecular gas mass traced by 12CO emission is identified as part of molecular clouds in 12CO surveys by cloud identification algorithms. This study indicates that CO emission not only traces star forming clouds, but also a significant diffuse molecular ISM component.

preprint2016arXiv

Extreme Scale-out SuperMUC Phase 2 - lessons learned

In spring 2015, the Leibniz Supercomputing Centre (Leibniz-Rechenzentrum, LRZ), installed their new Peta-Scale System SuperMUC Phase2. Selected users were invited for a 28 day extreme scale-out block operation during which they were allowed to use the full system for their applications. The following projects participated in the extreme scale-out workshop: BQCD (Quantum Physics), SeisSol (Geophysics, Seismics), GPI-2/GASPI (Toolkit for HPC), Seven-League Hydro (Astrophysics), ILBDC (Lattice Boltzmann CFD), Iphigenie (Molecular Dynamic), FLASH (Astrophysics), GADGET (Cosmological Dynamics), PSC (Plasma Physics), waLBerla (Lattice Boltzmann CFD), Musubi (Lattice Boltzmann CFD), Vertex3D (Stellar Astrophysics), CIAO (Combustion CFD), and LS1-Mardyn (Material Science). The projects were allowed to use the machine exclusively during the 28 day period, which corresponds to a total of 63.4 million core-hours, of which 43.8 million core-hours were used by the applications, resulting in a utilization of 69%. The top 3 users were using 15.2, 6.4, and 4.7 million core-hours, respectively.

preprint2016arXiv

Modelling Dust Evolution in Galaxies with a Multiphase, Inhomogeneous ISM

We develop a model of dust evolution in a multiphase, inhomogeneous ISM including dust growth and destruction processes. The physical conditions for grain evolution are taken from hydrodynamical simulations of giant molecular clouds in a Milky Way-like spiral galaxy. We improve the treatment of dust growth by accretion in the ISM to investigate the role of the temperature-dependent sticking coefficient and ion-grain interactions. From detailed observational data on the gas-phase Si abundances [Si/H]_{gas} measured in the local Galaxy, we derive a relation between the average [Si/H]_{gas} and the local gas density n(H) which we use as a critical constraint for the models. This relation requires a sticking coefficient that decreases with the gas temperature. The synthetic relation constructed from the spatial dust distribution reproduces the slope of -0.5 of the observed relation in cold clouds. This slope is steeper than that for the warm medium and is explained by the dust growth. We find that it occurs for all adopted values of the minimum grain size a_{min} from 1 to 5nm. For the classical cut-off of a_{min}=5 nm, the ion-grain interactions result in longer growth timescales and higher [Si/H]_{gas} than the observed values. For a_{min} below 3 nm, the ion-grain interactions enhance the growth rates, steepen the slope of [Si/H]_{gas}-n(H) relation and provide a better match to observations. The rates of dust re-formation in the ISM by far exceed the rates of dust production by stellar sources as expected from simple evolution models. After the cycle of matter in and out of dust reaches a steady state, the dust growth balances the destruction operating on similar timescales of 350 Myr.

preprint2012arXiv

Magnetic Field Amplification by Small-Scale Dynamo Action: Dependence on Turbulence Models and Reynolds and Prandtl Numbers

The small-scale dynamo is a process by which turbulent kinetic energy is converted into magnetic energy, and thus is expected to depend crucially on the nature of turbulence. In this work, we present a model for the small-scale dynamo that takes into account the slope of the turbulent velocity spectrum v(l) ~ l^theta, where l and v(l) are the size of a turbulent fluctuation and the typical velocity on that scale. The time evolution of the fluctuation component of the magnetic field, i.e., the small-scale field, is described by the Kazantsev equation. We solve this linear differential equation for its eigenvalues with the quantum-mechanical WKB-approximation. The validity of this method is estimated as a function of the magnetic Prandtl number Pm. We calculate the minimal magnetic Reynolds number for dynamo action, Rm_crit, using our model of the turbulent velocity correlation function. For Kolmogorov turbulence (theta=1/3), we find that the critical magnetic Reynolds number is approximately 110 and for Burgers turbulence (theta=1/2) approximately 2700. Furthermore, we derive that the growth rate of the small-scale magnetic field for a general type of turbulence is Gamma ~ Re^((1-theta)/(1+theta)) in the limit of infinite magnetic Prandtl numbers. For decreasing magnetic Prandtl number (down to Pm approximately larger than 10), the growth rate of the small-scale dynamo decreases. The details of this drop depend on the WKB-approximation, which becomes invalid for a magnetic Prandtl number of about unity.

preprint2012arXiv

The Small-Scale Dynamo and Non-Ideal MHD in Primordial Star Formation

We study the amplification of magnetic fields during the formation of primordial halos. The turbulence generated by gravitational infall motions during the formation of the first stars and galaxies can amplify magnetic fields very efficiently and on short timescales up to dynamically significant values. Using the Kazantsev theory, which describes the so-called small-scale dynamo - a magnetohydrodynamical process converting kinetic energy from turbulence into magnetic energy - we can then calculate the growth rate of the small-scale magnetic field. Our calculations are based on a detailed chemical network and we include non-ideal magnetohydrodynamical effects such as ambipolar diffusion and Ohmic dissipation. We follow the evolution of the magnetic field up to larger scales until saturation occurs on the Jeans scale. Assuming a weak magnetic seed field generated by the Biermann battery process, both Burgers and Kolmogorov turbulence lead to saturation within a rather small density range. Such fields are likely to become relevant after the formation of a protostellar disk and, thus, could influence the formation of the first stars and galaxies in the Universe.

preprint2011arXiv

Importance of the Initial Conditions for Star Formation - III: Statistical Properties of Embedded Protostellar Clusters

We investigate the formation of protostellar clusters during the collapse of dense molecular cloud cores with a focus on the evolution of potential and kinetic energy, the degree of substructure, and the early phase of mass segregation. Our study is based on a series of hydrodynamic simulations of dense cores, where we vary the initial density profile and the initial turbulent velocity. In the three-dimensional adaptive mesh refinement simulations, we follow the dynamical formation of filaments and protostars until a star formation efficiency of 20%. Despite the different initial configurations, the global ensemble of all protostars in a setup shows a similar energy evolution and forms sub-virial clusters with an energy ratio $E_\mathrm{kin}/|E_\mathrm{pot}|\sim0.2$. Concentrating on the innermost central region, the clusters show a roughly virialised energy balance. However, the region of virial balance only covers the innermost $\sim10-30%$ of all the protostars. In all simulations with multiple protostars, the total kinetic energy of the protostars is higher than the kinetic energy of the gas cloud, although the protostars only contain 20% of the total mass. The clusters vary significantly in size, mass, and number of protostars, and show different degrees of substructure and mass segregation. Flat density profiles and compressive turbulent modes produce more subclusters then centrally concentrated profiles and solenoidal turbulence. We find that dynamical relaxation and hence dynamical mass segregation is very efficient in all cases from the very beginning of the nascent cluster, i.e., during a phase when protostars are constantly forming and accreting.

preprint2011arXiv

Simulations on a Moving Mesh: The Clustered Formation of Population III Protostars

The cosmic dark ages ended a few hundred million years after the Big Bang, when the first stars began to fill the universe with new light. It has generally been argued that these stars formed in isolation and were extremely massive - perhaps 100 times as massive as the Sun. In a recent study, Clark and collaborators showed that this picture requires revision. They demonstrated that the accretion disks that build up around Population III stars are strongly susceptible to fragmentation and that the first stars should therefore form in clusters rather than in isolation. We here use a series of high-resolution hydrodynamical simulations performed with the moving mesh code AREPO to follow up on this proposal and to study the influence of environmental parameters on the level of fragmentation. We model the collapse of five independent minihalos from cosmological initial conditions, through the runaway condensation of their central gas clouds, to the formation of the first protostar, and beyond for a further 1000 years. During this latter accretion phase, we represent the optically thick regions of protostars by sink particles. Gas accumulates rapidly in the circumstellar disk around the first protostar, fragmenting vigorously to produce a small group of protostars. After an initial burst, gravitational instability recurs periodically, forming additional protostars with masses ranging from ~ 0.1 to 10 M_sun. Although the shape, multiplicity, and normalization of the protostellar mass function depend on the details of the sink-particle algorithm, fragmentation into protostars with diverse masses occurs in all cases, confirming earlier reports of Population III stars forming in clusters. Depending on the efficiency of later accretion and merging, Population III stars may enter the main sequence in clusters and with much more diverse masses than are commonly assumed.

preprint2011arXiv

The Delay of Population III Star Formation by Supersonic Streaming Velocities

It has recently been demonstrated that coherent relative streaming velocities of order 30 km / s between dark matter and gas permeated the universe on scales below a few Mpc directly after recombination. We here use a series of high-resolution moving-mesh calculations to show that these supersonic motions significantly influence the virialization of the gas in minihalos, and delay the formation of the first stars. As the gas streams into minihalos with bulk velocities around 1 km / s at z ~ 20, the additional momentum and energy input reduces the gas fractions and central densities of the halos, increasing the typical virial mass required for efficient cooling by a factor of three, and delaying Population III star formation by dz ~ 4. Since the distribution of the magnitude of the streaming velocities is narrowly peaked around a non-negligible value, this effect is important in most regions of the universe. As a consequence, the increased minimum halo mass implies a reduction of the absolute number of minihalos that can be expected to cool and form Population III stars by up to an order of magnitude. We further find that the streaming velocities increase the turbulent velocity dispersion of the minihalo gas, which could affect its ability to fragment and hence alter the mass function of the first stars.

preprint2011arXiv

The Turbulence Spectrum of Molecular Clouds in the Galactic Ring Survey: A Density-Dependent PCA Calibration

Turbulence plays a major role in the formation and evolution of molecular clouds. The problem is that turbulent velocities are convolved with the density of an observed region. To correct for this convolution, we investigate the relation between the turbulence spectrum of model clouds, and the statistics of their synthetic observations obtained from Principal Component Analysis (PCA). We apply PCA to spectral maps generated from simulated density and velocity fields, obtained from hydrodynamic simulations of supersonic turbulence, and from fractional Brownian motion fields with varying velocity, density spectra, and density dispersion. We examine the dependence of the slope of the PCA structure function, alpha_PCA, on intermittency, on the turbulence velocity (beta_v) and density (beta_n) spectral indexes, and on density dispersion. We find that PCA is insensitive to beta_n and to the log-density dispersion sigma_s, provided sigma_s < 2. For sigma_s>2, alpha_PCA increases with sigma_s due to the intermittent sampling of the velocity field by the density field. The PCA calibration also depends on intermittency. We derive a PCA calibration based on fBms with sigma_s<2 and apply it to 367 CO spectral maps of molecular clouds in the Galactic Ring Survey. The average slope of the PCA structure function, <alpha_PCA>=0.62\pm0.2, is consistent with the hydrodynamic simulations and leads to a turbulence velocity exponent <beta_v>=2.06\pm0.6 for a non-intermittent, low density dispersion flow. Accounting for intermittency and density dispersion, the coincidence between the PCA slope of the GRS clouds and the hydrodynamic simulations suggests beta_v~1.9, consistent with both Burgers and compressible intermittent turbulence.

preprint2010arXiv

Effects of Varying the Three-Body Molecular Hydrogen Formation Rate in Primordial Star Formation

The transformation of atomic hydrogen to molecular hydrogen through three-body reactions is a crucial stage in the collapse of primordial, metal-free halos, where the first generation of stars (Population III stars) in the Universe are formed. However, in the published literature, the rate coefficient for this reaction is uncertain by nearly an order of magnitude. We report on the results of both adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) simulations of the collapse of metal-free halos as a function of the value of this rate coefficient. For each simulation method, we have simulated a single halo three times, using three different values of the rate coefficient. We find that while variation between halo realizations may be greater than that caused by the three-body rate coefficient being used, both the accretion physics onto Population III protostars as well as the long-term stability of the disk and any potential fragmentation may depend strongly on this rate coefficient.

preprint2009arXiv

Clump morphology and evolution in MHD simulations of molecular cloud formation

Abridged: We study the properties of clumps formed in three-dimensional weakly magnetized magneto-hydrodynamic simulations of converging flows in the thermally bistable, warm neutral medium (WNM). We find that: (1) Similarly to the situation in the classical two-phase medium, cold, dense clumps form through dynamically-triggered thermal instability in the compressed layer between the convergent flows, and are often characterised by a sharp density jump at their boundaries though not always. (2) However, the clumps are bounded by phase-transition fronts rather than by contact discontinuities, and thus they grow in size and mass mainly by accretion of WNM material through their boundaries. (3) The clump boundaries generally consist of thin layers of thermally unstable gas, but these layers are often widened by the turbulence, and penetrate deep into the clumps. (4) The clumps are approximately in both ram and thermal pressure balance with their surroundings, a condition which causes their internal Mach numbers to be comparable to the bulk Mach number of the colliding WNM flows. (5) The clumps typically have mean temperatures 20 < T < 50 K, corresponding to the wide range of densities they contain (20 < n < 5000 pcc) under a nearly-isothermal equation of state. (6) The turbulent ram pressure fluctuations of the WNM induce density fluctuations that then serve as seeds for local gravitational collapse within the clumps. (7) The velocity and magnetic fields tend to be aligned with each other within the clumps, although both are significantly fluctuating, suggesting that the velocity tends to stretch and align the magnetic field with it. (8) The typical mean field strength in the clumps is a few times larger than that in the WNM. (9) The magnetic field strength has a mean value of B ~ 6 mu G ...

preprint2001arXiv

The Formation of Stellar Clusters: Mass Spectra from Turbulent Molecular Cloud Fragmentation

Turbulent fragmentation determines where and when protostellar cores form, and how they contract and grow in mass from the surrounding cloud material. This process is investigated, using numerical models of self-gravitating molecular cloud dynamics. Molecular cloud regions without turbulent driving sources, or where turbulence is driven on large scales, exhibit rapid and efficient star formation in a clustered mode, whereas interstellar turbulence that carries most energy on small scales results in isolated star formation with low efficiency. The clump mass spectrum of shock-generated density fluctuations in pure hydrodynamic, supersonic turbulence is not well fit by a power law, and it is too steep at the high-mass end to be in agreement with the observational data. When gravity is included in the turbulence models, local collapse occurs, and the spectrum extends towards larger masses as clumps merge together, a power-law description dN/dM ~ M^nu becomes possible with slope nu < -2. In the case of pure gravitational contraction, i.e. in regions without turbulent support, the clump mass spectrum is shallower with nu = -3/2. The mass spectrum of protostellar cores in regions without turbulent support and where turbulence is replenished on large-scales, however, is well described by a log-normal or by multiple power laws, similar to the stellar IMF at low and intermediate masses. In the case of small-scale turbulence, the core mass spectrum is too flat compared to the IMF for all masses.

preprint1997arXiv

GRAPESPH with Fully Periodic Boundary Conditions: Fragmentation of Molecular Clouds

A method of adapting smoothed particle hydrodynamics (SPH) with periodic boundary conditions for use with the special purpose device GRAPE is presented. GRAPE (GRAvity PipE) solves the Poisson and force equations for an N-body system by direct summation on a specially designed chip and in addition returns the neighbour list for each particle. Due to its design, GRAPE cannot treat periodic particle distributions directly. This limitation of GRAPESPH can be overcome by computing a correction force for each particle due to periodicity (Ewald correction) on the host computer using a PM-like method. This scheme is applied to study the fragmentation process in giant molecular clouds. Assuming a pure isothermal model, we follow the dynamical evolution in the interior of a molecular cloud starting from an Gaussian initial density distribution to the formation of selfgravitating clumps until most of the gas is consumed in these dense cores. Despite its simplicity, this model can reproduce some fundamental properties of observed molecular clouds, like a clump mass distribution of the form $dN/dm \propto m^n$, with $n ~ -1.5$.

preprint1995arXiv

Constraints on massive black holes as dark matter candidates using Galactic globular clusters

This work considers the idea of massive black holes being the constituents of the Galactic dark matter halo. It constrains the maximum black hole mass to $\mbh \sil 5 \times 10^4 \msun$ by examining their influence on the population of globular clusters in our Milky Way. In the adopted halo model, globular clusters are exposed to constant bombardment of halo objects on their orbits through the Galaxy and thus will steadily gain internal energy. Depending on the mass of these halo objects and the structural parameters of the globular clusters, they can be disrupted on time scales of a few billion years and below. These disruption time scales are calculated using a modification of the well known (classical) impulsive approximation and compared with direct N-body simulations of such encounter events to ensure the method works correctly. For a set of ten prototypical globular cluster models and black hole masses ranging from $10^3$ to $10^7\:\msun$, Monte-Carlo-simulations of $10\,000$ encounter histories over the period of 10 billion years were calculated each, at three different galactocentric distances \mbox{$R = 5$}, $10$ and $15\, \kpc$. These data were compared with the real globular cluster population in our Galaxy and used to obtain the above constraint of $\mbh \sil 5 \times 10^4 \msun$.