Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Entanglement in N-harmonium: bosons and fermions

C. L. Benavides-riveros, I. V. Toranzo, J. S. Dehesa

Published Jan 1, 2014• 25 citations
OpenAlex

Abstract

The entanglement of the N-harmonium system (i.e., a completely-integrable model of N particles where both the confinement and the two-particle interaction are har-monic) is shown to be analytically determined in terms of N and the relative interaction strength. First, the entanglement of the N-boson system via the von Neumann entropy is computed by use of the corresponding natural occupation numbers. Below a critical value Nc of N it is found that, for positive values of the coupling constant the entan-glement grows when the number of particles is increasing; opposite situation occurs for N> Nc. Then the entanglement of the N-fermion system is investigated by means of the linear entropy, which is a linearized form of the von Neumann entropy of the one-body reduced density matrix. We have shown that the entanglement grows when the number of particles is increasing in the region of negative coupling. For positively increasing values of the coupling the entanglement as a function of N varies according to various regimes. Moreover, in the strong coupling regime the entropy of the system decreases when N is growing. For fixed N, the entanglement is found (a) to decrease (increase) for negatively (positively) increasing values of the coupling constant, and (b) to grow when the energy is increasing. Finally, let us point out that the spatial and spin contributions to the total entanglement are of comparable size.

Related papers

Entanglement in N-harmonium: bosons and fermions

C L Benavides-Riveros, I V Toranzo, J S Dehesa

The ground-state entanglement of a single particle of the N-harmonium system (i.e., a completely-integrable model of $N$ particles where both the confinement and the two-particle interaction are harmonic) is shown to be analytically determined in terms of $N$ and the relative interaction strength. For bosons, we compute the von Neumann entropy of the one-body reduced density matrix by using the corresponding natural occupation numbers. There exists a critical number $N_c$ of particles so that below it, for positive values of the coupling constant, the entanglement grows when the number of particles is increasing; the opposite occurs for $N > N_c$. For fermions, we compute the one-body reduced density matrix for the closed-shell spinned case. In the strong coupling regime, the linear entropy of the system decreases when $N$ is growing. For fixed $N$, the entanglement is found (a) to decrease (increase) for negatively (positively) increasing values of the coupling constant, and (b) to grow when the energy is increasing. Moreover, the spatial and spin contributions to the total entanglement are found to be of comparable size.

arXiv:1404.4447•Journal of Physics B Atomic Molecular and Optical Physics•Sep 24, 2014•35 citations

Entanglement and bifurcation in the integrable dimer

Xi-Wen Hou, Jinghua Chen, Bambi Hu

In this Brief Report the properties of both dynamical and static entanglement in the integrable quantum dimer are studied in terms of the reduced-density linear entropy and von Neumann entropy with various coupling parameters, total boson numbers, and initial states. The mean entanglement, which is defined to be averaged over time, is used to describe the influence of the classical separatrix on the behavior of entanglement. It is shown that the mean entanglement exhibits a maximum near the position of the corresponding classical separatrix energy and that the static entanglement of the state with the largest eigenvalue of the quantum spectrum displays a maximum near the bifurcation point. For weak coupling and larger total boson number the maximum entanglement state is exactly at the position of the classical separatrix and bifurcation. In strong coupling all initial states have nearly the same mean entanglement.

Physical Review A•Mar 10, 2005•37 citations

Natural orbitals and occupation numbers for harmonium: Fermions versus bosons

Christian Schilling

For a quantum system of $N$ identical, harmonically interacting particles in a one-dimensional harmonic trap we calculate for the bosonic and fermionic ground state the corresponding one-particle reduced density operator ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\rho}}}_{1}$ analytically. In case of bosons ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\rho}}}_{1}$ is a Gibbs state for an effective harmonic oscillator. Hence the natural orbitals are Hermite functions and their occupation numbers obey a Boltzmann distribution. Intriguingly, for fermions with not too large couplings the natural orbitals coincide up to just a small error with the bosonic ones. In case of strong coupling this still holds qualitatively. Moreover, the decay of the decreasingly ordered fermionic natural occupation numbers is given by the bosonic one, but modified by an algebraic prefactor. Significant differences to bosons occur only for the largest occupation numbers. After all, the ``gap'' at the ``Fermi level'' decreases with increasing coupling strength but remains well pronounced even for strong interaction.

arXiv:1307.6858•Physical Review A•Oct 9, 2013•38 citations

Entanglement entropy in quantum many particle systems and their simulation via ansatz states

Thomas Barthel, Ulrich Schollwöck

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.

•Jan 1, 2009•0 citations

Entanglement of an impurity in a few-body one-dimensional ideal Bose system

Miguel Ángel García-March, A S Dehkharghani, Nikolaj Thomas Zinner

Abstract We study the correlation between an impurity and a small ensemble of bosonic particles in one dimension. Our study analyzes the one-body density matrix and calculates the corresponding von Neumann entanglement entropy as a function of the interaction strength between the impurity and the bosons when all particles have the same mass. We show that the entropy grows very fast for small and moderate interaction strength and then increases slowly toward the strongly interacting regime. Then we study the effect over the quantum correlations of a mass imbalance between the impurity and the bosons. In the strongly interacting case, we discover that when the impurity is much heavier than the bosons, then we have the least possible correlation. However, the entropy tops its maximum when the mass ratio is between 3 and 4 in the case where there are four bosonic particles and then falls off to its minimum for higher mass imbalance.

arXiv:1512.03367•Journal of Physics B Atomic Molecular and Optical Physics•Mar 17, 2016•18 citations

Natural Multiparticle Entanglement in a Fermi Gas

Christian Lunkes, Časlav Brukner, Vlatko Vedral

We investigate multipartite entanglement in a noninteracting fermion gas, as a function of fermion separation, starting from the many particle fermion density matrix. We prove that all multiparticle entanglement can be built only out of two-fermion entanglement. Although from the Pauli exclusion principle we would always expect entanglement to decrease with fermion distance, we surprisingly find the opposite effect for certain fermion configurations. The von Neumann entropy is found to be proportional to the volume for a large number of particles even when they are arbitrarily close to each other. We will illustrate our results using different configurations of two, three, and four fermions at zero temperature although all our results can be applied to any temperature and any number of particles.

arXiv:quant-ph/0502122•Physical Review Letters•Jul 13, 2005•41 citations
The entanglement of the N-harmonium system (i.e., a completely-integrable model of N particles where both the confinement and the two-particle interaction are har-monic) is shown to be analytically determined in terms of N and the relative interaction strength. First, the entanglement of the N-boson system via the von Neumann entropy is computed by use of the corresponding natural occupation numbers. Below a critical value Nc of N it is found that, for positive values of the coupling constant the entan-glement grows when the number of particles is increasing; opposite situation occurs for N> Nc. Then the entanglement of the N-fermion system is investigated by means of the linear entropy, which is a linearized form of the von Neumann entropy of the one-body reduced density matrix. We have shown that the entanglement grows when the number of particles is increasing in the region of negative coupling. For positively increasing values of the coupling the entanglement as a function of N varies according to various regimes. Moreover, in the strong coupling regime the entropy of the system decreases when N is growing. For fixed N, the entanglement is found (a) to decrease (increase) for negatively (positively) increasing values of the coupling constant, and (b) to grow when the energy is increasing. Finally, let us point out that the spatial and spin contributions to the total entanglement are of comparable size.
The ground-state entanglement of a single particle of the N-harmonium system (i.e., a completely-integrable model of NNN particles where both the confinement and the two-particle interaction are harmonic) is shown to be analytically determined in terms of NNN and the relative interaction strength. For bosons, we compute the von Neumann entropy of the one-body reduced density matrix by using the corresponding natural occupation numbers. There exists a critical number NcN_cNc​ of particles so that below it, for positive values of the coupling constant, the entanglement grows when the number of particles is increasing; the opposite occurs for N>NcN > N_cN>Nc​. For fermions, we compute the one-body reduced density matrix for the closed-shell spinned case. In the strong coupling regime, the linear entropy of the system decreases when NNN is growing. For fixed NNN, the entanglement is found (a) to decrease (increase) for negatively (positively) increasing values of the coupling constant, and (b) to grow when the energy is increasing. Moreover, the spatial and spin contributions to the total entanglement are found to be of comparable size.
In this Brief Report the properties of both dynamical and static entanglement in the integrable quantum dimer are studied in terms of the reduced-density linear entropy and von Neumann entropy with various coupling parameters, total boson numbers, and initial states. The mean entanglement, which is defined to be averaged over time, is used to describe the influence of the classical separatrix on the behavior of entanglement. It is shown that the mean entanglement exhibits a maximum near the position of the corresponding classical separatrix energy and that the static entanglement of the state with the largest eigenvalue of the quantum spectrum displays a maximum near the bifurcation point. For weak coupling and larger total boson number the maximum entanglement state is exactly at the position of the classical separatrix and bifurcation. In strong coupling all initial states have nearly the same mean entanglement.
For a quantum system of NNN identical, harmonically interacting particles in a one-dimensional harmonic trap we calculate for the bosonic and fermionic ground state the corresponding one-particle reduced density operator ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\rho}}}_{1}$ analytically. In case of bosons ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\rho}}}_{1}$ is a Gibbs state for an effective harmonic oscillator. Hence the natural orbitals are Hermite functions and their occupation numbers obey a Boltzmann distribution. Intriguingly, for fermions with not too large couplings the natural orbitals coincide up to just a small error with the bosonic ones. In case of strong coupling this still holds qualitatively. Moreover, the decay of the decreasingly ordered fermionic natural occupation numbers is given by the bosonic one, but modified by an algebraic prefactor. Significant differences to bosons occur only for the largest occupation numbers. After all, the ``gap'' at the ``Fermi level'' decreases with increasing coupling strength but remains well pronounced even for strong interaction.
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.
Abstract We study the correlation between an impurity and a small ensemble of bosonic particles in one dimension. Our study analyzes the one-body density matrix and calculates the corresponding von Neumann entanglement entropy as a function of the interaction strength between the impurity and the bosons when all particles have the same mass. We show that the entropy grows very fast for small and moderate interaction strength and then increases slowly toward the strongly interacting regime. Then we study the effect over the quantum correlations of a mass imbalance between the impurity and the bosons. In the strongly interacting case, we discover that when the impurity is much heavier than the bosons, then we have the least possible correlation. However, the entropy tops its maximum when the mass ratio is between 3 and 4 in the case where there are four bosonic particles and then falls off to its minimum for higher mass imbalance.
We investigate multipartite entanglement in a noninteracting fermion gas, as a function of fermion separation, starting from the many particle fermion density matrix. We prove that all multiparticle entanglement can be built only out of two-fermion entanglement. Although from the Pauli exclusion principle we would always expect entanglement to decrease with fermion distance, we surprisingly find the opposite effect for certain fermion configurations. The von Neumann entropy is found to be proportional to the volume for a large number of particles even when they are arbitrarily close to each other. We will illustrate our results using different configurations of two, three, and four fermions at zero temperature although all our results can be applied to any temperature and any number of particles.