[论文解读] Structure and thermochemistry of K$_2$Rb, KRb$_2$ and K$_2$Rb$_2$

Jason N. Byrd, John A. Montgomery+1 more|arXiv (Cornell University)|Mar 23, 2010

Cold Atom Physics and Bose-Einstein Condensates被引用 7

一句话总结

本研究采用CCSD(T)方法结合大基组和CBS外推，对K₂Rb、KRb₂和K₂Rb₂进行了高精度从头算结构与热化学计算。研究识别出KRb + KRb与K₂ + Rb₂之间无能垒的反应路径，其中D₂h和Cₛ对称性为全局最小值，并报告了相应的原子化能分别为11,131 cm⁻¹和11,133 cm⁻¹，为超冷分子化学提供了关键数据。

ABSTRACT

The formation and interaction of ultracold polar molecules is a topic of active research. Understanding possible reaction paths and molecular combinations requires accurate studies of the fragment and product energetics. We have calculated accurate gradient optimized ground state structures and zero point corrected atomization energies for the trimers and tetramers formed by the reaction of KRb with KRb and corresponding isolated atoms. The K${}_2$Rb andKRb${}_2$ trimers are found to have global minima at the $C_{2v}$ configurationwith atomization energies of 6065 and 5931 cm${}^{-1}$ while the tetramer is found to have two stable planar structures, of $D_{2h}$ and $C_s$ symmetry, which have atomization energies of 11131 cm${}^{-1}$ an 11133 cm${}^{-1}$, respectively. We have calculated the minimum energy reaction path for the reaction KRb+KRb to K$_2+$Rb$_2$ and found it to be barrierless.

研究动机与目标

• 确定与超冷分子碰撞相关的异核碱金属三聚体和四聚体（K₂Rb、KRb₂、K₂Rb₂）的基态结构与原子化能。
• 研究KRb + KRb → K₂ + Rb₂反应的反应路径与能量学，特别是其无能垒特性。
• 为模拟涉及KRb分子的超冷化学反应提供精确的热化学数据。
• 通过比较费米子与玻色子KRb体系中的反应速率，评估自旋统计在四体交换反应中的作用。

提出的方法

• 在CCSD(T)/def2-TZVPP和CCSD(T)/def2-QZVPP水平上进行从头算电子结构计算，通过Rb的ECP28MDF考虑核心-价相关效应与相对论效应。
• 对几何结构进行优化，并计算振动频率以确认为能量极小值，同时提取零点振动能（ZPE）校正。
• 采用显式相关CCSD(T)-F12b方法提高电子相关能的精度，并应用CBS外推法估算完全基组极限。
• 在CCSD(T)/def2-TZVPP水平上计算最小能量反应路径，单点能量在CCSD(T)-F12b水平和CBS极限下进行评估。
• 计算内禀反应坐标（IRC），以定位过渡态并确认D₂h与Cₛ极小值之间的反应路径。
• 最小化长程静电相互作用，以确定反应路径的初始双原子-双原子构型。

实验结果

研究问题

• RQ1 在高精度理论水平下，K₂Rb、KRb₂和K₂Rb₂的最稳定几何构型与原子化能分别是什么？
• RQ2 连接KRb + KRb与K₂ + Rb₂的反应路径上是否存在能垒？过渡态的性质如何？
• RQ3 自旋统计（费米子与玻色子KRb）如何影响超冷陷阱中四体交换反应的速率？
• RQ4 标量相对论效应与核心-价相关效应在这些异核碱金属团簇热化学中起什么作用？
• RQ5 长程静电相互作用在预测双原子-双原子缔合的初始构型方面表现如何？

主要发现

• K₂Rb与KRb₂三聚体的全局最小值位于C₂v对称性，ZPE校正后的原子化能分别为6065 cm⁻¹与5931 cm⁻¹。
• K₂Rb₂四聚体表现出两种稳定的平面结构：D₂h与Cₛ，其原子化能分别为11,131 cm⁻¹与11,133 cm⁻¹。
• 反应KRb + KRb → K₂ + Rb₂为无能垒过程，D₂h结构为产物侧的最小值。
• 该反应的过渡态位于D₂h解离能以上1167.3 cm⁻¹处，证实了反应路径的平滑性。
• CBS外推的CCSD(T)-F12b计算得到高度精确的原子化能，其中K₂Rb₂的Cₛ结构表现出最高稳定性，为11,133 cm⁻¹。
• 本研究证实，自旋统计会抑制费米子⁴⁰KRb的s波散射，这可能通过抑制四体损失而稳定陷阱，从而实现对交换反应的受控研究。