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 solvation model proposed by Fattebert and Gygi [1] and Scherlis et al. [2] is reformulated, overcoming some of the numerical limitations encountered and extending its range of applicability.We first recast the problem in terms of induced polarization charges that act as a direct mapping of the self-consistent continuum dielectric; this allows to define a functional form for the dielectric that is well behaved both in the high-density region of the nuclear charges and in the low-density region where the electronic wavefunctions decay into the solvent. Second, we outline an iterative procedure to solve the Poisson equation for the quantum fragment embedded in the solvent that does not require multi-grid algorithms, is trivially parallel, and can be applied to any Bravais crystallographic system. Last, we capture some of the non-electrostatic or cavitation terms via a combined use of the quantum volume and quantum surface kcal/mol, whereby larger discrepancies are mostly limited to self-dissociating species and strong hydrogen-bond forming compounds.

In catalysis science stability is as crucial as activity and selectivity. Understanding the degradation pathways occurring during operation and developing mitigation strategies will eventually improve catalyst design, thus facilitating the translation of basic science to technological applications. Herein, we reveal the unique and general degradation mechanism of metallic nanocatalysts during electrochemical CO2 reduction, exemplified by different sized copper nanocubes. We follow their morphological evolution during operation and correlate it with the electrocatalytic performance. In contrast with the most common coalescence and dissolution/precipitation mechanisms, we find a potential-driven nanoclustering to be the predominant degradation pathway. Grand-potential density functional theory calculations confirm the role of the negative potential applied to reduce CO2 as the main driving force for the clustering. This study offers a novel outlook on future investigations of stability and degradation reaction mechanisms of nanocatalysts in electrochemical CO2 reduction and, more generally, in electroreduction reactions.

We investigated the cathodic and anodic limits of six room-temperature ionic liquids (ILs) formed from a combination of two common cations-1-butyl-3-methylimidazolium (BMIM) and N,N-propylmethylpyrrolidinium (P13)-and three common anions-PF 6 , BF 4 and bis(trifluoromethylsulfonyl)imide (TFSI)-using an approach that combines molecular dynamics (MD) simulations and density functional theory (DFT) calculations. All inter-ion interactions were taken into account by explicitly modeling the entire liquid structure using classical MD, followed by DFT computations of the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energies. The relative cathodic and anodic limits of BMIM PF 6 , BMIM BF 4 , BMIM TFSI and P13 TFSI obtained from our approach are * To whom correspondence should be addressed

We discuss grand canonical simulations based on density-functional theory to study the thermodynamic properties of electrochemical interfaces of metallic electrodes in aqueous environments. Water is represented using implicit solvation, here via the self-consistent continuum solvation (SCCS) model, providing a charge-density dependent dielectric boundary. The electrochemical double layer is accounted for in terms of a phenomenological continuum description. It is shown that the experimental potentials of zero charge and interfacial capacitances can be reproduced for an optimized SCCS parameter set [ρ min = 0.0013, ρ max = 0.010 25]. By performing a detailed derivation and analysis of the interface energetics for selected electrochemical systems, we are able to relate the widely used approach of the computational hydrogen electrode (CHE) to a general grand canonical description of electrified interfaces. In particular, charge-neutral CHE results are shown to be an upper-boundary estimate for the grand canonical interfacial free energies. In order to demonstrate the differences between the CHE and full grand canonical calculations, we study the pristine (100), (110), and (111) surfaces for Pt, Au, Cu, and Ag, and H or Cl electrosorbed on Pt. The calculations support the known surface reconstructions in the aqueous solution for Pt and Au. Furthermore, the predicted potential-pH dependence of proton coverage, surface charge, and interfacial pseudocapacitance for Pt is found to be in close agreement with experimental or other theoretical data as well as the predicted equilibrium shapes for Pt nanoparticles. Finally, Cl is found to interact more strongly than H with the interfacial fields, leading to significantly altered interface energetics and structure upon explicit application of an electrode potential. This work underscores the strengths and eventual limits of the CHE approach and might guide further understanding of the thermodynamics of electrified interfaces.

The CO2 electro-reduction reaction (CORR) is a promising avenue to convert greenhouse gases into high-value fuels and chemicals, in addition to being an attractive method for storing intermittent renewable energy. Although polycrystalline Cu surfaces have long known to be unique in their capabilities of catalyzing the conversion of CO2 to higher-order C1 and C2 fuels, such as hydrocarbons (CH4, C2H4 etc.) and alcohols (CH3OH, C2H5OH), product selectivity remains a challenge. Rational design of more selective catalysts would greatly benefit from a mechanistic understanding of the complex, multi-proton and multi-electron conversion of CO2. In this study, we select three metal catalysts (Pt, Au, Cu) and apply in situ surface enhanced infrared absorption spectroscopy (SEIRAS) and ambient-pressure X-ray photoelectron spectroscopy (APXPS), coupled to density-functional theory (DFT) calculations, to get insight into the reaction pathway for the CORR. We present a comprehensive reaction mechanism for the CORR, and show that the preferential reaction pathway can be rationalized in terms of metal-carbon (M-C) and metaloxygen (M-O) affinity. We show that the final products are determined by the configuration of the initial intermediates, C-bound and O-bound, which can be obtained from CO2 and (H)CO3, respectively. C1 hydrocarbons are produced via OCH3,ad intermediates obtained from O-bound CO3,ad and require a catalyst with relatively high affinity for O-bound intermediates. Additionally, C2 hydrocarbon formation is suggested to result from the C-C coupling between C-bound COad and (H)COad, which requires an optimal affinity for the C-bound species, so that (H)COad can be further reduced without poisoning the catalyst surface. It is suggested that the formation of C1 alcohols (CH3OH) is the most challenging process to optimize, as stabilization of the O-bound species would both accelerate the formation of key-intermediates (OCH3,ad) but also simultaneously inhibit their desorption from the catalyst surface. Our findings pave the way

We present an implicit solvation approach where the interface between the quantum-mechanical solute and the surrounding environment is described by a fully continuous permittivity built up with atomic-centered "soft" spheres. This approach combines many of the advantages of the self-consistent continuum solvation model in handling solutes and surfaces in contact with complex dielectric environments or electrolytes in electronic-structure calculations. In addition it is able to describe accurately both neutral and charged systems. The continuous function, describing the variation of the permittivity, allows to compute analytically the nonelectrostatic contributions to the solvation free energy that are described in terms of the quantum surface. The whole methodology is computationally stable, provides consistent energies and forces, and keeps the computational efforts and runtimes comparable to those of standard vacuum calculations. The capabilitiy to treat arbitrary molecular or slab-like geometries as well as charged molecules is key to tackle electrolytes within mixed explicit/implicit frameworks. We show that, with given, fixed atomic radii, two parameters are sufficient to give a mean absolute error of only 1.12 kcal/mol with respect to the experimental aqueous solvation energies for a set of 274 neutral solutes. For charged systems, the same set of parameters provides solvation energies for a set of 60 anions and 52 cations with an error of 2.96 and 2.13 kcal/mol, respectively, improving upon previous literature values. To tackle elements not present in most solvation databases, a new benchmark scheme on wettability and contact angles is proposed for solid-liquid interfaces and applied to the investigation of the stable terminations of a CdS (112̅0) surface in an electrochemical medium.

The recently developed self-consistent continuum solvation model (SCCS) [O. Andreussi, I. Dabo, and N. Marzari, J. Chem. Phys. 136, 064102 (2012)] is applied here to charged species in aqueous solutions. Describing ions in solution represents a great challenge because of the large electrostatic interactions between the solute and the solvent. The SCCS model is tested over 106 monocharged species, both cations and anions, and we demonstrate its flexibility, notwithstanding its much reduced set of parameters, to describe charged species in solution. Remarkably low mean absolute errors are obtained with values of 2.27 and 5.54 kcal/mol for cations and anions, respectively. These results are comparable or better than the state of the art to describe solvation of charged species in water. Finally, differences of behavior between cations and anions are discussed.