Papers

A compact formulation of electrostatics is presented. Without using constitutive relations, including the aether relations, an expression for the potential energy of a charged region under a potential function is proposed, where the charge distribution is specified by the charge potential field. Assuming that during a virtual motion of the charge, the virtual work expended by the field is equal to minus the time derivative of the potential energy, an expression for the mechanical force functional is derived. The force functional contains the action of an asymmetric active stress field D_{i}E_{j}, the skew-symmetric part of which is the mechanical couple density.

We construct a finite discrete realization of the Smorodinsky–Winternitz I superintegrable system on a triangular region of the two-dimensional square lattice. The construction is based on a pair of commuting number operators with finite spectrum together with an associated ladder-operator structure. We show that the resulting model is maximally superintegrable and that its symmetry algebra admits a Hahn-algebra presentation. Its spectral problem is solved exactly in terms of the bivariate dual Hahn polynomials of Tratnik type. Finally, we show that the continuum limit recovers the continuous Smorodinsky–Winternitz I system together with its ladder operators, eigenfunctions and symmetry algebra, thereby establishing the present construction as a genuine finite discrete realization of the continuous model.

For every Lie–Rinehart algebra, we construct an $F$-manifold algebra on the direct sum of its base algebra and module. Contrary to the assertion in \cite[Proposition 13.3.26]{LodayVallette}, the resulting structure is generally not Poisson. We determine when powers of the positive-degree ideal in the associated symmetric Poisson algebra are Poisson ideals, and relate the Leibnizator to the Lie–Rinehart differential. For a finite projective module of constant rank, the trace of the Leibnizator recovers the anchor and yields a rigidity result for injective anchors. We conclude with algebraic and geometric examples.

We analyze the dynamics of two harmonic oscillators with nonlinear, non-Hermitian coupling between them. Specifically, we consider a synthetic antiferromagnet composed of two ferromagnetic layers with antiferromagnetic interactions between the two layers, leading to the emergence of two oscillation modes of the combined system. We analyze the response of this system to external driving fields. We determine the allowed steady states under various driving conditions. We then analyze the dynamics of small deviations away from the different steady states, which reveals the normal modes that can be probed using spectroscopic techniques in experiment. As expected, the nonlinear system exhibits multistability, in which multiple steady states can exist for the same driving conditions. We find that the dynamical response can be drastically different depending on which one of the two modes is driven and on the strength of the driving. Specifically, we can obtain level-repulsion or level-attraction patterns in the spectrum. In addition to the rich variety of spectra, we find that some steady states exhibit normal modes that contain multiple frequencies. Our results explain recent experiments on synthetic antiferromagnets and provide guidance for designing experiments that can explore previously unseen phenomena in these systems.

In this paper, we develop the thermodynamic formalism for intermittent full-branch interval maps with finitely many neutral fixed points and for a broad class of possibly discontinuous potentials. Under a hyperbolicity assumption on the potential, we prove the existence of a unique equilibrium state with positive entropy and exponential decay of correlations. Furthermore, we characterize the phase transitions for the one-parameter family $\beta\varphi$, with $\beta>0$, and prove that, at the critical parameter, multiple ergodic equilibrium states may coexist, at most one of which has positive entropy. Our approach relies on a detailed analysis of conformal measures.

We develop the first rigorous topological framework for an auxiliary boundary-integral model motivated by physical scaling identities for wave scattering from self-similar or fractal media. The model is periodic in logarithmic scale, and a multiplicative Bloch-Floquet-Zak transform fiberizes its interscale coupling. For sufficiently small, well-separated components, equilibrium densities define a computable finite-dimensional projected matrix, while Riesz projections select the corresponding invariant spectral subspaces of the full boundary symbol. Under suitable spectral-isolation and point-gap conditions, and after recentering the two families at their respective same-scale reference values, we prove that the determinant winding of the exact Riesz-reduced family agrees with that of the projected matrix family. For regular-simplex configurations, we derive explicit nonzero winding formulas and track the resulting local winding data across finite prefractal levels and along a geometric sequence of wavenumbers. We also formulate conditional winding data for multiple dilation centers and show that the Zak phase of the chiral Hermitianization equals $\pi$ times the point-gap winding modulo $2\pi$. Thus, the scale-periodic boundary model admits a rigorous topological reduction.

We construct rank-zero irregular vertex operators for the Neveu–Schwarz algebra as linear maps between irregular Verma modules of the same rank satisfying the usual superconformal commutation relations and an irregular asymptotic condition. Under the nondegeneracy assumption $\Lambda_{2p}\neq0$, we prove their existence and uniqueness. We then extend the decomposition theorem for the Neveu–Schwarz algebra to the irregular setting: the tensor product of the free-fermion Fock module with a rank-$p$ Neveu–Schwarz irregular Verma module decomposes into an infinite direct sum of tensor products of two Virasoro irregular Verma modules. We also prove a compatible decomposition of the irregular vertex operators into tensor products of Virasoro irregular vertex operators. Using suitable pairings and mode insertions, we derive bilinear differential equations for weighted sums of products of Virasoro irregular conformal blocks of types $(0,0,1)$ and $(0,2)$. After explicit parameter identifications and a gauge transformation in the Painlevé V case, the resulting bilinear differential operators agree with those appearing in the quantum Painlevé V and IV tau-function equations.

We compare several "pathological" notions of continuity that make every additive function continuous. Our goal is to determine which of these notions is best behaved and deserves to be called Cauchy-Hamel continuity. The conclusion is that a bilateral version of a previously introduced notion of Q-continuity seems the most promising candidate. Our work is relevant to the axiomatic foundations of mechanics, specifically to the problem of axiomatically characterizing the parallelogram rule for the composition of forces. It was noted by Darboux that a continuity assumption is needed for such a characterization of the parallelogram rule and that, absent such a continuity axiom, "exotic" physical models based on alternative composition rules may exist. An intriguing question is whether such exotic physical models still possess some weak, residual notion of continuity. The concepts studied in this paper offer a first hint of a response to this question.

In our previous paper, we described the proposal and validation of a new physical implementation-type simulator (spontaneous symmetry breaking machine (SSBM)) for a combinatorial optimization problem that utilizes a phenomenon dual to the model in which a pseudo-spin many-body system is created by spontaneous symmetry breaking. Furthermore, we performed numerical simulations with different initial fluctuations on a large-scale benchmark problem (K_{2000}) and reported that we identified a condition under which pseudo-spin patterns gradually converge during the solution search process to ultimately form a single pattern and that the cut value of this converged solution reached 99.7% of the known best. In this paper, we discuss the conditions under which the SSBM can exhibit behavior analogous to that of a continuous-variable-type Ising machine, and report that by adding an effect similar to physical annealing to these conditions, the SSBM can be made to search for the best known cuts in K_{2000}.

Uncertainty in nonlinear dynamical systems is often organized by transport structures that are not apparent from posterior geometry alone. We introduce Dynamical Uncertainty Geometry (DUG), a framework that separates three ingredients that are frequently conflated: posterior uncertainty, future observations, and model-induced transport. DUG indexes posterior credible sets and transport-conditioned subsets by a common enclosed-probability coordinate, allowing physically meaningful classes to be tracked across credibility levels rather than examined at a single threshold. The framework combines a mass-ranked credible filtration, a description of future experiments through their induced probability laws and Fisher geometry, and a labelled transport skeleton representing dynamically distinct outcomes. We establish stability results for the resulting mass-indexed persistence modules under simultaneous perturbations of posterior density and probability measure, including continuum, refinement, and finite atomic formulations. These results provide quantitative control of persistence computed from numerical approximations and weighted grids. Two benchmark problems illustrate the framework. In a short-arc orbit-determination setting, DUG identifies dynamically distinct return classes within a connected uncertainty region and highlights a difference between information-based and local Fisher-based observation-design criteria. In the Earth-Moon planar circular restricted three-body problem, DUG reveals multiple first-hit transport outcomes coexisting within a connected credible region and provides diagnostics for assessing their numerical resolution. Together, these examples show how topological summaries of posterior geometry, when coupled to transport labels and future experiments, yield a richer description of uncertainty than either posterior probabilities or dynamical classifications alone.

The number of non-Clifford gates needed to perform a task, or simply magic, is a resource for fault-tolerant quantum computation. Von Neumann algebras provide a formal mathematical structure to describe infinite-dimensional quantum systems, such as those in quantum field theory or quantum statistical mechanics. A particularly important class of von Neumann algebras is called Type III algebras, for which the standard notions of finite-dimensional systems such as density matrices and traces break down. In this work, we argue that Type III von Neumann algebras fundamentally require an infinite amount of magic. Specifically, we consider the thermodynamic limit of finite-dimensional quantum systems, such as lattice systems, and show that if the states in the thermodynamic limit possess only a bounded amount of magic, the resulting local von Neumann algebra cannot be of Type III. Our result has direct implications for the quantum simulations of quantum field theories, for which the algebra of a local subregion is known to be of Type III.

Real-to-complex spectral transitions mark the onset of amplification in non-Hermitian systems, but their thresholds are often treated as model-specific quantities. Here we develop a general, non-perturbative activated-mode principle that governs the real-to-complex threshold across broad classes of non-Hermitian systems. A central insight is that only a small Hilbert subspace is “activated" at the transition onset, which can be variationally determined through the competition between spectral detuning and mode-level projected non-Hermitian couplings. The result is a closed-form exceptional-activation condition for arbitrarily large “disturbances", rather than a perturbative estimate. We apply our framework to three contrasting illustrative problems, establishing (i) a closed-form threshold for critical non-Hermitian skin amplification at all system sizes; (ii) a new link between impurity tunneling threshold and exceptional point switching; and (iii) activation channel switching without underlying topological phase transition. Overall, our findings recast real-to-complex transitions as generic mode-selection problems independent of any specific symmetry.

We study a spin-$\frac12$ impurity coupled to the boundary of a strongly correlated spin-$s$ Takhtajan–Babujian chain, an integrable model whose low-energy physics is described by a perturbed $SU(2)_{2s}$ Wess–Zumino–Witten conformal field theory. While boundary conformal field theory determines the low-energy universality class of the weak-coupling regime, exact Bethe Ansatz methods reveal a sequence of boundary quantum phase transitions in which impurity-bound states emerge and reorganize the Hilbert space into multiple excitation towers built on distinct boundary configurations. This tower restructuring provides the organizing principle for a rich boundary phase diagram extending beyond the conventional Kondo regime. Weak antiferromagnetic coupling realizes the overscreened $2s$-channel Kondo universality class, whereas stronger couplings generate localized boundary modes and qualitatively new screening mechanisms. To describe the resulting thermodynamics, we develop a generalized thermodynamic Bethe Ansatz framework that captures the multi-tower structure across all regimes. The impurity entropy reproduces the boundary conformal field theory prediction in the overscreened Kondo regime but develops pronounced nonmonotonic temperature dependence once boundary-bound states appear, in quantitative agreement with large-scale finite-temperature matrix-product-operator simulations. Complementary dynamical calculations reveal sharp threshold features in the impurity spectral function that directly track the underlying tower structure. Together, boundary conformal field theory, exact Bethe Ansatz, generalized thermodynamic Bethe Ansatz, and tensor-network simulations provide a unified description of impurity screening, boundary-bound-state formation, and excitation-tower reconstruction in a correlated spin-$s$ chain.

We develop a probabilistic cavity-contiguity framework for proving replica-symmetric convergence of local marginals in mean-field Gibbs systems. The approach is based on cavity decompositions of the Hamiltonian together with direct comparison of probability measures through Radon-Nikodym derivatives and Hellinger-type estimates. At a conceptual level, the method separates the concentration of the relevant order parameters from the identification of the asymptotic cavity model and the comparison of the associated Gibbs measures. In contrast with interpolation-based approaches, the argument relies only weakly on the Gaussianity of the disorder and naturally accommodates concentration tools such as Poincaré and log-Sobolev inequalities. Rather than pursuing maximal generality, with the aim of making the method transparent, we implement the framework in a canonical example: the high-temperature Sherrington-Kirkpatrick model. In this setting, we prove that the marginal law of a fixed spin converges in total variation toward the effective one-dimensional cavity measure predicted by the replica method. Beyond the specific result for the Sherrington-Kirkpatrick model, the paper illustrates a broader cavity-contiguity methodology which is expected to extend naturally to other mean-field Gibbs systems, particularly Bayesian inference models with or without mismatch.

The min-plus process is a stochastic coagulation-annihilation-type process on the binary tree, of interest in mathematics, physics, and computer science as a tractable instance of max-type recursive distributional equations. We carry out large Monte Carlo simulations at effective tree depths up to $N=60$ that provide finite-depth corroboration of the Beta(2,1) stretched-exponential limit for its root value $X_{N}$ at $p=1/2$, on the asymmetric $\sqrt{N}$ side of the random-homogeneous-systems classification recently introduced by Chen, Duquesne, and Shi and by Morfe. Off criticality, our simulations confirm the sub-critical closed form $\mathbb{P}(X_{\infty}=1)=(1-2p)/(1-p)$ within Monte Carlo error and document a super-critical mean growth exceeding the elementary $(2p)^{N}$ lower bound at the depths we reach. For a Bernoulli($q$)-initial-condition variant, we identify an elementary closed-form identity at $p=1/2$ that pins down the order parameter $\mathbb{P}(X_{N}=0)=q$ exactly, locates the absorbing-state phase transition at $p_{c}=1/2$ in the operator-mixing probability rather than in the initial-zero density, and shows that the conditional law on positives deforms substantially with $q$. Our simulations use a level-wise recursion and an FFT-based precomputed leaf table which reduce the effective simulation depth while preserving the recursive tree law and may be useful for the simulation of related recursive equations on large trees.

Many nonlinear models across physics, chemistry, and biology exhibit multiple solutions for the same parameters, and capturing this entire solution set is essential for understanding pattern-forming systems. Yet existing learned surrogates are fundamentally single-valued: neural operators map each parameter to a single output, and physics-informed neural networks converge to one branch. We develop PatternFormer (PF), a large language model-based framework for learning the multiple solutions of nonlinear partial differential equations. By transforming unordered coexisting solutions into canonical sequences, PF produces structured solution sets in a single autoregressive pass, terminating automatically for finite families and enforcing physical residual constraints for unbounded ones. On nonlinear elliptic problems it recovers all solution branches in one inference step; on Gray–Scott it generates coexisting Turing patterns, including physically valid states absent from the reference data and beyond training. PF can also be sequentially fine-tuned across multistable systems, toward general foundation models for solution landscapes.

While the solution to the $1+1$ dimensional stochastic heat equation with multiplicative noise is closely related to the exponential of a Brownian motion, the two-dimensional picture exhibits an additional weak-to-strong disorder transition. In [CSZ '23], the critical two-dimensional stochastic heat flow (SHF) was constructed as the scaling limit of the partition function of $2+1$ dimensional directed polymers under the logarithmic intermediate-disorder scaling at criticality. The SHF is a random measure and, like many naturally occurring random measures, it is expected to exhibit rich intermittency. [CSZ '25] established that it is almost surely singular with respect to the Lebesgue measure. More recently, [GT '26] showed that the logarithm of the SHF averaged over small balls is asymptotically Gaussian, with both its mean and variance diverging as the ball radius tends to zero. In this paper we prove a sharp result quantifying the singularity of the support of the SHF as well as its intermittency. In particular, we show that, almost surely, for all small $\varepsilon>0$, up to a vanishing error, all the mass of the point-to-plane SHF in any domain is concentrated on ${1}/{\big(\varepsilon^2\log^{1/2+o(1)}(1/\varepsilon)\big)}$ balls of radius $\varepsilon$, each containing $\varepsilon^2{\log^{1/2+o(1)}(1/\varepsilon)}$ mass, thus precisely establishing its logarithmic fractal behavior. A key ingredient in the proof is a refined large-deviations theory, which allows access to conditional distributions, by taking advantage of the Gaussian-like behavior of the SHF at quasi-critical scales. A further useful observation that features prominently is that conditioning a Brownian motion on its endpoint being unusually large essentially induces a shift in the mean of its increments, and consequently, at small enough scales, their distributions do not alter significantly.

A nonlinear model of a memristor based on charged mobile vacancies is considered taking into account the electrostatic interaction between them. This interaction significantly affects the stationary (limiting) vacancy distributions formed under the action of the electric current flowing through the memristor, for which analytical expressions are obtained in this work. Between the regions with reduced and increased vacancy concentrations, an intermediate electrically neutral region is formed due to the electrostatic interaction. Interestingly, the limiting resistances in the “on” and “off” states of such a memristor do not depend on the strength of the electrostatic interaction, at least in the leading first order in the vacancy concentration.

Let $\mathcal{P}_{N}$ be the $(N+1)$-dimensional Hilbert space of analytic polynomials of degree at most $N$. This is the natural environment to define $SU(2)$ (Bloch) coherent states. Let $Q_{\rho }$ be the Husimi function of a density operator $\rho $ on $% \mathcal{P}_{N}$. We prove an isospectral version of Lieb-Solovej inequality: if $\rho ^{\downarrow }$ is obtained by placing the eigenvalues of $\rho $ in decreasing order along the monomial basis, then \begin{equation*} \int_{\mathbb{C}}\Phi (Q_{\rho }(z))\,dm(z)\leq \int_{\mathbb{C}}\Phi (Q_{\rho ^{\downarrow }}(z))\,dm(z) \end{equation*}% for every convex function $\Phi $ on $[0,1]$. Applying the corresponding reversed inequality to the concave function $\Phi (t)=-t\log t$ gives the Wehrl entropy. In the process it is show that the output state of $\rho $ under Lieb-Solovej's channel is majorized by the output of the state $\rho ^{\downarrow }$. As an application, among all measurable subsets of the sphere having a fixed area, spherical caps maximize every partial sum of the eigenvalues of Toeplitz operators with symbol $1_{\Omega}$. Equivalently, caps maximize all Ky Fan norms, leading to the sharp inequalities, previously known for $r=1$: \begin{equation*} \sum_{j=1}^{r}\lambda_{j}(\Omega) \leq r-\sum_{k=0}^{r-1}(r-k)\binom{N+1}{k} m(\Omega)^{k}\bigl(1-m(\Omega)\bigr)^{N+1-k}. \end{equation*} This implies isoperimetric inequalities for all Schatten sums of $T_{\Omega}$, obtained without using the spherical isoperimetric inequality.

We study the quantum geometry of doubly degenerate energy levels in a four-dimensional parameter space. We find that the scalar quantum metric $g$ and the Berry curvature $F$ obey $(\operatorname{tr} g)^2/16\geq\sqrt{\det g}\geq |2\Tr(F\wedge F)-\Tr F\wedge \Tr F|/24$. The first inequality characterizes the anisotropy in the metric. The second determinant inequality measures the self-duality of the traceless $SU(2)$ part of the curvature under Hodge star operation and the algebraic closedness of the inter-level polarization amplitudes under $SU(2)$ rotations in the doubly degenerate levels. The saturation of the determinant bound imposes a quaternionic Cauchy-Riemann equation, analogous to the complex analyticity imposed by ideal-band conditions in two-dimensional Chern insulators. As examples, four-band Dirac Hamiltonians automatically saturate the determinant bound and possess a topological zero in $\operatorname{tr}(F\wedge F)$. We compare the differences between Kramers degeneracy and ordinary $U(2)$ degeneracy. In addition to the non-Abelian geometric bound, the latter also obeys an independent first-Chern bound.

A rigorous framework for consistent study of junctions in a coordinate independent manner is presented. This is achieved by extending a theory of distributions in curved spacetimes with torsion and combining it with covariant formalism. In this way, one can derive general conditions on the differentiability of the joined manifold. Given a specific theory of gravity, in particular the Einstein-Cartan-Sciama-Kibble one, we evaluate the conditions for obtaining a smooth junction. These conditions can be successfully applied independently of the choice of coordinates used to describe the metric, thereby, evading the drawbacks of coordinate dependent Israel-Darmois framework.

Boundary conditions for a drift-kinetic model at the collisional presheath entrance with perpendicular incidence of the magnetic field to the wall are derived and numerically implemented. The drift-kinetic model for the plasma is based on the expansion of the ion distribution function on a Hermite-Laguerre basis, and the evolution of the resulting gyromoments. A linear-plasma-device geometry is considered. Comparison with simpler simulations with previously used ad hoc boundary conditions is presented. For the new set of boundary conditions, a significant increase of the plasma outflow to the wall is observed, leading to a significantly smaller plasma density in the whole volume of the device.

We study the stability of one-dimensional planar isotropic–nematic interfaces in the Landau–de Gennes model with anisotropic elastic constant $L$. Earlier work proved the instability for $L<0$ only under an extra condition. We remove this condition and prove that any non-negative minimizer of the reduced energy within the diagonal planar class is unstable under general one-dimensional perturbations throughout $-\frac{3}{2}<L<0$. This result shows that $L=0$ is the sharp endpoint of the instability range from the negative-$L$ side: the critical operator is non-negative at $L=0$, while instability holds over the full negative-$L$ range of the reduced problem. We also define the optimal stability index and correct a factor-of-two normalization inconsistency in a previously stated interface profile.

The statistics of current fluctuations has long been a central object of study in non-equilibrium physics. Most existing work has focused on one-time statistics, while their multi-time generalisation remains comparatively less explored. We address this gap by extending the fluctuating hydrodynamics framework of Macroscopic Fluctuation Theory (MFT) to study multi-time statistics in the non-stationary state of a diffusive system on an infinite line. For the simplest cases of a non-interacting lattice gas and hard-core Brownian point particles, we present explicit solution of the MFT leading to multi-time large-deviation statistics. For generic systems, the MFT is solved perturbatively, yielding explicit results for two-time correlations. These reveal that the connection between current fluctuations and fractional Brownian motion, previously observed for flat initial conditions, does not persist for step initial conditions. We independently verify these hydrodynamic results by solving the corresponding microscopic dynamics for the non-interacting gas and for the symmetric simple exclusion process. Additional confirmation comes from numerical simulations.

Reconstruction of 3D electromagnetic currents from sparse, multi-frequency, far-field radiation data is severely ill-posed because sparse sampling masks specific current Fourier modes and the missing spectrum enlarges the null space of the forward operator. Although structured-Hankel completion can recover missing spectrum by exploiting the finite rate of innovations structure of a compactly supported geometrically sparse source, it treats current components independently without guaranteeing compatibility with the model. This often renders non-physical solutions that violate model consistency. To resolve this, we propose a training-free Maxwell-model-consistent 3D structured-Hankel framework (MC-Hankel) that alternates joint low-rank Hankel completion with an exact, closed-form projection onto the true Maxwell synthesis subspace while enforcing conjugate symmetry and noise-aware data fidelity. Across tests on two source models, three sparse sampling rates (30%, 40%, and 50% of a Nyquist-sampled grid), and clean versus noisy (10 dB SNR additive white Gaussian noise) conditions, MC-Hankel consistently outperforms sub-sampled Fourier inversion, scale-normalized $\ell_1$-compressed sensing baseline, and the standard joint 3D annihilating filter-based low-rank Hankel matrix completion approach (ALOHA). Across all considered source-noise-sampling configurations, MC-Hankel improves mean peak SNR (PSNR) by 0.65-4.22 dB and tri-planer structural similarity (3D SSIM) by $0.0266$-$0.0703$, reduces relative full-volume $\ell_2$ reconstruction error by $7.4\%$-$38.9\%$, and brings non-physical model residuals down to machine precision over standard joint 3D ALOHA with merely 17% computational overhead. Exact signed-rank tests are underpowered for five paired trials, so statistical conclusions are reported together with paired bootstrap intervals, effect sizes, and small-sample limitations.

In this paper, we investigate the integrability of Lotka-Volterra (replicator) systems arising from interaction matrices generated from corresponding graph structures, continuing work started by Visomirski and Griffin [J. Phys. A., 58:015701, 2025] and Evripidou et al. [J. Phys. A., 55:325201, 2022] (among others). In particular, we define a new family of graphs, the spoked graphs, and show that all dynamics generated from this family are integrable. In reference to Evripidou et al. (2022), we define a new anti-cloning operator and show that its action on balanced tournament graphs (with odd vertex count) generates new graphs whose dynamics are integrable. Interestingly, we provide numerical evidence that this anti-cloning operation leads to chaotic behaviour when applied to other graph families (e.g., the directed cycles that generate the classically integrable Volterra lattice). This work completes a taxonomy of all integrable dynamics generated by directed graphs with up to six vertices started by Visomirski and Griffin (2025), and suggests several future directions of study on this topic.

The constrained equilibrium problem for the logarithmic potential was introduced by Rakhmanov (1996) as he realized that the asymptotic distribution of the zeros of discrete orthogonal polynomials in a compact interval could be described in terms of the equilibrium measure of this interval in a class of measures subject to a certain constraint. Other authors such as Dragnev and Saff (1997), and Kuijlaars and Van Assche (1999), extended this approach to more general settings. These problems have proven to be useful for describing asymptotic distributions in different settings. In the current paper, we consider constrained equilibrium problems in the Riesz setting, that is, for $s$-Riesz potentials in the hyperplane $\mathbb{R}^d,$ with $d\geq 1$ and $\max (0,d-2) < s < d$. Along with some general results, an illustrative example consisting of the solution of a constrained equilibrium problem in the unit ball is studied in detail. This is the main part of the paper. Finally, a number of numerical experiments show how the constrained equilibrium measure, the solution of the problem, may be discretized using the so-called 'constrained Leja points' introduced by Coroian and Dragnev (2001).

We use a variety of analytic and numerical techniques to approximate the geometry of individual dark matter halos when modeled as collections of ultralight scalar bosons (ULDM), followed by an analysis on the future directions and implications of this work. Ultimately, we aim to understand the morphology and creation of complex structures (like spirals) in galaxies from a geometric standpoint and to investigate the connection between dark matter and baryonic matter in galaxies. We find that superpositions of gravitational eigenmodes for ULDM in a fixed background potential can yield spiral density perturbations which rotate on the order of one hundred million years and appear within ten to twenty kpc of the halo center, both features comparable to realistic galaxies, and that this model provides a natural explanation for the density waves necessary to seed baryonic spirals.

We prove a conifold gap theorem for the topological recursion of toric mirror curves: for every local analytic family acquiring a generic one-node degeneration and every fixed genus at least two, the conifold-polarized free energy has one universal polar term, and the remainder is jointly holomorphic in the transverse and spectator parameters. The result covers separating and nonseparating nodes and allows general family deformations beyond the pure filling-fraction case.

Coherences between different energy levels are strongly constrained by thermodynamics. Here we ask a related question: if several coherence transfers can be optimized separately, can they also be optimized simultaneously by the same thermodynamic process? We show that this is in general a compatibility problem. Using an energy-resolved Stinespring representation of thermal operations, we express elementary coherence-transfer amplitudes as inner products of bath-weighted vectors associated with total-energy shells. Saturation of the corresponding Cauchy–Schwarz bounds requires these vectors to be collinear. Since all transfers have to arise from the same energy-preserving system–bath unitary, such saturation conditions cannot in general be imposed independently. Trace and Gibbs preservation lead to additional phase-closure conditions and to polygon inequalities that can rule out simultaneous saturation for fixed population-transfer data. We complement these analytic results with a rigorous computer-assisted semidefinite-program analysis. For the mixed-degenerate Hamiltonian \(H=\operatorname{diag}(0,1,1,3)\), machine-certified primal and dual bounds show that no covariant Gibbs-preserving channel can simultaneously attain three individually optimal directional coherence transfers. Since every thermal operation is covariant and Gibbs preserving, the same triple of \(\CGP\)-optimal values cannot be attained by a thermal operation. We also analyze auxiliary Choi-sector Gram matrices for mixed-degenerate ladders. Under specified directional rank assumptions, positivity imposes range-inclusion conditions, from which we obtain conditional non-genericity statements for physical joint saturation.