We report a search for astronomical neutrinos in the energy region from several GeV to TeV in the direction of the blazar TXS0506+056 using the Super-Kamiokande detector following the detection of a 100 TeV neutrino from the same location by the IceCube collaboration. Using Super-Kamiokande neutrino data across several data samples observed from April 1996 to February 2018 we have searched for both a total excess above known backgrounds across the entire period as well as localized excesses on smaller time scales in that interval. No significant excess nor significant variation in the observed event rate are found in the blazar direction. Upper limits are placed on the electron and muon neutrino fluxes at 90\% confidence level as $6.03 \times 10^{-7}$ and $4.52 \times 10^{-7}$ to $9.26 \times 10^{-10}$ [${\rm erg}/{\rm cm}^2/{\rm s}$], respectively.

Supernova detection is a major objective of the Super-Kamiokande (SK) experiment. In the next stage of SK (SK-Gd), gadolinium (Gd) sulfate will be added to the detector, which will improve the ability of the detector to identify neutrons. A core-collapse supernova will be preceded by an increasing flux of neutrinos and anti-neutrinos, from thermal and weak nuclear processes in the star, over a timescale of hours; some of which may be detected at SK-Gd. This could provide an early warning of an imminent core-collapse supernova, hours earlier than the detection of the neutrinos from core collapse. Electron anti-neutrino detection will rely on inverse beta decay events below the usual analysis energy threshold of SK, so Gd loading is vital to reduce backgrounds while maximising detection efficiency. Assuming normal neutrino mass ordering, more than 200 events could be detected in the final 12 hours before core collapse for a 15-25 solar mass star at around 200 pc, which is representative of the nearest red supergiant to Earth, $\mathrm{\alpha}$Ori (Betelgeuse). At a statistical false alarm rate of 1 per century, detection could be up to 10 hours before core collapse, and a pre-supernova star could be detected by SK-Gd up to 600 pc away. A pre-supernova alert could be provided to the astrophysics community following gadolinium loading.

Cassini/UVIS FUV observations of stellar occultations at Titan are well suited for probing its atmospheric composition and structure. However, due to instrument pointing motion, only five out of tens of observations have been analyzed. We present an innovative retrieval method that corrects for the effect of pointing motion by forward modeling the Cassini/UVIS instrument response function with the pointing motion value obtained from the SPICE C-kernel along the spectral dimension. To illustrate the methodology, an occultation observation made during flyby T52 is analyzed, when the Cassini spacecraft had insufficient attitude control. A high-resolution stellar model and an instrument response simulator that includes the position of the point source on the detector are used for the analysis of the pointing motion. The Markov Chain Monte-Carlo method is used to retrieve the line-of-sight abundance profiles of eleven species (CH4, C2H2, C2H4, C2H6, C4H2, C6H6, HCN, C2N2, HC3N, C6N2 and haze particles) in the spectral vector fitting process. We obtain tight constraints on all of the species aside from C2H6, C2N2 and C6N2, for which we only retrieved upper limits. This is the first time that the T52 occultation was used to derive abundances of major hydrocarbon and nitrile species in Titan's upper and middle atmosphere, as pointing motion prohibited prior analysis. With this new method, nearly all of the occultations obtained over the entire Cassini mission could yield reliable profiles of atmospheric composition, allowing exploration of Titan's upper atmosphere over seasons, latitudes, and longitudes.

Magnetic reconnection is a fundamental physical process in various astrophysical, space, and laboratory environments. Many pieces of evidence for magnetic reconnection have been uncovered. However, its specific processes that could be fragmented and turbulent have been short of direct observational evidence. Here, we present observations of a super-hot current sheet during SOL2017-09-10T X8.2-class solar flare that display the fragmented and turbulent nature of magnetic reconnection. As bilateral plasmas converge toward the current sheet, significant plasma heating and non-thermal motions are detected therein. Two oppositely directed outflow jets are intermittently expelled out of the fragmenting current sheet, whose intensity shows a power-law distribution in spatial frequency domain. The intensity and velocity of the sunward outflow jets also display a power-law distribution in temporal frequency domain. The length-to-width ratio of the current sheet is estimated to be larger than theoretical threshold of and thus ensures occurrence of tearing mode instability. The observations therefore suggest fragmented and turbulent magnetic reconnection occurring in the long stretching current sheet.

We report the results of a neutrino search in Super-Kamiokande for coincident signals with the first detected gravitational wave produced by a binary neutron star merger, GW170817, which was followed by a short gamma-ray burst, GRB170817A, and a kilonova/macronova. We searched for coincident neutrino events in the range from 3.5 MeV to $\sim$100 PeV, in a time window $\pm$500 seconds around the gravitational wave detection time, as well as during a 14-day period after the detection. No significant neutrino signal was observed for either time window. We calculated 90% confidence level upper limits on the neutrino fluence for GW170817. From the upward-going-muon events in the energy region above 1.6 GeV, the neutrino fluence limit is $16.0^{+0.7}_{-0.6}$ ($21.3^{+1.1}_{-0.8}$) cm$^{-2}$ for muon neutrinos (muon antineutrinos), with an error range of $\pm5^{\circ}$ around the zenith angle of NGC4993, and the energy spectrum is under the assumption of an index of $-2$. The fluence limit for neutrino energies less than 100 MeV, for which the emission mechanism would be different than for higher-energy neutrinos, is also calculated. It is $6.6 \times 10^7$ cm$^{-2}$ for anti-electron neutrinos under the assumption of a Fermi-Dirac spectrum with average energy of 20 MeV.

We present the results of a search in the Super-Kamiokande (SK) detector for excesses of neutrinos with energies above a few GeV that are in the direction of the track events reported in IceCube. Data from all SK phases (SK-I through SK-IV) were used, spanning a period from April 1996 to April 2016 and corresponding to an exposure of 225 kilotonne-years . We considered the 14 IceCube track events from a data set with 1347 livetime days taken from 2010 to 2014. We use Poisson counting to determine if there is an excess of neutrinos detected in SK in a 10 degree search cone (5 degrees for the highest energy data set) around the reconstructed direction of the IceCube event. No significant excess was found in any of the search directions we examined. We also looked for coincidences with a recently reported IceCube multiplet event. No events were detected within a $\pm$ 500 s time window around the first detected event, and no significant excess was seen from that direction over the lifetime of SK.

We report the results from a search in Super-Kamiokande for neutrino signals coincident with the first detected gravitational wave events, GW150914 and GW151226, using a neutrino energy range from 3.5 MeV to 100 PeV. We searched for coincident neutrino events within a time window of $\pm$500 seconds around the gravitational wave detection time. Four neutrino candidates are found for GW150914 and no candidates are found for GW151226. The remaining neutrino candidates are consistent with the expected background events. We calculated the 90\% confidence level upper limits on the combined neutrino fluence for both gravitational wave events, which depends on event energy and topologies. Considering the upward going muon data set (1.6 GeV - 100 PeV) the neutrino fluence limit for each gravitational wave event is 14 - 37 (19 - 50) cm$^{-2}$ for muon neutrinos (muon antineutrinos), depending on the zenith angle of the event. In the other data sets, the combined fluence limits for both gravitational wave events range from 2.4$\times 10^{4}$ to 7.0$\times 10^{9}$ cm$^{-2}$.

In this paper, we study the formation and early evolution of a limb coronal mass ejection (CME) and its associated shock wave that occurred on 2014 January 8. The extreme ultraviolet (EUV) images provided by the Atmospheric Imaging Assembly (AIA) on board \textit{Solar Dynamics Observatory} disclose that the CME first appears as a bubble-like structure. Subsequently, its expansion forms the CME and causes a quasi-circular EUV wave. Interestingly, both the CME and the wave front are clearly visible at all of the AIA EUV passbands. Through a detailed kinematical analysis, it is found that the expansion of the CME undergoes two phases: a first phase with a strong but transient lateral over-expansion followed by a second phase with a self-similar expansion. The temporal evolution of the expansion velocity coincides very well with the variation of the 25--50 keV hard X-ray flux of the associated flare, which indicates that magnetic reconnection most likely plays an important role in driving the expansion. Moreover, we find that, when the velocity of the CME reaches $\sim$600 km s$^{-1}$, the EUV wave starts to evolve into a shock wave, which is evidenced by the appearance of a type II radio burst. The shock's formation height is estimated to be $\sim$0.2$R_{sun}$, which is much lower than the height derived previously. Finally, we also study the thermal properties of the CME and the EUV wave. We find that the plasma in the CME leading front and the wave front has a temperature of $\sim$2 MK, while that in the CME core region and the flare region has a much higher temperature of $\ge$8 MK.

The relation between star formation rates and stellar masses, i.e. the galaxy main sequence, is a useful diagnostic of galaxy evolution. We present the distributions relative to the main sequence of 55 optically-selected PG and 12 near-IR-selected 2MASS quasars at z <= 0.5. We estimate the quasar host stellar masses from Hubble Space Telescope or ground-based AO photometry, and the star formation rates through the mid-infrared aromatic features and far-IR photometry. We find that PG quasar hosts more or less follow the main sequence defined by normal star-forming galaxies while 2MASS quasar hosts lie systematically above the main sequence. PG and 2MASS quasars with higher nuclear luminosities seem to have higher specific SFRs (sSFRs), although there is a large scatter. No trends are seen between sSFRs and SMBH masses, Eddington ratios or even morphology types (ellipticals, spirals and mergers). Our results could be placed in an evolutionary scenario with quasars emerging during the transition from ULIRGs/mergers to ellipticals. However, combined with results at higher redshift, they suggest that quasars can be widely triggered in normal galaxies as long as they contain abundant gas and have ongoing star formation.

Accurately predicting the arrival of coronal mass ejections (CMEs) at the Earth based on remote images is of critical significance in the study of space weather. In this paper, we make a statistical study of 21 Earth directed CMEs, exploring in particular the relationship between CME initial speeds and transit times. The initial speed of a CME is obtained by fitting the CME with the Graduated Cylindrical Shell model and is thus free of projection effects. We then use the drag force model to fit results of the transit time versus the initial speed. By adopting different drag regimes, i.e., the viscous, aerodynamics, and hybrid regimes, we get similar results, with the least mean estimation error of the hybrid model of 12.9 hours. CMEs with a propagation angle (the angle between the propagation direction and the Sun-Earth line) larger than its half angular width arrive at the Earth with an angular deviation caused by factors other than the radial solar wind drag. The drag force model cannot be well applied to such events. If we exclude these events in the sample, the prediction accuracy can be improved, i.e., the estimation error reduces to 6.8 hours. This work suggests that it is viable to predict the arrival time of CMEs at the Earth based on the initial parameters with a fairly good accuracy. Thus, it provides a method of space weather forecast of 1--5 days following the occurrence of CMEs.

Gamma ray burst (GRB) optical and X-ray afterglow luminosity is expected to correlate with the GRB isotropic equivalent kinetic energy of the outflow in the standard synchrotron model for GRB afterglows. Previous studies, using prompt GRB isotropic equivalent energy ($E_{\rm iso}$) as a proxy for isotropic equivalent kinetic energy, have generally confirmed a correlation between X-ray and optical afterglow luminosities. Assuming that GRB afterglow luminosity does not evolve strongly with redshift, we identify a strong Malmquist bias in GRB optical and X-ray afterglow luminosity data. We show that selection effects dominate the observed E$_{\rm iso}$ -- L$_{\rm opt,X}$ correlations, and have likely been underestimated in other studies. The bias is strongest for a subset of optically faint bursts $m>24$ at 24 hr with $z>2$. After removing this optical selection bias, the E$_{\rm iso}$ -- L$_{\rm opt,X}$ correlation for long GRBs is not statistically significant, but combining both long and short GRB luminosity data the correlation is significant. Using the median of the $E_{\rm iso}$ and $L_{\rm opt,X}$ distributions, we apply the synchrotron model assuming the same power law index for short and long GRBs, but different microphysical parameter distributions. Comparing the ratio of optical and X-ray luminosities, we find tentative evidence that the fraction of post-shock energy in magnetic fields, $\epsilon_B$, could be systematically higher in SGRBs compared to LGRBs.

Multi-frequency radio observations of the radio bridge of a powerful classical double radio source can be used to determine: the beam power of the jets emanating from the AGN; the total time the source will actively produce jets that power large-scale radio emission; the thermal pressure of the medium in the vicinity of the radio source; and the total mass, including dark matter, of the galaxy or cluster of galaxies traced by the ambient gas that surrounds the radio source. Empirical determinations of each of these quantities are obtained and analyzed for 22 radio sources. Typical beam powers are about $10^{45} {erg s}^{-1}$. The characteristic or total time the AGN will actively produce a collimated outflow is estimated. Typical total lifetimes are $\sim (10^7$ to $10^8$) years. Total masses, and mass-density profiles, similar to those of low-redshift clusters of galaxies are obtained. Thus, some clusters of galaxies, or cores of clusters, exist at redshifts of one to two. A new method of estimating the thermal pressure of the ambient gas in the vicinity of a powerful classical double radio source is presented. This new estimate is independent of synchrotron and inverse Compton aging arguments, and depends only upon the properties of the radio lobe and the shape of the radio bridge. Thermal pressures on the order of $10^{-10} {dyne cm}^{-2}$, typical of gas in low-redshift clusters of galaxies, are found for the environments of the sources studied here. It is shown that appreciable amounts of cosmic microwave background diminution (the Sunyaev-Zel'dovich effect) are expected from many of these clusters. This could be detected at high frequency where the emission from the radio sources is weak.

A sample of 20 powerful extended radio galaxies with redshifts between zero and two were used to determine constraints on global cosmological parameters. Data for six radio sources were obtained from the VLA archive, analyzed, and combined with the sample of 14 radio galaxies used previously by Guerra & Daly to determine cosmological parameters. The results are consistent with our previous results, and indicate that the current value of the mean mass density of the universe is significantly less than the critical value. A universe with $\Omega_m$ of unity is ruled out at 99.0% confidence, and the best fitting values of $\Omega_m$ in matter are $0.10^{+0.25}_{-0.10}$ and $-0.25^{+0.35}_{-0.25}$ assuming zero space curvature and zero cosmological constant, respectively. Note that identical results obtain when the low redshift bin, which includes Cygnus A, is excluded; these results are independent of whether the radio source Cygnus A is included. The method does not rely on a zero-redshift normalization. The radio properties of each source are also used to determine the density of the gas in the vicinity of the source, and the beam power of the source. The six new radio sources have physical characteristics similar to those found for the original 14 sources. The density of the gas around these radio sources is typical of gas in present day clusters of galaxies. The beam powers are typically about $10^{45} \hbox{erg s}^{-1}$.

A sample of 20 powerful extended radio galaxies with redshifts between zero and two were used to determine constraints on global cosmological parameters. Data for six radio sources were obtained from the VLA archive, analyzed, and combined with the sample of 14 radio galaxies used previously by Guerra & Daly to determine cosmological parameters. The new results are consistent with our previous results, and indicate that the current value of the mean mass density of the universe is significantly less than the critical value. A universe with $\Omega_m$ in matter of unity is ruled out at 99.0% confidence, and the best fitting values of $\Omega_m$ in matter are $0.10^{+0.25}_{-0.10}$ and $-0.25^{+0.35}_{-0.25}$ assuming zero space curvature and zero cosmological constant, respectively. The radio properties of each source are also used to determine the density of the gas in the vicinity of the source, and the beam power of the source. The six new radio sources have physical characteristics similar to those found for the original 14 sources. The density of the gas around these radio sources is typical of gas in present day clusters of galaxies. The beam powers are typically about $10^{45} {erg s}^{-1}$.

Multi-frequency radio observations of the radio bridge of powerful classical double radio sources can be used to determine: the beam power of the jets emanating from the AGN; the total time the source will actively produce jets that power large-scale radio emission; the thermal pressure of the medium in the vicinity of the radio source; and the total mass, including dark matter, of the galaxy or cluster of galaxies traced by the ambient gas that surrounds the radio source. The theoretical constructs that allow a determination of each of these quantities using radio observations are presented and discussed. Empirical determinations of each of these quantities are obtained and analyzed. A sample of 14 radio galaxies and 8 radio loud quasars with redshifts between zero and two for which there is enough radio information to be able to determine the physical parameters listed above was studied in detail. (abridged)

Fourteen classical double radio galaxies with redshifts between zero and two were used to determine the cosmological parameters $\Omega_m$, $\Omega_{\Lambda}$, and $\Omega_k$, where these are the normalized values of the mean mass density, cosmological constant, and space curvature at the present epoch. A low value of $\Omega_m$ is obtained, and $\Omega_m = 1$ is ruled out with 97.5% confidence. The low value of $\Omega_m$ determined using the radio source method described here is also indicated by several independent tests. Thus, it appears that either a cosmological constant, or space curvature, is significant at the present epoch. This means that the universe is undergoing, or has recently undergone, a transition away from a state of matter domination and into a state where either a cosmological constant or space curvature is determining the expansion rate of the universe. The low value of $\Omega_m$ presented here and by Guerra & Daly (1998) means that we can state with 97.5 % confidence that the universe will continue to expand forever.

A sample of FRII radio sources in rich clusters of galaxies, both at high (z$\sim$0.5) and low (z$\sim$0) redshift, has been constructed to study the effect of environment on radio sources. Comparisons are made between the properties of FRII sources in cluster and non-cluster environments, and between X-ray clusters with and without FRII sources. The principal results are the following: 1. Most low-redshift FRII sources in clusters appear to be similar to FRII sources in group or field environments in terms of radio power, optical properties of the host galaxy, and nonthermal pressure of the radio bridge. Most low-redshift clusters with FRII sources tend to lie at the lower end of cluster X-ray luminosity distribution, some having $L_x$ comparable to FRII sources in non-cluster environments. 2. High-redshift FRII sources are all quite similar to each other irrespective of their environments. The nonthermal pressures of their radio bridges appear to be similar to those of low-redshift FRII's. 3. The nonthermal pressures of the bridges of FRII sources appear to be similar to the thermal pressures of the ICM around them. This result would allow FRII sources to be used as probes of their gaseous environments. The current data imply that the evolution in the clustering strength around FRII sources toward high-redshift is likely to be closely linked to an evolution of the state of the intracluster medium.