455 research outputs found
Lyapunov exponents from geodesic spread in configuration space
The exact form of the Jacobi -- Levi-Civita (JLC) equation for geodesic
spread is here explicitly worked out at arbitrary dimension for the
configuration space manifold M_E = {q in R^N | V(q) < E} of a standard
Hamiltonian system, equipped with the Jacobi (or kinetic energy) metric g_J. As
the Hamiltonian flow corresponds to a geodesic flow on (M_E,g_J), the JLC
equation can be used to study the degree of instability of the Hamiltonian
flow. It is found that the solutions of the JLC equation are closely resembling
the solutions of the standard tangent dynamics equation which is used to
compute Lyapunov exponents. Therefore the instability exponents obtained
through the JLC equation are in perfect quantitative agreement with usual
Lyapunov exponents. This work completes a previous investigation that was
limited only to two-degrees of freedom systems.Comment: REVTEX file, 10 pages, 2 figure
Topological aspects of geometrical signatures of phase transitions
Certain geometric properties of submanifolds of configuration space are
numerically investigated for classical lattice phi^4 models in one and two
dimensions. Peculiar behaviors of the computed geometric quantities are found
only in the two-dimensional case, when a phase transition is present. The
observed phenomenology strongly supports, though in an indirect way, a recently
proposed topological conjecture about a topology change of the configuration
space submanifolds as counterpart of a phase transition.Comment: REVTEX file, 4 pages, 5 figure
Macroscopic detection of the strong stochasticity threshold in Fermi-Pasta-Ulam chains of oscillators
The largest Lyapunov exponent of a system composed by a heavy impurity
embedded in a chain of anharmonic nearest-neighbor Fermi-Pasta-Ulam oscillators
is numerically computed for various values of the impurity mass . A
crossover between weak and strong chaos is obtained at the same value
of the energy density (energy per degree of freedom)
for all the considered values of the impurity mass . The threshold \epsi
lon_{_T} coincides with the value of the energy density at which a
change of scaling of the relaxation time of the momentum autocorrelation
function of the impurity ocurrs and that was obtained in a previous work ~[M.
Romero-Bastida and E. Braun, Phys. Rev. E {\bf65}, 036228 (2002)]. The complete
Lyapunov spectrum does not depend significantly on the impurity mass . These
results suggest that the impurity does not contribute significantly to the
dynamical instability (chaos) of the chain and can be considered as a probe for
the dynamics of the system to which the impurity is coupled. Finally, it is
shown that the Kolmogorov-Sinai entropy of the chain has a crossover from weak
to strong chaos at the same value of the energy density that the crossover
value of largest Lyapunov exponent. Implications of this result
are discussed.Comment: 6 pages, 5 figures, revtex4 styl
Geometric dynamical observables in rare gas crystals
We present a detailed description of how a differential geometric approach to
Hamiltonian dynamics can be used for determining the existence of a crossover
between different dynamical regimes in a realistic system, a model of a rare
gas solid. Such a geometric approach allows to locate the energy threshold
between weakly and strongly chaotic regimes, and to estimate the largest
Lyapunov exponent. We show how standard mehods of classical statistical
mechanics, i.e. Monte Carlo simulations, can be used for our computational
purposes. Finally we consider a Lennard Jones crystal modeling solid Xenon. The
value of the energy threshold turns out to be in excellent agreement with the
numerical estimate based on the crossover between slow and fast relaxation to
equilibrium obtained in a previous work by molecular dynamics simulations.Comment: RevTeX, 19 pages, 6 PostScript figures, submitted to Phys. Rev.
Phase transitions as topology changes in configuration space: an exact result
The phase transition in the mean-field XY model is shown analytically to be
related to a topological change in its configuration space. Such a topology
change is completely described by means of Morse theory allowing a computation
of the Euler characteristic--of suitable submanifolds of configuration
space--which shows a sharp discontinuity at the phase transition point, also at
finite N. The present analytic result provides, with previous work, a new key
to a possible connection of topological changes in configuration space as the
origin of phase transitions in a variety of systems.Comment: REVTeX file, 5 pages, 1 PostScript figur
Geometric Approach to Lyapunov Analysis in Hamiltonian Dynamics
As is widely recognized in Lyapunov analysis, linearized Hamilton's equations
of motion have two marginal directions for which the Lyapunov exponents vanish.
Those directions are the tangent one to a Hamiltonian flow and the gradient one
of the Hamiltonian function. To separate out these two directions and to apply
Lyapunov analysis effectively in directions for which Lyapunov exponents are
not trivial, a geometric method is proposed for natural Hamiltonian systems, in
particular. In this geometric method, Hamiltonian flows of a natural
Hamiltonian system are regarded as geodesic flows on the cotangent bundle of a
Riemannian manifold with a suitable metric. Stability/instability of the
geodesic flows is then analyzed by linearized equations of motion which are
related to the Jacobi equations on the Riemannian manifold. On some geometric
setting on the cotangent bundle, it is shown that along a geodesic flow in
question, there exist Lyapunov vectors such that two of them are in the two
marginal directions and the others orthogonal to the marginal directions. It is
also pointed out that Lyapunov vectors with such properties can not be obtained
in general by the usual method which uses linearized Hamilton's equations of
motion. Furthermore, it is observed from numerical calculation for a model
system that Lyapunov exponents calculated in both methods, geometric and usual,
coincide with each other, independently of the choice of the methods.Comment: 22 pages, 14 figures, REVTeX
Geometric approach to chaos in the classical dynamics of abelian lattice gauge theory
A Riemannian geometrization of dynamics is used to study chaoticity in the
classical Hamiltonian dynamics of a U(1) lattice gauge theory. This approach
allows one to obtain analytical estimates of the largest Lyapunov exponent in
terms of time averages of geometric quantities. These estimates are compared
with the results of numerical simulations, and turn out to be very close to the
values extrapolated for very large lattice sizes even when the geometric
quantities are computed using small lattices. The scaling of the Lyapunov
exponent with the energy density is found to be well described by a quadratic
power law.Comment: REVTeX, 9 pages, 4 PostScript figures include
Geometry of dynamics, Lyapunov exponents and phase transitions
The Hamiltonian dynamics of classical planar Heisenberg model is numerically
investigated in two and three dimensions. By considering the dynamics as a
geodesic flow on a suitable Riemannian manifold, it is possible to analytically
estimate the largest Lyapunov exponent in terms of some curvature fluctuations.
The agreement between numerical and analytical values for Lyapunov exponents is
very good in a wide range of temperatures. Moreover, in the three dimensional
case, in correspondence with the second order phase transition, the curvature
fluctuations exibit a singular behaviour which is reproduced in an abstract
geometric model suggesting that the phase transition might correspond to a
change in the topology of the manifold whose geodesics are the motions of the
system.Comment: REVTeX, 10 pages, 5 PostScript figures, published versio
Topological conditions for discrete symmetry breaking and phase transitions
In the framework of a recently proposed topological approach to phase
transitions, some sufficient conditions ensuring the presence of the
spontaneous breaking of a Z_2 symmetry and of a symmetry-breaking phase
transition are introduced and discussed. A very simple model, which we refer to
as the hypercubic model, is introduced and solved. The main purpose of this
model is that of illustrating the content of the sufficient conditions, but it
is interesting also in itself due to its simplicity. Then some mean-field
models already known in the literature are discussed in the light of the
sufficient conditions introduced here
Kinematics in Kapteyn's Selected Area 76: Orbital Motions Within the Highly Substructured Anticenter Stream
We have measured the mean three-dimensional kinematics of stars in Kapteyn's
Selected Area (SA) 76 (l=209.3, b=26.4 degrees) that were selected to be
Anticenter Stream (ACS) members on the basis of their radial velocities, proper
motions, and location in the color-magnitude diagram. From a total of 31 stars
ascertained to be ACS members primarily from its main sequence turnoff, a mean
ACS radial velocity (derived from spectra obtained with the Hydra multi-object
spectrograph on the WIYN 3.5m telescope) of V_helio = 97.0 +/- 2.8 km/s was
determined, with an intrinsic velocity dispersion sigma_0 = 12.8 \pm 2.1 km/s.
The mean absolute proper motions of these 31 ACS members are mu_alpha cos
(delta) = -1.20 +/- 0.34 mas/yr and mu_delta = -0.78 \pm 0.36 mas/yr. At a
distance to the ACS of 10 \pm 3 kpc, these measured kinematical quantities
produce an orbit that deviates by ~30 degrees from the well-defined swath of
stellar overdensity constituting the Anticenter Stream in the western portion
of the Sloan Digital Sky Survey footprint. We explore possible explanations for
this, and suggest that our data in SA 76 are measuring the motion of a
kinematically cold sub-stream among the ACS debris that was likely a fragment
of the same infalling structure that created the larger ACS system. The ACS is
clearly separated spatially from the majority of claimed Monoceros ring
detections in this region of the sky; however, with the data in hand, we are
unable to either confirm or rule out an association between the ACS and the
poorly-understood Monoceros structure.Comment: Accepted to ApJ. 48 pages, 20 figures, preprint forma
- …