Researcher profile

D. Li

D. Li contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

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

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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

Building this graph slice

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

Published work

44 published item(s)

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

Evidence for particle acceleration approaching PeV energies in the W51 complex

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2025arXiv

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

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

preprint2025arXiv

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

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

preprint2025arXiv

The Nearby Evolved Stars Survey III: First data release of JCMT CO-line observations

Low- to intermediate-mass ($\sim$0.8$-$8 M$_\odot$) evolved stars contribute significantly to the chemical enrichment of the interstellar medium in the local Universe, making accurate mass-return estimates in their final stages crucial. The Nearby Evolved Stars Survey (NESS) is a large multi-telescope project targeting a volume-limited sample of $\sim$850 stars within 3 kpc in order to derive the dust and gas return rates in the Solar Neighbourhood, and to constrain the physics underlying these processes. We present an initial analysis of the CO-line observations, including detection statistics, carbon isotopic ratios, initial mass-loss rates, and gas-to-dust ratios. We describe a new data reduction pipeline to analyse the available NESS CO data from the JCMT, measuring line parameters and calculating empirical gas mass-loss rates. We present the first release of the available data on 485 sources, one of the largest homogeneous samples of CO data to date. Comparison with a large literature sample finds that high mass-loss rate and especially carbon-rich sources are over-represented in literature, while NESS is probing significantly more sources at low mass-loss rates, detecting 59 sources in CO for the first time and providing useful upper limits. CO line detection rates are 81% for the CO (2--1) line and 75% for CO (3--2). The majority (82%) of detected lines conform to the expected soft parabola shape, while eleven sources show a double wind. Calculated mass-loss rates show power-law relations with both the dust-production rates and expansion velocities up to $\sim 5 \times 10^{-6}$~\msunyr. Median gas-to-dust ratios of 250 and 680 are found for oxygen-rich and carbon-rich sources, respectively. Our analysis of CO observations in this first data release highlights the importance of our volume-limited approach in characterizing the local AGB population as a whole.

preprint2022arXiv

A repeating fast radio burst associated with a persistent radio source

The dispersive sweep of fast radio bursts (FRBs) has been used to probe the ionized baryon content of the intergalactic medium, which is assumed to dominate the total extragalactic dispersion. While the host galaxy contributions to dispersion measure (DM) appear to be small for most FRBs, in at least one case there is evidence for an extreme magneto-ionic local environment and a compact persistent radio source. Here we report the detection and localization of the repeating FRB 20190520B, which is co-located with a compact, persistent radio source and associated with a dwarf host galaxy of high specific star formation rate at a redshift $z=0.241\pm0.001$. The estimated host galaxy DM $\approx 903^{+72}_{-111}$ pc cm$^{-3}$, nearly an order of magnitude higher than the average of FRB host galaxies, far exceeds the DM contribution of the intergalactic medium. Caution is thus warranted in inferring redshifts for FRBs without accurate host galaxy identifications.

preprint2022arXiv

Arecibo and FAST Timing Follow-up of twelve Millisecond Pulsars Discovered in Commensal Radio Astronomy FAST Survey

We report the phase-connected timing ephemeris, polarization pulse profiles, Faraday rotation measurements, and Rotating-Vector-Model (RVM) fitting results of twelve millisecond pulsars (MSPs) discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) in the Commensal radio Astronomy FAST survey (CRAFTS). The timing campaigns were carried out with FAST and Arecibo over three years. Eleven of the twelve pulsars are in neutron star - white dwarf binary systems, with orbital periods between 2.4 and 100 d. Ten of them have spin periods, companion masses, and orbital eccentricities that are consistent with the theoretical expectations for MSP - Helium white dwarf (He WD) systems. The last binary pulsar (PSR J1912$-$0952) has a significantly smaller spin frequency and a smaller companion mass, the latter could be caused by a low orbital inclination for the system. Its orbital period of 29 days is well within the range of orbital periods where some MSP - He WD systems have shown anomalous eccentricities, however, the eccentricity of PSR J1912$-$0952 is typical of what one finds for the remaining MSP - He WD systems.

preprint2022arXiv

Detection of strong scattering close to the eclipse region of PSR B1957+20

We present the first measurement of pulse scattering close to the eclipse region of PSR B1957+20, which is in a compact binary system with a low-mass star. We measured pulse scattering time-scales up to 0.2 ms close to the eclipse and showed that it scales with the dispersion measure (DM) excess roughly as $τ\proptoΔ{\rm DM}^{2}$. Our observations provide the first evidence of strong scattering due to multi-path propagation effects in the eclipsing material. We show that Kolmogorov turbulence in the eclipsing material with an inner scale of $\sim100$ m and an outer scale of the size of the eclipse region can naturally explain the observation. Our results show that the eclipsing material in such systems can be highly turbulent and suggest that scattering is one of the main eclipsing mechanisms at around 1.4 GHz.

preprint2022arXiv

European Strategy for Particle Physics -- Accelerator R&D Roadmap

The 2020 update of the European Strategy for Particle Physics emphasised the importance of an intensified and well-coordinated programme of accelerator R&D, supporting the design and delivery of future particle accelerators in a timely, affordable and sustainable way. This report sets out a roadmap for European accelerator R&D for the next five to ten years, covering five topical areas identified in the Strategy update. The R&D objectives include: improvement of the performance and cost-performance of magnet and radio-frequency acceleration systems; investigations of the potential of laser / plasma acceleration and energy-recovery linac techniques; and development of new concepts for muon beams and muon colliders. The goal of the roadmap is to document the collective view of the field on the next steps for the R&D programme, and to provide the evidence base to support subsequent decisions on prioritisation, resourcing and implementation.

preprint2022arXiv

FAST observations of an extremely active episode of FRB 20201124A: I. Burst morphology

We report the properties of more than 600 bursts (including cluster-bursts) detected from the repeating fast radio burst (FRB) source FRB 20201124A with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) during an extremely active episode on UTC September 25-28, 2021, in a series of four papers. The observations were carried out in the band of 1.0 - 1.5 GHz by using the center beam of the L-band 19-beam receiver. We monitored the source in sixteen 1-hour sessions and one 3-hour session spanning 23 days. All the bursts were detected during the first four days. In this first paper of the series, we perform a detailed morphological study of 624 bursts using the 2-dimensional frequency-time ``waterfall&#39;&#39; plots, with a burst (or cluster-burst) defined as an emission episode during which the adjacent emission peaks have a separation shorter than 400 ms. The duration of a burst is therefore always longer than 1 ms, with the longest up to more than 120 ms. The emission spectra of the sub-bursts are typically narrow within the observing band with a characteristic width of $\sim$277 MHz. The center frequency distribution has a dominant peak at about 1091.9 MHz and a secondary weak peak around 1327.9 MHz. Most bursts show a frequency-downward-drifting pattern. Based on the drifting patterns, we classify the bursts into five main categories: downward drifting (263) bursts, upward drifting (3) bursts, complex (203), no drifting (35) bursts, and no evidence for drifting (121) bursts. Subtypes are introduced based on the emission frequency range in the band (low, middle, high and wide) as well as the number of components in one burst (1, 2, or multiple). We measured a varying scintillation bandwidth from about 0.5 MHz at 1.0 GHz to 1.4 MHz at 1.5 GHz with a spectral index of 3.0.

preprint2022arXiv

Magnetic Field Measurement and Analysis for the Muon g-2 Experiment at Fermilab

The Fermi National Accelerator Laboratory has measured the anomalous precession frequency $a^{}_μ= (g^{}_μ-2)/2$ of the muon to a combined precision of 0.46 parts per million with data collected during its first physics run in 2018. This paper documents the measurement of the magnetic field in the muon storage ring. The magnetic field is monitored by nuclear magnetic resonance systems and calibrated in terms of the equivalent proton spin precession frequency in a spherical water sample at 34.7$^\circ$C. The measured field is weighted by the muon distribution resulting in $\tildeω&#39;^{}_p$, the denominator in the ratio $ω^{}_a$/$\tildeω&#39;^{}_p$ that together with known fundamental constants yields $a^{}_μ$. The reported uncertainty on $\tildeω&#39;^{}_p$ for the Run-1 data set is 114 ppb consisting of uncertainty contributions from frequency extraction, calibration, mapping, tracking, and averaging of 56 ppb, and contributions from fast transient fields of 99 ppb.

preprint2022arXiv

Quantifying the Magnetic Structure of a Coronal Shock Producing a Type II Radio Burst

Type II radio bursts are thought to be produced by shock waves in the solar atmosphere. However, what magnetic conditions are needed for the generation of type II radio bursts is still a puzzling issue. Here, we quantify the magnetic structure of a coronal shock associated with a type II radio burst. Based on the multi-perspective extreme-ultraviolet observations, we reconstruct the three-dimensional (3D) shock surface. By using a magnetic field extrapolation model, we then derive the orientation of the magnetic field relative to the normal of the shock front ($θ_{\rm Bn}$) and Alfvén Mach number ($M_A$) on the shock front. Combining the radio observations from Nancay Radio Heliograph, we obtain the source region of the type II radio burst on the shock front. It is found that the radio burst is generated by a shock with $M_A \gtrsim 1.5$ and a bimodal distribution of $θ_{Bn}$. We also use the Rankine-Hugoniot relations to quantify the properties of the shock downstream. Our results provide a quantitative 3D magnetic structure condition of a coronal shock that produces a type II radio burst.

preprint2022arXiv

The discovery of a rotating radio transient J1918$-$0449 with intriguing emission properties with the five hundred meter aperture spherical radio telescope

In this study, we report on a detailed single pulse analysis of the radio emission from a rotating radio transient (RRAT) J1918$-$0449 which is the first RRAT discovered with the five hundred meter aperture spherical radio telescope (FAST). The sensitive observations were carried out on 30 April 2021 using the FAST with a central frequency of 1250 MHz and a short time resolution of 49.152 $μ$s, which forms a reliable basis to probe single pulse emission properties in detail. The source was successively observed for around 2 hours. A total of 83 dispersed bursts with significance above 6$σ$ are detected over 1.8 hours. The source&#39;s DM and rotational period are determined to be 116.1$\pm$0.4 \pcm \ and 2479.21$\pm$0.03 ms, respectively. The share of registered pulses from the total number of observed period is 3.12\%. No underlying emission is detected in the averaged off pulse profile. For bursts with fluence larger than 10 Jy ms, the pulse energy follows a power-law distribution with an index of $-3.1\pm0.4$, suggesting the existence of bright pulse emission. We find that the distribution of time between subsequent pulses is consistent with a stationary Poisson process and find no evidence of clustering over the 1.8 h observations, giving a mean burst rate of one burst every 66 s. Close inspection of the detected bright pulses reveals that 21 pulses exhibit well-defined quasi-periodicities. The subpulse drifting is present in non-successive rotations with periodicity of $2.51\pm0.06$ periods. Finally, possible physical mechanisms are discussed.

preprint2022arXiv

The Large Dispersion and Scattering of FRB 20190520B are Dominated by the Host Galaxy

The repeating FRB 20190520B is localized to a galaxy at $z=0.241$, much closer than expected given its dispersion measure $\rm DM=1205\pm4\ pc\ cm^{-3}$. Here we assess implications of the large DM and scattering observed from FRB 20190520B for the host galaxy&#39;s plasma properties. A sample of 75 bursts detected with the Five-hundred-meter Aperture Spherical radio Telescope shows scattering on two scales: a mean temporal delay $τ(1.41\ {\rm GHz})=10.9\pm1.5$ ms, which is attributed to the host galaxy, and a mean scintillation bandwidth $ν_{\rm d}(1.41\ {\rm GHz})=0.21\pm0.01$ MHz, which is attributed to the Milky Way. Balmer line measurements for the host imply an H$α$ emission measure (galaxy frame) $\rm EM_s=620$ pc cm$^{-6} \times (T/10^4 {\rm K})^{0.9}$, implying $\rm DM_{\rm Hα}$ of order the value inferred from the FRB DM budget, $\rm DM_h=1121^{+89}_{-138}$ pc cm$^{-3}$ for plasma temperatures greater than the typical value $10^4$ K. Combining $τ$ and $\rm DM_h$ yields a nominal constraint on the scattering amplification from the host galaxy $\tilde{F} G=1.5^{+0.8}_{-0.3}$ (pc$^2$ km)$^{-1/3}$, where $\tilde{F}$ describes turbulent density fluctuations and $G$ represents the geometric leverage to scattering that depends on the location of the scattering material. For a two-screen scattering geometry where $τ$ arises from the host galaxy and $Δν_{\rm d}$ from the Milky Way, the implied distance between the FRB source and dominant scattering material is $\lesssim100$ pc. The host galaxy scattering and DM contributions support a novel technique for estimating FRB redshifts using the $τ-\rm DM$ relation, and are consistent with previous findings that scattering of localized FRBs is largely dominated by plasma within host galaxies and the Milky Way.

preprint2022arXiv

Transverse Coronal-Loop Oscillations Induced by the Non-radial Eruption of a Magnetic Flux Rope

We investigate the transverse coronal-loop oscillations induced by the eruption of a prominence-carrying flux rope on 7 December 2012. The flux rope originating from NOAA Active Region (AR) 11621 was observed in EUV wavelengths by the SDO/AIA and in H$α$ line center by the ground-based telescope at the BBSO. The early evolution of the flux rope is divided into two steps: a slow rise phase at a speed of $\approx$230\,km\,s$^{-1}$ and a fast rise phase at a speed of $\approx$706\,km\,s$^{-1}$. The eruption generates a C5.8 flare and the onset of the fast rise is consistent with the HXR peak time of the flare. The embedded prominence has a lower speed of $\approx$452\,km\,s$^{-1}$. During the early eruption of the flux rope, the nearby coronal loops are disturbed and experience independent kink-mode oscillations in the horizontal and vertical directions. The oscillation in the horizontal direction has an initial amplitude of $\approx$3.1\,Mm, a period of $\approx$294\,seconds, and a damping time of $\approx$645\,seconds. It is most striking in 171\,Å and lasts for three to four cycles. The oscillations in the vertical directions are observed mainly in 171, 193, and 211\,Å. The initial amplitudes lie in the range of 3.4\,--\,5.2\,Mm, with an average value of 4.5\,Mm. The periods are between 407\,seconds and 441\,seconds, with an average value of 423\,seconds. The oscillations are damping and last for nearly four cycles. The damping times lie in the range of 570\,--\,1012\,seconds, with an average value of 741\,seconds. Assuming a semi-circular shape of the vertically oscillating loops, we calculate the loop lengths according to their heights. Using the observed periods, we carry out coronal seismology and estimate the internal Alfvén speeds (988\,--\,1145\,km\,s$^{-1}$) and the magnetic-field strengths (12\,--\,43\,G) of the oscillating loops.

preprint2021arXiv

A Demonstration of Improved Constraints on Primordial Gravitational Waves with Delensing

We present a constraint on the tensor-to-scalar ratio, $r$, derived from measurements of cosmic microwave background (CMB) polarization $B$-modes with &#34;delensing,&#34; whereby the uncertainty on $r$ contributed by the sample variance of the gravitational lensing $B$-modes is reduced by cross-correlating against a lensing $B$-mode template. This template is constructed by combining an estimate of the polarized CMB with a tracer of the projected large-scale structure. The large-scale-structure tracer used is a map of the cosmic infrared background derived from Planck satellite data, while the polarized CMB map comes from a combination of South Pole Telescope, BICEP/Keck, and Planck data. We expand the BICEP/Keck likelihood analysis framework to accept a lensing template and apply it to the BICEP/Keck data set collected through 2014 using the same parametric foreground modelling as in the previous analysis. From simulations, we find that the uncertainty on $r$ is reduced by $\sim10\%$, from $σ(r)$= 0.024 to 0.022, which can be compared with a $\sim26\%$ reduction obtained when using a perfect lensing template. Applying the technique to the real data, the constraint on $r$ is improved from $r_{0.05} < 0.090$ to $r_{0.05} < 0.082$ (95\% C.L.). This is the first demonstration of improvement in an $r$ constraint through delensing.

preprint2021arXiv

Long and Short Fast Radio Bursts are Different from Repeating and Non-repeating Transients

We collect 133 Fast Radio Bursts (FRBs), including 110 non-repeating and 23 repeating ones, and systematically investigate their observational properties. To check the frequency dependence of FRB classifications, we define our samples with a central frequency below/above 1GHz as subsample I/II. We find that there is a clear bimodal distribution of pulse width for the subsample I. And If we classify FRBs into short FRBs (\emph{s}FRBs) ($<$100 ms) and long FRBs (\emph{l}FRBs) ($>$100 ms) as done for short and long Gamma-Ray Bursts (GRBs), the \emph{s}FRBs at higher central frequency are commonly shorter than those at lower central frequency not only for non-repeating but also repeating \emph{s}FRBs. Secondly, we find that fluence and peak flux density are correlated with a power law relation of $F \varpropto S{^γ_{p,obs}}$ for both \emph{s}FRBs and \emph{l}FRBs whose distributions are obviously different. Thirdly, the \emph{l}FRBs with isotropic energies ranging from $10^{42}$ to $10^{44}$ erg are more energetic than the \emph{s}FRBs in the $F- DM_{EX}$ plane, indicating that they are two representative types. Finally, it is interestingly note that the peak flux density behaves an independence on the redshift when the distance of the FRBs becomes far enough, which is similar to the scenario of peak flux evolving with redshift in the field of GRBs. We predict that fainter FRBs at higher redshift of $z>2$ can be successfully detected by FAST and SKA in the near future.

preprint2021arXiv

The Nearby Evolved Stars Survey II: Constructing a volume-limited sample and first results from the James Clerk Maxwell Telescope

The Nearby Evolved Stars Survey (NESS) is a volume-complete sample of $\sim$850 Galactic evolved stars within 3\,kpc at (sub-)mm wavelengths, observed in the CO $J = $ (2$-$1) and (3$-$2) rotational lines, and the sub-mm continuum, using the James Clark Maxwell Telescope and Atacama Pathfinder Experiment. NESS consists of five tiers, based on distances and dust-production rate (DPR). We define a new metric for estimating the distances to evolved stars and compare its results to \emph{Gaia} EDR3. Replicating other studies, the most-evolved, highly enshrouded objects in the Galactic Plane dominate the dust returned by our sources, and we initially estimate a total DPR of $4.7\times 10^{-5}$ M$_\odot$ yr$^{-1}$ from our sample. Our sub-mm fluxes are systematically higher and spectral indices are typically shallower than dust models typically predict. The 450/850 $μ$m spectral indices are consistent with the blackbody Rayleigh--Jeans regime, suggesting a large fraction of evolved stars have unexpectedly large envelopes of cold dust.

preprint2020arXiv

An extensive spectroscopic time series of three Wolf-Rayet stars -- II. A search for wind asymmetries in the dust-forming WC7 binary WR137

We present the results of a four-month, spectroscopic campaign of the Wolf-Rayet dust-making binary, WR137. We detect only small-amplitude, random variability in the CIII5696 emission line and its integrated quantities (radial velocity, equivalent width, skewness, kurtosis) that can be explained by stochastic clumps in the wind of the WC star. We find no evidence of large-scale, periodic variations often associated with Corotating Interaction Regions that could have explained the observed intrinsic continuum polarization of this star. Our moderately high-resolution and high signal-to-noise average Keck spectrum shows narrow double-peak emission profiles in the Halpha, Hbeta, Hgamma, HeII6678 and HeII5876 lines. These peaks have a stable blue-to-red intensity ratio with a mean of 0.997 and a root-mean-square of 0.004, commensurate with the noise level; no variability is found during the entire observing period. We suggest that these profiles arise in a decretion disk around the O9 companion, which is thus an O9e star. The characteristics of the profiles are compatible with those of other Be/Oe stars. The presence of this disk can explain the constant component of the continuum polarization of this system, for which the angle is perpendicular to the plane of the orbit, implying that the rotation axis of the O9e star is aligned with that of the orbit. It remains to be explained why the disk is so stable within the strong ultraviolet radiation field of the O star. We present a binary evolutionary scenario that is compatible with the current stellar and system parameters.

preprint2020arXiv

An Impedance-Modulated Code-Division Microwave SQUID Multiplexer

Large arrays of cryogenic detectors, including transition-edge sensors (TESs) or magnetic micro-calorimeters (MMCs), are needed for future experiments across a wide range of applications. Complexities in integration and cryogenic wiring have driven efforts to develop cryogenic readout technologies with large multiplexing factors while maintaining minimal readout noise. One such example is the microwave SQUID multiplexer ($μ$mux), which couples an incoming TES or magnetic calorimeter signal to a unique GHz-frequency resonance that is modulated in frequency. Here, we present a hybrid scheme combining the microwave SQUID multiplexer with code division multiplexing: the impedance-modulated code-division multiplexer (Z-CDM), which may enable an order of magnitude increase in multiplexing factor particularly for low-bandwidth signal applications.

preprint2020arXiv

An Improved Measurement of the Secondary Cosmic Microwave Background Anisotropies from the SPT-SZ + SPTpol Surveys

We report new measurements of millimeter-wave power spectra in the angular multipole range $2000 \le \ell \le 11,000$ (angular scales $5^\prime \gtrsim θ\gtrsim 1^\prime$). By adding 95 and 150\,GHz data from the low-noise 500 deg$^2$ SPTpol survey to the SPT-SZ three-frequency 2540 deg$^2$ survey, we substantially reduce the uncertainties in these bands. These power spectra include contributions from the primary cosmic microwave background, cosmic infrared background, radio galaxies, and thermal and kinematic Sunyaev-Zel&#39;dovich (SZ) effects. The data favor a thermal SZ (tSZ) power at 143\,GHz of $D^{\rm tSZ}_{3000} = 3.42 \pm 0.54~ μ{\rm K}^2$ and a kinematic SZ (kSZ) power of $D^{\rm kSZ}_{3000} = 3.0 \pm 1.0~ μ{\rm K}^2$. This is the first measurement of kSZ power at $\ge 3\,σ$. We study the implications of the measured kSZ power for the epoch of reionization, finding the duration of reionization to be $Δz_{re} = 1.0^{+1.6}_{-0.7}$ ($Δz_{re}< 4.1$ at 95% confidence), when combined with our previously published tSZ bispectrum measurement.

preprint2020arXiv

An in-depth investigation of 11 pulsars discovered by FAST

We present timing solutions and analyses of 11 pulsars discovered by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). These pulsars were discovered using an ultra-wide bandwidth receiver in drift-scan observations made during the commissioning phase of FAST, and were then confirmed and timed using the 64-m Parkes Radio Telescope. Each pulsar has been observed over a span of at least one year. Highlighted discoveries include PSR J0344-0901, which displays mode-changing behaviour and may belong to the class of so-called `swooshing&#39; pulsars (alongside PSRs B0919+06 and B1859+07); PSR J0803-0942, whose emission is almost completely linearly polarised; and PSRs J1900-0134 and J1945+1211, whose well defined polarisation angle curves place stringent constraints on their emission geometry. We further discuss the detectability of these pulsars by earlier surveys, and highlight lessons learned from our work in carrying out confirmation and monitoring observations of pulsars discovered by a highly sensitive telescope, many of which may be applicable to next-generation pulsar surveys. This paper marks one of the first major releases of FAST-discovered pulsars, and paves the way for future discoveries anticipated from the Commensal Radio Astronomy FAST Survey (CRAFTS).

preprint2020arXiv

Constraints on Cosmological Parameters from the 500 deg$^2$ SPTpol Lensing Power Spectrum

We present cosmological constraints based on the cosmic microwave background (CMB) lensing potential power spectrum measurement from the recent 500 deg$^2$ SPTpol survey, the most precise CMB lensing measurement from the ground to date. We fit a flat $Λ$CDM model to the reconstructed lensing power spectrum alone and in addition with other data sets: baryon acoustic oscillations (BAO) as well as primary CMB spectra from Planck and SPTpol. The cosmological constraints based on SPTpol and Planck lensing band powers are in good agreement when analysed alone and in combination with Planck full-sky primary CMB data. With weak priors on the baryon density and other parameters, the CMB lensing data alone provide a 4\% constraint on $σ_8Ω_m^{0.25} = 0.0593 \pm 0.025$.. Jointly fitting with BAO data, we find $σ_8=0.779 \pm 0.023$, $Ω_m = 0.368^{+0.032}_{-0.037}$, and $H_0 = 72.0^{+2.1}_{-2.5}\,\text{km}\,\text{s}^{-1}\,\text{Mpc}^{-1} $, up to $2\,σ$ away from the central values preferred by Planck lensing + BAO. However, we recover good agreement between SPTpol and Planck when restricting the analysis to similar scales. We also consider single-parameter extensions to the flat $Λ$CDM model. The SPTpol lensing spectrum constrains the spatial curvature to be $Ω_K = -0.0007 \pm 0.0025$ and the sum of the neutrino masses to be $\sum m_ν < 0.23$ eV at 95\% C.L. (with Planck primary CMB and BAO data), in good agreement with the Planck lensing results. With the differences in the $S/N$ of the lensing modes and the angular scales covered in the lensing spectra, this analysis represents an important independent check on the full-sky Planck lensing measurement.

preprint2020arXiv

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

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

preprint2020arXiv

Fractional Polarisation of Extragalactic Sources in the 500-square-degree SPTpol Survey

We study the polarisation properties of extragalactic sources at 95 and 150 GHz in the SPTpol 500 deg$^2$ survey. We estimate the polarised power by stacking maps at known source positions, and correct for noise bias by subtracting the mean polarised power at random positions in the maps. We show that the method is unbiased using a set of simulated maps with similar noise properties to the real SPTpol maps. We find a flux-weighted mean-squared polarisation fraction $\langle p^2 \rangle= [8.9\pm1.1] \times 10^{-4}$ at 95 GHz and $[6.9\pm1.1] \times 10^{-4}$ at 150~GHz for the full sample. This is consistent with the values obtained for a sub-sample of active galactic nuclei. For dusty sources, we find 95 per cent upper limits of $\langle p^2 \rangle_{\rm 95}<16.9 \times 10^{-3}$ and $\langle p^2 \rangle_{\rm 150}<2.6 \times 10^{-3}$. We find no evidence that the polarisation fraction depends on the source flux or observing frequency. The 1-$σ$ upper limit on measured mean squared polarisation fraction at 150 GHz implies that extragalactic foregrounds will be subdominant to the CMB E and B mode polarisation power spectra out to at least $\ell\lesssim5700$ ($\ell\lesssim4700$) and $\ell\lesssim5300$ ($\ell\lesssim3600$), respectively at 95 (150) GHz.

preprint2020arXiv

Galaxy Clusters Selected via the Sunyaev-Zel&#39;dovich Effect in the SPTpol 100-Square-Degree Survey

We present a catalog of galaxy cluster candidates detected in 100 square degrees surveyed with the SPTpol receiver on the South Pole Telescope. The catalog contains 89 candidates detected with a signal-to-noise ratio greater than 4.6. The candidates are selected using the Sunyaev-Zel&#39;dovich effect at 95 and 150 GHz. Using both space- and ground-based optical and infrared telescopes, we have confirmed 81 candidates as galaxy clusters. We use these follow-up images and archival images to estimate photometric redshifts for 66 galaxy clusters and spectroscopic observations to obtain redshifts for 13 systems. An additional 2 galaxy clusters are confirmed using the overdensity of near-infrared galaxies only, and are presented without redshifts. We find that 15 candidates (18% of the total sample) are at redshift of $z \geq 1.0$, with a maximum confirmed redshift of $z_{\rm{max}} = 1.38 \pm 0.10$. We expect this catalog to contain every galaxy cluster with $M_{500c} > 2.6 \times 10^{14} M_\odot h^{-1}_{70}$ and $z > 0.25$ in the survey area. The mass threshold is approximately constant above $z = 0.25$, and the complete catalog has a median mass of approximately $ M_{500c} = 2.7 \times 10^{14} M_\odot h^{-1}_{70}$. Compared to previous SPT works, the increased depth of the millimeter-wave data (11.2 and 6.5 $μ$K-arcmin at 95 and 150 GHz, respectively) makes it possible to find more galaxy clusters at high redshift and lower mass.

preprint2020arXiv

IRIS Si IV Line Profiles at Flare Ribbons as Indications of Chromospheric Condensation

We present temporal variations of the Si IV line profiles at the flare ribbons in three solar flares observed by the Interface Region Imaging Spectrograph (IRIS). In the M1.1 flare on 2014 September 6 and the X1.6 flare on 2014 September 10, the Si IV line profiles evolve from wholly redshifted to red-wing enhanced with the flare development. However, in the B1.8 flare on 2016 December 2, the Si IV line profiles are wholly redshifted throughout the flare evolution. We fit the wholly redshifted line profiles with a single Gaussian function but the red-asymmetric ones with a double Gaussian function to deduce the corresponding Doppler velocities. In addition, we find that hard X-ray emission above 25 keV shows up in the two large flares, implying a nonthermal electron beam heating. In the microflare, there only appears weak hard X-ray emission up to 12 keV, indicative of a thermal heating mostly. We interpret the redshifts or red asymmetries of the Si IV line at the ribbons in the three flares as spectral manifestations of chromospheric condensation. We propose that whether the line appears to be wholly redshifted or red-asymmetric depends on the heating mechanisms and also on the propagation of the condensation.

preprint2020arXiv

Measurements of B-mode Polarization of the Cosmic Microwave Background from 500 Square Degrees of SPTpol Data

We report a B-mode power spectrum measurement from the cosmic microwave background (CMB) polarization anisotropy observations made using the SPTpol instrument on the South Pole Telescope. This work uses 500 deg$^2$ of SPTpol data, a five-fold increase over the last SPTpol B-mode release. As a result, the bandpower uncertainties have been reduced by more than a factor of two, and the measurement extends to lower multipoles: $52 < \ell < 2301$. Data from both 95 and 150 GHz are used, allowing for three cross-spectra: 95 GHz x 95 GHz, 95 GHz x 150 GHz, and 150 GHz x 150 GHz. B-mode power is detected at very high significance; we find $P(BB < 0) = 5.8 \times 10^{-71}$, corresponding to a $18.1 σ$ detection of power. An upper limit is set on the tensor-to-scalar ratio, $r < 0.44$ at 95% confidence (the expected $1 σ$ constraint on $r$ given the measurement uncertainties is 0.22). We find the measured B-mode power is consistent with the Planck best-fit $Λ$CDM model predictions. Scaling the predicted lensing B-mode power in this model by a factor Alens, the data prefer Alens = $1.17 \pm 0.13$. These data are currently the most precise measurements of B-mode power at $\ell > 320$.

preprint2020arXiv

Parkes transient events: I. Database of single pulses, initial results and missing FRBs

A large number of observations from the Parkes 64\,m-diameter radio telescope, recorded with high time resolution, are publicly available. We have re-processed all of the observations obtained during the first four years (from 1997 to 2001) of the Parkes Multibeam receiver system in order to identify transient events and have built a database that records the 568,736,756 pulse candidates generated during this search. We have discovered a new fast radio burst (FRB), FRB~010305, with a dispersion measure of 350$\pm$5\,\,cm$^{-3}\,$pc and explored why so few FRBs have been discovered in data prior to 2001. After accounting for the dispersion smearing across the channel bandwidth and the sky regions surveyed, the number of FRBs is found to be consistent with model predictions. We also present five single pulse candidates from unknown sources, but with Galactic dispersion measures. We extract a diverse range of sources from the database, which can be used, for example, as a training set of data for new software being developed to search for FRBs in the presence of radio frequency interference.

preprint2020arXiv

Transverse coronal loop oscillations excited by homologous circular-ribbon flares

We report our multiwavelength observations of two homologous circular-ribbon flares (CRFs) in active region 11991 on 2014 March 5, focusing on the transverse oscillations of an extreme-ultraviolet (EUV) loop excited by the flares. The transverse oscillations are of fast standing kink-mode. The first-stage oscillation triggered by the C2.8 flare is decayless with lower amplitudes (310$-$510 km). The periods (115$-$118 s) in different wavelengths are nearly the same, indicating coherent oscillations. The magnetic field of the loop is estimated to be 65$-$78 G. The second-stage oscillation triggered by the M1.0 flare is decaying with larger amplitudes (1250$-$1280 km). The periods decreases from 117 s in 211 Å to 70 s in 171 Å, implying a decrease of loop length or an implosion after a gradual expansion. The damping time, being 147$-$315 s, increases with the period, so that the values of $τ/P$ are close to each other in different wavelengths. The thickness of the inhomogeneous layer is estimated to be $\sim$0\farcs45 under the assumption of resonant absorption. This is the first observation of the excitation of two kink-mode loop oscillations by two sympathetic flares. The results are important for understanding of the excitation of kink oscillations of coronal loops and hence the energy balance in the solar corona. Our findings also validate the prevalence of significantly amplified amplitudes of oscillations by successive drivers.

preprint2019arXiv

Emergence of superconductivity in strongly correlated hole-dominated Fe1-xSe

Here we establish a more complete phase diagram for FeSe system, based on experimental results of nonstoichiometric Fe1-xSe single crystals that we have developed recently, as well as nearly stoichiometric FeSe single crystals. The electronic correlation is found to be strongly enhanced in hole-dominated Fe1-xSe, as compared with electron-dominated FeSe, from the magnetic susceptibility and electrical transport measurements in the normal state. A superconducting dome is found to emerge starting from the strongly correlated hole-dominated regime with electron doping, while the tetragonal-orthorhombic phase transition at ~90 K is observed only at higher electron-doping levels in the electron-dominated regime.

preprint2019arXiv

Probing the cold magnetized Universe with SPICA-POL (B-BOP)

SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A&#39; concept study as one of the three remaining candidates for ESA&#39;s fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars.

preprint2019arXiv

Quasi-periodic pulsations of gamma-ray emissions from a solar flare on 2017 September 06

We investigate quasi-periodic pulsations (QPPs) of high-energy nonthermal emissions from an X9.3 flare (SOL2017-Sep-06T11:53), the most powerful flare since the beginning of solar cycle 24. The QPPs are identified as a series of regular and repeating peaks in the light curves in the gamma- and hard X-ray (HXR) channels recorded by the Konus-Wind, as well as the radio and microwave fluxes measured by the CALLISTO radio spectrograph during the impulsive phase. The periods are determined from the global wavelet and Fourier power spectra, as 24-30 s in the HXR and microwave channels which are associated with nonthermal electrons, and ~20 s in the gamma-ray band related to nonthermal ions. Both nonthermal electrons and ions may be accelerated by repetitive magnetic reconnection during the impulsive phase. However, we could not rule out other mechanisms such as the MHD oscillation in a sausage mode. The QPP detected in this study is useful for understanding the particle acceleration and dynamic process in solar flares and also bridging the gap between stellar and solar flares since the energy realm of the X9.3 solar flare is almost compared with a typical stellar flare.

preprint2019arXiv

The SPTpol Extended Cluster Survey

We describe the observations and resultant galaxy cluster catalog from the 2770 deg$^2$ SPTpol Extended Cluster Survey (SPT-ECS). Clusters are identified via the Sunyaev-Zel&#39;dovich (SZ) effect, and confirmed with a combination of archival and targeted follow-up data, making particular use of data from the Dark Energy Survey (DES). With incomplete followup we have confirmed as clusters 244 of 266 candidates at a detection significance $ξ\ge 5$ and an additional 204 systems at $4<ξ<5$. The confirmed sample has a median mass of $M_{500c} \sim {4.4 \times 10^{14} M_\odot h_{70}^{-1}}$, a median redshift of $z=0.49$, and we have identified 44 strong gravitational lenses in the sample thus far. Radio data are used to characterize contamination to the SZ signal; the median contamination for confirmed clusters is predicted to be $\sim$1% of the SZ signal at the $ξ>4$ threshold, and $<4\%$ of clusters have a predicted contamination $>10\% $ of their measured SZ flux. We associate SZ-selected clusters, from both SPT-ECS and the SPT-SZ survey, with clusters from the DES redMaPPer sample, and find an offset distribution between the SZ center and central galaxy in general agreement with previous work, though with a larger fraction of clusters with significant offsets. Adopting a fixed Planck-like cosmology, we measure the optical richness-to-SZ-mass ($λ-M$) relation and find it to be 28% shallower than that from a weak-lensing analysis of the DES data---a difference significant at the 4 $σ$ level---with the relations intersecting at $λ=60$ . The SPT-ECS cluster sample will be particularly useful for studying the evolution of massive clusters and, in combination with DES lensing observations and the SPT-SZ cluster sample, will be an important component of future cosmological analyses.