Papers

Among the various "no-go theorems" in the foundations of quantum mechanics, Bell inequalities are among the few, which do not implicitly rely on the mathematical structure of quantum mechanics. Instead, Bell inequalities are commonly understood to rest on Kolmogorov's theory of probability and a "locality" assumption. Physical evidence supporting their violation is hence taken as evidence against at least one of these two assumptions. This work provides a comprehensive analysis of the physical debate on Bell inequalities from the point of view of Kolmogorov's theory of probability, focusing on the Bell-Clauser-Horn-Shimony-Holt inequality (BCHSH inequality). The central contribution is the motivation and mathematical study of the notion of contextuality in Bell experiments: in so-called contextual Bell models the detectors themselves behave randomly and the distribution of the random detector parameters depends on the chosen detector settings. It is shown how even "classical" Bell experiments motivate such a description. The need to account for contextuality was pointed out before by de la Peña, Cetto, Brody, Lochak, Bohm, Hiley and others. A broad class of contextual models is proven to satisfy the BCHSH inequality, thereby generalizing a previous result by Gill and Lambare. Nonetheless, even under strong locality requirements, contextual Bell models may violate the BCHSH inequality: based on work by Popescu and Rohrlich, we consider a "proof of concept" example exhibiting the maximum possible violation. Contrary to the much debated "loopholes" in empirical tests of Bell inequalities, contextuality thus reveals a conceptual misunderstanding in the physical conclusions commonly drawn from Bell inequality violations.

We present a direct simulation Monte Carlo (DSMC) method for the kinetic theory of ordered fluids recently introduced by the authors. The microscopic constituents of these fluids carry an order parameter belonging to a manifold. The method applies Strang splitting to couple a Monte Carlo treatment of the binary collisions with a symplectic integration of the mean-field Vlasov transport on the order-parameter manifold. We use the Nanbu-Babovsky approach for Maxwellian collision kernels, and Bird's acceptance-rejection method for general ones. The mean-field Vlasov force is reconstructed from the particle ensemble at each time step. We also include a Bhatnagar-Gross-Krook (BGK) / Andersen thermostat to control the system temperature. We formulate the method for a generic order-parameter manifold, and specialize it to two-dimensional calamitic (rod-like) molecules. Numerical tests validate the method. We confirm the relaxation to the Maxwellian distribution predicted by the H-theorem, and simulate the mean-field alignment driven by quadratic, Kuramoto, and Onsager interaction potentials. With the Onsager potential we successfully capture the isotropic-nematic phase transition at the correct critical temperature.

The Teukolsky Master Equation (TME) is a powerful technique used to study perturbations of a spacetime with gravitational, electromagnetic, and neutrino components. The well-known Kerr metric is a suitable candidate for this formalism because it is a Petrov type D metric. However, the Kerr metric represents the exterior spacetime of a perfect sphere and only black holes are described using this metric. For another astrophysical object, such as fast-spinning neutron stars, the Kerr metric is not valid due to the rotation induced deformations to the gravitational source. The main goal of the present work is to show a TME for a Kerr-like metric (KL) with a mass quadrupole moment parameter $q$. This new parameter will represent deformations from a perfect sphere; therefore, KL is a better candidate to study perturbations surrounding spinning stars.

We consider an interpolation with the parameter $p\in [0,1]$ between the discrete-time quantum ($p=0$) and its intrinsic random ($p=1$) walks on a finite connected graph. Its time evolution is defined by the convex combination of the Kraus (CPTP) maps of the quantum and random walks. The Kraus map of the random walk is represented by taking a kind of projection of that of the quantum walk. The time evolution of the interpolation is interpreted as the combination of the following two dynamics of $2$ walkers on the same graph correlated with each other: $2$ walkers move independently of each other (no-correlation) with probability $1-p$, while $2$ walkers always move into the same position (the maximal correlation) with probability $p$ at each time step. In this paper, we generalize the random walk to a new walk, namely, the correlated walk, which preserves the intrinsic property and relaxes the maximal correlation of the random walk. This correlation is determined by a partition of the arc sets in the underlying graph. We show that the eigenvalues of the interpolating walk between the correlated and quantum walks live in $\{ z\in \mathbb{C}\;|\;1-p\leq |z|\leq 1 \}$ and that the absorption state coincides with that of the correlated walk, which is characterized by the underlying graph structures.

We study the Cauchy problem for the Tzitzéica equation, which is an important integrable model arising in affine differential geometry and characterizing proper affine spheres. Recently, Huang, Wang and Zhu (Math. Ann. 395, 18 (2026)) reported the long-time asymptotic result for the Tzitzéica equation for the case of the purely continuous spectrum. Inspired by their work, we conduct in-depth research on the soliton resolution conjecture and asymptotic analysis of the Tzitzéica equation. Compared with the previous work, by developing the $\bar{\partial}$-nonlinear steepest descent method, we reveal that the solution of the Tzitzéica equation exhibits three different asymptotic regions depending on $\xi:=x/t$. For the region $\xi\in(-\infty,-1]\cup[1,\infty)$, where there are no stationary phase points, we rigorously prove that the solution of the Tzitzéica equation decays algebraically to zero with error $O(t^{-1})$. In the region $\xi\in(-1,1)$, the phase function $\theta(z)$ has two stationary phase points. The corresponding asymptotic approximations can be characterized with an $N$-soliton solution as well as an interaction term between soliton solutions and the dispersion term with a residual error of order $O(t^{-3/4})$. In the transition regions $\{(x,t):1-\varepsilon<|\xi|<1\}$, the corresponding asymptotic approximation can be characterized by the $N$-soliton solution together with the interaction terms between the soliton solution and the dispersion term, with a residual error of order $O(t^{-3/4})$. Our results provide a rigorous analysis of the long-time asymptotic behavior of solutions of the Tzitzéica equation in different regions and confirm the soliton resolution conjecture for the Tzitzéica equation.

We associate to finite block size data \[ \bm n=(n_d)_{d\ge1},\qquad n_1>0, \] a positive definite cyclotomic form and a Householder reflection representation. Writing \[ N=\operatorname{lcm}\{d:n_d\ne0\},\qquad \wt=\sum_d d^2n_d, \] the form on the cyclic group $C_N$ is \[ B_{\bm n}(g,h)=\frac1{\wt}\sum_{(g^{-1}h)^d=1}d^2n_d. \] Its labelled Fourier spectrum reconstructs the complete multiplicity sequence by Möbius inversion. If the $n_d$ are the Wedderburn multiplicities of a split semisimple algebra $S=\prod_d M_d(k)^{n_d}$, the divisibility widths are dimensions of characteristic sums of Wedderburn blocks, and the form together with $\dim_kS$ recovers the isomorphism type of $S$. The main theorem is an arithmetic rigidity result. For $N\ge2$, the reflection group is finite if and only if it preserves a full rank lattice, equivalently \[ \supp(\bm n)=\{1,N\},\qquad n_1=N^2n_N. \] In this case $B_{\bm n}=\frac12(I+J)$ and the group is $W(A_N)\cong S_{N+1}$. In semisimple algebra terms, finiteness is equivalent to \[ S\cong k^{N^2m}\times M_N(k)^m \] for some $m\ge1$. For $N=2$, every nonfinite case is a dense copy of $D_\infty$ in the corresponding orthogonal group. For contraction algebras of threefold flops, Toda's width formulas identify the total scale and primitive Fourier mode with $\dim A_{\con}$ and $\dim A_{\con}^{\ab}$. This yields infinitude for all irreducible flops of length at least three and, in length two, \[ \Gamma_f\text{ finite}\quad\Longleftrightarrow\quad \dim A_{\con}=2\dim A_{\con}^{\ab}. \] The final part gives a complementary BPS/Riemann–Hilbert realization of the same reflection matrices.

The para-Krawtchouk polynomials arose in the search for spin chains with perfect state transfer; the para-Racah, $q$-para-Racah and para-Bannai–Ito polynomials followed, obtained through singular truncations of families of the Askey scheme and used in turn to design spin chains. All are orthogonal on bi-lattices; the prefix goes back to the para-Krawtchouk case, whose spectrum is that of the parabose oscillator. The representations that these polynomials provide of the corresponding algebras — Hahn, Racah, $q$-Racah, Bannai–Ito and complementary Bannai–Ito — are shown to be reducible. Each decomposes into a direct sum of two irreducible modules of the same algebra, one attached to each of the two interlaced sub-lattices, and the parameters of the submodules are given. The para polynomials are thereby expressed in terms of the classical polynomials that the submodules carry, and the pairing $\lambda_n=\lambda_{N-n}$ of their eigenvalues is accounted for; the three features of the bi-lattice setting on which the reduction rests are identified. The properties of the para families are then presented in a unified way in the light of this decomposition, together with the constructions through which these families arise and the settings in which they have been put to use.

Stochastic quantization formulates $d$-dimensional Euclidean QFT as the equilibrium limit of a $(d+1)$-dimensional Langevin system. Using the supersymmetric formulation of stochastic systems, we apply techniques from Morse theory to study Langevin dynamics. We interpret stochastic quantization as the absolutization of a relative QFT by a bulk cohomological topological Equilibrium TFT ("EqmTFT"), and relate its enriched Neumann boundary data to the $d$-dimensional BV/BRST complex. Analytic continuation reveals connections to the A-model and non-relativistic and Carrollian limits.

The para-Bannai–Ito polynomials are orthogonal on a Bannai–Ito bi-lattice that is two-sided — it extends on both sides of a central value $c_0$ — and depend, besides the two grid parameters, on an isospectral deformation parameter. Their Jacobi matrix $J$ is invariant under the reversal of the sites at one value of that parameter only, the persymmetric point, and we determine its continuum limit there. Because the grid is two-sided, the natural object is the shifted operator $Q=J-c_0$, and the limit is taken about the centre of the band: $Q$ contracts to a first–order Dunkl operator, a supercharge, and $H=Q^2$ to the generalized Pöschl–Teller (Scarf I) Hamiltonian, with little $-1$ Jacobi eigenfunctions — supersymmetric quantum mechanics with reflections, the reflection entering the square root of the Hamiltonian. The second grid parameter adds a reflection term to the Hamiltonian itself, and the bi-lattice is, up to a shift, the spectrum of the limiting supercharge. The convergence of the eigenvectors is proved directly on the polynomials. On this grid the persymmetric model transfers states perfectly, under an arithmetic condition on the parameters, but never revives fractionally, revival being the work of the deformation; both properties pass to the limit. These results are for an even number of sites; for an odd number of sites the limit is instead the matrix supersymmetric pair of Pöschl–Teller Hamiltonians, with the reflection in the Hamiltonian as a Dunkl potential. The bispectral pair of the discrete model goes over to the Dunkl realization of the anticommutator algebra of the little $-1$ Jacobi polynomials.

We consider generalised Casimir operators associated with irreducible representations for finite-dimensional, semisimple, complex Lie algebras. Using Klimyk's formula and Okubo's formula, we derive a rational ordinary generating function for the Casimir operator eigenvalues for simple Lie algebras which does not rely on knowledge of the Weyl group. We also directly reproduce generating function expressions for the classical Lie algebras originally due to Perelomov and Popov.

The para–Racah Hamiltonian is the Jacobi matrix $J$ of the para–Racah polynomials, taken at the persymmetric value of their isospectral deformation parameter: a mirror–symmetric tridiagonal matrix whose spectrum is a quadratic bi–lattice, the squares $(s+a)^2$ and $(s+c)^2$ of two interlaced arithmetic progressions of unit step, and which carries out fractional revival when the grid parameters $a$ and $c$ satisfy a Diophantine condition. As the number of sites grows, $J$ contracts to the trigonometric Pöschl–Teller Hamiltonian on a finite interval. We establish this twice: on the recurrence relation, and on the polynomials themselves, for which it is the confluence Racah~$\to$~Jacobi. Both grid parameters survive: $c-a$ becomes a scalar centrifugal term at the interior point fixed by the mirror — the exact analogue of the para–Krawtchouk origin — while $a+c-1$ sets two confining singular walls at the ends, which turn the spectrum into the quadratic pair $(t+a)^2,(t+c)^2$ with Jacobi eigenfunctions $P_t^{(a+c-1,\,\mp(c-a))}$, the two sublattices converging to two families that differ only by the sign of the second Jacobi parameter. The mirror symmetry of $J$ becomes the reflection about the midpoint; it grades the two families and carries the fractional revival to the limit. Algebraically, the bispectral pair of the discrete model, which generates the Racah algebra $\mathfrak{r}_1$, becomes a pair of generators of the Jacobi algebra $\mathfrak{j}_1$ — one rung down the Askey–Wilson hierarchy — the dynamical algebra of the Pöschl–Teller Hamiltonian.

The para–Krawtchouk polynomials define, through their recurrence relation, a finite oscillator model that carries out fractional revival. We determine the continuum limit of this Hamiltonian, both from its recurrence relation and directly on the polynomials. The limit is the singular (isotonic) oscillator $-\partial_\eta^2+\eta^2+\tfrac{\gamma^2-1}{4\eta^2}$: the deformation parameter $\gamma$ survives as a central impurity of the couplings, decaying as the inverse square of the distance — the marginal rate, which the continuum scaling turns into a finite $1/\eta^2$ term. Using the decomposition of the para–Krawtchouk polynomials into two Hahn families, we prove the confluence, for every degree, of the eigenvectors on the two sublattices to generalized Hermite functions. The bispectral pair of the discrete model, which generates the Hahn algebra, becomes in the limit a pair of generators of $\su$, the dynamical algebra of the singular oscillator.

We give a polynomial-time algorithm to sample from the Gibbs measure of the Sherrington-Kirkpatrick (SK) model with $o_n(1)$ error in total-variation distance (TVD) at any inverse-temperature $\beta < 1$. The algorithm combines algorithmic stochastic localization (ASL) with rejection sampling over path-space via Jarzynski's equality (JE). The analysis extends the authors' prior $\beta < 1/2$ result [arXiv:2605.03718] by replacing all global regularity requirements in the stochastic differential equation (SDE) error analysis with local regularity around likely trajectories. The relaxed regularity is established using Celentano's proof of the local strong convexity of the TAP free energy [arXiv:2208.09550]. The analysis utilizes the cavity interpolation theory and free probability toolkit developed in the authors' previous result, where the former applies nearly verbatim and the latter applies supplemented with Lipschitz and $C^2$ extensions of various functions. The ASL and JE analysis arises from using the TAP free energy as an efficiently computable proxy for the actual free energy of the stochastically localized Gibbs measure [$ \S $ 3, arXiv:2605.03718]. We give a list of \vocab{desiderata} encapsulating the approximation and regularity properties required of the TAP free energy, relaxing those of [$ \S $ 2.5, arXiv:2605.03718] to only require local regularity. These generic desiderata are potentially applicable in other settings where a free energy surrogate exists, giving algorithmic sampling guarantees while bypassing the usual functional inequalities based approach.

The purpose of this article is to explicitly describe how the integration technique known as the change of variables is used in graduate-level Probability and Statistics courses. Throughout the article, we employ this integration technique seven times across various examples and results. We demonstrate that four applications of this technique are required to derive the probability density function (PDF) of the Student's t-distribution.

The structure of divergences in several nonpolynomial scalar quantum field theories are investigated. Various $n$-point Green functions for nonpolynomial theories are calculated and their renormalization is analyzed, with special regard for the nonpolynomial $\left|\phi\right|^{3}$ theory in $4$ spacetime dimensions.

We study critical Strichartz estimates for the resonant hyperbolic Schrödinger flow on the waveguide $\mathbb{R} \times \mathbb{T}^2$. For initial data localized to frequencies at most $N$, we prove that the unit-time $L^2\to L^{10/3}$ operator norm has sharp order $N^{1/5}$. We further track its dependence on the observation length, obtaining sharp estimates in the short- and long-time regimes and characterizing the frequency-dependent intervals on which the critical estimate remains lossless. The precise dependence on $T$ in the intermediate regime $N^{-1}<T<1$ remains open. The upper bounds are obtained by lifting the global endpoint estimate on the elliptic cylinder $\mathbb{R} \times \mathbb{T}$ to Hilbert-valued data, preserving its time summability, and combining this with periodic Bernstein and interpolation. The lower bounds exploit an exact integer null direction in the periodic frequency lattice together with a suitably scaled Euclidean wave packet. Finally, using a Bloch–Floquet transference argument, we recover the known elliptic unit-time estimate from compact torus theory and obtain a lossless elliptic estimate on frequency-dependent short time intervals.

This review surveys diverse types of random motion within a unified framework, introducing essential mathematical tools, including special functions and fractional calculus, in a pedagogical manner. We focus on generalized random motions in one-dimensional space described by Fokker-Planck like equations, whose fundamental solutions take a universal and simple form in the Laplace-Fourier domain in terms of a generic function $\varphi(s)$. By exploring the functional forms of $\varphi(s)$ consistent with physical constraints, we derive a broad class of stochastic dynamics. This includes Brownian and active motion, simple and anomalous diffusion, and processes featuring stochastic resetting or trapping mechanisms. For each case, we provide the corresponding kinetic equation, explicit solutions (where available), and mean-square displacements alongside their asymptotic behaviors. Finally, we discuss the physical interpretation of these motions through the lens of time-changed processes and subordinators.

We study planar openbooks obtained by lifting braids through branched covers of D^2, together with the quantum operations in the Ising representation. We give a criterion for the Stein fillings to be Calabi-Yau (CY). Among positive factorizations of a fixed monodromy, a CY one has minimal length and its filling minimizes \chi and b_2. Applied to Baykur's recent examples, this gives a planar contact 3-manifold with infinitely many non-homeomorphic CY fillings, the cover there having degree \ge 6. At degree 4 a single CY filling forces every filling to be CY, and in one case the filling is unique; at degree \le 3 it is unique, and CY under a mild condition. Admissible cuts of D^2 decompose the openbook into subopenbooks. Every positive factorization then localizes, and the CY condition holds exactly when it holds locally. The state space is then a direct sum of tensor products indexed by the compatible parity choices, with at most two qubits per factor when the pieces have degree $\le 4$, the same range in which the CY condition depends only on the monodromy. For one degree 4 cover, the points, lines and flags of the two-qubit doily are realized by its subopenbooks. Fifteen liftable braids there share the same quantum operation and Stein filling, and are separated only by the subopenbook systems they admit. A choice of tensor-product structure is thus carried by the lift, not by the braid group representation, suggesting a link between contact topology and quantum entanglement.

The most general strategy for estimating parameters of a quantum channel is to probe the channel many times sequentially with the help of coherent ancillary systems. However, finding the strategy that minimizes the mean-squared error (MSE) is a notoriously hard problem, requiring optimization over the initial probe state, intermediate evolutions, final measurement, and classical postprocessing. In this work, we study how the minimum attainable MSE can be bounded from below in this setting. We propose a lower bound for the MSE of any physically implementable sequential sensing strategy that applies to estimating parameters of any physical quantum evolution, including quantum channels probed over multiple time steps under temporally correlated noise and general temporally correlated quantum processes. We show that this lower bound can be formulated as a semidefinite program (SDP). Our bound generalizes the Nagaoka-Hayashi lower bound (that applies to parameter estimation tasks with separable measurements) to the regime where the quantum channel can be accessed sequentially multiple times. For regular finite-dimensional single-parameter sequential estimation tasks admitting an optimal strategy, we show that our SDP bound is tight and provide a procedure to obtain a locally optimal strategy from the SDP optimizer. Using our optimization bound, we compute the MSE lower bound and derive new protocols for several other physically motivated problems including error-corrected quantum sensing under an unknown correlated noise model, multiparameter estimation, and estimation with qudit systems.

Commuting-projector Hamiltonians with integer spectrum such as the 1+1d Ising ferromagnet and the 2+1d toric code may be regarded as the generators of non-on-site $U(1)$ symmetries. In this paper we study whether these symmetries can be made on-site in the presence of finite and infinite-dimensional ancillae. In the finite-dimensional case, we prove that the fractionalized excitations, such as domain walls and anyons, make it impossible to on-site these Hamiltonians. On the other hand, we show that when infinite dimensional ancillae (e.g., rotors) are allowed, both the 1+1d Ising ferromagnet and the 2+1d toric code can be disentangled, yielding on-site symmetries. This is not in contradiction with the non-trivial ground state order of these Hamiltonians, since the ancilla Hamiltonians are no longer bounded from below. Our results indicate agreement between lattice and quantum field theory anomalies in the presence of infinite-dimensional ancillae, while with finite-dimensional ancillae the obstruction to on-siteability may even be non-invertible.

We consider a fixed finite number of independent two-dimensional discrete Gaussian free fields on a box of side $n$, with zero boundary values, conditioned so that the vector-valued field stays outside a fixed ball on a macroscopic interior domain. We prove that the field, divided by $\log n$, converges in distribution in $L^2$ to twice the equilibrium potential of the domain multiplied by a uniform random unit vector. We also identify the exact capacity cost of the conditioning event and obtain alignment on a deterministic density-one set of sites, with no restriction on their mutual distance. The proof combines the norm repulsion theorem of Korzhenkova and Sepúlveda with Gaussian large deviations and the equality case in the capacity variational problem.

We describe a family of convex uniformly continuous functions $E^\lambda_F$, $\lambda>0$, on the set of states of a bipartite quantum system (consisting of finite-dimensional or infinite-dimensional subsystems), which monotonically increase and converge pointwise to the Entanglement of Formation (EoF) as $\lambda\to0$. These functions are "nonselective" entanglement monotones defined in a way close to the construction of the Moreau-Yosida regularization (the Moreau envelope) of a convex function on a convex set used in the modern convex analysis. So, we call the functions $E^\lambda_F$ the Moreau-Yosida approximations of the EoF and describe their equivalent definitions and basic properties. The semicontinuity bounds for the EoF (presented in [Lob.J.Math., 46(6), 2632-2658]) allow us to obtain easily computable upper bounds on the difference $E_F(\rho)-E^\lambda_F(\rho)$ for a given state $\rho$. These bounds give easily computable bounds on the rate of uniform convergence of the function $E^\lambda_F$ to the EoF as $\lambda\to0^+$ on the sets of states with bounded rank/energy of one of the marginal states. We also discuss sufficient conditions for the coincidence of $E_F(\rho)$ and $E^\lambda_F(\rho)$ at a given state $\rho$ for all $\lambda$ small enough and consider several classes of states for which such coincidence takes place. The conjectured selective LOCC-monotonicity of the functions $E^\lambda_F$ and a possible way to prove it are briefly discussed. Secondary parts of the article are written with the help of ChatGPT-5.6.

The Anderson tight-binding model is a fundamental model of quantum transport and localization in disordered media. Completing a problem left open by Bourgain and Kenig, this paper proves Anderson localization near the bottom of the spectrum for the lattice Anderson model with Bernoulli potential, in any dimension $d\ge 2$. The proof uses the multiscale framework of Fröhlich-Spencer and Bourgain-Kenig, and the main new ingredient is a probabilistic discrete unique continuation principle (PDUC) for the discrete Schrödinger equation. This PDUC is established via a bootstrap argument and a key probabilistic lemma proved by adaptively revealing the random potential.

Non-Abelian quantum codes connect quantum error correction, phases of matter, and computational resources. In this work, we develop a general framework for constructing non-Abelian quantum low-density parity-check (qLDPC) codes by gauging sheaf codes via cup products and use it to obtain families with constant encoding rate and linear distance. We resolve the coupled logical constraints using explicit representatives to characterize the full gauged code space. We provide a fundamental treatment of code distance based on the general Knill–Laflamme condition and combine expansion with cleaning to establish protection against arbitrary low-weight errors. We further construct an almost-good family whose entire code space exhibits long-range magic. Gauging and ungauging also enable logical Clifford measurements that prepare encoded magic states. These results extend good qLDPC codes beyond the Pauli stabilizer setting and provide a concrete foundation for exploring non-Abelian phases beyond geometric locality and pursuing the no low-energy trivial magic conjecture.

We determine the exact universal high-temperature separability threshold for Pauli Hamiltonians of bounded degree $\Delta\ge2$. If every coefficient in the Hamiltonian has magnitude at most one and each term has overlapping support with at most $\Delta$ other terms, the Gibbs state is a mixture of product Pauli eigenstates whenever \[ \beta \le z_\Delta := \operatorname{arctanh}\left[\max_{0\le x \le 1}x\left(\frac{1-x}{1+x}\right)^{\Delta-1} \right]. \] For every $\beta>z_\Delta$, a finite commuting Hamiltonian with maximum overlap degree at most $\Delta$ has an entangled Gibbs state. At any fixed $\beta<z_\Delta$ strictly below the threshold, a classical polynomial-time algorithm produces samples from a distribution over product Pauli eigenstates approximating the Gibbs state in trace distance.

We explicitly construct the additional integrals of motion, ensuring maximal superintegrability of the discrete-time rational Ruijsenaars-Schneider model. Using them, we investigate the algebraic aspects of superintegrability in both continuous- and discrete-time settings. In particular, we determine the complete structures of the polynomial symmetry algebras associated with both the rational Ruijsenaars-Schneider model and its discretization. We demonstrate that discretization leads to a nontrivial deformation of the continuous symmetry algebra with respect to the discretization parameter, thereby extending recent analogous results from the rational Calogero-Moser system to its relativistic generalization.

The classical Becker-Döring equations describe the formation of clusters by aggregation and fragmentation of monomers. If the total amount of mass is supercritical, larger and larger clusters are formed, leading to an asymptotic loss of mass in the long time limit. The two-component Becker-Döring system arises when clusters are built from two types of monomers. Here, we study the natural extension of the one-component system, where no energy or entropy is leaving or entering the system - the so called detailed balance assumption. We rigorously prove that all initial conditions admit a solution minimising the relative entropy as time approaches infinity. This shows, that the long time limit in the weak* topology is selected through the initial Type I and Type II masses. Furthermore, the proportion of lost Type I mass is determined by the limit point via the mixing ratio of Type I and Type II monomers that maximises the binding energy. The main difficulty of the two-component system is that no relative entropy is weak* continuous, which is crucial for the classical argument [2]. Instead, our proof is based on a discrete two-dimensional logarithmic Sobolev inequality, that bounds the entropy dissipation on the correct timescale. Our approach also improves current assumptions to determine the long time behaviour for the one-component system.

The two-component Becker-Döring model is a two-dimensional coagulation-fragmentation equation that describes the dynamics of clusters built from two different types of monomers growing and shrinking through monomer interaction. We show the existence and mass conservation of solutions and give a partial result regarding uniqueness. Furthermore, we prove under the assumption of detailed balance the entropy equation, followed by a detailed study of the equilibrium states and relative entropies. We can determine the minimum of all relative entropies under sequences of fixed Type I and Type II masses and fully describe the minimising sequences in the weak* topology.

For independent bond percolation on $\mathbb{Z}^d$ with parameter $p$, let $p_c^b(d)$ be the critical probability. We prove that for each $d \geq 9$, there is $\epsilon_d>0$ such that for each $p \in (p_c^b(d), p_c^b(d)+\epsilon_d)$, the complement of the infinite open cluster stochastically dominates a supercritical site percolation on $\mathbb{Z}^d$. This improves the previous results by Grimmett, Holroyd and Kozma 2014, and Bock, Damron, Newman and Sidoravicius 2020. Numerical estimates of $p_c^b(d)$ and $p_c^s(d)$ (the site critical probability) suggest that a similar stochastic domination result holds for all $d \geq 4$. For the Ising model on $\mathbb{Z}^d$ with inverse temperature $\beta$, let $\beta_c(d)$ be the critical inverse temperature. We prove that for each $d \geq 8$, there is $\epsilon_d>0$ such that for each $\beta \in [0,\beta_c(d)+\epsilon_d)$, the $+$ spins under the minus phase stochastically dominate a supercritical site percolation on $\mathbb{Z}^d$. This improves the previous result of Aizenman, Bricmont and Lebowitz 1987. Numerical estimates of $\beta_c(d)$ and $p_c^s(d)$ suggest that a similar stochastic domination result holds for all $d \geq 5$.

We study the $L^p$-boundedness of Isozaki–Kitada type modified wave operators associated with Schrödinger operators $P=-\partial_x^2 + V$ in one dimension for potentials $V$ including long-range potentials such as negative Coulomb-like ones $V(x) = -(1+|x|^2)^{-\mu/2}$ for $\mu>0$. For such potentials, if $ \mu \in (1, 2)$, we prove that the middle and the high energy parts are bounded on $L^p (\mathbb{R})$ for any $p \in [1, \infty]$ but the low energy part is unbounded except for $p=2$. Moreover, if $\mu \in (0, 1]$, we show that even the middle energy part is unbounded except for $p=2$. Actually, regarding the $L^p$-boundedness in the middle and the high energy regimes, we give a complete classification of slowly decaying potentials without imposing any sign condition.