GPAW: An open Python package for electronic structure calculations.

We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for the implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE), providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe-Salpeter Equation, variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn-Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support for graphics processing unit (GPU) acceleration has been achieved with minor modifications to the GPAW code thanks to the CuPy library. We end the review with an outlook, describing some future plans for GPAW.

Birth of the localized surface plasmon resonance in monolayer-protected gold nanoclusters.

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.

Time-dependent density-functional theory in the projector augmented-wave method.

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.