Molecular dynamics simulations of liquid-liquid interfaces in an electric field: The water-1,2-dichloroethane interface

Equivalence of Charge Imbalance and External Electric Fields during Free Energy Calculations of Membrane Electroporation

Electric fields across lipid membranes play important roles in physiology, medicine, and biotechnology, rationalizing the wide interest in modeling transmembrane potentials in molecular dynamics simulations. Transmembrane potentials have been implemented with external electric fields or by imposing charge imbalance between the two water compartments of a stacked double-membrane system. We compare the two methods in the context of membrane electroporation, which involves a large change of membrane structure and capacitance. We show that, given that Ewald electrostatics are defined with tinfoil boundary conditions, the two methods lead to (i) identical potentials of mean force (PMFs) of pore formation and expansion at various potentials, demonstrating that the two methods impose equivalent driving forces for large-scale transitions at membranes, and (ii) to identical polarization of water within thin water wires or open pores, suggesting that the two methods furthermore impose equivalent local electric fields. Without tinfoil boundary conditions, effects from external fields on pore formation are spuriously suppressed or even removed. Together, our study shows that both methods, external fields and charge imbalance, are well suitable for studying large-scale transitions of lipid membranes that involve changes of membrane capacitance. However, using charge imbalance is technically more challenging for maintaining a constant transmembrane potential since it requires updating of the charge imbalance as the membrane capacitance changes.

Fully periodic, computationally efficient constant potential molecular dynamics simulations of ionic liquid supercapacitors

Molecular dynamics (MD) simulations of complex electrochemical systems, such as ionic liquid supercapacitors, are increasingly including the constant potential method (CPM) to model conductive electrodes at a specified potential difference, but the inclusion of CPM can be computationally expensive. We demonstrate the computational savings available in CPM MD simulations of ionic liquid supercapacitors when the usual non-periodic slab geometry is replaced with fully periodic boundary conditions. We show how a doubled cell approach, previously used in non-CPM MD simulations of charged interfaces, can be used to enable fully periodic CPM MD simulations. Using either a doubled cell approach or a finite field approach previously reported by others, fully periodic CPM MD simulations produce comparable results to the traditional slab geometry simulations with a nearly double speedup in computational time. Indeed, these savings can offset the additional cost of the CPM algorithm, resulting in periodic CPM MD simulations that are computationally competitive with the non-periodic, fixed charge equivalent simulations for the ionic liquid supercapacitors studied here.

Electromagnetic bioeffects: a multiscale molecular simulation perspective

Electromagnetic bioeffects remain an enigma from both the experimental and theoretical perspectives despite the ubiquitous presence of related technologies in contemporary life. Multiscale computational modelling can provide valuable insights into biochemical systems and predict how they will be perturbed by external stimuli. At a microscopic level, it can be used to determine what (sub)molecular scale reactions various stimuli might induce; at a macroscopic level, it can be used to examine how these changes affect dynamic behaviour of essential molecules within the crowded biomolecular milieu in living tissues. In this review, we summarise and evaluate recent computational studies that examined the impact of externally applied electric and electromagnetic fields on biologically relevant molecular systems. First, we briefly outline the various methodological approaches that have been employed to study static and oscillating field effects across different time and length scales. The practical value of such modelling is then illustrated through representative case-studies that showcase the diverse effects of electric and electromagnetic field on the main physiological solvent - water, and the essential biomolecules - DNA, proteins, lipids, as well as some novel biomedically relevant nanomaterials. The implications and relevance of the theoretical multiscale modelling to practical applications in therapeutic medicine are also discussed. Finally, we summarise ongoing challenges and potential opportunities for theoretical modelling to advance the current understanding of electromagnetic bioeffects for their modulation and/or beneficial exploitation in biomedicine and industry.

Potential dependence of the ionic structure at the ionic liquid/water interface studied using MD simulation

The structure at the electrochemical liquid/liquid interface between water (W) and trioctylmethylammonium bis(nonafluorobutanesulfonyl)amide, a hydrophobic ionic liquid (IL), was studied using molecular dynamics (MD) simulation in which the interfacial potential difference was controlled. On the IL side of the IL/W interface, ionic multilayers were found in the number density distribution of IL ions whereas monolayer-thick charge accumulation was found in the charge density distribution. This suggests that the potential screening is completed within the first ionic layer and the effect of overlayers on the potential is marginal. The W side of the interface showed the diffuse electric double layer as expected, and unexpectedly unveiled a density depletion layer, indicating that the IL surface is hydrophobic enough to be repelled by water. The IL ions in the first ionic layer showed anisotropic orientation even at the potential of zero charge, in which the polar moieties were oriented to the W side and the non-polar moieties preferred parallel orientation to the interface. When an electric field is applied across the interface so that the IL ions are more accumulated, the non-polar moieties changed the parallel preference to more oriented to the IL side due to the dipolar nature of the IL ions. The ionic orientations at the IL/W interface were compared with those at other two IL interfaces, the vacuum and graphene interfaces of the IL. The parallel preference of the non-polar moieties was similar at the IL/graphene interface but different from the perpendicular orientation toward the vacuum side at the IL/vacuum interface. The comparison suggests that water behaves like a wall that repels IL ions like a solid electrode.