Source author record

Y. Huang

Y. Huang 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

218works
26topics
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

218 published item(s)

preprint2026arXiv

An Enigmatic PeVatron in an Area around HII Region G35.6$-$0.5

Identifying Galactic PeVatrons (PeV particle accelerators) from the ultra-high-energy (UHE, >100 TeV) $γ$-ray sources plays a crucial role in revealing the origin of Galactic cosmic rays. The UHE source 1LHAASO J1857+0203u is suggested to be associated with HESS J1858+020, which may be attributed to the possible PeVatron candidate supernova remnant (SNR) G35.6$-$0.4 or HII region G35.6$-$0.5. We perform detailed analysis on the very-high-energy and UHE $γ$-ray emissions towards this region with data from the Large High Altitude Air Shower Observatory (LHAASO). 1LHAASO J1857+0203u is detected with a significance of 11.6$σ$ above 100 TeV, indicating the presence of a PeVatron. It has an extension of $\sim 0.18^\circ$ with a power-law (PL) spectral index of $\sim$2.5 in 1-25 TeV and a point-like emission with a PL spectral index of $\sim$3.2 above 25 TeV. Using the archival CO and HI data, we identify some molecular and atomic clouds that may be associated with the TeV $γ$-ray emissions. Our modelling indicates that the TeV $γ$-ray emissions are unlikely to arise from the clouds illuminated by the protons that escaped from SNR G35.6$-$0.4. In the scenario that HII region G35.6$-$0.5 could accelerate particles to the UHE band, the observed GeV-TeV $γ$-ray emission could be well explained by a hadronic model with a PL spectral index of $\sim$2.0 and cutoff energy of $\sim$450 TeV. However, an evolved pulsar wind nebula origin cannot be ruled out.

preprint2026arXiv

An Ultrahigh-energy $γ$-ray Bubble Powered by a Super PeVatron

We report the detection of a $γ$-ray bubble spanning at least 100$\rm deg^2$ in ultra high energy (UHE) up to a few PeV in the direction of the star-forming region Cygnus X, implying the presence Super PeVatron(s) accelerating protons to at least 10 PeV. A log-parabola form with the photon index $Γ(E) = (2.71 \pm 0.02) + (0.11 \pm 0.02) \times \log_{10} (E/10 \ {\rm TeV})$ is found fitting the gamma-ray energy spectrum of the bubble well. UHE sources, `hot spots' correlated with very massive molecular clouds, and a quasi-spherical amorphous $γ$-ray emitter with a sharp central brightening are observed in the bubble. In the core of $\sim 0.5^{\circ}$, spatially associating with a region containing massive OB association (Cygnus OB2) and a microquasar (Cygnus X-3), as well as previously reported multi-TeV sources, an enhanced concentration of UHE $γ$-rays are observed with 2 photons at energies above 1 PeV. The general feature of the bubble, the morphology and the energy spectrum, are reasonably reproduced by the assumption of a particle accelerator in the core, continuously injecting protons into the ambient medium.

preprint2026arXiv

Constraining the Cosmic-ray Energy Based on Observations of Nearby Galaxy Clusters by LHAASO

Galaxy clusters act as reservoirs of high-energy cosmic rays (CRs). As CRs propagate through the intracluster medium, they generate diffuse $γ$-rays detectable by arrays such as LHAASO. These $γ$-rays result from proton-proton ($pp$) collisions of very high-energy cosmic rays (VHECRs) or inverse Compton (IC) scattering of positron-electron pairs created by $pγ$ interactions of ultra-high-energy cosmic rays (UHECRs). We analyzed diffuse $γ$-ray emission from the Coma, Perseus, and Virgo clusters using LHAASO data. Diffuse emission was modeled as a disk of radius $R_{500}$ for each cluster while accounting for point sources. No significant diffuse emission was detected, yielding 95\% confidence level (C.L.) upper limits on the $γ$-ray flux: for WCDA (1-25~TeV) and KM2A ($>25$~TeV), less than $(49.4, 13.7, 54.0)$ and $(1.34, 1.14, 0.40) \times 10^{-14}$~ph~cm$^{-2}$~s$^{-1}$ for Coma, Perseus, and Virgo, respectively. The $γ$-ray upper limits can be used to derive model-independent constraints on the integral energy of CRp above 10~TeV (corresponding to the LHAASO observational range $>1$~TeV under the $pp$ scenario) to be less than $(1.96, 0.59, 0.08) \times 10^{61}$~erg. The absence of detectable annuli/ring-like structures, indicative of cluster accretion or merging shocks, imposes further constraints on models in which the UHECRs are accelerated in the merging shocks of galaxy clusters.

preprint2026arXiv

Constraints on heavy decaying dark matter from 570 days of LHAASO observations

The Kilometer Square Array~(KM2A) of the Large High Altitude Air Shower Observatory (LHAASO) aims at surveying the northern gamma-ray sky at energies above 10 TeV with unprecedented sensitivity. Gamma-ray observations have long been one of the most powerful tools for dark matter searches, as e.g., high-energy gamma-rays could be produced by the decays of heavy dark matter particles. In this letter, we present the first dark matter analysis with LHAASO-KM2A, using the first 340~days of data from 1/2-KM2A and 230~days of data from 3/4-KM2A. Several regions of interest are used to search for a signal and account for the residual cosmic-ray background after gamma/hadron separation. We find no excess of dark matter signals, and thus place some of the strongest gamma-ray constraints on the lifetime of heavy dark matter particles with mass between 10^5 and 10^9~GeV. Our results with LHAASO are robust, and have important implications for dark matter interpretations of the diffuse astrophysical high-energy neutrino emission.

preprint2026arXiv

Deep view of Composite SNR CTA1 with LHAASO in $γ$-rays up to 300 TeV

The ultra-high-energy (UHE) gamma-ray source 1LHAASO J0007+7303u is positionally associated with the composite SNR CTA1 that is located at high Galactic Latitude $b\approx 10.5^\circ$. This provides a rare opportunity to spatially resolve the component of the pulsar wind nebula (PWN) and supernova remnant (SNR) at UHE. This paper conducted a dedicated data analysis of 1LHAASO J0007+7303u using the data collected from December 2019 to July 2023. This source is well detected with significances of 21$σ$ and 17$σ$ at 8$-$100 TeV and $>$100 TeV, respectively. The corresponding extensions are determined to be 0.23$^{\circ}\pm$0.03$^{\circ}$ and 0.17$^{\circ}\pm$0.03$^{\circ}$. The emission is proposed to originate from the relativistic electrons and positrons accelerated within the PWN of PSR J0007+7303. The energy spectrum is well described by a power-law with an exponential cutoff function $dN/dE = (42.4\pm4.1)(\frac{E}{20\rm\ TeV})^{-2.31\pm0.11}\exp(-\frac{E}{110\pm25\rm\ TeV})$ $\rm\ TeV^{-1}\ cm^{-2}\ s^{-1}$in the energy range from 8 TeV to 300 TeV, implying a steady-state parent electron spectrum $dN_e/dE_e\propto (\frac{E_e}{100\rm\ TeV})^{-3.13\pm0.16}\exp[(\frac{-E_e}{373\pm70\rm\ TeV})^2]$ at energies above $\approx 50 \rm\ TeV$. The cutoff energy of the electron spectrum is roughly equal to the expected current maximum energy of particles accelerated at the PWN terminal shock. Combining the X-ray and gamma-ray emission, the current space-averaged magnetic field can be limited to $\approx 4.5\rm\ μG$. To satisfy the multi-wavelength spectrum and the $γ$-ray extensions, the transport of relativistic particles within the PWN is likely dominated by the advection process under the free-expansion phase assumption.

preprint2026arXiv

Energy calibration of LHAASO-KM2A using the cosmic ray Moon shadow

We present a precise measurement of the westward, rigidity-dependent shift of the Moon's shadow using three and a half years of cosmic-ray data collected by the Kilometer Square Array (KM2A) of the Large High Altitude Air Shower Observatory (LHAASO). These measurements enable us to calibrate the detector energy response in the range 20-260 TeV, with results showing excellent agreement with the response derived from Monte Carlo (MC) simulations of the KM2A detector. We also measure a best-fit parameter $ε= 0.015 \pm 0.08$, corresponding to a 95% confidence interval of [-14%, +17%] for the energy-scale estimation. This result establishes the exceptional accuracy of the KM2A-MC in simulating the detector's response within this energy range.

preprint2026arXiv

Energy-Dependent Shifts of Medium-Scale Anisotropies in Very-High-Energy Cosmic Rays Observed by LHAASO-KM2A

Small deviations from isotropy in the arrival directions of Galactic cosmic rays serve as a unique probe of the local magnetic environment. In this Letter, we report observations of medium-scale anisotropies (MSA) at energies above 10 TeV using the LHAASO-KM2A array. Our analysis identifies four regions of excess and four regions of deficit, each spanning angular scales of approximately ten degrees. Crucially, we detect significant energy-dependent shifts in the centroids of two excess regions: Region B and the newly identified Region $\mathrm{\widetilde{D}}$. We also characterize the energy evolution of the fractional relative intensity across both excess and deficit regions. These findings imply that the observed anisotropies are shaped by the specific realization of the local turbulent magnetic field within the cosmic ray scattering length. Such energy-dependent behaviors impose strict constraints on local turbulence models and cosmic ray propagation theories.

preprint2026arXiv

Evidence for particle acceleration approaching PeV energies in the W51 complex

The $γ$-ray emission from the W51 complex is widely acknowledged to be attributed to the interaction between the cosmic rays (CRs) accelerated by the shock of supernova remnant (SNR) W51C and the dense molecular clouds in the adjacent star-forming region, W51B. However, the maximum acceleration capability of W51C for CRs remains elusive. Based on observations conducted with the Large High Altitude Air Shower Observatory (LHAASO), we report a significant detection of $γ$ rays emanating from the W51 complex, with energies from 2 TeV to 200 TeV. The LHAASO measurements, for the first time, extend the $γ$-ray emission from the W51 complex beyond 100 TeV and reveal a significant spectrum bending at tens of TeV. By combining the ``$π^0$-decay bump" featured data from Fermi-LAT, the broadband $γ$-ray spectrum of the W51 region can be well-characterized by a simple pp-collision model. The observed spectral bending feature suggests an exponential cutoff at $\sim400$~TeV or a power-law break at $\sim200$~TeV in the CR proton spectrum, most likely providing the first evidence of SNRs serving as CR accelerators approaching the PeV regime. Additionally, two young star clusters within W51B could also be theoretically viable to produce the most energetic $γ$ rays observed by LHAASO. Our findings strongly support the presence of extreme CR accelerators within the W51 complex and provide new insights into the origin of Galactic CRs.

preprint2026arXiv

Extended Very-High-Energy Gamma-Ray Emission Surrounding PSR J0622 + 3749 Observed by LHAASO-KM2A

We report the discovery of an extended very-high-energy (VHE) gamma-ray source around the location of the middle-aged (207.8 kyr) pulsar PSR J0622+3749 with the Large High Altitude Air Shower Observatory (LHAASO). The source is detected with a significance of $8.2σ$ for $E>25$~TeV assuming a Gaussian template. The best-fit location is (R.A., Dec.)$=(95^{\circ}\!.47\pm0^{\circ}\!.11,\,37^{\circ}\!.92 \pm0^{\circ}\!.09)$, and the extension is $0^{\circ}\!.40\pm0^{\circ}\!.07$. The energy spectrum can be described by a power-law spectrum with an index of ${-2.92 \pm 0.17_{\rm stat} \pm 0.02_{\rm sys} }$. No clear extended multi-wavelength counterpart of the LHAASO source has been found from the radio to sub-TeV bands. The LHAASO observations are consistent with the scenario that VHE electrons escaped from the pulsar, diffused in the interstellar medium, and scattered the interstellar radiation field. If interpreted as the pulsar halo scenario, the diffusion coefficient, inferred for electrons with median energies of $\sim160$~TeV, is consistent with those obtained from the extended halos around Geminga and Monogem and much smaller than that derived from cosmic ray secondaries. The LHAASO discovery of this source thus likely enriches the class of so-called pulsar halos and confirms that high-energy particles generally diffuse very slowly in the disturbed medium around pulsars.

preprint2026arXiv

LHAASO Detection of Ultra-High-Energy Gamma-Ray Emission toward the Giant Molecular Clouds

The $γ$-ray from Giant molecular clouds (GMCs) is regarded as the most ideal tool to perform in-situ measurement of cosmic ray (CR) density and spectra in our Galaxy. We report the first detection of $γ$-ray emissions in the very-high-energy (VHE) domain from the five nearby GMCs with a stacking analysis based on a 4.5-year $γ$-ray observation with the Large High Altitude Air Shower Observatory (LHAASO) experiment. The spectral energy distributions derived from the GMCs are consistent with the expected $γ$-ray flux produced via CR interacting with the ISM in the energy interval 1 - 100 $~\rm$ TeV. In addition, we investigate the presence of the CR spectral `knee' by introducing a spectral break in the $γ$-ray data. While no significant evidence for the CR knee is found, the current KM2A measurements from GMCs strongly favor a proton CR knee located above 0.9$~\rm$ PeV, which is consistent with the latest measurement of the CR spectrum by ground-based experiments.

preprint2026arXiv

Measurement of attenuation length of the muon content in extensive air showers from 0.3 to 30 PeV with LHAASO

The attenuation length of the muon content in extensive air showers provides important information regarding the generation and development of air showers. This information can be used not only to improve the description of such showers but also to test fundamental models of hadronic interactions. Using data from the LHAASO-KM2A experiment, the development of the muon content in high-energy air showers was studied. The attenuation length of muon content in the air showers was measured from experimental data in the energy range from 0.3 to 30 PeV using the constant intensity cut method. By comparing the attenuation length of the muon content with predictions from high-energy hadronic interaction models (QGSJET-II-04, SIBYLL 2.3d, and EPOS-LHC), it is evident that LHAASO results are significantly shorter than those predicted by the first two models (QGSJET-II-04 and SIBYLL 2.3d) but relatively close to those predicted by the third model (EPOS-LHC). Thus, the LHAASO data favor the EPOS-LHC model over the other two models. The three interaction models confirmed an increasing trend in the attenuation length as the cosmic-ray energy increases.

preprint2026arXiv

Measurement of Very-high-energy Diffuse Gamma-ray Emissions from the Galactic Plane with LHAASO-WCDA

The diffuse Galactic gamma-ray emission is a very important tool used to study the propagation and interaction of cosmic rays in the Milky Way. In this work, we report the measurements of the diffuse emission from the Galactic plane, covering Galactic longitudes from $15^{\circ}$ to $235^{\circ}$ and latitudes from $-5^{\circ}$ to $+5^{\circ}$, in an energy range of 1 TeV to 25 TeV, with the Water Cherenkov Detector Array (WCDA) of the Large High Altitude Air Shower Observatory (LHAASO). After masking the sky regions of known sources, the diffuse emission is detected with $24.6σ$ and $9.1σ$ significance in the inner Galactic plane and outer Galactic plane, respectively. The WCDA spectra in both regions can be well described by a power-law function, with spectral indices of $-2.67\pm0.05_{\rm stat}$ in the inner region and $-2.83\pm0.19_{\rm stat}$ in the outer region, respectively. Combined with the Square Kilometer Array (KM2A) measurements at higher energies, a clear softening of the spectrum is found in the inner region, with change of spectral indices by $\sim0.5$ at a break energy around $30$ TeV. The fluxes of the diffuse emission are higher by a factor of $1.5-2.7$ than the model prediction assuming local CR spectra and the gas column density, which are consistent with those measured by the KM2A. Along Galactic longitude, the spatial distribution of the diffuse emission shows deviation from that of the gas column density. The spectral shape of the diffuse emission are possibly variation in different longitude region. The WCDA measurements bridge the gap between the low-energy measurements by space detectors and the ultra-high-energy observations by LHAASO-KM2A and other experiments. These results suggest that improved modeling of the wide-band diffuse emission is required.

preprint2026arXiv

Transient Large-Scale Anisotropy in TeV Cosmic Rays due to an Interplanetary Coronal Mass Ejection

Large- or medium-scale cosmic ray anisotropy at TeV energies has not previously been confirmed to vary with time. Transient anisotropy changes have been observed below 150 GeV, especially near the passage of an interplanetary shock and coronal mass ejection containing a magnetic flux rope ejected by a solar storm, which can trigger a geomagnetic storm with practical consequences. In such events, cosmic rays provide remote sensing of the magnetic field properties. Here we report the observation of transient large-scale anisotropy in TeV cosmic ray ions using data from the Large High Altitude Air Shower Observatory (LHAASO). We analyze hourly skymaps of the transient cosmic ray intensity excess or deficit, the gradient of which indicates the direction and magnitude of transient large-scale anisotropy across the field of view. We observe enhanced anisotropy above typical hourly fluctuations with $>$5$σ$ significance during some hours of November 4, 2021, in separate data sets for four primary cosmic ray energy ranges of median energy from $E$=0.7 to 3.1 TeV. The gradient varies with energy as $E^γ$, where $γ\approx-0.5$. At a median energy $\leq$1.0 TeV, this gradient corresponds to a dipole anisotropy of at least 1\%, or possibly a weaker anisotropy of higher order. This new type of observation opens the opportunity to study interplanetary magnetic structures using air shower arrays around the world, complementing existing in situ and remote measurements of plasma properties.

preprint2025arXiv

Observation of the $γ$-ray Emission from W43 with LHAASO

In this paper, we report the detection of the very-high-energy (VHE, $ 100{\rm\ GeV} < E < 100{\rm\ TeV} $) and ultra-high-energy (UHE, $E > 100\rm\ TeV$) $γ$-ray emissions from the direction of the young star-forming region W43, observed by the Large High Altitude Air Shower Observation (LHAASO). The extended $γ$-ray source was detected with a significance of ${\sim}16\,σ$ by KM2A and ${\sim}17\,σ$ by WCDA, respectively. The angular extension of this $γ$-ray source is about 0.5 degrees, corresponding to a physical size of about 50 pc. We discuss the origin of the $γ$-ray emission and possible cosmic ray acceleration in the W43 region using multi-wavelength data. Our findings suggest that W43 is likely another young star cluster capable of accelerating cosmic rays (CRs) to at least several hundred TeV.

preprint2025arXiv

Study of Ultra-High-Energy Gamma-Ray Source 1LHAASO J0056+6346u and Its Possible Origins

We report a dedicated study of the newly discovered extended UHE $γ$-ray source 1LHAASO J0056+6346u. Analyzing 979 days of LHAASO-WCDA data and 1389 days of LHAASO-KM2A data, we observed a significant excess of $γ$-ray events with both WCDA and KM2A. Assuming a point power-law source with a fixed spectral index, the significance maps reveal excesses of ${\sim}12.65\,σ$, ${\sim}22.18\,σ$, and ${\sim}10.24\,σ$ in the energy ranges of 1--25 TeV, 25--100 TeV, and $> 100$ TeV, respectively. We use a 3D likelihood algorithm to derive the morphological and spectral parameters, and the source is detected with significances of $12.65\,σ$ by WCDA and $25.27\,σ$ by KM2A. The best-fit positions derived from WCDA and KM2A data are (R.A. = $13.96^\circ\pm0.09^\circ$, Decl. = $63.92^\circ\pm0.05^\circ$) and (R.A. = $14.00^\circ\pm0.05^\circ$, Decl. = $63.79^\circ\pm0.02^\circ$), respectively. The angular size ($r_{39}$) of 1LHAASO J0056+6346u is $0.34^\circ\pm0.04^\circ$ at 1--25 TeV and $0.24^\circ\pm0.02^\circ$ at $> 25$ TeV. The differential flux of this UHE $γ$-ray source can be described by an exponential cutoff power-law function: $(2.67\pm0.25) \times 10^{-15} (E/20\,\text{TeV})^{-1.97\pm0.10} e^{-E/(55.1\pm7.2)\,\text{TeV}} \,\text{TeV}^{-1}\,\text{cm}^{-2}\,\text{s}^{-1}$. To explore potential sources of $γ$-ray emission, we investigated the gas distribution around 1LHAASO J0056+6346u. 1LHAASO J0056+6346u is likely to be a TeV PWN powered by an unknown pulsar, which would naturally explain both its spatial and spectral properties. Another explanation is that this UHE $γ$-ray source might be associated with gas content illuminated by a nearby CR accelerator, possibly the SNR candidate G124.0+1.4.

preprint2023arXiv

Longitudinal and transverse spin transfer to $Λ$ and $\overlineΛ$ hyperons in polarized $p$+$p$ collisions at $\sqrt{s} = 200$ GeV

The longitudinal and transverse spin transfers to $Λ$ ($\overlineΛ$) hyperons in polarized proton-proton collisions are expected to be sensitive to the helicity and transversity distributions, respectively, of (anti-)strange quarks in the proton, and to the corresponding polarized fragmentation functions. We report improved measurements of the longitudinal spin transfer coefficient, $D_{LL}$, and the transverse spin transfer coefficient, $D_{TT}$, to $Λ$ and $\overlineΛ$ in polarized proton-proton collisions at $\sqrt{s}$ = 200 GeV by the STAR experiment at RHIC. The data set includes longitudinally polarized proton-proton collisions with an integrated luminosity of 52 pb$^{-1}$, and transversely polarized proton-proton collisions with a similar integrated luminosity. Both data sets have about twice the statistics of previous results and cover a kinematic range of $|η_{Λ(\overlineΛ)}|$ $<$ 1.2 and transverse momentum $p_{T,{Λ(\overlineΛ)}}$ up to 8 GeV/$c$. We also report the first measurements of the hyperon spin transfer coefficients $D_{LL}$ and $D_{TT}$ as a function of the fractional jet momentum $z$ carried by the hyperon, which can provide more direct constraints on the polarized fragmentation functions.

preprint2023arXiv

Model-based cross-correlation search for gravitational waves from the low-mass X-ray binary Scorpius X-1 in LIGO O3 data

We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to balance sensitivity with computing cost. The search covered a range of gravitational-wave frequencies from 25Hz to 1600Hz, as well as ranges in orbital speed, frequency and phase determined from observational constraints. No significant detection candidates were found, and upper limits were set as a function of frequency. The most stringent limits, between 100Hz and 200Hz, correspond to an amplitude h0 of about 1e-25 when marginalized isotropically over the unknown inclination angle of the neutron star's rotation axis, or less than 4e-26 assuming the optimal orientation. The sensitivity of this search is now probing amplitudes predicted by models of torque balance equilibrium. For the usual conservative model assuming accretion at the surface of the neutron star, our isotropically-marginalized upper limits are close to the predicted amplitude from about 70Hz to 100Hz; the limits assuming the neutron star spin is aligned with the most likely orbital angular momentum are below the conservative torque balance predictions from 40Hz to 200Hz. Assuming a broader range of accretion models, our direct limits on gravitational-wave amplitude delve into the relevant parameter space over a wide range of frequencies, to 500Hz or more.

preprint2023arXiv

Pattern of Global Spin Alignment of $ϕ$ and $K^{*0}$ mesons in Heavy-Ion Collisions

Notwithstanding decades of progress since Yukawa first developed a description of the force between nucleons in terms of meson exchange, a full understanding of the strong interaction remains a major challenge in modern science. One remaining difficulty arises from the non-perturbative nature of the strong force, which leads to the phenomenon of quark confinement at distances on the order of the size of the proton. Here we show that in relativistic heavy-ion collisions, where quarks and gluons are set free over an extended volume, two species of produced vector (spin-1) mesons, namely $ϕ$ and $K^{*0}$, emerge with a surprising pattern of global spin alignment. In particular, the global spin alignment for $ϕ$ is unexpectedly large, while that for $K^{*0}$ is consistent with zero. The observed spin-alignment pattern and magnitude for the $ϕ$ cannot be explained by conventional mechanisms, while a model with a connection to strong force fields, i.e. an effective proxy description within the Standard Model and Quantum Chromodynamics, accommodates the current data. This connection, if fully established, will open a potential new avenue for studying the behaviour of strong force fields.

preprint2022arXiv

Centrality and transverse momentum dependence of higher-order flow harmonics of identified hadrons in Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV

We present high-precision measurements of elliptic, triangular, and quadrangular flow $v_{2}$, $v_{3}$, and $v_{4}$, respectively, at midrapidity ($|η|<1.0$) for identified hadrons $π$, $p$, $K$, $φ$, $K_s$, $Λ$ as a function of centrality and transverse momentum in Au+Au collisions at the center-of-mass energy $\sqrt{s_{\rm NN}}=$ 200 GeV. We observe similar $v_{n}$ trends between light and strange mesons which indicates that the heavier strange quarks flow as strongly as the lighter up and down quarks. The number-of-constituent-quark scaling for $v_{2}$, $v_{3}$, and $v_{4}$ is found to hold within statistical uncertainty for 0-10$\%$, 10-40$\%$ and 40-80$\%$ collision centrality intervals. The results are compared to several viscous hydrodynamic calculations with varying initial conditions, and could serve as an additional constraint to the development of hydrodynamic models.

preprint2022arXiv

Deeply virtual Compton scattering cross section at high Bjorken $x_B$

We report high-precision measurements of the Deeply Virtual Compton Scattering (DVCS) cross section at high values of the Bjorken variable $x_B$. DVCS is sensitive to the Generalized Parton Distributions of the nucleon, which provide a three-dimensional description of its internal constituents. Using the exact analytic expression of the DVCS cross section for all possible polarization states of the initial and final electron and nucleon, and final state photon, we present the first experimental extraction of all four helicity-conserving Compton Form Factors (CFFs) of the nucleon as a function of $x_B$, while systematically including helicity flip amplitudes. In particular, the high accuracy of the present data demonstrates sensitivity to some very poorly known CFFs.

preprint2022arXiv

Development of a low-background neutron detector array

A low-background neutron detector array was developed to measure the cross section of the $^{13}$C($α$,n)$^{16}$O reaction, which is the neutron source for the $s$-process in AGB stars, in the Gamow window ($E_{c.m.}$ = 190 $\pm$ 40 keV) at the China Jinping Underground Laboratory (CJPL). The detector array consists of 24 $^{3}$He proportional counters embedded in a polyethylene cube. Due to the deep underground location and a borated polyethylene shield around the detector array, a low background of 4.5(2)/hour was achieved. The $^{51}$V(p, n)$^{51}$Cr reaction was used to determine the neutron detection efficiency of the array for neutrons with energy $E_n$ $<$ 1 MeV. Geant4 simulations, which were shown to well reproduce experimental results, were used to extrapolate the detection efficiency to higher energies for neutrons emitted in the $^{13}$C($α$,n) $^{16}$O reaction. The theoretical angular distributions of the $^{13}$C($α$,n)$^{16}$O reaction were shown to be important in estimating the uncertainties of the detection efficiency.

preprint2022arXiv

Differential measurements of jet substructure and partonic energy loss in Au$+$Au collisions at $\sqrt{s_{\rm{NN}}} =200$ GeV

The STAR collaboration presents jet substructure measurements related to both the momentum fraction and the opening angle within jets in \pp and \AuAu collisions at \sqrtsn $= 200$ GeV. The substructure observables include SoftDrop groomed momentum fraction (\zg), groomed jet radius (\rg), and subjet momentum fraction (\zsj) and opening angle (\tsj). The latter observable is introduced for the first time. Fully corrected subjet measurements are presented for \pp collisions and are compared to leading order Monte Carlo models. The subjet \tsj~distributions reflect the jets leading opening angle and are utilized as a proxy for the resolution scale of the medium in \AuAu collisions. We compare data from \AuAu collisions to those from \pp which are embedded in minimum-bias \AuAu events in order to include the effects of detector smearing and the heavy-ion collision underlying event. The subjet observables are shown to be more robust to the background than \zg~and \rg. We observe no significant modifications of the subjet observables within the two highest-energy, back-to-back jets, resulting in a distribution of opening angles and the splittings that are vacuum-like. We also report measurements of the differential di-jet momentum imbalance ($A_{\rm{J}}$) for jets of varying \tsj. We find no qualitative differences in energy loss signatures for varying angular scales in the range $0.1 < $ \tsj $ < 0.3$, leading to the possible interpretation that energy loss in this population of high momentum di-jet pairs, is due to soft medium-induced gluon radiation from a single color-charge as it traverses the medium.

preprint2022arXiv

Direct Observation of Coherent Longitudinal and Shear Acoustic Phonons in TaAs Using Ultrafast X-ray Diffraction

Using femtosecond time-resolved X-ray diffraction, we investigated optically excited coherent acoustic phonons in the Weyl semimetal TaAs. The low symmetry of the (112) surface probed in our experiment enables the simultaneous excitation of longitudinal and shear acoustic modes, whose dispersion closely matches our simulations. We observed an asymmetry in the spectral lineshape of the longitudinal mode that is notably absent from the shear mode, suggesting a time-dependent frequency chirp that is likely driven by photoinduced carrier diffusion. We argue on the basis of symmetry that these acoustic deformations can transiently alter the electronic structure near the Weyl points and support this with model calculations. Our study underscores the benefit of using off-axis crystal orientations when optically exciting acoustic deformations in topological semimetals, allowing one to transiently change their crystal and electronic structures.

preprint2022arXiv

Disappearance of partonic collectivity in $\sqrt{s_{NN}}$ = 3 GeV Au+Au collisions at RHIC

We report on the measurements of directed flow $v_1$ and elliptic flow $v_2$ for hadrons ($π^{\pm}$, $K^{\pm}$, $K_{S}^0$, $p$, $ϕ$, $Λ$ and $Ξ^{-}$) from Au+Au collisions at $\sqrt{s_{NN}}$ = 3\,GeV and $v_{2}$ for ($π^{\pm}$, $K^{\pm}$, $p$ and $\overline{p}$) at 27 and 54.4\,GeV with the STAR experiment. While at the two higher energy midcentral collisions the number-of-constituent-quark (NCQ) scaling holds, at 3\,GeV the $v_{2}$ at midrapidity is negative for all hadrons and the NCQ scaling is absent. In addition, the $v_1$ slopes at midrapidity for almost all observed hadrons are found to be positive, implying dominant repulsive baryonic interactions. The features of negative $v_2$ and positive $v_1$ slope at 3\,GeV can be reproduced with a baryonic mean-field in transport model calculations. These results imply that the medium in such collisions is likely characterized by baryonic interactions.

preprint2022arXiv

Evidence for Nonlinear Gluon Effects in QCD and their $A$ Dependence at STAR

The STAR Collaboration reports measurements of back-to-back azimuthal correlations of di-$π^0$s produced at forward pseudorapidities ($2.6<η<4.0$) in $p$+$p$, $p+$Al, and $p+$Au collisions at a center-of-mass energy of 200 GeV. We observe a clear suppression of the correlated yields of back-to-back $π^0$ pairs in $p+$Al and $p+$Au collisions compared to the $p$+$p$ data. The observed suppression of back-to-back pairs as a function of transverse momentum suggests nonlinear gluon dynamics arising at high parton densities. The larger suppression found in $p+$Au relative to $p+$Al collisions exhibits a dependence of the saturation scale, $Q_s^2$, on the mass number, $A$. A linear scaling of the suppression with $A^{1/3}$ is observed with a slope of $-0.09$ $\pm$ $0.01$.

preprint2022arXiv

Evolution of Accretion Modes between Spectral States Inferred from Spectral and Timing Analysis of Cygnus X-1 with Insight-HXMT observations

We execute a detailed spectral-timing study of Cygnus X-1 in the low/hard, intermediate and high/soft states with the Hard X-ray Modulation Telescope observations. The broad band energy spectra fit well with the "truncated disk model" with the inner boundary of the accretion disk stays within $\sim$10 \textit R$_{\rm g}$ and moves inward as the source softens. Through studying of PDS, rms and Fourier-frequency component resolved spectroscopy, we find that the X-ray variations are generated in two different regions for each state. We discover that the major role that contributes to the X-ray variation is the hot corona rather than the accretion disk. We suggest a scenario with different corona geometry for each state based on the truncated disk geometry in which the corona wraps up the disk to form a sandwich geometry in the low/hard state, and then gradually moves away from the disk in direction that is perpendicular to the disk until forming a jet-like geometry in the high/soft state.

preprint2022arXiv

Fabrication and installation of the Mu2e cryogenic distribution system

The muon-to-electron conversion (Mu2e) experiment at Fermilab will be used to search for the charged lepton flavor-violating conversion of muons to electrons in the field of an atomic nucleus. The Mu2e experiment is currently in the construction stage. The scope of this paper is the cryogenic distribution system and superconducting power leads for four superconducting solenoid magnets: Production Solenoid (PS), an Upstream and Downstream Transport Solenoids (TSu and TSd) and Detector Solenoid (DS). The design of the cryogenic distribution system and the fabrication of several sub-systems was reported previously. This paper reports on additional fabrication and installation progress that has been performed over the past two years. Lessons learned during fabrication and testing of the cryogenic distribution system components are described. In particular, the challenges and solutions implemented for aluminum welding are reported. A description of the process used to qualify the welding procedure and welders for welding the aluminium stabilized NbTi superconducting power leads is provided. Additionally, the progress made with regards to installing the power leads into the cryogenic Feedboxes is covered.

preprint2022arXiv

Global $Λ$-hyperon polarization in Au+Au collisions at $\sqrt{s_\mathrm{NN}}=3$ GeV

Global hyperon polarization, $\overline{P}_\mathrm{H}$, in Au+Au collisions over a large range of collision energy, $\sqrt{s_\mathrm{NN}}$, was recently measured and successfully reproduced by hydrodynamic and transport models with intense fluid vorticity of the quark-gluon plasma. While naïve extrapolation of data trends suggests a large $\overline{P}_\mathrm{H}$ as the collision energy is reduced, the behavior of $\overline{P}_\mathrm{H}$ at small $\sqrt{s_\mathrm{NN}}<7.7$ GeV is unknown. Operating the STAR experiment in fixed-target mode, we measured the polarization of $Λ$ hyperons along the direction of global angular momentum in Au+Au collisions at $\sqrt{s_\mathrm{NN}}=3$ GeV. The observation of substantial polarization of $4.91\pm0.81(\rm stat.)\pm0.15(\rm syst.)$% in these collisions may require a reexamination of the viscosity of any fluid created in the collision, of the thermalization timescale of rotational modes, and of hadronic mechanisms to produce global polarization.

preprint2022arXiv

GWTC-2.1: Deep Extended Catalog of Compact Binary Coalescences Observed by LIGO and Virgo During the First Half of the Third Observing Run

The second Gravitational-Wave Transient Catalog reported on 39 compact binary coalescences observed by the Advanced LIGO and Advanced Virgo detectors between 1 April 2019 15:00 UTC and 1 October 2019 15:00 UTC. We present GWTC-2.1, which reports on a deeper list of candidate events observed over the same period. We analyze the final version of the strain data over this period with improved calibration and better subtraction of excess noise, which has been publicly released. We employ three matched-filter search pipelines for candidate identification, and estimate the astrophysical probability for each candidate event. While GWTC-2 used a false alarm rate threshold of 2 per year, we include in GWTC-2.1, 1201 candidates that pass a false alarm rate threshold of 2 per day. We calculate the source properties of a subset of 44 high-significance candidates that have an astrophysical probability greater than 0.5. Of these candidates, 36 have been reported in GWTC-2. If the 8 additional high-significance candidates presented here are astrophysical, the mass range of events that are unambiguously identified as binary black holes (both objects $\geq 3M_\odot$) is increased compared to GWTC-2, with total masses from $\sim 14 M_\odot$ for GW190924_021846 to $\sim 182 M_\odot$ for GW190426_190642. The primary components of two new candidate events (GW190403_051519 and GW190426_190642) fall in the mass gap predicted by pair instability supernova theory. We also expand the population of binaries with significantly asymmetric mass ratios reported in GWTC-2 by an additional two events (the mass ratio is less than $0.65$ and $0.44$ at $90\%$ probability for GW190403_051519 and GW190917_114630 respectively), and find that 2 of the 8 new events have effective inspiral spins $χ_\mathrm{eff} > 0$ (at $90\%$ credibility), while no binary is consistent with $χ_\mathrm{eff} < 0$ at the same significance.

preprint2022arXiv

Incoherent transport across the strange metal regime of highly overdoped cuprates

Strange metals possess highly unconventional transport characteristics, such as a linear-in-temperature ($T$) resistivity, an inverse Hall angle that varies as $T^2$ and a linear-in-field ($H$) magnetoresistance. Identifying the origin of these collective anomalies has proved profoundly challenging, even in materials such as the hole-doped cuprates that possess a simple band structure. The prevailing dogma is that strange metallicity in the cuprates is tied to a quantum critical point at a doping $p*$ inside the superconducting dome. Here, we study the high-field in-plane magnetoresistance of two superconducting cuprate families at doping levels beyond $p*$. At all dopings, the magnetoresistance exhibits quadrature scaling and becomes linear at high $H/T$ ratios. Moreover, its magnitude is found to be much larger than predicted by conventional theory and insensitive to both impurity scattering and magnetic field orientation. These observations, coupled with analysis of the zero-field and Hall resistivities, suggest that despite having a single band, the cuprate strange metal phase hosts two charge sectors, one containing coherent quasiparticles, the other scale-invariant `Planckian' dissipators.

preprint2022arXiv

Joint User Scheduling and Beamforming Design for Multiuser MISO Downlink Systems

In multiuser communication systems, user scheduling and beamforming (US-BF) design are two fundamental problems that are usually studied separately in the existing literature. In this work, we focus on the joint US-BF design with the goal of maximizing the set cardinality of scheduled users, which is computationally challenging due to the non-convex objective function and the coupled constraints with discrete-continuous variables. To tackle these difficulties, a successive convex approximation based US-BF (SCA-USBF) optimization algorithm is firstly proposed. Then, inspired by wireless intelligent communication, a graph neural network based joint US-BF (J-USBF) learning algorithm is developed by combining the joint US and power allocation network model with the BF analytical solution. The effectiveness of SCA-USBF and J-USBF is verified by various numerical results, the latter achieves close performance and higher computational efficiency. Furthermore, the proposed J-USBF also enjoys the generalizability in dynamic wireless network scenarios.

preprint2022arXiv

Light Nuclei Collectivity from $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au Collisions at RHIC

In high-energy heavy-ion collisions, partonic collectivity is evidenced by the constituent quark number scaling of elliptic flow anisotropy for identified hadrons. A breaking of this scaling and dominance of baryonic interactions is found for identified hadron collective flow measurements in $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions. In this paper, we report measurements of the first- and second-order azimuthal anisotropic parameters, $v_1$ and $v_2$, of light nuclei ($d$, $t$, $^{3}$He, $^{4}$He) produced in $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions at the STAR experiment. An atomic mass number scaling is found in the measured $v_1$ slopes of light nuclei at mid-rapidity. For the measured $v_2$ magnitude, a strong rapidity dependence is observed. Unlike $v_2$ at higher collision energies, the $v_2$ values at mid-rapidity for all light nuclei are negative and no scaling is observed with the atomic mass number. Calculations by the Jet AA Microscopic Transport Model (JAM), with baryonic mean-field plus nucleon coalescence, are in good agreement with our observations, implying baryonic interactions dominate the collective dynamics in 3 GeV Au+Au collisions at RHIC.

preprint2022arXiv

Longitudinal double-spin asymmetry for inclusive jet and dijet production in polarized proton collisions at $\sqrt{s}=510$ GeV

We report measurements of the longitudinal double-spin asymmetry, $A_{LL}$, for inclusive jet and dijet production in polarized proton-proton collisions at midrapidity and center-of-mass energy $\sqrt{s}$ = 510 GeV, using the high luminosity data sample collected by the STAR experiment in 2013. These measurements complement and improve the precision of previous STAR measurements at the same center-of-mass energy that probe the polarized gluon distribution function at partonic momentum fraction 0.015 $\lesssim x \lesssim$ 0.25. The dijet asymmetries are separated into four jet-pair topologies, which provide further constraints on the $x$ dependence of the polarized gluon distribution function. These measurements are in agreement with previous STAR measurements and with predictions from current next-to-leading order global analyses. They provide more precise data at low dijet invariant mass that will better constraint the shape of the polarized gluon distribution function of the proton.

preprint2022arXiv

Measurement of cold nuclear matter effects for inclusive $J/ψ$ in $p$+Au collisions at $\sqrt{s_{_{\mathrm{NN}}}}$ = 200 GeV

Measurement by the STAR experiment at RHIC of the cold nuclear matter (CNM) effects experienced by inclusive $J/ψ$ at mid-rapidity in 0-100\% $p$+Au collisions at $\sqrt{s_{_{\mathrm{NN}}}}$ = 200 GeV is presented. Such effects are quantified utilizing the nuclear modification factor, $R_{p\mathrm{Au}}$, obtained by taking a ratio of $J/ψ$ yield in $p$+Au collisions to that in $p$+$p$ collisions scaled by the number of binary nucleon-nucleon collisions. The differential $J/ψ$ yield in both $p$+$p$ and $p$+Au collisions is measured through the dimuon decay channel, taking advantage of the trigger capability provided by the Muon Telescope Detector in the RHIC 2015 run. Consequently, the $J/ψ$ $R_{p\mathrm{Au}}$ is derived within the transverse momentum ($p_{\mathrm{T}}$) range of 0 to 10 GeV/$c$. A suppression of approximately 30% is observed for $p_{\mathrm{T}}<2$ GeV/$c$, while $J/ψ$ $R_{p\mathrm{Au}}$ becomes compatible with unity for $p_{\mathrm{T}}$ greater than 3 GeV/$c$, indicating the $J/ψ$ yield is minimally affected by the CNM effects at high $p_{\mathrm{T}}$. Comparison to a similar measurement from 0-20% central Au+Au collisions reveals that the observed strong $J/ψ$ suppression above 3 Gev/$c$ is mostly due to the hot medium effects, providing strong evidence for the formation of the quark-gluon plasma in these collisions. Several model calculations show qualitative agreement with the measured $J/ψ$ $R_{p\mathrm{Au}}$, while their agreement with the $J/ψ$ yields in $p$+$p$ and $p$+Au collisions is worse.

preprint2022arXiv

Measurement of inclusive electrons from open heavy-flavor hadron decays in $p$+$p$ collisions at $\sqrt{s} = 200$ GeV with the STAR detector

We report a new measurement of the production cross section for inclusive electrons from open heavy-flavor hadron decays as a function of transverse momentum ($p_{\rm T}$) at mid-rapidity ($|y|<$ 0.7) in $p$+$p$ collisions at $\sqrt{s} = 200$ GeV. The result is presented for 2.5 $<p_{\rm T}<$ 10 GeV/$c$ with an improved precision above 6 GeV/$c$ with respect to the previous measurements, providing more constraints on perturbative QCD calculations. Moreover, this measurement also provides a high-precision reference for measurements of nuclear modification factors for inclusive electrons from open-charm and -bottom hadron decays in heavy-ion collisions.

preprint2022arXiv

Measurements of ${}^3_Λ\rm{H}$ and ${}^4_Λ\rm{H}$ Lifetimes and Yields in Au+Au Collisions in the High Baryon Density Region

We report precision measurements of hypernuclei ${}^3_Λ\rm{H}$ and ${}^4_Λ\rm{H}$ lifetimes obtained from Au+Au collisions at \snn = 3.0\,GeV and 7.2\,GeV collected by the STAR experiment at RHIC, and the first measurement of ${}^3_Λ\rm{H}$ and ${}^4_Λ\rm{H}$ mid-rapidity yields in Au+Au collisions at \snn = 3.0\,GeV. ${}^3_Λ\rm{H}$ and ${}^4_Λ\rm{H}$, being the two simplest bound states composed of hyperons and nucleons, are cornerstones in the field of hypernuclear physics. Their lifetimes are measured to be $221\pm15(\rm stat.)\pm19(\rm syst.)$\,ps for ${}^3_Λ\rm{H}$ and $218\pm6(\rm stat.)\pm13(\rm syst.)$\,ps for ${}^4_Λ\rm{H}$. The $p_T$-integrated yields of ${}^3_Λ\rm{H}$ and ${}^4_Λ\rm{H}$ are presented in different centrality and rapidity intervals. It is observed that the shape of the rapidity distribution of ${}^4_Λ\rm{H}$ is different for 0--10\% and 10--50\% centrality collisions. Thermal model calculations, using the canonical ensemble for strangeness, describes the ${}^3_Λ\rm{H}$ yield well, while underestimating the ${}^4_Λ\rm{H}$ yield. Transport models, combining baryonic mean-field and coalescence (JAM) or utilizing dynamical cluster formation via baryonic interactions (PHQMD) for light nuclei and hypernuclei production, approximately describe the measured ${}^3_Λ\rm{H}$ and ${}^4_Λ\rm{H}$ yields. Our measurements provide means to precisely assess our understanding of the fundamental baryonic interactions with strange quarks, which can impact our understanding of more complicated systems involving hyperons, such as the interior of neutron stars or exotic hypernuclei.

preprint2022arXiv

Measurements of Proton High Order Cumulants in 3 GeV Au+Au Collisions and Implications for the QCD Critical Point

We report cumulants of the proton multiplicity distribution from dedicated fixed-target Au+Au collisions at 3.0 GeV, measured by the STAR experiment in the kinematic acceptance of rapidity ($y$) and transverse momentum ($p_{\rm T}$) within $-0.5 < y<0$ and $0.4 < p_{\rm T} <2.0 $ GeV/$c$. In the most central 0--5\% collisions, a proton cumulant ratio is measured to be $C_4/C_2=-0.85 \pm 0.09 ~(\rm stat.) \pm 0.82 ~(\rm syst.)$, which is less than unity, the Poisson baseline. The hadronic transport UrQMD model reproduces our $C_4/C_2$ in the measured acceptance. Compared to higher energy results and the transport model calculations, the suppression in $C_4/C_2$ is consistent with fluctuations driven by baryon number conservation and indicates an energy regime dominated by hadronic interactions. These data imply that the QCD critical region, if created in heavy-ion collisions, could only exist at energies higher than 3\,GeV.

preprint2022arXiv

Narrowband searches for continuous and long-duration transient gravitational waves from known pulsars in the LIGO-Virgo third observing run

Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational radiation is phase-locked to the electromagnetic emission. In the search presented here, we relax this assumption and allow the frequency and frequency time-derivative of the gravitational waves to vary in a small range around those inferred from electromagnetic observations. We find no evidence for continuous gravitational waves, and set upper limits on the strain amplitude for each target. These limits are more constraining for seven of the targets than the spin-down limit defined by ascribing all rotational energy loss to gravitational radiation. In an additional search we look in O3 data for long-duration (hours-months) transient gravitational waves in the aftermath of pulsar glitches for six targets with a total of nine glitches. We report two marginal outliers from this search, but find no clear evidence for such emission either. The resulting duration-dependent strain upper limits do not surpass indirect energy constraints for any of these targets.

preprint2022arXiv

Peculiar disk behaviors of the black hole candidate MAXI J1348-630 in the hard state observed by Insight-HXMT and Swift

We present a spectral study of the black hole candidate MAXI J1348-630 during its 2019 outburst, based on monitoring observations with Insight-HXMT and Swift. Throughout the outburst, the spectra are well fitted with power-law plus disk-blackbody components. In the soft-intermediate and soft states, we observed the canonical relation L ~ T_in^4 between disk luminosity L and peak colour temperature T_in, with a constant inner radius R_in (traditionally identified with the innermost stable circular orbit). At other stages of the outburst cycle, the behaviour is more unusual, inconsistent with the canonical outburst evolution of black hole transients. In particular, during the hard rise, the apparent inner radius is smaller than in the soft state (and increasing), and the peak colour temperature is higher (and decreasing). This anomalous behaviour is found even when we model the spectra with self-consistent Comptonization models, which take into account the up-scattering of photons from the disk component into the power-law component. To explain both those anomalous trends at the same time, we suggest that the hardening factor for the inner disk emission was larger than the canonical value of ~1.7 at the beginning of the outburst. A more physical trend of radii and temperature evolution requires a hardening factor evolving from ~3.5 at the beginning of the hard state to ~1.7 in the hard intermediate state. This could be evidence that the inner disk was in the process of condensing from the hot, optically thin medium and had not yet reached a sufficiently high optical depth for its emission spectrum to be described by the standard optically-thick disk solution.

preprint2022arXiv

Probing Strangeness Canonical Ensemble with $K^{-}$, $ϕ(1020)$ and $Ξ^{-}$ Production in Au+Au Collisions at ${\sqrt{s_{\rm NN}} = \rm{3\,GeV}}$

We report the first multi-differential measurements of strange hadrons of $K^{-}$, $ϕ$ and $Ξ^{-}$ yields as well as the ratios of $ϕ/K^-$ and $ϕ/Ξ^-$ in Au+Au collisions at ${\sqrt{s_{\rm NN}} = \rm{3\,GeV}}$ with the STAR experiment fixed target configuration at RHIC. The $ϕ$ mesons and $Ξ^{-}$ hyperons are measured through hadronic decay channels, $ϕ\rightarrow K^+K^-$ and $Ξ^-\rightarrow Λπ^-$. Collision centrality and rapidity dependence of the transverse momentum spectra for these strange hadrons are presented. The $4π$ yields and ratios are compared to thermal model and hadronic transport model predictions. At this collision energy, thermal model with grand canonical ensemble (GCE) under-predicts the $ϕ/K^-$ and $ϕ/Ξ^-$ ratios while the result of canonical ensemble (CE) calculations reproduce $ϕ/K^-$, with the correlation length $r_c \sim 2.7$\,fm, and $ϕ/Ξ^-$, $r_c \sim 4.2$\,fm, for the 0-10\% central collisions. Hadronic transport models including high mass resonance decays could also describe the ratios. While thermal calculations with GCE work well for strangeness production in high energy collisions, the change to CE at $\rm{3\,GeV}$ implies a rather different medium property at high baryon density.

preprint2022arXiv

Probing the Galactic Halo with RR Lyrae Stars. II. The Substructures of the Milky Way

We identify substructures of the Galactic halo using 3,003 type $ab$ RR Lyraes (RRab) with 6D position-velocity information from the SDSS, LAMOST, and Gaia EDR3. Based on the information, we define the separation of any two of the stars in the integrals of motion space and identify substructures by utilizing the friends-of-friends algorithm. We identify members belonging to several known substructures: the Sagittarius stream, the Gaia-Enceladus-Sausage (GES), the Sequoia, and the Helmi streams. In addition to these known substructures, there are three other substructures possibly associated with globular clusters NGC 5272, NGC 6656, and NGC 5024, respectively. Finally, we also find three remaining unknown substructures and one of them has large angular momentum and a mean metallicity $\rm -2.13\,dex$ which may be a new substructure. As for GES, we find that it accounts for a large part of substructures in the inner halo and the range of apocenter distance is from 10 to $34\,\rm kpc$, which suggests that the GES is mainly distributed in the inner halo. The near one-third proportion of the GES and the peak value $20\,\rm kpc$ of the apocenter distances suggest that GES could account for the break in the density profile of the Galactic halo at Galactocentric distance ${\sim}20-25\,\rm kpc$. The similarity of comparing the kinematic properties of Gaia-Enceladus-Sausage with the Hercules-Aquila Cloud and Virgo Overdensity suggests that the three substructures may have similar origins.

preprint2022arXiv

Probing the gluonic structure of the deuteron with $J/ψ$ photoproduction in d+Au ultra-peripheral collisions

Understanding gluon density distributions and how they are modified in nuclei are among the most important goals in nuclear physics. In recent years, diffractive vector meson production measured in ultra-peripheral collisions (UPCs) at heavy-ion colliders has provided a new tool for probing the gluon density. In this Letter, we report the first measurement of $J/ψ$ photoproduction off the deuteron in UPCs at the center-of-mass energy $\sqrt{s_{_{\rm NN}}}=200~\rm GeV$ in d$+$Au collisions. The differential cross section as a function of momentum transfer $-t$ is measured. In addition, data with a neutron tagged in the deuteron-going Zero-Degree Calorimeter is investigated for the first time, which is found to be consistent with the expectation of incoherent diffractive scattering at low momentum transfer. Theoretical predictions based on the Color Glass Condensate saturation model and the gluon shadowing model are compared with the data quantitatively. A better agreement with the saturation model has been observed. With the current measurement, the results are found to be directly sensitive to the gluon density distribution of the deuteron and the deuteron breakup, which provides insights into the nuclear gluonic structure.

preprint2022arXiv

Search for continuous gravitational waves from 20 accreting millisecond X-ray pulsars in O3 LIGO data

Results are presented of searches for continuous gravitational waves from 20 accreting millisecond X-ray pulsars with accurately measured spin frequencies and orbital parameters, using data from the third observing run of the Advanced LIGO and Advanced Virgo detectors. The search algorithm uses a hidden Markov model, where the transition probabilities allow the frequency to wander according to an unbiased random walk, while the $\mathcal{J}$-statistic maximum-likelihood matched filter tracks the binary orbital phase. Three narrow sub-bands are searched for each target, centered on harmonics of the measured spin frequency. The search yields 16 candidates, consistent with a false alarm probability of 30% per sub-band and target searched. These candidates, along with one candidate from an additional target-of-opportunity search done for SAX J1808.4$-$3658, which was in outburst during one month of the observing run, cannot be confidently associated with a known noise source. Additional follow-up does not provide convincing evidence that any are a true astrophysical signal. When all candidates are assumed non-astrophysical, upper limits are set on the maximum wave strain detectable at 95% confidence, $h_0^{95\%}$. The strictest constraint is $h_0^{95\%} = 4.7\times 10^{-26}$ from IGR J17062$-$6143. Constraints on the detectable wave strain from each target lead to constraints on neutron star ellipticity and $r$-mode amplitude, the strictest of which are $ε^{95\%} = 3.1\times 10^{-7}$ and $α^{95\%} = 1.8\times 10^{-5}$ respectively. This analysis is the most comprehensive and sensitive search of continuous gravitational waves from accreting millisecond X-ray pulsars to date.

preprint2022arXiv

Search of the Early O3 LIGO Data for Continuous Gravitational Waves from the Cassiopeia A and Vela Jr. Supernova Remnants

We present directed searches for continuous gravitational waves from the neutron stars in the Cassiopeia A (Cas A) and Vela Jr. supernova remnants. We carry out the searches in the LIGO data from the first six months of the third Advanced LIGO and Virgo observing run, using the Weave semi-coherent method, which sums matched-filter detection-statistic values over many time segments spanning the observation period. No gravitational wave signal is detected in the search band of 20--976 Hz for assumed source ages greater than 300 years for Cas A and greater than 700 years for Vela Jr. Estimates from simulated continuous wave signals indicate we achieve the most sensitive results to date across the explored parameter space volume, probing to strain magnitudes as low as ~$6.3\times10^{-26}$ for Cas A and ~$5.6\times10^{-26}$ for Vela Jr. at frequencies near 166 Hz at 95% efficiency.

preprint2022arXiv

The peculiar spectral evolution of the new X-ray transient MAXI J0637-430

We studied the transient Galactic black hole candidate MAXI J0637-430 with data from Insight-HXMT, Swift and XMM-Newton. The broad-band X-ray observations from Insight-HXMT help us constrain the power-law component. MAXI J0637-430 is located at unusually high Galactic latitude; if it belongs to the Galactic thick disk, we suggest a most likely distance <7 kpc. Compared with other black hole transients, MAXI J0637-430 is also unusual for other reasons: a fast transition to the thermal dominant state at the start of the outburst; a low peak temperature and luminosity (we estimate them at ~ 0.7 keV and <0.1 times Eddington, respectively); a short decline timescale; a low soft-to-hard transition luminosity (<0.01 times Eddington). We argue that such properties are consistent with a small binary separation, short binary period (P ~ 2 hr), and low-mass donor star (M2 ~ 0.2 M_sun). Moreover, spectral modelling shows that a single disk-blackbody component is not a good fit to the thermal emission. Soft spectral residuals, and deviations from the standard L_disk ~ T^4 in relation, suggest the need for a second thermal component. We propose and discuss various scenarios for such component, in addition to those presented in previous studies of this source. For example, a gap in the accretion disk between a hotter inner ring near the innermost stable orbit, and a cooler outer disk. Another possibility is that the second thermal component is the thermal plasma emission from an ionized outflow.

preprint2022arXiv

The Second Catalog of Interplanetary Network Localizations of Konus Short Duration Gamma-Ray Bursts

We present the catalog of Interplanetary Network (IPN) localizations for 199 short-duration gamma-ray bursts (sGRBs) detected by the Konus-Wind (KW) experiment between 2011 January 1 and 2021 August 31, which extends the initial sample of IPN localized KW sGRBs (arXiv:1301.3740) to 495 events. We present the most comprehensive IPN localization data on these events, including probability sky maps in HEALPix format.

preprint2022arXiv

Two-particle correlations on transverse rapidity in Au+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV at STAR

Two-particle correlation measurements projected onto two-dimensional, transverse rapidity coordinates ($y_{T1},y_{T2}$), allow access to dynamical properties of the QCD medium produced in relativistic heavy-ion collisions that angular correlation measurements are not sensitive to. We report non-identified charged-particle correlations for Au + Au minimum-bias collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV taken by the STAR experiment at the Relativistic Heavy-Ion Collider (RHIC). Correlations are presented as 2D functions of transverse rapidity for like-sign, unlike-sign and all charged-particle pairs, as well as for particle pairs whose relative azimuthal angles lie on the near-side, the away-side, or at all relative azimuth. The correlations are constructed using charged particles with transverse momentum $p_T \geq 0.15$ GeV/$c$, pseudorapidity from $-$1 to 1, and azimuthal angles from $-π$ to $π$. The significant correlation structures that are observed evolve smoothly with collision centrality. The major correlation features include a saddle shape plus a broad peak with maximum near $y_T \approx 3$, corresponding to $p_T \approx$ 1.5 GeV/$c$. The broad peak is observed in both like- and unlike-sign charge combinations and in near- and away-side relative azimuthal angles. The all-charge, all-azimuth correlation measurements are compared with the theoretical predictions of {\sc hijing} and {\sc epos}. The results indicate that the correlations for peripheral to mid-central collisions can be approximately described as a superposition of nucleon + nucleon collisions with minimal effects from the QCD medium. Strong medium effects are indicated in mid- to most-central collisions.

preprint2021arXiv

Accretion Torque Reversals in GRO J1008-57 Revealed by Insight-HXMT

GRO J1008-57, as a Be/X-ray transient pulsar, is considered to have the highest magnetic field in known neutron star X-ray binary systems. Observational data of the X-ray outbursts in GRO J1008-57 from 2017 to 2020 were collected by the Insight-HXMT satellite. In this work, the spin period of the neutron star in GRO J1008-57 was determined to be about 93.28 seconds in August 2017, 93.22 seconds in February 2018, 93.25 seconds in June 2019 and 93.14 seconds in June 2020. GRO J1008-57 evolved in the spin-up process with a mean rate of $-(2.10\pm 0.05)\times$10$^{-4}$ s/d from 2009 -- 2018, and turned into a spin down process with a rate of $(6.7\pm 0.6)\times$10$^{-5}$ s/d from Feb 2018 to June 2019. During the type II outburst of 2020, GRO J1008-57 had the spin-up torque again. During the torque reversals, the pulse profiles and continuum X-ray spectra did not change significantly, and the cyclotron resonant scattering feature around 80 keV was only detected during the outbursts in 2017 and 2020. Based on the observed mean spin-up rate, we estimated the inner accretion disk radius in GRO J1008-57 (about 1 - 2 times of the Alfvén radius) by comparing different accretion torque models of magnetic neutron stars. During the spin-down process, the magnetic torque should dominate over the matter accreting inflow torque, and we constrained the surface dipole magnetic field $B\geq 6\times 10^{12}$ G for the neutron star in GRO J1008-57, which is consistent with the magnetic field strength obtained by cyclotron line centroid energy.

preprint2021arXiv

All-sky search for long-duration gravitational-wave bursts in the third Advanced LIGO and Advanced Virgo run

After the detection of gravitational waves from compact binary coalescences, the search for transient gravitational-wave signals with less well-defined waveforms for which matched filtering is not well-suited is one of the frontiers for gravitational-wave astronomy. Broadly classified into "short" $ \lesssim 1~$\,s and "long" $ \gtrsim 1~$\,s duration signals, these signals are expected from a variety of astrophysical processes, including non-axisymmetric deformations in magnetars or eccentric binary black hole coalescences. In this work, we present a search for long-duration gravitational-wave transients from Advanced LIGO and Advanced Virgo's third observing run from April 2019 to March 2020. For this search, we use minimal assumptions for the sky location, event time, waveform morphology, and duration of the source. The search covers the range of $2~\text{--}~ 500$~s in duration and a frequency band of $24 - 2048$ Hz. We find no significant triggers within this parameter space; we report sensitivity limits on the signal strength of gravitational waves characterized by the root-sum-square amplitude $h_{\mathrm{rss}}$ as a function of waveform morphology. These $h_{\mathrm{rss}}$ limits improve upon the results from the second observing run by an average factor of 1.8.

preprint2021arXiv

All-sky search for short gravitational-wave bursts in the third Advanced LIGO and Advanced Virgo run

This paper presents the results of a search for generic short-duration gravitational-wave transients in data from the third observing run of Advanced LIGO and Advanced Virgo. Transients with durations of milliseconds to a few seconds in the 24--4096 Hz frequency band are targeted by the search, with no assumptions made regarding the incoming signal direction, polarization or morphology. Gravitational waves from compact binary coalescences that have been identified by other targeted analyses are detected, but no statistically significant evidence for other gravitational wave bursts is found. Sensitivities to a variety of signals are presented. These include updated upper limits on the source rate-density as a function of the characteristic frequency of the signal, which are roughly an order of magnitude better than previous upper limits. This search is sensitive to sources radiating as little as $\sim$10$^{-10} M_{\odot} c^2$ in gravitational waves at $\sim$70 Hz from a distance of 10~kpc, with 50\% detection efficiency at a false alarm rate of one per century. The sensitivity of this search to two plausible astrophysical sources is estimated: neutron star f-modes, which may be excited by pulsar glitches, as well as selected core-collapse supernova models.

preprint2021arXiv

All-sky, all-frequency directional search for persistent gravitational-waves from Advanced LIGO's and Advanced Virgo's first three observing runs

We present the first results from an all-sky all-frequency (ASAF) search for an anisotropic stochastic gravitational-wave background using the data from the first three observing runs of the Advanced LIGO and Advanced Virgo detectors. Upper limit maps on broadband anisotropies of a persistent stochastic background were published for all observing runs of the LIGO-Virgo detectors. However, a broadband analysis is likely to miss narrowband signals as the signal-to-noise ratio of a narrowband signal can be significantly reduced when combined with detector output from other frequencies. Data folding and the computationally efficient analysis pipeline, {\tt PyStoch}, enable us to perform the radiometer map-making at every frequency bin. We perform the search at 3072 {\tt{HEALPix}} equal area pixels uniformly tiling the sky and in every frequency bin of width $1/32$~Hz in the range $20-1726$~Hz, except for bins that are likely to contain instrumental artefacts and hence are notched. We do not find any statistically significant evidence for the existence of narrowband gravitational-wave signals in the analyzed frequency bins. Therefore, we place $95\%$ confidence upper limits on the gravitational-wave strain for each pixel-frequency pair, the limits are in the range $(0.030 - 9.6) \times10^{-24}$. In addition, we outline a method to identify candidate pixel-frequency pairs that could be followed up by a more sensitive (and potentially computationally expensive) search, e.g., a matched-filtering-based analysis, to look for fainter nearly monochromatic coherent signals. The ASAF analysis is inherently independent of models describing any spectral or spatial distribution of power. We demonstrate that the ASAF results can be appropriately combined over frequencies and sky directions to successfully recover the broadband directional and isotropic results.

preprint2021arXiv

Determining impact parameters of heavy-ion collisions at low-intermediate incident energies using deep learning with convolutional neural network

A deep learning based method with the convolutional neural network (CNN) algorithm for determining the impact parameters is developed using the constrained molecular dynamics model simulations, focusing on the heavy-ion collisions at the low-intermediate incident energies from several ten to one hundred MeV/nucleon in which the emissions of heavy fragments with the charge numbers larger than 3 become crucial. To make the CNN applicable in the task of the impact parameter determination at the present energy range, specific improvements are made in the input selection, the CNN construction and the CNN training. It is demonstrated from the comparisons of the deep CNN method and the conventional methods with the impact parameter-sensitive observables, that the deep CNN method shows better performance for determining the impact parameters, especially leading to the capability of providing better recognition of the central collision events. With a proper consideration of the experimental filter effect in both training and testing processes to keep consistency with the actual experiments, the good performance of the deep CNN method holds, and shows significantly better in terms of predicting the impact parameters and recognizing the central collision events, compared to that of the conventional methods, demonstrating the superiority of the present deep CNN method. The deep CNN method with the consideration of the filter effect is applied in the deduction of nuclear stopping power. Higher accuracy for the stopping power deduction is achieved benefitting from the better impact parameter determination using the deep CNN method, compared to using the the conventional methods. This result reveals the importance to select a reliable impact parameter determination method in the experimental deduction of the nuclear stopping power as well as other observables.

preprint2021arXiv

GECAM detection of a bright type-I X-ray burst from 4U 0614+09: confirmation its spin frequency at 413 Hz

One month after launching Gravitational wave high-energy Electromagnetic Counterpart All-sky Monitor (GECAM), a bright thermonuclear X-ray burst from 4U~0614+09, was observed on January 24, 2021. We report the time-resolved spectroscopy of the burst and a burst oscillation detection at 413 Hz with a fractional amplitude 3.4\% (rms). This coincides with the burst oscillation previously discovered with \textit{Swift}/BAT \citep{Strohmayer2008}, and therefore confirms the spin frequency of this source. This burst is the brightest one in the normal bursts (except the superburst) ever detected from 4U~0614+09, which leads to an upper limit of distance estimation as 3.1 kpc. The folded light curve during the burst oscillation shows a multi-peak structure, which is the first case observed during a single burst oscillation in nonpulsating sources. The multi-peak profile could be due to additional harmonics of the burst oscillation, which is corresponding to several brighter/fainter spots at the stellar surface.

preprint2021arXiv

Open data from the first and second observing runs of Advanced LIGO and Advanced Virgo

Advanced LIGO and Advanced Virgo are actively monitoring the sky and collecting gravitational-wave strain data with sufficient sensitivity to detect signals routinely. In this paper we describe the data recorded by these instruments during their first and second observing runs. The main data products are the gravitational-wave strain arrays, released as time series sampled at 16384 Hz. The datasets that include this strain measurement can be freely accessed through the Gravitational Wave Open Science Center at http://gw-openscience.org, together with data-quality information essential for the analysis of LIGO and Virgo data, documentation, tutorials, and supporting software.

preprint2021arXiv

Overview of KAGRA: Calibration, detector characterization, physical environmental monitors, and the geophysics interferometer

KAGRA is a newly built gravitational wave observatory, a laser interferometer with a 3 km arm length, located in Kamioka, Gifu, Japan. In this series of articles, we present an overview of the baseline KAGRA, for which we finished installing the designed configuration in 2019. This article describes the method of calibration (CAL) used for reconstructing gravitational wave signals from the detector outputs, as well as the characterization of the detector (DET). We also review the physical environmental monitors (PEM) system and the geophysics interferometer (GIF). Both are used for characterizing and evaluating the data quality of the gravitational wave channel. They play important roles in utilizing the detector output for gravitational wave searches. These characterization investigations will be even more important in the near future, once gravitational wave detection has been achieved, and in using KAGRA in the gravitational wave astronomy era.

preprint2021arXiv

Search for Gravitational Waves Associated with Gamma-Ray Bursts Detected by Fermi and Swift During the LIGO-Virgo Run O3b

We search for gravitational-wave signals associated with gamma-ray bursts detected by the Fermi and Swift satellites during the second half of the third observing run of Advanced LIGO and Advanced Virgo (1 November 2019 15:00 UTC-27 March 2020 17:00 UTC).We conduct two independent searches: a generic gravitational-wave transients search to analyze 86 gamma-ray bursts and an analysis to target binary mergers with at least one neutron star as short gamma-ray burst progenitors for 17 events. We find no significant evidence for gravitational-wave signals associated with any of these gamma-ray bursts. A weighted binomial test of the combined results finds no evidence for sub-threshold gravitational wave signals associated with this GRB ensemble either. We use several source types and signal morphologies during the searches, resulting in lower bounds on the estimated distance to each gamma-ray burst. Finally, we constrain the population of low luminosity short gamma-ray bursts using results from the first to the third observing runs of Advanced LIGO and Advanced Virgo. The resulting population is in accordance with the local binary neutron star merger rate.

preprint2021arXiv

Search for lensing signatures in the gravitational-wave observations from the first half of LIGO-Virgo's third observing run

We search for signatures of gravitational lensing in the gravitational-wave signals from compact binary coalescences detected by Advanced LIGO and Advanced Virgo during O3a, the first half of their third observing run. We study: 1) the expected rate of lensing at current detector sensitivity and the implications of a non-observation of strong lensing or a stochastic gravitational-wave background on the merger-rate density at high redshift; 2) how the interpretation of individual high-mass events would change if they were found to be lensed; 3) the possibility of multiple images due to strong lensing by galaxies or galaxy clusters; and 4) possible wave-optics effects due to point-mass microlenses. Several pairs of signals in the multiple-image analysis show similar parameters and, in this sense, are nominally consistent with the strong lensing hypothesis. However, taking into account population priors, selection effects, and the prior odds against lensing, these events do not provide sufficient evidence for lensing. Overall, we find no compelling evidence for lensing in the observed gravitational-wave signals from any of these analyses.

preprint2021arXiv

Search for subsolar-mass binaries in the first half of Advanced LIGO and Virgo's third observing run

We report on a search for compact binary coalescences where at least one binary component has a mass between 0.2 $M_\odot$ and 1.0 $M_\odot$ in Advanced LIGO and Advanced Virgo data collected between 1 April 2019 1500 UTC and 1 October 2019 1500 UTC. We extend previous analyses in two main ways: we include data from the Virgo detector and we allow for more unequal mass systems, with mass ratio $q \geq 0.1$. We do not report any gravitational-wave candidates. The most significant trigger has a false alarm rate of 0.14 $\mathrm{yr}^{-1}$. This implies an upper limit on the merger rate of subsolar binaries in the range $[220-24200] \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$, depending on the chirp mass of the binary. We use this upper limit to derive astrophysical constraints on two phenomenological models that could produce subsolar-mass compact objects. One is an isotropic distribution of equal-mass primordial black holes. Using this model, we find that the fraction of dark matter in primordial black holes is $f_\mathrm{PBH} \equiv Ω_\mathrm{PBH} / Ω_\mathrm{DM} \lesssim 6\%$. The other is a dissipative dark matter model, in which fermionic dark matter can collapse and form black holes. The upper limit on the fraction of dark matter black holes depends on the minimum mass of the black holes that can be formed: the most constraining result is obtained at $M_\mathrm{min}=1 M_\odot$, where $f_\mathrm{DBH} \equiv Ω_\mathrm{PBH} / Ω_\mathrm{DM} \lesssim 0.003\%$. These are the tightest limits on spinning subsolar-mass binaries to date.

preprint2020arXiv

A catalogue of H$α$ emission-line point sources in the vicinity fields of M31 and M33 from the LAMOST survey

We present a catalogue of 3,305 H$α$ emission-line point sources observed with the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) in the vicinity fields of M31 and M33 during September 2011 and January 2016. The catalogue contains 1,487 emission-line stars, 532 emission-line nebulae including 377 likely planetary nebulae (PNe), 83 H~{\textsc{ii}} regions candidates and 20 possible supernovae remnants (SNRs) and 1,286 unknown objects. Among them, 24 PN candidates, 19 H~{\sc ii} region candidates, 10 SNR candidates and 1 symbiotic star candidate are new discoveries. Radial velocities and fluxes estimated from the H$α$ line and those quantities of seven other major emission lines including H$β$, [O~{\textsc{iii}}]~$λ$4959, [O~{\textsc{iii}}]~$λ$5007, [N~{\textsc{ii}}]~$λ$6548, [N~{\textsc{ii}}]~$λ$6583, [S~{\textsc{ii}}]~$λ$6717 and [S~{\textsc{ii}}]~$λ$6731 lines of all the catalogued sources yielded from the LAMOST spectra are also presented in our catalogue. Our catalogue is an ideal starting point to study the chemistry properties and kinematics of M31 and M33.

preprint2020arXiv

A Joint Fermi-GBM and LIGO/Virgo Analysis of Compact Binary Mergers From the First and Second Gravitational-wave Observing Runs

We present results from offline searches of Fermi Gamma-ray Burst Monitor (GBM) data for gamma-ray transients coincident with the compact binary coalescences observed by the gravitational-wave (GW) detectors Advanced LIGO and Advanced Virgo during their first and second observing runs. In particular, we perform follow-up for both confirmed events and low significance candidates reported in the LIGO/Virgo catalog GWTC-1. We search for temporal coincidences between these GW signals and GBM triggered gamma-ray bursts (GRBs). We also use the GBM Untargeted and Targeted subthreshold searches to find coincident gamma-rays below the on-board triggering threshold. This work implements a refined statistical approach by incorporating GW astrophysical source probabilities and GBM visibilities of LIGO/Virgo sky localizations to search for cumulative signatures of coincident subthreshold gamma-rays. All search methods recover the short gamma-ray burst GRB 170817A occurring ~1.7 s after the binary neutron star merger GW170817. We also present results from a new search seeking GBM counterparts to LIGO single-interferometer triggers. This search finds a candidate joint event, but given the nature of the GBM signal and localization, as well as the high joint false alarm rate of $1.1 \times 10^{-6}$ Hz, we do not consider it an astrophysical association. We find no additional coincidences.

preprint2020arXiv

A large catalogue of molecular clouds with accurate distances within 4 kpc of the Galactic disk

We present a large, homogeneous catalogue of molecular clouds within 4 kpc from the Sun at low Galactic latitudes ($|b|$ $<$ 10\degr) with unprecedented accurate distance determinations. Based on the three-dimensional dust reddening map and estimates of colour excesses and distances of over 32 million stars presented in Chen et al, we have identified 567 dust/molecular clouds with a hierarchical structure identification method and obtained their distance estimates by a dust model fitting algorithm. The typical distance uncertainty is less than 5 per cent. As far as we know, this is the first large catalogue of molecular clouds in the Galactic plane with distances derived in a direct manner. The clouds are seen to lie along the Sagittarius, Local and Perseus Arms. In addition to the known structures, we propose the existence of a possible {\it spur}, with a pitch angle of about 34\degr, connecting the Local and the Sagittarius Arms in the fourth quadrant. We have also derived the physical properties of those molecular clouds. The distribution of cloud properties in different parameter spaces agrees grossly with the previous results. Our cloud sample is an ideal starting point to study the concentration of dust and gas in the solar vicinity and their star formation activities.

preprint2020arXiv

A search for prompt gamma-ray counterparts to fast radio bursts in the Insight-HXMT data

No robust detection of prompt electromagnetic counterparts to fast radio bursts (FRBs) has yet been obtained, in spite of several multi-wavelength searches carried out so far. Specifically, X/gamma-ray counterparts are predicted by some models. We planned on searching for prompt gamma-ray counterparts in the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) data, taking advantage of the unique combination of large effective area in the keV-MeV energy range and of sub-ms time resolution. We selected 39 FRBs that were promptly visible from the High-Energy (HE) instrument aboard Insight-HXMT. After calculating the expected arrival times at the location of the spacecraft, we searched for a significant excess in both individual and cumulative time profiles over a wide range of time resolutions, from several seconds down to sub-ms scales. Using the dispersion measures in excess of the Galactic terms, we estimated the upper limits on the redshifts. No convincing signal was found and for each FRB we constrained the gamma-ray isotropic-equivalent luminosity and the released energy as a function of emission timescale. For the nearest FRB source, the periodic repeater FRB180916.J0158+65, we find $L_{γ,iso}<5.5\times 10^{47}$ erg/s over 1 s, whereas $L_{γ,iso}<10^{49}-10^{51}$ erg/s for the bulk of FRBs. The same values scale up by a factor of ~100 for a ms-long emission. Even on a timescale comparable with that of the radio pulse itself no keV-MeV emission is observed. A systematic association with either long or short GRBs is ruled out with high confidence, except for subluminous events, as is the case for core-collapse of massive stars (long) or binary neutron star mergers (short) viewed off axis. Only giant flares from extra-galactic magnetars at least ten times more energetic than Galactic siblings are ruled out for the nearest FRB.

preprint2020arXiv

A variable ionized disk wind in the black-hole candidate EXO 1846-031

After 34 years, the black-hole candidate EXO 1846-031 went into outburst again in 2019. We investigate its spectral properties in the hard intermediate and the soft states with NuSTAR and Insight-HXMT. A reflection component has been detected in the two spectral states but possibly originating from different illumination spectra: in the intermediate state, the illuminating source is attributed to a hard coronal component, which has been commonly observed in other X-ray binaries, whereas in the soft state the reflection is probably produced by the disk self-irradiation. Both cases support EXO 1846-031 as a low inclination system of ~40 degrees. An absorption line is clearly detected at ~7.2 keV in the hard intermediate state, corresponding to a highly ionized disk wind (log ξ > 6.1) with a velocity up to 0.06c. Meanwhile, quasi-simultaneous radio emissions have been detected before and after the X-rays, implying the co-existence of disk winds and jets in this system. Additionally, the observed wind in this source is potentially driven by magnetic forces. The absorption line disappeared in the soft state and a narrow emission line appeared at ~6.7 keV on top of the reflection component, which may be evidence for disk winds, but data with the higher spectral resolution are required to examine this.

preprint2020arXiv

Application of the independent component analysis to the iKAGRA data

We apply the independent component analysis (ICA) to the real data from a gravitational wave detector for the first time. Specifically we use the iKAGRA data taken in April 2016, and calculate the correlations between the gravitational wave strain channel and 35 physical environmental channels. Using a couple of seismic channels which are found to be strongly correlated with the strain, we perform ICA. Injecting a sinusoidal continuous signal in the strain channel, we find that ICA recovers correct parameters with enhanced signal-to-noise ratio, which demonstrates usefulness of this method. Among the two implementations of ICA used here, we find the correlation method yields the optimal result for the case environmental noises act on the strain channel linearly.

preprint2020arXiv

Beam test results of IHEP-NDL Low Gain Avalanche Detectors(LGAD)

To meet the timing resolution requirement of up-coming High Luminosity LHC (HL-LHC), a new detector based on the Low-Gain Avalanche Detector(LGAD), High-Granularity Timing Detector (HGTD), is under intensive research in ATLAS. Two types of IHEP-NDL LGADs(BV60 and BV170) for this update is being developed by Institute of High Energy Physics (IHEP) of Chinese Academic of Sciences (CAS) cooperated with Novel Device Laboratory (NDL) of Beijing Normal University and they are now under detailed study. These detectors are tested with $5GeV$ electron beam at DESY. A SiPM detector is chosen as a reference detector to get the timing resolution of LGADs. The fluctuation of time difference between LGAD and SiPM is extracted by fitting with a Gaussian function. Constant fraction discriminator (CFD) method is used to mitigate the effect of time walk. The timing resolution of $41 \pm 1 ps$ and $63 \pm 1 ps$ are obtained for BV60 and BV170 respectively.

preprint2020arXiv

Discovery of delayed spin-up behavior following two large glitches in the Crab pulsar, and the statistics of such processes

Glitches correspond to sudden jumps of rotation frequency ($ν$) and its derivative ($\dotν$) of pulsars, the origin of which remains not well understood yet, partly because the jump processes of most glitches are not well time-resolved. There are three large glitches of the Crab pulsar, detected in 1989, 1996 and 2017, which were found to have delayed spin-up processes before the normal recovery processes. Here we report two additional glitches of the Crab pulsar occurred in 2004 and 2011 for which we discovered delayed spin up processes, and present refined parameters of the largest glitch occurred in 2017. The initial rising time of the glitch is determined as $<0.48$ hour. We also carried out a statistical study of these five glitches with observed spin-up processes. The two glitches occurred in 2004 and 2011 have delayed spin-up time scales ($τ_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $Δν$ vs. $|Δ{\dotν}|$ relation of these five glitches is similar to those with no detected delayed spin-up process, indicating that they are similar to the others in nature except that they have larger amplitudes. For these five glitches, the amplitudes of the delayed spin-up process ($|Δν_{\rm d1}|$) and recovery process ($Δν_{\rm d2}$), their time scales ($τ_{1}$, $τ_{2}$), and permanent changes in spin frequency ($Δν_{\rm p}$) and total frequency step ($Δν_{\rm g}$) have positive correlations. From these correlations, we suggest that the delayed spin-up processes are common for all glitches, but are too short and thus difficult to be detected for most glitches.

preprint2020arXiv

GW190412: Observation of a Binary-Black-Hole Coalescence with Asymmetric Masses

We report the observation of gravitational waves from a binary-black-hole coalescence during the first two weeks of LIGO's and Virgo's third observing run. The signal was recorded on April 12, 2019 at 05:30:44 UTC with a network signal-to-noise ratio of 19. The binary is different from observations during the first two observing runs most notably due to its asymmetric masses: a ~30 solar mass black hole merged with a ~8 solar mass black hole companion. The more massive black hole rotated with a dimensionless spin magnitude between 0.22 and 0.60 (90% probability). Asymmetric systems are predicted to emit gravitational waves with stronger contributions from higher multipoles, and indeed we find strong evidence for gravitational radiation beyond the leading quadrupolar order in the observed signal. A suite of tests performed on GW190412 indicates consistency with Einstein's general theory of relativity. While the mass ratio of this system differs from all previous detections, we show that it is consistent with the population model of stellar binary black holes inferred from the first two observing runs.

preprint2020arXiv

GW190425: Observation of a Compact Binary Coalescence with Total Mass $\sim 3.4 M_{\odot}$

On 2019 April 25, the LIGO Livingston detector observed a compact binary coalescence with signal-to-noise ratio 12.9. The Virgo detector was also taking data that did not contribute to detection due to a low signal-to-noise ratio, but were used for subsequent parameter estimation. The 90% credible intervals for the component masses range from 1.12 to 2.52 $M_{\odot}$ (1.45 to 1.88 $M_{\odot}$ if we restrict the dimensionless component spin magnitudes to be smaller than 0.05). These mass parameters are consistent with the individual binary components being neutron stars. However, both the source-frame chirp mass $1.44^{+0.02}_{-0.02} M_{\odot}$ and the total mass $3.4^{+0.3}_{-0.1}\,M_{\odot}$ of this system are significantly larger than those of any other known binary neutron star system. The possibility that one or both binary components of the system are black holes cannot be ruled out from gravitational-wave data. We discuss possible origins of the system based on its inconsistency with the known Galactic binary neutron star population. Under the assumption that the signal was produced by a binary neutron star coalescence, the local rate of neutron star mergers is updated to $250-2810 \text{Gpc}^{-3}\text{yr}^{-1}$.

preprint2020arXiv

GW190521: A Binary Black Hole Merger with a Total Mass of $150 ~ M_{\odot}$

On May 21, 2019 at 03:02:29 UTC Advanced LIGO and Advanced Virgo observed a short duration gravitational-wave signal, GW190521, with a three-detector network signal-to-noise ratio of 14.7, and an estimated false-alarm rate of 1 in 4900 yr using a search sensitive to generic transients. If GW190521 is from a quasicircular binary inspiral, then the detected signal is consistent with the merger of two black holes with masses of $85^{+21}_{-14} M_{\odot}$ and $66^{+17}_{-18} M_{\odot}$ (90 % credible intervals). We infer that the primary black hole mass lies within the gap produced by (pulsational) pair-instability supernova processes, and has only a 0.32 % probability of being below $65 M_{\odot}$. We calculate the mass of the remnant to be $142^{+28}_{-16} M_{\odot}$, which can be considered an intermediate mass black hole (IMBH). The luminosity distance of the source is $5.3^{+2.4}_{-2.6}$ Gpc, corresponding to a redshift of $0.82^{+0.28}_{-0.34}$. The inferred rate of mergers similar to GW190521 is $0.13^{+0.30}_{-0.11}\,\mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}$.

preprint2020arXiv

GW190814: Gravitational Waves from the Coalescence of a 23 M$_\odot$ Black Hole with a 2.6 M$_\odot$ Compact Object

We report the observation of a compact binary coalescence involving a 22.2 - 24.3 $M_{\odot}$ black hole and a compact object with a mass of 2.50 - 2.67 $M_{\odot}$ (all measurements quoted at the 90$\%$ credible level). The gravitational-wave signal, GW190814, was observed during LIGO's and Virgo's third observing run on August 14, 2019 at 21:10:39 UTC and has a signal-to-noise ratio of 25 in the three-detector network. The source was localized to 18.5 deg$^2$ at a distance of $241^{+41}_{-45}$ Mpc; no electromagnetic counterpart has been confirmed to date. The source has the most unequal mass ratio yet measured with gravitational waves, $0.112^{+0.008}_{-0.009}$, and its secondary component is either the lightest black hole or the heaviest neutron star ever discovered in a double compact-object system. The dimensionless spin of the primary black hole is tightly constrained to $\leq 0.07$. Tests of general relativity reveal no measurable deviations from the theory, and its prediction of higher-multipole emission is confirmed at high confidence. We estimate a merger rate density of 1-23 Gpc$^{-3}$ yr$^{-1}$ for the new class of binary coalescence sources that GW190814 represents. Astrophysical models predict that binaries with mass ratios similar to this event can form through several channels, but are unlikely to have formed in globular clusters. However, the combination of mass ratio, component masses, and the inferred merger rate for this event challenges all current models for the formation and mass distribution of compact-object binaries.

preprint2020arXiv

Insight-HXMT firm detection of the highest energy fundamental cyclotron resonance scattering feature in the spectrum of GRO J1008-57

We report on the observation of the accreting pulsar GRO J1008-57 performed by Insight-HXMT at the peak of the source's 2017 outburst. Pulsations are detected with a spin period of 93.283(1) s. The pulse profile shows double peaks at soft X-rays, and only one peak above 20 keV. The spectrum is well described by the phenomenological models of X-ray pulsars. A cyclotron resonant scattering feature is detected with very high statistical significance at a centroid energy of $E_{\rm cyc}=90.32_{-0.28}^{+0.32}$ keV, for the reference continuum and line models, HIGHECUT and GABS respectively. Detection is very robust with respect to different continuum models. The line energy is significantly higher than what suggested from previous observations, which provided very marginal evidence for the line. This establishes a new record for the centroid energy of a fundamental cyclotron resonant scattering feature observed in accreting pulsars. We also discuss the accretion regime of the source during the Insight-HXMT observation.

preprint2020arXiv

Insight-HXMT insight into switch of the accretion mode: the case of the X-ray pulsar 4U 1901+03

We use the In data collected during the 2019 outburst from X-ray pulsar 4U 1901+03 to complement the orbital parameters reported by Fermi/GBM. Using the Insight-HXMT, we examine the correlation between the derivative of the intrinsic spin frequency and bolometric flux based on accretion torque models. It was found that the pulse profiles significantly evolve during the outburst. The existence of two types of the profile's pattern discovered in the Insight-HXMT data indicates that this source experienced transition between a super-critical and a sub-critical accretion regime during its 2019 outburst. Based on the evolution of the pulse profiles and the torque model, we derive the distance to 4U 1901+03 as 12.4+-0.2 kpc.

preprint2020arXiv

Insight-HXMT observations of a possible fast transition from jet to wind dominated state during a huge flare of GRS~1915+105

We present the analysis of the brightest flare that was recorded in the \emph{Insight}-HMXT data set, in a broad energy range (2$-$200 keV) from the microquasar GRS~1915+105 during an unusual low-luminosity state. This flare was detected by \emph{Insight}-HXMT among a series of flares during 2 June 2019 UTC 16:37:06 to 20:11:36, with a 2-200 keV luminosity of 3.4$-$7.27$\times10^{38}$ erg s$^{-1}$. Basing on the broad-band spectral analysis, we find that the flare spectrum shows different behaviors during bright and faint epochs. The spectrum of the flare can be fitted with a model dominated by a power-law component. Additional components show up in the bright epoch with a hard tail and in the faint epoch with an absorption line $\sim$ 6.78 keV. The reflection component of the latter is consistent with an inner disk radius $\sim$ 5 times larger than that of the former. These results on the giant flare during the "unusual" low-luminosity state of GRS~1915+105 may suggest that the source experiences a possible fast transition from a jet-dominated state to a wind-dominated state. We speculate that the evolving accretion disk and the large-scale magnetic field may play important roles in this peculiar huge flare.

preprint2020arXiv

Joint Analysis of Energy and RMS Spectra from MAXI J1535-571 with Insight-HXMT

A new black hole X-ray binary (BHXRB) MAXI J1535-571 was discovered by MAXI during its outburst in 2017. Using observations taken by the first Chinese X-ray satellite, the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT), we perform a joint spectral analysis (2-150 keV) in both energy and time domains. The energy spectra provide the essential input for probing the intrinsic Quasi-Periodic Oscillation (QPO) fractional rms spectra (FRS). Our results show that during the intermediate state, the energy spectra are in general consistent with those reported by Swift/XRT and NuSTAR. However, the QPO FRS become harder and the FRS residuals may suggest the presence of either an additional power-law component in the energy spectrum or a turn-over in the intrinsic QPO FRS at high energies.

preprint2020arXiv

Mapping the Galactic disk with the LAMOST and Gaia Red clump sample. VI: An evidence for the long lived non-steady warp of non-gravitational scenarios

By combining LAMOST DR4 and Gaia DR2 common red clump stars with age and proper motion, we analyze the amplitude evolution of the stellar warp independently of any assumption with a simple model. The greatest height of the warp disk increases with Galactocentric distance in different populations and it is dependent on the age: the younger stellar populations exhibit stronger warp features than the old ones, accompanied with the warp amplitude $γ({\rm age})$ decreasing with age and its first derivative $\dotγ ({\rm age})$ is different from zero. The azimuth of line of nodes $ϕ_w$ is stable at $-$5 degree without clear time evolution, which perfectly confirms some of previous works. All these self-consistent evidences support that our Galactic warp should most likely be a long lived, but non-steady structure and not a transient one, which is supporting the warp is originated from gas infall onto the disk or other hypotheses that suppose that the warp mainly affects to the gas, and consequently younger populations tracing the gas are stronger than the older ones. In other words, the Galactic warp is induced by the non-gravitational interaction over the disk models.

preprint2020arXiv

Mapping the unoccupied state dispersions in Ta$_2$NiSe$_5$ with resonant inelastic x-ray scattering

The transition metal chalcogenide Ta$_2$NiSe$_5$ undergoes a second-order phase transition at $T_c=328$ K involving a small lattice distortion. Below $T_c$, a band gap at the center of its Brillouin zone increases up to about 0.35 eV. In this work, we study the electronic structure of Ta$_2$NiSe$_5$ in its low-temperature semiconducting phase, using resonant inelastic x-ray scattering (RIXS) at the Ni $L_3$-edge. In addition to a weak fluorescence response, we observe a collection of intense Raman-like peaks that we attribute to electron-hole excitations. Using density functional theory calculations of its electronic band structure, we identify the main Raman-like peaks as interband transitions between valence and conduction bands. By performing angle-dependent RIXS measurements, we uncover the dispersion of these electron-hole excitations that allows us to extract the low-energy boundary of the electron-hole continuum. From the dispersion of the valence band measured by angle-resolved photoemission spectroscopy, we derive the effective mass of the lowest unoccupied conduction band.

preprint2020arXiv

Model comparison from LIGO-Virgo data on GW170817's binary components and consequences for the merger remnant

GW170817 is the very first observation of gravitational waves originating from the coalescence of two compact objects in the mass range of neutron stars, accompanied by electromagnetic counterparts, and offers an opportunity to directly probe the internal structure of neutron stars. We perform Bayesian model selection on a wide range of theoretical predictions for the neutron star equation of state. For the binary neutron star hypothesis, we find that we cannot rule out the majority of theoretical models considered. In addition, the gravitational-wave data alone does not rule out the possibility that one or both objects were low-mass black holes. We discuss the possible outcomes in the case of a binary neutron star merger, finding that all scenarios from prompt collapse to long-lived or even stable remnants are possible. For long-lived remnants, we place an upper limit of 1.9 kHz on the rotation rate. If a black hole was formed any time after merger and the coalescing stars were slowly rotating, then the maximum baryonic mass of non-rotating neutron stars is at most 3.05 $M_\odot$, and three equations of state considered here can be ruled out. We obtain a tighter limit of 2.67 $M_\odot$ for the case that the merger results in a hypermassive neutron star.

preprint2020arXiv

Probing the galactic halo with RR Lyrae stars I: The catalog

We present a catalog of 5,290 RR Lyrae stars (RRLs) with metallicities estimated from spectra of the LAMOST Experiment for Galactic Understanding and Exploration (LEGUE) and the Sloan Extension for Galactic Understanding and Exploration (SEGUE) surveys. Nearly 70 per cent of them (3,642 objects) also have systemic radial velocities measured. Given the pulsating nature of RRLs, metallicity estimates are based on spectra of individual exposures, by matching them with the synthetic templates. The systemic radial velocities are measured by fitting the observed velocity as a function of phase assuming an empirical pulsating velocity template curve. Various tests show that our analyses yield metallicities with a typical precision of 0.20\,dex and systemic radial velocities with uncertainties ranging from 5 to 21\,km\,s$^{-1}$ (depending on the number of radial velocity measurements available for a given star). Based on the well calibrated near-infrared $PM_{W1}Z$ or $PM_{K_{\rm s}}Z$, and $M_{V}$-[Fe/H] relations, precise distances are derived for these RRLs. Finally, we include Gaia DR2 proper motions in our catalog. The catalog should be very useful for various Galactic studies, especially of the Galactic halo.

preprint2020arXiv

Properties and astrophysical implications of the 150 Msun binary black hole merger GW190521

The gravitational-wave signal GW190521 is consistent with a binary black hole merger source at redshift 0.8 with unusually high component masses, $85^{+21}_{-14}\,M_{\odot}$ and $66^{+17}_{-18}\,M_{\odot}$, compared to previously reported events, and shows mild evidence for spin-induced orbital precession. The primary falls in the mass gap predicted by (pulsational) pair-instability supernova theory, in the approximate range $65 - 120\,M_{\odot}$. The probability that at least one of the black holes in GW190521 is in that range is 99.0%. The final mass of the merger $(142^{+28}_{-16}\,M_{\odot})$ classifies it as an intermediate-mass black hole. Under the assumption of a quasi-circular binary black hole coalescence, we detail the physical properties of GW190521's source binary and its post-merger remnant, including component masses and spin vectors. Three different waveform models, as well as direct comparison to numerical solutions of general relativity, yield consistent estimates of these properties. Tests of strong-field general relativity targeting the merger-ringdown stages of coalescence indicate consistency of the observed signal with theoretical predictions. We estimate the merger rate of similar systems to be $0.13^{+0.30}_{-0.11}\,{\rm Gpc}^{-3}\,\rm{yr}^{-1}$. We discuss the astrophysical implications of GW190521 for stellar collapse, and for the possible formation of black holes in the pair-instability mass gap through various channels: via (multiple) stellar coalescence, or via hierarchical merger of lower-mass black holes in star clusters or in active galactic nuclei. We find it to be unlikely that GW190521 is a strongly lensed signal of a lower-mass black hole binary merger. We also discuss more exotic possible sources for GW190521, including a highly eccentric black hole binary, or a primordial black hole binary.

preprint2020arXiv

Radiation Campaign of HPK Prototype LGAD sensors for the High-Granularity Timing Detector (HGTD)

We report on the results of a radiation campaign with neutrons and protons of Low Gain Avalanche Detectors (LGAD) produced by Hamamatsu (HPK) as prototypes for the High-Granularity Timing Detector (HGTD) in ATLAS. Sensors with an active thickness of 50~$μ$m were irradiated in steps of roughly 2$\times$ up to a fluence of $3\times10^{15}~\mathrm{n_{eq}cm^{-2}}$. As a function of the fluence, the collected charge and time resolution of the irradiated sensors will be reported for operation at $-30^{\circ}$.

preprint2020arXiv

Search for Eccentric Binary Black Hole Mergers with Advanced LIGO and Advanced Virgo during their First and Second Observing Runs

When formed through dynamical interactions, stellar-mass binary black holes may retain eccentric orbits ($e>0.1$ at 10 Hz) detectable by ground-based gravitational-wave detectors. Eccentricity can therefore be used to differentiate dynamically-formed binaries from isolated binary black hole mergers. Current template-based gravitational-wave searches do not use waveform models associated to eccentric orbits, rendering the search less efficient to eccentric binary systems. Here we present results of a search for binary black hole mergers that inspiral in eccentric orbits using data from the first and second observing runs (O1 and O2) of Advanced LIGO and Advanced Virgo. The search uses minimal assumptions on the morphology of the transient gravitational waveform. We show that it is sensitive to binary mergers with a detection range that is weakly dependent on eccentricity for all bound systems. Our search did not identify any new binary merger candidates. We interpret these results in light of eccentric binary formation models.

preprint2020arXiv

Search for intermediate mass black hole binaries in the first and second observing runs of the Advanced LIGO and Virgo network

Gravitational wave astronomy has been firmly established with the detection of gravitational waves from the merger of ten stellar mass binary black holes and a neutron star binary. This paper reports on the all-sky search for gravitational waves from intermediate mass black hole binaries in the first and second observing runs of the Advanced LIGO and Virgo network. The search uses three independent algorithms: two based on matched filtering of the data with waveform templates of gravitational wave signals from compact binaries, and a third, model-independent algorithm that employs no signal model for the incoming signal. No intermediate mass black hole binary event was detected in this search. Consequently, we place upper limits on the merger rate density for a family of intermediate mass black hole binaries. In particular, we choose sources with total masses $M=m_1+m_2\in[120,800]$M$_\odot$ and mass ratios $q = m_2/m_1 \in[0.1,1.0]$. For the first time, this calculation is done using numerical relativity waveforms (which include higher modes) as models of the real emitted signal. We place a most stringent upper limit of $0.20$~Gpc$^{-3}$yr$^{-1}$ (in co-moving units at the 90% confidence level) for equal-mass binaries with individual masses $m_{1,2}=100$M$_\odot$ and dimensionless spins $χ_{1,2}= 0.8$ aligned with the orbital angular momentum of the binary. This improves by a factor of $\sim 5$ that reported after Advanced LIGO's first observing run.

preprint2020arXiv

Surface Crystallization of Monoatomic Pd Metallic Glasses

Crystallization from an amorphous atomic structure is usually seen as a spontaneous process in pursuit of a lower energy state, but for alloy systems it is often hard to elucidate because of the intrinsic structural and compositional complexity. Here, by means of electron beam irradiation, we found surface-limited, and thus size-dependent crystallization in a system of monoatomic Pd metallic glass, which is ascribed to the structural differences between the surface and the interior. The equilibrium thickness of the surface crystallization is controllable, presenting a promising approach to fabricate novel nanostructures. The investigation is believed to provide a general understanding of solid amorphous-to-crystalline phase transition from the nanoscale to the bulk size.

preprint2020arXiv

Switches between accretion structures during flares in 4U 1901+03

We report on our analysis of the 2019 outburst of the X-ray accreting pulsar 4U 1901+03 observed with Insight-HXMT and NICER. Both spectra and pulse profiles evolve significantly in the decaying phase of the outburst. Dozens of flares are observed throughout the outburst. They are more frequent and brighter at the outburst peak. We find that the flares, which have a duration from tens to hundreds of seconds, are generally brighter than the persistent emission by a factor of $\sim$ 1.5. The pulse profile shape during the flares can be significantly different than that of the persistent emission. In particular, a phase shift is clearly observed in many cases. We interpret these findings as direct evidence of changes of the pulsed beam pattern, due to transitions between the sub- and super-critical accretion regimes on a short time scale. We also observe that at comparable luminosities the flares' pulse profiles are rather similar to those of the persistent emission. This indicates that the accretion on the polar cap of the neutron star is mainly determined by the luminosity, i.e., the mass accretion rate.

preprint2020arXiv

The evolution of the broadband temporal features observed in the black-hole transient MAXI J1820+070 with Insight-HXMT

We study the evolution of the temporal properties of MAXI 1820+070 during the 2018 outburst in its hard state from MJD 58190 to 58289 with Insight-HXMT in a broad energy band 1-150 keV. We find different behaviors of the hardness ratio, the fractional rms and time lag before and after MJD 58257, suggesting a transition occurred at around this point. The observed time lags between the soft photons in the 1-5 keV band and the hard photons in higher energy bands, up to 150 keV, are frequency-dependent: the time lags in the low-frequency range, 2-10 mHz, are both soft and hard lags with a timescale of dozens of seconds but without a clear trend along the outburst; the time lags in the high-frequency range, 1-10 Hz, are only hard lags with a timescale of tens of milliseconds; first increase until around MJD 58257 and decrease after this date. The high-frequency time lags are significantly correlated to the photon index derived from the fit to the quasi-simultaneous NICER spectrum in the 1-10 keV band. This result is qualitatively consistent with a model in which the high-frequency time lags are produced by Comptonization in a jet.

preprint2020arXiv

The Integrated Pulse Profiles of Fast Radio Bursts

Multi-peaked features appear on the integrated pulse profiles of fast radio burst observed below 2.5 GHz and the instantaneous spectrum of many bursts observed between 4 and 8 GHz. The mechanism of pulse or spectrum shaping has attracted little attention. Here we show that these interference-like pulse profiles are mostly the instantaneous spectra near the source regions of fast radio bursts. The corresponding instantaneous spectra are coincident to the spectrum from a single electron passing through a tapered undulator. The multi-peaked spectrum observed between 4 and 8 GHz can also be explained consistently by this type of spectrum. The spectrum is invisible unless the particles in the radiation beam are bunched. The bunching effect is probably due to the acceleration of particles in the plasma wakefield.

preprint2020arXiv

Three-Dimensional Distribution of the Interstellar Dust in the Milky Way

We present a three-dimensional (3D) extinction map of the southern sky. The map covers the SkyMapper Southern Survey (SMSS) area of $\sim$ 14,000 ${\rm deg^{2}}$ and has spatial resolutions between 6.9 and 27 arcmin. Based on the multi-band photometry of SMSS, the Two Micron All Sky Survey, the Wide-Field Infrared Survey Explorer Survey and the Gaia mission, we have estimated values of the $r$-band extinction for $\sim$ 19 million stars with the spectral energy distribution (SED) analysis. Together with the distances calculated from the Gaia data release 2 (DR2) parallaxes, we have constructed a three-dimensional extinction map of the southern sky. By combining our 3D extinction map with those from the literature, we present an all-sky 3D extinction map, and use it to explore the 3D distribution of the Galactic dust grains. We use two different models, one consisting a single disk and another of two disks, to fit the 3D distribution of the Galactic dust grains. The data is better fitted by a two-disk model, yielding smaller values of the Bayesian Information Criterion (BIC). The best fit model has scale heights of 73 and 225 pc for the "thin" and "thick" dust disks, respectively.

preprint2019arXiv

$Insight$-HXMT Science Operations

The Hard X-ray Modulation Telescope ($Insight$-HXMT) was successfully launched on June 15th, 2017. It performs broad band X-ray scan survey of the Galactic Plane to detect new black holes and other objects in active states. It also observes X-ray binaries to study their X-ray variabilities. Here we will introduce the Science Operations of $Insight$-HXMT, which is responsible for collecting and evaluating observation proposals, scheduling observations, and monitoring the working status of the payloads.

preprint2019arXiv

An arm length stabilization system for KAGRA and future gravitational-wave detectors

Modern ground-based gravitational wave (GW) detectors require a complex interferometer configuration with multiple coupled optical cavities. Since achieving the resonances of the arm cavities is the most challenging among the lock acquisition processes, the scheme called arm length stabilization (ALS) had been employed for lock acquisition of the arm cavities. We designed a new type of the ALS, which is compatible with the interferometers having long arms like the next generation GW detectors. The features of the new ALS are that the control configuration is simpler than those of previous ones and that it is not necessary to lay optical fibers for the ALS along the kilometer-long arms of the detector. Along with simulations of its noise performance, an experimental test of the new ALS was performed utilizing a single arm cavity of KAGRA. This paper presents the first results of the test where we demonstrated that lock acquisition of the arm cavity was achieved using the new ALS and residual noise was measured to be $8.2\,\mathrm{Hz}$ in units of frequency, which is smaller than the linewidth of the arm cavity and thus low enough to lock the full interferometer of KAGRA in a repeatable and reliable manner.

preprint2019arXiv

An Optically Targeted Search for Gravitational Waves emitted by Core-Collapse Supernovae during the First and Second Observing Runs of Advanced LIGO and Advanced Virgo

We present the results from a search for gravitational-wave transients associated with core-collapse supernovae observed within a source distance of approximately 20 Mpc during the first and second observing runs of Advanced LIGO and Advanced Virgo. No significant gravitational-wave candidate was detected. We report the detection efficiencies as a function of the distance for waveforms derived from multidimensional numerical simulations and phenomenological extreme emission models. For neutrino-driven explosions the distance at which we reach 50% detection efficiency is approaching 5 kpc, and for magnetorotationally-driven explosions is up to 54 kpc. However, waveforms for extreme emission models are detectable up to 28 Mpc. For the first time, the gravitational-wave data enabled us to exclude part of the parameter spaces of two extreme emission models with confidence up to 83%, limited by coincident data coverage. Besides, using ad hoc harmonic signals windowed with Gaussian envelopes we constrained the gravitational-wave energy emitted during core-collapse at the levels of $4.27\times 10^{-4}\,M_\odot c^2$ and $1.28\times 10^{-1}\,M_\odot c^2$ for emissions at 235 Hz and 1304 Hz respectively. These constraints are two orders of magnitude more stringent than previously derived in the corresponding analysis using initial LIGO, initial Virgo and GEO 600 data.

preprint2019arXiv

Mapping the Galactic disk with the LAMOST and Gaia Red clump sample. II. 3D asymmetrical kinematics of mono-age populations in the disk between 6$-14$ \,kpc

We perform analysis of the three-dimensional kinematics of Milky Way disk stars in mono-age populations. We focus on stars between Galactocentric distances of $R=6$ and 14 \,kpc, selected from the combined LAMOST DR4 red clump giant stars and Gaia DR2 proper motion catalogue. We confirm the 3D asymmetrical motions of recent works and provide time tagging of the Galactic outer disk asymmetrical motions near the anticenter direction out to Galactocentric distances of 14\,kpc. Radial Galactocentric motions reach values up to 10 km s$^{-1}$, depending on the age of the population, and present a north-south asymmetry in the region corresponding to density and velocity substructures that were sensitive to the perturbations in the early 6 \,Gyr. After that time, the disk stars in this asymmetrical structure have become kinematically hotter, and are thus not sensitive to perturbations, and we find the structure is a relatively younger population. With quantitative analysis, we find stars both above and below the plane at $R\gtrsim 9$ kpc that exhibit bending mode motions of which the sensitive duration is around 8 \,Gyr. We speculate that the in-plane asymmetries might not be mainly caused by a fast rotating bar, intrinsically elliptical outer disk, secular expansion of the disk, or streams. Spiral arm dynamics, out-of-equilibrium models, minor mergers or others are important contributors. Vertical motions might be dominated by bending and breathing modes induced by complicated inner or external perturbers. It is likely that many of these mechanisms are coupled together.

preprint2019arXiv

Search for sub-solar mass ultracompact binaries in Advanced LIGO's second observing run

We present an Advanced LIGO and Advanced Virgo search for sub-solar mass ultracompact objects in data obtained during Advanced LIGO's second observing run. In contrast to a previous search of Advanced LIGO data from the first observing run, this search includes the effects of component spin on the gravitational waveform. We identify no viable gravitational wave candidates consistent with sub-solar mass ultracompact binaries with at least one component between 0.2 - 1.0 solar masses. We use the null result to constrain the binary merger rate of (0.2 solar mass, 0.2 solar mass) binaries to be less than 3.7 x 10^5 Gpc^-3 yr^-1 and the binary merger rate of (1.0 solar mass, 1.0 solar mass) binaries to be less than 5.2 x 10^3 Gpc^-3 yr^-1. Sub-solar mass ultracompact objects are not expected to form via known stellar evolution channels, though it has been suggested that primordial density fluctuations or particle dark matter with cooling mechanisms and/or nuclear interactions could form black holes with sub-solar masses. Assuming a particular primordial black hole formation model, we constrain a population of merging 0.2 solar mass black holes to account for less than 16% of the dark matter density and a population of merging 1.0 solar mass black holes to account for less than 2% of the dark matter density. We discuss how constraints on the merger rate and dark matter fraction may be extended to arbitrary black hole population models that predict sub-solar mass binaries.

preprint2019arXiv

Streaming controlled by meniscus shape

Surface waves called meniscus waves often appear in the systems that are close to the capillary length scale. Since the meniscus shape determines the form of the meniscus waves, the resulting streaming circulation has a structure distinct from that caused by other capillary-gravity waves recently reported in the literature. In the present study, we produce symmetric and antisymmetric meniscus shapes by controlling boundary wettability and excite meniscus waves by oscillating the meniscus vertically. The symmetric and antisymmetric configurations produce different surface capillary-gravity wave modes and streaming flow structures. The energy density of the streaming circulation increases at the rate of the fourth power of the forcing amplitude in both configurations. The flow symmetry of streaming circulation is retained under the symmetric meniscus, while it is lost under the antisymmetric meniscus. In our experiments, the streaming circulation primarily originates from the Stokes boundary layer beneath the meniscus and can be successfully explained using the existing streaming theory.

preprint2018arXiv

A variational principle for fluid sloshing with vorticity, dynamically coupled to vessel motion

A variational principle is derived for two-dimensional incompressible rotational fluid flow with a free surface in a moving vessel when both the vessel and fluid motion are to be determined. The fluid is represented by a stream function and the vessel motion is represented by a path in the planar Euclidean group. Novelties in the formulation include how the pressure boundary condition is treated, the introduction of a stream function into the Euler-Poincaré variations, the derivation of free surface variations, and how the equations for the vessel path in the Euclidean group, coupled to the fluid motion, are generated automatically.

preprint2018arXiv

Charge trapping and super-Poissonian noise centers in a cuprate high-temperature superconductor

The electronic properties of cuprate high temperature superconductors in their normal state are very two-dimensional: while transport in the ab plane is perfectly metallic, it is insulating along the c-axis, with ratios between the two exceeding 10^4. This anisotropy has been identified as one of the mysteries of the cuprates early on, and while widely different proposals exist for its microscopic origin, little is known empirically on the microscopic scale. Here, we elucidate the properties of the insulating layers with a newly developed scanning noise spectroscopy technique that can spatially map not only the current but also the current fluctuations in time. We discover atomic-scale noise centers that exhibit MHz current fluctuations 40 times the expectation from Poissonian noise, more than what has been observed in mesoscopic systems. Such behaviour can only happen in highly polarizable insulators and represents strong evidence for trapping of charge in the charge reservoir layers. Our measurements suggest a picture of metallic layers separated by polarizable insulators within a three-dimensional superconducting state.

preprint2016arXiv

A universal description of pseudorapidity distributions in both nucleus-nucleus and p-p collisions at now available energy regions

Investigations have shown that the collective motion not only appears in nucleus-nucleus but also in p-p collisions. The best tool for depicting such collective motion is relativistic hydrodynamics. In this paper, the collective motion is assumed obeying the hydro model which integrates the features of Landau and Hwa-Bjorken theory and is one of a very few analytically solvable models. The fluid is then supposed freezing out into charged particles from a time-like hypersurface with a fixed time of t_FO . The researches of present paper show that this part of charged particles together with leading particles, which, by conventional definition, carry on the quantum numbers of colliding nucleons and take away the most part of incident energy, can give a proper universal description to the pseudorapidity distributions of charged particles measured in both nucleus-nucleus and p-p collisions at now available energy regions.

preprint2016arXiv

Andreev reflection in s-type superconductor proximized 3D topological insulator

We investigate transport and shot noise in lateral N-TI-S contacts, where N is a normal metal, TI is a Bi-based three dimensional topological insulator (3D TI), and S is an s-type superconductor. In normal state, the devices are in the elastic diffusive transport regime, as demonstrated by a nearly universal value of the shot noise Fano factor $F_{\rm N}\approx1/3$ in magnetic field and in reference normal contact. In the absence of magnetic field, we identify the Andreev reflection (AR) regime, which gives rise to the effective charge doubling in shot noise measurements. Surprisingly, the Fano factor $F_{\rm AR}\approx0.22\pm0.02$ is considerably reduced in the AR regime compared to $F_{\rm N}$, in contrast to previous AR experiments in normal metals and semiconductors. We suggest that this effect is related to a finite thermal conduction of the proximized, superconducting TI owing to a residual density of states at low energies.

preprint2016arXiv

Epitaxial growth and thermodynamic stability of SrIrO3/SrTiO3 heterostructures

Obtaining high-quality thin films of 5d transition metal oxides is essential to explore the exotic semimetallic and topological phases predicted to arise from the combination of strong electron correlations and spin-orbit coupling. Here, we show that the transport properties of SrIrO3 thin films, grown by pulsed laser deposition, can be optimized by considering the effect of laser-induced modification of the SrIrO3 target surface. We further demonstrate that bare SrIrO3 thin films are subject to degradation in air and are highly sensitive to lithographic processing. A crystalline SrTiO3 cap layer deposited in-situ is effective in preserving the film quality, allowing us to measure metallic transport behavior in films with thicknesses down to 4 unit cells. In addition, the SrTiO3 encapsulation enables the fabrication of devices such as Hall bars without altering the film properties, allowing precise (magneto)transport measurements on micro- and nanoscale devices.

preprint2016arXiv

GaAs/AlGaAs Nanowire Photodetector

We demonstrate an efficient core-shell GaAs/AlGaAs nanowire photodetector operating at room temperature. The design of this nanoscale detector is based on a type-I heterostructure combined with a metal-semiconductor-metal (MSM) radial architecture, in which built-in electric fields at the semiconductor heterointerface and at the metal/semiconductor Schottky contact promote photogenerated charge separation, enhancing photosensitivity. The spectral photoconductive response shows that the nanowire supports resonant optical modes in the near-infrared region, which lead to large photocurrent density in agreement with the predictions of electromagnetic and transport computational models. The single nanowire photodetector shows remarkable peak photoresponsivity of 0.57 A/W, comparable to large-area planar GaAs photodetectors on the market, and a high detectivity of 7.2 10^10 cm\sqrt{Hz}/W at λ=855 nm. This is promising for the design of a new generation of highly sensitive single nanowire photodetectors by controlling optical mode confinement, bandgap, density of states, and electrode engineering.

preprint2016arXiv

Improved measurement of the absolute branching fraction of $D^{+}\rightarrow \bar K^0 μ^{+}ν_μ$

By analyzing 2.93 fb$^{-1}$ of data collected at $\sqrt s=3.773$ GeV with the BESIII detector, we measure the absolute branching fraction ${\mathcal B}(D^{+}\rightarrow\bar K^0μ^{+}ν_μ)=(8.72 \pm 0.07_{\rm stat.} \pm 0.18_{\rm sys.})\%$, which is consistent with previous measurements within uncertainties but with significantly improved precision. Combining the Particle Data Group values of ${\mathcal B}(D^0\to K^-μ^+ν_μ)$, ${\mathcal B}(D^{+}\rightarrow\bar K^0 e^{+}ν_{e})$, and the lifetimes of the $D^0$ and $D^+$ mesons with the value of ${\mathcal B}(D^{+}\rightarrow\bar K^0 μ^{+}ν_μ)$ measured in this work, we determine the following ratios of partial widths: $Γ(D^0\to K^-μ^+ν_μ)/Γ(D^{+}\rightarrow\bar K^0μ^{+}ν_μ)=0.963\pm0.044$ and $Γ(D^{+}\rightarrow\bar K^0 μ^{+}ν_μ)/Γ(D^{+}\rightarrow\bar K^0 e^{+}ν_{e})=0.988\pm0.033$.

preprint2016arXiv

Measurement of azimuthal asymmetries in inclusive charged dipion production in $e^+e^-$ annihilations at $\sqrt{s}$ = 3.65 GeV

We present a measurement of the azimuthal asymmetries of two charged pions in the inclusive process $e^+e^-\rightarrow ππX$ based on a data set of 62 $\rm{pb}^{-1}$ at the center-of-mass energy $\sqrt{s}=3.65$ GeV collected with the BESIII detector. These asymmetries can be attributed to the Collins fragmentation function. We observe a nonzero asymmetry, which increases with increasing pion momentum. As our energy scale is close to that of the existing semi-inclusive deep inelastic scattering experimental data, the measured asymmetries are important inputs for the global analysis of extracting the quark transversity distribution inside the nucleon and are valuable to explore the energy evolution of the spin-dependent fragmentation function.

preprint2016arXiv

Measurement of the $D_s^+ \to \ell^+ν_\ell$ branching fractions and the decay constant $f_{D_s^+}$

Using 482~pb$^{-1}$ of $e^+e^-$ collision data collected at a center-of-mass energy of $\sqrt{s} = 4.009$ GeV with the BESIII detector, we measure the branching fractions of the decays $D_s^+\toμ^+ν_μ$ and $D_s^+\toτ^+ν_τ$. By constraining the ratio of decay rates of $D_s^+$ to $τ^+ν_τ$ and to $μ^+ν_μ$ to the Standard Model prediction, the branching fractions are determined to be $\mathcal{B}(D_s^+ \to μ^+ν_μ) = (0.495 \pm 0.067 \pm 0.026)\%$ and $\mathcal{B}(D_s^+ \to τ^+ν_τ) = (4.83 \pm 0.65 \pm 0.26)\%$. Using these branching fractions, we obtain a value for the decay constant $f_{D_s^+}$ of $(241.0 \pm 16.3 \pm 6.5)~\text{MeV}$, where the first error is statistical and the second systematic.

preprint2016arXiv

Measurement of the absolute branching fraction of $D^{+}\rightarrow\bar K^0 e^{+}ν_{e}$ via $\bar K^0\toπ^0π^0$

By analyzing 2.93 fb$^{-1}$ data collected at the center-of-mass energy $\sqrt s=3.773$ GeV with the BESIII detector, we measure the absolute branching fraction of the semileptonic decay $D^+\rightarrow\bar K^0 e^{+}ν_{e}$ to be ${\mathcal B}(D^{+}\rightarrow\bar K^0 e^{+}ν_{e})=(8.59 \pm 0.14 \pm 0.21)\%$ using $\bar K^0\to K^0_S\to π^0π^0$, where the first uncertainty is statistical and the second systematic. Our result is consistent with previous measurements within uncertainties.

preprint2016arXiv

Measurement of the leptonic decay width of $J/ψ$ using initial state radiation

Using a data set of 2.93 fb$^{-1}$ taken at a center-of-mass energy of $\sqrt{s}$ = 3.773 GeV with the BESIII detector at the BEPCII collider, we measure the process $e^+e^-\rightarrow J/ψγ\rightarrow μ^+μ^-γ$ and determine the product of the branching fraction and the electronic width $\mathcal B_{μμ}\cdot Γ_{ee} = (333.4 \pm 2.5_{\rm stat} \pm 4.4_{\rm sys})$~eV. Using the earlier-published BESIII result for $\mathcal B_{μμ}$ = (5.973 $\pm$ 0.007$_{\rm stat}$ $\pm$ 0.037$_{\rm sys}$)\%, we derive the $J/ψ$ electronic width $Γ_{ee}$~= (5.58 $\pm$ 0.05$_{\rm stat}$ $\pm$ 0.08$_{\rm sys}$) keV.

preprint2016arXiv

Measurements of absolute hadronic branching fractions of $Λ_{c}^{+}$ baryon

Using $567\rm{pb}^{-1}$ of $e^+e^-$ collisions recorded at $\sqrt{s}=4.599\rm{GeV}$ with the BESIII detector, we report first measurements of absolute hadronic branching fractions of Cabibbo-favored decays of the $Λ_{c}^{+}$ baryon with a double-tag technique. A global least-square fitter is utilized to improve the measured precision. Among the measurements for twelve $Λ_{c}^{+}$ decay modes, the branching fraction for $Λ_{c}^{+} \rightarrow pK^-π^+$ is determined to be $(5.84\pm0.27\pm0.23)\%$, where the first uncertainty is statistical and the second is systematic. In addition, the measurements of the branching fractions of the other eleven Cabbibo-favored hadronic decay modes are significantly improved.

preprint2016arXiv

Numerical Study of a Particle Method for Gradient Flows

We study the numerical behaviour of a particle method for gradient flows involving linear and nonlinear diffusion. This method relies on the discretisation of the energy via non-overlapping balls centred at the particles. The resulting scheme preserves the gradient flow structure at the particle level, and enables us to obtain a gradient descent formulation after time discretisation. We give several simulations to illustrate the validity of this method, as well as a detailed study of one-dimensional aggregation-diffusion equations.

preprint2016arXiv

Observation of $e^{+}e^{-} \to η^{\prime} J/ψ$ at center-of-mass energies between 4.189 and 4.600 GeV

The process $e^{+}e^{-}\to η^{\prime} J/ψ$ is observed for the first time with a statistical significance of $8.6σ$ at center-of-mass energy $\sqrt{s} = 4.226$ GeV and $7.3σ$ at $\sqrt{s} = 4.258$ GeV using data samples collected with the BESIII detector. The Born cross sections are measured to be $(3.7 \pm 0.7 \pm 0.3)$ and $(3.9 \pm 0.8 \pm 0.3)$ pb at $\sqrt{s} = 4.226$ and $4.258$ GeV, respectively, where the first errors are statistical and the second systematic. Upper limits at the 90% confidence level of the Born cross sections are also reported at other 12 energy points.

preprint2016arXiv

Observation of $e^+e^- \rightarrow ωχ_{c1,2}$ near $\sqrt{s}$ = 4.42 and 4.6 GeV

Based on data samples collected with the BESIII detector operating at the BEPCII storage ring at center-of-mass energies $\sqrt{s} >$ 4.4 GeV, the processes $e^+e^- \rightarrow ωχ_{c1,2}$ are observed for the first time. With an integrated luminosity of $1074 pb^{-1}$ near $\sqrt{s} =$ 4.42 GeV, a significant $ωχ_{c2}$ signal is found, and the cross section is measured to be $(20.9 \pm 3.2 \pm 2.5)\pb$. With $567 pb^{-1}$ near $\sqrt{s} =$ 4.6 GeV, a clear $ωχ_{c1}$ signal is seen, and the cross section is measured to be $(9.5 \pm 2.1 \pm 1.3) \pb$, while evidence is found for an $ωχ_{c2}$ signal. The first errors are statistical and the second are systematic. Due to low luminosity or low cross section at other energies, no significant signals are observed. In the $ωχ_{c2}$ cross section, an enhancement is seen around $\sqrt{s} =$ 4.42 GeV. Fitting the cross section with a coherent sum of the $ψ(4415)$ Breit-Wigner function and a phase space term, the branching fraction $\mathcal{B}(ψ(4415)\toωχ_{c2})$ is obtained to be of the order of $10^{-3}$.

preprint2016arXiv

Observation of $h_{c}$ radiative decay $h_{c} \rightarrow γη'$ and evidence for $h_{c} \rightarrow γη$

A search for radiative decays of the $P$-wave spin singlet charmonium resonance $h_c$ is performed based on $4.48 \times 10^{8}$ $ψ'$ events collected with the BESIII detector operating at the BEPCII storage ring. Events of the reaction channels $h_{c} \rightarrow γη'$ and $γη$ are observed with a statistical significance of $8.4 σ$ and $4.0 σ$, respectively, for the first time. The branching fractions of $h_{c} \rightarrow γη'$ and $h_{c} \rightarrow γη$ are measured to be $\mathcal{B}(h_{c} \rightarrow γη')=(1.52 \pm 0.27 \pm 0.29)\times10^{-3}$ and $\mathcal{B}(h_{c} \rightarrow γη)=(4.7 \pm 1.5 \pm 1.4)\times10^{-4}$, respectively, where the first errors are statistical and the second are systematic uncertainties.

preprint2016arXiv

Observation of an anomalous line shape of the $η^{\prime}π^{+}π^{-}$ mass spectrum near the $p\bar{p}$ mass threshold in $J/ψ\rightarrowγη^{\prime}π^{+}π^{-}$

Using $1.09\times10^{9}$ $J/ψ$ events collected by the BESIII experiment in 2012, we study the $J/ψ\rightarrowγη^{\prime}π^{+}π^{-}$ process and observe a significant abrupt change in the slope of the $η^{\prime}π^{+}π^{-}$ invariant mass distribution at the proton-antiproton ($p\bar{p}$) mass threshold. We use two models to characterize the $η^{\prime}π^{+}π^{-}$ line shape around $1.85~\text{GeV}/c^{2}$: one which explicitly incorporates the opening of a decay threshold in the mass spectrum (Flatté formula), and another which is the coherent sum of two resonant amplitudes. Both fits show almost equally good agreement with data, and suggest the existence of either a broad state around $1.85~\text{GeV}/c^{2}$ with strong couplings to $p\bar{p}$ final states or a narrow state just below the $p\bar{p}$ mass threshold. Although we cannot distinguish between the fits, either one supports the existence of a $p\bar{p}$ molecule-like state or bound state with greater than $7σ$ significance.

preprint2016arXiv

Observation of pseudoscalar and tensor resonances in $J/ψ\rightarrow γϕϕ$

Based on a sample of $(1310.6 \pm 10.5) \times 10^{6}$ $J/ψ$ events collected with the BESIII detector operating at the BEPCII storage ring, a partial wave analysis of the decay $J/ψ\rightarrow γϕϕ$ is performed in order to study the intermediate states. Results of the partial wave analysis show that the structures are predominantly $0^{-+}$ states. The existence of the $η(2225)$ is confirmed, and its resonance parameters are measured. Two additional pseudoscalar states, the $η(2100)$ with a mass of $2050_{-24}^{+30}{}_{-26}^{+75}$~MeV/$c^{2}$ and a width of $250_{-30}^{+36}{}_{-164}^{+181}$~MeV/$c^{2}$ and the $X(2500)$ with a mass of $2470_{-19}^{+15}{}_{-23}^{+101}$~MeV/$c^{2}$ and a width of $230_{-35}^{+64}{}_{-33}^{+56}$~MeV/$c^{2}$, are observed. In addition to these three pseudoscalar states, the scalar state $f_0(2100)$, and three tensor states, the $f_2(2010)$, $f_2(2300)$ and $f_2(2340)$, are observed in the process $J/ψ\rightarrow γϕϕ$. The product branching fractions ${\cal B}$($J/ψ\rightarrow γX)\times{\cal B}(X\rightarrow ϕϕ)$ are reported.

preprint2016arXiv

Quantifying discordance in the 2015 Planck CMB spectrum

We examine the internal consistency of the Planck 2015 cosmic microwave background (CMB) temperature anisotropy power spectrum. We show that tension exists between cosmological constant cold dark matter (LCDM) model parameters inferred from multipoles l<1000 (roughly those accessible to Wilkinson Microwave Anisotropy Probe), and from l>=1000, particularly the CDM density, Omega_ch^2, which is discrepant at 2.5 sigma for a Planck-motivated prior on the optical depth, tau=0.07+/-0.02. We find some parameter tensions to be larger than previously reported because of inaccuracy in the code used by the Planck Collaboration to generate model spectra. The Planck l>=1000 constraints are also in tension with low-redshift data sets, including Planck's own measurement of the CMB lensing power spectrum (2.4 sigma), and the most precise baryon acoustic oscillation (BAO) scale determination (2.5 sigma). The Hubble constant predicted by Planck from l>=1000, H_0=64.1+/-1.7 km/s/Mpc, disagrees with the most precise local distance ladder measurement of 73.0+/-2.4 km/s/Mpc at the 3.0 sigma level, while the Planck value from l<1000, 69.7+/-1.7 km/s/Mpc, is consistent within 1 sigma. A discrepancy between the Planck and South Pole Telescope (SPT) high-multipole CMB spectra disfavors interpreting these tensions as evidence for new physics. We conclude that the parameters from the Planck high-multipole spectrum probably differ from the underlying values due to either an unlikely statistical fluctuation or unaccounted-for systematics persisting in the Planck data.

preprint2016arXiv

Quenched Magnon excitations by oxygen sublattice reconstruction in (SrCuO$_2$)$_n$/(SrTiO$_3$)$_2$ superlattices

The recently discovered structural reconstruction in the cuprate superlattice (SrCuO$_2$)$_n$/(SrTiO$_3$)$_2$ has been investigated across the critical value of $n=5$ using resonant inelastic x-ray scattering (RIXS). We find that at the critical value of $n$, the cuprate layer remains largely in the bulk-like two-dimensional structure with a minority of Cu plaquettes being reconstructed. The partial reconstruction leads to quenching of the magnons starting at the $Γ$-point due to the minority plaquettes acting as scattering points. Although comparable in relative abundance, the doped charge impurities in electron-doped cuprate superconductors do not show this quenching of magnetic excitations.

preprint2016arXiv

Search for the weak decay $η^{\prime}\to K^{\pm}π^{\mp}$ and precise measurement of the branching fraction ${\cal B}(J/ψ\toϕη^{\prime})$

We present the first search for the rare decay of $η^\prime$ into $K^{\pm}π^{\mp}$ in $J/ψ\toϕη^{\prime}$, using a sample of $1.3\times10^{9}$ $J/ψ$ events collected with the BESIII detector. No significant signal is observed, and the upper limit at the 90\% confidence level for the ratio $\frac{{\cal B}(η^{\prime}\to K^{\pm}π^{\mp})}{{\cal B}(η^{\prime}\toγπ^{+}π^{-})}$ is determined to be $1.3\times10^{-4}$. In addition, we report the measurement of the branching fraction of $J/ψ\toϕη^{\prime}$ to be $(5.10\pm0.03(\text{stat.})\pm0.32(\text{syst.}))\times10^{-4}$, which agrees with previous results from BESII.

preprint2016arXiv

Study of $D^{+} \to K^{-} π^+ e^+ ν_e$

We present an analysis of the decay $D^{+} \to K^{-} π^+ e^+ ν_e$ based on data collected by the BESIII experiment at the $ψ(3770)$ resonance. Using a nearly background-free sample of 18262 events, we measure the branching fraction $\mathcal{B}(D^{+} \to K^{-} π^+ e^+ ν_e) = (3.71 \pm 0.03 \pm 0.08)\%$. For $0.8<m_{Kπ}<1.0$ GeV/$c^{2}$ the partial branching fraction is $\mathcal{B}(D^{+} \to K^{-} π^+ e^+ ν_e)_{[0.8,1]} = (3.33 \pm 0.03 \pm 0.07)\%$. A partial wave analysis shows that the dominant $\bar K^{*}(892)^{0}$ component is accompanied by an \emph{S}-wave contribution accounting for $(6.05\pm0.22\pm0.18)\%$ of the total rate and that other components are negligible. The parameters of the $\bar K^{*}(892)^{0}$ resonance and of the form factors based on the spectroscopic pole dominance predictions are also measured. We also present a measurement of the $\bar K^{*}(892)^{0}$ helicity basis form factors in a model-independent way.

preprint2016arXiv

Study of $J/ψ\to p\bar{p}ϕ$ at BESIII

Using a data sample of $1.31 \times 10^{9}$ $J/ψ$ events accumulated with the BESIII detector, the decay $J/ψ\to p\bar{p}ϕ$ is studied via two decay modes, $ϕ\to K^{0}_{S}K^{0}_{L}$ and $ϕ\to K^{+}K^{-}$. The branching fraction of $J/ψ\to p\bar{p}ϕ$ is measured to be $\mathcal{B}(J/ψ\to p\bar{p}ϕ)=[5.23\pm0.06(\mbox{stat})\pm0.33(\mbox{syst})]\times10^{-5}$, which agrees well with a previously published measurement, but with a significantly improved precision. No evident enhancement near the $p\bar{p}$ mass threshold, denoted as $X(p\bar{p})$, is observed, and the upper limit on the branching fraction of $J/ψ\to X(p\bar{p})ϕ\to p\bar{p}ϕ$ is determined to be $\mathcal{B}(J/ψ\to X(p\bar{p})ϕ\to p\bar{p}ϕ)<2.1\times10^{-7}$ at the 90\% confidence level.

preprint2016arXiv

Study of $ψ$ decays to the $Ξ^{-}\barΞ^{+}$ and $Σ(1385)^{\mp}\barΣ(1385)^{\pm}$ final states

We study the decays of the charmonium resonances $J/ψ$ and $ψ(3686)$ to the final states $Ξ^{-}\barΞ^{+}$, $Σ(1385)^{\mp}\barΣ(1385)^{\pm}$ based on a single baryon tag method using data samples of $(223.7 \pm 1.4) \times 10^{6}$ $J/ψ$ and $(106.4 \pm 0.9) \times 10^{6}$ $ψ(3686)$ events collected with the BESIII detector at the BEPCII collider. The decay $ψ(3686)\rightarrowΣ(1385)^{\mp}\barΣ(1385)^{\pm}$ is observed for the first time, and the measurements of the other processes, including the branching fractions and angular distributions, are in good agreement with and much more precise than the previously published results. Additionally, the ratios $\frac{\cal{B}(ψ(3686)\rightarrowΞ^{-}\barΞ^{+})}{\cal{B}(J/ψ\rightarrowΞ^{-}\barΞ^{+})}$, $\frac{\cal{B}(ψ(3686)\rightarrowΣ(1385)^{-}\barΣ(1385)^{+})}{\cal{B}(J/ψ\rightarrowΣ(1385)^{-}\barΣ(1385)^{+})}$ and $\frac{\cal{B}(ψ(3686)\rightarrowΣ(1385)^{+}\barΣ(1385)^{-})}{\cal{B}(J/ψ\rightarrowΣ(1385)^{+}\barΣ(1385)^{-})}$ are determined.

preprint2015arXiv

An amplitude analysis of the $π^{0}π^{0}$ system produced in radiative $J/ψ$ decays

An amplitude analysis of the $π^{0}π^{0}$ system produced in radiative $J/ψ$ decays is presented. In particular, a piecewise function that describes the dynamics of the $π^{0}π^{0}$ system is determined as a function of $M_{π^{0}π^{0}}$ from an analysis of the $(1.311\pm0.011)\times10^{9}$ $J/ψ$ decays collected by the BESIII detector. The goal of this analysis is to provide a description of the scalar and tensor components of the $π^0π^0$ system while making minimal assumptions about the properties or number of poles in the amplitude. Such a model-independent description allows one to integrate these results with other related results from complementary reactions in the development of phenomenological models, which can then be used to directly fit experimental data to obtain parameters of interest. The branching fraction of $J/ψ\to γπ^{0}π^{0}$ is determined to be $(1.15\pm0.05)\times10^{-3}$, where the uncertainty is systematic only and the statistical uncertainty is negligible.

preprint2015arXiv

Computational design of high performance hybrid perovskite on silicon tandem solar cells

In this study, the optoelectronic properties of a monolithically integrated series-connected tandem solar cell are simulated. Following the large success of hybrid organic-inorganic perovskites, which have recently demonstrated large efficiencies with low production costs, we examine the possibility of using the same perovskites as absorbers in a tandem solar cell. The cell consists in a methylammonium mixed bromide-iodide lead perovskite, CH3NH3PbI3(1-x)Br3x (0 < x < 1), top sub-cell and a single-crystalline silicon bottom sub-cell. A Si-based tunnel junction connects the two sub-cells. Numerical simulations are based on a one-dimensional numerical drift-diffusion model. It is shown that a top cell absorbing material with 20% of bromide and a thickness in the 300-400 nm range affords current matching with the silicon bottom cell. Good interconnection between single cells is ensured by standard n and p doping of the silicon at 5.10^19cm-3 in the tunnel junction. A maximum efficiency of 27% is predicted for the tandem cell, exceeding the efficiencies of stand-alone silicon (17.3%) and perovskite cells (17.9%) taken for our simulations, and more importantly, that of the record crystalline Si cells.

preprint2015arXiv

Confirmation of a charged charmoniumlike state $Z_c(3885)^{\mp}$ in $e^+e^-\toπ^{\pm}(D\bar{D}^*)^\mp$ with double $D$ tag

We present a study of the process $e^+e^-\toπ^{\pm}(D\bar{D}^*)^{\mp}$ using data samples of 1092~pb$^{-1}$ at $\sqrt{s}=4.23$~GeV and 826~pb$^{-1}$ at $\sqrt{s}=4.26$~GeV collected with the BESIII detector at the BEPCII storage ring. With full reconstruction of the $D$ meson pair and the bachelor $π^{\pm}$ in the final state, we confirm the existence of the charged structure $Z_c(3885)^{\mp}$ in the $(D\bar{D}^*)^{\mp}$ system in the two isospin processes $e^+e^-\toπ^+D^0D^{*-}$ and $e^+e^-\toπ^+D^-D^{*0}$. By performing a simultaneous fit, the statistical significance of $Zc(3885)^{\mp}$ signal is determined to be greater than 10$σ$, and its pole mass and width are measured to be $M_{\rm{pole}}$=(3881.7$\pm$1.6(stat.)$\pm$1.6(syst.))~MeV/$c^2$ and $Γ_{\rm{pole}}$=(26.6$\pm$2.0(stat.)$\pm$2.1(syst.))~MeV, respectively. The Born cross section times the $(D\bar{D}^*)^{\mp}$ branching fraction ($σ(e^+e^-\toπ^{\pm}Z_{c}(3885)^{\mp}) \times Br(Z_{c}(3885)^{\mp}\to(D\bar{D}^*)^{\mp})$) is measured to be $(141.6\pm7.9(\text{stat.})\pm12.3(\text{syst.}))~\text{pb}$ at $\sqrt{s}=4.23$~GeV and $(108.4\pm6.9(\text{stat.})\pm8.8(\text{syst.}))~\text{pb}$ at $\sqrt{s}=4.26$~GeV. The polar angular distribution of the $π^{\pm}$-$Z_c(3885)^{\mp}$ system is consistent with the expectation of a quantum number assignment of $J^P=1^+$ for $Z_c(3885)^{\mp}$.

preprint2015arXiv

DA white dwarfs from the LSS-GAC survey DR1: the preliminary luminosity and mass functions and formation rate

Modern large-scale surveys have allowed the identification of large numbers of white dwarfs. However, these surveys are subject to complicated target selection algorithms, which make it almost impossible to quantify to what extent the observational biases affect the observed populations. The LAMOST (Large Sky Area Multi-Object Fiber Spectroscopic Telescope) Spectroscopic Survey of the Galactic anti-center (LSS-GAC) follows a well-defined set of criteria for selecting targets for observations. This advantage over previous surveys has been fully exploited here to identify a small yet well-characterised magnitude-limited sample of hydrogen-rich (DA) white dwarfs. We derive preliminary LSS-GAC DA white dwarf luminosity and mass functions. The space density and average formation rate of DA white dwarfs we derive are 0.83+/-0.16 x 10^{-3} pc^{-3} and 5.42 +/- 0.08 x 10^{-13} pc^{-3} yr^{-1}, respectively. Additionally, using an existing Monte Carlo population synthesis code we simulate the population of single DA white dwarfs in the Galactic anti-center, under various assumptions. The synthetic populations are passed through the LSS-GAC selection criteria, taking into account all possible observational biases. This allows us to perform a meaningful comparison of the observed and simulated distributions. We find that the LSS-GAC set of criteria is highly efficient in selecting white dwarfs for spectroscopic observations (80-85 per cent) and that, overall, our simulations reproduce well the observed luminosity function. However, they fail at reproducing an excess of massive white dwarfs present in the observed mass function. A plausible explanation for this is that a sizable fraction of massive white dwarfs in the Galaxy are the product of white dwarf-white dwarf mergers.

preprint2015arXiv

Dust-to-gas ratio, $X_{\rm CO}$ factor and CO-dark gas in the Galactic anticentre: an observational study

We investigate the correlation between extinction and H~{\sc i} and CO emission at intermediate and high Galactic latitudes ($|b|>10\degr$) within the footprint of the Xuyi Schmidt Telescope Photometric Survey of the Galactic anticentre (XSTPS-GAC) on small and large scales. In Paper I (Chen et al. 2014), we present a three-dimensional dust extinction map within the footprint of XSTPS-GAC, covering a sky area of over 6,000\,deg$^2$ at a spatial angular resolution of 6\,arcmin. In the current work, the map is combined with data from gas tracers, including H~{\sc i} data from the Galactic Arecibo L-band Feed Array H~{\sc i} survey and CO data from the Planck mission, to constrain the values of dust-to-gas ratio $DGR=A_V/N({\rm H})$ and CO-to-$\rm H_2$ conversion factor $X_{\rm CO}=N({\rm H_2})/W_{\rm CO}$ for the entire GAC footprint excluding the Galactic plane, as well as for selected star-forming regions (such as the Orion, Taurus and Perseus clouds) and a region of diffuse gas in the northern Galactic hemisphere. For the whole GAC footprint, we find $DGR=(4.15\pm0.01) \times 10^{-22}$\,$\rm mag\,cm^{2}$ and $X_{\rm CO}=(1.72 \pm 0.03) \times 10^{20}$\,$\rm cm^{-2}\,(K\,km\,s^{-1})^{-1}$. We have also investigated the distribution of "CO-dark" gas (DG) within the footprint of GAC and found a linear correlation between the DG column density and the $V$-band extinction: $N({\rm DG}) \simeq 2.2 \times 10^{21} (A_V - A^{c}_{V})\,\rm cm^{-2}$. The mass fraction of DG is found to be $f_{\rm DG}\sim 0.55$ toward the Galactic anticentre, which is respectively about 23 and 124 per cent of the atomic and CO-traced molecular gas in the same region. This result is consistent with the theoretical work of Papadopoulos et al. but much larger than that expected in the $\rm H_2$ cloud models by Wolfire et al.

preprint2015arXiv

Electron correlation effects and scattering rates in Fe$_{1+y}$Te$_{1-x}$Se$_x$ superconductor

Using angle-resolved photoemission spectroscopy we have studied the low-energy electronic structure and the Fermi surface topology of Fe$_{1+y}$Te$_{1-x}$Se$_x$ superconductors. Similar to the known iron pnictides we observe hole pockets at the center and electron pockets at the corner of the Brillouin zone (BZ). However, on a finer level, the electronic structure around the $Γ$- and $Z$-points in $k$-space is substantially different from other iron pnictides, in that we observe two hole pockets at the $Γ$-point, and more interestingly only one hole pocket is seen at the $Z$-point, whereas in $1111$-, $111$-, and $122$-type compounds, three hole pockets could be readily found at the zone center. Another major difference noted in the Fe$_{1+y}$Te$_{1-x}$Se$_x$ superconductors is that the top of innermost hole-like band moves away from the Fermi level to higher binding energy on going from $Γ$ to $Z$, quite opposite to the iron pnictides. The polarization dependence of the observed features was used to aid the attribution of the orbital character of the observed bands. Photon energy dependent measurements suggest a weak $k_z$ dispersion for the outer hole pocket and a moderate $k_z$ dispersion for the inner hole pocket. By evaluating the momentum and energy dependent spectral widths, the single-particle self-energy was extracted and interestingly this shows a pronounced non-Fermi liquid behaviour for these compounds. The experimental observations are discussed in context of electronic band structure calculations and models for the self-energy such as the spin-fermion model and the marginal-Fermi liquid.

preprint2015arXiv

Measurement of $\mathcal B(ψ(3770)\toγχ_{c1})$ and search for $ψ(3770)\toγχ_{c2}$

We report a measurement of the branching fraction of $ψ(3770)\toγχ_{c1}$ and search for the transition $ψ(3770)\rightarrow γχ_{c2}$ based on 2.92~fb$^{-1}$ of $e^+e^-$ data accumulated at $\sqrt{s}=3.773$~GeV with the BESIII detector at the BEPCII collider. The measured branching fraction of $ψ(3770) \rightarrow γχ_{c1}$ is $\mathcal B(ψ(3770) \rightarrow γχ_{c1})=(2.48 \pm 0.15 \pm 0.23) \times 10^{-3}$, which is the most precise measurement to date. The upper limit on the branching fraction of $ψ(3770)\rightarrow γχ_{c2}$ at a $90\%$ confidence level is $\mathcal B(ψ(3770) \rightarrow γχ_{c2})<0.64 \times 10^{-3}$. The corresponding partial widths are $Γ(ψ(3770) \to γχ_{c1}) =(67.5\pm 4.1\pm 6.7)$~keV and $Γ(ψ(3770) \to γχ_{c2}) < 17.4$~keV.

preprint2015arXiv

Measurement of $y_{CP}$ in $D^0-\overline{D}^0$ oscillation using quantum correlations in $e^+e^-\to D^0\overline{D}^0$ at $\sqrt{s}$ = 3.773\,GeV

We report a measurement of the parameter $y_{CP}$ in $\ensuremath{D^0}\xspace-\ensuremath{\overline{D}^{0}}\xspace$ oscillations performed by taking advantage of quantum coherence between pairs of $\ensuremath{D^0}\xspace \ensuremath{\overline{D}^{0}}\xspace$ mesons produced in $e^+e^-$ annihilations near threshold. In this work, doubly-tagged $\ensuremath{D^0}\xspace \ensuremath{\overline{D}^{0}}\xspace$ events, where one $D$ decays to a $CP$ eigenstate and the other $D$ decays in a semileptonic mode, are reconstructed using a data sample of 2.92\,fb$^{-1}$ collected with the BESIII detector at the center-of-mass energy of $\sqrt{s}$ = 3.773\,GeV. We obtain $y_{CP} = (-2.0\pm1.3\pm0.7)\%$, where the first uncertainty is statistical and the second is systematic. This result is compatible with the current world average.

preprint2015arXiv

Measurement of the $\mathrm e^+\mathrm e^-\rightarrow\mathrmπ^+\mathrmπ^-$ Cross Section between 600 and 900 MeV Using Initial State Radiation

We extract the $e^+e^-\rightarrow π^+π^-$ cross section in the energy range between 600 and 900 MeV, exploiting the method of initial state radiation. A data set with an integrated luminosity of 2.93 fb$^{-1}$ taken at a center-of-mass energy of 3.773 GeV with the BESIII detector at the BEPCII collider is used. The cross section is measured with a systematic uncertainty of 0.9%. We extract the pion form factor $|F_π|^2$ as well as the contribution of the measured cross section to the leading order hadronic vacuum polarization contribution to $(g-2)_μ$. We find this value to be $a_μ^{ππ,\rm LO}(600-900\;\rm MeV) = (368.2 \pm 2.5_{\rm stat} \pm 3.3_{\rm sys})\cdot 10^{-10}$.

preprint2015arXiv

Measurement of the $e^{+}e^{-} \to ηJ/ψ$ cross section and search for $e^{+}e^{-} \to π^{0} J/ψ$ at center-of-mass energies between 3.810 and 4.600~GeV

Using data samples collected with the BESIII detector operating at the BEPCII collider at center-of-mass energies from 3.810 to 4.600 GeV, we perform a study of $e^{+}e^{-} \to ηJ/ψ$ and $π^0 J/ψ$. Statistically significant signals of $e^{+}e^{-} \to ηJ/ψ$ are observed at $\sqrt{s}$ = 4.190, 4.210, 4.220, 4.230, 4.245, 4.260, 4.360 and 4.420 GeV, while no signals of $e^{+}e^{-} \to π^{0} J/ψ$ are observed. The measured energy-dependent Born cross section for $e^{+}e^{-} \to ηJ/ψ$ shows an enhancement around 4.2~GeV. The measurement is compatible with an earlier measurement by Belle, but with a significantly improved precision.

preprint2015arXiv

Measurement of the branching fraction for $ψ(3770) \to γχ_{c0}$

By analyzing a data set of 2.92 fb$^{-1}$ of $e^+e^-$ collision data taken at $\sqrt s= 3.773~\rm GeV$ and 106.41$\times 10^{6}$ $ψ(3686)$ decays taken at $\sqrt s= 3.686~\rm GeV$ with the BESIII detector at the BEPCII collider, we measure the branching fraction and the partial decay width for $ψ(3770)\toγχ_{c0}$ to be ${\mathcal B}(ψ(3770)\toγχ_{c0})=(6.88\pm0.28\pm0.67)\times 10^{-3}$ and $Γ[ψ(3770)\toγχ_{c0}]=(187\pm8\pm19)~\rm keV$, respectively. These are the most precise measurements to date.

preprint2015arXiv

Measurement of the branching fractions of $D_{s}^{+}\rightarrow η'X$ and $D_{s}^{+}\rightarrow η'ρ^{+}$ in $e^+e^-\to D^+_{s}D^-_{s}$

We study $D_{s}^{+}$ decays to final states involving the $η'$ with a 482$\,$pb$^{-1}$ data sample collected at $\sqrt{s}$ = 4.009$\,$GeV with the \mbox{BESIII} detector at the BEPCII collider. We measure the branching fractions $\mathcal{B}(D^+_{s}\rightarrow η'X)$ = (8.8$\pm$1.8$\pm$0.5)$\%$ and $\mathcal{B}(D_{s}^{+}\rightarrow η'ρ^{+})$ = ($5.8\pm1.4\pm0.4$)$\%$ where the first uncertainty is statistical and the second is systematic. In addition, we estimate an upper limit on the non-resonant branching ratio $\mathcal{B}(D_{s}^{+}\rightarrow η'π^+π^0)<5.1\%$ at the 90$\%$ confidence level. Our results are consistent with CLEO's recent measurements and help to resolve the disagreement between the theoretical prediction and CLEO's previous measurement of $\mathcal{B}(D_{s}^{+}\rightarrow η'ρ^{+})$.

preprint2015arXiv

Measurement of the form factors in the decay $D^+ \to ωe^+ ν_{e}$ and search for the decay $D^+ \to ϕe^+ ν_{e}$

Using 2.92 fb$^{-1}$ of electron-positron annihilation data collected at a center-of-mass energy of $\sqrt{s}= 3.773$ GeV with the BESIII detector, we present an improved measurement of the branching fraction $\mathcal{B}(D^+ \to ωe^+ ν_{e}) = (1.63\pm0.11\pm0.08)\times 10^{-3}$. The parameters defining the corresponding hadronic form factor ratios at zero momentum transfer are determined for the first time, we measure them to be $r_V = 1.24\pm0.09\pm0.06$ and $r_2 = 1.06\pm0.15 \pm 0.05$. The first and second uncertainties are statistical and systematic, respectively. We also search for the decay $D^+ \to ϕe^+ ν_{e}$. An improved upper limit $\mathcal{B}(D^+ \to ϕe^+ ν_{e}) < 1.3 \times 10^{-5}$ is set at 90\% confidence level.

preprint2015arXiv

Measurement of the Matrix Elements for the Decays $η\rightarrow π^{+}π^{-}π^0$ and $η/η^{\prime}\rightarrowπ^0π^0π^0$

Based on a sample of $1.31 \times 10^9$ $J/ψ$ events collected with the BESIII detector at the BEPCII collider, Dalitz plot analyses of selected 79,625 $η\rightarrowπ^{+}π^{-}π^0$ events, 33,908 $η\rightarrowπ^0π^0π^0$ events and 1,888 $η^{\prime}\rightarrowπ^0π^0π^0$ events are performed. The measured matrix elements of $η\rightarrowπ^+π^-π^0$ are in reasonable agreement with previous measurements. The Dalitz plot slope parameters of $η\rightarrowπ^0π^0π^0$ and $η^{\prime}\rightarrowπ^0π^0π^0$ are determined to be $-0.055 \pm 0.014 \pm 0.004$ and $-0.640 \pm 0.046 \pm 0.047$, respectively, where the first uncertainties are statistical and the second systematic. Both values are consistent with previous measurements, while the precision of the latter one is improved by a factor of three. Final state interactions are found to have an important role in those decays.

preprint2015arXiv

Measurement of the proton form factor by studying $e^{+} e^{-}\rightarrow p\bar{p}$

Using data samples collected with the BESIII detector at the BEPCII collider, we measure the Born cross section of $e^{+}e^{-}\rightarrow p\bar{p}$ at 12 center-of-mass energies from 2232.4 to 3671.0 MeV. The corresponding effective electromagnetic form factor of the proton is deduced under the assumption that the electric and magnetic form factors are equal $(|G_{E}|= |G_{M}|)$. In addition, the ratio of electric to magnetic form factors, $|G_{E}/G_{M}|$, and $|G_{M}|$ are extracted by fitting the polar angle distribution of the proton for the data samples with larger statistics, namely at $\sqrt{s}=$ 2232.4 and 2400.0 MeV and a combined sample at $\sqrt{s}$ = 3050.0, 3060.0 and 3080.0 MeV, respectively. The measured cross sections are in agreement with recent results from BaBar, improving the overall uncertainty by about 30\%. The $|G_{E}/G_{M}|$ ratios are close to unity and consistent with BaBar results in the same $q^{2}$ region, which indicates the data are consistent with the assumption that $|G_{E}|=|G_{M}|$ within uncertainties.

preprint2015arXiv

Measurements of $ψ(3686) \to K^{-} Λ\barΞ^{+} +c.c.$ and $ψ(3686) \to γK^{-} Λ\barΞ^{+} +c.c$

Using a sample of $1.06\times10^8\ \psip$ events produced in $e^+e^-$ collisions at $\sqrt{s}$ = 3.686 GeV and collected with the BESIII detector at the BEPCII collider, we present studies of the decays $\klx+c.c.$ and $\gklx+c.c.$. We observe two hyperons, $Ξ(1690)^-$ and $Ξ(1820)^-$, in the $K^-Λ$ invariant mass distribution in the decay $\klx+c.c.$ with significances of $4.9 σ$ and $6.2 σ$, respectively. The branching fractions of $\klx+c.c.$, $\ksx+c.c.$, $\psip\toγχ_{cJ}\to γK^- Λ\barΞ^+ +c.c.$ $(J=0,\ 1,\ 2)$, and $\psip\to Ξ(1690/1820)^{-} \barΞ^++c.c$ with subsequent decay $Ξ(1690/1820)^-\to K^-Λ$ are measured for the first time.

preprint2015arXiv

Observation and Spin-Parity Determination of the $X(1835)$ in $J/ψ\rightarrowγK^0_S K^0_Sη$

We report an observation of the process $J/ψ\rightarrowγX(1835)\rightarrowγK^0_S K^0_Sη$ at low $K^0_S K^0_S$ mass with a statistical significance larger than 12.9$σ$ using a data sample of $1.31 \times 10^{9}$ $J/ψ$ events collected with the BESIII detector. In this region of phase space the $K^0_S K^0_S$ system is dominantly produced through the $f_0(980)$. By performing a partial wave analysis, we determine the spin-parity of the $X(1835)$ to be $J^{PC}=0^{-+}$. The mass and width of the observed $X(1835)$ are $1844\pm9(\text{stat})^{+16}_{-25}(\text{syst})$~\MeV~and $192^{+20}_{-17}(\text{stat})^{+62}_{-43}(\text{syst})$~MeV, respectively, which are consistent with the results obtained by BESIII in the channel $J/ψ\rightarrowγπ^+π^-η^{\prime}$.

preprint2015arXiv

Observation of $η^{\prime}\toωe^{+} e^{-}$

Based on a sample of \etapr mesons produced in the radiative decay $J/ψ\toγη^{\prime}$ in $1.31\times 10^9$ $J/ψ$ events collected with the BESIII detector, the decay $η^{\prime}\toωe^{+} e^{-}$ is observed for the first time, with a statistical significance of $8σ$. The branching fraction is measured to be $\mathcal{B}(η^{\prime}\toωe^{+} e^{-})=(1.97\pm0.34(\text{stat})\pm0.17(\text{syst}))\times10^{-4}$, which is in agreement with theoretical predictions. The branching fraction of $η^{\prime}\toωγ$ is also measured to be $(2.55\pm0.03(\text{stat})\pm0.16(\text{syst}))\times10^{-2}$, which is the most precise measurement to date, and the relative branching fraction $\frac{\mathcal{B}(η^{\prime}\to ωe^{+}e^{-})}{\mathcal{B}(η^{\prime}\to ωγ)}$ is determined to be $(7.71\pm1.34(\text{stat})\pm0.54(\text{syst}))\times10^{-3}$.

preprint2015arXiv

Observation of a neutral charmoniumlike state $Z_c(4025)^0$ in $e^{+} e^{-} \to (D^{*} \bar{D}^{*})^{0} π^0$

We report a study of the process $e^{+} e^{-} \to (D^{*} \bar{D}^{*})^{0} π^0$ using $e^+e^-$ collision data samples with integrated luminosities of $1092 \rm{pb}^{-1}$ at $\sqrt{s}=4.23 \rm{GeV}$ and $826 \rm{pb}^{-1}$ at $\sqrt{s}=4.26 \rm{GeV}$ collected with the BESIII detector at the BEPCII storage ring. We observe a new neutral structure near the $(D^{*} \bar{D}^{*})^{0}$ mass threshold in the $π^0$ recoil mass spectrum, which we denote as $Z_{c}(4025)^{0}$. Assuming a Breit-Wigner line shape, its pole mass and pole width are determined to be $(4025.5^{+2.0}_{-4.7}\pm3.1) \rm{MeV}/c^2$ and $(23.0\pm 6.0\pm 1.0) \rm{MeV}$, respectively. The Born cross sections of $e^{+}e^{-}\to Z_{c}(4025)^{0} π^0\to (D^{*} \bar{D}^{*})^{0}π^0$ are measured to be $(61.6\pm8.2\pm9.0) \rm{pb}$ at $\sqrt{s}=4.23 \rm{GeV}$ and $(43.4\pm8.0\pm5.4) \rm{pb}$ at $\sqrt{s}=4.26 \rm{GeV}$. The first uncertainties are statistical and the second are systematic.

preprint2015arXiv

Observation of a Neutral Structure near the $D\bar{D}^{*}$ Mass Threshold in $e^{+}e^{-}\to (D \bar{D}^*)^0π^0$ at $\sqrt{s}$ = 4.226 and 4.257 GeV

A neutral structure in the $D\bar{D}^{*}$ system around the $D\bar{D}^{*}$ mass threshold is observed with a statistical significance greater than 10$σ$ in the processes $e^{+}e^{-}\rightarrow D^{+}D^{*-}π^{0}+c.c.$ and $e^{+}e^{-}\rightarrow D^{0}\bar{D}^{*0}π^{0}+c.c.$ at $\sqrt{s}$ = 4.226 and 4.257 GeV in the BESIII experiment. The structure is denoted as $Z_{c}(3885)^{0}$. Assuming the presence of a resonance, its pole mass and width are determined to be ($3885.7^{+4.3}_{-5.7}$(stat)$\pm 8.4$(syst)) MeV/$c^{2}$ and ($35^{+11}_{-12}$(stat)$ \pm 15$(syst)) MeV, respectively. The Born cross sections are measured to be $σ(e^{+}e^{-}\to Z_{c}(3885)^{0}π^{0}, Z_{c}(3885)^{0} \to D\bar{D}^{*})=(77 \pm 13$(stat)$\pm 17$(syst)) pb at 4.226 GeV and ($47 \pm 9$(stat)$ \pm 10$(syst)) pb at 4.257 GeV. The ratio of decay rates $\frac{\mathcal{B}({Z_{c}(3885)^{0} \to D^{+}D^{*-}+c.c.})}{\mathcal{B}({Z_{c}(3885)^{0} \to D^{0}\bar{D}^{*0}+c.c.})}$ is determined to be $0.96 \pm 0.18$(stat)$\pm 0.12$(syst), consistent with no isospin violation in the process $Z_{c}(3885)^{0}\to D\bar{D}^*$.

preprint2015arXiv

Observation of the $ψ(1^3D_2)$ state in $e^+e^-\toπ^+π^-γχ_{c1}$ at BESIII

We report the observation of the $X(3823)$ in the process $e^+e^-\to π^+π^-X(3823) \to π^+π^-γχ_{c1}$ with a statistical significance of $6.2σ$, in data samples at center-of-mass energies $\sqrt{s}=$4.230, 4.260, 4.360, 4.420 and 4.600~GeV collected with the BESIII detector at the BEPCII electron positron collider. The measured mass of the $X(3823)$ is $(3821.7\pm 1.3\pm 0.7)$~MeV/$c^2$, where the first error is statistical and the second systematic, and the width is less than $16$~MeV at the 90\% confidence level. The products of the Born cross sections for $e^+e^-\to π^+π^-X(3823)$ and the branching ratio $\mathcal{B}[X(3823)\to γχ_{c1,c2}]$ are also measured. These measurements are in good agreement with the assignment of the $X(3823)$ as the $ψ(1^3D_2)$ charmonium state.

preprint2015arXiv

Observation of the Dalitz Decay $η' \to γe^+e^-$

We report the first observation of the Dalitz decay $η' \to γe^+e^-$, based on a data sample of 1.31 billion $J/ψ$ events collected with the BESIII detector. The $η'$ mesons are produced via the $J/ψ\to γη'$ decay process. The ratio $Γ(η' \to γe^+ e^-)/Γ(η'\toγγ)$ is measured to be $(2.13\pm0.09(\text{stat.})\pm0.07(\text{sys.}))\times10^{-2}$. This corresponds to a branching fraction ${\cal B}(η' \to γe^+e^-)= (4.69 \pm0.20(\text{stat.})\pm0.23(\text{sys.}))\times10^{-4}$. The transition form factor is extracted and different expressions are compared to the measured dependence on the $e^+e^-$ invariant mass. The results are consistent with the prediction of the Vector Meson Dominance model.

preprint2015arXiv

Observation of the electromagnetic doubly OZI-suppressed decay $J/ψ\rightarrow ϕπ^{0}$

Using a sample of $1.31$ billion $J/ψ$ events accumulated with the BESIII detector at the BEPCII collider, we report the observation of the decay $J/ψ\rightarrow ϕπ^{0}$, which is the first evidence for a doubly Okubo-Zweig-Iizuka suppressed electromagnetic $J/ψ$ decay. A clear structure is observed in the $K^{+} K^{-}$ mass spectrum around 1.02 GeV/$c^2$, which can be attributed to interference between $J/ψ\rightarrow ϕπ^{0}$ and $J/ψ\rightarrow K^{+}K^{-}π^{0}$ decays. Due to this interference, two possible solutions are found. The corresponding measured values of the branching fraction of $J/ψ\to ϕπ^{0}$ are $[2.94 \pm 0.16\text{(stat.)} \pm 0.16\text{(syst.)}] \times 10^{-6}$ and $[1.24 \pm 0.33\text{(stat.)} \pm 0.30\text{(syst.)}] \times 10^{-7}$.

preprint2015arXiv

Observation of the isospin-violating decay $J/ψ\to ϕπ^{0}f_{0}(980)$

Using a sample of 1.31 billion $J/ψ$ events collected with the BESIII detector at the BEPCII collider, the decays $J/ψ\to ϕπ^{+}π^{-}π^{0}$ and $J/ψ\to ϕπ^{0}π^{0}π^{0}$ are investigated. The isospin violating decay $J/ψ\to ϕπ^{0} f_{0}(980)$ with $f_{0}(980) \to ππ$, is observed for the first time. The width of the $f_{0}(980)$ obtained from the dipion mass spectrum is found to be much smaller than the world average value. In the $π^{0} f_{0}(980)$ mass spectrum, there is evidence of $f_1(1285)$ production. By studying the decay $J/ψ\to ϕη'$, the branching fractions of $η' \to π^{+}π^{-}π^{0}$ and $η' \to π^{0}π^{0}π^{0}$, as well as their ratio, are also measured.

preprint2015arXiv

Observation of the Singly Cabibbo-Suppressed Decay $D^{+}\toωπ^{+}$ and Evidence for $D^{0}\toωπ^{0}$

Based on 2.93 fb$^{-1}$ $e^+e^-$ collision data taken at center-of-mass energy of 3.773 GeV by the BESIII detector, we report searches for the singly Cabibbo-suppressed decays $D^{+}\toωπ^{+}$ and $D^{0}\toωπ^{0}$. A double tag technique is used to measure the absolute branching fractions $\mathcal{B}(D^{+}\toωπ^{+})=(2.79\pm0.57\pm0.16)\times 10^{-4}$ and $\mathcal{B}(D^{0}\toωπ^{0})=(1.17\pm0.34\pm0.07)\times 10^{-4}$, with statistical significances of $5.5σ$ and $4.1σ$, respectively. We also present measurements of the absolute branching fractions for the related $ηπ$ decay modes. We find $\mathcal{B}(D^{+}\toηπ^{+})=(3.07\pm0.22\pm0.13)\times10^{-3}$ and $\mathcal{B}(D^{0}\toηπ^{0})=(0.65\pm0.09\pm0.04)\times10^{-3}$, which are consistent with the current world averages. The first and second uncertainties are statistical and systematic, respectively.

preprint2015arXiv

Precision measurement of the $D^{*0}$ decay branching fractions

Using 482 pb$^{-1}$ of data taken at $\sqrt{s}=4.009$ GeV, we measure the branching fractions of the decays of $D^{*0}$ into $D^0π^0$ and $D^0γ$ to be $\BR(D^{*0} \to D^0π^0)=(65.5\pm 0.8\pm 0.5)%$ and $\BR(D^{*0} \to D^0γ)=(34.5\pm 0.8\pm 0.5)%$ respectively, by assuming that the $D^{*0}$ decays only into these two modes. The ratio of the two branching fractions is $\BR(D^{*0} \to D^0π^0)/\BR(D^{*0} \to D^0γ) =1.90\pm 0.07\pm 0.05$, which is independent of the assumption made above. The first uncertainties are statistical and the second ones systematic. The precision is improved by a factor of three compared to the present world average values.

preprint2015arXiv

Precision measurement of the integrated luminosity of the data taken by BESIII at center of mass energies between 3.810 GeV and 4.600 GeV

From December 2011 to May 2014, about 5 $\rm fb^{-1}$ of data were taken with the BESIII detector at center-of-mass energies between 3.810 GeV and 4.600 GeV to study the charmoniumlike states and higher excited charmonium states. The time integrated luminosity of the collected data sample is measured to a precision of 1% by analyzing events produced by the large-angle Bhabha scattering process.

preprint2015arXiv

Search for $D^0\toγγ$ and improved measurement of the branching fraction for $D^0 \to π^0π^0$

Using $2.92$ fb$^{-1}$ of electron-positron annihilation data collected at $\sqrt{s} = 3.773$~GeV with the BESIII detector, we report the results of a search for the flavor-changing neutral current process $D^0\toγγ$ using a double-tag technique. We find no signal and set an upper limit at $90\%$ confidence level for the branching fraction of $B(D^0\toγγ) < 3.8\times10^{-6}$. We also investigate $D^0$-meson decay into two neutral pions, obtaining a branching fraction of $B(D^0\toπ^0π^0) = (8.24\pm0.21(\text{stat.})\pm0.30(\text{syst.}))\times10^{-4}$, the most precise measurement to date and consistent with the current world average

preprint2015arXiv

Search for $Z_c(3900)^\pm\toωπ^\pm$

The decay $Z_c(3900)^\pm\toωπ^\pm$ is searched for using data samples collected with the BESIII detector operating at the BEPCII storage ring at center-of-mass energies $\sqrt{s}=4.23$ and $4.26$~GeV. No significant signal for the $Z_c(3900)^\pm$ is found, and upper limits at the 90\% confidence level on the Born cross section for the process $e^+e^-\to Z_c(3900)^\pmπ^\mp\toωπ^+π^-$ are determined to be $0.26$ and $0.18$ pb at $\sqrt{s}=4.23$ and 4.26 GeV, respectively.

preprint2015arXiv

Search for the $Y(4140)$ via $e^{+} e^{-} \to γϕJ/ψ$ at $\sqrt{s}=$ 4.23, 4.26 and 4.36 GeV

Using data samples collected at center-of-mass energies $\sqrt{s} = 4.23$, 4.26, and 4.36 GeV with the BESIII detector operating at the BEPCII storage ring, we search for the production of the charmoniumlike state $Y(4140)$ through a radiative transition followed by its decay to $ϕJ/ψ$. No significant signal is observed and upper limits on $σ[e^{+} e^{-} \rightarrow γY(4140)] \cdot \mathcal{B}(Y(4140)\rightarrow ϕJ/ψ)$ at the $90\%$ confidence level are estimated as 0.35, 0.28, and 0.33 pb at $\sqrt{s} = 4.23$, 4.26, and 4.36 GeV, respectively.

preprint2015arXiv

Search for the $η_c(2S)/h_c\to p\bar{p}$ decays and measurements of the $χ_{cJ}\to p\bar{p}$ branching fractions

Using a sample of $1.06\times10^{8}~ψ(3686)$ events collected with the BESIII detector at BEPCII, the decays $η_{c}(2S)\to p\bar{p}$ and $h_{c}\to p\bar{p}$ are searched for, where $η_c(2S)$ and $h_c$ are reconstructed in the decay chains $ψ(3686)\toγη_{c}(2S)$, $η_{c}(2S)\to p\bar{p}$ and $ψ(3686)\toπ^{0}h_{c}$, $h_{c}\to p\bar{p}$, respectively. No significant signals are observed. The upper limits of the product branching fractions are determined to be $\mathcal{B}(ψ(3686)\toγη_c(2S))\times\mathcal{B}(η_{c}(2S)\to p\bar{p})<1.4\times10^{-6}$ and $\mathcal{B}(ψ(3686)\toπ^0h_c)\times\mathcal{B}(h_{c}\to p\bar{p})<1.3\times10^{-7}$ at the 90% C.L.. The branching fractions for $χ_{cJ}\to p\bar{p}$ $(J=0,~1,~2)$ are also measured to be $(24.5\pm0.8\pm1.3,~8.6\pm0.5\pm0.5,~8.4\pm0.5\pm0.5)\times10^{-5}$, which are the world's most precise measurements.

preprint2015arXiv

Search for the isospin violating decay $Y(4260)\rightarrow J/ψηπ^{0}$

Using data samples collected at center of mass energies of $\sqrt{s}$ = 4.009, 4.226, 4.257, 4.358, 4.416 and 4.599 GeV with the BESIII detector operating at the BEPCII storage ring, we search for the isospin violating decay $Y(4260)\rightarrow J/ψηπ^{0}$. No signal is observed, and upper limits on the cross section $σ(e^{+}e^{-}\rightarrow J/ψηπ^{0})$ at the 90\% confidence level are determined to be 3.6, 1.7, 2.4, 1.4, 0.9 and 1.9 pb, respectively.

preprint2015arXiv

Search for the weak decays $J/ψ\to D^{(*)}_{s}eν_{e}$

Using a sample of $2.25\times 10^8$ $J/ψ$ events collected with the BESIII detector at the BEPCII collider, we search for the $J/ψ$ semi-leptonic weak decay $J/ψ\to D^{-}_{s} e^{+}ν_{e}+c.c.$ with a much higher sensitivity than previous searches. We also perform the first search for $J/ψ\to D^{*-}_{s} e^{+}ν_{e}+c.c.$ No significant excess of a signal above background is observed in either channel. At the $90\%$ confidence level, the upper limits are determined to be $\mathcal{B}(J/ψ\to D^{-}_{s}e^{+}ν_{e}+c.c.)<1.3\times10^{-6}$ and $\mathcal{B}(J/ψ\to {D^{*}_{s}}^{-}e^{+}ν_{e}+c.c.)<1.8\times10^{-6}$, respectively. Both are consistent with Standard Model predictions.

preprint2015arXiv

Searches for isospin-violating transitions $χ_{c0,2} \rightarrow π^{0} η_{c}$

We present the first upper-limit measurement of the branching fractions of the isospin-violating transitions $χ_{c0,2} \rightarrow π^{0} η_{c}$. The measurements are performed using $106\times 10^{6}$ $ψ(3686)$ events accumulated with the BESIII detector at the BEPCII $e^{+}e^{-}$ collider at a center-of-mass energy corresponding to the $ψ(3686)$ mass. We obtained upper limits on the branching fractions at a 90\% confidence level of $B(χ_{c0} \rightarrow π^{0} η_{c}) < 1.6 \times 10^{-3}$ and $B(χ_{c2} \rightarrow π^{0} η_{c}) < 3.2 \times 10^{-3}$.

preprint2015arXiv

Study of $e^+e^-\toωχ_{cJ}$ at center-of-mass energies from 4.21 to 4.42 GeV

Based on data samples collected with the BESIII detector at the BEPCII collider at 9 center-of-mass energies from 4.21 to 4.42 GeV, we search for the production of $e^+e^-\to ωχ_{cJ}$ ($J$=0, 1, 2). The process $e^+e^-\to ωχ_{c0}$ is observed for the first time, and the Born cross sections at $\sqrt{s}$=4.23 and 4.26 GeV are measured to be $(55.4\pm 6.0\pm 5.9)$ and $(23.7\pm 5.3\pm 3.5)$ pb, respectively, where the first uncertainties are statistical and the second are systematic. The $ωχ_{c0}$ signals at the other 7 energies and $e^+e^-\to ωχ_{c1}$ and $ωχ_{c2}$ signals are not significant, and the upper limits on the cross sections are determined. By examining the $ωχ_{c0}$ cross section as a function of center-of-mass energy, we find that it is inconsistent with the line shape of the $Y(4260)$ observed in $e^+ e^-\toπ^+π^-J/ψ$. Assuming the $ωχ_{c0}$ signals come from a single resonance, we extract mass and width of the resonance to be $(4230\pm8\pm6)$ MeV/$c^2$ and $(38\pm12\pm2)$ MeV, respectively, and the statistical significance is more than $9σ$.

preprint2015arXiv

Study of $χ_{cJ}$ decaying into $ϕ$ $K^*(892)$ $\bar{K}$

Using a data sample of 106 million $ψ(3686)$ events collected with the BESIII detector operated at the BEPCII storage ring, we study for the first time the decay $χ_{cJ}\toϕK^{0}_S K^{\pm}π^{\mp}$ and $χ_{cJ}\toϕK^{+} K^{-}π^{0}$ in the E1 radiative transition $ψ(3686)\toγχ_{cJ}$. The decays are dominated by the three-body decay $χ_{cJ}\to ϕK^*(892)\bar{K}$. We measure branching fractions for this reaction via the neutral and charged $K^*(892)$ and find them consistent with each other within the expectation of isospin symmetry. In the $K\bar{K}π$ invariant mass distribution a structure near the $K^*(892)\bar{K}$ mass threshold is observed, and the corresponding mass and width are measured to be $1412\pm4(\mathrm{stat.})\pm8(\mathrm{sys.}) \mathrm{MeV}/c^2$ and $Γ$ = $84\pm12(\mathrm{stat.})\pm40(\mathrm{sys.}) \mathrm{MeV}$, respectively. The observed state favors an assignment to the $h_1(1380)$, considering its possible $J^{PC}$ and comparing its mass, width and decay mode to those reported in the Particle Data Group.

preprint2015arXiv

Study of decay dynamics and $CP$ asymmetry in $D^+ \to K^0_L e^+ ν_e$ decay

Using 2.92~fb$^{-1}$ of electron-positron annihilation data collected at $\sqrt{s} = 3.773$~GeV with the BESIII detector, we obtain the first measurements of the absolute branching fraction $\mathcal{B}(D^+ \to K^0_L e^+ ν_e) = (4.481 \pm 0.027(\mathrm{stat.}) \pm 0.103(\mathrm{sys.}))\%$ and the $CP$ asymmetry $A_{CP}^{D^+ \to K^0_L e^+ ν_e} = (-0.59 \pm 0.60(\mathrm{stat.}) \pm 1.48(\mathrm{sys.}))\%$. From the $D^+ \to K^0_L e^+ ν_e$ differential decay rate distribution, the product of the hadronic form factor and the magnitude of the CKM matrix element, $f_{+}^{K}(0)|V_{cs}|$, is determined to be $0.728 \pm 0.006(\mathrm{stat.}) \pm 0.011(\mathrm{sys.})$. Using $|V_{cs}|$ from the SM constrained fit with the measured $f_{+}^{K}(0)|V_{cs}|$, $f_{+}^{K}(0) = 0.748 \pm 0.007(\mathrm{stat.}) \pm 0.012(\mathrm{sys.})$ is obtained, and utilizing the unquenched LQCD calculation for $f_{+}^{K}(0)$, $|V_{cs}| = 0.975 \pm 0.008(\mathrm{stat.}) \pm 0.015(\mathrm{sys.}) \pm 0.025(\mathrm{LQCD})$.

preprint2015arXiv

Study of Dynamics of $D^0 \to K^- e^+ ν_{e}$ and $D^0\toπ^- e^+ ν_{e}$ Decays

In an analysis of a 2.92~fb$^{-1}$ data sample taken at 3.773~GeV with the BESIII detector operated at the BEPCII collider, we measure the absolute decay branching fractions to be $\mathcal B(D^0 \to K^-e^+ν_e)=(3.505\pm 0.014 \pm 0.033)\%$ and $\mathcal B(D^0 \to π^-e^+ν_e)=(0.295\pm 0.004\pm 0.003)\%$. From a study of the differential decay rates we obtain the products of hadronic form factor and the magnitude of the CKM matrix element $f_{+}^K(0)|V_{cs}|=0.7172\pm0.0025\pm 0.0035$ and $f_{+}^π(0)|V_{cd}|=0.1435\pm0.0018\pm 0.0009$. Combining these products with the values of $|V_{cs(d)}|$ from the SM constraint fit, we extract the hadronic form factors $f^K_+(0) = 0.7368\pm0.0026\pm 0.0036$ and $f^π_+(0) = 0.6372\pm0.0080\pm 0.0044$, and their ratio $f_+^π(0)/f_+^{K}(0)=0.8649\pm 0.0112\pm 0.0073$. These form factors and their ratio are used to test unquenched Lattice QCD calculations of the form factors and a light cone sum rule (LCSR) calculation of their ratio. The measured value of $f_+^{K(π)}(0) |V_{cs(d)}|$ and the lattice QCD value for $f^{K(π)}_+(0)$ are used to extract values of the CKM matrix elements of $|V_{cs}|=0.9601 \pm 0.0033 \pm 0.0047 \pm 0.0239$ and $|V_{cd}|=0.2155 \pm 0.0027 \pm 0.0014 \pm 0.0094$, where the third errors are due to the uncertainties in lattice QCD calculations of the form factors. Using the LCSR value for $f_+^π(0)/f_+^K(0)$, we determine the ratio $|V_{cd}|/|V_{cs}|=0.238\pm 0.004\pm 0.002\pm 0.011$, where the third error is from the uncertainty in the LCSR normalization. In addition, we measure form factor parameters for three different theoretical models that describe the weak hadronic charged currents for these two semileptonic decays. All of these measurements are the most precise to date.

preprint2015arXiv

The evolution of stellar metallicity gradients of the Milky Way disk from LSS-GAC main sequence turn-off stars: a two-phase disk formation history?

We use 297 042 main sequence turn-off stars selected from the LSS-GAC to determine the radial and vertical gradients of stellar metallicity of the Galactic disk in the anti-center direction. We determine ages of those turn-off stars by isochrone fitting and measure the temporal variations of metallicity gradients. Our results show that the gradients, both in the radial and vertical directions, exhibit significant spatial and temporal variations. The radial gradients yielded by stars of oldest ages (>11 Gyr) are essentially zero at all heights from the disk midplane, while those given by younger stars are always negative. The vertical gradients deduced from stars of oldest ages (>11Gyr) are negative and show only very weak variations with the Galactocentric distance in the disk plane, $R$, while those yielded by younger stars show strong variations with $R$. After being essentially flat at the earliest epochs of disk formation, the radial gradients steepen as age decreases, reaching a maxima (steepest) at age 7-8 Gyr, and then they flatten again. Similar temporal trends are also found for the vertical gradients. We infer that the assemblage of the Milky Way disk may have experienced at least two distinct phases. The earlier phase is probably related to a slow, pressure-supported collapse of gas, when the gas settles down to the disk mainly in the vertical direction. In the later phase, there are significant radial flows of gas in the disk, and the rate of gas inflow near the solar neighborhood reaches a maximum around a lookback time of 7-8 Gyr. The transition of the two phases occurs around a lookback time between 8 and 11 Gyr. The two phases may be responsible for the formation of the Milky Way thick and thin disks, respectively. And, as a consequence, we recommend that stellar age is a natural, physical criterion to distinguish thin and thick disk stars. ... (abridged)

preprint2015arXiv

The LAMOST Spectroscopic Survey of Globular Clusters in M 31 and M 33. I. Catalog and new identifications

We present a catalog of 908 objects observed with the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) in the vicinity fields of M31 and M33, targeted as globular clusters (GCs) and candidates. The targets include known GCs and candidates selected from the literature, as well as new candidates selected from the Sloan Digital Sky Survey (SDSS). Analysis shows that 356 of them are likely GCs of various degree of confidence, while the remaining ones turn out to be background galaxies and quasars, stars and HII regions in M31 or foreground Galactic stars. The 356 likely GCs include 298 bona fide GCs and 26 candidates known in the literature. Three candidates selected from the Revised Bologna Catalog of M31 GCs and candidates (RBC) and one possible cluster from Johnson et al. are confirmed to be bona fide clusters. We search for new GCs in the halo of the M31 amongst the new candidates selected from the SDSS photometry. Based on radial velocities yielded by LAMSOT spectra and visual examination of the SDSS images, we find 28 objects, 5 bona fide and 23 likely GCs. Amongst the five bona fide GCs, three have been recently discovered independently by others, the remaining 25 are our new identifications, including two bona fide ones. The new identified objects fall at projected distances ranging from 13 to 265 kpc from M31. Of the two newly discovered bona fide GCs, one is located near M33, probably a GC belonging to M33. The other bona fide GC falls on the Giant Stream with a projected distance of 78 kpc from M31. Of the 23 newly identified likely GCs, one has a projected distance of about 265 kpc from M31 and could be an intergalactic cluster.

preprint2015arXiv

Thermal and electromagnetic properties of Bi$_2$Sr$_2$CaCu$_2$O$_8$ intrinsic Josephson junction stacks studied via one-dimensional coupled sine-Gordon equations

We used one-dimensional coupled sine-Gordon equations combined with heat diffusion equations to numerically investigate the thermal and electromagnetic properties of a $300\,μ\mathrm{m}$ long intrinsic Josephson junction stack consisting of $N = 700$ junctions. The junctions in the stack are combined to $M$ segments where we assume that inside a segment all junctions behave identically. Most simulations are for $M = 20$. For not too high bath temperatures there is the appearence of a hot spot at high bias currents. In terms of electromagnetic properties, robust standing wave patterns appear in the current density and electric field distributions. These patterns come together with vortex/antivortex lines across the stack that correspond to $π$ kink states, discussed before in the literature for a homogeneous temperature distribution in the stack. We also discuss scaling of the thermal and electromagnetic properties with $M$, on the basis of simulations with $M$ between 10 and 350.

preprint2014arXiv

A three dimensional extinction map of the Galactic Anticentre from multi-band photometry

We present a three dimensional extinction map in $r$ band. The map has a spatial angular resolution, depending on latitude, between 3 -- 9\,arcmin and covers the entire XSTPS-GAC survey area of over 6,000\,$\rm deg^2$ for Galactic longitude $\rm 140 \leq$ $l$ $ \leq 220°$ and latitude $\rm -40\leq$ $b$ $ \leq 40°$. By cross-matching the photometric catalog of the Xuyi Schmidt Telescope Photometric Survey of the Galactic Anticentre (XSTPS-GAC) with those of 2MASS and WISE, we have built a multi-band photometric stellar sample of about 30 million stars and applied spectral energy distribution (SED) fitting to the sample. By combining photometric data from the optical to the near-infrared, we are able to break the degeneracy between the intrinsic stellar colours and the amounts of extinction by dust grains for stars with high photometric accuracy, and trace the extinction as a function of distance for low Galactic latitude and thus highly extincted regions. This has allowed us to derive the best-fit extinction and distance information of more than 13 million stars, which are used to construct the three dimensional extinction map. We have also applied a Rayleigh-Jeans colour excess (RJCE) method to the data using the 2MASS and WISE colour $(H-W2)$. The resulting RJCE extinction map is consistent with the integrated two dimensional map deduced using the best-fit SED algorithm. However for individual stars, the amounts of extinction yielded by the RJCE method suffer from larger errors than those given by the best-fit SED algorithm.

preprint2014arXiv

Amplitude Analysis of the $D^+ \to K^0_S π^+ π^0$ Dalitz Plot

We perform an analysis of the $D^+ \to K^0_S π^+ π^0$ Dalitz plot using a data set of 2.92 fb$^{-1}$ of $e^+e^-$ collisions at the $ψ(3770)$ mass accumulated by the BESIII Experiment, in which 166694 candidate events are selected with a background of 15.1%. The Dalitz plot is found to be well-represented by a combination of six quasi-two-body decay channels ($K^0_Sρ^+$, $K^0_Sρ(1450)^+$, $\overline{K}^{*0}π^+$, $\overline{K}_0(1430)^0π^+$, $\overline{K}(1680)^0π^+$, $\overlineκ^0π^+$) plus a small non-resonant component. Using the fit fractions from this analysis, partial branching ratios are updated with higher precision than previous measurements.

preprint2014arXiv

An integrated low phase noise radiation-pressure-driven optomechanical oscillator chipset

High-quality frequency references are the cornerstones in position, navigation and timing applications of both scientific and commercial domains. Optomechanical oscillators, with direct coupling to continuous-wave light and non-material-limited f Q product, are long regarded as a potential platform for frequency reference in radio-frequency-photonic architectures. However, one major challenge is the compatibility with standard CMOS fabrication processes while maintaining optomechanical high quality performance. Here we demonstrate the monolithic integration of photonic crystal optomechanical oscillators and on-chip high speed Ge detectors based on the silicon CMOS platform. With the generation of both high harmonics (up to 59th order) and subharmonics (down to 1/4), our chipset provides multiple frequency tones for applications in both frequency multipliers and dividers. The phase noise is measured down to -125 dBc/Hz at 10 kHz offset at ~ 400 μW dropped-in powers, one of the lowest noise optomechanical oscillators to date and in room-temperature and atmospheric non-vacuum operating conditions. These characteristics enable optomechanical oscillators as a frequency reference platform for radio-frequency-photonic information processing.

preprint2014arXiv

Coexistence of superconductivity and density wave in $\mathrm{Ba_2Ti_2Fe_2As_4O}$: an optical spectroscopy study

We performed an optical spectroscopy measurement on single crystals of $\mathrm{Ba_2Ti_2Fe_2As_4O}$, which is a newly discovered superconductor showing a coexistence of superconductivity and density wave order. The study reveals a spectral change related to the formation density wave energy gap below $T_{DW}\approx$125 K, leading to the removal of about half of the multiple Fermi surface sheets. The ratio of 2$Δ_{DW}$/$k_B T_{DW}\approx$ 11.9 is considerably larger than the mean-field value based on the weak-coupling BCS theory. At the lowest temperature in the superconducting state, we observed opening of superconducting energy gaps $Δ_1(0) =3.4$ meV and $Δ_2(0)=7.9$ meV. The properties of superconducting state in $\mathrm{Ba_2Ti_2Fe_2As_4O}$ are similar to that in $\mathrm{BaFe_{1.85}Co_{0.15}As_2}$.

preprint2014arXiv

Evidence for $e^+e^-\toγχ_{c1, 2}$ at center-of-mass energies from 4.009 to 4.360 GeV

Using data samples collected at center-of-mass energies of $\sqrt{s}$ = 4.009, 4.230, 4.260, and 4.360 GeV with the BESIII detector operating at the BEPCII collider, we perform a search for the process $e^+e^-\toγχ_{cJ}$ $(J = 0, 1, 2)$ and find evidence for $e^+e^-\toγχ_{c1}$ and $e^+e^-\toγχ_{c2}$ with statistical significances of 3.0$σ$ and 3.4$σ$, respectively. The Born cross sections $σ^{B}(e^+e^-\toγχ_{cJ})$, as well as their upper limits at the 90% confidence level are determined at each center-of-mass energy.

preprint2014arXiv

Exponential Convergence Towards Stationary States for the 1D Porous Medium Equation with Fractional Pressure

We analyse the asymptotic behaviour of solutions to the one dimensional fractional version of the porous medium equation introduced by Caffarelli and Vázquez, where the pressure is obtained as a Riesz potential associated to the density. We take advantage of the displacement convexity of the Riesz potential in one dimension to show a functional inequality involving the entropy, entropy dissipation, and the Euclidean transport distance. An argument by approximation shows that this functional inequality is enough to deduce the exponential convergence of solutions in self-similar variables to the unique steady states.

preprint2014arXiv

Geo-neutrinos and Earth Models

We present the current status of geo-neutrino measurements and their implications for radiogenic heating in the mantle. Earth models predict different levels of radiogenic heating and, therefore, different geo-neutrino fluxes from the mantle. Seismic tomography reveals features in the deep mantle possibly correlated with radiogenic heating and causing spatial variations in the mantle geo-neutrino flux at the Earth surface. An ocean-based observatory offers the greatest sensitivity to the mantle flux and potential for resolving Earth models and mantle features. Refinements to estimates of the geo-neutrino flux from continental crust reduce uncertainty in measurements of the mantle flux, especially measurements from land-based observatories. These refinements enable the resolution of Earth models using the combined measurements from multiple continental observatories.

preprint2014arXiv

HST and LAMOST discover a dual active galactic nucleus in J0038+4128

We report the discovery of a kiloparsec-scale dual active galactic nucleus (AGN) in J0038+4128. From the Hubble Space Telescope (HST) Wide Field Planetary Camera (WFPC2) images, we find two optical nuclei with a projection separation of 4.7 kpc (3.44 arcsec). The southern component (J0038+4128S) is spectroscopically observed with the HST Goddard High Resolution Spectrograph in the UV range and is found to be a Seyfert 1 galaxy with a broad Ly alpha emission line. The northern component (J0038+4128N) is spectroscopically observed during the Large Sky Area Multi-Object Fibre Spectroscopic Telescope (also named the Guoshoujing Telescope) pilot survey in the optical range. The observed line ratios as well as the consistency of redshift of the nucleus emission lines and the host galaxy's absorption lines indicate that J0038+4128N is a Seyfert 2 galaxy with narrow lines only. These results thus confirm that J0038+4128 is a Seyfert 1-Seyfert 2 AGN pair. The HST WFPC2 F336W/U-band image of J0038+4128 also reveals for the first time for a dual AGN system two pairs of bi-symmetric arms, as are expected from the numerical simulations of such system. Being one of a few confirmed kiloparsec-scale dual AGNs exhibiting a clear morphological structure of the host galaxies, J0038+4128 provides an unique opportunity to study the co-evolution of the host galaxies and their central supermassive black holes undergoing a merging process.

preprint2014arXiv

Itinerant effects and enhanced magnetic interactions in Bi-based multilayer cuprates

The cuprate high temperature superconductors exhibit a pronounced trend in which the superconducting transition temperature, $T_{\rm c}$, increases with the number of CuO$_2$ planes, $n$, in the crystal structure. We compare the magnetic excitation spectrum of Bi$_{2+x}$Sr$_{2-x}$CuO$_{6+δ}$ (Bi-2201) and Bi$_2$Sr$_2$Ca$_2$Cu$_3$O$_{10 + δ}$ (Bi-2223), with $n=1$ and $n=3$ respectively, using Cu $L_3$-edge resonant inelastic x-ray scattering (RIXS). Near the anti-nodal zone boundary we find the paramagnon energy in Bi-2223 is substantially higher than that in Bi-2201, indicating that multilayer cuprates host stronger effective magnetic exchange interactions, providing a possible explanation for the $T_{\rm c}$ vs.\ $n$ scaling. In contrast, the nodal direction exhibits very strongly damped, almost non-dispersive excitations. We argue that this implies that the magnetism in the doped cuprates is partially itinerant in nature.

preprint2014arXiv

Measurement of $χ_{cJ}$ decaying into $η^{\prime}K^+K^-$

Using $(106.41\pm 0.86) \times 10^{6}$ $ψ(3686)$ events collected with the BESIII detector at BEPCII, we study for the first time the decay $χ_{cJ}\toη'K^+K^-$, where $η'\toγρ^{0}$ and $η'\toηπ^+π^-$. A partial wave analysis in the covariant tensor amplitude formalism is performed for the decay $χ_{c1}\toη'K^+K^-$. Intermediate processes $χ_{c1}\toη' f_0(980)$, $χ_{c1}\toη' f_0(1710)$, $χ_{c1}\toη' f_2'(1525)$ and $χ_{c1}\to K^*_0(1430)^{\pm}K^{\mp}$ ($K^*_0(1430)^{\pm}\toη' K^{\pm}$) are observed with statistical significances larger than 5$σ$, and their branching fractions are measured.

preprint2014arXiv

Measurement of the $D\to K^-π^+$ strong phase difference in $ψ(3770)\to D^0\overline{D}{}^0$

We study $D^0\overline{D}{}^0$ pairs produced in $e^+e^-$ collisions at $\sqrt{s}=3.773$ GeV using a data sample of 2.92 fb$^{-1}$ collected with the BESIII detector. We measure the asymmetry $\mathcal{A}^{CP}_{Kπ}$ of the branching fractions of $D \to K^-π^+$ in $CP$-odd and $CP$-even eigenstates to be $(12.7\pm1.3\pm0.7)\times10^{-2}$. $\mathcal{A}^{CP}_{Kπ}$ can be used to extract the strong phase difference $δ_{Kπ}$ between the doubly Cabibbo-suppressed process $\overline{D}{}^{0}\to K^-π^+$ and the Cabibbo-favored process $D^0\to K^- π^+$. Using world-average values of external parameters, we obtain $\cosδ_{Kπ} = 1.02\pm0.11\pm0.06\pm0.01$. Here, the first and second uncertainties are statistical and systematic, respectively, while the third uncertainty arises from the external parameters. This is the most precise measurement of $δ_{Kπ}$ to date.

preprint2014arXiv

Measurement of the branching fraction for $ψ(3686)\toωK^+ K^-$

With $1.06\times 10^8$ $ψ(3686)$ events collected with the BESIII detector, the branching fraction of $ψ(3686) \to ωK^+ K^-$ is measured to be $(1.54 \pm 0.04 \pm 0.11) \times 10^{-4}$. This is the most precise result to date, due to the largest $ψ(3686)$ sample, improved signal reconstruction efficiency, good simulation of the detector performance, and a more accurate knowledge of the continuum contribution. Using the branching fraction of $J/ψ\to ωK^+ K^-$, the ratio $\mathcal{B}(ψ(3868) \to K^+K^-) / \mathcal{B}(J/ψ\to K^+K^-)$ is determined to be $(18.4 \pm 3.7)\,\%$. This constitutes a significantly improved test of the $12\,\%$ rule, with the uncertainty now dominated by the $J/ψ$ branching fraction.

preprint2014arXiv

Observation of $J/ψ\rightarrow p\bar{p}a_{0}(980)$ at BESIII

Using $2.25\times10^{8}$ $J/ψ$ events collected with the BESIII detector at the BEPCII storage rings, we observe for the first time the process $J/ψ\rightarrow p\bar{p}a_{0}(980)$, $a_{0}(980)\rightarrow π^{0}η$ with a significance of $6.5σ$ ($3.2σ$ including systematic uncertainties). The product branching fraction of $J/ψ\rightarrow p\bar{p}a_{0}(980)\rightarrow p\bar{p}π^{0}η$ is measured to be $(6.8\pm1.2\pm1.3)\times 10^{-5}$, where the first error is statistical and the second is systematic. This measurement provides information on the $a_{0}$ production near threshold coupling to $p\bar{p}$ and improves the understanding of the dynamics of $J/ψ$ decays to four body processes.

preprint2014arXiv

Observation of $η^{\prime}\toπ^{+}π^{-}π^{+}π^{-}$ and $η^{\prime}\toπ^{+}π^{-}π^{0}π^{0}$

Using a sample of $1.3\times 10^9$ $J/ψ$ events collected with the BESIII detector, we report the first observation of $η^{\prime}\toπ^{+}π^{-}π^{+}π^{-}$ and $η^{\prime}\toπ^{+}π^{-}π^{0}π^{0}$. The measured branching fractions are $\mathcal{B}$($η^{\prime}\toπ^{+}π^{-}π^{+}π^{-}$) = $(8.53\pm0.69({\rm stat.})\pm0.64({\rm syst.}))\times10^{-5}$ and $\mathcal{B}$($η^{\prime}\toπ^{+}π^{-}π^{0}π^{0}$) = $(1.82\pm0.35({\rm stat.})\pm0.18({\rm syst.}))\times10^{-4}$, which are consistent with theoretical predictions based on a combination of chiral perturbation theory and vector-meson dominance.

preprint2014arXiv

Observation of a charged charmoniumlike structure in $e^+e^- \to (D^{*} \bar{D}^{*})^{\pm} π^\mp$ at $\sqrt{s}=4.26$GeV

We study the process $e^+e^- \to (D^{*} \bar{D}^{*})^{\pm} π^\mp$ at a center-of-mass energy of 4.26GeV using a 827pb$^{-1}$ data sample obtained with the BESIII detector at the Beijing Electron Positron Collider. Based on a partial reconstruction technique, the Born cross section is measured to be $(137\pm9\pm15)$pb. We observe a structure near the $(D^{*} \bar{D}^{*})^{\pm}$ threshold in the $π^\mp$ recoil mass spectrum, which we denote as the $Z^{\pm}_c(4025)$. The measured mass and width of the structure are $(4026.3\pm2.6\pm3.7)$MeV/c$^2$ and $(24.8\pm5.6\pm7.7)$MeV, respectively. Its production ratio $\frac{σ(e^+e^-\to Z^{\pm}_c(4025)π^\mp \to (D^{*} \bar{D}^{*})^{\pm} π^\mp)}{σ(e^+e^-\to (D^{*} \bar{D}^{*})^{\pm} π^\mp)}$ is determined to be $0.65\pm0.09\pm0.06$. The first uncertainties are statistical and the second are systematic.

preprint2014arXiv

Observation of Electromagnetic Dalitz decays J/ψ\to P e^+e^-

Based on a sample of (225.3\pm2.8)\times 10^{6} J/ψevents collected with the BESIII detector, the electromagnetic Dalitz decays of J/ψ\to P e^+e^-(P=η'/η/π^0) are studied. By reconstructing the pseudoscalar mesons in various decay modes, the decays J/ψ\to η' e^+e^-, J/ψ\to ηe^+e^- and J/ψ\to π^0 e^+e^- are observed for the first time. The branching fractions are determined to be \mathcal{B}(J/ψ\to η' e^+e^-) = (5.81\pm0.16\pm0.31)\times10^{-5}, \mathcal{B}(J/ψ\to ηe^+e^-) = (1.16\pm0.07\pm0.06)\times10^{-5}, and \mathcal{B}(J/ψ\to π^0 e^+e^-)=(7.56\pm1.32\pm0.50)\times10^{-7}, where the first errors are statistical and the second ones systematic.

preprint2014arXiv

Precision Measurement of the Mass of the $τ$ Lepton

An energy scan near the $τ$ pair production threshold has been performed using the BESIII detector. About $24$ pb$^{-1}$ of data, distributed over four scan points, was collected. This analysis is based on $τ$ pair decays to $ee$, $eμ$, $eh$, $μμ$, $μh$, $hh$, $eρ$, $μρ$ and $πρ$ final states, where $h$ denotes a charged $π$ or $K$. The mass of the $τ$ lepton is measured from a maximum likelihood fit to the $τ$ pair production cross section data to be $m_τ = (1776.91\pm0.12 ^{+0.10}_{-0.13}$) MeV/$c^2$, which is currently the most precise value in a single measurement.

preprint2014arXiv

Precision measurements of $B(D^+ \rightarrow μ^+ ν_μ)$, the pseudoscalar decay constant $f_{D^+}$, and the quark mixing matrix element $|V_{\rm cd}|$

We report a measurement of the branching fraction $B(D^+ \rightarrow μ^+ ν_μ) = [3.71 \pm 0.19 (\rm stat) \pm 0.06 (\rm sys)]\times 10^{-4}$ based on 2.92 ${\rm fb^{-1}}$ of data accumulated at $\sqrt{s}=3.773$ GeV with the BESIII detector at the BEPCII collider. This measurement, in conjunction with the Cabibbo-Kobayashi-Maskawa matrix element $|V_{\rm cd}|$ determined from a global Standard Model fit, implies a value for the weak decay constant $f_{D^+}=(203.2 \pm 5.3 \pm 1.8)$ MeV. Additionally, using this branching fraction measurement together with a Lattice QCD prediction for $f_{D^+}$, we find $|V_{\rm cd}|=0.2210\pm 0.0058 \pm 0.0047$. In either case, these are the most precise results for these quantities to date.

preprint2014arXiv

Revealing multiple density wave orders in non-superconducting titanium oxypnictide Na$_2$Ti$_2$As$_2$O

We report an optical spectroscopy study on the single crystal of Na$_2$Ti$_2$As$_2$O, a sister compound of superconductor BaTi$_2$Sb$_2$O. The study reveals unexpectedly two density wave phase transitions. The first transition at 320 K results in the formation of a large energy gap and removes most part of the Fermi surfaces. But the compound remains metallic with residual itinerant carriers. Below 42 K, another density wave phase transition with smaller energy gap scale occurs and drives the compound into semiconducting ground state. These experiments thus enable us to shed light on the complex electronic structure in the titanium oxypnictides.

preprint2014arXiv

Search for C-parity violation in $J/ ψ\to γγ$ and $ γϕ$

Using $1.06\times10^8$ $ψ(3686)$ events recorded in $e^{+}e^{-}$ collisions at $\sqrt{s}=$ 3.686 GeV with the BESIII at the BEPCII collider, we present searches for C-parity violation in $J/ψ\to γγ$ and $ γϕ$ decays via $ψ(3686) \to J/ψπ^+π^-$. No significant signals are observed in either channel. Upper limits on the branching fractions are set to be $\mathcal{B}(J/ψ\to γγ) < 2.7 \times 10^{-7}$ and $\mathcal{B}(J/ψ\to γϕ) < 1.4 \times 10^{-6}$ at the 90\% confidence level. The former is one order of magnitude more stringent than the previous upper limit, and the latter represents the first limit on this decay channel.

preprint2014arXiv

Search for the radiative transitions $ψ(3770)\toγη_c$ and $γη_c(2S)$

By using a 2.92 fb$^{-1}$ data sample taken at $\sqrt{s} = 3.773$ GeV with the BESIII detector operating at the BEPCII collider, we search for the radiative transitions $ψ(3770)\toγη_c$ and $γη_c(2S)$ through the hadronic decays $η_c(η_c(2S))\to K^0_SK^\pmπ^\mp$. No significant excess of signal events above background is observed. We set upper limits at a 90% confidence level for the product branching fractions to be $\mathcal{B}(ψ(3770)\toγη_c)\times \mathcal{B}(η_c\to K^0_SK^\pmπ^\mp) < 1.6\times10^{-5}$ and $\mathcal{B}(ψ(3770)\toγη_c(2S))\times \mathcal{B}(η_c(2S)\to K^0_SK^\pmπ^\mp) < 5.6\times10^{-6}$. Combining our result with world-average values of $\mathcal{B}(η_c(η_c(2S))\to K^0_SK^\pmπ^\mp)$, we find the branching fractions $\mathcal{B}(ψ(3770)\toγη_c) < 6.8\times10^{-4}$ and $\mathcal{B}(ψ(3770)\toγη_c(2S)) < 2.0\times10^{-3}$ at a 90% confidence level.

preprint2014arXiv

Search for the rare decays $J/ψ\to D_s^-ρ^+$ and $J/ψ\to \overline{D}{}^0\overline{K}{}^{\ast0}$

A search for the rare decays of $J/ψ\to D_s^-ρ^+ + c.c.$ and $J/ψ\to \overline{D}{}^0\overline{K}{}^{\ast0} + c.c.$ is performed with a data sample of 225.3 million $J/ψ$ events collected with the BESIII detector. No evident signal is observed. Upper limits on the branching fractions are determined to be $\mathcal{B}(J/ψ\to D_s^-ρ^+ + c.c.) < 1.3 \times 10^{-5}$ and $\mathcal{B}(J/ψ\to \overline{D}{}^0\overline{K}{}^{\ast0} + c.c.) < 2.5 \times 10^{-6}$ at the $90\%$ confidence level.

preprint2014arXiv

Study of $e^+e^- \rightarrow p\bar{p}$ in the vicinity of $ψ(3770)$

Using 2917 $\rm{pb}^{-1}$ of data accumulated at 3.773~$\rm{GeV}$, 44.5~$\rm{pb}^{-1}$ of data accumulated at 3.65~$\rm{GeV}$ and data accumulated during a $ψ(3770)$ line-shape scan with the BESIII detector, the reaction $e^+e^-\rightarrow p\bar{p}$ is studied considering a possible interference between resonant and continuum amplitudes. The cross section of $e^+e^-\rightarrowψ(3770)\rightarrow p\bar{p}$, $σ(e^+e^-\rightarrowψ(3770)\rightarrow p\bar{p})$, is found to have two solutions, determined to be ($0.059\pm0.032\pm0.012$) pb with the phase angle $ϕ= (255.8\pm37.9\pm4.8)^\circ$ ($<$0.11 pb at the 90% confidence level), or $σ(e^+e^-\rightarrowψ(3770)\rightarrow p\bar{p}) = (2.57\pm0.12\pm0.12$) pb with $ϕ= (266.9\pm6.1\pm0.9)^\circ$ both of which agree with a destructive interference. Using the obtained cross section of $ψ(3770)\rightarrow p\bar{p}$, the cross section of $p\bar{p}\rightarrow ψ(3770)$, which is useful information for the future PANDA experiment, is estimated to be either ($9.8\pm5.7$) nb ($<17.2$ nb at 90% C.L.) or $(425.6\pm42.9)$ nb.

preprint2014arXiv

Study of J/psi->eta phi pi+pi- at BESIII

Based on a sample of $2.25\times 10^{8} J/ψ$ events taken with the BESIII detector at the BEPCII collider, we present the results of a study of the decay $J/ψ\to ηϕπ^{+}π^{-}$. The $Y(2175)$ resonance is observed in the invariant mass spectrum of $ϕf_{0}(980)$ with a statistical significance of greater than $10σ$. The corresponding mass and width are determined to be $M=2200\pm 6 \mathrm{(stat.)} \pm 5\mathrm{(syst.)} \mathrm{MeV}/c^{2}$ and $Γ=104\pm 15\mathrm{(stat.)}\pm 15\mathrm{(syst.)}$ MeV, respectively, and the product branching fraction is measured to be $\mathcal{B}(J/ψ\toηY(2175)$, $Y(2175)\to ϕf_{0}(980)$, $f_{0}(980)\to π^{+}π^{-})= (1.20\pm 0.14\mathrm{(stat.)}\pm 0.37 \mathrm{(syst.)})\times 10^{-4}$. The results are consistent within errors with those of previous experiments. We also measure the branching fraction of $J/ψ\to ϕf_1(1285)$ with $f_1(1285)\to ηπ^{+}π^{-}$ and set upper limits on the branching fractions for $J/ψ\to ϕη(1405)$/$ϕX(1835)$/$ϕX(1870)$ with $η(1405)$/$X(1835)$/$X(1870)\to ηπ^{+}π^{-}$ at the 90\% confidence level.

preprint2014arXiv

Systematic study of the symmetry energy coefficient in finite nuclei

The symmetry energy coefficients in finite nuclei have been studied systematically with a covariant density functional theory (DFT) and compared with the values calculated using several available mass tables. Due to the contamination of shell effect, the nuclear symmetry energy coefficients extracted from the binding energies have large fluctuations around the nuclei with double magic numbers. The size of this contamination is shown to be smaller for the nuclei with larger isospin value. After subtracting the shell effect with the Strutinsky method, the obtained nuclear symmetry energy coefficients with different isospin values are shown to decrease smoothly with the mass number $A$ and are subsequently fitted to the relation $\dfrac{4a_{\rm sym}}{A}=\dfrac{b_v}{A}-\dfrac{b_s}{A^{4/3}}$. The resultant volume $b_v$ and surface $b_s$ coefficients from axially deformed covariant DFT calculations are $121.73$ and $197.98$ MeV respectively. The ratio $b_s/b_v=1.63$ is in good agreement with the value derived from the previous calculations with the non-relativistic Skyrme energy functionals. The coefficients $b_v$ and $b_s$ corresponding to several available mass tables are also extracted. It is shown that there is a strong linear correlation between the volume $b_v$ and surface $b_s$ coefficients and the ratios $b_s/b_v$ are in between $1.6-2.0$ for all the cases.

preprint2013arXiv

LSS-GAC - A LAMOST Spectroscopic Survey of the Galactic Anti-center

As a major component of the LAMOST Galactic surveys, the LAMOST Spectroscopic Survey of the Galactic Anti-center (LSS-GAC) will survey a significant volume of the Galactic thin/thick disks and halo in a contiguous sky area of ~ 3,400sq.deg., centered on the Galactic anti-center (|b| <= 30°, 150 <= l <= 210°), and obtain λλ3800--9000 low resolution (R ~ 1,800) spectra for a statistically complete sample of >= 3M stars of all colors, uniformly and randomly selected from (r, g - r) and (r, r - i) Hess diagrams obtained from a CCD imaging photometric survey of ~ 5,400sq.deg. with the Xuyi 1.04/1.20 m Schmidt Telescope, ranging from r = 14.0 to a limiting magnitude of r = 17.8 (18.5 for limited fields). The survey will deliver spectral classification, radial velocity Vr and stellar parameters (effective temperature Teff, surface gravity log g and metallicity [Fe/H]) for millions of Galactic stars. Together with Gaia which will provide accurate distances and tangential velocities for a billion stars, the LSS-GAC will yield a unique dataset to study the stellar populations, chemical composition, kinematics and structure of the disks and their interface with the halo, identify streams of debris of tidally disrupted dwarf galaxies and clusters, probe the gravitational potential and dark matter distribution, map the 3D distribution of interstellar dust extinction, search for rare objects (e.g. extremely metal-poor or hyper-velocity stars), and ultimately advance our understanding of the assemblage of the Milky Way and other galaxies and the origin of regularity and diversity of their properties. ... (abridged)

preprint2013arXiv

Measurement of $η^\prime\rightarrow π^+π^- e^+e^-$ and $η^\prime\rightarrow π^+π^-μ^+μ^-$

Based on a sample of 225.3 million J/ψevents accumulated with the BESIII detector at the BEPCII, the decays of η' to pi+pi-l+l- are studied via J/ψto γη'. A clear η' signal is observed in the pi+pi-e+e- mass spectrum, and the branching fraction is measured to be \BR(η' to pi+pi-e+e-) = (2.11\pm0.12 (stat.)\pm0.15 (syst.))\times10^{-3}, which is in good agreement with theoretical predictions and the previous measurement, but is determined with much higher precision. No η' signal is found in the pi+ pi- mu+ mu- mass spectrum, and the upper limit is determined to be \BR(η' to pi+ pi- mu+ mu-)<2.9\times10^{-5} at the 90% confidence level.

preprint2013arXiv

Measurements of Baryon Pair Decays of $χ_{cJ}$ Mesons

Using 106 $\times 10^{6}$ $ψ^{\prime}$ decays collected with the BESIII detector at the BEPCII, three decays of $χ_{cJ}$ ($J=0,1,2$) with baryon pairs ($\llb$, $\ssb$, $\SSB$) in the final state have been studied. The branching fractions are measured to be $\cal{B}$$(χ_{c0,1,2}\rightarrowΛ\barΛ) =(33.3 \pm 2.0 \pm 2.6)\times 10^{-5}$, $(12.2 \pm 1.1 \pm 1.1)\times 10^{-5}$, $(20.8 \pm 1.6 \pm 2.3)\times 10^{-5}$; $\cal{B}$$(χ_{c0,1,2}\rightarrowΣ^{0}\barΣ^{0})$ = $(47.8 \pm 3.4 \pm 3.9)\times 10^{-5}$, $(3.8 \pm 1.0 \pm 0.5)\times 10^{-5}$, $(4.0 \pm 1.1 \pm 0.5) \times 10^{-5}$; and $\cal{B}$$(χ_{c0,1,2}\rightarrowΣ^{+}\barΣ^{-})$ = $(45.4 \pm 4.2 \pm 3.0)\times 10^{-5}$, $(5.4 \pm 1.5 \pm 0.5)\times 10^{-5}$, $(4.9 \pm 1.9 \pm 0.7)\times 10^{-5}$, where the first error is statistical and the second is systematic. Upper limits on the branching fractions for the decays of $χ_{c1,2}\rightarrowΣ^{0}\barΣ^{0}$, $Σ^{+}\barΣ^{-}$, are estimated to be $\cal{B}$$(χ_{c1}\rightarrowΣ^{0}\barΣ^{0}) < 6.2\times 10^{-5}$, $\cal{B}$$(χ_{c2}\rightarrowΣ^{0}\barΣ^{0}) < 6.5\times 10^{-5}$, $\cal{B}$$(χ_{c1}\rightarrowΣ^{+}\barΣ^{-}) < 8.7\times 10^{-5}$ and $\cal{B}$$(χ_{c2}\rightarrowΣ^{+}\barΣ^{-}) < 8.8\times 10^{-5}$ at the 90% confidence level.

preprint2013arXiv

Modeling the linewidth dependence of coherent terahertz emission from intrinsic Josephson junction stacks in the hot-spot regime

Recently it has been found that, when operated at large input power, the linewidth of terahertz radiation emitted from intrinsic Josephson junction stacks can be as narrow as some megahertz. In this high-bias regime a hot spot coexists with regions which are still superconducting. Surprisingly, the linewidth was found to decrease with increasing bath temperature. We present a simple model describing the dynamics of the stack in the presence of a hot spot by two parallel arrays of pointlike Josephson junctions and an additional shunt resistor in parallel. Heat diffusion is taken into account by thermally coupling all elements to a bath at temperature T_b. We present current-voltage characteristics of the coupled system and calculations of the linewidth of the radiation as a function of T_b. In the presence of a spatial gradient of the junction parameters critical current and resistance, the linewidth deceases with increasing T_b, similar to the experimental observation.

preprint2013arXiv

Nonlinear stability of flock solutions in second-order swarming models

In this paper we consider interacting particle systems which are frequently used to model collective behavior in animal swarms and other applications. We study the stability of orientationally aligned formations called flock solutions, one of the typical patterns emerging from such dynamics. We provide an analysis showing that the nonlinear stability of flocks in second-order models entirely depends on the linear stability of the first-order aggregation equation. Flocks are shown to be nonlinearly stable as a family of states under reasonable assumptions on the interaction potential. Furthermore, we numerically verify that commonly used potentials satisfy these hypotheses and investigate the nonlinear stability of flocks by an extensive case-study of uniform perturbations.

preprint2013arXiv

Observation of a charged charmoniumlike structure in e+e- to pi+pi-J/psi at \sqrt{s}=4.26 GeV

We study the process e+e- to pi+pi-J/psi at a center-of-mass energy of 4.260 GeV using a 525 pb^{-1} data sample collected with the BESIII detector operating at the Beijing Electron Positron Collider. The Born cross section is measured to be (62.9\pm 1.9\pm 3.7) pb, consistent with the production of the Y(4260). We observe a structure at around 3.9 GeV/c^2 in the π^\pm J/psi mass spectrum, which we refer to as the Z_c(3900). If interpreted as a new particle, it is unusual in that it carries an electric charge and couples to charmonium. A fit to the π^\pm J/psi invariant mass spectrum, neglecting interference, results in a mass of (3899.0\pm 3.6\pm 4.9) MeV/c^2 and a width of (46\pm 10\pm 20) MeV. Its production ratio is measured to be R=\frac{σ(e+e- to π^\pm Z_c(3900)^\mp to pi+pi-J/psi))} {σ(e+e- to pi+pi-J/psi)}=(21.5\pm 3.3\pm 7.5)%. In all measurements the first errors are statistical and the second are systematic.

preprint2013arXiv

Observation of a charged charmoniumlike structure Z_c(4020) and search for the Z_c(3900) in e+e- to pi+pi-h_c

We study e+e- --> pi+pi-h_c at center-of-mass energies from 3.90 GeV to 4.42 GeV using data samples collected with the BESIII detector operating at the Beijing Electron Positron Collider. The Born cross sections are measured at 13 energies, and are found to be of the same order of magnitude as those of e+e- --> pi+pi-J/psi but with a different line shape. In the π^\pm h_c mass spectrum, a distinct structure, referred to as Z_c(4020), is observed at 4.02 GeV/c^2. The Z_c(4020) carries an electric charge and couples to charmonium. A fit to the π^\pm h_c invariant mass spectrum, neglecting possible interferences, results in a mass of (4022.9\pm 0.8\pm 2.7) MeV/c^2 and a width of (7.9\pm 2.7\pm 2.6) MeV for the Z_c(4020), where the first errors are statistical and the second systematic. No significant Z_c(3900) signal is observed, and upper limits on the Z_c(3900) production cross sections in π^\pm h_c at center-of-mass energies of 4.23 and 4.26 GeV are set.

preprint2013arXiv

Observation of a structure at 1.84 GeV/c$^2$ in the $3(π^+π^-)$ mass spectrum in $J/ψ\rightarrow γ3(π^+π^-)$ decays

With a sample of 225.3 million $J/ψ$ events taken with the BESIII detector, the decay $J/ψ\rightarrow γ3(π^+π^-)$ is analyzed. A structure at 1.84 GeV/c$^2$ is observed in the $3(π^+π^-)$ invariant mass spectrum with a statistical significance of 7.6$σ$ . The mass and width are measured to be $M=1842.2\pm 4.2^{+7.1}_{-2.6}$ MeV/c$^2$ and $Γ=83\pm 14 \pm 11$ MeV. The product branching fraction is determined to be $B(J/ψ\rightarrowγX(1840))\times B(X(1840)\rightarrow 3(π^+π^-))=(2.44\pm0.36^{+0.60}_{-0.74})\times 10^{-5}$. No $η^\prime$ signals are observed in the $3(π^+π^-)$ invariant mass spectrum, and the upper limit of the branching fraction for the decay $η^\prime\rightarrow 3(π^+π^-)$ is set to be $3.1\times10^{-5}$ at a 90% confidence level.

preprint2013arXiv

Partial wave analysis of $J/ψ\to γηη$

Based on a sample of $2.25\times 10^{8}$ $J/ψ$ events collected with the BESIII detector at BEPCII, a full partial wave analysis on $J/ψ\toγηη$ was performed using the relativistic covariant tensor amplitude method. The results show that the dominant $0^{++}$ and $2^{++}$ components are from the $f_0(1710)$, $f_0(2100)$, $f_0(1500)$, $f_2'(1525)$, $f_2(1810)$ and $f_2(2340)$. The resonance parameters and branching fractions are also presented.

preprint2013arXiv

Partial wave analysis of $ψ(2S) \to p \bar{p}η$

Using a sample of $1.06 \times 10^{8}$ $ψ(2S)$ events collected with the BESIII detector at BEPCII, the decay $ψ(2S) \to p \bar{p}η$ is studied. A partial wave analysis determines that the intermediate state N(1535) with a mass of $1524\pm5^{+10}_{-4}$ MeV/$c^2$ and a width of $130^{+27+57}_{-24-10}$ MeV/$c^2$ is dominant in the decay; the product branching fraction is determined to be $B(ψ(2S) \to N(1535)\bar{p})\times B(N(1535)\to pη)+c.c. = (5.2\pm0.3^{+3.2}_{-1.2})\times 10^{-5}$. Furthermore, the branching fraction of $ψ(2S) \to ηp \bar{p}$ is measured to be $(6.4\pm0.2\pm0.6)\times 10^{-5}$.

preprint2013arXiv

Precision Measurements of ${\cal B}[ψ(3686) \to π^+π^- J/ψ]$ and ${\cal B}[J/ψ\to l^+l^-]$

Based on $(106.41 \pm 0.86)\times 10^{6}$ $ψ(3686)$ events collected with the BESIII detector at the BEPCII collider, the branching fractions of $ψ(3686) \to π^+π^- J/ψ$, $J/ψ\to e^+e^- $, and $J/ψ\to μ^+μ^-$ are measured. We obtain ${\cal B}[ψ(3686) \to π^+π^-J/ψ]=(34.98\pm 0.02\pm 0.45)%$, ${\cal B}[J/ψ\to e^+e^-] = (5.983 \pm 0.007 \pm 0.037)%$ and ${\cal B}[J/ψ\to μ^+μ^-] = (5.973 \pm 0.007 \pm 0.038)%$. The measurement of ${\cal B}[ψ(3686) \to π^{+}π^{-}J/ψ]$ confirms the CLEO-c measurement, and is apparently larger than the others. The measured $J/ψ$ leptonic decay branching fractions agree with previous experiments within one standard deviation. These results lead to ${\cal B}[J/ψ\to l^+l^-] = (5.978 \pm 0.005 \pm 0.040)%$ by averaging over the $e^{+}e^{-}$ and $μ^{+}μ^{-}$ channels and a ratio of ${\cal B}[J/ψ\to e^+e^-] / {\cal B}[J/ψ\to μ^+μ^-] = 1.0017 \pm 0.0017 \pm 0.0033$, which tests $e$-$μ$ universality at the four tenths of a percent level. All the measurements presented in this paper are the most precise in the world to date.

preprint2013arXiv

Scaling of electrical and thermal conductivities in an almost integrable chain

Many low-dimensional materials are well described by integrable one-dimensional models such as the Hubbard model of electrons or the Heisenberg model of spins. However, the small perturbations to these models required to describe real materials are expected to have singular effects on transport quantities: integrable models often support dissipationless transport, while weak non-integrable terms lead to finite conductivities. We use matrix-product-state methods to obtain quantitative values of spin/electrical and thermal conductivities in an almost integrable gapless chain (an XXZ spin chain with staggered fields, or equivalently a spinless fermion chain with staggered on-site potentials). The results at low temperatures validate a scaling theory based on bosonization.

preprint2013arXiv

Search for $η$ and $η^\prime\to π^+ e^- \barν_e +c.c.$ decays in $\jpsi \to ϕη$ and $ϕη^\prime$

Using a sample of 225.3 million $\jpsi$ events collected with the BESIII detector at the BEPCII $e^+e^-$ collider in 2009, searches for the decays of $η$ and $η^\prime\toπ^+ e^- \barν_e +c.c.$ in $\jpsi \to ϕη$ and $ϕη^\prime$ are performed. The $ϕ$ signals, which are reconstructed in $K^+K^-$ final states, are used to tag $η$ and $η^\prime$ semileptonic decays. No signals are observed for either $η$ or $η^\prime$, and upper limits at the 90% confidence level are determined to be $7.3\times 10^{-4}$ and $5.0\times 10^{-4}$ for the ratios $\frac{{\mathcal B}(η\to π^+ e^- \barν_e +c.c.)}{{\mathcal B}(η\to \pip\pim\piz)}$ and $\frac{{\mathcal B}(η^\prime\to π^+ e^-\barν_e +c.c.)}{{\mathcal B}(η^\prime \to \pip\pimη)}$, respectively. These are the first upper limit values determined for $η$ and $η^\prime$ semileptonic weak decays.

preprint2013arXiv

Search for Baryonic Decays of ψ(3770) and ψ(4040)

By analyzing data samples of 2.9 fb^{-1} collected at \sqrt s=3.773 GeV, 482 pb^{-1} collected at \sqrt s=4.009 GeV and 67 pb^{-1} collected at \sqrt s=3.542, 3.554, 3.561, 3.600 and 3.650 GeV with the BESIII detector at the BEPCII storage ring, we search for ψ(3770) and ψ(4040) decay to baryonic final states, including Λ\barΛπ^+π^-, Λ\barΛπ^0, Λ\barΛη, Σ^+ \barΣ^-, Σ^0 \barΣ^0, Ξ^-\barΞ^+ and Ξ^0\barΞ^0 decays. None are observed, and upper limits are set at the 90% confidence level.

preprint2013arXiv

Search for the Lepton Flavor Violation Process $J/ψ\to eμ$ at BESIII

We search for the lepton-flavor-violating decay of the $J/ψ$ into an electron and a muon using $(225.3\pm2.8)\times 10^{6}$ $J/ψ$ events collected with the BESIII detector at the BEPCII collider. Four candidate events are found in the signal region, consistent with background expectations. An upper limit on the branching fraction of $\mathcal{B}(J/ψ\to eμ)< 1.5 \times 10^{-7}$ (90% C.L.) is obtained.

preprint2013arXiv

Study of $J/ψ\rightarrow ωp \bar{p}$ at BESIII

The decay $J/ψ\rightarrow ωp \bar{p}$ has been studied, using $225.3\times 10^{6}$ $J/ψ$ events accumulated at BESIII. No significant enhancement near the $p\bar{p}$ invariant-mass threshold (denoted as $X(p\bar{p})$) is observed. The upper limit of the branching fraction $\mathcal{B}(J/ψ\rightarrow ωX(p\bar{p}) \rightarrow ωp \bar{p})$ is determined to be $3.9\times10^{-6}$ at the 95% confidence level. The branching fraction of $J/ψ\rightarrow ωp \bar{p}$ is measured to be $\mathcal{B}(J/ψ\rightarrow ωp \bar{p}) =(9.0 \pm 0.2\ (\text{stat.})\pm 0.9\ (\text{syst.})) \times 10^{-4}$.

preprint2013arXiv

Surface and Bulk Electronic Structure of the Strongly Correlated System SmB$_{6}$ and Implications for a Topological Kondo Insulator

Recent theoretical calculations and experimental results suggest that the strongly correlated material SmB$_{6}$ may be a realization of a topological Kondo insulator. We have performed an angle-resolved photoemission spectroscopy study on SmB$_{6}$ in order to elucidate elements of the electronic structure relevant to the possible occurrence of a topological Kondo insulator state. The obtained electronic structure in the whole three-dimensional momentum space reveals one electron-like 5d bulk band centred at the X point of the bulk Brillouin zone that is hybridized with strongly correlated f electrons, as well as the opening of a Kondo bandgap ($Δ_B$ $\sim$ 20 meV) at low temperature. In addition, we observe electron-like bands forming three Fermi surfaces at the center $\barΓ$ point and boundary $\bar{X}$ point of the surface Brillouin zone. These bands are not expected from calculations of the bulk electronic structure, and their observed dispersion characteristics are consistent with surface states. Our results suggest that the unusual low-temperature transport behavior of SmB$_{6}$ is likely to be related to the pronounced surface states sitting inside the band hybridisation gap and/or the presence of a topological Kondo insulating state.

preprint2012arXiv

Beaming Effect in Fermi Blazars

The γ-ray loud blazars (flat spectrum radio quasars--FSRQs and BL Lacertae objects-BLs) are very bright in the γ-ray bands, which is perhaps associated with a beaming effect. Therefore, one can expect that the γ-ray luminosity is correlated with the beaming factor. In this paper, we investigated the relation between the radio Doppler factors and the gamma-ray luminosities. Our analysis suggests that the γ-ray luminosity be strongly correlated with the factor of δ_R for the whole sample, FSRQs, and BLs. When the effect of a common redshift is excluded, the correlation still exists for the FSRQs sub-sample suggesting that the γ-rays are strongly beamed. However, the partial correlation analysis does not show a correlation for the small BL Lac sample.

preprint2012arXiv

Effect of disorder in the charge-density-wave compounds LaTe$_{1.95}$ and CeTe$_{1.95-x}$Se$_x$ (x=0 and 0.16) as revealed by optical spectroscopy

We present optical spectroscopy measurements on rare-earth ditelluride single crystals of LaTe$_{1.95}$ and CeTe$_{1.95-x}$Se$_x$ (x=0 and 0.16). The measurements reveal formation of charge density wave energy gaps at rather high energy levels, e.g. 2$Δ\sim$ 8500 \cm for LaTe$_{1.95}$, and 6800 \cm for CeTe$_{1.95}$. More strikingly, the study reveals that, different from the rare-earth tri-tellurides, the Te vacancies and disorder effect play a key role in the low-energy charge excitations of ditelluride systems. Although an eminent peak is observed between 800 and 1500 \cm in conductivity spectra for LaTe$_{1.95}$, and CeTe$_{1.95-x}$Se$_x$ (x=0. 0.16), our analysis indicates that it could not be attributed to the formation of a small energy gap, instead it could be well accounted for by the localization modified Drude model. Our study also indicates that the low-tempreature optical spectroscopic features are distinctly different from a semiconducting CDW state with entirely gapped Fermi surfaces.

preprint2012arXiv

Study of the near-threshold $ωϕ$ mass enhancement in doubly OZI suppressed $J/ψ\to γωϕ$ decays

A 2.25$\times10^8$ $\jpsi$ event sample accumulated with the BESIII detector is used to study the doubly OZI suppressed decay modes $\jpsi\toγ\of$, $ω\to\ppp$, $ϕ\to\kk$. A strong deviation ($>$ 30$σ$) from three-body $\jpsi\toγωϕ$ phase space is observed near the $ωϕ$ mass threshold that is consistent with a previous observation reported by the BESII experiment. A partial wave analysis with a tensor covariant amplitude that assumes that the enhancement is due to the presence of a resonance, the X(1810), is performed, and confirms that the spin-parity of the X(1810) is $0^{++}$. The mass and width of the X(1810) are determined to be $M=1795\pm7$(stat)$^{+13}_{-5}$(syst)$\pm$19(mod) MeV/$c^2$ and $Γ=95\pm10$(stat)$^{+21}_{-34}$(syst)$\pm$75(mod) MeV/$c^2$, respectively, and the product branching fraction is measured to be ${\cal B}(\jpsi\toγX(1810))\times{\cal B}(X(1810)\to\of)=(2.00\pm0.08$(stat)$^{+0.45}_{-1.00}$(syst)$\pm$1.30(mod))$\times10^{-4}$. %where the first error is the statistical error and the second the systematical error. These results are consistent within errors with those of the BESII experiment.

preprint2012arXiv

Surface adatom conductance filtering in scanning tunneling spectroscopy of Co-doped BaFe2As2 iron pnictide superconductors

We establish in a combination of ab initio theory and experiments that the tunneling process in scanning tunneling microscopy/spectroscopy on the A-122 iron pnictide superconductors - in this case BaFe$_{2-x}$Co$_x$As$_2$ - involve a strong adatom filtering of the differential conductance from the near-EF Fe3d states, which in turn originates from the top-most sub-surface Fe layer of the crystal. The calculations show that the dominance of surface Ba-related tunneling pathways leaves fingerprints found in the experimental differential conductance data, including large particle-hole asymmetry and an energy-dependent contrast inversion.

preprint2011arXiv

Equilibria of biological aggregations with nonlocal repulsive-attractive interactions

We consider the aggregation equation $ρ_{t}-\nabla\cdot(ρ\nabla K\astρ) =0$ in $\mathbb{R}^{n}$, where the interaction potential $K$ incorporates short-range Newtonian repulsion and long-range power-law attraction. We study the global well-posedness of solutions and investigate analytically and numerically the equilibrium solutions. We show that there exist unique equilibria supported on a ball of $\mathbb{R}^n$. By using the method of moving planes we prove that such equilibria are radially symmetric and monotone in the radial coordinate. We perform asymptotic studies for the limiting cases when the exponent of the power-law attraction approaches infinity and a Newtonian singularity, respectively. Numerical simulations suggest that equilibria studied here are global attractors for the dynamics of the aggregation model.

preprint2011arXiv

Hertz-level Measurement of the 40Ca+ 4s 2S1/2-3d 2D5/2 Clock Transition Frequency With Respect to the SI Second through GPS

We report a frequency measurement of the clock transition of a single ^40Ca^+ ion trapped and laser cooled in a miniature ring Paul trap with 10^-15 level uncertainty. In the measurement, we used an optical frequency comb referenced to a Hydrogen maser, which was calibrated to the SI second through the Global Positioning System (GPS). Two rounds of measurements were taken in May and June 2011, respectively. The frequency was measured to be 411 042 129 776 393.0(1.6) Hz with a fractional uncertainty of 3.9{\times}10^-15 in a total averaging time of > 2{\times}10^6 s within 32 days.

preprint2010arXiv

Mini review of Poincaré invariant quantum theory

We review the construction and applications of exactly Poincaré invariant quantum mechanical models of few-degree of freedom systems. We discuss the construction of dynamical representations of the Poincaré group on few-particle Hilbert spaces, the relation to quantum field theory, the formulation of cluster properties, and practical considerations related to the construction of realistic interactions and the solution of the dynamical equations. Selected applications illustrate the utility of this approach.

preprint2009arXiv

A high resolution, hard x-ray photoemission investigation of La_(2-2x)Sr_(1+2x)Mn_2O_7 (0.30<x<0.50): on microscopic phase separation and the surface electronic structure of a bilayered CMR manganite

Photoemission data taken with hard x-ray radiation on cleaved single crystals of the bilayered, colossal magnetoresistant manganite La_(2-2x)Sr_(1+2x)Mn_2O_7 (LSMO) with 0.30<x<0.50 are presented. Making use of the increased bulk-sensitivity upon hard x-ray excitation it is shown that the core level footprint of the electronic structure of the LSMO cleavage surface is identical to that of the bulk. Furthermore, by comparing the core level shift of the different elements as a function of doping level x, it is shown that microscopic phase separation is unlikely to occur for this particular manganite well above the Curie temperature.

preprint2009arXiv

Electronic structure studies of BaFe2As2 by angle-resolved photoemission spectroscopy

We report high resolution angle-resolved photoemission spectroscopy (ARPES) studies of the electronic structure of BaFe$_2$As$_2$, which is one of the parent compounds of the Fe-pnictide superconductors. ARPES measurements have been performed at 20 K and 300 K, corresponding to the orthorhombic antiferromagnetic phase and the tetragonal paramagnetic phase, respectively. Photon energies between 30 and 175 eV and polarizations parallel and perpendicular to the scattering plane have been used. Measurements of the Fermi surface yield two hole pockets at the $Γ$-point and an electron pocket at each of the X-points. The topology of the pockets has been concluded from the dispersion of the spectral weight as a function of binding energy. Changes in the spectral weight at the Fermi level upon variation of the polarization of the incident photons yield important information on the orbital character of the states near the Fermi level. No differences in the electronic structure between 20 and 300 K could be resolved. The results are compared with density functional theory band structure calculations for the tetragonal paramagnetic phase.

preprint2009arXiv

Optical properties of BaFe$_{2-x}$Co$_x$As$_2$

We present detailed temperature dependent optical data on BaFe$_{2-x}$Co$_{x}$As$_{2}$ (BCFA), with x = 0.14, between 4 meV and 6.5 eV. We analyze our spectra to determine the main optical parameters and show that in this material the interband conductivity already starts around 10 meV. We determine the superfluid density to be 2.2 10^{7}$ cm^{-2}, which places optimally doped BFCA close to the Uemura line. Our experimental data shows clear signs of a superconducting gap with 2$Δ_{1}$ = 6.2 $\pm$ 0.8 meV. In addition we show that the optical spectra are consistent with the presence of an additional band of strongly scattered carriers with a larger gap, 2$Δ_{2}$ = 14 $\pm$ 2 meV.

preprint2009arXiv

Orbital character variation of the Fermi surface and doping dependent changes of the dimensionality in BaFe2-xCoxAs2 from angle-resolved photoemission spectroscopy

From a combination of high resolution angle-resolved photoemission spectroscopy and density functional calculations, we show that BaFe2As2 possesses essentially two-dimensional electronic states, with a strong change of orbital character of two of the Gamma-centered Fermi surfaces as a function of kz. Upon Co doping, the electronic states in the vicinity of the Fermi level take on increasingly three-dimensional character. Both the orbital variation with kz and the more three-dimensional nature of the doped compounds have important consequences for the nesting conditions and thus possibly also for the appearance of antiferromagnetic and superconducting phases.

preprint2008arXiv

Critical Scaling of the Magnetization and Magnetostriction in the Weak Itinerant Ferromagnet UIr

The weak itinerant ferromagnet UIr is studied by magnetization and magnetostriction measurements. Critical behavior, which surprisingly extends up to several Tesla, is observed at the Curie temperature $T_C\simeq45$ K and is analyzed using Arrott and Maxwell relations. Critical exponents are found that do not match with any of the well-known universality classes. The low-temperature magnetization $M_s\simeq0.5$ $μ_B \cong const.$ below 3 T rises towards higher fields and converges asymptotically around 50 T with the magnetization at $T_C$. From the magnetostriction and magnetization data, we extract the uniaxial pressure dependences of $T_C$, using a new method presented here, and of $M_s$. These results should serve as a basis for understanding spin fluctuations in anisotropic itinerant ferromagnets.

preprint2006arXiv

Optical Photometrical Observations and Variability for Quasar 4C 29.45

We reported the result of long term optical variability of the blazar 4C 29.45 (QSO 1156+295, Ton 599), carried out optical photometric observations in Johnson V, Cousins RI passbands during April 1997 to March 2002 using the 1.56 meter telescope of the Shanghai Astronomical Observatory (SHAO) at Sheshan, China, compiled the post-1974 optical photometric data of the blazar by combining our new observations with the published optical data, and found maximum variations in different passbands: $Δ$U = 4.41 mag, $Δ$B = 5.55 mag, $Δ$V = 4.53 mag, $Δ$R = 5.80 mag, and $Δ$I = 5.34 mag. The average color indices are: U$-$B = $-0.54\pm$0.18 mag, B$-$V = 0.56$\pm$0.21 mag, B$-$R = 0.93$\pm$0.18 mag, B$-$I = 1.51$\pm$0.24 mag, V$-$R = 0.44$\pm$0.15 mag and V$-$I = 1.03$\pm$0.23 mag. The post-1974 data give us an excellent opportunity to search for the existence of possible periodicity in the light curve. In search for periodicity in the R passband light curve, we performed Jurkevich test and power spectral (Fourier) analysis methods, and CLEANest algorithms to remove false signals. We found possible periods of 3.55 and 1.58 years. The possible mechanism for the periodic variability was discussed.

preprint2005arXiv

Search for First-Generation Scalar Leptoquarks in $\bm{p \bar{p}}$ collisions at $\sqrt{s}$=1.96 TeV

We report on a search for pair production of first-generation scalar leptoquarks ($LQ$) in $p \bar{p}$ collisions at $\sqrt{s}$=1.96 TeV using an integrated luminosity of 203 $pb^{-1}$ collected at the Fermilab Tevatron collider by the CDF experiment. We observe no evidence for $LQ$ production in the topologies arising from $LQ \bar{LQ} \to eqeq$ and $LQ \bar{LQ} \to eq νq$, and derive 95% C.L. upper limits on the $LQ$ production cross section. %as a function of $β$, where $β$ is the branching fraction for $LQ \to eq$. The results are combined with those obtained from a separately reported CDF search in the topology arising from $LQ\bar{LQ} \to νq νq$ and 95% C.L. lower limits on the LQ mass as a function of $β= BR(LQ \to eq) $ are derived. The limits are 236, 205 and 145 GeV/c$^2$ for $β$ = 1, $β$ = 0.5 and $β$ = 0.1, respectively.

preprint2000arXiv

Artifactual log-periodicity in finite size data: Relevance for earthquake aftershocks

The recently proposed discrete scale invariance and its associated log-periodicity are an elaboration of the concept of scale invariance in which the system is scale invariant only under powers of specific values of the magnification factor. We report on the discovery of a novel mechanism for such log-periodicity relying solely on the manipulation of data. This ``synthetic'' scenario for log-periodicity relies on two steps: (1) the fact that approximately logarithmic sampling in time corresponds to uniform sampling in the logarithm of time; and (2) a low-pass-filtering step, as occurs in constructing cumulative functions, in maximum likelihood estimations, and in de-trending, reddens the noise and, in a finite sample, creates a maximum in the spectrum leading to a most probable frequency in the logarithm of time. We explore in detail this mechanism and present extensive numerical simulations. We use this insight to analyze the 27 best aftershock sequences studied by Kisslinger and Jones [1991] to search for traces of genuine log-periodic corrections to Omori's law, which states that the earthquake rate decays approximately as the inverse of the time since the last main shock. The observed log-periodicity is shown to almost entirely result from the ``synthetic scenario'' owing to the data analysis. From a statistical point of view, resolving the issue of the possible existence of log-periodicity in aftershocks will be very difficult as Omori's law describes a point process with a uniform sampling in the logarithm of the time. By construction, strong log-periodic fluctuations are thus created by this logarithmic sampling.

preprint2000arXiv

Re-examination of log-periodicity observed in the seismic precursors of the 1989 Loma Prieta earthquake

Based on several empirical evidence, a series of papers has advocated the concept that seismicity prior to a large earthquake can be understood in terms of the statistical physics of a critical phase transition. In this model, the cumulative Benioff strain (BS) increases as a power-law time-to-failure before the final event. This power law reflects a kind of scale invariance with respect to the distance to the critical point. A few years ago, on the basis of a fit of the cumulative BS released prior to the 1989 Loma Prieta earthquake, Sornette and Sammis [1995] proposed that this scale invariance could be partially broken into a discrete scale invariance (DSI). The observable consequence of DSI takes the form of log-periodic oscillations decorating the accelerating power law. They found that the quality of the fit and the predicted time of the event are significantly improved by the introduction of log-periodicity. Here, we present a battery of synthetic tests performed to quantify the statistical significance of this claim. We find that log-periodic oscillations with frequency and regularity similar to those of the Loma Prieta case are very likely to be generated by the interplay of the low pass filtering step due to the construction of cumulative functions together with the approximate power law acceleration. Thus, the single Loma Prieta case alone cannot support the initial claim and additional cases and further study are needed to increase the signal-to-noise ratio if any. The present study will be a useful methodological benchmark for future testing of additional events when the methodology and data to construct reliable Benioff strain function become available.