Researcher profile

Rafael de la Llave

Rafael de la Llave contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

11 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.