Escaping from nonhyperbolic chaotic attractors

The noise-induced escape process from a nonhyperbolic chaotic attractor is of physical and fundamental importance. We address this problem by uncovering the general mechanism of escape in the relevant low noise limit using the Hamiltonian theory of large fluctuations and by establishing the crucial role of the primary homoclinic tangency closest to the basin boundary in the dynamical process. In order to demonstrate that, we provide an unambiguous solution of the variational equations from the Hamiltonian theory. Our results are substantiated with the help of physical and dynamical paradigms, such as the Hénon and the Ikeda maps. It is further pointed out that our findings should be valid for driven flow systems and for experimental data.     

Exponential decay and scaling laws in noisy chaotic scattering

In this Letter we present a numerical study of the effect of noise on a chaotic scattering problem in open Hamiltonian systems. We use the second order Heun method for stochastic differential equations in order to integrate the equations of motion of a two-dimensional flow with additive white Gaussian noise. We use as a prototype model the paradigmatic Hénon-Heiles Hamiltonian with weak dissipation which is a well-known example of a system with escapes. We study the behavior of the scattering particles in the scattering region, finding an abrupt change of the decay law from algebraic to exponential due to the effects of noise. Moreover, we find a linear scaling law between the coefficient of the exponential law and the intensity of noise. These results are of a general nature in the sense that the same behavior appears when we choose as a model a two-dimensional discrete map with uniform noise (bounded in a particular interval and zero otherwise), showing the validity of the algorithm used. We believe the results of this work be useful for a better understanding of chaotic scattering in more realistic situations, where noise is presented.     

Scattering from correlated systems

A difficulty usually encountered in formulating the problem of scattering of identical particles from correlated systems is that the customary choice of an unperturbed Hamiltonian as the target Hamiltonian plus the kinetic energy of the projectile is not symmetric under particle exchange. This choice of unperturbed Hamiltonian leads to wavefunctions which, if they are antisymmetrized, are not orthonormal. In this paper an orthonormal, antisymmetrized set of basis states is constructed. These states are then used to construct a symmetric unperturbed Hamiltonian, so that a formal scattering equation with appropriate boundary conditions can be written. An expression for a T matrix describing nucleon-nucleus scattering can then be obtained. The formalism leads to a two-potential form for the T matrix, the first term of which describes the effect of the orthogonality of the scattering state and the negative energy states.     

Generalized variational principles,global weak solutions and behavior with random initial data for systems of conservation laws arising in adhesion particle dynamics

We study systems of conservation laws arising in two models of adhesion particle dynamics. The first is the system of free particles which stick under collision. The second is a system of gravitationally interacting particles which also stick under collision. In both cases, mass and momentum are conserved at the collisions, so the dynamics is described by 2×2 systems of conservations laws. We show that for these systems, global weak solutions can be constructed explicitly using the initial data by a procedure analogous to the Lax-Oleinik variational principle for scalar conservation laws. However, this weak solution is not unique among weak solutions satisfying the standard entropy condition. We also study a modified gravitational model in which, instead of momentum, some other weighted velocity is conserved at collisions. For this model, we prove both existence and uniqueness of global weak solutions. We then study the qualitative behavior of the solutions with random initial data. We show that for continuous but nowhere differentiable random initial velocities, all masses immediately concentrate on points even though they were continuously distributed initially, and the set of shock locations is dense.     

Chaotic and Arnold stripes in weakly chaotic Hamiltonian systems

The dynamics in weakly chaotic Hamiltonian systems strongly depends on initial conditions (ICs) and little can be affirmed about generic behaviors. Using two distinct Hamiltonian systems, namely one particle in an open rectangular billiard and four particles globally coupled on a discrete lattice, we show that in these models, the transition from integrable motion to weak chaos emerges via chaotic stripes as the nonlinear parameter is increased. The stripes represent intervals of initial conditions which generate chaotic trajectories and increase with the nonlinear parameter of the system. In the billiard case, the initial conditions are the injection angles. For higher-dimensional systems and small nonlinearities, the chaotic stripes are the initial condition inside which Arnold diffusion occurs.     

Logarithmic correction to the probability of capture for dissipatively perturbed Hamiltonian systems

Hamiltonian systems are analyzed with a double homoclinic orbit connecting a saddle to itself. Competing centers exist. A small dissipative perturbation causes the stable and unstable manifolds of the saddle point to break apart. The stable manifolds of the saddle point are the boundaries of the basin of attraction for the competing attractors. With small dissipation, the boundaries of the basins of attraction are known to be tightly wound and spiral-like. Small changes in the initial condition can alter the equilibrium to which the solution is attracted. Near the unperturbed homoclinic orbit, the boundary of the basin of attraction consists of a large sequence of nearly homoclinic orbits surrounded by close approaches to the saddle point. The slow passage through an unperturbed homoclinic orbit (separatrix) is determined by the change in the value of the Hamiltonian from one saddle approach to the next. The probability of capture can be asymptotically approximated using this change in the Hamiltonian. The well-known leading-order change of the Hamiltonian from one saddle approach to the next is due to the effect of the perturbation on the homoclinic orbit. A logarithmic correction to this change of the Hamiltonian is shown to be due to the effect of the perturbation on the saddle point itself. It is shown that the probability of capture can be significantly altered from the well-known leading-order probability for Hamiltonian systems with double homoclinic orbits of the twisted type, an example of which is the Hamiltonian system corresponding to primary resonance. Numerical integration of the perturbed Hamiltonian system is used to verify the accuracy of the analytic formulas for the change in the Hamiltonian from one saddle approach to the next. (c) 1995 American Institute of Physics.     

One-particle and few-particle billiards

We study the dynamics of one-particle and few-particle billiard systems in containers of various shapes. In few-particle systems, the particles collide elastically both against the boundary and against each other. In the one-particle case, we investigate the formation and destruction of resonance islands in (generalized) mushroom billiards, which are a recently discovered class of Hamiltonian systems with mixed regular-chaotic dynamics. In the few-particle case, we compare the dynamics in container geometries whose counterpart one-particle billiards are integrable, chaotic, and mixed. One of our findings is that two-, three-, and four-particle billiards confined to containers with integrable one-particle counterparts inherit some integrals of motion and exhibit a regular partition of phase space into ergodic components of positive measure. Therefore, the shape of a container matters not only for noninteracting particles but also for interacting particles.     

Slow passage through a transcritical bifurcation for Hamiltonian systems and the change in action due to a nonhyperbolic homoclinic orbit

One-degree of freedom conservative slowly varying Hamiltonian systems are analyzed in the case in which a saddle-center pair undergo a transcritical bifurcation. We analyze the case in which the method of averaging predicts the solution crosses the unperturbed homoclinic orbit at the precise time at which the transcritical bifurcation occurs. For the slow passage through the nonhyperbolic homoclinic orbit associated with a transcritical bifurcation, the solution consists of a large sequence of nonhyperbolic homoclinic orbits surrounded by autonomous nonlinear saddle approaches. The change in action is computed by matching these solutions to those obtained by averaging, valid before and after crossing the nonhyperbolic homoclinic orbit. For initial conditions near the stable manifold of the nonhyperbolic saddle point, one saddle approach has particularly small energy and instead satisfies a nonautonomous nonlinear equation, which provides a transition between nonhyperbolic homoclinic orbits, centers, and saddles. (c) 2000 American Institute of Physics.     

Macroscopic Diffusion from a Hamilton-like Dynamics

We introduce and analyze a model for the transport of particles or energy in extended lattice systems. The dynamics of the model acts on a discrete phase space at discrete times but has nonetheless some of the characteristic properties of Hamiltonian dynamics in a confined phase space: it is deterministic, periodic, reversible and conservative. Randomness enters the model as a way to model ignorance about initial conditions and interactions between the components of the system. The orbits of the particles are non-intersecting random loops. We prove, by a weak law of large number, the validity of a diffusion equation for the macroscopic observables of interest for times that are arbitrary large, but small compared to the minimal recurrence time of the dynamics.     

Quantum Dynamics in Noisy Backgrounds: from Sampling to Dissipation and Fluctuations

We investigate the dynamics of a quantum system coupled linearly to Gaussian white noise using functional methods. By performing the integration over the noisy field in the evolution operator, we get an equivalent non-Hermitian Hamiltonian, which evolves the quantum state with a dissipative dynamics. We also show that if the integration over the noisy field is done for the time evolution of the density matrix, a gain contribution from the fluctuations can be accessed in addition to the loss one from the non-hermitian Hamiltonian dynamics. We illustrate our study by computing analytically the effective non-Hermitian Hamiltonian, which we found to be the complex frequency harmonic oscillator, with a known evolution operator. It leads to space and time localisation, a common feature of noisy quantum systems in general applications.     

Vlasov moments, integrable systems and singular solutions

The Vlasov equation governs the evolution of the single-particle probability distribution function (PDF) for a system of particles interacting without dissipation. Its singular solutions correspond to the individual particle motions. The operation of taking the moments of the Vlasov equation is a Poisson map. The resulting Lie-Poisson Hamiltonian dynamics of the Vlasov moments is found to be integrable is several cases. For example, the dynamics for coasting beams in particle accelerators is associated by a hodograph transformation to the known integrable Benney shallow-water equation. After setting the context, the Letter focuses on geodesic Vlasov moment equations. Continuum closures of these equations at two different orders are found to be integrable systems whose singular solutions characterize the geodesic motion of the individual particles.     

Lyapunov Instability for a Hard-Disk Fluid in Equilibrium and Nonequilibrium Thermostated by Deterministic Scattering

We compute the full Lyapunov spectra for a hard-disk fluid under temperature gradient and under shear. The Lyapunov exponents are calculated using a recently developed formalism for systems with elastic hard collisions. The system is thermalized by deterministic and time-reversible scattering at the boundary, whereas the bulk dynamics remains Hamiltonian. This thermostating mechanism allows for energy fluctuations around a mean value which is reflected by only two vanishing Lyapunov exponents in equilibrium and nonequilibrium. In nonequilibrium steady states the phase-space volume is contracted on average, leading to a negative sum of the Lyapunov exponents. Since the system is driven inhomogeneously we do not expect the conjugate pairing rule to hold, which is indeed shown to be the case. Finally, the Kaplan–Yorke dimension and the Kolmogorov–Sinai entropy are calculated from the Lyapunov spectra.     

Comparison of several approaches to the relativistic dynamics of directly interacting particles

Several approaches to the relativistic dynamics of directly interacting particles are compared. The equivalence between constrained Hamiltonian relativistic systems and a priori Hamiltonian predictive ones is completely proved. Coordinate transformations are obtained to express these systems in the framework of noncovariant predictive mechanics. The world line condition for constrained Hamiltonian relativistic systems is analyzed and is proved to be also necessary in the predictive Hamiltonian framework.     

Can aerosols be trapped in open flows?

The fate of aerosols in open flows is relevant in a variety of physical contexts. Previous results are consistent with the assumption that such finite-size particles always escape in open chaotic advection. Here we show that a different behavior is possible. We analyze the dynamics of aerosols both in the absence and presence of gravitational effects, and both when the dynamics of the fluid particles is hyperbolic and nonhyperbolic. Permanent trapping of aerosols much heavier than the advecting fluid is shown to occur in all these cases. This phenomenon is determined by the occurrence of multiple vortices in the flow and is predicted to happen for realistic particle-fluid density ratios.     

Large Deviations for Countable to One Markov Systems

In this paper, we study large deviation properties for countable to one Markov systems associated to weak Gibbs measures for non-Hölder potentials. Furthermore, we establish multifractal large deviation laws for countable to one piecewise conformal Markov systems, which are derived systems constructed over hyperbolic regions for certain nonhyperbolic systems exhibiting intermittency. We apply our results to higher-dimensional number theoretical transformations.     

Bipartite entanglement induced by classically-constrained quantum dissipative dynamics

The properties of some complex many body systems can be modeled by introducing in the dissipative dynamics of each single component a set of kinetic constraints that depend on the state of the neighbor systems. Here, we characterize this kind of dynamics for two quantum systems whose independent dissipative evolutions are defined by a Lindblad equation. The constraints are introduced through a set of projectors that restrict the action of each single dissipative Lindblad channel to the state of the other system. Conditions that guarantee a classical interpretation of the kinetic constraints are found. The generation and evolution of entanglement is studied for two optical qubits systems. Classically constrained dissipation leads to a stationary state whose degree of entanglement depends on the initial state. Nevertheless, independently of the initial conditions, a maximal entangled state is generated when both systems are subjected to the action of local Hamiltonian fields that do not commutate with the constraints. The underlying physical mechanism is analyzed in detail.     

Extended Lagrangian formalisms for dyons and some applications to solid systems under external fields

We analyze the conditions of the electromagnetic potentials for systems with electric and magnetic charges and the Lagrangian theory with these potentials. The constructed Lagrangian function is valid for obtaining the field equations and the extended Lorentz force for dyonic charges for both relativistic particles in vacuum and non-relativistic entities in solids. In a second part, with the one-body Hamiltonian of independent particles in external fields, we explore some dual properties of the dyonic system under external fields. We analyze the possible diamagnetic (and ‘diaelectric’) response of magnetic monopoles under a weak and constant electromagnetic field and the theory of Landau levels in the case of magnetic charges under strong electromagnetic constant fields.     

Six-Bodies Calculations Using the Hyperspherical Harmonics Method

In this work we show results for light nuclear systems and small clusters of helium atoms using the hyperspherical harmonics basis. We use the basis without previous symmetrization or antisymmetrization of the state. After the diagonalization of the Hamiltonian matrix, the eigenvectors have well defined symmetry under particle permutation and the identification of the physical states is possible. We show results for systems composed up to six particles. As an example of a fermionic system, we consider a nucleon system interacting through the Volkov potential, used many times in the literature. For the case of bosons, we consider helium atoms interacting through a potential model which does not present a strong repulsion at short distances. We have used an attractive gaussian potential to reproduce the values of the dimer binding energy, the atom-atom scattering length, and the effective range obtained with one of the most widely used He–He interaction, the LM2M2 potential. In addition, we include a repulsive hypercentral three-body force to reproduce the trimer binding energy.