Source author record

Rafael de la Llave

Rafael de la Llave 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

33works
9topics
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

33 published item(s)

preprint2022arXiv

Lyapunov exponents everywhere and rigidity

In the present work we obtain rigidity results analysing the set of regular points, in the sense of Oseledec's Theorem. It is presented a study on the possibility of an Anosov diffeomorphisms having all Lyapunov exponents defined everywhere. We prove that this condition implies local rigidity of an Anosov automorphism of the torus $\mathbb{T}^d, d \geq 3,$ $C^1-$close to a linear automorphism diagonalizable over $\mathbb{R}$ and such that its characteristic polynomial is irreducible over $\mathbb{Q}.$

preprint2022arXiv

Nonconmutative coboundary equations over integrable systems

\def\G{\mathcal G} \def\M{\mathcal M} \def\cE{\mathcal E} We prove an analog of Livšic theorem for real-analytic families of cocycles over an integrable system with values in a Banach algebra $\G$ or a Lie group. Namely, we consider an integrable dynamical system $f:\M \equiv\torus^d \times [-1,1]^d\to \M$, $f(θ, I)=(θ+ I, I)$, and a real-analytic family of cocycles $η_\eps : \M \to \G$, indexed by a complex parameter $\eps$ in an open ball $\cE_ρ\in\CC$. We show that if $η_\eps$ has trivial periodic data, i.e., $$ η_\eps(f^{n-1}(p))\dots η_{\eps} (f(p))\cdot η_{\eps} (p)=Id $$ for each periodic point $p=f^n p$ and each $\eps \in \cE_ρ$, then there exists a real-analytic family of maps $ϕ_\eps: \M \to \G$ satisfying the coboundary equation $$ η_\eps(θ, I)=ϕ_\eps^{-1}\circ f(θ, I)\cdot ϕ_\eps (θ, I) $$ for all $(θ, I)\in \M$ and $\eps \in \cE_{ρ/2}$. We also show that if the coboundary equation above with an analytic left-hand side $η_\eps$ has a solution in the sense of formal power series in $\eps$, then it has an analytic solution.

preprint2022arXiv

Resonances and Phase Locking Phenomena for Foliation Preserving Torus Maps

It is well known for experts that resonances in nonlinear systems lead to new invariant objects that lead to new behaviors. The goal of this paper is to study the invariant sets generated by resonances under foliation preserving torus maps. That is torus which preserve a foliation of irrational lines $L_{θ_{0}}=\{θ_{0}+Ωt | t\in\mathbb{R}\}\subset\mathbb{T}^{d}$. Foliation preserving maps appear naturally as reparametrization of linear flows in the torus and also play an important role in several applications involving coupled oscillators, delay equations, resonators with moving walls, etc. The invariant objects we find here, lead to predictions on the behavior of these models. Since the results of this paper are meant to be applied for other problems, we have developed very quantitative results giving very explicit descriptions of the phenomena and the invariant objects that control them. The structure of the phase locking regions for foliation preserving maps is very different than for generic maps of the torus. Indeed, for the sake of completeness, we have developed similar analysis for the case of generic maps of the torus and shown that the objects that appear in foliation preserving maps are quantitatively and qualitatively different from those of generic torus maps. This has consequences in applications.

preprint2021arXiv

Melnikov method for non-conservative perturbations of the three-body problem

We consider the planar circular restricted three-body problem (PCRTBP), as a model for the motion of a spacecraft relative to the Earth-Moon system. We focus on the Lagrange equilibrium points $L_1$ and $L_2$. There are families of Lyapunov periodic orbits around either $L_1$ or $L_2$, forming Lyapunov manifolds. There also exist homoclinic orbits to the Lyapunov manifolds around either $L_1$ or $L_2$, as well as heteroclinic orbits between the Lyapunov manifold around $L_1$ and the one around $L_2$. The motion along the homoclinic/heteroclinic orbits can be described via the scattering map, which gives the future asymptotic of a homoclinic orbit as a function of the past asymptotic. In contrast with the more customary Melnikov theory, we do not need to assume that the asymptotic orbits have a special nature (periodic, quasi-periodic, etc.). We add a non-conservative, time-dependent perturbation, as a model for a thrust applied to the spacecraft for some duration of time, or for some other effect, such as solar radiation pressure. We compute the first order approximation of the perturbed scattering map, in terms of fast convergent integrals of the perturbation along homoclinic/heteroclinic orbits of the unperturbed system. As a possible application, this result can be used to determine the trajectory of the spacecraft upon using the thrust.

preprint2021arXiv

Persistence and Smooth Dependence on Parameters of Periodic Orbits in Functional Differential Equations Close to an ODE or an Evolutionary PDE

We consider functional differential equations(FDEs) which are perturbations of smooth ordinary differential equations(ODEs). The FDE can involve multiple state-dependent delays or distributed delays (forward or backward). We show that, under some mild assumptions, if the ODE has a nondegenerate periodic orbit, then the FDE has a smooth periodic orbit. Moreover, we get smooth dependence of the periodic orbit and its frequency on parameters with high regularity. The result also applies to FDEs which are perturbations of some evolutionary partial differential equations(PDEs). The proof consists in solving functional equations satisfied by the parameterization of the periodic orbit and the frequency using a fixed point approach. We do not need to consider the smoothness of the evolution or even the phase space of the FDEs.

preprint2021arXiv

Rapid and Accurate Methods for Computing Whiskered Tori and their Manifolds in Periodically Perturbed Planar Circular Restricted 3-Body Problems

When the planar circular restricted 3-body problem (RTBP) is periodically perturbed, families of unstable periodic orbits break up into whiskered tori, with most tori persisting into the perturbed system. In this study, we 1) develop a quasi-Newton method which simultaneously solves for the tori and their center, stable, and unstable directions; 2) implement continuation by both perturbation as well as rotation numbers; 3) compute Fourier-Taylor parameterizations of the stable and unstable manifolds; 4) regularize the equations of motion; and 5) globalize these manifolds. Our methodology improves on efficiency and accuracy compared to prior studies, and applies to a variety of periodic perturbations. We demonstrate the tools near resonances in the planar elliptic RTBP.

preprint2020arXiv

KAM estimates for the dissipative standard map

From the beginning of KAM theory, it was realized that its applicability to realistic problems depended on developing quantitative estimates on the sizes of the perturbations allowed. In this paper we present results on the existence of quasi-periodic solutions for conformally symplectic systems in non-perturbative regimes. We recall that, for conformally symplectic systems, finding the solution requires also to find a "drift parameter". We present a proof on the existence of solutions for values of the parameters which agree with more than three figures with the numerically conjectured optimal values. The verification of the estimates presented here is not completely rigorous since we do not control the round-off error. Nevertheless, running with different precision shows very little difference in the results. Given the high precision of the calculation and the simplicity of the estimates, this does not seem to affect the results. A full verification should be done implementing interval arithmetic. We make available the approximate solutions, the highly efficient algorithms to generate them (incorporating high precision based on the MPFR library) and the routines used to verify the applicability of the theorem.

preprint2020arXiv

KAM theory for some dissipative systems

Dissipative systems play a very important role in several physical models, most notably in Celestial Mechanics, where the dissipation drives the motion of natural and artificial satellites, leading them to migration of orbits, resonant states, etc. Hence the need to develop theories that ensure the existence of structures such as invariant tori or periodic orbits and device efficient computational methods. In this work we concentrate on the existence of invariant tori for the specific case of dissipative systems known as "conformally symplectic" systems, which have the property that they transform the symplectic form into a multiple of itself. To give explicit examples of conformally symplectic systems, we will present two different models: a discrete system, the standard map, and a continuous system, the spin-orbit problem. In both cases we will consider the conservative and dissipative versions, that will help to highlight the differences between the symplectic and conformally symplectic dynamics. For such dissipative systems we will present a KAM theorem in an a-posteriori format. The method of proof is based on extending geometric identities originally developed in [39] for the symplectic case. Besides leading to streamlined proofs of KAM theorem, this method provides a very efficient algorithm which has been implemented. Coupling an efficient numerical algorithm with an a-posteriori theorem, we have a very efficient way to provide rigorous estimates close to optimal. Indeed, the method gives a criterion (the Sobolev blow up criterion) that allows to compute numerically the breakdown. We will review this method as well as an extension of J. Greene's method and present the results in the conservative and dissipative standard maps. Computing close to the breakdown, allows to discover new mathematical phenomena such as the "bundle collapse mechanism".

preprint2020arXiv

Numerical computation of periodic orbits and isochrons for state-dependent delay perturbation of an ODE in the plane

We present algorithms and their implementation to compute limit cycles and their isochrons for state-dependent delay equations (SDDE's) which are perturbed from a planar differential equation with a limit cycle. Note that the space of solutions of an SDDE is infinite dimensional. We compute a two parameter family of solutions of the SDDE which converge to the solutions of the ODE as the perturbation goes to zero in a neighborhood of the limit cycle. The method we use formulates functional equations among periodic functions (or functions converging exponentially to periodic). The functional equations express that the functions solve the SDDE. Therefore, rather than evolving initial data and finding solutions of a certain shape, we consider spaces of functions with the desired shape and require that they are solutions. The mathematical theory of these invariance equations is developed in a companion paper, which develops "a posteriori" theorems. They show that, if there is a sufficiently approximate solution (with respect to some explicit condition numbers), then there is a true solution close to the approximate one. Since the numerical methods produce an approximate solution, and provide estimates of the condition numbers, we can make sure that the numerical solutions we consider approximate true solutions. In this paper, we choose a systematic way to approximate functions by a finite set of numbers (Taylor-Fourier series) and develop a toolkit of algorithms that implement the operators -- notably composition -- that enter into the theory. We also present several implementation results and present the results of running the algorithms and their implementation in some representative cases.

preprint2019arXiv

Existence of whiskered KAM tori of conformally symplectic systems

We study the existence of whiskered tori in a family $f_μ$ of conformally symplectic maps depending on parameters $μ$. Whiskered tori are tori on which the motion is a rotation, but they have as many expanding/contracting directions as allowed by the preservation of the geometric structure. Our main result is formulated in an "a-posteriori" format. We fix $ω$ satisfying Diophantine conditions. We assume that we are given 1) a value of the parameter $μ_0$, 2) an embedding of the torus $K_0$ into the phase space, approximately invariant under $f_{μ_0}$ in the sense that $f_{μ_0} \circ K_0 - K_0 \circ T_ω$ is small, 3) a splitting of the tangent space at the range of $K_0$, into three bundles which are approximately invariant under $D f_{μ_0}$ and such that the derivative satisfies "rate conditions" on each of the components. Then, if some non-degeneracy conditions are satisfied, we show that there is another parameter $μ_\infty$, an embedding $K_\infty$ and splittings close to the original ones which are invariant under $f_{μ_\infty}$. We also bound $|μ_\infty - μ_0|$, $\|K_\infty - K_0 \|$ and the distance of the initial and final splittings in terms of the initial error. The proof of the main theorem consists in describing an iterative process that takes advantage of cancellations coming from the geometry. Then, we show that the process converges to a true solution when started from an approximate enough solution. The iterative process leads to an efficient algorithm that is quite practical to implement. As an application, we study the singular problem of effects of small dissipation on whiskered tori. We develop formal expansions in the perturbative parameter and use them as input for the a-posteriori theorem. This allows to obtain lower bounds for the domain of analyticity of the tori as function of the perturbative parameter.

preprint2019arXiv

Expansions in the delay of quasi-periodic solutions for state dependent delay equations

We consider several models of State Dependent Delay Differential Equations (SDDEs), in which the delay is affected by a small parameter. This is a very singular perturbation since the nature of the equation changes. Under some conditions, we construct formal power series, which solve the SDDEs order by order. These series are quasi-periodic functions of time. This is very similar to the Lindstedt procedure in celestial mechanics. Truncations of these power series can be taken as input for a-posteriori theorems, that show that near the approximate solutions there are true solutions. In this way, we hope that one can construct a catalogue of solutions for SDDEs, bypassing the need of a systematic theory of existence and uniqueness for all initial conditions.

preprint2016arXiv

A Framework for the Numerical Computation and a Posteriori Verification of Invariant Objects of Evolution Equations

We develop a theoretical framework for computer-assisted proofs of the existence of invariant objects in semilinear PDEs. The invariant objects considered in this paper are equilibrium points, traveling waves, periodic orbits and invariant manifolds attached to fixed points or periodic orbits. The core of the study is writing down the invariance condition as a zero of an operator. These operators are in general not continuous, so one needs to smooth them by means of preconditioners before classical fixed point theorems can be applied. We develop in detail all the aspects of how to work with these objects: how to precondition the equations, how to work with the nonlinear terms, which function spaces can be useful, and how to work with them in a computationally rigorous way. In two companion papers, we present two different implementations of the tools developed in this paper to study periodic orbits.

preprint2016arXiv

An a posteriori KAM theorem for whiskered tori in Hamiltonian partial differential equations with applications to some ill-posed equations

The goal of this paper is to develop a KAM theory for tori with hyperbolic directions, which applies to Hamiltonian partial differential equations, even to some ill-posed ones. The main result has an \emph{a-posteriori} format, i.e., we show that if there is an approximate solution of an invariance equation which also satisfies some non-degeneracy conditions, then there is a true solution nearby. The a-posteriori format also has other automatic consequences (smooth dependence on parameters, bootstrap of regularity, etc.). The method of proof is based on an iterative method to solve a functional equation for the parameterization of the torus satisfying the invariance equations and for parametrization of directions invariant under the linearizatation. The iterative method does not use transformation theory or action-angle variables. It does not assume that the system is close to integrable. We first develop an abstract theorem. Then, we show how this abstract result applies to some concrete examples, including the scalar Boussinesq equation and the Boussinesq system so that we construct {\sl small amplitude} tori for the equations, which are even in the spatial variable. Note that the equations we use as examples are ill-posed.

preprint2016arXiv

Arnold diffusion in the planar elliptic restricted three-body problem: mechanism and numerical verification

We present a diffusion mechanism for time-dependent perturbations of autonomous Hamiltonian systems introduced in [25]. This mechanism is based on shadowing of pseudo-orbits generated by two dynamics: an `outer dynamics', given by homoclinic trajectories to a normally hyperbolic invariant manifold, and an `inner dynamics', given by the restriction to that manifold. On the inner dynamics the only assumption is that it preserves area. Unlike other approaches, [25] does not rely on the KAM theory and/or Aubry-Mather theory to establish the existence of diffusion. Moreover, it does not require to check twist conditions or non-degeneracy conditions near resonances. The conditions are explicit and can be checked by finite precision calculations in concrete systems. As an application, we study the planar elliptic restricted three-body problem. We present a rigorous theorem that shows that if some concrete calculations yield a non zero value, then for any sufficiently small, positive value of the eccentricity of the orbits of the main bodies, there are orbits of the infinitesimal body that exhibit a change of energy that is bigger than some fixed number, which is independent of the eccentricity. We verify numerically these calculations for values of the masses close to that of the Jupiter/Sun system. The numerical calculations are not completely rigorous, because we ignore issues of round-off error and do not estimate the truncations, but they are not delicate at all by the standard of numerical analysis. (Standard tests indicate that we get 7 or 8 figures of accuracy where 1 would be enough). The code of this verifications is available. We hope that some full computer assisted proofs will be obtained in a near future since there are packages (CAPD) designed for problems of this type.

preprint2016arXiv

Numerical Computations and Computer Assisted Proofs of Periodic Orbits of the Kuramoto-Sivashinsky Equation

We present numerical results and computer assisted proofs of the existence of periodic orbits for the Kuramoto-Sivashinky equation. These two results are based on writing down the existence of periodic orbits as zeros of functionals. This leads to the use of Newton's algorithm for the numerical computation of the solutions and, with some a posteriori analysis in combination with rigorous interval arithmetic, to the rigorous verification of the existence of solutions. This is a particular case of the methodology developed in [19] for several types of orbits. An independent implementation, covering overlapping but different ground, using different functional setups appears in [33].

preprint2015arXiv

A continuous family of equilibria in ferromagnetic media are ground states

We show that a foliation of equilibria (a continuous family of equilibria whose graph covers all the configuration space) in ferromagnetic models are ground states. The result we prove is very general, and it applies to models with long range interactions and many body. As an application, we consider several models of networks of interacting particles including models of Frenkel-Kontorova type on $\mathbb{Z}^d$ and one-dimensional quasi-periodic media. The result above is an analogue of several results in the calculus variations (fields of extremals) and in PDE's. Since the models we consider are discrete and long range, new proofs need to be given. We also note that the main hypothesis of our result (the existence of foliations of equilibria) is the conclusion (using KAM theory) of several recent papers. Hence, we obtain that the KAM solutions recently established are minimizers when the interaction is ferromagnetic (and transitive).

preprint2015arXiv

Domains of analyticity of Lindstedt expansions of KAM tori in dissipative perturbations of Hamiltonian systems

Many problems in Physics are described by dynamical systems that are conformally symplectic (e.g., mechanical systems with a friction proportional to the velocity, variational problems with a small discount or thermostated systems). Conformally symplectic systems are characterized by the property that they transform a symplectic form into a multiple of itself. The limit of small dissipation, which is the object of the present study, is particularly interesting. We provide all details for maps, but we present also the modifications needed to obtain a direct proof for the case of differential equations. We consider a family of conformally symplectic maps $f_{μ, ε}$ defined on a $2d$-dimensional symplectic manifold $\mathcal M$ with exact symplectic form $Ω$; we assume that $f_{μ,ε}$ satisfies $f_{μ,ε}^*Ω=λ(ε) Ω$. We assume that the family depends on a $d$-dimensional parameter $μ$ (called drift) and also on a small scalar parameter $ε$. Furthermore, we assume that the conformal factor $λ$ depends on $ε$, in such a way that for $ε=0$ we have $λ(0)=1$ (the symplectic case). We study the domains of analyticity in $ε$ near $ε=0$ of perturbative expansions (Lindstedt series) of the parameterization of the quasi--periodic orbits of frequency $ω$ (assumed to be Diophantine) and of the parameter $μ$. Notice that this is a singular perturbation, since any friction (no matter how small) reduces the set of quasi-periodic solutions in the system. We prove that the Lindstedt series are analytic in a domain in the complex $ε$ plane, which is obtained by taking from a ball centered at zero a sequence of smaller balls with center along smooth lines going through the origin. The radii of the excluded balls decrease faster than any power of the distance of the center to the origin.

preprint2015arXiv

Resonant equilibrium configurations in quasi-periodic media: KAM theory

We develop an a-posteriori KAM theory for the equilibrium equations for quasi-periodic solutions in a quasi-periodic Frenkel-Kontorova model when the frequency of the solutions resonates with the frequencies of the substratum. The KAM theory we develop is very different both in the methods and in the conclusions from the more customary KAM theory for Hamiltonian systems or from the KAM theory in quasi-periodic media for solutions with frequencies which are Diophantine with respect to the frequencies of the media. The main difficulty is that we cannot use transformations (as in the Hamiltonian case) nor Ward identities (as in the case of frequencies Diophantine with those of the media). The technique we use is to add an extra equation that ensures the linearization of the equilibrium equation factorizes. To solve the extra equation requires an extra counterterm. We compare this phenomenon with other phenomena in KAM theory. It seems that this technique could be used in several other problems. As a conclusion, we obtain that the perturbation expansions developed in the previous paper \cite{SuZL15} converge when the potentials are in a codimension one manifold in a space of potentials. The method of proof also leads to efficient (low storage requirements and low operation count) algorithms to compute the quasi-periodic solutions.

preprint2015arXiv

Resonant Equilibrium configurations in quasi-periodic media: perturbative expansions

We consider 1-D quasi-periodic Frenkel-Kontorova models. We study the existence of equilibria whose frequency (i.e. the inverse of the density of deposited material) is resonant with the frequencies of the substratum. We study perturbation theory for small potential. We show that there are perturbative expansions to all orders for the quasi-periodic equilibria with resonant frequencies. Under very general conditions, we show that there are at least two such perturbative expansions for equilibria for small values of the parameter. We also develop a dynamical interpretation of the equilibria in these quasi-periodic media. We show that equilibria are orbits of a dynamical system which has very unusual properties. We obtain results on the Lyapunov exponents of the dynamical systems, i.e. the phonon gap of the resonant quasi-periodic equilibria. We show that the equilibria can be pinned even if the gap is zero.

preprint2015arXiv

Response solutions for quasi-periodically forced, dissipative wave equations

We consider several models of nonlinear wave equations subject to very strong damping and quasi-periodic external forcing. This is a singular perturbation, since the damping is not the highest order term. We study the existence of response solutions (i.e., quasi-periodic solutions with the same frequency as the forcing). Under very general non-resonance conditions on the frequency, we show the existence of asymptotic expansions of the response solution; moreover, we prove that the response solution indeed exists and depends analytically on $\varepsilon$ (where $\varepsilon$ is the inverse of the coefficient multiplying the damping) for $\varepsilon$ in a complex domain, which in some cases includes disks tangent to the imaginary axis at the origin. In other models, we prove analyticity in cones of aperture $π/2$ and we conjecture it is optimal. These results have consequences for the asymptotic expansions of the response solutions considered in the literature. The proof of our results relies on reformulating the problem as a fixed point problem, constructing an approximate solution and studying the properties of iterations that converge to the solutions of the fixed point problem.

preprint2014arXiv

Barriers to Transport and Mixing in Volume-Preserving Maps with Nonzero Flux

In this paper we identify the geometric structures that restrict transport and mixing in perturbations of integrable volume-preserving systems with nonzero net flux. Unlike KAM tori, these objects cannot be continued to the tori present in the integrable system but are generated by resonance and have a contractible direction. We introduce a remarkably simple algorithm to analyze the behavior of these maps and obtain quantitative properties of the tori. In particular, we present assertions regarding the distribution of the escape times of the unbounded orbits, the abundance of tori, and the size of the resonant regions.

preprint2014arXiv

Construction of invariant whiskered tori by a parameterization method. Part II: Quasi-periodic and almost periodic breathers in coupled map lattices

We construct quasi-periodic and almost periodic solutions for coupled Hamiltonian systems on an infinite lattice which is translation invariant. The couplings can be long range, provided that they decay moderately fast with respect to the distance. For the solutions we construct, most of the sites are moving in a neighborhood of a hyperbolic fixed point, but there are oscillating sites clustered around a sequence of nodes. The amplitude of these oscillations does not need to tend to zero. In particular, the almost periodic solutions do not decay at infinity. We formulate an invariance equation. Solutions of this equation are embeddings of an invariant torus on which the motion is conjugate to a rotation. We show that, if there is an approximate solution of the invariance equation that satisfies some non-degeneracy conditions, there is a true solution close by. The proof of this \emph{a-posteriori} theorem is based on a Nash-Moser iteration, which does not use transformation theory. Simpler versions of the scheme were developed in E. Fontich, R. de la Llave,Y. Sire \emph{J. Differential. Equations.} {\bf 246}, 3136 (2009). One technical tool, important for our purposes, is the use of weighted spaces that capture the idea that the maps under consideration are local interactions. Using these weighted spaces, the estimates of iterative steps are similar to those in finite dimensional spaces. In particular, the estimates are independent of the number of nodes that get excited. Using these techniques, given two breathers, we can place them apart and obtain an approximate solution, which leads to a true solution nearby. By repeating the process infinitely often, we can get solutions with infinitely many frequencies which do not tend to zero at infinity.

preprint2013arXiv

Instability of high dimensional Hamiltonian Systems: Multiple resonances do not impede diffusion

We consider models given by Hamiltonians of the form $$H(I,ϕ,p,q,t;ε) = h(I) + \sum_{j = 1}^n \pm(\frac{1}{2} p_j^2 + V_j(q_j)) + εQ(I,ϕ,p,q,t;ε)$$ where $I,ϕ$ are d-dimensional actions and angles, $p,q$ are n-dimensional real conjugated variables, and $t$ is an angle. These are higher dimensional analogues, both in the center and hyperbolic directions, of the models studied in previous papers by the athors. All these models present the large gap problem. We show that, for $ε$ small enough, under regularity and explicit non-degeneracy conditions on the model, there are orbits whose action variables $I$ perform rather arbitrary excursions in a domain of size O(1). This domain includes resonance lines and, hence, large gaps among $d$-dimensional KAM tori. The method of proof follows closely the strategy of [DLS(=Delshams, Llave and Seara) 2006]. The main new phenomenon that appears when the dimension $d$ of the center directions is larger than one, is the existence of multiple resonances. We show that, since these multiple resonances happen in sets of codimension greater than one in the space of actions $I$, they can be contoured. This corresponds to the mechanism called diffusion across resonances in the Physics literature. The present paper, however, differs substantially from [DLS 2006]. On the technical details of the proofs, we have taken advantage of the theory of the scattering map [DLS 2008], not available when the above papers were written. We have analyzed the conditions imposed on the resonances in more detail. More precisely, we have found that there is a simple condition on the Melnikov potential which allows us to conclude that the resonances are crossed. In particular, this condition does not depend on the resonances. So that the results are new even when applied to the models in [DLS 2006].

preprint2013arXiv

Localized Stable Manifolds for Whiskered Tori in Coupled Map Lattices with Decaying Interaction

In this paper we consider lattice systems coupled by local interactions. We prove invariant manifold theorems for whiskered tori (we recall that whiskered tori are quasi-periodic solutions with exponentially contracting and expanding directions in the linearized system). The invariant manifolds we construct generalize the usual (strong) (un) stable manifolds and allow us to consider also non-resonant manifolds. We show that if the whiskered tori are localized near a collection of specific sites, then so are the invariant manifolds. We recall that the existence of localized whiskered tori has recently been proven for symplectic maps and flows in Fontich et. al. (submitted), but our results do not need that the systems are symplectic. For simplicity we will present first the main results for maps, but we will show tha the result for maps imply the results for flows. It is also true that the results for flows can be proved directly following the same ideas.

preprint2013arXiv

Perturbations of geodesic flows by recurrent dynamics

We consider a geodesic flow on a compact manifold endowed with a Riemannian (or Finsler, or Lorentz) metric satisfying some generic, explicit conditions. We couple the geodesic flow with a time-dependent potential, driven by an external flow on some other compact manifold. If the external flow satisfies some very general recurrence condition, and the potential satisfies some explicit conditions that are also very general, we show that the coupled system has orbits whose energy grows at a linear rate with respect to time. This growth rate is optimal. We also show the existence of symbolic dynamics. The existence of orbits whose energy grows unboundedly in time is related to Arnold's diffusion problem. A particular case of this phenomenon is obtained when the external dynamical system is quasi-periodic, of rationally independent frequency vector, not necessarily Diophantine, thus extending earlier results. We also recover `Mather's acceleration theorem' for time-periodic perturbations. Since the general class of perturbations that we consider in this paper are not characterized by a frequency of motion, we are not able to use KAM or Aubry-Mather or averaging theory. We devise a new mechanism to construct orbits whose energy evolves in a prescribed fashion. This mechanism intertwines the inner dynamics restricted to some normally hyperbolic manifold with the outer dynamics corresponding to two distinct choices of homoclinic excursions to that manifold. As this mechanism uses very coarse information on the dynamics, it is potentially of wider applicability.

preprint2013arXiv

The Analyticity Breakdown for Frenkel-Kontorova Models in Quasi-periodic Media: Numerical Explorations

We study numerically the "analyticity breakdown" transition in 1-dimensional quasi-periodic media. This transition corresponds physically to the transition between pinned down and sliding ground states. Mathematically, it corresponds to the solutions of a functional equation losing their analyticity properties. We implemented some recent numerical algorithms that are efficient and backed up by rigorous results so that we can compute with confidence even close to the breakdown. We have uncovered several phenomena that we believe deserve a theoretical explanation: A) The transition happens in a smooth surface. B) There are scaling relations near breakdown. C) The scaling near breakdown is very anisotropic. Derivatives in different directions blow up at different rates. Similar phenomena seem to happen in other KAM problems.

preprint2012arXiv

Invariant manifolds for analytic difference equations

We use a modification of the parameterization method to study invariant manifolds for difference equations. We establish existence, regularity, smooth dependence on parameters and study several singular limits, even if the difference equations do not define a dynamical system. This method also leads to efficient algorithms that we present with their implementations. The manifolds we consider include not only the classical strong stable and unstable manifolds but also manifolds associated to non-resonant spaces. When the difference equations are the Euler-Lagrange equations of a discrete variational we present sharper results. Note that, if the Legendre condition fails, the Euler-Lagrange equations can not be treated as a dynamical system. If the Legendre condition becomes singular, the dynamical system may be singular while the difference equation remains regular. We present numerical applications to several examples in the physics literature: the Frenkel-Kontorova model with long-range interactions and the Heisenberg model of spin chains with a perturbation. We also present extensions to finite differentiable difference equations.

preprint2012arXiv

Tracing KAM tori in presymplectic dynamical systems

We present a KAM theorem for presymplectic dynamical systems. The theorem has a " a posteriori " format. We show that given a Diophantine frequency $ω$ and a family of presymplectic mappings, if we find an embedded torus which is approximately invariant with rotation $ω$ such that the torus and the family of mappings satisfy some explicit non-degeneracy condition, then we can find an embedded torus and a value of the parameter close to to the original ones so that the torus is invariant under the map associated to the value of the parameter. Furthermore, we show that the dimension of the parameter space is reduced if we assume that the systems are exact.

preprint2011arXiv

A Comparison Principle for a Sobolev Gradient Semi-Flow

We consider gradient descent equations for energy functionals of the type S(u) = 1/2 < u(x), A(x)u(x) >_{L^2} + \int_Ω V(x,u) dx, where A is a uniformly elliptic operator of order 2, with smooth coefficients. The gradient descent equation for such a functional depends on the metric under consideration. We consider the steepest descent equation for S where the gradient is an element of the Sobolev space H^β, β\in (0,1), with a metric that depends on A and a positive number γ> \sup |V_{22}|. We prove a weak comparison principle for such a gradient flow. We extend our methods to the case where A is a fractional power of an elliptic operator. We provide an application to the Aubry-Mather theory for partial differential equations and pseudo-differential equations by finding plane-like minimizers of the energy functional

preprint2011arXiv

KAM theory for quasi-periodic equilibria in 1-D quasiperiodic media

We consider Frenkel-Kontorova models corresponding to 1 dimensional quasicrystals. We present a KAM theory for quasi-periodic equilibria. The theorem presented has an \emph{a-posteriori} format. We show that, given an approximate solution of the equilibrium equation, which satisfies some appropriate non-degeneracy conditions, then, there is a true solution nearby. This solution is locally unique. Such a-posteriori theorems can be used to validate numerical computations and also lead immediately to several consequences a) Existence to all orders of perturbative expansion and their convergence b) Bootstrap for regularity c) An efficient method to compute the breakdown of analyticity. Since the system does not admit an easy dynamical formulation, the method of proof is based on developing several identities. These identities also lead to very efficient algorithms.

preprint2011arXiv

Perturbation and Numerical Methods for Computing the Minimal Average Energy

We investigate the differentiability of minimal average energy associated to the functionals $S_\ep (u) = \int_{\mathbb{R}^d} (1/2)|\nabla u|^2 + \ep V(x,u)\, dx$, using numerical and perturbative methods. We use the Sobolev gradient descent method as a numerical tool to compute solutions of the Euler-Lagrange equations with some periodicity conditions; this is the cell problem in homogenization. We use these solutions to determine the average minimal energy as a function of the slope. We also obtain a representation of the solutions to the Euler-Lagrange equations as a Lindstedt series in the perturbation parameter $\ep$, and use this to confirm our numerical results. Additionally, we prove convergence of the Lindstedt series.

preprint2010arXiv

Computation of whiskered invariant tori and their associated manifolds: new fast algorithms

In this paper we present efficient algorithms for the computation of several invariant objects for Hamiltonian dynamics. More precisely, we consider KAM tori (i.e diffeomorphic copies of the torus such that the motion on them is conjugated to a rigid rotation) both Lagrangian tori (of maximal dimension) and whiskered tori (i.e. tori with hyperbolic directions which, together with the tangents to the torus and the symplectic conjugates span the whole tangent space). In the case of whiskered tori, we also present algorithms to compute the invariant splitting and the invariant manifolds associated to the splitting. We present them both for the case of discrete time and for differential equations. The algorithms are based on a Newton method to solve an appropriately chosen functional equation that expresses invariance. The algorithms are efficient: if we discretize the objects by $N$ elements, one step of the Newton method requires only O(N) storage and $O(N \ln(N))$ operations. Furthermore, if the object we consider is of dimension $\ell$, we only need to compute functions of $\ell$ variables, independently of what is the dimension of the phase space. The algorithms do not require that the system is presented in action-angle variables nor that it is close to integrable. The algorithms are backed up by rigorous \emph{a-posteriori} bounds which state that if the equations are solved with a small residual and some explicitly computable condition numbers are not too big, then, there is a true solution which is close to the computed one. The algorithms apply both to primary (i.e non-contractible) and secondary tori (i.e. contractible to a torus of lower dimension, such as islands). They have already been implemented. We will report on the technicalities of the implementation and the results of running them elsewhere.