富勒烯（LuNC@C76）结构优化结果与实际单晶结构相差较大求助

775242627

小白一枚，最近在进行LuCN@C76的结构优化时，计算的碳笼中的LuNC团簇结构与实际单晶测出的团簇角度相差较大，差了有10~20°左右。
泛函：PBE1PBE,B3LYP，C和N原子用的基组分别有6-31G*，def2-SVP和def2-TZVP，金属Lu一直使用的SDD赝势（大核版本），计算用的g16。在上面各种泛函和基组的组合下，最后优化出的碳笼中的LuNC团簇结构与实际单晶测出的团簇角度相差较大。目前有两个问题想请教一下各位大佬：
1、不知道我选择的泛函和基组是否合适？是否有其他泛函和基组推荐？或者换成g09？可能可以优化的结构更接近单晶的结果。
2、如果最后实在优化不好，是否可以直接采用晶体结构进行单点计算，分析LUMO-HOMO和静电势等？如果可以，应该采用什么泛函和基组好一些呢？
晶体结构如图所示，求大神指点迷津，万分感谢。

wzkchem5

可以改用小核赝势，或者做全电子相对论计算试试，毕竟你只有一个Lu，计算量增加不大；此外可以考虑改用范围分离泛函，因为你有大共轭体系，可能有SIE的影响
没有任何理由改用G09
这个分子用晶体结构问题不大，因为不存在容易受晶体堆积力影响的基团，也不存在氢原子等XRD解不出来或解不准的原子。但假如实验晶体有无序，那就要慎用了

775242627

wzkchem5 发表于 2023-5-3 15:34
可以改用小核赝势，或者做全电子相对论计算试试，毕竟你只有一个Lu，计算量增加不大；此外可以考虑改用范围 ...

好的，好的，我试一下，感谢老师的指导和帮助。

七尺贱

这个富勒烯没有别人算过嘛，一个新的富勒烯合成不久之后就会有计算组的人算，如果是早以前合成的，那么找找人家的坐标不就行了