G. Elastic fluctuations in the mesoscopic model
In order to obtain a full description of the model that includes the effects of the elastic fluctuations of the polymers, we have to split the contribution of the handles and the ssDNA, so that they can have independent lengths (see Fig. G.1).
Full model. Each element is treated indepentdently.
The total energy of the system is then written as,
where the explicit dependencies of the handles are given by Eq. F.2; either Eq. F.2 or Eq. F.5 give the dependency of the ssDNA; Eq. 3.17 gives the free energy of the duplex and Eq. 3.16 gives the energy of the bead in the optical trap.
The total extension of the system is the control parameter and it can be written as
is the diameter of the DNA duplex and is the diameter of the bead. The partition function can be calculated by summing over all the possible states of the system according to
where the integral is extended over all positive values of the extensions
and ; provided that Eq. G.2 is fulfilled. Now, the calculation of this integral is quite complex because the extensions in Eq. G.3 are restricted by Eq. G.2, so they are not independent variables. Our goal is to write the total free energy of the system in terms of independent variables that can be integrated separately. So we can define a new set of positively defined variables , , and (see Fig. G.1) that indicate the position of the ends of the elastic elements, which are related to the extensions according to the following relations:
and allow us to rewrite Eq. G.1 as
Now the total energy of the system can be calculated for different extensions of the elements without affecting the total extension of the system. For instance, we can give different values to that modify the extensions of and without changing the value of
. Now, the partition function can be calculated by integrating over the new variables of the elements according to
where the integral is extended over all values of
while keeping the value of
constant. Note that Eqs. G.4 define a set of limits of integration for the new variables that do not coincide with the ones used in Eq. G.6. Nevertheless, we extend the limits of integration of the newly defined variables over all the space. We can do this because the contribution to the integral vanishes for large extensions, due to the exponential dependence of the energy. The advantage of such approximation is that we will have to perform Gaussian integrals, which are easy to calculate when the limits of integration expand over all the space.
For each value of
and , there is always a combination of
that minimizes the total energy in Eq. G.5. It is convenient to define this global minimum as
because it allows to perform the integral by using a saddle-point approximation. If we define
as the vector that contains all the variables, we can perform a Taylor expansion of Eq. G.5 around the minimum
up to order
stands for the transposed vector and
is the Hessian matrix calculated at fixed
and . Formally, the Hessian matrix is a quadratic form that acts over the vector
and its transposed
and returns the change in
. Note that the variables of the function are , while
and are fixed parameters. Since the Taylor expansion has been performed around the minimum, the gradient (i.e., first derivative) of the function vanishes and the second term in the right side of Eq. G.7 does not contribute.
In order to obtain the Hessian matrix, we have to calculate all the second order derivatives of Eq. G.1. The first derivatives of elastic energy terms give forces and the second derivatives give stiffnesses. So for the first derivatives of Eq. G.5 we can write
are the FEC of the elastic elements and the explicit dependencies of
and have not been written. The components of the Hessian matrix will be given by
which can be computed from Eqs. G.8 according to
) are the stiffnesses of the elastic components (handles, ssDNA and optical trap) and the variables (
) have been written in terms of the extensions again. All the other components of the Hessian matrix (
) vanish because the first derivative does not depend on the the variable that is being derivated in second order. For instance, because
does not depend on . Note that all second derivatives are symmetric (
), as expected in a function that has continuous second partial derivatives. Now the Hessian matrix can be written as
and its determinant is equal to
which will be necessary to compute the integral. At this point it is important to remember that the components of the Hessian matrix (and so the determinant) depend on
and and the calculation of the stiffnesses (i.e., the second order derivatives) is performed at fixed values of
Now we recall the Taylor expansion written in Eq. G.7 and we introduce it into the calculation of the partition function in Eq. G.6 to obtain the following expression,
where the gradient of
vanishes because we are expanding around the minimum. Now we can perform the integral of the quadratic contribution according to
where we have applied a change of variable
that keeps the same integration limits and allows to perform the Gaussian integration by computing the determinant of the Hessian matrix.
So, in order to calculate the partition function of the whole model we have proceed as follows:
- Fix a value of
will be calculated.
- For this fixed value of
, we have to find the total energy for all the values of . It requires to solve the transcendental equation G.2 for fixed values of
and , find the extensions of all the elastic elements and introduce them into Eq. G.1.
- For each value of (at the same fixed value of
), we have to calculate the Hessian matrix according to Eqs. G.10 and compute its determinant, given by Eq. G.12.
- Once we have all the values of
from steps 2 and 3 respectively for all values of , we have to multiply them and sum over according to Eq. G.13.
- Go back to step 1 to calculate
at another value of
Once the partition function is calculated at all desired values of
, the equation of state is obtained after derivating it according to Eq. 3.27,
Figure G.2 shows the correction introduced in the calculation of the FDC when considering the elasticity of the elements that constitute the whole model. Note that the differences are less than 0.05 pN. Therefore the effort to compute all the second derivatives of the Hessian matrix does not compensate the improvement in the computation of the FDC, which is lower than the experimental resolution of the experiments. However, this effect is important when the length of the handles is longer and they no longer behave like rigid rods.
FDCs for the simplified and full models. Red curve shows the simplified model described in Sec. 3.4.1. Blue curve shows the full model described here in Appendix G.