Large-scale first-principles molecular dynamics for electrochemical systems with O(N) methods

A method for large-scale first-principles molecular dynamics (MD) simulations on electrochemical systems has been developed by combining the effective screening medium (ESM) method with O(N) density functional theory (DFT). This implementation has been significantly simplified by the introduction of neutral atom potentials, which minimizes the modifications to existing DFT code. In order to demonstrate ability of this implementation, it has been applied to an electrochemical system consisting of a H-Si(111) electrode, which is a candidate anode for high-capacity Li-ion secondary batteries, and a propylene carbonate (PC) solvent to simulate how PC molecules in the vicinity of the electrode surface respond to an imposed electric field. The large-scale MD simulation clearly demonstrates that the combination of the ESM and O(N) DFT methods provides a useful tool for first-principles investigation of complicated electrochemical systems such as high-capacity batteries.