We use a compilation of disc galaxy rotation curves to assess the role of the luminous component ("baryons") in the rotation curve diversity problem. As in earlier work, we find that rotation curve shape correlates with baryonic surface density: high surface density galaxies have rapidly-rising rotation curves consistent with cuspy cold dark matter halos; slowly-rising rotation curves (characteristic of galaxies with inner mass deficits or "cores") occur only in low surface density galaxies. The correlation, however, seems too weak to be the main driver of the diversity. In addition, dwarf galaxies exhibit a clear trend, from "cuspy" systems where baryons are unimportant in the inner mass budget to "cored" galaxies where baryons actually dominate. This trend constrains the various scenarios proposed to explain the diversity, such as (i) baryonic inflows and outflows during galaxy formation; (ii) dark matter self-interactions; (iii) variations in the baryonic mass structure coupled to rotation velocities through the "mass discrepancy-acceleration relation" (MDAR); or (iv) non-circular motions in gaseous discs. Together with analytical modeling and cosmological hydrodynamical simulations, our analysis shows that each of these scenarios has promising features, but none seems to fully account for the observed diversity. The MDAR, in particular, is inconsistent with the observed trend between rotation curve shape and baryonic importance; either the trend is caused by systematic errors in the data or the MDAR does not apply. The origin of the dwarf galaxy rotation curve diversity and its relation to the structure of cold dark matter halos remains an open issue.

We address the 'plane of satellites problem' by studying planar configurations around two disc galaxies with no late major mergers, formed in zoom-in hydro-simulations. Due to the current lack of good quality kinematic data for M31 satellites, we use only positional information. So far, positional analyses of simulations are unable to find planes as thin and populated as the observed ones. Moreover, they miss systematicity and detail in the plane-searching techniques, as well as in the study of the properties and quality of planes, both in simulations or real data. To fill this gap, i) we extend the 4-galaxy-normal density plot method (Pawlowski et al. 2013) in a way designed to efficiently identify the best quality planes (i.e., thin and populated) without imposing extra constraints on their properties, and ii), we apply it for the first time to simulations. Using zoom-in simulations allows us to mimic MW/M31-like systems regarding the number of satellites involved as well as the galactic disc mass and morphology, in view of possible disc effects. At all timesteps analyzed in both simulations we find satellite planar configurations that are compatible, along given time intervals, with all the spatial characteristics of observed planes identified using the same methodology. However, the fraction of co-orbiting satellites within them is in general low, suggesting time-varying satellite membership. We conclude that high-quality positional planes of satellites are not infrequent in LCDM-formed disc galaxies with a quiet assembly history. Detecting kinematically-coherent, time-persistent planes demands considering the full six-dimensional phase-space information of satellites.

We use N-body simulations to revisit the globular cluster (GC) ``timing problem'' in the Fornax dwarf spheroidal (dSph). In agreement with earlier work, we find that, due to dynamical friction, GCs sink to the center of dark matter halos with a cuspy inner density profile but ``stall'' at roughly 1/3 of the core radius ($r_{\rm core}$) in halos with constant-density cores. The timescales to sink or stall depend strongly on the mass of the GC and on the initial orbital radius, but are essentially the same for either cuspy (NFW) or cored halos normalized to have the same total mass within $r_{\rm core}$. Arguing against a cusp on the basis that GCs have not sunk to the center is thus no different from arguing against a core, unless all clusters are today at $\sim (1/3)\, r_{\rm core}$. This would imply a core radius exceeding $\sim 3$ kpc, much larger than seems plausible in any core-formation scenario. (The average projected distance of Fornax GCs is $\langle R_{\rm GC,Fnx}\rangle\sim 1$ kpc and its effective radius is $\sim 700$ pc.) A simpler explanation is that Fornax GCs have only been modestly affected by dynamical friction, as expected if clusters started orbiting at initial radii of order $\sim 1$-$2$ kpc, just outside Fornax's present-day half-light radius but well within the tidal radius imprinted by Galactic tides. This is not entirely unexpected. Fornax GCs are significantly older and more metal-poor than most Fornax stars, and such populations in dSphs tend to be more spatially extended than their younger and more metal-rich counterparts. Contrary to some earlier claims, our simulations further suggest that GCs do not truly ``stall'' at $\sim 0.3\, r_{\rm core}$, but rather continue decaying toward the center, albeit at reduced rates. We conclude that dismissing the presence of a cusp in Fornax based on the spatial distribution of its GC population is unwarranted.

We present a detailed characterization of planes of satellites in the Milky Way (MW) and M31 systems. To this end we introduce an extension to the '4-galaxy-normal density plot' method \citep{Pawlowski13}, by which plot over-densities signal the normal direction to predominant planes of satellites within a given sample. For a given over-density, the extension provides a \textit{collection} of planes, each including a different number of objects $N_{\rm sat}$. We apply this method to the position data of confirmed MW and M31 satellites and quantify the quality of planes through the outputs of a Tensor of Inertia plane-fitting technique. Plane quality is quantified in terms of population ($N_{\rm sat}$) and flattening (the short-to-long axis ratio $c/a$ or the rms thickness normal to the plane). Therefore, planes with the same population or flattening can be compared with each other allowing us to single-out best-quality planes. For the first time, we study the second-most predominant planar configuration of satellites in M31, singling out a plane with 18 satellite members that shows a quality comparable to the Great Plane of Andromeda (GPoA, with $N_{\rm sat}=19$) despite it being more affected by distance uncertainties. This structure is viewed nearly face-on from the MW and is approximately normal to the GPoA. Overall, we find planes of satellites around the MW and M31 with higher qualities than those previously reported with a given $N_{\rm sat}$. We also show that mass plays no role in determining a satellite's membership or not to the respective best-quality planes.

The significant diversity of rotation curve (RC) shapes in dwarf galaxies has recently emerged as a challenge to LCDM: in dark matter (DM) only simulations, DM halos have a universal cuspy density profile that results in self-similar RC shapes. We compare RC shapes of simulated galaxies from the NIHAO project with observed galaxies from the homogeneous SPARC dataset. The DM halos of the NIHAO galaxies can expand to form cores, with the degree of expansion depending on their stellar-to-halo mass ratio. By means of the $V_{\rm 2kpc}-V_{\rm Rlast}$ relation (where $V_{\rm Rlast}$ is the outermost measured rotation velocity), we show that both the average trend and the scatter in RC shapes of NIHAO galaxies are in reasonable agreement with SPARC: this represents a significant improvement compared to simulations that do not result in DM core formation, suggesting that halo expansion is a key process in matching the diversity of dwarf galaxy RCs. Note that NIHAO galaxies can reproduce even the extremely slowly rising RCs of IC 2574 and UGC 5750. Revealingly, the range where observed galaxies show the highest diversity corresponds to the range where core formation is most efficient in NIHAO simulations, 50$<V_{\rm Rlast}$\km s$^{-1}<$100. A few observed galaxies in this range cannot be matched by any NIHAO RC nor by simulations that predict a universal halo profile. Interestingly, the majority of these are starbursts or emission-line galaxies, with steep RCs and small effective radii. Such galaxies represent an interesting observational target providing new clues to the process/viability of cusp-core transformation, the relationship between starburst and inner potential well, and the nature of DM.

The spectral energy distributions (SEDs) of low-mass low-metallicity (dwarf) galaxies are a challenging piece of the puzzle of galaxy formation in the near Universe. These SEDs show some particular features in the submillimeter to far-infrared wavelength range compared to normal, larger, and metal-richer galaxies that cannot be explained by the current models. These can be summarized as: a broadening of the IR peak, which implies a warmer dust component; an excess of emission in the submm ($\sim$500 $\mu$m), that causes a flattening of the submm/FIR slope; and a very low intensity of PAH emission features. With the aim of explaining these emission patterns, the SEDs of a sample of 27 simulated dwarf galaxies were calculated using the GRASIL-3D radiation transfer code. This code has the particularity that it separately treats the radiative transfer through dust grains within molecular clouds and within the cirrus, the dense and diffuse components of the gas phase, respectively. The simulated galaxies have stellar masses ranging from 10$^6$-10$^9$ M$_\odot$, and were obtained from a single simulation run within a Local Group environment with initial conditions from the CLUES project. We report a study of the IRAS, Spitzer and Herschel bands luminosities, and of the SFRs, dust, and gas (HI and H$_2$) mass contents. We find a satisfactory agreement with observational data, with GRASIL-3D naturally reproducing the specific spectral features mentioned above. We conclude that the GRASIL-3D two-component dust model gives a physical interpretation of the emission of dwarf galaxies: molecular clouds and cirrus represent the warm and cold dust components, respectively, needed to reproduce observational data.

We compare the Baryonic Tully-Fisher relation (BTFR) of simulations and observations of galaxies ranging from dwarfs to spirals, using various measures of rotational velocity Vrot. We explore the BTFR when measuring Vrot at the flat part of the rotation curve, Vflat, at the extent of HI gas, Vlast, and using 20% (W20) and 50% (W50) of the width of HI line profiles. We also compare with the maximum circular velocity of the parent halo, Vmax, within dark matter only simulations. The different BTFRs increasingly diverge as galaxy mass decreases. Using Vlast one obtains a power law over four orders of magnitude in baryonic mass, with slope similar to the observed BTFR. Measuring Vflat gives similar results as Vlast when galaxies with rising rotation curves are excluded. However, higher rotation velocities would be found for low mass galaxies if the cold gas extended far enough for Vrot to reach a maximum. W20 gives a similar slope as Vlast but with slightly lower values of Vrot for low mass galaxies, although this may depend on the extent of the gas in your galaxy sample. W50 bends away from these other relations toward low velocities at low masses. By contrast, Vmax bends toward high velocities for low mass galaxies, as cold gas does not extend out to the radius at which halos reach Vmax. Our study highlights the need for careful comparisons between observations and models: one needs to be consistent about the particular method of measuring Vrot, and precise about the radius at which velocities are measured.

Using 22 hydrodynamical simulated galaxies in a LCDM cosmological context we recover not only the observed baryonic Tully-Fisher relation, but also the observed "mass discrepancy--acceleration" relation, which reflects the distribution of the main components of the galaxies throughout their disks. This implies that the simulations, which span the range 52 < V$_{\rm flat}$ < 222 km/s where V$_{\rm flat}$ is the circular velocity at the flat part of the rotation curve, and match galaxy scaling relations, are able to recover the observed relations between the distributions of stars, gas and dark matter over the radial range for which we have observational rotation curve data. Furthermore, we explicitly match the observed baryonic to halo mass relation for the first time with simulated galaxies. We discuss our results in the context of the baryon cycle that is inherent in these simulations, and with regards to the effect of baryonic processes on the distribution of dark matter.