Precise modeling of detector energy response is crucial for next-generation neutrino experiments which present computational challenges due to lack of analytical likelihoods. We propose a solution using neural likelihood estimation within the simulation-based inference framework. We develop two complementary neural density estimators that model likelihoods of calibration data: conditional normalizing flows and a transformer-based regressor. We adopt JUNO - a large neutrino experiment - as a case study. The energy response of JUNO depends on several parameters, all of which should be tuned, given their non-linear behavior and strong correlations in the calibration data. To this end, we integrate the modeled likelihoods with Bayesian nested sampling for parameter inference, achieving uncertainties limited only by statistics with near-zero systematic biases. The normalizing flows model enables unbinned likelihood analysis, while the transformer provides an efficient binned alternative. By providing both options, our framework offers flexibility to choose the most appropriate method for specific needs. Finally, our approach establishes a template for similar applications across experimental neutrino and broader particle physics.

In scalar-tensor Horndeski theories, nonsingular cosmological models - bounce and genesis - are problematic because of potential ghost and/or gradient instabilities. One way to get around this obstacle is to send the effective Planck mass to zero in the asymptotic past ("strong gravity in the past"). One may suspect that this feature is a signal of a strong coupling problem at early times. However, the classical treatment of the cosmological background is legitimate, provided that the strong coupling energy scale remains at all times much higher than the scale associated with the classical evolution. We construct various models of this sort, namely (i) bouncing Universe which proceeds through inflationary epoch to kination (expansion within general relativity, driven by massless scalar field); (ii) bouncing Universe with kination stage immediately after bounce; (iii) combination of genesis and bounce, with the Universe starting from flat space-time, then contracting and bouncing to the expansion epoch; (iv) "standard" genesis evading the strong coupling problem in the past. All these models are stable, and perturbations about the backgrounds are not superluminal.

We construct a concrete model of Horndeski bounce with strong gravity in the past. Within this model we show that the correct spectra of cosmological perturbations may be generated at early contracting epoch, with mild fine-tuning ensuring that the scalar spectral tilt $n_S$ and tensor-to-scalar ratio $r$ are consistent with observations. The smallness of $r$ is governed by the smallness of the scalar sound speed. Arbitrarily small values of $r$ are forbidden in our setup because of the strong coupling in the past. Nevertheless, we show that it is possible to generate perturbations in a controllable way, i.e. in the regime where the background evolution and perturbations are legitimately described within classical field theory and weakly coupled quantum theory.

In this note the finite formulation of quantum field theory, which is based on the system of the differential equations which are reminiscent of Callan-Symanzik equations is discussed. This system of equations was previously formulated on the bare language. We rederive these equations on a fully renormalized language. In this language, it was demonstrated for a simple $\phi^4$ toy model, that with the specific choice of renormalization conditions - on-shell scheme for renormalized mass - the class of such finite renormalization prescriptions is equivalent to the classical renormalization group equation written in Callan-Symanzik-Ovsyannikov form.

We study whether a violation of the null energy condition necessarily implies the presence of instabilities. We prove that this is the case in a large class of situations, including isotropic solids and fluids relevant for cosmology. On the other hand we present several counter-examples of consistent effective field theories possessing a stable background where the null energy condition is violated. Two necessary features of these counter-examples are the lack of isotropy of the background and the presence of superluminal modes. We argue that many of the properties of massive gravity can be understood by associating it to a solid at the edge of violating the null energy condition. We briefly analyze the difficulties of mimicking $\dot H>0$ in scalar tensor theories of gravity.

We present a new type of Scalar-Vector-Tensor (SVT) theories with higher derivatives of all the fields in the action, but with second order equations of motion. The higher derivative vector field is invariant under a U(1) gauge transformation and the Scalar-Tensor sector corresponds to Horndeski theory. We also present a subclass of these SVT theories with 8 free functions of $\pi$ and $X$ where the speed of the tensor and vector modes is exactly the same. In particular, the Horndeski functions $G_4(\pi,X)$ and $G_5(\pi)$ remain free, while the speed of the vector modes tracks the speed of the tensor modes. Additionally, the vector sector retains freedom through the four new functions. All the theories here shown are a generalization of the Kaluza-Klein reduction of 5D Horndeski theory, sharing the main properties in cosmology, but including new free scalar functions in the Lagrangian.

We review recent progress in the study of the anomalous baryon number non-conservation at high temperatures and in high energy collisions. Recent results on high temperature phase transitions are described, and applications to electroweak baryogenesis are considered. The current status of the problem of electroweak instanton-like processes at high energies is outlined. This paper is written on the occasion of Sakharov's 75th anniversary and will appear in the memorial volume of Uspekhi (Usp. Fiz. Nauk, volume 166, No 5, May 1996).

Evolution of cosmic domain walls (DWs) settles to the scaling solution, which is often assumed to be independent of initial conditions. However, lattice simulations performed in this work reveal a clear dependence of the scaling DW area on the initial configuration of the sourcing scalar field, specifically, its infrared (IR) properties. Namely, the DW area grows as one suppresses IR modes in the initial scalar field spectrum. This growth is saturated, when the area parameter $\xi$ commonly used in the literature reaches the value $\xi_{max} \approx 1.2$. The dependence of $\xi$ on IR modes is argued to be of non-physical origin: it is likely to be due to effects of the lattice boundary. Assuming that physically the memory of initial conditions is erased, one recognizes $\xi \approx 1.2$ obtained in the situation with maximally suppressed IR modes as a genuine universal value of the area parameter in the scaling regime. We demonstrate that ignorance about initial conditions may affect predictions for the energy density of gravitational waves by the factor five. The spectral shape of gravitational waves is also affected by the choice of initial conditions, most notably in the low-frequency part. Likewise, we revisit annihilation of DWs under the influence of a potential bias. It has been previously found in Ref. [19] that the annihilation happens significantly earlier compared to the estimate based on the simple balance between the potential bias and surface energy density. We further support this observation and show that the tendency towards an earlier annihilation gets even stronger upon removing IR modes in simulations.

We study the production of hypothetical vector portal mediators, dark photons $\gamma^\prime$, with masses in the range 0.4-3.5 GeV in the negatively charged pion-proton collisions $\pi^-p\rightarrow \gamma^\prime X$ via inelastic pion bremsstrahlung and QCD Drell-Yan-like process. In both cases we estimate the value of total dark photon production cross section and obtain the energy distribution for dark photons that could be produced in the NA64h experiment. We also present the mean energies of dark photons produced in the same way by the secondary pions with momenta typical for T2K, DUNE and SHiP experiments.

Based on $10.64~\mathrm{fb}^{-1}$ of $e^+e^-$ collision data taken at center-of-mass energies between 4.237 and 4.699 GeV with the BESIII detector, we study the leptonic $D^+_s$ decays using the $e^+e^-\to D^{*+}_{s} D^{*-}_{s}$ process. The branching fractions of $D_s^+\to\ell^+\nu_{\ell}\,(\ell=\mu,\tau)$ are measured to be $\mathcal{B}(D_s^+\to\mu^+\nu_\mu)=(0.547\pm0.026_{\rm stat}\pm0.016_{\rm syst})\%$ and $\mathcal{B}(D_s^+\to\tau^+\nu_\tau)=(5.60\pm0.16_{\rm stat}\pm0.20_{\rm syst})\%$, respectively. The product of the decay constant and Cabibbo-Kobayashi-Maskawa matrix element $|V_{cs}|$ is determined to be $f_{D_s^+}|V_{cs}|=(246.5\pm5.9_{\rm stat}\pm3.6_{\rm syst}\pm0.5_{\rm input})_{\mu\nu}~\mathrm{MeV}$ and $f_{D_s^+}|V_{cs}|=(252.7\pm3.6_{\rm stat}\pm4.5_{\rm syst}\pm0.6_{\rm input}))_{\tau \nu}~\mathrm{MeV}$, respectively. Taking the value of $|V_{cs}|$ from a global fit in the Standard Model, we obtain ${f_{D^+_s}}=(252.8\pm6.0_{\rm stat}\pm3.7_{\rm syst}\pm0.6_{\rm input})_{\mu\nu}$ MeV and ${f_{D^+_s}}=(259.2\pm3.6_{\rm stat}\pm4.5_{\rm syst}\pm0.6_{\rm input})_{\tau \nu}$ MeV, respectively. Conversely, taking the value for $f_{D_s^+}$ from the latest lattice quantum chromodynamics calculation, we obtain $|V_{cs}| =(0.986\pm0.023_{\rm stat}\pm0.014_{\rm syst}\pm0.003_{\rm input})_{\mu\nu}$ and $|V_{cs}| = (1.011\pm0.014_{\rm stat}\pm0.018_{\rm syst}\pm0.003_{\rm input})_{\tau \nu}$, respectively.

We give a mini-review of scalar field theories with second-derivative Lagrangians, whose field equations are second order. Some of these theories admit solutions violating the Null Energy Condition and having no obvious pathologies. We give a few examples of using these theories in cosmological setting and also in the context of the creation of a universe in the laboratory.

The NA62 experiment at the CERN SPS reports a study of a sample of $4 \times10^{9}$ tagged $\pi^0$ mesons from $K^+ \to \pi^+ \pi^0 (\gamma)$, searching for the decay of the $\pi^0$ to invisible particles. No signal is observed in excess of the expected background fluctuations. An upper limit of $4.4 \times10^{-9}$ is set on the branching ratio at 90% confidence level, improving on previous results by a factor of 60. This result can also be interpreted as a model-independent upper limit on the branching ratio for the decay $K^+ \to \pi^+ X$, where $X$ is a particle escaping detection with mass in the range 0.110-0.155 GeV$/c^2$ and rest lifetime greater than 100 ps. Model-dependent upper limits are obtained assuming $X$ to be an axion-like particle with dominant fermion couplings or a dark scalar mixing with the Standard Model Higgs boson.

We argue that any viable mechanism of gauge field localization should automatically imply charge universality on the brane. We study whether this condition is satisfied in the two known proposals aimed to localize vector field in flat bulk space. We construct a simple calculable model with confinement in the bulk and deconfinement on the brane, as in the Shifman--Dvali set up. We find that in our model the 4-dimensional Coulomb law is indeed reproduced on the brane due to the massless localized photon mode. The charge universality is enforced by the presence of ``confining strings''. On the other hand, charge universality condition is not satisfied in another, brane-induced localization mechanism when the number of extra dimensions d is larger than two. We demonstrate that in the non-Abelian case the gauge fields inside the brane are never four-dimensional and their self-interaction is strong at all distances of interest. Hence this mechanism does not work for d>2. At d=2 the charge universality is still a problem, but it holds automatically at d=1. At d=1, however, the bulk gauge fields are strongly coupled in the non-Abelian case.

Extensive air showers (EAS), produced by cosmic rays in the atmosphere, serve as probes of particle interactions, providing access to energies and kinematical regimes beyond the reach of laboratory experiments. Measurements from multiple cosmic-ray detectors indicate a significant, yet unexplained, discrepancy between the observed muon content in EAS and that predicted by state-of-the-art interaction models, suggesting a need for refinements in our understanding of fundamental physics. Here we show that a tiny, experimentally allowed, violation of the Lorentz invariance (LIV) may result in the suppression of the number of electrons in EAS, leaving the muon number intact and explaining both the ''muon excess'' and its energy dependence. On the other hand, we use the lack of a much stronger discrepancy between EAS data and simulations to obtain strict constraints on the LIV scale. Future experimental tests of this LIV scenario are outlined.

This article gives a concise overview of the development and current status of studies on healthy models of the early Universe without an initial singularity, namely the cosmological bounce and Genesis scenarios, constructed within a broad class of scalar-tensor theories, specifically Horndeski theories and their generalizations. The review focuses on the topics related to linear stability at the perturbation level over the non-singular background solutions: 1) the no-go theorem, valid for non-singular cosmologies within Horndeski theory, 2) the updates on possible approaches to evade the no-go theorem, 3) the role of disformal transformations relating the Horndeski subclasses with the generalized theories like DHOST, 4) the effects on stability caused by additional matter coupling and potential emergence of superluminal perturbation modes in the multi-component setting.

With the establishment and maturation of the experimental programs searching for new physics with sizeable couplings at the LHC, there is an increasing interest in the broader particle and astrophysics community for exploring the physics of light and feebly-interacting particles as a paradigm complementary to a New Physics sector at the TeV scale and beyond. FIPs@LHCb continues the successful series of the FIPs workshops, FIPs 2020 and FIPs 2022. The main focus of the workshop was to explore the LHCb potential to search for FIPs thanks to the new software trigger deployed during the recent upgrade. Equally important goals of the workshop were to update the available parameter space in the commonly used FIPs benchmarks by including recent results from the high energy physics community and to discuss recent theory progress necessary for a more accurate definition of observables related to FIP benchmarks. This document presents the summary of the talks presented at the workshops and the outcome of subsequent discussions.

From gamma-ray burst (GRB) 221009A, very high-energy photons were detected: >10 TeV with LHAASO and >100 TeV with Carpet-3. Such energetic photons are expected to be absorbed via electron-positron pair production on their way to the Earth. Their observation might be explained by new physics, including Lorentz invariance violation (LIV) or photon mixing with axion-like particles (ALPs). Here, we construct a joint fluence spectrum by combining flux measurements from both experiments, and fit it under these hypotheses. LIV can account for the Carpet-3 observation, providing a modest improvement over standard physics in the overall fit. ALP mixing improves the description of both LHAASO and Carpet-3 data, yielding a substantial enhancement in fit quality for a specific region of the ALP parameter space.

The NA62 experiment reports the branching ratio measurement BR$(K^+ \rightarrow \pi^+ \nu\bar{\nu}) = (10.6^{+4.0}_{-3.4} |_{\rm stat} \pm 0.9_{\rm syst}) \times 10 ^{-11}$ at 68% CL, based on the observation of 20 signal candidates with an expected background of 7.0 events from the total data sample collected at the CERN SPS during 2016-2018. This provides evidence for the very rare $K^+ \rightarrow \pi^+ \nu\bar{\nu}$ decay, observed with a significance of 3.4$\sigma$. The experiment achieves a single event sensitivity of $(0.839\pm 0.054)\times 10^{-11}$, corresponding to 10.0 events assuming the Standard Model branching ratio of $(8.4\pm1.0)\times10^{-11}$. This measurement is also used to set limits on BR($K^+ \to \pi^+ X$), where $X$ is a scalar or pseudo-scalar particle. Details are given of the analysis of the 2018 data sample, which corresponds to about 80% of the total data sample.

These notes are an introduction to General Relativity as a Quantum Effective Field Theory, following the material given in a short course on the subject at EPFL. The intent is to develop General Relativity starting from a quantum field theoretic viewpoint, and to introduce some of the techniques needed to understand the subject.

We study a "classical" bouncing scenario in beyond Horndeski theory. We give an example of spatially flat bouncing solution that is non-singular and stable throughout the whole evolution. The model is arranged in such a way that the scalar field driving the cosmological evolution initially behaves like full-fledged beyond Horndeski, whereas at late times it becomes a massless scalar field minimally coupled to gravity.