“不是，徐佑，你主攻的方向不是物理嘛，什么时候又研究起算法了？而且，材料模拟是很难完全彷真的，不是你随便编几个程序，就能做出来什么成果的。”
“徐教授，除了物理之外，数学、计算机等学科，也同样是我主攻的科目。不管怎样，先让我尝试一下吧。”
“行，数据在我电脑里，你去看就是，密码是你妈的生日。”
“好的，我知道了。”
徐佑并不知道，自己现在的计算机水平，到底到达了一个什么样的程度。
但至少，计算机专业的本科毕设课题，对于徐佑来说，已经可以轻松的解决了。
徐佑觉得，掌握好了材料学与信息学的基础，并能很好的将材料与模型联系到一起的话。
还是有希望通过模拟去解决材料结构的问题的。
徐佑顺利的解开了徐洲电脑的密码，并找到了相关的资料文件。
通过徐洲电脑里的各种信息，徐佑了解到了，徐佑在金属原子簇制备中的方法和数据。
铝原子簇和铁原子簇的制备，虽然总体上比较顺利，但还是有一些原子数量的原子簇，是没有成功制备出来的。
至于其他的金属原子簇，制备起来更是完全没有头绪。
这说明仅仅是同一种结构模型，还是完全不够的。
“想要得到一种原子簇的最稳定结构，需要保证原子簇的势能是最低时才行。”
很快，徐佑便找到了其中的关键。
当原子簇的势能比较高时，整体结构显得不稳定，原子簇很容易被破坏。
而找到势能与空间结构的关系，则成为了问题的重点。
“看样子，这是一个数学问题啊。”
影响到原子簇势能的因素有很多，而具体每个因素的影响方式，也是各不相同的。
想要构造一个精准的，能够表征原子簇的空间结构与其势能之间关系的函数，并不是一件容易的事情。
“先把几个影响因素罗列出来吧。”徐佑想道。
首先，最重要的一个因素，就是原子的个数。
整体的趋势上，原子的个数越多，原子簇的势能越大。
这也是为什么，当原子簇数量足够多时，原子簇材料非常难以制备的原因。
其次的因素，便是原子簇的结构方式。
徐佑通过表示原子簇中，各个原子的空间三维坐标，去描述原子簇的结构方式。
除了这些之外，徐佑也将其他的次要影响因素也罗列了出来。
通过徐洲电脑中的各项数据，徐佑进行着函数参数的调整和修正。
徐洲在闲暇时，看到徐佑认真的样子，心里不禁对徐佑夸赞起来。
“不愧是我徐洲的儿子啊，不仅天资过人，学习态度也非常认真。可惜时间有点短，要不然他还真说不定做出什么成果出来。”
这个时候，徐洲自然是不会对徐佑抱有什么期望，也只是把徐佑当成一个实习的本科生去看待。
徐佑在写满了一篇又一篇的草稿纸之后，终于得到了一个满意的函数关系式。
“这样一来，原子簇的势能就可以表示出来了。”
徐佑已经对之前的数据，进行了多次的验证，确保这个函数公式得出的结果，与真实值的误差是非常小的。
而通过这个函数公式，也解释了为什么个别原子个数的原子簇，是难以通过同样的结构去制备出来的。
完成了函数公式的设计后，徐佑带着公式找到了徐洲。
“嗯？这个是什么玩意？”
看着徐佑纸上的公式，徐洲竟有些摸不着头脑。
其实材料这门学科，与数学、物理、化学都是有很大关系的。
在平时的科研中，也会用到很多数学公式。
但徐洲还从来没见过这样的一条公式。
“徐教授，这是金属原子簇势能的影响公式。这里面，E表示原子间的作用势能，N表示原子个数，ε表示二聚体的势能阱深度，σ表示二聚体的核心距，rij为原子簇中任意两个原子i与原子j间的欧氏距离。”
徐佑给徐洲讲解着，公式中各个物理量的含义。
“虽然看不懂这个函数是怎么构造出来的，但确实很像那么回事儿啊。”徐洲心说道。
徐洲心里清楚，这些物理量，确实是影响势能的因素。
而徐佑构造出来的这个公式，也很有那种教材上标准公式的感觉。
即使徐洲还不能确定公式的正确与否，但却非常想和去验证一下。
“那我们挑几个数据，来验证一下吧。”徐洲说道。
“不用了，徐教授，我已经验证过了。至少目前的这些数据，都是没有问题的。”
听到这话，徐洲却有些不敢相信。
“真的都试过了？”
“不信的话，你再试一下吧。”
徐洲也不磨叽，连忙把自己的一些数据，代入验证了一番。
而当看到这一个个结果后。
徐洲也不得不承认，通过这个公式去计算，所得出来的结果确实是正确的。
即使个别的数据，有极小的误差，但也完全在允许的范围之内。
“可以啊，儿子！你到底是怎么想出这样一个公式的？”
“徐教授，在实验室里，我们尽量还是以师生相称。”
“说得也是。那徐佑同学，说说你的思路吧。”
徐洲也明白，在这样的氛围下，尽量还是忽略父子的身份更好一些。
随后，徐佑向徐洲讲述着，自己构造的一些思路。
徐洲自然还是能听懂大部分的，只是想要独立的创建出这样一个函数公式，远不是能听懂这么简单。
“怪不得有几个铝原子簇和铁原子簇是做不出来的。按照这个函数公式，这几个原子个数的原子簇，会出现势能反常的现象。”
虽说整体上，原子个数越多，势能越大。
但这个规律也不是绝对的。
解决了函数公式的问题后，接下来的问题仍然很棘手。
即使知道什么时候势能大，什么时候势能小。
但想要得到最小的势能，并不是通过简单的计算就能够确定的。
在这里面，有太多的影响因素。
徐佑还需要构造一个合适的算法，去计算出每个原子数对应的最小势能结构。</p>