Showing 1–50 of 52 results for author: Galli, E

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

Abstract: PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observati… ▽ More PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases. △ Less

Submitted 8 June, 2024; originally announced June 2024.

Euclid. II. The VIS Instrument

Authors: Euclid Collaboration, M. Cropper, A. Al-Bahlawan, J. Amiaux, S. Awan, R. Azzollini, K. Benson, M. Berthe, J. Boucher, E. Bozzo, C. Brockley-Blatt, G. P. Candini, C. Cara, R. A. Chaudery, R. E. Cole, P. Danto, J. Denniston, A. M. Di Giorgio, B. Dryer, J. Endicott, J. -P. Dubois, M. Farina, E. Galli, L. Genolet, J. P. D. Gow , et al. (403 additional authors not shown)

Abstract: This paper presents the specification, design, and development of the Visible Camera (VIS) on the ESA Euclid mission. VIS is a large optical-band imager with a field of view of 0.54 deg^2 sampled at 0.1" with an array of 609 Megapixels and spatial resolution of 0.18". It will be used to survey approximately 14,000 deg^2 of extragalactic sky to measure the distortion of galaxies in the redshift ran… ▽ More This paper presents the specification, design, and development of the Visible Camera (VIS) on the ESA Euclid mission. VIS is a large optical-band imager with a field of view of 0.54 deg^2 sampled at 0.1" with an array of 609 Megapixels and spatial resolution of 0.18". It will be used to survey approximately 14,000 deg^2 of extragalactic sky to measure the distortion of galaxies in the redshift range z=0.1-1.5 resulting from weak gravitational lensing, one of the two principal cosmology probes of Euclid. With photometric redshifts, the distribution of dark matter can be mapped in three dimensions, and, from how this has changed with look-back time, the nature of dark energy and theories of gravity can be constrained. The entire VIS focal plane will be transmitted to provide the largest images of the Universe from space to date, reaching m_AB>24.5 with S/N >10 in a single broad I_E~(r+i+z) band over a six year survey. The particularly challenging aspects of the instrument are the control and calibration of observational biases, which lead to stringent performance requirements and calibration regimes. With its combination of spatial resolution, calibration knowledge, depth, and area covering most of the extra-Galactic sky, VIS will also provide a legacy data set for many other fields. This paper discusses the rationale behind the VIS concept and describes the instrument design and development before reporting the pre-launch performance derived from ground calibrations and brief results from the in-orbit commissioning. VIS should reach fainter than m_AB=25 with S/N>10 for galaxies of full-width half-maximum of 0.3" in a 1.3" diameter aperture over the Wide Survey, and m_AB>26.4 for a Deep Survey that will cover more than 50 deg^2. The paper also describes how VIS works with the other Euclid components of survey, telescope, and science data processing to extract the cosmological information. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Spectrum of density, spin and pairing fluctuations of an attractive two-dimensional Fermi gas

Authors: Christian Apostoli, Patrick Kelly, Annette Lopez, Kaelyn Dauer, Gianluca Bertaina, Davide Emilio Galli, Ettore Vitali

Abstract: We leverage random phase approximation and unbiased auxiliary-field quantum Monte Carlo methods to compute dynamical correlations for a dilute homogeneous two-dimensional attractive Fermi gas. Our main purpose is to quantitatively study the collective excitations of the system to generate robust benchmark results and to shed light into fermionic superfluidity in the strongly correlated regime. In… ▽ More We leverage random phase approximation and unbiased auxiliary-field quantum Monte Carlo methods to compute dynamical correlations for a dilute homogeneous two-dimensional attractive Fermi gas. Our main purpose is to quantitatively study the collective excitations of the system to generate robust benchmark results and to shed light into fermionic superfluidity in the strongly correlated regime. In particular we are motivated by a recent paper suggesting that the Higgs mode can be detected in the spectrum of spin fluctuations. Despite the fact that we are somewhat limited by finite-size effects, our study pinpoints indeed a clear peak in the spin channel at low momentum, but a detailed analysis suggests that such a peak, though certainly interesting, does not correspond to the Higgs mode. We propose a different explanation for the shape of the spin structure factor. On the other hand, our results clearly show that the Higgs mode can be detected in the density channel at very small wave vectors, although very good resolution is necessary. △ Less

Submitted 5 July, 2024; v1 submitted 20 March, 2024; originally announced March 2024.

Comments: 10 pages, 5 figures

Quantum Monte Carlo and perturbative study of two-dimensional Bose-Fermi mixtures

Authors: Jacopo D'Alberto, Lorenzo Cardarelli, Davide Emilio Galli, Gianluca Bertaina, Pierbiagio Pieri

Abstract: We derive analytically the leading beyond-mean field contributions to the zero-temperature equation of state and to the fermionic quasi-particle residue and effective mass of a dilute Bose-Fermi mixture in two dimensions. In the repulsive case, we perform quantum Monte Carlo simulations for two representative bosonic concentrations and equal masses, extending a method for correcting finite-size ef… ▽ More We derive analytically the leading beyond-mean field contributions to the zero-temperature equation of state and to the fermionic quasi-particle residue and effective mass of a dilute Bose-Fermi mixture in two dimensions. In the repulsive case, we perform quantum Monte Carlo simulations for two representative bosonic concentrations and equal masses, extending a method for correcting finite-size effects in fermionic gases to Bose-Fermi mixtures. We find good agreement between analytic expressions and numerical results for weak interactions, while significant discrepancies appear in the regime close to mechanical instability, above which we provide evidence of phase separation of the bosonic component. △ Less

Submitted 8 May, 2024; v1 submitted 21 February, 2024; originally announced February 2024.

Comments: 16 pages, 9 figures. Added appendix and changed title

Journal ref: Phys. Rev. A 109, 053302 (2024)

N-electron valence perturbation theory with reference wavefunctions from quantum computing: application to the relative stability of hydroxide anion and hydroxyl radical

Authors: Alessandro Tammaro, Davide E. Galli, Julia E. Rice, Mario Motta

Abstract: Quantum simulations of the hydroxide anion and hydroxyl radical are reported, employing variational quantum algorithms for near-term quantum devices. The energy of each species is calculated along the dissociation curve, to obtain information about the stability of the molecular species being investigated. It is shown that simulations restricted to valence spaces incorrectly predict the hydroxyl r… ▽ More Quantum simulations of the hydroxide anion and hydroxyl radical are reported, employing variational quantum algorithms for near-term quantum devices. The energy of each species is calculated along the dissociation curve, to obtain information about the stability of the molecular species being investigated. It is shown that simulations restricted to valence spaces incorrectly predict the hydroxyl radical to be more stable than the hydroxide anion. Inclusion of dynamical electron correlation from non-valence orbitals is demonstrated, through the integration of the variational quantum eigensolver and quantum subspace expansion methods in the workflow of N-electron valence perturbation theory, and shown to correctly predict the hydroxide anion to be more stable than the hydroxyl radical, provided that basis sets with diffuse orbitals are also employed. Finally, we calculate the electron affinity of the hydroxyl radical using an aug-cc-pVQZ basis on IBM's quantum devices. △ Less

Submitted 16 October, 2022; v1 submitted 25 February, 2022; originally announced February 2022.

Comments: 12 pages, 6 figures

Journal ref: J. Phys. Chem. A 127, 817-827 (2023)

Quantum circuits for the preparation of spin eigenfunctions on quantum computers

Authors: Alessandro Carbone, Davide Emilio Galli, Mario Motta, Barbara Jones

Abstract: The application of quantum algorithms to the study of many-particle quantum systems requires the ability to prepare wavefunctions that are relevant in the behavior of the system under study. Hamiltonian symmetries are an important instrument, to classify relevant many-particle wavefunctions, and to improve the efficiency of numerical simulations. In this work, quantum circuits for the exact and ap… ▽ More The application of quantum algorithms to the study of many-particle quantum systems requires the ability to prepare wavefunctions that are relevant in the behavior of the system under study. Hamiltonian symmetries are an important instrument, to classify relevant many-particle wavefunctions, and to improve the efficiency of numerical simulations. In this work, quantum circuits for the exact and approximate preparation of total spin eigenfunctions on quantum computers are presented. Two different strategies are discussed and compared: exact recursive construction of total spin eigenfunctions based on the addition theorem of angular momentum, and heuristic approximation of total spin eigenfunctions based on the variational optimization of a suitable cost function. The construction of these quantum circuits is illustrated in detail, and the preparation of total spin eigenfunctions is demonstrated on IBM quantum devices, focusing on 3- and 5-spin systems on graphs with triangle connectivity. △ Less

Submitted 18 February, 2022; originally announced February 2022.

Comments: 15 pages, 14 figures, 6 tables

Journal ref: Symmetry 2022, 14(3), 624

Multi-class quantum classifiers with tensor network circuits for quantum phase recognition

Authors: Marco Lazzarin, Davide Emilio Galli, Enrico Prati

Abstract: Hybrid quantum-classical algorithms based on variational circuits are a promising approach to quantum machine learning problems for near-term devices, but the selection of the variational ansatz is an open issue. Recently, tensor network-inspired circuits have been proposed as a natural choice for such ansatz. Their employment on binary classification tasks provided encouraging results. However, t… ▽ More Hybrid quantum-classical algorithms based on variational circuits are a promising approach to quantum machine learning problems for near-term devices, but the selection of the variational ansatz is an open issue. Recently, tensor network-inspired circuits have been proposed as a natural choice for such ansatz. Their employment on binary classification tasks provided encouraging results. However, their effectiveness on more difficult tasks is still unknown. Here, we present numerical experiments on multi-class classifiers based on tree tensor network and multiscale entanglement renormalization ansatz circuits. We conducted experiments on image classification with the MNIST dataset and on quantum phase recognition with the XXZ model by Cirq and TensorFlow Quantum. In the former case, we reduced the number of classes to four to match the aimed output based on 2 qubits. The quantum data of the XXZ model consist of three classes of ground states prepared by a checkerboard circuit used for the ansatz of the variational quantum eigensolver, corresponding to three distinct quantum phases. Test accuracy turned out to be 59%-93% and 82%-96% respectively, depending on the model architecture and on the type of preprocessing. △ Less

Submitted 15 October, 2021; originally announced October 2021.

Solving Rubik's Cube via Quantum Mechanics and Deep Reinforcement Learning

Authors: Sebastiano Corli, Lorenzo Moro, Davide E. Galli, Enrico Prati

Abstract: Rubik's Cube is one of the most famous combinatorial puzzles involving nearly $4.3 \times 10^{19}$ possible configurations. Its mathematical description is expressed by the Rubik's group, whose elements define how its layers rotate. We develop a unitary representation of such group and a quantum formalism to describe the Cube from its geometrical constraints. Cubies are describedby single particle… ▽ More Rubik's Cube is one of the most famous combinatorial puzzles involving nearly $4.3 \times 10^{19}$ possible configurations. Its mathematical description is expressed by the Rubik's group, whose elements define how its layers rotate. We develop a unitary representation of such group and a quantum formalism to describe the Cube from its geometrical constraints. Cubies are describedby single particle states which turn out to behave like bosons for corners and fermions for edges, respectively. When in its solved configuration, the Cube, as a geometrical object, shows symmetrieswhich are broken when driven away from this configuration. For each of such symmetries, we build a Hamiltonian operator. When a Hamiltonian lies in its ground state, the respective symmetry of the Cube is preserved. When all such symmetries are preserved, the configuration of the Cube matches the solution of the game. To reach the ground state of all the Hamiltonian operators, we make use of a Deep Reinforcement Learning algorithm based on a Hamiltonian reward. The Cube is solved in four phases, all based on a respective Hamiltonian reward based on its spectrum, inspired by the Ising model. Embedding combinatorial problems into the quantum mechanics formalism suggests new possible algorithms and future implementations on quantum hardware. △ Less

Submitted 15 September, 2021; originally announced September 2021.

Comments: 26 pages, 48 references, 55 equations, 5 tables, 7 figures

Journal ref: Journal of Physics A: Mathematical and Theoretical, 2021

Evolution of static and dynamical density correlations in a one-dimensional soft-core gas from the Tonks-Girardeau limit to a clustering fluid

Authors: Martina Teruzzi, Christian Apostoli, Davide Pini, Davide Emilio Galli, Gianluca Bertaina

Abstract: Repulsive soft-core atomic systems may undergo clustering if their density is high enough that core overlap is unavoidable. In one-dimensional quantum systems, it has been shown that this instability triggers a transition from a Luttinger liquid to various cluster Luttinger liquids. Here, we focus on the Luttinger liquid regime and theoretically study the evolution of key observables related to de… ▽ More Repulsive soft-core atomic systems may undergo clustering if their density is high enough that core overlap is unavoidable. In one-dimensional quantum systems, it has been shown that this instability triggers a transition from a Luttinger liquid to various cluster Luttinger liquids. Here, we focus on the Luttinger liquid regime and theoretically study the evolution of key observables related to density fluctuations, that manifest a striking dependence on density. We tune the interaction so that the low-density regime corresponds to a Tonks-Girardeau gas, and show that as the density is increased the system departs more and more from Tonks-Girardeau behavior, displaying a much larger compressibility as well as rotonic excitations that finally drive the clustering transition. We compare various theoretical approaches, which are accurate in different regimes. Using quantum Monte Carlo methods and analytic continuation as a benchmark, we investigate the regime of validity of the mean-field Bogoliubov and the real-time multiconfiguration time-dependent Hartree-Fock approaches. Part of the behavior that we describe should be observable in ultracold Rydberg-dressed gases, provided that system losses are prevented. △ Less

Submitted 22 June, 2021; originally announced June 2021.

Comments: 9 pages, 8 figures

Quantum simulations of molecular systems with intrinsic atomic orbitals

Authors: Stefano Barison, Davide Emilio Galli, Mario Motta

Abstract: Quantum simulations of molecular systems on quantum computers often employ minimal basis sets of Gaussian orbitals. In comparison with more realistic basis sets, quantum simulations employing minimal basis sets require fewer qubits and quantum gates, but yield results of lower accuracy. A natural strategy to achieve more accurate results is to increase the basis set size, which in turn requires in… ▽ More Quantum simulations of molecular systems on quantum computers often employ minimal basis sets of Gaussian orbitals. In comparison with more realistic basis sets, quantum simulations employing minimal basis sets require fewer qubits and quantum gates, but yield results of lower accuracy. A natural strategy to achieve more accurate results is to increase the basis set size, which in turn requires increasing the number of qubits and quantum gates. Here we explore the use of intrinsic atomic orbitals (IAOs) in quantum simulations of molecules, to improve the accuracy of energies and properties at the same computational cost required by a minimal basis. We investigate ground-state energies and one- and two-body density operators in the framework of the variational quantum eigensolver, employing and comparing different Ansätze. We also demonstrate the use of this approach in the calculation of ground- and excited-states energies of small molecules by a combination of quantum algorithms, using IBM Quantum computers. △ Less

Submitted 16 August, 2022; v1 submitted 16 November, 2020; originally announced November 2020.

Comments: 14 pages, 14 figures. Published in Physical Review A

Journal ref: Phys. Rev. A 106, 022404 (2022)

Emergence of an Ising critical regime in the clustering of 1D soft matter revealed through string variables

Authors: Francesco Mambretti, Simone Molinelli, Davide Pini, Gianluca Bertaina, Davide Emilio Galli

Abstract: Soft matter systems are renowned for being able to display complex emerging phenomena such as clustering phases. Recently, a surprising quantum phase transition has been revealed in a one-dimensional (1D) system composed of bosons interacting via a pairwise soft potential in the continuum. It was shown that the spatial coordinates undergoing two-particle clustering could be mapped into quantum spi… ▽ More Soft matter systems are renowned for being able to display complex emerging phenomena such as clustering phases. Recently, a surprising quantum phase transition has been revealed in a one-dimensional (1D) system composed of bosons interacting via a pairwise soft potential in the continuum. It was shown that the spatial coordinates undergoing two-particle clustering could be mapped into quantum spin variables of a 1D transverse Ising model. In this work we investigate the manifestation of an analogous critical phenomenon in 1D classical fluids of soft particles in the continuum. In particular, we study the low-temperature behavior of three different classical models of 1D soft matter, whose inter-particle interactions allow for clustering. The same string variables highlight that, at the commensurate density for the two-particle cluster phase, the peculiar pairing of neighboring soft particles can be nontrivially mapped onto a 1D discrete classical Ising model. We also observe a related phenomenon, namely the presence of an anomalous peak in the low-temperature specific heat, thus indicating the emergence of Schottky phenomenology in a non-magnetic fluid. △ Less

Submitted 14 October, 2020; v1 submitted 31 July, 2020; originally announced July 2020.

Comments: 18 figures, 12 pages Accepted by Phys. Rev. E. Post-print version

Journal ref: Phys. Rev. E 102, 042134 (2020)

Crystal growth rates in supercooled atomic liquid mixtures

Authors: Alexander Schottelius, Francesco Mambretti, Anton Kalinin, Björn Beyersdorff, Andre Rothkirch, Claudia Goy, Jan Müller, Nikolaos Petridis, Maurizio Ritzer, Florian Trinter, José M. Fernández, Tiberio A. Ezquerra, Davide E. Galli, Robert E. Grisenti

Abstract: Crystallization is a fundamental process in materials science, providing the primary route for the realization of a wide range of new materials. Crystallization rates are also considered to be useful probes of glass-forming ability. At the microscopic level, crystallization is described by the classical crystal nucleation and growth theories, yet in general solid formation is a far more complex pr… ▽ More Crystallization is a fundamental process in materials science, providing the primary route for the realization of a wide range of new materials. Crystallization rates are also considered to be useful probes of glass-forming ability. At the microscopic level, crystallization is described by the classical crystal nucleation and growth theories, yet in general solid formation is a far more complex process. In particular, the observation of apparently different crystal growth regimes in many binary liquid mixtures greatly challenges our understanding of crystallization. Here, we study by experiments, theory and computer simulations the crystallization of supercooled mixtures of argon and krypton, showing that crystal growth rates in these systems can be reconciled with existing crystal growth models only by explicitly accounting for the non-ideality of the mixtures. Our results highlight the importance of thermodynamic aspects in describing the crystal growth kinetics, providing a substantial step towards a more sophisticated theory of crystal growth. △ Less

Submitted 8 March, 2020; originally announced March 2020.

Comments: 19 pages, 3 figures

Journal ref: Nature Materials 19, 512 (2020)

Dynamical structure factor of a fermionic supersolid on an optical lattice

Authors: Ettore Vitali, Patrick Kelly, Annette Lopez, Gianluca Bertaina, Davide Emilio Galli

Abstract: Interfacing unbiased quantum Monte Carlo simulations with state-of-art analytic continuation techniques, we obtain exact numerical results for dynamical density and spin correlations in the attractive Hubbard model, describing a spin-balanced two-dimensional cold Fermi gas on an optical lattice. We focus on half-filling, where on average one fermion occupies each lattice site, and the system displ… ▽ More Interfacing unbiased quantum Monte Carlo simulations with state-of-art analytic continuation techniques, we obtain exact numerical results for dynamical density and spin correlations in the attractive Hubbard model, describing a spin-balanced two-dimensional cold Fermi gas on an optical lattice. We focus on half-filling, where on average one fermion occupies each lattice site, and the system displays an intriguing supersolid phase: a superfluid with a checkerboard density modulation. The coexistence of $U(1)$ broken symmetry and the density modulations makes this regime very challenging and interesting for the calculation of dynamical properties. We compare our unbiased results with state-of-the-art Generalized Random Phase Approximation calculations: both approaches agree on a well-defined low-energy Nambu-Goldstone collective mode in the density correlations, while the higher energy structures appear to differ significantly. We also observe an interesting high-energy spin mode. We argue that our results provide a robust benchmark for Generalized Random Phase Approximation techniques, which are widely considered to be the method of choice for dynamical correlations in Fermi gases. Also, our calculations yield new physical insight in the high-energy behavior of the dynamical structure factor of the attractive Hubbard model, which is a well known prototype lattice model for superconductors and is a fertile field to target the observation of collective modes in strongly correlated systems. △ Less

Submitted 4 December, 2020; v1 submitted 12 December, 2019; originally announced December 2019.

Comments: 12 pages, 6 figures, post-print version

Journal ref: Phys. Rev. A 102, 053324 (2020)

Static density response of one-dimensional soft bosons across the clustering transition

Authors: Martina Teruzzi, Davide Pini, Stefano Rossotti, Davide Emilio Galli, Gianluca Bertaina

Abstract: One-dimensional bosons interacting via a soft-shoulder potential are investigated at zero temperature. The flatness of the potential at short distances introduces a typical length, such that, at relatively high densities and sufficiently strong interactions, clusters are formed, even in the presence of a completely repulsive potential. We evaluate the static density response function of this syste… ▽ More One-dimensional bosons interacting via a soft-shoulder potential are investigated at zero temperature. The flatness of the potential at short distances introduces a typical length, such that, at relatively high densities and sufficiently strong interactions, clusters are formed, even in the presence of a completely repulsive potential. We evaluate the static density response function of this system across the transition from the liquid to the cluster liquid phases. Such quantity reveals the density modulations induced by a weak periodic external potential, and is maximal at the clustering wavevector. It is known that this response function is proportional to the static structure factor in the classical regime at high temperature, while for this zero-temperature quantum systems, we extract it from the dynamical structure factor evaluated with quantum Monte Carlo methods. △ Less

Submitted 9 April, 2018; originally announced April 2018.

Comments: 7 pages, 3 figures. Proceedings of the "Recent Progress in Many-Body Theories" international conference, Pohang, South Korea, 2017

Probing quantum turbulence in $^4$He by quantum evaporation measurements

Authors: Ivan Amelio, Davide Emilio Galli, Luciano Reatto

Abstract: Theory of superfluid $^4$He shows that, due to strong correlations and backflow effects, the density profile of a vortex line has the character of a density modulation and it is not a simple rarefaction region as found in clouds of cold bosonic atoms. We find that the basic features of this density modulation are represented by a wave--packet of cylindrical symmetry in which rotons with positive g… ▽ More Theory of superfluid $^4$He shows that, due to strong correlations and backflow effects, the density profile of a vortex line has the character of a density modulation and it is not a simple rarefaction region as found in clouds of cold bosonic atoms. We find that the basic features of this density modulation are represented by a wave--packet of cylindrical symmetry in which rotons with positive group velocity have a dominant role: The vortex density modulation can be viewed as a cloud of virtual excitations, mainly rotons, sustained by the phase of the vortex wave function. This suggests that in a vortex reconnection some of these rotons become real so that a vortex tangle is predicted to be a source of non-thermal rotons. The presence of such vorticity induced rotons can be verified by measurements at low temperature of quantum evaporation of $^4$He atoms. We estimate the rate of evaporation and this turns out to be detectable by current instrumentation. Additional information on the microscopic processes in the decay of quantum turbulence will be obtained if quantum evaporation by high energy phonons should be detected. △ Less

Submitted 20 February, 2018; originally announced February 2018.

Comments: 9 pages, 7 figures

Journal ref: Phys. Rev. Lett. 121, 015302 (2018)

Quantum Critical Behavior of One-Dimensional Soft Bosons in the Continuum

Authors: Stefano Rossotti, Martina Teruzzi, Davide Pini, Davide Emilio Galli, Gianluca Bertaina

Abstract: We consider a zero-temperature one-dimensional system of bosons interacting via the soft-shoulder potential in the continuum, typical of dressed Rydberg gases. We employ quantum Monte Carlo simulations, which allow for the exact calculation of imaginary-time correlations, and a stochastic analytic continuation method, to extract the dynamical structure factor. At finite densities, in the weakly-in… ▽ More We consider a zero-temperature one-dimensional system of bosons interacting via the soft-shoulder potential in the continuum, typical of dressed Rydberg gases. We employ quantum Monte Carlo simulations, which allow for the exact calculation of imaginary-time correlations, and a stochastic analytic continuation method, to extract the dynamical structure factor. At finite densities, in the weakly-interacting homogeneous regime, a rotonic spectrum marks the tendency to clustering. With strong interactions, we indeed observe cluster liquid phases emerging, characterized by the spectrum of a composite harmonic chain. Luttinger theory has to be adapted by changing the reference lattice density field. In both the liquid and cluster liquid phases, we find convincing evidence of a secondary mode, which becomes gapless only at the transition. In that region, we also measure the central charge and observe its increase towards c = 3/2, as recently evaluated in a related extended Bose-Hubbard model, and we note a fast reduction of the Luttinger parameter. For 2-particle clusters, we then interpret such observations in terms of the compresence of a Luttinger liquid and a critical transverse Ising model, related to the instability of the reference lattice density field towards coalescence of sites, typical of potentials which are flat at short distances. Even in the absence of a true lattice, we are able to evaluate the spatial correlation function of a suitable pseudo-spin operator, which manifests ferromagnetic order in the cluster liquid phase, exponential decay in the liquid phase, and algebraic order at criticality. △ Less

Submitted 25 November, 2017; v1 submitted 27 July, 2017; originally announced July 2017.

Comments: 11 pages, 13 figures. Added one reference, cleared discussion on liquid phase, results unchanged

Journal ref: Phys. Rev. Lett. 119, 215301 (2017)

Facing the phase problem in Coherent Diffractive Imaging via Memetic Algorithms

Authors: Alessandro Colombo, Davide Emilio Galli, Liberato De Caro, Francesco Scattarella, Elvio Carlino

Abstract: Coherent Diffractive Imaging is a lensless technique that allows imaging of matter at a spatial resolution not limited by lens aberrations. This technique exploits the measured diffraction pattern of a coherent beam scattered by periodic and non-periodic objects to retrieve spatial information. The diffracted intensity, for weak-scattering objects, is proportional to the modulus of the Fourier Tra… ▽ More Coherent Diffractive Imaging is a lensless technique that allows imaging of matter at a spatial resolution not limited by lens aberrations. This technique exploits the measured diffraction pattern of a coherent beam scattered by periodic and non-periodic objects to retrieve spatial information. The diffracted intensity, for weak-scattering objects, is proportional to the modulus of the Fourier Transform of the object scattering function. Any phase information, needed to retrieve its scattering function, has to be retrieved by means of suitable algorithms. Here we present a new approach, based on a memetic algorithm, i.e. a hybrid genetic algorithm, to face the phase problem, which exploits the synergy of deterministic and stochastic optimization methods. The new approach has been tested on simulated data and applied to the phasing of transmission electron microscopy coherent electron diffraction data of a $\text{SrTiO}_\text{3}$ sample. We have been able to quantitatively retrieve the projected atomic potential, and also image the oxygen columns, which are not directly visible in the relevant high-resolution transmission electron microscopy images. Our approach proves to be a new powerful tool for the study of matter at atomic resolution and opens new perspectives in those applications in which effective phase retrieval is necessary. △ Less

Submitted 26 January, 2017; originally announced January 2017.

Comments: 11 pages, 8 figures. To be published on Sci. Rep. 7, 42236 (2017)

Statistical and computational intelligence approach to analytic continuation in Quantum Monte Carlo

Authors: G. Bertaina, D. E. Galli, E. Vitali

Abstract: The term analytic continuation emerges in many branches of Mathematics, Physics, and, more generally, applied Science. Generally speaking, in many situations, given some amount of information that could arise from experimental or numerical measurements, one is interested in extending the domain of such information, to infer the values of some variables which are central for the study of a given pr… ▽ More The term analytic continuation emerges in many branches of Mathematics, Physics, and, more generally, applied Science. Generally speaking, in many situations, given some amount of information that could arise from experimental or numerical measurements, one is interested in extending the domain of such information, to infer the values of some variables which are central for the study of a given problem. For example, focusing on Condensed Matter Physics, state-of-the-art methodologies to study strongly correlated quantum physical systems are able to yield accurate estimations of dynamical correlations in imaginary time. Those functions have to be extended to the whole complex plane, via analytic continuation, in order to infer real-time properties of those physical systems. In this Review, we will present the Genetic Inversion via Falsification of Theories method, which allowed us to compute dynamical properties of strongly interacting quantum many--body systems with very high accuracy. Even though the method arose in the realm of Condensed Matter Physics, it provides a very general framework to face analytic continuation problems that could emerge in several areas of applied Science. Here we provide a pedagogical review that elucidates the approach we have developed. △ Less

Submitted 6 March, 2017; v1 submitted 25 November, 2016; originally announced November 2016.

Comments: 18 pages, 3 figures. This is an Accepted Manuscript of an article published by Taylor & Francis in Advances in Physics X on 19 Feb 2017 (online), available online: http://www.tandfonline.com/doi/full/10.1080/23746149.2017.1288585

Journal ref: Advances in Physics: X 2, 302 (2017)

Dynamical structure factor of one-dimensional hard rods

Authors: M. Motta, E. Vitali, M. Rossi, D. E. Galli, G. Bertaina

Abstract: The zero-temperature dynamical structure factor $S(q,ω)$ of one-dimensional hard rods is computed using state-of-the-art quantum Monte Carlo and analytic continuation techniques, complemented by a Bethe Ansatz analysis. As the density increases, $S(q,ω)$ reveals a crossover from the Tonks-Girardeau gas to a quasi-solid regime, along which the low-energy properties are found in agreement with the n… ▽ More The zero-temperature dynamical structure factor $S(q,ω)$ of one-dimensional hard rods is computed using state-of-the-art quantum Monte Carlo and analytic continuation techniques, complemented by a Bethe Ansatz analysis. As the density increases, $S(q,ω)$ reveals a crossover from the Tonks-Girardeau gas to a quasi-solid regime, along which the low-energy properties are found in agreement with the nonlinear Luttinger liquid theory. Our quantitative estimate of $S(q,ω)$ extends beyond the low-energy limit and confirms a theoretical prediction regarding the behavior of $S(q,ω)$ at specific wavevectors $\mathcal{Q}_n=n 2 π/a$, where $a$ is the core radius, resulting from the interplay of the particle-hole boundaries of suitably rescaled ideal Fermi gases. We observe significant similarities between hard rods and one-dimensional $^4$He at high density, suggesting that the hard-rods model may provide an accurate description of dense one-dimensional liquids of quantum particles interacting through a strongly repulsive, finite-range potential. △ Less

Submitted 27 August, 2016; originally announced August 2016.

Comments: 13 pages, 9 figures

Journal ref: Phys. Rev. A 94, 043627 (2016)

Microscopic study of static and dynamical properties of dilute one-dimensional soft bosons

Authors: M. Teruzzi, D. E. Galli, G. Bertaina

Abstract: We study static properties and the dynamical structure factor of zero-temperature dilute bosons interacting via a soft-shoulder potential in one dimension. Our approach is fully microscopic and employs state-of-the-art quantum Monte Carlo and analytic continuation techniques. By increasing the interaction strength, our model reproduces the Lieb-Liniger gas, the Tonks-Girardeau and the Hard-Rods mo… ▽ More We study static properties and the dynamical structure factor of zero-temperature dilute bosons interacting via a soft-shoulder potential in one dimension. Our approach is fully microscopic and employs state-of-the-art quantum Monte Carlo and analytic continuation techniques. By increasing the interaction strength, our model reproduces the Lieb-Liniger gas, the Tonks-Girardeau and the Hard-Rods models. △ Less

Submitted 11 January, 2017; v1 submitted 18 July, 2016; originally announced July 2016.

Comments: 8 pages, 5 figures, published on J. Low. Temp. Phys, Special Issue QFS 2016. Updated with hard-rods regime results obtained with a better initial wavefunction

Linear response of one-dimensional liquid $^4$He to external perturbations

Authors: M. Motta, G. Bertaina, E. Vitali, D. E. Galli, M. Rossi

Abstract: We study the response of one-dimensional liquid $^4$He to weak perturbations relying on the dynamical structure factor, $S(q,ω)$, recently obtained via ab-initio techniques [Phys. Rev. Lett. 116, 135302 (2016)]. We evaluate the drag force, $F_v$, experienced by an impurity moving along the system with velocity $v$ and the static response function, $χ(q)$, describing the density modulations induced… ▽ More We study the response of one-dimensional liquid $^4$He to weak perturbations relying on the dynamical structure factor, $S(q,ω)$, recently obtained via ab-initio techniques [Phys. Rev. Lett. 116, 135302 (2016)]. We evaluate the drag force, $F_v$, experienced by an impurity moving along the system with velocity $v$ and the static response function, $χ(q)$, describing the density modulations induced by a periodic perturbation with wave vector $q$. △ Less

Submitted 25 November, 2016; v1 submitted 18 July, 2016; originally announced July 2016.

Comments: 7 pages, 2 figures

Roton excitations and the fluid-solid phase transition in superfluid 2D Yukawa Bosons

Authors: S. Molinelli, D. E. Galli, L. Reatto, M. Motta

Abstract: We compute several ground state properties and the dynamical structure factor of a 0-temperature system of Bosons interacting with the 2D screened Coulomb (2D-SC) potential. We resort to the exact shadow path integral ground state (SPIGS) quantum Monte Carlo method to compute the imaginary-time correlation function of the model, and to the genetic algorithm via falsification of theories (GIFT) to… ▽ More We compute several ground state properties and the dynamical structure factor of a 0-temperature system of Bosons interacting with the 2D screened Coulomb (2D-SC) potential. We resort to the exact shadow path integral ground state (SPIGS) quantum Monte Carlo method to compute the imaginary-time correlation function of the model, and to the genetic algorithm via falsification of theories (GIFT) to retrieve the dynamical structure factor. We provide a detailed comparison of ground state properties and collective excitations of 2D-SC and 4He atoms. The roton energy of the 2D-SC system is an increasing function of density, and not a decreasing one as in 4He. This result is in contrast with the view that the roton is the soft mode of the fluid-solid transition. We uncover a remarkable quasi-universality of backflow and of other properties when expressed in terms of the amount of short range order as quantified by the height of the first peak of the static structure factor. △ Less

Submitted 4 May, 2016; originally announced May 2016.

Comments: 15 pages, 10 figures, accepted for publication in J. Low Temp. Phys

Low-density phases of $^3$He monolayers adsorbed on graphite

Authors: Michele Ruggeri, Ettore Vitali, Davide Emilio Galli, Massimo Boninsegni, Saverio Moroni

Abstract: Quantum Monte Carlo simulations at zero temperature of a $^3$He monolayer adsorbed on graphite, either clean or preplated with $^4$He, unexpectedly point to a gas-liquid phase transition at a very low areal density of the order of 0.01Å$^{-2}$. This result stems from an essentially unbiased calculation of the ground state energy for an infinite, defect-free substrate which interacts with He atoms… ▽ More Quantum Monte Carlo simulations at zero temperature of a $^3$He monolayer adsorbed on graphite, either clean or preplated with $^4$He, unexpectedly point to a gas-liquid phase transition at a very low areal density of the order of 0.01Å$^{-2}$. This result stems from an essentially unbiased calculation of the ground state energy for an infinite, defect-free substrate which interacts with He atoms via a realistic potential, whereas the interaction between two He atoms includes two- and three-body terms. The sensitivity of the gas-liquid coexistence region on the model Hamiltonian employed is discussed. △ Less

Submitted 10 February, 2016; v1 submitted 12 December, 2015; originally announced December 2015.

Comments: 8 pages, 9 figures

Journal ref: Phys. Rev. B 93, 104102 (2016)

Imaginary Time Correlations for a High-Density two-dimensional Electron Gas

Authors: M. Motta, D. E. Galli, S. Moroni, E. Vitali

Abstract: We evaluate imaginary time density-density correlation functions for a two-dimensional homogeneous electron gas using the phaseless auxiliary field quantum Monte Carlo method. 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… ▽ More We evaluate imaginary time density-density correlation functions for a two-dimensional homogeneous electron gas using the phaseless auxiliary field quantum Monte Carlo method. 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 present simulations of a number up to 42 correlated fermions in the continuum, using up to 300 plane waves as basis set elements. We discuss the numerical stabilization techniques and the computational complexity of the methodology. 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 many-fermion systems. △ Less

Submitted 28 July, 2015; originally announced September 2015.

Comments: 17 pages, 8 figures; submitted to J. Chem. Phys

Mixing Effects in the Crystallization of Supercooled Quantum Binary Liquids

Authors: M. Kühnel, J. M. Fernández, F. Tramonto, G. Tejeda, E. Moreno, A. Kalinin, M. Nava, D. E. Galli, S. Montero, R. E. Grisenti

Abstract: By means of Raman spectroscopy of liquid microjets we have investigated the crystallization process of supercooled quantum liquid mixtures composed of parahydrogen (pH$_2$) diluted with small amounts of up to 5\% of either neon or orthodeuterium (oD$_2$), and of oD$_2$ diluted with either Ne or pH$_2$. We show that the introduction of Ne impurities affects the crystallization kinetics in both the… ▽ More By means of Raman spectroscopy of liquid microjets we have investigated the crystallization process of supercooled quantum liquid mixtures composed of parahydrogen (pH$_2$) diluted with small amounts of up to 5\% of either neon or orthodeuterium (oD$_2$), and of oD$_2$ diluted with either Ne or pH$_2$. We show that the introduction of Ne impurities affects the crystallization kinetics in both the pH$_2$-Ne and oD$_2$-Ne mixtures in terms of a significant reduction of the crystal growth rate, similarly to what found in our previous work on supercooled pH$_2$-oD$_2$ liquid mixtures [M. Kühnel et {\it al.}, Phys. Rev. B \textbf{89}, 180506(R) (2014)]. Our experimental results, in combination with path-integral simulations of the supercooled liquid mixtures, suggest in particular a correlation between the measured growth rates and the ratio of the effective particle sizes originating from quantum delocalization effects. We further show that the crystalline structure of the mixture is also affected to a large extent by the presence of the Ne impurities, which likely initiate the freezing process through the formation of Ne crystallites. △ Less

Submitted 28 April, 2015; originally announced April 2015.

Comments: 19 pages, 7 figures, submitted to J. Chem. Phys

One-dimensional liquid $^4$He: dynamical properties beyond Luttinger liquid theory

Authors: G. Bertaina, M. Motta, M. Rossi, E. Vitali, D. E. Galli

Abstract: We compute the zero-temperature dynamical structure factor of one-dimensional liquid $^4$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 quasi-solid regime. In the low-energy limit, the dynamical structure factor can be described b… ▽ More We compute the zero-temperature dynamical structure factor of one-dimensional liquid $^4$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 quasi-solid 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 {\it{particle-hole}} continuum typical of fermionic systems. At the low-energy boundary of such region and moderate densities, we find consistency, within statistical uncertainties, with predictions of a power-law structure by the recently-developed non-linear Luttinger liquid theory. In the quasi-solid 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. △ Less

Submitted 10 March, 2016; v1 submitted 22 December, 2014; originally announced December 2014.

Comments: 12 pages, 9 figures, containing Supplemental Material. Accepted for publication in Physical Review Letters. Completely revised text, stricter finite-size analysis, more sophisticate analytic continuations, new Fig. 4 showing comparison to non-linear Luttinger liquid theory for hard-rods

Journal ref: Phys. Rev. Lett. 116, 135302 (2016)

Condensed phase of Bose-Fermi mixtures with a pairing interaction

Authors: Andrea Guidini, Gianluca Bertaina, Davide Emilio Galli, Pierbiagio Pieri

Abstract: We study the condensed phase of a Bose-Fermi mixture with a tunable pairing interaction between bosons and fermions with many-body diagrammatic methods and fixed-node diffusion Quantum Monte Carlo simulations. A universal behavior of the condensate fraction and bosonic momentum distribution with respect to the boson concentration is found to hold in an extended range of boson-fermion couplings and… ▽ More We study the condensed phase of a Bose-Fermi mixture with a tunable pairing interaction between bosons and fermions with many-body diagrammatic methods and fixed-node diffusion Quantum Monte Carlo simulations. A universal behavior of the condensate fraction and bosonic momentum distribution with respect to the boson concentration is found to hold in an extended range of boson-fermion couplings and concentrations. For vanishing boson density, we prove that the bosonic condensate fraction reduces to the quasiparticle weight Z of the Fermi polaron studied in the context of polarized Fermi gases, unifying in this way two apparently unrelated quantities. △ Less

Submitted 9 February, 2015; v1 submitted 8 December, 2014; originally announced December 2014.

Comments: 5 pages, 4 figures; final version as published

Journal ref: Physical Review A 91, 023603 (2015)

Implementation of the Linear Method for the optimization of Jastrow-Feenberg and Backflow Correlations

Authors: M. Motta, G. Bertaina, D. E. Galli, E. Vitali

Abstract: We present a fully detailed and highly performing implementation of the Linear Method [J. Toulouse and C. J. Umrigar (2007)] to optimize Jastrow-Feenberg and Backflow Correlations in many-body wave-functions, which are widely used in condensed matter physics. We show that it is possible to implement such optimization scheme performing analytical derivatives of the wave-function with respect to the… ▽ More We present a fully detailed and highly performing implementation of the Linear Method [J. Toulouse and C. J. Umrigar (2007)] to optimize Jastrow-Feenberg and Backflow Correlations in many-body wave-functions, which are widely used in condensed matter physics. We show that it is possible to implement such optimization scheme performing analytical derivatives of the wave-function with respect to the variational parameters achieving the best possible complexity O(N^3) in the number of particles N. △ Less

Submitted 10 December, 2014; v1 submitted 1 December, 2014; originally announced December 2014.

Comments: submitted to the Comp. Phys. Comm

Path Integral Monte Carlo study confirms a highly ordered snowball in 4He nanodroplets doped with an Ar+ ion

Authors: Filippo Tramonto, Paolo Salvestrini, Marco Nava, Davide E. Galli

Abstract: By means of the exact Path Integral Monte Carlo method we have performed a detailed microscopic study of $^4$He nanodroplets doped with an argon ion, Ar$^+$, at $T=0.5$ K. We have computed density profiles, energies, dissociation energies and characterized the local order around the ion for nanodroplets with a number of 4He atoms ranging from 10 to 64 and also 128. We have found the formation of a… ▽ More By means of the exact Path Integral Monte Carlo method we have performed a detailed microscopic study of $^4$He nanodroplets doped with an argon ion, Ar$^+$, at $T=0.5$ K. We have computed density profiles, energies, dissociation energies and characterized the local order around the ion for nanodroplets with a number of 4He atoms ranging from 10 to 64 and also 128. We have found the formation of a stable solid structure around the ion, a "snowball", consisting of 3 concentric shells in which the 4He atoms are placed on at the vertices of platonic solids: the first inner shell is an icosahedron (12 atoms); the second one is a dodecahedron with 20 atoms placed on the faces of the icosahedron of the first shell; the third shell is again an icosahedron composed of 12 atoms placed on the faces of the dodecahedron of the second shell. The "magic numbers" implied by this structure, 12, 32 and 44 helium atoms, have been observed in a recent experimental study [Bartl et al, J. Phys. Chem. A 118, 2014] of these complexes; the dissociation energy curve computed in the present work shows jumps in correspondence with those found in the nanodroplets abundance distribution measured in that experiment, strengthening the agreement between theory and experiment. The same structures were predicted in Ref. [Galli et al, J. Phys. Chem. A 115, 2011] in a study regarding Na$^+$@$^4$He$_n$ when n>30; a comparison between Ar$^+$@$^4$He$_n$ and Na$^+$@$^4$He$_n$ complexes is also presented. △ Less

Submitted 17 October, 2014; v1 submitted 10 October, 2014; originally announced October 2014.

Comments: 10 pages, 4 figures

Observation of crystallization slowdown in supercooled para-hydrogen and ortho-deuterium quantum liquid mixtures

Authors: Matthias Kühnel, José M. Fernández, Filippo Tramonto, Guzmán Tejeda, Elena Moreno, Anton Kalinin, Marco Nava, Davide E. Galli, Salvador Montero, Robert E. Grisenti

Abstract: We report a quantitative experimental study of the crystallization kinetics of supercooled quantum liquid mixtures of para-hydrogen (pH$_2$) and ortho-deuterium (oD$_2$) by high spatial resolution Raman spectroscopy of liquid microjets. We show that in a wide range of compositions the crystallization rate of the isotopic mixtures is significantly reduced with respect to that of the pure substances… ▽ More We report a quantitative experimental study of the crystallization kinetics of supercooled quantum liquid mixtures of para-hydrogen (pH$_2$) and ortho-deuterium (oD$_2$) by high spatial resolution Raman spectroscopy of liquid microjets. We show that in a wide range of compositions the crystallization rate of the isotopic mixtures is significantly reduced with respect to that of the pure substances. To clarify this behavior we have performed path-integral simulations of the non-equilibrium pH$_2$-oD$_2$ liquid mixtures, revealing that differences in quantum delocalization between the two isotopic species translate into different effective particle sizes. Our results provide first experimental evidence for crystallization slowdown of quantum origin, offering a benchmark for theoretical studies of quantum behavior in supercooled liquids. △ Less

Submitted 10 October, 2014; originally announced October 2014.

Comments: 6 pages, 3 figures

Journal ref: Phys. Rev. B 89, 180201(R) (2014)

Quantum Monte Carlo study of a vortex in superfluid $^4$He and search for a vortex state in the solid

Authors: Davide E. Galli, Luciano Reatto, Maurizio Rossi

Abstract: We have performed a microscopic study of a straight quantized vortex line in three dimensions in condensed $^4$He at zero temperature using the Shadow Path Integral Ground State method and the fixed-phase approximation. We have characterized the energy and the local density profile around the vortex axis in superfluid $^4$He at several densities, ranging from below the equilibrium density up to th… ▽ More We have performed a microscopic study of a straight quantized vortex line in three dimensions in condensed $^4$He at zero temperature using the Shadow Path Integral Ground State method and the fixed-phase approximation. We have characterized the energy and the local density profile around the vortex axis in superfluid $^4$He at several densities, ranging from below the equilibrium density up to the overpressurized regime. For the Onsager-Feynman (OF) phase our results are exact and represent a benchmark for other theories. The inclusion of backflow correlations in the phase improves the description of the vortex with respect to the OF phase by a large reduction of the core energy of the topological excitation. At all densities the phase with backflow induces a partial filling of the vortex core and this filling slightly increases with density. The core size slightly decreases for increasing density and the density profile has well defined density dependent oscillations whose wave vector is closer to the wave vector of the main peak in the static density response function rather than to the roton wave vector. Our results can be applied to vortex rings of large radius $R$ and we find good agreement with the experimental value of the energy as function of $R$ without any free parameter. We have studied also $^4$He above the melting density in the solid phase using the same functional form for the phase as in the liquid. We found that off-diagonal properties of the solid are not qualitatively affected by the velocity field induced by the vortex phase, both with and without backflow correlations. Therefore we find evidence that a perfect $^4$He crystal is not a marginally stable quantum solid in which rotation would be able to induce off-diagonal long-range coherence. △ Less

Submitted 28 June, 2014; v1 submitted 29 May, 2014; originally announced May 2014.

Comments: 15 pages, 8 figures

Journal ref: Phys. Rev. B 89, 224516 (2014)

Imaginary Time Correlations and the phaseless Auxiliary Field Quantum Monte Carlo

Authors: M. Motta, D. E. Galli, S. Moroni, E. Vitali

Abstract: The phaseless Auxiliary Field Quantum Monte Carlo method provides a well established approximation scheme for accurate calculations of ground state energies of many-fermions systems. Here we apply the method to the calculation of imaginary time correlation functions. We give a detailed description of the technique and we test the quality of the results for static and dynamic properties against exa… ▽ More The phaseless Auxiliary Field Quantum Monte Carlo method provides a well established approximation scheme for accurate calculations of ground state energies of many-fermions systems. Here we apply the method to the calculation of imaginary time correlation functions. We give a detailed description of the technique and we test the quality of the results for static and dynamic properties against exact values for small systems. △ Less

Submitted 5 November, 2013; originally announced November 2013.

Comments: 13 pages, 6 figures; submitted to J. Chem. Phys

Quantum Monte Carlo study of the dynamic structure factor in the gas and crystal phase of hard-sphere bosons

Authors: Riccardo Rota, Filippo Tramonto, Davide Emilio Galli, Stefano Giorgini

Abstract: We investigate the dynamic structure factor of a system of Bose particles at zero temperature using quantum Monte Carlo methods. Interactions are modeled using a hard-sphere potential of size $a$ and simulations are performed for values of the gas parameter $na^3$ ranging from the dilute regime up to densities $n$ where the thermodynamically stable phase is a solid. With increasing density we obse… ▽ More We investigate the dynamic structure factor of a system of Bose particles at zero temperature using quantum Monte Carlo methods. Interactions are modeled using a hard-sphere potential of size $a$ and simulations are performed for values of the gas parameter $na^3$ ranging from the dilute regime up to densities $n$ where the thermodynamically stable phase is a solid. With increasing density we observe a crossover of the dispersion of elementary excitations from a Bogoliubov-like spectrum to a phonon--maxon--roton curve and the emergence of a broad multiphonon contribution accompanying the single-quasiparticle peak. In particular, for $na^3=0.2138$, which corresponds to superfluid $^4$He at equilibrium density, the extracted spectrum turns out to be in good agreement with the experimental energy--momentum dispersion relation in the roton region and for higher momenta. The behavior of the spectral function at the same density in the stable solid and metastable gas phase above the freezing point is also discussed. △ Less

Submitted 7 October, 2013; originally announced October 2013.

Journal ref: Phys. Rev. B 88, 214505 (2013)

Novel behaviors of monolayer quantum gases on Graphene, Graphane and Fluorographene

Authors: Luciano Reatto, Davide E. Galli, Marco Nava, Milton W. Cole

Abstract: This article discusses the behavior of submonolayer quantum films (He and H2) on graphene and newly discovered surfaces that are derived from graphene. Among these substrates are graphane (abbreviated GH), which has an H atom bonded to each C atom, and fluorographene (GF). For the graphite case, extensive experimental and theoretical investigations have revealed that the phase diagrams of the bose… ▽ More This article discusses the behavior of submonolayer quantum films (He and H2) on graphene and newly discovered surfaces that are derived from graphene. Among these substrates are graphane (abbreviated GH), which has an H atom bonded to each C atom, and fluorographene (GF). For the graphite case, extensive experimental and theoretical investigations have revealed that the phase diagrams of the bose gases 4He and para-H2 are qualitatively similar, differing primarily in a higher characteristic temperature scale for H2 than for He. The phase behaviors of these films on one side of pristine graphene, or both sides of free-standing graphene, are expected to be similar to those on graphite. We point out the possibility of novel phenomena in adsorption on graphene related to the large flexibility of the graphene sheet, to the non--negligible interaction between atoms adsorbed on opposite sides of the sheet and to the perturbation effect of the adsorbed layer on the Dirac electrons. In contrast, the behaviors predicted on GF and GH surfaces are very different from those on graphite, a result of the different corrugation, i.e., the lateral variation of the potential experienced by these gases. As a result of this novel geometry and potential, distinct properties are observed. For example, the 4He film's ground-state on graphite is a two-dimensional (2D) crystal commensurate with the substrate; on GF and GH, instead, it is predicted to be an anisotropic superfluid. On GF the anisotropy is so extreme that the roton excitations are very anisotropic, as if the bosons are moving in a multiconnected space along the bonds of a honeycomb lattice. Such a novel phase has not been predicted or observed previously on any substrate. Also, in the case of 3He the film's ground-state is a fluid, thus offering the possibility of studying an anisotropic Fermi fluid with a tunable density. △ Less

Submitted 19 September, 2013; originally announced September 2013.

Comments: 32 pages, 18 figures, accepted for publication in Journal of Physics: Condensed Matter as a topical review

Excitation spectrum in two-dimensional superfluid 4He

Authors: F. Arrigoni, E. Vitali, D. E. Galli, L. Reatto

Abstract: In this work we perform an ab-initio study of an ideal two-dimensional sample of 4He atoms, a model for 4He films adsorbed on several kinds of substrates. Starting from a realistic hamiltonian we face the microscopic study of the excitation phonon-roton spectrum of the system at zero temperature. Our approach relies on Path Integral Ground State Monte Carlo projection methods, allowing to evaluate… ▽ More In this work we perform an ab-initio study of an ideal two-dimensional sample of 4He atoms, a model for 4He films adsorbed on several kinds of substrates. Starting from a realistic hamiltonian we face the microscopic study of the excitation phonon-roton spectrum of the system at zero temperature. Our approach relies on Path Integral Ground State Monte Carlo projection methods, allowing to evaluate exactly the dynamical density correlation functions in imaginary time, and this gives access to the dynamical structure factor of the system S(q,omega), containing information about the excitation spectrum E(q), resulting in sharp peaks in S(q,omega). The actual evaluation of S(q,omega) requires the inversion of the Laplace transform in ill-posed conditions, which we face via the Genetic Inversion via Falsification of Theories technique. We explore the full density range from the region of spinodal decomposition to the freezing density, i.e. 0.0321 A^-2 - 0.0658 A^-2. In particular we follow the density dependence of the excitation spectrum, focusing on the low wave--vector behavior of E(q), the roton dispersion, the strength of single quasi--particle peak, Z(q), and the static density response function, chi(q). As the density increases, the dispersion E(q) at low wave--vector changes from a super-linear (anomalous dispersion) trend to a sub-linear (normal dispersion) one, anticipating the crystallization of the system; at the same time the maxon-roton structure, which is barely visible at low density, becomes well developed at high densities and the roton wave vector has a strong density dependence. Connection is made with recent inelastic neutron scattering results from highly ordered silica nanopores partially filled with 4He. △ Less

Submitted 17 May, 2013; v1 submitted 16 May, 2013; originally announced May 2013.

Comments: Accepted for publication: Special issue of FNT (Low Temp. Phys.) 75th anniversary of the discovery of superfluidity

Journal ref: Low Temp. Phys./Fiz. Niz. Temp. 39, 1021 (2013)

Dynamic structure factor for 3He in two-dimensions

Authors: M. Nava, D. E. Galli, S. Moroni, E. Vitali

Abstract: Recent neutron scattering experiments on 3He films have observed a zero-sound mode, its dispersion relation and its merging with -and possibly emerging from- the particle-hole continuum. Here we address the study of the excitations in the system via quantum Monte Carlo methods: we suggest a practical scheme to calculate imaginary time correlation functions for moderate-size fermionic systems. Comb… ▽ More Recent neutron scattering experiments on 3He films have observed a zero-sound mode, its dispersion relation and its merging with -and possibly emerging from- the particle-hole continuum. Here we address the study of the excitations in the system via quantum Monte Carlo methods: we suggest a practical scheme to calculate imaginary time correlation functions for moderate-size fermionic systems. Combined with an efficient method for analytic continuation, this scheme affords an extremely convincing description of the experimental findings. △ Less

Submitted 11 March, 2013; v1 submitted 7 February, 2013; originally announced February 2013.

Comments: 5 pages, 5 figures

Novel substrates for Helium adsorption: Graphane and Graphene-Fluoride

Authors: L. Reatto, M. Nava, D. E. Galli, C. Billman, J. O. Sofo, M. W. Cole

Abstract: The discovery of fullerenes has stimulated extensive exploration of the resulting behavior of adsorbed films. Our study addresses the planar substrates graphene-fluoride (GF) and graphane (GH) in comparison to graphene. We present initial results concerning the potential energy, energy bands and low density behavior of 4He and 3He films on such different surfaces. For example, while graphene prese… ▽ More The discovery of fullerenes has stimulated extensive exploration of the resulting behavior of adsorbed films. Our study addresses the planar substrates graphene-fluoride (GF) and graphane (GH) in comparison to graphene. We present initial results concerning the potential energy, energy bands and low density behavior of 4He and 3He films on such different surfaces. For example, while graphene presents an adsorption potential that is qualitatively similar to that on graphite, GF and GH yield potentials with different symmetry, a number of adsorption sites double that on graphene/graphite and a larger corrugation for the adatom. In the case of GF, the lowest energy band width is similar to that on graphite but the He atom has a significantly larger effective mass and the adsorption energy is about three time that on graphite. Implications concerning the monolayer phase diagram of 4He are explored with the exact path integral ground state method. A commensurate ordered state similar to the sqrt{3} x sqrt{3} R30^o state on graphite is found the be unstable both on GF and on GH. The ground states of submonolayer 4He on both GF and GH are superfluids with a Bose Einstein condensate fraction of about 10%. △ Less

Submitted 13 April, 2012; originally announced April 2012.

Comments: 6 pages, 3 figures, LT26 proceedings, accepted for publication in Journal of Physics: Conference Series

Study of solid 4He in two dimensions. The issue of zero-point defects and study of confined crystal

Authors: M. Rossi, L. Reatto, D. E. Galli

Abstract: Defects are believed to play a fundamental role in the supersolid state of 4He. We report on studies by exact Quantum Monte Carlo (QMC) simulations at zero temperature of the properties of solid 4He in presence of many vacancies, up to 30 in two dimensions (2D). In all studied cases the crystalline order is stable at least as long as the concentration of vacancies is below 2.5%. In the 2D system f… ▽ More Defects are believed to play a fundamental role in the supersolid state of 4He. We report on studies by exact Quantum Monte Carlo (QMC) simulations at zero temperature of the properties of solid 4He in presence of many vacancies, up to 30 in two dimensions (2D). In all studied cases the crystalline order is stable at least as long as the concentration of vacancies is below 2.5%. In the 2D system for a small number, n_v, of vacancies such defects can be identified in the crystalline lattice and are strongly correlated with an attractive interaction. On the contrary when n_v~10 vacancies in the relaxed system disappear and in their place one finds dislocations and a revival of the Bose-Einstein condensation. Thus, should zero-point motion defects be present in solid 4He, such defects would be dislocations and not vacancies, at least in 2D. In order to avoid using periodic boundary conditions we have studied the exact ground state of solid 4He confined in a circular region by an external potential. We find that defects tend to be localized in an interfacial region of width of about 15 A. Our computation allows to put as upper bound limit to zero--point defects the concentration 0.003 in the 2D system close to melting density. △ Less

Submitted 21 March, 2012; originally announced March 2012.

Comments: 17 pages, accepted for publication in J. Low Temp. Phys., Special Issue on Supersolids

Quantized vortices in two dimensional solid 4He

Authors: M. Rossi, D. E. Galli, P. Salvestrini, L. Reatto

Abstract: Diagonal and off-diagonal properties of 2D solid 4He systems doped with a quantized vortex have been investigated via the Shadow Path Integral Ground State method using the fixed-phase approach. The chosen approximate phase induces the standard Onsager-Feynman flow field. In this approximation the vortex acts as a static external potential and the resulting Hamiltonian can be treated exactly with… ▽ More Diagonal and off-diagonal properties of 2D solid 4He systems doped with a quantized vortex have been investigated via the Shadow Path Integral Ground State method using the fixed-phase approach. The chosen approximate phase induces the standard Onsager-Feynman flow field. In this approximation the vortex acts as a static external potential and the resulting Hamiltonian can be treated exactly with Quantum Monte Carlo methods. The vortex core is found to sit in an interstitial site and a very weak relaxation of the lattice positions away from the vortex core position has been observed. Also other properties like Bragg peaks in the static structure factor or the behavior of vacancies are very little affected by the presence of the vortex. We have computed also the one-body density matrix in perfect and defected 4He crystals finding that the vortex has no sensible effect on the off-diagonal long range tail of the density matrix. Within the assumed Onsager Feynman phase, we find that a quantized vortex cannot auto-sustain itself unless a condensate is already present like when dislocations are present. It remains to be investigated if backflow can change this conclusion. △ Less

Submitted 8 November, 2011; originally announced November 2011.

Comments: 4 pages, 3 figures, LT26 proceedings, accepted for publication in Journal of Physics: Conference Series

Equation of state of two--dimensional $^3$He at zero temperature

Authors: M. Nava, E. Vitali, A. Motta, D. E. Galli, S. Moroni

Abstract: We have performed a Quantum Monte Carlo study of a two-dimensional bulk sample of interacting 1/2-spin structureless fermions, a model of $^3$He adsorbed on a variety of preplated graphite substrates. We have computed the equation of state and the polarization energy using both the standard fixed-node approximate technique and a formally exact methodology, relying on bosonic imaginary-time correla… ▽ More We have performed a Quantum Monte Carlo study of a two-dimensional bulk sample of interacting 1/2-spin structureless fermions, a model of $^3$He adsorbed on a variety of preplated graphite substrates. We have computed the equation of state and the polarization energy using both the standard fixed-node approximate technique and a formally exact methodology, relying on bosonic imaginary-time correlation functions of operators suitably chosen in order to extract fermionic energies. As the density increases, the fixed-node approximation predicts a transition to an itinerant ferromagnetic fluid, whereas the unbiased methodology indicates that the paramagnetic fluid is the stable phase until crystallization takes place. We find that two-dimensional $^3$He at zero temperature crystallizes from the paramagnetic fluid at a density of 0.061 Å$^{-2}$ with a narrow coexistence region of about 0.002 Å$^{-2}$. Remarkably, the spin susceptibility turns out in very good agreement with experiments. △ Less

Submitted 12 January, 2012; v1 submitted 4 March, 2011; originally announced March 2011.

Comments: 7 pages, 7 figures

Exact ground state Monte Carlo method for Bosons without importance sampling

Authors: M. Rossi, M. Nava, L. Reatto, D. E. Galli

Abstract: Generally ``exact'' Quantum Monte Carlo computations for the ground state of many Bosons make use of importance sampling. The importance sampling is based, either on a guiding function or on an initial variational wave function. Here we investigate the need of importance sampling in the case of Path Integral Ground State (PIGS) Monte Carlo. PIGS is based on a discrete imaginary time evolution of… ▽ More Generally ``exact'' Quantum Monte Carlo computations for the ground state of many Bosons make use of importance sampling. The importance sampling is based, either on a guiding function or on an initial variational wave function. Here we investigate the need of importance sampling in the case of Path Integral Ground State (PIGS) Monte Carlo. PIGS is based on a discrete imaginary time evolution of an initial wave function with a non zero overlap with the ground state, that gives rise to a discrete path which is sampled via a Metropolis like algorithm. In principle the exact ground state is reached in the limit of an infinite imaginary time evolution, but actual computations are based on finite time evolutions and the question is whether such computations give unbiased exact results. We have studied bulk liquid and solid 4He with PIGS by considering as initial wave function a constant, i.e. the ground state of an ideal Bose gas. This implies that the evolution toward the ground state is driven only by the imaginary time propagator, i.e. there is no importance sampling. For both the phases we obtain results converging to those obtained by considering the best available variational wave function (the Shadow wave function) as initial wave function. Moreover we obtain the same results even by considering wave functions with the wrong correlations, for instance a wave function of a strongly localized Einstein crystal for the liquid phase. This convergence is true not only for diagonal properties such as the energy, the radial distribution function and the static structure factor, but also for off-diagonal ones, such as the one--body density matrix. From this analysis we conclude that zero temperature PIGS calculations can be as unbiased as those of finite temperature Path Integral Monte Carlo. △ Less

Submitted 20 October, 2009; v1 submitted 27 July, 2009; originally announced July 2009.

Comments: 11 pages, 10 figures

Journal ref: published J. Chem. Phys. 131, 154108 (2009)

Quantum dislocations: the fate of multiple vacancies in two dimensional solid 4He

Authors: M. Rossi, E. Vitali, D. E. Galli, L. Reatto

Abstract: Defects are believed to play a fundamental role in the supersolid state of 4He. We have studied solid 4He in two dimensions (2D) as function of the number of vacancies n_v, up to 30, inserted in the initial configuration at rho = 0.0765 A^-2, close to the melting density, with the exact zero temperature Shadow Path Integral Ground State method. The crystalline order is found to be stable also in p… ▽ More Defects are believed to play a fundamental role in the supersolid state of 4He. We have studied solid 4He in two dimensions (2D) as function of the number of vacancies n_v, up to 30, inserted in the initial configuration at rho = 0.0765 A^-2, close to the melting density, with the exact zero temperature Shadow Path Integral Ground State method. The crystalline order is found to be stable also in presence of many vacancies and we observe two completely different regimes. For small n_v, up to about 6, vacancies form a bound state and cause a decrease of the crystalline order. At larger n_v, the formation energy of an extra vacancy at fixed density decreases by one order of magnitude to about 0.6 K. In the equilibrated state it is no more possible to recognize vacancies because they mainly transform into quantum dislocations and crystalline order is found almost independent on how many vacancies have been inserted in the initial configuration. The one--body density matrix in this latter regime shows a non decaying large distance tail: dislocations, that in 2D are point defects, turn out to be mobile, their number is fluctuating, and they are able to induce exchanges of particles across the system mainly triggered by the dislocation cores. These results indicate that the notion of incommensurate versus commensurate state loses meaning for solid 4He in 2D, because the number of lattice sites becomes ill defined when the system is not commensurate. Crystalline order is found to be stable also in 3D in presence of up to 100 vacancies. △ Less

Submitted 18 March, 2010; v1 submitted 21 July, 2009; originally announced July 2009.

Journal ref: J. Phys.: Condens. Matter 22, 145401 (2010)

Ab-initio low-energy dynamics of superfluid and solid 4He

Authors: E. Vitali, M. Rossi, L. Reatto, D. E. Galli

Abstract: We have extracted information about real time dynamics of 4He 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… ▽ More We have extracted information about real time dynamics of 4He 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 4He we have recovered from exact QMC simulations sharp quasi-particle excitations with spectral functions displaying also the multiphonon branch. As further applications, we have studied the impuriton branch of one 3He atom in liquid 4He and the vacancy-wave excitations in hcp solid 4He finding a novel roton like feature. △ Less

Submitted 4 November, 2010; v1 submitted 27 May, 2009; originally announced May 2009.

Comments: 15 pages, 15 figures; Accepted for publication: Phys.Rev. B 82, (2010)

Solid 4He and the Supersolid Phase: from Theoretical Speculation to the Discovery of a New State of Matter? A Review of the Past and Present Status of Research

Authors: D. E. Galli, L. Reatto

Abstract: The possibility of a supersolid state of matter, i.e., a crystalline solid exhibiting superfluid properties, first appeared in theoretical studies about forty years ago. After a long period of little interest due to the lack of experimental evidence, it has attracted strong experimental and theoretical attention in the last few years since Kim and Chan (Penn State, USA) reported evidence for non… ▽ More The possibility of a supersolid state of matter, i.e., a crystalline solid exhibiting superfluid properties, first appeared in theoretical studies about forty years ago. After a long period of little interest due to the lack of experimental evidence, it has attracted strong experimental and theoretical attention in the last few years since Kim and Chan (Penn State, USA) reported evidence for nonclassical rotational inertia effects, a typical signature of superfluidity, in samples of solid 4He. Since this "first observation", other experimental groups have observed such effects in the response to the rotation of samples of crystalline helium, and it has become clear that the response of the solid is extremely sensitive to growth conditions, annealing processes, and 3He impurities. A peak in the specific heat in the same range of temperatures has been reported as well as anomalies in the elastic behaviour of solid 4He with a strong resemblance to the phenomena revealed by torsional oscillator experiments. Very recently, the observation of unusual mass transport in hcp solid 4He has also been reported, suggesting superflow. From the theoretical point of view, powerful simulation methods have been used to study solid 4He, but the interpretation of the data is still rather difficult; dealing with the question of supersolidity means that one has to face not only the problem of the coexistence of quantum coherence phenomena and crystalline order, exploring the realm of spontaneous symmetry breaking and quantum field theory, but also the problem of the role of disorder, i.e., how defects, such as vacancies, impurities, dislocations, and grain boundaries, participate in the phase transition mechanism. △ Less

Submitted 21 November, 2008; originally announced November 2008.

Comments: Published on J. Phys. Soc. Jpn., Vol.77, No.11, p.111010

Zero-temperature study of vacancies in solid 4He

Authors: M. Rossi, E. Vitali, D. E. Galli, L. Reatto

Abstract: We are investigating the properties of vacancies in solid 4He with the exact zero temperature SPIGS method. Our aim is to study the possibility of phase--separation between vacancies and the perfect crystal. We find a significant correlation among vacancies both in two and in three dimensional solid 4He. In 3D we have considered up to 5 vacancies in a hcp crystal at a density rho = 0.0293 AA^-3,… ▽ More We are investigating the properties of vacancies in solid 4He with the exact zero temperature SPIGS method. Our aim is to study the possibility of phase--separation between vacancies and the perfect crystal. We find a significant correlation among vacancies both in two and in three dimensional solid 4He. In 3D we have considered up to 5 vacancies in a hcp crystal at a density rho = 0.0293 AA^-3, close to the melting, and in 2D we have considered up to 6 vacancies in a triangular crystal at rho = 0.0765 AA^-2. In all the considered cases the correlation among vacancies seems to display an exponential decay that would suggest the presence of a bound state. The systematically sublinear dependence on the number of vacancies found for the activation energy strengthens the indication of an attractive interaction, even if we have direct evidence that vacancies do not form a single compact cluster. △ Less

Submitted 13 June, 2008; originally announced June 2008.

Comments: 4 pages, LT25 Conference proceedings

Zero-point vacancies in quantum solids

Authors: M. Rossi, E. Vitali, D. E. Galli, L. Reatto

Abstract: A Jastrow wave function (JWF) and a shadow wave function (SWF) describe a quantum solid with Bose--Einstein condensate; i.e. a supersolid. It is known that both JWF and SWF describe a quantum solid with also a finite equilibrium concentration of vacancies x_v. We outline a route for estimating x_v by exploiting the existing formal equivalence between the absolute square of the ground state wave… ▽ More A Jastrow wave function (JWF) and a shadow wave function (SWF) describe a quantum solid with Bose--Einstein condensate; i.e. a supersolid. It is known that both JWF and SWF describe a quantum solid with also a finite equilibrium concentration of vacancies x_v. We outline a route for estimating x_v by exploiting the existing formal equivalence between the absolute square of the ground state wave function and the Boltzmann weight of a classical solid. We compute x_v for the quantum solids described by JWF and SWF employing very accurate numerical techniques. For JWF we find a very small value for the zero point vacancy concentration, x_v=(1.4\pm0.1) x 10^-6. For SWF, which presently gives the best variational description of solid 4He, we find the significantly larger value x_v=(1.4\pm0.1) x 10^-3 at a density close to melting. We also study two and three vacancies. We find that there is a strong short range attraction but the vacancies do not form a bound state. △ Less

Submitted 13 June, 2008; originally announced June 2008.

Comments: 19 pages, submitted to J. Low Temp. Phys

Journal ref: J. Low Temp. Phys. 153 (2008) : 250-265

Path Integral Ground State study of 2D solid 4He

Authors: E. Vitali, M. Rossi, F. Tramonto, D. E. Galli, L. Reatto

Abstract: We have studied a two-dimensional triangular commensurate crystal of 4He with the exact T=0 K Path Integral Ground State (PIGS) Monte Carlo method. We have projected onto the true ground state both a Jastrow-Nosanow wave function, in which equilibrium positions are explicitly given and no Bose-Einstein (BEC) is present, and a translationally invariant shadow wave function, in which the solid pha… ▽ More We have studied a two-dimensional triangular commensurate crystal of 4He with the exact T=0 K Path Integral Ground State (PIGS) Monte Carlo method. We have projected onto the true ground state both a Jastrow-Nosanow wave function, in which equilibrium positions are explicitly given and no Bose-Einstein (BEC) is present, and a translationally invariant shadow wave function, in which the solid phase emerges through a spontaneously broken symmetry process and it has BEC. We find a remarkable convergence to the same properties, both the diagonal ones as well as the off-diagonal one-body density matrix rho_1. This supplies a strong evidence that no variational bias are present in the PIGS method. We find no BEC in the commensurate 2D 4He crystal at T=0 K, rho_1 shows an exponential decay in the large distance range. The structure found in rho_1 is due to virtual vacancy--interstitial pairs and this shows up in the momentum distribution. △ Less

Submitted 1 February, 2008; v1 submitted 11 January, 2008; originally announced January 2008.

Comments: 4 pages, 4 figures

Microscopic studies of solid 4He with Path Integral Projector Monte Carlo

Authors: M. Rossi, R. Rota, E. Vitali, D. E. Galli, L. Reatto

Abstract: We have investigated the ground state properties of solid $^4$He with the Shadow Path Integral Ground State method. This exact T=0 K projector method allows to describes quantum solids without introducing any a priori equilibrium position. We have found that the efficiency in computing off-diagonal properties in the solid phase sensibly improves when the direct sampling of permutations, in princ… ▽ More We have investigated the ground state properties of solid $^4$He with the Shadow Path Integral Ground State method. This exact T=0 K projector method allows to describes quantum solids without introducing any a priori equilibrium position. We have found that the efficiency in computing off-diagonal properties in the solid phase sensibly improves when the direct sampling of permutations, in principle not required, is introduced. We have computed the exact one-body density matrix (obdm) in large commensurate 4He crystal finding a decreasing condensate fraction with increasing imaginary time of projection, making our result not conclusive on the presence of Bose-Einstein condensation in bulk solid 4He. We can only give an upper bound of 2.5 times 10^-8 on the condensate fraction. We have exploited the SPIGS method to study also 4He crystal containing grain boundaries by computing the related surface energy and the obdm along these defects. We have found that also highly symmetrical grain boundaries have a finite condensate fraction. We have also derived a route for the estimation of the true equilibrium concentration of vacancies x_v in bulk T=0 K solid 4He, which is shown to be finite, x_v=0.0014(1) at the melting density, when computed with the variational shadow wave function technique. △ Less

Submitted 27 July, 2007; originally announced July 2007.

Comments: Submitted to Proceedings of the 14th International Conference on Recent Progress in Many Body Theories (RPMBT14)

Bounds for the Superfluid Fraction from Exact Quantum Monte Carlo Local Densities

Authors: D. E. Galli, L. Reatto, W. M. Saslow

Abstract: For solid 4He and solid p-H2, using the flow-energy-minimizing one-body phase function and exact T=0 K Monte Carlo calculations of the local density, we have calculated the phase function, the velocity profile and upper bounds for the superfluid fraction f_s. At the melting pressure for solid 4He we find that f_s < 0.20-0.21, about ten times what is observed. This strongly indicates that the the… ▽ More For solid 4He and solid p-H2, using the flow-energy-minimizing one-body phase function and exact T=0 K Monte Carlo calculations of the local density, we have calculated the phase function, the velocity profile and upper bounds for the superfluid fraction f_s. At the melting pressure for solid 4He we find that f_s < 0.20-0.21, about ten times what is observed. This strongly indicates that the theory for the calculation of these upper bounds needs substantial improvements. △ Less

Submitted 27 July, 2007; originally announced July 2007.

Comments: to be published in Phys. Rev. B (Brief Reports)

Two-body correlations and the superfluid fraction for nonuniform systems

Authors: Wayne M. Saslow, Davide E. Galli, Luciano Reatto

Abstract: We extend the one-body phase function upper bound on the superfluid fraction in a periodic solid (a spatially ordered supersolid) to include two-body phase correlations. The one-body current density is no longer proportional to the gradient of the one-body phase times the one-body density, but rather it depends also on two-body correlation functions. The equations that simultaneously determine t… ▽ More We extend the one-body phase function upper bound on the superfluid fraction in a periodic solid (a spatially ordered supersolid) to include two-body phase correlations. The one-body current density is no longer proportional to the gradient of the one-body phase times the one-body density, but rather it depends also on two-body correlation functions. The equations that simultaneously determine the one-body and two-body phase functions require a knowledge of one-, two-, and three-body correlation functions. The approach can also be extended to disordered solids. Fluids, with two-body densities and two-body phase functions that are translationally invariant, cannot take advantage of this additional degree of freedom to lower their energy. △ Less

Submitted 18 December, 2006; originally announced December 2006.

Comments: 13 pages