-
Suppression of polaron self-localization by correlations
Authors:
Lilith Zschetzsche,
Robert E. Zillich
Abstract:
We investigate self-localization of a polaron in a homogeneous Bose-Einstein condensate in one dimension. This effect, where an impurity is trapped by the deformation that it causes in the surrounding Bose gas, has been first predicted by mean field calculations, but has not been seen in experiments. We study the system in one dimension, where, according to the mean field approximation, the self-l…
▽ More
We investigate self-localization of a polaron in a homogeneous Bose-Einstein condensate in one dimension. This effect, where an impurity is trapped by the deformation that it causes in the surrounding Bose gas, has been first predicted by mean field calculations, but has not been seen in experiments. We study the system in one dimension, where, according to the mean field approximation, the self-localization effect is particularly robust, and present for arbitrarily weak impurity-boson interactions. We address the question whether self-localization is a real effect by developing a variational method which incorporates impurity-boson correlations non-perturbatively and solving the resulting inhomogeneous correlated polaron equations. We find that correlations inhibit self-localization except for very strongly repulsive or attractive impurity-boson interactions. Our prediction for the critical interaction strength for self-localization agrees with a sharp drop of the inverse effective mass to almost zero found in quantum Monte Carlo simulations of polarons in one dimension.
△ Less
Submitted 18 December, 2023;
originally announced December 2023.
-
Nonadiabatic Laser-Induced Alignment Dynamics of Molecules on a Surface
Authors:
Lorenz Kranabetter,
Henrik H. Kristensen,
Areg Ghazaryan,
Constant A. Schouder,
Adam S. Chatterley,
Paul Janssen,
Frank Jensen,
Robert E. Zillich,
Mikhail Lemeshko,
Henrik Stapelfeldt
Abstract:
We demonstrate that a sodium dimer, Na$_{2}$($1^3Σ_{u}^+$), residing on the surface of a helium nanodroplet, can be set into rotation by a nonresonant 1.0 ps infrared laser pulse. The time-dependent degree of alignment measured, exhibits a periodic, gradually decreasing structure that deviates qualitatively from that expected for gas phase dimers. Comparison to alignment dynamics calculated from t…
▽ More
We demonstrate that a sodium dimer, Na$_{2}$($1^3Σ_{u}^+$), residing on the surface of a helium nanodroplet, can be set into rotation by a nonresonant 1.0 ps infrared laser pulse. The time-dependent degree of alignment measured, exhibits a periodic, gradually decreasing structure that deviates qualitatively from that expected for gas phase dimers. Comparison to alignment dynamics calculated from the time-dependent rotational Schrödinger equation shows that the deviation is due to the alignment dependent interaction between the dimer and the droplet surface. This interaction confines the dimer to the tangential plane of the droplet surface at the point where it resides and is the reason that the observed alignment dynamics is also well-described by a 2D quantum rotor model.
△ Less
Submitted 29 August, 2023;
originally announced August 2023.
-
Molecular Dynamics Study of the Sonic Horizon of Microscopic Laval Nozzles
Authors:
Helmut Ortmayer,
Robert E. Zillich
Abstract:
A Laval nozzle can accelerate expanding gas above supersonic velocities, while cooling the gas in the process. This work investigates this process for microscopic Laval nozzles by means of non-equilibrium molecular dynamics simulations of statioary flow, using grand canonical Monte-Carlo particle reservoirs. We study the expansion of a simple fluid, a mono-atomic gas interacting via a Lennard-Jone…
▽ More
A Laval nozzle can accelerate expanding gas above supersonic velocities, while cooling the gas in the process. This work investigates this process for microscopic Laval nozzles by means of non-equilibrium molecular dynamics simulations of statioary flow, using grand canonical Monte-Carlo particle reservoirs. We study the expansion of a simple fluid, a mono-atomic gas interacting via a Lennard-Jones potential, through an idealized nozzle with atomically smooth walls. We obtain the thermodynamic state variables pressure, density, and temperature, but also the Knudsen number, speed of sound, velocity, and the corresponing Mach number of the expanding gas for nozzles of different sizes. We find that the temperature is well-defined in the sense that the each velocity components of the particles obey the Maxwell-Boltzmann distribution, but it is anisotropic, especially for small nozzles. The velocity auto-correlation function reveals a tendency towards condensation of the cooled supersonic gas, although the nozzles are too small for the formation of clusters. Overall we find that microscopic nozzles act qualitatively like macroscopic nozzles in that the particles are accelerated to supersonic speeds while their thermal motion relative to the stationary flow is cooled. We find that, like macroscopic Laval nozzles, microscopic nozzles also exhibit a sonic horizon, which is well-defined on a microscopic scale. The sonic horizon is positioned only slightly further downstream compared to isentropic expansion through macroscopic nozzles, where it is situated in the most narrow part. We analyze the sonic horizon by studying spacetime density correlations, i.e.\ how thermal fluctuations at two positions of the gas density are correlated in time and find that after the sonic horizon there are indeed no upstream correlations on a microscopic scale.
△ Less
Submitted 30 June, 2023;
originally announced June 2023.
-
Interaction quenches in Bose gases studied with a time-dependent hypernetted-chain Euler-Lagrange method
Authors:
Mathias Gartner,
David Miesbauer,
Michael Kobler,
Julia Freund,
Giuseppe Carleo,
Robert E. Zillich
Abstract:
We present a new variational method to study the dynamics of a closed bosonic many-body system, the time-dependent hypernetted-chain Euler-Lagrange method, tHNC . Based on the Jastrow ansatz, it accounts for quantum fluctuations in a non-perturbative way. tHNC scales well with the number of dimensions, as demonstrated by our results on one, two, and three dimensions. We apply the tHNC method to in…
▽ More
We present a new variational method to study the dynamics of a closed bosonic many-body system, the time-dependent hypernetted-chain Euler-Lagrange method, tHNC . Based on the Jastrow ansatz, it accounts for quantum fluctuations in a non-perturbative way. tHNC scales well with the number of dimensions, as demonstrated by our results on one, two, and three dimensions. We apply the tHNC method to interaction quenches, i.e. sudden changes of the interaction strength, in homogeneous Bose gases. When the quench is strong enough that the final state has roton excitations (as found and predicted for dipolar and Rydberg-dressed Bose-Einstein condensates, respectively), the pair distribution function exhibits stable oscillations. For validation, we compare tHNC results with time-dependent variational Monte Carlo results in one and two dimensions.
△ Less
Submitted 14 December, 2022;
originally announced December 2022.
-
Striped Ultradilute Liquid of Dipolar Bosons in Two Dimensions
Authors:
Clemens Staudinger,
Diana Hufnagl,
Ferran Mazzanti,
Robert E. Zillich
Abstract:
We investigate the phases of a Bose-Einstein condensate of dipolar atoms restricted to move in a two-dimensional plane. The dipole moments are all aligned in a direction tilted with respect to the plane normal. As a result of the attractive and repulsive components of the dipole-dipole interaction, the dipolar gas has a self-bound phase, which is stabilized by quantum fluctuations. Furthermore, ti…
▽ More
We investigate the phases of a Bose-Einstein condensate of dipolar atoms restricted to move in a two-dimensional plane. The dipole moments are all aligned in a direction tilted with respect to the plane normal. As a result of the attractive and repulsive components of the dipole-dipole interaction, the dipolar gas has a self-bound phase, which is stabilized by quantum fluctuations. Furthermore, tilting the dipoles tunes the anisotropy of the dipole-dipole interaction, which can trigger a spatial density modulation. In this work we study these two aspects and investigate the conditions for the formation of a self-bound and striped phase, which has been realized in experiments with dipolar droplets. We use a variational method based on the hypernetted-chain Euler-Lagrange optimization of a Jastrow-Feenberg ansatz for the many-body wave function to study the ground state properties. This method takes into account quantum fluctuations in a non-perturbative way and thus can be used also for strongly correlated systems.
△ Less
Submitted 20 March, 2023; v1 submitted 9 August, 2022;
originally announced August 2022.
-
Time-dependent variational Monte Carlo study of the dynamic response of bosons in an optical lattice
Authors:
Mathias Gartner,
Ferran Mazzanti,
Robert E. Zillich
Abstract:
We study the dynamics of a one-dimensional Bose gas at unit filling in both shallow and deep optical lattices and obtain the dynamic structure factor ${S(k,ω)}$ by monitoring the linear response to a weak probe pulse. We introduce a new procedure, based on the time-dependent variational Monte Carlo method (tVMC), which allows to evolve the system in real time, using as a variational model a Jastro…
▽ More
We study the dynamics of a one-dimensional Bose gas at unit filling in both shallow and deep optical lattices and obtain the dynamic structure factor ${S(k,ω)}$ by monitoring the linear response to a weak probe pulse. We introduce a new procedure, based on the time-dependent variational Monte Carlo method (tVMC), which allows to evolve the system in real time, using as a variational model a Jastrow-Feenberg wave function that includes pair correlations. Comparison with exact diagonalization results of ${S(k,ω)}$ obtained on a lattice in the Bose-Hubbard limit shows good agreement of the dispersion relation for sufficiently deep optical lattices, while for shallow lattices we observe the influence of higher Bloch bands. We also investigate non-linear response to strong pulses. From the power spectrum of the density fluctuations we obtain the excitation spectrum, albeit broadened, by higher harmonic generation after a strong pulse with a single low wave number. As a remarkable feature of our simulations we furthermore demonstrate that the full excitation spectrum can be retrieved from the power spectrum of the density fluctuations due to the stochastic noise inherent in any Monte Carlo method, without applying an actual perturbation.
△ Less
Submitted 27 May, 2022; v1 submitted 15 February, 2022;
originally announced February 2022.
-
Rotational coherence spectroscopy of molecules in helium nanodroplets: Reconciling the time and the frequency domains
Authors:
Adam S. Chatterley,
Lars Christiansen,
Constant A. Schouder,
Anders V. Jørgensen,
Benjamin Shepperson,
Igor N. Cherepanov,
Giacomo Bighin,
Robert E. Zillich,
Mikhail Lemeshko,
Henrik Stapelfeldt
Abstract:
Alignment of OCS, CS$_2$ and I$_2$ molecules embedded in helium nanodroplets is measured as a function of time following rotational excitation by a non-resonant, comparatively weak ps laser pulse. The distinct peaks in the power spectra, obtained by Fourier analysis, are used to determine the rotational, B, and centrifugal distortion, D, constants. For OCS, B and D match the values known from IR s…
▽ More
Alignment of OCS, CS$_2$ and I$_2$ molecules embedded in helium nanodroplets is measured as a function of time following rotational excitation by a non-resonant, comparatively weak ps laser pulse. The distinct peaks in the power spectra, obtained by Fourier analysis, are used to determine the rotational, B, and centrifugal distortion, D, constants. For OCS, B and D match the values known from IR spectroscopy. For CS$_2$ and I$_2$, they are the first experimental results reported. The alignment dynamics calculated from the gas-phase rotational Schrödinger equation, using the experimental in-droplet B and D values, agree in detail with the measurement for all three molecules. The rotational spectroscopy technique for molecules in helium droplets introduced here should apply to a range of molecules and complexes.
△ Less
Submitted 4 June, 2020;
originally announced June 2020.
-
A quantum heat engine based on dynamical materials design
Authors:
G. Tulzer,
M. Hoffmann,
R. E. Zillich
Abstract:
We propose a novel type of quantum heat engine based on the ultrafast dynamical control of the magnetic properties of a nano-scale working body. The working principle relies on nonlinear phononics, an example for dynamical materials design. We describe the general recipe for identifying candidate materials, and also propose Cr$_{2}$O$_{3}$ as a promising working body for a quantum Otto cycle. Usin…
▽ More
We propose a novel type of quantum heat engine based on the ultrafast dynamical control of the magnetic properties of a nano-scale working body. The working principle relies on nonlinear phononics, an example for dynamical materials design. We describe the general recipe for identifying candidate materials, and also propose Cr$_{2}$O$_{3}$ as a promising working body for a quantum Otto cycle. Using a spin Hamiltonian as a model for Cr$_{2}$O$_{3}$, we investigate the performance in terms of efficiency, output power, and quantum friction. To assess the assumptions underlying our effective spin Hamiltonian we also consider a working substance composed of several unit cells. We show that even without an implementation of transitionless driving, the quantum friction is very low compared to the total produced work and the energy cost of counterdiabatic driving is negligible. This is an advantage of the working substance, as experimentally hard-to-implement shortcuts to adiabaticity are not needed. Moreover, we discuss some remarkable thermodynamic features due to the quantumness of the proposed system such as a non-monotonic dependence of the efficiency on the temperature of the hot-bath. Finally, we explore the dependence of the performance on the system parameters for a generic model of this type of quantum heat engine and identify properties of the energy spectrum required for a well-performing quantum heat engine.
△ Less
Submitted 25 May, 2020; v1 submitted 13 September, 2019;
originally announced September 2019.
-
Rotons and Bose condensation in Rydberg-dressed Bose Gases
Authors:
Iran Seydi,
Saeed H. Abedinpour,
Robert E. Zillich,
Reza Asgari,
B. Tanatar
Abstract:
We investigate the ground-state properties and excitations of Rydberg-dressed bosons in both three and two dimensions, using the hypernetted-chain Euler-Lagrange approximation, which accounts for correlations and thus goes beyond the mean field approximation. The short-range behavior of the pair distribution function signals the instability of the homogeneous system towards the formation of drople…
▽ More
We investigate the ground-state properties and excitations of Rydberg-dressed bosons in both three and two dimensions, using the hypernetted-chain Euler-Lagrange approximation, which accounts for correlations and thus goes beyond the mean field approximation. The short-range behavior of the pair distribution function signals the instability of the homogeneous system towards the formation of droplet crystals at strong couplings and large soft-core radius. This tendency to spatial density modulation coexists with off-diagonal long-range order. The contribution of the correlation energy to the ground-state energy is significant at large coupling strengths and intermediate values of the soft-core radius while for a larger soft-core radius the ground-state energy is dominated by the mean-field (Hartree) energy. We have also performed path integral Monte Carlo simulations to verify the performance of our hypernetted-chain Euler-Lagrange results in three dimensions. In the homogeneous phase, the two approaches are in very good agreement. Moreover, Monte Carlo simulations predict a first-order quantum phase transition from a homogeneous superfluid phase to the quantum droplet phase with face-centered cubic symmetry for Rydberg-dressed bosons in three dimensions.
△ Less
Submitted 30 March, 2020; v1 submitted 5 May, 2019;
originally announced May 2019.
-
Self-bound Bose mixtures
Authors:
Clemens Staudinger,
Ferran Mazzanti,
Robert E. Zillich
Abstract:
Recent experiments confirmed that fluctuations beyond the mean-field approximation can lead to self-bound liquid droplets of ultra-dilute binary Bose mixtures. We proceed beyond the beyond-mean-field approximation, and study liquid Bose mixtures using the variational hypernetted-chain Euler Lagrange method, which accounts for correlations non-perturbatively. Focusing on the case of a mixture of un…
▽ More
Recent experiments confirmed that fluctuations beyond the mean-field approximation can lead to self-bound liquid droplets of ultra-dilute binary Bose mixtures. We proceed beyond the beyond-mean-field approximation, and study liquid Bose mixtures using the variational hypernetted-chain Euler Lagrange method, which accounts for correlations non-perturbatively. Focusing on the case of a mixture of uniform density, as realized inside large saturated droplets, we study the conditions for stability against evaporation of one of the components (both chemical potentials need to be negative) and against liquid-gas phase separation (spinodal instability), the latter being accompanied by a vanishing speed of sound. Dilute Bose mixtures are stable only in a narrow range near an optimal ratio $ρ_1/ρ_2$ and near the total energy minimum. Deviations from a universal dependence on the s-wave scattering lengths are significant despite the low density.
△ Less
Submitted 16 May, 2018;
originally announced May 2018.
-
Quantum phases of dipolar rotors on two-dimensional lattices
Authors:
Brendan P. Abolins,
Robert E. Zillich,
K. Birgitta Whaley
Abstract:
The quantum phase transitions of dipoles confined to the vertices of two dimensional (2D) lattices of square and triangular geometry is studied using path integral ground state quantum Monte Carlo (PIGS). We analyze the phase diagram as a function of the strength of both the dipolar interaction and a transverse electric field. The study reveals the existence of a class of orientational phases of q…
▽ More
The quantum phase transitions of dipoles confined to the vertices of two dimensional (2D) lattices of square and triangular geometry is studied using path integral ground state quantum Monte Carlo (PIGS). We analyze the phase diagram as a function of the strength of both the dipolar interaction and a transverse electric field. The study reveals the existence of a class of orientational phases of quantum dipolar rotors whose properties are determined by the ratios between the strength anisotropic dipole-dipole interaction, the strength of the applied transverse field, and the rotational constant. For the triangular lattice, the generic orientationally disordered phase found at zero and weak values of both dipolar interaction strength and applied field, is found to show a transition to a phase characterized by net polarization in the lattice plane as the strength of the dipole-dipole interaction is increased, independent of the strength of the applied transverse field, in addition to the expected transition to a transverse polarized phase as the electric field strength increases. The square lattice is also found to exhibit a transition from a disordered phase to an ordered phase as the dipole-dipole interaction strength is increased, as well as the expected transition to a transverse polarized phase as the electric field strength increases. In contrast to the situation with a triangular lattice, on square lattices the ordered phase at high dipole-dipole interaction strength possesses a striped ordering. The properties of these quantum dipolar rotor phases are dominated by the anisotropy of the interaction and provide useful models for developing quantum phases beyond the well-known paradigms of spin Hamiltonian models, realizing in particular a novel physical realization of a quantum rotor-like Hamiltonian that possesses an anisotropic long range interaction.
△ Less
Submitted 8 March, 2018; v1 submitted 6 March, 2018;
originally announced March 2018.
-
Laser-induced rotation of iodine molecules in He-nanodroplets: revivals and breaking-free
Authors:
Benjamin Shepperson,
Anders A. Søndergaard,
Lars Christiansen,
Jan Kaczmarczyk,
Robert E. Zillich,
Mikhail Lemeshko,
Henrik Stapelfeldt
Abstract:
Rotation of molecules embedded in He nanodroplets is explored by a combination of fs laser-induced alignment experiments and angulon quasiparticle theory. We demonstrate that at low fluence of the fs alignment pulse, the molecule and its solvation shell can be set into coherent collective rotation lasting long enough to form revivals. With increasing fluence, however, the revivals disappear -- ins…
▽ More
Rotation of molecules embedded in He nanodroplets is explored by a combination of fs laser-induced alignment experiments and angulon quasiparticle theory. We demonstrate that at low fluence of the fs alignment pulse, the molecule and its solvation shell can be set into coherent collective rotation lasting long enough to form revivals. With increasing fluence, however, the revivals disappear -- instead, rotational dynamics as rapid as for an isolated molecule is observed during the first few picoseconds. Classical calculations trace this phenomenon to transient decoupling of the molecule from its He shell. Our results open novel opportunities for studying non-equilibrium solute-solvent dynamics and quantum thermalization.
△ Less
Submitted 7 February, 2017;
originally announced February 2017.
-
The Multi-component Correlated Basis Function Method and its Application to Multilayered Dipolar Bose Gases
Authors:
Michael Rader,
Martin Hebenstreit,
Robert E. Zillich
Abstract:
We present a method for calculating the dynamics of a bosonic mixture, the multi-component correlated basis function (CBF) method. For single components, CBF results for the excitation energies agree quite well with experimental results, even for highly correlated systems like $^4$He, and recent systematic improvements of CBF achieve perfect agreement. We give a full derivation of multi-component…
▽ More
We present a method for calculating the dynamics of a bosonic mixture, the multi-component correlated basis function (CBF) method. For single components, CBF results for the excitation energies agree quite well with experimental results, even for highly correlated systems like $^4$He, and recent systematic improvements of CBF achieve perfect agreement. We give a full derivation of multi-component 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 8 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.
△ Less
Submitted 15 November, 2016;
originally announced November 2016.
-
Solvation of Mg in Helium-4: Are there Meta-stable Mg Dimers ?
Authors:
Eckhard Krotscheck,
Robert E. Zillich
Abstract:
Experiments on the formation of magnesium complexes in $^4$He nanodroplets were interpreted as the observation of the formation of weakly bound magnesium complexes. We present results for single Mg and Mg dimer solvation using the hypernetted chain / Euler-Lagrange method as well as path integral Monte Carlo simulations. We find that the phonon-mediated, indirect Mg-Mg interaction adds an oscillat…
▽ More
Experiments on the formation of magnesium complexes in $^4$He nanodroplets were interpreted as the observation of the formation of weakly bound magnesium complexes. We present results for single Mg and Mg dimer solvation using the hypernetted chain / Euler-Lagrange method as well as path integral Monte Carlo simulations. We find that the phonon-mediated, indirect Mg-Mg interaction adds an oscillatory component to the direct Mg-Mg interaction. We undertake a step-by-step examination of the ingredients of the calculation of the phonon-induced interaction, comparing the results of semi-analytic HNC-EL calculations for bulk and single impurity results with experiments as well as Monte Carlo data. We do not find evidence for a sufficiently strong secondary minimum in the effective Mg-Mg interaction to support a metastable state.
△ Less
Submitted 14 October, 2016;
originally announced October 2016.
-
Gapped spectrum in pair-superfluid bosons
Authors:
G. E. Astrakharchik,
R. E. Zillich,
F. Mazzanti,
J. Boronat
Abstract:
We study the ground state of a bilayer system of dipolar bosons with dipoles oriented by an external field perpendicularly to the two parallel planes. By decreasing the interlayer distance, for a fixed value of the strength of the dipolar interaction, the system undergoes a quantum phase transition from an atomic to a pair superfluid. We investigate the excitation spectrum across this transition b…
▽ More
We study the ground state of a bilayer system of dipolar bosons with dipoles oriented by an external field perpendicularly to the two parallel planes. By decreasing the interlayer distance, for a fixed value of the strength of the dipolar interaction, the system undergoes a quantum phase transition from an atomic to a pair superfluid. We investigate the excitation spectrum across this transition by using microscopic approaches. Quantum Monte Carlo methods are employed to obtain the static structure factors and intermediate scattering functions in imaginary time. The dynamic response is calculated using both the correlated basis functions (CBF) method and the approximate inversion of the Laplace transform of the quantum Monte Carlo imaginary time data. In the atomic phase, both density and spin excitations are gapless. However, in the pair-superfluid phase a gap opens in the excitation energy of the spin mode. For small separation between layers, the minimal spin excitation energy equals the binding energy of a dimer and is twice the gap value.
△ Less
Submitted 25 August, 2016;
originally announced August 2016.
-
Dipolar Bilayer with Antiparallel Polarization -- a Self-Bound Liquid
Authors:
Martin Hebenstreit,
Michael Rader,
Robert E. Zillich
Abstract:
Dipolar bilayers with antiparallel polarization, i.e. opposite polarization in the two layers, exhibit liquid-like rather than gas-like behavior. In particular, even without external pressure a self-bound liquid puddle of constant density will form. We investigate the symmetric case of two identical layers, corresponding to a two-component Bose system with equal partial densities. The zero-tempera…
▽ More
Dipolar bilayers with antiparallel polarization, i.e. opposite polarization in the two layers, exhibit liquid-like rather than gas-like behavior. In particular, even without external pressure a self-bound liquid puddle of constant density will form. We investigate the symmetric case of two identical layers, corresponding to a two-component Bose system with equal partial densities. The zero-temperature equation of state $E(ρ)/N$, where $ρ$ is the total density, has a minimum, with an equilibrium density that decreases with increasing distance between the layers. The attraction necessary for a self-bound liquid comes from the inter-layer dipole-dipole interaction that leads to a mediated intra-layer attraction. We investigate the regime of negative pressure towards the spinodal instability, where the bilayer is unstable against infinitesimal fluctuations of the total density, conformed by calculations of the speed of sound of total density fluctuations.
△ Less
Submitted 9 January, 2015;
originally announced January 2015.
-
Correlations in the low-density Fermi gas: Fermi-Liquid state, Dimerization, and BCS Pairing
Authors:
H. H. Fan,
E. Krotscheck,
T. Lichtenegger,
D. Mateo,
R. E. Zillich
Abstract:
We present ground state calculations for low-density Fermi gases described by two model interactions, an attractive square-well potential and a Lennard-Jones potential, of varying strength. We use the optimized Fermi-Hypernetted Chain integral equation method which has been proved to provide, in the density regimes of interest here, an accuracy better than one percent. We first examine the low-den…
▽ More
We present ground state calculations for low-density Fermi gases described by two model interactions, an attractive square-well potential and a Lennard-Jones potential, of varying strength. We use the optimized Fermi-Hypernetted Chain integral equation method which has been proved to provide, in the density regimes of interest here, an accuracy better than one percent. We first examine the low-density expansion of the energy and compare with the exact answer by Huang and Yang (H. Huang and C. N. Yang, {\em Phys. Rev.\/} {\bf 105}, 767 (1957)). It is shown that a locally correlated wave function of the Jastrow-Feenberg type does not recover the quadratic term in the expansion of the energy in powers of $\a0\KF$, where $\a0$ is the vacuum $s$-wave scattering length and $\KF$ the Fermi wave number. The problem is cured by adding second-order perturbation corrections in a correlated basis. Going to higher densities and/or more strongly coupled systems, we encounter an instability of the normal state of the system which is characterized by a divergence of the {\em in-medium\/} scattering length. We interpret this divergence as a phonon-exchange driven dimerization of the system, similar to what one has at zero density when the vacuum scattering length $\a0$ diverges. We then study, in the stable regime, the superfluid gap and its dependence on the density and the interaction strength. We identify two different corrections to low-density expansions: One is medium corrections to the pairing interaction, and the other one finite-range corrections. We show that the most important finite-range corrections are a direct manifestation of the many-body nature of the system.
△ Less
Submitted 3 November, 2016; v1 submitted 19 September, 2014;
originally announced September 2014.
-
Stability and excitations of a bilayer of strongly correlated dipolar Bosons
Authors:
Diana Hufnagl,
Robert E. Zillich
Abstract:
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…
▽ More
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 \emph{intra-layer} attraction of the dipoles, and for narrow layers that are close to each other dominated by \emph{inter-layer} 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 non-trivial multi-peak 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.
△ Less
Submitted 3 March, 2013; v1 submitted 2 January, 2013;
originally announced January 2013.
-
Excitations and Stripe Phase Formation in a 2D Dipolar Bose Gas with Tilted Polarization
Authors:
Adrian Macia,
Diana Hufnagl,
Ferran Mazzanti,
Jordi Boronat,
Robert E. Zillich
Abstract:
We present calculations of the ground state and excitations of an anisotropic dipolar Bose gas in two dimensions, realized by a non-perpendicular polarization with respect to the system plane. For sufficiently high density an increase of the polarization angle leads to a density instability of the gas phase in the direction where the anisotropic interaction is strongest. Using a dynamic many-body…
▽ More
We present calculations of the ground state and excitations of an anisotropic dipolar Bose gas in two dimensions, realized by a non-perpendicular polarization with respect to the system plane. For sufficiently high density an increase of the polarization angle leads to a density instability of the gas phase in the direction where the anisotropic interaction is strongest. Using a dynamic many-body theory, we calculate the dynamic structure function in the gas phase which shows the anisotropic dispersion of the excitations. We find that the energy of roton excitations in the strongly interacting direction decreases with increasing polarization angle and almost vanishes close to the instability. Exact path integral ground state Monte Carlo simulations show that this instability is indeed a quantum phase transition to a stripe phase, characterized by long-range order in the strongly interacting direction.
△ Less
Submitted 7 September, 2012;
originally announced September 2012.
-
Microscopic description of anisotropic low-density dipolar Bose gases in two dimensions
Authors:
A. Macia,
F. Mazzanti,
J. Boronat,
R. E. Zillich
Abstract:
A microscopic description of the zero energy two-body ground state and many-body static properties of anisotropic homogeneous gases of bosonic dipoles in two dimensions at low densities is presented and discussed. By changing the polarization angle with respect to the plane, we study the impact of the anisotropy, present in the dipole--dipole interaction, on the energy per particle, comparing the…
▽ More
A microscopic description of the zero energy two-body ground state and many-body static properties of anisotropic homogeneous gases of bosonic dipoles in two dimensions at low densities is presented and discussed. By changing the polarization angle with respect to the plane, we study the impact of the anisotropy, present in the dipole--dipole interaction, on the energy per particle, comparing the results with mean field predictions. We restrict the analysis to the regime where the interaction is always repulsive, although the strength of the repulsion depends on the orientation with respect to the polarization field. We present a series expansion of the solution of the zero energy two-body problem which allows us to find the scattering length of the interaction and to build a suitable Jastrow factor that we use as a trial wave function for both a variational and diffusion Monte Carlo simulation of the infinite system. We find that the anisotropy has an almost negligible impact on the ground state properties of the many-body system in the universal regime where the scattering length governs the physics of the system. We also show that scaling in the gas parameter persists in the dipolar case up to values where other isotropic interactions with the same scattering length yield different predictions.
△ Less
Submitted 30 May, 2011;
originally announced May 2011.
-
Extrapolated High-Order Propagators for Path Integral Monte Carlo Simulations
Authors:
Robert E. Zillich,
Johannes M. Mayrhofer,
Siu A. Chin
Abstract:
We present a new class of high-order imaginary time propagators for path-integral Monte Carlo simulations by subtracting lower order propagators. By requiring all terms of the extrapolated propagator be sampled uniformly, the subtraction only affects the potential part of the path integral. The negligible violation of positivity of the resulting path integral at small time steps has no discernab…
▽ More
We present a new class of high-order imaginary time propagators for path-integral Monte Carlo simulations by subtracting lower order propagators. By requiring all terms of the extrapolated propagator be sampled uniformly, the subtraction only affects the potential part of the path integral. The negligible violation of positivity of the resulting path integral at small time steps has no discernable affect on the accuracy of our method. Thus in principle arbitrarily high order algorithms can be devised for path-integral Monte Carlo simulations. We verify this claim is by showing that fourth, sixth, and eighth order convergence can indeed be achieved in solving for the ground state of strongly interacting quantum many-body systems such as bulk liquid $^4$He.
△ Less
Submitted 20 July, 2009;
originally announced July 2009.
-
Roton-Rotation Coupling of Acetylene in $^4$He
Authors:
R. E. Zillich,
K. B. Whaley
Abstract:
Rotational absorption spectra of acetylene in superfluid $^4$He are calculated using a combined correlated basis function-diffusion Monte Carlo method (CBF/DMC) as well as a numerically exact path integral correlation function approach (PICF). Both approaches predict anomalously large distortion constants in addition to reduced rotational constants, and the higher rotational excitations are seen…
▽ More
Rotational absorption spectra of acetylene in superfluid $^4$He are calculated using a combined correlated basis function-diffusion Monte Carlo method (CBF/DMC) as well as a numerically exact path integral correlation function approach (PICF). Both approaches predict anomalously large distortion constants in addition to reduced rotational constants, and the higher rotational excitations are seen to be accompanied by an additional peak and absorption band. The CBF analysis shows that these phenomena are due to strong coupling of the higher rotational states with the roton and maxon excitations of $^4$He, and the assocated divergence of the helium density of states in this region.
△ Less
Submitted 15 April, 2004;
originally announced April 2004.
-
Comparison of Rotational Energies and Rigidity of OCS-paraH_2 and OCS-4He complexes
Authors:
R. E. Zillich,
K. B. Whaley
Abstract:
We analyze the nature of the rotational energy level structure of the OCS-He and OCS-H_2 complexes with a comparison of exact calculations to several differentdynamical approximations. We compare with the clamped coordinate quasiadiabatic approximation that introduces an effective potential for each asymmetric rotor level, with an effective rotation Hamiltonian constructed from ground state aver…
▽ More
We analyze the nature of the rotational energy level structure of the OCS-He and OCS-H_2 complexes with a comparison of exact calculations to several differentdynamical approximations. We compare with the clamped coordinate quasiadiabatic approximation that introduces an effective potential for each asymmetric rotor level, with an effective rotation Hamiltonian constructed from ground state averages of the inverse of the inertial matrix, and investigate the usefulness of the Eckart condition to decouple rotations and vibrations of these weakly bound complexes between linear OCS and 4He or H_2. Comparison with exact results allows an assessment of the accuracies of the different approximate methods and indicates which approaches are suitable for larger clusters of OCS with 4He and with H_2. We find the OCS-H_2 complex is considerably more rigid than the OCS-4He complex, suggesting that semi-rigid models are useful for analysis of larger clusters of H_2 with OCS.
△ Less
Submitted 17 September, 2003;
originally announced September 2003.
-
OCS in small para-hydrogen clusters: energetics and structure with N=1-8 complexed hydrogen molecules
Authors:
F. Paesani,
R. E. Zillich,
K. B. Whaley
Abstract:
We determine the structure and energetics of complexes of the linear OCS molecule with small numbers of para-hydrogen molecules, N=1-8, using zero temperature quantum Monte Carlo methods. Ground state calculations are carried out with importance-sampled rigid body diffusion Monte Carlo (IS-RBDMC) and excited state calculations with the projection operator imaginary time spectral evolution (POITS…
▽ More
We determine the structure and energetics of complexes of the linear OCS molecule with small numbers of para-hydrogen molecules, N=1-8, using zero temperature quantum Monte Carlo methods. Ground state calculations are carried out with importance-sampled rigid body diffusion Monte Carlo (IS-RBDMC) and excited state calculations with the projection operator imaginary time spectral evolution (POITSE) methodology. The ground states are found to be highly structured, with a gradual build up of two axial rings as N increases to 8. Analysis of the azimuthal density correlations around the OCS molecule shows that these rings are quite delocalized for small N values, but become strongly localized for N \geq 5 . Excited state calculations are made for a range of total cluster angular momentum values and the rotational energy levels fitted to obtain effective rotational and distortion constants of the complexed OCS molecule as a function of cluster size N. Detailed analysis of these spectroscopic constants indicates that the complexes of OCS with para-hydrogen have an unusually rich variation in dynamical behavior, with sizes N=1-2 showing near rigid behavior, sizes N=3-4 showing extremely floppy behavior, and the larger sizes N=5-8 showing more rigid behavior again. The large values of the distortion constant D obtained for N=3-4 are rationalized in terms of the coupling between the OCS rotations and the "breathing" mode of the first, partially filled ring of para-hydrogen molecules.
△ Less
Submitted 16 September, 2003;
originally announced September 2003.
-
Quantum Rotation of HCN and DCN in Helium-4
Authors:
R. E. Zillich,
K. B. Whaley
Abstract:
We present calculations of rotational absorption spectra of the molecules HCN and DCN in superfluid helium-4, using a combination of the Diffusion Monte Carlo method for ground state properties and an analytic many-body method (Correlated Basis Function theory) for the excited states. Our results agree with the experimentally determined effective moment of inertia which has been obtained from th…
▽ More
We present calculations of rotational absorption spectra of the molecules HCN and DCN in superfluid helium-4, using a combination of the Diffusion Monte Carlo method for ground state properties and an analytic many-body method (Correlated Basis Function theory) for the excited states. Our results agree with the experimentally determined effective moment of inertia which has been obtained from the $J=0\to 1$ spectral transition. The correlated basis function analysis shows that, unlike heavy rotors such as OCS, the J=2 and higher rotational excitations of HCN and DCN have high enough energy to strongly couple to rotons, leading to large shifts of the lines and accordingly to anomalous large spectroscopic distortion constants, to the emergence of roton-maxon bands, and to secondary peaks in the absorption spectra for J=2 and J=3.
△ Less
Submitted 15 December, 2003; v1 submitted 10 September, 2003;
originally announced September 2003.