Difference between revisions of "DsDNA persistence length"

From OxDNA
(Persistence length of a double-stranded DNA)
m (Persistence length of a double-stranded DNA)
Line 11: Line 11:
 
Note that for calculating a persistence length of a dsDNA, one needs a large number of decorrelated states. To obtain the states (which will be saved into a trajectory file), run the simulation program using the prepared input_persistence file:
 
Note that for calculating a persistence length of a dsDNA, one needs a large number of decorrelated states. To obtain the states (which will be saved into a trajectory file), run the simulation program using the prepared input_persistence file:
  
<tt> oxDNA input_persistence </tt>.
+
<tt>oxDNA input_persistence</tt>
  
 
The program will run a molecular dynamics simulation at 23 °C and record the individual configurations. They are saved in trajectory.dat file. To analyze the data, use the python script dspl.py:
 
The program will run a molecular dynamics simulation at 23 °C and record the individual configurations. They are saved in trajectory.dat file. To analyze the data, use the python script dspl.py:
  
<tt> dspl.py trajectory.dat init.top 10 50 </tt>
+
<tt>dspl.py trajectory.dat init.top 10 50</tt>
  
 
This program will produce a table of correlations between helical vectors, <math> \langle {\bf n_k} \cdot {\bf n_0} \rangle </math>. Using an exponential fit to these data, one can find the persistence length.
 
This program will produce a table of correlations between helical vectors, <math> \langle {\bf n_k} \cdot {\bf n_0} \rangle </math>. Using an exponential fit to these data, one can find the persistence length.

Revision as of 12:31, 17 April 2012

Persistence length of a double-stranded DNA

The example shows how to calculate a persistence length of a double stranded DNA molecule. dsDNA persistence length. The persistence length in this example is calculated using the following formula (see [1] for details):

In the EXAMPLES/PERSISTENCE_LENGTH directory, you will find a setup for calculating the persistence length of a 202 base pairs long dsDNA. Note that for calculating a persistence length of a dsDNA, one needs a large number of decorrelated states. To obtain the states (which will be saved into a trajectory file), run the simulation program using the prepared input_persistence file:

oxDNA input_persistence

The program will run a molecular dynamics simulation at 23 °C and record the individual configurations. They are saved in trajectory.dat file. To analyze the data, use the python script dspl.py:

dspl.py trajectory.dat init.top 10 50

This program will produce a table of correlations between helical vectors, . Using an exponential fit to these data, one can find the persistence length.