From 59b59573a736057543bf9ea1ffa2c2c5ff2ded57 Mon Sep 17 00:00:00 2001
From: PabloPiaggi <ppiaggi@gmail.com>
Date: Wed, 2 May 2018 16:10:31 +0200
Subject: [PATCH] Added documentation for compute_pair_entropy_atom

---
 doc/src/Eqs/pair_entropy.jpg          | Bin 0 -> 8870 bytes
 doc/src/Eqs/pair_entropy.tex          |  10 +++
 doc/src/Eqs/pair_entropy2.jpg         | Bin 0 -> 8474 bytes
 doc/src/Eqs/pair_entropy2.tex         |  10 +++
 doc/src/Eqs/pair_entropy3.jpg         | Bin 0 -> 5759 bytes
 doc/src/Eqs/pair_entropy3.tex         |  10 +++
 doc/src/compute_pair_entropy_atom.txt | 121 ++++++++++++++++++++++++++
 src/compute_pair_entropy_atom.cpp     |  11 ++-
 8 files changed, 160 insertions(+), 2 deletions(-)
 create mode 100644 doc/src/Eqs/pair_entropy.jpg
 create mode 100644 doc/src/Eqs/pair_entropy.tex
 create mode 100644 doc/src/Eqs/pair_entropy2.jpg
 create mode 100644 doc/src/Eqs/pair_entropy2.tex
 create mode 100644 doc/src/Eqs/pair_entropy3.jpg
 create mode 100644 doc/src/Eqs/pair_entropy3.tex
 create mode 100644 doc/src/compute_pair_entropy_atom.txt

diff --git a/doc/src/Eqs/pair_entropy.jpg b/doc/src/Eqs/pair_entropy.jpg
new file mode 100644
index 0000000000000000000000000000000000000000..44161e6aa0e74a05544f8f29aab3ca942ef305a5
GIT binary patch
literal 8870
zcmb7o1yq&Y^7nIs?v(Cs5DwiO8tI0Ez@bCBLlluxN*Y8F=@t&9pwa>&Ab9|3MY<b4
zeBb~5-n;Joul22O*4ndwduE>5d(X4ij=7n-SpbMNRW(!r5C{O9fz1HmW(k)`Ls{8I
z&p=mILmU1_0pJn?w*<cB9-e+a25L%7X66=5e`Nn6TYFzGMLj+3f3^R+J{|{dJpsTN
z-+$}=pAuZ8gRlLq)#Ysov-k4x0{~FaE#?aJ^ZJ8d-C}Y(x4&@1AMA5Gz*~9WAMEru
zp8d<_Z@l}Ljgf)!?W~{xfQ|3;Z`khN@c;CCYYwm|`uMpxIs`EpK5+C4wDoad^8FKu
z+gLH_`8;rRad-G*^*{9gW&E#?{-gGD2)dmc04P833ifeva`t2TuMvm{2@3=N{EGwq
zzhVrc9Q;T2kJEo-4+;T5;s^jB$^XdgvH_s!2>?*d{3By~1pq|R0Du~{_xB0-+tY2h
z000BP1&9E0fEr)`SOG489}oe=0ck)1fCHL<E?@+h0|>w#a0Wa8A0QA21EPQ!AOT1Q
z(t#YH04M^=fNG!~Xae2?pMV}<0Qd%s1886YSOqqKUEmnF0D(YQAOa8>h#JHQVh8bp
zL_iWCIS?GA4Ke_kgYJWzK%O9fP#EYjC;^lT$^pFwm4WI&Euc@Je$WVL8ng&n2knE-
z!4NP$m>f(G<^T(VCBTYcEwB;T8ter228V#7!Oy{2;6iXE7zO?W9srMn7r-0fBM1P&
zgHS@4A^Z>th%!VMVgYf4_&_2c@sJEiA*3483h9N6K^7pJkTWPIlnlxQ<%dc`)uF~v
zd#E=w0-6ZTftEs1&@Sj0bP>7>y~ZHGpu^zBkiyWwFvD=d2*ilN$iOJZXvFBo7{^$}
zIL5@pq{QUJ6vtG@G{<zs48=^u%*U+3{D?V%xr}*)g@r|p#fv42rHf^Y<%bo6m5o)2
z)qypNwTg9$jfc&MErP9zZGr8H{Rle)yBxb6dldU8_5}_x4m*w%jxG)oCm1IQrwFGN
zXBcM{=K_}mmlIbO*9g}YHwrfkw;Hz_cNX^m4-bz8PZCcb&js%xUKU<0ULW2f-WfhA
zJ`cV!z7@Vd{&V~i{Ezt4_y+_81RMkk1QrB-1kVY|2s#Po2~G*g2n7f=2$6)5gt>$$
z!coF)A{-)iB1IwuQ3z26(L17HqD^8PVh&<uVq4+};#}eu;z{Bo5;78D5?vAxl0=en
zl75n3q*$aJq;OIsX*B5@(oWK4G6)$fnKGF@*<-RIvTm|fa!hhgat(4<@&xiK@~`B(
z6r>cQ6ebiw6uA^1C>AIolpK^Al<t&Cl=YMol;?No?<n4JxD$7$`p)Q`6DnFN1u6%s
zc&ZwzajJ7_CTcjf8+8hG6ZISol!k{!pC*vzB~3TYCM`LwG_5^tJZ(Mg3>}z`ht7a5
zl<p1PAl(r?1HBr(7kxH;C;cV^C4&NkE5i$hc7}CEGDbN@7seNi9gM%3D3}zP+?leN
zx|nvE>6q1-{h41ge`UU4;bbvpdBjq~GRKO`D#7Z=`hxW{>n<B3n+{tTTLs$;I}W=9
zyEA(xdoTM52PcOq$5W0bj$fQqoEn@VoE4n2Tm)S5T;5!TT;I7dxh1&WxbwNc@ql^6
zc$|51d4_qxyf9uD-hAE>J`6rdK2N?vzDa(3entKO{tEs@0ZIXFfky&O0^5RYf);`)
zg1v&*LSjPhLPbJoVNzjD;YY%)!uujTBK9J=BH!;4+*Q3BdAIrQo+z)VgXl}qDKRoJ
z9kHikpT(|V;xHdrHEdm+P5i!ip7^8$xrClXqC}r0hNO~Yq~v?aGbxyqpH#injx@it
zn{>JKFBwi52bm(7Wmy(kTiMsL3vx_y_vK#6Eyy#=+sYTpFDbApAQg%gekyV+x++#G
zZYv2X`6x9iohVBwhbw<nfvBjd#H)OTlfsSRIq-Q^7F8$J3e{aTF||;&kLnoeTIwn4
zlN$6I_8MgxyPD#f;hJ4q_*#ZqIa*8FJla0mtvX;GO`SBI*?a8wJnuE>0=gQyX}WWI
zoO<4RZTgt{diuHgYX-swVFta1WQI0|<%Y*bN=8XWXk#v8f8)<4L?%`yWhN)4DyFHX
z3uc05;bsHoH0Cbm%@$Y|CKkmON0xBQbjwvMF{>D>32RR4AnSev4Z<Dq!G_4j*5=)P
zjQghd%kN*>>e?3C9@%Nw<=O4pE8Az<Zy*(rFOcgFvJR;ZKOJQqQyqUf$vUMutvky*
zXE<-VD7)mi?76DD7Py|c>AID;-ME{&*LYxi-1lhlB=vOj?0UfPAn3t}7q3^0*Mhf{
zce?kkkCsodFUZ&0x7m;U_D~%3=kkBzzZf7FkQ;CsXcG7?h$P4}XfT*3I4<~Sh-%2&
zP-v)q=;tuzu&A)baK-S~5ugaWh|iI%k<pQ>57i!)M&U-eMGZX?dX)O;@UiLR)@b_Z
zhtbPV)Sgtt5XSh#Og@!<S`Z73b&35NcQ-CO?m8YB|0O{vAtT{3(LV9ZGvQ}h&u*SO
zKL45oOM01%ne3T7ks_Z`no5!ylDeFxlZJZ1^dkPnVY*Fve}-sAK_+gdf967#c2;vX
zTXu5xWsY;sc&<`zO&(odT;54OGJo`?!prIc`htXl^H(mfre3SPMip`tX1u|C6YyrO
z$fT(Et;E~1Vw&QF;;RzRl7&)((ylV`va)iz@}vrIg@46*rB&s(DwV3{YX0iiHIy|8
zwLq<Z?eDt#brbbk^`GBKysK_tZ^&yTYm7&MP{F9(Cg-MwX7lEe7R{E=t<tTHZ31m2
z@0s7{exUe})Q;O8-G0*%(sA_B>*MApr%%hDZ9bzr%{sq#>34nY*6Qx-QSIsMRqXBP
zlkI!oFV)}jMdC}-0Bis?C^pzQBsSFe74{W1EI!=)P4Zjoh|Ea)sQl=sF_p2N?;76+
z#_x@fOc+m0O<GMZOd+R!O}kI;%>>L`pdZd+%*M|V&t=Th%@_XQ{!z04Tllc3vN*V8
zxP)G|U*1^pUAbJ1UL#n`_{sRQ?3c)|_v@<b-+o*EUfuB8xY&%@BH7B@=Gt!9QP>&W
zHQQa?^Vz%EPduPLC^-~8{CuQ$G=J=Id~p(gN_|>-CUMq#ZhHRnBJdLHGV6-_s`Xmy
zdiKWS=H><sC<EZz^BDyC!@++D3i-n^pil@D6XW(U`%|!Rux|?%Cg$xCg@gNtZV^5n
zF8-g!%|`$a8;}DfgFqy}O&36Sdvk&O<@4`r#x2Lh!iHi%aKN`(1!4dMfr3F`3=k9x
zig&97-?jmaTN4seGG-id7GXtw3Q7^%J5;PnhTc)QY?6VZ2Btni8SffVqhmKS0QR4I
ztUPda`IcML;Kd)7o^a^YIcc*ZQ#OnK;pvinkIy5)DZSdJRRe}d-f9=pmd?sW${`xo
zb-^htuc0O%Q&dQJH{!zJO%s|hYD51R>IY9kn;pS{kPjN(%>nPjKYQv-rO@ia;=ZfG
zi`F%?b%#id4As0t{c2bwuHrlctN1P7B4J!(oM`Hl-$mS*Q)!98hxw(Xr8TP*6nEw*
zOB(*~*!~*ef3Uh?%793gGsC`t`PuOUuiS{lu7k}4Y-wl(_T|=R84gRA={iPWOIE8;
zjsvPb(?}PQDrTS)?$Ti4PC4N|!D|(yw+w~`rBG3#oeeJvSz;d*^z#<}M<;)p#-K{^
z3bbW@a70~ZG}5F=h)%W;hc~%ppVT5-Y~o{)qM`y0);`AuN4dlP&mSYd51X!T!&z!Q
zum_n8vFlECr#{b>X}6U?mQW46OY|8W%%beR-R%k<wTrMp{gu?7#pEgIvz_R%X}!T8
z^!A90wK&*hxL{a!{D7X$cWF`f9V*k$hr;GCF5_@m3a=9o={)IcPX$-@z;lnRdiC!O
zKqRM3-os<ICQVpdXpBx$imi47=SFlZ`?2OQ>MdH!G;7YR+xYW5mKi>+U@VScpE^?e
zyid(JjIUdk_?2{?To_qp7>Ob4;WvH*;PLaj^;sIVA-s<FX-6{mBt_O1XBk9x_Dk4Z
z2Uh2&W^$USyI#B0*NTM=_@X0Z2H$7kjoNkNfY)iu#{@+kOP~40#uS7!R6X{=bN6~L
zGmtv+P~)Dehhn{N(Y_w3u|5Uj9Byo+`{in*Z-$naU@T<m0tI7V;6LW``F_>2B{bIH
z5)-98PPSb!8p;w{ei*mi_|v)aot;zbm$|&nO^v*W<sgi#p=bgfw4bS5F?Cpa2sidj
zQrYsi*X8j6k5}5)kS!s0nc$BG^G$WJ4sNvNj03d4@?%>TO4WJVzm#p~pNt-j_5^Kw
z@V-|0J=I3p{~+WClS6-#jeGqwiVK~uyj@+*o=2NS=#Ta@p`<~np`=DiV*@|geR~!^
z<xUhgE2MsHx!P-3b?~*Fnxu*oV$OB`b^|z%BoudaZ<flI%vZ-7R1ly))KCTVbT+8I
zYv^imG*27W@M!6>5p-*Sn@l)-I>b<7W;8C&g&G<9@80{S*Uxued9r=n;(dJs#3QrC
z%7qsubH_}yuBORTHbndHLG+M}9dWASaa3Z*dYm}|-}PGG-LA0;3Co$)EFZ^!m`PM=
zy{#Q3%F*eh462#J_yVf8wHaF;a0U+9=DC1P7U%Q4=Yctxq%~Duil=@PRsGxqPa+aR
z$qEAR)eKZJF7-&W9|YqVlk8y}p^c7VDP8Ln=?3*1&o`3IDs`2TLf&`BuJ+ieG(~52
zQEL5$I*&P59OcAP2bFXvRGeCDmLH^E7L4cVOuR!?qqUZ178KL?4{Miww(q^spXSVG
zry9zEae*sVmSQg=Db!8`_J_A+L_^bNJ+mDWe=oYEe{vIQ_`uxs>M~ACRRhK795u=a
z-9KQX*qJIHm<*=N>$RDxU$2DdTD$K#%8Xf`MoE3kDGb~h<u<J|tFfC27LI-}XB%&C
z8L=2SEpFY$v6J2C)Y4BYDjRXHPxqdmS{{cvv)gJGGM8AKEk2p2q`a}=B`2s5mV4WW
zKH>Sk%~(gbp%u^S;sMc}ArbhvVV%wFT!LR0!qa*r+a9I@#+(l?p*}b6(2JGdff@Tw
z{ygDsK{xxK@c7K(WFV;{P}#okr#*6L`LpA5r{HtbAjaY2C!ex*xY~U*POqaX>|zND
zc}U8NlL{(NCwm1H5y7-C9LJWB(|d6(u09%lI7IJAVtWpQsz5RbKA~?ldk|;U(udWl
z>2E4g?Woz)R>#&7v8ENuC#jp>%`+>kl@(Q)2@j)Reo;~{)Hr1`vhg!2PfD?RGTl09
z|MG794grR3Y;ul2|HU1S5cYoTw$1cfZ?Z8&d-}$Y`B!fSWfzL#uAEO(wmTBnbZ>wv
zH(Ix^8ck}+2ZPO?B2Z6N1e)ClZj5mme=VPG)DkGbu@tOwiIjNTEKmm4ug&CD<r=6%
zp4K4VsMf>P?`+4H=3l3e)kx~zz3(VgcPFU3ag$Xjj&dt`8NZHDUg<6RIsPMk_(svW
z>6}3sgT`+<>E}`qTkkv3Q_AlG<e3D-TL!1s<`~i$pb3Q`hC2r7qoN5IFT~_YRFr+Y
zYbjVu=O;(=(w7Pdk9y*lbp_;oHmOv_l5D3MzfRtN?AyGWSu*E}=lpXBwXCg1IMPko
zt=%SOK>naG&+7qVzx*Nx0W;}8jafBclT2DOt-uTVz6t9liTlB^sKQAfz^3P?ZiA&k
ztgau1)v(ciSy_`N%^0Y*0vqXpZEg}-|9aIj?UBi=9n^I+e+DKR9L_jvN5|H^X}<qa
zgiFp5{l&Fxi90)RNBS;(Q|fMs-1Pam`B1|K=ZQd<9ecTQGM^9O46*Ik9Q?S}LDhm>
z?o}gFBh?hyjyUo8Cv~3g0iKp9kqB+oX^8;1GitWhGBluhPfk6;?AeKdmfO;@0rtaO
zxQgmW!`XE8jE;Ag`re#D%h5iQyd&m11>mncs11ku6`w_U_S6C~=B}V+*nYO?a?#B3
zk8CE1r^};)78$iR7wTt2OGCR&ePV=y%t{BS`VB$?H;_PFXMt>up5f6}#+>fH1v=-1
zsA*#mv-hEeB}tt9yPl?Q;(IBy4zJzy$qs)Rlt;MV0HcRLFL&!B*R@JxdVip$eZ!P!
z7c$-qu-Q{;OGE4hr>Di*rJclfa+BRHx4h%N2{{dqg#WbOGuU-k2s3(-c6Yw1XKy94
z^NSeTdu=h<{O)svg*RP@51Kl0Wzay*lz}XwTr^f7_Xo#1^<k{~@59)978l)wnDu#?
zpe2}Ei(k7VO43AHYXH=@)vyvb*)(?>VH!wPRvIUVT3T#ubm>$)4>P}R1FU$}leK<W
zV)in_jiYW>sIPLctC!wYqN^jM+=pE)=uX|U^>R^sQqn6q6(r&ayhyffc2+Pz?+f{b
zf=^BG4hE2I&Ssw;vZj_6B(XXWL>jQz45{rT($zd)p&HWke&icR2DNeO9cr{nDcWab
zGuo0{pW%sQFTV@S%}xqF?ch4Y3tJt(cDMmFoq`8tz?K=gK76sqrr0yI5AkiWF0aH~
zRESAH#)iKUW08x}5#-krGvaHmQMh-fw8I{4t+SNY;^8CPCRg<g2Ws{!SyO!>Il;r2
z{x^Vl?N;r!SjG;2_WedTI&Ak5q-7P#GXE0|lJO3$DdUS;{eH)j9HtzmSSH-p&N{*Z
zs+N0>^Wu1IFX6GSgRmhdkKyLK{2N96nDTwIlqshBl_l%?_-`au*R1yh61pz~Lt7l`
z)F+z~OZej+L{T526s`)pq<QrdXZBr{=U+52h#2aU3X{W$g<tqI8E(;<e!k=AhEVZ1
zJ$q^G@jbJtx+OOtr75~BQ7uH7D1k}y884E+uk^$5%$+&cv*l}R-zgeeoMAx)dy$h^
zw9~*HvO+aPJUTkz@`Ak-p+C^y{8dAqo{-a%H&roX8ON(}NxO|kH;yOd8BIT&3DE{O
zx;GLrj;@x33)(pf3u6ThHG1!*|C~d!^;jek6c8I+k9<mYl|32_=bb5M_X<;xbi;dm
zsXlPVvJGo0n=l*ft(_YAK{R73soO+J=P$S&*M8uv`SO>XX1Ye-uw-bZ<m+BfDQijS
z(e-M++0HcO@)M4JvQ?@CFt0i@2vjFg#W(#-^rZ??VI}~9&@@PjiN0vO0r2P1+|5H8
z9f%f!Y#iS1B>_HaZqL}!%sHreY~Kwa7*c+i-n)eLMdn9x$T0qz)N<8E#Z~pEXGr!$
zcm^cpdMJToGg=Px`?EvucN(V#HWr(^%#0X~?VJKS-Fy7h1x-v=@0Ty*J9CqSEO+-c
zQZwe*rF0j`>@54rE=PDZS6e`KL0y+G+{I}NKgB;_ZylkhPv`8GKAgcYQWCph;VJ1@
ztJ~p|jDl^YuRq0$v&<0e-=<0>c39#AcZ^rB$>oM?*s$6NoF5;PBCaab*SD@XJ?9;!
z=SOJq$kgjnRx-1uEMo_L7!W;15;3M+oJ?wL7KFrKYDO;3edt)Z*0=#s>ZLDbZvZSV
z`BNvi8z70@LSoE(^EtP<-S2f0^FG99rY7eU>A7MCuA$ggjnOlvz}4lmy^RC&2T#O3
zz4lM6hG(01gcwow6VLX1AIF7V3sY(!_IUkA`gRgiNgfQH;SN-xc25&L3I*eG&e+>r
zP90s{5H~<>U!Zu*T(;b5hFFS~_gNdx!8^^25!_k<Wkm&!oLwTF#&fkXhW2p_Yssgz
zCnh@m{3Z4?cZI4fdztbb;zY(NB|n6Be2zNsU6FebRoL^q2o3ZFoyR=qCM}75%j%gw
zz(8`zD0;MG>R35_P?_gJs*cO+ng8={e~ov62@<tHTIq$8hVzMn4Dqy3mvAH4c}}PA
zLu7ru+wsIoieu2+DB6>jJ{2C7HEhWVmU@Rc%Gu~IQAaDd^<D_OkoCo>XW&%H(M_Q~
zd_PU1otwTZKr@EhF6`fEPl$fJ#?fbIr}O-sBXxG4LTc>(X@cuwd@xUm-R+9wE%XR0
z33QODUlixvoeks~J6%aMl<!9IM>$Tv6<g^3ZaD68JRkitVd%9KgSy5jS50yHs~`tl
z*0B4XKA4RK+kTa)-y|uPX-u&vFnU(Bq%2+GBRrQ~Pa}RMoUTj#E;hZgsS;_6pWDw5
zv}@7QQ!qEhQx=pcEI{rH8?Vc2E?qx5OuBk**Krmie!I)P5FK+z@zqZ~1YSMF^?MX|
ztfJc;X;c}SeJOMX3ZD(mSnibLv3yLPl5G_pOw3QI1$p-QzCJuQF?Uzzlf!FcMMl9E
z62b4t3Rgkwy{C4+He|M*g{HfYibMp^VO@!79qu;e!5jvSR5sB^h+txU5-Z=1?YH5L
z0_ok&+~u(6TJB$}Cv$Z=`TBf+GHbi5S8edN_^=2syj1;BDv)(D9q_wk-MnP+<fBZR
zTbRV2%-xh+*ddIdbb0z~|EG3E6rzA*gR5H5axQ)I7o6P>V<LsPpSbh80z1vOOq-|2
z^&OrBZ&B6L`~w4yznd+t`0Er=mbKGUo_fhD`uDoE-;&fs^i@687&U_U{n|?DDtKBP
zTq9~I$q@O{T8LvWQQyNvh0KAviJ+a7Rfp#>e?9$sm5L*3iLy&2Zb)S8sC<!t=zS}<
zu@-G-srVg?%{LswVAIXmW9MB$>6GQpmkvM0CQ079^w3b)7Ygd>hiQdOKu)X8W`8ge
zUz6_HzCJ3%9BOr^-HhBW3Z)~5#ikva7LeQkVpV5G^3nns16;0nRclLK>Xq9U%(5JY
zpV1yAZ@Hye*j03oDwn^=K6yCr^WmUHqA_slSf}gip+CBLsD_@=MmU%LO$eF;L`4nz
z_~HN~PmceME3!O(D#vVc8{%M=`B|BGNua|e^PRPKqM1Md24_UJB*EL$6X@pTrHghJ
zcg^XS3KP|?35$*;c#+t=YB0Qq;L%*3dHorzc#7OdCD-luz9oB9iASj>NKLtj87AYp
z-EPA-HI#u+OR=VU42Jrqz>p<+w6?D(SVY*5NmSs*##cN~7Z?b_11Nh|wyhn~9yEJu
zAE7U)?&j!4h-GIDDxe^b+85xITyn44m+Fm7-@>a0YVdXpTl(@A<2nRJs`(eH?X>MR
zkTG89CfarGh(;n~tFtJwgZ&H6Bx_@Thf|Z;%}XTuz2!Bkg?GYZApHHDh~(-sbgH(g
zz!x_!>?`plmY)S<E74U)o9#yqo8n&Abh}F7HNFx1n{<JsUEr|*&;F!@N!}M2B(V>L
z-<%ef!<?hvEXA^@e(-C|XpE#`9BfJJwKgq72^eO%%u6mGY(J4&CMfnuUb`&lGhcX9
zh2oa8*OA#vDa<S)z;fnTDs#448oIK%W}cW=m{l?F4rsq~hp>HCBgB(n8^2`55@-AZ
zyU)`KMVuaTq?mlsEOvk6tKf1L5m$g9PXc?c_lHed>zpZlBno@HnxdsC1x0OiVkC}~
zSr5Egx*~lYNvhK7E+?zK<VcEK(`}*f4$^_;S6YT%5q$N&0mfZoUH7;J>$61}eaS7C
zKG~y|VO(Z&IZ}5Ce~>hf(CF5RPQ8k~0SZ~z{YT#oclnq^_ctQsUj^eqI6Wt8_ouu(
zS}3JG)gX+`aKRaplK!9(%S5ziH(^|t5sN+%#lA&lvux!t;mB#!mZ!&>B>eUjoY>P@
z6(rZ({-f#*PRrKguUajAPpAP|iK6pfcWWaXvuqJ1<~@tr(?)waiT$S@-pA+if}X=7
zO~;fFrnKVKXjC>{k)BOc`Cw1~Dp5y?5Dr<9K*{C3r3taFONJXD&dkSUvP@DgFyzZl
zpj$HVfyJZ8yie4K=m9YW=F@|vZzp>949w<M)HoUi9@dN0^G~UoeEm4g_K2m4n7QE>
znkDD<p?=s4TD<j*-yi-ybR^w+r%EkyCUEP4U)@mrKu|=d-Iccnn|gJxC;rrLg2ROA
zGYyROkbB7=uKQ5zCNj;-^NL?|IKxK*u0DBnZe_zW9ExuMo=KNy>?e!;$ZRhV5(f{y
z1Rr7Hxww_r+^w?VOt2x?95ZmdXx`^Cp*QtvFxzfw_P*bmGTo5NSC5FtmT-+aStrlj
zGubXj^}Y8-(>B1HTL`AjtYDTIUs!UX`hZ^dooOTE(0jKV_Ce6Io<Xgrv#y;%!;k0K
zby15GCJT+;c`Qjo-V{31giIH~&5P8Th`F5cXXf?}m;7}Ek1J7YF{tgq1jbL|42p7B
zsb-sdtO;AQju02lu~enE``Z%VTD41$HC5#zpsx+wegeYoZRcWvB2y=l)oxcg+21qJ
z@n~&?fKF!!-+kEC6{VRPt;t7YV>v4l<M8otn^;G(Ja-1Jm7ukHv?Pi;l`O@o*1GZF
zp*W{w!a}7gTu4|o-}h4!g<+BVGh=-c*Ui!=!%62nN6FX6^P&@iJiqGY?pKTO2@fp=
z_7e^1C^rhSDQ%gyk5Iln{F6(Jk1ZR0OCoiAo9z7Q&=SM>a5?pWZiBOF+Jl((h|^l&
z9atu4afUOcdFj=ofYSW(lKFbM`rginFjwEJm40XiNBVczEXOpPX@A3eYa7eF?-&={
zP`yX&ljsSRsx*Y0cTFvQt}UM>a*;ADX38_+sYkmdO~xmx8aqvNp;CFmF``ZLa8c>0
zN6dKZ97*Zkasz*DC;uzReeQ3u@u4h#D3sIqs5&&fc--aa25>;z_nc)zFAFk}%Hw_=
zai7=a3K41bHvYdeEp@8T9t4bqqI$C__Vz+Vi{U-9_bz@#r|5@|9axjU<+eJpgqL3c
zTY+uJ;1uhP@6WhXyx*<s8`iNIA<{_<&k9Z9zUjpiri$$pD5)V#oQGMJN#*?9%v{Vz
zGq~j}<$<H+EZ2YXp}a4xizmI(XxuN{1LvcEw-HQclX4Ki71b#W>x@a6w`9-Rb?YbL
zzMTy-Rr|=HD39lzW7cYpQUfe(XRlr{(-a11K4Our=}?M0*8KIB8v*mL_CNkscfIL8
zNIH~a8EzYT9DY3U9VM%{({&}T^0Z$Xky$;~SYsg*w7c9y$_k9>fFU|CE*5ay=xuyh
zSR!~Y{yj0QM`&!(V7>fox{WTkwzFkpO}NZw|I`p|=+egrt7$HeEiEs7cE+ihFja!M
zGFI&TFkIyomUicxm&E>Yt877n_Sd~6i{rHUwwgTu?1O63;KW57VeRC&@A?8yx)0OI
zbpK9||1XK*zZBJH!6ileWbwJ2Dxce+SO;{J8;*9qh2{8kS<gi|)MLvF<DDxeJ2cF8
zqG@0L7`irDpPlIRjQEzLtTCLBQYWN#b-wB&k$j$52uFi1V2Z(vQ8lGM#xAWRqW`y!
IO>d_E4~Hw}5C8xG

literal 0
HcmV?d00001

diff --git a/doc/src/Eqs/pair_entropy.tex b/doc/src/Eqs/pair_entropy.tex
new file mode 100644
index 0000000000..304c9d6138
--- /dev/null
+++ b/doc/src/Eqs/pair_entropy.tex
@@ -0,0 +1,10 @@
+\documentclass[12pt]{article}
+
+\begin{document}
+\thispagestyle{empty}
+
+$$
+   s_S^i=-2\pi\rho k_B \int\limits_0^{r_m} \left [ g(r) \ln g(r) - g(r) + 1 \right ] r^2 dr ,
+$$
+
+\end{document}
diff --git a/doc/src/Eqs/pair_entropy2.jpg b/doc/src/Eqs/pair_entropy2.jpg
new file mode 100644
index 0000000000000000000000000000000000000000..db3dcb0ca75f444a6640e704c61326d7c8c598a9
GIT binary patch
literal 8474
zcmbt&1ymf(^6%^}ySPh&!{Q#CKyVKp+#y(y#R3F}1cEG12m~htC&Ar=JHg!{1h)`8
z{K@y-d%k<#J^%C0JFjNC=GWELHPbacQ`LPpbGHKEJ(W|G13(}EtOwQufV(ve8bxVo
zQ}t(Ra*E2be;NP`Jm8+d_Z;Eq=K4%tibhZ0faZ_v-(+s_(pgenUHRYj|Gxg31nw;X
zz&QJVxBXuw7?xHqE$+2$?u)O5v#T2bfIROp^D8&!KRD|i6Fs;83s?QYuJ;dkFYo<>
zZT`j!f9d><kN(opdM158Dii=f!?O7g{QN)gzwUdl4!|W{-E6I`JZUtYtleIjyIRq_
z{NssxU(u+$I$7H~SpCuZ7yUme|GU%wadoruydN6?NIN-ux!T&iaHILJC*XO+%?<qX
zE>1Q6@-c{H@E_SfM*opHWdnfF82~_H|B*dU27ubP0Pt|`9~pfX0N@4#K+TASyQ|0F
zmhRmJ08jt}zy*i^Qh*Ae1DF90fCmr+L;wju7I+G%0a}1QU;<bGF8~DK3cLb*f!9C?
z5COyj?|~E`6UYGyffAqss0EsUuRu320DJ=`fq7sB*aUWgBj6Ws4FZ8sLD(Qd5GjZn
z!~kLg@qmOtVjx+NGUypdA7loy0Xc%)LB60MPy{F*lmhw)Dg>2->Oo&YeV|d$ENB(9
z4LSi`fgxZlFcFv%%n0TJ3xOrUN?<LpG1vy|0`>+6gQLM8z}eu>;2Q8(@BnxcyaL_@
zpFsc!CWHh+3*mqWL8KvS5Ce!c#1-NXK|&HC*^m-Q1EdEs4q1WhK`v2XD1<09C>$sv
zC<-XrC>AI#DE=ssC@ClfC^aaZDB~!rC`Ty2q1aFgC>vB5stDDC+CX1HL!gP!JZLqv
z3pxqig#LoTU?eamm>^66rVq1&`M@G!>9A5*8*CJ|4m(3dMI}XLLls3;Lp4WrLk&Sq
zM*WQ1iaLh6iF$#CiAIgagC>V&faZwy1}zb-2(1Ne3~dYT8XXUv0bLkf4c!vm3q1xs
z2fYD(1bq|z8Ur7L2}2Y^3&RfMHO2>w5{xd41&mWnOiVcDV@wT9Tg*Vr513_`y_l<*
zmskW?tXR@mMp*7x(OCIdZCJBdr`XuojMx&`2H0-c(b$F99oS3Q7dV7CoH&X&mN)@8
zsW>$_V>k!6=(r5HlDH<g-nfalmAE6gdwA%0jCj&`=6L>isd)8x(|Bk2g!tU}YWN8J
zNc<xFKKvg9s054zvILd{!33WOItbPYA%t{<(u5X-L4-MkU4)xNFd`-*MIt+*2%;}U
z!$e2K1jKyAI>esDsl?61D<lvSMiNC52a*_)3X&<3s|S=1Bp+Bk2!BxWVC=#9L-L0b
z53L>|AC^9xe0W7lLn=#ZPZ~#BOS(vgLdHs_LH3F)gRF~ekDQ2HgxrE0NnSxdM**f_
zrFcf+L-C1Xkm8Jzic+4^nKGHOgL03Ggi3<Sjw*qwg=(9ckXnq|mO6pDm3oJUm`0Mu
zf#w5EC(RKp1+4<DJMBl>VcKgr6I>hq23`tZq{E;SqO+z;plhc)qNk=;q4%XPrk`U#
zXAoj|!H~qz!*I^X#Hh;{%2><zgXtlYB9k{$G1CGwHnTXh3v)K}1PhEsh{c{Io#h)V
zn3bRP1#2qn2pgD9fX$XIoo$pI%Kn($kv*GzngfeNlEZ_em}8Zbgj1RG4QDOq0T(@&
z0aqMX57+NU{Er+S<vg0_Cg6U`{f4`N`-F#;$ATx7XMz`-SB^J;w~qIikB!fYFN1G}
zpO9aLKa{_n|5iXyz*V3`U|W!0&`dB*a9W5+NL?sWsP{4SvDD*$$4!qfg$0D&ge!y(
zML0z4MT$gzJYjlb^(5!Xx+q-KT=b*piWrTUnOK(Cia4#fxp=ntngpGMr9_^@mL!X$
zo#bc91F1(+u2R)f=hDK`e$s6+5E*$Hq|C4^fvmP{itLgcT+T+WSnf!kU*1Q)O#!N)
zq!6bttw^b8p;)MR^i=Sv-_uSdEG12)6s0v~R%KV^1{JW%Q<Zlr3#tsNj;gh4fSRJ(
zJGDi1CUqC}Mh%#Tx<;zT_h;PCe4q7b5^9=i7HR&{lG2LNn%8F5cGqs#!PPO+DbzXF
zmC=pYUD4yx^V1v9C)2mpuQNb3&@sp}I5U(rd~djE#BUU0G-b?W>}lL*LT2J%(rk)r
zYHnI-1~t<)D>Az=S2NExKYOnDJni|Bg|tPI#jd4<WrF3lm8ey`)t0rWb-eWt8&R8g
zHrp@6UnIWRvz4|@u|2j^u*<YNw^y^zx4(1HbtpxkA<Pi<js%YOj-5_aPM%Jq&TP&h
z&MPj$F7I8AT$Nn&UV>g4zpQg3x(|wj?#%9Q-B&%tJW@R_Ug^B5^u+gc^c?hJ^$Pde
z@|N??^+EBm@M-s@^?mKT>L=;<(I4dh+`m15E+9ByGf+OT;5Ei;``1Hn9=(Zw^D{^{
zs3Dj#I52qqt^C{K5S$R#km*p7(9AHDFx#-<aNh9b@ZU&FWPik?h{TASNQ=n+DDJ2a
zQFqbS(ZewUF&VM2SjX6@IPtiGc>H+p`1N-x?`jfg5|9Z$-<!VgOXN$;Ou|TVPg?n)
z{Gl$HJ~=k|CgnxSWU5qZX&OaZc-ncoW%^i#L`F#_WoAU?RhDhm%t!f;HQ9{WiJxGf
zJU)HT(aGt_70NBlBg>1(yUlmZUnzK2&{-&0SXe|+6jKZ?b}!!kZ1nlt7nv`0B^)Im
zOG!#2%78NWvYm3X@~H}?iuOvO%91LEs<djtYGe(l#;fM2_C@VVoqpYD{nPsP29bv9
zM$X3kCfcUdX5!|U7L1nQmb+H(*0VO}w!N=5U)S4B+vhv<Iwm?bI)}TIx_Z0ix;uI#
zds=%%dz<=%`|A6J`fCRS25JWR2djtphpLAKhHFLyN9w*k{?;)1WVB^WeC+GE%y{>N
z;>5tD>g4E@_SDR@(e%oU<;;&+huPyfkGbplzy;_6auIJaafxCndzodqbVXpLc~xe0
za7}Y<e%)eycjM*8&1Ud-?C*(N)LVr=cz!f(%WZ$#G2GeQb>6++3)#otPdi{fs5+E5
z96ZuH+B|kWzB`FLB|Xjm$@jDUO#N)>7vk6TIr4(^qTo{Kvgb<oYU}#d4eHH@TbA2~
z-%7t1?hto(ck_TW0KSjUAkZHU{zE8`KOBmJ0zrX6?}OQ&f{KoIUr=GN`w)eW@rUja
z7A6MPpUPbufQbh1f=oaleBiDVpt#RmAb;8XJI%P~FjO=s3OWRQ-ztU&fKcuaFa&mg
z3-tR$5)6W%Kw$tX8a@FbEjkhRLpYso(gZPyq{ahM9z7|SI0il;X>*Uns`1IYIRNcX
z7Ap>1U*@IQ{&m`O`dU9?RW4mEZhlOhu<K#DK5rNxP<vtZv(9`s@&S?JH}kd!qeL3&
zMp8tk+r_O>%NHJP0l$oO?*NT2DtzOAPW~GOGx{D^-ryIU;BYqUYz1Vx535i19k7#e
z)wS^OtF6GuT4sMF6Irs|9$H;Dc-OR>i$U*6$c7;uzpj<gdc-VisfjVuPnoHSzwG@n
z*;g>~$l=Y2HLqmy#BfFF9nd=6QEofE?sIk#70!a%aBtryuVK!%6egdrF66GUv^h(J
zQNR&T_l`X`2pgm5XIP$HSo$3>Vv@2;cI}BZH!)(ZrH_pnx6$K{T`Q?sMN7MnPAodA
zX<Cz#doCTT4tK3|s+-2~(AE{6lW7De8G>2b%N_XyDLm*lPT0c6Xe%stnKL*Ef&2@k
za09>T$g&Ro34tJ0*`Oc7rXyd;p*>*m%3-wAT^SzkH+e?w*S&CLKc6?}+^2!3zK1RG
z@i0wcAt<h=sUfnMX)xE&_o-{^_n`&339x!hh3AuwXu<Du=QB@cvr8>si~ezg7H;P*
zY>8ua6}IWX*o}}$E2bTBj(ExpELI6Qba*IRyrA^z;|yJp6`THw)!dSj?xMZ;)#7>i
zXQ!7nH69F3^70rB+P=HY@y*iKyk@Eb8@fisxb*a@6H4=~M3dopO<oit5ft8j{-Un$
z$CpEPogGMo8c#~`SRpOYJ%9QF6Slo(w}wmmqucLJosgGa>2-Gic81?euqtMgXF_~>
z*d4&NlF(t8q*3sTC;RF+Ea6)KzI>nmHdh6k{M1_2dJGr66l39UxIVot=5)}={vsqh
z^OOiFdIwBv9}H^Vs9oIl&t3VexS~=n?tkoGZ3aWDRWXJ*@*X>+3w>2Yyy5m1$-=4L
zp=F8)UVa{CylZHi*BXu^H2Bgz4}A<zo8+|y(InvyWKsCst?2X@OEF&aCt2OU=^95{
zFSj0W6uMv!=XYg|>o>leXxZcW9K%u@EzRBs$A|m+`aYSS^^uusy`N}<1C`1mK+a3X
z;gg0rgQ-Z_c}&}(9!s>OkFtWqP-5L#(c3eMj~lfVC#fH14!B!C7)J&7XcWYnZQ@gT
zotkC^Vj^XVZ{tKJxdTtg^!b~f)DEm08t6Q~jNFqZ=qWbcu36I+#@}dOv*e(Q9wbF^
z5-JEB;K+aR#&}Qwk|yMFZg$|9CtpiSwY27>uC*{>IeXL-gwxF*h+h>|zY`nr=8c%y
zQ)EVQz(j*$H4~%f>N~^kp5vk4?Vp(jJq^M!NHp=MGhK9Z0FaATxbBaz^qul~r!Wfa
zsm1Miwb0{6=M8KzgD&E;6J73;_YrX15&76VVC^b!mM1F@DW+09gwA9gUYst-gl}ZM
z%%%A3w^Y=&iazyGAypQI6}>yCRc&ZOmU!EcK>7{{9jKbrY?;lS8a^?!ed%rQC970d
zJG?qesT3VwrAmaZX*LM!^G37@mGC_c%aEiye6<~_O-=4BqGaNzy(GQH=~X+^vSe{y
zHS_*$ljkQgA5Hp-pW)BwZ^p;y*-C0yJ5ZRTVhUCt#GyeeQyuB1DAF4GsH+pFSSmiS
zM__*c)FQFasQ-m6dp4U(?z5CxOqY+0{(jOVp$11<Oy8597;>wI`!60_mvdF@vB|+=
zt7BxVnlT&_uNMWPww+2v#Or-#{z)BzL!Fe2&xIP*bUOUveo`CSH_a@mrCTP;`2^@!
zz5l7D%=<D~UD|hwzpz4X0P9VU46pBGUaxkxjI1X#fSl8!xqn!weIs$v8TVz0!*Nq{
zw{yQ?NLECqn7E)Rt6^lhsLX8CM6AG~)-x_5lzsJ{-Yy);j?7RJ>_*Q$th_hYlf&d}
z2Fk%H)-wI|0(Obq8$Hq19SX`%)`$356)93JW!{5)=A;S~g4yO+eK^Z&tcE$p7DS?D
zwhT=_Os^Tf>Op&Frdsz-osg=x0xg!nttMhZOYUP2f$e*V?P&_$A74RFa+};{*|*}(
z8d-*iB{p1bYoP)yG)wDs7vB79&xXPe_?Ycft1=3BIN{W5bYO;aAKCGSfn{8BC9H~p
zIFTVZe6X&1d0^*zU=&I}zmmr8^PYXbxA=5%Os&QZx%(Rpr2>tw#OQCrh6Fnt1UfvO
zkzVt*7Am^BPKd~>(c#+*{5-ujaTVQu3JqTwMZzgwdRb=;Vx7*5u>Q=qmT4~+3q~*u
z#}CgqluFCQ#6BD9ofPNtg*D%{pS4XA2b{$DO=ucek6yai)At3=GhRE8@(@tz%+pC~
z%h1V~;veobb15E<?TX#zJVsm^m)`*;PItgfl!}0e$E$%S{ia0=pLreavue@O5PR}2
zi*>Ojv=}t+7@)y;+q#t}Chzfu&-a-O!!xuI8FXxXjYV~-WK5xjWFhggk)KZJiH>3(
zVt1TtX&)1Bi)jCnoizP*G}qeSymCw28(@<^-(hzLu$~9~)TO=6$aKOkzA0Y99W@I5
z<?YubP<m48DK6_uJ-hwC_xW<;al)x0i3xv8Vu&@2TnI_d#O^%|{P2;cZagexHs60u
zAvj{MAmX#=anps}_l4SdZz>%Z#?FGrgcTJeyNE+`H=LOu`<=%++X0<8ZY^l0uZFaA
z59++@f0Gz3CLd$1;!yVOHaC}b?eCL%P~OVS<b4w~C+)J9e3X|N_Em^JOL4jOdGSn1
zl|Z&neUs<+LCvnlKJLNnzCu?&*j!in<C=rlo`oKW4%>Fnw)liX__E@%eF5tK_rU*G
zXYD)~Z1Gs!H~u}cmp7pP^ta1WF)5}g>oVDgON)tPx}5MS+%^Y|Q}=R#R2Vwyu`xBj
z41tQ?jbq)SZ(mw+^SEPo54Cl>(DOunF`-3=%b`!;($qM0`L9(zRS3nJA<dR94HV>{
z1x3%NCp>f-(ug*bwsvMX@(ys*QbX9J_;Z<v%eOQsHF_*n;w_ST)bb@Ggi#aZ_UDhv
zLuWmrS_sg$q|1+K;?q5{{7ip7^cRVG&OgI%Wyqqd3;N)^M*czqs|pN=a@Ce01IA4&
zg}#fDj{1vYF%CW`dxS*GZe1iMoL5&VTMrSfa+?pLrfMut6>+W@U=LYxC}bWlWYjeK
zR)r%FmW9nX&H>K1gwJu=H4G%<Ogrg_b}5sG&L%{@QQcm!EN4i&@sK(Hs4E-ifHR|{
z|AKkkmNGsul23jH_0zXP_b-@D*G}&b<rzGNX*F7@iRXQEPZ~vZdpa3si5@@tC}VS#
z`9{(29P@`5Hf}@8ne))Q6X(UGA3}#YqtrRQm!v<_f3n<a`-i0z)Qwk8<@a2xFR}7Q
z1=PoX=B!k}C1B}s^=8iDG0k6*m7@!AYMZn*?O5zzI(V?&B)VX{F!<}Y+;U%pn!q86
zHGZ@j42?#{#UBwm?<z7x$GHBKxh0F-63KF&BXrnt)$`DYYT%P=v2fW?0eqo@V(mcG
zH{7YS;?h*VRKX(HqggyYYMTv##La6CzkKitr)RGPwI{tWclJwV;GVY~bJ2uW{8htg
zj-0B{X$e<jblqxdi2Z?_W$U!Mqy2u1TC1TgYBV2FJPwZ4E9}qNv#5D&J9NEv%B}F<
zB~P$#>cAdFr7T>FU}K}Ou<$V^g_KLpri&7~en*j!IY*2GIRkg2k4rEoUMw;b_VUjx
zS<M`=&F@5~#cr6b2L%;L-DEplRsjZrqp83R8x}9pD#QJT3uBm1%r`Sbi}GmeXB5wa
zE&BqWm(+rVQjvi<cR;@q9hg88u!ZC@a42|=_qRwGoH&ur;z(HjMCF{fX5+z1X{agP
zoG<FNZm_3n;btL>hWCV<t$solWdz1EfAri@>3-_EKEnRW$AarM-PeUDX|<Tpc>8Zd
zmtmtbRgtVM{kuGR@(E%sX>IhXao*|2yb)&$2AwA<)+gtOdzF;!4FgV&^j#z>+d2e_
z6HcAl{<|u4${6Pd?Ah0^#w6>x6nuj=YrZ+WXY+5uZs@%OI@#wGe1iF(vxiN(ga>>o
zRg86_x%CsWF5W9#3Qxvih?KYAgV$W+)~K{}v<2gUn4Z$#|FjB<kP&B-_P>}{Yp<2(
z#ebul-T~_+2g}D>W*!-<gQ4~w;&|Qp4ee!nJ6Lw?f>F6z+E9U10u5~<lg8aDwf(~A
z0>6PRG2WolOoWiQ8dBWGCat|L3E_rw(UbpBy~UE$)as1d9Zk2R>KH0zdWOMqoR?zH
zcp{R}R5mlaXX9GV`s3SoCwCpvNHbSjtnkMw%K7>)_AC_4*90;^E|*lWIW$4+PTFlh
z6?+aVe+lc^Q2tuS@t|wMIHlr_f?o1vm9}Cpi5Z73bnaeU+%b!l!c&}>CTrXrE%&zL
zF;6l4u>BHB{n{BE5rXXaL`?Oaq7|N~O5%_%6mpPSHXP2FJ-bC;Qwm;e(lpQAU{@Zj
zMCZ&k5yX%nDnky)?4BAswaiD)AmMc_{c_;+IJli(US$spSLNd7a->>ls%E{Q-Y06B
z1I_w+F>XGKk&;5Ez7H;BZozT+4E_7Q4Tqz_cYvvuv9{+Dbyr%g^KMN^#9Dp}j)_s-
zshfuJLcAqMRK~4F*l8TJ&{CAOh)1p&UOLKedj}*{W*%n7Q5Qn45#j^0Bj+Q%F%z$E
zYce<Ywyw7}<aZ=~&ufLa;98cE&mJE<Zth1dD6IVO!CyK_aslm7cfwdv=WW55y5ZTp
z+GG5AT~P)Q(^sN65v?)upiUAV_ALj&Afg5RNNXQmn7g0KZ&~t-aQ*gF${H(vZoD}r
zw8=!{GnBU)-PhadVLCoFbyV8n_2w_r4wU!er@qbW4$dfiTMMN1tR{O@NGE;%URxom
z-Smts8-u?1Uu#UhalMhXq3~1Ivq8>|5Q_Mba>MR;3lF_+a!DKu%{8u^T<*2_Sx-W_
zC;?_<h9b=*;Xz4_ULZyx8l3f3n{c-Tx-`D;|5Vp1c7Kz7v(LPe5U%eY_4r0r#Kh#C
z+9k$S7M0B!!yHqx`<bzS=4#&J9UvP0BY;lp4&XsS34XBLF@ISei7blz)a>O6kMmo)
zx*oWlHfH(hR#UO2R#^~1dv*T)2d-FKa<_}Dq)N5c9=c4Jk7`o3))i;7!Kq|!{^ul4
zKZMqc!x1mav{QDpv`{=0ydez{9!Y$v>FOUEV&SE#DR!cbqymN^`(M6mzBX6ec%^P{
z>`{c|4MY`?Vp?$O?bX0nDw(eOby(`>`*WnyrsjN2r87lwo}<ob1Bbs(zze1R=cr+v
zyMU~vq`x{%AFTxRSPe@5bB?6bM^a|Gat*V>5b>MsyS~nf<H3fBsrj{h+yP84VLlF_
zfrBO2Hk%(fq3|r#HujVhb#p8=0;{*TX_DVpDrc7K?J|V7L}~RGDMux-9>InHjumP1
z`##m!PJTwoid5QVMM~#bZ>KI|715-eN!)tA`OBK0zN&@!3z+&wZx_F`M_FWADzQvM
zQ>s75!B?xJwzebDeyn2)^wR`>sA#?{jlZt5>LQCaxdUPn>s3!2u^gT<o&<0qPA8J>
z_a0;0lD(z!3qmPU;8ISdk`ne~Ex63OD2)HIS@-18W(lHybGbKdI6o++5-!3oG>XV-
z=E<v4sdz`mzLtE!*(aglD#e#qpdjjr#Hs4~jJsCtYByzhzkubXn7a%s@f4npE%E{d
zv@0&ZXdTtyk1EG~{feaKs{YCsH}*+LZaAZF_?U&ZQ}v+7V!}0YeFQJ40J)mXg3r}Z
z^sC`?`G(H(!*4}gdmdFAS4u<)tvDx?C|cf_KHO!7adcX;)wH(_z581@VU8x%@4~_W
z^cHL*h82Oc#S<orLf7Qv_qRNgBP|0bvjSsc`o^eeG8+U1g?Om;%XUFiRW<fLPLXTx
zN?mTA8vpz;$Rd;f=2WGTV<hVhc@7O_(RCavoXP%Hc#^z*BuwBzzaz_%tF&-;5+pg;
z&{*41mEWA7`iUdLq68~3%>~7GlpNn39f~j8xocFy?O-*0<}%{^%J8M-yv3|M{JCPX
z=y#=d5AkO=QW3h?t@PK?@*xsy+4+jJTS-?t>EEZmZWx9lj|Wy8tLoh|j8j(XySWl9
z^>0(n;<<MWETml9l_Iks&lo&?D^)Py^*cfZ88kBS!&<JDYb>?;Ki|wN2iR~oS+O`d
z!aGL5{l*286IH>tL3aT2E`kpyxq_oALniEq<$=EhPyJ2lMk|R>hqzyTH3ZtyP@Sv`
z0Xl!xQ1l)e%#P67iFy%Ds}T24JEIN><AHcNM@5T`V%||d#PUfJ!hKj;c(J8(u)0<=
z{mlqbT%T*ee;^~jCaiF3Qs_B=`yf~bE%05ZkBWD_CZeAFOYaXuV%|W*aGTmQ(sY%<
zVsEmYAJwGB>y8zOp-fsUx(*mErendC+<Fg%$aPd!#`k}3QvJRA^%vKcTRVCeWmw*v
zun;-&$hPvOvKtmNmto(Yp74sZj4=tQ5Ip~UCGcYwlJ@Mpsd}>4^Glh&<?Yw1K7sG-
zv*h_QL!6$whYd>hCgV6|3(v2YL%{RH`6?IJb+W&654BM{*~07PE{<QqR`W(y5zQfL
z!^g<IDUtEFU@9iCv@89a;@pjshq89!NH2Wv#Zue3owa$BM|+9mo1q!27hZP&g+W%a
z_1aj!%r<v%TT=)*r=oLT2)?(eGLft{zBaLj@Xp?FKO&U+k$}g5&9Q*3EonNZnUqzB
z*v8Y?jiiNVJaD*Ry7Q#A?m}y=VxdloJwJlVr74M;fMVBy>*e<YmqAlUTD3Cd`_Dn{
zH6#K&{;DFd*E_h85B*gyNxef>lBB|zQ*>2)&lgKW4{6CSKA3jQrF4_8AJV|5Hpo|k
zcbEp|HS3tm=W#jjfPKc_`uX0f38Vq7?yBd*ie$HY&wZZ=xGCtog?=k~$HUXmfm^^T
zU_BFC7?^}C_=r&i&$m6BVxMzF^;Ipq1j(kpY(;F1E^OlOSI%l(%fcy#PB8cDjy`vl
z^EpJ-jK*mpWNWVQh$Ll{!~-t4lAV6ZH=j?>S9PG*%}d_A4Jsq<{1jET{tQDaRO~)~
zAs;h&F-h2SAm`vWy==UHx~8MnY{*cf*OlR9db3$W$`4a9wud(%v38%X9Pzm)aB<sC
z2*wFBFiv#7>yys?fhc`&K(_HdhD)%PG1Mjbanv2aY2j@5iO`vQC)bPGtHw`O(|v1^
zq;EQ+r^>brkBGY2b^G5Nw102Y{=<t-lk}nAnYt*4`^!^xSM+rdXM;tTqeki;BKu!G
zzqZkz6A;K~X<f0}bNj+7ME0lG(x+)_sb;Sur@8uqd5n<u$Y)HELfoI^q)qcBzGX8h
Wq!-0;S;{EY{l|3vSJ36ZoBbc#ECSa6

literal 0
HcmV?d00001

diff --git a/doc/src/Eqs/pair_entropy2.tex b/doc/src/Eqs/pair_entropy2.tex
new file mode 100644
index 0000000000..c91004840f
--- /dev/null
+++ b/doc/src/Eqs/pair_entropy2.tex
@@ -0,0 +1,10 @@
+\documentclass[12pt]{article}
+
+\begin{document}
+\thispagestyle{empty}
+
+$$
+   g_m^i(r) = \frac{1}{4 \pi \rho r^2} \sum\limits_{j} \frac{1}{\sqrt{2 \pi \sigma^2}} e^{-(r-r_{ij})^2/(2\sigma^2)} ,
+$$
+
+\end{document}
diff --git a/doc/src/Eqs/pair_entropy3.jpg b/doc/src/Eqs/pair_entropy3.jpg
new file mode 100644
index 0000000000000000000000000000000000000000..0d5e0025b7200f9ca4255c48bd2bbabb291553a2
GIT binary patch
literal 5759
zcmb7GXH*o+w(iL=<R~CvkQ|kq<B)^EkfUT6GDDDzAVCo%X9URt3W_iwA`&Dh93)5v
zMFAzlA&P*!;hcNdx$E5@Z@s;`s=nIwb#>LQ?%v<U^u;1TsiUE-0f0aN*aB<;02kko
zC~Y+QmLbMKLt9V%asWV50t5pSYHy!F97a<WWnpQBy7c|q9Gn9DR16LE{@wrY_Gkhi
z%me`Ag8x1Ce|nJ4E&)yi))}EhIQiiM0RR+2z`VBu{Vwri0;YBJ{0lc-VjN+C1pmM#
zcKaL8|Hb(m@BYOx!JrAdf&u^$soQ^G$N#|pdM<$uTvNdXVqINAP{zKlfwvuSE~tRZ
zNDyL$GQ|11V!d20S^u*C7x~{e{jWCAC4{gy06_cth2pSo?t!TPT7l>d5fR{@v)DKK
zE5;!Dp?`e;O!|+{w+sNJ4*>v@_K(jo9{^h70f1rlA0Ovq0HBBi0Q^U%AYAa@GZDfC
z05AXvPyn<5Bftu90K9+@APPtUa)2_R4(I>|fC*p;*a1#}JKzoAfZIR>5CbFt4}dfv
z8z=xufC``ns0SK>7N8yI0{VbKU<8-|=72?D1=s*~fg|7)1Og#IWFQ(4BZv*e1rh*>
zf}}x;Aa#%)2m`VN*@N6bKA<2_1n4g40Vo4h04fF5fL?%Jfx1BNL8G7<&{xnJXb<!Y
z3;~mZX~8UDZt!)mG*|_!3pN4Ug5ALW;4p9;_#rqCTn2s$#)G@SgWw79B6tIQ2mv4@
z5PAqZL<k}cK|>56RuESR4iW`Pg5*HTAoY+oNIzs8vIyCPoIv4F8Yl`X1eJqoK~14f
zP=9C?G#Od|eFDWp-$BQrU!l9u-!L*5GfV&`3)6;Kz}#TBVF|DtSQV@p)(e|}t-y}p
za5z1j2QCHIf?L8p;NkFO_#=1&yc0eOUxFVZ5C}$u073y_fN($rA`%e!h^L6Rh%v+p
z;wKRa5gU;xkp__!kq=QUQ4UcpQ3uf&(GQ|iVoG8zVp(DXVrSw|;#A@a;x^)s#4E(7
zNGc=`QUPg#^gzZS^N{t(UgSJ-pM-?u8i@>v5eb&$4oM!#bCLm)uOug=)TI2RXi^)}
zAkv4V)uf%IGo<@uWMtfA%4Ak#fn*QKYRG!XK9l_<ry&<6*Cux+k0vi9$CHndZ&MIc
za8amG*inR0<WRh%_(-uyNleL2iKcX*jG`>0d__4$c}PV=B|>FD<xQ1LRZI1rYLyy6
z%}uRN?Mxj<T~6IYy+i|{;h;g&IMLjtsi5hlS)ql~^3ZD2deA<geMUP>yGut+Cr)QZ
z7eZG^_l9nf9zxGeuTAeopGx0IKS}@V3d<FhD=t^=U#Y(`cIBAiDuXhE3quk^1H%Nv
zFGdujI-@6JI%5mt0uz*ppUH^nHd8TEFVp5#+N*L`ovtQbZM-_m3})tM#xRF7mopDB
zAF{BrXtMaR<g@g!Y_ig`Dzkd9X0dj#uCdXuDY9YNve@3Tt)u8rDkv{h9_k%xmz|ki
zi#>?Flzo`}^cv4K(`&KU8m=vHAUUKtTsg8hx;b_^**Nt%BRK0gXSs;Eq`BO=a=H4s
zj=6cb&AAh~Tew$w7<jaK!g%U<=6T6@m3aMm%XmNW!TF^5Joz5+jqro{CHUR>3;90^
zfCVH4umX<+Mg?JlGJ-yWWr9;eq(Uk}!9sOHUxn#~^@L-ETZFf-b6&T)o_@Xm`tKVO
zH@t3C+?W%g7SR!j6=@UM6Xh3m5-k+{Bt|BtAr>vxD)y(ifVhizvG}wEjfB2LqC~gE
zxuldNPO@HdO^Q>>UaClHN}5*MP&!$9Kn5nGDibZ!E^{I)DH|x;D7zykB<CqtE4M1o
zBkv+#A-|+>O~FB-RAEsOrD(7CSaDH_UCBYIOzE34hqAMBmGTc2J{1p@rz+d3H&k(|
z&8o*}S#%`2Qw^e~sg|TRtWK?Ns$QV}S>u|9n?{|+uBL=$xMrsoOiNcQU2964McYZc
zMtfIBN+(k1oi3@av2KCxH$8qmoL-wgSYJm!Q-A&@*G-?BEe3#rwn3)Bf+3Hgzu{{m
zxRIezq0x7Y2qps4Z%ku+%edC~$VAm7)nv|;*EGnq+l<1@#;nHd*j&v#!+g==x<#bL
zpe2(f*0R+KVP$4jWp!w+Zk=ttVk2RbU^8jUV;f@o-j2!6%kIrBidzo1UfRR#&FyRL
z&m0UK${Y?IwH=Ebcb(8qxlS9-%FbENYc2{d87@Cu6<jl1SKSocGTqkPmE3dOH?e4J
z0rpQ1Esqk9V^0ImYR?NVbFT((B5!-|S3cA}o<8q<S$#u%NBso+68skZW&N}LcX7J7
zssK=cZ9r=vEg>lm1@Q*O2Yn4z3@!}*dE4yv%Mhv%pOB$Y{?PlOKf*M^D#M}SPT}1V
z>=7{$Un5l_OQS$hj#1sw9MN&nD|a;SJc&WZc*gt@dm}a@_TaAh-L^QExI1x6@tW~<
z3FHa5gsDWi#FBf^d)Rx!_r>n#-~XNDob=(rjR!do&XS#yKctAH<fU9ZbbUCSDw$fG
z22b-zn@m?qf099!5tgx(sh^3@LS-do9c16iexD<rQ<96!4a!~2)5~kk=gd#bKPzxA
zm?%^&Y$#$bx?gns$obJ&v2t;J2}{X?l3$OpkEcsDOYvpgWjW>W^5F9C6=oIvmC}_p
zRZLY6s?Mu@su!PNp1iA(s;Q}Eu1&22*9FzBJ+*l{@=WbnYrRl?X#;)3gXh5Wpy%r^
z>|aba>Na-2lzv&?#MM;POw*i%2jN5UyDjc5i>;Qeqpx&cb+^g2HNO^qUER*!UigOY
zO=<_SBd+7(ZP?qxPQT90F1N0w?pxh+Jr+Hm-Wk0c?$zxb=+o%y=~wB0JD@Po{$BR|
zs}Is2S_UNt@k0_r&3{Pz(L5|UjQ=S0v2{ddq-|7wv|~(ZtZQ6tyzi6tr@@Jv6Qh%+
zlhac+Q;XBi)2lOHGk?wo&z{cRnTO3MEl@7xd}jV!_J!|D!=mKko3CnLhrStqn_F^P
z+E@-)K3j?VPWC<L2iuRDRngV<HI22Ab?fz&4Zn@k&4ewgt)gw-?WP^&ouOTe-IYIa
ze_re*?=$XKABZ1x9~vHhKJq>~Jx=<``18q$^hy6O^ItzsZ=WH~^3M6r+kWf*p1<(E
zxVV@D&;Xc_pFyC@8hmL`$Yl+KLLpE%jF8MO3xb%4P!Mo9Aw>}*FO7gmNsy$Mor_L@
zgovOTazG#|;Nl&?LQrCm%j+K_g2E9H7zj*E7?z*}KoAHN1cDKPL9j~>2n>P3;0S<-
zikh8<nD&~83f&b^BMXPz2{BcFj+of^^ov=5=u*8Z0jDSB8YxWTEmI1W^j<;}3*l1u
z;FZjJlF5-u^;&~ez3Oo*S6fBRgq2>CnxTBs<CSJUzMu3+)2^j_Mv`bLSJ7q{J9CPc
zJBsfU$r_jAjJ(%PRd&;4IK@~iTIf32)-!^h*~>3D`sk8cTP`sUq^%&=D<ll7!#8n*
zTcu8eKKKayI^7J3%3d}XdaW)vMexNQ*%#cY{%ULeQ9u6aOLX5v+qrPnyo7MMuX=DP
zJ!ur@F(0?PkIv+@YFL#l!iDxpj;GOjbC@YPay_W6D8#>zk0hOM86TI;%Y~{@1AVv$
z+2<tF8m+Zk6pO5Gd1`C)eky~zUXc4Enej9_Rmj$)3`1Y>7+NaeuVa;wRoP6f$#gt7
z(JQ>$B^kLSxtHr#|H8YcSRi1HsQRtM^~PMxMtb`0^Rn!R)@*S<C;BzogAIJzABOx$
z2^}rh$lCD;o)dnMFj(oE=xe$j^fvFglc=ZJA6Atp5$Xm#q8-Dsn|^5ule*d{x)ATq
z_=nOwdgNbluV?1qR*Yrp>8ALiwr2xDrgMjVkurgMr+j|T*Q1gQd??O5n;l`P_qTMi
zHOOR=Ou6|hC-{kzy4JIWbfeI32y#Y8z>=6b!{Fh~Dyu08>U})lRLUzMgYnwuR{5~Q
zM*LuM>zyXvQhPr6io4x+IMzX9j0Rt-a2+kiQl6$H-&&G!5k)1kuxGw@N0>Q$9sRp#
zll!zMHOViG@>*;3(o#b_mo|&RhN{hHD$4lahFqVa_pc-S$Li&CgRBlhvhxG$BBw{O
zO<_1Yv@c)IV3n9p$657!P%nGO2R;L0U%klOu-Sa3k?LEMLct1VDdutWZMHg0rxJYe
z)s4(57l88<N|h5z6(*y|$JU*9P3Ks8Rc_Aho#6z;9d>XpPPd+inX#AjnfCWi#r<}U
zIh?!zl8;p=xIdFxZrLZ7s2XpzzQDYCT3w3PG|Bd>Hqp|S_d|U8IaPnK7i`mGZSkft
z81cN)<~WEosd}|Zn)3X+TKtigRe@7u<Za3s>-(PR>-62BKbx8yJQqqlqF4u~gNNsX
zeAIJA<8PLklz*mgTWirLudfuEdGWnVq%7gvTecA$WlU1W(t=qrJWtQG^-af2Q=u=_
z(8o~$3Cm0}{%oTt*;$X5(t^j1Ox=-J!dN@S$EMawJBYtRH=HN^aXy$-Z}Tcc)_0by
zj-?${+5Cl$Ls4=SZ13jo8<1%KqFUUzK9}J3#ARz|`GfZt$GKfs%gJ-uT<w~sYAPdL
zoY`_e-*$W}=382K&6v%4Tu!cpP+jrkoLhviS+B*rXhRCWGTR-><8_iyleD=LRAXY(
z{Gxp}OLw^LYoAk1U*@OoWK|n0gOlRb3|aBVuWe7Z<QylOD^0f);<ho%#p606?z<wi
zRrJ3)cS|s;X>s5PoI&cqdSTmH&8um<K^f~i?<K67Q_`wOY}2Vn$_7lZH#Q_R)vVMd
ze48+hPICOxeiVDHQ@y);8DsO(t;_NhVkFDDN#T9X*!mQkI^{H{_kOP#8uz!Ju`T%>
zVPe>>1(i50<Y~6JEG76B-91|2pw`s!k`?q2b(Ekx=)qBQDK<;5e+i;`<ll6rr6qt#
z(q61|nuqmg^GWEQM#-V(D_{B>3a__^hwb-B#Y44hgS{11*N`dtc|B`U9K)6TtFsy3
zkoN*=RX+uuS<=V9NqyR;o`}U5fwR@BM8Ns5mmi$^U7GN7_UEdYG^MTDs0|UWuA*IA
zb4IIi$Noo>K@vxg3M#ptb9~}>Z^4Bs)?qYcRdkLXZqj$PzB6YlC%quas_eQRm2J~+
z4#^G=?Zvlz<NLw!iY2N&uyr~-VcuNTM3#3?W{*tg`%4kYasPVtKSKLCbkk)ofJbnJ
zAHk(nUmJ4KCIj9#;))r(AFr#1VzQ7q&3tN)9IcIeWsXUGl-Rs>gjQZ_idoR>KRyrZ
zeZ@Za)^+k3k0U$nPi&DzsQ(ILNiQu<nsPR)ioU)&7X7gy|KzS;Xnl45IwP`D^O^Yz
zZq0!MB8F6d!vEE=8a$J<e9vD^3aj${=OK3$T3A~vkS+9QY%Ey(n`zaemso3UP2aR#
zx$XPg`n;1qlHL`PwrWbh{NGhia2-v%O^8^#L(^K49lolR8Ao(6)=tbzUXRJyEbunt
z60JDMSrW6Zt`$tPNLwP|$KIk=l+5l|)LJ|bO58TTG371z#;{a*Vdca_a$JxKUAlOm
zZeVC{TE5B=CCkizrgmFFd#Rt<BVThKHhgZ!LGoqTAW|S>R#(_N!C2dm4)p2iIH0U_
ztzMIPzgjP%BdpIO-!0sz7UwjE5AgbsG`FFL6?S<(TkORf_q2R-X(YXC<NZ>xjO^m+
zmoEmZ^^?CH4G!Brh?>#NE7}hEiAk!=6ouVAIL0=#_`GiAiR#DcEBv^X$nN)h_+apv
z?0ONR6O%DxXync8aRG#%uAS4(ei_pK-J<#V_pfO6;Ff*r=w}mQ#ow3LB#yS*%wL8x
z_0@_NKmAm`d~YBvri&Q?6L`L>wHu{r#`fm?;CIrN?|04%z`A@VI#Q|Zwe`t;rODO*
zT$$qW@%S%+=C^s7=RReVqOAMd47c3Ok9b^fJ~YeD;dZH7&~z@wgw^J0qD~NNPjXzs
zgCuEmsx1{i7J~26qiamEjP?RV#>xN#$DB1R5*5Py+r~2j7ooL4BdM$R47<$k6SGH%
z$H|`q7r=I7A;dZ4>Qz<z)&?%$RPwx#PgLkC#?<8pUtl6e+uK`{r$#74*Cq0eXs7`c
zQdxN0Y@xCk$Fh~PNirzu=r!t8rT)zI;|q=t%I(wch%y9+0*gOCRpg&PMnANp5Xd4`
z(}6(074FSaM^BH;XwB+oKQYxLVJOBz*Z<8~QzmCq)=#4Abm{#+>bMDN#C?pU{d{&?
z3TwU6B-*6+S~M8uaa=z0%533DXF^qtn5Uq^cBW+{Y326<mmWl$Ni$@}F6yoxt47&M
zVSl*U!b_KJ(Y2PXHrwJu*G)Y3c}-o4h^w()w~GJIdim*qr`;X(<Y*L6L&oh$<JwnI
zWC98!yRxed!umd+YRs=bK`5iO60gx9-K!sJ8|7tC$df1TPxR{|nmqkq7t9Gr|K991
z7NXx8Ocb!ndE-DQ>$%5I;!|$oBeqhBpdPmuZ$R={E|);M#+Ax7GH-dZy-a}r&UKAU
zaSS^qYv=l9&itU(TUOrUW>{2(>b(Z`n?h5Nn_!HpnhUw6JxTUodd&TkT*9}OEq=km
zoQ?XPD&tad+kC?(I%vTgG5eUwM(Xs_%)YEyADVZ?V<m3E#?57R;<In0B!}yDA6Qlm
zqSq4E0>-LjX=6jIh1BCZ<SC5~rVuZJPyPt)@osQ?^ucP3&wSKc(gj<C{uUGnv|jI!
zIpLb=__gyp%zZEAbz1!c?GG9^2DWCiUF74XEQ0SXD6<1+S_d{(m~aZ_{Ux20sYI#R
zloN|xme0u&AMFHQ(f6Xa`%q-!1@R0*Zw;ot+9sui2MnzBpR$OD#EuQp$5ZPGzRI0#
zuX-)ejvbmyF-8joLB<c}th(D{!)ZPiD;Ah4v{>@6%;rVYVtgJPrQFPVG+-f`%@kx*
zRYZENaV%%Yl$Z6sj)_i(GG|5vcBEorO=3MjB8d(oKb=7}bKi(fH>I`Y4}NV_>b%FL
zq;)6Au@t5hi?Y#eey2~?v6)PPGQ{g=+R0p#Ma=u%8IM(%WA>T8wg_3<eXK6fFZb@O
zR2RNN>Jcb@W5#kRJ4hLl*|~q}VdMV70dIAm5uagJnT5ZTW<rnyOY7R4B{j0k1T9m+
z=B|E^@V-;cpRt{1(^%hZlS^e-^6c}gS(rXbG|kBA!KPy>8c$;&88kCko>ZRxo%*p%
lt%?!s(CDg3>u#B*rrEW34np|ZXwFI$pyW`)B<^`J^FKD8lB@s#

literal 0
HcmV?d00001

diff --git a/doc/src/Eqs/pair_entropy3.tex b/doc/src/Eqs/pair_entropy3.tex
new file mode 100644
index 0000000000..b0b3c76358
--- /dev/null
+++ b/doc/src/Eqs/pair_entropy3.tex
@@ -0,0 +1,10 @@
+\documentclass[12pt]{article}
+
+\begin{document}
+\thispagestyle{empty}
+
+$$
+  \bar{s}_S^i  = \frac{\sum_j s_S^j + s_S^i}{N + 1} ,
+$$
+
+\end{document}
diff --git a/doc/src/compute_pair_entropy_atom.txt b/doc/src/compute_pair_entropy_atom.txt
new file mode 100644
index 0000000000..bed36ed4b6
--- /dev/null
+++ b/doc/src/compute_pair_entropy_atom.txt
@@ -0,0 +1,121 @@
+"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
+
+compute pentropy/atom command :h3
+
+[Syntax:]
+
+compute ID group-ID pentropy/atom sigma cutoff  :l
+compute ID group-ID pentropy/atom sigma cutoff avg yes/no cutoff2 ... :pre
+
+ID, group-ID are documented in "compute"_compute.html command :l
+pentropy/atom = style name of this compute command :l
+sigma = width of gaussians used in the g(r) smoothening :l
+cutoff = cutoff for the g(r) calculation :l
+one or more keyword/value pairs may be appended :l
+{avg} yes/no cutoff2
+  avg = {no} or {yes}
+    {no} = do not average the pair entropy over neighbors
+    {yes} = average the pair entropy over neighbors
+  cutoff2 = cutoff for the averaging over neighbors :pre
+:ule
+
+[Examples:]
+
+compute 1 all pentropy/atom 0.25 5. :pre
+compute 1 all pentropy/atom 0.25 5. avg yes 5. :pre
+
+[Description:]
+
+Define a computation that calculates the pair entropy fingerprint for
+each atom in the group. The fingerprint is useful to distinguish between
+ordered and disordered environments, for instance liquid and solid-like 
+environments, or glassy and crystalline-like environments. Some 
+ applications could be the identification of grain boundaries, a 
+ melt-solid interface, or a solid cluster emerging from the melt. 
+The advantage of this parameter over others is that no a priori 
+ information about the solid structure is required.
+
+This parameter for atom i is computed using the following formula from
+"(Piaggi)"_#Piaggi
+
+:c,image(Eqs/pair_entropy.jpg)
+
+where r is a distance, g(r) is the radial distribution function of atom 
+ i and rho is the density of the system. The g(r) computed for each 
+ atom i can be noisy and therefore it is smoothened using:
+
+:c,image(Eqs/pair_entropy2.jpg)
+
+where the sum in j goes through the neighbors of atom i, and sigma is a
+parameter to control the smoothening.
+
+The input parameters are {sigma} the smoothening parameter, and the
+  {cutoff} for the calculation of g(r). 
+
+If the keyword {avg} has the setting {yes}, then this compute also
+ averages the parameter over the neighbors  of atom i according to:
+
+:c,image(Eqs/pair_entropy3.jpg)
+
+where the sum j goes over the neighbors of atom i and N is the number
+ of neighbors. This procedure provides a sharper distinction between
+order and disorder environments. In this case the input parameter 
+ {cutoff2} is the cutoff for the averaging over the neighbors and 
+ must also be specified.
+
+If the {avg yes} option is used, the effective cutoff of the neighbor
+ list should be {cutoff}+{cutoff2} and therefore it might be necessary 
+ to increase the skin of the neighbor list with:
+
+neighbor skin bin :pre
+
+See "neighbor"_neighbor.html for details.
+
+The neighbor list needed to compute this quantity is constructed each
+time the calculation is performed (e.g. each time a snapshot of atoms
+is dumped).  Thus it can be inefficient to compute/dump this quantity
+too frequently or to have multiple compute/dump commands, each with a
+{centro/atom} style.
+
+[Output info:]
+
+By default, this compute calculates the pair entropy value for each
+atom as a per-atom vector, which can be accessed by any command that
+uses per-atom values from a compute as input.  See "Section
+6.15"_Section_howto.html#howto_15 for an overview of LAMMPS output
+options.
+
+The pair entropy values have units of the Boltzmann constant. They are 
+ always negative, and lower values (lower entropy) correspond to more
+ ordered environments.
+
+Here are typical input parameters for fcc aluminum (lattice 
+ constant 4.05 Angstroms),
+
+compute 1 all pentropy/atom 0.25 5.7 avg yes 3.7 :pre
+
+and for bcc sodium (lattice constant 4.23 Angstroms),
+
+compute 1 all pentropy/atom 0.25 7.3 avg yes 5.1 :pre
+
+[Restrictions:] none
+
+[Related commands:]
+
+"compute cna/atom"_compute_cna_atom.html
+"compute centro/atom"_compute_centro_atom.html
+
+[Default:]
+
+The default value for the optional keyword is avg = no.
+
+:line
+
+:link(Piaggi)
+[(Piaggi)] Piaggi and Parrinello, J Chem Phys, 147, 114112 (2017).
diff --git a/src/compute_pair_entropy_atom.cpp b/src/compute_pair_entropy_atom.cpp
index c4e411f761..97461b4eec 100644
--- a/src/compute_pair_entropy_atom.cpp
+++ b/src/compute_pair_entropy_atom.cpp
@@ -47,6 +47,13 @@ ComputePairEntropyAtom::ComputePairEntropyAtom(LAMMPS *lmp, int narg, char **arg
   if (narg < 5 || narg > 8)
     error->all(FLERR,"Illegal compute pentropy/atom command");
 
+  // Arguments are: sigma cutoff avg yes/no cutoff2
+  //   sigma is the gaussian width 
+  //   cutoff is the cutoff for the calculation of g(r) 
+  //   avg is optional and it means averaginf the pair entropy over the neighbors
+  //   the next argument should be yes or no 
+  //   cutoff2 is the cutoff for the averaging
+
   sigma = force->numeric(FLERR,arg[3]);
   cutoff = force->numeric(FLERR,arg[4]);
   if (cutoff < 0.0) error->all(FLERR,"Illegal compute pentropy/atom command; negative cutoff");
@@ -74,7 +81,7 @@ ComputePairEntropyAtom::ComputePairEntropyAtom(LAMMPS *lmp, int narg, char **arg
   nbin = static_cast<int>(cutoff / sigma) + 1;
   nmax = 0;
   maxneigh = 0;
-  deltabin = 2;
+  deltabin = 2; // 2 seems a good compromise between speed and good mollification
   deltar = sigma; 
   peratom_flag = 1;
   size_peratom_cols = 0;
@@ -234,7 +241,7 @@ void ComputePairEntropyAtom::compute_peratom()
           maxbin=bin +  deltabin;
           if (maxbin > (nbin-1)) maxbin=nbin-1;
           for(int k=minbin;k<maxbin+1;k++) {
-            double invNormKernel=invNormConstantBase/rbinsq[bin];
+            double invNormKernel=invNormConstantBase/rbinsq[k];
             double distance = r - rbin[k];
             gofr[k] += invNormKernel*exp(-distance*distance/sigmasq2) ; 
           }
-- 
GitLab