-
Adaptive energy reference for machine-learning models of the electronic density of states
Authors:
Wei Bin How,
Sanggyu Chong,
Federico Grasselli,
Kevin K. Huguenin-Dumittan,
Michele Ceriotti
Abstract:
The electronic density of states (DOS) provides information regarding the distribution of electronic states in a material, and can be used to approximate its optical and electronic properties and therefore guide computational material design. Given its usefulness and relative simplicity, it has been one of the first electronic properties used as target for machine-learning approaches going beyond…
▽ More
The electronic density of states (DOS) provides information regarding the distribution of electronic states in a material, and can be used to approximate its optical and electronic properties and therefore guide computational material design. Given its usefulness and relative simplicity, it has been one of the first electronic properties used as target for machine-learning approaches going beyond interatomic potentials. A subtle but important point, well-appreciated in the condensed matter community but usually overlooked in the construction of data-driven models, is that for bulk configurations the absolute energy reference of single-particle energy levels is ill-defined. Only energy differences matter, and quantities derived from the DOS are typically independent on the absolute alignment. We introduce an adaptive scheme that optimizes the energy reference of each structure as part of training, and show that it consistently improves the quality of ML models compared to traditional choices of energy reference, for different classes of materials and different model architectures. On a practical level, we trace the improved performance to the ability of this self-aligning scheme to match the most prominent features in the DOS. More broadly, we believe that this work highlights the importance of incorporating insights into the nature of the physical target into the definition of the architecture and of the appropriate figures of merit for machine-learning models, that translate in better transferability and overall performance.
△ Less
Submitted 2 July, 2024; v1 submitted 1 July, 2024;
originally announced July 2024.
-
Thermal transport of glasses via machine learning driven simulations
Authors:
Paolo Pegolo,
Federico Grasselli
Abstract:
Accessing the thermal transport properties of glasses is a major issue for the design of production strategies of glass industry, as well as for the plethora of applications and devices where glasses are employed. From the computational standpoint, the chemical and morphological complexity of glasses calls for atomistic simulations where the interatomic potentials are able to capture the variety o…
▽ More
Accessing the thermal transport properties of glasses is a major issue for the design of production strategies of glass industry, as well as for the plethora of applications and devices where glasses are employed. From the computational standpoint, the chemical and morphological complexity of glasses calls for atomistic simulations where the interatomic potentials are able to capture the variety of local environments, composition, and (dis)order that typically characterize glassy phases. Machine-learning potentials (MLPs) are emerging as a valid alternative to computationally expensive ab initio simulations, inevitably run on very small samples which cannot account for disorder at different scales, as well as to empirical force fields, fast but often reliable only in a narrow portion of the thermodynamic and composition phase diagrams. In this article, we make the point on the use of MLPs to compute the thermal conductivity of glasses, through a review of recent theoretical and computational tools and a series of numerical applications on vitreous silica and vitreous silicon, both pure and intercalated with lithium.
△ Less
Submitted 9 February, 2024;
originally announced February 2024.
-
Thermal conductivity of Li$_3$PS$_4$ solid electrolytes with ab initio accuracy
Authors:
Davide Tisi,
Federico Grasselli,
Lorenzo Gigli,
Michele Ceriotti
Abstract:
The vast amount of computational studies on electrical conduction in solid-state electrolytes is not mirrored by comparable efforts addressing thermal conduction, which has been scarcely investigated despite its relevance to thermal management and (over)heating of batteries. The reason for this lies in the complexity of the calculations: on one hand, the diffusion of ionic charge carriers makes la…
▽ More
The vast amount of computational studies on electrical conduction in solid-state electrolytes is not mirrored by comparable efforts addressing thermal conduction, which has been scarcely investigated despite its relevance to thermal management and (over)heating of batteries. The reason for this lies in the complexity of the calculations: on one hand, the diffusion of ionic charge carriers makes lattice methods formally unsuitable, due to the lack of equilibrium atomic positions needed for normal-mode expansion. On the other hand, the prohibitive cost of large-scale molecular dynamics (MD) simulations of heat transport in large systems at ab initio levels has hindered the use of MD-based methods. In this paper, we leverage recently developed machine-learning potentials targeting different ab initio functionals (PBEsol, r$^2$SCAN, PBE0) and a state-of-the-art formulation of the Green-Kubo theory of heat transport in multicomponent systems to compute the thermal conductivity of a promising solid-state electrolyte, Li$_3$PS$_4$, in all its polymorphs ($α$, $β$, and $γ$). By comparing MD estimates with lattice methods on the low-temperature, nondiffusive $γ$-Li$_3$PS$_4$, we highlight strong anharmonicities and negligible nuclear quantum effects, hence further justifying MD-based methods even for nondiffusive phases. Finally, for the ion-conducting $α$ and $β$ phases, where the multicomponent Green-Kubo MD approach is mandatory, our simulations indicate a weak temperature dependence of the thermal conductivity, a glass-like behavior due to the effective local disorder characterizing these Li-diffusing phases.
△ Less
Submitted 17 June, 2024; v1 submitted 23 January, 2024;
originally announced January 2024.
-
Mechanism of charge transport in lithium thiophosphate
Authors:
Lorenzo Gigli,
Davide Tisi,
Federico Grasselli,
Michele Ceriotti
Abstract:
Lithium ortho-thiophosphate (Li$_3$PS$_4$) has emerged as a promising candidate for solid-state-electrolyte batteries, thanks to its highly conductive phases, cheap components, and large electrochemical stability range. Nonetheless, the microscopic mechanisms of Li-ion transport in Li$_3$PS$_4$ are far to be fully understood, the role of PS$_4$ dynamics in charge transport still being controversia…
▽ More
Lithium ortho-thiophosphate (Li$_3$PS$_4$) has emerged as a promising candidate for solid-state-electrolyte batteries, thanks to its highly conductive phases, cheap components, and large electrochemical stability range. Nonetheless, the microscopic mechanisms of Li-ion transport in Li$_3$PS$_4$ are far to be fully understood, the role of PS$_4$ dynamics in charge transport still being controversial. In this work, we build machine learning potentials targeting state-of-the-art DFT references (PBEsol, r$^2$SCAN, and PBE0) to tackle this problem in all known phases of Li$_3$PS$_4$ ($α$, $β$ and $γ$), for large system sizes and timescales. We discuss the physical origin of the observed superionic behavior of Li$_3$PS$_4$: the activation of PS$_4$ flipping drives a structural transition to a highly conductive phase, characterized by an increase of Li-site availability and by a drastic reduction in the activation energy of Li-ion diffusion. We also rule out any paddle-wheel effects of PS$_4$ tetrahedra in the superionic phases -- previously claimed to enhance Li-ion diffusion -- due to the orders-of-magnitude difference between the rate of PS$_4$ flips and Li-ion hops at all temperatures below melting. We finally elucidate the role of inter-ionic dynamical correlations in charge transport, by highlighting the failure of the Nernst-Einstein approximation to estimate the electrical conductivity. Our results show a strong dependence on the target DFT reference, with PBE0 yielding the best quantitative agreement with experimental measurements not only for the electronic band-gap but also for the electrical conductivity of $β$- and $α$-Li$_3$PS$_4$.
△ Less
Submitted 10 January, 2024; v1 submitted 24 October, 2023;
originally announced October 2023.
-
Robustness of Local Predictions in Atomistic Machine Learning Models
Authors:
Sanggyu Chong,
Federico Grasselli,
Chiheb Ben Mahmoud,
Joe D. Morrow,
Volker L. Deringer,
Michele Ceriotti
Abstract:
Machine learning (ML) models for molecules and materials commonly rely on a decomposition of the global target quantity into local, atom-centered contributions. This approach is convenient from a computational perspective, enabling large-scale ML-driven simulations with a linear-scaling cost, and also allow for the identification and post-hoc interpretation of contributions from individual chemica…
▽ More
Machine learning (ML) models for molecules and materials commonly rely on a decomposition of the global target quantity into local, atom-centered contributions. This approach is convenient from a computational perspective, enabling large-scale ML-driven simulations with a linear-scaling cost, and also allow for the identification and post-hoc interpretation of contributions from individual chemical environments and motifs to complicated macroscopic properties. However, even though there exist practical justifications for these decompositions, only the global quantity is rigorously defined, and thus it is unclear to what extent the atomistic terms predicted by the model can be trusted. Here, we introduce a quantitative metric, which we call the local prediction rigidity (LPR), that allows one to assess how robust the locally decomposed predictions of ML models are. We investigate the dependence of LPR on the aspects of model training, particularly the composition of training dataset, for a range of different problems from simple toy models to real chemical systems. We present strategies to systematically enhance the LPR, which can be used to improve the robustness, interpretability, and transferability of atomistic ML models.
△ Less
Submitted 27 June, 2023;
originally announced June 2023.
-
Self-interaction and transport of solvated electrons in molten salts
Authors:
Paolo Pegolo,
Stefano Baroni,
Federico Grasselli
Abstract:
The dynamics of (few) electrons dissolved in an ionic fluid--as when a small amount of metal is added to a solution while upholding its electronic insulation--manifests interesting properties that can be ascribed to nontrivial topological features of particle transport (e.g., Thouless' pumps). In the adiabatic regime, the charge distribution and the dynamics of these dissolved electrons are unique…
▽ More
The dynamics of (few) electrons dissolved in an ionic fluid--as when a small amount of metal is added to a solution while upholding its electronic insulation--manifests interesting properties that can be ascribed to nontrivial topological features of particle transport (e.g., Thouless' pumps). In the adiabatic regime, the charge distribution and the dynamics of these dissolved electrons are uniquely determined by the nuclear configuration. Yet, their localization into effective potential wells and their diffusivity are dictated by how the self-interaction is modeled. In this article, we investigate the role of self-interaction in the description of localization and transport properties of dissolved electrons in non-stoichiometric molten salts. Although the account for the exact (Fock) exchange strongly localizes the dissolved electrons, decreasing their tunneling probability and diffusivity, we show that the dynamics of the ions and of the dissolved electrons are largely uncorrelated, irrespective of the degree to which the electron self-interaction is treated, and in accordance with topological arguments.
△ Less
Submitted 7 September, 2023; v1 submitted 17 May, 2023;
originally announced May 2023.
-
Effect of a Temperature Gradient on the Screening Properties of Ionic Fluids
Authors:
Andrea Grisafi,
Federico Grasselli
Abstract:
The electrostatic screening properties of ionic fluids are of paramount importance in countless physical processes. Yet, the behavior of ionic conductors out of thermal equilibrium has to date mainly been studied in the context of thermodiffusion phenomena by virtue of direct extensions of Debye-Hückel theories. We investigate how the static response of a symmetric ionic fluid is influenced by the…
▽ More
The electrostatic screening properties of ionic fluids are of paramount importance in countless physical processes. Yet, the behavior of ionic conductors out of thermal equilibrium has to date mainly been studied in the context of thermodiffusion phenomena by virtue of direct extensions of Debye-Hückel theories. We investigate how the static response of a symmetric ionic fluid is influenced by the presence of a thermal gradient by introducing a theory of electrostatic screening under a stationary temperature profile. By borrowing mathematical methods commonly used in the semiclassical approximation of quantum particles, we find analytical solutions to the asymptotic decay of the charge density which can be used to describe the non-equilibrium response of the system to external charge perturbations and for arbitrary ionic concentrations. Notably, a transition between monotonic and oscillatory screening regimes is observed as an effect of the temperature variation which generalizes known results of thermal equilibrium to out of equilibrium conditions. A final quantitative example on the screening of charged surfaces in aqueous electrolytes shows that the deviation from thermal equilibrium predicted by our solutions is generally larger than thermodiffusion effects, and should therefore be taken into account for a comprehensive description of the electrical double layer. Our findings pave the way to the rigorous treatment of non-equilibrium steady states in ionic systems with potential applications to the study of energy materials, nanostructured systems and waste-heat-recovery technologies.
△ Less
Submitted 18 April, 2023; v1 submitted 23 December, 2022;
originally announced December 2022.
-
Topology, oxidation states, and charge transport in ionic conductors
Authors:
Paolo Pegolo,
Stefano Baroni,
Federico Grasselli
Abstract:
Recent theoretical advances, based on a combination of concepts from Thouless' topological theory of adiabatic charge transport and a newly introduced gauge-invariance principle for transport coefficients, have permitted to connect (and reconcile) Faraday's picture of ionic transport -- whereby each atom carries a well-defined integer charge -- with a rigorous quantum description of the electronic…
▽ More
Recent theoretical advances, based on a combination of concepts from Thouless' topological theory of adiabatic charge transport and a newly introduced gauge-invariance principle for transport coefficients, have permitted to connect (and reconcile) Faraday's picture of ionic transport -- whereby each atom carries a well-defined integer charge -- with a rigorous quantum description of the electronic charge-density distribution, which hardly suggests its partition into well defined atomic contributions. In this paper we review these progress and in particular that, by relaxing some general topological conditions, charge may be transported in ionic conductors without any net ionic displacements. After reporting numerical experiments which corroborate these findings, we introduce a new connection between our topological picture and the well-known Marcus-Hush theory of electron transfer, which we are able to connect with the topology of adiabatic paths drawn by atomic trajectories. As a significant byproduct, the results reviewed here permit to classify different regimes of ionic transport according to the topological properties of the electronic structure of the conducting material. We finally report on a few recent applications to energy materials and planetary sciences.
△ Less
Submitted 15 June, 2022;
originally announced June 2022.
-
Predicting hot-electron free energies from ground-state data
Authors:
Chiheb Ben Mahmoud,
Federico Grasselli,
Michele Ceriotti
Abstract:
Machine-learning potentials are usually trained on the ground-state, Born-Oppenheimer energy surface, which depends exclusively on the atomic positions and not on the simulation temperature. This disregards the effect of thermally-excited electrons, that is important in metals, and essential to the description of warm dense matter. An accurate physical description of these effects requires that th…
▽ More
Machine-learning potentials are usually trained on the ground-state, Born-Oppenheimer energy surface, which depends exclusively on the atomic positions and not on the simulation temperature. This disregards the effect of thermally-excited electrons, that is important in metals, and essential to the description of warm dense matter. An accurate physical description of these effects requires that the nuclei move on a temperature-dependent electronic free energy. We propose a method to obtain machine-learning predictions of this free energy at an arbitrary electron temperature using exclusively training data from ground-state calculations, avoiding the need to train temperature-dependent potentials, and benchmark it on metallic liquid hydrogen at the conditions of the core of gas giants and brown dwarfs. This work demonstrates the advantages of hybrid schemes that use physical consideration to combine machine-learning predictions, providing a blueprint for the development of similar approaches that extend the reach of atomistic modelling by removing the barrier between physics and data-driven methodologies.
△ Less
Submitted 28 September, 2022; v1 submitted 11 May, 2022;
originally announced May 2022.
-
Investigating finite-size effects in molecular dynamics simulations of ion diffusion, heat transport, and thermal motion in superionic materials
Authors:
Federico Grasselli
Abstract:
The effects of the finite size of the simulation box in equilibrium molecular dynamics simulations are investigated for prototypical superionic conductors of different types, namely the fluorite-structure materials PbF2, CaF2, and UO2 (type-II), and the α phase of AgI (type I). Largely validated empirical force-fields are employed to run ns-long simulations and extract general trends for several p…
▽ More
The effects of the finite size of the simulation box in equilibrium molecular dynamics simulations are investigated for prototypical superionic conductors of different types, namely the fluorite-structure materials PbF2, CaF2, and UO2 (type-II), and the α phase of AgI (type I). Largely validated empirical force-fields are employed to run ns-long simulations and extract general trends for several properties, at increasing size and in a wide temperature range. This work shows that, for the considered type-II superionic conductors, the diffusivity dramatically depends on the system size and that the superionic regime is shifted to larger temperatures in smaller cells. Furthermore, only simulations of several hundred atoms are able to capture the experimentally-observed, characteristic change in the activation energy of the diffusion process, occurring at the order-disorder transition to the superionic regime. Finite-size effects on ion diffusion are instead much weaker in α-AgI. The thermal conductivity is found generally smaller for smaller cells, where the temperature-independent (Allen-Feldman) regime is also reached at significantly lower temperatures. The finite-size effects on the thermal motion of the non-mobile ions composing the solid matrix follow the simple law which holds for solids.
△ Less
Submitted 17 March, 2022;
originally announced March 2022.
-
Temperature- and vacancy-concentration-dependence of heat transport in Li$_3$ClO from multi-method numerical simulations
Authors:
Paolo Pegolo,
Stefano Baroni,
Federico Grasselli
Abstract:
Despite governing heat management in any realistic device, the microscopic mechanisms of heat transport in all-solid-state electrolytes are poorly known: existing calculations, all based on simplistic semi-empirical models, are unreliable for superionic conductors and largely overestimate their thermal conductivity. In this work, we deploy a combination of state-of-the-art methods to calculate the…
▽ More
Despite governing heat management in any realistic device, the microscopic mechanisms of heat transport in all-solid-state electrolytes are poorly known: existing calculations, all based on simplistic semi-empirical models, are unreliable for superionic conductors and largely overestimate their thermal conductivity. In this work, we deploy a combination of state-of-the-art methods to calculate the thermal conductivity of a prototypical Li-ion conductor, the Li$_3$ClO antiperovskite. By leveraging ab initio, machine-learning, and force-field descriptions of inter-atomic forces, we are able to reveal the massive role of anharmonic interactions and diffusive defects on the thermal conductivity and its temperature dependence, and to eventually embed their effects into a simple rationale which is likely applicable to a wide class of ionic conductors.
△ Less
Submitted 23 December, 2021; v1 submitted 29 November, 2021;
originally announced November 2021.
-
Invariance principles in the theory and computation of transport coefficients
Authors:
Federico Grasselli,
Stefano Baroni
Abstract:
In this work we elaborate on two recently discovered invariance principles, according to which transport coefficients are, to a large extent, independent of the microscopic definition of the densities and currents of the conserved quantities being transported (energy, momentum, mass, charge). The first such principle, gauge invariance, allows one to define a quantum adiabatic energy current from d…
▽ More
In this work we elaborate on two recently discovered invariance principles, according to which transport coefficients are, to a large extent, independent of the microscopic definition of the densities and currents of the conserved quantities being transported (energy, momentum, mass, charge). The first such principle, gauge invariance, allows one to define a quantum adiabatic energy current from density-functional theory, from which the heat conductivity can be uniquely defined and computed using equilibrium ab initio molecular dynamics. When combined with a novel topological definition of atomic oxidation states, gauge invariance also sheds new light onto the mechanisms of charge transport in ionic conductors. The second principle, convective invariance, allows one to extend the analysis to multi-component systems. These invariance principles can be combined with new spectral analysis methods for the current time series to be fed into the Green-Kubo formula to obtain accurate estimates of transport coefficients from relatively short molecular dynamics simulations.
△ Less
Submitted 18 July, 2021; v1 submitted 5 May, 2021;
originally announced May 2021.
-
Oxidation states, Thouless' pumps, and non-trivial ionic transport in non-stoichiometric electrolytes
Authors:
Paolo Pegolo,
Federico Grasselli,
Stefano Baroni
Abstract:
Thouless' quantization of adiabatic particle transport permits to associate an integer topological charge with each atom of an electronically gapped material. If these charges are additive and independent of atomic positions, they provide a rigorous definition of atomic oxidation states and atoms can be identified as integer-charge carriers in ionic conductors. Whenever these conditions are met, c…
▽ More
Thouless' quantization of adiabatic particle transport permits to associate an integer topological charge with each atom of an electronically gapped material. If these charges are additive and independent of atomic positions, they provide a rigorous definition of atomic oxidation states and atoms can be identified as integer-charge carriers in ionic conductors. Whenever these conditions are met, charge transport is necessarily convective, i.e. it cannot occur without substantial ionic flow, a transport regime that we dub trivial. We show that the topological requirements that allow these conditions to be broken are the same that would determine a Thouless' pump mechanism if the system were subject to a suitably defined time-periodic Hamiltonian. The occurrence of these requirements determines a non-trivial transport regime whereby charge can flow without any ionic convection, even in electronic insulators. These results are first demonstrated with a couple of simple molecular models that display a quantum pump mechanism upon introduction of a fictitious time dependence of the atomic positions along a closed loop in configuration space. We finally examine the impact of our findings on the transport properties of non-stoichiometric alkali-halide melts, where the same topological conditions that would induce a quantum pump mechanism along certain closed loops in configuration space also determine a non-trivial transport regime such that most of the total charge current results to be uncorrelated from the ionic ones.
△ Less
Submitted 30 September, 2020; v1 submitted 30 June, 2020;
originally announced June 2020.
-
Heat and charge transport in H$_2$O at ice-giant conditions from ab initio molecular dynamics simulations
Authors:
Federico Grasselli,
Lars Stixrude,
Stefano Baroni
Abstract:
The impact of the inner structure and thermal history of planets on their observable features, such as luminosity or magnetic field, crucially depends on the poorly known heat and charge transport properties of their internal layers. The thermal and electric conductivities of different phases of water (liquid, solid, and super-ionic) occurring in the interior of ice giant planets, such as Uranus o…
▽ More
The impact of the inner structure and thermal history of planets on their observable features, such as luminosity or magnetic field, crucially depends on the poorly known heat and charge transport properties of their internal layers. The thermal and electric conductivities of different phases of water (liquid, solid, and super-ionic) occurring in the interior of ice giant planets, such as Uranus or Neptune, are evaluated from equilibrium ab initio molecular dynamics, leveraging recent progresses in the theory and data analysis of transport in extended systems. The implications of our findings on the evolution models of the ice giants are briefly discussed
△ Less
Submitted 14 July, 2020; v1 submitted 27 March, 2020;
originally announced March 2020.
-
Classical and quantum dynamics of indirect excitons driven by surface acoustic waves
Authors:
Federico Grasselli,
Andrea Bertoni,
Guido Goldoni
Abstract:
We perform explicit time-dependent classical and quantum propagation of a spatially indirect exciton (SIX) driven by surface acoustic waves (SAWs) in a semiconductor heterostructure device. We model the SIX dynamics at different levels of description, from the Euler-Lagrange propagation of structureless classical particles to unitary Schrödinger propagation of an electron-hole wave packet in a mea…
▽ More
We perform explicit time-dependent classical and quantum propagation of a spatially indirect exciton (SIX) driven by surface acoustic waves (SAWs) in a semiconductor heterostructure device. We model the SIX dynamics at different levels of description, from the Euler-Lagrange propagation of structureless classical particles to unitary Schrödinger propagation of an electron-hole wave packet in a mean field and to the full quantum propagation of the two-particle complex. A recently proposed beyond mean-field self-energy approach, adding internal virtual transitions to the c.m. dynamics, has been generalized to time-dependent potentials and turns out to describe very well full quantum calculations, while being orders of magnitude numerically less demanding. We show that SAW-driven SIXs are a sensitive probe of scattering potentials in the devices originating, for example, from single impurities or metallic gates, due to competing length and energy scales between the SAW elastic potential, the scattering potential, and the internal electron-hole dynamic of the SIX. Comparison between different approximations allow us to show that internal correlation of the electron-hole pair is crucial in scattering from shallow impurities, where tunneling plays a major role. On the other hand, scattering from broad potentials, i.e., with length scales exceeding the SIX Bohr radius, is well described as the classical dynamics of a pointlike SIX. Recent experiments are discussed in light of our calculations.
△ Less
Submitted 20 February, 2019;
originally announced February 2019.
-
Topological quantisation and gauge invariance of charge transport in liquid insulators
Authors:
Federico Grasselli,
Stefano Baroni
Abstract:
According to the Green-Kubo theory of linear response, the conductivity of an electronically gapped liquid can be expressed in terms of the time correlations of the adiabatic charge flux, which is determined by the atomic velocities and Born effective charges. We show that topological quantisation of adiabatic charge transport and gauge invariance of transport coefficients allow one to rigorously…
▽ More
According to the Green-Kubo theory of linear response, the conductivity of an electronically gapped liquid can be expressed in terms of the time correlations of the adiabatic charge flux, which is determined by the atomic velocities and Born effective charges. We show that topological quantisation of adiabatic charge transport and gauge invariance of transport coefficients allow one to rigorously express the electrical conductivity of an insulating fluid in terms of integer-valued, scalar, and time-independent atomic oxidation numbers, instead of real-valued, tensor, and time-dependent Born charges.
△ Less
Submitted 8 May, 2019; v1 submitted 19 February, 2019;
originally announced February 2019.
-
Theory and numerical simulation of heat transport in multi-component systems
Authors:
Riccardo Bertossa,
Federico Grasselli,
Loris Ercole,
Stefano Baroni
Abstract:
The thermal conductivity of classical multi-component fluids is seemingly affected by the intrinsic arbitrariness in the definition of the atomic energies and it is ill-conditioned numerically, when evaluated from the Green-Kubo theory of linear response. To cope with these two problems we introduce two new concepts: a convective invariance principle for transport coefficients, in the first case,…
▽ More
The thermal conductivity of classical multi-component fluids is seemingly affected by the intrinsic arbitrariness in the definition of the atomic energies and it is ill-conditioned numerically, when evaluated from the Green-Kubo theory of linear response. To cope with these two problems we introduce two new concepts: a convective invariance principle for transport coefficients, in the first case, and multi-variate cepstral analysis, in the second. A combination of these two concepts allows one to substantially reduce the noise affecting the estimate of the thermal conductivity from equilibrium molecular dynamics, even for one-component systems.
△ Less
Submitted 15 May, 2019; v1 submitted 9 August, 2018;
originally announced August 2018.
-
Heat transport in insulators from ab initio Green-Kubo theory
Authors:
Stefano Baroni,
Riccardo Bertossa,
Loris Ercole,
Federico Grasselli,
Aris Marcolongo
Abstract:
The Green-Kubo theory of thermal transport has long be considered incompatible with modern simulation methods based on electronic-structure theory, because it is based on such concepts as energy density and current, which are ill-defined at the quantum-mechanical level. Besides, experience with classical simulations indicates that the estimate of heat-transport coefficients requires analysing mole…
▽ More
The Green-Kubo theory of thermal transport has long be considered incompatible with modern simulation methods based on electronic-structure theory, because it is based on such concepts as energy density and current, which are ill-defined at the quantum-mechanical level. Besides, experience with classical simulations indicates that the estimate of heat-transport coefficients requires analysing molecular trajectories that are more than one order of magnitude longer than deemed feasible using ab initio molecular dynamics. In this paper we report on recent theoretical advances that are allowing one to overcome these two obstacles. First, a general gauge invariance principle has been established, stating that thermal conductivity is insensitive to many details of the microscopic expression for the energy density and current from which it is derived, thus permitting to establish a rigorous expression for the energy flux from Density-Functional Theory, from which the conductivity can be computed in practice. Second, a novel data analysis method based on the statistical theory of time series has been proposed, which allows one to considerably reduce the simulation time required to achieve a target accuracy on the computed conductivity. These concepts are illustrated in detail, starting from a pedagogical introduction to the Green-Kubo theory of linear response and transport, and demonstrated with a few applications done with both classical and quantum-mechanical simulation methods.
△ Less
Submitted 16 August, 2018; v1 submitted 22 February, 2018;
originally announced February 2018.
-
Exact two-body quantum dynamics of an electron-hole pair in semiconductor coupled quantum wells: a time-dependent approach
Authors:
Federico Grasselli,
Andrea Bertoni,
Guido Goldoni
Abstract:
We simulate the time-dependent coherent dynamics of a spatially indirect exciton (an electron-hole pair with the two particles confined in different layers) in a GaAs coupled quantum well system. We use a unitary wave-packet propagation method taking into account in full the four degrees of freedom of the two particles in a two-dimensional system, including both the long-range Coulomb attraction a…
▽ More
We simulate the time-dependent coherent dynamics of a spatially indirect exciton (an electron-hole pair with the two particles confined in different layers) in a GaAs coupled quantum well system. We use a unitary wave-packet propagation method taking into account in full the four degrees of freedom of the two particles in a two-dimensional system, including both the long-range Coulomb attraction and arbitrary two-dimensional electrostatic potentials affecting the electron and/or the hole separately. The method has been implemented for massively parallel architectures to cope with the huge numerical problem, showing good scaling properties and allowing evolution for tens of picoseconds. We have investigated both transient time phenomena and asymptotic time transmission and reflection coefficients for potential profiles consisting of i) extended barriers and wells and ii) a single-slit geometry. We found clear signatures of the internal two-body dynamics, with transient phenomena in the picosecond time-scale which might be revealed by optical spectroscopy. Exact results have been compared with mean-field approaches which, neglecting dynamical correlations by construction, turn out to be inadequate to describe the electron-hole pair evolution in realistic experimental conditions.
△ Less
Submitted 28 April, 2016;
originally announced April 2016.
-
Space- and time-dependent quantum dynamics of spatially indirect excitons in semiconductor heterostructures
Authors:
Federico Grasselli,
Andrea Bertoni,
Guido Goldoni
Abstract:
We study the unitary propagation of a two-particle one-dimensional Schrödinger equation by means of the Split-Step Fourier method, to study the coherent evolution of a spatially indirect exciton (IX) in semiconductor heterostructures. The mutual Coulomb interaction of the electron-hole pair and the electrostatic potentials generated by external gates and acting on the two particles separately are…
▽ More
We study the unitary propagation of a two-particle one-dimensional Schrödinger equation by means of the Split-Step Fourier method, to study the coherent evolution of a spatially indirect exciton (IX) in semiconductor heterostructures. The mutual Coulomb interaction of the electron-hole pair and the electrostatic potentials generated by external gates and acting on the two particles separately are taken into account exactly in the two-particle dynamics. As relevant examples, step/downhill and barrier/well potential profiles are considered. The space- and time-dependent evolution during the scattering event as well as the asymptotic time behavior are analyzed. For typical parameters of GaAs-based devices the transmission or reflection of the pair turns out to be a complex two-particle process, due to comparable and competing Coulomb, electrostatic and kinetic energy scales. Depending on the intensity and anisotropy of the scattering potentials, the quantum evolution may result in excitation of the IX internal degrees of freedom, dissociation of the pair, or transmission in small periodic IX wavepackets due to dwelling of one particle in the barrier region. We discuss the occurrence of each process in the full parameter space of the scattering potentials and the relevance of our results for current excitronic technologies.
△ Less
Submitted 15 April, 2015;
originally announced April 2015.