AiiDA 1.0, a scalable computational infrastructure for automated reproducible workflows and data provenance.

The ever-growing availability of computing power and the sustained development of advanced computational methods have contributed much to recent scientific progress. These developments present new challenges driven by the sheer amount of calculations and data to manage. Next-generation exascale supercomputers will harden these challenges, such that automated and scalable solutions become crucial. In recent years, we have been developing AiiDA (aiida.net), a robust open-source high-throughput infrastructure addressing the challenges arising from the needs of automated workflow management and data provenance recording. Here, we introduce developments and capabilities required to reach sustained performance, with AiiDA supporting throughputs of tens of thousands processes/hour, while automatically preserving and storing the full data provenance in a relational database making it queryable and traversable, thus enabling high-performance data analytics. AiiDA's workflow language provides advanced automation, error handling features and a flexible plugin model to allow interfacing with external simulation software. The associated plugin registry enables seamless sharing of extensions, empowering a vibrant user community dedicated to making simulations more robust, user-friendly and reproducible.

Optomechanical Measurement of Thermal Transport in Two-Dimensional MoSe2 Lattices.

Nanomechanical resonators have emerged as sensors with exceptional sensitivities. These sensing capabilities open new possibilities in the studies of the thermodynamic properties in condensed matter. Here, we use mechanical sensing as a novel approach to measure the thermal properties of low-dimensional materials. We measure the temperature dependence of both the thermal conductivity and the specific heat capacity of a transition metal dichalcogenide monolayer down to cryogenic temperature, something that has not been achieved thus far with a single nanoscale object. These measurements show how heat is transported by phonons in two-dimensional systems. Both the thermal conductivity and the specific heat capacity measurements are consistent with predictions based on first-principles.

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.

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.

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.

Rippling ultrafast dynamics of suspended 2D monolayers, graphene.

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.

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.

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 Fourier's 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.