Guido Fratesi

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.

Advanced capabilities for materials modelling with Quantum ESPRESSO

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.Quantum EXPRESSO 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 pseudopotential and projector-augmented-wave approaches. Quantum EXPRESSO 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 their 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.

Analysis of methane-to-methanol conversion on clean and defective Rh surfaces

We investigate by density-functional theory simulations several elementary reactions associated to direct methane-to-methanol conversion on clean Rh(111) surfaces and on Rh adatoms on Rh(111). Energy barriers and reaction paths have been determined by the nudged elastic band method. The rate-limiting step in the process, C-O bond formation, has higher activation energy than the one for complete methane dehydrogenation. Our analysis enables us to understand the effect of defects on the reactivity and rules out Rh as candidate catalyst for methanol synthesis.

Ab initio calculation of core-valence-valence Auger spectra in closed shell systems

We propose an ab initio method to evaluate the core-valence-valence (CVV) Auger spectrum of systems with filled valence bands. The method is based on the Cini-Sawatzky theory, and aims at estimating the parameters by first-principles calculations in the framework of density-functional theory (DFT). Photoemission energies and the interaction energy for the two holes in the final state are evaluated by performing DFT simulations for the system with varied population of electronic levels. Transition matrix elements are taken from atomic results. The approach takes into account the non-sphericity of the density of states of the emitting atom, spin-orbit interaction in core and valence, and non quadratic terms in the total energy expansion with respect to fractional occupation numbers. It is tested on two benchmark systems, Zn and Cu metals, leading in both cases to L23M45M45 Auger peaks within 2 eV from the experimental ones. Detailed analysis is presented on the relative weight of the various contributions considered in our method, providing the basis for future development. Especially problematic is the evaluation of the hole-hole interaction for systems with broad valence bands: our method underestimates its value in Cu, while we obtain excellent results for this quantity in Zn.

High resolution NEXAFS of perylene and PTCDI: a surface science approach to molecular orbital analysis

We made use of synchrotron radiation to perform near edge X-ray absorption fine structure spectroscopy, NEXAFS, at the carbon K-edge of perylene and perylene-tetracarboxylic-diimide, PTCDI. Reference spectra measured for isolated molecules in the gas phase are compared with polarization dependent NEXAFS spectra measured on highly oriented thin films in order to study the symmetry of the molecular orbitals. The molecular overlayers are grown onto the rutile TiO2(110) surface for which the large anisotropic corrugation effectively drives the molecular orientation, while its dielectric nature prevents the rehybridization of the molecular orbitals. We employed density functional theory, DFT, calculations to disentangle the contribution of specific carbon atoms to the molecular density of states. Numerical simulations correctly predict the observed NEXAFS azimuthal dichroism of the σ* resonances above the ionization threshold, from which we determine the full geometric orientation of the overlayer molecules. A discrepancy observed for the spectral contribution of the imide carbon atom to the calculated unoccupied molecular orbitals has been explained in terms of initial state effects, as determined by Hartree-Fock corrections and in full agreement with the corresponding shift of the C 1s core level measured by X-ray photoelectron spectroscopy, XPS.

Electronic States of Silicene Allotropes on Ag(111)

Silicene, a honeycomb lattice of silicon, presents a particular case of allotropism on Ag(111). Silicene forms multiple structures with alike in-plane geometry but different out-of-plane atomic buckling and registry to the substrate. Angle-resolved photoemission and first-principles calculations show that these silicene structures, with (4×4), (√13×√13)R13.9°, and (2√3×2√3)R30° lattice periodicity, display similar electronic bands despite the structural differences. In all cases the interaction with the substrate modifies the electronic states, which significantly differ from those of free-standing silicene. Complex photoemission patterns arise from surface umklapp processes, varying according to the periodicity of the silicene allotropes.

Theoretical approaches in adsorption: alkali adatom investigations

We discuss how different properties at surfaces could require for different theoretical treatments within the first-principles density functional theory. Energies and structures are accurately determined by adopting the supercell geometry. Surface states are more conveniently described by the Green function embedding approach, which is able to take into account a truly semi-infinite solid and hence real continuous spectra. In this way a detailed analysis of discrete and resonant states is provided. We mainly describe the embedding method and provide examples to compare the two approaches. We focus next on the structural and electronic properties of alkali adatoms. The adsorption structure of Na/Cu(001) at low coverages is calculated within the supercell geometry motivated by the results of the Cambridge group on surface diffusion by 3He spin echo scattering. The dispersion, energy, effective mass, and width of surface (quantum well and image) states of alkali atoms on Cu(111) are worked out by the embedding approach and compared with experiments.

Femtomagnetism in graphene induced by core level excitation of organic adsorbates

We predict the induction or suppression of magnetism in the valence shell of physisorbed and chemisorbed organic molecules on graphene occurring on the femtosecond time scale as a result of core level excitations. For physisorbed molecules, where the interaction with graphene is dominated by van der Waals forces and the system is non-magnetic in the ground state, numerical simulations based on density functional theory show that the valence electrons relax towards a spin polarized configuration upon excitation of a core-level electron. The magnetism depends on efficient electron transfer from graphene on the femtosecond time scale. On the other hand, when graphene is covalently functionalized, the system is magnetic in the ground state showing two spin dependent mid gap states localized around the adsorption site. At variance with the physisorbed case upon core-level excitation, the LUMO of the molecule and the mid gap states of graphene hybridize and the relaxed valence shell is not magnetic anymore.

Image resonance in the many-body density of states at a metal surface

The electronic properties of a semi-infinite metal surface without a bulk gap are studied by a formalism that is able to account for the continuous spectrum of the system. The density of states at the surface is calculated within the GW approximation of many-body perturbation theory. We demonstrate the presence of an unoccupied surface resonance peaked at the position of the first image state. The resonance encompasses the whole Rydberg series of image states and cannot be resolved into individual peaks. Its origin is the shift in spectral weight when many-body correlation effects are taken into account.

Unusually Large Magnetic Anisotropy in Electrochemically Deposited Co-Rich Co–Pt Films

Co-rich Co-Pt films grown by electrodeposition from an amino-nitrite/citrate/glycine electrolyte onto Au(111) substrates apparently grow with a hexagonal structure, with its c-axis directed perpendicular to the surface. The films exhibit a perpendicular magnetic anisotropy (MCA) of the same order of magnitude as the shape anisotropy. Experimental estimates of the MCA result in a higher anisotropy than that reported for bulk materials of the same composition, but similar to values measured in films grown by vacuum methods at relatively high temperature, which partly consist of a high anisotropy, metastable orthorhombic Pmm2 phase. Comparison of valence band X-ray photoelectron spectroscopy measurements on electrodeposited films with density functional theory simulations of the electronic structure of the various reported Co(3)Pt structures support the notion that the films may consist of a mixture of the hexagonal and the Pmm2 structure.