Source author record

Mats Granath

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

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

11 published item(s)

preprint2022arXiv

Error-rate-agnostic decoding of topological stabilizer codes

Efficient high-performance decoding of topological stabilizer codes has the potential to crucially improve the balance between logical failure rates and the number and individual error rates of the constituent qubits. High-threshold maximum-likelihood decoders require an explicit error model for Pauli errors to decode a specific syndrome, whereas lower-threshold heuristic approaches such as minimum weight matching are "error agnostic". Here we consider an intermediate approach, formulating a decoder that depends on the bias, i.e., the relative probability of phase-flip to bit-flip errors, but is agnostic to error rate. Our decoder is based on counting the number and effective weight of the most likely error chains in each equivalence class of a given syndrome. We use Metropolis-based Monte Carlo sampling to explore the space of error chains and find unique chains, that are efficiently identified using a hash table. Using the error-rate invariance the decoder can sample chains effectively at an error rate which is higher than the physical error rate and without the need for "thermalization" between chains in different equivalence classes. Applied to the surface code and the XZZX code, the decoder matches maximum-likelihood decoders for moderate code sizes or low error rates. We anticipate that, because of the compressed information content per syndrome, it can be taken full advantage of in combination with machine-learning methods to extrapolate Monte Carlo-generated data.

preprint2022arXiv

The XYZ$^2$ hexagonal stabilizer code

We consider a topological stabilizer code on a honeycomb grid, the "XYZ$^2$" code. The code is inspired by the Kitaev honeycomb model and is a simple realization of a "matching code" discussed by Wootton [J. Phys. A: Math. Theor. 48, 215302 (2015)], with a specific implementation of the boundary. It utilizes weight-six ($XYZXYZ$) plaquette stabilizers and weight-two ($XX$) link stabilizers on a planar hexagonal grid composed of $2d^2$ qubits for code distance $d$, with weight-three stabilizers at the boundary, stabilizing one logical qubit. We study the properties of the code using maximum-likelihood decoding, assuming perfect stabilizer measurements. For pure $X$, $Y$, or $Z$ noise, we can solve for the logical failure rate analytically, giving a threshold of 50%. In contrast to the rotated surface code and the XZZX code, which have code distance $d^2$ only for pure $Y$ noise, here the code distance is $2d^2$ for both pure $Z$ and pure $Y$ noise. Thresholds for noise with finite $Z$ bias are similar to the XZZX code, but with markedly lower sub-threshold logical failure rates. The code possesses distinctive syndrome properties with unidirectional pairs of plaquette defects along the three directions of the triangular lattice for isolated errors, which may be useful for efficient matching-based or other approximate decoding.

preprint2019arXiv

Deep Q-learning decoder for depolarizing noise on the toric code

We present an AI-based decoding agent for quantum error correction of depolarizing noise on the toric code. The agent is trained using deep reinforcement learning (DRL), where an artificial neural network encodes the state-action Q-values of error-correcting $X$, $Y$, and $Z$ Pauli operations, occurring with probabilities $p_x$, $p_y$, and $p_z$, respectively. By learning to take advantage of the correlations between bit-flip and phase-flip errors, the decoder outperforms the minimum-weight-perfect-matching (MWPM) algorithm, achieving higher success rate and higher error threshold for depolarizing noise ($p_z = p_x = p_y$), for code distances $d\leq 9$. The decoder trained on depolarizing noise also has close to optimal performance for uncorrelated noise and provides functional but sub-optimal decoding for biased noise ($p_z \neq p_x = p_y$). We argue that the DRL-type decoder provides a promising framework for future practical error correction of topological codes, striking a balance between on-the-fly calculations, in the form of forward evaluation of a deep Q-network, and pre-training and information storage. The complete code, as well as ready-to-use decoders (pre-trained networks), can be found in the repository https://github.com/mats-granath/toric-RL-decoder.

preprint2019arXiv

Unsupervised learning using topological data augmentation

Unsupervised machine learning is a cornerstone of artificial intelligence as it provides algorithms capable of learning tasks, such as classification of data, without explicit human assistance. We present an unsupervised deep learning protocol for finding topological indices of quantum systems. The core of the proposed scheme is a 'topological data augmentation' procedure that uses seed objects to generate ensembles of topologically equivalent data. Such data, assigned with dummy labels, can then be used to train a neural network classifier for sorting arbitrary objects into topological equivalence classes. Our protocol is explicitly illustrated on 2-band insulators in 1d and 2d, characterized by a winding number and a Chern number respectively. By using the augmentation technique also in the classification step we can achieve accuracy arbitrarily close to 100% even for objects with indices outside the training regime.

preprint2016arXiv

Diffusion and localization of relative strategy scores in the Minority Game

We study the equilibrium distribution of relative strategy scores of agents in the asymmetric phase ($α\equiv P/N\gtrsim 1$) of the basic Minority Game using sign-payoff, with $N$ agents holding two strategies over $P$ histories. We formulate a statistical model that makes use of the gauge freedom with respect to the ordering of an agent's strategies to quantify the correlation between the attendance and the distribution of strategies. The relative score $x\in\mathbb{Z}$ of the two strategies of an agent is described in terms of a one dimensional random walk with asymmetric jump probabilities, leading either to a static and asymmetric exponential distribution centered at $x=0$ for fickle agents or to diffusion with a positive or negative drift for frozen agents. In terms of scaled coordinates $x/\sqrt{N}$ and $t/N$ the distributions are uniquely given by $α$ and in quantitative agreement with direct simulations of the game. As the model avoids the reformulation in terms of a constrained minimization problem it can be used for arbitrary payoff functions with little calculational effort and provides a transparent and simple formulation of the dynamics of the basic Minority Game in the asymmetric phase.

preprint2014arXiv

Signatures of coherent electronic quasiparticles in the paramagnetic Mott insulator

We show that the Mott insulating state of the half-filled paramagnetic Hubbard model in infinite dimensions contains electronic quasiparticles with very small quasiparticle weight at the inner edge of the Hubbard bands. We use a stochastic and non-perturbative quantum impurity solver based on calculating the impurity self energy as a sample average over a representative distribution of impurity models solved by exact diagonalization. Due to the natural parallelization of the method, millions of poles are readily generated for the self energy which allows to work with very small pole-broadening $η$. Solutions at small and large $η$ are qualitatively different, and only at $η\leq 0.001$ (in units of half bare band width) are quasiparticles found. Evaluated on the imaginary frequency axis we find that the small $η$ solutions agree within statistical error with results using continuous time quantum Monte Carlo.

preprint2013arXiv

Bilayer graphene spectral function in RPA and self-consistent GW

We calculate the single-particle spectral function for doped bilayer graphene in the low energy limit, described by two parabolic bands with zero band gap and long range Coulomb interaction. Calculations are done using thermal Green's functions in both the random phase approximation (RPA) and the fully self-consistent GW approximation. RPA (in line with previous studies) yields a spectral function which apart from the Landau quasiparticle peaks shows additional coherent features interpreted as plasmarons, i.e. composite electron-plasmon excitations. In GW the plasmaron becomes incoherent and peaks are replaced by much broader features. The deviation of the quasiparticle weight and mass renormalization from their non-interacting values is small which indicates that bilayer graphene is a weakly interacting system. The electron energy loss function, $Im[-ε^{-1}_q(ω)]$ shows a sharp plasmon mode in RPA which in GW approximation becomes less coherent and thus consistent with the weaker plasmaron features in the corresponding single-particle spectral function.

preprint2012arXiv

Distributional exact diagonalization formalism for quantum impurity models

We develop a method for calculating the self-energy of a quantum impurity coupled to a continuous bath by stochastically generating a distribution of finite Anderson models that are solved by exact diagonalization, using the noninteracting local spectral function as a probability distribution for the sampling. The method enables calculation of the full analytic self-energy and single-particle Green's function in the complex frequency plane, without analytic continuation, and can be used for both finite and zero temperature at arbitrary fillings. Results are in good agreement with imaginary frequency data from continuous-time quantum Monte Carlo calculations for the single impurity Anderson model and the two-orbital Hubbard model within dynamical mean field theory (DMFT) as well as real frequency data for self energy of the single band Hubbard model within DMFT using numerical renormalization group. The method should be applicable to a wide range of quantum impurity models and particularly useful when high-precision real frequency results are sought.

preprint2011arXiv

Discretized Thermal Green's Functions

We present a spectral weight conserving formalism for Fermionic thermal Green's functions that are discretized in imaginary time and thus periodic in imaginary ("Matsubara") frequency. The formalism requires a generalization of the Dyson equation and the Baym-Kadanoff-Luttinger-Ward functional for the free energy. A conformal transformation is used to analytically continue the periodized Matsubara Green's function to the continuous real axis in a way that conserves the discontinuity at t=0 of the corresponding real-time Green's function. For given discretization the method allows numerical Green's function calculations of very high precision and it appears to give a well controlled convergent approximation as we decrease the discretization interval. The ideas are tested on Dynamical Mean Field Theory calculations of the paramagnetic Hubbard model.

preprint2011arXiv

Lattice expansion from isotope substitution in the Iron based superconductors

We estimate the changes in lattice parameters due to iron isotope substitution in the iron-based high temperature superconductors using a variational calculation based on the anharmonicity of the Fe-As(Se) bond. For BaFe2As2 and 57Fe to 54$Fe substitution we find a c-axis expansion of 1 10^{-3}Å and in-plane a-axis contraction of 1 10^{-5}Å at 250K in good agreement with experimental values on (Ba,K)Fe2As2 for which a negative isotope exponent α=-\frac{d\ln T_c}{d\ln m}<0 was found [P.M. Shirage et al., Phys. Rev. Lett. 103, 257003 (2009)]. For FeSe in contrast we find an expansion for both a and c lattice parameters. We discuss the relevance of these isotope induced changes of lattice parameters to the isotope effect in light of the well established sensitivity to the Fe-As(Se) coordination.

preprint2010arXiv

The Dynamical Mean Field Theory phase space extension and critical properties of the finite temperature Mott transition

We consider the finite temperature metal-insulator transition in the half filled paramagnetic Hubbard model on the infinite dimensional Bethe lattice. A new method for calculating the Dynamical Mean Field Theory fixpoint surface in the phase diagram is presented and shown to be free from the convergence problems of standard forward recursion. The fixpoint equation is then analyzed using dynamical systems methods. On the fixpoint surface the eigenspectra of its Jacobian is used to characterize the hysteresis boundaries of the first order transition line and its second order critical end point. The critical point is shown to be a cusp catastrophe in the parameter space, opening a pitchfork bifurcation along the first order transition line, while the hysteresis boundaries are shown to be saddle-node bifurcations of two merging fixpoints. Using Landau theory the properties of the critical end point is determined and related to the critical eigenmode of the Jacobian. Our findings provide new insights into basic properties of this intensively studied transition.