±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 5024  |  »Ø¸´: 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µÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] Ò»Ö¾Ô¸±±¾©»¯¹¤´óѧ²ÄÁÏÓ뻯¹¤£¨085600£©296Çóµ÷¼Á +13 µ¾ÆÞС±à 2026-03-26 13/650 2026-03-29 08:20 by fmesaito
[¿¼ÑÐ] ±¾¿ÆË«·Ç²ÄÁÏ£¬¿ç¿¼Ò»Ö¾Ô¸»ªµç085801µçÆø£¬283Çóµ÷¼Á£¬ÈκÎרҵ¶¼¿ÉÒÔ +6 ֥ʿѩbaoo 2026-03-28 8/400 2026-03-29 08:16 by ËÉ»¨¸×1201
[¿¼ÑÐ] 321Çóµ÷¼Á +7 è±Óñ~~ 2026-03-25 8/400 2026-03-29 06:41 by 544594351
[¿¼ÑÐ] 316Çóµ÷¼Á +7 ½­´Ç666 2026-03-26 7/350 2026-03-28 21:28 by sanrepian
[¿¼ÑÐ] Çóµ÷¼Á +7 ÕùÈ¡¾Åµã˯ 2026-03-28 8/400 2026-03-28 21:07 by ÕùÈ¡¾Åµã˯
[¿¼ÑÐ] 083000ѧ˶274Çóµ÷¼Á +8 LiÀîÓã 2026-03-26 8/400 2026-03-28 20:33 by ¼ÓÓÍÏòδÀ´°¡
[¿¼ÑÐ] »¯Ñ§0703 µ÷¼Á 306·Ö Ò»Ö¾Ô¸211 +4 26ÒªÉϰ¶ 2026-03-28 4/200 2026-03-28 15:30 by 1018329917
[¿¼ÑÐ] 0703±¾¿ÆÖ£ÖÝ´óѧÇóµ÷¼Á +3 nhj_ 2026-03-25 3/150 2026-03-28 13:24 by Iveryant
[¿¼ÑÐ] 291Çóµ÷¼Á +15 hhhhxn.. 2026-03-23 21/1050 2026-03-28 11:26 by self2008
[¿¼ÑÐ] 081200-314 +3 LILIQQ 2026-03-27 4/200 2026-03-28 09:41 by ±£»¤µØÇòÄãÎÒ×öÆ
[¿¼ÑÐ] 352·Ö »¯¹¤Óë²ÄÁÏ +5 º£Äɰٴ¨Ly 2026-03-27 5/250 2026-03-28 03:39 by fmesaito
[¿¼ÑÐ] 086000µ÷¼Á +3 7901117076 2026-03-26 3/150 2026-03-27 21:34 by Jianing_Mi
[¿¼ÑÐ] 287Çóµ÷¼Á +10 land xuxu 2026-03-26 10/500 2026-03-27 15:33 by ÅÁ¶ûÂíÀ­ÌØ
[¿¼ÑÐ] 279 ·Ö Çóµ÷¼Á +4 ˯¸öºÃ¾õ_16 2026-03-24 4/200 2026-03-27 15:05 by ×íÔÚ·çÀï
[¿¼ÑÐ] ²ÄÁÏÇóµ÷¼Á +5 .m.. 2026-03-25 5/250 2026-03-27 11:08 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] 342Çóµ÷¼Á +3 ¼ÓÓÍaÀîzs 2026-03-26 3/150 2026-03-27 00:29 by wxiongid
[¿¼ÑÐ] 294·Ö080500²ÄÁÏ¿ÆÑ§Ó빤³ÌÇóµ÷¼Á +4 ÁøÏª±ß 2026-03-26 4/200 2026-03-26 21:14 by XPUÀîÇì
[¿¼ÑÐ] 321Çóµ÷¼Á +6 wasdssaa 2026-03-26 6/300 2026-03-26 20:57 by sanrepian
[¿¼ÑÐ] Ò»Ö¾Ô¸ÄϺ½ 335·Ö | 0856²ÄÁÏ»¯¹¤ | GPA 4.07 | ÓпÆÑо­Àú +6 cccchenso 2026-03-23 6/300 2026-03-25 22:25 by 544594351
[¿¼ÑÐ] һ־Ըɽ¶«´óѧҩѧѧ˶Çóµ÷¼Á +3 ¿ª¿ªÐÄÐÄû·³ÄÕ 2026-03-23 4/200 2026-03-24 00:06 by ¿ª¿ªÐÄÐÄû·³ÄÕ
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û