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:1611.04757 · Physical Review A

Multicomponent correlated-basis-function method and its application to multilayered dipolar Bose gases

Michael Rader, Martin Hebenstreit, Robert E. Zillich

Published Mar 24, 2017• 6 citations
arXivHTMLPDF

Abstract

We present a method for calculating the dynamics of a bosonic mixture, the multicomponent correlated-basis-function (CBF) method. For a single component, CBF results for the excitation energies agree quite well with experimental results, even for highly correlated systems like $^{4}\mathrm{He}$, and recent systematic improvements of CBF achieve perfect agreement. We give a full derivation of multicomponent CBF, and apply the method to a dipolar Bose gas cut into two-dimensional layers by a deep optical lattice, with coupling between layers due to the long-ranged dipole-dipole interaction. We consider the case of strong coupling, leading to large positive interlayer correlations. We calculate the spectrum for a system of eight layers and show that the strong coupling can lead to a simpler spectrum than in the uncoupled case, with a single peak carrying most of the spectral weight.

DOI
10.1103/physreva.95.033625

Related papers

Stability and excitations of a bilayer of strongly correlated dipolar bosons

Diana Hufnagl, Robert E. Zillich

We study correlation effects and excitations in a dipolar Bose gas bilayer which is modeled by a one-dimensional double-well trap that determines the width of an individual layer, the distance between the two layers, and the height of the barrier between them. For the ground-state calculations we use the hypernetted-chain Euler-Lagrange method, and for the calculation of the excitations we use the correlated basis function method. We observe instabilities both for wide, well-separated layers dominated by intralayer attraction of the dipoles, and for narrow layers that are close to each other dominated by interlayer attraction. The behavior of the pair distribution function leads to the interpretation that the monomer phase becomes unstable when pairing of two dipoles becomes energetically favorable between or within layers, respectively. In both cases we observe a tendency towards ``rotonization,'' i.e., the appearance of a soft mode with finite momentum in the excitation spectrum. The dynamic structure function is not simply characterized by a single excitation mode but has a nontrivial multipeak structure that is not captured by the Bijl-Feynman approximation. The dipole-dipole interaction between different layers leads to additional damping compared to the damping obtained for uncoupled layers.

arXiv:1301.0343•Physical Review A•Mar 25, 2013•16 citations

Diagrammatic Analysis of the Method of Correlated Basis Functions. III. Excitation Spectrum of a Weakly Interacting Bose Gas

Archana Bhattacharyya, Chia-Wei Woo

In the continuing effort to establish contacts between the method of correlated basis functions (CBF) and the diagrammatic perturbation theory, we carry out parallel calculations of the excitation spectrum of a weakly interacting Bose gas. The Hugenholtz---Pines (HP) theory is first used to obtain the exact liquid-structure function and the energy spectrum, each to three orders in the strength parameter $\ensuremath{\lambda}$ that characterizes the interaction between a pair of particles in the system. Next we construct sets of CBF using three different correlating factors, and in each case calculate the excitation energy of the system. The results are compared to those obtained from the HP theory. Complete agreement is found to all orders considered. In conjunction with a paper to follow, the analysis demonstrates the flexibility of the CBF method and the relative ease of application in comparison to other methods.

Physical Review A•Jan 1, 1973•10 citations

Dynamics of a Two-Dimensional System of Quantum Dipoles

Ferran Mazzanti, Robert E. Zillich, G. E. Astrakharchik +1 more

A detailed microscopic analysis of the dynamic structure function S(k,omega) of a two-dimensional Bose system of dipoles polarized along the direction perpendicular to the plane is presented and discussed. Starting from ground-state quantities obtained using a quantum diffusion Monte Carlo algorithm, the density-density response is evaluated in the context of the correlated basis functions (CBF) theory. CBF predicts a sharp peak and a multiexcitation component at higher energies produced by the decay of excitations. We discuss the structure of the phonon-roton peak and show that the Feynman and Bogoliubov predictions depart from the CBF result already at low densities. We finally discuss the emergence of a roton in the spectrum, but find the roton energy not low enough to make the system unstable under density fluctuations up to the highest density considered that is close to the freezing point.

Physical Review Letters•Mar 20, 2009•38 citations

Excitation Spectrum and Liquid-Structure Function of a Charged Bose Gas

A. Bhattacharyya, Chia-Wei Woo

The excitation spectrum of a charged Bose gas in the high-density limit is calculated to two orders in the expansion parameter ${r}_{s}^{\frac{3}{4}}$ by the method of correlated basis functions (CBF). To $O({r}_{s}^{\frac{3}{4}})$, the result differs from that obtained variationally, but agrees with that obtained by Ma and Woo via diagrammatic perturbation theory. The choice of the basis functions is, however, not unique. By choosing an alternate set of CBF which employs the ground-state eigenfunction as a correlating factor, we find it possible to extract from the excitation spectrum an expression for the liquid-structure function, also to two orders in ${r}_{s}^{\frac{3}{4}}$, which represents a new result. This new facet of CBF helps to demonstrate the flexibility of the method and the ease with which it can be applied to calculating properties of many-particle systems, model as well as realistic.

Physical Review A•Jan 1, 1973•11 citations

Correlation effects and collective excitations in bosonic bilayers: Role of quantum statistics, superfluidity, and the dimerization transition

Alexey V. Filinov

A two-component, two-dimensional (2D) dipolar bosonic system in the bilayer geometry is considered. By performing quantum Monte Carlo simulations in a wide range of layer spacings we analyze in detail the pair correlation functions, the static response function, and the kinetic and interaction energies. By reducing the layer spacing we observe a transition from weakly to strongly bound dimer states. The transition is accompanied by the onset of short-range correlations, suppression of the superfluid response, and rotonization of the excitation spectrum. A dispersion law and a dynamic structure factor for the in-phase (symmetric) and out-of-phase (antisymmetric) collective modes during the dimerization is studied in detail with the stochastic reconstruction method and the method of moments. The antisymmetric mode spectrum is most strongly influenced by suppression of the inlayer superfluidity (specified by the superfluid fraction ${\ensuremath{\gamma}}_{s}={\ensuremath{\rho}}_{s}/\ensuremath{\rho}$). In a pure superfluid (normal fluid) phase, only an acoustic [optical (gapped)] mode is recovered. In a partially superfluid phase, both are present simultaneously, and the dispersion splits into two branches corresponding to a normal and a superfluid component. The spectral weight of the acoustic mode scales linearly with ${\ensuremath{\gamma}}_{s}$. This weight transfers to the optical branch when ${\ensuremath{\gamma}}_{s}$ is reduced due to formation of dimer states. In summary, we demonstrate how the interlayer dimerization in dipolar bilayers can be uniquely identified by static and dynamic properties.

arXiv:1604.08110•Physical Review A•Jul 11, 2016•36 citations

Comparison between the correlated-basis-functions method and the density-functional theory for inhomogeneous systems

Zhao-Qing Zhang, Ruibao Tao, Chia-Wei Woo

As a hypothetical inhomogeneous quantum many-body model, we take a weakly interacting Bose gas made inhomogeneous by the presence of a weak external field. The interaction is assumed to be Fourier transformable and characterized by a strength parameter \ensuremath{\lambda}; this is also true for the external field, characterized by a strength parameter \ensuremath{\varepsilon}. The ground-state energy is calculated to order ${\ensuremath{\varepsilon}}^{2}$${\ensuremath{\lambda}}^{2}$ using three methods: diagrammatic perturbation theory, correlated basis functions (CBF), and density-functional theory (DFT). The diagrammatic perturbation theoretical calculation is new. It differs from our earlier approach for homogeneous Bose systems which evolved from the Hugenholtz-Pines theory. We are able to reproduce the same results in an efficient manner. For the inhomogeneous system, the results are used to measure CBF against DFT. It is found that in their simplest forms CBF yields results which are better than the renormalized DFT by Ebner and Saam. In particular, CBF is superior to DFT in the local-density approximation, even though the latter was designed to take care of systems with slow and small density variations such as the model considered here.

Physical review. B, Condensed matter•Feb 1, 1988•3 citations
We present a method for calculating the dynamics of a bosonic mixture, the multicomponent correlated-basis-function (CBF) method. For a single component, CBF results for the excitation energies agree quite well with experimental results, even for highly correlated systems like 4He^{4}\mathrm{He}4He, and recent systematic improvements of CBF achieve perfect agreement. We give a full derivation of multicomponent CBF, and apply the method to a dipolar Bose gas cut into two-dimensional layers by a deep optical lattice, with coupling between layers due to the long-ranged dipole-dipole interaction. We consider the case of strong coupling, leading to large positive interlayer correlations. We calculate the spectrum for a system of eight layers and show that the strong coupling can lead to a simpler spectrum than in the uncoupled case, with a single peak carrying most of the spectral weight.
We study correlation effects and excitations in a dipolar Bose gas bilayer which is modeled by a one-dimensional double-well trap that determines the width of an individual layer, the distance between the two layers, and the height of the barrier between them. For the ground-state calculations we use the hypernetted-chain Euler-Lagrange method, and for the calculation of the excitations we use the correlated basis function method. We observe instabilities both for wide, well-separated layers dominated by intralayer attraction of the dipoles, and for narrow layers that are close to each other dominated by interlayer attraction. The behavior of the pair distribution function leads to the interpretation that the monomer phase becomes unstable when pairing of two dipoles becomes energetically favorable between or within layers, respectively. In both cases we observe a tendency towards ``rotonization,'' i.e., the appearance of a soft mode with finite momentum in the excitation spectrum. The dynamic structure function is not simply characterized by a single excitation mode but has a nontrivial multipeak structure that is not captured by the Bijl-Feynman approximation. The dipole-dipole interaction between different layers leads to additional damping compared to the damping obtained for uncoupled layers.
In the continuing effort to establish contacts between the method of correlated basis functions (CBF) and the diagrammatic perturbation theory, we carry out parallel calculations of the excitation spectrum of a weakly interacting Bose gas. The Hugenholtz---Pines (HP) theory is first used to obtain the exact liquid-structure function and the energy spectrum, each to three orders in the strength parameter $\ensuremath{\lambda}$ that characterizes the interaction between a pair of particles in the system. Next we construct sets of CBF using three different correlating factors, and in each case calculate the excitation energy of the system. The results are compared to those obtained from the HP theory. Complete agreement is found to all orders considered. In conjunction with a paper to follow, the analysis demonstrates the flexibility of the CBF method and the relative ease of application in comparison to other methods.
A detailed microscopic analysis of the dynamic structure function S(k,omega) of a two-dimensional Bose system of dipoles polarized along the direction perpendicular to the plane is presented and discussed. Starting from ground-state quantities obtained using a quantum diffusion Monte Carlo algorithm, the density-density response is evaluated in the context of the correlated basis functions (CBF) theory. CBF predicts a sharp peak and a multiexcitation component at higher energies produced by the decay of excitations. We discuss the structure of the phonon-roton peak and show that the Feynman and Bogoliubov predictions depart from the CBF result already at low densities. We finally discuss the emergence of a roton in the spectrum, but find the roton energy not low enough to make the system unstable under density fluctuations up to the highest density considered that is close to the freezing point.
The excitation spectrum of a charged Bose gas in the high-density limit is calculated to two orders in the expansion parameter rs34{r}_{s}^{\frac{3}{4}}rs43​​ by the method of correlated basis functions (CBF). To O(rs34)O({r}_{s}^{\frac{3}{4}})O(rs43​​), the result differs from that obtained variationally, but agrees with that obtained by Ma and Woo via diagrammatic perturbation theory. The choice of the basis functions is, however, not unique. By choosing an alternate set of CBF which employs the ground-state eigenfunction as a correlating factor, we find it possible to extract from the excitation spectrum an expression for the liquid-structure function, also to two orders in rs34{r}_{s}^{\frac{3}{4}}rs43​​, which represents a new result. This new facet of CBF helps to demonstrate the flexibility of the method and the ease with which it can be applied to calculating properties of many-particle systems, model as well as realistic.
A two-component, two-dimensional (2D) dipolar bosonic system in the bilayer geometry is considered. By performing quantum Monte Carlo simulations in a wide range of layer spacings we analyze in detail the pair correlation functions, the static response function, and the kinetic and interaction energies. By reducing the layer spacing we observe a transition from weakly to strongly bound dimer states. The transition is accompanied by the onset of short-range correlations, suppression of the superfluid response, and rotonization of the excitation spectrum. A dispersion law and a dynamic structure factor for the in-phase (symmetric) and out-of-phase (antisymmetric) collective modes during the dimerization is studied in detail with the stochastic reconstruction method and the method of moments. The antisymmetric mode spectrum is most strongly influenced by suppression of the inlayer superfluidity (specified by the superfluid fraction ${\ensuremath{\gamma}}_{s}={\ensuremath{\rho}}_{s}/\ensuremath{\rho}$). In a pure superfluid (normal fluid) phase, only an acoustic [optical (gapped)] mode is recovered. In a partially superfluid phase, both are present simultaneously, and the dispersion splits into two branches corresponding to a normal and a superfluid component. The spectral weight of the acoustic mode scales linearly with ${\ensuremath{\gamma}}_{s}$. This weight transfers to the optical branch when ${\ensuremath{\gamma}}_{s}$ is reduced due to formation of dimer states. In summary, we demonstrate how the interlayer dimerization in dipolar bilayers can be uniquely identified by static and dynamic properties.
As a hypothetical inhomogeneous quantum many-body model, we take a weakly interacting Bose gas made inhomogeneous by the presence of a weak external field. The interaction is assumed to be Fourier transformable and characterized by a strength parameter \ensuremath{\lambda}; this is also true for the external field, characterized by a strength parameter \ensuremath{\varepsilon}. The ground-state energy is calculated to order ${\ensuremath{\varepsilon}}^{2}$${\ensuremath{\lambda}}^{2}$ using three methods: diagrammatic perturbation theory, correlated basis functions (CBF), and density-functional theory (DFT). The diagrammatic perturbation theoretical calculation is new. It differs from our earlier approach for homogeneous Bose systems which evolved from the Hugenholtz-Pines theory. We are able to reproduce the same results in an efficient manner. For the inhomogeneous system, the results are used to measure CBF against DFT. It is found that in their simplest forms CBF yields results which are better than the renormalized DFT by Ebner and Saam. In particular, CBF is superior to DFT in the local-density approximation, even though the latter was designed to take care of systems with slow and small density variations such as the model considered here.