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: A&A accepted

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

It has been suggested that the gravitational potential can have a significant role in suppressing the star formation in the nearby galaxies. To establish observational constrains on this scenario, we investigate the connection between the dynamics, through the circular velocity curves (CVCs) as a proxy of the inner gravitational potential, and star formation quenching in 215 non-active galaxies across Hubble sequence from the Calar Alto Legacy Integral Field Area (CALIFA) survey. Our results show that galaxies with similar CVCs tend to have a certain star-formation quenching pattern. To explore these findings in more details, we construct kpc-resolved relations of the equivalent width of the H$\alpha$ ($W_{{\rm H}\alpha}$) versus the amplitude ($V_c$) and shape ($\beta= d\ln V_c/ d\ln R$) of the circular velocity at given radius. We find that the $W_{{\rm H}\alpha}-V_c$ is a declining relationship, where the retired regions of the galaxies (the ones with $W_{{\rm H}\alpha}$ values below 3 \r{A}) tend to have higher $V_c$. Differently, $W_{{\rm H}\alpha}-\beta$ is a bi-modal relationship, characterised by two peaks: concentration of the star forming regions at a positive $\beta$ (rising CVC) and another one of the retired regions with a negative $\beta$ (declining CVC). Our results show that both the amplitude of the CVC, driven by the mass of the galaxies, and the shape of the CVC, reflecting the internal structure of the galaxies, play an important role in galaxy's quenching history., Comment: 16 pages, 12 figures, Accepted for publication in A&A

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

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

We use the results of a high resolution N-body simulation to investigate the role of the environment on the formation and evolution of galaxy-sized halos. Starting from a set of constrained initial conditions, we have produced a final configuration hosting a double cluster in one octant and a large void extending over two octants of the simulation box. We present results for two statistics: the relationship between 1-D velocity dispersion and mass and the probability distribution of the spin parameter $P(\lambda)$. The \svm relationship is well reproduced by the Truncated Isothermal Sphere (TIS) model introduced by Shapiro et al. (1999), although the slope is different from the original prediction. A series of \svm relationships for different values of the anisotropy parameter $\beta$, obtained using the theoretical predictions by Lokas and Mamon (2001) for NFW density profiles are found to be only marginally consistent with the data. Using some properties of the equilibrium TIS models, we construct subsamples of {\em fiducial} equilibrium TIS halos from each of the three subregions, and we study their properties. For these halos, we do find an environmental dependence of their properties, in particular of the spin parameter distribution $P(\lambda)$. We study in more detail the TIS model, and we find new relationships between the truncation radius and other structural parameters. No gravitationally bound halo is found having a radius larger than the critical value for gravithermal instability for TIS halos (\rt $\ge 34.2 r_{0}$, where $r_{0}$ is the core radius of the TIS solution). We do however find a dependence of this relationship on the environment, like for the $P(\lambda)$ statistics. These facts hint at a possible r\^{o}le of tidal fields at determining the statistical properties of halos., Comment: 12 pages, 14 figures. Accepted by MNRAS. Adopted an improved algorithm for halo finding and added a comparison with NFW models

The workshop `Astrophysical Dynamics 1999/2000' followed a homonymous advanced research course, and both activities were organized by me. In this opening paper of the proceedings book, I describe them and document their strong impact on the academic life of the local institutions. The advanced research course was open to graduate students, senior researchers, and motivated under-graduate students with good background in physics and mathematics. The course covered several multi-disciplinary issues of modern research on astrophysical dynamics, and thus also of interest to physicists, mathematicians and engineers. The major topic was gas dynamics, viewed in context with stellar dynamics and plasma physics. The course was complemented by parallel seminars on hot topics given by experts in such fields, and open to a wide scientific audience. In particular, I gave a friendly introduction to wavelets, which are becoming an increasingly powerful tool not only for processing signals and images but also for analysing fractals and turbulence, and which promise to have important applications to dynamical modelling of disc galaxies. The workshop was open to a wide scientific audience. The workshop with published proceedings book was, as a matter of fact, the innovative form of exam that I proposed for the advanced research course. The contributions were refereed and their quality is high on average, exceptionally high in a few cases. The advanced research course and the workshop all together produced great enthusiasm in the students and welcomed the participation of a hundred different people, which means an order of magnitude more than an average graduate course at Chalmers University of Technology and G\"oteborg University., Comment: opening paper; the proceedings book is in http://www.oso.chalmers.se/~romeo/PROCEEDINGS_BOOK_E

In this final report, I briefly reflect on two parallel teaching experiences as tutor. Besides, I briefly view such experiences in interaction with my research work, private life and new teaching position. In harmony with my conception of teaching, I avoid the standard formal style of reports and try an interactive dialogue with the reader., Comment: 9 pages (tex)

Modelling gravity is a fundamental problem that must be tackled in N-body simulations of stellar systems, and satisfactory solutions require a deep understanding of the dynamical effects of softening. In a previous paper (Romeo 1997), we have devised a method for exploring such effects, and we have focused on two applications that reveal the dynamical differences between the most representative types of softened gravity. In the present paper we show that our method can be applied in another, more fruitful, way: for developing new ideas about softening. Indeed, it opens a direct route to the discovery of optimal types of softened gravity for given dynamical requirements, and thus to the accomplishment of a physically consistent modelling of disc galaxies, even in the presence of a cold interstellar gaseous component and in situations that demand anisotropic resolution., Comment: A&A (in press). 7 pages (latex, A&A style), 2 postscript figures

Two questions that naturally arise in N-body simulations of stellar systems are: (1) How can we compare experiments that employ different types of softened gravity? (2) Given a particular type of softened gravity, which choices of the softening length optimize the faithfulness of the experiments to the Newtonian dynamics? We devise a method for exploring the dynamical effects of softening, which provides detailed answers in the case of 2-D simulations of disc galaxies and also solves important aspects of the 3-D problem. In the present paper we focus on two applications that reveal the dynamical differences between the most representative types of softened gravity, including certain anisotropic alternatives. Our method is potentially important not only for testing but also for developing new ideas about softening. Indeed, it opens a direct route to the discovery of optimal types of softened gravity for given dynamical requirements, and thus to the accomplishment of a physically consistent modelling., Comment: A&A (in press). 12 pages (latex, A&A style), 5 postscript figures