原创力文档- 1、本文档共3页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
lammps学习L-J势函数涉及各物理量单位体制.doc
LJ势函数单位体制
Thelj/cutstylescomputethestandard12/6Lennard-Jonespotential,givenby
E⼆4rf(-)12-(-)6]rqrrc
Rcisthecutoff.
unitslj
Forstylelj,allquantitiesareunitless.Withoutlossofgenerality,LAMMPSsetsthefundamentalquantitiesmass,sigma,
epsilon,andtheBoltzmannconstant=1.Themasses,distances,energiesyouspecifyaremultiplesofthesefundamental
values?Theformulasrelatingthereducedorunitlessquantity(withanasterisk)tothesamequantitywithunitsisalsogiven.
ThusyoucanusethemasssigmaepsilonvaluesforaspecificmaterialandconverttheresultsfromaunitlessLJ
simulationintophysicalquantities.
mass=massorm?
distance=sigma,wherex*⼆x/sigma?time=tau,wheret*=t(epsilon/m/sigmaA2)Al/2?energy=epsilon,whereE*⼆
E/epsilon?velocity=sigma/tau,wherev*=vtau/sigma?force=epsilon/sigma,wheref*⼆fsigma/epsilon?torque(⼒
矩)=epsilon,wheret*=t/epsilon?temperature=reducedLJtemperature,whereT*=TKb/epsilon?pressure=reducedLJ
pressure,whereP*⼆PsigmaA3/epsilon?
dynamicviscosity=reducedLJviscosity,whereeta*=etasigmaA3/epsilon/tau
charge=reducedLJcharge,whereq*=q/(4pipermOsigmaepsilon)A1/2
dipole=reducedLJdipole,momentwhere*mu⼆mu/(4pipermOsigmaA3epsilon)Al/2
electricfield=force/charge,whereE*⼆E(4pipermOsigmaepsilon)Al/2sigma/epsilon
density=mass/volume,whererho*=rhosigmaAdim
Lamps中的LJ单位问题
LJ单位
mass：=2.43x10*kg,即［mass相当于2-43x108kg
sigma：=lxl0_5m,即1sigma相当于lxl(T5m
epsilon：=8.41xlOl3N-m,即1epsilon相当于8.41xlO_13Nm
那么根据manual⾥⾯的dipole的解释，即：
dipole=reducedLJdipole,momentwhere*mu=mu/(4pipermOsigmaA3epsilon)Al/2
所以1dipole相当于(4pipermOsigmaA3epsilon)Al/2
即.y]47rx/?(?n/7Ox(lxlO~5m)3x8.41xlO_13Nm
根据您的解释，这⾥的pennO应该也是⼀个约化单位，即permO=l,但在计算的过程中要转化成实际单位，故应该此处的
permO=7
4^X10-7N/A2,也只有这样，dipo
文档评论（0）