The streaming instability is a leading candidate mechanism to explain the formation of planetesimals. Yet, the role of this instability in the driving of turbulence in protoplanetary disks, given its fundamental nature as a linear hydrodynamical instability, has so far not been investigated in detail. We study the turbulence that is induced by the streaming instability as well as its interaction with the vertical shear instability. For this purpose, we employ the FLASH Code to conduct two-dimensional axisymmetric global disk simulations spanning radii from $1$ au to $100$ au, including the mutual drag between gas and dust as well as the radial and vertical stellar gravity. If the streaming instability and the vertical shear instability start their growth at the same time, we find the turbulence in the dust mid-plane layer to be primarily driven by the streaming instability. It gives rise to vertical gas motions with a Mach number of up to ${\sim}10^{-2}$. The dust scale height is set in a self-regulatory manner to about $1\%$ of the gas scale height. In contrast, if the vertical shear instability is allowed to saturate before the dust is introduced into our simulations, then it continues to be the main source of the turbulence in the dust layer. The vertical shear instability induces turbulence with a Mach number of ${\sim}10^{-1}$ and thus impedes dust sedimentation. Nonetheless, we find the vertical shear instability and the streaming instability in combination to lead to radial dust concentration in long-lived accumulations which are significantly denser than those formed by the streaming instability alone. Thus, the vertical shear instability may promote planetesimal formation by creating weak overdensities that act as seeds for the streaming instability.

Magnetic fields are an elemental part of the interstellar medium in galaxies. However, their impact on gas dynamics and star formation in galaxies remains controversial. We use a suite of global magnetohydrodynamical simulations of isolated disc galaxies to study the influence of magnetic fields on the diffuse and dense gas in the discs. We find that the magnetic field acts in multiple ways. Stronger magnetised discs fragment earlier due to the shorter growth time of the Parker instability. Due to the Parker instability in the magnetised discs we also find cold ($T<50\,\mathrm{K}$) and dense ($n\sim10^3-10^4\,\mathrm{cm}^{-3}$) gas several hundred pc above/below the midplane without any form of stellar feedback. In addition, magnetic fields change the fragmentation pattern. While in the hydrodynamical case, the disc breaks up into ring-like structures, magnetised discs show the formation of filamentary entities that extent both in the azimuthal and radial direction. These kpc scale filaments become magnetically (super-)critical very quickly and allow for the rapid formation of massive giant molecular clouds. Our simulations suggest that major differences in the behaviour of star formation - due to a varying magnetisation - in galaxies could arise.

We present radiation-magnetohydrodynamic simulations aimed at studying evolutionary properties of H\,{\normalsize II} regions in turbulent, magnetised, and collapsing molecular clouds formed by converging flows in the warm neutral medium. We focus on the structure, dynamics and expansion laws of these regions. Once a massive star forms in our highly structured clouds, its ionising radiation eventually stops the accretion (through filaments) toward the massive star-forming regions. The new over-pressured H\,{\normalsize II} regions push away the dense gas, thus disrupting the more massive collapse centres. Also, because of the complex density structure in the cloud, the H\,{\normalsize II} regions expand in a hybrid manner: they virtually do not expand toward the densest regions (cores), while they expand according to the classical analytical result towards the rest of the cloud, and in an accelerated way, as a blister region, towards the diffuse medium. Thus, the ionised regions grow anisotropically, and the ionising stars generally appear off-centre of the regions. Finally, we find that the hypotheses assumed in standard H\,{\normalsize II}-region expansion models (fully embedded region, blister-type, or expansion in a density gradient) apply simultaneously in different parts of our simulated H\,{\normalsize II} regions, producing a net expansion law ($R \propto t^\alpha$, with $\alpha$ in the range of 0.93-1.47 and a mean value of $1.2 \pm 0.17$) that differs from any of those of the standard models.

Both observational and theoretical research over the past decade has demonstrated that the probability distribution function (PDF) of the gas density in turbulent molecular clouds is a key ingredient for understanding star formation. It has recently been argued that the PDF of molecular clouds is a pure power-law distribution. It has been claimed that the log-normal part is ruled out when using only the part of the PDF up/down to which it is complete, that is where the column density contours are still closed. By using the results from high-resolution magnetohydrodynamical simulations of molecular cloud formation and evolution, we find that the column density PDF is indeed composed of a log-normal and, if including self-gravity, a power-law part. We show that insufficient sampling of a molecular cloud results in closed contours that cut off the log-normal part. In contrast, systematically increasing the field of view and sampling the entire cloud yields a completeness limit at the lower column densities, which also recovers the log-normal part. This demonstrates that the field of view must be sufficiently large for the PDF to be complete down to its log-normal part, which has important implications for predictions of star-formation activity based on the PDF.

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

Observations show that galaxies and their interstellar media are pervaded by strong magnetic fields with energies in the diffuse component being at least comparable to the thermal and even as large or larger than the turbulent energy. Such strong magnetic fields prevent the formation of stars because patches of the interstellar medium are magnetically subcritical. Here we present the results from global numerical simulations of strongly magnetised and self-gravitating galactic discs, which show that the buoyancy of the magnetic field due to the Parker instability leads at first to the formation of giant filamentary regions. These filamentary structures become gravitationally unstable and fragment into $\sim10^5 M_{\odot}$ clouds that attract kpc long, coherent filamentary flows that build them into GMCs. Our results thus provide a solution to the long-standing problem of how the transition from sub- to supercritical regions in the interstellar medium proceeds.

The origin of cosmic magnetic fields is an unsolved problem and magnetogenesis could have occurred in the early Universe. We study the evolution of such primordial magnetic fields across the cosmological recombination epoch via 3D magnetohydrodynamic numerical simulations. We compute the effective or net heating rate of baryons due to decaying magnetic fields and its dependence on the magnetic field strength and spectral index. In the drag-dominated regime ($z \gtrsim 1500$), prior to recombination, we find no real heating is produced. Our simulations allow us to smoothly trace a new transition regime ($600 \lesssim z \lesssim 1500$), where magnetic energy decays, at first, into the kinetic energy of baryons. A turbulent velocity field is built up until it saturates, as the net heating rate rises from a low value at recombination to its peak towards the end of the transition regime. This is followed by a turbulent decay regime ($z \lesssim 600$) where magnetic energy dissipates via turbulent decay of both magnetic and velocity fields while net heating remains appreciable and declines slowly. Both the peak of the net heating rate and the onset of turbulent decay are delayed significantly beyond recombination, by up to 0.5 Myr (until $z\simeq 600-700$), for scale-invariant magnetic fields. We provide analytic approximations and present numerical results for a range of field strengths and spectral indices, illustrating the redshift-dependence of dissipation and net heating rates. These can be used to study cosmic microwave background constraints on primordial magnetic fields.

Deuterium fractionation is considered as an important process to infer the chemical ages of prestellar cores in filaments. We present here the first magneto-hydrodynamical simulations including a chemical network to study deuterium fractionation in magnetized and turbulent filaments and their substructures. The filaments typically show widespread deuterium fractionation with average values $\gtrsim0.01$. For individual cores of similar age, we observe the deuteration fraction to increase with time, but also to be independent of their average properties such as density, virial or mass-to-magnetic flux ratio. We further find a correlation of the deuteration fraction with core mass, average H$_2$ density and virial parameter only at late evolutionary stages of the filament and attribute this to the lifetime of the individual cores. Specifically, chemically old cores reveal higher deuteration fractions. Within the radial profiles of selected cores, we notice differences in the structure of the deuteration fraction or surface density, which we can attribute to their different turbulent properties. High deuteration fractions of the order $0.01-0.1$ may be reached within approximately $200$~kyrs, corresponding to two free-fall times, as defined for cylindrical systems, of the filaments

We explore the phenomenon of the recently discovered inverse transfer of energy from small to large scales in decaying magnetohydrodynamical turbulence by Brandenburg et al. (2015) even for nonhelical magnetic fields. For this investigation we mainly employ the Pencil-Code performing a parameter study, where we vary the Prandtl number, the kinematic viscosity and the initial spectrum. We find that in order to get a decay which exhibits this inverse transfer, large Reynolds numbers ($\mathcal{O}\sim 10^{3}$) are needed and low Prandtl numbers of the order unity $Pr = 1$ are preferred. Compared to helical MHD turbulence, though, the inverse transfer is much less efficient in transferring magnetic energy to larger scales than the well-known effect of the inverse cascade. Hence, applying the inverse transfer to the magnetic field evolution in the Early Universe, we question whether the nonhelical inverse transfer is effective enough to explain the observed void magnetic fields if a magneto- genesis scenario during the electroweak phase transition is assumed.

Astrochemistry plays a central role during the process of star formation, both in the primordial regime as well as in the present-day Universe. We revisit here the chemistry in both regimes, focusing first on the chemistry under close to primordial conditions, as observed in the so-called Caffau star SDSS J102915+172927, and subsequently discuss deuteration processes in present-day star-forming cores. In models of the high-redshift Universe, the chemistry is particularly relevant to determine the cooling, while it also serves as an important diagnostic in the case of present-day star formation.

Molecular clouds are the principle stellar nurseries of our universe, keeping them in the focus of both observational and theoretical studies. From observations, some of the key properties of molecular clouds are well known but many questions regarding their evolution and star formation activity remain open. While numerical simulations feature a large number and complexity of involved physical processes, this plenty of effects may hide the fundamentals that determine the evolution of molecular clouds and enable the formation of stars. Purely analytical models, on the other hand, tend to suffer from rough approximations or a lack of completeness, limiting their predictive power. In this paper, we present a model that incorporates central concepts of astrophysics as well as reliable results from recent simulations of molecular clouds and their evolutionary paths. Based on that, we construct a self-consistent semi-analytical framework that describes the formation, evolution and star formation activity of molecular clouds, including a number of feedback effects to account for the complex processes inside those objects. The final equation system is solved numerically but at much lower computational expense than, e.g., hydrodynamical descriptions of comparable systems. The model presented in this paper agrees well with a broad range of observational results, showing that molecular cloud evolution can be understood as an interplay between accretion, global collapse, star formation and stellar feedback.

Molecular clouds are to a great extent influenced by turbulent motions in the gas. Numerical and observational studies indicate that the star formation rate and efficiency crucially depend on the mixture of solenoidal and compressive modes in the turbulent acceleration field, which can be quantified by the turbulent driving parameter b. For purely solenoidal (divergence-free) driving previous studies showed that b=1/3 and for entirely compressive (curl-free) driving b=1. In this study, we determine the evolution of the turbulent driving parameter b in magnetohydrodynamical simulations of molecular cloud formation and evolution. The clouds form due to the convergence of two flows of warm neutral gas. We explore different scenarios by varying the magnitude of the initial turbulent perturbations in the flows. We show that the driving mode of the turbulence within the cloud strongly fluctuates with time and exhibits no clear correlation with typical cloud properties, such as the cloud mass and the (Alfven) Mach number. We specifically find that $b$ strongly varies from b=0.3 to b=0.8 on timescales t<5 Myr, where the timescale and range of variation can change from cloud to cloud. This rapid change of b from solenoidal to compressive driving is primarily associated with global contraction of the cloud and subsequent onset of star formation. We conclude that the effective turbulence driving parameter should be treated as a free parameter that can vary from solenoidal to compressive in both time and space.

Aims: In this paper we present a case study to investigate conditions necessary to detect a characteristic magnetic field substructure embedded in a large-scale field. A helical magnetic field with a surrounding hourglass shaped field is expected from theoretical predictions and self-consistent magnetohydrodynamical (MHD) simulations to be present in the specific case of protostellar outflows. Hence, such an outflow environment is the perfect for our study. Methodes: We present synthetic polarisation maps in the infrared and millimeter regime of protostellar outflows performed with the newly developed RT and polarisation code POLARIS. The code, as the first, includes a self-consistent description of various alignement mechanism like the imperfect Davis-Greenstein (IDG) and the radiative torque (RAT) alignment. We investigate for which effects the grain size distribution, and applied alignement mechanism have. Results: We find that the IDG mechanism cannot produce any measurable polarization degree (< 1 %) whereas RAT alignment produced polarization degrees of a few 1 %. Furthermore, we developed a method to identify the origin of the polarization. We show that the helical magnetic field in the outflow can only be observed close to the outflow axis and at its tip, whereas in the surrounding regions the hourglass field in the foreground dominates the polarization. Furthermore, the polarization degree in the outflow lobe is lower than in the surroundings in agreement with observations. We also find that the orientation of the polarization vector flips around a few 100 micron due to the transition from dichroic extinction to thermal re-emission. Hence, in order to avoid ambiguities when interpreting polarization data, we suggest to observed in the far-infrared and mm regime. Finally, we show that with ALMA it is possible to observe the polarization emerging from protostellar outflows.

High-mass stars are expected to form from dense prestellar cores. Their precise formation conditions are widely discussed, including their virial condition, which results in slow collapse for super-virial cores with strong support by turbulence or magnetic fields, or fast collapse for sub-virial sources. To disentangle their formation processes, measurements of the deuterium fractions are frequently employed to approximately estimate the ages of these cores and to obtain constraints on their dynamical evolution. We here present 3D magneto-hydrodynamical simulations including for the first time an accurate non-equilibrium chemical network with 21 gas-phase species plus dust grains and 213 reactions. With this network we model the deuteration process in fully depleted prestellar cores in great detail and determine its response to variations in the initial conditions. We explore the dependence on the initial gas column density, the turbulent Mach number, the mass-to-magnetic flux ratio and the distribution of the magnetic field, as well as the initial ortho-to-para ratio of H2. We find excellent agreement with recent observations of deuterium fractions in quiescent sources. Our results show that deuteration is rather efficient, even when assuming a conservative ortho-to-para ratio of 3 and highly sub-virial initial conditions, leading to large deuterium fractions already within roughly a free-fall time. We discuss the implications of our results and give an outlook to relevant future investigations.

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

We present synthetic ALMA observations of Keplerian, protostellar discs in the Class 0 stage studying the emission of molecular tracers like $^{13}$CO, C$^{18}$O, HCO$^+$, H$^{13}$CO$^+$, N$_2$H$^+$, and H$_2$CO. We model the emission of discs around low- and intermediate-mass protostars. We show that under optimal observing conditions ALMA is able to detect the discs already in the earliest stage of protostellar evolution, although the emission is often concentrated to the innermost 50 AU. Therefore, a resolution of a few 0.1" might be too low to detect Keplerian discs around Class 0 objects. We also demonstrate that under optimal conditions for edge-on discs Keplerian rotation signatures are recognisable, from which protostellar masses can be inferred. For this we here introduce a new approach, which allows us to determine protostellar masses with higher fidelity than before. Furthermore, we show that it is possible to reveal Keplerian rotation even for strongly inclined discs and that ALMA should be able to detect possible signs of fragmentation in face-on discs. In order to give some guidance for future ALMA observations, we investigate the influence of varying observing conditions and source distances. We show that it is possible to probe Keplerian rotation in inclined discs with an observing time of 2 h and a resolution of 0.1", even in the case of moderate weather conditions. Furthermore, we demonstrate that under optimal conditions, Keplerian discs around intermediate-mass protostars should be detectable up to kpc-distances.

We present radiation hydrodynamic simulations of collapsing protostellar cores with initial masses of 30, 100, and 200 M$_{\odot}$. We follow their gravitational collapse and the formation of a massive protostar and protostellar accretion disk. We employ a new hybrid radiative feedback method blending raytracing techniques with flux-limited diffusion for a more accurate treatment of the temperature and radiative force. In each case, the disk that forms becomes Toomre-unstable and develops spiral arms. This occurs between 0.35 and 0.55 freefall times and is accompanied by an increase in the accretion rate by a factor of 2-10. Although the disk becomes unstable, no other stars are formed. In the case of our 100 and 200 M$_{\odot}$ simulation, the star becomes highly super-Eddington and begins to drive bipolar outflow cavities that expand outwards. These radiatively-driven bubbles appear stable, and appear to be channeling gas back onto the protostellar accretion disk. Accretion proceeds strongly through the disk. After 81.4 kyr of evolution, our 30 M$_{\odot}$ simulation shows a star with a mass of 5.48 M$_{\odot}$ and a disk of mass 3.3 M$_{\odot}$, while our 100 M$_{\odot}$ simulation forms a 28.8 M$_{\odot}$ mass star with a 15.8 M$_{\odot}$ disk over the course of 41.6 kyr, and our 200 M$_{\odot}$ simulation forms a 43.7 M$_{\odot}$ star with an 18 M$_{\odot}$ disk in 21.9 kyr. In the absence of magnetic fields or other forms of feedback, the masses of the stars in our simulation do not appear limited by their own luminosities.

We use magnetohydrodynamical simulations of converging warm neutral medium flows to analyse the formation and global evolution of magnetised and turbulent molecular clouds subject to supernova feedback from massive stars. We show that supernova feedback alone fails to disrupt entire, gravitationally bound, molecular clouds, but is able to disperse small--sized (~10 pc) regions on timescales of less than 1 Myr. Efficient radiative cooling of the supernova remnant as well as strong compression of the surrounding gas result in non-persistent energy and momentum input from the supernovae. However, if the time between subsequent supernovae is short and they are clustered, large hot bubbles form that disperse larger regions of the parental cloud. On longer timescales, supernova feedback increases the amount of gas with moderate temperatures (T~300-3000 K). Despite its inability to disrupt molecular clouds, supernova feedback leaves a strong imprint on the star formation process. We find an overall reduction of the star formation efficiency by a factor of 2 and of the star formation rate by roughly factors of 2-4.

We present and test chemical models for three-dimensional hydrodynamical simulations of galaxies. We explore the effect of changing key parameters such as metallicity, radiation and non-equilibrium versus equilibrium metal cooling approximations on the transition between the gas phases in the interstellar medium. The microphysics is modelled by employing the public chemistry package KROME and the chemical networks have been tested to work in a wide range of densities and temperatures. We describe a simple H/He network following the formation of H$_2$, and a more sophisticated network which includes metals. Photochemistry, thermal processes, and different prescriptions for the H$_2$ catalysis on dust are presented and tested within a one-zone framework. The resulting network is made publicly available on the KROME webpage. We find that employing an accurate treatment of the dust-related processes induces a faster HI--H$_2$ transition. In addition, we show when the equilibrium assumption for metal cooling holds, and how a non-equilibrium approach affects the thermal evolution of the gas and the HII--HI transition. These models can be employed in any hydrodynamical code via an interface to KROME and can be applied to different problems including isolated galaxies, cosmological simulations of galaxy formation and evolution, supernova explosions in molecular clouds, and the modelling of star-forming regions. The metal network can be used for a comparison with observational data of CII 158 $\mu$m emission both for high-redshift as well as for local galaxies.

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

We present a GPU accelerated CUDA-C implementation of the Barnes Hut (BH) tree code for calculating the gravitational potential on octree adaptive meshes. The tree code algorithm is implemented within the FLASH4 adaptive mesh refinement (AMR) code framework and therefore fully MPI parallel. We describe the algorithm and present test results that demonstrate its accuracy and performance in comparison to the algorithms available in the current FLASH4 version. We use a MacLaurin spheroid to test the accuracy of our new implementation and use spherical, collapsing cloud cores with effective AMR to carry out performance tests also in comparison with previous gravity solvers. Depending on the setup and the GPU/CPU ratio, we find a speedup for the gravity unit of at least a factor of 3 and up to 60 in comparison to the gravity solvers implemented in the FLASH4 code. We find an overall speedup factor for full simulations of at least factor 1.6 up to a factor of 10

We analyse simulations of turbulent, magnetised molecular cloud cores focussing on the formation of Class 0 stage protostellar discs and the physical conditions in their surroundings. We show that for a wide range of initial conditions Keplerian discs are formed in the Class 0 stage already. Furthermore, we show that the accretion of mass and angular momentum in the surroundings of protostellar discs occurs in a highly anisotropic manner, by means of a few narrow accretion channels. The magnetic field structure in the vicinity of the discs is highly disordered, revealing field reversals up to distances of 1000 AU. These findings demonstrate that as soon as even mild turbulent motions are included, the classical disc formation scenario of a coherently rotating environment and a well-ordered magnetic field breaks down.

Molecular clouds, the birthplaces of stars in galaxies, form dynamically from the diffuse atomic gas of the interstellar medium (ISM). The ISM is also threaded by magnetic fields which have a large impact on its dynamics. In particular, star forming regions must be magnetically supercrit- ical in order to accomodate gas clumps which can collapse under their own weight. Based on a parameter study of three dimensional magneto-hydrodyamical (MHD) simulations, we show that the long-standing problem of how such supercritical regions are generated is still an open issue.

We improve previous calculations of the CMB spectral distortions due to the decay of primordial magnetic fields. We focus our studies on causally generated magnetic fields at the electroweak and QCD phase transitions. We also consider the decay of helical magnetic fields. We show that the decay of non-helical magnetic fields generated at either the electroweak or QCD scale produce $\mu$ and $y$-type distortions below $10^{-8}$ which are probably not detectable by a future PIXIE-like experiment. We show that magnetic fields generated at the electroweak scale must have a helicity fraction $f_*>10^{-4}$ in order to produce detectable $\mu$-type distortions. Hence a positive detection coming from the decay of magnetic fields would rule out non-helical primordial magnetic fields and provide a lower bound on the magnetic helicity.

In this letter we show that magnetic fields generated at the electroweak phase transition are most likely too weak to explain the void magnetic fields apparently observed today unless they have considerable helicity. We show that, in the simplest estimates, the helicity naturally produced in conjunction with the baryon asymmetry is too small to explain observations, which require a helicity fraction at least of order $10^{-14}-10^{-10}$ depending on the void fields constraint used. Therefore new mechanisms to generate primordial helicity are required if magnetic fields generated during the electroweak phase transition should explain the extragalactic fields.

We use magnetohydrodynamical simulations of converging flows to investigate the process of molecular cloud formation and evolution out of the magnetised ISM. Here, we investigate whether the observed subcritical HI clouds can become supercritical and hence allow the formation of stars within them. To do so, we vary the turbulent Mach number of the flows, as well as the initial magnetic field strength. We show that dense cores are able to build up under all conditions, but that star formation in these cores is either heavily delayed or completely suppressed if the initial field strength is B>3 microGauss. To probe the effect of magnetic diffusion, we introduce a tilting angle between the flows and the uniform background magnetic field, which mimics non--ideal MHD effects. Even with highly inclined flows, the formed cores are devoid of star formation, because no magnetically supercritical regions are build up. Hence we conclude, that the problem of how supercritical cloud cores are generated still persists.

We present an algorithm for solving the radiative transfer problem on massively parallel computers using adaptive mesh refinement and domain decomposition. The solver is based on the method of characteristics which requires an adaptive raytracer that integrates the equation of radiative transfer. The radiation field is split into local and global components which are handled separately to overcome the non-locality problem. The solver is implemented in the framework of the magneto-hydrodynamics code FLASH and is coupled by an operator splitting step. The goal is the study of radiation in the context of star formation simulations with a focus on early disc formation and evolution. This requires a proper treatment of radiation physics that covers both the optically thin as well as the optically thick regimes and the transition region in particular. We successfully show the accuracy and feasibility of our method in a series of standard radiative transfer problems and two 3D collapse simulations resembling the early stages of protostar and disc formation.

We present an interface between the (magneto-) hydrodynamics code PLUTO and the plasma simulation and spectral synthesis code CLOUDY. By combining these codes, we constructed a new photoionization hydrodynamics solver: The PLUTO-CLOUDY Interface (TPCI), which is well suited to simulate photoevaporative flows under strong irradiation. The code includes the electromagnetic spectrum from X-rays to the radio range and solves the photoionization and chemical network of the 30 lightest elements. TPCI follows an iterative numerical scheme: First, the equilibrium state of the medium is solved for a given radiation field by CLOUDY, resulting in a net radiative heating or cooling. In the second step, the latter influences the (magneto-) hydrodynamic evolution calculated by PLUTO. Here, we validated the one-dimensional version of the code on the basis of four test problems: Photoevaporation of a cool hydrogen cloud, cooling of coronal plasma, formation of a Stroemgren sphere, and the evaporating atmosphere of a hot Jupiter. This combination of an equilibrium photoionization solver with a general MHD code provides an advanced simulation tool applicable to a variety of astrophysical problems.

To understand the evolution of planetary systems, it is important to investigate planets in highly evolved stellar systems, and to explore the implications of their observed properties with respect to potential formation scenarios. Observations suggest the presence of giant planets in post-common-envelope binaries (PCEBs). A particularly well-studied system with planetary masses of 1.7 M_J and 7.0 M_J is NN Ser. We show here that a pure first-generation scenario where the planets form before the common envelope (CE) phase and the orbits evolve due to the changes in the gravitational potential is inconsistent with the current data. We propose a second-generation scenario where the planets are formed from the material that is ejected during the CE, which may naturally explain the observed planetary masses. In addition, hybrid scenarios where the planets form before the CE and evolve due to the accretion of the ejected gas appear as a realistic possibility.

We analyse simulations of turbulent, magnetised molecular cloud cores focussing on the formation of Class 0 stage protostellar discs and the physical conditions in their surroundings. We show that for a wide range of initial conditions Keplerian discs are formed in the Class 0 stage already. In particular, we show that even subsonic turbulent motions reduce the magnetic braking efficiency sufficiently in order to allow rotationally supported discs to form. We therefore suggest that already during the Class 0 stage the fraction of Keplerian discs is significantly higher than 50%, consistent with recent observational trends but significantly higher than predictions based on simulations with misaligned magnetic fields, demonstrating the importance of turbulent motions for the formation of Keplerian discs. We show that the accretion of mass and angular momentum in the surroundings of protostellar discs occurs in a highly anisotropic manner, by means of a few narrow accretion channels. The magnetic field structure in the vicinity of the discs is highly disordered, revealing field reversals up to distances of 1000 AU. These findings demonstrate that as soon as even mild turbulent motions are included, the classical disc formation scenario of a coherently rotating environment and a well-ordered magnetic field breaks down. Hence, it is highly questionable to assess the magnetic braking efficiency based on non-turbulent collapse simulation. We strongly suggest that, in addition to the global magnetic field properties, the small-scale accretion flow and detailed magnetic field structure have to be considered in order to assess the likelihood of Keplerian discs to be present.

Radiation feedback plays a crucial role in the process of star formation. In order to simulate the thermodynamic evolution of disks, filaments, and the molecular gas surrounding clusters of young stars, we require an efficient and accurate method for solving the radiation transfer problem. We describe the implementation of a hybrid radiation transport scheme in the adaptive grid-based FLASH general magnetohydrodynamics code. The hybrid scheme splits the radiative transport problem into a raytracing step and a diffusion step. The raytracer captures the first absorption event, as stars irradiate their environments, while the evolution of the diffuse component of the radiation field is handled by a flux-limited diffusion (FLD) solver. We demonstrate the accuracy of our method through a variety of benchmark tests including the irradiation of a static disk, subcritical and supercritical radiative shocks, and thermal energy equilibration. We also demonstrate the capability of our method for casting shadows and calculating gas and dust temperatures in the presence of multiple stellar sources. Our method enables radiation-hydrodynamic studies of young stellar objects, protostellar disks, and clustered star formation in magnetized, filamentary environments.

Here we present a simple, but nevertheless, instructive model for the star formation efficiency in turbulent molecular clouds. The model is based on the assumption of log-normal density distribution which reflects the turbulent nature of the interstellar medium (ISM). Together with the number count of cloud cores, which follows a Salpeter-like core mass function (CMF), and the minimum mass for the collapse of individual cloud cores, given by the local Jeans mass, we are able to derive the SFE for clouds as a function of their Jeans masses. We find a very generic power-law, SFE \propto (M_cloud/M_J)^{-0.26} and a maximum SFE_max \sim 1/3 for the Salpeter case. This result is independent of the turbulent Mach number but fairly sensitive to variations of the CMF.

Powerful jets and outflows are launched from the protostellar disks around newborn stars. These outflows carry enough mass and momentum to transform the structure of their parent molecular cloud and to potentially control star formation itself. Despite their importance, we have not been able to fully quantify the impact of jets and outflows during the formation of a star cluster. The main problem lies in limited computing power. We would have to resolve the magnetic jet-launching mechanism close to the protostar and at the same time follow the evolution of a parsec-size cloud for a million years. Current computer power and codes fall orders of magnitude short of achieving this. In order to overcome this problem, we implement a subgrid-scale (SGS) model for launching jets and outflows, which demonstrably converges and reproduces the mass, linear and angular momentum transfer, and the speed of real jets, with ~ 1000 times lower resolution than would be required without SGS model. We apply the new SGS model to turbulent, magnetized star cluster formation and show that jets and outflows (1) eject about 1/4 of their parent molecular clump in high-speed jets, quickly reaching distances of more than a parsec, (2) reduce the star formation rate by about a factor of two, and (3) lead to the formation of ~ 1.5 times as many stars compared to the no-outflow case. Most importantly, we find that jets and outflows reduce the average star mass by a factor of ~ 3 and may thus be essential for understanding the characteristic mass of the stellar initial mass function.

In this paper we show that the Universe is already strongly magnetized at very early epochs during cosmic evolution. Our calculations are based on the efficient amplification of weak magnetic seed fields, which are unavoidably present in the early Universe, by the turbulent small-scale dynamo. We identify two mechanisms for the generation of turbulence in the radiation dominated epoch where velocity fluctuations are produced by the primordial density perturbation and by possible first-order phase transitions at the electroweak or QCD scales. We show that all the necessities for the small-scale dynamo to work are fulfilled. Hence, this mechanism, operating due to primordial density perturbations, guarantees fields with comoving field strength $B_0\sim10^{-6}\varepsilon^{1/2}$ nG on scales up to $\lambda_c\sim0.1$ pc, where $\varepsilon$ is the saturation efficiency. The amplification of magnetic seed fields could be even larger if there are first-order phase transitions in the early Universe. Where, on scales up to $\lambda_c\sim100$ pc, the comoving field strength due to this mechanism will be $B_0\sim10^{-3}\varepsilon^{1/2}$ nG at the present time. Such fields, albeit on small scales, can play an important role in structure formation and could provide an explanation to the apparently observed magnetic fields in the voids of the large-scale structure.

The first supernova explosions are potentially relevant sources for the production of the first large-scale magnetic fields. For this reason we present a set of high resolution simulations studying the effect of supernova explosions on magnetized, primordial halos. We focus on the evolution of an initially small-scale magnetic field formed during the collapse of the halo. We vary the degree of magnetization, the halo mass, and the amount of explosion energy in order to account for expected variations as well as to infer systematical dependencies of the results on initial conditions. Our simulations suggest that core collapse supernovae with an explosion energy of 10^51 erg and more violent pair instability supernovae with 10^53 erg are able to disrupt halos with masses up to a about 10^6 and 10^7 M_sun, respectively. The peak of the magnetic field spectra shows a continuous shift towards smaller k-values, i.e. larger length scales, over time reaching values as low as k = 4. On small scales the magnetic energy decreases at the cost of the energy on large scales resulting in a well-ordered magnetic field with a strength up to ~ 10^-8 G depending on the initial conditions. The coherence length of the magnetic field inferred from the spectra reaches values up to 250 pc in agreement with those obtained from autocorrelation functions. We find the coherence length to be as large as $50\%$ of the radius of the supernova bubble. Extrapolating this relation to later stages we suggest that significantly strong magnetic fields with coherence lengths as large as 1.5 kpc could be created. We discuss possible implications of our results on processes like recollapse of the halo, first galaxy formation, and the magnetization of the intergalactic medium.

(Abridged) In this review we focus on the observations and theory of the formation of early disks and outflows, and their connections with the first phases of planet formation. Large rotationally supported circumstellar disks, although common around more evolved young stellar objects, are rarely detected during the earliest, "Class 0" phase; however, a few excellent candidates have been discovered recently around both low and high mass protostars. In this early phase, prominent outflows are ubiquitously observed; they are expected to be associated with at least small magnetized disks. Disk formation - once thought to be a simple consequence of the conservation of angular momentum during hydrodynamic core collapse - is far more subtle in magnetized gas. In this case, the rotation can be strongly magnetically braked. Indeed, both analytic arguments and numerical simulations have shown that disk formation is suppressed in the strict ideal magnetohydrodynamic (MHD) limit for the observed level of core magnetization. We review what is known about this "magnetic braking catastrophe", possible ways to resolve it, and the current status of early disk observations. Outflows are also intimately linked to disk formation; they are a natural product of magnetic fields and rotation and are important signposts of star formation. We review new developments on early outflow generation since PPV. The properties of early disks and outflows are a key component of planet formation in its early stages and we review these major connections.

In this paper, we study the general impact of stellar mass-ejection events in planetary orbits in post-common envelope binaries with circumbinary planets like those around NN Serpentis. We discuss a set of simple equations that determine upper and lower limits for orbital expansion and investigate the effect of initial eccentricity. We deduce the range of possible semi-major axes and initial eccentricity values of the planets prior to the common-envelope event. In addition to spherically-symmetric mass-ejection events, we consider planetary dynamics under the influence of an expanding disk. In order to have survived, we suggest that the present planets in NN Ser must have had semi-major axes $\,{}^>_{\sim}\, 10$ AU and high eccentricity values which is in conflict with current observations. Consequently, we argue that these planets were not formed together with their hosting stellar system, but rather originated from the fraction of matter of the envelope that remained bound to the binary. According to the cooling age of the white dwarf primary of $10^6$ yr, the planets around NN Ser might be the youngest known so far and open up a wide range of further study of second generation planet formation.

The formation of stars is usually accompanied by the launching of protostellar outflows. Observations with the Atacama Large Millimetre/sub-millimetre Array (ALMA) will soon revolutionalise our understanding of the morphologies and kinematics of these objects. In this paper, we present synthetic ALMA observations of protostellar outflows based on numerical magnetohydrodynamic collapse simulations. We find significant velocity gradients in our outflow models and a very prominent helical structure within the outflows. We speculate that the disk wind found in the ALMA Science Verification Data of HD 163296 presents a first instance of such an observation.

We present collapse simulations of strongly magnetised, turbulent molecular cloud cores with masses ranging from 2.6 to 1000 M_sun in order to study the influence of the initial conditions on the turbulence-induced disc formation mechanism proposed recently by Seifried et al. 2012. We find that Keplerian discs are formed in all cases independently of the core mass, the strength of turbulence, or the presence of global rotation. The discs appear within a few kyr after the formation of the protostar, are 50 - 150 AU in size, and have masses between 0.05 and a few 0.1 M_sun. During the formation of the discs the mass-to-flux ratio stays well below the critical value of 10 for Keplerian disc formation. Hence, flux-loss alone cannot explain the formation of Keplerian discs. The formation of rotationally supported discs at such early phases is rather due to the disordered magnetic field structure and due to turbulent motions in the surroundings of the discs, two effects lowering the classical magnetic braking efficiency. Binary systems occurring in the discs are mainly formed via the disc capturing mechanism rather than via disc fragmentation, which is largely suppressed by the presence of magnetic fields.

We study the influence of the choice of cooling function on the formation of molecular clouds in high-resolution three-dimensional simulations of converging flows. We directly compare the results obtained using the simple, parametrized cooling function introduced by Koyama & Inutsuka (2002) and used by a number of converging flow studies with the results of the detailed calculation of the non-equilibrium chemistry and thermal balance of the gas. We find that a number of the cloud properties, such as the mass and volume filling fractions of cold gas, are relatively insensitive to the choice of cooling function. On the other hand, the cloud morphology and the large-scale velocity distribution of the gas do strongly depend on the cooling function. We show that the differences that we see can largely be explained by differences in the way that Lyman-alpha cooling is treated in the two complementary approaches, and that a proper non-equilibrium treatment of the ionisation and recombination of the gas is necessary in order to model the high-temperature cooling correctly. We also investigate the properties of the dense clumps formed within the cloud. In agreement with previous models, we find that the majority of these clumps are not self-gravitating, suggesting that some form of large-scale collapse of the cloud may be required in order to produce gravitationally unstable clumps and hence stars. Overall, the physical properties of the dense clumps are similar in both simulations, suggesting that they do not depend strongly on the choice of cooling function. However, we do find a systematic difference of around 10K in the mean temperatures of the clumps produced by the two models.

The Square Kilometre Array (SKA) is the most ambitious radio telescope ever planned. With a collecting area of about a square kilometre, the SKA will be far superior in sensitivity and observing speed to all current radio facilities. The scientific capability promised by the SKA and its technological challenges provide an ideal base for interdisciplinary research, technology transfer, and collaboration between universities, research centres and industry. The SKA in the radio regime and the European Extreme Large Telescope (E-ELT) in the optical band are on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI) and have been recognised as the essential facilities for European research in astronomy. This "White Paper" outlines the German science and R&D interests in the SKA project and will provide the basis for future funding applications to secure German involvement in the Square Kilometre Array.

We explore the amplification of magnetic fields in the high-redshift Universe. For this purpose, we perform high-resolution cosmological simulations following the formation of primordial halos with \sim10^7 M_solar, revealing the presence of turbulent structures and complex morphologies at resolutions of at least 32 cells per Jeans length. Employing a turbulence subgrid-scale model, we quantify the amount of unresolved turbulence and show that the resulting turbulent viscosity has a significant impact on the gas morphology, suppressing the formation of low-mass clumps. We further demonstrate that such turbulence implies the efficient amplification of magnetic fields via the small-scale dynamo. We discuss the properties of the dynamo in the kinematic and non-linear regime, and explore the resulting magnetic field amplification during primordial star formation. We show that field strengths of \sim10^{-5} G can be expected at number densities of \sim5 cm^{-3}.

The first galaxies form due to gravitational collapse of primordial halos. During this collapse, weak magnetic seed fields get amplified exponentially by the small-scale dynamo - a process converting kinetic energy from turbulence into magnetic energy. We use the Kazantsev theory, which describes the small-scale dynamo analytically, to study magnetic field amplification for different turbulent velocity correlation functions. For incompressible turbulence (Kolmogorov turbulence), we find that the growth rate is proportional to the square root of the hydrodynamic Reynolds number, Re^(1/2). In the case of highly compressible turbulence (Burgers turbulence) the growth rate increases proportional to Re^(1/3). With a detailed chemical network we are able to follow the chemical evolution and determine the kinetic and magnetic viscosities (due to Ohmic and ambipolar diffusion) during the collapse of the halo. This way, we can calculate the growth rate of the small-scale dynamo quantitatively and predict the evolution of the small-scale magnetic field. As the magnetic energy is transported to larger scales on the local eddy-timescale, we obtain an estimate for the magnetic field on the Jeans scale. Even there, we find that equipartition with the kinetic energy is reached on small timescales. Dynamically relevant field structures can thus be expected already during the formation of the first objects in the Universe.

Stars form by the gravitational collapse of interstellar gas. The thermodynamic response of the gas can be characterized by an effective equation of state. It determines how gas heats up or cools as it gets compressed, and hence plays a key role in regulating the process of stellar birth on virtually all scales, ranging from individual star clusters up to the galaxy as a whole. We present a systematic study of the impact of thermodynamics on gravitational collapse in the context of high-redshift star formation, but argue that our findings are also relevant for present-day star formation in molecular clouds. We consider a polytropic equation of state, P = k rho^Gamma, with both sub-isothermal exponents Gamma < 1 and super-isothermal exponents Gamma > 1. We find significant differences between these two cases. For Gamma > 1, pressure gradients slow down the contraction and lead to the formation of a virialized, turbulent core. Weak magnetic fields are strongly tangled and efficiently amplified via the small-scale turbulent dynamo on timescales corresponding to the eddy-turnover time at the viscous scale. For Gamma < 1, on the other hand, pressure support is not sufficient for the formation of such a core. Gravitational contraction proceeds much more rapidly and the flow develops very strong shocks, creating a network of intersecting sheets and extended filaments. The resulting magnetic field lines are very coherent and exhibit a considerable degree of order. Nevertheless, even under these conditions we still find exponential growth of the magnetic energy density in the kinematic regime.

The formation of massive stars exceeding 10 solar masses usually results in large-scale molecular outflows. Numerical simulations, including ionization, of the formation of such stars show evidence for ionization-driven molecular outflows. We here examine whether the outflows seen in these models reproduce the observations. We compute synthetic ALMA and CARMA maps of CO emission lines of the outflows, and compare their signatures to existing single-dish and interferometric data. We find that the ionization-driven models can only reproduce weak outflows around high-mass star-forming regions. We argue that expanding H II regions probably do not represent the dominant mechanism for driving observed outflows. We suggest instead that observed outflows are driven by the collective action of the outflows from the many lower-mass stars that inevitably form around young massive stars in a cluster.

We present collapse simulations of strongly magnetised, 100 M_sun, turbulent cloud cores. Around the protostars formed during the collapse Keplerian discs with typical sizes of up to 100 AU build up in contrast to previous simulations neglecting turbulence. Analysing the condensations in which the discs form, we show that the magnetic flux loss is not sufficient to explain the build-up of Keplerian discs. The average magnetic field is strongly inclined to the disc which might reduce the magnetic braking efficiency. However, the main reason for the reduced magnetic braking efficiency is the highly disordered magnetic field in the surroundings of the discs. Furthermore, due to the lack of a coherently rotating structure in the turbulent environment of the disc no toroidal magnetic field necessary for angular momentum extraction can build up. Simultaneously the angular momentum inflow remains high due to local shear flows created by the turbulent motions. We suggest that the "magnetic braking catastrophe" is an artefact of the idealised non-turbulent initial conditions and that turbulence provides a natural mechanism to circumvent this problem.

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.

We present collapse simulations of 100 M_{\sun}, turbulent cloud cores threaded by a strong magnetic field. During the initial collapse phase filaments are generated which fragment quickly and form several protostars. Around these protostars Keplerian discs with typical sizes of up to 100 AU build up in contrast to previous simulations neglecting turbulence. We examine three mechanisms potentially responsible for lowering the magnetic braking efficiency and therefore allowing for the formation of Keplerian discs. Analysing the condensations in which the discs form, we show that the build-up of Keplerian discs is neither caused by magnetic flux loss due to turbulent reconnection nor by the misalignment of the magnetic field and the angular momentum. It is rather a consequence of the turbulent surroundings of the disc which exhibit no coherent rotation structure while strong local shear flows carry large amounts of angular momentum. We suggest that the "magnetic braking catastrophe", i.e. the formation of sub-Keplerian discs only, is an artefact of the idealised non-turbulent initial conditions and that turbulence provides a natural mechanism to circumvent this problem.

(abridged) The nature of early Class 0/I protostellar discs is not clearly understood. Early protostellar discs are needed to drive molecular outflows and jets observed in star forming regions, but there has been some debate to how they form. From a theoretical perspective, the consequences of disc and outflow generation are crucial to understanding the very nature of how stars are assembled. We have performed 3D ideal magnetohydrodynamic (MHD) simulations of collapsing Bonnor-Ebert spheres, employing sink particles with a radius of 3.2 AU alongside an AMR grid and using a cooling function to model radiative cooling of the gas. This has allowed us to explore 2-8x10^4 yr further into the evolution of an early Class 0 disc-outflow system than previous simulations. Our outflow is precessing, kinked, turbulent, contains internal shocks and has a scale of 0.1 pc end-to-end. We form a rotationally dominated disc with a radius of 100 AU embedded inside a transient, unstable, flattened, rotating core extending out to 2000 AU. The larger flattened structure launches a low speed wind (v_r < 1.5 km/s) dominated by B_phi, while the inner disc launches a centrifugally driven jet dominated by B_p with speeds up to 20$ km/s. From the inner disk, the value of dM/dt_out/dM/dt_in ~ 0.1, wheras in the outer core dM/dt_out/dM/dt_in ~ 1.0. The inner disc becomes unstable to a warping instability due to the magnetic structure of the outflow and warps to 30 deg with respect to the z-axis by the end of the simulation. The envelope is cleared out and is less massive than the disc. We measure star formation efficiencies of eta_core=0.63 (and growing), higher than theoretical predictions. This indicates that outflows are not as efficient at expelling envelope mass as some current models estimate. We discuss the relevance of our disc misalignment concerning the formation of mis-aligned hot Jupiters.

Here we summarize our recent results of high-resolution computer simulations on the turbulent amplification of weak magnetic seed fields showing that such fields will be exponentially amplified also during the gravitational collapse reminiscent to the situation during primordial star formation. The exponential magnetic field amplification is driven by the turbulent small-scale dynamo that can be only observed in computer simulations if the turbulent motions in the central core are sufficiently resolved. We find that the Jeans length, which determines the central core region, has to be resolved by at least 30 grid cells to capture the dynamo activity. We conclude from our studies that strong magnetic fields will be unavoidably created already during the formation of the first stars in the Universe, potentially influencing their evolution and mass distribution.

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.

We study the influence of initial conditions on the magnetic field amplification during the collapse of a magnetised gas cloud. We focus on the dependence of the growth and saturation level of the dynamo generated field on the turbulent properties of the collapsing cloud. In particular, we explore the effect of varying the initial strength and injection scale of turbulence and the initial uniform rotation of the collapsing magnetised cloud. In order to follow the evolution of the magnetic field in both the kinematic and the nonlinear regime, we choose an initial field strength of $\simeq 1\,\mkG$ with the magnetic to kinetic energy ratio, $E_{\rm m}/E_{\rm k} \sim 10^{-4}$. Both gravitational compression and the small-scale dynamo initially amplify the magnetic field. Further into the evolution, the dynamo-generated magnetic field saturates but the total magnetic field continues to grow because of compression. The saturation of the small-scale dynamo is marked by a change in the slope of $B/\rho^{2/3}$ and by a shift in the peak of the magnetic energy spectrum from small scales to larger scales. For the range of initial Mach numbers explored in this study, the dynamo growth rate increases as the Mach number increases from $v_{\rm rms}/c_{\rm s}\sim 0.2$ to 0.4 and then starts decreasing from $v_{\rm rms}/c_{\rm s}\sim 1.0$. We obtain saturation values of $E_{\rm m}/E_{\rm k} = 0.2 - 0.3$ for these runs. Simulations with different initial injection scales of turbulence also show saturation at similar levels. For runs with different initial rotation of the cloud, the magnetic energy saturates at $E_{\rm m}/E_{\rm k}\sim 0.2 - 0.4$ of the equipartition value. (Abridged)

Numerical simulations of outflows formed during the collapse of 100 M_sun cloud cores are presented. We derive a generalised criterion from MHD wind theory to analyse the launching mechanism of these outflows. The criterion is successfully applied to the whole outflow structure and cases with sub-Keplerian disc rotation. It allows us to decide whether an outflow is driven centrifugally or by the toroidal magnetic pressure. We show that quantities such as the magnetic field line inclination or the ratio of the toroidal to poloidal magnetic field alone are insufficient to determine the driving mechanism of outflows. By performing 12 runs we are able to study the influence of the initial conditions on the properties of outflows around massive protostars in detail. Our simulations reveal a strong effect of the magnetic field strength on the outflow morphology. In runs with weak fields or high rotational energies, well collimated, fast jets are observed whereas for strong fields poorly collimated, low-velocity outflows are found. We show that the occurrence of a fast jet is coupled to the existence of a Keplerian protostellar disc. Despite the very different morphologies all outflows are launched from the discs by centrifugal acceleration with the toroidal magnetic field increasingly contributing to the gas acceleration further away from the discs. The poor collimation of the outflows in runs with strong fields is a consequence of the weak hoop stresses. This in turn is caused by the slow build-up of a toroidal field due to sub-Keplerian disc rotation. The mass ejection/accretion ratios scatter around a mean of 0.3 in accordance with observational and analytical results. We suggest an evolutionary scenario for the earliest stage of massive star formation in which initially poorly collimated outflows develop which progressively get better collimated due to the generation of fast jets.

We study the mass-to-flux ratio (M/\Phi) of clumps and cores in simulations of supersonic, magnetohydrodynamical turbulence for different initial magnetic field strengths. We investigate whether the (M/\Phi)-ratio of core and envelope, R = (M/\Phi)_{core}/(M/\Phi)_{envelope} can be used to distinguish between theories of ambipolar diffusion and turbulence-regulated star formation. We analyse R for different Lines-of-Sight (LoS) in various sub-cubes of our simulation box. We find that, 1) the average and median values of |R| for different times and initial magnetic field strengths are typically greater, but close to unity, 2) the average and median values of |R| saturate at average values of |R| ~ 1 for smaller magnetic fields, 3) values of |R| < 1 for small magnetic fields in the envelope are caused by field reversals when turbulence twists the field lines such that field components in different directions average out. Finally, we propose two mechanisms for generating values |R| ~< 1 for the weak and strong magnetic field limit in the context of a turbulent model. First, in the weak field limit, the small-scale turbulent dynamo leads to a significantly increased flux in the core and we find |R| ~< 1. Second, in the strong field limit, field reversals in the envelope also lead to values |R| ~< 1. These reversals are less likely to occur in the core region where the velocity field is more coherent and the internal velocity dispersion is typically subsonic.

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.

We investigate the impact of different initial conditions for the initial density profile and the initial turbulence on the formation process of protostellar clusters. We study the collapse of dense molecula r cloud cores with three-dimensional adaptive mesh refinement simulations. We focus our discussion on the distribution of the gas among the protostellar objects in the turbulent dynamical cores. Despite the large variations in the initial configurations and the resulting gas and cluster morphology we find that all stellar clusters follow a very similar gas accretion behaviour. Once secondary protostars begin to form, the central region of a cluster is efficiently shielded from further accretion. Hence, objects located close to the centre are starved of material, as indicated by a strong decrease of the central acc retion rate. This Fragmentation Induced Starvation occurs not only in rotationally supported discs and filaments, but also in more spherically symmetric clusters with complex chaotic motions.

We discuss the amplification of magnetic fields by the small-scale dynamo, a process that could efficiently produce strong magnetic fields in the first galaxies. In addition, we derive constraints on the primordial field strength from the epoch of reionization.

Understanding the origin of high-mass stars is central to modern astrophysics. We shed light on this problem with simulations using a novel, adaptive-mesh, ray-tracing algorithm. These simulations consistently follow the gravitational collapse of a massive molecular cloud core, the subsequent build-up and fragmentation of the accretion disk surrounding the nascent star, and, for the first time, the interaction between its intense UV radiation field and the infalling material. We show that ionization feedback can neither stop protostellar mass growth nor suppress fragmentation. We discuss the effects of feedback by ionizing and non-ionizing radiation on the evolution of the stellar cluster. The accretion is not limited by radiative feedback but by the formation of low-mass companions in a process we call "fragmentation-induced starvation". This behavior consistently reproduces the observed relation between the most massive star and the total mass of stars in a cluster. We show that magnetic fields reduce the star formation rate and lead to the formation of more massive stars.

We study the growth rate and saturation level of the turbulent dynamo in magnetohydrodynamical simulations of turbulence, driven with solenoidal (divergence-free) or compressive (curl-free) forcing. For models with Mach numbers ranging from 0.02 to 20, we find significantly different magnetic field geometries, amplification rates, and saturation levels, decreasing strongly at the transition from subsonic to supersonic flows, due to the development of shocks. Both extreme types of turbulent forcing drive the dynamo, but solenoidal forcing is more efficient, because it produces more vorticity.

We present simulations of collapsing 100 M_\sun mass cores in the context of massive star formation. The effect of variable initial rotational and magnetic energies on the formation of massive stars is studied in detail. We focus on accretion rates and on the question under which conditions massive Keplerian disks can form in the very early evolutionary stage of massive protostars. For this purpose, we perform 12 simulations with different initial conditions extending over a wide range in parameter space. The equations of magnetohydrodynamics (MHD) are solved under the assumption of ideal MHD. We find that the formation of Keplerian disks in the very early stages is suppressed for a mass-to-flux ratio normalised to the critical value \mu below 10, in agreement with a series of low-mass star formation simulations. This is caused by very efficient magnetic braking resulting in a nearly instantaneous removal of angular momentum from the disk. For weak magnetic fields, corresponding to \mu > 10, large-scale, centrifugally supported disks build up with radii exceeding 100 AU. A stability analysis reveals that the disks are supported against gravitationally induced perturbations by the magnetic field and tend to form single stars rather than multiple objects. We find protostellar accretion rates of the order of a few 10^-4 M_\sun yr^-1 which, considering the large range covered by the initial conditions, vary only by a factor of ~ 3 between the different simulations. We attribute this fact to two competing effects of magnetic fields. On the one hand, magnetic braking enhances accretion by removing angular momentum from the disk thus lowering the centrifugal support against gravity. On the other hand, the combined effect of magnetic pressure and magnetic tension counteracts gravity by exerting an outward directed force on the gas in the disk thus reducing the accretion onto the protostars.

Ultracompact and hypercompact HII regions appear when a star with a mass larger than about 15 solar masses starts to ionize its own environment. Recent observations of time variability in these objects are one of the pieces of evidence that suggest that at least some of them harbor stars that are still accreting from an infalling neutral accretion flow that becomes ionized in its innermost part. We present an analysis of the properties of the HII regions formed in the 3D radiation-hydrodynamic simulations presented by Peters et al. as a function of time. Flickering of the HII regions is a natural outcome of this model. The radio-continuum fluxes of the simulated HII regions, as well as their flux and size variations are in agreement with the available observations. From the simulations, we estimate that a small but non-negligible fraction (~ 10 %) of observed HII regions should have detectable flux variations (larger than 10 %) on timescales of ~ 10 years, with positive variations being more likely to happen than negative variations. A novel result of these simulations is that negative flux changes do happen, in contrast to the simple expectation of ever growing HII regions. We also explore the temporal correlations between properties that are directly observed (flux and size) and other quantities like density and ionization rates.

We employ simulations of supersonic super-Alfvenic turbulence decay as a benchmark test problem to assess and compare the performance of nine astrophysical MHD methods actively used to model star formation. The set of nine codes includes: ENZO, FLASH, KT-MHD, LL-MHD, PLUTO, PPML, RAMSES, STAGGER, and ZEUS. We present a comprehensive set of statistical measures designed to quantify the effects of numerical dissipation in these MHD solvers. We compare power spectra for basic fields to determine the effective spectral bandwidth of the methods and rank them based on their relative effective Reynolds numbers. We also compare numerical dissipation for solenoidal and dilatational velocity components to check for possible impacts of the numerics on small-scale density statistics. Finally, we discuss convergence of various characteristics for the turbulence decay test and impacts of various components of numerical schemes on the accuracy of solutions. We show that the best performing codes employ a consistently high order of accuracy for spatial reconstruction of the evolved fields, transverse gradient interpolation, conservation law update step, and Lorentz force computation. The best results are achieved with divergence-free evolution of the magnetic field using the constrained transport method, and using little to no explicit artificial viscosity. Codes which fall short in one or more of these areas are still useful, but they must compensate higher numerical dissipation with higher numerical resolution. This paper is the largest, most comprehensive MHD code comparison on an application-like test problem to date. We hope this work will help developers improve their numerical algorithms while helping users to make informed choices in picking optimal applications for their specific astrophysical problems.

We present an analysis of the large-scale molecular cloud structure and of the stability of clumpy structures in nearby molecular clouds. In our recent work, we identified a structural transition in molecular clouds by studying the probability distributions of gas column densities in them. In this paper, we further examine the nature of this transition. The transition takes place at the visual extinction of A_V^tail = 2-4 mag, or equivalently, at \Sigma^tail = 40-80 Ms pc^{-2}. The clumps identified above this limit have wide ranges of masses and sizes, but a remarkably constant mean volume density of n = 10^3 cm^{-3}. This is 5-10 times larger than the density of the medium surrounding the clumps. By examining the stability of the clumps, we show that they are gravitationally unbound entities, and that the external pressure from the parental molecular cloud is a significant source of confining pressure for them. Then, the structural transition at A_V^tail may be linked to a transition between this population and the surrounding medium. The star formation rates in the clouds correlate strongly with the total mass in the clumps, i.e, with the mass above A_V^tail, dropping abruptly below that threshold. These results imply that the formation of pressure confined clumps introduces a prerequisite for star formation. Furthermore, they give a physically motivated explanation for the recently reported relation between the star formation rates and the amount of dense material in molecular clouds. Likewise, they give rise to a natural threshold for star formation at A_V^tail.

We present a detailed parameter study of collapsing turbulent cloud cores, varying the initial density profile and the initial turbulent velocity field. We systematically investigate the influence of different initial conditions on the star formation process, mainly focusing on the fragmentation, the number of formed stars, and the resulting mass distributions. Our study compares four different density profiles (uniform, Bonnor-Ebert type, $\rho\propto r^{-1.5}$, and $\rho\propto r^{-2}$), combined with six different supersonic turbulent velocity fields (compressive, mixed, and solenoidal, initialised with two different random seeds each) in three-dimensional simulations using the adaptive-mesh refinement, hydrodynamics code FLASH. The simulations show that density profiles with flat cores produce hundreds of low-mass stars, either distributed throughout the entire cloud or found in subclusters, depending on the initial turbulence. Concentrated density profiles always lead to the formation of one high-mass star in the centre of the cloud and, if at all, low-mass stars surrounding the central one. In uniform and Bonnor-Ebert type density distributions, compressive initial turbulence leads to local collapse about 25% earlier than solenoidal turbulence. However, central collapse in the steep power-law profiles is too fast for the turbulence to have any significant influence. We conclude that (I) the initial density profile and turbulence mainly determine the cloud evolution and the formation of clusters, (II) the initial mass function (IMF) is not universal for all setups, and (III) that massive stars are much less likely to form in flat density distributions. The IMFs obtained in the uniform and Bonnor-Ebert type density profiles are more consistent with the observed IMF, but shifted to lower masses.

Recent FERMI observations provide a lower limit of 10^{-15} G for the magnetic field strength in the intergalactic medium (IGM). This is consistent with theoretical expectations based on the Biermann battery effect, which predicts such IGM fields already at redshifts z~10. During gravitational collapse, such magnetic fields can be amplified by compression and by turbulence, giving rise to the small-scale dynamo. On scales below the Jeans length, the eddy turnover timescale is much shorter than the free-fall timescale, so that saturation can be reached during collapse. This scenario has been tested and confirmed with magneto-hydrodynamical simulations following the collapse of a turbulent, weakly magnetized cloud. Based on a spectral analysis, we confirm that turbulence is injected on the Jeans scale. For the power spectrum of the magnetic field, we obtain the Kazantsev slope which is characteristic for the small-scale dynamo. A calculation of the critical length scales for ambipolar diffusion and Ohmic dissipation shows that these scales are always small enough to allow significant amplification of the magnetic field by small-scale eddies. We discuss potential implications for the protostellar accretion disk, with particular focus on the magneto-rotational instability, which may change the morphology of the disk and reduce the accretion rate by a factor of a few.

Cosmic structure formation is characterized by the complex interplay between gravity, turbulence, and magnetic fields. The processes by which gravitational energy is converted into turbulent and magnetic energies, however, remain poorly understood. Here, we show with high-resolution, adaptive-mesh simulations that MHD turbulence is efficiently driven by extracting energy from the gravitational potential during the collapse of a dense gas cloud. Compressible motions generated during the contraction are converted into solenoidal, turbulent motions, leading to a natural energy ratio of E_sol/E_tot of approximately 2/3. We find that the energy injection scale of gravity-driven turbulence is close to the local Jeans scale. If small seeds of the magnetic field are present, they are amplified exponentially fast via the small-scale dynamo process. The magnetic field grows most efficiently on the smallest scales, for which the stretching, twisting, and folding of field lines, and the turbulent vortices are sufficiently resolved. We find that this scale corresponds to about 30 grid cells in the simulations. We thus suggest a new minimum resolution criterion of 30 cells per Jeans length in (magneto)hydrodynamical simulations of self-gravitating gas, in order to resolve turbulence on the Jeans scale, and to capture minimum dynamo amplification of the magnetic field. Due to numerical diffusion, however, any existing simulation today can at best provide lower limits on the physical growth rates. We conclude that a small, initial magnetic field can grow to dynamically important strength on time scales significantly shorter than the free-fall time of the cloud.

Massive stars disproportionately influence their surroundings. How they form has only started to become clear recently through radiation gas dynamical simulations. However, until now, no simulation has simultaneously included both magnetic fields and ionizing radiation. Here we present the results from the first radiation-magnetohydrodynamical (RMHD) simulation including ionization feedback, comparing an RMHD model of a 1000 M_sol rotating cloud to earlier radiation gas dynamical models with the same initial density and velocity distributions. We find that despite starting with a strongly supercritical mass to flux ratio, the magnetic field has three effects. First, the field offers locally support against gravitational collapse in the accretion flow, substantially reducing the amount of secondary fragmentation in comparison to the gas dynamical case. Second, the field drains angular momentum from the collapsing gas, further increasing the amount of material available for accretion by the central, massive, protostar, and thus increasing its final mass by about 50% from the purely gas dynamical case. Third, the field is wound up by the rotation of the flow, driving a tower flow. However, this flow never achieves the strength seen in low-mass star formation simulations for two reasons: gravitational fragmentation disrupts the circular flow in the central regions where the protostars form, and the expanding H II regions tend to further disrupt the field geometry. Therefore, outflows driven by ionization heating look likely to be more dynamically important in regions of massive star formation.

We investigate the formation and evolution of giant molecular clouds (GMCs) by the collision of convergent warm neutral medium (WNM) streams in the interstellar medium, in the presence of magnetic fields and ambipolar diffusion (AD), focusing on the evolution of the star formation rate (SFR) and efficiency (SFE), as well as of the mass-to-magnetic-flux ratio (M2FR) in the forming clouds. We find that: 1) Clouds formed by supercritical inflow streams proceed directly to collapse, while clouds formed by subcritical streams first contract and then re-expand, oscillating on the scale of tens of Myr. 2) Our suite of simulations with initial magnetic field strength of 2, 3, and 4 $\mu\G$ show that only supercritical or marginal critical streams lead to reasonable star forming rates. 3) The GMC's M2FR is a generally increasing function of time, whose growth rate depends on the details of how mass is added to the GMC from the WNM. 4) The M2FR is a highly fluctuating function of position in the clouds. 5) In our simulations, the SFE approaches stationarity, because mass is added to the GMC at a similar rate at which it converts mass to stars. In such an approximately stationary regime, the SFE provides a proxy of the supercritical mass fraction in the cloud. 6) We observe the occurrence of buoyancy of the low-M2FR regions within the gravitationally-contracting GMCs, so that the latter naturally segregate into a high-density, high-M2FR "core" and a low-density, low-M2FR "envelope", without the intervention of AD. (Abridged)

Massive stars influence their surroundings through radiation, winds, and supernova explosions far out of proportion to their small numbers. However, the physical processes that initiate and govern the birth of massive stars remain poorly understood. Two widely discussed models are monolithic collapse of molecular cloud cores and competitive accretion. To learn more about massive star formation, we perform simulations of the collapse of rotating, massive, cloud cores including radiative heating by both non-ionizing and ionizing radiation using the FLASH adaptive mesh refinement code. These simulations show fragmentation from gravitational instability in the enormously dense accretion flows required to build up massive stars. Secondary stars form rapidly in these flows and accrete mass that would have otherwise been consumed by the massive star in the center, in a process that we term fragmentation-induced starvation. This explains why massive stars are usually found as members of high-order stellar systems that themselves belong to large clusters containing stars of all masses. The radiative heating does not prevent fragmentation, but does lead to a higher Jeans mass, resulting in fewer and more massive stars than would form without the heating. This mechanism reproduces the observed relation between the total stellar mass in the cluster and the mass of the largest star. It predicts strong clumping and filamentary structure in the center of collapsing cores, as has recently been observed. We speculate that a similar mechanism will act during primordial star formation.

We implemented sink particles in the Adaptive Mesh Refinement (AMR) code FLASH to model the gravitational collapse and accretion in turbulent molecular clouds and cores. Sink particles are frequently used to measure properties of star formation in numerical simulations, such as the star formation rate and efficiency, and the mass distribution of stars. We show that only using a density threshold for sink particle creation is insufficient in case of supersonic flows, because the density can exceed the threshold in strong shocks that do not necessarily lead to local collapse. Additional physical collapse indicators have to be considered. We apply our AMR sink particle module to the formation of a star cluster, and compare it to a Smoothed Particle Hydrodynamics (SPH) code with sink particles. Our comparison shows encouraging agreement of gas and sink particle properties between the AMR and SPH code.

We explore the amplification of magnetic seed fields during the formation of the first stars and galaxies. During gravitational collapse, turbulence is created from accretion shocks, which may act to amplify weak magnetic fields in the protostellar cloud. Numerical simulations showed that such turbulence is sub-sonic in the first star-forming minihalos, and highly supersonic in the first galaxies with virial temperatures larger than 10^4 K. We investigate the magnetic field amplification during the collapse both for Kolmogorov and Burgers-type turbulence with a semi-analytic model that incorporates the effects of gravitational compression and small-scale dynamo amplification. We find that the magnetic field may be substantially amplified before the formation of a disk. On scales of 1/10 of the Jeans length, saturation occurs after ~10^8 yr. Although the saturation behaviour of the small-scale dynamo is still somewhat uncertain, we expect a saturation field strength of the order ~10^{-7} n^{0.5} G in the first star-forming halos, with n the number density in cgs units. In the first galaxies with higher turbulent velocities, the magnetic field strength may be increased by an order of magnitude, and saturation may occur after 10^6 to 10^7 yr. In the Kolmogorov case, the magnetic field strength on the integral scale (i.e. the scale with most magnetic power) is higher due to the characteristic power-law indices, but the difference is less than a factor of 2 in the saturated phase. Our results thus indicate that the precise scaling of the turbulent velocity with length scale is of minor importance. They further imply that magnetic fields will be significantly enhanced before the formation of a protostellar disk, where they may change the fragmentation properties of the gas and the accretion rate.

Cosmological hydrodynamical simulations of primordial star formation suggest that the gas within the first star-forming halos is turbulent. This has strong implications on the subsequent evolution, in particular on the generation of magnetic fields. Using high-resolution numerical simulations, we show that in the presence of turbulence, weak seed magnetic fields are exponentially amplified by the small-scale dynamo during the formation of the first stars. We conclude that strong magnetic fields are generated during the birth of the first stars in the universe, potentially modifying the mass distribution of these stars and influencing the subsequent cosmic evolution. We find that the presence of the small-scale turbulent dynamo can only be identified in numerical simulations in which the turbulent motions in the central core are resolved with at least 32 grid cells.

Massive stars influence the surrounding universe far out of proportion to their numbers through ionizing radiation, supernova explosions, and heavy element production. Their formation requires the collapse of massive interstellar gas clouds with very high accretion rates. We discuss results from the first three-dimensional simulations of the gravitational collapse of a massive, rotating molecular cloud core that include heating by both non-ionizing and ionizing radiation. Local gravitational instabilities in the accretion flow lead to the build-up of a small cluster of stars. These lower-mass companions subsequently compete with the high-mass star for the same common gas reservoir and limit its overall mass growth. This process is called fragmentation-induced starvation, and explains why massive stars are usually found as members of high-order stellar systems. These simulations also show that the HII regions forming around massive stars are initially trapped by the infalling gas, but soon begin to fluctuate rapidly. Over time, the same ultracompact HII region can expand anisotropically, contract again, and take on any of the observed morphological classes. The total lifetime of HII regions is given by the global accretion timescale, rather than their short internal sound-crossing time. This solves the so-called lifetime problem of ultracompact HII region. We conclude that the the most significant differences between the formation of low-mass and high-mass stars are all explained as the result of rapid accretion within a dense, gravitationally unstable flow.

The spatial morphology, spectral characteristics, and time variability of ultracompact H II regions provide strong constraints on the process of massive star formation. We have performed simulations of the gravitational collapse of rotating molecular cloud cores, including treatments of the propagation of ionizing and non-ionizing radiation. We here present synthetic radio continuum observations of H II regions from our collapse simulations, to investigate how well they agree with observation, and what we can learn about how massive star formation proceeds. We find that intermittent shielding by dense filaments in the gravitationally unstable accretion flow around the massive star leads to highly variable H II regions that do not grow monotonically, but rather flicker, growing and shrinking repeatedly. This behavior appears able to resolve the well-known lifetime problem. We find that multiple ionizing sources generally form, resulting in groups of ultracompact H II regions, consistent with observations. We confirm that our model reproduces the qualitative H II region morphologies found in surveys, with generally consistent relative frequencies. We also find that simulated spectral energy distributions (SEDs) from our model are consistent with the range of observed H II region SEDs, including both regions showing a normal transition from optically thick to optically thin emission, and those with intermediate spectral slopes. In our models, anomalous slopes are solely produced by inhomogeneities in the H II region, with no contribution from dust emission at millimeter or submillimeter wavelengths. We conclude that many observed characteristics of ultracompact H II regions appear consistent with massive star formation in fast, gravitationally unstable, accretion flows.

We discuss the formation of supermassive black holes in the early universe, and how to probe their subsequent evolution with the upcoming mm/sub-mm telescope ALMA. We first focus on the chemical and radiative conditions for black hole formation, in particular considering radiation trapping and molecular dissociation effects. We then turn our attention towards the magnetic properties in the halos where the first black holes form, and show that the presence of turbulence may lead to a magnetic dynamo, which could support the black hole formation process by providing an efficient means of transporting the angular momentum. We finally focus on observable properties of high-redshift black holes with respect to ALMA, and discuss how to distinguish between chemistry driven by the starburst and chemistry driven by X-rays from the black hole.

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.

We describe the first three-dimensional simulation of the gravitational collapse of a massive, rotating molecular cloud that includes heating by both non-ionizing and ionizing radiation. We find that as the first protostars gain sufficient mass to ionize the accretion flow, their H II regions are initially gravitationally trapped, but soon begin to rapidly fluctuate between trapped and extended states, in agreement with observations. Over time, the same ultracompact H II region can expand anisotropically, contract again, and take on any of the observed morphological classes. In their extended phases, expanding H II regions drive bipolar neutral outflows characteristic of high-mass star formation. The total lifetime of H II regions is given by the global accretion timescale, rather than their short internal sound-crossing time. The pressure of the hot, ionized gas does not terminate accretion. Instead the final stellar mass is set by fragmentation-induced starvation. Local gravitational instabilities in the accretion flow lead to the build-up of a small cluster of stars, all with relatively high masses due to heating from accretion radiation. These companions subsequently compete with the initial high-mass star for the same common gas reservoir and limit its mass growth. Our findings show that the most significant differences between the formation of low-mass and high-mass stars are all explained as the result of rapid accretion within a dense, gravitationally unstable, ionized flow.

We explore the effects of magnetic energy dissipation on the formation of the first stars. For this purpose, we follow the evolution of primordial chemistry in the presence of magnetic fields in the post-recombination universe until the formation of the first virialized halos. From the point of virialization, we follow the protostellar collapse up to densities of $\sim10^{12}$ cm$^{-3}$ in a one-zone model. In the intergalactic medium (IGM), comoving field strengths of $\gtrsim0.1$ nG lead to Jeans masses of $10^8 M_\odot$ or more and thus delay gravitational collapse in the first halos until they are sufficiently massive. During protostellar collapse, we find that the temperature minimum at densities of $\sim10^3$ cm$^{-3}$ does not change significantly, such that the characteristic mass scale for fragmentation is not affected. However, we find a significant temperature increase at higher densities for comoving field strengths of $\gtrsim0.1$ nG. This may delay gravitational collapse, in particular at densities of $\sim10^9$ cm$^{-3}$, where the proton abundance drops rapidly and the main contribution to the ambipolar diffusion resistivity is due to collisions with Li$^+$. After the formation of the protostar, the increased gas temperatures may enhance the protostellar accretion rate. Our model confirms that initial weak magnetic fields may be amplified considerably during gravitational collapse and become dynamically relevant. For instance, a comoving field strength above $10^{-5}$ nG will be amplified above the critical value for the onset of jets which can magnetize the IGM.

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 ...

Dark stars powered by dark matter annihilation have been proposed as the first luminous sources in the universe. These stars are believed to form in the central dark matter cusp of low-mass minihalos. Recent calculations indicate stellar masses up to \sim1000 solar masses and/or have very long lifetimes. The UV photons from these objects could therefore contribute significantly to cosmic reionization. Here we show that such dark star models would require a somewhat artificial reionization history, based on a double-reionization phase and a late star-burst near redshift $z\sim6$, in order to fulfill the WMAP constraint on the optical depth as well as the Gunn-Peterson constraint at $z\sim6$. This suggests that, if dark stars were common in the early universe, then models are preferred which predict a number of UV photons similar to conventional Pop. III stars. This excludes dark stars with 100 solar masses that enter a main-sequence phase and other models that lead to a strong increase in the number of UV photons. We also derive constraints for massive as well as light dark matter candidates from the observed X-ray, gamma-ray and neutrino background, considering dark matter profiles which have been steepened during the formation of dark stars. This increases the clumping factor at high redshift and gives rise to a higher dark matter annihilation rate in the early universe. We furthermore estimate the potential contribution from the annihilation products in the remnants of dark stars, which may provide a promising path to constrain such models further, but which is currently still uncertain.

Light dark matter annihilating into electron-positron pairs emits a significant amount of internal bremsstrahlung that may contribute to the cosmic gamma-ray background. The amount of emitted gamma-rays depends on the dark matter clumping factor. Recent calculations indicate that this value should be of order $10^6-10^7$. That allows us to calculate the expected gamma-ray background contribution from dark matter annihilation. We find that the light dark matter model can be ruled out if a constant thermally-averaged cross section is assumed (s-wave annihilation). For more massive dark matter candidates like neutralinos, however, cosmic constraints are weaker.

Magnetic fields in the early universe can significantly alter the thermal evolution and the ionization history during the dark ages. This is reflected in the 21 cm line of atomic hydrogen, which is coupled to the gas temperature through collisions at high redshifts, and through the Wouthuysen-Field effect at low redshifts. We present a semi-analytic model for star formation and the build-up of a Lyman alpha background in the presence of magnetic fields, and calculate the evolution of the mean 21 cm brightness temperature and its frequency gradient as a function of redshift. We further discuss the evolution of linear fluctuations in temperature and ionization in the presence of magnetic fields and calculate the effect on the 21 cm power spectrum. At high redshifts, the signal is increased compared to the non-magnetic case due to the additional heat input into the IGM from ambipolar diffusion and the decay of MHD turbulence. At lower redshifts, the formation of luminous objects and the build-up of a Lyman alpha background can be delayed by a redshift interval of 10 due to the strong increase of the filtering mass scale in the presence of magnetic fields. This tends to decrease the 21 cm signal compared to the zero-field case. In summary, we find that 21 cm observations may become a promising tool to constrain primordial magnetic fields.

We calculate the reionization history for different models of the stellar population and explore the effects of primordial magnetic fields, dark matter decay and dark matter annihilation on reionization. We find that stellar populations based on a Scalo-type initial mass function for Population II stars can be ruled out as sole sources for reionization, unless star formation efficiencies of more than 10% or very high photon escape fractions from the parental halo are adopted. When considering primordial magnetic fields, we find that the additional heat injection from ambipolar diffusion and decaying MHD turbulence has significant impact on the thermal evolution and the ionization history of the post-recombination universe and on structure formation. The magnetic Jeans mass changes the typical mass scale of the star forming halos, and depending on the adopted stellar model we derive upper limits to the magnetic field strength between 0.7 and $5 $nG (comoving). For dark matter annihilation, we find an upper limit to the thermally averaged mass-weighted cross section of $10^{-33} \mathrm{cm}^3\mathrm{/s/eV}$. For dark matter decay, our calculations yield a lower limit to the lifetime of dark matter particles of $3\times10^{23}$ s. These limits are in agreement with constraints from recombination and the X-ray background and provide an independent confirmation at a much later epoch.

Astrophysical jets are associated with the formation of young stars of all masses, stellar and massive black holes, and perhaps even with the formation of massive planets. Their role in the formation of planets, stars, and galaxies is increasingly appreciated and probably reflects a deep connection between the accretion flows - by which stars and black holes may be formed - and the efficiency by which magnetic torques can remove angular momentum from such flows. We compare the properties and physics of jets in both non-relativistic and relativistic systems and trace, by means of theoretical argument and numerical simulations, the physical connections between these different phenomena. We discuss the properties of jets from young stars and black holes, give some basic theoretical results that underpin the origin of jets in these systems, and then show results of recent simulations on jet production in collapsing star-forming cores as well as from jets around rotating Kerr black holes.

We present 3D radiation-gasdynamical simulations of an ionization front running into a dense clump. In our setup, a B0 star irradiates an overdensity which is at a distance of 10 pc and modelled as a supercritical 100 M_sol Bonnor-Ebert sphere. The radiation from the star heats up the gas and creates a shock front that expands into the interstellar medium. The shock compresses the clump material while the ionizing radiation heats it up. The outcome of this "cloud-crushing" process is a fully turbulent gas in the wake of the clump. In the end, the clump entirely dissolves. We propose that this mechanism is very efficient in creating short-living supersonic turbulence in the vicinity of massive stars.

We present results from our numerical simulations of collapsing massive molecular cloud cores. These numerical calculations show that massive stars assemble quickly with mass accretion rates exceeding 10^-3 Msol/yr and confirm that the mass accretion during the collapsing phase is much more efficient than predicted by selfsimilar collapse solutions, dM/dt ~ c^3/G. We find that during protostellar assembly out of a non-turbulent core, the mass accretion reaches 20 - 100 c^3/G. Furthermore, we explore the self-consistent structure of bipolar outflows that are produced in our three dimensional magnetized collapse simulations. These outflows produce cavities out of which radiation pressure can be released, thereby reducing the limitations on the final mass of massive stars formed by gravitational collapse. Additional enhancement of the mass accretion rate comes from accretion along filaments that are built up by supersonic turbulent motions. Our numerical calculations of collapsing turbulent cores result in mass accretion rates as high as 10^-2 Msol/yr.

{It has recently been proposed that giant molecular complexes form at the sites where streams of diffuse warm atomic gas collide at transonic velocities.} {We study the global statistics of molecular clouds formed by large scale colliding flows of warm neutral atomic interstellar gas under ideal MHD conditions. The flows deliver material as well as kinetic energy and trigger thermal instability leading eventually to gravitational collapse.} {We perform adaptive mesh refinement MHD simulations which, for the first time in this context, treat self-consistently cooling and self-gravity.} {The clouds formed in the simulations develop a highly inhomogeneous density and temperature structure, with cold dense filaments and clumps condensing from converging flows of warm atomic gas. In the clouds, the column density probability density distribution (PDF) peaks at $\sim 2 \times 10^{21} \psc$ and decays rapidly at higher values; the magnetic intensity correlates weakly with density from $n \sim 0.1$ to $10^4 \pcc$, and then varies roughly as $n^{1/2}$ for higher densities.} {The global statistical properties of such molecular clouds are reasonably consistent with observational determinations. Our numerical simulations suggest that molecular clouds formed by the moderately supersonic collision of warm atomic gas streams.}

Jets and outflows from young stellar objects are proposed candidates to drive supersonic turbulence in molecular clouds. Here, we present the results from multi-dimensional jet simulations where we investigate in detail the energy and momentum deposition from jets into their surrounding environment and quantify the character of the excited turbulence with velocity probability density functions. Our study include jet--clump interaction, transient jets, and magnetised jets. We find that collimated supersonic jets do not excite supersonic motions far from the vicinity of the jet. Supersonic fluctuations are damped quickly and do not spread into the parent cloud. Instead subsonic, non-compressional modes occupy most of the excited volume. This is a generic feature which can not be fully circumvented by overdense jets or magnetic fields. Nevertheless, jets are able to leave strong imprints in their cloud structure and can disrupt dense clumps. Our results question the ability of collimated jets to sustain supersonic turbulence in molecular clouds.

We present an investigation of massive star formation that results from the gravitational collapse of massive, magnetized molecular cloud cores. We investigate this by means of highly resolved, numerical simulations of initial magnetized Bonnor-Ebert-Spheres that undergo collapse and cooling. By comparing three different cases - an isothermal collapse, a collapse with radiative cooling, and a magnetized collapse - we show that massive stars assemble quickly with mass accretion rates exceeding 10^-3 Msol/yr. We confirm that the mass accretion during the collapsing phase is much more efficient than predicted by selfsimilar collapse solutions, i.e. dM/dt ~ c^3/G. We find that during protostellar assembly the mass accretion reaches 20 - 100 c^3/G. Furthermore, we determined the self-consistent structure of bipolar outflows that are produced in our three dimensional magnetized collapse simulations. These outflows produce cavities out of which radiation pressure can be released, thereby reducing the limitations on the final mass of massive stars formed by gravitational collapse. Moreover, we argue that the extraction of angular momentum by disk-threaded magnetic fields and/or by the appearance of bars with spiral arms significantly enhance the mass accretion rate, thereby helping the massive protostar to assemble more quickly.

We present a detailed computational study of the assembly of protostellar disks and massive stars in molecular clouds with supersonic turbulence. We follow the evolution of large scale filamentary structures in a cluster-forming clump down to protostellar length scales by means of very highly resolved, 3D adaptive mesh refined (AMR) simulations, and show how accretion disks and massive stars form in such environments. We find that an initially elongated cloud core which has a slight spin from oblique shocks collapses first to a filament and later develops a turbulent disk close to the center of the filament. The continued large scale flow that shocks with the filament maintains the high density and pressure within it. Material within the cooling filament undergoes gravitational collapse and an outside-in assembly of a massive protostar. Our simulations show that very high mass accretion rates of up to 10^-2 Msol/yr and high, supersonic, infall velocities result from such filamentary accretion. Accretion at these rates is higher by an order of magnitude than those found in semi-analytic studies, and can quench the radiation field of a growing massive young star.Our simulations include a comprehensive set of the important chemical and radiative processes such as cooling by molecular line emission, gas-dust interaction, and radiative diffusion in the optical thick regime, as well as H2 formation and dissociation. Therefore, we are able to probe, for the first time, the relevant physical phenomena on all scales from those characterizing the clump down to protostellar core.

In this work we investigate the effects of the primordial magnetic fields on cosmic microwave background anisotropies (CMB). Based on cosmological magnetohydrodynamic (MHD) simulations we calculate the CMB anisotropy spectra and polarization induced by fluid fluctuations (Alfv\'en modes) generated by primordial magnetic fields. The strongest effect on the CMB spectra comes from the transition epoch from a turbulent regime to a viscous regime. The balance between magnetic and kinetic energy until the onset of the viscous regime provides a one to one relation between the comoving coherence length $L$ and the comoving magnetic field strength $B$, such as $L \sim 30 (B/10^{-9}{\rm G})^3 \rm pc$. The resulting CMB temperature and polarization anisotropies are somewhat different from the ones previously obtained by using linear perturbation theory. Our calculation gives a constraint on the magnetic field strength in the intermediate scale of CMB observations. Upper limits are set by WMAP and BOOMERANG results for comoving magnetic field strength of $B < 28 \rm nG$ with a comoving coherence length of $L > 0.7 \rm Mpc$ for the most extreme case, or $B < 30 \rm nG$ and $L> 0.8 \rm Mpc$ for the most conservative case.

Star formation is usually accompanied by outflow phenomena. There is strong evidence that these outflows and jets are launched from the protostellar disk by magneto-rotational processes. Here, we report on our three dimensional, adaptive mesh, magneto-hydrodynamic simulations of collapsing, rotating, magnetized Bonnor-Ebert-Spheres whose properties are taken directly from observations. In contrast to the pure hydro case where no outflows are seen, our present simulations show an outflow from the protodisk surface at ~ AU and a jet at ~ 0.07 AU after a strong toroidal magnetic field build up. The large scale outflow, which extends up to ~ AU at the end of our simulation, is driven by toroidal magnetic pressure (spring), whereas the jet is powered by magneto-centrifugal force (fling). At the final stage of our simulation these winds are still confined within two respective shock fronts. Furthermore, we find that the jet-wind and the disk-anchored magnetic field extracts a considerable amount of angular momentum from the protostellar disk. The initial spin of our cloud core was chosen high enough to produce a binary system. We indeed find a close binary system (separation ~ 3 R_sol) which results from the fragmentation of an earlier formed ring structure. The magnetic field strength in these protostars reaches ~ 3 kG and becomes about 3 G at 1 AU from the center in agreement with recent observational results.

A large body of theoretical and computational work shows that jets - modelled as magnetized disk winds - exert an external torque on their underlying disks that can efficiently remove angular momentum and act as major drivers of disk accretion. These predictions have recently been confirmed in direct HST measurements of the jet rotation and angular momentum transport in low mass protostellar systems. We review the theory of disc winds and show that their physics is universal and scales to jets from both low and high mass star forming regions. This explains the observed properties of outflows in massive star forming regions, before the central massive star generates an ultracompact HII region. We also discuss the recent numerical studies on the formation of massive accretion disks and outflows through gravitational collapse, including our own work on 3D Adaptive Mesh simulations (using the FLASH code) of the hydromagnetic collapse of an initial rotating, and cooling Bonner-Ebert sphere. Magnetized collapse gives rise to outflows. Our own simulations show that both a jet-like disk wind on sub AU scales, and a larger scale molecular outflow occur (Banerjee and Pudritz 2005).

(abridged) A detailed examination of the evolution of stochastic magnetic fields between high cosmic temperatures and the present epoch is presented. A simple analytical model matching the results of the 3D MHD simulations allows for the prediction of present day magnetic field correlation lengths and energy. Our conclusions are multi fold. (a) Initial primordial fields with only a small amount of helicity are evolving into maximally helical fields. (b) There exists a correlation between the strength of the magnetic field, B, at the peak of it's spectrum and the location of the peak, given at the present epoch by: B ~ 5x10^{-12} (L/kpc) Gauss, where L is the correlation length determined by the initial magnetic field. (c) Concerning studies of generation of cosmic microwave background (CMBR) anisotropies due to primordial magnetic fields of B~10^{-9} Gauss on ~ 10 Mpc scales, such fields are not only impossible to generate in early causal magnetogenesis scenarios but also seemingly ruled out by distortions of the CMBR spectrum due to magnetic field dissipation on smaller scales and the overproduction of cluster magnetic fields. (d) The most promising detection possibility of CMBR distortions due to primordial magnetic fields may be on much smaller scales at higher multipoles l~10^6 where the signal is predicted to be the strongest. (e) It seems possible that magnetic fields in clusters of galaxies are entirely of primordial origin, without invoking dynamo amplification. Such fields would be of (pre-collapse) strength 10^{-12} - 10^{-11} Gauss with correlation lengths in the kpc range, and would also exist in voids of galaxies.

We present results of a detailed and fully non-linear numerical and analytical investigation of magnetic field evolution from the very earliest cosmic epochs to the present. We find that, under reasonable assumptions concerning the efficiency of a putative magnetogenesis era during cosmic phase transitions, surprisingly strong magnetic fields 10^{-13} - 10^{-11} Gauss, on comparatively small scales 100 pc - 10 kpc may survive to the present. Building on prior work on the evolution of magnetic fields during the course of gravitational collapse of a cluster, which indicates that pre-collapse fields of 4\times 10^{-12} Gauss extant on small scales may suffice to produce clusters with acceptable Faraday rotation measures, we question the widely hold view that cluster magnetic fields may not be entirely of primordial origin.

We present a detailed study of the collapse of molecular cloud cores using high resolution 3D adaptive mesh refinement (AMR) numerical simulations. In this first in a series of investigations our initial conditions consists of spherical molecular core obeying the hydrostatic Bonnor-Ebert-Profile with varying degrees of initial rotation. Our simulations cover both the formation of massive disks in which massive stars form as well as low mass disks. We use a customized version of the FLASH code whose AMR technique allows us to follow the formation of a protstellar disk and protostellar core(s) through more than ten orders in density increase while continuously resolving the local Jeans length (i.e. obeying the Truelove criterion, Truelove et al. (1997)). Our numerical simulations also incorporate the energy loss due to molecular line emission in order to obtain a more realistic picture of protostellar core and disk formation. Out initial states model system of mass 168 M_sol and 2.1 M_sol that will form high and low mass stars, respectively. We follow many features such as the development complex shock structures, and the possible fragmentation of the disk. We find that slowly rotating cores (Omega t_ff = 0.1) produce disks in which a strong bar develops but which does not fragment. Faster initial rotation rates (Omega t_ff = 0.2) result in the formation of a ring which may fragment into two star-forming cores. The size of the rings found in our simulated disks agree with the observations of similar systems.

It has been recently suggested that stable, supersymmetric B-balls formed in the early universe could not only be the dark matter at the present epoch, but also be responsible for baryogenesis by their partial evaporation at high temperatures. We reinvestigate the efficiency of B-ball baryogenesis and find it to be limited by the diffusion of baryon number away from the B-balls. Successful baryogenesis may only occur for B-balls with charges Q\simle 10^{20} - 5\times 10^{23}, which is close to the observational lower limits on the Q of a significant B-ball dark matter component. We also present some cosmological constraints on the abundances of larger B-balls in the early universe.