Search

We study the abundance, radial distribution, and orbits of luminous satellites in simulations of Milky Way-mass dark halos in the $\Lambda$CDM cosmology. We follow the evolution of a halo from the Aquarius project and the formation of satellites with the GALFORM semi-analytic model of galaxy formation, in which gas cools radiatively into halos before reionization and in halos that exceed a redshift-dependent ``critical'' virial mass after that. Subhalos are prone to disruption in the tidal field of the main halo, with the number of surviving self-bound subhalos increasing with resolution. Even in the highest resolution simulation (Aq-L1, with particle mass $m_{\rm p}\sim10^3\, M_\odot$), a substantial number of subhalos are disrupted but their galaxies may survive as ``orphans''. Whether or not a satellite becomes an orphan depends primarily on its time of infall. When orphans are included, the simulations yield a converged satellite stellar mass function across different resolution levels. The total number of luminous satellites is sensitive to the assumed redshift of reionization, but the shape of the satellite stellar mass function is robust, peaking at the stellar mass ($\sim 10^3\, M_\odot$) of a halo just above the critical threshold. Most orphans are found in the central regions of the main halo and make up roughly half of all satellites in Aq-L1. When orphans are taken into account there is no need to populate subhalos below the critical mass with satellites to fit the radial distribution of Milky Way satellites, as had been argued in recent work. Our model predicts that orphans dominate the ultra-faint population and that many more satellites with small apocentric radii should be detected in upcoming deep wide-field surveys., Comment: 16 pages, 10 figures, submitted to MNRAS

Kinematically-persistent planes of satellites (KPPs) are fixed sets of satellites co-orbiting around their host galaxy, whose orbital poles are conserved and clustered across long cosmic time intervals. They play the role of 'skeletons', ensuring the long-term durability of positional planes. We explore the physical processes behind their formation in terms of the dynamics of the local Cosmic Web (CW), characterized via the so-called Lagrangian Volumes (LVs) built up around two zoom-in, cosmological hydro-simulations of MW-mass disk galaxy + satellites systems, where three KPPs have been identified. By analyzing the LVs deformations in terms of the reduced Tensor of Inertia (TOI), we find an outstanding alignment between the LV principal directions and KPP satellites' orbital poles. The most compressive local mass flows (along the $\hat{e}_3$ eigenvector) are strong at early times, feeding the so-called $\hat{e}_3$-structure, while the smallest TOI axis rapidly decreases. The $\hat{e}_3$-structure collapse marks the end of this regime and is the timescale for the establishment of satellite orbital pole clustering when the Universe is $\lesssim$ 4 Gyr old. KPP proto-satellites aligned with $\hat{e}_3$ are those whose orbital poles are either aligned from early times, or have been successfully bent at $\hat{e}_3$-structure collapse. KPP satellites associated to $\hat{e}_1$ tend to have early trajectories already parallel to $\hat{e}_3$. We show that KPPs can arise as a result of the $\Lambda$CDM-predicted large-scale dynamics acting on particular sets of proto-satellites, the same dynamics that shape the local CW environment., Comment: 27 pages, 12 figures, 2 tables. Accepted for publication in The Astrophysical Journal

The Local Group (LG) of galaxies is dominated by the M31 and Milky Way (MW) pair, a configuration which suggests that the mass distribution in the LG and its surroundings should be highly anisotropic. We use the APOSTLE cosmological simulations to examine how this anisotropy manifests on the spatial distribution and kinematics of dwarf galaxies out to a distance of 3 Mpc from the MW. The simulations indicate a clear preference for dwarfs to be located close to the axis defined by the MW-M31 direction, even for dwarfs in the LG periphery (LGP; i.e., those at distances 1.25, varies with angular distance to M31, peaking in the anti-M31 direction and reaching a minimum behind M31. The combined M31-MW mass decelerates the local expansion; the LG "turnaround radius" (where =0) in APOSTLE is located at r ~ 1.25 Mpc from the LG barycentre and the pure Hubble flow (where ~ H_0*d) is not reached out to at least 3 Mpc. The predicted flow is very cold, with a barycentric dispersion of <40 km/s. A comparison of these predicted features with existing observations gives mixed results. There is clear observational evidence for an angular anisotropy in V_{rad} around the LGP, but little evidence for a preferred spatial distribution of LGP dwarfs along the MW-M31 axis. The observed local Hubble flow is also peculiar. Although the coldness of the flow is consistent with the simulations, it is significantly less decelerated: relative to the MW, on average, all galaxies beyond d~1.25 Mpc seem to be on a pure Hubble flow. We argue that these oddities may result from incompleteness and inhomogeneous sky coverage, but a full explanation may need to await the completion of deep all-sky surveys able to fill the gaps in our current inventory of nearby dwarfs., Comment: 10 pages, 5 figures, submitted to MNRAS. [Fig. 1 is a video that will play on compatible software (e.g. Okular and Adobe Acrobat, but not Preview or browser viewers).]

We use the APOSTLE suite of cosmological hydrodynamical simulations of the Local Group to examine the high speed tail of the local dark matter velocity distribution in simulated Milky Way analogues. The velocity distribution in the Solar neighborhood is well approximated by a generalized Maxwellian distribution sharply truncated at a well-defined maximum ``escape" speed. The truncated generalized Maxwellian distribution accurately models the local dark matter velocity distribution of all our Milky Way analogues, with no evidence for any separate extragalactic high-speed components. The local maximum speed is well approximated by the terminal velocity expected for particles able to reach the Solar neighborhood in a Hubble time from the farthest confines of the Local Group. This timing constraint means that the local dark matter velocity distribution is unlikely to contain any high-speed particles contributed by the Virgo Supercluster ``envelope", as argued in recent works. Particles in the Solar neighborhood with speeds close to the local maximum speed can reach well outside the virial radius of the Galaxy, and, in that sense, belong to the Local Group envelope posited in earlier work. The local manifestation of such envelope is thus not a distinct high-speed component, but rather simply the high-speed tail of the truncated Maxwellian distribution., Comment: 16 pages, 8 figures. To be submitted to JCAP

We use two zoom-in $\Lambda$CDM hydrodynamical simulations of massive disk galaxies to study the possible existence of fixed satellite groups showing a kinematically-coherent behaviour across evolution (angular momentum conservation and clustering). We identify three such groups in the two simulations, defining kinematically-coherent, time-persistent planes (KPPs) that last at least from virialization to $z=0$ (more than 7 Gyrs). This proves that orbital pole clustering is not necessarily set in at low redshift, representing a long-lived property of galaxy systems. KPPs are thin and oblate, represent $\sim25-40\%$ of the total number of satellites in the system, and are roughly perpendicular to their corresponding central disk galaxies during certain periods, consistently with Milky Way $z=0$ data. KPP satellite members are statistically distinguishable from satellites outside KPPs: they show higher specific orbital angular momenta, orbit more perpendicularly to the central disk galaxy, and have larger pericentric distances, than the latter. We numerically prove, for the first time, that KPPs and the best-quality positional planes share the same space configuration across time, such that KPPs act as `skeletons' preventing the latter of being washed out in short timescales. In one of the satellite-host systems, we witness the late capture of a massive dwarf galaxy endowed with its own satellite system, also organized into a KPP configuration prior to its capture. We briefly explore the consequences this event has on the host's KPP, and on the possible enhancement of the asymmetry in the number of satellites rotating in one sense or the opposite within the KPP., Comment: 24 pages, 8 figures, 2 tables. Accepted for publication in The Astrophysical Journal

Light-harvesting complexes in natural photosynthetic systems, such as those in purple bacteria, consist of photo-reactive chromophores embedded in densely packed "antenna" systems organized in well-defined nanostructures. In the case of purple bacteria, the chromophore antennas are composed of natural J-aggregates such as bacteriochlorophylls and carotenoids. Inspired by the molecular composition of such biological systems, we create a library of organic materials composed of densely packed J-aggregates in a polymeric matrix, in which the matrix mimics a protein scaffold. This library of organic materials shows polaritonic properties which can be tuned from the visible to the infrared by choice of the model molecule. Inspired by the molecular architecture of the light-harvesting complexes of \textit{Rhodospirillum molischianum} bacteria, we study the light-matter interactions of J-aggregate-based nanorings with similar dimensions to the analogous natural nanoscale architectures. Electromagnetic simulations show that these nanorings of J-aggregates can act as resonators, with subwavelength confinement of light while concentrating the electric field in specific regions. These results open the door to bio-inspired building blocks for all-organic metamaterials while offering a new perspective on light-matter interactions at the nanoscale in densely packed organic matter in biological organisms including photosynthetic organelles., Comment: Including supplementary information

Optical Tamm States (OTS) are confined optical modes that can occur at the interface between two highly reflective structures. However, due to the strong reflectance required, their implemen-tation with highly processable and metal-free flexible materials has proven challenging. Herein, we develop the first structure supporting OTS based only on organic polymeric materials, demon-strating a photonic platform based on non-critical, widely available, and easily processable mate-rials. The structures fabricated present large areas and consist of a narrowband multi-layered polymeric Distributed Bragg Reflector (DBR) followed by a thin film of J-aggregate molecular exci-tonic material that can act as a highly reflective surface within a narrowband range. We take ad-vantage of the narrowband spectral response of the DBR and of the reflective molecular layer to tune the OTS band by varying the periodicity of the multilayer, opening the door for the fabrica-tion of OTS structures based on lightweight integrable excitonic devices with cost-effective proce-dures., Comment: 15 pages, 8 (3+5) figures, supplementary information included

We use the APOSTLE Local Group (LG) cosmological hydro-simulations to examine the properties of "backsplash" galaxies, i.e, dwarfs which were within the virial boundaries of the Milky Way (MW) or M31 in the past, but are today outside their virial radius ($r_{200}$). More than half of all dwarfs between $1-2\,r_{200}$ of each primary are backsplash. More distant backsplash systems, i.e., those reaching distances well beyond $2\,r_{200}$, are typically close to apocentre of nearly radial orbits, and, therefore, essentially at rest relative to their primary. We use this result to investigate which LG dwarfs beyond $\sim500$ kpc of either primary could be a distant backsplash satellite of MW or M31. Tucana dSph, one of the few known quiescent LG field dwarfs, at $d_{\rm M31}\approx1350$ kpc and $d_{\rm MW}\approx880$ kpc, is a promising candidate. Tucana's radial velocity is consistent with being at rest relative to M31. Further, Tucana is located close to M33's orbital plane around M31, and simple orbit integrations indicate that Tucana may have been ejected during an early pericentric passage of M33 $\sim11$ Gyr ago, a timing which approximately coincides with Tucana's last episode of star formation. We suggest that Tucana may have been an early-infalling satellite of M31 or M33, providing a compelling explanation for its puzzling lack of gas and ongoing star formation despite its isolated nature. In this scenario, M33 should have completed some orbits around M31, a result that may help to explain the relative dearth of M33 satellite-candidates identified so far., Comment: 9 pages, 5 figures, Accepted for publication in MNRAS

We use the APOSTLE cosmological simulations to examine the role of the cosmic UV background in regulating star formation (SF) in low-mass LCDM halos. In agreement with earlier work, we find that after reionization SF proceeds mainly in halos whose mass exceeds a redshift-dependent ``critical'' mass, Mcrit, set by the structure of the halos and by the thermal pressure of UV-heated gas. Mcrit increases from ~10^8 Msun at z~10 to Mcrit ~10^9.7 Msun at z=0, roughly following the average mass growth of halos in that mass range. This implies that halos well above or below critical at present have remained so since early times. Halos of luminous dwarfs today were already above-critical and star-forming at high redshift, explaining naturally the ubiquitous presence of ancient stellar populations in dwarfs, regardless of luminosity. The SF history of systems close to the critical boundary is more complex. SF may cease or reignite in dwarfs whose host halo falls below or climbs above the critical boundary, suggesting an attractive explanation for the episodic nature of SF in some dwarfs. Also, some subcritical halos today may have been above critical in the past; these systems should at present make up a sizable population of faint field dwarfs lacking ongoing star formation. Although few such galaxies are currently known, the discovery of this population would provide strong support for our results. Our work indicates that, rather than stellar feedback, it is the ionizing UV background and mass accretion history what regulates SF in the faintest dwarfs., Comment: 12 pages, 12 figures. Accepted by MNRAS

We use $\Lambda$CDM cosmological hydrodynamical simulations to explore the kinematics of gaseous discs in late-type dwarf galaxies. We create high-resolution 21-cm 'observations' of simulated dwarfs produced in two variations of the EAGLE galaxy formation model: one where supernova-driven gas flows redistribute dark matter and form constant-density central 'cores', and another where the central 'cusps' survive intact. We 'observe' each galaxy along multiple sight lines and derive a rotation curve for each observation using a conventional tilted-ring approach to model the gas kinematics. We find that the modelling process introduces systematic discrepancies between the recovered rotation curve and the actual circular velocity curve driven primarily by (i) non-circular gas orbits within the discs; (ii) the finite thickness of gaseous discs, which leads to overlap of different radii in projection; and (iii) departures from dynamical equilibrium. Dwarfs with dark matter cusps often appear to have a core, whilst the inverse error is less common. These effects naturally reproduce an observed trend which other models struggle to explain: late-type dwarfs with more steeply-rising rotation curves appear to be dark matter-dominated in the inner regions, whereas the opposite seems to hold in galaxies with core-like rotation curves. We conclude that if similar effects affect the rotation curves of observed dwarfs, a late-type dwarf population in which all galaxies have sizeable dark matter cores is most likely incompatible with current measurements., Comment: 21 pages, 13 figures, 1 table. Accepted for publication in MNRAS. Appendices included as ancillary file

The abundance of the faintest galaxies provides insight into the nature of dark matter and the process of dwarf galaxy formation. In the LCDM scenario, low mass halos are so numerous that the efficiency of dwarf formation must decline sharply with decreasing halo mass in order to accommodate the relative scarcity of observed dwarfs and satellites in the Local Group. The nature of this decline contains important clues to the mechanisms regulating the onset of galaxy formation in the faintest systems. We explore here two possible models for the stellar mass ($M_*$)-halo mass ($M_{200}$) relation at the faint end, motivated by some of the latest LCDM cosmological hydrodynamical simulations. One model includes a sharp mass threshold below which no luminous galaxies form, as expected if galaxy formation proceeds only in systems above the Hydrogen-cooling limit. In the second model, $M_*$ scales as a steep power-law of $M_{200}$ with no explicit cutoff, as suggested by recent semianalytic work. Although both models predict satellite numbers around Milky Way-like galaxies consistent with current observations, they predict vastly different numbers of ultra-faint dwarfs and of satellites around isolated dwarf galaxies. Our results illustrate how the satellite mass function around dwarfs may be used to probe the $M_*$-$M_{200}$ relation at the faint end and to elucidate the mechanisms that determine which low-mass halos "light up" or remain dark in the LCDM scenario., Comment: 14 pages, 11 figures. Accepted version for publication in MNRAS

The dark matter (DM) content of the Fornax dwarf spheroidal galaxy inferred from its kinematics is substantially lower than expected from LCDM cosmological simulations. We use N-body simulations to examine whether this may be the result of Galactic tides. We find that, despite improved proper motions from the Gaia mission, the pericentric distance of Fornax remains poorly constrained, mainly because its largest velocity component is roughly anti-parallel to the solar motion. Translating Fornax's proper motion into a Galactocentric velocity is thus sensitively dependent on Fornax's assumed distance: the observed distance uncertainty, $\pm 8\%$, implies pericentric distances that vary between $r_{\rm peri}\sim 50$ and $r_{\rm peri}\sim 150$ kpc. Our simulations show that for $r_{\rm peri}$ in the lower range of that estimate, a LCDM subhalo with maximum circular velocity $V_{\rm max}=40$ km s$^{-1}$ (or virial mass $M_{200}\approx 10^{10} M_\odot$, as expected from LCDM) would be tidally stripped to $V_{\rm max} \sim 23$ km s$^{-1}$ over $10$ Gyr. This would reduce the DM mass within the Fornax stellar half-mass radius to about half its initial value, bringing it into agreement with observations. Tidal stripping affects mainly Fornax's DM halo; its stellar component is affected little, losing less than $5\%$ of its initial mass in the process. We also explore the effect of Galactic tides on the dynamical friction decay times of Fornax's population of globular clusters (GC) and find little evidence for substantial changes, compared with models run in isolation. A population of GCs with initial orbital radii between $1$ and $2$ kpc is consistent with the present-day spatial distribution of Fornax GCs, despite assuming a cuspy halo. Neither the DM content nor the spatial distribution of GCs seem inconsistent with a simple model where Fornax inhabits a tidally-stripped cuspy cold DM halo., Comment: 10 pages, 8 figures, 4 tables. Accepted for publication in MNRAS

We use the APOSTLE $\Lambda$CDM cosmological hydrodynamical simulations of the Local Group to study the recent accretion of massive satellites into the halo of Milky Way (MW)-sized galaxies. These systems are selected to be close analogues to the Large Magellanic Cloud (LMC), the most massive satellite of the MW. The simulations allow us to address, in a cosmological context, the impact of the Clouds on the MW, including the contribution of Magellanic satellites to the MW satellite population, and the constraints placed on the Galactic potential by the motion of the LMC. We show that LMC-like satellites are twice more common around Local Group-like primaries than around isolated halos of similar mass; these satellites come from large turnaround radii and are on highly eccentric orbits whose velocities at first pericentre are comparable with the primary's escape velocity. This implies $V_{\rm esc}^{\rm MW} (50 $ kpc$)\sim 365$ km/s, a strong constraint on Galactic potential models. LMC analogues contribute about 2 satellites with $M_*>10^5\, M_\odot$, having thus only a mild impact on the luminous satellite population of their hosts. At first pericentre, LMC-associated satellites are close to the LMC in position and velocity, and are distributed along the LMC's orbital plane. Their orbital angular momenta roughly align with the LMC's, but, interestingly, they may appear to "counter-rotate" the MW in some cases. These criteria refine earlier estimates of the LMC association of MW satellites: only the SMC, Hydrus1, Car3, Hor1, Tuc4, Ret2 and Phoenix2 are compatible with all criteria. Carina, Grus2, Hor2 and Fornax are less probable associates given their large LMC relative velocity., Comment: 17 pages, 14 figures, 1 table. Accepted for publication in MNRAS

We present a detailed characterization of planes of satellite galaxies in the Milky Way (MW) and M31. For a positional analysis, we introduce an extension to the `4-galaxy-normal density plot' method \citep[][P13]{Pawlowski13}. It finds the normal directions to the predominant planar configurations of satellites of a system, yielding for each a \textit{collection} of planes of increasing member satellites. This allows to quantify the quality of planes in terms of population ($N_{\rm sat}$) and spatial flattening ($c/a$). We apply this method to the latest data for confirmed MW and M31 satellite samples, with 46 and 34 satellites, respectively. New MW satellites form part of planes previously identified from the sample with $N_{\rm sat}=27$ studied in P13: we identify a new plane with $N_{\rm sat}=39$ as thin as the VPOS-3 ($c/a\sim 0.2$), and with roughly the same normal direction; so far the most populated plane that thin reported in the Local Group. We introduce a new method to determine, using kinematic data, the axis of maximum co-orbitation of MW satellites. Interestingly, this axis approximately coincides with the normal to the former plane: $\geq45\pm5\%$ of satellites co-orbit. In M31 we discover a plane with $N_{\rm sat}=18$ and $c/a\sim0.15$, i.e., quality comparable to the GPoA, and perpendicular to it. This structure is viewed face-on from the Sun making it susceptible to M31 satellite distance uncertainties. An estimation of the perpendicular velocity dispersion suggests it is dynamically unstable. Finally we find that mass is not a property determining a satellite's membership to good quality planes., Comment: 20 pages. 13 figures. 5 tables. Accepted for publication in MNRAS. Supercedes arXiv:1908.02298. Please cite the NEW version. This updated and extended version of the paper considers up-to-date samples of confirmed MW and M31 satellite galaxies and studies the co-orbitation of MW satellites

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., Comment: 26 pages, 13 figures, 2 tables. Accepted for publication in ApJ