The equilibrium polymerization of sulfur is investigated by Monte Carlo simulations. The potential energy model is based on density functional results for the cohesive energy, structural, and vibrational properties as well as reactivity of sulfur rings and chains ͓Part I, J. Chem. Phys. 118, 9257 ͑2003͔͒. Liquid samples of 2048 atoms are simulated at temperatures 450рTр850 K and Pϭ0 starting from monodisperse S 8 molecular compositions. Thermally activated bond breaking processes lead to an equilibrium population of unsaturated atoms that can change the local pattern of covalent bonds and allow the system to approach equilibrium. The concentration of unsaturated atoms and the kinetics of bond interchanges is determined by the energy ⌬E b required to break a covalent bond. Equilibrium with respect to the bond distribution is achieved for 15р⌬E b р21 kcal/mol over a wide temperature range (Tу450 K), within which polymerization occurs readily, with entropy from the bond distribution overcompensating the increase in enthalpy. There is a maximum in the polymerized fraction at temperature T max that depends on ⌬E b . This fraction decreases at higher temperature because broken bonds and short chains proliferate and, for T рT max , because entropy is less important than enthalpy. The molecular size distribution is described well by a Zimm-Schulz function, plus an isolated peak for S 8 . Large molecules are almost exclusively open chains. Rings tend to have fewer than 24 atoms, and only S 8 is present in significant concentrations at all T. The T dependence of the density and the dependence of polymerization fraction and degree on ⌬E b give estimates of the polymerization temperature T f ϭ450Ϯ20 K.