In [1]:
import matplotlib.pyplot as plt
import pandas as pd
import numpy as np
In [2]:
def exp_dsc(filename):
    df = pd.read_csv(filename,index_col=0,usecols=[0,1,2,3])
    #print(df)
    l_time = np.array(df.iloc[:,0])
    l_time = np.array([float(d) for d in l_time])
    l_temp = np.array(df.iloc[:,1])
    l_temp = np.array([float(d) for d in l_temp])
    l_dsc = np.array(df.iloc[:,2])
    l_dsc = np.array([float(d) for d in l_dsc])
    print(l_time)
    print(l_temp)
    print(l_dsc)
    #dsc_max = max(l_dsc)
    #l_dsc = l_dsc/dsc_max
    return df, l_time, l_temp, l_dsc
In [3]:
df0, l_time0, l_temp0, l_dsc0 = exp_dsc("20210701_Ag2S_No.1_stoichi.csv")

[0.0000000e+00 1.8336000e-02 3.6673000e-02 ... 9.3661667e+01 9.3680000e+01
 9.3698334e+01]
[28.534    28.533001 28.530001 ... 46.118    46.105999 46.094002]
[-0.174436 -0.175011 -0.175626 ... -0.287353 -0.287209 -0.287183]
In [4]:
df2, l_time2, l_temp2, l_dsc2 = exp_dsc("20220621_No.2_Ag2S0.7Te0.3_400deg.csv")

[0.00000000e+00 4.49440000e-02 8.98880000e-02 ... 2.22876801e+02
 2.22921738e+02            nan]
[21.344999 21.353001 21.355    ... 30.003    30.             nan]
[ 0.350611  0.344334  0.344072 ... -0.268154 -0.267945       nan]
In [5]:
df3, l_time3, l_temp3, l_dsc3 = exp_dsc("20220719_#3_Ag2S0.7Te0.3_1st_400degC.csv")

[0.00000000e+00 3.83220000e-02 7.66450000e-02 ... 1.91190033e+02
 1.91228363e+02            nan]
[27.093    27.051001 27.046    ... 40.005001 40.             nan]
[-0.353894 -0.352652 -0.352665 ... -0.333348 -0.333283       nan]
In [6]:
df17, l_time17, l_temp17, l_dsc17 = exp_dsc("20211123_No.17_Ag2S0.8Te0.2.csv")

[0.00000000e+00 2.83160000e-02 5.66320000e-02 ... 1.46026718e+02
 1.46055023e+02            nan]
[21.743    21.747999 21.75     ... 30.003    30.             nan]
[ 0.511172  0.507824  0.507144 ... -0.28904  -0.288896       nan]
In [7]:
df17_2, l_time17_2, l_temp17_2, l_dsc17_2 = exp_dsc("20230313_Ag2S0.8Te0.2_#17_400degC.csv")

[0.00000000e+00 4.16280000e-02 8.32570000e-02 ... 2.08516739e+02
 2.08558380e+02            nan]
[19.514999 19.52     19.521    ... 30.004    30.             nan]
[ 0.568375  0.558737  0.557834 ... -0.196329 -0.19595        nan]
In [8]:
df20, l_time20, l_temp20, l_dsc20 = exp_dsc("20230207_#20_Ag2S0.2Te0.8.csv")

[0.00000000e+00 3.66620000e-02 7.33240000e-02 ... 1.88443359e+02
 1.88480026e+02            nan]
[19.541    19.548    19.548    ... 40.016998 40.             nan]
[-0.231941 -0.233732 -0.234098 ... -0.353227 -0.352756       nan]
In [9]:
df21, l_time21, l_temp21, l_dsc21 = exp_dsc("20230207_#21_Ag2S0.3Te0.7.csv")

[0.00000000e+00 3.66620000e-02 7.33240000e-02 ... 1.84776688e+02
 1.84813339e+02            nan]
[24.753    16.240999 19.74     ... 40.014    40.             nan]
[-0.383555 -0.420161 -0.425915 ... -0.322807 -0.32333        nan]
In [10]:
df19, l_time19, l_temp19, l_dsc19 = exp_dsc("20221010_#19_Ag2S0.6Te0.4_1day_400degC.csv")

[0.00000000e+00 3.83030000e-02 7.66050000e-02 ... 1.93006729e+02
 1.93045029e+02            nan]
[22.382    22.386999 22.389    ... 40.005001 40.             nan]
[ 0.583415  0.576615  0.576144 ... -0.396542 -0.39649        nan]
In [11]:
df29, l_time29, l_temp29, l_dsc29 = exp_dsc("20220220_#29_Ag2S0.4Te0.6_annealed.csv")

[0.0000000e+00 1.6667000e-02 3.3333000e-02 ... 3.4183334e+01 3.4200001e+01
           nan]
[ 28.608999  28.603001  28.597    ... 202.625    202.628998        nan]
[-0.279283 -0.279401 -0.279663 ... -0.719859 -0.71948        nan]
In [12]:
df42, l_time42, l_temp42, l_dsc42 = exp_dsc("20220826_#42_Ag2S0.6Te0.4_400degC_1st.csv")

[0.00000000e+00 3.82950000e-02 7.65910000e-02 ... 1.90290100e+02
 1.90328384e+02            nan]
[21.285999 21.291    21.292999 ... 40.008999 40.             nan]
[-0.201102 -0.200004 -0.198604 ... -0.40312  -0.40312        nan]
In [13]:
df42_2, l_time42_2, l_temp42_2, l_dsc42_2 = exp_dsc("20220909_#42_Ag2S0.5Te0.5_400degC_2nd.csv")

[0.00000000e+00 3.83020000e-02 7.66050000e-02 ... 1.91856735e+02
 1.91895050e+02            nan]
[21.141001 21.146    21.148001 ... 40.005001 40.             nan]
[-0.250668 -0.249387 -0.247857 ... -0.335637 -0.33548        nan]
In [14]:
df42_3, l_time42_3, l_temp42_3, l_dsc42_3 = exp_dsc("20220930_#42_Ag2S0.6Te0.4_320degC.csv")

[0.00000000e+00 3.33270000e-02 6.66530000e-02 ... 1.64733353e+02
 1.64766678e+02            nan]
[22.388    22.395    22.396999 ... 40.       40.             nan]
[-0.230253 -0.227416 -0.226592 ... -0.363206 -0.363232       nan]
In [15]:
df43, l_time43, l_temp43, l_dsc43 = exp_dsc("20220826_#43_Ag2S0.5Te0.5_400degC_1st.csv")

[0.00000000e+00 3.83110000e-02 7.66210000e-02 ... 1.90173401e+02
 1.90211700e+02            nan]
[26.282    26.247    26.245001 ... 40.006001 40.             nan]
[-0.320244 -0.314359 -0.314006 ... -0.418343 -0.418055       nan]
In [16]:
df43_2, l_time43_2, l_temp43_2, l_dsc43_2 = exp_dsc("20220909_#43_Ag2S0.5Te0.5_400degC_2nd.csv")

[0.00000000e+00 3.49890000e-02 6.99780000e-02 ... 1.80263351e+02
 1.80298355e+02            nan]
[29.726999 29.667    29.663    ... 47.895    47.883999       nan]
[-0.346505 -0.341692 -0.341169 ... -0.377722 -0.377369       nan]
In [17]:
df43_3, l_time43_3, l_temp43_3, l_dsc43_3 = exp_dsc("20220930_#43_Ag2S0.5Te0.5_320degC.csv")

[0.00000000e+00 3.16570000e-02 6.33140000e-02 ... 1.62370026e+02
 1.62401688e+02            nan]
[23.893    23.896    23.896999 ... 40.012001 40.002998       nan]
[-0.299057 -0.296756 -0.296363 ... -0.310906 -0.311089       nan]
In [18]:
df43_4, l_time43_4, l_temp43_4, l_dsc43_4 = exp_dsc("20220220_#43_Ag2S0.5Te0.5_annealed.csv")

[0.0000000e+00 1.6667000e-02 3.3333000e-02 ... 4.4416668e+01 4.4433334e+01
           nan]
[ 19.443001  19.448999  19.452    ... 201.233002 201.233002        nan]
[-0.291054 -0.290413 -0.289223 ...  0.20823   0.207968       nan]
In [19]:
df50_1, l_time50_1, l_temp50_1, l_dsc50_1 = exp_dsc("20221004_#50_Ag2S0.5Te0.5_400degC.csv")

[0.00000000e+00 3.83190000e-02 7.66380000e-02 ... 1.91173386e+02
 1.91211700e+02            nan]
[34.916    34.827    34.818001 ... 40.007    40.             nan]
[-0.306734 -0.301725 -0.301268 ... -0.309598 -0.309533       nan]
In [20]:
df51_1, l_time51_1, l_temp51_1, l_dsc51_1 = exp_dsc("20221017_#51_Ag2S0.7Te0.3_400degC.csv")

[0.00000000e+00 3.83320000e-02 7.66650000e-02 ... 1.91623352e+02
 1.91661682e+02            nan]
[21.724001 21.731001 21.731001 ... 40.019001 40.007999       nan]
[-0.235825 -0.234805 -0.233445 ... -0.297658 -0.297632       nan]
In [21]:
df53_1, l_time53_1, l_temp53_1, l_dsc53_1 = exp_dsc("20221013_#53_Ag2S0.4Te0.6_400degC.csv")

[0.00000000e+00 3.82940000e-02 7.65880000e-02 ... 1.95223419e+02
 1.95261719e+02            nan]
[24.042    24.047001 24.049    ... 40.004002 40.             nan]
[-0.329451 -0.325292 -0.325214 ... -0.244339 -0.244273       nan]
In [22]:
df53_2, l_time53_2, l_temp53_2, l_dsc53_2 = exp_dsc("20221207_#53_Ag2S0.4Te0.6_400degC_new.csv")

[0.00000000e+00 3.83120000e-02 7.66250000e-02 ... 1.91906723e+02
 1.91945023e+02            nan]
[22.016001 22.02     22.02     ... 40.006001 40.             nan]
[-0.265787 -0.264714 -0.264688 ... -0.397287 -0.397327       nan]
In [23]:
df54_1, l_time54_1, l_temp54_1, l_dsc54_1 = exp_dsc("20221017_#54_Ag2S0.6Te0.4_400degC.csv")

[0.00000000e+00 3.66310000e-02 7.32620000e-02 ... 1.88210083e+02
 1.88246719e+02            nan]
[32.515999 32.400002 32.389999 ... 40.007    40.             nan]
[-0.286751 -0.293329 -0.293735 ... -0.291485 -0.291211       nan]
In [24]:
df57_1_1, l_time57_1_1, l_temp57_1_1, l_dsc57_1_1 = exp_dsc("20221117_#57-1_Ag2S0.6Te0.4_1st.csv")

[0.00000000e+00 3.83280000e-02 7.66560000e-02 ... 1.93556686e+02
 1.93595016e+02            nan]
[22.816    22.82     22.820999 ... 40.015999 40.000999       nan]
[ 0.570769  0.562425  0.56181  ... -0.353659 -0.35409        nan]
In [25]:
df57_1_2, l_time57_1_2, l_temp57_1_2, l_dsc57_1_2 = exp_dsc("20221118_#57-1_Ag2S0.6Te0.4_2nd.csv")

[0.00000000e+00 3.83240000e-02 7.66480000e-02 ... 1.93573364e+02
 1.93611679e+02            nan]
[22.665001 22.669001 22.671    ... 40.012001 40.001999       nan]
[-0.264348 -0.26338  -0.262753 ... -0.313051 -0.313182       nan]
In [26]:
df57_1_3, l_time57_1_3, l_temp57_1_3, l_dsc57_1_3 = exp_dsc("20220217_#57-1_Ag2S0.6Te0.4_annealed.csv")

[0.00000000e+00 1.99970000e-02 3.99940000e-02 ... 1.04743347e+02
 1.04763344e+02            nan]
[19.604 19.608 19.608 ... 30.003 30.       nan]
[ 0.49791   0.497531  0.497531 ... -0.269226 -0.269357       nan]
In [27]:
df57_2, l_time57_2, l_temp57_2, l_dsc57_2 = exp_dsc("20221117_#57-2_Ag2S0.6Te0.4_1st.csv")

[0.00000000e+00 3.83160000e-02 7.66320000e-02 ... 1.93073364e+02
 1.93111679e+02            nan]
[24.895    24.874001 24.872999 ... 40.005001 40.             nan]
[-0.292505 -0.292061 -0.292008 ... -0.348937 -0.349055       nan]
In [28]:
df58_1, l_time58_1, l_temp58_1, l_dsc58_1 = exp_dsc("20221110_#58-1_Ag2S0.6Te0.4_1st.csv")

[0.00000000e+00 3.83060000e-02 7.66110000e-02 ... 1.90723404e+02
 1.90761703e+02            nan]
[21.035    21.028999 21.028    ... 40.011002 40.001999       nan]
[ 0.554604  0.546574  0.546012 ... -0.306787 -0.306669       nan]
In [29]:
df58_2, l_time58_2, l_temp58_2, l_dsc58_2 = exp_dsc("20221110_#58-2_Ag2S0.6Te0.4_1st.csv")

[0.00000000e+00 3.66590000e-02 7.33190000e-02 ... 1.85460037e+02
 1.85496689e+02            nan]
[37.189999 37.062    37.053001 ... 40.009998 40.000999       nan]
[-0.28866  -0.295278 -0.296337 ... -0.291119 -0.291001       nan]
In [30]:
df58_2_2, l_time58_2_2, l_temp58_2_2, l_dsc58_2_2 = exp_dsc("20221111_#58-2_Ag2S0.6Te0.4_2nd.csv")

[0.00000000e+00 3.83190000e-02 7.66380000e-02 ... 1.90023376e+02
 1.90061691e+02            nan]
[20.583    20.589001 20.591    ... 40.009998 40.             nan]
[-0.266702 -0.265237 -0.26444  ... -0.370202 -0.370189       nan]
In [31]:
df59_1, l_time59_1, l_temp59_1, l_dsc59_1 = exp_dsc("20221123_#59_Ag2S0.6Te0.4_1st.csv")

[0.00000000e+00 3.83060000e-02 7.66110000e-02 ... 1.90723404e+02
 1.90761703e+02            nan]
[20.958    20.965    20.968    ... 40.007999 40.             nan]
[-0.231941 -0.229586 -0.229207 ... -0.349892 -0.349774       nan]
In [32]:
df64_1, l_time64_1, l_temp64_1, l_dsc64_1 = exp_dsc("20230822_#64_Ag2S0.7Te0.3_1st_dsc.csv")

[0.00000000e+00 4.16420000e-02 8.32840000e-02 ... 2.09000061e+02
 2.09041718e+02            nan]
[20.454    20.459    20.461    ... 30.006001 30.000999       nan]
[-0.458977 -0.463881 -0.464143 ... -0.260438 -0.260608       nan]

Background¶

In [44]:
def background(l_temp,l_dsc,i1,i2):
    t1,d1 = l_temp[i1],l_dsc[i1]
    t2,d2 = l_temp[i2],l_dsc[i2]
    slope = (d2-d1)/(t2-t1)
    l_db = []
    for i,t in enumerate(l_temp):
        back = slope*(t-t1)+d1
        l_db.append(l_dsc[i]-back)
    return np.array(l_db)
In [45]:
l_db42_2 = background(l_temp42_2,l_dsc42_2,363,597)
In [46]:
l_db43_3 = background(l_temp43_3,l_dsc43_3,707,828)
In [47]:
l_db43 = background(l_temp43,l_dsc43,500,597)
In [48]:
l_db21 = background(l_temp21,l_dsc21,581,1096)
In [49]:
l_db20 = background(l_temp20,l_dsc20,649,1300)
In [50]:
l_db50_1 = background(l_temp50_1,l_dsc50_1,500,750)
In [51]:
def DSC_picture(l_temp,l_dsc,l_lab,xmin,xmax,ymin,ymax,filename):
    dT = 273.15
    l_color = ["black","red", "blue","green","cyan","orange"]
    l_mrks = ["o","^", "*","d",">"]
    plt.clf() # initialization
    plt.figure(figsize=(8, 6))
    plt.axis([xmin,xmax,ymin,ymax])
    plt.tick_params(direction="in")
    plt.rcParams["font.family"] = "times new roman"
    plt.rcParams['mathtext.fontset'] = 'cm'
    plt.rcParams["axes.linewidth"] = 1.5
    plt.tick_params(direction="in")
    plt.rcParams["xtick.major.size"] = 10
    plt.rcParams["xtick.major.width"] = 1.5
    plt.rcParams["ytick.major.size"] = 10
    plt.rcParams["ytick.major.width"] = 1.5
    plt.rcParams["font.size"] = 25
    plt.tick_params(labelleft=False)
    #plt.ylabel("Heat flow, mW")
    plt.ylabel("Heat flow (arb. units)")
    plt.xlabel("Temperature (K)")
    l_st = ["--","-",":","-."]
    for i in range(len(l_dsc)):
        plt.plot(l_temp[i]+dT, l_dsc[i],color=l_color[i],marker="None",label=l_lab[i],linestyle="-")
    #plt.legend(loc='upper left')
    plt.savefig(filename,bbox_inches="tight")
    plt.show()
In [52]:
def DSC_picture_shifted(l_temp,l_dsc,l_ad,xmin,xmax,ymin,ymax,filename):
    dT = 273.15
    l_color = ["black","red","blue","green","cyan","purple","orange","brown"]
    #l_color = ["green","purple"]
    l_mrks = ["o","^", "*","d",">"]
    plt.clf() # initialization
    plt.figure(figsize=(8, 6))
    plt.axis([xmin,xmax,ymin,ymax])
    plt.tick_params(direction="in")
    plt.rcParams["font.family"] = "times new roman"
    plt.rcParams['mathtext.fontset'] = 'cm'
    plt.rcParams["axes.linewidth"] = 1.5
    plt.tick_params(direction="in")
    plt.rcParams["xtick.major.size"] = 10
    plt.rcParams["xtick.major.width"] = 1.5
    plt.rcParams["ytick.major.size"] = 10
    plt.rcParams["ytick.major.width"] = 1.5
    plt.rcParams["font.size"] = 25
    plt.tick_params(labelleft=False)
    #plt.ylabel("Heat flow, mW")
    plt.ylabel("Heat flow (arb. units)")
    plt.xlabel("Temperature (K)")
    for i in range(len(l_dsc)):
        #plt.plot(l_temp[i]+dT, l_dsc[i]+l_ad,color=l_color[i],marker="None",label=l_lab[i],linestyle=l_st[i])
        plt.plot(l_temp[i]+dT, l_dsc[i]-l_ad[i],color=l_color[i],marker="None",linestyle="-")
    #plt.legend(loc='upper right')
    plt.savefig(filename,bbox_inches="tight")
    plt.show()
In [53]:
N1,N2,N3,N4,N5,N6=1500,2000,2000,2000,2000,2000
l_temp = [l_temp2[:N2],l_temp42_2[:N3]]#,l_temp17_2[:N1],l_temp50_1[:N4],l_temp53_2[:N5],l_temp21[:N6]]
l_dsc =  [l_dsc2[:N2],l_db42_2[:N3]]#,l_dsc17_2[:N1],l_dsc50_1[:N4],l_dsc53_2[:N5],l_db21[:N6]]
l_ad = [1,2,3]#,3.5,4.5,6]
DSC_picture_shifted(l_temp,l_dsc,l_ad,300,600,-2.5,0,"DSC_analysis_picture_20230314_Ag2S1-xTex_x=0.2-0.4.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [42]:
N1,N2,N3,N4,N5,N6=1500,2000,2000,2000,2000,2000
l_temp = [l_temp50_1[:N4],l_temp53_2[:N5]]#,l_temp21[:N6]]
l_dsc =  [l_dsc50_1[:N4],l_dsc53_2[:N5]]#,l_db21[:N6]]
l_ad = [1,2,3.5]#,3.5,4.5,6]
DSC_picture_shifted(l_temp,l_dsc,l_ad,300,600,-3.5,-1,"DSC_analysis_picture_20230314_Ag2S1-xTex_x=0.5-0.7.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [43]:
N1,N2,N3,N4,N5,N6=1500,2000,2000,2000,2000,2000
l_temp = [l_temp17_2[:N1],l_temp2[:N2],l_temp42_2[:N3],l_temp50_1[:N4],l_temp53_2[:N5],l_temp21[:N6]]
l_dsc =  [l_dsc17_2[:N1],l_dsc2[:N2],l_db42_2[:N3],l_dsc50_1[:N4],l_dsc53_2[:N5],l_db21[:N6]]
l_ad = [1,2,3,3.5,4.5,6.5]
DSC_picture_shifted(l_temp,l_dsc,l_ad,300,600,-8,0.5,"DSC_analysis_picture_20230313_Ag2S1-xTex_x=0.2-0.7.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [44]:
N1,N2,N3,N4,N5,N6=1500,2000,2000,2000,2000,2000
l_temp = [l_temp42_2[:N3]]
l_dsc =  [l_db42_2[:N3]]
l_ad = [1]
DSC_picture_shifted(l_temp,l_dsc,l_ad,350,600,-2,0.,"DSC_analysis_picture_20230313_Ag2S1-xTex_x=0.2-0.7.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [45]:
N1,N2,N3,N4,N5,N6=1500,2000,2000,2000,2000,2000
l_temp = [l_temp42[:N3]]
l_dsc =  [l_dsc42[:N3]]
l_ad = [1]
DSC_picture_shifted(l_temp,l_dsc,l_ad,350,650,-2,0.,"DSC_analysis_picture_20230313_Ag2S1-xTex_x=0.2-0.7.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [46]:
N=2000
l_temp = [l_temp2[:N]]#,l_temp20[N:]]
l_dsc =  [l_dsc2[:N]]#,l_dsc20[N:]]
l_lab = ["Heating"]#,"Cooling"]
DSC_picture(l_temp,l_dsc,l_lab,350,650,-0.,1,"DSC_analysis_picture_20230220_Ag2S1-xTex_x=0.3_heating.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [47]:
N=2000
l_temp = [l_temp51_1[:N]]#,l_temp57_1_3[N:]]
l_dsc =  [l_dsc51_1[:N]]#,l_dsc57_1_3[N:]]
l_lab = ["Ag2S0.7Te0.3"]#,"Cooling"]
DSC_picture(l_temp,l_dsc,l_lab,350,550,-2,2,"test.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [48]:
N=1750
l_temp = [l_temp17_2[:N]]#,l_temp20[N:]]
l_dsc =  [l_dsc17_2[:N]]#,l_dsc20[N:]]
l_lab = ["Heating"]#,"Cooling"]
DSC_picture(l_temp,l_dsc,l_lab,320,700,-0.5,2,"DSC_analysis_picture_20230220_Ag2S1-xTex_x=0.2_heating.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [49]:
N=2000
l_temp = [l_temp21[:N]]#,l_temp21[N:]]
l_dsc =  [l_dsc21[:N]]#,l_dsc21[N:]]
l_lab = ["Heating","Cooling"]
DSC_picture(l_temp,l_dsc,l_lab,300,600,-1,1,"DSC_analysis_picture_20230216_Ag2S1-xTex_x=0.3_heating.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [50]:
N=2000
l_temp = [l_temp21[:N]]#,l_temp21[N:]]
l_dsc =  [l_db21[:N]]#,l_dsc21[N:]]
l_lab = ["Heating","Cooling"]
DSC_picture(l_temp,l_dsc,l_lab,350,600,-1,1,"DSC_analysis_picture_20230216_Ag2S1-xTex_x=0.3_heating.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [51]:
N=2000
l_temp = [l_temp20[:N]]#,l_temp21[N:]]
l_dsc =  [l_dsc20[:N]]#,l_dsc21[N:]]
l_lab = ["Heating","Cooling"]
DSC_picture(l_temp,l_dsc,l_lab,300,600,-2,0,"DSC_analysis_picture_20230314_Ag2S1-xTex_x=0.2_heating.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [52]:
N=2000
l_temp = [l_temp20[:N]]#,l_temp21[N:]]
l_dsc =  [l_db20[:N]]#,l_dsc21[N:]]
l_lab = ["Heating","Cooling"]
DSC_picture(l_temp,l_dsc,l_lab,300,700,-2,2,"DSC_analysis_picture_20230314_Ag2S1-xTex_x=0.2_heating.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [53]:
N=2000
l_temp = [l_temp42_2[:N]]#,l_temp42_2[N:]]
l_dsc =  [l_dsc42_2[:N]]#,l_dsc42_2[N:]]
l_lab = ["Heating","Cooling"]
DSC_picture(l_temp,l_dsc,l_lab,400,600,-0.25,4.,"DSC_analysis_picture_20230220_Ag2S1-xTex_x=0.4.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [135]:
N=2000
l_temp = [l_temp43_1[:N]]#,l_temp43[N:]]
l_dsc =  [l_db43_3[:N]]#,l_dsc43[N:]]
l_lab = ["heating","cooling"]
DSC_picture(l_temp,l_dsc,l_lab,320,500,-0.5,0.5,"DSC_analysis_picture_20230220_Ag2S1-xTex_x=0.5.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [140]:
N=2000
l_temp = [l_temp43[:N]]#,l_temp43[N:]]
l_dsc =  [l_dsc43[:N]]#,l_dsc43[N:]]
l_lab = ["heating","cooling"]
DSC_picture(l_temp,l_dsc,l_lab,320,500,-0.5,0.5,"DSC_analysis_picture_20230220_Ag2S1-xTex_x=0.5.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [142]:
N=2000
l_temp = [l_temp43[:N]]#,l_temp43[N:]]
l_dsc =  [l_db43[:N]]#,l_dsc43[N:]]
l_lab = ["heating","cooling"]
DSC_picture(l_temp,l_dsc,l_lab,320,500,-0.5,0.5,"DSC_analysis_picture_20230220_Ag2S1-xTex_x=0.5.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [62]:
N=2000
l_temp = [l_temp57_1_1[:N]]#,l_temp43[N:]]
l_dsc =  [l_dsc57_1_1[:N]]#,l_dsc43[N:]]
l_lab = ["heating","cooling"]
DSC_picture(l_temp,l_dsc,l_lab,300,600,-1,1,"test.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [54]:
N=1800
l_temp = [l_temp50_1[:N],l_temp50_1[N:]]
l_dsc =  [l_dsc50_1[:N],l_dsc50_1[N:]]
l_lab = ["Heating","Cooling"]
DSC_picture(l_temp,l_dsc,l_lab,340,650,-0.9,1,"DSC_analysis_picture_20221205_Ag2S1-xTex_x=0.5_#50_up-down.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [61]:
N=1800
l_temp = [l_temp50_1[:N]]
l_dsc =  [l_dsc50_1[:N]]
l_lab = ["Heating","Cooling"]
DSC_picture(l_temp,l_dsc,l_lab,300,650,-1,-0,"DSC_analysis_picture_20230309_Ag2S1-xTex_x=0.5_precipitation.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [60]:
N=1800
l_temp = [l_temp50_1[:N]]
l_dsc =  [l_db50_1[:N]]
l_lab = ["Heating","Cooling"]
DSC_picture(l_temp,l_dsc,l_lab,300,600,-0.5,0.5,"DSC_analysis_picture_20230309_Ag2S1-xTex_x=0.5_precipitation.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [59]:
N=1900
l_temp = [l_temp53_1[:N]]#,l_temp53_2[N:]]
l_dsc =  [l_dsc53_1[:N]]#,l_dsc53_2[N:]]
l_lab = ["Heating","Cooling"]
DSC_picture(l_temp,l_dsc,l_lab,300,600,-1.5,0.5,"DSC_analysis_picture_20230309_Ag2S1-xTex_x=0.6_precipitation.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [29]:
N=2000
l_temp = [l_temp57_1[:N],l_temp57_1[N:]]
l_dsc =  [l_dsc57_1[:N],l_dsc57_1[N:]]
l_lab = ["heating","cooling"]
DSC_picture(l_temp,l_dsc,l_lab,350,650,-1,2,"DSC_analysis_picture_20221118_Ag2S0.6Te0.4_#57-1_up_down.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [30]:
N=2000
l_temp = [l_temp57_1_1[:N],l_temp57_1_1[N:]]
l_dsc =  [l_dsc57_1_1[:N],l_dsc57_1_1[N:]]
l_lab = ["heating","cooling"]
DSC_picture(l_temp,l_dsc,l_lab,350,650,-1,2,"DSC_analysis_picture_20221118_Ag2S0.6Te0.4_#57-1_2nd_up_down.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [31]:
N=2000
l_temp = [l_temp57_2[:N],l_temp57_2[N:]]
l_dsc =  [l_dsc57_2[:N],l_dsc57_2[N:]]
l_lab = ["heating","cooling"]
DSC_picture(l_temp,l_dsc,l_lab,350,650,-1,2,"DSC_analysis_picture_20221118_Ag2S0.6Te0.4_#57-2_up_down.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [32]:
N=2000
l_temp = [l_temp58_1[:N],l_temp58_1[N:]]
l_dsc =  [l_dsc58_1[:N],l_dsc58_1[N:]]
l_lab = ["heating","cooling"]
DSC_picture(l_temp,l_dsc,l_lab,350,650,-1,2,"DSC_analysis_picture_20221111_Ag2S0.6Te0.4_#58-1_up_down.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [33]:
N=2000
l_temp = [l_temp58_2_2[:N],l_temp58_2_2[N:]]
l_dsc =  [l_dsc58_2_2[:N],l_dsc58_2_2[N:]]
l_lab = ["heating","cooling"]
DSC_picture(l_temp,l_dsc,l_lab,350,650,-1,2,"DSC_analysis_picture_20221111_Ag2S0.6Te0.4_#58-2_2nd_up_down.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [34]:
N=2000
l_temp = [l_temp59_1[:N],l_temp59_1[N:]]
l_dsc =  [l_dsc59_1[:N],l_dsc59_1[N:]]
l_lab = ["heating","cooling"]
DSC_picture(l_temp,l_dsc,l_lab,350,650,-1,2,"DSC_analysis_picture_20221124_Ag2S0.6Te0.4_#59_up_down.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image

Paper¶

In [56]:
N=1500
l_temp = [l_temp3[:N]]#,l_temp2[N:]]
l_dsc =  [l_dsc3[:N]]#,l_dsc2[N:]]
l_lab = ["Heating"]#,"cooling"]
DSC_picture(l_temp,l_dsc,l_lab,300,525,-2.5,2,"DSC_analysis_picture_20221205_Ag2S1-xTex_x=0.3_up.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [34]:
N=1200
l_temp = [l_temp2[N:]]
l_dsc =  [l_dsc2[N:]]
l_lab = ["cooling"]
DSC_picture(l_temp,l_dsc,l_lab,350,550,-0.4,0.2,"DSC_analysis_picture_20221221_Ag2S1-xTex_x=0.3_cooling.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [54]:
N=1850
l_temp = [l_temp57_1_1[:N]]
l_dsc =  [l_dsc57_1_1[:N]]
l_lab = ["heating"]
DSC_picture(l_temp,l_dsc,l_lab,300,600,-1,1,"DSC_analysis_picture_20221118_Ag2S0.6Te0.4_#57-1_up.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [55]:
N=2000
l_temp = [l_temp42_2[:N]]#,l_temp54_1[N:]]
l_dsc =  [l_dsc42_2[:N]]#,l_dsc54_1[N:]]
l_lab = ["heating","cooling"]
DSC_picture(l_temp,l_dsc,l_lab,350,600,-0.5,1,"DSC_analysis_picture_20221017_Ag2S1-xTex_x=0.4_up.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [35]:
N=2000
l_temp = [l_temp42_2[N:]]#,l_temp54_1[N:]]
l_dsc =  [l_dsc42_2[N:]]#,l_dsc54_1[N:]]
l_lab = ["cooling"]
DSC_picture(l_temp,l_dsc,l_lab,350,600,-0.5,1,"DSC_analysis_picture_20221221_Ag2S1-xTex_x=0.4_cooling.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [63]:
N1s,N1e,N2s,N2e = 100,1800,1900,-1
l_temp = [l_temp50_1[:N1e]]#, l_temp50_1[N2s:N2e]]
l_dsc =  [l_dsc50_1[:N1e]]#,  l_dsc50_1[N2s:N2e]]
l_lab = ["up","down"]
DSC_picture(l_temp,l_dsc,l_lab,300,650,-1,-0.2,"DSC_analysis_picture_20221011_Ag2S1-xTex_x=0.5_#50_up.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [65]:
N=1900
l_temp = [l_temp53_1[:N]]#,l_temp53_1[N:]]
l_dsc =  [l_dsc53_1[:N]]#,l_dsc53_1[N:]]
l_lab = ["Heating","Cooling"]
DSC_picture(l_temp,l_dsc,l_lab,300,700,-1.5,0.,"DSC_analysis_picture_20221205_Ag2S1-xTex_x=0.6_up.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [44]:
N=1900
l_temp = [l_temp53_1[N:]]#,l_temp53_1[N:]]
l_dsc =  [l_dsc53_1[N:]]#,l_dsc53_1[N:]]
l_lab = ["Cooling"]
DSC_picture(l_temp,l_dsc,l_lab,350,650,-0.5,-0.3,"DSC_analysis_picture_20221221_Ag2S1-xTex_x=0.6_#53_cooling.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [55]:
def Phase_transition_temperature_graph(cstat,l_x,l_temp,l_color,xmin,xmax,ymin,ymax,filename):
    #l_color = ["black","red", "blue","green","cyan","orange"]
    l_mrks = ["o","^", "*","d",">"]
    plt.clf() # initialization
    plt.figure(figsize=(8, 6))
    plt.axis([xmin,xmax,ymin,ymax])
    plt.tick_params(direction="in")
    plt.rcParams["font.family"] = "Times New roman"
    plt.rcParams['mathtext.fontset'] = 'cm'
    plt.rcParams["axes.linewidth"] = 1.5
    plt.tick_params(direction="in")
    plt.rcParams["xtick.major.size"] = 10
    plt.rcParams["xtick.major.width"] = 1.5
    plt.rcParams["ytick.major.size"] = 10
    plt.rcParams["ytick.major.width"] = 1.5
    plt.rcParams["font.size"] = 23
    #plt.tick_params(labelleft=False)
    #plt.ylabel("Heat flow, mW")
    plt.xlabel("$x$ in Ag$_{2}$S$_{1-x}$Te$_{x}$")
    plt.ylabel(r"Temperature (K)")
    l_st = ["--","-",":","-."]
    for i in range(len(l_temp)):
        if i !=4 and i !=5 and i != 10:
            plt.plot(l_x[i], l_temp[i],color=l_color[i],marker="o")#,label=l_lab[i])
        elif i == 10:
            plt.plot(l_x[i], l_temp[i],color=l_color[i],marker="o",linestyle = "--")#,label=l_lab[i])
    if cstat == 1:
        for i in range(len(l_temp)):
            plt.fill_between(l_x[i],l_temp[i],0,facecolor='blue', alpha=0.2)
            plt.fill_between(l_x[i],ymax,l_temp[i],facecolor='red', alpha=0.2)
            plt.fill_between([0.12,0.2,0.3,0.4,0.5,0.6,0.7,0.8],ymax,ymin,facecolor='red', alpha=0.1)
    elif cstat == 2:
        plt.fill_between(l_x[0],l_temp[0],0,facecolor='purple', alpha=1)
        #plt.fill_between(l_x[0][2:],ymax,l_temp[0][2:],facecolor='red', alpha=0.5)
        plt.fill_between(l_x[1],l_temp[1],l_temp[0][:3],facecolor='purple', alpha=1)
        plt.fill_between(l_x[1],ymax,l_temp[1],facecolor='red', alpha=1)
        plt.fill_between(l_x[2],l_temp[2],0,facecolor='blue', alpha=1)
        #plt.fill_between(l_x[2],ymax,l_temp[2],facecolor='gray', alpha=0.5)
        plt.fill_between(l_x[3],l_temp[3],0,facecolor='blue', alpha=1)
        plt.fill_between(l_x[3][:2],ymax,l_temp[3][:2],facecolor='red', alpha=1)
        #plt.fill_between(l_x[5],l_temp[5],0,facecolor='cyan', alpha=0.5)
        #plt.fill_between(l_x[5],ymax,l_temp[5],facecolor='brown', alpha=0.5)
        plt.fill_between(l_x[6],l_temp[6],0,facecolor='cyan', alpha=1)
        plt.fill_between(l_x[7],ymax,l_temp[7],facecolor='black', alpha=1)
        plt.fill_between(l_x[8],ymax,l_temp[8],facecolor='brown', alpha=1)
        plt.fill_between([0.7,0.8],ymax,ymin,facecolor='gray', alpha=1)
        plt.fill_between(l_x[9],l_temp[9],0,facecolor="yellow", alpha=1)
        plt.fill_between(l_x[9],ymax,l_temp[9],facecolor='brown', alpha=1)
        plt.fill_between(l_x[10],l_temp[10],0,facecolor='green', alpha=1)
        plt.fill_between(l_x[10],ymax,l_temp[10],facecolor='brown', alpha=1)
        plt.fill_between(l_x[11],l_temp[11],l_temp[2],facecolor='purple', alpha=1)
        plt.fill_between(l_x[11],ymax,l_temp[11],facecolor='red', alpha=1)
    #plt.legend(loc='upper left')
    plt.savefig(filename,bbox_inches="tight")
    plt.show()
In [57]:
l_x0 = [0.3,0.4,0.5,0.6,0.7]
l_htp_t0 = [384,438,537,557,560]
l_x1 = [0.3,0.4,0.5]
l_htp_t1 = [498,526,537]
l_x2 = [0.2,0.3]
l_htp_t2 = [373,384]
l_x3 = [0.,0.2]
l_htp_t3 = [450,373]
l_x4 = [0,0.03,0.06,0.09,0.12] #ref
l_htp_t4 = [450,429,374.5,346.3,310.7] #ref
l_x5 = [0.8,0.9,1.0] # ref
l_htp_t5 = [320.6,378.7,429] #ref
l_x6 = [0.3,0.4,0.5,0.6,0.7]
l_htp_t6 =  [384,438,537,557,560]
l_htp_t7 = [537,557]
l_x7 = [0.5,0.6]
l_htp_t8 = [557,560]
l_x8 = [0.6,0.7]
l_x9 = [0.8,1.0] 
l_htp_t9 = [335,429]
l_x10 = [0.7,0.8]
l_htp_t10 = [350,335]
l_x11 = [0.2,0.3]
l_htp_t11 = [440,498]
l_temp = [l_htp_t0,l_htp_t1,l_htp_t2,l_htp_t3,l_htp_t4,l_htp_t5,\
          l_htp_t6,l_htp_t7,l_htp_t8,l_htp_t9,l_htp_t10,l_htp_t11]#,l_htp_tp]#,l_htp_ta]
l_x = [l_x0,l_x1,l_x2,l_x3,l_x4,l_x5,l_x6,l_x7,l_x8,l_x9,l_x10,l_x11]
#l_color = ["red","red","red","red","gray", "gray","red","red","red","red","red","red"]
l_color = ["black"]*12
Phase_transition_temperature_graph(2,l_x,l_temp,l_color,-0.01,1.01,300,600,"20230117_phase-transition_temperatures.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [44]:
#l_x1 = [0.2,0.3,0.4,0.5,0.6,0.7]
#l_htp_t1 = [410,384,438,537,557,560]
l_x0 = [0.3,0.4,0.5,0.6]
l_htp_t0 = [384,438,537,557]
l_x2 = [0,0.03,0.06,0.09,0.12]
l_htp_t2 = [450,429,374.5,346.3,310.7]
l_x3 = [0.8,0.9,1.0]
l_htp_t3 = [320.6,378.7,429]
l_temp = [l_htp_t2,l_htp_t3]#,l_htp_tp]#,l_htp_ta]
l_x = [l_x2,l_x3]
l_color = ["black","black"]
Phase_transition_temperature_graph(1,l_x,l_temp,l_color,-0.01,1.01,300,600,"20230117_phase-transition_temperatures.jpg")

<Figure size 432x288 with 0 Axes>
No description has been provided for this image
In [ ]: