Dwarf Galaxies Research Articles

Forming off-center massive black hole binaries in dwarf galaxies through Jacobi capture

It is well established that black holes reside in the central regions of virtually all types of known galaxies. Recent observational and numerical studies however challenge this picture, suggesting that intermediate-mass black holes in dwarf galaxies may be found on orbits far from the center. In particular, constant-density cores minimize orbital energy losses due to dynamical friction, and allow black holes to settle on stable off-center orbits. Using controlled simulations, we study the dynamics of off-center black holes in dwarf galaxies with such cores. We propose a new scenario to describe off-center mergers of massive black holes, starting with a Jacobi capture. We focus on initially circular and co-planar black hole orbits and explore a large parameter space of black hole masses and orbital parameters. We find that Jacobi captures are a complex and chaotic phenomenon that occurs in about 13% of cases in this simplified setup, and we quantify how the likelihood of capture depends on the simulation parameters. We note that this percentage is likely an upper limit of the general case. Nevertheless, we show that Jacobi captures in cored dwarf galaxies can facilitate the formation of off-center black hole binaries, and that this process is sufficiently common to have a substantial effect on black hole growth. While our setup only allows for temporary captures, we expect dissipative forces from baryons and post-Newtonian corrections to maintain the captures over time and to lead to the formation of stable binary systems. This motivates future studies of the effectiveness of such dissipative forces, within stripped nuclei or globular clusters, in forming stable bound systems.

Bursting with Feedback: The Relationship between Feedback Model and Bursty Star Formation Histories in Dwarf Galaxies

Due to their inability to self-regulate, ultrafaint dwarfs are sensitive to prescriptions in subgrid physics models that converge and regulate at higher masses. We use high-resolution cosmological simulations to compare the effect of bursty star formation histories (SFHs) on dwarf galaxy structure for two different subgrid supernova (SN) feedback models, superbubble and blastwave, in dwarf galaxies with stellar masses from 5000 < M */M ⊙ < 109. We find that in the “MARVEL-ous Dwarfs” suite both feedback models produce cored galaxies and reproduce observed scaling relations for luminosity, mass, and size. Our sample accurately predicts the average stellar metallicity at higher masses, however low-mass dwarfs are metal poor relative to observed galaxies in the Local Group. We show that continuous bursty star formation and the resulting stellar feedback are able to create dark matter (DM) cores in the higher dwarf galaxy mass regime, while the majority of ultrafaint and classical dwarfs retain cuspy central DM density profiles. We find that the effective core formation peaks at M */M halo ≃ 5 × 10−3 for both feedback models. Both subgrid SN models yield bursty SFHs at higher masses; however, galaxies simulated with superbubble feedback reach maximum mean burstiness values at lower stellar mass fractions relative to blastwave feedback. As a result, core formation may be better predicted by stellar mass fraction than the burstiness of SFHs.

What drives the corpulence of galaxies?

Nearby dwarf galaxies display a variety of effective radii (sizes) at a given stellar mass, suggesting different evolution scenarios according to their final “stellar” size. The TNG hydrodynamical simulations present a bimodality in the z = 0 size–mass relation (SMRz0) of dwarf galaxies, at r1/2, ⋆ ∼ 450 pc. Using the TNG50 simulation, we explored the evolution of the most massive progenitors of dwarf galaxies (z = 0 log(M⋆/M⊙) between 8.4 and 9.2) that end up as central galaxies of their groups. We split these dwarfs into three classes of the SMRz0: “Normals” from the central spine of the main branch, and “Compacts” from the secondary branch as well as the lower envelope of the main branch. Both classes of Compacts see their stellar sizes decrease from z ∼ 1 onwards in contrast to Normals, while the sizes of the gas and dark matter (DM) components continue to increase (as for Normals). A detailed analysis reveals that Compacts live in poorer environments, and thus suffer fewer major mergers from z = 0.8 onwards, which otherwise would pump angular momentum into the gas, allowing strong gas inflows, producing inner star formation, and thus leading to the buildup of a stellar core. Compacts are predicted to be rounder and to have bluer cores. Compact dwarfs of similar sizes are observed in the GAMA survey, but the bimodality in size is less evident and the most compact dwarfs tend to be passive rather than star forming, as in TNG50. Our conclusions should therefore be confirmed with future cosmological hydrodynamical simulations.

Exploring the Evolution of a Dwarf Spheroidal Galaxy with Smoothed-particle Hydrodynamics Simulations. I. Stellar Feedback

A fundamental question regarding the evolution of dwarf spheroidal galaxies is the identification of the key physical mechanisms responsible for gas depletion. Here, we focus on the study of stellar feedback in isolated dwarf spheroidal galaxies by performing numerical simulations using a modified version of the smoothed-particle hydrodynamics code GADGET-3. The Milky Way satellite Leo II (PGC 34176) in the Local Group was considered as our default model dwarf galaxy. The parameter space for the stellar feedback models was explored to match observational constraints of Leo II, such as residual gas mass, total mass within the tidal radius, star formation history, final stellar mass, stellar ages, and metallicity. Additionally, we examined the impact of the binary fraction of stars, initial mass function, dark matter halo mass, and initial gas reservoir. Many simulations revealed recent star formation quenching due to stellar feedback. In general, the gas depletion, expected star formation history, total mass of stars, and total mass within the tidal radius were adequately reproduced in the simulations when compared to observational estimates. However, there were discrepancies in the distribution of stellar ages and metallicities, which suggested that the cosmic gas infall would play a more complex role in our dwarf spheroidal galaxy than captured by a monolithic infall scenario. Our results suggest that currently quenched dwarf galaxies may not necessarily need to evolve within clusters or groups and that stellar feedback alone could be a sufficient factor in shaping at least some of these galaxies as we observe them today.

A theoretical perspective on the almost dark galaxy Nube: exploring the fuzzy dark matter model

In recent astronomical observations, an almost dark galaxy, designated as Nube, has unveiled an intriguing anomaly in its stellar distribution. Specifically, Nube exhibits an exceptionally low central brightness, with the 2D half-light radius of its stars far exceeding the typical values found in dwarf galaxies, and even surpassing those observed in ultra-diffuse galaxies (UDGs). This phenomenon is difficult to explain within the framework of cold dark matter (CDM). Meanwhile, due to its ultralight particle mass, fuzzy dark matter (FDM) exhibits a de Broglie wavelength on the order of kiloparsecs under the typical velocities of galaxies. The interference between different modes of the FDM wave gives rise to fluctuations in the gravitational field, which can lead to the dynamical heating of stars within galaxies, resulting in an expansion of their spatial distribution. In this paper, we aim to interpret the anomalous stellar distribution observed in Nube as a consequence of the dynamical heating effect induced by FDM. Our findings suggest that a FDM particle mass around 1-2 × 10-23 eV can effectively account for this anomaly. And we propose that the FDM dynamical heating effect provides a new insight into understanding the formation of field UDGs.

The Variation in the Galaxy-wide Initial Mass Function for Low-mass Stars: Modeling and Observational Insights

The stellar initial mass function (IMF) characterizes the mass distribution of newly formed stars in various cosmic environments, serving as a fundamental assumption in astrophysical research. Recent findings challenge the prevalent notion of a universal and static IMF, proposing instead that the IMF’s shape is contingent upon the star formation environment. In this study, we analyze the galaxy-wide variation in the IMF for low-mass stars in both dwarf and massive galaxies with diverse observational methods. Despite systematic discrepancies between different approaches, an IMF model with a metallicity-dependent slope for the low-mass stars aligns with the majority of observations, indicating a high degree of uniformity in the star formation processes across the Universe. We also emphasize the need for a more comprehensive understanding of the variation in the low-mass IMF, considering measurement biases and factors beyond metallicity.

Anisotropies in the spatial distribution and kinematics of dwarf galaxies in the Local Group and beyond

ABSTRACT The mass distribution in the Local Group (LG), dominated by the Andromeda (M31) and Milky Way (MW) pair, is highly anisotropic. We use the APOSTLE simulations to examine how this anisotropy manifests on the spatial distribution and kinematics of dwarf galaxies out to a distance of $d_{\rm MW}\sim 3$ Mpc from the MW. The simulations indicate a preference for dwarfs to be located near the axis defined by the MW-M31 direction, even for dwarfs in the LG periphery (LGP; i.e. at distances $1.25\ \lt\ d_{\rm MW}/$Mpc $\lt\ 3$). The ‘Hubble flow’ in the periphery is also affected; at fixed $d_{\rm MW}$ the mean recession speed, $\langle V_{\rm rad} \rangle$, varies with angular distance to M31, peaking in the anti-M31 direction and reaching a minimum behind M31. The M31-MW mass decelerates the local expansion; the LG ‘turnaround radius’ (i.e. where $\langle V_{\rm rad} \rangle =0$) in APOSTLE is at $r \sim 1.25$ Mpc from the LG barycentre and the pure Hubble flow (where $\langle V_{\rm rad} \rangle \sim H_0*d$) is reached beyond $r\sim 3$ Mpc. The predicted flow is very cold, with a barycentric dispersion of $\lt 40$ km s−1. Comparing these predictions with observations yields mixed results. There is little evidence for a preferred alignment of dwarfs along the MW-M31 direction, but some evidence for an angular anisotropy in $\langle V_{\rm rad} \rangle$. Although the ‘coldness’ of the Hubble flow is consistent with the simulations, it is less decelerated: relative to the MW all galaxies beyond $d_{\rm MW} \sim 1.25$ Mpc seem to be already on a pure Hubble flow. We argue that these oddities may result at least partly from incompleteness and inhomogeneous sky coverage in our current inventory of nearby dwarfs.

Constraining the Low-mass End of the Black Hole Mass Function and the Active Fraction of the Intermediate-mass Black Holes

We investigate the black hole mass function (BHMF) and the Eddington ratio distribution function (ERDF), focusing on the intermediate-mass black holes (IMBHs) with masses down to M • ∼ 104 M ⊙. Based on the active galactic nuclei (AGNs) with a detected broad Hα emission line, we construct a sample of 14,242 AGNs at redshift z < 0.35, including 243 IMBHs with M • < 106 M ⊙. By jointly modeling the BHMF and ERDF via the maximum posterior estimation, we find that the BHMF peaks at ∼106 M ⊙ and exhibits a relatively constant value of 10−4 Mpc−3 dex−1 at the low-mass end. By comparing the derived BHMF of type 1 AGNs with the galaxy mass function based on the updated black hole mass–host galaxy stellar mass relation, we derive the active fraction. We also determine the active fraction for all AGNs using the upper and lower limit of the type 1 fraction. The active fraction decreases from 15%–40% for massive galaxies (M ⋆ > 1010 M ⊙) to lower than ∼2% for dwarf galaxies with M ⋆ ∼ 108 M ⊙. These results suggest that the black hole occupation fraction is expected to be ∼50% for low-mass galaxies (M ⋆ ∼ 108.5–109 M ⊙) if the duty cycle is similar between IMBHs and supermassive black holes.

Figuring Out Gas and Galaxies in Enzo (FOGGIE). VII. The (Dis)assembly of Stellar Halos

Over the next decade, the astronomical community will be commissioning multiple wide-field observatories well suited for studying stellar halos in both integrated light and resolved stars. In preparation for this, we use five high-resolution cosmological simulations of Milky Way–like galaxies from the FOGGIE suite to explore the properties and components of stellar halos. These simulations are run with high time (5 Myr) and stellar mass (1000 M ⊙) resolution to better model the properties and origins of low-density regions like stellar halos. We find that the FOGGIE stellar halos have masses, metallicity gradients, and surface brightness profiles that are consistent with observations. In agreement with other simulations, the FOGGIE stellar halos receive 30%–40% of their mass from in situ stars. However, this population is more centrally concentrated in the FOGGIE simulations and therefore does not contribute excess light to the halo outskirts. The remaining stars are accreted from ∼10–50 other galaxies, with the majority of the accreted mass originating in two to four galaxies. While the inner halo (r < 50 kpc) of each FOGGIE galaxy has a large number of contributors, the halo outskirts of three of the five galaxies are primarily made up of stars from only a few contributors. We predict that upcoming wide-field observatories, like the Nancy Grace Roman Space Telescope, will probe stellar halos around Milky Way–like galaxies out to ∼100 kpc in integrated light and will be able to distinguish the debris of dwarf galaxies with extended star formation histories from the underlying halo with resolved color–magnitude diagrams.

A census of new globular clusters in the Galactic bulge

Context. The number of known globular clusters in the Galactic bulge has been increasing steadily thanks to different new surveys. Aims. The aim of this study is to provide a census of the newly revealed globular clusters in the Galactic bulge, and analyze their characteristics. Methods. In recent years, many globular clusters have been discovered or identified. The stellar populations to which they belong are indicated in their original studies: they are mostly bulge clusters, with some identified as disk or halo members. We collected 41 new globular clusters revealed in the last decade and compared them to the known bulge clusters. Results. The new clusters are intrinsically faint with MV of around −6.0 mag. The distance to the Sun of the ensemble of well-known and new bulge clusters is compatible with the Galactocentric distance measurements from the Galactic black hole location. The ensemble sample shows metallicity peaks at [Fe/H] ∼ −1.08 ± 0.35 and −0.51 ± 0.25 dex, confirming previous findings. The age–metallicity relation of the new clusters younger than 10 Gyr is compatible with that of the ex situ samples of the dwarf galaxies Sagittarius, Canis Majoris, and Gaia-Enceladus-Sausage. The clusters with ages between 11.5 and 13.5 Gyr show no age–metallicity relation, because they are all old. This is compatible with their formation in situ in the early Galaxy.

Satellite group infall into the Milky Way: Exploring the Crater-Leo case with new HST proper motions

Context. Within Λ cold dark matter (ΛCDM) simulations, Milky Way-like galaxies accrete some of their satellite galaxies in groups of 3–5 members rather than individually, and this has been suggested as a possible mechanism driving the formation of satellite planes. Objects accreted in groups are expected to share similar specific total energy and angular momentum, and to also have identical orbital planes and directions. Aims. Looking at observatio ns of Milky Way satellites, the dwarf galaxies Leo II, IV, V, and Crater II, and the star cluster Crater 1 were proposed to be a vestige of group infall. The suggested ‘Crater-Leo group’ shows a coherent distance gradient and all these objects align along a great circle on the sky. We used proper motion data to investigate whether the phase-space distribution of the members of the proposed group are indeed consistent with group infall. Methods. To further investigate this possibility, we used Gaia Data Release 3 (DR3) and new Hubble Space Telescope (HST) proper motions – namely, (μα*, μδ) = (−0.1921 ± 0.0514, −0.0686 ± 0.0523) mas yr−1 for Leo IV and (μα*, μδ) = (0.1186 ± 0.1943, −0.1183 ± 0.1704) mas yr−1 for Leo V – to derive accurate orbital properties for the proposed group objects. In addition, we explored other possible members of this putative association. Results. Leo II, Leo IV, and Crater 1 show orbital properties consistent with those we predict from assuming group infall. However, our results suggest that Crater II was not accreted with the rest of the objects. If confirmed with increasingly accurate proper motions in the future, the Crater-Leo objects would appear to constitute the first identified case of a cosmologically expected, typical group infall event, as opposed to the highly hierarchical Magellanic Cloud system.

Modeling Ionized Gas in the Small Magellanic Cloud: The Wolf–Rayet Nebula N76

We present Cloudy modeling of infrared emission lines in the Wolf–Rayet (WR) nebula N76 caused by one of the most luminous and hottest WR stars in the low metallicity Small Magellanic Cloud. We use spatially resolved mid-infrared Spitzer/InfRared Spectrograph and far-infrared Herschel/PACS spectroscopy to establish the physical conditions of the ionized gas. The spatially resolved distribution of the emission allows us to constrain properties much more accurately than using spatially integrated quantities. We construct models with a range of constant hydrogen densities between nH = 4–10 cm−3 and a stellar wind-blown cavity of 10 pc, which reproduces the intensity and shape of most ionized gas emission lines, including the high ionization lines [O iv] and [Ne v], as well as [S iii], [S iv], [O iii], and [Ne iii]. Our models suggest that the majority of [Si ii] emission (91%) is produced at the edge of the H ii region around the transition between ionized and atomic gas while very little of the [C ii] (<5%) is associated with the ionized gas. The physical conditions of N76 are characterized by a hot HII region with a maximum electron temperature of T e ∼ 24,000 K, electron densities that range from n e ∼ 4 to 12 cm−3, and high ionization parameters of log(U)∼−1.15to−1.77. By analyzing a low-metallicity WR nebula with a single ionization source, this work gives valuable insights into the impact WR stars have on the galaxy-integrated ionized gas properties in nearby dwarf galaxies.

New H i Observations Toward the NGC 5055 Galaxy Group with FAST

We report a new high-sensitivity H i mapping observation of the NGC 5055 galaxy group over an area of 1.°5 × 0.°75 with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Our observation reveals that the warped H i disk of NGC 5055 is more extended than what was previously observed by WSRT, out to 23.′9 (61.7 kpc). The total H i mass of NGC 5055 is determined to be ∼1.1 × 1010 M ⊙. We identified three H i clouds with H i masses of the order of ∼107 M ⊙ at the southeastern edge of the H i disk, as well as a candidate high-velocity cloud with an H i mass of (1.2 ± 0.5) × 106 M ⊙ to the north of NGC 5055. The H i content of UGCA 337 is robustly detected for the first time by the FAST observations. It has a narrow H i linewidth of W 50 = 17.4 ± 3.8 km s−1 with a total H i mass of (3.5 ± 0.3) × 106 M ⊙. Comparing the gas content and g−r color of UGCA 337 with typical low-mass dwarf galaxies, UGCA 337 appears relatively gas-poor despite its blue color. This suggests that UGCA 337 may have undergone gas stripping in the past. We also analyzed the possible origin of the diffuse H i clouds located at the outskirts of NGC 5055, and speculate that they might be the remnant features of a merger event in the past.

VLBA Astrometry of the Galactic Double Neutron Stars PSR J0509+3801 and PSR J1930–1852: A Preliminary Transverse Velocity Distribution of Double Neutron Stars and its Implications

The mergers of double neutron star (DNS) systems are believed to drive the majority of short γ-ray bursts (SGRBs), while also serving as production sites of heavy r-process elements. Despite being key to (i) confirming the nature of the extragalactic SGRBs, (ii) addressing the poorly understood r-process enrichment in the ultrafaint dwarf galaxies (UFDGs), and (iii) probing the formation process of DNS systems, the space velocity distribution of DNSs is still poorly constrained, due to the small number of DNSs with well-determined astrometry. In this work, we determine new proper motions and parallaxes of two Galactic DNSs, PSR J0509+3801 and PSR J1930−1852, using the Very Long Baseline Array, and we estimate the transverse velocities v ⊥ of all 11 isolated Galactic DNSs having proper-motion measurements in a consistent manner. Our correlation analysis reveals that the DNS v ⊥ is tentatively correlated with three parameters: spin period, orbital eccentricity, and companion mass. With the preliminary v ⊥ distribution, we obtain the following findings. First, the refined v ⊥ distribution is confirmed to agree with the observed displacements of the localized SGRBs from their host galaxy birth sites. Second, we estimate that around 11% and 25% of DNSs remain gravitationally bound to UFDGs with escape velocities of 15 and 25 km s−1, respectively. Hence, the retained DNSs might indeed be responsible for the r-process enrichment confirmed so far in a few UFDGs. Finally, we discuss how a future ensemble of astrometrically determined DNSs may probe the multimodality of the v ⊥ distribution.

HSTPROMO Internal Proper-motion Kinematics of Dwarf Spheroidal Galaxies. I. Velocity Anisotropy and Dark Matter Cusp Slope of Draco

We analyze four epochs of Hubble Space Telescope imaging over 18 yr for the Draco dwarf spheroidal galaxy. We measure precise proper motions for hundreds of stars and combine these with existing line-of-sight (LOS) velocities. This provides the first radially resolved 3D velocity dispersion profiles for any dwarf galaxy. These constrain the intrinsic velocity anisotropy and resolve the mass–anisotropy degeneracy. We solve the Jeans equations in oblate axisymmetric geometry to infer the mass profile. We find the velocity dispersion to be radially anisotropic along the symmetry axis and tangentially anisotropic in the equatorial plane, with a globally averaged value βB¯=−0.20−0.53+0.28 , (where 1 – βB≡〈vtan2〉/〈vrad2〉 in 3D). The logarithmic dark matter (DM) density slope over the observed radial range, Γdark, is −0.83−0.37+0.32 , consistent with the inner cusp predicted in ΛCDM cosmology. As expected given Draco’s low mass and ancient star formation history, it does not appear to have been dissolved by baryonic processes. We rule out cores larger than 487, 717, and 942 pc at 1σ, 2σ, and 3σ confidence, respectively, thus imposing important constraints on the self-interacting DM cross section. Spherical models yield biased estimates for both the velocity anisotropy and the inferred slope. The circular velocity at our outermost data point (900 pc) is 24.19−2.97+6.31kms−1 . We infer a dynamical distance of 75.37−4.00+4.73 kpc and show that Draco has a modest LOS rotation, with v/σ=0.22±0.09 . Our results provide a new stringent test of the so-called “cusp–core” problem that can be readily extended to other dwarfs.

Lyα emission as a sensitive probe of feedback-regulated LyC escape from dwarf galaxies

ABSTRACT Lyα emission is an exceptionally informative tracer of the life cycle of evolving galaxies and the escape of ionizing photons. However, theoretical studies of Lyα emission are often limited by insufficient numerical resolution, incomplete sets of physical models, and poor line-of-sight (LOS) statistics. To overcome such limitations, we utilize here the novel pandora suite of high-resolution dwarf galaxy simulations that include a comprehensive set of state-of-the-art physical models for ionizing radiation, magnetic fields, supernova feedback, and cosmic rays. We post-process the simulations with the radiative transfer code rascas to generate synthetic observations and compare to the observed properties of Lyα emitters. Our simulated Lyα haloes are more extended than the spatial region from which the intrinsic emission emanates, and our spatially resolved maps of spectral parameters of the Lyα emission are very sensitive to the underlying spatial distribution and kinematics of neutral hydrogen. Lyα and LyC emissions display strongly varying signatures along different LOS depending on how each LOS intersects low-density channels generated by stellar feedback. Comparing galaxies simulated with different physics, we find the Lyα signatures to exhibit systematic offsets determined by the different levels of feedback strength and the clumpiness of the neutral gas. Despite this variance, and regardless of the different physics included in each model, we find universal correlations between Lyα observables and LyC escape fraction, demonstrating a robust connection between Lyα and LyC emission. Lyα observations from a large sample of dwarf galaxies should thus give strong constraints on their stellar feedback-regulated LyC escape and confirm their important role for the re-ionization of the Universe.

Overview and public data release of the augmented Auriga Project: cosmological simulations of dwarf and Milky Way-mass galaxies

ABSTRACT We present an extended suite of the Auriga cosmological gravo-magnetohydrodynamical ‘zoom-in’ simulations of 40 Milky Way-mass haloes and 26 dwarf galaxy–mass haloes run with the moving-mesh code arepo. Auriga adopts the Lambda cold dark matter cosmogony and includes a comprehensive galaxy formation physics model following the coupled cosmic evolution of dark matter, gas, stars, and supermassive black holes which has been shown to produce numerically well-converged galaxy properties for Milky Way-mass systems. We describe the first public data release of this augmented suite of Auriga simulations, which includes raw snapshots, group catalogues, merger trees, initial conditions, and supplementary data, as well as public analysis tools with worked examples of how to use the data. To demonstrate the value and robustness of the simulation predictions, we analyse a series of low-redshift global properties that compare well with many observed scaling relations, such as the Tully–Fisher relation, the star-forming (SF) main sequence, and H i gas fraction/disc thickness. Finally, we show that SF gas discs appear to build rotation and velocity dispersion rapidly for $z\gtrsim 3$ before they ‘settle’ into ever-increasing rotation-dispersion ratios ($V/\sigma$). This evolution appears to be in rough agreement with some kinematic measurements from H$\alpha$ observations, and demonstrates an application of how to utilize the released data.

Scalar field dark matter: impact of supernova-driven blowouts on the soliton structure of low-mass dark matter haloes

ABSTRACT We present the first study on the gravitational impact of supernova feedback in an isolated soliton and a spherically symmetric dwarf scalar field dark matter (SFDM) halo of virial mass $1\times 10^{10}\,\mathrm{M_\odot }$. We use a boson mass $m=10^{-22}\,\mathrm{eV\,c^{-2}}$ and a soliton core $r_\mathrm{ c} \approx 0.7$ kpc, comparable to typical half-light radii of Local Group dwarf galaxies. We simulate the rapid gas removal from the centre of the soliton by a concentric external time-dependent Hernquist potential. We explore two scenarios of feedback blowouts: (i) a massive single burst and (ii) multiple consecutive blowouts injecting the same total energy to the system, including various magnitudes for the blowouts in both scenarios. In all cases, we find one single blowout has a stronger effect on reducing the soliton central density. Feedback leads to central soliton densities that oscillate quasi-periodically for an isolated soliton and stochastically for an SFDM halo. The range in the density amplitude depends on the strength of the blowout; however, we observe typical variations of a factor of $\geqslant$2. One important consequence of the stochastic fluctuating densities is that, if we had no prior knowledge of the system evolution, we can only know the configuration profile at a specific time within some accuracy. By fitting soliton profiles at different times to our simulated structures, we found the (1$\sigma$) scatter of their time-dependent density profiles. For configurations within the 1$\sigma$ range, we find the inferred boson mass is typically less than 20 per cent different from the real value used in our simulations. Finally, we compare the observed dynamical masses of field dwarf galaxies in our Local Group with the implied range of viable solitons from our simulations and find good agreement.

The torsion of stellar streams and the overall shape of galactic gravity’s source

Context. Flat rotation curves, v(r), are naturally explained by elongated (prolate) dark matter (DM) distributions, and we have provided competitive fits to the SPARC database. To further probe the geometry of the halo, or the equivalent source of gravity in other formulations, one needs observables outside the galactic plane. Stellar streams, poetically analogous to airplane contrails, but caused by tidal dispersion of massive substructures such as satellite dwarf galaxies, would lie on their own plane (consistently with angular momentum conservation) should the DM-halo gravitational field be spherically symmetric. Tracks resembling entire orbits are seldom available because their periods are commensurable with Hubble time, with streams often presenting themselves as short segments. Aims. Therefore, we aim to establish stellar stream torsion, a local observable that measures the deviation from planarity in differential curve geometry, as a diagnostic providing sensitivity to aspherical DM distributions and ensuring the use of even relatively short streams. Methods. We performed small-scale simulations of tidally distorted star clusters to check that indeed a central force center produces negligible torsion, while distorted halos can generate it. Turning to observational data, we identified among the known streams those that are at the largest distance from the Galactic center, and that are likely not affected by the Magellanic clouds, as the most promising for the study, and by means of polynomial fits we extracted their differential torsion. Results. We find that the torsion of the few known streams that should be sensitive to most of the Milky Way’s DM halo is much larger than expected for a central spherical bulb alone. This is consistent with the nonsphericity of the halo. Conclusions. Future studies of stellar stream torsion with larger samples and further out of the galactic plane should be able to extract the ellipticity of the halo to see whether it is just a slight distortion of a spherical shape or whether it rather resembles a more elongated cigar.

Dark-matter-free Dwarf Galaxy Formation at the Tips of the Tentacles of Jellyfish Galaxies

When falling into a galaxy cluster, galaxies experience a loss of gas due to ram pressure stripping. In particular, disk galaxies lose gas from their disks, and very large tentacles of gas can be formed. Because of the morphology of these stripped galaxies, they have been referred to as jellyfish galaxies. It has been found that star formation is triggered not only in the disk, but also in the tentacles of such jellyfish galaxies. The observed star-forming regions located in the tentacles of those galaxies have been found to be as massive as 3 × 107 M ⊙ and with sizes >100 pc. Interestingly, these parameters in mass and size agree with those of dwarf galaxies. In this work, we make use of the state-of-the-art magnetohydrodynamic (MHD) cosmological simulation IllustrisTNG-50 to study massive jellyfish galaxies with long tentacles. We find that, in the tentacles of TNG-50 jellyfish galaxies, the star formation regions (gas+stars) formed could be as massive as ∼2 × 108 M ⊙. A particular star-forming region was analyzed. This region has a star formation rate of 0.04 M ⊙ yr−1, it is metal-rich, has an average age of 0.46 Gyr, and has a half-mass radius of ∼1 kpc, typical of standard dwarf galaxies. Most importantly, this region is gravitationally self-bound. Overall, we identify a new type of dwarf galaxy being born from the gas tentacles of jellyfish galaxies that, by construction, lacks a dark matter halo.