Many applications in chemistry, biology,
and energy storage/conversion
research rely on molecular simulations to provide fundamental insight
into structural and transport properties of materials with high ionic
concentrations. Whether the system is comprised entirely of ions,
like ionic liquids, or is a mixture of a polar solvent with a salt,
e.g., liquid electrolytes for battery applications, the presence of
ions in these materials results in strong local electric fields polarizing
solvent molecules and large ions. To predict properties of such systems
from molecular simulations often requires either explicit or mean-field
inclusion of the influence of polarization on electrostatic interactions.
In this manuscript, we review the pros and cons of different treatments
of polarization ranging from the mean-field approaches to the most
popular explicit polarization models in molecular dynamics simulations
of ionic materials. For each method, we discuss their advantages and
disadvantages and emphasize key assumptions as well as their adjustable
parameters. Strategies for the development of polarizable models are
presented with a specific focus on extracting atomic polarizabilities.
Finally, we compare simulations using polarizable and nonpolarizable
models for several classes of ionic systems, discussing the underlying
physics that each approach includes or ignores, implications for implementation
and computational efficiency, and the accuracy of properties predicted
by these methods compared to experiments.