arxiv: v1 [physics.atom-ph] 9 Jun 2014

Transcription

1 epl drft Finite Size Dependent Dispersion Potentils Between Atoms nd Ions Dissolved in Wter rxiv: v [physis.tom-ph] 9 Jun 204 P. Thiym, M. Boström,2,3, C. Persson,2,4, D. F. Prsons 5, I. Brevik 3 nd Bo E. Sernelius 6 Dept of Mterils Siene nd Engineering, Royl Institute of Tehnology, SE Stokholm, Sweden, EU 2 Centre for Mterils Siene nd Nnotehnology, University of Oslo, P. O. Box 048 Blindern, NO-036 Oslo, Norwy 3 Deprtment of Energy nd Proess Engineering, Norwegin University of Siene nd Tehnology, NO-749 Trondheim, Norwy 4 Deprtment of Physis, University of Oslo, P. O. Box 048 Blindern, NO-036 Oslo, Norwy 5 Deprtment of Applied Mthemtis, Austrlin Ntionl University, Cnberr, Austrli 6 Division of Theory nd Modeling, Deprtment of Physis, Chemistry nd Biology, Linköping University, SE Linköping, Sweden, EU PACS Cf Intertomi potentils nd fores PACS k Theory of quntized fields PACS L Quntum flututions, quntum noise, nd quntum jumps Abstrt A non-expnded theory is used for dispersion potentils between toms nd ions dissolved in medium. The first order dispersion intertion between two toms in n exited stte must ount for the ft tht the two toms re oupled vi the eletromgneti field nd must inlude effets from bkground medi, retrdtion nd finite size. We show tht finite size orretions when two prtiles re lose hnge the dispersion intertions in wter by severl orders of mgnitude. We onsider s four illustrtive exmples helium toms, krypton toms, phosphte ions, nd iodide ions. We demonstrte tht due to lrge nelltion effets retrdtion domintes the intertion for helium tom pirs in n isotropi exited stte down to the very smll tom-tom seprtions where finite size orretions re lso importnt. We im in the present letter t lulting the Csimir- Polder intertion energy between toms helium nd krypton) nd ions phosphte nd iodine) dissolved in wter nd to study wht effets the finite tomi size hs on this intertion. We lso give results for the resonne intertion energy between He toms. The lultions presented re for very smll intertomi distneslower thn ten Ångströms) t whih effets beyond dipole-dipole dispersion intertions re importnt. The intertion energy between ground stte helium nd krypton toms, in vuum, re well known from stndrd b initio lultions. They show potentil energy minimum min =2.97Å, U min =0.00 ev for He 2, nd min =4 Å, U min =0.06 ev for Kr 2, where is the tom-tom seprtion nd U is the intertion potentil). The minimum is the result of repulsive short-rnge eletro-stti prt nd n ttrtive long-rnge prt, due eletron orreltion. We re here onerned with the ttrtive long rnge prt of the potentil only. The formlism we use fouses on the dipole-dipole dispersion intertion. At very smll seprtion higher order multipolr ontributions ontribute nd when there is substntil wve-funtion overlp one hs to resort to quntum-hemistry lultions sine the eletrons re no longer onfined to their respetive prtile. In our formlism we find tht the intertion stys finite when the prtiles ome lose to eh other. We emphsize tht the reson is tht the prtile rdii hve been ssumed to be finite. There re reent developments in the literture relted to our pproh. Thus Kysylyhyn et l. [] hve studied the intertion between finite-size nnoprtile ould be n tom) nd the surfe of solid, mking use of the lol-field method. Assuming tht the nnoprtile hs finite nonliner polrizbility, n intertion potentil is derivedtht is repulsivet short rngesnd hs n ttrtive long-rnge til. Moreover, Przybytek et l. [2, 3] hve done high preision work on the retrded Csimir-Polder potentilbetweentwohetoms. When > 00Åretrdp-

2 P. Thiym et l. tion mkes the ditomi potentil swith from the London 6 dey to the Csimir-Polder 7 form. The retrdtion effet is found to inrese the size of the helium dimer bound stte by bout 5%. The nonretrded lultions yieldnvergevlueoffortheheliumboundsttetobe quite lrge, 47.Å, whih grees firly well with the mesured retrded) vlue of 52Å. It my lso be mentioned tht n extensive tretment of quntum eletrodynmil QED) intertions nd thermophysil properties of helium hve reently been given by Cenek et l. [4]. A reent pper by DiStsio et l. [5] deling with density funtionl quntum mehnil) theory of the intertion potentil between finite-sized quntum hrmoni osilltors expliitly tkes orbitl overlp into ount. They find results tht re similr to ours in wter. In both pprohes finite tomi size renders the intertion finite t zero tom-tom seprtion. Wht is evident in our theory is the removl of the divergenes t zero seprtion. In this letter we tret both the Csimir-Polder intertion nd the resonne intertion. We begin by giving brief introdution to resonne intertion. Ninhm nd o-workers [6] demonstrted dede go tht, due to some drsti pproximtions, the underlying theory of resonne intertions in free spe derived from perturbtive QED is inorret t lrge distnes. [6 9] For smll seprtions the QED result hs been onsidered to be orret. However, the QED formlism supposes point dipolr desription. One im of this letter is to demonstrte how the tking into ount of finite tomi size nd bkground medium fundmentlly lters the resonne intertion energy. In the formlism one tom is in its ground stte nd the other in n exited stte. At resonne the exittion swithes bk nd forth between the two toms. The intertion n be seprted into three brnhes, the x-, y-, nd z-brnhwhere the nme of the brnh denotes in wht diretion the osillting dipoles re pointing. We let the z-diretion be defined by the line joining the two toms. When there is no preferred diretion for the exittion ll three brnhes re tivted. We refer to this se s isotropilly exited tom pirs. From writing up the equtions of motion for the exited system it is strightforwrd to derive the zero temperture Green s funtion for two identil nd isotropi) toms dissolved in wter [6, 0]. The resonne ondition [6] n be obtined from the following ondition: α ω) 2 T2 ω) = 0, where T is the suseptibility tensor nd α ω) the exess polrisbility of the tom dissolved in wter. The exess polrisbilities α i)) nd tomi rdii ) for helium nd other noble gs toms dissolved in wter were derived s in severl ppers by Prsons nd Ninhm [, 2]. They were presented reently by Boström, Prsons, nd o-workers. [3] Dynmi polrisbilities of noble gs toms in vuum were lulted using Molpro [4] t oupled luster singles nd double CCSD) level of theory. An ug--pv6z bsis set [5, 6] ws used for seletion of noble gs toms He, Ne nd Ar while ug--pv5z ws used for Kr [7]). In this work we onsider helium nd krypton toms nd iodide nd phosphte ions. Polrisbilities, αi), in vuum were trnsformed to exess polrisbilities, α i), in wter vi the reltion for dieletri sphere embedded in dieletri medium [8,9], α i) = R 3 ε ε w )/ε +2ε w ). Here ε w is the dieletri funtion of wter nd R is the rdius of the tom. ε is the effetive dieletri funtion of the tomi sphere, estimted from the tomi polrisbility in vuum s ε i) = +4παi)/V, where V is the volume of the tomi sphere. In the se of two identil toms the bove resonne ondition n be seprted in one nti-symmetri nd one symmetri prt. Sine the exited symmetri stte hs muh shorter life time thn the nti-symmetri stte the system n be trpped in n exited ntisymmetri stte. The resonne intertion energy of this ntisymmetri stte n be evluted by simple expression for two dissolved toms exited in the j=x-, y-, or z-brnh, where z is in the diretion of the line onneting the two toms) in wter, U j ) = h d 2π ln[+α i)t jj i)]. ) The orresponding vn der Wls Csimir-Polder) intertion between isotropi) toms is given by the following expression: U CP ) = h 2 j=x,y,z d [ 2π ln α i) 2 T jj i) 2], 2) where is the distne between the two toms. Trditionlly one ssumes tht the intertion is so wek tht one my expnd the logrithm in Eqs.) nd 2) nd keep the lowest order ontribution only ln + x) x). In wht follows we nme this the expnded theory nd when the full logrithmi expression is used the non-expnded theory. We now present the theoretil frmework used to study how finite tomi size influenes the Csimir-Polder intertion nd the resonne intertion when two toms re nerehother. Thepolristionloudofreltomshs finite spred nd we onsider s n interesting se onsistent with the modeling of the exess polrisbility) sptil distribution of the tom following Gussin funtionforhelium inwterthe rdiusis0.60å). [20,2]We hve, following the formlism developed by Mhnty nd Ninhm [20, 2], derived the Green s funtion elements tht ount for retrdtion, bkground medi, nd fip-2

3 Finite Size Dependent Dispersion Potentils 00 0 He-He in wter U ev) CP Full CP Expnded Dipole-Qudrupole Qudrupole-Qudrupole 0-6 Csimir Polder Csimir Polder Finite tom size !Å) Fig. : Color online) The vn der Wls intertion between pirs of finite size helium toms in vuum. Note tht finite size effets keep the intertion finite when the prtiles ome lose together. We lso show the vn der Wls intertion between pirs of finite size helium toms when the theory is perturbtively expnded. We lso show the urves for dipolequdrupole nd qudrupole-qudrupole ontributions in the expnded theory of dispersion energies. [22] nite size in n urte wy. They re nd T xx i) [ = T yy i) = exp ) 2 2) 2] {[ ) 2 ) 2 + [ )] ) + erf 2 exp [ ) 2 ) ) ] 2 + [ )] ) erf exp 4 π exp [ 4 π 2 + ) 2 2) 2 ) 2 ]}, T zz i) [ = exp ) 2 2) 2] {[ ) ) ] 2 + [ )] ) erf 2 exp [ ) ] 2 [ + ) erf + ) ) [ 2 exp )] 2 + ) 2 2) 2 ) 2 ] ) exp ]} ) 2, 3) 4) where = 0 / εi) with 0 the veloity of light in vuum nd εi) the dieletri funtion of wter for imginry frequenies. The effet of the bkground medium is entirely ontined in the veloity of light in the medium nd the exess polrisbility. In some ppers ftor εi) multiplies the exess polrisbility nd divides the Green s funtion elements. The results would of ourse be the sme if we dpted tht nottion. When retrdtion is negleted we find the following re- 0-8 Fig. 2: Color online) Csimir-Polder ground-stte) intertion between two helium toms dissolved in wter. Compring theory with or without finite size dependene. sults: nd T NR T NR xx zz i) [ = 2 π 3 = π 3 i) = TNR [ yy i) πerf ) 2 ) e ) 2], πerf ) +2 ) [ + ) 2 ]e ] ) 2, 5) 6) ] Tre[ TNR i) = 4 / 2. 7) π 3 e2 We then find with either the x-brnh or z-brnh exited / 3 -dependene of the resonne intertion in the non-retrded limit. We observe tht the resonne intertion t lose ontt depends on the rdius 3. This should be ompred to the orresponding result for the ground stte Csimir-Polder intertion tht depends on the rdius 6. The nlytil symptotes re only vlid for tomi systems where the logrithm n be series expnded. It is importnt to note tht we re only onsidering the dispersive dipole intertion. We hve not inluded ontributions due to higher-order multipole ontributions. [22 24] We show the resulting dipole-dipole, dipole-qudrupole, nd qudrupole-qudrupole dispersion potentils [22] in n expnded theory for two helium toms in vuum in Fig. note tht the dipole-dipole intertionis ftor of2too smll in Ref. [22] nd the preftor in their Eq.34) should be 6/9π insted of 8/4π). It is ler tht dipole-dipole ontributions give n importnt ontribution t ll seprtions. As omprison we lso show the non-expnded result for the dipole-dipole dispersion potentil. Here we see tht non-expnded theory my ount prtilly for the observed short-rnge repulsion between two toms in vuum. The ft tht the toms re in wter will s we will see hnge this ft. In simultions tht use Lennrd-Jones potentil p-3

4 P. Thiym et l. He-He in wter -U retrded -U nonretrded -U retrded, finite size -U nonretrded, finite size U ev) 0 He-He in wter U zz U zz finite size -U xx -U xx finite size Fig. 3: Color online) Resonne intertion between isotropilly exited ditomi helium dissolved in wter. The effets of retrdtion nd finite size = 0.6Å) re onsidered. For distnes smller thn the tomi rdii retrdtion effets n be negleted. Fig. 4: Color online) Resonne intertion between x- or z- mode in n exited stte, with nd without finite size effets inluded. In this plot retrded nd non-retrded results re indistinguishble. U xx isttrtive, nditsmgnitudeisshown here for omprison with U zz. the intertion is mde up of n ttrtive dispersion ontribution nd repulsive ontribution. We here lulte the ttrtive prt in orret wy. A full quntum mehnil moleulr simultion ould be used to study how well our preditions represent the effet of wter moleules on the intertion energy between two toms t smll intertomi distnes. However, it should be observed tht while quntum hemistry lultion is n lterntive wy to model the intertion tht too hs its limittions. Quntum hemistry does not, for exmple, in n urte wy ount for the oupling of toms in n exited stte. We show in Fig.2 the fully retrded finite temperture Csimir-Polderground-stte) intertion between two helium toms dissolved in wter. Finite size orretions prevent the ttrtion from going to infinity when two toms ome very lose together. This is onsistent with the ground stte Csimir-Polder intertion between finite sizetomsinvuum found bymhntynd Ninhm. [2] In Fig.3 we onsider the effet of inluding retrdtion nd finite size on the resonne intertion between isotropilly exited ditomi helium dissolved in wter. Finite size effets soften the ttrtion t very smll seprtions but eqully importnt is the effet of retrdtion. One observes tht for isotropilly exited ditomi helium the non-retrded pproximtion whih neglets finite size is not pplible t ny seprtion. Both finite size dependene nd retrdtion must be onsidered in order to hve theory tht works for ll seprtions. For isotropilly exited tom pirs the first expnsion term in the logrithm nels out nd the leding non-retrded term is 6. The dominting term is retrded term 4. We ompre in Fig.4 the effet of finite size on the ttrtive x-mode exittion nd repulsive z-mode exittion. It is seen tht while finite size effets re oneptully importnt t very smll tom-tom seprtions these effets ply little role for seprtions beyond two tomi rdii. Retrdtion effets re negligible for eh of the ontributing terms. However, there is lrge nelltion effet between the different U j ) terms) when the ontributions re dded together nd retrdtion is importnt for the resonne intertion between isotropilly exited tom pirs s ws seen in Fig.3. One importnt point in this letter is tht while retrdtion redues the x, y, nd z exited stte resonne intertions / 3 / 4 ) the opposite is found for resonne intertion between isotropilly exited tom pirs / 6 / 4 ). Severe problems hve been noted with the epted theoretil expressions for the resonne intertion between identil tom pirs in n exited stte when the toms re fr prt [6, 7, 9]. Different reserh groups find very different results for the first order dispersion potentil for tom pirs in n exited stte t tht seprtion limit. They do however find the sme result when two toms or moleules re lose together. There re no debte onerning ground stte intertions, nor for nisotropi) resonne intertions in the non-retrded limit. However, our results for isotropilly exited tom pirs propose tht retrdtion due to strong nelltions between x-, y- nd z-omponents) ould be importnt lmost down to lose ontt. This my provide se where experiments my be ble to selet the more pproprite wy to lulte intertions between exited tom pirs. As we hve demonstrted in this letter finite size effets tht hve been ignored influene the resonne intertion between two toms for seprtions of the order of few tomi rdii. More substntil long rnge effets re expeted for toms, moleules nd ions with lrger Gussin rdii. We end by giving in Fig.5 nd Fig.6 s further exmples the resulting Csimir-Polder intertion for pir of krypton toms, pir of phosphte ions, nd pir of p-4

5 Finite Size Dependent Dispersion Potentils 0 Kr-Kr in wter Iodide-iodide nd phosphte-phosphte in wter zero size finite size 0-5 Iodide finite size Iodide zero size Phosphte finite size Phosphte zero size Fig. 5: Color online) Retrded Csimir-Polder intertion between pirs of krypton toms dissolved in wter. Solid urve is for zero-size toms nd dshed for toms of finite size. The finite rdius of the krypton toms is =.4Å. iodide ions. Here, substntilly longer rnged finite size effets re observed. MB nd CP knowledge support from VR Contrt No. C04850) nd the Reserh Counil of Norwy Projet: 22469). PT knowledges support from the Europen Commission; this publition reflets the views only of the uthors, nd the Commission nnot be held responsiblefornyuse whih mybe mde ofthe informtion ontined therein. MB lso thnks the Deprtment of Energy nd Proess Engineering NTNU, Norwy) for finnil support. REFERENCES [] D. Kysylyhyn, V. Pitnytsi nd V. Lozovski, Phys. Rev. E 88, ). [2] M. Przybytek, W. Cenek, J. Koms, G. Lh, B. Jeziorski nd K. Szlewiz, Phys. Rev. Lett. 04, ). [3] M. Przybytek, B. Jeziorski, W. Cenek, J. Koms, J. B. Mehl nd K. Szlewiz, Phys. Rev. Lett. 08, ). [4] W. Cenek, M. Przybytek, J. Koms, J. B. Mehl, B. Jeziorski nd K. Szlewiz, J. Chem. Phys. 36, ). [5] Jr R. A DiStsio, V. V Gobre nd A. Tkthenko,.J. Phys.: Condens. Mtter, ). [6] Boström M., Longdell J. J., Mithell D. J., nd Ninhm B. W., Eur. Phys. J. D ) 47. [7] Boström M., Brevik I., Sernelius B. E., Dou M., Persson C., nd Ninhm B. W., Phys. Rev. A ) [8] Hrtmn R. L. nd Leung P. T., Phys. Rev. B ) [9] Sherkunov Y., Phys. Rev. A ) [0] MLhln A. D., Mole. Phys ) 409. [] Prsons D.F. nd Ninhm B. W., J. Phys. Chem. A ) 4. Fig. 6: Color online) Retrded Csimir-Polder intertion between pirs of iodide ions I ) nd pirs of phosphte ions PO 3 4 ) dissolved in wter. The urves ompre theory with nd without finite size. The sizes for iodide nd phosphte ions re = 2.2Å nd = 2.40Å, respetively. [2] Prsons D.F. nd Ninhm B. W., Lngmuir ) 86. [3] Boström M., Ellingsen S. Å., Brevik I., Prsons D. F., Sernelius Bo E., Phys. Rev. A ) [4] Werner H.-J., et l., MOLPRO, version 2008., pkge of b initio progrms, 2008). [5] Woon D. E. nd Dunning T. H., Jr., J. Chem. Phys ) [6] Peterson K. A. nd Dunning T. H., Jr., J. Chem. Phys ) [7] Wilson A. K., Woon D. E., Peterson K. A., nd Dunning T. H., Jr., J. Chem. Phys ) [8] Boroudjerdi H., Kim Y.-W., Nji A., Netz R.R., Shlgberger X., nd Serr A., Phys. Rep ) 29. [9] Lndu L. D. nd Lifshitz E. M., Eletrodynmis of Continuous Medi, Pergmon Press, Oxford, 960). [20] Mhnty J. nd Ninhm B. W., Dispersion Fores Ademi, London, 976). [2] Mhnty J. nd Ninhm B. W., J. Chem. So. Frdy Trns. II 7 975) 9 975). [22] D. D. Rihrdson, J. Phys. A: Mth. Gen. 8, ). [23] Buhmnn S. Y., Dispersion Fores. Mrosopi Quntum Eletrodynmis nd Ground-stte Csimir, Csimir- Polder nd Vn der Wls Fores Springer-Verlg, Berlin nd Heidelberg, 203). [24] Duignn T., Prsons D. F., nd Ninhm B. W., J. Phys. Chem. B 7 203) 942. p-5