RESEARCH STARTER
Phase equilibria
Phase equilibria is a scientific concept that explores the conditions under which different phases of matter (solids, liquids, gases) coexist in a stable state. It is particularly relevant in geology, where it helps to understand the formation processes of rocks, especially igneous and metamorphic types. By analyzing the mineral assemblages found in these rocks, researchers can infer the temperature, pressure, and chemical composition conditions that existed during their formation.
The principles of phase equilibria are rooted in the concept of chemical equilibrium, which posits that a chemical reaction reaches a stable state when the chemical potentials of reactants and products are equal. This understanding allows scientists to predict how changes in temperature and pressure can affect rock formation. The Gibbs phase rule is a key tool in this area, which provides a mathematical relationship between the number of phases present, the number of chemical components, and the degrees of freedom in a system.
In practice, phase equilibria studies involve creating phase diagrams to graphically represent how different minerals behave under varying conditions. These diagrams are essential in both metamorphic and igneous petrology, helping geologists to infer the evolutionary history of rocks and the geological processes at play. The ultimate goal of phase equilibria studies is to deepen our understanding of how rocks form and the dynamic processes that shape the Earth's crust.
Authored By: Lowell, Gary R. 1 of 4
Published In: 2022 2 of 4
- Related Topics:
3 of 4
- Related Articles:A Barrovian Overprint Superimposed on Partly Preserved High-Pressure Metamorphic Rocks: Pressure–Temperature–Time Evolution in the Ordovician Corner Brook Collisional Complex (Western Newfoundland, Canada).;Advances in Igneous Petrology: Coupling Chemical Maps with Thermodynamic Models to Tackle Mush Crystallization Dynamics and Crystal–Melt Segregation.;Age of Skarn Formation and Meteoric Water Infiltration in the Willsboro–Lewis Wollastonite District (Adirondack Highlands, New York).;Identification of Neoarchean UHT (Ultrahigh‐Temperature) Metamorphism From Intermediate Migmatitic Granulites in the Yinshan Block, North China Craton.;Influence of the mineral gangue on pellets softening melting properties in blast furnace. Experimental study of phase equilibria during the melting of pre-reduced pellets.
4 of 4
Full Article
The mineral assemblages in most igneous and metamorphic rocks preserve a record of the chemical equilibrium related to the initial rock-forming process. Phase equilibria studies attempt to determine quantitatively the pressure and temperature conditions of rock formations from these mineral assemblages.
Chemical Equilibrium
The traditional methods of studying rock bodies are descriptive in nature and involve mapping large-scale outcrop features in the field and detailed microscopic observations of rock textures and mineralogy in the laboratory. These methods, which are often successful by themselves, are supplemented by a second, more theoretical approach wherein rocks are treated as chemical systems and the principles of phase equilibria are applied to determine the conditions of their origin. A full and complete description of a rock body is still required, but that is no longer the main goal of petrologic study. The principles of phase equilibria are simply the laws that govern the attainment of equilibrium of chemical reactions such as A + B = C + D, where A and B are known as reactants, and C and D are known as products.
Before exploring how these principles cast light on rock-forming processes, the concept of chemical equilibrium must first be developed. By analogy with gravitational potential, there must exist a similar tendency in chemical systems to lower their energy state through chemical reactions. Reasoning along these lines, J. Willard Gibbs introduced the term “chemical potential” to describe the flow of chemical components from one site (of high potential) to another (of lower potential) during reactions that lead a chemical system toward its lowest energy state. The total energy available to drive a chemical reaction must therefore be the sum of the chemical potentials of each component in the system multiplied by the number of moles of each component. The usual definition of chemical potential of a phase (or pure substance) is “the molar free energy” (or free energy per mole). This simple statement leads to a workable, three-part definition of chemical equilibrium, which is central to the understanding of phase equilibria.
First, if the chemical potential of a component is the same on either side of a reaction equation, the component can have no tendency to participate in the reaction. Second, in a multicomponent system consisting of several phases under uniform temperature (T) and pressure (P), equilibrium must prevail when the chemical potential of each component is the same in all phases in which the component is present. Third, the condition of equilibrium is one of maximum chemical stability. The second part of the definition is equivalent to saying that, for a given chemical reaction, the free energy of the reactants must equal the free energy of the products if a condition of chemical equilibrium prevails under fixed conditions of T and P. If either T or P changes, the system is no longer in equilibrium, but it will immediately adjust itself in such a way as to “moderate” the effect of the disturbing factors. The last statement is known as the moderation theorem, or Le Châtelier's principle. If the free energy on the product side of an equation is less than the free energy on the reactant side, the reaction will be spontaneous. If the opposite is true, no reaction is possible. Scientists are thus able to predict the result of any chemical reaction if the free energy of the reactants and products under the reaction conditions (T, P) is known.
Equilibrium State of a System
The most significant processes of rock formation—magmatism, metamorphism, and sedimentation—all involve large-scale flow of energy and movement of matter that produce an uneven distribution of chemical potential. Inevitably, the result must be chemical reactions tending to restore these natural systems to a state of equilibrium. The equilibrium state of a system is governed by its bulk composition (X), temperature (T), and pressure (P). For most geological processes, T and P change slowly relative to the rates of most chemical reactions, which means that most rock-forming reactions may be considered to take place under constant T and P, and, if X also remains constant, most such reactions should easily attain chemical equilibrium.
In geology, the major concern is not so much with the achievement of equilibrium but rather with the preservation of equilibrium mineral assemblages through hundreds of millions of years, which must follow before deep-seated rocks are finally exposed at the surface. Rocks formed at depth must clearly undergo significant reductions in T and P prior to exposure at the surface, and there are several mechanisms that may induce changes in X during this lengthy period. Geologists are acutely aware of the implications of the moderation theorem: Retrograde metamorphism, mineralogical inversions and exsolutions, hydrothermal alteration, and weathering are but a few of the processes that could trigger reequilibration in rock bodies before they are exposed for study. Fortunately, microscopic studies coupled with Gibbs's pioneering work in phase equilibria provide the means to discern whether a given mineralogical assemblage preserves a former equilibrium.
Gibbs Phase Rule
The “phase rule,” also called the Gibbs phase rule, was initially derived by Gibbs in 1878 from the mathematical formulas of thermodynamics. The phase rule, which is fully applicable to all chemical systems, expresses the relationship between the governing variables T, P, and X, and the number of phases that may coexist in a state of equilibrium. Usually, the phase rule is expressed in equation form as P + F = C + 2, where P = number of phases, C = number of chemical components, and F = degrees of freedom possessed by the system (normally T, P, X). Phases are chemically pure, physically separable subparts of the system and may be gases, liquids, or solids. In the formal sense, C is the minimum number of chemical entities needed to define completely the composition of each reactant and product phase in a given reaction.
Although the objective of phase rule applications is to determine F for major rock-forming reactions, a far simpler situation that could be experimentally verified in any high school laboratory or even in an ordinary kitchen may be considered. Pure water (H2O) boils at T = 100 degrees Celsius at sea level (P = 1 bar; atmospheric pressure). The effect of dissolving common salt (NaCl) in water is to raise the boiling temperature approximately 0.8 degrees Celsius for each mole percent of NaCl in the liquid phase. The steam given off by boiling is pure H2O, and, therefore, the salt concentration in the remaining liquid must progressively increase with temperature during boiling.
To apply the phase rule to this simple system, one first must tally up the participating phases: there is steam and there is liquid salt solution, and one must conclude that P = 2. Both pure water and salt are required to form these coexisting phases and, therefore, C = 2. The phase rule for this process (boiling), under conditions of fixed pressure (P = 1 bar), tells one that
F=C-P+2=2-2+2=2
Therefore, the system has two degrees of freedom and is said to be “divariant,” which means that, because one degree of freedom is utilized by fixing P = 1 bar, only one additional variable needs to be known to specify completely the state of the system. That may be either T (boiling temperature) or X (composition of the boiling solution); in other words, T and X are dependent variables at constant P. This T-X compositional dependence is easily determined for P = 1 bar by direct experiment. If the resulting T-X data were graphically plotted, the diagram would indicate, for example, that the boiling temperature for an 8 mole percent solution is close to 106 degrees Celsius. Conversely, if it were known only that the boiling temperature of a salt solution were 106 degrees Celsius, that would necessitate a solution concentration of 8 mole percent at P = 1 bar.
Application of the Phase Rule
To apply the phase rule to a reaction that has some geological significance, consider the appearance of diopside in siliceous dolomite during contact metamorphism by the reaction
CaMg(CO3)2(dolomite)+2SiO2(quartz)=CaMgSi2O6(diopside)+2CO2(gas)
The reaction involves three mineral phases and a fugitive gas phase, which is necessarily lost from the rock if diopside appears; P = 4. Note that the Ca:Mg (calcium-to-magnesium) ratio is the same in the reactant phase (dolomite) and the product phase (diopside); consequently, the minimum number of components needed to define the compositions of the four phases in the reaction is three. Therefore, C = 3. Substituting these values into the phrase rule, one obtains
F=C−P+2=3−4+2=1
and the reaction has one degree of freedom and is said to be univariant. For any given reaction pressure, the phases appearing in that reaction can coexist in equilibrium at one—and only one—temperature.
Univariant reactions are of great interest to petrologists because it is frequently possible to estimate the depth of a rock-forming reaction and, therefore, P, from field relationships. It is then a simple matter to estimate the reaction temperature from the P-T diagram for the univariant reaction involved. For example, suppose that field relationships lead to the conclusion that diopside-bearing contact metamorphic rocks formed by the reaction above at an estimated depth of 4 kilometers. The pressure equivalent of this depth is about 1,000 bars. The experimentally derived P-T diagram indicates a reaction temperature of about 450 degrees Celsius.
Mineralogical Phase Rule
The principles of phase equilibria were first applied to rocks in 1911 by Swiss geochemist V. M. Goldschmidt in his classic account of contact metamorphism in the Oslo area in Norway. Countless similar attempts have followed, with most authors concluding, as did Goldschmidt, that rocks, in general, record a state of chemical equilibrium governed by the temperature and pressure prevailing at their time of origin. The major generalization emerging from eighty years of such studies is that rocks with a large number of mineral phases tend to have a low number of degrees of freedom. Goldschmidt recognized that, at Oslo, metamorphism must be controlled mainly by T and P and that divariant equilibrium (F = 2) is the general case for isochemical rock-forming reactions. For F = 2, the phase rule reduces to P = C. Dubbed “the mineralogical phase rule,” this equation cannot be mathematically derived from thermodynamic principles. It simply reflects the common case in nature where rock-forming processes are approximately isochemical, and their phase equilibria are controlled by both T and P operating as independent variables. It follows that rocks recording univariant equilibria are more restricted in occurrence than those recording divariant equilibria. Similarly, the rarity of rocks recording invariant equilibria (F = 0) can be understood, as the phase rule requires that both a unique T and a unique P be maintained during their formation. Rocks of complex mineralogy are occasionally encountered in which the number of phases, P, exceeds C + 2 and, consequently, the value of F is negative. Phase rule departures, indicated by negative F values, are a sure sign of disequilibrium and serve as a reminder that not all mineral assemblages can be treated by the methods of phase equilibria.
Petrogenetic Grids
The temperature and pressure range in which any mineral may exist is limited, and it is the task of experimental petrology to ascertain these limits for rock-forming minerals. The data resulting from this experimental work are utilized to construct “phase diagrams.” Such diagrams show the effects of changing P, T, and X values on mineral stability fields and are, therefore, simply graphical expressions of the phase rule. In metamorphic petrology, the major use of phase equilibria data has been for construction of “petrogenetic grids.” A petrogenetic grid, initially conceived by Canadian petrologist N. L. Bowen in 1940, is a P-T diagram on which experimentally derived univariant reaction curves are plotted for a particular metamorphic rock type (for example, blueschists, marbles, calcsilicates, or pelites). The value of such “grids” lies in the fact that each natural equilibrium assemblage recognized in the field will fall within a definite P-T pigeonhole and thereby inform the field geologist immediately of the P-T conditions of the metamorphic terrane under study. This goal, so simple in concept, has proved elusive even after half a century of vigorous experimental, theoretical, and field effort. The problem lies in the lack of truly univariant reactions. For nearly seventy years, field geologists mapped isograds recording the “first appearance” of notable zone minerals such as biotite, garnet, staurolite, kyanite, and sillimanite under the impression that they represented the intersection of the ground surface with a plane of univariant equilibrium. Virtually all such “isograds” have proved to be the result of divariant equilibria and thus plot as a “band”—which may be rather wide—on a P-T diagram. This undesirable result has the effect of “smearing” grid boundaries and rendering them less useful as metamorphic indicators.
The general absence of univariant reactions in metamorphic rocks was eventually recognized because of theoretical and experimental advances in phase equilibria studies. The problem stems from the fact that most mineral phases participating in metamorphic reactions are solid solutions of variable composition, and the most common reactions lead to the release of a fluid, the composition of which may vary with time. Each of these effects introduces an additional degree of freedom in phase rule terms, and, as a result, virtually all important reactions are divariant.
In spite of these difficulties, petrogenetic grids, based on divariant and quasi-univariant equilibria, have gradually evolved for all major metamorphic rock types. These are not the simple, quantitative grids envisioned by Bowen, but they do provide quick, reliable, and fairly narrow estimates of P-T conditions for common metamorphic mineral assemblages. Modern grids, continually subject to refinement, are phase equilibria's greatest contribution to metamorphic petrology.
Since 2022, efficient open-source Gibbs free energy minimization tools such as Mineral Assemblage Gibbs Energy Minimization (MAGEMin) have revolutionized phase equilibria modeling by enabling rapid, high-resolution calculations of stable mineral assemblages and melt properties across complex multicomponent systems. These advances allow seamless integration with large-scale geodynamic simulations. This shift transforms phase equilibria from mainly interpretive tools for ancient rocks into powerful predictive engines for active Earth processes, significantly reducing uncertainties in P-T paths and strengthening quantitative links between petrology, tectonics, and mantle dynamics.
Application to Igneous Petrology
In the area of igneous petrology, phase equilibria methods and data have become indispensable. Natural rocks, spanning the compositional spectrum, are melted under strictly controlled laboratory conditions to determine solidus and liquidus temperatures at pressures ranging from 1 to 35,000 bars. The results of such experiments place tight constraints on the depths and temperatures of magma generation. They also permit the experimentalist to explore P-T effects on partial melting (anatexis) in terms of melt composition and refractory solid phases. The resulting phase diagrams, like metamorphic grids, permit petrologists to “see” deep into the crust and upper mantle and to test hypotheses dealing with the origin of magma.
For nearly a century, igneous petrologists have studied crystal-melt equilibria of simplified, synthetic melts as models for complex, natural magmas. The objective is to reduce the number of equilibrium phases by the elimination of minor components of real magmas. Studies of this type were introduced by Bowen at the Geophysical Laboratory of the Carnegie Institution in Washington, D.C. Through its many subsequent researchers, this laboratory published hundreds of phase diagrams and earned a reputation for meticulous and exhaustive experimental work.
Phase equilibrium studies have provided a rather complete understanding of two fundamentally different modes of magma crystallization. Equilibrium crystallization occurs when P-T-X conditions change so slowly that chemical reactions within the melt are able to maintain the state of chemical equilibrium. Conversely, fractional crystallization results when changes in P-T-X conditions outpace the compensating reactions. This disequilibrium process greatly influences the behavior of natural magmas and extends the range of melt compositions that can be derived from a given parent magma. This latter type of behavior, recognized through the early phase equilibria studies of the Geophysical Laboratory, is the major factor in explaining the compositional diversity of igneous rocks.
The relatively simple phase diagrams of synthetic systems unraveled the complexities of sequential crystallization, cast light on the mechanics of crystal nucleation, and exposed the crucial role that water plays in magmatic processes. Collectively, these diagrams are the foundation of modern igneous petrology.
Goal of Phase Equilibria Studies
The refined symbolic notation and elegant mathematical derivations of thermodynamics are likely to remain unappreciated by the majority of laypersons and geologists alike. It is precisely these formalisms, however, that place phase equilibria on a quantitative footing and permit calculation of mineral stability fields from compositional data. Future development in the area of phase equilibria will follow this theoretical line.
The qualitative form of phase equilibria is expressed in phase diagrams rather than equations. Such diagrams have been a major part of petrology since the 1950s. Historically, emphasis in phase equilibria studies has been on high-temperature igneous and metamorphic rocks, which are most likely to preserve former equilibrium mineral assemblages. This preservation is the fundamental prerequisite for any application of phase equilibria methods. For this reason, the phase equilibria approach has generally not been applied to sedimentary rocks, except for saline deposits formed by intense evaporation of seawater and record chemical equilibrium.
The phase diagrams and sophisticated calculations utilized in phase equilibria studies are often imposing, but that merely reflects the compositional complexity of natural rocks and minerals. What must be appreciated is that the goal of such studies is both simple and practical: to determine how rocks form. All processes taking place on or within the Earth (as well as all other solar system bodies) involve the flow of energy and mass. If scientists wish to advance beyond simply describing these processes—that is, to understand the chemical nature of the world—the phase equilibria approach must be employed.
Principal Terms
degree of freedom: the variance of a system; the least number of variables that must be fixed to define the state of a system in equilibrium, generally symbolized by F in the phase rule (P + F = C + 2), where P is the number of phases and C is the number of chemical components
equilibrium: the condition of a system at its lowest energy state compatible with the composition (X), temperature (T), and pressure (P) of the system; the smallest change in T, P, or X induces a state of disequilibrium that the system attempts to rectify
isochemical processes: processes that leave rock compositions unchanged; in thermodynamic terms, a system in which X remains constant even if T and P change
mole: the amount of pure substance that contains as many elementary units as there are atoms in 12 grams of the isotope carbon-12
phase: any part of a system—solid, liquid, or gaseous—that is physically distinct and mechanically separable from other parts of the system; a boundary surface separates adjacent phases
phase diagrams: graphical devices that show the stability limits of rocks or minerals in terms of the variables T, P, and X; the simplest and most widely used are P-T diagrams (X = constant) and T-X diagrams (P = constant)
system: any part of the universe (for example, a crystal, a given volume of rock, or an entire lithospheric plate) that is set aside for thermodynamic analysis; open systems permit energy and mass to enter and leave, while closed systems do not
thermodynamics: the science that treats transformations of heat into mechanical work and the flow of energy and mass from one system to another, based on the assumption that energy can neither be created nor destroyed (the first law of thermodynamics)
Bibliography
Aharonov, Einat. Solid-Fluid Interactions in Porous Media: Processes That Form Rocks. Massachusetts Institute of Technology, 1996.
Angrist, Stanley W., and Loren G. Hepler. Order and Chaos: Laws of Energy and Entropy. Basic Books, 1967.
Asimow, Paul. “A New, Fast Computational Tool for Magmatic Phase Equilibria.” Eos, 11 Aug. 2022, eos.org/editor-highlights/a-new-fast-computational-tool-for-magmatic-phase-equilibria. Accessed 18 May 2026.
Best, Myron G. Igneous and Metamorphic Petrology. 2nd ed., Blackwell Science, 2003.
Blatt, Harvey, et al. Petrology: Igneous, Sedimentary, and Metamorphic. 3rd ed., W. H. Freeman, 2005.
Bowen, Norman L. The Evolution of the Igneous Rocks. Dover, 1956.
Carey, Van P. Liquid Vapor Phase Change Phenomena. 2nd ed., Taylor & Francis, 2007.
Ciccioli, Andrea, and Leslie Glasser. “Complexities of One-Component Phase Diagrams.” Journal of Chemical Education, vol. 88, no. 5, 2011, pp. 586–91.
Ehlers, Ernest G. The Interpretation of Geological Phase Diagrams. W. H. Freeman, 1972.
Emsley, John. The Elements. 3rd ed., Oxford UP, 1998.
Fermi, Enrico. Thermodynamics. Dover, 1956.
Fyfe, W. S., et al. Metamorphic Reactions and Metamorphic Facies. Geological Society of America, 1958.
Greenwood, Norman Neill, and Allen Earnshaw. Chemistry of Elements. 2nd ed., Butterworth-Heinemann, 1997.
Grotzinger, John, et al. Understanding Earth. 5th ed., W. H. Freeman, 2006.
Hillert, Mats. Phase Equilibria, Phase Diagrams and Phase Transformations. 2nd ed., Cambridge UP, 2008.
"Phase Equilibria." LibreTexts, 20 July 2022, chem.libretexts.org/Bookshelves/Physical_and_Theoretical_Chemistry_Textbook_Maps/Physical_Chemistry_(LibreTexts)/23%3A_Phase_Equilibria. Accessed 18 May 2026.
Powell, Roger. Equilibrium Thermodynamics in Petrology: An Introduction. Harper & Row, 1978.
Van Ness, H. C. Understanding Thermodynamics. Dover Publications, 1983.
Full Article
The mineral assemblages in most igneous and metamorphic rocks preserve a record of the chemical equilibrium related to the initial rock-forming process. Phase equilibria studies attempt to determine quantitatively the pressure and temperature conditions of rock formations from these mineral assemblages.
Chemical Equilibrium
The traditional methods of studying rock bodies are descriptive in nature and involve mapping large-scale outcrop features in the field and detailed microscopic observations of rock textures and mineralogy in the laboratory. These methods, which are often successful by themselves, are supplemented by a second, more theoretical approach wherein rocks are treated as chemical systems and the principles of phase equilibria are applied to determine the conditions of their origin. A full and complete description of a rock body is still required, but that is no longer the main goal of petrologic study. The principles of phase equilibria are simply the laws that govern the attainment of equilibrium of chemical reactions such as A + B = C + D, where A and B are known as reactants, and C and D are known as products.
Before exploring how these principles cast light on rock-forming processes, the concept of chemical equilibrium must first be developed. By analogy with gravitational potential, there must exist a similar tendency in chemical systems to lower their energy state through chemical reactions. Reasoning along these lines, J. Willard Gibbs introduced the term “chemical potential” to describe the flow of chemical components from one site (of high potential) to another (of lower potential) during reactions that lead a chemical system toward its lowest energy state. The total energy available to drive a chemical reaction must therefore be the sum of the chemical potentials of each component in the system multiplied by the number of moles of each component. The usual definition of chemical potential of a phase (or pure substance) is “the molar free energy” (or free energy per mole). This simple statement leads to a workable, three-part definition of chemical equilibrium, which is central to the understanding of phase equilibria.
First, if the chemical potential of a component is the same on either side of a reaction equation, the component can have no tendency to participate in the reaction. Second, in a multicomponent system consisting of several phases under uniform temperature (T) and pressure (P), equilibrium must prevail when the chemical potential of each component is the same in all phases in which the component is present. Third, the condition of equilibrium is one of maximum chemical stability. The second part of the definition is equivalent to saying that, for a given chemical reaction, the free energy of the reactants must equal the free energy of the products if a condition of chemical equilibrium prevails under fixed conditions of T and P. If either T or P changes, the system is no longer in equilibrium, but it will immediately adjust itself in such a way as to “moderate” the effect of the disturbing factors. The last statement is known as the moderation theorem, or Le Châtelier's principle. If the free energy on the product side of an equation is less than the free energy on the reactant side, the reaction will be spontaneous. If the opposite is true, no reaction is possible. Scientists are thus able to predict the result of any chemical reaction if the free energy of the reactants and products under the reaction conditions (T, P) is known.
Equilibrium State of a System
The most significant processes of rock formation—magmatism, metamorphism, and sedimentation—all involve large-scale flow of energy and movement of matter that produce an uneven distribution of chemical potential. Inevitably, the result must be chemical reactions tending to restore these natural systems to a state of equilibrium. The equilibrium state of a system is governed by its bulk composition (X), temperature (T), and pressure (P). For most geological processes, T and P change slowly relative to the rates of most chemical reactions, which means that most rock-forming reactions may be considered to take place under constant T and P, and, if X also remains constant, most such reactions should easily attain chemical equilibrium.
In geology, the major concern is not so much with the achievement of equilibrium but rather with the preservation of equilibrium mineral assemblages through hundreds of millions of years, which must follow before deep-seated rocks are finally exposed at the surface. Rocks formed at depth must clearly undergo significant reductions in T and P prior to exposure at the surface, and there are several mechanisms that may induce changes in X during this lengthy period. Geologists are acutely aware of the implications of the moderation theorem: Retrograde metamorphism, mineralogical inversions and exsolutions, hydrothermal alteration, and weathering are but a few of the processes that could trigger reequilibration in rock bodies before they are exposed for study. Fortunately, microscopic studies coupled with Gibbs's pioneering work in phase equilibria provide the means to discern whether a given mineralogical assemblage preserves a former equilibrium.
Gibbs Phase Rule
The “phase rule,” also called the Gibbs phase rule, was initially derived by Gibbs in 1878 from the mathematical formulas of thermodynamics. The phase rule, which is fully applicable to all chemical systems, expresses the relationship between the governing variables T, P, and X, and the number of phases that may coexist in a state of equilibrium. Usually, the phase rule is expressed in equation form as P + F = C + 2, where P = number of phases, C = number of chemical components, and F = degrees of freedom possessed by the system (normally T, P, X). Phases are chemically pure, physically separable subparts of the system and may be gases, liquids, or solids. In the formal sense, C is the minimum number of chemical entities needed to define completely the composition of each reactant and product phase in a given reaction.
Although the objective of phase rule applications is to determine F for major rock-forming reactions, a far simpler situation that could be experimentally verified in any high school laboratory or even in an ordinary kitchen may be considered. Pure water (H2O) boils at T = 100 degrees Celsius at sea level (P = 1 bar; atmospheric pressure). The effect of dissolving common salt (NaCl) in water is to raise the boiling temperature approximately 0.8 degrees Celsius for each mole percent of NaCl in the liquid phase. The steam given off by boiling is pure H2O, and, therefore, the salt concentration in the remaining liquid must progressively increase with temperature during boiling.
To apply the phase rule to this simple system, one first must tally up the participating phases: there is steam and there is liquid salt solution, and one must conclude that P = 2. Both pure water and salt are required to form these coexisting phases and, therefore, C = 2. The phase rule for this process (boiling), under conditions of fixed pressure (P = 1 bar), tells one that
F=C-P+2=2-2+2=2
Therefore, the system has two degrees of freedom and is said to be “divariant,” which means that, because one degree of freedom is utilized by fixing P = 1 bar, only one additional variable needs to be known to specify completely the state of the system. That may be either T (boiling temperature) or X (composition of the boiling solution); in other words, T and X are dependent variables at constant P. This T-X compositional dependence is easily determined for P = 1 bar by direct experiment. If the resulting T-X data were graphically plotted, the diagram would indicate, for example, that the boiling temperature for an 8 mole percent solution is close to 106 degrees Celsius. Conversely, if it were known only that the boiling temperature of a salt solution were 106 degrees Celsius, that would necessitate a solution concentration of 8 mole percent at P = 1 bar.
Application of the Phase Rule
To apply the phase rule to a reaction that has some geological significance, consider the appearance of diopside in siliceous dolomite during contact metamorphism by the reaction
CaMg(CO3)2(dolomite)+2SiO2(quartz)=CaMgSi2O6(diopside)+2CO2(gas)
The reaction involves three mineral phases and a fugitive gas phase, which is necessarily lost from the rock if diopside appears; P = 4. Note that the Ca:Mg (calcium-to-magnesium) ratio is the same in the reactant phase (dolomite) and the product phase (diopside); consequently, the minimum number of components needed to define the compositions of the four phases in the reaction is three. Therefore, C = 3. Substituting these values into the phrase rule, one obtains
F=C−P+2=3−4+2=1
and the reaction has one degree of freedom and is said to be univariant. For any given reaction pressure, the phases appearing in that reaction can coexist in equilibrium at one—and only one—temperature.
Univariant reactions are of great interest to petrologists because it is frequently possible to estimate the depth of a rock-forming reaction and, therefore, P, from field relationships. It is then a simple matter to estimate the reaction temperature from the P-T diagram for the univariant reaction involved. For example, suppose that field relationships lead to the conclusion that diopside-bearing contact metamorphic rocks formed by the reaction above at an estimated depth of 4 kilometers. The pressure equivalent of this depth is about 1,000 bars. The experimentally derived P-T diagram indicates a reaction temperature of about 450 degrees Celsius.
Mineralogical Phase Rule
The principles of phase equilibria were first applied to rocks in 1911 by Swiss geochemist V. M. Goldschmidt in his classic account of contact metamorphism in the Oslo area in Norway. Countless similar attempts have followed, with most authors concluding, as did Goldschmidt, that rocks, in general, record a state of chemical equilibrium governed by the temperature and pressure prevailing at their time of origin. The major generalization emerging from eighty years of such studies is that rocks with a large number of mineral phases tend to have a low number of degrees of freedom. Goldschmidt recognized that, at Oslo, metamorphism must be controlled mainly by T and P and that divariant equilibrium (F = 2) is the general case for isochemical rock-forming reactions. For F = 2, the phase rule reduces to P = C. Dubbed “the mineralogical phase rule,” this equation cannot be mathematically derived from thermodynamic principles. It simply reflects the common case in nature where rock-forming processes are approximately isochemical, and their phase equilibria are controlled by both T and P operating as independent variables. It follows that rocks recording univariant equilibria are more restricted in occurrence than those recording divariant equilibria. Similarly, the rarity of rocks recording invariant equilibria (F = 0) can be understood, as the phase rule requires that both a unique T and a unique P be maintained during their formation. Rocks of complex mineralogy are occasionally encountered in which the number of phases, P, exceeds C + 2 and, consequently, the value of F is negative. Phase rule departures, indicated by negative F values, are a sure sign of disequilibrium and serve as a reminder that not all mineral assemblages can be treated by the methods of phase equilibria.
Petrogenetic Grids
The temperature and pressure range in which any mineral may exist is limited, and it is the task of experimental petrology to ascertain these limits for rock-forming minerals. The data resulting from this experimental work are utilized to construct “phase diagrams.” Such diagrams show the effects of changing P, T, and X values on mineral stability fields and are, therefore, simply graphical expressions of the phase rule. In metamorphic petrology, the major use of phase equilibria data has been for construction of “petrogenetic grids.” A petrogenetic grid, initially conceived by Canadian petrologist N. L. Bowen in 1940, is a P-T diagram on which experimentally derived univariant reaction curves are plotted for a particular metamorphic rock type (for example, blueschists, marbles, calcsilicates, or pelites). The value of such “grids” lies in the fact that each natural equilibrium assemblage recognized in the field will fall within a definite P-T pigeonhole and thereby inform the field geologist immediately of the P-T conditions of the metamorphic terrane under study. This goal, so simple in concept, has proved elusive even after half a century of vigorous experimental, theoretical, and field effort. The problem lies in the lack of truly univariant reactions. For nearly seventy years, field geologists mapped isograds recording the “first appearance” of notable zone minerals such as biotite, garnet, staurolite, kyanite, and sillimanite under the impression that they represented the intersection of the ground surface with a plane of univariant equilibrium. Virtually all such “isograds” have proved to be the result of divariant equilibria and thus plot as a “band”—which may be rather wide—on a P-T diagram. This undesirable result has the effect of “smearing” grid boundaries and rendering them less useful as metamorphic indicators.
The general absence of univariant reactions in metamorphic rocks was eventually recognized because of theoretical and experimental advances in phase equilibria studies. The problem stems from the fact that most mineral phases participating in metamorphic reactions are solid solutions of variable composition, and the most common reactions lead to the release of a fluid, the composition of which may vary with time. Each of these effects introduces an additional degree of freedom in phase rule terms, and, as a result, virtually all important reactions are divariant.
In spite of these difficulties, petrogenetic grids, based on divariant and quasi-univariant equilibria, have gradually evolved for all major metamorphic rock types. These are not the simple, quantitative grids envisioned by Bowen, but they do provide quick, reliable, and fairly narrow estimates of P-T conditions for common metamorphic mineral assemblages. Modern grids, continually subject to refinement, are phase equilibria's greatest contribution to metamorphic petrology.
Since 2022, efficient open-source Gibbs free energy minimization tools such as Mineral Assemblage Gibbs Energy Minimization (MAGEMin) have revolutionized phase equilibria modeling by enabling rapid, high-resolution calculations of stable mineral assemblages and melt properties across complex multicomponent systems. These advances allow seamless integration with large-scale geodynamic simulations. This shift transforms phase equilibria from mainly interpretive tools for ancient rocks into powerful predictive engines for active Earth processes, significantly reducing uncertainties in P-T paths and strengthening quantitative links between petrology, tectonics, and mantle dynamics.
Application to Igneous Petrology
In the area of igneous petrology, phase equilibria methods and data have become indispensable. Natural rocks, spanning the compositional spectrum, are melted under strictly controlled laboratory conditions to determine solidus and liquidus temperatures at pressures ranging from 1 to 35,000 bars. The results of such experiments place tight constraints on the depths and temperatures of magma generation. They also permit the experimentalist to explore P-T effects on partial melting (anatexis) in terms of melt composition and refractory solid phases. The resulting phase diagrams, like metamorphic grids, permit petrologists to “see” deep into the crust and upper mantle and to test hypotheses dealing with the origin of magma.
For nearly a century, igneous petrologists have studied crystal-melt equilibria of simplified, synthetic melts as models for complex, natural magmas. The objective is to reduce the number of equilibrium phases by the elimination of minor components of real magmas. Studies of this type were introduced by Bowen at the Geophysical Laboratory of the Carnegie Institution in Washington, D.C. Through its many subsequent researchers, this laboratory published hundreds of phase diagrams and earned a reputation for meticulous and exhaustive experimental work.
Phase equilibrium studies have provided a rather complete understanding of two fundamentally different modes of magma crystallization. Equilibrium crystallization occurs when P-T-X conditions change so slowly that chemical reactions within the melt are able to maintain the state of chemical equilibrium. Conversely, fractional crystallization results when changes in P-T-X conditions outpace the compensating reactions. This disequilibrium process greatly influences the behavior of natural magmas and extends the range of melt compositions that can be derived from a given parent magma. This latter type of behavior, recognized through the early phase equilibria studies of the Geophysical Laboratory, is the major factor in explaining the compositional diversity of igneous rocks.
The relatively simple phase diagrams of synthetic systems unraveled the complexities of sequential crystallization, cast light on the mechanics of crystal nucleation, and exposed the crucial role that water plays in magmatic processes. Collectively, these diagrams are the foundation of modern igneous petrology.
Goal of Phase Equilibria Studies
The refined symbolic notation and elegant mathematical derivations of thermodynamics are likely to remain unappreciated by the majority of laypersons and geologists alike. It is precisely these formalisms, however, that place phase equilibria on a quantitative footing and permit calculation of mineral stability fields from compositional data. Future development in the area of phase equilibria will follow this theoretical line.
The qualitative form of phase equilibria is expressed in phase diagrams rather than equations. Such diagrams have been a major part of petrology since the 1950s. Historically, emphasis in phase equilibria studies has been on high-temperature igneous and metamorphic rocks, which are most likely to preserve former equilibrium mineral assemblages. This preservation is the fundamental prerequisite for any application of phase equilibria methods. For this reason, the phase equilibria approach has generally not been applied to sedimentary rocks, except for saline deposits formed by intense evaporation of seawater and record chemical equilibrium.
The phase diagrams and sophisticated calculations utilized in phase equilibria studies are often imposing, but that merely reflects the compositional complexity of natural rocks and minerals. What must be appreciated is that the goal of such studies is both simple and practical: to determine how rocks form. All processes taking place on or within the Earth (as well as all other solar system bodies) involve the flow of energy and mass. If scientists wish to advance beyond simply describing these processes—that is, to understand the chemical nature of the world—the phase equilibria approach must be employed.
Principal Terms
degree of freedom: the variance of a system; the least number of variables that must be fixed to define the state of a system in equilibrium, generally symbolized by F in the phase rule (P + F = C + 2), where P is the number of phases and C is the number of chemical components
equilibrium: the condition of a system at its lowest energy state compatible with the composition (X), temperature (T), and pressure (P) of the system; the smallest change in T, P, or X induces a state of disequilibrium that the system attempts to rectify
isochemical processes: processes that leave rock compositions unchanged; in thermodynamic terms, a system in which X remains constant even if T and P change
mole: the amount of pure substance that contains as many elementary units as there are atoms in 12 grams of the isotope carbon-12
phase: any part of a system—solid, liquid, or gaseous—that is physically distinct and mechanically separable from other parts of the system; a boundary surface separates adjacent phases
phase diagrams: graphical devices that show the stability limits of rocks or minerals in terms of the variables T, P, and X; the simplest and most widely used are P-T diagrams (X = constant) and T-X diagrams (P = constant)
system: any part of the universe (for example, a crystal, a given volume of rock, or an entire lithospheric plate) that is set aside for thermodynamic analysis; open systems permit energy and mass to enter and leave, while closed systems do not
thermodynamics: the science that treats transformations of heat into mechanical work and the flow of energy and mass from one system to another, based on the assumption that energy can neither be created nor destroyed (the first law of thermodynamics)
Bibliography
Aharonov, Einat. Solid-Fluid Interactions in Porous Media: Processes That Form Rocks. Massachusetts Institute of Technology, 1996.
Angrist, Stanley W., and Loren G. Hepler. Order and Chaos: Laws of Energy and Entropy. Basic Books, 1967.
Asimow, Paul. “A New, Fast Computational Tool for Magmatic Phase Equilibria.” Eos, 11 Aug. 2022, eos.org/editor-highlights/a-new-fast-computational-tool-for-magmatic-phase-equilibria. Accessed 18 May 2026.
Best, Myron G. Igneous and Metamorphic Petrology. 2nd ed., Blackwell Science, 2003.
Blatt, Harvey, et al. Petrology: Igneous, Sedimentary, and Metamorphic. 3rd ed., W. H. Freeman, 2005.
Bowen, Norman L. The Evolution of the Igneous Rocks. Dover, 1956.
Carey, Van P. Liquid Vapor Phase Change Phenomena. 2nd ed., Taylor & Francis, 2007.
Ciccioli, Andrea, and Leslie Glasser. “Complexities of One-Component Phase Diagrams.” Journal of Chemical Education, vol. 88, no. 5, 2011, pp. 586–91.
Ehlers, Ernest G. The Interpretation of Geological Phase Diagrams. W. H. Freeman, 1972.
Emsley, John. The Elements. 3rd ed., Oxford UP, 1998.
Fermi, Enrico. Thermodynamics. Dover, 1956.
Fyfe, W. S., et al. Metamorphic Reactions and Metamorphic Facies. Geological Society of America, 1958.
Greenwood, Norman Neill, and Allen Earnshaw. Chemistry of Elements. 2nd ed., Butterworth-Heinemann, 1997.
Grotzinger, John, et al. Understanding Earth. 5th ed., W. H. Freeman, 2006.
Hillert, Mats. Phase Equilibria, Phase Diagrams and Phase Transformations. 2nd ed., Cambridge UP, 2008.
"Phase Equilibria." LibreTexts, 20 July 2022, chem.libretexts.org/Bookshelves/Physical_and_Theoretical_Chemistry_Textbook_Maps/Physical_Chemistry_(LibreTexts)/23%3A_Phase_Equilibria. Accessed 18 May 2026.
Powell, Roger. Equilibrium Thermodynamics in Petrology: An Introduction. Harper & Row, 1978.
Van Ness, H. C. Understanding Thermodynamics. Dover Publications, 1983.
More Like ThisRelated Articles
Related Articles (5)
Related Articles (5)
- A Barrovian Overprint Superimposed on Partly Preserved High-Pressure Metamorphic Rocks: Pressure–Temperature–Time Evolution in the Ordovician Corner Brook Collisional Complex (Western Newfoundland, Canada).Published In: Journal of Petrology, 2025, v. 66, n. 4. P. 1Authored By: Bissbort, Thilo; Willner, Arne P; Glodny, Johannes; Jöns, Niels; Staal, Cees R van; Schertl, Hans-PeterPublication Type: Academic Journal
- Advances in Igneous Petrology: Coupling Chemical Maps with Thermodynamic Models to Tackle Mush Crystallization Dynamics and Crystal–Melt Segregation.Published In: Journal of Petrology, 2025, v. 66, n. 5. P. 1Authored By: France, Lydéric; Charvet, Valentine; Toussaint, AurorePublication Type: Academic Journal
- Age of Skarn Formation and Meteoric Water Infiltration in the Willsboro–Lewis Wollastonite District (Adirondack Highlands, New York).Published In: Journal of Petrology, 2025, v. 66, n. 4. P. 1Authored By: Peck, William H; Regan, Sean P; Blum, Tyler B; Valley, John W; Timothy, Samuel CPublication Type: Academic Journal
- Identification of Neoarchean UHT (Ultrahigh‐Temperature) Metamorphism From Intermediate Migmatitic Granulites in the Yinshan Block, North China Craton.Published In: Journal of Metamorphic Geology, 2025, v. 43, n. 6. P. 591Authored By: Yang, Zizhen; Wei, Chunjing; Wang, Bin; Song, Shuguang; TIAN, WeiPublication Type: Academic Journal
- Influence of the mineral gangue on pellets softening melting properties in blast furnace. Experimental study of phase equilibria during the melting of pre-reduced pellets.Published In: Metallurgical Research & Technology, 2026, v. 123, n. 2. P. 1Authored By: Graz, Yann; Athoy, Evangeline; Oudich, Fayssal; Maurice, Yann; Husson, Aurore; Berthelemy, Nadine; Jactard, Jeremy; Basselin, Yves; Siboni, GeraldinePublication Type: Academic Journal