Source author record

Holger R. Dullin

Holger R. Dullin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

27works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

27 published item(s)

preprint2021arXiv

The Harmonic Lagrange Top and the Confluent Heun Equation

The harmonic Lagrange top is the Lagrange top plus a quadratic (harmonic) potential term. We describe the top in the space fixed frame using a global description with a Poisson structure on $T^*S^3$. This global description naturally leads to a rational parametrisation of the set of critical values of the energy-momentum map. We show that there are 4 different topological types for generic parameter values. The quantum mechanics of the harmonic Lagrange top is described by the most general confluent Heun equation (also known as the generalised spheroidal wave equation). We derive formulas for an infinite pentadiagonal symmetric matrix representing the Hamiltonian from which the spectrum is computed.

preprint2020arXiv

Relative equilibria of the 3-body problem in $\mathbb{R}^4$

The classical equations of the Newtonian 3-body problem do not only define the familiar 3-dimensional motions. The dimension of the motion may also be 4, and cannot be higher. We prove that in dimension 4, for three arbitrary positive masses, and for an arbitrary value (of rank 4) of the angular momentum, the energy possesses a minimum, which corresponds to a motion of relative equilibrium which is Lyapunov stable when considered as an equilibrium of the reduced problem. The nearby motions are nonsingular and bounded for all time. We also describe the full family of relative equilibria, and show that its image by the energy-momentum map presents cusps and other interesting features.

preprint2019arXiv

On the $ C^{8/3} $-Regularisation of Simultaneous Binary Collisions in the Collinear 4-Body Problem

The singularity at a simultaneous binary collision is explored in the collinear 4-body problem. It is known that any attempt to remove the singularity via block regularisation will result in a regularised flow that is no more than $ C^{8/3} $ differentiable with respect to initial conditions. Through a blow-up of the singularity, this loss of differentiability is investigated and a new proof of the $ C^{8/3} $ regularity is provided. In the process, it is revealed that the collision manifold consists of two manifolds of normally hyperbolic saddle singularities which are connected by a manifold of heteroclinics. By utilising recent work on transitions near such objects and their normal forms, an asymptotic series of the transition past the singularity is explicitly computed. It becomes remarkably apparent that the finite differentiability at $ 8/3 $ is due to the inability to construct a set of integrals local to the simultaneous binary collision. The finite differentiability is shown to be independent from a choice of initial condition or value of the masses.

preprint2019arXiv

Poisson Structure of the Three-Dimensional Euler Equations in Fourier Space

We derive a simple Poisson structure in the space of Fourier modes for the vorticity formulation of the Euler equations on a three-dimensional periodic domain. This allows us to analyse the structure of the Euler equations using a Hamiltonian framework. The Poisson structure is valid on the divergence free subspace only, and we show that using a projection operator it can be extended to be valid in the full space. We then restrict the simple Poisson structure to the divergence-free subspace on which the dynamics of the Euler equations take place, reducing the size of the system of ODEs by a third. The projected and the restricted Poisson structures are shown to have the helicity as a Casimir invariant. We conclude by showing that periodic shear flows in three dimensions are equilibria that correspond to singular points of the projected Poisson structure, and hence that the usual approach to study their nonlinear stability through the Energy-Casimir method fails.

preprint2019arXiv

Stability Theory of the 3-Dimensional Euler Equations

The Euler equations on a three-dimensional periodic domain have a family of shear flow steady states. We show that the linearised system around these steady states decomposes into subsystems equivalent to the linearisation of shear flows in a two-dimensional periodic domain. To do so, we derive a formulation of the dynamics of the vorticity Fourier modes on a periodic domain and linearise around the shear flows. The linearised system has a decomposition analogous to the two-dimensional problem, which can be significantly simplified. By appealing to previous results it is shown that some subset of the shear flows are spectrally stable, and another subset are spectrally unstable. For a dense set of parameter values the linearised operator has a nilpotent part, leading to linear instability. This is connected to the nonnormality of the linearised dynamics and the transition to turbulence. Finally we show that all shear flows in the family considered (even the linearly stable ones) are parametrically unstable.

preprint2019arXiv

Symmetry reduction of the 3-body problem in $\mathbb{R}^4$

The 3-body problem in $\mathbb{R}^4$ has 24 dimensions and is invariant under translations and rotations. We do the full symplectic symmetry reduction and obtain a reduced Hamiltonian in local symplectic coordinates on a reduced phase space with 8 dimensions. The Hamiltonian depends on two parameters $μ_1 > μ_2 \ge 0$, related to the conserved angular momentum. The limit $μ_2 \to 0$ corresponds to the 3-dimensional limit. We show that the reduced Hamiltonian has relative equilibria that are local minima and hence Lyapunov stable when $μ_2$ is sufficiently small. This proves the existence of balls of initial conditions of full dimension that do not contain any orbits that are unbounded.

preprint2018arXiv

Quantum Monodromy in the Isotropic 3-Dimensional Harmonic Oscillator

The isotropic harmonic oscillator in dimension 3 separates in several different coordinate systems. Separating in a particular coordinate system defines a system of three commuting operators, one of which is the Hamiltonian. We show that the joint spectrum of the Hamilton operator, the $z$ component of the angular momentum, and a quartic integral obtained from separation in prolate spheroidal coordinates has quantum monodromy for sufficiently large energies. This means that one cannot globally assign quantum numbers to the joint spectrum. The effect can be classically explained by showing that the corresponding Liouville integrable system has a non-degenerate focus-focus point, and hence Hamiltonian monodromy.

preprint2018arXiv

Scattering invariants in Euler's two-center problem

The problem of two fixed centers was introduced by Euler as early as in 1760. It plays an important role both in celestial mechanics and in the microscopic world. In the present paper we study the spatial problem in the case of arbitrary (both positive and negative) strengths of the centers. Combining techniques from scattering theory and Liouville integrability, we show that this spatial problem has topologically non-trivial scattering dynamics, which we identify as scattering monodromy. The approach that we introduce in this paper applies more generally to scattering systems that are integrable in the Liouville sense.

preprint2018arXiv

Symplectic classification of coupled angular momenta

The coupled angular momenta are a family of completely integrable systems that depend on three parameters and have a compact phase space. They correspond to the classical version of the coupling of two quantum angular momenta and they constitute one of the fundamental examples of so-called semitoric systems. Pelayo & Vu Ngoc have given a classification of semitoric systems in terms of five symplectic invariants. Three of these invariants have already been partially calculated in the literature for a certain parameter range, together with the linear terms of the so-called Taylor series invariant for a fixed choice of parameter values. In the present paper we complete the classification by calculating the polygon invariant, the height invariant, the twisting-index invariant, and the higher-order terms of the Taylor series invariant for the whole family of systems. We also analyse the explicit dependence of the coefficients of the Taylor series with respect to the three parameters of the system, in particular near the Hopf bifurcation where the focus-focus point becomes degenerate.

preprint2018arXiv

Taylor series and twisting-index invariants of coupled spin-oscillators

About six years ago, semitoric systems on 4-dimensional manifolds were classified by Pelayo & Vu Ngoc by means of five invariants. A standard example of such a system is the coupled spin-oscillator on $\mathbb{S}^2 \times \mathbb{R}^2$. Calculations of three of the five semitoric invariants of this system (namely the number of focus-focus singularities, the generalised semitoric polygon, and the height invariant) already appeared in the literature, but the so-called twisting index was not yet computed and, of the so-called Taylor series invariant, only the linear terms were known. In the present paper, we complete the list of invariants for the coupled spin-oscillator by calculating higher order terms of the Taylor series invariant and by computing the twisting index. Moreover, we prove that the Taylor series invariant has certain symmetry properties that make the even powers in one of the variables vanish and allow us to show superintegrability of the coupled spin-oscillator on the zero energy level.

preprint2018arXiv

Using the Geometric Phase to Optimise Planar Somersaults

We derive the equations of motion for the planar somersault, which consist of two additive terms. The first is the dynamic phase that is proportional to the angular momentum, and the second is the geometric phase that is independent of angular momentum and depends solely on the details of the shape change. Next, we import digitised footage of an elite athlete performing 3.5 forward somersaults off the 3m springboard, and use the data to validate our model. We show that reversing and reordering certain sections of the digitised dive can maximise the geometric phase without affecting the dynamic phase, thereby increasing the overall rotation achieved. Finally, we propose a theoretical planar somersault consisting of four shape changing states, where the optimisation lies in finding the shape change strategy that maximises the overall rotation of the dive. This is achieved by balancing the rotational contributions from the dynamic and geometric phases, in which we show the geometric phase plays a small but important role in the optimisation process.

preprint2016arXiv

Instability of Equilibria for the 2D Euler Equations on the torus

We consider the hydrodynamics of an incompressible fluid on a 2D periodic domain. There exists a family of stationary solutions with vorticity given by $Ω^*=α\cos (\mathbf{p} \cdot \mathbf{x} )+β\sin (\mathbf{p} \cdot \mathbf{x} )$. This situation can be approximated as a structure preserving finite dimensional Hamiltonian system by a truncation introduced by Zeitlin (1990,2005) or by the more standard Galerkin style finite element method. We use these two truncations to analyse the linear stability of these solutions and analytical and numerical results are compared. Following the methods used by Li (2000) the problem is divided into subsystems and we prove that most subsystems are linearly stable. We derive a sufficient condition for a subsystem to be linearly unstable and derive an explicit lower bound for the associated real eigenvalues independent of the truncation size $N$. Then we show that the corresponding eigenvectors are in $\ell^2$. This together with known stability results for the 2D periodic Euler equations allows us to conclude that most of these stationary solutions are nonlinearly unstable. We confirm our results with a numerical computation of the spectrum for a large, finite truncation. Finally we discuss the essential spectrum of the full problem as the limit of the truncated problem.

preprint2013arXiv

Semi-global symplectic invariants of the Euler top

We compute the semi-global symplectic invariants near the hyperbolic equilibrium points of the Euler top. The Birkhoff normal form at the hyperbolic point is computed using Lie series. The actions near the hyperbolic point are found using Frobenius expansion of its Picard-Fuchs equation. We show that the Birkhoff normal form can also be found by inverting the regular solution of the Picard-Fuchs equation. Composition of the singular action integral with the Birkhoff normal form gives the semi-global symplectic invariant. Finally, we discuss the convergence of these invariants and show that in a neighbourhood of the separatrix the pendulum is not symplectically equivalent to any Euler top.

preprint2012arXiv

The degenerate C. Neumann system I: symmetry reduction and convexity

The C. Neumann system describes a particle on the sphere S^n under the influence of a potential that is a quadratic form. We study the case that the quadratic form has l+1 distinct eigenvalues with multiplicity. Each group of m_σequal eigenvalues gives rise to an O(m_σ)-symmetry in configuration space. The combined symmetry group G is a direct product of l+1 such factors, and its cotangent lift has an Ad^*-equivariant Momentum mapping. Regular reduction leads to the Rosochatius system on S^l, which has the same form as the Neumann system albeit for an additional effective potential. To understand how the reduced systems fit together we use singular reduction to construct an embedding of the reduced Poisson space T^*{S^n}/G into R^{3l+3}$. The global geometry is described, in particular the bundle structure that appears as a result of the superintegrability of the system. We show how the reduced Neumann system separates in elliptical-spherical co-ordinates. We derive the action variables and frequencies as complete hyperelliptic integrals of genus l. Finally we prove a convexity result for the image of the Casimir mapping restricted to the energy surface.

preprint2012arXiv

The Equilateral Pentagon at Zero Angular Momentum: Maximal Rotation Through Optimal Deformation

A pentagon in the plane with fixed side-lengths has a two-dimensional shape space. Considering the pentagon as a mechanical system with point masses at the corners we answer the question of how much the pentagon can rotate with zero angular momentum. We show that the shape space of the equilateral pentagon has genus 4 and find a fundamental region by discrete symmetry reduction with respect to symmetry group D_5. The amount of rotation Δθfor a loop in shape space at zero angular momentum is interpreted as a geometric phase and is obtained as an integral of a function B over the region of shape space enclosed by the loop. With a simple variational argument we determine locally optimal loops as the zero contours of the function B. The resulting shape change is represented as a Fourier series, and the global maximum of Δθ\approx 45\degree is found for a loop around the regular pentagram. We also show that restricting allowed shapes to convex pentagons the optimal loop is the boundary of the convex region and gives Δθ\approx 19\degree.

preprint2012arXiv

The Lie-Poisson structure of the reduced n-body problem

The classical n-body problem in d-dimensional space is invariant under the Galilean symmetry group. We reduce by this symmetry group using the method of polynomial invariants. As a result we obtain a reduced system with a Lie-Poisson structure which is isomorphic to sp(2n-2), independently of d. The reduction preserves the natural form of the Hamiltonian as a sum of kinetic energy that depends on velocities only and a potential that depends on positions only. Hence we proceed to construct a Poisson integrator for the reduced n-body problem using a splitting method.

preprint2011arXiv

Semi-global symplectic invariants of the spherical pendulum

We explicitly compute the semi-global symplectic invariants near the focus-focus point of the spherical pendulum. A modified Birkhoff normal form procedure is presented to compute the expansion of the Hamiltonian near the unstable equilibrium point in Eliasson-variables. Combining this with explicit formulas for the action we find the semi-global symplectic invariants near the focus-focus point introduced by Vu Ngoc 2003. We also show that the Birkhoff normal form is the inverse of a complete elliptic integral over a vanishing cycle. To our knowledge this is the first time that semi-global symplectic invariants near a focus-focus point have been computed explicitly. We close with some remarks about the pendulum, for which the invariants can be related to theta functions in a beautiful way.

preprint2009arXiv

Dynamics near the p : -q Resonance

We study the dynamics near the truncated p : +/- q resonant Hamiltonian equilibrium for p, q coprime. The critical values of the momentum map of the Liouville integrable system are found. The three basic objects reduced period, rotation number, and non-trivial action for the leading order dynamics are computed in terms of complete hyperelliptic integrals. A relation between the three functions that can be interpreted as a decomposition of the rotation number into geometric and dynamic phase is found. Using this relation we show that the p : -q resonance has fractional monodromy. Finally we prove that near the origin of the 1 : -q resonance the twist vanishes.

preprint2008arXiv

Quadratic Volume-Preserving Maps: Invariant Circles and Bifurcations

We study the dynamics of the five-parameter quadratic family of volume-preserving diffeomorphisms of R^3. This family is the unfolded normal form for a bifurcation of a fixed point with a triple-one multiplier and also is the general form of a quadratic three-dimensional map with a quadratic inverse. Much of the nontrivial dynamics of this map occurs when its two fixed points are saddle-foci with intersecting two-dimensional stable and unstable manifolds that bound a spherical ``vortex-bubble''. We show that this occurs near a saddle-center-Neimark-Sacker (SCNS) bifurcation that also creates, at least in its normal form, an elliptic invariant circle. We develop a simple algorithm to accurately compute these elliptic invariant circles and their longitudinal and transverse rotation numbers and use it to study their bifurcations, classifying them by the resonances between the rotation numbers. In particular, rational values of the longitudinal rotation number are shown to give rise to a string of pearls that creates multiple copies of the original spherical structure for an iterate of the map.

preprint2008arXiv

The 1:+2 / 1:-2 resonance

On the linear level elliptic equilibria of Hamiltonian systems are mere superpositions of harmonic oscillators. Non-linear terms can produce instability, if the ratio of frequencies is rational and the Hamiltonian is indefinite. In this paper we study the frequency ratio 1/2 and its unfolding. In particular we show that for the indefinite case (1:-2) the frequency ratio map in a neighbourhood of the origin has a critical point, i.e. the twist condition is violated for one torus on every energy surface near the energy of the equilibrium. In contrast, we show that the frequency map itself is non-degenerate (i.e. the Kolmogorov non-degeneracy condition holds) for every torus in a neighbourhood of the equilibrium point. As a byproduct we are able to obtain another proof of fractional monodromy in the 1:-2 resonance.

preprint2006arXiv

Geodesic flow on three dimensional ellipsoids with equal semi-axes

Following on from our previous study of the geodesic flow on three dimensional ellipsoid with equal middle semi-axes, here we study the remaining cases: Ellipsoids with two sets of equal semi-axes with $SO(2) \times SO(2)$ symmetry, ellipsoids with equal larger or smaller semi-axes with SO(2) symmetry, and ellipsoids with three semi-axes coinciding with SO(3) symmetry. All of these cases are Liouville-integrable, and reduction of the symmetry leads to singular reduced systems on lower-dimensional ellipsoids. The critical values of the energy-momentum maps and their singular fibers are completely classified. In the cases with SO(2) symmetry there are corank 1 degenerate critical points; all other critical points are non-degenreate. We show that in the case with $SO(2) \times SO(2)$ symmetry three global action variables exist and the image of the energy surface under the energy-momentum map is a convex polyhedron. The case with SO(3) symmetry is non-commutatively integrable, and we show that the fibers over regular points of the energy-casimir map are $T^2$ bundles over $S^2$.

preprint2006arXiv

Geodesics on the Ellipsoid and Monodromy

The equations for geodesic flow on the ellipsoid are well known, and were first solved by Jacobi in 1838 by separating the variables of the Hamilton-Jacobi equation. In 1979 Moser investigated the case of the general ellipsoid with distinct semi-axes and described a set of integrals which weren't know classically. After reviewing the properties of geodesic flow on the three dimensional ellipsoid with distinct semi-axes, we investigate the three dimensional ellipsoid with the two middle semi-axes being equal, corresponding to a Hamiltonian invariant under rotations. The system is Liouville-integrable and thus the invariant manifolds corresponding to regular points of the energy momentum map are 3-dimensional tori. An analysis of the critical points of the energy momentum maps gives the bifurcation diagram. We find the fibres of the critical values of the energy momentum map, and carry out an analysis of the action variables. We show that the obstruction to the existence of single valued globally smooth action variables is monodromy.

preprint2005arXiv

Extended phase diagram of the Lorenz model

The parameter dependence of the various attractive solutions of the three variable nonlinear Lorenz model equations for thermal convection in Rayleigh-Bénard flow is studied. Its bifurcation structure has commonly been investigated as a function of r, the normalized Rayleigh number, at fixed Prandtl number σ. The present work extends the analysis to the entire (r,σ) parameter plane. An onion like periodic pattern is found which is due to the alternating stability of symmetric and non-symmetric periodic orbits. This periodic pattern is explained by considering non-trivial limits of large r and σ. In addition to the limit which was previously analyzed by Sparrow, we identify two more distinct asymptotic regimes in which either σ/r or σ^2/r is constant. In both limits the dynamics is approximately described by Airy functions whence the periodicity in parameter space can be calculated analytically. Furthermore, some observations about sequences of bifurcations and coexistence of attractors, periodic as well as chaotic, are reported.

preprint2001arXiv

An Integrable Shallow Water Equation with Linear and Nonlinear Dispersion

We study a class of 1+1 quadratically nonlinear water wave equations that combines the linear dispersion of the Korteweg-deVries (KdV) equation with the nonlinear/nonlocal dispersion of the Camassa-Holm (CH) equation, yet still preserves integrability via the inverse scattering transform (IST) method. This IST-integrable class of equations contains both the KdV equation and the CH equation as limiting cases. It arises as the compatibility condition for a second order isospectral eigenvalue problem and a first order equation for the evolution of its eigenfunctions. This integrable equation is shown to be a shallow water wave equation derived by asymptotic expansion at one order higher approximation than KdV. We compare its traveling wave solutions to KdV solitons.