Using molecular dynamics (MD), we have studied the mechanism of heat accommodation between carbon dioxide clusters and monomers for temperatures and cluster size conditions that exist in homogeneous condensing supersonic expansion plumes. The work was motivated by our meso-scale direct simulation Monte Carlo and Bhatnagar-Gross-Krook based condensation simulations where we found that the heat accommodation model plays a key role in the near-field of the nozzle expansion particularly as the degree of condensation increases [R. Kumar, Z. Li, and D. Levin, Phys. Fluids 23, 052001 (2011)]. The heat released by nucleation and condensation and the heat removed by cluster evaporation can be transferred or removed from either the kinetic or translational modes of the carbon dioxide monomers. The molecular dynamics results show that the time required for gas-cluster interactions to establish an equilibrium from an initial state of non-equilibrium is less than the time step used in meso-scale analyses [R. Kumar, Z. Li, and D. Levin, Phys. Fluids 23, 052001 (2011)]. Therefore, the good agreement obtained between the measured cluster and gas number density and gas temperature profiles with the meso-scale modeling using the second energy exchange mechanism is not fortuitous but is physically based. Our MD simulations also showed that a dynamic equilibrium is established by the gas-cluster interactions in which condensation and evaporation processes take place constantly to and from a cluster.