| ²é¿´: 3362 | »Ø¸´: 17 | |||
| µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû | |||
[½»Á÷]
¡¾ÇóÖú¡¿Çó²úÉú±äÐκóPOSCARµÄ defvector.f ³ÌÐòµÄÔ´Îļþ
|
|||
| ×î½ü¸ÕѧϰʹÓÃVASPÈí¼þ£¬¿´ºî²©µÄ˵Ã÷ÊéÖÐÌáµ½¼ÆË㵯ÐÔ³£ÊýÒªÓõ½defvector.f£¬ÇëÎÊÕâ¸ö³ÌÐòÔÚÄÄÀï¿ÉÒÔÕÒµ½ÄØ£¿Íû´ó¼Ò¸øÖ¸µ¼Ö¸µ¼£¬²»Ê¤¸Ð¼¤°¡ |
» ²ÂÄãϲ»¶
¿¼Ñл¯Ñ§Ñ§Ë¶µ÷¼Á£¬Ò»Ö¾Ô¸985
ÒѾÓÐ5È˻ظ´
286Çóµ÷¼Á
ÒѾÓÐ4È˻ظ´
328Çóµ÷¼Á£¬Ó¢ÓïÁù¼¶551£¬ÓпÆÑоÀú
ÒѾÓÐ4È˻ظ´
»úеר˶325£¬Ñ°ÕÒµ÷¼ÁԺУ
ÒѾÓÐ5È˻ظ´
²ÄÁÏר˶306Ó¢Ò»Êý¶þ
ÒѾÓÐ6È˻ظ´
»ï°éÃÇ£¬×£ÎÒÉúÈÕ¿ìÀÖ°É
ÒѾÓÐ26È˻ظ´
333Çóµ÷¼Á
ÒѾÓÐ7È˻ظ´
0854¿ØÖƹ¤³Ì 359Çóµ÷¼Á ¿É¿çרҵ
ÒѾÓÐ9È˻ظ´
Áº³ÉΰÀÏʦ¿ÎÌâ×é»¶ÓÄãµÄ¼ÓÈë
ÒѾÓÐ9È˻ظ´
»¯Ñ§µ÷¼Á0703
ÒѾÓÐ8È˻ظ´
» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:
ÇóÖúѹËõÓÀ¾Ã±äÐÎÎÊÌâ
ÒѾÓÐ6È˻ظ´
¾Û¼×»ù±ûÏ¡ËáôÇÒÒõ¥£¨HEMA£©¾ÛºÏ¼Ó¹¤³ÉÐͱäÐÎ
ÒѾÓÐ4È˻ظ´
²ÄÁϲúÉúËÜÐÔ±äÐÎÊDz»ÊÇʹÄÚÓ¦Á¦µÃµ½ÁËÊÍ·Å
ÒѾÓÐ14È˻ظ´
PA6¡¢PETµÄÈȱäÐÎζȷֱðÊǶàÉÙ£¿
ÒѾÓÐ7È˻ظ´
¸÷λÑо¿Ñ¹µçÌմɵĸßÊÖÃÇ ÇóÖú¹ØÓÚ·½ÐÎѹµçƬµÄ±äÐÎÐÎʽ
ÒѾÓÐ3È˻ظ´
Îļþ.fÈçºÎ±àÒë³É¿ÉÖ´ÐÐÎļþ
ÒѾÓÐ3È˻ظ´
¡¾ÇóÖú¡¿×ö±äÐνð¸Õ¹Ø½ÚÓ¦¸ÃÓÃÄÄÖÖÏð½º
ÒѾÓÐ8È˻ظ´
¡¾ÇóÖú¡¿¼Ó¹¤Ó²»¯ÂÊ»òÕßÓ¦±äÓëÀäÔþµÄ±äÐÎÁ¿ÓÐûÓйØÏµÄØ£¿
ÒѾÓÐ9È˻ظ´
¡¾ÇóÖú¡¿ÇëÎʽðÊô¡¢Ìմɺ͸߷Ö×ӵıäÐλúÖÆ·Ö±ðÊÇʲô£¿
ÒѾÓÐ5È˻ظ´
» ÇÀ½ð±ÒÀ²£¡»ØÌû¾Í¿ÉÒԵõ½:
Î÷ºþ´óѧ¹¤Ñ§ÔºÖÇÄÜ΢ÐÍ»úÆ÷ÈËʵÑéÊÒ²©Ê¿ºóÕÐÆ¸ÆôÊÂ
+4/236
µçÁ¦È«¹úÖØµãʵÑéÊÒ˫һÁ÷AÀ೤½Ñ§ÕßÍŶÓÕÐ2026ÄêÈ«ÈÕÖÆ²©Ê¿1Ãû£¨4ÔÂ3ÈÕ½ØÖ¹£©
+2/168
ÄϾ©Ò½¿Æ´óѧÉúֳҽѧÓë×Ó´ú½¡¿µ¹úÖØÊµÑéÊÒ-¶·Éºê¿ÎÌâ×é-ѧÊõ˶ʿÕÐÉú
+1/104
×ø±ê¹ãÖÝ£¬ÄÐÕ÷Å®£¬Ë«Ïò±¼¸°
+1/83
Call for paper: Frontiers in Plant Science ÖпÆÔº2Çø
+2/40
ÕÒ¹¤×÷¾ÑéÇóÖú
+1/36
ҽѧ327·Öר˶µ÷¼Á
+1/16
Çൺ¿Æ¼¼´óѧ0860 µ÷¼ÁÕÐÉú
+1/13
ʳƷ¿ÆÑ§Ó빤³Ì083200-297·Ö
+1/8
¸£ÖÝ´óѧ¹ú¼Ò¼¶È˲ÅÌÀÓýÐÀ½ÌÊÚ¿ÎÌâ×é2026ÄêÕÐÊÕ¡±ÉêÇ뿼ºËÖÆ¡°´¢ÄÜµç³Ø·½Ïò²©Ê¿Ñо¿Éú
+1/8
Çൺ¿Æ¼¼´óѧ¿É³ÖÐø¸ß·Ö×ÓÍÅ¶Ó ¿¼ÑÐÕÐÉú
+1/7
ÖйúÃñº½·ÉÐÐѧԺÏȽøÁ¬½ÓÍŶÓÕÐÊÕº½¿ÕÓ¡¢º½¿Õ»úеרҵ˶ʿÑо¿Éú£¨Ñ§Ë¶×¨Ë¶¾ù¿É£©
+1/7
26Ä격ʿÕÐÉú
+1/6
ÇëÎÊΪʲôÎÒµÄÌù×Ó×Ô¶¯µ½¡°ÒÑɾ³ýµÄÖ÷Ì⡱Àï¶øÃ»Óб»·¢²¼£¿
+1/5
·¢paper
+1/5
Ïã¸Û¿Æ¼¼´óѧ£¨¹ãÖÝ£© ÕÐÊÕÎÂÊÒÆøÌå¼à²âÓë¼õÅÅ·½ÏòµÄ²©Ê¿¡¢²©ºó¡¢Ë¶Ê¿¡¢RA
+1/4
Á¬·¢Á½Æª Nature£¡´óÑù±¾µ¥Ï¸°û²âÐò½âÎöÈÏÖª¡¢Ë¥ÀÏÓëÁ鳤Àà½ø»¯
+1/3
085404 270Çóµ÷¼Á£¬bÇøÇóµ¼Ê¦ÊÕÁô½ÓÊÜ¿çרҵ
+1/2
Î÷½»ÀûÎïÆÖ´óѧ/µª»¯ïØ»ù CMOS ¼¼ÊõµÄÎïÀíÇý¶¯ÓëÈ˹¤ÖÇÄÜÔöÇ¿·½·¨/Õв©Ê¿Ñо¿Éú
+1/2
284 Çóµ÷¼Á
+1/1
·çÓêͬÖÜ113(½ð±Ò+16): 2011-07-16 21:39:37
|
C%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% C >this simple program to get the primitive vectors after C $\delta$ strain, in order to calculate the independent C elastic constants of solids. C usage: C!!!!! Please first prepare the undeformed POSCAR in OLDPOS C >defvector.x C >type defvector.x > create new POSCAR in file fort.3 C%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% program defvector real*8 privect,strvect,delta,strten,strain,pos, alat dimension privect(3,3),strvect(3,3),strten(3,3),strain(6) dimension pos(50,3) character*10 bravlat, title, direct integer i,j,k,ntype, natomi, nn dimension natomi(10) C%%%%%%%%% Read the undeformed primitive vector and atomic postion %%%%%%% open(7,file='OLDPOS') C%% In first line of OLDPOS, please add the number C%% of the type of atoms after the title read(7,*) title, ntype read(7,*) alat do i=1,3 read(7,*) (privect(i,j),j=1,3) write(*,*) (privect(i,j),j=1,3) enddo read(7,*) (natomi(i),i=1,ntype) nn=0 do i =1, ntype nn=nn+natomi(i) enddo read(7,*) direct do i=1, nn read(7,*) (pos(i,j),j=1,3) enddo C%%%%%%%%% Read the amti of strain %%%%%%%%%%%%%%% read(*,*) delta C%%%%%%%%% Define the strain %%%%%%%%%%%%%% strain(1)=0.0 strain(2)=delta strain(3)=delta strain(4)=0.0 strain(5)=0.0 strain(6)=0.0 C%%%%%%%%% Define the strain tensor %%%%%%%%%%%%%%%%%%%%%%%% strten(1,1)=strain(1)+1.0 strten(1,2)=0.5*strain(6) strten(1,3)=0.5*strain(5) strten(2,1)=0.5*strain(6) strten(2,2)=strain(2)+1.0 strten(2,3)=0.5*strain(4) strten(3,1)=0.5*strain(5) strten(3,2)=0.5*strain(4) strten(3,3)=strain(3)+1.0 C%%%%%%%%% Transform the primitive vector to the new vector under strain%%%%% C strvect(i,j)=privect(i,j)*(I+strten(i,j)) do k=1,3 do i=1,3 strvect(i,k)=0.0 do j=1,3 strvect(i,k)=strvect(i,k)+privect(i,j)*strten(j,k) enddo enddo enddo C%%%%%%%% Write the new vector under strain%%%%%%%%%%%% do i=1,3 write(*,100)(strvect(i,j),j=1,3) enddo 100 format(3f20.15) C%%%%%%%%% Create the POSCAR for total energy calculation %%%%%%%%%%%%%%5 write(3,'(A10)') title write(3,'(f15.10)') alat do i=1,3 write(3,100)(strvect(i,j),j=1,3) enddo write(3,'(10I4)') (natomi(i), i=1,ntype) write(3,'(A6)') Direct do i=1, nn write(3,100) (pos(i,j),j=1,3) enddo C%%%%%%% end |
4Â¥2011-04-12 10:27:42
¡ï
sunyang1988(½ð±Ò+1): ллÌáʾ 2011-04-10 00:04:45
·çÓêͬÖÜ113(½ð±Ò+2): 2011-04-10 12:53:42
sunyang1988(½ð±Ò+1): ллÌáʾ 2011-04-10 00:04:45
·çÓêͬÖÜ113(½ð±Ò+2): 2011-04-10 12:53:42
|
´ò¿ªÖ¸ÄÏ,¼ô/ÌùÒ»ÏÂdefvector.f,±àÒë. »ò Ö±½Ó¸øºî²©Ê¿·¢ÓʼþË÷È¡. |
2Â¥2011-04-09 23:40:39
ufo991
Ìú³æ (ÕýʽдÊÖ)
- 1STÇ¿Ìû: 1
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 1148.3
- Ìû×Ó: 358
- ÔÚÏß: 104.6Сʱ
- ³æºÅ: 671678
¡ï
Сľ³æ: ½ð±Ò+0.5, ¸ø¸öºì°ü£¬Ð»Ð»»ØÌû
Сľ³æ: ½ð±Ò+0.5, ¸ø¸öºì°ü£¬Ð»Ð»»ØÌû
|
ΪºÎÎÒ±àÒëÄÜͨ¹ý£¬ÔËÐеÄʱºò³öÏÖinvalid number:incomprehensible list input apparent state: unit 7 named OLDPOS last format: list io lately reading direct formatted external IO »òÕßapparent state: unit 5 (unnamed) last format: list io lately reading direct formatted external IO list in: end of file µÄÎÊÌâÄØ£¿ |
5Â¥2012-05-30 10:31:31
6Â¥2012-06-07 06:00:01













»Ø¸´´ËÂ¥