Source author record

Steven Tomsovic

Steven Tomsovic 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

21works
11topics
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

21 published item(s)

preprint2022arXiv

A corrected Maslov index for complex saddle trajectories

Saddle point approximations, extremely important in a wide variety of physical contexts, require the analytical continuation of canonically conjugate quantities to complex variables in quantum mechanics. An important component of this approximation's implementation is arriving at the phase correction attributable to caustics, which involves determinantal prefactors. The common prescription of using the inverse of half a certain determinant's total accumulated phase sometimes leads to sign errors. The root of this problem is traced to the zeros of the determinants at complex times crossing the real time axis. Deformed complex time contours around the zeros can repair the sign errors that sometimes occur, but a much more practical way is given that links saddles back to associated real trajectories and avoids the necessity of locating the complex time zeros of the determinants.

preprint2022arXiv

Geometric determination of heteroclinic and unstable periodic orbit classical actions

Semiclassical sum rules, such as the Gutzwiller trace formula, depend on the properties of periodic, closed, or homoclinic (heteroclinic) orbits. The interferences embedded in such orbit sums are governed by classical action functions and Maslov indices. For chaotic systems, the relative actions of such orbits can be expressed in terms of phase space areas bounded by segments of stable and unstable manifolds, and Moser invariant curves. This also generates direct relations between periodic orbits and homoclinic (heteroclinic) orbit actions. Simpler, explicit approximate expressions following from the exact relations are given with error estimates. They arise from asymptotic scaling of certain bounded phase space areas. The actions of infinite subsets of periodic orbits are determined by their periods and the locations of the limiting homoclinic points on which they accumulate.

preprint2022arXiv

Semiclassical propagation of coherent states and wave packets: hidden saddles

Semiclassical methods are extremely important in the subjects of wave packet and coherent state dynamics. Unfortunately, these essentially saddle point approximations are considered nearly impossible to carry out in detail for systems with multiple degrees of freedom due to the difficulties of solving the resulting two-point boundary value problems. However, recent developments have extended the applicability to a broader range of systems and circumstances. The most important advances are first to generate a set of real reference trajectories using appropriately reduced dimensional spaces of initial conditions, and second to feed that set into a Newton-Raphson search scheme to locate the $exposed$ complex saddle trajectories. The arguments for this approach were based mostly on intuition and numerical verification. In this paper, the methods are put on a firmer theoretical foundation and then extended to incorporate saddles $hidden$ from Newton-Raphson searches initiated with real trajectories. This hidden class of saddles is relevant to tunneling-type processes, but a hidden saddle can sometimes contribute just as much as or more than an exposed one. The distinctions between hidden and exposed saddles clarifies the interpretation of what constitutes tunneling for wave packets and coherent states in the time domain.

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.

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.

preprint2016arXiv

Generalized Gaussian wave packet dynamics: Integrable and Chaotic Systems

The ultimate semiclassical wave packet propagation technique is a complex, time-dependent WBK method known as generalized Gaussian wave packet dynamics (GGWPD). It requires overcoming many technical difficulties in order to be carried out fully in practice. In its place roughly twenty years ago, linearized wave packet dynamics was generalized to methods that include sets of off-center, real trajectories for both classically integrable and chaotic dynamical systems that completely capture the dynamical transport. The connections between those methods and GGWPD are developed in a way that enables a far more practical implementation of GGWPD. The generally complex saddle point trajectories at its foundation are found using a multi-dimensional, Newton-Raphson root search method that begins with the set of off-center, real trajectories. This is possible because there is a one-to-one correspondence. The neighboring trajectories associated with each off-center, real trajectory form a path that crosses a unique saddle; there are exceptions which are straightforward to identify. The method is applied to the kicked rotor to demonstrate the accuracy improvement as a function of $\hbar$ that comes with using the saddle point trajectories.

preprint2015arXiv

A semiclassical reversibility paradox in simple chaotic systems

Using semiclassical methods, it is possible to construct very accurate approximations in the short wavelength limit of quantum dynamics that rely exclusively on classical dynamical input. For systems whose classical realization is strongly chaotic, there is an exceedingly short logarithmic Ehrenfest time scale, beyond which the quantum and classical dynamics of a system necessarily diverge, and yet the semiclassical construction remains valid far beyond that time. This fact leads to a paradox if one ponders the reversibility and predictability properties of quantum and classical mechanics. They behave very differently relative to each other, with classical dynamics being essentially irreversible/unpredictable, whereas quantum dynamics is reversible/stable. This begs the question, "how can an accurate approximation to a reversible/stable dynamics be constructed from an irreversible/unpredictable one?" The resolution of this incongruity depends on a couple key ingredients, a well-known, inherent, one-way structural stability of chaotic systems and an overlap integral not being amenable to the saddle point method.

preprint2015arXiv

Action differences between fixed points and accurate heteroclinic orbits

A general relation is derived for the action difference between two fixed points and a phase space area bounded by the irreducible component of a heteroclinic tangle. The determination of this area can require accurate calculation of heteroclinic orbits, which are important in a wide range of dynamical system problems. For very strongly chaotic systems initial deviations from a true orbit are magnified by a large exponential rate making direct computational methods fail quickly. Here, a method is developed that avoids direct calculation of the orbit by making use of the well-known stability property of the invariant unstable and stable manifolds. Under an area-preserving map, this property assures that any initial deviation from the stable (unstable) manifold collapses onto themselves under inverse (forward) iterations of the map. Using a set of judiciously chosen auxiliary points on the manifolds, long orbit segments can be calculated using the stable and unstable manifold intersections of the heteroclinic (homoclinic) tangle. Detailed calculations using the example of the kicked rotor are provided along with verification of the relation between action differences. The loop structure of the heteroclinic tangle is necessarily quite different from that of the turnstile for a homoclinic tangle, its analogous partner.

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.

preprint2015arXiv

Wave Transport in One-Dimensional Disordered Systems with Finite-Size Scatterers

We study the problem of wave transport in a one-dimensional disordered system, where the scatterers of the chain are $n$ barriers and wells with statistically independent intensities and with a spatial extension $ł_c$ which may contain an arbitrary number $δ/2π$ of wavelengths, where $δ= k l_c$. We analyze the average Landauer resistance and transmission coefficient of the chain as a function of $n$ and the phase parameter $δ$. For weak scatterers, we find: i) a regime, to be called I, associated with an exponential behavior of the resistance with $n$, ii) a regime, to be called II, for $δ$ in the vicinity of $π$, where the system is almost transparent and less localized, and iii) right in the middle of regime II, for $δ$ very close to $π$, the formation of a band gap, which becomes ever more conspicuous as $n$ increases. In regime II, both the average Landauer resistance and the transmission coefficient show an oscillatory behavior with $n$ and $δ$. These characteristics of the system are found analytically, some of them exactly and some others approximately. The agreement between theory and simulations is excellent, which suggests a strong motivation for the experimental study of these systems. We also present a qualitative discussion of the results.

preprint2012arXiv

Constructing acoustic timefronts using random matrix theory

In a recent letter [Europhys. Lett. 97, 34002 (2012)], random matrix theory is introduced for long-range acoustic propagation in the ocean. The theory is expressed in terms of unitary propagation matrices that represent the scattering between acoustic modes due to sound speed fluctuations induced by the ocean's internal waves. The scattering exhibits a power-law decay as a function of the differences in mode numbers thereby generating a power-law, banded, random unitary matrix ensemble. This work gives a more complete account of that approach and extends the methods to the construction of an ensemble of acoustic timefronts. The result is a very efficient method for studying the statistical properties of timefronts at various propagation ranges that agrees well with propagation based on the parabolic equation. It helps identify which information about the ocean environment survives in the timefronts and how to connect features of the data to the surviving environmental information. It also makes direct connections to methods used in other disordered wave guide contexts where the use of random matrix theory has a multi-decade history.

preprint2012arXiv

Universal Quantum Localizing Transition of a Partial Barrier in a Chaotic Sea

Generic 2D Hamiltonian systems possess partial barriers in their chaotic phase space that restrict classical transport. Quantum mechanically the transport is suppressed if Planck's constant h is large compared to the classical flux, h >> Phi, such that wave packets and states are localized. In contrast, classical transport is mimicked for h << Phi. Designing a quantum map with an isolated partial barrier of controllable flux Phi is the key to investigating the transition from this form of quantum localization to mimicking classical transport. It is observed that quantum transport follows a universal transition curve as a function of the expected scaling parameter Phi/h. We find this curve to be symmetric to Phi/h=1, having a width of two orders of magnitude in Phi/h, and exhibiting no quantized steps. We establish the relevance of local coupling, improving on previous random matrix models relying on global coupling. It turns out that a phenomenological 2x2-model gives an accurate analytical description of the transition curve.

preprint2012arXiv

Wave transport in one-dimensional disordered systems with finite-width potential steps

An amazingly simple model of correlated disorder is a one-dimensional chain of n potential steps with a fixed width lc and random heights. A theoretical analysis of the average transmission coefficient and Landauer resistance as functions of n and klc predicts two distinct regimes of behavior, one marked by extreme sensitivity and the other associated with exponential behavior of the resistance. The sensitivity arises in n and klc for klc approximately pi, where the system is nearly transparent. Numerical simulations match the predictions well, and they suggest a strong motivation for experimental study.

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.

preprint2011arXiv

Random matrix theory for underwater sound propagation

Ocean acoustic propagation can be formulated as a wave guide with a weakly random medium generating multiple scattering. Twenty years ago, this was recognized as a quantum chaos problem, and yet random matrix theory, one pillar of quantum or wave chaos studies, has never been introduced into the subject. The modes of the wave guide provide a representation for the propagation, which in the parabolic approximation is unitary. Scattering induced by the ocean's internal waves leads to a power-law random banded unitary matrix ensemble for long-range deep ocean acoustic propagation. The ensemble has similarities, but differs, from those introduced for studying the Anderson metal-insulator transition. The resulting long-range propagation ensemble statistics agree well with those of full wave propagation using the parabolic equation.

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.

preprint2010arXiv

What is measured in the scanning gate microscopy of a quantum point contact?

The conductance change due to a local perturbation in a phase-coherent nanostructure is calculated. The general expressions to first and second order in the perturbation are applied to the scanning gate microscopy of a two-dimensional electron gas containing a quantum point contact. The first-order correction depends on two scattering states with electrons incoming from opposite leads and is suppressed on a conductance plateau; it is significant in the step regions. On the plateaus, the dominant second-order term likewise depends on scattering states incoming from both sides. It is always negative, exhibits fringes, and has a spatial decay consistent with experiments.

preprint2009arXiv

Investigating quantum transport with an initial value representation of the semiclassical propagator

Quantized systems whose underlying classical dynamics possess an elaborate mixture of regular and chaotic motion can exhibit rather subtle long-time quantum transport phenomena. In a short wavelength regime where semiclassical theories are most relevant, such transport phenomena, being quintessentially interference based, are difficult to understand with the system's specific long-time classical dynamics. Fortunately, semiclassical methods applied to wave packet propagation can provide a natural approach to understanding the connections, even though they are known to break down progressively as time increases. This is due to the fact that some long-time transport properties can be deduced from intermediate-time behavior. Thus, these methods need only retain validity and be carried out on much shorter time scales than the transport phenomena themselves in order to be valuable. The initial value representation of the semiclassical propagator of Herman and Kluk [M. F. Herman and E. Kluk, Chem. Phys. 91, 27 (1984)] is heavily used in a number of molecular and atomic physics contexts, and is of interest here. It is known to be increasingly challenging to implement as the underlying classical chaos strengthens, and we ask whether it is possible to implement it well enough to extract the kind of intermediate time information that reflects wave packet localization at long times. Using a system of two coupled quartic oscillators, we focus on the localizing effects of transport barriers formed by stable and unstable manifolds in the chaotic sea and show that these effects can be captured with the Herman-Kluk propagator.

preprint2008arXiv

Semiclassical theory of non-local statistical measures: residual Coulomb interactions

In a recent letter [Phys. Rev. Lett. {\bf 100}, 164101 (2008)] and within the context of quantized chaotic billiards, random plane wave and semiclassical theoretical approaches were applied to an example of a relatively new class of statistical measures, i.e. measures involving both complete spatial integration and energy summation as essential ingredients. A quintessential example comes from the desire to understand the short-range approximation to the first order ground state contribution of the residual Coulomb interaction. Billiards, fully chaotic or otherwise, provide an ideal class of systems on which to focus as they have proven to be successful in modeling the single particle properties of a Landau-Fermi liquid in typical mesoscopic systems, i.e. closed or nearly closed quantum dots. It happens that both theoretical approaches give fully consistent results for measure averages, but that somewhat surprisingly for fully chaotic systems the semiclassical theory gives a much improved approximation for the fluctuations. Comparison of the theories highlights a couple of key shortcomings inherent in the random plane wave approach. This paper contains a complete account of the theoretical approaches, elucidates the two shortcomings of the oft-relied-upon random plane wave approach, and treats non-fully chaotic systems as well.

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.