â³µ¢j[KXuÀKXÌªq®Íw@åw¶YZp¤ä ã V

Molecular Dynamics Calculation for glasses
Hiroyuki Inoue
Institute of Industrial Science, The University of TokyoA.

äÌ¤ÌêðÐîµ½¢BÍ¶ßÉªq®Íw@ÍC1950 NãÌã¼É¶Ýo³êC1970 NãÉKXóÔÉßÄKp³B. ªq®Íw@ê½1)BùÉ»±ÅÍCVJKX\¢Ìîªq®Íw@ÍC½ÌnÌ^®ûö®ðl{Æà¾¦é SiO4 lÊÌªvZ@ãÅÄ»³ðÍIÉð«C±qÌOÕâzuðvIÉ

êÄ¢éBvZ@Ì\Íâè@Ì­WÆ¤ÉC

·é±ÆÉæèC»Ì±qzñÉæé¨«llXÈKXg¬É¨¯é´qzñÌÄ»â»â\¢îñð^¦évZ@V~Ì¨«ÖÌ¢ðSÛèÆµÄCKXÞ¿@ÅÌªìÉ¨¢Äàªq®Íw@Íp¢çêÄ«â³ÍjðÝèµC±qzuðîÉ^®ûö®

½B·ÅÉC{Ì Vol. 13 No. 1 Ìuj[VèéB±qÔÌÝìpÆnÌóÔi·x[Kðð­±ÆÉæèC±qÌzñâ±qÌUXÌvZÈwIAv[`vÌÁWÌÉZ®â^®ðlðÆµÄ¾é±ÆªÅ«éBµÌÌV~éªC»Ì¦é±qÌÍCÅåÅ 106 ÂöxÅC»ÌKXóÔÆÌÎäªÐî³êÄ¢éBNa2O­

Ôà 10|9 s öxÆÀçêÄ¢éB} 1 É¦µ[VªSÅÍ

B2O3 nÉ¨¯é_f 4 zÊÌ BO4 PÊÌg¬½æ¤ÉC

ÉÎ·éÏ»EñË´_fÌ¶¬CN2O­SiO2vZµC»ÌÍ Fi Æ±qÌÂ¬x©çCDté t Å±q i É­­Í Fi ðnÉ¨¯é Qn ÌªzCNa­B­O­Cl nâ Na­Si­O­N nÉ¨¯éªâ\¢PÊC¬AJøÊÈÇCåÉ 2 ¬ª¬nÉ¨¯é´q\¢â¬ÌÀè«ÉÂ¢ÄCÀÛðñíÉæ­Ä»Å«é±Æª¦³êÄ¢éB³çÉCÚ×Èè@ÍCQlÐÈÇ2)ðQÆ¢½¾­
ÆµÄCÈPÉî{IÈè@ðà¾µ½ãCä8505 sÚæîê 4_6_1§153_TEL 03_5452_6315
FAX 03_5452_6316
E­mail: inoueiis.u­tokyo.ac.jp

Fig. 1

The trajectories of particle i in the molecu­
lar dynamics calculation. r, [and F are thecoordinate, velocity and force of particle i,respectively.
63

NEW GLASS Vol. 19 No. 4 2004

bãÌ±q i ÌÊuÀWÆ^®GlM[ðð\µCr ÍÌd³ð\·p[^ÅßéB±êðCÎÛÆ·éSÄÌ±qÉÎµÄéBC ¨æÑ D Í»ê¼êoÉq_oÉqCos¢CÀWðXVµCÄÑVµ¢zuðp¢Éq_ldÉqÌp[^ÅÄCÍðvZ·éB±ÌvZÌJèÔµðs¤ÉÎµÄCp[^ZbgðßÄClXÈ±ÆÉæèCÝèðºÅÌÀèÈ±qÌzu¨¿ÉKpÅ«êÎæ¢ªCÀÛÉÍCÂXÌéBeCIâC»ÌóÔÉ¨¯é±qÌÚ®âU®ðménÉæÁÄCp[^ð²®·éìÆâV½±ÆªÅ«éBp¢çêéÊuÀWÌXVÌû

È|eVÖÌtÁªKvÅ@ÉÍCx@EMA@ªéBZÌâKçêép[^ÍC®aªzÖÌs[NÊXÈÇÌêCµ¦é±qªÀÝÌXP[uðÄ»·éæ¤ÉsöëÉæÁÄßéûéBp¢ÆärµÄCÉßÄ­È¢½ßÉC»Ìæ¤@C»\¢Ì´qÊuÅ|eVªÉ¬

Èî{ZÆÄÎêéÌÉ´qâCIðÆÈéæ¤Éèßéû@CªqO¹@ÉæÁÄzèµ½±qðzuµCî{ZÆ¯¶±qz¾çêéNX^Ì|eV©çp[uÌvJðüÍÉzu·é±ÆiüúI«^ðÝè·éû@ÖÆiñÅ«½BÜ½CLøEðjÉæèC­È¢î{ZàÌ±qÅ

d×CMorse |eVC3 Ì|eVà\ÊÌe¿ÌÈ¢óÔðìèo·±ÆªÅ«

ÈÇÌ±üàsíêÄ«½B³çÉCæê´éBzñÌüú«Ìe¿ð¬³­³¦é½ßªq®Íw@àp¢½vZàñ³êÄ¢é

ÉCî{ZÍÂ\ÈÍÍÅå«­µC½­ÌªC»ÝÌÆ±ëvZ×ª¢½ßC2 Ì|±qðzu·é±Æª]ÜêéBâNXeVC^ÈÇðÎÛÆ·éêÉÍC±Ìüú«ð 2ðp¢Äì¬µ½´qzuÌÉaÉp¢é³ÉÏXCáª½¢BªqO¹@ª»ÝÌ ab initio @Öé¢ÍCp¢È¢±ÆÉæèK

p·é±ÆªÅ«éBÜ½Cé¢ÍC2 Ì{3 Ì|eVé´qzñðêiWµÄ«½æ¤ÉCªq®Íw@ào±IÈèÉÛÂæ¤ÈS©ððÝ¯é±ÆàÂ\Åè@©çCño±IÈè@ÖÆiWµC»Ì¸éBxÌüãÆÎÛÌgåðÚwµÄ­WµÄ«Ä¢éB®(1)É¦µ½ 2 Ì|eVÍCIðC. |eVzèµÄ¢é½ßÉC±êçÌ|eVðíXÌÝèðÌÅCÀÛÌ¨¿Ì´qzñâ¨«ÌÄ»«ðÚIÆ·éÆCp¢é|eVÍCÅàdvÈvfÅ

éBã\IÈp¢½vZÌÎÛÆµÌKXÉÍCéöxÌCI«ªúÒÅ«ég¬nªKµÄ¢éB»ÌáÆµÄCÅà½­ÌñáÌéP2 Ì|eVÌ Born­Mayer­Huggins |eC_nðl¦éÆCSi­O Í¤LEVÍCÃdÍCdq_Ì½­CoÉq_CIªÆàÉ 50÷öxÆ©ÏàçêÄoÉqÝìpCoÉq_ldÉqÝìpÌ¢éB1980 NãÌã¼È~É»ÝÅàp¢ç©çÈéÈºÌ®(1)Å\³êéBêÄ¢é SiO2 ðÎÛÆµ½|eVªñq(r )

Zi Zj{Aij b expr
Cij Dij| 6| 8rr()

³êÄ¢éB

si{sj|rr1)(1)

TTAM3)

SiO44| CIÌüÍÉd×ða·é½ß
É¼z_d× e{ ðzuµCab initio Harrtree­±±ÅCr Í±q i Æ j ÌÔÌ£ÅCZ Í±Fock SCF MO vZðp¢ÄC3 íÌÏ`[qÌd×CAij Í|[OöqCb ÍËÍÌhÉ¨¯é|eVGlM[ðßÄ¢å«³ð\·p[^Cs ÍCIÌå«³éB±Ì|eVGlM[©çCÈºÌ64NEW GLASS Vol. 19 No. 4 2004®(2)Ìep[^ðÅ¬©æ@ÉæèÅK|eVÖÅCrij Í±q i Æ j ÌÔÌ»µC³çÉC¾çê½p[^ZbgÉÎ£C ujik ÍpxÌÝèlÅC|eVÌlµÄCSiO2 »ÌÄ»\Í©çÅIIÈpÍCÝèµ½pxÅ 0 ÆÈèCpxª¸êé[^ÌlðèµÄ¢éBSiO2 »\¢ÌÙÇlªå«­ÈéBf(rij, rik jÍCÌÅ¼É]Ú·xE³ÍÉæéiqèÌÏ»ð¿Øè£ ro ÉßÃ­Æ±Ì|eVÌÄ»µÄ¢éBlð 0 Éû©·é½ßÌàÌÅq(r )

Zi Zj{Aij b expr(s {sr |r)| ri

Cijj6(2)2)

BKS4)éBãLÈOÉà SiO2 ÌKXóÔâ»¿ÌV~[VÍCáª½­C½¬ªnâvZÉæÁÄ¾çêé¨«â¨«­»@\ÌðÍâªõwIXyNgÖÆiñÅ«Ä¢éBæÌ TTAM Æ¯lÉªqO¹vZðîÉµ4)d×½t»@6)

ÄC|eVp[^ðÝèµÄ¢éBãLÌ 3 ÂÌ|eVÌp[^ÍªqO¹@ÌvZÉÍCH4SiO4 ÈÇÌNX

üÍÌ±qÌzuÉæÁÄÏ»µÈ¢B±êÉ^ðp¢Ä¢éBÅIIÉÍCªqO¹@©çÎµÄCd×½t»@ÅÍCüÍÌ´qÆÌÔß½p[^ðp¢ÄÎp»ÌiqèÅÌdqÌÚ®ðl¦é±ÆÉæèCüÍÌÂÆÌÏe«¦ðßÄCp[^Ìlðè

«É·é|eVÆµÄúÒ³êéBµÄ¢éBp[^ÌlÍæÌ TTAM ÉßÜ½CSiO2 KXÖÌKpàñ³êÄ¢é7)B¢lðñµÄ¢éBÜ½CSi­O ¾¯Åªq®Íw@ðp¢ÄClXÈ\¢ðÍè@È­CNa­O CAl­O CP­O ÉÖµÅ¾çêéîñðµÄC3 ³ÌfðÄàLÚ³êÄ¢éB¾é±ÆªÅ«C»êªíXÌ¨«ÆeÕÉq±êçÍCCIß©çCSi â O Ìd×Ìlð¬³­µCr Ìlð¬³­·é±Æªéæ¤ÉÈêÎCÞ¿J­EÝvÌxc[ÌÐÆÂÉÈéÅë¤BÅC¤L«ðæèüêé±ÆªÅ«é±Æð¦µÄ¢éBoûÆà»ðîÉµÄ¢éªC
SiO2 ÌKXóÔâ¼Ì¬ªðYÁµ½K
Xg¬ÉL­p¢çêÄ¢éB³çÉCMorseðÁ¦é±ÆÅ¤L«ðæèüêé±ÆCIß|eVÉæéKX`¬ªq®Íw@ÉæévZÍCÀÝÌ¨¿ðÍí·éisimulatej±Æ¾¯ªÚIÅàsíêÄ¢éB3)D.

3 Ì|eV5)ë¤©BP»µ½|eVðp¢½vZÌÉC2 Ì|eVÉ Si­O­Si â O­Si­O ÌKXóÔðl¦éûÍÈ¢ÌÅë¤©BpÉË¶µ½ÈºÌ®(3)Ì|eVðtPÈCIßÌ|eVðp¢ÄàCÁµÄ¢éBSi Ép¢çêÄ¢½ÈºÌ®(3)ZÌ©ç}âµÄCKXEñ»¿óÔª¾ç

Å\³êé|eVÖð SiO2 ÉKpµêéBâp¬xª 1014 K/s ÈÇÆÉßÄ¬¢Ä¢éB±ÆÆCî{ZªÉßÄ¬³¢½ßÉCÕE

V3Bjik f(rij, rik )p(ujik, ujik jf(rij, rik )expjª¶¬·éæ¤Èå«Èhç¬ÌN±ém¦(r |r {r |r )lijªÉßÄ¬³¢½ßÉC»»ªN±çÈ¢Bl0ikÀÝÌKXÆÍCå«Èu½èª0for rij, rikr0, 0 for rij, rikr0p(ujik, ujik )(cos ujik|cos ujik)2(3)±êÍC±q j, i, k Ììépx ujik ÉË¶µ½

éªCKXóÔÉ¨¯élXÈÁ¥ðLµÄ¢éÌÅÍÈ¢©BÉßÄàIEÃTIÈàeÅéªC±±ÅÍCPê¬ªÌñ»¿óÔÌCI
65

NEW GLASS Vol. 19 No. 4 2004ÌzñðìÁÄÝæ¤BïÌIÉÍCî{ZÌÉzCIð 200ÂÆ»ÌzCIÌd×É©¤¾¯ÌACIðzu·éBúzuÉÍCðp¢CúzuÌe¿ð³¦é½ßÉC1 ÂÌðÉÎµÄ 10 ÂÌÙÈéúzu©çÌvZðs
ÈÁ½BêAÌvZÌÅÍCî{ZÌÌÏÍêèÆµÄC·x¾¯ð§äµC3000 K Å
10 ps Ûµ½ãÉC293 K ÜÅ 10 ps ÅâpµC293 K Å 5 ps ÛµCÅãÌÀW©çC³ÍâCIÌzuÌf[^ð¾½Bî{ZÌå«³ðÏ¦ÄCvZ·é±ÆÉæèC293
K É¨¢ÄC0 GPa Ì³Íð^¦éZTCYÆ»ÌÆ«Ìzu©ç\¢îñð¾½B|eVÍæÌ®(1)ÌEÓÌnßÌ 2 ÌÝðp¢½BFig. 2

The molar volume (a) and coordination
number around cations (b) of the structuralmodels for C4{A2|2Fig. 3

The molar volume (a) and coordination
number around cations (b) of the structuralmodels for C3{A2|3/2Ü¸C{4 ¿ÌzCI C4{ Æ|2 ¿ÌACI A2| Æ·é CA2 ÌP¬ªð©ÄÝæ¤BA­A Ì|eVÌ½­WÈÇÍÅèµÄCC­A Ì|eVÌ½­W¾¯Ï¦
é±ÆÅCzCIÌå«³ðÏ¦½êÌCIzñÆÌÏÌÏ»ðß½B} 2 ÉCA2 g¬É¨¯é C­A £ÆÌÏÌÖWÆ C4{ CIÌzÊÌÖWð¦µ½BzCI C4{ ÌCI¼aª¬³¢Æ«ÍC
ÙÚ 4 ÊÌ\¢ðÆèCCI¼aªå«­ÈéÉµ½ªÁÄCÌÏªÁ·éªCz
ÊÌÏ»ÍÙÆñÇÈ¢B±êªCC­A Ì£ª 1.85 Å ð´¦éÆCÌÏª}É¸­µCzÊªÁ·éB±êæèå«¢zCI C4{ ÅÍCÌÏÆzÊÌoûªÁ·éB±êÍCC4{ CIª¬³¢Æ«É CA4 ÌlÊÌ\¢ªÀèÅèCC4{ CIªå«­ÈéÆCCA4 ÌlÊÌ\¢ðÛÅ«È­ÈèCCIÌzñªå«­Ï»·éBCI«»ÌêÍCCA4 ÌlÊÌÌÉ CA6 ÌªÊÌªÀèÉÈéªC±ÌvZÅÍCzÊª 6 zÊÅêèÆÈéÌæÍÈ
¢BÉ±êð{3 ¿Cé¢ÍC{5 ¿ÌzÌzCIÌêÆ¯¶æ¤ÉzCIÌ¼aCIÆµ½êð} 3 Æ 4 É¦µ½B{4 ¿ª¬³¢êÉÍCzCIÌ¼aªÁµÄ66

NEW GLASS Vol. 19 No. 4 2004

Fig. 4The molar volume (a) and coordination
number around cations (b) of the structuralmodels for C5{A2|5/2.

Fig. 5The molar volume (a) and coordination
number around cations (b) of the structuralmodels for C2{A1|2.

àCzÊªêèÅCÌÏªÁ·éÌæÆl¦çêéBzÊ½ÊÌª\¢PÊÆÈÁ½ª

é±ÆC±ÌzÊð´¦éÆÌÏÆP âòQ ÍCÜ³µ­ Zacha­ñ»¿Ì´qzñÌò

zÊÌoûªÁ·é±Æªí©éBÜ½Criasen ÌKX`¬¥ð¦µÄ¢éB³çÉC{5 ¿ÌzCIÌêÍC4 zÊÉ¨¢ÄàR ÍACI¼að 1.30 Å Æµ½Æ«CCIòéBzCI¼aªå«¢Ìæ»ÌC3E4 zÊÌÕE¼aäizCIÅÍC3 ÂÌnÅÙÚ¯¶ACIÌ§xÆÆACIÌ¼aÌäjð´¦éÆCÌÏÁXüð¦µÄ¢é±Æªí©éBµ½ªÁÌûkªN«éÆ¢¤æ¤É C­A £ÅÄCÌÏªûkµ½ãÌ\¢ÍC[U³ê®Å«éBÀèÈÌæª½ACIÌÔÉzCIªÊu·éæ¤ÈC
ACIÌ[UÉæé\¢Æl¦çêéBÌ C4{A2|2 ÌêÌÌÏâzÊÆÙÆ±±ÜÅÌÊÍC
Pò

QòR
ò

ÉCACIÌd×ª|1 ¿ÌêðÝÄæ¤BzCIª{2 ¿Ìêi} 5jÍCæzCI¼aª¬³¢êCzCIÌñÇ¯¶Ï»ð¦µÄ¢éB³çÉ} 6 Æ 7ACIzÊÍ¬³­CACIÌìéÉ¦·æ¤ÉCzCIÌd×ð{3 ¿C{4zÊ½ÊÌð\¢PÊÆµ½\¢ÆÈéB¿Æâ·ÆCzÊª 6 â 8 ÅêèÆÈé±ÌzÊÍC1 ÂÌACIÌüèÉÌæªé±Æªí©éBµ½ªÁÄCACI2 ÂÌzCIªµÄ¢é±ÆÉÌ¿Ìá¢ð´¦ÄCzÊ½ÊÌð\¢P·éBÊÆ·é\¢ÆACIÌ[UÉæé\¢ÆªéCIQ ÌACIÌüèÉ 2 ÂzC¶ÝµCæÌò¼aÌÆ±ëÅå«­ÌÏªûkIªµÄ¢éóÔªÀèÈzÊð¦·µCACIÌ[UÉæé\¢ÆÈéB½ßÌdvÈvÅCzÊª¬³¢±ÆÍzCIÌ¼aªÁ·éÆC67NEW GLASS Vol. 19 No. 4 2004

Fig. 8Fig. 6

The molar volume (a) and coordination
number around cations (b) of the structuralmodels for C3{A1|3.

Average deviation from the coordinates ofions in the structural model for C4{A2|2 after
4000 time steps. The magnitude of the appliedelectric eld were ü : 1, ¢: 3, : 4 and ¤: 5GV/m.

VÖÌ«¿©çC|eVÌÉ¬ÆÈé£ð·­·éÆC[³ªó­Èé½ßÉCzCIÆACIÌªã­ÈèzÊ½ÊÌÌÀè«ªºªé±Æªl¦çêéB³çÉCzÊªÁ·é±ÆÉæèCACIÌüèÉ 3 ÂÌzCIªzu³êCd×Ì
å«ÈzCIÔÌ½­Ì½ßÉCV½ÈzCIÆACIÌªÅ«é±ÆÉæé|eVÌáºª¬³­Èé½ßÉCzÊ½Ê
ÌÌÀè«ªáºµCACIÔÌ½­Éæè\¢ª`¬³êéÆl¦çêéB±±ÜÅÌÊÅÍCzÊ½ÊÌð\¢PÊ
Æ·é\¢ªCACIª|2 ¿Æµ½Æ«É
Zachariasen ÌKX`¬¥ÉÎµCACI
ð|1 ¿Æ·éÆCzÊªÁ·éûüÉFig. 7

The molar volume (a) and coordination
number around cations (b) of the structuralmodels for C4{A1|4.

g£³êé±Æð¦µÄ¢éB³çÉC±ÌCIÌzuÉÎµÄCdêðóÁµÄCªóêCCIª®«o·dêÌlð²×½B}
8 É{4 ¿Æ|2 ¿ÌgÝí¹É¨¢Ä 1`5{4 ¿Æ|2 ¿C{3 ¿Æ|2 ¿ÈÇÌgÝ

GV/m ÌdêðóÁµÄ 4000 step oßµ½Oí¹ÉæéÆ±ëªå«¢±Æªí©éB|eãÅÌeCIÌÚ®£Ì½Ïð¦µ½B

68

NEW GLASS Vol. 19 No. 4 2004£ª 1.85 Å ÜÅÌzÊ½ÊÌð\¢PÊQl¶£Æ·é\¢ÌûªCACIÌ[UÉæé\¢1) A. Rahman, R. H. Fowler, A. H. Narten, J. Chem.

Éä×CdêÉÎ·éÏ«ª¢±Æªí©ÁPhys., 57 3010 (1972); L. V. Woodcock, C. A.

½BAmgell, P. Cheeseman, J. Chem. Phys., 65, 1565(1976).
2) ìYÇKCOãvOCåì©¨éCgRsE.

^EV~¨íèÉ

[VÉ¨¯é¨¿Èw[ªq®

Íw@ÆeJ@h¤§oÅi1996jEãªq®Íw@ÉÂ¢ÄC»Ìè@ÌTvðàc

¾µ½BÚ×ÍQlÐðQÆ³ê½¢B´qãcÔ|eVÌ¸§³ªvZÌLø«ðèqXi1990jEÍºYsC
gp\RªqV~µÄ¢éÆ·éÈçÎCÂ\ÈÍÍÅ¸xÌ¢±Æª]Üµ¢ªC»óÅÍeÕÈ±ÆÅÍÈ¢BÚIâpÌÝèÉæèClXÈxÌvZÌÉ©o·àÌªéÌ¾Æv¤Bá

ûC
gªqV~
ûCgRs[VhÖØ[i2003jE[^V~[Vh©

[Vªq®ÍwÀ±üåhC¶°i1990jEªcMCåVfñCgªqV~

[VüåhC¶°i1989jD3) S. Tsuneyuki, M. Tsukada, H. Aoki, Y. Matsui,Phys. Rev. Lett., 61, 869 (1988).

¦ÎCnZ}â@ÅÍì»¢ïÈg¬É¨¢Ä4) B. W. H. van Beest, G. J. Kramer, R. A. vanàCtâCðîµÄCKXEñ»¿ªì

Santen, Phys. Rev. Lett., 64, 1955 (1990).

»³êÄ¢éB»Ì´qzñÆÀè«ÌÖWÈÇÌm©Í­È¢Bã¼ÅÐîµ½àeÍCÜèÉP»µß¬Ä¢é©àµêÈ¢ªC±êð 2 ¬ªnÈÇÉg£µ½Æ«C³çÉL
¢KX`¬âÀè»Ìûªé©àµêÈ

¢BP»µ½©ç±»©¦Ä­é±ÆªÈ¢¾ë¤©B»±©çCÉßÄ¡GÈÀÝÌ¨¿âÞ¿ÖÌðð[ßé±ÆàìÉ¢êê
ÎCªq®Íw@ðp¢½vZ@À±ÌÂ\«5) F. H. Stllinger, T. A. Weber, Phys. Rev. B 315262 (1985); B. Feuston, S. H. Garofalini, J.
Chem. Phys., 89 5818 (1988).; P. Vashista, R. K.
Kalia, J. P. Rino, I. Ebbsj äo, Phys. Rev. B 41,12197 (1990).
6) A. K. Rappe, W. A. Goddard III, J. Chem. Phys.,95, 3358 (1991).
7) E. Demiralp, T. Ca ²gin, W. A. Goddard, Phys.
Rev. Lett., 82, 1708 (1999); C. Wang, N. Kuzuu,
Y. Tamai, J. Appl. Phys. 92, 4408 (2002).

Í³çÉLªéÆv¤B69