Nicola Marzari

QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials

QUANTUM ESPRESSO is an integrated suite of computer codes for electronic-structure calculations and materials modeling, based on density-functional theory, plane waves, and pseudopotentials (norm-conserving, ultrasoft, and projector-augmented wave). The acronym ESPRESSO stands for opEn Source Package for Research in Electronic Structure, Simulation, and Optimization. It is freely available to researchers around the world under the terms of the GNU General Public License. QUANTUM ESPRESSO builds upon newly-restructured electronic-structure codes that have been developed and tested by some of the original authors of novel electronic-structure algorithms and applied in the last twenty years by some of the leading materials modeling groups worldwide. Innovation and efficiency are still its main focus, with special attention paid to massively parallel architectures, and a great effort being devoted to user friendliness. QUANTUM ESPRESSO is evolving towards a distribution of independent and interoperable codes in the spirit of an open-source project, where researchers active in the field of electronic-structure calculations are encouraged to participate in the project by contributing their own codes or by implementing their own ideas into existing codes.

First-principles determination of the structural, vibrational and thermodynamic properties of diamond, graphite, and derivatives

The structural, dynamical, and thermodynamic properties of diamond, graphite and layered derivatives (graphene, rhombohedral graphite) are computed using a combination of density-functional theory total-energy calculations and density-functional perturbation theory lattice dynamics in the generalized gradient approximation. Overall, very good agreement is found for the structural properties and phonon dispersions, with the exception of the c/a ratio in graphite and the associated elastic constants and phonon dispersions. Both the C-33 elastic constant and the F to A phonon dispersions are brought to close agreement with available data once the experimental c/a is chosen for the calculations. The vibrational free energy and the thermal expansion, the temperature dependence of the elastic moduli and the specific heat are calculated using the quasiharmonic approximation. Graphite shows a distinctive in-plane negative thermal-expansion coefficient that reaches its lowest value around room temperature, in very good agreement with experiments. Thermal contraction in graphene is found to be three times as large; in both cases, bending acoustic modes are shown to be responsible for the contraction, in a direct manifestation of the membrane effect predicted by Lifshitz over 50 years ago. Stacking directly affects the bending modes, explaining the large numerical difference between the thermal-contraction coefficients in graphite and graphene, notwithstanding their common physical origin.

Thermal contraction and disordering of the Al(110) surface

Al(110) has been studied for temperatures up to 900 K via ensemble density-functional molecular dynamics. The strong anharmonicity displayed by this surface results in a negative coefficient of thermal expansion, where the first interlayer distance decreases with increasing temperature. Very shallow channels of oscillation for the second-layer atoms in the direction perpendicular to the surface support this anomalous contraction, and provide a novel mechanism for the formation of adatom-vacancy pairs, preliminary to the disordering and premelting transition. Such characteristic behaviour originates in the free-electron-gas bonding at a loosely packed surface. [S0031-9007(99)08925-5].

The shear mode of multilayer graphene

The quest for materials capable of realizing the next generation of electronic and photonic devices continues to fuel research on the electronic, optical and vibrational properties of graphene. Few-layer graphene (FLG) flakes with less than ten layers each show a distinctive band structure. Thus, there is an increasing interest in the physics and applications of FLGs. Raman spectroscopy is one of the most useful and versatile tools to probe graphene samples. Here, we uncover the interlayer shear mode of FLGs, ranging from bilayer graphene (BLG) to bulk graphite, and suggest that the corresponding Raman peak measures the interlayer coupling. This peak scales from ~43 cm(-1) in bulk graphite to ~31 cm(-1) in BLG. Its low energy makes it sensitive to near-Dirac point quasiparticles. Similar shear modes are expected in all layered materials, providing a direct probe of interlayer interactions.

Density functional theory in transition-metal chemistry: a self-consistent Hubbard U approach.

Transition-metal centers are the active sites for a broad variety of biological and inorganic chemical reactions. Notwithstanding this central importance, density-functional theory calculations based on generalized-gradient approximations often fail to describe energetics, multiplet structures, reaction barriers, and geometries around the active sites. We suggest here an alternative approach, derived from the Hubbard U correction to solid-state problems, that provides an excellent agreement with correlated-electron quantum chemistry calculations in test cases that range from the ground state of Fe2 and Fe2- to the addition elimination of molecular hydrogen on FeO+. The Hubbard U is determined with a novel self-consistent procedure based on a linear-response approach.

Large-Area Epitaxial Monolayer MoS2

Two-dimensional semiconductors such as MoS2 are an emerging material family with wide-ranging potential applications in electronics, optoelectronics, and energy harvesting. Large-area growth methods are needed to open the way to applications. Control over lattice orientation during growth remains a challenge. This is needed to minimize or even avoid the formation of grain boundaries, detrimental to electrical, optical, and mechanical properties of MoS2 and other 2D semiconductors. Here, we report on the growth of high-quality monolayer MoS2 with control over lattice orientation. We show that the monolayer film is composed of coalescing single islands with limited numbers of lattice orientation due to an epitaxial growth mechanism. Optical absorbance spectra acquired over large areas show significant absorbance in the high-energy part of the spectrum, indicating that MoS2 could also be interesting for harvesting this region of the solar spectrum and fabrication of UV-sensitive photodetectors. Even though the interaction between the growth substrate and MoS2 is strong enough to induce lattice alignment via van der Waals interaction, we can easily transfer the grown material and fabricate devices. Local potential mapping along channels in field-effect transistors shows that the single-crystal MoS2 grains in our film are well connected, with interfaces that do not degrade the electrical conductivity. This is also confirmed by the relatively large and length-independent mobility in devices with a channel length reaching 80 μm.

Reproducibility in density functional theory calculations of solids

A comparison of DFT methods Density functional theory (DFT) is now routinely used for simulating material properties. Many software packages are available, which makes it challenging to know which are the best to use for a specific calculation. Lejaeghere et al. compared the calculated values for the equation of states for 71 elemental crystals from 15 different widely used DFT codes employing 40 different potentials (see the Perspective by Skylaris). Although there were variations in the calculated values, most recent codes and methods converged toward a single value, with errors comparable to those of experiment. Science, this issue p. 10.1126/science.aad3000; see also p. 1394 A survey of recent density functional theory methods shows a convergence to more accurate property calculations. [Also see Perspective by Skylaris] INTRODUCTION The reproducibility of results is one of the underlying principles of science. An observation can only be accepted by the scientific community when it can be confirmed by independent studies. However, reproducibility does not come easily. Recent works have painfully exposed cases where previous conclusions were not upheld. The scrutiny of the scientific community has also turned to research involving computer programs, finding that reproducibility depends more strongly on implementation than commonly thought. These problems are especially relevant for property predictions of crystals and molecules, which hinge on precise computer implementations of the governing equation of quantum physics. RATIONALE This work focuses on density functional theory (DFT), a particularly popular quantum method for both academic and industrial applications. More than 15,000 DFT papers are published each year, and DFT is now increasingly used in an automated fashion to build large databases or apply multiscale techniques with limited human supervision. Therefore, the reproducibility of DFT results underlies the scientific credibility of a substantial fraction of current work in the natural and engineering sciences. A plethora of DFT computer codes are available, many of them differing considerably in their details of implementation, and each yielding a certain “precision” relative to other codes. How is one to decide for more than a few simple cases which code predicts the correct result, and which does not? We devised a procedure to assess the precision of DFT methods and used this to demonstrate reproducibility among many of the most widely used DFT codes. The essential part of this assessment is a pairwise comparison of a wide range of methods with respect to their predictions of the equations of state of the elemental crystals. This effort required the combined expertise of a large group of code developers and expert users. RESULTS We calculated equation-of-state data for four classes of DFT implementations, totaling 40 methods. Most codes agree very well, with pairwise differences that are comparable to those between different high-precision experiments. Even in the case of pseudization approaches, which largely depend on the atomic potentials used, a similar precision can be obtained as when using the full potential. The remaining deviations are due to subtle effects, such as specific numerical implementations or the treatment of relativistic terms. CONCLUSION Our work demonstrates that the precision of DFT implementations can be determined, even in the absence of one absolute reference code. Although this was not the case 5 to 10 years ago, most of the commonly used codes and methods are now found to predict essentially identical results. The established precision of DFT codes not only ensures the reproducibility of DFT predictions but also puts several past and future developments on a firmer footing. Any newly developed methodology can now be tested against the benchmark to verify whether it reaches the same level of precision. New DFT applications can be shown to have used a sufficiently precise method. Moreover, high-precision DFT calculations are essential for developing improvements to DFT methodology, such as new density functionals, which may further increase the predictive power of the simulations. Recent DFT methods yield reproducible results. Whereas older DFT implementations predict different values (red darts), codes have now evolved to mutual agreement (green darts). The scoreboard illustrates the good pairwise agreement of four classes of DFT implementations (horizontal direction) with all-electron results (vertical direction). Each number reflects the average difference between the equations of state for a given pair of methods, with the green-to-red color scheme showing the range from the best to the poorest agreement. The widespread popularity of density functional theory has given rise to an extensive range of dedicated codes for predicting molecular and crystalline properties. However, each code implements the formalism in a different way, raising questions about the reproducibility of such predictions. We report the results of a community-wide effort that compared 15 solid-state codes, using 40 different potentials or basis set types, to assess the quality of the Perdew-Burke-Ernzerhof equations of state for 71 elemental crystals. We conclude that predictions from recent codes and pseudopotentials agree very well, with pairwise differences that are comparable to those between different high-precision experiments. Older methods, however, have less precise agreement. Our benchmark provides a framework for users and developers to document the precision of new applications and methodological improvements.

Self-assembled quantum dots in a nanowire system for quantum photonics

Quantum dots embedded within nanowires represent one of the most promising technologies for applications in quantum photonics. Whereas the top-down fabrication of such structures remains a technological challenge, their bottom-up fabrication through self-assembly is a potentially more powerful strategy. However, present approaches often yield quantum dots with large optical linewidths, making reproducibility of their physical properties difficult. We present a versatile quantum-dot-in-nanowire system that reproducibly self-assembles in core-shell GaAs/AlGaAs nanowires. The quantum dots form at the apex of a GaAs/AlGaAs interface, are highly stable, and can be positioned with nanometre precision relative to the nanowire centre. Unusually, their emission is blue-shifted relative to the lowest energy continuum states of the GaAs core. Large-scale electronic structure calculations show that the origin of the optical transitions lies in quantum confinement due to Al-rich barriers. By emitting in the red and self-assembling on silicon substrates, these quantum dots could therefore become building blocks for solid-state lighting devices and third-generation solar cells.

Maximally-localized Wannier functions for disordered systems: Application to amorphous silicon

We use the maximally-localized Wannier function method to study bonding properties in amorphous silicon. This study represents, to our knowledge, the first application of the Wannier-function analysis to a disordered system. Our results show that, in the presence of disorder, this method is extremely helpful in providing an unambiguous picture of the bond distribution. In particular, defect configurations can be studied and characterized with a novel degree of accuracy that was not available before.

Surface energies, work functions, and surface relaxations of low index metallic surfaces from first principles

We study the relaxations, surface energies, and work functions of low-index metallic surfaces using pseudopotential plane-wave density-functional calculations within the generalized gradient approximation. We study here the (100), (110), and (111) surfaces of Al, Pd, Pt, and Au and the (0001) surface of Ti, chosen for their use as contact or lead materials in nanoscale devices. We consider clean, mostly nonreconstructed surfaces in the slab-supercell approximation. Particular attention is paid to the convergence of these quantities with respect to slab thickness; furthermore, different methodologies for the calculation of work functions and surfaces energies are compared. We find that the use of bulk references for calculations of surface energies and work functions can be detrimental to convergence unless numerical grids are closely matched, especially when surface relaxations are being considered. Our results and comparison show that calculated values often do not quantitatively match experimental values. This may be understandable for the surface relaxations and surface energies, where experimental values can have large error but even for the work functions, neither local nor semilocal functionals emerge as an accurate choice for every case.