硫系化合物半导体长期以来被用作光学信息存储介质，如DVD-RW和DVD-RAM等[1]。近年来，其也被广泛研究用于最有前景的下一代非易失性电子信息存储器，相变存储器(Phase Change Memories, PCM)[2-3]。相变存储器利用相变存储材料在低电导、低反射率的非晶态与高电导、高反射率的晶态之间快速相变的特性来实现信息存储。常用的相变存储材料主要集中于GeTe-Sb2Te3伪二相线上，包括位于两端的母体材料GeTe、Sb2Te3以及三元的Ge1Sb2Te4、Ge1Sb4Te7和Ge2Te2Sb5等[4-6]。其中，GeTe是一种性能优异的相变存储材料[7-8]。相较于Ge2Sb2Te5(GST)，GeTe具有较高的再结晶温度(189℃)和较高的熔点(700℃)[9]，其数据能在105℃左右长期保持，可用于较高温度下的信息存储。然而，随着现代科技的进一步发展，汽车电子和航空航天领域对信息存储提出了更高的挑战。汽车电子领域要求存储器工作温度在120℃左右，而航空航天领域则要求其能在150℃下工作10年。GeTe材料已无法满足这些行业的需求。亟需对GeTe相变存储材料进行性能改进，以满足信息存储领域日益提升的高温应用要求。

元素掺杂是一种有效的半导体改性方式。科研人员对此进行了研究，实验表明Se掺杂GeTe相变存储材料具有更优异的性质。Vinod等制备了不同Se掺杂量的(GeTe)1-xSex薄膜，指出Se掺杂能大幅提高GeTe的再结晶温度，增强其数据保持力[10]；Se掺杂同时能提高GeTe非晶的光学带隙，10%的Se掺杂可大幅提高GeTe的非晶电阻，使晶相与非晶相的电阻差异增加2个数量级[11]。由于GeTe中可以存在高达8%～10%的本征Ge空位[12-13]，Vinod和Sangunni推测掺入的Se原子可能占据Ge空位的位置，从而抑制了GeTe的结晶[11]。然而，这些实验研究主要针对Se掺杂后的物相组成、电学性质、光学性质和热动态性质等。Se掺杂对于体系的原子构型、成键情况等微观性质的影响以及性能增强的作用机制尚不明晰。这些微观性质需要通过第一性原理计算模拟来进一步理解，从而更好地为未来实验提供理论指导。此外，与GeTe同系列的Ge2Te2Sb5的Se掺杂实验表明，当GST1-xSex中Se含量x超过10%时，出现了奇特的相转变性质：掺杂后的GST材料在受热时直接由无定型态结晶为稳态六方，而跳过中间的亚稳面心立方相(face-centered cubic, fcc)[14]。而Wang等进行了Sb2Se掺杂Ge2Sb2Te5的研究，指出当Se原子含量超过9%时，亚稳面心立方相转变为稳定六方相的过程被抑制[15]。Se掺杂PCM对相转变过程的影响较大。对硫系化合物相变存储材料的Se掺杂进行计算研究，也可以为解释相变存储材料晶相演变机制提供线索。

本文采用第一性原理计算方法，对Se掺杂完美GeTe和含本征Ge空位的GeTe的几何构型、成键情况和电子性质进行了研究。结果表明，对于完美GeTe，Se优先取代Te原子。而对于存在本征Ge空位的GeTe晶体，Se掺杂原子不会占据Ge空位，而倾向于取代Ge空位最近邻Te原子。Se原子与Ge空位的相互吸引，形成了较强的Ge—Se键，抑制了GeTe相变存储材料的结晶，从而提高了其数据保持能力和服役温度。Se原子与本征Ge空位最近邻时，使两相的电子局域函数(Electron Localization Function,ELF)分布和体积更加接近，可能是导致奇特相转变性质的重要原因。

1 计算方法

本文计算基于维也纳大学开发的从头算模拟软件包(Vienna Ab-initio Simulation Package，VASP)[16]。体系采用了Ge(4s24p2)、Te(5s25p4)的标准PAW赝势。与先前的GeTe计算一致，本文使用了广义梯度近似GGA-PBE[17-19]来描述交换关联泛函。截断能设置为350 eV，k点设置为3×3×3的Gamma中心k点网格。截断能和k点的设置都进行了收敛性测试。自洽计算的收敛判据设置为电子步能量差小于10-5 eV，而结构弛豫收敛判据则为每个原子上的残余应力均小于0.1 eV/nm。

GeTe包含2个晶相，分别为稳态菱方相(rhombohedral)和高温亚稳面心立方相[20]。信息存储主要通过GeTe非晶态和稳态菱方相之间的快速转变来实现[21]。因此，本文主要研究Se掺杂对于菱方相的影响，并与Se掺杂的面心立方相进行了类比分析。GeTe菱方相是一种扭曲的氯化钠结构，可以看做由GeTe面心立方相中Ge原子亚结构和Te原子亚结构沿[111]方向错动一定距离形成的[22]。这种结构变形属于派尔斯变形(Peierls distortion)[23]。GeTe菱方相中形成了交替的3个Ge—Te长键和3个Ge—Te短键。从[111]方向观察，GeTe也可以被当做层状结构；其中，由Ge—Te短键连接的原子可视为处在同一层。基于GeTe菱方结构的惯用胞，构建了2×2×2的超胞，包含64个原子，结构如图1(a)所示。对GeTe面心立方相，同样采用了2×2×2的超胞，如图1(b)所示。所有构型均进行了结构弛豫。计算优化得到的单胞晶格数据如表1所示。对于无缺陷的菱方相和面心立方相，计算得到晶格常数和晶胞角度分别为a0=0.6086nm，α=88.14°和a0=0.601 9nm，α=90°与实验测量值及前人的计算值吻合[20,24]。由于GeTe中存在大量的本征Ge空位，在掺杂分析中必须考虑Ge空位的影响。对于存在本征Ge空位的GeTe体系的Se掺杂计算，分别构造了Ge31Te32(共63个原子)的菱方结构和面心立方结构，并对所有构型均进行了结构优化。

To ease the calculation, we define angular frequencies and quality factors of primary and secondary sides as:

图1 GeTe晶体结构 Fig.1 Crystalline structure of GeTe

表1 计算得到的单胞晶格数据

Table 1 Calculated crystal lattice datas of single cell

结 构a0/nmα/(°)V/nm3Eg/eV菱方相GeTe0.608688.140.225050.62菱方相Ge32Te31Se10.608387.940.224620.58菱方相Ge31Te31Se10.605887.770.221880.58菱方相Ge31Te320.606287.820.222340.59菱方相Ge31Se1Te320.608488.120.22486—面心立方相GeTe0.6019900.218080.38面心立方相Ge32Te31Se10.6003900.216340.40面心立方相Ge31Te31Se10.6024900.218520.52面心立方相Ge31Te320.5994900.215370.31面心立方相Ge31Se1Te320.6010900.21704—

注：晶格常数a0、晶胞角度α、惯用胞体积V以及带隙Eg。

掺杂形成能的计算式[25-26]为

(1)

本节研究了Se掺杂GeTe完美晶体。当掺杂Se原子时，可能的占位情况有：取代Ge原子、取代Te原子或占据间隙位置。由于Se间隙原子造成较大的晶格畸变，具有较高形成能，故不作比较。计算得到不同取代位置的形成能Ef如表2所示。

本文考虑了缺陷生长环境对于掺杂形成能的影响。当Se取代Te时，在富Te情况下，被取代的Te原子视为与体相Te单质中的一个Te原子结合；而富Ge情况下，被取代的Te原子则视为与体相Ge单质中的一个Ge原子结合，重新形成新的Ge—Te单元。故而在富Te情况下，Te原子化学势即为Te体相单质中单个原子的能量，而Ge原子化学势则为μGe=μGeTe-μTe。其中μGeTe为GeTe晶相中Ge、Te原子对的能量，μTe为Te原子化学势。反之，在富Ge情况下，Ge原子化学势为Ge体相单质中单个原子能量，而Te原子化学势为μTe=μGeTe-μGe。

2 结果与讨论

2.1 GeTe完美晶体的Se掺杂

式中：Efinal和Einitio分别为末态体系和初始体系的总能；μi为第i种原子的化学势，通常是其最稳定单质中单个原子的能量；Ni为缺陷形成后与形成前第i种原子的数量差。

从表2可以看出，对于GeTe的菱方相和面心立方相，无论在富Ge或富Te的情况下，掺杂的Se原子均优先取代Te原子。Se掺杂的形成能为负，表明Se原子容易掺杂到GeTe中。从表1中还可以看出，少量Se掺杂对于GeTe完美晶体菱方结构和面心立方结构的体积影响很小，造成体积的微量坍缩。主要由于Se的原子半径小于Te。此外，较低的形成能和较小的晶格变化，表明在GeTe中能实现高浓度的Se掺杂。这与实验观测一致，实验中Se掺杂量达到20%时，依然没有造成GeTe结构的破坏[11]。

表2 Se掺杂完美GeTe不同取代位置的形成能

Table 2 Formation energy of Se doping at different substitution positions in ideal GeTe eV

Se 掺 杂菱 方 相面心立方相富Ge富Te富Ge富Te取代Ge1.050.880.670.50取代Te-0.41-0.24-0.54-0.37

对GeTe完美晶体的Se掺杂进行了电子局域函数分析，如图2所示。由于Se、Te属于同主族元素，Se原子与Te原子具有十分相似的电子局域函数分布。对于GeTe菱方相，由于派尔斯形变，形成的Ge—Te短键键长RGe1-Te=0.286 nm，Ge—Te长键键长RGe2-Te=0.325 nm，与实验值十分接近[20]。当Se取代Te原子时，也呈现出长短键交替的情况，然而，其中Ge—Se短键变得更短(RGe1-Se=0.266 nm)，而Ge—Se长键变得更长(RGe2-Se=0.334 nm)。这是由于Se原子具有较强的电负性，与其同层最近邻Ge原子(Ge1)形成了较强的Ge—Se键。同时Se原子较小的原子半径，使得其与邻层最近邻Ge原子(Ge2)的键合变弱。Se掺杂使得体系中原子分层排布的现象更加明显。

此外，GeTe面心立方结构则具有6个等同的Ge—Te键，RGe-Te=0.301 nm，如图2(c)所示。Se取代Te原子之后，Se原子的电子局域函数依然呈球状分布，形成6个等同的Ge—Se键，RGe-Se=0.292 nm，键长较Ge—Te键变短。两相中，Se掺杂均形成了较强的Ge—Se键，对Ge原子的移动起到抑制作用。Ge—Se原子主要为共价键结合，详细解释见后续的电子态密度分析。

图2 完美GeTe晶体掺杂前后在 (100) 面上的ELF截面图(等高线间隔为0.14) Fig.2 ELF contour plots on (100) plane for Se doped ideal GeTe crystal (interval is 0.14)

2.2 含本征Ge空位的GeTe晶体的Se掺杂

本文2.1节计算基于GeTe的完美晶体模型。然而，实际上GeTe是一种非化学计量比的半导体，GeTe中存在高浓度的本征Ge空位。对GeTe进行掺杂研究时，必须考虑其中大量存在的本征Ge空位的影响。由于Ge、Se具有相近的离子半径(分别为67 pm和64 pm)，Vinod和Sangunni[11]推测Se原子可能占据Ge空位的位置。因此，本文研究了Se掺杂含Ge空位的GeTe体系(Ge31Te32，共63个原子)。对Se占据Ge空位和Se取代Te原子的情况作了重点比较。

根据掺入的Se原子与Ge空位距离的不同，构建了多种模型。计算得到的形成能如图3所示。当Se占据Ge空位时，Se原子与Ge空位距离即为0。对于菱方相，Ge空位最近邻Te原子有2种，分别为同层最近邻和邻层最近邻，用A1和A2表示。点间连线以便于观察。从图3可以看出，无论对于菱方相还是面心立方相，当存在Ge空位时，掺杂的Se原子依然不会占据Ge空位位置，而是倾向于取代Ge空位最近邻的Te原子。对于菱方相，则优先取代Ge空位同层最近邻的Te原子。从另一个角度来看，Se杂质原子与Ge空位具有一定互相吸引作用。这种吸引作用会抑制Ge空位的移动。

式中时刻的后验估为k 时刻的先验估计为k-1时刻的控制量，wk-1为k-1时刻的过程噪声，系统的控制量

图3 Se掺杂Ge31Te32菱方相和Ge31Te32面心立方相的形成能随Se原子与Ge空位之间距离的变化 Fig.3 Formation energy of Se doping in rhombohedral Ge31Te32 and fcc Ge31Te32 versus distance between doped Se atom and Ge vacancy

图4为菱方相中Se取代Ge空位同层最近邻Te原子时(A1构型)和邻层最近邻Te原子时(A2构型)的原子结构, 虚线环表示Ge空位，深色球、浅色球分别表示Ge、Te原子，1号原子为Se。形成能计算结果表明，Se原子更倾向于占据Ge空位同层最近邻的位置。A1与A2之间形成能差异在50 meV/supercell左右。由图4可以发现，Ge空位使体系中形成了较强的层间Ge—Te键，部分抵消了Ge空位带来的影响。A1构型的成键情况为：3个Ge—Se键(R1-2=R1-3=0.256 nm，R1-4=0.297 nm)，2个层间Ge—Te较强键(R5-6=R7-8=0.312 nm)。A2构型成键情况为：3个Ge—Se键(R1-2=0.282 nm，R1-3=R1-4=0.260 nm)，3个层间Ge—Te键(R5-6=R7-8=0.312 nm，R9-10=0.309 nm)。由于Se取代同层的最近邻Te原子(A1)相较于邻层的最近邻Te原子(A2)所引起的晶格畸变更小，因此更容易形成。

进一步对含Ge空位的菱方相和面心立方相的Se掺杂进行了电子局域函数分析，如图5所示。其中图5(a)～(c)为菱方相，图5(d)～(f)为面心立方相。从图5(a)中可以看出，对于菱方相，Ge空位对ELF的影响十分局域。Ge空位同层最近邻Te原子(Te1),与邻层的Ge原子(Ge1)形成较强的层间Ge—Te结合(RGe1-Te1=0.311 nm)。

图4 Se掺杂GeTe的2种构型 Fig.4 Two configurations of Se doped GeTe

图5 Se掺杂含Ge空位的GeTe体系的在(100)面上的ELF截面图(等高线间隔为0.12) Fig.5 ELF contour plots on (100) plane for Se doped GeTe systems with intrinsic Ge vacancy (interval is 0.12)

图5(b)为Se占据Ge空位同层最近邻时的情况，ELF分布与未掺杂体系十分接近，因此Se掺杂对GeTe菱方相电子分布的作用较小。而当Se占据Ge空位时，Se原子周围的电子分布类似于球形，形成近乎等同的6个Se—Te键，(RSe-Te1=0.309 nm，RSe-Te2=0.316 nm)。此时Se杂质原子对周围Te原子的电子局域函数分布造成较大的影响，导致体系能量升高，从而具有较高的形成能。掺杂的Se原子不会占据Ge空位。

[6] RAOUX S.Phase change materials[J].Annual Review of Materials Research,2009,39(1):25-48.

图7(a)、(b)分别为Ge32Te31Se1菱方相和面心立方相中，掺杂的Se原子与其最近邻Ge原子和Te原子的分波态密度图(PDOS)，虚线表示费米能级位置，设置为0 eV处。可以看出，Ge 4p电子和Te 4p电子具有很强的p-p耦合，形成Ge—Te共价键。而Se原子的PDOS分布与Te原子十分相似。掺杂的Se原子与Ge原子形成类似于Ge—Te键的Ge—Se共价键。Se原子的3p电子态密度分布离价带顶较远，处在较低能量位置。

另外，从表1晶胞数据中可以明显观察到，对于含Ge空位的菱方相，Se掺杂导致晶胞体积减小。而反常的是，Se掺杂含Ge空位的面心立方相导致晶格常数和体积增大。主要是Se掺杂后电子局域函数变化所造成的。Se掺杂后，GeTe菱方相和面心立方相将具有更相近的晶格常数和体积。这对于减小GeTe相变过程中的体积变化十分有益，有助于提高相变存储材料的稳定性。

不同体系的总态密度如图6所示，虚线表示费米能级位置，设置为0 eV处。其带隙变化值如表1所示。Ge空位导致费米能级向价带内移动，呈现出p型半导体的性质。对于GeTe菱方相，Se掺杂对态密度的影响较小，造成带隙的少量降低。3% Se掺杂使菱方结构Ge31Te32的带隙仅仅下降了1.7%。Se掺杂可使GeTe菱方晶体保持良好的导电性质。此外，对于含Ge空位的GeTe面心立方晶体，Se掺杂则会使导带底向能量更高处移动，显著增加其带隙。面心立方结构的Ge31Te31Se1相较于Ge31Te32，带隙增加了67.7%。这主要是由掺杂后晶格常数增大所引起的。

综上所述，采取经皮椎间孔镜治疗腰椎间盘突出症，不仅手术时间短，而且并发症少，可缩短术后康复时间，减轻疼痛，改善功能障碍，临床疗效显著。

图6 菱方相和面心立方相总态密度 Fig.6 Total density of states for rhombohedral phase and fcc

以上分析得出，对于GeTe实际晶体，Se掺杂元素并不会占据Ge空位的位置。掺入的Se原子更倾向于靠近Ge空位，二者具有一定吸引作用，从而限制了Ge空位的移动。GeTe相变存储材料的相转变与电场下Ge空位的定向移动关系密切[27-28]。Se掺杂对于Ge空位移动的抑制作用，是导致其再结晶温度升高的重要原因。此外，实验和计算指出，GeTe和Ge2Sb2Te5中大量存在的本征Ge空位并非以单一Ge空位的形式存在，而是形成多空位群或者小孔洞[29-30]。Se杂质原子可能聚集在空位群或孔洞处，从而对相变存储材料的动态相变过程起到更大的影响。Se掺杂含有Ge空位的GeTe体系，使得两晶相的相似度增加。这种性质也可能是实验中观察到Se掺杂PCM出现奇特相变现象的原因。

图7 菱方相Ge32Te31Se1和面心立方相Ge32Te31Se1Ge原子、Te原子和Se原子的分波态密度 Fig.7 Partial density of states for Ge atom, Te atom and Se atom in rhombohedral Ge32Te31Se1 and fcc Ge32Te31Se1

Se原子与Ge原子的p电子态密度耦合更加显著，因而Ge—Se共价键较Ge—Te共价键更强。

3 结 论

1) 对于GeTe完美晶体，Se掺杂优先取代Te原子。而对于含Ge空位的GeTe晶体，Se原子不会占据Ge空位，而更倾向于取代Ge空位最近邻的Te原子。对于菱方相，则倾向于取代Ge空位同层最近邻的Te原子。Se杂质原子与Ge空位存在一定的吸引作用，阻碍了Ge空位的移动，可以提高GeTe相变存储材料的再结晶温度，使其具有更好的数据保持力和更高的服役温度。

卫星传输技术主要是将新闻采集现场所收集到的视频和音频等内容通过卫星采集系统的处理，发射到与信息同步的卫星上，再通过卫星传送到电视台，该种技术所需要的成本较高，并且操作较为复杂。光纤传播技术主要的信息传输介质为光导纤维，该种技术有着较强的抗干扰能力，并且在传输过程中信息容量较大，但是，这种技术需要铺设光纤，受到了地理环境的影响。但是，4k技术的应用就有效解决了这些问题，其较高的灵活性以及高效性受到了电视台的青睐。

2) Se掺杂导致菱方相的体积收缩，而面心立方相的体积膨胀，减小了由菱方相到面心立方相转变时的体积差异，提高了相变存储材料的稳定性。

3) Se原子掺杂对GeTe菱方相的电子结构影响较小，使其带隙略微减小。故少量Se掺杂不会影响GeTe菱方相的导电性质。此外，Se掺杂使面心立方相的带隙增大。Se掺杂含Ge空位的面心立方晶体呈现类似于菱方相的电子局域函数分布，减小了两相的差异。为解释实验中Se掺杂PCM出现的奇特相转变性质提供了重要线索。

综上，Se掺杂的GeTe是一种性质优异的相变存储材料，具有较高的再结晶温度，适用于航空航天、汽车电子等较高温度的信息存储领域。对Se掺杂后的非晶态性质需要后期进一步探索。

致谢 感谢北京航空航天大学司晨老师在研究过程中给予的帮助和指导。

参考文献 (References)

[1] OVSHINSKY S R.Reversible electrical switching phenomena in disordered structures[J].Physical Review Letters,1968,21(20):1450-1453.

[2] RAOUX S,WELNIC W,IELMINI D.Phase change materials and their application to nonvolatile memories[J].Chemical Reviews,2010,110(1):240-267.

[8] BRUNS G,MERKELBACH P,SCHLOCKERMANN C,et al.Nanosecond switching in GeTe phase change memory cells[J].Applied Physics Letters,2009,95(4):043108.

1°clYf(A)⊆clY{g(A)}。若x∈A,由A为X的可数既约闭集及clcsX{A}={B∈csX:B⊆A},则{x}-∈clcsX{A}。又f=gη及g连续,f(x)=gη(x)=g({x}-)∈g(clcxX{A})⊆clY{g(A)},由x的任意性,f(A)⊆clY{g(A)},于是clYf(A)⊆clY{g(A)}。

[4] LENCER D,SALINGA M,GRABOWSKI B,et al.A map for phase-change materials[J].Nature Materials,2008,7(12):972-977.

[24] SUN Z,TIAN S,SA B.Investigation of the structure and properties of rhombohedral Cu-Ge-Te alloys by ab initio calculations[J].Intermetallics,2013,32(2):292-296.

对于GeTe面心立方相，Ge空位对电子局域函数的影响也较小。与众不同的是，当掺入的Se原子与Ge空位最近邻时，造成电子局域函数的较大变化。此时GeTe体系中的电子局域函数分布甚至接近于菱方相中的分布，部分Ge—Te键出现长短键分化，从而降低了体系的能量，形成能较低。Se掺杂模糊了GeTe面心立方相和菱方相的界线。而Se占据Ge空位时，对周围Te原子的影响较大，形成能较高，不易存在。

[7] PERNIOLA L,SOUSA V,FANTINI A,et al.Electrical behavior of phase-change memory cells based on GeTe[J].IEEE Electron Device Letters,2010,31(5):488-490.

表5反映的是社会保障与基本就业模型的逐个回归，其中城镇化率变量系数从最开始的不显著到最后变得显著，但其系数一直为正，与预期假设不一致；人均收入的系数一直为正，从第（3）个回归开始，其系数在1%的显著性水平上显著，但正负性未发生改变；社会保障财政支出的引入，其变量系数一直为负，符合预期假设，同时在1%的水平上显著；城乡消费差距的系数为正值，随着变量的引入未发生改变，与预期假设一致，但系数并不显著；总抚养比的系数符合前文的假设预期为负值；随着税收的引入，其系数为负，且在10%的水平上显著，与前文假设一致。

[3] WONG H S P,SALAHUDDIN S.Memory leads the way to better computing[J].Nature Nanotechnology,2015,10(3):191-194.

[9] CHONG T,SHI L,ZHAO R,et al.Phase change random access memory cell with superlattice-like structure[J].Applied Physics Letters,2006,88(12):122114.

[10] VINOD E M,SINGH A K,GANESAN R,et al.Effect of selenium addition on the GeTe phase change memory alloys[J].Journal of Alloys and Compounds,2012,537(19):127-132.

[11] VINOD E M,SANGUNNI K S.The effect of Se doping on spectroscopic and electrical properties of GeTe[J].Thin Solid Films,2014,550(1):569-574.

作者的父亲在上世纪90年代中期开始，就饱受阿尔茨海默氏病的折磨，病情逐渐加重后，子女们不能不请看护日夜照料他，并请假回老家陪伴他；连与父亲离婚多年的母亲，基于同情与责任，也回到老宅，替父亲整理堆积了几十年收藏与回忆的家，为他理发、剪指甲。

[12] KOLOBOV A V,TOMINAGA J,FONS P,et al.Local structure of crystallized GeTe films[J].Applied Physics Letters,2003,82(3):382-384.

两年历练，许峰比李莉现实，他率先戳破了童话，李莉再没有过多追问，许峰坦白了自己的欺骗后再没有多作解释。他们都明白，他们之间，这样是最好的结局。

[13] TONG F,MIAO X S,WU Y,et al.Effective method to identify the vacancies in crystalline GeTe[J].Applied Physics Letters,2010,97(26):261904.

[14] VINOD E M,RAMESH K,SANGUNNI K S.Structural transition and enhanced phase transition properties of Se doped Ge2Sb2Te5 alloys[J].Scientific Reports,2015,5:8050.

[15] WANG M,LU Y,SHEN X,et al.Effect of Sb2Se on phase change characteristics of Ge2Sb2Te5[J].CrystEngComm,2015,17(26):4871-4876.

[16] HAFNER J.Ab-initio simulations of materials using VASP:Density-functional theory and beyond[J].Journal of Computational Chemistry,2008,29(13):2044-2078.

[17] BLÖCHL P E.Projector augmented-wave method[J].Physical Review B,1994,50(24):17953-17979.

[18] PERDEW J P,WANG Y.Pair-distribution function and its coupling-constant average for the spin-polarized electron gas[J].Physical Review B,1992,46(20):12947-12954.

[19] MIAO N,SA B,ZHOU J,et al.Investigation on Ge5-xSbxTe5 phase-change materials by first-principles method[J].Applied Physics A,2010,99(4):961-964.

3.选聘信息平台建设可以促进职业经理人市场的发育和完善。职业经理人选聘信息平台建设是我国市场经济发展的必然结果，也是职业经理人市场逐渐走向成熟的标志。信息平台的建设有利于企业选聘职业经理人，职业经理人的行为也将被规范，我们的职业经理人市场将逐步走向正规，市场也会逐渐成熟和完善。

[20] GOLDAK J,BARRETT C S,INNES D,et al.Structure of alpha GeTe[J].The Journal of Chemical Physics,1966,44(9):3323-3325.

[21] LEVIN E M,BESSER M F,HANUS R.Electronic and thermal transport in GeTe:A versatile base for thermoelectric materials[J].Journal of Applied Physics,2013,114(8):083713.

棚架绿化主要是在构件不同的棚架下端种植攀援类植物，棚架绿化的植物布置与棚架的功能和结构密切相关。通常情况下，棚架功能分成两种，即经济性和观赏性，在对经济性植物进行种植时，可以选取葫芦、茑萝等植物，在对观赏性植物进行种植时，则需要选取一些开花观叶植物，进一步提高观赏效果[3]。在结构上，不同的结构选取的植物种类也有较大的差异。如果棚架是由砖石或者是混凝土结构构成的，则可以种植紫藤、凌霄等植物，如果棚架是由竹、绳等结构制成，则宜选用牵牛花、啤酒花等攀援植物[4]。

[22] CHATTOPADHYAY T,BOUCHERLE J.Neutron diffraction study on the structural phase transition in GeTe[J].Journal of Physics C:Solid State Physics,1987,20(10):1431.

[23] PEIERLS R E.Quantum theory of solids[M].Oxford：Oxford University Press,1955.

[5] LENCER D,SALINGA M,WUTTIG M.Design rules for phase-change materials in data storage applications[J].Advanced Materials,2011,23(18):2030-2058.

[25] ZHANG S,WEI S H,ZUNGER A,et al.Defect physics of the CuInSe2 chalcopyrite semiconductor[J].Physical Review B,1998,57(16):9642.

SAP患者在病情发展过程中由于大量呕吐及渗出，会导致机体有效循环血容量不足，从而损伤肠道屏障功能。肠道屏障功能受损后会加重机体肠道动力障碍，最终出现肠道细菌移位，内毒素入血，导致多器官功能障碍的发生［9］。近些年，液体复苏被认为是SAP急性期治疗的基础，充分、早期的液体复苏可改善患者的血流动力学，降低病死率［3］。因此，在SAP的标准治疗中，液体复苏已经得到广泛的认可，并成为治疗方案中重要组成部分。但是液体复苏的时机、复苏液体的种类、复苏液体量和速度以及复苏终点等都是目前临床治疗方案的研究点。因此，探求有效的液体复苏方式和改善胃肠功能障碍的治疗对SAP具有重要的临床意义。

[26] WEI S H.Overcoming the doping bottleneck in semiconductors[J].Computational Materials Science,2004,30(3):337-348.

[27] NAM S W,CHUNG H S,LO Y C,et al.Electrical wind force-driven and dislocation-templated amorphization in phase-change nanowires[J].Science,2012,336(6088):1561-1566.

[28] NUKALA P,AGARWAL R,QIAN X,et al.Direct observation of metal-insulator transition in single-crystalline germanium telluride nanowire memory devices prior to amorphization[J].Nano Letters,2014,14(4):2201-2209.

Belaruskali公司总经理I.Golovaty认为：“新合作备忘录的签署，昭示着我们的产品在中国拥有巨大的市场和广阔的发展前景。我们相信，通过对产品的进一步升级，添加锌和硫元素的肥料产品将会走出亚洲，受到世界各国的青睐。此外，中化集团是世界上最大的钾肥进口商之一，与之合作，不仅能够使双方互利互惠，实现共赢，还能够促进中国与白俄罗斯之间的经济贸易交流，为两国谋求共同发展注入新动能。”

[29] SUN Z,ZHOU J,BLOMQVIST A,et al.Formation of large voids in the amorphous phase-change memory Ge2Sb2Te5 alloy[J].Physical Review Letters,2009,102:075504.

ISA95标准定义了企业商业系统和控制系统之间的集成，主要可以分成三个层次，即企业功能部分，信息流部分和控制功能部分[5]。企业功能基于Pursue大学当初建立的CIM功能模型；信息流部分基于Pursue大学的数据流模型图和S88批次标准，包括产品定义、生产能力、生产计划和生产性能4种信息流；而其控制功能则基于Pursue和MESA的功能模型[6]。

[30] UPADHYAY M,ABHAYA S,MURUGAVEL S,et al.Experimental evidence for presence of voids in phase change memory material[J].RSC Advances,2014,4(8):3659-3668.