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

arXiv:1307.6858 · Physical Review A

Natural orbitals and occupation numbers for harmonium: Fermions versus bosons

Christian Schilling

Published Oct 9, 2013• 38 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physreva.88.042105

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 in N-harmonium: bosons and fermions

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

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.

•Jan 1, 2014•25 citations

Contactium: A strongly correlated model system

Jerzy Ciosłowski, Berthold‐Georg Englert, Martin-Isbjörn Trappe +1 more

At the limit of an infinite confinement strength ω, the ground state of a system that comprises two fermions or bosons in harmonic confinement interacting through the Fermi-Huang pseudopotential remains strongly correlated. A detailed analysis of the one-particle description of this "contactium" reveals several peculiarities that are not encountered in conventional model systems (such as the two-electron harmonium atom, ballium, and spherium) involving Coulombic interparticle interactions. First of all, none of the natural orbitals (NOs) {ψn(ω;r)} of the contactium is unoccupied, which implies nonzero collective occupancies for all the angular momenta. Second, the NOs and their non-ascendingly ordered occupation numbers {νn} turn out to be related to the eigenfunctions and eigenvalues of a zero-energy Schrödinger equation with an attractive Gaussian potential. This observation enables the derivation of their properties, such as the n-4/3 asymptotic decay of νn at the n→∞ limit (which differs from that of n-8/3 in the Coulombic systems), the independence of the confinement energy vn=⟨ψn(ω;r)|12ω2r2|ψn(ω;r)⟩ of n, and the n-2/3 asymptotic decay of the respective contribution νntn to the kinetic energy. Upon suitable scaling, the weakly occupied NOs of the contactium turn out to be virtually identical to those of the two-electron harmonium atom at the ω → ∞ limit, despite the entirely different interparticle interactions in these systems.

arXiv:2303.14982•The Journal of Chemical Physics•May 12, 2023•3 citations

Unitary Fermi gas in a harmonic trap

S. Y. Chang, G. F. Bertsch

We present an ab initio calculation of small numbers of trapped, strongly interacting fermions using the Green's function Monte Carlo method. The ground-state energy, density profile, and pairing gap are calculated for particle numbers $N=2--22$ using the parameter-free ``unitary'' interaction. Trial wave functions are taken in the form of correlated pairs in a harmonic oscillator basis. We find that the lowest energies are obtained with a minimum explicit pair correlation beyond that needed to exploit the degeneracy of oscillator states. We find that the energies can be well fitted by the expression ${a}_{\mathit{TF}}{E}_{\mathit{TF}}+\ensuremath{\Delta}\phantom{\rule{0.2em}{0ex}}\mathrm{mod}(N,2)$ where ${E}_{\mathit{TF}}$ is the Thomas-Fermi energy of a noninteracting gas in the trap and $\ensuremath{\Delta}$ is the pairing gap. There is no evidence of a shell correction energy in the systematics, but the density distributions show pronounced shell effects. We find the value $\ensuremath{\Delta}=0.7\ifmmode\pm\else\textpm\fi{}0.2\ensuremath{\omega}$ for the pairing gap. This is smaller than the value found for the uniform gas at a density corresponding to the central density of the trapped gas.

arXiv:physics/0703190•Physical Review A•Aug 2, 2007•74 citations

Ground-state properties of hard-core anyons in a harmonic potential

Yajiang Hao

Using an anyon-fermion mapping method, we investigate the ground-state properties of hard-core anyons confined in a one-dimensional harmonic trap. The concise analytical formula of the reduced one-body density matrix is obtained. Based on this formula, we evaluated the momentum distribution, the natural orbitals, and their occupation distributions for different statistical parameters. The occupation and occupation fraction of the lowest natural orbital versus anyon number are also displayed. It is shown that the ground-state properties of anyons interplay between bosons and fermions continuously. We can expect that the hard-core anyons of larger statistical parameters exhibit similar properties to the hard-core bosons although anyon systems satisfy specific fractional statistics.

arXiv:1604.00458•Physical Review A•Jun 21, 2016•21 citations

Natural orbitals and their occupation numbers for free anyons in the magnetic gauge

Jerzy Ciosłowski, Oliver M. Brown, Tomasz Maciążek

We investigate the properties of natural orbitals and their occupation numbers of the ground state of two noninteracting anyons characterized by the fractional statistics parameter α and confined in a harmonic trap. We work in the boson magnetic gauge where the anyons are modeled as composite bosons with magnetic flux quanta attached to their positions. We derive an asymptotic form of the weakly occupied natural orbitals, and show that their corresponding (ordered descendingly) occupation numbers decay according to the power law n−(4+2α) , where n is the index of the natural orbital. We find remarkable numerical agreement of the theory with the natural orbitals and their occupation numbers computed from the spectral decomposition of the system's wave function. We explain that the same results apply to the fermion magnetic gauge. Published by the American Physical Society 2024

arXiv:2309.14023•Physical Review A•Feb 9, 2024•1 citation
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 \ensuremathρ\ifmmode^\elseˆ\fi1{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\rho}}}_{1}\ensuremathρ\ifmmode^\elseˆ\fi​1​ analytically. In case of bosons \ensuremathρ\ifmmode^\elseˆ\fi1{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\rho}}}_{1}\ensuremathρ\ifmmode^\elseˆ\fi​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.
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.
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.
At the limit of an infinite confinement strength ω, the ground state of a system that comprises two fermions or bosons in harmonic confinement interacting through the Fermi-Huang pseudopotential remains strongly correlated. A detailed analysis of the one-particle description of this "contactium" reveals several peculiarities that are not encountered in conventional model systems (such as the two-electron harmonium atom, ballium, and spherium) involving Coulombic interparticle interactions. First of all, none of the natural orbitals (NOs) {ψn(ω;r)} of the contactium is unoccupied, which implies nonzero collective occupancies for all the angular momenta. Second, the NOs and their non-ascendingly ordered occupation numbers {νn} turn out to be related to the eigenfunctions and eigenvalues of a zero-energy Schrödinger equation with an attractive Gaussian potential. This observation enables the derivation of their properties, such as the n-4/3 asymptotic decay of νn at the n→∞ limit (which differs from that of n-8/3 in the Coulombic systems), the independence of the confinement energy vn=⟨ψn(ω;r)|12ω2r2|ψn(ω;r)⟩ of n, and the n-2/3 asymptotic decay of the respective contribution νntn to the kinetic energy. Upon suitable scaling, the weakly occupied NOs of the contactium turn out to be virtually identical to those of the two-electron harmonium atom at the ω → ∞ limit, despite the entirely different interparticle interactions in these systems.
We present an ab initio calculation of small numbers of trapped, strongly interacting fermions using the Green's function Monte Carlo method. The ground-state energy, density profile, and pairing gap are calculated for particle numbers N=2−−22N=2--22N=2−−22 using the parameter-free ``unitary'' interaction. Trial wave functions are taken in the form of correlated pairs in a harmonic oscillator basis. We find that the lowest energies are obtained with a minimum explicit pair correlation beyond that needed to exploit the degeneracy of oscillator states. We find that the energies can be well fitted by the expression ${a}_{\mathit{TF}}{E}_{\mathit{TF}}+\ensuremath{\Delta}\phantom{\rule{0.2em}{0ex}}\mathrm{mod}(N,2)$ where ETF{E}_{\mathit{TF}}ETF​ is the Thomas-Fermi energy of a noninteracting gas in the trap and $\ensuremath{\Delta}$ is the pairing gap. There is no evidence of a shell correction energy in the systematics, but the density distributions show pronounced shell effects. We find the value $\ensuremath{\Delta}=0.7\ifmmode\pm\else\textpm\fi{}0.2\ensuremath{\omega}$ for the pairing gap. This is smaller than the value found for the uniform gas at a density corresponding to the central density of the trapped gas.
Using an anyon-fermion mapping method, we investigate the ground-state properties of hard-core anyons confined in a one-dimensional harmonic trap. The concise analytical formula of the reduced one-body density matrix is obtained. Based on this formula, we evaluated the momentum distribution, the natural orbitals, and their occupation distributions for different statistical parameters. The occupation and occupation fraction of the lowest natural orbital versus anyon number are also displayed. It is shown that the ground-state properties of anyons interplay between bosons and fermions continuously. We can expect that the hard-core anyons of larger statistical parameters exhibit similar properties to the hard-core bosons although anyon systems satisfy specific fractional statistics.
We investigate the properties of natural orbitals and their occupation numbers of the ground state of two noninteracting anyons characterized by the fractional statistics parameter α and confined in a harmonic trap. We work in the boson magnetic gauge where the anyons are modeled as composite bosons with magnetic flux quanta attached to their positions. We derive an asymptotic form of the weakly occupied natural orbitals, and show that their corresponding (ordered descendingly) occupation numbers decay according to the power law n−(4+2α) , where n is the index of the natural orbital. We find remarkable numerical agreement of the theory with the natural orbitals and their occupation numbers computed from the spectral decomposition of the system's wave function. We explain that the same results apply to the fermion magnetic gauge. Published by the American Physical Society 2024