Endonuclease PvuII (1PVI) DNA - GATTACAGATTACA
CAP - Catabolite gene Activating Protein (1BER)
DNA - GATTACAGATTACAGATTACA Endonuclease PvuII bound to palindromic DNA recognition site CAGCTG (1PVI) DNA - GATTACAGATTACAGATTACA TBP - TATA box Binding Protein (1C9B)
CAP - Catabolite gene Activating Protein (1BER)
GCN4 - leucine zipper transcription factor bound to palindromic DNA recognition site ATGAC(G)TCAT (1YSA)
GCN4 - leucine zipper transcription factor bound to palindromic DNA recognition site ATGAC(G)TCAT (1YSA)
GCN4 - leucine zipper transcription factor bound to palindromic DNA recognition site ATGAC(G)TCAT (1YSA)
GCN4 - leucine zipper transcription factor bound to palindromic DNA recognition site ATGAC(G)TCAT (1YSA)
GCN4 - leucine zipper transcription factor bound to palindromic DNA recognition site ATGAC(G)TCAT (1YSA)
TBP - TATA box Binding Protein (1C9B)
 

° 

PressureCtrl

-

Set pressure control


CommandArgument DatatypeDefaultMinMax
Format 1:PressureCtrl SolventProbe,STRING - --
  Name = Solvent residue name, STRING---
  Density = Solvent density in g/ml,FLOAT 1.0--
   Axis = X | Y | Z | XY | XZ | YZ | XYZ STRINGXYZ- -
Format 2:PressureCtrlOff STRING-- -
Python:PressureCtrl(noname1,name,density=None,axis=None)
PressureCtrlOff()
Menu:Simulation > Pressure control
Related: TempCtrl , SolvDensity, FillCellWater , FillCellObj
Required:


The PressureCtrl command sets the pressure control mode.

Conformational changes of the simulated biomolecules lead to fluctuations in density. If the simulation box has a constant size, changes in density lead to changes in pressure. This is somewhat unrealistic, because proteins normally fold and function in a constant pressure environment.

If PressureCtrl is activated, YASARA uses the time averaged solvent density as a probe for the pressure in the box. The box is then rescaled to arrive at the target solvent density and thus pressure. Simulations are normally run at a pressure of 1 atm and 298 K, where the experimentally determined water density is 0.997 g/ml. The densities of other solvents are listed at the end.

Instructions on how to fill the cell with a certain solvent are provided at the FillCellWater and FillCellObj commands.

Since the solvent density depends on the temperature, it has to be adjusted for simulations run at temperatures other than 298K. As an example, the water density as a function of temperature is listed below:

Temperature [Celsius]Density [g/ml]
0 0.99984
10 0.99970
20 0.99821
30 0.99565
40 0.99222
50 0.98803
60 0.98320
70 0.97778
80 0.97182
90 0.96535
100 0.95840

Rescaling does not only involve the simulation cell, but includes also the atom coordinates. An exception are simulations in non-periodic cells with atoms outside the cell. In this case, only the cell is rescaled. Note that pressure control in non-periodic cells is of limited use : if the cell size increases, the surface tension causes the water molecules to detach from the cell wall and form a liquid droplet, that is not influenced by pressure control and cell scaling.

The Axis parameter defines along which axes the cell is rescaled to reach the requested water density. It is sometimes helpful to keep the length of certain cell axes fixed, for example when simulating objects like membranes that wrap around the periodic boundaries. Note however that the solvent density is not a reliable pressure indicator when there are no solvent molecules left in a certain plane through the cell.

Densities of various solvents at 298K in g/ml:

Acetic acid1.0429
Acetone 0.7856
Acetonitrile 0.7793
Benzene0.8729
Bromoform2.8761
Butyl acetate 0.8761
Butyl alcohol 0.8061
Butylamine0.7406
Carbon disulfide1.2556
Carbon tetrachloride 1.5833
Chlorobenzene 1.1007
Chloroform1.4800
Cumene0.8597
Cyclohexane 0.7731
Cyclohexanol 0.9604
Cyclohexanone0.9425
Cyclohexylamine0.8627
Diacetone alcohol 0.9342
Dibutylamine 0.7571
Dichloroethyl ether1.2130
Dichloromethane1.3182
Diethanolamine 1.0899
Diethylamine 0.7017
Diethylene glycol1.1150
Diethyl ether0.7080
Dimethylamine 0.6501
Dimethyl sulfoxide 1.0955
Ethanol0.7873
Ethanolamine1.0136
Ethyl acetate 0.8942
Ethylamine 0.6766
Ethylbenzene0.8654
Ethyl chloride0.8902
Ethylene glycol 1.1101
Ethyl formate 0.9165
Furan0.9348
Furfural1.1554
Furfuryl alcohol 1.1308
Glycerol1.2567
Heptane0.6816
Hexane 0.6563
Hexylene glycol 0.9182
Isobutyl acetate0.8695
Isobutyl alcohol0.7967
Isobutylamine0.7296
Isopentyl acetate 0.8666
Isopropyl acetate 0.8711
Isopropyl alcohol0.7827
Methanol0.7872
Methyl acetate 0.9273
Methylamine 0.6552
Morpholine0.9959
Nitrobenzene1.1985
Nitroethane 1.0427
Nitromethane 1.1286
Octane0.6986
Pentachloroethane1.6749
Pentane 0.6215
Pentyl acetate 0.8721
Piperidine0.8578
Propyl alcohol0.8020
Propyl acetate 0.8826
Propylamine 0.7139
Propylene glycol1.0327
Pyridine0.9786
Sulfolane 1.2660
Tetrachloroethylene 1.6130
Tetrahydrofuran0.8800
Toluene0.8647
Tributylamine 0.7748
Trichloroethylene 1.4578
Water0.997

Example 1:
PressureCtrl Off

Disable pressure control.


Example 2:
PressureCtrl SolventProbe,Name=HOH,Density=0.997

Measure water density and slowly resize simulation cell to reach 0.997 g/ml.


Example 3:
PressureCtrl SolventProbe,Name=DMS,Density=1.096,Axis=Y

Measure density of the dimethyl sulfoxide solvent and resize cell along the Y axis to reach the density 1.096 g/ml.