Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

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.

The potential association between the blazar TXS 0506+056 and the neutrino event IceCube-170922A provides a unique opportunity to study the possible physical connection between the high-energy photons and neutrinos. We explore the correlated electromagnetic and neutrino emissions of blazar TXS 0506+056 by a self-consistent leptonic-hadronic model, taking into account particle stochastic acceleration and all relevant radiative processes self-consistently. The electromagnetic and neutrino spectra of blazar TXS 0506+056 are reproduced by the proton synchrotron and hybrid leptonic-hadronic models based on the proton-photon interactions. It is found that the hybrid leptonic-hadronic model can be used to better explain the observed X-ray and $\gamma$-ray spectra of blazar TXS 0506+056 than the proton synchrotron model. Moreover, the predicted neutrino spectrum of the hybrid leptonic-hadronic model is closer to the observed one compared to the proton synchrotron model. We suggest that the hybrid leptonic-hadronic model is more favored if the neutrino event IceCube-170922A is associated with the blazar TXS 0506+056.

A good compromise between the resistive model and the PIC model is Aristotelian electrodynamics, which can include the back-reaction of the radiative photons onto particle motion and allow for a local dissipation where the force-free condition is violated. We study the dissipative pulsar magnetosphere with Aristotelian electrodynamics where particle acceleration is fully balanced by radiation. The expression for the current density is defined by introducing a pair multiplicity. The 3D structure of the pulsar magnetosphere is then presented by solving the time-dependent Maxwell equations using a pseudo-spectral algorithm. It is found that the dissipative magnetosphere approaches the force-free solution and the dissipative region is more restricted to the current sheet outside the light-cylinder (LC) as the pair multiplicity increases. The spatial extension of the dissipative region is self-consistently controlled by the pair multiplicity. Our simulations show the high magnetospheric dissipation outside the LC for the low pair multiplicity.

The primary advantage of moderately superheated bubble chamber detectors is their simultaneous sensitivity to nuclear recoils from WIMP dark matter and insensitivity to electron recoil backgrounds. A comprehensive analysis of PICO gamma calibration data demonstrates for the first time that electron recoils in C$_3$F$_8$ scale in accordance with a new nucleation mechanism, rather than one driven by a hot-spike as previously supposed. Using this semi-empirical model, bubble chamber nucleation thresholds may be tuned to be sensitive to lower energy nuclear recoils while maintaining excellent electron recoil rejection. The PICO-40L detector will exploit this model to achieve thermodynamic thresholds as low as 2.8 keV while being dominated by single-scatter events from coherent elastic neutrino-nucleus scattering of solar neutrinos. In one year of operation, PICO-40L can improve existing leading limits from PICO on spin-dependent WIMP-proton coupling by nearly an order of magnitude for WIMP masses greater than 3 GeV c$^{-2}$ and will have the ability to surpass all existing non-xenon bounds on spin-independent WIMP-nucleon coupling for WIMP masses from 3 to 40 GeV c$^{-2}$.

In this paper we report on the characterization of the Hamamatsu VUV4 (S/N: S13370-6152) Vacuum Ultra-Violet (VUV) sensitive Silicon Photo-Multipliers (SiPMs) as part of the development of a solution for the detection of liquid xenon scintillation light for the nEXO experiment. Various SiPM features, such as: dark noise, gain, correlated avalanches, direct crosstalk and Photon Detection Efficiency (PDE) were measured in a dedicated setup at TRIUMF. SiPMs were characterized in the range $163 \text{ } \text{K} \leq \text{T}\leq 233 \text{ } \text{K}$. At an over voltage of $3.1\pm0.2$ V and at $\text{T}=163 \text{ }\text{K}$ we report a number of Correlated Avalanches (CAs) per pulse in the $1 \upmu\text{s}$ interval following the trigger pulse of $0.161\pm0.005$. At the same settings the Dark-Noise (DN) rate is $0.137\pm0.002 \text{ Hz/mm}^{2}$. Both the number of CAs and the DN rate are within nEXO specifications. The PDE of the Hamamatsu VUV4 was measured for two different devices at $\text{T}=233 \text{ }\text{K}$ for a mean wavelength of $189\pm7\text{ nm}$. At $3.6\pm0.2$ V and $3.5\pm0.2$ V of over voltage we report a PDE of $13.4\pm2.6\text{ }\%$ and $11\pm2\%$, corresponding to a saturation PDE of $14.8\pm2.8\text{ }\%$ and $12.2\pm2.3\%$, respectively. Both values are well below the $24\text{ }\%$ saturation PDE advertised by Hamamatsu. More generally, the second device tested at $3.5\pm0.2$ V of over voltage is below the nEXO PDE requirement. The first one instead yields a PDE that is marginally close to meeting the nEXO specifications. This suggests that with modest improvements the Hamamatsu VUV4 MPPCs could be considered as an alternative to the FBK-LF SiPMs for the final design of the nEXO detector.

We study the gamma-ray emission patterns and light curves in dissipative pulsar magnetospheres. We produce the gamma-ray light curves by using the geometric method and the particle trajectory method. For the geometric method, assuming the gamma-ray emission originates in a finite-width layer along the last closed lines, we generate the gamma-ray light curves based on the uniform emissivity along the magnetic field lines in the comoving frame (CF). For the particle trajectory method, we consider the spatial distribution of conductivity $\sigma$ by assuming a very high conductivity within the light cylinder and a finite conductivity outside the light cylinder (LC). Assuming that all the $\gamma$-ray emission originates in the outer magnetosphere outside the LC, we generate the gamma-ray light curves by computing realistic particle trajectories and Lorentz factors, taking into account both the accelerating electric field and curvature radiation loss. Further, we compare the modeling light curves to the observed light curves at $>0.1\, \rm GeV$ energies for Vela pulsar. Our results show that the magnetosphere with the low $\sigma$ value is more preferred for the geometric method. However, the magnetosphere with a near force-free regime within the LC and a high $\sigma$ value outside the LC is more favored for the particle trajectory method. It is noted that the particle trajectory method uses the parallel electric fields that are self-consistent with the magnetic fields of the magnetosphere. We suggest that the results from the particle trajectory method are more supported on the physical ground.

Final results are reported from operation of the PICO-60 C$_3$F$_8$ dark matter detector, a bubble chamber filled with 52 kg of C$_3$F$_8$ located in the SNOLAB underground laboratory. The chamber was operated at thermodynamic thresholds as low as 1.2 keV without loss of stability. A new blind 1404-kg-day exposure at 2.45 keV threshold was acquired with approximately the same expected total background rate as the previous 1167-kg-day exposure at 3.3 keV. This increased exposure is enabled in part by a new optical tracking analysis to better identify events near detector walls, permitting a larger fiducial volume. These results set the most stringent direct-detection constraint to date on the WIMP-proton spin-dependent cross section at 2.5 $\times$ 10$^{-41}$ cm$^2$ for a 25 GeV WIMP, and improve on previous PICO results for 3-5 GeV WIMPs by an order of magnitude.

In this work, we propose a new source for gravitational wave (GW) radiation associated with the quantum chromodynamics (QCD) phase transition in the inner cores of neutron stars. The mechanism is based on the bubble dynamics during the first-order phase transition from nuclear matter to quark matter. We identify the characteristic frequency to be of order $\omega_c\sim 10^6~{\rm rad/s}$ for this kind of sources and the strain magnitude ($h\sim 10^{-24}$ for a neutron star at a distance of $0.1~{\rm Mpc}$) reachable by future GW detectors. The GW spectra are shown to be useful to check the transition nature at high baryon chemical potential as well as to constrain the radius and density of the inner cores, which are still indistinct up to now.

New results are reported from the operation of the PICO-60 dark matter detector, a bubble chamber filled with 52 kg of C$_3$F$_8$ located in the SNOLAB underground laboratory. As in previous PICO bubble chambers, PICO-60 C$_3$F$_8$ exhibits excellent electron recoil and alpha decay rejection, and the observed multiple-scattering neutron rate indicates a single-scatter neutron background of less than 1 event per month. A blind analysis of an efficiency-corrected 1167-kg-day exposure at a 3.3-keV thermodynamic threshold reveals no single-scattering nuclear recoil candidates, consistent with the predicted background. These results set the most stringent direct-detection constraint to date on the WIMP-proton spin-dependent cross section at 3.4 $\times$ 10$^{-41}$ cm$^2$ for a 30-GeV$\thinspace$c$^{-2}$ WIMP, more than one order of magnitude improvement from previous PICO results.

An oblique pulsar magnetosphere with a plasma conductivity is studied by using a pseudo-spectral method. In the pseudo-spectral method, the time-dependent Maxwell equations are solved, both electric and magnetic fields are expanded in terms of the vector spherical harmonic (VSH) functions in spherical geometry and the divergencelessness of magnetic field is analytically enforced by a projection method. The pulsar magnetospheres in infinite (i. e., force-free approximation) and finite conductivities are simulated and a family of solutions that smoothly transition from the Deutsch vacuum solution to the force-free solution are obtained. The $\sin^2\alpha$ dependence of the spin-down luminosity on the magnetic inclination angle $\alpha$ in which the full electric current density are taken into account is retrieved in the force-free regime.

A pseudo-spectral method with an absorbing outer boundary is used to solve a set of the time-dependent force-free equations. In the method, both electric and magnetic fields are expanded in terms of the vector spherical harmonic (VSH) functions in spherical geometry and the divergencelessness of magnetic field is analytically enforced by a projection method. Our simulations show that the Deutsch vacuum solution and the Michel monopole solution can be well reproduced by our pseudo-spectral code. Further the method is used to present the time-dependent simulation of the force-free pulsar magnetosphere for an aligned rotator. The simulations show that the current sheet in the equatorial plane can be resolved well, and the obtained spin-down luminosity in the steady state is in good agreement with the value given by Spitkovsky (2006).

The very hard $\gamma$-ray spectrum from distant blazars challenges the traditional synchrotron self-Compton (SSC) model, which may indicate that there is the contribution of an additional high-energy component beyond the SSC emission. In this paper, we study the possible origin of the hard $\gamma$-ray spectrum from distant blazars. We develop a model to explain the hard $\gamma$-ray spectrum from blazar 1ES 1101-232. In the model, the optical and X-ray radiation would come from the synchrotron radiation of primary electrons and secondary pairs, the GeV emission would be produced by the SSC process, however, the hard $\gamma$-ray spectrum would originate from the decay of neutral pion produced through proton-photon interactions with the synchrotron radiation photons within the jet. Our model can explain the observed SED of 1ES 1101-232 well, especially the very hard $\gamma$-ray spectrum. However, our model requires the very large proton power to efficiently produce the $\gamma$-ray through proton-photon interactions.

We try to infer the location of the GeV emission region for 21 flat spectrum radio quasars (FSRQs) with quasi-simultaneous spectral energy distributions (SEDs), in which the SEDs of 21 FSRQs are reproduced by the one-zone leptonic model including the synchrotron-self Compton (SSC) and external Compton (EC) processes. We suggest that the X-ray emission could be produced by the SSC process and the GeV emission could come from the EC process. The EC emission could originate from the inverse Compton (IC) scattering of photons from the broad line region (BLR) and accretion disk or dust torus by the same electron population, which mainly depend on the location of the $\gamma$-ray emission region. We propose a method to constrain the location of the GeV emission region based on the spectral shapes. When the GeV emission is located within the BLR, the IC scattering could occur at the Klein-Nishina regime and produce a broken/steep spectrum in the GeV energy band. When the GeV emission is produced outside the BLR, the IC scattering could take place at the Thomson regime and the GeV spectrum would have the same spectral index as the optical-infrared spectrum. We infer that the location of the GeV emission region is inside the BLR for 5 FSRQs and beyond the BLR for 16 FSRQs. Our results show that the ratio of the magnetic field and electron energy density is close to equipartition condition for 21 FSRQs.

We investigate the X-ray and $\gamma$-ray flares of Mrk 421 on 2008 June 6-15 using the synchrotron self- Compton(SSC) model with electron acceleration, in which an evident correlation between the X-ray and $\gamma$-ray bands appears, while no significant correlation between the optical and X-ray band is observed. We argue that the emission from Mrk 421 may originate from two different components. One is the steady component from the outer region that is mainly attributed to the optical band, in which the electrons are accelerated by first-order Fermi acceleration mechanism. We use a steady electron spectrum to produce the synchrotron self-Compton emission. The other is the variable component from the inner region, in which the electrons are accelerated by the stochastic acceleration process. We use the time-dependent SSC model to produce the emission from the variable component. We suggest that the flares are due to the hardening of the electron spectrum under the process of the stochastic acceleration, which leads to the hardening of the observed spectrum in the X-ray and $\gamma$-ray bands. Furthermore, we find that the energy densities of electrons and magnetic fields are near equipartition in both jet regions.