Quantization of generic chaotic 3D billiard with smooth boundary II: structure of high-lying eigenstates

This is the first survey of highly excited eigenstates of a chaotic 3D billiard. We introduce a strongly chaotic 3D billiard with a smooth boundary and we manage to calculate accurate eigenstates with sequential number (of a 48-fold desymmetrized billiard) about 45,000. Besides the brute-force calculation of 3D wavefunctions we propose and illustrate another two representations of eigenstates of quantum 3D billiards: (i) normal derivative of a wavefunction over the boundary surface, and (ii) ray - angular momentum representation. The majority of eigenstates is found to be more or less uniformly extended over the entire energy surface, as expected, but there is also a fraction of strongly localized - scarred eigenstates which are localized either (i) on to classical periodic orbits or (ii) on to planes which carry (2+2)-dim classically invariant manifolds, although the classical dynamics is strongly chaotic and non-diffusive.Comment: 12 pages in plain Latex (5 figures in PCL format available upon request) Submitted to Phys.Lett.

Super-diffusion in one-dimensional quantum lattice models

We identify a class of one-dimensional spin and fermionic lattice models which display diverging spin and charge diffusion constants, including several paradigmatic models of exactly solvable strongly correlated many-body dynamics such as the isotropic Heisenberg spin chains, the Fermi-Hubbard model, and the t-J model at the integrable point. Using the hydrodynamic transport theory, we derive an analytic lower bound on the spin and charge diffusion constants by calculating the curvature of the corresponding Drude weights at half filling, and demonstrate that for certain lattice models with isotropic interactions some of the Noether charges exhibit super-diffusive transport at finite temperature and half filling.Comment: 4 pages + appendices, v2 as publishe

Dynamical properties of a particle in a wave packet: scaling invariance and boundary crisis

Some dynamical properties present in a problem concerning the acceleration of particles in a wave packet are studied. The dynamics of the model is described in terms of a two-dimensional area preserving map. We show that the phase space is mixed in the sense that there are regular and chaotic regions coexisting. We use a connection with the standard map in order to find the position of the first invariant spanning curve which borders the chaotic sea. We find that the position of the first invariant spanning curve increases as a power of the control parameter with the exponent 2/3. The standard deviation of the kinetic energy of an ensemble of initial conditions obeys a power law as a function of time, and saturates after some crossover. Scaling formalism is used in order to characterize the chaotic region close to the transition from integrability to nonintegrability and a relationship between the power law exponents is derived. The formalism can be applied in many different systems with mixed phase space. Then, dissipation is introduced into the model and therefore the property of area preservation is broken, and consequently attractors are observed. We show that after a small change of the dissipation, the chaotic attractor as well as its basin of attraction are destroyed, thus leading the system to experience a boundary crisis. The transient after the crisis follows a power law with exponent -2.Comment: Chaos, Solitons & Fractals, 201

Engineering fidelity echoes in Bose-Hubbard Hamiltonians

We analyze the fidelity decay for a system of interacting bosons described by a Bose-Hubbard Hamiltonian. We find echoes associated with "non-universal" structures that dominate the energy landscape of the perturbation operator. Despite their classical origin, these echoes persist deep into the quantum (perturbative) regime and can be described by an improved random matrix modeling. In the opposite limit of strong perturbations (and high enough energies), classical considerations reveal the importance of self-trapping phenomena in the echo efficiency.Comment: 6 pages, use epl2.cls class, 5 figures Cross reference with nlin, quant-phy

Entanglement in a periodic quench

We consider a chain of free electrons with periodically switched dimerization and study the entanglement entropy of a segment with the remainder of the system. We show that it evolves in a stepwise manner towards a value proportional to the length of the segment and displays in general slow oscillations. For particular quench periods and full dimerization an explicit solution is given. Relations to equilibrium lattice models are pointed out.Comment: 19 pages, 12 figures, 5 references adde

Chaos and Complexity of quantum motion

The problem of characterizing complexity of quantum dynamics - in particular of locally interacting chains of quantum particles - will be reviewed and discussed from several different perspectives: (i) stability of motion against external perturbations and decoherence, (ii) efficiency of quantum simulation in terms of classical computation and entanglement production in operator spaces, (iii) quantum transport, relaxation to equilibrium and quantum mixing, and (iv) computation of quantum dynamical entropies. Discussions of all these criteria will be confronted with the established criteria of integrability or quantum chaos, and sometimes quite surprising conclusions are found. Some conjectures and interesting open problems in ergodic theory of the quantum many problem are suggested.Comment: 45 pages, 22 figures, final version, at press in J. Phys. A, special issue on Quantum Informatio

Is efficiency of classical simulations of quantum dynamics related to integrability?

Efficiency of time-evolution of quantum observables, and thermal states of quenched hamiltonians, is studied using time-dependent density matrix renormalization group method in a family of generic quantum spin chains which undergo a transition from integrable to non-integrable - quantum chaotic case as control parameters are varied. Quantum states (observables) are represented in terms of matrix-product-operators with rank D_\epsilon(t), such that evolution of a long chain is accurate within fidelity error \epsilon up to time t. We find that rank generally increases exponentially, D_\epsilon(t) \propto \exp(const t), unless the system is integrable in which case we find polynomial increase.Comment: 4 pages; v2. added paragraph discussing pure state

Can chaotic quantum energy levels statistics be characterized using information geometry and inference methods?

In this paper, we review our novel information geometrodynamical approach to chaos (IGAC) on curved statistical manifolds and we emphasize the usefulness of our information-geometrodynamical entropy (IGE) as an indicator of chaoticity in a simple application. Furthermore, knowing that integrable and chaotic quantum antiferromagnetic Ising chains are characterized by asymptotic logarithmic and linear growths of their operator space entanglement entropies, respectively, we apply our IGAC to present an alternative characterization of such systems. Remarkably, we show that in the former case the IGE exhibits asymptotic logarithmic growth while in the latter case the IGE exhibits asymptotic linear growth. At this stage of its development, IGAC remains an ambitious unifying information-geometric theoretical construct for the study of chaotic dynamics with several unsolved problems. However, based on our recent findings, we believe it could provide an interesting, innovative and potentially powerful way to study and understand the very important and challenging problems of classical and quantum chaos.Comment: 21 page

Estimating the number of neurons in multi-neuronal spike trains

A common way of studying the relationship between neural activity and behavior is through the analysis of neuronal spike trains that are recorded using one or more electrodes implanted in the brain. Each spike train typically contains spikes generated by multiple neurons. A natural question that arises is "what is the number of neurons $\nu$ generating the spike train?"; This article proposes a method-of-moments technique for estimating $\nu$. This technique estimates the noise nonparametrically using data from the silent region of the spike train and it applies to isolated spikes with a possibly small, but nonnegligible, presence of overlapping spikes. Conditions are established in which the resulting estimator for $\nu$ is shown to be strongly consistent. To gauge its finite sample performance, the technique is applied to simulated spike trains as well as to actual neuronal spike train data.Comment: Published in at http://dx.doi.org/10.1214/10-AOAS371 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org

Dual symplectic classical circuits: An exactly solvable model of many-body chaos

We propose a general exact method of calculating dynamical correlation functions in dual symplectic brick-wall circuits in one dimension. These are deterministic classical many-body dynamical systems which can be interpreted in terms of symplectic dynamics in two orthogonal (time and space) directions. In close analogy with quantum dual-unitary circuits, we prove that two-point dynamical correlation functions are non-vanishing only along the edges of the light cones. The dynamical correlations are exactly computable in terms of a one-site Markov transfer operator, which is generally of infinite dimensionality. We test our theory in a specific family of dual-symplectic circuits, describing the dynamics of a classical Floquet spin chain. Remarkably, for these models, the rotational symmetry leads to a transfer operator with a block diagonal form on the basis of spherical harmonics. This allows us to obtain analytical predictions for simple local observables. We demonstrate the validity of our theory by comparison with Montecarlo simulations, displaying excellent agreement for different choices of observables.Comment: 16 pages, 5 figure