Source author record

S. Pilati

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

16 published item(s)

preprint2022arXiv

Path-integral Monte Carlo worm algorithm for Bose systems with periodic boundary conditions

We provide a detailed description of the path-integral Monte Carlo worm algorithm used to exactly calculate the thermodynamics of Bose systems in the canonical ensemble. The algorithm is fully consistent with periodic boundary conditions, that are applied to simulate homogeneous phases of bulk systems, and it does not require any limitation in the length of the Monte Carlo moves realizing the sampling of the probability distribution function in the space of path configurations. The result is achieved adopting a representation of the path coordinates where only the initial point of each path is inside the simulation box, the remaining ones being free to span the entire space. Detailed balance can thereby be ensured for any update of the path configurations without the ambiguity of the selection of the periodic image of the particles involved. We benchmark the algorithm using the non-interacting Bose gas model for which exact results for the partition function at finite number of particles can be derived. Convergence issues and the approach to the thermodynamic limit are also addressed for interacting systems of hard spheres in the regime of high density.

preprint2020arXiv

Boosting Monte Carlo simulations of spin glasses using autoregressive neural networks

The autoregressive neural networks are emerging as a powerful computational tool to solve relevant problems in classical and quantum mechanics. One of their appealing functionalities is that, after they have learned a probability distribution from a dataset, they allow exact and efficient sampling of typical system configurations. Here we employ a neural autoregressive distribution estimator (NADE) to boost Markov chain Monte Carlo (MCMC) simulations of a paradigmatic classical model of spin-glass theory, namely the two-dimensional Edwards-Anderson Hamiltonian. We show that a NADE can be trained to accurately mimic the Boltzmann distribution using unsupervised learning from system configurations generated using standard MCMC algorithms. The trained NADE is then employed as smart proposal distribution for the Metropolis-Hastings algorithm. This allows us to perform efficient MCMC simulations, which provide unbiased results even if the expectation value corresponding to the probability distribution learned by the NADE is not exact. Notably, we implement a sequential tempering procedure, whereby a NADE trained at a higher temperature is iteratively employed as proposal distribution in a MCMC simulation run at a slightly lower temperature. This allows one to efficiently simulate the spin-glass model even in the low-temperature regime, avoiding the divergent correlation times that plague MCMC simulations driven by local-update algorithms. Furthermore, we show that the NADE-driven simulations quickly sample ground-state configurations, paving the way to their future utilization to tackle binary optimization problems.

preprint2020arXiv

Scalable neural networks for the efficient learning of disordered quantum systems

Supervised machine learning is emerging as a powerful computational tool to predict the properties of complex quantum systems at a limited computational cost. In this article, we quantify how accurately deep neural networks can learn the properties of disordered quantum systems as a function of the system size. We implement a scalable convolutional network that can address arbitrary system sizes. This network is compared with a recently introduced extensive convolutional architecture [K. Mills et al., Chem. Sci. 10, 4129 (2019)] and with conventional dense networks with all-to-all connectivity. The networks are trained to predict the exact ground-state energies of various disordered systems, namely a continuous-space single-particle Hamiltonian for cold-atoms in speckle disorder, and different setups of a quantum Ising chain with random couplings, including one with only short-range interactions and one augmented with a long-range term. In all testbeds we consider, the scalable network retains high accuracy as the system size increases. Furthermore, we demonstrate that the network scalability enables a transfer-learning protocol, whereby a pre-training performed on small systems drastically accelerates the learning of large-system properties, allowing reaching high accuracy with small training sets. In fact, with the scalable network one can even extrapolate to sizes larger than those included in the training set, accurately reproducing the results of state-of-the-art quantum Monte Carlo simulations.

preprint2020arXiv

Simulating disordered quantum systems via dense and sparse restricted Boltzmann machines

In recent years, generative artificial neural networks based on restricted Boltzmann machines (RBMs) have been successfully employed as accurate and flexible variational wave functions for clean quantum many-body systems. In this article we explore their use in simulations of disordered quantum spin models. The standard dense RBM with all-to-all inter-layer connectivity is not particularly appropriate for large disordered systems, since in such systems one cannot exploit translational invariance to reduce the amount of parameters to be optimized. To circumvent this problem, we implement sparse RBMs, whereby the visible spins are connected only to a subset of local hidden neurons, thus reducing the amount of parameters. We assess the performance of sparse RBMs as a function of the range of the allowed connections, and compare it with the one of dense RBMs. Benchmark results are provided for two sign-problem free Hamiltonians, namely pure and random quantum Ising chains. The RBM ansatzes are trained using the unsupervised learning scheme based on projective quantum Monte Carlo (PQMC) algorithms. We find that the sparse connectivity facilitates the training process and allows sparse RBMs to outperform the dense counterparts. Furthermore, the use of sparse RBMs as guiding functions for PQMC simulations allows us to perform PQMC simulations at a reduced computational cost, avoiding possible biases due to finite random-walker populations. We obtain unbiased predictions for the ground-state energies and the magnetization profiles with fixed boundary conditions, at the ferromagnetic quantum critical point. The magnetization profiles agree with the Fisher-de Gennes scaling relation for conformally invariant systems, including the scaling dimension predicted by the renormalization-group analysis.

preprint2016arXiv

Ferromagnetism in a Repulsive Atomic Fermi Gas with Correlated Disorder

We investigate the zero-temperature ferromagnetic behavior of a two-component repulsive Fermi gas in the presence of a correlated random field that represents an optical speckle pattern. The density is tuned so that the (noninteracting) Fermi energy is close to the mobility edge of the Anderson localization transition. We employ quantum Monte Carlo simulations to determine various ground-state properties, including the equation of state, the magnetic susceptibility, and the energy of an impurity immersed in a polarized Fermi gas (repulsive polaron). In the weakly interacting limit, the magnetic susceptibility is found to be suppressed by disorder. However, it rapidly increases with the interaction strength, and it diverges at a much weaker interaction strength compared to the clean gas. Both the transition from the paramagnetic phase to the partially ferromagnetic phase, and the one from the partially to the fully ferromagnetic phase are strongly favored by disorder, indicating a case of order induced by disorder.

preprint2015arXiv

Simulated quantum annealing of double-well and multi-well potentials

We analyze the performance of quantum annealing as a heuristic optimization method to find the absolute minimum of various continuous models, including landscapes with only two wells and also models with many competing minima and with disorder. The simulations performed using a projective quantum Monte Carlo (QMC) algorithm are compared with those based on the finite-temperature path-integral QMC technique and with classical annealing. We show that the projective QMC algorithm is more efficient than the finite-temperature QMC technique, and that both are inferior to classical annealing if this is performed with appropriate long-range moves. However, as the difficulty of the optimization problem increases, classical annealing looses efficiency, while the projective QMC algorithm keeps stable performance and is finally the most effective optimization tool. We discuss the implications of our results for the outstanding problem of testing the efficiency of adiabatic quantum computers using stochastic simulations performed on classical computers.

preprint2014arXiv

Ferromagnetism of a Repulsive Atomic Fermi Gas in an Optical Lattice: a Quantum Monte Carlo Study

Using continuous-space quantum Monte Carlo methods we investigate the zero-temperature ferromagnetic behavior of a two-component repulsive Fermi gas under the influence of periodic potentials that describe the effect of a simple-cubic optical lattice. Simulations are performed with balanced and with imbalanced components, including the case of a single impurity immersed in a polarized Fermi sea (repulsive polaron). For an intermediate density below half filling, we locate the transitions between the paramagnetic, and the partially and the fully ferromagnetic phases. As the intensity of the optical lattice increases, the ferromagnetic instability takes place at weaker interactions, indicating a possible route to observe ferromagnetism in experiments performed with ultracold atoms. We compare our findings with previous predictions based on the standard computational method used in material science, namely density functional theory, and with results based on tight-binding models.

preprint2013arXiv

Critical Temperature of Interacting Bose Gases in Periodic Potentials

The superfluid transition of a repulsive Bose gas in the presence of a sinusoidal potential which represents a simple-cubic optical lattice is investigate using quantum Monte Carlo simulations. At the average filling of one particle per well the critical temperature has a nonmonotonic dependence on the interaction strength, with an initial sharp increase and a rapid suppression at strong interactions in the vicinity of the Mott transition. In an optical lattice the positive shift of the transition is strongly enhanced compared to the homogenous gas. By varying the lattice filling we find a crossover from a regime where the optical lattice has the dominant effect to a regime where interactions dominate and the presence of the lattice potential becomes almost irrelevant.

preprint2012arXiv

Bosonic superfluid-insulator transition in continuous space

We investigate the zero-temperature phase diagram of interacting Bose gases in the presence of a simple cubic optical lattice, going beyond the regime where the mapping to the single-band Bose-Hubbard model is reliable. Our computational approach is a new hybrid quantum Monte Carlo method which combines algorithms used to simulate homogeneous quantum fluids in continuous space with those used for discrete lattice models of strongly correlated systems. We determine the critical interaction strength and optical lattice intensity where the superfluid-to-insulator transition takes place, considering also the regime of shallow optical lattices and strong inter-atomic interactions. The implications of our findings for the super-solid state of matter are discussed.

preprint2010arXiv

Dilute Bose gas with correlated disorder: A Path Integral Monte Carlo study

We investigate the thermodynamic properties of a dilute Bose gas in a correlated random potential using exact path integral Monte Carlo methods. The study is carried out in continuous space and disorder is produced in the simulations by a 3D speckle pattern with tunable intensity and correlation length. We calculate the shift of the superfluid transition temperature due to disorder and we highlight the role of quantum localization by comparing the critical chemical potential with the classical percolation threshold. The equation of state of the gas is determined in the regime of strong disorder, where superfluidity is suppressed and the normal phase exists down to very low temperatures. We find a $T^2$ dependence of the energy in agreement with the expected behavior in the Bose glass phase. We also discuss the major role played by the disorder correlation length and we make contact with a Hartree-Fock mean-field approach that holds valid if the correlation length is very large. The density profiles are analyzed as a function of temperature and interaction strength. Effects of localization and the depletion of the order parameter are emphasized in the comparison between local condensate and total density. At very low temperature we find that the energy and the particle distribution of the gas are very well described by the T=0 Gross-Pitaevskii theory even in the regime of very strong disorder.

preprint2010arXiv

Itinerant ferromagnetism of a repulsive atomic Fermi gas: a quantum Monte Carlo study

We investigate the phase diagram of a two-component repulsive Fermi gas at T=0 by means of quantum Monte Carlo simulations. For a given value of the positive s-wave scattering length, both purely repulsive and purely attractive model potentials are considered in order to analyze the limits of the universal regime where the details of interatomic forces can be neglected. The equation of state of both balanced and unbalanced systems is calculated as a function of the interaction strength and the critical density for the onset of ferromagnetism is determined. The energy per particle of the strongly polarized gas is calculated and parametrized in terms of the physical properties of repulsive polarons, which are relevant for the stability of the fully magnetized ferromagnetic state. Finally, we analyze the phase diagram in the polarization/interaction plane under the assumption that only phases with homogeneous magnetization can be produced.

preprint2009arXiv

Beliaev technique for a weakly interacting Bose gas

Aiming for simplicity of explicit equations and at the same time controllable accuracy of the theory we present results for all thermodynamic quantities and correlation functions for the weakly interacting Bose gas at short-to-intermediate distances obtained within an improved version of Beliaev's diagrammatic technique. With a small symmetry breaking term Beliaev's diagrammatic technique becomes regular in the infrared limit. Up to higher-order terms (for which we present order-of-magnitude estimates), the partition function and entropy of the system formally correspond to those of a non-interacting bosonic (pseudo-)Hamiltonian with a temperature dependent Bogoliubov-type dispersion relation. Away from the fluctuation region, this approach provides the most accurate--in fact, the best possible within the Bogoliubov-type pseudo-Hamiltonian framework--description of the system with controlled accuracy. It produces accurate answers for the off-diagonal correlation functions up to distances where the behaviour of correlators is controlled by generic hydrodynamic relations, and thus can be accurately extrapolated to arbitrarily large scales. In the fluctuation region, the non-perturbative contributions are given by universal (for all weakly interacting U(1) systems) constants and scaling functions, which can be obtained separately--by simulating classical U(1) models--and then used to extend the description of the weakly interacting Bose gas to the fluctuation region. The theory works in all spatial dimensions and we explicitly check its validity against first-principle Monte Carlo simulations for various thermodynamic properties and the single-particle density matrix.

preprint2009arXiv

Superfluid transition in a Bose gas with correlated disorder

The superfluid transition of a three-dimensional gas of hard-sphere bosons in a disordered medium is studied using quantum Monte Carlo methods. Simulations are performed in continuous space both in the canonical and in the grand-canonical ensemble. At fixed density we calculate the shift of the transition temperature as a function of the disorder strength, while at fixed temperature we determine both the critical chemical potential and the critical density separating normal and superfluid phases. In the regime of strong disorder the normal phase extends up to large values of the degeneracy parameter and the critical chemical potential exhibits a linear dependence in the intensity of the random potential. The role of interactions and disorder correlations is also discussed.

preprint2008arXiv

Critical temperature of interacting Bose gases in two and three dimensions

We calculate the superfluid transition temperature of homogeneous interacting Bose gases in three and two spatial dimensions using large-scale Path Integral Monte Carlo simulations (with up to $N=10^5$ particles). In 3D we investigate the limits of the universal critical behavior in terms of the scattering length alone by using different models for the interatomic potential. We find that this type of universality sets in at small values of the gas parameter $na^3 \lesssim 10^{-4}$. This value is different from the estimate $na^3 \lesssim 10^{-6}$ for the validity of the asymptotic expansion in the limit of vanishing $na^3$. In 2D we study the Berezinskii-Kosterlitz-Thouless transition of a gas with hard-core interactions. For this system we find good agreement with the classical lattice $|ψ|^4$ model up to very large densities. We also explain the origin of the existing discrepancy between previous studies of the same problem.

preprint2007arXiv

Phase separation in a polarized Fermi gas at zero temperature

We investigate the phase diagram of asymmetric two-component Fermi gases at zero temperature as a function of polarization and interaction strength. The equations of state of the uniform superfluid and normal phase are determined using quantum Monte Carlo simulations. We find three different mixed states, where the superfluid and the normal phase coexist in equilibrium, corresponding to phase separation between: (a) the polarized superfluid and the fully polarized normal gas, (b) the polarized superfluid and the partially polarized normal gas and (c) the unpolarized superfluid and the partially polarized normal gas.

preprint2006arXiv

Equation of state of an interacting Bose gas at finite temperature: a Path Integral Monte Carlo study

By using exact Path Integral Monte Carlo methods we calculate the equation of state of an interacting Bose gas as a function of temperature both below and above the superfluid transition. The universal character of the equation of state for dilute systems and low temperatures is investigated by modeling the interatomic interactions using different repulsive potentials corresponding to the same s-wave scattering length. The results obtained for the energy and the pressure are compared to the virial expansion for temperatures larger than the critical temperature. At very low temperatures we find agreement with the ground-state energy calculated using the diffusion Monte Carlo method.