No Arabic abstract
We utilize zoom-in cosmological simulations to study the nature of violent disc instability (VDI) in clumpy galaxies at high redshift, $z=1$--$5$. Our simulated galaxies are not in the ideal state assumed in Toomre instability, of linear fluctuations in an isolated, uniform, rotating disk. There, instability is characterised by a $Q$ parameter below unity, and lower when the disk is thick. Instead, the high-redshift discs are highly perturbed. Over long periods they consist of non-linear perturbations, compact massive clumps and extended structures, with new clumps forming in inter-clump regions. This is while the galaxy is subject to frequent external perturbances. We compute the local, two-component $Q$ parameter for gas and stars, smoothed on a $sim1~{rm kpc}$ scale to capture clumps of $10^{8-9}~{rm M}_odot$. The $Q<1$ regions are confined to collapsed clumps due to the high surface density there, while the inter-clump regions show $Q$ significantly higher than unity. Tracing the clumps back to their relatively smooth Lagrangian patches, we find that $Q$ prior to clump formation typically ranges from unity to a few. This is unlike the expectations from standard Toomre instability. We discuss possible mechanisms for high-$Q$ clump formation, e.g. rapid turbulence decay leading to small clumps that grow by mergers, non-axisymmetric instability, or clump formation induced by non-linear perturbations in the disk. Alternatively, the high-$Q$ non-linear VDI may be stimulated by the external perturbations such as mergers and counter-rotating streams. The high $Q$ may represent excessive compressive modes of turbulence, possibly induced by tidal interactions.
Many massive star forming disc galaxies in the redshift range 3 to 0.5 are observed to have a clumpy morphology showing giant clumps of size $sim$1 kpc and masses of about $10^7M_{odot}$ to $10^{10} M_{odot}$. The nature and fate of these giant clumps is still under debate. In this work we use 19 high-resolution simulations of disc galaxies from the NIHAO sample to study the formation and the evolution of clumps in the discs of high redshift galaxies. We use mock HST - CANDELS observations created with the radiative transfer code GRASIL-3D to carry out, for the first time, a quantitative comparison of the observed fraction of clumpy galaxies and its evolution with redshift with simulations. We find a good agreement between the observed clumpy fraction and the one of the NIHAO galaxies. We find that dust attenuation can suppress intrinsically bright clumps and enhance less luminous ones. In our galaxy sample we only find clumps in light (u-band) from young stars but not in stellar mass surface density maps. This means that the NIHAO sample does not show clumpy stellar discs but rather a clumpy light distribution originating from clumpy star formation events. The clumps found in the NIHAO sample match observed age/color gradients as a function of distance from the galaxy center but they show no sign of inward migration. Clumps in our simulations disperse on timescales of a about a hundred Myr and their contribution to bulge growth is negligible.
We have developed a novel Markov Chain Mote Carlo (MCMC) chemical painting technique to explore possible radial and vertical metallicity gradients for the thick disc progenitor. In our analysis we match an N-body simulation to the data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. We assume that the thick disc has a constant scale-height and has completed its formation at an early epoch, after which time radial mixing of its stars has taken place. Under these assumptions, we find that the initial radial metallicity gradient of the thick disc progenitor should not be negative, but either flat or even positive, to explain the current negative vertical metallicity gradient of the thick disc. Our study suggests that the thick disc was built-up in an inside-out and upside-down fashion, and older, smaller and thicker populations are more metal poor. In this case, star forming discs at different epochs of the thick disc formation are allowed to have different radial metallicity gradients, including a negative one, which helps to explain a variety of slopes observed in high redshift disc galaxies. This scenario helps to explain the positive slope of the metallicity-rotation velocity relation observed for the Galactic thick disc. On the other hand, radial mixing flattens the slope of an existing gradient.
The disc structure of the Milky Way is marked by a chemical dichotomy, with high-alpha and low-alpha abundance sequences, traditionally identified with the geometric thick and thin discs. This identification is aided by the old ages of the high-alpha stars, and lower average ages of the low-alpha ones. Recent large scale surveys such as APOGEE have provided a wealth of data on this chemical structure, including showing that an identification of chemical and geometric thick discs is not exact, but the origin of the chemical dichotomy has remained unclear. Here we demonstrate that a dichotomy arises naturally if the early gas-rich disc fragments, leading to some fraction of the star formation occurring in clumps of the type observed in high-redshift galaxies. These clumps have high star formation rate density. They, therefore, enrich rapidly, moving from the low-alpha to the high-alpha sequence, while more distributed star formation produces the low-alpha sequence. We demonstrate that this model produces a chemically-defined thick disc that has many of the properties of the Milky Ways thick disc. Because clump formation is common in high redshift galaxies, we predict that chemical bimodalities are common in massive galaxies.
Massive stellar clumps in high redshift galaxies interact and migrate to the center to form a bulge and exponential disk in <1 Gyr. Here we consider the fate of intermediate mass black holes (BHs) that might form by massive-star coalescence in the dense young clusters of these disk clumps. We find that the BHs move inward with the clumps and reach the inner few hundred parsecs in only a few orbit times. There they could merge into a supermassive BH by dynamical friction. The ratio of BH mass to stellar mass in the disk clumps is approximately preserved in the final ratio of BH to bulge mass. Because this ratio for individual clusters has been estimated to be ~10^{-3}, the observed BH-to-bulge mass ratio results. We also obtain a relation between BH mass and bulge velocity dispersion that is compatible with observations of present-day galaxies.
A strong signature of a circumstellar disc around a high-mass protostar has been inferred from high resolution methanol maser observations in NGC7538-IRS1 N. This interpretation has however been challenged with a bipolar outflow proposed as an alternative explanation. We compare the two proposed scenarios for best consistency with the observations. Using a newly developed formalism we model the optical depth of the maser emission at each observed point in the map and LOS velocity for the two scenarios. We find that if the emission is symmetric around a central peak in both space and LOS velocity then it has to arise from an edge-on disc in sufficiently fast differential rotation. Disc models successfully fit ~100 independent measurement points in position-velocity space with 4 free parameters to an overall accuracy of 3-4%. Solutions for Keplerian rotation require a central mass of at least 4 solar masses. Close to best-fitting models are obtained if Keplerian motion is assumed around a central mass equaling ~30 solar masses as inferred from other observations. In contrast we find that classical bipolar outflow models cannot fit the data, although could be applicable in other sources. Our results strongly favour the differentially rotating disc hypothesis to describe the main feature of the 12.2 (and 6.7) GHz methanol maser emission in NGC7538 IRS1 N. Furthermore, for Keplerian rotation around a ~30 solar masses protostar we predict the position and velocity at which tangentially amplified masers should be detected in high dynamic range observations. [abridged]