level 1
圣德太子666
楼主
求助,可能是哪里出问题了、
输入:
f1[b_] = GyromagneticRatio["1H"] (1 - δ₁ 10⁻⁶) b;
f2[b_] = GyromagneticRatio["1H"] (1 - δ₂ 10⁻⁶) b;
f3[b_] = GyromagneticRatio["1H"] (1 - δ₃ 10⁻⁶) b;
ave[b_] := (f1[b] + f2[b] + f3[b]) / 3.;
Ham0[b_] = (f1[b] - ave[b]) opI[1, "z"] + (f2[b] - ave[b]) opI[2, "z"] + (f3[b] - ave[b]) opI[3, "z"] + 2 π (J12 opI[1].opI[2] + J23 opI[1].opI[3]);
{ev1, vecs1} = CorrelatedEigensystems[Ham0[b] /. {J12 -> 25., J23 -> 5, δ₁ -> +0, δ₂ -> 10, δ₃ -> -10}, {b, 0.00, .25, 1/100}];
2026年03月18日 11点03分
1
输入:
f1[b_] = GyromagneticRatio["1H"] (1 - δ₁ 10⁻⁶) b;
f2[b_] = GyromagneticRatio["1H"] (1 - δ₂ 10⁻⁶) b;
f3[b_] = GyromagneticRatio["1H"] (1 - δ₃ 10⁻⁶) b;
ave[b_] := (f1[b] + f2[b] + f3[b]) / 3.;
Ham0[b_] = (f1[b] - ave[b]) opI[1, "z"] + (f2[b] - ave[b]) opI[2, "z"] + (f3[b] - ave[b]) opI[3, "z"] + 2 π (J12 opI[1].opI[2] + J23 opI[1].opI[3]);
{ev1, vecs1} = CorrelatedEigensystems[Ham0[b] /. {J12 -> 25., J23 -> 5, δ₁ -> +0, δ₂ -> 10, δ₃ -> -10}, {b, 0.00, .25, 1/100}];



