Source author record

Isaac Shlosman

Isaac Shlosman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2022arXiv

Modeling Evolution of Galactic Bars at Cosmic Dawn

We study evolution of galactic bars using suite of very high-resolution zoom-in cosmological simulations of galaxies at z ~ 9-2. Our models were chosen to lie within similar mass DM halos, log(Mvir/Mo) ~ 11.65 +- 0.05, at z=6, 4, and 2, in high and low overdensity environments. We apply two galactic wind feedback mechanisms for each model. All galaxies develop sub-kpc stellar bars differing in their properties. We find that (1) The high-z bars form in response to various perturbations: mergers, close flybys, cold accretion inflows along the cosmological filaments, etc.; (2) These bars account for large-mass fraction of galaxies; (3) Bars display large corotation-to-bar-size ratios, and are weaker compared to their low-redshift counterparts, by measuring their Fourier amplitudes, and are very gas-rich; (4) Their pattern speed does not exhibit monotonic decline with time due to braking against DM, as at low z; (5) Bar properties, including their stellar population (SFRs and metal enrichment) depend sensitively on prevailing feedback; (6) Finally, we find that bars can weaken substantially during cosmological evolution, becoming weak oval distortions -- hence bars are destroyed and reformed multiple times unlike their low-z counterparts. In all cases, bars in our simulations have been triggered by interactions. In summary, stellar bars appear to be not only contemporary phenomenon, but based on increased frequency of mergers, flybys and the strength of cold accretion flows at high z, we expect them to be ubiquitous at redshifts > 2 -- the epoch of rapid galaxy growth and larger stellar dispersion velocities.

preprint2021arXiv

Emergence of Galactic Morphologies at Cosmic Dawn: Input from Numerical Modeling

We employ high-resolution zoom-in cosmological simulations to analyze the emerging morphology of galaxies in dark matter halos at redshifts z > 2. We choose DM halos of similar masses of log (Mvir/Mo) ~11.65 +- 0.05 at the target redshifts of z_f = 6, 4 and 2. The rationale for this choice, among others, allows us to analyze how the different growth rate in these halos propagates down to galaxy scales. Halos were embedded in high or low overdensity regions, and two different versions of a galactic wind feedback have been employed. Our main results are: (1) Although our galaxies evolve in different epochs, their global parameters remain within a narrow range. Their morphology, kinematics and stellar populations differ substantially, yet all of them host sub-kpc stellar bars; (2) The SFRs appear higher for larger z_f, in tandem with their energy and momentum feedback; (3) The stellar kinematics allowed separation of bulge from the stellar spheroid. The existence of disk-like bulges has been revealed based on stellar surface density and photometry, but displayed a mixed disk-like and classical bulges based on their kinematics. The bulge-to-total mass ratios appear independent of the last merger time for all z_f. The stellar spheroid-to-total mass ratios of these galaxies lie in the range of ~0.5-0.8; (4) The synthetic redshifted, pixelized and PSF-degraded JWST images allow to detect stellar disks at all z_f. Some bars disappear in degraded images, but others remain visible; (5) Based on the kinematic decomposition, for stellar disks separated from bulges and spheroids. we observe that rotational support in disks depends on the feedback type, but increases with decreasing z_f; (6) Finally, the ALMA images detect disks at all z_f, but their spiral structure is only detectable in z_f=2 galaxies.

preprint2020arXiv

Direct collapse to supermassive black hole seeds: the critical conditions for suppression of $\rm H_2$ cooling

Observations of high-redshift quasars imply the presence of supermassive black holes already at z~ 7.5. An appealing and promising pathway to their formation is the direct collapse scenario of a primordial gas in atomic-cooling haloes at z ~ 10 - 20, when the $\rm H_2$ formation is inhibited by a strong background radiation field, whose intensity exceeds a critical value, $J_{\rm crit}$. To estimate $J_{\rm crit}$, typically, studies have assumed idealized spectra, with a fixed ratio of $\rm H_{2}$ photo-dissociation rate $k_{\rm H_2}$ to the $\rm H^-$ photo-detachment rate $k_{\rm H^-}$. This assumption, however, could be too narrow in scope as the nature of the background radiation field is not known precisely. In this work we argue that the critical condition for suppressing the $\rm H_2$ cooling in the collapsing gas could be described in a more general way by a combination of $k_{\rm H_2}$ and $k_{\rm H^-}$ parameters. By performing a series of cosmological zoom-in simulations with an encompassing set of $k_{\rm H_2}$ and $k_{\rm H^-}$, we examine the gas flow by following evolution of basic parameters of the accretion flow. We test under what conditions the gas evolution is dominated by $\rm H_{2}$ and/or atomic cooling. We confirm the existence of a critical curve in the $k_{\rm H_2}-k_{\rm H^-}$ plane, and provide an analytical fit to it. This curve depends on the conditions in the direct collapse, and reveals domains where the atomic cooling dominates over the molecular cooling. Furthermore, we have considered the effect of $\rm H_{2}$ self-shielding on the critical curve, by adopting three methods for the effective column density approximation in $\rm H_{2}$. We find that the estimate of the characteristic length-scale for shielding can be improved by using $λ_{\rm Jeans25}$, which is 0.25 times that of the local Jeans length.

preprint2016arXiv

AGN jet power and feedback characterised by Bondi accretion in brightest cluster galaxies

(Abridged) We propose a new method to estimate the Bondi (hot gas) accretion rates, \dot M_B, onto the supermassive black holes (SMBHs) at the centres of elliptical galaxies. It can be applied even if the Bondi radius is not well-resolved in X-ray observations. This method is based on two simple assumptions: (1) hot gas outside the Bondi radius is in nearly a hydrostatic equilibrium, and (2) the gas temperature near the galaxy centre is close to the virial temperature. We apply this method to 28 bright elliptical galaxies in nearby galaxy clusters (27 of them are the brightest cluster galaxies; BCGs). We find a correlation between the Bondi accretion rates and the power of jets associated with the SMBHs over four orders of magnitude in \dot M_B. For most galaxies, the accretion rates are large enough to account for the jet powers. Our results indicate that the feedback from the active galactic nuclei (AGN) correlates with the properties of the hot gas surrounding the SMBHs. We also find that more massive SMBHs in BCGs tend to have larger specific growth rates. This may explain the hyper masses (~ 10^10 M_sun) of some of the SMBHs. Comparison between the accretion rates and the X-ray luminosities of the AGN suggests that the AGN in the BCGs are extremely radiatively inefficient compared with X-ray binaries in the Milky Way, even when their Eddington accretion ratio, \dot M_B/\dot M_Edd, exceeds 0.01. Lastly, we find a tight correlation between the Bondi accretion rates and the X-ray luminosities of cool cores. Their relation is linear and the power generated by the Bondi accretion is large enough to compensate the radiative cooling of the cool cores. Although the `classical' Bondi accretion model is a greatly oversimplified one, the correlations we find here demonstrate that the accretion onto the SMBHs reflects broadly the properties of the Bondi accretion.

preprint2016arXiv

AGN Jet Power, Formation of X-ray Cavities, and FR I/II Dichotomy in Galaxy Clusters

We investigate the ability of jets in active galactic nuclei (AGNs) to break out of the ambient gas with sufficiently large advance velocities. Using observationally estimated jet power, we analyze 28 bright elliptical galaxies in nearby galaxy clusters. Because the gas density profiles in the innermost regions of galaxies have not been resolved so far, we consider two extreme cases for temperature and density profiles. We also follow two types of evolution for the jet cocoons: being driven by the pressure inside the cocoon (Fanaroff-Riley [FR] I type), and being driven by the jet momentum (FR II type). Our main result is that regardless of the assumed form of density profiles, jets with observed powers of <~ 10^{44} erg s^{-1} are not powerful enough to evolve as FR II sources. Instead, they evolve as FR I sources and appear to be decelerated below the buoyant velocities of the cocoons when jets were propagating through the central dense regions of the host galaxies. This explains the reason why FR I sources are more frequent than FR II sources in clusters. Furthermore, we predict the sizes of X-ray cavities from the observed jet powers and compare them with the observed ones --- they are consistent within a factor of two if the FR I type evolution is realized. Finally, we find that the jets with a power >~ 10^{44} erg s^{-1} are less affected by the ambient medium, and some of them, but not all, could serve as precursors of the FR II sources.

preprint2016arXiv

Direct Collapse to Supermassive Black Hole Seeds: Comparing the AMR and SPH Approaches

We provide detailed comparison between the AMR code Enzo-2.4 and the SPH/N- body code GADGET-3 in the context of isolated or cosmological direct baryonic collapse within dark matter (DM) halos to form supermassive black holes. Gas flow is examined by following evolution of basic parameters of accretion flows. Both codes show an overall agreement in the general features of the collapse, however, many subtle differences exist. For isolated models, the codes increase their spatial and mass resolutions at different pace, which leads to substantially earlier collapse in SPH than in AMR cases due to higher gravitational resolution in GADGET-3. In cosmological runs, the AMR develops a slightly higher baryonic resolution than SPH during halo growth via cold accretion permeated by mergers. Still, both codes agree in the buildup of DM and baryonic structures. However, with the onset of collapse, this difference in mass and spatial resolution is amplified, so evolution of SPH models begins to lag behind. Such a delay can have effect on formation/destruction rate of H2 due to UV background, and on basic properties of host halos. Finally, isolated non-cosmological models in spinning halos, with spin parameter λ ~ 0.01 - 0.07, show delayed collapse for greater λ, but pace of this increase is faster for AMR. Within our simulation setup, GADGET-3 requires significantly larger computational resources than Enzo- 2.4 during collapse, and needs similar resources, during the pre-collapse, cosmological structure formation phase. Yet it benefits from substantially higher gravitational force and hydrodynamic resolutions, except at the end of collapse.

preprint2016arXiv

The Baryon Cycle at High Redshifts: Effects of Galactic Winds on Galaxy Evolution in Overdense and Average Regions

We employ high-resolution cosmological zoom-in simulations focusing on a high-sigma peak and an average cosmological field at $z\sim 6-12$, in order to investigate the influence of environment and baryonic feedback on galaxy evolution in the reionization epoch. Strong feedback, e.g., galactic winds, caused by elevated star formation rates (SFRs) is expected to play an important role in this evolution. We compare different outflow prescriptions: (i) constant wind velocity (CW), (ii) variable wind scaling with galaxy properties (VW), and (iii) no outflows (NW). The overdensity leads to accelerated evolution of dark matter and baryonic structures, absent in the "normal" region, and to shallow galaxy stellar mass functions at the low-mass end. Although CW shows little dependence on both environments, the more physically motivated VW model does exhibit this effect. In addition, VW can reproduce the observed specific SFR (sSFR) and the sSFR-stellar mass relation, which CW and NW fail to satisfy simultaneously. Winds also differ substantially in affecting the state of the intergalactic medium (IGM). The difference lies in volume-filling factor of hot, high-metallicity gas which is near unity for CW, while it remains confined in massive filaments for VW, and locked up in galaxies for NW. Such gas is nearly absent in the normal region. Although all wind models suffer from deficiencies, the VW model seems to be promising in correlating the outflow properties to those of host galaxies. Further constraints on the state of the IGM at high-$z$ are needed to separate different wind models.

preprint2015arXiv

Observational Properties of Simulated Galaxies in Overdense and Average Regions at High Redshifts z= 6-12

We use high-resolution zoom-in cosmological simulations of galaxies of Romano-Diaz et al., post-processing them with a panchromatic three-dimensional radiation transfer code to obtain the galaxy UV luminosity function (LF) at z ~ 6-12. The galaxies are followed in a rare, heavily overdense region within a ~ 5-sigma density peak, which can host high-z quasars, and in an average density region, down to the stellar mass of M_star ~ 4* 10^7 Msun. We find that the overdense regions evolve at a substantially accelerated pace --- the most massive galaxy has grown to M_star ~ 8.4*10^10 Msun by z = 6.3, contains dust of M_dust~ 4.1*10^8 Msun, and is associated with a very high star formation rate, SFR ~ 745 Msun/yr.The attained SFR-M_star correlation results in the specific SFR slowly increasing with M_star. Most of the UV radiation in massive galaxies is absorbed by the dust, its escape fraction f_esc is low, increasing slowly with time. Galaxies in the average region have less dust, and agree with the observed UV LF. The LF of the overdense region is substantially higher, and contains much brighter galaxies. The massive galaxies are bright in the infrared (IR) due to the dust thermal emission, with L_IR~ 3.7*10^12 Lsun at z = 6.3, while L_IR < 10^11 Lsun for the low-mass galaxies. Therefore, ALMA can probe massive galaxies in the overdense region up to z ~ 10 with a reasonable integration time. The UV spectral properties of disky galaxies depend significantly upon the viewing angle.The stellar and dust masses of the most massive galaxy in the overdense region are comparable to those of the sub-millimetre galaxy (SMG) found by Riechers et al. at z = 6.3, while the modelled SFR and the sub-millimetre flux fall slightly below the observed one. Statistical significance of these similarities and differences will only become clear with the upcoming ALMA observations.

preprint2015arXiv

The young radio lobe of 3C 84: inferred gas properties in the central 10 parsec

We analyse the environment of the supermassive black hole (SMBH) in the centre of a massive elliptical galaxy NGC 1275 in the Perseus cluster, hosting the radio source 3C 84. We focus on the young radio lobe observed inside the estimated Bondi accretion radius. We discuss the momentum balance between the jet associated with the lobe and the surrounding gas. The results are compared with the proper motion of the radio lobe obtained with the very long baseline interferometry. We find that under assumption of a high-density environment >~ 100 cm^-3), the jet power must be comparable to the Eddington luminosity --- this is clearly inconsistent with the current moderate activity of 3C 84, which indicates instead that the jet is expanding in a very low density region (<~1 cm^-3), along the rotation axis of the accretion flow. The power required for the jet to expand in the low-density environment is comparable to the past average jet power estimated from the X-ray observations. We estimate the classical Bondi accretion rate, assuming that (1) gas accretion is spherically symmetric, (2) accretion is associated with the jet environment, and (3) the medium surrounding the jet is representative of the properties of the dominant accreting gas. We find that Bondi accretion is inconsistent with the estimated jet power. This means that either accretion of the cold gas in the NGC 1275 is more efficient than that of the hot gas, or the jets are powered by the SMBH spin.

preprint2014arXiv

Secular Damping of Stellar Bars in Spinning Dark Matter Halos

We demonstrate that growth of stellar bars in spinning dark matter halos is heavily suppressed in the secular phase of evolution, using numerical simulations of isolated galaxies. In a representative set of models, we show that for values of the cosmological spin parameter lambda > 0.03, bar growth (in strength and size) becomes increasingly quenched. Furthermore, slowdown of bar pattern speed weakens substantially with increasing `lambda,' until it ceases completely. The terminal structure of bars is affected as well, including extent and shape of their boxy/peanut bulges. The essence of this effect lies in the modified angular momentum exchange between the disk and the halo facilitated by the bar. For the first time we have demonstrated that a dark matter halo can emit and not purely absorb angular momentum. Although the halo as a whole is not found to emit, the net transfer of angular momentum from the disk to the halo is significantly reduced or completely eliminated. The paradigm shift implies that the accepted view that disks serve as sources of angular momentum and halos serve as sinks, must be revised. Halos with lambda > 0.03 are expected to form a substantial fraction, based on lognormal distribution of lambda. Dependence of secular bar evolution on halo spin, therefore, implies profound corollaries for the cosmological evolution of galactic disks.

preprint2014arXiv

The Gentle Growth of Galaxies at High Redshifts in Overdense Environments

We have explored prevailing modes of galaxy growth for redshifts z ~ 6-14, comparing substantially overdense and normal regions of the universe, using high-resolution zoom-in cosmological simulations. Such rare overdense regions have been projected to host high-z quasars. We demonstrate that galaxies in such environments grow predominantly by a smooth accretion from cosmological filaments which dominates the mass input from major, intermediate and minor mergers. We find that by z ~6, the accumulated galaxy mass fraction from mergers falls short by a factor of 10 of the cumulative accretion mass for galaxies in the overdense regions, and by a factor of 5 in the normal environments. Moreover, the rate of the stellar mass input from mergers also lies below that of an in-situ star formation (SF) rate. The fraction of stellar masses in galaxies contributed by mergers in overdense regions is ~12%, and ~33% in the normal regions, at these redshifts. Our median SF rates for ~few X 10^9 Mo galaxies agrees well with the recently estimated rates for z ~ 7 galaxies from Spitzer's SURF-UP survey. Finally, we find that the main difference between the normal and overdense regions lies in the amplified growth of massive galaxies in massive dark matter halos. This leads to the formation of >= 10^{10} Mo galaxies due to the ~100-fold increase in mass during the above time period. Such galaxies are basically absent in the normal regions at these redshifts.

preprint2013arXiv

Supermassive Black Hole Formation at High Redshifts via Direct Collapse: Physical Processes in the Early Stage

We use numerical simulations to explore whether direct collapse can lead to the formation of SMBH seeds at high-z. We follow the evolution of gas within DM halos of 2 x 10^8 Mo and 1 kpc. We adopt cosmological density profiles and j-distributions. Our goal is to understand how the collapsing flow overcomes the centrifugal barrier and whether it is subject to fragmentation. We find that the collapse leads either to a central runaway or to off-center fragmentation. A disk-like configuration is formed inside the centrifugal barrier. For more cuspy DM distribution, the gas collapses more and experiences a bar-like perturbation and a central runaway. We have followed this inflow down to ~10^{-4} pc. The flow remains isothermal and the specific angular momentum is efficiently transferred by gravitational torques in a cascade of nested bars. This cascade supports a self-similar, disk-like collapse. In the collapsing phase, virial supersonic turbulence develops and fragmentation is damped. For larger initial DM cores the timescales become longer. In models with more organized initial rotation, a torus forms and appears to be supported by turbulent motions. The evolution depends on the competition between two timescales, corresponding to the onset of the central runaway and off-center fragmentation. For less organized rotation, the torus is greatly weakened, the central accretion timescale is shortened, and off-center fragmentation is suppressed --- triggering the central runaway even in previously `stable' models. The resulting SMBH masses lie in the range 2 x 10^4 Mo - 2 x 10^6 Mo, much higher than for Population III remnants. We argue that the above upper limit appears to be more realistic mass. Corollaries of this model include a possible correlation between SMBH and DM halo masses, and similarity between the SMBH and halo mass functions, at time of formation.

preprint2012arXiv

Cosmological Evolution of Galaxies

I review the subject of the cosmological evolution of galaxies, including different aspects of growth in disk galaxies, by focussing on the angular momentum problem, mergers, and their by-products. I discuss the alternative to merger-driven growth -- cold accretion and related issues. In the follow-up, I review possible feedback mechanisms and their role in galaxy evolution. Special attention is paid to high-redshift galaxies and their properties. In the next step, I discuss a number of processes, gas- and stellar-dynamical, within the central kiloparsec of disk galaxies, and their effect on the larger spatial scales, as well as on the formation and fuelling of the seed black holes in galactic centres at high redshifts.

preprint2011arXiv

Galaxy Formation in Heavily Overdense Regions at z~10: the Prevalence of Disks in Massive Halos

Using a high-resolution cosmological numerical simulation, we have analyzed the evolution of galaxies at z~10 in a highly overdense region of the universe. These objects could represent the high redshift galaxies recently observed by the Hubble's WFC3, and be as well possible precursors of QSOs at z~6-7. To overcome the sampling and resolution problems in cosmological simulations, we have used the Constrained Realizations method. Our main result for z~10 shows the region of 3.5h^{-1}Mpc radius in comoving coordinates completely dominated by disk galaxies in the total mass range of >=10^9h^{-1}Mo. We have verified that the gaseous and stellar disks we identify are robust morphological features, capable of surviving the ongoing merger process at these redshifts. Below this mass range, we find a sharp decline in the disk fraction to negligible numbers. At this redshift, the disks appear to be gas-rich and the dark matter halos baryon-rich, by a factor of ~2-3 above the average fraction of baryons in the universe. The prevalence of disk galaxies in the high density peaks during the epoch of reionization is contrary to the morphology-density trend observed at low redshifts.

preprint2011arXiv

The Dark Side of QSO Formation at High Redshifts

Observed high-redshift QSOs, at z~6, may reside in massive dark matter (DM) halos of more than 10^{12} Msun and are thus expected to be surrounded by overdense regions. In a series of 10 constrained simulations, we have tested the environment of such QSOs. Comparing the computed overdensities with respect to the unconstrained simulations of regions empty of QSOs, assuming there is no bias between the DM and baryon distributions, and invoking an observationally-constrained duty-cycle for Lyman Break Galaxies, we have obtained the galaxy count number for the QSO environment. We find that a clear discrepancy exists between the computed and observed galaxy counts in the Kim et al. (2009) samples. Our simulations predict that on average eight z~6 galaxies per QSO field should have been observed, while Kim et al. detect on average four galaxies per QSO field compared to an average of three galaxies in a control sample (GOODS fields). While we cannot rule out a small number statistics for the observed fields to high confidence, the discrepancy suggests that galaxy formation in the QSO neighborhood proceeds differently than in the field. We also find that QSO halos are the most massive of the simulated volume at z~6 but this is no longer true at z~3. This implies that QSO halos, even in the case they are the most massive ones at high redshifts, do not evolve into most massive galaxy clusters at z=0.

preprint2010arXiv

Dark Matter Halos and Evolution of Bars in Disk Galaxies: Varying Gas Fraction and Gas Spatial Resolution

We conduct numerical experiments by evolving gaseous/stellar disks embedded in live dark matter halos aiming at quantifying the effect of gas spatial resolution and gas content on the bar evolution. Model sequences have been constructed using different resolution, and gas fraction has been varied along each sequence within fgas=0%-50%, keeping the disk and halo properties unchanged. We find that the spatial resolution becomes important with an increase in `fgas'. For the higher resolution model sequences, we observe a bimodal behavior in the bar evolution with respect to the gas fraction, especially during the secular phase of this evolution. The switch from the gas-poor to gas-rich behavior is abrupt and depends on the resolution used. The diverging evolution has been observed in nearly all basic parameters characterizing bars, such as the bar strength, central mass concentration, vertical buckling amplitude, size, etc. We find that the presence of the gas component severely limits the bar growth and affects its pattern speed evolution. Gas-poor models display rapidly decelerating bars, while gas-rich models exhibit bars with constant or even slowly accelerating tumbling. The gas-rich models have bar corotation (CR) radii within the disk at all times, in constrast with gas-poor and purely stellar disks. The CR-to-bar size ratio is less than 2 for gas rich-models. We have confirmed that the disk angular momentum within the CR remains unchanged in the gas-poor models, as long as the CR stays within the disk, but experiences a sharp drop before leveling off in the gas-rich models. Finally, we discuss a number of observed correlations between various parameters of simulated bars, e.g., bar sizes and gas fractions, bar strength and buckling amplitude, bar strength and its size, etc.

preprint2010arXiv

Dissecting Galaxy Formation: II. Comparing Substructure in Pure Dark Matter and Baryonic Models

We compare the substructure evolution in pure dark matter (DM) halos with those in the presence of baryons (PDM and BDM). The prime halos have been analyzed by Romano-Diaz et al (2009). Models have been evolved from identical initial conditions using Constrained Realizations, including star formation and feedback. A comprehensive catalog of subhalos has been compiled and properties of subhalos analyzed in the mass range of 10^8 Mo - 10^11 Mo. We find that subhalo mass functions are consistent with a single power law, M_sbh^{alpha}, but detect a nonnegligible shift between these functions, alpha -0.86 for the PDM, and -0.98 for the BDM. Overall, alpha const. in time with variations of +-15%. Second, we find that the radial mass distribution of subhalos can be approximated by a power law, R^{gamma} with a steepening around the radius of a maximal circular velocity, Rvmax, in the prime halos. Gamma ~-1.5 for the PDM and -1 for the BDM, inside Rvmax, and is steeper outside. We detect little spatial bias between the subhalo populations and the DM of the main halos. The subhalo population exhibits much less triaxiality with baryons, in tandem with the prime halo. Finally, we find that, counter-intuitively, the BDM population is depleted at a faster rate than the PDM one within the central 30kpc of the prime. Although the baryons provide a substantial glue to the subhalos, the main halos exhibit the same trend. This assures a more efficient tidal disruption of the BDM subhalos. This effect can be reversed for a more efficient feedback from stellar evolution and supermassive black holes, which will expel baryons from the center and decrease the concentration of the prime halo. We compare our results with via Lactea and Aquarius simulations and other published results.

preprint2009arXiv

Angular Momentum Transfer and Lack of Fragmentation in Self-Gravitating Accretion Flows

Rapid inflows associated with early galaxy formation lead to the accumulation of self-gravitating gas in the centers of proto-galaxies. Such gas accumulations are prone to non-axisymmetric instabilities, as in the well-known Maclaurin sequence of rotating ellipsoids, which are accompanied by a catastrophic loss of angular momentum (J). Self-gravitating gas is also intuitively associated with star formation. However, recent simulations of the infall process display highly turbulent continuous flows. We propose that J-transfer, which enables the inflow, also suppresses fragmentation. Inefficient J loss by the gas leads to decay of turbulence, triggering global instabilities and renewed turbulence driving. Flow regulated in this way is stable against fragmentation, whilst staying close to the instability threshold for bar formation -- thick self-gravitating disks are prone to global instabilities before they become unstable locally. On smaller scales, the fraction of gravitationally unstable matter swept up by shocks in such a flow is a small and decreasing function of the Mach number. We conclude counterintuitively that gas able to cool down to a small fraction of its virial temperature will not fragment as it collapses. This provides a venue for supermassive black holes to form via direct infall, without the intermediary stage of forming a star cluster. Some black holes could have formed or grown in massive halos at low redshifts. Thus the fragmentation is intimately related to J redistribution within the system: it is less dependent on the molecular and metal cooling but is conditioned by the ability of the flow to develop virial, supersonic turbulence.

preprint2009arXiv

Dark Matter Halos and Evolution of Bars in Disk Galaxies: Collisionless Models Revisited

We construct and evolve families of steady-state models of stellar disks embedded in live DM halos, in order to study the dynamical and secular phases of bar evolution. These models are tested against those published in the literature in order to extend them and include the gaseous component in the follow up paper. We are interested in the angular momentum (J) redistribution in the disk-halo system. We confirm the previous results and quantify for the first time the dual role that the DM halos play in the bar evolution: more centrally concentrated halos dilute the dynamical processes, such as spontaneous bar instability and vertical buckling instability, and slowdown the J transfer, while facilitating it in the secular phase. Within the corotation radius (Rcr), the disk J remains nearly constant, as long as Rcr stays within the disk -- a sign that the lost J to the outer disk and the halo is being compensated by an influx of fresh J due to the outward motion of Rcr. This is feasible as long as the bar slowdown dominates the loss of J inside Rcr. We find that in some models the bar pattern speed stalls for prolonged time periods when Rcr is located outside the disk. This phenomenon appears concurrent with the near absence of J transfer between the disk and the halo. Furthermore, we confirm that stellar bars generally display the corotation to bar size ratios in the range of ~1-1.4, but only between the times of the first buckling and Rcr leaving the disk. The corotation-to-disk size ratio emerges as an important dynamic discriminator between various stages of barred disk evolution. Finally, we analyze a number of correlations between the basic parameters of a barred disk and a halo, some already reported in the literature and some new.

preprint2008arXiv

Anatomy of the Bar Instability in Cuspy Dark Matter Halos

We examine the bar instability in galactic models with an exponential disk and a cuspy dark matter (DM) halo with a Navarro-Frenk-White (NFW) cosmological density profile. We construct equilibrium models from a 3-integral composite distribution function that are subject to the bar instability. We generate a sequence of models with a range of mass resolution from 1.8K to 18M particles in the disk and 10K to 100M particles in the halo along with a multi-mass model with an effective resolution of ~10^10 particles. We describe how mass resolution affects the bar instability, including its linear growth phase, the buckling instability, pattern speed decay through the resonant transfer of angular momentum to the DM halo, and the possible destruction of the halo cusp. Our higher resolution simulations show a converging spectrum of discrete resonance interactions between the bar and DM halo orbits. As the pattern speed decays, orbital resonances sweep through most of the DM halo phase space and widely distribute angular momentum among the halo particles. The halo does not develop a flat density core and preserves the cusp, except in the region dominated by gravitational softening. The formation of the bar increases the central stellar density and the DM is compressed adiabatically increasing the halo central density by 1.7X. Overall, the evolution of the bar displays a convergent behavior for halo particle numbers between 1M and 10M particles, when comparing bar growth, pattern speed evolution, the DM halo density profile and a nonlinear analysis of the orbital resonances. Higher resolution simulations clearly illustrate the importance of discrete resonances in transporting the angular momentum from the bar to the halo.

preprint2008arXiv

Bars in Cuspy Dark Halos

We examine the bar instability in models with an exponential disk and a cuspy NFW-like dark matter (DM) halo inspired by cosmological simulations. Bar evolution is studied as a function of numerical resolution in a sequence of models spanning 10K to 100M DM particles - including a multi-mass model with an effective resolution of 10G. The goal is to find convergence in dynamical behaviour. We characterize the bar growth, the buckling instability, pattern speed decay through resonant transfer of angular momentum, and possible destruction of the DM halo cusp. Overall, most characteristics converge in behaviour in detail for halos containing more than 10M particles. Notably, the formation of the bar does not destroy the density cusp in this case. These higher resolution simulations clearly illustrate the importance of discrete resonances in transporting angular momentum from the bar to the halo.

preprint2007arXiv

Gas Feedback on Stellar Bar Evolution

We analyze evolution of live disk-halo systems in the presence of various gas fractions, f_gas less than 8% in the disk. We addressed the issue of angular momentum (J) transfer from the gas to the bar and its effect on the bar evolution. We find that the weakening of the bar, reported in the literature, is not related to the J-exchange with the gas, but is caused by the vertical buckling instability in the gas-poor disks and by a steep heating of a stellar velocity dispersion by the central mass concentration (CMC) in the gas-rich disks. The gas has a profound effect on the onset of the buckling -- larger f_gas brings it forth due to the more massive CMCs. The former process leads to the well-known formation of the peanut-shaped bulges, while the latter results in the formation of progressively more elliptical bulges, for larger f_gas. The subsequent (secular) evolution of the bar differs -- the gas-poor models exhibit a growing bar while gas-rich models show a declining bar whose vertical swelling is driven by a secular resonance heating. The border line between the gas-poor and -rich models lies at f_gas ~ 3% in our models, but is model-dependent and will be affected by additional processes, like star formation and feedback from stellar evolution. The overall effect of the gas on the evolution of the bar is not in a direct J transfer to the stars, but in the loss of J by the gas and its influx to the center that increases the CMC. The more massive CMC damps the vertical buckling instability and depopulates orbits responsible for the appearance of peanut-shaped bulges. The action of resonant and non-resonant processes in gas-poor and gas-rich disks leads to a converging evolution in the vertical extent of the bar and its stellar dispersion velocities, and to a diverging evolution in the bulge properties.

preprint2006arXiv

The AGN Obscuring Torus -- End of the "Doughnut" Paradigm?

Unified schemes of active galactic nuclei (AGN) require an obscuring dusty torus around the central engine. The compact sizes (only a few pc) determined in recent high-resolution observations require that the obscuring matter be clumpy and located inside the region where the black-hole gravity dominates over the galactic bulge. This location is in line with the scenario depicting the torus as the region of the clumpy wind coming off the accretion disk in which the clouds are dusty and optically thick. We study here the outflow scenario within the framework of hydromagnetic disk winds, incorporating the cloud properties determined from detailed modeling of the IR emission from clumpy tori. We find that torus clouds were likely detected in recent water maser observations of NGC 3079. In the wind scenario, the AGN main dynamic channel for release of accreted mass seems to be switching at low luminosities from torus outflow to radio jets. The torus disappears when the bolometric luminosity decreases below about \E{42} erg/sec because the accretion onto the central black hole can no longer sustain the required cloud outflow rate. This disappearance seems to have been observed in both LINERs and radio galaxies. With further luminosity decrease, suppression of cloud outflow spreads radially inward from the disk's dusty, molecular region into its atomic, ionized zone, resulting in disappearance of the broad emission line region at lower luminosities, yet to be determined.

preprint2001arXiv

Nested and Single Bars in Seyfert and Non-Seyfert Galaxies

We analyze the observed properties of nested and single stellar bar systems in disk galaxies. The 112 galaxies in our sample comprise the largest matched Seyfert vs. non-Seyfert sample of nearby galaxies with complete near-infrared or optical imaging sensitive to lengthscales ranging from tens of pc to tens of kpc. We find that a significant fraction of the sample galaxies, 17% +/- 4%, has more than one bar, and that 28% +/- 5% of barred galaxies have nested bars. The bar fractions appear to be stable according to reasonable changes in our adopted bar criteria. For these nested bars, we detect a clear division in length between the large-scale (primary) bars and small-scale (secondary) bars, both in absolute and normalized (to the size of the galaxy) length. We argue that this bimodal distribution can be understood within the framework of the inner Lindblad resonances (ILRs), which are located where the gravitational potential of the innermost galaxy switches effectively from 3D to 2D. While primary bars are found to correlate with the host galaxy sizes, no such correlation is observed for the secondary bars. Moreover, we find that secondary bars differ morphologically from single bars. Overall, our matched Seyfert and non-Seyfert samples show a statistically significant excess of bars among the Seyfert galaxies at practically all lengthscales. We confirm our previous result that Seyfert galaxies always show a preponderance of "thick" bars compared to the bars in non-Seyfert galaxies. Finally, no correlation is observed between the presence of a bar and that of companion galaxies, even relatively bright ones.

preprint1999arXiv

A subarcsecond resolution near-infrared study of Seyfert and `normal' galaxies: II. Morphology

We present a detailed study of the bar fraction in the CfA sample of Seyfert galaxies, and in a carefully selected control sample of non-active galaxies, to investigate the relation between the presence of bars and of nuclear activity. To avoid the problems related to bar classification in the RC3, e.g., subjectivity, low resolution and contamination by dust, we have developed an objective bar classification method, which we conservatively apply to our new sub-arcsecond resolution near-infrared imaging data set (Peletier et al. 1999). We are able to use stringent criteria based on radial profiles of ellipticity and major axis position angle to determine the presence of a bar and its axial ratio. Concentrating on non-interacting galaxies in our sample for which morphological information can be obtained, we find that Seyfert hosts are barred more often (79% +/- 7.5%) than the non-active galaxies in our control sample (59% +/- 9%), a result which is at the 2.5 sigma significance level. The fraction of non-axisymmetric hosts becomes even larger when interacting galaxies are taken into account. We discuss the implications of this result for the fueling of central activity by large-scale bars. This paper improves on previous work by means of imaging at higher spatial resolution and by the use of a set of stringent criteria for bar presence, and confirms that the use of NIR is superior to optical imaging for detection of bars in disk galaxies.

preprint1996arXiv

Dynamics of Inner Galactic Disks: The Striking Case of M100

We investigate gas dynamics in the presence of a double inner Lindblad resonance within a barred disk galaxy. Using an example of a prominent spiral, M100, we reproduce the basic central morphology, including four dominant regions of star formation corresponding to the compression maxima in the gas. These active star forming sites delineate an inner boundary (so-called nuclear ring) of a rather broad oval detected in the near infrared. We find that inclusion of self-gravitational effects in the gas is necessary in order to understand its behavior in the vicinity of the resonances and its subsequent evolution. The self-gravity of the gas is also crucial to estimate the effect of a massive nuclear ring on periodic orbits in the stellar bar.