Computing Lyapunov spectra with continuous Gram-Schmidt orthonormalization

We present a straightforward and reliable continuous method for computing the full or a partial Lyapunov spectrum associated with a dynamical system specified by a set of differential equations. We do this by introducing a stability parameter beta>0 and augmenting the dynamical system with an orthonormal k-dimensional frame and a Lyapunov vector such that the frame is continuously Gram-Schmidt orthonormalized and at most linear growth of the dynamical variables is involved. We prove that the method is strongly stable when beta > -lambda_k where lambda_k is the k'th Lyapunov exponent in descending order and we show through examples how the method is implemented. It extends many previous results.Comment: 14 pages, 10 PS figures, ioplppt.sty, iopl12.sty, epsfig.sty 44 k

Phase transition in a class of non-linear random networks

We discuss the complex dynamics of a non-linear random networks model, as a function of the connectivity k between the elements of the network. We show that this class of networks exhibit an order-chaos phase transition for a critical connectivity k = 2. Also, we show that both, pairwise correlation and complexity measures are maximized in dynamically critical networks. These results are in good agreement with the previously reported studies on random Boolean networks and random threshold networks, and show once again that critical networks provide an optimal coordination of diverse behavior.Comment: 9 pages, 3 figures, revised versio

Amplitude death in coupled chaotic oscillators

Amplitude death can occur in chaotic dynamical systems with time-delay coupling, similar to the case of coupled limit cycles. The coupling leads to stabilization of fixed points of the subsystems. This phenomenon is quite general, and occurs for identical as well as nonidentical coupled chaotic systems. Using the Lorenz and R\"ossler chaotic oscillators to construct representative systems, various possible transitions from chaotic dynamics to fixed points are discussed.Comment: To be published in PR

Refining Finite-Time Lyapunov Exponent Ridges and the Challenges of Classifying Them

While more rigorous and sophisticated methods for identifying Lagrangian based coherent structures exist, the finite-time Lyapunov exponent (FTLE) field remains a straightforward and popular method for gaining some insight into transport by complex, time-dependent two-dimensional flows. In light of its enduring appeal, and in support of good practice, we begin by investigating the effects of discretization and noise on two numerical approaches for calculating the FTLE field. A practical method to extract and refine FTLE ridges in two-dimensional flows, which builds on previous methods, is then presented. Seeking to better ascertain the role of a FTLE ridge in flow transport, we adapt an existing classification scheme and provide a thorough treatment of the challenges of classifying the types of deformation represented by a FTLE ridge. As a practical demonstration, the methods are applied to an ocean surface velocity field data set generated by a numerical model. (C) 2015 AIP Publishing LLC.ONR N000141210665Center for Nonlinear Dynamic

Time-reversed symmetry and covariant Lyapunov vectors for simple particle models in and out of thermal equilibrium

Recently, a new algorithm for the computation of covariant Lyapunov vectors and of corresponding local Lyapunov exponents has become available. Here we study the properties of these still unfamiliar quantities for a number of simple models, including an harmonic oscillator coupled to a thermal gradient with a two-stage thermostat, which leaves the system ergodic and fully time reversible. We explicitly demonstrate how time-reversal invariance affects the perturbation vectors in tangent space and the associated local Lyapunov exponents. We also find that the local covariant exponents vary discontinuously along directions transverse to the phase flow.Comment: 13 pages, 11 figures submitted to Physical Review E, 201

Measuring Topological Chaos

The orbits of fluid particles in two dimensions effectively act as topological obstacles to material lines. A spacetime plot of the orbits of such particles can be regarded as a braid whose properties reflect the underlying dynamics. For a chaotic flow, the braid generated by the motion of three or more fluid particles is computed. A ``braiding exponent'' is then defined to characterize the complexity of the braid. This exponent is proportional to the usual Lyapunov exponent of the flow, associated with separation of nearby trajectories. Measuring chaos in this manner has several advantages, especially from the experimental viewpoint, since neither nearby trajectories nor derivatives of the velocity field are needed.Comment: 4 pages, 6 figures. RevTeX 4 with PSFrag macro

Quantum oscillations in graphene in the presence of disorder and interactions

Quantum oscillations in graphene is discussed. The effect of interactions are addressed by Kohn's theorem regarding de Haas-van Alphen oscillations, which states that electron-electron interactions cannot affect the oscillation frequencies as long as disorder is neglected and the system is sufficiently screened, which should be valid for chemical potentials not very close to the Dirac point. We determine the positions of Landau levels in the presence of potential disorder from exact transfer matrix and finite size diagonalization calculations. The positions are shown to be unshifted even for moderate disorder; stronger disorder, can, however, lead to shifts, but this also appears minimal even for disorder width as large as one-half of the bare hopping matrix element on the graphene lattice. Shubnikov-de Haas oscillations of the conductivity are calculated analytically within a self-consistent Born approximation of impurity scattering. The oscillatory part of the conductivity follows the widely invoked Lifshitz-Kosevich form when certain mass and frequency parameters are properly interpreted.Comment: Appendix A was removed, as the content of it is already contained in Ref. 17. Thanks to M. A. H. Vozmedian

Critical conductance of two-dimensional chiral systems with random magnetic flux

The zero temperature transport properties of two-dimensional lattice systems with static random magnetic flux per plaquette and zero mean are investigated numerically. We study the two-terminal conductance and its dependence on energy, sample size, and magnetic flux strength. The influence of boundary conditions and of the oddness of the number of sites in the transverse direction is also studied. We confirm the existence of a critical chiral state in the middle of the energy band and calculate the critical exponent nu=0.35 +/- 0.03 for the divergence of the localization length. The sample averaged scale independent critical conductance _c turns out to be a function of the amplitude of the flux fluctuations whereas the variance of the respective conductance distributions appears to be universal. All electronic states outside of the band center are found to be localized.Comment: to appear in Phys. Rev.

Advection of vector fields by chaotic flows

We have introduced a new transfer operator for chaotic flows whose leading eigenvalue yields the dynamo rate of the fast kinematic dynamo and applied cycle expansion of the Fredholm determinant of the new operator to evaluation of its spectrum. The theory hs been tested on a normal form model of the vector advecting dynamical flow. If the model is a simple map with constant time between two iterations, the dynamo rate is the same as the escape rate of scalar quantties. However, a spread in Poincar\'e section return times lifts the degeneracy of the vector and scalar advection rates, and leads to dynamo rates that dominate over the scalar advection rates. For sufficiently large time spreads we have even found repellers for which the magnetic field grows exponentially, even though the scalar densities are decaying exponentially.Comment: 12 pages, Latex. Ask for figures from [email protected]

The Structure on Invariant Measures of C1C^1 generic diffeomorphisms

Let $\Lambda$ be an isolated non-trival transitive set of a $C^1$ generic diffeomorphism f\in\Diff(M). We show that the space of invariant measures supported on $\Lambda$ coincides with the space of accumulation measures of time averages on one orbit. Moreover, the set of points having this property is residual in $\Lambda$ (which implies the set of irregular$^+$ points is also residual in $\Lambda$). As an application, we show that the non-uniform hyperbolicity of irregular$^+$ points in $\Lambda$ with totally 0 measure (resp., the non-uniform hyperbolicity of a generic subset in $\Lambda$) determines the uniform hyperbolicity of $\Lambda$