Lucio Mayer

Forming disc galaxies in ΛCDM simulations

We used fully cosmological, high resolution N-body + SPH simulations to follow the formation of disk galaxies with rotational velocities between 135 and 270 km/sec in a ΛCDM universe. The simulations include gas cooling, star formation, the effects of a uniform UV background and a physically motivated description of feedback from supernovae. The host dark matter halos have a spin and last major merger redshift typical of galaxy sized halos as measured in recent large scale N–Body simulations. The simulated galaxies form rotationally supported disks with realistic exponential scale lengths and fall on both the I-band and baryonic Tully Fisher relations. An extended stellar disk forms inside the Milky Way sized halo immediately after the last major merger. The combination of UV background and SN feedback drastically reduces the number of visible satellites orbiting inside a Milky Way sized halo, bringing it in fair agreement with observations. Our simulations predict that the average age of a primary galaxy’s stellar population decreases with mass, because feedback delays star formation in less massive galaxies. Galaxies have stellar masses and current star formation rates as a function of total mass that are in good agreement with observational data. We discuss how both high mass and force resolution and a realistic description of star formation and feedback are important ingredients to match the observed properties of galaxies.

Fundamental differences between SPH and grid methods

We have carried out a comparison study of hydrodynamical codes by investigating their performance in modelling interacting multiphase fluids. The two commonly used techniques of grid and smoothed particle hydrodynamics (SPH) show striking differences in their ability to model processes that are fundamentally important across many areas of astrophysics. Whilst Eulerian grid based methods are able to resolve and treat important dynamical instabilities, such as Kelvin-Helmholtz or Rayleigh-Taylor, these processes are poorly or not at all resolved by existing SPH techniques. We show that the reason for this is that SPH, at least in its standard implementation, introduces spurious pressure forces on particles in regions where there are steep density gradients. This results in a boundary gap of the size of an SPH smoothing kernel radius over which interactions are severely damped.

The Metamorphosis of Tidally Stirred Dwarf Galaxies

We present results from high-resolution N-body/SPH (smoothed particle hydrodynamic) simulations of rotationally supported dwarf irregular galaxies moving on bound orbits in the massive dark matter halo of the Milky Way. The dwarf models span a range in disk surface density and the masses and sizes of their dark halos are consistent with the predictions of cold dark matter cosmogonies. We show that the strong tidal field of the Milky Way determines severe mass loss in their halos and disks and induces bar and bending instabilities that transform low surface brightness dwarfs (LSBs) into dwarf spheroidals (dSphs) and high surface brightness dwarfs (HSBs) into dwarf ellipticals (dEs) in less than 10 Gyr. The final central velocity dispersions of the remnants are in the range 8-30 km s-1 and their final v/? falls to values less than 0.5, matching well the kinematics of early-type dwarfs. The transformation requires the orbital time of the dwarf to be 3-4 Gyr, which implies a halo as massive and extended as predicted by hierarchical models of galaxy formation to explain the origin of even the farthest dSph satellites of the Milky Way, Leo I, and Leo II. We show that only dwarfs with central dark matter densities as high as those of Draco and Ursa Minor can survive for 10 Gyr in the proximity of the Milky Way. A correlation between the central density and the distance of the dwarfs from the primary galaxy is indeed expected in hierarchical models, in which the densest objects should have small orbital times because of their early formation epochs. Part of the gas is stripped and part is funneled to the center because of the bar, generating one strong burst of star formation in HSBs and smaller, multiple bursts in LSBs. Therefore, the large variety of star formation histories observed in Local Group dSphs arises because different types of dIrr progenitors respond differently to the external perturbation of the Milky Way. Our evolutionary model naturally explains the morphology-density relation observed in the Local Group and in other nearby loose groups. Extended low surface brightness stellar and gaseous streams originate from LSBs and follow the orbit of the dwarfs for several gigayears. Because of their high velocities, unbound stars projected along the line of sight can lead to overestimating the mass-to-light ratio of the bound remnant by a factor 2, but this does not eliminate the need of extremely high dark matter contents in some of the dSphs.

The Secular Evolution of Disk Structural Parameters

We present a comprehensive series of simulations to study the secular evolution of disk galaxies expected in a ΛCDM universe. Our simulations are organized in a hierarchy of increasing complexity, ranging from rigid-halo collisionless simulations to fully live simulations with gas and star formation. Our goal is to examine which structural properties of disk galaxies may result from secular evolution rather than from hierarchical assembly. In the vertical direction, we find that various mechanisms lead to heating, the strongest of which is the buckling instability of a bar, which leads to peanut-shaped bulges; these can be recognized face-on even in the presence of gas. We find that bars are robust structures that survive buckling and require a large (~20% of the total mass of the disk) central mass concentration to be destroyed. This can occur in dissipative simulations, where bars induce strong gas inflows, but requires that radiative cooling overcome heating. We show how angular momentum redistribution leads to increasing central densities and disk scale lengths and to profile breaks at large radii. The breaks in these simulations are in excellent agreement with observations, even when the evolution is collisionless. Disk scale lengths increase even when the total disk angular momentum is conserved; thus, mapping halo angular momenta to scale lengths is nontrivial. A decomposition of the resulting profile into a bulge+disk gives structural parameters in reasonable agreement with observations although kinematics betrays their bar nature. These findings have important implications for galaxy formation models, which have so far ignored or introduced in a very simplified way the effects of nonaxisymmetric instabilities on the morphological evolution of disk galaxies.

Simultaneous ram pressure and tidal stripping; how dwarf spheroidals lost their gas

We perform high-resolution N-body+SPH (smoothed particle hydrodynamics) simulations of gas-rich dwarf galaxy satellites orbiting within a Milky Way-sized halo and study for the first time the combined effects of tides and ram pressure. The structure of the galaxy models and the orbital configurations are chosen in accordance with those expected in a Lambda cold dark matter (ACDM) universe. While tidal stirring of disky dwarfs produces objects whose stellar structure and kinematics resembles that of dwarf spheroidals after a few orbits, ram pressure stripping is needed to entirely remove their gas component. Gravitational tides can aid ram pressure stripping by diminishing the overall potential of the dwarf, but tides also induce bar formation which funnels gas inwards making subsequent stripping more difficult. This inflow is particularly effective when the gas can cool radiatively. Assuming a low density of the hot Galactic corona consistent with observational constraints, dwarfs with V peak 30 km s -1 lose most or all of their gas content only if a heating source keeps the gas extended, partially counteracting the bar-driven inflow. We show that the ionizing radiation from the cosmic ultraviolet (UV) background at z > 2 can provide the required heating. In these objects, most of the gas is removed or becomes ionized at the first pericenter passage, explaining the early truncation of the star formation observed in Draco and Ursa Minor. Galaxies on orbits with larger pericenters and/or falling into the Milky Way halo at lower redshift can retain significant amounts of the centrally concentrated gas. These dwarfs would continue to form stars over a longer period of time, especially close to pericenter passages, as observed in Fornax and other dwarf spheroidal galaxies (dSphs) of the Local Group. The stripped gas breaks up into individual clouds pressure confined by the outer gaseous medium that have masses, sizes and densities comparable to the H I clouds recently discovered around M31.

The Formation of a Realistic Disk Galaxy in Λ-dominated Cosmologies

We simulate the formation of a realistic disk galaxy within the hierarchical scenario of structure formation and study its internal properties to the present epoch. We use a set of smoothed particle hydrodynamic (SPH) simulations, with a high dynamical range and force resolution, that include cooling, star formation, supernovae (SNe) feedback and a redshift-dependent UV background. We compare results from a Λ cold dark matter (ΛCDM) simulation to a Λ warm dark matter (ΛWDM) (2 keV) simulation that forms significantly less small-scale structure. We show how high mass and force resolution in both the gas and dark-matter components play an important role in solving the angular momentum catastrophe claimed from previous simulations of galaxy formation within the hierarchical framework. Hence, a large disk forms without the need of strong energy injection, the z = 0 galaxies lie close to the I band Tully-Fisher relation, and the stellar material in the disk component has a final specific angular momentum equal to 40% and 90% of the dark halo in the ΛCDM and ΛWDM models, respectively. If rescaled to the Milky Way, the ΛCDM galaxy has an overabundance of satellites, with a total mass in the stellar halo 40% of that in the bulge+disk system. The ΛWDM galaxy has a drastically reduced satellite population and a negligible stellar spheroidal component. Encounters with satellites play only a minor role in disturbing the disk. Satellites possess a variety of star formation histories linked to mergers and pericentric passages along their orbit around the primary galaxy. In both cosmologies, the galactic halo retains most of the baryons accreted and builds up a hot gas phase with a substantial X-ray emission. Therefore, while we have been successful in creating a realistic stellar disk in a massive galaxy within the ΛCDM scenario, energy injection emerges as necessary ingredient to reduce the baryon fraction in galactic halos, independent of the cosmology adopted.

Formation of Giant Planets by Fragmentation of Protoplanetary Disks

The evolution of gravitationally unstable protoplanetary gaseous disks has been studied with the use of three-dimensional smoothed particle hydrodynamics simulations with unprecedented resolution. We have considered disks with initial masses and temperature profiles consistent with those inferred for the protosolar nebula and for other protoplanetary disks. We show that long-lasting, self-gravitating protoplanets arise after a few disk orbital periods if cooling is efficient enough to maintain the temperature close to 50 K. The resulting bodies have masses and orbital eccentricities similar to those of detected extrasolar planets.

Rapid Formation of Supermassive Black Hole Binaries in Galaxy Mergers with Gas

Supermassive black holes (SMBHs) are a ubiquitous component of the nuclei of galaxies. It is normally assumed that after the merger of two massive galaxies, a SMBH binary will form, shrink because of stellar or gas dynamical processes, and ultimately coalesce by emitting a burst of gravitational waves. However, so far it has not been possible to show how two SMBHs bind during a galaxy merger with gas because of the difficulty of modeling a wide range of spatial scales. Here we report hydrodynamical simulations that track the formation of a SMBH binary down to scales of a few light years after the collision between two spiral galaxies. A massive, turbulent, nuclear gaseous disk arises as a result of the galaxy merger. The black holes form an eccentric binary in the disk in less than 1 million years as a result of the gravitational drag from the gas rather than from the stars.

Clumps in the Outer Disk by Disk Instability: Why They are Initially Gas Giants and the Legacy of Disruption

Abstract We explore the initial conditions for fragments in the extended regions ( r ≳ 50 AU ) of gravitationally unstable disks. We combine analytic estimates for the fragmentation of spiral arms with 3D SPH simulations to show that initial fragment masses are in the gas giant regime. These initial fragments will have substantial angular momentum, and should form disks with radii of a few AU. We show that clumps will survive for multiple orbits before they undergo a second, rapid collapse due to H 2 dissociation and that it is possible to destroy bound clumps by transporting them into the inner disk. The consequences of disrupted clumps for planet formation, dust processing, and disk evolution are discussed. We argue that it is possible to produce Earth-mass cores in the outer disk during the earliest phases of disk evolution.

Supermassive black hole binaries in gaseous and stellar circumnuclear discs: orbital dynamics and gas accretion

The dynamics of two massive black holes in a rotationally supported nuclear disc of mass M disc = 10 8 M⊙ is explored using N-body/smoothed particle hydrodynamics simulations. Gas and star particles are copresent in the disc. Described by a Mestel profile, the disc has a vertical support provided by turbulence of the gas, and by stellar velocity dispersion. A primary black hole of mass 4 x 10 6 M⊙ is placed at the centre of the disc, while a secondary black hole is set initially on an eccentric corotating orbit in the disc plane. Its mass is in a 1:1, 1:4, and 1:10 ratio, relative to the primary. With this choice, we mimic the dynamics of black hole pairs released in the nuclear region at the end of a gas-rich galaxy merger. It is found that, under the action of dynamical friction, the two black holes form a close binary in ∼10 Myr. The inspiral process is insensitive to the mass fraction in stars and gas present in the disc and is accompanied by the circularization of the orbit. We detail the gaseous mass profile bound to each black hole that can lead to the formation of two small Keplerian discs, weighing ≈2 per cent of the black hole mass, and of size ∼0.01 pc. The mass of the tightly (loosely) bound particles increases (decreases) with time as the black holes spiral into closer and closer orbits. Double active galactic nucleus activity is expected to occur on an estimated time-scale of ≤10 Myr, comparable to the inspiral time-scale. The double nuclear point-like sources that may appear during dynamical evolution will have typical separations of ≤10 pc.