Deviations from the Schmidt-Kennicutt relations during early galaxy evolution

We utilize detailed time-varying models of the coupled evolution of stars and the HI, H_2, and CO-bright H_2 gas phases in galaxy-sized numerical simulations to explore the evolution of gas-rich and/or metal-poor systems, expected to be numerous in the Early Universe. The inclusion of the CO-bright H_2 gas phase, and the realistic rendering of star formation as an H_2-regulated process (and the new feedback processes that this entails) allows the most realistic tracking of strongly evolving galaxies, and much better comparison with observations. We find that while galaxies eventually settle into states conforming to Schmidt-Kennicutt (S-K) relations, significant and systematic deviations of their star formation rates (SFRs) from the latter occur, especially pronounced and prolonged for ... ...This indicates potentially serious limitations of (S-K)-type relations as reliable sub-grid elements of star formation physics in simulations of structure formation in the Early Universe. We anticipate that galaxies with marked deviations from the S-K relations will be found at high redshifts as unbiased inventories of total gas mass become possible with ALMA and the EVLA.Comment: 13 pages, 3 figures, accepted for publication in the Astrophysical Journa

Diagnostics of the molecular component of PDRs with mechanical heating. II: line intensities and ratios

CO observations in active galactic nuclei and star-bursts reveal high kinetic temperatures. Those environments are thought to be very turbulent due to dynamic phenomena such as outflows and high supernova rates. We investigate the effect of mechanical heating (MH) on atomic fine-structure and molecular lines, and their ratios. We use those ratios as a diagnostic to constrain the amount of MH in an object and also study its significance on estimating the H2 mass. Equilibrium PDRs models were used to compute the thermal and chemical balance for the clouds. The equilibria were solved for numerically using the optimized version of the Leiden PDR-XDR code. Large velocity gradient calculations were done as post-processing on the output of the PDR models using RADEX. High-J CO line ratios are very sensitive to MH. Emission becomes at least one order of magnitude brighter in clouds with n~10^5~cm^-3 and a star formation rate of 1 Solar Mass per year (corresponding to a MH rate of 2 * 10^-19 erg cm^-3 s^-1). Emission of low-J CO lines is not as sensitive to MH, but they do become brighter in response to MH. Generally, for all of the lines we considered, MH increases excitation temperatures and decreases the optical depth at the line centre. Hence line ratios are also affected, strongly in some cases. Ratios involving HCN are a good diagnostic for MH, such as HCN(1-0)/CO(1-0) and HCN(1-0)/HCO^+(1-0). Both ratios increase by a factor 3 or more for a MH equivalent to > 5 percent of the surface heating, as opposed to pure PDRs. The first major conclusion is that low-J to high-J intensity ratios will yield a good estimate of the MH rate (as opposed to only low-J ratios). The second one is that the MH rate should be taken into account when determining A_V or equivalently N_H, and consequently the cloud mass. Ignoring MH will also lead to large errors in density and radiation field estimates.Comment: 38 pages, to appear in A&

Face-on accretion onto a protoplanetary disc

Globular clusters (GCs) are known to harbor multiple stellar populations. To explain these observations Bastian et al. suggested a scenario in which a second population is formed by the accretion of enriched material onto the low-mass stars in the initial GC population. The idea is that the low-mass, pre-main sequence stars sweep up gas expelled by the massive stars of the same generation into their protoplanetary disc as they move through the GC core. We perform simulations with 2 different smoothed particle hydrodynamics codes to investigate if a low-mass star surrounded by a protoplanetary disc can accrete the amount of enriched material required in this scenario. We focus on the gas loading rate onto the disc and star as well as on the lifetime of the disc. We find that the gas loading rate is a factor of 2 smaller than the geometric rate, because the effective cross section of the disc is smaller than its surface area. The loading rate is consistent for both codes, irrespective of resolution. The disc gains mass in the high resolution runs, but loses angular momentum on a time scale of 10^4 yrs. Two effects determine the loss of (specific) angular momentum in our simulations: 1) continuous ram pressure stripping and 2) accretion of material with no azimuthal angular momentum. Our study and previous work suggest that the former, dominant process is mainly caused by numerical rather than physical effects, while the latter is not. The latter process causes the disc to become more compact, increasing the surface density profile at smaller radii. The disc size is determined in the first place by the ram pressure when the flow first hits the disc. Further evolution is governed by the decrease in the specific angular momentum of the disc. We conclude that the size and lifetime of the disc are probably not sufficient to accrete the amount of mass required in Bastian et al.'s scenario.Comment: Accepted for publication in A&A, 15 pages, 5 figures, 4 table

Molecular gas, CO, and star formation in galaxies: emergent empirical relations, feedback, and the evolution of very gas-rich systems

We use time-varying models of the coupled evolution of the HI, H_2 gas phases and stars in galaxy-sized numerical simulations to: a) test for the emergence of the Kennicutt-Schmidt (K-S) and the H_2-pressure relation, b) explore a realistic H_2-regulated star formation recipe which brings forth a neglected and potentially significant SF-regulating factor, and c) go beyond typical galactic environments (for which these galactic empirical relations are deduced) to explore the early evolution of very gas-rich galaxies. In this work we model low mass galaxies (M_{\rm baryon} \le 10^9 \msun), while incorporating an independent treatment of CO formation and destruction, the most important tracer molecule of H2 in galaxies, along with that for the H2 gas itself. We find that both the K-S and the H_2-pressure empirical relations can robustly emerge in galaxies after a dynamic equilibrium sets in between the various ISM states, the stellar component and its feedback. (abridged)Comment: 32 pages, 9 figures, accepted for publication in Ap

Various facets of intermolecular transfer of phase coherence by nuclear dipolar fields

It has long been recognized that dipolar fields can mediate intermolecular transfer of phase coherence from abundant solvent to sparse solute spins. Generally, the dipolar field has been considered while acting during prolonged free-precession delays. Recently, we have shown that transfer can also occur during suitable uninterrupted radio frequency pulse trains that are used for total correlation spectroscopy. Here, we will expand upon the latter work. First, analytical expressions for the evolution of the solvent magnetization under continuous irradiation and the influence of the dipolar field are derived. These expressions facilitate the simulations of the transfer process. Then, a pulse sequence for the retrieval of high-resolution spectra in inhomogeneous magnetic fields is described, along with another sequence to detect a transfer from an intermolecular double-quantum coherence. Finally, various schemes are discussed where the magnetization is modulated by a combination of multiple selective radio frequency pulses and pulsed field gradients in different directions. In these schemes, the magnetization is manipulated in such a way that the dipolar field, which originates from a single-spin species, can be decomposed into two components. Each component originates from a part of the magnetization that is modulated in a different direction. Both can independently, but simultaneously, mediate an intermolecular transfer of phase coherence.</p

Periodic bursts of Star Formation in Irregular Galaxies

We present N-body/SPH simulations of the evolution of an isolated dwarf galaxy including a detailed model for the ISM, star formation and stellar feedback. Depending on the strength of the feedback, the modelled dwarf galaxy shows periodic or quasi-periodic bursts of star formation of moderate strength. The period of the variations is related to the dynamical timescale, of the order of $1.5~10^8$ yr. We show that the results of these simulations are in good agreement with recent detailed observations of dwarf irregulars (dIrr) and that the peculiar kinematic and morphological properties of these objects,as revealed by high resolution HI studies, are fully reproduced. We discuss these results in the context of recent surveys of dwarf galaxies and point out that if the star formation pattern of our model galaxy is typical for dwarf irregulars this could explain the scatter of observed properties of dwarf galaxies. Specifically, we show that the time sampled distribution of the ratio between the instanteneous star formation rate (SFR) and the mean SFR is similar to that distribution in observed sample of dwarf galaxies.Comment: 11 pages, 6 figures, accepted for A&

Diagnostics of the Molecular Component of PDRs with Mechanical Heating

Context. Multitransition CO observations of galaxy centers have revealed that significant fractions of the dense circumnuclear gas have high kinetic temperatures, which are hard to explain by pure photon excitation, but may be caused by dissipation of turbulent energy. Aims. We aim to determine to what extent mechanical heating should be taken into account while modelling PDRs. To this end, the effect of dissipated turbulence on the thermal and chemical properties of PDRs is explored. Methods. Clouds are modelled as 1D semi-infinite slabs whose thermal and chemical equilibrium is solved for using the Leiden PDR-XDR code. Results. In a steady-state treatment, mechanical heating seems to play an important role in determining the kinetic temperature of the gas in molecular clouds. Particularly in high-energy environments such as starburst galaxies and galaxy centers, model gas temperatures are underestimated by at least a factor of two if mechanical heating is ignored. The models also show that CO, HCN and H2 O column densities increase as a function of mechanical heating. The HNC/HCN integrated column density ratio shows a decrease by a factor of at least two in high density regions with n \sim 105 cm-3, whereas that of HCN/HCO+ shows a strong dependence on mechanical heating for this same density range, with boosts of up to three orders of magnitude. Conclusions. The effects of mechanical heating cannot be ignored in studies of the molecular gas excitation whenever the ratio of the star formation rate to the gas density is close to, or exceeds, 7 \times 10-6 M yr-1 cm4.5 . If mechanical heating is not included, predicted column densities are underestimated, sometimes even by a few orders of magnitude. As a lower bound to its importance, we determined that it has non-negligible effects already when mechanical heating is as little as 1% of the UV heating in a PDR.Comment: 26 pages, 14 figures in the text and 13 figures as supplementary material. Accepted for publication in A&

Performance optimization and load-balancing modeling for superparametrization by 3D LES

In order to eliminate climate uncertainty w.r.t. cloud and convection parametrizations, superpramaterization (SP) [1] has emerged as one of the possible ways forward. We have implemented (regional) superparametrization of the ECMWF weather model OpenIFS [2] by cloud-resolving, three-dimensional large-eddy simulations. This setup, described in [3], contains a two-way coupling between a global meteorological model that resolves large-scale dynamics, with many local instances of the Dutch Atmospheric Large Eddy Simulation (DALES) [4], resolving cloud and boundary layer physics. The model is currently prohibitively expensive to run over climate or even seasonal time scales, and a global SP requires the allocation of millions of cores. In this paper, we study the performance and scaling behavior of the LES models and the coupling code and present our implemented optimizations. We mimic the observed load imbalance with a simple performance model and present strategies to improve hardware utilization in order to assess the feasibility of a world-covering superparametrization. We conclude that (quasi-)dynamical load-balancing can significantly reduce the runtime for such large-scale systems with wide variability in LES time-stepping speeds