In this work, we explore physical systems which support not only multipartite interparticle entanglement, but also intraparticle entanglement between different degrees of freedom of the constituent particles and entanglement between different degrees of freedom of different particles, i.e., algebraic entanglement. We derive a simple method for calculating the algebraic entanglement entropy between two of the particles' degrees of freedom from collective states of the whole ensemble. Our procedure makes use of underlying symmetries in these systems, in particular permutation symmetry of the particle indices, and shows a connection between the algebraic entanglement entropy in these systems and the irreducible representations of Lie groups which describe the particles' degrees of freedom. Namely, we use the direct sum over irreducible representations to diagonalize the reduced density matrices in a block-by-block manner, then utilize the multiplicity of these irreducible representations to reproduce the results from an exponentially-scaled Hilbert space in only polynomial complexity. We use this to explore a variety of systems where the constituent particles support two degrees of freedom each with two levels, such as atoms with two electronic states and two momentum states. Notably, these systems may be exactly simulated in a polynomial-scaled Hilbert space, yet they support an algebraic entanglement entropy that grows linearly with the particle number which typically requires an exponentially-scaled Hilbert space.
In this work, we explore physical systems which support not only multipartite interparticle entanglement, but also intraparticle entanglement between different degrees of freedom of the constituent particles and entanglement between different degrees of freedom of different particles, i.e., algebraic entanglement. We derive a simple method for calculating the algebraic entanglement entropy between two of the particles' degrees of freedom from collective states of the whole ensemble. Our procedure makes use of underlying symmetries in these systems, in particular permutation symmetry of the particle indices, and shows a connection between the algebraic entanglement entropy in these systems and the irreducible representations of Lie groups which describe the particles' degrees of freedom. Namely, we use the direct sum over irreducible representations to diagonalize the reduced density matrices in a block-by-block manner, then utilize the multiplicity of these irreducible representations to reproduce the results from an exponentially-scaled Hilbert space in only polynomial complexity. We use this to explore a variety of systems where the constituent particles support two degrees of freedom each with two levels, such as atoms with two electronic states and two momentum states. Notably, these systems may be exactly simulated in a polynomial-scaled Hilbert space, yet they support an algebraic entanglement entropy that grows linearly with the particle number which typically requires an exponentially-scaled Hilbert space.
In contrast to classical many-particle systems for which the number of degrees of freedom grows linearly with the system size, it grows exponentially for quantum many-particle systems. On the one hand, the resulting complexity of the quantum systems has a high potential for technological advances if it can be exploited directly. For example, in recent years, the fundamentals for the realization of quantum computers and their possible applications have been investigated intensely. On the other hand, the high complexity, especially in strongly correlated regimes, often makes the analytical or numerical analysis of underlying physical phenomena very difficult. This is why some systems of high technological relevance, like the high-temperature superconductors, are still unsatisfactorily understood. A main topic of this thesis is the development of efficient numerical methods for the simulation of strongly correlated quantum lattice models. For one-dimensional systems, the density-matrix renormalization-group (DMRG) is such a very successful method. The physical states of interest are approximated within a certain class of ansatz states. These ansatz states are designed in a way that the number of degrees of freedom are prevented from growing exponentially. With the DMRG method, ground states of interacting one-dimensional systems can often be determined with very high precision. This observation indicates already that ground states actually do not exhaust the available exponential number of degrees of freedom but can be described efficiently with much fewer parameters. The ansatz states used in DMRG are the so-called matrix product states. A very useful criterion for the construction of corresponding appropriate ansatz states for higher-dimensional systems is the quantum nonlocality as for example quantified by the entanglement entropy. The ansatz states should be able to reproduce the scaling properties of the quantum nonlocality in the physical states of interest with as few parameters as possible. The first part of the thesis, therefore, provides analytical and numerical analysis of the scaling of quantum nonlocality with the system size or time in different, physically relevant scenarios. For example, the scaling of Renyi entropies and their dependence on boundary conditions is derived within the 1+1-dimensional conformal field theory. Conjectures and analytical indications concerning the properties of entanglement entropy in critical fermionic and bosonic systems are confirmed numerically with high precision. For integrable models in the thermodynamic limit, general preconditions are derived under which subsystems converge to steady states. These steady states are non-thermal and retain information about the initial state. It is shown that the entanglement entropy in such steady states is extensive. For short times, the entanglement entropy grows typically linearly with time, causing an exponential increase in computation costs for the DMRG method. Consequently, time-dependent DMRG simulations usually have no direct access to the long time regime. The second part of the thesis focuses on the development and improvement of the above-mentioned numerical techniques. The time-dependent DMRG is complemented with an extrapolation technique for the evaluated observables. In this way, the problem of the entropy increase can be circumvented, allowing for a precise determination of spectral functions. The method is demonstrated using the example of the Heisenberg antiferromagnet and results are compared to Bethe-Ansatz data for T=0 and quantum Monte Carlo data for T>0. For the simulation of higher-dimensional systems, projected-entangled pair-states (PEPS) and the multiscale entanglement renormalization ansatz (MERA) turn out to be appropriate wavefunction classes. Initially, they could only be used for spin systems. Within this thesis, variants of those state classes are being suggested for the fermionic case. It is shown how the known algorithms for spin systems can be translated into corresponding algorithms for the fermionic systems. The occurring computation costs increase only by a marginal overhead. First numerical benchmarking tests proceeded successfully. Thus, new techniques are available which can for example be used to examine the two-dimensional Hubbard model. This model is a candidate for the description of high-temperature superconductivity.
Quantum Field Theory (QFT) represents a vast generalization of Quantum Mechanics (QM), as it deals with systems that have an infinite number of degrees of freedom. The Stone-von Neumann theorem, which establishes the equivalence of irreducible representations of the canonical commutation relations (CCR) in QM, does not extend to QFT. Consequently, QFT admits multiple inequivalent irreducible representations, leading to a much richer algebraic structure. This essay aims to explore the physics of QFT from the operator algebra perspective, particularly focusing on entanglement entropy. We discuss the role of von Neumann algebras of different types in QFT, describe the local operator algebra approach to QFT, and explain how entanglement entropy can be defined in terms of the algebra of observables. Additionally, we explore the benefits of this approach in concrete applications, specifically in quantum field theory on curved spacetime.
Analyzing the properties of entanglement in many-particle spin-1/2 systems is generally difficult because the system's Hilbert space grows exponentially with the number of constituent particles, N. Fortunately, it is still possible to investigate a many-particle entanglement when the state of the system possesses sufficient symmetry. In this paper, we present a practical method for efficiently computing various bipartite entanglement measures for states in the symmetric subspace and perform these calculations for $N\ensuremath{\sim}{10}^{3}.$ By considering all possible bipartite splits, we construct a picture of the multiscale entanglement in large symmetric systems. In particular, we characterize dynamically generated spin-squeezed states by comparing them to known reference states (e.g., Greenberger-Horne-Zeilinger and Dicke states), and families of states with near-maximal bipartite entropy. We quantify the trade-off between the degree of entanglement and its robustness to particle loss, emphasizing that substantial entanglement need not be fragile.
The density matrix formalism is a fundamental tool in studying various problems in quantum information processing. In the space of density matrices, the most well-known and physically relevant measures are the Hilbert-Schmidt ensemble and the Bures-Hall ensemble. In this work, we propose a generalized ensemble of density matrices, termed quantum interpolating ensemble, which is able to interpolate between these two seemingly unrelated ensembles. As a first step to understand the proposed ensemble, we derive the exact mean formulas of entanglement entropies over such an ensemble generalizing several recent results in the literature. We also derive some key properties of the corresponding orthogonal polynomials relevant to obtaining other statistical information of the entropies. Numerical results demonstrate the usefulness of the proposed ensemble in estimating the degree of entanglement of quantum states.
Although the quantum state of a many-body system (e.g, a system consisting of several atoms) is described in terms of its degrees of freedom, i.e., a large number of variables representing a point on the compact space of quantum states, there is usually a small number of parameters that can be actually measured experimentally. We show that the problem of deciding whether the underlying quantum state is entangled is equivalent to the problem of deciding whether the set of measured parameters are the moments of some measure defined over the compact space of states of the many-body system. By considering quantum entanglement of n identical atoms we arrive at the truncated moment problem defined for symmetric measures over a product of n copies of unit balls in R^d. We consider moments up to degree 2 since these are most readily available experimentally. We derive necessary and sufficient conditions for belonging to the moment cone, which can be expressed via a linear matrix inequality (LMI) of size at most 2d+2, which is independent of n. We can convert the LMIs into a set of explicit semialgebraic inequalities giving necessary and sufficient conditions for membership in the moment cone, and show that the two conditions approach each other in the limit of large n. The inequalities are derived via considering the dual cone of nonnegative polynomials, and its sum-of-squares (SOS) relaxation. We show that the SOS relaxation of the dual cone is asymptotically exact, and using symmetry reduction techniques (J. Pure Appl. Algebra 192, no. 1-3, 95, arXiv:1205.3102), it can be written as a small linear matrix inequality of size at most 2d+2, which is independent of n. For the cone of nonnegative polynomials we also prove an analogue of the half degree principle for globally nonnegative symmetric polynomials (J. Math. Anal. Appl. 284, no. 1, 174, J. Pure Appl. Algebra 216, no. 4, 850).
Area laws describe how entanglement entropy scales and thus provide important necessary conditions for efficient quantum many-body simulation, but they do not, by themselves, yield a direct measure of computational complexity. Here we introduce a complementary framework based on p-particle positivity conditions from reduced density matrix (RDM) theory. These conditions form a hierarchy of N-representability constraints for an RDM to correspond to a valid N-particle quantum system, becoming exact when the Hamiltonian can be expressed as a convex combination of positive semidefinite p-particle operators. We prove a general complexity bound: if a quantum system is solvable with level-p positivity independent of its size, then its entanglement complexity scales polynomially with order p. This theorem connects structural constraints on RDMs with computational tractability and provides a rigorous framework for certifying when many-body methods including RDM methods can efficiently simulate correlated quantum matter and materials.