Monte Carlo simulation of an inhomogeneous dielectric continuum model for B-DNA

Scott Grandison, Robert Penfold, Jean-Marc Vanden-Broeck

Research output: Contribution to journalArticlepeer-review

7 Citations (Scopus)


Thermodynamic and structural properties of the counterion atmosphere surrounding B-DNA are calculated by Monte Carlo simulation in a spatially inhomogeneous, but piecewise uniform, dielectric continuum cell model — the “barbarous” model. A boundary element formulation is implemented to study the sensitivity of these properties with respect to perturbations in the location of discontinuous dielectric boundaries relative to fixed and mobile charges. High concentrations are considered corresponding to the liquid crystalline hexagonally ordered phase of DNA. Primitive model results are verified against other simulation reports and a comparison of barbarous model predictions with experimental data is discussed. The internal energy, osmotic coefficient, radial distributions and the population ratio of counterions in the geometrically resolved major and minor grooves are all found to strongly depend on the dielectric boundary position. This suggests that a self-consistent development of the model should consider a free surface problem where the boundary is not specified a priori.
Original languageEnglish
Pages (from-to)3486-3495
Number of pages10
JournalPhysical Chemistry Chemical Physics
Issue number19
Publication statusPublished - 2005

Cite this