Znn3bq.jpeg
ÉÇÍ·´óѧº£Ñó¿ÆÑ§½ÓÊܵ÷¼Á
²é¿´: 2573  |  »Ø¸´: 5
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

gtssongchi

ľ³æ (СÓÐÃûÆø)

[ÇóÖú] fix nvt ζÈÉèÖÃΪ300K£¬¿ÉÊǽá¹ûÈ´ÏÔʾֻÓÐ150K£¬ÊÇʲôԭÒòÄØ£¿ ÒÑÓÐ2È˲ÎÓë

´ó¼ÒºÃ£¬lammps³õѧ£¬ÇëÎÊÎÒÔÚÉèÖÃÖгöÁËʲôÎÊÌ⣬µ¼ÖÂ×îºóÈÈÁ¦Ñ§Êä³öµÄζÈÖ»ÓÐ150K
#lammps in file for melt of Cu

echo                        screen

units                        metal
boundary                p p p
atom_style                atomic

lattice                        fcc        3.61        #½¨Á¢Ò»¸ö¾§¸ñ²ÎÊýΪ3.61°£µÄfcc¾§Ìå
region                        box        block        0 10 0 10 0 10       
create_box                1        box       
create_atoms                1        box       

pair_style                eam
pair_coeff                * * ./potentials/Cu_u3.eam        #¸ñʽȷ±£ÕýÈ·


timestep                0.003        #µ¥Î»ÊÇps
neighbor                2.0        bin
neigh_modify                delay 0 every 5 check yes

velocity                all        create        300  474747 dist gaussian

#variable                N        equal        abs{floor(x[100]+10000)}
#velocity                all        create        N        dist gaussian

#fix                        1        all  nvt  temp  300  300  100.0        # 100ÊÇñîºÏϵÊý£¬¸ø³ö³Ûԥʱ¼ä100¡Áʱ¼äµ¥Î»
fix                        1         all        nve
fix                        2        all  temp/berendsen 300.0 300.0 100.0


thermo                        100       
#thermo_style                ĬÈÏ
#dump                        1        all        custom  100  dump.lammpstrj  id  type  x  y  z       

run                        10000
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

brucefan

ר¼Ò¹ËÎÊ (ÖøÃûдÊÖ)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
gtssongchi: ½ð±Ò+15, ¡ï¡ï¡ï¡ï¡ï×î¼Ñ´ð°¸ 2017-04-20 21:56:31
ÒýÓûØÌû:
2Â¥: Originally posted by cl1994861 at 2017-04-20 16:23:07
³Úԥʱ¼äÊDz»ÊÇÉèÖõIJ»ºÏÀí£¬¹ÙÍømannualÀï¸øµÄÒ»°ãÊÇ1000¸öʱ¼ä²½

No, 100 here is too large in the "metal" unit system, where the unit of time is ps. Change it to 0.1.
Also, 10000 steps are too short to reach equilibrium.
4Â¥2017-04-20 17:22:06
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 6 ¸ö»Ø´ð

cl1994861

½ð³æ (СÓÐÃûÆø)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
gtssongchi: ½ð±Ò+5, ¡ï¡ï¡ïºÜÓаïÖú 2017-04-20 21:56:48
³Úԥʱ¼äÊDz»ÊÇÉèÖõIJ»ºÏÀí£¬¹ÙÍømannualÀï¸øµÄÒ»°ãÊÇ1000¸öʱ¼ä²½
2Â¥2017-04-20 16:23:07
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

gtssongchi

ľ³æ (СÓÐÃûÆø)

ÒýÓûØÌû:
2Â¥: Originally posted by cl1994861 at 2017-04-20 16:23:07
³Úԥʱ¼äÊDz»ÊÇÉèÖõIJ»ºÏÀí£¬¹ÙÍømannualÀï¸øµÄÒ»°ãÊÇ1000¸öʱ¼ä²½

ÎÒ·¢ÏÖÎÊÌâÔ­Òò¿ÉÄÜÊÇÅܵÃʱ¼ä²½²»¹»£¬ÕâÀïÖ®ÉèÖÃÁË10000²½£¬½«Ê±¼ä²½¸ÄΪ50000²½ºó£¬ºóÆÚζÈȷʵÎȶ¨ÔÚ300KÁË£¬Ö»ÊÇÎ񼆮æ¹Ö£¬ÎªÊ²Ã´lammps¸Õ¿ªÊ¼µÄζȻáÊÇ150KÄØ£¬Ëü´Ó150KÂýÂýÉÏÉýµ½300K£¬È»ºóÎȶ¨¡£
3Â¥2017-04-20 17:19:03
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

gtssongchi

ľ³æ (СÓÐÃûÆø)

ÒýÓûØÌû:
4Â¥: Originally posted by brucefan at 2017-04-20 17:22:06
No, 100 here is too large in the "metal" unit system, where the unit of time is ps. Change it to 0.1.
Also, 10000 steps are too short to reach equilibrium....

¸Ðл£¬¿ÉÄÜÊÇÎÒÀí½âÓÐÎ󣬻»³É0.1Ö®ºó£¬µÄÈ·½â¾öÁËÕâ¸öÎÊÌâ¡£
5Â¥2017-04-20 21:12:20
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] ¸´ÊÔµ÷¼Á +17 »ý¼«ÏòÉÏ£» 2026-04-10 19/950 2026-04-13 19:53 by ÏÐÈËÖÕÄÏɽ
[¿¼ÑÐ] ʳƷÓëÓªÑø£¨0955£©271Çóµ÷¼Á +9 Éý¸ñ°¢´ï 2026-04-12 9/450 2026-04-13 14:12 by maddjdld
[¿¼ÑÐ] Ò»Ö¾Ô¸ÏôóÉúÎïѧ332Çóµ÷¼Á +11 ³Ø³Ø³Ø³Ø³Ø³Ø 2026-04-08 11/550 2026-04-13 14:10 by ¿ÆÑÐÂÛ
[¿¼ÑÐ] 271Çóµ÷¼Á +23 2261744733 2026-04-11 23/1150 2026-04-13 08:48 by ѧzh
[¿¼ÑÐ] µç×ÓÐÅÏ¢270Çóµ÷¼Á +18 terminal469 2026-04-07 18/900 2026-04-12 16:23 by ajpv·çÀ×
[¿¼ÑÐ] 307Çóµ÷¼Á +10 tzq94092 2026-04-10 10/500 2026-04-12 08:18 by wise999
[¿¼ÑÐ] 280Çóµ÷¼Á +13 wzzzÍõ 2026-04-09 13/650 2026-04-12 00:31 by ÓÂÅʸ߷å0126
[¿¼ÑÐ] 296Çóµ÷¼Á +8 Íô£¡£¿£¡ 2026-04-09 8/400 2026-04-11 21:02 by ÄæË®³Ë·ç
[¿¼ÑÐ] 086003µ÷¼ÁÇóÖú +21 ËÕß®Íò 2026-04-09 22/1100 2026-04-11 20:25 by dongdian1
[¿¼ÑÐ] 283Çóµ÷¼Á +22 ÄǸöàà×Ó 2026-04-09 22/1100 2026-04-11 10:41 by ÄæË®³Ë·ç
[¿¼ÑÐ] 342µç×ÓÐÅϢר˶Çóµ÷¼Á +9 ÄãÈÃÎÒÔõôÀóÖ¦ 2026-04-10 10/500 2026-04-11 08:33 by zhq0425
[¿¼ÑÐ] Ò»Ö¾Ô¸»ªÄÏÀí¹¤´óѧ331·Ö²ÄÁÏÇóµ÷¼Á +9 ÌìÏÂww 2026-04-09 9/450 2026-04-10 22:58 by Ftglcn90
[¿¼ÑÐ] Ò»Ö¾Ô¸¾©Çø985£¬085401µç×ÓÐÅÏ¢£¬±¾¿Æµç×ÓÐÅÏ¢ +3 Ñô¹â¿ªÀʵÄÄк¢ 2026-04-10 3/150 2026-04-10 16:29 by sophia_93
[¿¼ÑÐ] ¿¼Ñе÷¼Á-²ÄÁÏÀà-284 +28 Ïë»»ÊÖ»ú²»Ïë½âÊ 2026-04-08 28/1400 2026-04-09 20:08 by µ¹Êý321?
[¿¼ÑÐ] 283µç×ÓÐÅÏ¢Çóµ÷¼Á +4 ÈýʯWL 2026-04-08 4/200 2026-04-09 10:21 by wp06
[¿¼ÑÐ] 085404£¬334·Ö£¬Çóµ÷¼Á +5 sunjie8888 2026-04-08 8/400 2026-04-09 07:26 by sunjie8888
[¿¼ÑÐ] Ò»Ö¾Ô¸»ªÄÏʦ·¶´óѧ0702ÎïÀíѧ305µ÷¼Á +4 Äî³£°² 2026-04-07 6/300 2026-04-08 22:53 by bljnqdcc
[¿¼ÑÐ] ¿¼ÑÐÇóµ÷¼Á +4 ö©??? 2026-04-08 4/200 2026-04-08 21:44 by ÍÁľ˶ʿÕÐÉú
[¿¼ÑÐ] ר˶085403£¬291·Ö£¬ÓÐÁ½Æª×¨Àû£¬Ò»¹úÒ»½± +3 ¹þ¼ªßä¹þ¼ªßä 2026-04-07 3/150 2026-04-07 18:21 by À¶ÔÆË¼Óê
[¿¼ÑÐ] 325 µ÷¼Á +6 QQСϺ 2026-04-07 6/300 2026-04-07 15:17 by Ccclqqq
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û