Nature 658, 658-665

X-ray liquidography decodes complex motions in azobenzene isomerization

용매 신호 90%를 걷어 내자 보인 아조벤젠의 첫 동작

빛을 받은 아조벤젠이 트랜스에서 시스로 바뀌는 과정을 펨토초 X선 산란으로 찍은 Nature 논문을 읽는다. 신호의 90% 넘게 차지하는 용매 응답을 부분공간 투영으로 지우고, 운동 조합을 BIC로 골라 두 중간체의 구조를 정했다.

키워드

물리화학

한 줄 요약과 선정 근거

빛을 받은 아조벤젠이 트랜스에서 시스로 바뀔 때 처음 0.4피코초 동안 일어나는 일은 N=N 결합을 축으로 한 회전이 아니었다. 두 페닐 고리가 C-N 결합을 축으로 90도 가까이 비틀리는 동작이 먼저 왔다. 연구진은 이 장면을 펨토초 X선 산란으로 찍었다. 측정 신호의 90% 넘게 차지하는 용매 응답은, 용매만 따로 잰 실험에서 그 응답이 사는 부분공간을 구한 뒤 투영으로 걷어 냈다.

Nature는 해설 기사(News & Views) “Beyond the solvent fog: watching molecules dance in the crowd”(Kathryn Froese, Germán Sciaini, 2026년 10월 12일)를 함께 실었고, KAIST와 국내외 언론도 보도했다. KAIST 화학과 이효철 교수 연구팀과 기초과학연구원(IBS) 첨단 반응동역학 연구단이 포항가속기연구소의 X선 자유전자 레이저(PAL-XFEL)에서 얻은 자료다. 고른 까닭은 세 가지다.

1. 교과서에 나오는 반응인데도 메커니즘을 두고 수십 년째 결론이 나지 않던 문제에 구조로 답했다.
2. 무거운 원자가 없는 유기 분자에도 이 측정법이 통한다는 것을 보여, 남이 쓸 수 있는 도구의 범위를 넓혔다.
3. 거대한 배경 속 작은 신호를 부분공간 투영으로 꺼내는 법, 서로 엮인 자료에서 유효 표본 수를 따로 재서 BIC를 쓰는 법처럼 ML 연구자가 가져갈 기법이 많다. 심사 기록과 코드도 공개돼 있다.

배경지식 세 가지

1. 아조벤젠의 광이성질화(photoisomerization): 페닐 고리 두 개를 N=N 다리로 이은 분자로, 평평한 트랜스 형태가 빛을 받으면 구부러진 시스 형태가 된다. 빛으로 켜고 끄는 분자 스위치로 널리 쓰인다. 이 논문은 400nm 빛으로 nπ* 들뜸(S0 → S1)을 일으켰다. 바뀌는 방식을 두고 N=N 결합을 축으로 고리가 도는 회전(rotation), 질소의 C-N-N 각이 평면 안에서 벌어지는 반전(inversion), 훌라 트위스트(hula twist) 같은 설명이 맞서 왔다(그림 1a).
2. X선 리퀴도그래피(time-resolved X-ray liquidography, TRXL): 펌프 레이저로 반응을 시작하고 정해진 시간 t 뒤에 X선 펄스로 용액의 산란 무늬를 찍는다. 레이저를 쏜 무늬에서 쏘지 않은 무늬를 빼면 차이 산란 곡선 ΔS(q, t)가 나온다. q는 산란 벡터의 크기다. 용액 속 분자는 아무 방향으로나 놓여 있어서 등방성 신호에는 원자 사이 거리의 분포만 담긴다. 구조 하나가 내는 산란은 아래의 데비 식(Debye equation)으로 계산하고, 식의 \(f_i\)는 원자의 산란 인자, \(r_{ij}\)는 두 원자 사이 거리다(보충 방법 6).

\[S(q) = \sum_i \sum_j f_i(q)\, f_j(q)\, \frac{\sin(q r_{ij})}{q r_{ij}}\]

3. 용질 항, 케이지 항, 용매 항: 차이 신호는 세 가지의 합이다. 분자 자체의 구조 변화(용질 항), 분자와 둘러싼 용매 사이 거리의 변화(케이지 항), 그리고 들뜬 분자가 내놓은 열로 용매 전체의 온도와 밀도가 바뀌어 생기는 신호(용매 항)다. 요오드나 금속처럼 전자가 많은 원자가 있으면 용질 항이 크지만, 탄소, 질소, 수소뿐인 아조벤젠에서는 용매 항이 전체의 90%를 넘는다(그림 2c). 그래서 지금까지 TRXL은 거의 무거운 원자가 든 분자에만 쓰였다(그림 1c).

무엇을 발견했나

세 단계의 순차 반응. 용매를 걷어 낸 자료를 특이값 분해(SVD)하자 의미 있는 성분이 셋 나왔다. 세 성분의 시간 변화를 공통 시간 상수를 쓰는 지수 함수의 합으로 함께 맞추니 0.4 ± 0.2ps, 2.7 ± 1.0ps, 14 ± 4ps가 나왔다(본문, 그림 3c). 연구진은 이를 첫 중간체 I1이 생기고(0.4ps), I1이 I2로 바뀌고(2.7ps), I2가 시스가 되는(14ps) 순차 모형으로 읽었다. 구조 맞춤에서 나온 분기 비율로는 I1의 100 ± 19%가 I2로, I2의 91 ± 7%가 시스로 갔다(본문, 보충 방법 12.3).

첫 동작은 고리의 비틀림이다. I1에서는 두 페닐 고리가 C-N 결합을 축으로 87도씩 서로 반대 방향으로 비틀렸다(∠CCNN = ±87.0°). 그런데 N=N 축의 이면각 ∠CNNC(180°)와 질소의 각 ∠CNN(113.4°)은 트랜스 그대로였다(본문, 그림 4a). 회전도 반전도 아닌 C-N 비틀림(C-N torsion, CNT)으로 시작한 셈이다. 이때 가운데 질소 다리가 자전거 페달처럼 맞물려 움직여서 고리가 용매를 크게 쓸고 지나가지 않는다. I2에서는 비틀림(−71.9°)이 남은 채 ∠CNNC가 115.5°까지 돌았다. 이 단계의 운동은 회전(ROT) 67%, 비틀림 33%로 나뉘었고, 마지막 I2 → 시스 단계에서야 반전(INV)이 39%로 크게 끼어들었다(보충 표 5). 비틀림이 먼저 오고 반전은 맨 나중이다.

(가) 구조 맞춤에 쓴 네 가지 기본 운동

(나) 단계별로 어떤 운동이 일했나

(가) 논문은 중간체 구조를 네 가지 기본 운동의 조합으로만 만들어 자료에 맞췄다. ROT는 N=N 결합을 축으로 한쪽 고리를 돌려 ∠CNNC를 바꾸고, INV는 질소에서 ∠CNN을 벌리며, NND는 N-N 결합 길이를, CNT는 고리를 C-N 결합 축으로 비틀어 ∠CCNN을 바꾼다. 평면으로 단순화한 개념도다. (나) 앞 구조에서 다음 구조로 넘어갈 때 각 운동을 하나씩 빼고 나머지를 다시 맞췄을 때 구조 오차가 얼마나 커지는지로 잰 기여(MoCA, 보충 표 5)다. 위 막대는 실험으로 정한 구조, 아래 막대는 같은 방법으로 분해한 양자화학 계산(CASPT2) 구조다. 같은 운동이 나오지만 순서가 다르다.

계산과는 순서가 달랐다. 연구진이 기준점으로 삼은 다중참조 양자화학 계산(XMS-CASPT2)의 구조도 같은 네 가지 운동으로 분해되지만, 반전이 일찍 들어오고 비틀림은 마지막에 커지는 정반대 순서였다(본문, 보충 표 5). 계산 구조를 실험 자료에 맞춰 다듬자 χν²가 I1에서 47%, I2에서 58% 줄었고 바닥 상태인 시스는 19%만 줄었다(Methods). 들뜬 상태에서 계산이 가장 많이 어긋났다는 뜻이다. 많은 이론 연구가 지지해 온 반전 보조 회전 메커니즘은 실험에서 본 ∠CCNN의 큰 변화를 설명하지 못한다(본문).

부피를 아끼는 경로. 용매의 점도를 바꿔도 반응 속도가 거의 변하지 않는다는 오래된 관찰은 고리를 크게 휘두르는 회전과 맞지 않았다. 모든 원자가 함께 조금씩 움직이는 기준틀에서 단계마다 용매를 쓸고 지나가는 부피는 20, 37, 114ų로, 질소를 고정하고 고리만 움직였을 때(104, 65, 115ų)보다 첫 단계에서 81% 작다(본문, 보충 표 7).

논문 그림 4. (a) 반응 경로를 따라 여섯 구조 변수를 실험 구조(빨강, TRXL)와 계산 구조(파랑, CASPT2)로 비교했다. 가장 큰 차이는 오른쪽 위의 ∠CCNN으로, 실험에서는 I1에서 이미 87도로 꺾이지만 계산에서는 I2까지 거의 0도다. (b) 네 구조를 위, 앞, 옆에서 본 모습으로, 옅은 그림은 앞 구조에서 넘어오는 중간 모습이다. 출처: Kim 외, Nature (2026), CC BY 4.0.

어떻게 했나

측정. 메탄올에 녹인 50mM 트랜스 아조벤젠을 400nm 레이저로 들뜨게 하고 15.01keV X선 펄스로 찍었다. 시간 지연은 −0.5ps부터 100ps까지 140개, 지연마다 레이저를 쏜 영상이 약 4,500장이고, 시간 분해능은 169 ± 5fs, q 범위는 1~7Å⁻¹다(Methods). 용매 응답의 기준을 얻으려고 순수 메탄올도 따로 쟀다.

용매 걷어 내기. 순수 메탄올 측정을 SVD해서 신호로 볼 만한 성분 7개를 골랐다. 시간 방향 자기상관이 0.5 이상이고 특이값이 잡음 바닥보다 뚜렷이 큰 성분이다(Methods). 그다음 아조벤젠 자료를 이 7차원 부분공간에 수직인 쪽으로 투영했다. 같은 연구팀이 2023년에 내놓은 ’수직 성분을 뽑는 투영(projection to extract the perpendicular component, PEPC)’이다. 순수 용매 자료에 같은 투영을 걸면 아무 무늬도 남지 않는다(본문, 확장 데이터 그림 2). 부작용도 적혀 있다. 성분별 차이 산란 곡선(species-associated difference scattering curve, SADS)의 낮은 q에서 보이는 출렁임은 투영이 만든 것이라 구조 정보로 읽으면 안 된다. 시간 0을 알려 주는 용매의 광 커 효과(optical Kerr effect) 신호도 함께 지워지기 때문에, 시간 0은 투영 전 자료에서 정했다(Methods).

논문 그림 2. (a) 원자료(맨 왼쪽)를 용질 항, 케이지 항, 용매 항으로 나눈 것이다. 가로축은 시간, 세로축은 q, 색은 q를 곱한 차이 산란 신호다. (b) 용매 부분공간에 수직인 쪽으로 투영한 뒤(PEPC)의 모습으로, 열 배로 키워 그렸다. (c) 각 항이 신호에서 차지하는 비율이다. 출처: Kim 외, Nature (2026), CC BY 4.0.

(가) 투영이 남기는 것과 지우는 것

(나) 투영한 공간에서 맞추기

(가) 차이 산란 곡선 하나를 q 점 511개짜리 벡터로 보면, 용매가 데워지며 내는 응답은 용매만 따로 잰 측정에서 SVD로 뽑은 성분 몇 개가 펼치는 부분공간(파란 평면) 안에 거의 다 들어간다. PEPC는 측정 벡터에서 이 평면에 수직인 성분(보라)만 남긴다. 용매 응답은 통째로 사라지지만, 용질 신호 가운데 평면 쪽으로 누운 부분(빨간 점선)도 함께 지워진다. (나) 그래서 논문은 투영 전의 깨끗한 용질 신호를 되살리려 하지 않고, 후보 구조로 계산한 곡선에도 같은 투영을 걸어 투영된 공간에서 χ²를 잰다(보충 방법 5). 실제 공간은 수백 차원이고 용매 성분은 등방성 신호에서 7개다. 3차원으로 줄여 그린 개념도다.

용매 항이 얼마나 압도적인지, 투영이 무엇을 함께 지우는지는 논문이 공개한 그림 2의 원천 자료로 계산해 볼 수 있다.

(가) 원자료를 채운 것은 용매 응답

(나) 투영이 용질 신호에서 지운 것

논문 그림 2의 원천 자료로 계산했다. (가) 원자료(등방성 차이 산란)를 이루는 세 항이 차지하는 비율이다. 시간마다 q를 곱한 곡선의 절댓값을 q = 1~7 Å⁻¹에서 적분한 넓이로 재고(논문 그림 2c와 같은 척도), 로그 간격의 시간 구간마다 더했다. 0.1~0.2 ps에서도 용매 항이 79%이고, 열이 용매로 퍼지는 수십 ps 뒤에는 92%가 된다. (나) 모형이 예측한 용질 항(q를 곱한 차이 산란, 세로 눈금은 구간마다 다름)을 두 시간 구간에서 평균해 투영 전과 후로 겹쳤다. 첫 중간체가 대부분인 이른 구간에서는 4~6 Å⁻¹의 무늬가 거의 그대로 남지만, 시스가 대부분인 늦은 구간에서는 1~2 Å⁻¹의 큰 성분이 대부분 지워진다. 0 ps 이후 전체로는 용질 항 넓이의 75%가 남는다. 투영 뒤 곡선이 거친 까닭은 투영에 쓴 용매 성분이 측정에서 왔기 때문이다(보충 방법 5).

오차 막대. 이웃한 시간 점을 평균하는 전처리가 잡음을 엮어 놓아 해석적 오차가 약 √2배 작게 나오므로, 시간 상수의 오차는 잡음을 새로 입힌 자료 500벌로 붓스트랩해 구했다(보충 방법 1).

구조 맞춤과 모형 선택. 중간체 구조는 네 가지 기본 운동(ROT, INV, NND, CNT)의 조합으로만 만들고, 부분집합 15개를 모두 맞춰 BIC로 골랐다(Methods). 계산한 곡선(데비 식과 분자동역학으로 구한 케이지 항)에도 같은 투영을 걸어 비교했고, 케이지 항은 구조가 바뀔 때마다 다시 계산해 수렴할 때까지 되풀이했다(보충 방법 5, 8, 9).

분야에서의 의미

• 탄소, 질소, 수소만으로 된 분자의 용액 속 구조 변화를 X선 산란으로 풀어낸 첫 사례다. 그림 1c는 지금까지 TRXL로 다룬 분자와 이제 다룰 수 있게 된 분자(스피로피란, 디아릴에텐, 벤조페논 등)를 나란히 보여 준다. 유기 광화학의 넓은 영역이 구조 측정의 대상이 된다.
• 메커니즘 논쟁을 ’각자 무엇을 봤나’로 정리했다. 회전을 지지한 분광 연구는 ∠CNNC 변화를, 반전을 지지한 연구는 마지막 단계의 ∠CNN 변화를 본 것일 수 있다(보충 노트 8). 연구진은 구조 변수 몇 개만 보고 메커니즘을 분류해 온 습관을 혼란의 원인으로 꼽는다.

ML 연구자를 위한 인사이트

1. 배경을 맞추지 말고 배경이 사는 부분공간을 지우고, 모형에도 같은 투영을 건다 (방법론 교훈, 직접 적용)

PEPC는 직교 투영 \(P_\perp = I - UU^\top\) 하나이고, \(U\)는 목표 신호가 없는 대조 측정에서 SVD로 뽑은 기저다. ML에서도 영공간 투영(INLP, Ravfogel 외, ACL 2020)과 LEACE(Belrose 외, NeurIPS 2023)가 표현에서 개념의 부분공간을 지울 때, 대조 PCA(Abid 외, Nature Communications 2018)가 배경 자료와 견줘 목표 자료에만 많은 방향을 찾을 때 같은 연산을 쓴다. 배울 점은 세 가지다.

• 지울 부분공간은 목표 신호가 없는 대조 측정에서 따로 얻는다. 목표 자료에서 배경을 추정하면 신호가 기저에 섞인다.
• 투영은 배경과 겹치는 신호도 지운다. 위 계산에서 용질 항 넓이의 약 4분의 1이 사라졌고, 시스 시기의 낮은 q 성분은 대부분 깎였다.
• 그래서 깨끗한 신호를 되살리려 하지 않고 예측에도 같은 투영을 걸어 같은 공간에서 비교한다.

블로그의 해석으로는, 선형 최소제곱에서 맞춤에 쓰는 내적과 같은 내적으로 배경 성분을 직교 투영해 지운다면, 자료와 모형을 함께 투영하는 것은 배경 계수를 아무 제약 없이 같이 맞추는 것과 답이 같다(프리슈 워 러벨 정리). 잃은 정보를 되찾으려면 배경 쪽에 제약이 있어야 한다. 개념을 지운 표현으로 다른 과제를 평가할 때도 같은 손실을 따져야 한다.

2. 후보 15개를 모두 맞추고, 유효 표본 수를 따로 재서 BIC를 쓴다 (방법론 교훈, 직접 적용)

연구진은 모든 운동 조합을 맞춘 뒤 \(\mathrm{BIC} = \chi^2 + (k+4)\ln N_{\mathrm{eff}}\)로 골랐다(\(k\)는 운동 매개변수 수, 4는 진동 폭, 크기, 기준선 두 개, 보충 방법 13). 눈여겨볼 것은 \(N_{\mathrm{eff}}\)다. q 점은 511개지만 이웃한 점은 거의 같은 정보를 담는다. 연구진은 곡선을 재현하는 데 필요한 체비쇼프 다항식의 차수를 다시 BIC로 정해 \(N_{\mathrm{eff}} = 110\)을 얻었고, 섀넌 표본화 정리로 어림한 60~120과 맞는다는 것을 보였다. χν²가 1보다 작게 나와도 과적합이 아니라는 설명도 여기서 나온다(Methods). 선택된 모형(I1은 CNT 하나, I2는 ROT + CNT)은 \(N_{\mathrm{eff}}\)를 23에서 511까지 바꿔도 그대로였다(보충 표 6). 시계열이나 증강한 표본처럼 서로 엮인 자료로 정보 기준을 계산할 때 따라 할 만한 점검이다.

3. 기여도는 하나씩 빼고 다시 맞춰서 잰다 (방법론 교훈, 영감 수준)

단계마다 운동별 기여를 잰 MoCA는 한 운동을 빼고 나머지를 다시 최적화했을 때 구조 오차(RMSD)가 얼마나 늘어나는지를 백분율로 나눈 값이다(보충 방법 15). 특징 하나를 빼고 다시 학습하는 drop-column 중요도와 같은 발상이고, 연구진도 이 증가분을 더하면 전체 개선과 같아지지 않는다고 밝힌다. 블로그의 해석으로는, 운동이 넷뿐이고 15개 조합을 이미 다 맞춰 두었으니 섀플리 값(Shapley value)을 정확히 계산할 수 있다. 섀플리 값은 더하면 전체 개선과 같아지고 빼는 순서에 기대지 않는다. 그림의 순서(CNT 먼저, INV 나중)가 섀플리 값으로도 그대로 나오는지 확인해 볼 만하다.

4. 1차원 거리 정보로 3차원 구조를 고르는 역문제 (새 문제와 데이터, 영감 수준)

심사 과정에서 자료에서 구조 변화를 바로 읽어 낼 수 없느냐는 질문을 받자, 연구진은 용액 산란이 원자 쌍 거리의 1차원 분포만 담으므로 자료가 아무리 좋아도 3차원 구조를 되찾으려면 구조 모형이 필요하다고 답했다(심사 기록). 이 논문에서 그 사전 지식은 사람이 고른 네 가지 운동이다. 학습된 분자 구조 사전 분포에 미분 가능한 데비 정방향 모형을 붙이면, 가장 잘 맞는 구조 한 개가 아니라 자료와 맞는 구조들의 분포를 얻을 수 있다. 원자 24개라 정방향 모형은 가볍고, SADS를 포함한 원천 자료와 구조 좌표, 코드가 공개돼 있다. 실험이 계산과 다른 순서를 보였으니, 양자화학 자료로 학습한 들뜬 상태 ML 퍼텐셜(아조벤젠 유도체를 다룬 Axelrod 외, Nature Communications 2022 등)을 채점하는 기준으로도 쓸 수 있다. 둘 다 논문이 시험한 것은 아니어서 영감 수준이다.

한계와 남은 질문

• 결론은 네 가지 운동 안에서 가장 단순한 설명이다. 고리의 비평면 각, C-C-N 각, C-N 거리 같은 다른 좌표가 끼어들 가능성은 연구진도 배제하지 않았다(본문). 최선 모형과 다음 모형의 BIC 차이는 2.2~5.9로(보충 표 6), 흔히 쓰는 기준으로 지지가 분명한 수준이지 매우 강한 수준은 아니다.
• 성분 수 셋은 경계선에서 고른 값이다. 세 번째와 네 번째 성분의 특이값과 자기상관(0.71, 0.69)이 거의 같아서, 최종본은 세 번째 성분의 시간 변화에만 오르내림이 뚜렷하다는 점을 근거로 삼았다(보충 방법 1, 심사 기록).
• 등방성 자료는 분자의 절대 방향을 알려 주지 않는다. 부피를 아끼는 기준틀을 택한 근거는 점도와 무관하다는 기존 결과다(본문).
• 판단이 들어간 곳이 있다. 시스의 진동 폭(σ²)을 자유롭게 두면 0.023Ų로 뜨거운 중간체보다 크게 나와, 연구진은 물리적으로 맞지 않는다며 0으로 고정했다(보충 방법 6).
• 들뜬 직후 곧장 트랜스로 돌아가는 분자는 구조 신호를 남기지 않는다. 분광 연구의 시스 생성 수율 0.32를 받아들이면 들뜬 분자의 약 65%가 이 길로 간다(보충 노트 10). 용매도 들뜸 파장도 한 가지씩만 다뤘다.
• 심사위원 연구실의 학생 둘이 공개 코드를 직접 돌려, 시스 구조의 오차 막대가 서로 다른 두 극소점을 섞어 계산됐다는 점과 \(N_{\mathrm{eff}}\)가 거의 용매 성분으로 정해진다는 점을 짚었고 최종본에 반영됐다(심사 기록).

열린 질문

1. 배경 부분공간을 투영으로 지울 때 함께 사라지는 신호를 배경 쪽에 물리적 제약을 걸어 되찾을 수 있을까. 그 이득은 신호와 배경이 얼마나 겹칠 때 커질까.
2. 사람이 고른 네 가지 운동 대신 학습된 구조 사전 분포를 쓰면, 같은 자료가 C-N 비틀림으로 시작하는 경로만 허락할까, 아니면 다른 경로도 자료와 맞을까.
3. 무거운 원자가 없는 다른 광스위치(스피로피란, 디아릴에텐)에서도 계산과 실험이 운동의 순서에서 갈릴까.

이해 확인

1. PEPC 투영을 거친 뒤 남는 신호를 바르게 설명한 것은?

1. 용매 항만 사라지고 용질 항은 손실 없이 그대로 남는다
2. 용매 항이 사라지고, 용질 항 가운데 용매 부분공간 쪽으로 놓인 부분도 함께 사라진다
3. 용매 항의 크기만 줄고 모양은 남는다
4. 케이지 항까지 모두 사라지고 용질 항만 남는다

PEPC는 순수 용매 측정에서 얻은 부분공간에 수직인 성분만 남기는 직교 투영이라, 용질 신호도 그 부분공간 쪽 성분은 잃는다(원천 자료로 계산하면 용질 항 넓이의 약 4분의 1). 그래서 계산한 곡선에도 같은 투영을 걸어 비교했다(보충 방법 5). 케이지 항은 투영 뒤에도 남는다(그림 2b).

2. 트랜스에서 첫 중간체 I1로 갈 때 일어난 구조 변화는?

1. N=N 결합을 축으로 고리가 돌아 ∠CNNC가 180°에서 크게 줄었다
2. 질소의 C-N-N 각이 평면 안에서 크게 벌어졌다
3. 두 고리가 C-N 결합을 축으로 87도씩 비틀렸고, ∠CNNC와 ∠CNN은 그대로였다
4. N=N 결합 길이가 크게 늘었다

BIC는 트랜스 → I1을 C-N 비틀림(CNT) 하나로 설명하는 모형을 골랐다. I1의 ∠CCNN은 ±87.0°이고 ∠CNNC(180°)와 ∠CNN(113.4°)은 트랜스와 같다(그림 4a). 회전은 I1 → I2에서 주도하고(67%), 반전은 마지막 단계에서야 커진다(39%).

3. 운동 조합을 고르는 BIC의 벌점에 q 점의 수 511 대신 N_eff = 110을 쓴 까닭은?

1. 계산 시간을 줄이려고 q 점을 솎아 냈기 때문이다
2. 이웃한 q 점이 거의 같은 정보를 담아서, 독립적인 정보의 수가 측정 점의 수보다 훨씬 적기 때문이다
3. 높은 q의 자료를 믿을 수 없어 버렸기 때문이다
4. 벌점을 키워 더 복잡한 모형이 뽑히게 하려고

산란 곡선은 측정 간격보다 훨씬 느리게 변해서 이웃한 q 점이 엮여 있다. 연구진은 곡선을 재현하는 데 필요한 다항식 차수로 독립 정보의 수를 110으로 잡았고, 섀넌 표본화 정리의 어림(60~120)과도 맞았다(보충 방법 13). N_eff가 작을수록 벌점이 작아지므로 4번은 방향이 거꾸로다.

논문 정보

• 논문: Jungmin Kim, Hosung Ki, Seonggon Lee, Alekos Segalina, Yunbeom Lee, Hyotcherl Ihee, “X-ray liquidography decodes complex motions in azobenzene isomerization”, Nature 658, 658-665. 2026년 9월 30일 온라인 게재. doi:10.1038/s41586-026-11068-4
• 해설 기사: Kathryn Froese, Germán Sciaini, “Beyond the solvent fog: watching molecules dance in the crowd”, Nature 658, 613, 2026년 10월 12일. doi:10.1038/d41586-026-02970-y
• PEPC 방법: Hosung Ki 외, “Projection to extract the perpendicular component (PEPC) method for extracting kinetics from time-resolved data”, Structural Dynamics 10, 034103 (2023). doi:10.1063/4.0000189
• 코드: Code Ocean capsule (전처리, PEPC, SVD와 동역학 분석, 분자동역학 케이지 항, SSPM, BMA, MoCA)
• 자료: 구조 좌표, 원천 자료, 심사 기록은 논문 페이지의 보충 자료