Computational thermodynamics

Last updated

Computational thermodynamics is the use of computers to simulate thermodynamic problems specific to materials science, particularly used in the construction of phase diagrams. [1] [2]

Contents

Several open and commercial programs exist to perform these operations. The concept of the technique is minimization of Gibbs free energy of the system; the success of this method is due not only to properly measuring thermodynamic properties, such as those in the list of thermodynamic properties, but also due to the extrapolation of the properties of metastable allotropes of the chemical elements.

History

The computational modeling of metal-based phase diagrams, which dates back to the beginning of the previous century mainly by Johannes van Laar and to the modeling of regular solutions, has evolved in more recent years to the CALPHAD (CALculation of PHAse Diagrams). [3] This has been pioneered by American metallurgist Larry Kaufman since the 1970s. [4] [5] [6]

Current state

Computational thermodynamics may be considered a part of materials informatics and is a cornerstone of the concepts behind the materials genome project. While crystallographic databases are used mainly as a reference source, thermodynamic databases represent one of the earliest examples of informatics, as these databases were integrated into thermochemical computations to map phase stability in binary and ternary alloys. [7] Many concepts and software used in computational thermodynamics are credited to the SGTE Group, a consortium devoted to the development of thermodynamic databases; the open elements database is freely available [8] based on the paper by Dinsdale. [9] This so-called "unary" system proves to be a common basis for the development of binary and multiple systems and is used by both commercial and open software in this field.

However, as stated in recent[ when? ] CALPHAD papers and meetings, such a Dinsdale/SGTE database will likely need to be corrected over time despite the utility in keeping a common base. In this case, most published assessments will likely have to be revised, similarly to rebuilding a house due to a severely broken foundation. This concept has also been depicted as an "inverted pyramid." [10] Merely extending the current approach (limited to temperatures above room temperature) is a complex task. [11] PyCalphad, a Python library, was designed to facilitate simple computational thermodynamics calculation using open source code. [12] In complex systems, computational methods such as CALPHAD are employed to model thermodynamic properties for each phase and simulate multicomponent phase behavior. [13] The application of CALPHAD to high pressures in some important applications, which are not restricted to one side of materials science like the Fe-C system, [14] confirms experimental results by using computational thermodynamic calculations of phase relations in the Fe–C system at high pressures. Other scientists even considered viscosity and other physical parameters, which are beyond the domain of thermodynamics. [15]

Future developments

There is still a gap between ab initio methods [16] and operative computational thermodynamics databases. In the past, a simplified approach introduced by the early works of Larry Kaufman, based on Miedema's Model, was employed to check the correctness of even the simplest binary systems. However, relating the two communities to Solid State Physics and Materials Science remains a challenge, [17] as it has been for many years. [18] Promising results from ab initio quantum mechanics molecular simulation packages like VASP are readily integrated in thermodynamic databases with approaches like Zentool. [19] A relatively easy way to collect data for intermetallic compounds is now possible by using Open Quantum Materials Database. A series of papers focused on the concept of Zentropy has been proposed by prof. Z.K. Liu and his research group has been recently proposed [20]

See also

Related Research Articles

<span class="mw-page-title-main">Computational chemistry</span> Branch of chemistry

Computational chemistry is a branch of chemistry that uses computer simulations to assist in solving chemical problems. It uses methods of theoretical chemistry incorporated into computer programs to calculate the structures and properties of molecules, groups of molecules, and solids. The importance of this subject stems from the fact that, with the exception of some relatively recent findings related to the hydrogen molecular ion, achieving an accurate quantum mechanical depiction of chemical systems analytically, or in a closed form, is not feasible. The complexity inherent in the many-body problem exacerbates the challenge of providing detailed descriptions of quantum mechanical systems. While computational results normally complement information obtained by chemical experiments, it can occasionally predict unobserved chemical phenomena.

<span class="mw-page-title-main">Phase diagram</span> Chart used to show conditions at which physical phases of a substance occur

A phase diagram in physical chemistry, engineering, mineralogy, and materials science is a type of chart used to show conditions at which thermodynamically distinct phases occur and coexist at equilibrium.

<span class="mw-page-title-main">Molecular dynamics</span> Computer simulations to discover and understand chemical properties

Molecular dynamics (MD) is a computer simulation method for analyzing the physical movements of atoms and molecules. The atoms and molecules are allowed to interact for a fixed period of time, giving a view of the dynamic "evolution" of the system. In the most common version, the trajectories of atoms and molecules are determined by numerically solving Newton's equations of motion for a system of interacting particles, where forces between the particles and their potential energies are often calculated using interatomic potentials or molecular mechanical force fields. The method is applied mostly in chemical physics, materials science, and biophysics.

<span class="mw-page-title-main">Lennard-Jones potential</span> Model of intermolecular interactions

In computational chemistry, molecular physics, and physical chemistry, the Lennard-Jones potential is an intermolecular pair potential. Out of all the intermolecular potentials, the Lennard-Jones potential is probably the one that has been the most extensively studied. It is considered an archetype model for simple yet realistic intermolecular interactions. The Lennard-Jones potential is often used as a building block in molecular models for more complex substances. Many studies of the idealized "Lennard-Jones substance" use the potential to understand the physical nature of matter.

<span class="mw-page-title-main">Molecular modelling</span> Discovering chemical properties by physical simulations

Molecular modelling encompasses all methods, theoretical and computational, used to model or mimic the behaviour of molecules. The methods are used in the fields of computational chemistry, drug design, computational biology and materials science to study molecular systems ranging from small chemical systems to large biological molecules and material assemblies. The simplest calculations can be performed by hand, but inevitably computers are required to perform molecular modelling of any reasonably sized system. The common feature of molecular modelling methods is the atomistic level description of the molecular systems. This may include treating atoms as the smallest individual unit, or explicitly modelling protons and neutrons with its quarks, anti-quarks and gluons and electrons with its photons.

<span class="mw-page-title-main">Force field (chemistry)</span> Concept on molecular modeling

In the context of chemistry, molecular physics, physical chemistry, and molecular modelling, a force field is a computational model that is used to describe the forces between atoms within molecules or between molecules as well as in crystals. Force fields are a variety of interatomic potentials. More precisely, the force field refers to the functional form and parameter sets used to calculate the potential energy of a system on the atomistic level. Force fields are usually used in molecular dynamics or Monte Carlo simulations. The parameters for a chosen energy function may be derived from classical laboratory experiment data, calculations in quantum mechanics, or both. Force fields utilize the same concept as force fields in classical physics, with the main difference being that the force field parameters in chemistry describe the energy landscape on the atomistic level. From a force field, the acting forces on every particle are derived as a gradient of the potential energy with respect to the particle coordinates.

<span class="mw-page-title-main">Scheil equation</span> Metallurgical equation for the redistribution of solutes during solidification of an alloy

In metallurgy, the Scheil-Gulliver equation describes solute redistribution during solidification of an alloy.

Car–Parrinello molecular dynamics or CPMD refers to either a method used in molecular dynamics or the computational chemistry software package used to implement this method.

<span class="mw-page-title-main">Geothermobarometry</span> History of rock pressure and temperature

Geothermobarometry is the methodology for estimating the pressure and temperature history of rocks. Geothermobarometry is a combination of geobarometry, where the pressure attained by a mineral assemblage is estimated, and geothermometry where the temperature attained by a mineral assemblage is estimated.

CALPHAD stands for CALculation of PHAse Diagrams, a methodology introduced in 1970 by Larry Kaufman. An equilibrium phase diagram is usually a diagram with axes for temperature and composition of a chemical system. It shows the regions where substances or solutions are stable and regions where two or more of them coexist. Phase diagrams are a very powerful tool for predicting the state of a system under different conditions and were initially a graphical method to rationalize experimental information on states of equilibrium. In complex systems, computational methods such as CALPHAD are employed to model thermodynamic properties for each phase and simulate multicomponent phase behavior. The CALPHAD approach is based on the fact that a phase diagram is a manifestation of the equilibrium thermodynamic properties of the system, which are the sum of the properties of the individual phases. It is thus possible to calculate a phase diagram by first assessing the thermodynamic properties of all the phases in a system.

<span class="mw-page-title-main">Spinodal</span>

In thermodynamics, the limit of local stability against phase separation with respect to small fluctuations is clearly defined by the condition that the second derivative of Gibbs free energy is zero.

Binary compounds of hydrogen are binary chemical compounds containing just hydrogen and one other chemical element. By convention all binary hydrogen compounds are called hydrides even when the hydrogen atom in it is not an anion. These hydrogen compounds can be grouped into several types.

Quantemol Ltd is based in University College London initiated by Professor Jonathan Tennyson FRS and Dr. Daniel Brown in 2004. The company initially developed a unique software tool, Quantemol-N, which provides full accessibility to the highly sophisticated UK molecular R-matrix codes, used to model electron polyatomic molecule interactions. Since then Quantemol has widened to further types of simulation, with plasmas and industrial plasma tools, in Quantemol-VT in 2013 and launched in 2016 a sustainable database Quantemol-DB, representing the chemical and radiative transport properties of a wide range of plasmas.

Geochemical modeling or theoretical geochemistry is the practice of using chemical thermodynamics, chemical kinetics, or both, to analyze the chemical reactions that affect geologic systems, commonly with the aid of a computer. It is used in high-temperature geochemistry to simulate reactions occurring deep in the Earth's interior, in magma, for instance, or to model low-temperature reactions in aqueous solutions near the Earth's surface, the subject of this article.

Miedema's model is a semi-empirical approach for estimating the heat of formation of solid or liquid metal alloys and compounds in the framework of thermodynamic calculations for metals and minerals. It was developed by the Dutch scientist Andries Rinse Miedema while working at Philips Natuurkundig Laboratorium. It may provide or confirm basic enthalpy data needed for the calculation of phase diagrams of metals, via CALPHAD or ab initio quantum chemistry methods. For a binary system composed by elements A and B, a generic Miedema Formula could be cast as where terms Phi and nwS are explained and reported below. For a binary system the physical picture could be simplified by considering a relatively simpler function of the difference of these three physical parameters resulting in a more complex form

<span class="mw-page-title-main">Mark S. Ghiorso</span> American geochemist

Mark S. Ghiorso is an American geochemist who resides in Seattle, Washington. He is best known for creating MELTS, a software tool for thermodynamic modeling of phase equilibria in magmatic systems.

Computational materials science and engineering uses modeling, simulation, theory, and informatics to understand materials. The main goals include discovering new materials, determining material behavior and mechanisms, explaining experiments, and exploring materials theories. It is analogous to computational chemistry and computational biology as an increasingly important subfield of materials science.

While chemically pure materials have a single melting point, chemical mixtures often partially melt at the solidus temperature (TS or Tsol), and fully melt at the higher liquidus temperature (TL or Tliq). The solidus is always less than or equal to the liquidus, but they need not coincide. If a gap exists between the solidus and liquidus it is called the freezing range, and within that gap, the substance consists of a mixture of solid and liquid phases (like a slurry). Such is the case, for example, with the olivine (forsterite-fayalite) system, which is common in Earth's mantle.

References

  1. Liu, Zi-Kui; Wang, Yi (2016). Computational Thermodynamics of Materials. Cambridge University Press. ISBN   9780521198967.
  2. Liu, Zi-Kui; Wang, Liu (2020). "ComputationalComputational thermodynamics and its applications". Acta Materialia. 200: 745. Bibcode:2020AcMat.200..745L. doi: 10.1016/j.actamat.2020.08.008 .
  3. Fabrichnaya, Olga; Saxena, Surendra K.; Richet, Pascal; Westrum, Edgar F. (14 March 2013). Thermodynamic Data, Models, and Phase Diagrams in Multicomponent Oxide Systems: An Assessment for Materials and Planetary Scientists Based on Calorimetric, Volumetric and Phase Equilibrium Data. Springer Science & Business Media. ISBN   9783662105047.
  4. L Kaufman and H Bernstein, Computer Calculation of Phase Diagrams, Academic Press N Y (1970) ISBN   0-12-402050-X [ page needed ]
  5. N Saunders and P Miodownik, Calphad, Pergamon Materials Series, Vol 1 Ed. R W Cahn (1998) ISBN   0-08-042129-6 [ page needed ]
  6. H L Lukas, S G Fries and B Sundman, Computational Thermodynamics, the Calphad Method, Cambridge University Press (2007) ISBN   0-521-86811-4 [ page needed ]
  7. K., Saxena, Surendra (1993). Thermodynamic Data on Oxides and Silicates : an Assessed Data Set Based on Thermochemistry and High Pressure Phase Equilibrium. Chatterjee, Nilanjan., Fei, Yingwei., Shen, Guoyin. Berlin, Heidelberg: Springer Berlin Heidelberg. ISBN   9783642783326. OCLC   840299125.{{cite book}}: CS1 maint: multiple names: authors list (link)
  8. http://www.crct.polymtl.ca/sgte/unary50.tdb%5B%5D%5B%5D
  9. Dinsdale, A.T. (1991). "SGTE data for pure elements". Calphad. 15 (4): 317–425. doi:10.1016/0364-5916(91)90030-N.
  10. "MICRESS® - the MICRostructure Evolution Simulation Software" (PDF).
  11. "Computational Materials Engineering" (PDF).
  12. Otis, Richard; Liu, Zi-Kui (2017). "Pycalphad: CALPHAD-based Computational Thermodynamics in Python". Journal of Open Research Software. 5: 1. doi: 10.5334/jors.140 .
  13. L., Lukas, H. (2007). Computational thermodynamics : the CALPHAD method. Fries, Suzana G., Sundman, Bo. Cambridge: Cambridge University Press. ISBN   978-0521868112. OCLC   663969016.{{cite book}}: CS1 maint: multiple names: authors list (link)
  14. Fei, Yingwei; Brosh, Eli (2014). "Experimental study and thermodynamic calculations of phase relations in the Fe–C system at high pressure". Earth and Planetary Science Letters. 408: 155–62. Bibcode:2014E&PSL.408..155F. doi:10.1016/j.epsl.2014.09.044.
  15. Zhang, Fan; Du, Yong; Liu, Shuhong; Jie, Wanqi (2015). "Modeling of the viscosity in the AL–Cu–Mg–Si system: Database construction". Calphad. 49: 79–86. doi:10.1016/j.calphad.2015.04.001.
  16. P. Turchi AB INITIO AND CALPHAD THERMODYNAMICS OF MATERIALS https://e-reports-ext.llnl.gov/pdf/306920.pdf
  17. J. A. Alonso and N. H. March Electrons in Metals and Alloys http://www.sciencedirect.com/science/book/9780120536207%5B%5D
  18. Proceedings of the International Symposium on Thermodynamics of Alloys - 1st Edition. 1 January 1981. ISBN   978-1-4832-2782-5 . Retrieved 1 July 2023.{{cite book}}: |website= ignored (help)[ full citation needed ][ page needed ]
  19. Manualzengen.cnrs.fr Archived 18 April 2016 at the Wayback Machine
  20. Liu, Zi-Kui (2023). "Thermodynamics and its prediction and CALPHAD modeling: Review, state of the art, and perspectives". Calphad. 82: 102580. arXiv: 2301.02132 . doi:10.1016/j.calphad.2023.102580. S2CID   259138637 via ScienceDirect.

University Courses on Computational Thermodynamics