| ²é¿´: 630 | »Ø¸´: 7 | |||
| µ±Ç°Ö÷ÌâÒѾ´æµµ¡£ | |||
xli6гæ (³õÈëÎÄ̳)
|
[½»Á÷]
¡¾ÇóÖú¡¿relax
|
||
|
Hi, I am the new guy for siesta, my current problem is how can I relax my model to get my bond length or structure? Hopefully, some one can help me or give some hint about command I need to use in siesta input file.Thanks a lot!Sample Text |
» ²ÂÄãϲ»¶
Çóµ÷¼Á
ÒѾÓÐ23È˻ظ´
301Çóµ÷¼Á
ÒѾÓÐ15È˻ظ´
304Çóµ÷¼Á£¨085602£¬¹ýËļ¶£¬Ò»Ö¾Ô¸985£©
ÒѾÓÐ17È˻ظ´
302·ÖÇóµ÷¼Á Ò»Ö¾Ô¸°²»Õ´óѧ085601
ÒѾÓÐ12È˻ظ´
288»·¾³×¨Ë¶,Çóµ÷²ÄÁÏ·½Ïò
ÒѾÓÐ23È˻ظ´
»·¾³×¨Ë¶µ÷¼Á
ÒѾÓÐ6È˻ظ´
22408 µ÷¼Á²ÄÁÏ
ÒѾÓÐ6È˻ظ´
285Çóµ÷¼Á
ÒѾÓÐ12È˻ظ´
Çóµ÷¼Á
ÒѾÓÐ6È˻ظ´
085600²ÄÁÏÓ뻯¹¤301·ÖÇóµ÷¼ÁԺУ
ÒѾÓÐ19È˻ظ´
xhsh
ľ³æ (СÓÐÃûÆø)
- 1STÇ¿Ìû: 1
- Ó¦Öú: 2 (Ó×¶ùÔ°)
- ½ð±Ò: 2389.2
- ºì»¨: 3
- Ìû×Ó: 287
- ÔÚÏß: 58.5Сʱ
- ³æºÅ: 501440
- ×¢²á: 2008-02-14
- ÐÔ±ð: GG
- רҵ: Äý¾Û̬ÎïÐÔ II £ºµç×ӽṹ
¡ï ¡ï
Сľ³æ(½ð±Ò+0.5):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû½»Á÷
zxzj05(½ð±Ò+1):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû½»Á÷ 1-23 14:45
Сľ³æ(½ð±Ò+0.5):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû½»Á÷
zxzj05(½ð±Ò+1):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû½»Á÷ 1-23 14:45
|
MD.TypeOfRun cg MD.NumCGsteps 80 MD.MaxCGDispl 0.1 Ang MD.MaxForceTol 0.04 eV/Ang You'd better carefully read the manual. This is the best way to learn siesta. |
2Â¥2010-01-22 08:21:26
xli6
гæ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 81.4
- Ìû×Ó: 17
- ÔÚÏß: 20.2Сʱ
- ³æºÅ: 941695
- ×¢²á: 2010-01-13
3Â¥2010-01-22 08:27:11
xli6
гæ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 81.4
- Ìû×Ó: 17
- ÔÚÏß: 20.2Сʱ
- ³æºÅ: 941695
- ×¢²á: 2010-01-13
| why I use this kind of command, it does not work? it look no calculation and I only have one unit cell(4 atoms). do I need any other options? ElectronicTemperature? usually what is the most important option for optimization? sorry because I can not read the manual at one time, so I ask too much! |
4Â¥2010-01-22 11:11:42
xhsh
ľ³æ (СÓÐÃûÆø)
- 1STÇ¿Ìû: 1
- Ó¦Öú: 2 (Ó×¶ùÔ°)
- ½ð±Ò: 2389.2
- ºì»¨: 3
- Ìû×Ó: 287
- ÔÚÏß: 58.5Сʱ
- ³æºÅ: 501440
- ×¢²á: 2008-02-14
- ÐÔ±ð: GG
- רҵ: Äý¾Û̬ÎïÐÔ II £ºµç×ӽṹ
¡ï
Сľ³æ(½ð±Ò+0.5):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû½»Á÷
Сľ³æ(½ð±Ò+0.5):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû½»Á÷
| You may post your input file here and let me check. |
5Â¥2010-01-22 16:43:24
fanchen021
Ìú¸Ëľ³æ (ÕýʽдÊÖ)
- Ó¦Öú: 7 (Ó×¶ùÔ°)
- ½ð±Ò: 7994.9
- É¢½ð: 1844
- ºì»¨: 4
- Ìû×Ó: 776
- ÔÚÏß: 569Сʱ
- ³æºÅ: 773188
- ×¢²á: 2009-05-18
- רҵ: ÀíÂۺͼÆË㻯ѧ
¡ï
Сľ³æ(½ð±Ò+0.5):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû½»Á÷
Сľ³æ(½ð±Ò+0.5):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû½»Á÷
| There are some examples for geometry optimization in SIESTA |
6Â¥2010-01-22 17:56:27
xli6
гæ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 81.4
- Ìû×Ó: 17
- ÔÚÏß: 20.2Сʱ
- ³æºÅ: 941695
- ×¢²á: 2010-01-13
|
two more questions, I appreciate someone who can help me. 1, In siesta, is boundary condition always assumed periodic? 2, when I am doing optimisation for my structure, do I need to care how much atoms I have since everything is periodic? Is there any difference between my model is one unit cell(4 atoms) or a couple of unit cells? Below are some options I used for optimisation, but it looks the whole structure is totally changed. I do not know why. Thanks a lot! MD.TypeOfRun CG MD.NumCGsteps 500 MD.MaxCGDispl 0.2 Bohr MD.MaxForceTol 0.04 eV/Ang SolutionMethod Diagon MeshCutoff 200 Ry |
7Â¥2010-01-23 07:10:35
xhsh
ľ³æ (СÓÐÃûÆø)
- 1STÇ¿Ìû: 1
- Ó¦Öú: 2 (Ó×¶ùÔ°)
- ½ð±Ò: 2389.2
- ºì»¨: 3
- Ìû×Ó: 287
- ÔÚÏß: 58.5Сʱ
- ³æºÅ: 501440
- ×¢²á: 2008-02-14
- ÐÔ±ð: GG
- רҵ: Äý¾Û̬ÎïÐÔ II £ºµç×ӽṹ
¡ï ¡ï
Сľ³æ(½ð±Ò+0.5):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû½»Á÷
zxzj05(½ð±Ò+1):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû½»Á÷ 1-23 14:46
Сľ³æ(½ð±Ò+0.5):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû½»Á÷
zxzj05(½ð±Ò+1):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû½»Á÷ 1-23 14:46
|
1. yes. 2. check your LatticeConstat and LatticeVectors. |
8Â¥2010-01-23 08:31:05














»Ø¸´´ËÂ¥