Leptoquark models are prime candidates for new physics (NP) explanations of the long-standing anomalies in semi-leptonic $B$ decays; $b\to c \tau \bar\nu$ (encoded in $R(D^{(\ast)})$) and $b\to s\ell\bar\ell (\ell=e,\mu)$ transitions. Furthermore, Belle II and NA62 reported weaker-than-expected limits on $B^+ \to K^+ \nu\bar\nu$ and $K^+ \to \pi^+ \nu \bar\nu$, respectively. While the $R(D^{(\ast)})$ and $b\to s\ell \bar\ell$ measurements can be explained with NP contributions at the $O(10\%)$ level, the neutrino channels suggest that the NP effect could be comparable in size to the Standard Model one. In this context, we consider the two types of leptoquark models with minimal sets of the couplings that can best describe the semi-leptonic $B$ anomalies and lead at the same time to effects in the neutrino modes, the singlet-triplet scalar leptoquark model ($S_1+S_3$) and the singlet vector leptoquark model ($U_1$). More specifically, the neutrino channels pose non-trivial constraints on the parameter space, and we find that large effects (i.e., accounting for the current central value) in $B\to K^{(*)}\nu\bar\nu$ are only possible in the $S_1+S_3$ setup, while both models can account for the central value of $K^+\to \pi^+\nu\bar\nu$.

We compute tree-level $n$-point scattering amplitudes in scalar field theories in terms of geometric invariants on a fibre bundle. All 0- and 2-derivative interactions are incorporated into a metric on this bundle. The on-shell amplitudes can be efficiently pieced together from covariant Feynman rules, and we present a general closed formula for obtaining the $n$-point amplitude in this way. The covariant Feynman rules themselves can be derived using a generalization of the normal coordinate expansion of the fibre bundle metric. We demonstrate the efficiency of this approach by computing the covariant Feynman rules up to $n=10$ points, from which one can obtain the full amplitudes using our general formula. The formalism offers a prototype for obtaining geometric amplitudes in theories with higher-derivative interactions, by passing from the fibre bundle to its jet bundles.

The Real Higgs Triplet model, known as the $\Delta$SM, is a minimal extension of the Standard Model (SM) obtained by adding a hypercharge 0 triplet ($\Delta$). This simple model is motivated by the multi-lepton anomalies and excesses in di-photon, $Z\gamma$, and $WW$ spectra at $\approx152\,\text{GeV}$. The model contains, in addition to the SM particle content, a $CP$-even neutral Higgs ($\Delta^0$) and a charged state ($\Delta^\pm$), which are quasi-degenerate in mass. Observing the charged scalar at the LHC and measuring its mass is very challenging, since it dominantly decays to $WZ$, $tb$, and $\tau\nu$. In this article, we consider the discovery prospects at future electron-positron colliders. Taking into account $e^+e^- \to \gamma^*,Z^* \to \Delta^\pm \Delta^\mp$ as the production mechanism and the dominant decay modes, we define three signal regions (SR) to study the charged Higgs properties: SR1: $\ge 3j + 1\ell$, SR2: $\ge 3\ell + \tau_{\text{had}}$, SR3: $\ge 4j + \tau_{\text{had}}$. For $m_{\Delta^\pm}\approx150\,\text{GeV}$, a $5\sigma$ significance can be achieved in SR1 with an integrated luminosity of less than $1\text{ fb}^{-1}$. SR2 is very clean with leptonic final states having low background and small systematic uncertainties. Furthermore, SR3 is crucial for reconstructing the charged scalar invariant mass, which can be measured with $\mathcal{O}(1)\,\text{GeV}$ accuracy with an integrated luminosity of $500\text{ fb}^{-1}$.

We calculate the two-mass three-loop contributions to the unpolarized and polarized massive operator matrix elements $\tilde{A}_{Qg}^{(3)}$ and $\Delta \tilde{A}_{Qg}^{(3)}$ in $x$-space for a general mass ratio by using a semi-analytic approach. We also compute Mellin moments up to $N = 2000 (3000)$ by an independent method, to which we compare the results in $x$-space. In the polarized case, we work in the Larin scheme. We present numerical results. The two-mass contributions amount to about $50 \%$ of the full \textcolor{blue}{$O(T_F^2)$} and \textcolor{blue}{$O(T_F^3)$} terms contributing to the operator matrix elements. The present result completes the calculation of all unpolarized and polarized massive three-loop operator matrix elements.

We present $100$ full-sky quasar spectrophotometric mock catalogs with smooth redshift evolution from $z=0$ to $z\sim 4$, tailored to analyze the Gaia-unWISE Quasar Catalog (Quaia). In particular, we apply a novel hierarchical nonlocal nonlinear bias scheme (Hicobian) to dark matter fields generated through Augmented Lagrangian Perturbation Theory on the lightcone (WebON code), calibrating the free parameters of the bias model on Abacus quasar HOD mock catalogs tuned to reproduce DESI Early Data Release observations in real and redshift space. After having obtained such accurate spectroscopic catalogs, we inject in the mocks the observational effects characterizing the Quaia catalog: (i) spectrophotometric redshift uncertainties, (ii) the angular selection function, and (iii) the redshift number counts distribution. We assess the accuracy of our catalogs by validating a number of summary statistics: the full-sky QSO maps, the redshift uncertainty distributions as a function of redshift, the redshift $n(z)$ distribution, the angular power spectra and their normalized covariance matrices, and the angular two-point correlation functions. We find excellent agreement between these metrics from the mocks and from the Quaia catalog. We publicly release the mock catalogs to the community.

A novel method to treat effects from evanescent operators in next-to-leading order (NLO) computations is introduced. The approach allows, besides further simplifications, to discard evanescent-to-physical mixing contributions in NLO calculations. The method is independent of the treatments of $\gamma_5$ and can therefore be combined with different renormalization schemes. We illustrate the utility of this result by reproducing literature results of two-loop anomalous dimension matrices for both $|\Delta F| = 1$ and $|\Delta F| = 2$ transitions.

We present a public catalogue of numerical-relativity binary-black-hole simulations. The catalogue contains datasets from 80 distinct configurations of precessing binary-black-hole systems, with mass ratios up to $m_2/m_1 = 8$, dimensionless spin magnitudes on the larger black hole up to $|\vec{S}_2|/m_2^2 = 0.8$ (the small black hole is non-spinning), and a range of five values of spin misalignment for each mass-ratio/spin combination. We discuss the physical properties of the configurations in our catalogue, and assess the accuracy of the initial configuration of each simulation and of the gravitational waveforms. We perform a careful analysis of the errors due to the finite resolution of our simulations and the finite distance from the source at which we extract the waveform data and provide a conservative estimate of the mismatch accuracy. We find that the upper limit on the mismatch uncertainty of our waveforms is $0.4\%$. In doing this we present a consistent approach to combining mismatch uncertainties from multiple error sources. We compare this release to previous catalogues and discuss how these new simulations complement the existing public datasets. In particular, this is the first catalogue to uniformly cover this parameter space of single-spin binaries and there was previously only sparse coverage of the precessing-binary parameter space for mass ratios $\gtrsim 5$. We discuss applications of these new data, and the most urgent directions for future simulation work. The public dataset can be accessed online at this https URL.

Measuring the Higgs boson couplings as precisely as possible is one of the major goals of the High Luminosity LHC. We show that the $(b\bar{b})(b\bar{b})$ final state in Higgs boson pair production can be exploited in the boosted regime to give constraints on the trilinear Higgs boson self-coupling. In these exclusive phase space regions, novel jet substructure techniques can be used to separate the signal from the large QCD and electroweak backgrounds. New developments on trigger and b-tagging strategies for the upcoming LHC runs are necessary in order to reconstruct the Higgs bosons in boosted final states, where the trilinear self-coupling sensitivity is reduced. We find that using our approach one can set a limit for $\lambda \leq 1.2$ at $95 \%$ CL after $3000~\mathrm{fb}^{-1}$. As the signal-to-background ratio is small we propose a data-driven side-band analysis to improve on the coupling measurement.

We consider 2-to-2 scatterings of Higgs bosons in a CP-conserving two-Higgs-doublet model (2HDM) and study the implication of maximizing the entanglement in the flavor space, where the two doublets $\Phi_a$, $a=1,2$, can be viewed as a qubit: $\Phi_1=|0\rangle$ and $\Phi_2=|1\rangle$. More specifically, we compute the scattering amplitudes for $\Phi_a \Phi_b \to \Phi_c \Phi_d$ and require the outgoing flavor entanglement to be maximal for a full product basis such as the computational basis, which consists of $\{|00\rangle,|01\rangle,|10\rangle,|11\rangle\}$. In the unbroken phase and turning off the gauge interactions, entanglement maximization results in the appearance of an $U(2)\times U(2)$ global symmetry among the quartic couplings, which in general is broken softly by the mass terms. Interestingly, once the Higgs bosons acquire vacuum expectation values, maximal entanglement enforces an exact $U(2) \times U(2)$ symmetry, which is spontaneously broken to $U(1)\times U(1)$. As a byproduct, this gives rise to Higgs alignment as well as to the existence of 6 massless Nambu-Goldstone bosons. The $U(2)\times U(2)$ symmetry can be gauged to lift the massless Goldstones, while maintaining maximal entanglement demands the presence of a discrete $\mathrm{Z}_2$ symmetry interchanging the two gauge sectors. The model is custodially invariant in the scalar sector, and the inclusion of fermions requires a mirror dark sector, related to the standard one by the $\mathrm{Z}_2$ symmetry.

The presence of scalar fields with non-minimal gravitational interactions of the form $\xi |\phi|^2 R$ may have important implications for the physics of the early universe. While many studies solve the dynamics of non-minimally coupled scalars in the Einstein frame, where gravity is simply described by the Einstein-Hilbert action, we instead propose a procedure to solve the dynamics directly in the original Jordan frame where the non-minimal couplings are maintained explicitly. Our algorithm can be applied to scenarios that include minimally coupled fields and an arbitrary number of non-minimally coupled scalars, with the expansion of the universe sourced by all fields present. This includes situations when the dynamics become fully inhomogeneous, fully non-linear (due to e.g.~backreaction or mode rescattering effects), and/or when the expansion of the universe is dominated by non-minimally coupled species. As an example, we study geometric preheating with a non-minimally coupled scalar spectator field when the inflaton oscillates following the end of inflation. In the future, our technique may be used to shed light on aspects of the equivalence of the Jordan and Einstein frames at the quantum level.

We consider mixed strong-electroweak corrections to Higgs production via gluon fusion, in which the Higgs boson couples to the top quark. Using the method of differential equations, we compute all of the master integrals that contribute to this process at two loops through $\mathcal{O}(\epsilon^2)$ in the dimensional regularization parameter $\epsilon = (d-4)/2$, keeping full analytic dependence on the top quark, Higgs, W, and Z boson masses. We present the results for these master integrals in terms of iterated integrals whose kernels depend on elliptic curves.

JWST observations have uncovered a new population of red, compact objects at high redshifts dubbed `Little Red Dots' (LRDs), which typically show broad emission lines and are thought to be dusty Active Galactic Nuclei (AGN). Some of their other features, however, challenge the AGN explanation, such as prominent Balmer breaks and extremely faint or even missing metal high-ionization lines, X-ray, or radio emission, including in deep stacks. Time variability is another, robust, test of AGN activity. Here, we exploit the $z=7.045$ multiply-imaged LRD A2744-QSO1, which offers a particularly unique test of variability due to lensing-induced time delays between the three images spanning 22 yr (2.7 yr in the rest-frame), to investigate its photometric and spectroscopic variability. We find the equivalent widths (EWs) of the broad H$\alpha$ and H$\beta$ lines, which are independent of magnification and other systematics, to exhibit significant variations, up to $18\pm3$ % for H$\alpha$ and up to $22\pm8$ % in H$\beta$, on a timescale of 875 d (2.4 yr) in the rest-frame. This suggests that A2744-QSO1 is indeed an AGN. We find no significant photometric variability beyond the limiting systematic uncertainties, so it currently cannot be determined whether the EW variations are due to line-flux or continuum variability. These results are consistent with a typical damped random walk (DRW) variability model for an AGN like A2744-QSO1 ($M_{\mathrm{BH}}=4\times10^7 \mathrm{M}_{\odot}$) given the sparse sampling of the light-curve with the available data. Our results therefore support the AGN interpretation of this LRD, and highlight the need for further photometric and spectroscopic monitoring in order to build a detailed and reliable light-curve.

We applied the global end-to-end model described in Paper~I of this series to perform a population synthesis of companions formed via disc instability (DI). By using initial conditions compatible with both observations and hydrodynamical simulations, and by studying a large range of primary masses (0.05-5 Msol), we can provide quantitative predictions of the outcome of DI. In the baseline population, we find that ~10 % of the discs fragment, and about half of these end up with a surviving companion after 100 Myr. 75\% of the companions are in the brown dwarf regime, 15 % are low-mass stars, and 10 % planets. At distances larger than ~100 au, DI produces planetary-mass companions on a low percent level. Inside of 100 AU, however, planetary-mass companions are very rare (low per mill level). The average companion mass is ~30 Mj scaling weakly with stellar mass. Most of the initial fragments do not survive on a Myr timescale; they either collide with other fragments or are ejected, resulting in a population of free-floating objects (about 1-2 per star). We also quantify several variant populations to critically assess some of our assumptions used in the baseline population. DI appears to be a key mechanism in the formation of distant companions with masses ranging from low-mass stars down to the planetary regime, contributing, however, only marginally to planetary mass objects inside of 100 AU. Our results are sensitive to a number of physical processes, which are not completely understood. Two of them, gas accretion and clump-clump collisions, are particularly important and need to be investigated further. Magnetic fields and heavy-element accretion have not been considered in our study, although they are also expected to affect the inferred population. We suggest acknowledging the importance of the gravito-turbulent phase, which most protoplanetary discs experience.

The cosmological redshift drift promises to be the first observable directly measuring the evolution of the cosmic expansion rate and should be detectable with upcoming surveys by the Square Kilometre Array and the Extremely Large Telescope. To prepare for these upcoming measurements we study the redshift drift in detail using the relativistic N-body code $\texttt{gevolution}$, focusing on inhomogeneity-induced fluctuations. Using a ray-tracer, we calculate the redshift drift directly from the light cone at two different time steps. To investigate observer-dependent biases we consider 10 different observers. We find that inhomogeneity-induced fluctuations in the redshift drift can in extreme cases be of the same order as the cosmic signal for $z\lesssim0.15$. By comparing our results to first-order perturbation theory, we find that the extreme outliers are due to peculiar motion in over-densities and can be described by first-order perturbation theory to percent precision. We calculate angular power spectra that fit very well with our predictions based on perturbation theory at linear scales and show a surprisingly large non-linear signal. This shows that redshift drift not only has the power to measure the background expansion, but could also deliver information about the velocity and acceleration fields in clusters.

The precision of fixed-order calculations on the dihadron production in electron-positron annihilation is paramount for probing QCD factorization and constraining non-perturbative inputs. This paper investigates the QCD corrections to the angular separation distribution $\theta_{12}$ between two observed hadrons, $H_1$ and $H_2$, in the process $e^+e^- \to H_1 H_2 + X$ up to $\mathcal{O}(\alpha_s^2)$, with particular emphasis on the intermediate region $\theta_{12} \in (0,\pi)$. The partonic processes at this accuracy consist of two sorts of contributions, the real-virtual and double-real corrections. Of them, the evaluation of four-body phase space integrals in the latter case is at the core of this study. To address them, we first employ the integration-by-parts (IBP) identities to reduce the number of independent integrals and then apply the differential equations (DE) method to recursively solve the resulting master integrals. In kinematic regions where the invariant mass of the unresolved partons vanishes, IBP coefficients can develop divergences. To this end, we resum higher-order terms in the dimensional regulator for each master integral based on the asymptotic behavior of the canonical DEs. After combining the real and virtual corrections with the counter terms from fragmentation function renormalization, we demonstrate that the pole terms in the final analytic expressions exactly cancel out in all partonic channels, thereby providing a non-trivial validation of collinear factorization at the next-to-leading order (NLO). Eventually, when presenting our analytic expressions of the finite partonic coefficients, we transform the transcendental functions resulting from the DE solutions into classical (poly)logarithmic functions, in order to facilitate the implementation in event generators.

We conduct the first full Bayesian inference analysis for LISA parameter estimation incorporating the effects of subdominant harmonics and spin-precession through a full time domain response. The substantial computational demands of using time domain waveforms for LISA are significantly mitigated by implementing a novel Python version of the IMRPhenomT family of waveform models and the LISA response with GPU acceleration. This time domain response alleviates the theoretical necessity of developing specific transfer functions to approximate the LISA response in the Fourier domain for each specific type of system and allows for the use of unequal arms configurations and realistic LISA orbits. Our analysis includes a series of zero-noise injections for a Massive Black Hole Binary with aligned and precessing spins. We investigate the impact of including subdominant harmonics, compare equal and unequal arm configurations, and analyze different Time-Delay-Interferometry (TDI) configurations. We utilize full and uniform priors, with a lower frequency cutoff of 0.1mHz, and a signal duration of approximately two months, sampled every 5 seconds. The sampler is initialized based on Fisher estimates. Our results demonstrate LISA capability to measure the two spin magnitudes and the primary spin tilt angle, alongside sky localization, with percent-level precision, while component masses are determined with sub-percent accuracy.

We present a three-site Pati-Salam gauge model able to explain the Standard Model flavor hierarchies while, at the same time, accommodating the recent experimental hints of lepton-flavor non-universality in $B$ decays. The model is consistent with low- and high-energy bounds, and predicts a rich spectrum of new states at TeV scale that could be probed in the near future by the high-$p_T$ experiments at the LHC.

There are two tensions related to the Cabibbo angle of the CKM matrix. First, the determinations of $V_{us}$ from $K_{\mu 2}$, $K_{\ell3}$, and $\tau$ decays disagree at the $3\sigma$ level. Second, using the average of these results in combination with $\beta$ decays (including super-allowed $\beta$ decays and neutron decay), a deficit in first-row CKM unitarity with a significance of again about $3\sigma$ is found. These discrepancies, known as the Cabibbo Angle anomaly, can in principle be solved by modifications of $W$ boson couplings to quarks. However, due to $SU(2)_L$ invariance, $Z$ couplings to quarks are also modified and flavour changing neutral currents can occur. In order to consistently assess the agreement of a new physics hypothesis with data, we perform a combined analysis for all dimension-six Standard Model Effective Field Theory operators that generate modified $W$ couplings to first and second generation quarks. We then study models with vector-like quarks, which are prime candidates for a corresponding UV completion as they can affect $W$-quark couplings at tree level, and we perform a global fit including flavour observables (in particular loop effects in $\Delta F=2$ processes). We find that the best fit can be obtained for the $SU(2)_L$ doublet vector-like quark $Q$ as it can generate right-handed $W$-$u$-$d$ and $W$-$u$-$s$ couplings as preferred by data.

Studying the impact of new-physics models on low-energy observables necessitates matching to effective field theories at the relevant mass thresholds. We introduce the first public version of Matchete, a computer tool for matching weakly-coupled models at one-loop order. It uses functional methods to directly compute all matching contributions in a manifestly gauge-covariant manner, while simplification methods eliminate redundant operators from the output. We sketch the workings of the program and provide examples of how to match simple Standard Model extensions. The package, documentation, and example notebooks are publicly available at this https URL