arxiv: v1 [math-ph] 21 Dec 2007

Transcription

1 Dynamic Phase ransitions in PV Systems ian Ma Deartment of Mathematics, Sichuan Uniersity, Chengdu, P. R. China Shouhong Wang Deartment of Mathematics, Indiana Uniersity, Bloomington, IN 4745 (Dated: February 4, 28) arxi: [math-h] 21 Dec 27 he main objectie of this article are two-fold. First, we introduce some general rinciles on hase transition dynamics, including a new dynamic transition classification scheme, and a Ginzburg- Landau theory for modeling equilibrium hase transitions. Second, aly the general rinciles and the recently deeloed dynamic transition theory to study dynamic hase transitions of PV systems. In articular, we establish a new time-deendent Ginzburg-Landau model, whose dynamic transition analysis is carried out. It is worth ointing out that the new dynamic transition theory, along with the dynamic classification scheme and new time-deendent Ginzburg Landau models for equilibrium hase transitions can be used in other hase transition roblems, including e.g. the ferromagnetism and suerfluidity, which will be reorted elsewhere. In addition, the analysis for the PV system in this article leads to a few hysical redications, which are otherwise unclear from the hysical oint of iew. PACS numbers: I. INRODUCION he main objecties of this article are 1) to introduce a new dynamic hase transition classification scheme for nonlinear systems, 2) to resent a new time-deendent Ginzburg-Landau model for general equilibrium hase transitions, and 3) to study the dynamic hase transitions for PV systems. As we know, there are seeral classification schemes for hase transitions. he one used most often is the classical Ehrenfest classification scheme, which is based on the order of differentiability of the energy functional. he new dynamic hase transition classification scheme we roose here is based on the dynamic transition theory deeloed recently by the authors [1, 2], which has been used to study numerous roblems in sciences and engineering, including in articular roblems in classical and geohysical fluid dynamics, biology and chemistry, and hase transitions. he new scheme introduced in this article classifies hase transitions into three categories: ye-i, ye-ii and ye-iii, corresonding resectiely to the continuous, the jum and mixed transitions in the mathematical dynamic transition theory. We shall see from the alications to PV systems that this new classification scheme is more transarent. In addition, it can be used to identify high-order (in the he work was suorted in art by the Office of Naal Research and by the National Science Foundation. URL: htt:// fluid Ehrenfest sense) transitions, which are usually hard to derie both theoretically and exerimentally. Second, based on the le Châtelier rincile and some general characteristics of seudo-gradient systems, we introduce a general rincile, leading to a unified aroach to derie Ginzburg-Landau tye of time-deendent models for equilibrium hase transitions. In iew of this new dynamic modeling and the dynamic transition theory, we shall see that the states after transition often include not only equilibria, but also some transient states, which are hysically imortant as well. hird, with the aforementioned rinciles in our disosal, we derie a new dynamic model for PV system, whose dynamic transition can be exlicitly obtained, leading to some secific hysical redictions. As we know, a P V system is a system comosed of one tye of molecules, and the interaction between molecules is goerned by the an der Waals law. he molecules generally hae a reulsie core and a shortrange attraction region outside the core. Such systems hae a number of hases: gas, liquid and solid, and a solid can aear in a few hases. he most tyical examle of a P V system is water. In general, the hase transitions of P V systems mainly refer to the gas-liquid, gas-solid and liquid-solid hase transitions. hese transitions are all first order in the Ehrenfest sense (i.e., discontinuous) and are accomanied by a latent heat and a change in density. A P -hase diagram of a tyical P V system is schematically illustrated by Figure 1, where oint A is the trile oint at which the gas, liquid, and solid hases can coexist. Point C is the Andrews critical oint at which the gas-liquid coexistence cure terminates [3]. he fact that the gas-liquid coexistence cure has a critical oint means that we can go continuously from a gaseous state to a liquid state without eer meeting an

2 2 obserable hase transition, if we choose the right ath. he critical oint in the liquid-solid coexistence cure has neer been obsered. It imlies that we must go through a first-order hase transition in going from the liquid to the solid state. solid A B liquid gas C In addition, in the gas-solid and liquid-solid transitions, there also exist metastable states. For the gas-solid case, the metastable states corresond to the suerheated solid and suercooled liquid, and for the liquid-solid case, the metastable states corresond to the suerheated solid and suercooled liquid; see Physical Conclusion V.2. his article is organized as follows. he new dynamic classification scheme and the Ginzburg-Landau model for equilibrium hase transitions are gien in Section II. Section III introduces a time-deendent model for PV systems, which is analyzed in Section IV. Section V roides hysical consequences of the analysis and some new hysical redictions.he abstract dynamical transition theory is recaitulated in the Aendix A for conenience. FIG. 1: Coexistence cures of a tyical P V system: A is the trile oint, C is the critical oint, and the dashed cure is a melting line with negatie sloe. II. GENERAL PRINCIPLES OF PHASE RANSIION DYNAMICS he hase diagram in the P V -lane is gien by Figure 2, where the dashed cures reresent lines of constant temerature. In the region of coexistence of hases, the isotherms (dashed lines) are always flat, indicating that in these regions the change in olume (density) occurs for constant ressure and temerature. B solid liquid - solid A B gas - liquid A gas - solid C A = < gas > FIG. 2: P V -hase diagram, the dashed line reresent isotherms. One of the main objecties of this article is to study the dynamic hase transitions beyond the Andrews oint C, and to make some hysical redications. For this urose, first, by both mathematical and hysical considerations, a roer form of the Gibbs free energy functional is gien, and then the time-deendent Ginzburg-Landau tye of model for PV systems is deried naturally using the general rincile as mentioned aboe. hen, using the dynamical transition theory, we derie a system of critical arameter equations to determine the Andrews oint C = ( C, C ) gien by (V.6). hen we show that the gas-liquid transition can be either firstorder, second order or third order, as demonstrated in Physical Conclusion V.1, if we go from a gas to a liquid by a ath without assing through the gas-liquid coexistence cure. c c c V In this section, we introduce a new hase dynamic transition classification scheme to classify hase transitions into three categories: ye-i, ye-ii and ye-iii, corresonding mathematically continuous, jum mixed transitions, resectiely. hen a new time-deendent Ginzburg-Landau theory for modeling equilibrium hase transitions in statistical hysics is deried based on the le Châtelier rincile and some mathematical insights on seudo-gradient systems. A. Dynamic Classification Scheme In sciences, nonlinear dissiatie systems are generally goerned by differential equations, which can be exressed in the following abstract form du dt = L u + G(u, ), u() = ϕ, (II.1) where u : [, ) X is the unknown function, R N (N 1) is the control arameter, X and X 1 are two Banach saces with X 1 X being a dense and comact inclusion, L = A + B and G(, ) : X 1 X are C r (r ) maings deending continuously on, L : X 1 X is a sectorial oerator, and A : X 1 X B : X 1 X a linear homeomorhism, a linear comact oerator. (II.2) In following, we introduce some basic and uniersal concets in nonlinear sciences. First, a state of the system (II.1) at is usually referred to as a comact inariant set Σ. In many alications, Σ is a singular oint or a eriodic orbit. A state Σ of (II.1) is stable if Σ is an attractor, otherwise Σ is called unstable. Second, we say that the system (II.1) has a hase transition from a state Σ at = if Σ is stable on <

3 3 (or on > ) and is unstable on > (or on < ). he critical arameter is called a critical oint. In other words, the hase transition corresonds to an exchange of stable states. he concet of hase transition originates from the statistical hysics and thermodynamics. In hysics and chemistry, hase means the homogeneous art in a heterogeneous system. Howeer, here the so called hase means the stable state in the systems of nonlinear sciences including hysics, chemistry, biology, ecology, economics, fluid dynamics and geohysical fluid dynamics, etc. Hence, here the content of hase transition has been endowed with more general significance. In fact, the hase transition dynamics introduced here can be alied to a wide ariety of toics inoling the uniersal critical henomena of state changes in nature in a unified mathematical iewoint and manner. hird, if the system (II.1) ossesses the gradient-tye structure, then the hase transitions are called equilibrium hase transition; otherwise they are called the nonequilibrium hase transitions. Fourth, classically, there are seeral ways to classify hase transitions. he one most used is the Ehrenfest classification scheme, which grous hase transitions based on the degree of non-analyticity inoled. First order hase transitions are also called discontinuous, and higher order hase transitions (n > 1) are called continuous. Here we introduce the following notion of dynamic classification scheme: Definition II.1 Let R N be a critical oint of (II.1), and (II.1) undergo a transition from state Σ 1 to Σ2. here are three tyes of hase transitions for (II.1) at =, deending on their dynamic roerties: continuous, jum, and mixed as gien in heorem A.1, which are called ye-i, ye-ii and ye-iii resectiely. he main characteristics of ye-ii hase transitions is that there is a ga between Σ 1 and Σ2 at the critical oint. In thermodynamics, the metastable states corresond in general to the suer-heated or suer-cooled states, which hae been found in many hysical henomena. In articular, ye-ii hase transitions are always accomanied with the latent heat to occur. In a ye-i hase transition, two states Σ 1 and Σ2 meet at, i.e., the system undergoes a continuous transition from Σ 1 to Σ2. In a ye-iii hase transition, there are at least two different stable states Σ 2 and Σ 3 at, and system undergoes a continuous transition to Σ 2 or a jum transition to Σ 3, deending on the fluctuations. It is clear that a ye-ii hase transition of gradienttye systems must be discontinuous or the zero order because there is a ga between Σ 1 and Σj (2 j K). For a ye-i hase transition, the energy is continuous, and consequently, it is an n-th order transition in the Ehrenfest sense for some n 2. A ye-iii hase transition is indefinite, for the transition from Σ 1 to Σ2 it may be continuous, i.e., df /d = df 2 + /d, and for the transitions from Σ 1 to Σ3 it may be discontinuous: df /d df 3 + /d at =. Finally, we recall the hysical meaning of the deriaties of F (u, ) on in statistical hysics. In thermodynamics, the arameter generally stands for temerature and ressure = (, ). he energy functional F (u,, ) is the thermodynamic otential, and u is the order arameter. he first order deriaties are F = S F = V the entroy, and the second order deriaties are the hase olume, 2 F 2 = C the heat caacity in constant ressure, 1 2 F = κ the comression coefficient, V 2 1 V 2 F = α the thermal exansion coefficient, for the P V system. hus, for the first order hase transition, the discontinuity of F/ at critical oint = (, ) imlies the discontinuity of hase olume: V = V 2 V 1 = F + F, (II.3) and the discontinuity of F/ imlies that there is a ga between both hase entroies: S = S 2 S 1 = F + + F (II.4) Physically, S is a non-measurable quantity, hence, we always use latent heat H = S to determine the first order hase transition, and H stand for the absorbing heat for transition from hase Σ 1 to hase Σ2. In the two hase coexistence situation, H and V are related by the Clausius-Claeyron equation d d = H V. For second order hase transitions the ariance of heat caacity (or secific heat), comression coefficient, and the thermal exansion coefficient at critical alue are imortant static roerties in equilibrium hase transition, which are measurable hysical quantities.

4 4 B. New Ginzburg-Landau Models for Equilibrium Phase ransitions In this subsection, we introduce the time-deendent Ginzburg-Landau model for equilibrium hase transitions. We start with thermodynamic otentials and the Ginzburg-Landau free energy. As we know, four thermodynamic otentials internal energy, the enthaly, the Helmholtz free energy and the Gibbs free energy are useful in the chemical thermodynamics of reactions and noncyclic rocesses. Consider a thermal system, its order arameter u changes in Ω R n (1 n 3). In this situation, the free energy of this system is of the form [ 1 H(u, ) = H + Ω 2 m i=1 ] µ i u i 2 + g(u, u, ) dx (II.5) where N 3 is an integer, u = (u 1,, u m ), µ i = µ i () >, and g(u, u, ) is a C r (r 2) function of (u, u) with the aylor exansion g(u, u, ) = α ijk u i D j u k + N α I u I +o( u N ) fx, I =1 (II.6) where I = (i 1,, i m ), i k are integer, I = m k=1 i k, the coefficients α ijk and α I continuously deend on, which are determined by the concrete hysical roblem, u I = u i1 1 uim m and fx the generalized work. hus, the study of thermal equilibrium hase transition for the static situation is referred to the steady state bifurcation of the system of ellitic equations δ H(u, ) =, δu u n Ω =, (or u Ω = ), where δ/δu is the ariational deriatie. A thermal system is controled by some arameter. When is for from the critical oint the system lies on a stable equilibrium state Σ 1, and when reaches or exceeds the state Σ 1 becomes unstable, and meanwhile the system will undergo a transition from Σ 1 to another stable state Σ 2. he basic rincile is that there often exists fluctuations in the system leading to a deiation from the equilibrium states, and the hase transition rocess is a dynamical behaior, which should be described by a time-deendent equation. o derie a general time-deendent model, first we recall that the classical le Châtelier rincile amounts to saying that for a stable equilibrium state of a system Σ, when the system deiates from Σ by a small erturbation or fluctuation, there will be a resuming force to retore this system to return to the stable state Σ. Second, we know that a stable equilibrium state of a thermal system must be the minimal alue oint of the thermodynamic otential. By the mathematical characterization of gradient systems and the le Châtelier rincile, for a system with thermodynamic otential H(u, ), the goerning equations are essentially determined by the functional H(u, ). When the order arameters (u 1,, u m ) are nonconsered ariables, i.e., the integers u i (x, t)dx = a i (t) constant. Ω then the time-deendent equations are gien by u i t = β δ i H(u, ) + Φ i (u, u, ), δu i u n Ω = (or u Ω = ), u(x, ) = ϕ(x), (II.7) for any 1 i m, where δ/δu i are the ariational deriatie, β i > and Φ i satisfy δ Φ i H(u, )dx =. (II.8) δu i Ω i he condition (II.8) is required by the Le Châtelier rincile. In the concrete roblem, the terms Φ i can be determined by hysical laws and (II.8). When the order arameters are the number density and the system has no material exchange with the external, then u j (1 j m) are consered, i.e., u j (x, t)dx = constant. (II.9) Ω his conseration law requires a continuous equation u j t = J j(u, ), (II.1) where J j (u, ) is the flux of comonent u j. In addition, J j satisfy J j = k j (µ j i j µ i ), (II.11) where µ l is the chemical otential of comonent u l, µ j i j µ i = δ δu j H(u, ) φ j (u, u, ), (II.12) and φ j (u, ) is a function deending on the other comonents u i (i j). When m = 1, i.e., the system consists of two comonents A and B, this term φ j =. hus, from (II.1)-(II.12) we obtain the dynamical equations as follows [ ] u j δ t = β j H(u, ) φ j (u, u, ), δu j u n u Ω =, n (II.13) Ω =, u(x, ) = ϕ(x),

5 5 for 1 j m, where β j > are constants, φ j satisfy δ φ j H(u, )dx =. (II.14) δu j Ω j If the order arameters (u 1,, u k ) are couled to the consered ariables (u k+1,, u m ), then the dynamical equations are u i t = β δ i H(u, ) + Φ i (u, u, ), δu [ i ] u j δ t = β j H(u, ) φ j (u, u, ), δu j u i n Ω = (or u i Ω = ), u j n u j Ω =, n Ω =, u(x, ) = ϕ(x). (II.15) for 1 i k and k + 1 j m. he model (II.15) gies a general form of the goerning equations to thermodynamic hase transitions. Hence, the dynamics of equilibrium hase transition in statistic hysics is based on the new Ginzburg-Landau formulation (II.15). Physically, the initial alue condition u() = ϕ in (II.15) stands for the fluctuation of system or erturbation from the external. Hence, ϕ is generally small. Howeer, we can not exclude the ossibility of a bigger noise ϕ. From conditions (II.8) and (II.14) it follows that a steady state solution u of (II.15) satisfies Φ i (u, u, ) = 1 i k, φ j (u, u, ) = k + 1 j m. (II.16) Hence a stable equilibrium state must reach the minimal alue of thermodynamic otential. In fact, u fulfills β i δ H(u, ) Φ i (u, u, ) =, δu i β j δ δu j H(u, ) φ j (u, u, ) =, u i n Ω = (or u i Ω = ), (II.17) u j n u j Ω =, n Ω =, for 1 i k and k + 1 j m. Multilying Φ i (u, u, ) and φ j (u, u, ) on the first and the second equations of (II.17) resectiely, and integrating them, then we infer from (II.8) and (II.14) that Φ 2 i (u, u, )dx =, Ω Ω i φ j (u, u, ) 2 dx =, j which imly that (II.16) holds true. III. IME-DEPENDEN MODEL FOR P V SYSEMS In this section, we use the general rinciles deried in the last section to derie a time-deendent Ginzburg- Landau model for PV systems, which will be used to carry out dynamic transition analysis and hysical redications in the next section. A. an der Waals equations and Gibbs energy he classical and the simlest equation of state which can exhibit many of the essential features of the gas-liquid hase transition is the an der Waals equation. It reads 3 ( b + R ) 2 + a ab =, (III.1) where is the molar olume, is the ressure, is the temerature, R is the uniersal gas constant, b is the reised constant of inherent olume, and a is the reised constant of attractie force between molecules. If we adot the molar density = 1/ to relace in (III.1), then the an der Waals equation becomes (b + R ) + a 2 ab 3 + =. (III.2) Now, we shall aly thermodynamic otentials to inestigate the hase transitions of P V systems, and we shall see later that the an der Waals equation can be deried as a Euler-Langrange equation for the minimizers of the Gibbs free energy for P V systems at gaseous states. Consider an isothermal-isoiestic rocess. he thermodynamic otential is taken to be the Gibbs free energy. In this case, the order arameters are the molar density and the entroy S, and the control arameters are the ressure and temerature. he general form of the Gibbs free energy for P V systems is gien as [ 1 G(, S,, ) = Ω 2 µ µ 2 S 2 + g(, S,, ) ] S α(,, ) dx, (III.3) where g and α are differentiable with resect to and S, Ω R 3 is the container, and α is the mechanical couling term in the Gibbs free energy, which can be exressed by α(,, ) = 1 2 b2, (III.4) where b = b(, ) deends continuously on and. Based on both the hysical and mathematical considerations, we take the aylor exansion of g(, S,, ) on and S as follows g = 1 2 α β2 S + β 2 S α α 3 4, (III.5)

6 6 where α i (1 i 3), β 1 and β 2 deend continuously on and, and α i = α i (, ) > i = 2, 3, β 1 = β 1 (, ) >. B. Dynamical equations for PV systems (III.6) In a P V system, the order arameter is u = (, S), = 1, S = S 1 S, where i and S i (i =, 1) reresent the density and entroy,, S are reference oints. Hence the conjugate ariables of and S are the ressure and the temerature. hus, by the standard model (II.7), we derie from (III.3)-(III.5) the following general form of the timedeendent equations goerning a P V system: t = µ 1 (α 1 + b) + α 2 2 α 3 3 2β 2 S +, (III.7) S t = µ 2 S β 1 S β Although the domain Ω deends on and, we can still take the Neumann boundary condition n Ω =, S n =. Ω (III.8) An imortant secial case for P V systems is that the ressure and temerature functions are homogeneous in Ω. hus we can assume that and S are indeendent of x Ω, and the free energy (III.3) with (III.4) and (III.5) can be exressed as G(, S,, ) = 1 2 α β 1S 2 + β 2 S α α b2 S. (III.9) From (III.9) we get the dynamical equations as d dt = (α 1 + b) + α 2 2 α 3 3 2β 2 S +, (III.1) ds dt = β 1S β Because β 1 > for all and, we can relace the second equation of (III.1) by S = β 1 1 ( β 2 2 ). (III.11) hen, (III.1) are equialent to the following equation d dt = (α 1 + b + 2β 1 1 β 2 ) + α 2 2 (α 3 2β 2 2β 1 1 )3 +. (III.12) It is clear that if α 1 =, 2β1 1 β 2 = R, α 2 = a, (α 3 2β 2β ) = ab, then the steady state equation of (III.12) is referred to the an der Waals equation. IV. PHASE RANSIION DYNAMICS FOR PV SYSEMS In this section we use (III.12) to discuss dynamical roerties of transitions for P V systems, and remark that similar results can also deried using the more general form of the time-deendent model (III.7). Let be a steady state solution of (III.12). We take the transformation = +. hen equation (III.12) becomes (dro the rime) where d dt = + a 2 2 a 3 3. = 2α 2 3a 3 2 α 1 b 2β 2 β 1 1, a 2 = α 2 3a 3, a 3 = α 3 2β 2 2β 1 1. (IV.1) In the P -lane, near a non-trile oint A = (, ) (see Figure 1), the critical arameter equation = (, ) = in < δ, < δ for some δ >, defines a continuous function = φ(), such that < if > φ(), = if = φ(), (IV.2) > if < φ(). he main result in this section is the following dynamic transition theorem for PV systems. heorem IV.1 Let = φ( ) and a 3 >. hen the system (IV.1) has a transition at (, ) = (, ), and the following assertions hold true: 1. If the coefficient a 2 = a 2 (, ) in (IV.1) is zero at (, ), i.e., a 2 (, ) =, then the transition is of ye-i, as schematically shown in Figure 3 and. 2. If a 2 (, ), then the transition is of ye-iii (i.e., the mixed tye), and the following assertions hold true: here are two transition solutions near (, ) as ± (, ) = 1 ( ) a 2 ± a 2 2 2a + 4a 3. (IV.3) 3 here is a saddle-node bifurcation at ( 1, 1 ), where 1 > and 1 < if φ ( ) >, and > 1 if φ ( ) <.

7 7 (c) For φ ( ) >, when a 2 (, ) > the transition diagrams are illustrated by Figure 4 -, where + is stable for all (, ) near (, ), and = is stable, a saddle for < < 1 and 1 < <, and is stable, = a saddle for <, > ; (d) When a 2 (, ) <, the transition diagrams are illustrated by Figure 5-, where is stable for all (, ) near (, ), and = is stable, + a saddle for < < 1, 1 < <, and + is stable, = a saddle for <, <. _ + _ 1 + _ 1 _ FIG. 4: ye-iii (mixed) transition for a 2 >. + 1 _ FIG. 3: Continuous transition for the case where φ ( ) > + Proof. By (IV.1), the theorem follows from the transition theorems (heorems A.2 and A.3) and the singularity searation theorem (heorem A.4). We omit the detailed routine analysis. 1 _ V. PHYSICAL CONCLUSIONS AND PREDICIONS o discuss the hysical significance of heorem IV.1, we recall the classical V -hase diagram gien by Figure 2. We take = 1/V to relace olume V, then the gasliquid coexistence cure in the -P lane is illustrated by Figure 6. According to the hysical exeriments, in the gasliquid coexistence region in the hase diagram, there exist metastable states. Mathematically seaking, the metastable states are the attractors which hae a small basin of attraction. In Figure 6, the dashed lines aa and bb reresent the metastable states, and the oints FIG. 5: ye-iii (mixed) transition for a 2 <. in aa corresond to suer-heated liquid, the oints in bb corresond to suer-cooled gas. he hase diagram shows that along the isothermal line > C where C = ( C, P C ) is the Andrews critical oint, when the ressure increases, the density aries continuously from gaseous to liquid states. Howeer along the isothermal line < C when the ressure increases to the density will undergo an abrut change, and a transition from gaseous state a to a liquid state b accom-

8 8 liquid b b < c C a = c > c Second, utting S = β 1 1 ( β 2 2 ) into (III.9), the free energy becomes G(,, ) = 1 2 β (α 1 + b + 2β1 1 β 2 ) α (α 3 2β 2 2β 1 1 )4. gas - liquid a gas he alues of G at equilibrium states of (IV.1) are gien by FIG. 6: he gas-liquid coexistence cure in lane: oint C is the Andrews critical oint, the dashed lines aa and bb stand for the metastable states. anied with an isothermal exothermal rocess to occur. Likewise, such rocesses also occur in the gas-solid and liquid-solid transitions. We now return to discuss heorem IV.1. he steady state solution = (, ) of (III.12) can be taken to reresent a desired state in inestigating different transition situation. For examle, for studying the gas-liquid transition we can take as the gas density near the transition temerature and ressure, and for the liquid-solid transition we take as the liquid density, i.e., the lower state density. For conenience, in the following we alway consider the gas-liquid transition, and take as the gas density. Let = φ( ), and φ ( ) >. hen heorem IV.1 imlies the following hysical conclusions: First, near (, ), there are three stable equilibrium states of (III.12) for each of the three cases: a 2 >, a 2 =, or a 2 < : a) If a 2 <, they are ϕ + = + + for < δ & < δ, ϕ = for < or <, ϕ = + for < & > ; b) If a 2 =, they are ϕ ± = + ± for < & >, ϕ = for > or <, where ± = ± /a 3 ; c) If a 2 <, they are ϕ + = + + for < & >, ϕ = for > or <, ϕ = + for < δ & < δ. (V.1) (V.2) (V.3) Here the state ϕ + reresents the liquid density, ϕ the real gas density, and ϕ the underlying state density. G( + ± ) =G( ) 2 (± ) 2 a 2 3 (± ) 3 + a 3 4 (± ) 4 =G( ) 4 (± ) 2 a 2 12 (± ) 3. (V.4) By = at (, ), it follows from (V.1)-(V.4) that for (, ) near (, ) we hae G(ϕ + ) < G(ϕ ) < G(ϕ ) for a 2 >, G(ϕ + ) = G(ϕ ) < G(ϕ ) for a 2 =, G(ϕ ) < G(ϕ ) < G(ϕ + ) for a 2 <. (V.5) hird, according to heorem IV.1 and the transition diagrams (Figures 3-5), there are two tyes of transition behaiors characterized by some functions Φ = (, ) near (, ) = (, ), which are called the transition functions. Here, for simlicity we fixed = and consider Φ = ( ) as function of. If a 2 >, the transition functions are { ( ) for >, Φ + ( ) = ( ) + + ( ) for <, { ( ) for Φ, ( ) = ( ) + for <, for some < 1. Φ + ( ) has a finite jum at = as shown in Figure 7, and Φ ( ) is continuous, but has a discontinuous deriatie at = ; see Figure 7. If a 2 =, the functions are { ( ) for Φ +, ( ) = ( ) + + ( ) for <, { ( ) for Φ, ( ) = ( ) + ( ) for <, and Φ + and Φ are continuous with discontinuous deriaties at = ; see Figure 8 and. If a 2 <, the functions are { ( ) for Φ +, ( ) = ( ) + + ( ) for <, { ( ) for >, Φ ( ) = ( ) + ( ) for <,

9 9 for some < 1, Φ + is continuous with a discontinuous deriatie at = as shown in Figure 9, and Φ is discontinuous at = as shown in Figure Fourth, based on hysical facts, near the gas-liquid transition oint (, ) the density is a decreasing function of, therefore the transition rocess Φ + is realistic, and Φ is unrealistic. From the transition function Φ + ( ), the gas-liquid transition can be well understood. + + * FIG. 8: he transition functions for a 2 = : the cure of = Φ +, the cure of = Φ. + _ + + FIG. 7: he transition functions for a 2 > : the cure of = Φ +, the cure of = Φ. Fifth, we now exlain the gas-liquid transition by using the transition function Φ + ( ). Let the two cures = (, ) = and a 2 = a 2 (, ) = (V.6) intersect at C = ( C, C ) in P -lane (see Figure 1), and the cure segment AB of = is diided into two arts AC and CB by the oint C such that a 2 (, ) > for (, ) AC, a 2 (, ) < for (, ) CB. We shall see that the oint C = ( C, C ) is the Andrews critical oint, and the cure segment AC is the gas-liquid coexistence cure, as shown in Figure 1. In fact, when (, ) AC, a 2 (, ) > and the transition function Φ + ( ) has a jum at = (see Figure 6): + +, + a 2 a 3 >. Hence the system undergoes a transition from a gaseous state to a liquid state with an abrut change in density. + _ * FIG. 9: he transition functions for a 2 < : the cure of = Φ +, the cure of = Φ. On the other hand, by (V.5) there is an energy ga between the gaseous and liquid states: E = G( + + ) G( ) < at =. his energy ga E stands for a latent heat, and E < shows that the transition from a gaseous state to a liquid state is an isothermal exothermal rocess, and from a liquid state to gaseous state is an isothermal endothermal rocess. When (, ) = ( C, C ), a 2 ( C, C ) = and the

10 be obsered by hysical exeriments. herefore we can derie the following hysical conclusion: 1 A = C FIG. 1: B a2 = transition function Φ + ( ) is continuous as shown in Figure 8. Near = C, Φ + ( ) = ( ) + /a 3 for < C. By the Landau mean field theory, we hae ( ) = α( C ) (α > a constant). (V.7) hus we infer from (V.4) and (V.7) that G(Φ + ( )) = G( ) α2 4a 3 ( C ) 2 for < C. he difference of the heat caacity at = C is C = C 2 2 ( G(Φ + ( C + )) G( ) ) = α2 2a 3 C >. Namely the heat caacity has a finite jum at = C, therefore the transition at = C is of the second order. When (, ) CB, a 2 (, ) < and the transition function Φ + ( ) is continuous as shown in Figure 9. Near =, Φ + ( ) = ( ) + 1 a 2 ( ) + o( ). Notice that (IV.3) imlies that + = a 2 + o( ). hen we deduce from (V.4) and (V.7) that if >, and if, G(Φ + ( )) = G( ( )) G(Φ + ( )) = G( ( )), α3 6 a 2 2 ( ) 3 + o( 3 ). Namely, the free energy G(Φ + ( )) is continuously differentiable u to the second order at =, and the transition is of the third order. It imlies that as (, ) CB the ye-i transition at (, ) can not Physical Conclusion V.1.he oint C = ( C, C ) satisfying (V.6) corresonds to the Andrews critical oint, at which the gas-liquid transition is of the second order. In addition, when (, ) =, a 2 (, ) >, the gas-liquid transition at (, ) is of the first order accomanied with a latent heat to occur, and when (, ) =, a 2 (, ) <, the transition at (, ) is of the third order. Sixth, as (, ) = and a 2 (, ) >, the gasliquid transition oint (, ) is in the range < 1 and 1 < ; see Figure 4-. In fact, in the region < 1 and 1 < <, the two stable states ϕ = (, ) and ϕ + = (, ) + + (, ) are attractors, each ossessing a small basin of attraction. herefore they corresond to metastable states, and ϕ can be considered as a suer cooled gas, while ϕ + can be considered a suer heated liquid. Finally, likewise, we can also discuss the gas-solid and liquid-solid transitions, and derie the following hysical conclusion: Physical Conclusion V.2. In the gas-solid and liquidsolid transitions, there also exist metastable states. For the gas-solid case, the metastable states corresond to the suerheated solid and suercooled liquid, and for the liquid-solid case, the metastable states corresond to the suerheated solid and suercooled liquid. APPENDIX A: RECAPIULAION OF HE DYNAMIC RANSIION HEORY OF NONLINEAR SYSEMS In this aendix we recall some basic elements of the dynamic transition theory deeloed by the authors [1, 2], which are used to carry out the dynamic transition analysis for the PV systems in this article. Let X and X 1 be two Banach saces, and X 1 X a comact and dense inclusion. In this chater, we always consider the following nonlinear eolution equations du dt = L u + G(u, ), u() = ϕ, (A.1) where u : [, ) X is unknown function, and R 1 is the system arameter. Assume that L : X 1 X is a arameterized linear comletely continuous field deending contiguously on R 1, which satisfies L = A + B A : X 1 X B : X 1 X a sectorial oerator, a linear homeomorhism, a linear comact oerator. (A.2)

11 11 In this case, we can define the fractional order saces X σ for σ R 1. hen we also assume that G(, ) : X α X is C r (r 1) bounded maing for some α < 1, deending continuously on R 1, and G(u, ) = o( u Xα ), R 1. (A.3) Hereafter we always assume the conditions (A.2) and (A.3), which reresent that the system (A.1) has a dissiatie structure. In the following we introduce the definition of transitions for (A.1). Definition A.1 We say that the system (A.1) has a transition of equilibrium from (u, ) = (, ) on > (or < ) if the following two conditions are satisfied: (1) when < (or > ), u = is locally asymtotically stable for (A.1); and (2) when > (or < ), there exists a neighborhood U X of u = indeendent of, such that for any ϕ U \ Γ the solution u (t, ϕ) of (A.1) satisfies that lim su u (t, ϕ) X δ() >, t lim δ(), where Γ is the stable manifold of u =, with codim Γ 1 in X for > (or < ). Obiously, the attractor bifurcation of (A.1) is a tye of transition. Howeer, bifurcation and transition are two different, but related concets. Definition A.1 defines the transition of (A.1) from a stable equilibrium oint to other states (not necessary equilibrium state). In general, we can define transitions from one attractor to another as follows. Definition A.2 Let Σ X be an inariant set of (A.1). We say that (A.1) has a transition of states from (Σ, ) on > (or < ) if the following conditions are satisfied: (1) when < (or > ), Σ is a local minimal attractor, and (2) when > (or < ), there exists a neighborhood U X of Σ indeendent of such that for any ϕ U \ (Γ Σ ), the solution u(t, ϕ) of (A.1) satisfies that lim su dist(u(t, ϕ), Σ ) δ() >, t lim δ() = δ, where Γ is the stable manifolds of Σ with codim Γ 1. Let the eigenalues (counting multilicity) of L be gien by {β j () C j = 1, 2, } Assume that < if <, Re β i () = if =, > if >, 1 i m, (A.4) Re β j ( ) < j m + 1. (A.5) he following theorem is a basic rincile of transitions from equilibrium states, which roides sufficient conditions and a basic classification for transitions of nonlinear dissiatie systems. his theorem is a direct consequence of the center manifold theorems and the stable manifold theorems; we omit the roof. heorem A.1 Let the conditions (A.4) and (A.5) hold true. hen, the system (A.1) must hae a transition from (u, ) = (, ), and there is a neighborhood U X of u = such that the transition is one of the following three tyes: (1) Continuous ransition: there exists an oen and dense set Ũ U such that for any ϕ Ũ, the solution u (t, ϕ) of (A.1) satisfies lim lim su u (t, ϕ) X =. t In articular, the attractor bifurcation of (A.1) at (, ) is a continuous transition. (2) Jum ransition: for any < < + ε with some ε >, there is an oen and dense set U U such that for any ϕ U, lim su u (t, ϕ) X δ >, t where δ > is indeendent of. his tye of transition is also called the discontinuous transition. (3) Mixed ransition: for any < < + ε with some ε >, U can be decomosed into two oen sets U 1 and U 2 (U i not necessarily connected): such that Ū = Ū 1 + Ū 2, U 1 U 2 =, lim lim su u(t, ϕ) X = ϕ U1, t lim su u(t, ϕ) X δ > ϕ U2. t he following theorem roides sufficient conditions for continuous transitions and gies local transition structure.

12 12 u = u = u = u = 1 2 FIG. 11: oological structure of the jum transition of (A.1) when k=odd and α > : < ;. Here the horizontal line reresents the center manifold. FIG. 12: oological structure of the continuous transition of (A.1) when k=odd and α < : ; >. We consider the transition of (A.1) from a simle critical eigenalue. Let the eigenalues β j () of L satisfy < if <, β 1 () = if =, (A.6) > if >, Reβ j ( ) < j 2, where β 1 () is a real eigenalue. Let e 1 () and e 1() be the eigenectors of L and L resectiely corresonding to β 1 () with L e 1 =, L e 1 =, < e 1, e 1 >= 1. Let Φ(x, ) be the center manifold function of (A.1) near =. We assume that < G(xe 1 + Φ(x, ), ), e 1 >= αx k + o( x k ), (A.7) where k 2 an integer and α a real number. We hae the following transition theorems. heorem A.2 Under the conditions (A.6) and (A.7), if k=odd and α in (A.7) then the following assertions hold true: (1) If α >, then (A.1) has a jum transition from (, ), and bifurcates on < to exactly two saddle oints 1 and 2 with the Morse index one, as shown in Figure 11. (2) If α <, then (A.1) has a continuous transition from (, ), which is an attractor bifurcation as shown in Figure 12. (3) he bifurcated singular oints 1 and 2 in the aboe cases can be exressed in the following form 1,2 = ± β 1 ()/α 1/k 1 e 1 () + o( β 1 1/k 1 ). heorem A.3 Under the conditions (A.6) and (A.7), if k=een and α, then we hae the following assertions: 1. (A.1) has a mixed transition from (, ). More recisely, there exists a neighborhood U X of u = such that U is searated into two disjoint oen sets U 1 and U 2 by the stable manifold Γ of u = satisfying the following roerties: U = U1 + U2 + Γ, the transition in U1 is jum, and (c) the transition in U2 is continuous. he local transition structure is as shown in Figure (A.1) bifurcates in U 2 to a unique singular oint on >, which is an attractor such that for any ϕ U 2, lim u(t, ϕ) t X =, where u(t, ϕ) is the solution of (A.1). 3. (A.1) bifurcates on < to a unique saddle oint with the Morse index one. 4. he bifurcated singular oint can be exressed as = (β 1 ()/α) 1/(k 1) e 1 + o( β 1 1/(k 1) ). We consider the equation (A.1) defined on the Hilbert saces X = H, X 1 = H 1. Let L = A + B. For L and G(, ) : H 1 H, we assume that A : H 1 H is symmetric, and < Au, u > H c u 2 H 1/2, (A.8) < Bu, u > H c u 2 H, (A.9) < Gu, u > H c 1 u H + c 2 u 2 H, (A.1) where > 2, c, c 1, c 2 > are constants.

13 13 heorem A.4 Assume the conditions (A.3), (A.4) and (A.8)-(A.1), then (A.1) has a transition at (u, ) = (, ), and the following assertions hold true: u = u = (1) If u = is an een-order nondegenerate singular oint of L +G at =, then (A.1) has a singular searation of singular oints at some (u 1, 1 ) H (, ). 1 2 u = (c) (2) If m = 1 and G satisfies (A.7) with α > if k=odd and α if k=een, then (A.1) has a saddlenode bifurcation at some singular oint (u 1, 1 ) with 1 <. FIG. 13: oological structure of the mixing transition of (A.1) when k=een and α : < ; = ; (c) >. Here U1 is the unstable domain, and U2 the stable domain. [1]. Ma and S. Wang, Bifurcation theory and alications, ol. 53 of World Scientific Series on Nonlinear Science. Series A: Monograhs and reatises, World Scientific Publishing Co. Pte. Ltd., Hackensack, NJ, 25. [2], Stability and Bifurcation of Nonlinear Eolution Equations, Science Press, 27. [3] L. E. Reichl, A modern course in statistical hysics, A Wiley-Interscience Publication, John Wiley & Sons Inc., New York, second ed., 1998.