The first direct measurement of gravitational waves by the LIGO and Virgo collaborations has opened up new avenues to explore our Universe. This white paper outlines the challenges and gains expected in gravitational-wave searches at frequencies above the LIGO/Virgo band. The scarcity of possible astrophysical sources in most of this frequency range provides a unique opportunity to discover physics beyond the Standard Model operating both in the early and late Universe, and we highlight some of the most promising of these sources. We review several detector concepts that have been proposed to take up this challenge, and compare their expected sensitivity with the signal strength predicted in various models. This report is the summary of a series of workshops on the topic of high-frequency gravitational wave detection, held in 2019 (ICTP, Trieste, Italy), 2021 (online) and 2023 (CERN, Geneva, Switzerland).

We investigate the effect of higher-dimensional marginal operators on the thermodynamics of cosmological phase transitions. Focusing on the Abelian Higgs model, we systematically match these operators, which arise at higher orders in the underlying high-temperature expansion of thermal effective field theory, and use field redefinitions to construct a complete, minimal, and gauge-invariant operator basis. We argue that for strong transitions, temporal gauge modes, which enhance the transition strength, should be treated on equal footing with spatial gauge modes in perturbation theory. Marginal operators are found to weaken the transition strength and induce significant uncertainties for strong transitions. For even stronger transitions that could potentially produce a gravitational wave background detectable by LISA, the validity of the high temperature expansion is uncertain, which may impact the applicability of effective theory techniques, including their use in non-perturbative lattice studies.

Scientific discoveries often hinge on synthesizing decades of research, a task that potentially outstrips human information processing capacities. Large language models (LLMs) offer a solution. LLMs trained on the vast scientific literature could potentially integrate noisy yet interrelated findings to forecast novel results better than human experts. To evaluate this possibility, we created BrainBench, a forward-looking benchmark for predicting neuroscience results. We find that LLMs surpass experts in predicting experimental outcomes. BrainGPT, an LLM we tuned on the neuroscience literature, performed better yet. Like human experts, when LLMs were confident in their predictions, they were more likely to be correct, which presages a future where humans and LLMs team together to make discoveries. Our approach is not neuroscience-specific and is transferable to other knowledge-intensive endeavors.

Differential geometry offers a powerful framework for optimising and characterising finite-time thermodynamic processes, both classical and quantum. Here, we start by a pedagogical introduction to the notion of thermodynamic length. We review and connect different frameworks where it emerges in the quantum regime: adiabatically driven closed systems, time-dependent Lindblad master equations, and discrete processes. A geometric lower bound on entropy production in finitetime is then presented, which represents a quantum generalisation of the original classical bound. Following this, we review and develop some general principles for the optimisation of thermodynamic processes in the linear-response regime. These include constant speed of control variation according to the thermodynamic metric, absence of quantum coherence, and optimality of small cycles around the point of maximal ratio between heat capacity and relaxation time for Carnot engines.

We investigate the impact and mitigation of extragalactic foregrounds for the CMB lensing power spectrum analysis of Atacama Cosmology Telescope (ACT) data release 6 (DR6) data. Two independent microwave sky simulations are used to test a range of mitigation strategies. We demonstrate that finding and then subtracting point sources, finding and then subtracting models of clusters, and using a profile bias-hardened lensing estimator, together reduce the fractional biases to well below statistical uncertainties, with the inferred lensing amplitude, $A_{\mathrm{lens}}$, biased by less than $0.2\sigma$. We also show that another method where a model for the cosmic infrared background (CIB) contribution is deprojected and high frequency data from Planck is included has similar performance. Other frequency-cleaned options do not perform as well, incurring either a large noise cost, or resulting in biased recovery of the lensing spectrum. In addition to these simulation-based tests, we also present null tests performed on the ACT DR6 data which test for sensitivity of our lensing spectrum estimation to differences in foreground levels between the two ACT frequencies used, while nulling the CMB lensing signal. These tests pass whether the nulling is performed at the map or bandpower level. The CIB-deprojected measurement performed on the DR6 data is consistent with our baseline measurement, implying contamination from the CIB is unlikely to significantly bias the DR6 lensing spectrum. This collection of tests gives confidence that the ACT DR6 lensing measurements and cosmological constraints presented in companion papers to this work are robust to extragalactic foregrounds.

The advent of the James Webb Space Telescope has revealed a wealth of new galaxies just a few hundred Myr after the Big Bang. Some of these galaxies exhibit unusual N/O ratios that are difficult to explain with stellar populations today. While Wolf-Rayet stars in multiple-burst populations, very massive or rapidly-rotating primordial stars, general relativistic explosions of metal-enriched supermassive stars, or the precursors of globular clusters can in principle account for the nitrogen excess in the galaxies GN-z11 and CEERS 1019, no known stars or supernovae can explain the far higher N/O ratio of 0.46 in GS 3073 at redshift $z =$ 5.55. Here we show that the extreme nitrogen abundances in GS 3073 can be produced by 1000 - 10,000 M$_{\odot}$ primordial (Pop III) stars. We find that these are the only candidates that can account for its large N/O ratios and its C/O and Ne/O ratios. GS 3073 is thus the first conclusive evidence in the fossil abundance record of the existence of supermassive Pop III stars at cosmic Dawn.

We show that the QCD axion couples to the electromagnetic kinetic term at one loop. The result is that if axions make up dark matter, they induce temporal variation of the fine structure constant $\alpha$, which is severely constrained. We recast these constraints on the QCD axion parameter space. We also discuss how to generalise our finding to axion-like particles, and the resulting constraints.

The Lunar Gravitational-wave Antenna (LGWA) is a proposed array of next-generation inertial sensors to monitor the response of the Moon to gravitational waves (GWs). Given the size of the Moon and the expected noise produced by the lunar seismic background, the LGWA would be able to observe GWs from about 1 mHz to 1 Hz. This would make the LGWA the missing link between space-borne detectors like LISA with peak sensitivities around a few millihertz and proposed future terrestrial detectors like Einstein Telescope or Cosmic Explorer. In this article, we provide a first comprehensive analysis of the LGWA science case including its multi-messenger aspects and lunar science with LGWA data. We also describe the scientific analyses of the Moon required to plan the LGWA mission.

Seven rocky planets orbit the nearby dwarf star TRAPPIST-1, providing a unique opportunity to search for atmospheres on small planets outside the Solar System (Gillon et al., 2017). Thanks to the recent launch of JWST, possible atmospheric constituents such as carbon dioxide (CO2) are now detectable (Morley et al., 2017, Lincowski et al., 2018}. Recent JWST observations of the innermost planet TRAPPIST-1 b showed that it is most probably a bare rock without any CO2 in its atmosphere (Greene et al., 2023). Here we report the detection of thermal emission from the dayside of TRAPPIST-1 c with the Mid-Infrared Instrument (MIRI) on JWST at 15 micron. We measure a planet-to-star flux ratio of fp/fs = 421 +/- 94 parts per million (ppm) which corresponds to an inferred dayside brightness temperature of 380 +/- 31 K. This high dayside temperature disfavours a thick, CO2-rich atmosphere on the planet. The data rule out cloud-free O2/CO2 mixtures with surface pressures ranging from 10 bar (with 10 ppm CO2) to 0.1 bar (pure CO2). A Venus-analogue atmosphere with sulfuric acid clouds is also disfavoured at 2.6 sigma confidence. Thinner atmospheres or bare-rock surfaces are consistent with our measured planet-to-star flux ratio. The absence of a thick, CO2-rich atmosphere on TRAPPIST-1 c suggests a relatively volatile-poor formation history, with less than 9.5 +7.5 -2.3 Earth oceans of water. If all planets in the system formed in the same way, this would indicate a limited reservoir of volatiles for the potentially habitable planets in the system.

The measurement of Higgs couplings constitute an important part of present Standard Model precision tests at colliders. In this article, we show that modifications of Higgs couplings induce energy-growing effects in specific amplitudes involving longitudinally polarized vector bosons, and we initiate a novel program to study these very modifications of Higgs couplings off-shell and at high-energy, rather than on the Higgs resonance. Our analysis suggests that these channels are complementary and, at times, competitive with familiar on-shell measurements; moreover these high-energy probes offer endless opportunities for refinements and improvements.

We address key points for an efficient implementation of likelihood codes for modern weak lensing large-scale structure surveys. Specifically, we focus on the joint weak lensing convergence power spectrum-bispectrum probe and we tackle the numerical challenges required by a realistic analysis. Under the assumption of (multivariate) Gaussian likelihoods, we have developed a high performance code that allows highly parallelised prediction of the binned tomographic observables and of their joint non-Gaussian covariance matrix accounting for terms up to the 6-point correlation function and super-sample effects. This performance allows us to qualitatively address several interesting scientific questions. We find that the bispectrum provides an improvement in terms of signal-to-noise ratio (S/N) of about 10% on top of the power spectrum, making it a non-negligible source of information for future surveys. Furthermore, we are capable to test the impact of theoretical uncertainties in the halo model used to build our observables; with presently allowed variations we conclude that the impact is negligible on the S/N. Finally, we consider data compression possibilities to optimise future analyses of the weak lensing bispectrum. We find that, ignoring systematics, 5 equipopulated redshift bins are enough to recover the information content of a Euclid-like survey, with negligible improvement when increasing to 10 bins. We also explore principal component analysis and dependence on the triangle shapes as ways to reduce the numerical complexity of the problem.

We compare multiple foreground-cleaning pipelines for estimating the tensor-to-scalar ratio, $r$, using simulated maps of the planned CMB-S4 experiment within the context of the South Pole Deep Patch. To evaluate robustness, we analyze bias and uncertainty on $r$ across various foreground suites using map-based simulations. The foreground-cleaning methods include: a parametric maximum likelihood approach applied to auto- and cross-power spectra between frequency maps; a map-based parametric maximum-likelihood method; and a harmonic-space internal linear combination using frequency maps. We summarize the conceptual basis of each method to highlight their similarities and differences. To better probe the impact of foreground residuals, we implement an iterative internal delensing step, leveraging a map-based pipeline to generate a lensing $B$-mode template from the Large Aperture Telescope frequency maps. Our results show that the performance of the three approaches is comparable for simple and intermediate-complexity foregrounds, with $\sigma(r)$ ranging from 3 to 5 $\times 10^{-4}$. However, biases at the $1-2\sigma$ level appear when analyzing more complex forms of foreground emission. By extending the baseline pipelines to marginalize over foreground residuals, we demonstrate that contamination can be reduced to within statistical uncertainties, albeit with a pipeline-dependent impact on $\sigma(r)$, which translates to a detection significance between 2 and 4$\sigma$ for an input value of $r = 0.003$. These findings suggest varying levels of maturity among the tested pipelines, with the auto- and cross-spectra-based approach demonstrating the best stability and overall performance. Moreover, given the extremely low noise levels, mutual validation of independent foreground-cleaning pipelines is essential to ensure the robustness of any potential detection.

Deep surveys of the CMB polarization have more information on the lensing signal than the quadratic estimators (QE) can capture. We showed in a recent work that a CMB lensing power spectrum built from a single optimized CMB lensing mass map, working in close analogy to state-of-the-art QE techniques, can result in an essentially optimal spectrum estimator at reasonable numerical cost. We extend this analysis here to account for real-life non-idealities including masking and realistic instrumental noise maps. As in the QE case, it is necessary to include small corrections to account for the estimator response to these anisotropies, which we demonstrate can be estimated easily from simulations. The realization-dependent debiasing of the spectrum remains robust, allowing unbiased recovery of the band powers even in cases where the statistical model used for the lensing map reconstruction is grossly wrong. This allows now robust and at the same time optimal CMB lensing constraints from CMB data, on all scales relevant for the inference of the neutrino mass, or other parameters of our cosmological model.

An important non-perturbative effect in quantum physics is the energy gap of superconductors, which is exponentially small in the coupling constant. A natural question is whether this effect can be incorporated in the theory of resurgence. In this paper we take some steps in this direction. We conjecture that the perturbative series for the ground state energy of a superconductor is factorially divergent, and that its leading Borel singularity is governed by the superconducting energy gap. We test this conjecture in detail in the attractive Gaudin-Yang model, an exactly solvable model in one dimension with a BCS-like ground state. In order to do this, we develop techniques to calculate the exact perturbative series of its ground state energy up to high order. We also argue that the Borel singularity is of the renormalon type, and we identify a class of diagrams leading to factorial growth. We give additional evidence for the conjecture in other models.

In light of the evidence for dynamical dark energy (DE) found from the most recent Dark Energy Spectroscopic Instrument (DESI) baryon acoustic oscillation (BAO) measurements, we perform a non-parametric, model-independent reconstruction of the DE density evolution. To do so, we develop and validate a new framework that reconstructs the DE density through a third-degree piece-wise polynomial interpolation, allowing for direct constraints on its redshift evolution without assuming any specific functional form. The strength of our approach resides in the choice of directly reconstructing the DE density, which provides a more straightforward relation to the distances measured by BAO than the equation of state parameter. We investigate the constraining power of cosmic microwave background (CMB) observations combined with supernovae (SNe) and BAO measurements. In agreement with results from other works, we find a preference for deviations from $\Lambda$CDM, with a significance of $2.4\sigma$ when using the Dark Energy Survey Year 5 (DESY5) SNe data, and $1.3\sigma$ with PantheonPlus. In all the cases we consider, the derived DE equation of state parameter presents evidence for phantom crossing. By investigating potential systematic effects in the low-redshift samples of DESY5 observations, we confirm that correcting for the offset in apparent magnitude with respect to PantheonPlus data, as suggested in previous studies, completely removes the tension. Furthermore, we assess the risk of potentially overfitting the data by changing the number of interpolation nodes. As expected, we find that with lesser nodes we get a smoother reconstructed behavior of the DE density, although with similar overall features. The pipeline developed in this work is ready to be used with future high-precision data to further investigate the evidence for a non-standard background evolution.

Weak gravitational lensing of the CMB has been established as a robust and powerful observable for precision cosmology. However, the impact of Galactic foregrounds, which has been studied less extensively than many other potential systematics, could in principle pose a problem for CMB lensing measurements. These foregrounds are inherently non-Gaussian and hence might mimic the characteristic signal that lensing estimators are designed to measure. We present an analysis that quantifies the level of contamination from Galactic dust in lensing measurements, focusing particularly on measurements with the Atacama Cosmology Telescope and the Simons Observatory. We employ a whole suite of foreground models and study the contamination of lensing measurements with both individual frequency channels and multifrequency combinations. We test the sensitivity of different estimators to the level of foreground non-Gaussianity, and the dependence on sky fraction and multipole range used. We find that Galactic foregrounds do not present a problem for the Atacama Cosmology Telescope experiment (the bias in the inferred CMB lensing power spectrum amplitude remains below $0.3\sigma$). For Simons Observatory, not all foreground models remain below this threshold. Although our results are conservative upper limits, they suggest that further work on characterizing dust biases and determining the impact of mitigation methods is well motivated, especially for the largest sky fractions.

We study the resurgent structure of the refined topological string partition function on a non-compact Calabi-Yau threefold, at large orders in the string coupling constant $g_s$ and fixed refinement parameter $\mathsf{b}$. For $\mathsf{b}\neq 1$, the Borel transform admits two families of simple poles, corresponding to integral periods rescaled by $\mathsf{b}$ and $1/\mathsf{b}$. We show that the corresponding Stokes automorphism is expressed in terms of a generalization of the non-compact quantum dilogarithm, and we conjecture that the Stokes constants are determined by the refined Donaldson-Thomas invariants counting spin-$j$ BPS states. This jump in the refined topological string partition function is a special case (unit five-brane charge) of a more general transformation property of wave functions on quantum twisted tori introduced in earlier work by two of the authors. We show that this property follows from the transformation of a suitable refined dual partition function across BPS rays, defined by extending the Moyal star product to the realm of contact geometry.

We show that the all-orders WKB periods of one-dimensional quantum mechanical oscillators are governed by the refined holomorphic anomaly equations of topological string theory. We analyze in detail the double-well potential and the cubic and quartic oscillators, and we calculate the WKB expansion of their quantum free energies by using the direct integration of the anomaly equations. We reproduce in this way all known results about the quantum periods of these models, which we express in terms of modular forms on the WKB curve. As an application of our results, we study the large order behavior of the WKB expansion in the case of the double well, which displays the double factorial growth typical of string theory.

We study lensing of gravitational waves by a black hole in the deep wave optics regime, i.e. when the wavelength is much larger than the black hole Schwarzschild radius. We apply it to triple systems, with a binary of stellar mass objects in the inspiraling phase orbiting around a central massive black hole. We describe the full polarisation structure of the wave and derive predictions for the polarisation modes of the scattered wave measured by the observer. We show that lensing in the wave optics regime is not helicity preserving, as opposed to lensing in the geometric optics regime. The amplitude of the total wave is modulated due to interference between the directly transmitted and lensed components. The relative amplitude of the modulation is fixed by the lensing geometry and can reach unity in the most favourable settings. This indicates that wave optics lensing is potentially detectable by LISA for sufficiently high SNR systems. Our findings show that in the wave optics regime it is necessary to go beyond the usual lensing description where the amplification factor is assumed to be the same for both helicity modes. While motivated by GW190521 and the AGN formation scenario, our results apply more broadly to stellar-mass binaries orbiting a third body described as a Schwarzschild black hole, with a period comparable to the GW observation time.