X-ray binaries exhibit a wide range of properties but there are few accepted methods to determine the nature of the compact object. Color-Color-Intensity diagrams have been suggested as a means of distinguishing between systems containing black holes from those containing neutron stars. However, this technique has been verified with data from only one instrument (RXTE/ASM) with a single set of X-ray colors defined using data available only in pre-determined energy bands. We test a selection of X-ray colors with a more sensitive instrument to determine the reliability of this method. We use data from the MAXI Gas Slit Camera, which allows users to specify energy-bands. We test X-ray colors that have been previously defined in the literature as well as ones that we define specifically in this paper. A representative set of systems are used to construct Color-Color-Intensity diagrams in each set of colors to determine which are best for separating different classes. For studying individual sources certain bands are more effective than others. For a specified energy range, the separation of soft states in black hole binaries was possible only where both soft and hard colors included information from the lowest energy band. We confirm that Color-Color-Intensity diagrams can distinguish between systems containing black holes or neutron stars in all X-ray colors tested; this suggests an universality in the accretion processes governing these different classes. We suggest possible physical processes driving different classes of X-ray binaries to different locations in Color-Color-Intensity diagrams.

Our understanding of how AGN feedback operates in galaxy clusters has improved in recent years owing to large efforts in multi-wavelength observations and hydrodynamical simulations. However, it is much less clear how feedback operates in galaxy groups, which have shallower gravitational potentials. In this work, using very deep VLA and new MeerKAT observations from the MIGHTEE survey, we compiled a sample of 247 X-ray selected galaxy groups detected in the COSMOS field. We have studied the relation between the X-ray emission of the intra-group medium and the 1.4 GHz radio emission of the central radio galaxy. For comparison, we have also built a control sample of 142 galaxy clusters using ROSAT and NVSS data. We find that clusters and groups follow the same correlation between X-ray and radio emission. Large radio galaxies hosted in the centres of groups and merging clusters increase the scatter of the distribution. Using statistical tests and Monte-Carlo simulations, we show that the correlation is not dominated by biases or selection effects. We also find that galaxy groups are more likely than clusters to host large radio galaxies, perhaps owing to the lower ambient gas density or a more efficient accretion mode. In these groups, radiative cooling of the ICM could be less suppressed by AGN heating. We conclude that the feedback processes that operate in galaxy clusters are also effective in groups.

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.

We investigate the stellar and dynamical mass profiles in the centres of 25 brightest cluster galaxies (BCGs) at redshifts of 0.05 $\leq z \leq$ 0.30. Our spectroscopy enables us to robustly measure the Gauss-Hermite higher order velocity moments $h_{3}$ and $h_{4}$, which we compare to measurements for massive early-type galaxies, and central group galaxies. We measure positive central values for $h_{4}$ for all the BCGs. We derive the stellar mass-to-light ratio ($\Upsilon_{\star \rm DYN}$), and velocity anisotropy ($\beta$) based on a Multi-Gaussian Expansion (MGE) and axisymmetric Jeans Anisotropic Methods (JAM, cylindrically- and spherically-aligned). We explicitly include a dark matter halo mass component, which is constrained by weak gravitational lensing measurements for these clusters. We find a strong correlation between anisotropy and velocity dispersion profile slope, with rising velocity dispersion profiles corresponding to tangential anisotropy and decreasing velocity dispersion profiles corresponding to radial anisotropy. The rising velocity dispersion profiles can also indicate a significant contribution from the intracluster light (ICL) to the total light (in projection) in the centre of the galaxy. For a small number of BCGs with rising velocity dispersion profiles, a variable stellar mass-to-light ratio can also account for the profile shape, instead of tangential anisotropy or a significant ICL contribution. We note that, for some BCGs, a variable $\beta_{z}(r)$ (from radial to tangential anisotropy) can improve the model fit to the observed kinematic profiles. The observed diversity in these properties illustrates that BCGs are not the homogeneous class of objects they are often assumed to be.

We present results from GMRT and Chandra observations of the NGC 1550 galaxy group. Although previously thought of as relaxed, we show evidence that gas sloshing and active galactic nucleus (AGN) heating have affected the structure of the system. The 610 and 235 MHz radio images show an asymmetric jet-lobe structure with a total size of $\sim$33 kpc, with a sharp kink at the base of the more extended western jet, and bending of the shorter eastern jet as it enters the lobe. The 235$-$610 MHz spectral index map shows that both radio lobes have steep spectral indices ($\alpha_{235}^{610}\geq-1.5$) indicating the presence of an old electron population. The X-ray images reveal an asymmetric structure in the hot gas correlated with the radio structure, as well as potential cavities coincident with the radio lobes, with rims and arms of gas that may have been uplifted by the cavity expansion. The X-ray residual map reveals an arc shaped structure to the east that resembles a sloshing cold front. Radio spectral analysis suggests a radiative age of about 33 Myr for the source, comparable to the sloshing timescale and dynamical estimates of the age of the lobes. An estimate of the mechanical energy required to inflate the cavities suggests that the AGN of NGC 1550 is capable of balancing radiative losses from the intragroup medium (IGM) and preventing excessive cooling, providing that the AGN jets are efficiently coupled to the IGM gas. In conclusion, we find evidence of sloshing motions from both radio and X-ray structures, suggesting that NGC 1550 was perturbed by a minor merger or infalling galaxy about 33 Myr ago.

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.

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.

The fate of cooling gas in the centers of galaxy clusters and groups is still not well understood, as is also the case for the complex processes of triggering star formation in central dominant galaxies (CDGs), re-heating of cooled gas by AGN, and the triggering/feeding of supermassive black hole outbursts. We present CO observations of the early type galaxy NGC 5044, which resides at the center of an X-ray bright group with a moderate cooling flow. For our analysis we combine CO(2-1) data from the 7m antennae of the Atacama Compact Array (ACA), and the ACA total power array (TP). We demonstrate, using the 7m array data, that we can recover the total flux inferred from IRAM 30m single dish observations, which corresponds to a total molecular mass of about 4x10^7 Msun. Most of the recovered flux is blueshifted with respect to the galaxy rest frame and is extended on kpc-scales, suggesting low filling factor dispersed clouds. We find 8 concentrations of molecular gas out to a radius of 10 arcsec (1.5 kpc), which we identify with giant molecular clouds. The total molecular gas mass is more centrally concentrated than the X-ray emitting gas, but extended in the north-east/south-west direction beyond the IRAM 30m beam. We also compare the spatial extent of the molecular gas to the Halpha emission: The CO emission coincides with the very bright Halpha region in the center. We do not detect CO emission in the fainter Halpha regions. Furthermore, we find two CO absorption features spatially located at the center of the galaxy, within 5 pc projected distance of the AGN, infalling at 255 and 265 km/s relative to the AGN. This indicates that the two giant molecular clouds seen in absorption are most likely within the sphere of influence of the supermassive black hole.

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.

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.

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.

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.

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.

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

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.

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.

Using the data products of the Chandra Galaxy Atlas (Kim et al. 2019a), we have investigated the radial profiles of the hot gas temperature in 60 early type galaxies. Considering the characteristic temperature and radius of the peak, dip, and break (when scaled by the gas temperature and virial radius of each galaxy), we propose a universal temperature profile of the hot halo in ETGs. In this scheme, the hot gas temperature peaks at RMAX = 35 +/- 25 kpc (or ~0.04 RVIR) and declines both inward and outward. The temperature dips (or breaks) at RMIN (or RBREAK) = 3 - 5 kpc (or ~0.006 RVIR). The mean slope between RMIN (RBREAK) and RMAX is 0.3 +/- 0.1. Allowing for selection effects and observational limits, we find that the universal temperature profile can describe the temperature profiles of 72% (possibly up to 82%) of our ETG sample. The remaining ETGs (18%) with irregular or monotonically declining profiles do not fit the universal profile and require another explanation. The temperature gradient inside RMIN (RBREAK) varies widely, indicating different degrees of additional heating at small radii. Investigating the nature of the hot core (HC with a negative gradient inside RMIN), we find that HC is most clearly visible in small galaxies. Searching for potential clues associated with stellar, AGN feedback, and gravitational heating, we find that HC may be related to recent star formation. But we see no clear evidence that AGN feedback and gravitational heating play any significant role for HC.

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.

We investigate correlated gravitational wave and neutrino signals from rotating core-collapse supernovae with simulations. Using an improved mode identification procedure based on mode function matching, we show that a linear quadrupolar mode of the core produces a dual imprint on gravitational waves and neutrinos in the early post-bounce phase of the supernova. The angular harmonics of the neutrino emission are consistent with the mode energy around the neutrinospheres, which points to a mechanism for the imprint on neutrinos. Thus, neutrinos carry information about the mode amplitude in the outer region of the core, whereas gravitational waves probe deeper in. We also find that the best-fit mode function has a frequency bounded above by $\sim 420$ Hz, and yet the mode's frequency in our simulations is $\sim 15\%$ higher, due to the use of Newtonian hydrodynamics and a widely used pseudo-Newtonian gravity approximation. This overestimation is particularly important for the analysis of gravitational wave detectability and asteroseismology, pointing to limitations of pseudo-Newtonian approaches for these purposes, possibly even resulting in excitation of incorrect modes. In addition, mode frequency matching (as opposed to mode function matching) could be resulting in mode misidentification in recent work. Lastly, we evaluate the prospects of a multimessenger detection of the mode using current technology. The detection of the imprint on neutrinos is most challenging, with a maximum detection distance of $\sim\!1$ kpc using the IceCube Neutrino Observatory. The maximum distance for detecting the complementary gravitational wave imprint is $\sim\!5$ kpc using Advanced LIGO at design sensitivity.

The past decade has welcomed the emergence of cosmic neutrinos as a new messenger to explore the most extreme environments of the universe. The discovery measurement of cosmic neutrinos, announced by IceCube in 2013, has opened a new window of observation that has already resulted in new fundamental information that holds the potential to answer key questions associated with the high-energy universe, including: what are the sources in the PeV sky and how do they drive particle acceleration; where are cosmic rays of extreme energies produced, and on which paths do they propagate through the universe; and are there signatures of new physics at TeV-PeV energies and above? The planned advancements in neutrino telescope arrays in the next decade, in conjunction with continued progress in broad multimessenger astrophysics, promise to elevate the cosmic neutrino field from the discovery to the precision era and to a survey of the sources in the neutrino sky. The planned detector upgrades to the IceCube Neutrino Observatory, culminating in IceCube-Gen2 (an envisaged $400M facility with anticipated operation in the next decade, described in this white paper) are the cornerstone that will drive the evolution of neutrino astrophysics measurements.

This paper presents the results from point-like neutrino source searches using ten years of IceCube data collected between Apr.~6, 2008 and Jul.~10, 2018. We evaluate the significance of an astrophysical signal from a point-like source looking for an excess of clustered neutrino events with energies typically above $\sim1\,$TeV among the background of atmospheric muons and neutrinos. We perform a full-sky scan, a search within a selected source catalog, a catalog population study, and three stacked Galactic catalog searches. The most significant point in the Northern hemisphere from scanning the sky is coincident with the Seyfert II galaxy NGC 1068, which was included in the source catalog search. The excess at the coordinates of NGC 1068 is inconsistent with background expectations at the level of $2.9\,\sigma$ after accounting for statistical trials. The combination of this result along with excesses observed at the coordinates of three other sources, including TXS 0506+056, suggests that, collectively, correlations with sources in the Northern catalog are inconsistent with background at 3.3$\,\sigma$ significance. These results, all based on searches for a cumulative neutrino signal integrated over the ten years of available data, motivate further study of these and similar sources, including time-dependent analyses, multimessenger correlations, and the possibility of stronger evidence with coming upgrades to the detector.

We present a combined radio/X-ray analysis of the poorly studied galaxy cluster Abell 2495 (z=0.07923) based on new EVLA and Chandra data. We also analyze and discuss Halpha emission and optical continuum data retrieved from the literature. We find an offset of 6 kpc between the cluster BCG (MCG+02-58-021) and the peak of the X-ray emission, suggesting that the cooling process is not taking place on the central galaxy nucleus. We propose that sloshing of the ICM could be responsible for this separation. Furthermore, we detect a second, 4 kpc offset between the peak of the Halpha emission and that of the X-ray emission. Optical images highlight the presence of a dust filament extending up to 6 kpc in the cluster BCG, and allow us to estimate a dust mass within the central 7 kpc of 1.7e+5 Msun. Exploiting the dust to gas ratio and the L_Halpha-M_mol relation, we argue that a significant amount (up to 10^9 Msun) of molecular gas should be present in the BCG of this cluster. We also investigate the presence of ICM depressions, finding two putative systems of cavities; the inner pair is characterized by t age = 18 Myr and P cav = 1.2e+43 erg/s, the outer one by t age = 53 Myr and P cav = 5.6e+42 erg/s. Their age difference appears to be consistent with the free-fall time of the central cooling gas and with the offset timescale estimated with the Halpha kinematic data, suggesting that sloshing is likely playing a key role in this environment. Furthermore, the cavities' power analysis shows that the AGN energy injection is able to sustain the feedback cycle, despite cooling being offset from the BCG nucleus.

The presence of a population of point sources in a dataset modifies the underlying neutrino-count statistics from the Poisson distribution. This deviation can be exactly quantified using the non-Poissonian template fitting technique, and in this work we present the first application this approach to the IceCube high-energy neutrino dataset. Using this method, we search in 7 years of IceCube data for point-source populations correlated with the disk of the Milky Way, the Fermi bubbles, the Schlegel, Finkbeiner, and Davis dust map, or with the isotropic extragalactic sky. No evidence for such a population is found in the data using this technique, and in the absence of a signal we establish constraints on population models with source count distribution functions that can be described by a power-law with a single break. The derived limits can be interpreted in the context of many possible source classes. In order to enhance the flexibility of the results, we publish the full posterior from our analysis, which can be used to establish limits on specific population models that would contribute to the observed IceCube neutrino flux.

RNO is the mid-scale discovery instrument designed to make the first observation of neutrinos from the cosmos at extreme energies, with sensitivity well beyond current instrument capabilities. This new observatory will be the largest ground-based neutrino telescope to date, enabling the measurement of neutrinos above $10^{16}$ eV, determining the nature of the astrophysical neutrino flux that has been measured by IceCube at higher energies, similarly extending the reach of multi-messenger astrophysics to the highest energies, and enabling investigations of fundamental physics at energies unreachable by particle accelerators on Earth.

We present an implementation of IceCube in the SNOwGLoBES package, which is used to calculate expected detection event rates resulting from supernova neutrinos. The SNOwGLoBES package is widely used to compare the sensitivity of different neutrino observatories, but currently does not include simulation files for IceCube. In this paper, we give a brief overview of the design process that went into this implementation.

Efficient treatment of systematic uncertainties that depend on a large number of nuisance parameters is a persistent difficulty in particle physics experiments. Where low-level effects are not amenable to simple parameterization or re-weighting, analyses often rely on discrete simulation sets to quantify the effects of nuisance parameters on key analysis observables. Such methods may become computationally untenable for analyses requiring high statistics Monte Carlo with a large number of nuisance degrees of freedom, especially in cases where these degrees of freedom parameterize the shape of a continuous distribution. In this paper we present a method for treating systematic uncertainties in a computationally efficient and comprehensive manner using a single simulation set with multiple and continuously varied nuisance parameters. This method is demonstrated for the case of the depth-dependent effective dust distribution within the IceCube Neutrino Telescope.

The advent of Multi-Messenger Astronomy has allowed for new types of source searches within the neutrino community. We present the results of the first search for GeV astrophysical neutrinos emitted from Compact Binary Mergers, i.e. binary black hole or binary neutron star mergers, detected by the LIGO and Virgo interferometers. We introduce a new approach that lowers the energy threshold of IceCube from roughly 10 GeV to <1 GeV. This method uses an innovative event selection of GeV neutrino events in IceCube and searches for a statistically significant increase in the amount of GeV-like events detected around the Compact Binary Merger time. We compare our results with constraints set by high-energy neutrino searches, and describe the complementarity of these low and high-energy searches.

In 2017, IceCube detected a high energy neutrino in coincidence with the blazar TXS 0506+056. In a follow up analysis of archival data, evidence for previous neutrino emission from the blazar during 2014-15 was found. In this analysis, we investigate whether other blazars might have had episodes of neutrino emission. We use a 5-year sample of muon-neutrino data and the same timedependent algorithm that was used for the archival analysis of TXS 0506+056, applying it to all of the objects in the northern sky listed in the Fermi 3LAC catalog. We fit for the most significant neutrino flare from each direction, without assuming any correlation with electromagnetic light curves. To search for an excess of significant flares which could indicate a small population of neutrino emitters, we apply a binomial test. This binomial test does not find a significant excess of flares relative to the results expected for background coincidences from atmospheric neutrinos. The final significance of the analysis is p=0.244 when TXS 0506+056 is excluded, consistent with the background-only hypothesis. We report the parameters of the individual fits for the most significant objects which contributed to the final result.