Source author record

Xian Chen

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

37works
16topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

37 published item(s)

preprint2022arXiv

A Higher Probability of Detecting Lensed Supermassive Black Hole Binaries by LISA

Gravitational lensing of gravitational waves (GWs) is a powerful probe of the matter distribution in the universe. Here we revisit the wave-optics effects induced by dark matter (DM) halos on the GW signals of merging massive black hole binaries (MBHBs), and we study the possibility of discerning these effects using the Laser Interferometer Space Antenna (LISA). In particular, we include the halos in the low-mass range of $\rm 10^5-10^8\, M_\odot$ since they are the most numerous according to the cold DM model. We simulate the lensed signals corresponding to a wide range of impact parameters, and we find distinguishable deviation from the standard best-fit GW templates even when the impact parameter is as large $y\simeq50$. Consequently, we estimate that over $(0.1-1.6)\%$ of the MBHBs in the mass range of $\rm 10^{5.0}-10^{6.5}\, M_\odot$ and the redshift range of $4-10$ should show detectable wave-optics effects. This probability is one order of magnitude higher than that derived in previous works. The uncertainty comes mainly from the mass function of the DM halos. Not detecting any signal during the LISA mission would imply that DM halos with $\rm 10^5-10^8\, M_\odot$ are less numerous than what the cold DM model predicts.

preprint2022arXiv

The Eccentric and Accelerating Stellar Binary Black Hole Mergers in Galactic Nuclei: Observing in Ground and Space Gravitational Wave Observatories

We study the stellar binary black holes (BBHs) inspiralling/merging in galactic nuclei based on our numerical method GNC. We find that $3-40\%$ of all new born BBHs will finally merge due to various dynamical effects. In a five year's mission, up to $10^4$, $10^5$, $\sim100$ of BBHs inspiralling/merging in galactic nuclei can be detected with SNR$>8$ in aLIGO, Einstein/DECIGO, TianQin/LISA/TaiJi, respectively. About tens are detectable in both LISA/TaiJi/TianQin and aLIGO. These BBHs have two unique characteristics: (1) Significant eccentricities. $1-3\%$, $2-7\%$, or $30-90\%$ of them is with $e_i>0.1$ when they enter into aLIGO, Einstein, or space observatories, respectively. Such high eccentricities provide a possible explanation for that of GW 190521. Most highly-eccentric BBHs are not detectable in LISA/Tianqin/TaiJi before entering into aLIGO/Einstein as their strain become significant only at $f_{\rm GW}\gtrsim0.1$ Hz. DECIGO become an ideal observatory to detect those events as it can fully cover the rising phase. (2) Up to $2\%$ of BBHs can inspiral/merge at distances $\lesssim10^3 r_{\rm SW}$ from the massive black hole (MBH), with significant accelerations, such that the Doppler phase drift of $\sim10-10^5$ of them can be detectable with SNR$>8$ in space observatories. The energy density of the gravitational wave backgrounds (GWB) contributed by these BBHs deviate from the powerlaw slope of $2/3$ at $f_{\rm GW}\lesssim 1$mHz. The high eccentricity, significant accelerations and different profile of GWB of these sources make them distinguishable, thus interesting for future GW detections and tests of relativities.

preprint2021arXiv

Excitation of gravitational wave modes by a center-of-mass velocity of the source

Most gravitational wave (GW) sources are moving relative to us. This motion is often closely related to the environment of the source and can thus provide crucial information about the formation of the source and its host. Recently, LIGO and Virgo detected for the first time the subdominant modes of GWs. We show that a motion of the center-of-mass of the source can affect these modes, where the effect is proportional to the velocity of the source. The effect on the GW modes in turn affects the overall frequency of the GW, thus leading to a phase shift. We study the impact of this effect on LIGO/Virgo detections and show that it is detectable for sources with high mass ratios and inclinations. This effect breaks the degeneracy between mass and Doppler shift in GW observations, and opens a new possibility of detecting the motion of a GW source even for constant velocities.

preprint2021arXiv

Milli-Hertz Gravitational Wave Background Produced by Quasi-Periodic Eruptions

Extreme-mass-ratio inspirals (EMRIs) are important targets for future space-borne gravitational-wave (GW) detectors, such as the Laser Interferometer Sapce Antenna (LISA). Recent works suggest that EMRI may reside in a population of newly discovered X-ray transients called "quasi-periodic eruptions" (QPEs). Here we follow this scenario and investigate the detectability of the five recently discovered QPEs by LISA. We consider two specific models in which the QPEs are made of either stellar-mass objects moving on circular orbits around massive black holes (MBHs) or white dwarfs (WDs) on eccentric orbits around MBHs. We find that in either case each QPE is too weak to be resolvable by LISA. However, if QPEs are made of eccentric WD-MBH binaries, they radiate GWs in a wide range of frequencies. The broad spectra overlap to form a background which, between $0.003-0.02$ Hz, exceeds the background known to exist due to other types of sources. Presence of this GW background in the LISA band could impact the future search for the seed black holes at high redshift as well as the stellar-mass binary black holes in the local universe.

preprint2020arXiv

A Shining Death of Unequal Supermassive Black Hole Binaries

In the $Λ$CDM scenario, small galaxies merge to produce larger entities. Since supermassive black holes (SMBHs) are found in galaxies of all sizes, SMBH binaries (SMBHBs) are generally expected to form during the amalgamation of galaxies. It is unclear what fraction of these binaries could eventually merge, but a general consensus is that initially the orbital decay is mediated by the surrounding gas and stars. In this Letter, we show that in active galactic nulcei (AGNs) the radiation field also causes the orbits of the accreting SMBHs to shrink. The corresponding mechanism, known as the "Poynting-Robertson drag" (PR drag), takes effect on a well-defined timescale $CT_{\rm Sal}$, where $T_{\rm Sal}$ is the Salpeter timescale of the AGN, presumably coincide with the primary SMBH, and $C=4ξ^{-1}ε^{-1}q^{1/3}(1+q)^{2/3}(1-ε)$ is a constant determined by the radiative efficiency $ε$, the mass ratio $q$ of the two black holes, and a parameter $ξ$ characterizing the size of the circum-secondary accretion disk. We find that when $q\lesssim$a few$\times10^{-5}$, the PR drag is more efficient in shrinking the binary than many other mechanisms, such as dynamical friction and type-I migration. Our finding points to a possible new channel for the coalescence of unequal SMBHBs and the clearing of intermediate-massive black holes in AGNs.

preprint2020arXiv

Degeneracy between mass and peculiar acceleration for the double white dwarfs in the LISA band

Mass and distance are fundamental quantities to measure in gravitational-wave (GW) astronomy. However, recent studies suggest that the measurement may be biased due to the acceleration of GW source. Here we develop an analytical method to quantify such a bias induced by a tertiary on a double white dwarf (DWD), since DWDs are the most common GW sources in the milli-Hertz band. We show that in a large parameter space the mass is degenerate with the peculiar acceleration, so that from the waveform we can only retrieve a mass of ${\cal M}(1+Γ)^{3/5}$, where ${\cal M}$ is the real chirp mass of the DWD and $Γ$ is a dimensionless factor proportional to the peculiar acceleration. Based on our analytical method, we conduct mock observation of DWDs by the Laser Interferometer Space Antenna (LISA). We find that in about $9\%$ of the cases the measured chirp mass is biased due to the presence of a tertiary by $(5-30)\%$. Even more extreme cases are found in about a dozen DWDs and they may be misclassified as double neutron stars, binary black holes, DWDs undergoing mass transfer, or even binaries containing lower-mass-gap objects and primordial black holes. The bias in mass also affects the measurement of distance, resulting in a seemingly over-density of DWDs within a heliocentric distance of $1$ kpc as well as beyond $100$ kpc. Our result highlights the necessity of modeling the astrophysical environments of GW sources to retrieve their correct physical parameters.

preprint2020arXiv

Fake massive black holes in the milli-Hertz gravitational-wave band

In gravitational wave (GW) astronomy accurate measurement of the source parameters, such as mass, relies on accurate waveform templates. Currently, the templates are developed assuming that the source, such as a binary black hole (BBH), is residing in a vacuum. However, astrophysical models predict that BBHs could form in gaseous environments, such as common envelops, stellar cores, and accretion disks of active galactic nuclei. Here we revisit the impact of gas on the GW waveforms of stellar-mass BBHs with a focus on the early inspiral phase when the GW frequency is around milli-Hertz. We show that for these BBHs, gas friction could dominate the dynamical evolution and hence duplicate chirp signals. The relevant hydrodynamical timescale, $τ_{\rm gas}$, could be much shorter than the GW radiation timescale, $τ_{\rm gw}$, in the above astrophysical scenarios. As a result, the observable chirp mass is higher than the real one by a factor of $(1+τ_{\rm gw}/τ_{\rm gas})^{3/5}$ if the gas effect is ignored in the data analysis. Such an error also results in an overestimation of the source distance by a factor of $(1+τ_{\rm gw}/τ_{\rm gas})$. By performing matched-filtering analysis in the milli-Hertz band, we prove that the gas-dominated signals are practically indistinguishable from the chirp signals of those more massive BBHs residing in a vacuum environment. Such fake massive objects in the milli-Hertz band, if not appropriately accounted for in the future, may alter our understanding of the formation, evolution, and detection of BBHs.

preprint2020arXiv

Impact of an Active Sgr A* on the Synthesis of Water and Organic Molecules Throughout the Milky Way

Sgr A*, the supermassive black hole (SMBH) in our Galaxy, is dormant today, but it should have gone through multiple gas-accretion episodes in the past billions of years to grow to its current mass of $4\times10^6\,M_\odot$. Each episode temporarily ignites the SMBH and turns the Galactic Center into an active galactic nucleus (AGN). Recently, we showed that the AGN could produce large amount of hard X-rays that can penetrate the dense interstellar medium in the Galactic plane. Here we further study the impact of the X-rays on the molecular chemistry in our Galaxy. We use a chemical reaction network to simulate the evolution of several molecular species including H2O, CH3OH, and H2CO, both in the gas phase and on the surface of dust grains. We find that the X-ray irradiation could significantly enhance the abundances of these species. The effect is the most significant in those young, high-density molecular clouds, and could be prominent at a Galactic distance of $8$ kpc or smaller. The imprint in the chemical abundance is visible even several million years after the AGN turns off.

preprint2020arXiv

Impact of Leakage for Electricity Generation by Pyroelectric Converter

Pyroelectric energy converter is a functional capacitor using pyroelectric material as the dielectric layer. Utilizing the first-order phase transformation of the material, the pyroelectric device can generate adequate electricity within small temperature fluctuations. However, most pyroelectric capacitors are leaking during energy conversion. In this paper, we analyze the thermodynamics of pyroelectric energy conversion with consideration of the electric leakage. Our thermodynamic model is verified by experiments using three phase-transforming ferroelectric materials with different pyroelectric properties and leakage behaviors. We demonstrate that the impact of leakage for electric generation is prominent, and sometimes may be confused with the actual power generation by pyroelectricity. We discover an ideal material candidate, (Ba,Ca)(Ti,Zr,Ce)O$_3$, which exhibits large pyroelectric current and extremely low leakage current. The pyroelectric converter made of this material generates 1.95 $μ$A/cm$^2$ pyroelectric current density and 0.2 J/cm$^3$ pyroelectric work density even after 1389 thermodynamic conversion cycles.

preprint2020arXiv

Measuring black hole masses from tidal disruption events and testing the $M_{\rm BH}-σ_*$ relation

Liu and collaborators recently proposed an elliptical accretion disk model for tidal disruption events (TDEs). They showed that the accretion disks of optical/UV TDEs are large and highly eccentric and suggested that the broad optical emission lines with complex and diverse profiles originate in the cool eccentric accretion disk of random inclination and orientation. In this paper, we calculate the radiation efficiency of the elliptical accretion disk and investigate the implications for the observations of TDEs. We compile observational data for the peak bolometric luminosity and total radiation energy after peak brightness of 18 TDE sources and compare these data to the predictions from the elliptical accretion disk model. Our results show that the observations are consistent with the theoretical predictions and that the majority of the orbital energy of the stellar debris is advected into the black hole (BH) without being converted into radiation. Furthermore, we derive the masses of the disrupted stars and the masses of the BHs of the TDEs. The BH masses obtained in this paper are also consistent with those calculated with the $M_{\rm BH} - σ_*$ relation. Our results provide an effective method for measuring the masses of BHs in large numbers of TDEs to be discovered in ongoing and next-generation sky surveys, regardless of whether the BHs are located at the centers of galactic nuclei or wander in disks and halos.

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.

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.

preprint2019arXiv

Impact of a Spinning Supermassive Black Hole on the Orbit and Gravitational Waves of a Nearby Compact Binary

Recent theoretical studies suggest that stellar-mass binary black holes (BBHs) would merge more efficiently due to the Kozai-Lidov mechanism if these binaries form in the vicinity of supermassive black holes (SMBHs). Since SMBHs are likely rotating rapidly, we continue our earlier study on the generalization of the Kozai-Lidov formalism to include the spin of the SMBH and study the evolution of a nearby BBH. We find that the eccentricity and orbital inclination of the BBH is significantly affected, because the spin (i) forces the orbital plane of the center-of-mass of the BBH around the SMBH to precess (the Lense-Thirring effect) and (ii) imposes an additional gravitomagnetic force on the BBH. As a result, the merger time of the BBH could be significantly different. We calculate the waveform from the BBH in one representative example and study its detectability by a milli-Hertz GW detector, such as the Laser Interferometer Space Antenna (LISA). We find that the signal is distinguishable from that in the case without spin. Our results imply that the BBHs in the LISA band could potentially be used to probe the spin of the SMBHs in galaxy centers.

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.

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.

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.

preprint2016arXiv

SDSS J0159+0105: A Radio-Quiet Quasar with a Centi-Parsec Supermassive Black Hole Binary Candidate

We report a candidate centi-parsec supermassive black hole binary (SMBHB) in the radio-quiet quasar SDSS J0159+0105 at z=0.217. With a modified lomb-scargle code GLSdeDRW and the auto-correlation analysis ACF, we detect two significant (at P>99%) periodic signals at ~741 day and ~1500 day from the 8.1-year Catalina V-band light curve of this quasar. The period ratio, which is close to 1:2, is typical of a black-hole binary system with a mass ratio of 0.05<q<0.8 according to recent numerical simulations. SDSS J0159+0105 has two SDSS spectroscopic observations separated by ~10 years. There is a significant change in the broad H-beta profile between the two epochs, which can be explained by a single broad-line region (BLR) around the binary system illuminated by the aforementioned mini-disks, or a stream of gas flowing from the circumbinary disk to one of the SMBHs. From the single BLR assumption and the orbital period t_orb ~1500 day, we estimate the total virial masses of M_SMBHB ~ 1.3x10^8 M_sun, the average distances of BLR of ~0.04pc (~50 light-day, with +/-0.3 dex uncertainty), and a SMBHB separation of d= (0.01pc)M_{8,tot}^{1/3} (T_rest/3.3yr)^{2/3} ~ 0.013 pc (15 light-day). Based on analytical work, the postulated circumbinary disk has an inner radius of 2d = 0.026 pc (30 light-day). SDSS J0159+0105 also displays unusual spectral energy distribution. The unique properties of SDSS J0159+0105 are consistent with it being a centi-parsec SMBHB. The GLSdeDRW code link: http://butler.lab.asu.edu/qso_period/

preprint2016arXiv

The Astropy Problem

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

preprint2016arXiv

The TDE ASASSN-14li and its host resolved at parsec scales with the EVN

We report European Very Long Baseline Interferometry Network (EVN) radio continuum observations of ASASSN-14li, one of the best studied tidal disruption events (TDEs) to date. At 1.7 GHz with ~12x6mas resolution, the emission is unresolved. At 5.0 GHz with ~3x2mas resolution, the radio emission shows an extended structure that can be modeled with two components: a core-like component and a fainter, possibly elongated source 4.3mas (~2pc) away. Our observations are not conclusive as to the nature of the components, but three scenarios are discussed. One possibility is a core-jet/outflow morphology, thus making of ASASSN-14li the first TDE jet/outflow directly imaged. For this case, the projected separation between the two components can only be explained by superluminal motion, rather than the lower velocities inferred from low-resolution radio observations. However, typical fast moving jets have brightness temperatures ~5 orders of magnitude higher than we find, thus making this scenario less likely. The second possibility is that we are imaging a non-relativistic jet from past AGN/TDE activity. In this case a past TDE is preferred given that the spatial extension and radio luminosity of the elongated component are consistent with the theoretical predictions for a TDE outflow. Alternatively, the two sources could indicate the presence of a binary black hole, which would then naturally explain the enhanced TDE rates of post-starburst galaxies. Future EVN observations will help us to distinguish between these scenarios.

preprint2015arXiv

Continuous-time Markov decision processes under the risk-sensitive average cost criterion

This paper studies continuous-time Markov decision processes under the risk-sensitive average cost criterion. The state space is a finite set, the action space is a Borel space, the cost and transition rates are bounded, and the risk-sensitivity coefficient can take arbitrary positive real numbers. Under the mild conditions, we develop a new approach to establish the existence of a solution to the risk-sensitive average cost optimality equation and obtain the existence of an optimal deterministic stationary policy.

preprint2015arXiv

Determination of the stretch tensor for structural transformations

The transformation stretch tensor plays an essential role in the evaluation of conditions of compatibility between phases and the use of the Cauchy-Born rule. This tensor is difficult to measure directly from experiment. We give an algorithm for the determination of the transformation stretch tensor from x-ray measurements of structure and lattice parameters. When evaluated on some traditional and emerging phase transformations the algorithm gives unexpected results.

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.

preprint2015arXiv

The Gamma-ray Afterglows of Tidal Disruption Events

A star wandering too close to a supermassive black hole (SMBH) will be tidally disrupted. Previous studies of such "tidal disruption event" (TDE) mostly focus on the stellar debris that are bound to the system, because they give rise to luminous flares. On the other hand, half of the stellar debris in principle are unbound and can stream to a great distance, but so far there is no clear evidence that this "unbound debris stream" (UDS) exists. Motivated by the fact that the circum-nuclear region around SMBHs is usually filled with dense molecular clouds (MCs), here we investigate the observational signatures resulting from the collision between an UDS and a MC, which is likely to happen hundreds of years after a TDE. We focus on $γ$-ray emission ($0.1-10^5$ GeV), which comes from the encounter of shock-accelerated cosmic rays with background protons and, more importantly, is not subject to extinction. We show that because of the high proton density inside a MC, the peak $γ$-ray luminosity, about $10^{39}~{\rm erg~s^{-1}}$, is at least $100$ times greater than that in the case without a MC (only with a smooth interstellar medium). The luminosity decays on a time-scale of decades, depending on the distance of the MC, and about a dozen of these "TDE afterglows" could be detected within a distance of about $16$ Mpc by the future Cherenkov Telescope Array. Without careful discrimination, these sources potentially could contaminate the searches for starburt galaxies, galactic nuclei containing millisecond pulsars or dark-matter annihilation signals.

preprint2015arXiv

The Midpoint Mixed Model with a Missingness Mechanism (M5): A Likelihood-Based Framework for Quantification of Mass Spectrometry Proteomics Data (Preprint)

Statistical models for proteomics data often estimate protein fold changes between two samples, A and B, as the average peptide intensity from sample A divided by the average peptide intensity from sample B. Such average intensity ratios fail to take full advantage of the experimental design which eliminates unseen confounding variables by processing peptides from both samples under identical conditions. Typically this structure is exploited through the estimation of a protein ratio as the median ratio of matched peptide intensities. This simple solution fails to account for a substantial missing data bias which has led to the development of more sophisticated average intensity models. Here we develop the first statistical model that accounts for nonignorable missingness while utilizing peptide level matched pairs across samples. Our simulation analysis shows that models which fail to utilize peptide level ratios, su.er astonishing losses to accuracy with basic ANOVA estimates having an average MSE 371% higher than median ratio estimates. In turn, median ratio estimates have an average MSE 35% higher than our model estimates. An analysis of breast cancer data reinforces these relationships and shows that our model is capable of increasing the number of proteins estimated by 22%.

preprint2015arXiv

Unbound Debris Streams and Remnants Resulting From the Tidal Disruptions of Stars by Supermassive Black Holes

The kinetic energy of a star in orbit about a supermassive black hole is a significant fraction of its rest mass energy when its periapse is comparable to its tidal radius. Upon its destruction, a fraction of this energy is extracted and injected into the stellar debris, half of which becomes unbound from the black hole, with the fastest material moving at $\sim 0.03 c$. In this paper, we present a formalism for determining the fate of these unbound debris streams (UDSs) as they depart from the black hole and interact with the surrounding gas. As the density and velocity varies along the length of a UDS, we find that hydrodynamical drag quickly shapes UDSs into loop-like structures, with the densest portions of the streams leading portions of lower density. As UDSs travel outwards, their drag against the ISM increases quadratically with distance, which causes UDSs to deposit their momentum and energy into the ambient medium before the surrounding shocked ISM has a chance to cool. This sudden injection of $\sim 10^{50}$ erg into the ambient medium generates a Sedov-like unbound debris remnant (UDR) that mimics supernova remnants (SNRs) in energetics and appearance, accelerates particles which will produce cosmic rays and synchrotron emission, and provides momentum feedback into the molecular clouds surrounding a black hole. We estimate that a few of these UDRs might be present within a couple degrees of the Galactic Center masquerading as SNRs, and that the UDR scenario is a plausible explanation for Sgr A East.

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.

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.

preprint2014arXiv

Markov jump processes in modeling coalescent with recombination

Genetic recombination is one of the most important mechanisms that can generate and maintain diversity, and recombination information plays an important role in population genetic studies. However, the phenomenon of recombination is extremely complex, and hence simulation methods are indispensable in the statistical inference of recombination. So far there are mainly two classes of simulation models practically in wide use: back-in-time models and spatially moving models. However, the statistical properties shared by the two classes of simulation models have not yet been theoretically studied. Based on our joint research with CAS-MPG Partner Institute for Computational Biology and with Beijing Jiaotong University, in this paper we provide for the first time a rigorous argument that the statistical properties of the two classes of simulation models are identical. That is, they share the same probability distribution on the space of ancestral recombination graphs (ARGs). As a consequence, our study provides a unified interpretation for the algorithms of simulating coalescent with recombination, and will facilitate the study of statistical inference on recombination.

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$^*$.

preprint2013arXiv

Enhanced off-center stellar tidal disruptions by supermassive black holes in merging galaxies

Off-center stellar tidal disruption flares have been suggested to be a powerful probe of recoiling supermassive black holes (SMBHs) out of galactic centers due to anisotropic gravitational wave radiations. However, off-center tidal flares can also be produced by SMBHs in merging galaxies. In this paper, we computed the tidal flare rates by dual SMBHs in two merging galaxies before the SMBHs become self-gravitationally bounded. We employ an analytical model to calculate the tidal loss-cone feeding rates for both SMBHs, taking into account two-body relaxation of stars, tidal perturbations by the companion galaxy, and chaotic stellar orbits in triaxial gravitational potential. We show that for typical SMBHs with mass 10^7 M_\sun, the loss-cone feeding rates are enhanced by mergers up to Γ~ 10^{-2} yr^{-1}, about two order of magnitude higher than those by single SMBHs in isolated galaxies and about four orders of magnitude higher than those by recoiling SMBHs. The enhancements are mainly due to tidal perturbations by the companion galaxy. We suggest that off-center tidal flares are overwhelmed by those from merging galaxies, making the identification of recoiling SMBHs challenging. Based on the calculated rates, we estimate the relative contributions of tidal flare events by single, binary, and dual SMBH systems during cosmic time. Our calculations show that the off-center tidal disruption flares by un-bound SMBHs in merging galaxies contribute a fraction comparable to that by single SMBHs in isolated galaxies. We conclude that off-center tidal disruptions are powerful tracers of the merging history of galaxies and SMBHs.

preprint2013arXiv

Study of the cofactor conditions: conditions of supercompatibility between phases

The cofactor conditions, introduced in James and Zhang, are conditions of compatibility between phases in martensitic materials. They consist of three subconditions: i) the condition that the middle principal stretch of the transformation stretch tensor $\mathbf U$ is unity ($λ_2 = 1$), ii) the condition $\mathbf a \cdot \mathbf U\, \cof (\mathbf U^2 - \mathbf I)\mathbf n = 0$, where the vectors $\mathbf a$ and $\mathbf n$ are certain vectors arising in the specification of the twin system, and iii) the inequality ${\rm tr} \mathbf U^2 + \det \mathbf U^2 -(1/4) |\mathbf a|^2 |\mathbf n|^2\ge 2$. Together, these conditions are necessary and sufficient for the equations of the crystallographic theory of martensite to be satisfied for the given twin system but for any volume fraction f of the twins, $0 \le f \le 1$. This contrasts sharply with the generic solutions of the crystallographic theory which have at most two such volume fractions for a given twin system of the form $f^*$ and $1-f^*$. In this paper we simplify the form of the cofactor conditions, we give their specific forms for various symmetries and twin types, we clarify the extent to which the satisfaction of the cofactor conditions for one twin system implies its satisfaction for other twin systems. In particular, we prove that the satisfaction of the cofactor conditions for either Type I or Type II twins implies that there are solutions of the crystallographic theory using these twins that have no elastic transition layer. We show that the latter further implies macroscopically curved, transition-layer-free austenite/martensite interfaces for Type I twins, and planar transition-layer-free interfaces for Type II twins which nevertheless permit significant flexibility of the martensite. We identify some real material systems nearly satisfying the cofactor conditions.

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.

preprint2012arXiv

Interaction of Recoiling Supermassive Black Holes with Stars in Galactic Nuclei

Supermassive black hole binaries (SMBHBs) are the products of frequent galaxy mergers. The coalescence of the SMBHBs is a distinct source of gravitational wave (GW) radiation. The detections of the strong GW radiation and their possible electromagnetic counterparts are essential. Numerical relativity suggests that the post-merger supermassive black hole (SMBH) gets a kick velocity up to 4000 km/s due to the anisotropic GW radiations. Here we investigate the dynamical co-evolution and interaction of the recoiling SMBHs and their galactic stellar environments with one million direct N-body simulations including the stellar tidal disruption by the recoiling SMBHs. Our results show that the accretion of disrupted stars does not significantly affect the SMBH dynamical evolution. We investigate the stellar tidal disruption rates as a function of the dynamical evolution of oscillating SMBHs in the galactic nuclei. Our simulations show that most of stellar tidal disruptions are contributed by the unbound stars and occur when the oscillating SMBHs pass through the galactic center. The averaged disruption rate is ~10^{-6} M_\odot yr^{-1}, which is about an order of magnitude lower than that by a stationary SMBH at similar galactic nuclei. Our results also show that a bound star cluster is around the oscillating SMBH of about ~ 0.7% the black hole mass. In addition, we discover a massive cloud of unbound stars following the oscillating SMBH. We also investigate the dependence of the results on the SMBH masses and density slopes of the galactic nuclei.

preprint2012arXiv

Is there an intermediate massive black hole in the Galactic center: Imprints on the stellar tidal-disruption rate

It has been suggested that an intermediate-massive black hole (IMBH) with mass $10^{3-5} M_\odot$ could fall into the galactic center (GC) and form a massive black hole binary (MBHB) with the central supermassive black hole, but current observations are not sensitive to constrain all mass and distance ranges. Motivated by the recent discovery that MBHBs could enhance the rate of tidal-disruption events (TDEs) of stellar objects, we investigate the prospect of using stellar-disruption rate to probe IMBHs in the GC. We incorporated the perturbation by an IMBH into the loss-cone theory and calculated the stellar-disruption rates in the GC. We found that an IMBH heavier than $2000 M_\odot$ could distinguishably enhance the stellar-disruption rate. By comparing observations of Sgr A* with the fall-back model for stellar debris, we suggested that the TDE rate in our Galaxy should not significantly exceed $0.002 yr^{-1}$, therefore a fraction of the parameter space for the IMBH, concentrating at the high-mass end, can already be excluded. To derive constraint in the remaining parameter space, it is crucial to observationally confirm or reject the stellar-disruption rate between $10^{-4}$ and $10^{-2} yr^{-1}$, and we discussed possible strategies to make such measurements.

preprint2011arXiv

Recoiling Supermassive Black Holes in Spin-flip Radio Galaxies

Numerical relativity simulations predict that coalescence of supermassive black hole (SMBH) binaries not only leads to a spin flip but also to a recoiling of the merger remnant SMBHs. In the literature, X-shaped radio sources are popularly suggested to be candidates for SMBH mergers with spin flip of jet-ejecting SMBHs. Here we investigate the spectral and spatial observational signatures of the recoiling SMBHs in radio sources undergoing black hole spin flip. Our results show that SMBHs in most spin-flip radio sources have mass ratio $q\ga 0.3$ with a minimum possible value $q_{\rm min} \simeq 0.05$. For major mergers, the remnant SMBHs can get a kick velocity as high as $2100 km s^{-1}$ in the direction within an angle $\la 40^\circ$ relative to the spin axes of remnant SMBHs, implying that recoiling quasars are biased to be with high Doppler-shifted broad emission lines while recoiling radio galaxies are biased to large apparent spatial off-center displacements. We also calculate the distribution functions of line-of-sight velocity and apparent spatial off-center for spin-flip radio sources with different apparent jet reorientation angles. Our results show that the larger the apparent jet reorientation angle is, the larger the Doppler-shifting recoiling velocity and apparent spatial off-center displacement will be. We investigate the effects of recoiling velocity on the dust torus in spin-flip radio sources and suggest that recoiling of SMBHs would lead to "dust poor" AGNs. Finally, we collect a sample of 19 X-shaped radio objects and for each object give the probability of detecting the predicted signatures of recoiling SMBH.

preprint2010arXiv

Tidal stellar disruptions by massive black hole pairs: II. Decaying binaries

Tidal stellar disruptions have traditionally been discussed as a probe of the single, massive black holes (MBHs) that are dormant in the nuclei of galaxies. In Chen et al. (2009), we used numerical scattering experiments to show that three-body interactions between bound stars in a stellar cusp and a non-evolving "hard" MBH binary will also produce a burst of tidal disruptions, caused by a combination of the secular "Kozai effect" and by close resonant encounters with the secondary hole. Here we derive basic analytical scalings of the stellar disruption rates with the system parameters, assess the relative importance of the Kozai and resonant encounter mechanisms as a function of time, discuss the impact of general relativistic (GR) and extended stellar cusp effects, and develop a hybrid model to self-consistently follow the shrinking of an MBH binary in a stellar background, including slingshot ejections and tidal disruptions. In the case of a fiducial binary with primary hole mass M_1=10^7\msun and mass ratio q=M_2/M_1=1/81, embedded in an isothermal cusp, we derive a stellar disruption rate \dot{N_*} ~ 0.2/yr lasting ~ 3X10^5 yr. This rate is 3 orders of magnitude larger than the corresponding value for a single MBH fed by two-body relaxation, confirming our previous findings. For q<<0.01, the Kozai/chaotic effect could be quenched due to GR/cusp effects by an order of magnitude, but even in this case the stellar-disruption rate is still two orders of magnitude larger than that given by standard relaxation processes around a single MBH. Our results suggest that >~10% of the tidal-disruption events may originate in MBH binaries.

preprint2009arXiv

Interruption of Tidal Disruption Flares By Supermassive Black Hole Binaries

Supermassive black hole binaries (SMBHBs) are products of galaxy mergers, and are important in testing Lambda cold dark matter cosmology and locating gravitational-wave-radiation sources. A unique electromagnetic signature of SMBHBs in galactic nuclei is essential in identifying the binaries in observations from the IR band through optical to X-ray. Recently, the flares in optical, UV, and X-ray caused by supermassive black holes (SMBHs) tidally disrupting nearby stars have been successfully used to observationally probe single SMBHs in normal galaxies. In this Letter, we investigate the accretion of the gaseous debris of a tidally disrupted star by a SMBHB. Using both stability analysis of three-body systems and numerical scattering experiments, we show that the accretion of stellar debris gas, which initially decays with time $\propto t^{-5/3}$, would stop at a time $T_{\rm tr} \simeq ηT_{\rm b}$. Here, $η\sim0.25$ and $T_{\rm b}$ is the orbital period of the SMBHB. After a period of interruption, the accretion recurs discretely at time $T_{\rm r} \simeq ξT_b$, where $ξ\sim 1$. Both $η$ and $ξ$ sensitively depend on the orbital parameters of the tidally disrupted star at the tidal radius and the orbit eccentricity of SMBHB. The interrupted accretion of the stellar debris gas gives rise to an interrupted tidal flare, which could be used to identify SMBHBs in non-active galaxies in the upcoming transient surveys.