-
Optimizing Unitary Coupled Cluster Wave Functions on Quantum Hardware: Error Bound and Resource-Efficient Optimizer
Authors:
Martin Plazanet,
Thomas Ayral
Abstract:
In this work, we study the projective quantum eigensolver (PQE) approach to optimizing unitary coupled cluster wave functions on quantum hardware, as introduced in arXiv:2102.00345. The projective quantum eigensolver is a hybrid quantum-classical algorithm which, by optimizing a unitary coupled cluster wave function, aims at computing the ground state of many-body systems. Instead of trying to min…
▽ More
In this work, we study the projective quantum eigensolver (PQE) approach to optimizing unitary coupled cluster wave functions on quantum hardware, as introduced in arXiv:2102.00345. The projective quantum eigensolver is a hybrid quantum-classical algorithm which, by optimizing a unitary coupled cluster wave function, aims at computing the ground state of many-body systems. Instead of trying to minimize the energy of the system like in the variational quantum eigensolver, PQE uses projections of the Schrödinger equation to efficiently bring the trial state closer to an eigenstate of the Hamiltonian. In this work, we provide a mathematical study of the algorithm, which allows us to obtain a number of interesting results. We first show that one can derive a bound relating off-diagonal coefficients (residues) of the Hamiltonian to the error of the algorithm. This bound not only gives a formal motivation to the projective approach to optimizing unitary coupled cluster wavefunctions, but it also allows us to formulate a well-informed convergence criterion for residue-based optimizers. We then introduce a mathematical study of the classical optimization itself, and show that, using our results, one can derive a residue-based optimizer for which we present numerical evidence of superiority over both the optimization introduced in arXiv:2102.00345 and VQE optimized using the Broyden Fletcher Goldfarb Shanno (BFGS) method.
△ Less
Submitted 19 October, 2024;
originally announced October 2024.
-
Compact fermionic quantum state preparation with a natural-orbitalizing variational quantum eigensolving scheme
Authors:
Pauline Besserve,
Michel Ferrero,
Thomas Ayral
Abstract:
Assemblies of strongly interacting fermions, whether in a condensed-matter or a quantum chemistry context, range amongst the most promising candidate systems for which quantum computing platforms could provide an advantage. Near-term quantum state preparation is typically realized by means of the variational quantum eigensolver (VQE) algorithm. One of the main challenges to a successful implementa…
▽ More
Assemblies of strongly interacting fermions, whether in a condensed-matter or a quantum chemistry context, range amongst the most promising candidate systems for which quantum computing platforms could provide an advantage. Near-term quantum state preparation is typically realized by means of the variational quantum eigensolver (VQE) algorithm. One of the main challenges to a successful implementation of VQE lies in the sensitivity to noise exhibited by deep variational circuits. On the other hand, sufficient depth must be allowed to be able to reach a good approximation to the target state. In this work, we present a refined VQE scheme that consists in topping VQE with state-informed updates of the elementary fermionic modes (spin-orbitals). These updates consist in moving to the natural-orbital basis of the current, converged variational state, a basis we argue eases the task of state preparation. We test the method on the Hubbard model in the presence of experimentally relevant noise levels. For a fixed circuit structure, the method is shown to enhance the capabilities of the circuit to reach a state close to the target state without incurring too much overhead from shot noise. Moreover, coupled with an adaptive VQE scheme that constructs the circuit on the fly, we evidence reduced requirements on the depth of the circuit as the orbitals get updated.
△ Less
Submitted 20 June, 2024;
originally announced June 2024.
-
Enabling large-depth simulation of noisy quantum circuits with positive tensor networks
Authors:
Ambroise Müller,
Thomas Ayral,
Corentin Bertrand
Abstract:
Matrix product density operators (MPDOs) are tensor network representations of locally purified density matrices where each physical degree of freedom is associated to an environment degree of freedom. MPDOs have interesting properties for mixed state representations: guaranteed positivity by construction, efficient conservation of the trace and computation of local observables. However, they have…
▽ More
Matrix product density operators (MPDOs) are tensor network representations of locally purified density matrices where each physical degree of freedom is associated to an environment degree of freedom. MPDOs have interesting properties for mixed state representations: guaranteed positivity by construction, efficient conservation of the trace and computation of local observables. However, they have been challenging to use for noisy quantum circuit simulation, as the application of noise increases the dimension of the environment Hilbert space, leading to an exponential growth of bond dimensions. MPDOs also lack a unique canonical form, due to the freedom in the choice of basis for the environment Hilbert space, which leads to a vast variation of bond dimensions.
In this work, we present a systematic way to reduce the bond dimensions of MPDOs by disentangling the purified state. We optimize the basis for the environment Hilbert space by performing density matrix renormalization group (DMRG)-like sweeps of local 2-qubit basis optimization. Interestingly, we find that targeting only the disentanglement of the purified state leads to a reduction of the environment dimension. In other words, a compact MPDO representation requires a low-entanglement purified state.
We apply our compression method to the emulation of noisy random quantum circuits. Our technique allows us to keep bounded bond dimensions, and thus bounded memory, contrary to previous works on MPDOs, while keeping reasonable truncation fidelities.
△ Less
Submitted 29 February, 2024;
originally announced March 2024.
-
Hubbard physics with Rydberg atoms: using a quantum spin simulator to simulate strong fermionic correlations
Authors:
Antoine Michel,
Loïc Henriet,
Christophe Domain,
Antoine Browaeys,
Thomas Ayral
Abstract:
We propose a hybrid quantum-classical method to investigate the equilibrium physics and the dynamics of strongly correlated fermionic models with spin-based quantum processors. Our proposal avoids the usual pitfalls of fermion-to-spin mappings thanks to a slave-spin method which allows to approximate the original Hamiltonian into a sum of self-correlated free-fermions and spin Hamiltonians. Taking…
▽ More
We propose a hybrid quantum-classical method to investigate the equilibrium physics and the dynamics of strongly correlated fermionic models with spin-based quantum processors. Our proposal avoids the usual pitfalls of fermion-to-spin mappings thanks to a slave-spin method which allows to approximate the original Hamiltonian into a sum of self-correlated free-fermions and spin Hamiltonians. Taking as an example a Rydberg-based analog quantum processor to solve the interacting spin model, we avoid the challenges of variational algorithms or Trotterization methods. We explore the robustness of the method to experimental imperfections by applying it to the half-filled, single-orbital Hubbard model on the square lattice in and out of equilibrium. We show, through realistic numerical simulations of current Rydberg processors, that the method yields quantitatively viable results even in the presence of imperfections: it allows to gain insights into equilibrium Mott physics as well as the dynamics under interaction quenches. This method thus paves the way to the investigation of physical regimes -- whether out-of-equilibrium, doped, or multiorbital -- that are difficult to explore with classical processors.
△ Less
Submitted 19 June, 2024; v1 submitted 13 December, 2023;
originally announced December 2023.
-
On the feasibility of performing quantum chemistry calculations on quantum computers
Authors:
Thibaud Louvet,
Thomas Ayral,
Xavier Waintal
Abstract:
Quantum chemistry is envisioned as an early and disruptive application for quantum computers. Yet, closer scrutiny of the proposed algorithms shows that there are considerable difficulties along the way. Here, we propose two criteria for evaluating two leading quantum approaches for finding the ground state of molecules. The first criterion applies to the variational quantum eigensolver (VQE) algo…
▽ More
Quantum chemistry is envisioned as an early and disruptive application for quantum computers. Yet, closer scrutiny of the proposed algorithms shows that there are considerable difficulties along the way. Here, we propose two criteria for evaluating two leading quantum approaches for finding the ground state of molecules. The first criterion applies to the variational quantum eigensolver (VQE) algorithm. It sets an upper bound to the level of imprecision/decoherence that can be tolerated in quantum hardware as a function of the targeted precision, the number of gates and the typical energy contribution from states populated by decoherence processes. We find that decoherence is highly detrimental to the accuracy of VQE and performing relevant chemistry calculations would require performances that are expected for fault-tolerant quantum computers, not mere noisy hardware, even with advanced error mitigation techniques. Physically, the sensitivity of VQE to decoherence originates from the fact that, in VQE, the spectrum of the studied molecule has no correlation with the spectrum of the quantum hardware used to perform the computation.
The second criterion applies to the quantum phase estimation (QPE) algorithm, which is often presented as the go-to replacement of VQE upon availability of (noiseless) fault-tolerant quantum computers. QPE requires an input state with a large enough overlap with the sought-after ground state. We provide a criterion to estimate quantitatively this overlap based on the energy and the energy variance of said input state. Using input states from a variety of state-of-the-art classical methods, we show that the scaling of this overlap with system size does display the standard orthogonality catastrophe, namely an exponential suppression with system size. This in turns leads to an exponentially reduced QPE success probability.
△ Less
Submitted 3 October, 2024; v1 submitted 5 June, 2023;
originally announced June 2023.
-
Combining Matrix Product States and Noisy Quantum Computers for Quantum Simulation
Authors:
Baptiste Anselme Martin,
Thomas Ayral,
François Jamet,
Marko J. Rančić,
Pascal Simon
Abstract:
Matrix Product States (MPS) and Operators (MPO) have been proven to be a powerful tool to study quantum many-body systems but are restricted to moderately entangled states as the number of parameters scales exponentially with the entanglement entropy. While MPS can efficiently find ground states of 1D systems, their capacities are limited when simulating their dynamics, where the entanglement can…
▽ More
Matrix Product States (MPS) and Operators (MPO) have been proven to be a powerful tool to study quantum many-body systems but are restricted to moderately entangled states as the number of parameters scales exponentially with the entanglement entropy. While MPS can efficiently find ground states of 1D systems, their capacities are limited when simulating their dynamics, where the entanglement can increase ballistically with time. On the other hand, quantum devices appear as a natural platform to encode and perform the time evolution of correlated many-body states. However, accessing the regime of long-time dynamics is hampered by quantum noise. In this study we use the best of worlds: the short-time dynamics is efficiently performed by MPSs, compiled into short-depth quantum circuits, and is performed further in time on a quantum computer thanks to efficient MPO-optimized quantum circuits. We quantify the capacities of this hybrid classical-quantum scheme in terms of fidelities taking into account a noise model. We show that using classical knowledge in the form of tensor networks provides a way to better use limited quantum resources and lowers drastically the noise requirements to reach a practical quantum advantage. Finally we successfully demonstrate our approach with an experimental realization of the technique. Combined with efficient circuit transpilation we simulate a 10-qubit system on an actual quantum device over a longer time scale than low-bond-dimension MPSs and purely quantum Trotter evolution.
△ Less
Submitted 8 January, 2024; v1 submitted 30 May, 2023;
originally announced May 2023.
-
Quantum computing with and for many-body physics
Authors:
Thomas Ayral,
Pauline Besserve,
Denis Lacroix,
Edgar Andres Ruiz Guzman
Abstract:
Quantum computing technologies are making steady progress. This has opened new opportunities for tackling problems whose complexity prevents their description on classical computers. A prototypical example of these complex problems are interacting quantum many-body systems: on the one hand, these systems are known to become rapidly prohibitive to describe using classical computers when their size…
▽ More
Quantum computing technologies are making steady progress. This has opened new opportunities for tackling problems whose complexity prevents their description on classical computers. A prototypical example of these complex problems are interacting quantum many-body systems: on the one hand, these systems are known to become rapidly prohibitive to describe using classical computers when their size increases. On the other hand, these systems are precisely those which are used in the laboratory to build quantum computing platforms. This arguably makes them one of the most promising early use cases of quantum computing. In this review, we explain how quantum many-body systems are used to build quantum processors, and how, in turn, current and future quantum processors can be used to describe large many-body systems of fermions such as electrons and nucleons. The review includes an introduction to analog and digital quantum devices, the mapping of Fermi systems and their Hamiltonians onto qubit registers, as well as an overview of methods to access their static and dynamical properties. We also highlight some aspects related to entanglement, and touch on the description, influence and processing of decoherence in quantum devices.
△ Less
Submitted 27 September, 2023; v1 submitted 8 March, 2023;
originally announced March 2023.
-
A density-matrix renormalization group algorithm for simulating quantum circuits with a finite fidelity
Authors:
Thomas Ayral,
Thibaud Louvet,
Yiqing Zhou,
Cyprien Lambert,
E. Miles Stoudenmire,
Xavier Waintal
Abstract:
We develop a density-matrix renormalization group (DMRG) algorithm for the simulation of quantum circuits. This algorithm can be seen as the extension of time-dependent DMRG from the usual situation of hermitian Hamiltonian matrices to quantum circuits defined by unitary matrices. For small circuit depths, the technique is exact and equivalent to other matrix product state (MPS) based techniques.…
▽ More
We develop a density-matrix renormalization group (DMRG) algorithm for the simulation of quantum circuits. This algorithm can be seen as the extension of time-dependent DMRG from the usual situation of hermitian Hamiltonian matrices to quantum circuits defined by unitary matrices. For small circuit depths, the technique is exact and equivalent to other matrix product state (MPS) based techniques. For larger depths, it becomes approximate in exchange for an exponential speed up in computational time. Like an actual quantum computer, the quality of the DMRG results is characterized by a finite fidelity. However, unlike a quantum computer, the fidelity depends strongly on the quantum circuit considered. For the most difficult possible circuit for this technique, the so-called "quantum supremacy" benchmark of Google Inc. , we find that the DMRG algorithm can generate bit strings of the same quality as the seminal Google experiment on a single computing core. For a more structured circuit used for combinatorial optimization (Quantum Approximate Optimization Algorithm or QAOA), we find a drastic improvement of the DMRG results with error rates dropping by a factor of 100 compared with random quantum circuits. Our results suggest that the current bottleneck of quantum computers is their fidelities rather than the number of qubits.
△ Less
Submitted 29 August, 2022; v1 submitted 12 July, 2022;
originally announced July 2022.
-
Open Source Variational Quantum Eigensolver Extension of the Quantum Learning Machine (QLM) for Quantum Chemistry
Authors:
Mohammad Haidar,
Marko J. Rančić,
Thomas Ayral,
Yvon Maday,
Jean-Philip Piquemal
Abstract:
Quantum Chemistry (QC) is one of the most promising applications of Quantum Computing. However, present quantum processing units (QPUs) are still subject to large errors. Therefore, noisy intermediate-scale quantum (NISQ) hardware is limited in terms of qubits counts and circuit depths. Specific algorithms such as Variational Quantum Eigensolvers (VQEs) can potentially overcome such issues. We int…
▽ More
Quantum Chemistry (QC) is one of the most promising applications of Quantum Computing. However, present quantum processing units (QPUs) are still subject to large errors. Therefore, noisy intermediate-scale quantum (NISQ) hardware is limited in terms of qubits counts and circuit depths. Specific algorithms such as Variational Quantum Eigensolvers (VQEs) can potentially overcome such issues. We introduce here a novel open-source QC package, denoted Open-VQE, providing tools for using and developing chemically-inspired adaptive methods derived from Unitary Coupled Cluster (UCC). It facilitates the development and testing of VQE algorithms. It is able to use the Atos Quantum Learning Machine (QLM), a general quantum programming framework enabling to write, optimize and simulate quantum computing programs. Along with Open-VQE, we introduce myQLM-Fermion, a new open-source module (that includes the key QLM ressources that are important for QC developments (fermionic second quantization tools etc...). The Open-VQE package extends therefore QLM to QC providing: (i) the functions to generate the different types of excitations beyond the commonly used UCCSD ans{ä}tz;(ii) a new implementation of the "adaptive derivative assembled pseudo-Trotter method" (ADAPT-VQE), written in simple class structure python codes. Interoperability with other major quantum programming frameworks is ensured thanks to myQLM, which allows users to easily build their own code and execute it on existing QPUs. The combined Open-VQE/myQLM-Fermion quantum simulator facilitates the implementation, tests and developments of variational quantum algorithms towards choosing the best compromise to run QC computations on present quantum computers while offering the possibility to test large molecules. We provide extensive benchmarks for several molecules associated to qubit counts ranging from 4 up to 24.
△ Less
Submitted 28 November, 2022; v1 submitted 17 June, 2022;
originally announced June 2022.
-
Unraveling correlated material properties with noisy quantum computers: Natural orbitalized variational quantum eigensolving of extended impurity models within a slave-boson approach
Authors:
Pauline Besserve,
Thomas Ayral
Abstract:
We propose a method for computing space-resolved correlation properties of the two-dimensional Hubbard model within a quantum-classical embedding strategy that uses a Noisy, Intermediate Scale Quantum (NISQ) computer to solve the embedded model. While previous approaches were limited to purely local, one-impurity embedded models, requiring at most four qubits and relatively shallow circuits, we so…
▽ More
We propose a method for computing space-resolved correlation properties of the two-dimensional Hubbard model within a quantum-classical embedding strategy that uses a Noisy, Intermediate Scale Quantum (NISQ) computer to solve the embedded model. While previous approaches were limited to purely local, one-impurity embedded models, requiring at most four qubits and relatively shallow circuits, we solve a two-impurity model requiring eight qubits with an advanced hybrid scheme on top of the Variational Quantum Eigensolver algorithm. This iterative scheme, dubbed Natural Orbitalization (NOization), gradually transforms the single-particle basis to the approximate Natural-Orbital basis, in which the ground state can be minimally expressed, at the cost of measuring the one-particle reduced density-matrix of the embedded problem. We show that this transformation tends to make the variational optimization of existing (but too deep) ansatz circuits faster and more accurate, and we propose an ansatz, the Multireference Excitation Preserving (MREP) ansatz, that achieves great expressivity without requiring a prohibitive gate count. The one-impurity version of the ansatz has only one parameter, making the ground state preparation a trivial step, which supports the optimal character of our approach. Within a Rotationally Invariant Slave Boson embedding scheme that requires a minimal number of bath sites and does not require computing the full Green's function, the NOization combined with the MREP ansatz allow us to compute accurate, space-resolved quasiparticle weights and static self-energies for the Hubbard model even in the presence of noise levels representative of current NISQ processors. This paves the way to a controlled solution of the Hubbard model with larger and larger embedded problems solved by quantum computers.
△ Less
Submitted 1 March, 2022; v1 submitted 24 August, 2021;
originally announced August 2021.
-
Benchmarking quantum co-processors in an application-centric, hardware-agnostic and scalable way
Authors:
Simon Martiel,
Thomas Ayral,
Cyril Allouche
Abstract:
Existing protocols for benchmarking current quantum co-processors fail to meet the usual standards for assessing the performance of High-Performance-Computing platforms. After a synthetic review of these protocols -- whether at the gate, circuit or application level -- we introduce a new benchmark, dubbed Atos Q-score (TM), that is application-centric, hardware-agnostic and scalable to quantum adv…
▽ More
Existing protocols for benchmarking current quantum co-processors fail to meet the usual standards for assessing the performance of High-Performance-Computing platforms. After a synthetic review of these protocols -- whether at the gate, circuit or application level -- we introduce a new benchmark, dubbed Atos Q-score (TM), that is application-centric, hardware-agnostic and scalable to quantum advantage processor sizes and beyond. The Q-score measures the maximum number of qubits that can be used effectively to solve the MaxCut combinatorial optimization problem with the Quantum Approximate Optimization Algorithm. We give a robust definition of the notion of effective performance by introducing an improved approximation ratio based on the scaling of random and optimal algorithms. We illustrate the behavior of Q-score using perfect and noisy simulations of quantum processors. Finally, we provide an open-source implementation of Q-score that makes it easy to compute the Q-score of any quantum hardware.
△ Less
Submitted 26 July, 2021; v1 submitted 25 February, 2021;
originally announced February 2021.
-
Quantum Divide and Compute: Exploring The Effect of Different Noise Sources
Authors:
Thomas Ayral,
François-Marie Le Régent,
Zain Saleem,
Yuri Alexeev,
Martin Suchara
Abstract:
Our recent work (Ayral et al., 2020 IEEE Computer Society Annual Symposium on VLSI (ISVLSI)) showed the first implementation of the Quantum Divide and Compute (QDC) method, which allows to break quantum circuits into smaller fragments with fewer qubits and shallower depth. QDC can thus deal with the limited number of qubits and short coherence times of noisy, intermediate-scale quantum processors.…
▽ More
Our recent work (Ayral et al., 2020 IEEE Computer Society Annual Symposium on VLSI (ISVLSI)) showed the first implementation of the Quantum Divide and Compute (QDC) method, which allows to break quantum circuits into smaller fragments with fewer qubits and shallower depth. QDC can thus deal with the limited number of qubits and short coherence times of noisy, intermediate-scale quantum processors. This article investigates the impact of different noise sources -- readout error, gate error and decoherence -- on the success probability of the QDC procedure. We perform detailed noise modeling on the Atos Quantum Learning Machine, allowing us to understand tradeoffs and formulate recommendations about which hardware noise sources should be preferentially optimized. We describe in detail the noise models we used to reproduce experimental runs on IBM's Johannesburg processor. This work also includes a detailed derivation of the equations used in the QDC procedure to compute the output distribution of the original quantum circuit from the output distribution of its fragments. Finally, we analyze the computational complexity of the QDC method for the circuit under study via tensor-network considerations, and elaborate on the relation the QDC method with tensor-network simulation methods.
△ Less
Submitted 7 February, 2021;
originally announced February 2021.
-
Temporal Stochastic Softmax for 3D CNNs: An Application in Facial Expression Recognition
Authors:
Théo Ayral,
Marco Pedersoli,
Simon Bacon,
Eric Granger
Abstract:
Training deep learning models for accurate spatiotemporal recognition of facial expressions in videos requires significant computational resources. For practical reasons, 3D Convolutional Neural Networks (3D CNNs) are usually trained with relatively short clips randomly extracted from videos. However, such uniform sampling is generally sub-optimal because equal importance is assigned to each tempo…
▽ More
Training deep learning models for accurate spatiotemporal recognition of facial expressions in videos requires significant computational resources. For practical reasons, 3D Convolutional Neural Networks (3D CNNs) are usually trained with relatively short clips randomly extracted from videos. However, such uniform sampling is generally sub-optimal because equal importance is assigned to each temporal clip. In this paper, we present a strategy for efficient video-based training of 3D CNNs. It relies on softmax temporal pooling and a weighted sampling mechanism to select the most relevant training clips. The proposed softmax strategy provides several advantages: a reduced computational complexity due to efficient clip sampling, and an improved accuracy since temporal weighting focuses on more relevant clips during both training and inference. Experimental results obtained with the proposed method on several facial expression recognition benchmarks show the benefits of focusing on more informative clips in training videos. In particular, our approach improves performance and computational cost by reducing the impact of inaccurate trimming and coarse annotation of videos, and heterogeneous distribution of visual information across time.
△ Less
Submitted 10 November, 2020;
originally announced November 2020.
-
Solving optimization problems with Rydberg analog quantum computers: Realistic requirements for quantum advantage using noisy simulation and classical benchmarks
Authors:
Michel Fabrice Serret,
Bertrand Marchand,
Thomas Ayral
Abstract:
Platforms of Rydberg atoms have been proposed as promising candidates to solve some combinatorial optimization problems. Here, we compute quantitative requirements on the system sizes and noise levels that these platforms must fulfill to reach quantum advantage in approximately solving the Unit-Disk Maximum Independent Set problem. Using noisy simulations of Rydberg platforms of up to 26 atoms int…
▽ More
Platforms of Rydberg atoms have been proposed as promising candidates to solve some combinatorial optimization problems. Here, we compute quantitative requirements on the system sizes and noise levels that these platforms must fulfill to reach quantum advantage in approximately solving the Unit-Disk Maximum Independent Set problem. Using noisy simulations of Rydberg platforms of up to 26 atoms interacting through realistic van der Waals interactions, we compute the average approximation ratio that can be attained with a simple quantum annealing-based heuristic within a fixed temporal computational budget. Based on estimates of the correlation lengths measured in the engineered quantum state, we extrapolate the results to large atom numbers and compare them to a simple classical approximation heuristic. We find that approximation ratios of at least $\approx 0.84$ are within reach for near-future noise levels. Not taking into account further classical and quantum algorithmic improvements, we estimate that quantum advantage could be reached by attaining a number of controlled atoms of $\sim8,000$ for a time budget of 2 seconds, and $\sim 1,000-1,200$ for a time budget of 0.2 seconds, provided the coherence levels of the system can be improved by a factor 10 while maintaining a constant repetition rate.
△ Less
Submitted 5 November, 2020; v1 submitted 19 June, 2020;
originally announced June 2020.
-
Quantum Divide and Compute: Hardware Demonstrations and Noisy Simulations
Authors:
Thomas Ayral,
François-Marie Le Régent,
Zain Saleem,
Yuri Alexeev,
Martin Suchara
Abstract:
Noisy, intermediate-scale quantum computers come with intrinsic limitations in terms of the number of qubits (circuit "width") and decoherence time (circuit "depth") they can have. Here, for the first time, we demonstrate a recently introduced method that breaks a circuit into smaller subcircuits or fragments, and thus makes it possible to run circuits that are either too wide or too deep for a gi…
▽ More
Noisy, intermediate-scale quantum computers come with intrinsic limitations in terms of the number of qubits (circuit "width") and decoherence time (circuit "depth") they can have. Here, for the first time, we demonstrate a recently introduced method that breaks a circuit into smaller subcircuits or fragments, and thus makes it possible to run circuits that are either too wide or too deep for a given quantum processor. We investigate the behavior of the method on one of IBM's 20-qubit superconducting quantum processors with various numbers of qubits and fragments. We build noise models that capture decoherence, readout error, and gate imperfections for this particular processor. We then carry out noisy simulations of the method in order to account for the observed experimental results. We find an agreement within 20% between the experimental and the simulated success probabilities, and we observe that recombining noisy fragments yields overall results that can outperform the results without fragmentation.
△ Less
Submitted 26 May, 2020;
originally announced May 2020.
-
Rotationally invariant slave-boson and density matrix embedding theory: A unified framework and a comparative study on the 1D and 2D Hubbard Model
Authors:
Tsung-Han Lee,
Thomas Ayral,
Yong-Xin Yao,
Nicola Lanata,
Gabriel Kotliar
Abstract:
We present detailed benchmark ground-state calculations of the one- and two-dimensional Hubbard model utilizing the cluster extensions of the rotationally invariant slave-boson (RISB) mean-field theory and the density matrix embedding theory (DMET). Our analysis shows that the overall accuracy and the performance of these two methods are very similar. Furthermore, we propose a unified computationa…
▽ More
We present detailed benchmark ground-state calculations of the one- and two-dimensional Hubbard model utilizing the cluster extensions of the rotationally invariant slave-boson (RISB) mean-field theory and the density matrix embedding theory (DMET). Our analysis shows that the overall accuracy and the performance of these two methods are very similar. Furthermore, we propose a unified computational framework that allows us to implement both of these techniques on the same footing. This provides us with a new line of interpretation and paves the ways for developing systematically new generalizations of these complementary approaches.
△ Less
Submitted 23 December, 2018;
originally announced December 2018.
-
Dynamical Mean Field Theory, Density-Matrix Embedding Theory and Rotationally Invariant Slave Bosons: a Unified Perspective
Authors:
Thomas Ayral,
Tsung-Han Lee,
Gabriel Kotliar
Abstract:
We present a unified perspective on Dynamical Mean Field Theory (DMFT), Density-Matrix Embedding Theory (DMET) and Rotationally Invariant Slave Bosons (RISB). We show that DMET can be regarded as a simplification of the RISB method where the quasiparticle weight is set to unity. This relation allows to easily transpose extensions of a given method to another: for instance, a temperature-dependent…
▽ More
We present a unified perspective on Dynamical Mean Field Theory (DMFT), Density-Matrix Embedding Theory (DMET) and Rotationally Invariant Slave Bosons (RISB). We show that DMET can be regarded as a simplification of the RISB method where the quasiparticle weight is set to unity. This relation allows to easily transpose extensions of a given method to another: for instance, a temperature-dependent version of RISB can be used to derive a temperature-dependent free-energy formula for DMET.
△ Less
Submitted 4 December, 2017; v1 submitted 21 October, 2017;
originally announced October 2017.
-
Chiral d-wave Superconductivity in a Triangular Surface Lattice Mediated by Long-range Interaction
Authors:
Xiaodong Cao,
Thomas Ayral,
Zhicheng Zhong,
Olivier Parcollet,
Dirk Manske,
Philipp Hansmann
Abstract:
Correlated ad-atom systems on the Si(111) surface have recently attracted an increased attention as strongly correlated systems with a rich phase diagram. We study these materials by a single band model on the triangular lattice including 1/r long-range interaction. Employing the recently proposed TRILEX method we find an unconventional superconducting phase of chiral d-wave symmetry in hole-doped…
▽ More
Correlated ad-atom systems on the Si(111) surface have recently attracted an increased attention as strongly correlated systems with a rich phase diagram. We study these materials by a single band model on the triangular lattice including 1/r long-range interaction. Employing the recently proposed TRILEX method we find an unconventional superconducting phase of chiral d-wave symmetry in hole-doped systems. The superconductivity is driven simultaneously by both charge and spin fluctuations and is strongly enhanced by the long-range tail of the interaction. We provide an analysis of the relevant collective bosonic modes and explain how in triangular symmetry both charge and spin channels contribute to the Cooper-pairing.
△ Less
Submitted 10 October, 2017;
originally announced October 2017.
-
Non-local Coulomb interactions on the triangular lattice in the high-doping regime: Spectra and charge dynamics from Extended Dynamical Mean Field Theory
Authors:
Sophie Chauvin,
Thomas Ayral,
Lucia Reining,
Silke Biermann
Abstract:
We explore the two-dimensional extended Hubbard model on the triangular lattice in the high doping regime. On-site and nearest-neighbour repulsive interactions are treated in a non-perturbative way by means of Extended Dynamical Mean Field Theory. We compute the low-temperature phase diagram, displaying a metallic phase and a symmetry-broken phase for strong intersite repulsions. We describe the c…
▽ More
We explore the two-dimensional extended Hubbard model on the triangular lattice in the high doping regime. On-site and nearest-neighbour repulsive interactions are treated in a non-perturbative way by means of Extended Dynamical Mean Field Theory. We compute the low-temperature phase diagram, displaying a metallic phase and a symmetry-broken phase for strong intersite repulsions. We describe the correlation effects on both single-particle and two-particle observables in the metallic phase. Whereas single-particle spectra feature a Hubbard satellite typical of strongly correlated systems, local susceptibilities remain close to their non-interacting limit, even for large on-site repulsions. We argue that this behaviour is typical of the strongly doped case. We also report a region in parameter space with negative static local screening.
△ Less
Submitted 22 September, 2017;
originally announced September 2017.
-
The Fierz convergence criterion: a controlled approach to strongly-interacting systems with small embedded clusters
Authors:
Thomas Ayral,
Jaksa Vucicevic,
Olivier Parcollet
Abstract:
We present an embedded-cluster method, based on the TRILEX formalism. It turns the Fierz ambiguity, inherent to approaches based on a bosonic decoupling of local fermionic interactions, into a convergence criterion. It is based on the approximation of the three-leg vertex by a coarse-grained vertex computed by solving a self-consistently determined multi-site effective impurity model. The computed…
▽ More
We present an embedded-cluster method, based on the TRILEX formalism. It turns the Fierz ambiguity, inherent to approaches based on a bosonic decoupling of local fermionic interactions, into a convergence criterion. It is based on the approximation of the three-leg vertex by a coarse-grained vertex computed by solving a self-consistently determined multi-site effective impurity model. The computed self-energies are, by construction, continuous functions of momentum. We show that, in three interaction and doping regimes of parameters of the two-dimensional Hubbard model, self-energies obtained with clusters of size four only are very close to numerically exact benchmark results. We show that the Fierz parameter, which parametrizes the freedom in the Hubbard-Stratonovich decoupling, can be used as a quality control parameter. By contrast, the $GW$+extended dynamical mean field theory approximation with four cluster sites is shown to yield good results only in the weak-coupling regime and for a particular decoupling. Finally, we show that the vertex has spatially nonlocal components only at low Matsubara frequencies.
△ Less
Submitted 27 July, 2017; v1 submitted 5 June, 2017;
originally announced June 2017.
-
TRILEX and $GW$+EDMFT approach to $d$-wave superconductivity in the Hubbard model
Authors:
Jaksa Vucicevic,
Thomas Ayral,
Olivier Parcollet
Abstract:
We generalize the recently introduced TRILEX approach (TRiply Irreducible Local EXpansion) to superconducting phases. The method treats simultaneously Mott and spin-fluctuation physics using an Eliashberg theory supplemented by local vertex corrections determined by a self-consistent quantum impurity model. We show that, in the two-dimensional Hubbard model, at strong coupling, TRILEX yields a…
▽ More
We generalize the recently introduced TRILEX approach (TRiply Irreducible Local EXpansion) to superconducting phases. The method treats simultaneously Mott and spin-fluctuation physics using an Eliashberg theory supplemented by local vertex corrections determined by a self-consistent quantum impurity model. We show that, in the two-dimensional Hubbard model, at strong coupling, TRILEX yields a $d$-wave superconducting dome as a function of doping. Contrary to the standard cluster dynamical mean field theory (DMFT) approaches, TRILEX can capture $d$-wave pairing using only a single-site effective impurity model. We also systematically explore the dependence of the superconducting temperature on the bare dispersion at weak coupling, which shows a clear link between strong antiferromagnetic (AF) correlations and the onset of superconductivity. We identify a combination of hopping amplitudes particularly favorable to superconductivity at intermediate doping. Finally, we study within $GW$+EDMFT the low-temperature $d$-wave superconducting phase at strong coupling in a region of parameter space with reduced AF fluctuations.
△ Less
Submitted 31 July, 2017; v1 submitted 23 May, 2017;
originally announced May 2017.
-
Influence of Fock exchange in combined many-body perturbation and dynamical mean field theory
Authors:
Thomas Ayral,
Silke Biermann,
Philipp Werner,
Lewin Boehnke
Abstract:
In electronic systems with long-range Coulomb interaction, the nonlocal Fock exchange term has a band-widening effect. While this effect is included in combined many-body perturbation theory and dynamical mean field theory schemes, it is not taken into account in standard extended DMFT (EDMFT) calculations. Here, we include this instantaneous term in both approaches and investigate its effect on t…
▽ More
In electronic systems with long-range Coulomb interaction, the nonlocal Fock exchange term has a band-widening effect. While this effect is included in combined many-body perturbation theory and dynamical mean field theory schemes, it is not taken into account in standard extended DMFT (EDMFT) calculations. Here, we include this instantaneous term in both approaches and investigate its effect on the phase diagram and dynamically screened interaction. We show that the largest deviations between previously presented EDMFT and $GW$+EDMFT results originate from the nonlocal Fock term, and that the quantitative differences are especially large in the strong-coupling limit. Furthermore, we show that the charge-ordering phase diagram obtained in $GW$+EDMFT methods for moderate interaction values is very similar to the one predicted by dual boson methods that include the fermion-boson or four-point vertex.
△ Less
Submitted 5 May, 2017; v1 submitted 26 January, 2017;
originally announced January 2017.
-
Large effects of subtle electronic correlations on the energetics of vacancies in alpha-Fe
Authors:
Pascal Delange,
Thomas Ayral,
Sergei I. Simak,
Michel Ferrero,
Olivier Parcollet,
Silke Biermann,
Leonid Pourovskii
Abstract:
We study the effect of electronic Coulomb correlations on the vacancy formation energy in paramagnetic alpha-Fe within ab initio dynamical mean-field theory. The calculated value for the formation energy is substantially lower than in standard density-functional calculations and in excellent agreement with experiment. The reduction is caused by an enhancement of electronic correlations at the near…
▽ More
We study the effect of electronic Coulomb correlations on the vacancy formation energy in paramagnetic alpha-Fe within ab initio dynamical mean-field theory. The calculated value for the formation energy is substantially lower than in standard density-functional calculations and in excellent agreement with experiment. The reduction is caused by an enhancement of electronic correlations at the nearest neighbors of the vacancy. This effect is explained by subtle changes in the corresponding spectral function of the d-electrons. The local lattice relaxations around the vacancy are substantially increased by many-body effects.
△ Less
Submitted 16 June, 2016;
originally announced June 2016.
-
Mott physics and collective modes: an atomic approximation of the four-particle irreducible functional
Authors:
Thomas Ayral,
Olivier Parcollet
Abstract:
We discuss a generalization of the dynamical mean field theory (DMFT) for strongly correlated systems close to a Mott transition based on a systematic approximation of the fully irreducible four-point vertex. It is an atomic-limit approximation of a functional of the one- and two-particle Green functions, built with the second Legendre transform of the free energy with respect to the two-particle…
▽ More
We discuss a generalization of the dynamical mean field theory (DMFT) for strongly correlated systems close to a Mott transition based on a systematic approximation of the fully irreducible four-point vertex. It is an atomic-limit approximation of a functional of the one- and two-particle Green functions, built with the second Legendre transform of the free energy with respect to the two-particle Green function. This functional is represented diagrammatically by four-particle irreducible (4PI) diagrams. Like the dynamical vertex approximation (D$Γ$A), the fully irreducible vertex is computed from a quantum impurity model whose bath is self-consistently determined by solving the parquet equations. However, in contrast with D$Γ$A and DMFT, the interaction term of the impurity model is also self-consistently determined. The method interpolates between the parquet approximation at weak coupling and the atomic limit, where it is exact. It is applicable to systems with short-range and long-range interactions.
△ Less
Submitted 4 August, 2016; v1 submitted 29 May, 2016;
originally announced May 2016.
-
Mott physics and spin fluctuations: a functional viewpoint
Authors:
Thomas Ayral,
Olivier Parcollet
Abstract:
We present a formalism for strongly correlated systems with fermions coupled to bosonic modes. We construct the three-particle irreducible functional $\mathcal{K}$ by successive Legendre transformations of the free energy of the system. We derive a closed set of equations for the fermionic and bosonic self-energies for a given $\mathcal{K}$. We then introduce a local approximation for…
▽ More
We present a formalism for strongly correlated systems with fermions coupled to bosonic modes. We construct the three-particle irreducible functional $\mathcal{K}$ by successive Legendre transformations of the free energy of the system. We derive a closed set of equations for the fermionic and bosonic self-energies for a given $\mathcal{K}$. We then introduce a local approximation for $\mathcal{K}$, which extends the idea of dynamical mean field theory (DMFT) approaches from two- to three-particle irreducibility. This approximation entails the locality of the three-leg electron-boson vertex $Λ(iω,iΩ)$, which is self-consistently computed using a quantum impurity model with dynamical charge and spin interactions. This local vertex is used to construct frequency- and momentum-dependent electronic self-energies and polarizations. By construction, the method interpolates between the spin-fluctuation or GW approximations at weak coupling and the atomic limit at strong coupling. We apply it to the Hubbard model on two-dimensional square and triangular lattices. We complement the results of Phys.Rev. B 92, 115109 by (i) showing that, at half-filling, as DMFT, the method describes the Fermi-liquid metallic state and the Mott insulator, separated by a first-order interacting-driven Mott transition at low temperatures, (ii) investigating the influence of frustration and (iii) discussing the influence of the bosonic decoupling channel.
△ Less
Submitted 31 March, 2016; v1 submitted 21 December, 2015;
originally announced December 2015.
-
Uncertainty principle for experimental measurements: Fast versus slow probes
Authors:
Philipp Hansmann,
Thomas Ayral,
Antonio Tejeda,
Silke Biermann
Abstract:
The result of a physical measurement depends on the timescale of the experimental probe. In solid-state systems, this simple quantum mechanical principle has far-reaching consequences: the interplay of several degrees of freedom close to charge, spin or orbital instabilities combined with the disparity of the time scales associated to their fluctuations can lead to seemingly contradictory experime…
▽ More
The result of a physical measurement depends on the timescale of the experimental probe. In solid-state systems, this simple quantum mechanical principle has far-reaching consequences: the interplay of several degrees of freedom close to charge, spin or orbital instabilities combined with the disparity of the time scales associated to their fluctuations can lead to seemingly contradictory experimental findings. A particularly striking example is provided by systems of adatoms adsorbed on semiconductor surfaces where different experiments -- angle-resolved photoemission, scanning tunneling microscopy and core-level spectroscopy -- suggest different ordering phenomena. Using most recent first principles many-body techniques, we resolve this puzzle by invoking the time scales of fluctuations when approaching the different instabilities. These findings suggest a re-interpretation of ordering phenomena and their fluctuations in a wide class of solid-state systems ranging from organic materials to high-temperature superconducting cuprates.
△ Less
Submitted 16 November, 2015;
originally announced November 2015.
-
TRIQS: A Toolbox for Research on Interacting Quantum Systems
Authors:
Olivier Parcollet,
Michel Ferrero,
Thomas Ayral,
Hartmut Hafermann,
Igor Krivenko,
Laura Messio,
Priyanka Seth
Abstract:
We present the TRIQS library, a Toolbox for Research on Interacting Quantum Systems. It is an open-source, computational physics library providing a framework for the quick development of applications in the field of many-body quantum physics, and in particular, strongly-correlated electronic systems. It supplies components to develop codes in a modern, concise and efficient way: e.g. Green's func…
▽ More
We present the TRIQS library, a Toolbox for Research on Interacting Quantum Systems. It is an open-source, computational physics library providing a framework for the quick development of applications in the field of many-body quantum physics, and in particular, strongly-correlated electronic systems. It supplies components to develop codes in a modern, concise and efficient way: e.g. Green's function containers, a generic Monte Carlo class, and simple interfaces to HDF5. TRIQS is a C++/Python library that can be used from either language. It is distributed under the GNU General Public License (GPLv3). State-of-the-art applications based on the library, such as modern quantum many-body solvers and interfaces between density-functional-theory codes and dynamical mean-field theory (DMFT) codes are distributed along with it.
△ Less
Submitted 21 July, 2015; v1 submitted 8 April, 2015;
originally announced April 2015.
-
Mott physics and spin fluctuations: a unified framework
Authors:
Thomas Ayral,
Olivier Parcollet
Abstract:
We present a formalism for strongly correlated electrons systems which consists in a local approximation of the dynamical three-leg interaction vertex. This vertex is self-consistently computed with a quantum impurity model with dynamical interactions in the charge and spin channels, similar to dynamical mean field theory (DMFT) approaches. The electronic self-energy and the polarization are both…
▽ More
We present a formalism for strongly correlated electrons systems which consists in a local approximation of the dynamical three-leg interaction vertex. This vertex is self-consistently computed with a quantum impurity model with dynamical interactions in the charge and spin channels, similar to dynamical mean field theory (DMFT) approaches. The electronic self-energy and the polarization are both frequency and momentum dependent. The method interpolates between the spin-fluctuation or GW approximations at weak coupling and the atomic limit at strong coupling. We apply the formalism to the Hubbard model on a two-dimensional square lattice and show that as interactions are increased towards the Mott insulating state, the local vertex acquires a strong frequency dependence, driving the system to a Mott transition, while at low enough temperatures the momentum-dependence of the self-energy is enhanced due to large spin fluctuations. Upon doping, we find a Fermi arc in the one-particle spectral function, which is one signature of the pseudo-gap state.
△ Less
Submitted 26 March, 2015;
originally announced March 2015.
-
Dynamical correlations and screened exchange on the experimental bench: spectral properties of the cobalt pnictide BaCo2As2
Authors:
Ambroise van Roekeghem,
Thomas Ayral,
Jan M. Tomczak,
Michele Casula,
Nan Xu,
Hong Ding,
Michel Ferrero,
Olivier Parcollet,
Hong Jiang,
Silke Biermann
Abstract:
Understanding the Fermi surface and low-energy excitations of iron or cobalt pnictides is crucial for assessing electronic instabilities such as magnetic or superconducting states. Here, we propose and implement a new approach to compute the low-energy properties of correlated electron materials, taking into account both screened exchange beyond the local density approximation and local dynamical…
▽ More
Understanding the Fermi surface and low-energy excitations of iron or cobalt pnictides is crucial for assessing electronic instabilities such as magnetic or superconducting states. Here, we propose and implement a new approach to compute the low-energy properties of correlated electron materials, taking into account both screened exchange beyond the local density approximation and local dynamical correlations. The scheme allows us to resolve the puzzle of BaCo2As2, for which standard electronic structure techniques predict a ferromagnetic instability not observed in nature.
△ Less
Submitted 27 December, 2014; v1 submitted 13 August, 2014;
originally announced August 2014.
-
Extended dynamical mean-field study of the Hubbard model with long range interactions
Authors:
Li Huang,
Thomas Ayral,
Silke Biermann,
Philipp Werner
Abstract:
Using extended dynamical mean-field theory and its combination with the $GW$ approximation, we compute the phase diagrams and local spectral functions of the single-band extended Hubbard model on the square and simple cubic lattices, considering long range interactions up to the third nearest neighbors. The longer range interactions shift the boundaries between the metallic, charge-ordered insulat…
▽ More
Using extended dynamical mean-field theory and its combination with the $GW$ approximation, we compute the phase diagrams and local spectral functions of the single-band extended Hubbard model on the square and simple cubic lattices, considering long range interactions up to the third nearest neighbors. The longer range interactions shift the boundaries between the metallic, charge-ordered insulating and Mott insulating phases, and lead to characteristic changes in the screening modes and local spectral functions. Momentum-dependent self-energy contributions enhance the correlation effects and thus compete with the additional screening effect from longer range Coulomb interactions. Our results suggest that the influence of longer range intersite interactions is significant, and that these effects deserve attention in realistic studies of correlated materials.
△ Less
Submitted 17 October, 2014; v1 submitted 28 April, 2014;
originally announced April 2014.
-
Long-range Coulomb interactions in surface systems: a first principles description within self-consistently combined GW and dynamical mean field theory
Authors:
Philipp Hansmann,
Thomas Ayral,
Loig Vaugier,
Philipp Werner,
Silke Biermann
Abstract:
Systems of adatoms on semiconductor surfaces display competing ground states and exotic spectral properties typical of two-dimensional correlated electron materials which are dominated by a complex interplay of spin and charge degrees of freedom. We report a fully ab initio derivation of low energy Hamiltonians for the adatom systems Si(111):X, with X=Sn, Si, C, Pb, that we solve within self-consi…
▽ More
Systems of adatoms on semiconductor surfaces display competing ground states and exotic spectral properties typical of two-dimensional correlated electron materials which are dominated by a complex interplay of spin and charge degrees of freedom. We report a fully ab initio derivation of low energy Hamiltonians for the adatom systems Si(111):X, with X=Sn, Si, C, Pb, that we solve within self-consistent combined GW and dynamical mean field theory ("GW+DMFT"). Calculated photoemission spectra are in agreement with available experimental data. We rationalize experimentally observed tendencies from Mott physics towards charge-ordering along the series as resulting from substantial long-range interactions.
△ Less
Submitted 18 January, 2013;
originally announced January 2013.
-
Screening and Non-local Correlations in the Extended Hubbard Model from Self-Consistent Combined GW and Dynamical Mean Field Theory
Authors:
Thomas Ayral,
Silke Biermann,
Philipp Werner
Abstract:
We describe a recent implementation of the combined GW and dynamical mean field (DMFT) method "GW+DMFT" for the two-dimensional Hubbard model with on-site and nearest-neighbor repulsion. We clarify the relation of the GW+DMFT scheme to alternative approaches in the literature, and discuss the corresponding approximations to the free energy functional of the model. Furthermore, we describe a numeri…
▽ More
We describe a recent implementation of the combined GW and dynamical mean field (DMFT) method "GW+DMFT" for the two-dimensional Hubbard model with on-site and nearest-neighbor repulsion. We clarify the relation of the GW+DMFT scheme to alternative approaches in the literature, and discuss the corresponding approximations to the free energy functional of the model. Furthermore, we describe a numerically exact technique for the self-consistent solution of the GW+DMFT equations, namely the hybridization expansion continuous-time algorithm for solving the dynamical impurity model that arises within the GW+DMFT scheme. We compute the low-temperature phase diagram of the extended Hubbard model, addressing the metal-insulator transition at small intersite interactions and the transition to a charge-ordered state for stronger intersite repulsions. Within the GW+DMFT framework, as in extended DMFT, the intersite repulsion translates into a frequency-dependence of the local effective interaction. We analyze this dependence by extracting a charateristic plasma frequency and show how it affects the local spectral function.
△ Less
Submitted 9 October, 2012;
originally announced October 2012.
-
Spectral Properties of Correlated Materials: Local Vertex and Non-Local Two-Particle Correlations from Combined GW and Dynamical Mean Field Theory
Authors:
Thomas Ayral,
Philipp Werner,
Silke Biermann
Abstract:
We present a fully self-consistent combined GW and dynamical mean field (GW+DMFT) study of the spectral properties of the extended two-dimensional Hubbard model. The inclusion of the local dynamical vertex stemming from the DMFT self-energy and polarization is shown to cure the problems of self-consistent GW in the description of spectral properties. We calculate the momentum-resolved spectral fun…
▽ More
We present a fully self-consistent combined GW and dynamical mean field (GW+DMFT) study of the spectral properties of the extended two-dimensional Hubbard model. The inclusion of the local dynamical vertex stemming from the DMFT self-energy and polarization is shown to cure the problems of self-consistent GW in the description of spectral properties. We calculate the momentum-resolved spectral functions, the two-particle polarization and electron loss spectra, and show that the inclusion of GW in extended DMFT leads to a narrowing of the quasi-particle width and more pronounced Hubbard bands in the metallic regime as one approaches the charge-ordering transition. Finally, the momentum-dependence introduced by GW into the extended DMFT description of collective modes is found to affect their shape, giving rise to dispersive plasmon-like long-wavelength and stripe modes.
△ Less
Submitted 24 May, 2012;
originally announced May 2012.