Michele Simoncelli
@cam.ac.uk
Patricia Crone Research Fellow, Gonville & Caius College, Cavendish Laboratory, Theory of Condensed Matter Group
University of Cambridge
RESEARCH, TEACHING, or OTHER INTERESTS
Condensed Matter Physics, Electronic, Optical and Magnetic Materials, Electrical and Electronic Engineering, Fluid Flow and Transfer Processes
Scopus Publications
Scholar Citations
Scholar h-index
Scholar i10-index
Scopus Publications
Angela F. Harper, Kamil Iwanowski, William C. Witt, Mike C. Payne, and Michele Simoncelli
American Physical Society (APS)
Amorphous alumina is employed ubiquitously as a high-dielectric-constant material in electronics, and its thermal-transport properties are of key relevance for heat management in electronic chips and devices. Experiments show that the thermal conductivity of alumina depends significantly on the synthesis process, indicating the need for a theoretical study to elucidate the atomistic origin of these variations. Here we employ first-principles simulations to characterize the atomistic structure, vibrational properties, and thermal conductivity of alumina at densities ranging from 2.28 g/cm3 to 3.49 g/cm3. Moreover, using an interatomic potential trained on first-principles data, we investigate how system size affects predictions of the thermal conductivity, showing that simulations containing 120 atoms can already reproduce the bulk limit of the conductivity. Finally, relying on the recently developed Wigner formulation of thermal transport, we shed light on the interplay between atomistic topological disorder and anharmonicity in the context of heat conduction, showing that the former dominates over the latter in determining the conductivity of alumina.
Michele Simoncelli, Francesco Mauri, and Nicola Marzari
Springer Science and Business Media LLC
AbstractPredicting the thermal conductivity of glasses from first principles has hitherto been a very complex problem. The established Allen-Feldman and Green-Kubo approaches employ approximations with limited validity—the former neglects anharmonicity, the latter misses the quantum Bose-Einstein statistics of vibrations—and require atomistic models that are very challenging for first-principles methods. Here, we present a protocol to determine from first principles the thermal conductivity κ(T) of glasses above the plateau (i.e., above the temperature-independent region appearing almost without exceptions in the κ(T) of all glasses at cryogenic temperatures). The protocol combines the Wigner formulation of thermal transport with convergence-acceleration techniques, and accounts comprehensively for the effects of structural disorder, anharmonicity, and Bose-Einstein statistics. We validate this approach in vitreous silica, showing that models containing less than 200 atoms can already reproduce κ(T) in the macroscopic limit. We discuss the effects of anharmonicity and the mechanisms determining the trend of κ(T) at high temperature, reproducing experiments at temperatures where radiative effects remain negligible.
Enrico Di Lucente, Michele Simoncelli, and Nicola Marzari
American Physical Society (APS)
Skutterudites are crystals with a cage-like structure that can be augmented with filler atoms ("rattlers"), usually leading to a reduction in thermal conductivity that can be exploited for thermoelectric applications. Here, we leverage the recently introduced Wigner formulation of thermal transport to elucidate the microscopic physics underlying heat conduction in skutterudites, showing that filler atoms can drive a crossover from the Boltzmann to the Wigner regimes of thermal transport, i.e., from particle-like conduction to wave-like tunnelling. At temperatures where the thermoelectric efficiency of skutterudites is largest, wave-like tunneling can become comparable to particle-like propagation. We define a Boltzmann deviation descriptor able to differentiate the two regimes and relate the competition between the two mechanisms to the materials' chemistry, providing a design strategy to select rattlers and identify optimal compositions for thermoelectric applications.
Michele Simoncelli, Nicola Marzari, and Francesco Mauri
American Physical Society (APS)
Two different heat-transport mechanisms are discussed in solids: in crystals, heat carriers propagate and scatter particle-like as described by Peierls’ formulation of the Boltzmann transport equation for phonon wavepackets. In glasses, instead, carriers behave wave-like, diffusing via a Zener-like tunneling between quasi-degenerate vibrational eigenstates, as described by the Allen-Feldman equation. Recently, it has been shown that these two conduction mechanisms emerge from a Wigner transport equation, which unifies and extends the Peierls-Boltzmann and Allen-Feldman formulations, allowing to describe also complex crystals where particle-like and wave-like conduction mechanisms coexist. Here, we discuss the theoretical foundations of such transport equation as is derived from the Wigner phase-space formulation of quantum mechanics, elucidating how the interplay between disorder, anharmonicity, and the quantum Bose-Einstein statistics of atomic vibrations determines thermal conductivity. This Wigner formulation argues for a preferential phase convention for the dynamical matrix in the reciprocal Bloch representation and related off-diagonal velocity operator’s elements; such convention is the only one yielding a conductivity which is invariant with respect to the non-unique choice of the crystal’s unit cell and is size-consistent. We rationalize the conditions determining the crossover from particle-like to wave-like heat conduction, showing that phonons below the Ioffe-Regel limit ( i.e. with a mean free path shorter than the interatomic spacing) contribute to heat transport due to their wave-like capability to interfere and tunnel. Finally, we show that the present approach overcomes the failures of the Peierls-Boltzmann formulation for crystals with ultralow or glass-like thermal conductivity, with case studies of materials for thermal barrier coatings and thermoelectric energy conversion.
Giovanni Caldarelli, Michele Simoncelli, Nicola Marzari, Francesco Mauri, and Lara Benfatto
American Physical Society (APS)
Recent progress in understanding thermal transport in complex crystals has highlighted the promi-nent role of heat conduction mediated by interband tunneling processes, which emerge between overlapping phonon bands ( i.e . with energy differences smaller than their broadenings). These processes have recently been described in different ways, relying on the Wigner or Green-Kubo formalism, leading to apparently different results which question the definition of the heat-current operator. Here, we implement a full quantum approach based on the Kubo formula, elucidating analogies and differences with the recently introduced Wigner or Green-Kubo formulations, and extending the description of thermal transport to the overdamped regime of atomic vibrations, where the phonon quasiparticle picture breaks down. We rely on first-principles calculations on complex crystals with ultralow conductivity to compare numerically the thermal conductivity obtained within the aforementioned approaches, showing that at least in the quasiparticle regime the differences are negligible for practical applications.
Michele Simoncelli, Nicola Marzari, and Andrea Cepellotti
American Physical Society (APS)
Heat conduction in dielectric crystals originates from the propagation of atomic vibrations, whose microscopic dynamics is well described by the linearized phonon Boltzmann transport equation. Recently, it was shown that thermal conductivity can be resolved exactly and in a closed form as a sum over relaxons, $\\mathit{i.e.}$ collective phonon excitations that are the eigenvectors of Boltzmann equation's scattering matrix [Cepellotti and Marzari, PRX $\\mathbf{6}$ (2016)]. Relaxons have a well-defined parity, and only odd relaxons contribute to the thermal conductivity. Here, we show that the complementary set of even relaxons determines another quantity --- the thermal viscosity --- that enters into the description of heat transport, and is especially relevant in the hydrodynamic regime, where dissipation of crystal momentum by Umklapp scattering phases out. We also show how the thermal conductivity and viscosity parametrize two novel viscous heat equations --- two coupled equations for the temperature and drift-velocity fields --- which represent the thermal counterpart of the Navier-Stokes equations of hydrodynamics in the linear, laminar regime. These viscous heat equations are derived from a coarse-graining of the linearized Boltzmann transport equation for phonons, and encompass both limits of Fourier's law and of second sound, taking place, respectively, in the regimes of strong or weak momentum dissipation. Last, we introduce the Fourier deviation number as a descriptor that captures the deviations from Fourier's law due to hydrodynamic effects. We showcase these findings in a test case of a complex-shaped device made of graphite, obtaining a remarkable agreement with the very recent experimental demonstration of hydrodynamic transport in this material. The present findings also suggest that hydrodynamic behavior can appear at room temperature in micrometer-sized diamond crystals.
Michele Simoncelli, Nicola Marzari, and Francesco Mauri
Springer Science and Business Media LLC
Michele Simoncelli, Nidhal Ganfoud, Assane Sene, Matthieu Haefele, Barbara Daffos, Pierre-Louis Taberna, Mathieu Salanne, Patrice Simon, and Benjamin Rotenberg
American Physical Society (APS)
Capacitive mixing (CapMix) and capacitive deionization (CDI) are currently developed as alternatives to membrane-based processes to harvest blue energy—from salinity gradients between river and sea water— and to desalinate water—using charge-discharge cycles of capacitors. Nanoporous electrodes increase the contact area with the electrolyte and hence, in principle, also the performance of the process. However, models to design and optimize devices should be used with caution when the size of the pores becomes comparable to that of ions and water molecules. Here, we address this issue by simulating realistic capacitors based on aqueous electrolytes and nanoporous carbide-derived carbon (CDC) electrodes, accounting for both their complex structure and their polarization by the electrolyte under applied voltage. We compute the capacitance for two salt concentrations and validate our simulations by comparison with cyclic voltammetry experiments. We discuss the predictions of Debye-Huckel and Poisson-Boltzmann theories, as well as modified Donnan models, and we show that the latter can be parametrized using the molecular simulation results at high concentration. This then allows us to extrapolate the capacitance and salt adsorption capacity at lower concentrations, which cannot be simulated, finding a reasonable agreement with the experimental capacitance. We analyze the solvation of ions and their confinement within the electrodes—microscopic properties that are much more difficult to obtain experimentally than the electrochemical response but very important to understand the mechanisms at play. We finally discuss the implications of our findings for CapMix and CDI, both from the modeling point of view and from the use of CDCs in these contexts.