Lara Ferrighi

Electronic structure calculations with GPAW: A real-space implementation of the projector augmented-wave method

Electronic structure calculations have become an indispensable tool in many areas of materials science and quantum chemistry. Even though the Kohn-Sham formulation of the density-functional theory (DFT) simplifies the many-body problem significantly, one is still confronted with several numerical challenges. In this article we present the projector augmented-wave (PAW) method as implemented in the GPAW program package (https://wiki.fysik.dtu.dk/gpaw) using a uniform real-space grid representation of the electronic wavefunctions. Compared to more traditional plane wave or localized basis set approaches, real-space grids offer several advantages, most notably good computational scalability and systematic convergence properties. However, as a unique feature GPAW also facilitates a localized atomic-orbital basis set in addition to the grid. The efficient atomic basis set is complementary to the more accurate grid, and the possibility to seamlessly switch between the two representations provides great flexibility. While DFT allows one to study ground state properties, time-dependent density-functional theory (TDDFT) provides access to the excited states. We have implemented the two common formulations of TDDFT, namely the linear-response and the time propagation schemes. Electron transport calculations under finite-bias conditions can be performed with GPAW using non-equilibrium Green functions and the localized basis set. In addition to the basic features of the real-space PAW method, we also describe the implementation of selected exchange-correlation functionals, parallelization schemes, ΔSCF-method, x-ray absorption spectra, and maximally localized Wannier orbitals.

2D−3D Transition for Cationic and Anionic Gold Clusters: A Kinetic Energy Density Functional Study

We present a density functional theory study of the energetics of isolated Au(n)+ (n = 5-10) and Au(n)- (n = 8-13) gold clusters. We compare our results to both theoretical and experimental values from the literature and find the use of meta-generalized gradient approximation (MGGA) functionals, in particular the M06-L functional, to be of importance in order to match experiment. The M06-L values suggest crossovers between 2D and 3D structures at n = 8 and 12 for cationic and anionic clusters, respectively. We suggest that the MGGAs stronger tendency toward 3D structures arises from their smaller gradient enhancement. Moreover, we show how MGGAs, in contrast to generalize gradient approximations with smaller gradient enhancements, avoid overestimating the bond energies by combining the information contained in the reduced gradient and the kinetic energy. This allows MGGAs to treat differently the exchange enhancement in the decaying density and bonding regions.

Density functional theory based screening of ternary alkali-transition metal borohydrides: A computational material design project

We present a computational screening study of ternary metal borohydrides for reversible hydrogen storage based on density functional theory. We investigate the stability and decomposition of alloys containing 1 alkali metal atom, Li, Na, or K (M(1)); and 1 alkali, alkaline earth or 3d/4d transition metal atom (M(2)) plus two to five (BH(4))(-) groups, i.e., M(1)M(2)(BH(4))(2-5), using a number of model structures with trigonal, tetrahedral, octahedral, and free coordination of the metal borohydride complexes. Of the over 700 investigated structures, about 20 were predicted to form potentially stable alloys with promising decomposition energies. The M(1)(Al/Mn/Fe)(BH(4))(4), (Li/Na)Zn(BH(4))(3), and (Na/K)(Ni/Co)(BH(4))(3) alloys are found to be the most promising, followed by selected M(1)(Nb/Rh)(BH(4))(4) alloys.

Direct observation of molecular preorganization for chirality transfer on a catalyst surface.

Scanning tunneling microscopy and theoretical calculations shed light on an asymmetric heterogeneous catalyst. The chemisorption of specific optically active compounds on metal surfaces can create catalytically active chirality transfer sites. However, the mechanism through which these sites bias the stereoselectivity of reactions (typically hydrogenations) is generally assumed to be so complex that continued progress in the area is uncertain. We show that the investigation of heterogeneous asymmetric induction with single-site resolution sufficient to distinguish stereochemical conformations at the submolecular level is finally accessible. A combination of scanning tunneling microscopy and density functional theory calculations reveals the stereodirecting forces governing preorganization into precise chiral modifier-substrate bimolecular surface complexes. The study shows that the chiral modifier induces prochiral switching on the surface and that different prochiral ratios prevail at different submolecular binding sites on the modifier at the reaction temperature.

Stereodirection of an α-Ketoester at Sub-molecular Sites on Chirally Modified Pt(111): Heterogeneous Asymmetric Catalysis

Chirally modified Pt catalysts are used in the heterogeneous asymmetric hydrogenation of α-ketoesters. Stereoinduction is believed to occur through the formation of chemisorbed modifier-substrate complexes. In this study, the formation of diastereomeric complexes by coadsorbed methyl 3,3,3-trifluoropyruvate, MTFP, and (R)-(+)-1-(1-naphthyl)ethylamine, (R)-NEA, on Pt(111) was studied using scanning tunneling microscopy and density functional theory methods. Individual complexes were imaged with sub-molecular resolution at 260 K and at room temperature. The calculations find that the most stable complex isolated in room-temperature experiments is formed by the minority rotamer of (R)-NEA and pro-S MTFP. The stereodirecting forces in this complex are identified as a combination of site-specific chemisorption of MTFP and multiple non-covalent attractive interactions between the carbonyl groups of MTFP and the amine and aromatic groups of (R)-NEA.

Self-consistent meta-generalized gradient approximation study of adsorption of aromatic molecules on noble metal surfaces

The adsorption of benzene, pyridine, and two nucleobases on the Au(111) surface has been investigated using a fully relaxed, self-consistent meta-generalized gradient approximation (meta-GGA) density functional theory setup with the M06-L functional. The meta-GGA based molecule-surface separations are shortened and the adsorption bond strengths of the molecules are greatly improved over the virtually non-interacting results obtained when using a plain GGA exchange-correlation functional. The nucleobases containing oxygen atoms show higher corrugation with adsorption site and orientation than the other aromatic molecules considered. The adsorption of pentacene is studied on Au, Ag, and Cu surfaces. In agreement with experiment, the adsorption energies are found to increase with decreasing nobleness, but the dependency is underestimated. We point out how the kinetic energy density can discriminate between covalent and non-covalent bonding regions of orbital overlap.

On-surface photo-dissociation of C–Br bonds: towards room temperature Ullmann coupling

The surface-assisted synthesis of gold-organometallic hybrids on the Au(111) surface both by thermo- and light-initiated dehalogenation of bromo-substituted tetracene is reported. Combined X-ray photoemission (XPS) and scanning tunneling microscopy (STM) data reveal a significant increase of the surface order when mild reaction conditions are combined with 405 nm light irradiation.

Control of the Intermolecular Coupling of Dibromotetracene on Cu(110) by the Sequential Activation of CBr and CH Bonds

Dibromotetracene molecules are deposited on the Cu(110) surface at room temperature. The complex evolution of this system has been monitored at different temperatures (i.e., 298, 523, 673, and 723 K) by means of a variety of complementary techniques that range from STM and temperature-programmed desorption (TPD) to high-resolution X-ray spectroscopy (XPS) and near-edge X-ray absorption fine structure spectroscopy (NEXAFS). State-of-the-art density-functional calculations were used to determine the chemical processes that take place on the surface. After deposition at room temperature, the organic molecules are transformed into organometallic monomers through debromination and carbon-radical binding to copper adatoms. Organometallic dimers, trimers, or small oligomers, which present copper-bridged molecules, are formed by increasing the temperature. Surprisingly, further heating to 673 K causes the formation of elongated chains along the Cu(110) close-packed rows as a consequence of radical-site migration to the thermodynamically more stable molecule heads. Finally, massive dehydrogenation occurs at the highest temperature followed by ring condensation to nanographenic patches. This study is a paradigmatic example of how intermolecular coupling can be modulated by the stepwise control of a simple parameter, such as temperature, through a sequence of domino reactions.

Theoretical Studies of Oxygen Reactivity of Free‐Standing and Supported Boron‐Doped Graphene

Graphene inertness towards chemical reactivity can be considered as an accepted postulate by the research community. This limit has been recently overcome by chemically and physically modifying graphene through non-metal doping or interfacing with acceptor/donor materials (metals or semiconductors). As a result, outstanding performances as catalytic, electrocatalytic, and photocatalytic material have been observed. In this critical Review we report computational work performed, by our group, on the reactivity of free-standing, metal- and semiconductor-supported B-doped graphene towards oxygen, which is at the basis of extremely important energy-related chemical processes, such as the oxygen reduction reaction. It appears that a combination of doping and interfacing approaches for the activation of graphene can open unconventional and unprecedented reaction paths, thus boosting the potential of modified graphene in many chemical applications.

Synthesis of graphene nanoribbons with a defined mixed edge-site sequence by surface assisted polymerization of (1,6)-dibromopyrene on Ag(110)

By a combination of scanning tunneling microscopy, X-ray spectroscopic techniques and density functional theory calculations, we prove the formation of extended patterns of parallel, graphene nanoribbons with alternate zig-zag and armchair edges and selected width by surface-assisted Ullmann coupling polymerization and dehydrogenation of 1,6-dibromopyrene (C16H8Br2). Besides the relevance of these nanostructures for their possible application in nanodevices, we demonstrate the peculiarity of halogenated pyrene derivatives for the formation of nanoribbons, in particular on Ag(110). These results open the possibility of tuning the shape and dimension of nanoribbons (and hence the correlated electronic properties) by choosing suitably tailored or on-purpose designed molecular precursors.