Finite Element Method of high-order accuracy for solving two Dimensional Elliptic boundary-value Problems of two and three identical atoms in a line Текст научной статьи по специальности «Физика»

#

Вестник РУДН. Серия МИФ

RUDN Journal of MIPh

http://journals.rudn.ru/miph

2018 Vol. 26 No. 3 226-243

UDC 519.632.4

DOI: 10.22363/2312-9735-2018-26-3-226-243

Finite Element Method of High-Order Accuracy for Solving Two Dimensional Elliptic Boundary-Value Problems of Two and Three Identical Atoms in a Line A. A. Gusev

We considered models of three identical atoms in a line with molecular pair interactions and diatomic molecule scattered by an atom or tunneling through potential barriers. The models are formulated as 2D elliptic boundary-value problems (BVPs) in the Jacobi and polar coordinates. The BVP in Jacobi coordinates solved by finite element method of high-order accuracy for discrete spectrums of models under consideration. To solve the scattering problems the BVP in polar coordinates are reduced by means of Kantorovich method to a system of second-order ordinary differential equations with respect to the radial variable using the expansion of the desired solutions in the set of angular basis functions that depend on the radial variable as a parameter. The efficiency of the elaborated method, algorithms and programs is demonstrated by benchmark calculations of the resonance scattering, metastable and bound states of the considered models and also by a comparison of results for bound states of the three atomic system in the framework of direct solving the BVP by FEM and Kantorovich reduction.

Key words and phrases: Elliptic Boundary-Value Problems, scattering problem, metastable and bound states, Kantorovich method

With the development of modern computing power, there are more possibilities for setting and numerically solving multidimensional boundary-value problems with high accuracy. For this, new numerical methods of high accuracy order are being developed. When reducing the boundary value problem to an algebraic one in the finite element method (FEM) of order p, one of the problems is the calculation of integrals on a finite element (we consider only simplicial finite elements) containing the products of two basis functions of Lagrange or Hermite interpolation polynomials of order p by the coefficients for the unknown functions [1-3]. The FEM schemes and algorithms for solving the d-dimensional BVP are presented in [4-6]. The algorithms for constructing 2-dimensional fully symmetric Gaussian quadratures are presented in [7] while d-dimensional ones will be published elsewhere.

The aim of this paper is to demonstrate the efficiency of the elaborated algorithms and program complexes KANTBP 4M, KANTBP 3, ODPEVP [8-10] by benchmark calculations of the resonance scattering below the dissociation threshold, metastable and bound states of the considered models and also by a comparison of results for bound states of the three atomic system in the framework of direct solving BVP by 2D FEM [3,7] and Kantorovich method (KM).

We elaborated also algorithms for calculating the asymptotic parametric angular functions, the effective potentials and the fundamental solutions of the SODEs in the form of expansions by inverse powers of radial variable [11] and apply them to the construction of the asymptotic states of the triatomic scattering problem, because in three-atomic systems at large values of the hyperradial variable the effective potentials of SODEs have

Received 27th May, 2018.

The author thank Profs. O. Chuluunbaatar, V. L. Derbov, V. P. Gerdt, S. I. Vinitsky and Dr. L. L. Hai for collaboration. The work was supported by the RFBR (grants No. 16-01-00080 and 18-51-18005).

Laboratory of Information Technologies Joint Institute for Nuclear Research 6, Joliot-Curie str., Dubna, Moscow region, Russia, 141980

1. Introduction

the form of expansions in inverse powers of hyperradius even for short-range potentials of pair interactions.

The paper is organized as follows. In Section 2 we formulate the 2D BVP for the dimer and trimer models. In Section 3 the 2D BVP in polar coordinates is reduced to the SODEs. In Sections 4, 5 and 6 we demonstrate the efficiency of the elaborated technique by test calculations of the resonance scattering of a diatomic molecule on an atom, the metastable and bound states of three atomic system, and the resonance tunneling of a diatomic molecule through a Gaussian barrier and the metastable states, respectively. In Conclusion the results and perspectives are discussed.

2. Setting of the Problem

Consider a model of three identical particles (trimer) on a straight line with the masses Mi = M and the coordinates xi G R1, i = 1,2, 3, of one-dimensional Euclidian space R1 coupled via the pair short-range potential T^(|Xi — Xj|), i,j = 1, 2, 3. Performing the change of variables corresponding to the cyclic permutation ) = (1, 2, 3) of

the permutation group S3 [12]:

_ __X q I X R 2cc-y ._ _ _ ,

Xj = X(aß) =Xß -xa, y7 = y(aß)j =--L, Xo = +X2 + X3),

T

we get three pairs (X(aß), y(aß)1) of the scaled Jacobi variables. In the center-of-mass reference frame of the configuration space x = (X1,X2,X3) G R3, i.e., on the hyperplane Q = {x = (X1,X2,c3) |X1 + X2 + X3 = 0} C R3 these Jacobi maps (x, y)T G Qxy ~ R2 are connected by the relations:

/Vj i X{aß)\ = (Xa

\y-rj \y(aßhj yy '

T ( ) = / cospia sin pla\ xy (lja) = sin pia cos ipla)

Va/ (1)

where the angles p-a = pap = pp-- = -pa- = -ppa = -p--p = (2^)/3, i.e. cos p-a = cos pap = cos = -1/2, sin p-a = sin pap = sin = \/3/2, are such that p-a + pap + Pp-( = 2^, pa- + ppa + p7p = -2^. With the variable parameter p in Txy(p) they and simply change the sign of pairs of Jacobi coordinates in Eq. (1),

( -X(a) = T y (p-a ±«) (X(P-r) ) . (2)

V-y(aPh) XJ \y(p-r)aj V J

In Fig. 1 (a) the hyperplane (x, y)T £ QXy is separated by three axes y(12)3, y(31)2 and V(23)i into six sectors that together with three orthogonal axes x12, x31, x23 describe six channels of the scattering problem for three identical particles. Our choice is determined by the parameterization of the pairs (y,x)t = (p cos p,p sinp)T, where the angle p £ [0, 2^) is counted counterclockwise from the axis y^12)3, for which p = 0.

As a result, we arrive at the Schrodinger equation for the wave function ^(y,x) corresponding to the total energy E of the triatomic system in Jacobi coordinates (x, y)T = (x-, y-)T £ vXy:

d2 d2 M ~ ~ \

- W2 - JX2 + M(E<X- rt - E)J y(y>x) = 0. (3)

(5(3X2)®©© (J®©

® @®\y\ k

y(i2)3'©@®

@©@

©0©'y,/21 \ ^ ©@® @@® ©@®

200

(a)

©A®

200-

100

JV®@ y JV@®

(I

-100

(b)

Figure 1. (a) The profiles of 2D potential functions of Be3 in Jacobi coordinates

(1) in the center-of-mass plane and the relative arrangement of particles in accordance with the region of the center-of-mass plane. The numbers of sectors are given in boxes. (b) The profiles of 2D potential functions of Be2 with barrier and the relative arrangement of particles and barrier. Here coordinates and potential functions are given in A and A-2, respectively

In the case of a diatomic molecule (dimer) with identical nuclei coupled via the pair potential, V(|xi — x2|), moving in the external potential field V6(|Xi — X31), i = 2,1, of the third atom having the infinite mass, the same equation (3) is valid for the variables

X = X3 = X2 — Xi, y = y3 = xi + X2,

with respect to the origin of the coordinate frame placed at the infinite-mass atom, X3 = 0.

Here the potential function for the trimer with the pair potentials (below this case is referred to as Task 2, for example, see Fig. 1 (a)),

V(X, y) = V(|x|) + V

i x — V3y \ v ( x + V3y \

V 2 ) V 2 )

(4)

or the potential function for a dimer in the field of a barrier potential (below this case is referred to as Task 3, see Fig. 1 (b))

V(X, y) = V(|x|) + v6

( | ^IH '(14*1 )•

(5)

is an even function with respect to the straight line X = 0 (i.e., x1 = x2), which allows one to consider the solutions of the problem in the half-plane X ^ 0. The use of the Dirichlet or Neumann boundary condition at x = 0 allows one to obtain the solutions, symmetric and antisymmetric with respect to the permutation of two particles. If the pair potential possesses a high peak at the pair collision point, then the solution of the problem in the vicinity of x = 0 is exponentially small and can be considered in the half-plane X ^ xm;n. In this case the Neumann or Dirichlet boundary condition imposed at Xmin yields only a minor contribution to the solution. The equation, describing the diatomic molecular subsystem (dimer) (below this case is referred to as Task 1), has the form

i d2 M ~ \

(—^ + f(v(|x|) — e)j 0(x)=0.

We assume the energy spectrum of the dimer to consist of the discrete part with a finite number n0^ 1 of bound states with the eigenfunctions <j(x), j = 1, no and the eigenvalues V = — |Vj |< 0, and the continuous part with the eigenvalues V > 0 and the corresponding eigenfunctions <i(x). In certain cases the eigenfunctions of the continuous spectrum are approximated by the eigenfunctions of pseudostates of the discrete spectrum ij > 0, j = 1 + no,... calculated in the sufficiently large but finite interval X £ [Xmin, Xmax].

The approach proposed below is illustrated by the example of Be2 with the reduced mass M/2 = 4.506 Da of the nuclei [13] and the molecular interaction approximated by the Morse potential

M ~ ~ ~ V (x) = -^V(x), V(x) = D {exp[-2(x — Xeq)a] — 2exp[-(x — Xeq)a]}. (7)

Here a = 2.96812 A-1 is the potential well width, Xeq = 2.47 A is the average distance between the nuclei, and DD = 1280 K, D = (M/h2)D = 236.51 A-2 (1 K=0.184766 A-2, 1 A-2=5.412262 K) is the potential well depth. This potential supports only five bound states corresponding to the covalent bounding of the Be2 dimer [14] having the energies ei = (M/h2) Vi, i = 1,...,n0 = 5. The parameter values are determined from the condition (V — V) / (2nhc) = 277.124 cm-1, 1 K/(2^h c)=0.69503476 cm-1. To solve the discrete spectrum problem we applied the seventh-order FEM using the Hermitian interpolation polynomials with double nodes [8]: £i = {—193.066, —119.392, —63.338, —24.904, —4.089}A-2, i = 1,...,no = 5. As an example, following [13], we use below the Gaussian barrier potentials Vb(xi) = D exp (—x2/a)

with D = 236.51 A-2 and a = 0.0523 A2. The values of parameters of the repulsive Gaussian barrier potential were estimated following the experimental observation of the quantum diffusion of hydrogen atoms on the copper surface [15].

3. Boundary-Value Problems

Using the change of variables x = p sinp, y = p cosp, we rewrite Eq. (3) in polar coordinates (p, p), Qp,p = (p £ (0, <x>),<p £ (0, 2^))

/18 8 1 \ d2

(—-pdpPdp + J2 p) — E) ^ p) = 0 p) = — + P2V(^ p), (8)

where for a trimer with pair potentials

V(p, p) = V(p| sin + V(p| sin(p — 2n/3)l) + V(p| sin(p — 4^/3)|), (9)

and for a dimer with pair potential in the external field of barrier potential:

V(p, p) = V(p| sin p|) + Vb(p| sin(p — ^/4)|) + Vb(p| sin(p + ^/4)|). (10)

where Vb(p| sin(p ±^/4)|) = D exp (—p2 sin2(p ±^/4)/a).

The solution of Eq. (8) is sought in the form of the Kantorovich expansion [16]

j max

^i o (p,p) = Y1 <j (p;p)Xji o (p). (11)

j=1

Here Xjio (p) are unknown matrix functions, j = 1,..., jmax = 2N. The angular basis functions <j(p;p) £ Tp ~ L2 (Q) in the interval Q = p £ [0,2^), which is divided into smax

subintervals QS = £ (2w(s — 1)/smax, 2ks/smax): = (J5=ai QS, are determined at each value of the parameter p £ (0, as the eigenfunctions corresponding to the real discrete eigenvalues £i(p) < e2(p) < • • • < £j (p) < • • • of the Sturm-Liouville problem for the equation

(A(p,p) — £j (p))0j (^;p) = 0, (12)

The functions p) = 0jp) have the parity (—1)p, p = 0,1 with respect to the

inversion of Jacobi coordinates (2), i.e., p) = (—1)p^jp)(^ ± w; p); here and below

= 0. These functions satisfy the orthogonality and completeness conditions

/ Smax ^max f

) = ^ / (^;p) = ^j^ (13)

max

= = ^ — ^ (14)

For the three problems under consideration the potential function V(^,p) depending on the parameter p can be defined as follows.

Task 1. The case of one pair potential in the intervals ^ £ (0, 2^a) = w/6, w/4 or w/2) V(^,p) = V(psin^).

Task 2. The case of three pair potentials, Eq. (9), in the interval ^ £ (0, = w/3); the potential has the symmetry of the dihedral point group [12].

Task 3. The case of one pair potential and two penetrable or almost impenetrable barriers, Eq. (10), in the interval ^ £ (0, = w/2) or in the intervals ^ £ (0, = w/4 — e) and ^ £ = w/4 — e,w/2), 0 < e < w/4.

The solutions symmetric with respect to the permutation of two particles satisfy the Neumann boundary condition at ^ = 0 and ^ = 2^a, while the antisymmetric ones satisfy the Dirichlet boundary condition at these points, i.e., at the aforementioned points X = Xa/3 = 0. If the pair potential possesses a high peak in the vicinity of the pair collision point, then the solution of the problem (8) will be considered in the half-plane Qv,p = (p £ (pmin, ro),^ £ [^min(p), — ^min(p)]) with the Neumann or Dirichlet boundary condition.

The potential function V(^,p) of the boundary-value problem (12) is symmetric under the reflection /a: ^ ^ (4s — 2)^a — ^ with respect to the lines ^ = (2s — 1)^a in each sector of the cycle, numbered by s = s max: /aV(^,p) = V(^,p). There-

fore, the set of eigenfunctions 0jp) is separated into two subsets, namely, the even and odd ^J=±1(^;p) ones: (<^;p) = ((4s — 2)^a — <^;p) = ±^J=±1(^;p). This fact allows separate calculation of gerade p) = ((4s — 2)^a — ^;p) or

ungerade = —((4s — 2)^a — ^;p) eigenfunctions in the reduced interval

^ £ [^min(p), subjecting them to Neumann or Dirichlet boundary condition at the boundary point ^ = of the interval, respectively. Below the parametric angular basis functions 0J=±1(^; p) with the numbers j = 1, • • •, n° are referred to as cluster states with £j (p) < 0, and those with j > n° + 1 as pseudostates with Sj (p) > 0, corresponding to the discrete and continuous spectrum of BVP for Eq. (6) at large values of the parameter p, respectively. To reveal the above structurization property, we introduce the linear

combinations of these functions ф]'3р) for the trimer and фГ'^(^; р) for the dimer:

ф]'3(р;р) = (±ф?=-1(^;р) + ф?=1(^;р)) /V2, (15)

фГ'^;р) = (±ф?=-1(^;р)+ф?=1(^;р))/V2.

The action of the parity operator Ia on these functions

1аф] (^;р) = ф3 (^;p), 1аф3 (^;р) = ф] (^;p), (iß)

consists in the permutation of states I7) о |ß) or | о | with respect to the lines ^ = (2s — Indeed, the parametric cluster functions at j = 1,..., n0 and at

large р have maximal in the vicinity of ^ = (2s — 2)^a and ^ = 2s^a, respectively, that correspond to the eigenfunctions of cluster states of the BVP for Eq. (6). In the particular case of sector 1, the dimer functions фГг'^((р; р) have maximal in the vicinity of ^ = 0 and ^ = i.e. for у > 0, v or for у < 0, v at large р, respectively while the trimer functions ф]'3р) have maximal in the vicinity of ^ = 0 and ^ = ^/3, i.e. for Ж12/У(12)3 > 0, v = 7, or for Ж31/2/(31)2 < 0t) = ^ at large р, respectively.

Using the above pairs of the basis functions (15), we rewrite the expansion (11) in the ^-representation

N

p) = EE ^ (^ p)*k (p) , (17)

where a = «,(/ or a = 7,^ for the unknown functions x]i0 (p) and (p) in the (7^)-representation, related to the functions (p) and (p) in the ((/«)-representation as

( ) Aï. (pA .(x'k (p)\ . win

^Uwj U.<rij • 4 = 72 l-i 0 ■ (i8)

The averaging of Eq. (8) with the basis functions in ^-representation (17) yields the system of coupled ODEs with determined by the hyperspherical parameterization of the two-dimensional configuration space

1 d p d + £*(P) F

--~Ä~P ~Ä~ +--2--F

p dp dp p2

X«0 (p) + (p)^ (p) = 0, (19)

3=1

Wij (p) = Hji (p) + --rpQji (p) + Qji (p) —. (20)

1 d _ , , ^ / n d

- TpPQ3^(P) + Q3^(P)Tp

Here the potential curves (terms) Sj (p) are eigenvalues of the BVP (12) and the effective potentials (EPs) Qij(p) = — Qji(p), Hij(p) = Hji(p) are expressed as integrals calculated in the reduced intervals ^ G [0, 2ipa] using the above u symmetry:

Qij (P) = )• Hij (p) = {dp^dp^j )■ (21)

For Task 3 the effective potentials Wij (p) = Wij (p) + Vij (p) are sums of Wij (p), calculated using the potential curves and the parametric basis functions of Task 1, and

the matrix elements V^ (p) of the barrier potentials

Vj(p) = <^|V6(p| sin(p - n/4)|) + V\pl sin(p + V4)|)|&)• (22)

As an example, we calculated with the required accuracy the parametric basis functions of BVP (12) and the effective potentials (21) for the models of Be2 dimer and Be3 trimer in collinear configuration using the FEM implemented in the program ODPEVP [10]. The results of calculation on the grid [1.8/p, pa] = {1.8/p(24)3/p(10)4/p(5)5/p(10)pa} for ipa = n/2 for Be2 dimer and = n/6 for Be3 trimer. Using the obtained result in the uncoupled (gM)-representation and the transformation matrix A from Eq. (18), we can rewrite the system of ODEs in the coupled (aft^representation with the effective potentials.

UM = (lrrr'(pl U?^ = A (Uiuj^(p) rT0f ^ A-1. (23)

\UiHrt (p) Ui3j3 (p)j \ 0 Ui9j9 (p)

In Section 4 one can see that the (7^)-representation provides the required compatibility of the solutions of Eqs. (19) with the asymptotic boundary conditions of the scattering problem on the full axis and its half-axis.

4. Asymptotic Expressions of Scattering, Metastable

and Bound States

The general solution Fj of the system of ODEs in the open channels iQ = 1,..., N0 is determined by a linear combination of the fundamental solutions x- (piap) and x%' (piop) calculated using Eqs. (19)-(21) following from [11] with the leading terms of the Hankel functions of the first and the second kind [17] H<(/^2(piop) and H(/2(piop) below the dissociation threshold at E < 0 in the form of incoming and outgoing waves

N

FjlS(P) = £ [x-(Pi„PK, + X+ (PioP)a+.

i'o = 1

at Pi0 = \/E — £?o > 0 in the open channels Na = m^xj ^ N. The scattering matrix

S?oii (E) or S?oii (p), where p = diag{p^o}^=1, is a diagonal matrix. In open channels it is defined as the matrix transforming the amplitudes of the incoming waves a- into those of the outgoing waves a+ [18]

N„

= (PK. (24)

i'o = l

The components of the radial asymptotic solutions Fj^ (p) of the scattering problem in the open channels i0 = 1,..., N0 have the form

N

FjH(p) = £ [x-(Pi.P)Sji'B + X+o(Pi0p)S?0i0(p)] , (25)

i'=l

while in the closed channel i0 = N0 + 1,...,N the asymptotic solutions Fjio(p) are determined by the fundamental solutions (ipip) with the leading term of Kelvin functions [17] ^1/2(Qiop) for the decaying waves

Ft (p) = 4*+, (^op). (26)

These asymptotic solutions F(p) = {Fio(p)}^=1 = {{Fjia(p)}j'L1}^=1 are used to have the conventional asymptotic boundary conditions for the components of the numerical solution Fjio (p) of the system of ODEs (19) at large p = pmax

F*(p) = Ft(p). dpF'*(p) = (p). (27)

The scattering problem (3)-(6) with the asymptotic boundary conditions (25) and (26) is reduced to a boundary-value problem for the set of close-coupled equations (19) with the boundary conditions at p = pmin and p = pm

Jmax-

F (Pmin) = 0, ^

dp

= R(pmax)F(Pmax) , (28)

p—Pmax

where ^.(pmax) is a N x N symmetric matrix function of E, F(p) = {Xia (p)}i^=i = {{Xjia(p)}j=1}il=1 is the required N x N numerical matrix solution.

These matrices and the N0 x N0 matrices S±1 = {Sf,'™ (E)}^°io=1 sought for in the open channels N0 = max j ^ N are calculated directly from (27) using the program

KANTBP 4M [8].

For metastable states the even and odd eigenfunctions obey the boundary conditions of the third kind (28), where the matrices ^(pmax) = diag(^(pmax)) depend on the complex energy eigenvalue E = Em = ReEm + «ImETO, ImETO < 0 sought for, and are expressed as [19]

, ipm, ReEm > ej | I-— -

Ki0i0 (pmax) = j ^, Re Em < \ ' Pm = yEm - £io' = 1 \J£t - E™ (29)

since the asymptotic solutions of this problem contain only outgoing waves in the open channels i0 = 1,..., N0 and closed ones i0 = N0 +1,..., N. In this case the eigenfunctions obey the orthogonality and normalization conditions

(Fm|Fm' ) (+ ^Pm')

pmax

y F^(p)F^(p)dp - S„

pmin

+ Cmm' — 0,

(30)

Cmm' FrTi(pmax)Fm/ (pmax)-

For bound states the even (g) and odd (w) eigenfunctions obey the boundary conditions (28), with the matrices K(^max) = diag(^(p max)) — 0. In this case the eigenfunctions obey the orthogonality and normalization conditions

Pmax

<Fro|Fro, > —J F^n(p)Fm, (p)dp — 5mm,. (31)

p

Taking the property (15) and (16) of the quasiangular parametric basis functions and the effective potentials (23) into account, we express the S-matrix in the (7 ft) representation on the full axis R1 via the matrix Sa (24), (25) calculated on the half-axis R+. The

matrix S is a unitary and symmetric scattering matrix

S = !t3>) , StS = SSt = I, (32)

VS37 S33j

consisting of the matrices Stt, S33 and ST3 = S^ of the dimension N0 x N0 determined by the relations Stt = S33 = (S+1 + S-1)/2, ST3 = SjT = (S+1 — S-1)/2, where S+1 = S g and S-1 = Su are the matrices from (25). Here I is the unit matrix with the dimension 2N0 x 2N0, Stt and S33 corresponds to the elastic scattering processes (in the considered case of 1D scattering it means reflection) of the dimer (a3)(or (7a)) on the atom 7 (or ft): 7 + (aft) ^7 + (aft), or (7a) + ft ^ (7a) + ft, and S3T and the matrices ST3 correspond to the inelastic rearrangement scattering processes (in the case of 1D scattering it means transmission) 7 + (aft) ^ (7a) + ft or (7a) + ft ^ 7 + (aft), for which the conventional relations between inelastic and elastic scattering below breakup threshold at E < 0 follow from (32) that provide conservation of the Wronskian, | S3T |2 = I — |Stt |2 and S3T = S^3, where I is the unit N0 x N0 matrix.

For the scattering of the dimer (7ft) by the potential barriers, considered on the full axis, the matrix S is the 2N0 x 2N0 scattering matrix (32) read as similar to [19]

S = (t R^) , StS = SSt = I, (33)

where I is the unit matrix with the dimension 2N0 x 2N0 consisting of the amplitudes of the reflected and transmitted waves Rv = Rv(E) and Tv = Tv(E), where v ^ indicates the direction of the incident wave propagation with respect to the y-axis, i.e., v and v for y > 0 and y < 0, respectively, and the N0 x N0 matrices R« = Rv(E) and Tv = R(E) are expressed as R^ = R^ = (S+i + S_i)/2, T^ = T^ = (S+i — S_i)/2.

For the scattering of the dimer (aft) on potential barriers similar relations determine the reflection R^ = Stt and R^ = S33, and transmission T^ = S3T and T^ = ST3 amplitudes. For the reflection coefficient |Rv|2 = R«Rv and the transmission coefficient |Tv|2 = TVTv the conventional relation below breakup threshold at E < 0 following from (32) and constant Wronskian, |Tv |2 = I — |Rv |2 is valid, where I is the unit N0 x N0 matrix.

5. Bound, Metastable and Scattering States of the trimer

For the considered models, the eigenvalues and the hyperradial components of 2D eigenfunctions of the BVP for the set of ODEs (19) with Dirichlet boundary conditions were calculated with the predetermined accuracy using the FEM implemented in the KANTBP 4M program [8].

The set of even (g) and odd (u) bound states of the trimer Be3 (Task 2) were calculated on the grid Qh = {pmin = 4.24(1)4.33(10)6.13(1)6.33(23)pmax = 11.39}, where in parentheses the number of fifth-order Hermitian elements [20] is indicated, for the number of equations N = 15 in the system (19).

The comparison of sets of total and binding energies of g and u bound states of the trimer Be3 calculated by KM and 2D FEM is presented in Table 1 and the corresponding eigenfunctions (11) are shown in Fig. 2. One can see that results obtained by KM and 2D FEM [3,7] are in agreement with an accuracy of the order 0.1A_2.

Table 1

The total Ei(*) and binding energies -El^(*)=Ei(*)—EtIsh of g and u bound states of the trimer Bes relative to the exact threshold Etlsh= - 193.06 (in A-2): Ei(K) and -Ebi(K) calculated by KM in polar coordinates using N=12 basis functions (g or u); Ei(F) and -Ebi(F) calculated by the 2D FEM using fifth order interpolation Lagrange polynomials on triangular finite elements in sector 1 of Jacobi

coordinates

i 1g 2g 3u 4g 5u 6g

-Ei(K) -EbdK) ) -Eb{F) 389.08 196.02 389.09 196.03 335.43 142.37 335.45 142.39 300.58 107.52 300.60 107.54 287.02 93.96 287.05 93.99 260.47 67.41 260.50 67.44 245.84 52.78 245.88 52.82

i 7u 8g 9u 10g 11u 12g

-Ei(K) -Eb(K) -Ei(F) -Eb(F) 227.66 34.60 227.70 34.64 225.39 32.33 225.42 32.36 215.37 22.31 215.41 22.35 204.85 11.79 204.89 11.83 198.21 5.15 198.24 5.18 193.86 0.80 193.91 0.85

45678 45678 45678 45678

y y y y

Figure 2. The density plots of the eigenfunctions tyf,'u(<p,p) displayed in sector 1 of (y, x)-plane (in A) of the gerade (g) and ungerade (u) bound states with energies of the Be3 trimer presented in Table 1. The negative, positive and near-zero values of the eigenfunctions are displayed by black, white and gray, respectively

These solutions, i.e. the real eigenvalues and the corresponding eigenfunctions, were used as an initial approximation in the continuous analogue of Newton's method [21] with additional condition (Fm|Fm/ ) = 0 to calculate the metastable states of the trimer Be3 on the same finite-element grid. The corresponding problem of the dimer scattering on an atom was solved on the same grid and N = 15.

The calculated complex energies of the metastable states E^ = Em = ReEm + «ImETO for the trimer are presented in Table 2.

Table 2

The sets of the first resonance energy values E at which the minimum of the transmission coefficient |£ab|2i is achieved, the number i of the threshold £j, the real and imaginary part of the complex energy eigenvalues Ev = Re Ef +1 Im Ef in A-2 of the even g and odd u metastable states of Be3 numbered by the index v calculated with N =15 equations (19)

E i \Sab\ii ReEf ImE^ ^ype

-193.066 1 thr

-189.676 1 1■10-7 -188.94 -4 10-2 1g

-164.654 1 4■10-6 -164.72 -1 10-2 1u

-156.882 1 3■10-6 -157.04 -2 10-2 2g

-140.545 1 1■10-4 -140.57 -5 10-3 2u

-132.485 1 1■10-6 -132.47 -4 10-3 3g

-124.256 1 1■10-3 -124.16 -6 10-3 3u

-120.638 1 1■10-6 -120.75 -8 10-2 4g

-119.392 2 thr

-113.248 2 0.10 -113.24 -2 10-2 5g

-89.319 2 9■10-4 -89.16 -3 10-2 6g

-77.271 2 0.77 -76.51 -4 10-6 4u

-70.309 2 0.35 -70.30 -2 10-3 7g

-63.385 2 0.41 -65.14 -3 10-4 5u

-63.338 3 thr

-42.858 3 0.06 -42.87 -6 10-3 8g

-29.396 3 0.13 -29.19 -4 10-2 9g

-24.899 3 0.19 -25.82 -1 10-3 6u

-24.904 4 thr

-6.799 4 0.40 -7.12 -1 10-3 10g

-4.089 5 thr

0 thr

These metastable states are responsible for resonance energies, corresponding to the minimal probability of inelastic scattering of the dimer by the atom, i.e., to the resonance quantum reflection from the potential well (Feshbach resonances, see Figs. 1 (a) and 3 (a).

As an example, in Fig. 4 we display the eigenfunctions of the scattering problem for gerade and ungerade states corresponding to the minimum of the transmission coefficient 1^7/312 = I^W|2 = -10-7 at E = -189.676, as well as the metastable state 1g from Table 2.

The isolines of the absolute values ( y, of the linear combinations (y, x) = (y± ^u(y,x))^\f2 demonstrates the effect of resonance reflection from the effective

Figure 3. The effective diagonal potentials Wjj (p) = £j (p)p-2 + Hjj (p) for the Be3 trimer (a) and effective diagonal potentials Wjj (p) = £j (p)p-2 + Hjj (p) + Vjj (p) for the tunneling problem of the dimer Be2 through the Gaussian barrier (b)

, | , | , | , | , | , | 4 5 6 7 8 9 10 456789 10

y y

Figure 4. Upper panel: the isolines of the absolute values (ip, p)\ of corresponding gerade (left-hand panel) and ungerade (right-hand panel) solutions in sector 1 of (y, x)-plane (in A) for the scattering of Be atom with the energy —E = 189.676 A-2 (relative to the three-body threshold) on the dimer Be2, corresponding to the 1g metastable state from Table 2. Lower panel: isolines of the absolute values (<p,p)\ of the linear combinations

^ (v, p) = (±^(ip,p) + ^g(<P,P))/V2

potential well. It can be seen from the figures that the shape of the wave functions of the gerade scattering states (Fig. 4 (a)) and metastable states (Fig. 5) are similar and they are localized in the vicinity of the potential well (Fig. 3 (a)). At the same time, for the same energy value E = -189.6 A-2, the wave function of ungerade scattering states (Fig. 4 b) is a typical nonresonant wave function.

Figure 5. The components \j and isolines of the absolute values || of the solution for the trimer Be3 in 1g metastable state with the real part of

energy eigenvalues Re E = -188.94 A-2, localized near the minimal of the trimer

potential

6. Metastable and Scattering States of the Dimer Tunneling

The metastable states of the dimer Be2 tunneling through Gaussian barriers of Task 3 were calculated for BVP calculated with N =15 equations in the system (19) with matrix elements of potential barrier on the finite element grid = {pmin =1.81(12)4.21(15)pmax = 7.51} with the fifth-order Hermitian elements (p = 5). The corresponding problem of a dimer tunneling through the barriers was solved on the same grid.

The corresponding algebraic eigenvalue problem for metastable states was solved using the above mentioned continuous analogue of Newton's method. As the initial approximation the real eigenvalues and the corresponding orthonormalized eigenfunctions (31) were used. They were found as a result of solving the bound state problem with = 0 on the grid Qp = {pmin =1.81(12)pmax = 4.21}. The complex values of energy of the metastable states E= Em = ReEm + zlrnEm for the dimer Be2 tunneling through the Gaussian barriers, are presented in Fig. 6.

These metastable states are responsible for the resonance values of energy, corresponding to the maximal transmission coefficient, i.e., the quantum transparency of the potential barriers (see Figs. 1 (b) and 3 (b)), i.e., the shape resonances. The position of peaks presented in Fig. 6. is seen to be in quantitative agreement with the real part Re(E), and the geometric halfwidth of the |T111 (E) peaks agrees by the order of magnitude with the imaginary part r = -2Im(E) of the complex energy eigenvalues E = Re(E) + zIm(E) of the metastable states. The obtained complex energy values corresponding to the resonance values of energy in the first open channel are in good agreement with the ones calculated analytically in the model of a rigid diatomic molecule with Morse potential tunneling trough the Gaussian barrier at the same values of parameters [13]. From Fig. 6 one can see that as the energy of the initial excited state increases, the transmission peaks demonstrate a shift towards higher energies, the set of peak positions keeping approximately the same as for the transitions from the ground state and the peaks just replace each other. For example, the left epure shows that the positions of the 13th and 14th peaks for transitions from the first state coincide with the positions of the 1st and

2nd peaks for the transitions from the second state, while the right epure shows that the positions of the 25th and 26th peaks for transitions from the first state coincide with the positions of the 13th and 14th peaks for transitions from the second state and with the positions of the 1st and 2nd peaks for the transitions from the third state.

1000 -900 -800 -700 -600 -500 -400 -300 -200 -100 0

Re E,K

Figure 6. The total probability |T|2(E) = Tji (lines) of penetration of the

dimer for the initial states i through the repulsive Gaussian potential barriers versus the total energy E = Re E counted from the main threshold E = 0. The values of the threshold energies E = £j, i = 1,..., 5 corresponding to the energies of ground and excited initial states are shown by arrows. The real Re E and imaginary (— 1) Im E part (with negative sign) of the metastable states energy

(circles)

7. Conclusion

The model for three atomic beryllium system in a straight line was formulated as a 2D boundary-value problem for the Schrodinger equation in Jacobi and polar coordinates. Using the Kantorovich expansions this problem has been reduced to the boundary-value problem for a set of second-order ordinary differential equations.

The efficiency of the elaborated method, algorithms and programs has been demonstrated by benchmark calculations of the resonance scattering, metastable and bound states of the considered models and also by a comparison of results for bound states of the three atomic system in the framework of direct solving BVP by FEM and Kantorovich reduction.

The effects of resonant quantum transmission of diatomic molecule through the potential barrier and reflection from the effective potential well of a three atomic system (see Figs. 3 (b) and 3 (a)), arising in the scattering process were revealed, that are generated by metastable states of the composite system (cluster + barrier or well) with complex energy eigenvalues below the dissociation threshold of dimer, corresponding to the shape and Feshbach resonances, respectively.

The elaborated method, algorithms and programs [8-10] for solving the three-atomic scattering problem as well as diatomic molecule tunneling through potential barrier can be applied to the further analysis of quantum transparency and reflection effects [13], quantum diffusion [22] and the resonance scattering in triatomic systems using modern theoretical and experimental results [23,24].

References

1. P. Ciarlet, The Finite Element Method for Elliptic Problems, North-holland Publ. Comp, Amsterdam, 1978. doi:10.1137/1.9780898719208.

2. V. G. Korneev, Schemes of the Finite Element Method for High Orders of Accuracy, Leningrad State University, Leningrad, 1977, in Russian.

3. A. A. Gusev, V. P. Gerdt, O. Chuluunbaatar, G. Chuluunbaatar, S. I. Vinitsky, V. L. Derbov, A. Gozdz, Symbolic-Numerical Algorithm for Generating Interpolation Multivariate Hermite Polynomials of High-Accuracy Finite Element Method, Lecture Notes in Computer Science 10490 (2017) 134-150. doi:10.1007/978-3-319-66320-3_11.

4. A. A. Gusev, S. I. Vinitsky, O. Chuluunbaatar, G. Chuluunbaatar, V. P. Gerdt, V. L. Derbov, A. Góódó, P. M. Krassovitskiy, High-Accuracy Finite Element Method: Benchmark Calculations, European Physics Journal - Web of Conferences 173 (2018) 03009.

5. A. A. Gusev, S. I. Vinitsky, O. Chuluunbaatar, G. Chuluunbaatar, V. P. Gerdt, V. L. Derbov, A. Gozdz, P. M. Krassovitskiy, Interpolation Hermite Polynomials For Finite Element Method, European Physics Journal - Web of Conferences 173 (2018) 03010.

6. A. A. Gusev, S. I. Vinitsky, O. Chuluunbaatar, G. Chuluunbaatar, V. P. Gerdt, V. L. Derbov, A. Gozdz, P. M. Krassovitskiy, Algorithm for calculating interpolation hermite polynomials for high-accuracy finite element method, in: Computer Algebra: International Conference Materials, Plekhanov Russian University of Economics, 2017, pp. 89-95.

7. A. A. Gusev, V. P. Gerdt, O. Chuluunbaatar, G. Chuluunbaatar, S. I. Vinitsky, V. L. Derbov, A. Gozdz, Symbolic-Numerical Algorithms for Solving the Parametric Self-Adjoint 2D Elliptic Boundary-Value Problem Using High-Accuracy Finite Element Method, Lecture Notes in Computer Science 10490 (2017) 151-166. doi:10.1007/978-3-319-66320-3-12.

8. A. A. Gusev, L. L. Hai, O. Chuluunbaatar, S. I. Vinitsky, Program KANTBP 4M for Solving Boundary-Value Problems for Systems of Ordinary Differential Equations of the Second Order.

URL http://wwwinfo.jinr.ru/programs/jinrlib/kantbp4m

9. A. A. Gusev, O. Chuluunbaatar, S. I. Vinitsky, A. G. Abrashkevich, KANTBP 3.0: New version of a Program for Computing Energy Levels, Reflection and Transmission Matrices, and Corresponding Wave Functions in the Coupled-Channel Adiabatic Approach, Computer Physics Communications 185 (2014) 3341-3343. doi:10.1016/j.cpc.2014.08.002.

10. O. Chuluunbaatar, A. A. Gusev, S. I. Vinitsky, A. G. Abrashkevich, ODPEVP: A Program for Computing Eigenvalues and Eigenfunctions and Their First Derivatives with Respect to the Parameter of the Parametric Self-Adjoined Sturm-Liouville Problem, Computer Physics Communications 180 (2009) 1358-1375. doi:10.1016/j.cpc.2009.04.017.

11. A. A. Gusev, O. Chuluunbaatar, S. I. Vinitsky, V. L. Derbov, Algorithms for Solving the Boundary-Value Problems for Atomic Trimers in Collinear Configuration using the Kantorovich Method, Bulletin of Peoples' Friendship University of Russia. Series: Mathematics. Information Sciences. Physics (4) (2016) 56-76.

12. J. F. Cornwell, Group Theory in Physics, Academic Press, New York, 1984.

13. P. M. Krassovitskiy, F. M. Pen'kov, Contribution of Resonance Tunneling of Molecule to Physical Observables, Journal of Physics B: Atomic, Molecular and Optical Physics 47 (22) (2014) 225210. doi:10.1088/0953-4075/47/22/225210.

14. J. Wang, G. Wang, J. Zhao, Density Functional Study of Beryllium Clusters, with Gradient Correction, Journal of Physics: Condensed Matter 13 (33) (2001) L753-L758. doi:10.1088/0953-8984/13/33/101.

15. L. J. Lauhon, W. Ho, Direct Observation of the Quantum Tunneling of Single Hydrogen Atoms with a Scanning Tunneling Microscope, Physical Review Letters 85 (2000) 4566-4569. doi:10.1103/PhysRevLett.85.4566.

16. L. V. Kantorovich, V. I. Krylov, Approximate Methods of Higher Analysis, Wiley, New York, 1964.

17. M. Abramovits, I. A. Stegun, Handbook of Mathematical Functions, Dover, New York, 1972.

18. R. G. Newton, Analytic Properties of Radial Wave Functions, Journal of Mathematical Physics 1 (1960) 319-348. doi:10.1063/1.1703680.

19. A. A. Gusev, S. I. Vinitsky, O. Chuluunbaatar, V. L. Derbov, A. Gozdz, P. M. Krassovitskiy, Metastable States of a Composite System Tunneling through Repulsive Barriers, Theoretical and Mathematical Physics 186 (2016) 21-40. doi:10.1134/S0040577916010037.

20. A. A. Gusev, O. Chuluunbaatar, S. I. Vinitsky, V. L. Derbov, A. Gozdz, L. L. Hai, V. A. Rostovtsev, Symbolic-Numerical Solution of Boundary-Value Problems with Self-Adjoint Second-Order Differential Equation using the Finite Element Method with Interpolation Hermite Polynomials, Lecture Notes in Computer Science 8660 (2014) 138-154. doi:10.1007/978-3-319-10515-4_11.

21. I. V. Puzynin, T. L. Boyadjiev, S. I. Vinitsky, E. V. Zemlyanaya, T. P. Puzynina, O. Chuluunbaatar, Methods of Computational Physics for Investigation of Models of Complex Physical Systems, Physics of Particles and Nuclei 38 (2007) 70-116. doi:10.1134/S1063779607010030.

22. E. Pijper, A. Fasolino, Quantum Surface Diffusion of Vibrationally Excited Molecular Dimers, Journal of Chemical Physics 126 (2007) 014708. doi:10.1063/1.2424699.

23. A. V. Mitin, Unusual Chemical Bonding in the Beryllium Dimer and its Twelve Vibrational Levels, Chemical Physics Letters 682 (2017) 30-33. doi:10.1016/j.cplett.2017.05.071.

24. J. M. Merritt, V. E. Bondybey, M. C. Heaven, Beryllium Dimer-Caught in the Act of Bonding, Science 324 (2009) 1548-1551. doi:10.1126/science.1174326.

УДК 519.632.4

Б01: 10.22363/2312-9735-2018-26-3-226-243

Метод конечных элементов высокого порядка точности для решения двухмерных эллиптических краевых задач двух и трёх тождественных атомов на прямой А. А. Гусев

Лаборатория информационных технологий Объединённый институт ядерных исследований ул. Жолио-Кюри, д. 6, г. Дубна, Московской обл., 141980, Россия

Рассмотрены модели трёх одинаковых атомов на прямой с парным молекулярным взаимодействием и рассеяние двухатомной молекулы на атоме или её туннелирования через потенциальные барьеры. Модели сформулированы в виде двумерных эллиптических краевых задач (КЗ) в координатах Якоби и полярных координатах. КЗ в координатах Якоби решаются методом конечных элементов высокого порядка точности для дискретного спектра рассматриваемых моделей. Для решения задач рассеяния КЗ в полярных координатах с помощью метода Канторовича сводится к системе обыкновенных дифференциальных уравнений второго порядка по радиальной переменной с использованием разложения искомых решений по набору угловых базисных функций, параметрически зависящих от радиальной переменной. Эффективность разработанного метода, алгоритмов и программ демонстрируется путём эталонных расчётов резонансного рассеяния, метастабильных и связанных состояний рассматриваемых моделей, а также путём сравнения результатов для связанных состояний трёх атомных систем в рамках прямого решения КЗ методом конечных элементов и редукции Канторовича.

Ключевые слова: эллиптические краевые задачи, задача рассеяния, метастабильные и связанные состояния, метода Канторовича

Литература

1. Ciarlet P. The Finite Element Method for Elliptic Problems. — Amsterdam: North-holland Publ. Comp, 1978. — DOI: 10.1137/1.9780898719208.

2. Корнеев В. Г. Схемы метода конечных элементов высоких порядков точности. — Ленинград: Издательство Ленинградского университета, 1977.

3. Symbolic-Numerical Algorithm for Generating Interpolation Multivariate Hermite Polynomials of High-Accuracy Finite Element Method / A. A. Gusev, V. P. Gerdt, O. Chuluunbaatar et al. // Lecture Notes in Computer Science. — 2017. — Vol. 10490. — Pp. 134-150. — DOI: 10.1007/978-3-319-66320-3_11.

4. High-Accuracy Finite Element Method: Benchmark Calculations / A. A. Gusev, S. I. Vinitsky, O. Chuluunbaatar et al. // European Physics Journal - Web of Conferences. — 2018. — Vol. 173. — P. 03009.

5. Interpolation Hermite Polynomials For Finite Element Method / A. A. Gusev, S. I. Vinitsky, O. Chuluunbaatar et al. // European Physics Journal - Web of Conferences. — 2018. — Vol. 173. — P. 03010.

6. Algorithm for Calculating Interpolation Hermite Polynomials for High-Accuracy Finite Element Method / A. A. Gusev, S. I. Vinitsky, O. Chuluunbaatar et al. // Computer Algebra: International Conference Materials / Plekhanov Russian University of Economics. — 2017. — Pp. 89-95.

7. Symbolic-Numerical Algorithms for Solving the Parametric Self-Adjoint 2D Elliptic Boundary-Value Problem Using High-Accuracy Finite Element Method / A. A. Gusev, V. P. Gerdt, O. Chuluunbaatar et al. // Lecture Notes in Computer Science. — 2017. — Vol. 10490. — Pp. 151-166. — DOI: 10.1007/978-3-319-66320-3_12.

8. Gusev A. A, Hai L. L., Chuluunbaatar O., Vinitsky S. I. Program KANTBP 4M for Solving Boundary-Value Problems for Systems of Ordinary Differential Equations of the Second Order. — http://wwwinfo.jinr.ru/programs/jinrlib/kantbp4m.

9. KANTBP 3.0: New version of a Program for Computing Energy Levels, Reflection and Transmission Matrices, and Corresponding Wave Functions in the CoupledChannel Adiabatic Approach / A. A. Gusev, O. Chuluunbaatar, S. I. Vinitsky, A. G. Abrashkevich // Computer Physics Communications. — 2014. — Vol. 185, issue 12. — Pp. 3341-3343. — DOI: 10.1016/j.cpc.2014.08.002.

10. ODPEVP: A Program for Computing Eigenvalues and Eigenfunctions and Their First Derivatives with Respect to the Parameter of the Parametric Self-Adjoined Sturm-Liouville Problem / O. Chuluunbaatar, A. A. Gusev, S. I. Vinitsky, A. G. Abrashkevich // Computer Physics Communications. — 2009. — Vol. 180, issue 8. — Pp. 1358-1375. — DOI: 10.1016/j.cpc.2009.04.017.

11. Алгоритмы решения краевых задач для атомных тримеров в коллинеарной конфигурации методом Канторовича / А. А. Гусев, О. Чулуунбаатар, С. И. Ви-ницкий, В. Л. Дербов // Вестник РУДН: Серия Математика. Информатика. Физика. — 2016. — № 4. — С. 56-76.

12. Cornwell J. F. Group Theory in Physics. — New York: Academic Press, 1984.

13. Krassovitskiy P. M., Pen'kov F. M. Contribution of Resonance Tunneling of Molecule to Physical Observables // Journal of Physics B: Atomic, Molecular and Optical Physics. — 2014. — Vol. 47, No 22. — P. 225210. — DOI: 10.1088/09534075/47/22/225210.

14. Wang J., Wang G., Zhao J. Density Functional Study of Beryllium Clusters, with Gradient Correction // Journal of Physics: Condensed Matter. — 2001. — Vol. 13, No 33. — Pp. L753-L758. — DOI: 10.1088/0953-8984/13/33/101.

15. Lauhon L. J., Ho W. Direct Observation of the Quantum Tunneling of Single Hydrogen Atoms with a Scanning Tunneling Microscope // Physical Review Letters. — 2000. — Vol. 85. — Pp. 4566-4569. — DOI: 10.1103/PhysRevLett.85.4566.

16. Kantorovich L. V., Krylov V. I. Approximate Methods of Higher Analysis. — New York: Wiley, 1964.

17. Абрамовиц М., Стиган И. Справочник по специальным функциям. — М.: Наука, 1979.

18. Newton R. G. Analytic Properties of Radial Wave Functions // Journal of Mathematical Physics. — 1960. — Vol. 1, issue 5. — Pp. 319-348. — DOI: 10.1063/1.1703680.

19. Метастабильные состояния составной системы при туннелировании через отталкивающие барьеры / А. А. Гусев, С. И. Виницкий, О. Чулуунбаатар и др. // Теоретическая математическая физика. — 2016. — Т. 186. — С. 27-50. — DOI: 10.4213/tmf8981.

20. Symbolic-Numerical Solution of Boundary-Value Problems with Self-Adjoint Second-Order Differential Equation using the Finite Element Method with Interpolation Hermite Polynomials / A. A. Gusev, O. Chuluunbaatar, S. I. Vinitsky et al. // Lecture Notes in Computer Science. — 2014. — Vol. 8660. — Pp. 138-154. — DOI: 10.1007/978-3-319-10515-4_11.

21. О методах вычислительной физики для исследования моделей сложных физических процессов / И. В. Пузынин, Т. Л. Бояджиев, С. И. Виницкий и др. // Физика элементарных частиц и атомного ядра. — 2007. — Т. 38. — С. 144-232.

22. Pijper E., Fasolino A. Quantum Surface Diffusion of Vibrationally Excited Molecular Dimers // Journal of Chemical Physics. — 2007. — Vol. 126, issue 1. — P. 014708. — DOI: 10.1063/1.2424699.

23. Mitin A. V. Unusual Chemical Bonding in the Beryllium Dimer and its Twelve Vibrational Levels // Chemical Physics Letters. — 2017. — Vol. 682. — Pp. 30-33. — DOI: 10.1016/j.cplett.2017.05.071.

24. Merritt J. M., Bondybey V. E., Heaven M. C. Beryllium Dimer-Caught in the Act of Bonding // Science. — 2009. — Vol. 324. — Pp. 1548-1551. — DOI: 10.1126/sci-ence.1174326.

© Gusev A. A., 2018

Для цитирования:

Gusev A. A. Finite Element Method of High-Order Accuracy for Solving Two Dimensional Elliptic Boundary-Value Problems of Two and Three Identical Atoms in a Line // RUDN Journal of Mathematics, Information Sciences and Physics. — 2018. — Vol. 26, No 3. — Pp. 226-243. — DOI: 10.22363/2312-9735-2018-26-3-226-243.

For citation:

Gusev A. A. Finite Element Method of High-Order Accuracy for Solving Two Dimensional Elliptic Boundary-Value Problems of Two and Three Identical Atoms in a Line, RUDN Journal of Mathematics, Information Sciences and Physics 26 (3) (2018) 226-243. DOI: 10.22363/2312-9735-2018-26-3-226-243.

Сведения об авторах:

Гусев Александр Александрович — кандидат физико-математических наук, старший научный сотрудник Лаборатории информационных технологий Объединённого института ядерных исследований (e-mail: [email protected], тел.: +7(49621)63536)

Information about the authors:

Gusev A. A. — Candidate of Physical and Mathematical Sciences, senior researcher of Laboratory of Information Technologies Joint Institute for Nuclear Research (e-mail: [email protected], phone: +7(49621)63536)