The problem of non-integrability of the circular restricted three-body problem is very classical and important in the theory of dynamical systems. It was partially solved by Poincaré in the nineteenth century: he showed that there exists no real-analytic first integral which depends analytically on the mass ratio of the second body to the total and is functionally independent of the Hamiltonian. When the mass of the second body becomes zero, the restricted three-body problem reduces to the two-body Kepler problem. We prove the non-integrability of the restricted three-body problem both in the planar and spatial cases for any non-zero mass of the second body. Our basic tool of the proofs is a technique developed here for determining whether perturbations of integrable systems which may be non-Hamiltonian are not meromorphically integrable near resonant periodic orbits such that the first integrals and commutative vector fields also depend meromorphically on the perturbation parameter. The technique is based on generalized versions due to Ayoul and Zung of the Morales–Ramis and Morales–Ramis–Simó theories. We emphasize that our results are not just applications of the theories.

The singularly perturbed Riccati equation is the first-order nonlinear ordinary differential equation $\hbar \partial _x f = af^2 + bf + c$ in the complex domain where $\hbar $ is a small complex parameter. We prove an existence and uniqueness theorem for exact solutions with prescribed asymptotics as $\hbar \to 0$ in a half-plane. These exact solutions are constructed using the Borel–Laplace method; that is, they are Borel summations of the formal divergent $\hbar $ -power series solutions. As an application, we prove existence and uniqueness of exact WKB solutions for the complex one-dimensional Schrödinger equation with a rational potential.

Networked dynamical systems, i.e., systems of dynamical units coupled via nontrivial interaction topologies, constitute models of broad classes of complex systems, ranging from gene regulatory and metabolic circuits in our cells to pandemics spreading across continents. Most of such systems are driven by irregular and distributed fluctuating input signals from the environment. Yet how networked dynamical systems collectively respond to such fluctuations depends on the location and type of driving signal, the interaction topology and several other factors and remains largely unknown to date. As a key example, modern electric power grids are undergoing a rapid and systematic transformation towards more sustainable systems, signified by high penetrations of renewable energy sources. These in turn introduce significant fluctuations in power input and thereby pose immediate challenges to the stable operation of power grid systems. How power grid systems dynamically respond to fluctuating power feed-in as well as other temporal changes is critical for ensuring a reliable operation of power grids yet not well understood. In this work, we systematically introduce a linear response theory (LRT) for fluctuation-driven networked dynamical systems. The derivations presented not only provide approximate analytical descriptions of the dynamical responses of networks, but more importantly, also allow to extract key qualitative features about spatio-temporally distributed response patterns. Specifically, we provide a general formulation of a LRT for perturbed networked dynamical systems, explicate how dynamic network response patterns arise from the solution of the linearised response dynamics, and emphasise the role of LRT in predicting and comprehending power grid responses on different temporal and spatial scales and to various types of disturbances. Understanding such patterns from a general, mathematical perspective enables to estimate network responses quickly and intuitively, and to develop guiding principles for, e.g., power grid operation, control and design.

We consider a spherical particle levitating above a liquid bath owing to the Leidenfrost effect, where the vapour of either the bath or sphere forms an insulating film whose pressure supports the sphere’s weight. Starting from a reduced formulation based on a lubrication-type approximation, we use matched asymptotics to describe the morphology of the vapour film assuming that the sphere is small relative to the capillary length (small Bond number) and that the densities of the bath and sphere are comparable. We find that this regime is comprised of two formally infinite sequences of distinguished limits which meet at an accumulation point, the limits being defined by the smallness of an intrinsic evaporation number relative to the Bond number. These sequences of limits reveal a surprisingly intricate evolution of the film morphology with increasing sphere size. Initially, the vapour film transitions from a featureless morphology, where the thickness profile is parabolic, to a neck–bubble morphology, which consists of a uniform pressure bubble bounded by a narrow and much thinner annular neck. Gravity effects then become important in the bubble leading to sequential formation of increasingly smaller neck–bubble pairs near the symmetry axis. This process terminates when the pairs closest to the symmetry axis become indistinguishable and merge. Subsequently, the inner section of that merger transitions into a uniform-thickness film that expands radially, gradually squishing increasingly larger neck–bubble pairs into a region of localised oscillations sandwiched between the uniform film and what remains of the bubble whose radial extent is presently comparable to the uniform film; the neck–bubble pairs farther from the axis remain essentially intact. Ultimately, the uniform film gobbles up the largest outermost bubble, whereby the morphology simplifies to a uniform film bounded by localised oscillations. Overall, the asymptotic analysis describes the continuous evolution of the vapour film from a neck–bubble morphology typical of a Leidenfrost drop levitating above a flat solid substrate to a uniform-film morphology which resembles that in the case of a large liquid drop levitating above a liquid bath.

We study the extended Stefan problem which includes constitutional supercooling for the solidification of a binary alloy in a finite spherical domain. We perform an asymptotic analysis in the limits of large Lewis number and small Stefan number which allows us to identify a number of spatio-temporal regimes signifying distinct behaviours in the solidification process, resulting in an intricate boundary layer structure. Our results generalise those present in the literature by considering all time regimes for the Stefan problem while also accounting for impurities and constitutional supercooling. These results also generalise recent work on the extended Stefan problem for finite planar domains to spherical domains, and we shall highlight key differences in the asymptotic solutions and the underlying boundary layer structure which result from this change in geometry. We compare our asymptotic solutions with both numerical simulations and real experimental data arising from the casting of molten metallurgical grade silicon through the water granulation process, with our analysis highlighting the role played by supercooling in the solidification of binary alloys appearing in such applications.

Modelling a normal–superconducting interface, we consider a semi-infinite wire whose edge is adjacent to a normal magnetic metal, assuming asymptotic convergence, away from the boundary, to the purely superconducting state. We obtain that the maximal current which can be carried by the interface diminishes in the small normal conductivity limit.

In this paper, we are concerned with unbounded solutions of the singularly perturbed forced pendulum equation in the presence of friction, namely

$$\varepsilon ^2u_\varepsilon ^{{\prime}{\prime}} + \sin u_\varepsilon = \varepsilon ^2\alpha (t)u_\varepsilon + \varepsilon ^2\beta (t)u_\varepsilon ^{\prime} \quad {\rm in}\;(-L,L).{\rm }$$

Applications of a WKBJ-type ‘ray ansatz’ to obtain asymptotic solutions of the Helmholtz equation in the high-frequency limit are now standard and underpin the construction of ‘geometrical optics’ ray diagrams in many electromagnetic, acoustic and elastic reflection, transmission and other scattering problems. These applications were subsequently extended by Keller to include other types of rays – called ‘diffracted’ rays – to provide an accessible and impressively accurate theory which is relevant in wide-ranging sets of circumstances. Friedlander and Keller then introduced a modified ray ansatz to extend yet further the scope of ray theory and its applicability to certain other classes of diffraction problems (tangential ray incidence upon an obstructing boundary, for instance) and did so by the inclusion of an extra term proportional to a power of the wave number within the exponent of the initial ansatz. Our purpose here is to generalise this further still by the inclusion of several such terms, ordered in a natural sequence in terms of strategically chosen fractional powers of the large wave number, and to derive a systematic sequence of boundary value problems for the coefficient phase functions that arise within this generalised exponent, as well as one for the leading-order amplitude occurring as a pre-exponential factor. One particular choice of fractional power is considered in detail, and waves with specified radially symmetric or planar wavefronts are then analysed, along with a boundary value problem typifying two-dimensional radiation whereby arbitrary phase and amplitude variations are specified on a prescribed boundary curve. This theory is then applied to the scattering of plane and cylindrical waves at curved boundaries with small-scale perturbations to their underlying profile.

In this paper, periodic steady-state of a liquid film flowing over a periodic uneven wall is investigated via a classical method. Specifically, we analyze a long-wave model that is valid at the near-critical Reynolds number. For the periodic wall surface, we construct an iteration scheme in terms of an integral form of the original steady-state problem. The uniform convergence of the scheme is proved so that we can derive the existence and the uniqueness as well as the asymptotic formula of the periodic solutions.

Let $\mathbf{S}:={{\{S(t)\}}_{t\ge 0}}$ be a ${{\text{C}}_{0}}$-semigroup of quasinilpotent operators (i.e., $\sigma (S(t))=\{0\}$ for each $t>0$). In dynamical systems theory the above quasinilpotency property is equivalent to a very strong concept of stability for the solutions of autonomous systems. This concept is frequently called superstability and weakens the classical finite time extinction property (roughly speaking, disappearing solutions). We show that under some assumptions, the quasinilpotency, or equivalently, the superstability property of a ${{\text{C}}_{0}}$-semigroup is preserved under the perturbations of its infinitesimal generator.

Entropy maximising spatial interaction models have been widely exploited in a range of disciplines and applications: from trade and migration flows to the spread of riots and the understanding of spatial patterns in archaeological sites of interest. When embedded into a dynamic system and framed in the context of a retail model, the dynamics of centre growth poses an interesting mathematical problem, with bifurcations and phase changes, which may be addressed analytically. In this paper, we present some analysis of the continuous retail model and the corresponding discrete version, which yields insights into the effect of space on the evolving system, and an understanding of why certain retail centres are more successful than others. The slowly developing growths and the fast explosive growths that are of particular concern are explained in detail.

In this paper, we study the stability in the Lyapunov sense of the equilibrium solutions of discrete or difference Hamiltonian systems in the plane. First, we perform a detailed study of linear Hamiltonian systems as a function of the parameters. In particular we analyze the regular and the degenerate cases. Next, we give a detailed study of the normal form associated with the linear Hamiltonian system. At the same time we obtain the conditions under which we can get stability (in linear approximation) of the equilibrium solution, classifying all the possible phase diagrams as a function of the parameters. After that, we study the stability of the equilibrium solutions of the first order difference system in the plane associated with mechanical Hamiltonian systems and Hamiltonian systems defined by cubic polynomials. Finally, we point out important differences with the continuous case.

We present a review on the accuracy of asymptotic models for the scattering problem of electromagnetic waves in domains with thin layer. These models appear as first order approximations of the electromagnetic field. They are obtained thanks to a multiscale expansion of the exact solution with respect to the thickness of the thin layer, that makes possible to replace the thin layer by approximate conditions. We present the advantages and the drawbacks of several approximations together with numerical validations and simulations. The main motivation of this work concerns the computation of electromagnetic field in biological cells. The main difficulty to compute the local electric field lies in the thinness of the membrane and in the high contrast between the electrical conductivities of the cytoplasm and of the membrane, which provides a specific behavior of the electromagnetic field at low frequencies.

A susceptible–exposed–infectious theoretical model describing Tasmanian devil population and disease dynamics is presented and mathematically analysed using a dynamical systems approach to determine its behaviour under a range of scenarios. The steady states of the system are calculated and their stability analysed. Closed forms for the bifurcation points between these steady states are found using the rate of removal of infected individuals as a bifurcation parameter. A small-amplitude Hopf region, in which the populations oscillate in time, is shown to be present and subjected to numerical analysis. The model is then studied in detail in relation to an unfolding parameter which describes the disease latent period. The model’s behaviour is found to be biologically reasonable for Tasmanian devils and potentially applicable to other species.

We improve several recent results in the asymptotic integration theory of nonlinear ordinary differential equations via a variant of the method devised by J. K. Hale and N. Onuchic The results are used for investigating the existence of positive solutions to certain reaction-diffusion equations.

The solutions of a perturbed linear ordinary differential equation (ODE) system are studied. Provided that some integrability and oddness conditions are satisfied, we show that they are asymptotically equivalent at t = ±∞ to the solutions of the unperturbed one. This fact is used to determine the existence of almost periodic or pseudo-almost periodic solutions of the perturbed system.

We analyze asymptotically a differential-difference equation that arises in a Markov-modulated fluid model. We use singular perturbation methods to analyze the problem with appropriate scalings of the two state variables. In particular, the ray method and asymptotic matching are used.

We consider the M/M/∞ queue with m primary servers and infinitely many secondary servers. All the servers are numbered and ordered. An arriving customer takes the lowest available server. We define the wasted spaces as the difference between the highest numbered occupied server and the total number of occupied servers. Letting ρ = λ0/μ be the ratio of arrival to service rates, we study the probability distribution of the wasted spaces asymptotically for ρ → ∞. We also give some numerical results and the tail behavior for ρ = O(1).

We study the semilinear differential equation u″ + F(t,u,u′)=0 on a half-line. Under different growth conditions on the function F, equations with globally defined solutions asymptotic to lines are characterized. Both fixed initial data and fixed asymptote are considered.

For a class of second order nonlinear differential equations, sufficient conditions are presented to ensure that some, respectively all solutions are asymptotic to lines.