We have extracted information about real time dynamics of H4e systems from noisy imaginary-time correlation functions f($\tau${}) computed via quantum Monte Carlo (QMC): production and falsification of model spectral functions s($\omega${}) are obtained via a survival-to-compatibility with f($\tau${}) evolutionary process, based on genetic algorithms. Statistical uncertainty in f($\tau${}) is promoted to be an asset via a sampling of equivalent f($\tau${}) within the noise, which give rise to independent evolutionary processes. In the case of pure superfluid H4e we have recovered from exact QMC simulations sharp quasiparticle excitations with spectral functions displaying also the multiphonon branch. As further applications, we have studied the impuriton branch of one H3e atom in liquid H4e and the vacancy-wave excitations in hcp solid H4e finding an unexpected rotonlike feature.

We determine the spin and charge orders in the ground state of the doped two-dimensional (2D) Hubbard model in its simplest form, namely with only nearest-neighbor hopping and on-site repulsion. At half-filling, the ground state is known to be an antiferromagnetic Mott insulator. Doping Mott insulators is believed to be relevant to the superconductivity observed in cuprates. A variety of candidates have been proposed for the ground state of the doped 2D Hubbard model. A recent work employing a combination of several state-of-the-art numerical many-body methods established the stripe order as the ground state near $1/8$ doping at strong interactions. In this paper, we apply one of these methods, the cutting-edge constrained-path auxiliary field quantum Monte Carlo (AFQMC) method with self-consistently optimized gauge constraints, to systematically study the model as a function of doping and interaction strength. With careful finite size scaling based on large-scale computations, we map out the ground state phase diagram in terms of its spin and charge order. We find that modulated antiferromagnetic order persists from near half-filling to about $1/5$ doping. At lower interaction strengths or larger doping, these ordered states are best described as spin-density waves, with essentially delocalized holes and modest oscillations in charge correlations. When the charge correlations are stronger (large interaction or small doping), they are best described as stripe states, with the holes more localized near the node in the antiferromagnetic spin order. In both cases, we find that the wavelength in the charge correlations is consistent with so-called filled stripes.

We address the calculation of dynamical correlation functions for many fermion systems at zero temperature, using the auxiliary-field quantum Monte Carlo method. The two-dimensional Hubbard hamiltonian is used as a model system. Although most of the calculations performed here are for cases where the sign problem is absent, the discussions are kept general for applications to physical problems when the sign problem does arise. We study the use of twisted boundary conditions to improve the extrapolation of the results to the thermodynamic limit. A strategy is proposed to drastically reduce finite size effects relying on a minimization among the twist angles. This approach is demonstrated by computing the charge gap at half filling. We obtain accurate results showing the scaling of the gap with the interaction strength $U$ in two dimensions, connecting to the scaling of the unrestricted Hartree-Fock method at small $U$ and Bethe ansatz exact result in one dimension at large $U$. An alternative algorithm is then proposed to compute dynamical Green functions and correlation functions which explicitly varies the number of particles during the random walks in the manifold of Slater determinants. In dilute systems, such as ultracold Fermi gases, this algorithm enables calculations with much more favorable complexity, with computational cost proportional to basis size or the number of lattice sites.

We compute the zero-temperature dynamical structure factor of one-dimensional liquid ${}^4\mathrm{He}$ by means of state-of-the-art quantum Monte Carlo and analytic continuation techniques. By increasing the density, the dynamical structure factor reveals a transition from a highly compressible critical liquid to a quasisolid regime. In the low-energy limit, the dynamical structure factor can be described by the quantum hydrodynamic Luttinger-liquid theory, with a Luttinger parameter spanning all possible values by increasing the density. At higher energies, our approach provides quantitative results beyond the Luttinger-liquid theory. In particular, as the density increases, the interplay between dimensionality and interaction makes the dynamical structure factor manifest a pseudo-particle-hole continuum typical of fermionic systems. At the low-energy boundary of such a region and moderate densities, we find consistency, within statistical uncertainties, with predictions of a power-law structure by the recently developed nonlinear Luttinger-liquid theory. In the quasisolid regime, we observe a novel behavior at intermediate momenta, which can be described by new analytical relations that we derive for the hard-rods model.

We evaluate imaginary time density-density correlation functions for two-dimensional homogeneous electron gases of up to 42 particles in the continuum using the phaseless auxiliary field quantum Monte Carlo method. We use periodic boundary conditions and up to 300 plane waves as basis set elements. We show that such methodology, once equipped with suitable numerical stabilization techniques necessary to deal with exponentials, products, and inversions of large matrices, gives access to the calculation of imaginary time correlation functions for medium-sized systems. We discuss the numerical stabilization techniques and the computational complexity of the methodology and we present the limitations related to the size of the systems on a quantitative basis. We perform the inverse Laplace transform of the obtained density-density correlation functions, assessing the ability of the phaseless auxiliary field quantum Monte Carlo method to evaluate dynamical properties of medium-sized homogeneous fermion systems.

We develop a self-consistent variant of the constrained path quantum Monte Carlo approach which ensures its independence of the trial wave function, and apply the method to compute ground-state correlations in the two-dimensional SU(3) Fermi-Hubbard Hamiltonian at $\frac{1}{3}$ filling, modeling fermions with three possible spin flavors moving on a square lattice with an average of one particle per site. We provide clear evidence of a quantum critical point separating a nonmagnetic uniform metallic phase from a regime where long-range “spin” order is present. This discovery of multiple successive transitions to magnetic states with regular, long-range alternation of the different flavors, whose symmetry changes as the interaction strength increases, significantly extends previous work in the Heisenberg limit to itinerant fermions. In addition to the rich quantum magnetism, this important physical system allows one to study integer filling and the associated Mott transition disentangled from nesting, in contrast to the usual SU(2) model, while preserving the square-lattice geometry. Our results also provide a significant step towards the interpretation of present and future experiments on fermionic alkaline-earth atoms, and other realizations of $\mathrm{SU}\left(N\right)$ physics.