Alessandro Coretti, Camille Bacon, Roxanne Berthin, Alessandra Serva, Laura Scalfi, Iurii Chubak, Kateryna Goloviznina, Matthieu Haefele, Abel Marin-Laflèche, Benjamin Rotenberg, Sara Bonella, Mathieu Salanne, Centre Européen de Calcul Atomique et Moléculaire (CECAM), École normale supérieure de Lyon (ENS de Lyon)-Commissariat à l'énergie atomique et aux énergies alternatives (CEA)-Centre National de la Recherche Scientifique (CNRS), Universität Wien, PHysicochimie des Electrolytes et Nanosystèmes InterfaciauX (PHENIX), Institut de Chimie du CNRS (INC)-Sorbonne Université (SU)-Centre National de la Recherche Scientifique (CNRS), Réseau sur le stockage électrochimique de l'énergie (RS2E), Aix Marseille Université (AMU)-Université Toulouse III - Paul Sabatier (UT3), Université de Toulouse (UT)-Université de Toulouse (UT)-Université de Haute-Alsace (UHA) Mulhouse - Colmar (Université de Haute-Alsace (UHA))-Collège de France (CdF (institution))-Université de Picardie Jules Verne (UPJV)-Ecole Nationale Supérieure de Chimie de Paris - Chimie ParisTech-PSL (ENSCP), Université Paris sciences et lettres (PSL)-Université Paris sciences et lettres (PSL)-Université de Pau et des Pays de l'Adour (UPPA)-Institut de Chimie du CNRS (INC)-Sorbonne Université (SU)-Centre National de la Recherche Scientifique (CNRS)-Institut National Polytechnique (Toulouse) (Toulouse INP), Université de Toulouse (UT)-Institut polytechnique de Grenoble - Grenoble Institute of Technology (Grenoble INP ), Université Grenoble Alpes (UGA)-Université Grenoble Alpes (UGA)-Nantes Université (Nantes Univ)-Université de Montpellier (UM)-Ecole Nationale Supérieure de Chimie de Montpellier (ENSCM), Université de Montpellier (UM), Laboratoire de Mathématiques et de leurs Applications [Pau] (LMAP), Université de Pau et des Pays de l'Adour (UPPA)-Centre National de la Recherche Scientifique (CNRS), Maison de la Simulation (MDLS), Université de Versailles Saint-Quentin-en-Yvelines (UVSQ)-Institut National de Recherche en Informatique et en Automatique (Inria)-Commissariat à l'énergie atomique et aux énergies alternatives (CEA)-Université Paris-Saclay-Centre National de la Recherche Scientifique (CNRS), Institut Universitaire de France (IUF), Ministère de l'Education nationale, de l’Enseignement supérieur et de la Recherche (M.E.N.E.S.R.), and European Project: 676629,H2020 Pilier Excellent Science,H2020-EINFRA-2015-1,EoCoE(2015)
International audience; Electrochemistry is central to many applications, ranging from biology to energy science. Studies now involve a wide range of techniques, both experimental and theoretical. Modeling and simulations methods, such as density functional theory or molecular dynamics, provide key information on the structural and dynamic properties of the systems. Of particular importance are polarization effects of the electrode/electrolyte interface, which are difficult to simulate accurately. Here, we show how these electrostatic interactions are taken into account in the framework of the Ewald summation method. We discuss, in particular, the formal setup for calculations that enforce periodic boundary conditions in two directions, a geometry that more closely reflects the characteristics of typical electrolyte/electrode systems and presents some differences with respect to the more common case of periodic boundary conditions in three dimensions. These formal developments are implemented and tested in MetalWalls, a molecular dynamics software that captures the polarization of the electrolyte and allows the simulation of electrodes maintained at a constant potential. We also discuss the technical aspects involved in the calculation of two sets of coupled degrees of freedom, namely the induced dipoles and the electrode charges. We validate the implementation, first on simple systems, then on the well-known interface between graphite electrodes and a room-temperature ionic liquid. We finally illustrate the capabilities of MetalWalls by studying the adsorption of a complex functionalized electrolyte on a graphite electrode.