最近，一项实验和计算模拟联合研究表明，ReaxFF可以真实地描述功能化的金纳米颗粒，从理论上理解和预测抗生素的释放。这为ReaxFF分子动力学模拟辅助开发可控的纳米粒子原位药物释放带来了希望。 作者使用ReaxFF对共轭壳聚糖-庆大霉素官能化的金纳米粒子在水溶液中的动力学行为进行模拟，研究了基质包埋对抗生素抗菌活性的结构和动力学影响及其释药机理。 模拟体系：约200000个原子（水分子：40000）。图中：柠檬酸盐（棕红色分子）、壳聚糖（浅蓝色分子）和庆大霉素（绿色分子）。 S. Monti, J. Jose, A. Sahajan, N. Kalarikkal, and S. Thomas, Structure and dynamics of gold nanoparticles decorated with chitosan–gentamicin conjugates: ReaxFF molecular dynamics simulations to disclose drug delivery, Phys. Chem. Chem. Phys. (2019)
有机聚合物具有较低的溶解性、较低的自放电速率、较高的机械强度、较高的柔韧性、优异的热稳定性和通用性，有可能成为锂离子电池的电极材料。 北京工业大学孙少瑞课题组计算了文献报道的11种锂离子电池的聚合物电极材料的 ΔEpoly，并与实验值线性拟合得到经验公式： 然后利用半经验公式对未知材料进行预测。考虑到材料的电位、容量和带隙，选择了三种理想的电极材料（16、17、23）。通过对这三种理想电极材料进一步研究发现，LUMO能量越低、杂原子/取代基的电负性越强、材料的电位越高，越有利于设计新型聚合物电极材料。 本文采用AMS软件中的BAND模块完成计算，对于一维周期性体系的计算，由于基组和周期性边界条件的差异，BAND的计算效率与精度，均高于平面波程序。 Huili Lu, Jun Yu, Ge Chen, Shaorui Sun, Theoretical screening of novel electrode materials for lithium–ion batteries from industrial polymers, Ionics (2019)
加氢脱硫工艺是化学工业的基石之一，从石油中将有害的硫脱除，生产出清洁的碳氢化合物。MoS2纳米岛的边缘是催化活性位点，反应压强5-160 bar，温度260-380 ℃。到目前为止，还不清楚这些苛刻的条件如何影响催化剂的结构。因此荷兰莱顿大学Rik V. Mom等人利用高压扫描隧道显微镜，成功地原位观察了该反应条件下活性MoS2模型催化剂。结果表明，MoS2纳米岛活性边缘区的硫、氢和烃覆盖情况随气体环境的变化而变化。通过与密度泛函理论计算结果进行比较，作者提出CH3SH脱硫过程中的优势边结构吸附了硫和CH3SH的混合物。 计算细节： 所有的DFT计算都使用AMS中的BAND模块完成，使用PBE泛函、Grimme van der Waals校正范德华作用，标量方法考虑相对论效应的影响，使用TZP基组描述价轨道，芯轨道被冻结在与自由原子相同的状态，自洽场收敛阈值为10−6 Hatree/bohr，几何优化收敛阈值为10−2 Hatree/nm，其他使用BAND的默认设置。 采用（4×4）MoS2元胞，构造出一维周期性条带。条带包含S型边缘和Mo型边缘。所有计算中S边缘都被完全覆盖（每个Mo边缘原子两个S边缘原子）。元胞长度是在50%S覆盖率优化得到的（1.248 nm）。对于包含金衬底的计算，Au（111）表面使用2层平板模型，晶格常数与MoS2条带相称。选择S-Mo-S-Au-Au堆叠，其中S-Au层间距为0.442 nm。K点总数为3。 过渡态计算：对于给定的反应坐标（通常是原子间距离），先通过线性过渡进行扫描，得到一系列限制性几何结构优化收敛后的能量最高点。对于该能量最高结构，计算与反应有关的原子及其邻近原子的“部分Hessian”。该部分Hessian方程的最负本征值对应过渡态特征，则用来搜索精确的过渡态（势能面鞍点）。一般需要反复多次重复“计算部分hessian”、“过渡态搜索”的过程，最终得到的部分Hessian有且仅有一个负本征值。 Rik V. Mom, Jaap N. Louwen, Joost W. M. Frenken & Irene M. N. Groot, In situ observations of an active MoS2model hydrodesulfurization catalyst, Nature Communications 10, Article number: 2546 (2019)
AMS新增功能 巨正则系综蒙特卡洛（GCMC） 内反应坐标（IRC）扫描 选模振动分析（Mode Selective Vibrational Analysis） 热力学性质（理想气体模型） 简正模PVDOS 分子动力学：分子碎片成分分析 分子动力学：分子枪（分子入射）、分子沉降（模拟过程中移除某些产物） PLUMED：分子动力学分析、多种自由能计算 Collective Variable-driven Hyperdynamics (CVHD) 注意：这些功能，在AMS所有计算模块中均能使用。例如DFTB的分子动力学模拟，也能像传统的ReaxFF分子动力学模拟那样，得到分子数量变化曲线。GCMC模拟在各个模块中均可使用。 GUI 聚合物建模 不同谱图的重叠比对 ADF 新QTAIM与Conceptual DFT AOResponse: 考虑旋轨耦合的极化率、Raman光谱 DRF梯度 BAND DFT-1/2方法得到更好的带隙 SM12溶剂化方法 DFTB GFN1-xTB方法：可以通过ADFView查看分子轨道 改善了SCC收敛性 MOPAC 并行化优化 ReaxFF ADFtrain: 力场优化训练集生成的辅助工具 OpenMP并行技术 COSMO-RS 聚合物热力学性质 新的蒸汽压预测方法
文献资料：Kairui Liu, Guangjin Hou, Jingbo Mao, Zhanwei Xu, Peifang Yan, Huixiang Li, Xinwen Guo, Shi Bai & Z. Conrad Zhang, Genesis of electron deficient Pt1(0) in PDMS-PEG aggregates, Nature Communications volume 10, Article number: 996 (2019) 虽然载体表面稳定单原子的报道已经屡见不鲜，但在液体介质中合成原位还原的稳定弱配位分散金属原子更具挑战性。本文报道了在液体PDMS-PEG中还原H2PtCl6，形成单核缺电子Pt1(0) （Pt1@PDMS-PEG），该还原反应通过紫外-可见光谱、远红外光谱和X-射线光电子能谱得到证实。 CO红外光谱、195Pt与13C核磁共振谱证实Pt1(0)以近似八面体结构(R1OR2)2Pt(0)Cl2H2存在，其中R1、R2 分别是H、C或Si基团。通过比较实验与理论计算的195Pt核磁共振谱，验证了弱配位结构(R1OR2)2Pt(0)Cl2H2以及缺电子Pt1(0) 。质子态H+与Cl−形成类似HCl的弱配位，用碱中和生成铂纳米粒子。Pt1@PDMS-PEG在石蜡氢化硅烷化中表现出超高的活性以及最终加成物选择性。 NMR位移计算使用ADF模块，利用ZORA方法考虑重元素Pt的相对论效应。
文献资料：Stephanie C. C. van der Lubbe, Francesco Zaccaria, Xiaobo Sun, and Célia Fonseca Guerra, Secondary Electrostatic Interaction Model Revised: Prediction Comes Mainly from Measuring Charge Accumulation in Hydrogen-Bonded Monomers, J. Am. Chem. Soc., Just Accepted Manuscript, 10.1021/jacs.8b13358 次静电相互作用（SEI）模型常用于预测和解释自组装体系的相对氢键强度。该模型将氢键的作用机制简化为相互作用的点电荷，但实验结果往往与该模型的预测一致。这一模型是如何预测的？本文研究了两个四重氢键异构体DDAA-AADD和DADA-ADAD。 结果表明，当质子给体D(电子供体)和质子A(电子受体)以DDAA的方式聚在一起时，前沿原子周围的电荷积累比DADA的方式要大。这种电荷积累增强D的正电性以及A的负电性，使得DDAA的静电作用和共价作用都得到增强。因此，该模型之所以有效，是因为它为氢键单体中的电荷积累提供了度量。这为超分子化学家们调控氢键强度、控制自组装体系性能提供了理论依据。 本文使用ADF中的EDA方法分析了氢键键能的构成，静电势分布、前线分子轨道等信息。
文献资料：Patrick M. J. Szell, Scott Zablotny & David L. Bryce, Halogen bonding as a supramolecular dynamics catalyst, Nature Communications volume 10, Article number: 916 (2019) 动力学过程对催化剂、酶、主-客复合物、分子机器等功能分子有重要影响。本文通过氘核磁共振弛豫实验，展示了2,3,5,6-四甲基吡嗪共晶中，卤键对甲基转动的催化作用。作者观察到，与纯的2,3,5,6-四甲基吡嗪共晶相比，卤键共晶中，甲基转动的活化能垒平均下降了56%，氢键共晶下降了36%。 密度泛函计算的结果表明，对位交叉构象不稳定，邻位交叉构象却非常稳定，共同降低了能垒，导致了卤键超强的催化作用。此外，计算结果还表明，卤键的催化能力可能是可调的，卤键供体作用越强，催化作用越强。 本文DFT计算使用AMS中的ADF模块完成，分子模型从晶体结构中提取。因为存在分子间弱相互作用，因此使用Grimme3 BJDAMP色散修正泛函（即-D3(BJ)类泛函），使用ZORA方法考虑相对论效应对重元素的影响。甲基转动过程，使用其中Linear Transit功能完成，每一步转动2.5°。
文献资料：shan-shan zhang, Fahri Alkan, Hai-Feng Su, Christine M. Aikens, Chen-Ho Tung, and Di Sun, [Ag48(C≡CtBu)20(CrO4)7]: An Atomically Precise Silver Nanocluster Co-Protected by Inorganic and Organic Ligands, J. Am. Chem. Soc., 2019 配体的精细选择对合成金属纳米团簇，具有重要意义。使用硫醇、炔基、膦或它们的组合，作为保护金属纳米团簇的配体，是最为广泛的。而无机氧阴离子在这一领域几乎被忽视。本文作者合成了第一个CrO42-/tBuC≡C-共配的Ag48纳米团簇(SD/Ag48，SD=SunDi)，并用单晶X射线衍射(SCXRD)对其结构进行了表征。 SD/Ag48的伪五重对称金属骨架呈核壳结构，由外部Ag25壳包围Ag23圆筒组成。在SD/Ag48表面，前所未见地出现了无机(CrO42-)和有机(TBUC≡C-)配体共存的现象。无机CrO42-阴离子在银纳米团簇的构建中起着重要的作用：(1)钝化Ag23内核；(2)连接核与壳；(3)保护Ag25壳层。该纳米团簇属于14电子超原子体系，从可见区到紫外区呈现连续的类分子吸收带。本工作不仅为银纳米团簇的合成提供了一种新的配体策略，而且也为研究CrO42-对纳米银团簇形状的控制提供了新的思路。 本文使用ADF计算了该结构的电子结构，以及紫外可见吸收光谱。使用ZORA方法中的标量相对论方法考虑相对论效应的影响。
文献资料：Shang-Fu Yuan Cong-Qiao Xu Jun Li Quan-Ming Wang, Ligand‐Protected “Golden Fullerene”: the Dipyridylamido Au32 8+ Nanocluster, Angew.Chem. Int. Ed., 10.1002/anie.201901478 通过吡啶胺和三苯基膦作为保护剂，合成了“金富勒烯”Au32团簇。单晶X射线结构分析表明，该金纳米团簇[Au32(Ph3P)8(dpa)6](SbF 6)2(Hdpa=2,2′-二吡啶胺)，具有S6群对称性，内核为Au328+。量子化学研究阐明了该簇结构的特殊稳定性源于金属-配体相互作用。 量子化学计算使用Amsterdam Modeling Suite(简称AMS)中的ADF模块完成。使用PBE泛函，为了节省计算量，Au元素使用TZP基组、所有其他元素都使用DZ基组，Au元素4f，P元素2P，C、N元素1s及其以内的电子冻结。一般而言，Au元素的相对论效应几乎是所有元素中最强烈的，相对论效应使用ZORA方法中的标量相对论（Scalar选项）。为了简化分子模型，计算使用PH3替代Ph3P。由于ADF中不包含S6群，因此使用其子群Ci群，优化之后得到的结构非常接近S6群。 本文也计算了原子电荷，包括Mulliken电荷布居、Hirshfeld、Voronoi以及多极导出电荷（MDC）。紫外-可见吸收谱使用LB94泛函计算得到，相对其他泛函，LB94具有较好的描述空轨道的能力。
文献资料：Jingxiu Xie, Pasi P. Paalanen, Tom W. van Deelen, Bert M. Weckhuysen, Manuel J. Louwerse & Krijn P. de Jong, Promoted cobalt metal catalysts suitable for the production of lower olefins from natural gas, Nature Communications, volume 10, Article number: 167 (2019) 近年来，由于天然气产量激增，化学原料需求向较轻的碳氢化合物，特别是甲烷转移。合成气(CO和H2)转化的效率，关键取决于催化剂对甲烷和二氧化碳产生的抑制能力。本文作者设计了一种Co/Mn/Na/S催化剂，该催化剂的水-气转换活性可以忽略不计，烃的产物谱也偏离了Anderson-Schulz-Flory分布。在240 °C和1 bar压强下，C2-C4烯烃选择性为54%。10 bar压强下，低碳烯烃和燃料的选择性分别为30%和59%。该催化剂由直径约10 nm的处于hcp金属相Co纳米颗粒组成。作者认为Na与S在Co表面作为电子促进剂，产生协同效应，从而提高了对低碳烯烃与燃料的选择性，同时也大大减少了甲烷和二氧化碳的生成。 文中使用AMS软件中的BAND模块计算了Na2S和Na2O在Co(0001)表面的成键结构。Na2O（图a）与Na2S（图b）在Co表面的结合方式非常相似。只不过O原子在一个亚表面钴原子上方，而S原子在一个空位上方。原包外的原子由较小的球表示，蓝色表示Co，橙色表示Na，黄色表示S，红色表示O。 由于此处研究的是表面吸附，因此使用了色散修正泛函rPBE-D3。其他参数： 基组：TZP与Small冻芯基组 Numerical Quality：Normal k-space：Good SCF收敛阈值：5 × 10−4 Hartree 结构优化梯度阈值：0.001 Hartree/Å 单胞的晶格常数也进行了再优化，得到了a = 2.43 (实验2.51 )和c = 3.91 (实验4.07 )。(0001)面厚度为6原子层，约为12Å ，底部两层原子被冻结，并使用最小基组SZ，大冻芯，Numerical Quality：设置为Basic。表面原包包含4 × 4个原子。为不同区域的原子，设置不同基组，参考中文教程：对同一种元素的不同原子指定不同基组 BAND作为新兴第一性原理计算软件，在材料化学领域具有非常鲜明的特点： 由于BAND中的表面模型为真实的二维模型，因此不再需要人为地添加真空层，也不需要进行偶极矫正，对表面吸附计算的可靠度与效率高于平面波程序。 擅长成键分析，例如表面成键的键能分析（pEDA）、共价键成键机理研究（NOCV），参考案例：硅表面醚吸附－计算表面化学与分子计算化学共同点的最佳案例（Angew. Chem. Int. Ed., 2017） […]
