We present numerical relativity simulations of cosmological scenarios in which the universe is smoothed and flattened by undergoing a phase of slow contraction and test their sensitivity to a wide range of initial conditions. Our numerical scheme enables the variation of all freely specifiable physical quantities that characterize the initial spatial hypersurface, such as the initial shear and spatial curvature contributions as well as the initial field and velocity distributions of the scalar that drives the cosmological evolution. In particular, we include initial conditions that are far outside the perturbative regime of the well-known attractor scaling solution. We complement our numerical results by analytically performing a complete dynamical systems analysis and show that the two approaches yield consistent results.

We explore the effect of a supermassive black hole (SMBH) on the density profile of a fuzzy dark matter (FDM) soliton core at the centre of a dark matter halo. We numerically solve the Schrödinger-Poisson equations, treating the black hole as a gravitational point mass, and demonstrate that this additional perturbing term has a 'squeezing' effect on the soliton density profile, decreasing the core radius and increasing the central density. In the limit of large black hole mass, the solution approaches one akin to the hydrogen atom, with radius inversely proportional to the black hole mass. By applying our analysis to two specific galaxies (M87 and the Milky Way) and pairing it with known observational limits on the amount of centrally concentrated dark matter, we obtain a constraint on the FDM particle mass, finding that the range 10 −22.2 eV m 10 −21.7 eV should be forbidden (taking into account additional factors concerning the life-time of the soliton in the vicinity of a black hole). ; Einstein Fellow matter (FDM), this scalar field is assumed to have a particle mass of ∼ 10 −22 eV and a de Broglie wavelength of λ = 1.2 m 10 −22 eV 100 km s −1 v

Relying on the dramatic increase in the number of stars with full 6D phase-space information provided by the Gaia Data Release 3, we resolve the distribution of the stellar halo around the Sun to uncover signatures of incomplete phase-mixing. We show that for the stars likely belonging to the last massive merger, the (vr, r) distribution contains a series of long and thin chevron-like overdensities. These phase-space sub-structures have been predicted to emerge following the dissolution of a satellite, when its tidal debris is given time to wind up, thin out and fold. Such chevrons have been spotted in external galaxies before, here we report the first detection in our own Milky Way. We also show that the observed angular momentum Lz distribution appears more prograde at high energies, possibly revealing the original orbital angular momentum of the in-falling galaxy. The energy distribution of the debris is strongly asymmetric with a peak at low E – which, we surmise, may be evidence of the dwarf’s rapid sinking – and riddled with wrinkles and bumps. We demonstrate that similar phase-space and (E, Lz) sub-structures are present in numerical simulations of galaxy interactions, both in bespoke N-body runs and in cosmological hydrodynamical zoom-in suites. The remnant traces of the progenitor’s disruption and the signatures of the on-going phase-mixing discovered here will not only help to constrain the properties of our Galaxy’s most important interaction, but also can be used as a novel tool to map out the Milky Way’s current gravitational potential and its perturbations.

In a galaxy merger, the stars tidally stripped from the satellite and accreted onto the host galaxy undergo phase mixing and form finely-grained structures in the phase space. However, these fragile structures may be destroyed in the subsequent galaxy evolution, in particular, by a rotating bar that appears well after the merger is completed. In this work, we investigate the survivability of phase-space structures in the presence of a bar. We find that a bar with amplitude and pattern speed similar to those of the Milky Way would blur and destroy a substantial amount of the substructure that consists of particles with pericentre radii comparable to the bar length. While this appears to be in tension with the recent discovery of phase-space chevrons in Gaia DR3 data, the most prominent chevrons in our simulations can still be recovered when applying the same analysis procedure as in observations. Moreover, the smoothing effect is less pronounced in the population of stars whose angular momenta have the opposite sign to the bar pattern speed.

Using data from the Gaia satellite’s Radial Velocity Spectrometer Data Release 3 (RVS, DR3), we find a new and robust feature in the phase space distribution of halo stars. It is a prominent ridge at constant energy and with angular momentum Lz > 0. We run test particle simulations of a stellar halo-like distribution of particles in a realistic Milky Way potential with a rotating bar. We observe similar structures generated in the simulations from the trapping of particles in resonances with the bar, particularly at the corotation resonance. Many of the orbits trapped at the resonances are halo-like, with large vertical excursions from the disc. The location of the observed structure in energy space is consistent with a bar pattern speed in the range Ωb ≈ 35 − 40 km s−1 kpc−1. Overall, the effect of the resonances is to give the inner stellar halo a mild, net spin in the direction of the bar’s rotation. As the distribution of the angular momentum becomes asymmetric, a population of stars with positive mean Lz and low vertical action is created. The variation of the average rotational velocity of the simulated stellar halo with radius is similar to the behaviour of metal-poor stars in data from the APOGEE survey. Though the effects of bar resonances have long been known in the Galactic disc, this is strong evidence that the bar can drive changes even in the diffuse and extended stellar halo through its resonances.