Electrostatic interactions between charged species play a prominent role in determining structures and states of physical system, leading to important technological and biological applications. In coarse-grained simulations, accurate description of electrostatic interactions is crucial in addressing physical phenomena at larger spatial and longer temporal scales.

In this thesis, we implement ENUF method, an abbreviation for Ewald summation based on non-uniform fast Fourier transform technique, into dissipative particle dynamics (DPD) scheme. With determined suitable parameters, the computational complexity of ENUF-DPD method is approximately described as O(N logN). The ENUF-DPD method is further validated by investigating dependence of polyelectrolyte conformations on charge fraction of polyelectrolyte and counterion valency of added salts, and studying of specific binding structures of dendrimers on amphiphilic membranes.

In coarse-grained simulations, electrostatic interactions are either explicitly calculated with suitable methods, or implicitly included in effective potentials. The effect of treatment fashion of electrostatic interactions on phase behavior of [BMIM][PF6] ionic liquid (IL) is systematically investigated. Our systematic analyses show that electrostatic interactions should be incorporated explicitly in development of effective potentials, as well as in coarse-grained simulations to improve reliability of simulation results.

Detailed image of microscopic structures and orientations of [BMIM][PF6] at graphene and vacuum interfaces are investigated by using atomistic simulations. Imidazolium rings and alkyl side chains of [BMIM] lie preferentially flat on graphene surface. At IL-vacuum interface, ionic groups pack closely together to form polar domains, leaving alkyl side chains populated at interface and imparting hydrophobic character. With the increase of IL filmthickness, orientations of [BMIM] change gradually from dominant flat distributions along graphene surface to orientations where imidazolium rings are either parallel or perpendicular to IL-vacuum interface with tilted angles. The interfacial spatial ionic structural heterogeneity formed by ionic groups also contributes to heterogeneous dynamics in interfacial regions.