From 33b919a5b64db5067e3e01a6faedcb6dca14663d Mon Sep 17 00:00:00 2001 From: athomps Date: Tue, 30 Aug 2011 18:01:33 +0000 Subject: [PATCH] Added equation for lj_cubic git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@6858 f3b2605a-c512-4ea7-a41b-209d697bcdaa --- doc/Eqs/pair_lj_cubic.jpg | Bin 0 -> 22768 bytes doc/Eqs/pair_lj_cubic.tex | 12 ++++++++++++ doc/pair_lj_cubic.html | 34 ++++++++++++++++++++++------------ doc/pair_lj_cubic.txt | 34 ++++++++++++++++++++++------------ 4 files changed, 56 insertions(+), 24 deletions(-) create mode 100644 doc/Eqs/pair_lj_cubic.jpg create mode 100644 doc/Eqs/pair_lj_cubic.tex diff --git a/doc/Eqs/pair_lj_cubic.jpg b/doc/Eqs/pair_lj_cubic.jpg new file mode 100644 index 0000000000000000000000000000000000000000..e52392b87e329e56e233512df0f2bdb0396444cd GIT binary patch literal 22768 zcmdSA2UJwgw=Q_P$vNko5dO6~@M z78-cR-~YaQ-@KW()_t>P?yM=g`s{N~o$A`Bs`lRBSHOKPgm+5>K0BZt(DS^Lr{V$9CU*8~eaP+kY0EjMl+s?t)-UF1IKv~e<+s6-- zw?O%Opr1D=L;qf<^Z^?P%Afzr&i|5^{>pa$k`e!uG1ON901Qk3z$LV`^>GFO0#HMM zHNf5(tcTzp0AQ1QI=FfQ0NDdjmbZ7YbpYQ*6jJPrJf70K#PUz@qWB?w?!ORG#0A9dT zF#7RJ^={eO#Sg{8*L1UUIz4s_AJjc$Zb=p0M!W+yB2wf5+#) z8hfh$&$xq7kN?LtXAuApM}p&n=RdBwW&%LRV*tQ;_8-@{2?2n(7XU`m?frcM{#765 z-^0JgJtaT~FasO_FCYMj01|*Kpa`e}4*?y(5HJI*0eip&@BrU(AP@$;0Ahf6;0^E= z$Oa0465u^h3p4<&Kqt@(3<0CScVGrs1lE9UU>`sL=fD*>A#or?5DEw#gayI{;fIJo zq#*JTRmdZVKEw?21mX2+DnYfNhEOZ0GxQlW1R4cRg1&nPNF(1!BcuWnh(KHDe88 z&0uX~onvETQ)BaBOJZwan_;_QhhisS=V8OKKVeT`Z(^U};NZ~V@Z;Ra(ZR9B@xzJ1 z$-=3@>A{)A*~UTP65z7nir}i@n&W!nM&f4RR^fKzPU8N=y}~2KHVVktmT^lZ23DkTjBvlN^xZk#du&klK<)kmiuKkWO^Hp(L#avWMj208 zOZklwK}Aj_MrBMDM3qC;MYTqaLCsCAN&S>MnYw{`mim&0iAIUWktUX=j%JePf|h|+ zk=BtmjuuWkLwiZbN~cQaPM1R0O1DgpNzX^GM;}O^Pd`Y1z(C0$%V5tC$I!sAzzAjJ zW7KC1Wh`O*%6P`a%%slb&6LH|&vd{{&8*1m#+=UF&AiJ(!6MJ%%975~!?MRp#j3>m zl=U6!04suxkxhfmpRI`P8`~8-54#b2Bs-janFF6gio=N`ouiKf!O6m@%^Aj7&AGsZ z$0fz(!j;K2%yq%d!)?O-io2Ej7Y{9uCQk@YHO~?+5w8O8Gu~p}={wkWr0%%i$-DEN z56UOb=fao6H_i{`m*98h&*PuCi*;A}uIJt2yK@4B0!jh_0@VWRg4BYK1S17o1&@R{ zh0KLggocG~?up-fdav}}qA<%^Q}$7AR=!Y?Q3+T1^Z?_5`h&O!-&Dy}ja0K$SJc?moYkt;5bEOUq3WMB zur#zZQZ(i?nKbP+D>V@hB_D=A9MB@vGSJG_+IYnG$mdbJHdOnecAECmW3I=Zk6U#B z9Zj7yon>7f-DkQTdf0lpdO3PK`oj8Q`ojhk1~vv&2Iq##hRKGDM!ZJ;Mt#O4#+JsF z#up|JOj1qOOodFtO-Ie>&0Nhsn&X%oo0prPS*TjPwb-^4w~Vozwc@b~vKq0bw|;8f zWkX_PYt#4y>xs#eswZe$9ou5tGdoSYTswrlihZX2zJsDey2GyHeaBSCpHBCkQk{M| z-*--P-gQxM$#6MvRdLOBMY?IY6}Vlv>$q3A-##^Y3irVEc;eCKN$%xKz4lDz+1qCbA1$A9Ux=@j??*pMFcA3a&+8xUzaAhTkQ0CkG!ASGA`9{i`Wnm^ z92@*IL@lH&6f@L5v@eV;>_ymmxKemg1SG;Pq7TLai-K)GSAYKg1^x^77h{q4B2y!O zzchK-9>o;(JZdvqJ-RxEIL0St?v?DTg4dX@U0;8T6^+e`y@_*(`x1XIJ|iBTV4v_M zQ8@8k;%$;s(zj%Z`3rYNMmPbEtYN!?7-PHRbLO^-|e{nqB~NQPKOK_-5t zf9Be|NAEsnab~^ALT9^V&*UiQz;hXMV{e9OQ(d^?y+rzfL zcG>pk4#AF!PPWdRE~>8NZv5`3?%SS_p0iKhpALGRdpG-R`WE|5`=iH|jZllay$E;;`3yY%<=3Au^xNrlPY zsRvU-)0)$xGmmE`W{qZl%vsK@{c!m4YyRmxav@;hYVr9J)>7Ou>2k&j<4W=BomKdn z#9G(-gY~Z)1{;f;_M7`#zFX++s2!r6jGrt&D}Ran>fBY^{k~_hx4rMZe{~RZNOqWe z#Cz0qta$ttVT#yB`XFyl5>Dw(D}Iap?mN>xTRHbQzq*J+(V^b|k@_=yX>$4VDiDo> zes_K6y8TA$X6e@B_SW0h$M&CC00aDnIynKrNd*89nE?RBcL2aA{#X3?PYE!*1j8>V zC>uj!q0ca+FsU)Ou|8v0;3VS)<2mA65Ev2a5g8NflbDk}A#)-3q==y$q$Z=Wpe?6E zGRQH$Vwz{U#~Q&l&B4$4l53eqhWFi_bAH3SO@dTH!S~ih9*H)HGfKor9!nd_bjq>c ze=UEcsISzlLi->{bzWUcBUAHC%kWXP_NESnu9BXIex^aM;f^uB$sJQQGi!4ni?NurG1QcYNanbM|zxc75z7=g$9>!Gpl#(sS4AhxhQaR-X#r48Lf9 zzX1C{{UD{_yCF28kkEs$x$xeI8d&D@7cV>_O<$@+2}aY$V8tL`ExrC6+Z0zApOBE4 zn4FZ9ob)CsB{4NV?R9$O+wcs(%%|_{vW&Aeb7XVx<}v0InQ>Ja+bNX?);&&BWWu zh$-i3-5Hr#jya;aiy!OrlM8)|4NJw#DJ#!cz1FPOA8v?mvTosToovtV^!zOS^=j9B zPkUeRfbd}Mu<*$JSQ3GSm_yc_M4dkQt#rn9PI7*Ifk6HGvvIj`wTa%jKD_yTd;3rS zt3s79gfNaUd$C?&8{i1yQsMr@8^SLqNFWR#vLV(WktAgzBOtq>5TdlEN~iuxi$Nzz z@4`^TxXet?qQ{!V_LT#lQ=KcCdxRJF&O^Q={%HXgK^LKBVGI!+(JZk|2?5DKsXiG3 zSp&J8`|ApK6g`w0l>a=CSAD70ss2M_R}-a${fJDP`7w`|+^6`H3jN1qBS@3Q2OZ0aRPzq!S{2lZyxFRGn)FMnQoG^SRq6hZwx$g_pNcoo> zQMge@(bF;QuL@sB$GXSq#)~J=CEO-%C4Ek=dy|?HkZO^pkk0xRc)Oa>lUexgMV4*0 zat=o>X6|<0aDHV$e4%HNZn0Y28b>!-og; zoDG-_TaANFWzA774y_M9a<*Z&9kkDMbaz&ErFBR2xPLP1RqvDR=N+IQBp!kep@xq> zZ;vd0nH?SfIx_b4+xYmz_o<25$sbb-(@Qfevukr3Kepy~7Je@7EgdW)R!&#X*Zyok zHo3MQZwK#G{G8t<+k3p9axj0ycl;8udGhGA<4p7%juQRTaixZyyK()u{uhMoL-Q~U zFmW;4umZ3pu#q?)aU<|F@EP$D1OtQ_L~g{YBy1!%WZ2}?6nvD@RI1d*G;Xxd>E6=U zG7K`VFkP}xvfg8R%w9d%Y9cP8abh*%V-g2a zgwnh+8nO;@u=~04bqdvr`AW&k5i0HvOjQ-txYhC0e`ySA=0Eh)QhP-CXidB7v5$^~ z&ZTabUWC4k{+YoCLw6%CqgCT<6GKyC(_yn0<}&7|7B!a6RxDO?)^BVypFo~;+J@Lk z+g;eVIQTh=I3k_k&K@p&F27vM-5lIm-It$cdl-9Cc#eC;duu$yde-mr!dKoG?bqfX z6d)de46F-s59SHp3@HFd3U%0Icv6H03-QZ<9)u4~4*wZp{bDt`KK6QC?ECg) z;k5It{13hbfhFyg>=ThcaSC!1s+3*p#o5Cs24N^+64Upy~be2(8Tb@c!$x8 zu@CwpWlT@ZT+A<+f3SE#To#AbiFJ(4fo+PNggu0PgCmLKgHwyMhs%TOgj)_;t2=le zc=dP}`11HM_~Qih1dastgxG|(gk408L}5g$#45y1>JQNFc4G~k^xx=J`rwFE-@Z)ehCptSt()ZJ2G6d9CDoZdF6!^ zWE3ANSt|RfBtNKAom9Kj;MO#L7^l_u==$*komAaTeJO)P(EBNw6qw$c+gc1+30P-5 zf!KQ4tvcvC4m&Hlw7Ds{fAz5RM0m&gaQlwR6mE5%r% zxVePlB(FD$sbpzeZ`(4Hvs`ji@>uh43)hPWOY6$ty^pT)uQ9Jv`5@H5(nQ{Z{}HPl zyOXG!?h}8X%7E=q%;$zL8)J0iMiUv+^s^Q7qKh3XGHZ>S?Axz?ZtRI3gdBZD9G$YB zJ-qP!lXcZ`z4nj&fB+a^{}Tcf029bJh=Ht)0pJJ(0f|5<$RW%Drw~#IFGLw)2?>Iv zLmI%oKY>y}#i6FqFc9B=h9WR%F_bYpFtRX)FitT!Fby!HKzzOnqH;YDi+5rD#umVK z#xB8L#bLs+#L2~3#%0BIz^%kZ;7Q@V!260%i|>fvNPt0LNKj4)C3GhIMkGa43ZmYZ z#1|wUBtJ>*NjJ%y$&Sc_DR3yVD1|6LQ#pWmw4BC>7N54C4n}{U{)VBCF^T@&?FJ~INdKYtW1m@MHh zx31c+<84%JX>70j^!#rpHT?bQoc|)@4{+&xHFhn2lk;!=j|t!b zWB>!e12Ro&;N15BUIIBl3or>Fz`nl=(F7TxS0D%U6@q|LLFJ$h&_w7*=oSVUh9WrQ zi!i1zaWLgDeKD&rH?f$&*`AIyiA|1efc+MG9)|_T2?vgI39=~ZxLbGH>W$)SuycTR0+7V_HF%iuaTa@6G z^p@(DA(M5L8yS93 z1d$|)WK_yZ+VI`=k@I9M-aWG&aE!+ zFPczRf3W}5Ug}?>uJX_d==JNDHyk&8w@?1|zps-h$XS93LZ}SVk1Mygr`#Y{nF;{c zDYv&*#kaTDB_P-UIqpwS|HV!JEsPEF*|{x$r6HB5=YNm?W4c`jo~s7Ag7`&S8{h-) z-~w0>a?s9zzw(zeh0udM=U*8B9zq!YUDk*2{#zLk*SvBC{t}{p%Yc0QKM7nz_~-mz zM`$7Rzn6gO6=M8%8HoN{_I44pV^9d>?*|Hgz=VN?@i$>(Vq#$7V1vXT$PD7(65!$B z;^Sgt;}PNE6A%y*65`+y5fc#-gCfD-O8!2J@z1kZ1i09^AeZ;QnQlJ;WVk>yR1X8f z4nWBu7-W#!0e}V61Pdf$!Ss)6e>H<*U}AwS0^0!=XeI@lkAVRN8;gYrz6j9qLO~tL zu*lhjm9Q!FY;oA1QHs1sDZ=GYuJ5JNpFnVm+WAD{;ZxJl($RBq^YGr`6BCz^l#-TF zd7!GMuA%wRz|hFp1f-zt9UPsUU0mIK{rm$0gMvd|Mn%WGdL0{=nwI`HBlBHWc5z8* zS$W0#%BqIOrskH`k8SOJ{R4wT!=FdKOioSD%+CFoUs&JR+}hsx`D=F%d2;&u?EC`t z=kl*!pf~xKTK}%u|D+cgs23FM15BL1dO@Ipe-$Ug#9|Z1CRftKv3*9tF7g7GQaPol zz88-}R3AZQ=QDv%%_+9dh5W19KbrmL6pQ@7rP;qL_8)pJ0E8gghLT~B0Sdr1x;Snl zg6n3ohd6@JUDvgTHVn2_fVU?pF?7OU8JI?j+!^cU9QMp(bagSN|VaO)>f!y&;hpX(o|NDRRATb zcqDAIAVwa;bwT`MFj2KBQPy&%_P2(R^@>Sz%--7D7oOazcb=wdP;;f~zjjFFz3@uG zFo=)ecDGpk+1%y%(Cy*A$Bo6({`=xM6@1>}=l|sy{p}Wr{2r$HT~6s1z+${kx?v95 z?6+n|K7Li?9`b5noNN&l+E3HXJ}t05Rl6W?KmXehopTWV$yzmOv=Y!-8Ml7aCUrPg(Tzr z%cOMv-k-XZvza7`)b$en?anNtt)%qATOeHfw@%JZ+Oxgxg=Co6%$I<3BiUvzG~WBs zqxiHx%%g~EWds~eMYEoiyjj2IPD1b_h+`wtO`g2K_`CJEd90ZGR?PZ#bf1>9SgXVw z;y_m8_4^b_nxyVe@Hah9$++mpsL?BQhgK{kDJ&-l=+bgEK zlZ&G#A8{NKn;R+!ahd0>d60qUddPx-p4XC>bI}Se1qt1lA39IawG_%th^ZfD9*Met zw8Kj*9|uwi27Xx;&9(I|kqj+Sb(hWk0lwbKuzQ@4^;;l0>=I3ed~*vdX@@f>aVnZf z+2$+ZoZRRC3ArLbtKbwYte(G?AL{R2x&^koY}AVn73tzh?`iJA___pk9=B3=J6*r; zSt~QwIha81_siqw^erT{64%0i$umrg(;j6kdkFbg?!;1BY#g`O1L8@2TOpo4wp~W5;-i$tP@5$clc31uBzU$MPFlQkfvZyl*6KJkKU?J#n zuVvYoJ9I{Btt-TcjXhXFb2m)oR#p z)YJyoYLpz)<}HA^w{rJGuzlOa-Lp}%GlqcIhi@HPGzsH$BWbA^bSzCqSP+kdqFBI1 zm550gwSUioNiX#cQ6uWbhKI#q`}XR2Sw;nM8Rx|A?1Fewo|p(-fKbo;s`ajE8{Q}gMf53Q+oHGSf#PDX-8jL|q!`l0lH;?KV$1#_-+ZzgVm zhM)aj#H$wj-B4r(JQPV`*gkcnI+)bxSky~2n5$+QYsTiGe5zegsYc!VSmnt!tXC^S z3}td&j3O&Rv!gWYOx-UOC~oe!-U3A8uRToL;$*c-yXHI=lfpU(WSe|MAl{`ft#}U0 zLT&-vTfnPl-e!0qx6*w*KK?pzeFj|v9rloiyJOCB{YXyN^e5*}G7XDav&G5QKyUtPtD=}S-e8u*T*M0PhzPC(xI4%kg zP@6;9HvKTHy&HDq;9+-8@h=GjH<{?ABXWL=Al(tpq~YgPabGnL_!a#1(yAMGP8!-x zB0?T6Z){u$iexs&m3D29wwhE=MSw1;HO*T6C9IcQ{@#8MgTh+D4KN+SyE0aZDp5%~ z?OPZ`*S4K8bn%^sOIT#w|2pR6IDRb#FLB=yLU9yCH472m0==?wSm7Is?$!La0KH~q z(8Qsq=*|LneWjIK%D^8W#`IWELvD-yGN9f4G*r_UZ8A3p-|(H<+P+@xKZo^CAwKA- zQpMU@!>%5NJxg1o>3`j9Me`tMU}2@Q(%fZ1Esg2RVtD(}j$^FcNS3Jzp_{d64J=2Ppnjb&Jnw5Tx!}3oW3Or?NZObq%&DxZD>1wQ? zA2GqgPeb^he~ojh}x^!L8&%uhE=AHGk#LN`gA?~?h1X7HI77+AHm!A(*253_H9fy|cd z@C1O)Q_#uOE`c0aI#6VMBYj`nZ^NFpO*v;_@A_;aORu0wFSo>fVgL8Fyoo#$3hzo0 z?FQQOHnivz+WAT$sH0y2(oNuQ@85Nv80ff>%(SQXGd+OWDX(}WiC%)Hv53o5iT#ay z$@M#!BMI)!2s2Ulq&!t(z>$X=VlQNOf_K?=d6QT_}`n5Mnph0m{+MUy{hp@}s77xW$*9fhZg!S^G z@y0?f4t%SEE>AzvC&CXQo8zxwdIg{&P7&%nklRW{`^ zT!Pd>ZF9@&%pQxq_Umu6S`B#Cdq;bheyyX3uEnF^g2P?VZKzH4*!IwvyVkBLkOy6S zzWFTENzaL>D}v0xo$9Ap8(Pw{;e?2f1ZwfHN0#;eCajh z6?MagvVpDT*$GH2!|NjEoa_3zY7Alz1;gIY>c$wy=?HyQ^<#60Sdu(!yOKfmoF`t_ z-2zTxQ(`Wj-3(z0NH9`KMC(Lc2p6?>HhoBC;$y?}IHArz*Ayv?KEadX@?DHrE?|-~ zKxnOPf4DBV`&4}V^$N#~hu<>BgX~bh_L*r7o*oO%wtU*M_V~xIbU!9?z9ntKc*eo! z8?Rdte6*-F@p1lgHCyI)u{~}dFFMd z-zsg;o=NEbn2l)z*?2`f!RUjg>7!cel@rdu_`_!eFg~Z~1R*|I)e|cC)2eh!7XGF5 zBT7A0jlMYELg;@{COoD=sc52Jtb3>{BabpG+_Smip67v=5k*fkC0akUzYQvnQ`CRp z$Vv4?2?x_82PJxri=shn4w!S-BTm;^`#f-r5mQa=sz-V2MDnvg9|q70c--}liF%{K zh5H$9RvcZ@zQkNL?sdAf@tb9S7z}cvySXic#``1i5Vf!3>k@XXoQMQzkKZj(O5%JY z=%jmeFAG6iKp^O25lx_nWBcHrOfsWW1ZNBqV*>UQkjagJigVHC>in z_O+VF7h}2Zkg78(tFt!ZPEgzps7io*zL~UXOiH{m2!fK}p|!eZ4EiQl`~7Qn15!D? zQAsOutv`hDly_!UM_oA*-+ZqB<2l5WD*2n9uk=rW=rL@XzWW2k``R;4OFEm{@O!5<39CM% z9F6I4u$1=ro#M3b@C%nuwaeMxF!f%^!6vX!!nZ)w=lvPWzoZn0=6?vC^-2nOknfWp8U>_8_S0B?cPO>I`0AI)VP?miS!%` zN;%>6aM507#`*~~U4^MDy;>&JeBN~4qO7aAl($GJ&~>6htiQ`5^Dnv<26bMQrLL>BoR7*=qZpst`2 zp9S&YGxT@(u~c(a>qdGHH#oI~qyJBE^M8h)|JiesShVGJPbVHiYXfXvOqg?Nm}SRw z)2h(7%m*m54G+)3)}aQW%wYQpEia)jYfPyUvtM%WvGUixavFbH-0j7JFcONvd@iSe zJT9<_low&SiiMX<47|qwNhtc6^+vtSD6fwHH+3BUj@Lr`4Zcz7%@_-Fr|4Q`r@;Gt z@4bbam+M_BdP%z;vy&I*`7Te_niJEn8yqZxJ)ILof+emvM>7rY2)e2f1Wa1PRWG{k z`!O}8w{tH3wkUWsVXcMeIB%c05r~+DQFc?H>>o{QANRQnFo)Y1dh|GiCaT>jPnwkv zq$8&-NlAY6s3%!-9ll_&KAtGgZiQ0mT~RAGd5S&}(%`nm5*Rc=Ndo zm+ggK)7cgY-`C>vvAdzaTQ@5g|Ff{ZI_xp3f1M>VgwPy5MGct zv&K!sQ*Qwe8x&uzLjBL7JdC~o7D72IB%Xhs%_GEono?VfMeEF&)G=3MV@u@vrsnIU z-tBz#46aD&HQ0)$EzA3e!CRmS9tRF|8-k!duY^^CFjYi_1p*d#j24+SPfxzzQqR<; zogE$B-;ejPIbCXzHQ|nRn!X^_my|vBPZteR~8*7=877>zH-884hONX7$ zkmLOL^gKd+3GIwqXk^1m^Ga5B2^u=%2?sH7SG}yP znLvlmh2ppc2yTIM7F1%<4@-%HI-s`>T^m>92g^!#8Td977#_^jv9+X|tAxeTTLc}^J@6f&UQ5NzLbrbPYMx1d z#)8mONH(uH=(@rSvqZgIa~WvAh?e7rg_l)oudU#kTBLS;N*~hrjr~zJJWvY^`%G^7 zSBb({*1DK=DBXn>v&9p;1`!Q1-wIm-Wim>_LjykFn7JH^^Us&@e*H#^mjja72o0NQJnVm#ni-Z{TUV2iTZ&%*dH=TmcztNl2;VFHjoH* zl*#(sAPK>8P_icSaTW4I%LzSc{fY{-(|FgLJD}BGI@K^aFuiJ?#dO6@;xfnpCU*~e zcSp+dTGWJBD$E-8()H3UD}hg~y#A&x&AMH@UYttnhFn|r26oVMdQl@HY%k*gz!F)W=rHgd4G@>w+HV0>+@5TcF>XtQqyfu05cFEMD(> z-9W8_10KaVb&+36aEfWED0u?G=V&XFqg6b2_(n4RZg$>a;!!4di7VBV*$Ve)yz(#F zK|XM_)dhSwP^xeS)wMF&$qvtMqzFiRxq3r)qphIsi}A^Z2`zvW8om)iz@yGMcw*)7 ze2KbS%g4P0n#T(B6-4kNA67X=8j&;vT|79XtfFnw=Q}$)u{kVjT!S4Vuf)P0A@;8` zeh$ejulB)7t4RLfn_8PL+swz}AjP}R^An=32irdmU-nxtL&4t6*oBeyc;t{oWlKi; zlN=mAbw|2dvt2!PK1}=N+wbwd(d=-Qs%lxAtH!s{Nwd}&DomSX#z@ojaB;7oM&Iie z;BqtG4Rv=Kec#3ab5`p)+?njT?xN#b=HG2NQ&_k-FGXA&SX@b=N;aGa`NNGb)DfL{ z9P_(xYUbunM)*{Sx|*JTf!}QyWppches?p)`jxGpSB-5}aR!E%+m5s$PBMqADO8c8 zXs6p0anw}`S}QX-j`d|{jkqgS1WUh8e4tmLdnoxMIS0Qxb!@GvG@8PU7zP@G*a)rg zcvoEp@iUJ{{>h^IPCb2-u9qEMGMw!#v1uj^yMEUOuJE_5`!Lozl}*WEFM=@0_oii) zjf~Q#1{N~{P1+*39^365)5moTPmRxK5!ISvuL!r;^C}rrC8pIK_xHE<`AtyA;Gj*O zQIDO}^z(YoLUQ0KW)VztuMQ5kwB{OIJKPy|EWeApL}f9WGO6-d=9_S+J1o@;*_`P# z6+*O&r`3SA& zNFsJ`iRvWVKOrkD@y0cg^UWQN<;?A48d$Ft+~hn97NdZDLyXECmdU83{|6=pYG-eYXQmDjLb(}L89_dXq=BOoP=uD)5k5uE|;@hXa=?B@_a9M8RnFoF#w zq!_=E|EkPgj}}DdZkOgcgE3ZnH5fGp7Trb!{mQ0J-QzNT|UaDQQh;+3EJPR~l z_UVxHDZ$}`lX$@lGcVpYg|}J<1Zw`ca*7*51#a@nv!fNhcdC@}3)FR0EbiT_CL~Q$ zYv7ef*_zl4$b>x}Joq6Z;Bh;Nu7)nrKC?^>QZHfx5v={tY8?;=)SB>=Gpa0)FHsHl zu9vR}`E?@D3FMBjX_P%F$>m_4`98Nrf~NEx4rv`uPAX_LIpnEuyYLRe=r3ULBf$+~ z|HcpJThyC$O5Ir1oh@auISk}Y=~MLf5U_uXOQUapmN}q6Ca18@TkA zbW-2!(YcgU_p^i=-UE7-&+AJ2Mp)xodBU?^G(oIm~@3VDhem z=(MEkKyy=_*Q;cHLQ(S2Svdowf^q%#PTBP45r^lW~qw@?J)S z!m;GNM0G=WM}SS{%63rM_RI*E=lyo>)y+R#fm`UTvofkpSUv2Je*@%G`&RDdY)ba6 zl7-!wMYC*l+yaEn8D$*f9yRKo22Mq`-iSU=0Y_qq&m8z>me5b0^3A@rX7gX|J35-< z6s1kx8YX|r_7HbJ(?n;tqwzu{3e@J7i`?SM^I1P+sG7oS=HEr@a_ei}2X@*N7DPA0 zTMteVsaFFR?cv|*y;~(hoboCx+Do(Zi$Vj5-sUxQgw_9PmAt2B)on;6=r)S4eK>W>s>wEU%^f5?P>g( zbPKGLB1sa&NigAJ4Gn6*o8P;BT!xu* zcN@F@NiEIIwb8ygFWAUFaVA)^NnQ{dN`79A#5fOyp|r*J_wkN?e->CgXtdF$K9fjl zjFy;+EpcI%jug}9zx%?>l!W5?je@5aP1w`xGI6088%ngo`n%>tmQR-x}M&wrDUhDl$yE4b-Iu>fc#ny5O^M#&6s zIL8Psj0vtGR@%bP4<$~1S$zEQYmQpN#DjrQuHo^p?3lj|px~$zD}1vk;|0c^2>~mF zh^%^4KsWj9U9{}``v>{=rt^3BW4bzF^$Ue@OT+hGZaFJ>s;qPXtF-9ec!mZ$6H+ZrZ|olC{9%;D-nyr z(o5A*o+;img$%F>4LFm(E^b6~)%C>5@t3;=QcJd5 zW+#@KXgdex2nrAs+EZ4e1=|+->#>{g7DPgg9tsjweAi3O*0;D<-0g~r1$~tw;t4Vj z1dTlo*9-#+EMc&UFSBzYI+5i3@l5QtZww4(x%@tSV67x#ObhOBdgM*wb#-q6{cIY3 zmgsxl*MSNwQ+Z$NUW(G``>?c}h|4eR)HZWO=yYj<1;KOojBR?c*t8wl9R9@1Gv7+NA9NHi?c{RcQHY^WhfgGuV&{6ZYNMA8eTxzjN2VTZ{LwM(m0X z4i_xXn`z-w_1Z4P4;+)j4d8Ez1(`M$gzhasbu(i_=hwwzhQdR#e3HjU>bxjI4!zh& zbhjTA^RfKH{^buF5mDYIn^x3w%ZEGEOZybi7riox?~i&$dPWr5YQJwYo1WeW7mr$N zGFS44r~OP4BUWH7AYaq7mS4~q14l}3VyA5!x~oo_ax4Y~YNYA5-8FhgC;a+wqmOY+ zc9=Z1@Wb6gRC50z!{#`eGP5_tD;8~!WKc%kL8J|k3XEkWBZkdhS|XVnpKfv-%5&?^ z&KAD>;yTe-xh+J*nG65Ya@7OZ`VeuL+kW`-I(R;~z;(jX`IFFdd8+Fyd72+6g|%-5 z5;;Y#C)y&M7mX`R+G3MUjoRXa{hf_koT;)+sj{M?=vO(Z%=8Q$C>?K5lv#f<|l3U-!7hLRsZq&ak&$B0{^jp9tpMsZb5wD=QI7EVYR2yLp1pJ8{&*-H``8w zuA^_*UXCAnT|Hp2w_Q7eoveTyLO>zm;A>jWGljPTR3?*st-A0O*jmvl9y+u-Vg}ae zJdpf~89E!m9i+4F1h3trUwKiq;(C5rv2f%!)Lh>g>(5Z_9<0y7)%W34rSdJE3h?7{ zFV4!~5PtTjZrteoaf#KGO~k7{oBuDC+PcS(uqd&D*UZ$@enQNmLgpH|gulUV{7fdje-hMt7#0etT;FF|*L7 zc>%`R|3!;y7OQ^$?^CAb0)4_BRx8IpX@~S|{BkYKrQJiZj(H}`$}UDK6=bE#GtgZ^ z6|oT$1&yuo-JsR()$DOzC6~{It(&2Fr#h837)JY_QI9l_zJ5~Q@Vd|9E0@5}o;OnZ zB}J@z4pn&W9JyWpYdZY zmRhW;rx&DKxT}kDRE0|@limg6teU;cSFZ(s?@QHI%m~d=wNG1-^6T7-l-eMd#(Jmr zI?KM<@Cl(#yaEtLfIJ#nnA=Q#878pF;4UO)F<9%WYg5W@IM1sja!+^nRbG+$17{AY? zr90EjtXZRWE+Mj8U1@2xnBs_Od?MTFtIau!#3cOm90VK+#Pp!sp2q3zonR3w8iy1x%)`?KtAS2Mr;V6)Q#l+_ zaq5+E3k32tP6d_^4=5n?VJNP-gIgf@iX0t;#EF7AgLYgbXvGHGeIYu($Z>6ydd5aq z!r*h7Dxb7lMv~`^($-w3o&m{nwW@{Jr>f)y&a$7UU9DF=6uq#u_qPDv*?1C26ecG| zP@|j)mQcyYh={^=OWdODon1Cp_E@D&!C&w1FMMQKWU)dj9G-&f`vG%;3U*qiohmhG zSA<7dT}ttVWOqnrX9Tk^f*jSq6-HaDBK90T;Sw!XNpe%f= zl-7?)_XA~q)8CyenRf}Kx$D@jMM^5CJL)~TM4OdmtXVo=kE@z*`b{KiW~y~PwC?64 zwAvz5`>0CP+dEgDw-vQ$bLBs*I+u0cxF=5+grS4dL#+1SP$P1+urt@YBv!Q=W7ROd zd&vjla2g8hY&L%$BqFNMtF#3aT;VY{SYW_Sm%iIS5W>__agHH4E;@2L(b;ul8QQWb zOLHF!B)I%J_eQ>jv>cy`R8T$SyEnn*l`;ISdwmol+*YN8?$Zs{ zL?6c+e;$wacRwcb`MXWmvN|nV?iRFU%B6Y_p7XIk2aEuK-#E@Bl=#bq{1Y&! z3cvZ*Er{|(R-dP$;*pzbuhizc`NEWg#okYX@mx5@dt0t* z?1Q+BQ8u!sXp;%pT3R>#tQ^%9aoD{%1s4`dRDm(L*I;0bd(pH9Tpmgu$jpw4+}M%! z(*H8c-hV3^d9rWTa}mr}-MTE7Oo$HCAg?JLSlIh1MPnHtR}A zH#B6=&kWq}FkZEy5e?a|Fbuv04zt|B%`G?h1XnUaBz??m$hz3`!P!ou4VT|G6mIk$ zT|3DYN{pIbauQRe9y@OIZ(n-xca{?8v``x;ylAGArv|-~;|-guHD^jm3#tt745oAFNjaeg8@@sCdpso*M~EL@Sk{baW6B zneT&QHt;wne}~?aSR7!D5f>7R*PfvP7xt&ei#AvkZU) zz{XssLanP$sY$y>U$rXV9Y4Py=hh7OU^{=upl#12vdPT5jCOf5+c4Q>?k<|U**6Z8 zAxkVhfwl2tO|p%1AJ0^*lI805tyrK2Sc22Ny++7Zbp zKA{07=XD)4h41EgYy{KNm}A|q2*@n_uAuGb)CIr#^wQVrZ0YMS4pi-}8!H|gtiDy~ zw_wmKuQJqKIE1C_T(Gpj4{C9eo4Tm{wCUZM+y_&F$cG2?g4op&XFDBiZ@3GHO9jxe zrK^PXjcQ}`8f8o~H|De8Yf0W`p$*eqAEQ|hJZMGa41$}pz=cMC#|l)>V7B^eu|!@M zddjqeY4!%K_#bujh(&&O5%wCpXY3l2YUJI@D0_HNM%Ve1iRirK*oE|I{jMz@_oS+yiqmoxe9KDE>Kx~dE#{-wnRPp;Ru;3ld;OE{ z0VaiSBeL=_2BL=u*rc?k6*MttA*yV^e^Y;=u{c#>ui zRG#?X+Bx@VsQQ18D@iJa2)T@mTnAkY9+#P-p&=xp5-PcjX$(cg7^Nf%p(xTkE`wad zFv1wu@yLBJ^UK{x?wXO&FhgfQzqQUf&syiNv({Pb{B!nNd+q)AcYnX@^ZvX)ulFvK zxY!?eSQtd)HVZCnr^6I*w>gmTzV;*`9=t7sWE&9u!-=j^_bzS1%e=12L~2R8n@&?W zu9$Y?RL}iuYik9iW2uXo$cOINFiiK5(D_)1PZ6}zX4HdTNA)RNX7rBuC1p`j<;Se2 z38S(ba9(*MtE+iybGYH~yB#vf+qs8~XiwBS*d!zIQ^8 z2M{(Q>e~d|g{>e!iTSc~?Kz^)vHId|Zt6_E zzcq?f&|>G?3_l9_*3i$beRF5c$4^f;&-&W&TJpPnVZQNK60go$$g3sDUr=%R@jrtd zLhXLKp$yYg?U_GcX1-Zxrge?`pr9fLwlm?W%-IvG)n)0cDCSGSbzlSC7YOk-XGSMw zaDluN>|C0%`q0zcKU^3%Ew9LkU@#)Vtdkg*wUMKGFg50q=x$0zd}fHQVUjSrPasvk zafj_1XRz(n0;*y*yX(8q%e)#YDb;L`tb-e=wCR~Nt)miuEq*DHlQ<_7zfK)d=C^J6 zHo+C#anf`edBzm8cf(!O;fxvxL%pgv+W!55h5CYpc=uWPP8C-#1mbFhpG{`C0dvCn zj|i=gUrpWcH9x+E=?jj`O&P#Hg$j-U-f}cdrN@`pH3brD!kr5o)QVQfmL)`+LKYzrH~2|>X>K(z6J0%cNkaoa})i1EESXg%1%Lm_nb77u$F_;1V-CgI7v4A z>`?G;p-0Nn?7l}%mt$0u*Tq#wx|DC=Z)n`Ml58M0A#OSY%tVmeP^s4N8O_(I zq!nwJ4=c6#!-}A5JotiX=b7r6K@p5ZLIg*q06L;&5coP2!F`Tr?8r3(Pa&R0EA{YP zmfOjmB7f;WpFD(z;_04}cr(BL`H^4k4$mFb-yFNq9@zIM zit*|ys*GlM+$isi-R%3*oBz6b2)(eS)xOr*8DbBQc){NuUb}nZ2g14s{uzq3!!9Os zur$xnT-II!KVZI5ns-fv*w>0qI)YM`z)D=Zetq$JWMkQ#42+xI@P&G}F1KVq-I$8_ z@f@VL_?b$n)b0lFnFaizZ&V^~}Yob+_0a<$cgj)MQn9>k5dK z(2=lC5Dv3o(p>MERiP_TM;8ivz-_PUtD787Rw1Mv(=<+i5gPl2mKR+}o9bm-f*~i5 zbh~m5%<3qHZ|_&HM~GmqR&wpQl>qO)BbHr)|MF=JptYNI&2RvIicw2*t6{5u_S^dL zxGevH?|r3r$$BK#YFnyrko3*9TljNa!q?at*MEcLy)k@8oJK7#qI@2yAa3ixSm6>5 zp})*N;~&VL9g8;fW__-CoE|uG2Luh@mT|-trx$RUQHaX zqjnF-fIx}3T|gW8kV*V{p#@~PHLVi&_DRgxx7`;+le|1ZJpRTXJZKJ&>c~l2hI8&; z-2Q+I{b<8BL*pV~=+F!XG zE%(o;kyc76wrvsaV~1`Z`PJNG5|wB&_ZVl(LZpS2K^^v~RV)}4rJ8i81$C(<9CjxO6c*{FuTALbBNj&%A|%plr?y6W87UD|Thn#y5F0)sIzxV`wp|KS z5**{{=st_R&&oRB!to)q4s@M)z_=drl;QGQ$hYcwB`v_4;28TW1$!Yk?kA8gEziOC zz>af7ImvMJ4=tZ}ya41;Lc)e&>nqs2P$3fCkB*-Cp5E-rD z=601KcPnPM^?uEkqQdDbK+5is$ob%eIxc^l0x2mSF6L^x0B$(bp>cJUUs|L7?;*k8 zgM$BdZ$6ocU9o2U7LvtJhKnx$x0*Azd~OsN!!s-!eB*rLPQgWi=|vEV^d|3daXI40xIXC)E76 zkpB;fo&NC9Hi8%~#7=hSsD3j+_C}bs$7i~#9|Rhx8qZ?v>{CW6kNc8qpObXK%gNrC zx`m+Y2qUz5)*!NA92&fX;xZd^?0?l&+2z+8EHgU=(eavjCv)OUix)ROfaqmd?ou&Y zM?v!g3%C&RI$=XceM4{WE)+^ZHVGSKBHa^o=3f}ND&+4BL>;cq@@>s{kTgTxJHDuV zij1c0WhnWNAggPyoH3~J0~zg;QA~Szau;$L?aIO?W%JzFzpiqWa{-K#()hyfSIXD< z5nNqj_PxT>mnRr1cQPzYRCo9E?+z=_4|P%7%u23DJR~nq@;wCx+hsV@kiwllWqo>s zRo&hENyjRSet#3Q))HBUlSnP;t@1%9t3=cIbzpo~>oJY+3L~?OVjF~mh?X;y2&(}f zCRWyA@bb=vj?a+ZRjE&Bl^&;A$UH}MDJpB1KZfb^Y*}wMs?qGOc9;@JsV`Qqmr&eF zT)Z@aB#u2!75UQE)~deG6eRwHp3dEloQDJl|HCE!Pq+M^PawsmFujq?@txkNpT+_{ zY2mzw+Ta|Wt8fiA+ta`y-J?amQ;MB&yK<#h@}sQ0PC>5r>A;lN?inyCHdP?y)=}5d z>JS4LgTZlb12{)Ee2P)!^>N`CrXVJx!-Ki9&JZVF^C#!8OodJ9$iD3#J@32b1oLnU z9VaTO*cT6J$f_OrEi^x{8uDP1jK$Vy<_(Y2 zvQn+TZwdA{9D>g%P_fBejO}P*sA}<*eea>l$y#LEG+&$_&O6ESpIx{ zFlfq?XN-1^k{oTudT|a+d6TZ|4n2KVPq1TnqYG7yo4pfv#N_X15$*T&!p6#UW_*`t zr}RfGtp?YD5R-=lATVBFf-|!tZ`Ki0ZYU@`BI4%dpUo4{r&lcw$t6CGcMOmL!&Bc_i!3fhVt7bfPIW{1{aGsSgxcw_f(_z!l>kj7*7Ks;Kr+mWgR4`| z^IVe(P!H15Auq#TL6+3+cD=fj=rQ$pqvsY7An=}^eB!ZZ&`h75Zb`>MrzujqFt_qi zX&=!jU%9!vK~}p9J2f76s&S#8ZT(urCeuzOl|wF5->gKS1AYV%ONis0zlEmk>TzYj zW)t4O_FHJwsMSM`qp>`|ih~(UjDHQYJ{u2`T-qxpdGtHe?t|&3?%YI4DfaZp+k6XH zSYGlyqdVa-zG37R6W6xG`~AaZ6G@{=z6eICcLntm<*{!Kg`8}Il&+}Np8wKKyr2?& zuCS1(2^B>JFO}PwEcMsfpL%~|c$m8m+Y4wrlc;86n9?G(LZg`)?5E|tvyaP^`71VkufAZ8`8W>@2qf&gZ~pG9TP0BJ<@7X zO|Gn{j8h?#qfaSbviegyNA(dv=>k7TkV)t0OxHZ&yl66^#^5zkF5z>rA$-@$pQ#eY#ZW{!?-%~Xw zI-kVTD1QVq!5Q-;?!vK?@?lXlye|5yT4kpFP5ak?c|X?fllHB(J?fj8kf z+{st~sm1}FZF5~zEM8_l9R92Up4X;QIcQL9yu7>uR@@UW_gYnRUC5>9ysn*b8I|#Z zxJB{IB56l8q@*rv`)`fr4Cpb31YBUbf8zUgK)lUt0fZeu7GV&}fuvSIDw{E(yC{2+ zt~R8?fJQtItWTbpM8wIr21*LiknpqD1t6Xym?1*(Z>yNUXnkN%CawJA&1 z4RxeM)cVJ~xa<_$hjK^bGPKlQR=%>S*P2T92mvx--Jt}f26uKxyYfNiT@qYMpv94{ z+er%%H|HwB=6i01+%j+)&0|- r_c +\end{eqnarray*} + + +\end{document} diff --git a/doc/pair_lj_cubic.html b/doc/pair_lj_cubic.html index 2fb0464c19..3df6c4f972 100644 --- a/doc/pair_lj_cubic.html +++ b/doc/pair_lj_cubic.html @@ -23,19 +23,29 @@ pair_coeff * * 1.0 0.8908987

Description:

The lj/cubic style computes a truncated LJ interaction potential whose -energy and force are continuous everywhere. This is -achieved by replacing the LJ function outside the inflection point with -a cubic function of distance, so that both the energy and force are -continuous at the inflection point, and go to zero at the -cutoff distance. The LJ potential inside the inflection point is -unchanged. The location of the inflection point rs is defined -by the LJ diameter, rs/sigma = (26/7)^1/6. The cutoff distance -is defined by rcut/rs = 67/48. +energy and force are continuous everywhere. +Inside the inflection point the interaction is identical to the +standard 12/6 Lennard-Jones potential. +The LJ function outside the inflection point is replaced +with a cubic function of distance. The energy, force and second +derivative are continuous at the inflection point. +The cubic coefficient A3 is chosen so +that both energy and force go to zero at the cutoff distance. +Outside the cutoff distance the energy and force are zero.

-

This potential is commonly used to study the shock compression +

+
+

The location of the inflection point rs is defined +by the LJ diameter, rs/sigma = (26/7)^1/6. The cutoff distance +is defined by rc/rs = 67/48. The analytic expression for the +the cubic coefficient +A3*rmin^3/epsilon = 27.93357 is given in the paper +Holian and Ravelo (Holian). +

+

This potential is commonly used to study the mechanical behavior of FCC solids, as in the paper by Holian and Ravelo (Holian).

-

The following coefficients must be defined for each pair of atoms +

The following coefficients must be defined for each pair of atom types via the pair_coeff command as in the example above, or in the data file or restart files read by the read_data or read_restart @@ -46,8 +56,8 @@ commands, or by mixing as described below:

Note that sigma is defined in the LJ formula as the zero-crossing distance for the potential, not as the energy minimum, which -is located at 2^(1/6)*sigma. In the above example, sigma = 0.8908987, -so the energy minimum is located at r = 1. +is located at rmin = 2^(1/6)*sigma. In the above example, sigma = 0.8908987, +so rmin = 1.


diff --git a/doc/pair_lj_cubic.txt b/doc/pair_lj_cubic.txt index f0a864c609..f1dcac30be 100644 --- a/doc/pair_lj_cubic.txt +++ b/doc/pair_lj_cubic.txt @@ -20,19 +20,29 @@ pair_coeff * * 1.0 0.8908987 :pre [Description:] The {lj/cubic} style computes a truncated LJ interaction potential whose -energy and force are continuous everywhere. This is -achieved by replacing the LJ function outside the inflection point with -a cubic function of distance, so that both the energy and force are -continuous at the inflection point, and go to zero at the -cutoff distance. The LJ potential inside the inflection point is -unchanged. The location of the inflection point rs is defined -by the LJ diameter, rs/sigma = (26/7)^1/6. The cutoff distance -is defined by rcut/rs = 67/48. +energy and force are continuous everywhere. +Inside the inflection point the interaction is identical to the +standard 12/6 "Lennard-Jones"_pair_lj.html potential. +The LJ function outside the inflection point is replaced +with a cubic function of distance. The energy, force and second +derivative are continuous at the inflection point. +The cubic coefficient A3 is chosen so +that both energy and force go to zero at the cutoff distance. +Outside the cutoff distance the energy and force are zero. -This potential is commonly used to study the shock compression +:c,image(Eqs/pair_lj_cubic.jpg) + +The location of the inflection point rs is defined +by the LJ diameter, rs/sigma = (26/7)^1/6. The cutoff distance +is defined by rc/rs = 67/48. The analytic expression for the +the cubic coefficient +A3*rmin^3/epsilon = 27.93357 is given in the paper +Holian and Ravelo "(Holian)"_#Holian. + +This potential is commonly used to study the mechanical behavior of FCC solids, as in the paper by Holian and Ravelo "(Holian)"_#Holian. -The following coefficients must be defined for each pair of atoms +The following coefficients must be defined for each pair of atom types via the "pair_coeff"_pair_coeff.html command as in the example above, or in the data file or restart files read by the "read_data"_read_data.html or "read_restart"_read_restart.html @@ -43,8 +53,8 @@ sigma (distance units) :ul Note that sigma is defined in the LJ formula as the zero-crossing distance for the potential, not as the energy minimum, which -is located at 2^(1/6)*sigma. In the above example, sigma = 0.8908987, -so the energy minimum is located at r = 1. +is located at rmin = 2^(1/6)*sigma. In the above example, sigma = 0.8908987, +so rmin = 1. :line