Showing 1–50 of 156 results for author: Martin-Delgado, M A

Implementing the Grover Algorithm in Homomorphic Encryption Schemes

Authors: Pablo Fernández, Miguel A. Martin-Delgado

Abstract: We apply quantum homomorphic encryption (QHE) schemes suitable for circuits with a polynomial number of $T/T^{\dagger}$-gates to Grover's algorithm, performing a simulation in Qiskit of a Grover circuit that contains 3 qubits. The $T/T^{\dagger}$ gate complexity of Grover's algorithm is also analysed in order to show that any Grover circuit can be evaluated homomorphically in an efficient manner.… ▽ More We apply quantum homomorphic encryption (QHE) schemes suitable for circuits with a polynomial number of $T/T^{\dagger}$-gates to Grover's algorithm, performing a simulation in Qiskit of a Grover circuit that contains 3 qubits. The $T/T^{\dagger}$ gate complexity of Grover's algorithm is also analysed in order to show that any Grover circuit can be evaluated homomorphically in an efficient manner. We discuss how to apply these QHE schemes to allow for the efficient homomorphic evaluation of any Grover circuit composed of $n$ qubits using $n-2$ extra ancilla qubits. We also show how the homomorphic evaluation of the special case where there is only one marked item can be implemented using an algorithm that makes the decryption process more efficient compared to the standard Grover algorithm. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Quantum Bayesian Inference with Renormalization for Gravitational Waves

Authors: Gabriel Escrig, Roberto Campos, Hong Qi, M. A. Martin-Delgado

Abstract: Advancements in gravitational-wave interferometers, particularly the next generation, are poised to profoundly impact gravitational wave astronomy and multimessenger astrophysics. A hybrid quantum algorithm is proposed to carry out quantum inference of parameters from compact binary coalescences detected in gravitational-wave interferometers. It performs quantum Bayesian Inference with Renormaliza… ▽ More Advancements in gravitational-wave interferometers, particularly the next generation, are poised to profoundly impact gravitational wave astronomy and multimessenger astrophysics. A hybrid quantum algorithm is proposed to carry out quantum inference of parameters from compact binary coalescences detected in gravitational-wave interferometers. It performs quantum Bayesian Inference with Renormalization and Downsampling (qBIRD). We choose binary black hole (BBH) mergers from LIGO observatories as the first case to test the algorithm, but its application can be extended to more general instances. The quantum algorithm is able to generate corner plots of relevant parameters such as chirp mass, mass ratio, spins, etc. by inference of simulated gravitational waves with known injected parameter values with zero noise, Gaussian noise and real data, thus recovering an accuracy equivalent to that of classical Markov Chain Monte Carlo inferences. The simulations are performed with sets of 2 and 4 parameters. These results enhance the possibilities to extend our capacity to track signals from coalescences over longer durations and at lower frequencies extending the accuracy and promptness of gravitational wave parameter estimation. △ Less

Submitted 29 February, 2024; originally announced March 2024.

Comments: (6 pages, 4 figures)

Randomized SearchRank: A Semiclassical Approach to a Quantum Search Engine

Authors: Sergio A. Ortega, Miguel A. Martin-Delgado

Abstract: The quantum SearchRank algorithm is a promising tool for a future quantum search engine based on PageRank quantization. However, this algorithm loses its functionality when the $N/M$ ratio between the network size $N$ and the number of marked nodes $M$ is sufficiently large. We propose a modification of the algorithm, replacing the underlying Szegedy quantum walk with a semiclassical walk. To main… ▽ More The quantum SearchRank algorithm is a promising tool for a future quantum search engine based on PageRank quantization. However, this algorithm loses its functionality when the $N/M$ ratio between the network size $N$ and the number of marked nodes $M$ is sufficiently large. We propose a modification of the algorithm, replacing the underlying Szegedy quantum walk with a semiclassical walk. To maintain the same time complexity as the quantum SearchRank algorithm we propose a simplification of the algorithm. This new algorithm is called Randomized SearchRank, since it corresponds to a quantum walk over a randomized mixed state. The performance of the SearchRank algorithms is first analyzed on an example network, and then statistically on a set of different networks of increasing size and different number of marked nodes. On the one hand, to test the search ability of the algorithms, it is computed how the probability of measuring the marked nodes decreases with $N/M$ for the quantum SearchRank, but remarkably it remains at a high value around $0.9$ for our semiclassical algorithms, solving the quantum SearchRank problem. The time complexity of the algorithms is also analyzed, obtaining a quadratic speedup with respect to the classical ones. On the other hand, the ranking functionality of the algorithms has been investigated, obtaining a good agreement with the classical PageRank distribution. Finally, the dependence of these algorithms on the intrinsic PageRank damping parameter has been clarified. Our results suggest that this parameter should be below a threshold so that the execution time does not increase drastically. △ Less

Submitted 3 January, 2024; originally announced January 2024.

Comments: RevTex 4.2, 14 color figures, 2 tables

Exotic Symmetry Breaking Properties of Self-Dual Fracton Spin Models

Authors: Giovanni Canossa, Lode Pollet, Miguel A. Martin-Delgado, Hao Song, Ke Liu

Abstract: Fracton codes host unconventional topological states of matter and are promising for fault-tolerant quantum computation due to their large coding space and strong resilience against decoherence and noise. In this work, we investigate the ground-state properties and phase transitions of two prototypical self-dual fracton spin models -- the tetrahedral Ising model and the fractal Ising model -- whic… ▽ More Fracton codes host unconventional topological states of matter and are promising for fault-tolerant quantum computation due to their large coding space and strong resilience against decoherence and noise. In this work, we investigate the ground-state properties and phase transitions of two prototypical self-dual fracton spin models -- the tetrahedral Ising model and the fractal Ising model -- which correspond to error-correction procedures for the representative fracton codes of type-I and type-II, the checkerboard code and the Haah's code, respectively, in the error-free limit. They are endowed with exotic symmetry-breaking properties that contrast sharply with the spontaneous breaking of global symmetries and deconfinement transition of gauge theories. To show these unconventional behaviors, which are associated with sub-dimensional symmetries, we construct and analyze the order parameters, correlators, and symmetry generators for both models. Notably, the tetrahedral Ising model acquires an extended semi-local ordering moment, while the fractal Ising model fits into a polynomial ring representation and leads to a fractal order parameter. Numerical studies combined with analytical tools show that both models experience a strong first-order phase transition with an anomalous $L^{-(D-1)}$ scaling, despite the fractal symmetry of the latter. Our work provides new understanding of sub-dimensional symmetry breaking and makes an important step for studying quantum-error-correction properties of the checkerboard and Haah's codes. △ Less

Submitted 28 March, 2024; v1 submitted 18 November, 2023; originally announced November 2023.

Comments: 13 pages, 6 figures. v2: published version

Journal ref: Phys. Rev. Research 6, 013304 (2024)

SQUWALS: A Szegedy QUantum WALks Simulator

Authors: Sergio A. Ortega, Miguel A. Martin-Delgado

Abstract: Szegedy's quantum walk is an algorithm for quantizing a general Markov chain. It has plenty of applications such as many variants of optimizations. In order to check its properties in an error-free environment, it is important to have a classical simulator. However, the current simulation algorithms require a great deal of memory due to the particular formulation of this quantum walk. In this pape… ▽ More Szegedy's quantum walk is an algorithm for quantizing a general Markov chain. It has plenty of applications such as many variants of optimizations. In order to check its properties in an error-free environment, it is important to have a classical simulator. However, the current simulation algorithms require a great deal of memory due to the particular formulation of this quantum walk. In this paper we propose a memory-saving algorithm that scales as $\mathcal{O}(N^2)$ with the size $N$ of the graph. We provide additional procedures for simulating Szegedy's quantum walk over mixed states and also the Semiclassical Szegedy walk. With these techniques we have built a classical simulator in Python called SQUWALS. We show that our simulator scales as $\mathcal{O}(N^2)$ in both time and memory resources. This package provides some high-level applications for algorithms based on Szegedy's quantum walk, as for example the quantum PageRank. △ Less

Submitted 26 July, 2023; originally announced July 2023.

Comments: RevTex 4.2, 16 pages, 9 color figures

Journal ref: Advanced Quantum Technologies, 2400022, (2024)

Discrete-time Semiclassical Szegedy Quantum Walks

Authors: Sergio A. Ortega, Miguel A. Martin-Delgado

Abstract: Quantum walks are promising tools based on classical random walks, with plenty of applications such as many variants of optimization. Here we introduce the semiclassical walks in discrete time, which are algorithms that combines classical and quantum dynamics. Specifically, a semiclassical walk can be understood as a classical walk where the transition matrix encodes the quantum evolution. We have… ▽ More Quantum walks are promising tools based on classical random walks, with plenty of applications such as many variants of optimization. Here we introduce the semiclassical walks in discrete time, which are algorithms that combines classical and quantum dynamics. Specifically, a semiclassical walk can be understood as a classical walk where the transition matrix encodes the quantum evolution. We have applied this algorithm to Szegedy's quantum walk, which can be applied to any arbitrary weighted graph. We first have solved the problem analytically on regular 1D cycles to show the performance of the semiclassical walks. Next, we have simulated our algorithm in a general inhomogeneous symmetric graph, finding that the inhomogeneity drives a symmetry breaking on the graph. Moreover, we show that this phenomenon is useful for the problem of ranking nodes in symmetric graphs, where the classical PageRank fails. We have demonstrated experimentally that the semiclassical walks can be applied on real quantum computers using the platform IBM Quantum. △ Less

Submitted 24 July, 2023; v1 submitted 31 March, 2023; originally announced March 2023.

Comments: RevTex 4.2, 16 pages, 11 color figures

Journal ref: Physica A: Statistical Mechanics and its Applications 625, 129021 (2023)

Homomorphic Encryption of the k=2 Bernstein-Vazirani Algorithm

Authors: Pablo Fernández, Miguel A. Martin-Delgado

Abstract: The nonrecursive Bernstein-Vazirani algorithm was the first quantum algorithm to show a superpolynomial improvement over the corresponding best classical algorithm. Here we define a class of circuits that solve a particular case of this problem for second-level recursion. This class of circuits simplifies the number of gates $T$ required to construct the oracle by making it grow linearly with the… ▽ More The nonrecursive Bernstein-Vazirani algorithm was the first quantum algorithm to show a superpolynomial improvement over the corresponding best classical algorithm. Here we define a class of circuits that solve a particular case of this problem for second-level recursion. This class of circuits simplifies the number of gates $T$ required to construct the oracle by making it grow linearly with the number of qubits in the problem. We find an application of this scheme to quantum homomorphic encryption (QHE) which is an important cryptographic technology useful for delegated quantum computation. It allows a remote server to perform quantum computations on encrypted quantum data, so that the server cannot know anything about the client's data. Liang developed QHE schemes with perfect security, $\mathcal{F}$-homomorphism, no interaction between server and client, and quasi-compactness bounded by $O(M)$ where M is the number of gates $T$ in the circuit. Precisely these schemes are suitable for circuits with a polynomial number of gates $T/T^{\dagger}$. Following these schemes, the simplified circuits we have constructed can be evaluated homomorphically in an efficient way. △ Less

Submitted 31 March, 2023; v1 submitted 30 March, 2023; originally announced March 2023.

Comments: Revtex file, color figures

Is the Majorana bound state a gravitationally neutral system?

Authors: Miguel A. Martin-Delgado

Abstract: An argument is provided that confronts the strong equivalence principle (SEP) of general relativity with the superposition principle (SP) of quantum mechanics. The result has implications on the possibility of synthesizing Majorana bound states in 1D systems or the possible violation of the strong equivalence principle by quantum effects. The argument is formulated by inducing a very special fine… ▽ More An argument is provided that confronts the strong equivalence principle (SEP) of general relativity with the superposition principle (SP) of quantum mechanics. The result has implications on the possibility of synthesizing Majorana bound states in 1D systems or the possible violation of the strong equivalence principle by quantum effects. The argument is formulated by inducing a very special fine tuning of coupling constants in the Kitaev Hamiltonian representing 1D topological superconductors. △ Less

Submitted 6 March, 2023; originally announced March 2023.

Comments: Revtex file. Color figures

Generalized Quantum PageRank Algorithm with Arbitrary Phase Rotations

Authors: Sergio A. Ortega, Miguel A. Martin-Delgado

Abstract: The quantization of the PageRank algorithm is a promising tool for a future quantum internet. Here we present a modification of the quantum PageRank introducing arbitrary phase rotations (APR) in the underlying Szegedy's quantum walk. We define three different APR schemes with only one phase as a degree of freedom. We have analyzed the behavior of these algorithms in a small generic graph, observi… ▽ More The quantization of the PageRank algorithm is a promising tool for a future quantum internet. Here we present a modification of the quantum PageRank introducing arbitrary phase rotations (APR) in the underlying Szegedy's quantum walk. We define three different APR schemes with only one phase as a degree of freedom. We have analyzed the behavior of these algorithms in a small generic graph, observing that a decrease of the phase reduces the standard deviation of the instantaneous PageRank, so the nodes of the network can be distinguished better. However, the algorithm takes more time to converge, so the phase can not be decreased arbitrarily. With these results we choose a concrete value for the phase to later apply the algorithm to complex scale-free graphs. In these networks, the original quantum PageRank is able to break the degeneracy of the residual nodes and detect secondary hubs that the classical algorithm suppresses. Nevertheless, not all of the detected secondary hubs are real according to the PageRank's definition. Some APR schemes can overcome this problem, restoring the degeneration of the residual nodes and highlighting the truly secondary hubs of the networks. Finally, we have studied the stability of the new algorithms. The original quantum algorithm was known to be more stable than the classical. We have found that one of our algorithms, whose PageRank distribution resembles the classical one, has a stability similar to the original quantum algorithm. △ Less

Submitted 31 January, 2023; v1 submitted 27 September, 2022; originally announced September 2022.

Comments: RevTex 4.2, 19 color figures, 4 tables

Journal ref: Physical Review Research 5, 013061 (2023)

Parameter Estimation of Gravitational Waves with a Quantum Metropolis Algorithm

Authors: Gabriel Escrig, Roberto Campos, Pablo A. M. Casares, M. A. Martin-Delgado

Abstract: After the first detection of a gravitational wave in 2015, the number of successes achieved by this innovative way of looking through the universe has not stopped growing. However, the current techniques for analyzing this type of events present a serious bottleneck due to the high computational power they require. In this article we explore how recent techniques based on quantum algorithms could… ▽ More After the first detection of a gravitational wave in 2015, the number of successes achieved by this innovative way of looking through the universe has not stopped growing. However, the current techniques for analyzing this type of events present a serious bottleneck due to the high computational power they require. In this article we explore how recent techniques based on quantum algorithms could surpass this obstacle. For this purpose, we propose a quantization of the classical algorithms used in the literature for the inference of gravitational wave parameters based on the well-known Quantum Walks technique applied to a Metropolis-Hastings algorithm. Finally, we develop a quantum environment on classical hardware, implementing a metric to compare quantum versus classical algorithms in a fair way. We further test all these developments in the real inference of several sets of parameters of all the events of the first detection period GWTC-1 and we find a polynomial advantage in the quantum algorithms, thus setting a first starting point for future algorithms. △ Less

Submitted 10 May, 2023; v1 submitted 10 August, 2022; originally announced August 2022.

Comments: RevTex 4.2, 3 color Figures and 1 Table

Journal ref: Classical and Quantum Gravity 40(4) 045001, 2023

Quantum Metropolis Solver: A Quantum Walks Approach to Optimization Problems

Authors: Roberto Campos, Pablo A M Casares, M A Martin-Delgado

Abstract: The efficient resolution of optimization problems is one of the key issues in today's industry. This task relies mainly on classical algorithms that present scalability problems and processing limitations. Quantum computing has emerged to challenge these types of problems. In this paper, we focus on the Metropolis-Hastings quantum algorithm that is based on quantum walks. We use this algorithm to… ▽ More The efficient resolution of optimization problems is one of the key issues in today's industry. This task relies mainly on classical algorithms that present scalability problems and processing limitations. Quantum computing has emerged to challenge these types of problems. In this paper, we focus on the Metropolis-Hastings quantum algorithm that is based on quantum walks. We use this algorithm to build a quantum software tool called Quantum Metropolis Solver (QMS). We validate QMS with the N-Queen problem to show a potential quantum advantage in an example that can be easily extrapolated to an Artificial Intelligence domain. We carry out different simulations to validate the performance of QMS and its configuration. △ Less

Submitted 13 July, 2022; originally announced July 2022.

Comments: RevTex 4.2, 6 color figures, 4 tables

Journal ref: Quantum Machine Intelligence 5, 2023

A General Relativistic Pendulum: Isochronous vs. Geodesic Motion

Authors: Miguel A. Martin-Delgado

Abstract: General Relativity (GR) is shown to be a complete theory with respect to the isochrony of the pendulum. This guarantees that time can be measured with a mechanical clock within the theory itself as a matter of principle. The proper and coordinate oscillation frequencies of a simple pendulum are computed as a function of 3 coordinate parameters: the distance of the fulcrum to the gravity source, it… ▽ More General Relativity (GR) is shown to be a complete theory with respect to the isochrony of the pendulum. This guarantees that time can be measured with a mechanical clock within the theory itself as a matter of principle. The proper and coordinate oscillation frequencies of a simple pendulum are computed as a function of 3 coordinate parameters: the distance of the fulcrum to the gravity source, its length and the Schwarzschild radius. Novel GR effects appear in the Schwarzschild coordinates as divergences, zeroes and other instabilities of the oscillation frequencies that are absent in Newtonian gravity. The anomalies in the proper and coordinate frequencies occur under extreme conditions: either the fulcrum or the pendulum mass or both remain inside the Event Horizon. △ Less

Submitted 5 May, 2022; originally announced May 2022.

Comments: RevTex file. Color figures

Reinforcement Learning Generation of 4-Qubits Entangled States

Authors: Sara Giordano, Miguel A. Martin-Delgado

Abstract: We have devised an artificial intelligence algorithm with machine reinforcement learning (Q-learning) to construct remarkable entangled states with 4 qubits. This way, the algorithm is able to generate representative states for some of the 49 true SLOCC classes of the four-qubit entanglement states. In particular, it is possible to reach at least one true SLOCC class for each of the nine entanglem… ▽ More We have devised an artificial intelligence algorithm with machine reinforcement learning (Q-learning) to construct remarkable entangled states with 4 qubits. This way, the algorithm is able to generate representative states for some of the 49 true SLOCC classes of the four-qubit entanglement states. In particular, it is possible to reach at least one true SLOCC class for each of the nine entanglement families. The quantum circuits synthesized by the algorithm may be useful for the experimental realization of these important classes of entangled states and to draw conclusions about the intrinsic properties of our universe. We introduce a graphical tool called the state-link graph (SLG) to represent the construction of the Quality matrix (Q-matrix) used by the algorithm to build a given objective state belonging to the corresponding entanglement class. This allows us to discover the necessary connections between specific entanglement features and the role of certain quantum gates that the algorithm needs to include in the quantum gate set of actions. The quantum circuits found are optimal by construction with respect to the quantum gate-set chosen. These SLGs make the algorithm simple, intuitive and a useful resource for the automated construction of entangled states with a low number of qubits. △ Less

Submitted 27 October, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Comments: Tex file, color figures

Journal ref: Physical Review Research, Vol. 4, No. 4, 1 October 2022 issue

Simulating key properties of lithium-ion batteries with a fault-tolerant quantum computer

Authors: Alain Delgado, Pablo A. M. Casares, Roberto dos Reis, Modjtaba Shokrian Zini, Roberto Campos, Norge Cruz-Hernández, Arne-Christian Voigt, Angus Lowe, Soran Jahangiri, M. A. Martin-Delgado, Jonathan E. Mueller, Juan Miguel Arrazola

Abstract: There is a pressing need to develop new rechargeable battery technologies that can offer higher energy storage, faster charging, and lower costs. Despite the success of existing methods for the simulation of battery materials, they can sometimes fall short of delivering accurate and reliable results. Quantum computing has been discussed as an avenue to overcome these issues, but only limited work… ▽ More There is a pressing need to develop new rechargeable battery technologies that can offer higher energy storage, faster charging, and lower costs. Despite the success of existing methods for the simulation of battery materials, they can sometimes fall short of delivering accurate and reliable results. Quantum computing has been discussed as an avenue to overcome these issues, but only limited work has been done to outline how they may impact battery simulations. In this work, we provide a detailed answer to the following question: how can a quantum computer be used to simulate key properties of a lithium-ion battery? Based on recently-introduced first-quantization techniques, we lay out an end-to-end quantum algorithm for calculating equilibrium cell voltages, ionic mobility, and thermal stability. These can be obtained from ground-state energies of materials, which is the core calculation executed by the quantum computer using qubitization-based quantum phase estimation. The algorithm includes explicit methods for preparing approximate ground states of periodic materials in first quantization. We bring these insights together to perform the first estimation of the resources required to implement a quantum algorithm for simulating a realistic cathode material, dilithium iron silicate. △ Less

Submitted 6 February, 2023; v1 submitted 25 April, 2022; originally announced April 2022.

Comments: 31 pages, 14 figures

Journal ref: Phys. Rev. A 106, 032428 (2022)

Robust nonequilibrium surface currents in the 3D Hofstadter model

Authors: Mark T. Mitchison, Ángel Rivas, Miguel A. Martin-Delgado

Abstract: Genuinely two-dimensional robust crosscurrents -- which flow against the natural direction of heat flux -- have been missing since the discovery of their one-dimensional counterpart. We provide a setup to realize them on a cubic three-dimensional (3D) lattice hosting a Hofstadter model coupled to two heat baths with different temperatures. We show that these currents exhibit dissipative robustness… ▽ More Genuinely two-dimensional robust crosscurrents -- which flow against the natural direction of heat flux -- have been missing since the discovery of their one-dimensional counterpart. We provide a setup to realize them on a cubic three-dimensional (3D) lattice hosting a Hofstadter model coupled to two heat baths with different temperatures. We show that these currents exhibit dissipative robustness: they are stable against the presence of impurities and tilting of the gauge field in certain nonequilibrium configurations. Moreover, we find protected boundary currents with genuinely 3D robustness, i.e. they are only stable if tunnelling can occur in all three spatial directions. The model also presents generic surface currents, which are robust for both bosonic and fermionic systems. We identify the underlying qualitative mechanism responsible for the robustness of the surface currents and the crucial role played by certain discrete symmetries. △ Less

Submitted 19 October, 2022; v1 submitted 12 January, 2022; originally announced January 2022.

Comments: Comments welcome. v1: 7 pages, 5 figures. v2: Final author version

Journal ref: Phys. Rev. Research 4, 043032 (2022)

Optimal Thresholds for Fracton Codes and Random Spin Models with Subsystem Symmetry

Authors: Hao Song, Janik Schönmeier-Kromer, Ke Liu, Oscar Viyuela, Lode Pollet, M. A. Martin-Delgado

Abstract: Fracton models provide examples of novel gapped quantum phases of matter that host intrinsically immobile excitations and therefore lie beyond the conventional notion of topological order. Here, we calculate optimal error thresholds for quantum error correcting codes based on fracton models. By mapping the error-correction process for bit-flip and phase-flip noises into novel statistical models wi… ▽ More Fracton models provide examples of novel gapped quantum phases of matter that host intrinsically immobile excitations and therefore lie beyond the conventional notion of topological order. Here, we calculate optimal error thresholds for quantum error correcting codes based on fracton models. By mapping the error-correction process for bit-flip and phase-flip noises into novel statistical models with Ising variables and random multi-body couplings, we obtain models that exhibit an unconventional subsystem symmetry instead of a more usual global symmetry. We perform large-scale parallel tempering Monte Carlo simulations to obtain disorder-temperature phase diagrams, which are then used to predict optimal error thresholds for the corresponding fracton code. Remarkably, we found that the X-cube fracton code displays a minimum error threshold ($7.5\%$) that is much higher than 3D topological codes such as the toric code ($3.3\%$), or the color code ($1.9\%$). This result, together with the predicted absence of glass order at the Nishimori line, shows great potential for fracton phases to be used as quantum memory platforms. △ Less

Submitted 4 December, 2022; v1 submitted 9 December, 2021; originally announced December 2021.

Comments: 6+10 pages, 4+4 figures, Fig 1 selected for PRL cover. v2: minor changes to introduction, more quantum error correction details added to supplementary materials, typos corrected, published version

Journal ref: Phys. Rev. Lett. 129, 230502 (2022)

TFermion: A non-Clifford gate cost assessment library of quantum phase estimation algorithms for quantum chemistry

Authors: Pablo A. M. Casares, Roberto Campos, M. A. Martin-Delgado

Abstract: Quantum Phase Estimation is one of the most useful quantum computing algorithms for quantum chemistry and as such, significant effort has been devoted to designing efficient implementations. In this article, we introduce TFermion, a library designed to estimate the T-gate cost of such algorithms, for an arbitrary molecule. As examples of usage, we estimate the T-gate cost of a few simple molecules… ▽ More Quantum Phase Estimation is one of the most useful quantum computing algorithms for quantum chemistry and as such, significant effort has been devoted to designing efficient implementations. In this article, we introduce TFermion, a library designed to estimate the T-gate cost of such algorithms, for an arbitrary molecule. As examples of usage, we estimate the T-gate cost of a few simple molecules and compare the same Taylorization algorithms using Gaussian and plane-wave basis. △ Less

Submitted 13 July, 2022; v1 submitted 12 October, 2021; originally announced October 2021.

Comments: 37 pages, 9 figures, version accepted for publication in Quantum

Journal ref: Quantum 6, 768 (2022)

Robust nonequilibrium edge currents with and without band topology

Authors: Mark T. Mitchison, Ángel Rivas, Miguel A. Martin-Delgado

Abstract: We study two-dimensional bosonic and fermionic lattice systems under nonequilibrium conditions corresponding to a sharp gradient of temperature imposed by two thermal baths. In particular, we consider a lattice model with broken time-reversal symmetry that exhibits both topologically trivial and nontrivial phases. Using a nonperturbative Green function approach, we characterize the nonequilibrium… ▽ More We study two-dimensional bosonic and fermionic lattice systems under nonequilibrium conditions corresponding to a sharp gradient of temperature imposed by two thermal baths. In particular, we consider a lattice model with broken time-reversal symmetry that exhibits both topologically trivial and nontrivial phases. Using a nonperturbative Green function approach, we characterize the nonequilibrium current distribution in different parameter regimes. For both bosonic and fermionic systems, we find chiral edge currents that are robust against coupling to reservoirs and to the presence of defects on the boundary or in the bulk. This robustness not only originates from topological effects at zero temperature but, remarkably, also persists as a result of dissipative symmetries in regimes where band topology plays no role. Chirality of the edge currents implies that energy locally flows against the temperature gradient without any external work input. In the fermionic case, there is also a regime with topologically protected boundary currents, which nonetheless do not circulate around all system edges. △ Less

Submitted 14 April, 2022; v1 submitted 10 June, 2021; originally announced June 2021.

Comments: v1: 5+4 pages, 4+2 figures. Comments welcome. v2: Final author version

Journal ref: Phys. Rev. Lett. 128, 120403 (2022)

QFold: Quantum Walks and Deep Learning to Solve Protein Folding

Authors: P A M Casares, Roberto Campos, M A Martin-Delgado

Abstract: Predicting the 3D structure of proteins is one of the most important problems in current biochemical research. In this article, we explain how to combine recent deep learning advances with the well known technique of quantum walks applied to a Metropolis algorithm. The result, QFold, is a fully scalable hybrid quantum algorithm that, in contrast to previous quantum approaches, does not require a l… ▽ More Predicting the 3D structure of proteins is one of the most important problems in current biochemical research. In this article, we explain how to combine recent deep learning advances with the well known technique of quantum walks applied to a Metropolis algorithm. The result, QFold, is a fully scalable hybrid quantum algorithm that, in contrast to previous quantum approaches, does not require a lattice model simplification and instead relies on the much more realistic assumption of parameterization in terms of torsion angles of the amino acids. We compare it with its classical analog for different annealing schedules and find a polynomial quantum advantage, and implement a minimal realization of the quantum Metropolis in the IBMQ Casablanca quantum system. △ Less

Submitted 9 March, 2022; v1 submitted 25 January, 2021; originally announced January 2021.

Comments: RevTex 4.1, 9 color figures, 2 tables

Journal ref: Quantum Science and Technology 7 025013, 2022

An Exclusion Principle for Sum-Free Quantum Numbers

Authors: Miguel A. Martin-Delgado

Abstract: A hypothetical exclusion principle for quantum particles is introduced that generalizes the exclusion and inclusion principles for fermions and bosons, respectively: the correlated exclusion principle. The sum-free condition for Schur numbers can be read off as a form of exclusion principle for quantum particles. Consequences of this interpretation are analysed within the framework of quantum many… ▽ More A hypothetical exclusion principle for quantum particles is introduced that generalizes the exclusion and inclusion principles for fermions and bosons, respectively: the correlated exclusion principle. The sum-free condition for Schur numbers can be read off as a form of exclusion principle for quantum particles. Consequences of this interpretation are analysed within the framework of quantum many-body systems. A particular instance of the correlated exclusion principle can be solved explicitly yielding a sequence of quantum numbers that exhibits a fractal structure and is a relative of the Thue-Thurston sequence. The corresponding algebra of creation and annihilation operators can be identified in terms of commutation and anticommutation relations of a restricted version of the hard-core boson algebra. △ Less

Submitted 30 September, 2020; originally announced September 2020.

Comments: Revtex 4.0 file. Color figures

Chiral Prime Concatenations

Authors: Miguel A. Martin-Delgado

Abstract: The notion of chiral prime concatenations is studied as a recursive construction of prime numbers starting from a seed set and with appropriate blocks to define the primality growth, generation by generation, either from the right or from the left. Several basic questions are addressed like whether chiral concatenation is a symmetric process, an endless process, as well as the calculation of large… ▽ More The notion of chiral prime concatenations is studied as a recursive construction of prime numbers starting from a seed set and with appropriate blocks to define the primality growth, generation by generation, either from the right or from the left. Several basic questions are addressed like whether chiral concatenation is a symmetric process, an endless process, as well as the calculation of largest chiral prime numbers. In particular, the largest left-concatenated prime number is constructed. It is a unique prime number of 24 digits. By introducing anomalous left-concatenations of primes we can surpass the limit of 24 digits for left-concatenated primes. It is conjectured that prime numbers are left chiral under anomalous concatenations. △ Less

Submitted 25 September, 2020; originally announced September 2020.

Comments: Revtex 4.0 file with figures

The New SI and the Fundamental Constants of Nature

Authors: Miguel A. Martin-Delgado

Abstract: The launch in 2019 of the new international system of units is an opportunity to highlight the key role that the fundamental laws of physics and chemistry play in our lives and in all the processes of basic research, industry and commerce. The main objective of these notes is to present the new SI in an accessible way for a wide audience. After reviewing the fundamental constants of nature and its… ▽ More The launch in 2019 of the new international system of units is an opportunity to highlight the key role that the fundamental laws of physics and chemistry play in our lives and in all the processes of basic research, industry and commerce. The main objective of these notes is to present the new SI in an accessible way for a wide audience. After reviewing the fundamental constants of nature and its universal laws, the new definitions of SI units are presented using, as a unifying principle, the discrete nature of energy, matter and information in these universal laws. The new SI system is here to stay: although the experimental realizations may change due to technological improvements, the definitions will remain unaffected. Quantum metrology is expected to be one of the driving forces to achieve new quantum technologies of the second generation. ----- La puesta en marcha en 2019 del nuevo sistema internacional de unidades es una oportunidad para resaltar el papel fundamental que las leyes fundamentales de la Física y la Química juegan en nuestra vida y en todos los procesos de la investigación fundamental, la industria y el comercio. El principal objetivo de estas notas es presentar el nuevo SI de forma accesible para una audiencia amplia. Tras repasar las constantes fundamentales de la naturaleza y sus leyes universales, se presentan las nuevas definiciones de las unidades SI utilizando como principio unificador la naturaleza discreta de la energía, la materia y la información en esas leyes universales. El nuevo sistema SI tiene vocación de futuro: aunque las realizaciones experimentales cambien por mejoras tecnolgicas, las definiciones permanecerán inalteradas. La Metrología cuántica está llamada a ser uno de las fuerzas motrices para conseguir nuevas tecnologías cuánticas de segunda generación. △ Less

Submitted 24 April, 2020; originally announced May 2020.

Comments: Revtex file, color figures. English version y en español

Journal ref: European Journal of Physics, Volume 41, Number 6 (2020)

Determination of the Semion Code Threshold using Neural Decoders

Authors: Santiago Varona, Miguel Angel Martin-Delgado

Abstract: We compute the error threshold for the semion code, the companion of the Kitaev toric code with the same gauge symmetry group $\mathbb{Z}_2$. The application of statistical mechanical mapping methods is highly discouraged for the semion code, since the code is non-Pauli and non-CSS. Thus, we use machine learning methods, taking advantage of the near-optimal performance of some neural network decod… ▽ More We compute the error threshold for the semion code, the companion of the Kitaev toric code with the same gauge symmetry group $\mathbb{Z}_2$. The application of statistical mechanical mapping methods is highly discouraged for the semion code, since the code is non-Pauli and non-CSS. Thus, we use machine learning methods, taking advantage of the near-optimal performance of some neural network decoders: multilayer perceptrons and convolutional neural networks (CNNs). We find the values $p_{\text {eff}}=9.5\%$ for uncorrelated bit-flip and phase-flip noise, and $p_{\text {eff}}=10.5\%$ for depolarizing noise. We contrast these values with a similar analysis of the Kitaev toric code on a hexagonal lattice with the same methods. For convolutional neural networks, we use the ResNet architecture, which allows us to implement very deep networks and results in better performance and scalability than the multilayer perceptron approach. We analyze and compare in detail both approaches and provide a clear argument favoring the CNN as the best suited numerical method for the semion code. △ Less

Submitted 12 September, 2020; v1 submitted 20 February, 2020; originally announced February 2020.

Comments: Minor corrections

Journal ref: Phys. Rev. A 102, 032411 (2020)

A quantum active learning algorithm for sampling against adversarial attacks

Authors: P. A. M. Casares, M. A. Martin-Delgado

Abstract: Adversarial attacks represent a serious menace for learning algorithms and may compromise the security of future autonomous systems. A theorem by Khoury and Hadfield-Menell (KH), provides sufficient conditions to guarantee the robustness of machine learning algorithms, but comes with a caveat: it is crucial to know the smallest distance among the classes of the corresponding classification problem… ▽ More Adversarial attacks represent a serious menace for learning algorithms and may compromise the security of future autonomous systems. A theorem by Khoury and Hadfield-Menell (KH), provides sufficient conditions to guarantee the robustness of machine learning algorithms, but comes with a caveat: it is crucial to know the smallest distance among the classes of the corresponding classification problem. We propose a theoretical framework that allows us to think of active learning as sampling the most promising new points to be classified, so that the minimum distance between classes can be found and the theorem KH used. Additionally, we introduce a quantum active learning algorithm that makes use of such framework and whose complexity is polylogarithmic in the dimension of the space, $m$, and the size of the initial training data $n$, provided the use of qRAMs; and polynomial in the precision, achieving an exponential speedup over the equivalent classical algorithm in $n$ and $m$. This algorithm may be nevertheless `dequantized' reducing the advantage to polynomial. △ Less

Submitted 19 July, 2020; v1 submitted 6 December, 2019; originally announced December 2019.

Comments: Contains an additional dequantization appendix E that does not appear in the published version

Journal ref: New Journal of Physics, 2020

A Quantum Interior-Point Predictor-Corrector Algorithm for Linear Programming

Authors: P. A. M. Casares, M. A. Martin-Delgado

Abstract: We introduce a new quantum optimization algorithm for dense Linear Programming problems, which can be seen as the quantization of the Interior Point Predictor-Corrector algorithm \cite{Predictor-Corrector} using a Quantum Linear System Algorithm \cite{DenseHHL}. The (worst case) work complexity of our method is, up to polylogarithmic factors, $O(L\sqrt{n}(n+m)\overline{||M||_F}\barκ^2ε^{-2})$ for… ▽ More We introduce a new quantum optimization algorithm for dense Linear Programming problems, which can be seen as the quantization of the Interior Point Predictor-Corrector algorithm \cite{Predictor-Corrector} using a Quantum Linear System Algorithm \cite{DenseHHL}. The (worst case) work complexity of our method is, up to polylogarithmic factors, $O(L\sqrt{n}(n+m)\overline{||M||_F}\barκ^2ε^{-2})$ for $n$ the number of variables in the cost function, $m$ the number of constraints, $ε^{-1}$ the target precision, $L$ the bit length of the input data, $\overline{||M||_F}$ an upper bound to the Frobenius norm of the linear systems of equations that appear, $||M||_F$, and $\barκ$ an upper bound to the condition number $κ$ of those systems of equations. This represents a quantum speed-up in the number $n$ of variables in the cost function with respect to the comparable classical Interior Point algorithms when the initial matrix of the problem $A$ is dense: if we substitute the quantum part of the algorithm by classical algorithms such as Conjugate Gradient Descent, that would mean the whole algorithm has complexity $O(L\sqrt{n}(n+m)^2\barκ \log(ε^{-1}))$, or with exact methods, at least $O(L\sqrt{n}(n+m)^{2.373})$. Also, in contrast with any Quantum Linear System Algorithm, the algorithm described in this article outputs a classical description of the solution vector, and the value of the optimal solution. △ Less

Submitted 6 July, 2020; v1 submitted 18 February, 2019; originally announced February 2019.

Comments: Revtex file, color figures. Minor changes and typo correction from previous version

Journal ref: Journal of Physics A: Mathematical and Theretical 2020

Quench dynamics and zero-energy modes: the case of the Creutz model

Authors: R. Jafari, Henrik Johannesson, A. Langari, M. A. Martin-Delgado

Abstract: In most lattice models, the closing of a band gap typically occurs at high-symmetry points in the Brillouin zone. Differently, in the Creutz model $-$ describing a system of spinless fermions hopping on a two-leg ladder pierced by a magnetic field $-$ the gap closing at the quantum phase transition between the two topologically nontrivial phases of the model can be moved by tuning the hopping ampl… ▽ More In most lattice models, the closing of a band gap typically occurs at high-symmetry points in the Brillouin zone. Differently, in the Creutz model $-$ describing a system of spinless fermions hopping on a two-leg ladder pierced by a magnetic field $-$ the gap closing at the quantum phase transition between the two topologically nontrivial phases of the model can be moved by tuning the hopping amplitudes. We take advantage of this property to examine the nonequilibrium dynamics of the model after a sudden quench of the magnetic flux through the plaquettes of the ladder. For a quench to one of the equilibrium quantum critical points we find that the revival period of the Loschmidt echo $-$ measuring the overlap between initial and time-evolved states $-$ is controlled by the gap closing zero-energy modes. In particular, and contrary to expectations, the revival period of the Loschmidt echo for a finite ladder does not scale linearly with size but exhibits jumps determined by the presence or absence of zero-energy modes. We further investigate the conditions for the appearance of dynamical quantum phase transitions in the model and find that, for a quench {\em to} an equilibrium critical point, such transitions occur only for ladders of sizes which host zero-energy modes. Exploiting concepts from quantum thermodynamics, we show that the average work and the irreversible work per lattice site exhibit a weak dependence on the size of the system after a quench {\em across} an equilibrium critical point, suggesting that quenching into a different phase induces effective correlations among the particles. △ Less

Submitted 11 February, 2019; v1 submitted 6 November, 2018; originally announced November 2018.

Journal ref: Phys. Rev. B 99, 054302 (2019)

Quantum Error Correction with the Semion Code

Authors: Guillaume Dauphinais, Laura Ortiz, Santiago Varona, Miguel Angel Martin-Delgado

Abstract: We present a full quantum error correcting procedure with the semion code: an off-shell extension of the double semion model. We construct open strings operators that recover the quantum memory from arbitrary errors and closed string operators that implement the basic logical operations for information processing. Physically, the new open string operators provide a detailed microscopic description… ▽ More We present a full quantum error correcting procedure with the semion code: an off-shell extension of the double semion model. We construct open strings operators that recover the quantum memory from arbitrary errors and closed string operators that implement the basic logical operations for information processing. Physically, the new open string operators provide a detailed microscopic description of the creation of semions at their endpoints. Remarkably, topological properties of the string operators are determined using fundamental properties of the Hamiltonian, namely the fact that it is composed of commuting local terms squaring to the identity. In all, the semion code is a topological code that, unlike previously studied topological codes, it is of non-CSS type and fits into the stabilizer formalism. This is in sharp contrast with previous attempts yielding non-commutative codes. △ Less

Submitted 27 May, 2019; v1 submitted 18 October, 2018; originally announced October 2018.

Comments: REVTeX 4 file, color figures

Journal ref: New J. Phys. 21, 053035 (2019)

Twisted Fracton Models in Three Dimensions

Authors: Hao Song, Abhinav Prem, Sheng-Jie Huang, M. A. Martin-Delgado

Abstract: We study novel three-dimensional gapped quantum phases of matter which support quasiparticles with restricted mobility, including immobile "fracton" excitations. So far, most existing fracton models may be instructively viewed as generalized Abelian lattice gauge theories. Here, by analogy with Dijkgraaf-Witten topological gauge theories, we discover a natural generalization of fracton models, obt… ▽ More We study novel three-dimensional gapped quantum phases of matter which support quasiparticles with restricted mobility, including immobile "fracton" excitations. So far, most existing fracton models may be instructively viewed as generalized Abelian lattice gauge theories. Here, by analogy with Dijkgraaf-Witten topological gauge theories, we discover a natural generalization of fracton models, obtained by twisting the gauge symmetries. Introducing generalized gauge transformation operators carrying an extra phase factor depending on local configurations, we construct a plethora of exactly solvable three-dimensional models, which we dub "twisted fracton models." A key result of our approach is to demonstrate the existence of rich non-Abelian fracton phases of distinct varieties in a three-dimensional system with finite-range interactions. For an accurate characterization of these novel phases, the notion of being inextricably non-Abelian is introduced for fractons and quasiparticles with one-dimensional mobility, referring to their new behavior of displaying braiding statistics that is, and remains, non-Abelian regardless of which quasiparticles with higher mobility are added to or removed from them. We also analyze these models by embedding them on a three-torus and computing their ground state degeneracies, which exhibit a surprising and novel dependence on the system size in the non-Abelian fracton phases. Moreover, as an important advance in the study of fracton order, we develop a general mathematical framework which systematically captures the fusion and braiding properties of fractons and other quasiparticles with restricted mobility. △ Less

Submitted 2 August, 2018; v1 submitted 17 May, 2018; originally announced May 2018.

Comments: revtex4 file, color figures. (v2) The notion of being inextricably non-Abelian is introduced for fractons and quasiparticles with one-dimensional mobility in characterizing non-Abelian fracton phases. Introduction revised. References added

Journal ref: Phys. Rev. B 99, 155118 (2019)

Topological Phases in Nodeless Tetragonal Superconductors

Authors: Santiago Varona, Laura Ortiz, Oscar Viyuela, Miguel Angel Martin-Delgado

Abstract: We compute the topological phase diagram of 2D tetragonal superconductors for the only possible nodeless pairing channels compatible with that crystal symmetry. Subject to a Zeeman field and spin-orbit coupling, we demonstrate that these superconductors show surprising topological features: non-trivial high Chern numbers, massive edge states, and zero-energy modes out of high symmetry points, even… ▽ More We compute the topological phase diagram of 2D tetragonal superconductors for the only possible nodeless pairing channels compatible with that crystal symmetry. Subject to a Zeeman field and spin-orbit coupling, we demonstrate that these superconductors show surprising topological features: non-trivial high Chern numbers, massive edge states, and zero-energy modes out of high symmetry points, even though the edge states remain topologically protected. Interestingly, one of these pairing symmetries, $d+id$, has been proposed to describe materials such as water-intercalated sodium cobaltates, bilayer silicene or highly doped monolayer graphene, which opens the way for further applications of our results. △ Less

Submitted 23 January, 2019; v1 submitted 30 March, 2018; originally announced April 2018.

Comments: Minor corrections, results unchanged

Journal ref: Journal of Physics: Condensed Matter, 30, 395602 (2018)

The effects of thermal and correlated noise on magnons in a quantum ferromagnet

Authors: Jan Jeske, Ángel Rivas, Muhammad H. Ahmed, Miguel A. Martin-Delgado, Jared H. Cole

Abstract: The dynamics and thermal equilibrium of spin waves (magnons) in a quantum ferromagnet as well as the macroscopic magnetisation are investigated. Thermal noise due to an interaction with lattice phonons and the effects of spatial correlations in the noise are considered. We first present a Markovian master equation approach with analytical solutions for any homogeneous spatial correlation function… ▽ More The dynamics and thermal equilibrium of spin waves (magnons) in a quantum ferromagnet as well as the macroscopic magnetisation are investigated. Thermal noise due to an interaction with lattice phonons and the effects of spatial correlations in the noise are considered. We first present a Markovian master equation approach with analytical solutions for any homogeneous spatial correlation function of the noise. We find that spatially correlated noise increases the decay rate of magnons with low wave vectors to their thermal equilibrium, which also leads to a faster decay of the ferromagnet's magnetisation to its steady-state value. For long correlation lengths and higher temperature we find that additionally there is a component of the magnetisation which decays very slowly, due to a reduced decay rate of fast magnons. This effect could be useful for fast and noise-protected quantum or classical information transfer and magnonics. We further compare ferromagnetic and antiferromagnetic behaviour in noisy environments and find qualitatively similar behaviour in Ohmic but fundamentally different behaviour in super-Ohmic environments. △ Less

Submitted 5 November, 2018; v1 submitted 27 February, 2018; originally announced February 2018.

Journal ref: New J. Phys. 20 093017 (2018)

Twins Percolation for Qubit Losses in Topological Color Codes

Authors: D. Vodola, D. Amaro, M. A. Martin-Delgado, M. Müller

Abstract: We establish and explore a new connection between quantum information theory and classical statistical mechanics by studying the problem of qubit losses in 2D topological color codes. We introduce a protocol to cope with qubit losses, which is based on the identification and removal of a twin qubit from the code, and which guarantees the recovery of a valid three-colorable and trivalent reconstruc… ▽ More We establish and explore a new connection between quantum information theory and classical statistical mechanics by studying the problem of qubit losses in 2D topological color codes. We introduce a protocol to cope with qubit losses, which is based on the identification and removal of a twin qubit from the code, and which guarantees the recovery of a valid three-colorable and trivalent reconstructed color code lattice. Moreover, we show that determining the corresponding qubit loss error threshold is equivalent to a new generalized classical percolation process. We numerically compute the associated qubit loss thresholds for two families of 2D color code and find that these are close to satisfying the fundamental limit $p_\text{fund} = 0.461 \pm 0.005$ close to the 50% as imposed by the no-cloning theorem. Our findings reveal a new connection between topological color codes and percolation theory, show high robustness of color codes against qubit loss, and are relevant for implementations of quantum error correction in various physical platforms. △ Less

Submitted 7 March, 2018; v1 submitted 13 February, 2018; originally announced February 2018.

Comments: 6+7 pages, 3+7 figures

Journal ref: Phys. Rev. Lett. 121, 060501 (2018)

Staircase to Higher-Order Topological Phase Transitions

Authors: P. Cats, A. Quelle, O. Viyuela, M. A. Martin-Delgado, C. Morais Smith

Abstract: We find a series of topological phase transitions of increasing order, beyond the more standard second-order phase transition in a one-dimensional topological superconductor. The jumps in the order of the transitions depend on the range of the pairing interaction, which is parametrized by an algebraic decay with exponent $α$. Remarkably, in the limit $α= 1$ the order of the topological transition… ▽ More We find a series of topological phase transitions of increasing order, beyond the more standard second-order phase transition in a one-dimensional topological superconductor. The jumps in the order of the transitions depend on the range of the pairing interaction, which is parametrized by an algebraic decay with exponent $α$. Remarkably, in the limit $α= 1$ the order of the topological transition becomes infinite. We compute the critical exponents for the series of higher-order transitions in exact form and find that they fulfill the hyperscaling relation. We also study the critical behaviour at the boundary of the system and discuss potential experimental platforms of magnetic atoms in superconductors. △ Less

Submitted 29 March, 2018; v1 submitted 16 October, 2017; originally announced October 2017.

Comments: 5+5pages, 7 figures. Accepted as a Rapid Communication

Journal ref: Phys. Rev. B 97, 121106 (2018)

Localization and oscillations of Majorana fermions in a two-dimensional electron gas coupled with $d$-wave superconductors

Authors: Laura Ortiz, Santiago Varona, Oscar Viyuela, Miguel Angel Martin-Delgado

Abstract: We study the localization and oscillation properties of the Majorana fermions that arise in a two-dimensional electron gas (2DEG) with spin-orbit coupling (SOC) and a Zeeman field coupled with a d-wave superconductor. Despite the angular dependence of the d-wave pairing, localization and oscillation properties are found to be similar to the ones seen in conventional s-wave superconductors. In addi… ▽ More We study the localization and oscillation properties of the Majorana fermions that arise in a two-dimensional electron gas (2DEG) with spin-orbit coupling (SOC) and a Zeeman field coupled with a d-wave superconductor. Despite the angular dependence of the d-wave pairing, localization and oscillation properties are found to be similar to the ones seen in conventional s-wave superconductors. In addition, we study a microscopic lattice version of the previous system that can be characterized by a topological invariant. We derive its real space representation that involves nearest and next-to-nearest-neighbors pairing. Finally, we show that the emerging chiral Majorana fermions are indeed robust against static disorder. This analysis has potential applications to quantum simulations and experiments in high-$T_c$ superconductors. △ Less

Submitted 15 February, 2018; v1 submitted 19 September, 2017; originally announced September 2017.

Comments: revtex4 file, color figures

Journal ref: Phys. Rev. B 97, 064501 (2018)

Chiral Topological Superconductors Enhanced by Long-Range Interactions

Authors: Oscar Viyuela, Liang Fu, Miguel Angel Martin-Delgado

Abstract: We study the phase diagram and edge states of a two-dimensional p-wave superconductor with long-range hopping and pairing amplitudes. New topological phases and quasiparticles different from the usual short-range model are obtained. When both hopping and pairing terms decay with the same exponent, one of the topological chiral phases with propagating Majorana edge states gets significantly enhance… ▽ More We study the phase diagram and edge states of a two-dimensional p-wave superconductor with long-range hopping and pairing amplitudes. New topological phases and quasiparticles different from the usual short-range model are obtained. When both hopping and pairing terms decay with the same exponent, one of the topological chiral phases with propagating Majorana edge states gets significantly enhanced by long-range couplings. On the other hand, when the long-range pairing amplitude decays more slowly than the hopping, we discover new topological phases where propagating Majorana fermions at each edge pair nonlocally and become gapped even in the thermodynamic limit. Remarkably, these nonlocal edge states are still robust, remain separated from the bulk, and are localized at both edges at the same time. The inclusion of long-range effects is potentially applicable to recent experiments with magnetic impurities and islands in 2D superconductors. △ Less

Submitted 11 January, 2018; v1 submitted 7 July, 2017; originally announced July 2017.

Comments: 5 pages, 3 figures. Close to published version

Journal ref: Phys. Rev. Lett. 120, 017001 (2018)

Observation of Majorization Principle for quantum algorithms via 3-D integrated photonic circuits

Authors: Fulvio Flamini, Niko Viggianiello, Taira Giordani, Marco Bentivegna, Nicoló Spagnolo, Andrea Crespi, Giacomo Corrielli, Roberto Osellame, Miguel Angel Martin-Delgado, Fabio Sciarrino

Abstract: The Majorization Principle is a fundamental statement governing the dynamics of information processing in optimal and efficient quantum algorithms. While quantum computation can be modeled to be reversible, due to the unitary evolution undergone by the system, these quantum algorithms are conjectured to obey a quantum arrow of time dictated by the Majorization Principle: the probability distributi… ▽ More The Majorization Principle is a fundamental statement governing the dynamics of information processing in optimal and efficient quantum algorithms. While quantum computation can be modeled to be reversible, due to the unitary evolution undergone by the system, these quantum algorithms are conjectured to obey a quantum arrow of time dictated by the Majorization Principle: the probability distribution associated to the outcomes gets ordered step-by-step until achieving the result of the computation. Here we report on the experimental observation of the effects of the Majorization Principle for two quantum algorithms, namely the quantum fast Fourier transform and a recently introduced validation protocol for the certification of genuine many-boson interference. The demonstration has been performed by employing integrated 3-D photonic circuits fabricated via femtosecond laser writing technique, which allows to monitor unambiguously the effects of majorization along the execution of the algorithms. The measured observables provide a strong indication that the Majorization Principle holds true for this wide class of quantum algorithms, thus paving the way for a general tool to design new optimal algorithms with a quantum speedup. △ Less

Submitted 3 August, 2016; originally announced August 2016.

Comments: 6 pages, 4 figures + 2 pages Supplementary Information

Observation of topological Uhlmann phases with superconducting qubits

Authors: O. Viyuela, A. Rivas, S. Gasparinetti, A. Wallraff, S. Filipp, M. A. Martin-Delgado

Abstract: Topological insulators and superconductors at finite temperature can be characterized by the topological Uhlmann phase. However, a direct experimental measurement of this invariant has remained elusive in condensed matter systems. Here, we report a measurement of the topological Uhlmann phase for a topological insulator simulated by a system of entangled qubits in the IBM Quantum Experience platfo… ▽ More Topological insulators and superconductors at finite temperature can be characterized by the topological Uhlmann phase. However, a direct experimental measurement of this invariant has remained elusive in condensed matter systems. Here, we report a measurement of the topological Uhlmann phase for a topological insulator simulated by a system of entangled qubits in the IBM Quantum Experience platform. By making use of ancilla states, otherwise unobservable phases carrying topological information about the system become accessible, enabling the experimental determination of a complete phase diagram including environmental effects. We employ a state-independent measurement protocol which does not involve prior knowledge of the system state. The proposed measurement scheme is extensible to interacting particles and topological models with a large number of bands. △ Less

Submitted 12 February, 2018; v1 submitted 29 July, 2016; originally announced July 2016.

Comments: RevTex4 file, color figures

Journal ref: npj Quantum Information 4, 10 (2018)

Topological Heat Transport and Symmetry-Protected Boson Currents

Authors: Ángel Rivas, Miguel A. Martin-Delgado

Abstract: The study of non-equilibrium properties in topological systems is of practical and fundamental importance. Here, we analyze the stationary properties of a two-dimensional bosonic Hofstadter lattice coupled to two thermal baths in the quantum open-system formalism. Novel phenomena appear like chiral edge heat currents that are the out-of-equilibrium counterparts of the zero-temperature edge current… ▽ More The study of non-equilibrium properties in topological systems is of practical and fundamental importance. Here, we analyze the stationary properties of a two-dimensional bosonic Hofstadter lattice coupled to two thermal baths in the quantum open-system formalism. Novel phenomena appear like chiral edge heat currents that are the out-of-equilibrium counterparts of the zero-temperature edge currents. They support a new concept of dissipative symmetry-protection, where a set of discrete symmetries protects topological heat currents, differing from the symmetry-protection devised in closed systems and zero-temperature. Remarkably, one of these currents flows opposite to the decreasing external temperature gradient. As the starting point, we consider the case of a single external reservoir already showing prominent results like thermal erasure effects and topological thermal currents. Our results are experimentally accessible with platforms like photonics systems and optical lattices. △ Less

Submitted 25 July, 2017; v1 submitted 24 June, 2016; originally announced June 2016.

Comments: RevTeX4 file, color figures

Journal ref: Sci. Rep. 7, 6350 (2017)

A bilayer Double Semion Model with Symmetry-Enriched Topological Order

Authors: L. Ortiz, M. A. Martin-Delgado

Abstract: We construct a new model of two-dimensional quantum spin systems that combines intrinsic topo- logical orders and a global symmetry called flavour symmetry. It is referred as the bilayer Doubled Semion model (bDS) and is an instance of symmetry-enriched topological order. A honeycomb bi- layer lattice is introduced to combine a Double Semion Topolgical Order with a global spin-flavour symmetry to… ▽ More We construct a new model of two-dimensional quantum spin systems that combines intrinsic topo- logical orders and a global symmetry called flavour symmetry. It is referred as the bilayer Doubled Semion model (bDS) and is an instance of symmetry-enriched topological order. A honeycomb bi- layer lattice is introduced to combine a Double Semion Topolgical Order with a global spin-flavour symmetry to get the fractionalization of its quasiparticles. The bDS model exhibits non-trival braid- ing self-statistics of excitations and its dual model constitutes a Symmetry-Protected Topological Order with novel edge states. This dual model gives rise to a bilayer Non-Trivial Paramagnet that is invariant under the flavour symmetry and the well-known spin flip symmetry. △ Less

Submitted 14 September, 2016; v1 submitted 28 April, 2016; originally announced April 2016.

Comments: revtex4 file, color figures

Journal ref: Annals of Physics 375 (2016) 193-226

Error tolerance of topological codes with independent bit-flip and measurement errors

Authors: Ruben S. Andrist, Helmut G. Katzgraber, H. Bombin, M. A. Martin-Delgado

Abstract: Topological quantum error correction codes are currently among the most promising candidates for efficiently dealing with the decoherence effects inherently present in quantum devices. Numerically, their theoretical error threshold can be calculated by mapping the underlying quantum problem to a related classical statistical-mechanical spin system with quenched disorder. Here, we present results f… ▽ More Topological quantum error correction codes are currently among the most promising candidates for efficiently dealing with the decoherence effects inherently present in quantum devices. Numerically, their theoretical error threshold can be calculated by mapping the underlying quantum problem to a related classical statistical-mechanical spin system with quenched disorder. Here, we present results for the general fault-tolerant regime, where we consider both qubit and measurement errors. However, unlike in previous studies, here we vary the strength of the different error sources independently. Our results highlight peculiar differences between toric and color codes. This study complements previous results published in New J. Phys. 13, 083006 (2011). △ Less

Submitted 12 July, 2016; v1 submitted 29 March, 2016; originally announced March 2016.

Comments: 6 pages, 3 figures

Journal ref: Phys. Rev. A 94, 012318 (2016)

Quantum Enhanced Energy Distribution for Information Heat Engines

Authors: J. M. Diaz de la Cruz, M. A. Martin-Delgado

Abstract: A new scenario for energy distribution, security and shareability is presented that assumes the availability of quantum information heat engines and a thermal bath. It is based on the convertibility between entropy and work in the presence of a thermal reservoir. Our approach to the informational content of physical systems that are distributed between users is complementary to the conventional pe… ▽ More A new scenario for energy distribution, security and shareability is presented that assumes the availability of quantum information heat engines and a thermal bath. It is based on the convertibility between entropy and work in the presence of a thermal reservoir. Our approach to the informational content of physical systems that are distributed between users is complementary to the conventional perspective of quantum communication. The latter places the value on the unpredictable content of the transmitted quantum states, while our interest focuses on their certainty. Some well-known results in quantum communication are reused in this context. Particularly, we describe a way to securely distribute quantum states to be used for unlocking energy from thermal sources. We also consider some multi-partite entangled and classically correlated states for a collaborative multi-user sharing of work extraction possibilities. Besides, the relation between the communication and work extraction capabilities is analyzed and written as an equation. △ Less

Submitted 1 August, 2016; v1 submitted 25 March, 2016; originally announced March 2016.

Comments: 10 pages, 8 figures, color. Section V has been added

Iterative Phase Optimisation of Elementary Quantum Error Correcting Codes

Authors: M. Müller, A. Rivas, E. A. Martínez, D. Nigg, P. Schindler, T. Monz, R. Blatt, M. A. Martin-Delgado

Abstract: Performing experiments on small-scale quantum computers is certainly a challenging endeavor. Many parameters need to be optimized to achieve high-fidelity operations. This can be done efficiently for operations acting on single qubits as errors can be fully characterized. For multi-qubit operations, though, this is no longer the case as in the most general case analyzing the effect of the operatio… ▽ More Performing experiments on small-scale quantum computers is certainly a challenging endeavor. Many parameters need to be optimized to achieve high-fidelity operations. This can be done efficiently for operations acting on single qubits as errors can be fully characterized. For multi-qubit operations, though, this is no longer the case as in the most general case analyzing the effect of the operation on the system requires a full state tomography for which resources scale exponentially with the system size. Furthermore, in recent experiments additional electronic levels beyond the two-level system encoding the qubit have been used to enhance the capabilities of quantum information processors, which additionally increases the number of parameters that need to be controlled. For the optimization of the experimental system for a given task (e.g.~a quantum algorithm), one has to find a satisfactory error model and also efficient observables to estimate the parameters of the model. In this manuscript we demonstrate a method to optimize the encoding procedure for a small quantum error correction code in the presence of unknown but constant phase shifts. The method, which we implement here on a small-scale linear ion-trap quantum computer, is readily applicable to other AMO platforms for quantum information processing. △ Less

Submitted 30 August, 2016; v1 submitted 1 March, 2016; originally announced March 2016.

Comments: 13 pages including appendices, 5 color figures, final version

Journal ref: Phys. Rev. X 6, 031030 (2016)

Topological Massive Dirac Edge Modes and Long-Range Superconducting Hamiltonians

Authors: O. Viyuela, D. Vodola, G. Pupillo, M. A. Martin-Delgado

Abstract: We discover novel topological effects in the one-dimensional Kitaev chain modified by long-range Hamiltonian deformations in the hopping and pairing terms. This class of models display symmetry-protected topological order measured by the Berry/Zak phase of the lower band eigenvector and the winding number of the Hamiltonians. For exponentially-decaying hopping amplitudes, the topological sector ca… ▽ More We discover novel topological effects in the one-dimensional Kitaev chain modified by long-range Hamiltonian deformations in the hopping and pairing terms. This class of models display symmetry-protected topological order measured by the Berry/Zak phase of the lower band eigenvector and the winding number of the Hamiltonians. For exponentially-decaying hopping amplitudes, the topological sector can be significantly augmented as the penetration length increases, something experimentally achievable. For power-law decaying superconducting pairings, the massless Majorana modes at the edges get paired together into a massive non-local Dirac fermion localised at both edges of the chain: a new topological quasiparticle that we call topological massive Dirac fermion. This topological phase has fractional topological numbers as a consequence of the long-range couplings. Possible applications to current experimental setups and topological quantum computation are also discussed. △ Less

Submitted 19 September, 2016; v1 submitted 16 November, 2015; originally announced November 2015.

Comments: RevTex4 file, color figures, close to published version

Journal ref: Phys. Rev. B 94, 125121 (2016)

Quantum Simulation of a Topological Mott Insulator with Rydberg Atoms in a Lieb Lattice

Authors: A. Dauphin, M. Müller, M. A. Martin-Delgado

Abstract: We propose a realistic scheme to quantum simulate the so-far experimentally unobserved topological Mott insulator phase -- an interaction-driven topological insulator -- using cold atoms in an optical Lieb lattice. To this end, we study a system of spinless fermions in a Lieb lattice, exhibiting repulsive nearest and next-to-nearest neighbor interactions, and derive the associated zero temperature… ▽ More We propose a realistic scheme to quantum simulate the so-far experimentally unobserved topological Mott insulator phase -- an interaction-driven topological insulator -- using cold atoms in an optical Lieb lattice. To this end, we study a system of spinless fermions in a Lieb lattice, exhibiting repulsive nearest and next-to-nearest neighbor interactions, and derive the associated zero temperature phase diagram within mean-field approximation. In particular, we analyze how the interactions can dynamically generate a charge density wave ordered, a nematic as well as a topologically non-trivial quantum anomalous Hall phase. We characterize the topology of the different phases by the Chern number and discuss the possibility of phase coexistence. Based on the identified phases, we propose a realistic implementation of this model using cold Rydberg-dressed atoms in an optical lattice. The scheme, which allows one to access in particular the topological Mott insulator phase, robustly and independently of its exact position in parameter space, merely requires global, always-on off-resonant laser coupling to Rydberg states and is feasible with state-of-the-art experimental techniques. △ Less

Submitted 18 October, 2015; v1 submitted 1 October, 2015; originally announced October 2015.

Comments: 8 pages, 6 color figures (incl. appendices)

Journal ref: Phys. Rev. A 93, 043611 (2016)

Symmetry-protected Topological Phases at Finite Temperature

Authors: O. Viyuela, A. Rivas, M. A. Martin-Delgado

Abstract: We have applied the recently developed theory of topological Uhlmann numbers to a representative model of a topological insulator in two dimensions, the Qi-Wu-Zhang model. We have found a stable symmetry-protected topological (SPT) phase under external thermal fluctuations in two-dimensions. A complete phase diagram for this model is computed as a function of temperature and coupling constants in… ▽ More We have applied the recently developed theory of topological Uhlmann numbers to a representative model of a topological insulator in two dimensions, the Qi-Wu-Zhang model. We have found a stable symmetry-protected topological (SPT) phase under external thermal fluctuations in two-dimensions. A complete phase diagram for this model is computed as a function of temperature and coupling constants in the original Hamiltonian. It shows the appearance of large stable phases of matter with topological properties compatible with thermal fluctuations or external noise and the existence of critical lines separating abruptly trivial phases from topological phases. These novel critical temperatures represent thermal topological phase transitions. The initial part of the paper comprises a self-contained explanation of the Uhlmann geometric phase needed to understand the topological properties that it may acquire when applied to topological insulators and superconductors. △ Less

Submitted 30 June, 2015; v1 submitted 4 February, 2015; originally announced February 2015.

Comments: Contribution to the focus issue on "Artificial Graphene". Edited by Maciej Lewenstein, Vittorio Pellegrini, Marco Polini and Mordechai (Moti) Segev

Journal ref: 2D Mater. 2 034006 (2015)

Quantum Google Algorithm: Construction and Application to Complex Networks

Authors: G. D. Paparo, M. Müller, F. Comellas, M. A. Martin-Delgado

Abstract: We review the main findings on the ranking capabilities of the recently proposed Quantum PageRank algorithm (G.D. Paparo et al., Sci. Rep. 2, 444 (2012) and G.D. Paparo et al., Sci. Rep. 3, 2773 (2013)) applied to large complex networks. The algorithm has been shown to identify unambiguously the underlying topology of the network and to be capable of clearly highlighting the structure of secondary… ▽ More We review the main findings on the ranking capabilities of the recently proposed Quantum PageRank algorithm (G.D. Paparo et al., Sci. Rep. 2, 444 (2012) and G.D. Paparo et al., Sci. Rep. 3, 2773 (2013)) applied to large complex networks. The algorithm has been shown to identify unambiguously the underlying topology of the network and to be capable of clearly highlighting the structure of secondary hubs of networks. Furthermore, it can resolve the degeneracy in importance of the low-lying part of the list of rankings. Examples of applications include real-world instances from the WWW, which typically display a scale-free network structure and models of hierarchical networks. The quantum algorithm has been shown to display an increased stability with respect to a variation of the damping parameter, present in the Google algorithm, and a more clearly pronounced power-law behaviour in the distribution of importance among the nodes, as compared to the classical algorithm. △ Less

Submitted 12 September, 2014; originally announced September 2014.

Journal ref: Eur. Phys. J. Plus (2014) 129: 150

Two-Dimensional Density-Matrix Topological Fermionic Phases: Topological Uhlmann Numbers

Authors: O. Viyuela, A. Rivas, M. A. Martin-Delgado

Abstract: We construct a topological invariant that classifies density matrices of symmetry-protected topological orders in two-dimensional fermionic systems. As it is constructed out of the previously introduced Uhlmann phase, we refer to it as the topological Uhlmann number ${\rm n}_{\rm U}$. With it, we study thermal topological phases in several two-dimensional models of topological insulators and super… ▽ More We construct a topological invariant that classifies density matrices of symmetry-protected topological orders in two-dimensional fermionic systems. As it is constructed out of the previously introduced Uhlmann phase, we refer to it as the topological Uhlmann number ${\rm n}_{\rm U}$. With it, we study thermal topological phases in several two-dimensional models of topological insulators and superconductors, computing phase diagrams where the temperature $T$ is on an equal footing with the coupling constants in the Hamiltonian. Moreover, we find novel thermal-topological transitions between two non-trivial phases in a model with high Chern numbers. At small temperature we recover the standard topological phases as the Uhlmann number approaches to the Chern number. △ Less

Submitted 25 August, 2014; v1 submitted 23 May, 2014; originally announced May 2014.

Comments: RevTex4 file, color figures. Close to published version

Journal ref: Phys. Rev. Lett. 113, 076408 (2014)

Experimental Quantum Computations on a Topologically Encoded Qubit

Authors: Daniel Nigg, Markus Mueller, Esteban A. Martinez, Philipp Schindler, Markus Hennrich, Thomas Monz, Miguel A. Martin-Delgado, Rainer Blatt

Abstract: The construction of a quantum computer remains a fundamental scientific and technological challenge, in particular due to unavoidable noise. Quantum states and operations can be protected from errors using protocols for fault-tolerant quantum computing (FTQC). Here we present a step towards this by implementing a quantum error correcting code, encoding one qubit in entangled states distributed ove… ▽ More The construction of a quantum computer remains a fundamental scientific and technological challenge, in particular due to unavoidable noise. Quantum states and operations can be protected from errors using protocols for fault-tolerant quantum computing (FTQC). Here we present a step towards this by implementing a quantum error correcting code, encoding one qubit in entangled states distributed over 7 trapped-ion qubits. We demonstrate the capability of the code to detect one bit flip, phase flip or a combined error of both, regardless on which of the qubits they occur. Furthermore, we apply combinations of the entire set of logical single-qubit Clifford gates on the encoded qubit to explore its computational capabilities. The implemented 7-qubit code is the first realization of a complete Calderbank-Shor-Steane (CSS) code and constitutes a central building block for FTQC schemes based on concatenated elementary quantum codes. It also represents the smallest fully functional instance of the color code, opening a route towards topological FTQC. △ Less

Submitted 21 March, 2014; originally announced March 2014.

Comments: 7 pages, 4 figures and appendix

Quantum Information Engines with Many-Body States attaining optimal Extractable Work with Quantum Control

Authors: J. M. Diazdelacruz, M. A. Martin-Delgado

Abstract: We introduce quantum information engines that extract work from quantum states and a single thermal reservoir. They may operate under three general conditions: i/ Unitarily Steered evolution (US); ii/ Irreversible Thermalization (IT) and iii/ Isothermal Relaxation (IR), and hence are called USITIR machines. They include novel engines without traditional feedback control mechanisms, as well as vers… ▽ More We introduce quantum information engines that extract work from quantum states and a single thermal reservoir. They may operate under three general conditions: i/ Unitarily Steered evolution (US); ii/ Irreversible Thermalization (IT) and iii/ Isothermal Relaxation (IR), and hence are called USITIR machines. They include novel engines without traditional feedback control mechanisms, as well as versions which also include them. Explicit constructions of USITIR engines are presented for one- and two-qubit states and their maximum extractable work is computed, which is optimal. Optimality is achieved when the notions of controllable thermalizability and density matrix controllability are fullfilled. Then, many-body extensions of USITIR engines are also analyzed and conditions for optimal work extraction are identified. When they are not met, we measure their lack of optimality by means of newly defined uncontrollable entropies, that are explicitly computed for some selected examples. This includes cases of distinguishable and indistinguishable particles. △ Less

Submitted 20 March, 2014; v1 submitted 27 January, 2014; originally announced January 2014.

Comments: RevTeX4 file, color figures, close to published version

Journal ref: Physical Review A, Vol. 89, No. 3 (2014)

Quantum speedup for active learning agents

Authors: Giuseppe Davide Paparo, Vedran Dunjko, Adi Makmal, Miguel Angel Martin-Delgado, Hans J. Briegel

Abstract: Can quantum mechanics help us in building intelligent robots and agents? One of the defining characteristics of intelligent behavior is the capacity to learn from experience. However, a major bottleneck for agents to learn in any real-life situation is the size and complexity of the corresponding task environment. Owing to, e.g., a large space of possible strategies, learning is typically slow. Ev… ▽ More Can quantum mechanics help us in building intelligent robots and agents? One of the defining characteristics of intelligent behavior is the capacity to learn from experience. However, a major bottleneck for agents to learn in any real-life situation is the size and complexity of the corresponding task environment. Owing to, e.g., a large space of possible strategies, learning is typically slow. Even for a moderate task environment, it may simply take too long to rationally respond to a given situation. If the environment is impatient, allowing only a certain time for a response, an agent may then be unable to cope with the situation and to learn at all. Here we show that quantum physics can help and provide a significant speed-up for active learning as a genuine problem of artificial intelligence. We introduce a large class of quantum learning agents for which we show a quadratic boost in their active learning efficiency over their classical analogues. This result will be particularly relevant for applications involving complex task environments. △ Less

Submitted 14 July, 2014; v1 submitted 20 January, 2014; originally announced January 2014.

Comments: Minor updates, 14 pages, 3 figures

Journal ref: Phys. Rev. X 4, 031002 (2014)

Efficient algorithm to compute the Berry conductivity

Authors: A. Dauphin, M. Müller, M. A. Martin-Delgado

Abstract: We propose and construct a numerical algorithm to calculate the Berry conductivity in topological band insulators. The method is applicable to cold atom systems as well as solid state setups, both for the insulating case where the Fermi energy lies in the gap between two bulk bands as well as in the metallic regime and interpolates smoothly between both regimes. The algorithm is gauge-invariant by… ▽ More We propose and construct a numerical algorithm to calculate the Berry conductivity in topological band insulators. The method is applicable to cold atom systems as well as solid state setups, both for the insulating case where the Fermi energy lies in the gap between two bulk bands as well as in the metallic regime and interpolates smoothly between both regimes. The algorithm is gauge-invariant by construction, efficient and yields the Berry conductivity with known and controllable statistical error bars. We apply the algorithm to several paradigmatic models in the field of topological insulators, including Haldane's model on the honeycomb lattice, the multi-band Hofstadter model and the BHZ model, which describes the 2D spin Hall effect observed in CdTe/HgTe/CdTe quantum well heterostructures. △ Less

Submitted 21 July, 2014; v1 submitted 8 January, 2014; originally announced January 2014.

Comments: Revtex4 file, color figures. Revised version published in New Journal of Physics

Journal ref: New J. Phys. 16, 073016 (2014)