Interactions of like-charged rods at low temperatures: Analytical theory vs. simulations

We investigate a system consisting of two like-charged infinitely long rods and neutralizing counterions at low temperatures, using both analytic theory and simulations. With some reasonable approximations we can analytically solve for several ground-state structures of the model, starting with states where all counterions are lined up in the gap between the rods, over planar configurations, where the counterions are divided up into a fraction which resides between the rods, and counterions which are located on the outer surfaces, up to configurations which cover the full rod surfaces. Using parallel tempering simulations, we are able to study the system over a wide range of temperatures. At low temperatures we find good agreement with our T = 0 results. At higher temperatures, the strong coupling (SC) theory delivers qualitatively better results. We furthermore demonstrate that for the SC theory and our ground-state approximations to yield quantitative agreement, three parameters are required to be large, the strong-coupling parameter Ξ, the Rouzina-Bloomfield parameter, and the ratio of the average distance of the counterions to the radius of the rods. In the case of the latter ratio being small, our T = 0 results show better agreement with the simulation data at very low temperatures.