Search

We explore the properties of Milky Way subhalos in self-interacting dark matter models for moderate cross sections of 1 to 5 cm$^2$g$^{-1}$ using high-resolution zoom-in N-body simulations. We include the gravitational potential of a baryonic disk and bulge matched to the Milky Way, which is critical for getting accurate predictions. The predicted number and distribution of subhalos within the host halo are similar for 1 and 5 cm$^2$g$^{-1}$ models, and they agree with observations of Milky Way satellite galaxies only if subhalos with peak circular velocity over all time > 4.5 km/s are able to form galaxies. We do not find distinctive signatures in the pericenter distribution of the subhalos that could help distinguish the models. Using an analytic model to extend the simulation results, we are able to show that subhalos in models with cross sections between 1 and 5 cm$^2$g$^{-1}$ are not dense enough to match the densest ultra-faint and classical dwarf spheroidal galaxies in the Milky Way. This motivates velocity-dependent cross sections with values larger than 5 cm$^2$g$^{-1}$ at the velocities relevant for the satellites such that core collapse would occur in some of the ultra-faint and classical dwarf spheroidals., Version accepted for publication by MNRAS. 19 pages, 16 figures, 2 tables, 5 appendices

Self-interacting dark matter (SIDM) models have received great attention over the past decade as solutions to the small-scale puzzles of astrophysics. Though there are different implementations of dark matter (DM) self-interactions in N-body codes of structure formation, there has not been a systematic study to compare the predictions of these different implementations. We investigate the implementation of dark matter self-interactions in two simulation codes:gizmo and arepo. We begin with identical initial conditions for an isolated 1010 M⊙ dark matter halo and investigate the evolution of the density and velocity dispersion profiles in gizmo and arepo for SIDM cross-section over mass of 1, 5, and 50 $\rm cm^2\, g^{-1}$. Our tests are restricted to the core expansion phase, where the core density decreases and core radius increases with time. We find better than 30 per cent agreement between the codes for the density profile in this phase of evolution, with the agreement improving at higher resolution. We find that varying code-specific SIDM parameters changes the central halo density by less than 10 per cent outside of the convergence radius. We argue that SIDM core formation is robust across the two different schemes and conclude that these codes can reliably differentiate between cross-sections of 1, 5, and 50 $\rm cm^2\, g^{-1}$, but finer distinctions would require further investigation.

We present a suite of baryonic cosmological zoom-in simulations of self-interacting dark matter (SIDM) haloes within the ``Feedback In Realistic Environment'' (FIRE) project. The three simulated haloes have virial masses of $\sim 10^{12}\, \text{M}_\odot$ at $z=0$, and we study velocity-independent self-interaction cross sections of 1 and 10 ${\rm cm^2 \, g^{-1}}$. We study star formation rates and the shape of dark matter density profiles of the parent haloes in both cold dark matter (CDM) and SIDM models. Galaxies formed in the SIDM haloes have higher star formation rates at $z\leq1$, resulting in more massive galaxies compared to the CDM simulations. While both CDM and SIDM simulations show diverse shape of the dark matter density profiles, the SIDM haloes can reach higher and more steep central densities within few kpcs compared to the CDM haloes. We identify a correlation between the build-up of the stars within the half-mass radii of the galaxies and the growth in the central dark matter densities. The thermalization process in the SIDM haloes is enhanced in the presence of a dense stellar component. Hence, SIDM haloes with highly concentrated baryonic profiles are predicted to have higher central dark matter densities than the CDM haloes. Overall, the SIDM haloes are more responsive to the presence of a massive baryonic distribution than their CDM counterparts., 10 pages, 5 figures. Accepted for publication in MNRAS

Fuzzy Dark Matter (FDM), consisting of ultralight bosons ($m_{\rm b} \sim 10^{-22}\ \rm eV$), is an intriguing alternative to Cold Dark Matter. Numerical simulations that solve the Schr\"odinger-Poisson (SP) equation show that FDM halos consist of a central solitonic core, which is the ground state of the SP equation, surrounded by an envelope of interfering excited states. These excited states also interfere with the soliton, causing it to oscillate and execute a confined random walk with respect to the halo center of mass. Using high-resolution numerical simulations of a $6.6 \times 10^9 M_{\odot}$ FDM halo with $m_{\rm b} = 8 \times 10^{-23}\ \rm eV$ in isolation, we demonstrate that the wobbling, oscillating soliton gravitationally perturbs nuclear objects, such as supermassive black holes or dense star clusters, causing them to diffuse outwards. In particular, we show that, on average, objects with mass $\lesssim 0.3 \%$ of the soliton mass ($M_{\rm sol}$) are expelled from the soliton in $\sim 3\ \rm Gyr$, after which they continue their outward diffusion due to gravitational interactions with the soliton and the halo granules. More massive objects ($\gtrsim 1 \% M_{\rm sol}$), while executing a random walk, remain largely confined to the soliton due to dynamical friction. We also present an effective treatment of the diffusion, based on kinetic theory, that accurately reproduces the outward motion of low mass objects and briefly discuss how the observed displacements of star clusters and active galactic nuclei from the centers of their host galaxies can be used to constrain FDM., Comment: 20 pages, 9 figures, minor revisions, accepted for publication in ApJ

In recent years, the scalar field dark matter (SFDM), also called ultralight bosonic dark matter, has received considerable attention due to the number of problems it might help to solve. Among these are the cusp-core problem and the abundance of small structures of the standard cold dark matter (CDM) model. In this paper we show that multi-state solutions of the low energy and weak gravitational field limit of field equations, interpreted as galactic halo density profiles, can provide a possible explanation to the anisotropic distribution of satellite galaxies observed in the Milky Way, M31 and Centaurus A, where satellites trajectories seem to concentrate on planes close to the poles of the galaxies instead of following homogeneously distributed trajectories. The core hypothesis is that multi-state solutions of the equations describing the dynamics of this dark matter candidate, namely, the Gross-Pitaevskii-Poisson equations, with monopolar and dipolar contributions, can possibly explain the anisotropy of satellite trajectories. In order to construct a proof of concept, we study the trajectories of a number of test particles traveling on top of the gravitational potential due to a multi-state halo with modes (1,0,0)+(2,1,0). The result is that particles accumulate asymptotically in time on planes passing close to the poles. Satellite galaxies are not test particles but interpreted as such, our results indicate that in the asymptotic time their trajectories do not distribute isotropically, instead they prefer to have orbital poles accumulating near the equatorial plane of the multistate halo. The concentration of orbital poles depends on whether the potential is monopolar or dipolar dominated., This version has been accepted for publication in Physical Review D

Using the phat-ELVIS suite of Milky Way-size halo simulations, we show that subhalo orbital pericenters, $r_{\rm peri}$, correlate with their dark matter halo structural properties. Specifically, at fixed maximum circular velocity, $V_{\rm max}$, subhalos with smaller $r_{\rm peri}$ are more concentrated (have smaller $r_{\rm max}$ values) and have lost more mass, with larger peak circular velocities, $V_{\rm peak}$, prior to infall. These trends provide information that can tighten constraints on the inferred $V_{\rm max}$ and $V_{\rm peak}$ values for known Milky Way satellites. We illustrate this using published pericenter estimates enabled by Gaia for the nine classical Milky Way dwarf spheroidal satellites. The two densest dSph satellites (Draco and Ursa Minor) have relatively small pericenters, and this pushes their inferred $r_{\rm max}$ and $V_{\rm max}$ values lower than they would have been without pericenter information. For Draco, we infer $V_{\rm max} = 23.5 \, \pm 3.3$ km s$^{-1}$ (compared to $27.3 \, \pm 7.1$ km s$^{-1}$ without pericenter information). Such a shift exacerbates the traditional Too Big to Fail problem. Draco's peak circular velocity range prior to infall narrows from $V_{\rm peak} = 21 - 49$ km s$^{-1}$ without pericenter information to $V_{\rm peak} = 25-37$ km s$^{-1}$ with the constraint. Over the full population of classical dwarf spheroidals, we find no correlation between $V_{\rm peak}$ and stellar mass today, indicative of a high level of stochasticity in galaxy formation at stellar masses below $\sim 10^7$ M$_\odot$. As proper motion measurements for dwarf satellites become more precise, they should enable useful priors on the expected structure and evolution of their host dark matter subhalos., 6 Pages, 4 figures, submitted to MNRAS

We present a suite of FIRE-2 cosmological zoom-in simulations of isolated field dwarf galaxies, all with masses of $M_\mathrm{halo} \approx 10^{10}\,$M$_\odot$ at $z=0$, across a range of dark matter models. For the first time, we compare how both self-interacting dark matter (SIDM) and/or warm dark matter (WDM) models affect the assembly histories as well as the central density structure in fully hydrodynamical simulations of dwarfs. Dwarfs with smaller stellar half-mass radii (r$_{1/2}, Comment: 16 pages, 12 figures. V2: matches version accepted by MNRAS

We compare a suite of four simulated dwarf galaxies formed in 10$^{10} M_{\odot}$ haloes of collisionless Cold Dark Matter (CDM) with galaxies simulated in the same haloes with an identical galaxy formation model but a non-zero cross-section for dark matter self-interactions. These cosmological zoom-in simulations are part of the Feedback In Realistic Environments (FIRE) project and utilize the FIRE-2 model for hydrodynamics and galaxy formation physics. We find the stellar masses of the galaxies formed in Self-Interacting Dark Matter (SIDM) with $\sigma/m= 1\, cm^2/g$ are very similar to those in CDM (spanning $M_{\star} \approx 10^{5.7 - 7.0} M_{\odot}$) and all runs lie on a similar stellar mass -- size relation. The logarithmic dark matter density slope ($\alpha=d\log \rho / d\log r$) in the central $250-500$ pc remains steeper than $\alpha= -0.8$ for the CDM-Hydro simulations with stellar mass $M_{\star} \sim 10^{6.6} M_{\odot}$ and core-like in the most massive galaxy. In contrast, every SIDM hydrodynamic simulation yields a flatter profile, with $\alpha >-0.4$. Moreover, the central density profiles predicted in SIDM runs without baryons are similar to the SIDM runs that include FIRE-2 baryonic physics. Thus, SIDM appears to be much more robust to the inclusion of (potentially uncertain) baryonic physics than CDM on this mass scale, suggesting SIDM will be easier to falsify than CDM using low-mass galaxies. Our FIRE simulations predict that galaxies less massive than $M_{\star} < 3 \times 10^6 M_{\odot}$ provide potentially ideal targets for discriminating models, with SIDM producing substantial cores in such tiny galaxies and CDM producing cusps., Comment: 10 Pages, 7 figures, submitted to MNRAS

In hierarchical models of structure formation, the first galaxies form in low-mass dark matter potential wells, probing the behavior of dark matter on kiloparsec (kpc) scales. Even though these objects are below the detection threshold of current telescopes, future missions will open an observational window into this emergent world. In this Letter we investigate how the first galaxies are assembled in a `fuzzy' dark matter (FDM) cosmology where dark matter is an ultralight $\sim 10^{-22}$~eV boson and the primordial stars are expected to form along dense dark matter filaments. Using a first-of-its-kind cosmological hydrodynamical simulation, we explore the interplay between baryonic physics and unique wavelike features inherent to FDM. In our simulation, the dark matter filaments show coherent interference patterns on the boson de Broglie scale and develop cylindrical soliton-like cores which are unstable under gravity and collapse into kpc-scale spherical solitons. Features of the dark matter distribution are largely unaffected by the baryonic feedback. On the contrary, the distributions of gas and stars, which do form along the entire filament, exhibit central cores imprinted by dark matter -- a smoking gun signature of FDM., Comment: 6 pages, 2 figures, accepted to Phys. Rev. Lett

Bose-Einstein Condensate Dark Matter (BECDM; also known as Fuzzy Dark Matter) is motivated by fundamental physics and has recently received significant attention as a serious alternative to the established Cold Dark Matter (CDM) model. We perform cosmological simulations of BECDM gravitationally coupled to baryons and investigate structure formation at high redshifts ($z \gtrsim 5$) for a boson mass $m=2.5\cdot 10^{-22}~{\rm eV}$, exploring the dynamical effects of its wavelike nature on the cosmic web and the formation of first galaxies. Our BECDM simulations are directly compared to CDM as well as to simulations where the dynamical quantum potential is ignored and only the initial suppression of the power spectrum is considered -- a Warm Dark Matter-like ("WDM") model often used as a proxy for BECDM. Our simulations confirm that "WDM" is a good approximation to BECDM on large cosmological scales even in the presence of the baryonic feedback. Similarities also exist on small scales, with primordial star formation happening both in isolated haloes and continuously along cosmic filaments; the latter effect is not present in CDM. Global star formation and metal enrichment in these first galaxies are delayed in BECDM/"WDM" compared to the CDM case: in BECDM/"WDM" first stars form at $z\sim 13$/$13.5$ while in CDM star formation starts at $z\sim 35$. The signature of BECDM interference, not present in "WDM", is seen in the evolved dark matter power spectrum: although the small scale structure is initially suppressed, power on kpc scales is added at lower redshifts. Our simulations lay the groundwork for realistic simulations of galaxy formation in BECDM., Comment: 19 pages, 12 figures, submitted to MNRAS

The quenching "maintenance'" and related "cooling flow" problems are important in galaxies from Milky Way mass through clusters. We investigate this in halos with masses $\sim 10^{12}-10^{14}\,{\rm M}_{\odot}$, using non-cosmological high-resolution hydrodynamic simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model. We specifically focus on physics present without AGN, and show that various proposed "non-AGN" solution mechanisms in the literature, including Type Ia supernovae, shocked AGB winds, other forms of stellar feedback (e.g. cosmic rays), magnetic fields, Spitzer-Braginskii conduction, or "morphological quenching" do not halt or substantially reduce cooling flows nor maintain "quenched" galaxies in this mass range. We show that stellar feedback (including cosmic rays from SNe) alters the balance of cold/warm gas and the rate at which the cooled gas within the galaxy turns into stars, but not the net baryonic inflow. If anything, outflowing metals and dense gas promote additional cooling. Conduction is important only in the most massive halos, as expected, but even at $\sim 10^{14}\,{\rm M}_{\odot}$ reduces inflow only by a factor $\sim 2$ (owing to saturation effects and anisotropic suppression). Changing the morphology of the galaxies only slightly alters their Toomre-$Q$ parameter, and has no effect on cooling (as expected), so has essentially no effect on cooling flows or maintaining quenching. This all supports the idea that additional physics, e.g., AGN feedback, must be important in massive galaxies., 16 pages, 12 figures

Building upon results of cosmological simulations of ultra-light scalar field dark matter (SFDM), we present a comprehensive model for the density profiles of SFDM haloes as a function of halo virial mass $M_{\rm h}$ and scalar field mass $m$. The central regions of SFDM haloes are dominated by solitons with characteristic densities that increase with increasing halo mass and asymptote to CDM-like profiles at large radii. For scalar field masses $m \sim 10^{-22}$ eV, consistent with large-scale structure observations, $M_{\rm h} \sim 10^{10} \,M_\odot$ haloes have lower core densities than their Cold Dark Matter (CDM) counterparts and this alleviates the Too Big to Fail problem (TBTF) in a regime where feedback is less effective. However, higher-mass SFDM haloes with $M_{\rm h} \sim 10^{11} \,M_\odot$ are denser than their CDM counterparts at small, observationally relevant radii. We use rotation curves of $V \sim 100$ km s$^{-1}$ galaxies from the SPARC database to show that SFDM exacerbates the cusp/core and central density problems seen in CDM at this scale. We conclude that if the conventional cosmological SFDM scaling relations are correct, then baryonic feedback is required to lower densities in SFDM haloes even more so than in CDM. This motivates cosmological and self-consistent hydrodynamic simulations of SFDM to determine whether central soliton structure can be altered by realistic feedback implementations., Comment: 13 pages, 8 figures, submitted to MNRAS

The quenching `maintenance' and `cooling flow' problems are important from the Milky Way through massive cluster elliptical galaxies. Previous work has shown that some source of energy beyond that from stars and pure magnetohydrodynamic processes is required, perhaps from AGN, but even the qualitative form of this energetic input remains uncertain. Different scenarios include thermal `heating,' direct wind or momentum injection, cosmic ray heating or pressure support, or turbulent `stirring' of the intra-cluster medium (ICM). We investigate these in $10^{12}-10^{14}\,{\rm M}_{\odot}$ halos using high-resolution non-cosmological simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, including simplified toy energy-injection models, where we arbitrarily vary the strength, injection scale, and physical form of the energy. We explore which scenarios can quench without violating observational constraints on energetics or ICM gas. We show that turbulent stirring in the central $\sim100\,$kpc, or cosmic-ray injection, can both maintain a stable low-SFR halo for $>$Gyr timescales with modest energy input, by providing a non-thermal pressure which stably lowers the core density and cooling rates. In both cases, associated thermal-heating processes are negligible. Turbulent stirring preserves cool-core features while mixing condensed core gas into the hotter halo and is by far the most energy efficient model. Pure thermal heating or nuclear isotropic momentum injection require vastly larger energy, are less efficient in lower-mass halos, easily over-heat cores, and require fine-tuning to avoid driving unphysical temperature gradients or gas expulsion from the halo center., 22 pages, 16 figures

Recent analysis of the rotation curves of a large sample of galaxies with very diverse stellar properties reveal a relation between the radial acceleration purely due to the baryonic matter and the one inferred directly from the observed rotation curves. Assuming the dark matter (DM) exists, this acceleration relation is tantamount to an acceleration relation between DM and baryons. This leads us to a universal maximum acceleration for all halos. Using the latter in DM profiles that predict inner cores implies that the central surface density $\mu_{DM} = \rho_s r_s$ must be a universal constant, as suggested by previous studies in selected galaxies, revealing a strong correlation between the density $\rho_s$ and scale $r_s$ parameters in each profile. We then explore the consequences of the constancy of $\mu_{DM}$ in the context of the ultra-light scalar field dark matter model (SFDM). We find that for this model $\mu_{DM} = 648 \, M_\odot {\rm pc}^{-2}$, and that the so-called WaveDM soliton profile should be an universal feature of the DM halos. Comparing with data from the Milky Way and Andromeda satellites, we find that they are consistent with a boson mass of the scalar field particle of the order of $10^{-21} \, {\rm eV}/c^2$, which puts the SFDM model in agreement with recent cosmological constraints., Comment: 6 pages, 1 figure

Recent high-quality observations of dwarf and low surface brightness (LSB) galaxies have shown that their dark matter (DM) halos prefer flat central density profiles. On the other hand the standard cold dark matter model simulations predict a more cuspy behavior. Feedback from star formation has been widely used to reconcile simulations with observations, this might be successful in field dwarf galaxies but its success in low mass galaxies remains uncertain. One model that have received much attention is the scalar field dark matter model. Here the dark matter is a self-interacting ultra light scalar field that forms a cosmological Bose-Einstein condensate, a mass of $10^{-22}$eV/c$^2$ is consistent with flat density profiles in the centers of dwarf spheroidal galaxies, reduces the abundance of small halos, might account for the rotation curves even to large radii in spiral galaxies and has an early galaxy formation. The next generation of telescopes will provide better constraints to the model that will help to distinguish this particular alternative to the standard model of cosmology shedding light into the nature of the mysterious dark matter., Comment: 6 pages, to appear in: Proceedings of the Fourteenth Marcel Grossman Meeting on General Relativity

We analyse the velocity dispersion for eight of the Milky Way dwarf spheroidal satellites in the context of finite temperature scalar field dark mater. In this model the finite temperature allows the scalar field to be in configurations that possess excited states, a feature that has proved to be necessary in order to explain the asymptotic rotational velocities found in low surface brightness (LSB) galaxies. In this work we show that excited states are not only important in large galaxies but also have visible effects in dwarf spheroidals. Additionally, we stress that contrary to previous works where the scalar field dark matter haloes are consider to be purely Bose-Einstein condensates, the inclusion of excited states in these halo configurations provides a consistent framework capable of describing LSBs and dwarf galaxies of different sizes without arriving to contradictions within the scalar field dark matter model. Using this new framework we find that the addition of excited states accounts very well for the raise in the velocity dispersion in Milky Way dwarf spheroidal galaxies improving the fit compared to the one obtained assuming all the DM to be in the form of a Bose Einstein Condensate., 7 pages, 2 figures. Accepted for publication in Phys. Rev. D

We investigate the gravitational constraints imposed to dark matter halos in the context of finite temperature scalar field dark matter. We find constraints to produce multiple images by dark matter only, we show that there are differences with respect to the full Bose Einstein condensate halo when the temperature of the scalar field in dark matter halos is taken into account. The non zero temperature allows the scalar field to be in excited states and recently, their existence has proved to be necessary to fit rotation curves of dark matter dominated galaxies of all sizes, it also explained the non universality of the halo density profiles. Therefore, we expect that combining observations of rotation curves and strong lensing systems can give us a clue to the nature of dark matter. Finally, we propose a method to identify the excited state of a strong lens halo, knowing various halo excited states can provide information of the scalar field dark matter halo evolution which can be tested using numerical simulations., Comment: 4 pages, 2 figures

Recent high-quality observations of low surface brightness (LSB) galaxies have shown that their dark matter (DM) halos prefer flat central density profiles. On the other hand, the standard cold dark matter model simulations predict a more cuspy behavior. One mechanism to reconcile the simulations with the observed data is the feedback from star formation, this might be successful in isolated dwarf galaxies but its success in LSB galaxies remains unclear. Additionally, including too much feedback in the simulations is a double-edged sword, in order to obtain a cored DM distribution from an initially cuspy one, the feedback recipes usually require to remove a large quantity of baryons from the center of galaxies, however, some feedback recipes produce twice more satellite galaxies of a given luminosity and with much smaller mass to light ratios from those that are observed. Therefore, one DM profile that produces cores naturally and that does not require large amounts of feedback would be preferable. We find both requirements to be satisfied in the scalar field dark matter model. Here, we consider that the dark matter is an auto-interacting real scalar field in a thermal bath at temperature T with an initial $Z_2$ symmetric potential, as the universe expands, the temperature drops so that the $Z_2$ symmetry is spontaneously broken and the field rolls down to a new minimum. We give an exact analytic solution to the Newtonian limit of this system and show that it can satisfy the two desired requirements and that the rotation curve profile is not longer universal., 11 pages, 3 figures, this version matches the one accepted for publication in the Astrophysical Journal

The cold dark matter (CDM) model has two unsolved issues: simulations overpredict the satellite abundance around the Milky Way (MW) and it disagrees with observations of the central densities of dwarf galaxies which prefer constant density (core) profiles.One alternative explanation known as the scalar field dark matter (SFDM) model, assumes that the dark matter is a scalar field of mass($\sim 10^{-22}$ eV/$c^2$); this model can reduce the overabundance issue due to the lack of halo formation below a mass scale of $\sim 10^8$M$_{\odot}$ and successfully fits the density distribution in dwarfs. One of the attractive features of the model is predicting core profiles in halos, although the determination of the core sizes is set by fitting the observational data. We perform \textit{N}-body simulations to explore the influence of tidal forces over a stellar distribution embedded in a SFDM halo orbiting a MW-like SFDM host halo with a disk. Our simulations intend to test the viability of SFDM as an alternative model by comparing the tidal effects that result in this paradigm with those obtained in CDM for similar mass halos. We found that galaxies in subhalos with core profiles and high central densities survive for 10 Gyr. The same occurs for galaxies in low density subhalos located far from the host disk influence, whereas satellites in low density DM halos and in tight orbits can eventually be stripped of stars. We conclude that SFDM shows consistency with results from CDM for dwarf galaxies, but naturally offer a possibility to solve the missing satellite problem., Comment: 17 pages, 6 figures, matches the accepted version in ApJ

Using high-resolution simulations with explicit treatment of stellar feedback physics based on the FIRE (Feedback In Realistic Environments) project, we study how galaxy formation and the interstellar medium (ISM) are affected by magnetic fields, anisotropic Spitzer–Braginskii conduction and viscosity, and sub-grid metal diffusion from unresolved turbulence. We consider controlled simulations of isolated (non-cosmological) galaxies but also a limited set of cosmological ‘zoom-in’ simulations. Although simulations have shown significant effects from these physics with weak or absent stellar feedback, the effects are much weaker than those of stellar feedback when the latter is modelled explicitly. The additional physics have no systematic effect on galactic star formation rates (SFRs). In contrast, removing stellar feedback leads to SFRs being overpredicted by factors of ∼10–100. Without feedback, neither galactic winds nor volume-filling hot-phase gas exist, and discs tend to runaway collapse to ultra-thin scaleheights with unphysically dense clumps congregating at the galactic centre. With stellar feedback, a multi-phase, turbulent medium with galactic fountains and winds is established. At currently achievable resolutions and for the investigated halo mass range 10^(10)–10^(13) M⊙, the additional physics investigated here (magnetohydrodynamic, conduction, viscosity, metal diffusion) have only weak (∼10 per cent-level) effects on regulating SFR and altering the balance of phases, outflows or the energy in ISM turbulence, consistent with simple equipartition arguments. We conclude that galactic star formation and the ISM are primarily governed by a combination of turbulence, gravitational instabilities and feedback. We add the caveat that active galactic nucleus feedback is not included in the present work.