Source author record

Robert E. Zillich

Robert E. Zillich appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.quant-gas cond-mat.stat-mech physics.atm-clus physics.chem-ph physics.comp-ph

Catalog footprint

What is connected

7works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Time-dependent variational Monte Carlo study of the dynamic response of bosons in an optical lattice

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.

preprint2020arXiv

Rotational coherence spectroscopy of molecules in helium nanodroplets: Reconciling the time and the frequency domains

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.

preprint2020arXiv

Rotons and Bose condensation in Rydberg-dressed Bose Gases

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.

preprint2015arXiv

Dipolar Bilayer with Antiparallel Polarization -- a Self-Bound Liquid

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.

preprint2013arXiv

Stability and excitations of a bilayer of strongly correlated dipolar Bosons

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.

preprint2012arXiv

Excitations and Stripe Phase Formation in a 2D Dipolar Bose Gas with Tilted Polarization

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.

preprint2009arXiv

Extrapolated High-Order Propagators for Path Integral Monte Carlo Simulations

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.