debyetools 1.4.2

Debye approximation implementation for the calculation of thermodynamic properties from ground-state atomistic simulations.

debyetools

Implementation of a tool for calculating self-consistent thermodynamic properties that can take into account all kinds of contributions to the free energy inluding explicit anharmonicity. The software presented here is based in the Debye approximation within the QHA using the crystal internal energetics parametrized at ground-state to project the thermodynamics properties at high temperatures.

Made by Javier Jofre: [e-mail] Please cite.

Requirements for Python module:

• numpy
• mpmath
• scipy

Requirements for Interface:

For the interface it will also be necesary:

• matplotlib
• PySide6

Installation

pip install --upgrade debyetools

Get started

To start getting familiar with the interface tProps you can download examples input files. The GUI can be launched by executing the interface script from the debyetools repository main folder:

python interface.py

Or you can launch inside python:

from debyetools.tpropsgui.gui import interface
interface()

Debye tools can also be used as a library. Example: heat capacity of Al fcc using 3rd order Birch-Murnaghan EOS

import numpy as np
import debyetools.potentials as potentials
from debyetools.ndeb import nDeb

# EOS parametrization
# =========================
EOS_parameters = [-3.607736520e+05, 9.929277050e-06, 7.729289055e+10, 4.604381753e+00]
EOS = potentials.BM()
EOS.fitEOS([0], [0], initial_parameters=EOS_parameters, fit=False)

# Other Contributions parametrization
# =========================
p_electronic = [3.8027342892e-01, -1.8875015171e-02, 5.3071034596e-04, -7.0100707467e-06]
mass = 0.026981500000000002
Tmelting = 933
p_defects = 8.46, 1.69, Tmelting, 0.1
p_anharmonicity = 0, 1
p_XS = 0, 0, 0
poissonsratio = 0.37

# F minimization using Slater approximaiton
# =========================
ndeb = nDeb(poissonsratio, mass, p_anharmonicity, EOS, p_electronic, p_defects, p_XS, mode='jjsl')
T_initial, T_final= 0.1, 1000
T = np.arange(T_initial, T_final, 10)
Pressure = 0
T, V = ndeb.min_G(T, EOS_parameters[0] * .9, P=Pressure)

# Evaluation of thermodynamic properties
# =========================
tprops_dict = ndeb.eval_props(T, V, P=Pressure)

To Do's:

• Improve error handling
• Add 'Compatible input files formats'
• Improve Documentation
• Add handling of anisotropic materials
• Prediction of explicit anharmonicity parameters