Researcher profile

Mouyuan Sun

Mouyuan Sun contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

13 published item(s)

preprint2026arXiv

Light curve modeling of the nearest neutron star candidate LAMOST J235456.73+335625.9

The discovery of heavy radioactive elements (e.g., $^{60}\mathrm{Fe}$) on Earth suggests that supernova explosions may have occurred near our planet within the past million years, potentially having a significant impact on the ecological environment. This finding has motivated the search for nearby neutron stars in the Solar neighborhood. In a recent study, a candidate for one of the closest neutron stars to Earth, LAMOST J235456.73+335625.9 (hereafter J2354), was reported. Based on dynamical mass measurements under different inclination angle assumptions, the inferred mass range for the unseen compact companion in the system is $1.4$--$1.6$ $M_{\odot}$. Hence, the unseen companion in J2354 is either a massive cold white dwarf or a neutron star. Here we model the flux variations of J2354 as a combination of ellipsoidal modulation and surface spots. We test both cold spot and hot spot models, setting the number of spots to two in each case, and constrain the spot properties through light curve fitting. In the cold spot scenario, the spots are mostly visible at phases $0.5$--$0.75$, whereas in the hot spot scenario, the spots appear predominantly at phases $0.25$--$0.5$. The hot spot model shows better agreement with the observed H$α$ phase variation than the cold spot model. Furthermore, the thermal radiation of a massive but cold white dwarf cannot produce the level of localized heating required to explain the hot spot unless additional heating mechanisms are involved; in contrast, a neutron star can naturally provide such heating through energetic winds. Our results are consistent with the neutron star interpretation of the compact object in J2354.

preprint2022arXiv

Evidence for quasar fast outflows being accelerated at the scale of tens of parsecs

Quasar outflows may play a crucial role in regulating the host galaxy, although the spatial scale of quasar outflows remain a major enigma, with their acceleration mechanism poorly understood. The kinematic information of outflow is the key to understanding its origin and acceleration mechanism. Here, we report the galactocentric distances of different outflow components for both a sample and an individual quasar. We find that the outflow distance increases with velocity, with a typical value from several parsecs to more than one hundred parsecs, providing direct evidence for an acceleration happening at a scale of the order of 10 parsecs. These outflows carry ~1% of the total quasar energy, while their kinematics are consistent with a dust driven model with a launching radius comparable to the scale of a dusty torus, indicating that the coupling between dust and quasar radiation may produce powerful feedback that is crucial to galaxy evolution.

preprint2022arXiv

Nucleosynthesis Contribution of Neutrino-dominated Accretion Flows to the Chemical Evolution of Active Galactic Nuclei

Recent observations of quasars show high line-flux ratios in their broad emission lines and the ratios appear to be independent of redshift up to $z \gtrsim 6$, which indicate that the broad-line regions of these early quasars are surprisingly metal-rich. Here, we revisit the chemical evolution of high-redshift quasars by adding a new ingredient, i.e., the neutrino-dominated accretion flows (NDAFs) with outflows, on top of the conventional core-collapse supernovae (CCSNe). In the presence of the chemical contribution from NDAFs with outflows, the total metal mass (i.e., the summation of the conventional CCSN and NDAFs with outflows) per CCSN depends weakly upon the mass of the progenitor star if the mass is in the range of $\sim 25-55~M_{\odot}$. We model the chemical evolution by adopting a improved open-box model with three typical initial mass functions (IMFs). We find that, with the additional chemical contribution from NDAFs with outflows, the quasar metallicity can be enriched more rapidly in the very early Universe ($z \sim 10$) and reaches higher saturation than the no-NDAF case at $z \sim 8$, after which they evolve slowly with redshift. The quasar metallicity can reach $\sim 20~Z_{\odot}$ ($Z_\odot$ denotes the metallicity of the Sun; and $\sim 20\%$ of which is produced by NDAF outflows) at $z \sim 8$ for the ``top-heavy'' IMF model in \citet{Toyouchi2022}, which readily explains the quasar observations on the super-solar metal abundance and redshift-independent evolution.

preprint2022arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: UV-Optical Accretion Disk Measurements with Hubble Space Telescope

We present accretion-disk structure measurements from UV-optical reverberation mapping observations of a sample of eight quasars at 0.24<z<0.85. Ultraviolet photometry comes from two cycles of Hubble Space Telescope monitoring, accompanied by multi-band optical monitoring by the Las Cumbres Observatory network and Liverpool Telescopes. The targets were selected from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project sample with reliable black-hole mass measurements from Hbeta reverberation mapping results. We measure significant lags between the UV and various optical griz bands using JAVELIN and CREAM methods. We use the significant lag results from both methods to fit the accretion-disk structure using a Markov chain Monte Carlo approach. We study the accretion disk as a function of disk normalization, temperature scaling, and efficiency. We find direct evidence for diffuse nebular emission from Balmer and FeII lines over discrete wavelength ranges. We also find that our best-fit disk color profile is broadly consistent with the Shakura \& Sunyaev disk model. We compare our UV-optical lags to the disk sizes inferred from optical-optical lags of the same quasars and find that our results are consistent with these quasars being drawn from a limited high-lag subset of the broader population. Our results are therefore broadly consistent with models that suggest longer disk lags in a subset of quasars, for example, due to a nonzero size of the ionizing corona and/or magnetic heating contributing to the disk response.

preprint2022arXiv

Thermal Equilibrium Solutions of Black Hole Accretion Flows: Outflows VS Advection

Observations and numerical simulations have shown that outflows generally exist in the accretion process. We revisit the thermal equilibrium solutions of black hole accretion flows by including the role of outflows. Our study focuses on the comparison of the cooling rate of outflows with that of advection. Our results show that, except for the inner region, outflows can dominate over advection in a wide range of the flow, which is in good agreement with previous numerical simulations. We argue that an advection-dominated inner region together with an outflow-dominated outer region should be a general radial distribution for both super-Eddington accretion flows and optically thin flows with low accretion rates.

preprint2022arXiv

X-ray emission of contact binary variables within 1 kpc

By assembling the largest sample to date of X-ray emitting EW-type binaries (EWXs), we carried out correlation analyses for the X-ray luminosity log$L_{\textrm{X}}$, and X-ray activity level log($L_{\textrm{X}}$/$L_{\textrm{bol}}$) versus the orbital period $P$ and effective temperature $T_{\rm eff}$. We find strong $P$-log$L_{\textrm{X}}$ and $P$-log($L_{\textrm{X}}$/$L_{\textrm{bol}}$) correlations for EWXs with $P$ < 0.44 days and we provide the linear parametrizations for these relations, on the basis of which the orbital period can be treated as a good predictor for log$L_{\textrm{X}}$ and log($L_{\textrm{X}}$/$L_{\textrm{bol}}$). The aforementioned binary stellar parameters are all correlated with log$L_{\textrm{X}}$, while only $T_{\rm eff}$ exhibits a strong correlation with log($L_{\textrm{X}}$/$L_{\textrm{bol}}$). Then, EWXs with higher temperature show lower X-ray activity level, which could indicate the thinning of the convective area related to the magnetic dynamo mechanism. The total X-ray luminosity of an EWX is essentially consistent with that of an X-ray saturated main sequence star with the same mass as its primary, which may imply that the primary star dominates the X-ray emission. The monotonically decreasing $P$-log($L_{\textrm{X}}$/$L_{\textrm{bol}}$) relation and the short orbital periods indicate that EWXs could all be in the X-ray saturated state, and they may inherit the changing trend of the saturated X-ray luminosities along with the mass shown by single stars. For EWXs, the orbital period, mass, and effective temperature increase in concordance. We demonstrate that the period $P=0.44$ days corresponds to the primary mass of $\sim1.1 \rm M_\odot$, beyond which the saturated X-ray luminosity of single stars will not continue to increase with mass. This explains the break in the positive $P$-log$L_{\textrm{X}}$ relation for EWXs with $P>0.44$ days.

preprint2021arXiv

An Extraordinary Response of Iron Emission to the Central Outburst in a Tidal Disruption Event Candidate

Understanding the origin of \feii\ emission is important because it is crucial to construct the main sequence of Active Galactic Nuclei (AGNs). Despite several decades of observational and theoretical effort, the location of the optical iron emitting region and the mechanism responsible for the positive correlation between the \feii\ strength and the black hole accretion rate remain open questions as yet. In this letter, we report the optical \feii\ response to the central outburst in PS1-10adi, a candidate tidal disruption event (TDE) taking place in an AGN at $z = 0.203$ that has aroused extensive attention. For the first time, we observe that the \feii\ response in the rising phase of its central luminosity is significantly more prominent than that in the decline phase, showing a hysteresis effect. We interpret this hysteresis effect as a consequence of the gradual sublimation of the dust grains situating at the inner surface of the torus into gas when the luminosity of the central engine increases. It is the iron element released from the sublimated dust that contributes evidently to the observed \feii\ emission. This interpretation, together with the weak response of the \hb\ emission as we observe, naturally explains the applicability of relative \feii\ strength as a tracer of the Eddington ratio. In addition, optical iron emission of this origin renders the \feii\ time lag a potential &#34;standard candle&#34; with cosmological implications.

preprint2021arXiv

The Disk Veiling Effect of the Black Hole Low-Mass X-ray Binary A0620-00

The optical light curves of quiescent black hole low-mass X-ray binaries often exhibit significant non-ellipsoidal variabilities, showing the photospheric radiation of the companion star is veiled by other source of optical emission. Assessing this &#34;veiling&#34; effect is critical to the black hole mass measurement. Here in this work, we carry out a strictly simultaneous spectroscopic and photometric campaign on the prototype of black hole low-mass X-ray binary A0620-00. We find that for each observation epoch, the extra optical flux beyond a pure ellipsoidal modulation is positively correlated with the fraction of veiling emission, indicating the accretion disk contributes most of the non-ellipsoidal variations. Meanwhile, we also obtain a K2V spectral classification of the companion, as well as the measurements of the companion&#39;s rotational velocity $v \sin i = 83.8\pm1.9$ km s$^{-1}$ and the mass ratio between the companion and the black hole $q=0.063\pm0.004$.

preprint2020arXiv

AGNs are not that cool: revisiting the intrinsic AGN far-infrared spectral energy distribution

We investigate the intrinsic spectral energy distribution (SED) of active galactic nuclei (AGNs) at infrared (IR) bands with 42 $z < 0.5$ optically luminous Palomar Green survey quasars through SED decomposition. We decompose the SEDs of the 42 quasars by combining an AGN IR template library Siebenmorgen2015 that covers a wide range of the AGN parameter space with three commonly used galaxy template libraries. We determine the median AGN SED from the best-fitting results. The far-IR (FIR) contribution of our median AGN SED is significantly smaller than that of Symeonidis et al. 2016, but roughly consistent with that of Lyu et al. 2017. The AGN IR SED becomes cooler with increasing bolometric luminosity, which might be due to that more luminous AGNs might have stronger radiative feedback to change torus structures and/or their tori might have higher metallicities. Our conclusions do not depend on the choice of galaxy template libraries. However, since the predicted polycyclic aromatic hydrocarbon (PAH) emission line flux is galaxy template-dependent, cautions should be taken on deriving galaxy FIR contribution from PAH fluxes.

preprint2020arXiv

Corona-Heated Accretion-disk Reprocessing (CHAR): A Physical Model to Decipher the Melody of AGN UV/optical Twinkling

Active galactic nuclei (AGNs) have long been observed to &#34;twinkle&#34; (i.e., their brightness varies with time) on timescales from days to years in the UV/optical bands. Such AGN UV/optical variability is essential for probing the physics of supermassive black holes (SMBHs), the accretion disk, and the broad-line region. Here we show that the temperature fluctuations of an AGN accretion disk, which is magnetically coupled with the corona, can account for observed high-quality AGN optical light curves. We calculate the temperature fluctuations by considering the gas physics of the accreted matter near the SMBH. We find that the resulting simulated AGN UV/optical light curves share the same statistical properties as the observed ones as long as the dimensionless viscosity parameter $α$, which is widely believed to be controlled by magnetohydrodynamic (MHD) turbulence in the accretion disk, is about $0.01$---$0.2$. Moreover, our model can simultaneously explain the larger-than-expected accretion disk sizes and the dependence of UV/optical variability upon wavelength for NGC 5548. Our model also has the potential to explain some other observational facts of AGN UV/optical variability, including the timescale-dependent bluer-when-brighter color variability and the dependence of UV/optical variability on AGN luminosity and black hole mass. Our results also demonstrate a promising way to infer the black-hole mass, the accretion rate, and the radiative efficiency, thereby facilitating understanding of the gas physics and MHD turbulence near the SMBH and its cosmic mass growth history by fitting the AGN UV/optical light curves in the era of time-domain astronomy.

preprint2020arXiv

EUCLIA. II. On the puzzling large UV to X-ray lags in Seyfert galaxies

Recent intense X-ray and UV monitoring campaigns with Swift have detected clear UV lags behind X-ray in several local active galactic nuclei (AGNs). The UV to X-ray lags are often larger (by a factor up to ~ 20) than expected if the UV variation is simply due to the X-ray reprocessing. We previously developed a model in which the UV/optical variations are attributed to disk turbulences, and the effect of large-scale turbulence is considered. Our model, which overcomes many severe challenges to the reprocessing scheme, can well explain the observed variations in NGC 5548, particularly the correlations and lags among the UV/optical bands. In this work, assuming the corona heating is associated with turbulences in the inner accretion disk, we extend our study to model the correlations and lags between the X-ray and UV/optical bands. We find that our model, without the need of light echoing, can well reproduce the observed UV to X-ray lags and the optical to UV lags simultaneously in four local Seyfert galaxies, including NGC 4151, NGC 4395, NGC 4593, and NGC 5548. In our scenario, relatively larger UV to X-ray lag is expected for AGN with smaller innermost disk radius and thus more compact corona. Interestingly, for these Seyfert galaxies studied in this work, sources with relatively larger UV to X-ray lags do have broader Fe Ka lines, indicative of relativistic broadening due to more compact corona and smaller innermost disk radius. If confirmed with more X-ray and UV monitoring campaigns, this interesting discovery would provide a new probe to the inner disk/corona.

preprint2019arXiv

Mining for Candidates of Galactic Stellar-mass Black Hole Binaries with LAMOST

We study the prospects of searching for black hole (BH) binary systems with a stellar-mass BH and a non-compact visible companion, by utilizing the spectroscopic data of Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST). We simulate the Galactic BH binary population and determine its optical visibility by considering the stellar synthetic population model and the distributions of binary orbital parameters. By convolving the visibility of BH binaries with the LAMOST detection sensitivity, we predict that $\gtrsim$ 400 candidate BH binaries can be found by the low-resolution, non-time-domain survey, and $\sim$ 50-350 candidates by the LAMOST ongoing medium-resolution, time-domain spectroscopic survey. Most of the candidates are short-period (0.2-2 days) binaries with M-, K-, G-, or F-type companions, in which $\sim$ 47% have a mass function (the lower limit of the BH mass) larger than 3 $M_{\odot}$. By complementing the LAMOST spectroscopic data with other photometric/spectroscopic surveys or follow-up observations, these candidates could be confirmed. Therefore, by exploring the LAMOST data, we can enlarge the sample of dynamically confirmed BH binaries significantly, which can improve our understanding of the mass distribution of BHs and the stellar evolution model.

preprint2019arXiv

Understanding Broad Mg II Variability in Quasars with Photoionization: Implications for Reverberation Mapping and Changing-Look Quasars

The broad Mg II line in quasars has distinct variability properties compared with broad Balmer lines: it is less variable, and usually does not display a &#34;breathing&#34; mode, the increase in the average cloud distance when luminosity increases. We demonstrate that these variability properties of Mg II can be reasonably well explained by simple Locally Optimally Emitting Cloud (LOC) photoionization models, confirming earlier photoionization results. In the fiducial LOC model, the Mg II-emitting gas is on average more distant from the ionizing source than the H$α$/H$β$ gas, and responds with a lower amplitude to continuum variations. If the broad-line region (BLR) is truncated at a physical radius of $\sim 0.3$ pc (for a $10^{8.5}M_{\odot}$ BH accreting at Eddington ratio of 0.1), most of the Mg II flux will always be emitted near this outer boundary and hence will not display breathing. These results indicate that reverberation mapping results on broad Mg II, while generally more difficult to obtain due to the lower line responsivity, can still be used to infer the Mg II BLR size and hence black hole mass. But it is possible that Mg II does not have a well defined intrinsic BLR size-luminosity relation for individual quasars, even though a global one for the general population may still exist. The dramatic changes in broad H$α$/H$β$ emission in the observationally-rare changing-look quasars are fully consistent with photoionization responses to extreme continuum variability, and the LOC model provides natural explanations for the persistence of broad Mg II in changing-look quasars defined on H$α$/H$β$, and the rare population of broad Mg II emitters in the spectra of massive inactive galaxies.