부분 편광과 Stokes 관측 재구성

측정된 빔에 하나의 타원이 필요하지는 않다

검출기는 여러 편광 상태로 준비한 빛을 평균할 수 있다. 각 실현을 하나의 Jones 벡터로 나타낼 수 있어도 모든 실현을 동시에 설명하는 벡터 하나가 존재하지 않을 수 있다. 그런 빔을 결정적 타원 하나로 그리면 관측은 맞아 보여도 물리적 설명이 틀릴 수 있다.

이 단원에서는 제곱 관측량을 평균하여 결맞음 행렬을 만들고 분석기 강도에서 Stokes 성분을 복원한다. 순수 상태·부분 편광·무편광·빛이 없는 경우를 구분한다. 또한 한 출력 상태가 유효하다는 것과 실제로 구성할 수 있는 광학 변환이라는 더 강한 조건을 분리한다.

확률 집단과 정규화 관측의 가정

앞 Jones 단원과 같은 오른손 횡방향 기저 및 $e^{-i\omega t}$ 관례를 쓴다. 유한 집단은 벡터 $\mathbf E_r$와 가중치 $w_r\ge0$로 구성하며 전체 가중치는 양수여야 한다. 확률을 $p_r=w_r/\sum_sw_s$로 정규화한다.0인 가중치는 허용하지만 모든 가중치가0인 경우에는 평균을 정의하지 않는다.

이것은 검출 평균을 위한 명시적 준비 분포 또는 표본 모델이다. 미시적 시간 변동·공간 상관·광자 통계를 모두 지정하지는 않는다. 같은 행렬 대수가 다른 통계 이론에도 나타나지만 고전적 이차 모멘트 계산만으로 양자광 모델을 검증한 것은 아니다.

정규화 강도는 복소장 제곱 크기의 집단 평균이다. 실제 SI 복사조도에는 매질 교정과 실수장 주기 평균 계수가 여전히 필요하다. 정규화 전후에서 영강도와 양의 강도의 무편광을 서로 대체하지 않는다.

관측량에 앞서 결맞음 행렬 구성하기

$$C=\sum_rp_r\mathbf E_r\mathbf E_r^\dagger =\begin{pmatrix}C_{xx}&C_{xy}\\C_{xy}^*&C_{yy}\end{pmatrix}$$

로 정의한다. dagger는 복소켤레 전치다. 임의 분석기 벡터 $\mathbf a$에 대해

$$\mathbf a^\dagger C\mathbf a =\sum_rp_r|\mathbf a^\dagger\mathbf E_r|^2\ge0$$

이므로 $C$는 구성 자체로 Hermitian 양의 준정부호 행렬이다. 대각합은 전체 정규화 강도다. 적합한 행렬을 나중에 잘라 물리적으로 보이게 만든 것이 아니라 실제 비음수 관측량을 평균해서 이 성질을 얻는다.

수평 장 $(1,0)$과 $(-1,0)$을 같은 비율로 섞으면 평균 장은0이지만 $C=\operatorname{diag}(1,0)$이다. 평균 장을 제곱하면 가짜 어둠을 예측하고, 실현 강도를 평균하면1이다. 진폭 평균과 검출 평균은 순서를 바꿀 수 없는 다른 연산이다.

Stokes 부호와 물리적 원뿔

정한 관례에서

$$S=(S_0,S_1,S_2,S_3) =(C_{xx}+C_{yy},C_{xx}-C_{yy},2\operatorname{Re}C_{xy},-2\operatorname{Im}C_{xy})$$

다. 역관계는

$$C=\frac12\begin{pmatrix}S_0+S_1&S_2-iS_3\\S_2+iS_3&S_0-S_1\end{pmatrix}$$

이며, $R=\sqrt{S_1^2+S_2^2+S_3^2}$라 두면 고유값은 $\lambda_\pm=(S_0\pm R)/2$다. 물리적 상태는 $S_0\ge0$, $R\le S_0$를 만족한다. 각 성분이 그럴듯한 범위에 있어도 이 원뿔 밖의 벡터는 양의 결맞음 행렬이 될 수 없다.

$(1,i)/\sqrt2$이면 $C_{xy}=-i/2$이므로 $S_3=+1$이다. 문헌이나 소프트웨어를 바꿀 때 시간 인자·원편광 관측 방향·이 부호를 함께 대조한다. 이름만 같다고 서로 다른 관례의 성분을 그대로 연결하지 않는다.

편광도와 정규화된 공

$S_0>0$에서 전체 편광도와 선형·원형 성분을

$$P=\frac{R}{S_0},\qquad P_L=\frac{\sqrt{S_1^2+S_2^2}}{S_0},\qquad P_C=\frac{S_3}{S_0}$$

로 정의한다. 여기서 $P_C$는 부호를 갖지만 $P$는 비음수다. 정규화 벡터 $(S_1,S_2,S_3)/S_0$는 단위 공 안에 있다. 영이 아닌 결정적 Jones 상태는 표면, 부분 편광은 내부, 무편광은 중심에 놓인다. 공의 중심과 빛이 없는 경우는 다르다.

$S_0=0$이면 물리적인 모든 Stokes 성분도0이며 정규화 상태와 편광도는 미정이다. 편의상 $P=0$을 넣으면 신호가 없는 것을 측정 가능한 무편광과 혼동한다. 원형이나 거의 무편광인 상태에서는 일부 방향 각도도 불안정하다. 이 경우에도 정규화하지 않은 행렬은 유용한 관측 자료로 유지한다.

편광 성분과 무편광 성분의 분해

$R>0$일 때 $C$의 주 고유벡터를 단위 벡터 $\mathbf u$라 두면

$$C=\frac{S_0-R}{2}I_2+R\mathbf u\mathbf u^\dagger$$

다. 두 강도는 $I_{\rm unpol}=S_0-R$, $I_{\rm pol}=R$이며 $P=I_{\rm pol}/S_0$를 설명한다. $R>0$에서는 주 방향을 전체 위상까지 제외하고 정할 수 있지만 두 고유값이 같아지면 그 방향도 사라진다.

이 분해가 유일한 미시적 집단을 찾아 주지는 않는다. 수평·수직의 같은 비율 혼합과 대각·반대각의 같은 비율 혼합은 모두 $C=I_2/2$다. 이차 편광 관측은 같아도 시간에 따라 준비된 과정까지 같다는 뜻은 아니다. 분해한 유효 편광 성분을 부분 편광 전체의 유일한 실제 전기 벡터로 그리지 않는다.

결정적 소자의 결맞음 행렬 변환

Jones 행렬 $J$에 대해

$$C'=JCJ^\dagger$$

다. 각 실현을 변환한 뒤 외적을 평균하면 바로 유도된다. 유니터리 지연판은 대각합과 고유값을 보존하므로 편광도도 유지한다. 손실 투영에서는 전체 강도와 고유값 비율이 일반적으로 모두 달라진다.

이상적 수평 분석기는 $C=I_2/2$를 $\operatorname{diag}(1/2,0)$으로 보낸다. 강도는1에서1/2로 감소하지만 투과분의 편광도는0에서1로 증가한다. 직교 성분을 제거하여 선택한 것이지 전력을 만들거나 입사 빔 전체를 같은 강도의 순수 상태로 바꾼 것은 아니다. 정규화 전 투과량과 조건부 정규화 출력을 함께 보고한다.

물리적 혼합으로 탈편광 구성하기

간단한 등방 채널은

$$\mathcal D_d(C)=dC+(1-d)\frac{\operatorname{Tr}C}{2}I_2, \qquad 0\le d\le1$$

다. 강도는 유지하고 편광 벡터를 $d$배로 줄인다. 항등 변환의 가중치를 $(1+3d)/4$, 세 Pauli 켤레 변환의 가중치를 각각 $(1-d)/4$로 두면 비음수 혼합으로 구성할 수 있다. 구면의 벡터를 임의로 줄인 그림보다 강한 물리적 구성 근거다.

채널의 위치도 중요하다. 구현에서는 이 혼합이 결정적 광학 소자보다 먼저 작용한다. 단위 강도의 수평 입력을 완전히 탈편광한 뒤 수평으로 거르면 강도1/2·편광도1이다. 먼저 수평으로 거른 뒤 탈편광하면 강도1·편광도0이다. 탈편광과 손실은 다른 작용이며 순서는 모델의 일부다.

Jones 변환에서 Mueller 행렬로

Hermitian 기저를 $(\sigma_0,\sigma_1,\sigma_2,\sigma_3)=(I_2,Z,X,Y)$로 두면 $C=\sum_bS_b\sigma_b/2$다. 변환 뒤 대각합을 취하여

$$S'_a=\sum_bM_{ab}S_b,\qquad M_{ab}=\tfrac12\operatorname{Tr}(\sigma_aJ\sigma_bJ^\dagger)$$

를 얻는다. Jones 행렬은 복소수여도 이 성분들은 실수다. Mueller 표현은 정한 제곱 관측량을 변환하고 Jones의 전체 위상은 잃는다. 무작위 소자 집단에서는 각 Mueller 변환 또는 변환한 결맞음 행렬을 평균해야 한다. Jones 행렬부터 평균하면 보통 다른 물리 과정이 된다.

위 등방 채널의 Mueller 행렬은 $\operatorname{diag}(1,d,d,d)$다. 탈편광이 먼저이면 결정적 Mueller 행렬의 오른쪽에 이 대각 행렬을 곱한다. 관측량의 공간으로 옮겨도 소자 순서의 의미가 사라지지 않는다.

여섯 분석기로 Stokes 복원하기

단위 분석기를 $H=(1,0)$, $V=(0,1)$, $D=(1,1)/\sqrt2$, $A=(1,-1)/\sqrt2$, $C_\pm=(1,\pm i)/\sqrt2$로 둔다. 강도 $I_a=\mathbf a^\dagger C\mathbf a$에서

$$S_0=I_H+I_V,\quad S_1=I_H-I_V,\quad S_2=I_D-I_A,\quad S_3=I_{C+}-I_{C-}$$

를 복원한다. 이상적으로 교정된 모델에서는 세 상보 쌍의 합이 같은 전체 강도여야 한다. 실제 실험에서 불일치하면 시간 변동·검출기 응답 차이·잡음·준비 모델의 오류를 점검한다. 서로 다른 교정 단위의 강도를 그대로 빼면 성분 해석도 틀릴 수 있다.

계산에서는 실제 실현마다 양의 분석기 강도를 구해 평균한다. 거의 어두운 채널에서 전개한 큰 항끼리 빼다가 가짜 음의 전력이 생기는 것을 피한다. 복원한 Stokes는 측정 결과로 보존하며, 잡음 때문에 물리 원뿔 밖으로 나간 자료를 몰래 안쪽으로 옮겨 원래 관측처럼 제시하지 않는다.

완전한 선형 스캔도 원편광 정보는 놓친다

각도 $\gamma$의 선형 분석기에는

$$I(\gamma)=\tfrac12[S_0+S_1\cos(2\gamma)+S_2\sin(2\gamma)]$$

가 성립하고 $S_3$ 항이 없다. 따라서 $(1,0,0,0)$, $(1,0,0,1)$, $(1,0,0,-1)$은 모두 상수1/2의 곡선을 만든다. 그러나 편광도는 각각0·1·1이다. 선형 분석기의 각도를 더 많이 측정하거나 잡음을 줄여도 빠진 원편광 성분이 저절로 생기지 않는다.

이 스캔의 대비는 일반적으로 전체 편광도가 아니라 선형 편광의 비율을 알려 준다. 교정된 파장판과 선형 분석기 등의 원편광 민감 관측으로 나머지 성분을 측정해야 한다. 이는 관측 사상의 식별 한계이며 더 큰 신경망으로 자동 해결되는 최적화 실패가 아니다.

부분 편광의 수치 예제

수평·수직 단위 Jones 상태를3:1로 혼합한다. 정규화한 행렬은 $\operatorname{diag}(3/4,1/4)$, Stokes는 $(1,1/2,0,0)$, 편광도는1/2다. 고유값이3/4와1/4이므로 하나의 결정적 Jones 벡터의 외적이 될 수 없다. 두 실현의 물리적 혼합과 진폭 합을 구별한다.

축이 $\theta=\pi/4$인 사분파장판을 통과시키면 정한 관례에서 수평 입력은 음의 $S_3$, 수직 입력은 양의 $S_3$인 원형 상태가 된다. 따라서 출력은 $(1,0,0,-1/2)$이고 강도와 편광도는 유지된다. 선형 방향의 선호가 원편광 방향의 선호로 옮겨간 것이다.

네 선형 분석기 값은 모두1/2지만 양·음 $S_3$ 분석기 값은1/4와3/4다. 선형 측정만 하면 남아 있는 편광을 놓친다. 이후 수평으로 거르면 강도1/2의 완전 수평 상태가 된다. 조건부 편광도가 높아졌어도 전력이 늘어난 것은 아니다.

유효 상태와 물리적 변환은 다른 조건이다

적합한 실수 $4\times4$ 행렬이 검사한 Stokes를 모두 유효한 출력으로 보낸다고 비음수 결정적 Jones 혼합으로 구현되는 것은 아니다. 이 모델 클래스에서는 전체 Stokes 원뿔을 보존하는 조건조차 더 약하다. 상태 제약만 검사하고 장치의 물리성을 인증하지 않는다.

반례는 결맞음 행렬의 전치다. 양의 행렬을 양의 행렬로 보내고 고정 관례에서 $M=\operatorname{diag}(1,1,1,-1)$에 해당한다. 그러나 비정규화 Choi 행렬

$$\mathcal C_\Phi=\sum_{i,j}|i\rangle\langle j|\otimes\Phi(|i\rangle\langle j|)$$

의 고유값은 $(-1,1,1,1)$이다. $C\mapsto\sum_rK_rCK_r^\dagger$ 꼴이면 Choi가 양의 준정부호이므로 이 전치는 무작위 Jones 형태가 될 수 없다. 이는 고정한 복소 선형 Jones 집단 모델에 대한 설명이지 모든 능동 소자나 위상 켤레 소자 클래스의 가능성을 부정하는 말은 아니다.

강도 보존·수동성·물리적 구현 가능성도 서로 다른 제약이다. 손실 Jones 변환은 대각합을 보존할 필요가 없다. 하나의 출력 상태가 유효하다는 근거만으로 학습한 광학 장치 전체가 타당하다고 결론 내리지 않는다.

연구 설계와 교정 및 AI 평가

관측에 어떤 정보가 있는지 정한 뒤 분석기 읽음값·결맞음 행렬·소자 사상 중 무엇을 예측할지 선택한다. 양의 원뿔 제약은 유효하지 않은 상태를 줄이지만, 빠진 $S_3$ 관측이나 알 수 없는 검출 이득·틀린 소자 모델을 해결하지는 않는다. 한 입력에서 유효한 행렬보다 실제 물리적 채널 구성이 있는 모델이 더 강한 근거를 갖는다.

무편광 근처에서는 여러 성분의 잡음을 노름으로 합치면서 추정 편광도가 양수로 치우칠 수 있다. 작은 양의 추정치를 모두 고유 편광으로 해석하기보다 불확실성과 교정 가정을 기록한다. 현재의 정확한 집단 모델은 그런 실험 오차 분포를 시뮬레이션하지 않는다.

지연판의 편광도 보존, 손실을 동반한 필터 선택, 영강도, 선형 관측의 모호성, 비가환 순서를 각각 대조군으로 둔다. 독립 실수장 적분·고유값 편광도·여섯 입력 상태의 Mueller 복원은 서로 다른 계산 부분을 검사한다. 수학 검사가 통과해도 독자 화면의 실제 시각화와 조작 검증은 별도로 필요하다.

해설 문제

문제 1 — 혼합인가 결맞은 중첩인가?

$H,V$의 같은 비율 비결맞은 혼합과 결정적 상태 $(H+V)/\sqrt2$를 비교하라.

해설: 전체 강도는 둘 다1이다. 혼합은 $C=I_2/2$이고 중첩은 비대각 성분1/2을 가진 순수 대각 편광이다. 대각 분석기는 혼합에서1/2, 결맞은 중첩에서1을 투과시킨다. 실현 관측의 평균을 진폭의 합으로 바꾸면 이 차이를 없애 버린다.

문제 2 — 제안된 Stokes의 물리성

$(1,0.8,0.8,0)$은 물리적 상태인가?

해설: 아니다. 편광 벡터 크기가 $\sqrt{1.28}>1$이므로 고유값 $(1-\sqrt{1.28})/2$가 음수다. 적절한 분석기에서는 음의 예측 강도가 나온다. 거절하거나 잡음이 있는 역추론 문제를 명시적으로 풀어야 하며 그대로 물리 상태라고 받아들이지 않는다.

문제 3 — 준비 비율에서 편광도 구하기

$H,V$를3:1 가중치로 준비했을 때 편광도와 편광 강도 비율을 구하라.

해설: 정규화한 $C=\operatorname{diag}(3/4,1/4)$에서 $R=1/2$다. 전체 강도1이므로 편광도와 편광 강도 비율은 모두1/2다. 양의 강도의 무편광 성분이 남아 있으며 작은 진폭의 결정적 타원 하나로 바꿀 수 없다.

문제 4 — 빠진 관측 알아보기

완벽한 선형 분석기 스캔이 상수1/2이면 무편광이라고 확정할 수 있는가?

해설: 없다. 순수 원편광도 같은 곡선을 만든다. 원편광 민감 분석기 쌍은 $S_3=0$과 $S_3=\pm1$을 구분한다. $S_3$에 무관한 관측을 아무리 많이 정확하게 반복해도 그 성분을 유일하게 정할 수 없다.

문제 5 — 전력 생성 없이 편광도 높이기

수평 분석기는 $C=I_2/2$를 어떻게 바꾸는가?

해설: 출력은 $\operatorname{diag}(1/2,0)$이며 강도1/2·편광도1이다. 직교 성분을 제거했기 때문에 투과 부분집단이 순수해졌다. 정규화된 방향만 보고하면 전체 투과량이 감소한 사실이 가려지므로 두 값을 함께 제시한다.

문제 6 — 출력이 유효하면 장치도 유효한가?

학습한 사상이 $S_3$만 뒤집고 Stokes 원뿔을 보존한다. 이것으로 무작위 Jones 구현을 인증할 수 있는가?

해설: 없다. 고정한 기저에서는 행렬 전치에 해당하며 Choi 고유값에−1이 들어 있다. 비음수 Jones 혼합에는 양의 준정부호 Choi가 필요하다. 개별 출력의 유효성보다 실제 장치 구성의 제약이 더 강하다는 반례다.

자료와 재현 가능한 근거

NIST StokesVector와 관례 문서는 사용한 편광 부호를 명시한다. NIST MuellerMatrix는 원뿔 보존과 더 강한 Jones 혼합 물리성 조건을 구별한다. TU Delft Interactive Optics는 광학 소자의 맥락을 제공한다. 위 집단 예제와 대수는 정한 모델에서 직접 유도했다.

stokesPolarizationExperiment는 비음수 유한 준비·구성한 Pauli 채널·결정적 소자를 사용한다. 원래 Stokes 입력은 원뿔 밖이면 거절한다. 정규화 편광도의 보정은 결맞음 양성을 확인한 뒤 명시한 부동소수점 오차 범위에만 적용하며 원래 행렬은 보존한다. 독립 생성기는 실제 분석기 장을 적분하고 여섯 물리 입력으로 Mueller 열을 복원하며 고유값으로 편광도를 구한다. 이 근거를 미실행 홈페이지 통합이나 설치하지 않은 NIST 라이브러리 검사로 확대하지 않는다.