From fb1509f2e79e06b5252fa5d5021486552c593b97 Mon Sep 17 00:00:00 2001 From: sjplimp Date: Tue, 10 Jan 2012 20:05:45 +0000 Subject: [PATCH] git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@7515 f3b2605a-c512-4ea7-a41b-209d697bcdaa --- doc/Eqs/angle_dipole_potential.jpg | Bin 0 -> 2616 bytes doc/Eqs/angle_dipole_potential.tex | 9 +++ doc/Eqs/angle_dipole_torque.jpg | Bin 0 -> 5110 bytes doc/Eqs/angle_dipole_torque.tex | 9 +++ doc/Section_commands.txt | 3 +- doc/angle_dipole.html | 105 +++++++++++++++++++++++++++++ doc/angle_dipole.txt | 98 +++++++++++++++++++++++++++ 7 files changed, 223 insertions(+), 1 deletion(-) create mode 100644 doc/Eqs/angle_dipole_potential.jpg create mode 100644 doc/Eqs/angle_dipole_potential.tex create mode 100644 doc/Eqs/angle_dipole_torque.jpg create mode 100644 doc/Eqs/angle_dipole_torque.tex create mode 100644 doc/angle_dipole.html create mode 100644 doc/angle_dipole.txt diff --git a/doc/Eqs/angle_dipole_potential.jpg b/doc/Eqs/angle_dipole_potential.jpg new file mode 100644 index 0000000000000000000000000000000000000000..b1285e895749cc236fa379eb6862f8fb03a44470 GIT binary patch literal 2616 zcmbuAX*Ao38pi(-5e}lEja@7u!O&MF0Qbpa`hN_kz zrNkOaMX9A!?ZgO%MpXK z`v9l2fH?pGb8~ZpA>Rc8f$%^Dc%k1B;^*fRfC~u=!-e2*5rmYe$XSVVaJZP9n1r-6 z5{VQMl}E|TprmAwGN)YtKQC|{I0Av-z-cca1%Lo9@Q)k+1;PX7hH`;;zq`6;01!79 z7YGdE0YmtHBz~m1A#fgHIjD%b=~;xlhJ*Jx(PRZP#~b0iKG>YxGV=JR(;0yG$HvRR zabDX0#e*bmyCjE#Ve{nob5Vwgo=^TQe>r-F@N0g-el0ydFO=rJyzdy;V$^9fbe?mJ zN%~V0H${FE_M7sZfMxzg@S=_-eUSoEYzs%l<7Y5>`=@|7p3kQa)$6BxEjRTI8=0vZ zlQJx~LJY2e3MQCxP6~AxSDY`nM&wgyl`RWaM5<@*RIB4GtdQD1%`f)n6lX7VHQJ$r zgtx;^FRkxsK8mP$qKQ4$-$O#~gru3H9zv%vae8A0A%Z{*3MUiUu>F(X-z{xVVdlAg zSJQF=3U#Nd;Qraj!LNf~5qrpTHw$8kvC)!-ON-E@0V{6S)!(~q_7GZx(dOWK{Tl{e zs#!)g;AJYsQn5z<^M77hejun4y{&0Vq$={4LMMx;g4vz!)6@lTW2wa|H2tr%VCjLx zGETTkNVk@YU$rZZfuUWvsWQM%CyRx#_&sGtx2Dd4jCim@PPk&vn$AAt!=|}&BpbQN zrz#pfqrGqAJGIs8;)<=9d(G|}*6X|!p>kIyO>K4I1?PSiL(?bKT%UQ(;&qJ^pX_t7 z^c)y{y`qI??k>9-Wy~{O-yTW`XEKHxN($2(qUPP%YtpxVb*clz{W%GqTFd@!Xr^2Y zW>&Jb$f(Erg%>I@v+#KH#Ytq?$N}aQ*z#Oj$t}{9JYp=3-R+9RU^W*Lp?jx*%JcJ@ zC!fb$`@aQWfK+s=)w&nIB&f#~>F|4PS^i1q6Xd$c2e6~Ny))h$}{gXr5W`YG}k$H5WNHwh_3{)8>iQ23`6IGlpUbR>JH5&O?Vasm$6etp! zf`HeJ24((W#MI)ojDC}bnFa`EKHoCR{w0Iz@D*{f8Qa02P;yL0gFtsy>R4HL#lHu~ zjux;sPAgCE^#v$6wZq8LfrlDO2IM%?^8o?(nC8ZKxZ^V#8qKq7kv#EEAHUR|H<1;! z@w0{T#cfes5)_(q&qSUmijFc2(c#-r63X9hU#@lWYp)bXm55^la|;XxcOM%nvrhBgxpZ~3G*QGhpm7{o1lzq2&i;Ebj0nW|qX zRMuNI!vXvC+T)%N(lS1sK0CF7yGnY?sJl7QpcK-$2Z}xG)(VwMJ<`Qp&am`~wHtSp zKGfD}&lj7)zn$MjHkf4&V?Blq&}!DjC0LJMRlzSI&)9qKl)BJn7zb*%zCj=h)^y2* zxOm>;#A)X;r^fyiz&81=%;BED~U(Kva@;RKuIH5td;?D*c-wHh8@-(cwdR zUpKZ*W^#-+uy5EmX*3-ENIj*jm%T}@EeqvPLk~jD2`mqO(nY-FOS9bo*BqDGI?b#Q z7%Y4)DX*jMXf`BW_R*`k!JMI#!x5X3@J?S(zM9#y2O;7tSCUJ0@s|JaUWHJtPY}jq zbY7Sy(n(9GXJ= zt0B`J)^DdU^Y&X;3vArwd|QX?8tXtqvW@w+r+{p-&JOioH!?Zi-RxnVy7IH23gHs?XEcpC2Mq`NY|syr}GD z`YUv0VhO9xp6x87oU($mS|8Uy42D8fR=>_0sUQY26Jp zY)4F_xx7fBbCjPU06=j7b9d)?!PRn(yFVdyr}z}OgKGxnEsaz{TqYjb_i zqvAlqtKvzB^+C_?M;e{DeoiE^YQZCbuY}05{p_Bxf?1{yyB_*;w1dQ#4sTqayoS6f z=DL=&K+oqgQLu>He@1`u+82UnKlr#-@4{%YK>5PLa;`6FdMrw&Z-ximRxr`If=GT% zkW)PIez+n-9n-p}lo?UT%OgR!+2N6-l-qQr+}cp@-QMJpFi1 z9NA>?i%a7bXMPXW$+f(e66taP5#|L;3LrHYY|`lup=7DY z$rB6;*g_{UD70A)Vb}4Om_n4z^}K?Llv<{QiwfTBm_Q{Tk>JTa= zyKPQzJYh(u;T9uiYW!-65RsJL>-5gi+utorqg@T1;M$kk7{S)R)cA^~+U@+|J(VJK*inR`mKJ_a%6}S^nRtDTQ-BeJ;BJ z6U{~a(}Pyv*@zt5RGh31VhD7@qXp}^}U{-P5)=GomZ zBPfk)+>V{l!$&Pl?!J@f9fkJY8*75Z@-$9OwL6;lBBw^h^L|5@%(OrVx;oQR!FNL& R^l^JB_z(I2!G+W5e*nm(<9h%A literal 0 HcmV?d00001 diff --git a/doc/Eqs/angle_dipole_potential.tex b/doc/Eqs/angle_dipole_potential.tex new file mode 100644 index 0000000000..8949835eb7 --- /dev/null +++ b/doc/Eqs/angle_dipole_potential.tex @@ -0,0 +1,9 @@ +\documentclass[12pt]{article} + +\begin{document} + +$$ + E = K (\cos\gamma - \cos\gamma_0)^2 +$$ + +\end{document} \ No newline at end of file diff --git a/doc/Eqs/angle_dipole_torque.jpg b/doc/Eqs/angle_dipole_torque.jpg new file mode 100644 index 0000000000000000000000000000000000000000..996a9df3cb3cd51f23678121f6e4765a7aed35e0 GIT binary patch literal 5110 zcmcIoXHXMdmkuS;r3OSuXbBJ@`hp75OOS*FDM3K`E9IpI3`IZ@LFv*7T~Htq5JPW4 zs&uhXB@{sd(m_E4>Fe&S^ZnUB-_Gv2GtWFT_s(si`USv)(KpZsFfcFx z(7zYpbP9k30GXJXnVEpU0SE-1Vc}$D`EA@B9BiDt+R+%odmG{r>+*@WPl9E;E4Y_NYBViwHHg#Xm%3&XYg^u=G>9IvafzT7Z#4nF@UVV6& z=e&FI1^t&x;d|G7-%IXh7uyti#=!L=m}#4EL`qQ$7i7tyymlge^P78b@0{RrgJH^o z%mjW$MHg}lzO*z%4{w__0pP{=#KrLTt`n6Lmg^H`4_fYI*=8GrZO`Jl$m)qi5)5PM zM3ZU5-+W~d7`;w0Y$eOTCx!%EN1u;LXXwsLnt_Wo7p<1-xw&z8{3|f;XDbxFRIXM` ziE`9@&NCFq2b`;pfwiQpMqe?V9SYH8-yr>oVQgTNjr{j?)o+AopSV;$2jYRX7yq(~ zYz!diTHA{^-k-5=EtGm0)%1>}Se#w*Q7 z)BVODm}wM54rsFTzUM?62$#Sx1`BOYCMlGO0)I|X=(#@UrCiuwpuqG}6~ zpJ|ClkwI&I+L4HmYsE2x8--?T)>;!I;+YS7%KQ6F#&=uC0z#)ASU;qzPL7(8zBPIk zXp^|0cxy5Kg0H;0RY72QvR@BSKMy3I3zEJNOBH`o%?)lDC8at;H4u5-T( zykCKq-?vGYIR*Fy+P#LE=O_fzHtOgFz1ceK>Kw7od5oljZqKGbYt6K;7Pd{tzIjhv zF@5syy4}gTA%swoql}wm7#f>n9kiBzNVYzDsc+}Bde%4~WTM9MuPJBqJMmIeIjK8! z;h)Dja@l-!Sc&gaL0zg<1NO#~A5>(=cYp z&zHFpjmfhatM+DI^~Hk1E=I2+5|WAHAlcn@;v?hd%%4y*_9~}<8}Cj5I|B`;fSear zez}O5#rI9}@_7yUP#A9ONxVf=xTdS!ZCT^uZlo)6G%hdCnRiS5;*T^e-l)p+(Ry9| zGxP5{mb?dfo(G>KyY4;qUkNC>5b5h}Z|vuBWuB{-uS?J-0BxNOkkBJ@r79ViisvK$ zzL~3eUG1;cR?7s`Gs#)^Zd2Y&Jv^*G)d6oLTk6b5l#mc-UQP$`p#kQ=S&e?LJ^6W_ z_}KBX^zah?xfbDwi%ykd)=Xcdl7P*_NHWAuu8$Q8{$xN6gm_Dm6MzFEJFUh%tUYvmv4UmJWFKR)@KmQ>J{DQy?;+&b^RgY_yvC0z}K zVQ6UMf(z%1s?qNhXqOf4f8G*0EPCRuHkOt5|m>L_cq^(xv`#pu!kH)U}p)`W? z5i(^O+b6ZHo(nz!wf=WnJHUO<#e!CFc34Bj$ZKnHs;p=c=t&s+k}$Wq^!zW*r!!k& zv*OlH)nd0&Bo!T4O&QdDx}dV#5VL6?lzz*CFnU?B1Tw{wl2&b@?L# zab-RsS}meTYQqpu)VZ5jwm30cZi%TXS zp_$i^&Swuobfp7*dY1b8K)0#Rsb7DcaH8^)fN=WMZ8$EKKJpbMHDP0_8Mgak`2{$Y z?z^3BC0?s_;S>;bZ*?IH_teV3fhX!3g0#Yl@#mH_ItC)E=YrF7%lFo27`u^os@m5XQBf*VFu^6SpW@kjKA(#jz_vE3L$jl|7_r`l#rp~>ad3+YQV&R_= zr}5GY_4Gn2kB+K4nxXaHkmy(!AA?2=xq~1d4d|Zk+A46*^+^0NP)Ak8QmCK?TU;j= zV^Sd?PO(-`zL_U8P;UOC(*NXZjAKwsCWu-&pHY^ZR&Hd-j-tbQNGv*ZZpLDJOmqnv zy-soaJW1&NaBugvCahJodn&Rm0wpjPpegPkBWRMarVS-*)CH>7UNui$dFy28C($pA zq1u&ygMz^OQXY~)rvM3gxYKXcI{#@nok_A*r;T<|P%hX0ib#a3WS_R`1-v`lWVU@d z|7TjKmvDwg)SJ$C1I}6BrCtblw#C=>edoJYfN=^pjv9nAIIlU;eIsGhbsTt*VgA`ru*cWnWJ66VkZ}yG@t6kTC0c zy;$De@%0rA;u4hn0W%y7I)E}*ZaUR?Pmt4KuXXC)gno$b#Ks*w3z!9oUIUfZx>r?~ ztEyS@Pucyfi>IO-4@VoE5_P6v3AZ_9`b2Z8lLJotU%-#G1NX}^tc88&2J@f&op3E` z8?B7B4NlhWxCf7hu|M!cZxYD}xNTziq<=`~17*tlP?;zxf#0P}Q*$`40{8{=h# zv=^{`W|k6%AyJG1h#H|IgniltB`Es{JuW$VTS5J$|`~DPodKYU78Xt&DM4v5o8^%ay(L+6|f3*s;*h3j?ui z#ySlugdiR!uBNDdG5)$_HFKHy3FYGi&IsYiI zfD{p~HPegLbeVf*m#{bTZ_rU26i!vky&hFRbr2g_C8)@Rm)crh9TdyH$Y;E00L@<{ z->0p(zH7d`R@Hj&|3#nhmx9G&mAq)jdT)j#eTNsKwIGR@wgv zTXK1c4v5UGgeEJS@32Co`z-(wB%TZ+L)XUjLIdL8*KlSG<;;J&?{&?GV{Gh>aFPVO zfxI%vSWPZVGv*d1Xls+o*eeZTG{7 z>uX5rrdPp0gY<_*?Ji|2Z!@->y~FcE5>B1x7iS?Z?P0_%spLVU*B&q0Q>SSPi&vtP zr#4`t(xbdBWnkSCuW@KOT_Iaar@jM{N(`nCh z=0UnqwR5AlmPZnJKN~n^s0v}1XiQLSa3H8ky<;>vbBK2k!?ZP_8)yYLUaRoFha33n zM05G5Z#S2)Ja(6}8bt4?0z!bd#-lS{OxVP@dzZqYN{>|J!1>7y2vhdaeDC7DM=Fe! z-0j45pShMqdUXq5EX#+bV*)67$C0QAg38zS7Mq=i-;suT5ShQ)w^>dlU1%b(z@0v~UHtC+9;2N!qQt znA*ctT$(StDQ+sGeAO=})-vWUNj9D95vj}_;PeljCJB>^)35M$w^RJ|Znf1Em!qU;`VgrW}^S4G=u%J6`76|fRy(1*d!CaJs^4|uewAoo;OM}<4{6dFwG3A_BB?SU) z@Eb+j&(}?|LgiC<7xq|D?FXljwh-LCZ_)&4G7GSD}zdSE43x7 z_d;}+g4=rRKK&Z0msM8sHrd!fvf?a;0rcgF9(?d|F+uU4c9O~)O5?pMlmp?cje3C{-} zli<$$eFcTqMqx9fg$5>Uv_IJd$5>`Q#b;!!m0?75F+?&Lz*2hkO%r#rox_Ux$sxYd zPAGZ7QZ14$V$^9Z9Ml{LwpWaf;+-QGe74zM(_01yA}M>GG6owK=s_WHpJ323`Wd*# zg1U;y6n`bBSClg{eRk=YpT+Uz*VE<0+a$x}Yyy6~u00xBOT4y>&uF%{)i$7S=V=el zq%us1Tx~I_o38k*z4DpTs0gmrqqYS-lR%_SVdMKbB``ZcBXw4m@lNPHi)UJ|FLcM+ zUb;ZcA`~c=@{y9+b6fxv7dG)o`G_sFl#Pk1>JYhp`dfnk%L)IMfdub7$Yi$iQjvs6 zJn1is?|7&{WMl%54#0O~#S*NX8$k6xgisxXZ$BZub`Jbn|MptJSxRlGHg{C4#>#*TZEu&5*54gyg9&0Qt^JxW|k!lWaOFojR#8MX1(FI3cm zDuqRy0=@%76fG5%a$}nBBI}PS+{%r2zATBE4ZsbjH`U1U1im*C4`VY?vmPz<=c8=%5yd6g Y*f9EGOxo +
LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Commands +
+ + + + + + +
+ +

angle_style dipole command +

+

Syntax: +

+
angle_style dipole 
+
+

Examples: +

+
angle_style dipole
+angle_coeff 6 2.1 180.0 
+
+

Description: +

+

The dipole angle style is used to control the orientation of a dipolar +atom within a molecule (Orsi). Specifically, the dipole angle +style restrains the orientation of a point dipole mu_j (embedded in atom +'j') with respect to a reference (bond) vector r_ij = r_i - r_j, where 'i' +is another atom of the same molecule (typically, 'i' and 'j' are also +covalently bonded). +

+

It is convenient to define an angle gamma between the 'free' vector mu_j +and the reference (bond) vector r_ij: +

+
+
+

The dipole angle style uses the potential: +

+
+
+

where K is a rigidity constant and gamma0 is an equilibrium (reference) +angle. +

+

The torque on the dipole can be obtained by differentiating the +potential using the 'chain rule' as in appendix C.3 of +(Allen): +

+
+
+

Example: if gamma0 is set to 0 degrees, the torque generated by +the potential will tend to align the dipole along the reference +direction defined by the (bond) vector r_ij (in other words, mu_j is +restrained to point towards atom 'i'). +

+

Note that the angle dipole potential does not give rise to any force, +because it does not depend on the distance between i and j (it only +depends on the angle between mu_j and r_ij). +

+

The following coefficients must be defined for each angle type via the +angle_coeff command as in the example above, or in +the data file or restart files read by the read_data +or read_restart commands: +

+
  • K (energy) +
  • gamma0 (degrees) +
+

Restrictions: +

+

This angle style can only be used if LAMMPS was built with the +"user-misc" package. See the Making LAMMPS +section for more info on packages. +

+

IMPORTANT: In the "Angles" section of the data file, the atom id 'j' +corresponding to the dipole to restrain must come before the atom id +of the reference atom 'i'. A third atom id 'k' must also be provided, +although 'k' is just a 'dummy' atom which can be any atom; it may be +useful to choose a convention (e.g., 'k'='i') and adhere to it. For +example, if id=1 for the dipolar atom to restrain, and id=2 for the +reference atom, the corresponding line in the "Angles" section of the +data file would read: X X 1 2 2 +

+

The "newton" command for intramolecular interactions must be "on" +(which is the default). +

+

This angle style should not be used with SHAKE. +

+

Related commands: +

+

angle_coeff, angle_hybrid +

+

Default: none +

+
+ + + +

(Orsi) Orsi & Essex, The ELBA force field for coarse-grain modeling of +lipid membranes, PloS ONE 6(12): e28637, 2011. +

+ + +

(Allen) Allen & Tildesley, Computer Simulation of Liquids, +Clarendon Press, Oxford, 1987. +

+ diff --git a/doc/angle_dipole.txt b/doc/angle_dipole.txt new file mode 100644 index 0000000000..7ead336f8a --- /dev/null +++ b/doc/angle_dipole.txt @@ -0,0 +1,98 @@ +"LAMMPS WWW Site"_lws - "LAMMPS Documentation"_ld - "LAMMPS Commands"_lc :c + +:link(lws,http://lammps.sandia.gov) +:link(ld,Manual.html) +:link(lc,Section_commands.html#comm) + +:line + +angle_style dipole command :h3 + +[Syntax:] + +angle_style dipole :pre + +[Examples:] + +angle_style dipole +angle_coeff 6 2.1 180.0 :pre + +[Description:] + +The {dipole} angle style is used to control the orientation of a dipolar +atom within a molecule "(Orsi)"_#Orsi. Specifically, the {dipole} angle +style restrains the orientation of a point dipole mu_j (embedded in atom +'j') with respect to a reference (bond) vector r_ij = r_i - r_j, where 'i' +is another atom of the same molecule (typically, 'i' and 'j' are also +covalently bonded). + +It is convenient to define an angle gamma between the 'free' vector mu_j +and the reference (bond) vector r_ij: + +:c,image(Eqs/angle_dipole_gamma.jpg) + +The {dipole} angle style uses the potential: + +:c,image(Eqs/angle_dipole_potential.jpg) + +where K is a rigidity constant and gamma0 is an equilibrium (reference) +angle. + +The torque on the dipole can be obtained by differentiating the +potential using the 'chain rule' as in appendix C.3 of +"(Allen)"_#Allen: + +:c,image(Eqs/angle_dipole_torque.jpg) + +Example: if gamma0 is set to 0 degrees, the torque generated by +the potential will tend to align the dipole along the reference +direction defined by the (bond) vector r_ij (in other words, mu_j is +restrained to point towards atom 'i'). + +Note that the angle dipole potential does not give rise to any force, +because it does not depend on the distance between i and j (it only +depends on the angle between mu_j and r_ij). + +The following coefficients must be defined for each angle type via the +"angle_coeff"_angle_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 commands: + +K (energy) +gamma0 (degrees) :ul + +[Restrictions:] + +This angle style can only be used if LAMMPS was built with the +"user-misc" package. See the "Making LAMMPS"_Section_start.html#2_3 +section for more info on packages. + +IMPORTANT: In the "Angles" section of the data file, the atom id 'j' +corresponding to the dipole to restrain must come before the atom id +of the reference atom 'i'. A third atom id 'k' must also be provided, +although 'k' is just a 'dummy' atom which can be any atom; it may be +useful to choose a convention (e.g., 'k'='i') and adhere to it. For +example, if id=1 for the dipolar atom to restrain, and id=2 for the +reference atom, the corresponding line in the "Angles" section of the +data file would read: X X 1 2 2 + +The "newton" command for intramolecular interactions must be "on" +(which is the default). + +This angle style should not be used with SHAKE. + +[Related commands:] + +"angle_coeff"_angle_coeff.html, "angle_hybrid"_angle_hybrid.html + +[Default:] none + +:line + +:link(Orsi) +[(Orsi)] Orsi & Essex, The ELBA force field for coarse-grain modeling of +lipid membranes, PloS ONE 6(12): e28637, 2011. + +:link(Allen) +[(Allen)] Allen & Tildesley, Computer Simulation of Liquids, +Clarendon Press, Oxford, 1987.