Molecular simulation is becoming a very powerful tool for studying dynamic phenomena in materials. The simulation yields information about interaction at length and time scales unattainable by experimental measurements and unpredictable by continuum theories. This is especially meaningful when referring to bonding between a polymer and a metal substrate. A very important characteristic of polymers is that their physical properties do not rely on the detailed chemical structure of the molecular chains but only on their flexibility, and accordingly they will be able to adopt different conformations. In this paper, a molecular simulation of the bonding between vinyl ester polymer and steel is presented. Four different polymers with increasing chain lengths have been studied. Atomic co-ordinates are adjusted in order to reduce the molecular energy. Conformational changes in the macromolecules have been followed to obtain the polymer pair correlation function. Radius of gyration and end-to-end distance distributions of the individual chains have been used as a quantitative measurement of their flexibility. There exists a correlation between flexibility of the molecular chains and the energy of adhesion between the polymer and the metal substrate. Close contacts between the two materials are established at certain points but every atom up to a certain distance from the interface contributes to the total value of the adhesion energy of the system. (C) Koninklijke Brill NV, Leiden, 2008