The issue of inheriting periodicity of an exact solution of a dynamic system by a difference scheme is considered. It is shown that some difference schemes (midpoint scheme, Kahan scheme) in some special cases provide approximate solutions of differential equations, which are periodic sequences. Such solutions are called periodic. A purely algebraic method for finding such solutions is developed. It is shown that midpoint scheme inherits periodicity not only in case of linear oscillator, but also in case of nonlinear oscillator, integrable into elliptic functions.

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.

We provide a coarse classification of all 8-dimensional Manin triples, that describe Poisson--Lie T-dualities between 4-dimensional group manifold solutions to supergravity equations. We find several such dualities and one Poisson--Lie triality.

This work derives exact solutions to the problem of interacting particle density evolution in relativistic and quasi-relativistic approximations for electromagnetic and gravitational interactions. Two types of radial symmetry for the initial density distribution are considered: spherical and cylindrical. It is shown that the relativistic effect delays the onset of the shock wave moment, and in some cases removes it entirely or, conversely, can facilitate it. The analysis of the system's dynamics is carried out within the Wigner-Vlasov formalism, which makes it possible to extend the obtained solutions to quantum systems, including those with gravitational interaction. The derived exact solutions can be directly used as a cross-check for modeling and optimizing nonlinear problems of beam dynamics with account for space charge, astrophysics, plasma physics, and quantum systems with a shock wave.

We study a five-leg scattering amplitude on the special Coulomb branch of planar N=4 super Yang-Mills theory. We reach this point of the moduli space of scalar vacuum expectation values by considering six-dimensional N=(1,1) super Yang-Mills theory and reducing it down to four space-time dimensions with extra-dimensional momenta being nonvanishing. This branch is characterized by massive external W-bosons and massless internal gluons propagating in loops. We analyze the five W-boson amplitude in the kinematics when their masses are much smaller than all Mandelstam-like invariants. This is what we dub the near mass-shell limit. We perform calculations to two-loop order in 't Hooft coupling, making use of recent advances in analytic calculations of required Feynman integrals. Our findings confirm exponentiation of infrared logarithms and enable us to conjecture a concise all-order expression for the amplitude in question. We further analyze its duality to the `square root' of a five-point correlation function of infinitely-heavy half-BPS operators, known as the decagon. By considering the near-null limit for inter-operators distances, we verify that the two objects coincide. This observation corroborates the novel Coulomb amplitudes/heavy correlator duality previously observed for four W-boson amplitudes and Sudakov form factors.

We apply the method of regions to the evaluation of dual conformal integrals with small off-shellness. In contrast to conventional approach, where the separation of regions is performed via dimensional regularization breaking the dual conformal invariance (DCI), we use a sufficiently generic combination of dimensional and analytic regularizations which preserves the DCI. Within this regularization (dubbed as DCI regularization), the contribution of each region becomes DCI. We show that our method dramatically simplifies the calculations. As a demonstration, we calculate the slightly off-shell DCI pentabox integral up to power corrections. The contributions of all 32 regions appear to be expressible in terms of products/ratios of $\Gamma$-functions multiplied by some powers of DCI cross-ratios. Therefore, after removing the regularization, we obtain the final expression in terms of cross-ratios logarithms only. We have checked that our result for pentabox integral numerically agrees with the result of the recent Belitsky\&Smirnov paper [arXiv:2508.14298] which has essentially more complicated form.

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.

Based on the $(2712.4\pm14.4)\times 10^{6}$ $\psi(3686)$ events collected with the BESIII detector, we present a high-precision study of the $\pi^+\pi^-$ mass spectrum in $\psi(3686)\rightarrow\pi^{+}\pi^{-}J/\psi$ decays. A clear resonance-like structure is observed near the $\pi^+\pi^-$ mass threshold for the first time. A fit with a Breit-Wigner function yields a mass of $285.6\pm 2.5~{\rm MeV}/c^2$ and a width of $16.3\pm 0.9~{\rm MeV}$ with a statistical significance exceeding 10$\sigma$. To interpret the data, we incorporate final-state interactions (FSI) within two theoretical frameworks: chiral perturbation theory (ChPT) and QCD multipole expansion (QCDME). ChPT describes the spectrum above 0.3 GeV/$c^2$ but fails to reproduce the threshold enhancement. In contrast, the QCDME model, assuming the $\psi(3686)$ is an admixture of S- and D-wave charmonium, reproduces the data well. The pronounced dip near 0.3 GeV/$c^2$ offers new insight into the interplay between chiral dynamics and low-energy QCD.

Context: The detection of the highest energy neutrino observed to date by KM3NeT, with an estimated energy of 220 PeV, opens up new possibilities for the study and identification of the astrophysical sources responsible for a diffuse flux of such ultra-high-energy neutrinos, among which gamma-ray bursts are longstanding candidates. Aims: Based on the event KM3-230213A, we derive constraints on the baryon loading and density of the surrounding environment in models of blastwaves in long-duration gamma-ray bursts. Methods: We compute the diffuse flux from gamma-ray burst blastwaves, either expanding in a constant density interstellar medium or developing in a radially decreasing density of a wind-like environment surrounding the gamma-ray burst progenitor star, by taking into account the expected neutrino spectra and luminosity function. We use a Poisson likelihood method to constrain the blastwave model parameters by calculating the expected number of neutrino events within the 90% confidence level energy range of KM3-230213A and by using the joint exposure of KM3NeT/ARCA, IceCube and Pierre Auger. Results: We constrain the baryon loading to be $\leq \{392, 131, 39, 13\}$ at 90% confidence level, which is inversely proportional to a varying interstellar medium particle density of $\{1, 3, 10, 30\}$ cm$^{-3}$. In the wind-like environment case, the baryon loading is $\leq \{20, 50, 100\}$ at 90% confidence level, which is proportional to the sixth power of a varying density parameter of $\{0.05, 0.06, 0.07\}$.

Turbulence is indispensable to redistribute nutrients for all life forms larger than microbial, on land and in the ocean. Yet, the development of deep-sea turbulence has not been studied in three dimensions (3D). As a disproportionate laboratory, an array of nearly 3000 high-resolution temperature sensors had been installed for three years on the flat 2500-m deep bottom of the Mediterranean Sea. The time series from the half-cubic hectometer 3D mooring-array allows for the creation of unique movies of deep-sea water motions. Although temperature differences are typically 0.001degrC, variable convection-turbulence is observed as expected from geothermal heating through the flat seafloor. During about 40% of the time, an additional turbulence, 3 times stronger in magnitude, is observed from slantwise advected warmer waters to pass in turbulent clouds. Besides turbulent clouds and seafloor heating, movies also reveal weakly turbulent interfacial-wave breakdown that commonly occurs in the open ocean far away from boundaries.

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 study the cross-section of heavy Higgs production at the LHC within the framework of the Constrained MSSM. It is not only enhanced by tan^2 beta but sometimes is also enhanced by the squark contribution. First, we consider the universal scenario within mSUGRA and find out that to get the desired enhancement one needs large negative values of A_0, which seems to be incompatible with the b->s gamma decay rate. To improve the situation, we release the unification requirement in the Higgs sector. Then it becomes possible to satisfy all requirements simultaneously and enhance the squark contribution. The latter can gain a factor of several units increasing the overall cross-section which, however, is still smaller than the cross-section of the associated H b bbar production. We consider also some other consequences of the chosen benchmark point.

The PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) has measured electrons from heavy flavor (charm and bottom) decays for 0.3 < p_T < 9 GeV/c at midrapidity (|y| < 0.35) in Au+Au collisions at sqrt(s_NN) = 200 GeV. The nuclear modification factor R_AA relative to p+p collisions shows a strong suppression in central Au+Au collisions, indicating substantial energy loss of heavy quarks in the medium produced at RHIC. A large azimuthal anisotropy, v_2, with respect to the reaction plane is observed for 0.5 < p_T < 5 GeV/c indicating non-zero heavy flavor elliptic flow. Both R_AA and v_2 show a p_T dependence different from those of neutral pions. A comparison to transport models which simultaneously describe R_AA(p_T) and v_2(p_T) suggests that the viscosity to entropy density ratio is close to the conjectured quantum lower bound, i.e., near a perfect fluid.

In the paper, the family of conformal four-point ladder diagrams in arbitrary space-time dimensions is considered. We use the representation obtained via explicit calculation using the operator approach and conformal quantum mechanics to study their properties, such as symmetries, loop and dimensional shift identities. In even integer dimensions, latter allows one to reduce the problem to two-dimensional case, where the notable factorization holds. Additionally, for a specific choice of propagator powers, we show that the representation can be written in the form of linear combinations of classical polylogarithms (with coefficients that are rational functions) and explore the structure of the resulting expressions.

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.

We study diffraction of twisted matter waves (electrons and light ions carrying orbital angular momentum $\ell/\hbar=0,\pm1,\pm2,\ldots$ by circular and triangular apertures. Within the scalar Kirchhoff-Fresnel framework, circular apertures preserve cylindrical symmetry and produce ringlike far-field profiles whose radii and widths depend on $|\ell|$ but are insensitive to its sign. In contrast, equilateral triangles break axial symmetry and yield structured patterns that encode both the magnitude and the sign of $\ell$. A transparent Fraunhofer mapping links detector coordinates to the Fourier plane, explaining the $(|\ell|+1)$-lobe rule and the sign-dependent rotation of the pattern. We validate these results for both ideal Bessel beams and localized Laguerre-Gaussian packets, and we cross-check them by split-step Fourier propagation of the time-dependent Schr"odinger equation. From these analyses we extract practical design rules (Fraunhofer distance, lattice pitch, detector sampling) relevant to OAM diagnostics with moderately relativistic electrons with $E_{\rm kin}\sim0.1$ to $5$ MeV and light ions with $E_{\rm kin}\sim0.1$ to $1$ MeV/u. Our results establish triangular diffraction as a simple, passive, and robust method for reading out the OAM content of structured quantum beams.

Parametric integration with hyperlogarithms so far has been successfully used in problems of high energy physics (HEP) and critical statics. In this work, for the first time, it is applied to a problem of critical dynamics, namely, a stochastic model of developed turbulence in high-dimensional spaces, which has a propagator that is non-standard with respect to the HEP: $(-i \omega + \nu k^2)^{-1}$. Adaptation of the hyperlogarithm method is carried out by choosing a proper renormalization scheme and considering an effective dimension of the space. Analytical calculation of the renormalization group functions is performed up to the fourth order of the perturbation theory, $\varepsilon$-expansion of the critical exponent $\omega$ responsible for the infrared stability of the fixed point is obtained.

In this note we address the discrepancy found by Hung, Myers and Smolkin between the holographic calculation of entanglement entropy (using the Jacobson-Myers functional for the holographic minimal surface) and the CFT trace anomaly calculation if one uses the Wald prescription to compute the entropy in six dimensions. As anticipated in our previous work [1] the discrepancy originates entirely from a total derivative term present in the trace anomaly in six dimensions.

The enclosed data release consists of a subset of the calibration data from the Majorana Demonstrator experiment. Each Majorana event is accompanied by raw Germanium detector waveforms, pulse shape discrimination cuts, and calibrated final energies, all shared in an HDF5 file format along with relevant metadata. This release is specifically designed to support the training and testing of Artificial Intelligence (AI) and Machine Learning (ML) algorithms upon our data. This document is structured as follows. Section I provides an overview of the dataset's content and format; Section II outlines the location of this dataset and the method for accessing it; Section III presents the NPML Machine Learning Challenge associated with this dataset; Section IV contains a disclaimer from the Majorana collaboration regarding the use of this dataset; Appendix A contains technical details of this data release. Please direct questions about the material provided within this release to liaobo77@ucsd.edu (A. Li).