Systems of conservation laws with third-order Hamiltonian structures

We investigate n-component systems of conservation laws that possess third-order Hamiltonian structures of differential-geometric type. The classification of such systems is reduced to the projective classification of linear congruences of lines in \(\mathbb {P}^{n+2}\) satisfying additional geometric constraints. Algebraically, the problem can be reformulated as follows: for a vector space W of dimension \(n+2\), classify n-tuples of skew-symmetric 2-forms \(A^{\alpha } \in \varLambda ^2(W)\) such that

$$\begin{aligned} \phi _{\beta \gamma }A^{\beta }\wedge A^{\gamma }=0, \end{aligned}$$

for some non-degenerate symmetric \(\phi \).

     

A Theory of Hydrodynamic Turbulence Based on Non-equilibrium Statistical Mechanics

In earlier papers, we have studied the turbulent flow exponents \(\zeta _p\), where \(\langle |\Delta \mathbf{v}|^p\rangle \sim \ell ^{\zeta _p}\) and \(\Delta \mathbf{v}\) is the contribution to the fluid velocity at small scale \(\ell \). Using ideas of non-equilibrium statistical mechanics we have found

$$\begin{aligned} \zeta _p={p\over 3}-{1\over \ln \kappa }\ln \Gamma \left( {p\over 3}+1\right) \end{aligned}$$

where \(1/\ln \kappa \) is experimentally \(\approx \,0.32\,\pm \,0.01\). The purpose of the present note is to propose a somewhat more physical derivation of the formula for \(\zeta _p\). We also present an estimate \(\approx \,100\) for the Reynolds number at the onset of turbulence.

     

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.

     

The Planar Ising Model and Total Positivity

A matrix is called totally positive (resp. totally nonnegative) if all its minors are positive (resp. nonnegative). Consider the Ising model with free boundary conditions and no external field on a planar graph G. Let \(a_1,\dots ,a_k,b_k,\dots ,b_1\) be vertices placed in a counterclockwise order on the outer face of G. We show that the \(k\times k\) matrix of the two-point spin correlation functions

$$\begin{aligned} M_{i,j} = \langle \sigma _{a_i} \sigma _{b_j} \rangle \end{aligned}$$

is totally nonnegative. Moreover, \(\det M > 0\) if and only if there exist k pairwise vertex-disjoint paths that connect \(a_i\) with \(b_i\). We also compute the scaling limit at criticality of the probability that there are k parallel and disjoint connections between \(a_i\) and \(b_i\) in the double random current model. Our results are based on a new distributional relation between double random currents and random alternating flows of Talaska [37].

     

Spherical Model on a Cayley Tree: Large Deviations

We study the spherical model of a ferromagnet on a Cayley tree and show that in the case of empty boundary conditions a ferromagnetic phase transition takes place at the critical temperature \(T_\mathrm{c} =\frac{6\sqrt{2}}{5}J\), where J is the interaction strength. For any temperature the equilibrium magnetization, \(m_n\), tends to zero in the thermodynamic limit, and the true order parameter is the renormalized magnetization \(r_n=n^{3/2}m_n\), where n is the number of generations in the Cayley tree. Below \(T_\mathrm{c}\), the equilibrium values of the order parameter are given by \(\pm \rho ^*\), where

$$\begin{aligned} \rho ^*=\frac{2\pi }{(\sqrt{2}-1)^2}\sqrt{1-\frac{T}{T_\mathrm{c}}}. \end{aligned}$$

One more notable temperature in the model is the penetration temperature

$$\begin{aligned} T_\mathrm{p}=\frac{J}{W_\mathrm{Cayley}(3/2)}\left( 1-\frac{1}{\sqrt{2}}\left( \frac{h}{2J}\right) ^2\right) . \end{aligned}$$

Below \(T_\mathrm{p}\) the influence of homogeneous boundary field of magnitude h penetrates throughout the tree. The main new technical result of the paper is a complete set of orthonormal eigenvectors for the discrete Laplace operator on a Cayley tree.

     

Free Energy of the Cauchy Directed Polymer Model at High Temperature

We study the Cauchy directed polymer model on \(\mathbb {Z}^{1+1}\), where the underlying random walk is in the domain of attraction to the 1-stable law. We show that, if the random walk satisfies certain regularity assumptions and its symmetrized version is recurrent, then the free energy is strictly negative at any inverse temperature \(\beta >0\). Moreover, under additional regularity assumptions on the random walk, we can identify the sharp asymptotics of the free energy in the high temperature limit, namely,

$$\begin{aligned} \lim \limits _{\beta \rightarrow 0}\beta ^{2}\log (-p(\beta ))=-c. \end{aligned}$$

     

A Mysterious Cluster Expansion Associated to the Expectation Value of the Permanent of 0–1 Matrices

We consider two ensembles of \(0-1\) \(n\times n\) matrices. The first is the set of all \(n\times n\) matrices with entries zeroes and ones such that all column sums and all row sums equal r, uniformly weighted. The second is the set of \(n \times n\) matrices with zero and one entries where the probability that any given entry is one is r / n, the probabilities of the set of individual entries being i.i.d.’s. Calling the two expectation values E and \(E_B\) respectively, we develop a formal relation

$$\begin{aligned} E({{\mathrm{perm}}}(A)) = E_B({{\mathrm{perm}}}(A)) e^{\sum _2 T_i}.\quad \quad \quad \quad \mathrm{(A1)} \end{aligned}$$

We use two well-known approximating ensembles to E, \(E_1\) and \(E_2\). Replacing E by either \(E_1\) or \(E_2\) we can evaluate all terms in (A1). For either \(E_1\) or \(E_2\) the terms \(T_i\) have amazing properties. We conjecture that all these properties hold also for E. We carry through a similar development treating \(E({{\mathrm{perm}}}_m(A))\), with m proportional to n, in place of \(E({{\mathrm{perm}}}(A))\).

     

A Lax pair structure for the half-wave maps equation

We consider the half-wave maps equation

$$\begin{aligned} \partial _t \vec {S} = \vec {S} \wedge |\nabla | \vec {S}, \end{aligned}$$

where \(\vec {S}= \vec {S}(t,x)\) takes values on the two-dimensional unit sphere \(\mathbb {S}^2\) and \(x \in \mathbb {R}\) (real line case) or \(x \in \mathbb {T}\) (periodic case). This an energy-critical Hamiltonian evolution equation recently introduced in Lenzmann and Schikorra (2017, arXiv:1702.05995v2), Zhou and Stone (Phys Lett A 379:2817–2825, 2015) which formally arises as an effective evolution equation in the classical and continuum limit of Haldane–Shastry quantum spin chains. We prove that the half-wave maps equation admits a Lax pair and we discuss some analytic consequences of this finding. As a variant of our arguments, we also obtain a Lax pair for the half-wave maps equation with target \(\mathbb {H}^2\) (hyperbolic plane).

     

The Largest Fragment of a Homogeneous Fragmentation Process

We show that in homogeneous fragmentation processes the largest fragment at time t has size

$$\begin{aligned} e^{-t \Phi '(\overline{p})}t^{-\frac{3}{2} (\log \Phi )'(\overline{p})+o(1)}, \end{aligned}$$

where \(\Phi \) is the Lévy exponent of the fragmentation process, and \(\overline{p}\) is the unique solution of the equation \((\log \Phi )'(\bar{p})=\frac{1}{1+\bar{p}}\). We argue that this result is in line with predictions arising from the classification of homogeneous fragmentation processes as logarithmically correlated random fields.

     

The Scaling Limit of the KPZ Equation in Space Dimension 3 and Higher

We study in the present article the Kardar–Parisi–Zhang (KPZ) equation

$$\begin{aligned} \partial _t h(t,x)=\nu \Delta h(t,x)+\lambda |\nabla h(t,x)|^2 +\sqrt{D}\, \eta (t,x), \qquad (t,x)\in \mathbb {R}_+\times \mathbb {R}^d \end{aligned}$$

in \(d\ge 3\) dimensions in the perturbative regime, i.e. for \(\lambda >0\) small enough and a smooth, bounded, integrable initial condition \(h_0=h(t=0,\cdot )\). The forcing term \(\eta \) in the right-hand side is a regularized space-time white noise. The exponential of h—its so-called Cole-Hopf transform—is known to satisfy a linear PDE with multiplicative noise. We prove a large-scale diffusive limit for the solution, in particular a time-integrated heat-kernel behavior for the covariance in a parabolic scaling. The proof is based on a rigorous implementation of K. Wilson’s renormalization group scheme. A double cluster/momentum-decoupling expansion allows for perturbative estimates of the bare resolvent of the Cole-Hopf linear PDE in the small-field region where the noise is not too large, following the broad lines of Iagolnitzer and Magnen (Commun Math Phys 162(1):85–121, 1994). Standard large deviation estimates for \(\eta \) make it possible to extend the above estimates to the large-field region. Finally, we show, by resumming all the by-products of the expansion, that the solution h may be written in the large-scale limit (after a suitable Galilei transformation) as a small perturbation of the solution of the underlying linear Edwards–Wilkinson model (\(\lambda =0\)) with renormalized coefficients \(\nu _{eff}=\nu +O(\lambda ^2),D_{eff}=D+O(\lambda ^2)\).

     

Kinetically modified non-minimal inflation with exponential frame function

We consider supersymmetric (SUSY) and non-SUSY models of chaotic inflation based on the \(\phi ^n\) potential with \(n=2\) or 4. We show that the coexistence of an exponential non-minimal coupling to gravity \(f_\mathcal{R}=\mathrm{e}^{c_\mathcal{R}\phi ^{p}}\) with a kinetic mixing of the form \(f_{\mathrm{K}}=c_{\mathrm{K}}f_\mathcal{R}^m\) can accommodate inflationary observables favored by the Planck and Bicep2/Keck Array results for \(p=1\) and 2, \(1\le m\le 15\) and \(2.6\times 10^{-3}\le r_{\mathcal {R}\mathrm{K}}=c_\mathcal{R}/c_{\mathrm{K}}^{p/2}\le 1,\) where the upper limit is not imposed for \(p=1\). Inflation is of hilltop type and it can be attained for subplanckian inflaton values with the corresponding effective theories retaining the perturbative unitarity up to the Planck scale. The supergravity embedding of these models is achieved employing two chiral gauge singlet supefields, a monomial superpotential and several (semi)logarithmic or semi-polynomial Kähler potentials.     

Unbounded Trace Orbits of Thue–Morse Hamiltonian

It is well known that, an energy is in the spectrum of Fibonacci Hamiltonian if and only if the corresponding trace orbit is bounded. However, it is not known whether the same result holds for the Thue–Morse Hamiltonian. In this paper, we give a negative answer to this question. More precisely, we construct two subsets \(\Sigma _{II}\) and \(\Sigma _{III}\) of the spectrum of the Thue–Morse Hamiltonian, both of which are dense and uncountable, such that each energy in \(\Sigma _{II}\cup \Sigma _{III}\) corresponds to an unbounded trace orbit. Exact estimates on the norm of the transfer matrices are also obtained for these energies: for \(E\in \Sigma _{II}\cup \Sigma _{III}, \) the norms of the transfer matrices behave like

$$\begin{aligned} e^{c_1\gamma \sqrt{n}}\le \Vert T_{ n}(E)\Vert \le e^{c_2\gamma \sqrt{n}}. \end{aligned}$$

However, two types of energies are quite different in the sense that each energy in \(\Sigma _{II}\) is associated with a two-sided pseudo-localized state, while each energy in \(\Sigma _{III}\) is associated with a one-sided pseudo-localized state. The difference is also reflected by the local dimensions of the spectral measure: the local dimension is 0 for energies in \(\Sigma _{II}\) and is larger than 1 for energies in \(\Sigma _{III}.\) As a comparison, we mention another known countable dense subset \(\Sigma _I\). Each energy in \(\Sigma _I\) corresponds to an eventually constant trace map and the associated eigenvector is an extended state. In summary, the Thue–Morse Hamiltonian exhibits “mixed spectral nature”.

     

Topological Expansion in the Complex Cubic Log–Gas Model: One-Cut Case

We prove the topological expansion for the cubic log–gas partition function

$$\begin{aligned} Z_N(t)= \int _\Gamma \cdots \int _\Gamma \prod _{1\le j<k\le N}(z_j-z_k)^2 \prod _{k=1}^Ne^{-N\left( -\frac{z^3}{3}+tz\right) }\mathrm{dz}_1\cdots \mathrm{dz}_N, \end{aligned}$$

where t is a complex parameter and \(\Gamma \) is an unbounded contour on the complex plane extending from \(e^{\pi \mathrm{i}}\infty \) to \(e^{\pi \mathrm{i}/3}\infty \). The complex cubic log–gas model exhibits two phase regions on the complex t-plane, with one cut and two cuts, separated by analytic critical arcs of the two types of phase transition: split of a cut and birth of a cut. The common point of the critical arcs is a tricritical point of the Painlevé I type. In the present paper we prove the topological expansion for \(\log Z_N(t)\) in the one-cut phase region. The proof is based on the Riemann–Hilbert approach to semiclassical asymptotic expansions for the associated orthogonal polynomials and the theory of S-curves and quadratic differentials.

     

Domain and Range of the Modified Wave Operator for Schrödinger Equations with a Critical Nonlinearity

We study the final problem for the nonlinear Schrödinger equation

$i{\partial }_{t}u+\frac{1}{2}\Delta u=\lambda|u|^{\frac{2}{n}}u,\quad (t,x)\in {\mathbf{R}}\times \mathbf{R}^{n},$

where\(\lambda \in{\bf R},n=1,2,3\). If the final data\(u_{+}\in {\bf H}^{0,\alpha }=\left\{ \phi \in {\bf L}^{2}:\left( 1+\left\vert x\right\vert \right) ^{\alpha }\phi \in {\bf L}^{2}\right\} \) with\(\frac{ n}{2} < \alpha < \min \left( n,2,1+\frac{2}{n}\right) \) and the norm\(\Vert \widehat{u_{+}}\Vert _{{\bf L}^{\infty }}\) is sufficiently small, then we prove the existence of the wave operator in L ². We also construct the modified scattering operator from H ^0,α to H ^0,δ with\(\frac{n}{2} < \delta < \alpha\).

     

Two Higgs doublet models and $$b\rightarrow s$$ exclusive decays

We computed the leading order Wilson coefficients relevant to all the exclusive \(b\rightarrow s\ell ^+\ell ^-\) decays in the framework of the two Higgs doublet model (2HDM) with a softly broken \(\mathbb {Z}_2\) symmetry by including the \(\mathcal {O}(m_b)\) corrections. We elucidate the issue of appropriate matching between the full and the effective theory when dealing with the (pseudo-)scalar operators for which keeping the external momenta different from zero is necessary. We then make a phenomenological analysis by using the measured \({\mathcal {B}}(B_s\rightarrow \mu ^+\mu ^-)\) and \({\mathcal {B}}(B\rightarrow K \mu ^+\mu ^-)_{\mathrm {high}-q^2}\), for which the hadronic uncertainties are well controlled, and we discuss their impact on various types of 2HDM. A brief discussion of the decays with \(\tau \)-leptons in the final state is provided too.     

The effects of magnetic field and hydrogen-like impurity on RbCl quantum pseudodot qubit

By using a variational method of Pekar type, we investigate the effects of the hydrogen-like impurity and magnetic field on the electron’s probability density (PD) and oscillating frequency (OF) of a RbCl quantum pseudodot qubit. Numerical results indicate that (1) the PD oscillates periodically; (2) the crest of the PD will decrease with increasing the cyclotron frequencies and the Coulombic impurity potential strength; (3) as the cyclotron frequency of the magnetic field and the strength of the Coulombic impurity potential increases, PD’s peaks will occur more frequently; (4) besides, Figs. 1b and 2b clearly show that in a single period the PD will decrease with increasing the cyclotron frequency and the Coulombic impurity potential strength when \( t > 1.8\;\text{fs} \); whereas the changing law is just the opposite when \( t < 1.8\;\text{fs} \); (5) the OF is an aggrandizing function of the strength of the Coulombic impurity potential, whereas it is a decaying one of the cyclotron frequencies of the magnetic field. The coherence of qubit is crucial to the investigations of quantum information and quantum computation, where the electron’s PD, the OF and the coherence time are the physical quantities representing the properties of coherence. Our research results fine that by changing the cyclotron frequency of the magnetic field and the strength of the Coulombic impurity potential one can adjust the electron’s PD and the OF.

Open image in new window

Fig. 1 The PD \( \text{Q}\left( {r,t} \right) \) versus the time \( t \) and the cyclotron frequency of the magnetic field \( \omega_{c} \) with \( \text{V}_{0} = 10.0\,\text{meV, r}_{0} = 1.0\,\text{nm, }\beta \text{ = 1.0}\,\text{meV} \cdot \text{nm} \) and \( x = y = z = 1.0\,\text{nm} \)

Open image in new window

Fig. 2 The PD \( \text{Q}\left( {r,t} \right) \) versus the time \( t \) and strength of the Coulombic impurity potential \( \beta \) with \( \text{V}_{0} = 10.0\,\text{meV, r}_{0} = 1.0\,\text{nm,} \, \omega_{c} \text{ = 2.0}\, \times \text{10}^{13}\,\text{Hz} \) and \( x = y = z = 1.0\,\text{nm} \)

     

$$\varvec{B^0\rightarrow K^{*0}\mu ^+\mu ^-}$$ decay in the aligned two-Higgs-doublet model

In the aligned two-Higgs-doublet model, we perform a complete one-loop computation of the short-distance Wilson coefficients \(C_{7,9,10}^{(\prime )}\), which are the most relevant ones for \(b\rightarrow s\ell ^+\ell ^-\) transitions. It is found that, when the model parameter \(\left| \varsigma _{u}\right| \) is much smaller than \(\left| \varsigma _{d}\right| \), the charged scalar contributes mainly to chirality-flipped \(C_{9,10}^\prime \), with the corresponding effects being proportional to \(\left| \varsigma _{d}\right| ^2\). Numerically, the charged-scalar effects fit into two categories: (A) \(C_{7,9,10}^\mathrm {H^\pm }\) are sizable, but \(C_{9,10}^{\prime \mathrm {H^\pm }}\simeq 0\), corresponding to the (large \(\left| \varsigma _{u}\right| \), small \(\left| \varsigma _{d}\right| \)) region; (B) \(C_7^\mathrm {H^\pm }\) and \(C_{9,10}^{\prime \mathrm {H^\pm }}\) are sizable, but \(C_{9,10}^\mathrm {H^\pm }\simeq 0\), corresponding to the (small \(\left| \varsigma _{u}\right| \), large \(\left| \varsigma _{d}\right| \)) region. Taking into account phenomenological constraints from the inclusive radiative decay \(B\rightarrow X_{s}{\gamma }\), as well as the latest model-independent global analysis of \(b\rightarrow s\ell ^+\ell ^-\) data, we obtain the much restricted parameter space of the model. We then study the impact of the allowed model parameters on the angular observables \(P_2\) and \(P_5'\) of \(B^0\rightarrow K^{*0}\mu ^+\mu ^-\) decay, and we find that \(P_5'\) could be increased significantly to be consistent with the experimental data in case B.     

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.