Bulk ionic screening lengths from extremely large-scale molecular dynamics simulations

Recent experiments have reported anomalously large screening lengths of interactions between charged surfaces confining concentrated electrolytes and ionic liquids. Termed underscreening, this effect was ascribed to bulk properties of dense ionic systems. Herein, we study bulk ionic screening with extremely large-scale molecular dynamics simulations, allowing us to assess the range of distances relevant to the experiments. Our results yield two screening lengths satisfying distinct scaling relations. However, with an accuracy of 10^-5k_BT in interionic potentials of mean force, we find no signs of underscreening, suggesting that other than bulk effects might be at play in the experiments.