분자는 빛을 받은 뒤 전자만 바뀌는 정지된 물체가 아니다. 원자핵이 움직이고 결합이 늘어나거나 비틀리면서 전자상태의 에너지와 성격도 함께 달라진다. 이 진동과 전자의 결합 운동, 즉 진동-전자(vibronic) 동역학은 내부전환, 계간전이, 에너지·전하이동처럼 시간에 따라 진행되는 광물리 현상의 중심에 있다.
2026년 8월 19일 공개된 PennyLane 데모 A quantum algorithm for vibronic dynamics는 이 문제를 양자회로로 옮기는 과정을 코드로 보여준다. 기반 방법은 Motlagh 등이 발표한 동료평가 논문 Quantum Algorithm for Vibronic Dynamics: Case Study on Singlet Fission Solar Cell Design, Quantum Science and Technology 10, 045048 (2025)이다.
핵심 판정부터 말하면, 이 데모는 실제 양자 프로세서에서 실제 분자를 계산한 결과가 아니다. 두 전자상태와 한 진동모드로 만든 작은 모형을 CPU 기반 lightning.qubit 상태벡터 시뮬레이터에서 실행했다. 가치는 다른 곳에 있다. 진동-전자 해밀토니안을 실공간 격자, Quantum Fourier Transform(QFT), Quantum Read-Only Memory(QROM), 가역 산술, 위상 기울기 상태와 2차 Trotter 전개로 연결해, 미래 내결함성 양자컴퓨터용 알고리즘의 내부 구조를 학습 가능한 형태로 공개했다.
그림 1. 전자상태 변화, 핵 파동묶음의 운동, 실공간 격자 인코딩을 연결한 개념 일러스트. 특정 분자의 구조식, 실제 양자회로 또는 계산 결과를 재현한 그림이 아니다.
한눈에 보는 검증 결과
| 항목 | 확인된 내용 |
|---|---|
| 공개물 | Emily Nobes, PennyLane 교육용 데모, 2026년 8월 19일 |
| 기반 연구 | Motlagh et al., Quantum Science and Technology 10, 045048 (2025) |
| 계산 방식 | 디지털 gate 기반 실시간 해밀토니안 시뮬레이션; VQE·QAOA·양자 어닐링이 아님 |
| 실행 장치 | CPU C++ 상태벡터 시뮬레이터 lightning.qubit |
| 데모 모형 | 전자상태 2개, 진동모드 1개, 모드당 격자점 4개, 선형 결합 |
| 전체 크기 | 23 simulated wires: system 3개 + 보조·정밀도 wire 20개 |
| 출력 | 두 diabatic 전자상태의 시간별 정확 확률(exact probability) |
| 직접 보여주지 않은 것 | 실제 분자 정확도, QPU 성능, noise·오류정정, 양자우위, OLED 물성 예측 |
| 연구적 위치 | 초기 내결함성 양자컴퓨팅을 겨냥한 알고리즘·자원 타당성 연구 |
이 글은 다음 근거 층을 구분한다.
- [원 데모] 공개 코드가 실제로 실행한 설정과 출력
- [기반 논문] 동료평가 알고리즘과 최신 arXiv v3가 제시한 자원 추정
- [계산 경계] 코드가 검증하지 않은 정확도·하드웨어·물리 조건
- [리뷰 해석] TADF·PhOLED 및 분자 역설계 연구로 연결할 때 필요한 후속 단계
1. 쉬운 설명: 여러 장의 지형도 위를 움직이는 파동묶음
Born–Oppenheimer(BO) 근사는 전자가 핵보다 훨씬 빠르게 움직인다는 차이를 이용한다. 각 핵 배치에서 전자상태를 먼저 계산하고, 핵은 그 상태가 만드는 잠재에너지표면(potential-energy surface, PES) 위를 움직인다고 본다. 평형 구조와 바닥상태 에너지를 계산할 때 매우 성공적인 접근이다.
빛을 받은 뒤에는 여러 전자상태의 표면이 가까워지거나 교차할 수 있다. 이를 여러 장의 유연한 지형도로 생각해 보자. 한 장은 한 전자상태이고, 지도 위 위치는 분자의 진동 좌표다. 핵의 파동묶음이 지형을 따라 움직이다 두 지도가 강하게 연결되는 구역을 만나면, 일부가 다른 지도에 해당하는 전자상태로 넘어갈 수 있다. 에너지 준위표는 몇 지점의 높이를 알려주지만, 어느 경로로 얼마나 빨리 이동하는지는 알려주지 않는다.
PennyLane 알고리즘은 이 움직임을 영화 프레임처럼 잘게 나눈다. 짧은 시간 $dt$마다 다음 두 일을 번갈아 수행한다.
- 운동량에 따라 핵 파동묶음을 이동시킨다.
- 핵의 위치와 전자상태에 따라 위상과 상태 간 결합을 갱신한다.
이 짧은 업데이트를 대칭 순서로 반복하는 것이 2차 Suzuki–Trotter product formula다.
PennyLane 소개문은 고전 계산이 동역학을 사실상 다룰 수 없다는 인상을 준다. 이는 범위가 너무 넓다. Surface hopping, Ehrenfest dynamics, multiconfiguration time-dependent Hartree(MCTDH), multilayer MCTDH와 tensor-network 등 여러 고전 방법이 이미 쓰인다. 더 정확한 문제 정의는 많은 진동모드와 전자상태가 강하게 상관된 전체 양자 동역학을 체계적으로 정확하게 전파할 때 비용이 빠르게 증가한다는 것이다.
2. 계산할 물리: KDC 진동-전자 해밀토니안
데모가 사용하는 Köppel–Domcke–Cederbaum(KDC) 해밀토니안은 전자상태와 핵 진동을 한 양자계로 적는다.
$$ H=\mathbb I_{\mathrm{el}}\otimes(T_{\mathrm{nuc}}+V_0)+\mathbf W(\mathbf Q). $$- $\mathbb I_{\mathrm{el}}$은 전자상태 공간의 항등연산자다.
- $T_{\mathrm{nuc}}$은 핵 진동의 운동에너지다.
- $V_0$는 기준 조화 퍼텐셜이다.
- $\mathbf W(\mathbf Q)$는 정상모드 좌표 $\mathbf Q$에 따라 변하는 전자상태 에너지와 상태 간 결합을 담은 diabatic potential matrix다.
두 전자상태 $i,j$ 사이의 퍼텐셜 블록은 기준 구조 주변에서 다음처럼 전개할 수 있다.
$$ W'_{ij}(\mathbf Q)=\lambda^{(i,j)} +\sum_r a_r^{(i,j)}Q_r +\sum_{r,r'}b_{rr'}^{(i,j)}Q_rQ_{r'}. $$$\lambda$는 좌표와 무관한 에너지 또는 결합, $a_rQ_r$은 선형 진동-전자 결합, $b_{rr'}Q_rQ_{r'}$은 이차 및 모드 간 결합을 나타낸다. 선형항까지 남기면 linear vibronic coupling(LVC), 이차항까지 포함하면 quadratic vibronic coupling(QVC) 모델이다. 여기서 $Q_r$은 원자 하나의 위치가 아니라 분자 전체가 함께 움직이는 무차원 정상모드 좌표다.
이 해밀토니안은 압축된 물리 모형이지 전자구조를 자동으로 만들어 주는 장치가 아니다. 전자상태, 진동수, 정상모드와 결합계수는 DFT/TDDFT 또는 다중참조 전자구조 계산과 diabatization으로 먼저 얻어야 한다. 기반 논문 최신 v3도 해밀토니안 구성 자체는 양자 동역학 알고리즘의 범위 밖이라고 명시한다.
3. 연속적인 핵 운동을 큐비트 격자에 넣는 법
각 진동모드를 $K=2^k$개의 위치 격자점으로 이산화하면 한 모드를 $k$개 큐비트에 저장할 수 있다. 계산기저 $|x\rangle$가 위치를 나타내며,
$$ Q|x\rangle=\Delta(x-K/2)|x\rangle, \qquad \Delta=\sqrt{2\pi/K} $$로 좌표를 정의한다. 전자상태가 $N$개이고 진동모드가 $M$개라면 물리계를 나타내는 주 레지스터는
$$ M\log_2K+\lceil\log_2N\rceil $$개의 큐비트를 사용한다.
중요한 정정이 하나 있다. 이 데모는 진동을 연속변수 qumode로 처리하는 광자 하이브리드 구현이 아니다. 전자상태 레지스터와 모든 진동 격자 레지스터를 모두 큐비트로 표현하고 lightning.qubit에서 계산한다. 코드의 electrons 변수도 전자 수가 아니라 diabatic electronic-state register를 뜻한다.
고전적인 primitive grid는 $N K^M$개의 복소 진폭을 명시적으로 저장해야 하므로 모드 수에 따라 메모리가 지수적으로 늘 수 있다. 양자상태는 이 진폭을 더 작은 큐비트 레지스터에 부호화할 수 있다. 하지만 이 압축이 곧 저비용 계산을 뜻하지는 않는다. 상태 준비, 계수 로딩, 양자 산술, Trotter 반복, 오류정정과 측정 비용이 남는다.
그림 2. KDC 진동-전자 모형을 전자상태·진동 격자 큐비트로 옮기고 대칭 2차 Trotter 단계로 전파하는 흐름. 하단은 PennyLane 코드가 실제 실행한 toy 범위를 표시한다. 모든 진동 좌표는 qumode가 아니라 큐비트 격자로 인코딩된다.
4. 한 번의 시간 진화를 회로로 조립하는 방법
목표는 $|\psi(t)\rangle=e^{-iHt}|\psi(0)\rangle$를 계산하는 것이다. 운동에너지 $T$와 퍼텐셜 $V$는 일반적으로 서로 교환하지 않기 때문에 지수연산을 정확히 두 조각으로 분리할 수 없다. 데모는 운동에너지 절반 단계, 퍼텐셜 조각의 정방향·역방향 전개, 다시 운동에너지 절반 단계를 적용하는 2차 Trotter 공식을 사용한다.
운동에너지: QFT로 운동량 공간에 다녀온다
위치 $Q$는 격자 계산기저에서 대각형이지만 운동량 $P$는 그렇지 않다. 회로는 각 모드에 QFT를 적용해 운동량기저로 이동하고, OutPoly로 $(x-K/2)^2$을 보조 레지스터에 계산한다. 이 값과 운동에너지 계수로 위상을 누적한 뒤 중간 계산을 되돌리고 역 QFT로 위치기저에 복귀한다.
퍼텐셜: 상태별 계수를 QROM으로 불러와 위상으로 바꾼다
퍼텐셜 계수는 전자상태에 따라 달라진다. LoadCoeffsKDC는 QROM을 사용해 현재 전자상태에 맞는 미리 계산된 비트열을 계수 레지스터로 적재한다. 여기서 QROM은 외부 데이터베이스에 양자적으로 접속하는 하드웨어 QRAM이 아니다. 고정된 coefficient table을 가역 회로로 조회하는 논리 연산이다.
선형항은 $Q_r$, 이차항은 $Q_rQ_{r'}$를 산술회로로 계산한다. 그 결과를 재사용 가능한 $b$-qubit 위상 기울기 상태
$$ |R_b\rangle=2^{-b/2}\sum_{y=0}^{2^b-1}e^{i2\pi y/2^b}|y\rangle $$에 더해 에너지와 시간에 해당하는 위상을 만든다. $b$는 회전의 고정소수점 정밀도를 정한다. 위상을 적용한 뒤 QROM과 산술 중간값을 uncompute해 보조 레지스터를 원상복구한다.
비대각 전자 결합: XOR fragment와 Clifford 대각화
전자상태 사이의 결합은 전자 부분에서 비대각 행렬을 만든다. 기반 알고리즘은 두 전자상태 인덱스의 XOR 차이가 같은 항을 하나의 fragment로 묶는다. 비트열이 한 자리에서 다르면 Hadamard로, 여러 자리에서 다르면 CNOT으로 차이를 한 자리에 모은 뒤 Hadamard로 block diagonalization한다. 이 과정은 Clifford gate만 사용한다. 각 fragment가 대각형이 된 뒤에야 위치별 퍼텐셜 위상을 효율적으로 적용할 수 있다.
5. [원 데모] 코드가 실제로 계산한 것
실행된 예제의 설정은 다음과 같다.
| 변수 | 값 | 의미 |
|---|---|---|
| electronic states | 2 | 1-qubit 전자상태 레지스터의 $\lvert 0\rangle$, $\lvert 1\rangle$ |
| vibrational modes | 1 | 정상모드 하나 |
| $k$ | 2 | 모드당 큐비트 수 |
| $K=2^k$ | 4 | 진동 좌표 격자점 수 |
| $b$ | 5 | 위상·계수 정밀도 비트; delta=0.04에서 계산 |
| $\omega$ | 1 | 예시 진동수 |
| $dt$ | 0.4 a.u. | 한 시간 단계 |
| steps | 10 | 총 4 a.u.까지 평가 |
| 퍼텐셜 차수 | $\alpha=1$ | 실제 실행은 선형 좌표항만 사용 |
| 계수 배열 | [[1.0, 0.0], [-1.3, 1.3]] |
실제 분자가 아닌 예시값 |
레지스터의 합은 23 wires다.
| 레지스터 | wires |
|---|---|
| 전자상태 | 1 |
| 진동 격자 | 2 |
| phase gradient | 5 |
| coefficients | 5 |
| scratch | 6 |
| cache | 4 |
| 합계 | 23 |
물리 모형 자체는 전자 1개와 진동 2개, 즉 3 wires에 부호화된다. 나머지 20개는 정밀도와 가역 산술을 위한 작업공간이다. 따라서 이를 “23-qubit 실제 분자 계산”이라고 부르면 부정확하다.
초기 핵 상태는 조화진동자 바닥상태를 흉내 낸 이산 Gaussian 파동묶음이다. 프로그램은 각 $t=0,\ldots,10$에서 회로를 처음부터 다시 상태벡터로 계산하고 전자 레지스터의 qp.probs를 반환한다. shots를 지정하지 않았으므로 유한 표본의 통계가 아니라 상태벡터에서 얻은 exact probabilities다. 페이지에 보고된 전체 스크립트 실행시간은 약 8분 9.7초다.
출력 그래프는 두 전자상태의 population이 진동하고 전달이 완전하지 않음을 보여준다. 두 곡선의 합이 1이라는 것은 unitary evolution과 두 상태 레지스터의 정규화에 따른 필요조건이다. 정확한 고전 격자 전파, MCTDH, Trotter 수렴 또는 실험과 비교하지 않았으므로 chemical accuracy의 검증은 아니다.
6. [계산 경계] 데모가 입증하지 않은 것
| 층 | 데모가 한 일 | 하지 않은 일 |
|---|---|---|
| 장치 | CPU 상태벡터 시뮬레이션 | 실제 QPU 또는 광자 칩 실행 |
| 측정 | exact electronic probabilities | finite-shot 오차와 측정 비용 검증 |
| 회로 | PennyLane 고수준 산술 연산 | 특정 하드웨어 native gate 컴파일 |
| 오차 | 작은 격자·고정 정밀도·Trotter 근사 사용 | 오차항별 수렴 연구와 총 error budget |
| 화학 | 임의 계수의 2-state·1-mode LVC toy | 실제 분자 parameterization과 고전 기준값 |
| OLED | 일반 spin-vibronic 확장 가능성 | SOC, singlet–triplet, ISC/RISC, host 환경 |
| 성능 주장 | 알고리즘 기능 설명 | 양자 가속 또는 양자우위 입증 |
실제 연구에서는 최소한 다음 오차를 분리해야 한다.
- 포함한 전자상태와 진동모드의 선택
- LVC/QVC 및 비조화 퍼텐셜의 절단
- 전자구조 parameterization과 diabatization
- 위치 격자 $K$와 범위
- Trotter 시간간격과 차수
- 고정소수점 위상 정밀도 $b$
- 상태 준비와 유한-shot 측정
- 하드웨어 noise, 오류정정과 논리→물리 자원 변환
7. [기반 논문] 실제 연구 규모의 자원 추정
최신 arXiv v3와 저널 논문은 모드당 $K=16$ 격자점과 2차 product formula를 사용해 다음 algorithmic/logical 자원을 추정한다.
| 모형 | states · modes | 전파 조건 | qubits | Toffoli gates |
|---|---|---|---|---|
| $(\mathrm{NO})_4$-Anth | 5 · 19 | 100 fs, 10% | 146 | $5.47\times10^6$ |
| $(\mathrm{NO})_4$-Anth | 5 · 19 | 100 fs, 1% | 146 | $1.73\times10^7$ |
| anthracene dimer | 6 · 21 | 100 fs, 1% | 154 | $2.76\times10^6$ |
| anthracene dimer | 6 · 21 | 500 fs, 1% | 154 | $3.54\times10^7$ |
| anthracene/$\mathrm{C}_{60}$, reduced | 4 · 11 | 100 fs, 1% | 113 | $6.62\times10^5$ |
| anthracene/$\mathrm{C}_{60}$, full | 4 · 246 | 100 fs, 1% | 1,053 | $2.66\times10^7$ |
이 표는 실행 결과가 아니다. 물리 큐비트 수, 오류정정 code distance, magic-state factory, gate cycle time과 wall-clock time을 포함하지 않는다. 154 logical/algorithmic qubits를 현재의 154개 물리 큐비트 장치에서 실행할 수 있다는 뜻이 아니다.
또한 버전을 섞어 읽으면 안 된다. arXiv v1은 full 246-mode anthracene/$\mathrm{C}_{60}$에 1,065 qubits와 $2.7\times10^9$ Toffoli를 제시했지만, v3는 empirical error extrapolation을 사용해 1,053 qubits와 $2.66\times10^7$ Toffoli로 갱신했다. 이는 QPU 성능 향상이 아니라 자원 산정법과 알고리즘 분석의 변경이다.
별도 PennyLane 자원 추정 데모는 희소성을 활용하지 않는 dense $(\mathrm{NO})_4$-Anth 상한으로 146 wires, $6.082\times10^7$ Toffoli, 총 $5.21\times10^8$ gates를 보고한다. 구조와 희소성을 어떻게 쓰는지가 feasibility를 크게 바꾼다는 뜻이다.
8. VQE와는 무엇이 다른가
| 구분 | VQE | 이번 알고리즘 |
|---|---|---|
| 질문 | 낮은 에너지 고유상태와 에너지는 무엇인가 | 광여기 뒤 상태 population은 시간에 따라 어디로 가는가 |
| 계산 | 매개변수화 상태의 에너지 최소화 | $e^{-iHt}$의 실시간 전파 |
| 출력 | 에너지와 정적 관측량 | $P_j(t)$, 전달 경로, 상관함수·스펙트럼 후보 |
| 회로 | ansatz, Pauli 측정, 고전 최적화 | QFT, QROM, 가역 산술, phase gradient, Trotter |
| 장치 지향 | NISQ 얕은 회로 연구가 많음 | 깊은 산술을 포함한 fault-tolerant 지향 |
VQE로 얻은 에너지는 동역학 해밀토니안의 입력 일부가 될 수 있다. 그러나 정적 에너지 하나가 진동이 유도하는 시간별 population transfer를 자동으로 주지는 않는다.
9. [리뷰 해석] TADF·PhOLED 연구에 연결되는 지점
Thermally activated delayed fluorescence(TADF)의 reverse intersystem crossing(RISC)은 $\Delta E_{\mathrm{ST}}$ 하나로 결정되지 않는다. singlet와 triplet의 state character, higher triplet, SOC, 특정 promoting mode, 주변 유전환경과 구조 변동이 함께 관여할 수 있다. KDC 형식은 원리상 퍼텐셜 행렬에 SOC 해밀토니안을 더해
$$ \mathbf W'(\mathbf Q)=\mathbf W(\mathbf Q)+\mathbf H_{\mathrm{SO}}(\mathbf Q) $$와 같은 spin-vibronic 모형으로 확장할 수 있다. 기반 논문도 TADF를 잠재 응용으로 언급한다.
| OLED 물리 | 진동-전자 모형의 대응 요소 |
|---|---|
| $S_1$, $T_1$, $T_n$ 에너지·성격 | diabatic electronic-state manifold |
| SOC | $\mathbf H_{\mathrm{SO}}$와 상태 간 결합 |
| promoting vibration | $a_r^{(i,j)}$, $b_{rr'}^{(i,j)}$ |
| ISC/RISC 경로 | singlet–triplet population transfer |
| 시간별 상태점유 | $P_{S_1}(t)$, $P_{T_1}(t)$, $P_{T_n}(t)$ |
| 스펙트럼 | 시간 상관함수의 Fourier transform |
그러나 이 연결은 리뷰의 확장 해석이지 현재 데모의 결과가 아니다. 데모에는 SOC도 triplet도 없다. 실제 TADF·PhOLED 박막에는 열적 진동분포, dephasing, host polarization, disorder, 분자간 모드와 개방계 relaxation이 중요하다. 닫힌 계 toy model의 coherent oscillation을 비가역적인 ISC/RISC rate로 바로 읽을 수 없다.
따라서 이 알고리즘은 가까운 시기의 high-throughput generator보다, 미래 fault-tolerant 환경에서 소수의 고가치 후보를 검증하는 동역학 커널로 보는 편이 타당하다.
10. OLED 역설계에 넣는다면: 제한된 검증형 workflow
현실적인 연결 순서는 다음과 같다.
- 생성모델 또는 데이터베이스에서 후보를 만든다.
- 합성가능성, 안정성, 기본 전자구조를 고전 방법으로 선별한다.
- DFT/TDDFT·다중참조 계산으로 상태, Hessian, 정상모드와 SOC를 얻는다.
- diabatization을 수행하고 작은 LVC/QVC spin-vibronic 해밀토니안을 만든다.
- exact grid, MCTDH 또는 tensor-network 기준선으로 먼저 검증한다.
- 동일 정확도 목표에서 Trotter·격자·정밀도 오차와 logical resource를 추정한다.
- 고전 기준이 어려워지는 제한된 후보에서만 양자 전파를 고려한다.
- population, mode sensitivity와 불확실성을 ML surrogate의 기작 중심 라벨로 되돌린다.
첫 PoC는 $S_1$, $T_1$과 한두 개 higher triplet, 5–10개 promoting mode 정도의 공개 해밀토니안으로 시작하는 편이 낫다. 모든 trajectory를 exact 또는 MCTDH 기준과 비교하고, $K$, $dt$, $b$를 독립적으로 변화시켜 오차를 분해해야 한다. 이 검증을 통과한 뒤에야 larger-state 또는 open-system 확장을 논의할 수 있다.
최종 평가
PennyLane 데모는 진동-전자 해밀토니안을 양자회로의 구체적인 구성요소로 번역한 좋은 구현 튜토리얼이다. QFT가 운동에너지에, QROM과 가역 산술이 상태·좌표 의존 퍼텐셜에, phase-gradient state가 시간진화 위상에, Clifford fragment diagonalization이 다전자상태 결합에 각각 어떤 역할을 하는지 코드 수준에서 볼 수 있다.
과학적 주장은 작게 유지해야 한다. 23-wire CPU 상태벡터 toy simulation이며, 실제 분자나 QPU, quantitative dynamics 또는 quantum advantage를 검증하지 않았다. 기반 논문의 수백 logical qubits와 수백만 Toffoli도 물리 하드웨어 비용이 아니라 fault-tolerant 알고리즘 수준의 추정이다.
그럼에도 연구 방향은 분명하다. OLED 분자설계가 정적 $\Delta E_{\mathrm{ST}}$와 oscillator strength만 보는 단계에서 벗어나, 상태별 population, higher-state mediation, SOC와 promoting mode의 시간동역학을 기작 중심 라벨로 사용하려면 이런 실시간 진동-전자 전파가 필요하다. 그 가능성을 현실적인 연구도구로 바꾸는 데에는 더 깊은 회로만큼이나 신뢰할 수 있는 해밀토니안 구성과 엄격한 고전 기준선이 중요하다.
근거 자료
- Emily Nobes, A quantum algorithm for vibronic dynamics, PennyLane Demos, 19 August 2026.
- D. Motlagh et al., Quantum Algorithm for Vibronic Dynamics: Case Study on Singlet Fission Solar Cell Design, Quantum Science and Technology 10, 045048 (2025).
- D. Motlagh et al., arXiv:2411.13669v3, current technical version.
- D. Dhawan, Quantifying resource requirements for vibronic dynamics simulation, PennyLane Demos, updated 27 May 2026.
- PennyLaneAI, official demo source at fixed commit.
- PennyLane,
lightning.qubitdocumentation.
검증 메모: 수치와 실행 조건은 PennyLane 공개 코드, 공식 문서, 동료평가 논문과 최신 arXiv v3를 대조했다. 데모 결과, 논문 자원 추정, OLED 적용에 관한 리뷰 해석을 서로 다른 증거 층으로 유지했다.