pH-Dependent Distribution of Functional Groups on Titanium-Based MXenes.

MXenes are a new rapidly developing class of two-dimensional materials with suitable properties for a broad range of applications. It has been shown that during synthesis of these materials the surfaces are usually functionalized by O, OH, and F and further suggested that controlling the surface allows controlling the material properties. However, a proper understanding of the surface structure is still missing, with a significant discrepancy between computational and experimental studies. Experiments consistently show formation of surfaces with mixed terminations, whereas computational studies point toward pure terminated surfaces. Here, we explain the formation of mixed functionalization on the surface of titanium-based two-dimensional carbides, Ti2C and Ti3C2, using a multiscale modeling scheme. Our scheme is based on calculating Gibbs free energy of formation by a combination of electronic structure calculations with cluster expansion and Monte Carlo simulations. Our calculations show formation of mixtures of O, OH, and F on the surface with the composition depending on pH, temperature, and the work function. On the other hand, our results also suggest a limited stable range of compositions, which challenges the paradigm of a high tunability of MXene properties.

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.

Plasmon-Induced Direct Hot-Carrier Transfer at Metal-Acceptor Interfaces.

Plasmon-induced hot-carrier transfer from a metal nanostructure to an acceptor is known to occur via two key mechanisms: (i) indirect transfer, where the hot carriers are produced in the metal nanostructure and subsequently transferred to the acceptor, and (ii) direct transfer, where the plasmons decay by directly exciting carriers from the metal to the acceptor. Unfortunately, an atomic-level understanding of the direct-transfer process, especially with regard to its quantification, remains elusive even though it is estimated to be more efficient compared to the indirect-transfer process. This is due to experimental challenges in separating direct from indirect transfer as both processes occur simultaneously at femtosecond time scales. Here, we employ time-dependent density-functional theory simulations to isolate and study the direct-transfer process at a model metal-acceptor (Ag147-Cd33Se33) interface. Our simulations show that, for a 10 fs Gaussian laser pulse tuned to the plasmon frequency, the plasmon formed in the Ag147-Cd33Se33 system decays within 10 fs and induces the direct transfer with a probability of about 40%. We decompose the direct-transfer process further and demonstrate that the direct injection of both electrons and holes into the acceptor, termed direct hot-electron transfer (DHET) and direct hot-hole transfer (DHHT), takes place with similar probabilities of about 20% each. Finally, effective strategies to control and tune the probabilities of DHET and DHHT processes are proposed. We envision our work to provide guidelines toward the design of metal-acceptor interfaces that enable more efficient plasmonic hot-carrier devices.

Alkali Postdeposition Treatment-Induced Changes of the Chemical and Electronic Structure of Cu(In,Ga)Se2 Thin-Film Solar Cell Absorbers: A First-Principle Perspective.

The effects of alkali postdeposition treatment (PDT) on the valence band structure of Cu(In,Ga)Se2 (CIGSe) thin-film solar cell absorbers are addressed from a first-principles perspective. In detail, experimentally derived hard X-ray photoelectron spectroscopy (HAXPES) data [ Handick , E. ; ACS Appl. Mater. Interfaces 2015 , 7 , 27414 - 27420 ] of the valence band structure of alkali-free and NaF/KF-PDT CIGSe are directly compared and fit by calculated density of states (DOS) of CuInSe2, its Cu-deficient counterpart CuIn5Se8, and different potentially formed secondary phases, such as KInSe2, InSe, and In2Se3. The DOSs are based on first-principles electronic structure calculations and weighted according to element-, symmetry-, and energy-dependent photoionization cross sections for the comparison to experimental data. The HAXPES spectra were recorded using photon energies ranging from 2 to 8 keV, allowing extraction of information from different sample depths. The analysis of the alkali-free CIGSe valence band structure reveals that it can best be described by a mixture of the DOS of CuInSe2 and CuIn5Se8, resulting in a stoichiometry slightly more Cu-rich than that of CuIn3Se5. The NaF/KF-PDT-induced changes in the HAXPES spectra for different alkali exposures are best reproduced by additional contributions from KInSe2, with some indications that the formation of a pronounced K-In-Se-type surface species might crucially depend on the amount of K available during PDT.

Kohn-Sham Decomposition in Real-Time Time-Dependent Density-Functional Theory: An Efficient Tool for Analyzing Plasmonic Excitations.

Electronic excitations can be efficiently analyzed in terms of the underlying Kohn-Sham (KS) electron-hole transitions. While such a decomposition is readily available in the linear-response time-dependent density-functional theory (TDDFT) approaches based on the Casida equations, a comparable analysis is less commonly conducted within the real-time-propagation TDDFT (RT-TDDFT). To improve this situation, we present here an implementation of a KS decomposition tool within the local-basis-set RT-TDDFT code in the free GPAW package. Our implementation is based on postprocessing of data that is readily available during time propagation, which is important for retaining the efficiency of the underlying RT-TDDFT to large systems. After benchmarking our implementation on small benzene derivatives by explicitly reconstructing the Casida eigenvectors from RT-TDDFT, we demonstrate the performance of the method by analyzing the plasmon resonances of icosahedral silver nanoparticles up to Ag561. The method provides a clear description of the splitting of the plasmon in small nanoparticles due to individual single-electron transitions as well as the formation of a distinct d-electron-screened plasmon resonance in larger nanoparticles.

Theory and Applications of Generalized Pipek-Mezey Wannier Functions.

The theory for the generation of Wannier functions within the generalized Pipek-Mezey approach (Lehtola, S.; Jónsson, H. J. Chem. Theory Comput. 2014, 10, 642) is presented and an implementation thereof is described. Results are shown for systems with periodicity in one, two, and three dimensions as well as isolated molecules. The generalized Pipek-Mezey Wannier functions (PMWF) are highly localized orbitals consistent with chemical intuition where a distinction is maintained between σ- and π-orbitals. The PMWF method is compared with the so-called maximally localized Wannier functions (MLWFs) that are frequently used for the analysis of condensed matter calculations. Whereas PMWFs maximize the localization criterion of Pipek and Mezey, MLWFs maximize that of Foster and Boys and have the disadvantage of mixing σ- and π-orbitals in many cases. The PMWF orbitals turn out to be as localized as the MLWF orbitals as evidenced by cross-comparison of the values of the PMWF and MLWF objective functions for the two types of orbitals. Our implementation in the atomic simulation environment (ASE) is compatible with various representations of the wave function, including real-space grids, plane waves, and linear combinations of atomic orbitals. The projector-augmented wave formalism for the representation of atomic core electrons is also supported. Results of calculations with the GPAW software are described here, but our implementation can also use output from other electronic structure software such as ABINIT, NWChem, and VASP.

Resistive Switching in All-Oxide Ferroelectric Tunnel Junctions with Ionic Interfaces.

Universal, giant and nonvolatile resistive switching is demonstrated for oxide tunnel junctions with ferroelectric PbZr0.2 Ti0.8 O3 , ferroelectric BaTiO3, and paraelectric SrTiO3 tunnel barriers. The effects are caused by reversible migration of oxygen vacancies between the tunnel barrier and bottom La2/3 Sr1/3 MnO3 electrode. The switching process, which is driven by large electric fields, is efficient down to a temperature of 5 K.

Quantized Evolution of the Plasmonic Response in a Stretched Nanorod.

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.

Impact of Ga-V Codoping on Interfacial Electron Transfer in Dye-Sensitized TiO2.

The improvement of charge transfer between an organic molecule and a semiconductor is an important and challenging goal in the fields of photovoltaics and photocatalysis. In this work, we present a time-dependent density functional theory investigation of the impact of Ga-V codoping of TiO2 on the excited-state electron injection from perylene-3-carboxylic acid. The doping is shown to raise the charge-transfer efficiency for the highest possible surface dye uptake by ∼16%. The strength of the effect depends on the dopant-pair-dye separation, dopant concentration, and distribution of Ga, V atoms in TiO2. The doping of the superficial level turns out to be more favorable than those in the bulk. The changes in electron injection dynamics are attributed to the modification of accepting semiconductor levels and hybridization profile between molecular and semiconductor states.

Nanoplasmonics simulations at the basis set limit through completeness-optimized, local numerical basis sets.

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.

Nonadiabatic Ehrenfest molecular dynamics within the projector augmented-wave method.

We derive equations for nonadiabatic Ehrenfest molecular dynamics within the projector augmented-wave (PAW) formalism. The discretization of the electrons is time-dependent as the augmentation functions depend on the positions of the nuclei. We describe the implementation of the Ehrenfest molecular dynamics equations within the real-space PAW method. We demonstrate the applicability of our method by studying the vibration of NaCl, the torsional rotation of H(2)C=NH(2)(+) in both the adiabatic and the nonadiabatic regimes, and the hydrogen bombardment of C(40)H(16).

All-electron time-dependent density functional theory with finite elements: time-propagation approach.

We present an all-electron method for time-dependent density functional theory which employs hierarchical nonuniform finite-element bases and the time-propagation approach. The method is capable of treating linear and nonlinear response of valence and core electrons to an external field. We also introduce (i) a preconditioner for the propagation equation, (ii) a stable way to implement absorbing boundary conditions, and (iii) a new kind of absorbing boundary condition inspired by perfectly matched layers.

Tunability of the optical absorption in small silver cluster-polymer hybrid systems.

We have calculated the absorption characteristics of different hybrid systems consisting of Ag, Ag(2), or Ag(3) atomic clusters and poly(methacrylic acid) using the time-dependent density-functional theory. The polymer is found to have an extensive structural-dependency on the spectral patterns of the hybrid systems relative to the bare clusters. The absorption spectrum can be "tuned" to the visible range for hybrid systems with an odd number of electrons per silver cluster, whereas for hybrid systems comprising an even number of electrons per silver cluster, the leading absorption edge can be shifted up to approximately 4.5 eV. The results give theoretical support to the experimental observations on the absorption in the visible range in metal cluster-polymer hybrid structures.

All-electron density functional theory and time-dependent density functional theory with high-order finite elements.

We present for static density functional theory and time-dependent density functional theory calculations an all-electron method which employs high-order hierarchical finite-element bases. Our mesh generation scheme, in which structured atomic meshes are merged to an unstructured molecular mesh, allows a highly nonuniform discretization of the space. Thus it is possible to represent the core and valence states using the same discretization scheme, i.e., no pseudopotentials or similar treatments are required. The nonuniform discretization also allows the use of large simulation cells, and therefore avoids any boundary effects.

Intrinsic n-type behavior in transparent conducting oxides: a comparative hybrid-functional study of In2O3, SnO2, and ZnO.

We present a comparative study of oxygen vacancies in In2O3, SnO2, and ZnO based on the hybrid-functional method within the density-functional theory (DFT). For In2O3 and SnO2, our results provide strong evidence of shallow donor states at oxygen vacancies. In comparison with the (semi)local exchange-correlation approximations in DFT, the hybrid-functional method strongly lowers the formation energy of the positive charge state and keeps that of the neutral state nearly intact. The trend is analyzed in terms of changes in lattice relaxation energies and in electron energy levels near the band gap. The existence of shallow donor states at oxygen vacancies and the consequent n-type conductivity are in line with experimental findings. The results invalidate some former theoretical interpretations based on standard DFT calculations.

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.

Photoabsorption spectra of small fullerenes and Si-heterofullerenes.

We study the spectral properties of two kinds of derivatives of the carbon fullerene C(60), small fullerenes and Si-heterofullerenes, by ab initio calculations. The principal method of study is the time-dependent density-functional theory in its full time-propagation form. C(20), C(28), C(32), C(36), and C(50), the most stable small fullerenes in the range of C(20)-C(50), are found to have characteristic features in their optical absorption spectra, originating from the geometry of the molecules in question. The comparison of measured and calculated absorption spectra is found to be a useful tool in differentiating between different, almost isoenergetic ground state structure candidates of small fullerenes. Substitutionally doped fullerenes are of interest due to their enhanced chemical reactivity. It is suggested that the doping degree can be obtained by studying the absorption spectra. For example, it is observed that the spectra gradually change when doping C(60) up to C(48)Si(12) so that absorption in the visible and near infrared regions increases.

Photoabsorption spectra of boron nitride fullerenelike structures.

Optical absorption spectra have been calculated for a series of boron nitride fullerenelike cage structures BnNn of sizes n=12-36. The method used is a real-time, real-space implementation of the time-dependent density-functional theory, involving the full time propagation of the time-dependent Kohn-Sham equations. The spectra are found to be a possible tool for distinguishing between different boron nitride fullerene species and isomers. The trends and differences in the spectra are found to be related to the general geometry of the molecules. Comparison between local-density and generalized-gradient approximations for electron exchange-correlation functionals shows that both of them produce essentially the same spectral characteristics.