1
|
Colognesi D, Bafile U, Guarini E, Neumann M. Microscopic self-dynamics in liquid Ne-D_{2} mixtures: Quantum features and itinerant oscillators reexamined. Phys Rev E 2024; 109:045104. [PMID: 38755808 DOI: 10.1103/physreve.109.045104] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2023] [Accepted: 03/20/2024] [Indexed: 05/18/2024]
Abstract
In this paper, we report the results of a centroid molecular dynamics (CMD) study of the canonical velocity autocorrelation functions (VACFs) in liquid Ne-D_{2} mixtures at a temperature of T=30K and in the full D_{2}-concentration range (0%≤x_{D_{2}}≤100%). This binary system was selected because of its moderate, although sizable, quantum effects which, as far as its equilibrium properties are concerned, are fully described by the path integral Monte Carlo (PIMC) simulations that have been also implemented. A comprehensive test of the VACF spectral moments carried out using three physical quantities (namely, mean kinetic energy, Einstein frequency, and mean-squared force) obtained from PIMC was performed revealing the potentialities, as well as the limitations, of the CMD approach to the single-particle dynamics in these low-T liquid mixtures. Additional physical information was extracted from the canonical VACFs by fitting their spectra via two distinct methods: the Levesque-Verlet model (LV, very flexible but highly heuristic) and the itinerant oscillator model (IO, based on the physical ground of a single particle rattling inside a short-lived diffusing pseudocage). Both provided good fits of the CMD outputs, with LV being always more adequate than IO in the case of the Ne VACFs, while, as for the D_{2} VACFs, the LV superiority is evident only at high x_{D_{2}} values. However, a peculiar and systematic effect was found after analyzing the IO-fitted parameters: the estimated pseudocage masses turned out to be at least one order of magnitude lower than the corresponding values inferred from the PIMC simulations. This outcome concerns both the Ne and the D_{2} rattling molecule and, as we also discovered, had already been observed (but promptly forgotten) in purely classical simulations of liquid Ar. The possible physical origins of this finding have been finally discussed in some detail, also in connection with the result of the more recent exponential expansion theory (EET), which manages to shed more light on the concept of single particles rattling inside short-lived pseudocages, ultimately demonstrating its untenability.
Collapse
Affiliation(s)
- Daniele Colognesi
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata "N. Carrara", via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| | - Ubaldo Bafile
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata "N. Carrara", via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| | - Eleonora Guarini
- Dipartimento di Fisica e Astronomia, Università di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
| | - Martin Neumann
- Fakultät für Physik der Universität Wien, Kolingasse 14-16, A-1090 Wien, Austria
| |
Collapse
|
2
|
Barocchi F, Guarini E. Universal Matsubara time decay of quantum autocorrelations for Boltzmann particles. Phys Rev E 2022; 106:044128. [PMID: 36397547 DOI: 10.1103/physreve.106.044128] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2022] [Accepted: 09/20/2022] [Indexed: 06/16/2023]
Abstract
The general properties of time dependent autocorrelations in many-body quantum systems are here analyzed at thermodynamic equilibrium in the Boltzmann canonical ensemble at temperature T, by means of the exponential expansion theory (EET). It is shown that the Kubo-Martin-Schwinger (KMS) symmetry applied to the exponential expansion of the correlation leads to the existence of two different sets of decay modes (channels) here indicated as "Matsubara modes" and "system modes," respectively. The Matsubara modes are a series of pure decay channels with time constants representing a direct action of the thermostat upon the correlation, with a characteristic principal decay time τ_{1}=ℏ/(2πk_{B}T), where ℏ and k_{B} are the Planck and Boltzmann constants, and T is the temperature. Moreover, the KMS condition implies that the amplitudes pertaining to the even and odd contribution of the system modes to the quantum correlation are not independent. These two properties are quantum mechanical in nature and "universal," in the sense that they are present for any autocorrelation of a quantum system at equilibrium at a temperature T. The Matsubara modes' contribution to the time behavior of a quantum correlation is limited to times of the order of τ_{1}, which however can be comparable with some of the characteristic decay times of the system modes. In addition, since the parameters representing the overall time behavior of the quantum correlation can be given in terms of the parameters of its Kubo transform, the EET representation turns out to be useful in calculations exploiting the outputs of some widespread quantum simulation methods. A discussion of the properties of these relations is described in detail with numerical examples. The case of the velocity autocorrelation function of para hydrogen at low temperature is also reported as a final example for a real system.
Collapse
Affiliation(s)
- Fabrizio Barocchi
- Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
| | - Eleonora Guarini
- Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
| |
Collapse
|
3
|
De Francesco A, Scaccia L, Formisano F, Guarini E, Bafile U, Maccarini M, Zhang Y, Nykypanchuck D, Alatas A, Cunsolo A. The damping of terahertz acoustic modes in aqueous nanoparticle suspensions. Sci Rep 2021; 11:20110. [PMID: 34635734 PMCID: PMC8505432 DOI: 10.1038/s41598-021-99503-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2021] [Accepted: 09/20/2021] [Indexed: 11/23/2022] Open
Abstract
In this work, we investigate the possibility of controlling the acoustic damping in a liquid when nanoparticles are suspended in it. To shed light on this topic, we performed Inelastic X-Ray Scattering (IXS) measurements of the terahertz collective dynamics of aqueous suspensions of nanospheres of various materials, size, and relative concentration, either charged or neutral. A Bayesian analysis of measured spectra indicates that the damping of the two acoustic modes of water increases upon nanoparticle immersion. This effect seems particularly pronounced for the longitudinal acoustic mode, which, whenever visible at all, rapidly damps off when increasing the exchanged wavevector. Results also indicate that the observed effect strongly depends on the material the immersed nanoparticles are made of.
Collapse
Affiliation(s)
- Alessio De Francesco
- CNR-IOM & INSIDE@ILL c/o Operative Group in Grenoble (OGG), 38042, Grenoble, France
- Institut Laue-Langevin (ILL), 38042, Grenoble, France
| | - Luisa Scaccia
- Dipartimento di Economia e Diritto, Università di Macerata, Via Crescimbeni 20, 62100, Macerata, Italy
| | - Ferdinando Formisano
- CNR-IOM & INSIDE@ILL c/o Operative Group in Grenoble (OGG), 38042, Grenoble, France
- Institut Laue-Langevin (ILL), 38042, Grenoble, France
| | - Eleonora Guarini
- Dipartimento di Fisica e Astronomia, Università di Firenze, via G. Sansone 1, 50019, Sesto Fiorentino, Italy
| | - Ubaldo Bafile
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata "Nello Carrara", via Madonna del Piano 10, 50019, Sesto Fiorentino, Italy
| | - Marco Maccarini
- Laboratoire TIMC/IMAG UMR CNRS 5525, Université Grenoble-Alpes, 38000, Grenoble, France
| | - Yugang Zhang
- Brookhaven National Laboratory, Center for Functional Nanomaterials, P.O. Box 5000, Upton, 11973, NY, USA
| | - Dmytro Nykypanchuck
- Brookhaven National Laboratory, Center for Functional Nanomaterials, P.O. Box 5000, Upton, 11973, NY, USA
| | - Ahmet Alatas
- Argonne National Laboratory, Advanced Photon Source, P.O. Box 5000, Upton, 11973, NY, USA
| | - Alessandro Cunsolo
- Department of Physics, University of Wisconsin at Madison, 1150 University Avenue, Madison, WI, USA.
| |
Collapse
|
4
|
De Francesco A, Bafile U, Cunsolo A, Scaccia L, Guarini E. Searching for a second excitation in the inelastic neutron scattering spectrum of a liquid metal: a Bayesian analysis. Sci Rep 2021; 11:13974. [PMID: 34234260 PMCID: PMC8263748 DOI: 10.1038/s41598-021-93452-w] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2021] [Accepted: 06/24/2021] [Indexed: 12/04/2022] Open
Abstract
When probed at nanometer and picosecond scales, the properties of a liquid present striking analogies with the ones of the corresponding solid, one of the most surprising is the ability of supporting shear wave propagation, as a rigid medium. Although this evidence is being reported by a growing number of terahertz scattering measurements, it remains an open question whether it is universal or rather typical of some liquids only. Furthermore, given its elusive signatures in the scattering signal, the detection of this effect appears as a typical case where an unintentional “bias of confirmation” can mislead experimentalists. We thus decided to use a Bayesian inference approach to achieve a probabilistically grounded and evidence-based lineshape modeling of the inelastic neutron scattering spectra from liquid silver, whose simulated density autocorrelations bear evidence of a shear mode propagation over very short distances. The result of our analysis indicates that the observation of any additional, non-longitudinal, acoustic modes in this simple system goes beyond the accuracy of the used scattering method.
Collapse
Affiliation(s)
- Alessio De Francesco
- CNR-IOM & INSIDE@ILL c/o Operative Group in Grenoble (OGG), F-38042, Grenoble, France.,Institut Laue-Langevin (ILL), F-38042, Grenoble, France
| | - Ubaldo Bafile
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata "Nello Carrara", Via Madonna del Piano 10, I-50019, Sesto Fiorentino, Italy
| | - Alessandro Cunsolo
- Department of Physics, University of Wisconsin at Madison, 1150 University Avenue, Madison, WI, USA
| | - Luisa Scaccia
- Dipartimento di Economia e Diritto, Università di Macerata, Via Crescimbeni 20, I-62100, Macerata, Italy
| | - Eleonora Guarini
- Dipartimento di Fisica e Astronomia, Università di Firenze, Via G. Sansone 1, I-50019, Sesto Fiorentino, Italy.
| |
Collapse
|
5
|
Bafile U, Neumann M, Colognesi D, Guarini E. Time dependence of quantum correlation functions. Phys Rev E 2020; 101:052110. [PMID: 32575332 DOI: 10.1103/physreve.101.052110] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2020] [Accepted: 04/21/2020] [Indexed: 01/02/2023]
Abstract
In the past few years, the exponential expansion analysis of time autocorrelation functions has provided profound insight into the leading microscopic processes driving the atomic-scale dynamics and has made it possible to highlight the presence and the role of various relaxation channels through which the fundamental correlation functions decay with time. Here we apply this method to the determination of the full time dependence of a correlation function c(t) in a quantum system at nonzero temperature, by making explicit its relationship with its Kubo transform c_{K}(t), which in some cases can be approximately computed with the presently available quantum simulation techniques. We obtain an exact expression for c(t) in terms of the exponential modes that describe the time behavior of c_{K}(t). The relative importance of the various modes in determining the overall shape of c(t) can then be studied in detail. This work extends to the full time domain the results of a previous paper [Guarini et al., Phys. Rev. Lett. 123, 135301 (2019)PRLTAO0031-900710.1103/PhysRevLett.123.135301], in which we employed the same method to calculate the zero time value of the velocity autocorrelation function, to obtain a microscopic description of the quantum mean kinetic energy in a fluid. In particular, we show that the decay constants and the frequencies of the dominant microscopic modes of c(t) are the same as those of c_{K}(t), but the dynamics of the quantum system also contains an additional term decaying on a time scale determined solely by temperature of the system.
Collapse
Affiliation(s)
- Ubaldo Bafile
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata "Nello Carrara", via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| | - Martin Neumann
- Fakultät für Physik der Universität Wien, Strudlhofgasse 4, A-1090 Wien, Austria
| | - Daniele Colognesi
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata "Nello Carrara", via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| | - Eleonora Guarini
- Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
| |
Collapse
|