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

Kelly Holley-Bockelmann

Kelly Holley-Bockelmann appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COastro-ph.HEastro-ph.IMgr-qchep-phastro-phastro-ph.SRhep-thphysics.ed-phphysics.soc-ph

Catalog footprint

What is connected

38works
11topics
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

38 published item(s)

preprint2022arXiv

Advancing the Landscape of Multimessenger Science in the Next Decade

The last decade has brought about a profound transformation in multimessenger science. Ten years ago, facilities had been built or were under construction that would eventually discover the nature of objects in our universe could be detected through multiple messengers. Nonetheless, multimessenger science was hardly more than a dream. The rewards for our foresight were finally realized through IceCube's discovery of the diffuse astrophysical neutrino flux, the first observation of gravitational waves by LIGO, and the first joint detections in gravitational waves and photons and in neutrinos and photons. Today we live in the dawn of the multimessenger era. The successes of the multimessenger campaigns of the last decade have pushed multimessenger science to the forefront of priority science areas in both the particle physics and the astrophysics communities. Multimessenger science provides new methods of testing fundamental theories about the nature of matter and energy, particularly in conditions that are not reproducible on Earth. This white paper will present the science and facilities that will provide opportunities for the particle physics community renew its commitment and maintain its leadership in multimessenger science.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Eccentricity evolution of massive black hole binaries from formation to coalescence

Coalescing supermassive black hole binaries (BHBs) are expected to be the loudest sources of gravitational waves (GWs) in the Universe. Detection rates for ground or space-based detectors based on cosmological simulations and semi-analytic models are highly uncertain. A major difficulty stems from the necessity to model the BHB from the scale of the merger to that of inspiral. Of particular relevance to the GW merger timescale is the binary eccentricity. Here we present a self-consistent numerical study of the eccentricity of BHBs formed in massive gas-free mergers from the early stages of the merger to the hardening phase, followed by a semi-analytical model down to coalescence. We find that the early eccentricity of the unbound black hole pair is largely determined by the initial orbit. It systematically decreases during the dynamical friction phase. The eccentricity at binary formation is affected by stochasticity and noise owing to encounters with stars, but preserves a strong correlation with the initial orbital eccentricity. Binding of the black holes is a phase characterised by strong perturbations, and we present a quantitative definition of the time of binary formation. During hardening the eccentricity increases in minor mergers, unless the binary is approximately circular, but remains largely unchanged in major mergers, in agreement with predictions from semi-analytical models based on isotropic scattering experiments. Coalescence times due to hardening and GW emission in gas-poor non-rotating ellipticals are <~0.5 Gyr for the large initial eccentricities (0.5 < e < 0.9) typical of galaxy mergers in cosmological simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

PhD bridge programmes as engines for access, diversity and inclusion

The lack of diversity in physics and astronomy PhD programmes is well known but has not improved despite decades of efforts. PhD bridge programmes provide an asset-based model to help overcome the societal and disciplinary obstacles to improving access and inclusion for students from underrepresented groups and are beginning to show some success. We describe several well-known PhD bridge programmes in the United States and discuss lessons learned from their experiences. Many of these lessons can be extended more broadly to physics and astronomy PhD programmes to increase access, diversity and inclusion.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Membership lists for 431 open clusters in Gaia DR2 using extreme deconvolution gaussian mixture models

Open clusters are groups of stars that form at the same time, making them an ideal laboratory to test theories of star formation, stellar evolution, and dynamics in the Milky Way disk. However, the utility of an open cluster can be limited by the accuracy and completeness of its known members. Here, we employ a "top-down" technique, {\it extreme deconvolution gaussian mixture models} (XDGMM), to extract and evaluate known open clusters from Gaia DR2 by fitting the distribution of stellar parallax and proper motion along a line-of-sight. Extreme deconvolution techniques can recover the intrinsic distribution of astrometric quantities, accounting for the full covariance matrix of the errors; this allows open cluster members to be identified even when presented with relatively uncertain measurement data. To date, open cluster studies have only applied extreme deconvolution to specialized searches for individual systems. We use XDGMM to characterize the open clusters reported by Ahumada et al. 2007 and are able to recover 420 of the 426 open clusters therein (98.1\%). Our membership list contains the overwhelming majority ($>95\%$) of previously known cluster members. We also identify a new, significant, and relatively faint cluster member population and validate their membership status using Gaia eDR3. We report the fortuitous discovery of 11 new open cluster candidates within the lines of sight we analyzed. We present our technique, its advantages and challenges, as well as publish our membership lists and updated cluster parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Getting Ready for LISA: The Data, Support and Preparation Needed to Maximize US Participation in Space-Based Gravitational Wave Science

The NASA LISA Study Team was tasked to study how NASA might support US scientists to participate and maximize the science return from the Laser Interferometer Space Antenna (LISA) mission. LISA is gravitational wave observatory led by ESA with NASA as a junior partner, and is scheduled to launch in 2034. Among our findings: LISA science productivity is greatly enhanced by a full-featured US science center and an open access data model. As other major missions have demonstrated, a science center acts as both a locus and an amplifier of research innovation, data analysis, user support, user training and user interaction. In its most basic function, a US Science Center could facilitate entry into LISA science by hosting a Data Processing Center and a portal for the US community to access LISA data products. However, an enhanced LISA Science Center could: support one of the parallel independent processing pipelines required for data product validation; stimulate the high level of research on data analysis that LISA demands; support users unfamiliar with a novel observatory; facilitate astrophysics and fundamental research; provide an interface into the subtleties of the instrument to validate extraordinary discoveries; train new users; and expand the research community through guest investigator, postdoc and student programs. Establishing a US LISA Science Center well before launch can have a beneficial impact on the participation of the broader astronomical community by providing training, hosting topical workshops, disseminating mock catalogs, software pipelines, and documentation. Past experience indicates that successful science centers are established several years before launch; this early adoption model may be especially relevant for a pioneering mission like LISA.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Inferring Parameters of GW170502: The Loudest Intermediate-mass Black Hole Trigger in LIGO's O1/O2 data

Gravitational wave (GW) measurements provide the most robust constraints of the mass of astrophysical black holes. Using state-of-the-art GW signal models and a unique parameter estimation technique, we infer the source parameters of the loudest marginal trigger, GW170502, found by LIGO from 2015 to 2017. If this trigger is assumed to be a binary black hole merger, we find it corresponds to a total mass in the source frame of $157^{+55}_{-41}~\rm{M}_\odot$ at redshift $z=1.37^{+0.93}_{-0.64}$. The primary and secondary black hole masses are constrained to $94^{+44}_{-28}~\rm{M}_{\odot}$ and $62^{+30}_{-25}~\rm{M}_{\odot}$ respectively, with 90\% confidence. Across all signal models, we find $\gtrsim 70\%$ probability for the effective spin parameter $χ_\mathrm{eff}>0.1$. Furthermore, we find that the inclusion of higher-order modes in the analysis narrows the confidence region for the primary black hole mass by 10\%, however, the evidence for these modes in the data remains negligible. The techniques outlined in this study could lead to robust inference of the physical parameters for all intermediate-mass black hole binary candidates $(\gtrsim100~\mathrm{M}_\odot)$ in the current GW network.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The role of gravitational recoil in the assembly of massive black hole seeds

When two black holes merge, the asymmetric emission of gravitational waves provides an impulse to the merged system; this gravitational wave recoil velocity can be up to 4000 km s$^{-1}$, easily fast enough for the black hole to escape its host galaxy. We combine semi-analytic modeling with cosmological zoom-in simulations of a Milky Way-type galaxy to investigate the role of black hole spin and gravitational recoil in the epoch of massive black hole seeding. We sample four different spin distributions (random, aligned, anti-aligned, and zero spin), and compare the resulting merger rates, occupation fractions, and MBH-host relations with what is expected by excluding the effect of recoil. The inclusion of gravitational recoil and MBH spin in the assembly of MBH seeds can reduce the final $z=5$ MBH mass by up to an order of magnitude. The MBH occupation fraction, however, remains effectively unaltered due to episodes of black hole formation following a recoil event. While electromagnetic detections of these events are unlikely, LISA is ideally suited to detect gravitational wave signals from such events.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Inward Bound: The incredible journey of massive black holes as they pair and merge; I. The effect of mass ratio in flattened rotating galactic nuclei

Understanding how supermassive black holes (SMBHs) pair and merge helps to inform predictions of off-center, dual, and binary AGN, and provides key insights into how SMBHs grow and co-evolve with their galaxy hosts. As the loudest known gravitational wave source, binary SMBH mergers also hold centerstage for the Laser Interferometer Space Antenna (LISA), a joint ESA/NASA gravitational wave observatory set to launch in 2034. Here, we continue our work to characterize SMBH binary formation and evolution through increasingly more realistic high resolution direct $N$-body simulations, focusing on the effect of SMBH mass ratio, orientation, and eccentricity within a rotating and flattened stellar host. During the dynamical friction phase, we found a prolonged orbital decay for retrograde SMBHs and swift pairing timescales for prograde SMBHs compared to their counterparts in non-rotating models, an effect that becomes more pronounced for smaller mass ratios $M_{\rm sec}/M_{\rm prim} = q$. During this pairing phase, the eccentricity dramatically increases for retrograde configurations, but as the binary forms, the orbital plane flips so that it is almost perfectly prograde, which stifles the rapid eccentricity growth. In prograde configurations, SMBH binaries form and remain at comparatively low eccentricities. As in our prior work, we note that the center of mass of a prograde SMBH binary itself settles into an orbit about the center of the galaxy. Since even the initially retrograde binaries flip their orbital plane, we expect few binaries in rotating systems to reside at rest in the dynamic center of the host galaxy, though this effect is smaller as $q$ decreases.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A First Look at Galaxy Flyby Interactions. II. Do Flybys matter?

In the second paper of this series, we present results from cosmological simulations on the demographics of flyby interactions to gauge their potential impact on galaxy evolution. In a previous paper, we demonstrated that flybys -- an interaction where two independent halos inter-penetrate but detach at a later time and do not merge -- occur much more frequently than previously believed. In particular, we found that the frequency of flybys increases at low redshift and is comparable to or even greater than the frequency of mergers for halos $\gtrsim 10^{11} M_\odot/h$. In this paper, we classify flybys according to their orbits and the level of perturbation exacted on both the halos involved. We find that the majority of flybys penetrate deeper than $\sim R_{half}$ of the primary and have an initial relative speed $\sim 1.6\times V_{vir}$ of the primary. The typical flyby mass-ratio is $\sim 0.1$ at high $z$ for all halos, while at low $z$, massive primary halos undergo flybys with small secondary halos. We estimate the perturbation from the flyby on both the primary and the secondary and find that a typical flyby is mostly non-perturbative for the primary halo. However, since a massive primary experiences so many flybys at any given time, they are nearly continually a victim of a perturbative event. In particular, we find flybys that cause $\sim 1\%$ change in the binding energy of a primary halo occurs $\gtrsim 1 $ Gyr$^{-1}$ for halos $> 10^{10} M_\odot/h$ for $z \lesssim 4$. Secondary halos, on the other hand, are highly perturbed by the typical encounter, experiencing a change in binding energy of nearly order unity. Our results imply that flybys can drive a significant part of galaxy transformation at moderate to lower redshifts ($z \lesssim 4$). We touch on implications for observational surveys, mass-to-light ratios, and galaxy assembly bias.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Galaxy Rotation and Rapid Supermassive Black Hole Binary Coalescence

During a galaxy merger, the supermassive black hole (SMBH) in each galaxy is thought to sink to the center of the potential and form a supermassive black hole binary; this binary can eject stars via 3-body scattering, bringing the SMBHs ever closer. In a static spherical galaxy model, the binary stalls at a separation of about a parsec after ejecting all the stars in its loss cone -- this is the well-known final parsec problem. However it has been shown that SMBH binaries in non-spherical galactic nuclei harden at a nearly constant rate until reaching the gravitational wave regime. Here we use a suite of direct N-body simulations to follow SMBH binary evolution in both corotating and counterrotating flattened galaxy models. For N larger than 500K, we find that the evolution of the SMBH binary is convergent, and is independent of the particle number. Rotation in general increases the hardening rate of SMBH binaries even more effectively than galaxy geometry alone. SMBH binary hardening rates are similar for co- and counterrotating galaxies. In the corotating case, the center of mass of SMBH binary settles into an orbit that is in a corotation resonance with the background rotating model, and the coalescence time is roughly few hundred Myr faster than a non-rotating flattened model. We find that counterrotation drives SMBHs to coalesce on a nearly radial orbit promptly after forming a hard binary. We discuss the implications for gravitational wave astronomy, hypervelocity star production, and the effect on the structure of the host galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Massive Black Hole Science with eLISA

The evolving Laser Interferometer Space Antenna (eLISA) will revolutionize our understanding of the formation and evolution of massive black holes along cosmic history by probing massive black hole binaries in the $10^3-10^7$ solar mass range out to redshift $z\gtrsim 10$. High signal-to-noise ratio detections of $\sim 10-100$ binary coalescences per year will allow accurate measurements of the parameters of individual binaries (such as their masses, spins and luminosity distance), and a deep understanding of the underlying cosmic massive black hole parent population. This wealth of unprecedented information can lead to breakthroughs in many areas of physics, including astrophysics, cosmology and fundamental physics. We review the current status of the field, recent progress and future challenges.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Probing Galactic Structure with the Spatial Correlation Function of SEGUE G-dwarf Stars

We measure the two-point correlation function of G-dwarf stars within 1-3 kpc of the Sun in multiple lines-of-sight using the Schlesinger et al. G-dwarf sample from the SDSS SEGUE survey. The shapes of the correlation functions along individual SEGUE lines-of-sight depend sensitively on both the stellar-density gradients and the survey geometry. We fit smooth disk galaxy models to our SEGUE clustering measurements, and obtain strong constraints on the thin- and thick-disk components of the Milky Way. Specifically, we constrain the values of the thin- and thick-disk scale heights with 3% and 2% precision, respectively, and the values of the thin- and thick-disk scale lengths with 20% and 8% precision, respectively. Moreover, we find that a two-disk model is unable to fully explain our clustering measurements, which exhibit an excess of clustering at small scales (< 50 pc). This suggests the presence of small-scale substructure in the disk system of the Milky Way.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Eleventh and Twelfth Data Releases of the Sloan Digital Sky Survey: Final Data from SDSS-III

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Voronoi Tessellation and Non-parametric Halo Concentration

We present and test TesseRACt, a non-parametric technique for recovering the concentration of simulated dark matter halos using Voronoi tessellation. TesseRACt is tested on idealized N-body halos that are axisymmetric, triaxial, and contain substructure and compared to traditional least-squares fitting as well as two non-parametric techniques that assume spherical symmetry. TesseRACt recovers halo concentrations within 0.3% of the true value regardless of whether the halo is spherical, axisymmetric, or triaxial. Traditional fitting and non-parametric techniques that assume spherical symmetry can return concentrations that are systematically off by as much as 10% from the true value for non-spherical halos. TesseRACt also performs significantly better when there is substructure present outside $0.5R_{200}$. Given that cosmological halos are rarely spherical and often contain substructure, we discuss implications for studies of halo concentration in cosmological N-body simulations including how choice of technique for measuring concentration might bias scaling relations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Bar Formation from Galaxy Flybys

Recently, both simulations and observations have revealed that flybys - fast, one-time interactions between two galaxy halos - are surprisingly common, nearing/comparable to galaxy mergers. Since these are rapid, transient events with the closest approach well outside the galaxy disk, it is unclear if flybys can transform the galaxy in a lasting way. We conduct collisionless N-body simulations of three co-planer flyby interactions between pure-disk galaxies to take a first look at the effects flybys have on disk structure, with particular focus on stellar bar formation. We find that some flybys are capable of inciting a bar with bars forming in both galaxies during our 1:1 interaction and in the secondary during our 10:1 interaction. The bars formed have ellipticities >0.5, sizes on the order of the host disk's scale length, and persist to the end of our simulations, ~5 Gyr after pericenter. The ability of flybys to incite bar formation implies that many processes associated with secular bar evolution may be more closely tied with interactions than previously though.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Classification of Stellar Orbits in Axisymmetric Galaxies

It is known that two supermassive black holes (SMBHs) cannot merge in a spherical galaxy within a Hubble time; an emerging picture is that galaxy geometry, rotation, and large potential perturbations may usher the SMBH binary through the critical three-body scattering phase and ultimately drive the SMBH to coalesce. We explore the orbital content within an N-body model of a mildly- flattened, non-rotating, SMBH-embedded elliptical galaxy. When used as the foundation for a study on the SMBH binary coalescence, the black holes bypassed the binary stalling often seen within spherical galaxies and merged on Gyr timescales (Khan et al. 2013). Using both frequency-mapping and angular momentum criteria, we identify a wealth of resonant orbits in the axisymmetric model, including saucers, that are absent from an otherwise identical spherical system and that can potentially interact with the binary. We quantified the set of orbits that could be scattered by the SMBH binary, and found that the axisymmetric model contained nearly seven times the number of these potential loss cone orbits compared to our equivalent spherical model. In this flattened model, the mass of these orbits is roughly 3 times of that of the SMBH, which is consistent with what the SMBH binary needs to scatter to transition into the gravitational wave regime.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Early Growth in a Perturbed Universe: Exploring Dark Matter Halo Populations in 2LPT and ZA Simulations

We study the structure and evolution of dark matter halos from z = 300 to z = 6 for two cosmological N-body simulation initialization techniques. While the second order Lagrangian perturbation theory (2LPT) and the Zel'dovich approximation (ZA) both produce accurate present day halo mass functions, earlier collapse of dense regions in 2LPT can result in larger mass halos at high redshift. We explore the differences in dark matter halo mass and concentration due to initialization method through three 2LPT and three ZA initialized cosmological simulations. We find that 2LPT induces more rapid halo growth, resulting in more massive halos compared to ZA. This effect is most pronounced for high mass halos and at high redshift. Halo concentration is, on average, largely similar between 2LPT and ZA, but retains differences when viewed as a function of halo mass. For both mass and concentration, the difference between typical individual halos can be very large, highlighting the shortcomings of ZA-initialized simulations for high-z halo population studies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Effects of Inclination on Measuring Velocity Dispersion and Implications for Black Holes

The relation of central black hole mass and stellar spheroid velocity dispersion (the M-$σ$ relation) is one of the best-known and tightest correlations linking black holes and their host galaxies. There has been much scrutiny concerning the difficulty of obtaining accurate black hole measurements, and rightly so; however, it has been taken for granted that measurements of velocity dispersion are essentially straightforward. We examine five disk galaxies from cosmological SPH simulations and find that line-of-sight effects due to galaxy orientation can affect the measured $σ$ by 30%, and consequently black hole mass predictions by up to 1.0 dex. Face-on orientations correspond to systematically lower velocity dispersion measurements, while more edge-on orientations give higher velocity dispersions, due to contamination by disk stars when measuring line of sight quantities. We caution observers that the uncertainty of velocity dispersion measurements is at least 20 km/s, and can be much larger for moderate inclinations. This effect may account for some of the scatter in the locally measured M-$σ$ relation, particularly at the low-mass end. We provide a method for correcting observed $σ_{\rm los}$ values for inclination effects based on observable quantities.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Expansion Techniques for Collisionless Stellar Dynamical Simulations

We present GPU implementations of two fast force calculation methods, based on series expansions of the Poisson equation. One is the Self-Consistent Field (SCF) method, which is a Fourier-like expansion of the density field in some basis set; the other is the Multipole Expansion (MEX) method, which is a Taylor-like expansion of the Green's function. MEX, which has been advocated in the past, has not gained as much popularity as SCF. Both are particle-field method and optimized for collisionless galactic dynamics, but while SCF is a "pure" expansion, MEX is an expansion in just the angular part; it is thus capable of capturing radial structure easily, where SCF needs a large number of radial terms. We show that despite the expansion bias, these methods are more accurate than direct techniques for the same number of particles. The performance of our GPU code, which we call ETICS, is profiled and compared to a CPU implementation. On the tested GPU hardware, a full force calculation for one million particles took ~ 0.1 seconds (depending on expansion cutoff), making simulations with as many as $10^8$ particles fast on a comparatively small number of nodes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hypervelocity Star Candidates in the SEGUE G & K Dwarf Sample

We present 20 candidate hypervelocity stars from the Sloan Extension for Galactic Understanding and Exploration (SEGUE) G and K dwarf samples. Previous searches for hypervelocity stars have only focused on large radial velocities; in this study we also use proper motions to select the candidates. We determine the hypervelocity likelihood of each candidate by means of Monte Carlo simulations, considering the significant errors often associated with high proper motion stars. We find that nearly half of the candidates exceed their escape velocities with at least 98% probability. Every candidate also has less than a 25% chance of being a high-velocity fluke within the SEGUE sample. Based on orbits calculated using the observed six-dimensional positions and velocities, few, if any, of these candidates originate from the Galactic center. If these candidates are truly hypervelocity stars, they were not ejected by interactions with the Milky Way's supermassive black hole. This calls for a more serious examination of alternative hypervelocity-star ejection scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ultramassive Black Hole Coalescence

Although supermassive black holes (SMBHs) correlate well with their host galaxies, there is an emerging view that outliers exist. Henize 2-10, NGC 4889, and NGC1277 are examples of SMBHs at least an order of magnitude more massive than their host galaxy suggests. The dynamical effects of such ultramassive central black holes is unclear. Here, we perform direct N-body simulations of mergers of galactic nuclei where one black hole is ultramassive to study the evolution of the remnant and the black hole dynamics in this extreme regime. We find that the merger remnant is axisymmetric near the center, while near the large SMBH influence radius, the galaxy is triaxial. The SMBH separation shrinks rapidly due to dynamical friction, and quickly forms a binary black hole; if we scale our model to the most massive estimate for the NGC1277 black hole, for example, the timescale for the SMBH separation to shrink from nearly a kiloparsec to less than a parsec is roughly 10 Myr. By the time the SMBHs form a hard binary, gravitational wave emission dominates, and the black holes coalesce in a mere few Myr. Curiously, these extremely massive binaries appear to nearly bypass the 3-body scattering evolutionary phase. Our study suggests that in this extreme case, SMBH coalescence is governed by dynamical friction followed nearly directly by gravitational wave emission, resulting in an rapid and efficient SMBH coalescence timescale. We discuss the implications for gravitational wave event rates and hypervelocity star production.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A First Look at Galaxy Flyby Interactions: Characterizing the Frequency of Flybys in a Cosmological Context

Hierarchical structure formation theory is based on the notion that mergers drive galaxy evolution, so a considerable framework of semi-analytic models and N-body simulations has been constructed to calculate how mergers transform a growing galaxy. However, galaxy mergers are only one type of major dynamical interaction between halos -- another class of encounter, a close flyby, has been largely ignored. We analyze a 50 Mpc/h, $1024^3$ collisionless cosmological simulation and find that the number of close flyby interactions is comparable to, or even surpasses, the number of mergers for halo masses $\ga 10^{11}\,{h^{-1} M_\odot}$ at $z \la 2$. Halo flybys occur so frequently to high mass halos that they are continually perturbed, unable to reach a dynamical equilibrium. We also find tentative evidence that at high redshift, $z \ga 14$, flybys are as frequent as mergers. Our results suggest that close halo flybys can play an important role in the evolution of the earliest dark matter halos and their galaxies, and can still influence galaxy evolution at the present epoch.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Can a Satellite Galaxy Merger Explain the Active Past of the Galactic Center?

Observations of the Galactic Center (GC) have accumulated a multitude of "forensic" evidence indicating that several million years ago the center of the Milky Way galaxy was teaming with starforming and accretion-powered activity -- this paints a rather different picture from the GC as we understand it today. We examine a possibility that this epoch of activity could have been triggered by the infall of a satellite galaxy into the Milky Way which began at the redshift of 10 and ended few million years ago with a merger of the Galactic supermassive black hole with an intermediate mass black hole brought in by the inspiralling satellite.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Supermassive Black Hole Binary Evolution in Axisymmetric Galaxies: the final parsec problem is not a problem

During a galaxy merger, the supermassive black hole (SMBH) in each galaxy is thought to sink to the center of the potential and form a supermassive black hole binary; this binary can eject stars via 3-body scattering, bringing the SMBHs ever closer. In a static spherical galaxy model, the binary stalls at a separation of about a parsec after ejecting all the stars in its loss cone -- this is the well-known final parsec problem. Earlier work has shown that the centrophilic orbits in triaxial galaxy models are key in refilling the loss cone at a high enough rate to prevent the black holes from stalling. However, the evolution of binary SMBHs has never been explored in axisymmetric galaxies, so it is not clear if the final parsec problem persists in these systems. Here we use a suite of direct N-body simulations to follow SMBH binary evolution in galaxy models with a range of ellipticity. For the first time, we show that mere axisymmetry can solve the final parsec problem; we find the the SMBH evolution is independent of N for an axis ratio of c/a=0.8, and that the SMBH binary separation reaches the gravitational radiation regime for c/a=0.75.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A First Look at Galaxy Flyby Interactions: I. Characterizing the Frequency of Flybys in a Cosmological Context

Hierarchical structure formation theory is based on the notion that mergers drive galaxy evolution, so a considerable framework of semi-analytic models and N-body simulations has been constructed to calculate how mergers transform a growing galaxy. However, galaxy mergers are only one type of major dynamical interaction between halos -- another class of encounter, a close flyby, has been largely ignored. We use cosmological N-body simulations to reconstruct the entire dynamical interaction history of dark matter halos. We present a careful method of identifying and tracking a dark matter halo which resolves the typical classes of anomalies that occur in N-body data. This technique allows us to robustly follow halos and several hierarchical levels of subhalos as they grow, dissolve, merge, and flyby one another -- thereby constructing both a census of the dynamical interactions in a volume and an archive of the dynamical evolution of an individual halo. In addition to a census of mergers, our tool characterizes the frequency of close flyby interactions in the Universe. We find that the number of close flyby interactions is comparable to, or even surpasses, the number of mergers for halo masses $\gtrsim 10^{11}\,\Msun$ at $z \lesssim 2$. Halo flybys occur so frequently to high mass halos that they are continually perturbed, unable to reach a dynamical equilibrium. In particular, we find that Milky Way type halos undergo a similar number of flybys as mergers irrespective of mass-ratio for $z\lesssim 2$. We also find tentative evidence that at high redshift, $z \gtrsim 14$, flybys are as frequent as mergers. Our results suggest that close halo flybys can play an important role in the evolution of the earliest dark matter halos and their galaxies, and can still influence galaxy evolution at the present epoch. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Identifying High Metallicity M Giants at Intragroup Distances with SDSS

Tidal stripping and three-body interactions with the central supermassive black hole may eject stars from the Milky Way. These stars would comprise a set of `intragroup' stars that trace the past history of interactions in our galactic neighborhood. Using the Sloan Digital Sky Survey DR7, we identify candidate solar metallicity red giant intragroup stars using color cuts that are designed to exclude nearby M and L dwarfs. We present 677 intragroup candidates that are selected between 300 kpc and 2 Mpc, and are either the reddest intragroup candidates (M7-M10) or are L dwarfs at larger distances than previously detected.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Kickstarting Reionization with the First Black Holes: the effects of second-order perturbation theory in pre-reionization volumes

We explore structure formation in the dark ages ($z\sim 30-6$) using two well-known methods for initializing cosmological $N$-body simulations. Overall, both the Zel'dovich approximation (\za) and second order Lagrangian perturbation theory (\lpt) are known to produce accurate present-day dark matter halo mass functions. However, since the \lpt method drives more rapid evolution of dense regions, it increases the occurrence of rare massive objects -- an effect that is most pronounced at high redshift. We find that \lpt produces more halos that could harbor Population III stars and their black hole remnants, and they produce them earlier. Although the differences between the \lpt and \za mass functions are nearly erased by $z=6$, this small boost to the number and mass of black holes more than doubles the reionized volume of the early Universe. We discuss the implications for reionization and massive black hole growth.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The HI environment of the M101 group

We present a wide (8.5x6.7 degree, 1050x825 kpc), deep (sigma(N_HI)=10^(16.8-17.5) cm^-2) neutral hydrogen (HI) map of the M101 galaxy group. We identify two new HI sources in the group environment, one an extremely low surface brightness (and hitherto unknown) dwarf galaxy, and the other a starless HI cloud, possibly primordial in origin. Our data show that M101's extended HI envelope (Huchtmeier & Witzel 1979) takes the form of a ~100 kpc long tidal loop or plume of HI extending to the southwest of the galaxy. The plume has an HI mass ~ 10^8 Msun and a peak column density of N_HI=5x10^17 cm^-2, and while it rotates with the main body of M101, it shows kinematic peculiarities suggestive of a warp or flaring out of the rotation plane of the galaxy. We also find two new HI clouds near the plume with masses ~ 10^7 Msun, similar to HI clouds seen in the M81/M82 group, and likely also tidal in nature. Comparing to deep optical imaging of the M101 group, neither the plume nor the clouds have any extended optical counterparts down to a limiting surface brightness of mu_B = 29.5. We also trace HI at intermediate velocities between M101 and NGC 5474, strengthening the case for a recent interaction between the two galaxies. The kinematically complex HI structure in the M101 group, coupled with the optical morphology of M101 and its companions, suggests that the group is in a dynamically active state that is likely common for galaxies in group environments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Lives of High Redshift Mergers

We present a comparative study of recent works on merger-timescales with dynamical friction and find a strong contrast between idealized/isolated mergers (Boylan-Kolchin et al. 2008) and mergers from a cosmological volume (Jiang et al. 2008). Our study measures the duration of mergers in a cosmological N-body simulation of dark matter, with emphasis on higher redshifts (z < 10) and a lower mass range. In our analysis we consider and compare two merger definitions; tidal disruption and coalescence. We find that the merger-time formula proposed by Jiang et al. (2008) describes our results well and conclude that cosmologically motivated merger-time formulae provide a more versatile and statistically robust approximation for practical applications such as semi-analytic/hybrid models.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Growing Massive Black Holes in a Local Group Environment: the Central Supermassive, Slowly Sinking, and Ejected Populations

We explore the growth of < 10^7 Msun black holes that reside at the centers of spiral and field dwarf galaxies in a Local Group type of environment. We use merger trees from a cosmological N-body simulation known as Via Lactea II (VL-2) as a framework to test two merger-driven semi-analytic recipes for black hole growth that include dynamical friction, tidal stripping, and gravitational wave recoil in over 20,000 merger tree realizations. First, we apply a Fundamental Plane limited (FPL) model to the growth of Sgr A*, which drives the central black hole to a maximum mass limited by the Black Hole Fundamental Plane after every merger. Next, we present a new model that allows for low-level Prolonged Gas Accretion (PGA) during the merger. We find that both models can generate a Sgr A* mass black hole. We predict a population of massive black holes in local field dwarf galaxies - if the VL-2 simulation is representative of the growth of the Local Group, we predict up to 35 massive black holes (< 10^6 Msun) in Local Group field dwarfs. We also predict that hundreds of < 10^5 Msun black holes fail to merge, and instead populate the Milky Way halo, with the most massive of them at roughly the virial radius. In addition, we find that there may be hundreds of massive black holes ejected from their hosts into the nearby intergalactic medium due to gravitational wave recoil. We discuss how the black hole population in the Local Group field dwarfs may help to constrain the growth mechanism for Sgr A*.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Origin of Neutral Hydrogen Clouds in Nearby Galaxy Groups: Exploring the Range Of Galaxy Interactions

We combine high resolution N-body simulations with deep observations of neutral hydrogen (HI) in nearby galaxy groups in order to explore two well-known theories of HI cloud formation: HI stripping by galaxy interactions and dark matter minihalos with embedded HI gas. This paper presents new data from three galaxy groups, Canes Venatici I, NGC 672, and NGC 45, and assembles data from our previous galaxy group campaign to generate a rich HI cloud archive to compare to our simulated data. We find no HI clouds in the Canes Venatici I, NGC 672, or NGC 45 galaxy groups. We conclude that HI clouds in our detection space are most likely to be generated through recent, strong galaxy interactions. We find no evidence of HI clouds associated with dark matter halos above M_HI = 10^6 M_Sun, within +/- 700 km/s of galaxies, and within 50 kpc projected distance of galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Balancing the Baryon Budget: The fraction of the IGM due to Galaxy Mergers

Observations indicate that roughly 60% of the baryons may exist in a Warm-Hot Intergalactic Medium (WHIM) at low redshifts. Following up on previous results showing that gas is released through galaxy mergers, we use a semi-analytic technique to estimate the fraction of gas mass lost from haloes solely due to mergers. We find that up to ~25% of the gas in a halo can unbind over the course of galaxy assembly. This process does not act preferentially on smaller mass haloes; bigger haloes \emph{always} release larger amounts of gas in a given volume of the Universe. However, if we include multi-phase gas accretion onto haloes, we find that only a few percent is unbound. We conclude that either non-gravitational processes may be in play to heat up the gas in the galaxies prior to unbinding by mergers or most of the baryons in the WHIM have never fallen into virialised dark matter haloes. We present a budget for stocking the WHIM compiled from recent work.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Gravitational Wave Signal from Assembling the Lightest Supermassive Black Holes

We calculate the gravitational wave signal from the growth of 10 million solar mass supermassive black holes (SMBH) from the remnants of Population III stars. The assembly of these lower mass black holes is particularly important because observing SMBHs in this mass range is one of the primary science goals for the Laser Interferometer Space Antenna (LISA), a planned NASA/ESA mission to detect gravitational waves. We use high resolution cosmological N-body simulations to track the merger history of the host dark matter halos, and model the growth of the SMBHs with a semi-analytic approach that combines dynamical friction, gas accretion, and feedback. We find that the most common source in the LISA band from our volume consists of mergers between intermediate mass black holes and SMBHs at redshifts less than 2. This type of high mass ratio merger has not been widely considered in the gravitational wave community; detection and characterization of this signal will likely require a different technique than is used for SMBH mergers or extreme mass ratio inspirals. We find that the event rate of this new LISA source depends on prescriptions for gas accretion onto the black hole as well as an accurate model of the dynamics on a galaxy scale; our best estimate yields about 40 sources with a signal-to-noise ratio greater than 30 occur within a volume like the Local Group during SMBH assembly -- extrapolated over the volume of the universe yields roughly 500 observed events over 10 years, although the accuracy of this rate is affected by cosmic variance.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

HI Clouds in the M81 Filament as Dark Matter Minihalos--A Phase-Space Mismatch

Cosmological galaxy formation models predict the existence of dark matter minihalos surrounding galaxies and in filaments connecting groups of galaxies. The more massive of these minihalos are predicted to host HI gas that should be detectable by current radio telescopes such as the GBT. We observed the region including the M81/M82 and NGC 2403 galaxy groups, searching for observational evidence of an HI component associated with dark matter halos within the "M81 Filament", using the Robert C. Byrd Green Bank Telescope (GBT). The map covers an 8.7 degree x 21.3 degree (480 kpc x 1.2 Mpc) region centered between the M81/M82 and NGC 2403 galaxy groups. Our observations cover a wide velocity range, from -890 to 1320 km/s, which spans much of the range predicted by cosmological N-body simulations for dark matter minihalo velocities. Our search is not complete in the velocity range -210 to 85 km/s, containing Galactic emission and the HVC Complex A. For an HI cloud at the distance of M81, with a size < 10 kpc, our average 5-sigma mass detection limit is 3.2 x 10^6 M_Sun, for a linewidth of 20 km/s. We compare our observations to two large cosmological N-body simulations and find that the simulation predicts a significantly greater number of detectable minihalos than are found in our observations, and that the simulated minihalos do not match the phase space of observed HI clouds. These results place strong constraints on the HI gas that can be associated with dark-matter halos. Our observations indicate that the majority of extragalactic HI clouds with a mass greater than 10^6 M_Sun are likely to be generated through tidal stripping caused by galaxy interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Interpretation of the 115 Day Periodic Modulation in the X-ray Flux of NGC 5408 X-1

We comment on the recent observation of a 115-day modulation in the X-ray flux of the ultraluminous X-ray source (ULX) NGC 5408 X-1, and in particular, the interpretation of this modulation as the orbital period. We suggest that this modulation may instead be due to a precessing jet, and is thus superorbital in nature. Comparing the properties of this ULX with those of the prototypical microquasar SS 433, we argue that NGC 5408 X-1 is very similar to SS 433: a hyper-accreting stellar mass black hole in a shorter-period binary. If the analogy holds, the 115-day modulation is best explained by the still poorly-understood physics of inner-disc/jet precession and a longer observing baseline would be able to reveal an intrinsic phase jitter that is associated with such a precession.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Numerical Simulations of Hot Halo Gas in Galaxy Mergers

Galaxy merger simulations have explored the behaviour of gas within the galactic disk, yet the dynamics of hot gas within the galaxy halo has been neglected. We report on the results of high-resolution hydrodynamic simulations of colliding galaxies with hot halo gas. We explore a range of mass ratios, gas fractions and orbital configurations to constrain the shocks and gas dynamics within the progenitor haloes. We find that : (i) A strong shock is produced in the galaxy haloes before the first passage, increasing the temperature of the gas by almost an order of magnitude to $T\sim 10^{6.3}$ K. (ii) The X-ray luminosity of the shock is strongly dependent on the gas fraction; it is $\gtrsim 10^{39}$ erg/s for halo gas fractions larger than 10%. (iii) The hot diffuse gas in the simulation produces X-ray luminosities as large as $10^{42}$ erg/s. This contributes to the total X-ray background in the Universe. (iv) We find an analytic fit to the maximum X-ray luminosity of the shock as a function of merger parameters. This fit can be used in semi-analytic recipes of galaxy formation to estimate the total X-ray emission from shocks in merging galaxies. (v) $\sim$ 10-20% of the initial gas mass is unbound from the galaxies for equal-mass mergers, while $3-5%$ of the gas mass is released for the 3:1 and 10:1 mergers. This unbound gas ends up far from the galaxy and can be a feasible mechanism to enrich the IGM with metals.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Search for Diffuse Neutral Hydrogen and HI Clouds in the NGC 2403 Group

We have observed the NGC 2403 group of galaxies using the Robert C. Byrd Green Bank Telescope (GBT) in a search for faint, extended neutral hydrogen clouds similar to the clouds found around the M81/M82 group, which is located approximately 250 kpc from the NGC 2403 group along the same filament of galaxies. For an HI cloud with a size < 10 kpc within 50 kpc of a group galaxy, our 7-sigma mass detection limit is 2.2 x 10^6 M_sun for a cloud with a linewidth of 20 km/s, over the velocity range from -890 to 1750 km/s. At this sensitivity level we detect 3 new HI clouds in the direction of the group, as well as the known galaxies. The mean velocity of the new clouds differs from that of the group galaxies by more than 250 km/s, but are in the range of Milky Way High Velocity Clouds (HVCs) in that direction. It is most likely that the clouds are part of the Milky Way HVC population. If HI clouds exist in the NGC 2403 group, their masses are less than 2.2 x 10^6 M_sun. We also compared our results to structures that are expected based on recent cosmological models, and found none of the predicted clouds. If NGC 2403 is surrounded by a population of dark matter halos similar to those proposed for the Milky Way in recent models, our observations imply that their HI content is less than 1% of their total mass.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Event Rate for Extreme Mass Ratio Burst Signals in the LISA Band

Stellar mass compact objects in short period orbits about a $10^{4.5}$--$10^{7.5}$ solar mass massive black hole (MBH) are thought to be a significant continuous-wave source of gravitational radiation for the ESA/NASA Laser Interferometer Space Antenna (LISA) gravitational wave detector. However, these extreme mass-ratio inspiral sources began in long-period, nearly parabolic orbits that have multiple close encounters with the MBH. The gravitational radiation emitted during these close encounters may be detectable by LISA as a gravitational wave burst if the characteristic passage timescale is less than $10^5$ seconds. Scaling a static, spherical model to the size and mass of the Milky Way bulge we estimate an event rate of ~ 15 per year for such burst signals, detectable by LISA with signal-to-noise greater than five, originating in our galaxy. When extended to include Virgo cluster galaxies our estimate increases to a gravitational wave burst rate of ~ 18. We conclude that these extreme mass-ratio burst sources may be a steady and significant source of gravitational radiation in the LISA data streams.

0 citations0 reviewsNo code linkedOpen access