位置加权的改进粒子群算法

针对基本粒子群算法具有后期收敛速度慢、容易陷入局部极值等缺陷,通过考虑粒子的位置之间的加权作用,对基本粒子群算法进行了改进,提出了一种位置加权的粒子群算法以减小搜索过程中的盲目性。测试函数结果表明,算法的收敛性以及收敛速度与粒子群算法位置加权因子有很大关系,通过选择合适的加权因子能有效提高算法的计算效率,算法适用于地球物理优化领域的波动方程反问题。

通过第一性原理预测内地核的成分与结构

内地核成分与结构的确定一直是地球深部研究的重要课题。目前地核的公认成分是铁和少量的镍。但由于地核密度低于纯粹的铁镍合金(固态内核2%～3%,液态外核6%～7%),其中必定掺杂有一定量的轻元素,其种类与浓度有待确定。除成分外,地核条件下铁的晶体结构也存在争议。根据地震学观测,声波沿地轴方向的传播速度比赤道平面方向快大约3%～4%。这意味着内地核是各向异性的;但在极端高压下,晶体结构中的原子应该按致密的密排六方结构(h.c.p)排列,而h.c.p结构对声波传输是高度各向同性的,这就需要确定地核条件下铁的晶体结构。根据第一性原理计算得到的高压下体系能量以及爱因斯坦谐振子模型,本项研究估算了给定结构的自由能以及掺杂轻元素后的影响。根据计算结果可以定性的分析得出,在高压OK下致密的h.c.p结构显然比疏松的体心立方(b.c.c)更稳定;而随着温度的升高,原子核的振动造成b.c.c结构的自由能比h.c.p结构下降得更快,因此在高温下b.c.c结构更稳定;掺杂轻元素后,这种优势变得更加明显,而3.6at.%的Si则恰好同时解释了2%～3%的密度缺失和b.c.c结构在内地核条件下的稳定性。因此我们建议内地核的基本结构与成分应为以体心立方结构存在的铁,掺杂约3.6at.%的硅元素,内地核温度至少在5500K以上。这一结论与其它更复杂的方法得到的结果一致。