A QM/MD Coupling Method to Model the Ion-Induced Polarization of Graphene

Joshua D. Elliott,Alessandro Troisi,Paola Carbone
DOI: https://doi.org/10.1021/acs.jctc.0c00239
2020-07-09
Journal of Chemical Theory and Computation
Abstract:We report a new Quantum Mechanical/Molecular Dynamics (QM/MD) simulation loop to model the coupling between the electron and atom dynamics in solid/liquid interfacial systems. The method can describe simultaneously both the quantum mechanical surface polarizability emerging from the proximity to the electrolyte and the electrolyte structure and dynamics. In the current setup, Density Functional Tight Binding calculations for the electronic structure calculations of the surface are coupled with classical molecular dynamics to simulate the electrolyte solution. The reduced computational cost of the QM part makes the coupling with a classical simulation engine computationally feasible and allows simulation of large systems for hundreds of nanoseconds. We tested the method by simulating both a noncharged graphene flake and a noncharged and charged infinite graphene sheet immersed in an NaCl electrolyte solution. We found that, when no bias is applied, ions preferentially remained in solution, and only cations are mildly attracted to the surface of the graphene. This preferential adsorption of cations vs anions seems to persist also when the surface is moderately charged and rules out any substantial ions/surface charge transfer.The Supporting Information is available free of charge at <a class="ext-link" href="/doi/10.1021/acs.jctc.0c00239?goto=supporting-info">https://pubs.acs.org/doi/10.1021/acs.jctc.0c00239</a>.Plots showing convergence of DFTB Mulliken charges and change in average computed DFTB Mulliken charge over classical MD trajectory with different time steps (<a class="ext-link" href="/doi/suppl/10.1021/acs.jctc.0c00239/suppl_file/ct0c00239_si_001.pdf">PDF</a>)Movies of time evolution of surface charge for neutral surfaces (<a class="ext-link" href="/doi/suppl/10.1021/acs.jctc.0c00239/suppl_file/ct0c00239_si_002.avi">AVI</a>)Movies of time evolution of surface charge for cathodic (negatively charged) surfaces (<a class="ext-link" href="/doi/suppl/10.1021/acs.jctc.0c00239/suppl_file/ct0c00239_si_003.avi">AVI</a>)Movies of time evolution of surface charge for anodic (positively charged) surfaces (<a class="ext-link" href="/doi/suppl/10.1021/acs.jctc.0c00239/suppl_file/ct0c00239_si_004.avi">AVI</a>)This article has not yet been cited by other publications.
chemistry, physical,physics, atomic, molecular & chemical
What problem does this paper attempt to address?