Quantum ESPRESSO is an integrated suite of open-source computer codes for quantum simulations of materials using state-of-the art electronic-structure techniques, based on density-functional theory, density-functional perturbation theory, and many-body perturbation theory, within the plane-wave pseudo-potential and projector-augmented-wave approaches. Quantum ESPRESSO owes its popularity to the wide variety of properties and processes it allows to simulate, to its performance on an increasingly broad array of hardware architectures, and to a community of researchers that rely on its capabilities as a core open-source development platform to implement theirs ideas. In this paper we describe recent extensions and improvements, covering new methodologies and property calculators, improved parallelization, code modularization, and extended interoperability both within the distribution and with external software.

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.

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.

We present a first-principles theoretical approach for evaluating the lattice thermal conductivity based on the exact solution of the Boltzmann transport equation. We use the variational principle and the conjugate gradient scheme, which provide us with an algorithm faster than the one previously used in literature and able to always converge to the exact solution. Three-phonon normal and umklapp collision, isotope scattering and border effects are rigorously treated in the calculation. Good agreement with experimental data for diamond is found. Moreover we show that by growing more enriched diamond samples it is possible to achieve values of thermal conductivity up to three times larger than the commonly observed in isotopically enriched diamond samples with 99.93% C 12 and 0.07 C 13 .

Thermal transport is a key feature for the operation of phase change memory devices which rest on a fast and reversible transformation between the crystalline and amorphous phases of chalcogenide alloys upon Joule heating. In this paper we report on the ab initio calculations of bulk thermal conductivity of the prototypical phase change compounds Ge 2 Sb 2 Te 5 and GeTe in their crystalline form. The related Sb 2 Te 3 compound is also investigated for the sake of comparison. Thermal conductivity is obtained from the solution of the Boltzmann transport equation with phonon scattering rates computed within density functional perturbation theory. The calculations show that the large spread in the experimental data on the lattice thermal conductivity of GeTe is due to a variable content of Ge vacancies which at concentrations realized experimentally can halve the bulk thermal conductivity with respect to the ideal crystal. We show that the very low thermal conductivity of hexagonal Ge 2 Sb 2 Te 5 of about 0.45 W m −1 K −1 measured experimentally is also resulting from disorder in the form of a random distribution of Ge/Sb atoms in one sublattice.

International audienceA promising route towards nanodevice applications relies on the association of graphene and transition metal dichalcogenides with hexagonal boron nitride (h-BN). Due to its insulating nature, h-BN has emerged as a natural substrate and gate dielectric for graphene-based electronic devices. However, some fundamental properties of bulk h-BN remain obscure. For example, the band structure and the position of the Fermi level have not been experimentally resolved. Here, we report a direct observation of parabolic dispersions of h-BN crystals using high-resolution angle-resolved photoemission spectroscopy (ARPES). We find that h-BN exfoliation on epitaxial graphene enables overcoming the technical difficulties of using ARPES with insulating materials. We show trigonal warping of the intensity maps at constant energy. The valence-band maxima are located around the K points, 2.5 eV below the Fermi level, thus confirming the residual p-type character of typical h-BN

Green tea catechins are polyphenols which are believed to provide health benefits; they are marketed as health supplements and are studied for their potential effects on a variety of medical conditions. However, their mechanisms of action and interaction with the environment at the molecular level are still not well-understood. Here, by means of atomistic simulations, we explore the structural properties of four green tea catechins, in the gas phase and water solution: specifically, (-)-epigallocatechin-3-gallate, which is the most abundant, (-)-epicatechin-3-gallate, (-)-epigallocatechin-3-O-(3-O-methyl)-gallate, and (-)-epigallocatechin. We characterize the free energy conformational landscapes of these catechins at ambient conditions, as a function of the torsional degrees of freedom of the pholyphenolic rings, determining the stable conformers and their connections. We show that these free energy landscapes are only subtly influenced by the interactions with the solvent and by the structural details of the polyphenolic rings. However, the number and position of the hydroxyl groups (or their sustituents) and the presence/absence of the galloyl moiety have significant impact on the selected catechin solvation shells and hydrogen bond capabilities, which are ultimately linked to their ability to interact with and affect the biological environment.

Bismuth is one of the rare materials in which second sound has been experimentally observed. Our exact calculations of thermal transport with the Boltzmann equation predict the occurrence of this Poiseuille phonon flow between ≈1.5 and ≈3.5  K, in a sample size of 3.86 and 9.06 mm, consistent with the experimental observations. Hydrodynamic heat flow characteristics are given for any temperature: heat wave propagation length, drift velocity, and Knudsen number. We discuss a gedanken experiment allowing us to assess the presence of a hydrodynamic regime in any bulk material.