BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Thomas J. Cox

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

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-phastro-ph.HEastro-ph.SR

Catalog footprint

What is connected

31works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

31 published item(s)

preprint2014arXiv

Simulating multiple merger pathways to the central kinematics of early-type galaxies

Two-dimensional integral field surveys such as ATLAS^3D are producing rich observational data sets yielding insights into galaxy formation. These new kinematic observations have highlighted the need to understand the evolutionary mechanisms leading to a spectrum of fast-rotators and slow-rotators in early-type galaxies. We address the formation of slow and fast rotators through a series of controlled, comprehensive hydrodynamical simulations sampling idealized galaxy merger scenarios constructed from model spiral galaxies. Idealized and controlled simulations of this sort complement the more 'realistic' cosmological simulations by isolating and analyzing the effects of specific parameters, as we do in this paper. We recreate minor and major binary mergers, binary merger trees with multiple progenitors, and multiple sequential mergers. Within each of these categories of formation history, we correlate progenitor gas fraction, mass ratio, orbital pericenter, orbital ellipticity, and spin with remnant kinematic properties. We create kinematic profiles of these 95 simulations comparable to ATLAS^3D data. By constructing remnant profiles of the projected specific angular momentum (lambda_R = <R|V|> / <sqrt(V^2+sigma^2)>, triaxiality, and measuring the incidences of kinematic twists and kinematically decoupled cores, we distinguish between varying formation scenarios. We find that binary mergers nearly always form fast rotators. Slow rotators can be formed from zero initial angular momentum configurations and gas-poor mergers, but are not as round as the ATLAS^3D galaxies. Remnants of binary merger trees are triaxial slow rotators. Sequential mergers form round slow rotators that most resemble the ATLAS^3D rotators.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Modeling Mid-Infrared Diagnostics of Obscured Quasars and Starbursts

We analyze the link between active galactic nuclei (AGN) and mid-infrared flux using dust radiative transfer calculations of starbursts realized in hydrodynamical simulations. Focusing on the effects of galaxy dust, we evaluate diagnostics commonly used to disentangle AGN and star formation in ultraluminous infrared galaxies (ULIRGs). We examine these quantities as a function of time, viewing angle, dust model, AGN spectrum, and AGN strength in merger simulations representing two possible extremes of the ULIRG population: one is a typical gas-rich merger at z ~ 0, and the other is characteristic of extremely obscured starbursts at z ~ 2 to 4. This highly obscured burst begins star-formation-dominated with significant PAH emission, and ends with a ~10^9 yr period of red near-IR colors. At coalescence, when the AGN is most luminous, dust obscures the near-infrared AGN signature, reduces the relative emission from polycyclic aromatic hydrocarbons (PAHs), and enhances the 9.7 micron absorption by silicate grains. Although generally consistent with previous interpretations, our results imply none of these indicators can unambiguously estimate the AGN luminosity fraction in all cases. Motivated by the simulations, we show that a combination of the extinction feature at 9.7 micron, the PAH strength, and a near-infrared slope can simultaneously constrain the AGN fraction and dust grain distribution for a wide range of obscuration. We find that this indicator, accessible to the James Webb Space Telescope, may estimate the AGN power as tightly as the hard X-ray flux alone, thereby providing a valuable future cross-check and constraint for large samples of distant ULIRGs.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Star Formation in Galaxy Mergers with Realistic Models of Stellar Feedback & the Interstellar Medium

We use simulations with realistic models for stellar feedback to study galaxy mergers. These high resolution (1 pc) simulations follow formation and destruction of individual GMCs and star clusters. The final starburst is dominated by in situ star formation, fueled by gas which flows inwards due to global torques. The resulting high gas density results in rapid star formation. The gas is self gravitating, and forms massive (~10^10 M_sun) GMCs and subsequent super-starclusters (masses up to 10^8 M_sun). However, in contrast to some recent simulations, the bulk of new stars which eventually form the central bulge are not born in superclusters which then sink to the center of the galaxy, because feedback efficiently disperses GMCs after they turn several percent of their mass into stars. Most of the mass that reaches the nucleus does so in the form of gas. The Kennicutt-Schmidt law emerges naturally as a consequence of feedback balancing gravitational collapse, independent of the small-scale star formation microphysics. The same mechanisms that drive this relation in isolated galaxies, in particular radiation pressure from IR photons, extend over seven decades in SFR to regulate star formation in the most extreme starbursts (densities >10^4 M_sun/pc^2). Feedback also drives super-winds with large mass loss rates; but a significant fraction of the wind material falls back onto the disks at later times, leading to higher post-starburst SFRs in the presence of stellar feedback. Strong AGN feedback is required to explain sharp cutoffs in star formation rate. We compare the predicted relic structure, mass profile, morphology, and efficiency of disk survival to simulations which do not explicitly resolve GMCs or feedback. Global galaxy properties are similar, but sub-galactic properties and star formation rates can differ significantly.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Spatial Extent and Distribution of Star Formation in 3D-HST Mergers at z~1.5

We present an analysis of the spatial distribution of star formation in a sample of 60 visually identified galaxy merger candidates at z>1. Our sample, drawn from the 3D-HST survey, is flux-limited and was selected to have high star formation rates based on fits of their broad-band, low spatial resolution spectral energy distributions. It includes plausible pre-merger (close pairs) and post-merger (single objects with tidal features) systems, with total stellar masses and star formation rates derived from multi-wavelength photometry. Here we use near-infrared slitless spectra from 3D-HST which produce Halpha or [OIII] emission line maps as proxies for star-formation maps. This provides a first comprehensive high-resolution, empirical picture of where star formation occurred in galaxy mergers at the epoch of peak cosmic star formation rate. We find that detectable star formation can occur in one or both galaxy centres, or in tidal tails. The most common case (58%) is that star formation is largely concentrated in a single, compact region, coincident with the centre of (one of) the merger components. No correlations between star formation morphology and redshift, total stellar mass, or star formation rate are found. A restricted set of hydrodynamical merger simulations between similarly massive and gas-rich objects implies that star formation should be detectable in both merger components, when the gas fractions of the individual components are the same. This suggests that z~1.5 mergers typically occur between galaxies whose gas fractions, masses, and/or star formation rates are distinctly different from one another.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On Galaxies and Homology

The definition of homology for single-component galaxies is clear, but for multi-component (luminous and dark matter) galaxies there is some ambiguity. We attempt to clarify the situation by carefully separating the different concepts of homology that have been used to date. We argue that the most useful definition is that a set of galaxies is homologous if they are the same in all respects up to a set of three dimensional scaling constants which may differ from one galaxy to the next. Noting that we are free to choose the dimensional constants, we find that a set of hydrodynamic simulated galaxy merger remnants is significantly closer to homologous when the dimensional length constant is taken to be the radius containing equal amounts of dark and baryonic matter rather than the usual observationally motivated choice of the baryonic half-mass radius. Once the correct dimensional scaling constants are used, the stellar velocity dispersion anisotropy is essentially the sole source of the variation in the kinematic structure of these simulated merger remnants. In order to facilitate the use of these scaling constants to analyse observed galaxies, we calculated the relationship between our preferred dimensional scaling constants and the typical observationally accessible quantities.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

K+A Galaxies as the Aftermath of Gas-Rich Mergers: Simulating the Evolution of Galaxies as Seen by Spectroscopic Surveys

Models of poststarburst (or "K+A") galaxies are constructed by combining fully three-dimensional hydrodynamic simulations of galaxy mergers with radiative transfer calculations of dust attenuation. Spectral line catalogs are generated automatically from moderate-resolution optical spectra calculated as a function of merger progress in each of a large suite of simulations. The mass, gas fraction, orbital parameters, and mass ratio of the merging galaxies are varied systematically, showing that the lifetime and properties of the K+A phase are strong functions of merger scenario. K+A durations are generally less than ~0.1-0.3 Gyr, significantly shorter than the commonly assumed 1 Gyr, which is obtained only in rare cases, owing to a wide variation in star formation histories resulting from different orbital and progenitor configurations. Combined with empirical merger rates, the model lifetimes predict rapidly-rising K+A fractions as a function of redshift that are consistent with results of large spectroscopic surveys, resolving tension between the observed K+A abundance and that predicted when one assumes the K+A duration is the lifetime of A stars (~1 Gyr). The effects of dust attenuation, viewing angle, and aperture bias on our models are analyzed. In some cases, the K+A features are longer-lived and more pronounced when AGN feedback removes dust from the center, uncovering the young stars formed during the burst. In this picture, the K+A phase begins during or shortly after the bright starburst/AGN phase in violent mergers, and thus offers a unique opportunity to study the effects of quasar and star formation feedback on the gas reservoir and evolution of the remnant. Analytic fitting formulae are provided for the estimates of K+A incidence as a function of merger scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Recoiling Black Holes in Merging Galaxies

Gravitational-wave (GW) recoil of merging supermassive black holes (SMBHs) may influence the co-evolution of SMBHs and their host galaxies. We examine this possibility using SPH/N-body simulations of gaseous galaxy mergers in which the merged BH receives a recoil kick. With our suite of over 200 merger simulations, we identify systematic trends in the behavior of recoiling BHs. Our main results are as follows. (1) While BHs kicked at nearly the central escape speed (vesc) are essentially "lost" to the galaxy, in gas rich mergers, BHs kicked with up to about 0.7 vesc may be confined to the central few kpc of the galaxy. (2) The inflow of cold gas during a gas-rich major merger may cause a rapid increase in central escape speed; in such cases recoil trajectories will depend on the timing of the BH merger relative to the change in vesc. (3) Recoil events generally reduce the lifetimes of bright active galactic nuclei (AGN) but may actually extend AGN lifetimes at lower luminosities. (4) Recoiling AGN may be observable via kinematic offsets (v > 500 km s^-1) or spatial offsets (R > 1 kpc) for lifetimes of up to about 10 - 100 Myr. (5) Rapidly-recoiling BHs may be up to about 5 times less massive than their stationary counterparts. These mass deficits lower the normalization of the M - sigma relation and contribute to both intrinsic and overall scatter. (6) Finally, the displacement of AGN feedback by a recoil event causes higher central star formation rates in the merger remnant, thereby extending the starburst phase of the merger and creating a denser, more massive stellar cusp.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Recoiling Black Holes in Merging Galaxies: Relationship to AGN Lifetimes and Merger Remnant Properties

Central supermassive black holes (SMBHs) are a ubiquitous feature of locally-observed galaxies, and ample evidence suggests that the growth of SMBHs and their host galaxies is closely linked. However, in the event of a merger, gravitational-wave (GW) recoil may displace a SMBH from its galactic center, or eject it entirely. To explore the consequences of this phenomenon, we use hydrodynamic simulations of gaseous galaxy mergers that include a range of BH recoil velocities. We have generated a suite of over 200 simulations with more than 60 merger models, enabling us to identify systematic trends in the behavior of recoiling BHs -- specifically (i) their dynamics, (ii) their observable signatures, and (iii) their effects on BH/galaxy co-evolution. (i) Recoiling BH trajectories depend heavily on the gas content of the host galaxy; maximal BH displacements from the center may vary by up to an order of magnitude between gas-rich and gas-poor mergers. In some cases, recoil trajectories also depend on the timing of the BH merger relative to the formation of the galaxy merger remnant. (ii) Recoiling BHs may be observable as offset active galactic nuclei (AGN) via either kinematic offsets (v > 800 km s^-1) or spatial offsets (R > 1 kpc) for lifetimes of about 1 - 100 Myr. In addition, recoil events affect the total AGN lifetime. GW recoil generally reduces the lifetimes of bright AGN, but may extend lower-luminosity AGN lifetimes. (iii) Rapidly-recoiling BHs may be up to about 5 times less massive than their stationary counterparts. These mass deficits lower the normalization of the M - sigma relation and contribute to both intrinsic and overall scatter. Furthermore, recoil events displace AGN feedback from the galactic center, which enhances central star formation rates. This results in longer starburst phases and higher central stellar densities in merger remnants.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The effects of a hot gaseous halo in galaxy major mergers

Cosmological hydrodynamical simulations as well as observations indicate that spiral galaxies are comprised of five different components: dark matter halo, stellar disc, stellar bulge, gaseous disc and gaseous halo. While the first four components have been extensively considered in numerical simulations of binary galaxy mergers, the effect of a hot gaseous halo has usually been neglected even though it can contain up to 80% of the total gas within the galaxy virial radius. We present a series of hydrodynamic simulations of major mergers of disc galaxies, that for the first time include a diffuse, rotating, hot gaseous halo. Through cooling and accretion, the hot halo can dissipate and refuel the cold gas disc before and after a merger. This cold gas can subsequently form stars, thus impacting the morphology and kinematics of the remnant. Simulations of isolated systems with total mass M~10^12Msun show a nearly constant star formation rate of ~5Msun/yr if the hot gaseous halo is included, while the star formation rate declines exponentially if it is neglected. We conduct a detailed study of the star formation efficiency during mergers and find that the presence of a hot gaseous halo reduces the starburst efficiency (e=0.5) compared to simulations without a hot halo (e=0.68). Moreover we find cases where the stellar mass of the merger remnant is lower than the sum of the stellar mass of the two progenitor galaxies when evolved in isolation. This suggests a revision to semi-analytic galaxy formation models which assume that a merger always leads to enhanced star formation. We show that adding the hot gas component has a significant effect on the kinematics and internal structure of the merger remnants, like an increased abundance of fast rotators and an r^(1/4) surface brightness profile at small scales.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The effects of a hot gaseous halo on disc thickening in galaxy minor mergers

We employ hydrodynamical simulations to study the effects of dissipational gas physics on the vertical heating and thickening of disc galaxies during minor mergers. For the first time we present a suite of simulations that includes a diffuse, rotating, cooling, hot gaseous halo, as predicted by cosmological hydrodynamical simulations as well as models of galaxy formation. We study the effect of this new gaseous component on the vertical structure of a Milky Way-like stellar disc during 1:10 and 1:5 mergers. For 1:10 mergers we find no increased final thin disc scale height compared to the isolated simulation, leading to the conclusion that thin discs can be present even after a 1:10 merger if a reasonable amount of hot gas is present. The reason for this is the accretion of new cold gas, leading to the formation of a massive new thin stellar disc that dominates the surface brightness profile. In a previous study, in which we included only cold gas in the disk, we showed that the presence of cold gas decreased the thickening by a minor merger relative to the no-gas case. Here, we show that the evolution of the scale height in the presences of a cooling hot halo is dominated by the formation of the new stellar disc. In this scenario, the thick disc is the old stellar disc that has been thickened in a minor merger at z>1, while the thin disc is the new stellar disc that reforms after this merger. In addition, we study the evolution of the scale height during a 1:5 merger and find that a thin disc can be present even after this merger, provided enough hot gas is available. The final scale height in our simulations depends on the mass of the hot gaseous halo, the efficiency of the winds and the merger mass ratio. We find post-merger values in the range 0.5<z0<1.0 kpc in good agreement with observational constraints by local galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Kennicutt-Schmidt Star Formation Relation at z~2

Recent observations of excited CO emission lines from z~2 disc galaxies have shed light on the Kennicutt-Schmidt relation at high-z via observed SFR-CO (J=2-1) and (J=3-2) relations. Here, we describe a novel methodology for utilising these observations of high-excitation CO to derive the underlying Schmidt (SFR-rho^N) relationship. To do this requires an understanding of the potential effects of differential CO excitation with SFR. If the most heavily star-forming galaxies have a larger fraction of their gas in highly excited CO states than the lower SFR galaxies, then the observed molecular SFR-CO^alpha index, alpha, will be less than the underlying (volumetric) Schmidt index, N. Utilising a combination of SPH models of galaxy evolution and molecular line radiative transfer, we present the first calculations of CO excitation in z~2 disc galaxies with the aim of developing a mapping between various observed SFR-CO relationships and the underlying Schmidt relation. We find that even in relatively luminous z~2 discs, differential excitation does indeed exist, resulting in alpha < N for highly excited CO lines. This means that an observed (e.g.) SFR-CO (J=3-2) relation does not map linearly to SFR-H2 relation. We utilise our model results to provide a mapping from alpha to N for the range of Schmidt indices N=1-2. By comparing to recent observational surveys, we find that the observed SFR-CO (J=2-1) and SFR-CO (J=3-2) relations suggest that an underlying SFR rho^1.5 relation describes z~2 disc galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Metallicity Evolution of Interacting Galaxies

Nuclear inflows of metal-poor interstellar gas triggered by galaxy interactions can account for the systematically lower central oxygen abundances observed in local interacting galaxies. Here, we investigate the metallicity evolution of a large set of simulations of colliding galaxies. Our models include cooling, star formation, feedback, and a new stochastic method for tracking the mass recycled back to the interstellar medium from stellar winds and supernovae. We study the influence of merger-induced inflows, enrichment, gas consumption, and galactic winds in determining the nuclear metallicity. The central metallicity is primarily a competition between the inflow of low-metallicity gas and enrichment from star formation. An average depression in the nuclear metallicity of ~0.07 is found for gas-poor disk-disk interactions. Gas-rich disk-disk interactions, on the other hand, typically have an enhancement in the central metallicity that is positively correlated with the gas content. The simulations fare reasonably well when compared to the observed mass-metallicity and separation-metallicity relationships, but further study is warranted.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Physical Model for z~2 Dust Obscured Galaxies

We present a physical model for the origin of z~2 Dust-Obscured Galaxies (DOGs), a class of high-redshift ULIRGs selected at 24 micron which are particularly optically faint (24/R>1000). By combining N-body/SPH simulations of high redshift galaxy evolution with 3D polychromatic dust radiative transfer models, we find that luminous DOGs (with F24 > 0.3 mJy at z~2 are well-modeled as extreme gas-rich mergers in massive (~5x10^12-10^13 Msun) halos, with elevated star formation rates (~500-1000 Msun/yr) and/or significant AGN growth (Mdot > 0.5 Msun/yr), whereas less luminous DOGs are more diverse in nature. At final coalescence, merger-driven DOGs transition from being starburst dominated to AGN dominated, evolving from a "bump" to a power-law shaped mid-IR (IRAC) spectral energy distribution (SED). After the DOG phase, the galaxy settles back to exhibiting a "bump" SED with bluer colors and lower star formation rates. While canonically power-law galaxies are associated with being AGN-dominated, we find that the power-law mid-IR SED can owe both to direct AGN contribution, as well as to a heavily dust obscured stellar bump at times that the galaxy is starburst dominated. Thus power-law galaxies can be either starburst or AGN dominated. Less luminous DOGs can be well-represented either by mergers, or by massive ($M_{\rm baryon} ~5x10^11 Msun) secularly evolving gas-rich disc galaxies (with SFR > 50 Msun/yr). By utilising similar models as those employed in the SMG formation study of Narayanan et al. (2010), we investigate the connection between DOGs and SMGs. We find that the most heavily star-forming merger driven DOGs can be selected as Submillimetre Galaxies (SMGs), while both merger-driven and secularly evolving DOGs typically satisfy the BzK selection criteria.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Discovery of a large-scale galaxy filament near a candidate intergalactic X-ray absorption system

We present an analysis of the large-scale galaxy distribution around two possible warm-hot intergalactic medium (WHIM) absorption systems reported along the Markarian 421 sightline. Using the Sloan Digital Sky Survey, we find a prominent galaxy filament at the redshift of the z=0.027 X-ray absorption line system. The filament exhibits a width of 3.2 Mpc and length of at least 20 Mpc, comparable to the size of WHIM filaments seen in cosmological simulations. No individual galaxies fall within 350 projected kpc so it is unlikely that the absorption is associated with gas in a galaxy halo or outflow. Another, lower-significance X-ray absorption system was reported in the same Chandra spectrum at z=0.011, but the large-scale structure in its vicinity is far weaker and may be a spurious alignment. By searching for similar galaxy structures in 140 random smoothed SDSS fields, we estimate a ~5-10% probability of the z=0.027 absorber-filament alignment occurring by chance. If these two systems are indeed physically associated, this would represent the first known coincidence between large-scale galaxy structure and a blind X-ray WHIM detection.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On Sizes, Kinematics, M/L Gradients, and Light Profiles of Massive Compact Galaxies at z~2

We present a detailed analysis of the structure and resolved stellar populations of simulated merger remnants, and compare them to observations of compact quiescent galaxies at z ~ 2. We find that major merging is a viable mechanism to produce systems of ~ 10^11 Msun and ~ 1 kpc size, provided the gas fraction at the time of final coalescence is high (~ 40%), and provided that the progenitors are compact star-forming galaxies, as expected at high redshift. Their integrated spectral energy distributions and velocity dispersions are in good agreement with the observations, and their position in the (v_{maj}/sigma, ellipticity) diagram traces the upper envelope of the distribution of lower redshift early-type galaxies. The simulated merger remnants show time- and sightline-dependent M/L ratio gradients that result from a superposition of radially dependent stellar age, stellar metallicity, and extinction. The median ratio of effective radius in rest-frame V-band light to that in mass surface density is ~ 2 during the quiescent remnant phase. This is typically expressed by a negative color gradient (i.e., red core), which we expect to correlate with the integrated color of the system. Finally, the simulations differ from the observations in their surface brightness profile shape. The simulated remnants are typically best fit by high (n >> 4) Sersic indices, whereas observed quiescent galaxies at z ~ 2 tend to be less cuspy (median n ~ 2.3). Limiting early star formation in the progenitors may be required to prevent the simulated merger remnants from having extended wings.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Orbital structure of merger remnants: Trends with gas fraction in 1:1 mergers

Since the violent relaxation in hierarchical merging is incomplete, elliptical galaxies retain a wealth of information about their formation pathways in their present-day orbital structure. A variety of evidence indicates that gas-rich major mergers play an important role in the formation of elliptical galaxies. We simulate 1:1 disk mergers at seven different initial gas fractions ranging from 0 to 40%, using the TreeSPH code Gadget-2. We classify the stellar orbits in each remnant and construct radial profiles of the orbital content, intrinsic shape, and orientation. The dissipationless remnants are typically prolate-triaxial, dominated by box orbits within r_c ~ 1.5Reff, and by tube orbits in their outer parts. As the gas fraction increases, the box orbits within r_c are increasingly replaced by a population of short axis tubes (z-tubes) with near zero net rotation, and the remnants become progressively more oblate and round. The long axis tube (x-tube) orbits are highly streaming and relatively insensitive to the gas fraction, implying that their angular momentum is retained from the dynamically cold initial conditions. Outside r_c, the orbital structure is essentially unchanged by the gas. The 15-20% gas remnants often display disk-like kinematically distinct cores (KDCs). These remnants show an interesting resemblance, in both their velocity maps and intrinsic orbital structure, to the KDC galaxy NGC4365 (van den Bosch et al. 2008). At 30-40% gas, the remnants are rapidly rotating, with sharp embedded disks on ~ 1Reff scales. We predict a characteristic, physically intuitive orbital structure for 1:1 disk merger remnants, with a distinct transition between 1 and 3Reff that will be readily observable with combined data from the 2D kinematics surveys SAURON and SMEAGOL.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Recoiling Black Holes in Merging Galaxies: Relationship to AGN Lifetimes, Starbursts, and the M-sigma Relation

Gravitational-wave (GW) recoil of merging supermassive black holes (SMBHs) may influence the co-evolution of SMBHs and their host galaxies. We examine this possibility using SPH/N-body simulations of gaseous galaxy mergers in which the merged BH receives a recoil kick. This enables us to follow recoiling BHs in self-consistent, evolving merger remnants. In contrast to recent studies on similar topics, we conduct a large parameter study, generating a suite of over 200 simulations with more than 60 merger models and a range of recoil velocities (vk). Our main results are as follows. (1) BHs kicked at nearly the central escape speed (vesc) may oscillate on large orbits for up to a Hubble time, but in gas-rich mergers, BHs kicked with up to ~ 0.7 vesc may be confined to the central few kpc of the galaxy, owing to gas drag and steep central potentials. (2) vesc in gas-rich mergers may increase rapidly during final coalescence, in which case trajectories may depend on the timing of the BH merger relative to the formation of the potential well. (3) Recoil events generally reduce the lifetimes of bright active galactic nuclei (AGN), but may actually extend AGN lifetimes at lower luminosities. (4) Kinematically-offset AGN (v > 800 km s^-1) may be observable for up to ~ 10 Myr either immediately after the recoil or during pericentric passages through a gas-rich remnant. (5) Spatially-offset AGN (R > 1 kpc) generally have low luminosities and lifetimes of ~ 1 - 100 Myr. (6) Rapidly-recoiling BHs may be up to ~ 5 times less massive than their stationary counterparts. This lowers the normalization of the M-sigma relation and contributes to both intrinsic and overall scatter. (7) Finally, the displacement of AGN feedback after a recoil event enhances central star formation rates, thereby extending the starburst phase of the merger and creating a denser stellar cusp. [Abridged.]

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

SDSS J1254+0846: A Binary Quasar Caught in the Act of Merging

We present the first luminous, spatially resolved binary quasar that clearly inhabits an ongoing galaxy merger. SDSS J125455.09+084653.9 and SDSS J125454.87+084652.1 (SDSS J1254+0846 hereafter) are two luminous z=0.44 radio quiet quasars, with a radial velocity difference of just 215 km/s, separated on the sky by 21 kpc in a disturbed host galaxy merger showing obvious tidal tails. The pair was targeted as part of a complete sample of binary quasar candidates with small transverse separations drawn from SDSS DR6 photometry. We present follow-up optical imaging which shows broad, symmetrical tidal arm features spanning some 75 kpc at the quasars' redshift. Numerical modeling suggests that the system consists of two massive disk galaxies prograde to their mutual orbit, caught during the first passage of an active merger. This demonstrates rapid black hole growth during the early stages of a merger between galaxies with pre-existing bulges. Neither of the two luminous nuclei show significant instrinsic absorption by gas or dust in our optical or X-ray observations, illustrating that not all merging quasars will be in an obscured, ultraluminous phase. We find that the Eddington ratio for the fainter component B is rather normal, while for the A component L/LEdd is quite (>3sigma) high compared to quasars of similar luminosity and redshift, possibly evidence for strong merger-triggered accretion. More such mergers should be identifiable at higher redshifts using binary quasars as tracers.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Discriminating Between the Physical Processes that Drive Spheroid Size Evolution

Massive galaxies at high-z have smaller effective radii than those today, but similar central densities. Their size growth therefore relates primarily to the evolving abundance of low-density material. Various models have been proposed to explain this evolution, which have different implications for galaxy, star, and BH formation. We compile observations of spheroid properties as a function of redshift and use them to test proposed models. Evolution in progenitor gas-richness with redshift gives rise to initial formation of smaller spheroids at high-z. These systems can then evolve in apparent or physical size via several channels: (1) equal-density 'dry' mergers, (2) later major or minor 'dry' mergers with less-dense galaxies, (3) adiabatic expansion, (4) evolution in stellar populations & mass-to-light-ratio gradients, (5) age-dependent bias in stellar mass estimators, (6) observational fitting/selection effects. If any one of these is tuned to explain observed size evolution, they make distinct predictions for evolution in other galaxy properties. Only model (2) is consistent with observations as a dominant effect. It is the only model which allows for an increase in M_BH/M_bulge with redshift. Still, the amount of merging needed is larger than that observed or predicted. We therefore compare cosmologically motivated simulations, in which all these effects occur, & show they are consistent with all the observational constraints. Effect (2), which builds up an extended low-density envelope, does dominate the evolution, but effects 1,3,4, & 6 each contribute ~20% to the size evolution (a net factor ~2). This naturally also predicts evolution in M_BH-sigma similar to that observed.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Recovering Stellar Population Properties and Redshifts from Broad-Band Photometry of Simulated Galaxies: Lessons for SED Modeling

We present a detailed analysis of our ability to determine stellar masses, ages, reddening and extinction values, and star formation rates of high-redshift galaxies by modeling broad-band SEDs with stellar population synthesis. In order to do so, we computed synthetic optical-to-NIR SEDs for model galaxies taken from hydrodynamical merger simulations placed at redshifts 1.5 < z < 3. Viewed under different angles and during different evolutionary phases, the simulations represent a wide variety of galaxy types (disks, mergers, spheroids). We show that simulated galaxies span a wide range in SEDs and color, comparable to these of observed galaxies. In all star-forming phases, dust attenuation has a large effect on colors, SEDs, and fluxes. The broad-band SEDs were then fed to a standard SED modeling procedure and resulting stellar population parameters were compared to their true values. Disk galaxies generally show a decent median correspondence between the true and estimated mass and age, but suffer from large uncertainties. During the merger itself, we find larger offsets (e.g., log M_recovered - log M_true = -0.13^{+0.10}_{-0.14}). E(B-V) values are generally recovered well, but the estimated total visual absorption Av is consistently too low, increasingly so for larger optical depths. Since the largest optical depths occur during the phases of most intense star formation, it is for the highest SFRs that we find the largest underestimates. The masses, ages, E(B-V), Av, and SFR of merger remnants (spheroids) are very well reproduced. We discuss possible biases in SED modeling results caused by mismatch between the true and template star formation history, dust distribution, metallicity variations and AGN contribution.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Resonant stripping as the origin of dwarf spheroidal galaxies

Dwarf spheroidal galaxies are the most dark matter dominated systems in the nearby Universe and their origin is one of the outstanding puzzles of how galaxies form. Dwarf spheroidals are poor in gas and stars, making them unusually faint, and those known as ultra-faint dwarfs have by far the lowest measured stellar content of any galaxy. Previous theories require that dwarf spheroidals orbit near giant galaxies like the Milky Way, but some dwarfs have been observed in the outskirts of the Local Group. Here we report simulations of encounters between dwarf disk galaxies and somewhat larger objects. We find that the encounters excite a process, which we term ``resonant stripping'', that can transform them into dwarf spheroidals. This effect is distinct from other mechanisms proposed to form dwarf spheroidals, including mergers, galaxy-galaxy harassment, or tidal and ram pressure stripping, because it is driven by gravitational resonances. It may account for the observed properties of dwarf spheroidals in the Local Group, including their morphologies and kinematics. Resonant stripping predicts that dwarf spheroidals should form through encounters, leaving detectable long stellar streams and tails.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Effects of Gas on Morphological Transformation in Mergers: Implications for Bulge and Disk Demographics

Transformation of disks into spheroids via mergers is a well-accepted element of galaxy formation models. However, recent simulations have shown that bulge formation is suppressed in increasingly gas-rich mergers. We investigate the global implications of these results in a cosmological framework, using independent approaches: empirical halo-occupation models (where galaxies are populated in halos according to observations) and semi-analytic models. In both, ignoring the effects of gas in mergers leads to the over-production of spheroids: low and intermediate-mass galaxies are predicted to be bulge-dominated (B/T~0.5 at <10^10 M_sun), with almost no bulgeless systems), even if they have avoided major mergers. Including the different physical behavior of gas in mergers immediately leads to a dramatic change: bulge formation is suppressed in low-mass galaxies, observed to be gas-rich (giving B/T~0.1 at <10^10 M_sun, with a number of bulgeless galaxies in good agreement with observations). Simulations and analytic models which neglect the similarity-breaking behavior of gas have difficulty reproducing the strong observed morphology-mass relation. However, the observed dependence of gas fractions on mass, combined with suppression of bulge formation in gas-rich mergers, naturally leads to the observed trends. Discrepancies between observations and models that ignore the role of gas increase with redshift; in models that treat gas properly, galaxies are predicted to be less bulge-dominated at high redshifts, in agreement with the observations. We discuss implications for the global bulge mass density and future observational tests.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Formation of High Redshift Submillimeter Galaxies

We describe a model for the formation of \zsim 2 Submillimeter Galaxies (SMGs) which simultaneously accounts for both average and bright SMGs while providing a reasonable match to their mean observed spectral energy distributions (SEDs). By coupling hydrodynamic simulations of galaxy mergers with the high resolution 3D polychromatic radiative transfer code Sunrise, we find that a mass sequence of merger models which use observational constraints as physical input naturally yield objects which exhibit black hole, bulge, and H2 gas masses similar to those observed in SMGs. The dominant drivers behind the 850 micron flux are the masses of the merging galaxies and the stellar birthcloud covering fraction. The most luminous (S850 ~ 15 mJy) sources are recovered by ~10^13 Msun 1:1 major mergers with a birthcloud covering fraction close to unity, whereas more average SMGs ~5-7 mJy) may be formed in lower mass halos ~5x10^12 Msun. These models demonstrate the need for high spatial resolution hydrodynamic and radiative transfer simulations in matching both the most luminous sources as well as the full SEDs of SMGs. While these models suggest a natural formation mechanism for SMGs, they do not attempt to match cosmological statistics of galaxy populations; future efforts along this line will help ascertain the robustness of these models.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The imprint of dissipation on the shapes of merger remnant LOSVDs

The properties of elliptical galaxies are broadly consistent with simulated remnants of gas-rich mergers between spirals, motivating more detailed studies of the imprint of this formation mechanism on the remnant distribution function. Gas has a strong impact on the non-Gaussian shapes of the line-of-sight velocity distributions (LOSVDs) of the merger remnant, owing to the embedded disk that forms out of the gas that retains its angular momentum during the merger, and the strong central mass concentration from the gas that falls to the center. The deviations from Gaussianity are parametrized by the Gauss-Hermite moments h_3 and h_4, which are related to the skewness and kurtosis of the LOSVDs. We quantify the dependence of the (h_3,h_4)-v/sigma relations on the initial gas fraction of the progenitor disks in 1:1 mergers, using Gadget-2 simulations including star formation, radiative cooling, and feedback from supernovae and AGN. For gas fractions f_gas < ~15% the overall correlation between h_3 and v/sigma is weak, consisting of a flat negatively correlated component arising from edge-on viewing angles plus a steep positively correlated part from face-on projections. The spread in v/sigma values decreases toward high positive h_4, and there is a trend toward lower h_4 as the gas fraction increases from 0 to 15%. For f_gas > ~20% the (h_3,4)- v/sigma distributions look quite different - there is a tight negative h_3- v/sigma correlation, and a wide spread in v/sigma values at all h_4, in better agreement with observations. Re-mergers of the high-f_gas remnants (dry mergers) produce slowly rotating systems with nearly Gaussian LOSVDs. We explain all of these trends in terms of the underlying orbit structure of the remnants, as molded by their dissipative formation histories.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Star-Forming Molecular Gas in High Redshift Submillimeter Galaxies

We present a model for the CO molecular line emission from high redshift Submillimeter Galaxies (SMGs). By combining hydrodynamic simulations of gas rich galaxy mergers with the polychromatic radiative transfer code, Sunrise, and the 3D non-LTE molecular line radiative transfer code, Turtlebeach, we show that if SMGs are typically a transient phase of major mergers, their observed compact CO spatial extents, broad line widths, and high excitation conditions (CO SED) are naturally explained. In this sense, SMGs can be understood as scaled-up analogs to local ULIRGs. We utilize these models to investigate the usage of CO as an indicator of physical conditions. We find that care must be taken when applying standard techniques. The usage of CO line widths as a dynamical mass estimator from SMGs can possibly overestimate the true enclosed mass by a factor ~1.5-2. At the same time, assumptions of line ratios of unity from CO J=3-2 (and higher lying lines) to CO (J=1-0) will oftentimes lead to underestimates of the inferred gas mass. We provide tests for these models by outlining predictions for experiments which are imminently feasible with the current generation of bolometer arrays and radio-wave spectrometers.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Dissipation and Extra Light in Galactic Nuclei: IV. Evolution in the Scaling Relations of Spheroids

We develop a model for the origins and redshift evolution of spheroid scaling relations. We consider spheroid sizes, velocity dispersions, masses, profile shapes (Sersic indices), and black hole (BH) masses, and their related scalings. Our approach combines advantages of observational constraints in halo occupation models and hydrodynamic merger simulations. This allows us to separate the relative roles of dissipation, dry mergers, formation time, and progenitor evolution, and identify their effects on scalings at each redshift. Dissipation is the most important factor determining spheroid sizes and fundamental plane (FP) scalings, and can account for the FP tilt and differences between disk and spheroid scalings. Because disks at high-z have higher gas fractions, mergers are more gas-rich, yielding more compact spheroids. This predicts mass-dependent evolution in spheroid sizes, in agreement with observations. This relates to subtle evolution in the FP, important to studies that assume a fixed intrinsic FP. This also predicts mild evolution in BH-host correlations, towards larger BHs at higher z. Dry mergers are significant, but only for massive systems which form early: they form compact, but undergo dry mergers (consistent with observations) such that their sizes at later times are similar to spheroids of similar mass formed more recently. We model descendants of observed compact high-z spheroids: most will become cores of BCGs, with sizes, velocity dispersions, and BH masses consistent with observations, but we identify a fraction that might survive to z=0 intact.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

How Do Disks Survive Mergers?

We develop a physical model for how galactic disks survive and/or are destroyed in interactions. Based on dynamical arguments, we show gas primarily loses angular momentum to internal torques in a merger. Gas within some characteristic radius (a function of the orbital parameters, mass ratio, and gas fraction of the merging galaxies), will quickly lose angular momentum to the stars sharing the perturbed disk, fall to the center and be consumed in a starburst. A similar analysis predicts where violent relaxation of the stellar disks is efficient. Our model allows us to predict the stellar and gas content that will survive to re-form a disk in the remnant, versus being violently relaxed or contributing to a starburst. We test this in hydrodynamic simulations and find good agreement as a function of mass ratio, orbital parameters, and gas fraction, in simulations spanning a wide range in these properties and others, including different prescriptions for gas physics and feedback. In an immediate sense, the amount of disk that re-forms can be understood in terms of well-understood gravitational physics, independent of details of ISM gas physics or feedback. This allows us to explicitly quantify the requirements for such feedback to (indirectly) enable disk survival, by changing the pre-merger gas content and distribution. The efficiency of disk destruction is a strong function of gas content: we show how and why sufficiently gas-rich major mergers can, under general conditions, yield systems with small bulges (B/T<0.2). We provide prescriptions for inclusion of our results in semi-analytic models.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Photometric Properties of the Most Massive High-Redshift Galaxies

We calculate the observable properties of the most massive high-redshift galaxies in the hierarchical formation scenario where stellar spheroid and supermassive black hole growth are fueled by gas-rich mergers. Combining high-resolution hydrodynamical simulations of the hierarchical formation of a z~6 quasar, stellar population synthesis models, template AGN spectra, prescriptions for interstellar and intergalactic absorption, and the response of modern telescopes, the photometric evolution of galaxies destined to host z~6 quasars are modeled at redshifts z~4-14. These massive galaxies, with enormous stellar masses of M_star ~10^11.5-10^12 M_sun. and star formation rates of SFR~10^3-10^4 M_sun yr^-1 at z>~7, satisfy a variety of photometric selection criteria based on Lyman-break techniques including V-band dropouts at z>~5, i-band dropouts at z>~6, and z-band dropouts at z>~7. The observability of the most massive high-redshift galaxies is assessed and compared with a wide range of existing and future photometric surveys including SDSS, GOODS/HUDF, NOAO WDFS, UKIDSS, the IRAC Shallow Survey, Pan-STARRS, LSST, and SNAP. Massive stellar spheroids descended from z~6 quasars will likely be detected at z~4 by existing surveys, but owing to their low number densities the discovery of quasar progenitor galaxies at z>7 will likely require future surveys of large portions of the sky (>~0.5%) at wavelengths lambda>1 micron. The detection of rare, star-bursting, massive galaxies at redshifts z>~6 would provide support for the hierarchical formation of the earliest quasars and characterize the primitive star-formation histories of the most luminous elliptical galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

An Upper Limit to the Degree of Evolution Between Supermassive Black Holes and their Host Galaxies

We describe a model-independent integral constraint which defines an upper limit to the allowed degree of evolution in the ratio of black hole (BH) mass to host galaxy luminosity or mass, as a function of redshift. Essentially, if the BH/host ratio is excessive at redshift z, then it would imply that the total mass density in BHs above some M_BH(min) is larger at that redshift than at z=0, which is impossible. This argument requires no knowledge of host or BH properties, only a lower limit to the observed luminosity density in the brightest galaxies at some z. We calculate this upper limit from a variety of luminosity and mass functions in different bands from redshifts z=0-2. We show that it is consistent with passive evolution of spheroid populations (with a fixed M_BH/M_host relation) in all cases, and provides tighter constraints than have generally been obtained previously, ruling out at >6 sigma observational and theoretical estimates suggesting that M_BH/M_host was significantly larger at high redshifts than locally, although relatively weak (factor ~2 by z=2) evolution is still allowed. We discuss a variety of possible 'loopholes' or changes in the BH/host populations and correlations, and show that they typically lower the upper limits and strengthen our conclusions.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

The Relation Between Quasar and Merging Galaxy Luminosity Functions and the Merger-Induced Star Formation Rate of the Universe

Using a model for self-regulated growth of black holes (BHs) in mergers involving gas-rich galaxies, we study the relationship between quasars and the population of merging galaxies and predict the merger-induced star formation rate density of the Universe. Mergers drive nuclear gas inflows, fueling starbursts and 'buried quasars' until accretion feedback expels the gas, rendering a briefly visible optical quasar. Star formation is shut down and accretion declines, leaving a passively evolving remnant with properties typical of red, elliptical galaxies. Based on evolution of these events in our simulations, we demonstrate that the observed statistics of merger rates, luminosity functions (LFs) and mass functions, SFR distributions, specific SFRs, quasar and quasar host galaxy LFs, and elliptical/red galaxy LFs are self-consistent and follow from one another as predicted by the merger hypothesis. We use our simulations to de-convolve both quasar and merging galaxy LFs to determine the birthrate of black holes of a given final mass and merger rates as a function of stellar mass. We use this to predict the merging galaxy LF in several observed wavebands, color-magnitude relations, mass functions, absolute and specific SFR distributions and SFR density, and quasar host galaxy LFs, as a function of redshift from z=0-6. We invert this and predict e.g. quasar LFs from observed merger LFs or SFR distributions. Our results agree well with observations, but idealized models of quasar lightcurves are ruled out by comparison of merger and quasar observations at >99.9% confidence. Using only observations of quasars, we estimate the contribution of mergers to the SFR density of the Universe even to high redshifts z~4.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

The Evolution of the M-sigma Relation

(Abridged) We examine the evolution of the black hole mass - stellar velocity dispersion (M-sigma) relation over cosmic time using simulations of galaxy mergers that include feedback from supermassive black hole growth. We consider mergers of galaxies varying the properties of the progenitors to match those expected at redshifts z=0-6. We find that the slope of the resulting M-sigma relation is the same at all redshifts considered. For the same feedback efficiency that reproduces the observed amplitude of the M-sigma relation at z=0, there is a weak redshift-dependence to the normalization that results from an increasing velocity dispersion for a given galactic stellar mass. We develop a formalism to connect redshift evolution in the M-sigma relation to the scatter in the local relation at z=0. We show that the scatter in the local relation places severe constraints on the redshift evolution of both the normalization and slope of the M-sigma relation. Furthermore, we demonstrate that cosmic downsizing introduces a black hole mass-dependent dispersion in the M-sigma relation and that the skewness of the distribution about the locally observed M-sigma relation is sensitive to redshift evolution in the normalization and slope. In principle, these various diagnostics provide a method for differentiating between theories for producing the M-sigma relation. In agreement with existing constraints, our simulations imply that hierarchical structure formation should produce the relation with small intrinsic scatter.

0 citations0 reviewsNo code linkedOpen access