Dynamical Evolution of Galaxies in Clusters

Tidal forces acting on galaxies in clusters lead to a strong dynamical evolution. In order to quantify the amount of evolution, I run self-consistent N-body simulations of disk galaxies for a variety of models in the hierarchically forming clusters. The tidal field along the galactic orbits is extracted from the simulations of cluster formation in the Omega_0=1; Omega_0=0.4; and Omega_0=0.4, Omega_Lambda=0.6 cosmological scenarios. For large spiral galaxies with the rotation speed of 250 km/s, tidal interactions truncate massive dark matter halos at 30 +- 6 kpc, and thicken stellar disks by a factor 2 to 3, increasing Toomre's parameter to Q > 2 and halting star formation. Low density galaxies, such as the dwarf spheroidals with the circular velocity of 20 km/s and the extended low surface brightness galaxies with the scale length of 10-15 kpc, are completely disrupted by tidal shocks. Their debris contribute to the diffuse intracluster light. The tidal effects are significant not only in the core but throughout the cluster and can be parametrized by the critical tidal density. The tidally-induced evolution results in the transformation of the infalling spirals into S0 galaxies and in the depletion of the LSB population. In the low Omega_0 cosmological models, clusters form earlier and produce stronger evolution of galaxies.Comment: accepted to Ap

Electron self-trapping on a nano-circle

We study the self-trapping of quasiparticles (electrons, holes, excitons, etc) in a molecular chain with the structure of a ring, taking into account the electron-phonon interaction and the radial and tangential deformations of the chain. A discrete system of equations is obtained and solved numerically. The analytical solutions for the wave function of a quasiparticle and for the molecule displacements that determine the distortion of the ring, are also obtained and solved in the continuum approximation. The numerical solutions of the system of discrete nonlinear equations reveals several regimes of quasiparticle localisation in the chain which depend on the values of the parameters of the system. It is shown that the transversal deformation of the chain favours the formation of a soliton.Comment: 43 pages 9 figure

Universal scaling of the elliptic flow data at RHIC

Recent PHOBOS measurements of the excitation function for the pseudo-rapidity dependence of elliptic flow in Au+Au collisions at RHIC, have posed a significant theoretical challenge. Here we show that these differential measurements, as well as the RHIC measurements on transverse momentum satisfy a universal scaling relation predicted by the Buda-Lund model, based on exact solutions of perfect fluid hydrodynamics. We also show that recently found transverse kinetic energy scaling of the elliptic flow is a special case of this universal scaling.Comment: 4 pages, 3 figures, 1 tabl

Constraints from Gravitational Recoil on the Growth of Supermassive Black Holes at High Redshift

Recent studies have shown that during their coalescence, binary supermassive black holes (SMBHs) experience a gravitational recoil with velocities of 100 km/s < v(kick) < 600 km/s. These velocities exceed the escape velocity v(esc) from typical dark matter (DM) halos at high-redshift (z>6), and therefore put constraints on scenarios in which early SMBHs grow at the centers of DM halos. Here we quantify these constraints for the most distant known SMBHs, with inferred masses in excess of 10^9 M(sun), powering the bright quasars discovered in the Sloan Digital Sky Survey at z>6. We assume that these SMBHs grew via a combination of accretion and mergers between pre-existing seed BHs in individual progenitor halos, and that mergers between progenitors with v(esc) < v(kick) disrupt the BH growth process. Our results suggest that under these assumptions, the z=6 SMBHs had a phase during which gained mass significantly more rapidly than under an Eddington-limited exponential growth rate.Comment: submitted to ApJ Letters, 5 emulateapj pages with 1 figur

The non-Gaussian tail of cosmic-shear statistics

Due to gravitational instability, an initially Gaussian density field develops non-Gaussian features as the Universe evolves. The most prominent non-Gaussian features are massive haloes, visible as clusters of galaxies. The distortion of high-redshift galaxy images due to the tidal gravitational field of the large-scale matter distribution, called cosmic shear, can be used to investigate the statistical properties of the LSS. In particular, non-Gaussian properties of the LSS will lead to a non-Gaussian distribution of cosmic-shear statistics. The aperture mass ($M_{\rm ap}$) statistics, recently introduced as a measure for cosmic shear, is particularly well suited for measuring these non-Gaussian properties. In this paper we calculate the highly non-Gaussian tail of the aperture mass probability distribution, assuming Press-Schechter theory for the halo abundance and the `universal' density profile of haloes as obtained from numerical simulations. We find that for values of $M_{\rm ap}$ much larger than its dispersion, this probability distribution is closely approximated by an exponential, rather than a Gaussian. We determine the amplitude and shape of this exponential for various cosmological models and aperture sizes, and show that wide-field imaging surveys can be used to distinguish between some of the currently most popular cosmogonies. Our study here is complementary to earlier cosmic-shear investigations which focussed more on two-point statistical properties.Comment: 9 pages, 5 figures, submitted to MNRA

Critical collapse and the primordial black hole initial mass function

It has normally been assumed that primordial black holes (PBHs) always form with mass approximately equal to the mass contained within the horizon at that time. Recent work studying the application of critical phenomena in gravitational collapse to PBH formation has shown that in fact, at a fixed time, PBHs with a range of masses are formed. When calculating the PBH initial mass function it is usually assumed that all PBHs form at the same horizon mass. It is not clear, however, that it is consistent to consider the spread in the mass of PBHs formed at a single horizon mass, whilst neglecting the range of horizon masses at which PBHs can form. We use the excursion set formalism to compute the PBH initial mass function, allowing for PBH formation at a range of horizon masses, for two forms of the density perturbation spectrum. First we examine power-law spectra with $n>1$, where PBHs form on small scales. We find that, in the limit where the number of PBHs formed is small enough to satisfy the observational constraints on their initial abundance, the mass function approaches that found by Niemeyer and Jedamzik under the assumption that all PBHs form at a single horizon mass. Second, we consider a flat perturbation spectrum with a spike at a scale corresponding to horizon mass $\sim 0.5 M_{\odot}$, and compare the resulting PBH mass function with that of the MACHOs (MAssive Compact Halo Objects) detected by microlensing observations. The predicted mass spectrum appears significantly wider than the steeply-falling spectrum found observationally.Comment: 8 pages RevTeX file with ten figures incorporated (uses RevTeX and epsf). Minor changes to dicussion onl

Massive perturbers and the efficient merger of binary massive black holes

We show that dynamical relaxation in the aftermath of a galactic merger and the ensuing formation and decay of a binary massive black hole (MBH), are dominated by massive perturbers (MPs) such as giant molecular clouds or clusters. MPs accelerate relaxation by orders of magnitude relative to 2-body stellar relaxation alone, and efficiently scatter stars into the binary MBH's orbit. The 3-body star-binary MBH interactions shrink the binary MBH to the point where energy losses from the emission of gravitational waves (GW) lead to rapid coalescence. We model this process based on observed and simulated MP distributions and take into account the decreased efficiency of the star-binary MBH interaction due to acceleration in the galactic potential. We show that mergers of gas-rich galactic nuclei lead to binary MBH coalescence well within the Hubble time. Moreover, lower-mass binary MBHs (<10^8 Msun) require only a few percent of the typical gas mass in a post-merger nucleus to coalesce in a Hubble time. The fate of a binary MBH in a gas poor galactic merger is less certain, although massive stellar structures (e.g. clusters, stellar rings) could likewise lead to efficient coalescence. These coalescence events are observable by their strong GW emission. MPs thus increase the cosmic rate of such GW events, lead to a higher mass deficit in the merged galactic core and suppress the formation of triple MBH systems and the resulting ejection of MBHs into intergalactic space.Comment: 14 pages, 4 figures, 3 tables. More detailed explanations and changes in structure. Section on hypervelocity stars moved to another paper (in preparation). Results and conclusions unchanged. Accepted to Ap

Cosmological Effects of Powerful AGN Outbursts in Galaxy Clusters: Insights from an XMM-Newton Observation of MS0735+7421

We report on the results of an analysis of XMM-Newton observations of MS0735+7421, the galaxy cluster which hosts the most energetic AGN outburst currently known. The previous Chandra image shows twin giant X-ray cavities (~200 kpc diameter) filled with radio emission and surrounded by a weak shock front. XMM data are consistent with these findings. The total energy in cavities and shock (~6 \times 10^{61} erg) is enough to quench the cooling flow and, since most of the energy is deposited outside the cooling region (~100 kpc), to heat the gas within 1 Mpc by ~1/4 keV per particle. The cluster exhibits an upward departure (factor ~2) from the mean L-T relation. The boost in emissivity produced by the ICM compression in the bright shells due to the cavity expansion may contribute to explain the high luminosity and high central gas mass fraction that we measure. The scaled temperature and metallicity profiles are in general agreement with those observed in relaxed clusters. Also, the quantities we measure are consistent with the observed M-T relation. We conclude that violent outbursts such as the one in MS0735+7421 do not cause dramatic instantaneous departures from cluster scaling relations (other than the L-T relation). However, if they are relatively common they may play a role in creating the global cluster properties.Comment: 69 pages, 30 figures, accepted for publication in ApJ Main Journa

Survival near ice sheet margins for some, but not all, North American trees

Temperate species experienced dramatic range reductions during the Last Glacial Maximum, yet refugial populations from which modern populations are descended have never been precisely located. Climate-based models identify only broad areas of potential habitat, traditional phylogeographic studies provide poor spatial resolution, and pollen records for temperate forest communities are difficult to interpret and do not provide species-level taxonomic resolution. Here we harness signals of range expansion from large genomic datasets, using a simulation-based framework to infer the precise latitude and longitude of glacial refugia in two widespread, codistributed hickories (Carya spp.) and to quantify uncertainty in these estimates. We show that one species likely expanded from close to ice sheet margins near the site of a previously described macrofossil for the genus, highlighting support for the controversial notion of northern microrefugia. In contrast, the expansion origin inferred for the second species is compatible with classic hypotheses of distant displacement into southern refugia. Our statistically rigorous, powerful approach demonstrates how refugia can be located from genomic data with high precision and accuracy, addressing fundamental questions about long-term responses to changing climates and providing statistical insight into longstanding questions that have previously been addressed primarily qualitatively.Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/148529/1/Bemmels2019.pdf85Description of Bemmels2019.pdf : Main articl

On the Cosmological Evolution of the Luminosity Function and the Accretion Rate of Quasars

We consider a class of models for the redshift evolution (between 0\lsim z \lsim 4) of the observed optical and X-ray quasar luminosity functions (LFs), with the following assumptions: (i) the mass-function of dark matter halos follows the Press-Schechter theory, (ii) the black hole (BH) mass scales linearly with the halo mass, (iii) quasars have a constant universal lifetime, and (iv) a thin accretion disk provides the optical luminosity of quasars, while the X-ray/optical flux ratio is calibrated from a sample of observed quasars. The mass accretion rate $\dot{M}$ onto quasar BHs is a free parameter of the models, that we constrain using the observed LFs. The accretion rate $\dot M$ inferred from either the optical or X-ray data under these assumptions generally decreases as a function of cosmic time from $z \simeq 4$ to $z \simeq 0$. We find that a comparable accretion rate is inferred from the X-ray and optical LF only if the X-ray/optical flux ratio decreases with BH mass. Near $z\simeq 0$, $\dot M$ drops to substantially sub-Eddington values at which advection-dominated accretion flows (ADAFs) exist. Such a decline of $\dot M$, possibly followed by a transition to radiatively inefficient ADAFs, could explain both the absence of bright quasars in the local universe and the faintness of accreting BHs at the centers of nearby galaxies. We argue that a decline of the accretion rate of the quasar population is indeed expected in cosmological structure formation models.Comment: Latex, 23 pages, 9 figures, accepted for publication in Ap