Molecular dynamics (MD) simulations of tetramyristoyl cardiolipin (TMCL) and tetraoleoyl cardiolipin (TOCL) were carried out with the newly developed CHARMM lipid force field (FF), C36, and with head group charges q= −1 and −2. The surface areas per lipid, AL, for q=-1 are 126 ± 0.1 Å2 for TOCL and 111 ± .1 Å2 for TMCL. These are 1.8 times than those of the diacyl equivalents: 63 Å2 for dimyristoylphosphatidylcholine (DMPC) at 328K, and 69 Å2 for dioleoylphosphatidylcholine (DOPC) at 310K. Area compressibility, Ka, of TOCL equals 340 ± 40 dyn/cm, approximately 50% higher than experimentally obtained for DOPC (and most diacyl lipids); an experimental value for Ka for cardiolipins is not presently available. The areas and compressibilities for TOCL from the present simulation differ substantially from those obtained by Dahlberg and Maliniak using the FF of Berger et al. under the same conditions (AL= 99 Å2; Ka = 1100 dyn/cm). The origin of the differences appears to be in the ion binding to the surfaces of the cardiolipin bilayer. Under FF of Berger et al., ions bind closer to the carbonyl group in the lipid chain region whereas C36 CHARMM FF results ion binding closer to the negatively charged phosphate groups of the head group. Deuterium order parameter measurements are underway to determine which FF yields areas more representative of the fluid state of cardiolipin bilayers.