Researcher profile

J. Liu

J. Liu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

hep-ex physics.ins-det nucl-ex cond-mat.str-el physics.optics physics.plasm-ph astro-ph.SR cond-mat.mtrl-sci Distributed, Parallel, and Cluster Computing Genomics Methodology Multiagent Systems nlin.PS physics.atom-ph physics.geo-ph physics.space-ph quant-ph Systems and Control

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author

90works

0followers

18topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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

Building this graph slice

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

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

Discovery of a new $γ$-ray source LHAASO J0341+5258 with emission up to 200TeV

We report the discovery of a new unidentified extended $γ$-ray source in the Galactic plane named LHAASO J0341+5258 with a pre-trial significance of 8.2 standard deviations above 25 TeV. The best fit position is R.A.$=55.34^{\circ}\pm0.11^{\circ}$ and Dec$=52.97^{\circ}\pm0.07^{\circ}$. The angular size of LHAASO J0341+5258 is $0.29^\circ \pm 0.06^\circ_{stat} \pm0.02^\circ_{sys}$. The flux above 25 TeV is about $20\%$ of the flux of Crab Nebula. Although a power-law fit of the spectrum from 10 TeV to 200 TeV with the photon index $α=2.98 \pm 0.19_{stat} \pm 0.02_{sys}$ is not excluded, the LHAASO data together with the flux upper limit at 10 GeV set by the Fermi LAT observation, indicate a noticeable steepening of an initially hard power-law spectrum %($α\leq 1.75$) spectrum with a cutoff at $\approx 50$ TeV. We briefly discuss the origin of UHE gamma-rays. The lack of an energetic pulsar and a young SNR inside or in the vicinity of LHAASO J0341+5258 challenge, but do not exclude both the leptonic and hadronic scenarios of gamma-ray production.

preprint2026arXiv

Discovery of the Ultra-high energy gamma-ray source LHAASO J2108+5157

We report the discovery of a UHE gamma-ray source, LHAASO J2108+5157, by analyzing the LHAASO-KM2A data of 308.33 live days. Significant excess of gamma-ray induced showers is observed in both energy bands of 25-100 TeV and $\gt$100 TeV with 9.5 sigma and 8.5 sigma, respectively. This source is not significantly favored as an extensive source with the angular extension smaller than the point-spread function of KM2A. The measured energy spectrum from 20 to 200 TeV can be approximately described by a power-law function with an index of -2.83$\pm$ 0.18stat. A harder spectrum is demanded at lower energies considering the flux upper limit set by Fermi-LAT observations. The position of the gamma-ray emission is correlated with a giant molecular cloud, which favors a hadronic origin. No obvious counterparts have been found, deeper multiwavelength observations will help to shed new light on this intriguing UHE source.

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

Exploring Lorentz Invariance Violation from Ultra-high-energy Gamma Rays Observed by LHAASO

Recently the LHAASO Collaboration published the detection of 12 ultra-high-energy gamma-ray sources above 100 TeV, with the highest energy photon reaching 1.4 PeV. The first detection of PeV gamma rays from astrophysical sources may provide a very sensitive probe of the effect of the Lorentz invariance violation (LIV), which results in decay of high-energy gamma rays in the superluminal scenario and hence a sharp cutoff of the energy spectrum. Two highest energy sources are studied in this work. No signature of the existence of LIV is found in their energy spectra, and the lower limits on the LIV energy scale are derived. Our results show that the first-order LIV energy scale should be higher than about 10^5 times the Planck scale M_{pl} and that the second-order LIV scale is >10^{-3}M_{pl}. Both limits improve by at least one order of magnitude the previous results.

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

Multimode Objective Lens for Momentum Microscopy and XPEEM: Experiments

A new type of objective lens has recently been proposed for use in X-ray photoemission electron microscopes (XPEEMs) and momentum microscopes. Adding a ring electrode concentric with the extractor allows the field in the gap between the sample and the extractor to be shaped. Forming a lens field in this gap reduces the field strength at the sample by up to an order of magnitude. This mitigates the risk of field emission, particularly for cleaved samples with sharp edges. A retarding field can redirect all slow electrons, thus eliminating the primary contribution to the space-charge interaction. Here we present the first experimental investigation of the new lens, examining its performance at photon energies ranging from the extreme ultraviolet produced by a high-harmonic generation (HHG)-based source to soft and hard X-rays at two synchrotron facilities. The gap lens in a region without electrodes enables large working distances up to 23 mm. Reduced aberrations allow for larger fields of view in both k-space and real-space imaging, with resolutions comparable to those of conventional cathode lenses. However, field strengths are an order of magnitude smaller. The zero-field mode enables the study of 3D structured objects and is therefore beneficial for small cleaved samples as well as for operando devices involving top electrodes. The repeller mode reduces space-charge effects, but results in a smaller k-field diameter. This reduction ranges from 10% at hard X-ray energies to 50% in the XUV range. The usable energy interval is also reduced by a factor of two. In time-of-flight XPEEM mode the raw data show a resolution of 250 nm, which can be improved to better than 100 nm through data processing.

preprint2026arXiv

Reconstruction of atmospheric neutrinos in DUNE's horizontal-drift far-detector module

This paper reports on the capabilities in reconstructing and identifying atmospheric neutrino interactions in one of the Deep Underground Neutrino Experiment's (DUNE) far detector modules, a liquid argon time projection chamber (LArTPC) with horizontal drift (FD-HD) of ionization electrons. The reconstruction is based upon the workflow developed for DUNE's long-baseline oscillation analysis, with some necessary machine-learning models' retraining and the addition of features relevant only to atmospheric neutrinos such as the neutrino direction reconstruction. Where relevant, the impact of the detection of the charged particles of the hadronic system is emphasized, and comparisons are carried out between the case when lepton-only information is considered in the reconstruction (as is the case for many neutrino oscillation experiments), versus when all particles identified in the LArTPC were included. Three neutrino direction reconstruction methods have been developed and studied for the atmospheric analyses: using lepton-only information, using all reconstructed particles, and using only correlations from reconstructed hits. The results indicate that incorporating more than just lepton information significantly improves the resolution of both neutrino direction and energy reconstruction. The angle reconstruction algorithms developed in this work result in no strong dependence on particle direction for reconstruction efficiencies or neutrino flavor identification. This comprehensive review of the reconstruction of atmospheric neutrinos in DUNE's FD-HD LArTPC is the first step towards developing a first neutrino oscillation sensitivity analysis, which will ready DUNE for its first measurements.

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

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.

preprint2022arXiv

A COHERENT constraint on leptophobic dark matter using CsI data

We use data from the COHERENT CsI[Na] scintillation detector to constrain sub-GeV leptophobic dark matter models. This detector was built to observe low-energy nuclear recoils from coherent elastic neutrino-nucleus scattering. These capabilities enable searches for dark matter particles produced at the Spallation Neutron Source mediated by a vector portal particle with masses between 2 and 400 MeV/c$^2$. No evidence for dark matter is observed and a limit on the mediator coupling to quarks is placed. This constraint improves upon previous results by two orders of magnitude. This newly explored parameter space probes the region where the dark matter relic abundance is explained by leptophobic dark matter when the mediator mass is roughly twice the dark matter mass. COHERENT sets the best constraint on leptophobic dark matter at these masses.

preprint2022arXiv

A new study of the $N=32$ and $N=34$ shell gap for Ti and V by the first high-precision MRTOF mass measurements at BigRIPS-SLOWRI

The atomic masses of $^{55}$Sc, $^{56,58}$Ti, and $^{56-59}$V have been determined using the high-precision multi-reflection time-of-flight technique. The radioisotopes have been produced at RIKEN's RIBF facility and delivered to the novel designed gas cell and multi-reflection system (ZD MRTOF), which has been recently commissioned downstream of the ZeroDegree spectrometer following the BigRIPS separator. For $^{56,58}$Ti and $^{56-59}$V the mass uncertainties have been reduced down to the order of $10\,\mathrm{keV}$, shedding new light on the $N=34$ shell effect in Ti and V isotopes by the first high-precision mass measurements of the critical species $^{58}$Ti and $^{59}$V. With the new precision achieved, we reveal the non-existence of the $N=34$ empirical two-neutron shell gaps for Ti and V, and the enhanced energy gap above the occupied $νp_{3/2}$ orbit is identified as a feature unique to Ca. We perform new Monte Carlo shell model calculations including the $νd_{5/2}$ and $νg_{9/2}$ orbits and compare the results with conventional shell model calculations, which exclude the $νg_{9/2}$ and the $νd_{5/2}$ orbits. The comparison indicates that the shell gap reduction in Ti is related to a partial occupation of the higher orbitals for the outer two valence neutrons at $N=34$.

preprint2022arXiv

All-sky search for gravitational wave emission from scalar boson clouds around spinning black holes in LIGO O3 data

This paper describes the first all-sky search for long-duration, quasi-monochromatic gravitational-wave signals emitted by ultralight scalar boson clouds around spinning black holes using data from the third observing run of Advanced LIGO. We analyze the frequency range from 20~Hz to 610~Hz, over a small frequency derivative range around zero, and use multiple frequency resolutions to be robust towards possible signal frequency wanderings. Outliers from this search are followed up using two different methods, one more suitable for nearly monochromatic signals, and the other more robust towards frequency fluctuations. We do not find any evidence for such signals and set upper limits on the signal strain amplitude, the most stringent being $\approx10^{-25}$ at around 130~Hz. We interpret these upper limits as both an "exclusion region" in the boson mass/black hole mass plane and the maximum detectable distance for a given boson mass, based on an assumption of the age of the black hole/boson cloud system.

preprint2022arXiv

An auto-locking ultra-stable laser with sub-hertz linewidth

We report in detail the design process and performance of an auto-locking ultra-stable laser with sub-hertz linewidth at the first time. The laser frequency is automatically stabilized to an optical reference cavity with a home-made controller, which is based on a combination of digital circuit and analog circuit. The digital circuit is used for diagnosing and manipulating the state of the ultra-stable laser, and the analog circuit is used for demodulating the discriminate signal and servo control. A method of searching the transmission signal in the closed-loop state instead of the open-loop state is proposed to reduce the locking time and improve the reliable of the auto-locking ultra-stable laser. The median time of 16.6s is obtained after 157 times of relocking, and the probability of less than 20 s is more than 86%. The median linewidth of 1.08 Hz is obtained, and the fractional frequency instability is less than 3.4E-15 at integration time between 0.1 and 100 s. The performance of this system demonstrates that will be used as an important subsystem to transfer the optical clock signal.

preprint2022arXiv

Antiferromagnetic Excitonic Insulator State in Sr3Ir2O7

Excitonic insulators are usually considered to form via the condensation of a soft charge mode of bound electron-hole pairs. This, however, presumes that the soft exciton is of spin-singlet character. Early theoretical considerations have also predicted a very distinct scenario, in which the condensation of magnetic excitons results in an antiferromagnetic excitonic insulator state. Here we report resonant inelastic x-ray scattering (RIXS) measurements of Sr3Ir2O7. By isolating the longitudinal component of the spectra, we identify a magnetic mode that is well-defined at the magnetic and structural Brillouin zone centers, but which merges with the electronic continuum in between these high-symmetry points and which decays upon heating concurrent with a decrease in the material's resistivity. We show that a bilayer Hubbard model, in which electron-hole pairs are bound by exchange interactions, consistently explains all the electronic and magnetic properties of Sr3Ir2O7 indicating that this material is a realization of the long-predicted antiferromagnetic excitonic insulators phase.

preprint2022arXiv

Cavity-mediated electron-photon pairs

Advancing quantum information, communication and sensing relies on the generation and control of quantum correlations in complementary degrees of freedom. Here, we demonstrate the preparation of electron-photon pair states using the phase-matched interaction of free electrons with the evanescent vacuum field of a photonic-chip-based optical microresonator. Spontaneous inelastic scattering produces intracavity photons coincident with energy-shifted electrons. Harnessing these pairs for correlation-enhanced imaging, we achieve a two-orders of magnitude contrast improvement in cavity-mode mapping by coincidence-gated electron spectroscopy. This parametric pair-state preparation will underpin the future development of free-electron quantum optics, providing a pathway to quantum-enhanced imaging, electron-photon entanglement, and heralded single-electron and Fock-state photon sources.

preprint2022arXiv

Constraints from LIGO O3 data on gravitational-wave emission due to r-modes in the glitching pulsar PSR J0537-6910

We present a search for continuous gravitational-wave emission due to r-modes in the pulsar PSR J0537-6910 using data from the LIGO-Virgo Collaboration observing run O3. PSR J0537-6910 is a young energetic X-ray pulsar and is the most frequent glitcher known. The inter-glitch braking index of the pulsar suggests that gravitational-wave emission due to r-mode oscillations may play an important role in the spin evolution of this pulsar. Theoretical models confirm this possibility and predict emission at a level that can be probed by ground-based detectors. In order to explore this scenario, we search for r-mode emission in the epochs between glitches by using a contemporaneous timing ephemeris obtained from NICER data. We do not detect any signals in the theoretically expected band of 86-97 Hz, and report upper limits on the amplitude of the gravitational waves. Our results improve on previous amplitude upper limits from r-modes in J0537-6910 by a factor of up to 3 and place stringent constraints on theoretical models for r-mode driven spin-down in PSR J0537-6910, especially for higher frequencies at which our results reach below the spin-down limit defined by energy conservation.

preprint2022arXiv

Controllable emergent spatial spin modulation in Sr2IrO4 by in situ shear strain

Symmetric anisotropic interaction can be ferromagnetic and antiferromagnetic at the same time but for different crystallographic axes. We show that inducing competition of anisotropic interactions of orthogonal irreducible representations represents a general route to obtain new exotic magnetic states. We demonstrate it here by observing the emergence of a continuously tunable 12-layer spatial spin modulation when distorting the square lattice planes in the quasi-2D antiferromagnetic Sr2IrO4 under in situ shear strain. This translation-symmetry-breaking phase is a result of an unusual strain activated anisotropic interaction which is at the 4th order and competing with the inherent quadratic anisotropic interaction. Such a mechanism of competing anisotropy is distinct from that among the ferromagnetic, antiferromagnetic, and/or the Dzyaloshinskii-Moriya interactions, and it could be widely applicable and highly controllable in low dimensional magnets.

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

Dirac Monopoles in Rogue Waves

We for the first time demonstrate that the widely existed nonlinear waves such as rogue waves, contain Dirac monopoles. We find that the density zeros of these nonlinear waves on an extended complex plane can constitute the Dirac virtual magnetic monopole fields with a quantized flux of elementary $π$. We then can explain the exotic property of ``appearing from nowhere and disappearing without a trace" of rogue waves by means of a Dirac monopole collision mechanism. The maximum amplification ratio and multiple phase steps of high-order rogue waves are found to be closely related to the number of their contained monopoles. Important implications of the intrinsic virtual magnetic monopole fields are discussed.

preprint2022arXiv

Effective Field Theory Analysis of CDMSlite Run 2 Data

CDMSlite Run 2 was a search for weakly interacting massive particles (WIMPs) with a cryogenic 600 g Ge detector operated in a high-voltage mode to optimize sensitivity to WIMPs of relatively low mass from 2 - 20 GeV/$c^2$. In this article, we present an effective field theory (EFT) analysis of the CDMSlite Run 2 data using an extended energy range and a comprehensive treatment of the expected background. A binned likelihood Bayesian analysis was performed on the recoil energy data, taking into account the parameters of the EFT interactions and optimizing the data selection with respect to the dominant background components. Energy regions within 5$σ$ of known activation peaks were removed from the analysis. The Bayesian evidences resulting from the different operator hypotheses show that the CDMSlite Run 2 data are consistent with the background-only models and do not allow for a signal interpretation assuming any additional EFT interaction. Consequently, upper limits on the WIMP mass and coupling-coefficient amplitudes and phases are presented for each EFT operator. These limits improve previous CDMSlite Run 2 bounds for WIMP masses above 5 GeV/$c^2$.

preprint2022arXiv

First joint observation by the underground gravitational-wave detector, KAGRA, with GEO600

We report the results of the first joint observation of the KAGRA detector with GEO600. KAGRA is a cryogenic and underground gravitational-wave detector consisting of a laser interferometer with three-kilometer arms, and located in Kamioka, Gifu, Japan. GEO600 is a British--German laser interferometer with 600 m arms, and located near Hannover, Germany. GEO600 and KAGRA performed a joint observing run from April 7 to 20, 2020. We present the results of the joint analysis of the GEO--KAGRA data for transient gravitational-wave signals, including the coalescence of neutron-star binaries and generic unmodeled transients. We also perform dedicated searches for binary coalescence signals and generic transients associated with gamma-ray burst events observed during the joint run. No gravitational-wave events were identified. We evaluate the minimum detectable amplitude for various types of transient signals and the spacetime volume for which the network is sensitive to binary neutron-star coalescences. We also place lower limits on the distances to the gamma-ray bursts analysed based on the non-detection of an associated gravitational-wave signal for several signal models, including binary coalescences. These analyses demonstrate the feasibility and utility of KAGRA as a member of the global gravitational-wave detector network.

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

Improving sensitivity of the ARIANNA detector by rejecting thermal noise with deep learning

The ARIANNA experiment is an Askaryan detector designed to record radio signals induced by neutrino interactions in the Antarctic ice. Because of the low neutrino flux at high energies ($E > 10^{16} $), the physics output is limited by statistics. Hence, an increase in sensitivity significantly improves the interpretation of data and offers the ability to probe new parameter spaces. The amplitudes of the trigger threshold are limited by the rate of triggering on unavoidable thermal noise fluctuations. We present a real-time thermal noise rejection algorithm that enables the trigger thresholds to be lowered, which increases the sensitivity to neutrinos by up to a factor of two (depending on energy) compared to the current ARIANNA capabilities. A deep learning discriminator, based on a Convolutional Neural Network (CNN), is implemented to identify and remove thermal events in real time. We describe a CNN trained on MC data that runs on the current ARIANNA microcomputer and retains 95 percent of the neutrino signal at a thermal noise rejection factor of $10^5$, compared to a template matching procedure which reaches only $10^2$ for the same signal efficiency. Then the results are verified in a lab measurement by feeding in generated neutrino-like signal pulses and thermal noise directly into the ARIANNA data acquisition system. Lastly, the same CNN is used to classify cosmic-rays events to make sure they are not rejected. The network classified 102 out of 104 cosmic-ray events as signal.

preprint2022arXiv

Measurement of the Coherent Elastic Neutrino-Nucleus Scattering Cross Section on CsI by COHERENT

We measured the cross section of coherent elastic neutrino-nucleus scattering (\cevns{}) using a CsI[Na] scintillating crystal in a high flux of neutrinos produced at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory. New data collected before detector decommissioning has more than doubled the dataset since the first observation of \cevns{}, achieved with this detector. Systematic uncertainties have also been reduced with an updated quenching model, allowing for improved precision. With these analysis improvements, the COHERENT collaboration determined the cross section to be $(165^{+30}_{-25})\times10^{-40}$~cm$^2$, consistent with the standard model, giving the most precise measurement of \cevns{} yet. The timing structure of the neutrino beam has been exploited to compare the \cevns{} cross section from scattering of different neutrino flavors. This result places leading constraints on neutrino non-standard interactions while testing lepton flavor universality and measures the weak mixing angle as $\sin^2θ_{W}=0.220^{+0.028}_{-0.026}$ at $Q^2\approx(50\text{ MeV})^2$

preprint2022arXiv

Measuring the Polarization Reconstruction Resolution of the ARIANNA Neutrino Detector with Cosmic Rays

The ARIANNA detector is designed to detect neutrinos with energies above $10^{17}$eV. Due to the similarities in generated radio signals, cosmic rays are often used as test beams for neutrino detectors. Some ARIANNA detector stations are equipped with antennas capable of detecting air showers. Since the radio emission properties of air showers are well understood, and the polarization of the radio signal can be predicted from the arrival direction, cosmic rays can be used as a proxy to assess the reconstruction capabilities of the ARIANNA neutrino detector. We report on dedicated efforts of reconstructing the polarization of cosmic-ray radio pulses. After correcting for difference in hardware, the two stations used in this study showed similar performance in terms of event rate and agreed with simulation. Subselecting high quality cosmic rays, the polarizations of these cosmic rays were reconstructed with a resolution of $2.5^{\circ}$ (68% containment), which agrees with the expected value obtained from simulation. A large fraction of this resolution originates from uncertainties in the predicted polarization because of the contribution of the subdominant Askaryan effect in addition to the dominant geomagnetic emission. Subselecting events with a zenith angle greater than $70^{\circ}$ removes most influence of the Askaryan emission, and, with limited statistics, we found the polarization uncertainty is reduced to $1.3^{\circ}$ (68% containment).

preprint2022arXiv

Monitoring the SNS basement neutron background with the MARS detector

We present the analysis and results of the first dataset collected with the MARS neutron detector deployed at the Oak Ridge National Laboratory Spallation Neutron Source (SNS) for the purpose of monitoring and characterizing the beam-related neutron (BRN) background for the COHERENT collaboration. MARS was positioned next to the COH-CsI coherent elastic neutrino-nucleus scattering detector in the SNS basement corridor. This is the basement location of closest proximity to the SNS target and thus, of highest neutrino flux, but it is also well shielded from the BRN flux by infill concrete and gravel. These data show the detector registered roughly one BRN per day. Using MARS' measured detection efficiency, the incoming BRN flux is estimated to be $1.20~\pm~0.56~\text{neutrons}/\text{m}^2/\text{MWh}$ for neutron energies above $\sim3.5$ MeV and up to a few tens of MeV. We compare our results with previous BRN measurements in the SNS basement corridor reported by other neutron detectors.

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

Physics Opportunities in the ORNL Spallation Neutron Source Second Target Station Era

The Oak Ridge National Laboratory (ORNL) Spallation Neutron Source (SNS) First Target Station (FTS), used by the COHERENT experiment, provides an intense and extremely high-quality source of pulsed stopped-pion neutrinos, with energies up to about 50 MeV. Upgrades to the SNS are planned, including a Second Target Station (STS), which will approximately double the expected neutrino flux while maintaining quality similar to the FTS source. Furthermore, additional space for ten-tonne scale detectors may be available. We describe here exciting opportunities for neutrino physics, other particle and nuclear physics, and detector development using the FTS and STS neutrino sources.

preprint2022arXiv

Scintillator Tile Batch Test of CEPC AHCAL

Hadron calorimeter (HCAL) is an essential sub-detector of the baseline detector system for Circular Electron Positron Collider (CEPC). We plan to build an Analog Hadron CALorimeter (AHCAL) prototype based on the Particle Flow Algorithm (PFA). The AHCAL of CEPC uses steel as absorber and scintillator tiles read out by Silicon Photo-Multipliers (SiPMs) as sensitive medium. The energy linearity and resolution of the calorimeter depends on the light yield uniformity of sensitive medium. It is essential to qualify the entire detector production in order to select scintillator tiles with light yield uniform within 10\%. An automated batch test platform has been designed with 144 channels, an automated 3D servo motor. The paper summarizes the tests performed on more than 15000 scintillator tiles. The measured light yield, corrected for the set-up response non-uniformity, is around 12.9 p.e. . About 91.6\% of scintillators (14219 pieces) are qualified within 10\% of light yield window.

preprint2022arXiv

Search for anisotropic gravitational-wave backgrounds using data from Advanced LIGO and Advanced Virgo's first three observing runs

We report results from searches for anisotropic stochastic gravitational-wave backgrounds using data from the first three observing runs of the Advanced LIGO and Advanced Virgo detectors. For the first time, we include Virgo data in our analysis and run our search with a new efficient pipeline called {\tt PyStoch} on data folded over one sidereal day. We use gravitational-wave radiometry (broadband and narrow band) to produce sky maps of stochastic gravitational-wave backgrounds and to search for gravitational waves from point sources. A spherical harmonic decomposition method is employed to look for gravitational-wave emission from spatially-extended sources. Neither technique found evidence of gravitational-wave signals. Hence we derive 95\% confidence-level upper limit sky maps on the gravitational-wave energy flux from broadband point sources, ranging from $F_{α, Θ} < {\rm (0.013 - 7.6)} \times 10^{-8} {\rm erg \, cm^{-2} \, s^{-1} \, Hz^{-1}},$ and on the (normalized) gravitational-wave energy density spectrum from extended sources, ranging from $Ω_{α, Θ} < {\rm (0.57 - 9.3)} \times 10^{-9} \, {\rm sr^{-1}}$, depending on direction ($Θ$) and spectral index ($α$). These limits improve upon previous limits by factors of $2.9 - 3.5$. We also set 95\% confidence level upper limits on the frequency-dependent strain amplitudes of quasimonochromatic gravitational waves coming from three interesting targets, Scorpius X-1, SN 1987A and the Galactic Center, with best upper limits range from $h_0 < {\rm (1.7-2.1)} \times 10^{-25},$ a factor of $\geq 2.0$ improvement compared to previous stochastic radiometer searches.

preprint2022arXiv

Search for continuous gravitational wave emission from the Milky Way center in O3 LIGO--Virgo data

We present a directed search for continuous gravitational wave (CW) signals emitted by spinning neutron stars located in the inner parsecs of the Galactic Center (GC). Compelling evidence for the presence of a numerous population of neutron stars has been reported in the literature, turning this region into a very interesting place to look for CWs. In this search, data from the full O3 LIGO--Virgo run in the detector frequency band $[10,2000]\rm~Hz$ have been used. No significant detection was found and 95$\%$ confidence level upper limits on the signal strain amplitude were computed, over the full search band, with the deepest limit of about $7.6\times 10^{-26}$ at $\simeq 142\rm~Hz$. These results are significantly more constraining than those reported in previous searches. We use these limits to put constraints on the fiducial neutron star ellipticity and r-mode amplitude. These limits can be also translated into constraints in the black hole mass -- boson mass plane for a hypothetical population of boson clouds around spinning black holes located in the GC.

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

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in the Second and Third LIGO-Virgo Observing Runs

We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both once and twice the rotation frequency (dual harmonic). No evidence of GWs was found so we present 95\% credible upper limits on the strain amplitudes $h_0$ for the single harmonic search along with limits on the pulsars' mass quadrupole moments $Q_{22}$ and ellipticities $\varepsilon$. Of the pulsars studied, 23 have strain amplitudes that are lower than the limits calculated from their electromagnetically measured spin-down rates. These pulsars include the millisecond pulsars J0437\textminus4715 and J0711\textminus6830 which have spin-down ratios of 0.87 and 0.57 respectively. For nine pulsars, their spin-down limits have been surpassed for the first time. For the Crab and Vela pulsars our limits are factors of $\sim 100$ and $\sim 20$ more constraining than their spin-down limits, respectively. For the dual harmonic searches, new limits are placed on the strain amplitudes $C_{21}$ and $C_{22}$. For 23 pulsars we also present limits on the emission amplitude assuming dipole radiation as predicted by Brans-Dicke theory.

preprint2022arXiv

Simulating the neutrino flux from the Spallation Neutron Source for the COHERENT experiment

The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory is a pulsed source of neutrons and, as a byproduct of this operation, an intense source of pulsed neutrinos via stopped-pion decay. The COHERENT collaboration uses this source to investigate coherent elastic neutrino-nucleus scattering and other physics with a suite of detectors. This work includes a description of our Geant4 simulation of neutrino production at the SNS and the flux calculation which informs the COHERENT studies. We estimate the uncertainty of this calculation at about 10% based on validation against available low-energy pion production data.

preprint2022arXiv

Snowmass 2021 White Paper Instrumentation Frontier 05 -- White Paper 1: MPGDs: Recent advances and current R&D

This paper will review the origins, development, and examples of new versions of Micro-Pattern Gas Detectors. The goal for MPGD development was the creation of detectors that could cost-effectively cover large areas while offering excellent position and timing resolution, and the ability to operate at high incident particle rates. The early MPGD developments culminated in the formation of the RD51 collaboration which has become the critical organization for the promotion of MPGDs and all aspects of their production, characterization, simulation, and uses in an expanding array of experimental configurations. For the Snowmass 2021 study, a number of Letters of Interest were received that illustrate ongoing developments and expansion of the use of MPGDs. In this paper, we highlight high precision timing, high rate application, trigger capability expansion of the SRS readout system, and a structure designed for low ion backflow.

preprint2022arXiv

The COHERENT Experimental Program

The COHERENT experiment located in Neutrino Alley at the Spallation Neutron Source (SNS), Oak Ridge National Laboratory (ORNL), has made the world's first two measurements of coherent elastic neutrino-nucleus scattering (CEvNS), on CsI and argon, using neutrinos produced at the SNS. The COHERENT collaboration continues to pursue CEvNS measurements on various targets as well as additional studies of inelastic neutrino-nucleus interactions, searches for accelerator-produced dark matter (DM) and physics beyond the Standard Model, using the uniquely high-quality and high-intensity neutrino source available at the SNS. This white paper describes primarily COHERENT's ongoing and near-future program at the SNS First Target Station (FTS). Opportunities enabled by the SNS Second Target Station (STS) for the study of neutrino physics and development of novel detector technologies are elaborated in a separate white paper.

preprint2022arXiv

The Proton Spin Structure Function $g_2$ and Generalized Polarizabilities in the Strong QCD Regime

The strong interaction is not well understood at low energy, or for interactions with low momentum transfer $Q^2$, but one of the clearest insights we have comes from Chiral Perturbation Theory ($χ$PT). This effective treatment gives testable predictions for the nucleonic generalized polarizabilities -- fundamental quantities describing the nucleon's response to an external field. We have measured the proton's generalized spin polarizabilities in the region where $χ$PT is expected to be valid. Our results include the first ever data for the transverse-longitudinal spin polarizability $δ_{LT}$, and also extend the coverage of the polarizability $\bar{d_2}$ to very low $Q^2$ for the first time. These results were extracted from moments of the structure function $g_2$, a quantity which characterizes the internal spin structure of the proton. Our experiment ran at Jefferson Lab using a polarized electron beam and a polarized solid ammonia (NH$_3$) target. The $δ_{LT}$ polarizability has remained a challenging quantity for $χ$PT to reproduce, despite its reduced sensitivity to higher resonance contributions; recent competing calculations still disagree with each other and also diverge from the measured neutron data at very low $Q^2$. Our proton results provide discriminating power between existing calculations, and will help provide a better understanding of this strong QCD regime.

preprint2022arXiv

Triboelectric Backgrounds to radio-based UHE Neutrino Exeperiments

The proposed IceCube-Gen2 (ICG2) seeks to instrument ~500 sq. km of Antarctic ice near the geographic South Pole with radio antennas, in order to observe the highest energy (E>1 EeV) neutrinos in the Universe. To this end, ICG2 will use the impulsive radio-frequency (RF) signal produced by neutrino interactions in polar ice caps. In such experiments, rare single event candidates must be unambiguously separated from background; to date, signal identification strategies primarily reject thermal noise and anthropogenic backgrounds. Here, we consider the possibility that fake neutrino signals may also be naturally generated via the 'triboelectric effect'. This broadly includes any process in which force applied at a boundary layer results in displacement of surface charge, generating a potential difference ΔV. Wind blowing over granular surfaces such as snow can induce such a ΔV, with subsequent discharge. Discharges over nanosecond-timescales can then lead to RF emissions at characteristic MHz-GHz frequencies. We find that such backgrounds are evident in the several neutrino experiments considered, and are generally characterized by: a) a threshold wind velocity which likely depends on the experimental signal trigger threshold and layout; for the experiments considered herein, this value is typically O(10 m/s), b) frequency spectra generally shifted to the low-end of the frequency regime to which current radio experiments are typically sensitive (100-200 MHz), c) for the strongest background signals, an apparent preference for discharges from above-surface structures, although the presence of more isotropic, lower amplitude triboelectric discharges cannot be excluded.

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

Characterization of a CMOS pixel sensor for charged particle tracking

A prototype of the CMOS pixel sensor named SUPIX-1 has been fabricated and tested in order to investigate the feasibility of a pixelated tracker for a proposed Higgs factory, namely, the Circular Electron-Positron Collider (CEPC). The sensor, taped out with a 180 nm CMOS Image Sensor (CIS) process, consists of nine different pixel arrays varying in pixel pitches, diode sizes and geometries in order to study the particle detection performance of enlarged pixels. The test was carried out with a Fe-55 radioactive source. Two soft X-ray peaks observed were used to calibrate the charge to voltage factor of the sensor. The pixel-wise equivalent noise charge, charge collection efficiency and signal-to-noise ratio were evaluated. A reconstruction method for clustering pixels of a signal has been developed and the cluster-wise performance was studied as well. The test results show that pixels with the area as large as of $21μ$m$\times 84μ$m have satisfactory noise level and charge collection performance, meeting gerenal requirements for a pixel sensor. This contribution demonstrates that the CMOS pixel sensor with enlarged pitches, using the CIS technology, can be used in tracking for upcoming collider detectors akin to the CEPC.

preprint2021arXiv

Development of Planar P-type Point Contact Germanium Detectors for Low-Mass Dark Matter Searches

The detection of low-energy deposition in the range of sub-eV through ionization using germanium (Ge) with a bandgap of $\sim$0.7 eV requires internal amplification of charge signal. This can be achieved through high electric field which accelerates charge carriers to generate more charge carriers. The minimum electric field required to generate internal charge amplification is derived for different temperatures. A point contact Ge detector provides extremely high electric field in proximity to the point contact. We show the development of a planar point contact detector and its performance. The field distribution is calculated for this planar point contact detector. We demonstrate the required electric field can be achieved with a point contact detector.

preprint2021arXiv

First Measurement of Coherent Elastic Neutrino-Nucleus Scattering on Argon

We report the first measurement of coherent elastic neutrino-nucleus scattering (\cevns) on argon using a liquid argon detector at the Oak Ridge National Laboratory Spallation Neutron Source. Two independent analyses prefer \cevns over the background-only null hypothesis with greater than $3σ$ significance. The measured cross section, averaged over the incident neutrino flux, is (2.2 $\pm$ 0.7) $\times$10$^{-39}$ cm$^2$ -- consistent with the standard model prediction. The neutron-number dependence of this result, together with that from our previous measurement on CsI, confirms the existence of the \cevns process and provides improved constraints on non-standard neutrino interactions.

preprint2021arXiv

Modeling Homophily in Dynamic Networks with Application to HIV Molecular Surveillance

This paper describes a novel approach to modeling homphily, i.e. the tendency of nodes that share (or differ in) certain attributes to be linked; we consider dynamic networks in which nodes can be added over time but not removed. Our application is to HIV genetic linkage analysis that has been used to investigate HIV transmission dynamics. In this setting, two HIV sequences from different persons with HIV (PWH) are said to be linked if the genetic distance between these sequences is less than a given threshold. Such linkage suggests that that the nodes representing the two infected PWH, are close to each other in a transmission network; such proximity would imply that either one of the infected people directly transmitted the virus to the other or indirectly transmitted it through a small number of intermediaries. These viral genetic linkage networks are dynamic in the sense that, over time, a group or cluster of genetically linked viral sequences may increase in size as new people are infected by those in the cluster either directly or through intermediaries. Our approach makes use of a logistic model to describe homophily with regard to demographic and behavioral characteristics that is we investigate whether similarities (or differences) between PWH in these characteristics impacts the probability that their sequences are be linked. Such analyses provide information about HIV transmission dynamics within a population.

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

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.

preprint2021arXiv

Semiclassical Trajectory Perspective of Glory Rescattering in Strong-field Photoelectron Holography

We investigate theoretically the photoelectron momentum distribution (PMD) of the ionized atoms irradiated by a linearly polarized intense laser, focusing on the holography interference patterns in PMD that carry important information of the initial wavefunction of a tunneled electron and its experienced atomic potential in rescattering. With the help of Dyson series and semiclassical propagator, we calculate the scattering amplitudes in the cylindrical coordinate representation. In contrast to conventional recognitions that the photoelectron holography is the interference of two branches of electron trajectories, however, we find strikingly that infinite semiclassical trajectories can be deflected by the combined Coulomb potential and laser field into the same final momentum: The initial momenta are found to be distributed on a ring-shape curve in the transverse momentum plane and the initial positions of these trajectories are perpendicular to their initial momentum vectors. For the zero final transverse momentum, the above ring-source trajectories degenerate into the point-source axial caustic trajectories (or Glory trajectories) and the quantum interference of these trajectories will dramatically alter the scattering amplitudes that is termed as Glory rescattering effect. With following Berry's spirit of uniform approximation for Glory scattering in optics, we can finally derive a uniform formulation of the rescattering amplitude in the Bessel functions for the strong-field photoelectron holography (SFPH) patterns. Our results are in good agreement with solutions of the time-dependent Schrödinger equation and can account for recent photoelectron holography experiments. Important applications of our theory are also discussed.

preprint2021arXiv

Timing performance of a multi-pad PICOSEC-Micromegas detector prototype

The multi-pad PICOSEC-Micromegas is an improved detector prototype with a segmented anode, consisting of 19 hexagonal pads. Detailed studies are performed with data collected in a muon beam over four representative pads. We demonstrate that such a device, scalable to a larger area, provides excellent time resolution and detection efficiency. As expected from earlier single-cell device studies, we measure a time resolution of approximately 25 picoseconds for charged particles hitting near the anode pad centers, and up to 30 picoseconds at the pad edges. Here, we study in detail the effect of drift gap thickness non-uniformity on the timing performance and evaluate impact position based corrections to obtain a uniform timing response over the full detector coverage.

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 TOSA/ROSA-Based Optical Transmitter (MTx+)/Transceiver (MTRx+) for High-Energy Physics Experiments

We present a dual-channel optical transmitter (MTx+)/transceiver (MTRx+) for the front-end readout electronics of high-energy physics experiments. MTx+ utilizes two Transmitter Optical Sub-Assemblies (TOSAs) and MTRx+ utilizes a TOSA and a Receiver Optical Sub-Assemblies (ROSA). Both MTx+ and MTRx+ receive multimode fibers with standard Lucent Connectors (LCs) as the optical interface and can be panel or board mounted to a motherboard with a standard Enhanced Small Form-factor Pluggable (SFP+) connector as the electrical interface. MTx+ and MTRx+ employ a dual-channel Vertical-Cavity Surface-Emitting Laser (VCSEL) driver ASIC called LOCld65, which brings the transmitting data rate up to 14 Gbps per channel. MTx+ and MTRx+ have been tested to survive 4.9 kGy(SiO2).

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

Binary Black Hole Population Properties Inferred from the First and Second Observing Runs of Advanced LIGO and Advanced Virgo

We present results on the mass, spin, and redshift distributions with phenomenological population models using the ten binary black hole mergers detected in the first and second observing runs completed by Advanced LIGO and Advanced Virgo. We constrain properties of the binary black hole (BBH) mass spectrum using models with a range of parameterizations of the BBH mass and spin distributions. We find that the mass distribution of the more massive black hole in such binaries is well approximated by models with no more than 1% of black holes more massive than $45\,M_\odot$, and a power law index of $α= {1.3}^{+1.4}_{-1.7}$ (90% credibility). We also show that BBHs are unlikely to be composed of black holes with large spins aligned to the orbital angular momentum. Modelling the evolution of the BBH merger rate with redshift, we show that it is at or increasing with redshift with 93% probability. Marginalizing over uncertainties in the BBH population, we find robust estimates of the BBH merger rate density of $R = {53.2}^{+55.8}_{-28.2}$ Gpc$^{-3}$ yr$^{-1}$ (90% credibility). As the BBH catalog grows in future observing runs, we expect that uncertainties in the population model parameters will shrink, potentially providing insights into the formation of black holes via supernovae, binary interactions of massive stars, stellar cluster dynamics, and the formation history of black holes across cosmic time.

preprint2020arXiv

Characterization of High-Purity Germanium Detectors with Amorphous Germanium Contacts in Cryogenic Liquids

For the first time, planar high-purity germanium detectors with thin amorphous germanium contacts were successfully operated directly in liquid nitrogen and liquid argon in a cryostat at the Max-Planck-Institut fuer Physics in Munich. The detectors were fabricated at the Lawrence Berkeley National Laboratory and the University of South Dakota, using crystals grown at the University of South Dakota. They survived long-distance transportation and multiple thermal cycles in both cryogenic liquids and showed reasonable leakage currents and spectroscopic performance. Also discussed are the pros and cons of using thin amorphous semiconductor materials as an alternative contact technology in large-scale germanium experiments searching for physics beyond the Standard Model.

preprint2020arXiv

COHERENT Collaboration data release from the first detection of coherent elastic neutrino-nucleus scattering on argon

Release of COHERENT collaboration data from the first detection of coherent elastic neutrino-nucleus scattering (CEvNS) on argon. This release corresponds with the results of "Analysis A" published in Akimov et al., arXiv:2003.10630 [nucl-ex]. Data is shared in a binned, text-based format representing both "signal" and "backgrounds" along with associated uncertainties such that the included data can be used to perform independent analyses. This document describes the contents of the data release as well as guidance on the use of the data. Included example code in C++ (ROOT) and Python show one possible use of the included data.

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

Helicon plasma in a magnetic shuttle

The definition of magnetic shuttle is introduced to describe the magnetic space enclosed by two tandem magnetic mirrors with the same field direction and high mirror ratio. Helicon plasma immersed in such a magnetic shuttle which can provide the confinement of charged particles is modeled using an electromagnetic solver. The perpendicular structure of wave field along this shuttle is given in terms of stream vector plots, showing significant change from midplane to ending throats, and the vector field rotates and forms a circular layer that separating plasma column radially into core and edge regions near the throats. The influences of driving frequency, plasma density and field strength on the wave field and power absorption are computed in detail. It is found that the wave magnitude and power absorption decrease for increased driving frequency and reduced field strength, and maximize around a certain level of plasma density. The axial standing-wave feature always exists, due to the interference between forward and reflected waves from ending magnetic mirrors, while the radial wave field structure largely stays the same. Distributions of wave energy density and power absorption density all show shrinking feature from midplane to ending throats, which is consistent with the nature of helicon mode that propagating along field lines. Theoretical analysis based on a simple magnetic shuttle and the governing equation of helicon waves shows consistency with computed results and previous studies.

preprint2020arXiv

Impact of Charge Trapping on the Energy Resolution of Ge Detectors for Rare-Event Physics Searches

Charge trapping degrades the energy resolution of germanium (Ge) detectors, which require to have increased experimental sensitivity in searching for dark matter and neutrinoless double-beta decay. We investigate the charge trapping processes utilizing nine planar detectors fabricated from USD-grown crystals with well-known net impurity levels. The charge collection efficiency as a function of charge trapping length is derived from the Shockley-Ramo theorem. Furthermore, we develop a model that correlates the energy resolution with the charge collection efficiency. This model is then applied to the experimental data. As a result, charge collection efficiency and charge trapping length are determined accordingly. Utilizing the Lax model (further developed by CDMS collaborators), the absolute impurity levels are determined for nine detectors. The knowledge of these parameters when combined with other traits such as the Fano factor serve as a reliable indicator of the intrinsic nature of charge trapping within the crystals. We demonstrate that electron trapping is more severe than hole trapping in a p-type detector and the charge collection efficiency depends on the absolute impurity level of the Ge crystal when an adequate bias voltage is applied to the detector. Negligible charge trapping is found when the absolute impurity level is less than 1.0$\times$10$^{11}/$cm$^{3}$ for collecting electrons and 2.0$\times$10$^{11}/$cm$^{3}$ for collecting holes.

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

Observation of a cross-section enhancement near mass threshold in $e^{+}e^{-}\rightarrowΛ\barΛ$

The process $e^{+}e^{-}\rightarrowΛ\barΛ$ is studied using data samples at $\sqrt{s}=2.2324$, 2.400, 2.800 and 3.080 GeV collected with the BESIII detector operating at the BEPCII collider. The Born cross section is measured at $\sqrt{s}$=2.2324 GeV, which is 1.0 MeV above the $Λ\barΛ$ mass threshold, to be $305\pm45^{+66}_{-36}$ pb, where the first uncertainty is statistical and the second systematic. The substantial cross section near threshold is significantly larger than that expected from theory, which predicts the cross section to vanish at threshold. The Born cross sections at $\sqrt{s}$=2.400, 2.800 and 3.080 GeV are measured and found to be consistent with previous experimental results, but with improved precision. Finally, the corresponding effective electromagnetic form factors of $Λ$ are deduced.

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

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in 2015-2017 LIGO Data

We present a search for gravitational waves from 221 pulsars with rotation frequencies $\gtrsim 10$ Hz. We use advanced LIGO data from its first and second observing runs spanning 2015-2017, which provides the highest-sensitivity gravitational-wave data so far obtained. In this search we target emission from both the $l = m = 2$ mass quadrupole mode, with a frequency at twice that of the pulsar's rotation, and from the $l = 2$, $m = 1$ mode, with a frequency at the pulsar rotation frequency. The search finds no evidence for gravitational-wave emission from any pulsar at either frequency. For the $l = m = 2$ mode search, we provide updated upper limits on the gravitational-wave amplitude, mass quadrupole moment, and fiducial ellipticity for 167 pulsars, and the first such limits for a further 55. For 20 young pulsars these results give limits that are below those inferred from the pulsars' spin-down. For the Crab and Vela pulsars our results constrain gravitational-wave emission to account for less than 0.017% and 0.18% of the spin-down luminosity, respectively. For the recycled millisecond pulsar J0711-6830 our limits are only a factor of 1.3 above the spin-down limit, assuming the canonical value of $10^{38}$ kg m$^2$ for the star's moment of inertia, and imply a gravitational-wave-derived upper limit on the star's ellipticity of $1.2\!\times\!10^{-8}$. We also place new limits on the emission amplitude at the rotation frequency of the pulsars.

preprint2020arXiv

Sensitivity and Performance of the Advanced LIGO Detectors in the Third Observing Run

On April 1st, 2019, the Advanced Laser Interferometer Gravitational-Wave Observatory (aLIGO), joined by the Advanced Virgo detector, began the third observing run, a year-long dedicated search for gravitational radiation. The LIGO detectors have achieved a higher duty cycle and greater sensitivity to gravitational waves than ever before, with LIGO Hanford achieving angle-averaged sensitivity to binary neutron star coalescences to a distance of 111 Mpc, and LIGO Livingston to 134 Mpc with duty factors of 74.6% and 77.0% respectively. The improvement in sensitivity and stability is a result of several upgrades to the detectors, including doubled intracavity power, the addition of an in-vacuum optical parametric oscillator for squeezed-light injection, replacement of core optics and end reaction masses, and installation of acoustic mode dampers. This paper explores the purposes behind these upgrades, and explains to the best of our knowledge the noise currently limiting the sensitivity of each detector.

preprint2020arXiv

Solar Flare Prediction Using Magnetic Field Diagnostics Above the Photosphere

In this article, we present the application of the weighted horizontal gradient of magnetic field (WGM) flare prediction method to 3-dimensional (3D) extrapolated magnetic configurations of 13 flaring solar active regions (ARs). The main aim is to identify an optimal height range, if any, in the interface region between the photosphere and lower corona, where the flare onset time prediction capability of WGM is best exploited. The optimal height is where flare prediction, by means of the WGM method, is achieved earlier than at the photospheric level. 3D magnetic structures, based on potential and non-linear force-free field extrapolations, are constructed to study a vertical range from the photosphere up to the low corona with a 45 km step size. The WGM method is applied as a function of height to all 13 flaring AR cases that are subject to certain selection criteria. We found that applying the WGM method between 1000 and 1800 km above the solar surface would improve the prediction of the flare onset time by around 2-8 hours.Certain caveats and an outlook for future work along these lines are also discussed.

preprint2020arXiv

Spontaneous Rotation of Ferrimagnetism Driven by Antiferromagnetic Spin Canting

Spin-reorientation phase transitions that involve the rotation of a crystal$'$s magnetization have been well characterized in distorted-perovskite oxides such as the orthoferrites. In these systems spin reorientation occurs due to competing rare-earth and transition metal anisotropies coupled via $f$-$d$ exchange. Here, we demonstrate an alternative paradigm for spin reorientation in distorted perovskites. We show that the $R_2\mathrm{CuMnMn_4O_{12}}$ (R = Y or Dy) triple A-site columnar-ordered quadruple perovskites have three ordered magnetic phases and up to two spin-reorientation phase transitions. Unlike the spin-reorientation phenomena in other distorted perovskites, these transitions are independent of rare-earth magnetism, but are instead driven by an instability towards antiferromagnetic spin canting likely originating in frustrated Heisenberg exchange interactions, and the competition between Dzyaloshinskii-Moriya and single-ion anisotropies.

preprint2020arXiv

The ELFIN Mission

The Electron Loss and Fields Investigation with a Spatio-Temporal Ambiguity-Resolving option (ELFIN-STAR, or simply: ELFIN) mission comprises two identical 3-Unit (3U) CubeSats on a polar (~93deg inclination), nearly circular, low-Earth (~450 km altitude) orbit. Launched on September 15, 2018, ELFIN is expected to have a >2.5 year lifetime. Its primary science objective is to resolve the mechanism of storm-time relativistic electron precipitation, for which electromagnetic ion cyclotron (EMIC) waves are a prime candidate. From its ionospheric vantage point, ELFIN uses its unique pitch-angle-resolving capability to determine whether measured relativistic electron pitch-angle and energy spectra within the loss cone bear the characteristic signatures of scattering by EMIC waves or whether such scattering may be due to other processes. Pairing identical ELFIN satellites with slowly-variable along-track separation allows disambiguation of spatial and temporal evolution of the precipitation over minutes-to-tens-of-minutes timescales, faster than the orbit period of a single low-altitude satellite (~90min). Each satellite carries an energetic particle detector for electrons (EPDE) that measures 50keV to 5MeV electrons with deltaE/E<40% and a fluxgate magnetometer (FGM) on a ~72cm boom that measures magnetic field waves (e.g., EMIC waves) in the range from DC to 5Hz Nyquist (nominally) with <0.3nT/sqrt(Hz) noise at 1Hz. The spinning satellites (T_spin~3s) are equipped with magnetorquers that permit spin-up/down and reorientation maneuvers. The spin axis is placed normal to the orbit plane, allowing full pitch-angle resolution twice per spin. An energetic particle detector for ions (EPDI) measures 250keV-5MeV ions, addressing secondary science. Funded initially by CalSpace and the University Nanosat Program, ELFIN was selected for flight with joint support from NSF and NASA between 2014 and 2018.

preprint2020arXiv

The Impact of the Charge Barrier Height on Germanium (Ge) Detectors with Amorphous-Ge Contacts for Light Dark Matter Searches

Germanium (Ge) detectors with ability of measuring a single electron-hole (e-h) pair are needed in searching for light dark matter (LDM) down to the MeV scale. We investigate the feasibility of Ge detectors with amorphous-Ge (a-Ge) contacts to achieve the sensitivity of measuring a single e-h pair, which requires extremely low leakage current. Three Ge detectors with a-Ge contacts are used to study the charge barrier height for blocking electrons and holes. Using the measured bulk leakage current and the D$\ddot{o}$hler-Brodsky model, we obtain the values for charge barrier height and the density of localized energy states near the Fermi energy level for the top and bottom contacts, respectively. We predict that the bulk leakage current is extremely small and can be neglected at helium temperature ($\sim$4 K). Thus, Ge detectors with a-Ge contacts possess the potential to measure a single e-h pair for detecting LDM particles.

preprint2020arXiv

The observation of the Crab Nebula with LHAASO-KM2A for the performance study

As a sub-array of the Large High Altitude Air Shower Observatory (LHAASO), KM2A is mainly designed to cover a large fraction of the northern sky to hunt for gamma-ray sources at energies above 10 TeV. Even though the detector construction is still underway, a half of the KM2A array has been operating stably since the end of 2019. In this paper, we present the pipeline of KM2A data analysis and the first observation on the Crab Nebula, a standard candle in very high energy gamma-ray astronomy. We detect gamma-ray signals from the Crab Nebula in both energy ranges of 10$-$100 TeV and $>$100 TeV with high significance, by analyzing the KM2A data of 136 live days between December 2019 and May 2020. With the observations, we test the detector performance including angular resolution, pointing accuracy and cosmic ray background rejection power. The energy spectrum of the Crab Nebula in the energy range 10-250 TeV fits well with a single power-law function dN/dE =(1.13$\pm$0.05$_{stat}$$\pm$0.08$_{sys}$)$\times$10$^{-14}$$\cdot$(E/20TeV)$^{-3.09\pm0.06_{stat}\pm0.02_{sys}}$ cm$^{-2}$ s$^{-1}$ TeV$^{-1}$. It is consistent with previous measurements by other experiments. This opens a new window of gamma-ray astronomy above 0.1 PeV through which ultrahigh-energy gamma-ray new phenomena, such as cosmic PeVatrons, might be discovered.

preprint2017arXiv

A Hybrid Observer for a Distributed Linear System with a Changing Neighbor Graph

A hybrid observer is described for estimating the state of an $m>0$ channel, $n$-dimensional, continuous-time, distributed linear system of the form $\dot{x} = Ax,\;y_i = C_ix,\;i\in\{1,2,\ldots, m\}$. The system's state $x$ is simultaneously estimated by $m$ agents assuming each agent $i$ senses $y_i$ and receives appropriately defined data from each of its current neighbors. Neighbor relations are characterized by a time-varying directed graph $\mathbb{N}(t)$ whose vertices correspond to agents and whose arcs depict neighbor relations. Agent $i$ updates its estimate $x_i$ of $x$ at "event times" $t_1,t_2,\ldots $ using a local observer and a local parameter estimator. The local observer is a continuous time linear system whose input is $y_i$ and whose output $w_i$ is an asymptotically correct estimate of $L_ix$ where $L_i$ a matrix with kernel equaling the unobservable space of $(C_i,A)$. The local parameter estimator is a recursive algorithm designed to estimate, prior to each event time $t_j$, a constant parameter $p_j$ which satisfies the linear equations $w_k(t_j-τ) = L_kp_j+μ_k(t_j-τ),\;k\in\{1,2,\ldots,m\}$, where $τ$ is a small positive constant and $μ_k$ is the state estimation error of local observer $k$. Agent $i$ accomplishes this by iterating its parameter estimator state $z_i$, $q$ times within the interval $[t_j-τ, t_j)$, and by making use of the state of each of its neighbors' parameter estimators at each iteration. The updated value of $x_i$ at event time $t_j$ is then $x_i(t_j) = e^{Aτ}z_i(q)$. Subject to the assumptions that (i) the neighbor graph $\mathbb{N}(t)$ is strongly connected for all time, (ii) the system whose state is to be estimated is jointly observable, (iii) $q$ is sufficiently large, it is shown that each estimate $x_i$ converges to $x$ exponentially fast as $t\rightarrow \infty$ at a rate which can be controlled.

preprint2017arXiv

Search for High-energy Neutrinos from Binary Neutron Star Merger GW170817 with ANTARES, IceCube, and the Pierre Auger Observatory

The Advanced LIGO and Advanced Virgo observatories recently discovered gravitational waves from a binary neutron star inspiral. A short gamma-ray burst (GRB) that followed the merger of this binary was also recorded by the Fermi Gamma-ray Burst Monitor (Fermi-GBM), and the Anticoincidence Shield for the Spectrometer for the International Gamma-Ray Astrophysics Laboratory (INTEGRAL), indicating particle acceleration by the source. The precise location of the event was determined by optical detections of emission following the merger. We searched for high-energy neutrinos from the merger in the GeV--EeV energy range using the ANTARES, IceCube, and Pierre Auger Observatories. No neutrinos directionally coincident with the source were detected within $\pm500$ s around the merger time. Additionally, no MeV neutrino burst signal was detected coincident with the merger. We further carried out an extended search in the direction of the source for high-energy neutrinos within the 14-day period following the merger, but found no evidence of emission. We used these results to probe dissipation mechanisms in relativistic outflows driven by the binary neutron star merger. The non-detection is consistent with model predictions of short GRBs observed at a large off-axis angle.

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}$.

J. Liu

Quick read

Decide how to stay connected

How to connect with this researcher

Open a focused conversation when the fit is right

See the researcher in context

Building this graph slice

90 published item(s)

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

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

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

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

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

Discovery of a new $γ$-ray source LHAASO J0341+5258 with emission up to 200TeV

Discovery of the Ultra-high energy gamma-ray source LHAASO J2108+5157

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

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

Evidence for particle acceleration approaching PeV energies in the W51 complex

Exploring Lorentz Invariance Violation from Ultra-high-energy Gamma Rays Observed by LHAASO

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

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

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

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

Multimode Objective Lens for Momentum Microscopy and XPEEM: Experiments

Reconstruction of atmospheric neutrinos in DUNE&#39;s horizontal-drift far-detector module

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

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

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

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

A COHERENT constraint on leptophobic dark matter using CsI data

A new study of the $N=32$ and $N=34$ shell gap for Ti and V by the first high-precision MRTOF mass measurements at BigRIPS-SLOWRI

All-sky search for gravitational wave emission from scalar boson clouds around spinning black holes in LIGO O3 data

An auto-locking ultra-stable laser with sub-hertz linewidth

Antiferromagnetic Excitonic Insulator State in Sr3Ir2O7

Cavity-mediated electron-photon pairs

Constraints from LIGO O3 data on gravitational-wave emission due to r-modes in the glitching pulsar PSR J0537-6910

Controllable emergent spatial spin modulation in Sr2IrO4 by in situ shear strain

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

Dirac Monopoles in Rogue Waves

Effective Field Theory Analysis of CDMSlite Run 2 Data

First joint observation by the underground gravitational-wave detector, KAGRA, with GEO600

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

Improving sensitivity of the ARIANNA detector by rejecting thermal noise with deep learning

Measurement of the Coherent Elastic Neutrino-Nucleus Scattering Cross Section on CsI by COHERENT

Measuring the Polarization Reconstruction Resolution of the ARIANNA Neutrino Detector with Cosmic Rays

Monitoring the SNS basement neutron background with the MARS detector

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

Physics Opportunities in the ORNL Spallation Neutron Source Second Target Station Era

Scintillator Tile Batch Test of CEPC AHCAL

Search for anisotropic gravitational-wave backgrounds using data from Advanced LIGO and Advanced Virgo&#39;s first three observing runs

Search for continuous gravitational wave emission from the Milky Way center in O3 LIGO--Virgo data

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

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

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in the Second and Third LIGO-Virgo Observing Runs

Simulating the neutrino flux from the Spallation Neutron Source for the COHERENT experiment

Snowmass 2021 White Paper Instrumentation Frontier 05 -- White Paper 1: MPGDs: Recent advances and current R&D

The COHERENT Experimental Program

The Proton Spin Structure Function $g_2$ and Generalized Polarizabilities in the Strong QCD Regime

Triboelectric Backgrounds to radio-based UHE Neutrino Exeperiments

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

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

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

Characterization of a CMOS pixel sensor for charged particle tracking

Development of Planar P-type Point Contact Germanium Detectors for Low-Mass Dark Matter Searches

First Measurement of Coherent Elastic Neutrino-Nucleus Scattering on Argon

Modeling Homophily in Dynamic Networks with Application to HIV Molecular Surveillance

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

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

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

Search for subsolar-mass binaries in the first half of Advanced LIGO and Virgo&#39;s third observing run

Semiclassical Trajectory Perspective of Glory Rescattering in Strong-field Photoelectron Holography

Timing performance of a multi-pad PICOSEC-Micromegas detector prototype

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

A TOSA/ROSA-Based Optical Transmitter (MTx+)/Transceiver (MTRx+) for High-Energy Physics Experiments

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

Binary Black Hole Population Properties Inferred from the First and Second Observing Runs of Advanced LIGO and Advanced Virgo

Characterization of High-Purity Germanium Detectors with Amorphous Germanium Contacts in Cryogenic Liquids

COHERENT Collaboration data release from the first detection of coherent elastic neutrino-nucleus scattering on argon

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

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

Reconstruction of atmospheric neutrinos in DUNE's horizontal-drift far-detector module

Search for anisotropic gravitational-wave backgrounds using data from Advanced LIGO and Advanced Virgo's first three observing runs

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

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

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

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