Researcher profile

Feng Yuan

Feng Yuan contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
54works
0followers
14topics
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

54 published item(s)

preprint2025arXiv

Radio signatures of AGN-wind-driven shocks in elliptical galaxies: From simulations to observations

We investigate the synchrotron emission signatures of shocks driven by active galactic nucleus (AGN) wind in elliptical galaxies based on our two-dimensional axisymmetric hydrodynamic MACER numerical simulations. Using these simulation data, we calculate the synchrotron radiation produced by nonthermal electrons accelerated at shocks, adopting reasonable assumptions for the magnetic field and relativistic electron distribution (derived from diffusive shock acceleration theory), and predict the resulting observational signatures. In our fiducial model, shocks driven by AGN winds produce synchrotron emission with luminosities of approximately $10^{29}\,\mathrm{erg\,s^{-1}\,Hz^{-1}}$ in the radio band (0.5-5 GHz), with spectral indices of $α\approx -0.4$ to $-0.6$ during the strongest shock phases, gradually steepening to about $-0.8$ to $-1.4$ as the electron population ages. Spatially, the emission is initially concentrated in regions of strong shocks, later expanding into more extended, diffuse structures. We also apply our model to the dwarf elliptical galaxy Messier 32 (M32), and find remarkable consistency between our simulated emission and the observed nuclear radio source, suggesting that this radio component likely originates from hot-wind-driven shocks. Our results indicate that AGN winds not only influence galaxy gas dynamics through mechanical energy input but also yield direct observational evidence via nonthermal radiation. With the advent of next-generation radio facilities such as the FAST Core Array, SKA, and ngVLA, these emission signatures serve as important probes for detecting and characterizing AGN feedback.

preprint2024arXiv

Dynamical excitations of one-dimensional Fulde-Ferrell pairing Fermi superfluid

We theoretically investigate a one-dimensional Fulde-Ferrell Fermi superfluid at a finite effective Zeeman field $h$, and study entire dynamical excitations related to density perturbation. By calculating the density dynamic structure factor, we find anisotropic dynamical excitations in both collective modes and single-particle excitations. Along the direction of centre-of-mass momentum $p$, there are two obvious gapless collective modes with different speed. The lower collective modes is from the usual gauge symmetry breaking and has a larger speed than the one in the negative direction of $p$. The higher one is due to the direction spontaneous symmetry breaking of centre-of-mass momentum $p$, and separates two kinds of single-particle excitations in the positive $p$ direction. However, this higher mode disappears in the opposite direction of $p$, where two single-particle excitations overlap with each other. These signals of dynamical excitations can do help to distinguish Fulde-Ferrell superfluid from the conventional Bardeen-Cooper-Schrieffer superfluid in the future experiment.

preprint2024arXiv

Proton's gluon GPDs at large skewness and gravitational form factors from near threshold heavy quarkonium photo-production

We study the exclusive near-threshold photo-production of heavy quarkonium in the framework of the generalized parton distribution (GPD) factorization, taking the $J/ψ$ production as an example. Due to the threshold kinematics, the Compton-like amplitudes are related to gluon GPDs at large skewness $ξ$, distinct from the common kinematics in asymptotic high energy where the skewness is typically small. We discuss the nature of large-$ξ$ expansion of these amplitudes in terms of the moments of gluon GPDs in the large-$ξ$ limit. Based on that, we propose several ways to extract the first few moments of the gluon GPDs from these amplitudes, with the leading ones corresponding to the gluonic gravitational or energy-momentum tensor form factors (GFFs). We apply these methods to analyze the recent near-threshold $J/ψ$ production measurements by the $J/ψ$ 007 experiment and GlueX collaboration, and find that the $ξ$-scaling of the measured differential cross sections is consistent with the asymptotic behavior. However, the current data are not accurate enough yet for a complete determination of the gluonic GFFs, and therefore we consider some prospects for better extractions in the future.

preprint2022arXiv

Azimuthal Angular Asymmetry of Soft Gluon Radiation in Jet Production

We investigate the impact of soft gluon resummation on the azimuthal angle correlation between the total and relative momenta of two energetic final state particles (jets). We show that the initial and final state radiations induce sizable $\cos(ϕ)$ and $\cos (2ϕ)$ asymmetries in single jet and dijet events, respectively. We numerically evaluate the magnitude of these asymmetries for a number of processes in collider experiments, including diffractive dijet and dilepton production in ultraperipheral $pA$ and $AA$ collisions, inclusive and diffractive dijet production at the EIC and inclusive dijet production in $pp$ collisions at the LHC. In particular, the $\cos (2ϕ)$ asymmetry of perturbative origin can dominate over the primordial asymmetry due to the linearly polarized gluon distribution.

preprint2022arXiv

Formation of episodic jets and associated flares from black hole accretion systems

Episodic ejections of blobs (episodic jets) are widely observed in black hole sources and usually associated with flares. In this paper, by performing and analyzing three dimensional general relativity magnetohydrodynamical numerical simulations of accretion flows, we investigate their physical mechanisms. We find that magnetic reconnection occurs in the accretion flow, likely due to the turbulent motion and differential rotation of the accretion flow, resulting in flares and formation of flux ropes. Flux ropes formed inside of 10-15 gravitational radii are found to mainly stay within the accretion flow, while flux ropes formed beyond this radius are ejected outward by magnetic forces and form the episodic jets. These results confirm the basic scenario proposed in Yuan et al.(2009). Moreover, our simulations find that the predicted velocity of the ejected blobs is in good consistency with observations of Sgr A*, M81, and M87. The whole processes are found to occur quasi-periodically, with the period being the orbital time at the radius where the flux rope is formed. The predicted period of flares and ejections is consistent with those found from the light curves or image of Sgr A*, M87, and PKS 1510-089. The possible applications to protostellar accretion systems are discussed.

preprint2022arXiv

Large-scale Dynamics of Winds Driven by Line Force from a Thin Accretion Disk

Winds play a significant role in active galactic nuclei feedback process. Previous simulations studying winds only focus on a small dynamical range. Therefore, it is unknown how far the winds can go and what the properties of the winds will be if they can move to large radii. We perform simulations to study the large scale dynamics of winds driven by line force. We find that the properties of the winds depend on both black hole mass ($M_{BH}$) and accretion disk luminosity. When the accretion disk luminosity is $0.6L_{edd}$ ($L_{edd}$ being Eddington luminosity), independent of $M_{BH}$, the winds have kinetic energy flux exceeding $1\% L_{edd}$ and can escape from the black hole potential. For the case with the accretion disk luminosity equaling 0.3$L_{edd}$, the strength of the winds decreases with the decrease of $M_{BH}$. If $M_{BH}$ decreases from $10^9$ to $10^6$ solar mass ($M_\odot$), the winds kinetic energy flux decreases from $\sim 0.01 L_{edd}$ to $ \sim 10^{-6} L_{edd}$. In case of $M_{BH}\geq 10^7 M_\odot$, winds can escape from black hole potential. In the case of $M_{BH}=10^6 M_\odot$, the winds can not escape. We find that for the ultra-fast winds observed in hard X-ray bands (\citealt{Gofford et al. 2015}), the observed dependence of the mass flux and the kinetic energy flux on accretion disk luminosity can be well produced by line force driven winds model. We also find that the properties of the ultra-fast winds observed in soft X-ray bands can be explained by the line force driven winds model.

preprint2022arXiv

Magnetic Spirals in Accretion Flows Originated from Misaligned Magnetic Field

Misalignment between rotation and magnetic field has been suggested to be one type of physical mechanisms which can easen the effects of magnetic braking during collapse of cloud cores leading to formation of protostellar disks. However, its essential factors are poorly understood. Therefore, we perform a more detailed analysis of the physics involved. We analyze existing simulation data to measure the system torques, mass accretion rates and Toomre Q parameters. We also examine the presence of shocks in the system. While advective torques are generally the strongest, we find that magnetic and gravitational torques can play substantial roles in how angular momentum is transferred during the disk formation process. Magnetic torques can shape the accretion flows, creating two-armed magnetized inflow spirals aligned with the magnetic field. We find evidence of an accretion shock that is aligned according to the spiral structure of the system. Inclusion of ambipolar diffusion as explored in this work has shown a slight influence in the small scale structures but not in the main morphology. We discuss potential candidate systems where some of these phenomena could be present.

preprint2022arXiv

Magnetically driven accretion disc winds: the role of gas thermodynamics and comparison to ultra-fast outflows

Winds are commonly observed in luminous active galactic nuclei (AGNs). A plausible model of those winds is magnetohydrodynamic (MHD) disc winds. In the case of disc winds from a thin accretion disc, isothermal or adiabatic assumption is usually adopted in such MHD models. In this work we perform two-dimensional MHD simulations implementing different thermal treatments (isothermal, adiabatic and radiative) to study their effects on winds from a thin accretion disc. We find that both the isothermal model and the adiabatic model overestimate the temperature, underestimate the power of disc winds, and cannot predict the local structure of the winds, compared to the results obtained by solving the energy equation with radiative cooling and heating. Based on the model with radiative cooling and heating, the ionization parameter, the column density and the velocity of the disc winds have been compared to the observed ultrafast outflows (UFOs). We find that in our simulations the UFOs can only be produced inside hundreds of Schwarzschild radius. At much larger radii, no UFOs are found. Thus, the pure MHD winds cannot interpret all the observed UFOs.

preprint2022arXiv

Observational Constraints on Direct Electron Heating in the Hot Accretion Flows in Sgr A* and M87*

An important parameter in the theory of hot accretion flows around black holes is $δ$, which describes the fraction of ``viscously'' dissipated energy in the accretion flow that goes directly into heating electrons. For a given mass accretion rate, the radiative efficiency of a hot accretion flow is determined by $δ$. Unfortunately, the value of $δ$ is hard to determine from first principles. The recent Event Horizon Telescope Collaboration (EHTC) results on M87* and Sgr A* provide us with a different way of constraining $δ$. By combining the mass accretion rates in M87* and Sgr A* estimated by the EHTC with the measured bolometric luminosities of the two sources, we derive good constraints on the radiative efficiencies of the respective accretion flows. In parallel, we use a theoretical model of hot magnetically arrested disks (MAD) to calculate the expected radiative efficiency as a function of $δ$ (and accretion rate). By comparing the EHTC-derived radiative efficiencies with the theoretical results from MAD models, we find that Sgr A* requires $δ\ga 0.3$. %with the most likely value being $δ\sim 0.5$. A similar comparison in the case of M87* gives inconclusive results as there is still a large uncertainty in the accretion rate in this source.

preprint2022arXiv

Pressure-induced superconductivity in flat-band Kagome compounds Pd$_3$P$_2$(S$_{1-x}$Se$_x$)$_8$

We performed high-pressure transport studies on the flat-band Kagome compounds, Pd$_3$P$_2$(S$_{1-x}$Se$_x$)$_8$ ($x$ = 0, 0.25), with a diamond anvil cell. For both compounds, the resistivity exhibits an insulating behavior with pressure up to 17 GPa. With pressure above 20 GPa, a metallic behavior is observed at high temperatures in Pd$_3$P$_2$S$_8$, and superconductivity emerges at low temperatures. The onset temperature of superconducting transition $T_{\rm C}$ rises monotonically from 2 K to 4.8 K and does not saturate with pressure up to 43 GPa. For the Se-doped compound Pd$_3$P$_2$(S$_{0.75}$Se$_{0.25}$)$_8$, the $T_{\rm C}$ is about 1.5 K higher than that of the undoped one over the whole pressure range, and reaches 6.4 K at 43 GPa. The upper critical field with field applied along the $c$ axis at typical pressures is about 50$\%$ of the Pauli limit, suggesting a 3D superconductivity. The Hall coefficient in the metallic phase is low and exhibits a peaked behavior at about 30 K, which suggests either a multi-band electronic structure or an electron correlation effect in the system.

preprint2022arXiv

The Accretion flow in M87 is really MAD

The supermassive black holes in most galaxies in the universe are powered by hot accretion flows. Both theoretical analysis and numerical simulations have indicated that, depending on the degree of magnetization, black hole hot accretion flow is divided into two modes, namely SANE (standard and normal evolution) and MAD (magnetically arrested disk). It has been an important question which mode the hot accretion flows in individual sources should belong to in reality, SANE or MAD. This issue has been investigated in some previous works but they all suffer from various uncertainties. By using the measured rotation measure values in the prototype low-luminosity active galactic nuclei in {M87} at 2, 5, and 8 GHz along the jet at various distances from the black hole, combined with three dimensional general relativity magnetohydrodynamical numerical simulations of SANE and MAD, we show in this paper that the predicted rotation measure values by MAD are well consistent with observations, while the SANE model overestimates the rotation measure by over two orders of magnitude thus is ruled out.

preprint2022arXiv

The Photon Ring in M87*

We report measurements of the gravitationally lensed secondary image -- the first in an infinite series of so-called "photon rings" -- around the supermassive black hole M87* via simultaneous modeling and imaging of the 2017 Event Horizon Telescope (EHT) observations. The inferred ring size remains constant across the seven days of the 2017 EHT observing campaign and is consistent with theoretical expectations, providing clear evidence that such measurements probe spacetime and a striking confirmation of the models underlying the first set of EHT results. The residual diffuse emission evolves on timescales comparable to one week. We are able to detect with high significance a southwestern extension consistent with that expected from the base of a jet that is rapidly rotating in the clockwise direction. This result adds further support to the identification of the jet in M87* with a black hole spin-driven outflow, launched via the Blandford-Znajek process. We present three revised estimates for the mass of M87* based on identifying the modeled thin ring component with the bright ringlike features seen in simulated images, one of which is only weakly sensitive to the astrophysics of the emission region. All three estimates agree with each other and previously reported values. Our strongest mass constraint combines information from both the ring and the diffuse emission region, which together imply a mass-to-distance ratio of $4.20^{+0.12}_{-0.06}~μ{\rm as}$ and a corresponding black hole mass of $(7.13\pm0.39)\times10^9M_\odot$, where the error on the latter is now dominated by the systematic uncertainty arising from the uncertain distance to M87*.

preprint2021arXiv

EIC Physics from An All-Silicon Tracking Detector

The proposed electron-ion collider has a rich physics program to study the internal structure of protons and heavy nuclei. This program will impose strict requirements on detector design. This paper explores how these requirements can be satisfied using an all-silicon tracking detector, by consideration of three representative probes: heavy flavor hadrons, jets, and exclusive vector mesons.

preprint2021arXiv

Evidence for A Hot Wind from High-resolution X-ray Spectroscopic Observation of the Low-luminosity Active Galactic Nucleus in NGC 7213

Super-massive black holes (SMBHs) spend most of their lifetime accreting at a rate well below the Eddington limit, manifesting themselves as low-luminosity active galactic nuclei (LLAGNs). The prevalence of a hot wind from LLAGNs is a generic prediction by theories and numerical simulations of black hole accretion and is recently becoming a crucial ingredient of AGN kinetic feedback in cosmological simulations of galaxy evolution. However, direct observational evidence for this hot wind is still scarce. In this work, we identify significant Fe XXVI Ly$α$ and Fe XXV K$α$ emission lines from high-resolution Chandra grating spectra of the LLAGN in NGC\,7213, a nearby Sa galaxy hosting a $\sim10^8\rm~M_\odot$ SMBH, confirming previous work. We find that these lines exhibit a blueshifted line-of-sight velocity of $\sim1100\rm~km s^{-1}$ and a high XXVI Ly$α$ to XXV K$α$ flux ratio implying for a $\sim16$ keV hot plasma. By confronting these spectral features with synthetic X-ray spectra based on our custom magnetohydrodynamical simulations, we find that the high-velocity, hot plasma is naturally explained by the putative hot wind driven by the hot accretion flow powering this LLAGN. Alternative plausible origins of this hot plasma, including stellar activities, AGN photoionization and the hot accretion flow itself, are quantitatively disfavored. The inferred kinetic energy and momentum carried by the wind can serve as strong feedback to the environment. We compare NGC\,7213 to M81*, in which strong evidence for a hot wind was recently presented, and discuss implications on the universality and detectability of hot winds from LLAGNs.

preprint2021arXiv

Near Threshold Heavy Quarkonium Photoproduction at Large Momentum Transfer

Perturbative QCD is applied to investigate the near threshold heavy quarkonium photoproduction at large momentum transfer. We take into account the contributions from the leading three-quark Fock states of the nucleon. The dominant contribution comes from the three-quark Fock state with one unit quark orbital angular momentum (OAM) whereas that from zero quark OAM is suppressed at the threshold. From our analysis, we also show that there is no direct connection between the near threshold heavy quarkonium photoproduction and the gluonic gravitational form factors of the nucleon. Based on the comparison between our result and recent GlueX data of $J/ψ$ photoproduction, we make predictions for $ψ'$ and $Υ$ (1S,2S) states which can be tested in future experiments.

preprint2021arXiv

Probing gluon helicity with heavy flavor at the EIC

We propose a new measurement of the heavy flavor hadron double spin asymmetry in deep-inelastic scattering at a future Electron-Ion Collider (EIC) to constrain the polarized gluon distribution function inside the proton. Statistical projection on $D^0$ meson double spin asymmetry is calculated with an EIC central detector using an all-silicon tracker and vertexing subsystem. A first impact study was done by interpreting pseudo-data at next-to-leading order in QCD. The sensitivity of the experimental observable in constraining gluon helicity distribution in a wide range of parton momentum fraction $x$ has been investigated considering different beam energy configurations. This measurement complements the inclusive spin-dependent structure function measurement and provides an opportunity to constrain the gluon helicity distribution in the moderate $x$ region.

preprint2021arXiv

Understanding Patterns of Users Who Repost Censored Posts on Weibo

In this study, we focus on understanding patterns of users whose repost contents would later be censored on Weibo, a counterpart of Twitter in China as a social media platform. Little is known about the way regulations and censorship work in this indigenous platform and what role it plays in affecting users' expression of ideas. We collect over a million reposts from over 18,000 users and investigate the patterns of users whose reposts contain content that is no longer visible to the public, from the perspective of user location, device, gender, social capital as well as certified status. We find that user characteristics play an important role in affecting behaviors on Weibo.

preprint2020arXiv

Discovery of oscillations above 200 keV in a black hole X-ray binary with Insight-HXMT

Low-frequency quasi-periodic oscillations (LFQPOs) are commonly found in black hole X-ray binaries, and their origin is still under debate. The properties of LFQPOs at high energies (above 30 keV) are closely related to the nature of the accretion flow in the innermost regions, and thus play a crucial role in critically testing various theoretical models. The Hard X-ray Modulation Telescope (Insight-HXMT) is capable of detecting emissions above 30 keV, and is therefore an ideal instrument to do so. Here we report the discovery of LFQPOs above 200 keV in the new black hole MAXI J1820+070 in the X-ray hard state, which allows us to understand the behaviours of LFQPOs at hundreds of kiloelectronvolts. The phase lag of the LFQPO is constant around zero below 30 keV, and becomes a soft lag (that is, the high-energy photons arrive first) above 30 keV. The soft lag gradually increases with energy and reaches ~0.9s in the 150-200 keV band. The detection at energies above 200 keV, the large soft lag and the energy-related behaviors of the LFQPO pose a great challenge for most currently existing models, but suggest that the LFQPO probably originates from the precession of a small-scale jet.

preprint2020arXiv

Dynamical structure factors of a two-dimensional Fermi superfluid within random phase approximation

Based on random phase approximation (RPA), we numerically calculate dynamical structure factors of a balanced two-dimensional (2D) Fermi superfluid, and discuss their energy, momentum and interaction strength dependence in the 2D BEC-BCS crossover. At a small transferred momentum, a stable Higgs mode is observed in the unitary 2D Fermi superfluid gas where the particle-hole symmetry is not satisfied. Stronger interaction strength will make the visibility of the dispersion of Higgs mode harder to be observed. We also discuss the dimension effect and find that the signal of the Higgs mode in two dimension is more obvious than that in 3D case. At a large transferred momentum regime, stronger interaction strength will induce the weight of the molecules excitation increasing, while in verse the atomic one decreasing, which shows the pairing information of Fermi superfluid. The theoretical results qualitatively agree with the corresponding Quantum Monte Carlo data.

preprint2020arXiv

Factorization and its Breaking in Dijet Single Transverse Spin Asymmetries in $pp$ Collisions

We study factorization in single transverse spin asymmetries for dijet production in proton-proton collisions, by considering soft gluon radiation at one-loop order. We show that the associated transverse momentum dependent (TMD) factorization is valid at the leading logarithmic level. At next-to-leading-logarithmic (NLL) accuracy, however, we find that soft gluon radiation generates terms in the single transverse spin dependent cross section that differ from those known for the unpolarized case. As a consequence, these terms cannot be organized in terms of a spin independent soft factor in the factorization formula. We present leading logarithmic predictions for the single transverse spin dijet asymmetry for $pp$ collisions at RHIC, based on quark Sivers functions constrained by semi-inclusive deep inelastic scattering data. We hope that our results will contribute to a better understanding of TMD factorization breaking effects at NLL accuracy and beyond.

preprint2020arXiv

From haloes to galaxies -- II. The fundamental relations in star formation and quenching

Star formation and quenching are two of the most important processes in galaxy formation and evolution. We explore in the local Universe the interrelationships among key integrated galaxy properties, including stellar mass $M_*$, star formation rate (SFR), specific SFR (sSFR), molecular gas mass $M_{\rm H_2}$, star formation efficiency (SFE) of the molecular gas and molecular gas to stellar mass ratio $μ$. We aim to identify the most fundamental scaling relations among these key galaxy properties and their interrelationships. We show the integrated $M_{\rm H_2}$-SFR, SFR-$M_*$ and $M_{\rm H_2}$-$M_*$ relation can be simply transformed from the $μ$-sSFR, SFE-$μ$ and SFE-sSFR relation, respectively. The transformation, in principle, can increase or decrease the scatter of each relation. Interestingly, we find the latter three relations all have significantly smaller scatter than the former three corresponding relations. We show the probability to achieve the observed small scatter by accident is extremely close to zero. This suggests that the smaller scatters of the latter three relations are driven by a more fundamental physical connection among these quantities. We then show the large scatters in the former relations are due to their systematic dependence on other galaxy properties, and on star formation and quenching process. We propose the sSFR-$μ$-SFE relation as the Fundamental Formation Relation (FFR), which governs the star formation and quenching process, and provides a simple framework to study galaxy evolution. Other scaling relations, including integrated Kennicutt-Schmidt law, star-forming main sequence and molecular gas main sequence, can all be derived from the FFR.

preprint2020arXiv

Large-scale dynamics of winds originated from black hole accretion flows: (I) Hydrodynamics

Winds from black hole accretion flows are ubiquitous. Previous works mainly focus on the launching of wind in the accretion flow scale. It still remains unclear how far the winds can propagate outward and what is their large-scale dynamics. As the first paper of this series, we study the large-scale dynamics of thermal wind beyond accretion scales via analytical and numerical methods. Boundary conditions, which are crucial to our problem, are analyzed and presented based on the small-scale simulations combined with observations of winds. Both black hole and galaxy potential are taken into account. For winds originated from hot accretion flows, we find that the wind can reach to large scales. The radial profiles of velocity, density, and temperature can be approximated by $v_r\approx v_{r0}, ρ\approx ρ_{0}(r/r_0)^{-2}$, and $T\approx T_0 (r/r_0)^{-2(γ-1)}$, where $v_{r0}, ρ_0, T_0$ are the velocity, density, and temperature of winds at the boundary $r_0(\equiv 10^3 r_g)$, $γ$ is the polytropic index. During the outward propagation, the enthalpy and the rotational energy compensate the increase of gravitational potential. For thin disks, we find that because the Bernoulli parameter is smaller, winds cannot propagate as far as the hot winds, but stop at a certain radius where the Bernoulli parameter is equal to the potential energy. Before the winds stop, the profiles of dynamical quantities can also be approximated by the above relations. In this case the rotational energy alone compensates the increase of the potential energy.

preprint2020arXiv

Large-scale dynamics of winds originated from black hole accretion flows: (II) Magnetohydrodynamics

Winds from black hole accretion disks are essential ingredients in understanding the coevolution between the supermassive black hole and its host galaxy. The great difference of dynamical ranges from small-scale accretion disk simulations to large-scale or cosmological simulations places barriers to track wind kinematics. In the first paper of this series, we have studied the dynamics of disk winds from the outer edge of the accretion disk toward galaxy scales in the hydrodynamical framework. In this paper, we further incorporate magnetic fields to understand the wind dynamics by adopting one-dimensional magnetohydrodynamical (MHD) model, with boundary conditions set for hot accretion flows. The geometry of poloidal magnetic field is prescribed as a straight line with an angle $θ=45^\circ$ from the rotational axis, and the strength satisfies the divergence free condition. The wind solution is achieved through requesting gas to pass through the slow, Alfvén and fast magneto-sonic points smoothly. Physical quantities are found to show a power-law dependence on cylindrical radius $R$ beyond the fast magneto-sonic point, for which $ρ\propto R^{-2}, v_{\rm p}\propto {\rm const.}, v_{\rm ϕ}\propto R^{-1}, B_{\rm ϕ}\propto R^{-1},$ and $ β\propto ρ^{γ-1}$. The magnetization of wind is dominant in determining the wind properties. The wind is accelerated to a greater terminal velocity with strong magnetization ($v_{\rm Ap0}>1$) compared to the hydrodynamical case, which the magnetic pressure gradient dominates and the centrifugal potential converts to the kinetic energy. The dependance of wind physical quantities on magnetization, temperature, field line angular velocity, and adiabatic index is also discussed.

preprint2020arXiv

MAMO: Memory-Augmented Meta-Optimization for Cold-start Recommendation

A common challenge for most current recommender systems is the cold-start problem. Due to the lack of user-item interactions, the fine-tuned recommender systems are unable to handle situations with new users or new items. Recently, some works introduce the meta-optimization idea into the recommendation scenarios, i.e. predicting the user preference by only a few of past interacted items. The core idea is learning a global sharing initialization parameter for all users and then learning the local parameters for each user separately. However, most meta-learning based recommendation approaches adopt model-agnostic meta-learning for parameter initialization, where the global sharing parameter may lead the model into local optima for some users. In this paper, we design two memory matrices that can store task-specific memories and feature-specific memories. Specifically, the feature-specific memories are used to guide the model with personalized parameter initialization, while the task-specific memories are used to guide the model fast predicting the user preference. And we adopt a meta-optimization approach for optimizing the proposed method. We test the model on two widely used recommendation datasets and consider four cold-start situations. The experimental results show the effectiveness of the proposed methods.

preprint2020arXiv

Proton spin after 30 years: what we know and what we don't?

More than three decades has passed since the European Muon Collaboration published the first surprising result on the spin structure of the proton. Much theoretical and experimental progress has been made in understanding the origins of the proton spin. In this review, we will discuss what we have learned so far, what are still missing, and what we shall expect to learn from the upcoming experiments including JLab 12 GeV and Electron-Ion Collider. In particular, we focus on first principles calculations and experimental measurements of the total gluon helicity $ΔG$, and quark and gluon orbital angular momenta.

preprint2020arXiv

Resolving the Nuclear Radio Emission from M32 with Very Large Array

The Local Group dwarf elliptical galaxy M32 hosts one of the nearest and most under-luminous super-massive black holes (SMBHs) ever known, offering a rare opportunity to study the physics of accreting SMBHs at the most quiescent state. Recent Very Large Array (VLA) observations have detected a radio source at the nucleus of M32, which is suggested to be the radio counterpart of the SMBH. To further investigate the radio properties of this nuclear source, we have conducted follow-up, high-resolution VLA observations in four epochs between 2015--2017, each with dual frequencies. At 6 GHz, the nuclear source is resolved under an angular resolution of $\sim$0\farcs4, exhibiting a coreless, slightly lopsided morphology with a detectable extent of $\sim$2.5 \arcsec ($\sim$10 parsec). No significant variability can be found among the four epochs. At 15 GHz, no significant emission can be detected within the same region, pointing to a steep intrinsic radio spectrum (with a 3\,$σ$ upper limit of -1.46 for the spectral index). We discuss possible scenarios for the nature of this nuclear source and conclude that a stellar origin, in particular planetary nebulae, X-ray binaries, supernova remnants or diffuse ionized gas powered by massive stars, can be ruled out.Instead, the observed radio properties can be explained by synchrotron radiation from a hypothetical wind driven by the weakly accreting SMBH.

preprint2020arXiv

Transverse spin sum rule of the proton

The transversely-polarized state of a proton with arbitrary momentum is not an eigenstate of transverse angular momentum operator. The latter does not commute with the QCD Hamiltonian. However, the expectation value of the transverse angular momentum in the state is well-defined and grows proportionally to the energy of the particle. The transverse spin content of the proton is analyzed in terms of the QCD angular momentum structure. In particular, we reconfirm that the generalized parton distributions $H+E$ provide transverse orbital angular momentum densities of quarks and gluons in the infinite momentum frame.

preprint2019arXiv

Can We Observe Jet $P_T$-broadening in Heavy-Ion Collisions at the LHC?

We investigate the effect of PT-broadening on jet substructure observables in heavy-ion collisions at the LHC. As an example, we focus on the opening angle of the two branches that satisfy the soft drop grooming condition in a highly energetic jet. The medium modification of the angular distribution can provide important information on the jet transport properties of hot QCD matter. In addition, we take into account a change of the overall fraction of quark and gluon jets in heavy-ion collisions. We comment on the comparison to a recent measurement from the ALICE Collaboration.

preprint2019arXiv

Gluonic Probe for the Short Range Correlation in Nucleus

We investigate the gluonic probe to the nucleon-nucleon short range correlation (SRC) in nucleus through heavy flavor production in deep inelastic scattering (DIS). The relevant EMC effects of $F_2^{c\bar c}$ structure function will provide a universality test of the SRCs which have been extensively studied in the quark-channel. These SRCs can also be studied through the sub-threshold production of heavy flavor in $eA$ collisions at the intermediate energy electron-ion collider, including open Charm and $J/ψ$ ($Υ$) production.

preprint2019arXiv

Medium-Induced Transverse Momentum Broadening via Forward Dijet Correlations

Dijet azimuthal angle correlation is arguably one of the most direct probes of the medium-induced broadening effects. The evidence for such broadening, however, is not yet clearly observed within the precision of current mid-rapidity measurements at RHIC and the LHC. We show that the dijet correlation in forward rapidity from the future LHC RUN3, aided by forward detector upgrades, can reveal this broadening thanks to the steeper jet spectra, suppressed vacuum radiations and lower underlying event background, with a sensitivity comparable to that of the future high-luminosity Au+Au run at RHIC. Dijet correlation measurements at the two facilities together can provide powerful constraints on the temperature dependence of medium transport properties.

preprint2019arXiv

On the Role of Hot Feedback Mode in Active Galactic Nuclei Feedback in an Elliptical Galaxy

Depending on the value of the accretion rate, black hole accretion is divided into cold and hot modes. The two modes have distinctly different physics and correspond to two feedback modes. Most previous feedback works either focus only on one mode, or the accretion physics is not always properly adopted. Here, by performing hydrodynamical numerical simulations of AGN feedback in an elliptical galaxy, we show that including both is important, and gives different results from including just one or the other. We specifically focus on the wind and radiation (but neglecting the jet) feedback in the hot mode, to explore whether this particular mode of feedback can play any important feedback role. For this aim, we have run two test models. In one model, we always adopt the cold mode no matter what value the accretion rate is; in another model, we turn off the AGN once it enters into the hot mode. We have calculated the AGN light curves, black hole growth, star formation, and AGN energy duty-cycle; and compared the results with the model in which the two modes are correctly included. Important differences are found. For example, if we were to adopt only the cold mode, the total mass of newly formed stars would become two orders of magnitude smaller, and the fraction of energy ejected within the high accretion regime (i.e., $L_{\rm BH} > 2\% L_{\rm Edd}$) would be too small to be consistent with observations. We have also investigated the respective roles of wind and radiation in the hot mode.

preprint2019arXiv

QCD Resummation in Hard Diffractive Dijet Production at the Electron-Ion Collider

Diffractive dijet production at the electron-ion collider (EIC) has been proposed to study the gluon Wigner distribution at small-$x$. We investigate the soft gluon radiation associated with the final state jets and an all order resummation formula is derived. We show that the soft gluon resummation plays an important role to describe E791 data on $π$-induced diffractive dijet production at Fermilab. Predictions for the EIC are presented, and we emphasize that the soft gluon resummation is an important aspect to explore the nucleon/nucleus tomography through these processes.

preprint2019arXiv

Sub-threshold $J/ψ$ and $Υ$ Production in $γA$ Collisions

We study sub-threshold heavy quarkonium ($J/ψ$ and $Υ$) photo-productions in $γA$ collisions as an independent test of the universality of the nucleon-nucleon short range correlation (SRC) in nuclear scattering processes. Just below the $γp$ threshold, the cross section is dominated by the mean field contribution of nucleons inside the nucleus. The SRC contributions start to dominate at lower photon energies, depending on the fraction of the SRC pairs in the target nucleus. We give an estimate of the cross sections in the sub-threshold region both for $J/ψ$ and $Υ$. This may be helpful for future measurements at JLab as well as at the Electron-Ion Collider in the U.S., and especially in China.

preprint2011arXiv

An Observation Concerning the Process Dependence of the Sivers Functions

The $k_\perp$-moment of a quark's Sivers function is known to be related to the corresponding twist-three quark-gluon correlation function $T_{q,F}(x, x)$. The two functions have been extracted from data for single-spin asymmetries in semi-inclusive deep inelastic scattering and in single-inclusive hadron production in $pp$ collisions, respectively. Performing a consistent comparison of the extracted functions, we find that they show a "sign mismatch": while the magnitude of the functions is roughly consistent, the $k_\perp$-moment of the Sivers function has opposite sign from that of $T_{q,F}(x, x)$, both for up and for down quarks. Barring any inconsistencies in our theoretical understanding of the Sivers functions and their process dependence, the implication of this mismatch is that either, the Sivers effect is not dominantly responsible for the observed single-spin asymmetries in $pp$ collisions or, the current semi-inclusive lepton scattering data do not sufficiently constrain the $k_\perp$-moment of the quark Sivers functions. Both possibilities strengthen the case for further experimental investigations of single-spin asymmetries in high-energy $pp$ and $ep$ scattering.

preprint2011arXiv

kt-factorization for Hard Processes in Nuclei

Two widely proposed kt-dependent gluon distributions in the small-x saturation regime are investigated using two particle back-to-back correlations in high energy scattering processes. The Weizsacker-Williams gluon distribution, interpreted as the number density of gluon inside the nucleus, is studied in the quark-antiquark jet correlation in deep inelastic scattering. On the other hand, the unintegrated gluon distribution, defined as the Fourier transform of the color-dipole cross section, is probed in the direct photon-jet correlation in pA collisions. Dijet-correlation in pA collisions depends on both gluon distributions through combination and convolution in the large Nc limit.

preprint2011arXiv

QCD Resummation for Single Spin Asymmetries

We study the transverse momentum dependent factorization for single spin asymmetries in Drell-Yan and semi-inclusive deep inelastic scattering processes at one-loop order. The next-to-leading order hard factors are calculated in the Ji-Ma-Yuan factorization scheme. We further derive the QCD resummation formalisms for these observables following the Collins-Soper-Sterman method. The results are expressed in terms of the collinear correlation functions from initial and/or final state hadrons coupled with the Sudakov form factor containing all order soft-gluon resummation effects. The scheme independent coefficients are calculated up to one-loop order.

preprint2011arXiv

Single Spin Asymmetry Scaling in the Forward Rapidity Region at RHIC

We investigate the scaling properties in inclusive hadron production and the associated single transverse spin asymmetry in the forward rapidity region at RHIC. We find that the spin-averaged experimental data in both $pp$ and $dAu$ collisions demonstrates a transverse-momentum-dependent geometric scaling. We introduce the transverse momentum dependent Collins fragmentation function to study the scaling of the single transverse spin asymmetries. The general feature of the scaling analysis is consistent with the experimental observations, in particular, for the transverse momentum dependence of the spin asymmetries at RHIC energy. We further propose to probe the saturation scale of nucleus by measuring the spin asymmetry normalized by that in $pp$ scattering at low transverse momentum.

preprint2011arXiv

Universality of Unintegrated Gluon Distributions at small x

We systematically study dijet production in various processes in the small-$x$ limit and establish an effective $k_t$-factorization for hard processes in a system with dilute probes scattering on a dense target. We find that the well-known Weizsäcker-Williams gluon distribution can be directly probed in the quark-antiquark jet correlation in deep inelastic scattering and the dipole gluon distribution can be directly measured in the direct photon-jet correlation in $pA$ collisions. In the large-$N_c$ limit, the unintegrated gluon distributions involved in other different dijet channels in pA collisions are shown to be related to two widely proposed ones: the Weizsäcker-Williams gluon distribution and the dipole gluon distribution.

preprint2010arXiv

Accretion and ejection in Sgr A*

We review our current understanding to the accretion and ejection processes in Sgr A*. Roughly speaking, they correspond to the quiescent and flare states of the source respectively. The high-resolution {\it Chandra} observations to the gas at the Bondi radius combined with the Bondi accretion theory, the spectral energy distribution from radio to X-ray, and the radio polarization provide us strict constraints and abundant information to the theory of accretion. We review these observational results and describe how the advection-dominated accretion flow model explains these observations. Recently more attentions have been paid to flares in Sgr A*. Many simultaneous multi-wavelength campaigns have been conducted, aiming at uncovering the nature of flares. The main observational properties of flares are briefly reviewed. Especially, the time lag between the peaks of flare at two radio frequencies strongly indicates that the flare is associated with ejection of radio-emitting blobs from the underlying accretion flow. Such kind of episodic jets is distinctive from the continuous jets and are quite common in black hole systems. We introduce the magnetohydrodynamical model for the formation of episodic jets recently proposed based on the analogy with the theory of coronal mass ejection in the Sun. We point out that the various observational appearances of flares should be explained in the framework of this model, since ejection and flare originate from the same physical process.

preprint2010arXiv

Initial and Final State Interaction Effects in Small-x Quark Distributions

We study the initial and final state interaction effects in the transverse momentum dependent parton distributions in the small-$x$ saturation region. In particular, we discuss the quark distributions in the semi-inclusive deep inelastic scattering, Drell-Yan lepton pair production and dijet-correlation processes in $pA$ collisions. We calculate the quark distributions in the scalar-QED model and then extend to the color glass condensate formalism in QCD. The quark distributions are found universal between the DIS and Drell-Yan processes. On the other hand, the quark distribution from the $qq'\to qq'$ channel contribution to the dijet-correlation process is not universal. However, we find that it can be related to the quark distribution in DIS process by a convolution with the normalized unintegrated gluon distribution in the color glass condensate formalism in the large $N_c$ limit.

preprint2010arXiv

On the convective instability of hot radiative accretion flows

How many fraction of gas available at the outer boundary can finally fall onto the black hole is an important question. It determines the observational appearance of accretion flows, and is also related with the evolution of black hole mass and spin. Previous two-dimensional hydrodynamical simulations of hot accretion flows find that the flow is convectively unstable because of its inward increase of entropy. As a result, the mass accretion rate decreases inward, i.e., only a small fraction of accretion gas can fall onto the black hole, while the rest circulates in the convective eddies or lost in convective outflows. Radiation is usually neglected in these simulations. In many cases, however, radiative cooling is important. In the regime of the luminous hot accretion flow (LHAF), radiative cooling is even stronger than the viscous dissipation. In the one dimensional case, this implies that the inward increase of entropy will become slower or the entropy even decreases inward in the case of an LHAF. We therefore expect that convective instability becomes weaker or completely disappears when radiative cooling is important. To examine the validity of this expectation, in this paper we perform two-dimensional hydrodynamical simulations of hot accretion flows with strong radiative cooling. We find that compared to the case of negligible radiation, convection only becomes slightly weaker. Even an LHAF is still strongly convectively unstable, its radial profile of accretion rate correspondingly changes little. We find the reason is that the entropy still increases inward in the two-dimensional case.

preprint2010arXiv

Test of the Universality of Naive-time-reversal-odd Fragmentation Functions

We investigate the ''spontaneous'' hyperon transverse polarization in $e^+e^-$ annihilation and semi-inclusive deep inelastic scattering processes as a test of the universality of the naive-time-reversal-odd transverse momentum dependent fragmentation functions. We find that universality implies definite sign relations among various observables. This provides a unique opportunity to study initial/final state interaction effects in the fragmentation process and test the associated factorization.

preprint2006arXiv

Modeling the hard states of XTE J1550--564 during its 2000 outburst

We study hard states of the black-hole binary XTE J1550--564 during its 2000 outburst. In order to explain those states at their highest luminosities, $L\sim 10%$ of the Eddington luminosity, $L_{\rm E}$, we propose a specific hot accretion flow model. We point out that the highest values of the hard-state $L$ are substantially above the $L$ an advection-dominated accretion flow (ADAF) can produce, $\sim 0.4α^2 L_{\rm E}$, which is only $\sim (3$--$4)%L_{\rm E}$ even for $α$ as high as 0.3. On the other hand, we successfully explain the hard states with $L\sim (4$--$10)%$ using the luminous hot accretion flow (LHAF) model. As $10%L_{\rm E}$ is also roughly the highest luminosity an LHAF can produce, such an agreement between the predicted and observed highest luminosities provides by itself strong support for this model. Then, we study multi-waveband spectral variability during the 2000 outburst. In addition to the primary maxima in the optical light curves, secondary maxima were detected after the transition from the very high state to the hard state. We show that the secondary maxima are well modeled by synchrotron emission from a jet formed during the state transition. We argue that the absence of the corresponding secondary peak in the X-ray light curve indicates that the X-ray jet emission, regardless of its radiative process, synchrotron or its Comptonization, is not important in the hard state compared to the emission from the accretion flow.

preprint2005arXiv

Single-Transverse Spin Asymmetries: From DIS to Hadronic Collisions

We study single-spin asymmetries in semi-inclusive deep inelastic scattering with transversely polarized target. Based on the QCD factorization approach, we consider Sivers and Collins contributions to the asymmetries. We fit simple parameterizations for the Sivers and Collins functions to the recent HERMES data, and compare to results from COMPASS. Using the fitted parameterizations for the Sivers functions, we predict the single transverse spin asymmetries for various processes in $pp$ collisions at RHIC, including the Drell-Yan process and angular correlations in dijet and jet-plus-photon production. These asymmetries are found to be sizable at forward rapidities

preprint2005arXiv

Transverse Momentum Distribution Through Soft-Gluon Resummation in Effective Field Theory

We study resummation of transverse-momentum-related large logarithms generated from soft-gluon radiations in soft-collinear effective field theory. The anomalous dimensions of the effective quark and gluon currents, an important ingredient for the resummation, are calculated to two-loop order. the result at next-to-leading-log reproduces that obtained using the standard method for deep-inelastic scattering, Drell-Yan process, and Higgs production through gluon-gluon fusion. We comment on the extension of the calculation to next-to-next-to-leading logarithms.

preprint2004arXiv

Factorization of Large-x Quark Distributions in a Hadron

We present a factorization formula for valence quark distributions in a hadron in x-->1 limit. For the example of pion, we arrive at the form of factorization by analyzing momentum flow in the leading and high-order Feynman diagrams. The result confirms the well-known 1-x scaling rule to all orders in perturbation theory, providing the non-perturbative matrix elements for the infrared-divergence factors. We comment on re-summation of perturbative single and double logarithms in 1-x.

preprint2004arXiv

QCD Factorization for Semi-Inclusive Deep-Inelastic Scattering at Low Transverse Momentum

We demonstrate a factorization formula for semi-inclusive deep-inelastic scattering with hadrons in the current fragmentation region detected at low transverse momentum. To facilitate the factorization, we introduce the transverse-momentum dependent parton distributions and fragmentation functions with gauge links slightly off the light-cone, and with soft-gluon radiations subtracted. We verify the factorization to one-loop order in perturbative quantum chromodynamics and argue that it is valid to all orders in perturbation theory.

preprint2003arXiv

Quark Imaging in the Proton Via Quantum Phase-Space Distributions

We develop the concept of quantum phase-space (Wigner) distributions for quarks and gluons in the proton. To appreciate their physical content, we analyze the contraints from special relativity on the interpretation of elastic form factors, and examine the physics of the Feynman parton distributions in the proton's rest frame. We relate the quark Wigner functions to the transverse-momentum dependent parton distributions and generalized parton distributions, emphasizing the physical role of the skewness parameter. We show that the Wigner functions allow to visualize quantum quarks and gluons using the language of the classical phase space. We present two examples of the quark Wigner distributions and point out some model-independent features.

preprint2002arXiv

Parton Distributions in Light-Cone Gauge: Where Are the Final-State Interactions?

We show that the final-state interaction effects in the single target spin asymmetry discovered by Brodsky et al. can be reproduced by either a standard light-cone gauge definition of the parton distributions with a prescription of the light-cone singularities consistent with the light-cone gauge link, or a modified light-cone gauge definition with a gauge link involving the gauge potential at the spatial infinity.

preprint2001arXiv

Polarizations of J/ψand ψ' in hadroproduction at Tevatron in the k_t factorization approach

We present a calculation for the polarizations of $J/ψ$ and $ψ'$ produced in the hadron collisions at the Fermilab Tevatron. Various color octet channels including ${}^1S_0^{(8)}$, ${}^3P_J^{(8)}$, and ${}^3S_1^{(8)}$ as well as contributions from $χ_{cJ}$ decays are considered in the $k_t$ factorization approach. We find that in a rather wide range of the transverse momenta of $J/ψ$ and $ψ'$, the production rates could be dominated by the ${}^1S_0^{(8)}$ channel, and the predicted polarizations from the ${}^1S_0^{(8)}$ channel and $χ_{cJ}$ feeddown contributions are roughly compatible with the preliminary CDF data. This might provide a possible release from the conflict between the NRQCD collinear parton model calculations and the CDF data.

preprint2000arXiv

Inelastic electroproduction of $η_c$ at ep colliders

Using the nonrelativistic QCD factorization formalism, we calculate the electroproduction cross sections of $η_c$ in ep collisions, including the contribution from both the transverse photon and the longitudinal photon. For this process the color-singlet contribution vanishes up to the next to leading order perturbative QCD calculations. The dominant contribution comes from the color-octet ${}^1S_0^{(8)}$ subprocess. The nonperturbative color-octet matrix element of ${}^1S_0^{(8)}$ of $η_c$ is related to that of ${}^3S_1^{(8)}$ of $J/ψ$ by heavy quark spin symmetry, and the latter can be determined from the direct production of $J/ψ$ at large transverse momentum at the Fermilab Tevatron. The measurement of this process at DESY HERA can be viewed as another independent test for the color-octet production mechanism

preprint2000arXiv

χ_{cJ} polarizations at the Fermilab Tevatron

We propose the measurement of $χ_{cJ}$ polarizations at high energy hadron colliders to study heavy quarkonium production mechanism. We find that the color-singlet model in the $k_t$ factorization approach predicts very different behavior for $χ_{cJ}$ polarizations at the Fermilab Tevatron compared with the NRQCD predictions in the collinear parton model. In the color-singlet $k_t$ factorization approach, for both $χ_{c1}$ and $χ_{c2}$ productions, the helicity $h=0$ states dominate over other helicity states at large $p_T$. These properties are very useful in distinguishing between the two production mechanisms which are related to the interesting issue of $Jψ$ and $ψ'$ polarizations, and may provide a crucial test for the $k_t$ factorization approach.

preprint1999arXiv

Diffractive charm jet production at hadron colliders in the two-gluon exchange model

We present a calculation of diffractive charm jet production at hadron colliders in perturbative QCD based on the two-gluon exchange model. Differing from the previous calculations, we abandon the use of the effective color-singlet vitual gluon simplification, and use the real gluon in the calculations for the partonic process. We find the final result is free of the linear singularities due to the small transverse momenta of the exchanged two gluons. In the leading logarithmic approximation (LLA) in QCD, this process is related to the off-diagonal gluon density in the proton. As a result, this process may provide a wide window for testing the two-gluon exchange model, and may be particularly useful in studying the small $x$ physics. In comparison, the diffractive bottom jet production is also discussed, and is found to be important only for large transverse momentum of the jet.

preprint1999arXiv

Inelastic $J/ψ$ production in polarized photon-hadron collisions

Presented here is a calculation of inelastic $J/ψ$ production in polarized photon-hadron collisions under the framework of NRQCD factorization formalism. We consider the photoproduction of $\jpsi$ in the energy range relevant to HERA. The Weizsäcker-Williams approximation is adopted in the evaluation of the cross sections for $ep$ collisions. We found that this process can give another independent test for the color-octet mechanism, and the different features for the two color-octet processes may provide further informations on the mechanism for inelastic $\jpsi$ photoproduction. And the discrepancy on the production asymmetry $A$ between various sets of polarized gluon distribution functions is also found to be distinctive.