def delta_P(T, rho):
HI = hi(T, rho)
PHI = progonka(T, rho, 1, 1)
#print("Hi_n = ", HI[N], "len = ", len(HI))
return 8 / (3 * pi**4) * (2 / pi)**(1 / 3) * (
2**(7 / 6) * 3**(2 / 3) * pi**(-5 / 3) * theta(T)**(1 / 2) *
volume(rho, 1)**(2 / 3) *
PHI[0]**2)**(-4 / 3) * (HI[N] * integral_1_2(PHI[N]) + igrek(PHI[N]))
f2 = -teta(V) / T
f3 = -(T_a * sigma**(2 / 3) / T)**(1 / 2)
f4 = exp((gamma_0 - 2 / 3) * (B**2 + D**2) / B *
arctan(B * log(sigma) / (B**2 + D * (log(sigma) + D))))
f5 = B * log(sigma) / (B**2 + D * (log(sigma) + D))
f6 = B * log(sigma)
f7 = B**2 + D * (log(sigma) + D)
f8 = (T_a * sigma**(2 / 3) / T)**(1 / 2)
df8_V = T_a / T * 2 / 3 * sigma**(-1 / 3) * (-V_0 / V**2)
df3_V = -1 / 2 * (f8)**(-1 / 2) * df8_V
df5_V = (B * V / V_0 * (-V_0 / V**2) * f7 - f6 * D * V / V_0 *
(-V_0 / V**2)) / (f7**2)
darctanf5_V = 1 / (1 + f5**2) * df5_V
df4_V = f4 * (gamma_0 - 2 / 3) * (B**2 + D**2) / B * darctanf5_V
df2_V = -teta_0 / T * (2 / 3 * (V_0 / V)**(-1 / 3) *
(-V_0 / V**2) * f4 + sigma**(2 / 3) * df4_V)
df1_V = f1 * (df2_V + df3_V)
return -3 * R * T / (1 + f1) * df1_V
def teta(V):
sigma = V_0 / V
return teta_0 * sigma**(2 / 3) * exp(
(gamma_0 - 2 / 3) * (B**2 + D**2) / B * arctan(B * log(sigma) /
(B**2 + D *
(log(sigma) + D))))
print(ionic_contribution_pressure(1, volume(1, 1)))
def P(T, rho, Atom_weight, z):
PHI = progonka(T, rho, Atom_weight, z)
#return (P_e(T, rho) + theta(T) / volume(rho))
return 32 / (3 * pi**3) * (2 / pi)**(2 / 3) * (
2**(7 / 6) * 3**(2 / 3) * pi**(-5 / 3) * theta(T)**(1 / 2) *
volume(rho, 1)**(2 / 3) * (PHI[0])**2)**(-5 / 3) * integral_3_2(PHI[N])
def state_function_pressure(T, rho, Atom_weight):
#bilinear_interpolation_pressure(T, rho)
return 2.942 * 10**4 * bilinear_interpolation_pressure(
T, rho) + ionic_contribution_pressure(T, volume(
rho, Atom_weight)) + cold_curve_pressure(volume(rho, Atom_weight))
def S(T, rho):
return 0.9648 * 10**2 / Atom_weight * (S_e(T, rho) + 3 / 2 * log(
1836 * Atom_weight * theta(T) * volume(rho, 1)**
(2 / 3) / 2 / pi, e) + 5 / 3)
sigma_c = V_0_c / V
summ_1 = 0
summ_2 = 0
for i in range(1, 4):
summ_1 += a[i] / i * (sigma_c**(-i / 3) - 1)
for i in range(1, 3):
summ_2 += b[i] / i * (sigma_c**(i / 3) - 1)
print(summ_1)
print(summ_2)
return a[0] * V_0_c * log(
sigma_c) - 3 * V_0_c * summ_1 + 3 * V_0_c * summ_2
def cold_curve_pressure(V):
sigma_c = V_0_c / V
summ_1 = 0
summ_2 = 0
for i in range(1, 4):
summ_1 += a[i] / i * (i / 3 * (V_0_c / V)**(-i / 3 - 1) * V_0_c / V**2)
for i in range(1, 3):
summ_2 += b[i] / i * (i / 3 * (V_0_c / V)**(i / 3 - 1) *
(-V_0_c / V**2))
return -(a[0] * V_0_c * V / V_0_c *
(-V_0_c / V**2) - 3 * V_0_c * summ_1 + 3 * V_0_c * summ_2)
print("energy = ", cold_curve_energy(volume(1, 1)))