Introduction to Focus Issue: Statistical mechanics and billiard-type dynamical systems

Dynamical systems of the billiard type are of fundamental importance for the description of numerous phenomena observed in many different fields of research, including statistical mechanics, Hamiltonian dynamics, nonlinear physics, and many others. This Focus Issue presents the recent progress in this area with contributions from the mathematical as well as physical stand point. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4730155]

A generalized Jaynes-Cummings Hamiltonian and supersymmetric shape invariance

A class of shape-invariant bound-state problems which represent two-level systems are introduced. It is shown that the coupled-channel Hamiltonians obtained correspond to the generalization of the Jaynes-Cummings Hamiltonian.

Critical exponents for a transition from integrability to non-integrability via localization of invariant tori in the Hamiltonian system

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

On the number of critical periods for planar polynomial systems

In this paper we get some lower bounds for the number of critical periods of families of centers which are perturbations of the linear one. We give a method which lets us prove that there are planar polynomial centers of degree l with at least 2[(l - 2)/2] critical periods as well as study concrete families of potential, reversible and Lienard centers. This last case is studied in more detail and we prove that the number of critical periods obtained with our approach does not. increases with the order of the perturbation. (C) 2007 Elsevier Ltd. All rights reserved.

Affine Toda systems coupled to matter fields

We investigate higher grading integrable generalizations of the affine Toda systems, where the flat connections defining the models take values in eigensubspaces of an integral gradation of an affine Kac-Moody algebra, with grades varying from l to -l (l > 1). The corresponding target space possesses nontrivial vacua and soliton configurations, which can be interpreted as particles of the theory, on the same footing as those associated to fundamental fields. The models can also be formulated by a hamiltonian reduction procedure from the so-called two-loop WZNW models. We construct the general solution and show the classes corresponding to the solitons. Some of the particles and solitons become massive when the conformal symmetry is spontaneously broken by a mechanism with an intriguing topological character and leading to a very simple mass formula. The massive fields associated to nonzero grade generators obey field equations of the Dirac type and may be regarded as matter fields. A special class of models is remarkable. These theories possess a U(1 ) Noether current, which, after a special gauge fixing of the conformal symmetry, is proportional to a topological current. This leads to the confinement of the matter field inside the solitons, which can be regarded as a one-dimensional bag model for QCD. These models are also relevant to the study of electron self-localization in (quasi-)one-dimensional electron-phonon systems.

Field-theoretic approach for systems of composite hadrons

Systems containing simultaneously hadrons and their constituents are most easily described by treating composite hadron field operators on the same kinematical footing as the constituent ones. Introduction of a unitary transformation allows redescription of hadrons by elementary-particle field operators. Transformation of the microscopic Hamiltonian leads to effective Hamiltonians describing all possible processes involving hadrons and their constituents.

Some comments on the bi(tri)-Hamiltonian structure of generalized AKNS and DNLS hierarchies

We give the correct prescriptions for the terms involving ∂ -1 xδ(x - y), in the Hamiltonian structures of the AKNS and DNLS systems, necessary for the Jacobi identities to hold. We establish that the sl(2) and sl(3) AKNS systems are tri-Hamiltonians and construct two compatible Hamiltonian structures for the sl(n) AKNS system. We give a method for the derivation of the recursion operator for the sl(n + 1) DNLS system, and apply it explicitly to the sl(2) case, showing that such a system is tri-Hamiltonian. © 1998 Elsevier Science B.V.

Analytical functions for the calculation of hyperspherical potential curves of atomic systems

We present angular basis functions for the Schrödinger equation of two-electron systems in hyperspherical coordinates. By using the hyperspherical adiabatic approach, the wave functions of two-electron systems are expanded in analytical functions, which generalizes the Jacobi polynomials. We show that these functions, obtained by selecting the diagonal terms of the angular equation, allow efficient diagonalization of the Hamiltonian for all values of the hyperspherical radius. The method is applied to the determination of the 1S e energy levels of the Li + and we show that the precision can be improved in a systematic and controllable way. ©2000 The American Physical Society.

Reversible Hamiltonian Liapunov center theorem

We study the existence of periodic solutions in the neighbourhood of symmetric (partially) elliptic equilibria in purely reversible Hamiltonian vector fields. These are Hamiltonian vector fields with an involutory reversing symmetry R. We contrast the cases where R acts symplectically and anti-symplectically. In case R acts anti-symplectically, generically purely imaginary eigenvalues are isolated, and the equilibrium is contained in a local two-dimensional invariant manifold containing symmetric periodic solutions encircling the equilibrium point. In case R acts symplectically, generically purely imaginary eigenvalues are doubly degenerate, and the equilibrium is contained in two two-dimensional invariant manifolds containing nonsymmetric periodic solutions encircling the equilibrium point. In addition, there exists a three-dimensional invariant surface containing a two-parameter family of symmetric periodic solutions.

Scaling invariance of the diffusion coefficient in a family of two-dimensional Hamiltonian mappings

We consider a family of two-dimensional nonlinear area-preserving mappings that generalize the Chirikov standard map and model a variety of periodically forced systems. The action variable diffuses in increments whose phase is controlled by a negative power of the action and hence effectively uncorrelated for small actions, leading to a chaotic sea in phase space. For larger values of the action the phase space is mixed and contains a family of elliptic islands centered on periodic orbits and invariant Kolmogorov-Arnold-Moser (KAM) curves. The transport of particles along the phase space is considered by starting an ensemble of particles with a very low action and letting them evolve in the phase until they reach a certain height h. For chaotic orbits below the periodic islands, the survival probability for the particles to reach h is characterized by an exponential function, well modeled by the solution of the diffusion equation. On the other hand, when h reaches the position of periodic islands, the diffusion slows markedly. We show that the diffusion coefficient is scaling invariant with respect to the control parameter of the mapping when h reaches the position of the lowest KAM island. © 2013 American Physical Society.

Multi-planet extrasolar systems - detection and dynamics

20 years after the discovery of the first planets outside our solar system, the current exoplanetary population includes more than 700 confirmed planets around main sequence stars. Approximately 50% belong to multiple-planet systems in very diverse dynamical configurations, from two-planet hierarchical systems to multiple resonances that could only have been attained as the consequence of a smooth large-scale orbital migration. The first part of this paper reviews the main detection techniques employed for the detection and orbital characterization of multiple-planet systems, from the (now) classical radial velocity (RV) method to the use of transit time variations (TTV) for the identification of additional planetary bodies orbiting the same star. In the second part we discuss the dynamical evolution of multi-planet systems due to their mutual gravitational interactions. We analyze possible modes of motion for hierarchical, secular or resonant configurations, and what stability criteria can be defined in each case. In some cases, the dynamics can be well approximated by simple analytical expressions for the Hamiltonian function, while other configurations can only be studied with semi-analytical or numerical tools. In particular, we show how mean-motion resonances can generate complex structures in the phase space where different libration islands and circulation domains are separated by chaotic layers. In all cases we use real exoplanetary systems as working examples.

Dynamics of the 3:1 resonant planetary systems

Many of the discovered exoplanetary systems are involved inside mean-motion resonances. In this work we focus on the dynamics of the 3:1 mean-motion resonant planetary systems. Our main purpose is to understand the dynamics in the vicinity of the apsidal corotation resonance (ACR) which are stationary solutions of the resonant problem. We apply the semi-analytical method (Michtchenko et al., 2006) to construct the averaged three-body Hamiltonian of a planetary system near a 3:1 resonance. Then we obtain the families of ACR, composed of symmetric and asymmetric solutions. Using the symmetric stable solutions we observe the law of structures (Ferraz-Mello,1988), for different mass ratio of the planets. We also study the evolution of the frequencies of σ1, resonant angle, and Δω, the secular angle. The resonant domains outside the immediate vicinity of ACR are studied using dynamical maps techniques. We compared the results obtained to planetary systems near a 3:1 MMR, namely 55 Cnc b-c, HD 60532 b-c and Kepler 20 b-c.

Study of the Influence of Localized Vibrational Modes in Charge Transport Properties at Nanoscale Systems.

In molecular and atomic devices the interaction between electrons and ionic vibrations has an important role in electronic transport. The electron-phonon coupling can cause the loss of the electron's phase coherence, the opening of new conductance channels and the suppression of purely elastic ones. From the technological viewpoint phonons might restrict the efficiency of electronic devices by energy dissipation, causing heating, power loss and instability. The state of the art in electron transport calculations consists in combining ab initio calculations via Density Functional Theory (DFT) with Non-Equilibrium Green's Function formalism (NEGF). In order to include electron-phonon interactions, one needs in principle to include a self-energy scattering term in the open system Hamiltonian which takes into account the effect of the phonons over the electrons and vice versa. Nevertheless this term could be obtained approximately by perturbative methods. In the First Born Approximation one considers only the first order terms of the electronic Green's function expansion. In the Self-Consistent Born Approximation, the interaction self-energy is calculated with the perturbed electronic Green's function in a self-consistent way. In this work we describe how to incorporate the electron-phonon interaction to the SMEAGOL program (Spin and Molecular Electronics in Atomically Generated Orbital Landscapes), an ab initio code for electronic transport based on the combination of DFT + NEGF. This provides a tool for calculating the transport properties of materials' specific system, particularly in molecular electronics. Preliminary results will be presented, showing the effects produced by considering the electron-phonon interaction in nanoscale devices.

Quantum Correlations in Field Theory and Integrable Systems

In this thesis we will investigate some properties of one-dimensional quantum systems. From a theoretical point of view quantum models in one dimension are particularly interesting because they are strongly interacting, since particles cannot avoid each other in their motion, and you we can never ignore collisions. Yet, integrable models often generate new and non-trivial solutions, which could not be found perturbatively. In this dissertation we shall focus on two important aspects of integrable one- dimensional models: Their entanglement properties at equilibrium and their dynamical correlators after a quantum quench. The first part of the thesis will be therefore devoted to the study of the entanglement entropy in one- dimensional integrable systems, with a special focus on the XYZ spin-1/2 chain, which, in addition to being integrable, is also an interacting model. We will derive its Renyi entropies in the thermodynamic limit and its behaviour in different phases and for different values of the mass-gap will be analysed. In the second part of the thesis we will instead study the dynamics of correlators after a quantum quench , which represent a powerful tool to measure how perturbations and signals propagate through a quantum chain. The emphasis will be on the Transverse Field Ising Chain and the O(3) non-linear sigma model, which will be both studied by means of a semi-classical approach. Moreover in the last chapter we will demonstrate a general result about the dynamics of correlation functions of local observables after a quantum quench in integrable systems. In particular we will show that if there are not long-range interactions in the final Hamiltonian, then the dynamics of the model (non equal- time correlations) is described by the same statistical ensemble that describes its statical properties (equal-time correlations).

Real-time simulation of large open quantum spin systems driven by dissipation

We consider a large quantum system with spins 12 whose dynamics is driven entirely by measurements of the total spin of spin pairs. This gives rise to a dissipative coupling to the environment. When one averages over the measurement results, the corresponding real-time path integral does not suffer from a sign problem. Using an efficient cluster algorithm, we study the real-time evolution from an initial antiferromagnetic state of the two-dimensional Heisenberg model, which is driven to a disordered phase, not by a Hamiltonian, but by sporadic measurements or by continuous Lindblad evolution.