Papers

In this article, we study a generalisation of Rényi's car-parking problem in which isosceles trapeziums are sequentially deposited on a strip. We derive an explicit formula for the parking constant in terms of the lengths of the two bases, recovering the classical rectangular and triangular models as special cases. We further obtain quantitative finite-size asymptotics for both the expected number and the variance of deposited trapeziums, with convergence rates that depend explicitly on the geometry of the deposited particle. In particular, although the recursive construction involves two different substrate geometries, their variances have the same leading asymptotic density. Finally, we show that the parking constant depends non-monotonically on the ratio of the two base lengths and possesses a unique minimiser, so that the least efficient shape is a genuine trapezium rather than a triangle.

We study three identical bosons with contact attraction of strength mu on the square lattice at total quasimomentum K = 0 and K = pi. First, at K = 0, for all sufficiently large mu there are exactly two eigenvalues below the essential spectrum. The trimer branch follows from Kato perturbation theory at the simple isolated eigenvalue 3 of the pair-contact operator, uniformly in K; the atom-dimer branch follows from a Feshbach-Schur reduction whose limiting equation coincides with the second Birman-Schwinger (Faddeev) eigenvalue of the finite-rank principal part. The spectral gap between the two branches grows linearly in mu. Second, at K = pi the trimer branch is unchanged up to O(mu^-2), the two-particle threshold equals -mu + 4 exactly, and the gap to it is 2mu - 2 plus lower order corrections. The second bound state disappears because the exchange operator changes sign, S_pi = -S_0: the s-channel eigenvalue that exceeds 1 at K = 0 becomes negative, while the only positive atom-dimer channel at K = pi is the odd p-channel, bounded by 2 - 4/pi < 1. Consequently the operator has exactly one eigenvalue below the threshold minus 36/(mu - 10). Third, the ground state is exponentially localized in the relative coordinates. The decay exponent is positive unconditionally; under a non-cancellation hypothesis it is bounded above by log(3mu) plus lower order terms, while finite-volume computations indicate the slowest tail is the atom-dimer dissociation channel with decay rate one half log(2mu). A dictionary between the principal-part method, Kato theory, and the atom-dimer threshold matrix is given, together with the three-dimensional constants for comparison.

We study an infinite family of strongly coupled four-dimensional $\mathcal N=2$ superconformal field theories (SCFTs) distinguished by a trivial Higgs branch, known as the $(A_2,D_{3m+1})$ Argyre-Douglas theories. We propose that their associated vertex operator algebras (VOAs) are the doublet algebras $\mathcal A(4m+2)$, an infinite family of non-rational vertex operator superalgebras with a remarkably simple strong generating set of only three fields. We provide several highly nontrivial checks of this proposal. In particular, we reproduce the four-dimensional conformal anomalies $a$ and $c$ from the VOA and analytically prove the exact equality between the Schur index of the SCFT and the supercharacter of the VOA. A key ingredient is a diagonal-gauging realization of the $(A_2,D_{3m+1})$ theories in terms of two simpler building blocks, which makes the Schur-index computation tractable. Remarkably, we find that the resulting Schur index of the $(A_2,D_{3m+1})$ theory coincides with that of $\mathcal N=4$ $SU(2)$ Super-Yang-Mills theory, up to an overall prefactor and an appropriate identification of fugacities. We also discuss a generalization to the two-parameter family $(A_{2s},D_{(2s+1)m+1})$, whose members likewise have trivial Higgs branches and admit diagonal-gauging realizations. A particularly interesting subfamily is $(A_{2s},D_{2s+2})$, for which the four-dimensional conformal anomalies coincide, $a=c$. Our results reveal a systematic connection between Higgsless SCFTs, diagonal gauging, and strongly finite but non-rational vertex operator algebras.

In positive degree, equivariant Gromov-Witten invariants of Calabi-Yau fivefolds are expected to admit an interpretation in terms of M2-branes, while we conjecture that constant maps are governed by the corresponding supergravity index. We make the first expectation precise by proposing a modular interpretation of M2-branes supported on several smooth curves meeting at an $n$-fold point. Together, the two pictures yield conjectural formulas for the Gromov-Witten invariants, which we translate into closed formulas for pointed and unpointed quintuple Hodge integrals. We show that these conjectures imply a K-theoretic Gromov-Witten/Pairs correspondence for local curves in degree one and link the generating series of constant maps to Donaldson-Thomas theory of points. We also prove the conjectures in two limits of the equivariant parameters. Along the way, we obtain new closed formulas for certain triple Hodge integrals.

Metastability is typically formulated for systems with finitely many metastable states, most often in positively recurrent settings. In this article, we extend the resolvent framework for metastability to Markov processes with countably many metastable states, thereby encompassing null-recurrent and transient dynamics. For a family of processes on locally compact Polish spaces, we prove, under a mild boundary regularity assumption, that the asymptotic flatness of microscopic resolvent solutions, supplemented by two compactness conditions, is equivalent to convergence in law of the projected trace processes to a limiting Markov chain and to the negligibility of the time spent outside the metastable sets. We apply this framework to two non-compact stochastic systems. First, we study a condensing inclusion process on a countably infinite, uniformly locally finite graph, in a setting where the process may be null recurrent or transient. We prove that the condensate location converges to a weighted random walk on the underlying graph, while the time spent away from the fully condensed configurations is negligible. Second, we study a small-noise one-dimensional Langevin diffusion with countably many stable equilibria, without assuming ergodicity. In the Eyring-Kramers time scale associated with the minimal energy barrier, we establish local equilibration inside each well, convergence of the well-index process to an explicit nearest-neighbor Markov chain on $\mathbb Z$, and negligibility of inter-well excursions. Together, these results broaden the scope of resolvent-based metastability theory beyond finite metastable state spaces and positive recurrence.

In this work we prove a Payne-Weinberger type inequality for quantum dot Dirac operators defined on bounded and simply connected planar domains. This is a sharp upper bound for their first positive eigenvalue depending only on the isoperimetric deficit of the domain. To this end, we use a recently studied connection with the so-called $\overline\partial$-Robin Laplacian and we establish an analogous inequality for its first eigenvalue, relying on the corresponding inequality for the Robin Laplacian.

We consider eigenvalue counting functions of Robin Laplace operators on cuboids where the Robin parameter and the spectral cut-off $\lambda$ are coupled. Our main focus is the critical regime, in which the Robin parameter is proportional to $\sqrt{\lambda}$. We obtain a two-term asymptotic expansion for the counting functions in this coupled setting. In the critical regime, the second term in the asymptotic expansion depends non-trivially on the proportionality constant and interpolates continuously between the corresponding second terms for the Dirichlet and Neumann Laplacians. Outside the critical regime, one recovers the corresponding Dirichlet or Neumann asymptotics. We also establish a Pólya-type inequality for the counting function whenever the ratio of the Robin parameter and $\sqrt{\lambda}$ exceeds a dimension-dependent threshold.

Heat transfer in a dense gas between two parallel plates is studied for assessing the quantitative difference of two variants of the Enskog equation, i.e., the original Enskog equation and the Enskog equation with a modified Enskog factor recently proposed in Phys. Rev. E 111, 065108 (2025). The main advantage of the latter is that the H theorem, which the former lacks, has been established. The influence of the modification is assessed by comparing the macroscopic quantities in detail. It is found that the impact of this modification is quite limited and the two variants reasonably agree as a whole.

This paper studies the category $\mathcal{C}$ of grading-restricted $C_1$-cofinite generalized modules for the vertex operator algebra associated to the $W$-algebra $W(2,2)$. We first show that $\mathcal{C}$ is the same as the category of finite length modules whose simple composition factors are the irreducible highest weight $W(2,2)$–modules $L[r]$ of highest weight $(\frac{1-r}{2}, 0)$ for $r \in \mathbb{Z}_{> 0}$, and hence $\mathcal{C}$ carries a braided tensor category structure. Then we prove the fusion rules for the simple objects $L[r]$ are governed by the $\mathfrak{sl}_2$ Clebsch–Gordan rule. In particular, we prove \[ L[r]\boxtimes L[s]\cong \bigoplus_{i=0}^{\min\{r,s\}-1} L[r+s-1-2i]. \] Using the fusion rules and a recent result of Etingof–Penneys, we establish the rigidity of $\mathcal{C}$. We also show the semisimple subcategory generated by the simple objects $L[r]$ is tensor equivalent to the category Rep $\mathfrak{sl}_2$ of finite dimensional $\mathfrak{sl}_2$-modules.

We prove the existence and uniqueness of a local-in-time solution to the stochastic Cahn-Hilliard equation in space dimensions 1, 2, 3, 4 and 5. Our method relies on the Da Prato-Debussche trick for $d=4$ and paracontrolled distributions for $d=5.$

The perihelion of comet 1P/Halley in 2061 is an excellent chance to revisit this object of outstanding scientific and cultural relevance. During its 1986 approach to the Sun, it was targeted by several flyby missions. Due to its retrograde, highly-inclined orbit, the relative velocities during the encounter were large, limiting the scientific return due to the short time spent inside the coma. A rendezvous trajectory would overcome this limitation, but the design is challenging due to the limitations of current propulsion technology. Given the lead times of spacecraft development and the long duration of the interplanetary transfer to the comet, it is imperative to start mission planning as soon as possible. We present a novel rendezvous strategy, combining unpowered Jupiter and Saturn gravity-assists with deep-space low-thrust arcs. It minimizes launch energy and propellant budget, constraining the arrival to occur before the onset of high activity. The double flyby strategy reduces the otherwise prohibitive cost of the plane change maneuver. Closed expressions for the optimal flyby geometry, together with an explicit low-thrust transcription technique, reduce the number of design parameter to three, improving computational efficiency. Crucially, this is the first rendezvous mission concept achievable with well-proven technology (standard radioisotope thermoelectric generators and a Hall-effect thruster) and compatible with several existing launchers. We describe the trajectory optimization strategy and perform a comprehensive exploration of the design space. Finally, we present two promising proof-of-concept trajectories in detail.

We investigate inverse scattering at a fixed energy for semiclassical Schrödinger operators with smooth potentials. For compactly supported potentials satisfying a natural virial condition, we prove that $ \liminf_{h\to0} \|S_1(\lambda,h)-S_2(\lambda,h)\|<\sqrt{2} $ implies $V_1=V_2$, provided the fixed energy $\lambda$ lies above both potentials. We then consider radial short-range repulsive potentials. Under a monotonicity assumption on the radial force, we show that the semiclassical differential cross section at a single fixed energy, up to $o(1)$ as $h\to0$, determines the potential throughout the classically accessible region. Finally, we obtain an analogous rigidity result for nontrapping compactly supported perturbations of the Euclidean metric under a strict convexity assumption.

Consider the polynomial $f=x^N+y^N-z^N$ where $x,\,y$ and $z$ are positive integers and $N \ge 3$ is an integer. By Fermat's Last Theorem, $f$ is never zero so that its reciprocal, $1/f$, has no singularities. We therefore study the finite sum of the reciprocal: $S(m,N)=\sum_{x=1}^m\sum_{y=1}^m\sum_{z=1}^{m}\frac{1}{f}$. The terms $1/f$ can be positive, negative and their magnitude is less than unity. A key observation is that $S(m,N)$ can be split into two convenient parts: a dominant contribution $D(m,N)$ that has a simple analytical expression and a remainder $R(m,N)$ which is more complicated but negligible compared to $D(m,N)$. Therefore, $S(m,N)$ is almost identical to $D(m,N)$. The analytical expression for $D(m,N)$ is $(2\,m-1)\,H_m^{(N)}$ where $H_m^{(N)}=\sum_{x=1}^m\frac{1}{x^N}$ approaches quickly the Riemann zeta function $\zeta(N)$ as $m$ increases. Therefore, the original sum $S(m,N)$ has a simple expression: it is basically linear in $m$ with slope equal to $2\,\zeta(N)$. Its linear behavior is not an asymptotic result; plots of $S(m,N)$ vs. $m$ for different $N$ show a straight line starting at $m=1$. $S(m,N)$ deviates slightly from a straight line over a small interval $8\le m\le 12$ for the case $N=3$. This slight deviation is due to Fermat near misses where $x^3+y^3-z^3=\pm 1$ (for $z\ne x$ and $z\ne y$); these create a jump in the remainder $R(m,3)$ at $m=9$. We make a numerical and analytical study of the remainder $R(m,N)$. From the numerical analysis, $R(m,N)$ converges for $N\ge 4$ but it was harder to tell whether $N=3$ converged. An analytical study based on a comparison of $R(m,N)$ to its Cauchy principal value integral, shows that $R(m,3)$ likely diverges logarithmically. It also shows that $R(m,N)$ converges for $N\ge 4$ in agreement with the numerical analysis. We discuss in the conclusion some interesting questions for future investigation.

We establish an identity relating the spatial derivative of a two-point free-energy correlation to two fundamental geometric quantities: the annealed probability that two independently sampled polymers, in the same environment, meet before reaching the prescribed terminal level, and the exit-point location of a single polymer. Our result is inspired in part by recent integration-by-parts work of Gu and Quastel for the KPZ equation. A crucial step in their continuous setting relies on an ingenious application of Itô's formula, which has no direct analogue in our semi-discrete setting. Instead, we exploit intrinsic symmetries of the polymer model together with the memoryless property. We also provide two applications: (i) we recover the horizontal Burke property by proving that the anchored stationary horizontal free-energy profile is a two-sided Brownian motion; (ii) in the zero-temperature limit, we obtain the corresponding identity for Brownian last-passage percolation.

The exact bound-state wavefunctions and the corresponding energy equation are calculated for a new generalized harmonic oscillator problem describing a spin-1/2 fermion in 3+1-dimensions, involving scalar, vector, and tensor couplings acting simultaneously within a particular plane of motion. For the scalar and vector coupling, singular harmonic oscillator shapes are considered, such that the singular term is needed to allow analytical solutions for the wavefunctions while preserving binding under adequate conditions. For the tensor sector,the Dirac oscillator potential is employed, which adds another independent binding mechanism to the problem. The exact bound-state solutions are computed by specifically tuning coefficients for an appropriate pair of Ansätze for the radial functions, which leads to wavefunctions in terms of generalized Laguerre polynomials. Although the energy equation cannot provide a general expression for the energy spectrum, specific constraints on the quantum numbers as simple functions of the external potential parameters can be derived for it, which determines the conditions for bound solutions to exist, and of what type: particle, antiparticle or both. It is shown that the results can be simply mapped to the spherically symmetric analogue problem, and this is used to show that the general result encompass several previous particular cases of spherically symmetric harmonic oscillator problems in the Dirac equation available in the literature.

Continuous-time quantum walks on a lattice spread ballistically and converge to a limiting distribution for the rescaled position. On the hard-wall half line the boundary reflects the walker but does not change the bulk dispersion, so the ballistic front remains set by the maximal group velocity. We ask how close one can get to this front when the initial state is constrained to occupy only the first $M$ sites. We show that optimizing the moment-generating function of the limiting-velocity distribution over this finite-support class reduces to the principal eigenvalue of an explicit $M\times M$ Hermitian matrix, which yields the optimal state by direct diagonalization. For large $M$, a near-front scaling limit produces a continuum description that controls the entire peak region. In particular, the maximal mean drift approaches the front with an inverse-square deficit in $M$ and a constant prefactor $\pi^2/4$. Numerical results verify both the finite-$M$ spectral formulation and the predicted scaling.

The idea that information can be localized is pervasive in physics. In thermodynamics, entropy flows from one reservoir to another, and in quantum gravity, the entropy content of some spatial regions is bounded. In quantum information theory, information is transmitted from one place to another, and entanglement can provide an advantage when there are constraints on such transmission. Shannon entropy quantifies localizable information in the sense that marginalization defines an outer measure on the set of parts of a classical multipartite system. In contrast, von Neumann entropy does not quantify localizable information in this sense. However, a variant of von Neumann entropy, which originates in noncommutative geometry, does. We prove that this adjusted von Neumann entropy is the minimum quantum entropy that quantifies localizable information. We also prove that this is the unique quantum entropy that characterizes Bell states in terms of redundant information, directly generalizing the classical case. We work with finite-dimensional $C^*$-algebras throughout, modeling finite quantum systems that may have superselection sectors.

We define the notion of superintegrability of a Hamiltonian system on a stratified symplectic space. We focus on spin Calogero-Moser-Sutherland (sCMS) systems, where the phase space is a stratified symplectic space obtained by the Hamiltonian reduction of the cotangent bundle over a compact Lie group, and demonstrate that the sCMS systems for $SU(3)$ are superintegrable.

Metadynamics is a class of enhanced-sampling methods that is widely used in molecular modeling. Here, we consider metadynamics with a one-dimensional collective variable and explore the limitations of using Gaussian hills in light of existing convergence results. We reduce the corresponding evolution problem to a replicator-type differential equation and analyze its long-time behavior. For a periodic collective variable, we prove the convergence of metadynamics. However, the situation is fundamentally different when the collective variable is defined on a finite interval. In this case, the stationary solution only exists in the weak sense as a finite atomic measure. The solution to the differential equation weakly converges to it. Consequently, metadynamics is ineffective due to the absence of a clear quasi-stationary state. A quasi-stationary, transient solution can only be captured in the "INTERVAL" framework if the free energy outside the finite interval remains nearly constant across sufficiently large regions compared to $\sigma$.

In this paper, we present a model for classical and quantum reference frames for $C^*$-algebras by utilizing Hilbert modules and equivariant correspondences of $C^*$-dynamical systems. We show the classical model provides a commutative baseline for the quantum model, and we discuss how the quantum model encodes the Hilbert space operational formulation of quantum reference frames. We also include various explicit examples for both classical and quantum reference frames.

We investigate how a fractional oscillator reacts to external forces when the dynamic law includes a damping term that is explicitly fractional. We show that the corresponding three-parametric fractional differential equation (in the Caputo sense) admits formal solution. Two of the fractional parameters are associated with the intrinsic dissipation mechanism derived from the fractional nature of the oscillator. The third fractional parameter characterizes the damping term against which the oscillator is driven by the external force. In the absence of fractional damping, the two-parametric fractional dynamic law can be associated with the Newton-Scott-Blair model for viscoelastic materials, in which position and external force are linked to the material's strain and stress, respectively. In general, the Laplace transform of the external force defines the profile of the three-parametric solution; specific examples include the absence of external forces as well as constant, sinusoidal and stepped external forces. In the case where the three fractional parameters converge to the same value, it is observed that, under the action of a constant driving force, the system reaches a fixed position after a transient period. With a sinusoidal driving force, the system exhibits persistent oscillatory behavior. These results suggest that the external driving force strongly influences the long-term behavior of the fractionally damped oscillator.

A collection and clarification is given of the key mathematical details of the method of applying the convolution theorem as part of the solution to the two-dimensional magnetic inverse problem: reconstructing current density from magnetic field data. Particular emphasis is given to obtaining a closed-form solution for the Fourier transform of the convolution kernel that features in the Biot-Savart law. Having done this, the Fourier-transformed current density is obtained, via the convolution theorem, as a function of Fourier-transformed magnetic field data and other variables.

We give a variational derivation of the central configuration equations in pair space, where the relative position vectors between pairs of bodies serve as the primary variables. Vector Lagrange multipliers enforce the linear triangle relations required for these pair vectors to be realizable in the plane. For arbitrary $N$, we show that the constrained Hessian decomposes as $H_{\mathcal C}=L^\Delta+\widetilde L$, where the gap Laplacian $L^\Delta$ is positive semidefinite and the transverse Laplacian $\widetilde L$ is signed and contains all possible negative directions. For non-collinear planar four-body central configurations, the signed part has rank two. This yields a new proof of the known bound that the Morse index is at most two. We then characterize positive definiteness and degeneracy of the Hessian by the eigenvalues of an explicit $2\times2$ matrix or, equivalently, by two generalized eigenvalues of a pencil of symmetric matrices. We illustrate the general criterion using the equal-mass square central configuration. For reflection-symmetric configurations, the generalized eigenvalue problem decomposes into two independent smaller problems. This reduction applies to kite and isosceles trapezoidal configurations. For the trapezoid, however, the reflection acts non-orthogonally on the pair coordinates, giving a different, non-orthogonal reduction. Applying the reflection-symmetry reduction to the rhombus family, we prove that the Hessian is positive definite modulo the rotational zero mode; equivalently, every rhombus central configuration is nondegenerate modulo rotations and has Morse index zero.

Three seemingly distinct extended-body propagators devised by Ansoldi et al., Erbin et al., and Stanford et al.\ respectively arise in completely different settings, though each can be considered a realization of the Polyakov path integral. This paper proves that the Ansoldi propagator corresponds to the characteristic function of the random variable with the probability density function $f(A)=C\,A\,(\alpha^{2}+A^{2})^{-D/2}\mathbf{1}_{A>0}$ in ambient dimension $D$, of scale $\alpha=\sqrt{b/2}$. The exponent is fixed by the dimension and by nothing else in the Ansoldi/geometric comparison; the on-shell two-point amplitude of Erbin et al. and JT gravity agree with its $D=4$ member in the senses made precise below. After an extensive literature review, the author believes that this observation and constructed area measure is novel. The area measure will, however, emerge as the natural candidate for the extension of the earlier non-extended body path space measures to two dimensional simplicial complexes. We further show that the area measure of this paper is exactly the law of the area of a triangle, two of whose sides are drawn from the measure obtained in the companion manuscripts for the one-dimensional Polyakov integral. We close with an argument, rigorous in part, which reduces the measure of tubular worldsheets between two arbitrary boundary curves to that of this triangle.

We develop a Caffarelli-Silvestre extension theory for the fully fractional Schrödinger operator $\mathcal{L}^s=(\partial_t-i\Delta_x)^s$, $0<s<1$, nonlocal in space and time, whose Cauchy problem prescribes a past history rather than data at a single time. Every extension theory so far rests on positivity or sectoriality of the generator; here the semigroup is unitary and the symbol changes sign across the characteristic paraboloid, so neither is at hand. We construct the extension nonetheless. Its Poisson kernel, computed explicitly, is oscillatory rather than positive; its Dirichlet-to-Neumann map is $\mathcal{L}^s$; its normalising constant is the Caffarelli-Silvestre constant times the phase $e^{i\frac{\pi s}{2}}$; and the theory undergoes a transition at $s=\frac12$. We then identify the intrinsic Hilbert space of histories, on which the Poisson lifting is an isometry. This closes a circle. A lifted history is precisely an initial datum for the singular Schrödinger equation with nonlinear Neumann interaction of our companion paper, whose well-posedness theory therefore transfers to $\mathcal{L}^su=\mu|u|^{p-1}u$ with prescribed past. With our earlier work on the Bessel operator on a half-line, the three papers form a single programme, of which the present one is the closing step.

We consider beta-Dyson Brownian motion, with $\beta >= 1$, started from a deterministic configuration with uniformly bounded support. Let $\mu_t$ be the semicircular free-convolution flow issued from the initial empirical measure, and set $S_t = supp(\mu_t)$. For every fixed $T, \epsilon > 0$, with probability at least $1 - C \exp(-(log n)^2)$, every particle remains within an epsilon-neighborhood of $S_t$ for all $0 <= t <= T$. The result holds from time zero, requires no regularity assumption at the initial spectral edges, and applies to multi-cut supports with macroscopic interior gaps. A key ingredient is a deterministic local resolvent exclusion principle: an $o((n \eta)^(-1))$ comparison of Stieltjes transforms on a complex disc above a real point separated from the reference support excludes eigenvalues from the corresponding real interval. This gives a model-independent mechanism for converting local resolvent estimates into spectral confinement.

In 1966, McKean asked whether the entropy production of the Boltzmann equation must be monotone decreasing in time. We show that this is not the case even in the space-homogeneous setting for the Boltzmann collision operator with a constant angular cross section and the kinetic parameter $\gamma \in [0,1]$. This recovers the classical case of hard spheres and the simplest case of Maxwell molecules. Our examples are radially symmetric mixtures of Maxwellians.

Relative periodic orbits (RPOs) are solutions of the three-body problem that are periodic in a uniformly rotating reference frame and, in general, quasi-periodic in inertial coordinates. We present a numerical procedure for computing and continuing one-parameter families of RPOs of the planar Newtonian three-body problem. The method exploits consecutive syzygies, understood here as configurations in which the three bodies are aligned and their velocities satisfy the corresponding symmetry conditions. Matching the positions and momenta at two consecutive alignments reduces the computation of RPOs to a low-dimensional nonlinear problem. Its solutions are then numerically continued, and linear stability is determined from the nontrivial eigenvalues of the rotated monodromy matrix after removing the neutral directions associated with conserved quantities and continuous symmetries. The procedure is applied to several mass distributions and initial configurations, producing families of Poincaré, Hill, and binary-type solutions. These families exhibit transitions from nearly circular to highly eccentric motion, changes of stability near resonances and turning points, and absolute periodic solutions when the rotation angle is a rational multiple of 2{\pi}. In the Hill families, the continuation connects satellite configurations with circumstellar motion as the smallest body loses its gravitational binding to the intermediate body. Circumbinary and circumstellar configurations are also obtained in the binary regime. The results illustrate the dynamical diversity of RPOs and provide coherent three-body motions that can be used as prescribed trajectories in restricted four-body models.

Additive functionals $\overline{J}_t=\frac{1}{t}\int_0^tU(X_s)\circ dX_s$ of Stratonovich-type recently attracted much attention in the context of inference of thermodynamic properties of complex systems from observations $U$ of individual fluctuating paths $(X_s)_{0\le s\le t}$, whereby $X_0$ is initiated from some general measure. Concentration results on $\overline{J}_t$, albeit desirable, are virtually nonexistent. They turn out to be significantly more challenging to prove than for classical Lebesgue-type functionals $\overline{\rho}_t=\frac{1}{t}\int_0^t V(X_s)ds$ because the tilt deforms the full second-order structure of the Feynman-Kac generator instead of contributing an additive potential. This renders the generator generally non-self-adjoint even under detailed balance. We overcome this by working with a symmetrized Dirichlet form with a new effective potential that now couples the observable to the non-equilibrium character of the dynamics. We prove concentration inequalities for $\overline{J}_t$ for any bounded, sufficiently smooth vector-valued function $U$ of a general geometrically ergodic diffusion process $X_s$, including explicit sub-gamma and Bernstein-type inequalities, and we obtain explicit upper bounds on ${\rm Var}(\overline{J}_t)$. Strikingly, under detailed balance the concentration of $\overline{J}_t$ is distinctively sub-Gaussian at all times and all deviations, with a variance proxy fixed by the noise alone and independent of the spectral gap, which has no analog for $\overline{\rho}_t$.

We study regularisation and irreducibility properties of Gaussian quantum Markov semigroups (GQMSs) acting on continuous-variable quantum systems. We first identify a natural notion of regularity for operators in this setting, which allows us to formulate and characterise the smoothing effects of GQMSs in terms of algebraic conditions involving the drift and quantum diffusion matrices. These conditions establish a connection with the controllability theory of quantum linear systems and with the structure of decoherence-free subsystems. We then characterise irreducibility through several equivalent algebraic criteria: one formulated in terms of the drift and quantum diffusion matrices, one in terms of the operators appearing in the generalised GKLS representation of the generator, and a third given by a quantum analogue of Hörmander's condition. A central and somewhat surprising consequence is that, in contrast with the classical case, irreducibility is strictly stronger than conditions ensuring regularisation. Our results provide an algebraic framework for analysing these properties and lay the groundwork for a broader study of reducible Gaussian quantum Markov semigroups and, more generally, more general relevant quantum Markov semigroups on continuous-variable systems.