deploy: 124098c

feos-org · Oct 24, 2024 · 59dbd96 · 59dbd96
1 parent d8aba92
commit 59dbd96
Show file tree

Hide file tree

Showing 810 changed files with 2,369 additions and 49 deletions.
diff --git a/_sources/theory/eos/critical_points.md.txt b/_sources/theory/eos/critical_points.md.txt
@@ -0,0 +1,30 @@
+# Stability and critical points
+
+The implementation of critical points in $\text{FeO}_\text{s}$ follows the algorithm by [Michelsen and Mollerup](https://tie-tech.com/new-book-release/). A necessary condition for stability is the positive-definiteness of the quadratic form ([Heidemann and Khalil 1980](https://doi.org/10.1002/aic.690260510))
+
+$$\sum_{ij}\left(\frac{\partial^2 A}{\partial N_i\partial N_j}\right)_{T,V}\Delta N_i\Delta N_j$$
+
+The **spinodal** or limit of stability consists of the points for which the quadratic form is positive semi-definite. Following Michelsen and Mollerup, the matrix $M$ can be defined as
+
+$$M_{ij}=\sqrt{z_iz_j}\left(\frac{\partial^2\beta A}{\partial N_i\partial N_j}\right)$$
+
+with the molar compositon $z_i$. Further, the variable $s$ is introduced that acts on the mole numbers $N_i$ via
+
+$$N_i=z_i+su_i\sqrt{z_i}$$
+
+with $u_i$ the elements of the eigenvector of $M$ corresponding to the smallest eigenvector $\lambda_1$. Then, the limit of stability can be expressed as
+
+$$c_1=\left.\frac{\partial^2\beta A}{\partial s^2}\right|_{s=0}=\sum_{ij}u_iu_jM_{ij}=\lambda_1=0$$
+
+A **critical point** is defined as a stable point on the limit of stability. This leads to the second criticality condition
+
+$$c_2=\left.\frac{\partial^3\beta A}{\partial s^3}\right|_{s=0}=0$$
+
+The derivatives of the Helmholtz energy can be calculated efficiently in a single evaluation using [generalized hyper-dual numbers](https://doi.org/10.3389/fceng.2021.758090). The following methods of `State` are available to determine spinodal or critical points for different specifications:
+
+||specified|unkonwns|equations|
+|-|-|-|-|
+|`spinodal`|$T,N_i$|$\rho$|$c_1(T,\rho,N_i)=0$|
+|`critical_point`|$N_i$|$T,\rho$|$c_1(T,\rho,N_i)=0$<br/>$c_2(T,\rho,N_i)=0$|
+|`critical_point_binary_t`|$T$|$\rho_1,\rho_2$|$c_1(T,\rho_1,\rho_2)=0$<br/>$c_2(T,\rho_1,\rho_2)=0$|
+|`critical_point_binary_p`|$p$|$T,\rho_1,\rho_2$|$c_1(T,\rho_1,\rho_2)=0$<br/>$c_2(T,\rho_1,\rho_2)=0$<br/>$p(T,\rho_1,\rho_2)=p$|
diff --git a/_sources/theory/eos/index.md.txt b/_sources/theory/eos/index.md.txt
@@ -6,6 +6,7 @@ This section explains the thermodynamic principles and algorithms used for equat
    :maxdepth: 1
 
    properties
+   critical_points
 ```
 
 It is currently still under construction. You can help by [contributing](https://github.com/feos-org/feos/issues/70).
diff --git a/api/dft.html b/api/dft.html
@@ -1301,6 +1301,7 @@
 <ul>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../theory/eos/index.html">Equations of state</a><input class="toctree-checkbox" id="toctree-checkbox-142" name="toctree-checkbox-142" role="switch" type="checkbox"/><label for="toctree-checkbox-142"><div class="visually-hidden">Toggle navigation of Equations of state</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>
 <li class="toctree-l2"><a class="reference internal" href="../theory/eos/properties.html">Properties</a></li>
+<li class="toctree-l2"><a class="reference internal" href="../theory/eos/critical_points.html">Stability and critical points</a></li>
 </ul>
 </li>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../theory/dft/index.html">Classical density functional theory</a><input class="toctree-checkbox" id="toctree-checkbox-143" name="toctree-checkbox-143" role="switch" type="checkbox"/><label for="toctree-checkbox-143"><div class="visually-hidden">Toggle navigation of Classical density functional theory</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>

diff --git a/api/dippr.html b/api/dippr.html
@@ -1301,6 +1301,7 @@
 <ul>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../theory/eos/index.html">Equations of state</a><input class="toctree-checkbox" id="toctree-checkbox-142" name="toctree-checkbox-142" role="switch" type="checkbox"/><label for="toctree-checkbox-142"><div class="visually-hidden">Toggle navigation of Equations of state</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>
 <li class="toctree-l2"><a class="reference internal" href="../theory/eos/properties.html">Properties</a></li>
+<li class="toctree-l2"><a class="reference internal" href="../theory/eos/critical_points.html">Stability and critical points</a></li>
 </ul>
 </li>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../theory/dft/index.html">Classical density functional theory</a><input class="toctree-checkbox" id="toctree-checkbox-143" name="toctree-checkbox-143" role="switch" type="checkbox"/><label for="toctree-checkbox-143"><div class="visually-hidden">Toggle navigation of Classical density functional theory</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>

diff --git a/api/eos.html b/api/eos.html
@@ -1301,6 +1301,7 @@
 <ul>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../theory/eos/index.html">Equations of state</a><input class="toctree-checkbox" id="toctree-checkbox-142" name="toctree-checkbox-142" role="switch" type="checkbox"/><label for="toctree-checkbox-142"><div class="visually-hidden">Toggle navigation of Equations of state</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>
 <li class="toctree-l2"><a class="reference internal" href="../theory/eos/properties.html">Properties</a></li>
+<li class="toctree-l2"><a class="reference internal" href="../theory/eos/critical_points.html">Stability and critical points</a></li>
 </ul>
 </li>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../theory/dft/index.html">Classical density functional theory</a><input class="toctree-checkbox" id="toctree-checkbox-143" name="toctree-checkbox-143" role="switch" type="checkbox"/><label for="toctree-checkbox-143"><div class="visually-hidden">Toggle navigation of Classical density functional theory</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>

diff --git a/api/epcsaft.html b/api/epcsaft.html
@@ -1301,6 +1301,7 @@
 <ul>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../theory/eos/index.html">Equations of state</a><input class="toctree-checkbox" id="toctree-checkbox-142" name="toctree-checkbox-142" role="switch" type="checkbox"/><label for="toctree-checkbox-142"><div class="visually-hidden">Toggle navigation of Equations of state</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>
 <li class="toctree-l2"><a class="reference internal" href="../theory/eos/properties.html">Properties</a></li>
+<li class="toctree-l2"><a class="reference internal" href="../theory/eos/critical_points.html">Stability and critical points</a></li>
 </ul>
 </li>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../theory/dft/index.html">Classical density functional theory</a><input class="toctree-checkbox" id="toctree-checkbox-143" name="toctree-checkbox-143" role="switch" type="checkbox"/><label for="toctree-checkbox-143"><div class="visually-hidden">Toggle navigation of Classical density functional theory</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>

diff --git a/api/gc_pcsaft.html b/api/gc_pcsaft.html
@@ -1301,6 +1301,7 @@
 <ul>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../theory/eos/index.html">Equations of state</a><input class="toctree-checkbox" id="toctree-checkbox-142" name="toctree-checkbox-142" role="switch" type="checkbox"/><label for="toctree-checkbox-142"><div class="visually-hidden">Toggle navigation of Equations of state</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>
 <li class="toctree-l2"><a class="reference internal" href="../theory/eos/properties.html">Properties</a></li>
+<li class="toctree-l2"><a class="reference internal" href="../theory/eos/critical_points.html">Stability and critical points</a></li>
 </ul>
 </li>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../theory/dft/index.html">Classical density functional theory</a><input class="toctree-checkbox" id="toctree-checkbox-143" name="toctree-checkbox-143" role="switch" type="checkbox"/><label for="toctree-checkbox-143"><div class="visually-hidden">Toggle navigation of Classical density functional theory</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>

diff --git a/api/generated/feos.cubic.BinaryRecord.__init__.html b/api/generated/feos.cubic.BinaryRecord.__init__.html
@@ -1301,6 +1301,7 @@
 <ul>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../../theory/eos/index.html">Equations of state</a><input class="toctree-checkbox" id="toctree-checkbox-142" name="toctree-checkbox-142" role="switch" type="checkbox"/><label for="toctree-checkbox-142"><div class="visually-hidden">Toggle navigation of Equations of state</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>
 <li class="toctree-l2"><a class="reference internal" href="../../theory/eos/properties.html">Properties</a></li>
+<li class="toctree-l2"><a class="reference internal" href="../../theory/eos/critical_points.html">Stability and critical points</a></li>
 </ul>
 </li>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../../theory/dft/index.html">Classical density functional theory</a><input class="toctree-checkbox" id="toctree-checkbox-143" name="toctree-checkbox-143" role="switch" type="checkbox"/><label for="toctree-checkbox-143"><div class="visually-hidden">Toggle navigation of Classical density functional theory</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>

diff --git a/api/generated/feos.cubic.BinaryRecord.from_json.html b/api/generated/feos.cubic.BinaryRecord.from_json.html
@@ -1301,6 +1301,7 @@
 <ul>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../../theory/eos/index.html">Equations of state</a><input class="toctree-checkbox" id="toctree-checkbox-142" name="toctree-checkbox-142" role="switch" type="checkbox"/><label for="toctree-checkbox-142"><div class="visually-hidden">Toggle navigation of Equations of state</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>
 <li class="toctree-l2"><a class="reference internal" href="../../theory/eos/properties.html">Properties</a></li>
+<li class="toctree-l2"><a class="reference internal" href="../../theory/eos/critical_points.html">Stability and critical points</a></li>
 </ul>
 </li>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../../theory/dft/index.html">Classical density functional theory</a><input class="toctree-checkbox" id="toctree-checkbox-143" name="toctree-checkbox-143" role="switch" type="checkbox"/><label for="toctree-checkbox-143"><div class="visually-hidden">Toggle navigation of Classical density functional theory</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>

diff --git a/api/generated/feos.cubic.BinaryRecord.from_json_str.html b/api/generated/feos.cubic.BinaryRecord.from_json_str.html
@@ -1301,6 +1301,7 @@
 <ul>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../../theory/eos/index.html">Equations of state</a><input class="toctree-checkbox" id="toctree-checkbox-142" name="toctree-checkbox-142" role="switch" type="checkbox"/><label for="toctree-checkbox-142"><div class="visually-hidden">Toggle navigation of Equations of state</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>
 <li class="toctree-l2"><a class="reference internal" href="../../theory/eos/properties.html">Properties</a></li>
+<li class="toctree-l2"><a class="reference internal" href="../../theory/eos/critical_points.html">Stability and critical points</a></li>
 </ul>
 </li>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../../theory/dft/index.html">Classical density functional theory</a><input class="toctree-checkbox" id="toctree-checkbox-143" name="toctree-checkbox-143" role="switch" type="checkbox"/><label for="toctree-checkbox-143"><div class="visually-hidden">Toggle navigation of Classical density functional theory</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>

diff --git a/api/generated/feos.cubic.BinaryRecord.html b/api/generated/feos.cubic.BinaryRecord.html
@@ -1301,6 +1301,7 @@
 <ul>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../../theory/eos/index.html">Equations of state</a><input class="toctree-checkbox" id="toctree-checkbox-142" name="toctree-checkbox-142" role="switch" type="checkbox"/><label for="toctree-checkbox-142"><div class="visually-hidden">Toggle navigation of Equations of state</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>
 <li class="toctree-l2"><a class="reference internal" href="../../theory/eos/properties.html">Properties</a></li>
+<li class="toctree-l2"><a class="reference internal" href="../../theory/eos/critical_points.html">Stability and critical points</a></li>
 </ul>
 </li>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../../theory/dft/index.html">Classical density functional theory</a><input class="toctree-checkbox" id="toctree-checkbox-143" name="toctree-checkbox-143" role="switch" type="checkbox"/><label for="toctree-checkbox-143"><div class="visually-hidden">Toggle navigation of Classical density functional theory</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>

diff --git a/api/generated/feos.cubic.BinaryRecord.to_json_str.html b/api/generated/feos.cubic.BinaryRecord.to_json_str.html
@@ -1301,6 +1301,7 @@
 <ul>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../../theory/eos/index.html">Equations of state</a><input class="toctree-checkbox" id="toctree-checkbox-142" name="toctree-checkbox-142" role="switch" type="checkbox"/><label for="toctree-checkbox-142"><div class="visually-hidden">Toggle navigation of Equations of state</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>
 <li class="toctree-l2"><a class="reference internal" href="../../theory/eos/properties.html">Properties</a></li>
+<li class="toctree-l2"><a class="reference internal" href="../../theory/eos/critical_points.html">Stability and critical points</a></li>
 </ul>
 </li>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../../theory/dft/index.html">Classical density functional theory</a><input class="toctree-checkbox" id="toctree-checkbox-143" name="toctree-checkbox-143" role="switch" type="checkbox"/><label for="toctree-checkbox-143"><div class="visually-hidden">Toggle navigation of Classical density functional theory</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>

diff --git a/api/generated/feos.cubic.ChemicalRecord.__init__.html b/api/generated/feos.cubic.ChemicalRecord.__init__.html
@@ -1301,6 +1301,7 @@
 <ul>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../../theory/eos/index.html">Equations of state</a><input class="toctree-checkbox" id="toctree-checkbox-142" name="toctree-checkbox-142" role="switch" type="checkbox"/><label for="toctree-checkbox-142"><div class="visually-hidden">Toggle navigation of Equations of state</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>
 <li class="toctree-l2"><a class="reference internal" href="../../theory/eos/properties.html">Properties</a></li>
+<li class="toctree-l2"><a class="reference internal" href="../../theory/eos/critical_points.html">Stability and critical points</a></li>
 </ul>
 </li>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../../theory/dft/index.html">Classical density functional theory</a><input class="toctree-checkbox" id="toctree-checkbox-143" name="toctree-checkbox-143" role="switch" type="checkbox"/><label for="toctree-checkbox-143"><div class="visually-hidden">Toggle navigation of Classical density functional theory</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>

diff --git a/api/generated/feos.cubic.ChemicalRecord.from_json_str.html b/api/generated/feos.cubic.ChemicalRecord.from_json_str.html
@@ -1301,6 +1301,7 @@
 <ul>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../../theory/eos/index.html">Equations of state</a><input class="toctree-checkbox" id="toctree-checkbox-142" name="toctree-checkbox-142" role="switch" type="checkbox"/><label for="toctree-checkbox-142"><div class="visually-hidden">Toggle navigation of Equations of state</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>
 <li class="toctree-l2"><a class="reference internal" href="../../theory/eos/properties.html">Properties</a></li>
+<li class="toctree-l2"><a class="reference internal" href="../../theory/eos/critical_points.html">Stability and critical points</a></li>
 </ul>
 </li>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../../theory/dft/index.html">Classical density functional theory</a><input class="toctree-checkbox" id="toctree-checkbox-143" name="toctree-checkbox-143" role="switch" type="checkbox"/><label for="toctree-checkbox-143"><div class="visually-hidden">Toggle navigation of Classical density functional theory</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>

diff --git a/api/generated/feos.cubic.ChemicalRecord.from_smiles.html b/api/generated/feos.cubic.ChemicalRecord.from_smiles.html
@@ -1301,6 +1301,7 @@
 <ul>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../../theory/eos/index.html">Equations of state</a><input class="toctree-checkbox" id="toctree-checkbox-142" name="toctree-checkbox-142" role="switch" type="checkbox"/><label for="toctree-checkbox-142"><div class="visually-hidden">Toggle navigation of Equations of state</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>
 <li class="toctree-l2"><a class="reference internal" href="../../theory/eos/properties.html">Properties</a></li>
+<li class="toctree-l2"><a class="reference internal" href="../../theory/eos/critical_points.html">Stability and critical points</a></li>
 </ul>
 </li>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../../theory/dft/index.html">Classical density functional theory</a><input class="toctree-checkbox" id="toctree-checkbox-143" name="toctree-checkbox-143" role="switch" type="checkbox"/><label for="toctree-checkbox-143"><div class="visually-hidden">Toggle navigation of Classical density functional theory</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>

diff --git a/api/generated/feos.cubic.ChemicalRecord.html b/api/generated/feos.cubic.ChemicalRecord.html
@@ -1301,6 +1301,7 @@
 <ul>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../../theory/eos/index.html">Equations of state</a><input class="toctree-checkbox" id="toctree-checkbox-142" name="toctree-checkbox-142" role="switch" type="checkbox"/><label for="toctree-checkbox-142"><div class="visually-hidden">Toggle navigation of Equations of state</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>
 <li class="toctree-l2"><a class="reference internal" href="../../theory/eos/properties.html">Properties</a></li>
+<li class="toctree-l2"><a class="reference internal" href="../../theory/eos/critical_points.html">Stability and critical points</a></li>
 </ul>
 </li>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../../theory/dft/index.html">Classical density functional theory</a><input class="toctree-checkbox" id="toctree-checkbox-143" name="toctree-checkbox-143" role="switch" type="checkbox"/><label for="toctree-checkbox-143"><div class="visually-hidden">Toggle navigation of Classical density functional theory</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>

diff --git a/api/generated/feos.cubic.ChemicalRecord.to_json_str.html b/api/generated/feos.cubic.ChemicalRecord.to_json_str.html
@@ -1301,6 +1301,7 @@
 <ul>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../../theory/eos/index.html">Equations of state</a><input class="toctree-checkbox" id="toctree-checkbox-142" name="toctree-checkbox-142" role="switch" type="checkbox"/><label for="toctree-checkbox-142"><div class="visually-hidden">Toggle navigation of Equations of state</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>
 <li class="toctree-l2"><a class="reference internal" href="../../theory/eos/properties.html">Properties</a></li>
+<li class="toctree-l2"><a class="reference internal" href="../../theory/eos/critical_points.html">Stability and critical points</a></li>
 </ul>
 </li>
 <li class="toctree-l1 has-children"><a class="reference internal" href="../../theory/dft/index.html">Classical density functional theory</a><input class="toctree-checkbox" id="toctree-checkbox-143" name="toctree-checkbox-143" role="switch" type="checkbox"/><label for="toctree-checkbox-143"><div class="visually-hidden">Toggle navigation of Classical density functional theory</div><i class="icon"><svg><use href="#svg-arrow-right"></use></svg></i></label><ul>