Source author record

Feng Yuan

Feng Yuan 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

147works
23topics
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

147 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.

preprint2016arXiv

Exploring the accretion model of M87 and 3C 84 with the Faraday rotation measure observations

Low-luminosity active galactic nuclei (LLAGNs) are believed to be powered by an accretion-jet model, consisting of an inner advection-dominated accretion flow (ADAF), an outer truncated standard thin disk, and a jet. But model degeneracy still exists in this framework. For example, the X-ray emission can originate from either the ADAF or jet. The aim of the present work is to check these models with the Faraday rotation measure (RM) observations recently detected for two LLAGNs, M87 and 3C 84, in the sub-mm band. For M87, we find that the RM predicted by the model in which the X-ray emission originates from the ADAF is larger than the observed upper limit of RM by over two orders of magnitude, while the model in which the X-ray emission originates from the jet predicts a RM lower than the observed upper limit. For 3C 84, the sub-mm emission is found to be dominated by the jet component, while the Faraday screen is attributed to the ADAFs. This scenario can naturally explain the observed {\it external} origin of the RM and why RM is found to be stable during a two-year interval although the sub-mm emission increases at the same period.

preprint2016arXiv

Hydrodynamical numerical simulation of wind production from black hole hot accretion flows at very large radii

Previous works show strong winds exist in hot accretion flows around black holes. Those works focus only on the region close to the black hole thus it is unknown whether or where the wind production stops at large radii. In this paper, we investigate this problem by hydrodynamical simulations. We take into account gravities of both the black hole and the nuclear star clusters. For the latter, we assume that the velocity dispersion of stars is a constant and its gravitational potential $\propto σ^2 \ln (r)$, where $σ$ is the velocity dispersion of stars and $r$ is the distance from the center of the galaxy. We focus on the region where the gravitational potential is dominated by the star cluster. We find, same as the accretion flow at small radii, the mass inflow rate decreases inward and the flow is convectively unstable. However, trajectory analysis shows that there is very few wind launched from the flow. Our result, combined with the results of Yuan et al. (2015), indicates that the mass flux of wind launched from hot accretion flow $\dot{M}_{\rm wind}=\dot{M}_{\rm BH}(r/20r_s)$, with $r\la R_A\equiv GM_{\rm BH}/σ^2$. Here $\dot{M}_{\rm BH}$ is accretion rate at black hole horizon. $R_A$ is similar to Bondi radius. We argue that the inward decrease of inflow rate is not because of mass loss via wind, but because of convective motion. The disappearance of wind outside of $R_A$ must be because of the change of the gravitational potential, but the exact reason remains to be probed.

preprint2016arXiv

Interpreting the radio/X-ray correlation of black hole sources based on the accretion-jet model

Two types of correlations between the radio and X-ray luminosities ($L_R$ and $L_X$) have been found in black hole X-ray binaries. For some sources, they follow the `original' type of correlation which is described by a single power-law. Later it was found that some other sources follow a different correlation consisting of three power-law branches, with each branch having different power-law indexes. In this work, we explain these two types of correlation under the coupled accretion--jet model. We attribute the difference between these two types of sources to the different value of viscosity parameter $α$. One possible reason for different $α$ is the different configuration of magnetic field in the accretion material coming from the companion stars. For the `single power-law' sources, their $α$ is high; so their accretion is always in the mode of advection-dominated accretion flow (ADAF) for the whole range of X-ray luminosity. For those `hybrid power-law' sources, the value of $α$ is small so their accretion mode changes from an ADAF to a luminous hot accretion flow, and eventually to two-phase accretion as the accretion rate increases. Because the dependence of radiative efficiency on the mass accretion rate is different for these three accretion modes, different power-law indexes in the $L_R--L_X$ correlation are expected. Constraints on the ratio of the mass loss rate into the jet and the mass accretion rate in the accretion flow are obtained, which can be tested in future by radiative magnetohydrodynamic numerical simulations of jet formation.

preprint2016arXiv

Long Range Correlation in Higgs Boson Plus Two Jets Production at the LHC

We study Higgs boson plus two high energy jets production at the LHC in the kinematics where the two jets are well separated in rapidity. The partonic processes are dominated by the t-channel weak boson fusion (WBF) and gluon fusion (GF) contributions. We derive the associated QCD resummation formalism for the correlation analysis where the total transverse momentum q_\perp of the Higgs boson and two jets is small. Because of different color structures, the resummation results lead to distinguished behaviors: the WBF contribution peaks at relative low q_\perp while all GF channel contributions are strongly de-correlated and spread to a much wider q_\perp range. By applying a kinematic cut on q_\perp, one can effectively increase the WBF signal to the GF background by a significant factor. This greatly strengthens the ability to investigate the WBF channel in Higgs boson production and study the couplings of Higgs to electroweak bosons.

preprint2016arXiv

Magneto-hydrodynamical Numerical simulation of wind production from black hole hot accretion flows at very large radii

Numerical simulations of black hole hot accretion flows have shown the existence of strong wind. Those works focus only on the region close to black hole thus it is unknown whether or where the wind production stops at large radii. To address this question, Bu et al. (2016) have performed hydrodynamic (HD) simulations by taking into account the gravitational potential of both the black hole and the nuclear star clusters. The latter is assumed to be $\propto σ^2 \ln(r)$, with $σ$ being the velocity dispersion of stars and $r$ be the distance from the center of the galaxy. It was found that when the gravity is dominated by nuclear stars, i.e., outside of radius $R_A\equiv GM_{\rm BH}/σ^2$, winds can no longer be produced. That work, however, neglects the magnetic field, which is believed to play a crucial dynamical role in the accretion and thus must be taken into account. In this paper, we revisit this problem by performing magneto-hydrodynamical (MHD) simulations. We confirm the result of Bu et al. (2016), namely wind can't be produced at the region of $R>R_A$. Our result, combined with the results of Yuan et al. (2015), indicates that the formula describing the mass flux of wind $\dot{M}_{\rm wind}=\dot{M}_{\rm BH}(r/20r_s)$ can only be applied to the region where the black hole potential is dominant. Here $\dot{M}_{\rm BH}$ is the mass accretion rate at the black hole horizon and the value of $R_A$ is similar to the Bondi radius.

preprint2016arXiv

On the evolution of a rogue wave along the orthogonal direction of the ($t,x$)-plane

The localization characters of the first-order rogue wave (RW) solution $u$ of the Kundu-Eckhaus equation is studied in this paper. We discover a full process of the evolution for the contour line with height $c^2+d$ along the orthogonal direction of the ($t,x$)-plane for a first-order RW $|u|^2$: A point at height $9c^2$ generates a convex curve for $3c^2\leq d<8c^2$, whereas it becomes a concave curve for $0<d<3c^2$, next it reduces to a hyperbola on asymptotic plane (i.e. equivalently $d=0$), and the two branches of the hyperbola become two separate convex curves when $-c^2<d<0$, and finally they reduce to two separate points at $d=-c^2$. Using the contour line method, the length, width, and area of the RW at height $c^2+d (0<d<8c^2)$ , i.e. above the asymptotic plane, are defined. We study the evolutions of three above-mentioned localization characters on $d$ through analytical and visual methods. The phase difference between the Kundu-Eckhaus and the nonlinear Schrodinger equation is also given by an explicit formula.

preprint2016arXiv

Probing the Small-$x$ Gluon Tomography in Correlated Hard Diffractive Dijet Production in DIS

We investigate the close connection between the quantum phase space Wigner distribution of small-$x$ gluons and the color dipole scattering amplitude, and propose to study it experimentally in the hard diffractive dijet production at the planned electron-ion collider. The angular correlation between the nucleon recoiled momentum and the dijet transverse momentum will probe the nontrivial correlation in the phase space Wigner distribution. This experimental study will not only provide us with three-dimensional tomographic pictures of gluons inside high energy proton, but also give a unique and interesting signal for the small-$x$ dynamics with QCD evolution effects.

preprint2016arXiv

Probing Transverse Momentum Broadening in Heavy Ion Collisions

We study the dijet azimuthal de-correlation in relativistic heavy ion collisions as an important probe of the transverse momentum broadening effects of a high energy jet traversing the quark-gluon plasma. We take into account both the soft gluon radiation in vacuum associated with the Sudakov logarithms and the jet P_T-broadening effects in the QCD medium. We find that the Sudakov effects are dominant at the LHC, while the medium effects can play an important role at RHIC energies. This explains why the LHC experiments have not yet observed sizable P_T-broadening effects in the measurement of dijet azimuthal correlations in heavy ion collisions. Future investigations at RHIC will provide a unique opportunity to study the P_T-broadening effects and help to pin down the underlying mechanism for jet energy loss in a hot and dense medium.

preprint2016arXiv

Single Spin Asymmetry in Forward $pA$ Collisions

We compute the transverse single spin asymmetry in light hadron production $p^\uparrow p\to hX$ and $p^\uparrow A\to hX$ including the gluon saturation effect in the unpolarized nucleon/nucleus. In the forward (large-$x_F$) region, the dominant contribution comes from the so-called derivative term associated with the soft gluonic pole. This leads to the cancellation of nuclear effects in $A_N$ which can be tested at RHIC. We also show that the soft fermionic pole disappears in the saturation environment.

preprint2016arXiv

Sudakov Resummations in Mueller-Navelet Dijet Production

In high energy hadron-hadron collisions, dijet production with large rapidity separation proposed by Mueller and Navelet, is one of the most interesting processes which can help us to directly access the well-known Balitsky-Fadin-Kuraev-Lipatov evolution dynamics. The objective of this work is to study the Sudakov resummation of Mueller-Navelet jets. Through the one-loop calculation, Sudakov type logarithms are obtained for this process when the produced dijets are almost back-to-back. These results could play an important role in the phenomenological study of dijet correlations with large rapidity separation at the LHC.

preprint2016arXiv

Thick Accretion Disk Model for Ultraluminous Supersoft Sources

We propose a geometrically thick, super-Eddington accretion disk model, where an optically thick wind is not necessary, to understand ultraluminous supersoft sources (ULSs). For high mass accretion rates $\dot M \ga 30\dot M_{\rm Edd}$ and not small inclination angles $θ\ga 25^{\circ}$, where $\dot M_{\rm Edd}$ is the Eddington accretion rate, the hard photons from the hot inner region may be shaded by the geometrically thick inner disk, and therefore only the soft photons from the outer thin disk and the outer photosphere of the thick disk can reach the observer. Our model can naturally explain the approximate relation between the typical thermal radius and the thermal temperature, $R_{\rm bb} \propto T_{\rm bb}^{-2}$. Moreover, the thick disk model can unify ULSs and normal ultraluminous X-ray sources, where the different observational characteristics are probably related to the inclination angle and the mass accretion rate. By comparing our model with the optically thick outflow model, we find that less mass accretion rate is required in our model.

preprint2016arXiv

Two-dimensional inflow-wind solution of black hole accretion with an evenly symmetric magnetic field

We solve the two-dimensional magnetohydrodynamic (MHD) equations of black hole accretion with the presence of magnetic field. The field includes a turbulent component, whose role is represented by the viscosity, and a large-scale ordered component. The latter is further assumed to be evenly symmetric with the equatorial plane. The equations are solved in the $r-θ$ plane of a spherical coordinate by assuming time-steady and radially self-similar. An inflow-wind solution is found. Around the equatorial plane, the gas is inflowing; while above and below the equatorial plane at a certain critical $θ$ angle, $θ\sim 47^{\circ}$, the inflow changes its direction of radial motion and becomes wind. The driving forces are analyzed and found to be the centrifugal force and the gradient of gas and magnetic pressure. The properties of wind are also calculated. The specific angular momentum of wind is found to be significantly larger than that of inflow, thus wind can transfer angular momentum outward. These analytical results are compared to those obtained by the trajectory analysis based on MHD numerical simulation data and good agreements are found.

preprint2015arXiv

Detection of a Compact Nuclear Radio Source in the Local Group Elliptical Galaxy M32

The Local Group compact elliptical galaxy M32 hosts one of the nearest candidate super-massive black holes (SMBHs), which has a previously suggested X-ray counterpart. Based on sensitive observations taken with the {\it Karl G. Jansky} Very Large Array (VLA), we detect for the first time a compact radio source coincident with the nucleus of M32, which exhibits an integrated flux density of $\sim$$47.3\pm6.1$ $μ$Jy at 6.6 GHz. We discuss several possibilities for the nature of this source, favoring an origin of the long-sought radio emission from the central SMBH, for which we also revisit the X-ray properties based on recently acquired {\sl Chandra} and {\sl XMM-Newton} data. Our VLA observations also discover radio emission from three previously known optical planetary nebulae in the inner region of M32.

preprint2015arXiv

Extraction of Quark Transversity Distribution and Collins Fragmentation Functions with QCD Evolution

We study the transverse momentum dependent (TMD) evolution of the Collins azimuthal asymmetries in $e^+e^-$ annihilations and semi-inclusive hadron production in deep inelastic scattering (SIDIS) processes. All the relevant coefficients are calculated up to the next-to-leading logarithmic (NLL) order accuracy. By applying the TMD evolution at the approximate NLL order in the Collins-Soper-Sterman (CSS) formalism, we extract transversity distributions for $u$ and $d$ quarks and Collins fragmentation functions from current experimental data by a global analysis of the Collins asymmetries in back-to-back di-hadron productions in $e^+e^-$ annihilations measured by BELLE and BABAR Collaborations and SIDIS data from HERMES, COMPASS, and JLab HALL A experiments. The impact of the evolution effects and the relevant theoretical uncertainties are discussed. We further discuss the TMD interpretation for our results, and illustrate the unpolarized quark distribution, transversity distribution, unpolarized quark fragmentation and Collins fragmentation functions depending on the transverse momentum and the hard momentum scale. We make detailed predictions for future experiments and discuss their impact.

preprint2015arXiv

Implementing the exact kinematical constraint in the saturation formalism

We revisit the issue of the large negative next-to-leading order (NLO) cross section for single inclusive hadron production in $pA$ collisions in the saturation formalism. By implementing the exact kinematical constraint in the modified dipole splitting functions, two additional positive NLO correction terms are obtained. In the asymptotic large $k_\perp$ limit, we analytically show that these two terms become as large as the negative NLO contributions found in our previous calculation. Furthermore, the numerical results demonstrate that the applicable regime of the saturation formalism can be extended to a larger $k_\perp$ window, where the exact matching between the saturation formalism (in the asymptotic $k_\perp$ regime) and the collinear factorization calculations will have to be performed separately. In addition, after significantly improving the numerical accuracy of the NLO correction, we obtain excellent agreement with the LHC and RHIC data for forward hadron productions.

preprint2015arXiv

Magneto-thermal Disk Wind from Protoplanetary Disks

Global evolution and dispersal of protoplanetary disks (PPDs) is governed by disk angular momentum transport and mass-loss processes. Recent numerical studies suggest that angular momentum transport in the inner region of PPDs is largely driven by magnetized disk wind, yet the wind mass-loss rate remains unconstrained. On the other hand, disk mass loss has conventionally been attributed to photoevaporation, where external heating on the disk surface drives a thermal wind. We unify the two scenarios by developing a 1D model of magnetized disk winds with a simple treatment of thermodynamics as a proxy for external heating. The wind properties largely depend on 1) the magnetic field strength at the wind base, characterized by the poloidal Alfvén speed $v_{Ap}$, 2) the sound speed $c_s$ near the wind base, and 3) how rapidly poloidal field lines diverge (achieve $R^{-2}$ scaling). When $v_{Ap}\gg c_s$, corotation is enforced near the wind base, resulting in centrifugal acceleration. Otherwise, the wind is accelerated mainly by the pressure of the toroidal magnetic field. In both cases, the dominant role played by magnetic forces likely yields wind outflow rates that well exceed purely hydrodynamical mechanisms. For typical PPD accretion-rate and wind-launching conditions, we expect $v_{Ap}$ to be comparable to $c_s$ at the wind base. The resulting wind is heavily loaded, with total wind mass loss rate likely reaching a considerable fraction of wind-driven accretion rate. Implications for modeling global disk evolution and planet formation are also discussed.

preprint2015arXiv

Nucleon tensor charge from Collins azimuthal asymmetry measurements

We investigate the nucleon tensor charge from current experiments by a combined analysis of the Collins asymmetries in two hadron production in $e^+e^-$ annihilations and semi-inclusive hadron production in deep inelastic scattering processes. The transverse momentum dependent evolution is taken into account, for the first time, in the global fit of the Collins fragmentation functions and the quark transversity distributions at the approximate next-to-leading logarithmic order. We obtain the nucleon tensor charge contribution from up and down quarks as $δu=+0.30^{+0.12}_{-0.08}$ and $δd=-0.20_{-0.11}^{+0.28}$ at 90\% confidence level for momentum fraction $0.0065 \le x_B \le 0.35$ and $Q^2=$ 10 GeV$^2$.

preprint2015arXiv

Numerical Simulation of Hot Accretion Flows (III): Revisiting wind properties using trajectory approach

Previous MHD simulations have shown that wind must exist in black hole hot accretion flows. In this paper, we continue our study by investigating the detailed properties of wind, such as mass flux and poloidal speed, and the mechanism of wind production. For this aim, we make use of a three dimensional GRMHD simulation of hot accretion flows around a Schwarzschild black hole. The simulation is designed so that the magnetic flux is not accumulated significantly around the black hole. To distinguish real wind from turbulent outflows, we track the trajectories of the virtual Largrangian particles from simulation data. We find two types of real outflows, i.e., a quasi-relativistic jet close to the axis and a sub-relativistic wind subtending a much larger solid angle. Most of the wind originates from the surface layer of the accretion flow. The poloidal wind speed almost remains constant once they are produced, but the flux-weighted wind speed roughly follows $v_{\rm p, wind}(r)\approx 0.25 v_k(r)$. The mass flux of jet is much lower but the speed is much higher, $v_{\rm p,jet}\sim (0.3-0.4) c$. Consequently, both the energy and momentum fluxes of the wind are much larger than those of the jet. We find that the wind is produced and accelerated primarily by the combination of centrifugal force and magnetic pressure gradient, while the jet is mainly accelerated by magnetic pressure gradient. Finally, we find that the wind production efficiency $ε_{\rm wind}\equiv\dot{E}_{\rm wind}/\dot{M}_{\rm BH}c^2\sim 1/1000$, in good agreement with the value required from large-scale galaxy simulations with AGN feedback.

preprint2015arXiv

Scheme dependence and Transverse Momentum Distribution interpretation of Collins-Soper-Sterman resummation

Following an earlier derivation by Catani-de Florian-Grazzini (2000) on the scheme dependence in the Collins-Soper-Sterman (CSS) resummation formalism in hard scattering processes, we investigate the scheme dependence of the Transverse Momentum Distributions (TMDs) and their applications. By adopting a universal $C$-coefficient function associated with the integrated parton distributions, the difference between various TMD schemes can be attributed to a perturbative calculable function depending on the hard momentum scale. We further apply several TMD schemes to the Drell-Yan process of lepton pair production in hadronic collisions, and find that the constrained non-perturbative form factors in different schemes are remarkably consistent with each other and with that of the standard CSS formalism.

preprint2015arXiv

Spin Decomposition of Electron in QED

We perform a systematic study on the spin decomposition of an electron in QED at one-loop order. It is found that the electron orbital angular momentum defined in Jaffe-Manohar and Ji spin sum rules agrees with each other, and the so-called potential angular momentum vanishes at this order. The calculations are performed in both dimensional regularization and Pauli-Villars regularization for the ultraviolet divergences, and they lead to consistent results. We further investigate the calculations in terms of light-front wave functions, and find a missing contribution from the instantaneous interaction in light-front quantization. This clarifies the confusing issues raised recently in the literature on the spin decomposition of an electron, and will help to consolidate the spin physics program for nucleons in QCD.

preprint2015arXiv

Statistics of X-ray flares of Sagittarius A*: evidence for solar-like self-organized criticality phenomenon

X-ray flares have routinely been observed from the supermassive black hole, Sagittarius A$^\star$ (Sgr A$^\star$), at our Galactic center. The nature of these flares remains largely unclear, despite of many theoretical models. In this paper, we study the statistical properties of the Sgr A$^\star$ X-ray flares, by fitting the count rate (CR) distribution and the structure function (SF) of the light curve with a Markov Chain Monte Carlo (MCMC) method. With the 3 million second \textit{Chandra} observations accumulated in the Sgr A$^\star$ X-ray Visionary Project, we construct the theoretical light curves through Monte Carlo simulations. We find that the $2-8$ keV X-ray light curve can be decomposed into a quiescent component with a constant count rate of $\sim6\times10^{-3}~$count s$^{-1}$ and a flare component with a power-law fluence distribution $dN/dE\propto E^{-α_{\rm E}}$ with $α_{\rm E}=1.65\pm0.17$. The duration-fluence correlation can also be modelled as a power-law $T\propto E^{α_{\rm ET}}$ with $α_{\rm ET} < 0.55$ ($95\%$ confidence). These statistical properties are consistent with the theoretical prediction of the self-organized criticality (SOC) system with the spatial dimension $S = 3$. We suggest that the X-ray flares represent plasmoid ejections driven by magnetic reconnection (similar to solar flares) in the accretion flow onto the black hole.

preprint2015arXiv

The Accretion Wind Model of the Fermi Bubbles (II): Radiation

In a previous work, we have shown that the formation of the Fermi bubbles can be due to the interaction between winds launched from the hot accretion flow in Sgr A* and the interstellar medium (ISM). In that work, we focus only on the morphology. In this paper we continue our study by calculating the gamma-ray radiation. Some cosmic ray protons (CRp) and electrons must be contained in the winds, which are likely formed by physical processes such as magnetic reconnection. We have performed MHD simulations to study the spatial distribution of CRp, considering the advection and diffusion of CRp in the presence of magnetic field. We find that a permeated zone is formed just outside of the contact discontinuity between winds and ISM, where the collisions between CRp and thermal nuclei mainly occur. The decay of neutral pions generated in the collisions, combined with the inverse Compton scattering of background soft photons by the secondary leptons generated in the collisions and primary CR electrons can well explain the observed gamma-ray spectral energy distribution. Other features such as the uniform surface brightness along the latitude and the boundary width of the bubbles are also explained. The advantage of this accretion wind model is that the adopted wind properties come from the detailed small scale MHD numerical simulation of accretion flows and the value of mass accretion rate has independent observational evidences. The success of the model suggests that we may seriously consider the possibility that cavities and bubbles observed in other contexts such as galaxy clusters may be formed by winds rather than jets.

preprint2015arXiv

The RHIC SPIN Program: Achievements and Future Opportunities

Time and again, spin has been a key element in the exploration of fundamental physics. Spin-dependent observables have often revealed deficits in the assumed theoretical framework and have led to novel developments and concepts. Spin is exploited in many parity-violating experiments searching for physics beyond the Standard Model or studying the nature of nucleon-nucleon forces. The RHIC spin program plays a special role in this grand scheme: it uses spin to study how a complex many-body system such as the proton arises from the dynamics of QCD. Many exciting results from RHIC spin have emerged to date, most of them from RHIC running after the 2007 Long Range Plan. In this document we present highlights from the RHIC program to date and lay out the roadmap for the significant advances that are possible with future RHIC running.

preprint2015arXiv

Transverse Momentum Resummation for Dijet Correlation in Hadronic Collisions

We study the transverse momentum resummation for dijet correlation in hadron collisions based on the Collins-Soper-Sterman formalism. The complete one-loop calculations are carried out in the collinear factorization framework for the differential cross sections at low imbalance transverse momentum between the two jets. Important cross checks are performed to demonstrate that the soft divergences cancelled out between different diagrams, and in particular, those associated with final state jets. The leading and sub-leading logarithms are identified. All order resummation is derived following the transverse momentum dependent factorization at this order. Its phenomenological applications are also presented.

preprint2015arXiv

Transverse Spin Structure of the Nucleon

We review the current status and future perspectives of theory and experiments of transverse spin phenomena in high-energy scattering processes off nucleon targets and related issues in nucleon structure and QCD. Systematic exploration of transverse spin effects requires measurements in polarized deep-inelastic scattering, polarized pp collisions, and e+e- annihilations. Sophisticated QCD-based techniques are also needed to analyze the experimental data sets.

preprint2015arXiv

Universal Non-perturbative Functions for SIDIS and Drell-Yan Processes

We update the well-known BLNY fit to the low transverse momentum Drell-Yan lepton pair productions in hadronic collisions, by considering the constraints from the semi-inclusive hadron production in deep inelastic scattering (SIDIS) from HERMES and COMPASS experiments. We follow the Collins-Soper-Sterman (CSS) formalism with the b_*-prescription. A universal non-perturbative form factor associated with the transverse momentum dependent quark distributions is found in the analysis with a new functional form different from that of BLNY. This releases the tension between the BLNY fit to the Drell-Yan data with the SIDIS data from HERMES/COMPASS in the CSS resummation formalism.

preprint2014arXiv

AGN feedback in an isolated elliptical galaxy: the effect of strong radiative feedback in the kinetic mode

Based on two-dimensional high resolution hydrodynamic numerical simulation, we study the mechanical and radiative feedback effects from the central AGN on the cosmological evolution of an isolated elliptical galaxy. Physical processes such as star formation and supernovae are considered. The inner boundary of the simulation domain is carefully chosen so that the fiducial Bondi radius is resolved and the accretion rate of the black hole is determined self-consistently. In analogy to previous works, we assume that the specific angular momentum of the galaxy is low. It is well-known that when the accretion rates are high and low, the central AGNs will be in cold and hot accretion modes, which correspond to the radiative and kinetic feedback modes, respectively. The emitted spectrum from the hot accretion flows is harder than that from the cold accretion flows, which could result in a higher Compton temperature accompanied by a more efficient radiative heating, according to previous theoretical works. Such a difference of the Compton temperature between the two feedback modes, the focus of this study, has been neglected in previous works. Significant differences in the kinetic feedback mode are found as a result of the stronger Compton heating and accretion becomes more chaotic. More importantly, if we constrain models to correctly predict black hole growth and AGN duty cycle after cosmological evolution, we find that the favored model parameters are constrained: mechanical feedback efficiency diminishes with decreasing luminosity (the maximum efficiency being $\simeq 10^{-3.5}$) and X-ray Compton temperature increases with decreasing luminosity, although models with fixed mechanical efficiency and Compton temperature can be found that are satisfactory as well. We conclude that radiative feedback in the kinetic mode is much more important than previously thought.

preprint2014arXiv

Comment on "Next-to-leading order forward hadron production in the small-x regime: rapidity factorization" arXiv:1403.5221 by Kang et al

In a recent paper (arXiv:1403.5221), Kang et al.proposed the so-called "rapidity factorization" for the single inclusive forward hadron production in pA collisions. We point out that the leading small-x logarithm was mis-identified in this paper, and hence the newly added next-to-leading order correction term is unjustified and should be absent in view of the small-x factorization.

preprint2014arXiv

Confronting the jet model of Sgr A* with the Faraday rotation measure observations

Sgr A* is probably the supermassive black hole being investigated most extensively due to its proximity. Several theoretical models for its steady state emission have been proposed in the past two decades. Both the radiative-inefficient accretion flow and the jet model have been shown to well explain the observed spectral energy distribution. Faraday rotation measure (RM) has been unambiguously measured at submillimeter wavelength, but has only been tested against the accretion flow model. Here we first calculate the RM based on the jet model and find that the predicted value is two orders of magnitude lower than the measured value. We then include an additional contribution from the accretion flow in front of the jet and show that the measured RM may be reconciled with the model under some tight constraints. The main constraint is that the inclination angle should be greater than $\sim 73^{\circ}$. But this requirement is not consistent with an existing observational estimate of the inclination angle.

preprint2014arXiv

Correlation between the photon index and X-ray luminosity of black hole X-ray binaries and active galactic nuclei: observations and interpretation

We investigate the observed correlation between the 2--10 keV X-ray luminosity (in unit of the Eddington luminosity; $l_X \equiv L_X/L_{Edd}$) and the photon index ($Γ$) of the X-ray spectrum for both black hole X-ray binaries (BHBs) and active galactic nuclei (AGNs). We construct a large sample, with $10^{-9} < l_X < 10^{-1}$. We find that $Γ$ is positively and negatively correlated with $l_X$ when $l_X > 10^{-3}$ and $10^{-6.5} < l_X < 10^{-3}$ respectively, while $Γ$ is nearly a constant when $l_X < 10^{-6.5}$. We explain the above correlation in the framework of a coupled hot accretion flow -- jet model. The radio emission always come from the jet while the X-ray emission comes from the accretion flow and jet when $l_X$ is above and below $10^{-6.5}$, respectively. More specifically, we assume that with the increase of mass accretion rate, the hot accretion flow develops into a clumpy and further a disc -- corona two-phase structure because of thermal instability. We argue that such kind of two-phase accretion flow can explain the observed positive correlation.

preprint2014arXiv

Does the circularization radius exist or not for low angular momentum accretion?

If the specific angular momentum of accretion gas at large radius is small compared to the local Keplerian value, one usually believes that there exists a "circularization radius" beyond which the angular momentum of accretion flow is almost a constant while within which a disk is formed and the angular momentum roughly follows the Keplerian distribution. In this paper, we perform numerical simulations to study whether the picture above is correct in the context of hot accretion flow. We find that for a steady accretion flow, the "circularization radius" does not exist and the angular momentum profile will be smooth throughout the flow. However, for transient accretion systems, such as the tidal disruption of a star by a black hole, a "turning point" should exist in the radial profile of the angular momentum, which is conceptually similar to the "circularization radius". At this radius, the viscous timescale equals the life time of the accretion event. The specific angular momentum is close to Keplerian within this radius, while beyond this radius the angular momentum is roughly constant.

preprint2014arXiv

Fermi Bubbles Inflated by Winds Launched from the Hot Accretion Flow in Sgr A*

A pair of giant gamma-ray bubbles have been revealed by the {\it Fermi} LAT. In this paper we investigate their formation mechanism. Observations have indicated that the activity of the supermassive black hole located at the Galactic center, Sgr A*, was much stronger than the present time. Specifically, one possibility is that while Sgr A* was also in the hot accretion regime, the accretion rate should be $10^3-10^4$ times higher during the past $\sim 10^7$ yr. On the other hand, recent MHD numerical simulations of hot accretion flows have unambiguously shown the existence of strong winds and obtained their properties. Based on these knowledge, by performing three-dimensional hydrodynamical simulations, we show in this paper that the Fermi bubbles could be inflated by winds launched from the ``past' hot accretion flow in Sgr A*. In our model, the active phase of Sgr A* is required to last for about 10 million years and it was quenched no more than 0.2 million years ago. The Central Molecular Zone (CMZ) is included and it collimates the wind orientation towards the Galactic poles. Viscosity suppresses the Rayleigh-Taylor and Kelvin-Helmholtz instabilities and results in the smoothness of the bubble edge. The main observational features of the bubbles can be well explained. Specifically, the {\it ROSAT} X-ray features are interpreted by the shocked interstellar medium and the interaction region between winds and CMZ gas. The thermal pressure and temperature obtained in our model are in good consistency with the recent {\it Suzaku} observations.

preprint2014arXiv

Hot Accretion Flows Around Black Holes

Black hole accretion flows can be divided into two broad classes: cold and hot. Cold accretion flows, which consist of cool optically thick gas, are found at relatively high mass accretion rates. Prominent examples are the standard thin disk, which occurs at a fraction of the Eddington mass accretion rate, and the slim disk at super-Eddington rates. These accretion flows are responsible for luminous systems such as active galactic nuclei radiating at or close to the Eddington luminosity and black hole X-ray binaries in the soft state. Hot accretion flows, the topic of this review, are virially hot and optically thin. They occur at lower mass accretion rates, and are described by models such as the advection-dominated accretion flow and luminous hot accretion flow. Because of energy advection, the radiative efficiency of these flows is in general lower than that of a standard thin accretion disk. Moreover, the efficiency decreases with decreasing mass accretion rate. Observations show that hot accretion flows are associated with jets. In addition, theoretical arguments suggest that hot flows should produce strong winds. Hot accretion flows are believed to be present in low-luminosity active galactic nuclei and in black hole X-ray binaries in the hard and quiescent states. The prototype is Sgr A*, the ultra-low-luminosity supermassive black hole at our Galactic center. The jet, wind and radiation from a supermassive black hole with a hot accretion flow can interact with the external interstellar medium and modify the evolution of the host galaxy. Details of this "maintenance-mode feedback" could, in principle, be worked out through theoretical studies and numerical simulations of hot accretion flows.

preprint2014arXiv

Kinematical correlations for Higgs boson plus High Pt Jet Production at Hadron Colliders

We investigate the effect of QCD resummation to kinematical correlations in the Higgs boson plus high transverse momentum (Pt) jet events produced at hadron colliders. We show that at the complete one-loop order, the Collins-Soper-Sterman resummation formalism can be applied to derive the Sudakov form factor, which is found to be independent of jet-finding algorithm. We compare the singular behavior of resummation calculation to fixed order prediction in the case that Higgs boson and high Pt jet are produced nearly back-to-back in their transverse momenta, and find a perfect agreement. The phenomenological importance of the resummation effect at the LHC is also demonstrated.

preprint2014arXiv

Matching Collinear and Small-x Factorization Calculations for Inclusive Hadron Production in pA Collisions

We construct a theoretical framework to match the formulas for forward inclusive hadron productions in pA collisions in the small-x saturation formalism and collinear factorization. The small-x calculation can be viewed as a power series in $Q_s^2/k_\perp^2$, in which the collinear factorization result corresponds to the leading term. At high transverse momentum, the subleading correction terms are insignficant, whereas at low $p_\perp$, the power corrections become important and the small-x resummation is essential to describe the differential cross section. We show that the familiar collinear factorization calculation can smoothly match the results from small-x factorization at next-to-leading order in $α_s$ when we use exact kinematics, as opposed to the approximate kinematics in previous work. With this matching, we can describe the experimental data from RHIC very well at high $p_\perp$.

preprint2014arXiv

Soft Factor Subtraction and Transverse Momentum Dependent Parton Distributions on Lattice

We study the transverse momentum dependent (TMD) parton distributions in the newly proposed quasi-parton distribution function framework in Euclidean space. A soft factor subtraction is found to be essential to make the TMDs calculable on lattice. We show that the quasi-TMDs with the associated soft factor subtraction can be applied in hard QCD scattering processes such as Drell-Yan lepton pair production in hadronic collisions. This allows future lattice calculations to provide information on the non-perturbative inputs and energy evolutions for the TMDs. Extension to the generalized parton distributions and quantum phase space Wigner distributions will lead to a complete nucleon tomography on lattice.

preprint2014arXiv

Soft Gluon Resummations in Dijet Azimuthal Angular Correlations at the Collider

We derive all order soft gluon resummation in dijet azimuthal angular correlation in hadronic collisions at the next-to-leading logarithmic level. The relevant coefficients for the Sudakov resummation factor, the soft and hard factors, are calculated. The theory predictions agree well with the experimental data from D0 Collaboration at the Tevatron. This provides a benchmark calculation for the transverse momentum dependent QCD resummation for jet productions in hadron collisions and can be readily applied at the CERN LHC.

preprint2014arXiv

Two dimensional numerical simulations of Supercritical Accretion Flows revisited

We study the dynamics of super-Eddington accretion flows by performing two-dimensional radiation-hydrodynamic simulations. Compared with previous works, in this paper we include the $T_{θϕ}$ component of the viscous stress and consider various values of the viscous parameter $α$. We find that when $T_{θϕ}$ is included, the rotational speed of the high-latitude flow decreases, while the density increases and decreases at the high and low latitudes, respectively. We calculate the radial profiles of inflow and outflow rates. We find that the inflow rate decreases inward, following a power law form of $\dot{M}_{\rm in}\propto r^s$. The value of $s$ depends on the magnitude of $α$ and is within the range of $\sim 0.4-1.0$. Correspondingly, the radial profile of density becomes flatter compared with the case of a constant $\dot{M}(r)$. We find that the density profile can be described by $ρ(r)\propto r^{-p}$, and the value of $p$ is almost same for a wide range of $α$ ranging from $α=0.1$ to $0.005$. The inward decrease of inflow accretion rate is very similar to hot accretion flows, which is attributed to the mass loss in outflows. To study the origin of outflow, we analyze the convective stability of slim disk. We find that depending on the value of $α$, the flow is marginally stable (when $α$ is small) or unstable (when $α$ is large). This is different from the case of hydrodynamical hot accretion flow where radiation is dynamically unimportant and the flow is always convectively unstable. We speculate that the reason for the difference is because radiation can stabilize convection. The origin of outflow is thus likely because of the joint function of convection and radiation, but further investigation is required.

preprint2013arXiv

Electron correlation and impurity-induced quasiparticle resonance states in cuprate superconductors

We theoretically study the quasiparticle resonance states around a nonmagnetic impurity in cuprate superconductors based on t-t'-J-U model. The purpose of introducing the Coulomb repulsive interaction U is to partially impose the double occupancy constraint by employing the Gutzwiller projected mean-field approximation. We determine the spatial variation of the order parameter and the local density of states (LDOS) by self-consistently solving Bogoliubov-de Gennes equations. We find that in the large U limit, a zero-energy resonance peak in the LDOS indeed appears for the impurity potential in the unitary limit, at the same time the asymmetric superconducting coherence peaks are strongly suppressed. As U decreases the electron double occupancy d is permitted and gradually increases, leading to the decreasing of order parameter. In particular the above zero-energy resonance peak begins to evolve into a double-peaked structure since a critical value $d_{c}$. These important feathers are qualitatively agreement with the scanning tunneling spectroscopy experiments, and uncover the essential role played by the electron correlation in cuprate superconductors.

preprint2013arXiv

Energy Evolution for the Sivers Asymmetries in Hard Processes

We investigate the energy evolution of the azimuthal spin asymmetries in semi-inclusive hadron production in deep inelastic scattering (SIDIS) and Drell-Yan lepton pair production in pp collisions. The scale dependence is evaluated by applying an approximate solution to the Collins-Soper-Sterman (CSS) evolution equation at one-loop order which is adequate for moderate Q^2 variations. This describes well the unpolarized cross sections for SIDIS and Drell-Yan process in the $Q^2$ range of 2.4-100GeV^2. A combined analysis of the Sivers asymmetries in SIDIS from HERMES and COMPASS experiments, and the predictions for the Drell-Yan process at RHIC at \sqrt{S}=200GeV are presented. We further extend to the Collins asymmetries and find, for the first time, a consistent description for HERMES/COMPASS and BELLE experiments with the evolution effects. We emphasize an important test of the evolution effects by studying di-hadron azimuthal asymmetry in e^+e^- annihilation at moderate energy range, such as at BEPC at \sqrt{S}=4.6GeV.

preprint2013arXiv

Global Neutrino Heating in Hyperaccretion Flows

The neutrino-dominated accretion flow (NDAF) with accretion rates \dot{M} = 0.01 - 10 M_{\sun} s^{-1} is a plausible candidate for the central engine of gamma-ray bursts (GRBs). This hyperaccretion disk is optically thin to neutrinos in the radial direction, therefore the neutrinos produced at one radius can travel for a long distance in the disk. Those neutrinos can thus be absorbed with certain probability by the disk matter at the other radius and heat the disk there. The effect of this "global neutrino heating" has been ignored in previous works and is the focus of this paper. We find that around the "ignition" radius r_{ign}, the global neutrino heating rate could be comparable to or even larger than the local viscous heating rate thus must be an important process. Two possible consequences are in order if the "global neutrino heating" is taken into account: i) the temperature of the disk is slightly raised and the "ignition" radius r_{ign} slightly shifts to a larger radius, both lead to the increasing of the total neutrino flux; ii) what is more interesting is that, the temperature of the ADAF just beyond r_{ign} may be raised above the virial temperature thus the accretion will be suppressed. In this case, the activity of the black hole is expected to oscillate between an active and inactive phases. The timescale of the active phases is estimated to be \sim 1 second. If the timescale of the inactive phase is comparable to or less than this value, this intermittent activity may explain the slow variability component of the GRBs. Self-consistent global calculations of NDAFs with the "global neutrino heating" included are required in the future to more precisely evaluate this effect.

preprint2013arXiv

Heavy Quarkonium Production at Low Pt in NRQCD with Soft Gluon Resummation

We extend the non-relativistic QCD (NRQCD) prediction for the production of heavy quarkonium with low transverse momentum in hadronic collisions by taking into account effects from all order soft gluon resummation. Following the Collins-Soper-Sterman formalism, we resum the most singular terms in the partonic subprocesses. The theoretical predictions of $J/ψ$ and $Υ$ productions are compared to the experimental data from the fixed target experiments (E866) and the collider experiments (RHIC, Tevatron, LHC). The associated non-perturbative Sudakov form factor for the gluon distributions is found to be different from the previous assumption of rescaling the quark form factor by the ratio of color factors. This conclusion should be further checked by future experiments on Higgs boson and/or di-photon production in $pp$ collisions. We also comment on the implication of our results on determining the color-octet matrix elements associated with the $J/ψ$ and $Υ$ productions in the NRQCD factorization formalism.

preprint2013arXiv

Low Transverse Momentum Heavy Quark Pair Production to Probe Gluon Tomography

We derive the transverse momentum dependent (TMD) factorization for heavy quark pair production in deep inelastic scattering, where the total transverse momentum is much smaller than the invariant mass of the pair. The factorization is demonstrated at one-loop order, in both Ji-Ma-Yuan and Collins-11 schemes for the TMD definitions, and the hard factors are calculated accordingly. Our result provides a solid theoretical foundation for the phenomenological investigations of the gluon TMDs in this process, and can be extended to other similar hard processes, including dijet (di-hadron) production in DIS.

preprint2013arXiv

On the role of initial and boundary conditions in numerical simulations of accretion flows

We study the effects of initial and boundary conditions, taking two-dimensional hydrodynamical numerical simulations of hot accretion flow as an example. The initial conditions considered include a rotating torus, a solution expanded from the one-dimensional global solution of hot accretion flows, injected gas with various angular momentum distributions, and the gas from a large-scale numerical simulation. Special attention is paid to the radial profiles of the mass accretion rate and density. Both can be described by a power-law function, $\dot{M}\propto r^s$ and $ρ\propto r^{-p}$. We find that if the angular momentum is not very low, the value of $s$ is not sensitive to the initial condition and lies within a narrow range, $0.47\la s \la 0.55$. However, the value of $p$ is more sensitive to the initial condition and lies in the range $0.48\la p \la 0.8$. The diversity of the density profile is because different initial conditions give different radial profiles of radial velocity due to the different angular momentum of the initial conditions. When the angular momentum of the accretion flow is very low, the inflow rate is constant with radius. Taking the torus model as an example, we have also investigated the effects of inner and outer boundary conditions by considering the widely adopted "outflow" boundary condition and the "mass flux conservation" condition. We find that the results are not sensitive to these two boundary conditions.

preprint2013arXiv

Radiation driven outflow in active galactic nuclei: the feedback effects of scattered and reprocessed photons

We perform time-dependent, 2DHD numerical simulations to study the dynamics of a slowly rotating accretion flow from sub-pc to pc scales under the irradiation from the central AGN. Compared to previous work, we improve the calculation of the radiative force due to X-rays. More importantly, in addition to radiative pressure and radiative heating/cooling directly from the central AGN, in the momentum equation we also include the force due to the scattered and reprocessed photons. We find that the accretion flow properties change significantly due to this "re-radiation" effect. The inflow rate at the inner boundary is reduced, while the outflow rate at the outer boundary is enhanced by about one order of magnitude. This effect is more significant when the density at the outer boundary is higher. The properties of outflows such as velocity, momentum and energy fluxes, and the ratio of outflow rate and the accretion rate, are calculated. We find that the efficiency of transferring the radiation power into the kinetic power of outflow is typically $10^{-3}$, far below the value of $\sim 0.05$ which is assumed in some cosmological simulations. The effect of the temperature of the gas at the outer boundary ($T_0$) is investigated. When $T_0$ is high, the emitted luminosity of the accretion flow oscillates. This is because in this case the gas around the Bondi radius can be more easily heated to be above the virial temperature due to its high internal energy. Another question we hope to address is the so-called "sub-Eddington" puzzle. Observationally, the luminosity of almost all AGNs are sub-Eddington, while theoretically the luminosity of an accretion flow can easily be super-Eddington. We find that even when the re-radiation effect is included and outflow does become much stronger, the luminosity, while reduced, can still be super-Eddington.

preprint2013arXiv

Sudakov Double Logarithms Resummation in Hard Processes in Small-x Saturation Formalism

In this manuscript, we present a complete study of the Sudakov double logarithms resummation for various hard processes in eA and pA collisions in the small-x saturation formalism. We first employ a couple of slightly different formalisms to perform the one-loop analysis of the Higgs boson production process in pA collisions, and demonstrate that Sudakov-type logarithms arise as the leading correction and that they can be systematically resummed in addition to the usual small-x resummation. We further study the Sudakov double logarithms for other processes such as heavy quark pair production and back-to-back dijet production in eA and pA collisions through detailed calculation of the corresponding one-loop diagrams. As the most important contribution from the one-loop correction, the Sudakov factor should play an important role in the phenomenological study of saturation physics in the pA programs at RHIC and the LHC.

preprint2013arXiv

The RHIC Spin Program: Achievements and Future Opportunities

This document summarizes recent achievements of the RHIC spin program and their impact on our understanding of the nucleon's spin structure, i.e. the individual parton (quark and gluon) contributions to the helicity structure of the nucleon and to understand the origin of the transverse spin phenomena. Open questions are identified and a suite of future measurements with polarized beams at RHIC to address them is laid out. Machine and detector requirements and upgrades are briefly discussed.

preprint2013arXiv

TMD Evolution: Matching SIDIS to Drell-Yan and W/Z Boson Production

We examine the QCD evolution for the transverse momentum dependent observables in hard processes of semi-inclusive hadron production in deep inelastic scattering and Drell-Yan lepton pair production in $pp$ collisions, including the spin-average cross sections and Sivers single transverse spin asymmetries. We show that the evolution equations derived by a direct integral of the Collins-Soper-Sterman evolution kernel from low to high Q can describe well the transverse momentum distribution of the unpolarized cross sections in the Q^2 range from 2 to 100 GeV^2. In addition, the matching is established between our evolution and the Collins-Soper-Sterman resummation with b*-prescription and Konychev-Nodalsky parameterization of the non-perturbative form factors, which are formulated to describe the Drell-Yan lepton pair and W/Z boson production in hadronic collisions. With these results, we present the predictions for the Sivers single transverse spin asymmetries in Drell-Yan lepton pair production and $W^\pm$ boson production in polarized pp and $π^- p$ collisions for several proposed experiments. We emphasize that these experiments will not only provide crucial test of the sign change of the Sivers asymmetry, but also provide important opportunities to study the QCD evolution effects.

preprint2013arXiv

Universal Suppression of Heavy Quarkonium Production in pA Collisions at Low Transverse Momentum

The nuclear suppression of heavy quarkonium production at low transverse momentum in pA collisions in high energy scatterings is investigated in the small-x factorization formalism. A universal suppression is found in the large Nc limit between the two formalisms to describe the heavy quarkonium production: the non-relativistic QCD (NRQCD) and the color-evaporation model (CEM). This provides an important probe to the saturation momentum at small-x in big nucleus. We also comment on the phenomenological applications of our results.

preprint2012arXiv

Back-to-Back Correlations of Di-hadrons in dAu Collisions at RHIC

We perform a complete theoretical analysis of the azimuthal angular correlation of two-hadron productions in the forward $dAu$ collisions at RHIC in the saturation formalism, and obtain a very good agreement with the experimental data. It is demonstrated that the suppression and broadening of the away side peak provide a unique signal for the onset of the saturation mechanism at small-$x$ in a large nucleus. We emphasize that future experiments of di-hadron correlations in $pA$ collisions at both RHIC and LHC, and in $eA$ collisions at the planned electron-ion collider, shall provide us with a thorough study and understanding of the strong interaction dynamics in the saturation regime.

preprint2012arXiv

Episodic Jets as the Central Engine of Gamma-Ray Bursts

Most Gamma-ray bursts (GRBs) have erratic light curves, which demand that the GRB central engine launches an episodic outflow. Recent Fermi observations of some GRBs indicate a lack of the thermal photosphere component as predicted by the baryonic fireball model, which suggests a magnetic origin of GRBs. In view that powerful episodic jets have been observed along with continuous jets in other astrophysical black hole systems, here we propose an intrinsically episodic, magnetically-dominated jet model for GRB central engine. Accumulation and eruption of free magnetic energy in the corona of a differentially-rotating, turbulent accretion flow around a hyperaccreting black hole lead to ejections of episodic, magnetically dominated plasma blobs. These blobs are accelerated magnetically, collide with each other at large radii, trigger rapid magnetic reconnection and turbulence, efficient particle acceleration and radiation, and power the observed episodic prompt gamma-ray emission from GRBs.

preprint2012arXiv

Inclusive Hadron Productions in pA Collisions

We calculate inclusive hadron productions in pA collisions in the small-x saturation formalism at one-loop order. The differential cross section is written into a factorization form in the coordinate space at the next-to-leading order, while the naive form of the convolution in the transverse momentum space does not hold. The rapidity divergence with small-x dipole gluon distribution of the nucleus is factorized into the energy evolution of the dipole gluon distribution function, which is known as the Balitsky-Kovchegov equation. Furthermore, the collinear divergences associated with the incoming parton distribution of the nucleon and the outgoing fragmentation function of the final state hadron are factorized into the splittings of the associated parton distribution and fragmentation functions, which allows us to reproduce the well-known DGLAP equation. The hard coefficient function, which is finite and free of divergence of any kind, is evaluated at one-loop order.

preprint2012arXiv

NLO QCD Corrections to $B_c$-to-Charmonium Form Factors

The $B_c(^1S_0)$ meson to S-wave Charmonia transition form factors are calculated in next-to-leading order(NLO) accuracy of Quantum Chromodynamics(QCD). Our results indicate that the higher order corrections to these form factors are remarkable, and hence are important to the phenomenological study of the corresponding processes. For the convenience of comparison and use, the relevant expressions in asymptotic form at the limit of $m_c\rightarrow0$ for the radiative corrections are presented.

preprint2012arXiv

Numerical Simulation of Hot Accretion Flows (I): A Large Radial Dynamical Range and the Density Profile of Accretion Flow

Numerical simulations of hot accretion flow have shown that the mass accretion rate decreases with decreasing radius; consequently the density profile of accretion flow becomes flatter compared to the case of a constant accretion rate. This result has important theoretical and observational implications. However, because of technical difficulties, the radial dynamic range in almost all previous simulations usually spans at most two orders of magnitude. This small dynamical range, combined with the effects of boundary conditions, makes the simulation results suspectable. Especially, the radial profiles of density and accretion rate may not be precise enough to be used to compare with observations. In this paper we present a "two-zone" approach to expand the radial dynamical range from two to four orders of magnitude. We confirm previous results and find that from $r_s$ to $ 10^4r_s$ the radial profiles of accretion rate and density can be well described by $\dot{M}(r)\propto r^s$ and $ρ\propto r^{-p}$. The values of (s, p) are (0.48, 0.65) and (0.4, 0.85), for viscous parameter $α=0.001$ and 0.01, respectively. We have looked up numerical simulation works in the literature and found that the values of $s$ and $p$ are all similar, no matter a magnetic field is included or not and what kind of initial conditions are adopted. The density profile we obtain is in good quantitative agreement with that obtained from the detailed observations and modeling to Sgr A* and NGC 3115. The origin of such a accretion rate profile will be investigated in a subsequent paper.

preprint2012arXiv

Numerical Simulation of Hot Accretion Flows (II): Nature, Origin, and Properties of Outflow and Their Possible Observational Applications

Previous hydrodynamical (HD) and magnetohydrodynamical (MHD) numerical simulations of hot accretion flows have shown that the mass accretion rate decreases with decreasing radius. Two models have been proposed to explain this result. In the ADIOS model, the inward decrease of accretion rate is because of the loss of gas in the outflow. In the CDAF model, the gas is assumed to be locked in convective eddies, which results in the inward decrease of the accretion rate. We investigate the nature of inward decrease of accretion rate using HD and MHD simulations. We calculate various properties of inflow and outflow, including the mass flux, radial and rotational velocities, temperature, and the Bernoulli parameter ($Be$). Systematic and significant differences between inflow and outflow are found. These results suggest that the inflow and outflow are not dominated by convective turbulence, but are systematic inward and outward motion. We have also analyzed the convective stability of MHD accretion flow and found that they are convectively stable. These results indicate that the ADIOS scenario is favored. The different properties of inflow and outflow also suggest that the mechanisms of producing outflow in HD and MHD flows are buoyancy and centrifugal force associated with the magnetic field, respectively. The latter mechanism is similar to the Blandford & Payne mechanism. We also study the effect of initial conditions in the simulations. We find that the value of $Be$ is mainly determined by the value of $Be$ of the initial condition. We discuss some possible observational applications of our outflow model. These observations include the Fermi bubble observed in the Galaxy center, and winds widely observed in AGNs and black hole X-ray binaries.

preprint2012arXiv

On the linearly polarized gluon distributions in the color dipole model

We show that the linearly polarized gluon distributions appear in the color dipole model as we derive the full cross sections of the DIS dijet production and the Drell-Yan dijet ($γ^*$ jet correlation) process. Together with the normal Weizsäcker-Williams gluon distribution, the linearly polarized one will contribute to the DIS dijet production cross section as the coefficient of the $\cos(2 Δϕ) $ term in the correlation limit. We also derive the exact results for the cross section of the Drell-Yan dijet process, and find that the linearly polarized dipole gluon distribution which is identical to the normal dipole gluon distribution involves in the cross section. The results obtained in this paper agree with the previous transverse momentum dependent factorization study. We further derive the small-$x$ evolution of these linearly polarized gluon distributions and find that they rise as $x$ gets small at high energy.

preprint2012arXiv

Outflow from Hot Accretion Flows

Numerical simulations of hot accretion flows have shown that the mass accretion rate decreases with decreasing radius. Two models have been proposed to explain this result. In the adiabatic inflow-outflow solution (ADIOS), it is thought to be due to the loss of gas in outflows. In the convection-dominated accretion flow (CDAF) model, it is explained as because that the gas is locked in convective eddies. In this paper we use hydrodynamical (HD) and magnetohydrodynamical (MHD) simulations to investigate which one is physical. We calculate and compare various properties of inflow (gas with an inward velocity) and outflow (gas with an outward velocity). Systematic and significant differences are found. For example, for HD flows, the temperature of outflow is higher than inflow; while for MHD flows, the specific angular momentum of outflow is much higher than inflow. We have also analyzed the convective stability of MHD accretion flow and found that they are stable. These results suggest that systematic inward and outward motion must exist, i.e., the ADIOS model is favored. The different properties of inflow and outflow also suggest that the mechanisms of producing outflow in HD and MHD flows are buoyancy associated with the convection and the centrifugal force associated with the angular momentum transport mediated by the magnetic field, respectively. The latter mechanism is similar to the Blandford & Payne mechanism but no large-scale open magnetic field is required here. Possible observational applications are briefly discussed.

preprint2012arXiv

Probing Parton Orbital Angular Momentum in Longitudinally Polarized Nucleon

While the total orbital angular momentum (OAM) of a definite quark flavor in a longitudinally-polarized nucleon can be obtained through a sum rule involving twist-two generalized parton distribution (GPDs), its distribution as a function of parton momentum in light-front coordinates is more complicated to define and measure because it involves intrinsically twist-three effects. In this paper, we consider two different parton OAM distributions. The first is manifestly gauge invariant, and its moments are local operators and calculable in lattice QCD. We show that it can potentially be measured through twist-three GPDs. The second is the much-debated canonical OAM distribution natural in free-field theory and light-cone gauge. We show the latter in light-cone gauge can also be related to twist-three GPDs as well as quantum phase-space Wigner distributions, both being measurable in high-energy experiments.

preprint2012arXiv

Proton Spin Structure from Measurable Parton Distributions

We present a systematic study of the proton spin structure in terms of measurable parton distributions. For a transversely-polarizedproton, we derive a polarization sum rule from the leading generalized parton distributions appearing in hard exclusive processes. For a longitudinally-polarized proton, we obtain a helicity decomposition from well-known quark and gluon helicity distributions and orbital angular-momentum contributions. The latter is shown to be related to measurable subleading generalized parton distributions and quantum-phase space Wigner distributions.

preprint2012arXiv

Sudakov Resummation in Small-x Saturation Formalism

Through an explicit calculation of massive scalar particle (e.g., Higgs boson) production in high energy pA collisions up to one-loop order, we demonstrate that the Sudakov-type logarithms in hard processes in small-x saturation formalism can be systematically separated from the typical small-x logarithms. The generic feature of the Sudakov logarithms and all order resummation is derived. We further comment on the phenomenological implications and extension to other hard processes in the small-x calculations.

preprint2012arXiv

The effects of viscosity on the circumplanetary disks

The effects of viscosity on the circumplanetary disks residing in the vicinity of protoplanets are investigated through two-dimensional hydrodynamical simulations with the shearing sheet model. We find that viscosity can affect properties of the circumplanetary disk considerably when the mass of the protoplanet is $M_p \lesssim 33M_\oplus$, where $M_\oplus$ is the Earth mass. However, effects of viscosity on the circumplanetary disk are negligibly small when the mass of the protoplanet $M_p \gtrsim 33M_\oplus$. We find that when $M_p \lesssim 33M_\oplus$, viscosity can disrupt the spiral structure of the gas around the planet considerably and make the gas smoothly distributed, which makes the torques exerted on the protoplanet weaker. Thus, viscosity can make the migration speed of a protoplanet lower. After including viscosity, size of the circumplanetary disk can be decreased by a factor of $\gtrsim 20%$. Viscosity helps to transport gas into the circumplanetary disk from the differentially rotating circumstellar disk. The mass of the circumplanetary disk can be increased by a factor of 50% after viscosity is taken into account when $M_p \lesssim 33M_\oplus$. Effects of viscosity on the formation of planets and satellites are briefly discussed.

preprint2012arXiv

The gamma-ray emitting region of the jet in Cyg X-3

We study models of the gamma-ray emission of Cyg X-3 observed by Fermi. We calculate the average X-ray spectrum during the gamma-ray active periods. Then, we calculate spectra from Compton scattering of a photon beam into a given direction by isotropic relativistic electrons with a power-law distribution, both based on the Klein-Nishina cross section and in the Thomson limit. Applying the results to scattering of stellar blackbody radiation in the inner jet of Cyg X-3, we find that a low-energy break in the electron distribution at a Lorentz factor of ~ 300--1000 is required by the shape of the observed X-ray/gamma-ray spectrum in order to avoid overproducing the observed X-ray flux. The electrons giving rise to the observed \g-rays are efficiently cooled by Compton scattering, and the power-law index of the acceleration process is ~ 2.5--3. The bulk Lorentz factor of the jet and the kinetic power before the dissipation region depend on the fraction of the dissipation power supplied to the electrons; if it is ~ 1/2, the Lorentz factor is ~ 2.5, and the kinetic power is ~ 10^38 erg/s, which represents a firm lower limit on the jet power, and is comparable to the bolometric luminosity of Cyg X-3. Most of the power supplied to the electrons is radiated. The broad band spectrum constrains the synchrotron and self-Compton emission from the gamma-ray emitting electrons, which requires the magnetic field to be relatively weak, with the magnetic energy density < a few times 10^-3 of that in the electrons. The actual value of the magnetic field strength can be inferred from a future simultaneous measurement of the IR and gamma-ray fluxes.

preprint2012arXiv

The Radiative Efficiency of Hot Accretion Flows

Two significant progresses have been made in the past years on our understanding of hot accretion flows. One is that only a small fraction of accretion flow available at the outer boundary can finally falls onto the black hole while most of them is lost in outflow. Another one is that electrons may directly receive a large fraction of the viscously dissipated energy in the accretion flow, i.e, $δ\sim 0.1-0.5$. The radiative efficiency of hot accretion flow when these two progresses are taken into account has not been systematically studied and is the subject of the present paper. We consider two regimes of hot accretion model. One is the advection dominated accretion flows (ADAFs) which lie on low accretion rate regime, $\la 10α^2\ledd/c^2$; another being the luminous hot accretion flows (LHAFs) which lie above this accretion rate. For the latter, we assume that the accretion flow will has a two-phase structure above a certain accretion rate, and a simplification is adopted in our calculation of the dynamics. Our results indicate that the radiative efficiency of hot accretion flow increases with the accretion rate and is highly enhanced by the direct viscous heating to electrons compared to the previous case of $δ\ll 1$. When the accretion rate is high, the radiative efficiency of hot accretion flow is comparable to that of the standard thin disk. Fitting formulae of radiative efficiency as a function of accretion rate for various $δ$ values are presented.

preprint2012arXiv

Transverse Polarization of the Nucleon in Parton Picture

The proton's transverse polarization structure is examined in terms of the Lorentz-covariant Pauli-Lubanski vector in QCD. We find that there are contributions from leading, subleading, and next-to-subleading partonic contributions in the light-front system of coordinates. The subleading and next-to-subleading contributions are related to the leading one through Lorentz symmetry. And the leading contribution obeys a simple partonic angular momentum sum rule that gives a clear physical interpretation to a relation known previously.

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

Gluon Distribution Functions and Higgs Boson Production at Moderate Transverse Momentum

We investigate the gluon distribution functions and their contributions to the Higgs boson production in pp collisions in the transverse momentum dependent factorization formalism. In addition to the usual azimuthal symmetric transverse momentum dependent gluon distribution, we find that the azimuthal correlated gluon distribution also contributes to the Higgs boson production. This explains recent findings on the additional contribution in the transverse momentum resummation for the Higgs boson production as compared to that for electroweak boson production processes. We further examine the small-x naive k_t-factorization in the dilute region and find that the azimuthal correlated gluon distribution contribution is consistently taken into account, and the result agrees with the transverse momentum dependent factorization formalism. We comment on the possible breakdown of the naive k_t-factorization in the dense medium region, due to the unique behaviors for the gluon distributions.

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

Magnetothermal and magnetorotational instabilities in hot accretion flows

In a hot, dilute, magnetized accretion flow, the electron mean-free path can be much greater than the Larmor radius, thus thermal conduction is anisotropic and along magnetic field lines. In this case, if the temperature decreases outward, the flow may be subject to a buoyancy instability (the magnetothermal instability, or MTI). The MTI amplifies the magnetic field, and aligns field lines with the radial direction. If the accretion flow is differentially rotating, the magnetorotational instability (MRI) may also be present. Using two-dimensional, time-dependent magnetohydrodynamic simulations, we investigate the interaction between these two instabilities. We use global simulations that span over two orders of magnitude in radius, centered on the region around the Bondi radius where the infall time of gas is longer than the growth time of both the MTI and MRI. Significant amplification of the magnetic field is produced by both instabilities, although we find that the MTI primarily amplifies the radial component, and the MRI primarily the toroidal component, of the field, respectively. Most importantly, we find that if the MTI can amplify the magnetic energy by a factor $F_t$, and the MRI by a factor $F_r$, then when the MTI and MRI are both present, the magnetic energy can be amplified by a factor of $F_t \cdot F_r$. We therefore conclude that amplification of the magnetic energy by the MTI and MRI operates independently. We also find that the MTI contributes to the transport of angular momentum, because radial motions induced by the MTI increase the Maxwell (by amplifying the magnetic field) and Reynolds stresses. Finally, we find that thermal conduction decreases the slope of the radial temperature profile. The increased temperature near the Bondi radius decreases the mass accretion rate.

preprint2011arXiv

Observational evidence for AGNs feedback at parsec scale

In a hot accretion flow, the radiation from the innermost region of the flow propagates outward and heats the electrons at large radii via Compton scattering. It has been shown in previous works that if the radiation is strong enough, $L\ga 2%L_{\rm Edd}$, the electrons at the Bondi radius ($r_B\sim 10^5 r_s$) will be heated to be above the virial temperature thus the accretion will be stopped. The accretion will recover after the gas cools down. This results in the oscillation of the black hole activity. In this paper we show that this mechanism is the origin of the intermittent activity of some compact young radio sources. Such intermittency is required to explain the population of these sources. We calculate the timescales of the black hole oscillation and find that the durations of active and inactive phases are $3\times 10^4 (0.1/α)(M/10^8\msun) (L/2%L_{\rm Edd})^{-1/2} {\rm yr}$ and $10^5(α/0.1)(M/10^8\msun) {\rm yr}$, respectively, consistent with those required to explain observations. Such kind of feedback occurring at parsec scale should be common in low-luminosity AGNs and should be considered when we consider their matter and energy output.

preprint2011arXiv

One-loop Factorization for Inclusive Hadron Production in $pA$ Collisions in the Saturation Formalism

We demonstrate the QCD factorization for inclusive hadron production in $pA$ collisions in the saturation formalism at one-loop order, with explicit calculation of both real and virtual gluon radiation diagrams. The collinear divergences associated with the incoming parton distribution of the nucleon and the outgoing fragmentation function of the final state hadron, as well as the rapidity divergence with small-$x$ dipole gluon distribution of the nucleus are factorized into the splittings of the associated parton distribution and fragmentation functions and the energy evolution of the dipole gluon distribution function. The hard coefficient function is evaluated at one-loop order, and contains no divergence.

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

Revisiting the Thermal Stability of Radiation-dominated Thin Disks

The standard thin disk model predicts that when the accretion rate is over a small fraction of the Eddington rate, which corresponds to $L \ga 0.06 L_{Edd}$, the inner region of the disk is radiation-pressure-dominated and thermally unstable. However, observations of the high/soft state of black hole X-ray binaries with luminosity well within this regime ($0.01L_{Edd} \la L \la 0.5L_{Edd}$) indicate that the disk has very little variability, i.e., quite stable. Recent radiation magnetohydrodynamic simulations of a vertically stratified shearing box have confirmed the absence of the thermal instability. In this paper, we revisit the thermal stability by linear analysis, taking into account the role of magnetic field in the accretion flow. By assuming that the field responses negatively to a positive temperature perturbation, we find that the threshold of accretion rate above which the disk becomes thermally unstable increases significantly compared with the case of not considering the role of magnetic field. This accounts for the stability of the observed sources with high luminosities. Our model also presents a possible explanation as to why only GRS 1915+105 seems to show thermally unstable behavior. This peculiar source holds the highest accretion rate (or luminosity) among the known high state sources, which is well above the accretion rate threshold of the instability.

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

Advection-dominated accretion, jets and the spectral energy distribution of LINERs

Low-luminosity active galactic nuclei (LLAGNs) represent the bulk of the AGN population in the present-day universe and they trace low-level accreting supermassive black holes. The observational properties of LLAGNs suggest that their central engines are intrinsically different from those of more luminous AGNs. It has been suggested that accretion in LLAGNs occurs via an advection-dominated accretion flow (ADAF) associated with strong jets. In order to probe the accretion physics in LLAGNs as a class, we model the multiwavelength spectral energy distributions (SEDs) of 24 LINERs (taken from a recent compilation by Eracleous et al.) with a coupled accretion-jet model. The accretion flow is modeled as an inner ADAF outside of which there is a truncated standard thin disk. These SEDs include radio, near-IR to near-UV HST data, and Chandra X-ray data. We find that the radio emission is severely underpredicted by ADAF models but can be explained by the relativistic jet. The origin of the X-ray radiation in most sources can be explained by three distinct scenarios: the X-rays can be dominated by emission from the ADAF, the jet, or from both components contributing at similar levels. From the model fits, we estimate important parameters of the central engine of LINERs, such as the mass accretion rate -- relevant for studies of the feeding of AGNs -- and the mass-loss rate in the jet and the jet power -- relevant for studies of the kinetic feedback from jets.

preprint2010arXiv

Dihadron Azimuthal Correlation from Collins Effect in Unpolarized Hadron Collisions

We study the dihadron azimuthal correlation produced nearly back-to-back in unpolarized hadron collisions, arising from the product of two Collins fragmentation functions. Using the latest Collins fragmentation functions extracted from the global analysis of available experimental data, we make predictions for the azimuthal correlation of two-pion production in $pp$ collisions at RHIC energies. We find that the correlation is sizable in the mid-rapidity region for moderate jet transverse momentum.

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

Modeling the hard states of three black hole candidates

Simultaneous multiwavelength observations were recently performed for three black hole candidates --- SWIFT J1753.5-0127, GRO J1655-40 and XTE J1720-318. In this paper we test the accretion-jet model originally proposed to XTE J1118+480 by investigating the hard state of these three sources using this model. The accretion flow in the model is composed of an inner hot accretion flow and an outer truncated thin disk. We find that the model satisfactorily explains the spectrum ranging from radio to X-rays, with the radio and X-ray spectra dominated by the synchrotron and thermal Comptonization emissions in the jet and the hot accretion flow respectively, while the infrared and optical being the sum of the emissions from the jet, hot accretion flow, and the truncated thin disk. Similar to the case of XTE J1118+480, the model can also explain, although only qualitatively in some cases, the observed timing features including QPO, and positive and negative time lags between the optical and X-ray emissions detected in SWIFT J1753.5-0127. The origin of the ejection events detected in XTE J1720-318 is also briefly discussed.

preprint2010arXiv

Non-Universality of Transverse Momentum Dependent Parton Distributions at Small-x

We study the universality issue of the transverse momentum dependent parton distributions at small-x, by comparing the initial/final state interaction effects in dijet-correlation in pA collisions with those in deep inelastic lepton nucleus scattering. We demonstrate the non-universality by performing an explicit calculation in a particular model where the multiple gauge boson exchange contributions are summed up to all orders. We comment on the implications of our results on the theoretical interpretation of di-hadron correlation in dA collisions in terms of the saturation phenomena in deep inelastic lepton nucleus scattering.

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

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

On the Origin of Ultraviolet Emission and the Accretion Model of Low-luminosity AGNs

Low-luminosity active galactic nuclei (LLAGNs) are generally believed to be powered by an inner radiatively inefficient, advection-dominated accretion flow (ADAF), an outer truncated thin disk, and a jet. Maoz (2007) recently challenged this picture based on the observation that the strength of ultraviolet emission relative to the X-ray and radio bands does not depart from empirical trends defined by more luminous sources. He advocates that AGNs across all luminosities have essentially the same accretion and radiative processes, which in luminous sources are described by a standard optically thick, geometrically thin disk. We calculate ADAF models and demonstrate that they can successfully fit the observed spectral energy distributions of the LLAGNs in Maoz's sample. Our model naturally accommodates the radio and X-ray emission, and the ultraviolet flux is well explained by a combination of the first-order Compton scattering in the ADAF, synchrotron emission in the jet, and black body emission in the truncated thin disk. It is premature to dismiss the ADAF model for LLAGNs. The UV data can be fit equally well using a standard thin disk, but an additional corona and jet would be required to account for the X-ray and radio emission. We argue that there are strong theoretical reasons to prefer the ADAF model over the thin disk scenario. We discuss testable predictions that can potentially discriminate between the two accretion models.

preprint2010arXiv

Phases of Augmented Hadronic Light-Front Wave Functions

It is an important question whether the final/initial state gluonic interactions which lead to naive-time-reversal-odd single-spin asymmetries and diffraction at leading twist can be associated in a definite way with the light-front wave function hadronic eigensolutions of QCD. We use light-front time-ordered perturbation theory to obtain augmented light-front wave functions which contain an imaginary phase which depends on the choice of advanced or retarded boundary condition for the gauge potential in light-cone gauge. We apply this formalism to the wave functions of the valence Fock states of nucleons and pions, and show how this illuminates the factorization properties of naive-time-reversal-odd transverse momentum dependent observables which arise from rescattering. In particular, one calculates the identical leading-twist Sivers function from the overlap of augmented light-front wavefunctions that one obtains from explicit calculations of the single-spin asymmetry in semi-inclusive deep inelastic lepton-polarized nucleon scattering where the required phases come from the final-state rescattering of the struck quark with the nucleon spectators.

preprint2010arXiv

Quark splitting in non-trivial θ-vacuum

Quark splitting in non-trivial $θ$-vacuum with a given helicity is investigated in pQCD with a modified quark propagator. We found that the quark splitting functions were modified by the presence of a topologically non-trivial QCD background field, though there is no explicit helicity flip associated with the radiative processes. The interaction with the topological non-trivial field leads to the degeneracy of the quark splitting functions for left- and right-handed quarks. Such degeneracy can lead to imbalance of left- and right-handed quarks in quark jet showers. We also discuss phenomenological consequences of such imbalance if there exists non-trivial topological gluon field configuration in heavy-ion collisions.

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.

preprint2010arXiv

Twist-three Fragmentation Function Contribution to the Single Spin Asymmetry in pp Collisions

We study the twist-three fragmentation function contribution to the single transverse spin asymmetries in inclusive hadron production in pp collisions, pp->h+X. In particular, we evaluate the so-called derivative contribution which dominates the spin asymmetry in the forward direction of the polarized proton. With certain parametrizations for the twist-three fragmentation function, we estimate its contribution to the asymmetry of pi0 production at RHIC energy. We find that the contribution is sizable and might be responsible for the big difference between the asymmetries in eta and pi0 productions observed by the STAR collaboration at RHIC.

preprint2009arXiv

Electron heating and acceleration by magnetic reconnection in hot accretion flows

Both analytical and numerical works show that magnetic reconnection must occur in hot accretion flows. This process will effectively heat and accelerate electrons. In this paper we use the numerical hybrid simulation of magnetic reconnection plus test-electron method to investigate the electron acceleration and heating due to magnetic reconnection in hot accretion flows. We consider fiducial values of density, temperature, and magnetic parameter $β_e$ (defined as the ratio of the electron pressure to the magnetic pressure) of the accretion flow as $n_{0} \sim 10^{6} {\rm cm^{-3}}$, $T_{e}^0\sim 2\times 10^9 {\rm K}$, and $β_e=1$. We find that electrons are heated to a higher temperature $T_{e}=5\times 10^9$K, and a fraction $η\sim 8%$ of electrons are accelerated into a broken power-law distribution, $dN(γ)\propto γ^{-p}$, with $p\approx 1.5$ and 4 below and above $\sim 1$ MeV, respectively. We also investigate the effect of varying $β$ and $n_0$. We find that when $β_e$ is smaller or $n_0$ is larger, i.e, the magnetic field is stronger, $T_e$, $η$, and $p$ all become larger.

preprint2009arXiv

Monte Carlo simulations of global Compton cooling in inner regions of hot accretion flows

Hot accretion flows such as advection-dominated accretion flows are generally optically thin in the radial direction. Thus photons generated at some radii can cool or heat electrons at other radii via Compton scattering. Such global Compton scattering has previously been shown to be important for the dynamics of accretion flows. Here, we extend previous treatments of this problem by using accurate global general relativistic Monte Carlo simulations. We focus on an inner region of the accretion flow (R < 600R_g), for which we obtain a global self-consistent solution. As compared to the initial, not self-consistent solution, the final solution has both the cooling rate and the electron temperature significantly reduced at radii >=10 gravitational radii. On the other hand, the radiation spectrum of the self-consistent solution has the shape similar to that of the initial iteration, except for the high-energy cut-off being at an energy lower by a factor of ~2 and the bolometric luminosity decreased by a factor of ~2. We also compare the global Compton scattering model with local models in spherical and slab geometry. We find that the slab model approximates the global model significantly better than the spherical one. Still, neither local model gives a good approximation to the radial profile of the cooling rate, and the differences can be up to two orders of magnitude. The local slab model underestimates the cooling rate at outer regions whereas it overestimates that rate at inner regions. We compare our modelling results to observed hard-state spectra of black-hole binaries and find an overall good agreement provided any disc outflow is weak. We find that general-relativistic effects in flows which dynamics is modified by global Comptonization is crucial in approaching this agreement.

preprint2009arXiv

Reprocessed emission from warped accretion discs induced by the Bardeen-Petterson effect

The broad Balmer emission-line profiles resulting from the reprocessing of UV/X-ray radiation from a warped accretion disc induced by the Bardeen-Petterson effect are studied. We adopt a thin warped disc geometry and a central ring-like illuminating source in our model. We compute the steady-state shape of the warped disc numerically, and then use it in the calculation of the line profile. We find that, from the outer radius to the inner radius of the disc, the warp is twisted by an angle of $\simπ$ before being flattened efficiently into the equatorial plane. The profiles obtained depend weakly on the illuminating source radius in the range from $3r_{g}$ to $10r_g$, but depend strongly on this radius when it approaches the marginally stable orbit of an extreme Kerr black hole. Double- or triplet-peaked line profiles are present in most cases when the illuminating source radius is low. The triplet-peaked line profiles observed from the Sloan Digital Sky Survey may be a {"}signature" of a warped disc.

preprint2009arXiv

Steady state solution of warped accretion discs

We consider a thin accretion disc warped due to the Bardeen-Petterson effect, presenting both analytical and numerical solutions for the situation that the two viscosity coefficients vary with radius as power law, with the two power law indices not necessarily equal. The analytical solutions are compared with numerical ones, showing that our new analytical solution is more accurate than previous one, which overestimates the inclination changing in the outer disc. Our new analytical solution is appropriate for moderately warped discs, while for extremely misaligned disc, only numerical solution is appropriate.

preprint2009arXiv

Transverse momentum dependent quark distributions and polarized Drell-Yan processes

We study the spin-dependent quark distributions at large transverse momentum. We derive their transverse momentum behaviors in the collinear factorization approach in this region. We further calculate the angular distribution of the Drell-Yan lepton pair production with polarized beams and present the results in terms of the collinear twist-three quark-gluon correlation functions. In the intermediate transverse momentum region, we find that the two pproaches: the collinear factorization and the transverse momentum dependent factorization approaches are consistent in the description of the lepton pair angular distributions.

preprint2008arXiv

Origin of the X-ray Emission in the Nuclei of FR Is

We investigate the X-ray origin in FRIs using the multi-waveband high resolution data of eight FR I sources, which have very low Eddington ratios. We fit their multi-waveband spectrum using a coupled accretion-jet model. We find that X-ray emission in the source with the highest L_X (~1.8*10^-4 L_Edd) is from the advection-dominated accretion flow (ADAF). Four sources with moderate L_X(~several*10^-6 L_Edd) are complicated. The X-ray emission of one FR I is from the jet, and the other three is from the sum of the jet and ADAF. The X-ray emission in the three least luminous sources (L_X<1.0*10^-6L_Edd) is dominated by the jet. These results roughly support the predictions of Yuan and Cui(2005) where they predict that when the X-ray luminosity of the system is below a critical value, the X-radiation will not be dominated by the emission from the ADAF any longer, but by the jet. We also find that the accretion rates in four sources must be higher than the Bondi rates, which implies that other fuel supply (e.g., stellar winds) inside the Bondi radius should be important.

preprint2008arXiv

The Influences of Outflow on the Dynamics of Inflow

Both numerical simulations and observations indicate that in an advection-dominated accretion flow most of the accretion material supplied at the outer boundary will not reach the inner boundary. Rather, they are lost via outflow. Previously, the influence of outflow on the dynamics of inflow is taken into account only by adopting a radius-dependent mass accretion rate $\dot{M}=\dot{M}_0 (r/r_{\rm out})^s$ with $s>0$. In this paper, based on a 1.5 dimensional description to the accretion flow, we investigate this problem in more detail by considering the interchange of mass, radial and azimuthal momentum, and the energy between the outflow and inflow. The physical quantities of the outflow is parameterized based on our current understandings to the properties of outflow mainly from numerical simulations of accretion flows. Our results indicate that under reasonable assumptions to the properties of outflow, the main influence of outflow has been properly included by adopting $\dot{M}=\dot{M}_0 (r/r_{\rm out})^s$.

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.

preprint2003arXiv

Sivers Function in the MIT Bag Model

The Sivers function, an asymmetric transverse-momentum distribution of the quarks in a transversely polarized nucleon, is calculated in the MIT bag model. The bag quark wave functions contain both $S$-wave and $P$-wave components, and their interference leads to nonvanishing Sivers function in the presence of the final state interactions. We approximate these interactions through one-gluon exchange. An estimate of another transverse momentum dependent distribution $h_1^\perp$ is also performed in the same model.

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.

preprint2002arXiv

Scale Dependence of Hadronic Wave Functions and Parton Densities

We study how the components of hadronic wave functions in light-cone quantization depend on the ultraviolet cut-off by relating them in a systematic way to the matrix elements of a class of quark-gluon operators between the QCD vacuum and the hadrons. From this, we derive an infinite set of scale-evolution equations for the individual contributions to parton distributions from the Fock expansion. When summed over all the contributions, we recover the well-known DGLAP equation.

preprint2001arXiv

Enhanced Minijet Production in $A-A$ Collisions from Gluons with Large Transverse Momenta

We find supersaturation for the intrinsic gluon distribution of nuclei, i.e. the low $x$ unintegrated nuclear gluon distribution peaks at intermediate transverse momenta $k_t=Q_s$ and vanishes at zero $k_t$. Taking into account the intrinsic transverse momenta of the gluons and the saturation scale $p_s$ of the produced gluons, we calculate the minijet cross section arising from gluon gluon scattering for RHIC energies. For central collisions at $\sqrt{s}=200GeV$ the saturation scale $p_s \approx 1.4$ GeV is larger due to intrinsic $k_t$ effects and increases with energy. Our theoretical results on charged particle multiplicity agree very well with the recent experimental data from RHIC.

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

A New Parameter In Accretion Disk Model

Taking optically thin accretion flows as an example, we investigate the dynamics and the emergent spectra of accretion flows with different outer boundary conditions (OBCs) and find that OBC plays an important role in accretion disk model. This is because the accretion equations describing the behavior of accretion flows are a set of {\em differential} equations, therefore, accretion is intrinsically an initial-value problem. We argue that optically thick accretion flow should also show OBC-dependent behavior. The result means that we should seriously consider the initial physical state of the accretion flow such as its angular momentum and its temperature. An application example to Sgr A$^*$ is presented.

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

Luminous hot accretion disks

We find a new two-temperature hot branch of equilibrium solutions for stationary accretion disks around black holes. In units of Eddington accretion rate defined as $10L_{\rm Edd}/c^2$, the accretion rates to which these solutions correspond are within the range $\dot{m}_1 \la \dot{m} \la 1$, here $\dot{m}_1$ is the critical rate of advecti$. In these solutions, the energy loss rate of the ions by Coulomb energy transfer between the ions and electrons is larger than the viscously heating rate and it is the advective heating together with the viscous dissipation that balances the Coulomb cooling of ions. When $\dot{m}_1 \la \dot{m} \la \dot{m}_2$, where $\dot{m}_2 \sim 5 \dot{m}_1 < 1$, the accretion flow remains hot throughout the disk. When $\dot{m}_2 \la \dot{m}\la 1$, Coulomb interaction will cool the inner region of the disk within a certain radius ($r_{\rm tr} \sim$ several -- tens of Schwarzschild radii or larger depending on the accretion rate and the outer boundary condition) and the disk will collapse onto the equatorial plane and form an optically thick cold annulus. Compared to ADAF, these hot solutions are much more luminous because of the high accretion rate and efficiency, therefore, we name them luminous hot accretion disks.

preprint2000arXiv

Pseudogap effects on the c-axis charge dynamics in copper oxide materials

The c-axis charge dynamics of copper oxide materials in the underdoped and optimally doped regimes has been studied by considering the incoherent interlayer hopping. It is shown that the c-axis charge dynamics for the chain copper oxide materials is mainly governed by the scattering from the in-plane fluctuation, and the c-axis charge dynamics for the no-chain copper oxide materials is dominated by the scattering from the in-plane fluctuation incorporating with the interlayer disorder, which would be suppressed when the holon pseudogap opens at low temperatures and lower doping levels, leading to the crossovers to the semiconducting-like range in the c-axis resistivity and the temperature linear to the nonlinear range in the in-plane resistivity.

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.