主要作者与研究机构 本研究题为《s-四嗪在其水合团簇和稀溶液中最低¹(n, π*)激发的氢键和溶剂化显色位移》，由来自南京大学化学系、理论与计算化学研究所及介观化学重点实验室的Daiqian Xie，以及澳大利亚La Trobe大学生物化学系和Cytopia有限公司的Jun Zeng共同完成。该研究成果发表于2004年的《Journal of Computational Chemistry》期刊第25卷第12期。
研究背景与目的 本研究聚焦于一个重要的科学问题：含氮杂环（azine）及其衍生物（如核酸碱基）与水分子间的氢键相互作用。这种特定的溶质-溶剂相互作用在极性介质和蛋白质等许多重要体系中普遍存在，并对生色团的电子光谱产生显著的溶剂化显色位移（solvatochromatic shift），在溶液和蛋白质的电子转移过程中扮演着关键角色。此前，作者团队已开发了一套结合计算机模拟与静电相互作用分析的方法，用于计算周围介质对分子电子性质的影响，并成功应用于一系列单嗪、二嗪和三嗪在水中的最低¹(n, π)激发态溶剂位移的研究。本研究的目的，正是为了完成该系列研究的最后一块拼图：探究s-四嗪（s-tetrazine）这一具有四个氮原子的特殊嗪类分子，其基态与激发态与水分子间特定的氢键相互作用，并阐释其稀水溶液中最低¹(n, π)吸收和荧光谱带中心的溶剂位移。传统的观点认为，基态(S₀)中氮原子的孤对电子轨道与水分子形成的氢键在(n, π*)激发态(S₁)会因电子被激发而消除，从而导致大的吸收蓝移和较小的荧光位移。然而，近期对吡啶-水团簇的高精度计算则提出了激发态氢键可能依然存在的新机制，这使得对嗪类-水氢键在不同电子态下行为的研究变得更为复杂和必要。
详细研究流程 该研究采用了一套系统的计算化学流程，主要包含三个核心步骤：液相模拟、电子结构计算和溶剂位移计算。 首先进行研究对象的构建与气相电子结构计算。研究对象为孤立的s-四嗪分子及其与水分子形成的团簇。研究人员采用GAMESS程序，在观测到的基态平衡几何构型下，使用三zeta价层加极化（TZVP）基组，通过完全活性空间自洽场（CASSCF）方法，计算了s-四嗪分子基态(S₀)和第一¹(n, π*)激发态(S₁)的性质。计算选定的活性空间能很好地描述该激发态，计算得到的激发能为2.38 eV，与实验观测值（2.25-2.35 eV）高度吻合。基于CASSCF波函数，研究进一步在分子范德华壳层外的特定区域计算了静电势（ESP），并通过拟合这些静电势，构建了描述分子电荷分布的模型。这个模型特别之处在于，它包含三项组分来描述库仑相互作用：原子点电荷、位于每个氮原子上的原子点偶极（用于显式描述局域s-p杂化对分子偶极矩的贡献），以及八个名为X9-X16的额外浮动点电荷和偶极（用于模拟π环的电荷密度）。 其次，开展液相蒙特卡洛模拟。为了研究s-四嗪在稀水溶液中的行为，研究人员构建了一个包含1个s-四嗪分子和102个TIP3P水分子模型的体系，在298K和1个大气压条件下，进行了等温等压（NPT）系综的刚体蒙特卡洛模拟。溶液体系采用周期性截角八面体边界条件。模拟中使用了改进的Kollman势函数形式来描述分子间的成对加和势能，其中伦纳德-琼斯（Lennard-Jones）参数在基态和激发态保持一致，库仑相互作用则采用了上述由从头算拟合得到的三组分电荷模型。模拟经过至少1亿次移动的充分平衡后，进行了2千万次移动的采样，以获取体系的热力学性质、径向分布函数（RDF）以及用于后续溶剂位移计算的溶液构型快照。 最后，执行溶剂位移计算。该研究使用了作者团队早前开发的溶剂位移方法。此方法从液相蒙特卡洛模拟中采样的构型出发，将溶剂位移推导为生色团在初始态和最终态下与溶剂静电相互作用的变化。在此步骤中，溶剂分子被显式处理，并使用了可极化的溶质分子模型。为了计算垂直激发的溶剂位移，还使用了球边界条件和Friedman镜像电荷技术来精确处理长程静电相互作用。该计算充分考虑了溶质在基态和激发态下的极化率张量差异，这些极化率是通过HONDO程序在自洽场（SCF）水平上计算得到的。
主要研究结果与论证 研究结果清晰地揭示了s-四嗪水合团簇和稀溶液的氢键结构与光谱位移特征。 在s-四嗪-水团簇层面，计算结果预测了两种稳定的异构体。第一种是线性（linear）氢键构型，水分子位于s-四嗪分子平面内，与四个氮原子之一形成直线的O-H···N氢键。该构型在基态的相互作用能为-6.59 kcal/mol，在激发态为-6.17 kcal/mol，二者非常接近，表明在离域的(n, π*)电子激发过程中，氢键结构变化很小。第二种是分叉（bifurcated）氢键构型，水分子同时与两个相邻的氮原子形成氢键，具有C₂ᵥ对称性。此构型在基态是能量为-5.32 kcal/mol的稳定态，但在激发态则变为一个过渡态，其键合强度下降至-4.04 kcal/mol。这些构型得到了MP2/cc-pVDZ级别的从头算优化的验证，尽管高精度计算显示，线性结构的最优构型中水分子会垂直于四嗪平面，类似于吡啶-水复合物的结构，且分叉结构在激发态也可能是稳定的。此外，由于基态与激发态团簇的能量最低点几何结构相似，从弗兰克-康登（Franck-Condon）区域出发的分子动力学轨迹表现出简单的周期性运动，这与其它嗪类-水团簇的动力学行为一致。 在s-四嗪-水稀溶液层面，径向分布函数和热力学数据揭示了更丰富的溶剂化结构。基态模拟显示，受线性与分叉两种氢键构型竞争影响，氮-氢（N-H）的径向分布函数峰比其它嗪类溶液更宽。平均每个s-四嗪溶质通过氢键吸引约3.0个水分子。在所有采样构型中，57%的构型仅包含线性氢键，而43%的构型则同时包含线性和分叉氢键。激发态模拟同样观察到氢键特征，但平均氢键水分子数降为2.4个，氢键作用被削弱。76%的构型以线性氢键为主，22%包含分叉氢键。基于这些溶液构型计算出的s₀→s₁吸收谱带中心的溶剂位移为1130 cm⁻¹，显著小于其它嗪类化合物（2000-4100 cm⁻¹）；而s₁→s₀荧光位移为560 cm⁻¹，与其它嗪类相当。这一较小的吸收位移结果，是由于s-四嗪的¹(n, π*)激发是高度离域的，其基态和激发态的偶极矩均为零，因此溶剂位移主要由特定的氢键相互作用主导，而非介电溶剂化效应。
研究结论与价值 本研究最终完成了对吡啶、二嗪、三嗪和s-四嗪这一系列稀水溶液中最低¹(n, π*)激发的氢键作用和溶剂位移的系统性调查。该研究证实，生色团在不同电子态下的氢键细节对电子光谱的溶剂位移起着决定性作用。研究揭示，与其它嗪类不同，s-四嗪-水体系中观察到了线性正统氢键和分叉氢键共存的独特现象。在线性构型中，氢键从基态到激发态得以良好保留；而在分叉构型中，其在激发态可能变为过渡态。这些发现反驳了激发态氢键必然断裂的简单化传统观念，凸显了不同嗪类杂环上氮原子数量和位置对氢键网络结构的精细调控。
研究亮点 本研究的亮点主要体现在以下几个方面：第一，研究对象独特：s-四嗪是唯一具有四个氮原子的嗪类，其独特性质填补了该系列研究的空白。第二，发现新颖现象：首次在s-四嗪-水团簇和溶液中明确提出了线性氢键与分叉氢键共存的模型，并发现分叉氢键在激发态转变为过渡态的特殊行为，这丰富了我们对氢键动态学的理解。第三，方法具有系统性：研究将经验性蒙特卡洛模拟、从头算电子结构计算和显式溶剂位移计算方法有机整合，成功将宏观光谱现象（溶剂位移）与微观分子间相互作用（氢键构型、能量、布居数）联系起来，展现了计算化学在解析复杂溶液体系中的强大能力。