Abstract

This work explores the use of continuous thermodynamic integration in field-theoretic simulations of a symmetric diblock copolymer melt. Free energies of the lamellar and disorder phases are evaluated by thermodynamic integration from a reference state (an Einstein crystal, λ = 0) to a diblock copolymer (λ = 1). This is followed by integration over the interaction parameter, χb, to locate the order–disorder transition (ODT). Then the equilibrium lamellar spacing and free energy cost of stretching and compressing lamellae are examined. The ODT, lamellar spacing, and compression modulus are consistent with previous calculations, though found faster and more precisely. The above quantities do not depend on simulation box size, suggesting that finite-size effects are small and simulating two lamellar periods is sufficient to accurately evaluate bulk behavior. Furthermore, the statistical uncertainty in the ODT increases quickly with system size, suggesting that small systems may lead to more precise results.