arxiv: v1 [math-ph] 22 Jul 2016

Transcription

1 A NEW APPROACH FOR THE STRONG UNIQUE CONTINUATION OF ELECTROMAGNETIC SCHRÖDINGER OPERATOR WITH COMPLEX-VALUED COEFFICIENT XIAOJUN LU 1 XIAOFEN LV 2 arxiv: v1 [math-ph] 22 Jul Department of Mathematics & Jiangsu Key Laboratory of Engineering Mechanics, Southeast University, , Nanjing, China 2. Jiangsu Testing Center for Quality of Construction Engineering Co., Ltd, , Nanjing, China Abstract. This paper mainly addresses the strong unique continuation property for the electromagnetic Schrödinger operator with complex-valued coefficients. Appropriate multipliers with physical backgrounds have been introduced to prove a priori estimates. Moreover, its application in an exact controllability problem has been shown, in which case, the boundary value determines the interior value completely. Résumé. Dans cet article, on considère essentiellement la propiété fortement unique pour l opérateur électromagnétique de Schrödinger avec les coéfficients de complexe. La méthode de multiplicateur a été introduite pour démontrer les estimations à priori. En plus, cette théorie est appliquée dans le problème de contrôlabilité exacte, où la valeur sur la frontière déterminera la valeur intérieure complètement. 1. Introduction Nowadays, quantum studies, especially multiphoton entanglement and interferometry, are attracting many scientists attention, either theoretically or practically[16]. A few world-famous high-tech companies, such as Apple, Microsoft, etc. are developing new generation of high-performance computers based on the quantum mechanics phenomena. In our paper, we discuss an important complex-valued operator in this research field. Let Ax be the vector potential of the magnetic field B, that is, B = A. Clearly, B = div rota = 0. From one of Maxwell s equationsµ is magnetic permeability E = µ B/ t = 0, we deduce that E = φ, where the scalar φ represents the electric potential. We choose an appropriate Lagrangian for the non-relativistic charged particle in the electromagnetic field q is the electric charge of the particle, and v is its velocity, m is mass, L = mv 2 /2 qφ+qv A. Particularly, the canonical Corresponding author: Xiaojun Lu, Department of Mathematics & Jiangsu Key Laboratory of Engineering Mechanics, Southeast University, , Nanjing, China. AMS 2010 subject classification: 35J10, 35J25, 81V10. Key words and phrases: electromagnetic Schrödinger operator, strong unique continuation, multiplier method. 1

2 momentum is specified by the vector p = v L = mv + qa. Next, we define the classical Hamiltonian by Legendre transform, H p v L = p qa 2 /2m+qφ. In quantum mechanics, when p is replaced by i, is the Planck constant, we have the following operator 1 P i +qa 2 /2m+qφ : H H, where H and H are corresponding function spaces. Lots of literature is devoted to the research of this kind of operator[6, 7, 11, 18]. Let Ω R N be an open, connected and bounded domain. From the structure of operator P, we define the corresponding simplified operators 2 H A i +Ax : L 2 Ω L 2 Ω N, 3 H 2 A i +Ax2 : L 2 Ω L 2 Ω, where A C 1 Ω is a real-valued potential vector. The corresponding derivative of the magnetic potential A is as follows, 1 a 1 1 a 2, 1 a N DA = 2 a 1 2 a 2 2 a N..., N a 1 N a 2 N a N where i a j a j / x i,i,j = 1,,N. In addition, one defines the following N N anti-symmetric matrix Ξ A given by ξ 11 ξ 12, ξ 1N ξ Ξ A DA T DA T = 21 ξ 22 ξ 2N... ξ N1 ξ N2 ξ NN with ξ jk j a k k a j, k,j = 1,,N. In quantum mechanics, Ξ A 0 stands for the case without magnetic field, i.e. B = rota = 0. Once the magnetic field exists, then Ξ A 0. Consequently, Ξ A serves as a test matrix for the magnetic field. Interested readers can refer to [5, 14, 15] for more details concerned with the vector operator H A and self-adjoint operator HA 2. In such a manner, 1 is simplified as 4 H 2 A φx : L 2 Ω L 2 Ω, where the complex-valued function φ L Ω. In this paper, we focus on the strong unique continuation propertysucp for the electromagnetic Schrödinger operator 4. In the following, we introduce a few important definitions. 2

3 Definition 1.1. A function u L 2 loc Ω is said to vanish of infinite order at x 0 Ω if for any sufficiently small R > 0, one has 5 u 2 dx = OR M, for every M N +. x x 0 <R Definition 1.2. We say that the operator 4 has SUCP if every solution ω of the equation H 2 A ω = φω, which vanishes of infinite order at x 0 is identically zero in a neighborhood of x 0. So far, the strong unique continuation problem for second order elliptic operators is well-understood. In the case of Ω = R 2, Carleman proved the SUCP of the elliptic equation with bounded coefficients and V L loc R2 6 u = W u+vu by introducing a weighted L 2 -estimate, the so-called Carleman estimate [4]. For the space dimension N 3 with bounded coefficients, N. Aronszajn, A. Krzywicki and J. Szarski proved the SUCP by means of Carleman type inequalities, namely, observability inequalities. Afterwards, D. Jerison, C. E. Kenig, C. D. Sogge treated the equation 6 with singular potentials V L N/2 loc RN and W L R N, N 3, by the approach of L p L q Carleman estimate involving sharp exponents[10, 11, 17]. Afterwards, N. Garofalo and F. H. Lin gave a new proof for the SUCP of the elliptic operator u = Vu with bounded potential by applying a variational method in [9]. There is a large body of work on SUCP for 6 with real-valued coefficients. In this paper, we investigate the complex-valued case. As a matter of fact, the operator H 2 A can be decomposed into 7 H 2 Aω = ω +ia ω +i Aω+AA T ω. In [12, 13], K. Kurata proved the SUCP for 4 with AA T K loc N Ω, where K loc N Ω denotes the Kato class. When the potential A L Ω N, in effect, it does not belong to the Kato class. As a result, we can not deduce corresponding results directly from K. Kurata s work. In this manuscript, we intend to provide a new approach of SUCP for 4 with complex-valued coefficients by developing new multipliers. At the moment one is ready to state the main results. Theorem 1.3. For N 2, let the complex-valued ω H 2 B 1 be a solution of the problem 8 ω +ia ω +i Aω+AA T ω = φxω in B 1, where B 1 is a unit ball B 1 Ω, A C 1 Ω is a real-valued potential vector and the complex-valued function φ L R N. If ω vanishes of infinite order at x 0 B 1, then ω 0 in B 1. By virtue of Theorem 1.3, one is able to prove the following statement for a mixed boundary value problem which is of great importance in the discussion of exact controllability through boundary control [14]. 3

4 Corollary 1.4. Assume that Ω is a bounded, open and connected domain in R N with the boundary Γ C 2, A C 1 Ω is a real-valued potential vector and the complexvalued function φ L R N. Let ω H 2 Ω be the solution of the mixed boundary problem ω +ia ω +i Aω+AA T ω = φxω in Ω, Then ω is identically 0 in Ω. ω = ω/ ν = 0 on Γ. Remark 1.5. Theorem 1.3 demonstrates, the asymptotic behavior of the solution ω at an interior point x 0 determines the interior value of ω in B 1. In contrast with Theorem 1.3, Corollary 1.4 indicates, the behavior of solution ω on the boundary determines the interior value of ω in Ω. The rest of the paper is organized as follows. First and foremost, in Section 2, we introduce some useful quantities and their particular properties. Next, we give an important comparison lemma and a frequency function. By carefully estimating the derivative of the frequency function, we reach the conclusion in the final analysis. In Section 3, as an important application in exact controllability, we prove Corollary 1.4 in detail. 2. Proof of the main theorem: A new multiplier method First, we introduce several quantities which will serve as useful tools for our purposes. For every r 0,1, we define the following two quantities 9 Φr ω 2 ds x, where is centered at the origin with radius r, denotes its sphere, ds x stands for the N-1-dimensional Hausdorff measure on the sphere. 10 Ψr H A ω 2 φ R ω 2 dv x, where φ R denotes the real part of φ. Actually, we have Lemma 2.1. By virtue of divergence theorem, the following identity holds, 11 Re ω 2 ia ω 2 x/rds x = 2 H A ω 2 +2φ R ω 2 dv x. Proof. On the one hand, HA 2 ω 2 dv x = ω 2 ia ω 2 x/rds x + A H A ω 2 dv x 12 = ω 2 ia ω 2 x/rds x + A iω ω +iω ω +A ω 2 dv x. 4

5 On the other hand, 13 Since H 2 A ω 2 dv x = ω ω ω ω 2 ω 2 +i A ω 2 dv x + iωa ω +iωa ω +AA T ω 2 dv x. A ω 2 = ω Aω +ωa ω = ω Aω +ωa ω, then by combining 12 and 13, we have Re ω 2 ia ω 2 x/rds x = ω ω ω ω 2 ω dv 2 x = 2 ω 2 +2iωA ω 2iωA ω 2AA T ω dv 2 x + ω ω +iωa ω +iω Aω +AA T ω dv 2 x + ω ω iωa ω iω Aω +AA T ω dv 2 x = 2 H A ω 2 +2φ R ω dv 2 x. Next we calculate the derivatives of Φr and Ψr with respect to r. Lemma 2.2. The derivatives of Φr and Ψr with respect to r are presented as follows, 14 Φ r = N 1Φr/r+2Ψr. 15 Ψ r = N 2Ψr/r+N 2/r φ R ω 2 dv x +2/r Re Aω H A ωdv x +2/r Re x ω φωdv x +2/r Re ωxda T H A ω T dv x +2 ν i ω +Aω 2 ds x 2 Re Aω νh A ω νds x φ R ω 2 ds x. Remark shows that Φr and Ψr are closed related with each other. This relation is very important for our discussion. 5

6 Proof. First we consider the derivative of Φr with respect to r. Indeed, we have Φ r = B 1 ωry 2 r N 1 r ds y = B 1 ω yω+ω ω yr N 1 + ω 2 N 1r N 2 ds y = ω x/rω +ω ω x/r ds x +N 1/r ω 2 ds x = N 1/rΦr+Re ω 2 ia ω 2 x/rds x = N 1Φr/r+2Ψr. As for the derivative of Ψr with respect to r, via the divergence theorem, we have Ψ r = H A ω 2 ds x φ R ω 2 ds x = 1/r H A ω 2 x x/rds x φ R ω 2 ds x 16 = 1/r div H A ω 2 xdv x φ R ω 2 ds x = N/r H A ω 2 dv x +1/r x H A ω 2 dv x φ R ω 2 ds x. B } r {{} I 6

7 Now we treat the term I carefully. I = 1/r x j j i k ω +a k ωi k ω +a k ω dv x j,k = j,k = j,k = j,k j,k + j,k = j,k j,k j,k + j,k 2/r Re x j j i k ω +a k ωi k ω +a k ω dv x 2/r Re ix j j k ωi k ω +a k ωdv x + j,k 2/r Re ix j j ωi k ω +a k ων k ds x j,k 2/r Re ix j j ω k i k ω +a k ωdv x + j,k 2/r Re a k x j j ωi k ω +a k ωdv x 2/r Re x j i j ω +a j ωi k ω +a k ων k ds x 2/r Re x j a j ωi k ω +a k ων k ds x j,k 2/r Re ix j j ω k i k ω +a k ωdv x + j,k 2/r Re a k x j j ωi k ω +a k ωdv x 2/r Re x j j a k ωi k ω +a k ωdv x 2/r Re i k x j j ωi k ω +a k ωdv x 2/r Re x j ω j a k i k ω +a k ωdv x 2/r Re i k x j j ωi k ω +a k ωdv x 2/r Re x j ω j a k i k ω +a k ωdv x = 2 ν H A ω 2 ds x 2 Re Aω νh A ω νds x 2/r H A ω 2 dv x +2/r Re Aω H A ωdv x +2/r Re x ω H 2 A ωdv x +2/r Re ωxda T H A ω T dv x = 2 ν H A ω 2 ds x 2 Re Aω νh A ω νds x 2/r H A ω 2 dv x +2/r Re Aω H A ωdv x +2/r Re x ω φωdv x +2/r Re 7 ωxda T H A ω T dv x.

8 Finally, keeping in mind the definition of Ψr, we reach the conclusion immediately. Next we show an important comparison lemma. Lemma 2.4. There exists an r 0 0,1 such that for every r 0,r 0, we have 17 ω 2 dv x r ω 2 ds x. Proof. On the one hand, H 2 A ω 2 r 2 x 2 dv x = ω 2 H 2 A r2 x 2 dv x + ω 2 r 2 x 2 / ν iha ds x = ω 2N 2 +2iA x i Ar 2 x 2 +AA T r 2 x 2 dv x On the other hand, 2r ω 2 ds x. H 2 A ω 2 r 2 x 2 dv x = ω ω ω ω 2 ω 2 r 2 x 2 dv x + AA T ω 2 r 2 x 2 dv x + i A ω 2 +iωa ω +iωa ω r 2 x 2 dv x = 2 H A ω 2 +2φ R ω r 2 2 x 2 dv x + AA T ω 2 r 2 x 2 dv x + i A ω 2 +iωa ω +iωa ω r 2 x 2 dv x. As a result, 2N ω 2 +2 H A ω 2 r 2 x 2 2φ R ω 2 r 2 x 2 dv x = 2r ω 2 ds x. When φ R L > 0, then we choose r 0 0,1/2 such that r 2 0 N 1/ φr L. It follows immediately that ω 2 dv x N ω 2 φ R ω 2 r 2 x 2 dv x r ω 2 ds x. When φ R L = 0, then it is evident ω 2 dv x r/n ω 2 ds x. 8

9 Assume that there exists a small r 1 0,1 such that 18 Φr 0 for r 0,r 1. Define the frequency function 19 r rψr/φr, r 0,r 1. Let r min{r 0,r 1 }, and we set 20 ℶ r { } r 0,r : r > 1. With the above definitions, we have the following inequality for the frequency function. Lemma 2.5. Under the assumptions 18-20, there exists a positive constant τ = τn,φ which is independent of r such that r is estimated in a uniform fashion, r rτ. Proof. Actually, from 17-20, we have H A ω 2 dv x > 1/r 2 φ R L ω 2 dv x. Indeed, H A ω 2 dv x = H A ω 2 φ R ω 2 dv x + φ R ω 2 dv x > 1/r ω 2 ds x + φ R ω 2 dv x 1/r 2 ω 2 dv x + φ R ω 2 dv x 1/r 2 φ R L ω 2 dv x. 9

10 This indicates the integral H A ω 2 dv x is the dominating part in Ψr. By calculating r with respect to r, we have the following identity, r = r Ψ r/ψr+1/r Φ r/φr = r 2 ν i ω +Aω 2 ds x /Re x/r ωωds x B }{{ r } II 2 Re x/r ωωds x / ω 2 ds x B }{{ r } II { + r N 2/r +2/r Re φ R ω 2 dv x B }{{ r } III x ω φωdv x B } r {{} V 2 Re Aω νh A ω νds x B } r {{} VII +2/r Re Aω H A ωdv x B } r {{} IV +2/r Re ωxda T H A ω T dv x B } r {{} VI φ R ω 2 ds x }{{} V III }/{1/2 Br x/r ω 2 ds x }. 10

11 We estimate each term respectively. For II-II, we apply Hölder s inequality and obtain II II = 2 ν i ω +Aω 2 ds x /Re x/r ωωds x 2 Re x/r ωωds x / ω 2 ds x = 2 ν i ω +Aω 2 ds x /Re x/r ω iaωωds x 2 Re x/r ω iaωωds x / ω 2 ds x 2 ν i ω +Aω 2 ds x / x/r ω iaω 2 ds Br x ω 2 ds Br x 0. 2 x/r ω iaω 2 ds x ω 2 ds x / ω 2 ds x In addition, we have Lemma 2.6. There exists a constant C φ independent of r such that ν i ω +Aω 2 ds x /Ψr C φ/r. Proof. Indeed, by multiplying Hx H A ω to HA 2 ω = φω and integrating by parts, we have the following identity, 1/2 ω/ νiha 2 Hx ν ds x = 1/2 Hx H A ω 2 dv x Im φω Hx H A ω dv x Re H A ωdh T H T A ωdv x Re ωh A ωξ A H T dv x. Since H A ω 2 dv x 1/r 2 φ R L ω 2 dv x, by choosing Hx x/r, 11

12 we have the following estimate, 1/2 ω/ νiha 2 dsx N/2r H A ω 2 dv x +1/2 φ L ω 2 + H A ω 2 dv x = +1/r H A ω 2 dv x +1/2max Ξ A F ω 2 + H A ω 2 dv x N +2/2r+1/2 φ L +max Ξ A F H A ω 2 dv x }{{} α +1/2 φ L +max Ξ A F ω 2 dv x }{{} β α+βr 2 /1 r 2 φ R L H Aω 2 dv x. where Ξ A F denotes the Frobenius norm of the test matrix Ξ A. Since Ψr = H A ω 2 φ R ω 2 dv x 1 2r 2 φ R L /1 r 2 φ R L H A ω 2 dv x, therefore, ν i ω +Aω 2 ds x /Ψr 1 r 2 φ R L /1 2r 2 φ R L 2α+2r 2 /1 r 2 φ R L β. The conclusion follows immediately. Taking Lemma 2.6 into account and noticing the fact r > 1, for the term III, we have III/Ψr = N 2 φ R ω 2 dv x /rψr N 2 φ R L Φr/Ψr rn 2 φ R L. 12

13 For the term IV, IV/Ψr = 2 Re Aω H A ωdv x /rψr 2 Aω H A ω dv x /rψr = 1/2ǫr Br Aω 2 dv x +2ǫ/r Br H A ω 2 dv x /Ψr 12ǫr Br AA T ω 2 dv x +2ǫ/r Br H A ω 2 φ R ω 2 dv x +2ǫ/r Br φ R ω 2 dv x /Ψr AA T L Φr/2ǫΨr+2ǫ/r+2ǫ φ R L Φr/Ψr. Let ǫ = r/2, since r > 1, then IV/Ψr AA T L +1+r 2 φ R L. For the term V, V/Ψr = 2 Re x ω φωdv x /rψr = 2 Re x ω iaω φωdv x /rψr = φ L /2ǫr ω 2 dv x +2/rǫ φ L x H A ω 2 dv x /Ψr φ L /2ǫr ω 2 dv x +2rǫ φ L H A ω 2 dv x /Ψr φ L /2ǫr ω 2 dv x +2rǫ φ L H A ω 2 φ R ω 2 dv x +2rǫ φ L φ R ω 2 dv x /Ψr r φ L /2ǫ+2rǫ φ L +2r 3 ǫ φ 2 L. Let ǫ = r/2, then V/Ψr φ L 1+r 2 +r 4 φ L. 13

14 For the term VI, V I/Ψr = 2 Re ωxda T H A ω T dv x /rψr 2 ωx 2 DA T H A ω T 2 dv x /rψr 2 ωx 2 DA T F H A ω 2 dv x /rψr 2max DA T F ω H A ω dv x /Ψr = max DA T F /2ǫ ω 2 dv x +2ǫmax DA T F { H A ω 2 φ R ω 2 dv x + φ R ω 2 dv x }/Ψr max DA T F r 2 /2ǫ+2ǫ+2ǫr 2 φ R L. Let ǫ = r/2, then VI/Ψr max DA T F 2r +r 3 φ R L, where DA T F denotes the Frobenius norm of DA T. Let ǫ = r/2. By Schwartz s inequality, we have the estimate below for the term VII, V II/Ψr = 2 Re Aω νh A ω νds x /Ψr 2 Aω ν H A ω ν ds x /Ψr 1/2ǫ Br Aω 2 ds x +2ǫ Br H A ω ν 2 ds x /Ψr AA T L ω 2 ds x /2ǫΨr+2ǫ H A ω ν 2 ds x /Ψr r AA T L /2ǫ+2ǫ H A ω ν 2 ds x /Ψr AA T L +C φ, where C φ is from Lemma 2.6. For the last term VIII, a simple calculation leads to V III/Ψr = φ R ω 2 ds x /Ψr r φ R L. 14

15 Fromtheaboveestimates, weconcludethatthereexists apositiveconstant τ = τn,φ which is independent of r such that r rτ. It follows that expτr r is monotonously increasing on 0,r, that is to say, expτr r expτr r. Keeping in mind the case 1, we know that, r is bounded on 0,r. Since Φ r = N 1/rΦr+2Ψr, then logφr/r N 1 = 2Ψr/Φr = 2 r/r Cτ/r. We integrate from γ to 2γ, then log2 1 N Φ2γ/Φγ Cτlog2. It follows that Φ2γ 2 Cτ+N 1 Φγ. Finally, integrating with respect to γ gives ω 2 dv x 2 Cτ+N 1 ω 2 dv x. B 2γ B γ Since B 1 is connected, then our theorem follows immediately. Remark 2.7. It is of great interest to explore the strong unique continuation for a variety of Schrödinger operators with singular or nonlinear potentials by the multiplier method. More results will be available in sequential papers. Set 3. Proof of Corollary 1.4 In this section, we show an important application of Theorem 1.3 in [14]. Proof of Corollary 1.4: Let B be an arbitrarily small open ball such that and define Γ B. Ω 1 Ω B, ω in Ω; ω 1 0 in B\Ω. It is sufficient to verify that ω 1 H 2. Denote by ωj 1, ω1 jk the extension by zero to Ω 1 of the derivatives j ω, j k ω, j,k = 1,,N. Then ω j, ω jk L 2 Ω 1 and it is necessary to demonstrate that, for ζ DΩ 1, ω 1 j ζdx = ωj 1 ζdx, Ω 1 Ω 15 1

16 and ωj 1 k ζdx = ωjkζdx. 1 Ω 1 Ω 1 Indeed, since ωj 1 = ω1 jk 0 outside of Ω, ζ 0 on Γ\Γ B and ω = ω/ ν 0 on Γ B, we have ω 1 Ω 1 j ζdx = ω Ω jζdx = ωζν Γ jdγ Ω jωζdx and = Γ B ωζν jdγ Ω jωζdx = Ω jωζdx = Ω 1 ω 1 j ζdx, Ω 1 ω 1 j kζdx = Ω jω k ζdx = Γ jωζν k dγ Ω k j ωζdx = Γ B jωζν k dγ Ω k j ωζdx = Ω k j ωζdx = Ω 1 ω 1 jk ζdx. Thus, the result is concluded due to the connectness of Ω. Acknowledgement This project is partially supported by Natural Science Foundation of Jiangsu Province BK , National Natural Science Foundation of ChinaNSFC , , , the Scientific Research Foundation for the Returned Overseas Chinese Scholars, Open Research Fund Program of Jiangsu Key Laboratory of Engineering Mechanics, Southeast University LEM16B06, Fundamental Research Funds for the Central Universities on the Field Research of Commercialization of Marriage between China and Vietnam No. 2014B In particular, the authors also express their deep gratitude to the referees for their careful reading and useful remarks. References [1] R. Adams, Sobolev Spaces, Academic Press, New York, [2] N. Aronszajn, A unique continuation theorem for solutions of elliptic partial differential equations or inequalities of second order. J. Math. Pures Appl. 9, , [3] N. Aronszajn, A. Krzywicki and J. Szarski, A unique continuation theorem for exterior differential forms on Riemannian manifold. Ark. Mat., 41996, [4] T. Carleman, Sur un problème d unicité pour les systèmes d équations aux dérivéés partielles à deux variables indépendantes. Ark. Mat., 26B1939, 1-9. [5] L. Fanelli, Dirpersive equations in quantum mechanics, doctoral thesis, University of Rome, preprint. [6] L. Fanelli, Electromagnetic Schrödinger flow: multiplier methods for dispersion, Proc. Journées EDP, Port D Albret, GDR 2434CNRS, [7] M. Goldberg, Strichartz estimates for Schrödinger operators with a non-smooth magnetic potential, Discrete and Continuous Dynmaical Systems, , [8] N.GarofaloandF.H. Lin, Monotonicitypropertiesofvariationalintegrals,A p weightsandunique continuation, Indiana Univ. Math. J., , [9] N. Garofalo and F. H. Lin, Unique continuation for elliptic operators: A geometric variational approach, Comm. Pure Appl. Math., , [10] D. Jerison and C. E. Kenig, Unique continuation and absence of positive eigenvalues for Schrödinger operator. Ann. Math ,

17 [11] H. Koch and D. Tataru, Recent results on unique continuation for second order elliptic equations. Carleman estimates and applications to uniqueness and control theory, Progr. Nonlinear Differential Equations Appl , [12] K. Kurata, A unique continuation theorem for the Schrödinger equation with singular magnetic field, Proc. AMS, , [13] K. Kurata, Local boundedness and continuity for weak solutions of ib 2 u+vu = 0, Math. Z., , [14] X. Lu, Z. Tu and X. Lv, On the exact controllability of hyperbolic magnetic Schrödinger equations, Nonlinear Analysis , [15] N. Arrizabalaga and M. Zubeldia, Unique continuation for magnetic Schrödinger operator with singular potentials, arxiv: v1, [16] J. Pan, Z. Chen, C. Lu, H. Weinfurter, A. Zeilinger and M. Zukowski, Multiphoton entanglement and interferometry, Reviews of Modern Physics , [17] C. D. Sogge, Strong uniqueness theorems for second order elliptic differential equations. Amer. J. Math., , [18] Michael E. Taylor, Partial Differential Equations II, Qualitative Studies of Linear Equations, Springer Verlag,