Stable determination of the first order perturbation of the biharmonic operator from partial data

기술 요약: 부분 데이터를 이용한 비조화 연산자의 1차 섭동의 안정적 결정

문제 정의
본 논문은 경계가 매끄러운 유계 영역 $\Omega \subseteq \mathbb{R}^n$ ($n \ge 3$)에서 섭동된 비조화 연산자(biharmonic operator)에 대한 역 경계값 문제를 다룬다. 연산자는 $L_{A,q} = \Delta^2 + A(x) \cdot D + q(x)$로 정의되며, 여기서 $D = i^{-1}\nabla$, $A$는 $W^{1,\infty}(\Omega, \mathbb{C}^n)$에 속하는 벡터장, $q$는 $L^\infty(\Omega, \mathbb{C})$에 속하는 퍼텐셜이다. 본 연구는 부분 데이터(partial data) 설정에 초점을 맞춘다. 즉, 임의의 작은 열린 부분 집합 $\gamma_1, \gamma_2 \subset \partial\Omega$에서 측정된 데이터를 통해 계수 $A$와 $q$를 복원하는 것이다.

이 문제는 계수들이 경계 근방 $\omega_0$에서 사전적으로 알려져 있다고 가정한다. 데이터는 $\gamma_1$에서 지지되는 경계 데이터 $(f, g)$로부터 $(\partial_\nu u, \partial_\nu(\Delta u))$를 측정하는 부분 디리클레-투-뉴먼(Dirichlet-to-Neumann, DtN) 맵 $\Lambda_{\gamma_1, \gamma_2}^{A,q}$에 의해 제공된다. 주요 목표는 1차 섭동 $A$와 0차 섭동 $q$ 모두에 대한 안정성 추정치(stability estimates)(측정 오차가 계수 복원에 미치는 영향을 정량화하는 것)를 확립하는 것이다.

방법론
저자들은 고차 타원형 연산자의 역 문제에서 표준적인 기법인 복소 기하 광학(Complex Geometric Optics, CGO) 해법을 사용한다. 증명 전략은 다음과 같은 몇 가지 핵심적인 기술적 구성 요소들을 포함한다:

1. CGO 해의 구성: 저자들은 $u(x) = e^{ix\cdot\zeta/h}(a(x) + r(x))$ 형태를 갖는 수반 방정식(adjoint equation)과 순방향 방정식(forward equation)의 해 $u_1, u_2$를 구성한다. 여기서 $\zeta$는 $\zeta \cdot \zeta = 0$을 만족하는 복소 벡터이고, $h$는 준클래식(semiclassical) 파라미터이며, $r(x)$는 $h \to 0$일 때 소멸하는 나머지 항(remainder term)이다.
2. 적분 항등식 및 그린 정리(Green's Formula): 두 연산자의 차이에 그린 정리를 적용하고 부분 DtN 맵을 활용함으로써, 계수의 차이($A_1 - A_2$ 및 $q_1 - q_2$)와 DtN 맵의 차이를 연결하는 적분 항등식을 유도한다.
3. 국소 정량적 유일 연속성(Local Quantitative Unique Continuation): 본 연구의 핵심적인 혁신은 균일한 경계 조건을 가진 비균질 비조화 방정식의 해에 대한 국소 정량적 유일 연속성 추정치를 도출하는 것이다. 이는 준클래식 비조화 연산자 $h^4\Delta^2$에 대한 **카를레만 추정치(Carleman estimates)**를 사용하여 달성된다. 이 단계는 저자들이 전체 경계가 아닌 측정 부분 집합 $\gamma_2$ 상에서의 해 노름을 제어할 수 있게 해주며, 이는 부분 데이터 설정에서 필수적이다.
4. 벡터장의 분해: 외미분 $dA$로부터 벡터장$A$를 복원하기 위해, 저자들은 헬름홀츠 유형의 분해$A = A_{sol} + \nabla \phi$를 활용한다. 즉,$dA$의 유계(bound)를 사용하여 솔레노이달(solenoidal) 부분 $A_{sol}$을 추정하고, $\phi$에 대한 경계값 문제를 해결한 뒤 소볼레프 임베딩 및 보간 부등식을 적용하여 퍼텐셜 부분 $\nabla \phi$를 추정한다.
5. 파라미터 최적화: CGO 나머지 항과 카를레만 추정치의 지수 가중치로부터 발생하는 오차 항들 사이의 균형을 맞추기 위해, 준클래식 파라미터 $h$와 푸리에 영역에서의 주파수 컷오프 $\rho$를 최적화함으로써 안정성 추정치를 도출한다.

주요 기여 및 결과
본 논문은 역 문제가 갖는 특성인 로그 유형(logarithmic-type) 및 이중 로그 유형(double-logarithmic-type) 안정성 추정치를 확립한다.

• 1차 섭동($A$)에 대한 안정성: 정리 1.1은 벡터장 $A$에 대한 로그 유형 안정성 추정치를 확립한다. 구체적으로, $A_1 - A_2$의 $L^\infty$ 노름은 DtN 차이의 거듭제곱과 DtN 차이의 음의 로그 거듭제곱의 합으로 유계된다:
  $$\|A_1 - A_2\|_{L^\infty} \lesssim \|\Lambda_1 - \Lambda_2\|^{\mu_1} + |\log \|\Lambda_1 - \Lambda_2\||^{-\mu_2}$$
  이 결과는 기존의 유일성 결과들을 보완하며 정량적인 안정성 경계치를 제공한다.

• 1차 항이 존재하는 경우 0차 섭동($q$)에 대한 안정성: 정리 1.2는 1차 항 $A$가 존재할 때 퍼텐셜 $q$에 대한 로그-로그 유형 안정성 추정치를 확립한다:
  $$\|q_1 - q_2\|_{H^{-1}} \lesssim |\log |\log \|\Lambda_1 - \Lambda_2\|||^{-\mu'}$$
  따로 lệ 1.3은 $q$에 대한 추가적인 정칙성 가정이 있을 때 이를 $L^\infty$ 추정치로 확장한다.

• 1차 섭동이 없는 경우 0차 섭동($q$)에 대한 안정성: 정리 1.4는 1차 섭동($A=0$)이 없을 때 $q$에 대한 안정성이 로그-로그 유형에서 로그 유형으로 개선됨을 보여준다:
  $$\|q_1 - q_2\|_{H^{-1}} \lesssim \|\Lambda_1 - \Lambda_2\|^{2/3} + |\log \|\Lambda_1 - \Lambda_2\||^{-2/(n+2)}$$
  이는 비조화 연산자의 1차 항이 유발하는 특정한 난이도를 강조한다.

의의 및 주장
저자들은 본 연구가 부분 데이터를 사용하는 비조화 연산자의 1차 섭동에 대한 최초의 안정성 결과를 제공한다고 주장한다. 유일성에 대해서는 이미 기존 연구([42])에서 확립되었으나, 안정성 분석은 훨씬 더 까다롭다.

본 논문은 두 가지 주요 기술적 참신함을 강조한다:

1. 국소 유일 연속성: 카를레만 추정치를 사용하여 고차 타원형 연산자(특히 비조화 연산자)에 대한 국소 정량적 유일 연속성 추정치를 도출하였다. 저자들은 이것이 특정 부분 데이터 맥락에서 고차 연산자에 대해 이루어진 첫 번째 결과임을 명시한다.
2. 부분 데이터 처리: 특정 기하학적 구성(예: 평평한 접근 불가능 영역)이나 전체 경계 데이터에 의존하지 않고, 경계의 임의의 작은 열린 부분 집합으로부터의 측정만을 사용하여 안정성 추정치를 도출할 수 있다.

저자들은 자신의 결과가 기존의 유일성 정리들을 보완하는 동시에, 슈뢰딩거 및 자기 슈뢰딩거 연산자에 사용된 방법론과 유사성을 가지면서도 비조화 연산자의 구조 및 1차 섭동의 게이지 불변성 결여가 초래하는 독특한 과제들을 해결했음을 밝히며, 이를 통해 고차 역 문제의 안정성 분석을 진전시켰다고 위치시킨다.