Cosmological Simulations Of Galaxy Formation Research Articles

A new model for including galactic winds in simulations of galaxy formation – I. Introducing the Physically Evolved Winds (PhEW) model

ABSTRACT The propagation and evolution of cold galactic winds in galactic haloes is crucial to galaxy formation models. However, modelling of this process in hydrodynamic simulations of galaxy formation is oversimplified owing to a lack of numerical resolution and often neglects critical physical processes such as hydrodynamic instabilities and thermal conduction. We propose an analytic model, Physically Evolved Winds, that calculates the evolution of individual clouds moving supersonically through a uniform ambient medium. Our model reproduces predictions from very high resolution cloud-crushing simulations that include isotropic thermal conduction over a wide range of physical conditions. We discuss the implementation of this model into cosmological hydrodynamic simulations of galaxy formation as a subgrid prescription to model galactic winds more robustly both physically and numerically.

The AREPO Public Code Release

We introduce the public version of the cosmological magnetohydrodynamical moving-mesh simulation code Arepo. This version contains a finite-volume magnetohydrodynamics algorithm on an unstructured, dynamic Voronoi tessellation coupled to a tree-particle-mesh algorithm for the Poisson equation either on a Newtonian or cosmologically expanding spacetime. Time-integration is performed adopting local timestep constraints for each cell individually, solving the fluxes only across active interfaces, and calculating gravitational forces only between active particles, using an operator-splitting approach. This allows simulations with high dynamic range to be performed efficiently. Arepo is a massively distributed-memory parallel code, using the Message Passing Interface (MPI) communication standard and employing a dynamical work-load and memory balancing scheme to allow optimal use of multi-node parallel computers. The employed parallelization algorithms of Arepo are deterministic and produce binary-identical results when re-run on the same machine and with the same number of MPI ranks. A simple primordial cooling and star formation model is included as an example of sub-resolution models commonly used in simulations of galaxy formation. Arepo also contains a suite of computationally inexpensive test problems, ranging from idealized tests for automated code verification to scaled-down versions of cosmological galaxy formation simulations, and is extensively documented in order to assist adoption of the code by new scientific users.

Cosmic variance of z > 7 galaxies: prediction from bluetides

ABSTRACT In the coming decade, a new generation of telescopes, including JWST and WFIRST, will probe the period of the formation of first galaxies and quasars, and open up the last frontier for structure formation. Recent simulations and observations have suggested that these galaxies are strongly clustered (with large-scale bias ≳6), and therefore have significant cosmic variance. In this work, we use bluetides, the largest volume cosmological simulation of galaxy formation, to directly estimate the cosmic variance for current and upcoming surveys. Given its resolution and volume, bluetides can probe the bias and cosmic variance of z > 7 galaxies between magnitude MUV ∼ −16 and MUV ∼ −22 over survey areas ∼0.1 arcmin2 to ∼10 deg2. Within this regime, the cosmic variance decreases with survey area/ volume as a power law with exponents between ∼−0.25 and ∼−0.45. For the planned 10 deg2 field of WFIRST, the cosmic variance is between $3{{\ \rm per\ cent}}$ and $10{{\ \rm per\ cent}}$. Upcoming JWST medium/ deep surveys with areas up to A ∼ 100 arcmin2 will have cosmic variance ranging from ${\sim}20\,\mathrm{ to}\,50{{\ \rm per\ cent}}$. Lensed surveys have the highest cosmic variance ${\gtrsim}40{{\ \rm per\ cent}}$; the cosmic variance of MUV ≲ −16 galaxies is ${\lesssim}100{{\ \rm per\ cent}}$ up to z ∼ 11. At higher redshifts such as z ∼ 12 (14), effective volumes of ≳ (8 Mpc h−1)3 (≳(12 Mpc h−1)3) are required to limit the cosmic variance to within $100{{\ \rm per\ cent}}$. Finally, we find that cosmic variance is larger than Poisson variance and forms the dominant component of the overall uncertainty in all current and upcoming surveys. We present our calculations in the form of simple fitting functions and an online cosmic variance calculator (CV_AT_COSMIC_DAWN) that we publicly release.

GW190412 as a Third-generation Black Hole Merger from a Super Star Cluster

Abstract We explore the possibility that GW190412, a binary black hole merger with a non-equal-mass ratio and significantly spinning primary, was formed through repeated black hole mergers in a dense super star cluster. Using a combination of semianalytic prescriptions for the remnant spin and recoil kick of black hole mergers, we show that the mass ratio and spin of GW190412 are consistent with a binary black hole whose primary component has undergone two successive mergers from a population of black holes in a high-metallicity environment. We then explore the production of GW190412-like analogs in the CMC Cluster Catalog, a grid of 148 N-body star cluster models, as well as a new model, behemoth, with nearly 107 particles and initial conditions taken from a cosmological MHD simulation of galaxy formation. We show that, if the spins of black holes born from stars are small, the production of binaries with GW190412-like masses and spins is dominated by massive super star clusters with high metallicities and large central escape speeds. While many are observed in the local universe, our results suggest that a careful treatment of these massive clusters, many of which may have been disrupted before the present day, is necessary to characterize the production of unique gravitational-wave events produced through dynamics.

QSO obscuration at high redshift (z ≳ 7): predictions from the bluetides simulation

ABSTRACT High-$z$ AGNs hosted in gas-rich galaxies are expected to grow through significantly obscured accretion phases. This may limit or bias their observability. In this work, we use bluetides, a large volume cosmological simulation of galaxy formation to examine quasar obscuration for the highest redshift ($z$ ≥ 7) supermassive black holes residing in the centre of galaxies. We find that for the bright quasars, most of the high-column density gas ($\rm {\gt} 90 {\rm {per\ cent}}$) resides in the innermost regions of the host galaxy (typically within <10 ckpc), while the gas in the outskirts is a minor contributor to the NH. The brightest quasars can have large angular variations in galactic obscuration, over 2 orders of magnitude (ranging from column density $N_\mathrm{H} \sim 10^{21.5 \!-\! 24}\, \rm {cm}^{-2}$), where the lines of sight with the lowest obscuration are those formed via strong gas outflows driven by AGN feedback. The obscured fraction P(NH > 1023 cm−2) typically ranges from 0.6 to 1.0 for increasing LX (with $L_\mathrm{ X} \gt 10^{43} \, \rm {erg\, s}^{-1}$), with no clear trend of redshift evolution. Due to the angular variation in NH, all relations between NH and LX, MBH, and galaxy host properties (global M*, $M_{\rm H_2}$, and star formation rate) show appreciable scatter. The dust optical depth in the UV band τUV has tight positive correlation with NH. Our dust-extincted UV luminosity function (UVLF) is about 1.5 dex lower than the intrinsic UVLF, implying that more than 99 per cent of the $z$ ∼ 7 AGNs are heavily dust extincted and therefore would be missed by the UV-band observation.

Origin of giant stellar clumps in high-redshift galaxies

ABSTRACT We examine the nature of kpc-scale clumps seen in high-redshift galaxies using a suite of cosmological simulations of galaxy formation. We identify rest-frame UV clumps in mock HST images smoothed to 500 pc resolution, and compare them with the intrinsic 3D clumps of young stars identified in the simulations with 100 pc resolution. According to this comparison for the progenitors of Milky Way-sized galaxies probed by our simulations, we expect that the stellar masses of the observed clumps are overestimated by as much as an order of magnitude, and that the sizes of these clumps are also overestimated by factor of several, due to a combination of spatial resolution and projection. The masses of young stars contributing most of the UV emission can also be overestimated by factor of a few. We find that most clumps of young stars present in a simulation at one time dissolve on a timescale shorter than ∼150 Myr. Some clumps with dense cores can last longer but eventually disperse. Most of the clumps are not bound structures, with virial parameter αvir > 1. We find similar results for clumps identified in mock maps of H α emission measure. We examine the predictions for effective clump sizes from the linear theory of gravitational perturbations and conclude that they are inconsistent with being formed by global disc instabilities. Instead, the observed clumps represent random projections of multiple compact star-forming regions.

The Milky Way’s stellar streams and globular clusters do not align in a Vast Polar Structure

ABSTRACT There is increasing evidence that a substantial fraction of Milky Way satellite galaxies align in a rotationally supported plane of satellites, a rare configuration in cosmological simulations of galaxy formation. It has been suggested that other Milky Way substructures (namely young halo globular clusters and stellar/gaseous streams) similarly tend to align with this plane, accordingly dubbed the Vast Polar Structure (VPOS). Using systemic proper motions inferred from Gaia data, we find that globular cluster orbital poles are not clustered in the VPOS direction, though the population with the highest VPOS membership fraction is the young halo clusters (∼30 per cent). We additionally provide a current census of stellar streams, including new streams discovered using the Dark Energy Survey and Gaia data sets, and find that stellar stream normals are also not clustered in the direction of the VPOS normal. We also find that, based on orbit modelling, there is a likely association between NGC 3201 and the Gjöll stellar stream and that, based on its orbital pole, NGC 4147 is likely not a Sagittarius globular cluster. That the Milky Way’s accreted globular clusters and streams do not align in the same planar configuration as its satellites suggests that the plane of satellites is either a particularly stable orbital configuration or a population of recently accreted satellites. Neither of these explanations is particularly likely in light of other recent studies, leaving the plane of satellites problem as one of the more consequential open problems in galaxy formation and cosmology.

R-process enrichment in ultrafaint dwarf galaxies

ABSTRACT We study the enrichment and mixing of r-process elements in ultrafaint dwarf galaxies (UFDs). We assume that r-process elements are produced by neutron-star mergers (NSMs), and examine multiple models with different natal kick velocities and explosion energies. To this end, we perform cosmological simulations of galaxy formation to follow mixing of the dispersed r-process elements driven by star formation and the associated stellar feedback in progenitors of UFDs. We show that the observed europium abundance in Reticulum II is reproduced by our inner explosion model where an NSM is triggered at the centre of the galaxy, whereas the relatively low abundance in Tucana III is reproduced if an NSM occurs near the virial radius of the progenitor galaxy. The latter case is realized only if the neutron-star binary has a large natal kick velocity and travels over a long distance of a kiloparsec before merger. In both the inner and outer explosion cases, it is necessary for the progenitor galaxy to sustain prolonged star formation over a few hundred million years after the NSM, so that the dispersed r-process elements are well mixed within the interstellar medium. Short-duration star formation results in inefficient mixing, and then a large variation is imprinted in the stellar europium abundances, which is inconsistent with the observations of Reticulum II and Tucana III.

Cosmological simulations of galaxy formation

Over recent decades, cosmological simulations of galaxy formation have been instrumental in advancing our understanding of structure and galaxy formation in the Universe. These simulations follow the nonlinear evolution of galaxies, modelling a variety of physical processes over an enormous range of time and length scales. A better understanding of the relevant physical processes, improved numerical methods and increased computing power have led to simulations that can reproduce a large number of the observed galaxy properties. Modern simulations model dark matter, dark energy and ordinary matter in an expanding space-time starting from well-defined initial conditions. The modelling of ordinary matter is most challenging due to the large array of physical processes affecting this component. Cosmological simulations have also proven useful to study alternative cosmological models and their impact on the galaxy population. This Technical Review presents a concise overview of the methodology of cosmological simulations of galaxy formation and their different applications. Cosmological computer simulations of galaxy formation emerged as the primary tool to study structure formation in the Universe. This Technical Review describes the main techniques and ingredients of such simulations and their application to develop and constrain galaxy formation theories.

The dust-to-gas and dust-to-metal ratio in galaxies from z = 0 to 6

ABSTRACT We present predictions for the evolution of the galaxy dust-to-gas ratio (DGR) and dust-to-metal ratio (DTM) from z = 0 → 6, using a model for the production, growth, and destruction of dust grains implemented into the simba cosmological hydrodynamic galaxy formation simulation. In our model, dust forms in stellar ejecta, grows by the accretion of metals, and is destroyed by thermal sputtering and supernovae. Our simulation reproduces the observed dust mass function at z = 0, but modestly underpredicts the mass function by ∼×3 at z ∼ 1–2. The z = 0 DGR versus metallicity relationship shows a tight positive correlation for star-forming galaxies, while it is uncorrelated for quenched systems. There is little evolution in the DGR–metallicity relationship between z = 0 and 6. We use machine learning techniques to search for the galaxy physical properties that best correlate with the DGR and DTM. We find that the DGR is primarily correlated with the gas-phase metallicity, though correlations with the depletion time-scale, stellar mass, and gas fraction are non-negligible. We provide a crude fitting relationship for DGR and DTM versus the gas-phase metallicity, along with a public code package that estimates the DGR and DTM given a set of galaxy physical properties.

Quantifying inhomogeneities in the HI distributions of simulated galaxies

The NIHAO cosmological simulations form a collection of a hundred high-resolution galaxies. We used these simulations to test the impact of stellar feedback on the morphology of the HI distribution in galaxies. We ran a subsample of twenty of the galaxies with different parameterizations of stellar feedback, looking for differences in the HI spatial distribution and morphology. We found that different feedback models do leave a signature in HI, and can potentially be compared with current and future observations. These findings can help inform future modeling efforts in the parameterization of stellar feedback in cosmological simulations of galaxy formation and evolution.

Mergers, starbursts, and quenching in the simba simulation

ABSTRACT We use the simba cosmological galaxy formation simulation to investigate the relationship between major mergers ($\lesssim$4:1), starbursts, and galaxy quenching. Mergers are identified via sudden jumps in stellar mass M* well above that expected from in situ star formation, while quenching is defined as going from specific star formation rate (sSFR) $\gt t_{\rm H}^{-1}$ to $\lt 0.2t_{\rm H}^{-1}$, where tH is the Hubble time. At z ≈ 0–3, mergers show ∼2–3× higher SFR than a mass-matched sample of star-forming galaxies, but globally represent $\lesssim 1{{\ \rm per\ cent}}$ of the cosmic SF budget. At low masses, the increase in SFR in mergers is mostly attributed to an increase in the H2 content, but for $M_*\gtrsim 10^{10.5} \,\mathrm{ M}_{\odot }$ mergers also show an elevated star formation efficiency suggesting denser gas within merging galaxies. The merger rate for star-forming galaxies shows a rapid increase with redshift, ∝(1 + z)3.5, but the quenching rate evolves much more slowly, ∝(1 + z)0.9; there are insufficient mergers to explain the quenching rate at $z\lesssim 1.5$. simba first quenches galaxies at $z\gtrsim 3$, with a number density in good agreement with observations. The quenching time-scales τq are strongly bimodal, with ‘slow’ quenchings (τq ∼ 0.1tH) dominating overall, but ‘fast’ quenchings (τq ∼ 0.01tH) dominating in M* ∼ 1010–1010.5 M$\odot$ galaxies, likely induced by simba’s jet-mode black hole feedback. The delay time distribution between mergers and quenching events suggests no physical connection to either fast or slow quenching. Hence, simba predicts that major mergers induce starbursts, but are unrelated to quenching in either fast or slow mode.

Predictions for the abundance of high-redshift galaxies in a fuzzy dark matter universe

ABSTRACT During the last decades, rapid progress has been made in measurements of the rest-frame ultraviolet (UV) luminosity function (LF) for high-redshift galaxies (z ≥ 6). The faint-end of the galaxy LF at these redshifts provides powerful constraints on different dark matter (DM) models that suppress small-scale structure formation. In this work we perform full hydrodynamical cosmological simulations of galaxy formation using an alternative DM model composed of extremely light bosonic particles (m ∼ 10−22 eV), also known as fuzzy dark matter (FDM), and examine the predictions for the galaxy stellar mass function and LF at z ≥ 6 for a range of FDM masses. We find that for FDM models with bosonic mass m = 5 × 10−22 eV, the number density of galaxies with stellar mass $\rm M_* \sim 10^7 M_{\odot }$ is suppressed by $\sim 40\, {\rm per\, cent}$ at z = 9, $\sim 20\, {\rm per\, cent}$ at z = 5, and the UV LFs within magnitude range of −16 < MUV < −14 is suppressed by $\sim 60\, {\rm per\, cent}$ at z = 9, $\sim 20\, {\rm per\, cent}$ at z = 5 comparing to the cold dark matter counterpart simulation. Comparing our predictions with current measurements of the faint-end LFs (−18 ≤ MUV ≤ −14), we find that FDM models with m22 < 5 × 10−22 are ruled out at 3σ confidence level. We expect that future LF measurements by James Webb Space Telescope, which will extend down to MUV ∼ −13 for z ≲ 10, with a survey volume that is comparable to the Hubble Ultra Deep Field would have the capability to constrain FDM models to m ≳ 10−21 eV.

Uncovering the birth of the Milky Way through accurate stellar ages with Gaia

Knowledge of ages for stars formed over a galaxy's lifetime is fundamental to understand its formation and evolution. However, stellar ages are difficult to obtain since they cannot be measured from observations, being comparison with stellar models (Soderblom 2010) required. Alternatively, age distributions can be derived applying the robust technique of colour-magnitude diagram fitting (Gallart et al. 2005), till now mainly employed to study nearby galaxies. The new distances to individual Milky Way stars from the Gaia mission (Brown et al. 2018) have allowed us to use this technique to derive ages from a thick disk colour-magnitude diagram, and from the enigmatic, two-sequenced colour-magnitude diagram of the kinematically hot local halo (Babusiaux et al. 2018), which blue-sequence has been linked to a major accretion event (Haywood et al. 2018, Helmi et al. 2018). Because accurate ages were lacking, the time of the merger and its role on our Galaxy's early evolution remained unclear. We show that the stars in both halo sequences share identical age distributions, and are older than the bulk of thick disc stars. The sharp halo age cut 10 Gyr ago can be identified with the accretion of Gaia-Enceladus. Along with state-of-the-art cosmological simulations of galaxy formation (Brook et al. 2012), these robust ages allow us to order the early sequence of events that shaped our Galaxy, identifying the red-sequence as the first stars formed within the Milky Way progenitor which, because of their kinematics, can be described as its long sought in-situ halo.

The Fundamental Metallicity Relation Emerges from the Local Anti-correlation between Star Formation Rate and Gas-phase Metallicity that Exists in Disk Galaxies

Abstract The fundamental metallicity relation (FMR) states that galaxies of the same stellar mass but larger star formation rate (SFR) tend to have smaller gas-phase metallicity ( ). It is thought to be fundamental because it naturally arises from the stochastic feeding of star formation from external metal-poor gas accretion, a process extremely elusive to observe but essential according the cosmological simulations of galaxy formation. In this Letter, we show how the FMR emerges from the local anti-correlation between SFR surface density and Z g recently observed to exist in disk galaxies. We analytically derive the global FMR from the local law, and then show that both relations agree quantitatively when considering the star-forming galaxies of the Mapping Nearby Galaxies at APO (MaNGA) survey. Thus, understanding the FMR becomes equivalent to understanding the origin of the anti-correlation between SFR and metallicity followed by the set of star-forming regions of any typical galaxy. The correspondence between local and global laws is not specific to the FMR, so that a number of local relations should exist that are associated with known global relations.

Fast and inefficient star formation due to short-lived molecular clouds and rapid feedback.

The physics of star formation and the deposition of mass, momentum, and energy into the interstellar medium by massive stars (‘feedback’) are the main uncertainties in modern cosmological simulations of galaxy formation and evolution1, 2. These processes determine the properties of galaxies3, 4, but are poorly understood on the ≲100 pc scale of individual giant molecular clouds (GMCs)5, 6 resolved in modern galaxy formation simulations7, 8. The key question is why the timescale for depleting molecular gas through star formation in galaxies (tdep ≈ 2 Gyr)9, 10 exceeds the dynamical timescale of GMCs by two orders of magnitude11. Either most of a GMC’s mass is converted into stars over many dynamical times12, or only a small fraction turns into stars before the GMC is dispersed on a dynamical timescale13, 14. Here we report our observation that molecular gas and star formation are spatially decorrelated on GMC scales in the nearby flocculent spiral galaxy NGC300, contrary to their tight correlation on galactic scales5. We demonstrate that this de-correlation implies rapid evolutionary cycling between GMCs, star formation, and feedback. We apply a novel statistical method15, 16 to quantify the evolutionary timeline and find that star formation is regulated by efficient stellar feedback, driving GMC dispersal on short timescales (~1.5 Myr) due to radiation and stellar winds, prior to supernova explosions. This feedback limits GMC lifetimes to about one dynamical timescale (~10 Myr), with integrated star formation efficiencies of only 2–3%. Our findings reveal that galaxies consist of building blocks undergoing vigorous, feedback-driven lifecycles, that vary with the galactic environment and collectively define how galaxies form stars. Systematic applications of this multi-scale analysis to large galaxy samples will provide key input for a predictive, bottom-up theory of galaxy formation and evolution.

An observational test for star formation prescriptions in cosmological hydrodynamical simulations

State-of-the-art cosmological hydrodynamical simulations of galaxy formation have reached the point at which their outcomes result in galaxies with ever more realism. Still, the employed sub-grid models include several free parameters such as the density threshold, $n$, to localize the star-forming gas. In this work, we investigate the possibilities to utilize the observed clustered nature of star formation (SF) in order to refine SF prescriptions and constrain the density threshold parameter. To this end, we measure the clustering strength, correlation length and power-law index of the two-point correlation function of young ($\tau<50$ Myr) stellar particles and compare our results to observations from the HST Legacy Extragalactic UV Survey (LEGUS). Our simulations reveal a clear trend of larger clustering signal and power-law index and lower correlation length as the SF threshold increases with only mild dependence on galaxy properties such as stellar mass or specific star formation rate. In conclusion, we find that the observed clustering of SF is inconsistent with a low threshold for SF ($n<1$ cm$^{-3}$) and strongly favours a high value for the density threshold of SF ($n>10$ cm$^{-3}$), as for example employed in the NIHAO project.

Simba: Cosmological simulations with black hole growth and feedback

We introduce the Simba simulations, the next generation of the Mufasa cosmological galaxy formation simulations run with Gizmo's meshless finite mass hydrodynamics. Simba includes updates to Mufasa's sub-resolution star formation and feedback prescriptions, and introduces black hole growth via the torque-limited accretion model of Angl\'es-Alc\'azar et al. (2017) from cold gas and Bondi accretion from hot gas, along with black hole feedback via kinetic bipolar outflows and X-ray energy. Ejection velocities are taken to be ~10^3 km/s at high Eddington ratios, increasing to ~8000 km/s at Eddington ratios below 2%, with a constant momentum input of 20L/c. Simba further includes an on-the-fly dust production, growth, and destruction model. Our Simba run with (100 Mpc/h)^3 and 1024^3 gas elements reproduces numerous observables, including galaxy stellar mass functions at z=0-6, the stellar mass--star formation rate main sequence, HI and H2 fractions, the mass-metallicity relation at z=0 and z=2, star-forming galaxy sizes, hot gas fractions in massive halos, and z=0 galaxy dust properties. However, Simba also yields an insufficiently sharp truncation of the z=0 mass function, and too-large sizes for low-mass quenched galaxies. We show that Simba's jet feedback is primarily responsible for quenching massive galaxies.

NIHAO XX: the impact of the star formation threshold on the cusp–core transformation of cold dark matter haloes

We use cosmological hydrodynamical galaxy formation simulations from the NIHAO project to investigate the impact of the threshold for star formation on the response of the dark matter (DM) halo to baryonic processes. The fiducial NIHAO threshold, $n=10\, {\rm cm}^{-3}$, results in strong expansion of the DM halo in galaxies with stellar masses in the range $10^{7.5} < M_{star} < 10^{9.5} M_{\odot}$. We find that lower thresholds such as $n=0.1$ (as employed by the EAGLE/APOSTLE and Illustris/AURIGA projects) do not result in significant halo expansion at any mass scale. Halo expansion driven by supernova feedback requires significant fluctuations in the local gas fraction on sub-dynamical times (i.e., < 50 Myr at galaxy half-light radii), which are themselves caused by variability in the star formation rate. At one per cent of the virial radius, simulations with $n=10$ have gas fractions of $\simeq 0.2$ and variations of $\simeq 0.1$, while $n=0.1$ simulations have order of magnitude lower gas fractions and hence do not expand the halo. The observed DM circular velocities of nearby dwarf galaxies are inconsistent with CDM simulations with $n=0.1$ and $n=1$, but in reasonable agreement with $n=10$. Star formation rates are more variable for higher $n$, lower galaxy masses, and when star formation is measured on shorter time scales. For example, simulations with $n=10$ have up to 0.4 dex higher scatter in specific star formation rates than simulations with $n=0.1$. Thus observationally constraining the sub-grid model for star formation, and hence the nature of DM, should be possible in the near future.

The effect of rotation on the thermal instability of stratified galactic atmospheres – II. The formation of high-velocity clouds

Whether High Velocity Clouds (HVCs) can form by condensation of the hot ($T \sim 10^6 \, {\rm K}$) Galactic corona as a consequence of thermal instabilities has been controversial. Here we re-examine this problem and we suggest that rotation of the corona might be a missing key ingredient. We do this by studying the evolution of the models of rotating galactic coronae presented in Sormani et al. (2018) under the presence of cooling and thermal conduction. We combine a linear stability analysis with the results of local and global hydrodynamical simulations. We find that condensations are likely to occur in regions where the corona has substantial rotational support. Under reasonably general assumptions on the rotation profile of the corona, the locations where condensations are expected are in remarkable agreement with the observed location of the major non-magellanic HVCs complexes in our Galaxy (namely, at distances $< 15 \, \rm kpc$ from the Sun and within $30^\circ$ from the disc plane). We conclude that HVCs can form by thermal instabilities provided that (i) the corona is rotating substantially in the inner ($R < 50 \, \rm kpc$) parts, as suggested by current observational data and predicted by cosmological simulations of galaxy formation; (ii) close to the disc the corona is well-represented by a nearly-equilibrium stratified rotating structure (as opposed to a fast cooling flow). Our results also suggest that a better understanding of the disc-halo interface, including supernova feedback, is critical to understand the origin of HVCs.