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

Light-Triggered Reversible Supracolloidal Self-Assembly of Precision Gold Nanoclusters

Monolayer thiol-protected noble metal nanoclusters are attractive nanoscale building blocks for well-defined colloidal superstructures. However, achieving facile reversible self-assembly of nanoclusters using external stimuli is still in its infancy. Herein, we report the synthesis and photon-assisted reversible self-assembly of thiolated azobenzene-stapled Au₂₅ nanoclusters. Photoactivation of functionalized nanoclusters in dichloromethane by irradiating ultraviolet light at 345 nm results in a visual change and formation of disc-like colloidal superstructures (d ∼ 100-1000 nm). The superstructures readily disassemble into individual nanoclusters upon irradiating with visible light at 435 nm. Systematic changes in both the electronic absorption bands and nuclear magnetic resonance spectra of chromophores in solution suggest that the photoisomerization of surface ligands drives the self-assembly. High-resolution transmission electron microscopy, electron tomographic reconstruction, dynamic light scattering, and small-angle X-ray powder diffraction show that the disc-like superstructures contain densely packed nanoclusters. Long-range self-assembly and disassembly under ultraviolet and visible light, respectively, demonstrate reversible photoswitching in nanoclusters.

The role of an interface in stabilizing reaction intermediates for hydrogen evolution in aprotic electrolytes

By combining idealized experiments with realistic quantum mechanical simulations of an interface, we investigate electro-reduction reactions of HF, water and methanesulfonic acid (MSA) on the single crystal (111) facets of Au, Pt, Ir and Cu in organic aprotic electrolytes, 1 M LiPF₆ in EC/EMC 3:7W (LP57), the aprotic electrolyte commonly used in Li-ion batteries, 1 M LiClO₄ in EC/EMC 3:7W and 0.2 M TBAPF₆ in 3 : 7 EC/EMC. In our previous work, we have established that LiF formation, accompanied by H₂ evolution, is caused by a reduction of HF impurities and requires the presence of Li at the interface, which catalyzes the HF dissociation. In the present paper, we find that the measured potential of the electrochemical response for these reduction reactions correlates with the work function of the electrode surfaces and that the work function determines the potential for Li⁺ adsorption. The reaction path is investigated further by electrochemical simulations suggesting that the overpotential of the reaction is related to stabilizing the active structure of the interface having adsorbed Li⁺. Li⁺ is needed to facilitate the dissociation of HF which is the source of protons. Further experiments on other proton sources, water and methanesulfonic acid, show that if the hydrogen evolution involves negatively charged intermediates, F^- or HO^-, a cation at the interface can stabilize them and facilitate the reaction kinetics. When the proton source is already significantly dissociated (in the case of a strong acid), there is no negatively charged intermediate and thus the hydrogen evolution can proceed at much lower overpotentials. This reveals a situation where the overpotential for electrocatalysis is related to stabilizing the active structure of the interface, facilitating the reaction rather than providing the reaction energy.

Steric Effects Control Dry Friction of H- and F-Terminated Carbon Surfaces

A stable passivation of surface dangling bonds underlies the outstanding friction properties of diamond and diamond-like carbon (DLC) coatings in boundary lubrication. While hydrogen is the simplest termination of a carbon dangling bond, fluorine can also be used as a monoatomic termination, providing an even higher chemical stability. However, whether and under which conditions a substitution of hydrogen with fluorine can be beneficial to friction is still an open question. Moreover, which of the chemical differences between C-H and C-F bonds are responsible for the change in friction has not been unequivocally understood yet. In order to shed light on this problem, we develop a density functional theory-based, nonreactive force field that describes the relevant properties of hydrogen- and fluorine-terminated diamond and DLC tribological interfaces. Molecular dynamics and nudged elastic band simulations reveal that the frictional stress at such interfaces correlates with the corrugation of the contact potential energy, thus ruling out a significant role of the mass of the terminating species on friction. Furthermore, the corrugation of the contact potential energy is almost exclusively determined by steric factors, while electrostatic interactions only play a minor role. In particular, friction between atomically flat diamond surfaces is controlled by the density of terminations, by the C-H and C-F bond lengths, and by the H and F atomic radii. For sliding DLC/DLC interfaces, the intrinsic atomic-scale surface roughness plays an additional role. While surface fluorination decreases the friction of incommensurate diamond contacts, it can negatively affect the friction performance of carbon surfaces that are disordered and not atomically flat. This work provides a general framework to understand the impact of chemical structure of surfaces on friction and to generate design rules for optimally terminated low-friction systems.

Oxygen evolution reaction: a perspective on a decade of atomic scale simulations

Multiple strategies to overcome the intrinsic limitations of the oxygen evolution reaction (OER) have been proposed by numerous research groups. Despite the substantial efforts, the driving force required for water oxidation is largely making the reaction inefficient. In the present work, we collected published studies involving DFT calculations for the OER, with the purpose to understand why the progress made so far, for lowering the overpotential of the reaction, is relatively small. The data revealed that the universal scaling relationship between HO* and HOO* intermediates is still present and robust, despite the variety in methods and structures used for calculating the binding energies of the intermediates. On the other hand, the data did not show a clear trend line regarding the O* binding. Our analysis suggested that trends in doped semiconducting oxides behave very differently from those in other oxides. This points towards a computational challenge in describing doped oxides in a realistic manner. We propose a way to overcome these computational challenges, which can be applied to simulations corresponding to doped semiconductors in general.

Understanding the Electron-Doping Mechanism in Potassium-Intercalated Single-Walled Carbon Nanotubes

Single-walled carbon nanotubes (SWCNTs) can be doped with potassium, similar to graphite, leading to intercalation compounds. These binary systems exhibit a clear metallic character. However, the entire picture of how electron doping (e-doping) modifies the SWCNTs' vibrational spectra as a function of their diameter, chirality, and metallicity is still elusive. Herein, we present a detailed study of the intercalation and solid state reduction of metallic and semiconducting enriched HiPco SWCNTs. We performed a combined experimental and theoretical study of the evolution of their Raman response with potassium exposure, focusing specifically on their radial breathing mode (RBM). We found the charge donated from the potassium atoms occupies antibonding π orbitals of the SWCNTs, weakening their C-C bonds, and reducing the RBM frequency. This RBM downshift with increasing doping level is quasi-linear with a steplike behavior when the Fermi level crosses a van Hove singularity for semiconducting species. Moreover, this weakening of the C-C bonds is greater with decreasing curvature, or increasing diameter. Overall, this suggests the RBM downshift with e-doping is proportional to both the SWCNT's integrated density of states (DOS) ϱ(ε) and diameter d. We have provided a precise and complete description of the complex electron doping mechanism in SWCNTs up to a charge density of -18 me/C, far beyond that achievable by standard gate voltage studies, not being the highest doping possible, but high enough to track the effects of doping in SWCNTs based on their excitation energy, diameter, band gap energy, chiral angle, and metallicity. This work is highly relevant to tuning the electronic properties of SWCNTs for applications in nanoelectronics, plasmonics, and thermoelectricity.

QuantumATK: an integrated platform of electronic and atomic-scale modelling tools

QuantumATK is an integrated set of atomic-scale modelling tools developed since 2003 by professional software engineers in collaboration with academic researchers. While different aspects and individual modules of the platform have been previously presented, the purpose of this paper is to give a general overview of the platform. The QuantumATK simulation engines enable electronic-structure calculations using density functional theory or tight-binding model Hamiltonians, and also offers bonded or reactive empirical force fields in many different parametrizations. Density functional theory is implemented using either a plane-wave basis or expansion of electronic states in a linear combination of atomic orbitals. The platform includes a long list of advanced modules, including Green's-function methods for electron transport simulations and surface calculations, first-principles electron-phonon and electron-photon couplings, simulation of atomic-scale heat transport, ion dynamics, spintronics, optical properties of materials, static polarization, and more. Seamless integration of the different simulation engines into a common platform allows for easy combination of different simulation methods into complex workflows. Besides giving a general overview and presenting a number of implementation details not previously published, we also present four different application examples. These are calculations of the phonon-limited mobility of Cu, Ag and Au, electron transport in a gated 2D device, multi-model simulation of lithium ion drift through a battery cathode in an external electric field, and electronic-structure calculations of the composition-dependent band gap of SiGe alloys.

Large-Z limit in atoms and solids from first principles

We present density functional theory (DFT) calculations of atomic ionization potentials and lattice constants of simple solids from low atomic numbers Z to the large-Z limit. We compare different kinetic energy functional approximations [Kohn-Sham (KS) vs simple orbital-free functionals] and, in the case of orbital free, also different methods for including the nuclear potential (all-electron with the projector augmented wave method vs local pseudopotentials). For both ionization potentials and lattice constants, all-electron orbital-free DFT does yield the general trend of KS DFT for moderate values of the atomic number Z. For large values of Z, all-electron orbital-free DFT deviates from the KS DFT results. Local pseudopotentials give a better qualitative description by adding shell oscillations to the orbital-free DFT model. We show that both all-electron orbital-free DFT and KS DFT have a finite value for nonrelativistic lattice constants in the large-Z limit.

Single-Shot Multi-keV X-Ray Absorption Spectroscopy Using an Ultrashort Laser-Wakefield Accelerator Source

Single-shot absorption measurements have been performed using the multi-keV x rays generated by a laser-wakefield accelerator. A 200 TW laser was used to drive a laser-wakefield accelerator in a mode which produced broadband electron beams with a maximum energy above 1 GeV and a broad divergence of ≈15  mrad FWHM. Betatron oscillations of these electrons generated 1.2±0.2×10^{6}  photons/eV in the 5 keV region, with a signal-to-noise ratio of approximately 300∶1. This was sufficient to allow high-resolution x-ray absorption near-edge structure measurements at the K edge of a titanium sample in a single shot. We demonstrate that this source is capable of single-shot, simultaneous measurements of both the electron and ion distributions in matter heated to eV temperatures by comparison with density functional theory simulations. The unique combination of a high-flux, large bandwidth, few femtosecond duration x-ray pulse synchronized to a high-power laser will enable key advances in the study of ultrafast energetic processes such as electron-ion equilibration.

Experimental and DFT Studies of Au Deposition Over WO3/g-C3N4 Z-Scheme Heterojunction

A typical Z-scheme system is composed of two photocatalysts which generate two sets of charge carriers and split water into H₂ and O₂ at different locations. Scientists are struggling to enhance the efficiencies of these systems by maximizing their light absorption, engineering more stable redox couples, and discovering new O₂ and H₂ evolutions co-catalysts. In this work, Au decorated WO₃/g-C₃N₄ Z-scheme nanocomposites are fabricated via wet-chemical and photo-deposition methods. The nanocomposites are utilized in photocatalysis for H₂ production and 2,4-dichlorophenol (2,4-DCP) degradation. It is investigated that the optimized 4Au/6% WO₃/CN nanocomposite is highly efficient for production of 69.9 and 307.3 µmol h^-1 g^-1 H₂ gas, respectively, under visible-light (λ > 420 nm) and UV-visible illumination. Further, the fabricated 4Au/6% WO₃/CN nanocomposite is significant (i.e., 100% degradation in 2 h) for 2,4-DCP degradation under visible light and highly stable in photocatalysis. A significant 4.17% quantum efficiency is recorded for H₂ production at wavelength 420 nm. This enhanced performance is attributed to the improved charge separation and the surface plasmon resonance effect of Au nanoparticles. Solid-state density functional theory simulations are performed to countercheck and validate our experimental data. Positive surface formation energy, high charge transfer, and strong non-bonding interaction via electrostatic forces confirm the stability of 4Au/6% WO₃/CN interface.

Polarizable Embedding with a Transferable H2O Potential Function II: Application to (H2O)n Clusters and Liquid Water

The incorporation of polarization in multiscale quantum-mechanics/molecular-mechanics (QM/MM) simulations is important for a variety of applications, for example, charge-transfer reactions. A recently developed formalism based on a density functional theory description of the QM region and a potential energy function for H₂O molecules that includes quadrupole as well as dipole polarizability of the MM region is used to simulate liquid water and water clusters. Analysis of the energy, atomic forces, MM polarization, and structure is presented. A quantitative assessment of the QM/MM-MM/MM interaction energy differences of all possible QM/MM configurations of (H₂O)_n clusters shows that the interquartile range of the distributions of the QM/MM binding energies is never more than 20 meV/molecule higher or lower than the binding energies produced with either of the single-model results. Comparing these interaction energy differences with the QM/MM induction differences show that they are not systematically caused by the induced MM moments of our polarizable embedding scheme. Optimized hexamer geometries as well as the liquid water structure are shown to be improved in comparison with results obtained using point-charge based embedding models neglecting polarization.

Polarizable Embedding with a Transferable H2O Potential Function I: Formulation and Tests on Dimer

The incorporation of mutual polarization in multiscale simulations where different regions of the system are treated at different level of theory is important in studies of, for example, electronic excitations and charge transfer processes. We present here an energy functional for describing a quantum mechanics/molecular mechanics (QM/MM) scheme that includes reciprocal polarization between the two subsystems. The inclusion of polarization alleviates shortcomings inherent in electrostatic embedding QM/MM models based on point-charge force fields. A density functional theory (DFT) description of the QM subsystem is coupled to a single center multipole expansion (SCME) description of H₂O molecules in the MM subsystem that includes anisotropic dipole and quadrupole polarizability as well as static multipoles up to and including the hexadecapole. The energy functional and the coupling scheme is general and can be extended to arbitrary order in terms of both the static and induced moments. Tests of the energy surface for the H₂O dimer show that the QM/MM results lie in between the pure DFT and pure SCME values. The consistency of the many-body contributions to the energy and analytical forces is demonstrated for an H₂O pentamer.

Classifying the Electronic and Optical Properties of Janus Monolayers

Inspired by the recent synthesis of monolayer MoSSe, we conduct a first-principles high-throughput investigation of 216 MXY Janus monolayers consisting of a middle layer of metal atoms (M) sandwiched between different types of chalcogen, halogen, or pnictogen atoms (X,Y). Using density functional theory and many-body perturbation theory, we perform an exhaustive computational characterization of the 70 most stable semiconducting monolayers. These are found to exhibit diverse and fascinating properties including finite out-of-plane dipoles, giant Rashba-splittings, direct and indirect band gaps ranging from 0.7 to 3.0 eV, large exciton binding energies, and very strong light-matter interactions. The data have been generated using the workflow behind the Computational 2D Materials Database and are freely available online. Our work expands the class of known Janus monolayers and points to several potentially synthesizable structures, which could be interesting candidates for valley- or optoelectronic applications or for generating out-of-plane electric fields to control charge transfer, charge separation, or band alignments in van der Waals heterostructures.

Descriptors for Electrolyte-Renormalized Oxidative Stability of Solvents in Lithium-Ion Batteries

Electrolyte stability against oxidation is one of the important factors limiting the development of high energy density batteries. The highest occupied molecular orbital (HOMO) level of solvent molecules has been traditionally used for understanding trends in their oxidative stability, but this assumes a noninteracting environment. However, solvent HOMO levels are renormalized because of molecules in their solvation shells. In this work, we first demonstrate an inexpensive and accurate method to determine the HOMO level of the solvent followed by simple descriptors for renormalization of the HOMO level due to different electrolyte components. The descriptors are based on Gutmann donor and acceptor numbers of the solvent and other components. The method uses fast generalized gradient approximation-level density functional theory calculations compared to previously used expensive, experimental data-dependent methods. This method can be used to screen for unexplored stable solvents among the large number of known organic compounds to design novel high-voltage stable electrolytes.

Nonstoichiometric Phases of Two-Dimensional Transition-Metal Dichalcogenides: From Chalcogen Vacancies to Pure Metal Membranes

Two-dimensional (2D) membranes consisting of a single layer of Mo atoms were recently manufactured [ Adv. Mater. 2018 , 30 , 1707281 ] from MoSe₂ sheets by sputtering Se atoms using an electron beam in a transmission electron microscope. This is an unexpected result as formation of Mo clusters should energetically be more favorable. To get microscopic insights into the energetics of realistic Mo membranes and nonstoichiometric phases of transition-metal dichalcogenides (TMDs) M_aX_b, where M = Mo and W and X = S, Se, and Te, we carry out first-principles calculations and demonstrate that the membranes, which can be referred to as metallic quantum dots embedded into a semiconducting matrix, can be stabilized by charge transfer. We also show that an ideal neutral 2D Mo or W sheet is not flat but a corrugated structure, with a square lattice being the lowest-energy configuration. We further demonstrate that several intermediate nonstoichiometric phases of TMDs are possible as they have lower formation energies than pure metal membranes. Among them, the orthorhombic metallic 2D M₄X₄ phase is particularly stable. Finally, we study the properties of this phase in detail and discuss how it can be manufactured by the top-down approaches.

Optical excitations of chlorophyll a and b monomers and dimers

A necessary first step in the development of technologies such as artificial photosynthesis is understanding the photoexcitation process within the basic building blocks of naturally occurring light harvesting complexes (LHCs). The most important of these building blocks in biological LHCs such as LHC II from green plants are the chlorophyll a (Chl a) and chlorophyll b (Chl b) chromophores dispersed throughout the protein matrix. However, efforts to describe such systems are still hampered by the lack of computationally efficient and accurate methods that are able to describe optical absorption in large biomolecules. In this work, we employ a highly efficient linear combination of atomic orbitals (LCAOs) to represent the Kohn-Sham (KS) wave functions at the density functional theory (DFT) level and perform time-dependent density functional theory (TDDFT) calculations in either the reciprocal space and frequency domain (LCAO-TDDFT-k-ω) or real space and time domain (LCAO-TDDFT-r-t) of the optical absorption spectra of Chl a and b monomers and dimers. We find that our LCAO-TDDFT-k-ω and LCAO-TDDFT-r-t calculations reproduce results obtained with a plane-wave (PW) representation of the KS wave functions (PW-TDDFT-k-ω) but with a significant reduction in computational effort. Moreover, by applying the Gritsenko, van Leeuwen, van Lenthe, and Baerends solid and correlation derivative discontinuity correction Δ_x to the KS eigenenergies, with both LCAO-TDDFT-k-ω and LCAO-TDDFT-r-t methods, we are able to semiquantitatively reproduce the experimentally measured photoinduced dissociation results. This work opens the path to first principles calculations of optical excitations in macromolecular systems.

Synthesis and Doping Strategies to Improve the Photoelectrochemical Water Oxidation Activity of BiVO4 Photoanodes

BiVO4 is one of the most investigated and most promising metal oxide based photoanode materials for photoelectrochemical (PEC) water oxidation. Although it has several advantages (suitable band gap around 2.4 eV, suitable valence-band position for water oxidation, low toxicity, high abundance), it suffers from slow charge-carrier transport properties, high surface recombination, and limited water-oxidation activity. In the present work, we review the synthesis and doping strategies that we developed in the last years to improve the PEC performance of BiVO4 photoanodes. Strategies ranging from single anion doping or cation doping to anion and cation co-doping will be presented for fluoride and molybdenum as anion and cation dopants, respectively. One major result is that co-doping allows combining the most important PEC specific benefits of each type of dopant, i.e. an increased charge-injection efficiency in case of fluoride as well as an increased charge-separation efficiency in case of molybdenum.

Fast nonadiabatic dynamics of many-body quantum systems

Modeling many-body quantum systems with strong interactions is one of the core challenges of modern physics. A range of methods has been developed to approach this task, each with its own idiosyncrasies, approximations, and realm of applicability. However, there remain many problems that are intractable for existing methods. In particular, many approaches face a huge computational barrier when modeling large numbers of coupled electrons and ions at finite temperature. Here, we address this shortfall with a new approach to modeling many-body quantum systems. On the basis of the Bohmian trajectory formalism, our new method treats the full particle dynamics with a considerable increase in computational speed. As a result, we are able to perform large-scale simulations of coupled electron-ion systems without using the adiabatic Born-Oppenheimer approximation.

Calculating spin crossover temperatures by a first-principles LDA+U scheme with parameter U evaluated from GW

The prediction of spin crossover (SCO) temperatures (T_1/2) depends sensitively on the description of local Coulomb correlation. Due to its balance between accuracy and computational cost, local density approximation combined with Hubbard U model (LDA+U) is an appealing tool for this purpose. Despite its accurate performance on energetic properties, such as spin adiabatic energy difference, it is well-known that the LDA+U approach would lose its predictive power if U is tuned to achieve close agreement with experiment for a certain property. On the other hand, a static U value cannot account for changes in the electronic structure. Here, we propose a framework to derive dynamical U (U_dyn) values for iron(ii) complexes from the many-body GW calculations. By performing model calculations on a series of compounds with varying ligand fields, we show that the U values determined in this way are local environment dependent, and the resulting LDA+U_dyn method could reproduce their experimental ground spin states. We present applications to selected SCO complexes illustrating that U_dyn considerably overcomes some of the drawbacks of employing a constant U in the calculation of thermochemical quantities. Using the described calculation procedure, the T_1/2 values are predicted with a small mean absolute error of 176 K with respect to experiment.

Permethylation Introduces Destructive Quantum Interference in Saturated Silanes

The single-molecule conductance of silanes is suppressed due to destructive quantum interference in conformations with cisoid dihedral angles along the molecular backbone. Yet, despite the structural similarity, σ-interference effects have not been observed in alkanes. Here we report that the methyl substituents used in silanes are a prerequisite for σ-interference in these systems. Through density functional theory calculations, we find that the destructive interference is not evident to the same extent in nonmethylated silanes. We find the same is true in alkanes as the transmission is significantly suppressed in permethylated cyclic and bicyclic alkanes. Using scanning tunneling microscope break-junction method we determine the single-molecule conductance of functionalized cyclohexane and bicyclo[2.2.2]octane that are found to be higher than that of equivalent permethylated silanes. Rather than the difference between carbon and silicon atoms in the molecular backbones, our calculations reveal that it is primarily the difference between hydrogen and methyl substituents that result in the different electron transport properties of nonmethylated alkanes and permethylated silanes. Chemical substituents play an important role in determining the single-molecule conductance of saturated molecules, and this must be considered when we improve and expand the chemical design of insulating organic molecules.

Hierarchy in the Halogen Activation During Surface-Promoted Ullmann Coupling

Within the collection of surface-supported reactions currently accessible for the production of extended molecular nanostructures under ultra-high vacuum, Ullmann coupling has been the most successful in the controlled formation of covalent single C-C bonds. Particularly advanced control of this synthetic tool has been obtained by means of hierarchical reactivity, commonly achieved by the use of different halogen atoms that consequently display distinct activation temperatures. Here we report on the site-selective reactivity of certain carbon-halogen bonds. We use precursor molecules halogenated with bromine atoms at two non-equivalent carbon atoms and found that the Ullmann coupling occurs on Au(111) with a remarkable predilection for one of the positions. Experimental evidence is provided by means of scanning tunneling microscopy and core level photoemission spectroscopy, and a rationalized understanding of the observed preference is obtained from density functional theory calculations.

A method for structure prediction of metal-ligand interfaces of hybrid nanoparticles

Hybrid metal nanoparticles, consisting of a nano-crystalline metal core and a protecting shell of organic ligand molecules, have applications in diverse areas such as biolabeling, catalysis, nanomedicine, and solar energy. Despite a rapidly growing database of experimentally determined atom-precise nanoparticle structures and their properties, there has been no successful, systematic way to predict the atomistic structure of the metal-ligand interface. Here, we devise and validate a general method to predict the structure of the metal-ligand interface of ligand-stabilized gold and silver nanoparticles, based on information about local chemical environments of atoms in experimental data. In addition to predicting realistic interface structures, our method is useful for investigations on the steric effects at the metal-ligand interface, as well as for predicting isomers and intermediate structures induced by thermal dynamics or interactions with the environment. Our method is applicable to other hybrid nanomaterials once a suitable set of reference structures is available.

Atomistic structure prediction of the metal-ligand interface of hybrid nanoparticles remains challenging. Here the authors present an algorithm to predict the structure of the metal-ligand interface of ligand-stabilized gold and silver nanoparticles, guided by experimental data on local chemical environments.

Structural, vibrational and electronic properties of SnMBO4 (M = Al, Ga): a predictive hybrid DFT study

We propose two new members of the mullite-type family, SnAlBO₄ and SnGaBO₄, and carry out an in-depth study of their crystal properties using the hybrid method PW1PW. Both are isostructural to PbMBO₄ (M  =  Fe, Mn, Al, Ga), which show axial negative linear compressibility (ANLC), among other interesting features. We find that, although Sn²⁺ is susceptible of being oxidized by oxygen, a suitable range of experimental parameters exists in which the compounds could be synthesized. We observe absence of ANLC below 20 GPa and explain it by the small space occupied by the lone electron pairs, as indicated by the small length of the corresponding Liebau Density Vectors. In agreement with this fact, the structures present a low number of negative mode-Grüneisen parameters, which may also suggest lack of negative thermal expansion. The electronic properties show a remarkable anisotropic behaviour, with a strong dependence of the absorption spectra on light polarization direction.

van der Waals exchange-correlation functionals over bulk and surface properties of transition metals

The correct description of van der Waals (vdW) interaction forces is required for accurately describing dispersion bonded systems. Several approaches have been proposed to include London dispersion in density functional theory exchange-correlation functionals, where the family of so-called van der Waals (vdW-DF) exchange-correlation functionals have shown a better performance than local or semi local exchange-correlation functionals for describing molecular adsorption on metals. Despite the numerous benchmarks performed with these functionals, their performance in predicting bulk properties of transition metals has hitherto not been investigated in detail. We have therefore tested five vdW-DF exchange-correlation functionals, vdW-DF2, optPBE-vdW, BEEF-vdW, optB88-vdW and C09 _x -vdW to assess their performance in the prediction of lattice constants, bulk moduli, cohesive energies and surface energies of bulk Ni, Cu, Rh, Pd, Ag, Ir, Pt and Au (in fcc crystal structure). These transition metals are commonly used for benchmarking density functionals because they are important for applications in catalysis. The results are compared with experimental data and the PBE exchange-correlation functional. We found that both the optB88-vdW and the C09 _x -vdW exchange-correlation functionals estimate all properties with high accuracy, in better agreement with experimental data than PBE and other considered vdW functionals. The C09 _x -vdW functional clearly outperforms all other exchange-correlation functionals for surface energies for the (1 1 1) termination of different metals. We have also evaluated the interatomic electron density emerging from different functionals, and concluded that the observed differences are a result of the predicted lattice parameter, rather than a direct consequence of the functional form. Plane-wave and real-space grid-based expansions of the electron density are also compared, revealing good agreement between the two approaches for lattice parameters, cohesive energies, and surface energies, but more severe differences in bulk moduli. On the basis of our results, we recommend using the C09 _x -vdW for studying bulk properties and surface energies of transition metals.

Crystal field coefficients for yttrium analogues of rare-earth/transition-metal magnets using density-functional theory in the projector-augmented wave formalism

We present a method of calculating crystal field coefficients of rare-earth/transition-metal (RE-TM) magnets within density-functional theory (DFT). The principal idea of the method is to calculate the crystal field potential of the yttrium analogue ('Y-analogue') of the RE-TM magnet, i.e. the material where the lanthanide elements have been substituted with yttrium. The advantage of dealing with Y-analogues is that the methodological and conceptual difficulties associated with treating the highly-localized 4f  electrons in DFT are avoided, whilst the nominal valence electronic structure principally responsible for the crystal field is preserved. In order to correctly describe the crystal field potential in the core region of the atoms we use the projector-augmented wave formalism of DFT, which allows the reconstruction of the full charge density and electrostatic potential. The Y-analogue crystal field potentials are combined with radial 4f  charge densities obtained in self-interaction-corrected calculations on the lanthanides to obtain crystal field coefficients. We demonstrate our method on a test set of ten materials comprising nine RE-TM magnets and elemental Tb. We show that the calculated easy directions of magnetization agree with experimental observations, including a correct description of the anisotropy within the basal plane of Tb and NdCo₅. We further show that the Y-analogue calculations generally agree quantitatively with previous calculations using the open-core approximation to treat the 4f  electrons, and argue that our simple approach may be useful for large-scale computational screening of new magnetic materials.

Strong plasmon-molecule coupling at the nanoscale revealed by first-principles modeling

Strong light-matter interactions in both the single-emitter and collective strong coupling regimes attract significant attention due to emerging applications in quantum and nonlinear optics as well as opportunities for modifying material-related properties. Exploration of these phenomena is theoretically demanding, as polaritons exist at the intersection between quantum optics, solid state physics, and quantum chemistry. Fortunately, nanoscale polaritons can be realized in small plasmon-molecule systems, enabling treatment with ab initio methods. Here, we show that time-dependent density-functional theory calculations access the physics of nanoscale plasmon-molecule hybrids and predict vacuum Rabi splitting. By considering a system comprising a few-hundred-atom aluminum nanoparticle interacting with benzene molecules, we show that cavity quantum electrodynamics holds down to resonators of a few cubic nanometers in size, yielding a single-molecule coupling strength exceeding 200 meV due to a massive vacuum field of 4.5 V · nm⁻¹. In a broader perspective, ab initio methods enable parameter-free in-depth studies of polaritonic systems for emerging applications.

Light-matter interaction is described by many different approaches and approximations depending on the coupling strength. Here the authors model a plasmonic system of aluminum nanoparticle interacting with benzene molecules using TDDFT and show its validity to a strong plasmon-molecule coupling regime.

The GW Compendium: A Practical Guide to Theoretical Photoemission Spectroscopy

The GW approximation in electronic structure theory has become a widespread tool for predicting electronic excitations in chemical compounds and materials. In the realm of theoretical spectroscopy, the GW method provides access to charged excitations as measured in direct or inverse photoemission spectroscopy. The number of GW calculations in the past two decades has exploded with increased computing power and modern codes. The success of GW can be attributed to many factors: favorable scaling with respect to system size, a formal interpretation for charged excitation energies, the importance of dynamical screening in real systems, and its practical combination with other theories. In this review, we provide an overview of these formal and practical considerations. We expand, in detail, on the choices presented to the scientist performing GW calculations for the first time. We also give an introduction to the many-body theory behind GW, a review of modern applications like molecules and surfaces, and a perspective on methods which go beyond conventional GW calculations. This review addresses chemists, physicists and material scientists with an interest in theoretical spectroscopy. It is intended for newcomers to GW calculations but can also serve as an alternative perspective for experts and an up-to-date source of computational techniques.

Optical signatures of pentacene in soft rare-gas environments

Acenes and pentacene (Pc), in particular, are promising candidates for organic dyes with interesting properties important for solar light to energy conversion. We present a combined experimental and computational study of Pc in an ultracold environment that allows for high resolution optical spectroscopy. The spectra and their vibrational substructure are interpreted with the help of density functional theory calculations. While there are only slight changes within superfluid helium as compared to vacuum, the neon surface shows more prominent effects. Additional vibrational coupling by neon modes leads to broadening as well as the emergence of new features, like the otherwise symmetry forbidden out-of-plane butterfly mode.

Cavity Control of Excitons in Two-Dimensional Materials

We propose a robust and efficient way of controlling the optical spectra of two-dimensional materials and van der Waals heterostructures by quantum cavity embedding. The cavity light-matter coupling leads to the formation of exciton-polaritons, a superposition of photons and excitons. Our first-principles study demonstrates a reordering and mixing of bright and dark excitons spectral features and in the case of a type II van-der-Waals heterostructure an inversion of intra- and interlayer excitonic resonances. We further show that the cavity light-matter coupling strongly depends on the dielectric environment and can be controlled by encapsulating the active two-dimensional (2D) crystal in another dielectric material. Our theoretical calculations are based on a newly developed nonperturbative many-body framework to solve the coupled electron-photon Schrödinger equation in a quantum-electrodynamical extension of the Bethe-Salpeter approach. This approach enables the ab initio simulations of exciton-polariton states and their dispersion from weak to strong cavity light-matter coupling regimes. Our method is then extended to treat van der Waals heterostructures and encapsulated 2D materials using a simplified Mott-Wannier description of the excitons that can be applied to very large systems beyond reach for fully ab initio approaches.

Plasmon Excitations in Mixed Metallic Nanoarrays

Features of the surface plasmon from macroscopic materials emerge in molecular systems, but differentiating collective excitations from single-particle excitations in molecular systems remains elusive. The rich interactions between single-particle electron-hole and collective electron excitations produce phenomena related to the chemical physics aspects within the atomic array. We study the plasmonic properties of atomic arrays of noble (Au, Ag, and Cu) and transition-metal (Pd, Pt) homonuclear chains using time-dependent density functional theory and their Kohn-Sham transition contributions. The response to the electromagnetic radiation is related to both the geometry-dependent confinement of sp-valence electrons and the energy position of d-electrons in the different atomic species and the hybridization between d and sp electrons. It is possible to tune the position of the plasmon resonance, split it into several peaks, and eventually achieve broadband absorption of radiation. Arrays of mixed noble and transition-metal chains may have strongly attenuated plasmonic behavior. The collective nature of the excitations is ascertained using their Kohn-Sham transition contributions. To manipulate the plasmonic response and achieve the desired properties for broad applications, it is vital to understand the origins of these phenomena in atomic chains and their arrays.

Atomically Precise, Thiolated Copper-Hydride Nanoclusters as Single-Site Hydrogenation Catalysts for Ketones in Mild Conditions

Copper-hydrides are known catalysts for several technologically important reactions such as hydrogenation of CO, hydroamination of alkenes and alkynes, and chemoselective hydrogenation of unsaturated ketones to unsaturated alcohols. Stabilizing copper-based particles by ligand chemistry to nanometer scale is an appealing route to make active catalysts with optimized material economy; however, it has been long believed that the ligand-metal interface, particularly if sulfur-containing thiols are used as stabilizing agent, may poison the catalyst. We report here a discovery of an ambient-stable thiolate-protected copper-hydride nanocluster [Cu₂₅H₁₀(SPhCl₂)₁₈]^3- that readily catalyzes hydrogenation of ketones to alcohols in mild conditions. A full experimental and theoretical characterization of its atomic and electronic structure shows that the 10 hydrides are instrumental for the stability of the nanocluster and are in an active role being continuously consumed and replenished in the hydrogenation reaction. Density functional theory computations suggest, backed up by the experimental evidence, that the hydrogenation takes place only around a single site of the 10 hydride locations, rendering the [Cu₂₅H₁₀(SPhCl₂)₁₈]^3- one of the first nanocatalysts whose structure and catalytic functions are characterized fully to atomic precision. Understanding of a working catalyst at the atomistic level helps to optimize its properties and provides fundamental insights into the controversial issue of how a stable, ligand-passivated, metal-containing nanocluster can be at the same time an active catalyst.

Catalysis-Hub.org, an open electronic structure database for surface reactions

We present a new open repository for chemical reactions on catalytic surfaces, available at https://www.catalysis-hub.org . The featured database for surface reactions contains more than 100,000 chemisorption and reaction energies obtained from electronic structure calculations, and is continuously being updated with new datasets. In addition to providing quantum-mechanical results for a broad range of reactions and surfaces from different publications, the database features a systematic, large-scale study of chemical adsorption and hydrogenation on bimetallic alloy surfaces. The database contains reaction specific information, such as the surface composition and reaction energy for each reaction, as well as the surface geometries and calculational parameters, essential for data reproducibility. By providing direct access via the web-interface as well as a Python API, we seek to accelerate the discovery of catalytic materials for sustainable energy applications by enabling researchers to efficiently use the data as a basis for new calculations and model generation.

Silicon Substitution in Nanotubes and Graphene via Intermittent Vacancies

The chemical and electrical properties of single-walled carbon nanotubes (SWCNTs) and graphene can be modified by the presence of covalently bound impurities. Although this can be achieved by introducing chemical additives during synthesis, it often hinders growth and leads to limited crystallite size and quality. Here, through the simultaneous formation of vacancies with low-energy argon plasma and the thermal activation of adatom diffusion by laser irradiation, silicon impurities are incorporated into the lattice of both materials. After an exposure of ∼1 ion/nm², we find Si-substitution densities of 0.15 nm^-2 in graphene and 0.05 nm^-2 in nanotubes, as revealed by atomically resolved scanning transmission electron microscopy. In good agreement with predictions of Ar irradiation effects in SWCNTs, we find Si incorporated in both mono- and divacancies, with ∼2/3 being of the first type. Controlled inclusion of impurities in the quasi-1D and -2D carbon lattices may prove useful for applications such as gas sensing, and a similar approach might also be used to substitute other elements with migration barriers lower than that of carbon.

R-NEB: Accelerated Nudged Elastic Band Calculations by Use of Reflection Symmetry

Many activated processes in materials science, physics, and chemistry, e.g. diffusion processes, have initial and final states related by symmetry. Identification of minimum energy paths in such systems with methods such as nudged elastic band (NEB) can gain substantial speed up if the symmetry is exploited. The identification of minimum energy paths and transition states for such processes constitute a large fraction of the CPU-usage within computational materials science; much of which is in essence redundant due to symmetry. Paths with a reflection symmetry can be calculated using about half the computational resources, and the activation energy can, for some transitions, be estimated with high precision with a speed up factor equal to the number of images used in a standard NEB calculation. We present the formal properties required for a system to guarantee a reflection symmetric minimum energy path and an implementation to prepare and effectively speed up nudged elastic band calculations through symmetry considerations. Five examples are given to show the versatility and effectiveness of the method and to validate the implementation. The method is implemented in the open source package Atomic Simulation Environment (ASE) and contains internal methods to identify symmetry relations between the given end point configurations.

Low-Scaling Algorithm for Nudged Elastic Band Calculations Using a Surrogate Machine Learning Model

We present the incorporation of a surrogate Gaussian process regression (GPR) atomistic model to greatly accelerate the rate of convergence of classical nudged elastic band (NEB) calculations. In our surrogate model approach, the cost of converging the elastic band no longer scales with the number of moving images on the path. This provides a far more efficient and robust transition state search. In contrast to a conventional NEB calculation, the algorithm presented here eliminates any need for manipulating the number of images to obtain a converged result. This is achieved by inventing a new convergence criteria that exploits the probabilistic nature of the GPR to use uncertainty estimates of all images in combination with the force in the saddle point in the target model potential. Our method is an order of magnitude faster in terms of function evaluations than the conventional NEB method with no accuracy loss for the converged energy barrier values.