{VERSION 5 0 "APPLE_PPC_MAC" "5.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 1 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 0 1 }{CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 1 } {PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }1 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Maple Output" -1 11 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }3 3 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Maple Output" -1 12 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }1 3 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Maple Plot" 0 13 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }} {SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "E_t:=z->z*p0-sqrt(z) *H0;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$E_tGf*6#%\"zG6\"6$%)operato rG%&arrowGF(,&*&9$\"\"\"%#p0GF/F/*&-%%sqrtG6#F.F/%#H0GF/!\"\"F(F(F(" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "Ev:=z*E_t(z-1)-z*E_t(z);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#EvG,&*&%\"zG\"\"\",&*&,&F'F(F(!\" \"F(%#p0GF(F(*&-%%sqrtG6#F+F(%#H0GF(F,F(F(*&F'F(,&*&F'F(F-F(F(*&-F06#F 'F(F2F(F,F(F," }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "H0:=poq*sq rt(z)*p0;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#H0G*(%$poqG\"\"\"-%%sq rtG6#%\"zGF'%#p0GF'" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "p0:= -Ec/(1-poq)/z;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#p0G,$*(%#EcG\"\" \",&F(F(%$poqG!\"\"F+%\"zGF+F+" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "assume(z>1);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "ser ies(subs(z=1/x+1,normal(Ev)),x,2);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# +)%\"xG*(%#EcG\"\"\",&%$poqG#F'\"\"#F'!\"\"F',&F,F'F)F'F,\"\"!,&F%#F,F +*(F&F',&F0F'*&#\"\"$\"\")F'F)F'F'F'F-F,F'F'-%\"OG6#F'F+" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "restart;\n#del:=(r-r0)/r0" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 56 "E_total:=del->12*A*exp(-p*de l)-B*sqrt(12*exp(-q*del)^2);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%(E_t otalGf*6#%$delG6\"6$%)operatorG%&arrowGF(,&*&%\"AG\"\"\"-%$expG6#,$*&% \"pGF/9$F/!\"\"F/\"#7*&%\"BGF/-%%sqrtG6#,$*$)-F16#,$*&%\"qGF/F6F/F7\" \"#F/F8F/F7F(F(F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 79 "BulkMo dulus:=77;\nLat_const:=4.032;\nEc:=3.39;\nLat_number:=4;\n#poq:=2.89; \nq:=2.0;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%,BulkModulusG\"#x" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%*Lat_constG$\"%KS!\"$" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#>%#EcG$\"$R$!\"#" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%+Lat_numberG\"\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"qG$\" #?!\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 105 "delta:=0.001;\n eVtoGPa:=160.218;\nOmega:=Lat_const^3/Lat_number;\ndE:=delta^2/2*9*Ome ga*BulkModulus/eVtoGPa;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&deltaG$ \"\"\"!\"$" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%(eVtoGPaG$\"'=-;!\"$" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&OmegaG$\"+>!3(Q;!\")" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#dEG$\"+d$)*Ra$!#9" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 64 "E1:=E_total(-0.001);\nE2:=evalf(E_total(0));\nE3:= E_total(+0.001);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#E1G,&*&%\"AG\" \"\"-%$expG6#,$%\"pG$F(!\"$F(\"#7*&$\"+]n.rM!\"*F(%\"BGF(!\"\"" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%#E2G,&%\"AG$\"#7\"\"!*&$\"+:;5kM!\"* \"\"\"%\"BGF.!\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#E3G,&*&%\"AG \"\"\"-%$expG6#,$%\"pG$!\"\"!\"$F(\"#7*&$\"+N.=dM!\"*F(%\"BGF(F/" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "A:=solve(E2=-Ec,A);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"AG,&%\"BG$\"+Y8v')G!#5$\"++++DGF)! \"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 107 "find_poq:=proc(tmp )\nglobal B,p,q,poq;\nB:='B';\npoq:=tmp;\np:=q*poq;\nB:=fsolve(E1=E3,B );\nevalf(E1-E2-dE);\nend;" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%)find_ poqGf*6#%$tmpG6\"F(F(C'>%\"BG.F+>%$poqG9$>%\"pG*&%\"qG\"\"\"F.F4>F+-%' fsolveG6$/%#E1G%#E3GF+-%&evalfG6#,(F:F4%#E2G!\"\"%#dEGFAF(6&F+F1F3F.F( " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 64 "poq_min:=2;\npoq_max:=1 0;\nplot('find_poq(x)',x=poq_min..poq_max);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%(poq_minG\"\"#" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%( poq_maxG\"#5" }}{PARA 13 "" 1 "" {GLPLOT2D 329 246 246 {PLOTDATA 2 "6% -%'CURVESG6$7S7$$\"\"#\"\"!$!+d$)p(=#!#97$$\"+LsPu@!\"*$!+d$)pp?F-7$$ \"+K>5EBF1$!+d$)*o'>F-7$$\"+BAt'\\#F1$!+d$)H^=F-7$$\"+kN\\oEF1$!+d$)fM \"\\Zg#\\%F1$!*d$)*y \\F-7$$\"+Fn*Gn%F1$!*d$)pv$F-7$$\"+2xiD[F1$!*d$)4s#F-7$$\"+Y,H.]F1$!*d $)f^\"F-7$$\"+2:bg^F1$!)d$)>XF-7$$\"+Y@4L`F1$\")V;!=(F-7$$\"+O;R(\\&F1 $\"*V;S$=F-7$$\"+<4#)ocF1$\"*V;]*HF-7$$\"+7lCEeF1$\"*V;?1%F-7$$\"+%G^g *fF1$\"*V;I@&F-7$$\"+>2VshF1$\"*V;!4kF-7$$\"+O&pfK'F1$\"*V;!\\uF-7$$\" +kcz\"\\'F1$\"*V;Id)F-7$$\"+\"G5Jm'F1$\"*V;St*F-7$$\"+6#32$oF1$\"+V;5( 3\"F-7$$\"+Ey'G*pF1$\"+V;?(>\"F-7$$\"+J%=H<(F1$\"+V;?>8F-7$$\"+3>qMtF1 $\"+V;!*G9F-7$$\"+,.W2vF1$\"+V;!fa\"F-7$$\"+fp'Rm(F1$\"+V;5_;F-7$$\"+T >4NyF1$\"+V;?oF- 7$$\"+od'*G$)F1$\"+V;!H5#F-7$$\"+EcB,&)F1$\"+V;q>AF-7$$\"+v>:n')F1$\"+ V;?KBF-7$$\"+0a#o$))F1$\"+V;?ZCF-7$$\"+`Q40!*F1$\"+V;IhDF-7$$\"+\"3:(f \"*F1$\"+V;5mEF-7$$\"+e%GpL*F1$\"+V;?'y#F-7$$\"+:-V&\\*F1$\"+V;!Q*GF-7 $$\"+ZhUk'*F1$\"+V;I3IF-7$$\"+ " 0 "" {MPLTEXT 1 0 46 " poq:=fsolve('find_poq(x)',x,poq_min..poq_max);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$poqG$\"+]t3F_!\"*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "A;B;p;q;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"*LVHo'! #5" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\":g7x$fmvCaE657!#C" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\":r*p*3ST(3quTX5!#B" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"#?!\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "16 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 } {PAGENUMBERS 0 1 2 33 1 1 }