Source author record

Roman Orus

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

37works
16topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

37 published item(s)

preprint2026arXiv

Fast Tensorization of Neural Networks via Slice-wise Feature Distillation

We propose a scalable tensorization framework for neural network compression based on slice-wise feature distillation. Unlike conventional tensor decomposition methods that rely on costly global finetuning, our approach decomposes the network into slices consisting of either individual layers or blocks (e.g., convolutional layers or MLPs), or small groups of consecutive layers, and tensorizes each slice independently to reproduce the intermediate representations of the original pretrained model. This modular strategy improves accuracy recovery, reduces data requirements, and enables efficient parallel optimization. Experiments on ResNet-34 show significant gains over conventional global tensorization, achieving near-lossless compression at moderate compression rates with faster optimization. Results on GPT-2 XL further demonstrate the scalability of the method and its applicability to large-scale models, particularly in distributed settings.

preprint2026arXiv

Quantum-enhanced Large Language Models on Quantum Hardware via Cayley Unitary Adapters

Large language models (LLMs) have transformed artificial intelligence, yet classical architectures impose a fundamental constraint: every trainable parameter demands classical memory that scales unfavourably with model size. Quantum computing offers a qualitatively different pathway, but practical demonstrations on real hardware have remained elusive for models of practical relevance. Here we show that Cayley-parameterised unitary adapters -- quantum circuit blocks inserted into the frozen projection layers of pre-trained LLMs and executed on a 156-qubit IBM Quantum System Two superconducting processor -- improve the perplexity of Llama 3.1 8B, an 8-billion-parameter model in widespread use, by 1.4% with only 6,000 additional parameters and end-to-end inference validated on real Quantum Processing Unit (QPU). A systematic study on SmolLM2 (135M parameters), chosen for its tractability, reveals monotonically improving perplexity with unitary block dimension, 83% recovery of compression-induced degradation, and correct answers to questions that both classical baselines fail -- with a sharp noise-expressivity phase transition identifying the concrete path to quantum utility at larger qubit scales.

preprint2024arXiv

Quantum artificial vision for defect detection in manufacturing

In this paper we consider several algorithms for quantum computer vision using Noisy Intermediate-Scale Quantum (NISQ) devices, and benchmark them for a real problem against their classical counterparts. Specifically, we consider two approaches: a quantum Support Vector Machine (QSVM) on a universal gate-based quantum computer, and QBoost on a quantum annealer. The quantum vision systems are benchmarked for an unbalanced dataset of images where the aim is to detect defects in manufactured car pieces. We see that the quantum algorithms outperform their classical counterparts in several ways, with QBoost allowing for larger problems to be analyzed with present-day quantum annealers. Data preprocessing, including dimensionality reduction and contrast enhancement, is also discussed, as well as hyperparameter tuning in QBoost. To the best of our knowledge, this is the first implementation of quantum computer vision systems for a problem of industrial relevance in a manufacturing production line.

preprint2024arXiv

Variational Quantum and Quantum-Inspired Clustering

Here we present a quantum algorithm for clustering data based on a variational quantum circuit. The algorithm allows to classify data into many clusters, and can easily be implemented in few-qubit Noisy Intermediate-Scale Quantum (NISQ) devices. The idea of the algorithm relies on reducing the clustering problem to an optimization, and then solving it via a Variational Quantum Eigensolver (VQE) combined with non-orthogonal qubit states. In practice, the method uses maximally-orthogonal states of the target Hilbert space instead of the usual computational basis, allowing for a large number of clusters to be considered even with few qubits. We benchmark the algorithm with numerical simulations using real datasets, showing excellent performance even with one single qubit. Moreover, a tensor network simulation of the algorithm implements, by construction, a quantum-inspired clustering algorithm that can run on current classical hardware.

preprint2022arXiv

Financial Index Tracking via Quantum Computing with Cardinality Constraints

In this work, we demonstrate how to apply non-linear cardinality constraints, important for real-world asset management, to quantum portfolio optimization. This enables us to tackle non-convex portfolio optimization problems using quantum annealing that would otherwise be challenging for classical algorithms. Being able to use cardinality constraints for portfolio optimization opens the doors to new applications for creating innovative portfolios and exchange-traded-funds (ETFs). We apply the methodology to the practical problem of enhanced index tracking and are able to construct smaller portfolios that significantly outperform the risk profile of the target index whilst retaining high degrees of tracking.

preprint2022arXiv

Language Design as Information Renormalization

Here we consider some well-known facts in syntax from a physics perspective, allowing us to establish equivalences between both fields with many consequences. Mainly, we observe that the operation MERGE, put forward by N. Chomsky in 1995, can be interpreted as a physical information coarse-graining. Thus, MERGE in linguistics entails information renormalization in physics, according to different time scales. We make this point mathematically formal in terms of language models. In this setting, MERGE amounts to a probability tensor implementing a coarse-graining, akin to a probabilistic context-free grammar. The probability vectors of meaningful sentences are given by stochastic tensor networks (TN) built from diagonal tensors and which are mostly loop-free, such as Tree Tensor Networks and Matrix Product States, thus being computationally very efficient to manipulate. We show that this implies the polynomially-decaying (long-range) correlations experimentally observed in language, and also provides arguments in favour of certain types of neural networks for language processing. Moreover, we show how to obtain such language models from quantum states that can be efficiently prepared on a quantum computer, and use this to find bounds on the perplexity of the probability distribution of words in a sentence. Implications of our results are discussed across several ambits.

preprint2022arXiv

Quantum-Inspired Tensor Neural Networks for Partial Differential Equations

Partial Differential Equations (PDEs) are used to model a variety of dynamical systems in science and engineering. Recent advances in deep learning have enabled us to solve them in a higher dimension by addressing the curse of dimensionality in new ways. However, deep learning methods are constrained by training time and memory. To tackle these shortcomings, we implement Tensor Neural Networks (TNN), a quantum-inspired neural network architecture that leverages Tensor Network ideas to improve upon deep learning approaches. We demonstrate that TNN provide significant parameter savings while attaining the same accuracy as compared to the classical Dense Neural Network (DNN). In addition, we also show how TNN can be trained faster than DNN for the same accuracy. We benchmark TNN by applying them to solve parabolic PDEs, specifically the Black-Scholes-Barenblatt equation, widely used in financial pricing theory, empirically showing the advantages of TNN over DNN. Further examples, such as the Hamilton-Jacobi-Bellman equation, are also discussed.

preprint2021arXiv

Dynamic Portfolio Optimization with Real Datasets Using Quantum Processors and Quantum-Inspired Tensor Networks

In this paper we tackle the problem of dynamic portfolio optimization, i.e., determining the optimal trading trajectory for an investment portfolio of assets over a period of time, taking into account transaction costs and other possible constraints. This problem is central to quantitative finance. After a detailed introduction to the problem, we implement a number of quantum and quantum-inspired algorithms on different hardware platforms to solve its discrete formulation using real data from daily prices over 8 years of 52 assets, and do a detailed comparison of the obtained Sharpe ratios, profits and computing times. In particular, we implement classical solvers (Gekko, exhaustive), D-Wave Hybrid quantum annealing, two different approaches based on Variational Quantum Eigensolvers on IBM-Q (one of them brand-new and tailored to the problem), and for the first time in this context also a quantum-inspired optimizer based on Tensor Networks. In order to fit the data into each specific hardware platform, we also consider doing a preprocessing based on clustering of assets. From our comparison, we conclude that D-Wave Hybrid and Tensor Networks are able to handle the largest systems, where we do calculations up to 1272 fully-connected qubits for demonstrative purposes. Finally, we also discuss how to mathematically implement other possible real-life constraints, as well as several ideas to further improve the performance of the studied methods.

preprint2020arXiv

Spin-$\frac{1}{2}$ kagome Heisenberg antiferromagnet with strong breathing anisotropy

We study the zero-temperature phase diagram of the spin-$\frac{1}{2}$ Heisenberg model with breathing anisotropy (i.e., with different coupling strength on the upward and downward triangles) on the kagome lattice. Our study relies on large scale tensor network simulations based on infinite projected entangled-pair state and infinite projected entangled-simplex state methods adapted to the kagome lattice. Our energy analysis suggests that the U(1) algebraic quantum spin-liquid (QSL) ground-state of the isotropic Heisenberg model is stable up to very large breathing anisotropy until it breaks down to a critical lattice-nematic phase that breaks rotational symmetry in real space through a first-order quantum phase transition. Our results also provide further insight into the recent experiment on vanadium oxyfluoride compounds which has been shown to be relevant platforms for realizing QSL in the presence of breathing anisotropy.

preprint2020arXiv

Topological order on the Bloch sphere

A Bloch sphere is the geometrical representation of an arbitrary two-dimensional Hilbert space. Possible classes of entanglement and separability for the pure and mixed states on the Bloch sphere were suggested by [M. Boyer, R. Liss, T. Mor, PRA 95, 032308 (2017)]. Here we construct a Bloch sphere for the Hilbert space spanned by one of the ground states of Kitaev's toric code model and one of its closest product states. We prove that this sphere contains only one separable state, thus belonging to the fourth class suggested by the said paper. We furthermore study the topological order of the pure states on its surface and conclude that, according to conventional definitions, only one state (the toric code ground state) seems to present non-trivial topological order. We conjecture that most of the states on this Bloch sphere are neither ``trivial'' states (namely, they cannot be generated from a product state using a trivial circuit) nor topologically ordered. In addition, we show that the whole setting can be understood in terms of Grover rotations with gauge symmetry, akin to the quantum search algorithm.

preprint2019arXiv

A programming guide for tensor networks with global $SU(2)$ symmetry

This paper is a manual with tips and tricks for programming tensor network algorithms with global $SU(2)$ symmetry. We focus on practical details that are many times overlooked when it comes to implementing the basic building blocks of codes, such as useful data structures to store the tensors, practical ways of manipulating them, and so forth. Here we do not restrict ourselves to any specific tensor network method, but keep always in mind that the implementation should scale well for simulations of higher-dimensional systems using, e.g., Projected Entangled Pair States, where tensors with many indices may show up. To this end, the structural tensors (or intertwiners) that arise in the usual decomposition of $SU(2)$-symmetric tensors are never explicitly stored throughout the simulation. Instead, we store and manipulate the corresponding fusion trees - an algebraic specification of the symmetry constraints on the tensor - in order to implement basic $SU(2)$-symmetric tensor operations.

preprint2019arXiv

Topological $\mathbb{Z}_2$ RVB quantum spin liquid on the ruby lattice

We construct a short-range resonating valence-bond state (RVB) on the ruby lattice, using projected entangled-pair states (PEPS) with bond dimension $D=3$. By introducing non-local moves to the dimer patterns on the torus, we distinguish four distinct sectors in the space of dimer coverings, which is a signature of the topological nature of the RVB wave function. Furthermore, by calculating the reduced density matrix of a bipartition of the RVB state on an infinite cylinder and exploring its entanglement entropy, we confirm the topological nature of the RVB wave function by obtaining non-zero topological contribution, $γ=-\rm{ln}\ 2$, consistent with that of a $\mathbb{Z}_2$ topological quantum spin liquid. We also calculate the ground-state energy of the spin-$\frac{1}{2}$ antiferromagnetic Heisenberg model on the ruby lattice and compare it with the RVB energy. Finally, we construct a quantum-dimer model for the ruby lattice and discuss it as a possible parent Hamiltonian for the RVB wave function.

preprint2016arXiv

Spin-S Kagome quantum antiferromagnets in a field with tensor networks

Spin-$S$ Heisenberg quantum antiferromagnets on the Kagome lattice offer, when placed in a magnetic field, a fantastic playground to observe exotic phases of matter with (magnetic analogs of) superfluid, charge, bond or nematic orders, or a coexistence of several of the latter. In this context, we have obtained the (zero temperature) phase diagrams up to $S=2$ directly in the thermodynamic limit thanks to infinite Projected Entangled Pair States (iPEPS), a tensor network numerical tool. We find incompressible phases characterized by a magnetization plateau vs field and stabilized by spontaneous breaking of point group or lattice translation symmetry(ies). The nature of such phases may be semi-classical, as the plateaus at $\frac{1}{3}$th, $(1-\frac{2}{9S})$th and $(1-\frac{1}{9S})$th of the saturated magnetization (the latter followed by a macroscopic magnetization jump), or fully quantum as the spin-$\frac{1}{2}$ $\frac{1}{9}$-plateau exhibiting coexistence of charge and bond orders. Upon restoration of the spin rotation $U(1)$ symmetry a finite compressibility appears, although lattice symmetry breaking persists. For integer spin values we also identify spin gapped phases at low enough field, such as the $S=2$ (topologically trivial) spin liquid with no symmetry breaking, neither spin nor lattice.

preprint2016arXiv

Symmetry-protected intermediate trivial phases in quantum spin chains

Symmetry-protected trivial (SPt) phases of matter are the product-state analogue of symmetry-protected topological (SPT) phases. This means, SPt phases can be adiabatically connected to a product state by some path that preserves the protecting symmetry. Moreover, SPt and SPT phases can be adiabatically connected to each other when interaction terms that break the symmetries protecting the SPT order are added in the Hamiltonian. It is also known that spin-1 SPT phases in quantum spin chains can emerge as effective intermediate phases of spin-2 Hamiltonians. In this paper we show that a similar scenario is also valid for SPt phases. More precisely, we show that for a given spin-2 quantum chain, effective intermediate spin-1 SPt phases emerge in some regions of the phase diagram, these also being adiabatically connected to non-trivial intermediate SPT phases. We characterize the phase diagram of our model by studying quantities such as the entanglement entropy, symmetry-related order parameters, and 1-site fidelities. Our numerical analysis uses Matrix Product States (MPS) and the infinite Time-Evolving Block Decimation (iTEBD) method to approximate ground states of the system in the thermodynamic limit. Moreover, we provide a field theory description of the possible quantum phase transitions between the SPt phases. Together with the numerical results, such a description shows that the transitions may be described by Conformal Field Theories (CFT) with central charge c=1. Our results are in agreement, and further generalize, those in [Y. Fuji, F. Pollmann, M. Oshikawa, Phys. Rev. Lett. 114, 177204 (2015)].

preprint2016arXiv

Systematic construction of spin liquids on the square lattice from tensor networks with SU(2) symmetry

We elaborate a simple classification scheme of all rank-5 SU(2)-spin rotational symmetric tensors according to i) the on-site physical spin-$S$, (ii) the local Hilbert space $V^{\otimes 4}$ of the four virtual (composite) spins attached to each site and (iii) the irreducible representations of the $C_{4v}$ point group of the square lattice. We apply our scheme to draw a complete list of all SU(2)-symmetric translationally and rotationally-invariant Projected Entangled Pair States (PEPS) with bond dimension $D\leqslant 6$. All known SU(2)-symmetric PEPS on the square lattice are recovered and simple generalizations are provided in some cases. More generally, to each of our symmetry class can be associated a $({\cal D}-1)$-dimensional manifold of spin liquids (potentially) preserving lattice symmetries and defined in terms of ${\cal D}$ independent tensors of a given bond dimension $D$. In addition, generic (low-dimensional) families of PEPS explicitly breaking either (i) particular point-group lattice symmetries (lattice nematics) or (ii) time reversal symmetry (chiral spin liquids) or (iii) SU(2)-spin rotation symmetry down to $U(1)$ (spin nematics or Néel antiferromagnets) can also be constructed. We apply this framework to search for new topological chiral spin liquids characterized by well-defined chiral edge modes, as revealed by their entanglement spectrum. In particular, we show how the symmetrization of a double-layer PEPS leads to a chiral topological state with a gapless edge described by a SU(2)$_2$ Wess-Zumino-Witten model.

preprint2016arXiv

The spin-1/2 Kagome XXZ model in a field: competition between lattice nematic and solid orders

We study numerically the spin-1/2 XXZ model in a field on an infinite Kagome lattice. We use different algorithms based on infinite Projected Entangled Pair States (iPEPS) for this, namely: (i) with simplex tensors and 9-site unit cell, and (ii) coarse-graining three spins in the Kagome lattice and mapping it to a square-lattice model with nearest-neighbor interactions, with usual PEPS tensors, 6- and 12-site unit cells. Similarly to our previous calculation at the SU(2)-symmetric point (Heisenberg Hamiltonian), for any anisotropy from the Ising limit to the XY limit, we also observe the emergence of magnetization plateaus as a function of the magnetic field, at $m_z = \frac{1}{3}$ using 6- 9- and 12-site PEPS unit cells, and at $m_z = \frac{1}{9}, \frac{5}{9}$ and $\frac{7}{9}$ using a 9-site PEPS unit cell, the later set-up being able to accommodate $\sqrt{3} \times \sqrt{3}$ solid order. We also find that, at $m_z = \frac{1}{3}$, (lattice) nematic and $\sqrt{3} \times \sqrt{3}$ VBC-order states are degenerate within the accuracy of the 9-site simplex-method, for all anisotropy. The 6- and 12-site coarse-grained PEPS methods produce almost-degenerate nematic and $1 \times 2$ VBC-Solid orders. Within our accuracy, the 6-site coarse-grained PEPS method gives slightly lower energies, which can be explained by the larger amount of entanglement this approach can handle, even when the PEPS unit-cell is not commensurate with the expected ground state. Furthermore, we do not observe chiral spin liquid behaviors at and close to the XY point, as has been recently proposed. Our results are the first tensor network investigations of the XXZ spin chain in a field, and reveal the subtle competition between nearby magnetic orders in numerical simulations of frustrated quantum antiferromagnets, as well as the delicate interplay between energy optimization and symmetry in tensor networks.

preprint2015arXiv

All spin-1 topological phases in a single spin-2 chain

Here we study the emergence of different Symmetry-Protected Topological (SPT) phases in a spin-2 quantum chain. We consider a Heisenberg-like model with bilinear, biquadratic, bicubic, and biquartic nearest-neighbor interactions, as well as uniaxial anisotropy. We show that this model contains four different effective spin-1 SPT phases, corresponding to different representations of the $(\mathbb{Z}_2 \times \mathbb{Z}_2) + T$ symmetry group, where $\mathbb{Z}_2$ is some $π$-rotation in the spin internal space and $T$ is time-reversal. One of these phases is equivalent to the usual spin-1 Haldane phase, while the other three are different but also typical of spin-1 systems. The model also exhibits an $SO(5)$-Haldane phase. Moreover, we also find that the transitions between the different effective spin-1 SPT phases are continuous, and can be described by a $c=2$ conformal field theory. At such transitions, indirect evidence suggests a possible effective field theory of four massless Majorana fermions. The results are obtained by approximating the ground state of the system in the thermodynamic limit using Matrix Product States via the infinite Time Evolving Block Decimation method, as well as by effective field theory considerations. Our findings show, for the first time, that different large effective spin-1 SPT phases separated by continuous quantum phase transitions can be stabilized in a simple quantum spin chain.

preprint2015arXiv

The iPEPS algorithm, improved: fast full update and gauge fixing

The infinite Projected Entangled Pair States (iPEPS) algorithm [J. Jordan et al, PRL 101, 250602 (2008)] has become a useful tool in the calculation of ground state properties of 2d quantum lattice systems in the thermodynamic limit. Despite its many successful implementations, the method has some limitations in its present formulation which hinder its application to some highly-entangled systems. The purpose of this paper is to unravel some of these issues, in turn enhancing the stability and efficiency of iPEPS methods. For this, we first introduce the fast full update scheme, where effective environment and iPEPS tensors are both simultaneously updated (or evolved) throughout time. As we shall show, this implies two crucial advantages: (i) dramatic computational savings, and (ii) improved overall stability. Besides, we extend the application of the \emph{local gauge fixing}, successfully implemented for finite-size PEPS [M. Lubasch, J. Ignacio Cirac, M.-C. Bañuls, PRB 90, 064425 (2014)], to the iPEPS algorithm. We see that the gauge fixing not only further improves the stability of the method, but also accelerates the convergence of the alternating least squares sweeping in the (either "full" or "fast full") tensor update scheme. The improvement in terms of computational cost and stability of the resulting "improved" iPEPS algorithm is benchmarked by studying the ground state properties of the quantum Heisenberg and transverse-field Ising models on an infinite square lattice.

preprint2014arXiv

A Practical Introduction to Tensor Networks: Matrix Product States and Projected Entangled Pair States

This is a partly non-technical introduction to selected topics on tensor network methods, based on several lectures and introductory seminars given on the subject. It should be a good place for newcomers to get familiarized with some of the key ideas in the field, specially regarding the numerics. After a very general introduction we motivate the concept of tensor network and provide several examples. We then move on to explain some basics about Matrix Product States (MPS) and Projected Entangled Pair States (PEPS). Selected details on some of the associated numerical methods for 1d and 2d quantum lattice systems are also discussed.

preprint2014arXiv

Advances on Tensor Network Theory: Symmetries, Fermions, Entanglement, and Holography

This is a short review on selected theory developments on Tensor Network (TN) states for strongly correlated systems. Specifically, we briefly review the effect of symmetries in TN states, fermionic TNs, the calculation of entanglement Hamiltonians from Projected Entangled Pair States (PEPS), and the relation between the Multi-scale Entanglement Renormalization Ansatz (MERA) and the AdS/CFT or gauge/gravity duality. We stress the role played by entanglement in the emergence of several physical properties and objects through the TN language. Some recent results along these lines are also discussed.

preprint2014arXiv

Topological Minimally Entangled States via Geometric Measure

Here we show how the Minimally Entangled States (MES) of a 2d system with topological order can be identified using the geometric measure of entanglement. We show this by minimizing this measure for the doubled semion, doubled Fibonacci and toric code models on a torus with non-trivial topological partitions. Our calculations are done either quasi-exactly for small system sizes, or using the tensor network approach in [R. Orus, T.-C. Wei, O. Buerschaper, A. Garcia-Saez, arXiv:1406.0585] for large sizes. As a byproduct of our methods, we see that the minimisation of the geometric entanglement can also determine the number of Abelian quasiparticle excitations in a given model. The results in this paper provide a very efficient and accurate way of extracting the full topological information of a 2d quantum lattice model from the multipartite entanglement structure of its ground states.

preprint2014arXiv

Topological transitions from multipartite entanglement with tensor networks: a procedure for sharper and faster characterization

Topological order in a 2d quantum matter can be determined by the topological contribution to the entanglement Rényi entropies. However, when close to a quantum phase transition, its calculation becomes cumbersome. Here we show how topological phase transitions in 2d systems can be much better assessed by multipartite entanglement, as measured by the topological geometric entanglement of blocks. Specifically, we present an efficient tensor network algorithm based on Projected Entangled Pair States to compute this quantity for a torus partitioned into cylinders, and then use this method to find sharp evidence of topological phase transitions in 2d systems with a string-tension perturbation. When compared to tensor network methods for Rényi entropies, our approach produces almost perfect accuracies close to criticality and, on top, is orders of magnitude faster. The method can be adapted to deal with any topological state of the system, including minimally entangled ground states. It also allows to extract the critical exponent of the correlation length, and shows that there is no continuous entanglement-loss along renormalization group flows in topological phases.

preprint2013arXiv

Geometric Entanglement in Topologically Ordered States

Here we investigate the connection between topological order and the geometric entanglement, as measured by the logarithm of the overlap between a given state and its closest product state of blocks. We do this for a variety of topologically-ordered systems such as the toric code, double semion, color code, and quantum double models. As happens for the entanglement entropy, we find that for sufficiently large block sizes the geometric entanglement is, up to possible sub-leading corrections, the sum of two contributions: a bulk contribution obeying a boundary law times the number of blocks, and a contribution quantifying the underlying pattern of long-range entanglement of the topologically-ordered state. This topological contribution is also present in the case of single-spin blocks in most cases, and constitutes an alternative characterisation of topological order for these quantum states based on a multipartite entanglement measure. In particular, we see that the topological term for the 2D color code is twice as much as the one for the toric code, in accordance with recent renormalization group arguments. Motivated by this, we also derive a general formalism to obtain upper- and lower-bounds to the geometric entanglement of states with a non-Abelian symmetry, and which we use to analyse quantum double models. Furthermore, we also analyse the robustness of the topological contribution using renormalization and perturbation theory arguments, as well as a numerical estimation for small systems. Some of our results rely on the ability to disentangle single sites from the system, which is always possible in our framework. Aditionally, we relate our results to the relative entropy of entanglement in topological systems, and discuss some tensor network numerical approaches that could be used to extract the topological contribution for large systems beyond exactly-solvable models.

preprint2012arXiv

Comment on "Topological quantum phase transitions of attractive spinless fermions in a honeycomb lattice" by Poletti D. et al

In a recent letter [D. Poletti et al., EPL 93, 37008 (2011)] a model of attractive spinless fermions on the honeycomb lattice at half filling has been studied by mean-field theory, where distinct homogenous phases at rather large attraction strength $V>3.36$, separated by (topological) phase transitions, have been predicted. In this comment we argue that without additional interactions the ground states in these phases are not stable against phase separation. We determine the onset of phase separation at half filling $V_{ps}\approx 1.7$ by means of infinite projected entangled-pair states (iPEPS) and exact diagonalization.

preprint2012arXiv

Exploring corner transfer matrices and corner tensors for the classical simulation of quantum lattice systems

In this paper we explore the practical use of the corner transfer matrix and its higher-dimensional generalization, the corner tensor, to develop tensor network algorithms for the classical simulation of quantum lattice systems of infinite size. This exploration is done mainly in one and two spatial dimensions (1d and 2d). We describe a number of numerical algorithms based on corner matri- ces and tensors to approximate different ground state properties of these systems. The proposed methods make also use of matrix product operators and projected entangled pair operators, and naturally preserve spatial symmetries of the system such as translation invariance. In order to assess the validity of our algorithms, we provide preliminary benchmarking calculations for the spin-1/2 quantum Ising model in a transverse field in both 1d and 2d. Our methods are a plausible alternative to other well-established tensor network approaches such as iDMRG and iTEBD in 1d, and iPEPS and TERG in 2d. The computational complexity of the proposed algorithms is also considered and, in 2d, important differences are found depending on the chosen simulation scheme. We also discuss further possibilities, such as 3d quantum lattice systems, periodic boundary conditions, and real time evolution. This discussion leads us to reinterpret the standard iTEBD and iPEPS algorithms in terms of corner transfer matrices and corner tensors. Our paper also offers a perspective on many properties of the corner transfer matrix and its higher-dimensional generalizations in the light of novel tensor network methods.

preprint2012arXiv

Topological Geometric Entanglement

Here we show the connection between topological order and the geometric entanglement, as measured by the logarithm of the overlap between a given state and its closest product state of blocks, for the topological universality class of the toric code model. As happens for the entanglement entropy, we find that for large block sizes the geometric entanglement is, up to possible subleading corrections, the sum of two contributions: a non-universal bulk contribution obeying a boundary law times the number of blocks, and a universal contribution quantifying the underlying pattern of long-range entanglement of a topologically-ordered state.

preprint2011arXiv

Effective field theory for the SO(n) bilinear-biquadratic spin chain

We present a low-energy effective field theory to describe the SO(n) bilinear-biquadratic spin chain. We start with n=6 and construct the effective theory by using six Majorana fermions. After determining various correlation functions we characterize the phases and establish the relation between the effective theories for SO(6) and SO(5). Together with the known results for n=3 and 4, a reduction mechanism is proposed to understand the ground state for arbitrary SO(n). Also, we provide a generalization of the Lieb-Schultz-Mattis theorem for SO(n). The implications of our results for entanglement and correlation functions are discussed.

preprint2011arXiv

Entanglement and SU(n) symmetry in one-dimensional valence bond solid states

Here we evaluate the many-body entanglement properties of a generalized SU(n) valence bond solid state on a chain. Our results follow from a derivation of the transfer matrix of the system which, in combination with symmetry properties, allows for a new, elegant and straightforward evaluation of different entanglement measures. In particular, the geometric entanglement per block, correlation length, von Neumann and Rényi entropies of a block, localizable entanglement and entanglement length are obtained in a very simple way. All our results are in agreement with previous derivations for the SU(2) case.

preprint2011arXiv

Geometric entanglement of one-dimensional systems: bounds and scalings in the thermodynamic limit

In this paper the geometric entanglement (GE) of systems in one spatial dimension (1D) and in the thermodynamic limit is analyzed focusing on two aspects. First, we reexamine the calculation of the GE for translation-invariant matrix product states (MPSs) in the limit of infinite system size. We obtain a lower bound to the GE which collapses to an equality under certain sufficient conditions that are fulfilled by many physical systems, such as those having unbroken space (P) or space-time (PT) inversion symmetry. Our analysis justifies the validity of several derivations carried out in previous works. Second, we derive scaling laws for the GE per site of infinite-size 1D systems with correlation length $ξ\gg 1$. In the case of MPSs, we combine this with the theory of finite-entanglement scaling, allowing to understand the scaling of the GE per site with the MPS bond dimension at conformally invariant quantum critical points.

preprint2011arXiv

Intermediate Haldane phase in spin-2 quantum chains with uniaxial anisotropy

We provide evidence of an intermediate Haldane phase in a spin-2 quantum chain. By combining effective field theory and numerical approaches, we show that the phase diagram of the proposed model includes SO(5) Haldane, intermediate Haldane, and large-D phases. We determine the characteristic properties of these phases, including edge states, string order parameters, and degeneracies in the entanglement spectrum. The symmetries responsible for the degeneracy patterns observed in the entanglement spectrum are also discussed.

preprint2011arXiv

Phase diagram of the SO(n) bilinear-biquadratic chain from many-body entanglement

Here we investigate the phase diagram of the SO(n) bilinear-biquadratic quantum spin chain by studying the global quantum correlations of the ground state. We consider the cases of n=3,4 and 5 and focus on the geometric entanglement in the thermodynamic limit. Apart from capturing all the known phase transitions, our analysis shows a number of novel distinctive behaviors in the phase diagrams which we conjecture to be general and valid for arbitrary n. In particular, we provide an intuitive argument in favor of an infinite entanglement length in the system at a purely-biquadratic point. Our results are also compared to other methods, such as fidelity diagrams.

preprint2010arXiv

Finite-Size Geometric Entanglement from Tensor Network Algorithms

The global geometric entanglement is studied in the context of newly-developed tensor network algorithms for finite systems. For one-dimensional quantum spin systems it is found that, at criticality, the leading finite-size correction to the global geometric entanglement per site behaves as $b/n$, where $n$ is the size of the system and $b$ a given coefficient. Our conclusion is based on the computation of the geometric entanglement per spin for the quantum Ising model in a transverse magnetic field and for the spin-1/2 XXZ model. We also discuss the possibility of coefficient $b$ being universal.

preprint2010arXiv

Simulation of strongly correlated fermions in two spatial dimensions with fermionic Projected Entangled-Pair States

We explain how to implement, in the context of projected entangled-pair states (PEPS), the general procedure of fermionization of a tensor network introduced in [P. Corboz, G. Vidal, Phys. Rev. B 80, 165129 (2009)]. The resulting fermionic PEPS, similar to previous proposals, can be used to study the ground state of interacting fermions on a two-dimensional lattice. As in the bosonic case, the cost of simulations depends on the amount of entanglement in the ground state and not directly on the strength of interactions. The present formulation of fermionic PEPS leads to a straightforward numerical implementation that allowed us to recycle much of the code for bosonic PEPS. We demonstrate that fermionic PEPS are a useful variational ansatz for interacting fermion systems by computing approximations to the ground state of several models on an infinite lattice. For a model of interacting spinless fermions, ground state energies lower than Hartree-Fock results are obtained, shifting the boundary between the metal and charge-density wave phases. For the t-J model, energies comparable with those of a specialized Gutzwiller-projected ansatz are also obtained.

preprint2010arXiv

Visualizing elusive phase transitions with geometric entanglement

We show that by examining the global geometric entanglement it is possible to identify "elusive" or hard to detect quantum phase transitions. We analyze several one-dimensional quantum spin chains and demonstrate the existence of non-analyticities in the geometric entanglement, in particular across a Kosterlitz-Thouless transition and across a transition for a gapped deformed Affleck-Kennedy-Lieb-Tasaki chain. The observed non-analyticities can be understood and classified in connection to the nature of the transitions, and are in sharp contrast to the analytic behavior of all the two-body reduced density operators and their derived entanglement measures.

preprint2009arXiv

Numerical study of the hard-core Bose-Hubbard Model on an Infinite Square Lattice

We present a study of the hard-core Bose-Hubbard model at zero temperature on an infinite square lattice using the infinite Projected Entangled Pair State algorithm [Jordan et al., Phys. Rev. Lett. 101, 250602 (2008)]. Throughout the whole phase diagram our values for the ground state energy, particle density and condensate fraction accurately reproduce those previously obtained by other methods. We also explore ground state entanglement, compute two-point correlators and conduct a fidelity-based analysis of the phase diagram. Furthermore, for illustrative purposes we simulate the response of the system when a perturbation is suddenly added to the Hamiltonian.

preprint2009arXiv

Simulation of two dimensional quantum systems on an infinite lattice revisited: corner transfer matrix for tensor contraction

An extension of the projected entangled-pair states (PEPS) algorithm to infinite systems, known as the iPEPS algorithm, was recently proposed to compute the ground state of quantum systems on an infinite two-dimensional lattice. Here we investigate a modification of the iPEPS algorithm, where the environment is computed using the corner transfer matrix renormalization group (CTMRG) method, instead of using one-dimensional transfer matrix methods as in the original proposal. We describe a variant of the CTMRG that addresses different directions of the lattice independently, and use it combined with imaginary time evolution to compute the ground state of the two-dimensional quantum Ising model. Near criticality, the modified iPEPS algorithm is seen to provide a better estimation of the order parameter and correlators.

preprint2005arXiv

Entanglement and majorization in (1+1)-dimensional quantum systems

Motivated by the idea of entanglement loss along Renormalization Group flows, analytical majorization relations are proven for the ground state of (1+1)-dimensional conformal field theories. For any of these theories, majorization is proven to hold in the spectrum of the reduced density matrices in a bipartite system when changing the size L of one of the subsystems. Continuous majorization along uniparametric flows is also proven as long as part of the conformal structure is preserved under the deformation and some monotonicity conditions hold as well. As particular examples of our derivations, we study the cases of the XX, Heisenberg and XY quantum spin chains. Our results provide in a rigorous way explicit proves for all the majorization conjectures raised by Latorre, Lutken, Rico, Vidal and Kitaev in previous papers on quantum spin chains.