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