先按你的四段格式要求整理：标题翻译、关键词、摘要专业翻译、再给出基于正文的学术解读。中文标题：芳基偶氮吡唑聚集体中的分子激子：一项量子化学研究

《Scientific Reports》：Molecular excitons in arylazopyrazole aggregates: a quantum chemical study

　　

分子光致开关的聚集可能影响其功能性。从基础层面看，单体在聚集态中的相互作用会导致激子态（exciton states）的形成，而这些激子态进一步支配由此类光致开关构成材料中的能量转移与电荷转移过程。本研究利用量子化学计算，考察了芳基偶氮吡唑（arylazopyrazole）聚集体的激子态。芳基偶氮吡唑是过去十年中作为偶氮苯（azobenzene）替代体系而受到广泛关注的一类光致开关。研究人员对包含多达32个芳基偶氮吡唑单体的聚集体开展了簇模型激发态计算，同时还对其晶体结构实施了周期性计算。研究获得并分析了激子态的组成、激子分裂以及从单体到聚集体的光谱位移，从而为真实芳基偶氮吡唑聚集体的电子态与吸收光谱提供了定量认识。

该论文发表于《Scientific Reports》，围绕芳基偶氮吡唑（arylazopyrazole, AAP）聚集态中的分子激子行为展开，核心问题是：当光致开关分子在晶体、单分子层、自组装体系或其他超分子环境中彼此靠近时，分子间相互作用会诱导形成集体电子激发态，即分子激子（molecular exciton）。这种聚集诱导的激子结构会直接影响材料的吸收光谱、能量转移与电荷转移过程，因此也是理解聚集态光响应功能的关键。既往研究已经对偶氮苯（azobenzene）聚集体的激子性质进行了较充分探索，但对于结构相近、近年来迅速兴起的芳基偶氮吡唑体系，尤其是在真实晶体构型中的激子形成规律、光谱位移模式以及激子局域/离域特征，相关认识仍然有限。芳基偶氮吡唑因具有优良的光稳态组成与较长的顺式异构体热半衰期而成为重要光开关候选，因此有必要从电子结构层面对其聚集态行为进行系统研究。

研究人员以甲氧基取代的methoxy-AAP1晶体实验结构为基础，构建了从二聚体到三十二聚体的一系列模型，包括沿不同晶轴方向截取的二聚体、四聚体、一维线性八聚体、二维十六聚体以及三维晶体簇模型，并结合周期性晶体计算，系统分析了不同聚集方向、聚集维度与分子排布方式对吸收带位置、激子分裂大小以及激子组成的影响。研究表明，芳基偶氮吡唑聚集体中的最低激发态主要来源于单体的 \(S₁\)（\(n\pi^*\)）与 \(S₂\)（\(\pi\pi^*\)）态，且总体上以局域激发（local excitation, LE）为主，电荷转移（charge transfer, CT）成分较小，因而其本质主要属于Frenkel激子。与此同时，不同堆积模式可分别诱导H-聚集（导致蓝移）或J-聚集（导致红移）特征，真实晶体中则常表现为二者竞争与耦合后的综合结果。该研究因而为理解AAP聚集体和晶体的电子光谱本质提供了量化依据，也为今后设计基于此类光开关的响应材料提供了理论支撑。

在方法上，研究人员首先对AAP1与methoxy-AAP1单体进行了密度泛函理论（DFT）几何优化与频率分析；随后采用长程校正含时密度泛函理论（time-dependent density functional theory, TD-DFT），以 \(\omega\)B97X-D 为主要方法，对从 \(N=2\) 至 \(N=32\) 的簇模型实施激发态计算，并用ADC(2)及其自旋标度变体对单体和二聚体结果进行对照评估。研究还基于实验晶胞开展了VASP周期性DFT/TD-DFT计算，并借助跃迁密度矩阵分数（fraction of transition density matrix, FTDM）与自然跃迁轨道（natural transition orbital, NTO）分析激子组成、局域化程度及CT贡献。样本来源为Nagai等人报道的methoxy-AAP1实验晶体结构。

在单体结果部分，研究首先比较了母体AAP1与甲氧基取代衍生物methoxy-AAP1的吸收性质。结果显示，甲氧基取代使 \(\pi\pi^*\) 吸收带发生明显红移，说明给电子取代基改变了电子分布并稳定了相关激发态。研究人员进一步发现，双\(\zeta\)与三\(\zeta\)基组会影响绝对峰位，但AAP1与methoxy-AAP1之间的相对红移幅度基本稳定。对于methoxy-AAP1，几何构型同样影响吸收带位置：采用实验几何时，光谱相对于优化构型表现出一定蓝移，说明分子非平面性会改变激发态能量。NTO分析证明，主要吸收带对应 \(S₀\rightarrow S₂\) 的 \(\pi\pi^*\) 跃迁，而 \(S₀\rightarrow S₁\) 为暗态 \(n\pi^*\) 跃迁；在methoxy-AAP1中，氧原子参与相关NTO分布，与其红移现象相一致。进一步的方法对比表明，CAM-B3LYP与 \(\omega\)B97X-D 结果接近，B3LYP偏红移，CIS则明显蓝移，而ADC(2)对主吸收峰的预测与 \(\omega\)B97X-D 高度接近，支持后者作为后续大体系“主力”方法的合理性。

在“Dimers”部分，研究人员考察了沿晶体A、B、C方向提取的三个最近邻二聚体。结果显示，不同相对取向导致截然不同的激子耦合与光谱位移模式。dimer A相对于单体几乎无位移，说明其激子耦合与范德华位移效应近乎相互抵消；dimer B出现明显蓝移，符合H-聚集特征；dimer C则表现出轻微红移，属于弱J-聚集特征。无论采用TD-\(\omega\)B97X-D还是ADC(2)，二聚体激子分裂大小规律均一致，即dimer B最大，dimer A次之，dimer C最小。研究进一步指出，尽管dimer A总光谱位移接近零，但其激子分裂并不小，这说明“峰位几乎不动”并不等同于“分子间耦合很弱”，而是多种效应抵消的结果。FTDM分析显示，二聚体中的 \(n\pi^*\) 态较为局域，\(\pi\pi^*\) 态则更为离域，但整体仍以LE为主；除dimer B的 \(S₅\) 呈显著CT态外，其余低激发态的CT贡献均很小。这一结果说明AAP二聚体中的主要光谱行为仍由Frenkel型激子主导，而非电荷转移激发主导。

在“Larger aggregates”部分，研究将体系扩展到四聚体、八聚体、十六聚体和三十二聚体，以揭示聚集维度和空间方向对激子性质的系统影响。对于四聚体，tetramer AB和晶胞模型表现为蓝移，而tetramer AC则略微红移，这与其内部分别包含较强H型排列或较多J型排列有关。FTDM结果表明，\(\pi\pi^*\) 态比 \(n\pi^*\) 态更加离域，但其对角元分布并不对称，显示真实晶体构型中的分子取向会破坏理想化模型中“每个单体等权贡献”的情形。

对于 \(N=8\) 体系，线性八聚体的蓝移最为显著，是全文所有模型中最大蓝移者，表明沿特定方向的连续H型耦合能够最大程度抬升明亮 \(\pi\pi^*\) 激子态能量。相比之下，同样含8个分子的三维crystal-211模型仅呈中等蓝移，说明三维堆积中不同方向相互作用的竞争削弱了纯H聚集效应。FTDM分析显示，线性八聚体最亮态离域程度较高，而crystal-211中的 \(n\pi^*\) 态则更局域。

对于 \(N=16\) 体系，二维hexadecamer与三维crystal-221均表现出较大的激子分裂，但其峰位与 \(N=8\) 体系相比并未持续单向蓝移。研究人员指出，沿C方向增加的第二层分子引入了J-聚集特征，这会抵消原有H-聚集导致的蓝移，因此hexadecamer虽仍保持蓝移，但比线性八聚体有所回落。crystal-221同样同时存在AB平面中的H-聚集和C方向上的J-聚集，因此其最亮态排序和光谱响应都体现出混合聚集的结果。FTDM矩阵继续表明，多聚体中的激子总体仍是局域激发构成，CT贡献很小。

在最大簇模型crystal-222（\(N=32\)）中，研究得到全文最大的 \(\pi\pi^*\) 激子分裂0.63 eV，说明随着体系尺寸增大，激子带宽明显扩展。尽管如此，其单体到聚集体的光谱蓝移仅为0.09 eV，并未随尺寸无限增大而持续增强。研究人员据此指出，在真实三维晶体中，H-聚集与J-聚集效应持续并存，最终导致适中的净光谱位移而非极端蓝移。FTDM空间映射进一步显示，最亮 \(\pi\pi^*\) 激子态主要集中在聚集体中部区域，表明在有限尺寸三维簇模型中，边界分子对激子态贡献较小，而内部区域承担了更主要的离域激子成分。对一系列三维晶体模型的比较还表明，吸收强度随分子数增加而增强，但峰位蓝移在 \(N=16\) 达到较高值后，在 \(N=32\) 时反而下降，进一步支持“聚集尺度增大并不简单对应更强蓝移”的结论。

在“Periodic calculations”部分，研究人员又采用周期性DFT框架直接计算体相晶体吸收性质，以检验簇模型结论的稳健性。研究比较了Green–Kubo（GK）、Tamm–Dancoff近似（TDA）和完整TD-DFT三种周期性方案，并与Gaussian中的簇模型TD-PBE结果对照。结果显示，不同计算框架对绝对峰位预测差异较大，但周期性TD-DFT给出的单体到晶体蓝移约为0.11 eV，与簇模型TD-PBE结果约0.12 eV高度一致，也与实验中单体溶液到晶体薄膜约0.07 eV的蓝移趋势相符。相比之下，GK与TDA方法得到的峰位位移接近于零，且对能带截断较敏感，因此研究人员认为周期性TD-DFT结果更具稳健性。相关能带与分波电荷密度分析还表明，价带顶部与 \(\pi\) 轨道、n轨道有关，导带底部则具有典型 \(\pi^*\) 分布，这与文中关于 \(n\pi^*\) 和 \(\pi\pi^*\) 激发来源的解释相互印证。

讨论与总结部分强调，本研究系统建立了芳基偶氮吡唑聚集体中激子形成、光谱位移和空间局域化之间的联系。首先，AAP聚集体低激发态以LE主导，CT成分通常很小，因此采用Frenkel激子框架理解其主吸收行为是合理的。其次，\(n\pi^*\) 激子较 \(\pi\pi^*\) 激子更局域，而 \(\pi\pi^*\) 激子在较大聚集体中展现出更明显的离域扩展。再次，光谱位移方向并不由聚集本身单独决定，而是依赖具体堆积模式：H型排列倾向于蓝移，J型排列倾向于红移，真实晶体常体现二者竞争后的净结果。最后，大尺寸簇模型与周期性计算均支持methoxy-AAP1从单体到晶体存在约0.1 eV量级的蓝移，这一理论结果与实验薄膜光谱趋势相近。

研究结论可概括翻译如下：研究人员基于第一性原理量子化学计算，系统研究了由methoxy-AAP1构成的芳基偶氮吡唑聚集体的吸收光谱与激子态。聚集体结构取自实验测定晶体构型，主要采用长程校正TD-DFT计算了从二聚体到三十二聚体的激发态，并对二聚体辅以ADC(2)系列方法验证，同时对晶体结构进行了周期性DFT与TD-DFT计算。FTDM分析表明，这些聚集体中由单体 \(S₁\)（\(n\pi^*\)）与 \(S₂\)（\(\pi\pi^*\)）态演化而来的低激发态主要由局域激发构成，CT贡献较小，因此本质上属于Frenkel激子；其中 \(n\pi^*\) 激子比 \(\pi\pi^*\) 激子更局域。所有模型中最大的 \(\pi\pi^*\) 激子分裂为0.63 eV，出现在crystal-222模型；二聚体中最大分裂为0.24 eV，出现在H型的dimer B。不同聚集方式可导致蓝移或红移，其中最大蓝移0.23 eV出现在线性八聚体，最大红移仅?0.03 eV出现在以J型特征为主的tetramer AC。周期性TD-DFT估计单体到晶体的位移约为0.1 eV，接近实验单体到薄膜的约0.07 eV位移。总体而言，该研究为芳基偶氮吡唑聚集体与晶体的电子光谱性质提供了新的定量认识。