Search

Current star formation models are based on the local structure of the interstellar medium (ISM), yet the details on how the small-scale physics propagates up to global galactic-scale properties are still under debate. To investigate this we use {\small VINTERGATAN}, a high-resolution (20 pc) cosmological zoom-in simulation of a Milky Way-like galaxy. We study how the velocity dispersion and density structure of the ISM on 50-100 pc scales evolve with redshift, and quantify their impact on the star formation efficiency per free-fall timescale, $\epsilon_{\rm ff}$. During starbursts the ISM can reach velocity dispersions as high as $\sim 50$ km s$^{-1}$ for the densest and coldest gas, most noticeable during the last major merger event ($1.3 < z < 1.5$). After a merger-dominated phase ($11$ to 0.3\% at $z<1$. The mass-weighted average shows a universal $\langle \epsilon_{\rm ff} \rangle \approx 1\%$, caused by an almost invariant virial parameter distribution in star forming clouds. Changes in their density and turbulence levels are coupled so the kinetic-to-gravitational energy ratio remains close to constant. Finally, we show that a \textit{theoretically} motivated instantaneous $\epsilon_{\rm ff}$ is intrinsically different to its \textit{observational} estimates adopting tracers of star formation e.g. H$\alpha$. Since the physics underlying star formation can be lost on short ($\sim$ 10 Myr) timescales, caution must be taken when constraining star formation models from observational estimates of $\epsilon_{\rm ff}$., Comment: 14 pages, 6 figures, comments welcome

We analyze the evolution of the radial profiles and the azimuthal variations of the stellar metallicities from the Vintergatan simulation of a Milky Way-like galaxy. We find that negative gradients exist as soon as the disk settles at high redshift, and are maintained throughout the long term evolution of the galaxy, including during major merger events. The inside-out growth of the disk and an overall outward radial migration tend to flatten these gradients in time. Major merger events only have a moderate and short-lived imprint on the [Fe/H] distributions with almost no radial dependence. The reason lies in the timescale for enrichment in Fe being significantly longer than the duration of the starbursts episodes, themselves slower than dynamical mixing during typical interactions. It results that signatures of major mergers become undetectable in [Fe/H] only a few Myr after pericenter passages. We note that considering other tracers like the warm interstellar medium, or monitoring the evolution of the metallicity gradient as a single value instead of a radial full profile could lead to different interpretations, and warn against an oversimplification of this complex problem., Comment: submitted to A&A

The clumpy nature of gas-rich galaxies at cosmic noon raises the question of universality of the scaling relations and average properties of the star-forming structures. Using controlled simulations of disk galaxies and varying only the gas fraction, we show that the influence of the galactic environments (large-scale turbulence, tides, shear) contributes, together with the different regime of instabilities, to setting a diversity of physical conditions for the formation and evolution of gas clumps from low to high gas fractions. However, the distributions of gas clumps at all gas fractions follow similar scaling relations as Larson's, suggesting the universality of median properties. Yet, we find that the scatter around these relations significantly increases with the gas fraction, allowing for the presence of massive, large, and highly turbulent clouds in gas-rich disks in addition to a more classical population of clouds. Clumps with an excess of mass for their size are slightly denser, more centrally concentrated, and host more abundant and faster star formation. We find that the star formation activity (rate, efficiency, depletion time) correlates much more strongly with the excess of mass than with the mass itself. Our results suggest the existence of universal scaling relations for gas clumps but with redshift-dependent scatters, which calls for deeper and more complete census of the populations of star-forming clumps and young stellar clusters at cosmic noon and beyond., Comment: A&A in press

The specific angular momenta ($j\equiv J/M$) of stars ($j_{\star}$), gas ($j_{\mathrm{gas}}$), baryons as a whole ($j_{\mathrm{b}}$) and dark matter haloes ($j_{\mathrm{h}}$) contain clues of vital importance about how galaxies form and evolve. Using one of the largest samples of disc galaxies (S0-BCD) with high-quality rotation curves and near-infrared surface photometry, we perform a detailed comparative analysis of $j$ that stretches across a variety of galaxy properties. Our analysis imposes tight constraints on the "retained" fractions of specific angular momentum ($j_{\star}/j_{\mathrm{h}}$, $j_{\mathrm{HI}}/j_{\mathrm{h}}$ and $j_{\mathrm{b}}/j_{\mathrm{h}}$), as well as on their systematic trends with mass fraction and galaxy morphology, thus on how well specific angular momentum is conserved in the process of disc galaxy formation and evolution. In particular, one of the most innovative results of our analysis is the finding that galaxies with larger baryon fractions have also retained larger fractions of their specific angular momentum. Furthermore, our analysis demonstrates how challenging it is to characterize barred galaxies from a gravitational instability point of view. This is true not only for the popular Efstathiou, Lake & Negroponte bar instability criterion, which fails to separate barred from non-barred galaxies in about 55% of the cases, but also for the mass-weighted Toomre parameter of atomic gas, $\langle Q_{\mathrm{HI}}\rangle$, which succeeds in separating barred from non-barred galaxies, but only in a statistical sense., Comment: MNRAS, in press

Supersonic gas turbulence is a ubiquitous property of the interstellar medium. The level of turbulence, quantified by the gas velocity dispersion ($\sigma_{\rm g}$), is observed to increase with the star formation rate (SFR) of a galaxy, but it is yet not established whether this trend is driven by stellar feedback or gravitational instabilities. In this work we carry out hydrodynamical simulations of entire disc galaxies, with different gas fractions, to understand the origins of the SFR-$\sigma_{\rm g}$ relation. We show that disc galaxies reach the same levels of turbulence regardless of the presence of stellar feedback processes, and argue that this is an outcome of the way disc galaxies regulate their gravitational stability. The simulations match the SFR-$\sigma_{\rm g}$ relation up to SFRs of the order of tens of M$_\odot$ yr$^{-1}$ and $\sigma_{\rm g}\sim 50$ km s$^{-1}$ in neutral hydrogen and molecular gas, but fail to reach the very large values ($>100$ km s$^{-1}$) reported in the literature for rapidly star-forming galaxies. We demonstrate that such high values of $\sigma_{\rm g}$ can be explained by 1) insufficient beam smearing corrections in observations, and 2) stellar feedback being coupled to the ionised gas phase traced by recombination lines. Given that the observed SFR-$\sigma_{\rm g}$ relation is composed of highly heterogeneous data, with $\sigma_{\rm g}$ at high SFRs almost exclusively being derived from H$\alpha$ observations of high redshift galaxies with complex morphologies, we caution against analytical models that attempt to explain the SFR-$\sigma_{\rm g}$ relation without accounting for these effects., Comment: Accepted for publication in MNRAS. 18 pages, 14 figures

The morphology of gas-rich disc galaxies at redshift ~1-3 is dominated by a few massive clumps. The process of formation or assembly of these clumps and their relation to molecular clouds in contemporary spiral galaxies are still unknown. Using simulations of isolated disc galaxies, we study how the structure of the interstellar medium and the stability regime of the discs change when varying the gas fraction. In all galaxies, the stellar component is the main driver of instabilities. However, the molecular gas plays a non-negligible role in the inter-clumps medium of gas-rich cases, and thus in the assembly of the massive clumps. At scales smaller than a few 100 pc, the Toomre-like disc instabilities are replaced by another regime, specially in the gas-rich galaxies. We find that galaxies at low gas fraction (10%) stand apart from discs with more gas, which all share similar properties on virtually all the aspects we explore. For gas fractions below approximately 20%, the clump-scale regime of instabilities disappears, only leaving the large-scale disc-driven regime. When translating the change of gas fraction to the cosmic evolution of galaxies, this transition marks the end of the clumpy phase of disc galaxies, and allows for the onset of spiral structures, as commonly found in the local Universe., Comment: 19 pages, MNRAS accepted

The stellar-to-halo mass relation (SHMR) is not only one of the main sources of information we have on the connection between galaxies and their dark matter haloes, but also an important indicator of the performance of galaxy formation models. Here we use one of the largest sample of galaxies with both high-quality rotation curves and near-infrared surface photometry, and perform a detailed comparative analysis of the SHMR. Our analysis shows that there are significant statistical differences between popular forms of the SHMR, and illustrates the predictive power of a new physically motivated scaling relation, which connects the stellar mass fraction ($M_{\star}/M_{\mathrm{h}}$) to the stellar specific angular momentum ($j_{\star}$) and the stellar radial velocity dispersion ($\sigma_{\star}$) via disc gravitational instability. Making use of such a relation, we demonstrate (i) how challenging it is to reproduce the efficiency of galaxy formation even for state-of-the-art cosmological hydrodynamical simulations, and (ii) that the evolution of the stellar mass fraction is regulated by disc gravitational instability: when $M_{\star}/M_{\mathrm{h}}$ varies, $j_{\star}$ and $\sigma_{\star}$ also vary as predicted by our scaling relation, thus erasing the memory of such evolution. This implies that the process of disc gravitational instability is intriguingly uniform across disc galaxies of all morphological types: from lenticulars to blue compact dwarfs. In particular, the cosmic variance of Toomre's $Q$ is 0.2 dex, a universal value for both stars and atomic gas., Comment: MNRAS, in press. Fig. 4, Sect. 4 and the last itemized paragraph of Sect. 5 are new, and demonstrate the importance of our results for the simulation community

We present a new set of galaxy scaling relations for the relative mass content of atomic gas, molecular gas and stars. Such relations are driven by disc gravitational instability, and originate from the low galaxy-to-galaxy variance of Toomre's $Q$ stability parameter. We test such relations using more than 100 galaxies, from massive spirals to dwarf irregulars, thus spanning several orders of magnitude in stellar mass ($M_{\star}\approx10^{6\mbox{-}11}\,\mbox{M}_{\odot}$) and atomic gas mass ($M_{\mathrm{HI}}\approx10^{7\mbox{-}10.5}\,\mbox{M}_{\odot}$). Such tests demonstrate (i) that our scaling relations are physically motivated and tightly constrained, (ii) that the mass-averaged gravitational instability properties of galaxy discs are remarkably uniform across the sequence Sa-dIrr, and (iii) that specific angular momentum plays an important role in such a scenario. Besides providing new insights into a very important topic in galaxy evolution, this work provides a simple formula (Eq. 5) that one can use for generating other galaxy relations driven by disc instability. We explain how to do that, mention a few possible applications, and stress the importance of testing our approach further., Comment: MNRAS, in press. Changed from letter to paper format, and revised significantly

Observations find a median star formation efficiency per free-fall time in Milky Way Giant Molecular Clouds (GMCs) on the order of $\epsilon_{\rm ff}\sim 1\%$ with dispersions of $\sim0.5\,{\rm dex}$. The origin of this scatter in $\epsilon_{\rm ff}$ is still debated and difficult to reproduce with analytical models. We track the formation, evolution and destruction of GMCs in a hydrodynamical simulation of a Milky Way-like galaxy and by deriving cloud properties in an observationally motivated way, measure the distribution of star formation efficiencies which are in excellent agreement with observations. We find no significant link between $\epsilon_{\rm ff}$ and any measured global property of GMCs (e.g. gas mass, velocity dispersion). Instead, a wide range of efficiencies exist in the entire parameter space. From the cloud evolutionary tracks, we find that each cloud follow a \emph{unique} evolutionary path which gives rise to wide diversity in all properties. We argue that it is this diversity in cloud properties, above all else, that results in the dispersion of $\epsilon_{\rm ff}$., Comment: 10 pages, 6 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

We introduce a new diagnostic for exploring the link between angular momentum and local gravitational instability in galaxy discs. Our diagnostic incorporates the latest developments in disc instability research, is fully consistent with approximations that are widely used for measuring the stellar specific angular momentum, $j_{\star}=J_{\star}/M_{\star}$, and is also very simple. We show that such a disc instability diagnostic hardly correlates with $j_{\star}$ or $M_{\star}$, and is remarkably constant across spiral galaxies of any given type (Sa$\!-\!$Sd), stellar mass ($M_{\star}=10^{9.5}\!-\!10^{11.5}\,\mathrm{M}_{\odot}$) and velocity dispersion anisotropy ($\sigma_{z\star}/\sigma_{R\star}=0\!-\!1$). The fact that $M_{\star}$ is tightly correlated with star formation rate ($\mathrm{SFR}$), molecular gas mass ($M_{\mathrm{mol}}$), metallicity ($12+\log\mathrm{O/H}$) and other fundamental galaxy properties thus implies that nearby star-forming spirals self-regulate to a quasi-universal disc stability level. This proves the existence of the self-regulation process postulated by several star formation models, but also raises important caveats., Comment: MNRAS Letters, in press. Minor revision to match the accepted version (added Fig. 1, Sect. 3.2, the final paragraph of Sect. 4, references and clarifications)

The stellar velocity dispersion, $\sigma$, is a quantity of crucial importance for spiral galaxies, where it enters fundamental dynamical processes such as gravitational instability and disc heating. Here we analyse a sample of 34 nearby spirals from the Calar Alto Legacy Integral Field Area (CALIFA) spectroscopic survey, deproject the line-of-sight $\sigma$ to $\sigma_{R}$ and present reliable radial profiles of $\sigma_{R}$ as well as accurate measurements of $\langle\sigma_{R}\rangle$, the radial average of $\sigma_{R}$ over one effective (half-light) radius. We show that there is a trend for $\sigma_{R}$ to increase with decreasing $R$, that $\langle\sigma_{R}\rangle$ correlates with stellar mass ($M_{\star}$) and tested correlations with other galaxy properties. The most significant and strongest correlation is the one with $M_{\star}$: $\langle\sigma_{R}\rangle \propto M_{\star}^{0.5}$. This tight scaling relation is applicable to spiral galaxies of type Sa $\mbox{--}$ Sd and stellar mass $M_{\star}\approx10^{9.5}\mbox{--}10^{11.5}\ \mbox{M}_{\odot}$. Simple models that relate $\sigma_{R}$ to the stellar surface density and disc scale length roughly reproduce that scaling, but overestimate $\langle\sigma_{R}\rangle$ significantly., Comment: MNRAS, in press

Using hydrodynamical simulations of entire galactic discs similar to the Milky Way, reaching 4.6pc resolution, we study the origins of observed physical properties of giant molecular clouds (GMCs). We find that efficient stellar feedback is a necessary ingredient in order to develop a realistic interstellar medium (ISM), leading to molecular cloud masses, sizes, velocity dispersions and virial parameters in excellent agreement with Milky Way observations. GMC scaling relations observed in the Milky Way, such as the mass-size ($M$--$R$), velocity dispersion-size ($\sigma$--$R$), and the $\sigma$--$R\Sigma$ relations, are reproduced in a feedback driven ISM when observed in projection, with $M\propto R^{2.3}$ and $\sigma\propto R^{0.56}$. When analysed in 3D, GMC scaling relations steepen significantly, indicating potential limitations of our understanding of molecular cloud 3D structure from observations. Furthermore, we demonstrate how a GMC population's underlying distribution of virial parameters can strongly influence the scatter in derived scaling relations. Finally, we show that GMCs with nearly identical global properties exist in different evolutionary stages, where a majority of clouds being either gravitationally bound or expanding, but with a significant fraction being compressed by external ISM pressure, at all times., Comment: 12 pages, 8 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

We introduce simulations aimed at assessing how well weak gravitational lensing of 21cm radiation from the Epoch of Reionization ($z \sim 8$) can be measured by an SKA-like radio telescope. A simulation pipeline has been implemented to study the performance of lensing reconstruction techniques. We show how well the lensing signal can be reconstructed using the three-dimensional quadratic lensing estimator in Fourier space assuming different survey strategies. The numerical code introduced in this work is capable of dealing with issues that can not be treated analytically such as the discreteness of visibility measurements and the inclusion of a realistic model for the antennae distribution. This paves the way for future numerical studies implementing more realistic reionization models, foreground subtraction schemes, and testing the performance of lensing estimators that take into account the non-Gaussian distribution of HI after reionization. If multiple frequency channels covering $z \sim 7-11.6$ are combined, Phase 1 of SKA-Low should be able to obtain good quality images of the lensing potential with a total resolution of $\sim 1.6$ arcmin. The SKA-Low Phase 2 should be capable of providing images with high-fidelity even using data from $z\sim 7.7 - 8.3$. We perform tests aimed at evaluating the numerical implementation of the mapping reconstruction. We also discuss the possibility of measuring an accurate lensing power spectrum. Combining data from $z \sim 7-11.6$ using the SKA2-Low telescope model, we find constraints comparable to sample variance in the range $L<1000$, even for survey areas as small as $25\mbox{ deg}^2$., Comment: 23 pages, 18 figures. Accepted for pubblication in MNRAS

We have proposed a method for measuring weak lensing using the Lyman-alpha forest. Here we estimate the noise expected in weak lensing maps and power spectra for different sets of observational parameters. We find that surveys of the size and quality of the ones being done today and ones planned for the future will be able to measure the lensing power spectrum at a source redshift of z~2.5 with high precision and even be able to image the distribution of foreground matter with high fidelity on degree scales. For example, we predict that Lyman-alpha forest lensing measurement from the Dark Energy Spectroscopic Instrument survey should yield the mass fluctuation amplitude with statistical errors of 1.5%. By dividing the redshift range into multiple bins some tomographic lensing information should be accessible as well. This would allow for cosmological lensing measurements at higher redshift than are accessible with galaxy shear surveys and correspondingly better constraints on the evolution of dark energy at relatively early times., Comment: 8 pages, 8 figures, submitted to MNRAS

The angular positions of quasars are deflected by the gravitational lensing effect of foreground matter. The Lyman-alpha forest seen in the spectra of these quasars is therefore also lensed. We propose that the signature of weak gravitational lensing of the forest could be measured using similar techniques that have been applied to the lensed Cosmic Microwave Background, and which have also been proposed for application to spectral data from 21cm radio telescopes. As with 21cm data, the forest has the advantage of spectral information, potentially yielding many lensed "slices" at different redshifts. We perform an illustrative idealized test, generating a high resolution angular grid of quasars (of order arcminute separation), and lensing the Lyman-alphaforest spectra at redshifts z=2-3 using a foreground density field. We find that standard quadratic estimators can be used to reconstruct images of the foreground mass distribution at z~1. There currently exists a wealth of Lya forest data from quasar and galaxy spectral surveys, with smaller sightline separations expected in the future. Lyman-alpha forest lensing is sensitive to the foreground mass distribution at redshifts intermediate between CMB lensing and galaxy shear, and avoids the difficulties of shape measurement associated with the latter. With further refinement and application of mass reconstruction techniques, weak gravitational lensing of the high redshift Lya forest may become a useful new cosmological probe., Comment: 9 pages, 7 figures, submitted to MNRAS

The velocity dispersion of cold interstellar gas, sigma, is one of the quantities that most radically affect the onset of gravitational instabilities in galaxy discs, and the quantity that is most drastically approximated in stability analyses. Here we analyse the stability of a large sample of nearby star-forming spirals treating molecular gas, atomic gas and stars as three distinct components, and using radial profiles of sigma_CO and sigma_HI derived from HERACLES and THINGS observations. We show that the radial variations of sigma_CO and sigma_HI have a weak effect on the local stability level of galaxy discs, which remains remarkably flat and well above unity, but is low enough to ensure (marginal) instability against non-axisymmetric perturbations and gas dissipation. More importantly, the radial variation of sigma_CO has a strong impact on the size of the regions over which gravitational instabilities develop, and results in a characteristic instability scale that is one order of magnitude larger than the Toomre length of molecular gas. Disc instabilities are driven, in fact, by the self-gravity of stars at kpc scales. This is true across the entire optical disc of every galaxy in the sample, with few exceptions. In the linear phase of the disc instability process, stars and molecular gas are strongly coupled, and it is such a coupling that ultimately triggers local gravitational collapse/fragmentation in the molecular gas., Comment: MNRAS, in press. Moderate revision to match the accepted version

We study the impact of stellar feedback in shaping the density and velocity structure of neutral hydrogen (HI) in disc galaxies. For our analysis, we carry out $\sim 4.6$pc resolution $N$-body+adaptive mesh refinement (AMR) hydrodynamic simulations of isolated galaxies, set up to mimic a Milky Way (MW), and a Large and Small Magellanic Cloud (LMC, SMC). We quantify the density and velocity structure of the interstellar medium using power spectra and compare the simulated galaxies to observed HI in local spiral galaxies from THINGS (The HI Nearby Galaxy Survey). Our models with stellar feedback give an excellent match to the observed THINGS HI density power spectra. We find that kinetic energy power spectra in feedback regulated galaxies, regardless of galaxy mass and size, show scalings in excellent agreement with super-sonic turbulence ($E(k)\propto k^{-2}$) on scales below the thickness of the HI layer. We show that feedback influences the gas density field, and drives gas turbulence, up to large (kpc) scales. This is in stark contrast to density fields generated by large scale gravity-only driven turbulence. We conclude that the neutral gas content of galaxies carries signatures of stellar feedback on all scales., Comment: 19 pages, 13 figures, 2 tables, accepted for publication in Monthly Notices of the Royal Astronomical Society

The mass-metallicity relation shows that the galaxies with the lowest mass have the lowest metallicities. As most dwarf galaxies are in group environments, interaction effects such as tides could contribute to this trend. We perform a series of smoothed particle hydrodynamics (SPH) simulations of dwarf galaxies in external tidal fields to examine the effects of tides on their metallicities and metallicity gradients. In our simulated galaxies, gravitational instabilities drive gas inwards and produce centralized star formation and a significant metallicity gradient. Strong tides can contribute to these instabilities, but their primary effect is to strip the outer low-metallicity gas, producing a truncated gas disk with a large metallicity. This suggests that the role of tides on the mass-metallicity relation is to move dwarf galaxies to higher metallicities., Comment: Accepted to ApJ

We explore the role that gravitational instability plays in NGC 1068, a nearby Seyfert galaxy that exhibits unusually vigorous starburst activity. For this purpose, we use the Romeo-Falstad disc instability diagnostics and data from BIMA SONG, SDSS and SAURON. Our analysis illustrates that NGC 1068 is a gravitationally unstable "monster". Its starburst disc is subject to unusually powerful instabilities. Several processes, including AGN/stellar feedback, try to quench such instabilities from inside out by depressing the surface density of molecular gas across the central kpc, but they do not succeed. Gravitational instability "wins" because it is driven by the stars via their much higher surface density. In this process, stars and molecular gas are strongly coupled, and it is such a coupling that ultimately triggers local gravitational collapse/fragmentation in the molecular gas., Comment: MNRAS, in press. Moderate revision to match the accepted version

We present an unprecedented measurement of the disc stability and local instability scales in the luminous infrared Seyfert 1 host, NGC7469, based on ALMA observations of dense gas tracers and with a synthesized beam of 165 x 132 pc. While we confirm that non-circular motions are not significant in redistributing the dense interstellar gas in this galaxy, we find compelling evidence that the dense gas is a suitable tracer for studying the origin of its intensely high-mass star forming ring-like structure. Our derived disc stability parameter accounts for a thick disc structure and its value falls below unity at the radii in which intense star formation is found. Furthermore, we derive the characteristic instability scale and find a striking agreement between our measured scale of ~ 180 pc, and the typical sizes of individual complexes of young and massive star clusters seen in high-resolution images., Comment: Accepted for publication in ApJ Letters

The late-type spiral galaxy NGC 6946 is a prime example of molecular gas dynamics driven by "bars within bars". Here we use data from the BIMA SONG and HERACLES surveys to analyse the structure and stability of its molecular disc. Our radial profiles exhibit a clear transition at distance R ~ 1 kpc from the galaxy centre. In particular, the surface density profile breaks at R ~ 0.8 kpc and is well fitted by a double exponential distribution with scale lengths R_1 ~ 200 pc and R_2 ~ 3 kpc, while the 1D velocity dispersion sigma decreases steeply in the central kpc and is approximately constant at larger radii. The fact that we derive and use the full radial profile of sigma rather than a constant value is perhaps the most novel feature of our stability analysis. We show that the profile of the Q stability parameter traced by CO emission is remarkably flat and well above unity, while the characteristic instability wavelength exhibits clear signatures of the nuclear starburst and inner bar within bar. We also show that CO-dark molecular gas, stars and other factors can play a significant role in the stability scenario of NGC 6946. Our results provide strong evidence that gravitational instability, radial inflow and disc heating have driven the formation of the inner structures and the dynamics of molecular gas in the central kpc., Comment: MNRAS, in press

Gravitational instabilities play an important role in galaxy evolution and in shaping the interstellar medium (ISM). The ISM is observed to be highly turbulent, meaning that observables like the gas surface density and velocity dispersion depend on the size of the region over which they are measured. In this work we investigate, using simulations of Milky Way-like disc galaxies with a resolution of $\sim 9$ pc, the nature of turbulence in the ISM and how this affects the gravitational stability of galaxies. By accounting for the measured average turbulent scalings of the density and velocity fields in the stability analysis, we can more robustly characterize the average level of stability of the galaxies as a function of scale, and in a straightforward manner identify scales prone to fragmentation. Furthermore, we find that the stability of a disc with feedback-driven turbulence can be well described by a "Toomre-like" $Q$ stability criterion on all scales, whereas the classical $Q$ can formally lose its meaning on small scales if violent disc instabilities occur in models lacking pressure support from stellar feedback., Comment: 11 pages, 5 figures, submitted to MNRAS

Gravitational instabilities play a primary role in shaping the clumpy structure and powering the star formation activity of gas-rich high-redshift galaxies. Here we analyse the stability of such systems, focusing on the size and mass ranges of unstable regions in the disc. Our analysis takes into account the mass-size and linewidth-size scaling relations observed in molecular gas, originally discovered by Larson. We show that such relations can have a strong impact on the size and mass of star-forming clumps, as well as on the stability properties of the disc at all observable scales, making the classical Toomre parameter a highly unreliable indicator of gravitational instability. For instance, a disc with Q=1 can be far from marginal instability, while a disc with Q<<1 can be marginally unstable. Our work raises an important caveat: if clumpy discs at high redshift have scale-dependent surface densities and velocity dispersions, as implied by the observed clump scaling relations, then we cannot thoroughly understand their stability and star formation properties unless we perform multi-scale observations. This will soon be possible thanks to dedicated ALMA surveys, which will explore the physical properties of super-giant molecular clouds at the peak of cosmic star formation and beyond., Comment: MNRAS, in press

In this paper, we propose a Q stability parameter that is more realistic than those commonly used, and is easy to evaluate [see Eq. (19)]. Using our Q_N parameter, you can take into account several stellar and/or gaseous components as well as the stabilizing effect of disc thickness, you can predict which component dominates the local stability level, and you can do all that simply and accurately. To illustrate the strength of Q_N, we analyse the stability of a large sample of spirals from The HI Nearby Galaxy Survey (THINGS), treating stars, HI and H_2 as three distinct components. Our analysis shows that H_2 plays a significant role in disc (in)stability even at distances as large as half the optical radius. This is an important aspect of the problem, which was missed by previous (two-component) analyses of THINGS spirals. We also show that HI plays a negligible role up to the edge of the optical disc; and that the stability level of THINGS spirals is, on average, remarkably flat and well above unity., Comment: MNRAS, in press

We investigate the gravitational instability of galactic discs, treating stars and cold interstellar gas as two distinct components, and taking into account the phenomenology of turbulence in the interstellar medium (ISM), i.e. the Larson-type scaling relations observed in the molecular and atomic gas. Besides deriving general properties of such systems, we analyse a large sample of galaxies from The HI Nearby Galaxy Survey (THINGS), and show in detail how interstellar turbulence affects disc instability in star-forming spirals. We find that turbulence has a significant effect on both the inner and the outer regions of the disc. In particular, it drives the inner gas disc to a regime of transition between two instability phases and makes the outer disc more prone to star-dominated instabilities., Comment: MNRAS, in press. Moderate revision to match the accepted version

The Wang-Silk approximation, 1/Q ~ 1/Q_stars + 1/Q_gas, is frequently used for estimating the effective Q parameter in two-component discs of stars and gas. Here we analyse this approximation in detail, and show how its accuracy depends on the radial velocity dispersions and Toomre parameters of the two components. We then propose a much more accurate but still simple approximation for the effective Q parameter, which further takes into account the stabilizing effect of disc thickness. Our effective Q parameter is a natural generalization of Toomre's Q, and as such can be used in a wide variety of contexts, e.g. for predicting star formation thresholds in galaxies or for measuring the stability level of galactic discs at low and high redshifts., Comment: MNRAS, in press

We explore the gravitational instability of clumpy and turbulent gas discs, taking into account the Larson-type scaling laws observed in giant molecular clouds (GMCs) and HI, as well as more general scaling relations. This degree of freedom is of special interest in view of the coming high-z ISM surveys, and is thus potentially important for understanding the dynamical effects of turbulence at all epochs of galaxy evolution. Our analysis shows that turbulence has a deep impact on the gravitational instability of the disc. It excites a rich variety of stability regimes, several of which have no classical counterpart. Among other diagnostics, we provide two useful tools for observers and simulators: (1) the stability map of turbulence, which illustrates our stability scenario and relates it to the phenomenology of interstellar turbulence: GMC/HI observations, simulations and models; (2) a Toomre-like stability criterion, $Q\geq\bar{Q}$, which applies to a large class of clumpy/turbulent discs. We make specific predictions about GMC and cold-HI turbulence, and point out the implications of our analysis for high-z galaxy surveys., Comment: MNRAS, in press. Moderate revision to match the accepted version

Observations of turbulent velocity dispersions in the HI component of galactic disks show a characteristic floor in galaxies with low star formation rates and within individual galaxies the dispersion profiles decline with radius. We carry out several high resolution adaptive mesh simulations of gaseous disks embedded within dark matter haloes to explore the roles of cooling, star-formation, feedback, shearing motions and baryon fraction in driving turbulent motions. In all simulations the disk slowly cools until gravitational and thermal instabilities give rise to a multi-phase medium in which a large population of dense self-gravitating cold clouds are embedded within a warm gaseous phase that forms through shock heating. The diffuse gas is highly turbulent and is an outcome of large scale driving of global non-axisymmetric modes as well as cloud-cloud tidal interactions and merging. At low star-formation rates these processes alone can explain the observed HI velocity dispersion profiles and the characteristic value of ~10 km/s observed within a wide range of disk galaxies. Supernovae feedback creates a significant hot gaseous phase and is an important driver of turbulence in galaxies with a star-formation rate per unit area >10^-3 M_sun/yr/kpc^2., Comment: 18 pages, 23 figures, MNRAS accepted. Typos and minor errors corrected. A version with high-resolution figures can be found at http://www-theorie.physik.unizh.ch/~agertz/DISK/

Melott has made pioneering studies of the effects of particle discreteness in N-body simulations, a fundamental point that needs careful thought and analysis since all such simulations suffer from numerical noise arising from the use of finite-mass particles. Melott (arXiv:0804.0589) claims that the conclusions of our paper (arXiv:0804.0294) are essentially equivalent to those of his earlier work. Melott is wrong: he has jumped onto one of our conclusions and interpreted that in his own way. Here we point out the whys and the wherefores.

The effects of particle discreteness in N-body simulations of Lambda Cold Dark Matter (LambdaCDM) are still an intensively debated issue. In this paper we explore such effects, taking into account the scatter caused by the randomness of the initial conditions, and focusing on the statistical properties of the cosmological density field. For this purpose, we run large sets of LambdaCDM simulations and analyse them using a large variety of diagnostics, including new and powerful wavelet statistics. Among other facts, we point out (1) that dynamical evolution does not propagate discreteness noise up from the small scales at which it is introduced, and (2) that one should aim to satisfy the condition epsilon ~ 2d, where "epsilon" is the force resolution and "d" is the interparticle distance. We clarify what such a condition means, and how to implement it in modern cosmological codes., Comment: ApJ, in press. Minor changes to match the accepted version

Wavelets are a new and powerful mathematical tool, whose most celebrated applications are data compression and de-noising. In Paper I (Romeo, Horellou & Bergh 2003, astro-ph/0302343), we have shown that wavelets can be used for removing noise efficiently from cosmological, galaxy and plasma N-body simulations. The expected two-orders-of-magnitude higher performance means, in terms of the well-known Moore's law, an advance of more than one decade in the future. In this paper, we describe a wavelet add-on code designed for such an application. Our code can be included in common grid-based N-body codes, is written in Fortran, is portable and available on request from the first author. The code can also be applied for removing noise from standard data, such as signals and images., Comment: Mon. Not. R. Astron. Soc., in press. The interested reader is strongly recommended to ignore the low-resolution Figs 10 and 11, and to download the full-resolution paper (800 kb) from http://www.oso.chalmers.se/~romeo/Paper_VII.ps.gz

Noise is a problem of major concern for N-body simulations of structure formation in the early Universe, of galaxies and plasmas. Here for the first time we use wavelets to remove noise from N-body simulations of disc galaxies, and show that they become equivalent to simulations with two orders of magnitude more particles. We expect a comparable improvement in performance for cosmological and plasma simulations. Our wavelet code will be described in a following paper, and will then be available on request., Comment: Mon. Not. R. Astron. Soc., in press. The interested reader is strongly recommended to ignore the low-resolution Fig. 3 (and Fig. 4), and to download the full-resolution paper (700 kb) from http://www.oso.chalmers.se/~romeo/Paper_VI.ps.gz