a
ɝh;B @ s d dl Z d dlZd dlZe je jeZe jeZeejvrNej e d dl
m Z d dl
mZ d dlmZ d dlZd dlZd dlZd dlmZ d dlZdd Z dd Z dd Z d d
Z dS ) N) Structurec C s d}zd| d j d d| d j d d| d j d d | d
j d d | d j d d
| d j d d| d j d d| d j d d| d| d j d d}W n8 ty } z td| d}W Y d }~n
d }~0 0 |S )Nz@$$\kappa_L=A\frac{M V^{\frac{1}{3}} \theta_a^3}{\gamma^2 T n} $$zA
The number of atoms of the crystal structure is Number of Atomsr z=.
The volume of the crystal structure is Volume (Å3)uO Å$^3$.
The total atomic mass of the crystal structure is the total atomic mass (amu)zG amu.
The Bulk modulus of the crystal structure is Bulk modulus (GPa)zH GPa.
The Shear modulus of the crystal structure is Shear modulus (GPa)zI GPa.
The sound velocity of the crystal structure is Speed of sound (m s-1)zV m/s.
The Acoustic Debye Temperature of the crystal structure is Acoustic Debye Temperature (K)uM K.
The Grüneisen parameter of the crystal structure is Grüneisen parameterzJ.
The formula of the lattice thermal conductivity is: zE.
The calculated lattice thermal conductivity is Kappa_Slack (W m-1 K-1)u W/(m·K).
z!Error in display_results_kappap: zError displaying results)iloc Exceptionstwrite)dfformulatemplatee r +/var/www/html/AI4Kappa/Pages/CustomKappa.pydisplay_results_kappap s6
r c C st d}d| d j d d| d j d d| d j d d | d
j d d | d j d d
| d| d j d d}|S )NzG$$\kappa_L=\frac{G \upsilon_s V^{\frac{1}{3}}}{N T} \cdot e^{-\gamma}$$z4
The volume of the crystal structure is r r uF Å$^3$.
The Bulk modulus of the crystal structure is r zD GPa.
The Shear modulus of the crystal structure is r zE GPa.
The sound velocity of the crystal structure is r zI m/s.
The Grunisen parameter of the crystal structure is r
zF.
The formula of the lattice thermal conductivity is: zA.
The calculated lattice thermal conductivity is Kappa_cal (W m-1 K-1)u W/(m·K).
)r )r r r r r r display_results_ai4kappa. s" r c C s | dkrg d}ng d}|S )NKappaP)
r Density (g cm-3)r r r r -Sound velocity of the transverse wave (m s-1),Sound velocity of the longitude wave (m s-1)r
Poisson ratior
r r )
r r r r r r r r r r r
r r r )methodlsr r r display_columns: s
r c 2 C s6 t d tjd} t jj}tj| d}t jdddgdd}t jj rt
tj|d }t |d
krxt d d S t
d t d
i }|D ]}tj|}t
d| d t d\}} }
zFtttj||} | d} | d u st| ttfs d} W n> tyL }
z$td| d|
d} W Y d }
~
n
d }
~
0 0 |( t jdd| d dd}W d n1 s0 Y | ( t jdd| d dd}W d n1 s0 Y d }|d urh|d urh| d urhz |d }|d }| d }|d | d | d }|| d }|d krNt d | d! d }n|| }t |d"d" d"t |d" d" }t d"d| d#k rt d | d$ d }nTdd%| d"d"d| }d | k rd&k sn t d | d'|d(d) d }W nl t y$ t d | d* d }Y nD tyf }
z*t d | d+t |
d }W Y d }
~
n
d }
~
0 0 |
Z | d,| }d-| t jvr|t jd-| < |t jd-| krd.| t jv rt jd.| = |t jd-| < t j!d/d0| t jd.| d1d2d}zL|" r4t|}|t jd.| < n$|}d.| t jv rXt jd.| = W n" t#y| t d3 d }Y n0 |d urt
d4|d( |" rt
d5 n
t
d6 W d n1 s0 Y i }|d ur||d< |d ur||d< |d ur||d/< |r|||< |||< qt $d7} | r i } i } |% D ]\}}g } g } d|vrp| &d8 n|d d kr| &d9 d|vr| &d: n|d d kr| &d; d/|vr| &d/ n|d/ d kr| &d< | r| | |< | rJ| | |< qJ| s| rt d= | rZt
d> | % D ]&\}}t
d?| dd@| q2| rt
dA | % D ]&\}}!t
d?| dd@|! qrd S z"t'|}"t(j)|"j*g dBdC}#dDD ] }$|$|"jv r|"|$ |#|$< q|% D ]H\}}tj+|d }%|%|#j*v r|% D ]\}&}'|'|#j,|%|&f< qq|#d/ }(t-.|#})|dkr|t-/|)|(}*t-0|*d%|(}+t-1|+},t2}-nt-/|)|(}*t-3|*},t4}-t
d t dE t
dF t 5|,j,d d t6|f t
d t dG |7 D ]}tj+|d }%t 8dH| t9tj||}.t j
|.ddI z&|,j,|%g }/|-|/}0t j:|0ddI W n@ ty }
z&t dJ| dt |
W Y d }
~
n
d }
~
0 0 W d n1 s0 Y qW n< ty }
z"t dKt |
W Y d }
~
d S d }
~
0 0 t;| t<| n
t =dL dM}1t j:|1ddI d S )NNzCustom Kappa Calculator.modelzSelect calculation method:r PINKT)
horizontalz*.cif z3Maximum 5 files allowed. Please upload fewer files.z---zCustom Parameters Inputz**Parameters for z:** r g @zError reading density for z: r Zbulk_zEnter value...)keyvalue placeholderr Zshear_g eAi g ?r z Warning for uW : Cannot calculate Grüneisen parameter - Sound velocity of the transverse wave is zero g|=u? : Cannot calculate Grüneisen parameter - Invalid Poisson ratio
u# : Calculated Grüneisen parameter (z.4fz#) is out of reasonable range (0-10)u@ : Cannot calculate Grüneisen parameter - Division by zero erroru+ : Error calculating Grüneisen parameter - _Zprev_bulk_shear_Z
user_grun_r
Z grun_input_ z&Enter value or leave empty for defaultzPlease enter a valid numberzCurrent value: z(User defined)z(Default calculated)Z Calculatez Bulk modulusz Bulk modulus must be positivez
Shear modulusz Shear modulus must be positiveu% Grüneisen parameter must be positivez$Please correct the following errors:zMissing parameters:z- z, zInvalid parameters:)r r r r r r r
)indexcolumns)r r r r ResultszCombined Results:z Crystal Structure InformationzStructure details for )unsafe_allow_htmlz Error displaying results for z&An error occurred during calculation: z9Please upload CIF files (maximum 5) in the sidebar first.aK We strive to have clear documentation and examples to help everyone with using our calculator.
We will happily fix issues in the documentation and examples should you find any,
however, we will not be able to offer extensive user support and training, except for our collaborators.
)>r titleospathabspath
session_state
root_dir_pathjoinradiouploaded_filesgloblenerrorr subheaderbasenamer1 foget_crystalline_datar from_fileget
isinstanceintfloatr
print number_inputwarningpowabsZeroDivisionErrorstr
text_inputstrip
ValueErrorbuttonitemsappendget_dir_crystalline_datapd DataFramer0 splitextloccalk cal_Debye_T cal_gammacal_A cal_K_Slackr by_MTPr dataframer keysexpanderget_crystalline_contentmarkdown del_cif_file
del_temp_fileinfo)2 sour_pathr9
model_pathr Z cif_filesZ file_paramscif_file file_namecol1col2Zcol3Z structure_dfdensityr ZbulkZshearZdefault_gruneisenZbulk_paZshear_paZ
density_siVlVtapoissonZcurrent_bulk_shearZ
user_inputZ gruneisenparams calculateZmissing_paramsZinvalid_paramsmissinginvaliderrors
all_cry_df
whole_info_dfcol base_nameparamr( custom_gammaDebye_dfgamma_dfA_dffinal_dfZ display_func cry_contentZ file_resultsr declarationr r r appG sx
&&
*
X
r ) r5 r= sysr6 dirnamer7 __file__ current_dir
parent_dirrU streamlit_scripts.file_opfile_oprB streamlit_scripts.chang_model chang_modelcm streamlit_scripts.calculate_K calculate_Kr[ predict streamlitr pandasrW
pymatgen.corer numpynpr r r r r r r r s$