熱膨張解析

//============== for Magmetostriction ====================================================================== 

//ME-０// Load Wave-----------------------------------------------------------------------------------

1. Macro PPMS_ME0_Load_Wave_()

2. LoadWave/J/D/W/K=0/L={0,1,0,0,0}

3. End

//ME-1// Set Folder Pass-----------------------------------------------------------------------------------

1. Macro PPMS_ME1_SetFolderPass()

2. //-●●●check●●●-// Folder Pass

3. string path2Kp = "root:PPMS:'yyyymmdd_Sampname':B100_para_L100:B_sweep:'2K':sage:"

4. string/G path2Kp

5. Setdatafolder path2Kp

6. End 

//ME-2// Prepare dL_L0--------------------------------------------------------------------------------------

1. Macro PPMS_ME2_Prepare_dL_L0()

3. //unit conversion------------------------------------------------

4. Duplicate/O deltaL_1E_6cm_, dL_samp

5. Duplicate/O deltaL_1E_6cm_, dL_L0

6. 'dL_samp'=('deltaL_1E_6cm_')*10^(-8) //[m]

8. //change Sample length-------------------------------------------

9. //-●●●check●●●-// 

10. Variable/G 'L_0'

11. 'L_0'=#####*10^(-3); // Sample Length mm 10^(-3)[m] at T(room)

12. //calc-----------------------------------------------------------

13. 'dL_L0'=('dL_samp')/'L_0'

14. End

//ME-3// Display ΔL/L0 vs B(T) ----------------------------------------------------------------------------

1. Macro PPMS_ME3_Display_Bsweep_()

3. //-●●●check●●●-// Window Name

4. string WZ2K = "La05zp_2K"

5. string/G WZ2K

7. Execute "Display/N=" + WZ2K + " dL_L0 vs Field_T_"

9. //Modify

10. Execute "ModifyGraph/W=" + WZ2K + " tick=2,noLabel=0,mirror=1,fSize=20,axThick=2,standoff=0,font=\"Arial\""

11. Execute "ModifyGraph/W=" + WZ2K + " width=340.157,height=255.118"

12. //Left

13. Execute "Label/W=" + WZ2K + " left \"(1/L\\B0\\M)*ΔL(\\um)\""

14. Execute "ModifyGraph/W=" + WZ2K + " manTick(left)={0,0.5,-6,2},manMinor(left)={1,0};DelayUpdate"

15. Execute "SetAxis/W=" + WZ2K + " left 0,400;DelayUpdate"

16. //Bottom

17. Execute "Label/W=" + WZ2K + " bottom \"\\f02B\\f]0 (T)\""

18. Execute "ModifyGraph/W=" + WZ2K + " manTick(bottom)={0,4,0,0},manMinor(bottom)={1,0};DelayUpdate"

19. Execute "SetAxis/W=" + WZ2K + " bottom 0,12"

20. End

//ME-4// Differentiate to obtain magnetostriction ------------------------------------------------------

1. Macro PPMS_ME4_smth_Dif_()

2. 
   

3. //-●●●check●●●-// Folder Pass

4. string path2Kp = "root:PPMS:'yyyymmdd_Sampname':B100_para_L100:B_sweep:'2K':sage:"

5. string/G path2Kp

6. SetdataFolder path2Kp

7. NewdataFolder/O DIF

8. Execute "SetdataFolder " + path2Kp + "DIF"

9. NewdataFolder/O smth31

10. NewdataFolder/O smth81

11. NewdataFolder/O smth121

12. NewdataFolder/O smth221

13. NewdataFolder/O smth451

14. NewdataFolder/O DIF_Cut

15. 
    

16. //smoothing_31

17. Execute "SetdataFolder " + path2Kp + "DIF:smth31"

18. Duplicate/O :::dL_L0,dL_L0_smth31;DelayUpdate

19. Duplicate/O :::Field_T_,Field_T_smth31;DelayUpdate

20. Smooth/S=2 31, dL_L0_smth31,Field_T_smth31;DelayUpdate

21. //Differentiate

22. Differentiate/METH=1 dL_L0_smth31/X=Field_T_smth31/D=dL_L0_smth31_DIF;DelayUpdate

23. //graph

24. //-●●●check●●●-// Window Name

25. string WZ2KD = "La06zp_2K_dif"

26. string/G WZ2KD

28. Execute "Display/N=" + WZ2KD + "  dL_L0_smth31_DIF vs Field_T_smth31"

29. //Modify

30. Execute "ModifyGraph/W=" + WZ2KD + " tick=2,noLabel=0,mirror=1,fSize=20,axThick=2,standoff=0,font=\"Arial\""

31. Execute "ModifyGraph/W=" + WZ2KD + " width=340.157,height=255.118"

32. //Left

33. Execute "Label/W=" + WZ2KD + " left \"(1/L0)*(dΔL/dB) \\u\""

34. Execute "ModifyGraph/W=" + WZ2KD + " manTick(left)={0,1,-7,0},manMinor(left)={1,0};DelayUpdate"

35. Execute "SetAxis/W=" + WZ2KD + " left -2e-07,2e-07;DelayUpdate"

36. //Bottom

37. Execute "Label/W=" + WZ2KD + " bottom \"\\f02B\\f]0 (T)\""

38. Execute "ModifyGraph/W=" + WZ2KD + " manTick(bottom)={0,4,0,0},manMinor(bottom)={1,0};DelayUpdate"

39. Execute "SetAxis/W=" + WZ2KD + " bottom 0,12"

41. //smoothing_81

42. Execute "SetdataFolder " + path2Kp + "DIF:smth81"

43. Duplicate/O :::dL_L0,dL_L0_smth81;DelayUpdate

44. Duplicate/O :::Field_T_,Field_T_smth81;DelayUpdate

45. Smooth/S=2 81, dL_L0_smth81,Field_T_smth81;DelayUpdate

46. //Differentiate

47. Differentiate/METH=1 dL_L0_smth81/X=Field_T_smth81/D=dL_L0_smth81_DIF;DelayUpdate

48. //graph

49. Execute "AppendToGraph/W=" + WZ2KD + " dL_L0_smth81_DIF vs Field_T_smth81"

50. 
    

51. //smoothing_121

52. Execute "SetdataFolder " + path2Kp + "DIF:smth121"

53. Duplicate/O :::dL_L0,dL_L0_smth121;DelayUpdate

54. Duplicate/O :::Field_T_,Field_T_smth121;DelayUpdate

55. Smooth/S=2 121, dL_L0_smth121,Field_T_smth121;DelayUpdate

56. //Differentiate

57. Differentiate/METH=1 dL_L0_smth121/X=Field_T_smth121/D=dL_L0_smth121_DIF;DelayUpdate

58. //graph

59. Execute "AppendToGraph/W=" + WZ2KD + " dL_L0_smth121_DIF vs Field_T_smth121"

61. //smoothing_221

62. Execute "SetdataFolder " + path2Kp + "DIF:smth221"

63. Duplicate/O :::dL_L0,dL_L0_smth221;DelayUpdate

64. Duplicate/O :::Field_T_,Field_T_smth221;DelayUpdate

65. Smooth/S=2 221, dL_L0_smth221,Field_T_smth221;DelayUpdate

66. //Differentiate

67. Differentiate/METH=1 dL_L0_smth221/X=Field_T_smth221/D=dL_L0_smth221_DIF;DelayUpdate

68. //graph

69. Execute "AppendToGraph/W=" + WZ2KD + " dL_L0_smth221_DIF vs Field_T_smth221"

70. 
    

71. //smoothing_451

72. Execute "SetdataFolder " + path2Kp + "DIF:smth451"

73. Duplicate/O :::dL_L0,dL_L0_smth451;DelayUpdate

74. Duplicate/O :::Field_T_,Field_T_smth451;DelayUpdate

75. Smooth/S=2 451, dL_L0_smth451,Field_T_smth451;DelayUpdate

76. //Differentiate

77. Differentiate/METH=1 dL_L0_smth451/X=Field_T_smth451/D=dL_L0_smth451_DIF;DelayUpdate

78. //graph

79. Execute "AppendToGraph/W=" + WZ2KD + " dL_L0_smth451_DIF vs Field_T_smth451"

81. //ten_sen

82. Execute "ModifyGraph/W=" + WZ2KD + " mode=3,marker=19,msize=1,rgb(dL_L0_smth31_DIF)=(52428,1,1),rgb(dL_L0_smth81_DIF)=(52428,34958,1),rgb(dL_L0_smth121_DIF)=(26205,52428,1),rgb(dL_L0_smth221_DIF)=(1,34817,52428),rgb(dL_L0_smth451_DIF)=(1,4,52428)"

83. Execute "DrawText/W=" + WZ2KD + " 0.5,0.78,\"\\Z20 \\f02T\\f]0 = 2K_sage\""

84. Execute "DrawText/W=" + WZ2KD + " 0.5,0.78,\"\\Z20 LaMgNi4 #5\""

85. Execute "DrawText/W=" + WZ2KD + " 0.5,0.78,\"\\Z20 \\f02L\\f]0 || [100] || \\f02B\\f]0 ||  [100]\""

86. End

87. 
    

//ME-5// hukusei ------------------------------------------------------

1. Macro PPMS_ME5_DIF_Cut()

2. //-●●●check●●●-// Folder Pass

3. string path2Kp = "root:PPMS:'yyyymmdd_Sampname':B100_para_L100:B_sweep:'2K':sage:"

4. string/G path2Kp

5. Execute "SetdataFolder " + path2Kp + "DIF:DIF_Cut"

6. NewdataFolder/O smth31

7. NewdataFolder/O smth81

8. NewdataFolder/O smth121

9. NewdataFolder/O smth221

10. NewdataFolder/O smth451

11. 
    

12. //Duplicate and move wave.

13. Execute "SetdataFolder " + path3Kp + "DIF:smth31"

14. Duplicate/O dL_L0_smth31_DIF, ::DIF_cut:smth31:dL_L0_smth31_DIF

15. Duplicate/O Field_T_smth31, ::DIF_cut:smth31:Field_T_smth31

17. Execute "SetdataFolder " + path3Kp + "DIF:smth81"

18. Duplicate/O dL_L0_smth81_DIF, ::DIF_cut:smth81:dL_L0_smth81_DIF

19. Duplicate/O Field_T_smth81, ::DIF_cut:smth81:Field_T_smth81

20. 
    

21. Execute "SetdataFolder " + path3Kp + "DIF:smth121"

22. Duplicate/O dL_L0_smth121_DIF, ::DIF_cut:smth121:dL_L0_smth121_DIF

23. Duplicate/O Field_T_smth121, ::DIF_cut:smth121:Field_T_smth121

24. //-●●●check●●●-// Specify range smth221[##1##,##1##]

25. Execute "SetdataFolder " + path3Kp + "DIF:smth221"

26. Duplicate/O/R=(##1##,##1##) dL_L0_smth221_DIF, ::DIF_cut:smth221:dL_L0_smth221_DIF

27. Duplicate/O/R=(##1##,##1##) Field_T_smth221, ::DIF_cut:smth221:Field_T_smth221

28. //-●●●check●●●-// Specify range smth451[##2##,##2##]

29. Execute "SetdataFolder " + path3Kp + "DIF:smth451"

30. Duplicate/O/R=(##2##,##2##) dL_L0_smth451_DIF, ::DIF_cut:smth451:dL_L0_smth451_DIF

31. Duplicate/O/R=(##2##,##2##) Field_T_smth451, ::DIF_cut:smth451:Field_T_smth451

32. //-●●●check●●●-// Specify range smth###

33. Concatenate {::DIF_Cut:smth451:dL_L0_smth451_DIF,::DIF_Cut:smth221:dL_L0_smth221_DIF},::DIF_Cut:dL_L0_smth00_DIF

34. Concatenate {::DIF_Cut:smth451:Field_T_smth451,::DIF_Cut:smth221:Field_T_smth221},::DIF_Cut:Field_T_smth00

35. Execute "AppendToGraph/W=" + WZ3KD + " ::DIF_Cut:dL_L0_smth00_DIF vs ::DIF_Cut:Field_T_smth00"

36. End

//ME-6// graph_DIF_Cut ------------------------------------------------------

1. Macro PPMS_ME6_graph()

2. //-●●●check●●●-// Folder Pass

3. string path2Kp = "root:PPMS:'yyyymmdd_Sampname':B100_para_L100:B_sweep:'2K':sage:"

4. string/G path2Kp

5. Execute "SetdataFolder " + path2Kp + "DIF:DIF_Cut"

6. 
   

7. //-●●●check●●●-// Window Name

8. string WZ2KDC = "La07zp_2K_dif"

9. string/G WZ2KDC

10. 
    

11. //-●●●check●●●-// Wave Name

12. Execute "Display/N=" + WZ2KDC + "  dL_L0_smth00_DIF vs Field_T_smth00"

13. //Modify

14. Execute "ModifyGraph/W=" + WZ2KDC + " tick=2,noLabel=0,mirror=1,fSize=20,axThick=2,standoff=0,font=\"Arial\""

15. Execute "ModifyGraph/W=" + WZ2KDC + " width=340.157,height=255.118"

16. //Left

17. Execute "Label/W=" + WZ2KDC + " left \"(1/L0)*(dΔL/dB) \\u\""

18. Execute "ModifyGraph/W=" + WZ2KDC + " manTick(left)={0,1,-7,0},manMinor(left)={1,0};DelayUpdate"

19. Execute "SetAxis/W=" + WZ2KDC + " left -2e-07,2e-07;DelayUpdate"

20. //Bottom

21. Execute "Label/W=" + WZ2KDC + " bottom \"\\f02B\\f]0 (T)\""

22. Execute "ModifyGraph/W=" + WZ2KDC + " manTick(bottom)={0,4,0,0},manMinor(bottom)={1,0};DelayUpdate"

23. Execute "SetAxis/W=" + WZ2KDC + " bottom 0,12"

24. End

//ME-7// Curve Fit to find the slope-----------------------------------------

1. Macro PPMS_ME7_CurveFit()

2. //-●●●check●●●-// Folder Pass

3. string path2Kp = "root:PPMS:'yyyymmdd_Sampname':B100_para_L100:B_sweep:'2K':sage:"

4. string/G path2Kp

5. SetdataFolder path2Kp

6. NewdataFolder/O FIT

7. Execute "SetdataFolder " + path2Kp + "FIT"

9. //Curve Fit

10. CurveFit/X=1/H="10"/TBOX=768 line ::DIF:DIF_Cut:dL_L0_smth00_DIF[pcsr(A),pcsr(B)] /X=::DIF:DIF_Cut:Field_T_smth00 /D

12. //The line is made into a function.

13. Duplicate/O W_coef, X_fit;DelayUpdate

14. X_fit[0]=0;

15. X_fit[1]=4;

16. Duplicate/O W_coef, Y_fit;DelayUpdate

17. Variable/G 'AAA'

18. 'AAA'=W_coef[1]

19. 'Y_fit'= 'AAA'*'X_fit'

20. Print/D W_coef[1]

21. //graph

22. Execute "AppendToGraph/W=" + WZ2KDC + " Y_fit vs X_fit"

23. End

//ME-8// quadrupole susceptibility-----------------------------------------

1. Macro  PPMS_ME8_katamuki()

2. string pathpb = "root:PPMS:'20230315_LaMgNi4':B100_para_L100:B_sweep:"

3. string/G pathpb

4. Execute "SetdataFolder " + pathpb + "katamuki"

5. //-●●●check●●●-// Folder Pass

6. Execute "Duplicate/O " + pathpb + "'2K':sage:FIT:W_coef, " + pathpb + "katamuki:wave0"

7. Execute "SetdataFolder " + pathpb + "katamuki"

9. Make/O/N=8 katamuki

10. Make/O/N=8 Temp

11. //-●●●check●●●-// katamuki[#]

12. katamuki[#]=wave0[1]

13. Temp[0]=2;Temp[1]=3;Temp[2]=5;Temp[3]=7;Temp[4]=10;Temp[5]=15;Temp[6]=20;Temp[7]=30

14. End