단일 입자 회절 영상 복원을 위한 EMC 재구성 알고리즘

초록

EMC 알고리즘은 입자 회전이 알려지지 않은 수많은 저광자 2D 회절 패턴으로부터 3차원 회절 강도를 효율적으로 추정한다. 확률적 최대우도(M) 단계와 3D 모델을 중복된 투시 표현으로 전·후 변환하는 확장(E)·압축(C) 단계가 순환하며, 몇 차례 반복 후 차이 지도(diff‑map) 알고리즘에 넘겨 입자 구조를 복원한다. 시뮬레이션을 통해 입자 복잡도, 측정 수, 프레임당 광자 수가 재구성 품질에 미치는 영향을 정량화하고, 향후 XFEL 실험에 필요한 데이터 처리 규모를 예측한다.

상세 분석

EMC(Expansion‑Maximization‑Compression) 알고리즘은 단일 입자 회절 영상(single‑particle diffraction imaging, SPI)에서 가장 큰 난제인 ‘방향성 불확실성’을 통계적 기대값(maximization)과 기하학적 변환(expansion/compression)으로 동시에 해결한다. 먼저 3차원 회절 강도 I(q) 를 초기화하고, 이를 2차원 투시 이미지들의 집합에 매핑하기 위해 ‘확장(E)’ 단계에서 각 가능한 회전 R에 대해 I(q) 를 회전된 좌표계에 투사한다. 이렇게 얻어진 중복된 2D 강도 집합은 실제 측정된 포톤 수 n_i와 포아송 확률 모델을 이용해 ‘최대우도(M)’ 단계에서 각 회전 가중치 w_i(R)를 업데이트한다. 여기서 w_i(R)∝exp