We present the discoveries of WISEA J041451.67-585456.7 and WISEA J181006.18-101000.5, two low-temperature (1200$-$1400 K), high proper motion T-type subdwarfs. Both objects were discovered via their high proper motion ($>$0.5 arcsec yr$^{-1}$); WISEA J181006.18-101000.5 as part of the NEOWISE proper motion survey and WISEA J041451.67-585456.7 as part of the citizen science project Backyard Worlds; Planet 9. We have confirmed both as brown dwarfs with follow-up near-infrared spectroscopy. Their spectra and near-infrared colors are unique amongst known brown dwarfs, with some colors consistent with L-type brown dwarfs and other colors resembling those of the latest-type T dwarfs. While no forward model consistently reproduces the features seen in their near-infrared spectra, the closest matches suggest very low metallicities ([Fe/H] $\leq$ -1), making these objects likely the first examples of extreme subdwarfs of the T spectral class (esdT). WISEA J041451.67-585456.7 and WISEA J181006.18-101000.5 are found to be part of a small population of objects that occupy the "substellar transition zone," and have the lowest masses and effective temperatures of all objects in this group.

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 report the results of a ${\sim}4$-year direct imaging survey of 104 stars to resolve and characterize circumstellar debris disks in scattered light as part of the Gemini Planet Imager Exoplanet Survey. We targeted nearby (${\lesssim}150$ pc), young (${\lesssim}500$ Myr) stars with high infrared excesses ($L_{\mathrm{IR}} / L_\star > 10^{-5}$), including 38 with previously resolved disks. Observations were made using the Gemini Planet Imager high-contrast integral field spectrograph in $H$-band (1.6 $\mu$m) coronagraphic polarimetry mode to measure both polarized and total intensities. We resolved 26 debris disks and three protoplanetary/transitional disks. Seven debris disks were resolved in scattered light for the first time, including newly presented HD 117214 and HD 156623, and we quantified basic morphologies of five of them using radiative transfer models. All of our detected debris disks but HD 156623 have dust-poor inner holes, and their scattered-light radii are generally larger than corresponding radii measured from resolved thermal emission and those inferred from spectral energy distributions. To assess sensitivity, we report contrasts and consider causes of non-detections. Detections were strongly correlated with high IR excess and high inclination, although polarimetry outperformed total intensity angular differential imaging for detecting low inclination disks (${\lesssim} 70 \deg$). Based on post-survey statistics, we improved upon our pre-survey target prioritization metric predicting polarimetric disk detectability. We also examined scattered-light disks in the contexts of gas, far-IR, and millimeter detections. Comparing $H$-band and ALMA fluxes for two disks revealed tentative evidence for differing grain properties. Finally, we found no preference for debris disks to be detected in scattered light if wide-separation substellar companions were present.

HR4796A hosts a well-studied debris disk with a long history due to its high fractional luminosity and favorable inclination lending itself well to both unresolved and resolved observations. We present new J- and K1-band images of the resolved debris disk HR4796A taken in the polarimetric mode of the Gemini Planet Imager (GPI). The polarized intensity features a strongly forward scattered brightness distribution and is undetected at the far side of the disk. The total intensity is detected at all scattering angles and also exhibits a strong forward scattering peak. We use a forward modelled geometric disk in order to extract geometric parameters, polarized fraction and total intensity scattering phase functions for these data as well as H-band data previously taken by GPI. We find the polarized phase function becomes increasingly more forward scattering as wavelength increases. We fit Mie and distribution of hollow spheres grain (DHS) models to the extracted functions. We find that while it is possible to describe generate a satisfactory model for the total intensity using a DHS model, but not with a Mie model. We find that no single grain population of DHS or Mie grains of arbitrary composition can simultaneously reproduce the polarized fraction and total intensity scattering phase functions, indicating the need for more sophisticated grain models.

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.

Despite their astrophysical relevance, nuclear pasta phases are relatively unstudied at high temperatures. We present molecular dynamics simulations of symmetric nuclear matter with several topologies of `lasagna' at a range of temperatures to study the pasta-uniform transition. Using the Minkowski functionals we quantify trends in the occupied volume, surface area, mean breadth, and Euler characteristic. The amplitude of surface displacements of the pasta increase with temperature which produce short lived topological defects such as holes and filaments near melting, resulting in power laws for increasing surface curvature with temperature. We calculate the static structure factor and report the shear viscosity and thermal conductivity of pasta, finding that the shear viscosity is minimized at the melting temperature. These results may have implications for the thermoelastic properties of nuclear pasta and finite temperature corrections to the equation of state at pasta densities.

Quantifying the evolution of stellar extreme ultraviolet (EUV, 100 -- 1000 $\overset{\circ}{A}$) emission is critical for assessing the evolution of planetary atmospheres and the habitability of M dwarf systems. Previous studies from the HAbitable Zones and M dwarf Activity across Time (HAZMAT) program showed the far- and near-UV (FUV, NUV) emission from M stars at various stages of a stellar lifetime through photometric measurements from the Galaxy Evolution Explorer (GALEX). The results revealed increased levels of short-wavelength emission that remain elevated for hundreds of millions of years. The trend for EUV flux as a function of age could not be determined empirically because absorption by the interstellar medium prevents access to the EUV wavelengths for the vast majority of stars. In this paper, we model the evolution of EUV flux from early M stars to address this observational gap. We present synthetic spectra spanning EUV to infrared wavelengths of 0.4 $\pm$ 0.05 M$_{\odot}$ stars at five distinct ages between 10 and 5000 Myr, computed with the PHOENIX atmosphere code and guided by the GALEX photometry. We model a range of EUV fluxes spanning two orders of magnitude, consistent with the observed spread in X-ray, FUV, and NUV flux at each epoch. Our results show that the stellar EUV emission from young M stars is 100 times stronger than field age M stars, and decreases as t$^{-1}$ after remaining constant for a few hundred million years. This decline stems from changes in the chromospheric temperature structure, which steadily shifts outward with time. Our models reconstruct the full spectrally and temporally resolved history of an M star's UV radiation, including the unobservable EUV radiation, which drives planetary atmospheric escape, directly impacting a planet's potential for habitability.

Cold brown dwarfs are excellent analogs of widely separated, gas giant exoplanets, and provide insight into the potential atmospheric chemistry and physics we may encounter in objects discovered by future direct imaging surveys. We present a low resolution R $\sim$ 300 $M$-band spectroscopic sequence of seven brown dwarfs with effective temperatures between 750 K and 250 K along with Jupiter. These spectra reveal disequilibrium abundances of carbon monoxide (CO) produced by atmospheric quenching. We use the eddy diffusion coefficient (K$_{zz}$) to estimate the strength of vertical mixing in each object. The K$_{zz}$ values of cooler gaseous objects are close to their theoretical maximum and warmer objects show weaker mixing, likely due to less efficient convective mixing in primarily radiative layers. The CO-derived K$_{zz}$ values imply that disequilibrium phosphine (PH$_{3}$) should be easily observable in all of the brown dwarfs, but none as yet show any evidence for PH$_{3}$ absorption. We find that ammonia is relatively insensitive to atmospheric quenching at these effective temperatures. We are able to improve the fit to WISE 0855's $M$-band spectrum by including both CO and water clouds in the atmospheric model.

We present the discovery of WISEA J083011.95+283716.0, the first Y dwarf candidate identified through the Backyard Worlds: Planet 9 citizen science project. We identified this object as a red, fast-moving source with a faint $W2$ detection in multi-epoch \textit{AllWISE} and unWISE images. We have characterized this object with Spitzer Space Telescope and \textit{Hubble Space Telescope} follow-up imaging. With mid-infrared detections in \textit{Spitzer}'s \emph{ch1} and \emph{ch2} bands and flux upper limits in Hubble Space Telescope $F105W$ and $F125W$ filters, we find that this object is both very faint and has extremely red colors ($ch1-ch2 = 3.25\pm0.23$ mag, $F125W-ch2 \geq 9.36$ mag), consistent with a T$_{eff}\sim300$ K source, as estimated from the known Y dwarf population. A preliminary parallax provides a distance of $11.1^{+2.0}_{-1.5}$ pc, leading to a slightly warmer temperature of $\sim350$ K. The extreme faintness and red Hubble Space Telescope and Spitzer Space Telescope colors of this object suggest it may be a link between the broader Y dwarf population and the coldest known brown dwarf WISE J0855$-$0714, and highlight our limited knowledge of the true spread of Y dwarf colors. We also present four additional Backyard Worlds: Planet 9 late-T brown dwarf discoveries within 30 pc.

We present new $H$-band scattered light images of the HD 32297 edge-on debris disk obtained with the Gemini Planet Imager (GPI). The disk is detected in total and polarized intensity down to a projected angular separation of 0.15", or 20au. On the other hand, the large scale swept-back halo remains undetected, likely a consequence of its markedly blue color relative to the parent body belt. We analyze the curvature of the disk spine and estimate a radius of $\approx$100au for the parent body belt, smaller than past scattered light studies but consistent with thermal emission maps of the system. We employ three different flux-preserving post-processing methods to suppress the residual starlight and evaluate the surface brightness and polarization profile along the disk spine. Unlike past studies of the system, our high fidelity images reveal the disk to be highly symmetric and devoid of morphological and surface brightness perturbations. We find the dust scattering properties of the system to be consistent with those observed in other debris disks, with the exception of HR 4796. Finally, we find no direct evidence for the presence of a planetary-mass object in the system.

We present a revision to the astrometric calibration of the Gemini Planet Imager (GPI), an instrument designed to achieve the high contrast at small angular separations necessary to image substellar and planetary-mass companions around nearby, young stars. We identified several issues with the GPI Data Reduction Pipeline (DRP) that significantly affected the determination of angle of north in reduced GPI images. As well as introducing a small error in position angle measurements for targets observed at small zenith distances, this error led to a significant error in the previous astrometric calibration that has affected all subsequent astrometric measurements. We present a detailed description of these issues, and how they were corrected. We reduced GPI observations of calibration binaries taken periodically since the instrument was commissioned in 2014 using an updated version of the DRP. These measurements were compared to observations obtained with the NIRC2 instrument on Keck II, an instrument with an excellent astrometric calibration, allowing us to derive an updated plate scale and north offset angle for GPI. This revised astrometric calibration should be used to calibrate all measurements obtained with GPI for the purposes of precision astrometry.

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.

We present a study of the HD 165054 astrometric calibration field that has been periodically observed with the Gemini Planet Imager. HD 165054 is a bright star within Baade's Window, a region of the galactic plane with relatively low extinction from interstellar dust. HD 165054 was selected as a calibrator target due to the high number density of stars within this region ($\sim 3$ stars per square arcsecond with $H<22$), necessary because of the small field-of-view of the Gemini Planet Imager. Using nine epochs spanning over five years, we have fit a standard five-parameter astrometric model to the astrometry of seven background stars within close proximity to HD 165054 (angular separation $< 2$ arcsec). We achieved a proper motion precision of $\sim 0.3$ mas/yr, and constrained the parallax of each star to be $\lesssim 1$ mas. Our measured proper motions and parallax limits are consistent with the background stars being a part of the galactic bulge. Using these measurements we find no evidence of any systematic trend of either the plate scale or the north angle offset of GPI between 2014 and 2019. We compared our model describing the motions of the seven background stars to observations of the same field in 2014 and 2018 obtained with Keck/NIRC2, an instrument with an excellent astrometric calibration. We find that predicted position of the background sources is consistent with that measured by NIRC2, within the uncertainties of the calibration of the two instruments. In the future, we will use this field as a standard astrometric calibrator for the upgrade of GPI and potentially for other high-contrast imagers.

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.

An accurate modelling of baryonic feedback effects is required to exploit the full potential of future weak-lensing surveys such as Euclid or LSST. In this second paper in a series of two, we combine Euclid-like mock data of the cosmic shear power spectrum with an eROSITA X-ray mock of the cluster gas fraction to run a combined likelihood analysis including both cosmological and baryonic parameters. Following the first paper of this series, the baryonic effects (based on the baryonic correction model of Schneider et al. 2019) are included in both the tomographic power spectrum and the covariance matrix. However, this time we assume the more realistic case of a $\Lambda$CDM cosmology with massive neutrinos, and we consider several extensions of the currently favoured cosmological model. For the standard $\Lambda$CDM case, we show that including X-ray data reduces the uncertainties on the sum of the neutrino mass by $\sim30$ percent, while there is only a mild improvement on other parameters such as $\Omega_m$ and $\sigma_8$. As extensions of $\Lambda$CDM, we consider the cases of a dynamical dark energy model (wCDM), a $f(R)$ gravity model (fRCDM), and a mixed dark matter model ($\Lambda$MDM) with both a cold and a warm/hot dark matter component. We find that combining weak lensing with X-ray data only leads to a mild improvement of the constraints on the additional parameters of wCDM, while the improvement is more substantial for both fRCDM and $\Lambda$MDM. Ignoring baryonic effects in the analysis pipeline leads to significant false-detections of either phantom dark energy or a light subdominant dark matter component. Overall we conclude that for all cosmologies considered, a general parametrisation of baryonic effects is both necessary and sufficient to obtain tight constraints on cosmological parameters.

Baryonic feedback effects lead to a suppression of the weak lensing angular power spectrum on small scales. The poorly constrained shape and amplitude of this suppression is an important source of uncertainties for upcoming cosmological weak lensing surveys such as Euclid or LSST. In this first paper in a series of two, we use simulations to build a Euclid-like tomographic mock data-set for the cosmic shear power spectrum and the corresponding covariance matrix, which are both corrected for baryonic effects following the baryonification method of Schneider et al. (2019). In addition, we develop an emulator to obtain fast predictions of the baryonic power suppression, allowing us to perform a likelihood inference analysis for a standard $\Lambda$CDM cosmology with both cosmological and astrophysical parameters. Our main findings are the following: (i) ignoring baryonic effects leads to a greater than 5$\sigma$ bias on the cosmological parameters $\Omega_m$ and $\sigma_8$; (ii) restricting the analysis to the largest scales, that are mostly unaffected by baryons, makes the bias disappear, but results in a blow-up of the $\Omega_m$-$\sigma_8$ contour area by more than a factor of 10; (iii) ignoring baryonic effects on the covariance matrix does not significantly affect cosmological parameter estimates; (iv) while the baryonic suppression is mildly cosmology dependent, this effect does not noticeably modify the posterior contours. Overall, we conclude that including baryonic uncertainties in terms of nuisance parameters results in unbiased and surprisingly tight constraints on cosmology.

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.

The core collapse of a massive star results in the formation of a proto-neutron star (PNS). If enough material is accreted onto a PNS it will become gravitationally unstable and further collapse into a black-hole (BH). We perform a systematic study of failing core-collapse supernovae in spherical symmetry for a wide range of presupernova progenitor stars and equations of state (EOSs) of nuclear matter. We analyze how variations in progenitor structure and the EOS of dense matter above nuclear saturation density affect the PNS evolution and subsequent BH formation. Comparisons of core-collapse for a given progenitor star and different EOSs show that the path traced by the PNS in mass-entropy phase space $M_{\mathrm{grav}}^{\mathrm{PNS}}-\tilde{s}$ is well correlated with the progenitor compactness and almost EOS independent, apart from the final endpoint. Furthermore, BH formation occurs, to a very good approximation, soon after the PNS overcomes the maximum \textit{gravitational} mass supported by a hot NS with constant entropy equal to $\tilde{s}$. These results show a path to constraining the temperature dependence of the EOS through the detection of neutrinos from a failed galactic supernova.

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.

We search for evidence of the cause of the exoplanet radius gap, i.e. the dearth of planets with radii near $1.8\ R_\oplus$. If the cause was photoevaporation, the radius gap should trend with proxies for the early-life high-energy emission of planet-hosting stars. If, alternatively, the cause was core-powered mass loss, no such trends should exist. Critically, spurious trends between the radius gap and stellar properties arise from an underlying correlation with instellation. After accounting for this underlying correlation, we find no trends remain between the radius gap and stellar mass or present-day stellar activity as measured by near-UV emission. We dismiss the nondetection of a radius gap trend with near-UV emission because present-day near-UV emission is unlikely to trace early-life high-energy emission, but we provide a catalog of GALEX near-UV and far-UV emission measurements for general use. We interpret the nondetection of a radius gap trend with stellar mass by simulating photoevaporation with mass-dependent evolution of stellar high-energy emission. The simulation produces an undetectable trend between the radius gap and stellar mass under realistic sources of error. We conclude that no evidence, from this analysis or others in the literature, currently exists that clearly favors either photoevaporation or core powered mass loss as the primary cause of the exoplanet radius gap. However, repeating this analysis once the body of well-characterized $< 4\ R_\oplus$ planets has roughly doubled could confirm or rule out photoevaporation.

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.

We present the first polarimetric detection of the inner disk component around the pre-main sequence B9.5 star HD 141569A. Gemini Planet Imager H-band (1.65 micron) polarimetric differential imaging reveals the highest signal-to-noise ratio detection of this ring yet attained and traces structure inwards to 0.25" (28 AU at a distance of 111 pc). The radial polarized intensity image shows the east side of the disk, peaking in intensity at 0.40" (44 AU) and extending out to 0.9" (100 AU). There is a spiral arm-like enhancement to the south, reminiscent of the known spiral structures on the outer rings of the disk. The location of the spiral arm is coincident with 12CO J=3-2 emission detected by ALMA, and hints at a dynamically active inner circumstellar region. Our observations also show a portion of the middle dusty ring at ~220 AU known from previous observations of this system. We fit the polarized H-band emission with a continuum radiative transfer Mie model. Our best-fit model favors an optically thin disk with a minimum dust grain size close to the blow-out size for this system: evidence of on-going dust production in the inner reaches of the disk. The thermal emission from this model accounts for virtually all of the far-infrared and millimeter flux from the entire HD 141569A disk, in agreement with the lack of ALMA continuum and CO emission beyond ~100 AU. A remaining 8-30 micron thermal excess a factor of ~2 above our model argues for a yet-unresolved warm innermost 5-15 AU component of the disk.

We present new observations of the planet beta Pictoris b from 2018 with GPI, the first GPI observations following conjunction. Based on these new measurements, we perform a joint orbit fit to the available relative astrometry from ground-based imaging, the Hipparcos Intermediate Astrometric Data (IAD), and the Gaia DR2 position, and demonstrate how to incorporate the IAD into direct imaging orbit fits. We find a mass consistent with predictions of hot-start evolutionary models and previous works following similar methods, though with larger uncertainties: 12.8 [+5.3, -3.2] M_Jup. Our eccentricity determination of 0.12 [+0.04, -0.03] disfavors circular orbits. We consider orbit fits to several different imaging datasets, and find generally similar posteriors on the mass for each combination of imaging data. Our analysis underscores the importance of performing joint fits to the absolute and relative astrometry simultaneously, given the strong covariance between orbital elements. Time of conjunction is well constrained within 2.8 days of 2017 September 13, with the star behind the planet's Hill sphere between 2017 April 11 and 2018 February 16 (+/- 18 days). Following the recent radial velocity detection of a second planet in the system, beta Pic c, we perform additional two-planet fits combining relative astrometry, absolute astrometry, and stellar radial velocities. These joint fits find a significantly smaller mass for the imaged planet beta Pic b, of 8.0 +/- 2.6 M_Jup, in a somewhat more circular orbit. We expect future ground-based observations to further constrain the visual orbit and mass of the planet in advance of the release of Gaia DR4.

We present the first spatially resolved scattered-light images of four debris disks around members of the Scorpius-Centaurus (Sco-Cen) OB Association with high-contrast imaging and polarimetry using the Gemini Planet Imager (GPI). All four disks are resolved for the first time in polarized light and one disk is also detected in total intensity. The three disks imaged around HD 111161, HD 143675, and HD 145560 are symmetric in both morphology and brightness distribution. The three systems span a range of inclinations and radial extents. The disk imaged around HD 98363 shows indications of asymmetries in morphology and brightness distribution, with some structural similarities to the HD 106906 planet-disk system. Uniquely, HD 98363 has a wide co-moving stellar companion Wray 15-788 with a recently resolved disk with very different morphological properties. HD 98363 A/B is the first binary debris disk system with two spatially resolved disks. All four targets have been observed with ALMA, and their continuum fluxes range from one non-detection to one of the brightest disks in the region. With the new results, a total of 15 A/F-stars in Sco-Cen have resolved scattered light debris disks, and approximately half of these systems exhibit some form of asymmetry. Combining the GPI disk structure results with information from the literature on millimeter fluxes and imaged planets reveals a diversity of disk properties in this young population. Overall, the four newly resolved disks contribute to the census of disk structures measured around A/F-stars at this important stage in the development of planetary systems.

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.

This white paper proposes leveraging high-quality Gaia data available to the worldwide scientific community and complement it with support from Canadian-related facilities to place Canada as a leader in the fields of stellar associations and exoplanet science, and to train Canadian highly qualified personnel through graduate and post-graduate research grants. Gaia has sparked a new era in the study of stellar kinematics by measuring precise distances and proper motions for 1.3 billion stars. These data have already generated more than 1700 scientific papers and are guaranteed to remain the source of many more papers for the upcoming decades. More than 900 new age-calibrated young low-mass stars have already been discovered as a direct consequence of the second Gaia data release. Some of these may already be host stars to known exoplanet systems or may become so with the progress of the TESS mission that is expected to discover 10,000 nearby transiting exoplanets in the upcoming decade. This places Canada in a strategic position to leverage Gaia data because it has access to several high-resolution spectrometers on 1-4 m class telescopes (e.g. The ESPaDOnS, SPIRou and NIRPS), that would allow to quickly characterize this large number of low-mass stars and their exoplanet systems. This white paper describes the opportunity in such scientific projects that could place Canada as a leader in the fields of stellar associations and exoplanets.

We report the discovery of WISE2150-7520AB (W2150AB): a widely separated (~ 341 AU) very low mass L1 + T8 co-moving system. The system consists of the previously known L1 primary 2MASS J21501592-7520367 and a newly discovered T8 secondary found at position 21:50:18.99 -75:20:54.6 (MJD=57947) using Wide-field Infrared Survey Explorer (WISE) data via the Backyard Worlds: Planet 9 citizen science project. We present Spitzer ch1 and ch2 photometry (ch1-ch2= 1.41 +/-0.04 mag) of the secondary and FIRE prism spectra of both components. The sources show no peculiar spectral or photometric signatures indicating that each component is likely field age. Using all observed data and the Gaia DR2 parallax of 41.3593 +/- 0.2799 mas for W2150A we deduce fundamental parameters of log(Lbol/Lsun)=-3.69 +/- 0.01, Teff=2118 +/- 62 K, and an estimated mass=72 +/- 12 MJup for the L1 and log(Lbol/Lsun)=-5.64 +/- 0.02, Teff=719 +/- 61 K, and an estimated mass=34 +/- 22 MJup for the T8. At a physical separation of ~341 AU this system has Ebin = 10^41 erg making it the lowest binding energy system of any pair with Mtot < 0.1 Msun not associated with a young cluster. It is equivalent in estimated mass ratio, Ebin, and physical separation to the ~ 2 Myr M7.25 + M8.25 binary brown dwarf 2MASS J11011926-7732383AB (2M1101AB) found in the Chameleon star forming region. W2150AB is the widest companion system yet observed in the field where the primary is an L dwarf or later.

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.

We present a revision to the visual orbit of the young, directly-imaged exoplanet 51 Eridani b using four years of observations with the Gemini Planet Imager. The relative astrometry is consistent with an eccentric ($e=0.53_{-0.13}^{+0.09}$) orbit at an intermediate inclination ($i=136_{-11}^{+10}$\,deg), although circular orbits cannot be excluded due to the complex shape of the multidimensional posterior distribution. We find a semi-major axis of $11.1_{-1.3}^{+4.2}$\,au and a period of $28.1_{-4.9}^{+17.2}$\,yr, assuming a mass of 1.75\,M$_{\odot}$ for the host star. We find consistent values with a recent analysis of VLT/SPHERE data covering a similar baseline. We investigated the potential of using absolute astrometry of the host star to obtain a dynamical mass constraint for the planet. The astrometric acceleration of 51~Eri derived from a comparison of the {\it Hipparcos} and {\it Gaia} catalogues was found to be inconsistent at the 2--3$\sigma$ level with the predicted reflex motion induced by the orbiting planet. Potential sources of this inconsistency include a combination of random and systematic errors between the two astrometric catalogs or the signature of an additional companion within the system interior to current detection limits. We also explored the potential of using {\it Gaia} astrometry alone for a dynamical mass measurement of the planet by simulating {\it Gaia} measurements of the motion of the photocenter of the system over the course of the extended eight-year mission. We find that such a measurement is only possible ($>98$\% probability) given the most optimistic predictions for the {\it Gaia} scan astrometric uncertainties for bright stars, and a high mass for the planet ($\gtrsim3.6$\,M$_{\rm Jup}$).

The $\sim500$\, Myr A2IV star HR 1645 has one of the most significant low-amplitude accelerations of nearby early-type stars measured from a comparison of the {\it Hipparcos} and {\it Gaia} astrometric catalogues. This signal is consistent with either a stellar companion with a moderate mass ratio ($q\sim0.5$) on a short period ($P<1$\,yr), or a substellar companion at a separation wide enough to be resolved with ground-based high contrast imaging instruments; long-period equal mass ratio stellar companions that are also consistent with the measured acceleration are excluded with previous imaging observations. The small but significant amplitude of the acceleration made HR 1645 a promising candidate for targeted searches for brown dwarf and planetary-mass companions around nearby, young stars. In this paper we explore the origin of the astrometric acceleration by modelling the signal induced by a wide-orbit M8 companion discovered with the Gemini Planet Imager, as well as the effects of an inner short-period spectroscopic companion discovered a century ago but not since followed-up. We present the first constraints on the orbit of the inner companion, and demonstrate that it is a plausible cause of the astrometric acceleration. This result demonstrates the importance of vetting of targets with measured astrometric acceleration for short-period stellar companions prior to conducting targeted direct imaging surveys for wide-orbit substellar companions.

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.

Upcoming weak-lensing surveys have the potential to become leading cosmological probes provided all systematic effects are under control. Recently, the ejection of gas due to feedback energy from active galactic nuclei (AGN) has been identified as major source of uncertainty, challenging the success of future weak-lensing probes in terms of cosmology. In this paper we investigate the effects of baryons on the number of weak-lensing peaks in the convergence field. Our analysis is based on full-sky convergence maps constructed via light-cones from $N$-body simulations, and we rely on the baryonic correction model of Schneider et al. (2019) to model the baryonic effects on the density field. As a result we find that the baryonic effects strongly depend on the Gaussian smoothing applied to the convergence map. For a DES-like survey setup, a smoothing of $\theta_k\gtrsim8$ arcmin is sufficient to keep the baryon signal below the expected statistical error. Smaller smoothing scales lead to a significant suppression of high peaks (with signal-to-noise above 2), while lower peaks are not affected. The situation is more severe for a Euclid-like setup, where a smoothing of $\theta_k\gtrsim16$ arcmin is required to keep the baryonic suppression signal below the statistical error. Smaller smoothing scales require a full modelling of baryonic effects since both low and high peaks are strongly affected by baryonic feedback.

The discovery of habitable zone (HZ) planets around low-mass stars has highlighted the need for a comprehensive understanding of the radiation environments in which such planets reside. Of particular importance is knowledge of the far-ultraviolet (FUV) radiation, as low-mass stars are typically much more active than solar-type stars and the proximity of their HZs can be one tenth the distance. The vast majority of the flux emitted by low-mass stars at FUV wavelengths occurs in the Lyman-$\alpha$ line at 1216 Angstroms. However, measuring a low-mass star's Lyman-$\alpha$ emission directly is almost always impossible because of the contaminating effects of interstellar hydrogen and geocoronal airglow. We observed Ross 825 (K3) and Ross 1044 (M0), two stars with exceptional radial velocities, with the STIS spectrograph aboard the Hubble Space Telescope (HST). Their radial velocities resulted in significant line shifts, allowing for a more complete view of their Lyman-$\alpha$ line profiles. We provide an updated relation between effective temperature and Lyman-$\alpha$ flux using Gaia DR2 astrometry as well as updated, model-independent relationships between Lyman-$\alpha$ flux and UV flux measurements from the Galaxy Evolution Explorer (GALEX) for low-mass stars. These new relations, in combination with GALEX's considerable spatial coverage, provide substantial predictive power for the Lyman-$\alpha$ environments for thousands of nearby, low-mass stars.