[논문 리뷰] Cis--Trans Rotational Isomerism of Seleno-, Thio-, and Formic Acids and Their Dimers: Chemical Kinetics under Interstellar Conditions

Judith Würmel, John M. Simmie|arXiv (Cornell University)|2026. 01. 22.

Astrophysics and Star Formation Studies인용 수 0

한 줄 요약

이 논문은 CVT/SCT/TST를 사용하고 ISPE 보정을 적용하여 10 K에서 300 K까지 포름산(formic acid) 및 그 thio- 및 seleno- 유도체의 cis–trans 회전 이성질화에 대한 기상상 단계의 속도 상수와 에너지 장벽 데이터를 제공합니다(단량체 및 이량체).

ABSTRACT

Tunnelling reactions of molecules embedded on cryogenic noble-gas matrices are being used in fundamental studies of how reactivity varies with the nature of the supposedly inert matrix as well as pointers to the chemistry occurring in the interstellar medium on ice-grains. To these ends we present chemical kinetic rate constants for the extit{cis} to extit{trans} isomerisation of seleno-, thio- and monomeric formic acids and that of their three dimeric species, based on multidimensional calculations in the gas-phase, from 10~K to 300~K as a guide to the matrix reactions.

연구 동기 및 목표

• 포름산 및 유사체의 cis–trans 회전 이성질화가 항성간의 냉각 조건에서 어떻게 이해되는지에 대한 동기를 부여한다.
• cryogenic 매트릭스 실험 및 ISM 화학의 해석을 안내하기 위한 정확한 기상상 속도 상수와 장벽 데이터를 제공한다.
• 단량체와 이량체 형태가 H, S, Se 유사체의 터널링 지배적 이성질화에 미치는 영향을 평가한다.

제안 방법

• B2PLYP-D3BJ/def2-TZVP로의 기하학 최적화 및 고유 진동수 계산과 고전적 및 비조화 보정.
• IRC 경로에 대한 DLPNO-CCSD(T)/def2-TZVP/C를 통한 고수준 단일점 에너지 계산.
• CVT를 SCT, QRC 및 ISPE 보정으로 보완한 표준 변이된 전이상 이론으로 반응 속도 상수 계산.
• 저수준 IRC(B2PLYP-D3BJ/def2-TZVP)에서의 2단계 접근 및 IRC를 따라가는 고수준 CCSD(T) 보정.
• BSSE 보정 영향 완화 없이 B3LYP-D3BJ/def2TZVP(및 해당되는 경우 상위 DLPNO-CCSD(T))로 이량체 에너지 계산.

실험 결과

연구 질문

• RQ1 가스상에서 HC(O)OH, HC(O)SH, HC(O)SeH의 cis→trans 이성질화의 장벽과 반응 엔thal피는 무엇인가?
• RQ2 10 K에서 300 K까지의 가스상 조건에서 단량체 및 이량체 포름산의 회전 이성질화가 온도에 따라 어떻게 변하는가?
• RQ3 O, S, Se 계열에서 양자 터널링 효과(SCT, QRC)가 저온에서 이성질화에 지배적인가?
• RQ4 Computed gas-phase rate constants가 매트릭스 고립 및 ISM 관측과 cis–trans 전환과 비교될 수 있는가?
• RQ5 이러한 동역학이 성간 얼음과 입자에서 cis 형태의 생존 및 탐지에 어떤 함의를 가지는가?

주요 결과

Reaction	E‡ (kJ mol−1)	ΔrH (kJ mol−1)
cis → trans HC(O)OH	31.84	−17.99
cis → trans HC(O)SH	35.56	−4.06
cis → trans HC(O)SeH	30.86	−1.38
tc1 → tt2 HC(O)OH	38.80	−13.79
tc1 → tt2 HC(O)SH	34.29	−9.67
tc1 → tt2 HC(O)SeH	34.98	−0.41
tc4 → tt3 HC(O)OH	34.24	−13.95
tc4 → tt3 HC(O)SH	39.61	−3.40
tc4 → tt3 HC(O)SeH	34.64	−0.82
cc5 → tc3 HC(O)OH	36.51	−15.53
cc5 → tc3 HC(O)SH	41.08	−3.11
cc5 → tc3 HC(O)SeH	29.00	−0.51

• Cis→trans 장벽은 HC(O)OH: 31.84 kJ/mol로 ΔrH = −17.99 kJ/mol; HC(O)SH: 35.56 kJ/mol로 ΔrH = −4.06 kJ/mol; HC(O)SeH: 30.86 kJ/mol로 ΔrH = −1.38 kJ/mol 입니다.
• 이량체 전이는 tc1→tt2 및 tc4→tt3에서 더 높은 장벽을 보이고, cc5→tc3은 전이상 구성이 상대적으로 더 폭넓고 결합이 다르게 이루어집니다.
• 10 K에서 모든 종에 대해 터널링이 이성질화의 지배적 요인; 단량체 HC(O)OH의 cis→trans 반감기는 약 7.6 ms로 예측되나 HC(O)SH 및 HC(O)SeH는 상당히 더 긴 편(어떤 이량체의 경우 ~136 시간 등).
• 본 연구의 속도 상수는 일반적으로 이전의 가스상 추정치보다 더 큰 경향을 보임(포름산의 경우 최대 ~3배, 티오-/셀레노- 유도체의 경우 최대 ~36배까지), 포함된 터널링 보정 반영.
• cis-HCOOH는 냉각 조건에서 빠른 이성질화로 인해 ISM에서 실질적으로 관측 불가일 가능성이 크며( Keq ≫ 10^15 at low T), cis- 포름산의 관측은 대체 형성 경로나 매트릭스/표면 효과를 시사합니다.
• 포름산 이량체는 터널링 지배적 역전 현상을 보이며 BSSE 보정의 질적 경향에 미치는 영향이 미미하며, 황(S) 및 셀레늄(Se) 이량체는 산소 유사체에 비해 몇 제곱의 차이로 느립니다.