Catalog footprint

What is connected

110works
22topics
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

110 published item(s)

preprint2022arXiv

Automated Design of Pulse Sequences for Magnetic Resonance Fingerprinting using Physics-Inspired Optimization

Magnetic Resonance Fingerprinting (MRF) is a method to extract quantitative tissue properties such as T1 and T2 relaxation rates from arbitrary pulse sequences using conventional magnetic resonance imaging hardware. MRF pulse sequences have thousands of tunable parameters which can be chosen to maximize precision and minimize scan time. Here we perform de novo automated design of MRF pulse sequences by applying physics-inspired optimization heuristics. Our experimental data suggests systematic errors dominate over random errors in MRF scans under clinically-relevant conditions of high undersampling. Thus, in contrast to prior optimization efforts, which focused on statistical error models, we use a cost function based on explicit first-principles simulation of systematic errors arising from Fourier undersampling and phase variation. The resulting pulse sequences display features qualitatively different from previously used MRF pulse sequences and achieve fourfold shorter scan time than prior human-designed sequences of equivalent precision in T1 and T2. Furthermore, the optimization algorithm has discovered the existence of MRF pulse sequences with intrinsic robustness against shading artifacts due to phase variation.

preprint2022arXiv

Independent quality assessment of a commercial quantum random number generator

We reverse-engineer, test and analyse hardware and firmware of the commercial quantum-optical random number generator Quantis from ID Quantique. We show that > 99% of its output data originates in physically random processes: random timing of photon absorption in a semiconductor material, and random growth of avalanche owing to impact ionisation. Under a strong assumption that these processes correspond to a measurement of an initially pure state of the components, our analysis implies the unpredictability of the generated randomness. We have also found minor non-random contributions from imperfections in detector electronics and an internal processing algorithm, specific to this particular device. Our work shows that the design quality of a commercial quantum-optical randomness source can be verified without cooperation of the manufacturer and without access to the engineering documentation.

preprint2021arXiv

Prospects of Quantum Computing for Molecular Sciences

Molecular science is governed by the dynamics of electrons, atomic nuclei, and their interaction with electromagnetic fields. A reliable physicochemical understanding of these processes is crucial for the design and synthesis of chemicals and materials of economic value. Although some problems in this field are adequately addressed by classical mechanics, many require an explicit quantum mechanical description. Such quantum problems represented by exponentially large wave function should naturally benefit from quantum computation on a number of logical qubits that scales only linearly with system size. In this perspective, we focus on the potential of quantum computing for solving relevant problems in the molecular sciences -- molecular physics, chemistry, biochemistry, and materials science.

preprint2020arXiv

Efficient Quantum Walk Circuits for Metropolis-Hastings Algorithm

We present a detailed circuit implementation of Szegedy's quantization of the Metropolis-Hastings walk. This quantum walk is usually defined with respect to an oracle. We find that a direct implementation of this oracle requires costly arithmetic operations and thus reformulate the quantum walk in a way that circumvents the implementation of that specific oracle and which closely follows the classical Metropolis-Hastings walk. We also present heuristic quantum algorithms that use the quantum walk in the context of discrete optimization problems and numerically study their performances. Our numerical results indicate polynomial quantum speedups in heuristic settings.

preprint2020arXiv

Evaluation of synthetic and experimental training data in supervised machine learning applied to charge state detection of quantum dots

Automated tuning of gate-defined quantum dots is a requirement for large scale semiconductor based qubit initialisation. An essential step of these tuning procedures is charge state detection based on charge stability diagrams. Using supervised machine learning to perform this task requires a large dataset for models to train on. In order to avoid hand labelling experimental data, synthetic data has been explored as an alternative. While providing a significant increase in the size of the training dataset compared to using experimental data, using synthetic data means that classifiers are trained on data sourced from a different distribution than the experimental data that is part of the tuning process. Here we evaluate the prediction accuracy of a range of machine learning models trained on simulated and experimental data and their ability to generalise to experimental charge stability diagrams in two dimensional electron gas and nanowire devices. We find that classifiers perform best on either purely experimental or a combination of synthetic and experimental training data, and that adding common experimental noise signatures to the synthetic data does not dramatically improve the classification accuracy. These results suggest that experimental training data as well as realistic quantum dot simulations and noise models are essential in charge state detection using supervised machine learning.

preprint2020arXiv

Topological Thouless Pumping of Ultracold Fermions

A gas of electrons in a one-dimensional periodic potential can be transported even in the absence of a voltage bias if the potential is modulated slowly and periodically in time. Remarkably, the transferred charge per cycle is only sensitive to the topology of the path in parameter space. Although this so-called Thouless charge pump has first been proposed more than thirty years ago, it has not yet been realized. Here we report the first demonstration of topological Thouless pumping using ultracold atoms in a dynamically controlled optical superlattice. We observe a shift of the atomic cloud as a result of pumping and extract the topological invariance of the pumping process from this shift. We demonstrate the topological nature of the Thouless pump by varying the topology of the pumping path and verify that the topological pump indeed works in the quantum region by varying speed and temperature.

preprint2016arXiv

Accelerated nuclear quantum effects sampling with open path integrals

We numericaly demonstrate that, in double well models, the autocorrelation time of open path integral Monte Carlo simulations can be much smaller compared to standard ones using ring polymers. We also provide an intuitive explanation based on the role of instantons as transition states of the path integral pseudodynamics. Therefore we propose that, in all cases when the ground state approximation to the finite temperature partition function holds, open path integral simulations can be used to accelerate the sampling in realistic simulations aimed to explore nuclear quantum effects.

preprint2016arXiv

Bind: a Partitioned Global Workflow Parallel Programming Model

High Performance Computing is notorious for its long and expensive software development cycle. To address this challenge, we present Bind: a "partitioned global workflow" parallel programming model for C++ applications that enables quick prototyping and agile development cycles for high performance computing software targeting heterogeneous distributed many-core architectures. We present applications of Bind to Linear Algebra and MapReduce algorithms alongside with performance results.

preprint2016arXiv

Density redistribution effects in fermionic optical lattices

We simulate a one dimensional fermionic optical lattice to analyse heating due to non-adiabatic lattice loading. Our simulations reveal that, similar to the bosonic case, density redistribution effects are the major cause of heating in harmonic traps. We suggest protocols to modulate the local density distribution during the process of lattice loading, in order to reduce the excess energy. Our numerical results confirm that linear interpolation of the trapping potential and/or the interaction strength is an efficient method of doing so, bearing practical applications relevant to experiments.

preprint2016arXiv

Experimental Signatures of the Inverted Phase in InAs/GaSb Coupled Quantum Wells

Transport measurements are performed on InAs/GaSb double quantum wells at zero and finite magnetic fields applied parallel and perpendicular to the quantum wells. We investigate a sample in the inverted regime where electrons and holes coexist, and compare it with another sample in the non-inverted semiconducting regime. Activated behavior in conjunction with a strong suppression of the resistance peak at the charge neutrality point in a parallel magnetic field attest to the topological hybridization gap between electron and hole bands in the inverted sample. We observe an unconventional Landau level spectrum with energy gaps modulated by the magnetic field applied perpendicular to the quantum wells. This is caused by strong spin-orbit interaction provided jointly by the InAs and the GaSb quantum wells.

preprint2016arXiv

First order topological phase transition of the Haldane--Hubbard model

We study the interplay of topological band structure and conventional magnetic long-range order in spinful Haldane model with onsite repulsive interaction. Using the dynamical cluster approximation with clusters of up to 24 sites we find evidence of a first order phase transition from a Chern insulator at weak coupling to a topologically trivial antiferromagnetic insulator at strong coupling. These results call into question a previously found intermediate state with coexisting topological character and antiferromagnetic long-range order. Experimentally measurable signatures of the first order transition include hysteretic behavior of the double occupancy, single-particle excitation gap and nearest neighbor spin-spin correlations. This first order transition is contrasted with a continuous phase transition from the conventional band insulator to the antiferromagnetic insulator in the ionic Hubbard model on the honeycomb lattice.

preprint2016arXiv

Fulde-Ferrell-Larkin-Ovchinnikov pairing as leading instability on the square lattice

We study attractively interacting spin-1/2 fermions on the square lattice subject to a spin population imbalance. Using unbiased diagrammatic Monte Carlo simulations we find an extended region in the parameter space where the Fermi liquid is unstable towards formation of Cooper pairs with non-zero center-of-mass momentum, known as the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state. In contrast to earlier mean-field and quasi-classical studies we provide quantitative and well-controlled predictions on the existence and location of the relevant Fermi-liquid instabilities. The highest temperature where the FFLO instability can be observed is about half of the superfluid transition temperature in the unpolarized system.

preprint2016arXiv

High-Temperature Series Expansion for Spin-1/2 Heisenberg Models

We present a high-temperature series expansion code for spin-1/2 Heisenberg models on arbitrary lattices. As an example we demonstrate how to use the application for an anisotropic triangular lattice with two independent couplings J1 and J2 and calculate the high-temperature series of the magnetic susceptibility and the static structure factor up to 12th and 10th order, respectively. We show how to extract effective coupling constants for the triangular Heisenberg model from experimental data on Cs2CuBr4.

preprint2016arXiv

Local spin operators for fermion simulations

Digital quantum simulation of fermionic systems is important in the context of chemistry and physics. Simulating fermionic models on general purpose quantum computers requires imposing a fermionic algebra on spins. The previously studied Jordan-Wigner and Bravyi-Kitaev transformations are two techniques for accomplishing this task. Here we re-examine an auxiliary fermion construction which maps fermionic operators to local operators on spins. The local simulation is performed by relaxing the requirement that the number of spins should match the number of fermionic modes. Instead, auxiliary modes are introduced to enable non-consecutive fermionic couplings to be simulated with constant low-rank tensor products on spins. We connect the auxiliary fermion construction to other topological models and give examples of the construction.

preprint2016arXiv

MoTe2: A Type-II Weyl Topological Metal

Based on the ab initio calculations, we show that MoTe2, in its low-temperature orthorhombic structure characterized by an X-ray diffraction study at 100 K, realizes 4 type-II Weyl points between the N-th and N+1-th bands, where N is the total number of valence electrons per unit cell. Other WPs and nodal lines between different other bands also appear close to the Fermi level due to a complex topological band structure. We predict a series of strain-driven topological phase transitions in this compound, opening a wide range of possible experimental realizations of different topological semimetal phases. Crucially, with no strain, the number of observable surface Fermi arcs in this material is 2 - the smallest number of arcs consistent with time-reversal symmetry.

preprint2016arXiv

Optimizing spin-orbit splittings in InSb Majorana nanowires

Semiconductor-superconductor heterostructures represent a promising platform for the detection of Majorana zero modes and subsequently the processing of quantum information using their exotic non-Abelian statistics. Theoretical modeling of such low-dimensional semiconductors is generally based on phenomenological effective models. However, a more microscopic understanding of their band structure and, especially, of the spin-orbit coupling of electrons in these devices is important for optimizing their parameters for applications in quantum computing. In this paper, we approach this problem by first obtaining a highly accurate effective tight-binding model from ab initio calculations in the bulk. This model is symmetrized and correctly reproduces both the band structure and the wavefunction character. It is then used to determine Dresselhaus and Rashba spin-orbit splittings induced by finite size effects and external electric field in zincblende InSb one- and two-dimensional nanostructures as a function of growth parameters. The method presented here enables reliable simulations of realistic devices, including those used to realize exotic topological states.

preprint2016arXiv

Smooth gauge and Wannier functions for topological band structures in arbitrary dimensions

The construction of exponentially localized Wannier functions for a set of bands requires a choice of Bloch-like functions that span the space of the bands in question, and are smooth and periodic functions of k in the entire Brillouin zone. For bands with nontrivial topology, such smooth Bloch functions can only be chosen such that they do not respect the symmetries that protect the topology. This symmetry breaking is a necessary, but not sufficient condition for smoothness, and, in general, finding smooth Bloch functions for topological bands is a complicated task. We present a generic technique for finding smooth Bloch functions and constructing exponentially localized Wannier functions in the presence of nontrivial topology, given that the net Chern number of the bands in question vanishes. The technique is verified against known results in the Kane-Mele model. It is then applied to the topological insulator Bi2Se3, where the topological state is protected by two symmetries: time reversal and inversion. The resultant exponentially localized Wannier functions break both these symmetries. Finally, we illustrate how the calculation of the Chern-Simons orbital magnetoelectric response is facilitated by the proposed smooth gauge construction.

preprint2016arXiv

Stochastic series expansion simulation of the $t$-$V$ model

We present an algorithm for the efficient simulation of the half-filled spinless $t$-$V$ model on bipartite lattices, which combines the stochastic series expansion method with determinantal quantum Monte Carlo techniques widely used in fermionic simulations. The algorithm scales linearly in the inverse temperature, cubically with the system size and is free from the time-discretization error. We use it to map out the finite temperature phase diagram of the spinless $t$-$V$ model on the honeycomb lattice and observe a suppression of the critical temperature of the charge density wave phase in the vicinity of a fermionic quantum critical point.

preprint2016arXiv

Superfluidity and Density Order in a Bilayer Extended Hubbard Model

We use cluster dynamical mean field theory to study the phase diagram of the square lattice bilayer Hubbard model with an interlayer interaction. The layers are populated by two-component fermions, and the densities in both layers and the strength of the interactions are varied. We find that an attractive interlayer interaction can induce a checkerboard density ordered phase and superfluid phases, with either interlayer or intralayer pairing. Remarkably, the latter phase does not require an intralayer interaction to be present: it can be attributed to an induced attractive interaction caused by density fluctuations in the other layer.

preprint2016arXiv

The enigmatic 12/5 fractional quantum Hall effect

We numerically study the fractional quantum Hall effect at filling factors $ν=12/5$ and 13/5 (the particle-hole conjugate of 12/5) in high-quality two-dimensional GaAs heterostructures via exact diagonalization including finite well width and Landau level mixing. We find that Landau level mixing suppresses $ν=13/5$ fractional quantum Hall effect relative to $ν=12/5$. By contrast, we find both $ν=2/5$ and (its particle-hole conjugate) $ν=3/5$ fractional quantum Hall effects in the lowest Landau level to be robust under Landau level mixing and finite well-width corrections. Our results provide a possible explanation for the experimental absence of the 13/5 fractional quantum Hall state as caused by Landau level mixing effects.

preprint2016arXiv

Topological phase transitions in the repulsively interacting Haldane-Hubbard model

Using dynamical mean-field theory and exact diagonalization we study the phase diagram of the repulsive Haldane-Hubbard model, varying the interaction strength and the sublattice potential difference. In addition to the quantum Hall phase with Chern number $C=2$ and the band insulator with $C=0$ present already in the noninteracting model, the system also exhibits a $C=0$ Mott insulating phase, and a $C=1$ quantum Hall phase. We explain the latter phase by a spontaneous symmetry breaking where one of the spin-components is in the Hall state and the other in the band insulating state.

preprint2016arXiv

Topological Phases in InAs$_{1-x}$Sb$_x$: From Novel Topological Semimetal to Majorana Wire

Superconductor proximitized one-dimensional semiconductor nanowires with strong spin-orbit interaction (SOI) are at this time the most promising candidates for the realization of topological quantum information processing. In current experiments the SOI originates predominantly from extrinsic fields, induced by finite size effects and applied gate voltages. The dependence of the topological transition in these devices on microscopic details makes scaling to a large number of devices difficult unless a material with dominant intrinsic bulk SOI is used. Here we show that wires made of certain ordered alloys InAs$_{1-x}$Sb$_x$ have spin-splittings up to 20 times larger than those reached in pristine InSb wires. In particular, we show this for a stable ordered CuPt-structure at $x = 0.5$, which has an inverted band ordering and realizes a novel type of a topological semimetal with triple degeneracy points in the bulk spectrum that produce topological surface Fermi arcs. Experimentally achievable strains can drive this compound either into a topological insulator phase, or restore the normal band ordering making the CuPt-ordered InAs$_{0.5}$Sb$_{0.5}$ a semiconductor with a large intrinsic linear in $k$ bulk spin splitting.

preprint2015arXiv

Accuracy of downfolding based on the constrained random phase approximation

We study the reliability of the constrained random phase approximation (cRPA) method for the calculation of low-energy effective Hamiltonians by considering multi-orbital lattice models with one strongly correlated "target" band and two weakly correlated "screening" bands. The full multi-orbital system and the effective model are solved within dynamical mean field theory (DMFT) in a consistent way. By comparing the quasi-particle weights for the correlated bands, we examine how accurately the effective model describes the low-energy properties of the multi-band system. We show that the violation of the Pauli principle in the cRPA method leads to overscreening effects when the inter-orbital interaction is small. This problem can be overcome by using a variant of the cRPA method which restores the Pauli principle.

preprint2015arXiv

An Efficient Matrix Product Operator Representation of the Quantum-Chemical Hamiltonian

We describe how to efficiently construct the quantum chemical Hamiltonian operator in matrix product form. We present its implementation as a density matrix renormalization group (DMRG) algorithm for quantum chemical applications in a purely matrix product based framework. Existing implementations of DMRG for quantum chemistry are based on the traditional formulation of the method, which was developed from a viewpoint of Hilbert space decimation and attained a higher performance compared to straightforward implementations of matrix product based DMRG. The latter variationally optimizes a class of ansatz states known as matrix product states (MPS), where operators are correspondingly represented as matrix product operators (MPO). The MPO construction scheme presented here eliminates the previous performance disadvantages while retaining the additional flexibility provided by a matrix product approach; for example, the specification of expectation values becomes an input parameter. In this way, MPOs for different symmetries - abelian and non-abelian - and different relativistic and non-relativistic models may be solved by an otherwise unmodified program.

preprint2015arXiv

Effective models of doped quantum ladders of non-Abelian anyons

Quantum spin models have been studied extensively in one and higher dimensions. Furthermore, these systems have been doped with holes to study $t$--$J$ models of $SU(2)$ spin-1/2. Their anyonic counterparts can be built from non-Abelian anyons, such as Fibonacci anyons described by $SU(2)_3$ theories, which are quantum deformations of the $SU(2)$ algebra. Inspired by the physics of $SU(2)$ spins, several works have explored ladders of Fibonacci anyons and also one-dimensional (1D) $t$--$J$ models. Here we aim to explore the combined effects of extended dimensionality and doping by studying ladders composed of coupled chains of interacting itinerant Fibonacci anyons. We show analytically that in the limit of strong rung couplings these models can be mapped onto effective 1D models. These effective models can either be gapped models of hole pairs, or gapless models described by $t$--$J$ (or modified $t$--$J$--$V$) chains of Fibonacci anyons, whose spectrum exhibits a fractionalization into charge and anyon degrees of freedom. By using exact diagonalizations for two-leg and three-leg ladders, we show that indeed the doped ladders show exactly the same behavior as that of $t$--$J$ chains. In the strong ferromagnetic rung limit, we can obtain a new model that hosts two different kinds of Fibonacci particles - which we denote as the heavy $τ$'s and light $τ$'s. These two particle types carry the same (non-Abelian) topological charge but different (Abelian) electric charges. Once again, we map the two-dimensional ladder onto an effective chain carrying these heavy and light $τ$'s. We perform a finite size scaling analysis to show the appearance of gapless modes for certain anyon densities whereas a topological gapped phase is suggested for another density regime.

preprint2015arXiv

Efficient continuous-time quantum Monte Carlo algorithm for fermionic lattice models

Efficient continuous time quantum Monte Carlo (CT-QMC) algorithms that do not suffer from time discretization errors have become the state-of-the-art for most discrete quantum models. They have not been widely used yet for fermionic quantum lattice models, such as the Hubbard model, due to a suboptimal scaling of $O(β^3)$ with inverse temperature $β$, compared to the linear scaling of discrete time algorithms. Here we present a CT-QMC algorithms for fermionic lattice models that matches the scaling of discrete-time methods but is more efficient and free of time discretization errors. This provides an efficient simulation scheme that is free from the systematic errors opening an avenue to more precise studies of large systems at low temperatures.

preprint2015arXiv

Efficient Continuous-time Quantum Monte Carlo Method for the Ground State of Correlated Fermions

We present the ground state extension of the efficient quantum Monte Carlo algorithm for lattice fermions of arXiv:1411.0683. Based on continuous-time expansion of imaginary-time projection operator, the algorithm is free of systematic error and scales \emph{linearly} with projection time and interaction strength. Compared to the conventional quantum Monte Carlo methods for lattice fermions, this approach has greater flexibility and is easier to combine with powerful machinery such as histogram reweighting and extended ensemble simulation techniques. We discuss the implementation of the continuous-time projection in detail using the spinless $t-V$ model as an example and compare the numerical results with exact diagonalization, density-matrix-renormalization-group and infinite projected entangled-pair states calculations. Finally we use the method to study the fermionic quantum critical point of spinless fermions on a honeycomb lattice and confirm previous results concerning its critical exponents.

preprint2015arXiv

Entanglement as a resource in adiabatic quantum optimization

We explore the role of entanglement in adiabatic quantum optimization by performing approximate simulations of the real-time evolution of a quantum system while limiting the amount of entanglement. To classically simulate the time evolution of the system with a limited amount of entanglement, we represent the quantum state using matrix-product states and projected entangled-pair states. We show that the probability of finding the ground state of an Ising spin glass on either a planar or non-planar two-dimensional graph increases rapidly as the amount of entanglement in the state is increased. Furthermore, we propose evolution in complex time as a way to improve simulated adiabatic evolution and mimic the effects of thermal cooling of the quantum annealer.

preprint2015arXiv

Fate of the Kondo Effect and Impurity Quantum Phase Transitions Through the Lens of Fidelity Susceptibility

The Kondo effect is an ubiquitous phenomenon appearing at low temperature in quantum confined systems coupled to a continuous bath. Efforts in understanding and controlling it have triggered important developments across several disciplines of condensed matter physics. A recurring pattern in these studies is that the suppression of the Kondo effect often results in intriguing physical phenomena such as impurity quantum phase transitions or non-Fermi-liquid behavior. We show that the fidelity susceptibility is a sensitive indicator for such phenomena because it quantifies the sensitivity of the system's state with respect to its coupling to the bath. We demonstrate the power of fidelity susceptibility approach by using it to identify the crossover and quantum phase transitions in the one and two impurity Anderson models.

preprint2015arXiv

Fidelity susceptibility made simple: A unified quantum Monte Carlo approach

The fidelity susceptibility is a general purpose probe of phase transitions. With its origin in quantum information and in the differential geometry perspective of quantum states, the fidelity susceptibility can indicate the presence of a phase transition without prior knowledge of the local order parameter, as well as reveal the universal properties of a critical point. The wide applicability of the fidelity susceptibility to quantum many-body systems is, however, hindered by the limited computational tools to evaluate it. We present a generic, efficient, and elegant approach to compute the fidelity susceptibility of correlated fermions, bosons, and quantum spin systems in a broad range of quantum Monte Carlo methods. It can be applied both to the ground-state and non-zero temperature cases. The Monte Carlo estimator has a simple yet universal form, which can be efficiently evaluated in simulations. We demonstrate the power of this approach with applications to the Bose-Hubbard model, the spin-$1/2$ XXZ model, and use it to examine the hypothetical intermediate spin-liquid phase in the Hubbard model on the honeycomb lattice.

preprint2015arXiv

From local to global ground states in Ising spin glasses

We consider whether it is possible to find ground states of frustrated spin systems by solving them locally. Using spin glass physics and Imry-Ma arguments in addition to numerical benchmarks we quantify the power of such local solution methods and show that for the average low-dimensional spin glass problem outside the spin- glass phase the exact ground state can be found in polynomial time. In the second part we present a heuristic, general-purpose hierarchical approach which for spin glasses on chimera graphs and lattices in two and three dimensions outperforms, to our knowledge, any other solver currently around, with significantly better scaling performance than simulated annealing.

preprint2015arXiv

Heavy tails in the distribution of time-to-solution for classical and quantum annealing

For many optimization algorithms the time-to-solution depends not only on the problem size but also on the specific problem instance and may vary by many orders of magnitude. It is then necessary to investigate the full distribution and especially its tail. Here we analyze the distributions of annealing times for simulated annealing and simulated quantum annealing (by path integral quantum Monte Carlo) for random Ising spin glass instances. We find power-law distributions with very heavy tails, corresponding to extremely hard instances, but far broader distributions - and thus worse performance for hard instances - for simulated quantum annealing than for simulated annealing. Fast, non-adiabatic, annealing schedules can improve the performance of simulated quantum annealing for very hard instances by many orders of magnitude.

preprint2015arXiv

Minimizing nonadiabaticities in optical-lattice loading

In the quest to reach lower temperatures of ultra-cold gases in optical lattice experiments, non-adiabaticites during lattice loading are one of the limiting factors that prevent the same low temperatures to be reached as in experiments without lattice. Simulating the loading of a bosonic quantum gas into a one-dimensional optical lattice with and without a trap, we find that the redistribution of atomic density inside a global confining potential is by far the dominant source of heating. Based on these results we propose to adjust the trapping potential during loading to minimize changes to the density distribution. Our simulations confirm that a very simple linear interpolation of the trapping potential during loading already significantly decreases the heating of a quantum gas and we discuss how loading protocols minimizing density redistributions can be designed.

preprint2015arXiv

Negative sign problem in continuous-time quantum Monte Carlo: optimal choice of single-particle basis for impurity problems

The negative sign problem in quantum Monte Carlo (QMC) simulations of cluster impurity problems is the major bottleneck in cluster dynamical mean field calculations. In this paper we systematically investigate the dependence of the sign problem on the single-particle basis. We explore both the hybridization-expansion and the interaction-expansion variants of continuous-time QMC for three-site and four-site impurity models with baths that are diagonal in the orbital degrees of freedom. We find that the sign problem in these models can be substantially reduced by using a non-trivial single-particle basis. Such bases can be generated by diagonalizing a subset of the intracluster hoppings.

preprint2015arXiv

p-wave Superfluidity by Spin-Nematic Fermi Surface Deformation

We study attractively interacting fermions on a square lattice with dispersion relations exhibiting strong spin-dependent anisotropy. The resulting Fermi surface mismatch suppresses the s-wave BCS-type instability, clearing the way for unconventional types of order. Unbiased sampling of the Feynman diagrammatic series using Diagrammatic Monte Carlo methods reveals a rich phase diagram in the regime of intermediate coupling strength. Instead of a proposed Cooper-pair Bose metal phase [A. E. Feiguin and M. P. A. Fisher, Phys. Rev. Lett. 103, 025303 (2009)] we find an incommensurate density wave at strong anisotropy and two different p-wave superfluid states with unconventional symmetry at intermediate anisotropy.

preprint2015arXiv

Pair Correlations in Doped Hubbard Ladders

Hubbard ladders are an important stepping stone to the physics of the two-dimensional Hubbard model. While many of their properties are accessible to numerical and analytical techniques, the question of whether weakly hole-doped Hubbard ladders are dominated by superconducting or charge-density-wave correlations has so far eluded a definitive answer. In particular, previous numerical simulations of Hubbard ladders have seen a much faster decay of superconducting correlations than expected based on analytical arguments. We revisit this question using a state-of-the-art implementation of the density matrix renormalization group algorithm that allows us to simulate larger system sizes with higher accuracy than before. Performing careful extrapolations of the results, we obtain improved estimates for the Luttinger liquid parameter and the correlation functions at long distances. Our results confirm that, as suggested by analytical considerations, superconducting correlations become dominant in the limit of very small doping.

preprint2015arXiv

Partial Reinitialisation for Optimisers

Heuristic optimisers which search for an optimal configuration of variables relative to an objective function often get stuck in local optima where the algorithm is unable to find further improvement. The standard approach to circumvent this problem involves periodically restarting the algorithm from random initial configurations when no further improvement can be found. We propose a method of partial reinitialization, whereby, in an attempt to find a better solution, only sub-sets of variables are re-initialised rather than the whole configuration. Much of the information gained from previous runs is hence retained. This leads to significant improvements in the quality of the solution found in a given time for a variety of optimisation problems in machine learning.

preprint2015arXiv

Phase diagram of pyrochlore iridates: all-in-all-out magnetic ordering and non-Fermi liquid properties

We study the prototype 5d pyrochlore iridate Y2Ir2O7 from first principles using the local density approximation and dynamical mean-field theory (LDA+DMFT). We map out the phase diagram in the space of temperature, onsite Coulomb repulsion U, and filling. Consistent with experiments, we find that an all-in-all-out ordered insulating phase is stable for realistic values of U. The trigonal crystal field enhances the magnetic moments of the jeff=1/2 doublet, and strong inter-band correlations are induced by the Coulomb interaction, which indicates that a three-band description is important. We demonstrate a substantial band narrowing in the paramagnetic metallic phase and non-Fermi liquid behavior in the electron/hole doped system originating from long-lived quasi-spin moments induced by nearly flat bands.

preprint2015arXiv

Probing for quantum speedup in spin glass problems with planted solutions

The availability of quantum annealing devices with hundreds of qubits has made the experimental demonstration of a quantum speedup for optimization problems a coveted, albeit elusive goal. Going beyond earlier studies of random Ising problems, here we introduce a method to construct a set of frustrated Ising-model optimization problems with tunable hardness. We study the performance of a D-Wave Two device (DW2) with up to 503 qubits on these problems and compare it to a suite of classical algorithms, including a highly optimized algorithm designed to compete directly with the DW2. The problems are generated around predetermined ground-state configurations, called planted solutions, which makes them particularly suitable for benchmarking purposes. The problem set exhibits properties familiar from constraint satisfaction (SAT) problems, such as a peak in the typical hardness of the problems, determined by a tunable clause density parameter. We bound the hardness regime where the DW2 device either does not or might exhibit a quantum speedup for our problem set. While we do not find evidence for a speedup for the hardest and most frustrated problems in our problem set, we cannot rule out that a speedup might exist for some of the easier, less frustrated problems. Our empirical findings pertain to the specific D-Wave processor and problem set we studied and leave open the possibility that future processors might exhibit a quantum speedup on the same problem set.

preprint2015arXiv

Solving strongly correlated electron models on a quantum computer

One of the main applications of future quantum computers will be the simulation of quantum models. While the evolution of a quantum state under a Hamiltonian is straightforward (if sometimes expensive), using quantum computers to determine the ground state phase diagram of a quantum model and the properties of its phases is more involved. Using the Hubbard model as a prototypical example, we here show all the steps necessary to determine its phase diagram and ground state properties on a quantum computer. In particular, we discuss strategies for efficiently determining and preparing the ground state of the Hubbard model starting from various mean-field states with broken symmetry. We present an efficient procedure to prepare arbitrary Slater determinants as initial states and present the complete set of quantum circuits needed to evolve from these to the ground state of the Hubbard model. We show that, using efficient nesting of the various terms each time step in the evolution can be performed with just $\mathcal{O}(N)$ gates and $\mathcal{O}(\log N)$ circuit depth. We give explicit circuits to measure arbitrary local observables and static and dynamic correlation functions, both in the time and frequency domain. We further present efficient non-destructive approaches to measurement that avoid the need to re-prepare the ground state after each measurement and that quadratically reduce the measurement error.

preprint2015arXiv

Split orthogonal group: A guiding principle for sign-problem-free fermionic simulations

We present a guiding principle for designing fermionic Hamiltonians and quantum Monte Carlo (QMC) methods that are free from the infamous sign problem by exploiting the Lie groups and Lie algebras that appear naturally in the Monte Carlo weight of fermionic QMC simulations. Specifically, rigorous mathematical constraints on the determinants involving matrices that lie in the split orthogonal group provide a guideline for sign-free simulations of fermionic models on bipartite lattices. This guiding principle not only unifies the recent solutions of the sign problem based on the continuous-time quantum Monte Carlo methods and the Majorana representation, but also suggests new efficient algorithms to simulate physical systems that were previously prohibitive because of the sign problem.

preprint2015arXiv

The effect of quenched bond disorder on first-order phase transitions

We investigate the effect of quenched bond disorder on the two-dimensional three-color Ashkin-Teller model, which undergoes a first-order phase transition in the absence of impurities. This is one of the simplest and striking models in which quantitative numerical simulations can be carried out to investigate emergent criticality due to disorder rounding of first-order transition. Utilizing extensive cluster Monte Carlo simulations on large lattice sizes of up to $128 \times 128$ spins, each of which is represented by three colors taking values $\pm 1$, we show that the rounding of the first-order phase transition is an emergent criticality. We further calculate the correlation length critical exponent, $ν$, and the magnetization critical exponent, $β$, from finite size scaling analysis. We find that the critical exponents, $ν$ and $β$, change as the strength of disorder or the four-spin coupling varies, and we show that the critical exponents appear not to be in the Ising universality class. We know of no analytical approaches that can explain our non-perturbative results. However our results should inspire further work on this important problem, either numerical or analytical.

preprint2015arXiv

Thermodynamics of the Hubbard model on stacked honeycomb and square lattices

We present a numerical study of the Hubbard model on simply stacked honeycomb and square lattices, motivated by a recent experimental realization of such models with ultracold atoms in optical lattices. We perform simulations with different interlayer coupling and interaction strengths and obtain Néel transition temperatures and entropies. We provide data for the equation of state to enable comparisons of experiments and theory. We find an enhancement of the short-range correlations in the anisotropic lattices compared to the isotropic cubic lattice, in parameter regimes suitable for the interaction driven adiabatic cooling.

preprint2015arXiv

Type-II Weyl Semimetals

Fermions in nature come in several types: Dirac, Majorana and Weyl are theoretically thought to form a complete list. Even though Majorana and Weyl fermions have for decades remained experimentally elusive, condensed matter has recently emerged as fertile ground for their discovery as low energy excitations of realistic materials. Here we show the existence of yet another particle - a new type of Weyl fermion - that emerges at the boundary between electron and hole pockets in a new type of Weyl semimetal phase of matter. This fermion was missed by Weyl in 1929 due to its breaking of the stringent Lorentz symmetry of high-energy physics. Lorentz invariance however is not present in condensed matter physics, and we predict that an established material, WTe$_2$, is an example of this novel type of topological semimetal hosting the new particle as a low energy excitation around a type-2 Weyl node. This node, although still a protected crossing, has an open, finite-density of states Fermi surface, likely resulting in a plethora physical properties very different from those of standard point-like Fermi surface Weyl points.

preprint2015arXiv

Understanding Quantum Tunneling through Quantum Monte Carlo Simulations

The tunneling between the two ground states of an Ising ferromagnet is a typical example of many-body tunneling processes between two local minima, as they occur during quantum annealing. Performing quantum Monte Carlo (QMC) simulations we find that the QMC tunneling rate displays the same scaling with system size, as the rate of incoherent tunneling. The scaling in both cases is $O(Δ^2)$, where $Δ$ is the tunneling splitting. An important consequence is that QMC simulations can be used to predict the performance of a quantum annealer for tunneling through a barrier. Furthermore, by using open instead of periodic boundary conditions in imaginary time, equivalent to a projector QMC algorithm, we obtain a quadratic speedup for QMC, and achieve linear scaling in $Δ$. We provide a physical understanding of these results and their range of applicability based on an instanton picture.

preprint2014arXiv

Competing states in the t-J model: uniform d-wave state versus stripe state

Variational studies of the t-J model on the square lattice based on infinite projected-entangled pair states (iPEPS) confirm an extremely close competition between a uniform d-wave superconducting state and different stripe states. The site-centered stripe with an in-phase d-wave order has an equal or only slightly lower energy than the stripe with anti-phase d-wave order. The optimal stripe filling is not constant but increases with J/t. A nematic anisotropy reduces the pairing amplitude and the energies of stripe phases are lowered relative to the uniform state with increasing nematicity.

preprint2014arXiv

Defining and detecting quantum speedup

The development of small-scale digital and analog quantum devices raises the question of how to fairly assess and compare the computational power of classical and quantum devices, and of how to detect quantum speedup. Here we show how to define and measure quantum speedup in various scenarios, and how to avoid pitfalls that might mask or fake quantum speedup. We illustrate our discussion with data from a randomized benchmark test on a D-Wave Two device with up to 503 qubits. Comparing the performance of the device on random spin glass instances with limited precision to simulated classical and quantum annealers, we find no evidence of quantum speedup when the entire data set is considered, and obtain inconclusive results when comparing subsets of instances on an instance-by-instance basis. Our results for one particular benchmark do not rule out the possibility of speedup for other classes of problems and illustrate that quantum speedup is elusive and can depend on the question posed.

preprint2014arXiv

Fermionic Quantum Critical Point of Spinless Fermions on a Honeycomb Lattice

Spinless fermions on a honeycomb lattice provide a minimal realization of lattice Dirac fermions. Repulsive interactions between nearest neighbors drive a quantum phase transition from a Dirac semimetal to a charge-density-wave state through a fermionic quantum critical point, where the coupling of Ising order parameter to the Dirac fermions at low energy drastically affects the quantum critical behavior. Encouraged by a recently discovery of absence of the fermion sign problem in this model, we study the fermionic quantum critical point using the continuous time quantum Monte Carlo method with worm sampling technique. We estimate the transition point $V/t= 1.356(1)$ with the critical exponents $ν=0.80(3)$ and $η=0.302(7)$. Compatible results for the transition point are also obtained with infinite projected entangled-pair states.

preprint2014arXiv

Gate count estimates for performing quantum chemistry on small quantum computers

As quantum computing technology improves and quantum computers with a small but non-trivial number of N > 100 qubits appear feasible in the near future the question of possible applications of small quantum computers gains importance. One frequently mentioned application is Feynman's original proposal of simulating quantum systems, and in particular the electronic structure of molecules and materials. In this paper, we analyze the computational requirements for one of the standard algorithms to perform quantum chemistry on a quantum computer. We focus on the quantum resources required to find the ground state of a molecule twice as large as what current classical computers can solve exactly. We find that while such a problem requires about a ten-fold increase in the number of qubits over current technology, the required increase in the number of gates that can be coherently executed is many orders of magnitude larger. This suggests that for quantum computation to become useful for quantum chemistry problems, drastic algorithmic improvements will be needed.

preprint2014arXiv

Hybridization expansion Monte Carlo simulation of multi-orbital quantum impurity problems: matrix product formalism and improved Monte Carlo sampling

We explore two complementary modifications of the hybridization-expansion continuous-time Monte Carlo method, aiming at large multi-orbital quantum impurity problems. One idea is to compute the imaginary-time propagation using a matrix product states representation. We show that bond dimensions considerably smaller than the dimension of the Hilbert space are sufficient to obtain accurate results, and that this approach scales polynomially, rather than exponentially with the number of orbitals. Based on scaling analyses, we conclude that a matrix product state implementation will outperform the exact-diagonalization based method for quantum impurity problems with more than 12 orbitals. The second idea is an improved Monte Carlo sampling scheme which is applicable to all variants of the hybridization expansion method. We show that this so-called sliding window sampling scheme speeds up the simulation by at least an order of magnitude for a broad range of model parameters, with the largest improvements at low temperature.

preprint2014arXiv

Matrix Product State applications for the ALPS project

The density-matrix renormalization group method has become a standard computational approach to the low-energy physics as well as dynamics of low-dimensional quantum systems. In this paper, we present a new set of applications, available as part of the ALPS package, that provide an efficient and flexible implementation of these methods based on a matrix-product state (MPS) representation. Our applications implement, within the same framework, algorithms to variationally find the ground state and low-lying excited states as well as simulate the time evolution of arbitrary one-dimensional and two-dimensional models. Implementing the conservation of quantum numbers for generic Abelian symmetries, we achieve performance competitive with the best codes in the community. Example results are provided for (i) a model of itinerant fermions in one dimension and (ii) a model of quantum magnetism.

preprint2014arXiv

Probing the stability of the spin liquid phases in the Kitaev-Heisenberg model using tensor network algorithms

We study the extent of the spin liquid phases in the Kitaev-Heisenberg model using infinite Projected Entangled-Pair States tensor network ansatz wave functions directly in the thermodynamic limit. To assess the accuracy of the ansatz wave functions we perform benchmarks against exact results for the Kitaev model and find very good agreement for various observables. In the case of the Kitaev-Heisenberg model we confirm the existence of 6 different phases: Néel, stripy, ferromagnetic, zigzag and two spin liquid phases. We find finite extents for both spin liquid phases and discontinuous phase transitions connecting them to symmetry-broken phases.

preprint2014arXiv

Quantum versus Classical Annealing of Ising Spin Glasses

The strongest evidence for superiority of quantum annealing on spin glass problems has come from comparing simulated quantum annealing using quantum Monte Carlo (QMC) methods to simulated classical annealing [G. Santoro et al., Science 295, 2427(2002)]. Motivated by experiments on programmable quantum annealing devices we revisit the question of when quantum speedup may be expected for Ising spin glass problems. We find that even though a better scaling compared to simulated classical annealing can be achieved for QMC simulations, this advantage is due to time discretization and measurements which are not possible on a physical quantum annealing device. QMC simulations in the physically relevant continuous time limit, on the other hand, do not show superiority. Our results imply that care has to be taken when using QMC simulations to assess quantum speedup potential and are consistent with recent arguments that no quantum speedup should be expected for two-dimensional spin glass problems.

preprint2014arXiv

Real time evolution at finite temperatures with operator space matrix product states

We propose a method to simulate the real time evolution of one dimensional quantum many-body systems at finite temperature by expressing both the density matrices and the observables as matrix product states. This allows the calculation of expectation values and correlation functions as scalar products in operator space. The simulations of density matrices in inverse temperature and the local operators in the Heisenberg picture are independent and result in a grid of expectation values for all intermediate temperatures and times. Simulations can be performed using real arithmetics with only polynomial growth of computational resources in inverse temperature and time for integrable systems. The method is illustrated for the XXZ model and the single impurity Anderson model.

preprint2014arXiv

Reexamining classical and quantum models for the D-Wave One processor

We revisit the evidence for quantum annealing in the D-Wave One device (DW1) based on the study of random Ising instances. Using the probability distributions of finding the ground states of such instances, previous work found agreement with both simulated quantum annealing (SQA) and a classical rotor model. Thus the DW1 ground state success probabilities are consistent with both models, and a different measure is needed to distinguish the data and the models. Here we consider measures that account for ground state degeneracy and the distributions of excited states, and present evidence that for these new measures neither SQA nor the classical rotor model correlate perfectly with the DW1 experiments. We thus provide evidence that SQA and the classical rotor model, both of which are classically efficient algorithms, do not satisfactorily explain all the DW1 data. A complete model for the DW1 remains an open problem. Using the same criteria we find that, on the other hand, SQA and the classical rotor model correlate closely with each other. To explain this we show that the rotor model can be derived as the semiclassical limit of the spin-coherent states path integral. We also find differences in which set of ground states is found by each method, though this feature is sensitive to calibration errors of the DW1 device and to simulation parameters.

preprint2014arXiv

Renyi Entanglement Entropy of Interacting Fermions Calculated Using Continuous-Time Quantum Monte Carlo Method

We present a new algorithm for calculating the Renyi entanglement entropy of interacting fermions using the continuous-time quantum Monte Carlo method. The algorithm only samples interaction correction of the entanglement entropy, which by design ensures efficient calculation of weakly interacting systems. Combined with Monte Carlo reweighting, the algorithm also performs well for systems with strong interactions. We demonstrate the potential of this method by studying the quantum entanglement signatures of the charge-density-wave transition of interacting fermions on a square lattice.

preprint2014arXiv

Spontaneous emissions and thermalization of cold bosons in optical lattices

We study the thermalization of excitations generated by spontaneous emission events for cold bosons in an optical lattice. Computing the dynamics described by the many-body master equation, we characterize equilibration timescales in different parameter regimes. For simple observables, we find regimes in which the system relaxes rapidly to values in agreement with a thermal distribution, and others where thermalization does not occur on typical experimental timescales. Because spontaneous emissions lead effectively to a local quantum quench, this behavior is strongly dependent on the low-energy spectrum of the Hamiltonian, and undergoes a qualitative change at the Mott Insulator-superfluid transition point. These results have important implications for the understanding of thermalization after localized quenches in isolated quantum gases, as well as the characterization of heating in experiments.

preprint2014arXiv

The Phase Diagram of the $ν=5/2$ Fractional Quantum Hall Effect: Effects of Landau Level Mixing and Non-Zero Width

Interesting non-Abelian states, e.g., the Moore-Read Pfaffian and the anti-Pfaffian, offer candidate descriptions of the $ν= 5/2$ fractional quantum Hall state. But the significant controversy surrounding the nature of the $ν= 5/2$ state has been hampered by the fact that the competition between these and other states is affected by small parameter changes. To study the phase diagram of the $ν= 5/2$ state we numerically diagonalize a comprehensive effective Hamiltonian describing the fractional quantum Hall effect of electrons under realistic conditions in GaAs semiconductors. The effective Hamiltonian takes Landau level mixing into account to lowest-order perturbatively in $κ$, the ratio of the Coulomb energy scale to the cyclotron gap. We also incorporate non-zero width $w$ of the quantum well and sub-band mixing. We find the ground state in both the torus and spherical geometries as a function of $κ$ and $w$. To sort out the non-trivial competition between candidate ground states we analyze the following 4 criteria: its overlap with trial wave functions; the magnitude of energy gaps; the sign of the expectation value of an order parameter for particle-hole symmetry breaking; and the entanglement spectrum. We conclude that the ground state is in the universality class of the Moore-Read Pfaffian state, rather than the anti-Pfaffian, for $κ< {κ_c}(w)$, where ${κ_c}(w)$ is a $w$-dependent critical value $0.6 \lesssim{κ_c}(w)\lesssim 1$. We observe that both Landau level mixing and non-zero width suppress the excitation gap, but Landau level mixing has a larger effect in this regard. Our findings have important implications for the identification of non-Abelian fractional quantum Hall states.

preprint2014arXiv

The Trotter Step Size Required for Accurate Quantum Simulation of Quantum Chemistry

The simulation of molecules is a widely anticipated application of quantum computers. However, recent studies \cite{WBCH13a,HWBT14a} have cast a shadow on this hope by revealing that the complexity in gate count of such simulations increases with the number of spin orbitals $N$ as $N^8$, which becomes prohibitive even for molecules of modest size $N\sim 100$. This study was partly based on a scaling analysis of the Trotter step required for an ensemble of random artificial molecules. Here, we revisit this analysis and find instead that the scaling is closer to $N^6$ in worst case for real model molecules we have studied, indicating that the random ensemble fails to accurately capture the statistical properties of real-world molecules. Actual scaling may be significantly better than this due to averaging effects. We then present an alternative simulation scheme and show that it can sometimes outperform existing schemes, but that this possibility depends crucially on the details of the simulated molecule. We obtain further improvements using a version of the coalescing scheme of \cite{WBCH13a}; this scheme is based on using different Trotter steps for different terms. The method we use to bound the complexity of simulating a given molecule is efficient, in contrast to the approach of \cite{WBCH13a,HWBT14a} which relied on exponentially costly classical exact simulation.

preprint2014arXiv

Thermalization of strongly interacting bosons after spontaneous emissions in optical lattices

We study the out-of-equilibrium dynamics of bosonic atoms in a 1D optical lattice, after the ground-state is excited by a single spontaneous emission event, i.e. after an absorption and re-emission of a lattice photon. This is an important fundamental source of decoherence for current experiments, and understanding the resulting dynamics and changes in the many-body state is important for controlling heating in quantum simulators. Previously it was found that in the superfluid regime, simple observables relax to values that can be described by a thermal distribution on experimental time-scales, and that this breaks down for strong interactions (in the Mott insulator regime). Here we expand on this result, investigating the relaxation of the momentum distribution as a function of time, and discussing the relationship to eigenstate thermalization. For the strongly interacting limit, we provide an analytical analysis for the behavior of the system, based on an effective low-energy Hamiltonian in which the dynamics can be understood based on correlated doublon-holon pairs.

preprint2014arXiv

Topological Origin of the Fermion Sign Problem

Monte Carlo simulations are a powerful tool for elucidating the properties of complex systems across many disciplines. Not requiring any a priori knowledge, they are particularly well suited for exploring new phenomena. However, when applied to fermionic quantum systems, quantum Monte Carlo (QMC) algorithms suffer from the so-called "negative sign problem", which causes the computational effort to grow exponentially with problem size. Here we demonstrate that the fermion sign problem originates in topological properties of the configurations. In particular, we show that in the widely used auxiliary field approaches the negative sign of a configuration is a geometric phase that is the imaginary time counterpart of the Aharonov-Anandan phase, and reduces to a Berry phase in the adiabatic limit. This provides an intriguing connection between QMC simulations and classification of topological states. Our results shed clarify the controversially debated origin of the sign problem in fermionic lattice models.

preprint2014arXiv

Topological Phase Transition in the Hofstadter-Hubbard Model

We study the interplay between topological and conventional long range order of attractive fermions in a time reversal symmetric Hofstadter lattice using quantum Monte Carlo simulations, focussing on the case of one-third flux quantum per plaquette. At half-filling, the system is unstable towards s-wave pairing and charge-density-wave order at infinitesimally small interactions. At one-third-filling, the noninteracting system is a topological insulator, and a nonzero critical interaction strength is needed to drive a transition from the quantum spin Hall insulator to a superfluid. We probe the topological signature of the phase transition by threading a magnetic flux through a cylinder and observe quantized topological charge pumping.

preprint2013arXiv

A supersymmetric multicritical point in a model of lattice fermions

We study a model of spinless fermions with infinite nearest-neighbor repulsion on the square ladder which has microscopic supersymmetry. It has been conjectured that in the continuum the model is described by the superconformal minimal model with central charge c=3/2. Thus far it has not been possible to confirm this conjecture due to strong finite-size corrections in numerical data. We trace the origin of these corrections to the presence of unusual marginal operators that break Lorentz invariance, but preserve part of the supersymmetry. By relying mostly on entanglement entropy calculations with the density-matrix renormalization group, we are able to reduce finite-size effects significantly. This allows us to unambiguously determine the continuum theory of the model. We also study perturbations of the model and establish that the supersymmetric model is a multicritical point. Our work underlines the power of entanglement entropy as a probe of the phases of quantum many-body systems.

preprint2013arXiv

Anyonic quantum spin chains: Spin-1 generalizations and topological stability

There are many interesting parallels between systems of interacting non-Abelian anyons and quantum magnetism, occuring in ordinary SU(2) quantum magnets. Here we consider theories of so-called su(2)_k anyons, well-known deformations of SU(2), in which only the first k+1 angular momenta of SU(2) occur. In this manuscript, we discuss in particular anyonic generalizations of ordinary SU(2) spin chains with an emphasis on anyonic spin S=1 chains. We find that the overall phase diagrams for these anyonic spin-1 chains closely mirror the phase diagram of the ordinary bilinear-biquadratic spin-1 chain including anyonic generalizations of the Haldane phase, the AKLT construction, and supersymmetric quantum critical points. A novel feature of the anyonic spin-1 chains is an additional topological symmetry that protects the gapless phases. Distinctions further arise in the form of an even/odd effect in the deformation parameter k when considering su(2)_k anyonic theories with k>4, as well as for the special case of the su(2)_4 theory for which the spin-1 representation plays a special role. We also address anyonic generalizations of spin-1/2 chains with a focus on the topological protection provided for their gapless ground states. Finally, we put our results into context of earlier generalizations of SU(2) quantum spin chains, in particular so-called (fused) Temperley-Lieb chains.

preprint2013arXiv

Comment on: "Classical signature of quantum annealing"

In a recent preprint (arXiv:1305.4904) entitled "Classical signature of quantum annealing" Smolin and Smith point out that a bimodal distribution presented in (arXiv:1304.4595) for the success probability in the D-Wave device does not in itself provide sufficient evidence for quantum annealing, by presenting a classical model that also exhibits bimodality. Here we analyze their model and in addition present a similar model derived from the semi-classical limit of quantum spin dynamics, which also exhibits a bimodal distribution. We find that in both cases the correlations between the success probabilities of these classical models and the D-Wave device are weak compared to the correlations between a simulated quantum annealer and the D-Wave device. Indeed, the evidence for quantum annealing presented in arXiv:1304.4595 is not limited to the bimodality, but relies in addition on the success probability correlations between the D-Wave device and the simulated quantum annealer. The Smolin-Smith model and our semi-classical spin model both fail this correlation test.

preprint2013arXiv

Direct measurement of topological invariants in optical lattices

We propose an experimental technique for classifying the topology of band structures realized in optical lattices, based on a generalization of topological charge pumping in quantum Hall systems to cold atom in optical lattices. Time-of-flight measurement along one spatial direction combined with in situ detection along the transverse direction provide a direct measure of the system's Chern number, as we illustrate by calculations for the Hofstadter lattice. Based on an analogy with Wannier functions techniques of topological band theory, the method is very general and also allows the measurement of other topological invariants, such as the $Z_2$ topological invariant of time-reversal symmetric insulators.

preprint2013arXiv

Ferromagnetism of the Repulsive Atomic Fermi Gas: three-body recombination and domain formation

The simplest model for itinerant ferromagnetism, the Stoner model, has so far eluded experimental observation in repulsive ultracold fermions due to rapid three-body recombination at large scattering lengths. Here we show that a ferromagnetic phase can be stabilised by imposing a moderate optical lattice. The reduced kinetic energy drop upon formation of a polarized phase in an optical lattice extends the ferromagnetic phase to smaller scattering lengths where three-body recombination is small enough to permit experimental detection of the phase. We also show, using time dependent density functional theory, that in such a setup ferromagnetic domains emerge rapidly from a paramagnetic initial state.

preprint2013arXiv

Fibonacci topological order from quantum nets

We analyze a model of quantum nets and show it has non-abelian topological order of doubled Fibonacci type. The ground state has the same topological behavior as that of the corresponding string-net model, but our Hamiltonian can be defined on any lattice, has less complicated interactions, and its excitations are dynamical, not fixed. This Hamiltonian includes terms acting on the spins around a face, around a vertex, and special "Jones-Wenzl" terms that serve to couple long loops together. We provide strong evidence for a gap by exact diagonalization, completing the list of ingredients necessary for topological order.

preprint2013arXiv

Quantum annealing with more than one hundred qubits

Quantum technology is maturing to the point where quantum devices, such as quantum communication systems, quantum random number generators and quantum simulators, may be built with capabilities exceeding classical computers. A quantum annealer, in particular, solves hard optimisation problems by evolving a known initial configuration at non-zero temperature towards the ground state of a Hamiltonian encoding a given problem. Here, we present results from experiments on a 108 qubit D-Wave One device based on superconducting flux qubits. The strong correlations between the device and a simulated quantum annealer, in contrast with weak correlations between the device and classical annealing or classical spin dynamics, demonstrate that the device performs quantum annealing. We find additional evidence for quantum annealing in the form of small-gap avoided level crossings characterizing the hard problems. To assess the computational power of the device we compare it to optimised classical algorithms.

preprint2013arXiv

Seeing Hofstadter's Butterfly in Atomic Fermi Gases

We propose a novel way to detect the fractal energy spectrum of the Hofstadter model from the density distributions of ultracold fermions in an external trap. At low temperature, the local compressibility is proportional to the density of states of the system which reveals the fractal energy spectrum. However, thermal broadening and noises in the real experimental situation inevitably smear out fine features in the density distribution. To overcome this difficulty, we use the maximum entropy method to extract the density of states directly from the noisy thermal density distributions. Simulations show that one is able to restore the core feature of the Hofstadter's butterfly spectrum with current experimental techniques. By further reducing the noise or the temperature, one can refine the resolution and observe fine structures of the butterfly spectrum.

preprint2013arXiv

Thermodynamics and magnetic properties of the anisotropic 3D Hubbard model

We study the 3D Hubbard model with anisotropic nearest neighbor tunneling amplitudes using the dynamical cluster approximation and compare the results with a quantum simulation experiment using ultracold fermions in an optical lattice, focussing on magnetic correlations. We find that the short-range spin correlations are significantly enhanced in the direction with stronger tunneling amplitudes. Our results agree with the experimental observations and show that the experimental temperature is lower than the strong tunneling amplitude. We characterize the system by examining the spin correlations beyond neighboring sites and determine the distribution of density, entropy and spin correlation in the trapped system. We furthermore investigate the dependence of the critical entropy at the Néel transition on anisotropy.

preprint2013arXiv

Topological charge pumping in a one-dimensional optical lattice

A topological charge pump [1] transfers charge in a quantized fashion. The quantization is stable against the detailed form of the pumping protocols and external noises and shares the same topological origin as the quantum Hall effect. We propose an experiment setup to realize topological charge pumping of cold atoms in a one-dimensional optical lattice. The quantization of the pumped charge is confirmed by first-principle simulations of the dynamics of uniform and trapped systems. Quantum effects are shown to be crucial for the topological protection of the charge quantization. Finite-temperature and non-adiabatic effect on the experimental observables are discussed. Realization of such a topological charge pump servers as a firm step towards exploring topological states and non-equilibrium dynamics using cold atoms.

preprint2013arXiv

True randomness from realistic quantum devices

Even if the output of a Random Number Generator (RNG) is perfectly uniformly distributed, it may be correlated to pre-existing information and therefore be predictable. Statistical tests are thus not sufficient to guarantee that an RNG is usable for applications, e.g., in cryptography or gambling, where unpredictability is important. To enable such applications a stronger notion of randomness, termed "true randomness", is required, which includes independence from prior information. Quantum systems are particularly suitable for true randomness generation, as their unpredictability can be proved based on physical principles. Practical implementations of Quantum RNGs (QRNGs) are however always subject to noise, i.e., influences which are not fully controlled. This reduces the quality of the raw randomness generated by the device, making it necessary to post-process it. Here we provide a framework to analyse realistic QRNGs and to determine the post-processing that is necessary to turn their raw output into true randomness.

preprint2012arXiv

Bond disorder induced criticality of the three-color Ashkin-Teller model

An intriguing result of statistical mechanics is that a first-order phase transition can be rounded by disorder coupled to energy-like variables. In fact, even more intriguing is that the rounding may manifest itself as a critical point, quantum or classical. In general, it is not known, however, what universality classes, if any, such criticalities belong to. In order to shed light on this question we examine in detail the disordered three-color Ashkin-Teller model by Monte Carlo methods. Extensive analyses indicate that the critical exponents define a new universality class. We show that the rounding of the first-order transition of the pure model due to the impurities is manifested as criticality. However, the magnetization critical exponent, (β), and the correlation length critical exponent, (ν), are found to vary with disorder and the four-spin coupling strength, and we conclusively rule out that the model belongs to the universality class of the two-dimensional Ising model.

preprint2012arXiv

Double transfer through Dirac points in a tunable honeycomb optical lattice

We report on Bloch-Zener oscillations of an ultracold Fermi gas in a tunable honeycomb lattice. The quasi-momentum distribution of the atoms is measured after sequentially passing through two Dirac points. We observe a double-peak feature in the transferred fraction to the second band, both as a function of the band gap at the Dirac points and the quasi-momentum of the trajectory. Our results are in good agreement with a simple analytical model based on two successive Landau-Zener transitions. Owing to the variation of the potential gradient over the cloud size, coherent Stückelberg oscillations are not visible in our measurements. This effect of the harmonic confinement is confirmed by a numerical simulation of the dynamics of a trapped 2D system.

preprint2012arXiv

Effect of thermal fluctuations in topological p-wave superconductors

We study the effect of thermal fluctuations on the topological stability of chiral p-wave superconductors. We consider two models of superconductors: spinless and spinful with a focus on topological properties and Majorana zero-energy modes. We show that proliferation of vortex-antivortex pairs above the Kosterlitz-Thouless temperature T_KT drives the transition from a thermal Quantum Hall insulator to a thermal metal/insulator, and dramatically modifies the ground-state degeneracy splitting. Therefore, in order to utilize 2D chiral p-wave superconductors for topological quantum computing, the temperature should be much smaller than T_KT. Within the spinful chiral p-wave model, we also investigate the interplay between half-quantum vortices carrying Majorana zero-energy modes and full-quantum vortices having trivial topological charge, and discuss topological properties of half-quantum vortices in the background of proliferating full-quantum vortices.

preprint2012arXiv

From the Cooper problem to canted supersolids in Bose-Fermi mixtures

We calculate the phase diagram of the Bose-Fermi Hubbard model on the 3d cubic lattice at fermionic half filling and bosonic unit filling by means of single-site dynamical mean-field theory. For fast bosons, this is equivalent to the Cooper problem in which the bosons can induce s-wave pairing between the fermions. We also find miscible superfluid and canted supersolid phases depending on the interspecies coupling strength. In contrast, slow bosons favor fermionic charge density wave structures for attractive fermionic interactions. These competing instabilities lead to a rich phase diagram within reach of cold gas experiments.

preprint2012arXiv

One-dimensional itinerant interacting non-Abelian anyons

We construct models of interacting itinerant non-Abelian anyons moving along one-dimensional chains. We focus on itinerant Ising (Majorana) and Fibonacci anyons, which are, respectively, related to SU(2)_2 and SU(2)_3 anyons and also, respectively, describe quasiparticles of the Moore-Read and Z_3-Read-Rezayi fractional quantum Hall states. Following the derivation of the electronic large-U effective Hubbard model, we derive effective anyonic t-J models for the low-energy sectors. Solving these models by exact diagonalization, we find a fractionalization of the anyons into charge and (neutral) anyonic degrees of freedom -- a generalization of spin-charge separation of electrons which occurs in Luttinger liquids. A detailed description of the excitation spectrum can be performed by combining spectra for charge and anyonic sectors. The anyonic sector is the one of a squeezed chain of localized interacting anyons, and hence is described by the same conformal field theory (CFT), with central charge c=1/2 for Ising anyons and c=7/10 or c=4/5 for Fibonacci anyons with antiferromagnetic or ferromagnetic coupling, respectively. The charge sector is the spectrum of a chain of hardcore bosons subject to phase shifts which coincide with the momenta of the combined anyonic eigenstates, revealing a subtle coupling between charge and anyonic excitations at the microscopic level (which we also find to be present in Luttinger liquids), despite the macroscopic fractionalization. The combined central charge extracted from the entanglement entropy between segments of the chain is shown to be 1+c, where c is the central charge of the underlying CFT of the localized anyon (squeezed) chain.

preprint2012arXiv

Three-sublattice order in the SU(3) Heisenberg model on the square and triangular lattice

We present a numerical study of the SU(3) Heisenberg model of three-flavor fermions on the triangular and square lattice by means of the density-matrix renormalization group (DMRG) and infinite projected entangled-pair states (iPEPS). For the triangular lattice we confirm that the ground state has a three-sublattice order with a finite ordered moment which is compatible with the result from linear flavor wave theory (LFWT). The same type of order has recently been predicted also for the square lattice [PRL 105, 265301 (2010)] from LFWT and exact diagonalization. However, for this case the ordered moment cannot be computed based on LFWT due to divergent fluctuations. Our numerical study clearly supports this three-sublattice order, with an ordered moment of m=0.2-0.4 in the thermodynamic limit.

preprint2012arXiv

Translation invariance, topology, and protection of criticality in chains of interacting anyons

Using finite size scaling arguments, the critical properties of a chain of interacting anyons can be extracted from the low energy spectrum of a finite system. In Phys. Rev. Lett. 98, 160409 (2007), Feiguin et al. showed that an antiferromagnetic (AFM) chain of Fibonacci anyons on a torus is in the same universality class as the tricritical Ising model, and that criticality is protected by a topological symmetry. In the present paper we first review the graphical formalism for the study of anyons on the disc and demonstrate how this formalism may be consistently extended to the study of systems on surfaces of higher genus. We then employ this graphical formalism to study finite rings of interacting anyons on both the disc and the torus, and show that analysis on the disc necessarily yields an energy spectrum which is a subset of that which is obtained on the torus. For a critical Hamiltonian, one may extract from this subset the scaling dimensions of the local scaling operators which respect the topological symmetry of the system. Related considerations are also shown to apply for open chains.

preprint2011arXiv

Continuous-time Monte Carlo methods for quantum impurity models

Quantum impurity models describe an atom or molecule embedded in a host material with which it can exchange electrons. They are basic to nanoscience as representations of quantum dots and molecular conductors and play an increasingly important role in the theory of "correlated electron" materials as auxiliary problems whose solution gives the "dynamical mean field" approximation to the self energy and local correlation functions. These applications require a method of solution which provides access to both high and low energy scales and is effective for wide classes of physically realistic models. The continuous-time quantum Monte Carlo algorithms reviewed in this article meet this challenge. We present derivations and descriptions of the algorithms in enough detail to allow other workers to write their own implementations, discuss the strengths and weaknesses of the methods, summarize the problems to which the new methods have been successfully applied and outline prospects for future applications.

preprint2011arXiv

Dynamical mean-field theory for bosons

We discuss the recently developed bosonic dynamical mean-field (B-DMFT) framework, which maps a bosonic lattice model onto the selfconsistent solution of a bosonic impurity model with coupling to a reservoir of normal and condensed bosons. The effective impurity action is derived in several ways: (i) as an approximation to the kinetic energy functional of the lattice problem, (ii) using a cavity approach, and (iii) by using an effective medium approach based on adding a one-loop correction to the selfconsistently defined condensate. To solve the impurity problem, we use a continuous-time Monte Carlo algorithm based on a sampling of a perturbation expansion in the hybridization functions and the condensate wave function. As applications of the formalism we present finite temperature B-DMFT phase diagrams for the bosonic Hubbard model on a 3d cubic and 2d square lattice, the condensate order parameter as a function of chemical potential, critical exponents for the condensate, the approach to the weakly interacting Bose gas regime for weak repulsions, and the kinetic energy as a function of temperature.

preprint2011arXiv

Fractionalization of itinerant anyons in one dimensional chains

We construct models of interacting itinerant non-Abelian anyons moving along one-dimensional chains, focusing in particular on itinerant Ising anyon chains, and derive effective anyonic t-J models for the low energy sectors. Solving these models by exact diagonalization, we find a fractionalization of the anyons into charge and (non-Abelian) anyonic degrees of freedom -- a generalization of spin-charge separation of electrons which occurs in Luttinger liquids. A detailed description of the excitation spectrum by combining spectra for charge and anyonic sectors requires a subtle coupling between charge and anyonic excitations at the microscopic level (which we also find to be present in Luttinger liquids), despite the macroscopic fractionalization.

preprint2011arXiv

Galois Conjugates of Topological Phases

Galois conjugation relates unitary conformal field theories (CFTs) and topological quantum field theories (TQFTs) to their non-unitary counterparts. Here we investigate Galois conjugates of quantum double models, such as the Levin-Wen model. While these Galois conjugated Hamiltonians are typically non-Hermitian, we find that their ground state wave functions still obey a generalized version of the usual code property (local operators do not act on the ground state manifold) and hence enjoy a generalized topological protection. The key question addressed in this paper is whether such non-unitary topological phases can also appear as the ground states of Hermitian Hamiltonians. Specific attempts at constructing Hermitian Hamiltonians with these ground states lead to a loss of the code property and topological protection of the degenerate ground states. Beyond this we rigorously prove that no local change of basis can transform the ground states of the Galois conjugated doubled Fibonacci theory into the ground states of a topological model whose Hermitian Hamiltonian satisfies Lieb-Robinson bounds. These include all gapped local or quasi-local Hamiltonians. A similar statement holds for many other non-unitary TQFTs. One consequence is that the "Gaffnian" wave function cannot be the ground state of a gapped fractional quantum Hall state.

preprint2011arXiv

Multi-orbital Kondo physics of Co in Cu hosts

We investigate the electronic structure of cobalt atoms on a copper surface and in a copper host by combining density functional calculations with a numerically exact continuous-time quantum Monte Carlo treatment of the five-orbital impurity problem. In both cases we find low energy resonances in the density of states of all five Co $d$-orbitals. The corresponding self-energies indicate the formation of a Fermi liquid state at low temperatures. Our calculations yield the characteristic energy scale -- the Kondo temperature -- for both systems in good agreement with experiments. We quantify the charge fluctuations in both geometries and suggest that Co in Cu must be described by an Anderson impurity model rather than by a model assuming frozen impurity valency at low energies. We show that fluctuations of the orbital degrees of freedom are crucial for explaining the Kondo temperatures obtained in our calculations and measured in experiments.

preprint2011arXiv

Quantum spin ladders of non-Abelian anyons

Quantum ladder models, consisting of coupled chains, form intriguing systems bridging one and two dimensions and have been well studied in the context of quantum magnets and fermionic systems. Here we consider ladder systems made of more exotic quantum mechanical degrees of freedom, so-called non-Abelian anyons, which can be thought of as certain quantum deformations of ordinary SU(2) spins. Such non-Abelian anyons occur as quasiparticle excitations in topological quantum fluids, including p_x + i p_y superconductors, certain fractional quantum Hall states, and rotating Bose-Einstein condensates. Here we use a combination of exact diagonalization and conformal field theory to determine the phase diagrams of ladders with up to four chains. We discuss how phenomena familiar from ordinary SU(2) spin ladders are generalized in their anyonic counterparts, such as gapless and gapped phases, odd/even effects with the ladder width, and elementary `magnon' excitations. Other features are entirely due to the topological nature of the anyonic degrees of freedom.

preprint2011arXiv

Simultaneous dimerization and SU(4) symmetry breaking of 4-color fermions on the square lattice

Using infinite projected entangled pair states (iPEPS), exact diagonalization, and flavor-wave theory, we show that the SU(4) Heisenberg model undergoes a spontaneous dimerization on the square lattice, in contrast to its SU(2) and SU(3) counterparts, which develop Neel and three-sublattice stripe-like long-range order. Since the ground state of a dimer is not a singlet for SU(4) but a 6-dimensional irrep, this leaves the door open for further symmetry breaking. We provide evidence that, unlike in SU(4) ladders, where dimers pair up to form singlet plaquettes, here the SU(4) symmetry is additionally broken, leading to a gapless spectrum in spite of the broken translational symmetry.

preprint2011arXiv

Spectral properties of the three-dimensional Hubbard model

We present momentum resolved single-particle spectra for the three-dimensional Hubbard model for the paramagnetic and antiferromagnetically ordered phase obtained within the dynamical cluster approximation. The effective cluster problem is solved by continuous-time Quantum Monte Carlo simulations. The absence of a time discretization error and the ability to perform Monte Carlo measurements directly in Matsubara frequencies enable us to analytically continue the self-energies by maximum entropy, which is essential to obtain momentum resolved spectral functions for the N'eel state. We investigate the dependence on temperature and interaction strength and the effect of magnetic frustration introduced by a next-nearest neighbor hopping. One particular question we address here is the influence of the frustrating interaction on the metal insulator transition of the three-dimensional Hubbard model.

preprint2011arXiv

Topological Phases: An Expedition off Lattice

Motivated by the goal to give the simplest possible microscopic foundation for a broad class of topological phases, we study quantum mechanical lattice models where the topology of the lattice is one of the dynamical variables. However, a fluctuating geometry can remove the separation between the system size and the range of local interactions, which is important for topological protection and ultimately the stability of a topological phase. In particular, it can open the door to a pathology, which has been studied in the context of quantum gravity and goes by the name of `baby universe', Here we discuss three distinct approaches to suppressing these pathological fluctuations. We complement this discussion by applying Cheeger's theory relating the geometry of manifolds to their vibrational modes to study the spectra of Hamiltonians. In particular, we present a detailed study of the statistical properties of loop gas and string net models on fluctuating lattices, both analytically and numerically.

preprint2010arXiv

Complete-Graph Tensor Network States: A New Fermionic Wave Function Ansatz for Molecules

We present a new class of tensor network states that are specifically designed to capture the electron correlation of a molecule of arbitrary structure. In this ansatz, the electronic wave function is represented by a Complete-Graph Tensor Network (CGTN) ansatz which implements an efficient reduction of the number of variational parameters by breaking down the complexity of the high-dimensional coefficient tensor of a full-configuration-interaction (FCI) wave function. We demonstrate that CGTN states approximate ground states of molecules accurately by comparison of the CGTN and FCI expansion coefficients. The CGTN parametrization is not biased towards any reference configuration in contrast to many standard quantum chemical methods. This feature allows one to obtain accurate relative energies between CGTN states which is central to molecular physics and chemistry. We discuss the implications for quantum chemistry and focus on the spin-state problem. Our CGTN approach is applied to the energy splitting of states of different spin for methylene and the strongly correlated ozone molecule at a transition state structure. The parameters of the tensor network ansatz are variationally optimized by means of a parallel-tempering Monte Carlo algorithm.

preprint2010arXiv

Dynamical mean field solution of the Bose-Hubbard model

We present the effective action and self-consistency equations for the bosonic dynamical mean field (B-DMFT) approximation to the bosonic Hubbard model and show that it provides remarkably accurate phase diagrams and correlation functions. To solve the bosonic dynamical mean field equations we use a continuous-time Monte Carlo method for bosonic impurity models based on a diagrammatic expansion in the hybridization and condensate coupling. This method is readily generalized to bosonic mixtures, spinful bosons, and Bose-Fermi mixtures.

preprint2010arXiv

Measuring the equation of state of trapped ultracold bosonic systems in an optical lattice with in-situ density imaging

We analyze quantitatively how imaging techniques with single-site resolution allow to measure thermodynamical properties that cannot be inferred from time-of-light images for the trapped Bose-Hubbard model. If the normal state extends over a sufficiently large range, the chemical potential and the temperature can be extracted from a single shot, provided the sample is in thermodynamic equilibrium. When the normal state is too narrow, temperature is low but can still be extracted using the fluctuation-dissipation theorem over the entire trap range as long as the local density approximation remains valid, as was recently suggested by Qi Zhou and Tin-Lun Ho [arXiv:0908.3015]. However, for typical present-day experiments, the number of samples needed is of the order of 1000 in order to get the temperature at least $10 \%$ accurate, but it is possible to reduce the variance by 2 orders of magnitude if the density-density correlation length is short, which is the case for the Bose-Hubbard model. Our results provide further evidence that cold gases in an optical lattices can be viewed as quantum analog computers.

preprint2010arXiv

Microscopic models of interacting Yang-Lee anyons

Collective states of interacting non-Abelian anyons have recently been studied mostly in the context of certain fractional quantum Hall states, such as the Moore-Read state proposed to describe the physics of the quantum Hall plateau at filling fraction v = 5/2. In this manuscript, we further expand this line of research and present non-unitary generalizations of interacting anyon models. In particular, we introduce the notion of Yang-Lee anyons, discuss their relation to the so-called `Gaffnian' quantum Hall wave function, and describe an elementary model for their interactions. A one-dimensional version of this model -- a non-unitary generalization of the original golden chain model -- can be fully understood in terms of an exact algebraic solution and numerical diagonalization. We discuss the gapless theories of these chain models for general su(2)_k anyonic theories and their Galois conjugates. We further introduce and solve a one-dimensional version of the Levin-Wen model for non-unitary Yang-Lee anyons.

preprint2010arXiv

Mutual information in classical spin models

The total many-body correlations present in finite temperature classical spin systems are studied using the concept of mutual information. As opposed to zero-temperature quantum phase transitions, the total correlations are not maximal at the phase transition, but reach a maximum in the high temperature paramagnetic phase. The Shannon and Renyi mutual information in both Ising and Potts models in 2 dimensions are calculated numerically by combining matrix product states algorithms and Monte Carlo sampling techniques.

preprint2010arXiv

Optimized broad-histogram simulations for strong first-order phase transitions: Droplet transitions in the large-Q Potts model

The numerical simulation of strongly first-order phase transitions has remained a notoriously difficult problem even for classical systems due to the exponentially suppressed (thermal) equilibration in the vicinity of such a transition. In the absence of efficient update techniques, a common approach to improve equilibration in Monte Carlo simulations is to broaden the sampled statistical ensemble beyond the bimodal distribution of the canonical ensemble. Here we show how a recently developed feedback algorithm can systematically optimize such broad-histogram ensembles and significantly speed up equilibration in comparison with other extended ensemble techniques such as flat-histogram, multicanonical or Wang-Landau sampling. As a prototypical example of a strong first-order transition we simulate the two-dimensional Potts model with up to Q=250 different states on large systems. The optimized histogram develops a distinct multipeak structure, thereby resolving entropic barriers and their associated phase transitions in the phase coexistence region such as droplet nucleation and annihilation or droplet-strip transitions for systems with periodic boundary conditions. We characterize the efficiency of the optimized histogram sampling by measuring round-trip times tau(N,Q) across the phase transition for samples of size N spins. While we find power-law scaling of tau vs. N for small Q \lesssim 50 and N \lesssim 40^2, we observe a crossover to exponential scaling for larger Q. These results demonstrate that despite the ensemble optimization broad-histogram simulations cannot fully eliminate the supercritical slowing down at strongly first-order transitions.

preprint2010arXiv

Quantitative Determination of Temperature in the Approach to Magnetic Order of Ultracold Fermions in an Optical Lattice

We perform a quantitative simulation of the repulsive Fermi-Hubbard model using an ultracold gas trapped in an optical lattice. The entropy of the system is determined by comparing accurate measurements of the equilibrium double occupancy with theoretical calculations over a wide range of parameters. We demonstrate the applicability of both high-temperature series and dynamical mean-field theory to obtain quantitative agreement with the experimental data. The reliability of the entropy determination is confirmed by a comprehensive analysis of all systematic errors. In the center of the Mott insulating cloud we obtain an entropy per atom as low as 0.77kB which is about twice as large as the entropy at the Neel transition. The corresponding temperature depends on the atom number and for small fillings reaches values on the order of the tunneling energy.

preprint2010arXiv

Simulation of anyons with tensor network algorithms

Interacting systems of anyons pose a unique challenge to condensed matter simulations due to their non-trivial exchange statistics. These systems are of great interest as they have the potential for robust universal quantum computation, but numerical tools for studying them are as yet limited. We show how existing tensor network algorithms may be adapted for use with systems of anyons, and demonstrate this process for the 1-D Multi-scale Entanglement Renormalisation Ansatz (MERA). We apply the MERA to infinite chains of interacting Fibonacci anyons, computing their scaling dimensions and local scaling operators. The scaling dimensions obtained are seen to be in agreement with conformal field theory. The techniques developed are applicable to any tensor network algorithm, and the ability to adapt these ansaetze for use on anyonic systems opens the door for numerical simulation of large systems of free and interacting anyons in one and two dimensions.

preprint2010arXiv

Topological Order and Quantum Criticality

In this chapter we discuss aspects of the quantum critical behavior that occurs at a quantum phase transition separating a topological phase from a conventionally ordered one. We concentrate on a family of quantum lattice models, namely certain deformations of the toric code model, that exhibit continuous quantum phase transitions. One such deformation leads to a Lorentz-invariant transition in the 3D Ising universality class. An alternative deformation gives rise to a so-called conformal quantum critical point where equal-time correlations become conformally invariant and can be related to those of the 2D Ising model. We study the behavior of several physical observables, such as non-local operators and entanglement entropies, that can be used to characterize these quantum phase transitions. Finally, we briefly consider the role of thermal fluctuations and related phase transitions, before closing with a short overview of field theoretical descriptions of these quantum critical points.

preprint2010arXiv

Two-dimensional quantum liquids from interacting non-Abelian anyons

A set of localized, non-Abelian anyons - such as vortices in a p_x + i p_y superconductor or quasiholes in certain quantum Hall states - gives rise to a macroscopic degeneracy. Such a degeneracy is split in the presence of interactions between the anyons. Here we show that in two spatial dimensions this splitting selects a unique collective state as ground state of the interacting many-body system. This collective state can be a novel gapped quantum liquid nucleated inside the original parent liquid (of which the anyons are excitations). This physics is of relevance for any quantum Hall plateau realizing a non-Abelian quantum Hall state when moving off the center of the plateau.

preprint2008arXiv

Influence of the trap shape on the superfluid-Mott transition in ultracold atomic gases

The coexistence of superfluid and Mott insulator, due to the quadratic confinement potential in current optical lattice experiments, makes the accurate detection of the superfluid-Mott transition difficult. Studying alternative trapping potentials which are experimentally realizable and have a flatter center, we find that the transition can be better resolved, but at the cost of a more difficult tuning of the particle filling. When mapping out the phase diagram using local probes and the local density approximation we find that the smoother gradient of the parabolic trap is advantageous.

preprint2008arXiv

Temperature changes when adiabatically ramping up an optical lattice

When atoms are loaded into an optical lattice, the process of gradually turning on the lattice is almost adiabatic. In this paper we investigate how the temperature changes when going from the gapless superfluid phase to the gapped Mott phase along isentropic lines. To do so we calculate the entropy in the single-band Bose-Hubbard model for various densities, interaction strengths and temperatures in one and two dimensions for homogeneous and trapped systems. Our theory is able to reproduce the experimentally observed visibilities and therefore strongly supports that current experiments remain in the quantum regime for all considered lattice depths with low temperatures and minimal heating.

preprint2007arXiv

Absence of a structural glass phase in a monoatomic model liquid predicted to undergo an ideal glass transition

We study numerically a monodisperse model of interacting classical particles predicted to exhibit a static liquid-glass transition. Using a dynamical Monte Carlo method we show that the model does not freeze into a glassy phase at low temperatures. Instead, depending on the choice of the hard-core radius for the particles the system either collapses trivially or a polycrystalline hexagonal structure emerges.

preprint2007arXiv

Performance analysis of continuous-time solvers for quantum impurity models

Impurity solvers play an essential role in the numerical investigation of strongly correlated electrons systems within the "dynamical mean field" approximation. Recently, a new class of continuous-time solvers has been developed, based on a diagrammatic expansion of the partition function in either the interactions or the impurity-bath hybridization. We investigate the performance of these two complementary approaches and compare them to the well-established Hirsch-Fye method. The results show that the continuous-time methods, and in particular the version which expands in the hybridization, provide substantial gains in computational efficiency.

preprint2007arXiv

Quantum Monte-Carlo study of a two-species boson Hubbard model

We consider a two-species hard-core boson Hubbard model for a supersolid, where the two types of bosons represent vacancies and interstitials doped into a commensurate crystal. The on-site inter-species interaction may create bound states of vacancies and interstitials facilitating vacancy condensation at lower energies than in a single-species model, as suggested in an earlier mean field study. Here we carry out quantum Monte Carlo simulation to study possible supersolid phases of the model, corresponding to superfluid phases of the vacancies or interstitials. At low temperatures, we find three distinct superfluid phases. The extent of the phases and the nature of the phase transitions are discussed in comparison to mean-field theory.

preprint2006arXiv

Dimer-Quadrupolar Quantum Phase Transition in the Quasi-One-Dimensional Heisenberg Model with Biquadratic Interaction

The quasi-one-dimensional S=1 Heisenberg antiferromagnet with a biquadratic term is investigated at zero temperature by quantum Monte Carlo simulation. As the magnitude of the inter-chain coupling is increased, the system undergoes a phase transition from a spontaneously dimerized phase to a Néel ordered or spin nematic phase. The numerical results suggest the possibility of an unconventional second-order transition in which the symmetry group characterizing one phase is not a subgroup of the other.

preprint2004arXiv

Quantum spin chains with site dissipation

We use Monte Carlo simulations to study chains of Ising- and XY-spins with dissipation coupling to the site variables. The phase diagram and critical exponents of the dissipative Ising chain in a transverse magnetic field have been computed previously, and here we consider a universal ratio of susceptibilities. We furthermore present the phase diagram and exponents of the dissipative XY-chain, which exhibits a second order phase transition. All our results compare well with the predictions from a dissipative $ϕ^4$ field theory.

preprint1998arXiv

Quantum Monte Carlo Simulation of the Trellis Lattice Heisenberg Model for SrCu$_2$O$_3$ and CaV$_2$O$_5$

We study the spin-1/2 trellis lattice Heisenberg model, a coupled spin ladder system, both by perturbation around the dimer limit and by quantum Monte Carlo simulations. We discuss the influence of the inter-ladder coupling on the spin gap and the dispersion, and present results for the temperature dependence of the uniform susceptibility. The latter was found to be parameterized well by a mean-field type scaling ansatz. Finally we discuss fits of experimental measurements on SrCu$_2$O$_3$ and CaV$_2$O$_5$ to our results.

preprint1998arXiv

Universal critical temperature for Kosterlitz-Thouless transitions in bilayer quantum magnets

Recent experiments show that double layer quantum Hall systems may have a ground state with canted antiferromagnetic order. In the experimentally accessible vicinity of a quantum critical point, the order vanishes at a temperature T_{KT} = κH, where H is the magnetic field and κis a universal number determined by the interactions and Berry phases of the thermal excitations. We present quantum Monte Carlo simulations on a model spin system which support the universality of κand determine its numerical value. This allows experimental tests of an intrinsically quantum-mechanical universal quantity, which is not also a property of a higher dimensional classical critical point.