Finite electric displacement simulations of polar ionic solid-electrolyte interfaces: Application to NaCl(111)/aqueous NaCl solution

Abstract

Tasker type III polar terminations of ionic crystals carry a net surface charge as well as a dipole moment and are fundamentally unstable. In contact with electrolytes, such polar surfaces can be stabilized by adsorption of counter ions from solution to form electric double layers (EDLs). In a previous work (J. Chem. Phys 147, 104702 (2017)) we reported on a classical force field based molecular dynamics study of a prototype model system namely a NaCl(111) slab interfaced with an aqueous NaCl solution on both sides. A serious hurdle in the simulation is that the finite width of the slab admits an electric field in the solid perturbing the theoretical charge balance at the interface of semi-infinite systems (half the surface charge density for NaCl(111)). It was demonstrated that the application of a finite macroscopic field $E$ cancelling the internal electric field can recover the correct charge compensation at the interface. In the present work, we expand this method by applying a conjugate electric displacement field $D$. The benefits of using $D$ instead of $E$ as the control variable are two fold: it does not only speed up the convergence of the polarization in the simulation but also leads to a succinct expression for the biasing displacement field involving only structural parameters which are known in advance. This makes it feasible to study the charge compensating phenomenon of this prototype system with density functional theory based molecular dynamics (DFTMD), as shown in this work.

arXiv comment13 pages, 8 figures

References in corpus (14)

• Canonical sampling through velocity-rescaling
• Separable Dual Space Gaussian Pseudo-potentials
• Maximally-localized generalized Wannier functions for composite energy bands
• Maximally localized Wannier functions: Theory and applications
• The Quantum-Mechanical Position Operator in Extended Systems
• Self-consistent hybrid functional for condensed systems
• Maximally-localized Wannier functions for disordered systems: application to amorphous silicon
• Electric displacement as the fundamental variable in electronic-structure calculations
• Recent Progress in Molecular Simulation of Aqueous Electrolytes: Force Fields, Chemical Potentials and Solubility
• Computing the dielectric constant of liquid water at constant dielectric displacement
• Finite Field Methods for the Supercell Modelling of Charged Insulator-Electrolyte Interfaces
• First-principles modeling of ferroelectric capacitors via constrained-D calculations
• Note: On the dielectric constant of nanoconfined water
• Charge compensation at the interface between the polar NaCl(111) surface and a NaCl aqueous solution

Cited by in corpus (8)

• Dielectric properties of nano-confined water: a canonical thermopotentiostat approach
• Coupling of surface chemistry and electric double layer at TiO${}_2$ electrochemical interfaces
• Computational Amperometry of Nanoscale Capacitors in Molecular Simulations
• Finite field formalism for bulk electrolyte solutions
• Stabilization of AgI's polar surfaces by the aqueous environment, and its implications for ice formation
• Learning classical density functionals for ionic fluids
• Macroscopic surface charges from microscopic simulations
• A theory for the stabilization of polar crystal surfaces by a liquid environment