The power spectrum response function of the large-scale structure of the Universe describes how the evolved power spectrum is modified by a small change in initial power through non-linear mode coupling of gravitational evolution. It was previously found that the response function for the coupling from small to large scales is strongly suppressed in amplitude, especially at late times, compared to predictions from perturbation theory (PT) based on the single-stream approximation. One obvious explanation for this is that PT fails to describe the dynamics beyond shell-crossing. We test this idea by comparing measurements in $N$-body simulations to prescriptions based on PT but augmented with adaptive smoothing to account for the formation of non-linear structures of various sizes in the multi-stream regime. We first start with one-dimensional (1D) cosmology, where the Zel'dovich approximation provides the exact solution in the single stream regime. Similarly to the 3D case, the response function of the large-scale modes exhibits a strong suppression in amplitude at small scales which cannot be explained by the Zel'dovich solution alone. However, by performing adaptive smoothing of initial conditions to identify haloes of different sizes and solving approximately post-collapse dynamics in the 3-streams regime, agreement between theory and simulations drastically improves. We extend our analyses to the 3D case using PINOCCHIO algorithm, in which similar adaptive smoothing is implemented on the Lagrangian PT fields to identify haloes and is combined with a spherical halo prescription to account for post-collapse dynamics. Again, a suppression is found in the coupling between small- and large-scale modes and the agreement with simulations is improved.

Fuzzy dark matter of ultra-light axions has gained attention, largely in light of the galactic scale problems associated with cold dark matter. But the large de Broglie wavelength, believed to possibly alleviate these problems, also leads to fluctuations that place constraints on ultra-light axions. We adapt and extend a method, previously devised to describe the effect of gaseous fluctuations on cold dark matter cusps, in order to determine the imprints of ultra-light axion haloes on the motion of classical test particles. We first evaluate the effect of fluctuations in a statistically homogeneous medium of classical particles, then in a similar system of ultra light axions. In the first case, one recovers the classical two body relaxation time (and diffusion coefficients) from white noise density fluctuations. In the second situation, the fluctuations are not born of discreteness noise but from the finite de Broglie wavelength; correlation therefore exists over this scale, while white noise is retained on larger scales, elucidating the correspondence with classical relaxation. The resulting density power spectra and correlation functions are compared with those inferred from numerical simulations, and the relaxation time arising from the associated potential fluctuations is evaluated. We then apply our results to estimate the heating of disks embedded in axion dark haloes. We find that this implies an axion mass $m \ga 2 \times 10^{-22} {\rm eV}$. We finally apply our model to the case of the central cluster of Eridanus II, confirming that far stronger constraints on $m$ may in principle be obtained, and discussing the limitations associated with the assumptions leading to these.

In the mean field limit, isolated gravitational systems often evolve towards a steady state through a violent relaxation phase. One question is to understand the nature of this relaxation phase, in particular the role of radial instabilities in the establishment/destruction of the steady profile. Here, through a detailed phase-space analysis based both on a spherical Vlasov solver, a shell code and a $N$-body code, we revisit the evolution of collisionless self-gravitating spherical systems with initial power-law density profiles $\rho(r) \propto r^n$, $0 \leq n \leq -1.5$, and Gaussian velocity dispersion. Two sub-classes of models are considered, with initial virial ratios $\eta=0.5$ ("warm") and $\eta=0.1$ ("cool"). Thanks to the numerical techniques used and the high resolution of the simulations, our numerical analyses are able, for the first time, to show the clear separation between two or three well known dynamical phases: (i) the establishment of a spherical quasi-steady state through a violent relaxation phase during which the phase-space density displays a smooth spiral structure presenting a morphology consistent with predictions from self-similar dynamics, (ii) a quasi-steady state phase during which radial instabilities can take place at small scales and destroy the spiral structure but do not change quantitatively the properties of the phase-space distribution at the coarse grained level and (iii) relaxation to non spherical state due to radial orbit instabilities for $n \leq -1$ in the cool case.

We study how migration affects stars of a galaxy with a thin stellar disc and thicker stellar components. The simulated galaxy has a strong bar and lasting spiral arms. We find that the amplitude of the churning (change in angular momentum) is similar for thin and thick components, and of limited amplitude, and that stars of all components can be trapped at the corotation of the bar. At the exception of those stars trapped at the corotation, we find that stars that are far from their initial guiding radius are more likely so due to blurring rather than churning effects. We compare the simulation to orbits integration with a fixed gravitational potential rotating at a constant speed. In the latter case, stars trapped at corotation are churned periodically outside and inside the corotation radius, with a zero net average. However, as the bar speed of the simulated galaxy decreases and its corotation radius increases, stars trapped at corotation for several Gyrs can be churned outwards on average. We study the location of extreme migrators (stars experimenting the largest churning) and find that extreme migrators come from regions on the leading side of the effective potential local maxima.

Radio interferometers have the ability to precisely localize and better characterize the properties of sources. This ability is having a powerful impact on the study of fast radio transients, where a few milliseconds of data is enough to pinpoint a source at cosmological distances. However, recording interferometric data at millisecond cadence produces a terabyte-per-hour data stream that strains networks, computing systems, and archives. This challenge mirrors that of other domains of science, where the science scope is limited by the computational architecture as much as the physical processes at play. Here, we present a solution to this problem in the context of radio transients: realfast, a commensal, fast transient search system at the Jansky Very Large Array. Realfast uses a novel architecture to distribute fast-sampled interferometric data to a 32-node, 64-GPU cluster for real-time imaging and transient detection. By detecting transients in situ, we can trigger the recording of data for those rare, brief instants when the event occurs and reduce the recorded data volume by a factor of 1000. This makes it possible to commensally search a data stream that would otherwise be impossible to record. This system will search for millisecond transients in more than 1000 hours of data per year, potentially localizing several Fast Radio Bursts, pulsars, and other sources of impulsive radio emission. We describe the science scope for realfast, the system design, expected outcomes, and ways real-time analysis can help in other fields of astrophysics.

In some galaxies, the stars orbiting the supermassive black hole take the form of an eccentric nuclear disk, in which every star is on a coherent, apsidally-aligned orbit. The most famous example of an eccentric nuclear disk is the double nucleus of Andromeda, and there is strong evidence for many more in the local universe. Despite their apparent ubiquity however, a dynamical explanation for their longevity has remained a mystery: differential precession should wipe out large-scale apsidal-alignment on a short timescale. Here we identify a new dynamical mechanism which stabilizes eccentric nuclear disks, and explain for first time the negative eccentricity gradient seen in the Andromeda nucleus. The stabilizing mechanism drives oscillations of the eccentricity vectors of individual orbits, both in direction (about the mean body of the disk) and in magnitude. Combined with the negative eccentricity gradient, the eccentricity oscillations push some stars near the inner edge of the disk extremely close to the black hole, potentially leading to tidal disruption events. Order of magnitude calculations predict extremely high rates in recently-formed eccentric nuclear disks ($\sim0.1 - 1$ ${\rm yr}^{-1} {\rm gal}^{-1}$). Unless the stellar disks are replenished, these rates should decrease with time as the disk depletes in mass. If eccentric nuclear disks form during gas-rich major mergers, this may explain the preferential occurrence of tidal disruption events in recently-merged and post-merger (E+A/K+A) galaxies.

Radial stellar migration in galactic discs has received much attention in studies of galactic dynamics and chemical evolution, but remains a dynamical phenomenon that needs to be fully quantified. In this work, using a Tree-SPH simulation of an Sb-type disc galaxy, we quantify the effects of blurring (epicyclic excursions) and churning (change of guiding radius). We quantify migration (either blurring or churning) both in terms of flux (the number of migrators passing at a given radius), and by estimating the population of migrators at a given radius at the end of the simulation compared to non-migrators, but also by giving the distance over which the migration is effective at all radii. We confirm that the corotation of the bar is the main source of migrators by churning in a bar-dominated galaxy, its intensity being directly linked to the episode of a strong bar, in the first 1-3 Gyr of the simulation. We show that within the outer Lindblad resonance (OLR), migration is strongly dominated by churning, while blurring gains progressively more importance towards the outer disc and at later times. Most importantly, we show that the OLR limits the exchange of angular momentum, separating the disc in two distinct parts with minimal or null exchange, except in the transition zone, which is delimited by the position of the OLR at the epoch of the formation of the bar, and at the final epoch. We discuss the consequences of these findings for our understanding of the structure of the Milky Way disc. Because the Sun is situated slightly outside the OLR, we suggest that the solar vicinity may have experienced very limited churning from the inner disc.

By means of idealized, dissipationless N-body simulations which follow the formation and subsequent buckling of a stellar bar, we study the characteristics of boxy/peanut-shaped bulges and compare them with the properties of the stellar populations in the Milky Way bulge. The main results of our modeling, valid for the general family of boxy/peanut shaped bulges, are the following: (i) because of the spatial redistribution in the disk initiated at the epoch of bar formation, stars from the innermost regions to the outer Lindblad resonance of the stellar bar are mapped into a boxy bulge; (ii) the contribution of stars to the local bulge density depends on their birth radius: stars born in the innermost disk tend to dominate the innermost regions of the boxy bulge, while stars originating closer to the OLR are preferentially found in the outer regions of the boxy/peanut structure; (iii) stellar birth radii are imprinted in the bulge kinematics, the larger the birth radii of stars ending up in the bulge, the greater their rotational support and the higher their line-of- sight velocity dispersions (but note that this last trend depends on the bar viewing angle); (iv) the higher the classical bulge-over-disk ratio, the larger its fractional contribution of stars at large vertical distance from the galaxy mid-plane. (ABRIDGED) On the basis of their chemical and kinematic characteristics, the results of our modeling suggests that the populations A, B and C, as defined by the ARGOS survey, can be associated, respectively, with the inner thin disk, to the young thick and to the old thick disk, following the nomenclature recently suggested for stars in the solar neighborhood by Haywood et al. (2013).

Recent work in galaxy formation has enlightened the important role of baryon physics, to solve the main problems encountered by the standard theory at the galactic scale, such as the galaxy stellar mass functions, or the missing satellites problem. The present work aims at investigating in particular the role of the cold and dense molecular phase, which could play a role of gas reservoir in the outer galaxy discs, with low star formation efficiency. Through TreeSPH simulations, implementing the cooling to low temperatures, and the inclusion of the molecular hydrogen component, several feedback efficiencies are studied, and results on the gas morphology and star formation are obtained. It is shown that molecular hydrogen allows some slow star formation (with gas depletion times of about 5 Gyr) to occur in the outer parts of the discs. This dense and quiescent phase might be a way to store a significant fraction of dark baryons, in a relatively long time-scale, in the complete baryonic cycle, connecting the galaxy discs to hot gaseous haloes and to the cosmic filaments.

We propose a unified single-field description of the galactic Dark Matter and various uniform scalar fields for the inflation and cosmological constant. The two types of effects could originate from a fluid of both spatially and temporally varying Vacuum Energy if the vacuum has an uneven pressure caused by a photon-like vector field (of perhaps an unstable massive boson). We propose a most general Lagrangian with a {\bf N}on-{\bf u}niform Cosmological Constant for this vacuum fluid (dubbed as a Nu-Lambda fluid), working within the framework of Einsteinian gravity. This theory includes a continuous spectrum of plausible dark energy theories and gravity theories, e.g., inflation, quintessence, k-essence, f(R), Generalized Einstein-Aether f(K), MOND, TeVeS, BSTV etc. theories. It also suggests new models such as a certain f(K+R) model, which suggests intriguing corrections to MOND depending of redshift and density. Some specific constructions of the Nu-Lambda fluid (e.g., Zhao's V-$\Lambda$ model) closely resemble the $\Lambda$CDM cosmology on large scale, but fit galaxy rotation curves as good as MOND. Perturbed Einstein Equations in a simple $f(K_4)$ model are solvable and show effects of a DM coupled to DE. Incorporating the perturbation equations here into standard simulations for cosmological structure growth offers a chance to falsify examples of the Nu-Lambda theories.

The MOdified Newtonian Dynamics (MOND) is presented here, as well as a theory that can be linked to it: the theory of the Aether, a four-vector field breaking Lorentz invariance. The form of its Lagrangian is studied, then basic equations of the theory are rederived in a detailed way, and calculated for different metrics, exploring the impact of non-zero spatial terms of the Aether. A brief attempt of making the Aether Lagrangian depend on a scalar field is presented. An analytic solving of a galaxy model with an external field is described, which highlights the MONDian external field effect that breaks the strong equivalence principle.