In this paper, we build from previous work (Bustard et al. 2018) and present simulations of recent (within the past Gyr), magnetized, cosmic ray driven outflows from the Large Magellanic Cloud (LMC), including our first attempts to explicitly use the derived star formation history of the LMC to seed outflow generation. We run a parameter set of simulations for different LMC gas masses and cosmic ray transport treatments, and we make preliminary comparisons to published outflow flux estimates, neutral and ionized hydrogen observations, and Faraday rotation measure maps. We additionally report on the gas mass that becomes unbound from the LMC disk and swept by ram pressure into the Trailing Magellanic Stream. We find that, even for our largest outburst, the mass contribution to the Stream is still quite small, as much of the outflow-turned-halo gas is shielded on the LMC's far-side due to the LMC's primarily face-on infall through the Milky Way halo over the past Gyr. On the LMC's near-side, past outflows have fought an uphill battle against ram pressure, with near-side halo mass being at least a factor of a few smaller than the far-side. Absorption line studies probing only the LMC foreground, then, may be severely underestimating the total mass of the LMC halo formed by outflows.

The Laser Interferometer Space Antenna (LISA) will enable Galactic gravitational wave (GW) astronomy by individually resolving $ > 10^4$ signals from double white dwarf (DWD) binaries throughout the Milky Way. In this work we assess for the first time the potential of LISA data to map the Galactic stellar bar and spiral arms, since GWs are unaffected by stellar crowding and dust extinction unlike optical observations of the bulge region. To achieve this goal we combine a realistic population of Galactic DWDs with a high-resolution N-Body simulation a galaxy in good agreement with the Milky Way. We then model GW signals from our synthetic DWD population and reconstruct the structure of the simulated Galaxy from mock LISA observations. Our results show that while the low signal contrast between the background disc and the spiral arms hampers our ability to characterise the spiral structure, the stellar bar will instead clearly appear in the GW map of the bulge. The bar length and bar width derived from these synthetic observations are underestimated, respectively within $1\sigma$ and at a level greater than $2\sigma$, but the resulting axis ratio agrees to well within $1\sigma$, while the viewing angle is recovered to within one degree. These are competitive constraints compared to those from electromagnetic tracers, and they are obtained with a completely independent method. We therefore foresee that the synergistic use of GWs and electromagnetic tracers will be a powerful strategy to map the bar and the bulge of the Milky Way.

Warps are common features in both stellar and gaseous disks of nearby spiral galaxies with the latter usually easier to detect. Several theories have been proposed in the literature to explain their formation and prevalence, including tidal interactions with external galaxies. Observational correlations also suggest the importance of tides for warp formation. Here, we use the TNG100 run from the magnetohydrodynamical cosmological simulation suite IllustrisTNG to investigate the connection between interactions and the formation of gas warps. We find that in the sample of well-resolved gas-rich spiral galaxies ($10^{10}\lesssim\mathrm{M_{*}/M_{\odot}}\lesssim10^{11}$ at $z=0$) from the simulation TNG100-1, about $16\%$ possess the characteristic S-shaped warp. Around one third of these objects have their vertical morphology induced by interactions with other galaxies. Half of these interactions end with the perturber absorbed by the host. Warps formed in interactions are more asymmetrical than the remaining sample, however after the interaction the asymmetry decreases with time. We find that warps induced by interactions survive on average for $<1$ Gyr. The angle between the orbital angular momentum of the perturber and the angular momentum of the host's disk that most likely leads to warp formation is around 45 degrees. While our main goal is to investigate tidally induced warps, we find that during interactions in addition to tides, new gas that is accreted from infalling satellites also can contribute to warp formation.

We use the Auriga cosmological simulations of Milky Way (MW)-mass galaxies and their surroundings to study the satellite populations of dwarf galaxies in $\Lambda$CDM. As expected from prior work, the number of satellites above a fixed stellar mass is a strong function of the mass of the primary dwarf. For galaxies as luminous as the Large Magellanic Cloud (LMC), and for halos as massive as expected for the LMC (determined by its rotation speed), the simulations predict about 3 satellites with stellar masses exceeding $M_*>10^5\, M_\odot$. If the LMC is on its first pericentric passage, then these satellites should be near the LMC and should have orbital angular momenta roughly coincident with that of the LMC. We use 3D positions and velocities from the 2nd data release of the Gaia mission to revisit which of the "classical" MW dwarf spheroidals could plausibly be LMC satellites. The new proper motions of the Fornax and Carina dwarf spheroidals place them on orbits closely aligned with the orbital plane of the Magellanic Clouds, hinting at a potential Magellanic association. Together with the Small Magellanic Cloud (SMC), this result raises to $3$ the number of LMC satellites with $M_*>10^5\, M_\odot$, as expected from simulations. This also fills the 12-mag luminosity gap between the SMC and the ultra-faints Hyi1, Car2, Hor1, and Car3, the few ultra-faint satellites confirmed to have orbits consistent with a Magellanic origin.

We want to study how the velocity segregation and the radial profile of the velocity dispersion depend on the prominence of the brightest cluster galaxies (BCGs). We divide a sample of 102 clusters and groups of galaxies into four bins of magnitude gap between the two brightest cluster members. We then compute the velocity segregation in bins of absolute and relative magnitudes. Moreover, for each bin of magnitude gap we compute the radial profile of the velocity dispersion. When using absolute magnitudes, the segregation in velocity is limited to the two brightest bins and no significant difference is found for different magnitude gaps. However, when we use relative magnitudes, a trend appears in the brightest bin: the larger the magnitude gap, the larger the velocity segregation. We also show that this trend is mainly due to the presence, in the brightest bin, of satellite galaxies in systems with small magnitude gaps: in fact, if we study separately central galaxies and satellites, this trend is mitigated and central galaxies are more segregated than satellites for any magnitude gap. A similar result is found in the radial velocity dispersion profiles: a trend is visible in central regions (where the BCGs dominate) but, if we analyse the profile using satellites alone, the trend disappears. In the latter case, the shape of the velocity dispersion profile in the centre of systems with different magnitude gaps show three types of behaviours: systems with the smallest magnitude gaps have an almost flat profile from the centre to the external regions; systems with the largest magnitude gaps show a monothonical growth from the low values of the central part to the flat ones in the external regions; finally, systems with $1.0 < \Delta m_{12} \le 1.5$ show a profile that peaks in the centres and then decreases towards the external regions. We suggest that two mechanisms could be respons....

About 20% of stars in the solar vicinity are in the Hercules stream, a bundle of stars that move together with a velocity distinct from the Sun. Its origin is still uncertain. Here, we explore the possibility that Hercules is made of trojans, stars captured at L4, one the Lagrangian points of the stellar bar. Using GALAKOS--a high-resolution N-body simulation of the Galactic disk--we follow the motions of stars in the co-rotating frame of the bar and confirm previous studies on Hercules being formed by stars in co-rotation resonance with the bar. Unlike previous work, we demonstrate that the retrograde nature of trojan orbits causes the asymmetry in the radial velocity distribution, typical of Hercules in the solar vicinity. We show that trojans remain at capture for only a finite amount of time, before escaping L4 without being captured again. We anticipate that in the kinematic plane the Hercules stream will de-populate along the bar major axis and be visible at azimuthal angles behind the solar vicinity with a peak towards L4. This test can exclude the OLR origin of the Hercules stream and be validated by Gaia DR3 and DR4.

Recent studies have shown that the passage of a massive satellite through the disk of a spiral galaxy can induce vertical wobbles in the disk and produce features such as in-plane rings and phase-space spirals. Here we analyze a high-resolution N-body simulation of a live stellar disk perturbed by the recent passage of a massive dwarf galaxy that induces such perturbations. We study the implications of the phase-space structures for the estimate of the matter density through traditional Jeans modelling. The dwarf satellite excites rapid time-variations in the potential, leading to a significant bias of the local matter surface density determined through such a method. In particular, while the Jeans modelling gives reasonable estimates in the most overdense regions of the disk, we show that it tends to overestimate the dynamical surface density in underdense regions. In these regions, the phase-space spiral is indeed more marked in the surface of section of height vs. vertical velocity. This prediction can be verified with future Gaia data releases. Our finding is highly relevant for future attempts at determining the dynamical surface density of the outer Milky Way disk as a function of radius. The outer disk of the Milky Way is indeed heavily perturbed, and \textit{Gaia} DR2 data have clearly shown that such phase-space perturbations are even present locally. While our results show that traditional Jeans modelling should give reliable results in overdense regions of the disk, the important biases in underdense regions call for the development of non-equilibrium methods to estimate the dynamical matter density locally and in the outer disk.

We use two unsupervised machine learning algorithms, Shared Nearest Neighbor (SNN) and \textit{EnLink}, as a single approach to identify stellar groupings in the Orion star-forming complex as an application to the 5-dimensional astrometric data from Gaia DR2. The algorithms present two different ways to limit user bias when evaluating the relative weights among the astrometric parameters, automatically determined by the machine and through a standard procedure by monitoring several outcome measures. Both algorithms complement each other and produce similar stellar groups. Because SNN groups have a much smaller spread in proper motions compared to \textit{EnLink}, we use \textit{EnLink}, which requires no input, as a first pass tool for group identification and validation. We then used the SNN algorithm to dissect the Orion star-forming complex. We identify 21 spatially- and kinematically-coherent groups in the Orion complex, nine of which previously unknown. The groups show a wide distribution of distances extending as far as about 150 pc in front of the star-forming Orion clouds, to about 50 pc beyond them where we find, unexpectedly, three groups. Our results expose to view the wealth of sub-structure in the OB association, within and beyond the classical Blaauw Orion OBI sub-groups. A full characterization of the new groups is of the essence as it offers the potential to unveil how star formation proceeds globally in large complexes such as Orion. The data and code that generated the groups is provided in this Letter.

Extending for over 200 degrees across the sky, the Magellanic Stream together with its Leading Arm is the most spectacular example of a gaseous stream in the local Universe. The Stream is an interwoven tail of filaments trailing the Magellanic Clouds as they orbit the Milky Way. Thought to be created by tidal forces, ram pressure, and halo interactions, the Stream is a benchmark for dynamical models of the Magellanic System, a case study for gas accretion and dwarf-galaxy accretion onto galaxies, a probe of the outer halo, and the bearer of more gas mass than all other Galactic high velocity clouds combined. If it survives to reach the Galactic disk, it may maintain or even elevate the Galactic star-formation rate. In this white paper, we emphasize the Stream's importance for many areas of Galactic astronomy, summarize key unanswered questions, and identify future observations and simulations needed to resolve them. We stress the importance of ultraviolet, optical, and radio spectroscopy, and the need for computational models that capture full particle and radiation treatments within an MHD environment.

Fossil groups (FGs) are galaxy aggregates with an extended and luminous X-ray halo, which are dominated by a very massive early-type galaxy and lack of L* objects. FGs are indeed characterized by a large magnitude gap between their central and surrounding galaxies. This is explained by either speculating that FGs are failed groups that formed without bright satellite galaxies and did not suffer any major merger, or by suggesting that FGs are very old systems that had enough time to exhaust their bright satellite galaxies through multiple major mergers. Since major mergers leave signatures in the stellar populations of the resulting galaxy, we study the stellar population parameters of the brightest central galaxies (BCGs) of FGs as a benchmark against which the formation and evolution scenarios of FGs can be compared. We present long-slit spectroscopic observations along different axes of NGC 6482 and NGC 7556, which are the BCGs of two nearby FGs. The measurements include spatially resolved stellar kinematics and radial profiles of line-strength indices, which we converted into stellar population parameters using single stellar-population models. NGC 6482 and NGC 7556 are very massive and large galaxies and host a centrally concentrated stellar population, which is significantly younger and more metal rich than the rest of the galaxy. The age gradients of both galaxies are somewhat larger than those of the other FG BCGs studied so far, whereas their metallicity gradients are similarly negative and shallow. They have negligible gradients of alpha-element abundance ratio. The measured metallicity gradients are less steep than those predicted for massive galaxies that formed monolithically and evolved without experiencing any major merger. We conclude that the observed FGs formed through major mergers rather than being failed groups that lacked bright satellite galaxies from the beginning.

We present a new precision measurement of gas-phase abundances of S, O, N, Si, Fe, P, Al, Ca as well as molecular hydrogen (H_2) in the Leading Arm (region II, LAII) of the Magellanic Stream (MS) towards the Seyfert galaxy NGC 3783. The results are based on high-quality archival ultraviolet/optical/radio data from various different instruments (HST/STIS, FUSE, AAT, GBT, GB140ft, ATCA). Our study updates previous results from lower-resolution data and provides for the first time a self-consistent component model of the complex multi-phase absorber, delivering important constraints on the nature and origin of LAII. We derive a uniform, moderate alpha abundance in the two main absorber groups at +245 and +190 km s^-1 of alpha/H = 0.30 pm 0.05 solar, a low nitrogen abundance of N/H = 0.05 pm 0.01 solar, and a high dust content with substantial dust depletion values for Si, Fe, Al, and Ca. These alpha, N, and dust abundances in LAII are similar to those observed in the Small Magellanic Cloud (SMC). From the analysis of the H_2 absorption, we determine a high thermal pressure of P/k = 1680 K cm^-3 in LAII, in line with the idea that LAII is located in the inner Milky Way halo at a z-height of <20 kpc where it hydrodynamically interacts with the ambient hot coronal gas. Our study supports a scenario, in which LAII stems from the break-up of a metal- and dust-enriched progenitor cloud that was recently (200-500 Myr ago) stripped from the SMC.

We present a hydrodynamical model of M33 and its recent interaction with M31. This scenario was previously proposed in the literature in order to explain the distorted gaseous and stellar disks of M33, as well as the increased star formation rate in both objects around 2 Gyr ago. We used an orbit integration scheme to find which estimate of the transverse velocity of M31 more favors the interaction scenario and then tried to reproduce it in our simulations. M33 was modeled as a stellar and gaseous disk embedded in a live dark-matter halo, while M31 was approximated only with a live dark halo. In the simulations the two galaxies passed each other with a pericenter distance of 37 kpc. Tides excited a two-armed spiral structure in the M33 disk, which is found to be the predominant spiral signal in the observed galaxy and has been long known as a feature easily induced by tidal interactions. We found that the gaseous warp produced by the interaction did not resemble enough the observed one and we performed an additional simulation including the hot gas halo of M31 to show that this feature can be properly reproduced by tidal forces and ram pressure stripping acting simultaneously on the gaseous disk. In addition to the spiral arms, tidal forces produced the stellar stream similar to the observed one and triggered a star formation burst at similar radii as it is observed.

Galactic winds and fountains driven by supernova-heated gas play an integral role in re-distributing gas in galaxies, depositing metals in the circumgalactic medium (CGM), and quenching star formation. The interplay between these outflows and ram pressure stripping due to the galaxy's motion through an ambient medium may enhance these effects by converting fountain flows into expelled gas. In this paper, we present controlled, 3D simulations of ram pressure stripping combined with thermally driven, local outflows from clustered supernovae in an isolated disk galaxy modeled on the Large Magellanic Cloud (LMC), a dwarf satellite of the Milky Way on its first infall. Observational evidence of local outflows emanating from supergiant shells in the LMC and a trailing filament of HI gas originating from these regions - with no obvious Leading Arm counterpart - may represent a perfect example of this process. Our simulations present a proof-of-concept that ram pressure can convert fountain flows into expelled gas. We find that fountains launched near the peak star formation time of the LMC can comprise part of the LMC filament in the Trailing Stream, but with lower column densities than observed. Larger, more numerous outflows from the LMC may be possible and may contribute more mass, but higher inertia gas will lengthen the timescale for this gas to be swept away by ram pressure. Given the high resolution observations, increased knowledge of star formation histories, and growing evidence of multiphase, ionized outflows, the LMC is an ideal test-bed for future wind models.

We have developed a novel technique based on a clustering algorithm which searches for kinematically- and chemically-clustered stars in the APOGEE DR12 Cannon data. As compared to classical chemical tagging, the kinematic information included in our methodology allows us to identify stars that are members of known globular clusters with greater confidence. We apply our algorithm to the entire APOGEE catalog of 150,615 stars whose chemical abundances are derived by the Cannon. Our methodology found anti-correlations between the elements Al and Mg, Na and O, and C and N previously identified in the optical spectra in globular clusters, even though we omit these elements in our algorithm. Our algorithm identifies globular clusters without a priori knowledge of their locations in the sky. Thus, not only does this technique promise to discover new globular clusters, but it also allows us to identify candidate streams of kinematically- and chemically-clustered stars in the Milky Way.

One of the main theories for explaining the formation of spiral arms in galaxies is the stationary density wave theory. This theory predicts the existence of an age gradient across the arms. We use the stellar cluster catalogues of the galaxies NGC 1566, M51a, and NGC 628 from the Legacy Extragalactic UV Survey (LEGUS) program. In order to test for the possible existence of an age sequence across the spiral arms, we quantified the azimuthal offset between star clusters of different ages in our target galaxies. We found that NGC 1566, a grand-design spiral galaxy with bisymmetric arms and a strong bar, shows a significant age gradient across the spiral arms that appears to be consistent with the prediction of the stationary density wave theory. In contrast, M51a with its two well-defined spiral arms and a weaker bar does not show an age gradient across the arms. In addition, a comparison with non LEGUS star cluster catalogues for M51a yields similar results. We believe that the spiral structure of M51a is not the result of a stationary density wave with a fixed pattern speed. Instead, tidal interactions could be the dominant mechanism for the formation of spiral arms. We also found no offset in the azimuthal distribution of star clusters with different ages across the weak spiral arms of NGC 628.

The Magellanic Stream and Leading Arm of HI that stretches from the Large and Small Magellanic Clouds (LMC and SMC) and over 200 degrees of the Southern sky is thought to be formed from multiple encounters between the LMC and SMC. In this scenario, most of the gas in the Stream and Leading Arm is stripped from the SMC, yet recent observations have shown a bifurcation of the Trailing Arm that reveals LMC origins for some of the gas. Absorption measurements in the Stream also reveal an order of magnitude more gas than in current tidal models. We present hydrodynamical simulations of the multiple encounters between the LMC and SMC at their first pass around the Milky Way, assuming that the Clouds were more extended and gas rich in the past. Our models create filamentary structures of gas in the Trailing Stream from both the LMC and SMC. While the SMC trailing filament matches the observed Stream location, the LMC filament is offset. In addition, the total observed mass of the Stream in these models is underestimated of a factor of four when the ionized component is accounted for. Our results suggest that there should also be gas stripped from both the LMC and SMC in the Leading Arm, mirroring the bifurcation in the Trailing Stream. This prediction is consistent with recent measurements of spatial variation in chemical abundances in the Leading Arm, which show that gas from multiple sources is present, although the nature is still uncertain.

The Leading Arm (LA) of the Magellanic Stream is a vast debris field of H I clouds connecting the Milky Way and the Magellanic Clouds. It represents an example of active gas accretion onto the Galaxy. Previously only one chemical abundance measurement had been made in the LA. Here we present chemical abundance measurements using Hubble Space Telescope/Cosmic Origins Spectrograph Green Bank Telescope spectra of four sightlines passing through the LA, and three nearby sightlines that may trace outer fragments of the LA. We find low oxygen abundances, ranging from 4.0(+4.0,-2.0) percent solar to 12.6(+6.2,-4.1) percent solar, in the confirmed LA directions, with the lowest values found in the region known as LA III, farthest from the LMC. These abundances are substantially lower than the single previous measurement, S/H=35+/-7 percent solar (Lu et al. 1998), but are in agreement with those reported in the SMC filament of the trailing Stream, supporting a common origin in the SMC (not the LMC) for the majority of the LA and the trailing Stream. This provides important constraints for models of the formation of the Magellanic System. Finally, the HVCs in two of the three nearby sightlines show H I columns, kinematics, and oxygen abundances consistent with LA membership. This suggests that the LA is larger than traditionally thought, extending at least 20 degrees further to the Galactic northwest.

It is thought that fossil systems are relics of structure formation in the primitive Universe. They are galaxy aggregations that have assembled their mass at high redshift with few or no subsequent accretion. Observationally these systems are selected by large magnitude gaps between their 1st and 2nd ranked galaxies. Nevertheless, there is still debate over whether or not this observational criterium selects dynamically evolved ancient systems. We have studied the properties of the nearby fossil group RXJ075243.6+455653 in order to understand the mass assembly of this system. Deep spectroscopic observations allow us to construct the galaxy luminosity function (LF) of RXJ075243.6+455653 down to M*+ 6. The analysis of the faint-end of the LF in groups and clusters provides valuable information about the mass assembly of the system. In addition, we have analyzed the nearby large-scale structure around this group. We identified 26 group members within r200=0.9 Mpc. The LF of the group shows a flat faint-end slope ( -1.08 +/- 0.33). This low density of dwarf galaxies is confirmed by the low value of the dwarf-to-giant ratio (DGR = 0.99 +/- 0.49) for this system. Both the lack of dwarf galaxies and the low luminosity of the BGG suggests that RXJ075243.6+455653 still has to accrete mass from its nearby environment. This mass accretion will be achieved because it is the dominant structure of a rich environment formed by several groups of galaxies (15) within 7 Mpc from the group center and with +/- 1000$ km/s. RXJ075243.6+455653 is a group of galaxies that has not yet completed the process of its mass assembly. This new mass accretion will change the fossil state of the group. This group is an example of a galaxy aggregation selected by a large magnitude gap but still in the process of the accretion of its mass (Abridged).

Using the Illustris cosmological simulation, we investigate the origin of fossil groups in the $M_{200}=10^{13}-10^{13.5}M_{\odot}/h$ mass regime. We examine the formation of the two primary features of fossil groups: the large magnitude gap between their two brightest galaxies, and their exceptionally luminous brightest group galaxy (BGG). For fossils and non-fossils identified at $z=0$, we find no difference in their halo mass assembly at early times, departing from previous studies. However, we do find a significant difference in the recent accretion history of fossil and non-fossil halos; in particular, fossil groups show a lack of recent accretion and have in majority assembled 80\% of their $M_{200}(z=0)$ mass before $z\sim 0.4$. For fossils, massive satellite galaxies accreted during this period have enough time to merge with the BGG by the present day, producing a more massive central galaxy; and, the lack of recent group accretion prevents replenishment of the bright satellite population, allowing for a large magnitude gap to develop within the past few Gyr. We thus find that the origin of the magnitude gap and overmassive BGG of fossils in Illustris depends on the recent accretion history of the groups and merger history of the BGG after their collapse at $z\sim1$. This indicates that selecting galaxy groups by their magnitude gap does not guarantee obtaining neither early-forming galaxy systems nor undisturbed central galaxies.

RoadMapping is a dynamical modeling machinery developed to constrain the Milky Way's (MW) gravitational potential by simultaneously fitting an axisymmetric parametrized potential and an action-based orbit distribution function (DF) to discrete 6D phase-space measurements of stars in the Galactic disk. In this work we demonstrate RoadMapping's robustness in the presence of spiral arms by modeling data drawn from an N-body simulation snapshot of a disk-dominated galaxy of MW mass with strong spiral arms (but no bar), exploring survey volumes with radii 500pc<=r_max<=5kpc. The potential constraints are very robust, even though we use a simple action-based DF, the quasi-isothermal DF (qDF). The best-fit RoadMapping model always recovers the correct gravitational forces where most of the stars that entered the analysis are located, even for small volumes. For data from large survey volumes, RoadMapping finds axisymmetric models that average well over the spiral arms. Unsurprisingly, the models are slightly biased by the excess of stars in the spiral arms. Gravitational potential models derived from survey volumes with at least r_max=3kpc can be reliably extrapolated to larger volumes. However, a large radial survey extent, r_max~5kpc, is needed to correctly recover the halo scale length. In general, the recovery and extrapolability of potentials inferred from data sets which were drawn from inter-arm regions appear to be better than those of data sets drawn from spiral arms. Our analysis implies that building axisymmetric models for the Galaxy with upcoming Gaia data will lead to sensible and robust approximations of the MW's potential.

Magellanic spiral galaxies are a class of one-armed systems that often exhibit an offset stellar bar, and are rarely found around massive spiral galaxies. Using a set of N-body and hydrodynamic simulations we consider a dwarf-dwarf galaxy interaction as the driving mechanism for the formation of this peculiar class of systems. We investigate here the relation between the dynamical, stellar and gaseous disk center and the bar. In all our simulations the bar center always coincides with the dynamical center, while the stellar disk becomes highly asymmetric during the encounter causing the photometric center of the Magellanic galaxy disk to become mismatched with both the bar and the dynamical center. The disk asymmetries persist for almost 2 Gyrs, the time that it takes for the disk to be re-centered with the bar, and well after the companion has passed. This explains the nature of the offset bar found in many Magellanic-type galaxies, including the Large Magellanic Cloud (LMC) and NGC 3906. In particular, these results, once applied to the LMC, suggest that the dynamical center should reside in the bar center instead of the HI center as previously assumed, pointing to a variation in the current estimate of the north component of the LMC proper motion.

We employ high-resolution N-body simulations of isolated spiral galaxy models, from low-amplitude, multi-armed galaxies to Milky Way-like disks, to estimate the vertical action of ensembles of stars in an axisymmetric potential. In the multi-armed galaxy the low-amplitude arms represent tiny perturbations of the potential, hence the vertical action for a set of stars is conserved, although after several orbital periods of revolution the conservation degrades significantly. For a Milky Way-like galaxy with vigorous spiral activity and the formation of a bar, our results show that the potential is far from steady, implying that the action is not a constant of motion. Furthermore, because of the presence of high-amplitude arms and the bar, considerable in-plane and vertical heating occurs that forces stars to deviate from near-circular orbits, reducing the degree at which the actions are conserved for individual stars, in agreement with previous results, but also for ensembles of stars. If confirmed, this result has several implications, including the assertion that the thick disk of our Galaxy forms by radial migration of stars, under the assumption of the conservation of the action describing the vertical motion of stars.

We use numerical simulations to examine the effects of radial migration on the vertical structure of galaxy disks. The simulations follow three exponential disks of different mass but similar circular velocity, radial scalelength, and (constant) scale height. The disks develop different non-axisymmetric patterns, ranging from feeble, long-lived multiple arms to strong, rapidly-evolving few-armed spirals. These fluctuations induce radial migration through secular changes in the angular momentum of disk particles, mixing the disk radially and blurring pre-existing gradients. Migration affects primarily stars with small vertical excursions, regardless of spiral pattern. This "provenance bias" largely determines the vertical structure of migrating stars: inward migrators thin down as they move in, whereas outward migrators do not thicken up but rather preserve the disk scale height at destination. Migrators of equal birth radius thus develop a strong scale-height gradient, not by flaring out as commonly assumed, but by thinning down as they spread inward. Similar gradients have been observed for low-[$\alpha$/Fe] mono-abundance populations (MAPs) in the Galaxy but our results argue against interpreting them as a consequence of radial migration. This is because outward migration does not lead to thickening, implying that the maximum scale height of any population should reflect its value at birth. In contrast, Galactic MAPs have scale heights that increase monotonically outwards, reaching values that greatly exceed those at their presumed birth radii. Given the strong vertical bias affecting migration, a proper assessment of the importance of radial migration in the Galaxy should take carefully into account the strong radial dependence of the scale heights of the various stellar populations.

We use a set of high-resolution N-body simulations of the Galactic disk to study its interactions with the population of satellites predicted cosmologically. One simulation illustrates that multiple passages of massive satellites with different velocities through the disk generate a wobble, having the appearance of rings in face-on projections of the stellar disk. They also produce flares in the disk outer parts and gradually heat the disk through bending waves. A different numerical experiment shows that an individual satellite as massive as the Sagittarius dwarf galaxy passing through the disk will drive coupled horizontal and vertical oscillations of stars in underdense regions, with small significant associated heating. This experiment shows that vertical excursions of stars in these low-density regions can exceed 1 kpc in the Solar neighborhood, resembling the coherent vertical oscillations recently detected locally. They can also induce non-zero vertical streaming motions as large as 10-20 km s$^{-1}$, consistent with recent observations in the Galactic disk. This phenomenon appears as a local ring, with modest associated disk heating.

We present a versatile family of model galactic outflows including non-uniform mass and energy source distributions, a gravitational potential from an extended mass source, and radiative losses. The model easily produces steady-state wind solutions for a range of mass-loading factors, energy-loading factors, galaxy mass and galaxy radius. We find that, with radiative losses included, highly mass-loaded winds must be driven at high central temperatures, whereas low mass-loaded winds can be driven at low temperatures just above the peak of the cooling curve, meaning radiative losses can drastically affect the wind solution even for low mass-loading factors. By including radiative losses, we are able to show that subsonic flows can be ignored as a possible mechanism for expelling mass and energy from a galaxy compared to the more efficient transonic solutions. Specifically, the transonic solutions with low mass-loading and high energy-loading are the most efficient. Our model also produces low-temperature, high-velocity winds that could explain the prevalence of low-temperature material in observed outflows. Finally, we show that our model, unlike the well-known Chevalier & Clegg (1985) model, can reproduce the observed linear relationship between wind X-ray luminosity and star formation rate (SFR) over a large range of SFR from $1-1000$ M$_{\odot}$/yr assuming the wind mass-loading factor is higher for low-mass, and hence, low-SFR galaxies. We also constrain the allowed mass-loading factors that can fit the observed X-ray luminosity vs. SFR trend, further suggesting an inverse relationship between mass-loading and SFR as explored in advanced numerical simulations.

The Magellanic Clouds are surrounded by an extended network of gaseous structures. Chief among these is the Magellanic Stream, an interwoven tail of filaments trailing the Clouds in their orbit around the Milky Way. When considered in tandem with its Leading Arm, the Stream stretches over 200 degrees on the sky. Thought to represent the result of tidal interactions between the Clouds and ram-pressure forces exerted by the Galactic corona, its kinematic properties reflect the dynamical history of the closest pair of dwarf galaxies to the Milky Way. The Stream is a benchmark for hydrodynamical simulations of accreting gas and cloud/corona interactions. If the Stream survives these interactions and arrives safely in the Galactic disk, its cargo of over a billion solar masses of gas has the potential to maintain or elevate the Galactic star formation rate. In this article, we review the current state of knowledge of the Stream, including its chemical composition, physical conditions, origin, and fate. We also review the dynamics of the Magellanic System, including the proper motions and orbital history of the Large and Small Magellanic Clouds, the first-passage and second-passage scenarios, and the evidence for a Magellanic Group of galaxies.

Fossil groups are expected to be the final product of galaxy merging within galaxy groups. In simulations, they are predicted to assemble their mass at high redshift. This early formation allows for the innermost $M^\ast$ galaxies to merge into a massive central galaxy. Then, they are expected to maintain their fossil status because of the few interactions with the large-scale structure. In this context, the magnitude gap between the two brightest galaxies of the system is considered a good indicator of its dynamical status. As a consequence, the systems with the largest gaps should be dynamically relaxed. In order to examine the dynamical status of these systems, we systematically analyze, for the first time, the presence of galaxy substructures in a sample of 12 spectroscopically-confirmed fossil systems with redshift $z \le 0.25$. We apply a number of tests in order to investigate the substructure in fossil systems in the two-dimensional space of projected positions out to $R_{200}$. Moreover, for a subsample of 5 systems with at least 30 spectroscopically-confirmed members we also analyze the substructure in the velocity and in the three-dimensional velocity-position spaces. Additionally, we look for signs of recent mergers in the regions around the central galaxies. We find that an important fraction of fossil systems show substructure. The fraction depends critically on the adopted test, since each test is more sensible to a particular type of substructure. Our interpretation of the results is that fossil systems are not, in general, as relaxed as expected from simulations. Our sample of 12 spectroscopically-confirmed fossil systems need to be extended in order to compute an accurate fraction, but our conclusion is that it is similar to the fraction of substructure detected in non-fossil clusters. THIS ABSTRACT IS TRUNCATED.

We present the first pointed X-ray observations of 10 candidate fossil galaxy groups and clusters. With these Suzaku observations, we determine global temperatures and bolometric X-ray luminosities of the intracluster medium (ICM) out to $r_{500}$ for six systems in our sample. The remaining four systems show signs of significant contamination from non-ICM sources. For the six objects with successfully determined $r_{500}$ properties, we measure global temperatures in the range $2.8 \leq T_{\mathrm{X}} \leq 5.3 \ \mathrm{keV}$, bolometric X-ray luminosities of $0.8 \times 10^{44} \ \leq L_{\mathrm{X,bol}} \leq 7.7\times 10^{44} \ \mathrm{erg} \ \mathrm{s}^{-1}$, and estimate masses, as derived from $T_{\mathrm{X}}$, of $M_{500} > 10^{14} \ \mathrm{M}_{\odot}$. Fossil cluster scaling relations are constructed for a sample that combines our Suzaku observed fossils with fossils in the literature. Using measurements of global X-ray luminosity, temperature, optical luminosity, and velocity dispersion, scaling relations for the fossil sample are then compared with a control sample of non-fossil systems. We find the fits of our fossil cluster scaling relations are consistent with the relations for normal groups and clusters, indicating fossil clusters have global ICM X-ray properties similar to those of comparable mass non-fossil systems.

We use test-particle integrations to show that epicyclic motions excited by a pericentre passage of a dwarf galaxy could account for bulk vertical velocity streaming motions recently observed in the Galactic stellar disc near the Sun. We use fixed potential test-particle integrations to isolate the role of phase wrapping of epicyclic perturbations from bending and breathing waves or modes, which require self-gravity to oscillate. Perturbations from a fairly massive Sagittarius dwarf galaxy, $M_d \sim 2.5 \times 10^{10} M_\odot$, are required to account for the size of the observed streaming motions from its orbital pericentre approximately a Gyr ago. A previous passage of the dwarf through the Galactic disc approximately 2.2 Gyr ago (with a then more massive dwarf galaxy) is less effective. If phase wrapping of epicyclic perturbations is responsible for stellar streaming motions in the Galactic disc, then there should be variations in velocity gradients on scales of a few kpc in the vicinity of the Sun.

Using N-body simulations we study the tidal evolution of initially disky dwarf galaxies orbiting a Milky Way-like host, a process known to lead to the formation of dwarf spheroidal galaxies. We focus on the effect of the orientation of the dwarf galaxy disk's angular momentum with respect to the orbital one and find very strong dependence of the evolution on this parameter. We consider four different orientations: the exactly prograde, the exactly retrograde and two intermediate ones. Tidal evolution is strongest for the exactly prograde and weakest for the exactly retrograde orbit. In the prograde case the stellar component forms a strong bar and remains prolate until the end of the simulation, while its rotation is very quickly replaced by random motions of the stars. In the retrograde case the dwarf remains oblate, does not form a bar and loses rotation very slowly. In the two cases of intermediate orientation of the disk, the evolution is between the two extremes, suggesting a monotonic dependence on the inclination. We interpret the results in terms of the resonance between the angular velocity of the stars in the dwarf and its orbital motion by comparing the measurements from simulations to semi-analytic predictions. We conclude that resonant effects are the most important mechanism underlying the tidal evolution of disky dwarf galaxies.

A model based on disk-stability criteria to determine the number of spiral arms of a general disk galaxy with an exponential disk, a bulge and a dark halo described by a Hernquist model is presented. The multifold rotational symmetry of the spiral structure can be evaluated analytically once the structural properties of a galaxy, such as the circular speed curve, and the disk surface brightness, are known. By changing the disk mass, these models are aimed at varying the critical length scale parameter of the disk and lead to a different spiral morphology in agreement with prior models. Previous studies based on the swing amplification and disk stability have been applied to constrain the mass-to-light ratio in disk galaxies. This formalism provides an analytic expression to estimate the number of arms expected by swing amplification making its application straight-forward to large surveys. It can be applied to predict the number of arms in the Milky Way as a function of radius and to constrain the mass-to-light ratio in disk galaxies for which photometric and kinematic measurements are available, like in the DiskMass survey. Hence, the halo contribution to the total mass in the inner parts of disk galaxies can be inferred in light of the ongoing and forthcoming surveys.

We use mid-infrared 3.6 and 4.5microns imaging of NGC 3906 from the Spitzer Survey of Stellar Structure in Galaxies (S4G) to understand the nature of an unusual offset between its stellar bar and the photometric center of an otherwise regular, circular outer stellar disk. We measure an offset of ~720 pc between the center of the stellar bar and photometric center of the stellar disk; the bar center coincides with the kinematic center of the disk determined from previous HI observations. Although the undisturbed shape of the disk suggests that NGC 3906 has not undergone a significant merger event in its recent history, the most plausible explanation for the observed offset is an interaction. Given the relatively isolated nature of NGC 3906 this interaction could be with dark matter sub structure in the galaxy's halo or from a recent interaction with a fast moving neighbor which remains to be identified. Simulations aimed at reproducing the observed offset between the stellar bar / kinematic center of the system and the photometric center of the disk are necessary to confirm this hypothesis and constrain the interaction history of the galaxy.

In nature we observe galaxy aggregations that span a wide range of magnitude gaps between the two first-ranked galaxies of a system ($\Delta m_{12}$). There are systems with gaps close to zero (e.g., the Coma cluster), and at the other extreme of the distribution, the largest gaps are found among the so-called fossil systems. Fossil and non-fossil systems could have different galaxy populations that should be reflected in their luminosity functions. In this work we study, for the first time, the dependence of the luminosity function parameters on $\Delta m_{12}$ using data obtained by the fossil group origins (FOGO) project. We constructed a hybrid luminosity function for 102 groups and clusters at $z \le 0.25$. We stacked all the individual luminosity functions, dividing them into bins of $\Delta m_{12}$, and studied their best-fit Schechter parameters. We additionally computed a relative luminosity function, expressed as a function of the central galaxy luminosity, which boosts our capacity to detect differences, especially at the bright end. We find trends as a function of $\Delta m_{12}$ at both the bright and faint ends of the luminosity function. In particular, at the bright end, the larger the magnitude gap, the fainter the characteristic magnitude $M^\ast$. We also find differences at the faint end. In this region, the larger the gap, the flatter the faint-end slope $\alpha$. The differences found at the bright end support a dissipationless, dynamical friction-driven merging model for the growth of the central galaxy in group- and cluster-sized halos. The differences in the faint end cannot be explained by this mechanism. Other processes, such as enhanced tidal disruption due to early infall and/or prevalence of eccentric orbits, may play a role. However, a larger sample of systems with $\Delta m_{12} > 1.5$ is needed to establish the differences at the faint end.

This is a brief rebuttal to arXiv:1502.03821, which claims to provide the first observational proof of dark matter interior to the solar circle. We point out that this result is not new, and can be traced back at least a quarter century.

We report the discovery of a giant stellar tidal stream in the halo of NGC 4631, a nearby edge-on spiral galaxy interacting with the spiral NGC 4656, in deep images taken with a 40-cm aperture robotic telescope. The stream has two components: a bridge-like feature extended between NGC 4631 and NGC 4656 (stream_SE) and an overdensity with extended features on the opposite side of the NGC 4631 disk (stream_NW). Together, these features extend more than 85 kpc and display a clear (g-r) colour gradient. The orientation of stream_SE relative to the orientations of NGC 4631 and NGC 4656 is not consistent with an origin from interaction between these two spirals, and is more likely debris from a satellite encounter. The stellar tidal features can be qualitatively reproduced in an N-body model of the tidal disruption of a single, massive dwarf satellite on a moderately eccentric orbit (e=0.6) around NGC 4631 over $\sim$ 3.5 Gyr, with a dynamical mass ratio (m1:m2) of ~40. Both modelling and inferences from the morphology of the streams indicate these are not associated with the complex HI tidal features observed between both spirals, which likely originate from a more recent, gas-rich accretion event. The detailed structure of stream_NW suggests it may contain the progenitor of the stream, in agreement with the N-body model. In addition, stream_NW is roughly aligned with two very faint dwarf spheroidal candidates. The system of dwarf galaxies and the tidal stream around NGC 4631 can provide an additional interesting case for exploring the anisotropy distribution of satellite galaxies recently reported in Local Group spiral galaxies by means of future follow-up observations.

We study the radial migration of stars driven by recurring multi-arm spiral features in an exponential disk embedded in a dark matter halo. The spiral perturbations redistribute angular momentum within the disk and lead to substantial radial displacements of individual stars, in a manner that largely preserves the circularity of their orbits and that results, after 5 Gyr (~40 full rotations at the disk scalelength), in little radial heating and no appreciable changes to the vertical or radial structure of the disk. Our results clarify a number of issues related to the spatial distribution and kinematics of migrators. In particular, we find that migrators are a heavily biased subset of stars with preferentially low vertical velocity dispersions. This "provenance bias" for migrators is not surprising in hindsight, for stars with small vertical excursions spend more time near the disk plane and thus respond more readily to non-axisymmetric perturbations. We also find that the vertical velocity dispersion of outward migrators always decreases, whereas the opposite holds for inward migrators. To first order, newly arrived migrators simply replace stars that have migrated off to other radii, thus inheriting the vertical bias of the latter. Extreme migrators might therefore be recognized, if present, by the unexpectedly small amplitude of their vertical excursions. Our results show that migration, understood as changes in angular momentum that preserve circularity, can affect strongly the thin disk, but cast doubts on models that envision the Galactic thick disk as a relic of radial migration.

This study is part of the FOssil Groups Origin (FOGO) project which aims at carrying out a systematic and multiwavelength study of a large sample of fossil systems. Here we focus on the relation between the optical luminosity (Lopt) and X-ray luminosity (Lx). Out of a sample of 28 candidate fossil systems, we consider a sample of 12 systems whose fossil classification has been confirmed by a companion study. They are compared with the complementary sample of 16 systems whose fossil nature is not confirmed and with a subsample of 102 galaxy systems from the RASS-SDSS galaxy cluster survey. Fossil and normal systems span the same redshift range 0<z<0.5 and have the same Lx distribution. For each fossil system, the Lx in the 0.1-2.4 keV band is computed using data from the ROSAT All Sky Survey. For each fossil and normal system we homogeneously compute Lopt in the r-band within the characteristic cluster radius, using data from the SDSS DR7. We sample the Lx-Lopt relation over two orders of magnitude in Lx. Our analysis shows that fossil systems are not statistically distinguishable from the normal systems both through the 2D KS test and the fit of the Lx-Lopt relation. The optical luminosity of the galaxy system does strongly correlate with the X-ray luminosity of the hot gas component, independently of whether the system is fossil or not. We conclude that our results are consistent with the classical "merging scenario" of the brightest galaxy formed via merger/cannibalism of other group galaxies, with conservation of the optical light. We find no evidence for a peculiar state of the hot intracluster medium.

(Abridged) Fossil systems are group- or cluster-sized objects whose luminosity is dominated by a very massive central galaxy. In the current cold dark matter scenario, these objects formed hierarchically at an early epoch of the Universe and then slowly evolved until present day. That is the reason why they are called {\it fossils}. We started an extensive observational program to characterize a sample of 34 fossil group candidates spanning a broad range of physical properties. Deep $r-$band images were taken for each candidate and optical spectroscopic observations were obtained for $\sim$ 1200 galaxies. This new dataset was completed with SDSS DR7 archival data to obtain robust cluster membership and global properties of each fossil group candidate. For each system, we recomputed the magnitude gaps between the two brightest galaxies ($\Delta m_{12}$) and the first and fourth ranked galaxies ($\Delta m_{14}$) within 0.5 $R_{{\rm 200}}$. We consider fossil systems those with $\Delta m_{12} \ge 2$ mag or $\Delta m_{14} \ge 2.5$ mag within the errors. We find that 15 candidates turned out to be fossil systems. Their observational properties agree with those of non-fossil systems. Both follow the same correlations, but fossils are always extreme cases. In particular, they host the brightest central galaxies and the fraction of total galaxy light enclosed in the central galaxy is larger in fossil than in non-fossil systems. Finally, we confirm the existence of genuine fossil clusters. Combining our results with others in the literature, we favor the merging scenario in which fossil systems formed due to mergers of $L^\ast$ galaxies. The large magnitude gap is a consequence of the extreme merger ratio within fossil systems and therefore it is an evolutionary effect. Moreover, we suggest that at least one candidate in our sample could represent a transitional fossil stage.

We have identified a low surface brightness stellar stream from visual inspection of SDSS imaging for the edge-on, spiral galaxy NGC5387. A blue overdensity was also identified in SDSS coincident with the stream intersection with the NGC5387 disk. The overdensity was also detected in the GALEX Deep Imaging Survey and found to contribute 38% of the total FUV integrated flux from NGC5387, which suggests that the region is actively forming stars. Deeper imaging was acquired with the Vatican Advanced Technology Telescope (VATT) in the B, V, and R filters that confirmed the presence of both the stellar stream and the blue overdensity. Analysis of the VATT photometry indicates the stellar stream is red in color, B-V = 0.7, and has a stellar mass of 6x10^8 M$_{\odot}$, which implies a merger ratio of 1:50. Assessment of the stream morphology suggests that the merger event has a circular radius, R~16 kpc, the stream formed in ~400 Myr, and the progenitor had a total mass of ~2x10 M$_{\odot}$. Spectroscopy from LBT+MODS1 was used to determine that the blue overdensity is at the same redshift as NGC5387, consists of young stellar populations (~10 Myr), is metal-poor (12 + log(O/H) = 8.03), and forming stars at an enhanced rate (~1-3 M$_{\odot}$/year) given its total stellar mass (2x10^7 M$_{\odot}$). Several interpretations are posited to explain these observational data, of which the most likely are (i) that the blue overdensity is a region of enhanced star formation in the outer disk of NGC5387 induced by the minor accretion event, and (ii) that the blue overdensity is the progenitor of the stellar stream undergoing a period of enhanced star formation as a result of its interaction with NGC5387. Confirmation and theoretical exploration of these scenarios are presented in a companion paper.

Recent observations have shown the presence of extra-solar planets in Galactic open stellar clusters, as in the Praesepe (M44). These systems provide a favorable environment for planetary formation due to the high heavy-element content exhibited by the majority of their population. The large stellar density, and corresponding high close-encounter event rate, may induce strong perturbations of planetary orbits with large semimajor axes. Here we present a set of N-body simulations implementing a novel scheme to treat the tidal effects of external stellar perturbers on planetary orbit eccentricity and inclination. By simulating five nearby open clusters we determine the rate of occurrence of bodies extracted from their parent stellar system by quasi-impulsive tidal interactions. We find that the specific free-floating planet production rate (total number of free-floating planets per unit of time, normalized by the total number of stars) is proportional to the stellar density of the cluster, with a constant of proportionality equal to (23 +/- 5)10^-6 pc^3 Myr^-1. For the Pleiades (M45) we predict that about 26% of stars should have lost their planets. This raises the exciting possibility of directly observing these wandering planets with the James Webb Space Telescope in the NIR band. Assuming a surface temperature of the planet of 500 K, a free-floating planet of Jupiter size inside the Pleiades would have a specific flux @4.4 micron of approximately 400 nJy, which would lead to a very clear detection (S/N of order 100) in only one hour of integration.

The causes of spiral structure in galaxies remain uncertain. Leaving aside the grand bisymmetric spirals with their own well-known complications, here we consider the possibility that multi-armed spiral features originate from density inhomogeneities orbiting within disks. Using high-resolution N-body simulations, we follow the motions of stars under the influence of gravity, and show that mass concentrations with properties similar to those of giant molecular clouds can induce the development of spiral arms through a process termed swing amplification. However, unlike in earlier work, we demonstrate that the eventual response of the disk can be highly non-linear, significantly modifying the formation and longevity of the resulting patterns. Contrary to expectations, ragged spiral structures can thus survive at least in a statistical sense long after the original perturbing influence has been removed.

(Abridged) Fossil systems are galaxy associations dominated by a relatively isolated, bright elliptical galaxy, surrounded by a group of smaller galaxies lacking L* objects. We analyzed the near-infrared photometric and structural properties of a sample of 20 BGGs present in FGs in order to better understand their formation mechanisms. Their surface-brightness distribution was fitted to a Sersic profile using the GASP2D algorithm. Then, the standard scaling relations were derived for the first time for these galaxies and compared with those of normal ellipticals and brightest cluster galaxies in non-fossil systems. The BGGs presented in this study represent a subset of the most massive galaxies in the Universe. We found that their ellipticity profiles are continuously increasing with the galactocentric radius. Our fossil BCGs follow closely the fundamental plane described by normal ellipticals. However, they depart from both the log \sigma_0 vs. log L_{K_{s}} and log r_{\rm e} vs. log L_{K_{s}} relations described by intermediate mass ellipticals. This occurs in the sense that our BGGs have larger effective radii and smaller velocity dispersions than those predicted by these relations. We also found that more elliptical galaxies systematically deviate from the previous relations while more rounder object do not. No similar correlation was found with the Sersic index. The derived scaling relations can be interpreted in terms of the formation scenario of the BGGs. Because our BGGs follow the fundamental plane tilt but they have larger effective radii than expected for intermediate mass ellipticals, we suggest that they only went through dissipational mergers in a early stage of their evolution and then assembled the bulk of their mass through subsequent dry mergers, contrary to previous claims that BGGs in FGs were formed mainly by the merging of gas-rich galaxies.

We report on 2 fossil groups of galaxies at z=0.425 and 0.372 discovered in the Cosmic Evolution Survey (COSMOS) area. Selected as X-ray extended sources, they have total masses (M_200) of 1.9(+/-0.41)E13 and 9.5(+/-0.42)E13 M_sun, respectively, as obtained from a recent X-ray luminosity-mass scaling relation. The lower mass system appears isolated, whereas the other sits in a well-known large-scale structure (LSS) populated by 27 other X-ray emitting groups. The identification as fossil is based on the i-band photometry of all the galaxies with a photo-z consistent with that of the group at the 2-sigma confidence level and within a projected group-centric distance equal to 0.5R_200, and i_AB<=22.5-mag limited spectroscopy. Both fossil groups exhibit high stellar-to-total mass ratios compared to all the X-ray selected groups of similar mass at 0.3<=z<=0.5 in the COSMOS. At variance with the composite galaxy stellar mass functions (GSMFs) of similarly massive systems, both fossil group GSMFs are dominated by passively evolving galaxies down to M^stars~1E10 M_sun (according to the galaxy broad-band spectral energy distributions). The relative lack of star-forming galaxies with 1E10<=M^stars<=1E11 M_sun is confirmed by the galaxy distribution in the b-r vs i color-magnitude diagram. Hence, the 2 fossil groups appear as more mature than the coeval, similarly massive groups. Their overall star formation activity ended rapidly after an accelerated build up of the total stellar mass; no significant infall of galaxies with M^stars>=1E10 M_sun took place in the last 3 to 6 Gyr. This similarity holds although the 2 fossil groups are embedded in two very different density environments of the LSS, which suggests that their galaxy populations were shaped by processes that do not depend on the LSS. However, their progenitors may do so. ...

The most accepted scenario for the origin of fossil groups (FGs) is that they are galaxy associations in which the merging rate was fast and efficient. These systems have assembled half of their mass at early epoch of the Universe, subsequently growing by minor mergers. They could contain a fossil record of the galaxy structure formation. We have started a project in order to characterize a large sample of FGs. In this paper we present the analysis of the fossil system RX J105453.3+552102. Optical deep images were used for studying the properties of the brightest group galaxy and for computing the photometric luminosity function of the group. We have also performed a detail dynamical analysis of the system based on redshift data for 116 galaxies. This galaxy system is located at z=0.47, and shows a quite large line-of-sight velocity dispersion \sigma_{v}~1000 km/s. Assuming the dynamical equilibrium, we estimated a virial mass of M ~ 10^{15} h_{70} M_{\odot}. No evidence of substructure was found within 1.4 Mpc radius. We found a statistically significant departure from Gaussianity of the group members velocities in the most external regions of the group. This could indicate the presence of galaxies in radial orbits in the external region of the group. We also found that the photometrical luminosity function is bimodal, showing a lack of M_{r} ~ -19.5 galaxies. The brightest group galaxy shows low Sersic parameter (n~2) and a small peculiar velocity. Indeed, our accurate photometry shows that the difference between the brightest and the second brightest galaxies is 1.9 mag in the r-band, while the classical definition of FGs is based on a magnitude gap of 2. We conclude that this fossil system does not follow the empirical definition of FGs. Nevertheless, it is a massive, old and undisturbed galaxy system with little infall of L^{*} galaxies since its initial collapse.

When a spinning system experiences a transient gravitational encounter with an external perturber, a quasi-resonance occurs if the spin frequency of the victim matches the peak orbital frequency of the perturber. Such encounters are responsible for the formation of long tails and bridges of stars during galaxy collisions. For high-speed encounters, the resulting velocity perturbations can be described within the impulse approximation. The traditional impulse approximation, however, does not distinguish between prograde and retrograde encounters, and therefore completely misses the resonant response. Here, using perturbation theory, we compute the effects of quasi-resonant phenomena on stars orbiting within a disk. Explicit expressions are derived for the velocity and energy change to the stars induced by tidal forces from an external gravitational perturber passing either on a straight line or parabolic orbit. Comparisons with numerical restricted three-body calculations illustrate the applicability of our analysis.

We present deep XMM observations and ESO WFI optical imaging of two X-ray-selected fossil group candidates, RXCJ0216.7-4749 and RXCJ2315.7-0222. Using the X-ray data, we derive total mass profiles under the hydrostatic equilibrium assumption. The central regions of RXCJ0216.7-4749 are found to be dominated by an X-ray bright AGN, and although we derive a mass profile, uncertainties are large and the constraints are significantly weakened due to the presence of the central source. The total mass profile of RXCJ2315.7-0222 is of high quality, being measured in fifteen bins from [0.075 - 0.75]R500 and containing three data points interior to 30 kpc, allowing comprehensive investigation of its properties. We probe several mass models based on the standard NFW profile or on the Sersic-like model recently suggested by high-resolution N-body simulations. We find that the addition of a stellar component due to the presence of the central galaxy is necessary for a good analytical model fit. In all mass profile models fitted, the mass concentration is not especially high compared to non-fossil systems. In addition, the modification of the dark matter halo by adiabatic contraction slightly improves the fit. However, our result depends critically on the choice of IMF used to convert galaxy luminosity to mass, which leads to a degeneracy between the central slope of the dark matter profile and the normalisation of the stellar component. While we argue on the basis of the range of M_*/L_R ratios that lower M_*/L_R ratios are preferred on physical grounds and that adiabatic contraction has thus operated in this system, better theoretical and observational convergence on this problem is needed to make further progess.

We employ numerical simulations and simple analytical estimates to argue that dark matter substructures orbiting in the inner regions of the Galaxy can be efficiently destroyed by disk shocking, a dynamical process known to affect globular star clusters. We carry out a set of fiducial high-resolution collisionless simulations in which we adiabatically grow a disk, allowing us to examine the impact of the disk on the substructure abundance. We also track the orbits of dark matter satellites in the high-resolution Aquarius simulations and analytically estimate the cumulative halo and disk shocking effect. Our calculations indicate that the presence of a disk with only 10% of the total Milky Way mass can significantly alter the mass function of substructures in the inner parts of halos. This has important implications especially for the relatively small number of satellites seen within ~30 kpc of the Milky Way center, where disk shocking is expected to reduce the substructure abundance by a factor of ~2 at 10^9 M$_{\odot}$ and ~3 at 10^7 M$_{\odot}$. The most massive subhalos with 10^10 M$_{\odot}$ survive even in the presence of the disk. This suggests that there is no inner missing satellite problem, and calls into question whether these substructures can produce transient features in disks, like multi-armed spiral patterns. Also, the depletion of dark matter substructures through shocking on the baryonic structures of the disk and central bulge may aggravate the problem to fully account for the observed flux anomalies in gravitational lens systems, and significantly reduces the dark matter annihilation signal expected from nearby substructures in the inner halo.

Dwarf spheroidal galaxies are the most dark matter dominated systems in the nearby Universe and their origin is one of the outstanding puzzles of how galaxies form. Dwarf spheroidals are poor in gas and stars, making them unusually faint, and those known as ultra-faint dwarfs have by far the lowest measured stellar content of any galaxy. Previous theories require that dwarf spheroidals orbit near giant galaxies like the Milky Way, but some dwarfs have been observed in the outskirts of the Local Group. Here we report simulations of encounters between dwarf disk galaxies and somewhat larger objects. We find that the encounters excite a process, which we term ``resonant stripping'', that can transform them into dwarf spheroidals. This effect is distinct from other mechanisms proposed to form dwarf spheroidals, including mergers, galaxy-galaxy harassment, or tidal and ram pressure stripping, because it is driven by gravitational resonances. It may account for the observed properties of dwarf spheroidals in the Local Group, including their morphologies and kinematics. Resonant stripping predicts that dwarf spheroidals should form through encounters, leaving detectable long stellar streams and tails.

The Magellanic Clouds were the largest members of a group of dwarf galaxies that entered the Milky Way (MW) halo at late times. This group, dominated by the LMC, contained ~4% of the mass of the Milky Way prior to its accretion and tidal disruption, but ~70% of the known dwarfs orbiting the MW. Our theory addresses many outstanding problems in galaxy formation associated with dwarf galaxies. First, it can explain the planar orbital configuration populated by some dSphs in the MW. Second, it provides a mechanism for lighting up a subset of dwarf galaxies to reproduce the cumulative circular velocity distribution of the satellites in the MW. Finally, our model predicts that most dwarfs will be found in association with other dwarfs. The recent discovery of Leo V (Belokurov et al. 2008), a dwarf spheroidal companion of Leo IV, and the nearby dwarf associations supports our hypothesis.

We consider the possibility that the Magellanic Clouds were the largest members of a group of dwarf galaxies that entered the Milky Way (MW) halo at late times. Seven of the eleven brightest satellites of the MW may have been part of this system. The proximity of some dwarfs to the plane of the orbit of the Large Magellanic Cloud (LMC) has been used to argue that they formed from tidal debris from the LMC and Small Magellanic Cloud (SMC). Instead, they may owe to the tidal breakup of the Magellanic Group. This can explain the association of many of the dwarf galaxies in the Local Group with the LMC system. It provides a mechanism for lighting up dwarf galaxies and reproduces the bright end of the cumulative circular velocity distribution of the satellites in the MW without invoking a stripping scenario for the subhalos to match the satellite distribution expected according to CDM theory. Finally, our model predicts that other isolated dwarfs will be found to have companions. Evidence for this prediction is provided by nearby, recently discovered dwarf associations.

We compare the observed merger rate of galaxies over cosmic time and the frequency of collisional ring galaxies (CRGs), with analytic models and halo merger and collision rates from a large cosmological simulation. In the Lambda cold dark matter (LCDM) model we find that the cosmic {\it merger fraction} does not evolve strongly between 0.2<z<2, implying that the observed decrease of the cosmic star formation rate since z~1 might not be tied to a disappearing population of major mergers. Halos hosting massive galaxies undergo on average ~2 mergers from z~2 up to present day, reflecting the late assembly time for the massive systems and the related downsizing problem. The cosmic {\it merger rate} declines with redshift: at the present time it is a factor of 10 lower than at z~2, in reasonable agreement with the current available data. The rate of CRG formation derived from the interactions between halo progenitors up to z=2 is found to be a good tracer of the cosmic merger rate. In the LCDM model the rate of CRGs as well as the merger rate do not scale as (1+z)^m, as suggested by previous models. Our predictions of cosmic merger and CRG rates may be applied to forthcoming surveys such as GOODS and zCOSMOS.

We use a numerical simulation of a loose group containing a Milky Way halo to probe that in the hierarchical universe the Magellanic Clouds and some dSphs have been accreted into the Milky Way halo from a late infalling group of dwarfs. Our simulations show that the tidal breakup of the Magellanic group occurs before it enters the Milky Way halo. Only half of the satellites contributed from the group are predicted to be inside the Milky Way virial radius. Half of its subhalos survive outside the current virial radius in the form of satellites, whereas the remaining material contributes to the diffuse Milky Way halo. At z~0 the disrupted group contributes less than 10% to the Milky Way halo mass but 20% of the brightest dwarf galaxies of the Milky Way have been part of this group. This scenario points out that some dSphs might have been form away from giant spirals and been accreted already as spheroids, by a late infall group in contrast with the classical picture of tidal stripping of dSph formation models. This would naturally explain several peculiarities of the local dSph: why Draco and the other luminous dSphs exist compared to other ultra-faint satellite galaxies, the location of Tucana and Cetus in the outskirts of the Local Group and the mismatch in metallicity between the stellar halo of the Milky Way and the dwarf galaxies that many have suspected dissolved to build it.

We investigate the origin and evolution of fossil groups in a concordance LCDM cosmological simulation. We consider haloes with masses between $(1-5)\times10^{13} \hMsun$ and study the physical mechanisms that lead to the formation of the large gap in magnitude between the brightest and the second most bright group member, which is typical for these fossil systems. Fossil groups are found to have high dark matter concentrations, which we can relate to their early formation time. The large magnitude-gaps arise after the groups have build up half of their final mass, due to merging of massive group members. We show that the existence of fossil systems is primarily driven by the relatively early infall of massive satellites, and that we do not find a strong environmental dependence for these systems. In addition, we find tentative evidence for fossil group satellites falling in on orbits with typically lower angular momentum, which might lead to a more efficient merger onto the host. We find a population of groups at higher redshifts that go through a ``fossil phase'': a stage where they show a large magnitude-gap, which is terminated by renewed infall from their environment.

We use a large cosmological N-body simulation to study the origin of possible correlations between the merging history and spin of cold dark matter halos. In particular, we examine claims that remnants of major mergers tend to have higher-than-average spins, and find that the effect is driven largely by unrelaxed systems: equilibrium dark matter halos show no significant correlation between spin and merger history. Out-of-equilibrium halos have, on average, higher spins than relaxed systems, suggesting that the virialization process leads to a net decrease in the value of the spin parameter. We find that this decrease is driven by the internal redistribution of mass and angular momentum that occurs during virialization, a process that is especially efficient during major mergers, when high angular momentum material is pushed beyond the virial radius of the remnant. Since such redistribution likely affects the angular momentum of baryons and dark matter unevenly, our findings question the common practice of identifying the specific angular momentum content of a halo with that of its embedded luminous component. Further work is needed to elucidate the true relation between the angular momentum content of baryons and dark matter in galaxy systems assembled hierarchically.

We use large cosmological N-body simulations to study the subhalo population in galaxy group sized halos. In particular, we look for fossil group candidates with typical masses ~10-25% of Virgo cluster but with an order of magnitude less substructure. We examine recent claims that the earliest systems to form are deficient enough in substructure to explain the luminosity function found in fossil groups. Although our simulations show a correlation between the halo formation time and the number of subhalos, the maximum suppression of subhalos is a factor of 2-2.5, whereas a factor of 6 is required to match fossil groups and galaxies. While the number of subhalos depends weakly on the formation time, the slope of the halo substructure velocity function does not. The satellite population within Cold Dark Matter (CDM) halos is self-similar at scales between galaxies and galaxy clusters regardless of mass, whereas current observations show a break in self-similarity at a mass scale corresponding to group of galaxies.

The processes are investigated by which gas loses its angular momentum during the protogalactic collapse phase, leading to disk galaxies that are too compact with respect to the observations. High-resolution N-body/SPH simulations in a cosmological context are presented including cold gas and dark matter. A halo with quiet merging activity since z~3.8 and with a high spin parameter is analysed that should be an ideal candidate for the formation of an extended galactic disk. We show that the gas and the dark matter have similar specific angular momenta until a merger event occurs at z~2 with a mass ratio of 5:1. All the gas involved in the merger loses a substantial fraction of its specific angular momentum due to tidal torques and falls quickly into the center. Dynamical friction plays a minor role,in contrast to previous claims. In fact, after this event a new extended disk begins to form from gas that was not involved in the 5:1 merger event and that falls in subsequently. We argue that the angular momentum problem of disk galaxy formation is a merger problem: in cold dark matter cosmology substantial mergers with mass ratios of 1:1 to 6:1 are expected to occur in almost all galaxies. We suggest that energetic feedback processes could in principle solve this problem, however only if the heating occurs at the time or shortly before the last substantial merger event. Good candidates for such a coordinated feedback would be a merger-triggered star burst or central black hole heating. If a large fraction of the low angular momentum gas would be ejected as a result of these processes, late-type galaxies could form with a dominant extended disk component, resulting from late infall, a small bulge-to-disk ratio and a low baryon fraction, in agreement with observations.

We use a set of twelve high-resolution N-body/hydrodynamical simulations in the $\Lambda$CDM cosmology to investigate the origin and formation rate of fossil groups (FGs), which are X-ray bright galaxy groups dominated by a large elliptical galaxy, with the second brightest galaxy being at least two magnitudes fainter. The simulations invoke star formation, chemical evolution with non-instantaneous recycling, metal dependent radiative cooling, strong star burst driven galactic super winds, effects of a meta-galactic UV field and full stellar population synthesis. We find an interesting correlation between the magnitude gap between the first and second brightest galaxy and the formation time of the group. It is found that FGs have assembled half of their final dark matter mass already at $z\ga1$, and subsequently typically grow by minor merging only, wheras non-FGs on average form later. The early assembly of FGs leaves sufficient time for galaxies of $L \sim L_*$ to merge into the central one by dynamical friction, resulting in the large magnitude gap at $z=0$. A fraction of 33$\pm$16% of the groups simulated are found to be fossil, whereas the observational estimate is $\sim$10-20%. The FGs are found to be X-ray over-luminous relative to non-FGs of the same optical luminosity, in qualitative agreement with observations. Finally, from a dynamical friction analysis is found that only because infall of $L \sim L_*$ galaxies happens along filaments with small impact parameters do FGs exist at all.

The importance of angular momentum in regulating the sizes of galactic disks and by this their star formation history is highlighted. Tidal torques and accretion of satellites in principle provide enough angular momentum to form disks with sizes that are in agreement with observations. However three major problems have been identified that challenge cold dark matter theory and affect models of galaxy evolution: (1) too much angular momentum is transferred from the gas to the dark halos during infall, leading to disks with scale lengths that are too small, (2) bulgeless disks require more specific angular momentum than is generated cosmologically even if gas would not lose angular momentum during infall, (3) gravitational torques and hierarchical merging produce a specific angular momentum distribution that does not match the distribution required to form exponential disks naturally; some gas has exceptionally high angular momentum, leading to extended outer disks while another large gas fraction will contain very little specific angular momentum and is expected to fall into the galactic center, forming a massive and dominant bulge component. Any self-consistent theory of galaxy formation will require to provide solutions to these questions. Selective mass loss of low-angular-momentum gas in an early phase of galaxy evolution currently seems to be the most promising scenario. Such a process would have a strong affect on the early protogalactic evolution phase, the origin and evolution of galactic morphologies and link central properties of galaxies like the origin of central massive black holes with their global structure.

The specific angular momentum of Cold Dark Matter (CDM) halos in a $\Lambda$CDM universe is investigated. Their dimensionless specific angular momentum $\lambda'=\frac{j}{\sqrt{2}V_{vir} R{vir}}$ with $V_{vir}$ and $R_{vir}$ the virial velocity and virial radius, respectively depends strongly on their merging histories. We investigate a set of $\Lambda$CDM simulations and explore the specific angular momentum content of halos formed through various merging histories. Halos with a quiet merging history, dominated by minor mergers and accretion until the present epoch, acquire by tidal torques on average only 2% to 3% of the angular momentum required for their rotational support ($\lambda'=0.02$). This is in conflict with observational data for a sample of late-type bulgeless galaxies which indicates that those galaxies reside in dark halos with exceptionally high values of $\lambda' \approx 0.06-0.07$. Minor mergers and accretion preserve or slowly increase the specific angular momentum of dark halos with time. This mechanism is however not efficient enough in order to explain the observed spin values for late-type dwarf galaxies. Energetic feedback processes have been invoked to solve the problem that gas loses a large fraction of its specific angular momentum during infall. Under the assumption that dark halos hosting bulgeless galaxies acquire their mass via quiescent accretion, our results indicate yet another serious problem: the specific angular momentum gained during the formation of these objects is not large enough to explain their observed rotational properties,even if no angular momentum would be lost during gas infall.

The Cold Dark Matter (CDM) theory predicts a wealth of substructure within dark halos. These predictions match observations of galaxy clusters like the nearby Virgo cluster. However, CDM has a "small scale crisis" since galaxies dominate the halo with little substructure while the model predicts that galaxies should be scaled versions of galaxy clusters with abundant substructure. Compared to CDM predictions, the Milky Way and Andromeda are "missing" objects with velocity dispersions $\sigma \geq 10$ km s$^{-1}$. The energy scale of these missing satellites is low enough that stellar winds and supernovae might remove gas and suppress the formation of their luminous stellar components. Here, we show that the small scale crisis persists in fossil groups that have masses of up to 40% of the nearby Virgo cluster of galaxies. Fossil groups are missing satellites with luminosities that occur at the predicted frequency in the Virgo cluster. Moreover, the "missing galaxies" in fossil groups are nearly as luminous as the Milky Way with a velocity dispersion $\sigma \leq$ 150 km s$^{-1}$.

Self-interacting dark matter was proposed by Spergel & Steinhardt (2000) to alleviate two conflicts between Cold Dark Matter (CDM) models and observations. Firstly, CDM N-body simulations predict dark matter halo density profiles that diverge at the centre in disagreement with the constant density cores observed in late-type dwarf and Low Surface Brightness (LSB) galaxies. Secondly, N-body simulations predict an overabundance of subhalos in the Galactic halo. Using a simple semi-analytical argument we show that weakly self-interacting dark matter models, which can produce halo cores of the sizes observed in dark matter dominated galaxies, are unable to reconcile the number of satellites in the Galactic halo with the observed number of dwarf galaxies in the Local Group.

We propose a new numerical procedure to simulate a single dark halo of any size and mass in a hierarchical framework coupling the extended Press-Schechter formalism (EPSF) to N-body simulations. The procedure consists of assigning cosmological initial conditions to the particles of a single halo with a EPSF technique and following only the dynamical evolution using a serial N-body code. The computational box is fixed with a side of $0.5 h^{-1}$ Mpc. This allows to simulate galaxy cluster halos using appropriate scaling relations, to ensure savings in computing time and code speed. The code can describe the properties of halos composed of collisionless or collisional dark matter. For collisionless Cold Dark Matter (CDM) particles the NFW profile is reproduced for galactic halos as well as galaxy cluster halos. Using this numerical technique we study some characteristics of halos assumed to be isolated or placed in a cosmological context in presence of weak self-interacting dark matter: the soft core formation and the core collapse. The self-interacting dark matter cross section per unit mass is assumed to be inversely proportional to the particle collision velocity: $\sigma/m_{x} \propto 1/v$.

We present high resolution Halpha rotation curves of 4 late-type dwarf galaxies and 2 low surface brightness galaxies (LSB) for which accurate HI rotation curves are available from the literature. Observations are carried out at Telescopio Nazionale Galileo (TNG). For LSB F583-1 an innovative dispersing element was used, the Volume Phase Holographic (VPH) with a dispersion of about 0.35 A/pxl. We find good agreement between the Halpha data and the HI observations and conclude that the HI data for these galaxies suffer very little from beam smearing. We show that the optical rotation curves of these dark matter dominated galaxies are best fitted by the Burkert profile. In the centers of galaxies, where the N-body simulations predict cuspy cores and fast rising rotation curves, our data seem to be in better agreement with the presence of soft cores.

This work is focused on the preliminary results of the observations of H\alpha rotation curves for some of the late type dwarf and LSB galaxies carried out at the TNG telescope. In light of the observational data and of the N-body simulations we have performed recently, we discuss some of the implications on the nature of the DM particles and the formation of the dark haloes.

We present high spatial resolution H\alpha rotation curves of late-type dwarf and LSB galaxies. From our analysis we find good agreement between our H\alpha data and the HI observations taken from the literature, concluding that the HI rotation curves for these galaxies suffer very little from beam smearing. A preliminary analysis of our data rules out the CDM model in the inner regions of these galaxies.

(Abridged) Using high-resolution cosmological N-body simulations, we study how the density profiles of dark matter halos are affected by the filtering of the density power spectrum below a given scale length and by the introduction of a thermal velocity dispersion. In the warm dark matter (WDM) scenario, both the free-streaming scale, R_f, and the velocity dispersion, v_w, are determined by the mass m_w of the WDM particle. We find that v_w is too small to affect the density profiles of WDM halos. Down to the resolution attained in our simulations, there is not any significant difference in the density profiles and concentrations of halos obtained in simulations with and without the inclusion of v_w. The density profiles of halos with masses down to ~0.01 the filtering mass M_f can be described by the NFW shape; significant soft cores are not formed. Nevertheless, the concentrations of these halos are lower than those of the CDM counterparts and are approximately independent of mass. The lower concentrations of WDM halos with respect to their CDM counterparts can be accounted for their late formation epoch. Overall, our results point to a series of advantages of a WDM model over the CDM one. In addition to solving the substructure problem, a WDM model with R_f~0.16 Mpc (m_w~0.75 kev; flat cosmology with Omega_L=h=0.7) also predicts concentrations, a Tully-Fisher relation, and formation epochs for small halos which seems to be in better agreement with observations, relative to CDM predictions.

We investigate the density profiles and evolution of weakly self-interacting cold dark matter haloes using a numerical code based on the collisional Boltzmann equation. This approach is alternative to N-body techniques in following the dynamical evolution of haloes in the cosmological context and taking into account particle self-interaction. The physical case with a cross section inversely proportional to the relative velocity of the colliding particles is modelled with an unprecedented resolution, spanning five orders of magnitude on the radius for each halo. The modelled haloes cover a mass range from dwarf galaxies to galaxy clusters. We find that for \sigma v_{100} \approx 10^{-24} cm^2/GeV, where \sigma is the cross section per unit mass and v_{100} is the collision velocity in units of 100 km/s, soft cores in good agreement with observations on galactic as well as on galaxy cluster scales are obtained. Remarkably, the observed nearly invariance of the halo central density with mass is reproduced.

One of the predictions of the standard CDM is that dark haloes have centrally divergent density profiles. An extensive body of rotation curve observations of dwarf and low surface brightness galaxies shows the dark haloes of those systems to be characterized by soft constant density central cores. Several physical processes have been proposed to produce soft cores in dark haloes, each one with different scaling properties. With the aim of discriminating among them we have examined the rotation curves of dark matter dominated dwarf and low surface brightness galaxies and the inner mass profiles of two clusters of galaxies lacking a central cD galaxy and with evidence of soft cores in the centre. The core radii and central densities of these haloes scale in a well defined manner with the depth of their potential wells, as measured through the maximum circular velocity. As a result of our analysis we identify self-interacting CDM as a viable solution to the core problem, where a non-singular isothermal core is formed in the halo center surrounded by a Navarro, Frenk, & White profile in the outer parts. We show that this particular physical situation predicts core radii in agreement with observations. Furthermore, using the observed scalings, we derive an expression for the minimum cross section (\sigma) which has an explicit dependence with the halo dispersion velocity (v). If m_x is the mass of the dark matter particle: \sigma/m_x ~4 10^-25 (v/100 km s^-1)^-1 cm^2/Gev.

Using observational data for a sample of dark matter dominated galaxies and two cluster of galaxies, we have found that the central halo density does not depend on its mass, and the core radius is roughly proportional to the maximum rotation velocity. A good agreement with these scaling laws is obtained for CDM halos whose dense inner parts were expanded by gravothermal instabilities if the particles efficiently self-interact only in these parts. We find that the particle cross-section is inversely proportional to the velocity dispersion.

Within the framework of the cold dark matter (CDM) cosmogony, a central cusp in the density profiles of virialized dark haloes is predicted. This prediction disagrees with the soft inner halo mass distribution inferred from observations of dwarf and low surface brightness galaxies, and some clusters of galaxies. By analysing data for some of these objects, we find that the halo central density is nearly independent of the mass from galactic to galaxy cluster scales with an average value of around 0.02 M_sun/pc^3. We show that soft cores can be produced in the CDM haloes by introducing a lower cut-off in the power spectra of fluctuations and assuming high orbital thermal energies during halo formation. However, the scale invariance of the halo central density is not reproduced in these cases. The introduction of self-interaction in the CDM particles offers the most attractive alternative to the core problem. We propose gravothermal expansion as a possible mechanism to produce soft cores in the CDM haloes with self-interacting particles. A global thermodynamical equilibrium can explain the central density scale invariance. We find a minimum cross section capable of establishing isothermal cores in agreement with the observed shallow cores. If \sigma and m_x are the cross section and mass of the dark matter particle, and v is the halo velocity dispersion, then (\sigma/m_x) \~ 4 10^{-25} (100 km s^{-1}/v) cm^2/GeV.