Source author record

Danail Obreschkow

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

48works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

48 published item(s)

preprint2022arXiv

An empirical measurement of the Halo Mass Function from the combination of GAMA DR4, SDSS DR12, and REFLEX II data

We construct the halo mass function (HMF) from the GAMA galaxy group catalogue over the mass range 10^12.7M_sol to 10^15.5M_sol, and find good agreement with the expectation from LambdaCDM. In comparison to previous studies, this result extends the mass range over which the HMF has now been measured over by an order of magnitude. We combine the GAMA DR4 HMF with similar data from the SDSS DR12 and REFLEX II surveys, and fit a four-parameter Murray-Robotham-Power (MRP) function, valid at z~0.1, yielding: a density normalisation of: log10 (phi Mpc^3)=-3.96[+0.55,-0.82], a high mass turn-over of: log10(M/M_sol)=14.13[+0.43,-0.40], a low mass power law slope of: alpha=-1.68[+0.21,-0.24] , and a high mass softening parameter of: beta= 0.63[+0.25,-0.11]. If we fold in the constraint on Omega_M from Planck 2018 Cosmology, we are able to reduce these uncertainties further, but this relies on the assumption that the power-law trend can be extrapolated from 10^12.7M_sol to zero mass. Throughout, we highlight the effort needed to improve on our HMF measurement: improved halo mass estimates that do not rely on calibration to simulations; reduced halo mass uncertainties needed to mitigate the strong Eddington Bias that arises from the steepness of the HMF low mass slope; and deeper wider area spectroscopic surveys. To our halo mass limit of 10^12.7 M_sol, we are directly resolving (`seeing') 41+/-5 per cent of the total mass density, i.e. Omega_[M>12.7]=0.128+/-0.016, opening the door for the direct construction of 3D dark matter mass maps at Mpc resolution.

preprint2022arXiv

Galaxy And Mass Assembly (GAMA): Data Release 4 and the z < 0.1 total and z < 0.08 morphological galaxy stellar mass functions

In Galaxy And Mass Assembly Data Release 4 (GAMA DR4), we make available our full spectroscopic redshift sample. This includes 248682 galaxy spectra, and, in combination with earlier surveys, results in 330542 redshifts across five sky regions covering ~250deg^2. The redshift density, is the highest available over such a sustained area, has exceptionally high completeness (95 per cent to r_KIDS=19.65mag), and is well suited for the study of galaxy mergers, galaxy groups, and the low redshift (z<0.25) galaxy population. DR4 includes 32 value-added tables or Data Management Units (DMUs) that provide a number of measured and derived data products including GALEX, ESO KiDS, ESO VIKING, WISE and Herschel Space Observatory imaging. Within this release, we provide visual morphologies for 15330 galaxies to z<0.08, photometric redshift estimates for all 18million objects to r_KIDS~25mag, and stellar velocity dispersions for 111830 galaxies. We conclude by deriving the total galaxy stellar mass function (GSMF) and its sub-division by morphological class (elliptical, compact-bulge and disc, diffuse-bulge and disc, and disc only). This extends our previous measurement of the total GSMF down to 10^6.75 M_sol h^-2_70 and we find a total stellar mass density of rho_*=(2.97+/-0.04)x10^8 M_sol h_70 Mpc^-3 or Omega_*=(2.17+/-0.03)x10^-3 h^-1_70. We conclude that at z<0.1, the Universe has converted 4.9+/-0.1 per cent of the baryonic mass implied by Big Bang Nucleosynthesis into stars that are gravitationally bound within the galaxy population.

preprint2022arXiv

LADUMA: Discovery of a luminous OH megamaser at $z > 0.5$

In the local Universe, OH megamasers (OHMs) are detected almost exclusively in infrared-luminous galaxies, with a prevalence that increases with IR luminosity, suggesting that they trace gas-rich galaxy mergers. Given the proximity of the rest frequencies of OH and the hyperfine transition of neutral atomic hydrogen (HI), radio surveys to probe the cosmic evolution of HI in galaxies also offer exciting prospects for exploiting OHMs to probe the cosmic history of gas-rich mergers. Using observations for the Looking At the Distant Universe with the MeerKAT Array (LADUMA) deep HI survey, we report the first untargeted detection of an OHM at $z > 0.5$, LADUMA J033046.20$-$275518.1 (nicknamed "Nkalakatha"). The host system, WISEA J033046.26$-$275518.3, is an infrared-luminous radio galaxy whose optical redshift $z \approx 0.52$ confirms the MeerKAT emission line detection as OH at a redshift $z_{\rm OH} = 0.5225 \pm 0.0001$ rather than HI at lower redshift. The detected spectral line has 18.4$σ$ peak significance, a width of $459 \pm 59\,{\rm km\,s^{-1}}$, and an integrated luminosity of $(6.31 \pm 0.18\,{\rm [statistical]}\,\pm 0.31\,{\rm [systematic]}) \times 10^3\,L_\odot$, placing it among the most luminous OHMs known. The galaxy's far-infrared luminosity $L_{\rm FIR} = (1.576 \pm 0.013) \times 10^{12}\,L_\odot$ marks it as an ultra-luminous infrared galaxy; its ratio of OH and infrared luminosities is similar to those for lower-redshift OHMs. A comparison between optical and OH redshifts offers a slight indication of an OH outflow. This detection represents the first step towards a systematic exploitation of OHMs as a tracer of galaxy growth at high redshifts.

preprint2022arXiv

Spin transfer from dark matter to gas during halo formation

In the protogalactic density field, diffuse gas and collision-less cold dark matter (DM) are often assumed sufficiently mixed that both components experience identical tidal torques. However, haloes in cosmological simulations consistently end up with a higher specific angular momentum (sAM) in gas, even in simulations without radiative cooling and galaxy formation physics. We refine this result by analysing the spin distributions of gas and DM in $\sim$50,000 well-resolved haloes in a non-radiative cosmological simulation from the SURFS suite. The sAM of the halo gas on average ends up $\sim$40\% above that of the DM. This can be pinned down to an excess AM in the inner halo ($<$50\% virial radius), paralleled by a more coherent rotation pattern in the gas. We uncover the leading driver for this AM difference through a series of control simulations of a collapsing ellipsoidal top-hat, where gas and DM are initially well mixed. These runs reveal that the pressurised inner gas shells collapse more slowly, causing the DM ellipsoid to spin ahead of the gas ellipsoid. The arising torque generally transfers AM from the DM to the gas. The amount of AM transferred via this mode depends on the initial spin, the initial axes ratios and the collapse factor. These quantities can be combined in a single dimensionless parameter, which robustly predicts the AM transfer of the ellipsoidal collapse. This simplistic model can quantitatively explain the average AM excess of the gas found in the more complex non-radiative cosmological simulation.

preprint2022arXiv

Stellar Masses of Clumps in Gas-rich, Turbulent Disk Galaxies

In this paper we use HST/WFC3 observations of 6 galaxies from the DYNAMO survey, combined with stellar population modelling of the SED, to determine the stellar masses of DYNAMO clumps. The DYNAMO sample has been shown to have properties similar to $z\approx1.5$ turbulent, clumpy disks. DYNAMO sample clump masses offer a useful comparison for studies of $z>1$ in that the galaxies have the same properties, yet the observational biases are significantly different. Using DYNAMO we can more easily probe rest-frame near-IR wavelengths and also probe finer spatial scales. We find that the stellar mass of DYNAMO clumps is typically $10^{7}-10^8 \mathrm{M}_\odot$. We employ a technique that makes non-parametric corrections in removal of light from nearby clumps, and carries out a locally determined disk subtraction. The process of disk subtraction is the dominant effect, and can alter clump masses at the 0.3~dex level. Using these masses, we investigate the stellar mass function of clumps in DYNAMO galaxies. DYNAMO stellar mass functions follow a declining power law with slope $α\approx -1.4$, which is slightly shallower than, but similar to what is observed in $z>1$ lensed galaxies. We compare DYNAMO clump masses to results of simulations. The masses and galactocentric position of clumps in DYNAMO galaxies are more similar to long-lived clumps in simulations. Similar to recent DYNAMO results on the stellar population gradients, these results are consistent with simulations that do not employ strong "early" radiative feedback prescriptions.

preprint2022arXiv

The variation of the gas content of galaxy groups and pairs compared to isolated galaxies

We measure how the atomic gas (HI) fraction ($f_{HI}={\rm \frac{M_{HI}}{M_{*}}}$) of groups and pairs taken as single units vary with average stellar mass ($\langle {\rm M_*} \rangle$) and average star-formation rate ($\langle {\rm SFR} \rangle$), compared to isolated galaxies. The HI 21 cm emission observation are from (i) archival ALFALFA survey data covering three fields from the GAMA survey (provides environmental and galaxy properties), and (ii) DINGO pilot survey data of one of those fields. The mean $f_{HI}$ for different units (groups/pairs/isolated galaxies) are measured in regions of the log($\langle {\rm M_*} \rangle$) -- log($\langle {\rm SFR} \rangle$) plane, relative to the z $\sim 0$ star-forming main sequence (SFMS) of individual galaxies, by stacking $f_{HI}$ spectra of individual units. For ALFALFA, $f_{HI}$ spectra of units are measured by extracting HI spectra over the full groups/pair areas and dividing by the total stellar mass of member galaxies. For DINGO, $f_{HI}$ spectra of units are measured by co-adding HI spectra of individual member galaxies, followed by division by their total stellar mass. For all units the mean $f_{HI}$ decreases as we move to higher $\langle {\rm M_*} \rangle$ along the SFMS, and as we move from above the SFMS to below it at any $\langle {\rm M_*} \rangle$. From the DINGO-based study, mean $f_{HI}$ in groups appears to be lower compared to isolated galaxies for all $\langle {\rm M_*} \rangle$ along the SFMS. From the ALFALFA-based study we find substantially higher mean $f_{HI}$ in groups compared to isolated galaxies (values for pairs being intermediate) for ${\langle{\rm M_*}\rangle}\lesssim10^{9.5}~{\rm M_{\odot}}$, indicating the presence of substantial amounts of HI not associated with cataloged member galaxies in low mass groups.

preprint2022arXiv

xGASS: The connection between angular momentum, mass and atomic gas fraction in nearby galaxies

We use a sample of 559 disc galaxies extracted from the eXtended GALEX Arecibo SDSS Survey (xGASS) to study the connection between baryonic angular momentum, mass and atomic gas fraction in the local Universe. Baryonic angular momenta are determined by combining HI and H$_{2}$ integrated profiles with two-dimensional stellar mass surface density profiles. In line with previous work, we confirm that specific angular momentum and atomic gas fraction are tightly correlated, but we find a larger scatter than previously observed. This is most likely due to the wider range of galaxy properties covered by our sample. We compare our findings with the predictions of the analytical stability model developed by Obreschkow et al. and find that, while the model provides a very good first-order approximation for the connection between baryonic angular momentum, mass and gas fraction, it does not fully match our data. Specifically, we find that at fixed baryonic mass, the dependence of specific angular momentum on gas fraction is significantly weaker, and at fixed gas fraction, the slope of the angular momentum vs. mass relation is shallower than what was predicted by the model. The reasons behind this tension remain unclear, but we speculate that multiple factors may simultaneously play a role, all related to the fact that the model is not able to encapsulate the full diversity of galaxy properties in our sample.

preprint2021arXiv

GAMA/DEVILS: Constraining the cosmic star-formation history from improved measurements of the 0.3-2.2 micron Extragalactic Background Light

We present a revised measurement of the optical extragalactic background light (EBL), based on the contribution of resolved galaxies to the integrated galaxy light (IGL). The cosmic optical background radiation (COB), encodes the light generated by star-formation, and provides a wealth of information about the cosmic star formation history (CSFH). We combine wide and deep galaxy number counts from the Galaxy And Mass Assembly survey (GAMA) and Deep Extragalactic VIsible Legacy Survey (DEVILS), along with the Hubble Space Telescope (HST) archive and other deep survey datasets, in 9 multi-wavelength filters to measure the COB in the range from 0.35 micron to 2.2 micron. We derive the luminosity density in each band independently and show good agreement with recent and complementary estimates of the optical-EBL from very high-energy (VHE) experiments. Our error analysis suggests that the IGL and Gamma-ray measurements are now fully consistent to within ~10%, suggesting little need for any additional source of diffuse light beyond the known galaxy population. We use our revised IGL measurements to constrain the cosmic star-formation history, and place amplitude constraints on a number of recent estimates. As a consistency check, we can now demonstrate convincingly, that the CSFH, stellar mass growth, and the optical-EBL provide a fully consistent picture of galaxy evolution. We conclude that the peak of star-formation rate lies in the range 0.066-0.076 Msol/yr/Mpc^3 at a lookback time of 9.1 to 10.9 Gyrs.

preprint2021arXiv

K-Band Imaging of the Nearby, Clumpy Turbulent Disk Galaxy DYNAMO G04-1

We present a case study of stellar clumps in G04-1, a clumpy, turbulent disk galaxy located at $z$ = 0.13 from the DYNAMO sample, using adaptive optics enabled K-band imaging ($\sim2.25$ kpc/arcsec) with Keck/NIRC2. We identify 15 stellar clumps in G04-1 with a range of masses from $3.6 \times 10^{6}$ to $2.7 \times 10^{8}\ \rm M_{\odot}$, and with a median mass of $\sim2.9 \times 10^{7}\ \rm M_{\odot}$. Note that these masses decrease by about one-half when we apply a light correction for the underlying stellar disk. A majority (12 of 15) of clumps observed in the $K_{P}$-band imaging have associated components in H$α$ maps ($\sim2.75$ kpc/arcsec; $<$R$_{clump}> \sim$500 pc) and appear co-located ($\overline{Δx} \sim 0.1$ arcsec). Using Hubble Space Telescope WFC/ACS observations with the F336W and F467M filters, we also find evidence of radial trends in clump stellar properties: clumps closer to the centre of G04-1 are more massive (consistent with observations at high-$z$) and appear more red, suggesting they may be more evolved. Using our high-resolution data, we construct a star forming main sequence for G04-1 in terms of spatially-resolved quantities and find that all regions (both clump and intra-clump) within the galaxy are experiencing an enhanced mode of star formation routinely observed in galaxies at high-$z$. In comparison to recent simulations, our observation of a number of clumps with masses $10^{7}-10^{8}\ \rm M_{\odot}$ is not consistent with strong radiative feedback in this galaxy.

preprint2021arXiv

The SAMI Galaxy Survey: the third and final data release

We have entered a new era where integral-field spectroscopic surveys of galaxies are sufficiently large to adequately sample large-scale structure over a cosmologically significant volume. This was the primary design goal of the SAMI Galaxy Survey. Here, in Data Release 3 (DR3), we release data for the full sample of 3068 unique galaxies observed. This includes the SAMI cluster sample of 888 unique galaxies for the first time. For each galaxy, there are two primary spectral cubes covering the blue (370-570nm) and red (630-740nm) optical wavelength ranges at spectral resolving power of R=1808 and 4304 respectively. For each primary cube, we also provide three spatially binned spectral cubes and a set of standardized aperture spectra. For each galaxy, we include complete 2D maps from parameterized fitting to the emission-line and absorption-line spectral data. These maps provide information on the gas ionization and kinematics, stellar kinematics and populations, and more. All data are available online through Australian Astronomical Optics (AAO) Data Central.

preprint2020arXiv

Angular momentum-related probe of cold gas deficiencies

Recent studies of neutral atomic hydrogen (HI) in nearby galaxies found that all field disk galaxies are HI saturated, in that they carry roughly as much HI as permitted before this gas becomes gravitationally unstable. By taking this HI saturation for granted, the atomic gas fraction $f_{\rm atm}$ of galactic disks can be predicted as a function of the stability parameter $q=jσ/(GM)$, where $M$ and $j$ are the baryonic mass and specific angular momentum of the disk and $σ$ is the HI velocity dispersion Obreschkow et al. 2016. The log-ratio $Δf_q$ between this predictor and the observed atomic fraction can be seen as a physically motivated `HI deficiency'. While field disk galaxies have $Δf_q \approx0$, objects subject to environmental removal of HI are expected to have $Δf_q>0$. Within this framework, we revisit the HI deficiencies of satellite galaxies in the Virgo cluster and in clusters of the EAGLE simulation. We find that observed and simulated cluster galaxies are HI deficient and that $Δf_q$ slightly increases when getting closer to the cluster centres. The $Δf_q$ values are similar to traditional HI deficiency estimators, but $Δf_q$ is more directly comparable between observations and simulations than morphology-based deficiency estimators. By tracking the simulated HI deficient cluster galaxies back in time, we confirm that $Δf_q\approx0$ until the galaxies first enter a halo with $M_{\rm halo}>10^{13} {\rm M_{\odot}}$, at which moment they quickly lose HI by environmental effects. Finally, we use the simulation to investigate the links between $Δf_q$ and quenching of star formation.

preprint2020arXiv

Characterising the Structure of Halo Merger Trees Using a Single Parameter: The Tree Entropy

Linking the properties of galaxies to the assembly history of their dark matter haloes is a central aim of galaxy evolution theory. This paper introduces a dimensionless parameter $s\in[0,1]$, the "tree entropy", to parametrise the geometry of a halo's entire mass assembly hierarchy, building on a generalisation of Shannon's information entropy. By construction, the minimum entropy ($s=0$) corresponds to smoothly assembled haloes without any mergers. In contrast, the highest entropy ($s=1$) represents haloes grown purely by equal-mass binary mergers. Using simulated merger trees extracted from the cosmological $N$-body simulation SURFS, we compute the natural distribution of $s$, a skewed bell curve peaking near $s=0.4$. This distribution exhibits weak dependences on halo mass $M$ and redshift $z$, which can be reduced to a single dependence on the relative peak height $δ_{\rm c}/σ(M,z)$ in the matter perturbation field. By exploring the correlations between $s$ and global galaxy properties generated by the SHARK semi-analytic model, we find that $s$ contains a significant amount of information on the morphology of galaxies $-$ in fact more information than the spin, concentration and assembly time of the halo. Therefore, the tree entropy provides an information-rich link between galaxies and their dark matter haloes.

preprint2020arXiv

From rest-frame luminosity functions to observer-frame colour distributions: tackling the next challenge in cosmological simulations

Galaxy spectral energy distributions (SEDs) remain among the most challenging yet informative quantities to reproduce in simulations due to the large and complex mixture of physical processes that shape the radiation output of a galaxy. With the increasing number of surveys utilising broadband colours as part of their target selection criteria, the production of realistic SEDs in simulations is necessary for assisting in survey design and interpretation of observations. The recent success in reproducing the observed luminosity functions (LF) from far-UV to far-IR, using the state-of-the-art semi-analytic model \shark\ and the SED generator \prospect, represents a critical step towards better galaxy colour predictions. We show that with \shark\ and \prospect\ we can closely reproduce the optical colour distributions observed in the panchromatic GAMA survey. The treatment of feedback, star formation, central-satellite interactions and radiation re-processing by dust are critical for this achievement. The first three processes create a bimodal distribution, while dust attenuation defines the location and shape of the blue and red populations. While a naive comparison between observation and simulations displays the known issue of over-quenching of satellite galaxies, the introduction of empirically-motivated observational errors and classification from the same group finder used in GAMA greatly reduces this tension. The introduction of random re-assignment of $\sim 15\%$ of centrals/satellites as satellites/centrals on the simulation classification closely resembles the outcome of the group finder, providing a computationally less intensive method to compare simulations with observations.

preprint2020arXiv

Multi-wavelength consensus of large-scale linear bias

We model the large-scale linear galaxy bias $b_g(x,z)$ as a function of redshift $z$ and observed absolute magnitude threshold $x$ for broadband continuum emission from the far infrared to ultra-violet, as well as for prominent emission lines, such as the H$α$, H$β$, Lya and [OII] lines. The modelling relies on the semi-analytic galaxy formation model GALFORM, run on the state-of-the-art $N$-body simulation SURFS with the Planck 2015 cosmology. We find that both the differential bias at observed absolute magnitude $x$ and the cumulative bias for magnitudes brighter than $x$ can be fitted with a five-parameter model: $b_g(x,z)=a + b(1+z)^e(1 + \exp{[(x-c)d]})$. We also find that the bias for the continuum bands follows a very similar form regardless of wavelength due to the mixing of star-forming and quiescent galaxies in a magnitude limited survey. Differences in bias only become apparent when an additional colour separation is included, which suggest extensions to this work could look at different colours at fixed magnitude limits. We test our fitting formula against observations, finding reasonable agreement with some measurements within $1σ$ statistical uncertainties, and highlighting areas of improvement. We provide the fitting parameters for various continuum bands, emission lines and intrinsic galaxy properties, enabling a quick estimation of the linear bias in any typical survey of large-scale structure.

preprint2020arXiv

Stellar angular momentum distribution linked to galaxy morphology

We study the spatially-resolved stellar specific angular momentum $j_*$ in a high-quality sample of 24 CALIFA galaxies covering a broad range of visual morphology, accounting for stellar velocity and velocity dispersion. The shape of the spaxel-wise probability density function of normalised $s=j_*/j_{*mean}$, PDF($s$), deviates significantly from the near-universal initial distribution expected of baryons in a dark matter halo and can be explained by the expected baryonic effects in galaxy formation that remove and redistribute angular momentum. Further we find that the observed shape of the PDF($s$) correlates significantly with photometric morphology, where late-type galaxies have PDF($s$) that is similar to a normal distribution, whereas early types have a strongly-skewed PDF($s$) resulting from an excess of low-angular momentum material. Galaxies that are known to host pseudobulges (bulge Sersic index $n_b <2.2$) tend to have less skewed bulge PDF($s$), with skewness $(b_{1rb})\lesssim0.8$. The PDF($s$) encodes both kinematic and photometric information and appears to be a robust tracer of morphology. Its use is motivated by the desire to move away from traditional component-based classifications which are subject to observer bias, to classification on a galaxy's fundamental (stellar mass, angular momentum) properties. In future, PDF($s$) may also be useful as a kinematic decomposition tool.

preprint2020arXiv

The influence of angular momentum and environment on the HI gas of late-type galaxies

We use high-resolution HI data from the WHISP survey to study the HI and angular momentum properties of a sample of 114 late-type galaxies. We explore the specific baryonic angular momentum -- baryonic mass ($j_b - M_b$) relation, and find that an unbroken power law of the form $j_b \propto M_b^{0.55 \pm 0.02}$ fits the data well, with an intrinsic scatter of $\sim 0.13 \pm 0.01$ dex. We revisit the relation between the atomic gas fraction, $f_{atm}$, and the integrated atomic stability parameter $q$ (the $f_{atm} - q$ relation), originally introduced by Obreschkow et al., and probe this parameter space by populating it with galaxies from different environments, in order to study the influence of the environment on their $j_b$, $f_{atm}$ and $q$ values. We find evidence that galaxies with close neighbours show a larger intrinsic scatter about the $f_{atm} - q$ relation compared to galaxies without close-neighbours. We also find enhanced SFR among the deviating galaxies with close neighbours. In addition, we use the bulge-to-total (B/T) ratio as a morphology proxy, and find a general trend of decreasing B/T values with increasing disc stability and HI fraction in the $f_{atm} - q$ plane, indicating a fundamental link between mass, specific angular momentum, gas fraction and morphology of galaxies.

preprint2020arXiv

The physical drivers of the atomic hydrogen-halo mass relation

We use SHARK, a semi-analytic galaxy formation model, to investigate the physical processes involved in dictating the shape, scatter and evolution of the HI-halo mass relation at $0\leq z \leq 2$. We compare SHARK with HI clustering and spectral stacking of the HI-halo mass relation derived from observations finding excellent agreement with the former and a deficiency of HI in SHARK at $M_{\rm vir}\approx 10^{12-13} M_{\odot}$ in the latter, but otherwise great agreement below and above that mass threshold. In SHARK, we find that the HI mass increases with the halo mass up to a critical mass of $\approx 10^{11.8} M_{\odot}$; between $\sim 10^{11.8}-10^{13}M_{\odot}$, the scatter in the relation increases by 0.7 dex and the HI mass decreases with the halo mass on average; at $M_{\rm vir} \geq 10^{13} M_{\odot}$, the HI content continues to increase with halo mass. We find that the critical halo mass of $\approx 10^{12} M_{\odot}$ is largely set by feedback from Active Galactic Nuclei (AGN), and the exact shape and scatter of the HI-halo mass relation around that mass is extremely sensitive to how AGN feedback is modelled, with other physical processes playing a less significant role. We determine the main secondary parameters responsible for the scatter of the HI-halo mass relation, namely the halo spin parameter at $M_{\rm vir}\leq 10^{11.8} M_{\odot}$, and the fractional contribution from substructure to the total halo mass for $M_{\rm vir}\geq 10^{13} M_{\odot}$. The scatter at $10^{11.8}<M_{\rm vir}<10^{13} M_{\odot}$ is best described by the black-hole-to-stellar mass ratio of the central galaxy, reflecting the AGN feedback relevance. We present a numerical model to populate dark matter-only simulations with HI at $0\leq z \leq 2$ based solely on halo parameters that are measurable in such simulations.

preprint2019arXiv

Fade to grey: systematic variation of the galaxy attenuation curves with galaxy properties in EAGLE

We present a simple model for galaxy attenuation by distilling SKIRT radiative transfer calculations for ~100,000 EAGLE galaxies at redshifts z=2-0. Our model adapts the two component screen model of Charlot & Fall (2000), parametrising the optical depth and slope of the ISM screen using the average dust surface density, $Σ_{\rm dust}$. We recover relatively tight relations between these parameters for the EAGLE sample, but also provide the scatter in these parameter owing to the morphological variation and orientation of galaxies. We also find that these relations are nearly independent of redshift in the EAGLE model. By pairing our model with an empirical prescription for birth clouds below the resolution scale of the simulation, we reproduce the observed relation between attenuation slope and optical depth for the first time in a cosmological simulation. We demonstrate that this result is remarkably independent of the attenuation properties assumed for birth cloud screen, merely requiring a boosted attenuation for infant stars. We present this model with a view to interpreting observations, as well as processing semi-analytic models and other hydrodynamic simulations.

preprint2019arXiv

The SAMI Galaxy Survey: The contribution of different kinematic classes to the stellar mass function of nearby galaxies

We use the complete Sydney-AAO Multi-object Integral field spectrograph (SAMI) Galaxy Survey to determine the contribution of slow rotators, as well as different types of fast rotators, to the stellar mass function of galaxies in the local Universe. We use stellar kinematics not only to discriminate between fast and slow rotators, but also to distinguish between dynamically cold systems (i.e., consistent with intrinsic axis ratios$<0.3$) and systems including a prominent dispersion-supported bulge. We show that fast rotators account for more than $80\%$ of the stellar mass budget of nearby galaxies, confirming that their number density overwhelms that of slow rotators at almost all masses from $10^{9}$ to $10^{11.5}{\rm M_\odot}$. Most importantly, dynamically cold disks contribute to at least $25\%$ of the stellar mass budget of the local Universe, significantly higher than what is estimated from visual morphology alone. For stellar masses up to $10^{10.5}{\rm M_\odot}$, this class makes up $>=30\%$ of the galaxy population in each stellar mass bin. The fact that many galaxies that are visually classified as having two-components have stellar spin consistent with dynamically cold disks suggests that the inner component is either rotationally-dominated (e.g., bar, pseudo-bulge) or has little effect on the global stellar kinematics of galaxies.

preprint2016arXiv

Angular Momentum of Dwarf Galaxies

We present measurements of baryonic mass Mb and specific angular momentum (sAM) jb in 14 rotating dwarf Irregular (dIrr) galaxies from the LITTLE THINGS sample. These measurements, based on 21cm kinematic data from the Very Large Array and stellar mass maps from the Spitzer Space Telescope, extend previous AM measurements by more than two orders of magnitude in Mb. The dwarf galaxies show systematically higher jb values than expected from the jb~Mb^{2/3} scaling of spiral galaxies, representative of a scale-free galaxy formation scenario. This offset can be explained by decreasing baryon mass fractions fM=Mb Mdyn (where Mdyn is the dynamical mass) with decreasing Mb (for Mb<10^{11}Msun). We find that the sAM of neutral atomic hydrogen HI alone is about 2.5 times higher than that of the stars. The M-j relation of HI is significantly steeper than that of the stars, as a direct consequence of the systematic variation of the HI fraction with Mb.

preprint2016arXiv

Angular Momentum Regulates Atomic Gas Fractions of Galactic Disks

We show that the mass fraction f_atm = 1.35*MHI/M of neutral atomic gas (HI and He) in isolated local disk galaxies of baryonic mass M is well described by a straightforward stability model for flat exponential disks. In the outer disk parts, where gas at the characteristic dispersion of the warm neutral medium is stable in the sense of Toomre (1964), the disk consists of neutral atomic gas; conversely the inner part where this medium would be Toomre-unstable, is dominated by stars and molecules. Within this model, f_atm only depends on a global stability parameter q=j*sigma/(GM), where j is the baryonic specific angular momentum of the disk and sigma the velocity dispersion of the atomic gas. The analytically derived first-order solution f_atm = min{1,2.5q^1.12} provides a good fit to all plausible rotation curves. This model, with no free parameters, agrees remarkably well (+-0.2 dex) with measurements of f_atm in isolated local disk galaxies, even with galaxies that are extremely HI-rich or HI-poor for their mass. The finding that f_atm increases monotonically with q for pure stability reasons offers a powerful intuitive explanation for the mean variation of f_atm with M: in a cold dark matter universe galaxies are expected to follow j~M^(2/3), which implies the average scaling q~M^(-1/3) and hence f_atm~M^(-0.37), in agreement with observations.

preprint2016arXiv

DYNAMO-HST Survey: Clumps in Nearby Massive Turbulent Disks and the Effects of Clump Clustering on Kiloparsec Scale Measurements of Clumps

We present $\sim$100 pc resolution Hubble Space Telescope H$α$ images of 10 galaxies from the DYnamics of Newly-Assembled Massive Objects (DYNAMO) survey of low-$z$ turbulent disk galaxies, and use these to undertake the first detailed systematic study of the effects of resolution and clump clustering on observations of clumps in turbulent disks. In the DYNAMO-{\em HST} sample we measure clump diameters spanning the range $d_{clump} \sim 100-800$~pc, and individual clump star formation rates as high as $\sim5$~M$_{\odot}$~yr$^{-1}$. DYNAMO clumps have very high SFR surface densities, $Σ_{SFR}\sim 15$~M$_{\odot}$~yr$^{-1}$~kpc$^{-2}$, $\sim100\times$ higher than in H{\sc ii} regions of nearby spirals. Indeed, SFR surface density provides a simple dividing line between massive star forming clumps and local star forming regions, where massive star forming clumps have $Σ_{SFR}> 0.5$~M$_{\odot}$~yr$^{-1}$~kpc$^{-2}$. When degraded to match the observations of galaxies in $z\sim 1-3$ surveys, DYNAMO galaxies are similar in morphology and measured clump properties to clumpy galaxies observed in the high-$z$ Universe. Emission peaks in the simulated high-redshift maps typically correspond to multiple clumps in full resolution images. This clustering of clumps systematically increases the apparent size and SFR of clumps in 1~kpc resolution maps, and decreases the measured SFR surface density of clumps by as much as a factor of 20$\times$. From these results we can infer that clump clustering is likely to strongly effect the measured properties of clumps in high-$z$ galaxies, which commonly have kiloparsec scale resolution.

preprint2016arXiv

Spurious Small-Scale Structure & Discreteness-Driven Relaxation in Cosmological Simulations

There is strong evidence that cosmological N-body simulations dominated by Warm Dark Matter (WDM) contain spurious or unphysical haloes, most readily apparent as regularly spaced low-mass haloes strung along filaments. We show that spurious haloes are a feature of traditional N-body simulations of cosmological structure formation models, including WDM and Cold Dark Matter (CDM) models, in which gravitational collapse proceeds in an initially anisotropic fashion, and arises naturally as a consequence of discreteness-driven relaxation. We demonstrate this using controlled N-body simulations of plane-symmetric collapse and show that spurious haloes are seeded at shell crossing by localised velocity perturbations induced by the discrete nature of the density field, and that their characteristic separation should be approximately the mean inter-particle separation of the N-body simulation, which is fixed by the mass resolution within the volume. Using cosmological N-body simulations in which particles are split into two collisionless components of fixed mass ratio, we find that the spatial distribution of the two components show signatures of discreteness-driven relaxation in their spatial distribution on both large and small scales. Adopting a spline kernel gravitational softening that is of order the comoving mean inter-particle separation helps to suppress the effect of discreteness-driven relaxation, but cannot eliminate it completely. These results provide further motivation for recent developments of new algorithms, which include, for example, revisions of the traditional N-body approach by means of spatially adaptive anistropric gravitational softenings or explicit solutions for the evolution of dark matter in phase space.

preprint2015arXiv

Connecting the Baryons: Multiwavelength Data for SKA HI Surveys

The science achievable with SKA HI surveys will be greatly increased through the combination of HI data with that at other wavelengths. These multiwavelength datasets will enable studies to move beyond an understanding of HI gas in isolation to instead understand HI as an integral part of the highly complex baryonic processes that drive galaxy evolution. As they evolve, galaxies experience a host of environmental and feedback influences, many of which can radically impact their gas content. Important processes include: accretion (hot and cold mode, mergers), depletion (star formation, galactic winds, AGN), phase changes (ionised/atomic/molecular), and environmental effects (ram pressure stripping, tidal effects, strangulation). Governing all of these to various extents is the underlying dark matter distribution. In turn, the result of these processes can significantly alter the baryonic states in which material is finally observed (stellar populations, dust, chemistry) and its morphology (galaxy type, bulge/disk ratio, bars, warps, radial profile). To fully understand the evolution of HI and the role it plays in galactic evolution requires the ability to quantify each of these separate processes, and hence to coordinate SKA HI surveys with extensive multi-band photometric and spectroscopic campaigns. In addition, multiwavelength data is essential for statistical methods of HI analysis such as HI stacking and intensity mapping cross-correlations. In this chapter, we examine some of the principal science motivations for acquiring multiwavelength data to match that from the extragalactic SKA HI surveys, and review the currently planned capacity to achieve this (eg. LSST, Euclid, W-FIRST, SPICA, ALMA, and 4MOST).

preprint2015arXiv

Extended Tully-Fisher Relations using HI Stacking

We present a new technique for the statistical evaluation of the Tully-Fisher relation (TFR) using spectral line stacking. This technique has the potential to extend TFR observations to lower masses and higher redshifts than possible through a galaxy-by-galaxy analysis. It further avoids the need for individual galaxy inclination measurements. To quantify the properties of stacked HI emission lines, we consider a simplistic model of galactic disks with analytically expressible line profiles. Using this model, we compare the widths of stacked profiles with those of individual galaxies. We then follow the same procedure using more realistic mock galaxies drawn from the S3-SAX model (a derivative of the Millennium simulation). Remarkably, when stacking the apparent HI lines of galaxies with similar absolute magnitude and random inclinations, the width of the stack is very similar to the width of the deprojected (= corrected for inclination) and dedispersed (= after removal of velocity dispersion) input lines. Therefore, the ratio between the widths of the stack and the deprojected/dedispersed input lines is approximately constant - about 0.93 - with very little dependence on the gas dispersion, galaxy mass, galaxy morphology, and shape of the rotation curve. Finally, we apply our technique to construct a stacked TFR using HIPASS data which already has a well defined TFR based on individual detections. We obtain a B-band TFR with a slope of $-8.5\pm0.4$ and a K-band relation with a slope of $-11.7\pm0.6$ for the HIPASS data set which is consistent with the existing results.

preprint2015arXiv

Galaxy And Mass Assembly (GAMA): Trends in galaxy colours, morphology, and stellar populations with large scale structure, group, and pair environments

We explore trends in galaxy properties with Mpc-scale structures using catalogues of environment and large scale structure from the Galaxy And Mass Assembly (GAMA) survey. Existing GAMA catalogues of large scale structure, group and pair membership allow us to construct galaxy stellar mass functions for different environmental types. To avoid simply extracting the known underlying correlations between galaxy properties and stellar mass, we create a mass matched sample of galaxies with stellar masses between $9.5 \leq \log{M_*/h^{-2} M_{\odot}} \leq 11$ for each environmental population. Using these samples, we show that mass normalised galaxies in different large scale environments have similar energy outputs, $u-r$ colours, luminosities, and morphologies. Extending our analysis to group and pair environments, we show galaxies that are not in groups or pairs exhibit similar characteristics to each other regardless of broader environment. For our mass controlled sample, we fail to see a strong dependence of Sérsic index or galaxy luminosity on halo mass, but do find that it correlates very strongly with colour. Repeating our analysis for galaxies that have not been mass controlled introduces and amplifies trends in the properties of galaxies in pairs, groups, and large scale structure, indicating that stellar mass is the most important predictor of the galaxy properties we examine, as opposed to environmental classifications.

preprint2014arXiv

100 deg$^2$ Mock Galaxy Cone for HI Surveys with the Early SKA

We distribute an easy-to-use mock catalog of galaxies with detailed neutral atomic hydrogen (HI) and auxiliary molecular and optical properties. The catalog covers a field of 10-by-10 degrees and a redshift range of z=0-1.2. It contains galaxies with 21cm peak flux densities down to 1uJy and is, within this flux limit, complete for HI masses above 10^8 solar masses. Five random realisations of the catalog in ASCII format (~4GB/file) and subtables with HI flux limits of 10u Jy (~500MB/file) and 100uJy$ (~30MB/file) can be downloaded at http://ict.icrar.org/store/staff/do/s3sax.

preprint2014arXiv

Cuboidal Dice and Gibbs Distributions

What are the face-probabilities of a cuboidal die, i.e. a die with different side-lengths? This paper introduces a model for these probabilities based on a Gibbs distribution. Experimental data produced in this work and drawn from the literature support the Gibbs model. The experiments also reveal that the physical conditions, such as the quality of the surface onto which the dice are dropped, can affect the face-probabilities. In the Gibbs model, those variations are condensed in a single parameter, adjustable to the physical conditions.

preprint2014arXiv

Extreme gas fractions in clumpy, turbulent disk galaxies at z~0.1

In this letter we report the discovery of CO fluxes, suggesting very high gas fractions in three disk galaxies seen in the nearby Universe (z ~ 0.1). These galaxies were investigated as part of the DYnamics of Newly Assembled Massive Objects (DYNAMO) survey. High-resolution Hubble Space Telescope imaging of these objects reveals the presence of large star forming clumps in the bodies of the galaxies, while spatially resolved spectroscopy of redshifted Halpha reveals the presence of high dispersion rotating disks. The internal dynamical state of these galaxies resembles that of disk systems seen at much higher redshifts (1 < z < 3). Using CO(1-0) observations made with the Plateau de Bure Interferometer, we find gas fractions of 20-30% and depletion times of tdep ~ 0.5 Gyr (assuming a Milky Way-like CO conversion factor). These properties are unlike those expected for low- redshift galaxies of comparable specific star formation rate, but they are normal for their high-z counterparts. DYNAMO galaxies break the degeneracy between gas fraction and redshift, and we show that the depletion time per specific star formation rate for galaxies is closely tied to gas fraction, independent of redshift. We also show that the gas dynamics of two of our local targets corresponds to those expected from unstable disks, again resembling the dynamics of high-z disks. These results provide evidence that DYNAMO galaxies are local analogues to the clumpy, turbulent disks, which are often found at high redshift.

preprint2014arXiv

Fundamental Mass-Spin-Morphology Relation of Spiral Galaxies

This work present high-precision measurements of the specific baryon angular momentum jb, contained in stars, atomic gas, and molecular gas, out to ~10 scale radii, in 16 nearby spiral galaxies of the THINGS sample. The accuracy of these measurements improves on existing studies by an order of magnitude, leading to the discovery of a strong correlation between the baryon mass Mb, jb, and the bulge mass fraction B/T, fitted by B/T=-(0.34+-0.03)*log(jb/Mb/[1e-7 kpc km/s/Msun])-(0.04+-0.01) on the full sample range of B/T=0.00-0.32 and Mb/Msun=1e9-1e11. The corresponding relation for the stellar quantities Ms and js is identical within the uncertainties. These M-j-B/T relations likely originate from the proportionality between j/M and the surface density of the disk that dictates its stability against (pseudo-)bulge formation. Using a CDM model, we can approximately explain classical scaling relations, such as the fundamental plane of spiral galaxies, the Tully-Fisher relation, and the mass-size relation, in terms of the M-j-B/T relation. These results advocate the use of mass and angular momentum as the most fundamental quantities of spiral galaxies.

preprint2014arXiv

Galaxy and Mass Assembly (GAMA): Fine filaments of galaxies detected within voids

Based on data from the Galaxy and Mass Assembly (GAMA) survey, we report on the discovery of structures that we refer to as `tendrils' of galaxies: coherent, thin chains of galaxies that are rooted in filaments and terminate in neighbouring filaments or voids. On average, tendrils contain 6 galaxies and span 10 $h^{-1}$ Mpc. We use the so-called line correlation function to prove that tendrils represent real structures rather than accidental alignments. We show that voids found in the SDSS-DR7 survey that overlap with GAMA regions contain a large number of galaxies, primarily belonging to tendrils. This implies that void sizes are strongly dependent on the number density and sensitivity limits of a survey. We caution that galaxies in low density regions, that may be defined as `void galaxies,' will have local galaxy number densities that depend on such observational limits and are likely higher than can be directly measured.

preprint2014arXiv

Surface wave dynamics in orbital shaken cylindrical containers

Be it to aerate a glass of wine before tasting, to accelerate a chemical reaction or to cultivate cells in suspension, the "swirling" (or orbital shaking) of a container ensures good mixing and gas exchange in an efficient and simple way. Despite being used in a large range of applications this intuitive motion is far from being understood and presents a richness of patterns and behaviors which has not yet been reported. The present research charts the evolution of the waves with the operating parameters identifying a large variety of patterns, ranging from single and multiple crested waves to breaking waves. Free surface and velocity fields measurements are compared to a potential sloshing model, highlighting the existence of various flow regimes. Our research assesses the importance of the modal response of the shaken liquids, laying the foundations for a rigorous mixing optimization of the orbital agitation in its applications. Copyright (2014) American Institute of Physics. This article may be downloaded for personal use only. Any other use requires prior permission of the author and the American Institute of Physics. The following article appeared in Physics of Fluids 26, 052104 (2014) and may be found at http://dx.doi.org/10.1063/1.4874612

preprint2014arXiv

Variation of galactic cold gas reservoirs with stellar mass

The stellar and neutral hydrogen (HI) mass functions at z~0 are fundamental benchmarks for current models of galaxy evolution. A natural extension of these benchmarks is the two-dimensional distribution of galaxies in the plane spanned by stellar and HI mass, which provides a more stringent test of simulations, as it requires the HI to be located in galaxies of the correct stellar mass. Combining HI data from the ALFALFA survey, with optical data from SDSS, we find a distinct envelope in the HI-to-stellar mass distribution, corresponding to an upper limit in the HI fraction that varies monotonically over five orders of magnitude in stellar mass. This upper envelope in HI fraction does not favour the existence of a significant population of dark galaxies with large amounts of gas but no corresponding stellar population. The envelope shows a break at a stellar mass of ~10^9 Msun, which is not reproduced by modern models of galaxy populations tracing both stellar and gas masses. The discrepancy between observations and models suggests a mass dependence in gas storage and consumption missing in current galaxy evolution prescriptions. The break coincides with the transition from galaxies with predominantly irregular morphology at low masses to regular disks at high masses, as well as the transition from cold to hot accretion of gas in simulations.

preprint2013arXiv

Confronting Cold Dark Matter Predictions with Observed Galaxy Rotations

The rich statistics of galaxy rotations as captured by the velocity function (VF) provides invaluable constraints on galactic baryon physics and the nature of dark matter (DM). However, the comparison of observed galaxy rotations against cosmological models is prone to subtle caveats that can easily lead to misinterpretations. Our analysis reveals full statistical consistency between ~5000 galaxy rotations, observed in line-of-sight projection, and predictions based on the standard cosmological model (LCDM) at the mass-resolution of the Millennium simulation (HI line-based circular velocities above ~50 km/s). Explicitly, the HI linewidths in the HI Parkes All Sky Survey (HIPASS) are found consistent with those in S3-SAX, a post-processed semi-analytic model for the Millennium simulation. Previously found anomalies in the VF can be plausibly attributed to (1) the mass-limit of the Millennium simulation, (2) confused sources in HIPASS, (3) inaccurate inclination measurements for optically faint sources, and (4) the non-detectability of gas-poor early-type galaxies. These issues can be bypassed by comparing observations and models using linewidth source counts rather than VFs. We investigate if and how well such source counts can constrain the temperature of DM.

preprint2013arXiv

Evolution of the Milky Way in Semi-Analytic Models: Detecting Cold Gas at z=3 with ALMA and SKA

We forecast the abilities of the Atacama Large Millimeter/submillimeter Array (ALMA) and the Square Kilometer Array (SKA) to detect CO and HI emission lines in galaxies at redshift z=3. A particular focus is set on Milky Way (MW) progenitors at z=3 for their detection within 24 h constitutes a key science goal of ALMA. The analysis relies on a semi-analytic model, which permits the construction of a MW progenitor sample by backtracking the cosmic history of all simulated present-day galaxies similar to the real MW. Results: (i) ALMA can best observe a MW at z=3 by looking at CO(3-2) emission. The probability of detecting a random model MW at 3-sigma in 24 h using 75 km/s channels is roughly 50%, and these odds can be increased by co-adding the CO(3-2) and CO(4-3) lines. These lines fall into ALMA band 3, which therefore represents the optimal choice towards MW detections at z=3. (ii) Higher CO transitions contained in the ALMA bands geq6 will be invisible, unless the considered MW progenitor coincidentally hosts a major starburst or an active black hole. (iii) The high-frequency array of SKA, fitted with 28.8 GHz receivers, would be a powerful instrument for observing CO(1-0) at z=3, able to detect nearly all simulated MWs in 24 h. (iv) HI detections in MWs at z=3 using the low-frequency array of SKA will be impossible in any reasonable observing time. (v) SKA will nonetheless be a supreme ha survey instrument through its enormous instantaneous field-of-view (FoV). A one year pointed HI survey with an assumed FoV of 410 sqdeg would reveal at least 10^5 galaxies at z=2.95-3.05. (vi) If the positions and redshifts of those galaxies are known from an optical/infrared spectroscopic survey, stacking allows the detection of HI at z=3 in less than 24 h.

preprint2013arXiv

Precise Tully-Fisher Relations without Galaxy Inclinations

Power-law relations between tracers of baryonic mass and rotational velocities of disk galaxies, so-called Tully-Fisher relations (TFRs), offer a wealth of applications in galaxy evolution and cosmology. However, measurements of rotational velocities require galaxy inclinations, which are difficult to measure, thus limiting the range of TFR studies. This work introduces a maximum likelihood estimation (MLE) method for recovering the TFR in galaxy samples with limited or no information on inclinations. The robustness and accuracy of this method is demonstrated using virtual and real galaxy samples. Intriguingly, the MLE reliably recovers the TFR of all test samples, even without using any inclination measurements - that is, assuming a sin(i)-distribution for galaxy inclinations. Explicitly, this 'inclination-free MLE' recovers the three TFR parameters (zero-point, slope, scatter) with statistical errors only about 1.5-times larger than the best estimates based on perfectly known galaxy inclinations with zero uncertainty. Thus, given realistic uncertainties, the inclination-free MLE is highly competitive. If inclination measurements have mean errors larger than 10 degrees, it is better not to use any inclinations, than to consider the inclination measurements to be exact. The inclination-free MLE opens interesting perspectives for future HI surveys by the SKA and its pathfinders.

preprint2013arXiv

The busy function: a new analytic function for describing the integrated 21-cm spectral profile of galaxies

Accurate parametrization of galaxies detected in the 21-cm HI emission is of fundamental importance to the measurement of commonly used indicators of galaxy evolution, including the Tully-Fisher relation and the HI mass function. Here, we propose a new analytic function, named the 'busy function', that can be used to accurately describe the characteristic double-horn HI profile of many galaxies. The busy function is a continuous, differentiable function that consists of only two basic functions, the error function, erf(x), and a polynomial, |x|^n, of degree n >= 2. We present the basic properties of the busy function and illustrate its great flexibility in fitting a wide range of HI profiles from the Gaussian profiles of dwarf galaxies to the broad, asymmetric double-horn profiles of spiral galaxies. Applications of the busy function include the accurate and efficient parametrization of observed HI spectra of galaxies and the construction of spectral templates for simulations and matched filtering algorithms. We demonstrate the busy function's power by automatically fitting it to the HI spectra of 1000 galaxies from the HIPASS Bright Galaxy Catalog, using our own C/C++ implementation, and comparing the resulting parameters with the catalogued ones. We also present two methods for determining the uncertainties of observational parameters derived from the fit.

preprint2013arXiv

The Quest for the Most Spherical Bubble

We describe a recently realized experiment producing the most spherical cavitation bubbles today. The bubbles grow inside a liquid from a point-plasma generated by a nanosecond laser pulse. Unlike in previous studies, the laser is focussed by a parabolic mirror, resulting in a plasma of unprecedented symmetry. The ensuing bubbles are sufficiently spherical that the hydrostatic pressure gradient caused by gravity becomes the dominant source of asymmetry in the collapse and rebound of the cavitation bubbles. To avoid this natural source of asymmetry, the whole experiment is therefore performed in microgravity conditions (ESA, 53rd and 56th parabolic flight campaign). Cavitation bubbles were observed in microgravity (~0g), where their collapse and rebound remain spherical, and in normal gravity (1g) to hyper-gravity (1.8g), where a gravity-driven jet appears. Here, we describe the experimental setup and technical results, and overview the science data. A selection of high-quality shadowgraphy movies and time-resolved pressure data is published online.

preprint2012arXiv

A Robust Measure of Cosmic Structure beyond the Power-Spectrum: Cosmic Filaments and the Temperature of Dark Matter

We discover that the mass of dark matter particles mDM is imprinted in phase-correlations of the cosmic density field more significantly than in the 2-point correlation. In particular, phase-correlations trace mDM out to scales about five times larger than the 2-point correlation. This result relies on a new estimator l(r) of pure phase-information in Fourier space, which can be interpreted as a parameter-free and scale-invariant tracer of filament-like structure. Based on simulated density fields we show how mDM can, in principle, be measured using l(r), given a suitably reconstructed density field.

preprint2012arXiv

Analytical Approximations for the Collapse of an Empty Spherical Bubble

The Rayleigh equation 3/2 R'+RR"+p/rho=0 with initial conditions R(0)=Rmax, R'(0)=0 models the collapse of an empty spherical bubble of radius R(T) in an ideal, infinite liquid with far-field pressure p and density rho. The solution for r=R/Rmax as a function of time t=T/Tcollapse, where R(Tcollapse)=0, is independent of Rmax, p, and rho. While no closed-form expression for r(t) is known we find that s(t)=(1-t^2)^(2/5) approximates r(t) with an error below 1%. A systematic development in orders of t^2 further yields the 0.001%-approximation r*(t)=s(t)[1-a Li(2.21,t^2)], where a=-0.01832099 is a constant and Li is the polylogarithm. The usefulness of these approximations is demonstrated by comparison to high-precision cavitation data obtained in microgravity.

preprint2012arXiv

Energy partition at the collapse of spherical cavitation bubbles

Spherically collapsing cavitation bubbles produce a shock wave followed by a rebound bubble. Here we present a systematic investigation of the energy partition between the rebound and the shock. Highly spherical cavitation bubbles are produced in microgravity, which suppress the buoyant pressure gradient that otherwise deteriorates the sphericity of the bubbles. We measure the radius of the rebound bubble and estimate the shock energy as a function of the initial bubble radius (2-5.6 mm) and the liquid pressure (10-80 kPa). Those measurements uncover a systematic pressure dependence of the energy partition between rebound and shock. We demonstrate that these observations agree with a physical model relying on a first-order approximation of the liquid compressibility and an adiabatic treatment of the non-condensable gas inside the bubble. Using this model we find that the energy partition between rebound and shock is dictated by a single non-dimensional parameter $ξ= Δpγ^6/[{p_{g0}}^{1/γ} (ρc^2)^{1-1/γ}]$, where $Δp=p_\infty-p_v$ is the driving pressure, $p_{\infty}$ is the static pressure in the liquid, $p_v$ is the vapor pressure, $p_{g0}$ is the pressure of the non-condensable gas at the maximal bubble radius, $γ$ is the adiabatic index of the non-condensable gas, $ρ$ is the liquid density, and $c$ is the speed of sound in the liquid.

preprint2011arXiv

"Oenodynamic": Hydrodynamic of wine swirling

A crucial step in wine tasting is the so called swirling, necessary to release the bouquet of the wine: a gentle circular movement of the glass generates a wave propagating along the glass walls, enhancing oxygenation and mixing. Although being used in a large variety of other applications (e.g. cells cultures in orbital shaken bioreactors) this motion is not yet well understood. In this fluid dynamics video we show the large variety of waves shapes generated by this simple movement, and we identify a group of dimensionless parameters governing the flow.

preprint2011arXiv

A Universal Scaling Law for Jets of Collapsing Bubbles

Cavitation bubbles collapsing and rebounding in a pressure gradient grad(p) form a "micro-jet" enveloped by a "vapor jet". This letter presents unprecedented observations of the vapor jets formed in a uniform gravity-induced grad(p), modulated aboard parabolic flights. The data uncovers that the normalized jet volume is independent of the liquid density and viscosity and proportional to zeta=grad(p)*R0/p, where R0 is the maximal bubble radius and p is the driving pressure. A derivation inspired by "Kelvin-Blake" considerations confirms this law and reveals its negligible dependence of surface tension. We further conjecture that the jet only pierces the bubble boundary if zeta>0.0004.

preprint2011arXiv

Bubble-Jets in Gravitational Fields (APS DFD Video V060)

We show visualizations of the gravity-induced jets formed by spherical bubbles collapsing in liquids subjected to normal gravity, micro-gravity, and hyper-gravity. These observations demonstrate that gravity can have a significant effect on cavitation bubbles. An analysis of the gravity-induced jets uncovers a scaling law between the size of bubble-induced jets and the non-dimensional parameter grad(p)*R0/p, where R0 is the maximal bubble radius and p is the driving pressure. This scaling law applies to any jet formed in any uniform pressure gradient grad(p).

preprint2011arXiv

Confined Shocks inside Isolated Liquid Volumes -- A New Path of Erosion?

The unique confinement of shock waves inside isolated liquid volumes amplifies the density of shock-liquid interactions. We investigate this universal principle through an interdisciplinary study of shock-induced cavitation inside liquid volumes, isolated in 2 and 3 dimensions. By combining high-speed visualizations of ideal water drops realized in microgravity with smoothed particle simulations we evidence strong shock-induced cavitation at the focus of the confined shocks. We extend this analysis to ground-observations of jets and drops using an analytic model, and argue that cavitation caused by trapped shocks offers a distinct mechanism of erosion in high-speed impacts (>100 m/s).

preprint2010arXiv

A Virtual Sky with Extragalactic HI and CO Lines for the SKA and ALMA

We present a sky simulation of the atomic HI emission line and the first ten CO rotational emission lines of molecular gas in galaxies beyond the Milky Way. The simulated sky field has a comoving diameter of 500/h Mpc, hence the actual field-of-view depends on the (user-defined) maximal redshift zmax; e.g. for zmax=10, the field of view yields ~4x4 sqdeg. For all galaxies, we estimate the line fluxes, line profiles, and angular sizes of the HI and CO emission lines. The galaxy sample is complete for galaxies with cold hydrogen masses above 10^8 Msun. This sky simulation builds on a semi-analytic model of the cosmic evolution of galaxies in a Lambda-cold dark matter (LCDM) cosmology. The evolving CDM-distribution was adopted from the Millennium Simulation, an N-body CDM-simulation in a cubic box with a side length of 500/h Mpc. This side length limits the coherence scale of our sky simulation: it is long enough to allow the extraction of the baryon acoustic oscillations (BAOs) in the galaxy power spectrum, yet the position and amplitude of the first acoustic peak will be imperfectly defined. This sky simulation is a tangible aid to the design and operation of future telescopes, such the SKA, the LMT, and ALMA. The results presented in this paper have been restricted to a graphical representation of the simulated sky and fundamental dN/dz-analyzes for peak flux density limited and total flux limited surveys of HI and CO. A key prediction is that HI will be harder to detect at redshifts z>2 than predicted by a no-evolution model. The future verification or falsification of this prediction will allow us to qualify the semi-analytic models.

preprint2010arXiv

Techniques for Generating Centimetric Drops in Microgravity and Application to Cavitation Studies

This paper describes the techniques and physical parameters used to produce stable centimetric water drops in microgravity, and to study single cavitation bubbles inside such drops (Parabolic Flight Campaigns, European Space Agency ESA). While the main scientific results have been presented in a previous paper, we shall herein provide the necessary technical background, with potential applications to other experiments. First, we present an original method to produce and capture large stable drops in microgravity. This technique succeeded in generating quasi-spherical water drops with volumes up to 8 ml, despite the residual g-jitter. We find that the equilibrium of the drops is essentially dictated by the ratio between the drop volume and the contact surface used to capture the drop, and formulate a simple stability criterion. In a second part, we present a setup for creating and studying single cavitation bubbles inside those drops. In addition, we analyze the influence of the bubble size and position on the drop behaviour after collapse, i.e. jets and surface perturbations.

preprint2009arXiv

Compactness of Cold Gas in High-Redshift Galaxies

Galaxies in the early Universe were more compact and contained more molecular gas than today. In this paper, we revisit the relation between these empirical findings, and we quantitatively predict the cosmic evolution of the surface densities of atomic (HI) and molecular (H2) hydrogen in regular galaxies. Our method uses a pressure-based model for the H2/HI-ratio of the Interstellar Medium, applied to ~3*10^7 virtual galaxies in the Millennium Simulation. We predict that, on average, the HI-surface density of these galaxies saturates at Sigma_HI<10 Msun/pc^2 at all redshifts (z), while H2-surface densities evolve dramatically as Sigma_H2(1+z)^2.4. This scaling is dominated by a (1+z)^2 surface brightness scaling originating from the (1+z)^-1 size scaling of galaxies at high z. Current measurements of Sigma_H2 at high z, derived from CO-observations, tend to have even higher values, which can be quantitatively explained by a selection bias towards merging systems. However, despite the consistency between our high-z predictions and the sparse empirical data, we emphasize that the empirical data potentially suffer from serious selection biases and that the semi-analytic models remain in many regards uncertain. As a case study, we investigate the cosmic evolution of simulated galaxies, which resemble the Milky Way at z=0. We explicitly predict their HI- and H2-distribution at z=1.5, corresponding to the CO-detected galaxy BzK-21000, and at z=3, corresponding to the primary science goal of the Atacama Large Millimeter/submillimeter Array (ALMA).