VMD-L Mailing List
From: zeynab hosseini (hosseinizeynab93_at_gmail.com)
Date: Sun Jun 13 2021 - 15:06:41 CDT
- Next message: Peter Freddolino: "Re: wraping creates extra bonds!"
- Previous message: zeynab hoseyni: "wraping creates extra bonds!"
- Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]
Hi All,
I'm simulating the bio-molecule (bmol) adsorption on carbon nano-tube (CNT)
solvated in a water box. The CNT is fixed along the z-axis during the
simulation and the height of the water box is the same as the height of the
CNT, in order to avoid entering water molecules into the CNT. Everything is
OK during the simulation (including minimization, equilibration and the NVT
simulation for the data collection). Since sometimes during the simulation
bmol leaves the primary water box (periodic boundary condition is applied),
it is needed to wrap bmol into the primary water box. I used to run the
following commands to wrap bmol into the cell with no problem:
package require pbctools
set nfiles 40
set str ./reduced.psf
mol new $str type psf waitfor all
for {set j 4} {$j <= $nfiles} {incr j 4} {
animate read dcd cnt-asp_red_$j.dcd waitfor all
pbc wrap -centersel "resname CNT" -center com -compound residue -all
set sel [atomselect top "all"]
animate write dcd cnt-asp_red_wrap_$j.dcd sel $sel waitfor all
animate delete all
}
But with this new setup (to consider the same height for the CNT and
simulation box) I see that after wraping the bmol into the cell, some extra
bonds appear on the top and bottom of CNT as it is seen in snapshots
uploaded in the following links:
https://urldefense.com/v3/__https://drive.google.com/file/d/1HsnlXz8bjkhPJkgAcj-nWvwj3WdyDdmm/view?usp=sharing__;!!DZ3fjg!qG6SeEicgqEtVNm1AjfMvi3LHPE9L_-cbGHd77iZSpr6qIdTOdyjCH5dVkro9AlrnA$
https://urldefense.com/v3/__https://drive.google.com/file/d/1zVU_o6GxR7zI6f4Ps-yIzghdMRJltRtu/view?usp=sharing__;!!DZ3fjg!qG6SeEicgqEtVNm1AjfMvi3LHPE9L_-cbGHd77iZSpr6qIdTOdyjCH5dVkpybl3Kqw$
https://urldefense.com/v3/__https://drive.google.com/file/d/19p2ykduU3RqBixjo2C4hkMknoSn9yp18/view?usp=sharing__;!!DZ3fjg!qG6SeEicgqEtVNm1AjfMvi3LHPE9L_-cbGHd77iZSpr6qIdTOdyjCH5dVkqOWG2edA$
I followed several ways to get rid of those bonds but still unsuccessful.
Would be thankful if anyone guides me.
All the best,
Zeynab
- Next message: Peter Freddolino: "Re: wraping creates extra bonds!"
- Previous message: zeynab hoseyni: "wraping creates extra bonds!"
- Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]