Daily Papers

We use a simple method that leads to the integrals involved in obtaining the traveling wave solutions of wave equations with one and two exponential nonlinearities. When the constant term in the integrand is zero, implicit solutions in terms of hypergeometric functions are obtained while when that term is nonzero all the basic traveling wave solutions of Liouville, Tzitzeica and their variants, as well as sine/sinh-Gordon equations with important applications in the phenomenology of nonlinear physics and dynamical systems are found through a detailed study of the corresponding elliptic equations

We study a one dimensional nonlinear Klein-Gordon equation with a potential, a localized quadratic nonlinearity, and a non-localized quartic nonlinearity. We assume that the linearized operator has a single discrete eigenvalue below the continuous spectrum, corresponding to an ''internal mode''. This eigenvalue generates localized, time-periodic solutions for the linear equation. Assuming the natural Fermi Golden Rule, we give a global description of the dynamics of the amplitude of the internal mode through the radiation damping mechanism, and prove scattering for the radiation. Our approach is based on the distorted Fourier transform and a collection of refined dispersive decay and smoothing estimates. The corresponding cubic problem is closely related to the question of global-in-space dispersive asymptotics for perturbations of kink solutions in the classical ϕ^4 model, which remains a challenging open question. From the perspective of nonlinear estimates, a quartic (non-localized) interaction is essentially sharp for our analysis, and closing the estimates with such a low degree nonlinear term requires a careful and delicate treatment. As a direct application of our analysis, we prove a result on the asymptotic stability of kinks for classes of scalar field models that can be seen as perturbations of the ϕ^4 model.

In the quantum inverse scattering method, the sine-Gordon model is solved on a lattice, and its continuum and infinite-volume limits become uniform only after the length of the box is rescaled. We show that this rescaling, supplemented by an on-shell condition, defines a renormalization scheme. The whole dependence on the cut-off is absorbed into a factor Z_L multiplying the length of the box, with an exponent fixed by the scaling dimension of the vertex operator; neither the field nor the coupling β is renormalized, no subtraction scale is introduced, and the mass parameter receives a finite renormalization, fixed by requiring that the lightest breather has the mass of the boson in the linearized theory. The soliton and breather masses follow from the eigenvalues of the monodromy operator as functions of parameters that do not run, and reduce to the classical and semiclassical results as βrightarrow0. In particular, the mass ratios are produced directly from the regularized construction, while the relation between the lattice parameters and the physical scale is absorbed into Z_L and never needs to be computed. The scheme is a reparameterization of the theory which makes explicit the interplay between integrability and renormalization. We compare it with Coleman's normal ordering and with conformal perturbation theory, and we argue that it extends to ultralocal integrable theories with a massive spectrum, possibly including asymptotically free models that admit an ultralocal lattice regularization.

This paper derives the second-order positive Burgers' equation from the standard Burgers' hierarchy to explore its complete integrability and exact analytical solutions. We construct this higher-order nonlinear evolution equation by systematically applying the recursion operator to the classical Burgers' equation. Expanding on this structural framework, we then derive the explicit third-order equation and use Complete Bell Polynomials to generalize the n-th order hierarchy. Through the Cole-Hopf transformation, the second-order nonlinear equation rigorously maps to the linear third-order dispersion equation. We establish complete integrability by explicitly formulating the scalar Lax pair (zero-curvature representation), which allows us to directly derive the associated differential and algebraic Darboux transformations. To systematically classify explicit stationary, time-dependent traveling wave, and self-similar profiles, we apply several techniques to the linear domain: separation of variables, the Hirota perturbation method, traveling wave reduction via the Complete Discrimination System for Polynomial Method (CDSPM), and Lie similarity reduction. Crucially, our analysis demonstrates the exact finite truncation of the Hirota perturbation series. We conclude by outlining how this equation impacts the theoretical understanding of fluid dynamics, nonlinear transport phenomena, and higher-order wave propagation.

In this manuscript we set up the direct and inverse scattering problems for step-like traveling-wave solutions of the nonlinear Schrödinger equation. Specifically, we consider initial data u(x,0) satisfying u(x,0)to u_0^ell(x) as xto-infty and u(x,0)to u_0^r(x) as xto+infty, where u_0^ell(x) and u_0^r(x) are elliptic traveling waves. Under suitable assumptions on the initial data we formulate the direct scattering problem and establish analytic properties of the scattering data. We then formulate the inverse problem as a Riemann--Hilbert problem and prove its solvability. Finally, we observe that this Riemann--Hilbert formulation is a special case of the one arising for full soliton-gas initial data.

We emphasize two connections, one well known and another less known, between the dissipative nonlinear second order differential equations and the Abel equations which in its first kind form have only cubic and quadratic terms. Then, employing an old integrability criterion due to Chiellini, we introduce the corresponding integrable dissipative equations. For illustration, we present the cases of some integrable dissipative Fisher, nonlinear pendulum, and Burgers-Huxley type equations which are obtained in this way and can be of interest in applications. We also show how to obtain Abel solutions directly from the factorization of second-order nonlinear equations

A one-parameter family of Emden-Fowler equations defined by Lampariello's parameter p which, upon using Thomas-Fermi boundary conditions, turns into a set of generalized Thomas-Fermi equations comprising the standard Thomas-Fermi equation for p=1 is studied in this paper. The entire family is shown to be non integrable by reduction to the corresponding Abel equations whose invariants do not satisfy a known integrability condition. We also discuss the equivalent dynamical system of equations for the standard Thomas-Fermi equation and perform its phase-plane analysis. The results of the latter analysis are similar for the whole class.

We present an approach to the bright soliton solution of the NLS equation from the standpoint of introducing a constant potential term in the equation. We discuss a `nongauge' bright soliton for which both the envelope and the phase depend only on the traveling variable. We also construct a family of generalized NLS equations with solitonic sech^p solutions in the traveling variable and find an exact equivalence with other nonlinear equations, such as the Korteveg-de Vries and Benjamin-Bona-Mahony equations when p=2

We first reformulate and expand with several novel findings some of the basic results in the integrability of Abel equations. Next, these results are applied to Vein's Abel equation whose solutions are expressed in terms of the third order hyperbolic functions and a phase space analysis of the corresponding nonlinear oscillator is also provided

In this paper the effect of a small dissipation on waves is included to find exact solutions to the modified BBM equation. Using Lyapunov functions and dynamical systems theory, we prove that when viscosity is added to the BBM equation, in certain regions there still exist bounded traveling wave solutions in the form of solitary waves, periodic, and elliptic functions. By using the canonical form of Abel equation, the polynomial Appell invariant make the equation integrable in terms of Weierstrass wp functions. We will use a general formalism based on Ince's transformation to write the general solution of dissipative BBM in terms of wp functions, from which all the other known solutions can be obtained via simplifying assumptions. Using ODE analysis we show that the traveling wave speed is a bifurcation parameter that makes transition between different classes of waves.

We introduce a special type of dissipative Ermakov-Pinney equations of the form v_{ζζ}+g(v)v_ζ+h(v)=0, where h(v)=h_0(v)+cv^{-3} and the nonlinear dissipation g(v) is based on the corresponding Chiellini integrable Abel equation. When h_0(v) is a linear function, h_0(v)=λ^2v, general solutions are obtained following the Abel equation route. Based on particular solutions, we also provide general solutions containing a factor with the phase of the Milne type. In addition, the same kinds of general solutions are constructed for the cases of higher-order Reid nonlinearities. The Chiellini dissipative function is actually a dissipation-gain function because it can be negative on some intervals. We also examine the nonlinear case h_0(v)=Ω_0^2(v-v^2) and show that it leads to an integrable hyperelliptic case

In this work, we generalize the reaction-diffusion equation in statistical physics, Schr\"odinger equation in quantum mechanics, Helmholtz equation in paraxial optics into the neural partial differential equations (NPDE), which can be considered as the fundamental equations in the field of artificial intelligence research. We take finite difference method to discretize NPDE for finding numerical solution, and the basic building blocks of deep neural network architecture, including multi-layer perceptron, convolutional neural network and recurrent neural networks, are generated. The learning strategies, such as Adaptive moment estimation, L-BFGS, pseudoinverse learning algorithms and partial differential equation constrained optimization, are also presented. We believe it is of significance that presented clear physical image of interpretable deep neural networks, which makes it be possible for applying to analog computing device design, and pave the road to physical artificial intelligence.

Existing model reduction techniques for high-dimensional models of conservative partial differential equations (PDEs) encounter computational bottlenecks when dealing with systems featuring non-polynomial nonlinearities. This work presents a nonlinear model reduction method that employs lifting variable transformations to derive structure-preserving quadratic reduced-order models for conservative PDEs with general nonlinearities. We present an energy-quadratization strategy that defines the auxiliary variable in terms of the nonlinear term in the energy expression to derive an equivalent quadratic lifted system with quadratic system energy. The proposed strategy combined with proper orthogonal decomposition model reduction yields quadratic reduced-order models that conserve the quadratized lifted energy exactly in high dimensions. We demonstrate the proposed model reduction approach on four nonlinear conservative PDEs: the one-dimensional wave equation with exponential nonlinearity, the two-dimensional sine-Gordon equation, the two-dimensional Klein-Gordon equation with parametric dependence, and the two-dimensional Klein-Gordon-Zakharov equations. The numerical results show that the proposed lifting approach is competitive with the state-of-the-art structure-preserving hyper-reduction method in terms of both accuracy and computational efficiency in the online stage while providing significant computational gains in the offline stage.

In this paper, we use a traveling wave reduction or a so-called spatial approximation to comprehensively investigate periodic and solitary wave solutions of the modified Benjamin, Bona & Mahony equation (BBM) to include both dissipative and dispersive effects of viscous boundary layers. Under certain circumstances that depend on the traveling wave velocity, classes of periodic and solitary wave like solutions are obtained in terms of Jacobi elliptic functions. An ad-hoc theory based on the dissipative term is presented, in which we have found a set of solutions in terms of an implicit function. Using dynamical systems theory we prove that the solutions of BBMv experience a transcritical bifurcation for a certain velocity of the traveling wave. Finally, we present qualitative numerical results.

External potentials break the translational invariance of nonlinear Schrödinger (NLS) equations and generally convert a continuous family of standing wave solutions into isolated pinned states. We study this mechanism for the prototypical case of a kink in the NLS equation with competing focusing cubic and defocusing quintic nonlinearities. The Lyapunov--Schmidt reduction method is used to prove that the location s_0 of the pinned kinks is selected near roots of an effective potential V_{rm eff}(s). We then determine how the translational zero eigenvalue of the Jacobian operator splits under the perturbation and prove that the sign of V_{rm eff}'(s_0) distinguishes eigendirections of positive and negative energy. We further show that if the Jacobian operator admits a negative eigenvalue, then the stability problem contains a real unstable eigenvalue. The derivation of the asymptotic expansion for the unstable eigenvalue is nonstandard because the zero eigenvalue is embedded in the continuous spectrum of the linearized operator. Finally, we show that the persistence of kinks lacks parity symmetry even if the external potential is odd. Numerical continuation, spectral computations with domains adapted to the slowly decaying eigenfunctions, and direct dynamical simulations of the time evolution corroborate the analytical predictions for different representative examples.

In this paper we review the physical relevance of a Korteweg-de Vries (KdV) equation with higher-order dispersion terms which is used in the applied sciences and engineering. We also present exact traveling wave solutions to this generalized KdV equation using an elliptic function method which can be readily applied to any scalar evolution or wave equation with polynomial terms involving only odd derivatives. We show that the generalized KdV equation still supports hump-shaped solitary waves as well as cnoidal wave solutions provided that the coefficients satisfy specific algebraic constraints. Analytical solutions in closed form serve as benchmarks for numerical solvers or comparison with experimental data. They often correspond to homoclinic orbits in the phase space and serve as separatrices of stable and unstable regions. Some of the solutions presented in this paper correct, complement, and illustrate results previously reported in the literature, while others are novel.

A review of three-dimensional waves on deep-water is presented. Three forms of three dimensionality, namely oblique, forced and spontaneous type, are identified. An alternative formulation for these three-dimensional waves is given through cubic nonlinear Schrödinger equation. The periodic solutions of the cubic nonlinear Schrödinger equation are found using Weierstrass elliptic wp functions. It is shown that the classification of solutions depends on the boundary conditions, wavenumber and frequency. For certain parameters, Weierstrass wp functions are reduced to periodic, hyperbolic or Jacobi elliptic functions. It is demonstrated that some of these solutions do not have any physical significance. An analytical solution of cubic nonlinear Schrödinger equation with wind forcing is also obtained which results in how groups of waves are generated on the surface of deep water in the ocean. In this case the dependency on the energy-transfer parameter, from wind to waves, make either the groups of wave to grow initially and eventually dissipate, or simply decay or grow in time.

Extending our previous work (Rosu, Mancas, Chen, Ann.Phys. 343 (2014) 87-102), we define supersymmetric partner potentials through a particular Riccati solution of the form F(x)=(x-c)^2-1, where c is a real shift parameter, and work out the quartic double-well family of one-parameter isospectral potentials obtained by using the corresponding general Riccati solution. For these parametric double well potentials, we study how the localization properties of the two wells depend on the parameter of the potentials for various values of the shifting parameter. We also consider the supersymmetric parametric family of the first double-well potential in the Razavy chain of double well potentials corresponding to F(x)=(1/2)sinh 2x-2(1+sqrt 2)sinh 2x/[(1+sqrt 2) cosh 2x+1], both unshifted and shifted, to test and compare the localization properties

This study introduces amangkurat, an open-source Python library designed for the robust numerical simulation of relativistic scalar field dynamics governed by the nonlinear Klein-Gordon equation in (1+1)D spacetime. The software implements a hybrid computational strategy that couples Fourier pseudo-spectral spatial discretization with a symplectic Størmer-Verlet temporal integrator, ensuring both exponential spatial convergence for smooth solutions and long-term preservation of Hamiltonian structure. To optimize performance, the solver incorporates adaptive timestepping based on Courant-Friedrichs-Lewy (CFL) stability criteria and utilizes Just-In-Time (JIT) compilation for parallelized force computation. The library's capabilities are validated across four canonical physical regimes: dispersive linear wave propagation, static topological kink preservation in phi-fourth theory, integrable breather dynamics in the sine-Gordon model, and non-integrable kink-antikink collisions. Beyond standard numerical validation, this work establishes a multi-faceted analysis framework employing information-theoretic entropy metrics (Shannon, Rényi, and Tsallis), kernel density estimation, and phase space reconstruction to quantify the distinct phenomenological signatures of these regimes. Statistical hypothesis testing confirms that these scenarios represent statistically distinguishable dynamical populations. Benchmarks on standard workstation hardware demonstrate that the implementation achieves high computational efficiency, making it a viable platform for exploratory research and education in nonlinear field theory.

We develop a spacetime neural network method with second order optimization for solving quantum dynamics from the high dimensional Schr\"{o}dinger equation. In contrast to the standard iterative first order optimization and the time-dependent variational principle, our approach utilizes the implicit mid-point method and generates the solution for all spatial and temporal values simultaneously after optimization. We demonstrate the method in the Schr\"{o}dinger equation with a self-normalized autoregressive spacetime neural network construction. Future explorations for solving different high dimensional differential equations are discussed.

The coupled cubic nonlinear Schrödinger (CNLS) equations are used to study modulational instabilities of a pair of nonlinearly interacting two-dimensional waves in deep water. It has been shown that the full dynamics of these interacting waves gives rise to localized large-amplitude wavepackets (wave focusing). In this short letter we attempt to verify this result numerically using a Fourier spectral method for the CNLS equations.

We investigate whether the geometry of a curved conductor--dielectric interface can modify the effective mass and propagation of scalar surface plasmons (SSPs). In the hydrodynamic limit, we describe SSPs as density perturbations in a semiclassical jellium electron gas. Using a thin--layer dimensional reduction procedure, we derive a massive Klein--Gordon equation for SSPs on the curved interface with a geometric potential. We find that under specific conditions, the geometric potential allows SSP propagation with a vanishing effective mass. We then apply the formalism to a well known analogue black hole configuration, determining the corresponding quasinormal mode frequencies. Our results show that interface curvature provides a geometric handle on plasmon propagation and localization, suggesting new possibilities for plasmon transport and analogue gravity experiments with surface plasmons.

We investigate soliton formation and interactions in the non-integrable Gardner-Ostrovsky equation (alias the rotation-modified Gardner equation) with anomalous (positive) dispersion. This equation is relevant to the description of wave dynamics in plasmas and other dispersive media. We show that a special class of solitons, characterized by zero total mass and non-monotonic asymptotics, can emerge from pulse-like initial perturbations. Depending on the initial conditions, these solitons may form regular amplitude-ordered trains, irregular nonstationary ensembles of interacting localized structures, or stationary propagating multi-soliton complexes. We further demonstrate that soliton interactions in the Gardner-Ostrovsky equation are inherently inelastic, leading to the emergence of a dominant "soliton-champion" in closed systems. For example, under periodic boundary conditions, only the largest-amplitude soliton ultimately survives, while all smaller solitons are gradually eliminated through successive interactions. Despite their inelastic interactions, Gardner-Ostrovsky solitons exhibit remarkable robustness.

We provide conditions for existence of hyperbolic, unbounded periodic and elliptic solutions in terms of Weierstrass wp functions of both third and fifth-order KdV--BBM (Korteweg-de Vries--Benjamin, Bona \& Mahony) regularized long wave equation. An analysis for the initial value problem is developed together with a local and global well-posedness theory for the third-order KdV--BBM equation. Traveling wave reduction is used together with zero boundary conditions to yield solitons and periodic unbounded solutions, while for nonzero boundary conditions we find solutions in terms of Weierstrass elliptic wp functions. For the fifth-order KdV--BBM equation we show that a parameter γ=1{12}, for which the equation has a Hamiltonian, represents a restriction for which there are constraint curves that never intersect a region of unbounded solitary waves, which in turn shows that only dark or bright solitons and no unbounded solutions exist. Motivated by the lack of a Hamiltonian structure for γneq1{12} we develop H^k bounds, and we show for the non Hamiltonian system that dark and bright solitons coexist together with unbounded periodic solutions. For nonzero boundary conditions, due to the complexity of the nonlinear algebraic system of coefficients of the elliptic equation we construct Weierstrass solutions for a particular set of parameters only.

We present a model for chi^{(2)} waveguides accounting for three modes, two of which make an avoided crossing at the second harmonic wavelength. We introduce two linearly coupled pure modes and adjust the coupling to replicate the waveguide dispersion near the avoided crossing. Analysis of the nonlinear system reveals continuous wave (CW) solutions across much of the parameter-space and prevalence of its modulational instability. We also predict the existence of the avoided-crossing solitons, and study peculiarities of their dynamics and spectral properties, which include formation of a pedestal in the pulse tails and associated pronounced spectral peaks. Mapping these solitons onto the linear dispersion diagrams, we make connections between their existence and CW existence and stability. We also simulate the two-color soliton generation from a single frequency pump pulse to back up its formation and stability properties.

We present a new derivation of the exact compact oscillon in the signum-Gordon model using Fourier methods. This approach decomposes the oscillon into an infinite, uniformly convergent series of discrete modes consisting of a fundamental frequency and its odd higher harmonics, in contrast to standard analytic potentials, where the corresponding infinite perturbative series fails to converge. Crucially, our oscillon evades the non-existence theorem for spatially localized, small-amplitude periodic states. We investigate the accuracy of truncating the Fourier series and find that the primary mode alone provides a remarkably accurate approximation of the exact solution, with agreement improving significantly when the first two non-vanishing harmonic modes are included. Furthermore, we generalize this framework to a swaying compact oscillon with a periodically moving boundary. We show that describing such an oscillon requires even-order modes, whose coupling with the odd source terms induces a dynamic modulation that generates the swaying motion, even though the spatial Fourier profiles themselves are time-independent.

The nonlinear pseudo-oscillator recently tackled by Gadella and Lara is mapped to an Emden-Fowler (EF) equation that is written as an autonomous two-dimensional ODE system for which we provide the phase-space analysis and the parametric solution. Through an invariant transformation we find periodic solutions to a certain class of EF equations that pass an integrability condition. We show that this condition is necessary to have periodic solutions and via the ODE analysis we also find the sufficient condition for periodic orbits. EF equations that do not pass integrability conditions can be made integrable via an invariant transformation which also allows us to construct periodic solutions to them. Two other nonlinear equations, a zero-frequency Ermakov equation and a positive power Emden-Fowler equation are discussed in the same context

Traveling wave solutions to Kawahara equation (KE), transmission line (TL), and Korteweg-de Vries (KdV) equation are found by using an elliptic function method which is more general than the tanh-method. The method works by assuming that a polynomial ansatz satisfies a Weierstrass equation, and has two advantages: first, it reduces the number of terms in the ansatz by an order of two, and second, it uses Weierstrass functions which satisfy an elliptic equation for the dependent variable instead of the hyperbolic tangent functions which only satisfy the Riccati equation with constant coefficients. When the polynomial ansatz in the traveling wave variable is of first order, the equation reduces to the KdV equation with only a cubic dispersion term, while for the KE which includes a fifth order dispersion term the polynomial ansatz must necessary be of quadratic type. By solving the elliptic equation with coefficients that depend on the boundary conditions, velocity of the traveling waves, nonlinear strength, and dispersion coefficients, in the case of KdV equation we find the well-known solitary waves (solitons) for zero boundary conditions, as well as wave-trains of cnoidal waves for nonzero boundary conditions. Both solutions are either compressive (bright) or rarefactive (dark), and either propagate to the left or right with arbitrary velocity.

In this work, we consider the generalized Benjamin-Bona-Mahony equation $partial_t u+partial_x u+partial_x( |u|^pu)-partial_t partial_x^{2}u=0, quad(t,x) in R times R, with p>4. This equation has the traveling wave solutions \phi_{c}(x-ct), for any frequency c>1. It has been proved by Souganidis and Strauss Strauss-1990 that, there exists a number c_{0}(p)>1, such that solitary waves \phi_{c}(x-ct) with 1<c<c_{0}(p) is orbitally unstable, while for c>c_{0}(p), \phi_{c}(x-ct) is orbitally stable. The linear exponential instability in the former case was further proved by Pego and Weinstein Pego-1991-eigenvalue. In this paper, we prove the orbital instability in the critical case c=c_{0}(p)$.

We present an open-source Python implementation of an idealized high-order pseudo-spectral solver for the one-dimensional nonlinear Schr\"odinger equation (NLSE). The solver combines Fourier spectral spatial discretization with an adaptive eighth-order Dormand-Prince time integration scheme to achieve machine-precision conservation of mass and near-perfect preservation of momentum and energy for smooth solutions. The implementation accurately reproduces fundamental NLSE phenomena including soliton collisions with analytically predicted phase shifts, Akhmediev breather dynamics, and the development of modulation instability from noisy initial conditions. Four canonical test cases validate the numerical scheme: single soliton propagation, two-soliton elastic collision, breather evolution, and noise-seeded modulation instability. The solver employs a 2/3 dealiasing rule with exponential filtering to prevent aliasing errors from the cubic nonlinearity. Statistical analysis using Shannon, R\'enyi, and Tsallis entropies quantifies the spatio-temporal complexity of solutions, while phase space representations reveal the underlying coherence structure. The implementation prioritizes code transparency and educational accessibility over computational performance, providing a valuable pedagogical tool for exploring nonlinear wave dynamics. Complete source code, documentation, and example configurations are freely available, enabling reproducible computational experiments across diverse physical contexts where the NLSE governs wave evolution, including nonlinear optics, Bose-Einstein condensates, and ocean surface waves.

Reid's m'th-order generalized Ermakov systems of nonlinear coupling constant alpha are equivalent to an integrable Emden-Fowler equation. The standard Ermakov-Lewis invariant is discussed from this perspective, and a closed formula for the invariant is obtained for the higher-order Reid systems (m\geq 3). We also discuss the parametric solutions of these systems of equations through the integration of the Emden-Fowler equation and present an example of a dynamical system for which the invariant is equivalent to the total energy

This thesis addresses a fundamental problem in deformation quantization: the difficulty of calculating the star-exponential, the symbol of the evolution operator, due to convergence issues. Inspired by the formalism that connects the star-exponential with the quantum propagator for bosonic systems, this work develops the analogous extension for the fermionic case. A rigorous method, based on Grassmann variables and coherent states, is constructed to obtain a closed-form expression for the fermionic star-exponential from its associated propagator. As a primary application, a fermionic version of the Feynman-Kac formula is derived within this formalism, allowing for the calculation of the ground state energy directly in phase space. Finally, the method is validated by successfully applying it to the simple and driven harmonic oscillators, where it is demonstrated that a simplified ("naive") approach (with an ad-hoc "remediation") is a valid weak-coupling limit of the rigorous ("meticulous") formalism, thereby providing a new and powerful computational tool for the study of fermionic systems.

In this paper we use a traveling wave reduction or a so-called spatial approximation to comprehensively investigate the periodic solutions of the complex cubic-quintic Ginzburg-Landau equation. The primary tools used here are Hopf bifurcation theory and perturbation theory. Explicit results are obtained for the post-bifurcation periodic orbits and their stability. Generalized and degenerate Hopf bifurcations are also briefly considered to track the emergence of global structure such as homoclinic orbits.

It is shown that in the complex trajectory representation of quantum mechanics, the Born's Psi^{\star}Ψprobability density can be obtained from the imaginary part of the velocity field of particles on the real axis. Extending this probability axiom to the complex plane, we first attempt to find a probability density by solving an appropriate conservation equation. The characteristic curves of this conservation equation are found to be the same as the complex paths of particles in the new representation. The boundary condition in this case is that the extended probability density should agree with the quantum probability rule along the real line. For the simple, time-independent, one-dimensional problems worked out here, we find that a conserved probability density can be derived from the velocity field of particles, except in regions where the trajectories were previously suspected to be nonviable. An alternative method to find this probability density in terms of a trajectory integral, which is easier to implement on a computer and useful for single particle solutions, is also presented. Most importantly, we show, by using the complex extension of Schrodinger equation, that the desired conservation equation can be derived from this definition of probability density.

This article deals with the following fractional (p,q)-Choquard equation with exponential growth of the form: $varepsilon^{ps}(-Delta)_{p}^{s}u+varepsilon^{qs}(-Delta)_q^su+ Z(x)(|u|^{p-2}u+|u|^{q-2}u)=varepsilon^{mu-N}[|x|^{-mu}*F(u)]f(u) in R^N, where s\in (0,1), \varepsilon>0 is a parameter, 2\leq p=N{s}<q, and 0<\mu<N. The nonlinear function f has an exponential growth at infinity and the continuous potential function Z satisfies suitable natural conditions. With the help of the Ljusternik-Schnirelmann category theory and variational methods, the multiplicity and concentration of positive solutions are obtained for \varepsilon>0$ small enough. In a certain sense, we generalize some previously known results.

In this paper, we study the long time dynamics of small perturbations of static domain walls in weighted Sobolev spaces for the one dimensional Schrödinger map with constant Dzyaloshinskii-Moriya interaction. The perturbation equation contains derivative nonlinearities, and its linearized operator has a threshold resonance. The results are twofold. (i) We prove the asymptotic stability of the domain wall modulo the translation and rotation symmetries. More precisely, the corresponding modulation parameters converge as time tends to infinity, while the radiation decays at the sharp rate t^{-1/2}. (ii) We prove that the radiation exhibits modified scattering, with an explicit logarithmic phase correction generated by the long range cubic interaction. The proof develops geometric reduction to overcome the loss of derivatives and exploits structural cancellations in the nonlinear interactions.

By using ZEUS cluster at Embry-Riddle Aeronautical University we perform extensive numerical simulations based on a two-dimensional Fourier spectral method Fourier spatial discretization and an explicit scheme for time differencing) to find the range of existence of the spatiotemporal solitons of the two-dimensional complex Ginzburg-Landau equation with cubic and quintic nonlinearities. We start from the parameters used by Akhmediev {\it et. al.} and slowly vary them one by one to determine the regimes where solitons exist as stable/unstable structures. We present eight classes of dissipative solitons from which six are known (stationary, pulsating, vortex spinning, filament, exploding, creeping) and two are novel (creeping-vortex propellers and spinning "bean-shaped" solitons). By running lengthy simulations for the different parameters of the equation, we find ranges of existence of stable structures (stationary, pulsating, circular vortex spinning, organized exploding), and unstable structures (elliptic vortex spinning that leads to filament, disorganized exploding, creeping). Moreover, by varying even the two initial conditions together with vorticity, we find a richer behavior in the form of creeping-vortex propellers, and spinning "bean-shaped" solitons. Each class differentiates from the other by distinctive features of their energy evolution, shape of initial conditions, as well as domain of existence of parameters.

Singularity Theory is used to comprehensively investigate the bifurcations of the steady-states of the traveling wave ODEs of the cubic-quintic Ginzburg-Landau equation (CGLE). These correspond to plane waves of the PDE. In addition to the most general situation, we also derive the degeneracy conditions on the eight coefficients of the CGLE under which the equation for the steady states assumes each of the possible quartic (the quartic fold and an unnamed form), cubic (the pitchfork and the winged cusp), and quadratic (four possible cases) normal forms for singularities of codimension up to three. Since the actual governing equations are employed, all results are globally valid, and not just of local applicability. In each case, the recognition problem for the unfolded singularity is treated. The transition varieties, i.e. the hysteresis, isola, and double limit curves are presented for each normal form. For both the most general case, as well as for various combinations of coefficients relevant to the particular cases, the bifurcations curves are mapped out in the various regions of parameter space delimited by these varieties. The multiplicities and interactions of the plane wave solutions are then comprehensively deduced from the bifurcation plots in each regime, and include features such as regimes of hysteresis among co-existing states, domains featuring more than one interval of hysteresis, and isola behavior featuring dynamics unrelated to the primary solution branch in limited ranges of parameter space.

In this paper, Liouville soliton surfaces based on some soliton solutions of the Liouville equation are constructed and displayed graphically, including some of those corresponding to Darboux-transformed counterparts. We find that the Liouville soliton surfaces are centroaffine surfaces of Tzitzeica type and their centroaffine invariant can be expressed in terms of the Hamiltonian. The traveling wave solutions to Liouville equation from which these soliton surfaces stem are also obtained through a modified variation of parameters method which is shown to lead to elliptic functions solution method.

We establish the nonlinear stability of solitary waves (solitons) and periodic traveling wave solutions (cnoidal waves) for a Korteweg-de Vries (KdV) equation which includes a fifth order dispersive term. The traveling wave solutions which yield solitons for zero boundary conditions and wave-trains of cnoidal waves for nonzero boundary conditions are analyzed using stability theorems, which rely on the positivity properties of the Fourier transforms. We show that all families of solutions considered here are (orbitally) stable.

We investigate the KdV-Burgers and Gardner equations with dissipation and external perturbation terms by the approach of dynamical systems and Shil'nikov's analysis. The stability of the equilibrium point is considered, and Hopf bifurcations are investigated after a certain scaling that reduces the parameter space of a three-mode dynamical system which now depends only on two parameters. The Hopf curve divides the two-dimensional space into two regions. On the left region the equilibrium point is stable leading to dissapative periodic orbits. While changing the bifurcation parameter given by the velocity of the traveling waves, the equilibrium point becomes unstable and a unique stable limit cycle bifurcates from the origin. This limit cycle is the result of a supercritical Hopf bifurcation which is proved using the Lyapunov coefficient together with the Routh-Hurwitz criterion. On the right side of the Hopf curve, in the case of the KdV-Burgers, we find homoclinic chaos by using Shil'nikov's theorem which requires the construction of a homoclinic orbit, while for the Gardner equation the supercritical Hopf bifurcation leads only to a stable periodic orbit.

This paper discusses the use of the Riemann-Silberstein vector to solve the source-free Maxwell's equations and obtains novel analytical solutions. The solving process naturally leads to the spinor form of the source-free Maxwell's equations. Several powerful theorems are established to solve this spinor form equation. The Waveguide Solution Theorem provides an elegant way to solve waveguide problems, while The Schrodinger Solution Theorem connects the Maxwell's equations with the two-dimensional Schrodinger equation. By utilizing The Schrodinger Solution Theorem, a precise formula for spatiotemporal diffraction of the Maxwell's equations is derived, which allows for the reconstruction of electromagnetic waves throughout space and time based on the field distribution on the diffraction screen. And finally, by studying some relevant contour integrals, two additional simple but beautiful mathematical theorems that are necessarily satisfied by electromagnetic field in any source-free region are derived.

In this paper, we propose an alternating numerical scheme for constructing quasi-periodic solutions to a class of nonlinear partial differential equations (PDEs) with bounded perturbations, including the one-dimensional nonlinear Schrödinger (NLS) equations and nonlinear wave (NLW) equations under periodic and Dirichlet boundary conditions. A key observation concerning the resonant Q-equations allows us to establish directly a local diffeomorphism between the amplitudes and the drifted frequencies. This renders the proof considerably cleaner, since it avoids the Birkhoff normal form and the associated coordinate transformations required in classical KAM theory for PDEs.Unlike in the finite-dimensional setting, the lattice boxes must also include the spatial index, so that the multiscale induction has to be carried out in one additional dimension. This introduces two difficulties. First, the measure estimate for the excluded set of ``bad'' frequencies must account for the spatial index. For both the NLS and NLW equations, the linear frequencies and their differences are integers or close to integers; this allows the spatial dependence to be controlled through the modulus of the temporal index, so that the excluded measure remains under control. Second, a singular set may contain sites with the same temporal index but widely separated spatial indices, which obstructs the ``inversion implies localization'' process. For the one-dimensional NLS and NLW equations, however, such sites occur only in pairs with opposite spatial indices. They can therefore be separated by small boxes, after which the global inverse is obtained from the local inverses via the resolvent identity. Numerical experiments illustrate the scheme for nonlinear traveling and standing waves, exhibiting both the computed quasi-periodic solutions and the convergence of the iteration.

We show that some Emden-Fowler (EF) equations encountered in astrophysics and cosmology belong to two EF integrable classes of the type d^2z/dχ^2=A χ^{-λ-2}z^n for λ=(n-1)/2 (class one), and λ=n+1 (class two). We find their corresponding invariants which reduce them to first order nonlinear ordinary differential equations. Using particular solutions of such EF equations, the two classes are set in the autonomous nonlinear oscillator form d^2ν/dt^2 +ad ν/dt +b(ν-ν^n)=0, where the coefficients a,b depend only on λ, n. For both classes, we write closed-form solutions in parametric form. The illustrative examples from astrophysics and general relativity correspond to two n=2 cases from class one and two, and one n=5 case from class one, all of them yielding Weierstrass elliptic solutions. It is also noticed that when n=2, EF equations can be studied using the Painlevé reduction method, since they are a particular case of equations of the type d^2z/dχ^2=F(χ)z^2, where F(χ) is the Kustaanheimo-Qvist function

In this paper, we prove the existence of stationary pulse solutions for a bulk-surface reaction-diffusion system modeling cell polarization in two-dimensional domains. Using a singular perturbation method, we construct pulse solutions whose centers lie near non-degenerate critical points of a geometry-dependent potential expressed in terms of the Neumann Green's function, as formally predicted in watanabe2026. The main analytical difficulty arises from the half-Laplacian induced by the logarithmic singularity of the Green's function, resulting in algebraically rather than exponentially decaying inner solutions. Our analysis may also serve as a useful guide for applying singular perturbation methods to other bulk-surface type systems.

We derive a quantum master equation in Lindblad form for the open quantum Rabi model in the perturbative deep strong coupling regime, where the spectrum is close to degenerate and harmonic. The Lindblad form is obtained by combining secular and singular-coupling approximations, working in a Schrieffer-Wolff transformed frame. We benchmark our equation against a Bloch-Redfield one, analyzing the time-evolution and steady state properties of several observables and find excellent agreement. Analytical results are obtained in limiting cases.

Using a variational formulation for partial differential equations (PDEs) combined with numerical simulations on ordinary differential equations (ODEs), we find two categories (pulses and snakes) of dissipative solitons, and analyze the dependence of both their shape and stability on the physical parameters of the cubic-quintic Ginzburg-Landau equation (CGLE). In contrast to the regular solitary waves investigated in numerous integrable and non-integrable systems over the last three decades, these dissipative solitons are not stationary in time. Rather, they are spatially confined pulse-type structures whose envelopes exhibit complicated temporal dynamics. Numerical simulations reveal very interesting bifurcations sequences as the parameters of the CGLE are varied. Our predictions on the variation of the soliton amplitude, width, position, speed and phase of the solutions using the variational formulation agree with simulation results.

We study four-point correlation functions of degenerated fields in the NS sector in Super-Liouville field theory. We find integral expressions for these functions using the BPZ equation, and study some superconformal properties of these solutions. Finally, we present the general form for three-point correlation functions.

The transition of an impulsively excited kink oscillation of a solar coronal loop to an oscillation with a stationary amplitude, i.e., the damping pattern, is determined using the low-dimensional self-oscillation model. In the model, the decayless kink oscillations are sustained by the interaction of the oscillating loop with an external quasi-steady flow. The analytical solution is based on the assumption that the combined effect of the effective dissipation, for example, by resonant absorption, and interaction with an external flow, is weak. The effect is characterised by a dimensionless coupling parameter. The damping pattern is found to depend upon the initial amplitude and the coupling parameter. The approximate expression shows a good agreement with a numerical solution of the self-oscillation equation. The plausibility of the established damping pattern is demonstrated by an observational example. Notably, the damping pattern is not exponential, and the characteristic decay time is different from the time determined by the traditionally used exponential damping fit. Implications of this finding for seismology of the solar coronal plasmas are discussed. In particular, it is suggested that a very rapid, in less than the oscillation period, decay of the oscillation to the stationary level, achieved for larger values of the coupling parameter, can explain the relative rareness of the kink oscillation events.

We introduce spatiotemporal solitons of the two-dimensional complex Ginzburg-Landau equation (2D CCQGLE) with cubic and quintic nonlinearities in which asymmetry between space-time variables is included. The 2D CCQGLE is solved by a powerful Fourier spectral method, i.e., a Fourier spatial discretization and an explicit scheme for time differencing. Varying the system's parameters, and using different initial conditions, numerical simulations reveal 2D solitons in the form of stationary, pulsating and exploding solitons which possess very distinctive properties. For certain regions of parameters, we have also found stable coherent structures in the form of spinning (vortex) solitons which exist as a result of a competition between focusing nonlinearities and spreading while propagating through medium.

This work presents a non-perturbative quantum field theoretic computation of spinor tunnelling through a potential barrier. Within position space potential formalism, we identify coupled recursion relations forming a self-consistent integral equation. Our method, enabling both analytical and numerical approaches, transforms into a formalism admitting both numerical framework for arbitrary potentials and analytical exact results for select barriers. For the rectangular barrier, the QFT-derived amplitudes reproduce the relativistic quantum mechanical (RQM) results, thus establishing a foundation for future calculation of beyond-RQM field-theoretic corrections.

Earlier, in the CINT program at Los Alamos National Laboratory, we focused ultrafast mode-locked lasers on the tip-sample junction of a scanning tunneling microscope to generate currents at hundreds of harmonics of the laser pulse repetition frequency. Each harmonic has a signal-to-noise ratio of 20 dB with a 10-dB linewidth of only 3 Hz. Now we model closed quantum nanocircuits with rectangular, triangular, or delta-function barrier, shunted by a beryllium filament for quasi-coherent electron transport over mean-free paths as great as 68 nm. The time-independent Schrödinger equation is solved with the boundary conditions that the wavefunction and its derivative are continuous at both connections. These four boundary conditions are used to form a four-by-four complex matrix equation with only zeros in the right-hand column vector which is required to have a non-trivial solution with each of the closed nanocircuits. Each model has four parameters: (1) the barrier length, (2) the height and shape of the barrier, (3) the length of the pre-barrier, and (4) the electron energy. Any three of these may be specified and then the fourth is varied to bring the determinant to zero to find the solutions on lines or surfaces in the space defined by the four parameters. First, we use a simplistic model having a rectangular barrier. The second model has a triangular barrier as a first approximation to field emission, and we are considering applying this approach for a self-contained nanoscale extension of our earlier effort to generate the harmonics at Los Alamos. The third model has a delta-function barrier, and the fourth model is an extension of the first one where the width of the rectangular barrier is varied inversely with its height.

We prove that all excited states -sin(nx), ngeq2, of the De Gregorio equation on the torus are nonlinearly unstable in the sense of Hadamard. More precisely, perturbations arbitrarily small in H^2 leave a fixed H^1 neighborhood of the state. For every ngeq2, the linearized operator has a real eigenvalue λ_nin(0.56,0.6). To prove this, we transform the eigenvalue problem into a Fuchsian equation on the complex disk and obtain a matching condition for the eigenvalue. This condition is verified rigorously using interval arithmetic. The nonlinear instability follows from a bootstrap argument in the spirit of Guo, Hallstrom, and Spirn.

In this paper, we investigate the quantum behaviour near the classical singularity of a black hole. We describe the black-hole interior by the Kantowski-Sachs metric, with a string propagating in this region. This approach is motivated by the fact that, near the classical singularity, spacetime curvature reaches Planck or string scales, necessitating a description in terms of string-theoretic degrees of freedom. The string dynamics is formulated via the Polyakov action, and the combined system is canonically quantized to derive the corresponding Wheeler-DeWitt equation. Although the equation is separable in the radial and azimuthal string sectors, the remaining geometric and string degrees of freedom are nontrivially coupled. To address this, we employ a perturbative approach and obtain the wave function up to first order. At the zeroth-order, the angular part is described by the general Heun function, while the geometric part is expressed in terms of the Bessel function. The first-order wave function is governed by a source term containing the zeroth-order solution, and we obtain its explicit expression in terms of the complementary function and the particular integral. We confirm regularity of the wave function at both zeroth and first order, at the singularity as well as at the horizon. Regularity at the horizon imposes a discretization condition on one eigenvalue, while another eigenvalue is quantized due to the azimuthal periodicity of the wave function. Importantly, both zeroth- and first-order solutions satisfy the DeWitt boundary condition, implying resolution of the black hole singularity within this quantum framework.

The two-body problem of variational electrodynamics possesses differential-delay equations of motion with state-dependent delays of neutral type and solutions that can have velocity discontinuities on countable sets. From a periodic orbit possessing some mild properties at breaking points, we define a synchronization function in R x R3, which is further used to construct two bounded oscillatory functions vanishing at breaking points and whose first derivatives are continuous and defined everywhere but at breaking points. The oscillatory functions are associated with a PDE identity in H2(R3), and we postulate ordering conditions for the PDE identity to define a Fredholm-Schroedinger operator with an O(1/r2) spin-orbit forcing term belonging to L2(R3). As an application, we introduce the Chemical Principle criterion to select orbits with asymptotically vanishing far-fields and estimate the Bohr radius parameter of the Fredholm-Schroedinger PDE using the boundary-layer of orbits chosen according to the Chemical Principle criterion. Last, working backward, we derive a Chemical Principle-like condition from the ordering conditions.

Schrödinger operators with nonlocal δ-interaction are studied with the use of the Lax-Phillips scattering theory methods. The condition of applicability of the Lax-Phillips approach in terms of non-cyclic functions is established. Two formulas for the S-matrix are obtained. The first one deals with the Krein-Naimark resolvent formula and the Weyl-Titchmarsh function, whereas the second one is based on modified reflection and transmission coefficients. The S-matrix S(z) is analytical in the lower half-plane C_- when the Schrödinger operator with nonlocal δ-interaction is positive self-adjoint. Otherwise, S(z) is a meromorphic matrix-valued function in C_- and its properties are closely related to the properties of the corresponding Schrödinger operator. Examples of S-matrices are given.

Let G be a Lie group equipped with commuting left and right linear representations on a vector space V, and let GJoin V be the associated centered semidirect product. We develop a reduction scheme for covariant field theories on principal bundles with structure group GJoin V from both Lagrangian and Hamiltonian perspectives, obtaining the reduced Euler-Poincaré and Lie-Poisson field equations decomposed into their mathfrak g- and V-components, where mathfrak g is the Lie algebra of G. We also derive the reduced multisymplectic form formula, as well as the Noether and energy-momentum balance laws. Our reduction results are compared with those for principal bundles with standard semidirect-product structure groups, and the role of the centered variables is investigated: they retain the original two-sided geometric data and separate the left and right contributions via the heart and diamond operators. The framework is then applied to a two-sided matrix microstructural model and to the second-order jet group on mathbb R^r, G_r^2. For the matrix model, the conjugation representation decomposes Mat(3) into scalar, vector, and symmetric-traceless sectors, clarifying its precise relation to director and alignment-tensor variables without identifying it with standard nematodynamics. For the one-dimensional jet group on a (1+1)-dimensional Minkowski base, we obtain a semilinear symmetric-hyperbolic system and an exact traveling conversion front with explicit energy-momentum exchange between the linear- and quadratic-jet sectors. Its hyperbolic-target interpretation yields a conserved positive energy for the linearized perturbations. A nonlinear numerical experiment illustrates the reduced evolution and reconstruction of the group-valued field, with consistency assessed through mesh refinement and an independent wave-map computation.

Despite the recent advances in the field of computational Schr\"odinger Bridges (SB), most existing SB solvers are still heavy-weighted and require complex optimization of several neural networks. It turns out that there is no principal solver which plays the role of simple-yet-effective baseline for SB just like, e.g., k-means method in clustering, logistic regression in classification or Sinkhorn algorithm in discrete optimal transport. We address this issue and propose a novel fast and simple SB solver. Our development is a smart combination of two ideas which recently appeared in the field: (a) parameterization of the Schr\"odinger potentials with sum-exp quadratic functions and (b) viewing the log-Schr\"odinger potentials as the energy functions. We show that combined together these ideas yield a lightweight, simulation-free and theoretically justified SB solver with a simple straightforward optimization objective. As a result, it allows solving SB in moderate dimensions in a matter of minutes on CPU without a painful hyperparameter selection. Our light solver resembles the Gaussian mixture model which is widely used for density estimation. Inspired by this similarity, we also prove an important theoretical result showing that our light solver is a universal approximator of SBs. Furthemore, we conduct the analysis of the generalization error of our light solver. The code for our solver can be found at https://github.com/ngushchin/LightSB

We obtain (small-parameter) well-posedness for the (space-time periodic) Phi^4 equation in the full subcritical regime in the context of regularity structures based on multi-indices. As opposed to Hairer's more extrinsic tree-based setting, due to the intrinsic description encoded by multi-indices, it is not possible to obtain a solution theory via the standard fixed-point argument. Instead, we develop a more intrinsic approach for existence using a variant of the continuity method from classical PDE theory based on a priori estimates for a new `robust' formulation of the equation. This formulation also allows us to obtain uniqueness of solutions and continuity of the solution map in the model norm even at the limit of vanishing regularisation scale. Since our proof relies on the structure of the nonlinearity in only a mild way, we expect the same ideas to be sufficient to treat a more general class of equations.

The Korteweg-de Vries (KdV) equation serves as a foundational model in nonlinear wave physics, describing the balance between dispersive spreading and nonlinear steepening that gives rise to solitons. This article introduces sangkuriang, an open-source Python library for solving this equation using Fourier pseudo-spectral spatial discretization coupled with adaptive high-order time integration. The implementation leverages just-in-time (JIT) compilation for computational efficiency while maintaining accessibility for instructional purposes. Validation encompasses progressively complex scenarios including isolated soliton propagation, symmetric two-wave configurations, overtaking collisions between waves of differing amplitudes, and three-body interactions. Conservation of the classical invariants is monitored throughout, with deviations remaining small across all test cases. Measured soliton velocities conform closely to theoretical predictions based on the amplitude-velocity relationship characteristic of integrable systems. Complementary diagnostics drawn from information theory and recurrence analysis confirm that computed solutions preserve the regular phase-space structure expected for completely integrable dynamics. The solver outputs data in standard scientific formats compatible with common analysis tools and generates visualizations of spatiotemporal wave evolution. By combining numerical accuracy with practical accessibility on modest computational resources, sangkuriang offers a platform suitable for both classroom demonstrations of nonlinear wave phenomena and exploratory research into soliton dynamics.

We consider Proca field perturbations in a five-dimensional Schwarzschild-anti-de Sitter (Schwarzschild-AdS_{5}) black hole geometry. Using the vector spherical harmonic (VSH) method, we show that the Proca field decomposes into scalar-type and vector-type components according to their tensorial behavior on the three-sphere. Two degrees of freedom of the field are described by scalar-type components, which are coupled due to the mass term, while the remaining two degrees of freedom are described by a vector-type component, which decouples completely. Motivated by the Frolov-Krtouš-Kubizňák-Santos (FKKS) ansatz in the limit of zero spin, we use a field transformation to decouple the scalar-type components at the expense of introducing a complex separation parameter β. This parameter can be determined analytically, and its values correspond to two distinct polarizations of the scalar-type sector: "electromagnetic" and "non-electromagnetic", denoted by β_{+} and β_{-}, respectively. In the scalar-type sector, the radial differential equation for each polarization is a Fuchsian differential equation with five singularities, whereas in the vector-type sector, the radial equation has four singularities. By means of the isomonodromy method, we reformulate the boundary value problem in terms of the initial conditions of the Painlevé VI τ function and, using a series expansion of the τ function, we compute the scalar-type and vector-type quasinormal modes (QNMs) in the small horizon limit. Our results are in overall very good agreement with those obtained via the numerical integration method. This shows that the isomonodromy method is a reliable method to compute quasinormal modes in the small horizon limit with high accuracy.

An implicit Euler finite-volume scheme for a nonlocal cross-diffusion system on the multidimensional torus is analyzed. The equations describe the dynamics of population species with repulsive or attractive interactions. The numerical scheme is based on a generalized Scharfetter-Gummel discretization of the nonlocal flux term. For merely integrable kernel functions, the scheme preserves the positivity, total mass, and entropy structure. The existence of a discrete solution and its convergence to a solution to the continuous problem, as the mesh size tends to zero, are shown. A key difficulty is the degeneracy of the generalized Bernoulli function in the Scharfetter-Gummel approximation. This issue is overcome by proving a uniform estimate for the discrete Fisher information, which requires both the Boltzmann and Rao entropy inequalities. Numerical simulations illustrate the features of the scheme in one and two space dimensions.

We study the exact controllability of the cubic nonlinear Schrödinger equation in the exterior Ω=R^nsetminusΘ of a non-trapping obstacle Θ, with internal controls supported outside a large ball, and the corresponding boundary control problem on Ω_0=B_{R_0}setminusΘ with Dirichlet controls acting only on the outer sphere partial B_{R_0}. Using the local smoothing effect of Burq, Gérard and Tzvetkov [9], we prove observability inequalities for the linear Schrödinger equation in the whole scale of Sobolev spaces H^σ_D(Ω), σin[-2,2], associated with the Dirichlet Laplacian; they require only the non-trapping assumption (no star-shapedness of Θ), and their proofs use no normal boundary traces. Combined with Strichartz-type estimates and a perturbation argument, this yields local exact controllability in H^1_0(Ω) and in H^2(Ω)cap H^1_0(Ω), in dimensions two and three, around the zero solution and, for internal controls and under a unique continuation assumption, around nontrivial trajectories. The same method applies to the quintic equation, which is energy-critical in dimension three. Several open problems are discussed.

We study the parameter dependence and transversality properties of trajectories of the controlled ODE $ U_x = P, qquad P_x = - f(U) - βP + tilde α(x)U, arising from traveling-wave solutions of a bistable reaction-diffusion equation with control. Optimal profiles minimizing the L^1-norm of the control \tilde α are determined by the unstable and stable manifolds Γ_u, Γ_s, the critical curve P^*(U) = Uf(U), and the associated function F(U) = dP^*{dU} + f(U){P^*(U)}. We establish generic transversality properties for these objects with respect to perturbations of the speed β and the nonlinear source f, together with quantitative estimates on the first- and second-order dependence of the invariant manifolds on these parameters. We also analyze the loss of regularity at the critical speed β^{**}, where the unstable manifold approaches the strongly stable manifold of the equilibrium (U^*,0)$. These results provide the ODE ingredients used in the analysis of the structure and regularity of optimal profiles in paper [S. Bianchini and C. Trifone. Traveling profiles and control cost for a PDE describing the evolution of invasive species. 2026. arXiv: 2606.23979 [math.OC]].

We present an analysis of the Rayleigh-Plesset equation for a three dimensional vacuous bubble in water. In the simplest case when the effects of surface tension are neglected, the known parametric solutions for the radius and time evolution of the bubble in terms of a hypergeometric function are briefly reviewed. By including the surface tension, we show the connection between the Rayleigh-Plesset equation and Abel's equation, and obtain the parametric rational Weierstrass periodic solutions following the Abel route. In the same Abel approach, we also provide a discussion of the nonintegrable case of nonzero viscosity for which we perform a numerical integration

Numerical simulations of the complex cubic-quintic Ginzburg-Landau equation (CCQGLE), a canonical equation governing the weakly nonlinear behavior of dissipative systems in a wide variety of disciplines, reveal five entirely novel classes of pulse or solitary waves solutions, viz. pulsating, creeping, snaking, erupting, and chaotical solitons. Here, we develop a theoretical framework for analyzing the full spatio-temporal structure of one class of dissipative solution (snaking soliton) of the CCQGLE using the variational approximation technique and the dynamical systems theory. The qualitative behavior of the snaking soliton is investigated using the numerical simulations of (a) the full nonlinear complex partial differential equation and (b) a system of three ordinary differential equations resulting from the variational approximation.

We study the Dirichlet problem for the weighted Schr\"odinger operator \[-\Delta u +Vu = \lambda \rho u,\] where rho is a positive weighting function and V is a potential. Such equations appear naturally in conformal geometry and in the composite membrane problem. Our primary goal is to establish concavity estimates for the principle eigenfunction with respect to conformal connections. Doing so, we obtain new bounds on the fundamental gap problem, which is the difference between the first and second eigenvalues. In particular, we partially resolve a conjecture of Nguyen, Stancu and Wei [IMRN 2022] on the fundamental gap of horoconvex domains. In addition, we obtain a power convexity estimate for solutions to the torsion problem in spherical geometry on convex domains which are not too large.

We consider using deep neural networks to solve time-dependent partial differential equations (PDEs), where multi-scale processing is crucial for modeling complex, time-evolving dynamics. While the U-Net architecture with skip connections is commonly used by prior studies to enable multi-scale processing, our analysis shows that the need for features to evolve across layers results in temporally misaligned features in skip connections, which limits the model's performance. To address this limitation, we propose SineNet, consisting of multiple sequentially connected U-shaped network blocks, referred to as waves. In SineNet, high-resolution features are evolved progressively through multiple stages, thereby reducing the amount of misalignment within each stage. We furthermore analyze the role of skip connections in enabling both parallel and sequential processing of multi-scale information. Our method is rigorously tested on multiple PDE datasets, including the Navier-Stokes equations and shallow water equations, showcasing the advantages of our proposed approach over conventional U-Nets with a comparable parameter budget. We further demonstrate that increasing the number of waves in SineNet while maintaining the same number of parameters leads to a monotonically improved performance. The results highlight the effectiveness of SineNet and the potential of our approach in advancing the state-of-the-art in neural PDE solver design. Our code is available as part of AIRS (https://github.com/divelab/AIRS).

In this paper, a PDE-ODE model with discontinuity in the flux as well as a flux constraint is analyzed. A modified Riemann solution is proposed and the existence of a weak solution to the Cauchy problem is rigorously investigated using the wavefront tracking scheme.

We find and analyse solutions of Einstein's equations in arbitrary d dimensions and in the presence of a scalar field with a Liouville potential coupled to a Maxwell field. We consider spacetimes of cylindrical symmetry or again subspaces of dimension d-2 with constant curvature and analyse in detail the field equations and manifest their symmetries. The field equations of the full system are shown to reduce to a single or couple of ODE's which can be used to solve analytically or numerically the theory for the symmetry at hand. Further solutions can also be generated by a solution generating technique akin to the EM duality in the absence of a cosmological constant. We then find and analyse explicit solutions including black holes and gravitating solitons for the case of four dimensional relativity and the higher-dimensional oxydised 5-dimensional spacetime. The general solution is obtained for a certain relation between couplings in the case of cylindrical symmetry.

Using the supersymmetric method of random matrix theory within the Heidelberg approach framework we provide statistical description of stationary intensity sampled in locations inside an open wave-chaotic cavity, assuming that the time-reversal invariance inside the cavity is fully broken. In particular, we show that when incoming waves are fed via a finite number M of open channels the probability density {cal P}(I) for the single-point intensity I decays as a power law for large intensities: {cal P}(I)sim I^{-(M+2)}, provided there is no internal losses. This behaviour is in marked difference with the Rayleigh law {cal P}(I)sim exp(-I/I) which turns out to be valid only in the limit Mto infty. We also find the joint probability density of intensities I_1, ldots, I_L in L>1 observation points, and then extract the corresponding statistics for the maximal intensity in the observation pattern. For Lto infty the resulting limiting extreme value statistics (EVS) turns out to be different from the classical EVS distributions.

The eight-component relativistic wave equation for spin-1{2} particles of Robson and Staudte (FV1{2}), a Feshbach--Villars linearization of the Feynman--Gell-Mann equation, reproduces the Dirac hydrogenic spectrum with a doubled solution space and an indefinite inner product. Its exact Stark shifts must equal the Dirac ones; whether its indefinite-metric perturbation theory delivers them is not guaranteed. We compute the first- and second-order Stark shifts of the n=1--3 levels of hydrogenic ions from H to U^{91+} in the Schrödinger, Dirac and FV1{2} theories, evaluating all spectral sums, continuum included, by the Dalgarno--Lewis method. The FV1{2} and Dirac shifts coincide level by level, to working precision at first order (closed-form roots) and to 45 significant figures at second order. Both ingredients are needed: without the explicit spin--field coupling the first-order roots become complex, and without the negative-norm states the (Zα)^2 coefficient is wrong. Norm sign and parity are locked in the doubled space: a potential lifting the 2s_{1/2}--2p_{1/2} degeneracy needs its own spin--gradient term, and shifts placed by hand on the κ-labeled states give a spurious imaginary linear Stark shift. The relativistic linear Stark effect is smaller by finite factors, but a second-moment sum rule shows the non-relativistic linear Stark strength is conserved, the missing part reappearing in singular (Zα)^{-2} second-order terms. For levels with a fine-structure partner, an exact intra-shell/regular split of the second-order shift gives the (Zα)^0 coefficients analytically and relates the π^2 terms of the (Zα)^2 coefficients to the singular ones. The Dirac weak/strong-field crossover for n=2 of hydrogen lies near 2.8~kV cm^{-1}, and the perturbative window set by the fine structure grows as Z^5.

We study the bi-parameter local linearization of the one-dimensional nonlinear stochastic wave equation driven by a Gaussian noise, which is white in time and has a spatially homogeneous covariance structure of Riesz-kernel type. We establish that the second-order increments of the solution can be approximated by those of the corresponding linearized wave equation, modulated by the diffusion coefficient. These findings extend the previous results of Huang et al. HOO2024, which addressed the case of space-time white noise. As applications, we analyze the quadratic variation of the solution and construct a consistent estimator for the diffusion parameter.

We address effects of spin-orbit coupling (SOC), phenomenologically added to a two-component Bose-Einstein condensate composed of particles moving by Levy flights, in one- and two-dimensional (1D and 2D) settings. The corresponding system of coupled Gross-Pitaevskii equations includes fractional kinetic-energy operators, characterized by the Levy index, \alpha < 2 (the normal kinetic energy corresponds to \alpha = 2). The SOC terms, with strength \lambda, produce strong effects in the 2D case: they create families of stable solitons of the semi-vortex (SV) and mixed-mode (MM) types in the interval of 1 < \alpha < 2, where the supercritical collapse does not admit the existence of stable solitons in the absence of the SOC. At \lambda --> 0, amplitudes of these solitons vanish as (\lambda)^{1/(\alpha - 1)}.

This paper presents the first numerical realization of Suzuki's Weil-Quadratic-Form operator, a candidate for the Hilbert--Pólya program linking spectral positivity to the Riemann Hypothesis (RH). Suzuki's 2026 construction was purely theoretical; here, the operator is instantiated via P1 finite-element discretization and Richardson extrapolation. Key results include: (R1) In the prime-free regime, the spectrum follows a closed Archimedean law A_k(a) = log(1/a) + log(k-2) + B_0 + O(a), with B_0 = log q - 2log 2, confirmed to 30-digit precision. (R2) A Mellin double-pole argument proves the head coefficient B(ν) and shows B_0 depends only on the conductor q, independent of the Archimedean parameter. (R2b) The degree d of an L-function appears directly as the logarithmic slope of the spectrum. (R3) Total spectral intensity follows the prime number theorem, S(a) sim (2a)^3/6. (R4) Nontrivial zeros are not eigenvalues but occur in the explicit-formula error term of the prime symbol. (R5) The best-match line σ^*(a) descends toward the critical line. (R6) Weil's positivity criterion is realized in operator form: bounded residual growth corresponds to all zeros on the line, while an injected off-line zero causes exponential blow-up. (R7) The lowest eigenvalue λ_1(a) is strictly positive, decays superexponentially, and passes smoothly through the first prime threshold. (R8) The characteristic function W(a,0;z) is computed for the first time, with all zeros confirmed real. (R9) Indirect traces of GUE statistics appear in the moment structure, even where direct detection is blocked. The authors emphasize that this work does not prove RH. All results are Archimedean and universal, with significance lying in the faithful numerical realization of classical identities rather than new arithmetic.

Starting from the Mellin-Barnes integral representation of a Feynman integral depending on set of kinematic variables z_i, we derive a system of partial differential equations w.r.t.\ new variables x_j, which parameterize the differentiable constraints z_i=y_i(x_j). In our algorithm, the powers of propagators can be considered as arbitrary parameters. Our algorithm can also be used for the reduction of multiple hypergeometric sums to sums of lower dimension, finding special values and reduction equations of hypergeometric functions in a singular locus of continuous variables, or finding systems of partial differential equations for master integrals with arbitrary powers of propagators. As an illustration, we produce a differential equation of fourth order in one variable for the one-loop two-point Feynman diagram with two different masses and arbitrary propagator powers.

The large-D membrane paradigm provides a perturbative, non-gravitational description of black hole dynamics. While the local second law of thermodynamics has been established for neutral membranes, the explicit determination of the local entropy production equation within the charged large-D membrane framework remains unsolved. In this paper, we compute the subleading-order correction to the divergence of the membrane velocity for a charged black hole in a flat background. We evaluate the corresponding scalar Einstein-Maxwell constraint equation at the physical event horizon via the non-affine null Raychaudhuri equation, and relate the physical horizon expressions to the auxiliary membrane kinematics, demonstrating that the membrane velocity produces a non-negative divergence at the first non-vanishing order in the membrane worldvolume, which is driven by two quadratic forms: the membrane shear and a novel contribution generated purely by the Maxwell field. Using this result, we construct the membrane entropy current and prove that its divergence is non-negative through this order, thereby establishing the local second law for the charged large-D membrane paradigm.

We study the dynamics of an inextensible, closed interface subject to bending forces and immersed in a two-dimensional and incompressible Stokes fluid. We formulate the problem as a boundary integral equation in terms of the tangent angle and demonstrate the well-posedness in suitable time-weighted spaces of the resulting nonlinear and nonlocal system. The solution is furthermore shown to be smooth for positive times. Numerical computations are performed to initiate the study of the long-time behavior of the interface.

We obtain closed-form solutions of several inhomogeneous Lienard equations by the factorization method. The two factorization conditions involved in the method are turned into a system of first-order differential equations containing the forcing term. In this way, one can find the forcing terms that lead to integrable cases. Because of the reduction of order feature of factorization, the solutions are simultaneously solutions of first-order differential equations with polynomial nonlinearities. The illustrative examples of Lienard solutions obtained in this way generically have rational parts, and consequently display singularities.

In this paper, both the pointwise space-time structures of the solution for a three-dimensional model describing gas dynamics in thermal non-equilibrium and the corresponding background stabilities within the L^p framework are investigated under the small initial perturbation assumption and the irrotationality condition. Our main tools are the modified Green's function according to the irrotationality condition, a priori energy method and a refined decay estimates of products composed of the entropy function and high-order derivatives.

Traveling waves of activity in neural tissue can be halted by spatial inhomogeneity in synaptic coupling. We study an integrate-and-fire network in which each neuron fires once and the coupling decays exponentially. For this model the leading-edge firing map reduces exactly to a scalar equation for the wave speed in space, with a slow unstable and a fast stable homogeneous speed, c_1 and c_2, and a periodic modulation of the coupling enters this equation pointwise. Positive periodic waves terminate in folds. For slowly varying modulation the fold amplitude approaches a plateau set by the trough of the modulation, where the local bottleneck speed is c_1 c_2; the approach to the plateau is set by the local geometry of the trough. For rapidly varying modulation the fold amplitude grows linearly with frequency, with a slope given by a full-amplitude average over the integrated modulation profile; a weak-ripple truncation of that average overestimates the slope about 2.5-fold at the default coupling. A perturbative recursion in the amplitude gives the speed profile explicitly. We distinguish loss of the periodic wave from failure of a particular launch, and the excitatory regime εle 1 from its sign-changing extension. The results are checked against direct integration of the reduced equation and first-spike firing-map simulations of the network.

I examine some analytical properties of a nonlinear reaction-diffusion system that has been used to model the propagation of a wildfire. I establish global-in-time existence and uniqueness of bounded mild solutions to the Cauchy problem for this system given bounded initial data. In particular, this shows that the model does not allow for thermal blow-up. If the initial temperature and fuel density also satisfy certain integrability conditions, the L^2-norms of these global solutions are uniformly bounded in time. Additionally, I use a bootstrap argument to show that small initial temperatures give rise to solutions that decay to zero as time goes to infinity, proving the existence of initial states that do not develop into travelling combustion waves.

We study Hamiltonian elliptic systems with prescribed bilinear mass int_{R^N}uv=a>0. We establish a Gagliardo--Nirenberg inequality in the crossed gradient pairing and bilinear mass, on scaling-invariant component-bounded classes. It identifies the critical curve 1/p+1/q=N/(N+2) and the restricted energy trichotomy. For power nonlinearities, we establish uniqueness up to common translations and nondegeneracy of positive profiles throughout the Sobolev-subcritical hyperbola. Exact scaling then classifies positive normalised solutions and identifies the unique critical mass. In dimension two, we develop a bilinear exponential Gagliardo--Nirenberg estimate with a sharp gradient threshold and a truncated scalar refinement with the optimal quartic leading coefficient. Adapting variational ideas from [Cassani-Tarsi, Calc.Var.PDE (2015)] we combine a mass-normalising quotient, reduction over the full negative fibres and compactness below the concentration threshold to obtain least-energy positive radial solutions for pure exponential nonlinearities at every prescribed mass below the cubic limiting mass. Local perturbative branches and a small-frequency exponential branch yield explicit mass-response formulas.

We study a version of adversarial classification where an adversary is empowered to corrupt data inputs up to some distance varepsilon, using tools from variational analysis. In particular, we describe necessary conditions associated with the optimal classifier subject to such an adversary. Using the necessary conditions, we derive a geometric evolution equation which can be used to track the change in classification boundaries as varepsilon varies. This evolution equation may be described as an uncoupled system of differential equations in one dimension, or as a mean curvature type equation in higher dimension. In one dimension, and under mild assumptions on the data distribution, we rigorously prove that one can use the initial value problem starting from varepsilon=0, which is simply the Bayes classifier, in order to solve for the global minimizer of the adversarial problem for small values of varepsilon. In higher dimensions we provide a similar result, albeit conditional to the existence of regular solutions of the initial value problem. In the process of proving our main results we obtain a result of independent interest connecting the original adversarial problem with an optimal transport problem under no assumptions on whether classes are balanced or not. Numerical examples illustrating these ideas are also presented.

The Redfield equation with time-dependent coefficients (or TCL2 equation) gives one of the simplest perturbative descriptions of open quantum system dynamics derived from a microscopic theory. Its status as an approximation is nevertheless subtle: while perturbation theory fixes the reduced dynamics order-by-order in the system-environment coupling, solving the Redfield equation yields a particular resummation that need not preserve complete positivity. Here we compare Redfield evolution with an exponentiated (Magnus) resummation of the same second-order perturbative map, which is completely positive and trace preserving (CPTP) for a broad class of physical conditions. We also construct a simple CPTP resolvent resummation, illustrating the non-uniqueness of resummations of the same perturbative map. The Redfield and Magnus prescriptions agree to second order but differ through time-ordering corrections at higher orders. We exhibit a weak-coupling criterion which guarantees that Redfield evolution, even when non-Markovian, remains close to a CPTP map. We illustrate the criterion in two examples, focusing on finite-dimensional (qubit) systems. Our results clarify when, and why, Redfield evolution can remain reliable even when it is highly non-Lindbladian.

We investigate in more detail the holographic model of a superconductor recently found by Hartnoll, Herzog, and Horowitz [Phys. Rev. Lett. 101, 031601], which is constructed from a condensate of a charged scalar field in AdS_4-Schwarzschild background. By analytically studying the perturbation of the gravitational system near the critical temperature T_c, we obtain the superconducting coherence length proportional to 1/1-T/T_c via AdS/CFT correspondence. By adding a small external homogeneous magnetic field to the system, we find that a stationary diamagnetic current proportional to the square of the order parameter is induced by the magnetic field. These results agree with Ginzburg-Landau theory and strongly support the idea that a superconductor can be described by a charged scalar field on a black hole via AdS/CFT duality.

For the case of nonlinear second-order differential equations with a constant coefficient of the first derivative term and polynomial nonlinearities, the factorization conditions of Rosu and Cornejo-Perez are approached in two ways: (i) by commuting the subindices of the factorization functions in the two factorization conditions and (ii) by leaving invariant only the first factorization condition achieved by using monomials or polynomial sequences. For the first case, the factorization brackets commute and the generated equations are only equations of Ermakov-Pinney type. The second modification is non commuting, leading to nonlinear equations with different nonlinear force terms, but the same first-order part as the initially factored equation. It is illustrated for monomials with the examples of the generalized Fisher and FitzHugh-Nagumo initial equations. A polynomial sequence example is also included.

We prove strong ill-posedness for the two-dimensional incompressible Euler equation in the endpoint velocity space W^{3,1}_σ(mathbb R^2). Given any compactly supported smooth divergence-free background and any δ>0, we construct a single compactly supported divergence-free datum within δ of the background in W^{3,1}. The perturbation can be localized inside any prescribed nonempty open neighborhood of the background support while remaining disjoint from that support. The associated unique global finite-energy solution u satisfies {rm curl} uin C([0,infty);W^{2,1}(mathbb R^2)), but \[ \operatorname*{ess\,sup}_{0<t<T}\|u(t)\|_{W^{3,1}}=\infty \qquadfor every T>0. \] The proof combines localized hyperbolic insertion, the Bonami--Poornima nonmultiplier mechanism, and many-packet dilution. Moment cancellation and positive material separation control the nonlocal Biot--Savart interactions, allowing a nested construction inside one fixed compact region. Without the compact-support requirement on the limiting datum, a complementary far-field argument yields a dense G_δ set of data whose solutions leave W^{3,1} at every positive rational time and on a residual set of times, with unbounded essential supremum on every nonempty open time interval.

Modern deep neural networks treat parameters as static floating-point matrices stored in physical memory, incurring Von Neumann memory bottlenecks and representation collapse. We formulate Orbital Error Dynamics (OED), an analytical framework wherein synaptic weights are not stored masses (O(W)), but transient topological resonances (O(1)) derived procedurally from the complex quadratic polynomial map z_{n+1} = z_n^2 + c. We introduce the Bent Sine Wave Hypothesis, demonstrating that non-equilibrium living systems emerge when harmonic waves curl inward through environmental drag toward the cardioid cusp (c = 1/4). We define the Observer Horizon Geometry in parameter space, identifying interior resonance shoulder loci X_upper = (0.25, +0.18) and X_lower = (0.25, -0.18) between the fixed-point basin and the true boundary at c = 0.25 +/- 0.50i. To escape non-convex stagnation without loss zeroing, we introduce a heavy-tailed Biomimetic Perturbed Jump Operator (Omega_tunneling) inspired by mammalian fertilization zinc sparks. We further couple an enteric-cranial Dual-Brain architecture shielded by adaptive CD4+ regulatory immune gating (M_CD4), and project the 4-nucleotide genetic basis (A, T, C, G) across quadrants in C. Multi-seed empirical validation on the Two-Moons manifold (5 seeds, 80/20 train/test split, 32x32 grid, zero test-time updates, zero label leakage) demonstrates that procedural parameterization from a 24-byte coordinate seed achieves 77.67% +/- 5.35% clean test accuracy (within an 8.00-point paired difference of an unconstrained gradient baseline at 85.67% +/- 5.35%, 95% CI: [-1.07%, 17.07%]) and 71.33% +/- 3.80% under distribution shift (N(1.2, 0.4)), alongside conceptual equivalence with an analog optical co-processor.

Deep Learning using the eponymous deep neural networks (DNNs) has become an attractive approach towards various data-based problems of theoretical physics in the past decade. There has been a clear trend to deeper architectures containing increasingly more powerful and involved layers. Contrarily, Taylor coefficients of DNNs still appear mainly in the light of interpretability studies, where they are computed at most to first order. However, especially in theoretical physics numerous problems benefit from accessing higher orders, as well. This gap motivates a general formulation of neural network (NN) Taylor expansions. Restricting our analysis to multilayer perceptrons (MLPs) and introducing quantities we refer to as propagators and vertices, both depending on the MLP's weights and biases, we establish a graph-theoretical approach. Similarly to Feynman rules in quantum field theories, we can systematically assign diagrams containing propagators and vertices to the corresponding partial derivative. Examining this approach for S-wave scattering lengths of shallow potentials, we observe NNs to adapt their derivatives mainly to the leading order of the target function's Taylor expansion. To circumvent this problem, we propose an iterative NN perturbation theory. During each iteration we eliminate the leading order, such that the next-to-leading order can be faithfully learned during the subsequent iteration. After performing two iterations, we find that the first- and second-order Born terms are correctly adapted during the respective iterations. Finally, we combine both results to find a proxy that acts as a machine-learned second-order Born approximation.

It is known that barotropic FRW equations written in the conformal time variable can be reduced to simple linear equations for an exponential function involving the conformal Hubble rate. Here, we show that an interesting class of barotropic universes can be obtained in the linear limit of a special type of nonlinear dissipative Ermakov-Pinney equations with the nonlinear dissipation built from Chiellini's integrability condition. These cosmologies, which evolutionary are similar to the standard ones, correspond to barotropic fluids with adiabatic indices rescaled by a particular factor and have amplitudes of the scale factors inverse proportional to the adiabatic index

One way to model the strange metal phase of certain materials is via a holographic description in terms of probe D-branes in a Lifshitz spacetime, characterised by a dynamical exponent z. The background geometry is dual to a strongly-interacting quantum critical theory while the probe D-branes are dual to a finite density of charge carriers that can exhibit the characteristic properties of strange metals. We compute holographically the low-frequency and low-momentum form of the charge density and current retarded Green's functions in these systems for massless charge carriers. The results reveal a quasi-particle excitation when z<2, which in analogy with Landau Fermi liquids we call zero sound. The real part of the dispersion relation depends on momentum k linearly, while the imaginary part goes as k^2/z. When z is greater than or equal to 2 the zero sound is not a well-defined quasi-particle. We also compute the frequency-dependent conductivity in arbitrary spacetime dimensions. Using that as a measure of the charge current spectral function, we find that the zero sound appears only when the spectral function consists of a single delta function at zero frequency.

We investigate Einstein-Gauss-Bonnet (EGB) 4D massive gravity coupled to nonlinear electrodynamics (NED) in an Anti-de-Sitter (AdS) background and find an exact magnetically charged black hole solution. The metric function was analyzed for different values of massive gravity parameters. We verified the first law of black hole thermodynamics and the generalized Smarr formula, treating the cosmological constant as thermodynamic pressure. Vacuum polarization is defined as the conjugate to the NED parameter. To analyze the local stability of the black hole, we compute specific heat. We investigated the Van der Waals-like/reentrant phase transition of the black holes and estimated the critical points. We observed small/large black hole (SBH/LBH) and small/intermediate/large black hole (SBH/IBH/LBH) phase transitions. The Joule-Thomson coefficient, inversion, and isenthalpic curves are discussed. Finally, the minimum inversion temperature and the corresponding event horizon radius are obtained using numerical techniques.

We consider a Schr\"odinger-Poisson system involving a general nonlinearity at critical growth and we prove the existence of positive solutions. The Ambrosetti-Rabinowitz condition is not required. We also study the asymptotics of solutions with respect to a parameter.

We establish a conditional theory of genuinely non-perturbative global dynamics for the cutoff Boltzmann equation with hard potentials on the three-dimensional torus. We prove that any global mild solution satisfying \[ \sup_{t\ge0}\int_{\mathbb R^3}W(v)\sup_{x\in\mathbb T^3}F(t,x,v)\,dv<\infty, \] for a suitable stretched-exponential weight \(W\), enjoys uniform weighted \(L^\infty\) control, is unique in the dissipative mild class, and converges exponentially to the global Maxwellian. In particular, these conclusions hold for arbitrary-amplitude solutions admitting a uniform stretched-exponential velocity upper bound. The key mechanism is a *tail-to-coercivity principle*: uniform velocity-tail control generates positive lower bounds for the local density and nonlinear collision frequency, without imposing macroscopic coercivity a priori. A double-Duhamel positivity argument also generates a Gaussian lower bound at arbitrarily short positive times, allowing initial data with local vacuum. We establish a polynomially weighted counterpart under an additional spatially uniform lower velocity profile, obtaining analogous weighted \(L^\infty\) bounds, uniqueness, and exponential relaxation. Our results provide a low-regularity alternative to conditional theories requiring smoothness and pointwise macroscopic bounds. They reveal how a single velocity-tail condition controls large-data Boltzmann dynamics and suggest a route toward problems with physical boundary conditions.

Singular regular points often arise in differential equations describing physical phenomena such as fluid dynamics, electromagnetism, and gravitation. Traditional numerical techniques often fail or become unstable near these points, requiring the use of semi-analytical tools, such as series expansions and perturbative methods, in combination with numerical algorithms; or to invoke more sophisticated methods. In this work, we take an alternative route and leverage the power of machine learning to exploit Physics Informed Neural Networks (PINNs) as a modern approach to solving ordinary differential equations with singular points. PINNs utilize deep learning architectures to approximate solutions by embedding the differential equations into the loss function of the neural network. We discuss the advantages of PINNs in handling singularities, particularly their ability to bypass traditional grid-based methods and provide smooth approximations across irregular regions. Techniques for enhancing the accuracy of PINNs near singular points, such as adaptive loss weighting, are used in order to achieve high efficiency in the training of the network. We exemplify our results by studying four differential equations of interest in mathematics and gravitation -- the Legendre equation, the hypergeometric equation, the solution for black hole space-times in theories of Lorentz violating gravity, and the spherical accretion of a perfect fluid in a Schwarzschild geometry.

We discuss the numerical solution of initial value problems for varepsilon^2,varphi''+a(x),varphi=0 in the highly oscillatory regime, i.e., with a(x)>0 and 0<varepsilonll 1. We analyze and implement an approximate solution based on the well-known WKB-ansatz. The resulting approximation error is of magnitude O(varepsilon^{N}) where N refers to the truncation order of the underlying asymptotic series. When the optimal truncation order N_{opt} is chosen, the error behaves like O(varepsilon^{-2}exp(-cvarepsilon^{-1})) with some c>0.

It is well established that the spectral analysis of canonically quantized four-dimensional Seiberg-Witten curves can be systematically studied via the Nekrasov-Shatashvili functions. In this paper, we explore another aspect of the relation between N=2 supersymmetric gauge theories in four dimensions and operator theory. Specifically, we study an example of an integral operator associated with Painlev\'e equations and whose spectral traces are related to correlation functions of the 2d Ising model. This operator does not correspond to a canonically quantized Seiberg-Witten curve, but its kernel can nevertheless be interpreted as the density matrix of an ideal Fermi gas. Adopting the approach of Tracy and Widom, we provide an explicit expression for its eigenfunctions via an O(2) matrix model. We then show that these eigenfunctions are computed by surface defects in SU(2) super Yang-Mills in the self-dual phase of the Omega-background. Our result also yields a strong coupling expression for such defects which resums the instanton expansion. Even though we focus on one concrete example, we expect these results to hold for a larger class of operators arising in the context of isomonodromic deformation equations.

We calculate the Rayleigh-Schrodinger perturbation series for the ground-state energy and the effective mass of the three-dimensional Frohlich polaron through fourth order in the coupling constant alpha. The third- and fourth-order coefficients of the effective mass are new. Because the Frohlich interaction is a form-bounded perturbation of an isolated nondegenerate ground state at fixed total momentum, the series has a nonzero radius of convergence, and we analyze it with the ratio and Pade methods appropriate to convergent series. The coefficients are strikingly regular and are described by a simple power-law singularity at alpha near 8.5 with an exponent near 1.4. The couplings at which successive truncations of the inverse mass vanish, the first of which is the textbook breakdown value alpha = 6, are slowly converging estimates of this radius. Pade approximants built from the weak-coupling coefficients alone agree with diagrammatic Monte Carlo results for the mass to better than one percent for alpha up to 5. Exact results on the analyticity of the ground state and on the strong-coupling asymptotics imply that the coefficients cannot keep a constant sign, although all known coefficients do, so that the apparent singularity must be a complex-conjugate pair close to the positive real axis. It marks the crossover to strong coupling: the weak-coupling series locates this crossover accurately but cannot be continued through it.

We study complex saddle points of the worldline path integral for massive charged particles in two-dimensional maximally symmetric backgrounds threaded by homogeneous electric fields. Within the worldline formalism, we show that the semiclassical saddles generically require complexifying both the target spacetime and the proper-time parameter. The resulting solutions comprise two slices of a single complex family, distinguished by the direction of proper-time evolution. The one evolving in imaginary proper time carries an action that typically controls the exponential suppression of non-perturbative effects, while that evolving in real proper time may encode additional complex phases responsible for quantum interference between saddles. Treating this family as a whole, we give a prescription for constructing and classifying all topologically inequivalent instanton sectors, and we find that the worldline action depends only on their topological data, namely winding numbers around the compact directions and certain isolated punctures. Moreover, in Anti-de Sitter spacetime, the qualitative structure of the solutions is governed by the dimensionless ratio R |eE|/m (R being the curvature radius), which distinguishes sub-extremal, extremal, and super-extremal regimes. By contrast, no analogous distinction occurs in de Sitter---where the instantons are insensitive to this ratio---or in Minkowski spacetime, which lacks a curvature scale altogether. We also recast the action as a sum over residues of the relevant (non-conserved) conjugate momentum, making transparent which analytic data control each contribution. Finally, we discuss the subtleties that arise when solutions built in various local patches are lifted to global (Anti-)de Sitter space. These impose further necessary conditions on which saddles can in principle contribute to the one-loop path integral.

Ap\'{e}ry-type (inverse) binomial series have appeared prominently in the calculations of Feynman integrals in recent years. In our previous work, we showed that a few large classes of the non-alternating Ap\'ery-type (inverse) central binomial series can be evaluated using colored multiple zeta values of level four (i.e., special values of multiple polylogarithms at fourth roots of unity) by expressing them in terms of iterated integrals. In this sequel, we shall prove that for several classes of the alternating versions we need to raise the level to eight. Our main idea is to adopt hyperbolic trigonometric 1-forms to replace the ordinary trigonometric ones used in the non-alternating setting.

In this paper we consider a particular class of solutions of the Rayleigh-Boltzmann equation, known in the nonlinear setting as homoenergetic solutions, which have the form gleft( x,v,t right) =fleft( v-Lleft( tright)x,tright) where the matrix L(t) describes a shear flow deformation. We began this analysis in [22] where we rigorously proved the existence of a stationary non-equilibrium solution and established the different behaviour of the solutions for small and large values of the shear parameter, for cut-off collision kernels with homogeneity parameter 0leq gamma <1, including Maxwell molecules and hard potentials. In this paper, we concentrate in the case where the deformation term dominates the collision term for large times (hyperbolic-dominated regime). This occurs for collision kernels with gamma < 0 and in particular we focus on gamma in (-1,0). In such a hyperbolic-dominated regime, it appears challenging to provide a clear description of the long-term asymptotics of the solutions. Here we present a formal analysis of the long-time asymptotics for the distribution of velocities and provide the explicit form for the asymptotic profile. Additionally, we discuss the different asymptotic behaviour expected in the case of homogeneity gamma < -1. Furthermore, we provide a probabilistic interpretation describing a stochastic process consisting in a combination of collisions and shear flows. The tagged particle velocity {v(t)}_{tgeq 0} is a Markov process that arises from the combination of free flights in a shear flow along with random jumps caused by collisions.

The dynamics of Einstein-conformally coupled Higgs field (EccH) system is investigated near the initial singularities in the presence of Friedman-Robertson--Walker symmetries. We solve the field equations asymptotically up to fourth order near the singularities analytically, and determine the solutions numerically as well. We found all the asymptotic, power series singular solutions, which are (1) solutions with a scalar polynomial curvature singularity but the Higgs field is bounded (`Small Bang'), or (2) solutions with a Milne type singularity with bounded spacetime curvature and Higgs field, or (3) solutions with a scalar polynomial curvature singularity and diverging Higgs field (`Big Bang'). Thus, in the present EccH model there is a new kind of physical spacetime singularity (`Small Bang'). We also show that, in a neighbourhood of the singularity in these solutions, the Higgs sector does not have any symmetry breaking instantaneous vacuum state, and hence then the Brout-Englert-Higgs mechanism does not work. The large scale behaviour of the solutions is investigated numerically as well. In particular, the numerical calculations indicate that there are singular solutions that cannot be approximated by power series.

In this work, we prove existence and uniqueness of a bounded viscosity solution for the Cauchy problem of degenerate parabolic equations with variable exponent coefficients. We construct the solution directly using the stochastic representation, then verify it satisfies the Cauchy problem. The corresponding SDE, on the other hand, allows the drift and diffusion coefficients to respond nonlinearly to the current state through the state-dependent variable exponents, and thus, extends the expressive power of classical SDEs to better capture complex dynamics. To validate our theoretical framework, we conduct comprehensive numerical experiments comparing finite difference solutions (Crank-Nicolson on logarithmic grids) with Monte Carlo simulations of the SDE.

We propose a novel neural algorithm for the fundamental problem of computing the entropic optimal transport (EOT) plan between continuous probability distributions which are accessible by samples. Our algorithm is based on the saddle point reformulation of the dynamic version of EOT which is known as the Schr\"odinger Bridge problem. In contrast to the prior methods for large-scale EOT, our algorithm is end-to-end and consists of a single learning step, has fast inference procedure, and allows handling small values of the entropy regularization coefficient which is of particular importance in some applied problems. Empirically, we show the performance of the method on several large-scale EOT tasks. https://github.com/ngushchin/EntropicNeuralOptimalTransport

We implement physics-informed neural networks (PINNs) to solve the time-independent Schr\"odinger equation for three canonical one-dimensional quantum potentials: an infinite square well, a finite square well, and a finite barrier. The PINN models incorporate trial wavefunctions that exactly satisfy boundary conditions (Dirichlet zeros at domain boundaries), and they optimize a loss functional combining the PDE residual with a normalization constraint. For the infinite well, the ground-state energy is known (E = pi^2 in dimensionless units) and held fixed in training, whereas for the finite well and barrier, the eigenenergy is treated as a trainable parameter. We use fully-connected neural networks with smooth activation functions to represent the wavefunction and demonstrate that PINNs can learn the ground-state eigenfunctions and eigenvalues for these quantum systems. The results show that the PINN-predicted wavefunctions closely match analytical solutions or expected behaviors, and the learned eigenenergies converge to known values. We present training logs and convergence of the energy parameter, as well as figures comparing the PINN solutions to exact results. The discussion addresses the performance of PINNs relative to traditional numerical methods, highlighting challenges such as convergence to the correct eigenvalue, sensitivity to initialization, and the difficulty of modeling discontinuous potentials. We also discuss the importance of the normalization term to resolve the scaling ambiguity of the wavefunction. Finally, we conclude that PINNs are a viable approach for quantum eigenvalue problems, and we outline future directions including extensions to higher-dimensional and time-dependent Schr\"odinger equations.

In the present work, we study the dynamical evolution of an homogeneous and anisotropic, noncommutative (NC) Bianchi I (BI) model coupled to a radiation perfect fluid. Our first motivation is determining if the present model tends to an homogeneous and isotropic NC Friedmann-Robertson-Walker (FRW) model, during its evolution. In order to simplify our task, we use the Misner parametrization of the BI metric. In terms of that parametrization the BI metric has three metric functions: the scale factor a(t) and the two parameters beta_pm (t), which measure the spatial anisotropy of the model. Our second motivation is trying to describe the present accelerated expansion of the universe using noncommutativity (NCTY). The NCTY is introduced by two nontrivial Poisson brackets between some geometrical as well as matter variables of the model. We recover the description in terms of commutative variables by introducing some variables transformations that depend on the NC parameter. Using those variables transformations, we rewrite the total NC Hamiltonian of the model in terms of commutative variables. From the resulting Hamiltonian, we obtain the dynamical equations for a generic perfect fluid. In order to solve these equations, we restrict our attention to a model where the perfect fluid is radiation. We solve, numerically, these equations and compare the NC solutions to the corresponding commutative ones. The comparison shows that the NC model may be considered as a possible candidate for describing the accelerated expansion of the universe. Finally, we obtain estimates for the NC parameter and compare the main results of the NC BI model coupled to radiation with the same NC BI model coupled to other perfect fluids. As our main result, we show that the solutions, after some time, produce an isotropic universe.

To model high-risk traffic densities a cellular automaton model is constructed, exhibiting Ising model-like properties. The attraction-repulsion forces between vehicles are evaluated as coupling, and the coupling function p(1-p)=g/8, where p is the initial distribution of cells with state 1, and g is the Moore neighbor count, which has roots cos^2(π/8) and sin^2(π/8) for g=1. This is achieved without the use of any trigonometric functions in the code. From the roots, the tangent polynomial tan^2(x) + tan(x) emerges. Kramers-Wannier duality is recovered and it is conjectured that the 1/2 difference between the roots serves as a fixed point for a projection mechanism from the 2-dimensional Ising model onto a 1-dimensional Ising chain through the Gudermannian function. The sigmoid evaluated at the proposed fixed point is substituted into the derivative of the logistic coupling function, σ(1/2)(1-σ(1/2)), yielding a numerical approximation to the three-dimensional Ising inverse critical coupling. Finally, the results are linked to risk densities in traffic and vehicle types, accounting for the amplification of fatal accidents.

In the context of supersymmetric quantum mechanics, we define a potential through a particular Riccati solution of the composition form, F(f(x)), and obtain a generalized Mielnik construction of one-parameter isospectral potentials when we use the general Riccati solution. Some examples for special cases of F and f are given to illustrate the method. An interesting result is obtained in the case of a parametric double well potential generated by this method, for which it is shown that the parameter of the potential controls the heights of the localization probability in the two wells, and for certain values of the parameter the height of the localization probability can be higher in the smaller well

We study the one-dimensional nonlocal Kirchhoff type bifurcation problem related to logistic equation of population dynamics. We establish the precise asymptotic formulas for bifurcation curve lambda = lambda(alpha) as alpha to infty in L^2-framework, where alpha:= Vert u_lambda Vert_2.

According to the generally accepted phase diagram of QCD, at low temperature and high baryon number density the chiral phase transition of QCD is of first order and the co-existence of the Nambu-Goldstone phase and the Wigner phase should appear. This is in conflict with the usual claim that the quark gap equation has no Wigner solution in the case of nonzero current quark mass. In this paper we analyze the reason why the Wigner solution does not exist in the usual treatment and try to propose a new approach to discuss this question. As a first step, we adopt a modified Nambu-Jona-Lasinio (NJL) model to study the Wigner solution at finite current quark mass. We then generalize this approach to the case of finite chemical potential and discuss partial restoration of chiral symmetry at finite chemical potential and compare our results with those in the normal NJL model.

Traveling wavetrains in generalized two-species predator-prey models and two-component reaction-diffusion equations are considered. The stability of the fixed points of the traveling wave ODEs (in the usual "spatial" variable) is considered. For general functional forms of the nonlinear prey birthrate/prey deathrate or reaction terms, a Hopf bifurcation is shown to occur at two different critical values of the traveling wave speed. The post-bifurcation dynamics is investigated for five different functional forms of the nonlinearities. In cases where the bifurcation is supercritical, the post-bifurcation behaviour yields stable periodic orbits of the traveling-wave ODEs in the spatial variable. These correspond to stable periodic wavetrains of the full PDEs. Subcritical Hopf bifurcations yield more complex post-bifurcation dynamics in the PDE wavetrains. In special cases where the subcritical bifurcation marks the end of the regime of stability, the post-bifurcation behavior in the spatial ODEs is chaotic, corresponding to wavetrains of the original PDEs which are spatially coherent, but have chaotic temporal dynamics. All the models are integrated numerically to investigate the post-bifurcation dynamics and chaotic regimes are characterized by computing power spectra, autocorrelation functions, and fractal dimensions.

We investigate the effects of quantum geometry on the motion of magnon wave packets in one-dimensional quantum spin chains. In contrast to fermionic systems, where transport phenomena driven by the quantum metric have been extensively studied, we show that the dynamics in bosonic magnon systems are strongly affected by the symplectic structure and the pseudo-unitarity of the transformation of the magnon Bogoliubov-de Gennes Hamiltonian. A chiral soliton lattice provides a setting in which the symplectic quantum metric strongly influences the longitudinal magnon dynamics. We derive the semiclassical response of a magnon wave packet in such a system to a weak, slowly varying perturbation and demonstrate that the wave packet velocity acquires a geometric contribution proportional to the momentum derivative of the symplectic quantum metric, which dominates over the conventional group velocity for a nearly flat band. Solving the resulting equations of motion, we obtain a bounded trajectory in real time with a phase-space area determined by the symplectic quantum metric.

We characterise the dynamics of neuronal activity, in terms of field theory, using neural units placed on a 2D-lattice modelling the cortical surface. The electrical activity of neuronal units was analysed with the aim of deriving a neural field model with a simple functional form that still able to predict or reproduce empirical findings. Each neural unit was modelled using a neural mass and the accompanying field theory was derived in the continuum limit. The field theory comprised coupled (real) Klein-Gordon fields, where predictions of the model fall within the range of experimental findings. These predictions included the frequency spectrum of electric activity measured from the cortex, which was derived using an equipartition of energy over eigenfunctions of the neural fields. Moreover, the neural field model was invariant, within a set of parameters, to the dynamical system used to model each neuronal mass. Specifically, topologically equivalent dynamical systems resulted in the same neural field model when connected in a lattice; indicating that the fields derived could be read as a canonical cortical field theory. We specifically investigated non-dispersive fields that provide a structure for the coding (or representation) of afferent information. Further elaboration of the ensuing neural field theory, including the effect of dispersive forces, could be of importance in the understanding of the cortical processing of information.

We provide a rigorous derivation of Boltzmann's kinetic equation from the hard sphere system for rarefied gas, which is valid for arbitrarily long times, as long as the solution to the Boltzmann equation exists. This extends Lanford's landmark theorem (1975), which justifies this derivation for a sufficiently short time. In a companion paper (arXiv:2503.01800), we connect this derivation to existing literature on hydrodynamic limits. This completes the resolution of Hilbert's Sixth Problem pertaining to the derivation of fluid equations from Newton's laws, in the case of a rarefied, hard sphere gas. The general strategy follows the paradigm introduced by the first two authors for the long-time derivation of the wave kinetic equation in wave turbulence theory. This is based on propagating a long-time cumulant ansatz, which keeps memory of the full collision history of the relevant particles, by an important partial time expansion. The heart of the matter is proving the smallness of these cumulants in L^1, which can be reduced to combinatorial properties for the associated diagrams which we call molecules. These properties are then proved by devising an elaborate cutting algorithm, which is a major novelty of this work.

We study two-point functions of chiral primary operators in N=4 Super Yang-Mills in the presence of a half-BPS surface defect. Such a defect preserves only the diagonal combination of transverse rotations and of a global U(1) R-symmetry, so that the correlator depends on one additional cross-ratio and splits into sectors of definite angular charge. We solve the superconformal Ward identities and construct the bulk and defect superconformal blocks in all sectors, including the charged ones, which have no analogue for defects preserving transverse rotations. The solution contains a holomorphic sector, which we fix non-perturbatively by mapping it to a chiral algebra correlator. We extend this map to defects breaking transverse rotations, where it requires a redefinition of the modes of the two-dimensional current algebra, and use it to fix the charged sectors as well. We then use the defect dispersion relation to reconstruct the full correlator of the dimension-two chiral primary at leading order at large N and strong coupling. In the charged sectors long multiplets are absent and the answer is conjecturally exact in both N and the coupling; expanding it at weak coupling, we find perfect agreement with known results. From these correlators we extract infinitely many new pieces of defect CFT data. Finally, we reproduce the strong-coupling answer by an independent Witten diagram bootstrap on the rotation-invariant locus.

Fourier Neural Operators (FNO) learn solution operators of partial differential equations by parameterizing global convolutions in the complex Fourier domain. For real-valued PDE solutions, the complex FFT carries representational redundancy through conjugate symmetry. We introduce the Hartley Neural Operator (HNO), the exact real-valued mirror of FNO: it replaces the FFT with the purely real Discrete Hartley Transform and learns a single real multiplier per retained spectral mode, with no complex arithmetic. Because the real Hartley spectrum is not halved by conjugate symmetry, HNO retains twice as many frequency corners as FNO but one real weight where FNO carries a complex pair, so the two operators are iso-parametric at equal width and differ only in spectral basis. Our central thesis is that the best basis is a property of the operator. Self-adjoint elliptic operators (Poisson, biharmonic) have real, symmetric Green's functions that the real Hartley multiplier diagonalizes exactly, and HNO is favored there. Time-dependent operators carry phase, from oscillation in the wave equation to transport in advection, Burgers, and Navier-Stokes, which a real diagonal multiplier cannot represent, so FNO is favored there, and increasingly so with the operator's phase content, leaving the phaseless heat equation as the borderline case. Training both operators identically and benchmarking across PDE classes, initial-condition families, and boundary conditions, we find an elliptic-versus-time-dependent split that is monotone in operator phase content and matches the Green's-function theory we develop. Rather than a universal winner, our findings give a predictive rule: match the spectral basis to the symmetry of the solution operator.

By introducing a dynamical version of the second fundamental form, we generalize a recent result of Filip-Fisher-Lowe to the setting of magnetic systems. Namely, we show that a real-analytic negatively s-curved magnetic system on a closed real-analytic manifold has only finitely many closed totally s-magnetic hypersurfaces, unless the magnetic 2-form is trivial and the underlying metric is hyperbolic.

We address the challenge of sound propagation simulations in 3D virtual rooms with moving sources, which have applications in virtual/augmented reality, game audio, and spatial computing. Solutions to the wave equation can describe wave phenomena such as diffraction and interference. However, simulating them using conventional numerical discretization methods with hundreds of source and receiver positions is intractable, making stimulating a sound field with moving sources impractical. To overcome this limitation, we propose using deep operator networks to approximate linear wave-equation operators. This enables the rapid prediction of sound propagation in realistic 3D acoustic scenes with moving sources, achieving millisecond-scale computations. By learning a compact surrogate model, we avoid the offline calculation and storage of impulse responses for all relevant source/listener pairs. Our experiments, including various complex scene geometries, show good agreement with reference solutions, with root mean squared errors ranging from 0.02 Pa to 0.10 Pa. Notably, our method signifies a paradigm shift as no prior machine learning approach has achieved precise predictions of complete wave fields within realistic domains. We anticipate that our findings will drive further exploration of deep neural operator methods, advancing research in immersive user experiences within virtual environments.

We establish lower bounds on the number of resonances for three-dimensional semiclassical Dirac operators with mass-type potentials. For a pair Dl_j=Dl_0+mbfβv_j, j=1,2, satisfying the exterior analyticity, decay and non-degeneracy assumptions of the paper, we associate a relative distribution μ to the corresponding mass functions. We prove that if E_0>1 is a boundary point of suppμ, then every complex neighbourhood of E_0, and likewise of -E_0, contains in total at least Chbar^{-3} resonances of the two operators, counted with multiplicity. The proof transfers an analytic wave-front singularity of μ to the relative classical density of states and combines this with Khochman's local trace formula and a relative semiclassical trace formula. As a consequence, if v is a single bounded mass-type potential, analytic with power decay in an exterior sector, with inf(1+v)>0 and sup v>0, then Dl=Dl_0+mbfβv has at least Chbar^{-3} resonances in every complex neighbourhood of E_0=1+sup v and of -E_0, for all sufficiently small hbar. The exponent agrees with Khochman's upper bound, the two resonance families are related by spectral symmetry, and E_0 is necessarily a trapping energy. We also prove the relative semiclassical trace formula needed in the argument for general Hermitian matrix-valued perturbations.