Universal Probability Distribution for the Wave Function of a Quantum System Entangled with its Environment

A quantum system (with Hilbert space \({\mathcal {H}_{1}}\)) entangled with its environment (with Hilbert space \({\mathcal {H}_{2}}\)) is usually not attributed to a wave function but only to a reduced density matrix \({\rho_{1}}\). Nevertheless, there is a precise way of attributing to it a random wave function \({\psi_{1}}\), called its conditional wave function, whose probability distribution \({\mu_{1}}\) depends on the entangled wave function \({\psi \in \mathcal {H}_{1} \otimes \mathcal {H}_{2}}\) in the Hilbert space of system and environment together. It also depends on a choice of orthonormal basis of \({\mathcal {H}_{2}}\) but in relevant cases, as we show, not very much. We prove several universality (or typicality) results about \({\mu_{1}}\), e.g., that if the environment is sufficiently large then for every orthonormal basis of \({\mathcal {H}_{2}}\), most entangled states \({\psi}\) with given reduced density matrix \({\rho_{1}}\) are such that \({\mu_{1}}\) is close to one of the so-called GAP (Gaussian adjusted projected) measures, \({GAP(\rho_{1})}\). We also show that, for most entangled states \({\psi}\) from a microcanonical subspace (spanned by the eigenvectors of the Hamiltonian with energies in a narrow interval \({[E, E+ \delta E]}\)) and most orthonormal bases of \({\mathcal {H}_{2}}\), \({\mu_{1}}\) is close to \({GAP(\rm {tr}_{2} \rho_{mc})}\) with \({\rho_{mc}}\) the normalized projection to the microcanonical subspace. In particular, if the coupling between the system and the environment is weak, then \({\mu_{1}}\) is close to \({GAP(\rho_\beta)}\) with \({\rho_\beta}\) the canonical density matrix on \({\mathcal {H}_{1}}\) at inverse temperature \({\beta=\beta(E)}\). This provides the mathematical justification of our claim in Goldstein et al. (J Stat Phys 125: 1193–1221, 2006) that GAP measures describe the thermal equilibrium distribution of the wave function.     

Light deflection and Gauss–Bonnet theorem: definition of total deflection angle and its applications

In this paper, we re-examine the light deflection in the Schwarzschild and the Schwarzschild–de Sitter spacetime. First, supposing a static and spherically symmetric spacetime, we propose the definition of the total deflection angle \(\alpha \) of the light ray by constructing a quadrilateral \(\varSigma ^4\) on the optical reference geometry \({\mathscr {M}}^\mathrm{opt}\) determined by the optical metric \(\bar{g}_{ij}\). On the basis of the definition of the total deflection angle \(\alpha \) and the Gauss–Bonnet theorem, we derive two formulas to calculate the total deflection angle \(\alpha \); (1) the angular formula that uses four angles determined on the optical reference geometry \({\mathscr {M}}^\mathrm{opt}\) or the curved \((r, \phi )\) subspace \({\mathscr {M}}^\mathrm{sub}\) being a slice of constant time t and (2) the integral formula on the optical reference geometry \({\mathscr {M}}^\mathrm{opt}\) which is the areal integral of the Gaussian curvature K in the area of a quadrilateral \(\varSigma ^4\) and the line integral of the geodesic curvature \(\kappa _g\) along the curve \(C_{\varGamma }\). As the curve \(C_{\varGamma }\), we introduce the unperturbed reference line that is the null geodesic \(\varGamma \) on the background spacetime such as the Minkowski or the de Sitter spacetime, and is obtained by projecting \(\varGamma \) vertically onto the curved \((r, \phi )\) subspace \({\mathscr {M}}^\mathrm{sub}\). We demonstrate that the two formulas give the same total deflection angle \(\alpha \) for the Schwarzschild and the Schwarzschild–de Sitter spacetime. In particular, in the Schwarzschild case, the result coincides with Epstein–Shapiro’s formula when the source S and the receiver R of the light ray are located at infinity. In addition, in the Schwarzschild–de Sitter case, there appear order \({\mathscr {O}}(\varLambda m)\) terms in addition to the Schwarzschild-like part, while order \({\mathscr {O}}(\varLambda )\) terms disappear.     

On the Convexity of the KdV Hamiltonian

Motivated by perturbation theory, we prove that the nonlinear part \({H^{*}}\) of the KdV Hamiltonian \({H^{kdv}}\), when expressed in action variables \({I = (I_{n})_{n \geqslant 1}}\), extends to a real analytic function on the positive quadrant \({\ell^{2}_{+}(\mathbb{N})}\) of \({\ell^{2}(\mathbb{N})}\) and is strictly concave near \({0}\). As a consequence, the differential of \({H^{*}}\) defines a local diffeomorphism near 0 of \({\ell_{\mathbb{C}}^{2}(\mathbb{N})}\). Furthermore, we prove that the Fourier-Lebesgue spaces \({\mathcal{F}\mathcal{L}^{s,p}}\) with \({-1/2 \leqslant s \leqslant 0}\) and \({2 \leqslant p < \infty}\), admit global KdV-Birkhoff coordinates. In particular, it means that \({\ell^{2}_+(\mathbb{N})}\) is the space of action variables of the underlying phase space \({\mathcal{F}\mathcal{L}^{-1/2,4}}\) and that the KdV equation is globally in time \({C^{0}}\)-well-posed on \({\mathcal{F}\mathcal{L}^{-1/2,4}}\).     

The Structure of the Kac–Wang–Yan Algebra

The Lie algebra \({\mathcal{D}}\) of regular differential operators on the circle has a universal central extension \({\hat{\mathcal{D}}}\). The invariant subalgebra \({\hat{\mathcal{D}}^+}\) under an involution preserving the principal gradation was introduced by Kac, Wang, and Yan. The vacuum \({\hat{\mathcal{D}}^+}\)-module with central charge \({c \in \mathbb{C}}\), and its irreducible quotient \({\mathcal{V}_c}\), possess vertex algebra structures, and \({\mathcal{V}_c}\) has a nontrivial structure if and only if \({c \in \frac{1}{2}\mathbb{Z}}\). We show that for each integer \({n > 0}\), \({\mathcal{V}_{n/2}}\) and \({\mathcal{V}_{-n}}\) are \({\mathcal{W}}\)-algebras of types \({\mathcal{W}(2, 4,\dots,2n)}\) and \({\mathcal{W}(2, 4,\dots, 2n^2 + 4n)}\), respectively. These results are formal consequences of Weyl’s first and second fundamental theorems of invariant theory for the orthogonal group \({{\rm O}(n)}\) and the symplectic group \({{\rm Sp}(2n)}\), respectively. Based on Sergeev’s theorems on the invariant theory of \({{\rm Osp}(1, 2n)}\) we conjecture that \({\mathcal{V}_{-n+1/2}}\) is of type \({\mathcal{W}(2, 4,\dots, 4n^2 + 8n + 2)}\), and we prove this for \({n = 1}\). As an application, we show that invariant subalgebras of \({\beta\gamma}\)-systems and free fermion algebras under arbitrary reductive group actions are strongly finitely generated.     

The Exoticness and Realisability of Twisted Haagerup–Izumi Modular Data

The quantum double of the Haagerup subfactor, the first irreducible finite depth subfactor with index above 4, is the most obvious candidate for exotic modular data. We show that its modular data \({\mathcal{D}{\rm Hg}}\) fits into a family \({\mathcal{D}^\omega {\rm Hg}_{2n+1}}\) , where n ≥  0 and \({\omega\in \mathbb{Z}_{2n+1}}\) . We show \({\mathcal{D}^0 {\rm Hg}_{2n+1}}\) is related to the subfactors Izumi hypothetically associates to the cyclic groups \({\mathbb{Z}_{2n+1}}\) . Their modular data comes equipped with canonical and dual canonical modular invariants; we compute the corresponding alpha-inductions, etc. In addition, we show there are (respectively) 1, 2, 0 subfactors of Izumi type \({\mathbb{Z}_7, \mathbb{Z}_9}\) and \({\mathbb{Z}_3^2}\) , and find numerical evidence for 2, 1, 1, 1, 2 subfactors of Izumi type \({\mathbb{Z}_{11},\mathbb{Z}_{13},\mathbb{Z}_{15},\mathbb{Z}_{17},\mathbb{Z}_{19}}\) (previously, Izumi had shown uniqueness for \({\mathbb{Z}_3}\) and \({\mathbb{Z}_5}\)), and we identify their modular data. We explain how \({\mathcal{D}{\rm Hg}}\) (more generally \({\mathcal{D}^\omega {\rm Hg}_{2n+1}}\)) is a graft of the quantum double \({\mathcal{D} Sym(3)}\) (resp. the twisted double \({\mathcal{D}^\omega D_{2n+1}}\)) by affine so(13) (resp. so\({(4n^2+4n+5)}\)) at level 2. We discuss the vertex operator algebra (or conformal field theory) realisation of the modular data \({\mathcal{D}^\omega {\rm Hg}_{2n+1}}\) . For example we show there are exactly 2 possible character vectors (giving graded dimensions of all modules) for the Haagerup VOA at central charge c = 8. It seems unlikely that any of this twisted Haagerup-Izumi modular data can be regarded as exotic, in any reasonable sense.     

Time Delay for the Dirac Equation

We consider time delay for the Dirac equation. A new method to calculate the asymptotics of the expectation values of the operator \({\int\limits_{0} ^{\infty}{\rm e}^{iH_{0}t}\zeta(\frac{\vert x\vert }{R}) {\rm e}^{-iH_{0}t}{\rm d}t}\), as \({R \rightarrow \infty}\), is presented. Here, H₀ is the free Dirac operator and \({\zeta\left(t\right)}\) is such that \({\zeta\left(t\right) = 1}\) for \({0 \leq t \leq 1}\) and \({\zeta\left(t\right) = 0}\) for \({t > 1}\). This approach allows us to obtain the time delay operator \({\delta \mathcal{T}\left(f\right)}\) for initial states f in \({\mathcal{H} _{2}^{3/2+\varepsilon}(\mathbb{R}^{3};\mathbb{C}^{4})}\), \({\varepsilon > 0}\), the Sobolev space of order \({3/2+\varepsilon}\) and weight 2. The relation between the time delay operator \({\delta\mathcal{T}\left(f\right)}\) and the Eisenbud–Wigner time delay operator is given. In addition, the relation between the averaged time delay and the spectral shift function is presented.     

Constraining the electric charges of some astronomical bodies in Reissner–Nordström spacetimes and generic <Emphasis Type="Italic">r</Emphasis><Superscript>−2</Superscript>-type power-law potentials from orbital motions

We put independent model dynamical constraints on the net electric charge Q of some astronomical and astrophysical objects by assuming that their exterior spacetimes are described by the Reissner-Nordström, metric, which induces an additional potential \({U_{\rm RN} \propto Q^2 r^{-2}}\). From the current bounds \({\Delta \dot \varpi}\) on any anomalies in the secular perihelion rate \({\dot \varpi}\) of Mercury and the Earth–mercury ranging Δρ, we have \({\left|Q_{\odot}\right| \lesssim 1-0.4 \times 10^{18}\ {\rm C}}\). Such constraints are ~60–200 times tighter than those recently inferred in literature. For the Earth, the perigee precession of the Moon, determined with the Lunar Laser Ranging technique, and the intersatellite ranging Δρ for the GRACE mission yield \({\left|Q_{\oplus} \right| \lesssim 5-0.4 \times 10^{14}\ {\rm C}}\). The periastron rate of the double pulsar PSR J0737-3039A/B system allows to infer \({\left|Q_{\rm NS} \right| \lesssim 5\times 10^{19}\ {\rm C}}\). According to the perinigricon precession of the main sequence S2 star in Sgr A^*, the electric charge carried by the compact object hosted in the Galactic Center may be as large as \({\left|Q_{\bullet} \right| \lesssim 4\times 10^{27} \ {\rm C}}\). Our results extend to other hypothetical power-law interactions inducing extra-potentials \({U_{\rm pert} = \Psi r^{-2}}\) as well. It turns out that the terrestrial GRACE mission yields the tightest constraint on the parameter \({\Psi}\), assumed as a universal constant, amounting to \({|\Psi| \lesssim 5\times 10^{9}\ {\rm m^4\ s^{-2}}}\).     

Spin density matrix of the $$\omega $$ in the reaction $${\bar{p}p}\,\rightarrow \,\omega \pi ^0$$p¯p→π0

The spin density matrix of the \(\omega \) has been determined for the reaction \({\bar{p}p}\,\rightarrow \,\omega \pi ^0\) with unpolarized in-flight data measured by the Crystal Barrel LEAR experiment at CERN. The two main decay modes of the \(\omega \) into \(\pi ^0 \gamma \) and \(\pi ^+ \pi ^- \pi ^0\) have been separately analyzed for various \({\bar{p}}\)momenta between 600 and 1940 MeV/c. The results obtained with the usual method by extracting the matrix elements via the \(\omega \) decay angular distributions and with the more sophisticated method via a full partial wave analysis are in good agreement. A strong spin alignment of the \(\omega \) is clearly visible in this energy regime and all individual spin density matrix elements exhibit an oscillatory dependence on the production angle. In addition, the largest contributing orbital angular momentum of the \({\bar{p}p~}\)system has been identified for the different beam momenta. It increases from \(L^{max}_{{\bar{p}p~}}\) \(=\) 2 at 600 MeV/c to \(L^{max}_{{\bar{p}p~}}\) \(=\) 5 at 1940 MeV/c.     

Stochastic Averaging Principle for Spatial Birth-and-Death Evolutions in the Continuum

We study a spatial birth-and-death process on the phase space of locally finite configurations \({\varGamma }^+ \times {\varGamma }^-\) over \({\mathbb {R}}^d\). Dynamics is described by an non-equilibrium evolution of states obtained from the Fokker-Planck equation and associated with the Markov operator \(L^+(\gamma ^-) + \frac{1}{\varepsilon }L^-\), \(\varepsilon > 0\). Here \(L^-\) describes the environment process on \({\varGamma }^-\) and \(L^+(\gamma ^-)\) describes the system process on \({\varGamma }^+\), where \(\gamma ^-\) indicates that the corresponding birth-and-death rates depend on another locally finite configuration \(\gamma ^- \in {\varGamma }^-\). We prove that, for a certain class of birth-and-death rates, the corresponding Fokker-Planck equation is well-posed, i.e. there exists a unique evolution of states \(\mu _t^{\varepsilon }\) on \({\varGamma }^+ \times {\varGamma }^-\). Moreover, we give a sufficient condition such that the environment is ergodic with exponential rate. Let \(\mu _{\mathrm {inv}}\) be the invariant measure for the environment process on \({\varGamma }^-\). In the main part of this work we establish the stochastic averaging principle, i.e. we prove that the marginal of \(\mu _t^{\varepsilon }\) onto \({\varGamma }^+\) converges weakly to an evolution of states on \({\varGamma }^+\) associated with the averaged Markov birth-and-death operator \({\overline{L}} = \int _{{\varGamma }^-}L^+(\gamma ^-)d \mu _{\mathrm {inv}}(\gamma ^-)\).     

On a Pre-Lie Algebra Defined by Insertion of Rooted Trees

Let K be a field of characteristic zero. For \({n \in \mathbb{N}^{*}}\) , let \({\mathcal{T}^{\,\prime}_{n}}\) be the vector space of non-planar rooted trees with n vertices (Foissy in Bull Sci Math 126, no. 3, 193–239; no. 4, 249–288, 2002). Let \({\vartriangleright}\) be the left pre-Lie product of insertion of a tree inside another defined on infinitesimal characters of the graded Hopf algebra \({\mathcal{H}}\) introduced by Calaque, Ebrahimi-Fard and Manchon. Let \({\mathcal{T}^{\,\prime}=\oplus_{n\geq 2}\mathcal{T}^{\,\prime}_{n}}\) . In this work, we first prove that \({(\mathcal{T}^{\,\prime}, \vartriangleright)}\) a pre-Lie algebra generated by the two ladders E ₁ and E ₂ where E ₁ is the ladder with one edge and E ₂ is the ladder with two edges. Second, we prove that \({(\mathcal{T}^{\,\prime}, \vartriangleright)}\) is not a free pre-Lie algebra, and we exhibit a family of relations.     

Likelihood analysis of the pMSSM11 in light of LHC 13-TeV data

We use MasterCode to perform a frequentist analysis of the constraints on a phenomenological MSSM model with 11 parameters, the pMSSM11, including constraints from \(\sim 36\)/fb of LHC data at 13 TeV and PICO, XENON1T and PandaX-II searches for dark matter scattering, as well as previous accelerator and astrophysical measurements, presenting fits both with and without the \((g-2)_\mu \) constraint. The pMSSM11 is specified by the following parameters: 3 gaugino masses \(M_{1,2,3}\), a common mass for the first-and second-generation squarks \(m_{\tilde{q}}\) and a distinct third-generation squark mass \(m_{\tilde{q}_3}\), a common mass for the first-and second-generation sleptons \(m_{\tilde{\ell }}\) and a distinct third-generation slepton mass \(m_{\tilde{\tau }}\), a common trilinear mixing parameter A, the Higgs mixing parameter \(\mu \), the pseudoscalar Higgs mass \(M_A\) and \(\tan \beta \). In the fit including \((g-2)_\mu \), a Bino-like \(\tilde{\chi }^0_{1}\) is preferred, whereas a Higgsino-like \(\tilde{\chi }^0_{1}\) is mildly favoured when the \((g-2)_\mu \) constraint is dropped. We identify the mechanisms that operate in different regions of the pMSSM11 parameter space to bring the relic density of the lightest neutralino, \(\tilde{\chi }^0_{1}\), into the range indicated by cosmological data. In the fit including \((g-2)_\mu \), coannihilations with \(\tilde{\chi }^0_{2}\) and the Wino-like \(\tilde{\chi }^\pm _{1}\) or with nearly-degenerate first- and second-generation sleptons are active, whereas coannihilations with the \(\tilde{\chi }^0_{2}\) and the Higgsino-like \(\tilde{\chi }^\pm _{1}\) or with first- and second-generation squarks may be important when the \((g-2)_\mu \) constraint is dropped. In the two cases, we present \(\chi ^2\) functions in two-dimensional mass planes as well as their one-dimensional profile projections and best-fit spectra. Prospects remain for discovering strongly-interacting sparticles at the LHC, in both the scenarios with and without the \((g-2)_\mu \) constraint, as well as for discovering electroweakly-interacting sparticles at a future linear \(e^+ e^-\) collider such as the ILC or CLIC.     

SU(3) symmetry breaking in charmed baryon decays

We explore the breaking effects of the SU(3) flavor symmetry in the singly Cabibbo-suppressed anti-triplet charmed baryon decays of \(\mathbf{B}_c\rightarrow \mathbf{B}_n M\), with \(\mathbf{B}_c=(\Xi _c^0,\Xi _c^+,\Lambda _c^+)\) and \(\mathbf{B}_n(M)\) the baryon (pseudo-scalar) octets. We find that these breaking effects can be used to account for the experimental data on the decay branching ratios of \({\mathcal {B}}(\Lambda _c^+\rightarrow \Sigma ^{0} K^{+},\Lambda ^{0} K^{+})\) and \(R'_{K/\pi }={\mathcal {B}}(\Xi ^0_c \rightarrow \Xi ^- K^+)\)/\({\mathcal {B}}(\Xi ^0_c \rightarrow \Xi ^- \pi ^+)\). In addition, we obtain that \({\mathcal {B}}(\Xi _{c}^{0} \rightarrow \Xi ^{-} K^{+},\Sigma ^{-} \pi ^{+})=(4.6 \pm 1.7,12.8 \pm 3.1)\times 10^{-4}\), \({\mathcal {B}}(\Xi _c^0\rightarrow pK^-,\Sigma ^+\pi ^-)=(3.0 \pm 1.0, 5.2 \pm 1.6)\times 10^{-4}\) and \({\mathcal {B}}(\Xi _c^+\rightarrow \Sigma ^{0(+)} \pi ^{+(0)})=(10.3 \pm 1.7)\times 10^{-4}\), which all receive significant contributions from the breaking effects, and can be tested by the BESIII and LHCb experiments.     

Noncommutative products of Euclidean spaces

We present natural families of coordinate algebras on noncommutative products of Euclidean spaces \({\mathbb {R}}^{N_1} \times _{\mathcal {R}} {\mathbb {R}}^{N_2}\). These coordinate algebras are quadratic ones associated with an \(\mathcal {R}\)-matrix which is involutive and satisfies the Yang–Baxter equations. As a consequence, they enjoy a list of nice properties, being regular of finite global dimension. Notably, we have eight-dimensional noncommutative euclidean spaces \({\mathbb {R}}^{4} \times _{\mathcal {R}} {\mathbb {R}}^{4}\). Among these, particularly well behaved ones have deformation parameter \(\mathbf{u} \in {\mathbb {S}}^2\). Quotients include seven spheres \({\mathbb {S}}^{7}_\mathbf{u}\) as well as noncommutative quaternionic tori \({\mathbb {T}}^{{\mathbb {H}}}_\mathbf{u} = {\mathbb {S}}^3 \times _\mathbf{u} {\mathbb {S}}^3\). There is invariance for an action of \({{\mathrm{SU}}}(2) \times {{\mathrm{SU}}}(2)\) on the torus \({\mathbb {T}}^{{\mathbb {H}}}_\mathbf{u}\) in parallel with the action of \(\mathrm{U}(1) \times \mathrm{U}(1)\) on a ‘complex’ noncommutative torus \({\mathbb {T}}^2_\theta \) which allows one to construct quaternionic toric noncommutative manifolds. Additional classes of solutions are disjoint from the classical case.     

Harmonic Pinnacles in the Discrete Gaussian Model

The 2D Discrete Gaussian model gives each height function \({\eta : {\mathbb{Z}^2\to\mathbb{Z}}}\) a probability proportional to \({\exp(-\beta \mathcal{H}(\eta))}\), where \({\beta}\) is the inverse-temperature and \({\mathcal{H}(\eta) = \sum_{x\sim y}(\eta_x-\eta_y)^2}\) sums over nearest-neighbor bonds. We consider the model at large fixed \({\beta}\), where it is flat unlike its continuous analog (the Discrete Gaussian Free Field). We first establish that the maximum height in an \({L\times L}\) box with 0 boundary conditions concentrates on two integers M, M + 1 with \({M\sim \sqrt{(1/2\pi\beta)\log L\log\log L}}\). The key is a large deviation estimate for the height at the origin in \({\mathbb{Z}^{2}}\), dominated by “harmonic pinnacles”, integer approximations of a harmonic variational problem. Second, in this model conditioned on \({\eta\geq 0}\) (a floor), the average height rises, and in fact the height of almost all sites concentrates on levels H, H + 1 where \({H\sim M/\sqrt{2}}\). This in particular pins down the asymptotics, and corrects the order, in results of Bricmont et al. (J. Stat. Phys. 42(5–6):743–798, 1986), where it was argued that the maximum and the height of the surface above a floor are both of order \({\sqrt{\log L}}\). Finally, our methods extend to other classical surface models (e.g., restricted SOS), featuring connections to p-harmonic analysis and alternating sign matrices.     

Quantum Electrodynamics of Atomic Resonances

A simple model of an atom interacting with the quantized electromagnetic field is studied. The atom has a finite mass m, finitely many excited states and an electric dipole moment, \({\vec{d}_0 = -\lambda_{0} \vec{d}}\), where \({\| d^{i}\| = 1, i = 1, 2, 3,}\) and \({\lambda_0}\) is proportional to the elementary electric charge. The interaction of the atom with the radiation field is described with the help of the Ritz Hamiltonian, \({-\vec{d}_0 \cdot \vec{E}}\), where \({\vec{E}}\) is the electric field, cut off at large frequencies. A mathematical study of the Lamb shift, the decay channels and the life times of the excited states of the atom is presented. It is rigorously proven that these quantities are analytic functions of the momentum \({\vec{p}}\) of the atom and of the coupling constant \({\lambda_0}\), provided \({\vert\vec{p} \vert < mc}\) and \({\vert \Im \vec{p} \vert}\) and \({\vert \lambda_{0} \vert}\) are sufficiently small. The proof relies on a somewhat novel inductive construction involving a sequence of ‘smooth Feshbach–Schur maps’ applied to a complex dilatation of the original Hamiltonian, which yields an algorithm for the calculation of resonance energies that converges super-exponentially fast.     

Scaling Limit of Symmetric Random Walk in High-Contrast Periodic Environment

It is shown that the deterministic infinite trigonometric products

$$\begin{aligned} \prod _{n\in \mathbb {N}}\left[ 1- p +p\cos \left( \textstyle n^{-s}_{_{}}t\right) \right] =: {\text{ Cl }_{p;s}^{}}(t) \end{aligned}$$

with parameters \( p\in (0,1]\ \& \ s>\frac{1}{2}\), and variable \(t\in \mathbb {R}\), are inverse Fourier transforms of the probability distributions for certain random series \(\Omega _{p}^\zeta (s)\) taking values in the real \(\omega \) line; i.e. the \({\text{ Cl }_{p;s}^{}}(t)\) are characteristic functions of the \(\Omega _{p}^\zeta (s)\). The special case \(p=1=s\) yields the familiar random harmonic series, while in general \(\Omega _{p}^\zeta (s)\) is a “random Riemann-\(\zeta \) function,” a notion which will be explained and illustrated—and connected to the Riemann hypothesis. It will be shown that \(\Omega _{p}^\zeta (s)\) is a very regular random variable, having a probability density function (PDF) on the \(\omega \) line which is a Schwartz function. More precisely, an elementary proof is given that there exists some \(K_{p;s}^{}>0\), and a function \(F_{p;s}^{}(|t|)\) bounded by \(|F_{p;s}^{}(|t|)|\!\le \! \exp \big (K_{p;s}^{} |t|^{1/(s+1)})\), and \(C_{p;s}^{}\!:=\!-\frac{1}{s}\int _0^\infty \ln |{1-p+p\cos \xi }|\frac{1}{\xi ^{1+1/s}}\mathrm{{d}}\xi \), such that

$$\begin{aligned} \forall \,t\in \mathbb {R}:\quad {\text{ Cl }_{p;s}^{}}(t) = \exp \bigl ({- C_{p;s}^{} \,|t|^{1/s}\bigr )F_{p;s}^{}(|t|)}; \end{aligned}$$

the regularity of \(\Omega _{p}^\zeta (s)\) follows. Incidentally, this theorem confirms a surmise by Benoit Cloitre, that \(\ln {\text{ Cl }_{{{1}/{3}};2}^{}}(t) \sim -C\sqrt{t}\; \left( t\rightarrow \infty \right) \) for some \(C>0\). Graphical evidence suggests that \({\text{ Cl }_{{{1}/{3}};2}^{}}(t)\) is an empirically unpredictable (chaotic) function of t. This is reflected in the rich structure of the pertinent PDF (the Fourier transform of \({\text{ Cl }_{{{1}/{3}};2}^{}}\)), and illustrated by random sampling of the Riemann-\(\zeta \) walks, whose branching rules allow the build-up of fractal-like structures.

     

Singularity of the Velocity Distribution Function in Molecular Velocity Space

We study the boundary singularity of the solutions to the Boltzmann equation in the kinetic theory. The solution has a jump discontinuity in the microscopic velocity \({\zeta}\) on the boundary and a secondary singularity of logarithmic type around the velocity tangential to the boundary, \({\zeta_{n} \sim 0_{-}}\), where \({\zeta_{n}}\) is the component of molecular velocity normal to the boundary, pointing to the gas. We demonstrate this secondary singularity by obtaining an asymptotic formula for the derivative of the solution on the boundary with respect to \({\zeta_{n}}\) that diverges logarithmically when \({\zeta_{n} \sim 0_{-}}\). Our study is for the thermal transpiration problem between two plates for the hard sphere gases with sufficiently large Knudsen number and with the diffuse reflection boundary condition. The solution is constructed and its singularity is studied by an iteration procedure.     

Likelihood analysis of the minimal AMSB model

We perform a likelihood analysis of the minimal anomaly-mediated supersymmetry-breaking (mAMSB) model using constraints from cosmology and accelerator experiments. We find that either a wino-like or a Higgsino-like neutralino LSP, \(\tilde{\chi }^0_{1}\), may provide the cold dark matter (DM), both with similar likelihoods. The upper limit on the DM density from Planck and other experiments enforces \(m_{\tilde{\chi }^0_{1}} \lesssim 3 \,\, \mathrm {TeV}\) after the inclusion of Sommerfeld enhancement in its annihilations. If most of the cold DM density is provided by the \(\tilde{\chi }^0_{1}\), the measured value of the Higgs mass favours a limited range of \(\tan \beta \sim 5\) (and also for \(\tan \beta \sim 45\) if \(\mu > 0\)) but the scalar mass \(m_0\) is poorly constrained. In the wino-LSP case, \(m_{3/2}\) is constrained to about \(900\,\, \mathrm {TeV}\) and \(m_{\tilde{\chi }^0_{1}}\) to \(2.9\pm 0.1\,\, \mathrm {TeV}\), whereas in the Higgsino-LSP case \(m_{3/2}\) has just a lower limit \(\gtrsim 650\,\, \mathrm {TeV}\) (\(\gtrsim 480\,\, \mathrm {TeV}\)) and \(m_{\tilde{\chi }^0_{1}}\) is constrained to \(1.12 ~(1.13) \pm 0.02\,\, \mathrm {TeV}\) in the \(\mu >0\) (\(\mu <0\)) scenario. In neither case can the anomalous magnetic moment of the muon, \((g-2)_\mu \), be improved significantly relative to its Standard Model (SM) value, nor do flavour measurements constrain the model significantly, and there are poor prospects for discovering supersymmetric particles at the LHC, though there are some prospects for direct DM detection. On the other hand, if the \(\tilde{\chi }^0_{1}\) contributes only a fraction of the cold DM density, future LHC Open image in new window -based searches for gluinos, squarks and heavier chargino and neutralino states as well as disappearing track searches in the wino-like LSP region will be relevant, and interference effects enable \(\mathrm{BR}(B_{s, d} \rightarrow \mu ^+\mu ^-)\) to agree with the data better than in the SM in the case of wino-like DM with \(\mu > 0\).