An axion-like field comprising $\sim 10\%$ of the energy density of the universe near matter-radiation equality is a candidate to resolve the Hubble tension; this is the "early dark energy" (EDE) model. However, as shown in Hill et al. (2020), the model fails to simultaneously resolve the Hubble tension and maintain a good fit to both cosmic microwave background (CMB) and large-scale structure (LSS) data. Here, we use redshift-space galaxy clustering data to sharpen constraints on the EDE model. We perform the first EDE analysis using the full-shape power spectrum likelihood from the Baryon Oscillation Spectroscopic Survey (BOSS), based on the effective field theory (EFT) of LSS. The inclusion of this likelihood in the EDE analysis yields a $25\%$ tighter error bar on $H_0$ compared to primary CMB data alone, yielding $H_0 = 68.54^{+0.52}_{-0.95}$ km/s/Mpc ($68\%$ CL). In addition, we constrain the maximum fractional energy density contribution of the EDE to $f_{\rm EDE} < 0.072$ ($95\%$ CL). We explicitly demonstrate that the EFT BOSS likelihood yields much stronger constraints on EDE than the standard BOSS likelihood. Including further information from photometric LSS surveys,the constraints narrow by an additional $20\%$, yielding $H_0 = 68.73^{+0.42}_{-0.69}$ km/s/Mpc ($68\%$ CL) and $f_{\rm EDE}<0.053$ ($95\%$ CL). These bounds are obtained without including local-universe $H_0$ data, which is in strong tension with the CMB and LSS, even in the EDE model. We also refute claims that MCMC analyses of EDE that omit SH0ES from the combined dataset yield misleading posteriors. Finally, we demonstrate that upcoming Euclid/DESI-like spectroscopic galaxy surveys can greatly improve the EDE constraints. We conclude that current data preclude the EDE model as a resolution of the Hubble tension, and that future LSS surveys can close the remaining parameter space of this model.

We present a search for high-energy $\gamma$-ray emission from 566 Active Galactic Nuclei at redshift $z > 0.2$, from the 2WHSP catalog of high-synchrotron peaked BL Lac objects with eight years of Fermi-LAT data. We focus on a redshift range where electromagnetic cascade emission induced by ultra-high-energy cosmic rays can be distinguished from leptonic emission based on the spectral properties of the sources. Our analysis leads to the detection of 160 sources above $\approx$ $5\sigma$ (TS $\geq 25$) in the 1 - 300 GeV energy range. By discriminating significant sources based on their $\gamma$-ray fluxes, variability properties, and photon index in the Fermi-LAT energy range, and modeling the expected hadronic signal in the TeV regime, we select a list of promising sources as potential candidate ultra-high-energy cosmic-ray emitters for follow-up observations by Imaging Atmospheric Cherenkov Telescopes.

Current cosmological data exhibit a tension between inferences of the Hubble constant, $H_0$, derived from early and late-universe measurements. One proposed solution is to introduce a new component in the early universe, which initially acts as "early dark energy" (EDE), thus decreasing the physical size of the sound horizon imprinted in the cosmic microwave background (CMB) and increasing the inferred $H_0$. Previous EDE analyses have shown this model can relax the $H_0$ tension, but the CMB-preferred value of the density fluctuation amplitude, $\sigma_8$, increases in EDE as compared to $\Lambda$CDM, increasing tension with large-scale structure (LSS) data. We show that the EDE model fit to CMB and SH0ES data yields scale-dependent changes in the matter power spectrum compared to $\Lambda$CDM, including $10\%$ more power at $k = 1~h$/Mpc. Motivated by this observation, we reanalyze the EDE scenario, considering LSS data in detail. We also update previous analyses by including $Planck$ 2018 CMB likelihoods, and perform the first search for EDE in $Planck$ data alone, which yields no evidence for EDE. We consider several data set combinations involving the primary CMB, CMB lensing, SNIa, BAO, RSD, weak lensing, galaxy clustering, and local distance-ladder data (SH0ES). While the EDE component is weakly detected (3$\sigma$) when including the SH0ES data and excluding most LSS data, this drops below 2$\sigma$ when further LSS data are included. Further, this result is in tension with strong constraints imposed on EDE by CMB and LSS data without SH0ES, which show no evidence for this model. We also show that physical priors on the fundamental scalar field parameters further weaken evidence for EDE. We conclude that the EDE scenario is, at best, no more likely to be concordant with all current cosmological data sets than $\Lambda$CDM, and appears unlikely to resolve the $H_0$ tension.

Strong gravitational lensing is a promising probe of the substructure of dark matter halos. Deep learning methods have the potential to accurately identify images containing substructure, and differentiate WIMP dark matter from other well motivated models, including vortex substructure of dark matter condensates and superfluids. This is crucial in future efforts to identify the true nature of dark matter. We implement, for the first time, a classification approach to identifying dark matter substructure based on simulated strong lensing images with different substructure. Utilizing convolutional neural networks trained on sets of simulated images, we demonstrate the feasibility of deep neural networks to reliably distinguish among different types of dark matter substructure. With thousands of strong lensing images anticipated with the coming launch of LSST, we expect that supervised and unsupervised deep learning models will play a crucial role in determining the nature of dark matter.

The recent detection of gravitational waves and electromagnetic counterparts from the double neutron star merger event GW+EM170817, supports the standard paradigm of short gamma-ray bursts (SGRBs) and kilonovae/macronovae. It is important to reveal the nature of the compact remnant left after the merger, either a black hole or neutron star, and their physical link to the origin of the long-lasting emission observed in SGRBs. The diversity of the merger remnants may also lead to different kinds of transients that can be detected in future. Here we study the high-energy emission from the long-lasting central engine left after the coalescence, under certain assumptions. In particular, we consider the X-ray emission from a remnant disk and the non-thermal nebular emission from disk-driven outflows or pulsar winds. We demonstrate that late-time X-ray and high-frequency radio emission can provide useful constraints on properties of the hidden compact remnants and their connections to long-lasting SGRB emission, and we discuss the detectability of nearby merger events through late-time observations at $\sim30-100$ d after the coalescence. We also investigate the GeV-TeV gamma-ray emission that occurs in the presence of long-lasting central engines, and show the importance of external inverse-Compton radiation due to up-scattering of X-ray photons by relativistic electrons in the jet. We also search for high-energy gamma-rays from GW170817 in the Fermi-LAT data, and report upper limits on such long-lasting emission. Finally, we consider the implications of GW+EM170817 and discuss the constraints placed by X-ray and high-frequency radio observations.

We present a targeted search for blazar flux-correlated high-energy ($\varepsilon_\nu > 1$ TeV) neutrinos from six bright northern blazars, using the public database of northern-hemisphere neutrinos detected during "IC40" 40-string operations of the IceCube neutrino observatory (April 2008 to May 2009). Our six targeted blazars are subjects of long-term monitoring campaigns by the VERITAS TeV gamma-ray observatory. We use the publicly-available VERITAS lightcurves to identify periods of excess and flaring emission. These predefined intervals serve as our "active temporal windows" in a search for an excess of neutrinos, relative to Poisson fluctuations of the near-isotropic atmospheric neutrino background which dominates at these energies. After defining the parameters of an optimized search, we confirm the expected Poisson behavior with Monte Carlo simulations prior to testing for excess neutrinos in the actual data. We make two searches: One for excess neutrinos associated with the bright flares of Mrk 421 that occurred during the IC40 run, and one for excess neutrinos associated with the brightest emission periods of five other blazars (Mrk 501, 1ES 0805+524, 1ES 1218+304, 3C66A, and W Comae), all significantly fainter than the Mrk 421 flares. We find no significant excess of neutrinos from the preselected blazar directions during the selected temporal windows. We derive 90%-confidence upper limits on the number of expected flux-associated neutrinos from each search. These limits are consistent with previous point-source searches and Fermi GeV flux-correlated searches. Our upper limits are sufficiently close to the physically-interesting regime that we anticipate future analyses using already-collected data will either constrain models or yield discovery of the first blazar-associated high-energy neutrinos.