Table of contents

This paper is devoted to the numerical computation and continuation of families of heteroclinic connections between hyperbolic periodic orbits (POs) of a Hamiltonian system. We describe a method that requires the numerical continuation of a nonlinear system that involves the initial conditions of the two POs, the linear approximations of the corresponding manifolds and a point in a given Poincaré section where the unstable and stable manifolds match. The method is applied to compute families of heteroclinic orbits between planar Lyapunov POs around the collinear equilibrium points of the restricted three-body problem in different scenarios. In one of them, for the Sun–Jupiter mass parameter, we provide energy ranges for which the transition between different resonances is possible.

We characterize the family ${\mathcal G}_F$ of homogeneous iterated function systems of similarities on $\mathbb{R}^d$ that satisfy the strong separation condition and generate the same self-similar set F. We also show that ${\mathcal G}_F$ as a semigroup is finitely generated.

This paper presents a systematic existence and uniqueness theory of weak measure solutions for systems of non-local interaction PDEs with two species, which are the PDE counterpart of systems of deterministic interacting particles with two species. The main motivations behind those models arise in cell biology, pedestrian movements, and opinion formation. In case of symmetrizable systems (i.e. with cross-interaction potentials one multiple of the other), we provide a complete existence and uniqueness theory within (a suitable generalization of) the Wasserstein gradient flow theory in Ambrosio et al (2008 Gradient Flows in Metric Spaces and in the Space of Probability Measures (Lectures in Mathematics ETH Zürich) 2nd edn (Basel: Birkhäuser)) and Carrillo et al (2011 Duke Math. J.156 229–71), which allows the consideration of interaction potentials with a discontinuous gradient at the origin. In the general case of non-symmetrizable systems, we provide an existence result for measure solutions which uses a semi-implicit version of the Jordan–Kinderlehrer–Otto (JKO) scheme (Jordan et al 1998 SIAM J. Math. Anal.29 1–17), which holds in a reasonable non-smooth setting for the interaction potentials. Uniqueness in the non-symmetrizable case is proven for C² potentials using a variant of the method of characteristics.

A practical method is described for computing the unique generator of the algebra of first integrals associated with a large class of Hopf-zero singularity. The set of all volume-preserving classical normal forms of this singularity is introduced via a Lie algebra description. This is a maximal vector space of classical normal forms with first integral; this is whence our approach works. Systems with a nonzero condition on their quadratic parts are considered. The algebra of all first integrals for any such system has a unique (modulo scalar multiplication) generator. The infinite level volume-preserving parametric normal forms of any nondegenerate perturbation within the Lie algebra of any such system is computed, where it can have rich dynamics. The associated unique generator of the algebra of first integrals are derived. The symmetry group of the infinite level normal forms are also discussed. Some necessary formulas are derived and applied to appropriately modified Rössler and generalized Kuramoto–Sivashinsky equations to demonstrate the applicability of our theoretical results. An approach (introduced by Iooss and Lombardi) is applied to find an optimal truncation for the first level normal forms of these examples with exponentially small remainders. The numerically suggested radius of convergence (for the first integral) associated with a hypernormalization step is discussed for the truncated first level normal forms of the examples. This is achieved by an efficient implementation of the results using Maple.

We consider population dynamics models dN/dt = f(N(t_τ)) − d(N(t)) with an increasing fecundity function f and any mortality function d which can be quadratic, as in the logistic equation, or have a different form provided that the equation has at most one positive equilibrium. Here the delay in the production term can be distributed and unbounded. It is demonstrated that the positive equilibrium is globally attractive if it exists, otherwise all positive solutions tend to zero. Moreover, we demonstrate that solutions of the equation are intrinsically non-oscillatory: once the initial function is less/greater than the equilibrium K > 0, so is the solution for any positive time value. The assumptions on f, d and the delay are rather nonrestrictive, and several examples demonstrate that none of them can be omitted.

We prove that a Hamiltonian system $H\in C^2(M,\mathbb{R})$ is globally hyperbolic if any of the following statements hold: H is robustly topologically stable; H is stably shadowable; H is stably expansive; and H has the stable weak specification property. Moreover, we prove that, for a C²-generic Hamiltonian H, the union of the partially hyperbolic regular energy hypersurfaces and the closed elliptic orbits, forms a dense subset of M. As a consequence, any robustly transitive regular energy hypersurface of a C²-Hamiltonian is partially hyperbolic. Finally, we prove that stable weakly-shadowable regular energy hypersurfaces are partially hyperbolic.

A brake orbit is an orbit that starts with zero initial velocity. The purpose of this paper is to find periodic brake orbits in the isosceles three-body problem. We use various shooting arguments; we follow a curve of brake initial conditions under the flow until it reaches a suitable surface, and show that the image curve includes a point that corresponds to a periodic orbit. As a result, we prove the existence of six types of periodic brake orbits.

Near a charged surface, counterions of different valences and sizes cluster; and their concentration profiles stratify. At a distance from such a surface larger than the Debye length, the electric field is screened by counterions. Both recent studies using a variational mean-field approach that includes ionic size effects and Monte Carlo simulations suggest that counterion stratification is determined by the ionic valence-to-volume ratios. Central in the mean-field approach is a free-energy functional of ionic concentrations in which the ionic size effects are included through the entropic effect of solvent molecules. The corresponding equilibrium conditions define the generalized Boltzmann distributions relating the ionic concentrations to the electrostatic potential. This paper presents a detailed analysis and numerical calculations for such a free-energy functional to understand the dependence of the ionic charge density on the electrostatic potential through the generalized Boltzmann distributions, the role of ionic valence-to-volume ratios in the counterion stratification and the modification of Debye length due to the effect of ionic sizes.

The full text of this article is available in the PDF provided.

The full text of this article is available in the PDF provided.