±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 5022  |  »Ø¸´: 18
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

Çà´ºåÐÒ£°¡

гæ (СÓÐÃûÆø)

[ÇóÖú] ¼ÆËãʯīϩÈȵ¼ÂÊÆ½ºâ̬£¬ÇóÖú£¡ ÒÑÓÐ2È˲ÎÓë

ÈçÌ⣬ÎÒÏëÓÃlammps¼ÆËãʯīϩµÄƽºâ̬Èȵ¼ÂÊ£¬ÎÒÓÃms½¨Á¢16*16*1µÄʯīϩģÐÍ£¬±ß½çÌõ¼þÈç¹ûÊÇppp£¬ÔòÔËÐÐûÓÐÎÊÌ⣬Èç¹û¸ÄΪpps£¬Ôò×ÜÊdzöÏÖÈçÏ´íÎó£º
     cannot use neighbor bins -box size << cutoff   (../neighbor.cpp:1737)
°´ÕÕlammps    errorÖеÄÌáʾ£¬ÏÔʾ

Cannot use neighbor bins - box size << cutoffToo many neighbor bins will be created. This typically happens when the simulation box is very small in some dimension, compared to the neighbor cutoff. Use the ¡°nsq¡± style instead of ¡°bin¡± style.


ÎÒ³¢ÊÔ׎«ÎÒµÄboxÉèÖôóÒ»µã£¬´Ó16*16*1¸ÄΪ50*50*1£¬¿ÉÊÇ´íÎóÒÀÈ»´æÔÚ¡£Êƺ¯Êý²ÉÓõÄÊÇtersoffÊÆº¯Êý¡£


²¿·ÖinÎļþÈçÏ£º

units       metal
variable    T equal 300          # in K
variable    V equal vol          # in A^3
variable    dt equal 0.0005      # in ps
variable    p equal 100000     # correlation length
variable    s equal 1      # sample interval
variable    d equal $p*$s   # dump interval

# convert from LAMMPS real units to SI   # donot change this part

variable    kB equal 1.3806504e-23    # [J/K] Boltzmann
variable    eV2J equal 1.60217646e-19
variable    A2m equal 1.0e-10
variable    ps2s equal 1.0e-12
variable    convert equal ${eV2J}*${eV2J}/${ps2s}/${A2m}

# setup problem

#dimension    3
#boundary     p p p
#lattice      fcc 5.376 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1

#region       box block 0 4 0 4 0 4
#create_box   1 box
#create_atoms 1 box
#mass         1 39.948
#pair_style   lj/cut 13.0
#pair_coeff   * * 0.010325 3.405
#timestep     ${dt}
#thermo       $d


# for reading data file

# read structure from data file

dimension    3
boundary     p p s
read_data         data.graphite

#potential function

pair_style  tersoff
pair_coeff * * SiC.tersoff Si
timestep     ${dt}
thermo       $d
»Ø¸´´ËÂ¥

» ÊÕ¼±¾ÌûµÄÌÔÌûר¼­ÍƼö

·Ö×Ó¶¯Á¦Ñ§Ä£Äâ

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹ØÉ̼ÒÍÆ¼ö: (ÎÒÒ²ÒªÔÚÕâÀïÍÆ¹ã)

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

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

brucefan

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

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

ÒýÓûØÌû:
7Â¥: Originally posted by brucefan at 2016-05-13 19:51:50
Hello,

I think you are using the Green-Kubo method. The problem of lammps is that the heat current (also called heat flux)
formula implemented in lammps is only correct of two-body potentials  ...

Sorry for a typo:
only correct of two-body potentials
should be
only correct for two-body potentials
8Â¥2016-05-13 19:53:50
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 19 ¸ö»Ø´ð

brucefan

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

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

¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
ÄãÈ·¶¨ÄãµÄÎļþÊÇÄ£ÄâʯīϩµÄ£¿Ôõô¿´ÉÏÈ¥ÊÇÄ£Äâ벵ģ¿¾§¸ñÊÇfcc£¬ÖÊÁ¿ÊÇ40£¿Õâ²»ÊÇʯīϩ°É£¿

µ±È»£¬¼´Ê¹ÄãÕâÀïÈ«¶¼¸ã¶ÔÁË£¬ÄãÒ²µÃ²»µ½ÕýÈ·½á¹û£¬ÒòΪlammpsÔÚÕâ¸öÎÊÌâÉÏÓÐÑÏÖØ´íÎó¡£¾ßÌåµÄϸ½ÚÎÒÒÔºóÓõçÄÔÔÙ˵¡£

[ ·¢×ÔÊÖ»ú°æ http://muchong.com/3g ]
3Â¥2016-05-13 18:17:37
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

brucefan

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

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

²»ºÃÒâ˼°¡£¬ÎÒ¿´´íÁË¡£ÄãÕâÊDz¿·ÖÎļþ£¬¶øÇÒÎÒ˵µÄÄDz¿·Ö×¢Ê͵ôÁË¡£

[ ·¢×ÔÊÖ»ú°æ http://muchong.com/3g ]
4Â¥2016-05-13 18:19:05
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

Çà´ºåÐÒ£°¡

гæ (СÓÐÃûÆø)

ÒýÓûØÌû:
4Â¥: Originally posted by brucefan at 2016-05-13 18:19:05
²»ºÃÒâ˼°¡£¬ÎÒ¿´´íÁË¡£ÄãÕâÊDz¿·ÖÎļþ£¬¶øÇÒÎÒ˵µÄÄDz¿·Ö×¢Ê͵ôÁË¡£

àŰ¡£¬ÎÒÊǽ¨µÄʯīϩ£¬ÓÃdata¶Á³öÀ´µÄ¡£²»ÖªµÀΪʲô³öÏÖÕâ¸ö´íÎó¡£
5Â¥2016-05-13 19:08:45
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 394Çóµ÷¼Á +3 ºÃʶàÄ¥¾²ºò¼ÑÒ 2026-03-26 5/250 2026-03-28 14:24 by ÌÆãå¶ù
[¿¼ÑÐ] 311Çóµ÷¼Á +4 ¶¬Ê®Èý 2026-03-24 4/200 2026-03-28 13:17 by ÌÆãå¶ù
[¿¼ÑÐ] 349Çóµ÷¼Á +6 ½Ü˹ËþÀï˹ 2026-03-21 6/300 2026-03-28 13:12 by ÌÆãå¶ù
[¿¼ÑÐ] ²ÄÁÏ277Çóµ÷¼Á +7 min3 2026-03-24 7/350 2026-03-28 11:39 by xuxiang
[¿¼ÑÐ] 0856£¬²ÄÁÏÓ뻯¹¤321·ÖÇóµ÷¼Á +12 ´ó²öС×Ó 2026-03-27 13/650 2026-03-28 10:56 by self2008
[¿¼ÑÐ] 085602 307·Ö Çóµ÷¼Á +7 ²»ÖªµÀ½Ðʲô£¡ 2026-03-26 7/350 2026-03-28 09:57 by ÉñÂí¶¼²»¶®
[¿¼ÑÐ] 081200-314 +3 LILIQQ 2026-03-27 4/200 2026-03-28 09:41 by ±£»¤µØÇòÄãÎÒ×öÆ
[¿¼ÑÐ] 265Çóµ÷¼Á11408 +3 ÁõС¹lu 2026-03-27 3/150 2026-03-27 20:53 by nihaoar
[¿¼ÑÐ] 0856µ÷¼Á +5 ÇóÇóÈÃÎÒÓÐÊé¶Á° 2026-03-26 6/300 2026-03-27 15:12 by caszguilin
[¿¼ÑÐ] 269ר˶Çóµ÷¼Á +10 ½ð¶÷±´ 2026-03-21 10/500 2026-03-27 15:10 by caszguilin
[¿¼ÑÐ] 314Çóµ÷¼Á +3 ÏªÔÆçæ 2026-03-26 3/150 2026-03-27 11:20 by sanrepian
[¿¼ÑÐ] ²ÄÁÏѧ˶333Çóµ÷¼Á +8 ±±µÀÏï 2026-03-24 8/400 2026-03-27 10:18 by ÎÒÊÇС¿µ
[¿¼ÑÐ] »¯Ñ§¹¤³Ì085602 305·ÖÇóµ÷¼Á +17 RichLi_ 2026-03-25 17/850 2026-03-26 19:44 by plmuchong
[¿¼ÑÐ] ÉúÎïѧѧ˶£¬Ò»Ö¾Ô¸ºþÄÏ´óѧ£¬³õÊԳɼ¨338 +4 YYYYYNNNNN 2026-03-26 4/200 2026-03-26 19:00 by macy2011
[¿¼ÑÐ] ¿¼Ñе÷¼Á +6 À´ºÃÔËÀ´À´À´ 2026-03-21 7/350 2026-03-25 22:43 by 418490947
[¿¼ÑÐ] ÉϺ£µçÁ¦´óѧ²ÄÁÏ·À»¤ÓëвÄÁÏÖØµãʵÑéÊÒÕÐÊÕµ÷¼ÁÑо¿Éú£¨²ÄÁÏ¡¢»¯Ñ§¡¢µç»¯Ñ§£¬»·¾³£© +4 ÎÒ°®Ñ§µç³Ø 2026-03-23 4/200 2026-03-25 00:59 by 1027_324
[¿¼ÑÐ] Ò»Ö¾Ô¸±±»¯315 Çóµ÷¼Á +3 akrrain 2026-03-24 3/150 2026-03-24 19:35 by ÁËÁËÁËÁË¡£¡£
[¿¼ÑÐ] 276Çóµ÷¼Á +3 YNRYG 2026-03-21 4/200 2026-03-23 08:31 by ×íÔÚ·çÀï
[¿¼ÑÐ] 308Çóµ÷¼Á +3 īīĮ 2026-03-21 3/150 2026-03-22 16:54 by i_cooler
[¿¼ÑÐ] 275Çóµ÷¼Á +6 shansx 2026-03-22 8/400 2026-03-22 15:27 by barlinike
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û