My Community

I know that in order to calculate vacuum potential of a system, one should turn on LVTOT tag and then get a LOCPOT file as result. I also know that you can use a script like vtotav.f to read the LOCPOT file and find the vacuum potential for surfaces. But I have no idea what I should do to get vacuum potential of carbon nanotube system from LOCPOT file. vtotav.f seems only work for surfaces, not carbon nanotubes with 2-dimention vacuum spaces. Any suggestions on this?
Thank you.

of course vtotav works for 1D periodc structures as well.

Thanks. Then, how do I answer the question'Which direction to keep? (1-3 --- 1=X,2=Y,3=Z)'? Suppose the periodic condition is applied in z direction.

then either give 1 or 2, one of the 2 directions perpendicular to the nanotube.