A Practical Approach to Large Scale Electronic Structure Calculations in Electrolyte Solutions via Continuum-Embedded Linear-Scaling DFT

29 January 2020, Version 1
This content is a preprint and has not undergone peer review at the time of posting.

Abstract

We present the implementation of a hybrid continuum-atomistic model for including the effects of surrounding electrolyte in large-scale density functional theory (DFT) calculations within the ONETEP linear-scaling DFT code, which allows the simulation of large complex systems such as electrochemical interfaces. The model represents the electrolyte ions as a scalar field and the solvent as a polarisable dielectric continuum, both surrounding the quantum solute. The overall energy expression is a grand canonical functional incorporating the electron kinetic and exchange correlation energies, the total electrostatic energy, entropy and chemical potentials of surrounding electrolyte, osmotic pressure, and the effects of cavitation, dispersion and repulsion. The DFT calculation is performed fully self-consistently in the electrolyte model, allowing the quantum mechanical system and the surrounding continuum environment to interact and mutually polarize. A bespoke parallel Poisson-Boltzmann solver library, DL_MG, deals with the electrostatic problem, solving a generalized Poisson-Boltzmann equation. Our model supports open boundary conditions, which allows the treatment of molecules, entire biomolecules or larger nanoparticle assemblies in electrolyte. We have also implemented the model for periodic boundary conditions, allowing the treatment of extended systems such as electrode surfaces in contact with electrolyte. A key feature of the model is the use of solute-size and solvation-shell-aware accessibility functions that prevent the unphysical accumulation of electrolyte charge near the quantum solute boundary. The model has a small number of parameters: here we demonstrate their calibration against experimental mean activity coefficients. We also present an exemplar simulation of a 1634-atom model of the interface between a graphite anode and LiPF6 electrolyte in ethylene carbonate solvent. We compare the cases where Li atoms are intercalated at opposite edges of the graphite slab and in solution, demonstrating a potential application of the model in simulations of fundamental processes in Li-ion batteries.

Keywords

DFT+PB
continuum solvent and electrolyte models
electrode-electrolyte interface
Poisson Boltzmann

Comments

Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.