We present here the analysis of giant micropulses from the Vela pulsar. A total of 4187 giant micropulses with peak flux density $>$2.5 Jy were detected during almost 4 hours of observations carried out with the Yunnan 40-m radio telescope at 6800 MHz. Nine of the giant micropulses arrived approximately 3 to 4 ms earlier than the peak of average pulse profile, longer than that at lower frequencies. The remaining giant micropulses were clustered into three distributions which correspond to three main emission regions, including four occurring on the trailing edge of averaged profile.

We use K-giant stars selected from the LAMOST DR5 to study the variation of the rotational velocity of the galactic halo at different space positions. Modelling the rotational velocity distribution with both the halo and disk components, we find that the rotational velocity of the halo population decreases almost linearly with increasing vertical distance to the galactic disk plane, $Z$, at fixed galactocentric radius, $R$. The samples are separated into two parts with $6<R<12$ kpc and $12<R<20$ kpc. We derive that the decreasing rates along $Z$ for the two subsamples are $-3.07\pm0.63$ and $-1.89\pm0.37$ km s$^{-1}$ kpc$^{-1}$, respectively. Compared with the TNG simulations, we suggest that this trend is probably caused by the interaction between the disk and halo. The results from the simulations show that only the oblate halo can provide a decreasing rotational velocity with an increasing $Z$. This indicates that the Galactic halo is oblate with galactocentric radius $R<20$ kpc. On the other hand, the flaring of the disk component (mainly the thick disk) is clearly traced by this study, with $R$ between 12 and 20 kpc, the disk can vertically extend to $6\sim10$ kpc above the disk plane. What is more interesting is that, we find the Gaia-Enceladus-Sausage (GES) component has a significant contribution only in the halo with $R<12$ kpc, i.e. a fraction of 23$-$47\%. While in the outer subsample, the contribution is too low to be well constrained.

We have performed microwave diagnostics of the magnetic field strengths in solar flare loops based on the theory of gyrosynchrotron emission. From Nobeyama Radioheliograph observations of three flare events at 17 and 34 GHz, we obtained the degree of circular polarization and the spectral index of microwave flux density, which were then used to map the magnetic field strengths in post-flare loops. Our results show that the magnetic field strength typically decreases from ~800 G near the loop footpoints to ~100 G at a height of 10--25 Mm. Comparison of our results with magnetic field modeling using a flux rope insertion method is also discussed. Our study demonstrates the potential of microwave imaging observations, even at only two frequencies, in diagnosing the coronal magnetic field of flaring regions.

Observations of the time-dependent cosmic-ray Sun shadow have been proven as a valuable diagnostic for the assessment of solar magnetic field models. In this paper, seven years of IceCube data are compared to solar activity and solar magnetic field models. A quantitative comparison of solar magnetic field models with IceCube data on the event rate level is performed for the first time. Additionally, a first energy-dependent analysis is presented and compared to recent predictions. We use seven years of IceCube data for the Moon and the Sun and compare them to simulations on data rate level. The simulations are performed for the geometrical shadow hypothesis for the Moon and the Sun and for a cosmic-ray propagation model governed by the solar magnetic field for the case of the Sun. We find that a linearly decreasing relationship between Sun shadow strength and solar activity is preferred over a constant relationship at the 6.4sigma level. We test two commonly used models of the coronal magnetic field, both combined with a Parker spiral, by modeling cosmic-ray propagation in the solar magnetic field. Both models predict a weakening of the shadow in times of high solar activity as it is also visible in the data. We find tensions with the data on the order of $3\sigma$ for both models, assuming only statistical uncertainties. The magnetic field model CSSS fits the data slightly better than the PFSS model. This is generally consistent with what is found previously by the Tibet AS-gamma Experiment, a deviation of the data from the two models is, however, not significant at this point. Regarding the energy dependence of the Sun shadow, we find indications that the shadowing effect increases with energy during times of high solar activity, in agreement with theoretical predictions.

The authors of this white paper met on 16-17 January 2020 at the New Jersey Institute of Technology, Newark, NJ, for a 2-day workshop that brought together a group of heliophysicists, data providers, expert modelers, and computer/data scientists. Their objective was to discuss critical developments and prospects of the application of machine and/or deep learning techniques for data analysis, modeling and forecasting in Heliophysics, and to shape a strategy for further developments in the field. The workshop combined a set of plenary sessions featuring invited introductory talks interleaved with a set of open discussion sessions. The outcome of the discussion is encapsulated in this white paper that also features a top-level list of recommendations agreed by participants.

The He I 10830 \AA\ triplet is a very informative indicator of chromospheric activities as the helium is the second most abundant element in the solar atmosphere. Taking advantage of the high resolution of the 1.6 m Goode Solar Telescope (GST) at Big Bear Solar Observatory (BBSO), previous observations have shown clear evidence of the enhanced absorption, instead of typically-observed emission, for two M-class flares. In this study, we analyze the evolution of the He I 10830 10830 \AA\ emission in numerical models and compare it with observations. The models represent the RADYN simulation results obtained from the F-CHROMA database. We consider the models with the injected electron spectra parameters close to observational estimates for the 2013-August-17 flare event ($\delta=8$, $E_c = \{15,20\}$ keV, $F=\{1\times 10^{11}, 3\times{}10^{11}\}$ erg cm$^{-2}$) in detail, as well as other available models. The modeling results agree well with observations, in the sense of both the maximum intensity decrease (-17.1%, compared to the observed value of -13.7%) and the trend of temporal variation (initial absorption phase followed by the emission). All models demonstrate the increased number densities and decreased ratio of the upper and lower level populations of He I 10830 10830 \AA\ transition in the initial phase, which enhances the opacity and forms an absorption feature. Models suggest that the temperatures and free electron densities at heights of 1.3-1.5 Mm should be larger than $\sim 10^4$ K and $6\times 10^{11}$ cm$^{-3}$ thresholds for the line to start being in emission.

Filaments are very common physical phenomena on the Sun and are often taken as important proxies of solar magnetic activities. The study of filaments has become a hot topic in the space weather research. For a more comprehensive understanding of filaments, especially for an understanding of solar activities of multiple solar cycles, it is necessary to perform a combined multifeature analysis by constructing a data set of multiple solar cycle data. To achieve this goal, we constructed a centennial data set that covers the H$\alpha$ data from five observatories around the world. During the data set construction, we encountered varieties of problems, such as data fusion, accurate determination of the solar edge, classifying data by quality, dynamic threshold, and so on, which arose mainly due to multiple sources and a large time span of data. But fortunately, these problems were well solved. The data set includes seven types of data products and eight types of feature parameters with which we can implement the functions of data searching and statistical analyses. It has the characteristics of better continuity and highly complementary to space observation data, especially in the wavelengths not covered by space observations, and covers many solar cycles (including more than 60 yr of high-cadence data). We expect that this new comprehensive data set as well as the tools will help researchers to significantly speed up their search for features or events of interest, for either statistical or case study purposes, and possibly help them get a better and more comprehensive understanding of solar filament mechanisms.

The Tianlai Cylinder Pathfinder is a radio interferometer array designed to test techniques for 21 cm intensity mapping in the post-reionization Universe, with the ultimate aim of mapping the large scale structure and measuring cosmological parameters such as the dark energy equation of state. Each of its three parallel cylinder reflectors is oriented in the north-south direction, and the array has a large field of view. As the Earth rotates, the northern sky is observed by drift scanning. The array is located in Hongliuxia, a radio-quiet site in Xinjiang, and saw its first light in September 2016. In this first data analysis paper for the Tianlai cylinder array, we discuss the sub-system qualification tests, and present basic system performance obtained from preliminary analysis of the commissioning observations during 2016-2018. We show typical interferometric visibility data, from which we derive the actual beam profile in the east-west direction and the frequency band-pass response. We describe also the calibration process to determine the complex gains for the array elements, either using bright astronomical point sources, or an artificial on site calibrator source, and discuss the instrument response stability, crucial for transit interferometry. Based on this analysis, we find a system temperature of about 90 K, and we also estimate the sensitivity of the array.

GWAC will have been built an integrated FOV of 5,000 $degree^2$ and have already built 1,800 square $degree^2$. The limit magnitude of a 10-second exposure image in the moonless night is 16R. In each observation night, GWAC produces about 0.7TB of raw data, and the data processing pipeline generates millions of single frame alerts. We describe the GWAC Data Processing and Management System (GPMS), including hardware architecture, database, detection-filtering-validation of transient candidates, data archiving, and user interfaces for the check of transient and the monitor of the system. GPMS combines general technology and software in astronomy and computer field, and use some advanced technologies such as deep learning. Practical results show that GPMS can fully meet the scientific data processing requirement of GWAC. It can online accomplish the detection, filtering and validation of millions of transient candidates, and feedback the final results to the astronomer in real-time. During the observation from October of 2018 to December of 2019, we have already found 102 transients.

Adopting the Standard Halo Model (SHM) of an isotropic Maxwellian velocity distribution for dark matter (DM) particles in the Galaxy, the most stringent current constraints on their spin-dependent scattering cross-section with nucleons come from the IceCube neutrino observatory and the PICO-60 C$_3$F$_8$ superheated bubble chamber experiments. The former is sensitive to high energy neutrinos from the self-annihilation of DM particles captured in the Sun, while the latter looks for nuclear recoil events from DM scattering off nucleons. Although slower DM particles are more likely to be captured by the Sun, the faster ones are more likely to be detected by PICO. Recent N-body simulations suggest significant deviations from the SHM for the smooth halo component of the DM, while observations hint at a dominant fraction of the local DM being in substructures. We use the method of Ferrer et al. (2015) to exploit the complementarity between the two approaches and derive conservative constraints on DM-nucleon scattering. Our results constrain $\sigma_{\mathrm{SD}} \lesssim 3 \times 10^{-39} \mathrm{cm}^2$ (6 $ \times 10^{-38} \mathrm{cm}^2$) at $\gtrsim 90\%$ C.L. for a DM particle of mass 1~TeV annihilating into $\tau^+ \tau^-$ ($b\bar{b}$) with a local density of $\rho_{\mathrm{DM}} = 0.3~\mathrm{ GeV/cm}^3$. The constraints scale inversely with $\rho_{\mathrm{DM}}$ and are independent of the DM velocity distribution.

Fast radio bursts (FRBs) are short pulses observed in radio band from cosmological distances, some of which emit repeating bursts. The physical origins of these mysterious events have been subject to wide speculations and heated debates. One class of models invoke soft gamma-ray repeaters (SGRs), or magnetars, as the sources of FRBs. Magnetars are rotating neutron stars with extremely strong magnetic field and can sporadically emit bursts from X-ray ($\sim$keV) to soft gamma-ray ($\sim$sub-MeV) with duration from $10^{-2}$~s to $10^2$~s. However, even though some bright radio bursts have been observed from some magnetars, no FRB-like events had been detected to be associated with any magnetar burst, including one giant flare, and no radio burst has been associated with any X-ray event from any magnetar. Therefore, there is still no observational evidence for magnetar-FRB association up to today. Recently, a pair of FRB-like bursts (FRB~200428 hereafter) separated by 30 milliseconds (ms) were detected from the general direction of the Galactic magnetar SGR~J1935+2154. Here we report the detection of a non-thermal X-ray burst in the 1-250 keV energy band with the Insight-HXMT satellite, which we identify as emitted from SGR~J1935+2154. The burst showed two hard peaks with a separation of $\sim30$ ms, consistent with the separation between the two bursts in FRB~200428. The delay time between the double radio and X-ray peaks is $\sim8.57$ s, fully consistent with the dispersion delay of FRB~200428. We thus identify the non-thermal X-ray burst is associated with FRB~200428 whose high energy counterpart is the two hard peaks in X-ray. Our results suggest that the non-thermal X-ray burst and FRB~200428 share the same physical origin in an explosive event from SGR~J1935+2154.

In the era of big data astronomy, next generation telescopes and large sky surveys produce data sets at the TB or even PB level. Due to their large data volumes, these astronomical data sets are extremely difficult to transfer and analyze using personal computers or small clusters. In order to offer better access to data, data centers now generally provide online science platforms that enable analysis close to the data. The Chinese Virtual Observatory (China-VO) is one of the member projects in the International Virtual Observatory Alliance and it is dedicated to providing a research and education environment where globally distributed astronomy archives are simple to find, access, and interoperate. In this study, we summarize highlights of the work conducted at the China-VO, as well the experiences and lessons learned during the full life-cycle management of astronomical data. Finally, We discuss the challenges and future trends for astronomical science platforms.

We propose a new machine learning approach to Stokes inversion based on a convolutional neural network (CNN) and the Milne-Eddington (ME) method. The Stokes measurements used in this study were taken by the Near InfraRed Imaging Spectropolarimeter (NIRIS) on the 1.6 m Goode Solar Telescope (GST) at the Big Bear Solar Observatory. By learning the latent patterns in the training data prepared by the physics-based ME tool, the proposed CNN method is able to infer vector magnetic fields from the Stokes profiles of GST/NIRIS. Experimental results show that our CNN method produces smoother and cleaner magnetic maps than the widely used ME method. Furthermore, the CNN method is 4~6 times faster than the ME method, and is able to produce vector magnetic fields in near real-time, which is essential to space weather forecasting. Specifically, it takes ~50 seconds for the CNN method to process an image of 720 x 720 pixels comprising Stokes profiles of GST/NIRIS. Finally, the CNN-inferred results are highly correlated to the ME-calculated results and are closer to the ME's results with the Pearson product-moment correlation coefficient (PPMCC) being closer to 1 on average than those from other machine learning algorithms such as multiple support vector regression and multilayer perceptrons (MLP). In particular, the CNN method outperforms the current best machine learning method (MLP) by 2.6% on average in PPMCC according to our experimental study. Thus, the proposed physics-assisted deep learning-based CNN tool can be considered as an alternative, efficient method for Stokes inversion for high resolution polarimetric observations obtained by GST/NIRIS.

We consider an $f(Q,T)$ type gravity model in which the scalar non-metricity $Q_{\alpha \mu \nu}$ of the space-time is expressed in its standard Weyl form, and it is fully determined by a vector field $w_{\mu}$. The field equations of the theory are obtained under the assumption of the vanishing of the total scalar curvature, a condition which is added into the gravitational action via a Lagrange multiplier. The gravitational field equations are obtained from a variational principle, and they explicitly depend on the scalar nonmetricity and on the Lagrange multiplier. The covariant divergence of the matter energy-momentum tensor is also determined, and it follows that the nonmetricity-matter coupling leads to the nonconservation of the energy and momentum. The energy and momentum balance equations are explicitly calculated, and the expressions of the energy source term and of the extra force are found. We investigate the cosmological implications of the theory, and we obtain the cosmological evolution equations for a flat, homogeneous and isotropic geometry, which generalize the Friedmann equations of standard general relativity. We consider several cosmological models by imposing some simple functional forms of the function $f(Q,T)$, and we compare the predictions of the theory with the standard $\Lambda$CDM model.

The CsI detectors of the High Energy X-ray Telescope of the Hard X-ray Modulation Telescope (HXMT/CsI) can be used for gamma-ray all sky monitoring and searching for the electromagnetic counterpart of gravitational wave source. The instrumental responses are mainly obtained by Monte Carlo simulation with the Geant4 tool and the mass model of both the satellite and all the payloads, which is updated and tested with the Crab pulse emission in various incident directions. Both the Energy-Channel relationship and the energy resolution are calibrated in two working modes (Normal-Gain mode & Low-Gain Mode) with the different detection energy ranges. The simulative spectral analyses show that HXMT/CsI can constrain the spectral parameters much better in the high energy band than that in the low energy band. The joint spectral analyses are performed to ten bright GRBs observed simultaneously with HXMT/CsI and other instruments (Fermi/GBM, Swift/BAT, Konus-Wind), and the results show that the GRB flux given by HXMT/CsI is systematically higher by $7.0\pm8.8\%$ than those given by the other instruments. The HXMT/CsI-Fermi/GBM joint fittings also show that the high energy spectral parameter can be constrained much better as the HXMT/CsI data are used in the joint fittings.

In the present paper, it is assumed that there exist two dark matter particles: superheavy dark matter particles (SHDM), whose mass $\sim$ inflaton mass, and light fermion dark matter (DM) particles which are the ultrahigh energy (UHE) products of its decay. The Earth will be taken as a detector to directly search for the UHE DM particles. These upward-going particles, which pass through the Earth and air and interact with nuclei, can be detected by the fluorescence detectors (FD) of the Pierre Auger observatory (Auger), via fluorescent photons due to the development of an EAS. The event rates of UHE DM particles measured at the FD of Auger are evaluated in the incoming energy range between 1 EeV and 1 ZeV when the different lifetimes of decay of SHDM. And the UHE DM fluxes are estimated at the FD of Auger when the different lifetimes of decay of SHDM. Finally, according to the Auger data from 2008 to 2019, the upper limit for UHE DM fluxes is estimated at 90\% C.L. at the FD of Auger.

We report a simultaneous 44 and 95 GHz class I methanol maser survey toward 144 sources from the 95 GHz class I methanol maser catalog. The observations were made with the three telescopes of the Korean very long baseline interferometry network operating in single-dish mode. The detection rates are 89% at 44 GHz and 77% at 95 GHz. There are 106 new discoveries at 44 GHz. Comparing the previous 95 GHz detections with new observations of the same transitions made using the Purple Mountain Observatory 13.7 m radio telescope shows no clear evidence of variability on a timescale of six years. Emission from the 44 and 95 GHz transitions shows strong correlations in peak velocity, peak flux density, and integrated flux density, indicating that they are likely cospatial. We found that the peak flux density ratio Spk,95/Spk,44 decreases as the 44 GHz peak flux density increases. We found that some class I methanol masers in our sample could be associated with infrared dark clouds, while others are associated with H II regions, indicating that some sources occur at an early stage of high-mass star formation, while others are located toward more evolved sources.

We present the results of a ~250 square degrees CO mapping (+26d<l<+50d and -5d<b<+5d) toward the Aquila Rift region at a spatial resolution of ~50" and a grid spacing of 30". The high dynamic range CO maps with a spectral resolution of ~0.2km/s display highly structured molecular cloud (MC) morphologies with valuable velocity information, revealing complex spatial and dynamical features of the local molecular gas. In combination with the MWISP CO data and the Gaia DR2, distances of the main MC structures in the local ISM are well determined toward the Aquila Rift. We find that the total MC mass within 1 kpc is about >4.1x10^5 Msun in the whole region. In fact, the mass of the molecular gas is dominated by the W40 giant molecular cloud (GMC) at ~474 pc (~1.4x10^5 Msun) and the GMC complex G036.0+01.0 at ~560-670 pc (~2.0x10^5 Msun), while the MCs at ~220-260pc have gas masses of 10^2-10^3 Msun. Interestingly, an ~80pc long filamentary MC G044.0-02.5 at a distance of ~404 pc shows a systematic velocity gradient along and perpendicular to the major axis of the filament. The HI gas with the enhanced emission has the similar spatial morphologies and velocity features compared to the corresponding CO structure, indicating that the large-scale converging HI flows are probably responsible for the formation of the MC. Meanwhile, the long filamentary MC consists of many sub-filaments with the lengths ranging from ~0.5 pc to several pc, as well as prevalent networks of filaments in other large-scale local MCs.

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

The GPS1 catalog was released in 2017. It delivered precise proper motions for around 350 million sources across three-fourths of the sky down to a magnitude of $r\sim20$\,mag. In this study, we present GPS1+ the extension GPS1 catalog down to $r\sim22.5$\,mag, based on {\it Gaia} DR2, PS1, SDSS and 2MASS astrometry. The GPS1+ totally provides proper motions for $\sim$400 million sources with a characteristic systematic error of less than 0.1\masyr. This catalog is divided into two sub-samples, i.e., the primary and secondary parts. The primary $\sim$264 million sources have either or both of the {\it Gaia} and SDSS astrometry, with a typical precision of 2.0-5.0 \masyr. In this part, $\sim$160 million sources have {\it Gaia} proper motions, we provide another new proper motion for each of them by building a Bayesian model. Relative to {\it Gaia}'s values, the precision is improved by $\sim$0.1\,dex on average at the faint end; $\sim$50 million sources are the objects whose proper motions are missing in {\it Gaia} DR2, we provide their proper motion with a precision of $\sim$4.5\masyr; the remaining $\sim$54 million faint sources are beyond {\it Gaia} detecting capability, we provide their proper motions for the first time with a precision of 7.0 \masyr. However, the secondary $\sim$136 million sources only have PS1 astrometry, the average precision is worse than 15.0 \masyr. All the proper motions have been validated using QSOs and the existing {\it Gaia} proper motions. The catalog will be released on-line and available via the VO-TAP Service.

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

Significant progress in the classification of Fermi unassociated sources , has led to an increasing number of blazars are being found. The optical spectrum is effectively used to classify the blazars into two groups such as BL Lacs and flat spectrum radio quasars (FSRQs). However, the accurate classification of the blazars without optical spectrum information, i.e., blazars of uncertain type (BCUs), remains a significant challenge. In this paper, we present a principal component analysis (PCA) and machine learning hybrid blazars classification method. The method, based on the data from Fermi LAT 3FGL Catalog, first used the PCA to extract the primary features of the BCUs and then used a machine learning algorithm to further classify the BCUs. Experimental results indicate that the that the use of PCA algorithms significantly improved the classification. More importantly, comparison with the Fermi LAT 4FGL Catalog, which contains the spectral classification of those BCUs in the Fermi-LAT 3FGL Catalog, reveals that the proposed classification method in the study exhibits higher accuracy than currently established methods; specifically, 151 out of 171 BL Lacs and 19 out of 24 FSRQs are correctly classified.

In the present paper, it is assumed that there exist two species of dark matter: a heavy dark matter particle (HDM) with the mass of O(TeV) which is generated in early universe and a lighter dark matter particle (LDM) which is a relativistic product due to the decay of HDM. HDMs, captured by the earth, decay to high energy LDMs, and these particles can be measured by km$^3$ neutrino telescopes, like the IceCube detector. A $Z^{\prime}$ portal dark matter model is taken for LDMs to interact with nuclei via a neutral current interaction mediated by a heavy gauge boson $Z^{\prime}$. With the different lifetimes of decay of HDMs and Z$^{\prime}$ masses, the event rates of LDMs, measured by IceCube, are evaluated in the energy range between 1 TeV and 100 TeV. According to the IceCube data, the upper limit for LDM fluxes is estimated at 90\% C.L. at IceCube. Finally, it is proved that LDMs could be directly detected in the energy range betwen O(1TeV) and O(10TeV) at IceCube with $m_{Z^{\prime}} \lesssim 500 GeV$ and $\tau_{\phi} \lesssim 10^{21}$ s.

In order to precisely measure the cosmological parameters and answer the fundamental questions in cosmology, it is necessary to develop new, powerful cosmological probes, in addition to the proposed next-generation optical survey projects. The neutral hydrogen 21 cm sky surveys will provide a promising tool to study the late-universe evolution, helping shed light on the nature of dark energy. The Square Kilometre Array is the largest radio telescope in the world to be constructed in the near future, and it will push the 21 cm cosmology into a new era and greatly promote the development of cosmology in the forthcoming decades.

We present a survey of molecular outflows across the dark cloud complex in the Cygnus region, based on 46.75 deg^2 field of CO isotopologues data from Milky Way Imaging Scroll Painting (MWISP) survey. A supervised machine learning algorithm, Support Vector Machine (SVM), is introduced to accelerate our visual assessment of outflow features in the data cube of 12CO and 13CO J = 1-0 emission. A total of 130 outflow candidates are identified, of which 77 show bipolar structures and 118 are new detections. Spatially, these outflows are located inside dense molecular clouds and some of them are found in clusters or in elongated linear structures tracing the underlying gas filament morphology. Along the line of sight, 97, 31, and 2 candidates reside in the Local, Perseus, and Outer arm, respectively. Young stellar objects as outflow drivers are found near most outflows, while 36 candidates show no associated source. The clusters of outflows that we detect are inhomogeneous in their properties; nevertheless, we show that the outflows cannot inject turbulent energy on cloud scales. Instead, at best, they are restricted to affecting the so called "clump" and "core" scales, and this only on short (~0.3 Myr) estimated timescales. Combined with outflow samples in the literature, our work shows a tight outflow mass-size correlation.

Using the IceCube Neutrino Observatory, we search for high-energy neutrino emission coincident with compact binary mergers observed by the LIGO and Virgo gravitational wave (GW) detectors during their first and second observing runs. We present results from two searches targeting emission coincident with the sky localization of each gravitational wave event within a 1000 second time window centered around the reported merger time. One search uses a model-independent unbinned maximum likelihood analysis, which uses neutrino data from IceCube to search for point-like neutrino sources consistent with the sky localization of GW events. The other uses the Low-Latency Algorithm for Multi-messenger Astrophysics, which incorporates astrophysical priors through a Bayesian framework and includes LIGO-Virgo detector characteristics to determine the association between the GW source and the neutrinos. No significant neutrino coincidence is seen by either search during the first two observing runs of the LIGO-Virgo detectors. We set upper limits on the time-integrated neutrino emission within the 1000 second window for each of the 11 GW events. These limits range from 0.02-0.7 $\mathrm{GeV~cm^{-2}}$. We also set limits on the total isotropic equivalent energy, $E_{\mathrm{iso}}$, emitted in high-energy neutrinos by each GW event. These limits range from 1.7 $\times$ 10$^{51}$ - 1.8 $\times$ 10$^{55}$ erg. We conclude with an outlook for LIGO-Virgo observing run O3, during which both analyses are running in real time.

During the first three flights of the Antarctic Impulsive Transient Antenna (ANITA) experiment, the collaboration detected several neutrino candidates. Two of these candidate events were consistent with an ultra-high-energy up-going air shower and compatible with a tau neutrino interpretation. A third neutrino candidate event was detected in a search for Askaryan radiation in the Antarctic ice, although it is also consistent with the background expectation. The inferred emergence angle of the first two events is in tension with IceCube and ANITA limits on isotropic cosmogenic neutrino fluxes. Here, we test the hypothesis that these events are astrophysical in origin, possibly caused by a point source in the reconstructed direction. Given that any ultra-high-energy tau neutrino flux traversing the Earth should be accompanied by a secondary flux in the TeV-PeV range, we search for these secondary counterparts in seven years of IceCube data using three complementary approaches. In the absence of any significant detection, we set upper limits on the neutrino flux from potential point sources. We compare these limits to ANITA's sensitivity in the same direction and show that an astrophysical explanation of these anomalous events under standard model assumptions is severely constrained regardless of source spectrum.

Glitches correspond to sudden jumps of rotation frequency ($\nu$) and its derivative ($\dot{\nu}$) 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 ($\tau_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $\Delta{\nu}$ vs. $|\Delta{\dot\nu}|$ 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 ($|\Delta{\nu}_{\rm d1}|$) and recovery process ($\Delta{\nu}_{\rm d2}$), their time scales ($\tau_{1}$, $\tau_{2}$), and permanent changes in spin frequency ($\Delta{\nu}_{\rm p}$) and total frequency step ($\Delta{\nu}_{\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.

Glitches correspond to sudden jumps of rotation frequency ($\nu$) and its derivative ($\dot{\nu}$) 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 ($\tau_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $\Delta{\nu}$ vs. $|\Delta{\dot\nu}|$ 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 ($|\Delta{\nu}_{\rm d1}|$) and recovery process ($\Delta{\nu}_{\rm d2}$), their time scales ($\tau_{1}$, $\tau_{2}$), and permanent changes in spin frequency ($\Delta{\nu}_{\rm p}$) and total frequency step ($\Delta{\nu}_{\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.

Pulsar wind nebulae (PWNe) are main gamma-ray emitters in the Galactic plane. They are diffuse nebulae that emit non-thermal radiation. Pulsar winds, relativistic magnetized outflows from the central star, shocked in the ambient medium produce a multiwavelength emission from the radio through gamma rays. Although the leptonic scenario is able to explain most PWNe emission, a hadronic contribution cannot be excluded. A possible hadronic contribution to the high-energy gamma-ray emission inevitably leads to the production of neutrinos. Using 9.5 years of all-sky IceCube data, we report results from a stacking analysis to search for neutrino emission from 35 PWNe that are high-energy gamma-ray emitters. In the absence of any significant correlation, we set upper limits on the total neutrino emission from those PWNe and constraints on hadronic spectral components.

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

We present two searches for IceCube neutrino events coincident with 28 fast radio bursts (FRBs) and one repeating FRB. The first improves upon a previous IceCube analysis -- searching for spatial and temporal correlation of events with FRBs at energies greater than roughly 50 GeV -- by increasing the effective area by an order of magnitude. The second is a search for temporal correlation of MeV neutrino events with FRBs. No significant correlation is found in either search, therefore, we set upper limits on the time-integrated neutrino flux emitted by FRBs for a range of emission timescales less than one day. These are the first limits on FRB neutrino emission at the MeV scale, and the limits set at higher energies are an order-of-magnitude improvement over those set by any neutrino telescope.

The third-generation ground-based gravitational-wave (GW) detector, Cosmic Explorer (CE), is scheduled to start its observation in the 2030s. In this paper, we make a forecast for cosmological parameter estimation with gravitational-wave standard siren observation from the CE. We use the simulated GW standard siren data of CE to constrain the $\Lambda$CDM, $w$CDM and CPL models. We combine the simulated GW data with the current cosmological electromagnetic observations including the latest cosmic microwave background anisotropies data from Planck, the optical baryon acoustic oscillation measurements, and the type Ia supernovae observation (Pantheon compilation) to do the analysis. We find that the future standard siren observation from CE will improve the cosmological parameter estimation to a great extent, since the future GW standard siren data can well break the degeneracies generated by the optical observations between various cosmological parameters. We also find that the CE's constraining capability on the cosmological parameters is slightly better than that of the same-type GW detector, the Einstein Telescope. In addition, the synergy between the GW standard siren observation from CE and the 21 cm emission observation from SKA is also discussed.

We performed a multi-wavelength observation toward LkHa 101 embedded cluster and its adjacent 85arcmin*60arcmin region. The LkHa 101 embedded cluster is the first and only one significant cluster in California molecular cloud (CMC). These observations have revealed that the LkHa 101 embedded cluster is just located at the projected intersectional region of two filaments. One filament is the highest-density section of the CMC, the other is a new identified filament with a low-density gas emission. Toward the projected intersection, we find the bridging features connecting the two filaments in velocity, and identify a V-shape gas structure. These agree with the scenario that the two filaments are colliding with each other. Using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we measured that the RRL velocity of the LkHa 101 H II region is 0.5 km/s, which is related to the velocity component of the CMC filament. Moreover, there are some YSOs distributed outside the intersectional region. We suggest that the cloud-cloud collision together with the fragmentation of the main filament may play an important role in the YSOs formation of the cluster.

We present the calibration of the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) X-ray satellite, which can be used to perform timing and spectral studies of bright X-ray sources. Insight-HXMT carries three main payloads onboard: the High Energy X-ray telescope (HE), the Medium Energy X-ray telescope (ME) and the Low Energy X-ray telescope (LE). In orbit, the radioactive sources, activated lines, the fluorescence lines and celestial sources are used to calibrate the energy scale and energy resolution of the payloads. The Crab nebular is adopted as the primary effective area calibrator and empirical functions are constructed to modify the simulated effective areas of the three payloads respectively. The systematic errors of HE, compared to the model of the Crab nebular, are less than 2% in 28--120 keV and 2%--10% above 120 keV. The systematic errors of ME are less than 1.5% in 10--35 keV. The systematic errors of LE are less than 1% in 1--7 keV except the Si K--edge (1.839 keV, up to 1.5%) and less than 2% in 7--10 keV.

We present the results of the first combined dark matter search targeting the Galactic Centre using the ANTARES and IceCube neutrino telescopes. For dark matter particles with masses from 50 to 1000 GeV, the sensitivities on the self-annihilation cross section set by ANTARES and IceCube are comparable, making this mass range particularly interesting for a joint analysis. Dark matter self-annihilation through the $\tau^+\tau^-$, $\mu^+\mu^-$, $b\bar{b}$ and $W^+W^-$ channels is considered for both the Navarro-Frenk-White and Burkert halo profiles. In the combination of 2,101.6 days of ANTARES data and 1,007 days of IceCube data, no excess over the expected background is observed. Limits on the thermally-averaged dark matter annihilation cross section $\langle\sigma_A\upsilon\rangle$ are set. These limits present an improvement of up to a factor of two in the studied dark matter mass range with respect to the individual limits published by both collaborations. When considering dark matter particles with a mass of 200 GeV annihilating through the $\tau^+\tau^-$ channel, the value obtained for the limit is $7.44 \times 10^{-24} \text{cm}^{3}\text{s}^{-1}$ for the Navarro-Frenk-White halo profile. For the purpose of this joint analysis, the model parameters and the likelihood are unified, providing a benchmark for forthcoming dark matter searches performed by neutrino telescopes.

We present results from the Nuclear Spectroscopic Telescope Array (NuSTAR) observations of the new black hole X-ray binary candidate MAXI J1631-479 at two epochs during its 2018-2019 outburst, which caught the source in a disk dominant state and a power-law dominant state. Strong relativistic disk reflection features are clearly detected, displaying significant variations in the shape and strength of the broad iron emission line between the two states. Spectral modeling of the reflection spectra reveals that the inner radius of the optically-thick accretion disk evolves from $<1.9$ $r_{\rm g}$ to $12\pm1$ $r_{\rm g}$ (statistical errors at 90% confidence level) from the disk dominant to the power-law dominant state. Assuming in the former case that the inner disk radius is consistent with being at the ISCO, we estimate a black hole spin of $a^*>0.94$. Given that the bolometric luminosity is similar in the two states, our results indicate that the disk truncation observed in MAXI J1631-479 in the power-law dominant state is unlikely to be driven by a global variation in the accretion rate. We propose that it may instead arise from local instabilities in the inner edge of the accretion disk at high accretion rates. In addition, we find an absorption feature in the spectra centered at $7.33\pm0.03$ keV during the disk dominant state, which is evidence for a rare case that an extremely fast disk wind ($v_{\rm out}=0.067^{+0.001}_{-0.004}~c$) is observed in a low-inclination black hole binary, with the viewing angle of $29\pm1^{\circ}$ as determined by the reflection modeling.

The WorldWide Telescope(WWT) is a scientific visualization platform which can browse deep space images, star catalogs, and planetary remote sensing data from different observation facilities in a three-dimensional virtual scene. First launched and then open-sourced by Microsoft Research, the WWT is now managed by the American Astronomical Society (AAS). Hierarchical Progressive Survey (HiPS) is an astronomical data release scheme proposed by Centre de Donn\'ees astronomiques de Strasbourg (CDS) and has been accepted as a recommendation by International Virtual Observatory Alliance (IVOA). The HiPS solution has been adopted widely by many astronomical institutions for data release. Since WWT selected Hierarchical Triangular Mesh (HTM) as the standard for data visualization in the early stage of development, data released by HiPS cannot be visualized in WWT, which significantly limits the application of WWT. This paper introduces the implementation method for HiPS dataset visualization in WWT, and introduces HiPS data projection, mesh rendering, and data index implementation in WWT. Taking Chang'E-2 lunar probe data as an example, this paper introduces how to convert planetary remote sensing data into a HiPS dataset and integrate it into WWT. This paper also compares the efficiency and memory consumption of WWT loading its native data and HiPS data, and illustrates the application of HiPS in scientific data visualization and science education in WWT.

We present results from NuSTAR and XMM-Newton observations of the new black hole X-ray binary MAXI J1820+070 at low accretion rates (below 1% of the Eddington luminosity). We detect a narrow Fe K$\alpha$ emission line, in contrast to the broad and asymmetric Fe K$\alpha$ line profiles commonly present in black hole binaries at high accretion rates. The narrow line, with weak relativistic broadening, indicates that the Fe K$\alpha$ line is produced at a large disk radius. Fitting with disk reflection models assuming standard disk emissivity finds a large disk truncation radius (a few tens to a few hundreds of gravitational radii, depending on the disk inclination). In addition, we detect a quasi-periodic oscillation (QPO) varying in frequency between $11.6\pm0.2$~mHz and $2.8\pm0.1$~mHz. The very low QPO frequencies suggest a large size for the optically-thin Comptonization region according to the Lense-Thirring precession model, supporting that the accretion disk recedes from the ISCO and is replaced by advection-dominated accretion flow at low accretion rates. We also discuss the possibility of an alternative accretion geometry that the narrow Fe K$\alpha$ line is produced by a lamppost corona with a large height illuminating the disk.

The measurement of diffuse PeV gamma-ray emission from the Galactic plane would provide information about the energy spectrum and propagation of Galactic cosmic rays, and the detection of a point-like source of PeV gamma rays would be strong evidence for a Galactic source capable of accelerating cosmic rays up to at least a few PeV. This paper presents several un-binned maximum likelihood searches for PeV gamma rays in the Southern Hemisphere using 5 years of data from the IceTop air shower surface detector and the in-ice array of the IceCube Observatory. The combination of both detectors takes advantage of the low muon content and deep shower maximum of gamma-ray air showers, and provides excellent sensitivity to gamma rays between $\sim$0.6 PeV and 100 PeV. Our measurements of point-like and diffuse Galactic emission of PeV gamma rays are consistent with background, so we constrain the angle-integrated diffuse gamma-ray flux from the Galactic Plane at 2 PeV to $2.61 \times 10^{-19}$ cm$^{-2}$ s$^{-1}$ TeV$^{-1}$ at 90% confidence, assuming an E$^{-3}$ spectrum, and we estimate 90% upper limits on point-like emission at 2 PeV between 10$^{-21}$ - 10$^{-20}$ cm$^{-2}$ s$^{-1}$ TeV$^{-1}$ for an E$^{-2}$ spectrum, depending on declination. Furthermore, we exclude unbroken power-law emission up to 2 PeV for several TeV gamma-ray sources observed by H.E.S.S., and calculate upper limits on the energy cutoffs of these sources at 90% confidence. We also find no PeV gamma rays correlated with neutrinos from IceCube's high-energy starting event sample. These are currently the strongest constraints on PeV gamma-ray emission.

With its large field of view, GWAC can record hundreds of meteors every day. These meteors are valuable treasures for some meteor research groups. It is therefore very important to accurately find all of these meteors. To address the challenge of precisely distinguishing meteors from other elongated objects in a GWAC-like sky survey system, we design and implement a meteor candidate recognition algorithm, including the recognizing and morphology analysis of the light curves of the meteor candidates. Although the algorithm may filter out some real meteors, it can provide a sample of meteor with high confidence. After processing the images of Mini-GWAC taken in two months, we detect 109,000 elongated objects in which more than 90 percent of objects are not meteor. Among the elongated objects, about 5.9% objects are identified as meteors with high confidence, after the filters based upon an existence in a single frame, a single peak in the light curves, and a slow variation of the light curves.

The ground-based wide-angle camera array (GWAC) generates millions of single frame alerts per night. After the complicated and elaborate filters by multiple methods, a couple of dozens of candidates are still needed to be confirmed by follow-up observations in real-time. In order to free scientists from the complex and high-intensity follow-up tasks, we developed a Real-time Automatic transient Validation System (RAVS), and introduce here its system architecture, data processing flow, database schema, automatic follow-up control flow, and mobile message notification solution. This system is capable of automatically carrying out all operations in real-time without human intervention, including the validation of transient candidates, the adaptive light-curve sampling for identified targets in multi-band, and the pushing of observation results to the mobile client. The running of RAVS shows that an M-type stellar flare event can be well sampled by RAVS without a significant loss of the details, while the observing time is only less than one-third of the time coverage. Because the control logic of RAVS is designed to be independent of the telescope hardware, RAVS can be conveniently transplanted to other telescopes, especially the follow-up system of SVOM. Some future improvements are presented for the adaptive light-curve sampling, after taking into account both the brightness of sources and the evolution trends of the corresponding light-curves.

Many asteroids in the main belt are continuously pushed by Yarkovsky effect into regions of different mean motion resonances (MMRs) and then ejected out. They are considered as the principal source of near-Earth objects. We investigate in this paper the effects of the 7/3 MMR with Jupiter (J7/3 MMR) on the transportation of asteroids from Koronis and Eos families that reside respectively on the inner and outer side of the resonance. The fraction of asteroids that make successful crossing through the resonance and the escaping rate from the resonance are found to depend on the Yarkovsky drifting rate, the initial inclination and the migrating direction. The excitation of eccentricity and inclination due to the combined influence from both the resonance and Yarkovsky effect is discussed. Only the eccentricity can be pumped up considerably, and it is attributed mainly to the resonance. In the observational data, family members are also found in the resonance and on the opposite side of the resonance with respect to the corresponding family centre. The existence of these family members is explained using our results of numerical simulations. Finally, the replenishment of asteroids in the J7/3 MMR and the transportation of asteroids by it are discussed.

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

Based on a linearized model of the Yarkovsky effect, we investigate in this paper the dependence of the semimajor axis drift $\Delta a$ of a celestial body on its size, spinning obliquity, initial orbit and thermal parameters on its surface. With appropriate simplification and approximation, we obtain the analytical solutions to the perturbation equations for the motion of asteroids influenced by the Yarkovsky effect, and they are then verified by numerical simulations of the full equations of motion. These solutions present explicitly the dependencies of $\Delta a$ on the thermal and dynamical parameters of the asteroid. With these analytical formulae for $\Delta a$, we investigate the combined seasonal and diurnal Yarkovsky effects. The critical points where the migration direction reverses are calculated and the consequent selective effects according to the size and rotation state of asteroids are discussed. %Solely the Yarkovsky effect is found to be able to produce some ring structure in the aged circumstellar debris disk. Finally, we apply the analytical formulae to calculate the migration of Eos family members. The space distribution of asteroids is well reproduced. Our calculations suggest that statistically the orientations of spin axes of family members satisfy a random-obliquity distribution, and the rotation rate $\omega_{\rm rot}$ of asteroid depends on its size $R$ by $\omega_{\rm rot}\propto R^{-1}$.

The recent observation of interstellar objects, 1I/Oumuamua and 2I/Borisov cross the solar system opened new opportunities for planetary science and planetary defense. As the first confirmed objects originating outside of the solar system, there are myriads of origin questions to explore and discuss, including where they came from, how did they get here and what are they composed of. Besides, there is a need to be cognizant especially if such interstellar objects pass by the Earth of potential dangers of impact. Specifically, in the case of Oumuamua, which was detected after its perihelion, passed by the Earth at around 0.2 AU, with an estimated excess speed of 60 km/s relative to the Earth. Without enough forewarning time, a collision with such high-speed objects can pose a catastrophic danger to all life Earth. Such challenges underscore the importance of detection and exploration systems to study these interstellar visitors. The detection system can include a spacecraft constellation with zenith-pointing telescope spacecraft. After an event is detected, a spacecraft swarm can be deployed from Earth to flyby past the visitor. The flyby can then be designed to perform a proximity operation of interest. This work aims to develop algorithms to design these swarm missions through the IDEAS (Integrated Design Engineering & Automation of Swarms) architecture. Specifically, we develop automated algorithms to design an Earth-based detection constellation and a spacecraft swarm that generates detailed surface maps of the visitor during the rendezvous, along with their heliocentric cruise trajectories.

We report on the first measurement of the astrophysical neutrino flux using particle showers (cascades) in IceCube data from 2010 -- 2015. Assuming standard oscillations, the astrophysical neutrinos in this dedicated cascade sample are dominated ($\sim 90 \%$) by electron and tau flavors. The flux, observed in the energy range from $16\,\mathrm{TeV} $ to $2.6\,\mathrm{PeV}$, is consistent with a single power-law as expected from Fermi-type acceleration of high energy particles at astrophysical sources. We find the flux spectral index to be $\gamma=2.53\pm0.07$ and a flux normalization for each neutrino flavor of $\phi_{astro} = 1.66^{+0.25}_{-0.27}$ at $E_{0} = 100\, \mathrm{TeV}$. This flux of electron and tau neutrinos is in agreement with IceCube muon neutrino results and with all-neutrino flavor results. Results from fits assuming more complex neutrino flux models suggest a flux softening at high energies and a flux hardening at low energies (p-value $\ge 0.06$).

M82 hosts two well-known ultraluminous X-ray sources (ULXs). X-1, an intermediate-mass black hole (IMBH) candidate, and X-2, an ultraluminous X-ray pulsar (ULXP). Here we present a broadband X-ray spectral analysis of both sources based on ten observations made simultaneously with Chandra and NuSTAR. Chandra provides the high spatial resolution to resolve the crowded field in the 0.5--8 keV band, and NuSTAR provides the sensitive hard X-ray spectral data, extending the bandpass of our study above 10 keV. The observations, taken in the period 2015--2016, cover a period of flaring from X-1, allowing us to study the spectral evolution of this source with luminosity. During four of these observations, X-2 was found to be at a low flux level, allowing an unambiguous view of the emission from X-1. We find that the broadband X-ray emission from X-1 is consistent with that seen in other ULXs observed in detail with NuSTAR, with a spectrum that includes a broadened disk-like component and a high-energy tail. We find that the luminosity of the disk scales with inner disk temperature as L~T^-3/2 contrary to expectations of a standard accretion disk and previous results. These findings rule out a thermal state for sub-Eddington accretion and therefore do not support M82 X-1 as an IMBH candidate. We also find evidence that the neutral column density of the material in the line of sight increases with L$_X$, perhaps due to an increased mass outflow with accretion rate. For X-2, we do not find any significant spectral evolution, but we find the spectral parameters of the phase-averaged broadband emission are consistent with the pulsed emission at the highest X-ray luminosities.

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

A search for point-like and extended sources of cosmic neutrinos using data collected by the ANTARES and IceCube neutrino telescopes is presented. The data set consists of all the track-like and shower-like events pointing in the direction of the Southern Sky included in the nine-year ANTARES point-source analysis, combined with the through-going track-like events used in the seven-year IceCube point-source search. The advantageous field of view of ANTARES and the large size of IceCube are exploited to improve the sensitivity in the Southern Sky by a factor $\sim$2 compared to both individual analyses. In this work, the Southern Sky is scanned for possible excesses of spatial clustering, and the positions of preselected candidate sources are investigated. In addition, special focus is given to the region around the Galactic Centre, whereby a dedicated search at the location of SgrA* is performed, and to the location of the supernova remnant RXJ 1713.7-3946. No significant evidence for cosmic neutrino sources is found and upper limits on the flux from the various searches are presented.

The observation of the transient sky through a multitude of astrophysical messengers hasled to several scientific breakthroughs these last two decades thanks to the fast evolution ofthe observational techniques and strategies employed by the astronomers. Now, it requiresto be able to coordinate multi-wavelength and multi-messenger follow-up campaign withinstruments both in space and on ground jointly capable of scanning a large fraction of thesky with a high imaging cadency and duty cycle. In the optical domain, the key challengeof the wide field of view telescopes covering tens to hundreds of square degrees is to dealwith the detection, the identification and the classification of hundreds to thousands of opticaltransient (OT) candidates every night in a reasonable amount of time. In the last decade, newautomated tools based on machine learning approaches have been developed to perform thosetasks with a low computing time and a high classification efficiency. In this paper, we presentan efficient classification method using Convolutional Neural Networks (CNN) to discard anybogus falsely detected in astrophysical images in the optical domain. We designed this toolto improve the performances of the OT detection pipeline of the Ground Wide field AngleCameras (GWAC) telescopes, a network of robotic telescopes aiming at monitoring the opticaltransient sky down to R=16 with a 15 seconds imaging cadency. We applied our trainedCNN classifier on a sample of 1472 GWAC OT candidates detected by the real-time detectionpipeline. It yields a good classification performance with 94% of well classified event and afalse positive rate of 4%.

Swift J1858.6-0814 was discovered by Swift-BAT on October 25, 2018. Here we report on the first follow-up NuSTAR observation of the source, which shows variability spanning two orders of magnitude in count rate on timescales of ~10-100 s. The power-spectrum of the source does not show any quasi-periodic oscillations or periodicity, but has a large fractional rms amplitude of 147%$\pm3$%, exhibiting a number of large flares throughout the observation. The hardness ratio (defined as $R_{10-79 \rm keV}/R_{3-10 \rm keV}$) of the flares tends to be soft, while the source spans a range of hardness ratios during non-flaring periods. The X-ray spectrum of the source shows strong reflection features, which become more narrow and peaked during the non-flaring intervals. We fit an absorbed relativistic reflection model to the source spectra to place physical constraints on the system. Most notably, we find that the source exhibits a large and varying intrinsic absorbing column density ($N_{\rm H}=1.4-4.2\times10^{23}$ cm$^{-2}$). This large intrinsic absorption is further supported by the energy spectra extracted from two flares observed simultaneously by NuSTAR and NICER. We find that the inner accretion disk of the source has a low inclination, $i<29^{\circ}$ ( 3$\sigma$ upper-limit), while the iron abundance in the disk is close to solar, $A_{\rm Fe}=1.0\pm0.3$. We set a 90% confidence upper limit on the inner radius of the accretion disk of $r_{\rm in}<8 \ r_{\rm ISCO}$, and, by fixing $r_{\rm in}$ to be at $r_{\rm ISCO}$, a 90% confidence lower-limit on the spin of the black hole of $a^{*}>0.0$. Lastly, we compare the properties of Swift J1858.6$-$0814 to those of V404 Cygni and V4641 Sgr, which both show rapid flaring and a strong and variable absorption.

Cosmic-ray interactions with the solar atmosphere are expected to produce particle showers which in turn produce neutrinos from weak decays of mesons. These solar atmospheric neutrinos (SA$\nu$s) have never been observed experimentally. A detection would be an important step in understanding cosmic-ray propagation in the inner solar system and the dynamics of solar magnetic fields. SA$\nu$s also represent an irreducible background to solar dark matter searches and a detection would allow precise characterization of this background. Here, we present the first experimental search based on seven years of data collected from May 2010 to May 2017 in the austral winter with the IceCube Neutrino Observatory. An unbinned likelihood analysis is performed for events reconstructed within 5 degrees of the center of the Sun. No evidence for a SA$\nu$ flux is observed. After inclusion of systematic uncertainties, we set a 90\% upper limit of $1.02^{+0.20}_{-0.18}\cdot10^{-13}$~$\mathrm{GeV^{-1}cm^{-2}s^{-1}}$ at 1 TeV.

We study the 2018 outburst of Aql X-1 via the monitor of all sky X-ray image (MAXI) data. We show that the outburst starting in February 2018 is a member of short-low class in the frame of outburst duration and the peak count rate although the outburst morphology is slightly different from the other fast-rise-exponential-decay (FRED) type outbursts with a milder rising stage. We study the partial accretion in the weak propeller stage of Aql X-1 via the MAXI data of the 2018 outburst. We report on the spectral analysis of 3 observations of Aquila X-1 obtained by Insight - hard X-ray modulation telescope (Insight-HXMT) during the late decay stage of the 2018 outburst. We discuss that the data taken by Insight-HXMT is just after the transition to the weak propeller stage. Our analysis shows the necessity of a comptonization component to take into account the existence of an electron cloud resulting photons partly up-scattered.

In this paper we describe the first results of a compact imaging air-Cherenkov telescope, IceAct, operating in coincidence with the IceCube Neutrino Observatory (IceCube) at the geographic South Pole. An array of IceAct telescopes (referred to as the IceAct project) is under consideration as part of the IceCube-Gen2 extension to IceCube. Surface detectors in general will be a powerful tool in IceCube-Gen2 for distinguishing astrophysical neutrinos from the dominant backgrounds of cosmic-ray induced atmospheric muons and neutrinos: the IceTop array is already in place as part of IceCube, but has a high energy threshold. Although the duty cycle will be lower for the IceAct telescopes than the present IceTop tanks, the IceAct telescopes may prove to be more effective at lowering the detection threshold for air showers. Additionally, small imaging air-Cherenkov telescopes in combination with IceTop, the deep IceCube detector or other future detector systems might improve measurements of the composition of the cosmic ray energy spectrum. In this paper we present measurements of a first 7-pixel imaging air Cherenkov telescope demonstrator, proving the capability of this technology to measure air showers at the South Pole in coincidence with IceTop and the deep IceCube detector.

The current and upcoming generation of Very Large Volume Neutrino Telescopes---collecting unprecedented quantities of neutrino events---can be used to explore subtle effects in oscillation physics, such as (but not restricted to) the neutrino mass ordering. The sensitivity of an experiment to these effects can be estimated from Monte Carlo simulations. With the high number of events that will be collected, there is a trade-off between the computational expense of running such simulations and the inherent statistical uncertainty in the determined values. In such a scenario, it becomes impractical to produce and use adequately-sized sets of simulated events with traditional methods, such as Monte Carlo weighting. In this work we present a staged approach to the generation of binned event distributions in order to overcome these challenges. By combining multiple integration and smoothing techniques which address limited statistics from simulation it arrives at reliable analysis results using modest computational resources.

The distribution of galaxies within the local universe is characterized by anisotropic features. Observatories searching for the production sites of astrophysical neutrinos can take advantage of these features to establish directional correlations between a neutrino dataset and overdensities in the galaxy distribution in the sky. The results of two correlation searches between a seven-year time-integrated neutrino dataset from the IceCube Neutrino Observatory, and the 2MASS Redshift Survey (2MRS) catalog are presented here. The first analysis searches for neutrinos produced via interactions between diffuse intergalactic Ultra-High Energy Cosmic Rays (UHECRs) and the matter contained within galaxies. The second analysis searches for low-luminosity sources within the local universe, which would produce subthreshold multiplets in the IceCube dataset that directionally correlate with galaxy distribution. No significant correlations were observed in either analyses. Constraints are presented on the flux of neutrinos originating within the local universe through diffuse intergalactic UHECR interactions, as well as on the density of standard candle sources of neutrinos at low luminosities.

Many Galactic sources of gamma rays, such as supernova remnants, are expected to produce neutrinos with a typical energy cutoff well below 100 TeV. For the IceCube Neutrino Observatory located at the South Pole, the southern sky, containing the inner part of the Galactic plane and the Galactic Center, is a particularly challenging region at these energies, because of the large background of atmospheric muons. In this paper, we present recent advancements in data selection strategies for track-like muon neutrino events with energies below 100 TeV from the southern sky. The strategies utilize the outer detector regions as veto and features of the signal pattern to reduce the background of atmospheric muons to a level which, for the first time, allows IceCube searching for point-like sources of neutrinos in the southern sky at energies between 100 GeV and several TeV in the muon neutrino charged current channel. No significant clustering of neutrinos above background expectation was observed in four years of data recorded with the completed IceCube detector. Upper limits on the neutrino flux for a number of spectral hypotheses are reported for a list of astrophysical objects in the southern hemisphere.