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

Pau Amaro-Seoane

Pau Amaro-Seoane appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAgr-qcastro-ph.COastro-ph.HEastro-phastro-ph.SRastro-ph.EPastro-ph.IMhep-th

Catalog footprint

What is connected

50works
9topics
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

50 published item(s)

preprint2020arXiv

Improved gravitational radiation time-scales: significance for LISA and LIGO-Virgo sources

We present a revised version of Peters' (1964) time-scale for the gravitational-wave (GW) induced decay of two point masses. The new formula includes the effects of the first-order post-Newtonian perturbation and additionally provides a simple fit to account for the Newtonian self-consistent evolution of the eccentricity. The revised time-scale is found by multiplying Peters' estimate by two factors, $R(e_0)= 8^{1-\sqrt{1-e_0}}$ and $Q_{\rm f}(p_0) = \exp \left(2.5 (r_{\rm S}/p_0) \right)$, where $e_0$ and $p_0$ are the initial eccentricity and periapsis, respectively, and $r_{\rm S}$ the Schwarzschild radius of the system. Their use can correct errors of a factor of 1-10 that arise from using the original Peters' formula. We apply the revised time-scales to a set of typical sources for existing ground-based laser interferometers and for the future Laser Interferometer Space Antenna (LISA), at the onset of their GW driven decay. We argue that our more accurate model for the orbital evolution will affect current event- and detection-rate estimates for mergers of compact object binaries, with stronger deviations for eccentric LISA sources, such as extreme and intermediate mass-ratio inspirals. We propose the correction factors $R$ and $Q_{\rm f}$ as a simple prescription to quantify decay time-scales more accurately in future population synthesis models. We also suggest that the corrected time-scale may be used as a computationally efficient alternative to numerical integration in other applications that include the modelling of radiation reaction for eccentric sources.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Phase shift of gravitational waves induced by aberration

The velocity of a gravitational wave (GW) source provides crucial information about its formation and evolution processes. Previous studies considered the Doppler effect on the phase of GWs as a potential signature of a time-dependent velocity of the source. However, the Doppler shift only accounts for the time component of the wave vector, and in principle motion also affects the spatial components. The latter effect, known as ``aberration'' for light, is analyzed in this paper for GWs and applied to the waveform modeling of an accelerating source. We show that the additional aberrational phase shift could be detectable in two astrophysical scenarios, namely, a recoiling binary black hole (BBH) due to GW radiation and a BBH in a triple system. Our results suggest that adding the aberrational phase shift in the waveform templates could significantly enhance the detectability of moving sources.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Rates of Stellar Tidal Disruption

Tidal disruption events occur rarely in any individual galaxy. Over the last decade, however, time-domain surveys have begun to accumulate statistical samples of these flares. What dynamical processes are responsible for feeding stars to supermassive black holes? At what rate are stars tidally disrupted in realistic galactic nuclei? What may we learn about supermassive black holes and broader astrophysical questions by estimating tidal disruption event rates from observational samples of flares? These are the questions we aim to address in this Chapter, which summarizes current theoretical knowledge about rates of stellar tidal disruption, and compares theoretical predictions to the current state of observations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Detecting the Beaming Effect of Gravitational Waves

The models currently used in the detection of gravitational waves (GWs) either do not consider a relative motion between the center-of-mass of the source and the observer, or usually only consider its effect on the frequencies of GWs. However, it is known for light waves that a relative motion not only changes the frequencies but also the brightness of the source, the latter of which is called the `beaming effect'. Here we investigate such an effect for GWs and find that the observed amplitude of a GW signal, unlike the behavior of light, is not a monotonic function of the relative velocity and responds differently to the two GW polarizations. We attribute the difference to a rotation of the wave-vector, as well as a reorientation of the GW polarizations. We find that even for velocities as small as $0.25\%$ of the speed of light, ignoring the aforementioned beaming effect could induce a systematic error that is larger than the designated calibration accuracy of LIGO. This error could lead to an incorrect estimation of the distance and orbital inclination of a GW source, or result in a spurious signal that appears to be incompatible with general relativity.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Making Bright Giants Invisible At The Galactic Centre

Current observations of the Galactic Center (GC) seem to display a core-like distribution of bright stars from $\sim 5"$ inwards. On the other hand, we observe young, massive stars at the GC, with roughly 20-50\% of them in a disc, mostly in the region where the bright giants appear to be lacking. In a previous publication we put the idea forward that the missing stars are deeply connected to the presence of this disc. The progenitor of the stellar disc is very likely to have been a gaseous disc that at some point fragmented and triggered star formation. This caused the appearance of overdensity regions in the disc that had high enough densities to ensure stripping large giants of their atmospheres and thus rendering them very faint. In this paper we use a stellar evolution code to derive the properties that a red giant would display in a colour-magnitude diagram, as well as a non-linearity factor required for a correct estimate of the mass loss. We find that in a very short timescale, the red giants (RGs) leave their standard evolutionary track. The non-linearity factor has values that not only depend on the properties of the clumps, but also on the physical conditions the giant stars, as we predicted analytically. According to our results, envelope stripping works, moving stars on a short timescale from the giant branch to the white dwarf stage, thus rendering them invisible to observations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Our Supermassive Black Hole Rivaled the Sun in the Ancient X-ray Sky

Sagittarius A* (SgrA*) lying in the Galactic Centre $8$ kpc from Earth, hosts the closest supermassive black hole known to us. It is now inactive, but there is evidence indicating that about six million years ago it underwent a powerful outburst where the luminosity could have approached the Eddington limit. Motivated by the fact that in extragalaxies the supermassive black holes with similar masses and near-Eddington luminosities are usually strong X-ray emitters, we calculate here the X-ray luminosity of SgrA*. For that, we assume that the outburst was due to accretion of gas or the tidal disruption of a star. We show that these cases could precipitate on Earth a hard X-ray (i.e. $hν>2~{\rm keV}$) flux comparable to that from the current quiescent sun. The flux in harder energy band $20~{\rm keV}<hν<100~{\rm keV}$, however, surpasses that from an X-class solar flare, and the irradiation timescale is also much longer, ranging from weeks to $10^5$ years depending on the outburst scenario. In the solar system gas giants will suffer the biggest impact in their atmospheres. Lower-mass planets such as Earth receive a level of radiation that might have played a role in the evolution of their primitive atmospheres, so that a detailed study of the consequences deserves further investigation. Planetary systems closer to SgrA* receive higher irradiance levels, making them more likely uninhabitable.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Relativistic mergers of black hole binaries have large, similar masses, low spins and are circular

Gravitational waves are a prediction of general relativity, and with ground-based detectors now running in their advanced configuration, we will soon be able to measure them directly for the first time. Binaries of stellar-mass black holes are among the most interesting sources for these detectors. Unfortunately, the many different parameters associated with the problem make it difficult to promptly produce a large set of waveforms for the search in the data stream. To reduce the number of templates to develop, one must restrict some of the physical parameters to a certain range of values predicted by either (electromagnetic) observations or theoretical modeling. In this work we show that "hyperstellar" black holes (HSBs) with masses $30 \lesssim M_{\rm BH}/M_{\odot} \lesssim 100$, i.e black holes significantly larger than the nominal $10\,M_{\odot}$, will have an associated low value for the spin, i.e. $a<0.5$. We prove that this is true regardless of the formation channel, and that when two HSBs build a binary, each of the spin magnitudes is also low, and the binary members have similar masses. We also address the distribution of the eccentricities of HSB binaries in dense stellar systems using a large suite of three-body scattering experiments that include binary-single interactions and long-lived hierarchical systems with a highly accurate integrator, including relativistic corrections up to ${\cal O}(1/c^5)$. We find that most sources in the detector band will have nearly zero eccentricities. This correlation between large, similar masses, low spin and low eccentricity will help to accelerate the searches for gravitational-wave signals.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Retrograde binaries of massive black holes in circum-binary accretion discs

We explore the hardening of a massive black hole binary embedded in a circum-binary gas disc when the binary and the gas are coplanar and the gas is counter-rotating. The secondary black hole, revolving in the direction opposite to the gas, experiences a drag from gas-dynamical friction and from direct accretion of part of it. Using two-dimensional (2D) hydrodynamical grid simulations we investigate the effect of changing the accretion prescriptions on the dynamics of the secondary black hole which in turn affect the binary hardening and eccentricity evolution. We find that realistic accretion prescriptions lead to results that differ from those inferred assuming accretion of all the gas within the Roche Lobe of the secondary black hole. Different accretion prescriptions result in different disc's surface densities which alter the black hole's dynamics back. Full 3D SPH realizations of a number of representative cases, run over a shorter interval of time, validate the general trends observed in the less computationally demanding 2D simulations. Initially circular black hole binaries increase only slightly their eccentricity which then oscillates around small values (<0.1) while they harden. By contrast, initially eccentric binaries become more and more eccentric. A semi-analytical model describing the black hole's dynamics under accretion only explores the late evolution stages of the binary in an otherwise unperturbed retrograde disc to illustrate how eccentricity evolves with time in relation to the shape of the underlying surface density distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Probing dark matter crests with white dwarfs and IMBHs

White dwarfs (WDs) are the most promising captors of dark matter (DM) particles in the crests that are expected to build up in the cores of dense stellar clusters. The DM particles could reach sufficient densities in WD cores to liberate energy through self-annihilation. The extinction associated with our Galactic Centre, the most promising region where to look for such effects, makes it impossible to detect the potential associated luminosity of the DM-burning WDs. However, in smaller stellar systems which are close enough to us and not heavily extincted, such as $ω-$Cen, we may be able to detect DM-burning WDs. We investigate the prospects of detection of DM-burning WDs in a stellar cluster harbouring an IMBH, which leads to higher densities of DM at the centre as compared with clusters without one. We calculate the capture rate of WIMPs by a WD around an IMBH and estimate the luminosity that a WD would emit depending on its distance to the center of the cluster. Direct-summation $N-$body simulations of $ω-$Cen yield a non-negligible number of WDs in the range of radii of interest. We apply our assumption to published HST/ACS observations of stars in the center of $ω-$Cen to search for DM burning WDs and, although we are not able to identify any evident candidate because of crowding and incompleteness, we proof that their bunching up at high luminosities would be unique. We predict that DM burning will lead to a truncation of the cooling sequence at the faint end. The detection of DM burning in future observations of dense stellar clusters, such as globular clusters or ultra-compact dwarf galaxies could allow us to probe different models of DM distributions and characteristics such as the DM particle scattering cross section on nucleons. On the other hand, if DM-burning WDs really exist, their number and properties could give hints to the existence of IMBHs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stability of Gas Clouds in Galactic Nuclei: An Extended Virial Theorem

Cold gas entering the central $1$ to $10^2$ pc of a galaxy fragments and condenses into clouds. The stability of the clouds determines whether they will be turned into stars or can be delivered to the central supermassive black hole (SMBH) to turn on an active galactic nucleus (AGN). The conventional criteria to assess the stability of these clouds, such as the Jeans criterion and Roche (or tidal) limit, are insufficient here, because they assume the dominance of self-gravity in binding a cloud, and neglect external agents, such as pressure and tidal forces, which are common in galactic nuclei. We formulate a new scheme for judging this stability. We first revisit the conventional Virial theorem, taking into account an external pressure, to identify the correct range of masses that lead to stable clouds. We then extend the theorem to include an external tidal field, crucial for the stability in the region of interest -- in dense star clusters, around SMBHs. We apply our extended Virial theorem to find the correct solutions to practical problems that until now were controversial, namely, the stability of the gas clumps in AGN tori, the circum-nuclear disk in the Galactic Center, and the central molecular zone of the Milky Way. The masses we derive for these structures are orders of magnitude smaller than the commonly-used Virial masses (equivalent to the Jeans mass). Moreover, we prove that these clumps are stable, contrary to what one would naively deduce from the Roche (tidal) limit.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Super massive black holes in star forming gaseous circumnuclear discs

Using N-body/SPH simulations we study the evolution of the separation of a pair of SMBHs embedded in a star forming circumnuclear disk (CND). This type of disk is expected to be formed in the central kilo parsec of the remnant of gas-rich galaxy mergers. Our simulations indicate that orbital decay of the SMBHs occurs more quickly when the mean density of the CND is higher, due to increased dynamical friction. However, in simulations where the CND is fragmented in high density gaseous clumps (clumpy CND), the orbits of the SMBHs are erratically perturbed by the gravitational interaction with these clumps, delaying, in some cases, the orbital decay of the SMBHs. The densities of these gaseous clumps in our simulations and in recent studies of clumpy CNDs are significantly higher than the observed density of molecular clouds in isolated galaxies or ULIRGs, thus, we expect that SMBH orbits are perturbed less in real CNDs than in the simulated CNDs of this study and other recent studies. We also find that the migration timescale has a weak dependence on the star formation rate of the CND. Furthermore, the migration timescale of a SMBH pair in a star-forming clumpy CND is at most a factor three longer than the migration timescale of a pair of SMBHs in a CND modeled with more simple gas physics. Therefore, we estimate that the migration timescale of the SMBHs in a clumpy CND is on the order of $10^7$ yrs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A rapid evolving region in the Galactic Center: Why S-stars thermalize and more massive stars are missing

The existence of "S-stars" within a distance of 1" from SgrA$^*$ contradicts our understanding of star formation, due to the forbiddingly violent environment. A suggested possibility is that they form far and have been brought in by some fast dynamical process, since they are young. Nonetheless, all conjectured mechanisms either fail to reproduce their eccentricities --without violating their young age-- or cannot explain the problem of "inverse mass segregation": The fact that lighter stars (the S-stars) are closer to SgrA$^*$ and more massive ones, Wolf-Rayet (WR) and O-stars, are farther out. In this Letter we propose that the responsible for both, the distribution of the eccentricities and the paucity of massive stars, is the Kozai-Lidov-{\em like} resonance induced by a sub-parsec disk recently discovered in the Galactic center. Considering that the disk probably extended to smaller radius in the past, we show that in as short as (a few) $10^6$ years, the stars populating the innermost 1" region would redistribute in angular-momentum space and recover the observed "super-thermal" distribution. Meanwhile, WR and O-stars in the same region intermittently attain ample eccentricities that will lead to their tidal disruptions by the central massive black hole. Our results provide new evidences that SgrA$^*$ was powered several millions years ago by an accretion disk as well as by tidal stellar disruptions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Descope of the ALIA mission

The present work reports on a feasibility study commissioned by the Chinese Academy of Sciences of China to explore various possible mission options to detect gravitational waves in space alternative to that of the eLISA/LISA mission concept. Based on the relative merits assigned to science and technological viability, a few representative mission options descoped from the ALIA mission are considered. A semi-analytic Monte Carlo simulation is carried out to understand the cosmic black hole merger histories starting from intermediate mass black holes at high redshift as well as the possible scientific merits of the mission options considered in probing the light seed black holes and their coevolution with galaxies in early Universe. The study indicates that, by choosing the armlength of the interferometer to be three million kilometers and shifting the sensitivity floor to around one-hundredth Hz, together with a very moderate improvement on the position noise budget, there are certain mission options capable of exploring light seed, intermediate mass black hole binaries at high redshift that are not readily accessible to eLISA/LISA, and yet the technological requirements seem to within reach in the next few decades for China.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Disruption of a Red Giant Star by a Supermassive Black Hole and the Case of PS1-10jh

The development of a new generation of theoretical models for tidal disruptions is timely, as increasingly diverse events are being captured in surveys of the transient sky. Recently, Gezari et al. reported a discovery of a new class of tidal disruption events: the disruption of a helium-rich stellar core, thought to be a remnant of a red giant (RG) star. Motivated by this discovery and in anticipation of others, we consider tidal interaction of an RG star with a supermassive black hole (SMBH) which leads to the stripping of the stellar envelope and subsequent inspiral of the compact core toward the black hole. Once the stellar envelope is removed the inspiral of the core is driven by tidal heating as well as the emission of gravitational radiation until the core either falls into the SMBH or is tidally disrupted. In the case of tidal disruption candidate PS1-10jh we find that there is a set of orbital solutions at high eccentricities in which the tidally stripped hydrogen envelope is accreted by the SMBH before the helium core is disrupted. This places the RG core in a portion of parameter space where strong tidal heating can lift the degeneracy of the compact remnant and disrupt it before it reaches the tidal radius. We consider how this sequence of events explains the puzzling absence of the hydrogen emission lines from the spectrum of PS1-10jh and gives rise to its other observational features.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rapid Eccentricity Oscillations and the Mergers of Compact Objects in Hierarchical Triples

Kozai-Lidov (KL) oscillations can accelerate compact object mergers via gravitational wave (GW) radiation by driving the inner binaries of hierarchical triples to high eccentricities. We perform direct three-body integrations of high mass ratio compact object triple systems using Fewbody including post-Newtonian terms. We find that the inner binary undergoes rapid eccentricity oscillations (REOs) on the timescale of the outer orbital period which drive it to higher eccentricities than secular theory would otherwise predict, resulting in substantially reduced merger times. For a uniform distribution of tertiary eccentricity ($e_2$), ~40% of systems merge within ~1-2 eccentric KL timescales whereas secular theory predicts that only ~20% of such systems merge that rapidly. This discrepancy becomes especially pronounced at low $e_2$, with secular theory overpredicting the merger time by many orders of magnitude. We show that a non-negligible fraction of systems have eccentricity > 0.8 when they merge, in contrast to predictions from secular theory. Our results are applicable to high mass ratio triple systems containing black holes or neutron stars. In objects in which tidal effects are important, such as white dwarfs, stars, and planets, REOs can reduce the tidal circularization timescale by an order of magnitude and bring the components of the inner binary into closer orbits than would be possible in the secular approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Relativistic mergers of compact binaries in clusters: The fingerprint of the spin

Dense stellar systems such as globular clusters and dense nuclear clusters are the breeding ground of sources of gravitational waves for the advanced detectors LIGO and Virgo. These systems deserve a close study to estimate rates and parameter distribution. This is not an easy task, since the evolution of a dense stellar cluster involves the integration of $N$ bodies with high resolution in time and space and including hard binaries and their encounters and, in the case of gravitational waves (GWs), one needs to take into account important relativistic corrections. In this work we present the first implementation of the effect of spin in mergers in a direct-summation code, NBODY6. We employ non-spinning post-Newtonian corrections to the Newtonian accelerations up to 3.5 post-Newtonian (PN) order as well as the spin-orbit coupling up to next-to-lowest order and the lowest order spin-spin coupling. We integrate spin precession and add a consistent treatment of mergers. We analyse the implementation by running a set of two-body experiments and then we run a set of 500 simulations of a relativistic stellar cluster. In spite of the large number of mergers in our tests, the application of the algorithm is robust. We find in particular the formation of a runaway star whose spin decays with the mass it wins, independently of the initial value of the spins of the stars. More remarkably, the subset of compact objects that do not undergo many mergers, and hence represent a more realistic system, has a correlation between the final absolute spin and the initial choice for the initial distribution, which could provide us with information about the evolution of spins in dense clusters once the first detections have started.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Research Update on Extreme-Mass-Ratio Inspirals

The inspirals of stellar-mass mass compact objects into massive black holes in the centres of galaxies are one of the most important sources of gravitational radiation for space-based detectors like LISA or eLISA. These extreme-mass-ratio inspirals (EMRIs) will enable an ambitious research program with implications for astrophysics, cosmology, and fundamental physics. This article is a summary of the talks delivered at the plenary session on EMRIs at the 10th International LISA Symposium. It contains research updates on the following topics: astrophysics of EMRIs; EMRI science potential; and EMRI modeling.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sculpting the Stellar Cusp in the Galactic Center

Observations of the innermost parsec surrounding Sgr A$^*$ ---the supermassive black hole in the center of our Galaxy--- have revealed a diversity of structures whose existence and characteristics apparently defy the fundamental principles of dynamics. In this article, we review the challenges to the dynamics theories that have been brought forth in the past two decades by the observations of the Galactic center (GC). We outline the theoretical framework that has been developed to reconcile the discrepancies between the theoretical predictions and the observational results. In particular, we highlight the role of the recently discovered sub-parsec stellar disk in determining the dynamics and resolving the inconsistencies. We also discuss the implications for the recent activity of Sgr A$^*$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sowing the seeds of massive black holes in small galaxies: Young clusters as the building blocks of Ultra-Compact-Dwarf Galaxies

Interacting galaxies often have complexes of hundreds of young stellar clusters of individual masses $\sim 10^{4-6}~M_\odot$ in regions that are a few hundred parsecs across. These cluster complexes interact dynamically, and their coalescence is a candidate for the origin of some ultracompact dwarf galaxies (UCDs). Individual clusters with short relaxation times are candidates for the production of intermediate-mass black holes of a few hundred solar masses, via runaway stellar collisions prior to the first supernovae in a cluster. It is therefore possible that a cluster complex hosts multiple intermediate-mass black holes that may be ejected from their individual clusters due to mergers or binary processes, but bound to the complex as a whole. Here we explore the dynamical interaction between initially free-flying massive black holes and clusters in an evolving cluster complex. We find that, after hitting some clusters, it is plausible that the massive black hole will be captured in an ultracompact dwarf forming near the center of the complex. In the process, the hole typically triggers electromagnetic flares via stellar disruptions, and is also likely to be a prominent source of gravitational radiation for the advanced ground-based detectors LIGO and VIRGO. We also discuss other implications of this scenario, notably that the central black hole could be considerably larger than expected in other formation scenarios for ultracompact dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Tidal disruptions in circumbinary disks. II: Observational signatures in the reverberation spectra

Supermassive Binary Black Holes (SMBBHs) with sub-pc separations form in the course of galaxy mergers, if both galaxies harbour massive black holes. Clear observational evidence for them however still eludes us. We propose a novel method of identifying these systems by means of reverberation mapping their circumbinary disk after a tidal disruption event has ionized it. The tidal disruption of a star at the secondary leads to strong asymmetries in the disk response. We model the shape of the velocity--delay maps for various toy disk models and more realistic gas distributions obtained by SPH simulations. The emissivity of the ionized disk is calculated with {\em Cloudy}. We find peculiar asymmetries in the maps for off center ionizing sources that may help us constrain geometrical parameters of a circumbinary disk such as semi-major axis and orbital phase of the secondary, as well as help strengthen the observational evidence for sub-parsec SMBBHs as such.

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

Hybrid methods in planetesimal dynamics: Formation of protoplanetary systems and the mill condition

The formation and evolution of protoplanetary discs remains a challenge from both a theoretical and numerical standpoint. In this work we first perform a series of tests of our new hybrid algorithm presented in Glaschke, Amaro-Seoane and Spurzem 2011 (henceforth Paper I) that combines the advantages of high accuracy of direct-summation N-body methods with a statistical description for the planetesimal disc based on Fokker-Planck techniques. We then address the formation of planets, with a focus on the formation of protoplanets out of planetesimals. We find that the evolution of the system is driven by encounters as well as direct collisions and requires a careful modelling of the evolution of the velocity dispersion and the size distribution over a large range of sizes. The simulations show no termination of the protoplanetary accretion due to gap formation, since the distribution of the planetesimals is only subjected to small fluctuations. We also show that these features are weakly correlated with the positions of the protoplanets. The exploration of different impact strengths indicates that fragmentation mainly controls the overall mass loss, which is less pronounced during the early runaway growth. We prove that the fragmentation in combination with the effective removal of collisional fragments by gas drag sets an universal upper limit of the protoplanetary mass as a function of the distance to the host star, which we refer to as the mill condition.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Investigating the retention of intermediate-mass black holes in star clusters using N-body simulations

Contrary to supermassive and stellar-mass black holes (SBHs), the existence of intermediate-mass black holes (IMBHs) with masses ranging between 10^{2-5} Msun has not yet been confirmed. The main problem in the detection is that the innermost stellar kinematics of globular clusters (GCs) or small galaxies, the possible natural loci to IMBHs, are very difficult to resolve. However, if IMBHs reside in the centre of GCs, a possibility is that they interact dynamically with their environment. A binary formed with the IMBH and a compact object of the GC would naturally lead to a prominent source of gravitational radiation, detectable with future observatories. We use N-body simulations to study the evolution of GCs containing an IMBH and calculate the gravitational radiation emitted from dynamically formed IMBH-SBH binaries and the possibility that the IMBH escapes the GC after an IMBH-SBH merger. We run for the first time direct-summation integrations of GCs with an IMBH including the dynamical evolution of the IMBH with the stellar system and relativistic effects, such as energy loss in gravitational waves (GWs) and periapsis shift, and gravitational recoil. We find in one of our models an intermediate mass-ratio inspiral (IMRI), which leads to a merger with a recoiling velocity higher than the escape velocity of the GC. The GWs emitted fall in the range of frequencies that a LISA-like observatory could detect, like the European eLISA or in mission options considered in the recent preliminary mission study conducted in China. The merger has an impact on the global dynamics of the cluster, as an important heating source is removed when the merged system leaves the GC. The detection of one IMRI would constitute a test of GR, as well as an irrefutable proof of the existence of IMBHs.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Mergers of multi-metallic globular clusters: The role of dynamics

Hubble Space Telescope observations of globular clusters (GCs) in the Antennae galaxy show clusters of clusters, or regions in the galaxy that span hundreds of parsecs, where many of the GCs are doomed to collide, and eventually merge. Several such objects appear likely to present a significant range in ages, hence possibly metallicities, and their merger could plausibly lead to multi-metallic GCs. Here we explore this process with direct-summation N-body simulations with GPU hardware. Our results reveal that colliding GCs with different metallicities and ages can produce a GC with multiplicity and occupation fractions not unlike those observed in multi-metallic clusters. In our simulations, the merged clusters have a phase with a larger amount of flattening than average, as a consequence of rapid rotation- thus suggesting that relatively recent mergers may play a role in producing highly flattened, multi-metallic clusters. We additionally explore the role of the King parameter of the cluster in the occupation fractions with a set of 160 direct-summation simulations and find that for equal size clusters the King parameter of the progenitor clusters determines the occupation fractions in the merger product, while in unequal size mergers the size of the clusters dominates the distribution of stars in the new GC. In particular, we find that the observed distribution of populations in Omega Cen can be described to some extent with our dynamical models.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The fragmenting past of the disk at the Galactic Center : The culprit for the missing red giants

Since 1996 we have known that the Galactic Center (GC) displays a core-like distribution of red giant branch (RGB) stars starting at ~ 10", which poses a theoretical problem, because the GC should have formed a segregated cusp of old stars. This issue has been addressed invoking stellar collisions, massive black hole binaries, and infalling star clusters, which can explain it to some extent. Another observational fact, key to the work presented here, is the presence of a stellar disk at the GC. We postulate that the reason for the missing stars in the RGB is closely intertwined with the disk formation, which initially was gaseous and went through a fragmentation phase to form the stars. Using simple analytical estimates, we prove that during fragmentation the disk developed regions with densities much higher than a homogeneous gaseous disk, i.e. "clumps", which were optically thick, and hence contracted slowly. Stars in the GC interacted with them and in the case of RGB stars, the clumps were dense enough to totally remove their outer envelopes after a relatively low number of impacts. Giant stars in the horizontal branch (HB), however, have much denser envelopes. Hence, the fragmentation phase of the disk must have had a lower impact in their distribution, because it was more difficult to remove their envelopes. We predict that future deeper observations of the GC should reveal less depletion of HB stars and that the released dense cores of RGB stars will still be populating the GC.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Blocking low-eccentricity EMRIs: A statistical direct-summation N-body study of the Schwarzschild barrier

The capture of a compact object in a galactic nucleus by a massive black hole (MBH), an extreme-mass ratio inspiral (EMRI), is the best way to map space and time around it. Recent work on stellar dynamics has demonstrated that there seems to be a complot in phase space acting on low-eccentricity captures, since their rates decrease significantly by the presence of a blockade in the rate at which orbital angular momenta change takes place. This so-called "Schwarzschild barrier" is a result of the impact of relativistic precession on to the stellar potential torques, and thus it affects the enhancement on lower-eccentricity EMRIs that one would expect from resonant relaxation. We confirm and quantify the existence of this barrier using a statistical sample of 2,500 direct-summation N-body simulations using both a post-Newtonian and also for the first time in a direct-summation code a geodesic approximation for the relativistic orbits. The existence of the barrier prevents low-eccentricity EMRIs from approaching the central MBH, but high-eccentricity EMRIs, which have been wrongly classified as "direct plunges" until recently, ignore the presence of the barrier, because they are driven by two-body relaxation. Hence, since the rates are significantly affected in the case of low-eccentricity EMRIs, we predict that a LISA-like observatory such as eLISA will predominantly detect high-eccentricity EMRIs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dynamics of dense star-gas systems: BHs and their precursors

This thesis embraces several aspects of theoretical stellar dynamics in clusters, both analytically and numerically. We try to elucidate the phenomena currently observed in all types of galaxies, including AGNs and quasars, some of the most powerful objects in the universe. The interactions between the stellar system and the central black hole give rise to a lot of interesting phenomena. The scheme we employ enables a study of clean-cut aspects without any noise that particle methods suffer from. We study the most important physical processes that are readily available in the evolution of a spherical cluster, like self-gravity, two-body relaxation etc, the interaction with a central black hole and the role of a mass spectrum. Not only embark we upon this subject, but we set about an analysis on super-massive stars. How these stars could power the quasar activity by star accretion and energy flows is one of the questions that arises. We undertake other questions, such as the uncertain evolution of such an object and its interaction with the surrounding stellar system. This is of crucial importance in astrophysics, for these objects could be regarded as super-massive black holes progenitors.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

eLISA: Astrophysics and cosmology in the millihertz regime

This document introduces the exciting and fundamentally new science and astronomy that the European New Gravitational Wave Observatory (NGO) mission (derived from the previous LISA proposal) will deliver. The mission (which we will refer to by its informal name "eLISA") will survey for the first time the low-frequency gravitational wave band (about 0.1 mHz to 1 Hz), with sufficient sensitivity to detect interesting individual astrophysical sources out to z = 15. The eLISA mission will discover and study a variety of cosmic events and systems with high sensitivity: coalescences of massive black holes binaries, brought together by galaxy mergers; mergers of earlier, less-massive black holes during the epoch of hierarchical galaxy and black-hole growth; stellar-mass black holes and compact stars in orbits just skimming the horizons of massive black holes in galactic nuclei of the present era; extremely compact white dwarf binaries in our Galaxy, a rich source of information about binary evolution and about future Type Ia supernovae; and possibly most interesting of all, the uncertain and unpredicted sources, for example relics of inflation and of the symmetry-breaking epoch directly after the Big Bang. eLISA's measurements will allow detailed studies of these signals with high signal-to-noise ratio, addressing most of the key scientific questions raised by ESA's Cosmic Vision programme in the areas of astrophysics and cosmology. They will also provide stringent tests of general relativity in the strong-field dynamical regime, which cannot be probed in any other way. This document not only describes the science but also gives an overview on the mission design and orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Evolution of binary black holes in self gravitating discs: dissecting the torques

We study the interplay between gas accretion and gravity torques in changing a binary elements and its total angular momentum (L) budget. Especially, we analyse the physical origin of the gravity torques (T_g) and their location within the disc. We analyse 3D SPH simulations of the evolution of initially quasi-circular massive black hole binaries (BHBs) residing in the central hollow of massive self-gravitating circumbinary discs. We use different thermodynamics within the cavity and for the numerical size of the black holes to show that (i) the BHB eccentricity growth found previously is a general result, independent of the accretion and the adopted thermodynamics; (ii) the semi-major axis decay depends both on the T_g and on the interplay with the disc-binary L-transfer due to accretion; (iii) the spectral structure of the T_g is predominately caused by disc edge overdensities and spiral arms developing in the body of the disc and, in general, does not reflect directly the period of the binary; (iv) the net T_g changes sign across the BHB corotation radius. We quantify the relative importance of the two, which appear to depend on the thermodynamical properties of the instreaming gas, and which is crucial in assessing the disc-binary L-transfer; (v) the net torque manifests as a purely kinematic (non-resonant) effect as it stems from the cavity, where the material flows in and out in highly eccentric orbits. Both accretion onto the black holes and the interaction with gas streams inside the cavity must be taken into account to assess the fate of the BHB. Moreover, the total torque exerted by the disc affects L(BHB) by changing all the elements (mass, mass ratio, eccentricity, semimajor axis) of the BHB. Common prescriptions equating tidal torque to semi-major axis shrinking might therefore be poor approximations for real astrophysical systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Fake plunges are very eccentric real EMRIs in disguise ... they dominate the rates and are blissfully ignorant of angular momentum barriers

The capture of a compact object in a galactic nucleus by a massive black hole (MBH) is the best way to map space and time around it. Compact objects such as stellar black holes on a capture orbit with a very high eccentricity have been wrongly assumed to be lost for the system after an intense burst of radiation, which has been described as a "direct plunge". We prove that these very eccentric capture orbits spend actually a similar number of cycles in a LISA-like detector as those with lower eccentricities if the central MBH is spinning. Although the rates are higher for high-eccentricity EMRIs, the spin also enhances the rates of lower-eccentricity EMRIs. This last kind have received more attention because of the fact that high-eccentricity EMRIs were thought to be direct plunges and thus negligible. On the other hand, recent work on stellar dynamics has demonstrated that there seems to be a complot in phase space acting on these lower-eccentricity captures, since their rates decrease significantly by the presence of a blockade in the rate at which orbital angular momenta change takes place. This so-called "Schwarzschild barrier" is a result of the impact of relativistic precession on to the stellar potential torques, and thus it affects the enhancement on lower-eccentricity EMRIs that one would expect from resonant relaxation. We confirm and quantify the existence of this barrier using a statitical sample of 2,500 direct-summation N-body simulations using both a post-Newtonian but also, and for the first time, a geodesic approximation for the relativistic orbits. The existence of the barrier prevents "traditional EMRIs" from approaching the central MBH, but if the central MBH is spinning the rate will be anyway dominated by highly-eccentric extreme-mass ratio inspirals, which insolently ignore the presence of the barrier, because they are driven by two-body relaxation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Low-frequency gravitational-wave science with eLISA/NGO

We review the expected science performance of the New Gravitational-Wave Observatory (NGO, a.k.a. eLISA), a mission under study by the European Space Agency for launch in the early 2020s. eLISA will survey the low-frequency gravitational-wave sky (from 0.1 mHz to 1 Hz), detecting and characterizing a broad variety of systems and events throughout the Universe, including the coalescences of massive black holes brought together by galaxy mergers; the inspirals of stellar-mass black holes and compact stars into central galactic black holes; several millions of ultracompact binaries, both detached and mass transferring, in the Galaxy; and possibly unforeseen sources such as the relic gravitational-wave radiation from the early Universe. eLISA's high signal-to-noise measurements will provide new insight into the structure and history of the Universe, and they will test general relativity in its strong-field dynamical regime.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The role of the supermassive black hole spin in the estimation of the EMRI event rate

A common result to all EMRI investigations on rates is that the possibility that a compact object merges with the MBH after only one intense burst of GWs is much more likely than a slow adiabatic inspiral, an EMRI. The later is referred to as a "plunge" because the compact object dives into the MBH. The event rates for plunges are orders of magnitude larger than slow inspirals. On the other hand, nature MBH's are most likely Kerr and the magnitude of the spin has been sized up to be high. We calculate the number of periapsis passages that a compact object set on to an extremely radial orbit goes through before being actually swallowed by the Kerr MBH and we then translate it into an event rate for a LISA-like observatory, such as the proposed ESA mission eLISA/NGO. We prove that a "plunging" compact object is conceptually indistinguishable from an adiabatic, slow inspiral; plunges spend on average up to hundred of thousands of cycles in the bandwidth of the detector for a two years mission. This has an important impact on the event rate, enhancing in some cases significantly, depending on the spin of the MBH and the inclination. Moreover, it has been recently proved that the production of low-eccentricity EMRIs is severely blocked by the presence of a blockade in the rate at which orbital angular momenta change takes place. This is the result of relativistic precession on to the stellar potential torques and hence affects EMRIs originating via resonant relaxation at distances of about $\sim 10^{-2}$ pc from the MBH. Since high-eccentricity EMRIs are a result of two-body relaxation, they are not affected by this phenomenon. Therefore we predict that eLISA EMRI event rates will be dominated by high-eccentricity binaries, as we present here.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Tidal disruptions of separated binaries in galactic nuclei

Several galaxies have exhibited X-ray flares that are consistent with the tidal disruption of a star by a central supermassive black hole. In theoretical treatments of this process it is usually assumed that the star was initially on a nearly parabolic orbit relative to the black hole. Such an assumption leads in the simplest approximation to a $t^{-5/3}$ decay of the bolometric luminosity and this is indeed consistent with the relatively poorly sampled light curves of such flares. We point out that there is another regime in which the decay would be different: if a binary is tidally separated and the star that remains close to the hole is eventually tidally disrupted from a moderate eccentricity orbit, the decay is slower, typically $\sim t^{-1.2}$. As a result, careful sampling of the light curves of such flares could distinguish between these processes and yield insight into the dynamics of binaries as well as single stars in galactic centres. We explore this process using three-body simulations and analytic treatments and discuss the consequences for present-day X-ray detections and future gravitational wave observations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Gravitational Waves Notes, Issue #5 : "The Capra research programme for capture of small compact objects by massive black holes"

GW Notes was born from the need for a journal where the distinct communities involved in gravitation wave research might gather. While these three communities - astrophysics, general relativity and data analysis - have made significant collaborative progress over recent years, we believe that it is indispensable to future advancement that they draw closer, and that they speak a common idiom. In this 6th GW Note (since we started numbering with #0), we present the work of Jonathan Thornburg, who has been fully-refereed, on the Capra research programme for capture of small compact objects by massive black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hybrid methods in planetesimal dynamics (I) : Description of a new composite algorithm

The formation and evolution of protoplanetary systems, the breeding grounds of planet formation, is a complex dynamical problem that involves many orders of magnitudes. To serve this purpose, we present a new hybrid algorithm that combines a Fokker-Planck approach with the advantages of a pure direct-summation N-body scheme, with a very accurate integration of close encounters for the orbital evolution of the larger bodies with a statistical model, envisaged to simulate the very large number of smaller planetesimals in the disc. Direct-summation techniques have been historically developped for the study of dense stellar systems such as open and globular clusters and, within some limits imposed by the number of stars, of galactic nuclei. The number of modifications to adapt direct-summation N-body techniques to planetary dynamics is not undemanding and requires modifications. These include the way close encounters are treated, as well as the selection process for the "neighbour radius" of the particles and the extended Hermite scheme, used for the very first time in this work, as well as the implementation of a central potential, drag forces and the adjustment of the regularisation treatment. For the statistical description of the planetesimal disc we employ a Fokker-Planck approach. We include dynamical friction, high- and low-speed encounters, the role of distant encounters as well as gas and collisional damping and then generalise the model to inhomogenous discs. We then describe the combination of the two techniques to address the whole problem of planetesimal dynamics in a realistic way via a transition mass to integrate the evolution of the particles according to their masses.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The butterfly effect in the extreme-mass ratio inspiral problem

Measurements of gravitational waves from the inspiral of a stellar-mass compact object into a massive black hole are unique probes to test General Relativity (GR) and MBH properties, as well as the stellar distribution about these holes in galactic nuclei. Current data analysis techniques can provide us with parameter estimation with very narrow errors. However, an EMRI is not a two-body problem, since other stellar bodies orbiting nearby will influence the capture orbit. Any deviation from the isolated inspiral will induce a small, though observable deviation from the idealised waveform which could be misinterpreted as a failure of GR. Based on conservative analysis of mass segregation in a Milky Way like nucleus, we estimate that the possibility that another star has a semi-major axis comparable to that of the EMRI is non-negligible, although probably very small. This star introduces an observable perturbation in the orbit in the case in which we consider only loss of energy via gravitational radiation. When considering the two first-order non-dissipative post-Newtonian contributions (the periapsis shift of the orbit), the evolution of the orbital elements of the EMRI turns out to be chaotic in nature. The implications of this study are twofold. From the one side, the application to testing GR and measuring MBHs parameters with the detection of EMRIs in galactic nuclei with a millihertz mission will be even more challenging than believed. From the other side, this behaviour could in principle be used as a signature of mass segregation in galactic nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The impact of realistic models of mass segregation on the event rate of extreme-mass ratio inspirals and cusp re-growth

One of the most interesting sources of gravitational waves (GWs) for LISA is the inspiral of compact objects on to a massive black hole (MBH), commonly referred to as an "extreme-mass ratio inspiral" (EMRI). The small object, typically a stellar black hole (bh), emits significant amounts of GW along each orbit in the detector bandwidth. The slowly, adiabatic inspiral of these sources will allow us to map space-time around MBHs in detail, as well as to test our current conception of gravitation in the strong regime. The event rate of this kind of source has been addressed many times in the literature and the numbers reported fluctuate by orders of magnitude. On the other hand, recent observations of the Galactic center revealed a dearth of giant stars inside the inner parsec relative to the numbers theoretically expected for a fully relaxed stellar cusp. The possibility of unrelaxed nuclei (or, equivalently, with no or only a very shallow cusp) adds substantial uncertainty to the estimates. Having this timely question in mind, we run a significant number of direct-summation $N-$body simulations with up to half a million particles to calibrate a much faster orbit-averaged Fokker-Planck code. We then investigate the regime of strong mass segregation (SMS) for models with two different stellar mass components. We show that, under quite generic initial conditions, the time required for the growth of a relaxed, mass segregated stellar cusp is shorter than a Hubble time for MBHs with $M_\bullet \lesssim 5 \times 10^6 M_\odot$ (i.e. nuclei in the range of LISA). SMS has a significant impact boosting the EMRI rates by a factor of $\sim 10$ for our fiducial models of Milky Way type galactic nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Constraining scalar fields with stellar kinematics and collisional dark matter

The existence and detection of scalar fields could provide solutions to long-standing puzzles about the nature of dark matter, the dark compact objects at the centre of most galaxies, and other phenomena. Yet, self-interacting scalar fields are very poorly constrained by astronomical observations, leading to great uncertainties in estimates of the mass $m_ϕ$ and the self-interacting coupling constant $λ$ of these fields. To counter this, we have systematically employed available astronomical observations to develop new constraints, considerably restricting this parameter space. In particular, by exploiting precise observations of stellar dynamics at the centre of our Galaxy and assuming that these dynamics can be explained by a single boson star, we determine an upper limit for the boson star compactness and impose significant limits on the values of the properties of possible scalar fields. Requiring the scalar field particle to follow a collisional dark matter model further narrows these constraints. Most importantly, we find that if a scalar dark matter particle does exist, then it cannot account for both the dark-matter halos and the existence of dark compact objects in galactic nuclei

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Detection of IMBHs with ground-based gravitational wave observatories: A biography of a binary of black holes, from birth to death

Even though the existence of intermediate-mass black holes (IMBHs, black holes with masses ranging between $10^{2-4}\,M_{\odot}$) has not yet been corroborated observationally, these objects are of high interest for astrophysics. Our understanding of the formation and evolution of supermassive black holes (SMBHs), as well as galaxy evolution modeling and cosmography would dramatically change if an IMBH were to be observed. From a point of view of traditional photon-based astronomy, {which relies on the monitoring of innermost stellar kinematics}, the {\em direct} detection of an IMBH seems to be rather far in the future. However, the prospect of the detection and characterization of an IMBH has good chances in lower-frequency gravitational-wave (GW) astrophysics using ground-based detectors such as LIGO, Virgo and the future Einstein Telescope (ET). We present an analysis of the signal of a system of a binary of IMBHs (BBH from now onwards) based on a waveform model obtained with numerical relativity simulations coupled with post-Newtonian calculations at the highest available order. IMBH binaries with total masses between $200-20000\,M_\odot$ would produce significant signal-to-noise ratios (SNRs) in advanced LIGO and Virgo and the ET. We have computed the expected event rate of IMBH binary coalescences for different configurations of the binary, finding interesting values that depend on the spin of the IMBHs. The prospects for IMBH detection and characterization with ground-based GW observatories would not only provide us with a robust test of general relativity, but would also corroborate the existence of these systems. Such detections should allow astrophysicists to probe the stellar environments of IMBHs and their formation processes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Gravitational Waves Notes, Issue #2 : "A probe of spacetime and astrophysics: EMRIs"

GW Notes was born from the need for a journal where the distinct communities involved in gravitation wave research might gather. While these three communities - Astrophysics, General Relativity and Data Analysis - have made significant collaborative progress over recent years, we believe that it is indispensable to future advancement that they draw closer, and that they speak a common idiom. For this GW Notes issue we have approached Nicolás Yunes (Princeton University) to extend in high detail his recent work on EMRI waveforms for our highlight article.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Gravitational Waves Notes, Issue #3 : "Stellar cusps in galactic nuclei - How stars distribute around a massive black hole"

GW Notes was born from the need for a journal where the distinct communities involved in gravitation wave research might gather. While these three communities - Astrophysics, General Relativity and Data Analysis - have made significant collaborative progress over recent years, we believe that it is indispensable to future advancement that they draw closer, and that they speak a common idiom. For this GW Notes issue we have approached Miguel Preto (Heidelberg University) to expand a recent work on how stars distribute around massive black holes for our highlight article.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Gravitational Waves Notes, Issue #4 : "A Roadmap to Fundamental Physics from LISA EMRI Observations"

GW Notes was born from the need for a journal where the distinct communities involved in gravitation wave research might gather. While these three communities - Astrophysics, General Relativity and Data Analysis - have made significant collaborative progress over recent years, we believe that it is indispensable to future advancement that they draw closer, and that they speak a common idiom. For this GW Notes issue we have approached Carlos F. Sopuerta to write the highlight article of this issue on his thoughts on fundamental physics with LISA.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Higher order moment models of dense stellar systems: Applications to the modeling of the stellar velocity distribution function

Dense stellar systems such as globular clusters, galactic nuclei and nuclear star clusters are ideal loci to study stellar dynamics due to the very high densities reached, usually a million times higher than in the solar neighborhood; they are unique laboratories to study processes related to relaxation. There are a number of different techniques to model the global evolution of such a system. In statistical models we assume that relaxation is the result of a large number of two-body gravitational encounters with a net local effect. We present two moment models that are based on the collisional Boltzmann equation. By taking moments of the Boltzmann equation one obtains an infinite set of differential moment equations where the equation for the moment of order $n$ contains moments of order $n+1$. In our models we assume spherical symmetry but we do not require dynamical equilibrium. We truncate the infinite set of moment equations at order $n=4$ for the first model and at order $n=5$ for the second model. The collisional terms on the right-hand side of the moment equations account for two-body relaxation and are computed by means of the Rosenbluth potentials. We complete the set of moment equations with closure relations which constrain the degree of anisotropy of our model by expressing moments of order $n+1$ by moments of order $n$. The accuracy of this approach relies on the number of moments included from the infinite series. Since both models include fourth order moments we can study mechanisms in more detail that increase or decrease the number of high velocity stars. The resulting model allows us to derive a velocity distribution function, with unprecedented accuracy, compared to previous moment models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Relativistic encounters in dense stellar systems

Two coalescing black holes (BHs) represent a conspicuous source of gravitational waves (GWs). The merger involves 17 parameters in the general case of Kerr BHs, so that a successful identification and parameter extraction of the information encoded in the waves will provide us with a detailed description of the physics of BHs. A search based on matched-filtering for characterization and parameter extraction requires the development of some $10^{15}$ waveforms. If a third additional BH perturbed the system, the waveforms would not be applicable, and we would need to increase the number of templates required for a valid detection. In this letter, we calculate the probability that more than two BHs interact in the regime of strong relativity in a dense stellar cluster. We determine the physical properties necessary in a stellar system for three black holes to have a close encounter in this regime and also for an existing binary of two BHs to have a strong interaction with a third hole. In both cases the event rate is negligible. While dense stellar systems such as galactic nuclei, globular clusters and nuclear stellar clusters are the breeding grounds for the sources of gravitational waves that ground-based detectors like Advanced LIGO and Advanced VIRGO will be exploring, the analysis of the waveforms in full general relativity needs only to evaluate the two-body problem. This reduces the number of templates of waveforms to create by orders of magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Binaries of massive black holes in rotating clusters: Dynamics, gravitational waves, detection and the role of eccentricity

The dynamical evolution of binaries of intermediate-massive black holes (IMBHs, massive black holes with a mass ranging between $10^2$ and $10^4 M_{\odot}$) in stellar clusters has recently received an increasing amount of attention. This is at least partially due to the fact that if the binary is hard enough to evolve to the phase at which it will start emitting gravitational waves (GWs) efficiently, there is a good probability that it will be detectable by future space-borne detectors like LISA. We study this evolution in the presence of rotation in the cluster. The eccentricity is strongly connected to the initial IMBHs velocities, and values of $\sim 0.7$ up to 0.9 are reached for low initial velocities, while almost circular orbits result if the initial velocities are increased. A Monte Carlo study indicates that these sources will be detectable by a detector such as LISA with median signal to noise ratios of between 10 and 20 over a three year period, although some events had signal to noise ratios of 300 or greater. Furthermore, one should also be able to estimate the chirp-mass with median fractional errors of $10^{-4}$, reduced mass on the order of $10^{-3}$ and luminosity distance on the order of $10^{-1}$. Finally, these sources will have a median angular resolution in the LISA detector of about 3 square degrees, putting events firmly in the field of view of future electromagnetic detectors such as LSST.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Gravitational waves from eccentric intermediate-mass black hole binaries

If binary intermediate-mass black holes (IMBHs; with masses between 100 and $10^4 \Msun$) form in dense stellar clusters, their inspiral will be detectable with the planned Laser Interferometer Space Antenna (LISA) out to several Gpc. Here we present a study of the dynamical evolution of such binaries using a combination of direct $N$-body techniques (when the binaries are well separated) and three-body relativistic scattering experiments (when the binaries are tight enough that interactions with stars occur one at a time). We find that for reasonable IMBH masses there is only a mild effect on the structure of the surrounding cluster even though the binary binding energy can exceed the binding energy of the cluster. We demonstrate that, contrary to standard assumptions, the eccentricity in the LISA band can be in {\em some} cases as large as $\sim 0.2 - 0.3$ and that it induces a measurable phase difference from circular binaries in the last year before merger. We also show that, even though energy input from the binary decreases the density of the core and slows down interactions, the total time to coalescence is short enough (typically less than a hundred million years) that such mergers will be unique snapshots of clustered star formation.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On strong mass segregation around a massive black hole: Implications for lower-frequency gravitational-wave astrophysics

We present, for the first time, a clear $N$-body realization of the {\it strong mass segregation} solution for the stellar distribution around a massive black hole. We compare our $N$-body results with those obtained by solving the orbit-averaged Fokker-Planck (FP) equation in energy space. The $N$-body segregation is slightly stronger than in the FP solution, but both confirm the {\it robustness} of the regime of strong segregation when the number fraction of heavy stars is a (realistically) small fraction of the total population. In view of recent observations revealing a dearth of giant stars in the sub-parsec region of the Milky Way, we show that the time scales associated with cusp re-growth are not longer than $(0.1-0.25) \times T_{rlx}(r_h)$. These time scales are shorter than a Hubble time for black holes masses $\mbul \lesssim 4 \times 10^6 M_\odot$ and we conclude that quasi-steady, mass segregated, stellar cusps may be common around MBHs in this mass range. Since EMRI rates scale as $\mbul^{-α}$, with $α\in [1\4,1]$, a good fraction of these events should originate from strongly segregated stellar cusps.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Triplets of supermassive black holes: Astrophysics, Gravitational Waves and Detection

Supermassive black holes (SMBHs) found in the centers of many galaxies have been recognized to play a fundamental active role in the cosmological structure formation process. In hierarchical formation scenarios, SMBHs are expected to form binaries following the merger of their host galaxies. If these binaries do not coalesce before the merger with a third galaxy, the formation of a black hole triple system is possible. Numerical simulations of the dynamics of triples within galaxy cores exhibit phases of very high eccentricity (as high as $e \sim 0.99$). During these phases, intense bursts of gravitational radiation can be emitted at orbital periapsis. This produces a gravitational wave signal at frequencies substantially higher than the orbital frequency. The likelihood of detection of these bursts with pulsar timing and the Laser Interferometer Space Antenna ({\it LISA}) is estimated using several population models of SMBHs with masses $\gtrsim 10^7 {\rm M_\odot}$. Assuming a fraction of binaries $\ge 0.1$ in triple system, we find that few to few dozens of these bursts will produce residuals $>1$ ns, within the sensitivity range of forthcoming pulsar timing arrays (PTAs). However, most of such bursts will be washed out in the underlying confusion noise produced by all the other 'standard' SMBH binaries emitting in the same frequency window. A detailed data analysis study would be required to assess resolvability of such sources. Implementing a basic resolvability criterion, we find that the chance of catching a resolvable burst at a one nanosecond precision level is 2-50%, depending on the adopted SMBH evolution model. On the other hand, the probability of detecting bursts produced by massive binaries (masses $\gtrsim 10^7\msun$) with {\it LISA} is negligible.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Dense gas-star systems: Super-massive stars evolution

In the 60s and 70s super-massive central objects (from now onwards SMOs) were thought to be the main source of active galactic nuclei (AGNs) characteristics (luminosities of $L \approx 10^{12} L_{\odot}$). The release of gravitational binding energy by the accretion of material on to an SMO in the range of $10^7 - 10^9 M_{\odot}$ has been suggested to be the primary powerhouse (Lynden-Bell 1969). That rather exotic idea in early time has become common sense nowadays. Not only our own galaxy harbours a few million-solar mass black hole (Genzel 2001) but also many of other non-active galaxies show kinematic and gas-dynamic evidence of these objects (Magorrian et al. 1998) The concept of central super-massive stars (SMSs henceforth) (${\cal M} \ge 5 \times 10^4 M_{\odot}$, where ${\cal M}$ is the mass of the SMS) embedded in dense stellar systems was suggested as a possible explanation for high- energy emissions phenomena occurring in AGNs and quasars (Vilkoviski 1976, Hara 1978), such as X-ray emissions (Bahcall and Ostriker, 1975). SMSs and super-massive black holes (SMBHs) are two possibilities to explain the nature of SMOs, and SMSs may be an intermediate step towards the formation of SMBHs (Rees 1984). In this paper we give the equations that describe the dynamics of such a dense star-gas system which are the basis for the code that will be used in a prochain future to simulate this scenario. We also briefly draw the mathematical fundamentals of the code

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

The loss-cone problem in dense nuclei

We address the classical problem of star accretion onto a supermassive central gaseous object in a galactic nucleus. The resulting supermassive central gas-star object is assumed to be located at the centre of a dense stellar system for which we use a simplified model consisting of a Plummer model with an embedded density cusp using stellar point masses. From the number of stars belonging to the loss-cone, which plunge onto the central object on elongated orbits from outside, we estimate the accretion rate taking into account a possible anisotropy of the surrounding stellar distribution. The total heating rate in the supermassive star due to the loss-cone stars plunging onto it is estimated. This semi-analytical study, revisiting and expanding classical paper's work, is a starting point of future work on a more detailed study of early evolutionary phases of galactic nuclei. It merits closer examination, because it is one of the key features for the link between cosmology and galaxy formation.

0 citations0 reviewsNo code linkedOpen access