Search

Density functional theory with linear combination of atomic orbitals (LCAO) basis sets is useful for studying large atomic systems, especially when it comes to computationally highly demanding time-dependent dynamics. We have implemented the Ehrenfest molecular dynamics (ED) method with the approximate approach of Tomfohr and Sankey within the projector augmented-wave code GPAW. We apply this method to small molecules as well as larger periodic systems, and elucidate its limits, advantages, and disadvantages in comparison to the existing implementation of Ehrenfest dynamics with a real-space grid representation. For modest atomic velocities, LCAO-ED shows satisfactory accuracy at a much reduced computational cost. This method will be particularly useful for modeling ion irradiation processes that require large amounts of vacuum in the simulation cell.

Alloyed metal nanoparticles are a promising platform for plasmonically enabled hot-carrier generation, which can be used to drive photochemical reactions. Although the non-plasmonic component in these systems has been investigated for its potential to enhance catalytic activity, its capacity to affect the photochemical process favorably has been underexplored by comparison. Here, we study the impact of surface alloy species and concentration on hot-carrier generation in Ag nanoparticles. By first-principles simulations, we photoexcite the localized surface plasmon, allow it to dephase, and calculate spatially and energetically resolved hot-carrier distributions. We show that the presence of non-noble species in the topmost surface layer drastically enhances hot-hole generation at the surface at the expense of hot-hole generation in the bulk, due to the additional d-type states that are introduced to the surface. The energy of the generated holes can be tuned by choice of the alloyant, with systematic trends across the d-band block. Already low surface alloy concentrations have a large impact, with a saturation of the enhancement effect typically close to 75% of a monolayer. Hot-electron generation at the surface is hindered slightly by alloying but here an judicious choice of the alloy composition allows one to strike a balance between hot electrons and holes. In this context, it is also important to consider that increasing the alloy concentration broadens the localized surface plasmon resonance, and thus decreases hot-carrier generation overall. Our work underscores the promise of utilizing multicomponent nanoparticles to achieve enhanced control over plasmonic catalysis, and provides guidelines for how hot-carrier distributions can be tailored by designing the electronic structure of the surface through alloying., Comment: 9 pages, 6 figures

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 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 (BSE), 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 of GPU acceleration has been achieved with minor modifications of the GPAW code thanks to the CuPy library. We end the review with an outlook describing some future plans for GPAW.

Bimetallic nanostructures combining plasmonic and catalytic metals are promising for tailoring and enhancing plasmonic hot-carrier generation utilized in plasmonic catalysis. In this work, we study the plasmonic hot-carrier generation in noble metal nanoparticles (Ag, Au, Cu) with single-atom dopants (Ag, Au, Cu, Pd, Pt) with first-principles time-dependent density-functional theory calculations. Our results show that the local hot-carrier generation at the dopant atom is significantly altered by the dopant element while the plasmonic response of the nanoparticle as a whole is not significantly affected. In particular, hot holes at the dopant atom originate from the discrete d-electron states of the dopant. The energies of these d-electron states, and hence those of the hot holes, depend on the dopant element, which opens up the possibility to tune hot-carrier generation with suitable dopants., Comment: 6 pages, 4 figures

While direct hot-carrier transfer can increase photo-catalytic activity, it is difficult to discern experimentally and competes with several other mechanisms. To shed light on these aspects, here, we model from first principles hot-carrier generation across the interface between plasmonic nanoparticles and a CO molecule. The hot-electron transfer probability depends non-monotonically on the nanoparticle-molecule distance and can be effective at long distances, well outside the region of chemisorption; hot-hole transfer on the other hand is limited to shorter distances. These observations can be explained by the energetic alignment between molecular and nanoparticle states as well as the excitation frequency. The hybridization of the molecular orbitals is the key predictor for hot-carrier transfer in these systems, emphasizing the need to include the effects of ground state hybridization for accurate predictions. Finally, we show a non-trivial dependence of the hot-carrier distribution on the excitation energy, which could be exploited when optimizing photo-catalytic systems., Comment: 11 pages, 5 figures

Pd nanoalloys show great potential as hysteresis-free, reliable hydrogen sensors. Here, a multi-scale modeling approach is employed to determine optimal conditions for optical hydrogen sensing using the Pd-Au-H system. Changes in hydrogen pressure translate to changes in hydrogen content and eventually the optical spectrum. At the single particle level, the shift of the plasmon peak position with hydrogen concentration (i.e., the "optical" sensitivity) is approximately constant at 180 nm/c_H for nanodisk diameters >~ 100 nm. For smaller particles, the optical sensitivity is negative and increases with decreasing diameter, due to the emergence of a second peak originating from coupling between a localized surface plasmon and interband transitions. In addition to tracking peak position, the onset of extinction as well as extinction at fixed wavelengths is considered. We carefully compare the simulation results with experimental data and assess the potential sources for discrepancies. Invariably, the results suggest that there is an upper bound for the optical sensitivity that cannot be overcome by engineering composition and/or geometry. While the alloy composition has a limited impact on optical sensitivity, it can strongly affect H uptake and consequently the "thermodynamic" sensitivity and the detection limit. Here, it is shown how the latter can be improved by compositional engineering and even substantially enhanced via the formation of an ordered phase that can be synthesized at higher hydrogen partial pressures., Comment: 14 pages, 8 figures

Strong light-matter interactions facilitate not only emerging applications in quantum and non-linear optics but also modifications of materials properties. In particular the latter possibility has spurred the development of advanced theoretical techniques that can accurately capture both quantum optical and quantum chemical degrees of freedom. These methods are, however, computationally very demanding, which limits their application range. Here, we demonstrate that the optical spectra of nanoparticle-molecule assemblies, including strong coupling effects, can be predicted with good accuracy using a subsystem approach, in which the response functions of the different units are coupled only at the dipolar level. We demonstrate this approach by comparison with previous time-dependent density functional theory calculations for fully coupled systems of Al nanoparticles and benzene molecules. While the present study only considers few-particle systems, the approach can be readily extended to much larger systems and to include explicit optical-cavity modes., Comment: 9 pages, 4 figures

Electronic circular dichroism (ECD) is a powerful spectroscopical method for investigating chiral properties at the molecular level. ECD calculations with the commonly used linear-response time-dependent density functional theory (LR-TDDFT) framework can be prohibitively costly for large systems. To alleviate this problem, we present here an ECD implementation for the projector augmented-wave method in the real-time-propagation TDDFT (RT-TDDFT) framework in the open-source GPAW code. Our implementation supports both local atomic basis set and real-space finite-difference representations of wave functions. We benchmark our implementation against an existing LR-TDDFT implementation in GPAW for small chiral molecules. We then demonstrate the efficiency of our local atomic basis set implementation for a large hybrid nanocluster.

Metal nanoparticles are attractive for plasmon-enhanced generation of hot carriers, which may be harnessed in photochemical reactions. In this work, we analyze the coherent femtosecond dynamics of photon absorption, plasmon formation, and subsequent hot-carrier generation through plasmon dephasing using first-principles simulations. We predict the energetic and spatial hot-carrier distributions in small metal nanoparticles and show that the distribution of hot electrons is very sensitive to the local structure. Our results show that surface sites exhibit enhanced hot-electron generation in comparison to the bulk of the nanoparticle. While the details of the distribution depend on particle size and shape, as a general trend lower-coordinated surface sites such as corners, edges, and {100} facets exhibit a higher proportion of hot electrons than higher-coordinated surface sites such as {111} facets or the core sites. The present results thereby demonstrate how hot carriers could be tailored by careful design of atomic-scale structures in nanoscale systems., Comment: 10 pages, 4 figures

Strong light-matter interactions in both the single-emitter and collective strong coupling regimes attract significant attention due to emerging quantum and nonlinear optics applications, as well as opportunities for modifying material-related properties. Further exploration of these phenomena requires an appropriate theoretical methodology, which is demanding since polaritons are at the intersection between quantum optics, solid state physics and quantum chemistry. Fortunately, however, nanoscale polaritons can be realized in small plasmon-molecule systems, which in principle allows treating them using ab initio methods, although this has not been demonstrated to date. Here, we show that time-dependent density-functional theory (TDDFT) calculations can access the physics of nanoscale plasmon-molecule hybrids and predict vacuum Rabi splitting in a system comprising a few-hundred-atom aluminum nanoparticle interacting with one or several benzene molecules. We show that the cavity quantum electrodynamics approach holds down to resonators on the order of a few cubic nanometers, yielding a single-molecule coupling strength exceeding 200 meV due to a massive vacuum field value of 4.5 V/nm. In a broader perspective, our approach enables parameter-free in-depth studies of polaritonic systems, including ground state, chemical and thermodynamic modifications of the molecules in the strong-coupling regime, which may find important use in emerging applications such as cavity enhanced catalysis., Comment: 6 pages, 5 figures

We study the plasmonic properties of arrays of atomic chains which comprise noble (Cu, Ag, and Au) and transition (Pd, Pt) metal atoms using time-dependent density-functional theory. We show that the response to the electromagnetic radiation is related to both physics, the geometry-dependent confinement of sp-valence electrons, and chemistry, the energy position of d-electrons in the different atomic species and the hybridization between d and sp electrons. As a result it is possible to tune the position of the surface plasmon resonance, split it to several peaks, and eventually achieve broadband absorption of radiation. Mixing the arrays with transition metals can strongly attenuate the plasmonic behaviour. We analyze the origin of these phenomena and show that they arise from rich interactions between single-particle electron-hole and collective electron excitations. The tunability of the plasmonic response of arrays of atomic chains, which can be realized on solid surfaces, opens wide possibilities for their applications. In the present study we obtain guidelines how the desired properties can be achieved.

The real-time-propagation formulation of time-dependent density-functional theory (RT-TDDFT) is an efficient method for modeling the optical response of molecules and nanoparticles. Compared to the widely adopted linear-response TDDFT approaches based on, e.g., the Casida equations, RT-TDDFT appears, however, lacking efficient analysis methods. This applies in particular to a decomposition of the response in the basis of the underlying single-electron states. In this work, we overcome this limitation by developing an analysis method for obtaining the Kohn-Sham electron-hole decomposition in RT-TDDFT. We demonstrate the equivalence between the developed method and the Casida approach by a benchmark on small benzene derivatives. Then, we use the method for analyzing the plasmonic response of icosahedral silver nanoparticles up to Ag$_{561}$. Based on the analysis, we conclude that in small nanoparticles individual single-electron transitions can split the plasmon into multiple resonances due to strong single-electron-plasmon coupling whereas in larger nanoparticles a distinct plasmon resonance is formed., Comment: 11 pages, 3 figures

We present an approach for generating local numerical basis sets of improving accuracy for first-principles nanoplasmonics simulations within time-dependent density functional theory. The method is demonstrated for copper, silver, and gold nanoparticles that are of experimental interest but computationally demanding due to the semi-core d-electrons that affect their plasmonic response. The basis sets are constructed by augmenting numerical atomic orbital basis sets by truncated Gaussian-type orbitals generated by the completeness-optimization scheme, which is applied to the photoabsorption spectra of homoatomic metal atom dimers. We obtain basis sets of improving accuracy up to the complete basis set limit and demonstrate that the performance of the basis sets transfers to simulations of larger nanoparticles and nanoalloys as well as to calculations with various exchange-correlation functionals. This work promotes the use of the local basis set approach of controllable accuracy in first-principles nanoplasmonics simulations and beyond., Comment: 11 pages, 6 figures

Quantum aspects, such as electron tunneling between closely separated metallic nanoparticles, are crucial for understanding the plasmonic response of nanoscale systems. We explore quantum effects on the response of the conductively coupled metallic nanoparticle dimer. This is realized by stretching a nanorod, which leads to the formation of a narrowing atomic contact between the two nanorod ends. Based on first-principles time-dependent density-functional-theory calculations, we find a discontinuous evolution of the plasmonic response as the nanorod is stretched. This is especially pronounced for the intensity of the main charge-transfer plasmon mode. We show the correlation between the observed discontinuities and the discrete nature of the conduction channels supported by the formed atomic-sized junction., Comment: Main text: 6 pages, 2 figures; Supplemental Material: 5 pages, 4 figures

We observe using ab initio methods that localized surface plasmon resonances in icosahedral silver nanoparticles enter the asymptotic region already between diameters of 1 and 2 nm, converging close to the classical quasistatic limit around 3.4 eV. We base the observation on time-dependent density-functional theory simulations of the icosahedral silver clusters Ag$_{55}$ (1.06 nm), Ag$_{147}$ (1.60 nm), Ag$_{309}$ (2.14 nm), and Ag$_{561}$ (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 for the electron structure code GPAW within the scope of this work. We obtain good agreement with experimental data and modeled results, including photoemission and plasmon resonance. Moreover, we can extrapolate the ab initio results to the classical quasistatically modeled icosahedral clusters., Comment: 8 pages, 5 figures