Source author record

Stephan Haas

Stephan Haas 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

52works
14topics
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

52 published item(s)

preprint2024arXiv

SNS Junctions along the BCS-BEC Crossover

We present a theory of SNS junctions, a normal metal sandwiched between two superconductors, along the crossover from the BCS to the BEC regime. We calculate the Josephson current as a function of the chemical potential relative to the band edge in the superconducting region, $μ_S$, where the BEC phase is indicated by $μ_S <0$. The chemical potential relative to the band edge in the normal metal, $μ_N$, allows us to tune the junction between the SNS case ($μ_N>0$) and the SIS case, where the superconductors are separated by a tunneling barrier. We find that there are Andreev levels in the BEC regime, as long as there is sufficient density of states in the normal region, i.e. when $μ_N>Δ$, where $Δ$ is the amplitude of the superconducting order parameter. For 1D SNS junctions, we find the Josephson current $I_S$ carried by these Andreev levels to be a function of the ratio $Δ/Δ_d$, where $Δ_d$ is the Andreev level spacing. At zero temperature, the Josephson current has a maximum on the BCS side of the transition where $Δ$ is maximal. At finite temperature, however, we find that the maximum moves to the BEC side of the crossover. We identify the mechanism for this phenomenon to be the decrease in the number of Andreev levels at the BCS-BEC crossover, accompanied by an increase in excitation energy to the unoccupied levels, making it less likely that these states are thermally occupied. Thereby, at finite temperature, the Josephson current is more strongly reduced on the BCS side of the crossover, resulting in a maximal Josephson current at the BCS-BEC crossover.

preprint2022arXiv

Collision Dynamics of Bose-Einstein Condensates

We study the collision dynamics of two Bose-Einstein condensates with their dynamical wave functions modeled by a set of coupled, time-dependent Gross-Pitaevskii equations. Beginning with an effective one-dimensional system, we identify regimes characterized by the relationship between inter- and intra-atomic interactions and the initial configuration of the system, akin to the equilibrium phase diagram of two interacting Bose condensates. We consider a dynamical setup in which two wave packets are initially at rest with a small separation about the center of an asymmetric harmonic trap. Upon release, we observe a rapid approach to dynamical equilibrium in the limits of very large and very small inter-particle repulsion, characterized by periodic transmission or reflection of the condensates as distinguishable units, whereas the intermediate, critical regime is characterized by extended transient dynamics, density fracturing, and dynamical mixing. We therefore have reason to believe that non-trivial behavior may exist during the collisions of Bose-Einstein condensates as a result of their non-linear interactions, and these effects may be observable in a suitable laboratory environment.

preprint2022arXiv

Emulating Quantum Dynamics with Neural Networks via Knowledge Distillation

High-fidelity quantum dynamics emulators can be used to predict the time evolution of complex physical systems. Here, we introduce an efficient training framework for constructing machine learning-based emulators. Our approach is based on the idea of knowledge distillation and uses elements of curriculum learning. It works by constructing a set of simple, but rich-in-physics training examples (a curriculum). These examples are used by the emulator to learn the general rules describing the time evolution of a quantum system (knowledge distillation). The goal is not only to obtain high-quality predictions, but also to examine the process of how the emulator learns the physics of the underlying problem. This allows us to discover new facts about the physical system, detect symmetries, and measure relative importance of the contributing physical processes. We illustrate this approach by training an artificial neural network to predict the time evolution of quantum wave packages propagating through a potential landscape. We focus on the question of how the emulator learns the rules of quantum dynamics from the curriculum of simple training examples and to which extent it can generalize the acquired knowledge to solve more challenging cases.

preprint2021arXiv

Probing Majorana Modes via Local Spin Dynamics

We investigate Majorana modes in a quantum spin chain with bond-dependent exchange interactions by studying its dynamics. Specifically, we consider two-time correlations for the Kitaev-Heisenberg (KH) Hamiltonian close to the so-called Kitaev critical point. Here, the model coincides with a phase boundary of two uncoupled instances of Kitaev's model for p-wave superconductors, together supporting a degenerate ground state characterized by multiple Majorana modes. In this regime, the real-time dynamics of local spins reveal a set of strong zero modes, corresponding to a set of protruding frequencies in the two-time correlation function. We derive perturbative interactions that map the KH spin chain onto the topological regime of Kitaev's fermionic model, thus opening up a bulk gap whilst retaining almost degenerate modes in the mesoscopic regime, i.e., for finite system sizes. This showcases the emergence of Majorana modes in a chain of effective dimers. Here, the binding energy within each unit cell competes with the inter-dimer coupling to generate a finite size energy gap, in analogy with local energy terms in the transverse-field Ising model. These modes give rise to long coherence times of local spins located at the system edges. By breaking the local symmetry in each dimer, one can also observe a second class of Majorana modes in terms of a beating frequency in the two-time correlations function of the edge spin. Furthermore, we develop a scenario for realizing these model predictions in ion-trap quantum simulators with collective addressing of the ions.

preprint2020arXiv

Detection of Topology via Entanglement Oscillations

We introduce a characterization of topological order based on bulk oscillations of the entanglement entropy and the definition of an `entanglement gap', showing that it is generally applicable to pure and disordered quantum systems. Using exact diagonalization and the strong disorder renormalization group method, we demonstrate that this approach gives results in agreement with the use of traditional topological invariants, especially in cases where topological order is known to persist in the presence of off-diagonal bond disorder. The entanglement gap is then used to analyze classes of quantum systems with alternating bond types, allowing us to construct their topological phase diagrams. The validity of these phase diagrams is verified in limiting cases of dominant bond types, where the solution is known.

preprint2019arXiv

Induced superconducting pair correlations in a quasicrystal coupled to a BCS superconductor

In this report, we describe the proximity effect which arises when a quasicrystal is placed in contact with a superconductor. We consider the simplest known model of a quasicrystal, the 1D Fibonacci chain, for which all states are known to be critical. A hybrid ring is made by connecting a finite piece of such a Fibonacci chain to a BCS-type superconductor. Solving the resulting Bogoliubov-de Gennes equations self-consistently, we show that the proximity effect in the Fibonacci chain is long ranged, and that the induced superconducting order carries information on topological properties, namely, the winding numbers of edge states of the Fibonacci chain.

preprint2019arXiv

Localized Plasmons in Topological Insulators

We investigate in a fully quantum-mechanical manner how the many-body excitation spectrum of topological insulators is affected by the presence of long-range Coulomb interactions. In the one-dimensional Su-Schrieffer-Heeger model and its mirror-symmetric variant strongly localized plasmonic excitations are observed which originate from topologically non-trivial single-particle states. These \textit{"topological plasmons"} inherit some of the characteristics of their constituent topological single-particle states, but they are not equally well protected against disorder due to the admixture of non-topological bulk single-particle states in the polarization function. The strength of the effective Coulomb interactions is also shown to have strong effects on the plasmonic modes. Furthermore, we show how external modifications via dielectric screening and applied electric fields with distinct symmetries can be used to study topological plasmons, thus allowing for experimental verification of our atomistic predictions.

preprint2015arXiv

Small quench dynamics as a probe for trapped ultracold atoms

Finite systems of bosons and/or fermions described by the Hubbard model can be realized using ultracold atoms confined in optical lattices. The ground states of these systems often exhibit a coexistence of compressible superfluid and incompressible Mott insulating regimes. We analyze such systems by studying the out-of-equilibrium dynamics following a weak sudden quench of the trapping potential. In particular, we show how the temporal variance of the site occupations reveals the location of spatial boundaries between compressible and incompressible regions. The feasibility of this approach is demonstrated for several models using numerical simulations. We first consider integrable systems, hard-core bosons (spinless fermions) confined by a harmonic potential, where space separated Mott and superfluid phases coexist. Then, we analyze a nonintegrable system, a $J-V-V'$ model with coexisting charge density wave and superfluid phases. We find that the temporal variance of the site occupations is a more effective measure than other standard indicators of phase boundaries such as a local compressibility. Based on these examples, we argue that analyzing temporal fluctuations is a valuable experimental tool for exploring phase boundaries in trapped atom systems.

preprint2014arXiv

Driven dipole oscillations and the lowest energy excitations of strongly interacting lattice bosons in a harmonic trap

We show that the analysis of the time evolution of the occupation of site and momentum modes of harmonically trapped lattice hard-core bosons, under driven dipole oscillations, allows one to determine the energy of the lowest one-particle excitations of the system in equilibrium. The analytic solution of a single particle in the absence of a lattice is used to identify which function of those time-dependent observables is best fit for the analysis, as well as to relate the dynamic response of the system to its single-particle spectrum. In the presence of the lattice and of multiple particles, a much richer and informative dynamical response is observed under the drive.

preprint2014arXiv

Quantification and Control of non-Markovian Evolution in Finite Quantum Systems via Feedback

We consider the unitary time evolution of continuous quantum mechanical systems confined to a cavity in contact with a finite bath of variable size. Measures for Markovianity for such finite system-bath configurations are developed in terms of Hilbert-Schmidt distances of time evolving wave packets. The relevant time scales are identified, which characterize pseudo-Markovian transient behavior, boundary scattering induced non-Markovian oscillations at intermediate times, and non-Markovian rephasing events at long time scales. It is shown how these time scales can be controlled by tunable parameters such as the bath size and the strength of the system-bath coupling.

preprint2014arXiv

Two-dimensional percolation with multiple seeds

We study non-uniform percolation in a two-dimensional cluster growth model with multiple seeds. With increasing concentration of seeds, the percolation threshold is found to increase monotonically, while the exponents for correlation length, order parameter, and average cluster size, keep invariant. The scaling law for an infinite square lattice keeps working for any nonzero concentration of seeds. Abnormal finite-size scaling behaviours happen at low concentration of seeds.

preprint2013arXiv

Bond Disorder in Even-Leg Heisenberg Ladders

Random bond disorder in antiferromagnetic spin-1/2 Heisenberg ladders is investigated using the Stochastic Series Expansion Quantum Monte Carlo method. We find that the effects of individual bond impurities vary strongly, depending on their position on legs and rungs. We initially focus on how the distribution of local bond energies depends on the impurity concentration. Then we study how the phase diagram of even-leg ladders is affected by random bond doping. We observe Bose glass phases in two regimes ($h' < h \lesssim h_{c1}$ and $h" < h <h_{c2}$) and a Bose-Einstein Condensate in-between. Their presence are discussed in the relation to the local bond energies.

preprint2013arXiv

Crossover physics in the non-equilibrium dynamics of quenched quantum impurity systems

A general framework is proposed to tackle analytically local quantum quenches in integrable impurity systems, combining a mapping onto a boundary problem with the form factor approach to boundary-condition-changing operators introduced in Phys. Rev. Lett. 80, 4370 (1998). We discuss how to compute exactly two central quantities of interest: the Loschmidt echo and the distribution of the work done during the quantum quench. Our results display an interesting crossover physics characterized by the energy scale T_b of the impurity corresponding to the Kondo temperature. We discuss in detail the non-interacting case as a paradigm and benchmark for more complicated integrable impurity models, and check our results using numerical methods.

preprint2013arXiv

Fast realization of a spatially correlated percolation model

We propose two schemes to achieve fast realizations of spatially correlated percolation models. The schemes are shown to be efficient in complementary regimes of correlation phase space. They are combined with a generalized Newman-Ziff algorithm to numerically determine the percolation thresholds of two-dimensional lattices in the presence of correlations. It is found that the spatial correlations affect only a relatively small part of phase space.

preprint2013arXiv

Non-equilibrium Approach to Thermodynamics using Jarzynski Equality and Diagonal Entropy

We combine the formalisms of diagonal entropy and Jarzynski Equality to study the thermodynamic properties of closed quantum systems. Applying this approach to a quantum harmonic oscillator, the diagonal entropy offers a notion of temperature for closed systems away from equilibrium, and allows computing free-energy profiles. We also apply this approach to a hard-core boson lattice model, and discuss measures how to estimate temperature, entropy and measure work distribution functions. This technique offers a path to investigate the non-equilibrium quantum dynamics by means of performing work in a series of quenches.

preprint2013arXiv

Non-Markovian Equilibration Controlled by Symmetry Breaking

We study the effects of symmetry breaking on non-Markovian dynamics in various system-bath arrangements. It is shown that by breaking certain symmetries features signaling non-Markovian time evolution disappear within a finite time t_{g}. We demonstrate numerically that the scaling of t_{g} with the symmetry breaking strength is different for various types of symmetry. We provide a mathematical explanation for these differences related to the spectrum of the total system-bath Hamiltonian and provide arguments that the scaling properties of t_{g} should be universal.

preprint2013arXiv

Possible Nematic Order Driven by Magnetic Fluctuations in Iron Pnictides

In this paper, instabilities of the isotropic metallic phase in iron pnictides are investigated. The relevant quartic fermionic interaction terms in the model are identified using phase space arguments. Using the functional integral formalism, a Hubbard-Stratonovich transformation is used to decouple these quartic terms. This procedure introduces several bosonic fields which describe the low-energy collective modes of the system. By studying the behavior of these collective modes, a possible instability is found in the forward scattering channel of the isotropic phase driven by magnetic fluctuations. Using mean field analysis, we obtain a static and homogeneous ground state. This ground state is metallic, but the electron Fermi pockets are distorted unequally at different pockets in momentum space. This results in a desirable nematic ordering which breaks the lattice C4 symmetry but preserves translational symmetry and may explain several experimental observations.

preprint2013arXiv

Scalable universal holonomic quantum computation realized with an adiabatic quantum data bus and potential implementation using superconducting flux qubits

In this paper we examine the use of an adiabatic quantum data transfer protocol to build a universal quantum computer. Single qubit gates are realized by using a bus protocol to transfer qubits of information down a spin chain with a unitary twist. This twist arises from altered couplings on the chain corresponding to unitary rotations performed on one region of the chain. We show how a controlled NOT gate can be realized by using a control qubit with Ising type coupling. The method discussed here can be extended to non-adiabatic quantum bus protocols. We also examine the potential of realizing such a quantum computer by using superconducting flux qubits.

preprint2012arXiv

Correlations in Quantum Spin Ladders with Site and Bond Dilution

We investigate the effects of quenched disorder, in the form of site and bond dilution, on the physics of the $S=1/2$ antiferromagnetic Heisenberg model on even-leg ladders. Site dilution is found to prune rung singlets and thus create localized moments which interact via a random, unfrustrated network of effective couplings, realizing a random-exchange Heisenberg model (REHM) in one spatial dimension. This system exhibits a power-law diverging correlation length as the temperature decreases. Contrary to previous claims, we observe that the scaling exponent is non-universal, i.e., doping dependent. This finding can be explained by the discrete nature of the values taken by the effective exchange couplings in the doped ladder. Bond dilution on even-leg ladders leads to a more complex evolution with doping of correlations, which are weakly enhanced in 2-leg ladders, and are even suppressed for low dilution in the case of 4-leg and 6-leg ladders. We clarify the different aspects of correlation enhancement and suppression due to bond dilution by isolating the contributions of rung-bond dilution and leg-bond dilution.

preprint2012arXiv

Demonstration of Jarzynski's Equality in Open Quantum Systems Using a Step-Wise Pulling Protocol

We present a generalization of Jarzynski's Equality, applicable to quantum systems, relating discretized mechanical work and free-energy changes. The theory is based on a step-wise pulling protocol. We find that work distribution functions can be constructed from fluctuations of a reaction coordinate along a reaction pathway in the step-wise pulling protocol. We also propose two sets of equations to determine the two possible optimal pathways that provide the most significant contributions to free-energy changes. We find that the transitions along these most optimal pathways, satisfying both sets of equations, follow the principle of detailed balance. We then test the theory by explicitly computing the free-energy changes for a one-dimensional quantum harmonic oscillator. This approach suggests a feasible way of measuring the fluctuations to experimentally test Jarzynski's Equality in many-body systems, such as Bose-Einstein condensates.

preprint2012arXiv

Equilibration times in clean and noisy systems

We study the equilibration dynamics of closed finite quantum systems and address the question of the time needed for the system to equilibrate. In particular we focus on the scaling of the equilibration time T_{\mathrm{eq}} with the system size L . For clean systems we give general arguments predicting T_{\mathrm{eq}}=O(L^{0}) for clustering initial states, while for small quenches around a critical point we find T_{\mathrm{eq}}=O(L^ζ) where ζis the dynamical critical exponent. We then analyze noisy systems where exponentially large time scales are known to exist. Specifically we consider the tight-binding model with diagonal impurities and give numerical evidence that in this case T_{\mathrm{eq}}\sim Be^{CL^ψ} where B,C, ψare observable dependent constants. Finally, we consider another noisy system whose evolution dynamics is randomly sampled from a circular unitary ensemble. Here, we are able to prove analytically that T_{\mathrm{eq}}=O(1), thus showing that noise alone is not sufficient for slow equilibration dynamics.

preprint2012arXiv

Excitonic Instabilities and Insulating States in Bilayer Graphene

The competing ground states of bilayer graphene are studied by applying renormalization group techniques to a bilayer honeycomb lattice with nearest neighbor hopping. In the absence of interactions, the Fermi surface of this model at half-filling consists of two nodal points with momenta $\mathbf{K}$, $\mathbf{K}'$, where the conduction band and valence band touch each other, yielding a semi-metal. Since near these two points the energy dispersion is quadratic with perfect particle-hole symmetry, excitonic instabilities are inevitable if inter-band interactions are present. Using a perturbative renormalization group analysis up to the one-loop level, we find different competing ordered ground states, including ferromagnetism, superconductivity, spin and charge density wave states with ordering vector $\mathbf{Q}=\mathbf{K}-\mathbf{K}'$, and excitonic insulator states. In addition, two states with valley symmetry breaking are found in the excitonic insulating and ferromagnetic phases. This analysis strongly suggests that the ground state of bilayer graphene should be gapped, and with the exception of superconductivity, all other possible ground states are insulating.

preprint2012arXiv

Quantum Hall States in Graphene from Strain-Induced Nonuniform Magnetic Fields

We examine strain-induced quantized Landau levels in graphene. Specifically, arc-bend strains are found to cause nonuniform pseudomagnetic fields. Using an effective Dirac model which describes the low-energy physics around the nodal points, we show that several of the key qualitative properties of graphene in a strain-induced pseudomagnetic field are different compared to the case of an externally applied physical magnetic field. We discuss how using different strain strengths allows us to spatially separate the two components of the pseudospinor on the different sublattices of graphene. These results are checked against a tight-binding calculation on the graphene honeycomb lattice, which is found to exhibit all the features described. Furthermore, we find that introducing a Hubbard repulsion on the mean-field level induces a measurable polarization difference between the A and the B sublattices, which provides an independent experimental test of the theory presented here.

preprint2012arXiv

Using the J1-J2 Quantum Spin Chain as an Adiabatic Quantum Data Bus

This paper investigates numerically a phenomenon which can be used to transport a single q-bit down a J1-J2 Heisenberg spin chain using a quantum adiabatic process. The motivation for investigating such processes comes from the idea that this method of transport could potentially be used as a means of sending data to various parts of a quantum computer made of artificial spins, and that this method could take advantage of the easily prepared ground state at the so called Majumdar-Ghosh point. We examine several annealing protocols for this process and find similar result for all of them. The annealing process works well up to a critical frustration threshold.

preprint2011arXiv

Bose glass and Mott glass of quasiparticles in a doped quantum magnet

The low-temperature states of bosonic fluids exhibit fundamental quantum effects at the macroscopic scale: the best-known examples are Bose-Einstein condensation (BEC) and superfluidity, which have been tested experimentally in a variety of different systems. When bosons are interacting, disorder can destroy condensation leading to a so-called Bose glass. This phase has been very elusive to experiments due to the absence of any broken symmetry and of a finite energy gap in the spectrum. Here we report the observation of a Bose glass of field-induced magnetic quasiparticles in a doped quantum magnet (Br-doped dichloro-tetrakis-thiourea-Nickel, DTN). The physics of DTN in a magnetic field is equivalent to that of a lattice gas of bosons in the grand-canonical ensemble; Br-doping introduces disorder in the hoppings and interaction strengths, leading to localization of the bosons into a Bose glass down to zero field, where it acquires the nature of an incompressible Mott glass. The transition from the Bose glass (corresponding to a gapless spin liquid) to the BEC (corresponding to a magnetically ordered phase) is marked by a novel, universal exponent governing the scaling on the critical temperature with the applied field, in excellent agreement with theoretical predictions. Our study represents the first, quantitative account of the universal features of disordered bosons in the grand-canonical ensemble.

preprint2011arXiv

Defect Induced Resonances and Magnetic Patterns in Graphene

We investigate the effects of point and line defects in monolayer graphene within the framework of the Hubbard model, using a self-consistent mean field theory. These defects are found to induce characteristic patterns into the electronic density of states and cause non-uniform distributions of magnetic moments in the vicinity of the impurity sites. Specifically, defect induced resonance bound states in the local density of states are observed at energies close to the Dirac points. The magnitudes of the frequencies of these resonance states are shown to decrease with the strength of the scattering potential, whereas their amplitudes decay algebraically with increasing distance from the defect. For the case of defect clusters, we observe that with increasing defect cluster size the local magnetic moments in the vicinity of the cluster center are strongly enhanced. Furthermore, non-trivial impurity induced magnetic patterns are observed in the presence of line defects: zigzag line defects are found to introduce stronger-amplitude magnetic patterns than armchair line defects. When the scattering strength of these topological defects is increased, the induced patterns of magnetic moments become more strongly localized.

preprint2011arXiv

Entanglement and its evolution after a quench in the presence of an energy current

We study the Ising spin chain with a Dzyaloshinskii-Moriya interaction focusing on the static and dynamic properties of the entanglement entropy, following a quantum quench. We show that the effects of the additional anisotropic interaction on the phase diagram and on the dynamics of the system are captured by the properties of the entanglement entropy. In particular, the model provides a way to study the quench dynamics in a system with an energy current. We consider quenches starting from an initial excited state of the Ising spin chain, and we analyze the effects of different initial conditions.

preprint2011arXiv

Propagation of Disturbances in Degenerate Quantum Systems

Disturbances in gapless quantum many-body models are known to travel an unlimited distance throughout the system. Here, we explore this phenomenon in finite clusters with degenerate ground states. The specific model studied here is the one-dimensional J1-J2 Heisenberg Hamiltonian at and close to the Majumdar-Ghosh point. Both open and periodic boundary conditions are considered. Quenches are performed using a local magnetic field. The degenerate Majumdar-Ghosh ground state allows disturbances which carry quantum entanglement to propagate throughout the system, and thus dephase the entire system within the degenerate subspace. These disturbances can also carry polarization, but not energy, as all energy is stored locally. The local evolution of the part of the system where energy is stored drives the rest of the system through long-range entanglement. We also examine approximations for the ground state of this Hamiltonian in the strong field limit, and study how couplings away from the Majumdar-Ghosh point affect the propagation of disturbances. We find that even in the case of approximate degeneracy, a disturbance can be propagated throughout a finite system.

preprint2011arXiv

Statistical properties of random matrix product states

We study the set of random matrix product states (RMPS) introduced in arXiv:0908.3877 as a tool to explore foundational aspects of quantum statistical mechanics. In the present work, we provide an accurate numerical and analytical investigation of the properties of RMPS. We calculate the average state of the ensemble in the non-homogeneous case, and numerically check the validity of this result. We also suggest using RMPS as a tool to approximate properties of general quantum random states. The numerical simulations presented here support the accuracy and efficiency of this approximation. These results suggest that any generalized canonical state can be approximated with high probability by the reduced density matrix of a random MPS, if the average MPS coincide with the associated microcanonical ensemble.

preprint2011arXiv

Visualization of nano-plasmons in graphene

We study localized plasmons at the nanoscale (nano-plasmons) in graphene. The collective excitations of induced charge density modulations in graphene are drastically changed in the vicinity of a single impurity compared to graphene's bulk behavior. The dispersion of nano-plasmons depends on the number of electrons and the sign, strength and size of the impurity potential. Due to this rich parameter space the calculated dispersions are intrinsically multidimensional requiring an advanced visualization tool for their efficient analysis, which can be achieved with parallel rendering. To overcome the problem of analyzing thousands of very complex spatial patterns of nano-plasmonic modes, we take a combined visual and quantitative approach to investigate the excitations on the two-dimensional graphene lattice. Our visual and quantitative analysis shows that impurities trigger the formation of localized plasmonic excitations of various symmetries. We visually identify dipolar, quadrupolar and radial modes, and quantify the spatial distributions of induced charges.

preprint2010arXiv

Impurity assisted nanoscale localization of plasmonic excitations in graphene

The plasmon modes of pristine and impurity doped graphene are calculated, using a real-space theory which determines the non-local dielectric response within the random phase approximation. A full diagonalization of the polarization operator is performed, allowing the extraction of all its poles. It is demonstrated how impurities induce the formation of localized modes which are absent in pristine graphene. The dependence of the spatial modulations over few lattice sites and frequencies of the localized plasmons on the electronic filling and impurity strength is discussed. Furthermore, it is shown that the chemical potential and impurity strength can be tuned to control target features of the localized modes. These predictions can be tested by scanning tunneling microscopy experiments.

preprint2010arXiv

Local quenches in frustrated quantum spin chains: global vs. subsystem equilibration

We study the equilibration behavior following local quenches, using frustrated quantum spin chains as an example of interacting closed quantum systems. Specifically, we examine the statistics of the time series of the Loschmidt echo, the trace distance of the time-evolved local density matrix to its average state, and the local magnetization. Depending on the quench parameters, the equilibration statistics of these quantities show features of good or poor equilibration, indicated by Gaussian, exponential or bistable distribution functions. These universal functions provide valuable tools to characterize the various time-evolution responses and give insight into the plethora of equilibration phenomena in complex quantum systems.

preprint2010arXiv

Magnetic Bose glass phases of coupled antiferromagnetic dimers with site dilution

We numerically investigate the phase diagram of two-dimensional site-diluted coupled dimer systems in an external magnetic field. We show that this phase diagram is characterized by the presence of an extended Bose glass, not accessible to mean-field approximation, and stemming from the localization of two distinct species of bosonic quasiparticles appearing in the ground state. On the one hand, non-magnetic impurities doped into the dimer-singlet phase of a weakly coupled dimer system are known to free up local magnetic moments. The deviations of these local moments from full polarization along the field can be mapped onto a gas of bosonic quasiparticles, which undergo condensation in zero and very weak magnetic fields, corresponding to transverse long-range antiferromagnetic order. An increasing magnetic field lowers the density of such quasiparticles to a critical value at which a quantum phase transition occurs, corresponding to the quasiparticle localization on clusters of local magnets (dimers, trimers, etc.) and to the onset of a Bose glass. Strong finite-size quantum fluctuations hinder further depletion of quasiparticles from such clusters, and thus lead to the appearance of pseudo-plateaus in the magnetization curve of the system. On the other hand, site dilution hinders the field-induced Bose-Einstein condensation of triplet quasiparticles on the intact dimers, and it introduces instead a Bose glass of triplets. A thorough numerical investigation of the phase diagram for a planar system of coupled dimers shows that the two above-mentioned Bose glass phases are continuously connected, and they overlap in a finite region of parameter space, thus featuring a two-species Bose glass. The quantum phase transition from Bose glass to magnetic order in two dimensions is marked by novel universal exponents.

preprint2009arXiv

Adaptation and Performance of the Cartesian Coordinates Fast Multipole Method for Nanomagnetic Simulations

An implementation of the fast multiple method (FMM) is performed for magnetic systems with long-ranged dipolar interactions. Expansion in spherical harmonics of the original FMM is replaced by expansion of polynomials in Cartesian coordinates, which is considerably simpler. Under open boundary conditions, an expression for multipole moments of point dipoles in a cell is derived. These make the program appropriate for nanomagnetic simulations, including magnetic nanoparticles and ferrofluids. The performance is optimized in terms of cell size and parameter set (expansion order and opening angle) and the trade off between computing time and accuracy is quantitatively studied. A rule of thumb is proposed to decide the appropriate average number of dipoles in the smallest cells, and an optimal choice of parameter set is suggested. Finally, the superiority of Cartesian coordinate FMM is demonstrated by comparison to spherical harmonics FMM and FFT.

preprint2009arXiv

Phase Diagram of Magnetization Reversal Processes in Nanorings

We investigate magnetization reversal processes of magnetic nanorings. Using a recently developed efficient cartesian coordinates fast multipole method, magnetization switching phase diagrams are constructed for such structures. We find that the dominant switching mechanisms strongly depend on the relative spatial dimensions of the rings, i.e. height, width and radius. Furthermore, they are found to depend on a characteristic exchange length scale, which is determined by the competing short-range ferromagnetic and long-range dipolar interactions in these systems. Detailed simulations allow us to identify three novel switching mechanisms, which should be useful for building high density storage systems.

preprint2009arXiv

Plasmonic Excitations in Tight-Binding Nanostructures

We explore the collective electromagnetic response in atomic clusters of various sizes and geometries. Our aim is to understand, and hence to control, their dielectric response, based on a fully quantum-mechanical description which captures accurately their relevant collective modes. The electronic energy levels and wave functions, calculated within the tight-binding model, are used to determine the non-local dielectric response function. It is found that the system shape, the electron filling and the driving frequency of the external electric field strongly control the resonance properties of the collective excitations in the frequency and spatial domains. Furthermore, it is shown that one can design spatially localized collective excitations by properly tailoring the nanostructure geometry.

preprint2009arXiv

Revealing Novel Quantum Phases in Quantum Antiferromagnets on Random Lattices

Quantum magnets represent an ideal playground for the controlled realization of novel quantum phases and of quantum phase transitions. The Hamiltonian of the system can be indeed manipulated by applying a magnetic field or pressure on the sample. When doping the system with non-magnetic impurities, novel inhomogeneous phases emerge from the interplay between geometric randomness and quantum fluctuations. In this paper we review our recent work on quantum phase transitions and novel quantum phases realized in disordered quantum magnets. The system inhomogeneity is found to strongly affect phase transitions by changing their universality class, giving the transition a novel, quantum percolative nature. Such transitions connect conventionally ordered phases to unconventional, quantum disordered ones - quantum Griffiths phases, magnetic Bose glass phases - exhibiting gapless spectra associated with low-energy localized excitations.

preprint2009arXiv

Universal phase diagram of disordered bosons from a doped quantum magnet

A quantum particle cannot in general diffuse through a disordered medium because of its wavelike nature, but interacting particles can escape localization by collectively percolating through the system. For bosonic particles this phenomenon corresponds to a quantum transition from a localized insulator phase - the Bose glass - to a superfluid phase, in which particles condense into an extended state. Here, we construct a universal phase diagram of disordered bosons in doped quantum magnets for which bosonic quasi-particles are represented by magnetized states (spin triplets) of the quantum spins, condensing into a magnetically ordered state. The appearance of a Bose glass leads to strong measurable signatures in the onset of superfluidity of the spin-triplet gas, exhibiting a complex crossover from low-temperature quantum percolation to a conventional condensation transition at intermediate temperatures.

preprint2006arXiv

Bottom-up approach to high-temperature superconductivity

Since the discovery of high-temperature superconductivity in the cuprates a theoretical understanding of their phase diagram has remained one of the major outstanding problems in condensed matter physics. Here we propose an effective low-energy Hamiltonian which produces both d-wave density wave (dDW) and d-wave superconducting (dSC) solutions within the BCS mean-field theory. This model predicts that (a) the observed pseudogap phase is a dDW state, (b) the superconducting phase is a d-wave BCS state, and (c) in the underdoped regime there is a gossamer superconducting state, i.e. dSC in coexistence with dDW. Moreover, this theory naturally explains the Uemura relation, the reduction of the quasiparticle density of states at the Fermi level, and the salient features in the tunneling conductivity measured in underdoped Bi2212.

preprint2006arXiv

Mott glass in site-diluted S=1 antiferromagnets with single-ion anisotropy

The interplay between site dilution and quantum fluctuations in S=1 Heisenberg antiferromagnets on the square lattice is investigated using quantum Monte Carlo simulations. Quantum fluctuations are tuned by a single-ion anisotropy D. In the clean limit, a sufficiently large D>D_c = 5.65(2) J forces each spin into its m_S=0 state, and thus destabilizes antiferromagnetic order. In the presence of site dilution, quantum fluctuations are found to destroy Néel order before the percolation threshold of the lattice is reached, if D exceeds a critical value D^* = 2.3(2) J. This mechanism opens up an extended quantum-disordered Mott glass phase on the percolated lattice, characterized by a gapless spectrum and vanishing uniform susceptibility.

preprint2006arXiv

Scaling Behavior of Entanglement in Two- and Three-Dimensional Free Fermions

Exactly solving a spinless fermionic system in two and three dimensions, we investigate the scaling behavior of the block entropy in critical and non-critical phases. The scaling of the block entropy crucially depends on the nature of the excitation spectrum of the system and on the topology of the Fermi surface. Noticeably, in the critical phases the scaling violates the area law and acquires a logarithmic correction \emph{only} when a well defined Fermi surface exists in the system. When the area law is violated, we accurately verify a conjecture for the prefactor of the logarithmic correction, proposed by D. Gioev and I. Klich [quant-ph/0504151].

preprint2006arXiv

Small $T_1^{-1}$ coherence peak near $T_c$ in unconventional superconductors

It is usually believed that a coherence peak just below T$_{c}$ in the nuclear spin lattice relaxation rate T$_{1}^{-1}$ in superconducting materials is a signature of conventional s-wave pairing. In this paper we demonstrate that any unconventional superconductor obeying BCS pure-case weak-coupling theory should show a small T$_{1}^{-1}$ coherence peak near T$_{c}$, generally with a height between 3 and 15 percent greater than the normal state T$_{1}^{-1}$ at T$_{c}$. It is largely due to impurity effects that this peak has not commonly been observed.

preprint2005arXiv

Bose-Einstein condensation vs. localization of bosonic quasiparticles in disordered weakly-coupled dimer antiferromagnets

We investigate the field-induced insulator-to-superfluid transition of bosonic quasiparticles in $S=1/2$ weakly-coupled dimer antiferromagnets. In presence of realistic disorder due to site dilution of the magnetic lattice, we show that the system displays an extended Bose-glass phase characterized by the localization of the hard-core quasiparticles.

preprint2005arXiv

Quantum localization in bilayer Heisenberg antiferromagnets with site dilution

The field-induced antiferromagnetic ordering in systems of weakly coupled S=1/2 dimers at zero temperature can be described as a Bose-Einstein condensation of triplet quasiparticles (singlet quasiholes) in the ground state. For the case of a Heisenberg bilayer, it is here shown how the above picture is altered in the presence of site dilution of the magnetic lattice. Geometric randomness leads to quantum localization of the quasiparticles/quasiholes and to an extended Bose-glass phase in a realistic disordered model. This localization phenomenon drives the system towards a quantum-disordered phase well before the classical geometric percolation threshold is reached.

preprint2005arXiv

Quantum percolation in two-dimensional antiferromagnets

The interplay of geometric randomness and strong quantum fluctuations is an exciting topic in quantum many-body physics, leading to the emergence of novel quantum phases in strongly correlated electron systems. Recent investigations have focused on the case of homogeneous site and bond dilution in the quantum antiferromagnet on the square lattice, reporting a classical geometric percolation transition between magnetic order and disorder. In this study we show how inhomogeneous bond dilution leads to percolative quantum phase transitions, which we have studied extensively by quantum Monte Carlo simulations. Quantum percolation introduces a new class of two-dimensional spin liquids, characterized by an infinite percolating network with vanishing antiferromagnetic order parameter.

preprint2002arXiv

Quantum Antiferromagnetism in Quasicrystals

The antiferromagnetic Heisenberg model is studied on a two-dimensional bipartite quasiperiodic lattice. The distribution of local staggered magnetic moments is determined on finite square approximants with up to 1393 sites, using the Stochastic Series Expansion Quantum Monte Carlo method. A non-trivial inhomogeneous ground state is found. For a given local coordination number, the values of the magnetic moments are spread out, reflecting the fact that no two sites in a quasicrystal are identical. A hierarchical structure in the values of the moments is observed which arises from the self-similarity of the quasiperiodic lattice. Furthermore, the computed spin structure factor shows antiferromagnetic modulations that can be measured in neutron scattering and nuclear magnetic resonance experiments. This generic model is a first step towards understanding magnetic quasicrystals such as the recently discovered Zn-Mg-Ho icosahedral structure.

preprint2000arXiv

Temperature dependence of impurity bound states in d-wave superconductors

We study the evolution with temperature of quasiparticle bound states around non-magnetic impurities in d-wave superconductors. The associated local density of states has a fourfold symmetry which has recently been observed in Zn-doped Bi2212 using scanning tunneling microscopy (STM). From the corresponding Bogoliubov-de Gennes equation we find that with increasing temperature the magnitude of the bound state energy increases and the amplitude of the fourfold contribution in the spinor wave functions decreases. In the pseudogap regime above T_c the fourfold angular dependence of the local tunneling conductance persists as long as the superconducting fluctuations are sufficiently strong to support a finite local order parameter. Once the gap function vanishes completely, the angular structure of the bound state wave function becomes featureless. These effects should be observable in STM studies of impurity doped high-temperature superconductors.

preprint1999arXiv

Phase diagram of three-leg ladders at strong coupling along the rungs

A phase diagram of the t-J three-leg ladder as a function of hole dopping is derived in the limit where the coupling parameters along the rungs, $t_{\perp}$ and $J_{\perp}$, are taken to be much larger than those along the legs, $t_{||}$ and $J_{||}$ At large exchange coupling along the rungs, $J_{\perp}/t_{\perp}> 3 / \sqrt{2}$, there is a transition from a low-dopping Luttinger liquid phase into a Luther-Emery liquid at a critical hole concentration $n_{crit}\approx 1/3$. In the opposite case, $J_{\perp}/t_{\perp}< 3 / \sqrt{2}$, there as a sequence of three Luttinger liquid phases (LLI, LLII and LLIII) as a function of hole dopping.

preprint1997arXiv

Hole Doping Evolution of the Quasiparticle Band in Models of Strongly Correlated Electrons for the High-T_c Cuprates

Quantum Monte Carlo (QMC) and Maximum Entropy (ME) techniques are used to study the spectral function $A({\bf p},ω)$ of the one band Hubbard model in strong coupling including a next-nearest-neighbor electronic hopping with amplitude $t'/t= -0.35$. These values of parameters are chosen to improve the comparison of the Hubbard model with angle-resolved photoemission (ARPES) data for $Sr_2 Cu O_2 Cl_2$. A narrow quasiparticle (q.p.) band is observed in the QMC analysis at the temperature of the simulation $T=t/3$, both at and away from half-filling. Such a narrow band produces a large accumulation of weight in the density of states at the top of the valence band. As the electronic density $< n >$ decreases further away from half-filling, the chemical potential travels through this energy window with a large number of states, and by $< n > \sim 0.70$ it has crossed it entirely. The region near momentum $(0,π)$ and $(π,0)$ in the spectral function is more sensitive to doping than momenta along the diagonal from $(0,0)$ to $(π,π)$. The evolution with hole density of the quasiparticle dispersion contains some of the features observed in recent ARPES data in the underdoped regime. For sufficiently large hole densities the ``flat'' bands at $(π,0)$ cross the Fermi energy, a prediction that could be tested with ARPES techniques applied to overdoped cuprates. The population of the q.p. band introduces a {\it hidden} density in the system which produces interesting consequences when the quasiparticles are assumed to interact through antiferromagnetic fluctuations and studied with the BCS gap equation formalism. In particular, a region of extended s-wave is found to compete with d-wave in the overdoped regime, i.e. when the chemical potential has almost entirely crossed the q.p.

preprint1996arXiv

Photoemission Spectra in t-J Ladders with Two Legs

Photoemission spectra for the isotropic two-leg t-J ladder are calculated at various hole-doping levels using exact diagonalization techniques. Low-energy sharp features caused by short-range antiferromagnetic correlations are observed at finite doping levels close to half-filling, above the naive Fermi momentum. These features should be observable in angle-resolved photoemission experiments. In addition, the formation of a d-wave pairing condensate as the ratio J/t is increased leads to dynamically generated spectral weight for momenta close to $k_F$ where the $d_{x^2-y^2 }$-order parameter is large.

preprint1995arXiv

Magnetic Raman Scattering in Two-Dimensional Spin-1/2 Heisenberg Antiferromagnets: Spectral Shape Anomaly and Magnetostrictive Effects

We calculate the Raman spectrum of the two-dimensional (2D) spin-1/2 Heisenberg antiferromagnet by exact diagonalization and quantum Monte Carlo techniques on clusters of up to 144 sites and, on a 16-site cluster, by considering the phonon-magnon interaction which leads to random fluctuations of the exchange integral. Results are in good agreement with experiments on various high-T_c precursors, such as La_2CuO_4 and YBa_2Cu_3O_{6.2}. In particular, our calculations reproduce the broad lineshape of the two-magnon peak, the asymmetry about its maximum, the existence of spectral weight at high energies, and the observation of nominally forbidden A_{1g} scattering.

preprint1993arXiv

Random Exchange Disorder in the Spin-1/2 XXZ Chain

The one-dimensional XXZ model is studied in the presence of disorder in the Heisenberg Exchange Integral. Recent predictions obtained from renormalization group calculations are investigated numerically using a Lanczos algorithm on chains of up to 18 sites. It is found that in the presence of strong X-Y-symmetric random exchange couplings, a ``random singlet'' phase with quasi-long-range order in the spin-spin correlations persists. As the planar anisotropy is varied, the full zero-temperature phase diagram is obtained and compared with predictions of Doty and Fisher [Phys. Rev. B {\bf 45 }, 2167 (1992)].