Skip to main content

Solutions of the Schrödinger equation in a Hilbert space

Abstract

Necessary and sufficient conditions for the existence of a solution of a boundary-value problem for the Schrödinger equation are obtained in the linear and nonlinear cases. Analytic solutions are represented using the generalized Green operator.

Introduction

The Schrödinger equation is the subject of numerous publications, and it is impossible to analyze all of them in detail. For this reason, we only briefly describe the methods and ideas that underlie the approach proposed in this paper for the investigation of the linear and a weakly nonlinear Schrödinger equation with different boundary conditions.

In this work, we develop constructive methods of analysis of linear and weakly nonlinear boundary-value problems, which occupy a central place in the qualitative theory of differential equations. The specific feature of these problems is that the operator of the linear part of the equation does not have an inverse. This does not allow one to use the traditional methods based on the principles of contracting mappings and a fixed point. These problems include the most complicated and inadequately studied problems known as critical (or resonance) problems [14]. Therefore, for the investigation of periodic problems for the Schrödinger equation, we develop the technique of generalized inverse operators [58] for the original linear operator in Banach and Hilbert spaces.

On the other hand, we use the notion of a strong generalized solution of an operator equation developed in [9]. The origins of this approach go back to the works of Weil and Sobolev. Using the process of completion, one can introduce the concept of a strong pseudoinverse operator for an arbitrary linear bounded operator and thus relax the requirement that the range of its values be closed. In this way, one can prove the existence of solutions of different types for the linear Schrödinger equation with arbitrary inhomogeneities. Thus, one may say that, in a certain sense, the Schrödinger equation is always solvable. There are three possible types of solutions: classical generalized solutions, strong generalized solution, and strong pseudosolutions [10].

For the analysis of a weakly nonlinear Schrödinger equation, we develop the ideas of the Lyapunov-Schmidt method and efficient methods of perturbation theory, namely the Vishik-Lyusternik method [11]. The combination of different approaches allows us to take a different look at the Schrödinger equation with a constant unbounded operator in the linear part and obtain all its solutions by using the generalized Green operator of this problem constructed in this work. Possible generalizations are discussed in the final part of the paper. By an example of the abstract van der Pol equation, we illustrate the results that can be obtained by using the proposed method.

Auxiliary result (linear case)

Statement of the problem

Consider the following boundary-value problem for the Schrödinger equation in a Hilbert space H T :

d φ ( t ) d t =i H 0 φ(t)+f(t),t[0;w],
(1)
φ(0)φ(w)=αD,
(2)

where H T =HH, H is Hilbert space and vector-function f(t) is integrable; for simplicity, the unbounded operator H 0 has the following form [12] for each t[0;w]:

H 0 =i ( 0 T T 0 ) = ( T 0 0 T ) ( 0 I I 0 ) =i ( 0 I I 0 ) ( T 0 0 T ) .

In a more general case, the operator H 0 has the form

H 0 =iJ ( T 0 0 T ) =i ( T 0 0 T ) J,J= J = J 1 ,

where T is a strongly positive self-adjoint operator in the Hilbert space H. Since the operator T is closed, the domain D(T) of the operator T is a Hilbert space with scalar product (Tu,Tu). The operator H 0 is self-adjoint in the domain D=D(T)D(T) with product

( u , v , u , v ) H T = ( T u , T u ) H + ( T v , T v ) H ,

the infinitesimal generator of a strongly continuous evolution semigroup has the form

U ( t ) : = U ( t , 0 ) = ( cos t T sin t T sin t T cos t T ) , U n ( t ) = ( cos n t T sin n t T sin n t T cos n t T ) = U ( n t ) ,

U n (t)=1, nN (nonexpanding group), φ(t)= ( φ 1 ( t ) , φ 2 ( t ) ) T , α= ( α 1 , α 2 ) T , and f(t)= ( f 1 ( t ) , f 2 ( t ) ) T . The mild solutions of equation (1) can be represented in the form

φ(t)=U(t)c+ 0 t U(t) U 1 (τ)f(τ)dτ,

for any element c H T . Substituting this in condition (2), we conclude that the solvability of the boundary-value problem (1), (2) is equivalent to the solvability of the following operator equation:

( I U ( w ) ) c=g,
(3)

where g=α+U(w) 0 w U 1 (τ)f(τ)dτ. Consider the case where the set of values of IU(w) is closed R(IU(w))= R ( I U ( w ) ) ¯ . Since U n (w)=U(wn)=1 for all nN, we can conclude [13] that the operator system (3) is solvable if and only if

U 0 (w)g=0,
(4)

where

U 0 (w)= lim n k = 0 n U k ( w ) n = lim n k = 0 n U ( k w ) n ,

is the orthoprojector that projects the space H T onto the subspace 1σ(U(w)). Under this condition, the solutions of (3) have the form

c= U 0 (w) c ¯ + ( k = 0 ( μ 1 ) k { l = 0 μ l 1 ( U ( w ) U 0 ( w ) ) l } k + 1 U 0 ( w ) ) g,

for 0<μ1< 1 R μ ( U ( w ) ) , and any c ¯ H T . Then we can formulate the first result as a lemma.

Lemma 1 Suppose that the operator IU(w) has a closed image R(IU(w))= R ( I U ( w ) ) ¯ .

  1. 1.

    Solutions of the boundary-value problem (1), (2) exist if and only if

    U 0 (w) ( α + 0 w U 1 ( τ ) f ( τ ) d τ ) =0.
    (5)
  2. 2.

    Under condition (5), solutions of (1), (2) have the form

    φ(t, c ¯ )=U(t) U 0 (w) c ¯ + ( G [ f , α ] ) (t),
    (6)

    where

    ( G [ f , α ] ) ( t ) = U ( t ) k = 0 ( μ 1 ) k { l = 0 μ l 1 ( U ( w ) U 0 ( w ) ) l } k + 1 × ( α + 0 w U ( w ) U 1 ( τ ) f ( τ ) d τ ) U ( t ) U 0 ( w ) ( α + 0 w U ( w ) U 1 ( τ ) f ( τ ) d τ ) + 0 t U ( t ) U 1 ( τ ) f ( τ ) d τ ,

    is the generalized Green operator of the boundary-value problem (1), (2) for 0<μ1<1/ R μ (U(w)).

We now show that the condition R(IU(w))= R ( I U ( w ) ) ¯ of Lemma 1 can be omitted and, in different senses, the boundary-value problem (1), (2) is always solvable.

  1. (1)

    Classical generalized solutions.

Consider the case where the set of values of IU(w) is closed (R(IU(w))= R ( I U ( w ) ) ¯ ). Then [5]gR(IU(w)) if and only if P N ( ( I U ( w ) ) ) g=0, and the set of solutions of (3) has the form [5]c=G[g]+ U 0 (w) c ¯ , c ¯ H T , where [5, 13]

G[g]= ( I U ( w ) ) + g= ( I ( U ( w ) U 0 ( w ) ) 1 U 0 ( w ) ) g

is the generalized Green operator (or it has the form of a convergent series).

  1. (2)

    Strong generalized solutions. Consider the case where R(IU(w)) R ( I U ( w ) ) ¯ and g R ( I U ( w ) ) ¯ . We show that the operator IU(w) can be extended to I U ( w ) ¯ in such a way that R( I U ( w ) ¯ ) is closed.

Since the operator IU(w) is bounded, the following representation of H T in the form of a direct sum is true:

H T =N ( I U ( w ) ) X, H T = R ( I U ( w ) ) ¯ Y,

where X=N ( I U ( w ) ) = R ( I U ( w ) ) ¯ and Y= R ( I U ( w ) ) ¯ =N(IU(w)). Let E= H T /N(IU(w)) be the quotient space of H T and let P R ( I U ( w ) ) ¯ and P N ( I U ( w ) ) be the orthoprojectors onto R ( I U ( w ) ) ¯ and N(IU(w)), respectively. Then the operator

I U (w)= P R ( I U ( w ) ) ¯ ( I U ( w ) ) j 1 p:XR ( I U ( w ) ) R ( I U ( w ) ) ¯

is linear, continuous, and injective. Here,

p:XE= H T /N ( I U ( w ) ) andj: H T E,

are a continuous bijection and a projection, respectively. The triple ( H T ,E,j) is a locally trivial bundle with typical fiber H 1 = P N ( I U ( w ) ) H T [14]. In this case [[9], p.26,29], we can define a strong generalized solution of the equation

( I U ( w ) ) x=g,xX.
(7)

We complete the space X with the norm x X ¯ = ( I U ( w ) ) x F , where F= R ( I U ( w ) ) ¯ [9]. Then the extended operator

I U ( w ) ¯ : X ¯ R ( I U ( w ) ) ¯ ,X X ¯

is a homeomorphism of X ¯ and R ( I U ( w ) ) ¯ . By the construction of a strong generalized solution [9], the equation

( I U ( w ) ¯ ) ξ ¯ =g,

has a unique solution ( I U ( w ) ¯ ) 1 g, which is called the generalized solution of equation (7).

Remark 1 It should be noted that there exist the following extensions of spaces and the corresponding operators:

p ¯ : X ¯ E ¯ , j ¯ : H ¯ T E ¯ , P ¯ X = P X ¯ : H ¯ T X ¯ , G ¯ : R ( I U ( w ) ) ¯ X ¯ ,

where

H ¯ T = N ( I U ( w ) ) X ¯ ; p ¯ ( x ) = p ( x ) , x X ; j ¯ ( x ) = j ( x ) , x H T , P ¯ X ( x ) = P X ( x ) , x H T ( P X = P X 2 = P X ) ; G ¯ [ g ] = G [ g ] , g R ( I U ( w ) ) .

Then the operator I U ( w ) ¯ =( I U ( w ) ¯ ) P X ¯ : H ¯ T H T is an extension of IU(w), and ( I U ( w ) ) ¯ c=(IU(w))c for all c H T .

  1. (3)

    Strong pseudosolutions.

Consider an element g R ( I U ( w ) ) ¯ . This condition is equivalent to P N ( ( I U ( w ) ) ) g0. In this case, there are elements of H ¯ T that minimize the norm ( I U ( w ) ) ¯ ξ g H T :

ξ= ( I U ( w ) ) ¯ 1 g+ P N ( I U ( w ) ) c ¯ , c ¯ H T .

These elements are called strong pseudosolutions by analogy with [5].

We now formulate the full theorem on solvability.

Theorem 1 The boundary-value problem (1), (2) is always solvable.

  1. (1)

    (a) Classical or strong generalized solutions of (1), (2) exist if and only if

    U 0 (w) ( α + 0 w U 1 ( τ ) f ( τ ) d τ ) =0.
    (8)

    If (α+ 0 w U 1 (τ)f(τ)dτ)R(IU(w)), then solutions of (1), (2) are classical.

  2. (b)

    Under assumption (8), the solutions of (1), (2) have the form

    φ(t, c ¯ )=U(t) U 0 (w) c ¯ + ( G [ f , α ] ¯ ) (t),

    where ( G [ f , α ] ¯ )(t) is an extension of the operator (G[f,α])(t).

  3. (2)

    (a) Strong pseudosolutions exist if and only if

    U 0 (w) ( α + 0 w U 1 ( τ ) f ( τ ) d τ ) 0.
    (9)
  4. (b)

    Under assumption (9), the strong pseudosolutions of (1), (2) have the form

    φ(t, c ¯ )=U(t) U 0 (w) c ¯ + ( G [ f , α ] ¯ ) (t),

    where

    ( G [ f , α ] ¯ ) ( t ) = U ( t ) G ¯ [ g ] + 0 t U ( t ) U 1 ( τ ) f ( τ ) d τ = U ( t ) ( I U ( w ) ) ¯ 1 g + 0 t U ( t ) U 1 ( τ ) f ( τ ) d τ .

1 Main result (nonlinear case)

1.1 Modification of the Lyapunov-Schmidt method

In the Hilbert space H T defined above, we consider the boundary-value problem

d φ ( t ) d t =i H 0 φ(t)+εZ ( φ ( t ) , t , ε ) +f(t),
(10)
φ(0,ε)φ(w,ε)=α.
(11)

We seek a generalized solution φ(t,ε) of the boundary-value problem (10), (11) that becomes one of the solutions of the generating equation (1), (2) φ 0 (t, c ¯ ) in the form (6) for ε=0.

To find a necessary condition for the operator function Z(φ,t,ε), we impose the joint constraints

Z(,,)C ( [ 0 ; w ] , H T ) ×C[0, ε 0 ]×C [ φ φ 0 q ] ,

where q is a positive constant.

The main idea of the next results was used in [15] for the investigation of bounded solutions.

Let us show that this problem can be solved with the use of the following operator equation for generating amplitudes:

F( c ¯ )= U 0 (w) 0 w U 1 (τ)Z ( φ 0 ( τ , c ¯ ) , τ , 0 ) dτ=0.
(12)

Theorem 2 (Necessary condition)

Suppose that the nonlinear boundary-value problem (10), (11) has a generalized solution φ(,ε) that becomes one of the solutions φ 0 (t, c ¯ ) of the generating equation (1), (2) with constant c ¯ = c 0 and φ(t,0)= φ 0 (t, c 0 ) for ε=0. Then this constant must satisfy the equation for generating amplitudes (12).

Proof If the boundary-value problem (10), (11) has classical generalized solutions, then, by Lemma 1, the following solvability condition must be satisfied:

U 0 (w) ( α + 0 w U 1 ( τ ) { f ( τ ) + ε Z ( φ ( τ , ε ) , τ , ε ) } d τ ) =0.
(13)

By using condition (5), we establish that condition (13) is equivalent to the following:

U 0 (w) 0 w U 1 (τ)Z ( φ ( τ , ε ) , τ , ε ) dτ=0.

Since φ(t,ε) φ 0 (t, c 0 ) as ε0, we finally obtain [by using the continuity of the operator function Z(φ,t,ε)] the required assertion.

To find a sufficient condition for the existence of solutions of the boundary-value problem (10), (11), we additionally assume that the operator function Z(φ,t,ε) is strongly differentiable in a neighborhood of the generating solution (Z(,t,ε) C 1 [φ φ 0 q]).

This problem can be solved with the use of the operator

B 0 = d F ( c ¯ ) d c ¯ | c ¯ = c 0 = U 0 (w) 0 w U 1 (t) A 1 (t)dt:HH,

where A 1 (t)= Z 1 (v,t,ε) | v = φ 0 , ε = 0 (Fréchet derivative). □

Theorem 3 (Sufficient condition)

Suppose that the operator B 0 satisfies the following conditions:

  1. (1)

    The operator B 0 is Moore-Penrose pseudoinvertible;

  2. (2)

    P N ( B 0 ) U 0 (w)=0.

Then, for an arbitrary element c= c 0 H T satisfying the equation for generating amplitudes (12), there exists at least one solution of (10), (11).

This solution can be found by using the following iterative process:

v ¯ k + 1 ( t , ε ) = ε G [ Z ( φ 0 ( τ , c 0 ) + v k , τ , ε ) , α ] ( t ) , c k = B 0 + U 0 ( w ) 0 w U 1 ( τ ) { A 1 ( τ ) v ¯ k ( τ , ε ) + R ( v k ( τ , ε ) , τ , ε ) } d τ , R ( v k ( t , ε ) , t , ε ) = Z ( φ 0 ( t , c 0 ) + v k ( t , ε ) , t , ε ) Z ( φ 0 ( t , c 0 ) , t , 0 ) A 1 ( t ) v k ( t , ε ) , R ( 0 , t , 0 ) = 0 , R x 1 ( 0 , t , 0 ) = 0 , v k + 1 ( t , ε ) = U ( t ) U 0 ( w ) c k + v ¯ k + 1 ( t , ε ) , φ k ( t , ε ) = φ 0 ( t , c 0 ) + v k ( t , ε ) , k = 0 , 1 , 2 , , v 0 ( t , ε ) = 0 , φ ( t , ε ) = lim k φ k ( t , ε ) .

1.2 Relationship between necessary and sufficient conditions

First, we formulate the following assertion:

Corollary Suppose that a functional F( c ¯ ) has the Fréchet derivative F ( 1 ) ( c ¯ ) for each element c 0 of the Hilbert space H satisfying the equation for generating constants (12). If F 1 ( c ¯ ) has a bounded inverse, then the boundary-value problem (10), (11) has a unique solution for each c 0 .

Remark 2 If the assumptions of the corollary are satisfied, then it follows from its proof that the operators B 0 and F ( 1 ) ( c 0 ) are equal. Since the operator F ( 1 ) ( c ¯ ) is invertible, it follows that assumptions 1 and 2 of Theorem 3 are necessarily satisfied for the operator B 0 . In this case, the boundary-value problem (10), (11) has a unique bounded solution for each c 0 H T satisfying (12). Therefore, the invertibility condition for the operator F 1 ( c ¯ ) expresses the relationship between the necessary and sufficient conditions. In the finite-dimensional case, the condition of invertibility of the operator F ( 1 ) ( c ¯ ) is equivalent to the condition of simplicity of the root c 0 of the equation for generating amplitudes [5].

In this way, we modify the well-known Lyapunov-Schmidt method. It should be emphasized that Theorems 2 and 3 give us a condition for the chaotic behavior of (10) and (11) [16].

1.3 Example

We now illustrate the obtained assertion. Consider the following differential equation in a separable Hilbert space H:

y ¨ (t)+Ty(t)=ε ( 1 y ( t ) 2 ) y ˙ (t),
(14)
y(0)=y(w), y ˙ (0)= y ˙ (w),
(15)

where T is an unbounded operator with compact T 1 . Then there exists an orthonormal basis e i H such that y(t)= i = 1 c i (t) e i and Ty(t)= i = 1 λ i c i (t) e i , λ i . In this case, the operator system (10), (11) for the boundary-value problem (14), (15) is equivalent to the following countable system of ordinary differential equations ( c k (t)= x k (t)):

x ˙ k ( t ) = λ k y k ( t ) , k = 1 , 2 , , y ˙ k ( t ) = λ k x k ( t ) + ε λ k ( 1 j = 1 x j 2 ( t ) ) y k ( t ) ,
(16)
x k (0)= x k (w), y k (0)= y k (w).
(17)

We find the solutions of these equations in the space W 2 1 ([0;w]) that, for ε=0, turn into one of the solutions of the generating equation. Consider the critical case λ i =4 π 2 i 2 / w 2 , iN. Let w=2π. In this case, the set of all periodic solutions of (16), (17) has the form

x k ( t ) = cos ( k t ) c 1 k + sin ( k t ) c 2 k , y k ( t ) = sin ( k t ) c 1 k + cos ( k t ) c 2 k

for all pairs of constants c 1 k , c 2 k R, kN. The equation for generating amplitudes (12) is equivalent in this case to the following countable systems of algebraic nonlinear equations:

( c 1 k ) 3 + 2 j = 1 , j k ( c 1 k ( c 1 j ) 2 + c 1 k ( c 2 j ) 2 ) + c 1 k ( c 2 k ) 2 4 c 1 k = 0 , ( c 2 k ) 3 + 2 j = 1 , j k ( c 2 k ( c 1 j ) 2 + c 2 k ( c 2 j ) 2 ) + ( c 1 k ) 2 c 2 k 4 c 2 k = 0 , k N .

Then we can obtain the next result.

Theorem 4 (Necessary condition for the van der Pol equation)

Suppose that the boundary-value problem (16), (17) has a bounded solution φ(,ε) that becomes one of the solutions of the generating equations with pairs of constants ( c 1 k , c 2 k ), kN. Then only a finite number of these pairs are not equal to zero. Moreover, if ( c 1 k i , c 2 k i )(0,0), i= 1 , N ¯ , then these constants lie on an N-dimensional torus in the infinite-dimensional space of constants:

( c 1 k i ) 2 + ( c 2 k i ) 2 = ( 2 2 N 1 ) 2 ,i= 1 , N ¯ .

Remark Similarly, we can study the Schrödinger equation with a variable operator and more general boundary conditions (as noted in the introduction).

Consider the differential Schrödinger equation

d φ ( t ) d t =iH(t)φ(t)+f(t),tJ
(18)

in a Hilbert space H with the boundary condition

Qφ()=α,
(19)

where, for each tJR, the unbounded operator H(t) has the form H(t)= H 0 +V(t), H 0 = H 0 is an unbounded self-adjoint operator with domain D=D( H 0 )H, and the mapping tV(t) is strongly continuous. The operator Q is linear and bounded and acts from the Hilbert space H to H 1 . As in [12], we define the operator-valued function

V ˜ (t)= e i t H 0 V(t) e i t H 0 .

In this case, V ˜ (t) admits the Dyson representation [[12], p.311]; denote its propagator by U ˜ (t,s). If U(t,s)= e i t H 0 U ˜ (t,s) e i s H 0 , then ψ s (t)=U(t,s)ψ is a weak solution of (14) with the condition φ s (s)=ψ in the sense that, for any ηD( H 0 ), the function (η, ψ s (t)) is differentiable and

d d t ( η , ψ s ( t ) ) =i ( H 0 η , ψ s ( t ) ) i ( V ( t ) η , ψ s ( t ) ) + ( f ( t ) , ψ s ( t ) ) ,tJ.

A detailed study of the boundary-value problem (18), (19) will be given in a separate paper.

References

  1. Vainberg MM, Trenogin SF: The Branching Theory of Solutions of Nonlinear Equations. Nauka, Moscow; 1969.

    Google Scholar 

  2. Malkin IG: Some Problems of the Theory of Nonlinear Oscillations. Nauka, Moscow; 1956. (in Russian)

    MATH  Google Scholar 

  3. Grebenikov EA, Ryabov YA: Constructive Methods of Analysis of Nonlinear Systems. Nauka, Moscow; 1979. (in Russian)

    MATH  Google Scholar 

  4. Boichuk AA, Korostil IA, Fec̆kan M: Bifurcation conditions for a solution of an abstract wave equation. Differ. Equ. 2007, 43: 481-487.

    Article  MathSciNet  Google Scholar 

  5. Boichuk AA, Samoilenko AM: Generalized Inverse Operators and Fredholm Boundary-Value Problems. VSP, Utrecht; 2004.

    Book  MATH  Google Scholar 

  6. Moore EH: On the reciprocal of the general algebraic matrix (abstract). Bull. Am. Math. Soc. 1920, 26: 394-395.

    Google Scholar 

  7. Penrose RA: A generalized inverse for matrices. Proc. Camb. Philos. Soc. 1955, 51: 406-413. 10.1017/S0305004100030401

    Article  MATH  MathSciNet  Google Scholar 

  8. Pokutnyi AA: Linear normally-resolvable equations in Banach spaces. J. Comput. Appl. Math. 2012, 107: 146-153. (in Russian)

    Google Scholar 

  9. Klyushin DA, Lyashko SI, Nomirovskii DA, Semenov VV: Generalized Solutions of Operator Equations and Extreme Elements. Springer, Berlin; 2012.

    Book  Google Scholar 

  10. Tikhonov AN, Arsenin VY: Methods for Solving Ill-Posed Problems. Nauka, Moscow; 1979. (in Russian)

    Google Scholar 

  11. Vishik MI, Lyusternik LA: The solution of some perturbation problems for matrices and selfadjoint or non-selfadjoint differential equations. Russ. Math. Surv. 1960, 15: 3-80.

    MathSciNet  Google Scholar 

  12. Reed M, Simon B 2. In Methods of Modern Mathematical Physics. Mir, Moscow; 1979.

    Google Scholar 

  13. Biletskyi BA, Boichuk AA, Pokutnyi AA: Periodic problems of difference equations and ergodic theory. Abstr. Appl. Anal. 2011., 12: Article ID 928587

    Google Scholar 

  14. Atiyah, MF: K-Theory. New York, Amsterdam (1967)

  15. Boichuk AA, Pokutnyi AA: Bounded solutions of weakly nonlinear differential equations in a Banach space. Nonlinear Oscil. 2008, 11: 158-167. 10.1007/s11072-008-0021-9

    Article  MATH  MathSciNet  Google Scholar 

  16. Chueshov ID: Introduction to the Theory of Infinite-Dimensional Dissipative Systems. Acta, Kharkiv; 2002.

    MATH  Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Alexander Boichuk.

Additional information

Competing interests

The authors did not provide this information.

Authors’ contributions

The authors did not provide this information.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Boichuk, A., Pokutnyi, O. Solutions of the Schrödinger equation in a Hilbert space. Bound Value Probl 2014, 4 (2014). https://doi.org/10.1186/1687-2770-2014-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2014-4

Keywords