Andrea Cepellotti

Thermal conductivity of graphene and graphite: collective excitations and mean free paths.

We characterize the thermal conductivity of graphite, monolayer graphene, graphane, fluorographane, and bilayer graphene, solving exactly the Boltzmann transport equation for phonons, with phonon-phonon collision rates obtained from density functional perturbation theory. For graphite, the results are found to be in excellent agreement with experiments; notably, the thermal conductivity is 1 order of magnitude larger than what found by solving the Boltzmann equation in the single mode approximation, commonly used to describe heat transport. For graphene, we point out that a meaningful value of intrinsic thermal conductivity at room temperature can be obtained only for sample sizes of the order of 1 mm, something not considered previously. This unusual requirement is because collective phonon excitations, and not single phonons, are the main heat carriers in these materials; these excitations are characterized by mean free paths of the order of hundreds of micrometers. As a result, even Fouriers law becomes questionable in typical sample sizes, because its statistical nature makes it applicable only in the thermodynamic limit to systems larger than a few mean free paths. Finally, we discuss the effects of isotopic disorder, strain, and chemical functionalization on thermal performance. Only chemical functionalization is found to play an important role, decreasing the conductivity by a factor of 2 in hydrogenated graphene, and by 1 order of magnitude in fluorogenated graphene.

Phonon hydrodynamics in two-dimensional materials

The conduction of heat in two dimensions displays a wealth of fascinating phenomena of key relevance to the scientific understanding and technological applications of graphene and related materials. Here, we use density-functional perturbation theory and an exact, variational solution of the Boltzmann transport equation to study fully from first-principles phonon transport and heat conductivity in graphene, boron nitride, molybdenum disulphide and the functionalized derivatives graphane and fluorographene. In all these materials, and at variance with typical three-dimensional solids, normal processes keep dominating over Umklapp scattering well-above cryogenic conditions, extending to room temperature and more. As a result, novel regimes emerge, with Poiseuille and Ziman hydrodynamics, hitherto typically confined to ultra-low temperatures, characterizing transport at ordinary conditions. Most remarkably, several of these two-dimensional materials admit wave-like heat diffusion, with second sound present at room temperature and above in graphene, boron nitride and graphane.

Two-dimensional materials from high-throughput computational exfoliation of experimentally known compounds

Two-dimensional (2D) materials have emerged as promising candidates for next-generation electronic and optoelectronic applications. Yet, only a few dozen 2D materials have been successfully synthesized or exfoliated. Here, we search for 2D materials that can be easily exfoliated from their parent compounds. Starting from 108,423 unique, experimentally known 3D compounds, we identify a subset of 5,619 compounds that appear layered according to robust geometric and bonding criteria. High-throughput calculations using van der Waals density functional theory, validated against experimental structural data and calculated random phase approximation binding energies, further allowed the identification of 1,825 compounds that are either easily or potentially exfoliable. In particular, the subset of 1,036 easily exfoliable cases provides novel structural prototypes and simple ternary compounds as well as a large portfolio of materials to search from for optimal properties. For a subset of 258 compounds, we explore vibrational, electronic, magnetic and topological properties, identifying 56 ferromagnetic and antiferromagnetic systems, including half-metals and half-semiconductors.The largest available database of potentially exfoliable 2D materials has been obtained via high-throughput calculations using van der Waals density functional theory.

Thermal Transport in Crystals as a Kinetic Theory of Relaxons

Thermal conductivity in dielectric crystals is the result of the relaxation of lattice vibrations described by the phonon Boltzmann transport equation. Remarkably, an exact microscopic definition of the heat carriers and their relaxation times is still missing: phonons, typically regarded as the relevant excitations for thermal transport, cannot be identified as the heat carriers when most scattering events conserve momentum and do not dissipate heat flux. This is the case for two-dimensional or layered materials at room temperature, or three-dimensional crystals at cryogenic temperatures. In this work we show that the eigenvectors of the scattering matrix in the Boltzmann equation define collective phonon excitations, termed here relaxons. These excitations have well defined relaxation times, directly related to heat flux dissipation, and provide an exact description of thermal transport as a kinetic theory of the relaxon gas. We show why Matthiessens rule is violated, and construct a procedure for obtaining the mean free paths and relaxation times of the relaxons. These considerations are general, and would apply also to other semiclassical transport models, such as the electronic Boltzmann equation. For heat transport, they remain relevant even in conventional crystals like silicon, but are of the utmost importance in the case of two-dimensional materials, where they can revise by several orders of magnitude the relevant time- and length-scales for thermal transport in the hydrodynamic regime.

Boltzmann Transport in Nanostructures as a Friction Effect

Surface scattering is the key limiting factor to thermal transport in dielectric crystals as the length scales are reduced or when temperature is lowered. To explain this phenomenon, it is commonly assumed that the mean free paths of heat carriers are bound by the crystal size and that thermal conductivity is reduced in a manner proportional to such mean free paths. We show here that these conclusions rely on simplifying assumptions and approximated transport models. Instead, starting from the linearized Boltzmann transport equation in the relaxon basis, we show how the problem can be reduced to a set of decoupled linear differential equations. Then, the heat flow can be interpreted as a hydrodynamic phenomenon with the relaxon gas being slowed down in proximity of a surface by friction effects, similar to the flux of a viscous fluid in a pipe. As an example, we study a ribbon and a trench of monolayer molybdenum disulfide, describing the procedure to reconstruct the temperature and thermal conductivity profile in the sample interior and showing how to estimate the effect of nanostructuring. The approach is general and could be extended to other transport carriers, such as electrons, or extended to materials of higher dimensionality and to different geometries, such as thin films.

Transport waves as crystal excitations

We introduce the concept of transport waves by showing that the linearized Boltzmann transport equation admits excitations in the form of waves that have well defined dispersion relations and decay times. Crucially, these waves do not represent single-particle excitations, but are collective excitations of the equilibrium distribution functions. We study in detail the case of thermal transport, where relaxons are found in the long-wavelength limit, and second sound is reinterpreted as the excitation of one or several temperature waves at finite frequencies. Graphene is studied numerically, finding decay times of the order of microseconds. The derivation, obtained by a spectral representation of the Boltzmann equation, holds in principle for any crystal or semiclassical transport theory and is particularly relevant when transport takes place in the hydrodynamic regime.

Rippling ultrafast dynamics of suspended 2D monolayers, graphene

Significance Rippling is an intrinsic feature of 2D materials, responsible for their structural stability, transport properties, and electron–hole charge redistribution. Modulating these ripples in a controlled manner not only provides a better understanding of their structural properties, but also has potential impact for applications. Here, we examine graphene monolayer as a prototypical 2D material. An ultrafast attenuation of the ripples intrinsically present in the graphene plane is followed by a significant enhancement of the rippling effect on a longer time scale, as driven by the successive excitation of in-plane and out-of-plane phonon modes. The methodology described is of a general nature and is suitable for the investigation of other 2D materials, where we expect to observe similar rippling effects. Here, using ultrafast electron crystallography (UEC), we report the observation of rippling dynamics in suspended monolayer graphene, the prototypical and most-studied 2D material. The high scattering cross-section for electron/matter interaction, the atomic-scale spatial resolution, and the ultrafast temporal resolution of UEC represent the key elements that make this technique a unique tool for the dynamic investigation of 2D materials, and nanostructures in general. We find that, at early time after the ultrafast optical excitation, graphene undergoes a lattice expansion on a time scale of 5 ps, which is due to the excitation of short-wavelength in-plane acoustic phonon modes that stretch the graphene plane. On a longer time scale, a slower thermal contraction with a time constant of 50 ps is observed and associated with the excitation of out-of-plane phonon modes, which drive the lattice toward thermal equilibrium with the well-known negative thermal expansion coefficient of graphene. From our results and first-principles lattice dynamics and out-of-equilibrium relaxation calculations, we quantitatively elucidate the deformation dynamics of the graphene unit cell.

A posteriori metadata from automated provenance tracking: integration of AiiDA and TCOD

In order to make results of computational scientific research findable, accessible, interoperable and re-usable, it is necessary to decorate them with standardised metadata. However, there are a number of technical and practical challenges that make this process difficult to achieve in practice. Here the implementation of a protocol is presented to tag crystal structures with their computed properties, without the need of human intervention to curate the data. This protocol leverages the capabilities of AiiDA, an open-source platform to manage and automate scientific computational workflows, and the TCOD, an open-access database storing computed materials properties using a well-defined and exhaustive ontology. Based on these, the complete procedure to deposit computed data in the TCOD database is automated. All relevant metadata are extracted from the full provenance information that AiiDA tracks and stores automatically while managing the calculations. Such a protocol also enables reproducibility of scientific data in the field of computational materials science. As a proof of concept, the AiiDA–TCOD interface is used to deposit 170 theoretical structures together with their computed properties and their full provenance graphs, consisting in over 4600 AiiDA nodes.