a fast solution method for three-dimensional many-particle

26

Upload: trinhnhan

Post on 05-Feb-2017

234 views

Category:

Documents


2 download

TRANSCRIPT

Page 1: A Fast Solution Method for Three-Dimensional Many-Particle

!#"$%&' (()("+*,(.-/*01234567& 4(8%9-&:)

;=<?>A@CBAD!EA<9FAGH.B?BAH.IJ>LKAMNGOQPSR(TU@WVYXZT[P\F^]Y_`H.DC@a_Jb%KAMNcd@[e^PQBfCFNghR(Hi_JbjH._`DCH._iFAKUklR!HWXnmYMj_J@WVdBAHUF?Kpo<A_`DCH.BqG3MNrUPQBADaH._.FNcd@CsNH._`ItM^uaklBve^H3]#HiPxwJB9F?klBNeLGy<?R%kl_YrNMAz,H.R%klDklB[IJP

|9H.AkaX#~\BNXIJPSIJ<AI`H#@C_0m,@CR(?<AI$kWIJPQ@aBNklOkWBNe tYAAOQPSHWeklIJ>?H.R%klI`P.X|>?HYBAPQuUH._+XZPSIZb3@W|9H.AkaXdklI4t#<NXIJPQBFAt#<NXIJPSB9FA|9H.?kCXW.[F9:r[t

!219(1-

tsN@a<ABNeAkW_Jb HiOQH.R(HiB^I,R(H.I`>A@[e1@C_4XZ@aOQu[PQBAD3I`>A_JH.H&e^PSR(H.BNXPQ@aBNklOOQPQB?HWkl_dH.OkCXIJP&PQIb3A_J@asAOQHiR%XIJ>NkWI=PSB^uU@aOQuUHkOkW_JDCHYB[<AR1sNH._9@W9NkW_JI`POSHWXH.RsNHWeAe^HWePSB(k#s?PQBNe^Hi_=PXPQB[IJ_`@[e^<NHepMh|>AHd?_J@l@XZHWeR(H.IJ>?@[e_JH.OSPQHWX@aB(kWB(PQI`H._$kWIJPSuUHdXZ@aOQ<AI`PQ@aB(XIJ_$kWIJH.Dab1PSB3Vd>AP$>3R%kWIJ_JPSuUHWI`@C_R1<?OQIJPSAOQP.klI`PQ@CB3PXNH._@C_`R(HWeVdPQI`>3IJ>AH,kCXZIR1<AOSIJPS@aOQHhR!H.IJ>A@[eMt4Xhk#_JHXZ<AOSII`>AHR(HiIJ>A@[e0PX.klNkWsAOQH@l4XZ@aOQu[PQBADLA_J@asAOQHiR%X#VdPQI`><ABAT[BA@WVdBNX<NXPQBADL@CBAOSb( R!H.R(@C_`bLklBNe!ZL@aHi_$kIJPQ@aBNXMcnHXZ<AOSI$Xkl_JHdDaPQuUH.B0@C_A_J@asAOQHiR%XVdPSIJ>>^<ABe^_JHWe?X@lNkl_JI`POSHWX9PQBVd>?P$>1 !J.¡a¢$M

£ #) ¤¥!v¦hsN@C<ABeAkl_`bH.OQHiR(H.B[IdR(H.IJ>?@[ep§pkCXZIR1<?OQIJPSN@COQHR(H.I`>A@[ep§9R%klB[bU\NkW_JI`[email protected]§OSPQBAHkl_dH.OkCXZI`PPSIbp§CPSIJH._+klI`PQuUH0X@COQ<?IJPQ@aBLXZIJ_+klI`H.DCb

Mf+IJ@Vd>A@CR¨@a_J_JHXZN@CBNe[H.BNHXZ>?@C<AOesNHkCeAe^_`HWXJXHWepFAHR%kWPQOpD&wJ_&©YIJPiklRMª<AIJH.?kCX&M«HWe[<

Page 2: A Fast Solution Method for Three-Dimensional Many-Particle

~ B3I`>APXNklNH._iFVHdPQB[IJ_`@[e^<NHYk#kCXI,s@a<ABNeAkW_JbH.OQHiR(H.B[IpR(H.I`>A@[e(\jdg, N@C_XZ@aOQu[PQBAD4IJ>A_`H.He^PQR(HiBNXZPS@CBNkWO=OSPQBAHkl_0H.OkaXZI`PPSIZbLA_J@asAOQHiR%XIJ>NkWI3PQB[uU@COQuUHk Okl_`DCH(B[<ARsHi_@l4Nkl_`IJPOQHWXH.RsHeAe^HWe3PQB!ksAPQBNe[H._.M HY_JHHi_hIJ@0IJ>?@[email protected] XkaX=R%kWB^bU\kl_JI`POSHA_`@CsAOSH.R%XklBNeIJ@0IJ>AHBAH.V R(H.IJ>?@[ekaX#Eg= rNEkaXZINPSBAHWkl_4ghOkaXZIJP0klB[bNkl_`IJPOQH?_J@Cs?OQH.RrU@COSuUH._.MEA@a_ R%kWB^bU\kl_JI`POSH A_`@CsAOSH.R%X&FEg, r @aR1s?PQBAHWXIJ>AH_`@CsA<XZIJB?HWXJX%@l!@aB^uUH.B[IJPS@CBNkWO

jdg=X3VdPQI`>6X<ANH._JPS@C_Hi_Z@C_`R%klBHUM6|9@e?klIJHF#Eg= r >NkCXsHiH.B6kWAAOQPSHWeIJ@ kWBNklOSbAXPX@lA_`@CsAOSH.R%XI`>NklIdPSB^uU@aOQuUHI`H.BNX,@lNkW_JIJPOQHX @aB%e^HWXT^I`@C%&@CR(A<?IJH._+XhklBNe(>[<ABNe^_`HWeAXh@lkl_ZIJP&OQHWX@aBqXZR%kWOQO9NkW_$klOSOQH.O@aR(A<AI`H._$X&Myr@COQ<?IJPQ@aB@lA_`@CsAOSH.R%XVdPQI`>IJ>A@a<NXJkWBNeAX@lNkW_JI`POSHWXPXHWkCXPQsAOSH#VdPSIJ>OkW_JDCHYNkW_$klOSOQH.O@aR(A<AI`H._$X&M|>AH0kls?PQOQPSIbIJ@X@COQuUH4XZ< >A_J@asAOQHiR%XR%k.b3e[_$kR%klI`P.kWOQOQb!$>NklBADaH3IJ>AHH.Oe%@ldR!P_`@CR(Hi >klBAPX&M#~\BI`>APXd_JH.Dkl_+epF9VHsNH.OSPQH.uUH0I`>NklI0PQBLIJ>AHBAHWkl_,<AIJ<?_JHuakl_`PQ@a<NXPX`XZ<AHWX@CBH._JB?HWe1VdPQI`>(IJ>?H#@WuUH._+klOQON_`HWXZN@CBXZH0@l @CR!@XZPSIJHYR%klI`H._JPklOXkaXZ>APSB9F, C"!4VdPQOQOsH!XZHiIJIJOSHWeklBe IJ>AH@[&<NX:@Wd_JHXZHWkW_$$>VdPQOQOsH!XZ>?P IJHe IJ@%A_J@asAOQHiR%X@CBN&H._JBAHe VdPSIJ>%HWe^DaH0H$#HWI$X&F^PSB^I`H._JOklR(PSBNkl_XZI`_JHWX`X#e[PXZI`_JPSsA<AIJPS@CB9F?XZP&%.HH$#HWI+XFAH.I+CM~\BL@C_$e[H._IJ@ @CR(kl_JHEg, r%VdPSIJ>q@aB^uUH.B[IJPS@CBNkWOjdg,XFNOQH.IY<NXd_JHW.kWOQOIJ>NkWI0@aB^uUH.B[

IJPQ@aBNklOjdg,X,DCPQuUH4_JPXHIJ@e[H.BNXZHkWOQDCHisA_$kWPA_J@asAOQHiR%X,klBNepFAPSB%IJ>AH0jdg,8@aR(R<ABAPQIbaFCPQIPX<NXZI`@CR%kW_Jb0IJ@X@COQuUHIJ>A@XZHdA_`@CsAOSH.R%X9<NXZPSBADe^PQ_`HWIR(H.I`>A@[eAXXZ<N$>kCX'(6e[HW@CR!@XZPSIJPQ@aB3@C_@CIJ>?H._$Xde^Hi_JPQuUHe_J@aR ]4kl<NX`XZPkWB H.OSPQR(PSBNklI`PQ@aB9M~ B%I`>APXd.kaXZHUF?k?_J@Cs?OQH.R/VdPQI`>'<?BAT^B?@lVdBX_JH*)U<APQ_`HWX(,+YR(H.R@C_`b%klBNe!Z-. @CNH._$kWIJPS@CBNX&M1|>?HWXZH3_`H*)<?PQ_JHiR(H.B[I$XkW_JH3A_+kCI`P.kWOQOSb<ABNkCiH.AI+klsAOSHy@a_dR%klB[bU\NkW_JI`POSH#A_`@CsAOSH.R%Xh@C_YVd>AP$> /. .¡*0WM~\BI`>AH(OkCXZI e^HW.kae^HUFhPSIJH._+klIJPSuUHXZ@aOQ<AI`PQ@CB XZIJ_+klI`H.DCPSHWX>Nk.uUH1sNHW&@CR(HPQBN&_JHWkaXZPQB?DCOQb_`HW@aDl

BAP1%.He s^b IJ>AH%jdg, @CR!R1<AB?PQI\bUML~\B Nkl_`IJP&<AOkW_.FhH$#HWI`PQuUH(PSIJH._+klI`PQuUH!XZ@COS<AIJPS@CB XZI`_$kWIJH.DaPQHWX@C_R%kWB^bU\kl_JI`POSHA_J@asAOQHiR%X(>Nk.uUHsNH.H.Be^HiuUH.OQ@aHe s[b32h>NkWB^ |PQH.B klBe GPSR ` 4*5^+MEg, r klOX@ _JH.OSPQHWX3@CB kWB PSIJH._+klIJPSuUHLXZ@COS<AIJPS@CB XIJ_$kWIJH.Dab6IJ>NkWI.F@C_%R%kWB^bUNkl_`IJP&OQH(A_`@Cs^OQH.R%X&F_`H*)<?PQ_JHX:@CB?OQb kWB !ZL0R(H.R(@a_JbqkWBNe ( 0@CNH._+klI`PQ@CBX§#PSB XZ<?R(R%kl_`bUFhEg= rPXkWB!ZL#R!H.IJ>A@[eM|>APX0e^_$kWR%klI`P:PSR(A_`@luUHiR1H.BUIdPSBNH._Z@a_JR%kWBNHUFpPSBq@CR!Nkl_`PXZ@aBIJ@

Page 3: A Fast Solution Method for Three-Dimensional Many-Particle

@CB[uUH.B[IJPQ@aBNklO?jdg,XFU&@CR(HWXp_J@aR IV@0X@C<A_+HWXMEPQ_+XZIiF^kaX,k#_J<AOSHUFCR%kWB^bUNkl_`IJP&OQH,A_J@asAOQHiR%XDCPQuUH,_`PXZH,I`@VH.OSO @CBe^PQI`PQ@CB?HWeklOSDCH.s?_$klP=A_`@CsAOSH.R%XA@C_Vd>AP$>H$#HW&IJPQuUHA_`HW@CBe^PQI`PQ@CB?H._$X.klBsH4@aBNXZI`_J<NI`HWepMtXhk_JHWX<AOQI,I`>AHdB^<?R1sNH._@l_JH*)U<APQ_`HWe:PSIJH._+klI`PQ@CBX PX<NX<NklOSOQb3XR%klOSOAkWBNe3u^PS_ZIJ<NkWOQOQb3PSBNe^H.NH.Be^H.B[I=@W klBNe(I`>AH._JH&@C_JHYIJ>?H@CNH._+klIJPS@CBL@a<AB^IPXh_JHe^<NHWe3_J@CR!Z-.IJ@(,+$M,rUHW@aBNepFUkWIhHWka >PQIJHi_$klI`PQ@aB9FlR klIJ_`PQU\uUHIJ@a_pR<AOQI`PQAOSP.kWIJPQ@aBPX9NH._Z@a_JR(Hey<XZPQB?DIJ>AHkCXZIR1<AOSIJPS@aOQHhR!H.IJ>A@[e\cd@CT[>AOQPSB9F[ C §A]Y_JH.HiBADUkW_$e3klBe3cn@aT^>AOSPQB9F? CUC9Vd>?P$>3@aR(A<AI`HWXIJ>AHA_J@[e^<I<NXPQBAD@CBAOSb( dR(HiR(@C_`bklBNe !ZL@aHi_$klI`PQ@aBNXM0|>AH_JHWkae^H._Y<AB^kWR(PQOSPkl_VdPQIJ>!IJ>AH#kaXZIR1<AOSIJPS@aOQHdR(HiIJ>A@[e%\Eh PXh_JH&H._J_`HWe1IJ@ cn@aT^>?OQPQBLJC $F9]Y_JH.HiBADUkW_$e(kWBNecn@aT^>AOSPQBLJCa$Fp]Y_JHiH.BADUkW_$e ` aC+FNklBNe%]#_`H.H.BADkl_+e(klBNe!cn@aT^>AOSPQBLJ4U$M#H.uUH._`IJ>AHiOQHWX`XFPQB @C_+e^H._3IJ@q<ABNe[H._$XI$klBeIJ>?PXNkWHi_.FI`>AH%_JHWkae^H._3BAH.HeAXIJ@qT^B?@lV/uUH._`b OSPQIJI`OQH%kWsN@C<AIIJ>AHEhF[klBNe1I`>APX,PQB^@a_JR%kWIJPS@CBPX=?_J@Wu^Pe^HWe3PQB rHWIJPS@CB!^M,EA@C_dB?@lV0F[PQI,PXXZ< PSH.B^I,I`@T[BA@WVIJ>NkWI#I`>AHEh PXdk3R(H.I`>A@[e@C_4XZ@aOQu[PQBAD3I`>AH0R%klB[bU\sN@[e^b(HiOQHWI`_J@XZI$kWIJP.X=A_`@CsAOSH.RMtYOQIJ>?@C<ADa>(IJ>AHYPe[HWk@lkWAAOSb^PSBAD0IJ>AH0E I`@PSIJH._+klIJPSuUHdX@COQ<?IJPQ@aB(@lPQB[IJH.Da_$kWOH )<NkWIJPS@CBNX

eAklI`HWXsNkC$T I`@(IJ>AH!XZHiR(PQBNkWO9NklNH._0@WYcd@CT[>AOQPSBJC $FIJ>AHi_JHkl_`H3DC@U@[e_JHkCXZ@aBNXVd>^bLIJ>AHEh >NkaX:[email protected]_JH*)U<AH.B[IJOSbkl?AOQPSHWe'I`@ I`>A_JH.H&e^PSR([email protected]`PsN@C<?BNeAkl_`buaklOS<[email protected] XFNHWXZNHWPklOSOQb(I`>A@UXHkl_JPXZPSBAD%PQB OSPQBAHkl_YH.OkaXZIJPPQIb%[email protected][BA@COeAX B[<AR1sNH._+X>^bUe^_J@[e^b[BNkWR(P.X&MEA@a_#I`>A_JH.H&e^PSR(H.BNXPQ@aBNklOs@a<ABNeAkW_JbU\uakWOQ<AH3A_`@CsAOSH.R%X&F^IJ>?H3Eh8PX<NXHWePSB^_JH*)U<AH.B[IJOSbsNHW.kl<XZHdPQIPXR(HWklB?PQBADW<AOU@CBAOSb0P IJ>AHA_J@asAOQHiR7XP1%.HH.?H.HWeAXh !J.¡*-%.¡*0.$M~ Bkl_JI`P<?Okl_iF9VHkl_`H(k.V kl_JH3@Wd@aBAOQb@aBAH1_`HWXZHkl_$$>DC_`@C<AqIJ>klI0<NXZHXyI`>AH!EhIJ@LXZ@aOQuUHs@a<ABNeAkW_JbPSB^I`H.DC_+klO,H*)U<NklI`PQ@aBNX0@C_`_JHWX@aBNe^PSBADIJ@LIJ>A_`H.H e^PQR(HiBNXZPS@CBNkWOkWAOkaH X0H*)U<NklI`PQ@CBkWs@a_$X F4 5[$M~\B &@CB[IJ_$kaXZIiFIJ>?H(sA_JHklTU\H.uUHiB N@CPSB^I0PSBIV@Le^PQR!H.BNXZPS@CBNXPX (`.¡ + $FaklBepFC@C_,<?_J_JHiB^Ie^HWXT^I`@C@aR(A<AI`H._$X&F.IJ>APXhe^P #H._`H.BNHPXXZPSDCBAP .klB[I.MEA@C_JI`<ABNkWIJH.OSbUFIJ>APXdXZPQI`<NklI`PQ@aB [email protected]<NXH0_JHWHiB^I`OQb%]#_`H.H.BADkl_+ekWBNe cd@CT[>AOSPQBqJ4UdPQR!A_J@WuUHWeIJ>AHLIJ>A_`H.H e^PQR(HiBNXZPS@CBNkWOdEh klBNe _`HWe^<N&HWePSI$X3sA_JHWkWTH.uUH.B N@CPSB^IIJ@ (`.¡4++§#@W@C<A_+XZH0<A_JI`>AH._dA_`@CDC_`HWXJXYPQBLe^HWXT^I`@C&@CR(A<?IJH._+XVdPQOQO9R klTUHI`>AH0IJ>A_`H.H e^PQR([email protected](@C_`HkWIJIJ_+kCI`PQuUHUM|>AH4_JHWkaXZ@CB%Vd>[b1I`>AH4Eh PXhBA@aI_JH*)U<AH.B[IJOSb3kl?AOQPSHWeI`@XZ@aOQ<AI`PQ@CB@lsN@C<ABNe?kl_JbUuCkWOQ<AH3?_J@Cs?OQH.R%X@aIJ>AHi_IJ>NkWBIJ>A@XZH3Da@luUHi_JBAHWe(s[b klAOkCH X#H*)U<NklI`PQ@aB PXYe^<AH

!

Page 4: A Fast Solution Method for Three-Dimensional Many-Particle

IJ@I`>AHdkCI,IJ>klI kWAOkaH X H )<NkWIJPS@CB(PX=e[PQ_JHIJOSb:_`H.OkWIJHWe3IJ@IJ>AHYR%kWB^bU\sN@[e^b3HiOQHWI`_J@[email protected]>APSOQH(@aIJ>AH._Nkl_`IJPklOYe^P #pHi_JH.B[IJPklO,H*)U<NklI`PQ@CBXkl_`H%BA@CIiM YH.uUH._JI`>AH.OSHWXJX0XZHiuUH._$kWO=_`HXZHWkW_$$>AH._$X>kuUH4sHiH.BLklsAOSHIJ@3@WuUH._+@CR(HYIJ>?PX@CsXZI$kaOQHM,t#R(@CB?DIJ>?@UXZH0VH&PQIJHYIJ>?HXIJ<Ne[PQHWX@l]Y_JH.H.B[sNkW<ARH.IhkWO\MJ "CU,@aBNH._`BAHWe:VdPSIJ>IJ>AHIZV@l e^PQR!H.BNXZPS@CBNkWOUsAPQ>kl_JR!@CBAP=H*)U<NklI`PQ@aB9F@l 2H.PS_$H klBNe kl?PQH._J4 ,@aBNH._`BAHWeVdPSIJ> k.u^PSH._ XpH*)<klIJPS@CBNXh@WpOSPQBAHkl_hH.OkCXIJP&PQIbUFklBNe @WrklBADklBAP klBNe @` "&@CBNHi_JBAHWeVdPQI`>IJ>?_JH.H e^PSR(H.BNXPQ@CBklOhH*)U<NkWIJPQ@aBNX:@[email protected][BA@aOeAX\\B[<ARsHi_$X>^b?e^_`@[e^b^BklR(P.X1rUIJ@aTUHWX H*)<klIJPS@CBNX$$M |>A@UXH _JHXZHWkW_$$>AH._$X0R(@[e^P NHWeH.PQI`>AH._I`>AH0PQB[IJH.Da_$kWOpH*)U<NklI`PQ@aBNXIJ@!kC.&@CR(R(@[e?klIJHYIJ>?HEZ]Y_JH.HiB^sNkW<AR H.IYkWO\MSFp "CUd@C_IJ>AHEh8IJ@!kC.&@CR(R(@[e?klIJH4IJ>AHPQB[IJH.Da_$klOH*)U<NklI`PQ@CBX: 2H.PQ_+H0kWBNe kl?PQH._iFp " §r[klB?DUklB?PklBNe @AFY 4 +M rU< > R(@[e^P .kWIJPQ@aBNX3kl_`H%BA@CB^IJ_`PQu[PkWOdklBe@lI`H.B A_J@asAOQHiR XZNHWP CM ~ BNkl_`IJP&<AOkW_.FhIJ>?H kWAA_`@UkC$> @W]Y_`H.H.B[sNkl<AR ` "C.klBAB?@CI3sNH%H.^I`H.BNe^He IJ@qIJ>A_`H.He^PQR(HiBNXZPS@CBNX&FIJ>?HuUHi_$XZPS@CB@WI`>AHLEh e^H.uUHiOQ@CNHWeqs^b r[klB?DUklB?PkWBNe @ J "U3.kWBABA@aIsH0HWkaXZPSOQbH.^I`H.BNe^He1IJ@IJ>A_`H.H e^PQR!H.BNXZPS@CBNkWOAOQPSBAHWkW_dH.OkaXZIJPPQIb3sNHW.kW<NXZH0PSIPXOSPQR(PSIJHWeI`@3_JPSDCPeNkl_`IJP&OQHWX&F#kWBNe IJ>?HuUH._$XPQ@CB6@WIJ>AHqEh e^HiuUH.OQ@aHes[b 2H.PS_$HLkWBNe 4klAPSH._J4 PXXZ@CR!H.Vd>NklIdPSBAH PQHiB^IiM|>?H3A_JPSBNPQklO_JHWX<AOQI+X#_JH.N@a_JIJHe PQBIJ>?PX#NkWNH._0kl_`HkCX#@COSOQ@WVYXM0EPQ_$XI.FpVHHWXI$kWsAOQPXZ>AHWe

k1_JH.OklI`PQ@CBXZ>APS sNH.IVH.HiBsN@a<ABNeAkW_JbPQB[IJHiDC_$kWOH )<NkWIJPS@CBNXd@lIJ>?_JH.H e^PSR(H.BNXPQ@CBklOpOSPQBAHkl_YH.OkaX IJP&PQIb(kWBNe IJ>?H3EhM|>APXd_JHiOklI`PQ@aBNXZ>APSkWOQOQ@WVYXd@CBAH3I`@(&@CRsAPQBAHyklB[bjdg,@a_JR<AOkWIJPS@CBVdPQIJ> klB[bEh@C_`R1<AOklI`PQ@CBM rUHW@aBNepFskCXZHe @CB IJ>?H(BAH.V _`H.OkWIJPQ@aBNXZ>?PQ9FVH!e^H.uUH.OS@CNHWeEg, rkykl_$kWOQOSH.OpmCm @[e^HYIJ>klIdiklBsHYHWkaXZPQOSb3R(@[e^P NHWeIJ@ka.@CR!R(@[eAklI`HsNH.I`IJH._jdg=8@C_YE@a_JR<AOkWIJPQ@aBNXM|>?H_`H.R%kWPQBNe^Hi_:@WI`>APXklNH._PX:@a_JDklBAP&%.HWe kCX@COSOQ@WVYXM ~ B6rUHIJPS@CB [F=VH e[HWXJ_`PQsNH

>A@WV7IJ@Lkae^@CAIIJ>AHE@a_XZ@COS<AIJPS@CBq@l PQB[IJH.Da_$klOH*)<klIJPS@CBNX&MEPQ_+XZI.FN@C_0IJ>A@XZH(<?B^klR(PSO Pkl_VdPQI`>IJ>AH EhF9VH!XZ<AR(R kl_JP&%.HPSI$X0sNkCXP1Pe^HWkCX&M|>AH.B VH(H.[AOklPQBq>A@WV IJ@Hi^I`H.BNeIJ>[email protected]?H._$kWOQP1%iHWe(R%klB[bsN@[e^b(H.OSHWI`_J@[email protected] XklBee^H.R(@aBNXZI`_$klI`HIJ>klI0k@CRsAPQBklIJPS@CB@WYIJ>A@XZH(A_`@CsAOSH.R%XiklBsNH(<NXHWe I`@ _JHiA_JHWXH.B[IIJ>AH!sAPQ>NkW_JR(@aBAPTUHi_JBAH.OM~ BrUHWI`PQ@CB!^FAVHY@C_JR<AOklIJH4jdg, H*)U<NklI`PQ@CBXM, @UXId@lIJ>AH4R%klI`H._JPklONA_JHXZH.B[IJHWe3PQBrUHIJPS@CB,!

5

Page 5: A Fast Solution Method for Three-Dimensional Many-Particle

PX,XZI+klBNeAkW_$epF[H.?H.AI@a_,BAH.V Hi^A_`HWXJXPQ@aBNX@C_,IJ>AH<?BNeAklR!H.B^I+klONX@COS<AIJPS@CBNXhIJ>klIdkWOQOS@lV <NX,IJ@_JH.OklIJHYe[PXJ&_JH.I`P1%.HWesN@C<?BNeAkl_`b1PSB^I`H.DC_+klO[H*)U<NklI`PQ@CBX @WIJ>A_`H.H e^PQR(HiBNXZPS@CBNkWO^OSPQBAHWkW_HiOkCXIJPPQIb0IJ@k(@CRsAPSBNklI`PQ@CBL@l XZPS^I`H.H.B%DCHiBAH._$kWOQP&%.HWeR klB[bsN@[e^b H.OSHWI`_J@[email protected] XM~\BrUHIJPS@CB5AFAVHye^HWXJ&_JPQsNHI`>AHPSR(AOQHiR1H.B[I$kWIJP@CB@lI`>AH0PQI`H._$kWIJPQuUHYX@COS<AIJPS@CBLXZI`_$klI`H.DCbM~\B rUHIJPS@CB F^VHH.^?OklPSBI`>AHdsNkCXPPe^HWkCXsHi>APQBNe0IJ>AHdkl_$kWOQOSH.O^uUHi_$XZPS@CB3@W9Eg= rMU~\B1rUHW&IJPQ@aB ^FCVHdA_JHXZH.B[IIJHWXIA_`@CsAOSH.R%XhPQB(Vd>?P$>1I`>AHB[<AR1sNH._,@WNkl_`IJP&OQHWX,uCkW_JPSHWX _`@CR IJ@ !*5 ![M~\B rUHIJPS@CB F^VHXZ<AR(R kl_JP&%.Hd_JHXZ<AOSI$XdklBesA_`PQH b!e^PXJ<NX`X <?IJ<A_`H:V@C_`TpM

$

~ B0I`>AH,R%klB[bsN@[e^b0H.OQHIJ_`@UXZI+klI`[email protected].@CBAH,PX9DCPSuUH.B0kd<ABAPSI<AsNHd@aB^I+klPSBAPQBAD$>Nkl_`DCHWX

klIYN@UXZPSIJPS@CBNX kWBNe kaXZTUHWe(IJ@!e^H.I`H._JR!PQBAHdI`>AH@C_$&HWX kCI`PQBAD@CB%I`>A@UXH3$>Nkl_`DCHWX&M|>AH@C_$&HWXdkl_`HiklO&<AOkWIJHWe(PQB%I`H._JR%Xh@W,IJ>AH0N@CI`H.B[IJPkWOXdkaX

`l

klBNe

= !"$#&%$' ")(#

"* + " *

§ ,.-/0/0/1-23/

6>APQOSH:@a_$e^PSBNkl_`PQOSbI`>AHXZ@aOQ<AI`PQ@CB@W,IJ>[email protected]_`H*)<?PQ_JHX#( +@CNH._$kWIJPS@CBNX&F@C_YOkW_JDCH4F@CBAH .klBq@aR(A<AI`HIJ>?H:@C_+HWXdVdPQI`> IJ>?H3EhPSB !9@aHi_$klI`PQ@aBNX#@CB?OQbUMqH@C<?IJOQPSBAH0IJ>AHEh PSB I`>AH1B?H.^I0NkW_$kWDC_$kWA>9MEA@a_3e^H.I+klPSOXF9I`>AH1_`HWkCe[H._0PX:_`HH._`_JHWe I`@cd@CT[>AOQPSB ` C +F]#_`H.H.BADkl_+eklBecd@CT[>AOQPSB3JCa$FA]Y_JH.HiBADUkW_$e3JCCU$F[klBe]Y_JH.HiBADUkW_$e0klBNecn@aT^>?OQPQB` 4 +M|>?HE H.[AOS@CPQI+XpIJ>AH,VH.OQO\T[BA@WVdBPe[HWkIJ>NkWII`>AH,N@CIJHiB^I`PklO4H.OePQBe^<NHWeYs^b0kYOS<NXZI`H._

@l$>Nkl_`DCHWXY.kWB sNHklA?_J@W^PSR%klI`HWe1s[b%IJ>AH0N@CI`H.B^I`PkWO NHiOe(@lIJ>?H.PQ_dI`@CI$kWOh$>NkW_JDCH@C_Y>APQDa>AH._@C_$e[H._,R(@CR(HiB^I+XFakCX=OS@CBADkaX=IJ>?H#I+kl_JDaH.Id$>NkW_JDCHX=kW_JH4XZ< PSH.B^I`OQbkl_ k.Vdk.b_J@aR IJ>AH4XZ@a<A_$H >NkW_JDaHWXM,~\B3IJ>?HEhFWIJ>?HOS<NXZI`H._$Xkl_JHd@C_`R(HWeu[Pk#>?PQH._+kl_$$>AP.klOAkl_JI`PQI`PQ@CB?PQBAD@WIJ>AHd<?BAPQI<AsNH0PQB[IJ@3*-CF 54-WFN*-WF^H.I+CMhPe^HiB^I`P.kWOp&<AsAP4H.OQOX=Vd>AP >!PQB[IJH._+kCI klR(@aBADI`>AH.R%XH.OQuUHWXhsNkaXZHWe@CBI`>AH1R kC$>APQBAHi_Jb%@l IJ>AHR1<?OQIJPSN@COQH klBNe XZ@COSPeqXZA>?H._JP.klO,>NkW_JR(@aBAP(OQ@[.kWO YH.[NklBNXPQ@aBNXM

Page 6: A Fast Solution Method for Three-Dimensional Many-Particle

|>AH@CR!A<AI$kWIJPS@CB ?_J@[H.HeAXPQB I`>A_JHiHXZI$kWDCHWX&MEPQ_$XI.FpI`>AHR<AOSIJPQN@aOQH@UH PQHiB^I+X: >NkW_JDaHUFe^[email protected]$CM0kW_JHd&@CR(A<?IJHWe#@C_ klOQO?H.OSOX&M rHW@CBepFlI`>AHYH.OSO H.OSOlPSB^I`H._$kaIJPS@CBNXkl_`Hd@aR(A<AI`HWeXZ@I`>NklIIJ>AHOQ@[.kWO^Hi^NkWBNXZPS@CB @UH PQH.B[I$Xp@WNIJ>?HdXZR klOQOSHWXZIH.OSOX9kl_`H_`H.OkWIJHWeIJ@4IJ>AH,R<AOQI`PQN@COSH@UH PQH.B[I$X&M|>?PQ_$eFI`>AHN@CI`H.B[IJPkWOXdkW_JH@aR(A<AI`HWe1s[b kaeAe^PSBAD1I`>[email protected][IJPklOXYe^<AHIJ@IJ>AHOQ@[.kWO9H.^klBNXPQ@CBX#kWBNeBAHkl_Js[b%$>Nkl_`DCHWX&M,tXdk3_JHWX<AOQIIJ>AHXZ@aOQ<AI`PQ@aB PXHi^A_`HWXJXHWe%kCX

, %! #&%

A !" (#

"* + " *

H._`_J@a_ !

klBNe

% ! #&%

_JH.R%kWPQBAPSBADqIJH._`R%X / 5^

~ B%I`>AHWXZHH*)U<NklI`PQ@CBXF Xdkl_`HI`>AHOS@[.klO9H.[NkWBNXZPS@CB&@UH %&PQH.B[I$X&F XYkl_`HIJ>?HX@COSPe XA>AH._`P.kWO>Nkl_`R(@CBAP.XF3PXI`>AH&H.B^I`H._h@W9IJ>?HXR%klOSOQHWXIhHiOQO?@CB[I$kWPQBAPSBAD F[klBeI`>AH#I`PQOe^HYXZPSDCB!e^H.BA@aIJHWXIJ>NkWII`>AH XZHW@aBNe XZ<AR PSBNOQ<e^HWX3BAHWkl_`s^b $>NkW_JDCHX @aBAOQbM |>AHq<AANH._%OQPSR(PQI3@WIJ>AHOS@[.klOH.^klBNXPQ@CBqPX 6W + sNHW.kW<NXZHIJ>APX#PXYIJ>AH3B[<ARsHi_d@ldI`>AHXZ@aOQPeXZA>AHi_JPiklO>Nkl_`R(@CBAP.X#@W@C_$e[H._OSHWXJXh@C_dH )<NkWONIJ@M,|>AHH._`_J@a_IJHi_JRFCVd>?P$>(.kWB(sNHsN@C<ABe^HWe QFCPX@CB[IJ_`@COQOSHWes^bL$>A@U@UXPQBAD%k!XZ< PQHiB^I`OQbOkl_`DCH M #@aIJHI`>NklI4R(@UXI#@WI`>AH3E R%ka >?PQBAH._`b1PX#BAH.He^HWe@C_@aR(A<AI`PQBAD X qH.uaklOS<NklI`PQ@CB @W!kWBNe6 5[#@C_3DaPQuUH.B X:PXkXPQR(AOSH@CR!A<AI$kWIJPS@CBHi_Z@C_`R(HWe1N@aPQB[IZVdPXZHUF[@CBAH >kl_JDaHklIYk3IJPQR!HUM~\BL@C_$e[H._#I`@1Hi^I`H.BNe(IJ>AHEh8IJ@R%kWB^bU\sN@[e^b%A_`@CsAOSH.R%XkCXJX@[PkWIJHWe VdPSIJ>ve^PX`_JH.I`P1%.He

PQB[IJH.Da_$klO9H )<NkWIJPS@CBNXd@l,I`>A_JHiHe[PQR(H.BXZPQ@aBNklOOQPSBAHWkl_dHiOkCXIJPPQIbVH0PQB[IJ_`@[e^<[email protected]?H._$kWOQP1%kIJPQ@aBNX#@W,IJ>AH3R klB[bsN@[e^b%H.OQHW&IJ_J@XZI+klIJP.XA_`@CsAOSH.RMY|>AH _$XZI4DCH.BAHi_$klOSP1%WkWIJPS@CB PXn_`H.OkWIJHWe%IJ@IJ>AH1kaI0IJ>NkWI@COSOQ@[.kWIJPS@CB jdg,X0kl_`H(H.BNe^@WVHWe'VdPSIJ>IV@ N@CPSB^I XZH.I+XM"!#B?H(N@CPQB[IXZH.IPX@C_JR!HWe1s[bI`>AH0BA@[eAkWOhN@CPQB[I$XdVd>AP >Lkl_`H3kWOXZ@<NXHWe kaXnI`>AH@COSOQ@[.kWIJPQ@aB N@CPSB^I+XMd|>?H@CI`>AH._@aPQB[IXH.IYPXn@a_JR(HWes^b%IJ>?HPQB[IJHiDC_$kWIJPS@CB N@aPQB[I$XM4|@!e^PXZIJPSBADC<?PXZ>sNH.IVH.H.B%IJ>AHXZHIV@(XH.I$XVH(_`H.DUkW_$e I`>AH%@aOQOQ@[iklIJPS@CBN@aPQB[I$XkaXIJ>AH!I$kl_`DCH.IN@CPQB[I$X3kWBNe I`>AH(PQB[IJHiDC_$kWIJPS@CBN@CPSB^I+XkaXIJ>AH XZ@C<?_$H@aPQB[I$X&M#~\B%Vd>NkWI#@COSOQ@WVYXFAVH3e^H.B?@CIJHIJ>AHI$kW_JDCHiI#N@CPSB^I+Xns[b 6klBNe%IJ>AH XZ@a<A_$H@aPQB[I$Xs[b$#=M!|>AH(XHW@aBNe DCHiBAH._$kWOQP&%WklI`PQ@CBPX_JH.OklIJHe'I`@I`>AH3kCI0IJ>klII`>AH<ABNe?klR(H.B[I$kWOXZ@COS<AIJPS@CBNXd@l OQPQB?HWkl_dH.OkCXIJP&PQIb1PSB^uU@aOQuUHYIJ>AHsAPQ>kl_JR!@CBAP * %# * TUH._JB?H.OpI`@1Vd>AP >FAVdPQIJ>?@C<AI

Page 7: A Fast Solution Method for Three-Dimensional Many-Particle

R(@[e^P .klI`PQ@CBXFI`>AHEh .kWBABA@CI!sHvklA?OQPQHepM ~ B I`>APX1_`H.DUkW_$epFVH @asNXZHi_JuUH I`>NklI@CB?OQbR(PQBA@a_dR(@[e^P .klI`PQ@aBNX#kW_JH_JH*)U<APQ_`HWe(PQBL@C_$e[H._IJ@!H.^I`H.BNe(IJ>?H3EhI`@(A_`@CsAOSH.R%XPQBLVd>AP$> P IJ>AH! >kl_JDaHWXkWI:I`>AH!XZ@C<?_$H3N@CPSB^I+X0kl_JHA_`HWXJ_`PQsNHWe PSBIJH._`R%X#@l k(<ABNI`PQ@CB # klBNePQP ,I`>AH@C_+HWXYklIYIJ>AHI$kl_`DCH.IYN@CPSB^I+Xkl_`HA_`HWXJ_`PQsNHWe(PQBLIJH._`R%X=@Wk1OSPQBAHWkW_dR%kl@CB

XZ@IJ>NkWI

2-2 p % ! #&%

N2-2 _`H.R%kWPQBAPSBAD I`H._JR X /

~ I%PX @as^u[PQ@a<NX(IJ>NkWI PX%kN@CPQB[IVdPXH@aR(A<AI+klI`PQ@CB IJ>klIL.klB sNHPQR!AOQH.RH.BUI`HWe@CBNHIJ>AHLEh PXXZ<?AAOQHiR(H.B[IJHWeVdPQI`> k «A_`HA_J@[HXJXZPSBADh?_J@[HWe[<A_JH(@C_3HiuCkWOQ<NkWIJPS@CB @l # klBe(k ª@XZI\A_`@[HWX`XZPQB?DAA_`@[HWe^<A_`H#@C_,H.uakWOQ<NkWIJPQ@aB(@l -2 \M #@CI`H#I`>NklIVdPQIJ>IJ>AHPQB[IJ_J@[e[<NIJPS@CB%@lI`>AH0DCH.BAHi_$klOSP1%.He R klB[bsN@[e^b%H.OQHW&IJ_J@XZI+klIJP.XA_`@CsAOSH.R VH klR(?OQP b!klB@CsNXH._JuaklI`PQ@aB(@l,]Y_JHiH.BADUkW_$eJ4*5^,IJ>klII`>AH4Eh¥PXkXH.NkW_$klI`PQ@aB(@luakl_`PkWsAOQHWX,IJHW$>ABAP )<?H@C_I`>AH>NkW_JR(@aBAPdTUH._`BAH.O * "# * % M,E_`@CR I`>APX,Hi_$XZNHW&IJPQuUHFVH#Hi^I`H.BNe3IJ>AH0E¥IJ@IJ>AHTUH._JBAHiO1 #h * %# * % Mt4XdklB%H.?klR(AOSHUFCOSH.I<NXY@aBNXZPe^H._>?@lV IJ>?H:DaH.BAH._+klOSP1%WkWIJPQ@aB(VH4>Nk.uUHPQB[IJ_`@[e^<NHWe%.kWB%sH

klAAOSPQHWe%IJ@(I`>AHs?PQ>NkW_JR(@aBAPTUH._`BAH.O\M ~ BLIJ>AH0@COSOQ@WVdPQBADH.[A_JHXJXZPS@CB9FNVH<NXZH3I`>AH3X<AsNX`_JPSAI$XIJ@!e^H.BA@aIJH3mkl_`IJHWXPklBv@[email protected][I$X,@l klBNe #

* %# *

* %# * #

%

* %# * % #

+

* %# * + #

-

* %# * - #

* %# * # # #

/

~ BI`>APX_`H.A_JHXZH.B[I$kWIJPQ@aBIJ>AH3>kl_JR!@CBAPTUH._JBAHiO9PX0e^PXIJPSOQOQHe _J@aR I`>AH1s?PQ>NkW_JR(@aBAP0@CBAHkWBNe

Page 8: A Fast Solution Method for Three-Dimensional Many-Particle

kCX,k:_`HWXZ<AOSIhIJ>?HdR%klB[bsN@[e^b3A_`@CsAOSH.R kaXJXZ@[&PklI`HWe1VdPSIJ>3IJ>?H#I`>A_JH.H&e^PSR(H.BNXPQ@aBNklO[sAPQ>kl_JR!@CBAPH*)<klIJPS@CB%PX_`HWe^<N&HWe(IJ@NuUH0DCHiBAH._$kWOQP&%.HWe1R%kWB^bUs@[e[b(H.OQHIJ_`@UXZI+klI`P.X,[email protected] XM~\Bv@CBN&OQ<NXPQ@CBL@lIJ>?PXYXZHWI`PQ@aB9FVH@asNXZH._`uUHI`>NklIiF9XZI`_JP&IJOQb XZNHWklT[PQB?DAFpI`>AH3Eh8PXYBA@CI

klB(!Z R(H.IJ>?@[e:sNHW.kW<NXZHPQB3@C_+e^H._IJ@@aAIJPSR(P1%iHhIJ>AHdNH._@C_`R%klBN&H=I`>AHkl_JI`PQI`PQ@CB?PQBAD0e[H.AIJ>XZ>A@a<AOe sNH!kOS@CDUkW_JPSIJ>AR(Py<ABIJPS@CBq@l0M H!k.uU@CPe'I`>APXPX`XZ<AHPQBqIJ>?PXNklNH._ klBNeqTUH.H._JHH._`_JPSBADI`@1I`>AHEhkaXdklBq( R(H.IJ>?@[ep§@a_Ye^H.I$kWPQOXdXZH.Ht#OS<A_J<` "+M

*

|>AHA_JPSBNPSNklO@as[wJHWI`PQuUH%PQB IJ>APX(XZHIJPS@CB6PXIJ@@C_`R1<AOklI`H IJ>?Hs@a<ABNeAkW_Jb PQB[IJH.Da_$kWOkWBNe@C_`_JHWX@aBNe^PSBADqklOSDCH.sA_+klPH*)U<NkWIJPQ@aBNX:@WR klB[bNkl_`[email protected] XyPSBk@a_JRX<APQI+klsAOSH1@C_IJ>AH!Eh kWAAOSP.kWIJPQ@aB9MqH%kC.&@CR(AOSPX>IJ>APXy@as[wJHWI`PQuUH3s[bq@CRsAPQB?PQBAD(VH.OQO\HXZI$kWsAOQPXZ>AHejdg= A_J@[&HWe^<A_`HWX1j_JH.sAs?Pk%klBNe #@aR(PQB?DC<AH %F=C*5^0VdPQI`>BAH.V _`H.A_`HWXZH.B[I$kWIJPS@CBNX#@C_IJ>AH<ABNeAkWR(H.B[I$klOXZ@aOQ<AI`PQ@CBXn@WhI`>A_JH.H&e^PSR(H.BNXPQ@aBNklONOQPSBAHWkW_dH.OkaXZIJPPQIbUM

!#"%$&!'

~ B IJ>APX(XZHIJPS@CB9FVH@a_JR<AOkWIJH%IJ>AHLsN@C<ABNe?kl_Jb6PSB^I`H.DC_+klOYH*)U<NklI`PQ@CBX @a_%kXZPSR(AOSH(R%klB[bUNkl_`IJP&OQHA_`@CsAOSH.R I`>NklIYPQB[uU@COSuUHWXkWB PQB NBAPQI`H0sAPQBNe[H._dVdPQI`>)( Pe^HiB^I`P.kWOpNkW_JIJPOQHX=NH._ZHIJOSbs@aBNe^HWeqIJ@%IJ>AH3s?PQBNe^Hi_.M|>APXA_`@CsAOSH.R PX0XZH.OSHWIJHe @C_e^HiR(@CBNXIJ_+klIJPS@CBA<?_JN@UXHWXM!|>AH._JHkl_JHYBA@ @CBN&H.AIJ<klOA@asNXZI+kCOSHWX=I`@:Hi^I`H.BNe^PSBAD0@C<A_kl?A_J@kC$>(IJ@[email protected] XI`>NklIPQB[uU@COQuUHNBAPSIJHsAPQBNe[H._$X&FNe^P #H._`H.B^INkl_`IJPOQHWX&FAkWBNe e^His@aBNe^PSBADAM|>?HA_`@CsAOSH.R @W,PQB[IJHi_JHWXIPX@C_`R1<?OklI`HWe3PQB I`>A_JHiH0XZI`H.NXMEPS_$XZIiFAVH@C_`R1<?OklI`H0klBvkl<AU

PQOQPkl_`b3s@a<ABNeAkW_JbU\uakWOQ<AH0A_`@CsAOSH.R @a_dk NB?PQIJHYsN@[e^b%VdPQI`>*( H.RsNHWeAe^HWekl_JI`POSHWXM,rUHW@aBNepFVH@C_`R1<AOklI`H:I`>AH0sN@C<ABNe?kl_JbLPQB[IJH.Da_$kWO9H*)<klIJPS@CBNX@C_YIJ>?H NBAPSIJH0sN@[e^bkWBNe PQB[IJ_`@[e^<NH0BAHiV_JH.A_`HWXZHiB^I+klIJPS@CBNX@C_YIJ>?Hy<?BNeAklR!H.B^I+klOX@COQ<?IJPQ@aBNXMYEPSBNklOSOQbUF[VH e^H._JPSuUHIJ>?HsN@C<ABeAkl_`b PQB[IJH&DC_$kWO,H*)<klIJPS@CBNX#@a_:I`>AHPQB NB?PQIJH3sN@[e^bM(EA_J@aR BA@WV7@aB9FVH3<NXZHIJ>?H(X<AsNXJ&_JPQ?I$XIJ@Le^H.B?@CIJHmkl_JI`HWXZPklB!@[email protected][I$XhkWBNe(X<ANH._$X`_JPSAI$XhIJ@e^HiBA@CI`HkWIJIJ_`PQsA<?IJHWX&MhEA@C_,N@UXPQIJPS@CBuUHWI`@C_$X&FCVH<NXZHPSB^I`H._$$>NkWBADCHWkWsAOQb3I`>AHPSBNe^P&PklOkWBNe%e^PQ_`HWI@CB[uUH.B[IJPS@CBNX&§@a_dH.?klR(AOSHUFCI`>AHBA@aI$klI`PQ@aBNX

Page 9: A Fast Solution Method for Three-Dimensional Many-Particle

Γ

Ω

n n

Ωm

m

EPQDC<?_JH4,t sAPSBNe^H._Y@aB^I+klPSBAPQBAD3Pe^H.B[IJP.klONNkW_JIJPOQHXM

klBNe kW_JH0H*)U<APQuakWOQH.B[I.Mm,@aBNXZPe^H._4k BAPQI`H:I`>A_JHiHe[PQR(H.BXZPQ@aBNklONOSPQBAHWkW_dH.OkaXZI`P#sN@[e^b%@C_JR!HWe1s[b%ks?PQBNe^Hi_ kWBNe

( H.RsNHWeAe^HWeqPe^H.B[IJP.klO,NkW_JI`POSHWX NH._ZHW&IJOQbqsN@CBNe[HWe IJ@qIJ>AH%sAPSBNe^H._%EPQD?M W$M He^H.BA@aIJH I`>AHX<A_ZkaH @WIJ>APXsN@[e^b s^bWFI`>AH PSB^I`H._ZkaH%s@a<ABNeAkW_JPSHWX1s[b FdklBe IJ>AH.PS_<ABAPQ@aB9F

#&% F9s^bMkWR(H X@CBXZI$kWB^I+X@ldI`>AH3sAPSBNe^H._kl_`He^H.BA@aIJHWes[b qkWBNe kWBNeIJ>AHY@a_J_`HWXZN@CBe^PQBAD 2@aPX`XZ@CB Xh_$kWIJPQ@0s[bpM,EA@C_hI`>AHdNkl_`IJPOQHWX&FCIJ>?HY@C_`_JHWX@aBNe^PSBAD@CBNXI$kWB^I+Xkl_JH4R%kl_`TUHWeVdPQI`>(IJ>AHykCXZI`H._JPXZT!XZ<ANH._+XJ_`PQAIiM,|>AHsN@[e^b(PX=OS@UkCe[HWe1s[b(ke^PXAOkaH.R(H.B[I NHiOe

A_JHWX`_JPSsNHWe1@CBWM|>?He[PXZ?OkC&H.R(H.B[I NH.Oe#PQB0I`>AHhsN@[e^bUF $F.PXpIJ>AH,X@COS<AIJPS@CB0@l?IJ>AHh@aOQOQ@WVdPSBADdsN@C<ABNe?kl_JbU

uCkWOQ<AH0A_`@CsAOSH.R

"' " ' " " =, ¡ -

"' " ! ' " " h ¡ # % - C

, " -

N ¡ klBNe # N ¡ - 4$ /

~ BLU MQMQM de[H.BA@CI`HWX=I`>AHwJ<AR(!klOQ@aBADI`>AHYBA@C_`R%klONuUHWI`@C_iF +F[kCXe[H$NBAHWePQB!EPQDa<A_JHCM

Page 10: A Fast Solution Method for Three-Dimensional Many-Particle

g)U<NklI`PQ@aBNX: kl_`H:H )<APSuaklOQHiB^I,[email protected] @l( PSB^I`H.DC_+klO9H*)U<NklI`PQ@aBNXM,|>AHN_+XZIH*)<klIJPS@CBqPQBqIJ>APXXZb?XZI`H.R8PX0e^H NBAHWe@CBqIJ>AHsAPSBNe^H._sN@a<ABNeAkW_JbUF FhklBe'HWka >q@l IJ>AH_JH.R%kWPQBAPSBADH )<NkWIJPS@CBNXPXde^H$NBAHe @aB I`>AH0s@a<ABNeAkW_Jb%@l=kNkl_`IJP&OQH

= " - #, " #he # " - #h " #=e #

( " - #h # " #he # ( " - #h # " #h `.¡U

klBNe

" -# " # e # ( " - # # " # e # - -$ -/1/0/1- ( / `CW

~ BLJi¡kWBNe Jal klBNe(PQBIJ>?H#_`HWXZI@lIJ>?H#NkWNH._dPQI,PXhPSR(AOQPSHWe:I`>NklI PX,klB%kW_JsAPSIJ_+kl_JbN@CPSB^I@CB I`>AHXZ<A_kCH0@aB Vd>AP >(I`>AH0PQB[IJH.Da_$kWO9H*)<klIJPS@CB%PXde[H$NBAHWeM|>?H3TUH._JBAHiOX ( " - #hYklBe " - #hPQB[IJ_J@[e[<NHWe PSBJi¡YkW_JH3I`>AH:<ABeAklR(HiB^I+klOXZ@aOQ<^

IJPQ@aBNX@l IJ>A_`H.H e^PQR(HiBNXZPS@CBNkWOOSPQBAHkl_dH.OkCXZI`PPSIbUMXZ<NkWOQOSb1I`>AH.b%kW_JH0H.[A_JHWX`XZHWe%kaX

( " - # ,

J * %# * ! 5 " $ " " * %# * + `.U

klBNe " - #h8 h` * # * - $ " #h 6 " " $ # " J$& #h

!h` * # * ¢ $ " " $ J$$ #h - ` !

Vd>AH._JH " PX1IJ>?HqXZHW@aBNe6_$klB?T Pe[H.B^I`PQIbI`H.BNXZ@a_.M |>?HTUH._JBAHiOX( " - # klBNe " -#

PQB[IJ_J@[e[<NHWe(PQBq`Cl .kWB sNH@asAI$kWPQBAHWe1_`@CRJi klBNeqJ !U@aBNHkWBNe %kl_`H_`H.AOkaHWe(s[b klBNe FA_JHWXNHWIJPSuUH.OQbM|>?H:TUHi_JBAH.OX ( " - #,klBNe " - # .kWBsNHklOXZ@H.[A_JHWX`XZHWe%kaX

( " - #h, " * %# * ,

* %# * " `$5[

klBNe " -# , " * %# * #

* %# * # " - `

Page 11: A Fast Solution Method for Three-Dimensional Many-Particle

Vd>AH._JH " ,= J ! 5 N " " - ` U

h h` - `la

" , h` J $ " " [J " " - `.U

klBNe ,

h` / ` U

~ BvJ .$+FI`>AH0XZb[Rs@aO e^H.BA@aIJHWX,IJ>AHkl_JI`PkWOe^H._JPSuaklIJPSuUHdVdPQI`>(_JHWXNHWII`@ M #@aIJHI`>NklIPQB`$5^ klBe'` ,IJ>AH0IJHi_JR%X,e^H.NH.Be^H.B[I@CB kl_JHXH.Nkl_+klI`HWe1_J@CR/IJ>?@UXZHe^HiHiBNe^H.B[Id@CB #H.A&H.AI=@a_#I`>AH0>Nkl_`R(@CB?P#TUH._`BAH.O * %# * % M

Y@lVqOQHiI9<NXHiuCkWOQ<NkWIJHYJi¡klBNeJal<ABNe^H._9I`>AH@CBe^PQI`PQ@CBXI`>NklI I`H.BNeAX9IJ@dPSB NBAPSIbkWBNe

, " " FVd>?H._JH " PXk@aBNXZI+klB[IIJHiBNXZ@a_.M,|9@3e^@3I`>APXhVHN_+XZId&@CBNXPe^H._ klB%kW<A[PQOQPkl_`[email protected]/@C_YIJ>AH NBAPQI`H0sAPQBNe[H.__J@aR¥Vd>?P$>I`>AH0Nkl_`IJPOQHWXdkW_JH_JH.R(@WuUHWeklBe PQB%Vd>AP$>%IJ>AHOQPQB?HWkl_ e^PXAOkaH.R(H.B[I NHiOe ,= " " PX,_JHWkWOQP1%iHWes[b(A_`HWXJ&_JPQs?PQBAD , " " @CB%IJ>AHH.B^I`PQ_`HsN@C<ABNe?kl_Jb M|>AH @C_`_JHWX@aBNe^PSBAD%XZI`_JHWX`XdklBNe I`[email protected]_`H3e^H.B?@CIJHes[b

" kWBNe # Fp_JHWXNHWIJPSuUH.OQb§ # " " M0|>?H3PQB[IJH.Da_$klOH*)U<NklI`PQ@aB @C_I`>AHkl<A[PQOSPkW_Jb%A_J@asAOQHiRNkl_+klOQOSH.OX`.¡U

" -# " #=e # " - # " #he #

( " - #h # " #he # ( " - #h # " #= / C¡U

!#BHq\C¡U3PX3X<AsAI`_$kC&IJHWe_J@CR J.¡U$FdIJ>?HPSB^I`H.DC_+klOX@luUHi_ %e^PX`klANHWkW_%kCX (I`H.BNeAX3IJ@PQB NB?PQIb(kWBNeVH0@CsAI+klPQB

" - #h " #he # ( " - #h # " #he # ¡ - ^W

Vd>AH._JH

, = klBe # = # # 2/ CU

a

Page 12: A Fast Solution Method for Three-Dimensional Many-Particle

~ B%I`H._JR%X,@lIJ>AH0NH._`IJ<A_`sNklI`PQ@CB,NH.OeAXF9`ClPX_`H.Vd_JPSIJI`H.B(PQB%I`>AH@C_JR

" - #h " #he # ( " - #h # " #he #¥ 4!U

VdPQIJ>

, " " " " -#= pe # " ( " - #h$ #=e # / 5[

~ Bq4! klBe \*5^PQIdPXPQR(AOSPQHWeI`>NklI ¥_`<ABNX=_`@CR80IJ@ (M~\BI`>AH3_JH.R klPQBe^H._Y@lIJ>AHNkWHi_.F9VHkW_JH@aBNH._`BAHWe VdPQI`>qXZ@aOQu[PQBAD(PSB^I`H.DC_+klOhH )<NkWIJPS@CBNX

\^W4klBNe \"!$M |>AH3R(@UXI$>NkWOQOQHiBADCPSBAD(NkW_JI0@lIJ>?PX#I+kCXZTPX#IJ@ XZ@COSuUH%\[lYVd>AP$>iklBqsH >NkW_$kaIJH._`P1%.HekaX#k!@C<A?OQHWe(A_`@CsAOSH.R @C_ ( Nkl_`IJP&OQHWX&MEA@C_4IJ>APXnA_`@CsAOSH.R I`>AH<NXH@W,IJ>AHEh PXd_JPSIJP.klOMh~\BL@CB[IJ_+kCXZIiF9\"!_JH.A_`HWXZHiB^I+X ( _JHiOklI`PQuUH.OSb3XR%klOSOp<ABN&@C<AAOSHWe(A_J@asAOQHiR%X@C_YVd>AP$>%k@CB[uUH.B[IJPQ@aBNklOjdg, PXdkCe[H*)<klIJHM

% " ! "%$&!'

~ BIJ>APX9XZHW&IJPQ@aB9F[\^WklBe1\"!kW_JH @CB[uUH._JI`HWe#PSB^I`@kWOQDCHisA_$kWPhH*)U<NklI`PQ@CBX@aOQOQ@WVdPSBADXI$kWBNeAkl_+eNBAPQI`H0H.OQH.R!H.B^I e^PX`_JH.I`P1%WkWIJPS@CB(I`HW >?BAP )U<AHWX&M|>?H BAPQI`H:HiOQH.R(HiB^Ih_`H.A_`HWXZH.B[I$kWIJPS@CB%@l PXskCXZHe@CBLIJ>AH0R%kW

N, !)' p /

~ BL\ (PXhIJ>?H#N@UXPQI`PQ@CB!uUHWIJ@a_,PQBIJ>AHYR%kCXIJH._,H.OSH.R(H.B[IXZkCHUF PX=IJ>?H#N@UXPQI`PQ@CB!uUHWIJ@a_@lYIJ>?H(BA@[e^H @l4IJ>?H1H.OSH.R(H.B[IpF,klBNe p4PX0IJ>AH%X>NklNH1<?BNIJPS@CB@WYIJ>AH!BA@[e^H MLEA@C_PXZ@aNkl_+klR(HiIJ_JP NBAPSIJHHiOQH.R(HiB^I+XFCI`>AH0R%klNX p kWBNe # d pR!PQR(P\ $F

p p !)' p klBe # p !

2' # p / U

6>AH.Bq\ klBNeq\UdkW_JHXZ<?sNXZI`PQIJ<?IJHWe%PQB^lVH0@asAI$kWPQB !

)' " -#=p " p #=[e

!)' ( " -#=p # " p #hpUe% ¡ / a

i

Page 13: A Fast Solution Method for Three-Dimensional Many-Particle

~ BL\UC,I`>AH#PSB^I`H.DC_+klOX=kW_JHY@WuUH._,IJ>AHYR%kCXIJH._,H.OSH.R(H.B[IXZkCH klBNe #,pPXhIJ>AH4KUkC@asAPkWB@C_`_JHWX@aBNe^PSBAD I`@ $M #<AR!H._JP.klO9PSB^I`H.DC_+klIJPS@CB@W\ad@COQOS@WVHWe s[bq@COSOQ@[.kWIJPS@CB klI0IJ>AHBA@[eAklON@CPSB^I+X b[PQH.OeAX=I`>AHklOQDaH.sA_+klPH*)U<NkWIJPQ@aB

!

2' ' " - #= " #h !

2' ' ( " - #= # " #= 6¡/ CU

~ B CU_J<ABNXY@luUHi_IJ>AH3PSB^I`H.DC_+klI`PQ@CBN@aPQB[I$X@l IJ>AH3R%kaXZI`H._#HiOQH.R(HiB^I klBNe PXyI`>AH3PQB[IJH&DC_$kWIJPS@CB VH.PQDa>^I klI MY@aBLXZ<AsNXIJPSIJ<AI`PQ@CBL@l J 5l+ +Fh\a .klB%sNH_`H.Vd_JPSIJI`H.B%kCX

" !

)' ' * # * " #h #

!2' '

* %#= * " $ #, " #h " !

2' ' * %#= * # " #h

!)' '

* %#, * # " " #, ¡ / 4U

t#DUkWPQB9FWBA@CI`HI`>NklIPQB\"IJ>AH,I`H._JR%Xe^H.NH.BNe[H.B^I9@aB kW_JHXZH.kl_$kWIJHWe_J@aR I`>A@UXHde^HiHiBNe^H.B[I@CB"#qH.A&H.AI=@a_#I`>AH0>Nkl_`R(@CB?P#TUH._`BAH.O * %# * % M|>?He[PXJ&_JH.I`P1%.HWe3@C_JR X=@ln\4!UdkWBNe 5[dkW_JHXZPSR(PQOkl_,IJ@L\a

!

)' ' " - #= " #h !

)' ' ( " - #, # " #h - !C¡U

klBNe

, " " " !)' ' " - #= #h

" !)' ' ( " - #, #, / !^W

|>AH%VdPe^H(>klI3PQR!AOQPSHWX#IJ>klI3IJ>?H X<AR%X0PQB[uU@COQuUH3I`>AH(BA@[e[HWX1kWBNePQB[IJH.Da_$kWIJPQ@aB @aPQB[I$XI`>NklIsHiOQ@CB?D1I`@k1XPQBADaOQH0Nkl_`IJPOQHUM

!

Page 14: A Fast Solution Method for Three-Dimensional Many-Particle

tdIpI`>APX9XI$kWDCHUFCPSIpPXp<NXZH<?OUIJ@d_JHiVd_JPQI`HY\4UNIJ>A_`@C<ADa>%![lNPQBkYXZb[R1sN@aOQP9R%kWIJ_JPS@C_`RM9|9@IJ>APX HiBNeVHn@a_JR e^PSR(H.BNXPQ@aBNklOAB?@[eAklOuUHWIJ@a_$X FF^kWBNe%@C_`DUkWBAP1%.He1@CB!k:NkW_JI`POSH\s[bUNkl_`IJP&OQHskCXZPXF?XZ@0IJ>NkWI,HWkC$>3uUHWI`@C_PX@C_`R(HWe:s[b ( sAOS@[ T?XM,m,@CBNXH*)<?H.B^I`OQb0IJ>?H#@aHi_$klI`@C_+XkCI`PQBAD0@CB klBNe4kl_`H_JHiA_JHWXH.B[IJHWe3kaX /R%kWIJ_JPHWXp@C_`R(HWes[b ( ( sAOS@[ T?XM,|>[<NXVH0_JH.Vd_`PQI`H3\" kCX

!CU

klBNeq!C¡UdkWBNe !^WdkaX : / !4!U

g)U<NkWIJPQ@aBNX0!CUdkWBNe !4!Ud.kWB%sH@aR1s?PQBAHWePQB[IJ@

% % ! 5[

@C_4XZPQR!AOQb / !

t#OX@3PQIdPX<XZH<AO9I`@3_JH.Vd_`PQI`H1 ! PSB(IJ>AH0_`HWXZPe^<NklOp@a_JR

, ¡ / !U

|>?H0XZI`_J<NI`<A_JH4@lIJ>?H#R%kWIJ_`PHWXPQB !C klBNe !4!_JH HWI$X,IJ>AHYA>[bAXP.Xh@WIJ>AH?_J@Cs?OQH.RM|>AH R%kWIJ_`PHWX kWBNe kl_`He^HiBNXZHLklBeIJ>?H.PQ_3@*# e^PklDC@aBNklO0sAOS@[ T?X_JHiA_JHWXH.B[I3NklPS_ZVdPXZHNkl_`IJP&OQH,PQB[IJH._+kCI`PQ@aBNXpVd>APQOSHI`>AHe^PkWDC@CBklOAsAOS@[ T?X9_JH.A_`HWXZHiB^I9NkW_JI`POSHXH.O PQB[IJHi_$kC&IJPQ@aBNXM|>AH@ #Ne[PklDa@CBNkWOsAOS@[ T?X1sNHW@aR(H B?H.DCOSPQDCPSsAOQH(@a_VH.OQOXH.NkW_$klI`HWekl_JI`POSHWXFYkWBNe IJ>AHiPQ_%e^H.k.b_$klI`HWX@COQOS@WV e^PS_JHW&IJOQb_`@CRJi klBNeqJ!$M|>AH0R%kWIJ_JPHWX kWBNe kW_JH0sAOQ@[$TUe^[email protected]<NXH3IJ>AHib_`H.A_`HWXZH.B[I ( <ABN&@C<AAOSHWe?_J@Cs?OQH.R%X&MdEA@C_Pe^H.B[IJP.klO9NkW_JI`POSHWXFAI`>AH3e^PklDa@CBNkWOsAOQ@[$TAXPSB%HWkC$>%R%klI`_JPQkW_JH0IJ>?HX`klR(HM

, *

~ BIJ>APXhXZHIJPS@CB9FCVH@a<AIJOSPQBAH k#A_J@[&HWe^<A_`H@C_,X@COQu[PQB?D ! <NXPQBAD( 9R(HiR(@C_`bYklBNe(( @CNH._$kWIJPS@CBNX&M,|>APXpA_`@[HWe^<?_JHd&@CBNXPXZI+X@WAI`>A_JHiH&@CR(N@CB?H.B^I+X$hP e^H@CR(N@XZPQI`PQ@aB:@W kWOQ@CB?D

5

Page 15: A Fast Solution Method for Three-Dimensional Many-Particle

IJ>AH(_`HWXZPe^<NkWOdG_Jb[OQ@Wu XuUHWI`@C_$X F F + F - FhH.I+CMQF, PSP YR%kWIJ_JPSuUHWI`@C_R<AOSIJPQ?OQPiklIJPS@CBVdPQIJ>%I`>AHEhFAkWBNe PQPQPhA_`HW@aBNe^PQI`PQ@aBAPQBAD?M|>?HdXZ@aOQ<AI`PQ@aB3uUHWIJ@a_ PX9@CsAI+klPQB?HWey@aOQOQ@WVdPSBAD#I`>AHDaH.BAH._+klOSP1%.HWeYR(PQBAPSR1<?R _JHXZPe[<NklO?klOSDC@l

_JPQI`>AR8]cYgrA@ldr[kakCeLklBNerU>^<?OQI %` C +M#|>?PXYklOQDa@C_`PQIJ>?R _JH*)U<APS_JHWXd@aBNXZI`_J<N&IJPQ@aB @WIJ>AH_JHWXPe^<NkWOAG_`b^OS@lu XNuUHWI`@C_$X&M|>AH IJ>3G_`b^OS@lu4uUHWIJ@a_ PX9&@CBNXIJ_J<IJHWe0s^b0PSBNe^<NI`PQ@aB

% VdPQI`>

/

E@aOQOS@lVdPSBAD%!*5^+FAIJ>APXR1<AOSIJPSAOQP.klI`PQ@CB3PXNH._Z@a_JR(HWe!kCX

% % - !a

VdPQIJ> % % % /

!#sXZH._`uUH3IJ>NkWI:I`>AH3IJHi_JR % .klBqsNH(&@CR(A<?IJHWe PSB( 4XZPSBNH3sN@CI`> klBNe kl_JHsAOQ@[$T e^PklDC@aBNklOM #@WV @aR(HWX#I`>AH!_J<N&PklO XZIJHi 'klAAOSP.kWIJPS@CBq@ldI`>AH!Eh I`@ !a$M%|9@H.^?OklPSBIJ>?PXXIJH. OSH.I0<NXDC@qsNkC$TqIJ@ \"kWBNe_JH.AOkCH!IJ>AH._`H(IJ>?H(BA@[eAkWOYe^PXAOkaH.R(H.B[I$XVdPQIJ>qIJ>?H1BA@[e?klOduaklOS<AHWX@l % klBe IJ>AHBA@[eAkWOdIJ_+kCI`PQ@CBXyVdPSIJ>qVdPQI`>IJ>AHBA@[eAkWOuaklOS<AHWX@l % ML|>AH1@C<?_XZ<AR%XPSB 40>kuUHIJ>?H%XJkWR(HXZI`_J<NI`<A_JH klBNe I`>AH._`H@C_`H(PQI0PXXZ<PQH.B[IYIJ@%H.[AOkWPQB>A@WVIJ@Lkl?AOQbI`>AHEhI`@I`>AHN_+XZI0@CBAHM3|@L@aBABAHW&IIJ>APX0XZ<AR¨VdPQIJ>IJ>AH(DaH.BAH._+klOSP1%.HWeR%kWB^bU\sN@[e^bqH.OSHWIJ_`@UXI$klI`P.X0A_J@asAOQHiRe^H NBAHWe PSB rUHW&IJPQ@aB ^FhVH(Pe[H.B^I`P bIJ>AH%&@COQOS@[.klI`PQ@aBN@aPQB[I$X VdPSIJ> IJ>AH!I$kl_`DCH.I3N@aPQB[I$XF,IJ>?H(PQB[IJH.Da_$kWIJPQ@aBN@CPSB^I+X #= 4VdPQIJ>IJ>AH%X@C<A_+H(N@CPSB^I+XFhI`>AH(H.[A_JHXJXZPS@CBPSBIJ>?H%<A_`OQbsA_+kC$TUH.I$XVdPQI`>IJ>?H(DCH.B?H._$kWOQP1%iHWe $>NkW_JDCH #h

\F[IJ>AHX<AR VdPQI`> I`>[email protected][IJPklO $FNkWBNe(IJ>AH@CNH._$kWIJ@a_ 0 VdPQIJ>%I`>AH0DCH.BAHi_ZklOQP&%.HWe@a_$H3@aHi_$klI`@C_ $2-2 M !YsNXZHi_JuUH3IJ>klI0IJ>AHDCHiBAH._$kWOQP&%.HWe $>NkW_JDCHPX0k%XZH@CBNe_$klB?T I`H.BNXZ@a_3kWBNe IJ>?H._JH@a_JHPQI>NkCX0IJ@ sNH(_`H.A_`HWXZH.B[IJHe's[bBAPQB?H1@a_$e^PSBNkl_`b$>Nkl_`DCHWX&M |>[<NXIJ>AH _$XZI XZ<AR PX@aR(A<AI`HWes[b1PSB^uU@aT^PSBADIJ>?H0Eh¥B?PQBAHYIJPSR(HWXMEA@C_IJ>AHY_JHiR%klPSBAPQBADyXZ<AR%X&FIJ>AH0E PX=PSB^uU@aTUHWe3IJ>A_`H.HIJPSR(HWX@a_dIJ>AHyXZHW@aBNe%klBNe(I`>APQ_+e(XZ<AR XFAkWBNe(@CBN&H#@C_dI`>AHOkaXZIXZ<ARF?XZ@IJ>NkWIdIJ>AH0I`@CI$kWOB[<ARsNH._,@l,I`PQR(HWXhI`>AHEh PXPQB[uU@CTUHWe3PXdXPQ[IJH.HiB9M

Page 16: A Fast Solution Method for Three-Dimensional Many-Particle

EA@a_OkW_JDaH0klOSDCH.sA_+klPnA_`@CsAOSH.R%X&FCOQPSTUHdIJ>A@XZH@C_`_JHWXN@CBNe^PSBADI`@R%kWB^bUNkl_`IJP&OQHdA_`@CsAOSH.R%XFA_JHW&@CBNe^PSIJPS@CBAPSBAD3PXdPQR(N@a_JI$kWB^IsNHW.kl<XZHPSIYklOSOQ@WVYXn@aBAH0IJ@_JHWe[<NH0IJ>AHB^<?R1sNH._,@l _JH*)U<APQ_`HWePQIJHi_$klI`PQ@aBNXM 2h_JHW&@CBNe^PSIJPS@CBAPSBAD PSB^uU@aOQuUHWX0R1<?OQIJPSAOQP.klI`PQ@aB@l0sN@aIJ> XZPe^HWX1@W! 3s[b k R%k.IJ_JPS(FIJ>AHLA_JHW&@CBNe^PSIJPS@CBAH._iF=Vd>?P$>PX1kaX(OS@UXZH%I`@ % kCX1N@UX`XZPSsAOQHUM qHL >A@U@XZHIJ@sHqsAOS@[ TUe[PklDa@CBNkWO\FVdPSIJ> HWkC$>6sAOQ@[$T sNH.PSBAD H )<NkWOIJ@ I`>AH PSB^uUHi_$XZHL@l IJ>AH @C_`_JHWXN@CBNePQBAD!e^PkWDC@aBNklOsAOQ@[$T%@l MEA@a_JR%kWOQOQbF^IJ>?PXA_`HW@CBe^PQI`PQ@CB?H._dPXdT^B?@lVdBvkCXdsAOQ@[$TUKUka@CsAPMdEA@C_R%klB[bNkW_JIJPOQH,A_`@CsAOSH.R%X&F.IJ>AHdsAOS@[ TU KUkC@asAPAA_`HW@CBe^PQI`PQ@CB?H._>kCX,kOSHWkl_A>[bAXP.kWO^R(HklBAPSBADPQIdH.?kCI`OQb_JH.?_JHWXH.B^I+X=BA@aB^PQB[IJH._+kCI`PQBAD3NkW_JI`POSHWXM|>AH._JH&@C_JH0I`>AH0R%klI`_JPQ _JH.?_JHWXH.B^I+XIJ>AH3PSB^I`H._ZNkl_`IJPOQHPQB[IJH._+kCI`PQ@CBX#@CBAOSbUM #HiI$klPSOX@lPSBN@C_`N@C_$kWIJPSBAD I`>AH3A_JH@CBNe[PQIJPS@CBAHi_PQB[IJ@IJ>AH3]4c#gr!.kWB sNH@C<ABNe%PQBr[kakCe kWBNe'r>^<AOSI %3` a+M

|>AH0DC@UkWOX@C_d@a<A_dNkl_+klOSOQH.OPQR(AOSH.RH.BUI$kWIJPQ@aB1@WhI`>AH0Eh8kWOQDa@C_JPSIJ>AR VH._`HH PQH.Bb3kWBNeXJ.kWOkls?PQOQPSIbkWBNe1A_`@Wu^PXZPQ@aB@a_,R(HWXZ>kaeAkl?IJPQu[PQIbPQB!kXZH@CBNe3A>NkaXZH@WkWOQDC@a_JPSIJ>AR/kWBNe(@[e^H_JH$NB?H.R(H.B[I.Md|>?HWXZH3Da@UklOX#VHi_JH3kWIJI+klPQB?HWe s^b%sA<?PQOe[PQBAD(I`>AHkl_$kWOQOSH.Oh&@[e^H3@CBI`>AH(r[.kWOkls?OQH#b^BklR(P #PXIJ_`PQsA<AI`HWe(t#_`_$k.br Yt:eAklI+kR klBNkWDCH.R(HiB^I,PQB[_$kCXIJ_`<NIJ<?_JH\g,e^V kl_$e?X#kWBNej_J@WVdBAHUF94[§9g,e^VdkW_$eAX&F4UC$MY|>AHr Yt A_`@Wu^Pe^HWXBA@a_JR%kWO9kl_`_$k.bLXZH.R klB[IJPiX @C_4klB[bkl_J_+k.b@lN@as[wJHWI+Xh Vd>AH._`H,HWkC$>0R%k.bd>Nk.uUHkde[P #H._JHiB^IXIJ_J<IJ<A_`HCNVd>AP >PXe[PXZI`_JPSsA<AIJHekC_`@UX`XR1<AOSIJPSAOQHA_J@[HXJXZ@a_$XdkWBNe1Vd>AP >R%k.bH.[NkWBNe(klBe(@CB[IJ_+kCIVdPQI`>1A_`HWXZHi_JuaklI`PQ@CB!@lOQ@[iklOQPSIbklBNe VdPSIJ>A@a<AI0@a^b[PQBADIJ@@COQOS@WVkCeAkWAIJPSuUHR(HXZ>%_JH$BAH.R(H.B[I.M |>AH1r #t7V kCX#@a_JPQDaPQBNkWOQOSbe^H.uUH.OS@CNHWe#I`@YXZ<AAN@a_JIh>A[kaeAklAI`PQuUHNBAPSIJHhHiOQH.R(HiB^IAR!H.IJ>A@[e?X\g,e^V kl_+eAXkWBNej_`@WVdBAHUFU "[§g=e[Vdkl_+eAXF94UC@a_dVd>AP$>(PQId?_J@Wu^Pe^HWXeAklI+kR%kWBNklDaH.R(H.B[I@a_dIJ>AH0DCHi@CR(H.I`_JbeAkWI$k3@lIJ>AHkCeAkWAIJPSBADNBAPQI`HH.OSH.R(H.BUIR(@[e^H.OM|>APX9_JH*)U<APS_JHWXXZ<A?@a_JIh@C_H %&PQH.B[Ip>NkWBNe^OSPQBAD@We^b^BklR(Pkl_J_+k.bAX#@W XZR%kWOQO@Cs[wJHWI+XnVd>?P$>PXdIJ>AH<ABeAklR(HiB^I+klOeAkWI$k(R%klBklDCHiR(H.B[[email protected] @a_IJ>AHNkl_+klOQOSH.O9E PSR(AOQHiR1H.BUI+klI`PQ@CBM|>?HPQR!AOQH.RH.B[I$kWIJPQ@aB#@WNIJ>?H#r YtPXpsNkCXHWe<AN@CB3<XZH@WXZkCH NOSOQPQB?DY<A_JuUHX9klBNe3XH.^

Page 17: A Fast Solution Method for Three-Dimensional Many-Particle

EPSDC<A_`H 2hkW_JIJPS@CBAPSBAD@lIJ>AH@aR(A<AI+klIJPS@CBNkWOe^@CR klPQB!sNkCXHWe@aBk3XZkCH NOSOQPQB?D<A_JuUHM

kl_$kWIJPS@CB@l kaeAe^_JHXJX#R%kWBNklDaH.R(H.B[IPQBNe[H.^PSBAD^h_`@CR¨XZI`@C_+klDCHR%kWBNklDaH.R(H.B[I 2hkl_+kCXZ>kl_0kWBNej_J@WVdBAHUF[ "C+MnrNkCH& NOSOQPSBAD&<A_JuUHWXhr[klDklB9FN4 5[R%klIJ>AHIJ>A_`H.H e^PQR(HiBNXZPS@CBNkWOCA>[bAXP.kWOXZNkaHI`@%k@aBAH e^[email protected]_$klOkW<AIJ>?@C_$X>Nk.uUH0A_`H.u^PS@C<NXOQb(<XZHWeLXZNkaH NOQOSPQBAD<A_JuUHX=@C_PSR(AOQHiR1H.BUI+klI`PQ@CB!@l,NkW_$kWOQOQHiOI`_JH.H@[e^HXZr[kWOQR(@aBLklBNe kl_`_JH.B9F "C[§ rklOQR!@CBklBNe kW_J_`H.B9F, "4!^§rPQBADaH._.F, " $M(t/XZPSR(AOSHYPQOQsNH._`I X#XNkCH& NOSOQPSBADL<A_JuUHPX0XZ>?@lVdBqPQBEPQDa<A_JH[M(~\I_JHXZH.RsAOQHXdk B?HW T[OkaHVdPSIJ>qIJ>AH!<A_`uUH3kaX:IJ>?H.kl_`_JPQHi_0klBNeqIJ>AHHiOQOX4kCXIJ>AHsHkCeAX&MtX`XZPSDCBAR(HiB^I,@lVHiPQDC>[I$XI`@3IJ>AHyH.OQOX=H.BklsAOSHWXdXZPSR(AOQHYkWBNeH %PSH.B[I=kl_JI`PQI`PQ@CB?PQBAD3@WIJ>AHI`>A_JHiHe[PQR(H.BXZPQ@aBNklOXZNkaH0VdPQIJ>!A_JHWXH._JuakWIJPQ@aB(@lPQI+X=DaH.@CR(HiIJ_JP.klO?OQ@[.kWOQPQI`PQHWX,PQBIJ>AH0@aBAHe^PQR(HiBNXZPS@CBNkWO,XZNkaHUM|>?H(R%kl?APQBAD%I`@qXZIJ@a_$kWDCHLe^PXIJ_`PQsA<AI`HWeqkC_`@UXJX3k%Nkl_+klOSOQH.O@aR(A<AI`H._kl_$$>APSIJHWI`<A_JHkWBNe ka.HWX`XZPSBAD1@WIJ>AHeAkWI$k1PXkCi@CR(?OQPX>AHWe(IJ>?_J@C<?DC>q<NXZH3@Wdk>APQH._+kl_+ >AP.klOPQBNe^Hi XNkC&H3Vd>AH._JHIJ>AH3TUHibAXdkW_JHe^H._`PQuUHWe1_`@CR I`>AHN@XZPQI`PQ@aBNX@lIJ>AH&H.OQOX@CBIJ>?HXZNkaHNOQOSPQBAD!<A_`uUHUM0t#OQO9kl_JI`PQI`PQ@CB?PQBAD?F9e^PXIJ_`PQsA<AI`PQ@aB9F9klBe @aR(R<ABAPS.klI`PQ@aB%kl_`H3A_J@Wu[Pe^HWe s[b%IJ>AHr YtI`_$kWBNXZNkW_JH.B[IJOSb(IJ@IJ>?H:<XZH._iM|>?H &<A_J_`H.B^I@[e[H e[@UHWXBA@CI3PSR(AOSH.R1HiBUIR(HWX>kaeAklAI`PQu[PQIbqsA<AII`>AH r #t A_`@lu[Pe[HWX

IJ>AHeAkWI$k1R klBNkWDCH.R(HiB^IdPSB^_$kaXZIJ_`<NI`<A_JH0@C_0IJ>APXn_`H$NBAHiR(H.B[IdPQBk(OklIJHi_#A>kCXZH@l @[e[He^HuUH.OQ@aAR(H.B[I.M

~ BIJ>APXhXZHIJPS@CB9FWVH e^HWXJ&_JPQsNH klAAOSP.kWIJPS@CB@WEg, r0IJ@0knkWR(PQOSb#@[email protected] XpIJ>NkWIhPQB[uU@COSuUHPe^H.B[IJP.klO?XZ?>AH._JP.klO?uU@CPe?XH.RsHeAe^HWe0PQBklBPQB NBAPQI`HdsAPQBe^H._.M~\B3HWkC$>3A_J@asAOQHiRFCIJ>?HYXZA>AHi_JHWX

W

Page 18: A Fast Solution Method for Three-Dimensional Many-Particle

@C_JR k XPQR(AOSH(&<AsAP!kl_J_+k.bPQB Vd>AP > IJ>?HB[<ARsNH._0@lXA>AH._`HWXPXsNH.IVHiH.B -¨kWBNe - ! 5 ! klBNe I`>AHXNkCPSBADqsHiIZVH.H.B IV@qBAHWkW_JHWXI(XA>AH._`HWX1PX3H*)U<NklOdI`@ ^M _+kCe^PSPEPQD?M!$MY|>AHA_JPSBNPSNklO9@as[wJHWIJPSuUHPQBLIJ>APX#XHWI`PQ@CB%PXnI`@%e^H.R(@aBNXZI`_$kWIJH0IJ>NkWI#VH3.kWBX@COQuUH0I`>[email protected] X)U<AP$T^OSbklBNeLkC.&<A_$kWIJH.OSbUM

EPQDa<A_JH ! tXZPSR(AOQH4<As?PkW_J_$k.b(@W3XZ?>AH._JHXM

~\BL@C<A_0&@CR(A<?I$klI`PQ@aBNXFAVH<NXHWe IJ>AHsN@C<ABNe?kl_JbLH.OQHiR(H.B[I=R!HWXZ>LXZ>?@lVdBPSBEPSDC<A_`H5?M#~\IklOQOS@WVHWe<NX9IJ@4@CR!A<AIJH XZIJ_`HWXJXHWX9BAHWkl_IJ>?HX<A_ZkaH@WNIJ>?H=uU@aPeAXpI`@#VdPSIJ>APSB kC.&<A_$kabUM Hkl_J_`PQuUHWe0kWIhIJ>APXh@aBNOQ<XZPQ@aBs[b@CR(kl_JPSBADdIJ>?HYjdg= XZ@COS<AIJPS@CBNX9@C_IZV@uU@CPeAX9VdPQIJ> e^H.I+klPQOSHWeNBAPQI`H3H.OQHiR(H.B[IdXZ@aOQ<AI`PQ@aBNXM3|>?H3R(HWXZ>LPXn@a_JR(Hes[b BA@[e^HPX@CNkW_$kWR(H.IJ_`P0H.OQH.R!H.B^I+XkWBNeIJ>AHB^<ARsNH._p@We^H.DC_`H.HWXp@W_`H.HWe^@aR NH._,XA>AH._`H=PX=a ^M6PQIJ>?PQB0HWkC$>HiOQH.R(HiB^IiF PSB^I`H.DC_+klI`PQ@CBPX=NH._@C_`R(HWe3<NXZPSBAD 5 51]4kW<NXJXPklBvklBNe%XPQBADa<AOkW_dPQB[IJH.Da_$klI`PQ@aB%XJ$>AH.R(HX NH F C +M!Y<A_uUH._$XPQ@CB@lIJ>?HYEh @[e^HdV kCX,e^H.uUHiOQ@CNHWe0s^b3rUPSBADCHi_dJ4 Vd>A@AFUIJ@4PQR(A_`@WuUHhIJ>AH

Hi_Z@C_`R%klBH#@WIJ>AHYR(H.I`>A@[epFU_JHiA_JHWXH.B[IJHWe3IJ>AH4XZ?>AH._JP.klO>Nkl_`R(@CB?P.XVdPSIJ>(B?@C_JR klOQP&%WklI`PQ@aBNXe^P #H._`H.B^Ih_`@CR]#_`H.H.BADkl_+e JCCU$M H4<NXZHWeIJ>AH4Eh¥Hi^NkWBNXZPS@CB%e[H.AIJ> .¡3klBe1IJ>AHNkl_`IJPQI`PQ@aBAPQBAD e^H.AI`>VdkCX4 >[email protected]<N$> IJ>klII`>AH3B^<?R1sNH._d@Wd$>NkW_JDCHX#PQBI`>AH3XR%klOSOQHWXIYH.OQOXVdkCX3&OQ@UXH1I`@PX0OS@UXZH!IJ@ @aBAH%>^<?BNe^_JHepM |>?PX:B[<ARsNH._0VdkCXe^H.I`H._JR(PSBAHWe H.[NH._JPSR(H.B[I$kWOQObPQB%@C_+e^H._YIJ@sNklOklBNH0I`>AHklR(@a<AB^Id@W@aR(A<AI+klIJPS@CBNXPSB^uU@aOQuUHWe3PQBLIJ>AHX<AR%XPQBq !4ZrUPSBADCH._iF 4 §9]Y_JH.H.B?DUkl_+e kWBNe cd@CT[>AOQPSB9Fp4 $MEPQ_$XI.FVHqklBNkWOQb %.HWe6IJ>AHqDCH.@aR(H.I`_JPQHX1IJ>NkWIL.klB sNH XZ@aOQuUHWe6@CB k e^HXZT[IJ@C @aR(A<AI`H._

i

Page 19: A Fast Solution Method for Three-Dimensional Many-Particle

EPQDa<A_JH 5 ,|>AH3X<A_ZkaHe^PX`_JH.I`P1%WkWIJPS@CB(@W,kXZA>?H._JHUM

@CB?H1!4! %YA_J@[HXJXZ@a_#@Wr 2=H )<APSANHWe1VdPSIJ>q.a1 s%cYt +M#@aBAH@lIJ>?HA_J@asAOQHiR%X_JH.N@C_`IJHWe0>AH._`Hd.kWBsNHdX@COQuUHe:@CBIJ>APX&@CR(A<?IJH._9VdPSIJ>3X@COSuUH._$X9IJ>NkWI_`H*)U<APQ_`Hn!Z + pR!H.R(@C_`bIJ>?H1A_`@CsAOSH.R8IJ>klIPSB^uU@aOQuUHWXdIV@vXZA>AHi_JHWX@CB?OQbqPXuUH._JbvOQ@XZHIJ@LIJ>AHR(H.R(@a_JbL.kWNkCPSIbOQPQR!PQI.M~\B@aB^I`_$kaXZI.FNVHVH._`H3kWsAOQHIJ@!XZ@COSuUHA_`@CsAOSH.R%XVdPSIJ> <?IJ@ 5%XA>AH._`HWXklBe 5AFSCi<ABAT[BA@WVdBNXn<XZPQB?D!ZL,R(H.R(@a_JbPSIJH._+klI`PQuUH0X@COQ<?IJPQ@aBLXZIJ_+klI`H.DCPSHWXMm,@aR(A<AI+klI`PQ@CBklO[_JHWX<AOQI+X9kl_JH XZ<AR!R%kl_`P1%.HWePSB3EPSDC<A_`H Vd>AP >@aB^I+klPSBNX9k#OQ@aDlOQ@aD#AOS@CI9@C_

IJ>AHYm2= IJPSR(H,kCXkn<ABNI`PQ@aB@lvMh|>AH _J@XJXhXb^RsN@COX9&@C_J_`HWXZN@aBNeI`@#I`>AH#!ZL9PQI`H._$kWIJPSuUHXZ@COS<AIJPS@CB XZI`_$klI`H.DCbqIJ>klI0I$kWTUHWX0kCe^uaklB[I$kWDCH(@WYIJ>AHE klBe'I`>AH(&PQ_$&<AOkW_0XZb[R1sN@aOX#I`@klB(,+ XZIJ_+klI`H.DCbUM |>AHLOklI`IJH._XIJ_+klIJHiDCb VdkaXPSR(AOQHiR1H.B[IJHe s^b X<AAA_`HWXJXPQBAD >APQH._+kl_+ >AP.klONkl_`IJPQI`PQ@aBAPQBADkWBNe &@CBNXH*)<?H.B^I`OQbH.OQPSR(PQBklIJPSBAD0IJ>AHHiOQO H.OSONPQB[IJH._+kCI`PQ@CBX=VdPQI`>APQB%I`>AHEhM|>^<NXIJ>AHY(,+XZIJ_+klI`H.DCbPX9BA@CIhI`>AHR!@UXZIH PQHiB^I9@aBAHPSBPQI+X&OkCX`XhXZPSBNHPQIPQB[uU@COQuUHXpXZ@aR(H@lIJ>AHE@WuUH._J>AHkCepM #H.uUH._`IJ>AH.OSHWXJX,VH3HWXIJPSR%klI`H:I`>NklIiFpPQBLIJ>AH3?_JHWXH.B^I4@aR(A<AI+klIJPS@CBNX&FIJ>AH%NH._Z@a_JR%kWBNH(@W0IJ>?H'!Z+XZI`_$klI`H.DCb VdkCX3VdPQI`>APQB a¡ [email protected]@aAIJPSR%klOM|>AH_JH.OklIJPSuUH

+ H._`_J@C_Y@lIJ>?HE R%kWIJ_`PQU\uUHW&IJ@C_,R<AOQI`PQAOSPS.klI`PQ@CBX V kCX(`.¡ 0.$M

|>AHB[<ARsHi_d@l,NH._@C_JR!HWe(R%klI`_JPQUuUHWIJ@a_R<AOQI`PQAOSP&klI`PQ@CBX=VdkaX#H*)U<NklOIJ@ @a_kWOQO.kCXHWXM

Page 20: A Fast Solution Method for Three-Dimensional Many-Particle

g=ka >!B^<ARsNH._h@WR%klI`_JPQUuUHWIJ@a_R<AOSIJPQ?OQPiklIJP@aBNXh@a_J_JHXZN@CBNe?XIJ@3k0_JHWXPe^<klOpuUHIJ@a_hVd>A@UXH + B?@C_JR VdkaX'i¡ @l IJ>AH + B?@C_JR @WIJ>AH_`PQDC>[IZ>NkWBNe XPe^HuUHIJ@a_.M cnHWX<AOQI+X1@lI`>AHWXZH@CR(?<AI$kWIJPQ@aBNXX<ADCDaHWXZIIJ>klI3IJ>?H(sA_JHklTU\H.uUHiBN@CPQB[IPXA_J@asNkls?OQb OS@UXH(IJ@ NuUH!IJ>A@a<NXJkWBNeklBNe I`>NklI@a_ ( 5LIJ>AH%(,+4XZI`_$kWIJH.DabsH@CR(HX#PQR(?_$kC&IJPiklOXZPQBH3PQI0I$kWTUHWX0klsN@C<?IIJH.BLeAk.b?XIJ@@aR(AOQHiIJHYIJ>AHwJ@CsM,~ IYkWANHWkl_+XnI`>NklIdI`>AH3( XZI`_$klI`H.DCb!Hi_Z@C_`R%X=@aAIJPSR%klOSOQbsHiIZVH.H.B ( 0kWBNe ( 5kWBNe1I`>AHikCXZH4@l ( kWANHWkl_+X=IJ@sNHBA@CIOkW_JDCH4H.BA@C<?DC>9M~ BI`>APXY_JH.Dkl_$eF@aBAHXZ>?@C<AOe'_`HWklOSP1%.HIJ>NkWI:PSB @a_$e^H._0I`@%klI`I$klPSB I`>AH(!Z YNH._@C_JR klBNHFXZ>A@a<AOesNH3Okl_`DCH3sA<?IBA@CI0I`@U@Okl_`DCH3PQB@a_$e^[email protected]@CPe<NXPQBAD%IJ>?H3>Nkl_+ee[PXZTLR(H.R(@a_Jb@C_NklDaPQBAD[$FVd>AP > R%k.b XZPSDCBAP .klB[IJOSbe[H.IJH._`PQ@a_$klI`H%IJ>AHLHi_Z@C_`R%klBHUM|>?PX3PXJX<AH%sNHW@CR!HWX@Cs[u^PS@C<NXPSBI`>AH0BAH.[IYXZH.Id@W@aR(A<AI+klIJPS@CBNX&M

104

105

103

104

105

106

Number of unknowns

CP

U ti

me

(sec

onds

)

EPQDa<A_JH ,m,@CR(A<?I$klI`PQ@aBNklON_JHXZ<AOSI$X@C_^F[U0klBe 5 XZA>AHi_JHWXh@Cs?I$klPSBAHWe3VdPQI`>1@aBAH#?_J@[HWX`XZ@a_@lr 2,^Mm,_`@UXJXHWX,@C_`_JHWXN@CBNe3IJ@I`>AHZLhXIJ_+klIJHiDCb3kWBNe3PS_$OSHWX&@C_J_`HWXZN@aBNeIJ@0IJ>?HZ+XZIJ_+klI`H.DCbUM

|>?H NkW_$klOSOQH.OuUH._$XPQ@aB@WI`>AH&@[e^H%VdkCX3IJHXZIJHe @aB6kWB6~jd r 2, &@CR(A<?IJH._VdPQIJ> !4!! %A_`@[HWXJX@C_+X#HWkC$>%H*)U<APQANHWe%VdPQI`> iC( sLcYtMN|>AH0A_`@CsAOSH.R%XIJHXZIJHe @aB I`>APX

Page 21: A Fast Solution Method for Three-Dimensional Many-Particle

@CR(?<AIJH._YPQB[uU@COSuUHWeI`>AH!.kCXHWX#@C_ 5?F=i F,^ LklBNe ! 5 ! XA>AH._`HWX0klBNeqklOSO,@CR!A<AI$kWIJPS@CBNXVH._`HNH._Z@a_JR(Hey<XZPQB?DIJ>AHY( XIJ_+klIJHiDCb@aBAOQbM=m,@CR!A<AI$kWIJPS@CBNkWO^_`HWXZ<AOSI$Xkl_`HYXZ<AR(R kl_JP&%.HWePQB!EPQDa<A_JH :Vd>AP >&@CB[I$klPSBNXk0OQ@aDlOQ@CDAOQ@aIh@C_I`>AH0m2= IJPQR!HNH._,PQI`H._$kWIJPQ@aB(kaXk<ABNI`PQ@CB@lM #@aIJHI`>NklIdI`>AHNH._Z@a_JR%kWBNHYPX,R1<N$>3sNH.IJI`H._=@a_dIJ>AH4R(HWe^PS<AR1 XZP1%&HY.kaXZHWX&F ( . klBNe( [ [FAIJ>NkWB@C_ ( *5(kWBNe)( ! 5 !^MEA@C_ ( 5?FI`>AH0@U@a_#NH._@C_`R%klBN&HPXe^<AHI`@(IJ>?H:kC&I#I`>NklIYIJ>AH @CR!R1<AB?P.kWIJPQ@aB1I`PQR(Hy@CBNXIJPSIJ<AI`HWX#k!XZPQDaBAP .klB[I#N@a_JIJPS@CB@W,IJ>AHIJ@CI+klO9I`PQR(HUFaP\M«HMI`>AHA_`@CsAOSH.R/PX=B?@CIdOkW_JDaH:HiBA@C<ADa>9M,EA@C_ ( ! 5 !^FAI`>AHN@U@C_YNH._Z@a_JR%kWBNHPX0e^<AH(I`@ I`>AHkaII`>NklIIJ>AHA_`@[HWXJX@C_3I`>NklI e^PXIJ_JPSsA<AI`HWXIJ>AH!eAklI+k>NkaeIJ@ _`HWXZ@a_JI0IJ@LIJ>AH>Nkl_+ee[PXZTqR(H.R!@C_JbFPM«HUMIJ>AH(?_J@Cs?OQH.R8PXyI`@U@OkW_JDCH3@C_IJ>AHv ?_J@[HWX`XZ@a_$XM H!e^PeqBA@CI@CR(?<AIJHY@aR(AOSH.IJH XZ@COS<AIJPS@CBNX@a_IJ>?H@aBNXZPe^H._`HWeA_J@asAOQHiR%X@CB?OQbe^<AHYIJ@0OSPQR(PSIJHWekC.HXJXhIJ@klOQO, ?_J@[HWX`XZ@a_$Xd@l,I`>AH@CR(?<AIJH._iM

105

106

102

103

104

Number of unknowns

CP

U ti

me

per

itera

tion

(sec

onds

)

EPQDa<A_JH m,@CR!A<AI$kWIJPS@CBNkWO#_`HWXZ<AOSI$X@a_(IJ>?H !XZI`_$kWIJH.Dab@C_ *5AFi F^ klBNe ! 5 !XZA>AHi_JHWX@asAI$kWPQBAHWe%VdPQI`> ?_J@[HWX`XZ@a_$Xd@ldr 2,[M

[

Page 22: A Fast Solution Method for Three-Dimensional Many-Particle

* *

~ B I`>APXNkWNH._.FVH%PQB[IJ_J@[e[<NHWe Eg, rNF,kWB !Z R(H.I`>A@[e'@a_(kWBNklOSbAXPX3@l4Okl_`DCH XJiklOQHIJ>A_`H.H e^PQR(HiBNXZPS@CBNkWOdR(P_`@CR(HJ>NklBAP.klOYkWBNklOSbAXPX1@W@aR(N@UXPQIJHR klIJHi_JPkWOX&M|>APX(R(H.I`>A@[ePX0XZ<APSI$kWsAOQH@a_kWBNklOSb %iPQBAD(?_J@Cs?OQH.R%XYIJ>NkWI:PSB^uU@aOQuUH0>[<ABNe^_`HWeAX0klBNeqIJ>A@a<NXJkWBNeAX:@WdNkW_JI`POSHWXVdPQIJ>?@C<AIPSR(N@UXPQBADkWB^bLHWXJXH.B^I`PkWO9_JHWXIJ_JPIJPS@CBNXd@CBLIJ>AHiPQ_YDCH.@aR(H.IJ_`bUMd|>AHA_JPSBNPSNklOPe^Hk@WEg, r PXyI`@L@CRsAPQB?H3k%@CB[uUH.B[IJPS@CBNkWO,jdg, VdPSIJ> klBqPQI`H._$kWIJPQuUH XZ@aOQ<AI`PQ@CB XZI`_$kWIJH.Dab I`>NklIPQBNOS<Ne^HWXIJ>AH4EhM4qH#PSR(AOSH.R(H&B^I`HWeIJ>APX Pe^HWkPQB!k:V k.b3IJ>NkWIkWOQOS@lVYX,<NXhIJ@ @CRsAPQB?HdklB[bjdg= VdPSIJ>!klB[bEhklBNe(e[H.uUH.OQ@aNHWe3kNkl_+klOQOSH.ON@[e[H#I`>NklI .klBsNH#<A?DC_$kae^[email protected]`IJH._jdg=8@C_YE@a_JR<AOkWIJPQ@aBNX=sNHW&@CR(H0k.uaklPSOkWsAOQHUM|9@ k H._JI+klPSB6H.^I`H.B[I.FI`>APX1klNH._L.klB sNH u[PQHiVHWe kCX%k N_+XZI%_JH.N@a_JI IJ>AHV@C_JT PX

PQBA_`@CDa_JHWX`X1kWBNe XZH.uUH._+klO,[email protected][I$X#I`@ IJ>?H &@[e^HLkl_JH <A_J_`H.B[IJOQb<?BNe^H._3V kbM|>A@UXHPQBNOS<Ne^H%PQB@C_`@a_$kWIJPQ@aB @WIJ>AHLPQR(A_`@WuUHWeEh ?_J@CN@XZHWe6s^b6]Y_JHiH.BADUkW_$e6kWBNe cd@CT[>AOSPQBJ "U$FR(@[e^HiOQPQB?D@l3R!@C_JHL@aR(AOSH.DCHi@CR(H.I`_JPSHWXF=<?_JIJ>?H._1kl_$kWOQOSH.OQP&%WklI`PQ@CB @lIJ>?Hq@[e^HUFe^H.uUH.OS@CAR(HiB^I@WOS@[.klOdR!HWXZ>_`H$NBAHiR(H.B[IIJHW$>ABAP )<?HWXFklBNe A_JHW&@CBNe^PSIJPS@CBAPSBADAM%!n4&@C<A_+XZHUFIJ>AH(R!H.IJ>A@[eqBAH.HWeAX0IJ@sNH(I`HWXZI`HWe @CB R(@a_JH!@CR!AOQH.LDCH.@aR(H.I`_JPQHXyI`>NklB I`>A@UXH &@CBNXPe^H._`HWePQBIJ>?PX#NkWNH._.M YH.uUH._JI`>AH.OSHWXJX&FAVH3sNH.OQPSH.uUH0IJ>NkWI0IJ>AH<ABeAklR(HiB^I+klOXd@l Eg, rvkl_`HXZ@C<?BNeklBNe1PSI,VdPQOQONsNHW@aR(HYk:<?BAP )U<AH4klBNe1<XZH<AOIJ@U@aOp@lR!P_`@CR(Hi >klBAP.klO[kWBNklOSbAXPX=@W@CR!@XZPSIJHR%klI`H._JPklOX&MEPQBNkWOQOQbFhVH%BA@CI`H%IJ>NkWI1I`>AH%DCH.BAHi_$klOSP1%.He R%klB[bU\sN@[e^bHiOQHWI`_J@XZI$kWIJP.X:A_`@CsAOSH.R.klB sH

H.^I`H.BNe^He1IJ@IJ>AH0TUH._`BAH.OX=IJ>klI4 >NkW_$kaIJH._`P1%.HYOQ@aBADl_$kWBADCHe[PXZOS@[.klI`PQ@aB PQB[IJHi_$kC&IJPQ@aBNX#cn@[e[PQB9F 4UC+F9XZ@(I`>NklIYIJ>?HPQB[IJHi_$kC&IJPQ@aBNX@l e^PXOQ@[[email protected]!H.B^I+X.kWB sNH@aR(A<AI`HWe1PQBq( _$klI`>AH._@C_$e[PQBNkW_JPQOSb(( + @CNH._+klI`PQ@CBXM !#_`PQDaPQBNkWOQOQbFCIJ>AH4EhVdkCXklA?OQPQHe1IJ@IJ>APX=&OkCX`X=@[email protected] X=s[bkWBAD!klBNe NHCr[kl_J4 $M

|>APXhV@a_JT3V kCXhDCH.B?H._J@a<NXZOSbXZ<AAN@a_JIJHe1s^b3I`>AH kWIJPQ@aBNklO9r[&PQH.BN&HYE@a<ABNeAkWIJPS@CBIJ>A_`@C<ADa>(IJ>AH

a

Page 23: A Fast Solution Method for Three-Dimensional Many-Particle

klI`PQ@CBklOh]Y_$kWBNe'm,>NkWOQOSH.BADCHXA_J@aDC_$kWRDC_$kWB^I4gmdrC 5Ua¡[M

3rNM[t#OQ<A_`<9F]#_`H.H.BADkl_+e X1sN@[e^b!klOQDa@C_`PQIJ>?R PXBA@aI@a_$e^H._ QFAr[~tKAM^rP\M^m,@aR(A<AIiM Fpa*!`C J "U$M

3mYMAtM?j_JH.sAs?Pk0klBNe%KAM #@CR(PSBADC<?H %UF $M

lBNe%HWe^B9M9m,@CR(?<AI$kWIJPQ@aBNklO H >NkWBAPiX2h<AsAOSP.kWIJPQ@aBNXFNW]#_+k.V #PSOQOj@U@CTm,@CR(klB[bUF#H.V ;=@a_JTFpC4[M

! MnmYMg,e^VdkW_$eAX&F t kl_$kWOQOSH.O#PSB^_$kaXZI`_J<NI`<A_JH%@C_%X`.klOklsAOSHkaeAklAI`PQuUH,NBAPQI`H H.OSH.R(H.B[IR(H.IJ>?@[eAXkWBNe(PQI$XdkWAAOSP.kWIJPQ@aB(IJ@3OSHWkCXIdX )<NkW_JHWX @aOQOQ@[iklIJPS@CB SF 2h> #PXJXZHi_JI$kWIJPS@CB9Fm,@CR(A<?I$klI`PQ@aBNklO4klBNe t#AAOSPQHWe klIJ>?H.R%klI`P.X&F,|>AHYBAPQuUH._+XZPSIZb@W|9H.AkaX3klI!t#<NXIJPQBJ "a$M

5 MmYMg,e^VdkW_$eAX0kWBNeKAMmYMj_J@WVdBAHF |>AH3X`.klOklsAOSH3e^PXZIJ_`PQsA<?IJHWeLe^b^BklR(P:kW_J_+kbLkWBNePQI$XYkl?AOQP.klI`PQ@CBLIJ@%kNkl_+klOQOSH.O9>A^ kCe?klAI`PQuUHNBAPSIJHH.OQH.R!H.B^I@[e[HQF !

#" $

%'& )( * '$

,+ .- $ NFAr[klB[I$k3EAHUF FN "[M

tMh]#_`H.H.B[sNkl<?RF M=]Y_`H.H.BADkl_$e klBNetMk.bU@?F !YB I`>AH(B[<AR(H._`P.kWOhX@COQ<?IJPQ@aB @WYIJ>AHsAPQ>NkW_JR(@aBAPYH*)<klIJPS@CB%PQB%I`>AH0AOkWBAHQF/0$1

3254 9FN^ .a J4a$M

M#]Y_JH.H.B?DUkl_+epF768$ 9 : <; = #"7 > @? =A 7BC

>

~Z| 2h_JHWX`XFpmkWR1sA_`Pe^DaHUFAtMQFN aC^M

M]Y_JHiH.BADUkW_$epF ENkaXZIklOSDC@a_JPQI`>AR%XA@C_OkaXJXZP.klOaA>[bAXP.X QF!BC > ED14GFF"C¡" ^$5` 4*5^+M

Mn]Y_JHiH.BADUkW_$e6kWBNez3M,cd@CT[>AOSPQB9F t kCXZI!klOSDC@a_JPQI`>AR%X:@a_(Nkl_`IJPOQHLXZPSR1<?OklI`PQ@aBNX SFIH C J$1

1K F !a !*5UJCa$M

"!

Page 24: A Fast Solution Method for Three-Dimensional Many-Particle

M]Y_JH.H.B?DUkl_+e kWBNez3Mpcd@CT[>AOSPQB9F t7BAHiV uUHi_$XZPS@CBL@lIJ>AHkaXZI0R<AOQI`PQN@COSHR(H.I`>A@[e@C_IJ>AH klAOkCHLH*)<klIJPS@CB PQB I`>A_JHiH e[PQR(H.BXZPQ@aBNX SF #H.NkW_JI`R(H.B[I3@l3m,@aR(A<AI`H._(r[PSH.BNHUF;kWOQH#BAPSuUH._$XPQIbUFUcdHWXHWkl_+ >%cdH.N@C_`I#;t g, :mdr c#c$aC J4U$M

.¡ M 4kCXZ>?PQB9F t#BNkWOQb?XZPX @W@CR!@XZPSIJHYR%klI`H._JPklOX QFH >$ F 9F 5U[ ¡ ` C"!+M

C MjM NP\F,]MMkWBkWBNe MtMklB?DAF t7BAH.V R(HiIJ>A@[e@a_H.uaklOS<NklI`PQBADvXZPSBADC<AOkl_PQB[IJH.Da_$klOXPQB!XZI`_JHWX`X,klBNkWOQb?XZPX@lX@COSPeAXhs[bIJ>?HYe^PQ_`HWIhsN@C<?BNeAkl_`b1HiOQH.R(HiB^IR(H.IJ>?@[e QF > H '$

! ! D F?C¡UU C¡"CqJ a $M

. GM kWsN@C_$X3kWBNeEM|MGy@a_$XR(H.bUH._0kWBNeEM|MNH.PSDC>[IJ@CB klBNe KAM 6>APQI`HUF 2h_JHW@aBNe^PIJPQ@aBAHWepFkaeAklAI`PQuUHUF9R<AOQI`PQN@COSH kC.HiO Hi_$klI`HWe(PQI`H._$kWIJPQuUHR(H.I`>A@[eAXn@a_:I`>A_JHiHe[PQR(H.BXZPQ@aBNklON_$XIZT^PSBNePQB[IJH.Da_$klO,H*)U<NkWIJPQ@aBNX:@[email protected][IJPklOI`>AH.@a_Jb QFJB1 H 0BC 0B C FFp$5l`*![%J " 5[$M

! M 2hkW_$kaXZ>NkW_0klBNeKAMpmYMj_J@WVdBAHF rUb?XZI`H.R H.BADaPQBAHiH._JPSBAD@a_#>?PQDC>LNH._Z@a_JR%kWBNH@CRA<AIJPSBAD0X@lIVdkW_JH IJ>AH #t #t] PQB^_+kCXZI`_J<N&IJ<A_`H@C_hPSR(AOQHiR1H.BUI+klI`PQ@CB0@Wpkl_$kWOQOSH.OXZIJ_`<NI`<A_JHWe kCeAkWAIJPSuUH%R(HWXZ> _JH$NB?H.R(H.B[I SFI`@kWANHWkl_!PQB B ;

$

9 $ F~t7uU@COS<AR(HWXYPQB klIJ>?H.R%klI`P.X4klBNe PSI$X0t#A?OQPiklIJPS@CBNX&F

rUA_JPSBADCHi_dz=Hi_JOkWDAF #HiV6;=@a_JTF9 "UM

$5tM2M 2H.PQ_+HkWBNe KAMdtM M kWAPQH._iF t XZNHW&IJ_$kWO0R1<AOSIJPS@aOQH%R(H.I`>A@[e@C_%H PSH.B^IXZ@COS<AIJPS@CB%@lOkl_`DCH XJ.kWOQHsN@a<ABNeAkW_Jb(H.OSH.R(H.B[IR(@[e^H.OX PSB(H.OkCXZI`@UXI$klI`P.X QF H $

! ! K F 5U¡a¡4 5U¡"! 5qJ4 $M

M2h>NkWB^|>?PQH.B klBe rM,GyPSRF$

:

1

C $ M !#U@C_+e%#BAPQuUHi_$XZPSIb 2h_JHWX`XF #H.V6;,@C_`TpF94*5AM

3]MKAMcd@[e^PQB9F |9@WVdkl_+e'_+klAPe'HiuCkWOQ<NkWIJPS@CBq@lYIJ>?H1H.OkCXIJPPSB^I`H._$kaIJPS@CBNX4klR(@aBAD%IJ>A_`H.He^PQR(HiBNXZPS@CBNkWO9e^PXZOQ@[iklIJPS@CBNX QFNX<AsAR(PSIJIJHeI`@ 0$

` "C+M

*5

Page 25: A Fast Solution Method for Three-Dimensional Many-Particle

l 0z3Mcn@aT^>AOSPQB9F cYklAPe X@COQ<?IJPQ@aB @WYPQB[IJH.Da_$kWOhH*)U<NklI`PQ@[email protected]@CIJHiB^I`PklO,I`>AH.@C_`b SFH 0 C 80$1

14 9F9i C¡ULJ a $M

. 0;3Mr[kakCe%klBeM MNr[$>^<AOSI %UF ]4c#gr ?t DaH.BAH._+klOSP1%.HWeR(PQBAPSR(kWO^_`HWXZPe^<NkWO9klOSDC@a_JPQI`>AR@C_XZ@COSu^PSBADYBA@CBNXb^R!R(H.IJ_`PpOSPQBAHWkW_XZb?XZI`H.R%X SF B CH BC B G C Fa 4J aU$M

Mpr[klDklB9FB ? ! 0 ; MprU?_JPQB?DCH._dz,H._`OklD?F#HiV ;,@C_`TpF " 5AM

C¡ tMrMr[kWBADUkWBAP,kWBNe ]M9 @AF t#B Z nklOQDa@C_`PQIJ>?R @a_rUI`@CTUHWXYklBNe kWAOkaH3PQB[IJHi_$kC$IJPQ@aBNXd@lNkl_`IJP&OQHWX 0$1

8? =

F "C¡ C¡[.¡qJ "U$M

^MrNM kl_J_`H.BqklBeKAMpG3M9r[klOSR(@CB9F tXZI`_J@aA>^b?XZP.klO d\sN@[e^bqXPQR<AOkWIJPQ@aBNX<NXPQBAD>APQHi_Zkl_$$>AP.klO,IJ_`H.H!eAklI+kXZIJ_`<NI`<A_JHWX QF~ B BC . C ! C[Fh PSBABAHWkW@aOQPXFh F~g,g,gm,@CR(A<?I.Mr@[CM 2h_`HWXJX&F W¡. J4a$M

C MrNM kl_`_JH.B5klBNe KAM0G3Mr[klOSR(@CB9F tNkl_+klOSOQH.O3>NkaXZ>AHe7@[IZIJ_`H.H n\sN@[e^b5klOSDC@l_JPQI`>AR SF#~\B BC C ! 4![FkaXZ>APSBADCI`@CB9F 0mYF#~g,g,g m,@CR(A<?I.MrU@[CM 2h_JHWX`XF. ^!J "4!U$M

4! KAMUG3MrUPSBADCHi_.F 2hkl_+klOQOSH.O[PQR(AOSH.RH.B^I`kWIJPQ@aB0@lIJ>AHkaXZI,R<AOQI`PQN@COSHhR(H.I`>A@[e:VdPSIJ>3NH._JPS@[e^Ps@a<ABNeAkW_Jbv@CBNe[PQIJPS@CBNX QF/

)(

IH $ K F94 ^ J4 $M

50;3M4M kWBADklBNe3c0M NHar[kl_iF ZLklOSDC@a_JPQI`>AR @a_he[PXZOS@[.klI`PQ@aBe^b^BklR(P.X SF!J$

C F9$5 + "! J4 $M

Page 26: A Fast Solution Method for Three-Dimensional Many-Particle

EPQDa<A_JHCM,t sAPSBNe^H._Y&@CB[I$klPSBAPQB?D3Pe^H.B[IJP.klONkl_`IJPOQHWX&M

EPQDa<A_JH^M2kl_JI`PQ@aBAPQBAD@lIJ>?H&@CR(A<?I$klI`PQ@aBNklOe^@aR%klPSB(sNkaXZHWe%@CBk3XNkC&H OQOQPSBAD<A_`uUHUM

EPQDa<A_JH !^M,tXPQR(AOSH&<AsAP0kW_J_$k.b(@W,XZ?>AH._JHXM

EPQDa<A_JH 5AM,|>?H NBAPSIJHHiOQH.R(HiB^I,e[PXJ&_JH.I`P1%WkWIJPQ@aB%@lIJ>AH0sN@C<ABeAkl_`b @l kXZA>AHi_JHUM

EPQDa<A_JH Mm,@CR!A<AI$kWIJPS@CBNkWO^_`HWXZ<AOSI$X@a_,^F[UUFUkWBNe *50XZ?>AH._JHX@CsAI+klPSBAHWeVdPQI`>3@CBAHd?_J@[HWX`XZ@a_@l9r 2,[M=m,PQ_+OQHX9@C_`_JHWX@aBNe:I`@#I`>AHd( XZI`_$klI`H.DCbkWBNeX )<NkW_JHWX@a_J_JHXZN@CBNeI`@#I`>AHd(,+XZIJ_+klI`H.DCbUM

EPQDa<A_JH ^M6m,@CR!A<AI$kWIJPS@CBNkWOd_JHWX<AOQI+X @a_3IJ>AHq( XZI`_$kWIJH.Dab @a_ *5AFi Fd^ ^FYklBNe ! 5 !

XZA>AHi_JHWX@asAI$kWPQBAHWe%VdPQI`> ?_J@[HWX`XZ@a_$Xd@ldr 2,[M