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.

Noble metal nanoparticles stabilized by organic ligands are important for applications in assembly, site-specific bioconjugate labelling and sensing, drug delivery and medical therapy, molecular recognition and molecular electronics, and catalysis. Here we report crystal structures and theoretical analysis of three Ag 44 (SR) 30 and three Au 12 Ag 32 (SR) 30 intermetallic nanoclusters stabilized with fluorinated arylthiols (SR ¼ SPhF, SPhF 2 or SPhCF 3 ). The nanocluster forms a Keplerate solid of concentric icosahedral and dodecahedral atom shells, protected by six Ag 2 (SR) 5 units. Positive counterions in the crystal indicate a high negative charge of 4 À per nanoparticle, and density functional theory calculations explain the stability as an 18-electron superatom shell closure in the metal core. Highly featured optical absorption spectra in the ultraviolet-visible region are analysed using time-dependent density functional perturbation theory. This work forms a basis for further understanding, engineering and controlling of stability as well as electronic and optical properties of these novel nanomaterials.

Gold nanoclusters protected by a thiolate monolayer (MPC) are widely studied for their potential applications in site-specific bioconjugate labeling, sensing, drug delivery, and molecular electronics. Several MPCs with 1-2 nm metal cores are currently known to have a well-defined molecular structure, and they serve as an important link between molecularly dispersed gold and colloidal gold to understand the size-dependent electronic and optical properties. Here, we show by using an ab initio method together with atomistic models for experimentally observed thiolate-stabilized gold clusters how collective electronic excitations change when the gold core of the MPC grows from 1.5 to 2.0 nm. A strong localized surface plasmon resonance (LSPR) develops at 540 nm (2.3 eV) in a cluster with a 2.0 nm metal core. The protecting molecular layer enhances the LSPR, while in a smaller cluster with 1.5 nm gold core, the plasmon-like resonance at 540 nm is confined in the metal core by the molecular layer. Our results demonstrate a threshold size for the emergence of LSPR in these systems and help to develop understanding of the effect of the molecular overlayer on plasmonic properties of MPCs enabling engineering of their properties for plasmonic applications.

We observe using ab initio methods that localized surface plasmon resonances in icosahedral silver nanoparticles enter the asymptotic region already between diameters of 1-2 nm, converging close to the classical quasistatic limit around 3.4 eV. We base the observation on time-dependent densityfunctional theory simulations of the icosahedral silver clusters Ag55 (1.06 nm), Ag147 (1.60 nm), Ag309 (2.14 nm), and Ag561 (2.68 nm). The simulation method combines the adiabatic GLLB-SC exchange-correlation functional with real time propagation in an atomic orbital basis set using the projector augmented wave method. The method has been implemented to the electron structure code GPAW within the scope of this work. We obtain good agreement with experimental data and modelled results, including photoemission and plasmon resonance. Moreover we can extrapolate the ab initio results to the classical quasistatically modelled icosahedral clusters.

Through a systematic structural search we found an allotrope of carbon with Cmmm symmetry which we predict to be more stable than graphite for pressures above 10 GPa. This material, which we refer to as Z-carbon, is formed by pure sp(3) bonds and it provides an explanation to several features in experimental x-ray diffraction and Raman spectra of graphite under pressure. The transition from graphite to Z-carbon can occur through simple sliding and buckling of graphene sheets. Our calculations predict that Z-carbon is a transparent wide band-gap semiconductor with a hardness comparable to diamond.

Articles you may be interested inMixed time-dependent density-functional theory/classical trajectory surface hopping study of oxirane photochemistry J. Chem. Phys. 129, 124108 (2008) We present the implementation of the time-dependent density-functional theory both in linear-response and in time-propagation formalisms using the projector augmented-wave method in real-space grids. The two technically very different methods are compared in the linear-response regime where we found perfect agreement in the calculated photoabsorption spectra. We discuss the strengths and weaknesses of the two methods as well as their convergence properties. We demonstrate different applications of the methods by calculating excitation energies and excited state Born-Oppenheimer potential surfaces for a set of atoms and molecules with the linear-response method and by calculating nonlinear emission spectra using the time-propagation method.

Determining the structures of nanoparticles at atomic resolution is vital to understand their structure–property correlations. Large metal nanoparticles with core diameter beyond 2 nm have, to date, eluded characterization by single-crystal X-ray analysis. Here we report the chemical syntheses and structures of two giant thiolated Ag nanoparticles containing 136 and 374 Ag atoms (that is, up to 3 nm core diameter). As the largest thiolated metal nanoparticles crystallographically determined so far, these Ag nanoparticles enter the truly metallic regime with the emergence of surface plasmon resonance. As miniatures of fivefold twinned nanostructures, these structures demonstrate a subtle distortion within fivefold twinned nanostructures of face-centred cubic metals. The Ag nanoparticles reported in this work serve as excellent models to understand the detailed structure distortion within twinned metal nanostructures and also how silver nanoparticles can span from the molecular to the metallic regime.

We predict and analyze density-functional theory (DFT)-based structures for the recently isolated Au(40)(SR)(24) cluster. Combining structural information extracted from ligand-exchange reactions, circular dichroism and transmission electron microscopy leads us to propose two families of low-energy structures that have a chiral Au-S framework on the surface. These families have a common geometrical motif where a nonchiral Au(26) bi-icosahedral cluster core is protected by 6 RS-Au-SR and 4 RS-Au-SR-Au-SR oligomeric units, analogously to the "Divide and Protect" motif of known clusters Au(25)(SR)(18)(-/0), Au(38)(SR)(24) and Au(102)(SR)(44). The strongly prolate shape of the proposed Au(26) core is supported by transmission electron microscopy. Density-of-state-analysis shows that the electronic structure of Au(40)(SR)(24) can be interpreted in terms of a dimer of two 8-electron superatoms, where the 8 shell electrons are localized at the two icosahedral halves of the metal core. The calculated optical and chiroptical characteristics of the optimal chiral structure are in a fair agreement with the reported data for Au(40)(SR)(24).