-
Quantum effects and temperature dependence of the free energy barrier and shear frequency in bilayer graphene
Authors:
Jean Paul Nery,
Lorenzo Monacelli,
Francesco Mauri
Abstract:
Little is known about the quantum effects of friction, and first-principles studies are virtually non-existent. In multilayer graphene, the order of the stacking of the layers plays a crucial role in the electronic properties and the manifestation of superconductivity, and friction affects how layers slide to change the stacking. Here, we study the effects of friction in bilayer graphene by focusi…
▽ More
Little is known about the quantum effects of friction, and first-principles studies are virtually non-existent. In multilayer graphene, the order of the stacking of the layers plays a crucial role in the electronic properties and the manifestation of superconductivity, and friction affects how layers slide to change the stacking. Here, we study the effects of friction in bilayer graphene by focusing on quantum fluctuations in the free energy barrier that layers have to overcome to slide with respect to each other. We evaluate the free energy at the barrier in bilayer graphene through the self-consistent harmonic approximation, accounting for quantum and thermal ionic fluctuations. We obtain a very large reduction of the barrier of more than 30%. Temperature effects are weaker, although high temperatures induce a significant reduction of the barrier as well. Our approach is general and paves the way for systematically accounting for nuclear quantum and thermal effects in free energy barriers of other macroscopic systems.
△ Less
Submitted 5 July, 2024;
originally announced July 2024.
-
Beyond Gaussian fluctuations of quantum anharmonic nuclei
Authors:
Antonio Siciliano,
Lorenzo Monacelli,
Francesco Mauri
Abstract:
The Self-Consistent Harmonic Approximation (SCHA) describes atoms in solids, including quantum fluctuations and anharmonic effects, in a non-perturbative way. It computes ionic free energy variationally, constraining the atomic quantum-thermal fluctuations to be Gaussian. Consequently, the entropy is analytical; there is no need for thermodynamic integration or heavy diagonalization to include fin…
▽ More
The Self-Consistent Harmonic Approximation (SCHA) describes atoms in solids, including quantum fluctuations and anharmonic effects, in a non-perturbative way. It computes ionic free energy variationally, constraining the atomic quantum-thermal fluctuations to be Gaussian. Consequently, the entropy is analytical; there is no need for thermodynamic integration or heavy diagonalization to include finite temperature effects. In addition, as the probability distribution is fixed, SCHA solves all the equations with Monte Carlo integration without employing Metropolis sampling of the quantum phase space. Unfortunately, the Gaussian approximation breaks down for rotational modes and tunneling effects. We show how to describe these non-Gaussian fluctuations using the quantum variational principle at finite temperatures, keeping the main advantage of SCHA: direct access to free energy. Our method, nonlinear SCHA (NLSCHA), employs an invertible nonlinear transformation to map Cartesian coordinates into an auxiliary manifold parametrized by a finite set of variables. So, we adopt a Gaussian \textit{ansatz} for the density matrix in this new coordinate system. The nonlinearity of the mapping ensures that NLSCHA enlarges the SCHA variational subspace, and its invertibility conserves the information encoded in the density matrix. We evaluate the entropy in the auxiliary space, where it has a simple analytical form. As in the SCHA, the variational principle allows for optimizing free parameters to minimize free energy. Finally, we show that, for the first time, NLSCHA gives direct access to the entropy of a crystal with non-Gaussian degrees of freedom.
△ Less
Submitted 4 July, 2024;
originally announced July 2024.
-
Simulating anharmonic crystals: Lights and shadows of first-principles approaches
Authors:
Lorenzo Monacelli
Abstract:
Understanding and simulating the thermodynamic and dynamical properties of materials affected by strong ionic anharmonicity is a central challenge in material science. Two powerful methodologies have emerged as frontrunners in this exploration: the Self-Consistent Harmonic Approximation (SCHA) with its counterpart, Self-Consistent Phonons (SCP), and the Temperature-Dependent Effective Potential (T…
▽ More
Understanding and simulating the thermodynamic and dynamical properties of materials affected by strong ionic anharmonicity is a central challenge in material science. Two powerful methodologies have emerged as frontrunners in this exploration: the Self-Consistent Harmonic Approximation (SCHA) with its counterpart, Self-Consistent Phonons (SCP), and the Temperature-Dependent Effective Potential (TDEP). Despite their widespread use, the limitations of these methods are often overlooked, and the absence of a formal derivation for the TDEP has led to significant missteps that hamper its predictive accuracy.
Here, the theoretical foundations of both SCHA and TDEP are revised, unveiling their profound interconnection. This work introduces formal derivation for TDEP, dispelling misconceptions about how to improve its accuracy and refining best practices. The perturbative limit of these methods uncovers that both SCHA and TDEP fall short in reproducing the lowest-order perturbative free energy. While TDEP excels in providing an exact treatment of the static susceptibility, its dynamical extension violates the perturbative regime by overcounting anharmonicity. Conversely, SCHA delivers an approximate static susceptibility but maintains correct static and dynamic response functions within the perturbative limit.
This work introduces a corrective strategy for TDEP's dynamical extension to address these shortcomings. All the claims provided here are benchmarked against the exact (numerical) solution on a simple one-dimensional anharmonic potential.
△ Less
Submitted 3 July, 2024;
originally announced July 2024.
-
Giant splitting of the hydrogen rotational eigenenergies in the C$_2$ filled ice
Authors:
Simone Di Cataldo,
Maria Rescigno,
Lorenzo Monacelli,
Umbertoluca Ranieri,
Richard Gaal,
Stefan Klotz,
Jacques Ollivier,
Michael Marek Koza,
Cristiano De Michele,
Livia Eleonora Bove
Abstract:
Hydrogen hydrates present a rich phase diagram influenced by both pressure and temperature, with the so-called C$_2$ phase emerging prominently above 2.5 GPa. In this phase, hydrogen molecules are densely packed within a cubic ice-like lattice and the interaction with the surrounding water molecules profoundly affects their quantum rotational dynamics. Herein, we delve into this intricate interpla…
▽ More
Hydrogen hydrates present a rich phase diagram influenced by both pressure and temperature, with the so-called C$_2$ phase emerging prominently above 2.5 GPa. In this phase, hydrogen molecules are densely packed within a cubic ice-like lattice and the interaction with the surrounding water molecules profoundly affects their quantum rotational dynamics. Herein, we delve into this intricate interplay by directly solving the Schrödinger's equation for a quantum H$_2$ rotor in the C$_2$ crystal field at finite temperature, generated through Density Functional Theory. Our calculations reveal a giant energy splitting relative to the magnetic quantum number of $\pm$3.2 meV for $l=1$. Employing inelastic neutron scattering, we experimentally measure the energy levels of H$_2$ within the C$_2$ phase at 6.0 and 3.4 GPa and low temperatures, finding remarkable agreement with our theoretical predictions. These findings underscore the pivotal role of hydrogen--water interactions in dictating the rotational behavior of the hydrogen molecules within the C$_2$ phase and indicate heightened induced-dipole interactions compared to other hydrogen hydrates.
△ Less
Submitted 23 May, 2024;
originally announced May 2024.
-
Quantum effects in the H-bond symmetrization and in the thermodynamic properties of high pressure ice
Authors:
Marco Cherubini,
Lorenzo Monacelli,
Bingjia Yang,
Roberto Car,
Michele Casula,
Francesco Mauri
Abstract:
We investigate the structural and thermodynamic properties of high-pressure ice by incorporating quantum anharmonicity at a non-perturbative level. Quantum fluctuations reduce the critical pressure of the phase transition between phase VIII (with asymmetric H-bonds) and phase X (with symmetric H-bonds) by 65 GPa from its classical value of 116 GPa at 0K. Moreover, quantum effects make it temperatu…
▽ More
We investigate the structural and thermodynamic properties of high-pressure ice by incorporating quantum anharmonicity at a non-perturbative level. Quantum fluctuations reduce the critical pressure of the phase transition between phase VIII (with asymmetric H-bonds) and phase X (with symmetric H-bonds) by 65 GPa from its classical value of 116 GPa at 0K. Moreover, quantum effects make it temperature-independent over a wide temperature range (0K-300K), in agreement with experimental estimates obtained through vibrational spectroscopy and in striking contrast to the strong temperature dependence found in the classical approximation. The equation of state shows fingerprints of the transition in accordance with experimental evidence. Additionally, we demonstrate that, within our approach, proton disorder in phase VII has a negligible impact on the occurrence of phase X. Finally, we reproduce with high accuracy the 10 GPa isotope shift due to the hydrogen-to-deuterium substitution.
△ Less
Submitted 14 March, 2024;
originally announced March 2024.
-
Electronic structure and lattice dynamics of 1T-VSe$_2$: origin of the 3D-CDW
Authors:
Josu Diego,
D. Subires,
A. H. Said,
D. A. Chaney,
A. Korshunov,
G. Garbarino,
F. Diekmann,
K. Mahatha,
V. Pardo,
J. Strempfer,
Pablo J. Bereciartua Perez,
S. Francoual,
C. Popescu,
M. Tallarida,
J. Dai,
Raffaello Bianco,
Lorenzo Monacelli,
Matteo Calandra,
A. Bosak,
Francesco Mauri,
K. Rossnagel,
Adolfo O. Fumega,
Ion Errea,
S. Blanco-Canosa
Abstract:
In order to characterize in detail the charge density wave (CDW) transition of 1$T$-VSe$_2$, its electronic structure and lattice dynamics are comprehensively studied by means of x-ray diffraction, angle resolved photoemission (ARPES), diffuse and inelastic x-ray scattering (IXS), and state-of-the-art first principles density functional theory calculations. Resonant elastic x-ray scattering (REXS)…
▽ More
In order to characterize in detail the charge density wave (CDW) transition of 1$T$-VSe$_2$, its electronic structure and lattice dynamics are comprehensively studied by means of x-ray diffraction, angle resolved photoemission (ARPES), diffuse and inelastic x-ray scattering (IXS), and state-of-the-art first principles density functional theory calculations. Resonant elastic x-ray scattering (REXS) does not show any resonant enhancement at either V or Se K-edges, indicating that the CDW peak describes a purely structural modulation of the electronic ordering. ARPES identifies (i) a pseudogap at T$>$T$_{CDW}$, which leads to a depletion of the density of states in the $ML-M'L'$ plane at T$<$T$_{CDW}$, and (ii) anomalies in the electronic dispersion reflecting a sizable impact of phonons on it. A diffuse scattering precursor, characteristic of soft phonons, is observed at room temperature (RT) and leads to the full collapse of the low-energy phonon ($ω_1$) with propagation vector (0.25 0 -0.3) r.l.u. We show that the frequency and linewidth of this mode are anisotropic in momentum space, reflecting the momentum dependence of the electron-phonon interaction (EPI), hence demonstrating that the origin of the CDW is, to a much larger extent, due to the momentum dependence EPI with a small contribution from nesting. The pressure dependence of the $ω_1$ soft mode remains nearly constant up to 13 GPa at RT, with only a modest softening before the transition to the high-pressure monoclinic $C2/m$ phase. The wide set of experimental data are well captured by our state-of-the art first-principles anharmonic calculations with the inclusion of van der Waals (vdW) corrections in the exchange-correlation functional. The description of the electronics and dynamics of VSe$_2$ reported here adds important pieces of information to the understanding of the electronic modulations of TMDs.
△ Less
Submitted 28 July, 2023;
originally announced July 2023.
-
Magnetostriction-Driven Muon Localization in an Antiferromagnetic Oxide
Authors:
Pietro Bonfà,
Ifeanyi John Onuorah,
Franz Lang,
Iurii Timrov,
Lorenzo Monacelli,
Chennan Wang,
Xiao Sun,
Oleg Petracic,
Giovanni Pizzi,
Nicola Marzari,
Stephen J. Blundell,
Roberto De Renzi
Abstract:
Magnetostriction drives a rhombohedral distortion in the cubic rock salt antiferromagnet MnO at the Néel temperature $T_{N}=118$ K. As an unexpected consequence we show that this distortion acts to localize the site of an implanted muon due to the accompanying redistribution of electron density. This lifts the degeneracy between equivalent sites, resulting in a single observed muon precession freq…
▽ More
Magnetostriction drives a rhombohedral distortion in the cubic rock salt antiferromagnet MnO at the Néel temperature $T_{N}=118$ K. As an unexpected consequence we show that this distortion acts to localize the site of an implanted muon due to the accompanying redistribution of electron density. This lifts the degeneracy between equivalent sites, resulting in a single observed muon precession frequency. Above $T_{N}$, the muon instead becomes delocalized around a network of equivalent sites. Our first-principles simulations based on Hubbard-corrected density-functional theory and molecular dynamics are consistent with our experimental data and help to resolve a long-standing puzzle regarding muon data on MnO, as well as having wider applicability to other magnetic oxides.
△ Less
Submitted 28 January, 2024; v1 submitted 20 May, 2023;
originally announced May 2023.
-
Large impact of phonon lineshapes on the superconductivity of solid hydrogen
Authors:
Đorđe Dangić,
Lorenzo Monacelli,
Raffaello Bianco,
Francesco Mauri,
Ion Errea
Abstract:
Phonon anharmonicity plays a crucial role in determining the stability and vibrational properties of high-pressure hydrides. Furthermore, strong anharmonicity can render phonon quasiparticle picture obsolete questioning standard approaches for modeling superconductivity in these material systems. In this work, we show the effects of non-Lorentzian phonon lineshapes on the superconductivity of high…
▽ More
Phonon anharmonicity plays a crucial role in determining the stability and vibrational properties of high-pressure hydrides. Furthermore, strong anharmonicity can render phonon quasiparticle picture obsolete questioning standard approaches for modeling superconductivity in these material systems. In this work, we show the effects of non-Lorentzian phonon lineshapes on the superconductivity of high-pressure solid hydrogen. We calculate the superconducting critical temperature T$_\mathrm{C}$ \emph{ab initio} considering the full phonon spectral function and show that it overall enhances the T$_\mathrm{C}$ estimate. The anharmonicity-induced phonon softening exhibited in spectral functions increases the estimate of the critical temperature, while the broadening of phonon lines due to phonon-phonon interaction decreases it. Our calculations also reveal that superconductivity emerges in hydrogen in the $Cmca-12$ molecular phase VI at pressures between 450 and 500 GPa and explain the disagreement between the previous theoretical results and experiments.
△ Less
Submitted 22 March, 2024; v1 submitted 14 March, 2023;
originally announced March 2023.
-
First-principles thermodynamics of CsSnI3
Authors:
Lorenzo Monacelli,
Nicola Marzari
Abstract:
CsSnI3 is a promising eco-friendly solution for energy harvesting technologies. It exists at room temperature in either a black perovskite polymorph or a yellow 1D double-chain, which irreversibly deteriorates in the air. In this work, we unveil the relative thermodynamic stability between the two structures with a first-principles sampling of the CsSnI3 finite-temperature phase diagram, discoveri…
▽ More
CsSnI3 is a promising eco-friendly solution for energy harvesting technologies. It exists at room temperature in either a black perovskite polymorph or a yellow 1D double-chain, which irreversibly deteriorates in the air. In this work, we unveil the relative thermodynamic stability between the two structures with a first-principles sampling of the CsSnI3 finite-temperature phase diagram, discovering how it is driven by anomalously large quantum and anharmonic ionic fluctuations. Thanks to a comprehensive treatment of anharmonicity, the simulations deliver a remarkable agreement with known experimental data for the transition temperatures of the orthorhombic, rhombohedral, and cubic perovskite structures and the thermal expansion coefficient. We disclose how the perovskite polymorphs are the ground state above \SI{270}{\kelvin} and discover an abnormal decrease in heat capacity upon heating in the cubic black perovskite. Our results also significantly downplay the Cs+ rattling modes' contribution to mechanical instability. The remarkable agreement with experiments validates our methodology, which can be systematically applied to all metal halides.
△ Less
Submitted 24 January, 2023;
originally announced January 2023.
-
Wigner Gaussian dynamics: simulating the anharmonic and quantum ionic motion
Authors:
Antonio Siciliano,
Lorenzo Monacelli,
Giovanni Caldarelli,
Francesco Mauri
Abstract:
The atomic motion controls important features of materials, such as thermal transport, phase transitions, and vibrational spectra. However, the simulation of ionic dynamics is exceptionally challenging when quantum fluctuations are relevant (e.g., at low temperatures or with light atoms) and the energy landscape is anharmonic. In this work, we formulate the Time-Dependent Self-Consistent Harmonic…
▽ More
The atomic motion controls important features of materials, such as thermal transport, phase transitions, and vibrational spectra. However, the simulation of ionic dynamics is exceptionally challenging when quantum fluctuations are relevant (e.g., at low temperatures or with light atoms) and the energy landscape is anharmonic. In this work, we formulate the Time-Dependent Self-Consistent Harmonic Approximation (TDSCHA) in the Wigner framework, paving the way for the efficient computation of the nuclear motion in systems with sizable quantum and thermal anharmonic fluctuations. Besides the improved numerical efficiency, the Wigner formalism unveils the classical limit of TDSCHA and provides a link with the many-body perturbation theory of Feynman diagrams.
We further extend the method to account for the non-linear couplings between phonons and photons, responsible, e.g., for a nonvanishing Raman signal in high-symmetry Raman inactive crystals, firstly discussed by Rasetti and Fermi. We benchmark the method in phase III of high-pressure hydrogen ab initio. The nonlinear photon-phonon coupling reshapes the IR spectra and explains the high-frequency shoulder of the H2 vibron observed in experiments.
The Wigner TDSCHA is computationally cheap and derived from first principles: it is unbiased by assumptions on the phonon-phonon and phonon-photon scattering and does not depend on empirical parameters. Therefore, the method can be adopted in unsupervised high-throughput calculations.
△ Less
Submitted 24 May, 2023; v1 submitted 20 January, 2023;
originally announced January 2023.
-
Temperature-dependent anharmonic phonons in quantum paraelectric KTaO$_3$ by first principles and machine-learned force fields
Authors:
Luigi Ranalli,
Carla Verdi,
Lorenzo Monacelli,
Matteo Calandra,
Georg Kresse,
Cesare Franchini
Abstract:
Understanding collective phenomena in quantum materials from first principles is a promising route toward engineering materials properties on demand and designing new functionalities. This work examines the quantum paraelectric state, an elusive state of matter characterized by the smooth saturation of the ferroelectric instability at low temperature due to quantum fluctuations associated with anh…
▽ More
Understanding collective phenomena in quantum materials from first principles is a promising route toward engineering materials properties on demand and designing new functionalities. This work examines the quantum paraelectric state, an elusive state of matter characterized by the smooth saturation of the ferroelectric instability at low temperature due to quantum fluctuations associated with anharmonic phonon effects. The temperature-dependent evolution of the soft ferroelectric phonon mode in the quantum paraelectric KTaO$_3$ in the range 0-300 K is modelled by combining density functional theory (DFT) calculations with the stochastic self-consistent harmonic approximation assisted by an on-the-fly machine-learned force field. The calculated data show that including anharmonic terms is essential to stabilize the spurious imaginary ferroelectric phonon predicted by DFT, in agreement with experiments. Augmenting the DFT workflow with machine-learned force fields allows for efficient stochastic sampling of the configurational space using large supercells in a broad and dense temperature range, inaccessible by conventional ab initio protocols. This work proposes a robust computational workflow capable of accounting for collective behaviors involving different degrees of freedom and occurring at large time/length scales, paving the way for precise modeling and control of quantum effects in materials.
△ Less
Submitted 24 September, 2022;
originally announced September 2022.
-
Quantum phase diagram of high-pressure hydrogen
Authors:
Lorenzo Monacelli,
Michele Casula,
Kosuke Nakano,
Sandro Sorella,
Francesco Mauri
Abstract:
The interplay between electron correlation and nuclear quantum effects makes our understanding of elemental hydrogen a formidable challenge. Here, we present the phase diagram of hydrogen and deuterium at low temperatures and high-pressure ($P > 300$ GPa by accounting for highly accurate electronic and nuclear enthalpies. We evaluated internal electronic energies by diffusion quantum Monte Carlo,…
▽ More
The interplay between electron correlation and nuclear quantum effects makes our understanding of elemental hydrogen a formidable challenge. Here, we present the phase diagram of hydrogen and deuterium at low temperatures and high-pressure ($P > 300$ GPa by accounting for highly accurate electronic and nuclear enthalpies. We evaluated internal electronic energies by diffusion quantum Monte Carlo, while nuclear quantum motion and anharmonicity have been included by the stochastic self-consistent harmonic approximation. Our results show that the long-sought atomic metallic hydrogen, predicted to host room-temperature superconductivity, forms at $577\pm 10$ GPa ($640\pm 14$ GPa in deuterium). Indeed, anharmonicity pushes the stability of this phase towards pressures much larger than previous theoretical estimates or attained experimental values. Before atomization, molecular hydrogen transforms from a conductive phase III to another metallic structure that is still molecular (phase VI) at $422\pm 40$ GPa ($442\pm30$ GPa in deuterium). We predict clear-cut signatures in optical spectroscopy and DC conductivity that can be used experimentally to distinguish between the two structural transitions. According to our findings, the experimental evidence of metallic hydrogen has so far been limited to molecular phases.
△ Less
Submitted 14 February, 2022; v1 submitted 11 February, 2022;
originally announced February 2022.
-
Understanding Anharmonic Effects on Hydrogen Desorption Characteristics of Mg$_n$H$_{2n}$ Nanoclusters by ab initio trained Deep Neural Network
Authors:
Andrea Pedrielli,
Paolo E. Trevisanutto,
Lorenzo Monacelli,
Giovanni Garberoglio,
Nicola M. Pugno,
Simone Taioli
Abstract:
Magnesium hydride (MgH$_2$) has been widely studied for effective hydrogen storage. However, its bulk desorption temperature (553 K) is deemed too high for practical applications. Besides doping, a strategy to decrease such reaction energy for releasing hydrogen is the use of MgH$_2$-based nanoparticles (NPs). Here, we investigate first the thermodynamic properties of Mg$_n$H$_{2n}$ NPs ($n<10$) f…
▽ More
Magnesium hydride (MgH$_2$) has been widely studied for effective hydrogen storage. However, its bulk desorption temperature (553 K) is deemed too high for practical applications. Besides doping, a strategy to decrease such reaction energy for releasing hydrogen is the use of MgH$_2$-based nanoparticles (NPs). Here, we investigate first the thermodynamic properties of Mg$_n$H$_{2n}$ NPs ($n<10$) from first-principles, in particular by assessing the anharmonic effects on the enthalpy, entropy and thermal expansion by means of the Stochastic Self Consistent Harmonic Approximation (SSCHA). The latter method goes beyond previous approaches, typically based on molecular mechanics and the quasi-harmonic approximation, allowing the ab initio calculation of the fully-anharmonic free energy. We find an almost linear dependence on temperature of the interatomic bond lengths - with a relative variation of few percent over 300K -, alongside with a bond distance decrease of the Mg-H bonds. In order to increase the size of NPs toward experiments of hydrogen desorption from MgH$_2$ we devise a computationally effective Machine Learning model trained to accurately determine the forces and total energies (i.e. the potential energy surfaces), integrating the latter with the SSCHA model to fully include the anharmonic effects. We find a significative decrease of the H-desorption temperature for sub-nanometric clusters Mg$_n$H$_{2n}$ with $n \leq 10$, with a non-negligible, although little effect due to anharmonicities (up to 10%).
△ Less
Submitted 27 November, 2021;
originally announced November 2021.
-
The microscopic origin of anomalous properties of ice relies on the strong quantum anharmonic regime of atomic vibrations
Authors:
Marco Cherubini,
Lorenzo Monacelli,
Francesco Mauri
Abstract:
Water ice is a unique material presenting intriguing physical properties, like negative thermal expansion and anomalous volume isotope effect (VIE). They arise from the interplay between weak hydrogen bonds and nuclear quantum fluctuations, making theoretical calculations challenging. Here, we employ the stochastic self-consistent harmonic approximation (SSCHA) to investigate how thermal and quant…
▽ More
Water ice is a unique material presenting intriguing physical properties, like negative thermal expansion and anomalous volume isotope effect (VIE). They arise from the interplay between weak hydrogen bonds and nuclear quantum fluctuations, making theoretical calculations challenging. Here, we employ the stochastic self-consistent harmonic approximation (SSCHA) to investigate how thermal and quantum fluctuations affect the physical properties of ice XI ab initio. Regarding the anomalous VIE, our work reveals that quantum effects on hydrogen are so strong to be in a nonlinear regime: when progressively increasing the mass of hydrogen from protium to infinity (classical limit), the volume firstly expands and then contracts, with a maximum slightly above the mass of tritium. We observe an anharmonic renormalization of about 10% in the bending and stretching phonon frequencies probed in IR and Raman experiments. For the first time, we report an accurate comparison of the low energy phonon dispersion with the experimental data, possible only thanks to high-level accuracy in the electronic correlation and nuclear quantum and thermal fluctuations, paving the way for the study of thermal transport in ice from first principles and the simulation of ice under pressure.
△ Less
Submitted 17 February, 2022; v1 submitted 7 July, 2021;
originally announced July 2021.
-
The Stochastic Self-Consistent Harmonic Approximation: Calculating Vibrational Properties of Materials with Full Quantum and Anharmonic Effects
Authors:
Lorenzo Monacelli,
Raffaello Bianco,
Marco Cherubini,
Matteo Calandra,
Ion Errea,
Francesco Mauri
Abstract:
The efficient and accurate calculation of how ionic quantum and thermal fluctuations impact the free energy of a crystal, its atomic structure, and phonon spectrum is one of the main challenges of solid state physics, especially when strong anharmonicy invalidates any perturbative approach. To tackle this problem, we present the implementation on a modular Python code of the stochastic self-consis…
▽ More
The efficient and accurate calculation of how ionic quantum and thermal fluctuations impact the free energy of a crystal, its atomic structure, and phonon spectrum is one of the main challenges of solid state physics, especially when strong anharmonicy invalidates any perturbative approach. To tackle this problem, we present the implementation on a modular Python code of the stochastic self-consistent harmonic approximation method. This technique rigorously describes the full thermodyamics of crystals accounting for nuclear quantum and thermal anharmonic fluctuations. The approach requires the evaluation of the Born-Oppenheimer energy, as well as its derivatives with respect to ionic positions (forces) and cell parameters (stress tensor) in supercells, which can be provided, for instance, by first principles density-functional-theory codes. The method performs crystal geometry relaxation on the quantum free energy landscape, optimizing the free energy with respect to all degrees of freedom of the crystal structure. It can be used to determine the phase diagram of any crystal at finite temperature. It enables the calculation of phase boundaries for both first-order and second-order phase transitions from the Hessian of the free energy. Finally, the code can also compute the anharmonic phonon spectra, including the phonon linewidths, as well as phonon spectral functions. We review the theoretical framework of the stochastic self-consistent harmonic approximation and its dynamical extension, making particular emphasis on the physical interpretation of the variables present in the theory that can enlighten the comparison with any other anharmonic theory. A modular and flexible Python environment is used for the implementation, which allows for a clean interaction with other packages. We briefly present a toy-model calculation to illustrate the potential of the code.
△ Less
Submitted 5 March, 2021;
originally announced March 2021.
-
Time-Dependent Self Consistent Harmonic Approximation: Anharmonic nuclear quantum dynamics and time correlation functions
Authors:
Lorenzo Monacelli,
Francesco Mauri
Abstract:
Most material properties of great physical interest are directly related to nuclear dynamics, e.g. the ionic thermal conductivity, Raman/IR vibrational spectra, inelastic X-ray, and Neutron scattering. A theory able to compute from first principles these properties, accounting for the anharmonicity and quantum fluctuations in the nuclear energy landscape that can be implemented in systems with hun…
▽ More
Most material properties of great physical interest are directly related to nuclear dynamics, e.g. the ionic thermal conductivity, Raman/IR vibrational spectra, inelastic X-ray, and Neutron scattering. A theory able to compute from first principles these properties, accounting for the anharmonicity and quantum fluctuations in the nuclear energy landscape that can be implemented in systems with hundreds of atoms is missing. Here, we derive an approximate theory for the quantum time evolution of lattice vibrations at finite temperature. This theory introduces the time dynamics in the Self-Consistent Harmonic Approximation (SCHA) and shares with the static case the same computational cost. It is nonempirical, as pure states evolve according to the Dirac least action principle and the dynamics of the thermal ensemble conserves both energy and entropy. The static SCHA is recovered as a stationary solution of the dynamical equations. We apply perturbation theory around the static SCHA solution and derive an algorithm to compute efficiently quantum dynamical response functions. Thanks to this new algorithm, we have access to the response function of any general external time-dependent perturbation, enabling the simulation of phonon spectra without following any perturbative expansion of the nuclear potential or empirical methods. We benchmark the algorithm on the IR and Raman spectroscopy of high-pressure hydrogen phase III, with a simulation cell of 96 atoms. Our work also explores the nonlinear regime of the dynamical nuclear motion, providing a paradigm to simulate the interaction with intense or multiple probes, as in pump-probe spectroscopy, or chemical reactions involving light atoms, as the proton transfer in biomolecules
△ Less
Submitted 5 March, 2021; v1 submitted 30 November, 2020;
originally announced November 2020.
-
Phonon collapse and van der Waals melting of the 3D charge density wave of VSe$_2$
Authors:
Josu Diego,
A. H. Said,
S. K. Mahatha,
Raffaello Bianco,
Lorenzo Monacelli,
Matteo Calandra,
Francesco Mauri,
K. Rossnagel,
Ion Errea,
S. Blanco-Canosa
Abstract:
Among transition metal dichalcogenides (TMDs), VSe$_2$ is considered to develop a purely 3-dimensional (3D) charge-density wave (CDW) at T$_{CDW}$=110 K. Here, by means of high resolution inelastic x-ray scattering (IXS), we show that the CDW transition is driven by the collapse of an acoustic mode at the critical wavevector \textit{q}$_{CDW}$= (2.25 0 0.7) r.l.u. and critical temperature T…
▽ More
Among transition metal dichalcogenides (TMDs), VSe$_2$ is considered to develop a purely 3-dimensional (3D) charge-density wave (CDW) at T$_{CDW}$=110 K. Here, by means of high resolution inelastic x-ray scattering (IXS), we show that the CDW transition is driven by the collapse of an acoustic mode at the critical wavevector \textit{q}$_{CDW}$= (2.25 0 0.7) r.l.u. and critical temperature T$_{CDW}$=110 K. The softening of this mode starts to be pronounced for temperatures below 2$\times$ T$_{CDW}$ and expands over a rather wide region of the Brillouin zone, suggesting a large contribution of the electron-phonon interaction to the CDW formation. This interpretation is supported by our first principles calculations that determine a large momentum-dependence of the electron-phonon interaction, peaking at the CDW wavevector, in the presence of nesting. Fully anharmonic {\it ab initio} calculations confirm the softening of one acoustic branch at \textit{q}$_{CDW}$ as responsible for the CDW formation and show that van der Waals interactions are crucial to melt the CDW. Our work also highlights the important role of out-of-plane interactions to describe 3D CDWs in TMDs.
△ Less
Submitted 16 July, 2020;
originally announced July 2020.
-
Theory of the thickness dependence of the charge density wave transition in 1T-TiTe$_2$
Authors:
Jianqiang Sky Zhou,
Raffaello Bianco,
Lorenzo Monacelli,
Ion Errea,
Francesco Mauri,
Matteo Calandra
Abstract:
Most metallic transition metal dichalcogenides undergo charge density wave (CDW) instabilities with similar or identical ordering vectors in bulk and in single layer, albeit with different critical temperatures. Metallic 1T-TiTe$_2$ is a remarkable exception as it shows no evidence of charge density wave formation in bulk, but it displays a stable $2\times2$ reconstruction in single-layer form. Th…
▽ More
Most metallic transition metal dichalcogenides undergo charge density wave (CDW) instabilities with similar or identical ordering vectors in bulk and in single layer, albeit with different critical temperatures. Metallic 1T-TiTe$_2$ is a remarkable exception as it shows no evidence of charge density wave formation in bulk, but it displays a stable $2\times2$ reconstruction in single-layer form. The mechanism for this 3D-2D crossover of the transition is still unclear, although strain from the substrate and the exchange interaction have been pointed out as possible formation mechanisms. Here, by performing non-perturbative anharmonic calculations with gradient corrected and hybrid functionals, we explain the thickness behaviour of the transition in 1T-TiTe$_2$. We demonstrate that the occurrence of the CDW in single-layer TiTe$_2$ occurs from the interplay of non-perturbative anharmonicity and an exchange enhancement of the electron-phonon interaction, larger in the single layer than in the bulk. Finally, we study the electronic and structural properties of the single-layer CDW phase and provide a complete description of its electronic structure, phonon dispersion as well as infrared and Raman active phonon modes.
△ Less
Submitted 13 July, 2020;
originally announced July 2020.
-
Bending rigidity, sound propagation and ripples in flat graphene
Authors:
Unai Aseginolaza,
Josu Diego,
Tommaso Cea,
Raffaello Bianco,
Lorenzo Monacelli,
Francesco Libbi,
Matteo Calandra,
Aitor Bergara,
Francesco Mauri,
Ion Errea
Abstract:
Despite many of the applications of graphene rely on its uneven stiffness and high thermal conductivity, the mechanical properties of graphene, and in general of all 2D materials, are still elusive. The harmonic theory predicts a quadratic dispersion for the flexural acoustic vibrational mode, which leads the unphysical result that long wavelength in-plane acoustic modes decay before vibrating one…
▽ More
Despite many of the applications of graphene rely on its uneven stiffness and high thermal conductivity, the mechanical properties of graphene, and in general of all 2D materials, are still elusive. The harmonic theory predicts a quadratic dispersion for the flexural acoustic vibrational mode, which leads the unphysical result that long wavelength in-plane acoustic modes decay before vibrating one period, preventing the propagation of sound. The robustness of the quadratic dispersion has been questioned by arguing that the anharmonic phonon-phonon interaction linearizes it. However, this implies a divergent bending rigidity in the long wavelength limit not reproduced experimentally. Here we show that rotational symmetry protects the quadratic flexural dispersion against phonon-phonon interactions and that, consequently, the bending stiffness is non-divergent irrespective of the temperature. Our non-perturbative anharmonic calculations also determine that sound propagation coexists with a quadratic dispersion. We also show that the temperature dependence of the height fluctuations of the membrane, known as ripples, is fully determined by thermal or quantum fluctuations, but without the anharmonic suppression of their amplitude previously assumed. The universality of our conclusions reconcile experimental evidence and theory not just in graphene, but all 2D materials.
△ Less
Submitted 7 June, 2023; v1 submitted 25 May, 2020;
originally announced May 2020.
-
Weak Dimensionality Dependence and Dominant Role of Ionic Fluctuations in the Charge-Density-Wave Transition of NbSe$_2$
Authors:
Raffaello Bianco,
Lorenzo Monacelli,
Matteo Calandra,
Francesco Mauri,
Ion Errea
Abstract:
Contradictory experiments have been reported about the dimensionality effect on the charge-density-wave transition in 2H NbSe$_2$. While scanning tunnelling experiments on single layers grown by molecular beam epitaxy measure a charge-density-wave transition temperature in the monolayer similar to the bulk, around 33 K, Raman experiments on exfoliated samples observe a large enhancement of the tra…
▽ More
Contradictory experiments have been reported about the dimensionality effect on the charge-density-wave transition in 2H NbSe$_2$. While scanning tunnelling experiments on single layers grown by molecular beam epitaxy measure a charge-density-wave transition temperature in the monolayer similar to the bulk, around 33 K, Raman experiments on exfoliated samples observe a large enhancement of the transition temperature up to 145 K. By employing a non-perturbative approach to deal with anharmonicity, we calculate from first principles the temperature dependence of the phonon spectra both for bulk and monolayer. In both cases, the charge-density-wave transition temperature is estimated as the temperature at which the phonon energy of the mode driving the structural instability vanishes. The obtained transition temperature in the bulk is around 59 K, in rather good agreement with experiments, and it is just slightly increased in the single-layer limit to 73 K, showing the weak dependence of the transition on dimensionality. Environmental factors could motivate the disagreement between the transition temperatures reported by experiments. Our analysis also demonstrates the predominance of ionic fluctuations over electronic ones in the melting of the charge-density-wave order.
△ Less
Submitted 1 September, 2020; v1 submitted 17 April, 2020;
originally announced April 2020.
-
Black metal hydrogen above 360 GPa driven by proton quantum fluctuations
Authors:
Lorenzo Monacelli,
Ion Errea,
Matteo Calandra,
Francesco Mauri
Abstract:
Production of metallic hydrogen is one of the top three open quests of physics. Recent low-temperature experiments report different metallization pressures, varying from 360GPa to 490GPa.
In this work, we simulate structural properties, vibrational Raman, IR and optical spectra of hydrogen phase III accounting for proton quantum effects. We demonstrate that nuclear quantum fluctuations downshift…
▽ More
Production of metallic hydrogen is one of the top three open quests of physics. Recent low-temperature experiments report different metallization pressures, varying from 360GPa to 490GPa.
In this work, we simulate structural properties, vibrational Raman, IR and optical spectra of hydrogen phase III accounting for proton quantum effects. We demonstrate that nuclear quantum fluctuations downshift the vibron frequencies by 25%, introduce a broad line-shape in the Raman spectra, and reduce the optical gap by 3eV. We show that hydrogen metallization occurs at 380GPa in phase III due to band overlap, in good agreement with transport data.
By simulating the optical properties, we predict this state to be a peculiar black metal, transparent in the IR. The transparent window closes at 450GPa, but the reflectivity remains low, discarding phase III as the shiny metal observed at 490GPa. We predict the conductivity onset to increase by 70GPa and the transparent window to increase by 1.3eV when replacing hydrogen by deuterium at 0K, underlining that metallization is driven by quantum fluctuations and is thus isotope dependent.
We show how hydrogen acquires metallic features (conductivity and brightness) at different pressures, explaining the apparent contradictions in existing experimental scenarios.
△ Less
Submitted 28 January, 2020; v1 submitted 11 December, 2019;
originally announced December 2019.
-
Anharmonic melting of the charge density wave in single-layer TiSe2
Authors:
Jianqiang Sky Zhou,
Lorenzo Monacelli,
Raffaello Bianco,
Ion Errea,
Francesco Mauri,
Matteo Calandra
Abstract:
Low dimensional systems with a vanishing band-gap and a large electron-hole interaction have been proposed to be unstable towards exciton formation. As the exciton binding energy increases in low dimension, conventional wisdom suggests that excitonic insulators should be more stable in 2D than in 3D. Here we study the effects of the electron-hole interaction and anharmonicity in single-layer TiSe2…
▽ More
Low dimensional systems with a vanishing band-gap and a large electron-hole interaction have been proposed to be unstable towards exciton formation. As the exciton binding energy increases in low dimension, conventional wisdom suggests that excitonic insulators should be more stable in 2D than in 3D. Here we study the effects of the electron-hole interaction and anharmonicity in single-layer TiSe2. We find that, contrary to the bulk case and to the generally accepted picture, the electron-hole exchange interaction is much smaller in 2D than in 3D and it has negligible effects on phonon spectra. By calculating anharmonic phonon spectra within the stochastic self-consistent harmonic approximation, we obtain TCDW = 440K for an isolated and undoped single-layer and TCDW = 364K for an electron-doping n = 4.6 x 10^13 cm^{-2} , close to the experimental result of 200-280K on supported samples. Our work demonstrates that anharmonicity and doping melt the charge density wave in single-layer TiSe2.
△ Less
Submitted 28 October, 2019;
originally announced October 2019.
-
Anomalous high-temperature superconductivity in YH$_6$
Authors:
Ivan A. Troyan,
Dmitrii V. Semenok,
Alexander G. Kvashnin,
Andrey V. Sadakov,
Oleg A. Sobolevskiy,
Vladimir M. Pudalov,
Anna G. Ivanova,
Vitali B. Prakapenka,
Eran Greenberg,
Alexander G. Gavriliuk,
Viktor V. Struzhkin,
Aitor Bergara,
Ion Errea,
Raffaello Bianco,
Matteo Calandra,
Francesco Mauri,
Lorenzo Monacelli,
Ryosuke Akashi,
Artem R. Oganov
Abstract:
Pressure-stabilized hydrides are a new rapidly growing class of high-temperature superconductors which is believed to be described within the conventional phonon-mediated mechanism of coupling. Here we report the synthesis of yttrium hexahydride Im3m-YH$_6$ that demonstrates the superconducting transition with T$_c$ = 224 K at 166 GPa, much lower than the theoretically predicted (>270 K). The meas…
▽ More
Pressure-stabilized hydrides are a new rapidly growing class of high-temperature superconductors which is believed to be described within the conventional phonon-mediated mechanism of coupling. Here we report the synthesis of yttrium hexahydride Im3m-YH$_6$ that demonstrates the superconducting transition with T$_c$ = 224 K at 166 GPa, much lower than the theoretically predicted (>270 K). The measured upper critical magnetic field B$_c$$_2$(0) of YH$_6$ was found to be 116-158 T, which is 2-2.5 times larger than the calculated value. A pronounced shift of T$_c$ in yttrium deuteride YD$_6$ with the isotope coefficient 0.4 supports the phonon-assisted superconductivity. Current-voltage measurements showed that the critical current I$_c$ and its density J$_c$ may exceed 1.75 A and 3500 A/mm$^2$ at 0 K, respectively, which is comparable with the parameters of commercial superconductors, such as NbTi and YBCO. The superconducting density functional theory (SCDFT) and anharmonic calculations suggest unusually large impact of the Coulomb repulsion in this compound. The results indicate notable departures of the superconducting properties of the discovered YH$_6$ from the conventional Migdal-Eliashberg and Bardeen-Cooper-Schrieffer theories.
△ Less
Submitted 12 June, 2020; v1 submitted 5 August, 2019;
originally announced August 2019.
-
Quantum Crystal Structure in the 250 K Superconducting Lanthanum Hydride
Authors:
Ion Errea,
Francesco Belli,
Lorenzo Monacelli,
Antonio Sanna,
Takashi Koretsune,
Terumasa Tadano,
Raffaello Bianco,
Matteo Calandra,
Ryotaro Arita,
Francesco Mauri,
José A. Flores-Livas
Abstract:
The discovery of superconductivity at 200 K in the hydrogen sulfide system at large pressures [1] was a clear demonstration that hydrogen-rich materials can be high-temperature superconductors. The recent synthesis of LaH$_{10}$ with a superconducting critical temperature (T$_{\text{c}}$) of 250 K [2,3] places these materials at the verge of reaching the long-dreamed room-temperature superconducti…
▽ More
The discovery of superconductivity at 200 K in the hydrogen sulfide system at large pressures [1] was a clear demonstration that hydrogen-rich materials can be high-temperature superconductors. The recent synthesis of LaH$_{10}$ with a superconducting critical temperature (T$_{\text{c}}$) of 250 K [2,3] places these materials at the verge of reaching the long-dreamed room-temperature superconductivity. Electrical and x-ray diffraction measurements determined a weakly pressure-dependent T$_{\text{c}}$ for LaH$_{10}$ between 137 and 218 gigapascals in a structure with a face-centered cubic (fcc) arrangement of La atoms [3]. Here we show that quantum atomic fluctuations stabilize in all this pressure range a high-symmetry Fm-3m crystal structure consistent with experiments, which has a colossal electron-phonon coupling of $λ\sim3.5$. Even if ab initio classical calculations neglecting quantum atomic vibrations predict this structure to distort below 230 GPa yielding a complex energy landscape with many local minima, the inclusion of quantum effects simplifies the energy landscape evidencing the Fm-3m as the true ground state. The agreement between the calculated and experimental T$_{\text{c}}$ values further supports this phase as responsible for the 250 K superconductivity. The relevance of quantum fluctuations in the energy landscape found here questions many of the crystal structure predictions made for hydrides within a classical approach that at the moment guide the experimental quest for room-temperature superconductivity [4,5,6]. Furthermore, quantum effects reveal crucial to sustain solids with extraordinary electron-phonon coupling that may otherwise be unstable [7].
△ Less
Submitted 27 July, 2019;
originally announced July 2019.
-
Strong anharmonicity and high thermoelectric efficiency in high temperature SnS from first-principles
Authors:
Unai Aseginolaza,
Raffaello Bianco,
Lorenzo Monacelli,
Lorenzo Paulatto,
Matteo Calandra,
Francesco Mauri,
Aitor Bergara,
Ion Errea
Abstract:
SnS and SnSe are isoelectronic materials with a common phase diagram. Recently, SnSe was found to be the most efficient intrinsic thermoelectric material in its high-temperature Cmcm phase above 800 K. Making use of first-principles calculations, here we show that the electronic and vibrational properties of both materials are very similar in this phase and, consequently, SnS is also expected to h…
▽ More
SnS and SnSe are isoelectronic materials with a common phase diagram. Recently, SnSe was found to be the most efficient intrinsic thermoelectric material in its high-temperature Cmcm phase above 800 K. Making use of first-principles calculations, here we show that the electronic and vibrational properties of both materials are very similar in this phase and, consequently, SnS is also expected to have a high thermoelectric figure of merit at high temperature in its Cmcm phase. In fact, the electronic power factor and lattice thermal conductivity are comparable for both materials, which ensures a similar figure of merit. As in the case of SnSe, the vibrational properties of SnS in the Cmcm phase are far from trivial and are dominated by huge anharmonic effects. Its phonon spectra are strongly renormalized by anharmonicity and the spectral functions of some particular in-plane modes depict anomalous non-lorentzian profiles. Finally, we show that non-perturbative anharmonic effects in the third-order force-constants are crucial in the calculation of the lattice thermal conductivity. Our results motivate new experiments in the high temperature regime to measure the figure of merit of SnS.
△ Less
Submitted 5 June, 2019;
originally announced June 2019.
-
Quantum effects in muon spin spectroscopy within the stochastic self-consistent harmonic approximation
Authors:
Ifeanyi John Onuorah,
Pietro Bonfà,
Roberto De Renzi,
Lorenzo Monacelli,
Francesco Mauri,
Matteo Calandra,
Ion Errea
Abstract:
In muon spin rotation experiments the positive implanted muon vibrates with large zero point amplitude by virtue of its light mass. Quantum mechanical calculations of the host material usually treat the muon as a point impurity, ignoring this large vibrational amplitude. As a first order correction, the muon zero point motion is usually described within the harmonic approximation, despite the larg…
▽ More
In muon spin rotation experiments the positive implanted muon vibrates with large zero point amplitude by virtue of its light mass. Quantum mechanical calculations of the host material usually treat the muon as a point impurity, ignoring this large vibrational amplitude. As a first order correction, the muon zero point motion is usually described within the harmonic approximation, despite the large anharmonicity of the crystal potential. Here we apply the stochastic self-consistent harmonic approximation, a quantum variational method devised to include strong anharmonic effects in total energy and vibrational frequency calculations, in order to overcome these limitations and provide an accurate ab initio description of the quantum nature of the muon. We applied this full quantum treatment to the calculation of the muon contact hyperfine field in textbook-case metallic systems, such as Fe, Ni, Co including MnSi and MnGe, significantly improving agreement with experiments. Our results show that muon vibrational frequencies are strongly renormalized by anharmonicity. Finally, in contrast to the harmonic approximation, we show that including quantum anharmonic fluctuations, the muon stabilizes at the octahedral site in bcc Fe.
△ Less
Submitted 19 July, 2019; v1 submitted 26 April, 2019;
originally announced April 2019.
-
Quantum enhancement of charge density wave in NbS$_2$ in the 2D limit
Authors:
Raffaello Bianco,
Ion Errea,
Lorenzo Monacelli,
Matteo Calandra,
Francesco Mauri
Abstract:
At ambient pressure, bulk 2H-NbS$_2$ displays no charge density wave instability at odds with the isostructural and isoelectronic compounds 2H-NbSe$_2$, 2H-TaS$_2$ and 2H-TaSe$_2$, and in disagreement with harmonic calculations. Contradictory experimental results have been reported in supported single layers, as 1H-NbS$_2$ on Au(111) does not display a charge density wave, while 1H-NbS$_2$ on 6H-S…
▽ More
At ambient pressure, bulk 2H-NbS$_2$ displays no charge density wave instability at odds with the isostructural and isoelectronic compounds 2H-NbSe$_2$, 2H-TaS$_2$ and 2H-TaSe$_2$, and in disagreement with harmonic calculations. Contradictory experimental results have been reported in supported single layers, as 1H-NbS$_2$ on Au(111) does not display a charge density wave, while 1H-NbS$_2$ on 6H-SiC(0001) endures a $3\times 3$ reconstruction. Here, by carrying out quantum anharmonic calculations from first-principles, we evaluate the temperature dependence of phonon spectra in NbS$_2$ bulk and single layer as a function of pressure/strain. For bulk 2H-NbS$_2$, we find excellent agreement with inelastic X-ray spectra and demonstrate the removal of charge ordering due to anharmonicity. In the 2D limit, we find an enhanced tendency toward charge density wave order. Freestanding 1H-NbS$_2$ undergoes a $3\times3$ reconstruction, in agreement with data on 6H-SiC(0001) supported samples. Moreover, as strains smaller than $0.5\%$ in the lattice parameter are enough to completely remove the $3\times3$ superstructure, deposition of 1H-NbS$_2$ on flexible substrates or a small charge transfer via field-effect could lead to devices with dynamical switching on/off of charge order.
△ Less
Submitted 24 April, 2019; v1 submitted 2 February, 2019;
originally announced February 2019.
-
Phonon Collapse and Second-Order Phase Transition in Thermoelectric SnSe
Authors:
Unai Aseginolaza,
Raffaello Bianco,
Lorenzo Monacelli,
Lorenzo Paulatto,
Matteo Calandra,
Francesco Mauri,
Aitor Bergara,
Ion Errea
Abstract:
Since 2014 the layered semiconductor SnSe in the high-temperature Cmcm phase is known to be the most efficient thermoelectric material. Making use of first-principles calculations we show that its vibrational and thermal transport properties are determined by huge non-perturbative anharmonic effects. We show that the transition from the Cmcm phase to the low-symmetry Pnma is a second-order phase t…
▽ More
Since 2014 the layered semiconductor SnSe in the high-temperature Cmcm phase is known to be the most efficient thermoelectric material. Making use of first-principles calculations we show that its vibrational and thermal transport properties are determined by huge non-perturbative anharmonic effects. We show that the transition from the Cmcm phase to the low-symmetry Pnma is a second-order phase transition driven by the collapse of a zone border phonon, whose frequency vanishes at the transition temperature. Our calculations show that the spectral function of the in-plane vibrational modes are strongly anomalous with shoulders and double-peak structures. We calculate the lattice thermal conductivity obtaining good agreement with experiments only when non-perturbative anharmonic scattering is included. Our results suggest that the good thermoelectric efficiency of SnSe is strongly affected by the non-perturbative anharmonicity. We expect similar effects to occur in other thermoelectric materials.
△ Less
Submitted 20 July, 2018;
originally announced July 2018.
-
Pressure and stress tensor of complex anharmonic crystals within the stochastic self-consistent harmonic approximation
Authors:
Lorenzo Monacelli,
Ion Errea,
Matteo Calandra,
Francesco Mauri
Abstract:
The self-consisted harmonic approximation (SCHA) allows the computation of free energy of anharmonic crystals considering both quantum and thermal fluctuations. Recently, a stochastic implementation of the SCHA has been developed, tailored for applications that use total energy and forces computed from first principles. In this work, we extend the applicability of the stochastic SCHA to complex cr…
▽ More
The self-consisted harmonic approximation (SCHA) allows the computation of free energy of anharmonic crystals considering both quantum and thermal fluctuations. Recently, a stochastic implementation of the SCHA has been developed, tailored for applications that use total energy and forces computed from first principles. In this work, we extend the applicability of the stochastic SCHA to complex crystals with many degrees of freedom, with the optimisation of both the lattice vectors and the atomic positions. To this goal, we provide an expression for the evaluation of the pressure and stress tensor within the stochastic SCHA formalism. Moreover, we develop a more robust free energy minimisation algorithm, which allows us to perform the SCHA variational minimisation very efficiently in systems having a broad spectrum of phonon frequencies and many degrees of freedom. We test and illustrate the new approach with an application to the phase XI of water ice using density-functional theory. We find that the SCHA reproduces extremely well the experimental thermal expansion of ice in the whole temperature range between 0 K and 270 K, in contrast with the results obtained within the quasi-harmonic approximation, that underestimates the effect by about 25 %.
△ Less
Submitted 18 April, 2018;
originally announced April 2018.
-
Entropy evaluation sheds light on ecosystem complexity
Authors:
Mattia Miotto,
Lorenzo Monacelli
Abstract:
Preserving biodiversity and ecosystem stability is a challenge that can be pursued through modern statistical mechanics modeling. Here we introduce a variational maximum entropy-based algorithm to evaluate the entropy in a minimal ecosystem on a lattice in which two species struggle for survival. The method quantitatively reproduces the scale-free law of the prey shoals size, where the simpler mea…
▽ More
Preserving biodiversity and ecosystem stability is a challenge that can be pursued through modern statistical mechanics modeling. Here we introduce a variational maximum entropy-based algorithm to evaluate the entropy in a minimal ecosystem on a lattice in which two species struggle for survival. The method quantitatively reproduces the scale-free law of the prey shoals size, where the simpler mean-field approach fails: the direct near neighbor correlations are found to be the fundamental ingredient describing the system self-organized behavior. Furthermore, entropy allows the measurement of structural ordering, that is found to be a key ingredient in characterizing two different coexistence behaviors, one where predators form localized patches in a sea of preys and another where species display more complex patterns. The general nature of the introduced method paves the way for its application in many other systems of interest.
△ Less
Submitted 10 October, 2018; v1 submitted 30 March, 2018;
originally announced March 2018.