We perform extensive molecular dynamics simulations of a charged polymer in a good solvent in the regime where the chain is collapsed. We analyze the dependence of the gyration radius R g on the reduced Bjerrum length l B and find two different regimes. In the first one, called a weak electrostatic regime,, which is consistent only with the predictions of the counterion-fluctuation theory. In the second one, called a strong electrostatic regime, we find R g ∼ l −1=5
B. To explain the novel regime we modify the counterion-fluctuation theory. DOI: 10.1103/PhysRevLett.117.147801 Introduction.-The conformational states of a flexible neutral polymer in different solvents are well known. It is extended in a good solvent due to favorable excluded volume interactions with the solvent molecules and collapses into a compact globule in a bad solvent [1][2][3]. In contrast, a flexible polyelectrolyte (PE)-charged polymer in the presence of counterions-undergoes an extended to collapsed transition in both good and bad solvents. Unlike neutral polymers, the conformations of a PE depend not only on the solvent quality, but also crucially on the interplay between electrostatic energy and translational entropy of counterions [4,5]. The strength of the electrostatic interactions depends on the charge density along the PE, which is quantified by the dimensionless Bjerrum length l B . For small charge density, counterions are dispersed away from the PE, and the chain is in an extended necklace conformation when in a good or theta solvent [3], and is collapsed into a compact globule in a bad solvent [3,6]. With increasing charge density, the PE attains an extended conformation, regardless of the solvent quality and counterions begin to condense onto the PE, renormalizing its charge density [3,7,8]. Further increase of the PE charge density results in an effective attraction between similarly charged monomers of the PE and it collapses into a globule conformation, independent of the solvent quality [4][5][6][9][10][11][12][13][14].