Cotunneling-Assisted Redistribution of Nonpassivity in a Photosynthetic Junction

기술 요약: 광합성 접합부 내 비수동성(Nonpassivity)의 코터널링 보조 재분배

문제 정의
본 논문은 비평형 분자 접합부로 모델링된 광계 II 반응 중심(PSIIRC)의 열역학적 특성을 조사한다. PSIIRC는 전하 이동 네트워크로서 기능하는 것으로 알려져 있으나, 유한한 비평형 시스템으로서의 거동은 다체(many-body) 에너지 준위 전반에 걸친 복잡한 인구 분포 재분배를 수반한다. 핵심적인 문제는 상관된 다전자 수송 과정, 특히 **코터널링(cotunneling)**이 시스템의 정상 상태(steady state)의 **수동성(passivity)**에 어떻게 영향을 미치는지를 이해하는 것이다. 수동성이란 에너지 준위에 따른 상태 인구(population)의 정렬을 의미하며, 비수동적(능동적) 상태는 낮은 에너지 상태가 높은 에너지 상태보다 인구가 높거나 같아야 한다는 조건을 위반함으로써 추출 가능한 일을 포함한다. 저자들은 쿨롱 상호작용과 코터널링 경로가 이러한 인구 정렬을 어떻게 재구성하고, 기존의 전류나 평균 에너지와 같은 전통적인 수송 관측량에서 드러나는 것 이상으로 시스템의 일 수행 능력을 어떻게 변화시키는지 규명하고자 한다.

방법론
본 연구는 산소 발생 복합체(OEC)를 소스(source)로, 하류의 플라스토퀴논(PQ) 풀을 드레인(drain)으로 취급하여 PSIIRC를 미시적인 소스-시스템-드레인 기술법으로 모델링한다.

• 모델 구축: 시스템은 기저 상태, 6개의 비국소화된 엑시톤 상태, 3개의 중간 전하 이동 상태, 라디칼 쌍 상태($|\alpha_1\rangle$), 말단 양이온 상태($|\beta_1\rangle$), 그리고 코터널링에 의해 보조되는 이중 환원 수용체 상태($|\alpha_2\rangle$)를 포함하는 13개 상태의 유한 양자 시스템으로 모델링된다. $|\alpha_2\rangle$ 상태는 수용체 측의 축소된 전자 부공간으로부터 미시적으로 구축되었으며, 구체적으로 터미널 퀴논 $Q_B$($|Q_B^{2-}\rangle$, $|Q_A^- Q_B^-\rangle$, $|Pheo_{D1}^- Q_B^-\rangle$)를 포함하는 2전자 구성(two-electron configurations)을 나타낸다.
• 역학(Dynamics): 축소된 역학은 나카지마-즈완직(Nakajima-Zwanzig) 투영 연산자 프레임워크를 사용하여 정식화된다.
  • 광 여기 및 완화는 2차 커널로 처리된다.
  • 엑시톤 매니폴드 내의 진동 재분배는 수정된 레드필드(Redfield) 처리를 사용한다.
  • 엑시톤 매니폴드와 전하 이동 매니폴드 사이의 내부 전이는 포르스터(Förster) 영역에서 처리된다.
  • 리저버(reservoir)와의 순차적 전자 교환은 2차까지 계산된다.
  • 코터널링 과정(상관된 2전자 전이)은 분자-리저버 결합의 4차까지 유지된다. 이 속도($\omega_{xy}$)는 쿨롱 상호작용 에너지 $U$에 의존하는 정규화된 유효 속도이다.
• 열역학적 정량화: 정상 상태 인구는 다음을 계산하는 데 사용된다:
  • 에르고트로피(Ergotropy, $\mathcal{E}$): 실제 비평형 정상 상태와 그 수동적 대응 상태(에너지를 최소화하도록 인구가 재정렬된 상태) 사이의 에너지 차이로 정의되는 추출 가능한 일.
  • 에너지 커패시턴스(Energetic Capacitance, $C$): 동일한 인구 스펙트럼에 대해 수동 상태와 안티-패시브(anti-passive) 상태(에너지를 최대화하도록 인구가 재정렬된 상태) 사이의 총 에너지 범위.
  • 평균 에너지($\langle E \rangle$): 정상 상태의 평균 에너지.

주요 결과

1. 비수동성의 재분배: 코터널링은 추출 가능한 비평형 에너지를 유의미하게 재분배한다. 낮은 쿨롱 상호작용($U$)에서, 비수동적 에너지는 $D1$ 라디칼 쌍($|\alpha_1\rangle$), 상관된 수용체 상태($|\alpha_2\rangle$), 그리고 $D2$ 분지 특성이 강한 엑시톤 상태($X_2, X_5$)에 의해 공유된다. 높은 $U$에서, $D1$ 라디칼 쌍의 기여는 부차적이 되는 반면, $|\alpha_2\rangle$ 상태와 $D2$ 가중 엑시톤이 에르고트로피의 주요 기여자로 남는다.
2. 열역학적 매끄러움 대 미시적 재구성: 본 연구는 미시적인 인구 정렬이 불연속적인 변화를 겪더라도 열역학적 관측량(에르고트로피 및 커패시턴스)은 매끄럽고 연속적인 거동을 보일 수 있음을 발견했다. 즉, "인구 교차(population crossings)"(두 상태의 상대적 인구 크기가 바뀌는 현상)가 수동 상태의 순열을 변경할 수 있지만, 지배적인 에너지 가중 인구 차이가 구조적으로 변하지 않는 한 총 에르고트로피에는 급격한 변화를 일으키지 않는다.
3. 분자 섹터의 구별된 역할: 분석 결과, 운동학적 참여도, 정상 상태 인구 크기, 그리고 열역학적 중요성은 서로 구별되는 속성임이 밝혀졌다. 예를 들어, $|\alpha_2\rangle$ 상태는 인구가 작고 코터널링 속도가 증가함에 따라 더욱 감소할 수 있음에도 불구하고, 그 에너지 가중 편차가 에르고트로피의 지배적인 항으로 남을 수 있다. 반대로, 인구가 매우 높은 상태라도 그 점유율이 이미 수동 정렬에 가깝다면 추출 가능한 일에 거의 기여하지 않을 수 있다.
4. 조각별 선형 관계(Piecewise-Linear Relations): 본 논문은 에르고트로피, 에너지 커패시턴스, 그리고 평균 에너지 사이의 조각별 선형 관계를 식별한다. 이러한 선형 구간은 고정된 인구 정렬에 대응한다. 인구 교차에서 꺾임(kink)이나 끊김이 발생하며, 이는 미시적 재구성을 나타내는 열역학적 징후 역할을 한다.
5. 상호작용 의존적 응답:
   • 낮은 $U$: 전자 방출 속도가 증가함에 따라 에르고트로피는 증가하는 반면, 터미널 전이와 관련된 일은 감소한다. 이는 터미널 작업 손실을 보상하기 위해 인구 재분배가 전역적 비수동성을 강화하는 네트워크 응답을 시사한다.
   • 높은 $U$: 방출 속도가 증가함에 따라 에르고트로피와 터미널 작업 모두 감소한다. 그러나 에너지 커패시턴스는 계속 증가하는데, 이는 실현된 정상 상태는 더 수동적으로 변하더라도, 기저의 인구 스펙트럼은 그 잠재적인 에너지 범위를 유지하거나 오히려 증가시킨다는 것을 나타낸다.

의의 및 주장
본 논문은 PSIIRC에 대한 완전한 열역학적 기술이 단순히 전통적인 수송 관측량에 의존하는 것이 아니라, 전체 운동학적 네트워크에 의해 생성된 인구 정렬을 반드시 고려해야 한다고 주장한다.

• 진단 도구로서의 수동성: 수동성 분석은 전류, 평균 에너지 또는 개별 인구 곡선에서는 숨겨져 있는 정상 상태의 재구성을 드러낸다. 이는 인구 크기만으로는 추론할 수 없는, 어디에 열역학적 중요성이 존재하는지를 식별하는 분자적 기준을 제공한다.
• 코터널링의 역할: 코터널링은 $D1$ 전하 분리 부분과 $D2$ 엑시톤 부분에 분산된 분자 인구의 열역학적 중요성을 재편하는 것으로 나타났다. 이는 강한 쿨롱 상호작용 하에서 추출 가능한 에너지의 균형을 생산적인 $D1$ 라디칼 쌍에서 터미널 수용체 및 $D2$ 엑시톤 매니폴드로 이동시킨다.
• 방법론적 기여: 본 연구는 에르고트로피와 에너지 커패시턴스의 결합이 실현된 비수동성의 감소와 인구 스펙트럼의 가용 에너지 범위 감소를 어떻게 구분할 수 있는지 보여준다. 에르고트로피와 에너지 커별 사이의 조각별 선형 관계는 코터널링에 의해 유도되는 연속적인 미시적 인구 재구조화의 직접적인 열역학적 징후 역할을 한다.

저자들은 이러한 현상이 색소 부위 파라미터화의 합리적인 변동에 대해 견고하며, 쿨롱 상호작용, 다체 수송, 그리고 광합성 접합부의 비평형 구동 사이의 상호작용에 내재된 현상이라고 결론짓는다.