No Arabic abstract
We use a high-temperature chemical network to derive the molecular abundances in axisymmetric accretion disk models around active galactic nuclei (AGNs) within 100 pc using simple radial and vertical density and temperature distributions motivated by more detailed physical models. We explore the effects of X-ray irradiation and cosmic ray ionization on the spatial distribution of the molecular abundances of CO, CN, CS, HCN, HCO+, HC3N, C2H, and c-C3H2 using a variety of plausible disk structures. These simple models have molecular regions with a layer of X-ray dominated regions, a midplane without the strong influence of X-rays, and a high-temperature region in the inner portion with moderate X-ray flux where families of polyynes (C$_{rm n}$H$_{2}$) and cyanopolyynes can be enhanced. For the high midplane density disks we explore, we find that cosmic rays produced by supernovae do not significantly affect the regions unless the star formation efficiency significantly exceeds that of the Milky Way. We highlight molecular abundance observations and ratios that may distinguish among theoretical models of the density distribution in AGN disks. Finally, we assess the importance of the shock crossing time and the accretion time relative to the formation time for various chemical species. Vertical column densities are tabulated for a number of molecular species at both the characteristic shock crossing time and steady state. Although we do not attempt to fit any particular system or set of observations, we discuss our models and results in the context of the nearby AGN NGC 1068.
Protoplanetary disks are the target of many chemical studies (both observational and theoretical) as they contain the building material for planets. Their large vertical and radial gradients in density and temperature make them challenging objects for chemical models. In the outer part of these disks, the large densities and low temperatures provide a particular environment where the binding of species onto the dust grains can be very efficient and can affect the gas-phase chemical composition. We attempt to quantify to what extent the vertical abundance profiles and the integrated column densities of molecules predicted by a detailed gas-grain code are affected by the treatment of the molecular hydrogen physisorption at the surface of the grains. We performed three different models using the Nautilus gas-grain code. One model uses a H2 binding energy on the surface of water (440 K) and produces strong sticking of H2. Another model uses a small binding energy of 23 K (as if there were already a monolayer of H2), and the sticking of H$_2$ is almost negligible. Finally, the remaining model is an intermediate solution known as the encounter desorption mechanism. We show that the efficiency of molecular hydrogen binding (and thus its abundance at the surface of the grains) can have a quantitative effect on the predicted column densities in the gas phase of major species such as CO, CS, CN, and HCN.
Large-scale cosmological simulations of galaxy formation currently do not resolve the densities at which molecular hydrogen forms, implying that the atomic-to-molecular transition must be modeled either on the fly or in postprocessing. We present an improved postprocessing framework to estimate the abundance of atomic and molecular hydrogen and apply it to the IllustrisTNG simulations. We compare five different models for the atomic-to-molecular transition, including empirical, simulation-based, and theoretical prescriptions. Most of these models rely on the surface density of neutral hydrogen and the ultraviolet (UV) flux in the Lyman-Werner band as input parameters. Computing these quantities on the kiloparsec scales resolved by the simulations emerges as the main challenge. We show that the commonly used Jeans length approximation to the column density of a system can be biased and exhibits large cell-to-cell scatter. Instead, we propose to compute all surface quantities in face-on projections and perform the modeling in two dimensions. In general, the two methods agree on average, but their predictions diverge for individual galaxies and for models based on the observed midplane pressure of galaxies. We model the UV radiation from young stars by assuming a constant escape fraction and optically thin propagation throughout the galaxy. With these improvements, we find that the five models for the atomic-to-molecular transition roughly agree on average but that the details of the modeling matter for individual galaxies and the spatial distribution of molecular hydrogen. We emphasize that the estimated molecular fractions are approximate due to the significant systematic uncertainties.
We perform collisionless N-body simulations to investigate the evolution of the structural and kinematical properties of simulated thick disks induced by the growth of an embedded thin disk. The thick disks used in the present study originate from cosmologically-common 5:1 encounters between initially-thin primary disk galaxies and infalling satellites. The growing thin disks are modeled as static gravitational potentials and we explore a variety of growing-disk parameters that are likely to influence the response of thick disks. We find that the final thick-disk properties depend strongly on the total mass and radial scale-length of the growing thin disk, and much less sensitively on its growth timescale and vertical scale-height as well as the initial sense of thick-disk rotation. Overall, the growth of an embedded thin disk can cause a substantial contraction in both the radial and vertical direction, resulting in a significant decrease in the scale-lengths and scale-heights of thick disks. Kinematically, a growing thin disk can induce a notable increase in the mean rotation and velocity dispersions of thick-disk stars. We conclude that the reformation of a thin disk via gas accretion may play a significant role in setting the structure and kinematics of thick disks, and thus it is an important ingredient in models of thick-disk formation.
We study outflows driven by Active Galactic Nuclei (AGNs) using high- resolution simulations of idealized z=2 isolated disk galaxies. Episodic accretion events lead to outflows with velocities >1000 km/s and mass outflow rates up to the star formation rate (several tens of Msun/yr). Outflowing winds escape perpendicular to the disk with wide opening angles, and are typically asymmetric (i.e. unipolar) because dense gas above or below the AGN in the resolved disk inhibits outflow. Owing to rapid variability in the accretion rates, outflowing gas may be detectable even when the AGN is effectively off. The highest velocity outflows are sometimes, but not always, concentrated within 2-3 kpc of the galactic center during the peak accretion. With our purely thermal AGN feedback model -- standard in previous literature -- the outflowing material is mostly hot (10^6 K) and diffuse (nH<10^(-2) cm-3), but includes a cold component entrained in the hot wind. Despite the powerful bursts and high outflow rates, AGN feedback has little effect on the dense gas in the galaxy disk. Thus AGN-driven outflows in our simulations do not cause rapid quenching of star-formation, although they may remove significant amounts of gas of long (>Gyr) timescales.
Passive early-type galaxies dominate cluster cores at z $lesssim$1.5. At higher redshift, cluster core galaxies are observed to have still on-going star-formation, fuelled by cold molecular gas. We measure the molecular gas reservoir of the central region around the radio-loud AGN in the cluster CARLA J1103+3449 at z=1.44 with NOEMA. The AGN synchrotron emission dominates the continuum emission at 94.48 GHz, and we measure its flux at the AGN position and at the position of two radio jets. Combining our measurements with published results over the range 4.71 GHz-94.5 GHz, we obtain a flat spectral index $alpha = 0.14 pm 0.03$ for the AGN core emission, and a steeper index $alpha = 1.43 pm 0.04$ and $alpha = 1.15 pm 0.04$ at positions close to the western and eastern lobe, respectively. The total spectral index is $alpha = 0.92 pm 0.02$ over the range 73.8 MHz-94.5 GHz. We detect two CO(2-1) emission lines, both blue-shifted with respect to the AGN. Their emission corresponds to two regions, ~17 kpc south-east and ~14 kpc south-west of the AGN, not associated with galaxies. In these two regions, we find a total massive molecular gas reservoir of $M_{gas}$ = 3.9 $pm$ 0.4 $10^{10} M_{odot}$, which dominates (~ 60%) the central total molecular gas reservoir. These results can be explained by massive cool gas flows in the center of the cluster. The AGN early-type host is not yet quenched; its star formation rate is consistent with being on the main sequence of star-forming galaxies in the field (SFR~30-140 $M_{odot}$/yr), and the cluster core molecular gas reservoir is expected to feed the AGN and the host star-formation before quiescence. The other cluster confirmed members show star formation rates at ~2 $sigma$ below the field main sequence at similar redshifts and do not have molecular gas masses larger than galaxies of similar stellar mass in the field.