Observations show a prevalence of high-redshift galaxies with large stellar masses and predominantly passive stellar populations. A variety of processes have been suggested that could reduce the star ...formation in such galaxies to observed levels, including quasar mode feedback, virial shock heating, or galactic winds driven by stellar feedback. However, the main quenching mechanisms have yet to be identified. Here we study the origin of star formation quenching using Argo, a cosmological, hydrodynamical zoom-in simulation that follows the evolution of a massive galaxy at z ~ 2. This simulation adopts the same subgrid recipes of the Eris simulations, which have been shown to form realistic disc galaxies, and, in one version, adopts also a mass and spatial resolution identical to Eris. The resulting galaxy has properties consistent with those of observed, massive (M* ~ 10... M...) galaxies at z ~ 2 and with abundance matching predictions. Our models do not include active galactic nuclei (AGN) feedback indicating that supermassive black holes likely play a subordinate role in determining masses and sizes of massive galaxies at high-z. The specific star formation rate (sSFR) of the simulated galaxy matches the observed M*-sSFR relation at early times. This period of smooth stellar mass growth comes to a sudden halt at z = 3.5 when the sSFR drops by almost an order of magnitude within a few hundred Myr. The suppression is initiated by a levelling off and a subsequent reduction of the cool gas accretion rate on to the galaxy, and not by feedback processes. This 'cosmological starvation' occurs as the parent dark matter halo switches from a fast collapsing mode to a slow accretion mode. Additional mechanisms, such as perhaps radio mode feedback from an AGN, are needed to quench any residual star formation of the galaxy and to maintain a low sSFR until the present time. (ProQuest: ... denotes formulae/symbols omitted.)
We analyze the formation and evolution of the stellar components in "Eris," a 120 pc resolution cosmological hydrodynamic simulation of a late-type spiral galaxy. The simulation includes the effects ...of a uniform UV background, a delayed-radiative-cooling scheme for supernova feedback, and a star formation recipe based on a high gas density threshold. It allows a detailed study of the relative contributions of "in-situ" (within the main host) and "ex-situ" (within satellite galaxies) star formation to each major Galactic component in a close Milky Way analog. We investigate these two star-formation channels as a function of galactocentric distance, along different lines of sight above and along the disk plane, and as a function of cosmic time. We find that: (1) approximately 70% of today's stars formed in-situ; (2) more than two thirds of the ex-situ stars formed within satellites after infall; (3) the majority of ex-situ stars are found today in the disk and in the bulge; (4) the stellar halo is dominated by ex-situ stars, whereas in-situ stars dominate the mass profile at distances <, ~5 kpc from the center at high latitudes; and (5) approximately 25% of the inner, r <, ~ 20 kpc, halo is composed of in-situ stars that have been displaced from their original birth sites during Eris' early assembly history.
We study the effect of sub-grid physics, galaxy mass, structural parameters and resolution on the fragmentation of gas-rich galaxy discs into massive star-forming clumps. The initial conditions are ...set up with the aid of the ARGO cosmological hydrodynamical simulation. Blast-wave feedback does not suppress fragmentation, but reduces both the number of clumps and the duration of the unstable phase. Once formed, bound clumps cannot be destroyed by our feedback model. Widespread fragmentation is promoted by high gas fractions and low halo concentrations. Yet giant clumps M > 108 M⊙ lasting several hundred Myr are rare and mainly produced by clump–clump mergers. They occur in massive discs with maximum rotational velocities V
max > 250 km s−1 at z ∼ 2, at the high-mass end of the observed galaxy population at those redshifts. The typical gaseous and stellar masses of clumps in all runs are in the range ∼107–108 M⊙ for galaxies with disc mass in the range 1010–8 × 1010 M⊙. Clumps sizes are usually in the range 100–400 pc, in agreement with recent clump observations in lensed high-z galaxies. We argue that many of the giant clumps identified in observations are not due to in situ fragmentation, or are the result of blending of smaller structures owing to insufficient resolution. Using an analytical model describing local collapse inside spiral arms, we can predict the characteristic gaseous masses of clumps at the onset of fragmentation (∼3–5 × 107 M⊙) quite accurately, while the conventional Toomre mass overestimates them. Due to their moderate masses, clumps which migrate to the centre have marginal effect on bulge growth.
We carry out simulations of gravitationally unstable disks using smoothed particle hydrodynamics (SPH) and the novel Lagrangian meshless finite mass (MFM) scheme in the GIZMO code. Our aim is to ...understand the cause of the nonconvergence of the cooling boundary for fragmentation reported in the literature. We run SPH simulations with two different artificial viscosity implementations and compare them with MFM, which does not employ any artificial viscosity. With MFM we demonstrate convergence of the critical cooling timescale for fragmentation at . Nonconvergence persists in SPH codes. We show how the nonconvergence problem is caused by artificial fragmentation triggered by excessive dissipation of angular momentum in domains with large velocity derivatives. With increased resolution, such domains become more prominent. Vorticity lags behind density, due to numerical viscous dissipation in these regions, promoting collapse with longer cooling times. Such effect is shown to be dominant over the competing tendency of artificial viscosity to diminish with increasing resolution. When the initial conditions are first relaxed for several orbits, the flow is more regular, with lower shear and vorticity in nonaxisymmetric regions, aiding convergence. Yet MFM is the only method that converges exactly. Our findings are of general interest, as numerical dissipation via artificial viscosity or advection errors can also occur in grid-based codes. Indeed, for the FARGO code values of significantly higher than our converged estimate have been reported in the literature. Finally, we discuss implications for giant planet formation via disk instability.
Satellites of giant planets have been thought to form in gaseous circumplanetary disks (CPDs) during the late planet-formation phase, but it was unknown whether or not smaller-mass planets such as ...the ice giants could form such disks, and thus moons, there. We combined radiative hydrodynamical simulations with satellite population synthesis to investigate the question in the case of Uranus and Neptune. For both ice giants we found that a gaseous CPD is created at the end of their formation. The population synthesis confirmed that Uranian-like, icy, prograde satellite system could form in these CPDs within a couple of 105 yr. This means that Neptune could have a Uranian-like moon system originally that was wiped away by the capture of Triton. Furthermore, the current moons of Uranus can be reproduced by our model without the need for planet-planet impact to create a debris disk for the moons to grow. These results highlight that even ice giants-among the most common mass category of exoplanets-can also form satellites, opening a way to a potentially much larger population of exomoons than previously thought.
ABSTRACT Supermassive black holes (SMBHs) are ubiquitous in galaxies with a sizable mass. It is expected that a pair of SMBHs originally in the nuclei of two merging galaxies would form a binary and ...eventually coalesce via a burst of gravitational waves. So far, theoretical models and simulations, focusing only on limited phases of the orbital decay of SMBHs under idealized conditions of the galaxy hosts, have been unable to directly predict the SMBH merger timescale from ab-initio galaxy formation theory. The predicted SMBH merger timescales are long, of order Gyrs, which could be problematic for future gravitational wave (GW) searches. Here, we present the first multi-scale ΛCDM cosmological simulation that follows the orbital decay of a pair of SMBHs in a merger of two typical massive galaxies at , all the way to the final coalescence driven by GW emission. The two SMBHs, with masses , settle quickly in the nucleus of the merger remnant. The remnant is triaxial and extremely dense due to the dissipative nature of the merger and the intrinsic compactness of galaxies at high redshift. Such properties naturally allow a very efficient hardening of the SMBH binary. The SMBH merger occurs in only ∼10 Myr after the galactic cores have merged, which is two orders of magnitude smaller than the Hubble time.
We present a semi-analytical population synthesis model of protoplanetary clumps formed by disk instability at radial distances of 80-120 au. Various clump density profiles, initial mass functions, ...protoplanetary disk models, stellar masses, and gap opening criteria are considered. When we use more realistic gap opening criteria, we find that gaps open only rarely, which strongly affects clump survival rates and their physical properties (mass, radius, and radial distance). The inferred surviving population is then shifted toward less massive clumps at smaller radial distances. We also find that populations of surviving clumps are very sensitive to the model assumptions and used parameters. Depending on the chosen parameters, the protoplanets occupy a mass range between 0.01 and 16 MJ and may either orbit close to the central star or as far out as 75 au, with a sweet spot at 10-30 au for the massive ones. However, in all of the cases we consider, we find that massive giant planets at very large radial distances are rare, in qualitative agreement with current direct imaging surveys. We conclude that caution should be taken in deriving population synthesis models as well as when comparing the models' results with observations.
We study the orbital decay of a pair of massive black holes (BHs) with masses and 107 , using hydrodynamical simulations of circumnuclear disks (CNDs) with the alternating presence of sub-grid ...physics, such as radiative cooling, star formation, supernova feedback, BH accretion, and BH feedback. In the absence of such processes, the orbit of the secondary BH decays over timescales of to the center of the CND, where the primary BH resides. When strong dissipation operates in CNDs, fragmentation into massive objects the size of giant molecular clouds with densities in the range 104-107 amu cm−3 occurs, causing stochastic torques and hits that can eject the secondary BH from the midplane. Outside the plane, the low-density medium provides only weak drag, and the BH return is governed by inefficient dynamical friction. In rare cases, clump-BH interactions can lead to a faster decay. Feedback processes lead to outflows, but do not significantly change the overall density of the CND midplane. However, with a spherically distributed BH feedback, a hot bubble is generated behind the secondary, which almost shuts off dynamical friction. We dub this phenomenon "wake evacuation." It leads to delays in the decay, possibly of . We discuss the non-trivial implications on the discovery space of the eLISA telescope. Our results suggest that the largest uncertainty in predicting BH merger rates lies in the potentially wide variety of galaxy host systems, with different degrees of gas dissipation and heating, yielding decay timescales from to .
ABSTRACT
Observations of hyper-luminous quasars at z>6 reveal the rapid growth of supermassive black holes (SMBHs ${\gt}10^9 \,\rm M_{\odot }$) whose origin is still difficult to explain. Their ...progenitors may have formed as remnants of massive, metal-free stars (light seeds), via stellar collisions (medium-weight seeds) and/or massive gas clouds direct collapse (heavy seeds). In this work, we investigate for the first time the relative role of these three seed populations in the formation of z>6 SMBHs within an Eddington-limited gas accretion scenario. To this aim, we implement in our semi-analytical data-constrained model a statistical description of the spatial fluctuations of Lyman–Werner (LW) photodissociating radiation and of metal/dust enrichment. This allows us to set the physical conditions for black hole seeds formation, exploring their relative birth rate in a highly biased region of the Universe at z>6. We find that the inclusion of medium-weight seeds does not qualitatively change the growth history of the first SMBHs: although less massive seeds (${\lt}10^3\, \rm M_\odot$) form at a higher rate, the mass growth of a ${\sim}10^9\, \rm M_\odot$ SMBH at z<15 is driven by efficient gas accretion (at a sub-Eddington rate) on to its heavy progenitors ($10^5\, \rm M_\odot$). This conclusion holds independently of the critical level of LW radiation and even when medium-weight seeds are allowed to form in higher metallicity galaxies, via the so-called supercompetitive accretion scenario. Our study suggests that the genealogy of z∼6 SMBHs is characterized by a rich variety of BH progenitors, which represent only a small fraction (${\lt} 10{-}20{{\ \rm per\ cent}}$) of all the BHs that seed galaxies at z>15.