Theoretical investigation of finite size effects at DNA melting

We investigated how the finiteness of the length of a sequence affects the phase transition that takes place at the DNA melting temperature. For this purpose, we modified the transfer integral method to adapt it to the calculation of both extensive (partition function, entropy, specific heat, etc.) and nonextensive (order parameter and average separation between paired bases) thermodynamic quantities of finite sequences with open boundary conditions, and applied the modified procedure to two different dynamical models. We characterized in some detail the three effects that take place when the length of the sequence is decreased, namely, (i) the decrease of the critical temperature, (ii) the decrease of the peak values of all quantities that diverge at the thermodynamic limit but remain finite for finite sequences, like the specific heat and the correlation length, and (iii) the broadening of the temperature range over which the transition affects the dynamics of the system. We also performed a finite size scaling analysis of the two models and showed that the singular part of the free energy can indeed be expressed in terms of a homogeneous function. However, Josephson's identity is satisfied for none of the investigated models, so that the derivation of the characteristic exponents which appear, for example, in the expression of the specific heat requires some care.