In Press
Luukko, P.J.J., et al. Strong quantum scarring by local impurities. Sci. Rep. (In Press). Publisher's VersionAbstract

We discover and characterize strong quantum scars, or eigenstates resembling classical periodic orbits, in two-dimensional quantum wells perturbed by local impurities. These scars are not explained by ordinary scar theory, which would require the existence of short, moderately unstable periodic orbits in the perturbed system. Instead, they are supported by classical resonances in the unperturbed system and the resulting quantum near-degeneracy. Even in the case of a large number of randomly scattered impurities, the scars prefer distinct orientations that extremize the overlap with the impurities. We demonstrate that these preferred orientations can be used for highly efficient transport of quantum wave packets across the perturbed potential landscape. Assisted by the scars, wave-packet recurrences are significantly stronger than in the unperturbed system. Together with the controllability of the preferred orientations, this property may be very useful for quantum transport applications.

Yang, Y., Kolesov, G., Kocia, L. & Heller, E.J. Graphene Terahertz Absorption. arxiv arXiv:1705.06267 (2017). Publisher's VersionAbstract
The unique terahertz properties of graphene has been identified for novel optoelectronic applications. In a graphene sample with bias voltage added, there is an enhanced absorption in the far infrared region and a diminished absorption in the infrared region. The strength of enhancement(diminishment) increases with the gate voltage, and the enhancement compensates the diminishment. We find that it is the coherence length of electrons in graphene that allows pure electronic transitions between states differing by small momentums and makes intraband transition possible, is responsible for the far infrared enhancement. Phonon assisted processes are not necessary and would not in any case contribute to a sum rule. This naturally leads to results obeying the general sum-rule in optical absorptions. Our prediction of the strength of enhancement(diminishment) in terms of the bias agrees with experiments. This is the first direct calculation we are aware of, since the prior phonon assisted model for indirect transition should not obey a sum rule.
Yang, Y., Kolesov, G., Kocia, L. & Heller, E.J. Reassessing Graphene Absorption and Emission Spectroscopy. arxiv arXiv:1704.07500 (2017). Publisher's VersionAbstract

We present a new paradigm for understanding optical absorption and hot electron dynamics experiments in graphene. Our analysis pivots on assigning proper importance to phonon assisted indirect processes and bleaching of direct processes. We show indirect processes figure in the excess absorption in the UV region. Experiments which were thought to indicate ultrafast relaxation of electrons and holes, reaching a thermal distribution from an extremely non-thermal one in under 5-10 fs, instead are explained by the nascent electron and hole distributions produced by indirect transitions. These need no relaxation or ad-hoc energy removal to agree with the observed emission spectra and fast pulsed absorption spectra. The fast emission following pulsed absorption is dominated by phonon assisted processes, which vastly outnumber direct ones and are always available, connecting any electron with any hole any time. Calculations are given, including explicitly calculating the magnitude of indirect processes, supporting these views.

Vardhan, S., Tomasi, G.D., Heyl, M., Heller, E.J. & Pollmann, F. Characterizing time-irreversibility in disordered fermionic systems by the effect of local perturbations. arxiv preprint arXiv:1702.01622 (2017). Publisher's VersionAbstract

We study the effects of local perturbations on the dynamics of disordered fermionic systems in order to characterize time-irreversibility. We focus on three different systems, the non-interacting Anderson and Aubry-Andr\'e-Harper (AAH-) models, and the interacting spinless disordered t-V chain. First, we consider the effect on the full many-body wave-functions by measuring the Loschmidt echo (LE). We show that in the extended/ergodic phase the LE decays exponentially fast with time, while in the localized phase the decay is algebraic. We demonstrate that the exponent of the decay of the LE in the localized phase diverges proportionally to the single-particle localization length as we approach the metal-insulator transition in the AAH model. Second, we probe different phases of disordered systems by studying the time expectation value of local observables evolved with two Hamiltonians that differ by a spatially local perturbation. Remarkably, we find that many-body localized systems could lose memory of the initial state in the long-time limit, in contrast to the non-interacting localized phase where some memory is always preserved.

Pittman, S.M., Tannenbaum, E. & Heller, E.J. Dynamical tunneling versus fast diffusion for a non-convex Hamiltonian. J. Chem. Phys. 145, 054303 (2016). Publisher's VersionAbstract

This paper attempts to resolve the issue of the nature of the 0.01-0.1 cm−1 peak splittings observed in high-resolution IR spectra of polyatomic molecules. One hypothesis is that these splittings are caused by dynamical tunneling, a quantum-mechanical phenomenon whereby energy flows between two disconnected regions of phase-space across dynamical barriers. However, a competing classical mechanism for energy flow is Arnol’d diffusion, which connects different regions of phase-space by a resonance network known as the Arnol’d web. The speed of diffusion is bounded by the Nekhoroshev theorem, which guarantees stability on exponentially long time scales if the Hamiltonian is steep. Here we consider a non-convex Hamiltonian that contains the characteristics of a molecular Hamiltonian, but does not satisfy the Nekhoroshev theorem. The diffusion along the Arnol’d web is expected to be fast for a non-convex Hamiltonian. While fast diffusion is an unlikely competitor for longtime energy flow in molecules, we show how dynamical tunneling dominates compared to fast diffusion in the nearly integrable regime for a non-convex Hamiltonian, as well as present a new kind of dynamical tunneling.

Heller, E.J., et al. Theory of Graphene Raman Scattering. ACS Nano 10, 2, 2803–2818 (2016). Publisher's VersionAbstract

Raman spectroscopy plays a key role in studies of graphene and related carbon systems. Graphene is perhaps the most promising material of recent times for many novel applications, including electronics. In this paper, the traditional and well established Kramers-Heisenberg-Dirac (KHD) Raman scattering theory (1925-1927) is extended to crystalline graphene for the first time. It demands different phonon production mechanisms and phonon energies than does the popular "double resonance" Raman scattering model. The latter has never been compared to KHD. Within KHD, phonons are produced instantly along with electrons and holes, in what we term an electron-hole-phonon triplet, which does not suffer Pauli blocking. A new mechanism for double phonon production we name "transition sliding" explains the brightness of the 2D mode and other overtones, as a result of linear (Dirac cone) electron dispersion. Direct evidence for sliding resides in hole doping experiments performed in 2011 \cite{chenCrommie}. Whole ranges of electronic transitions are permitted and may even constructively interfere for the same laser energy and phonon q, explaining the dispersion, bandwidth, and strength of many two phonon Raman bands. Graphene's entire Raman spectrum, including dispersive and fixed bands, missing bands not forbidden by symmetries, weak bands, overtone bands, Stokes anti-Stokes anomalies, individual bandwidths, trends with doping, and D-2D band spacing anomalies emerge naturally and directly in KHD theory.

Luukko, P.J.J., et al. Strong quantum scarring by local impurities. arXiv preprint arXiv:1511.04198 (2015). Publisher's Version PDF
Baker, S.M. & Heller, E.J. The Degree of Ergodicity of Ortho-and Para-aminobenzonitrile in an Electric Field. The Journal of Physical Chemistry A (2015). Publisher's VersionAbstract

We study the dynamics of the two molecules ortho-aminobenzonitrile (OABN) and para-aminobenzonitrile (PABN). They are structural isomers, with differing asymmetries and dipole moments. In this paper, we show that the dynamics of the system strongly depends on the region of phase space of the initial rotational state, the asymmetry of the molecule, and on the direction of the dipole. We also show that the ergodicity of the system varies gradually with energy, except where the rotational energy of the initial state is much less than the Stark interaction. In this regime, the projection of the dipole along the lab-frame z-axis varies linearly with increasing energy and follows the microcanonical ergodic estimate. Both molecules are far from full chaos for total angular momentum quanta J ∈ [0,45]. However, the initial rotational states in OABN access much more of the available phase space than in PABN. We show that this is a likely cause for the experimental discrepancies in molecular beam deflection experiments.

Kocia, L. & Heller, E.J. Directed HK propagator. The Journal of Chemical Physics 143, 124102 (2015). Publisher's Version PDF
Heller, E.J., Yang, Y. & Kocia, L. Raman Scattering in Carbon Nanosystems: Solving Polyacetylene. ACS Central Science (2015). Publisher's VersionAbstract

Polyacetylene has been a paradigm conjugated organic conductor since well before other conjugated carbon systems such as nanotubes and graphene became front and center. It is widely acknowledged that Raman spectroscopy of these systems is extremely important to characterize them and understand their internal quantum behavior. Here we show, for the first time, what information the Raman spectrum of polyacetylene contains, by solving the 35-year-old mystery of its spectrum. Our methods have immediate and clear implications for other conjugated carbon systems. By relaxing the nearly universal approximation of ignoring the nuclear coordinate dependence of the transition moment (Condon approximation), we find the reasons for its unusual spectroscopic features. When the Kramers–Heisenberg–Dirac Raman scattering theory is fully applied, incorporating this nuclear coordinate dependence, and also the energy and momentum dependence of the electronic and phonon band structure, then unusual line shapes, growth, and dispersion of the bands are explained and very well matched by theory.

Li, Z., Krems, R.V. & Heller, E.J. Collision dynamics of polyatomic molecules containing carbon rings at low temperatures. The Journal of chemical physics 141, 104317 (2014). Publisher's VersionAbstract

We explore the collision dynamics of complex hydrocarbon molecules (benzene, coronene, adamantane, and anthracene) containing carbon rings in a cold buffer gas of 3He. For benzene, we present a comparative analysis of the fully classical and fully quantum calculations of elastic and inelastic scattering cross sections at collision energies between 1 and 10 cm−1. The quantum calculations are performed using the time-independent coupled channel approach and the coupled-states approximation. We show that the coupled-states approximation is accurate at collision energies between 1 and 20 cm−1. For the classical dynamics calculations, we develop an approach exploiting the rigidity of the carbon rings and including low-energy vibrational modes without holonomic constraints. Our results illustrate the effect of the molecular shape and the vibrational degrees of freedom on the formation of long-lived resonance states that lead to low-temperature clustering.

Kocia, L. & Heller, E.J. Communication: HK propagator uniformized along a one-dimensional manifold in weakly anharmonic systems. The Journal of chemical physics 141, 181102 (2014). Publisher's Version PDF
Liu, B. & Heller, E.J. Multiple Scattering and Plasmon Resonance in the Intermediate Regime. arXiv preprint arXiv:1403.4310 (2014). Publisher's VersionAbstract

The collective excitation of the conduction electrons in subwavelength structures gives rise to the Localized Surface Plasmon(LSP). The system consisting of two such LSPs, known as the dimer system,is of fundamental interest and is being actively investigated in the literature. Three regimes have been previously identified and they are the photonic regime, the strong coupling regime and the quantum tunneling regime. In this Letter, we propose a new regime for this intriguing systems, the intermediate regime. In this new regime, the quasistatic approximation, which is widely used to study such LSP systems, fails to capture the main physics: the multiple scattering of the electromagnetic waves between the two LSPs, which significantly modifies the properties of the resonant modes in the system. This intermediate regime provides a new route to explore in plasmonics, where controlling both the excited plasmon modes and the damping rates are of paramount significance.

Liu, B. & Heller, E.J. Shaping Electromagnetic Fields. arXiv preprint arXiv:1405.2807 (2014). Publisher's VersionAbstract

The ability to control electromagnetic fields on the subwavelength scale could open exciting new venues in many fields of science. Transformation optics provides one way to attain such control through the local variation of the permittivity and permeability of a material. Here, we demonstrate another way to shape electromagnetic fields, taking advantage of the enormous size of the configuration space in combinatorial problems and the resonant scattering properties of metallic nanoparticles. Our design does not require the engineering of a material's electromagnetic properties and has relevance to the design of more flexible platforms for probing light-matter interaction and many body physics.

Kotimäki, V., Räsänen, E., Hennig, H. & Heller, E.J. Fractal dynamics in chaotic quantum transport. Physical Review E 88, 022913 (2013). Publisher's VersionAbstract

Despite several experiments on chaotic quantum transport in two-dimensional systems such as semiconductor quantum dots, corresponding quantum simulations within a real-space model have been out of reach so far. Here we carry out quantum transport calculations in real space and real time for a two-dimensional stadium cavity that shows chaotic dynamics. By applying a large set of magnetic fields we obtain a complete picture of magnetoconductance that indicates fractal scaling. In the calculations of the fractality we use detrended fluctuation analysis—a widely used method in time-series analysis—and show its usefulness in the interpretation of the conductance curves. Comparison with a standard method to extract the fractal dimension leads to consistent results that in turn qualitatively agree with the previous experimental data.

Räsänen, E. & Heller, E.J. Optimal control of quantum revival. The European Physical Journal B 86, (2013). Publisher's VersionAbstract


Increasing fidelity is the ultimate challenge of quantum information technology. In addition to decoherence and dissipation, fidelity is affected by internal imperfections such as impurities in the system. Here we show that the quality of quantum revival, i.e., periodic recurrence in the time evolution, can be restored almost completely by coupling the distorted system to an external field obtained from quantum optimal control theory. We demonstrate the procedure with wave-packet calculations in both one- and two-dimensional quantum wells, and analyze the required physical characteristics of the control field. Our results generally show that the inherent dynamics of a quantum system can be idealized at an extremely low cost.


Mason, D.J., Borunda, M.F. & Heller, E.J. Quantum flux and reverse engineering of quantum wave functions. EPL (Europhysics Letters) 102, 60005 (2013). Publisher's VersionAbstract

An interpretation of the probability flux is given, based on a derivation of its eigenstates and relating them to coherent-state projections on a quantum wave function. An extended definition of the flux operator is obtained using coherent states. We present a "processed Husimi" representation, which makes decisions using many Husimi projections at each location. The processed Husimi representation reverse engineers or deconstructs the wave function, yielding the underlying classical ray structure. Our approach makes possible interpreting the dynamics of systems where the probability flux is uniformly zero or strongly misleading. The new technique is demonstrated by the calculation of particle flow maps of the classical dynamics underlying a quantum wave function in simple model systems such as a circular billiard with and without a magnetic field.

Liu, B. & Heller, E.J. Stability of Branched Flow from a Quantum Point Contact. Physical Review Letters 111, 236804 (2013). Full TextAbstract
In classically chaotic systems, small differences in initial conditions are exponentially magnified over time. However, it was observed experimentally that the (necessarily quantum) ‘‘branched flow’’ pattern of electron flux from a quantum point contact (QPC) traveling over a random background potential in two- dimensional electron gases remains substantially invariant to large changes in initial conditions. Since such a potential is classically chaotic and unstable to changes in initial conditions, it was conjectured that the origin of the observed stability is purely quantum mechanical, with no classical analog. In this Letter, we show that the observed stability is a result of the physics of the quantum point contact and the nature of the experiment. We show that the same stability can indeed be reproduced classically, or quantum mechanically. In addition, we explore the stability of the branched flow with regards to changes in the eigenmodes of the quantum point contact.
Borunda, M.F., Hennig, H. & Heller, E.J. Ballistic versus diffusive transport in graphene. Physical Review B 88, 125415 (2013). Publisher's VersionAbstract

We investigate the transport of electrons in disordered and pristine graphene devices. Fano shot noise, a standard metric to assess the mechanism for electronic transport in mesoscopic devices, has been shown to produce almost the same magnitude (1/3) in ballistic and diffusive graphene devices and is therefore of limited applicability. We consider a two-terminal geometry where the graphene flake is contacted by narrow metallic leads. We propose that the dependence of the conductance on the position of one of the leads, a conductance profile, can give us insight into the charge flow, which can in turn be used to analyze the transport mechanism. Moreover, we simulate scanning probe microscopy (SPM) measurements for the same devices, which can visualize the flow of charge inside the device, thus complementing the transport calculations. From our simulations, we find that both the conductance profile and SPM measurements are excellent tools to assess the transport mechanism differentiating ballistic and diffusive graphene systems.

Kocia, L. & Heller, E.J. Generalized dephasing relation for fidelity and application as an efficient propagator. The Journal of Chemical Physics 139, 124110 (2013). Publisher's VersionAbstract

The dephasing relation (DR), a linearization of semiclassical fidelity, is generalized to include the overlap of “off-diagonal” elements. The accuracy of the formulation is tested in integrable and chaotic systems and its scaling with dimensionality is studied in a Caldeira-Leggett model with many degrees of freedom. It is shown that the DR is often in very good agreement with numerically analytic quantum results and frequently outperforms an alternative semiclassical treatment. Most importantly, since there is no computationally expensive prefactor, and Monte Carlo Metropolis sampling is used to facilitate the calculation, the DR is found to scale remarkably well with increasing dimension. We further demonstrate that a propagator based on the DR can include more quantum coherence and outperform other popular linearized semiclassical methods, such as forward-backward semiclassical dynamics (FBSD) and the linearized semiclassical initial value representation (LSC-IVR).