From: kanchi subbarao rao (ksubbu85_at_gmail.com)
Date: Tue May 31 2011 - 05:10:32 CDT
---------- Forwarded message ----------
From: kanchi subbarao rao <ksubbu85_at_gmail.com>
Date: Tue, May 31, 2011 at 11:24 AM
Subject: Fwd: hi
To: jordi_at_ks.uiuc.edu
---------- Forwarded message ----------
From: kanchi subbarao rao <ksubbu85_at_gmail.com>
Date: Tue, May 31, 2011 at 11:23 AM
Subject: hi
To: jordi_at_ks.uiuc.edu
hi jordi,
i am doing nanotubes tutorial, in that to build the nano tubes
there is given a make_arm.tcl file.
that is given below
# Written by Fangqiang Zhu, 6/3/2002
# Make a m*m nanotube with length n.
# Write to $output.psf and $output.pdb
proc make_arm {m n output} {
# Bond length
set a 1.3750
set PI 3.141592654
set tt [expr 2*$PI/(3*$m)]
set r [expr $a/(2*sin($PI/3/$m))]
set fd_top [open temp_armchair.top w]
# Version number
puts $fd_top "0 0"
puts $fd_top "MASS 1 CA 12.01100 C"
puts $fd_top "AUTO ANGLES DIHE"
puts $fd_top "RESI ARM 0.00"
set fd_pgn [open temp_armchair.pgn w]
puts $fd_pgn "topology temp_armchair.top"
puts $fd_pgn "segment NT {residue 1 ARM}"
# Define the atom names for each carbon
set ind 0
for {set i 0} {$i<$n} {incr i} {
for {set j 0} {$j<$m} {incr j} {
for {set k 0} {$k<2} {incr k} {
if {$ind > 99} {
set name($i,$j,$k) C$ind
} elseif {$ind > 9} {
set name($i,$j,$k) C0$ind
} else {
set name($i,$j,$k) C00$ind
}
incr ind
}
}
}
for {set i 0} {$i<$n} {incr i} {
set z [expr ($i-($n-1)/2.0)*sqrt(3)/2*$a]
set t0 [expr ($i%2==0)?0:1.5*$tt]
for {set j 0} {$j<$m} {incr j} {
puts $fd_top "ATOM $name($i,$j,0) CA 0.00"
set t [expr $t0+(3*$j-0.5)*$tt]
puts $fd_pgn "coord NT 1 $name($i,$j,0) {[expr $r*cos($t)] [expr
$r*sin($t)] $z}"
puts $fd_top "ATOM $name($i,$j,1) CA 0.00"
set t [expr $t0+(3*$j+0.5)*$tt]
puts $fd_pgn "coord NT 1 $name($i,$j,1) {[expr $r*cos($t)] [expr
$r*sin($t)] $z}"
puts $fd_top "BOND $name($i,$j,0) $name($i,$j,1)"
if {$i%2 == 1} {
set ii [expr $i-1]
set jj [expr ($j+1)%$m]
puts $fd_top "BOND $name($ii,$j,1) $name($i,$j,0)"
puts $fd_top "BOND $name($ii,$jj,0) $name($i,$j,1)"
} elseif {$i > 0} {
set ii [expr $i-1]
set jj [expr ($j-1)%$m]
puts $fd_top "BOND $name($ii,$jj,1) $name($i,$j,0)"
puts $fd_top "BOND $name($ii,$j,0) $name($i,$j,1)"
}
}
}
puts $fd_pgn "writepsf $output.psf"
puts $fd_pgn "writepdb $output.pdb"
close $fd_top
close $fd_pgn
exec psfgen temp_armchair.pgn
exec rm temp_armchair.top
exec rm temp_armchair.pgn
}
while building nano tube in tcl-console by the commands
*make arm 6 12 nanotube*
*
*
i am getting the error *couldn't execute "psfgen": no such file or
directory *
i can not debug this error shall you helpme?
thanking you
subbarao kanchi.
This archive was generated by hypermail 2.1.6 : Mon Dec 31 2012 - 23:20:21 CST