The computational study of chemical reactions in complex, wet environments is critical for applications in many fields. It is often essential to study chemical reactions in the presence of applied electrochemical potentials, taking into account the non-trivial electrostatic screening coming from the solvent and the electrolytes. As a consequence, the electrostatic potential has to be found by solving the generalized Poisson and the Poisson-Boltzmann equations for neutral and ionic solutions, respectively. In the present work, solvers for both problems have been developed. A preconditioned conjugate gradient method has been implemented for the solution of the generalized Poisson equation and the linear regime of the Poisson-Boltzmann, allowing to solve iteratively the minimization problem with some ten iterations of the ordinary Poisson equation solver. In addition, a self-consistent procedure enables us to solve the non-linear Poisson-Boltzmann problem. Both solvers exhibit very high accuracy and parallel efficiency and allow for the treatment of periodic, free, and slab boundary conditions. The solver has been integrated into the BigDFT and Quantum-ESPRESSO electronic-structure packages and will be released as an independent program, suitable for integration in other codes. C 2016 AIP Publishing LLC. [http://dx.

The behavior of metals across a broad frequency range from microwave to ultraviolet frequencies is of interest in plasmonics, nanophotonics, and metamaterials. Depending on the frequency, losses of collective excitations in metals can be predominantly classical resistive effects or Landau damping. In this context, we present first-principles calculations that capture all of the significant microscopic mechanisms underlying surface plasmon decay and predict the initial excited carrier distributions so generated. Specifically, we include ab initio predictions of phonon-assisted optical excitations in metals, which are critical to bridging the frequency range between resistive losses at low frequencies and direct interband transitions at high frequencies. In the commonly used plasmonic materials, gold, silver, copper, and aluminum, we find that resistive losses compete with phonon-assisted carrier generation below the interband threshold, but hot carrier generation via direct transitions dominates above threshold. Finally, we predict energy-dependent lifetimes and mean free paths of hot carriers, accounting for electron-electron and electron-phonon scattering, to provide insight toward transport of plasmonically generated carriers at the nanoscale.

Decay of surface plasmons to hot carriers finds a wide variety of applications in energy conversion, photocatalysis and photodetection. However, a detailed theoretical description of plasmonic hot-carrier generation in real materials has remained incomplete. Here we report predictions for the prompt distributions of excited ‘hot’ electrons and holes generated by plasmon decay, before inelastic relaxation, using a quantized plasmon model with detailed electronic structure. We find that carrier energy distributions are sensitive to the electronic band structure of the metal: gold and copper produce holes hotter than electrons by 1–2 eV, while silver and aluminium distribute energies more equitably between electrons and holes. Momentum-direction distributions for hot carriers are anisotropic, dominated by the plasmon polarization for aluminium and by the crystal orientation for noble metals. We show that in thin metallic films intraband transitions can alter the carrier distributions, producing hotter electrons in gold, but interband transitions remain dominant.

Energy and environmental concerns demand development of more efficient and selective electrodes for electrochemical reduction of CO2 to form fuels and chemicals. Since Cu is the only pure metal exhibiting reduction to form hydrocarbon chemicals, we focus here on the Cu (111) electrode. We present a methodology for density functional theory calculations to obtain accurate onset electrochemical potentials with explicit constant electrochemical potential and pH effects using implicit solvation. We predict the atomistic mechanisms underlying electrochemical reduction of CO, finding that (1) at acidic pH, the C1 pathway proceeds through COH to CHOH to form CH4 while C2 (C3) pathways are kinetically blocked; (2) at neutral pH, the C1 and C2 (C3) pathways share the COH common intermediate, where the branch to C-C coupling is realized by a novel CO-COH pathway; and (3) at high pH, early C-C coupling through adsorbed CO dimerization dominates, suppressing the C1 pathways by kinetics, thereby boosting selectivity for multi-carbon products.

Nanoscale localization of electromagnetic fields near metallic nanostructures underpins the fundamentals and applications of plasmonics. The unavoidable energy loss from plasmon decay, initially seen as a detriment, has now expanded the scope of plasmonic applications to exploit the generated hot carriers. However, quantitative understanding of the spatial localization of these hot carriers, akin to electromagnetic near-field maps, has been elusive. Here we spatially map hot-electron-driven reduction chemistry with 15 nm resolution as a function of time and electromagnetic field polarization for different plasmonic nanostructures. We combine experiments employing a six-electron photo-recycling process that modify the terminal group of a self-assembled monolayer on plasmonic silver nanoantennas, with theoretical predictions from first-principles calculations of non-equilibrium hot-carrier transport in these systems. The resulting localization of reactive regions, determined by hot-carrier transport from high-field regions, paves the way for improving efficiency in hot-carrier extraction science and nanoscale regio-selective surface chemistry.

Delivering the full benefits of first principles calculations to battery materials demands the development of accurate and computationally-efficient electronic structure methods that incorporate the effects of the electrolyte environment and electrode potential.Realistic electrochemical interfaces containing polar surfaces are beyond the regime of validity of existing continuum solvation theories developed for molecules, due to the presence of significantly stronger electric fields. We present an ab initio theory of the nonlinear dielectric and ionic response of solvent environments within the framework of joint densityfunctional theory, with precisely the same optimizable parameters as conventional polarizable continuum models. We demonstrate that the resulting nonlinear theory agrees with the standard linear models for organic molecules and metallic surfaces under typical operating conditions. However, we find that the saturation effects in the rotational response of polar solvent molecules, inherent to our nonlinear theory, are crucial for a qualitatively correct description of the ionic surfaces typical of the solid electrolyte interface. arXiv:1301.6189v1 [cond-mat.mtrl-sci]

First-principles calculations combining density-functional theory and continuum solvation models enable realistic theoretical modeling and design of electrochemical systems. When a reaction proceeds in such systems, the number of electrons in the portion of the system treated quantum mechanically changes continuously, with a balancing charge appearing in the continuum electrolyte. A grand-canonical ensemble of electrons at a chemical potential set by the electrode potential is therefore the ideal description of such systems that directly mimics the experimental condition. We present two distinct algorithms, a self-consistent field method (GC-SCF) and a direct variational free energy minimization method using auxiliary Hamiltonians (GC-AuxH), to solve the Kohn-Sham equations of electronic density-functional theory directly in the grand canonical ensemble at fixed potential. Both methods substantially improve performance compared to a sequence of conventional fixed-number calculations targeting the desired potential, with the GC-AuxH method additionally exhibiting reliable and smooth exponential convergence of the grand free energy. Finally, we apply grand-canonical DFT to the under-potential deposition of copper on platinum from chloride-containing electrolytes and show that chloride desorption, not partial copper monolayer formation, is responsible for the second voltammetric peak.Density-functional theory (DFT) enables theoretical elucidation of reaction mechanisms at complex catalyst surfaces, making it now possible to design efficient heterogeneous catalysts for various industrial applications from first principles, for example for high-temperature gas-phase transformation of hydrocarbons to a variety of valuable chemical products. 1,2 The extension of this predictive power to electrocatalysis would be highly valuable for an even broader class of technological problems, including a cornerstone of future technology for renewable energy: converting solar energy to chemical fuels by electrochemical water splitting and carbon dioxide reduction. 3 Accurately describing electrochemical phenomena, however, presents two additional challenges.First, the electrolyte, typically consisting of ions in a liquid solvent, strongly affects the energetics of structures and reactions at the interface. Treating liquids directly in DFT requires expensive molecular dynamics to sample the thermodynamic phase space of atomic configurations. Historically, a number of continuum solvation models that empirically capture liquid effects have enabled theoretical design of liquid-phase catalysts. 4,5 More recently, empirical solvation models suitable for solid-liquid interfaces, 6-8 joint density-functional theory (JDFT) for efficiently treating liquids with atomic-scale structure, 9 and minimally-empirical solvation models derived from JDFT, 10,11 have made great strides towards a) Electronic mail: sundar@rpi.edu b) Electronic mail: wag@wag.caltech.edu c) Electronic mail: taa2@cornell.edu reliable yet efficient treatment of electrochemica...

Density-functional theory (DFT) has revolutionized computational prediction of atomic-scale properties from first principles in physics, chemistry and materials science. Continuing development of new methods is necessary for accurate predictions of new classes of materials and properties, and for connecting to nano- and mesoscale properties using coarse-grained theories. JDFTx is a fully-featured open-source electronic DFT software designed specifically to facilitate rapid development of new theories, models and algorithms. Using an algebraic formulation as an abstraction layer, compact C++11 code automatically performs well on diverse hardware including GPUs (Graphics Processing Units). This code hosts the development of joint density-functional theory (JDFT) that combines electronic DFT with classical DFT and continuum models of liquids for first-principles calculations of solvated and electrochemical systems. In addition, the modular nature of the code makes it easy to extend and interface with, facilitating the development of multi-scale toolkits that connect to ab initio calculations, e.g. photo-excited carrier dynamics combining electron and phonon calculations with electromagnetic simulations.