We present novel 3D multi-scale SPH simulations of gas-rich galaxy mergers between the most massive galaxies at $z sim 8 - 10$, designed to scrutinize the direct collapse formation scenario for massive black hole seeds proposed in citet{mayer+10}. Th
e simulations achieve a resolution of 0.1 pc, and include both metallicity-dependent optically-thin cooling and a model for thermal balance at high optical depth. We consider different formulations of the SPH hydrodynamical equations, including thermal and metal diffusion. When the two merging galaxy cores collide, gas infall produces a compact, optically thick nuclear disk with densities exceeding $10^{-10}$ g cm$^3$. The disk rapidly accretes higher angular momentum gas from its surroundings reaching $sim 5$ pc and a mass of $gtrsim 10^9$ $M_{odot}$ in only a few $10^4$ yr. Outside $gtrsim 2$ pc it fragments into massive clumps. Instead, supersonic turbulence prevents fragmentation in the inner parsec region, which remains warm ($sim 3000-6000$ K) and develops strong non-axisymmetric modes that cause prominent radial gas inflows ($> 10^4$ $M_{odot}$ yr$^{-1}$), forming an ultra-dense massive disky core. Angular momentum transport by non-axisymmetric modes should continue below our spatial resolution limit, quickly turning the disky core into a supermassive protostar which can collapse directly into a massive black hole of mass $10^8-10^9$ $M_{odot}$ via the relativistic radial instability. Such a cold direct collapse explains naturally the early emergence of high-z QSOs. Its telltale signature would be a burst of gravitational waves in the frequency range $10^{-4} - 10^{-1}$ Hz, possibly detectable by the planned eLISA interferometer.
The Hubble sequence is a common classification scheme for the structure of galaxies. Despite the tremendous usefulness of this diagnostic, we still do not fully understand when, where, and how this morphological ordering was put in place. Here, we in
vestigate the morphological evolution of a sample of 22 high redshift ($zgeq3$) galaxies extracted from the Argo simulation. Argo is a cosmological zoom-in simulation of a group-sized halo and its environment. It adopts the same high resolution ($sim10^4$ M$_odot$, $sim100$ pc) and sub-grid physical model that was used in the Eris simulation but probes a sub-volume almost ten times bigger with as many as 45 million gas and star particles in the zoom-in region. Argo follows the early assembly of galaxies with a broad range of stellar masses ($log M_{star}/{rm M}_{odot}sim8-11$ at $zsimeq3$), while resolving properly their structural properties. We recover a diversity of morphologies, including late-type/irregular disc galaxies with flat rotation curves, spheroid dominated early-type discs, and a massive elliptical galaxy, already established at $zsim3$. We identify major mergers as the main trigger for the formation of bulges and the steepening of the circular velocity curves. Minor mergers and non-axisymmetric perturbations (stellar bars) drive the bulge growth in some cases. The specific angular momenta of the simulated disc components fairly match the values inferred from nearby galaxies of similar $M_{star}$ once the expected redshift evolution of disc sizes is accounted for. We conclude that morphological transformations of high redshift galaxies of intermediate mass are likely triggered by processes similar to those at low redshift and result in an early build-up of the Hubble sequence.
We study the dynamics of massive black hole pairs in clumpy gaseous circumnuclear disks. We track the orbital decay of the light, secondary black hole $M_{bullet2}$ orbiting around the more massive primary at the center of the disk, using $N$-body/sm
oothed particle hydrodynamic simulations. We find that the gravitational interaction of $M_{bullet2}$ with massive clumps $M_{rm cl}$ erratically perturbs the otherwise smooth orbital decay. In close encounters with massive clumps, gravitational slingshots can kick the secondary black hole out of the disk plane. The black hole moving on an inclined orbit then experiences the weaker dynamical friction of the stellar background, resulting in a longer orbital decay timescale. Interactions between clumps can also favor orbital decay when the black hole is captured by a massive clump which is segregating toward the center of the disk. The stochastic behavior of the black hole orbit emerges mainly when the ratio $M_{bullet2}/M_{rm cl}$ falls below unity, with decay timescales ranging from $sim1$ to $sim50$ Myr. This suggests that describing the cold clumpy phase of the inter-stellar medium in self-consistent simulations of galaxy mergers, albeit so far neglected, is important to predict the black hole dynamics in galaxy merger remnants.
