²é¿´: 3209  |  »Ø¸´: 1

renmadou

гæ (³õÈëÎÄ̳)

[ÇóÖú] lammpsÄ£ÄâÒºÌåË®±¨´í

´úÂëÈçÏ£¬Ò»¸öTIP3P-CharmmÄ£Ä⣬´®ÐÐûÎÊÌ⣬²¢Ðоͱ¨´í£¬ÊÇÎÒû¼û¹ýµÄ´í¡£


###  simulation of H2O with tip3p flexible bond and angle
###  using the parameter in lammps manual.

dimension     3
units              real
atom_style    full
boundary      p  p  p

read_data     dataTIP3P-C.Ag
mass 1 15.9994
mass 2 1.00797
pair_style    lj/charmm/coul/long 10.0  10.5  10.5
bond_style    harmonic
angle_style   harmonic
kspace_style  pppm  0.0001

###  O-O interactions,   epsion,  sigma
pair_coeff          1   1   0.1521   3.1507

###  O-H interactions,   epsion,  sigma
pair_coeff          1   2   0.0836   1.7753

###  H-H interactions,   epsion,  sigma
pair_coeff          2   2   0.0460   0.4000


### bond coefficient, K, r0
bond_coeff     1    450    0.9572

##  angle coefficient, K, theta
angle_coeff    1    55      104.52


dump  dump_npt   all  xyz 1000  dump_npt.xyz

timestep          1.0

velocity  all  create 298.0  328093  rot  yes  dist  gaussian
fix       1   all  npt  temp  298.0  298.0  100.0  iso 1.0   1.0  1000

thermo_style      custom step pe ke etotal temp press density vol evdwl ecoul ebond eangle
thermo            200

run   100000
write_data data.*


###  now, NVT

undump  dump_npt
dump    dump_nvt   all  xyz 1000  dump_nvt.xyz
unfix  1
fix    1  nvt  temp   298.0   298.0   100.0

group  Oatom  type  1
group Hatom type 2
compute  mymsd   all  msd  com  yes

compute 2 all vacf
fix vac all vector 1 c_2[4]
variable diff equal dt*trap(f_vac)

compute  OOrdf   Oatom  rdf  100
compute HHrdf Hatom rdf 100
compute OHrdf all rdf 100 1 2
compute HOrdf all rdf 100 2 1

fix 2  all  ave/time  200  1  200  c_mymsd  file  my.msd  mode vector
fix 3  all  ave/time  200  1  200   v_diff   file my.vac mode scalar
fix 4  Oatom  ave/time  200  1  200  c_OOrdf  file  OO.rdf  mode vector
fix 5  Hatom ave/time  200  1  200  c_HHrdf  file  HH.rdf  mode vector
fix 6  all ave/time 200 1 200 c_OHrdf file OH.rdf mode vector
fix 7  all ave/time 200 1 200 c_HOrdf file HO.rdf mode vector

run    1000000

write_restart restart.*

write_data data.*

run       1000000

write_restart restart.*

write_data data.*
run       1000000

write_restart restart.*

write_data data.*


³ö´í´úÂ룺
=====================================================================================
=   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
=   EXIT CODE: 11
=   CLEANING UP REMAINING PROCESSES
=   YOU CAN IGNORE THE BELOW CLEANUP MESSAGES
==================================================== =================================
APPLICATION TERMINATED WITH THE EXIT STRING: Segmentation fault (signal 11)

lammpsÄ£ÄâÒºÌåË®±¨´í
111.PNG
»Ø¸´´ËÂ¥
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

Broad_Fan

гæ (³õÈëÎÄ̳)

ÎÒÒ²³öÏÖÁËͬÑùµÄÎÊÌâ £¬ÇëÎÊÂ¥Ö÷½â¾öÁËÂð
2Â¥2018-04-12 22:44:04
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ renmadou µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[»ù½ðÉêÇë] Ã÷ÌìÓ¦¸Ã¿É²éÁË£¡£¿ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[»ù½ðÉêÇë] ÓÐûÓдóÉñ°ïÎÒ¿´¿´»ù½ð´úÂë 28+4 1234567wang 2026-08-24 10/500 2026-08-25 19:15 by lfy8008
[»ù½ðÉêÇë] ijЩ»ú¹¹£¬ÒÔЧÂʵÍΪÈÙ£¬ÒÔЧÂʵÍ×÷Ϊ´æÔڸР+9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[»ù½ðÉêÇë] 2026¹ú×ÔÈ»º¯ÆÀ·Ñµ½ÕË +20 ÑòÑü°å 2026-08-21 23/1150 2026-08-25 16:34 by zsna
[»ù½ðÉêÇë] ½ñÈÕ²»·Å°ñ£¿Íø´«¹ú×ÔȻԤ¼Æ 8 Ô 27 Èտɲé½á¹û +17 ҽѧÀÏÄк¢ 2026-08-20 22/1100 2026-08-25 15:36 by ҽѧÀÏÄк¢
[»ù½ðÉêÇë] ûÓÐÈκÎÏûÏ¢-ÊDz»ÊǾÍÁ¹ÁË +9 ͼÀ²Í¼À² 2026-08-24 10/500 2026-08-25 11:59 by ÄϺ£Ð¡¸ç
[½Ìʦ֮¼Ò] µ¼Ê¦Í²ۣºÎÒÔõô̯ÉÏÁËÕâô¸ö¼«Æ·Ñо¿Éú£¡ +3 ËÕ¶«ÆÂ¶þÊÀ 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[»ù½ðÉêÇë] 2026ÄêµÄ¹ú¼ÒÉç¿Æ»ù½ðÏîĿͨѶÆÀÉóµÄйæÔòÓëж¯Ïò¡¢ÐÂÌôÕ½ +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[»ù½ðÉêÇë] ½ñÌì»ù½ð»á³ö½á¹ûÂð£¿20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[»ù½ðÉêÇë] ÎÒÃæÉÏÍêµ°ÁË +13 ÇÒÌý»¢Ð¥ 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[»ù½ðÉêÇë] ·¶½øÖоÙÒ»ÎĵÄÖÐÐÄ˼Ïë +7 Ñ׻ƹóëÐ 2026-08-22 8/400 2026-08-25 08:48 by ZJTJZ
[»ù½ðÉêÇë] Ö»ÓÐÿÄêÕâÖÖʱºòÀ´¹ä¹äСľ³æ +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[»ù½ðÉêÇë] filecode£¬4¸öjtjcÁË +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by ¹þ¹þ¸ò£¿
[»ù½ðÉêÇë] ½¨Òé»ù½ð·¢²¼Ìáǰ¸ø³öÃ÷È·µÄʱ¼äµã +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[»ù½ðÉêÇë] ÈÃÎÒÖÐÒ»¸öÃæÉϰɣ¡ +13 ´óƼ1987 2026-08-20 16/800 2026-08-24 10:23 by ̫ɵÁË
[»ù½ðÉêÇë] ¿ÆÑй¶ùÌ«ÄÑÁË +18 ÎÒ4´ó°×²Ë 2026-08-20 19/950 2026-08-24 09:47 by ¿­¶÷¹ãÊ¢´ß»¯·ÖÎ
[½Ìʦ֮¼Ò] Ìø²ÛºóÔÚÑÐÏîÄ¿Ôõô°ì£¿ +5 ¼òµ¥»¯xn 2026-08-22 10/500 2026-08-23 12:38 by ¼òµ¥»¯xn
[»ù½ðÉêÇë] ʱ¼ä´Á½ñÌ죬20ºÅ±äÁË +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[»ù½ðÉêÇë] ¿´À´½ñÌì²»»á·Å°ñÁË£¿ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[»ù½ðÉêÇë] Ó¦¸ÃÊÇÏÂÖÜÈý26ÈÕ¹«²¼Á˰ɣ¿ +4 ¹þ¹þ¸ò£¿ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û