calphy 1.3.11__py3-none-any.whl → 1.3.13__py3-none-any.whl

calphy/__init__.py

@@ -4,7 +4,7 @@ from calphy.solid import Solid
 from calphy.alchemy import Alchemy
 from calphy.routines import MeltingTemp
 
-__version__ = "1.3.11"
+__version__ = "1.3.13"
 
 def addtest(a,b):
     return a+b

calphy/input.py

@@ -40,7 +40,7 @@ from pyscal3.core import structure_dict, element_dict, _make_crystal
 from ase.io import read, write
 import shutil
 
-__version__ = "1.3.11"
+__version__ = "1.3.13"
 
 def _check_equal(val):
     if not (val[0]==val[1]==val[2]):

calphy/phase.py

@@ -598,8 +598,12 @@ class Phase:
         
         #now we can check if it converted
         file = os.path.join(self.simfolder, "avg.dat")
-        lx, ly, lz, ipress = np.loadtxt(file, usecols=(1, 2, 3, 4), unpack=True)
-        lxpc = ipress
+        #we have to clean the data, so as just the last block is selected
+        lx, ly, lz, lxpc = np.loadtxt(file, usecols=(1, 2, 3, 4), unpack=True)
+        lx = lx[-ncount+1:]
+        ly = ly[-ncount+1:]
+        lz = lx[-ncount+1:]
+        lxpc = lxpc[-ncount+1:]
         mean = np.mean(lxpc)
         std = np.std(lxpc)
         volatom = np.mean((lx*ly*lz)/self.natoms)
@@ -612,16 +616,20 @@ class Phase:
             ncount = int(self.calc.md.n_small_steps)//int(self.calc.md.n_every_steps*self.calc.md.n_repeat_steps)
 
         file = os.path.join(self.simfolder, "avg.dat")
-        lx, ly, lz, ipress = np.loadtxt(file, usecols=(1, 2, 3, 4), unpack=True)        
-        lxpc = ipress
+        lx, ly, lz, lxpc = np.loadtxt(file, usecols=(1, 2, 3, 4), unpack=True)        
+        lx = lx[-ncount+1:]
+        ly = ly[-ncount+1:]
+        lz = lx[-ncount+1:]
+        lxpc = lxpc[-ncount+1:]
+
         mean = np.mean(lxpc)
         std = np.std(lxpc)
         volatom = np.mean((lx*ly*lz)/self.natoms)
 
         self.calc._pressure = mean
-        self.lx = np.round(np.mean(lx[-ncount+1:]), decimals=3)
-        self.ly = np.round(np.mean(ly[-ncount+1:]), decimals=3)
-        self.lz = np.round(np.mean(lz[-ncount+1:]), decimals=3)
+        self.lx = np.round(np.mean(lx), decimals=3)
+        self.ly = np.round(np.mean(ly), decimals=3)
+        self.lz = np.round(np.mean(lz), decimals=3)
         self.volatom = volatom
         self.vol = self.lx*self.ly*self.lz
         self.rho = self.natoms/(self.lx*self.ly*self.lz)

calphy/postprocessing.py

@@ -1,6 +1,8 @@
 import os
 import numpy as np
 import yaml
+import matplotlib.pyplot as plt
+import warnings
 
 def read_report(folder):
     """
@@ -145,4 +147,83 @@ def gather_results(mainfolder):
                 datadict['error_code'][-1] = _extract_error(errfile)
 
     df = pd.DataFrame(data=datadict)
-    return df
+    return df
+
+def find_transition_temperature(folder1, folder2, fit_order=4, plot=True):
+    """
+    Find transition temperature where free energy of two phases are equal.
+
+    Parameters
+    ----------
+    folder1: string
+        directory with temperature scale calculation
+
+    folder2: string
+        directory with temperature scale calculation
+
+    fit_order: int, optional
+        default 4. Order for polynomial fit of temperature vs free energy
+    
+    plot: bool, optional
+        default True. Plot the results.
+    """
+    file1 = os.path.join(folder1, 'temperature_sweep.dat')
+    file2 = os.path.join(folder2, 'temperature_sweep.dat')
+    if not os.path.exists(file1):
+        raise FileNotFoundError(f'{file1} does not exist')
+    if not os.path.exists(file2):
+        raise FileNotFoundError(f'{file2} does not exist')
+
+    t1, f1 = np.loadtxt(file1, unpack=True, usecols=(0,1))
+    t2, f2 = np.loadtxt(file2, unpack=True, usecols=(0,1))
+
+    #do some fitting to determine temps
+    t1min = np.min(t1)
+    t2min = np.min(t2)
+    t1max = np.max(t1)
+    t2max = np.max(t2)
+
+    tmin = np.min([t1min, t2min])
+    tmax = np.max([t1max, t2max])
+
+    #warn about extrapolation
+    if not t1min == t2min:
+        warnings.warn(f'free energy is being extrapolated!')
+    if not t1max == t2max:
+        warnings.warn(f'free energy is being extrapolated!')
+
+    #now fit
+    f1fit = np.polyfit(t1, f1, fit_order)
+    f2fit = np.polyfit(t2, f2, fit_order)
+
+    #reevaluate over the new range
+    fit_t = np.arange(tmin, tmax+1, 1)
+    fit_f1 = np.polyval(f1fit, fit_t)
+    fit_f2 = np.polyval(f2fit, fit_t)
+
+    #now evaluate the intersection temp
+    arg = np.argsort(np.abs(fit_f1-fit_f2))[0]
+    transition_temp = fit_t[arg]
+
+    #warn if the temperature is shady
+    if np.abs(transition_temp-tmin) < 1E-3:
+        warnings.warn('It is likely there is no intersection of free energies')
+    elif np.abs(transition_temp-tmax) < 1E-3:
+        warnings.warn('It is likely there is no intersection of free energies')
+
+    #plot
+    if plot:
+        c1lo = '#ef9a9a'
+        c1hi = '#b71c1c'
+        c2lo = '#90caf9'
+        c2hi = '#0d47a1'
+
+        plt.plot(fit_t, fit_f1, color=c1lo, label=f'{folder1} fit')
+        plt.plot(fit_t, fit_f2, color=c2lo, label=f'{folder2} fit')
+        plt.plot(t1, f1, color=c1hi, label=folder1, ls='dashed')
+        plt.plot(t2, f2, color=c2hi, label=folder2, ls='dashed')
+        plt.axvline(transition_temp, ls='dashed', c='#37474f')
+        plt.ylabel('Free energy (eV/atom)')
+        plt.xlabel('Temperature (K)')
+        plt.legend(frameon=False)
+    return transition_temp

{calphy-1.3.11.dist-info → calphy-1.3.13.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: calphy
-Version: 1.3.11
+Version: 1.3.13
 Summary: free energy calculation for python
 Home-page: https://github.com/ICAMS/calphy
 Author: Sarath Menon, Yury Lysogorskiy, Ralf Drautz

{calphy-1.3.11.dist-info → calphy-1.3.13.dist-info}/RECORD

@@ -1,25 +1,25 @@
-calphy/__init__.py,sha256=yzftsRGAExKS7j5P49gh0fpg-c43MIsUCXKiUgGnCCg,234
+calphy/__init__.py,sha256=O46MqskwAlD2DmoXB8prjLeuCUBlzHuKsv-SaqjTohY,234
 calphy/alchemy.py,sha256=epv_vJoAMVbj-S1TuPSbBGr4w_a9qEr4EIxKlgwdlSo,12893
 calphy/clitools.py,sha256=oDqaw0s-LJ7tAdW_Njk2SFljVx6K34Rs8sdMz48SNSI,4125
 calphy/composition_transformation.py,sha256=Hh240-iVGC8ATlJ3nQK757qVZIBrE2aw7dsF46mi3oo,15353
 calphy/errors.py,sha256=KN47RWTLbg1H_NZMrhCiJCbqjqJScJ1pgQAuzj1-l84,1268
 calphy/helpers.py,sha256=nj8g53H6ScohQtyV7Enzms8Rv_89PRrDDY1IHjFSZdQ,8292
-calphy/input.py,sha256=R-M0C-7i1ux1IWop0UzlqAKoYIre7cxVuFLspcUwsPk,29400
+calphy/input.py,sha256=836ifhhVHm6Nd9jMWLjMUYKukk_pcZwvoGcLzQXIWuU,29400
 calphy/integrators.py,sha256=hJYmbznsgY-x-eZm7ng8gW0qsGbopERl5gGDnw5uLms,21719
 calphy/kernel.py,sha256=rd_-EfCiBhQjkxcVaoLtVJB2_qDgS-g-dQ0BZBTJQ_A,6190
 calphy/liquid.py,sha256=a5NTAjc3VsrPBQoEMGe9O1WXfv1vxGoB0xPguf0fOyM,13762
-calphy/phase.py,sha256=FOm5_eP_8XbCYG3_UJoAToD3KKxlXcbSoEfQncZXTeA,44967
+calphy/phase.py,sha256=NBkOGkG3U9Uo-dbupFWndOaaCiiSRoD8WHUpunYZty0,45190
 calphy/phase_diagram.py,sha256=2EwmT_qkT5BEP6Sx0_rm09kPP2RWh5JY4sogIBK_AhA,12992
-calphy/postprocessing.py,sha256=hFsc5Kr1FOgfjicjm23Dy19BUbZU9GxKpbCa8_U6gbU,4947
+calphy/postprocessing.py,sha256=CztgekCvQ8lZPVfZ_hMAwlyvwetjSfADjTca8aVTftY,7602
 calphy/queuekernel.py,sha256=4GMIYnjMiAPipoLNKP5noYcfeEOI_vCqm84zgokk7Xw,5321
 calphy/routines.py,sha256=W6OZEPv6HEG11EYsoIbum2WVrO3Ly36i5UWsYQ4oBdQ,17473
 calphy/scheduler.py,sha256=IN8ogDedpTbZZsdpOj-hYZI05gWoD4bN7mHOb1z77Vo,8579
 calphy/solid.py,sha256=ZU_c4LqASoLzhjXX6eQY_8k_FNO9wjmMVTiOCKNsmis,19896
 calphy/splines.py,sha256=BGwUVz_qXQxUzpUCuZo6CsELcd5JVNWzI-Ttcz22G_E,61627
 calphy/utils.py,sha256=0UpsYoxjS5N-iGs-cdm0YDMkLF8IHvKO3smXDHrj3eg,3818
-calphy-1.3.11.dist-info/LICENSE,sha256=XIHGB5RZLIhOjjoO1bPf0II-qDbjhP5Cv5HJMRE9v1g,16651
-calphy-1.3.11.dist-info/METADATA,sha256=Pw5F6YJyyaGZlJG4kPU_MMf8UoilMH16JPa4iU7TFvY,4216
-calphy-1.3.11.dist-info/WHEEL,sha256=OVMc5UfuAQiSplgO0_WdW7vXVGAt9Hdd6qtN4HotdyA,91
-calphy-1.3.11.dist-info/entry_points.txt,sha256=W9qq254koyWnAgo1jtfQP9bO5Q7sgZrzc8BMnfo3vf4,386
-calphy-1.3.11.dist-info/top_level.txt,sha256=w871dhMqPwgjjbifBWdkT9_aOnK1ek4Odrh8UnSG3PE,7
-calphy-1.3.11.dist-info/RECORD,,
+calphy-1.3.13.dist-info/LICENSE,sha256=XIHGB5RZLIhOjjoO1bPf0II-qDbjhP5Cv5HJMRE9v1g,16651
+calphy-1.3.13.dist-info/METADATA,sha256=JCQSN5W690zkuPLQ9vqmc8mSOr56K6GOq8EiqEaymWo,4216
+calphy-1.3.13.dist-info/WHEEL,sha256=PZUExdf71Ui_so67QXpySuHtCi3-J3wvF4ORK6k_S8U,91
+calphy-1.3.13.dist-info/entry_points.txt,sha256=W9qq254koyWnAgo1jtfQP9bO5Q7sgZrzc8BMnfo3vf4,386
+calphy-1.3.13.dist-info/top_level.txt,sha256=w871dhMqPwgjjbifBWdkT9_aOnK1ek4Odrh8UnSG3PE,7
+calphy-1.3.13.dist-info/RECORD,,

{calphy-1.3.11.dist-info → calphy-1.3.13.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (75.2.0)
+Generator: setuptools (75.6.0)
 Root-Is-Purelib: true
 Tag: py3-none-any