116
J. Liu et al.
(*Define the parameters*)
(*constants*)
h = 6.62606896*10ˆ(-34)/(2 Pi);
\[Nu] = 7.29*10ˆ12; (* unit:Hz *)
\[Omega] = \[Nu]*2 Pi;
Tm = 350;(* unit:K *)
Subscript[f, -2] = 3.737;
Subscript[f, -1] = 1.699;
Subscript[f, 2] = 0.479;
m = 63.546/(6.02*10ˆ23)*10ˆ(-3);(* unit:Kg *)
a = 0.361*10ˆ(-9);(* unit:m *)
(*parameters*)
Subscript[\[Alpha], 200] = 0.175;
Subscript[\[Alpha], 111] = Sqrt[3]/2*Subscript[\[Alpha], 200];
Subscript[\[Beta], 111] = 4;
Subscript[\[Beta], 200] = 4;
Subscript[\[Rho], 111] = (3/2 + 3/6)/(Sqrt[6]/2*Sqrt[2]/2*1ˆ2);
Subscript[\[Rho], 200] = 1/1;
Subscript[k, 111] = Sqrt[2]*2 Pi/a;
Subscript[k, 200] = 2 Pi/a;
(*Debye-Waller factor*)
J[T_] := 2 Subscript[f, -2] (T/Tm) + 1/6 (T/Tm) -
Subscript[f, 2]/360*(T/Tm)ˆ3;
Subscript[\[Sigma], 111] = 4 h/(m*\[Omega])*J[273 + 800];
Subscript[\[Sigma], 200] = 4 h/(m*\[Omega])*J[273 + 800];
(**correlation function**)
k0 = 1;
k1 = Sqrt[2] k0;
k2 = k0;
Subscript[C, 2 _ 111] =
Exp[-Subscript[\[Sigma], 111]
Subscript[k,
111]ˆ2/(2 Subscript[\[Rho], 111] Subscript[\[Beta],
111])] Exp[-(k - k1)ˆ2/(2 Subscript[\[Alpha], 111]ˆ2)];
Subscript[C, 2 _ 200] =
Exp[-Subscript[\[Sigma], 200]
Subscript[k,
200]ˆ2/(2 Subscript[\[Rho], 200] Subscript[\[Beta],
200])] Exp[-(k - k2)ˆ2/(2 Subscript[\[Alpha], 200]ˆ2)];
Subscript[C, 2] = Subscript[C, 2 _ 111] + Subscript[C, 2 _ 200];
p1 = Plot[Subscript[C, 2 _ 111], {k, 0, 3}, PlotRange -> 1.5,
PlotStyle -> Blue];
p2 = Plot[Subscript[C, 2 _ 200], {k, 0, 3}, PlotRange -> 1.5,
PlotStyle -> Blue];
p3 = Plot[Subscript[C, 2], {k, 0, 3}, PlotRange -> 1.5,
PlotStyle -> Red];
(*two-mode PFC model [35]*)
Q1 = Sqrt[2/1];
Cwu[k_] := -r - Bx (1 - (k)ˆ2)ˆ2 ((Q1ˆ2 - (k)ˆ2)ˆ2 + R1/Bx);