Source author record

Arul Lakshminarayan

Arul Lakshminarayan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

30 published item(s)

preprint2023arXiv

Probing Dynamical Sensitivity of a Non-KAM System Through Out-of-Time-Order Correlators

Non-KAM (Kolmogorov-Arnold-Moser) systems, when perturbed by weak time-dependent fields, offer a fast route to classical chaos through an abrupt breaking of invariant phase space tori. In this work, we employ out-of-time-order correlators (OTOCs) to study the dynamical sensitivity of a perturbed non-KAM system in the quantum limit as the parameter that characterizes the $\textit{resonance}$ condition is slowly varied. For this purpose, we consider a quantized kicked harmonic oscillator (KHO) model, which displays stochastic webs resembling Arnold's diffusion that facilitate large-scale diffusion in the phase space. Although the Lyapunov exponent of the KHO at resonances remains close to zero in the weak perturbative regime, making the system weakly chaotic in the conventional sense, the classical phase space undergoes significant structural changes. Motivated by this, we study the OTOCs when the system is in resonance and contrast the results with the non-resonant case. At resonances, we observe that the long-time dynamics of the OTOCs are sensitive to these structural changes, where they grow quadratically as opposed to linear or stagnant growth at non-resonances. On the other hand, our findings suggest that the short-time dynamics remain relatively more stable and show the exponential growth found in the literature for unstable fixed points. The numerical results are backed by analytical expressions derived for a few special cases. We will then extend our findings concerning the non-resonant cases to a broad class of near-integrable KAM systems.

preprint2022arXiv

Construction and local equivalence of dual-unitary operators: from dynamical maps to quantum combinatorial designs

While quantum circuits built from two-particle dual-unitary (maximally entangled) operators serve as minimal models of typically nonintegrable many-body systems, the construction and characterization of dual-unitary operators themselves are only partially understood. A nonlinear map on the space of unitary operators was proposed in PRL.~125, 070501 (2020) that results in operators being arbitrarily close to dual unitaries. Here we study the map analytically for the two-qubit case describing the basins of attraction, fixed points, and rates of approach to dual unitaries. A subset of dual-unitary operators having maximum entangling power are 2-unitary operators or perfect tensors, and are equivalent to four-party absolutely maximally entangled states. It is known that they only exist if the local dimension is larger than $d=2$. We use the nonlinear map, and introduce stochastic variants of it, to construct explicit examples of new dual and 2-unitary operators. A necessary criterion for their local unitary equivalence to distinguish classes is also introduced and used to display various concrete results and a conjecture in $d=3$. It is known that orthogonal Latin squares provide a ``classical combinatorial design" for constructing permutations that are 2-unitary. We extend the underlying design from classical to genuine quantum ones for general dual-unitary operators and give an example of what might be the smallest sized genuinely quantum design of a 2-unitary in $d=4$.

preprint2022arXiv

Out-of-time-order correlators of nonlocal block-spin and random observables in integrable and nonintegrable spin chains

Out-of-time-order correlators (OTOC) in the Ising Floquet system, that can be both integrable and nonintegrable is studied. Instead of localized spin observables, we study contiguous symmetric blocks of spins or random operators localized on these blocks as observables. We find only power-law growth of OTOC in both integrable and nonintegrable regimes. In the non-integrable regime, beyond the scrambling time, there is an exponential saturation of the OTOC to values consistent with random matrix theory. This motivates the use of "pre-scrambled" random block operators as observables. A pure exponential saturation of OTOC in both integrable and nonintegrable system is observed, without a scrambling phase. Averaging over random observables from the Gaussian unitary ensemble, the OTOC is found to be exactly same as the operator entanglement entropy, whose exponential saturation has been observed in previous studies of such spin-chains.

preprint2021arXiv

From dual-unitary to quantum Bernoulli circuits: Role of the entangling power in constructing a quantum ergodic hierarchy

Deterministic classical dynamical systems have an ergodic hierarchy, from ergodic through mixing, to Bernoulli systems that are "as random as a coin-toss". Dual-unitary circuits have been recently introduced as solvable models of many-body nonintegrable quantum chaotic systems having a hierarchy of ergodic properties. We extend this to include the apex of a putative quantum ergodic hierarchy which is Bernoulli, in the sense that correlations of single and two-particle observables vanish at space-time separated points. We derive a condition based on the entangling power $e_p(U)$ of the basic two-particle unitary building block, $U$, of the circuit, that guarantees mixing, and when maximized, corresponds to Bernoulli circuits. Additionally we show, both analytically and numerically, how local-averaging over random realizations of the single-particle unitaries, $u_i$ and $v_i$ such that the building block is $U^\prime = (u_1 \otimes u_2 ) U (v_1 \otimes v_2 )$ leads to an identification of the average mixing rate as being determined predominantly by the entangling power $e_p(U)$. Finally we provide several, both analytical and numerical, ways to construct dual-unitary operators covering the entire possible range of entangling power. We construct a coupled quantum cat map which is dual-unitary for all local dimensions and a 2-unitary or perfect tensor for odd local dimensions, and can be used to build Bernoulli circuits.

preprint2021arXiv

Quantum walks with quantum chaotic coins: Of the Loschmidt echo, classical limit and thermalization

Coined discrete-time quantum walks are studied using simple deterministic dynamical systems as coins whose classical limit can range from being integrable to chaotic. It is shown that a Loschmidt echo like fidelity plays a central role and when the coin is chaotic this is approximately the characteristic function of a classical random walker. Thus the classical binomial distribution arises as a limit of the quantum walk and the walker exhibits diffusive growth before eventually becoming ballistic. The coin-walker entanglement growth is shown to be logarithmic in time as in the case of many-body localization and coupled kicked rotors, and saturates to a value that depends on the relative coin and walker space dimensions. In a coin dominated scenario, the chaos can thermalize the quantum walk to typical random states such that the entanglement saturates at the Haar averaged Page value, unlike in a walker dominated case when atypical states seem to be produced.

preprint2021arXiv

Thirty-six entangled officers of Euler: Quantum solution to a classically impossible problem

The negative solution to the famous problem of $36$ officers of Euler implies that there are no two orthogonal Latin squares of order six. We show that the problem has a solution, provided the officers are entangled, and construct orthogonal quantum Latin squares of this size. As a consequence, we find an example of the long-elusive Absolutely Maximally Entangled state AME$(4,6)$ of four subsystems with six levels each, equivalently a $2$-unitary matrix of size $36$, which maximizes the entangling power among all bipartite unitary gates of this dimension, or a perfect tensor with four indices, each running from one to six. This special state deserves the appellation golden AME state as the golden ratio appears prominently in its elements. This result allows us to construct a pure nonadditive quhex quantum error detection code $(\!(3,6,2)\!)_6$, which saturates the Singleton bound and allows one to encode a $6$-level state into a triplet of such states.

preprint2020arXiv

Creating ensembles of dual unitary and maximally entangling quantum evolutions

Maximally entangled bipartite unitary operators or gates find various applications from quantum information to being building blocks of minimal models of many-body quantum chaos, and have been referred to as "dual unitaries". Dual unitary operators that can create the maximum average entanglement when acting on product states have to satisfy additional constraints. These have been called "2-unitaries" and are examples of perfect tensors that can be used to construct absolutely maximally entangled states of four parties. Hitherto, no systematic method exists, in any local dimension, which result in the formation of such special classes of unitary operators. We outline an iterative protocol, a nonlinear map on the space of unitary operators, that creates ensembles whose members are arbitrarily close to being dual unitaries, while for qutrits and ququads we find that a slightly modified protocol yields a plethora of 2-unitaries. We further characterize the dual unitary operators via their entangling power and the 2-unitaries via the distribution of entanglement created from unentangled states.

preprint2020arXiv

Probing the randomness of ergodic states: extreme-value statistics in the ergodic and many-body-localized phases

The extreme-value statistics of the entanglement spectrum in disordered spin chains possessing a many-body localization transition is examined. It is expected that eigenstates in the metallic or ergodic phase, behave as random states and hence the eigenvalues of the reduced density matrix, commonly referred to as the entanglement spectrum, are expected to follow the eigenvalue statistics of a trace normalized Wishart ensemble. In particular, the density of eigenvalues is supposed to follow the universal Marchenko-Pastur distribution. We find deviations in the tails both for the disordered XXZ with total $S_z$ conserved in the half-filled sector as well as in a model that breaks this conservation. A sensitive measure of deviations is provided by the largest eigenvalue, which in the case of the Wishart ensemble after appropriate shift and scaling follows the universal Tracy-Widom distribution. We show that for the models considered, in the metallic phase, the largest eigenvalue of the reduced density matrix of eigenvector, instead follows the generalized extreme-value statistics bordering on the Fisher-Tipett-Gumbel distribution indicating that the correlations between eigenvalues are much weaker compared to the Wishart ensemble. We show by means of distributions conditional on the high entropy and normalized participation ratio of eigenstates that the conditional entanglement spectrum still follows generalized extreme value distribution. In the deeply localized phase we find heavy tailed distributions and Lévy stable laws in an appropriately scaled function of the largest and second largest eigenvalues. The scaling is motivated by a recently developed perturbation theory of weakly coupled chaotic systems.

preprint2019arXiv

Entanglement production by interaction quenches of quantum chaotic subsystems

The entanglement production in bipartite quantum systems is studied for initially unentangled product eigenstates of the subsystems, which are assumed to be quantum chaotic. Based on a perturbative computation of the Schmidt eigenvalues of the reduced density matrix, explicit expressions for the time-dependence of entanglement entropies, including the von Neumann entropy, are given. An appropriate re-scaling of time and the entropies by their saturation values leads a universal curve, independent of the interaction. The extension to the non-perturbative regime is performed using a recursively embedded perturbation theory to produce the full transition and the saturation values. The analytical results are found to be in good agreement with numerical results for random matrix computations and a dynamical system given by a pair of coupled kicked rotors.

preprint2019arXiv

Scrambling in strongly chaotic weakly coupled bipartite systems: Universality beyond the Ehrenfest time-scale

Out-of-time-order correlators (OTOC), vigorously being explored as a measure of quantum chaos and information scrambling, is studied here in the natural and simplest multi-particle context of bipartite systems. We show that two strongly chaotic and weakly interacting subsystems display two distinct phases in the growth of OTOC.The first is dominated by intra-subsystem scrambling, when an exponential growth with a positive Lyapunov exponent is observed till the Ehrenfest time. This phase is essentially independent of the interaction, while the second phase is an interaction dominated exponential approach to saturation that is universal and described by a random matrix model. This simple random matrix model of weakly interacting strongly chaotic bipartite systems, previously employed for studying entanglement and spectral transitions, is approximately analytically solvable for its OTOC. The example of two coupled kicked rotors is used to demonstrate the two phases, and the extent to which the random matrix model is applicable. That the two phases correspond to delocalization in the subsystems followed by inter-subsystem mixing is seen via the participation ratio in phase-space. We also point out that the second, universal, phase alone exists when the observables are in a sense already scrambled. Thus, while the post-Ehrenfest time OTOC growth is in general not well-understood, the case of strongly chaotic and weakly coupled systems presents an, perhaps important, exception.

preprint2016arXiv

Entanglement and localization transitions in eigenstates of interacting chaotic systems

The entanglement and localization in eigenstates of strongly chaotic subsystems are studied as a function of their interaction strength. Excellent measures for this purpose are the von-Neumann entropy, Havrda-Charv{\' a}t-Tsallis entropies, and the averaged inverse participation ratio. All the entropies are shown to follow a remarkably simple exponential form, which describes a universal and rapid transition to nearly maximal entanglement for increasing interaction strength. An unexpectedly exact relationship between the subsystem averaged inverse participation ratio and purity is derived that infers the transition in the localization as well.

preprint2015arXiv

Local entanglement structure across a many-body localization transition

Local entanglement between pairs of spins, as measured by concurrence, is investigated in a disordered spin model that displays a transition from an ergodic to a many-body localized phase in excited states. It is shown that the concurrence vanishes in the ergodic phase and becomes nonzero and increases in the many-body localized phase. This happen to be correlated with the transition in the spectral statistics from Wigner to Poissonian distribution. A scaling form is found to exist in the second derivative of the concurrence with the disorder strength. It also displays a critical value for the localization transition that is close to what is known in the literature from other measures. An exponential decay of concurrence with distance between spins is observed in the localized phase. Nearest neighbor spin concurrence in this phase is also found to be strongly correlated with the disorder configuration of onsite fields: nearly similar fields implying larger entanglement.

preprint2015arXiv

Real eigenvalues of non-Gaussian random matrices and their products

We study the properties of the eigenvalues of real random matrices and their products. It is known that when the matrix elements are Gaussian-distributed independent random variables, the fraction of real eigenvalues tends to unity as the number of matrices in the product increases. Here we present numerical evidence that this phenomenon is robust with respect to the probability distribution of matrix elements, and is therefore a general property that merits detailed investigation. Since the elements of the product matrix are no longer distributed as those of the single matrix nor they remain independent random variables, we study the role of these two factors in detail. We study numerically the properties of the Hadamard (or Schur) product of matrices and also the product of matrices whose entries are independent but have the same marginal distribution as that of normal products of matrices, and find that under repeated multiplication, the probability of all eigenvalues to be real increases in both cases, but saturates to a constant below unity showing that the correlations amongst the matrix elements are responsible for the approach to one. To investigate the role of the non-normal nature of the probability distributions, we present a thorough analytical treatment of the $2 \times 2$ single matrix for several standard distributions. Within the class of smooth distributions with zero mean and finite variance, our results indicate that the Gaussian distribution has the maximum probability of real eigenvalues, but the Cauchy distribution characterised by infinite variance is found to have a larger probability of real eigenvalues than the normal. We also find that for the two-dimensional single matrices, the probability of real eigenvalues lies in the range [5/8,7/8].

preprint2015arXiv

Records in the classical and quantum standard map

Record statistics is the study of how new highs or lows are created and sustained in any dynamical process. The study of the highest or lowest records constitute the study of extreme values. This paper represents an exploration of record statistics for certain aspects of the classical and quantum standard map. For instance the momentum square or energy records is shown to behave like that of records in random walks when the classical standard map is in a regime of hard chaos. However different power laws is observed for the mixed phase space regimes. The presence of accelerator modes are well-known to create anomalous diffusion and we notice here that the record statistics is very sensitive to their presence. We also discuss records in random vectors and use it to analyze the {\it quantum} standard map via records in their eigenfunction intensities, reviewing some recent results along the way.

preprint2015arXiv

Universal scaling of spectral fluctuation transitions for interacting chaotic systems

The spectral properties of interacting strongly chaotic systems are investigated for growing interaction strength. A very sensitive transition from Poisson statistics to that of random matrix theory is found. We introduce a new random matrix ensemble modeling this dynamical symmetry breaking transition which turns out to be universal and depends on a single scaling parameter only. Coupled kicked rotors, a dynamical systems paradigm for such transitions, are compared with this ensemble and excellent agreement is found for the nearest-neighbor-spacing distribution. It turns out that this transition is described quite accurately using perturbation theory.

preprint2014arXiv

Diagonal unitary entangling gates and contradiagonal quantum states

Nonlocal properties of an ensemble of diagonal random unitary matrices of order $N^2$ are investigated. The average Schmidt strength of such a bipartite diagonal quantum gate is shown to scale as $\log N$, in contrast to the $\log N^2$ behavior characteristic to random unitary gates. Entangling power of a diagonal gate $U$ is related to the von Neumann entropy of an auxiliary quantum state $ρ=AA^{\dagger}/N^2$, where the square matrix $A$ is obtained by reshaping the vector of diagonal elements of $U$ of length $N^2$ into a square matrix of order $N$. This fact provides a motivation to study the ensemble of non-hermitian unimodular matrices $A$, with all entries of the same modulus and random phases and the ensemble of quantum states $ρ$, such that all their diagonal entries are equal to $1/N$. Such a state is contradiagonal with respect to the computational basis, in sense that among all unitary equivalent states it maximizes the entropy copied to the environment due to the coarse graining process. The first four moments of the squared singular values of the unimodular ensemble are derived, based on which we conjecture a connection to a recently studied combinatorial object called the "Borel triangle". This allows us to find exactly the mean von Neumann entropy for random phase density matrices and the average entanglement for the corresponding ensemble of bipartite pure states.

preprint2013arXiv

On the number of real eigenvalues of products of random matrices and an application to quantum entanglement

The probability that there are $k$ real eigenvalues for an $n$ dimensional real random matrix is known. Here we study this for the case of products of independent random matrices. Relating the problem of the probability that the product of two real 2 dimensional random matrices has real eigenvalues to an issue of optimal quantum entanglement, this is fully analytically solved. It is shown that in $π/4$ fraction of such products the eigenvalues are real. Being greater than the corresponding known probability ($1/\sqrt{2}$) for a single matrix, it is shown numerically that the probability that {\it all} eigenvalues of a product of random matrices are real tends to unity as the number of matrices in the product increases indefinitely. Some other numerical explorations, including the expected number of real eigenvalues is also presented, where an exponential approach of the expected number to the dimension of the matrix seems to hold.

preprint2013arXiv

Quenching and generation of random states in a kicked Ising model

The kicked Ising model with both a pulsed transverse and a continuous longitudinal field is studied numerically. Starting from a large transverse field and a state that is nearly an eigenstate, the pulsed transverse field is quenched with a simultaneous enhancement of the longitudinal field. The generation of multipartite entanglement is observed along with a phenomenon akin to quantum resonance when the entanglement does not evolve for certain values of the pulse duration. Away from the resonance, the longitudinal field can drive the entanglement to near maximum values that is shown to agree well with those of random states. Further evidence is presented that the time evolved states obtained do have some statistical properties of such random states. For contrast the case when the fields have a steady value is also discussed.

preprint2013arXiv

Record statistics in random vectors and quantum chaos

The record statistics of complex random states are analytically calculated, and shown that the probability of a record intensity is a Bernoulli process. The correlation due to normalization leads to a probability distribution of the records that is non-universal but tends to the Gumbel distribution asymptotically. The quantum standard map is used to study these statistics for the effect of correlations apart from normalization. It is seen that in the mixed phase space regime the number of intensity records is a power law in the dimensionality of the state as opposed to the logarithmic growth for random states.

preprint2013arXiv

Using Partial Transpose and Realignment to generate Local Unitary Invariants

Motivated by link transformations of lattice gauge theory, a method for generating local unitary invariants, especially for a system of qubits, has been pointed out in an earlier work [M. S. Williamson {\it et. al.}, Phys. Rev. A {\bf 83}, 062308 (2011)]. This paper first points the equivalence of the so constructed transformations to the combined operations of partial transpose and realignment. This allows construction of local unitary invariants of any system, with subsystems of arbitrary dimensions. Some properties of the resulting operators and consequences for pure tripartite higher dimensional states are briefly discussed.

preprint2012arXiv

Entanglement signatures for the dimerization transition in the Majumdar-Ghosh model

The transition from a gapless liquid to a gapped dimerized ground state that occurs in the frustrated antiferromagnetic Majumdar-Ghosh (or $J_1-J_2$ Heisenberg) model is revisited from the point of view of entanglement. We study the evolution of entanglement spectra, a "projected subspace" block entropy, and concurrence in the Schmidt vectors through the transition. The standard tool of Schmidt decomposition along with the existence of the unique MG point where the ground states are degenerate and known exactly, suggests the projection into two orthogonal subspaces that is useful even away from this point. Of these, one is a dominant five dimensional subspace containing the complete state at the MG point and the other contributes marginally, albeit with increasing weight as the number of spins is increased. We find that the marginally contributing subspace has a minimum von Neumann entropy in the vicinity of the dimerization transition. Entanglement content between pairs of spins in the Schmidt vectors, studied via concurrence, shows that those belonging to the dominant five dimensional subspace display a clear progress towards dimerization, with the concurrence vanishing on odd/even sublattices, again in the vicinity of the dimerization, and maximizing in the even/odd sublattices at the MG point. In contrast, study of the Schmidt vectors in the marginally contributing subspace, as well as in the projection of the ground state in this space, display pair concurrence which decrease on both the sublattices as the MG point is approached. The robustness of these observations indicate their possible usefulness in the study of models that have similar transitions, and have hitherto been difficult to study using standard entanglement signatures.

preprint2011arXiv

Entanglement optimizing mixtures of two-qubit states

Entanglement in incoherent mixtures of pure states of two qubits is considered via the concurrence measure. A set of pure states is optimal if the concurrence for any mixture of them is the weighted sum of the concurrences of the generating states. When two or three pure real states are mixed it is shown that 28.5% and 5.12% of the cases respectively, are optimal. Conditions that are obeyed by the pure states generating such optimally entangled mixtures are derived. For four or more pure states it is shown that there are no such sets of real states. The implications of these on superposition of two or more dimerized states is discussed. A corollary of these results also show in how many cases rebit concurrence can be the same as that of qubit concurrence.

preprint2011arXiv

Entanglement transitions in random definite particle states

Entanglement within qubits are studied for the subspace of definite particle states or definite number of up spins. A transition from an algebraic decay of entanglement within two qubits with the total number $N$ of qubits, to an exponential one when the number of particles is increased from two to three is studied in detail. In particular the probability that the concurrence is non-zero is calculated using statistical methods and shown to agree with numerical simulations. Further entanglement within a block of $m$ qubits is studied using the log-negativity measure which indicates that a transition from algebraic to exponential decay occurs when the number of particles exceeds $m$. Several algebraic exponents for the decay of the log-negativity are analytically calculated. The transition is shown to be possibly connected with the changes in the density of states of the reduced density matrix, which has a divergence at the zero eigenvalue when the entanglement decays algebraically.

preprint2011arXiv

Localised zero-energy modes in the Kitaev model with vacancy-disorder

We study the effects of vacancy disorder on the Kitaev model defined on a hexagonal lattice. We show that the vacancy disorder induces a zero-mode that is localized at the defect site. We derive analytical forms for these localized wave functions in both the gapped and gapless phases of the Kitaev model. We conjecture that the vacancy disorder can be utilized as a probe of the quantum phase transition (from the gapped to gapless phases) in this model. The behavior of the Inverse Participation Ratio (IPR) in the gapless phase and across the transition is also studied numerically. Comments are made about the behavior of site-site entanglement in the single particle states for the case of a single vacancy.

preprint2011arXiv

On the Kolmogorov-Sinai entropy of many-body Hamiltonian systems

The Kolmogorov-Sinai (K-S) entropy is a central measure of complexity and chaos. Its calculation for many-body systems is an interesting and important challenge. In this paper, the evaluation is formulated by considering $N$-dimensional symplectic maps and deriving a transfer matrix formalism for the stability problem. This approach makes explicit a duality relation that is exactly analogous to one found in a generalized Anderson tight-binding model, and leads to a formally exact expression for the finite-time K-S entropy. Within this formalism there is a hierarchy of approximations, the final one being a diagonal approximation that only makes use of instantaneous Hessians of the potential to find the K-S entropy. By way of a non-trivial illustration, the K-S entropy of $N$ identically coupled kicked rotors (standard maps) is investigated. The validity of the various approximations with kicking strength, particle number, and time are elucidated. An analytic formula for the K-S entropy within the diagonal approximation is derived and its range of validity is also explored.

preprint2010arXiv

Fluctuations in classical sum rules

Classical sum rules arise in a wide variety of physical contexts. Asymptotic expressions have been derived for many of these sum rules in the limit of long orbital period (or large action). Although sum rule convergence may well be exponentially rapid for chaotic systems in a global sense with time, individual contributions to the sums may fluctuate with a width which diverges in time. Our interest is in the global convergence of sum rules as well as their local fluctuations. It turns out that a simple version of a lazy baker map gives an ideal system in which classical sum rules, their corrections, and their fluctuations can be worked out analytically. This is worked out in detail for the Hannay-Ozorio sum rule. In this particular case the rate of convergence of the sum rule is found to be governed by the Pollicott-Ruelle resonances, and both local and global boundaries for which the sum rule may converge are given. In addition, the width of the fluctuations is considered and worked out analytically, and it is shown to have an interesting dependence on the location of the region over which the sum rule is applied. It is also found that as the region of application is decreased in size the fluctuations grow. This suggests a way of controlling the length scale of the fluctuations by considering a time dependent phase space volume, which for the lazy baker map decreases exponentially rapidly with time.

preprint2007arXiv

Extreme statistics of complex random and quantum chaotic states

An exact analytical description of extreme intensity statistics in complex random states is derived. These states have the statistical properties of the Gaussian and Circular Unitary Ensemble eigenstates of random matrix theory. Although the components are correlated by the normalization constraint, it is still possible to derive compact formulae for all values of the dimensionality N. The maximum intensity result slowly approaches the Gumbel distribution even though the variables are bounded, whereas the minimum intensity result rapidly approaches the Weibull distribution. Since random matrix theory is conjectured to be applicable to chaotic quantum systems, we calculate the extreme eigenfunction statistics for the standard map with parameters at which its classical map is fully chaotic. The statistical behaviors are consistent with the finite-N formulae.

preprint2007arXiv

Fluctuations of finite-time stability exponents in the standard map and the detection of small islands

Some statistical properties of finite-time stability exponents in the standard map can be estimated analytically. The mean exponent averaged over the entire phase space behaves quite differently from all the other cumulants. Whereas the mean carries information about the strength of the interaction, and only indirect information about dynamical correlations, the higher cumulants carry information about dynamical correlations and essentially no information about the interaction strength. In particular, the variance and higher cumulants of the exponent are very sensitive to dynamical correlations and easily detect the presence of very small islands of regular motion via their anomalous time-scalings. The average of the stability matrix' inverse trace is even more sensitive to the presence of small islands and has a seemingly fractal behavior in the standard map parameter. The usual accelerator modes and the small islands created through double saddle node bifurcations, which come halfway between the positions in interaction strength of the usual accelerator modes, are clearly visible in the variance, whose time scaling is capable of detecting the presence of islands as small as 0.01% of the phase space. We study these quantities with a local approximation to the trace of the stability matrix which significantly simplifies the numerical calculations as well as allows for generalization of these methods to higher dimensions. We also discuss the nature of this local approximation in some detail.

preprint2007arXiv

The modular multiplication operator and the quantized bakers maps

The modular multiplication operator, a central subroutine in Shor's factoring algorithm, is shown to be a coherent superposition of two quantum bakers maps when the multiplier is 2. The classical limit of the maps being completely chaotic, it is shown that there exist perturbations that push the modular multiplication operator into regimes of generic quantum chaos with spectral fluctuations that are those of random matrices. For the initial state of relevance to Shor's algorithm we study fidelity decay due to phase and bit-flip errors in a single qubit and show exponential decay with shoulders at multiples or half-multiples of the order. A simple model is used to gain some understanding of this behavior.

preprint1995arXiv

On the Quantum Baker's Map and its Unusual Traces

The quantum baker's map is the quantization of a simple classically chaotic system, and has many generic features that have been studied over the last few years. While there exists a semiclassical theory of this map, a more rigorous study of the same revealed some unexpected features which indicated that correction terms of the order of log(h) had to be included in the periodic orbit sum. Such singular semiclassical behaviour was also found in the simplest traces of the quantum map. In this note we study the quantum mechanics of a baker's map which is obtained by reflecting the classical map about its edges, in an effort to understand and circumvent these anomalies. This leads to a real quantum map with traces that follow the usual Gutzwiller-Tabor like semiclassical formulae. We develop the relevant semiclassical periodic orbit sum for this map which is closely related to that of the usual baker's map, with the important difference that the propagators leading to this sum have no anomalous traces.