BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Paul C. Clark

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

ResearcherUnclaimed source record
astro-ph.GAastro-ph.SRastro-ph.COastro-phastro-ph.HEastro-ph.IM

Catalog footprint

What is connected

48works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

48 published item(s)

preprint2023arXiv

Towards the impact of GMC collisions on the star formation rate

Collisions between giant molecular clouds (GMCs) are one of the pathways for massive star formation, due to the high densities created. However the enhancement of the star formation rate (SFR) is not well constrained. In this study we perform a parameter study of cloud-cloud collisions, and investigate how the resulting SFR depends on the details of set-up. Our parameter study explores variations in: collision speed; magnetic field inclination (with respect to the collisional axis); and resolution, as defined by the number of cells per Jeans length. In all our collision simulations we find a factor of 2-3 increase in the SFR compared to our no collision simulation, with star formation beginning sooner with a) high collisional velocities, b) parallel orientation between the magnetic field and collision axis, c) and lower resolution. The mean virial parameter of high density (and thus possible star-forming) gas increases with collisional velocity, but has little variation with magnetic field inclination. The alignment of the velocity and magnetic field remains uniform in low density environments but becomes more perpendicular with increasing density, indicating the compression of the magnetic field by collapsing gas. Comparing the trends in the SFR with other GMC collision studies, we find good agreement with studies that account for the gravitational boundedness of the gas in their star formation algorithm, but not with those that simply form stars above a prescribed density threshold. This suggests that the latter approach should be used with caution when modelling star formation on resolved cloud scales.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Fragmentation induced starvation in Population III star formation: a resolution study

The Population III initial mass function (IMF) is currently unknown, but recent studies agree that fragmentation of primordial gas gives a broader IMF than the initially suggested singular star per halo. In this study we introduce sink particle mergers into Arepo, to perform the first resolution study for primordial star formation simulations and present the first Population III simulations to run up to densities of 10-6g cm-3 for hundreds of years after the formation of sink particles. The total number of sinks formed increases with increasing sink particle creation density, without achieving numerical convergence. The total mass in sinks remains invariant to the maximum resolution and is safe to estimate using low resolution studies. This results in an IMF that shifts towards lower masses with increasing resolution. Greater numbers of sinks cause increased fragmentation-induced starvation of the most massive sink, yielding lower accretion rates, masses and ionising photons emitted per second. The lack of convergence up to densities 2 orders of magnitudes higher than all relevant chemical reactions suggests that the number of sinks will continue to grow with increasing resolution until H2 is fully dissociated and the collapse becomes almost adiabatic at 10-4g cm-3. These results imply that many Population III studies utilising sink particles have produced IMFs which have overestimated the masses of primordial stars, and underestimated the number of stars formed. In the highest resolution runs, sinks with masses capable of surviving until the present day had an ejection fraction of 0.21.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

High contrast and resolution near infrared photometry of the core of R136

We present the sharpest and deepest near infrared photometric analysis of the core of R136, a newly formed massive star cluster at the centre of the 30 Doradus star forming region in the Large Magellanic Cloud. We used the extreme adaptive optics of the SPHERE focal instrument implemented on the ESO Very Large Telescope and operated in its IRDIS imaging mode, for the second time with longer exposure time in the H- and K filters. Our aim was to (i) increase the number of resolved sources in the core of R136, and (ii) to compare with the first epoch to classify the properties of the detected common sources between the two epochs. Within the field of view (FOV) of 10.8"x12.1" (2.7pc x3.0pc), we detected 1499 sources in both H and K filters, for which 76% of these sources have visual companions closer than 0.2". The larger number of detected sources, enabled us to better sample the mass function (MF). The MF slopes are estimated at ages of 1, 1.5 and 2 Myr, at different radii, and for different mass ranges. The MF slopes for the mass range of 10-300 solar-mass are about 0.3 dex steeper than the mass range of 3-300 solar-mass, for the whole FOV and different radii. Comparing the JHK colours of 790 sources common in between the two epochs, 67% of detected sources in the outer region (r >3") are not consistent with evolutionary models at 1-2 Myr and with extinctions similar to the average cluster value, suggesting an origin from ongoing star formation within 30 Doradus, unrelated to R136.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Formation sites of Population III star formation: The effects of different levels of rotation and turbulence on the fragmentation behavior of primordial gas

We use the moving-mesh code AREPO to investigate the effects of different levels of rotation and turbulence on the fragmentation of primordial gas and the formation of Population III stars. We consider 9 different combinations of turbulence and rotation and carry out 5 different realizations of each setup, yielding one of the largest sets of simulations of Population III star formation ever performed. We find that fragmentation in Population III star-forming systems is a highly chaotic process and show that the outcomes of individual realizations of the same initial conditions often vary significantly. However, some general trends are apparent. Increasing the turbulent energy promotes fragmentation, while increasing the rotational energy inhibits fragmentation. Within the 1000 yr period that we simulate, runs including turbulence yield flat protostellar mass functions while purely rotational runs show a more top-heavy distribution. The masses of the individual protostars are distributed over a wide range from a few $10^{-3} \, {\rm M_{\odot}}$ to several tens of ${\rm M_\odot}$. The total mass growth rate of the stellar systems remains high throughout the simulations and depends only weakly on the degree of rotation and turbulence. Mergers between protostars are common, but predictions of the merger fraction are highly sensitive to the criterion used to decide whether two protostars should merge. Previous studies of Population III star formation have often considered only one realization per set of initial conditions. However, our results demonstrate that robust trends can only be reliably identified by considering averages over a larger sample of runs.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Simulations of the star-forming molecular gas in an interacting M51-like galaxy

We present here the first of a series of papers aimed at better understanding the evolution and properties of giant molecular clouds (GMCs) in a galactic context. We perform high resolution, three-dimensional {\sc arepo} simulations of an interacting galaxy inspired by the well-observed M51 galaxy. Our fiducial simulations include a non-equilibrium, time-dependent, chemical network that follows the evolution of atomic and molecular hydrogen as well as carbon and oxygen self-consistently. Our calculations also treat gas self-gravity and subsequent star formation (described by sink particles), and coupled supernova feedback. In the densest parts of the simulated interstellar medium (ISM) we reach sub-parsec resolution, granting us the ability to resolve individual GMCs and their formation and destruction self-consistently throughout the galaxy. In this initial work we focus on the general properties of the ISM with a particular focus on the cold star-forming gas. We discuss the role of the interaction with the companion galaxy in generating cold molecular gas and controlling stellar birth. We find that while the interaction drives large-scale gas flows and induces spiral arms in the galaxy, it is of secondary importance in determining gas fractions in the different ISM phases and the overall star-formation rate. The behaviour of the gas on small GMC scales instead is mostly controlled by the self-regulating property of the ISM driven by coupled feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Cloud Factory I: Generating resolved filamentary molecular clouds from galactic-scale forces

We introduce a new suite of simulations, "The Cloud Factory", which self-consistently forms molecular cloud complexes at high enough resolution to resolve internal substructure (up to 0.25 Msol in mass) all while including galactic-scale forces. We use a version of the Arepo code modified to include a detailed treatment of the physics of the cold molecular ISM, and an analytical galactic gravitational potential for computational efficiency. The simulations have nested levels of resolution, with the lowest layer tied to tracer particles injected into individual cloud complexes. These tracer refinement regions are embedded in the larger simulation so continue to experience forces from outside the cloud. This allows the simulations to act as a laboratory for testing the effect of galactic environment on star formation. Here we introduce our method and investigate the effect of galactic environment on filamentary clouds. We find that cloud complexes formed after a clustered burst of feedback, have shorter lengths and are less likely to fragment compared to quiescent clouds (e.g. the Musca filament) or those dominated by the galactic potential (e.g. Nessie). Spiral arms and differential rotation preferentially align filaments, but strong feedback randomises them. Long filaments formed within the cloud complexes are necessarily coherent with low internal velocity gradients, which has implications for the formation of filamentary star-clusters. Cloud complexes formed in regions dominated by supernova feedback have fewer star-forming cores, and these are more widely distributed. These differences show galactic-scale forces can have a significant impact on star formation within molecular clouds.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Launching cosmic-ray-driven outflows from the magnetized interstellar medium

We present a hydrodynamical simulation of the turbulent, magnetized, supernova (SN)-driven interstellar medium (ISM) in a stratified box that dynamically couples the injection and evolution of cosmic rays (CRs) and a self-consistent evolution of the chemical composition. CRs are treated as a relativistic fluid in the advection-diffusion approximation. The thermodynamic evolution of the gas is computed using a chemical network that follows the abundances of H+, H, H2, CO, C+, and free electrons and includes (self-)shielding of the gas and dust. We find that CRs perceptibly thicken the disk with the heights of 90% (70%) enclosed mass reaching ~1.5 kpc (~0.2 kpc). The simulations indicate that CRs alone can launch and sustain strong outflows of atomic and ionized gas with mass loading factors of order unity, even in solar neighborhood conditions and with a CR energy injection per SN of 10^50 erg, 10% of the fiducial thermal energy of an SN. The CR-driven outflows have moderate launching velocities close to the midplane (~100 km/s) and are denser (ρ~1e-24 - 1e-26 g/cm^3), smoother, and colder than the (thermal) SN-driven winds. The simulations support the importance of CRs for setting the vertical structure of the disk as well as the driving of winds.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Star formation and molecular hydrogen in dwarf galaxies: a non-equilibrium view

We study the connection of star formation to atomic (HI) and molecular hydrogen (H$_2$) in isolated, low metallicity dwarf galaxies with high-resolution ($m_{\rm gas}$ = 4 M$_\odot$, $N_{\rm ngb}$ = 100) SPH simulations. The model includes self-gravity, non-equilibrium cooling, shielding from an interstellar radiation field, the chemistry of H$_2$ formation, H$_2$-independent star formation, supernova feedback and metal enrichment. We find that the H$_2$ mass fraction is sensitive to the adopted dust-to-gas ratio and the strength of the interstellar radiation field, while the star formation rate is not. Star formation is regulated by stellar feedback, keeping the gas out of thermal equilibrium for densities $n <$ 1 cm$^{-3}$. Because of the long chemical timescales, the H$_2$ mass remains out of chemical equilibrium throughout the simulation. Star formation is well-correlated with cold ( T $\leqslant$ 100 K ) gas, but this dense and cold gas - the reservoir for star formation - is dominated by HI, not H$_2$. In addition, a significant fraction of H$_2$ resides in a diffuse, warm phase, which is not star-forming. The ISM is dominated by warm gas (100 K $<$ T $\leqslant 3\times 10^4$ K) both in mass and in volume. The scale height of the gaseous disc increases with radius while the cold gas is always confined to a thin layer in the mid-plane. The cold gas fraction is regulated by feedback at small radii and by the assumed radiation field at large radii. The decreasing cold gas fractions result in a rapid increase in depletion time (up to 100 Gyrs) for total gas surface densities $Σ_{\rm HI+H_2} \lesssim$ 10 M$_\odot$pc$^{-2}$, in agreement with observations of dwarf galaxies in the Kennicutt-Schmidt plane.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The IMF as a function of supersonic turbulence

Recent studies seem to suggest that the stellar initial mass function (IMF) in early-type galaxies might be different from a classical Kroupa or Chabrier IMF, i.e. contain a larger fraction of the total mass in low-mass stars. From a theoretical point of view, supersonic turbulence has been the subject of interest in many analytical theories proposing a strong correlation with the characteristic mass of the core mass function (CMF) in star forming regions, and as a consequence with the stellar IMF. Performing two suites of smoothed particles hydrodynamics (SPH) simulations with different mass resolutions, we aim at testing the effects of variations in the turbulent properties of a dense, star forming molecular cloud on the shape of the system mass function in different density regimes. While analytical theories predict a shift of the peak of the CMF towards lower masses with increasing velocity dispersion of the cloud, we observe in the low-density regime the opposite trend, with high Mach numbers giving rise to a top-heavy mass distribution. For the high-density regime we do not find any trend correlating the Mach number with the characteristic mass of the resulting IMF, implying that the dynamics of protostellar accretion discs and fragmentation on small scales is not strongly affected by turbulence driven at the scale of the cloud. Furthermore, we suggest that a significant fraction of dense cores are disrupted by turbulence before stars can be formed in their interior through gravitational collapse. Although this particular study has limitations in its numerical resolution, we suggest that our results, along with those from other studies, cast doubt on the turbulent fragmentation models on the IMF that simply map the CMF to the IMF.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The role of cosmic ray pressure in accelerating galactic outflows

We study the formation of galactic outflows from supernova explosions (SNe) with the moving-mesh code AREPO in a stratified column of gas with a surface density similar to the Milky Way disk at the solar circle. We compare different simulation models for SNe placement and energy feedback, including cosmic rays (CR), and find that models that place SNe in dense gas and account for CR diffusion are able to drive outflows with similar mass loading as obtained from a random placement of SNe with no CRs. Despite this similarity, CR-driven outflows differ in several other key properties including their overall clumpiness and velocity. Moreover, the forces driving these outflows originate in different sources of pressure, with the CR diffusion model relying on non-thermal pressure gradients to create an outflow driven by internal pressure and the random-placement model depending on kinetic pressure gradients to propel a ballistic outflow. CRs therefore appear to be non-negligible physics in the formation of outflows from the interstellar medium.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Using CO line ratios to trace the physical properties of molecular clouds

The carbon monoxide (CO) rotational transition lines are the most common tracers of molecular gas within giant molecular clouds (MCs). We study the ratio ($R_{2-1/1-0}$) between CO's first two emission lines and examine what information it provides about the physical properties of the cloud. To study $R_{2-1/1-0}$ we perform smooth particle hydrodynamic simulations with time dependent chemistry (using GADGET-2), along with post-process radiative transfer calculations on an adaptive grid (using RADMC-3D) to create synthetic emission maps of a MC. $R_{2-1/1-0}$ has a bimodal distribution that is a consequence of the excitation properties of each line, given that $J=1$ reaches local thermal equilibrium (LTE) while $J=2$ is still sub-thermally excited in the considered clouds. The bimodality of $R_{2-1/1-0}$ serves as a tracer of the physical properties of different regions of the cloud and it helps constrain local temperatures, densities and opacities. Additionally this bimodal structure shows an important portion of the CO emission comes from diffuse regions of the cloud, suggesting that the commonly used conversion factor of $R_{2-1/1-0}\sim 0.7$ between both lines may need to be studied further.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Does the CO-to-H2 conversion factor depend on the star formation rate?

We present a series of numerical simulations that explore how the `X-factor', $X_{CO}$ -- the conversion factor between the observed integrated CO emission and the column density of molecular hydrogen -- varies with the environmental conditions in which a molecular cloud is placed. Our investigation is centred around two environmental conditions in particular: the cosmic ray ionisation rate (CRIR) and the strength of the interstellar radiation field (ISRF). Since both these properties of the interstellar medium have their origins in massive stars, we make the assumption in this paper that both the strength of the ISRF and the CRIR scale linearly with the local star formation rate (SFR). The cloud modelling in this study first involves running numerical simulations that capture the cloud dynamics, as well as the time-dependent chemistry, and ISM heating and cooling. These simulations are then post-processed with a line radiative transfer code to create synthetic 12CO (1-0) emission maps from which $X_{CO}$ can be calculated. We find that for 1e4 solar mass virialised clouds with mean density 100 cm$^{-3}$, $X_{CO}$ is only weakly dependent on the local SFR, varying by a factor of a few over two orders of magnitude in SFR. In contrast, we find that for similar clouds but with masses of 1e5 solar masses, the X-factor will vary by an order of magnitude over the same range in SFR, implying that extra-galactic star formation laws should be viewed with caution. However, for denser ($10^4$ cm$^{-3}$), super-virial clouds such as those found at the centre of the Milky Way, the X-factor is once again independent of the local SFR.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hunting for the signatures of molecular cloud formation

In order to understand how molecular clouds form in the Galactic interstellar medium, we would like to be able to map the structure and kinematics of the gas flows responsible for forming them. However, doing so is observationally challenging. CO, the workhorse molecule for studies of molecular clouds, traces only relatively dense gas and hence only allows us to study those portions of the clouds that have already assembled. Numerical simulations suggest that the inflowing gas that forms these clouds is largely composed of CO-dark H2. These same simulations allow us to explore the usefulness of different tracers of this CO-dark molecular material, and we use them here to show that the [CII] fine structure line is potentially a very powerful tracer of this gas and should be readily detectable using modern instrumentation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Impact of supernova and cosmic-ray driving on the surface brightness of the galactic halo in soft X-rays

The halo of the Milky Way contains a hot plasma with a surface brightness in soft X-rays of the order $10^{-12}$erg cm$^{-2}$ s$^{-1}$ deg$^{-2}$. The origin of this gas is unclear, but so far numerical models of galactic star formation have failed to reproduce such a large surface brightness by several orders of magnitude. In this paper, we analyze simulations of the turbulent, magnetized, multi-phase interstellar medium including thermal feedback by supernova explosions as well as cosmic-ray feedback. We include a time-dependent chemical network, self-shielding by gas and dust, and self-gravity. Pure thermal feedback alone is sufficient to produce the observed surface brightness, although it is very sensitive to the supernova rate. Cosmic rays suppress this sensitivity and reduce the surface brightness because they drive cooler outflows. Self-gravity has by far the largest effect because it accumulates the diffuse gas in the disk in dense clumps and filaments, so that supernovae exploding in voids can eject a large amount of hot gas into the halo. This can boost the surface brightness by several orders of magnitude. Although our simulations do not reach a steady state, all simulations produce surface brightness values of the same order of magnitude as the observations, with the exact value depending sensitively on the simulation parameters. We conclude that star formation feedback alone is sufficient to explain the origin of the hot halo gas, but measurements of the surface brightness alone do not provide useful diagnostics for the study of galactic star formation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Is atomic carbon a good tracer of molecular gas in metal-poor galaxies?

Carbon monoxide (CO) is widely used as a tracer of molecular hydrogen (H2) in metal-rich galaxies, but is known to become ineffective in low metallicity dwarf galaxies. Atomic carbon has been suggested as a superior tracer of H2 in these metal-poor systems, but its suitability remains unproven. To help us to assess how well atomic carbon traces H2 at low metallicity, we have performed a series of numerical simulations of turbulent molecular clouds that cover a wide range of different metallicities. Our simulations demonstrate that in star-forming clouds, the conversion factor between [CI] emission and H2 mass, $X_{\rm CI}$, scales approximately as $X_{\rm CI} \propto Z^{-1}$. We recover a similar scaling for the CO-to-H2 conversion factor, $X_{\rm CO}$, but find that at this point in the evolution of the clouds, $X_{\rm CO}$ is consistently smaller than $X_{\rm CI}$, by a factor of a few or more. We have also examined how $X_{\rm CI}$ and $X_{\rm CO}$ evolve with time. We find that $X_{\rm CI}$ does not vary strongly with time, demonstrating that atomic carbon remains a good tracer of H2 in metal-poor systems even at times significantly before the onset of star formation. On the other hand, $X_{\rm CO}$ varies very strongly with time in metal-poor clouds, showing that CO does not trace H2 well in starless clouds at low metallicity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Modelling [CI] emission from turbulent molecular clouds

We use detailed numerical simulations of the coupled chemical, thermal and dynamical evolution of the gas in a turbulent molecular cloud to study the usefulness of the [CI] 609 micron and 370 micron fine structure emission lines as tracers of cloud structure. Emission from these lines is observed throughout molecular clouds, and yet the question of what we can learn from them about the physics of the clouds remains largely unexplored. We show that the fact that [CI] emission is widespread within molecular clouds is a simple consequence of the fact that the clouds are dominated by turbulent motions. Turbulence creates large density inhomogeneities, allowing radiation to penetrate deeply into the clouds. As a result, [CI] emitting gas is found throughout the cloud, rather than being concentrated at the edges. We examine how well we can use [CI] emission to trace the structure of the cloud, and show that the integrated intensity of the 609 micron line traces column density accurately over a wide range of visual extinctions. For extinctions greater than a few, [CI] and 13CO both perform well, but [CI] becomes a superior tracer of column densities for visual extinctions A_V <= 3 We have also studied the distribution of [CI] excitation temperatures in the gas, and show that these are typically smaller than the kinetic temperature, indicating that most of the carbon atoms are not in local thermodynamic equilibrium. We discuss how best to determine T_ex from observations of the [CI] lines, and how to use these values to estimate the column density of neutral atomic carbon. We show that even in the best case, we tend to systematically underestimate the atomic carbon content of the gas. Our results suggest that observationally-derived estimates of the atomic carbon content of real GMCs could be in error by as much as a factor of two.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Star formation efficiencies of molecular clouds in a galactic center environment

We use the Arepo moving mesh code to simulate the evolution of molecular clouds exposed to a harsh environment similar to that found in the galactic center (GC), in an effort to understand why the star formation efficiency (SFE) of clouds in this environment is so small. Our simulations include a simplified treatment of time-dependent chemistry and account for the highly non-isothermal nature of the gas and the dust. We model clouds with a total mass of 1.3x10^5 M_{sun} and explore the effects of varying the mean cloud density and the virial parameter, alpha = E_{kin}/|E_{pot}|. We vary the latter from alpha = 0.5 to alpha = 8.0, and so many of the clouds that we simulate are gravitationally unbound. We expose our model clouds to an interstellar radiation field (ISRF) and cosmic ray flux (CRF) that are both a factor of 1000 higher than the values found in the solar neighbourhood. As a reference, we also run simulations with local solar neighbourhood values of the ISRF and the CRF in order to better constrain the effects of the extreme conditions in the GC on the SFE. Despite the harsh environment and the large turbulent velocity dispersions adopted, we find that all of the simulated clouds form stars within less than a gravitational free-fall time. Increasing the virial parameter from alpha = 0.5 to alpha = 8.0 decreases the SFE by a factor ~4-10, while increasing the ISRF/CRF by a factor of 1000 decreases the SFE again by a factor ~2-6. However, even in our most unbound clouds, the SFE remains higher than that inferred for real GC clouds. We therefore conclude that high levels of turbulence and strong external heating are not enough by themselves to lead to a persistently low SFE at the center of the Galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Synthetic observations of molecular clouds in a galactic center environment: I. Studying maps of column density and integrated intensity

We run numerical simulations of molecular clouds (MCs), adopting properties similar to those found in the Central Molecular Zone (CMZ) of the Milky Way. For this, we employ the moving mesh code Arepo and perform simulations which account for a simplified treatment of time-dependent chemistry and the non-isothermal nature of gas and dust. We perform simulations using an initial density of n_0 = 10^3 cm^{-3} and a mass of 1.3x10^5 M_sun. Furthermore, we vary the virial parameter, defined as the ratio of kinetic and potential energy, alpha = E_{kin} / |E_{pot}|, by adjusting the velocity dispersion. We set it to alpha = 0.5, 2.0 and 8.0, in order to analyze the impact of the kinetic energy on our results. We account for the extreme conditions in the CMZ and increase both the interstellar radiation field (ISRF) and the cosmic-ray flux (CRF) by a factor of 1000 compared to the values found in the solar neighbourhood. We use the radiative transfer code RADMC-3D to compute synthetic images in various diagnostic lines. These are [CII] at 158 micron, [OI] (145 micron), [OI] (63 micron), 12CO (J = 1 -> 0) and 13CO (J = 1 -> 0) at 2600 micron and 2720 micron, respectively. When alpha is large, the turbulence disperses much of the gas in the cloud, reducing its mean density and allowing the ISRF to penetrate more deeply into the cloud's interior. This significantly alters the chemical composition of the cloud, leading to the dissociation of a significant amount of the molecular gas. On the other hand, when alpha is small, the cloud remains compact, allowing more of the molecular gas to survive. We show that in each case the atomic tracers accurately reflect most of the physical properties of both the H2 and the total gas of the cloud and that they provide a useful alternative to molecular lines when studying the ISM in the CMZ.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The role of 3-body H$_2$ formation in the fragmentation of primordial gas

It has been shown that the behaviour of primordial gas collapsing in a dark matter minihalo can depend on the adopted choice of 3-body H$_2$ formation rate. The uncertainties in this rate span two orders of magnitude in the current literature, and so it remains a source of uncertainty in our knowledge of population III star formation. Here we investigate how the amount of fragmentation in primordial gas depends on the adopted 3-body rate. We present the results of calculations that follow the chemical and thermal evolution of primordial gas as it collapses in two dark matter minihalos. Our results on the effect of 3-body rate on the evolution until the first protostar forms agree well with previous studies. However, our modified version of GADGET-2 SPH also includes sink particles, which allows us to follow the initial evolution of the accretion disc that builds up on the centre of each halo, and capture the fragmentation in gas as well as its dependence on the adopted 3-body H$_2$ formation rate. We find that the fragmentation behaviour of the gas is only marginally effected by the choice of 3-body rate co-efficient, and that halo-to-halo differences are of equal importance in affecting the final mass distribution of stars.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The SILCC (SImulating the LifeCycle of molecular Clouds) project - II. Dynamical evolution of the supernova-driven ISM and the launching of outflows

The SILCC project (SImulating the Life-Cycle of molecular Clouds) aims at a more self-consistent understanding of the interstellar medium (ISM) on small scales and its link to galaxy evolution. We present three-dimensional (magneto)hydrodynamic simulations of the ISM in a vertically stratified box including self-gravity, an external potential due to the stellar component of the galactic disc, and stellar feedback in the form of an interstellar radiation field and supernovae (SNe). The cooling of the gas is based on a chemical network that follows the abundances of H+, H, H2, C+, and CO and takes shielding into account consistently. We vary the SN feedback by comparing different SN rates, clustering and different positioning, in particular SNe in density peaks and at random positions, which has a major impact on the dynamics. Only for random SN positions the energy is injected in sufficiently low-density environments to reduce energy losses and enhance the effective kinetic coupling of the SNe with the gas. This leads to more realistic velocity dispersions (σ_HI ~ 0.8σ_(300-8000K) ~ 10-20km/s, σ_Hα~ 0.6σ_(8000-3e5K) ~ 20-30km/s), and strong outflows with mass loading factors of up to 10 even for solar neighbourhood conditions. Clustered SNe abet the onset of outflows compared to individual SNe but do not influence the net outflow rate. The outflows do not contain any molecular gas and are mainly composed of atomic hydrogen. The bulk of the outflowing mass is dense (ρ~ 1e-25-1e-24g/cc) and slow (v ~ 20-40km/s) but there is a high-velocity tail of up to v ~ 500km/s with ρ~ 1e-28-1e-27g/cc.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A new approach to determine optically thick H2 cooling and its effect on primordial star formation

We present a new method for estimating the H2 cooling rate in the optically thick regime in simulations of primordial star formation. Our new approach is based on the TreeCol algorithm, which projects matter distributions onto a spherical grid to create maps of column densities for each fluid element in the computational domain. We have improved this algorithm by using the relative gas velocities, to weight the individual matter contributions with the relative spectral line overlaps, in order to properly account for the Doppler effect. We compare our new method to the widely used Sobolev approximation, which yields an estimate for the column density based on the local velocity gradient and the thermal velocity. This approach generally underestimates the photon escape probability, because it neglects the density gradient and the actual shape of the cloud. We present a correction factor for the true line overlap in the Sobolev approximation and a new method based on local quantities, which fits the exact results reasonably well during the collapse of the cloud, with the error in the cooling rates always being less than 10%. Analytical fitting formulae fail at determining the photon escape probability after formation of the first protostar (error of 40%) because they are based on the assumption of spherical symmetry and therefore break down once a protostellar accretion disc has formed. Our method yields lower temperatures and hence promotes fragmentation for densities above 10^{10}/ccm at a distance of 200AU from the first protostar. Since the overall accretion rates are hardly affected by the cooling implementation, we expect Pop III stars to have lower masses in our simulations, compared to the results of previous simulations that used the Sobolev approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

CO-dark gas and molecular filaments in Milky Way type galaxies

We use the moving mesh code AREPO coupled to a time-dependent chemical network to investigate the formation and destruction of molecular gas in simulated spiral galaxies. This allows us to determine the characteristics of the gas that is not traced by CO emission. Our extremely high resolution AREPO simulations allow us to capture the chemical evolution of the disc, without recourse to a parameterised `clumping factor'. We calculate H2 and CO column densities through our simulated disc galaxies, and estimate the CO emission and CO-H2 conversion factor. We find that in conditions akin to those in the local interstellar medium, around 42% of the total molecular mass should be in CO-dark regions, in reasonable agreement with observational estimates. This fraction is almost insensitive to the CO integrated intensity threshold used to discriminate between CO-bright and CO-dark gas, as long as this threshold is less than 10 K km/s. The CO-dark molecular gas primarily resides in extremely long (>100 pc) filaments that are stretched between spiral arms by galactic shear. Only the centres of these filaments are bright in CO, suggesting that filamentary molecular clouds observed in the Milky Way may only be small parts of much larger structures. The CO-dark molecular gas mainly exists in a partially molecular phase which accounts for a significant fraction of the total disc mass budget. The dark gas fraction is higher in simulations with higher ambient UV fields or lower surface densities, implying that external galaxies with these conditions might have a greater proportion of dark gas.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Interpreting the sub-linear Kennicutt-Schmidt relationship: The case for diffuse molecular gas

Recent statistical analysis of two extragalactic observational surveys strongly indicate a sublinear Kennicutt-Schmidt (KS) relationship between the star formation rate (Sigsfr) and molecular gas surface density (Sigmol). Here, we consider the consequences of these results in the context of common assumptions, as well as observational support for a linear relationship between Sigsfr and the surface density of dense gas. If the CO traced gas depletion time (tau_mol) is constant, and if CO only traces star forming giant molecular clouds (GMCs), then the physical properties of each GMC must vary, such as the volume densities or star formation rates. Another possibility is that the conversion between CO luminosity and Sigmol, the XCO factor, differs from cloud-to-cloud. A more straightforward explanation is that CO permeates the hierarchical ISM, including the filaments and lower density regions within which GMCs are embedded. A number of independent observational results support this description, with the diffuse gas comprising at least 30% of the total molecular content. The CO bright diffuse gas can explain the sublinear KS relationship, and consequently leads to an increasing tau_mol with Sigmol. If Sigsfr linearly correlates with the dense gas surface density, a sublinear KS relationship indicates that the fraction of diffuse gas fdiff grows with Sigmol. In galaxies where Sigmol falls towards the outer disk, this description suggests that fdiff also decreases radially.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On column density thresholds and the star formation rate

We present the results of a numerical study designed to address the question of whether there is a column density threshold for star formation within molecular clouds. We have simulated a large number of different clouds, with volume and column densities spanning a wide range of different values, using a state-of-the-art model for the coupled chemical, thermal and dynamical evolution of the gas. We show that star formation is only possible in regions where the mean (area-averaged) column density exceeds $10^{21} \: {\rm cm^{-2}}$. Within the clouds, we also show that there is a good correlation between the mass of gas above a K-band extinction $A_{\rm K} = 0.8$ and the star formation rate (SFR), in agreement with recent observational work. Previously, this relationship has been explained in terms of a correlation between the SFR and the mass in dense gas. However, we find that this correlation is weaker and more time-dependent than that between the SFR and the column density. In support of previous studies, we argue that dust shielding is the key process: the true correlation is one between the SFR and the mass in cold, well-shielded gas, and the latter correlates better with the column density than the volume density.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Statistical properties of dark matter mini-haloes at z >= 15

Understanding the formation of the first objects in the universe critically depends on knowing whether the properties of small dark matter structures at high-redshift (z > 15) are different from their more massive lower-redshift counterparts. To clarify this point, we performed a high-resolution N-body simulation of a cosmological volume 1 Mpc/h comoving on a side, reaching the highest mass resolution to date in this regime. We make precision measurements of various physical properties that characterize dark matter haloes (such as the virial ratio, spin parameter, shape, and formation times, etc.) for the high-redshift (z > 15) dark matter mini-haloes we find in our simulation, and compare them to literature results and a moderate-resolution comparison run within a cube of side-length 100 Mpc/h. We find that dark matter haloes at high-redshift have a log-normal distribution of the dimensionless spin parameter centered around λ $\sim$ 0.03, similar to their more massive counterparts. They tend to have a small ratio of the length of the shortest axis to the longest axis (sphericity), and are highly prolate. In fact, haloes of given mass that formed recently are the least spherical, have the highest virial ratios, and have the highest spins. Interestingly, the formation times of our mini-halos depend only very weakly on mass, in contrast to more massive objects. This is expected from the slope of the linear power spectrum of density perturbations at this scale, but despite this difference, dark matter structures at high-redshift share many properties with their much more massive counterparts observed at later times.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Origin and Universality of the Stellar Initial Mass Function

We review current theories for the origin of the Stellar Initial Mass Function (IMF) with particular focus on the extent to which the IMF can be considered universal across various environments. To place the issue in an observational context, we summarize the techniques used to determine the IMF for different stellar populations, the uncertainties affecting the results, and the evidence for systematic departures from universality under extreme circumstances. We next consider theories for the formation of prestellar cores by turbulent fragmentation and the possible impact of various thermal, hydrodynamic and magneto-hydrodynamic instabilities. We address the conversion of prestellar cores into stars and evaluate the roles played by different processes: competitive accretion, dynamical fragmentation, ejection and starvation, filament fragmentation and filamentary accretion flows, disk formation and fragmentation, critical scales imposed by thermodynamics, and magnetic braking. We present explanations for the characteristic shapes of the Present-Day Prestellar Core Mass Function and the IMF and consider what significance can be attached to their apparent similarity. Substantial computational advances have occurred in recent years, and we review the numerical simulations that have been performed to predict the IMF directly and discuss the influence of dynamics, time-dependent phenomena, and initial conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Indications of a sub-linear and non-universal Kennicutt-Schmidt relationship

We estimate the parameters of the Kennicutt-Schmidt (KS) relationship, linking the star formation rate (Sigma_SFR) to the molecular gas surface density (Sigma_mol), in the STING sample of nearby disk galaxies using a hierarchical Bayesian method. This method rigorously treats measurement uncertainties, and provides accurate parameter estimates for both individual galaxies and the entire population. Assuming standard conversion factors to estimate Sigma_SFR and Sigma_mol from the observations, we find that the KS parameters vary between galaxies, indicating that no universal relationship holds for all galaxies. The KS slope of the whole population is 0.76, with the 2sigma range extending from 0.58 to 0.94. These results imply that the molecular gas depletion time is not constant, but varies from galaxy to galaxy, and increases with the molecular gas surface density. Therefore, other galactic properties besides just Sigma_mol affect Sigma_SFR, such as the gas fraction or stellar mass. The non-universality of the KS relationship indicates that a comprehensive theory of star formation must take into account additional physical processes that may vary from galaxy to galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Formation and evolution of primordial protostellar systems

We investigate the formation of the first stars at the end of the cosmic dark ages with a suite of three-dimensional, moving mesh simulations that directly resolve the collapse of the gas beyond the formation of the first protostar at the centre of a dark matter minihalo. The simulations cover more than 25 orders of magnitude in density and have a maximum spatial resolution of 0.05 R_sun, which extends well below the radius of individual protostars and captures their interaction with the surrounding gas. In analogy to previous studies that employed sink particles, we find that the Keplerian disc around the primary protostar fragments into a number of secondary protostars, which is facilitated by H2 collisional dissociation cooling and collision-induced emission. The further evolution of the protostellar system is characterized by strong gravitational torques that transfer angular momentum between the secondary protostars formed in the disc and the surrounding gas. This leads to the migration of about half of the secondary protostars to the centre of the cloud in a free-fall time, where they merge with the primary protostar and enhance its growth to about five times the mass of the second most massive protostar. By the same token, a fraction of the protostars obtain angular momentum from other protostars via N-body interactions and migrate to higher orbits. On average, only every third protostar survives until the end of the simulation. However, the number of protostars present at any given time increases monotonically, suggesting that the system will continue to grow beyond the limited period of time simulated here.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

How long does it take to form a molecular cloud?

We present the first numerical simulations that self-consistently follow the formation of dense molecular clouds in colliding flows. Our calculations include a time-dependent model for the H2 and CO chemistry that runs alongside a detailed treatment of the dominant heating and cooling processes in the ISM. We adopt initial conditions characteristic of the warm neutral medium and study two different flow velocities - a slow flow with v = 6.8 km/s and a fast flow with v = 13.6 km/s. The clouds formed by the collision of these flows form stars, with star formation beginning after 16 Myr in the case of the slower flow, but after only 4.4 Myr in the case of the faster flow. In both flows, the formation of CO-dominated regions occurs only around 2 Myr before the onset of star formation. Prior to this, the clouds produce very little emission in the J = 1 -> 0 transition line of CO, and would probably not be identified as molecular clouds in observational surveys. In contrast, our models show that H2-dominated regions can form much earlier, with the timing depending on the details of the flow. In the case of the slow flow, small pockets of gas become fully molecular around 10 Myr before star formation begins, while in the fast flow, the first H2-dominated regions occur around 3 Myr before the first prestellar cores form. Our results are consistent with models of molecular cloud formation in which the clouds are dominated by "dark" molecular gas for a considerable proportion of their assembly history.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the effects of rotation during the formation of population III protostars

It has been suggested that turbulent motions are responsible for the transport of angular momentum during the formation of Population III stars, however the exact details of this process have never been studied. We report the results from three dimensional SPH simulations of a rotating self-gravitating primordial molecular cloud, in which the initial velocity of solid-body rotation has been changed. We also examine the build-up of the discs that form in these idealized calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the formation of very metal-poor stars: The case of SDSS J1029151+172927

The formation of stars is a key process in the early universe with far reaching consequences for further cosmic evolution. While stars forming from truly primordial gas are thought to be considerably more massive than our Sun, stars in the universe today have typical masses below one solar mass. The physical origin of this transition and the conditions under which it occurs are highly debated. There are two competing models, one based on metal-line cooling as the primary agent and one based on dust cooling. The recent discovery of the extremely metal poor star SDSS J1029151+172927 provides a unique opportunity to distinguish between these two models. Based on simple thermodynamic considerations we argue that SDSS J1029151+172927 was more likely formed as a result of dust continuum cooling rather than cooling by metal lines. We conclude that the masses of extremely metal-poor stars are determined by dust-induced fragmentation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Star formation in metal-poor gas clouds

Observations of molecular clouds in metal-poor environments typically find that they have much higher star formation rates than one would expect based on their observed CO luminosities and the molecular gas masses that are inferred from them. This finding can be understood if one assumes that the conversion factor between CO luminosity and H2 mass is much larger in these low metallicity systems than in nearby molecular clouds. However, it is unclear whether this is the only factor at work, or whether the star formation rate of the clouds is directly sensitive to the metallicity of the gas. To investigate this, we have performed numerical simulations of the coupled dynamical, chemical and thermal evolution of model clouds with metallicities ranging from 0.01 Z_solar to Z_solar. We find that the star formation rate in our model clouds has little sensitivity to the metallicity. Reducing the metallicity of the gas by two orders of magnitude delays the onset of star formation in the clouds by no more than a cloud free-fall time and reduces the time-averaged star formation rate by at most a factor of two. On the other hand, the chemical state of the clouds is highly sensitive to the metallicity, and at the lowest metallicities, the clouds are completely dominated by atomic gas. Our results confirm that the CO-to-H2 conversion factor in these systems depends strongly on the metallicity, but also show that the precise value is highly time-dependent, as the integrated CO luminosity of the most metal-poor clouds is dominated by emission from short-lived gravitationally collapsing regions. Finally, we find evidence that the star formation rate per unit H2 mass increases with decreasing metallicity, owing to the much smaller H2 fractions present in our low metallicity clouds.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The stellar IMF at very low metallicities

The theory for the formation of the first population of stars (Pop III) predicts an initial mass function (IMF) dominated by high-mass stars, in contrast to the present-day IMF, which tends to yield mostly stars with masses less than 1 M_Sol. The leading theory for the transition in the characteristic stellar mass predicts that the cause is the extra cooling provided by increasing metallicity. In particular, dust can overtake H_2 as the leading coolant at very high densities. The aim of this work is to determine the influence of dust cooling on the fragmentation of very low metallicity gas. To investigate this, we make use of high-resolution hydrodynamic simulations with sink particles to replace contracting protostars, and analyze the collapse and further fragmentation of star-forming clouds. We follow the thermodynamic response of the gas by solving the full thermal energy equation, and also track the behavior of the dust temperature and the chemical evolution of the gas. We model four clouds with different metallicities (10^{-4}, 10^{-5}, 10^{-6} Z_Sol, and 0), and determine the properties of each cloud at the point at which it undergoes gravitational fragmentation. We find evidence for fragmentation in all four cases, and hence conclude that there is no critical metallicity below which fragmentation is impossible. Nevertheless, there is a clear change in the behavior of the clouds at Z = 10^{-5} Z_Sol, caused by the fact that at this metallicity, fragmentation takes longer to occur than accretion, leading to a flat mass function at lower metallicities.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Protostellar outflows with Smoothed Particle Magnetohydrodynamics (SPMHD)

The protostellar collapse of a molecular cloud core is usually accompanied by outflow phenomena. The latter are thought to be driven by magnetorotational processes from the central parts of the protostellar disc. While several 3D AMR/nested grid studies of outflow phenomena in collapsing magnetically supercritical dense cores have been reported in the literature, so far no such simulation has been performed using the Smoothed Particle Hydrodynamics (SPH) method. This is mainly due to intrinsic numerical difficulties in handling magnetohydrodynamics within SPH, which only recently were partly resolved. In this work, we use an approach where we evolve the magnetic field via the induction equation, augmented with stability correction and divergence cleaning schemes. We consider the collapse of a rotating core of one solar mass, threaded by a weak magnetic field initially parallel to the rotation axis so that the core is magnetically supercritical. We show, that Smoothed Particle Magnetohydrodynamics (SPMHD) is able to handle the magnetorotational processes connected with outflow phenomena, and to produce meaningful results which are in good agreement with findings reported in the literature. Especially, our numerical scheme allows for a quantitative analysis of the evolution of the ratio of the toroidal to the poloidal magnetic field, which we performed in this work.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The effect of dust cooling on low-metallicity star-forming clouds

The theory for the formation of the first population of stars (Pop III) predicts a IMF composed predominantly of high-mass stars, in contrast to the present-day IMF, which tends to yield stars with masses less than 1 M_Solar. The leading theory for the transition in the characteristic stellar mass predicts that the cause is the extra cooling provided by increasing metallicity and in particular the cooling provided at high densities by dust. The aim of this work is to test whether dust cooling can lead to fragmentation and be responsible for this transition. To investigate this, we make use of high-resolution hydrodynamic simulations. We follow the thermodynamic evolution of the gas by solving the full thermal energy equation, and also track the evolution of the dust temperature and the chemical evolution of the gas. We model clouds with different metallicities, and determine the properties of the cloud at the point at which it undergoes gravitational fragmentation. We follow the further collapse to scales of an AU when we replace very dense, gravitationally bound, and collapsing regions by a simple and nongaseous object, a sink particle. Our results suggest that for metallicities as small as 10^{-5}Z_Solar, dust cooling produces low-mass fragments and hence can potentially enable the formation of low mass stars. We conclude that dust cooling affects the fragmentation of low-metallicity gas clouds and plays an important role in shaping the stellar IMF even at these very low metallicities. We find that the characteristic fragment mass increases with decreasing metallicity, but find no evidence for a sudden transition in the behaviour of the IMF within the range of metallicites examined in our present study.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Effects of Accretion Luminosity upon Fragmentation in the Early Universe

We introduce a prescription for the luminosity from accreting protostars into smoothed particle hydrodynamics simulation, and apply the method to simulations of five primordial minihalos generated from cosmological initial conditions. We find that accretion luminosity delays fragmentation within the halos, but does not prevent it. In halos that slowly form a low number of protostars, the accretion luminosity can reduce the number of fragments that are formed before the protostars start ionising their surroundings. However, halos that rapidly form many protostars become dominated by dynamical processes, and the effect of accretion luminosity becomes negligible. Generally the fragmentation found in the halos is highly dependent on the initial conditions. Accretion luminosity does not substantially affect the accretion rates experienced by the protostars, and is far less important than dynamical interactions, which can lead to ejections that effectively terminate the accretion. We find that the accretion rates onto the inner regions of the disks (20 AU) around the protostars are highly variable, in contrast to the constant or smoothly decreasing accretion rates currently used in models of the pre-main sequence evolution of Population III stars.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Formation and Fragmentation of Disks around Primordial Protostars

The very first stars to form in the Universe heralded an end to the cosmic dark ages and introduced new physical processes that shaped early cosmic evolution. Until now, it was thought that these stars lived short, solitary lives, with only one extremely massive star, or possibly a very wide binary system, forming in each dark matter minihalo. Here we describe numerical simulations that show that these stars were, to the contrary, often members of tight multiple systems. Our results show that the disks that formed around the first young stars were unstable to gravitational fragmentation, possibly producing small binary and higher-order systems that had separations as small as the distance between the Earth and the Sun.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

TreeCol: a novel approach to estimating column densities in astrophysical simulations

We present TreeCol, a new and efficient tree-based scheme to calculate column densities in numerical simulations. Knowing the column density in any direction at any location in space is a prerequisite for modelling the propagation of radiation through the computational domain. TreeCol therefore forms the basis for a fast, approximate method for modelling the attenuation of radiation within large numerical simulations. It constructs a HEALPix sphere at any desired location and accumulates the column density by walking the tree and by adding up the contributions from all tree nodes whose line of sight contributes to the pixel under consideration. In particular when combined with widely-used tree-based gravity solvers the new scheme requires little additional computational cost. In a simulation with $N$ resolution elements, the computational cost of TreeCol scales as $N \log N$, instead of the $N^{5/3}$ scaling of most other radiative transfer schemes. TreeCol is naturally adaptable to arbitrary density distributions and is easy to implement and to parallelize. We discuss its accuracy and performance characteristics for the examples of a spherical protostellar core and for the turbulent interstellar medium. We find that the column density estimates provided by TreeCol are on average accurate to better than 10 percent. In another application, we compute the dust temperatures for solar neighborhood conditions and compare with the result of a full-fledged Monte Carlo radiation-transfer calculation. We find that both methods give very similar answers. We conclude that TreeCol provides a fast, easy to use, and sufficiently accurate method of calculating column densities that comes with little additional computational cost when combined with an existing tree-based gravity solver.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A quantification of the non-spherical geometry and accretion of collapsing cores

We present the first detailed classification of the structures of Class 0 cores in a high resolution simulation of a giant molecular cloud. The simulated cloud contains 10^4 solar masses and produces over 350 cores which allows for meaningful statistics. Cores are classified into three types according to how much they depart from spherical symmetry. We find that three quarters of the cores are better described as irregular filaments than as spheres. Recent Herschel results have shown that cores are formed within a network of filaments, which we find has had a significant impact on the resulting core geometries. We show that the column densities and ram pressure seen by the protostar are not uniform and generally peak along the axes of the filament. The angular momentum vector of the material in the cores varies both in magnitude and direction, which will cause the rotation vector of the central source to fluctuate during the collapse of the core. In the case of the more massive stars, accretion from the environment outside the original core volume is even more important than that from the core itself. This additional gas is primarily accreted onto the cores along the dense filaments in which the cores are embedded, and the sections of the surfaces of the cores which do not coincide with a filament have very little additional material passing through them. The assumption of spherical symmetry cannot be applied to the majority of collapsing cores, and is never a good description of how stars accrete gas from outside the original core radius. This has ramifications for our understanding of collapsing cores, in particular their line profiles, the effect of radiation upon them and their ability to fragment.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Gravitational fragmentation in turbulent primordial gas and the initial mass function of Population III stars

We report results from numerical simulations of star formation in the early universe that focus on the dynamical behavior of metal-free gas under different initial and environmental conditions. In particular we investigate the role of turbulence, which is thought to ubiquitously accompany the collapse of high-redshift halos. We distinguish between two main cases: the birth of Population III.1 stars - those which form in the pristine halos unaffected by prior star formation - and the formation of Population III.2 stars - those forming in halos where the gas is still metal free but has an increased ionization fraction. This latter case can arise either from exposure to the intense UV radiation of stellar sources in neighboring halos, or from the high virial temperatures associated with the formation of massive halos, that is, those with masses greater than 1e8 solar masses. We find that turbulent primordial gas is highly susceptible to fragmentation in both cases, even for turbulence in the subsonic regime, i.e. for rms velocity dispersions as low as 20 % of the sound speed. Contrary to our original expectations, fragmentation is more vigorous and more widespread in pristine halos compared to pre-ionized ones. We therefore predict Pop III.1 stars to be on average of somewhat lower mass, and form in larger groups, than Pop III.2 stars. We find that fragment masses cover over two orders of magnitude, indicating that the resulting Population III initial mass function was significantly extended in mass as well. This prompts the need for a large, high-resolution study of the formation of dark matter minihalos that is capable of resolving the turbulent flows in the gas at the moment when the baryons become self-gravitating. This would help determine which, if any, of the initial conditions presented in our study are realized in nature.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Modeling Collapse and Accretion in Turbulent Gas Clouds: Implementation and Comparison of Sink Particles in AMR and SPH

We implemented sink particles in the adaptive mesh refinement (AMR) hydrodynamics code FLASH. Sink particles are created in regions of local gravitational collapse, and their trajectories and accretion can be followed over many dynamical times. We perform a series of tests including the time integration of circular and elliptical orbits, the collapse of a Bonnor-Ebert sphere and a rotating, fragmenting cloud core. We compare the collapse of a highly unstable singular isothermal sphere to the theory by Shu (1977), and show that the sink particle accretion rate is in excellent agreement with the theoretical prediction. To model eccentric orbits and close encounters of sink particles accurately, we show that a very small timestep is often required, for which we implemented subcycling of the N-body system. We emphasize that a sole density threshold for sink particle creation is insufficient in supersonic flows, if the density threshold is below the opacity limit. In that case, the density can exceed the threshold in strong shocks that do not necessarily lead to local collapse. Additional checks for bound state, gravitational potential minimum, Jeans instability and converging flows are absolutely necessary for a meaningful creation of sink particles. We apply our new sink particle module for FLASH to the formation of a stellar cluster, and compare to a smoothed particle hydrodynamics (SPH) code with sink particles. Our comparison shows encouraging agreement of gas properties, indicated by column density distributions and radial profiles, and of sink particle formation times and positions. We find excellent agreement in the number of sink particles formed, and in their accretion and mass distributions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Protostellar collapse and fragmentation using an MHD GADGET

Although the influence of magnetic fields is regarded as vital in the star formation process, only a few magnetohydrodynamics (MHD) simulations have been performed on this subject within the smoothed particle hydrodynamics (SPH) method. This is largely due to the unsatisfactory treatment of non-vanishing divergence of the magnetic field. Recently smoothed particle magnetohydrodynamics (SPMHD) simulations based on Euler potentials have proven to be successful in treating MHD collapse and fragmentation problems, however these methods are known to have some intrinsical difficulties. We have performed SPMHD simulations based on a traditional approach evolving the magnetic field itself using the induction equation. To account for the numerical divergence, we have chosen an approach that subtracts the effects of numerical divergence from the force equation, and additionally we employ artificial magnetic dissipation as a regularization scheme. We apply this realization of SPMHD to a widely known setup, a variation of the 'Boss & Bodenheimer standard isothermal test case', to study the impact of the magnetic fields on collapse and fragmentation. In our simulations, we concentrate on setups, where the initial magnetic field is parallel to the rotation axis. We examine different field strengths and compare our results to other findings reported in the literature. We are able to confirm specific results found elsewhere, namely the delayed onset of star formation for strong fields, accompanied by the tendency to form only single stars. We also find that the 'magnetic cushioning effect', where the magnetic field is wound up to form a 'cushion' between the binary, aids binary fragmentation in a case, where previously only formation of a single protostar was expected.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stellar and substellar initial mass function: a model that implements gravoturbulent fragmentation and accretion

In this work, we derive the stellar initial mass function (IMF) from the superposition of mass distributions of dense cores, generated through gravoturbulent fragmentation of unstable clumps in molecular clouds (MCs) and growing through competitive accretion. MCs are formed by the turbulent cascade in the interstellar medium at scales L from 100 down to ~0.1 pc. Their internal turbulence is essentially supersonic and creates clumps with a lognormal distribution of densities n. Our model is based on the assumption of a power-law relationship between clump mass and clump density: n~m^x, where x is a scale-free parameter. Gravitationally unstable clumps are assumed to undergo isothermal fragmentation and produce protostellar cores with a lognormal mass distribution, centred around the clump Jeans mass. Masses of individual cores are then assumed to grow further through competitive accretion until the rest of the gas within the clump is being exhausted. The observed IMF is best reproduced for a choice of x=0.25, for a characteristic star formation timescale of ~5 Myr, and for a low star formation efficiency of ~10 %.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Effect of Environment on Massive Star Formation

In this contribution we review our recent numerical work discussing the essential role of the local cluster environment in assembling massive stars. First we show that massive stars are formed from low mass pre-stellar cores and become massive due to accretion. Proto-stars that benefit from this accretion are those situated at the centre of a cluster's potential well, which is the focal point of the contraction of the cluster gas. Given that most of the mass which makes up a massive star in this model comes from the cluster environment rather than the core, it is important to model the molecular cloud environment accurately. Preliminary results of a simulation which accurately treats the chemistry and time-dependent thermodynamics of a molecular cloud show quantitatively similar star formation to previous models, but allow a true comparison to be made between simulation and observations. This method can also be applied to cases with varying metallicities allowing star formation in primordial gas to be studied. In general, these numerical studies of clustered star formation yield IMFs which are compatible with the Salpeter mass function. The only possible exception to this is in low density unbound regions of molecular clouds which lack very low and high mass stars.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The efficiency of star formation in clustered and distributed regions

We investigate the formation of both clustered and distributed populations of young stars in a single molecular cloud. We present a numerical simulation of a 10,000 solar mass elongated, turbulent, molecular cloud and the formation of over 2500 stars. The stars form both in stellar clusters and in a distributed mode which is determined by the local gravitational binding of the cloud. A density gradient along the major axis of the cloud produces bound regions that form stellar clusters and unbound regions that form a more distributed population. The initial mass function also depends on the local gravitational binding of the cloud with bound regions forming full IMFs whereas in the unbound, distributed regions the stellar masses cluster around the local Jeans mass and lack both the high-mass and the low-mass stars. The overall efficiency of star formation is ~ 15 % in the cloud when the calculation is terminated, but varies from less than 1 % in the the regions of distributed star formation to ~ 40 % in regions containing large stellar clusters. Considering that large scale surveys are likely to catch clouds at all evolutionary stages, estimates of the (time-averaged) star formation efficiency for the giant molecular cloud reported here is only ~ 4 %. This would lead to the erroneous conclusion of 'slow' star formation when in fact it is occurring on a dynamical timescale.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Fragmentation in Molecular Clouds and its connection to the IMF

We present an analysis of star-forming gas cores in an SPH simulation of a Giant Molecular Cloud. We identify cores using their deep potential wells. This yields a smoother distribution with clearer boundaries than density. Additionally, this gives an indication of future collapse, as bound potential cores (p-cores) represent the earliest stages of fragmentation in molecular clouds. We find that the mass function of the p-cores resembles the stellar IMF and the observed clump mass function, although p-core masses (~0.7 Msol) are smaller than typical density clumps. The bound p-cores are generally subsonic, have internal substructure, and are only quasi-spherical. We see no evidence of massive bound cores supported by turbulence. We trace the evolution of the p-cores forward in time, and investigate the connection between the original p-core mass and the stellar mass that formed from it. We find that there is a poor correlation, with considerable scatter suggesting accretion onto the core is dependent on more factors than just the initial core mass. During the accretion process the p-cores accrete from beyond the region first bound, highlighting the importance of the core environment to its subsequent evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Formation of Stellar Clusters and the Importance of Thermodynamics for Fragmentation

We discuss results from numerical simulations of star cluster formation in the turbulent interstellar medium (ISM). The thermodynamic behavior of the star-forming gas plays a crucial role in fragmentation and determines the stellar mass function as well as the dynamic properties of the nascent stellar cluster. This holds for star formation in molecular clouds in the solar neighborhood as well as for the formation of the very first stars in the early universe. The thermodynamic state of the ISM is a result of the balance between heating and cooling processes, which in turn are determined by atomic and molecular physics and by chemical abundances. Features in the effective equation of state of the gas, such as a transition from a cooling to a heating regime, define a characteristic mass scale for fragmentation and so set the peak of the initial mass function of stars (IMF). As it is based on fundamental physical quantities and constants, this is an attractive approach to explain the apparent universality of the IMF in the solar neighborhood as well as the transition from purely primordial high-mass star formation to the more normal low-mass mode observed today.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Star formation in unbound giant molecular clouds: the origin of OB associations?

We investigate the formation of star clusters in an unbound GMC, where the supporting kinetic energy is twice as large as the cloud's self-gravity. This cloud manages to form a series of star clusters and disperse, all within roughly 2 crossing times (10 Myr), supporting recent claims that star formation is a rapid process. Simple assumptions about the nature of the star formation occurring in the clusters allows us to place an estimate for the star formation efficiency at about 5 to 10 %, consistent with observations. We also propose that unbound clouds can act as a mechanism for forming OB associations. The clusters that form in the cloud behave as OB subgroups. These clusters are naturally expanding from one another due to unbound nature of the flows that create them. The properties of the cloud we present here are are consistent with those of classic OB associations.

0 citations0 reviewsNo code linkedOpen access