2 dimensional calculation

Queries about input and output files, running specific calculations, etc.


Moderators: Global Moderator, Moderator

Post Reply
Message
Author
acharyachethan

2 dimensional calculation

#1 Post by acharyachethan » Tue Feb 28, 2006 5:50 am

Hi,

Can I perform ISIF=2 on a 2 dimensional unit cell like graphite. I want to consider only the first layer of graphite. I tried running the program, but I got an input/outut error 137:value not recognized. I also got an error that says The real-space-projection scheme for the treatment of the nonlocal pseudopotentials has been switched on -- but on file POTCAR I have not found any entries signalling me that you have ever performed a real-space optimization of all nonlocal projectors!

I can run 3 dimensional graphite without any problem. I changed the POSCAR file accordingly. Any suggestion would be appreciated!

Thanks!
Last edited by acharyachethan on Tue Feb 28, 2006 5:50 am, edited 1 time in total.

admin
Administrator
Administrator
Posts: 2921
Joined: Tue Aug 03, 2004 8:18 am
License Nr.: 458

2 dimensional calculation

#2 Post by admin » Tue Feb 28, 2006 3:19 pm

1) it is not conflicting to use ISIF=2 for 2D unit cells. However, for your purpose it is more appropriate to use the possibility to choose 'Selective dynamics' in POSCAR and specify the atoms which are allowed to relax: after each ion position a triplet of 'T' or 'F' defines whether this ion is relaxed along the respective x,y, or z co-ordinate.
(please have a look into the online handbook,
http://cms.mpi.univie.ac.at/vasp/vasp/node63.html
on how to use the 'Selective Dynamics' feature)

2) unless turned off explicitely (LREAL=.False.), VASP switches to real space optimisation automatically if the number of atoms in the unit cell exceeds 16. in order to save computing time for large systems
Last edited by admin on Tue Feb 28, 2006 3:19 pm, edited 1 time in total.

acharyachethan

2 dimensional calculation

#3 Post by acharyachethan » Tue Feb 28, 2006 9:48 pm

I can choose Selectice dynamics to fix the coordinates of the bottom carbon atoms in graphite. But, then I want to use a unit cell with only x and y axis, without z axis, with lattice vectors and lattice constants for only x and y axis. Can I do that?
Last edited by acharyachethan on Tue Feb 28, 2006 9:48 pm, edited 1 time in total.

admin
Administrator
Administrator
Posts: 2921
Joined: Tue Aug 03, 2004 8:18 am
License Nr.: 458

2 dimensional calculation

#4 Post by admin » Wed Mar 01, 2006 10:08 am

no, VASP uses 3-dimensional boxes for all unit cells by default
Last edited by admin on Wed Mar 01, 2006 10:08 am, edited 1 time in total.

admin
Administrator
Administrator
Posts: 2921
Joined: Tue Aug 03, 2004 8:18 am
License Nr.: 458

2 dimensional calculation

#5 Post by admin » Wed Mar 01, 2006 10:13 am

no, VASP uses 3-dimensional boxes for all unit cells by default. Surfaces are modelled by 3D repeated slabs (2D-slabs separated by vacuum of sufficient thickness along z)
Last edited by admin on Wed Mar 01, 2006 10:13 am, edited 1 time in total.

Post Reply