From 1e51172f944a0d314cb1c2b53f1363b25d022187 Mon Sep 17 00:00:00 2001
From: sjplimp
Date: Fri, 23 Aug 2013 21:01:41 +0000
Subject: [PATCH] git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@10686
f3b2605a-c512-4ea7-a41b-209d697bcdaa
---
doc/Eqs/fix_wall_lj1043.jpg | Bin 0 -> 10996 bytes
doc/Eqs/fix_wall_lj1043.tex | 12 ++++++++++++
doc/fix_wall.html | 28 +++++++++++++++++++++-------
doc/fix_wall.txt | 26 +++++++++++++++++++-------
4 files changed, 52 insertions(+), 14 deletions(-)
create mode 100644 doc/Eqs/fix_wall_lj1043.jpg
create mode 100644 doc/Eqs/fix_wall_lj1043.tex
diff --git a/doc/Eqs/fix_wall_lj1043.jpg b/doc/Eqs/fix_wall_lj1043.jpg
new file mode 100644
index 0000000000000000000000000000000000000000..303743bf3e84dd0631cc4f1976aea6eaeaf90a26
GIT binary patch
literal 10996
zcmeHtXH-WNlr;gMZ-i(b)SKnl9G;{j)9q(
zm6er>mV=9fg^P)WmF1=pKuiemyfcA=LkYNP0Wbk@0JyhJ+&1#Bxcw|1?j0P$+tT+G
z0319V+&j2;hzRZy;FAJwYu~}e!zZA;2V|!bl;uz-yiZM|ZS3%ML`eLFlW!TPsMx!Y
z&XK95n-z@F~hukAv3i3xI5*Gouv+r4mvl+<2QQ`D{V^c65_vgwy@4A{q3IT9MD3j8}
z9W^u4|6K<)Yj@v@umB%C6b2QPk3Bguq8+RG=0^`9#N2DIrmSW|b-923Rkc1|C*LEM
zs~c1%DjC0$h;7twX1lW-Sj{gw=|*ZuvsX|1iS9|1iWtSJuU6tb^${0765z&=Yc!Gu^nJrINnH{ZEkcj?e(i
z>`S;;E87s$z-ojnqaRp%gM&kW2jO8@v$~fV>ofOBWqHlhSG)MxICb<;yo!>Pmf+VC
zE&e8#B@ee3)`}-@mYVt)aT&u+kOHC76?x9940W}(7>n-WUpCNuy8iIiWOn2ErPeEP
zZ>sfi7Ys$gqraWaySUVXmcvt=pv?H&6!sMX$!YMIWy@KmWO+XBwpD)iFiuM@HNLa-uArxk1#>
ziQNh#jcNNAIcLd91{-2o39t92KUXefiTH_W`nFr|Gh<`wnzjWT4iFvOA$sLrzsB9e
zpD$QjW94Y5I5RU6=yn^J7sd3E5A2~2V3@^e5;BA=dr49t8|wq;hV<*P?`>t;uSEJX
zL!?T-q2cYrT&__#>u2qIJ{J{*R)=_y1i1vJAwo`jgb4RcXKfYcFr7cF)O^v3k%`)wgv-R~?oz8?&IRa_iN9D(zooh|HP7wXtaX
zhO$h}lQrLK&QvUe%-(0n&@Y$AAK3Hqt;dXtw7LutdD0Sze@OWOVmwUk_ZY-+(kFM+
z92axF@eXZ{HXp6lTUZ@t+g(lW6LFce_&g`Tp6#%A7;#H^^RK3
z0}F*;c4P<~^P^8#Q{<{L9;CoIaOImgrS@G~8=)B$CEujO3vF?R*6WzjnyuzDgiaPInHx|@{^j&r$;PLxC5XgXRrZ;LvGYC5=EO8BV=8{3
z=FDQTjQgQ$Nzx4<0UN-6yhZ!spKSVM#n!^VY=nn_fj-08uwG_`2uPZ~clG=?Zk;|6
zeBRBG@w^cJ^Pn!Q-72nIsc)GxM_A1sChnGh2izCT-M%v`e|e&Q%VoM|;c4#WFHWEK
zqC8zB>y4vq87nW&s`F5Cw4oh~8Lg%c$NbP`Bh>3kcAH3WdS@mIk8;A9M{L83JxM8w
z52!Cj`uN3&@CJ{EYb&R*J@T=_x)K6ns9;Xud#+&oV`Hj#b|~Imi(qdeUF01((ty$T
zda<;}IlotB82;6&N3=u(oL8c$znPi3n(4Zsz!#suL?|#(5w=XY!l{Hu8fu<17?w2_
z#+<JZ$(Lbb_`Kix_(5O`lo5Ek(w{-L9IaXdxb3*M|s6K9xA-
zM|ry;ET%koy1gks`H2e@6lS9bIZ}*f9xuMg+mQb2rR@P>L($H
z!otpO_k+FZBb~3W@Hd7!JQ6w~i{3^&nSbk9+$#D8(DSewt7LI%roSiC@=+#T;hLh9
zEY}H|k;93aV*<^IIvHd|kjOfpN-
zGI7yJ^T;liZ3c0g+F$}_$Iw5d>U?ZK!c}Hgy5X-LIA`|5_z?W92_$vv^ifvA%Uakvo1T#U(k12zJJMpvE2arn~mfV5!UQwPF+b-x@R?oF>M^M4_lPtY~yNlHI$IO(XVM$99X6LM{96|=2%*5$NA6kFOx{$1Xb_r+?i<<=UsCa}X_1zvayF_2$W-8RX`
zfKK}tC2|O?vBaniZ!*8$s4pGiexs32rW=8$k5+*HKToDlT1~0h)VV{V~Xw;6-t<
zM+s|nePIP|i^CJq6OVP=M~Ba=&%~Ta5}QL@|9HeNGdMfXtz%EPaCfwsVp5x-Pb1A8
z#r0IbHYHPKLE}!*m-a@XD08`LobQg@mD4AN)o&{FTeT;0sLv+a_c&IFavMB-wbFR_RJjYwk`!OPK>j=wlKH;7Xg!6fTUAWhz5!-A-qhNHa=p
zK_!2?W{S>0GOWF-Vfk!GMRI#rO5p@QFO8$iUH(ule@4Gac+xkX*GnmE?kIoQPTXt!
zg?EW{$6sQ-MK;!b7-8KPVf{w*`I6hnj7g;33Nb|_uO6#hZx|yMHka;Tn(`xO6=5Ie
zc>de53d-hB^>eZ((0`N3i^+SBCwR3WJPWd&_if%rL!ZydZ|HesFOyv0!AiSvw`~VZ
zmP^ET2%XtJ9Q%MwUWeiGt>9oYXOJ|VGp0J6O|==ZR%JHy&yEC>XG?SU@X0L@#7ire
zHA-xW?o8dBo(DY=Nxl!VL2hq*sg)PPeq+1hE*xhX
ztnYp^wtH-z@Vrg1{#4=wcUw~o20{;U4oW;OwT~|8nt512aM{fy8)cW>_Pcr6a1^1r
zHCsL$r=KBlm3adQ?i5+fca-LKm2Xb}SP&4$==wQ~pY#x1r$69XUZ4AWLz({BL>{?k
zsDUqJOX%hL9LZP)R~px(akZ+MIfibCxJkHazA_^TN(}v4Ux_(!tdD`(oIyR>^~^9Ko1?De6(v
zXS;eIZvgL2ZT?8Lzb&s9z}JG384{II$br|&FEuL#@FVRe*Lr8B#F&DT`e96Dpr0DB
zd{OJNy8|i;gdB^ypwL+o-jHzkJ#eILdx;&@GOhF-^F?=b;GkV>%ejpFS%l)!_G>Re
zCdZmhnWi?i%Px@qg58=I#gqw%K4c=7?uxX=wUbAIa?Q1?WBCT4-(Fn(Vy8KrjcJjP
zu>eF)wfTCi832`%cEuEG}Pg?ImqR{g3{a>mm5ay+Ko3kbe&ykNl3D38*?O
z5ouS#3V!?finHnyS8M)Xa1%m?#8-+XE!Seec|H&UagLAe`MOVEe-pUO9JJl(TIR%1WIu@k`
z55DC0T9ea;lYceI?g!eQF73N=%z3t5p}iER;g{`BU$nIUJS<7pD0i5*1Gn3r+^Oj_
zG`RQq_(Cv^r?X0;e479S-Qa1vTb1?%R_>g}&*k(%SA=T$$s*GGM0qSGs(&1fcW^AD
zHI0&FIBgQO(J}qO+|8t`QaEwW!o`9O1dmpj&PZJY7x+eXW}TgaC6v-}M+<_TQ+6R;
zKRy1vzeVm0YGBAQlFCV33lsbpayp(+`BdR4(6OCH1Lm(5xwx7XZbI&4qbvTVuHln~
z_KMX^K{$};BJoMb61LxZOkJxi0VW*T*$vNn>%`AO`(-~zl?0giw(qdTizA^Y7W$F`_9m4l=%K@xQ|V
zKy|U>RA=}9|KtCzPN;jo4(ry!g~sH^wJNfUv}pSu+e&E!!XWz2cTme)Qho)CYu(c$
z5*~boX~%azP5!z!*H&*GXEF&5Zx(=64?hUMqs1IJ-H<&0IeoKXSsHn48$KiX;)wf@
z{l^U;&vo};FfW0F<58V1j1@IpRYtD`9BF&8c>s6DZ6Z+pPC
zWB6FWG!7vyS$xk%H)4VxbRDhSIM?IAvXai*4@wt)Ydv|)JD@dA9Xcy&3jCS=ZLA^t!CDe+n*;gzNOZ@So(ItJy5bt~{Ef_#K4`1!ZLnC6
z!b{ZB_6O*Au2f$io1xU_D|Qeh^It$yJI|#XK>0E{cW#d_ZhOX7Wupg>ZzZ-j0G>pzMy-Ib95NYZ
z$Y=YPqG>jQelF|<0P-E;8A?hTy$Ivq9^Csbv&t}QPf<~V?PZ&hxHsabcK1uYJ5LIf
z>`f68iMF}Nj`R!_=B4u!y1c2G%he({4Rw7mh1zNyj8DBoJVs1)>!{@XwR8TsrANuD
zM~H+WMY0KZ_&DiRl@*Io4t2~GYr3<%0U|%CCcS?H`7q=ANmPN+)L^yePwtd5vz))7
zLcs>JtT~@j_|Tc6amr(}6o%4sSwW5LaljuI96WTpvtUhFKj(5g&u!Hwv$St@`>Yc`
z>zh8k3JXT&PF3zHXH|;vqxeu5=$u9MZcUKoL`hf2Uq`R($`<@8&vwOe50-dM%!COY
zeWaOc=pP65Ia7;lb03n5ech>eQn{#Cd@M#P?$d2jB;C4gmn}T_jHAi$N<^1Yb^egx
z?qE&C1Zs!Vf@ON^l`pm3JBw|rs~ySHMW1;A@Ymj$9j-j9ohcbAL<#Bzq*1=Lhz}Uy
zwN=XTWF{1wZ(iX~92-oa-sLAKwuG@0rqjPT@KdXVjDt~s{V5w_agXHe&d$%nZoIV;
zo4Su5YeqSbL&{fXz>ZRw_-9&6-Q=M{{hi;MUxp%>dRt_N
zzR2s<_ixExiEH<|n+awdO|8&sfd{)jH1cpP
zX|@pU1&1hdQ6MR_Gv1i0S3Zgx6Z8=+u$&MfX9J==Ip+Ho9O$P$-0NAgo-1Tx8B?S4
zrX}YV_wKYm;xqx%y06(XsB|=NO*ysxVYeYK6EShmwWr@cS;=l&<6vZ#8W4Z@cB}0+Y&sA;dR;QR;`Zb5286EZ97dl&D#ON
zj~b6)E{K}*WK|b$z5uhSzcYG^|0+MPxlN5=P3X4{$sZaMtYoiu7}hUsy8(=m2Bp@<
zV^T<50L%eJIhrzl+aYR0&tiO64&;Tkme=PyK~;;qv1d=4H^!>#(lLz1TX)7{oJ2kd
z{i@=Bb4kYxU3g`Dc_~^NS@j7n9Y*z%q4x%mLIwQ(yWT9jDY{L!U3X<(H)!xRgYyFv
zqlDdY*zc}q;-6YzyoNQS^ogDxpAGbqDk%29S0{uyK4}>~X2tKrrekPo@YJ>aB-p}A
zS{dc&+QR%H1dTZBV!AUOs#O~+{G7?S^59``%e`)T?c$08mVm~fKCTPtw^CDDs9PQ6
z)aeh8A3~ff)
zqT7pSGs~;shP0Dv8gy1QTllN*$A`>%e9=R{-U5~RbNYr+M{sbGe0UIC_5?m4vvf)G
zO#}Z+&i9h5Vr_*X+Rnss(WK
zbn$j=Je8!P+9jJ?_2UBitaU3D?vP>Ew*9QI+^#`5s9MBE->cMx>mus?FAk2PV`7cw
z43XMN^rnhgegYk&BOTfQ*v!T#qOE$1xC(?j@nki1Y|`^VQF(GwU%MQ;hK35C4)1W#
zfD**{sf1%#zqgnKT=I@YQWelCHGea8=nAdz7KZjjixL<*=q{h$&jXKcmLtCm)~nPV
zR9_TOQ!1Ds?jN`u9Obzt^K&d&O)bv^1g>)=I-H`Y${HFR)-yDrc?{n!WPVQhzk4xo
zdEc_*qDNXS5hW;seYKK)Qlvs1n3wj;6V*DAC?;jgueBw$+ouw;`mKEgy|T+a5|wxV
zc-4JPg>|Gu@TiYM%S&F$5|#vi&4-U^NvsZf1C0y!W!&mBI432hcP~_9)P5ju-v186
zn@q&Gd;(7GwZ8~t%RNE7ee?Og#j|ynr_A9{kY}iDye366iVdSr)_AyYud-(NZrnc+RR<7#nq>r6N
zv2A5~Kom^O&dmlLMhdW7J!dZ7U0vfRKA(OjprBEcBcSS@Kc?0TQxyWSBQxjWwVV+Tm3H02eBhW2KsRf`#L?-j4+DE{|``vy77l
z?#$HWxs*(0_`@X!rLNM7a{$>;&;Im}gD5#(>qP51K#upziH8yq{f
zWg8msMp$ToJzmO2<_R=ox`ua*!}11D(X;5e(Bei79j3QF$BKsZjWB9=Y6N2pPEWXu
zTz5AoK6Ct^)5`zIsnK}B_kt
zKzxEU^Qmh-oSagRG9VqFfm(;tua<-q%D4mYwi{F3tZ+0h`6?`tM^~Gvv#nA3yVJxw
zI$u5ZuJv1VAl3JxW5gj-_o@?$U8mUYQIxXGV~7oC3sU45w3s(L1#w^~6*^KVqfG
z?tZsgvQk*k_ANwIi}&}oM7eP7C_QkR^F7J#hk!IjJW9}E$z;-8reSFD2u_DT)tnoK
zbWHLf{2+ho;QXmkHatO}D5YC}iE5^gD@tF9D@uv0aQQY@tS(U}qfDRJVfIMcmEdK-
zWbQaQ&M>*8zw?Qy{e*-cZY7={vPij_S7k;vtRM1o(ov^lEjV$dZ_&NVYv0~7N$u5WAR^0%kDi%0_0ifNrv3@PI%A9dcbUY}#;^c*{Mc9pP?
zqxVZvCKPB9giEZyv@L>DG7(YHlhI!kT44f4w-V3Jv}aqftJjk{RPy_-`^;eTrO3IG
zIuYMHQnnwwP$Gk@8E*U;OsHdawT~&)tsF1i>uec=bVZ!s0Dw&;N@ek_EkJLn-Binb
zykB-CS2pO7bK9+6ka_jIVOz(+0^M|H0Wgqd;gU)qp;P_W#&J2yj++h;LT&bQqVM{<
zGKjhQ2H>P#&XU^Qp8Lz^#M0o^PrHYfy^_RXG#ydK7
zYV8$4=#c107W}N2YumfC+;ugqRZh>(UP~l@mwkU56ICR)$UOp@=-4$X&Ye7-t*CT4
zg6RoW2{2E&pM?ia9GumDxADBpO-?}=?19oC^9zSK{`M|Ag7efyXJ9esBOW=ic7AWR
z;pxcT_9i>lE6@|I&+U30SC5`mPFHfHzB(1xQyIq(bgU-yPg@iR5IEk0i|;s}@{lV#
zk-bV%9X9|qz>&wr0y!{7^Oi44wZwu5LB+$FD;7dBuIlozPIgcI8$e$QEmjeI?Vek|
z!CUc$443_TheYL1>GBv0Hn^)`X1no*^IhRAx*FF4=GRiuuFG48HIXbxNK%{bbLsih
zm~aFHn07V%tRI+N!S9U|lAy#l|ACmFRyaiITOmT3{13TVvo~MUv}deX(Mt5_l!aP(
z`R~cOnCjs<*Uc>syAcmJcGnxg0?3Yku+EO(i&k|AbS8GJ^-!`A`%r9{YkG@6><3ma
zFIk&5v|E;mUf&6{|BA;VEXK&Dz%&e0F)zWKNN(@(yHD7)c_hirujBPDv7*%*pBlU(
zMvtsWOWvTt-P?>lC8Dlj*ln`LN9pUP!(Iu*WXXH1DbE2XU}h0wdYk`k;%VR3UPDe=l+hZosWHF30o;{H
z;^v!M9c$aJ!T-_ohQ#NllBAcM37?fq+g7?GkXv_)gALM*y(?
zo^A7otK~^;n=;AI;8=*HQch*qdV#i36P!lkmyj#E&?c1MxesrL2_k!<>~%yFW}et#
zAslN*665_=$Ml>hET9QmD`p@3XlVumg@OS3OPDCt4Tjl((b8iMJ_nYf76mwoO?5{e
z6V(i@FIRpn
z5=}Hc;VCPN?u6$xmFG;7ezQTZrJlp$!v4WWi=yasnt6KeJPTuC4Hm+IStj%;
zVGQn>THOm0(qe~!KoOv5mhsMXrz)&k)dgau0w$I${Wkc$epis`ZN>H{RzwOYPNDV{
zH`>4q$q_cxby5=r*05mItjRudN1pA?QQGV5o@RJP`r8t=~4Fg4H3
z=1Yb>Hc^4G6#q7b4R$fPbXH0WznuuS#QP$TLVkOYh}ZcDRS({zgyF662rz({W`2s-!?o!
z!8Y)Q$Tk&7b>tPqa5X|)+TV#eRS2sXv(0Os{^vxT5`|*9XAtlA
z&?3t914)lfWXI5WNwfXDd<*5yd#_m6DXZ{5t13HdA;4v;m_t|9?U$9Yey?po{0EJpUZ
zY-{j*XARE^x7o_XH8MSX=2DWezt23=huc|`gmIpBp{2q^xoe9`c|MD=)h)u+1y==e
zM>`s@18~sNKBPVeEWs(GL51dmOV1a$kV3laJz`l#bBp&4wt6L{Wm{|G>9&+Hgj?gb
z4sGf=p>3mW1HKXObBNje9+G}%vay6gq<;X2Q)xU$T*)eW-VabuE<$ys#GUTfx
zek*S5E`Q*enrXwPlg%AW*ZYNlHnik&?uGs&SG6R8YcQLq|3F<$`Di)shFefkex8d`
z!t!6b-)`8&ryoz>9980SjEU`tlfFyxX?L`YlF(_l@s!^NkB0p_dH-j39^Ficm+P5f
zrSkh_2@BPnm=_l>U%omxTS88pwv9g%t|4NKNeaXuYwpQd-
yLb?l1{d9^!mCRKAAZ?4oD|$DS*wZ>pGGss!WLqVDIN@z_`Onw>ujPlE$^QYP=y_HE
literal 0
HcmV?d00001
diff --git a/doc/Eqs/fix_wall_lj1043.tex b/doc/Eqs/fix_wall_lj1043.tex
new file mode 100644
index 0000000000..8ef981ca97
--- /dev/null
+++ b/doc/Eqs/fix_wall_lj1043.tex
@@ -0,0 +1,12 @@
+\documentstyle[12pt]{article}
+
+\begin{document}
+
+$$
+ E = 2 \pi \epsilon \left[ \frac{2}{5} \left(\frac{\sigma}{r}\right)^{10} -
+ \left(\frac{\sigma}{r}\right)^4 -
+ \frac{\sqrt(2)\sigma^3}{3\left(r+\left(0.61/\sqrt(2)\right)\sigma\right)^3}\right]
+ \qquad r < r_c
+$$
+
+\end{document}
\ No newline at end of file
diff --git a/doc/fix_wall.html b/doc/fix_wall.html
index bac3e55692..7b3762f5ab 100644
--- a/doc/fix_wall.html
+++ b/doc/fix_wall.html
@@ -13,6 +13,8 @@
fix wall/lj126 command
+fix wall/lj1043 command
+
fix wall/colloid command
fix wall/harmonic command
@@ -23,7 +25,7 @@
- ID, group-ID are documented in fix command
-
- style = wall/lj93 or wall/lj126 or wall/colloid or wall/harmonic
+
- style = wall/lj93 or wall/lj126 or wall/lj1043 or wall/colloid or wall/harmonic
- one or more face/arg pairs may be appended
@@ -78,6 +80,10 @@ wall-particle interactions depends on the style.
+For style wall/lj1043, the energy E is given by the 10/4/3 potential:
+
+
+
For style wall/colloid, the energy E is given by an integrated form
of the pair_style colloid potential:
@@ -115,16 +121,18 @@ can specify formulas with various mathematical functions, and include
box parameters and timestep and elapsed time. Thus it is easy to
specify a time-dependent wall position. See examples below.
-For the wall/lj93 and wall/lj126 styles, epsilon and sigma are
-the usual Lennard-Jones parameters, which determine the strength and
-size of the particle as it interacts with the wall. Epsilon has
-energy units. Note that this epsilon and sigma may be different
-than any epsilon or sigma values defined for a pair style that
-computes particle-particle interactions.
+
For the wall/lj93 and wall/lj126 and wall/lj1043 styles,
+epsilon and sigma are the usual Lennard-Jones parameters, which
+determine the strength and size of the particle as it interacts with
+the wall. Epsilon has energy units. Note that this epsilon and
+sigma may be different than any epsilon or sigma values defined
+for a pair style that computes particle-particle interactions.
The wall/lj93 interaction is derived by integrating over a 3d
half-lattice of Lennard-Jones 12/6 particles. The wall/lj126
interaction is effectively a harder, more repulsive wall interaction.
+The wall/lj1043 interaction is yet a different form of wall
+interaction, described in Magda et al in (Magda).
For the wall/colloid style, R is the radius of the colloid
particle, D is the distance from the surface of the colloid particle
@@ -292,4 +300,10 @@ being minimized), you MUST enable the fix_modify
The option defaults units = lattice, fld = no, and pbc = no.
+
+
+
+
+(Magda) Magda, Tirrell, Davis, J Chem Phys, 83, 1888-1901 (1985).
+