Courant-sharp eigenvalues of Neumann 2-rep-tiles

We find the Courant-sharp Neumann eigenvalues of the Laplacian on some 2-rep-tile domains. In \(\mathbb {R}^{2}\), the domains we consider are the isosceles right triangle and the rectangle with edge ratio \(\sqrt{2}\) (also known as the A4 paper). In \(\mathbb {R}^{n}\), the domains are boxes which generalize the mentioned planar rectangle. The symmetries of those domains reveal a special structure of their eigenfunctions, which we call folding\unfolding. This structure affects the nodal set of the eigenfunctions, which, in turn, allows to derive necessary conditions for Courant-sharpness. In addition, the eigenvalues of these domains are arranged as a lattice which allows for a comparison between the nodal count and the spectral position. The Courant-sharpness of most eigenvalues is ruled out using those methods. In addition, this analysis allows to estimate the nodal deficiency—the difference between the spectral position and the nodal count.     

Equilibration in the Kac Model Using the GTW Metric $$d_2$$

We use the Fourier based Gabetta–Toscani–Wennberg metric \(d_2\) to study the rate of convergence to equilibrium for the Kac model in 1 dimension. We take the initial velocity distribution of the particles to be a Borel probability measure \(\mu \) on \(\mathbb {R}^n\) that is symmetric in all its variables, has mean \(\vec {0}\) and finite second moment. Let \(\mu _t(dv)\) denote the Kac-evolved distribution at time t, and let \(R_\mu \) be the angular average of \(\mu \). We give an upper bound to \(d_2(\mu _t, R_\mu )\) of the form \(\min \left\{ B e^{-\frac{4 \lambda _1}{n+3}t}, d_2(\mu ,R_\mu )\right\} ,\) where \(\lambda _1 = \frac{n+2}{2(n-1)}\) is the gap of the Kac model in \(L^2\) and B depends only on the second moment of \(\mu \). We also construct a family of Schwartz probability densities \(\{f_0^{(n)}: \mathbb {R}^n\rightarrow \mathbb {R}\}\) with finite second moments that shows practically no decrease in \(d_2(f_0(t), R_{f_0})\) for time at least \(\frac{1}{2\lambda }\) with \(\lambda \) the rate of the Kac operator. We also present a propagation of chaos result for the partially thermostated Kac model in Tossounian and Vaidyanathan (J Math Phys 56(8):083301, 2015).     

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.     

Eigenvalue Order Statistics for Random Schrödinger Operators with Doubly-Exponential Tails

We consider random Schrödinger operators of the form \({\Delta+\xi}\), where \({\Delta}\) is the lattice Laplacian on \({\mathbb{Z}^{d}}\) and \({\xi}\) is an i.i.d. random field, and study the extreme order statistics of the Dirichlet eigenvalues for this operator restricted to large but finite subsets of \({\mathbb{Z}^{d}}\). We show that, for \({\xi}\) with a doubly-exponential type of upper tail, the upper extreme order statistics of the eigenvalues falls into the Gumbel max-order class, and the corresponding eigenfunctions are exponentially localized in regions where \({\xi}\) takes large, and properly arranged, values. The picture we prove is thus closely connected with the phenomenon of Anderson localization at the spectral edge. Notwithstanding, our approach is largely independent of existing methods for proofs of Anderson localization and it is based on studying individual eigenvalue/eigenfunction pairs and characterizing the regions where the leading eigenfunctions put most of their mass.     

Critical Two-Point Function for Long-Range <Emphasis Type="Italic">O</Emphasis>(<Emphasis Type="Italic">n</Emphasis>) Models Below the Upper Critical Dimension

We consider the n-component \(|\varphi |^4\) lattice spin model (\(n \ge 1\)) and the weakly self-avoiding walk (\(n=0\)) on \(\mathbb Z^d\), in dimensions \(d=1,2,3\). We study long-range models based on the fractional Laplacian, with spin-spin interactions or walk step probabilities decaying with distance r as \(r^{-(d+\alpha )}\) with \(\alpha \in (0,2)\). The upper critical dimension is \(d_c=2\alpha \). For \(\varepsilon >0\), and \(\alpha = \frac{1}{2} (d+\varepsilon )\), the dimension \(d=d_c-\varepsilon \) is below the upper critical dimension. For small \(\varepsilon \), weak coupling, and all integers \(n \ge 0\), we prove that the two-point function at the critical point decays with distance as \(r^{-(d-\alpha )}\). This “sticking” of the critical exponent at its mean-field value was first predicted in the physics literature in 1972. Our proof is based on a rigorous renormalisation group method. The treatment of observables differs from that used in recent work on the nearest-neighbour 4-dimensional case, via our use of a cluster expansion.     

Truncated Linear Statistics Associated with the Eigenvalues of Random Matrices II. Partial Sums over Proper Time Delays for Chaotic Quantum Dots

Invariant ensembles of random matrices are characterized by the distribution of their eigenvalues \(\{\lambda _1,\ldots ,\lambda _N\}\). We study the distribution of truncated linear statistics of the form \(\tilde{L}=\sum _{i=1}^p f(\lambda _i)\) with \(p<N\). This problem has been considered by us in a previous paper when the p eigenvalues are further constrained to be the largest ones (or the smallest). In this second paper we consider the same problem without this restriction which leads to a rather different analysis. We introduce a new ensemble which is related, but not equivalent, to the “thinned ensembles” introduced by Bohigas and Pato. This question is motivated by the study of partial sums of proper time delays in chaotic quantum dots, which are characteristic times of the scattering process. Using the Coulomb gas technique, we derive the large deviation function for \(\tilde{L}\). Large deviations of linear statistics \(L=\sum _{i=1}^N f(\lambda _i)\) are usually dominated by the energy of the Coulomb gas, which scales as \(\sim N^2\), implying that the relative fluctuations are of order 1 / N. For the truncated linear statistics considered here, there is a whole region (including the typical fluctuations region), where the energy of the Coulomb gas is frozen and the large deviation function is purely controlled by an entropic effect. Because the entropy scales as \(\sim N\), the relative fluctuations are of order \(1/\sqrt{N}\). Our analysis relies on the mapping on a problem of p fictitious non-interacting fermions in N energy levels, which can exhibit both positive and negative effective (absolute) temperatures. We determine the large deviation function characterizing the distribution of the truncated linear statistics, and show that, for the case considered here (\(f(\lambda )=1/\lambda \)), the corresponding phase diagram is separated in three different phases.     

Universality Conjecture and Results for a Model of Several Coupled Positive-Definite Matrices

The paper contains two main parts: in the first part, we analyze the general case of \({p \geq 2}\) matrices coupled in a chain subject to Cauchy interaction. Similarly to the Itzykson-Zuber interaction model, the eigenvalues of the Cauchy chain form a multi level determinantal point process. We first compute all correlations functions in terms of Cauchy biorthogonal polynomials and locate them as specific entries of a \({(p+1) \times (p+1)}\) matrix valued solution of a Riemann–Hilbert problem. In the second part, we fix the external potentials as classical Laguerre weights. We then derive strong asymptotics for the Cauchy biorthogonal polynomials when the support of the equilibrium measures contains the origin. As a result, we obtain a new family of universality classes for multi-level random determinantal point fields, which include the Bessel\({_{\nu}}\) universality for 1-level and the Meijer-G universality for 2-level. Our analysis uses the Deift-Zhou nonlinear steepest descent method and the explicit construction of a \({(p+1) \times (p+1)}\) origin parametrix in terms of Meijer G-functions. The solution of the full Riemann–Hilbert problem is derived rigorously only for p = 3 but the general framework of the proof can be extended to the Cauchy chain of arbitrary length p.     

Universality for 1d Random Band Matrices: Sigma-Model Approximation

The paper continues the development of the rigorous supersymmetric transfer matrix approach to the random band matrices started in (J Stat Phys 164:1233–1260, 2016; Commun Math Phys 351:1009–1044, 2017). We consider random Hermitian block band matrices consisting of \(W\times W\) random Gaussian blocks (parametrized by \(j,k \in \Lambda =[1,n]^d\cap \mathbb {Z}^d\)) with a fixed entry’s variance \(J_{jk}=\delta _{j,k}W^{-1}+\beta \Delta _{j,k}W^{-2}\), \(\beta >0\) in each block. Taking the limit \(W\rightarrow \infty \) with fixed n and \(\beta \), we derive the sigma-model approximation of the second correlation function similar to Efetov’s one. Then, considering the limit \(\beta , n\rightarrow \infty \), we prove that in the dimension \(d=1\) the behaviour of the sigma-model approximation in the bulk of the spectrum, as \(\beta \gg n\), is determined by the classical Wigner–Dyson statistics.     

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}}\).     

Universal Components of Random Nodal Sets

We give, as L grows to infinity, an explicit lower bound of order \({L^{\frac{n}{m}}}\) for the expected Betti numbers of the vanishing locus of a random linear combination of eigenvectors of P with eigenvalues below L. Here, P denotes an elliptic self-adjoint pseudo-differential operator of order \({m > 0}\), bounded from below and acting on the sections of a Riemannian line bundle over a smooth closed n-dimensional manifold M equipped with some Lebesgue measure. In fact, for every closed hypersurface \({\Sigma}\) of \({\mathbb{R}^n}\), we prove that there exists a positive constant \({p_\Sigma}\) depending only on \({\Sigma}\), such that for every large enough L and every \({x \in M}\), a component diffeomorphic to \({\Sigma}\) appears with probability at least \({p_\Sigma}\) in the vanishing locus of a random section and in the ball of radius \({L^{-\frac{1}{m}}}\) centered at x. These results apply in particular to Laplace–Beltrami and Dirichlet-to-Neumann operators.     

Non-Real Eigenvalues for {{\mathcal{PT}}}-Symmetric Double Wells

We study small, \({{\mathcal{PT}}}\)-symmetric perturbations of self-adjoint double-well Schrödinger operators in dimension \({n\ge 1}\). We prove that the eigenvalues stay real for a very small perturbation, then bifurcate to the complex plane as the perturbation gets stronger.     

Compound Poisson Law for Hitting Times to Periodic Orbits in Two-Dimensional Hyperbolic Systems

We show that a compound Poisson distribution holds for scaled exceedances of observables \(\phi \) uniquely maximized at a periodic point \(\zeta \) in a variety of two-dimensional hyperbolic dynamical systems with singularities \((M,T,\mu )\), including the billiard maps of Sinai dispersing billiards in both the finite and infinite horizon case. The observable we consider is of form \(\phi (z)=-\ln d(z,\zeta )\) where d is a metric defined in terms of the stable and unstable foliation. The compound Poisson process we obtain is a Pólya-Aeppli distibution of index \(\theta \). We calculate \(\theta \) in terms of the derivative of the map T. Furthermore if we define \(M_n=\max \{\phi ,\ldots ,\phi \circ T^n\}\) and \(u_n (\tau )\) by \(\lim _{n\rightarrow \infty } n\mu (\phi >u_n (\tau ) )=\tau \) the maximal process satisfies an extreme value law of form \(\mu (M_n \le u_n)=e^{-\theta \tau }\). These results generalize to a broader class of functions maximized at \(\zeta \), though the formulas regarding the parameters in the distribution need to be modified.     

Hölder Regularity of the Solutions of the Cohomological Equation for Roth Type Interval Exchange Maps

There is a general method for constructing a soliton hierarchy from a splitting \({{L_{\pm}}}\) of a loop group as positive and negative sub-groups together with a commuting linearly independent sequence in the positive Lie algebra \({\mathcal{L}_{+}}\). Many known soliton hierarchies can be constructed this way. The formal inverse scattering associates to each f in the negative subgroup \({L_-}\) a solution \({u_{f}}\) of the hierarchy. When there is a 2 co-cycle of the Lie algebra that vanishes on both sub-algebras, Wilson constructed a tau function \({\tau_{f}}\) for each element \({f \in L_-}\). In this paper, we give integral formulas for variations of \({\ln\tau_{f}}\) and second partials of \({\ln\tau_{f}}\), discuss whether we can recover solutions \({u_{f}}\) from \({\tau_{f}}\), and give a general construction of actions of the positive half of the Virasoro algebra on tau functions. We write down formulas relating tau functions and formal inverse scattering solutions and the Virasoro vector fields for the \({GL(n,\mathbb{C})}\)-hierarchy.     

Neutron Pairing Correlations in an {\alpha}–{n}–{n} Three-Cluster Model of the {^{6}}He Nucleus

An \({\alpha}\)–n–n three-cluster model of the \({^6}\)He nucleus is studied by solving the Faddeev equations, where the cluster potential between \({\alpha}\) and n takes into account the Pauli exclusion correction, using the Fish-Bone Optical Model (Schmid in Z Phys A 297:105, 1980). The resulting binding energy of the ground state (\({0^+}\)) is 0.831 MeV and the resonance energy of the first excited state (\({2^+}\)), 0.60–i0.012 MeV, is extracted from the three-cluster break-up threshold. These theoretical values are in reasonable agreement with the experimental data: 0.973 MeV and 0.824–i0.056 MeV, respectively. In order to investigate the structure of these states, we calculate the angle density matrix for the \({\angle n_1 \alpha n_2}\) angle in the triangle formed by the three clusters. The angle density matrix of the ground state has two peaks and the configuration of \({0^+}\) wave function corresponding to the peaks constitutes a mixture of an acute-angled triangle structure and an obtuse-angled one. This finding is consistent with the former result from a variational approach (Hagino and Sagawa in Phys Rev C 72:044321, 2005). On the other hand, in the case of \({2^+}\) state only a single peak is obtained.     

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.     

Dynamical system approach to scalar–vector–tensor cosmology

Using scalar–vector–tensor Brans Dicke (VBD) gravity (Ghaffarnejad in Gen Relativ Gravit 40:2229, 2008; Gen Relativ Gravit 41:2941, 2009) in presence of self interaction BD potential \(V(\phi )\) and perfect fluid matter field action we solve corresponding field equations via dynamical system approach for flat Friedmann Robertson Walker metric (FRW). We obtained three type critical points for \(\Lambda CDM\) vacuum de Sitter era where stability of our solutions are depended to choose particular values of BD parameter \(\omega \). One of these fixed points is supported by a constant potential which is stable for \(\omega <0\) and behaves as saddle (quasi stable) for \(\omega \ge 0\). Two other ones are supported by a linear potential \(V(\phi )\sim \phi \) which one of them is stable for \(\omega =0.27647\). For a fixed value of \(\omega \) there is at least 2 out of 3 critical points reaching to a unique critical point. Namely for \(\omega =-0.16856(-0.56038)\) the second (third) critical point become unique with the first critical point. In dust and radiation eras we obtained one critical point which never become unique fixed point. In the latter case coordinates of fixed points are also depended to \(\omega \). To determine stability of our solutions we calculate eigenvalues of Jacobi matrix of 4D phase space dynamical field equations for de Sitter, dust and radiation eras. We should point also potentials which support dust and radiation eras must be similar to \(V(\phi )\sim \phi ^{-\frac{1}{2}}\) and \(V(\phi )\sim \phi ^{-1}\) respectively. In short our study predicts that radiation and dust eras of our VBD–FRW cosmology transmit to stable de Sitter state via non-constant potential (effective variable cosmological parameter) by choosing \(\omega =0.27647\).     

Aharonov and Bohm versus Welsh eigenvalues

We consider a class of two-dimensional Schrödinger operator with a singular interaction of the \(\delta \) type and a fixed strength \(\beta \) supported by an infinite family of concentric, equidistantly spaced circles, and discuss what happens below the essential spectrum when the system is amended by an Aharonov–Bohm flux \(\alpha \in [0,\frac{1}{2}]\) in the center. It is shown that if \(\beta \ne 0\), there is a critical value \(\alpha _{\mathrm {crit}}\in (0,\frac{1}{2})\) such that the discrete spectrum has an accumulation point when \(\alpha <\alpha _{\mathrm {crit}}\), while for \(\alpha \ge \alpha _{\mathrm {crit}}\) the number of eigenvalues is at most finite, in particular, the discrete spectrum is empty for any fixed \(\alpha \in (0,\frac{1}{2})\) and \(|\beta |\) small enough.     

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))\).

     

Exotic tetraquark states with the $$qq\bar{Q}\bar{Q}$$ configuration

In this work, we study systematically the mass splittings of the \(qq\bar{Q}\bar{Q}\) (\(q=u\), d, s and \(Q=c\), b) tetraquark states with the color-magnetic interaction by considering color mixing effects and estimate roughly their masses. We find that the color mixing effect is relatively important for the \(J^P=0^+\) states and possible stable tetraquarks exist in the \(nn\bar{Q}\bar{Q}\) (\(n=u\), d) and \(ns\bar{Q}\bar{Q}\) systems either with \(J=0\) or with \(J=1\). Possible decay patterns of the tetraquarks are briefly discussed.     

CGC/saturation approach: re-visiting the problem of odd harmonics in angular correlations

In this paper we demonstrate that the selection of events with different multiplicities of produced particles, leads to the violation of the azimuthal angular symmetry, \(\phi \rightarrow \pi - \phi \). We find for LHC and lower energies, that this violation can be so large for the events with multiplicities \(n \ge 2 \bar{n}\), where \(\bar{n}\) is the mean multiplicity, that it leads to almost no suppression of \(v_n\), with odd n. However, this can only occur if the typical size of the dipole in DIS with a nuclear target is small, or \(Q^2 \,>\,Q^2_s\left( A; Y_{\mathrm{min}},b\right) \), where \(Q_s\) is the saturation momentum of the nucleus at \(Y = Y_{\mathrm{min}}\). In the case of large sizes of dipoles, when \(Q^2 \,<\,Q^2_s\left( A; Y_{\mathrm{min}},b\right) \), we show that \(v_n =0\) for odd n. Hadron-nucleus scattering is discussed.