Correlated exponential functions in high precision calculations for diatomic molecules. Abstract
|
|
- Jonas Dixon
- 6 years ago
- Views:
Transcription
1 Version 3.1 Correlated exponential functions in high precision calculations for diatomic molecules Krzysztof Pachucki Faculty of Physics University of Warsaw Hoża Warsaw Poland Abstract Various properties of the general two-center two-electron integral over the explicitly correlated exponential function are analyzed for the potential use in high precision calculations for diatomic molecules. A compact one dimensional integral representation is found which is suited for the numerical evaluation. Together with recurrence relations it makes possible the calculation of the two-center two-electron integral with arbitrary powers of electron distances. Alternative approach via the Taylor series in the internuclear distance is also investigated. Although numerically slower it can be used in cases when recurrences lose stability. Separate analysis is devoted to molecular integrals with integer powers of interelectronic distances r 1 and the vanishing corresponding nonlinear parameter. Several methods of their evaluation are proposed. PACS numbers: ac Ge krp@fuw.edu.pl 1
2 I. INTRODUCTION In order to achieve high accuracy for the nonrelativistic energy as well as for relativistic and quantum electrodynamics QED corrections in molecular systems well optimized basis functions have to be applied. The explicitly correlated Gaussian ECG functions do not satisfy the cusp condition therefore their use in the evaluation of higher order QED corrections 1 is problematic. Nevertheless the recent calculations of the leading Oα 3 QED effects performed for H using ECG functions lead to the most accurate to date theoretical predictions of about 10 3 cm 1 uncertainty for dissociation energies. To obtain even more accurate results and to include higher order QED effects we aim to use the basis of explicitly correlated exponential functions of the form φ = e w 1 r 1 u 3 r 1A u r 1B w r A w 3 r B r n 1 1 r n 1A rn 3 1B rn 4 A rn 5 B. 1 The notation for interparticle distances is explained after Eq. 4. The use of such exponential functions in molecular applications is quite limited. Starting from the pioneering work of Kolos and Roothan 3 the Neumann expansion of r 1 1 at ω = 0 in spherical oblate coordinates 4 has been applied most often see the most recent review by Harris in Ref. 5 and the collection of works in Ref. 6. This expansion and the more general one of e ω r 1 r 1 1 in terms of spheroidal functions 4 has been applied by several authors 7 11 in their accurate calculations for the hydrogen molecule. An alternative approach to perform integrals with exponential function in Eq. 1 has originally been proposed in Ref. 1. Authors have obtained an analytic although quite complex formula for the general four-body integral with exponential functions. By taking the inverse Laplace transform see Eq. 6 in one of the nonlinear parameters one can in principle obtain the general two-center two-electron integral in Eq.. Due to the very complicated analytic structure of the four-body formula this inverse Laplace transform has not been applied so far. In our recent work 13 we have reformulated the problem of the calculation of the inverse Laplace transform into the solution of some differential equations. From these differential equations satisfied by the master integral fr see Eq. below one derived analytic recursion relations for integrals with positive integer powers of interparticle distances assuming that the nonlinear parameter corresponding to r 1 vanishes i.e. for standard molecular integrals. Using these recursions one obtained analytic results for integrals with James-Coolidge and extended Heitler-London basis functions. The application
3 of analytic formulas has been demonstrated by the calculation of Born-Oppenheimer energies for H in Ref. 14 with accuracy of au and for HeH + in Ref. 15 with accuracy of about 10 1 au. While James-Coolidge and Heitler-London basis functions work very well for the ground states they are not equally good for excited states where arbitrary values of nonlinear parameters are needed. This is the main subject of this work to develop a computational technique for general two-center two-electron integrals which can be used for the calculation of nonrelativistic energies and also for relativistic and QED effects in two-electron diatomic molecule. Finally we think that various properties of the general twocenter two-body integral derived here can be applied for the calculation of nonrelativistic energies of an arbitrary diatomic molecule. This subject has recently been pursued also by Lesiuk and Moszyński in Ref. 16. The general structure of this work is the following. In Sec. II we derive the differential equation for the master integral fr in Eq.. In Sec. III we find a solution in terms of one-dimensional integral over elementary functions and additionally consider few special cases. In Sec. IV the Taylor series of fr in the internuclear distance r is derived and this shows the general analytic properties of the master integral f. In Sec. V integrals with positive powers of interparticle distances are obtained by differentiation of the master integral with respect to the corresponding nonlinear parameter. In Sec. VI we consider standard molecular integrals i.e. with the vanishing nonlinear parameter corresponding to r 1. In Sec. VII we present a summary with possible further applications. II. THE MASTER INTEGRAL We introduce here notation and definitions following the previous work in Ref. 13 and obtain the fourth order differential equation which is satisfied by the master integral fr. Matrix elements of the nonrelativistic Hamiltonian between functions of the type 1 can be expressed in terms of fr and its derivatives with respect to nonlinear parameters. The master two-electron two-center integral fr is defined by fr = r d 3 r 1 4 π d 3 r 4 π e w 1 r 1 r 1 e u3 r1a r 1A e u r1b r 1B e w ra and the related class of integrals with the nonnegative integer n is d 3 r 1 d 3 r 1 e u 3 r 1A e u r1b fr n = r 4 π 4 π r1 1 n r 1A r 1B r A e w ra r A e w3 rb r B e w3 rb r B 3 3
4 where 1 are indices of the electrons A B that of the nuclei and r = r AB is the distance between the nuclei. The notation for nonlinear parameters comes from the general four-body integral g defined by gu 1 = d 3 ρ 1 4 π d 3 ρ 4 π d 3 ρ 3 4 π e w 1 ρ 1 w ρ w 3 ρ 3 u 1 ρ 3 u ρ 31 u 3 ρ 1 ρ 3 ρ 31 ρ 1 ρ 1 ρ ρ 3 4 with ρ 1 = r 1 ρ = r A ρ 3 = r B. Function g is related to f by a Laplace transform namely d 3 r gt = 4 π fr e t r = dr fr e t r 5 r and the opposite relation is the inverse Laplace transform fr = 1 π i i +ɛ i +ɛ 0 dt e t r gt. 6 It was shown in Ref. 13 that the function g satisfies the first order differential equation in any of it s parameter α = u i w i : or equivalently σ g α + 1 σ α g + P α = 0 7 σ α σ g + P α = 0 8 with the inhomogeneous term P α presented in Appendix A and σ being the sixth order polynomial in six variables σ = u 1 u w 3 + u u 3 w 1 + u 1 u 3 w + w 1 w w 3 + u 1 w 1 u 1 + w 1 u u 3 w w 3 +u w u + w u 1 u 3 w 1 w 3 + u 3 w 3 u 3 + w 3 u u 1 w 1 w. 9 For example Eq. 7 for α = w 1 reads σ g + 1 σ g + P w 1 u 1 ; w u ; w 3 u 3 = 0 10 where the inhomogeneous term P w 1 u 1 ; w u ; w 3 u 3 is given by Eq. A. The solution of this differential equation is presented in the ingenious work of Fromm and Hill 1 by direct integration of Eq. in the momentum representation. A more compact formula was obtained by Harris 17. We do not present here their results as its explicit form is quite long and we will not need it. 4
5 To obtain the fundamental differential equation for the master integral f we use the differential equation 7 in variable t = u 1 σ g t + 1 σ t g + P t w 1; u 3 w 3 ; w u = 0 11 and perform the inverse Laplace transform in t. This differential equation using the following new parameters which are adapted to the symmetry of the problem w = w + x w 3 = w x u = u y u 3 = u + y 1 takes the form d σ 4 d r r d d r + σ d d r r d d r + σ 0 r fr = F r 13 where σ = σ 0 + t σ + t 4 σ 4 14 σ 4 = w 1 σ = w 4 1 w 1 u + w + x + y + 16 u w x y = w w 1 σ + σ 0 σ 0 = w 1 u + w x y u w + x y u w x + y u + w + x + y +16 w x u y u x w y u w x y = w 1 σ 0 + σ 00 and F r = F u1 r is presented in Appendix A. The differential Eq. 13 is supplemented by the boundary conditions namely fr vanishes at small and large r. Taking the occasion we present here for later use a differential equation which is obtained from the inverse Laplace transform of Eq σ 0 1 σ +σ 0 fr+ +σ f r+ w 1 +w 1 f 4 r = F w1 r 15 where F w1 r is given in Appendix A and that for an arbitrary parameter α = u w x y 1 σ 0 α + σ 1 σ 0 fr + α α + σ f r + w 1 α α f 4 r = F α r. 16 III. POISSON REPRESENTATION We derive here the integral representation for solutions of the homogeneous and the inhomogeneous differential equation 13 which is the analog of the Poisson representation 5
6 for Bessel functions 18. σ is a quadratic polynomial in t so it has four zeros ±t a ±t b. Assuming t b > t a > 0 the four solutions of the homogeneous equation d σ 4 d r r d d r + σ d d r r d d r + σ 0 r fr = 0 17 are tb e t r σt dt ta t b e t r σt dt Let us prove it as an example for the last term d σ 4 d r r d d r + σ tb = = = 1 t a tb t a tb ta t a d d r r d d r + σ 0 r dt d σ 4 σt d r r d d r + σ e t r σt dt tb dt t a d d r r d dt σt σ 4 t d dt t + σ t d dt t + σ 0 t a d dt tb t a e t r σt d r + σ 0 r e t r e t r σt dt. 18 dt σt d σt dt + d dt σt e t r 19 where the last equation holds because σt = σ 0 + σ t + σ 4 t 4. One integrates by parts boundary terms vanish because σt a = σt b = 0 and one obtains = 1 tb d dt e t r d σt + σt dt dt = 1 t a tb t a dt e t r 1 σt 1 σt e t r 1 1 = 0. 0 The solution of the inhomogeneous fundamental differential equation 13 is obtained by noting that the integration contour in the inverse Laplace transform can be deformed to encircle all branch cuts on the left side of the complex plain fr = 1 π i 0 dt e t r gt + i ɛ gt i ɛ. 1 The analytic properties of g have been analyzed by Fromm and Hill in 1. The function gt has four overlapping branch cuts similarly to the function P in variable t on the negative real axis starting at t 1 = u 3 + w = u + y + w + x t = u + w 3 = u y + w x t 3 = u 3 + w 1 + w 3 = u + y + w x + w 1 t 4 = u + w 1 + w = u y + w + x + w 1 6
7 correspondingly. The fact that all t i with i = are not positive comes from the requirement that the integral in Eq. is finite for positive values of u i and w i. Fromm and Hill found the imaginary part of g for the particular ordering of t i. Their result is generalized here to an arbitrary ordering of t i and takes the form 1 1 gt + i ɛ gt i ɛ = π i θt1 t ln β 00 + θt t ln β 33 σ θt 3 t ln β 31 β 33 θt 4 t ln β 01 β 00 3 where β ij = σ γij σ + γij 4 and cf. Ref.1 with a 1 = u 3 a = w a 3 = t a 3 = w 3 a 31 = u a 1 = w 1 γ 00 = u w 1 w 3 + u u 3 + w1 w 3 + w 1 t + u + w3 + u w1 w + w3 γ 33 = u 3 w 1 w + u + u 3 + w1 w + w 1 t + u 3 + w + u 3 w1 + w w3 γ 01 = t w w 3 t u 3 + w w 3 + w t u + w3 + t w1 + w + w3 γ 31 = t w w 3 + t u 3 + w w 3 w t u + w3 + t w1 + w + w3. 5 There is some arbitrariness in the form of Eq. 3 as products of β can be expressed in several ways 1 17 β µµ β µν = β νν β νµ = β ρρ β ρσ β µσ β νσ β ρσ = β σσ 6 for {µ ν ρ σ} being an arbitrary permutation of {0 1 3}. One can verify that Img in Eq. 3 satisfies the differential equation 11 with ImP as an inhomogeneous term. As a result fr can be expressed in terms of the one-dimensional integral fr = 0 dt e t r 1 θt1 t ln β 00 + θt t ln β 33 σ θt 3 t ln β 31 β 33 θt 4 t ln β 01 β This is a main result obtained in this work. All properties of the master integral f can in principle be obtained from this integral representation. Some of them however can be more easily obtained from the differential equation 13. The rest of the paper will be devoted to various properties and various ways of numerical calculation of fr in Eq. 7 and its 7
8 derivatives with respect to nonlinear parameters. This form of Eq. 7 is suited for the direct numerical quadrature once we know all the singularities on the integration path and this is presented in Appendix B. It would be interesting to investigate various expansions of fr in Eq. 7 in particular the expansion around w 1 = 0 is studied in Sec. VI. Here we obtain the master integral f at w 1 = 0 in a simple form which is convenient for the numerical evaluation. The σ polynomial reads σ w1 =0 = u u 3 + w w 3 u w u 3 w 3 t u u 3 w w 3 = 16 w x u yu x w yu w x y + 16 t u w x y 8 and the differential equation becomes u w x y d d r r d + 16 r w x u y u x w y u w x y fr = F r d r 9 w1 =0 where the inhomogeneous term F r using Eq.A5 is F r w1 =0 = u w x u w y u x y + w x y F 1 + u w x u w y + u x y w x y F u w x + u w y + u x y + w x y F 3 + u w x + u w y u x y w x y F 4 30 and where F 1± = Ei r u e r u w x+y ± Ei r w e r u+w+x y F ± = Ei r w e r u+w x+y ± Ei r u e r u w+x y F 3± = Ei r x + Ei r y Ei r x + y + ln F 4± = u w x + y x y u + w r u+w+x+y e ±Ei r u + w e r u+w+x+y u w x + y r u+w x y Ei r x + Ei r y Ei r x + y + ln e x y u + w ±Ei r u + w e r u+w x y 31 with Ei being the exponential integral function. The solution of this differential equation assuming p > 0 in terms of solutions of the homogeneous equation was presented in Ref fr = I 0 p r 16 u w x y r dr F r K 0 p r + K 0 p r 8 r 0 dr F r I 0 p r 3
9 where p = w x u y w y u x u w x y 33 u w x y and I 0 K 0 are modified Bessel functions. A more convenient however is the integral form of Eq. 7 fr = u w x y u w+x y ln β 00 + ln β 31 β 33 u w+x+y ln β 33 u w x+y e t r ln β 01 β 00 σ dt 34 with the following γ ij coefficients γ 00 γ 33 = 4 u w x + u w y u x y w x y = 4 u w x + u w y + u x y + w x y Γβ 31 β 33 = 16 u w x y/t + u + w x + y + 4 u w x u w y + u x y + w x y Γβ 01 β 00 = 16 u w x y/t + u + w + x y 4 u w x u w y u x y w x y 35 where Γ is defined in Eq. B3. In the particular case of exchange integral w = u x = y σ 00 vanishes and the master integral becomes fr = 1 dt ln t + w x e t r 8 w x w+x t w x t + dt ln t w + x e t r w x t + w + x t dt ln t + w x e t r t + w + x. 36 t w An equivalent integral representation obtained from the differential equation 9 was presented in Ref. 13 but this form is more convenient for the numerical evaluation. IV. TAYLOR EXPANSION IN r Here we find the Taylor expansion in the internuclear distance r of the master integral f. As it was noticed by Fromm and Hill in 1 this expansion is absolutely convergent therefore it is another way to calculate the master integral as well as its derivatives with respect to u i and w i. Following Ref. 1 we find the initial terms and the recurrence relation for subsequent terms of the Taylor series. 9
10 The expansion in small r of the master integral has the form fr = r k f 1 k lnr + f k. 37 k=1 The leading term can be obtained by noting that the limit r 0 corresponds to the one-center integral 1 1 d 3 r 16 π 1 d 3 e α r 1 β r γ r 1 r r1 r r 1 = 1 π γ α + γ ln β + γ + L 1 α + β α + γ + L 1 β + α 38 β + γ where L is the dilogarthmic function 18. Let us introduce for the later use the following X i symbols X 0 = 1 π w u + w1 ln w + w 1 u + w1 X 1 = ln u + w w + w1 X = ln u + w + L 1 u + w + L 1 u + w 1 u + w w + w 1 X r = 1 lnr u + w 1 w + w 1 + γ E. 39 where γ E is the Euler constant. The leading terms in small r expansion using Eq. 38 are f 1 1 = 0 f 1 = X and f i 0 = f i 1 = 0. The next terms of this expansion f 1 f are obtained from the large t asymptotics of the function g g = X 0 t + g1 ln t + g t 3 + o t Since gt satisfies the differential equation 11 with the inhomogenous term P t = P t w 1 u 3 w 3 w u P t = w 1 the 1/t 3 coefficients are + lnu + u 3 + w 1 + lnw 1 + w + w 3 lnt + o t 1 4 g 1 = g = lnu + u 3 + w 1 + lnw 1 + w + w
11 As noticed by Fromm and Hill 1 they are related to f i coefficients by f 1 = g1 = ψ3 g1 + g 44 f where ψ is the Euler ψ-function 18 and ψ3 = γ E +1+1/. The next terms in r expansion can be obtained from the fundamental differential equation 13. As the inhomogeneous term F r has a similar expansion F r = k= 1 the recurrence relations are the following r k F 1 k lnr + F k 45 f 1 k+1 = F 1 k σ 0 f 1 k 3 σ k 1 f 1 k 1 46 σ 4 k + 1 k k 1 f k+1 = F k σ 0 f k 3 σ k 1 f k 1 σ k 1 f 1 k 1 σ 4 k k + 1 k f 1 k+1 σ 4 k + 1 k k 1 The resulting first terms in the expansion of the master integral f in powers of r are fr = r X 0 + r X r 3 +r 3 w 1 X 0 w u + w + x + y X 0 w X 1 u X u + w 6 + x y x y u w X 0 + u X 1 + w X 3 w 1 3 w1 7 +r 4 54 w 1 + w 1 u + w + 3 u w 3 x y 11 u 11 w 9 x 9 y X r w u + w + x + y + Or Let s consider now the special case of w 1 = 0. The recurrence relation takes the form: k+1 = F 1 k σ 00 f 1 k 1 σ 0 k + 1 f 1 k+1 = F k σ 00 f k 1 σ 0 k + 1 k + 1 f 1 k+1 48 f 11
12 and the first terms in r-expansion are 1 fr w1 =0 = r 1 u ln + u w w ln + w + r u +r 3 u 3 w 3 + x y u 18 u + x y 18 w + w u + x + y u x y 1 w 18 w u w x y +r 4 1 +Or 5. w + u 18 In order to obtain the Taylor expansion of fr n 3 + γ + 1 ln4 r u w 4 + w + x + y w x y ln 1 + u 1 u 18 u w ln 1 + w u 3 + γ + 1 ln4 r u w u + w + x + y 49 fr n = k=1 r k f 1 k n lnr + f k n 50 one expands Eq. 15 in power series in r and obtains the following recursions for coefficients f 1 k+ n = 1 σ 00 f 1 k σ 0 k + k + 1 n + σ 0 n 1 f 1 k n +σ n 1 k + k + 1 f 1 k+ n +n 1 k + 4 k + 3 k + k + 1 f 1 k+4 n +n 3 n n 1 k + k + 1 f 1 k+ n 4 f k+ n = 1 σ 0 k + k + 1 F w1 kn 1 σ 0 k + 3 f 1 k+ n σ 00 f k n σ 0 n 1 f k n σ n 1 k + k + 1 f n 3 n n 1 k+ n + k + 3 f 1 k+ n 51 k + k + 1 f 1 k+ n 4 + k + 3 f k+ n 4 n 1 k + 4 k + 3 k + k + 1 f k+4 n + k + 5 k + 5 k + 5 f 1 k+4 n 5 where f i i 0 = f. The usefulness of the calculation of fr and fr n via Taylor series k k in r and the above recursions needs to be verified numerically. Nevertheless the expansion terms say a lot about analytic properties of f as a function of w 1 x y u and w. In particular the expansion in Eq. 47 suggests that f is an entire function in x y at fixed r u w w 1. 1
13 V. GENERAL RECURSION RELATIONS In order to use the exponentially correlated basis in molecular calculations one has to be able to calculate derivatives of the master integral f with respect to nonlinear parameters. Differential equations 1315 and 16 are used below to express the derivative of fr with respect to an arbitrary parameter α = w 1 x y u w in terms of fr f r f r and f r. This is done as follows. From Eq. 13 the fourth derivative f 4 r can be expressed in terms of lower derivatives of fr and we use this in Eqs. 15 and 16 to eliminate fourthth derivative. Next differentiate these equations with respect to r and again eliminate fourth derivative. Doing this differentiation and elimination three times one obtains four equations for four unknowns: the derivative of f with respect to w 1 and its three further derivatives with respect to r. The solution of these linear equations for α = w 1 is σ 0 σ 0 σ w 1 σ δ σ 0 fr = r δ σ 0 f r + w 1 and for α = u w x y σ 0 4 w 1 σ 0 + w 1 δ σ 0 fr α = r δ σ 0 α f r w 1 σ 0 w 1 σ 0 σ σ σ 0 α σ σ σ 0 σ α f r + r f r σ + w 1 σ 0 w 1 fr + r f r σ f r + r f r α α σ 0 fr + r f r where... denotes inhomogeneous terms presented in Eqs. C3C5 and where δ = 4 w 1 σ 0 σ = w 1 4 u w 1 4 w w 1 4 x w 1 4 y. 55 The solution for higher order derivatives with respect to r can be represented in the matrix form. Let us introduce a symbol f to denote fr f fr r =. 56 f r f r Then f α = Âα f + F α 57 13
14 where Âα and F α are presented in Appendix C. These vector differential equations with respect to the nonlinear parameter α can be used as recursion relations to express the integral with any powers of electron distances in terms of the master integral fr and it s first three derivatives with respect to r. These recursions become unstable when w 1 0 δ 0 or σ 0 0. These cases require a separate analysis. In the next Section we consider recursion relations for w 1 = 0 which has applications for standard molecular integrals. The particular case of σ 0 = 0 happens for an exchange integral where w = u and x = y. Here one can use recurrence only for f r to generate the analytic expressions for derivatives with respect to w u and x y. They are quite compact as there will be only three independent parameters and can be expressed in terms of fr and f r. The derivatives with respect to remaining nonlinear parameters can be calculated afterwards by setting w = u x = y. VI. RECURSIONS FOR STANDARD MOLECULAR INTEGRALS We describe in this Section the evaluation of standard molecular integrals where w 1 = 0. Recursion relations for the integral fr n in Eq. 3 have been presented in Ref. 13. Here we rederive them using a more compact notation. In particular cases close to spurious singularities we propose to use these recursions to derive an analytical expression for individual integrals which are then implemented in the numerical codes. Let us take Eqs. 13 and 15 and differentiate them with respect to w 1 n and n 1 times respectively at w 1 = 0 r σ 00 fr n + σ 0 f r n + r σ 0 f r n + n 1 n r σ 0 fr n + n 1 n σ f r n + n 1 n r σ f r n + n 1 n f 3 r n + n 1 n r f 4 r n + n 3 n n 1 n f r n 4 + n 3 n n 1 n r f r n 4 = F u1 r n 58 σ 00 fr n + σ 0 f r n + n 1 σ 0 fr n + n 1 σ f r n + n 1 f 4 r n + n 3 n n 1 f r n 4 = F w1 r n These are two linear equations for three unknowns fr n f r n f r n. The third equation is obtained by elimination of f r n and further differentiation with respect to r. The 14
15 solution of these three equations for fr n and f r n are fr n = 1 σ 00 n n 1 σ0 fr n + n 1 r σ 0 f r n + n 1 σ f r n + n 1 r σ f 3 r n +n 1 n + f 4 r n + n 1 r f 5 r n +4 n 3 n n 1 f r n 4 + n 3 n n 1 r f 3 r n 4 + F w1 r n 1 + r F w 1 r n 1 F u 1 r n 60 f r n = 1 σ 0 n 1 r σ0 fr n + n 3 n n 1 n f r n 4 +n 1 n σ f r n + n 3 n n 1 r f r n 4 where +n 1 r σ f r n + n 1 n f 3 r n + n 1 r f 4 r n +r F w1 r n 1 F u1 r n 61 F X r n = 1 n n w n 1 F X r 6 w1 =0 for X = w 1 u 1. Equation 60 allows one to obtain integral fr n with an arbitrary power n 1 of r 1 in terms of fr for example fr 0 = fr fr 1 = r3 4 h 0r u h 0 r w j 0 r x j 0 r y fr = 1 r σ0 f r + σ f r + r σ f 3 r + 4 f 4 r + r f 5 r σ 00 + F w1 r 1 + r F w 1 r 1 F u 1 r fr 3 = r5 3 h 1 r u h 1 r w j 1 r x j 1 r y h 0 r u h 0 r w j 0 r x j 0 r y 4 +h 0 r w h r u j 0 r x j 0 r y + h 0 r u h r w j 0 r x j 0 r y +h 0 r u h 0 r w j r x j 0 r y + h 0 r u h 0 r w j 0 r x j r y 63 where j n and h n are modified spherical Bessel functions 1 j n x = x n d n sinhx x dx x 1 h n x = x n d n exp x. 64 x dx x One notices that fr n for odd n can be expressed in terms of j n and h n only and their numerical evaluation is straightforward. In contrast numerical evaluation of the formulas 15
16 for even n is much more difficult. If σ 0 is not small then derivatives of fr can be reduced to fr and f r using f n r = 1 r σ 0 σ 00 n f n 3 r+r f n r +σ 0 n 1 f n 1 r F n If σ 0 is small then the inverse recursion f n 1 r = r σ00 f n+1 r + n + σ 0 f n+ r + r σ 0 f n+3 r F u n+1 n + 1 σ 1 r 00 is stable and can be used to obtain all the derivatives including the function fr itself. Alternatively f n r can be obtained directly using the integral form of Eq. 34. If σ 00 is small but not σ 0 then these formulas can be rewritten using Eq. 61 to the form without σ 00 in the denominator for example fr = 1 σ 0 σ0 f r r σ 0 f 1 r σ fr r σ f r 3 f r r f r where f n r is defined by f n r = r 0 u 1 r f n+1 r dr 68 and can be obtained from Eqs. 34 or 66. It may lead however to other type of numerical instabilities if the branch point starts at 0 see Eq. and in this case one can build the explicit table of integrals as in Ref When both σ 0 0 and σ 00 0 then parameters x y are small and in this case the Neumann expansion 5 can be applied. What remains are the derivative of f with respect to α = u w x y at w 1 = 0. We adapt here the derivation of corresponding formulas from Ref. 13. One takes fundamental differential equations in u 1 and α at w 1 = 0 r σ 00 fr + σ 0 f r + r σ 0 f r = F u1 r 69 1 σ 00 α fr + σ fr 00 α + 1 σ 0 α f f r r + σ 0 α = F αr 70 differentiates the first equation with respect to α eliminates f r/α and differentiate resulting equation again with respect to r. The obtained equation for the derivative of f using σ 0 /α = σ 0 /α is f α = 1 α fr + r 1 σ 00 σ 00 α 1 f r + G αr 71 α σ 00 16
17 where and G α r = 1 α F u α 1 r + r d F α r dr w 1 =0 G w r = 1 u y w x u x w y F 1+ F + + F 3+ + F 4+ w +u + w x y u w x y F 1+ + F + +u w x + y u w x y F 3+ + F 4+ G u r = 1 u y w x u x w y F 1+ F + + F 3+ + F 4+ u +u + w x y u w x y F 1+ F + +u w x + y u w x y F 3+ + F 4+ G x r = 1 u y w x u x w y F 1+ + F + F 3+ F 4+ x +u + w x y u w x y F 3+ + F 4+ +u w x + y u w x y F 1+ + F + G y r = 1 u y w x u x w y F 1+ + F + F 3+ F 4+ y 7 +u + w x y u w x y F 3+ F 4+ +u w x + y u w x y F 1+ F Here σ 0 from the denominator can be removed when necessary as in the case of fr n by introducing f n r. The formula 71 allows one to obtain integrals with arbitrary powers of electronic distances. It requires however further investigations of various numerical instabilities when the denominator approaches 0. In particular σ 00 vanishes for exchange integrals. In this case w = u x = y and one can generate a table of integrals avoiding recursions by Taylor expansion of the fundamental solution in Eq. 7 in w 1 w u and x y. Alternatively one can take recursion relations presented in Eq. 61 and integrate it once over r. VII. SUMMARY The evaluation of two-center two-electron integrals is a difficult task. Depending on a physical problem different approaches are employed. For the calculation of nonrelativistic energies of a two-electron diatomic molecule the best way is by the use of Kolos-Wolniewicz 17
18 basis functions namely w 1 = 0 and arbitrary nonlinear parameters u w x y. The master integral can be calculated numerically using Eqs. 34 and 37. Integrals with powers of interparticle distances can be obtained by using analytic recursions Eq and 71 followed by numerical evaluation of obtained formulas. For the calculation of QED effects in a two-electron diatomic molecule the best way is to use the most general basis functions with nonvanishing w 1 u w x y. The master integral can be evaluated according to Eq. 7 and integrals with additional positive powers of interparticle distance can be obtained through recursion relations Eq. 57. Integrals with negative powers can be obtained by numerical integration with respect to the corresponding nonlinear parameter. In order to demonstrate correctness of obtained formulas for twocenter two-electron integrals Table I presents a simple examplary calculations of the ground electronic Σ + g state of the H molecule. TABLE I: Examplary calculations of the Born-Oppenheimer energy at the equilibrium distance R = 1.4 of the ground Σ + g state of H molecule obtained using functions of the form Eq. 1 with parameters w 1 = 0.5 u = 1 w = 1 y = 0.15 x = 0.15 defined in Eqs. 1 N is a number of basis functions and n 1 + n + n 3 + n 4 + n 5 Ω. The result with is accurate to all digits and was obtained in Ω N E BO a.u An additional approach to evaluate two-center integrals is by the Taylor expansion in small internuclear distance r. Its applicability has to be verified numerically. Nevertheless it serves at present as a simple check of various formulas for recursion relations and the master integral. The advantage of calculations by Taylor series is that expansion terms are quite simple. They are related to the known one-center helium-like integrals and can be obtained recursively A similar Taylor expansion can in principle be obtained for many-body three-linked integrals at most three odd powers of interparticle integrals. 18
19 This would open a window for high precision results in many electron diatomic molecules. Acknowledgments I wish to thank Janek Dereziński for useful information about properties of Bessel functions. This work was supported by NCN grant 01/04/A/ST/ K. Pachucki Phys. Rev. A J. Komasa K. Piszczatowski G. Lach M. Przybytek B. Jeziorski and K. Pachucki JCTC W. Kolos and C. C. J. Roothaan Rev. Mod. Phys ; ibid C. Flammer Spheroidal Wave Functions Stanford University Press Stanford F.E. Harris Int. J. Quant. Chem Recent Advances in Computational Chemistry: Molecular Integrals over Slater Orbitals Ed. Telhat Ozdogan and Maria Belen Ruiz Transworld Research Network 008 India. 7 W. Ko los J. Chem. Phys L. Wolniewicz J. Chem. Phys J. S. Sims and S. A. Hagstrom J. Chem. Phys H. Nakatsuji Phys. Rev. Lett A. Ferron and P. Serra J. Chem. Theor. Comp D. M. Fromm and R. N. Hill Phys. Rev. A K. Pachucki Phys. Rev. A K. Pachucki Phys. Rev. A K. Pachucki Phys. Rev. A M. Lesiuk and R. Moszyński arxiv: quant-ph 17 F.E. Harris Phys. Rev. A NIST Handbook of Mathematical Functions F.W.J. Olver D.W. Lozier R.F. Boivert and C.W. Clark Cambridge University Press
20 Appendix A: Inhomogeneous terms The first order differential equation 8 satisfied by the general four body integral defined by Eq. 3 involves the inhomogeneous term P α for α = u i w i. They are related to each other by P w1 = P w 1 u 1 ; w u ; w 3 u 3 = P w 1 u 1 ; w 3 u 3 ; w u P u1 = P u 1 w 1 ; w u ; u 3 w 3 P w = P w u ; w 3 u 3 ; w 1 u 1 P u = P u w ; w 3 u 3 ; u 1 w 1 P w3 = P w 3 u 3 ; w 1 u 1 ; w u P u3 = P u 3 w 3 ; u 1 w 1 ; w u A1 The expression for P was obtained in 13 and is the following = P w 1 u 1 ; w u ; w 3 u 3 u 1 w 1 u 1 + w u 3 u 1 + u 3 w u 1 + u 3 + w ln u + u 3 + w 1 u 1 + u + w 1 + w u 1 w 1 u 1 + u 3 w + u 1 u 3 + w u 1 + u 3 + w ln w 1 + w + w 3 u 1 + u 3 + w 1 + w 3 u 1 w1 + u w u 3 w3 + w 1 w u 1 + u w3 u 1 + u + w 3 ln w 1 w + w 3 w 1 + w + w 3 u 1 + u + w 1 + w u 1 w1 u w + u 3 w3 + w 1 w 3 u 1 + u 3 w u 1 + u 3 + w ln w 1 + w w 3 w 1 + w + w 3 u 1 + u 3 + w 1 + w 3 + u u + w 1 u 1 + u 3 w u 3 u 1 + u w3 u 1 + u 3 + w ln u + u 3 w 1 u + u 3 + w 1 u 1 + u + w 1 + w + u 3 u 3 + w 1 u 1 + u w3 u u 1 + u 3 w u 1 + u + w 3 ln u u 3 w 1 u + u 3 + w 1 u 1 + u 3 + w 1 + w 3 w 1 w u 1 u + w3 + w 3 u 1 u 3 + w u + u 3 + w 1 ln w 1 w w 3 w 1 + w + w 3 u + u 3 + w + w 3 w 1 u u 1 + u 3 w + u 3 u 1 + u w3 w 1 + w + w 3 ln u u 3 + w 1 u + u 3 + w 1 u + u 3 + w + w 3 A 0
21 In the case the nonlinear parameter α is the combination of w i u i as in Eq. 1 then P α is given by P w = P w + P w3 P x = P w P w3 P u = P u + P u3 P y = P u3 P u A3 The inhomogeneous term in the fundamental differential equation 13 is the inverse Laplace transform of P In particular F u1 F u1 r = F α = 1 π i is given by { 1 w 1 r + w 1 + u + w x + y r i +ɛ i +ɛ dt e t r P α u1 =t e r u+w+w 1 x+y A4 1 w 1 r + u + w + x + y r u+w+x+y e r w + 1 u w x + y + u w x y + x y w u Ei r u + w 1 e r u w x+y Ei r w + w 1 e r u w x+y w + 1 u + w + x + y + u w x + y + x y u + w Ei r x + y Ei r w 1 x Ei r w 1 y } x + y w1 + u w 1 + w ln e r u+w+x+y + Ei r u + w e r u+w+x+y u + w w 1 x w 1 y { } + x x y y A5 1
22 where Ei is the exponential integral function and F w1 is { r u + w + x + y + 1 F w1 r = w 1 w 1 w 1 r 3 u + w 1 w + w 1 w 1 x w 1 w 1 y + u + w u + x u + y u + w w + x w + y + r u + w 1 w + w 1 u + x w + x x + y u + y w + y x + y u + w + x + y w 1 x w 1 y r u+w+x+y e r u + w + w1 x + y r 3 + w 1 u + w 1 + w 1 w + w 1 + w 1 w 1 x + w 1 w 1 + y + w 1 r + 1 w x y + u x + w x u y w y u w + w 1 + u + w x + y r u w u + x u y u w w + x w y + u + w 1 w + w 1 } u + x w + x x + y u y w y x y + + e r u+w+w 1 x+y w 1 x w 1 + y { } + x x y y A6 One notes that F w1 involves only the exponential and rational functions. Appendix B: Numerical evaluation of the master integral For the numerical evaluation of the master integral in Eq. 7 one should find all singularities on the integration path 1. For this one decomposes σ γij into products σ γ00 = u u 3 + w 1 u + u 3 + w 1 t u w 3 t + u + w 3 w 1 w + w 3 w 1 + w + w 3 σ γ33 = u u 3 w 1 u + u 3 + w 1 t + u 3 + w t + u 3 + w w 1 + w w 3 w 1 + w + w 3 σ γ31 = t + u 3 w t + u 3 + w t + u w 3 w 1 + w w 3 t + u + w 3 w 1 w + w 3 σ γ01 = t + u 3 w t + u 3 + w t + u w 3 w 1 + w w 3 t + u + w 3 w 1 w + w 3 B1 and observes that β 00 vanishes on the negative axis at t = t β 33 at t = t 1 β 31 at t = t and β 01 at t = t 1. All singularities at other points t u i w j t+u i w j and t u i +w j cancel
23 out among various ln β s because the original integral Eq. is finite at these points and g is analytic. The cancellation can be achieved explicitly by combining the corresponding β s under the common logarithm. For this one has to consider all possible orderings of t i. Because of the symmetries 1 and A B one can assume without loosing generality that t 1 > t and t 3 > t 4. We remain with six possible orderings of t i t 1 t t 3 t 4 t 1 t 3 t t 4 t 1 t 3 t 4 t t 3 t 1 t t 4 t 3 t 1 t 4 t t 3 t 4 t 1 t and six corresponding representations of the master integral f fr = = = = = = t3 t4 e t r ln β 00 + ln β 00 β 33 + ln β 00 /β 31 ln β 01 β 31 dt t 1 t t 3 t 4 σ t t4 e t r ln β 00 ln β 0 + ln β 00 /β 31 ln β 01 β 31 dt t 1 t 3 t t 4 σ t4 t e t r ln β 00 ln β 0 ln β 00 β 01 β 0 ln β 01 β 31 dt t 1 t 3 t 4 t σ t1 t t4 e t r ln β 00 β 0 ln β 0 + ln β 00 /β 31 ln β 01 β 31 dt t 3 t 1 t t 4 σ t1 t4 t e t r ln β 00 β 0 ln β 0 ln β 00 β 01 β 0 ln β 01 β 31 dt t 3 t 1 t 4 t σ t4 t1 t e ln β 00 β 0 ln β00 t r β 01 β 0 ln β 00 β 01 β 0 ln β 01 β 31 dt t σ t t3 t3 t 3 t 4 t 1 B In order to eliminate singularities algebraically the products and ratios of β can be combined together 17 where σ γ σ + γ γ σ γ + σ = Γ = σ + γ γ γ + γ σ Γ σ + Γ B3 B4 3
24 As a result no singularity appears on the integration path which is demonstrated by the following decompositions σ Γ β 01 β 31 = t u 3 w t + u 3 w t u 3 + w t + u 3 + w t u w 3 t + u w 3 t u + w 3 t + u + w 3 /4 t σ Γ β 00 /β 31 = u u 3 + w 1 u + u 3 + w 1 t u 3 w t + u 3 w t u w 3 w 1 + w w 3 t u + w 3 w 1 + w + w 3 /t u w 1 w σ Γ β 00 β 0 = u + u 3 + w 1 u + u 3 + w 1 t + u 3 w t + u 3 + w t u + w 3 t + u + w 3 w 1 w + w 3 w 1 + w + w 3 /t + u 3 + w 1 + w 3 σ Γ β 00 β 01 = u u 3 + w 1 u + u 3 + w 1 t u 3 + w t + u 3 + w t + u w 3 w 1 + w w 3 t + u + w 3 w 1 + w + w 3 /t + u + w 1 + w B5 Table II presents numerical results for f for several internuclear distances r and nonlinear parameters w 1 x y u and w. For r < 0.01 f is evaluated from the Taylor series up to r 8 and for r > 0.01 from the integral representation in Eq. B. TABLE II: Master integral f for selected values of internuclear distances r and nonlinear parameters w 1 u w u 3 w 3 f f f
25 Appendix C: Exact form of general recursion relations The recursion relations for the general master integral as given by Eq. 57 involves the 4 by 4 matrix A α. For α = x y u w A is the following 1 δ r 4 α 4 A α = 1 δ r σ 0 4 w 1 δ δ σ 0 α σ 0 w 1 σα α σ 0 r w 1 σα σ 0 r σ α r 4 σ α δ α r σ α 1 4 δ σ r δ w1 α w1 δ α δ α σ α r σα r 4 δ α C1 where and the inhomogeneous term F α is F α = 1 δ σ σ σ 0 α σ α σ σ α σ w1 + 1 δ r σ σ α = σ w 1 w 1 σ 0 α + w r 1 δ σ 0 α w 1 σ α + σ 0 3 r α rα w 1 w 1 δ σ 0 + r α + σ 4 r 3 α σ 0 σ α rα + w 1 σ σ 0 3 r 3 3 r α w 1 δ σ 0 r σ 8 w 1 r r w 1 r σ r σ r 10 w 1 r σ σ r + r σ r r r w 1 F u1 r r + 5 w 1 σ σ 0 r + r w 1 σ σ 0 3 r 3 3 r 3 F α r C C3 For α = w 1 matrix A is A w1 = 1 δ r σ 0 δ w r σ 0 4 w 1 δ r 4 r σ w1 δ δ σ 0 σ 0 δ r 4 σ w1 r σ w1 1 4 δ δ σ w 3 1 δ w 1 w 1 σ r δ σ w 3 1 r δ w 1 σ r σw 1 δ w 1 w 1 σw 1 σ 0 r w 1 σw 1 σ 0 δ r δ w 1 δ δ r δ w 1 4 r δ w 1 C4 5
26 and the inhomogeneous term 1 δ 4 σ 0 4 w1 σ 0 + F w1 = 1 σ σ w 1 σ δ σ σ σ w 1 r + 1 δ δ + 1 w1 3 4 w1 δ σ 0 w 1 δ σ 0 + r w 1 δ σ 0 + w 1 w r 1 r w 1 σ σ 0 3 r 3 + w 1 σ σ 0 4 r 3 3 r r + w 3 1 w r σ 0 w1 δ σ 0 r σ 8 w 1 r r w 1 r σ r σ r 10 w 1 r σ σ r + r σ r 4 r 3 σ w 1 r + σ 3 r r r w 1 r + 5 w 1 σ σ 0 r + r w 1 σ σ 0 3 r 3 3 r 3 F w1 r F u1 r C5 6
Accurate adiabatic correction in the hydrogen molecule. Abstract
Version 1.2 Accurate adiabatic correction in the hydrogen molecule Krzysztof Pachucki Faculty of Physics, University of Warsaw, Pasteura 5, 02-093 Warsaw, Poland Jacek Komasa Faculty of Chemistry, Adam
More informationarxiv: v1 [physics.chem-ph] 24 Apr 2017
Relativistic corrections for the ground electronic state of molecular hydrogen Mariusz Puchalski and Jacek Komasa Faculty of Chemistry, Adam Mickiewicz University, Umultowska 89b, 61-614 Poznań, Poland
More informationSchrödinger equation solved for the hydrogen molecule with unprecedented accuracy. and Jacek Komasa
Schrödinger equation solved for the hydrogen molecule with unprecedented accuracy Krzysztof Pachucki 1, a) 2, b) and Jacek Komasa 1) Faculty of Physics, University of Warsaw, Pasteura 5, 02-093 Warsaw,
More informationLowest vibrational states of 4 He 3 He + : Non-Born-Oppenheimer calculations
Lowest vibrational states of He 3 He + : Non-Born-Oppenheimer calculations Monika Stanke,, * Dariusz Kȩdziera, Sergiy Bubin, Marcin Molski, 3 and Ludwik Adamowicz, Department of Chemistry, University of
More informationHigh Precision Variational Calculations for the Born-Oppenheimer Energies of the Ground State of the Hydrogen Molecule. Abstract
High Precision Variational Calculations for the Born-Oppenheimer Energies of the Ground State of the Hydrogen Molecule James S. Sims National Institute of Standards and Technology Gaithersburg, MD 20899-8911
More informationHelium fine structure theory for determination of α
Journal of Physics: Conference Series Helium fine structure theory for determination of α To cite this article: Krzysztof Pachucki and Vladimir A Yerokhin 0 J. Phys.: Conf. Ser. 64 0007 View the article
More informationOn the derivatives 2 P ν (z)/ ν 2 and Q ν (z)/ ν of the Legendre functions with respect to their degrees
INTEGRAL TRANSFORMS AND SPECIAL FUNCTIONS, 07 VOL. 8, NO. 9, 65 66 https://doi.org/0.080/06569.07.339039 Research Article On the derivatives / and / of the Legendre functions with respect to their degrees
More informationThe 3 dimensional Schrödinger Equation
Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum
More informationTwo special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p
LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p. 259-268. Fourier-Bessel and Fourier-Legendre series. p. 453-460. Boundary value problems in other coordinate system.
More informationPower series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0
Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can
More informationOpen problems about many-body Dirac operators
Open problems about many-body Dirac operators Jan Dereziński Dep. Math. Meth. in Phys., Faculty of Physics, University of Warsaw, Hoża 74, 00-682 Warszawa, Poland email: Jan.Derezinski@fuw.edu.pl October
More informationThe hydrogen molecule in magnetic fields: The ground states of the Σ manifold of the parallel configuration
The hydrogen molecule in magnetic fields: The ground states of the Σ manifold of the parallel configuration T.Detmer, P. Schmelcher, F. K. Diakonos and L. S. Cederbaum Theoretische Chemie, Physikalisch
More informationAdditional material: Linear Differential Equations
Chapter 5 Additional material: Linear Differential Equations 5.1 Introduction The material in this chapter is not formally part of the LTCC course. It is included for completeness as it contains proofs
More information2. The Schrödinger equation for one-particle problems. 5. Atoms and the periodic table of chemical elements
1 Historical introduction The Schrödinger equation for one-particle problems 3 Mathematical tools for quantum chemistry 4 The postulates of quantum mechanics 5 Atoms and the periodic table of chemical
More informationNuclear spin-spin coupling in HD, HT, and DT
Nuclear spin-spin coupling in HD, HT, and DT Mariusz Puchalski, Jacek Komasa, and Krzysztof Pachucki 2 Faculty of Chemistry, Adam Mickiewicz University, Umultowska 89b, 6-64 Poznań, Poland 2 Faculty of
More informationSymbolic Calculation of Two-Center Overlap Integrals Over Slater-Type Orbitals
Journal of the Chinese Chemical Society, 2004, 51, 243-252 243 Symbolic Calculation of Two-Center Overlap Integrals Over Slater-Type Orbitals Sedat Gümü and Telhat Özdoan Department of Physics, Amasya
More informationHigher-order C n dispersion coefficients for hydrogen
Higher-order C n dispersion coefficients for hydrogen J Mitroy* and M W J Bromley Faculty of Technology, Charles Darwin University, Darwin NT 0909, Australia Received 2 November 2004; published 11 March
More informationNew series expansions of the Gauss hypergeometric function
New series expansions of the Gauss hypergeometric function José L. López and Nico. M. Temme 2 Departamento de Matemática e Informática, Universidad Pública de Navarra, 36-Pamplona, Spain. e-mail: jl.lopez@unavarra.es.
More informationSolving ground eigenvalue and eigenfunction of spheroidal wave equation at low frequency by supersymmetric quantum mechanics method
Chin. Phys. B Vol. 0, No. (0) 00304 Solving ground eigenvalue eigenfunction of spheroidal wave equation at low frequency by supersymmetric quantum mechanics method Tang Wen-Lin( ) Tian Gui-Hua( ) School
More informationSeries solutions to a second order linear differential equation with regular singular points
Physics 6C Fall 0 Series solutions to a second order linear differential equation with regular singular points Consider the second-order linear differential equation, d y dx + p(x) dy x dx + q(x) y = 0,
More informationEvaluation of Coulomb and exchange integrals for higher excited states of helium atom by using spherical harmonics series
J Math Chem () 5:86 DOI.7/s9--9997-6 ORIGINAL PAPER Evaluation of Coulomb and exchange integrals for higher excited states of helium atom by using spherical harmonics series Artit Hutem Sutee Boonchui
More informationContents. I Basic Methods 13
Preface xiii 1 Introduction 1 I Basic Methods 13 2 Convergent and Divergent Series 15 2.1 Introduction... 15 2.1.1 Power series: First steps... 15 2.1.2 Further practical aspects... 17 2.2 Differential
More informationp-adic Feynman s path integrals
p-adic Feynman s path integrals G.S. Djordjević, B. Dragovich and Lj. Nešić Abstract The Feynman path integral method plays even more important role in p-adic and adelic quantum mechanics than in ordinary
More information1 The Quantum Anharmonic Oscillator
1 The Quantum Anharmonic Oscillator Perturbation theory based on Feynman diagrams can be used to calculate observables in Quantum Electrodynamics, like the anomalous magnetic moment of the electron, and
More informationQualification Exam: Mathematical Methods
Qualification Exam: Mathematical Methods Name:, QEID#41534189: August, 218 Qualification Exam QEID#41534189 2 1 Mathematical Methods I Problem 1. ID:MM-1-2 Solve the differential equation dy + y = sin
More informationWeek 5-6: Lectures The Charged Scalar Field
Notes for Phys. 610, 2011. These summaries are meant to be informal, and are subject to revision, elaboration and correction. They will be based on material covered in class, but may differ from it by
More informationStability conditions for hydrogen-antihydrogen-like quasimolecules
Purdue University Purdue e-pubs Birck and NCN Publications Birck Nanotechnology Center May 8 Stability conditions for hydrogen-antihydrogen-like quasimolecules Alejandro Ferrón Universidad Nacional de
More informationarxiv:quant-ph/ v5 10 Feb 2003
Quantum entanglement of identical particles Yu Shi Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, United Kingdom and Theory of
More informationarxiv: v1 [cond-mat.stat-mech] 14 Apr 2009
arxiv:0904.2126v1 [cond-mat.stat-mech] 14 Apr 2009 Critical exponents for Gaussian fixed point of renormalization Witold Haliniak 1 and Wojciech Wislicki 2 1 University of Warsaw, Faculty of Mathematics,
More informationMethod of Frobenius. General Considerations. L. Nielsen, Ph.D. Dierential Equations, Fall Department of Mathematics, Creighton University
Method of Frobenius General Considerations L. Nielsen, Ph.D. Department of Mathematics, Creighton University Dierential Equations, Fall 2008 Outline 1 The Dierential Equation and Assumptions 2 3 Main Theorem
More information5.4 Bessel s Equation. Bessel Functions
SEC 54 Bessel s Equation Bessel Functions J (x) 87 # with y dy>dt, etc, constant A, B, C, D, K, and t 5 HYPERGEOMETRIC ODE At B (t t )(t t ), t t, can be reduced to the hypergeometric equation with independent
More informationarxiv: v1 [hep-ph] 30 Dec 2015
June 3, 8 Derivation of functional equations for Feynman integrals from algebraic relations arxiv:5.94v [hep-ph] 3 Dec 5 O.V. Tarasov II. Institut für Theoretische Physik, Universität Hamburg, Luruper
More information37. f(t) sin 2t cos 2t 38. f(t) cos 2 t. 39. f(t) sin(4t 5) 40.
28 CHAPTER 7 THE LAPLACE TRANSFORM EXERCISES 7 In Problems 8 use Definition 7 to find {f(t)} 2 3 4 5 6 7 8 9 f (t),, f (t) 4,, f (t) t,, f (t) 2t,, f (t) sin t,, f (t), cos t, t t t 2 t 2 t t t t t t t
More informationarxiv:quant-ph/ v1 11 Nov 1995
Quantization of electromagnetic fields in a circular cylindrical cavity arxiv:quant-ph/950v Nov 995 K. Kakazu Department of Physics, University of the Ryukyus, Okinawa 903-0, Japan Y. S. Kim Department
More informationChem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer. Lecture 9, February 8, 2006
Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer Lecture 9, February 8, 2006 The Harmonic Oscillator Consider a diatomic molecule. Such a molecule
More informationA note on the Lipkin model in arbitrary fermion number
Prog. Theor. Exp. Phys. 017, 081D01 9 pages) DOI: 10.1093/ptep/ptx105 Letter A note on the Lipkin model in arbitrary fermion number Yasuhiko Tsue 1,,, Constança Providência 1,, João da Providência 1,,
More informationSecond Quantization Method for Bosons
Second Quantization Method for Bosons Hartree-Fock-based methods cannot describe the effects of the classical image potential (cf. fig. 1) because HF is a mean-field theory. DFF-LDA is not able either
More informationarxiv: v1 [math-ph] 15 May 2012
All solutions of the n = 5 Lane Emden equation Patryk Mach arxiv:05.3480v [math-ph] 5 May 0 M. Smoluchowski Institute of Physics, Jagiellonian University, Reymonta 4, 30-059 Kraków, Poland Abstract All
More informationIntroduction to Mathematical Physics
Introduction to Mathematical Physics Methods and Concepts Second Edition Chun Wa Wong Department of Physics and Astronomy University of California Los Angeles OXFORD UNIVERSITY PRESS Contents 1 Vectors
More informationEnergy spectrum for a short-range 1/r singular potential with a nonorbital barrier using the asymptotic iteration method
Energy spectrum for a short-range 1/r singular potential with a nonorbital barrier using the asymptotic iteration method A. J. Sous 1 and A. D. Alhaidari 1 Al-Quds Open University, Tulkarm, Palestine Saudi
More informationLecture 10. The Dirac equation. WS2010/11: Introduction to Nuclear and Particle Physics
Lecture 10 The Dirac equation WS2010/11: Introduction to Nuclear and Particle Physics The Dirac equation The Dirac equation is a relativistic quantum mechanical wave equation formulated by British physicist
More informationMA22S3 Summary Sheet: Ordinary Differential Equations
MA22S3 Summary Sheet: Ordinary Differential Equations December 14, 2017 Kreyszig s textbook is a suitable guide for this part of the module. Contents 1 Terminology 1 2 First order separable 2 2.1 Separable
More informationGround State Calculations of Confined Hydrogen Molecule H Using Variational Monte Carlo Method
Ground State Calculations of Confined Hydrogen Molecule H Using Variational Monte Carlo Method S. B. Doma 1), F. N. ElGammal 2) and A. A. Amer 2) 1) Mathematics Department, Faculty of Science, Alexandria
More informationLaplace Transform of Spherical Bessel Functions
Laplace Transform of Spherical Bessel Functions arxiv:math-ph/01000v 18 Jan 00 A. Ludu Department of Chemistry and Physics, Northwestern State University, Natchitoches, LA 71497 R. F. O Connell Department
More informationBessel Functions Michael Taylor. Lecture Notes for Math 524
Bessel Functions Michael Taylor Lecture Notes for Math 54 Contents 1. Introduction. Conversion to first order systems 3. The Bessel functions J ν 4. The Bessel functions Y ν 5. Relations between J ν and
More information221B Lecture Notes Many-Body Problems II Molecular Physics
1 Molecules 221B Lecture Notes Many-Body Problems II Molecular Physics In this lecture note, we discuss molecules. I cannot go into much details given I myself am not familiar enough with chemistry. But
More informationPolyexponentials. Khristo N. Boyadzhiev Ohio Northern University Departnment of Mathematics Ada, OH
Polyexponentials Khristo N. Boyadzhiev Ohio Northern University Departnment of Mathematics Ada, OH 45810 k-boyadzhiev@onu.edu 1. Introduction. The polylogarithmic function [15] (1.1) and the more general
More informationChem120a : Exam 3 (Chem Bio) Solutions
Chem10a : Exam 3 (Chem Bio) Solutions November 7, 006 Problem 1 This problem will basically involve us doing two Hückel calculations: one for the linear geometry, and one for the triangular geometry. We
More informationIntroductions to ExpIntegralEi
Introductions to ExpIntegralEi Introduction to the exponential integrals General The exponential-type integrals have a long history. After the early developments of differential calculus, mathematicians
More informationStructural Aspects of Numerical Loop Calculus
Structural Aspects of Numerical Loop Calculus Do we need it? What is it about? Can we handle it? Giampiero Passarino Dipartimento di Fisica Teorica, Università di Torino, Italy INFN, Sezione di Torino,
More informationHigher-order recoil corrections in Helium and Helium-like ions
Higher-order recoil corrections in Helium and Helium-like ions Vojtěch Patkóš Department of Chemical Physics and Optics Charles University in Prague 19th May 2017 Motivation Motivation Hydrogenic systems
More informationChemistry 3502/4502. Final Exam Part I. May 14, 2005
Advocacy chit Chemistry 350/450 Final Exam Part I May 4, 005. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle
More informationHighly accurate evaluation of the few-body auxiliary functions and four-body integrals
INSTITUTE OF PHYSICSPUBLISHING JOURNAL OFPHYSICSB: ATOMIC, MOLECULAR AND OPTICAL PHYSICS J. Phys. B: At. Mol. Opt. Phys. 36 (2003) 1857 1867 PII: S0953-4075(03)58918-0 Highly accurate evaluation of the
More informationDensity Functional Theory
Density Functional Theory March 26, 2009 ? DENSITY FUNCTIONAL THEORY is a method to successfully describe the behavior of atomic and molecular systems and is used for instance for: structural prediction
More informationd 1 µ 2 Θ = 0. (4.1) consider first the case of m = 0 where there is no azimuthal dependence on the angle φ.
4 Legendre Functions In order to investigate the solutions of Legendre s differential equation d ( µ ) dθ ] ] + l(l + ) m dµ dµ µ Θ = 0. (4.) consider first the case of m = 0 where there is no azimuthal
More informationConsidering our result for the sum and product of analytic functions, this means that for (a 0, a 1,..., a N ) C N+1, the polynomial.
Lecture 3 Usual complex functions MATH-GA 245.00 Complex Variables Polynomials. Construction f : z z is analytic on all of C since its real and imaginary parts satisfy the Cauchy-Riemann relations and
More informationLecture 10. Born-Oppenheimer approximation LCAO-MO application to H + The potential energy surface MOs for diatomic molecules. NC State University
Chemistry 431 Lecture 10 Diatomic molecules Born-Oppenheimer approximation LCAO-MO application to H + 2 The potential energy surface MOs for diatomic molecules NC State University Born-Oppenheimer approximation
More informationPhysics 139B Solutions to Homework Set 4 Fall 2009
Physics 139B Solutions to Homework Set 4 Fall 9 1. Liboff, problem 1.16 on page 594 595. Consider an atom whose electrons are L S coupled so that the good quantum numbers are j l s m j and eigenstates
More informationLaplace s Method for Ordinary Differential Equations
Physics 2400 Spring 2016 Laplace s Method for Ordinary Differential Equations Lecture notes by M. G. Rozman Last modified: March 31, 2016 We can [1, pp. 124-8], [2, h. VIII], [3, Apps. a-b], [4, h. 18],
More information4 Power Series Solutions: Frobenius Method
4 Power Series Solutions: Frobenius Method Now the ODE adventure takes us to series solutions for ODEs, a technique A & W, that is often viable, valuable and informative. These can be readily applied Sec.
More informationH 2 in the minimal basis
H 2 in the minimal basis Alston J. Misquitta Centre for Condensed Matter and Materials Physics Queen Mary, University of London January 27, 2016 Overview H 2 : The 1-electron basis. The two-electron basis
More informationProject: Vibration of Diatomic Molecules
Project: Vibration of Diatomic Molecules Objective: Find the vibration energies and the corresponding vibrational wavefunctions of diatomic molecules H 2 and I 2 using the Morse potential. equired Numerical
More informationMethods of Fermion Monte Carlo
Methods of Fermion Monte Carlo Malvin H. Kalos kalos@llnl.gov Institute for Nuclear Theory University of Washington Seattle, USA July, 2013 LLNL-PRES-XXXXXX This work was performed under the auspices of
More informationLecture X: External fields and generation of gravitational waves
Lecture X: External fields and generation of gravitational waves Christopher M. Hirata Caltech M/C 350-17, Pasadena CA 91125, USA (Dated: November 12, 2012) I. OVEVIEW Having examined weak field gravity
More informationExact solutions of the radial Schrödinger equation for some physical potentials
arxiv:quant-ph/070141v1 14 Feb 007 Exact solutions of the radial Schrödinger equation for some physical potentials Sameer M. Ikhdair and Ramazan Sever Department of Physics, Near East University, Nicosia,
More informationHartree-Fock-Roothan Self-Consistent Field Method
Hartree-Fock-Roothan Self-Consistent Field Method 1. Helium Here is a summary of the derivation of the Hartree-Fock equations presented in class. First consider the ground state of He and start with with
More informationBy Soroku TOH. (Read Dec. 16, 1939 by M. KOTANI).
Quantum-Mechanical Treatment of Helium-Hydride Molecule Ion HeH+. By Soroku TOH. (Read Dec. 16, 1939 by M. KOTANI). 1. Helium-hydride molecule ion HeH+, whose existence and process of formation was investigated
More informationPauli Deformation APPENDIX Y
APPENDIX Y Two molecules, when isolated say at infinite distance, are independent and the wave function of the total system might be taken as a product of the wave functions for the individual molecules.
More informationUniversal Associated Legendre Polynomials and Some Useful Definite Integrals
Commun. Theor. Phys. 66 0) 158 Vol. 66, No., August 1, 0 Universal Associated Legendre Polynomials and Some Useful Definite Integrals Chang-Yuan Chen í ), 1, Yuan You ), 1 Fa-Lin Lu öß ), 1 Dong-Sheng
More informationRelevant sections from AMATH 351 Course Notes (Wainwright): Relevant sections from AMATH 351 Course Notes (Poulin and Ingalls):
Lecture 5 Series solutions to DEs Relevant sections from AMATH 35 Course Notes (Wainwright):.4. Relevant sections from AMATH 35 Course Notes (Poulin and Ingalls): 2.-2.3 As mentioned earlier in this course,
More informationBäcklund transformation and special solutions for Drinfeld Sokolov Satsuma Hirota system of coupled equations
arxiv:nlin/0102001v1 [nlin.si] 1 Feb 2001 Bäcklund transformation and special solutions for Drinfeld Sokolov Satsuma Hirota system of coupled equations Ayşe Karasu (Kalkanli) and S Yu Sakovich Department
More informationIsotropic harmonic oscillator
Isotropic harmonic oscillator 1 Isotropic harmonic oscillator The hamiltonian of the isotropic harmonic oscillator is H = h m + 1 mω r (1) = [ h d m dρ + 1 ] m ω ρ, () ρ=x,y,z a sum of three one-dimensional
More information2 Electronic structure theory
Electronic structure theory. Generalities.. Born-Oppenheimer approximation revisited In Sec..3 (lecture 3) the Born-Oppenheimer approximation was introduced (see also, for instance, [Tannor.]). We are
More informationHIROSHI NAKATSUJI, HIROYUKI NAKASHIMA
How Does the Free Complement Wave Function Become Accurate and Exact Finally for the Hydrogen Atom Starting From the Slater and Gaussian Initial Functions and for the Helium Atom on the Cusp Conditions?
More informationThe Sommerfeld Polynomial Method: Harmonic Oscillator Example
Chemistry 460 Fall 2017 Dr. Jean M. Standard October 2, 2017 The Sommerfeld Polynomial Method: Harmonic Oscillator Example Scaling the Harmonic Oscillator Equation Recall the basic definitions of the harmonic
More informationExplicit solution of a class of quartic Thue equations
ACTA ARITHMETICA LXIV.3 (1993) Explicit solution of a class of quartic Thue equations by Nikos Tzanakis (Iraklion) 1. Introduction. In this paper we deal with the efficient solution of a certain interesting
More informationwhich implies that we can take solutions which are simultaneous eigen functions of
Module 1 : Quantum Mechanics Chapter 6 : Quantum mechanics in 3-D Quantum mechanics in 3-D For most physical systems, the dynamics is in 3-D. The solutions to the general 3-d problem are quite complicated,
More informationCHAPTER 11 MOLECULAR ORBITAL THEORY
CHAPTER 11 MOLECULAR ORBITAL THEORY Molecular orbital theory is a conceptual extension of the orbital model, which was so successfully applied to atomic structure. As was once playfuly remarked, a molecue
More informationTwo-Body Problem. Central Potential. 1D Motion
Two-Body Problem. Central Potential. D Motion The simplest non-trivial dynamical problem is the problem of two particles. The equations of motion read. m r = F 2, () We already know that the center of
More informationChemistry 593: The Semi-Classical Limit David Ronis McGill University
Chemistry 593: The Semi-Classical Limit David Ronis McGill University. The Semi-Classical Limit: Quantum Corrections Here we will work out what the leading order contribution to the canonical partition
More informationSOME CONGRUENCES ASSOCIATED WITH THE EQUATION X α = X β IN CERTAIN FINITE SEMIGROUPS
SOME CONGRUENCES ASSOCIATED WITH THE EQUATION X α = X β IN CERTAIN FINITE SEMIGROUPS THOMAS W. MÜLLER Abstract. Let H be a finite group, T n the symmetric semigroup of degree n, and let α, β be integers
More informationarxiv:physics/ v1 [physics.plasm-ph] 13 Jan 2005
PPPL 2467 (Aug. 1987 Phys. Rev. Lett. 59(16, 1817 1820 (Oct. 1987 DIFFERENTIAL FORM OF THE COLLISION INTEGRAL FOR A RELATIVISTIC PLASMA arxiv:physics/0501067v1 [physics.plasm-ph] 13 Jan 2005 Bastiaan J.
More informationPHY 407 QUANTUM MECHANICS Fall 05 Problem set 1 Due Sep
Problem set 1 Due Sep 15 2005 1. Let V be the set of all complex valued functions of a real variable θ, that are periodic with period 2π. That is u(θ + 2π) = u(θ), for all u V. (1) (i) Show that this V
More informationI. Perturbation Theory and the Problem of Degeneracy[?,?,?]
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Chemistry 5.76 Spring 19 THE VAN VLECK TRANSFORMATION IN PERTURBATION THEORY 1 Although frequently it is desirable to carry a perturbation treatment to second or third
More informationA new class of sum rules for products of Bessel functions. G. Bevilacqua, 1, a) V. Biancalana, 1 Y. Dancheva, 1 L. Moi, 1 2, b) and T.
A new class of sum rules for products of Bessel functions G. Bevilacqua,, a V. Biancalana, Y. Dancheva, L. Moi,, b and T. Mansour CNISM and Dipartimento di Fisica, Università di Siena, Via Roma 56, 5300
More informationCHEM3023: Spins, Atoms and Molecules
CHEM3023: Spins, Atoms and Molecules Lecture 3 The Born-Oppenheimer approximation C.-K. Skylaris Learning outcomes Separate molecular Hamiltonians to electronic and nuclear parts according to the Born-Oppenheimer
More informationarxiv: v2 [math.ca] 2 Sep 2017
A note on the asymptotics of the modified Bessel functions on the Stokes lines arxiv:1708.09656v2 [math.ca] 2 Sep 2017 R. B. Paris Division of Computing and Mathematics, University of Abertay Dundee, Dundee
More informationInstantaneous Multiphoton Ionization Rate and Initial Distribution of Electron Momenta. Abstract
Instantaneous Multiphoton Ionization Rate and Initial Distribution of Electron Momenta Denys I. Bondar Department of Physics and Astronomy, University of Waterloo, 00 University Avenue West, arxiv:0805.1890v
More informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More informationMOLECULAR STRUCTURE. The general molecular Schrödinger equation, apart from electron spin effects, is. nn ee en
MOLECULAR STRUCTURE The Born-Oppenheimer Approximation The general molecular Schrödinger equation, apart from electron spin effects, is ( ) + + V + V + V =E nn ee en T T ψ ψ n e where the operators in
More informationEnergy spectrum inverse problem of q-deformed harmonic oscillator and WBK approximation
Journal of Physics: Conference Series PAPER OPEN ACCESS Energy spectrum inverse problem of q-deformed harmonic oscillator and WBK approximation To cite this article: Nguyen Anh Sang et al 06 J. Phys.:
More informationHarmonic Oscillator Eigenvalues and Eigenfunctions
Chemistry 46 Fall 217 Dr. Jean M. Standard October 4, 217 Harmonic Oscillator Eigenvalues and Eigenfunctions The Quantum Mechanical Harmonic Oscillator The quantum mechanical harmonic oscillator in one
More informationA General Expression for Symmetry Factors of Feynman Diagrams. Abstract
A General Expression for Symmetry Factors of Feynman Diagrams C.D. Palmer a and M.E. Carrington b,c a Department of Mathematics, Brandon University, Brandon, Manitoba, R7A 6A9 Canada b Department of Physics,
More informationHelium fine structure Ingvar Lindgren (mod )
Helium fine structure Ingvar Lindgren 2004.09.20 (mod. 2005.01.23) The leading contributions to the helium fine structure beyond the first-order relativistic contribution were first derived by Araki [1]
More informationGeneral Exam Part II, Fall 1998 Quantum Mechanics Solutions
General Exam Part II, Fall 1998 Quantum Mechanics Solutions Leo C. Stein Problem 1 Consider a particle of charge q and mass m confined to the x-y plane and subject to a harmonic oscillator potential V
More informationBayero Journal of Pure and Applied Sciences, 3(1): Received: September, 2009 Accepted: May, 2010
Bajopas Volume 3 Number June Bayero Journal of Pure and Applied Sciences, 3(: - 7 Received: September, 9 Accepted: May, VARIAIONAL QUANUM MON CARLO CALCULAION OF H GROUND SA NRG OF HDROGN MOLCUL Suleiman,
More informationRegularization Physics 230A, Spring 2007, Hitoshi Murayama
Regularization Physics 3A, Spring 7, Hitoshi Murayama Introduction In quantum field theories, we encounter many apparent divergences. Of course all physical quantities are finite, and therefore divergences
More informationElectron States of Diatomic Molecules
IISER Pune March 2018 Hamiltonian for a Diatomic Molecule The hamiltonian for a diatomic molecule can be considered to be made up of three terms Ĥ = ˆT N + ˆT el + ˆV where ˆT N is the kinetic energy operator
More informationOn the Weyl symbol of the resolvent of the harmonic oscillator
On the Weyl symbol of the resolvent of the harmonic oscillator Jan Dereziński, Maciej Karczmarczyk Department of Mathematical Methods in Physics, Faculty of Physics University of Warsaw, Pasteura 5, -93,
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 21 Quantum Mechanics in Three Dimensions Lecture 21 Physics 342 Quantum Mechanics I Monday, March 22nd, 21 We are used to the temporal separation that gives, for example, the timeindependent
More information