Mass Generation via the Phase Transition of the Higgs Field

: The commonly quoted bistable Higgs potential is not a proper description of the Higgs ﬁeld because, among other technical reasons, one of its stable states acquires a negative expectation value in vacuum. We rely on formal catastrophe theory to derive the form of the Higgs potential that admits only one positive mean value in vacuum. No symmetry is broken during the ensuing phase transition that assigns mass to the Higgs ﬁeld; only gauge redundancy is “broken” by the appearance of phase in the massive state, but this redundancy is not a true symmetry of the massless ﬁeld. Furthermore, a secondary, certainly amusing conclusion, is that, in its high-energy state, the ﬁeld oscillates about its potential minimum between positive and negative masses, but it is doubtful that such evanescent states can survive below the critical temperature of 159.5 GeV, where the known particles were actually created.


Introduction
Bistable potential wells possessing two minima separated by an energy barrier are quite common in the natural sciences [1][2][3][4][5][6].Despite their frequent use in descriptions of discontinuous transitions occurring in physical, chemical, and biological systems and their intimate connections to catastrophe theory, bifurcation theory, singularity theory, structural stability, and phase transitions [4-9], the ensuing dynamical evolution is not understood in virtually all cases, to the point that some famous accounts of transitions are not only technically unphysical, but they are also visibly preposterous.The deeper reason for such absurdities is the lack of temporal variables in Landau's phase-transition theory and in Thom's catastrophe theory.These theories apply only to gradient systems [1,[4][5][6], and the notion of time-dependent phenomena is added ad hoc by describing arbitrarily drawn paths in the control parameter space of the cusp and higher elementary catastrophes.
For instance, Landau's phenomenological theory of second-order phase transitions predicts the appearance of two minima of equal depth past the critical point, although we know from experiments that only one stable state exists below the critical temperature T c .To work around this problem, the theory postulates, against the odds, that an evolving system will arbitrarily choose to settle into one of these states.Even in this hypothetical scenario, the model remains unphysical because these states continue to evolve and change their mean values as the temperature T < T c is lowered toward absolute zero.So, no matter which minimum the system "chooses", it finds itself out of equilibrium all the time; thus, the system has to evolve again and again trying to catch up with the ever-changing equilibrium state.In contemporary parlance, such a situation is described by the metaphor "moving the goalposts" which has a negative connotation.Furthermore, it is well-known that an infinitesimal linear perturbation wipes out entirely Landau's second-order phase transition [6], which means that such transitions should not occur in nature, or that the theory is wrong.All these absurdities come to life because of the assumption that the system finds itself at a local maximum of the potential as T crosses to just below T c , where it sees two new minima opening up as the control parameter becomes nonzero (negative, to be specific [1,[5][6][7]).This assumption places the system in an unstable initial state, a practice that is heavily at odds with basic physics and with the stable stationary states that we describe in Section 2 below.It is also mathematically puzzling how the initial potential minimum at T ≥ T c changes directly to an isolated maximum without passing through a degenerate inflection point (this procedure builds a priori a discontinuity in the second derivative of the potential [1]).
The stability of the Higgs potential in particle physics [9][10][11][12][13] is another case in point.The descriptions of how the Higgs field acquires mass are cursory and nonsensical at their roots: At high energies, the massless Higgs field is supposed to be stripped of any and all features, yet it is hypothesized to have "some high symmetry" (zero weak isospin?) 1  supporting an even potential function.This is the symmetry that will be broken in the final stationary state, when the unstable field will conveniently ignore the negative-minimum state and it will choose against the odds to settle into the other available state of positive vacuum expectation value (VEV).But how can such a symmetry be broken when the potential continues to be an even function, just as prior to the transition?And how can the system ever settle into either one of the low-energy states, when these states are not really stationary but continue to move the goalposts (dotted curves in Figure 1) to different VEVs all the time?
For m = 0, three unphysical features are observed: (a) The system suddenly finds itself at a local maximum.(b) Two global minima are available, the one at φ 0 < 0 with negative VEV. (c) The stable minima continue to relocate to φ 0 = ±m/2 (along the dotted curves) as m increases, throwing the system out of equilibrium all the time and preventing its settling to a specific VEV, irrespective of which side it chooses to evolve.These features appear because the perturbation (−m 2 φ 2 /2 + bφ) of the cusp catastrophe germ (φ 4 ) has been overconstrained by setting b ≡ 0.
All of the above descriptions should have been taken with a large grain of salt because, after all, an infinitesimal linear perturbation at T = T c eliminates the second-order phase transition altogether.This occurs because Landau's assumption of a "higher symmetry" in the initial state [1] arbitrarily alters the perturbation 2 (−m 2 φ 2 /2 + bφ) attached to Thom's cusp catastrophe germ (φ 4 ; [4]); as a result, one control parameter is eliminated (b = 0; [6]) and the drawn ad hoc evolutionary path {m, b = 0} in the control parameter space (m, b) becomes incorrect and irreparable-even if an infinitesimal b = 0 perturbation is brought back in.The reason for this structural instability is that m and b are related along the transition path; thus, the value of b cannot be chosen independently.The proof is given in Section 2.1 below using polynomial theory.
The resulting overconstrained (b = 0) potential with one remaining control parameter, V(φ) = φ 4 − m 2 φ 2 /2, is illustrated in Figure 1.The phase-transition path highlighted by the dotted curves is unphysical for the reasons discussed above; thus, naturally occurring phase transitions (of first and second order) require a different mathematical approach.We undertake this task for the Higgs field in Section 2, and we discuss our results for the various types of phase transitions in Section 3. Finally, we summarize our conclusions in Section 4. For the sake of completeness of the methodology, the two higher-order elementary catastrophes (the swallowtail and the butterfly) are also analyzed in this work, and their results are collected in Appendix A.

Derivation of the Higgs Potential from Catastrophe Theory
In cosmology and particle physics, the scalar Higgs field is massless and featureless at the very high energies occurring right after the big bang [9,10,[14][15][16][17].When the universe cools down to a critical temperature of T c = 159.5 ± 1.5 GeV [16,17], the electroweak phase transition takes place [9][10][11][16][17][18][19][20].Lattice monte-carlo simulations indicate that the cross-over of the Higgs field is smooth but fast, lasting for only ∆T ∼ 5 GeV [16,17], during which the field settles down to a nonzero (positive) VEV of v = 246.22GeV, where it has remained until the present time.This value of the Higgs VEV is a natural constant [20], and it is responsible for the corresponding particle, the Higgs boson, acquiring its observed mass (125.25 GeV; [21-24]).
There are two methods by which we can derive the scalar Higgs potential at all temperatures and observe the phase transition to the massive Higgs boson.The first derivation is more tedious and requires more steps, but it is also transparent in justifying the various assumptions being made; it further shows that the Higgs potential obeys Thom's theorem [4] for the cusp catastrophe.The second derivation is an astute shortcut, but it is opaque and reveals no details; this formulation hides the influence of catastrophe theory, so it could have been carried out at the time that Landau [1] presented his phase-transition theory.We summarize both methods below.

Method 1: Relying on Catastrophe Theory and Stable Isolated States
For the Higgs potential V(φ) to generally exhibit three isolated extrema, its derivative V must have the form where a, b, c are interrelated control parameters to be constrained below.Then, V (φ 0 ) = 0 gives the extrema φ 0 = −a, −b, −c.Integrating Equation (1), we find that where the integration constant has been dropped.In the neighborhood of the critical point of the germ V = φ 4 , the Taylor expansion does not have a cubic term or terms higher than O(φ 4 ).These terms are eliminated by Thom's inhomogeneous linear transformation and his nonlinear transformation, respectively [4][5][6].Thus, we must set in which case we obtain the canonical form of the cusp catastrophe 3 with the extrema located at φ 0 = −a, −b, (a + b).Note that if we arbitrarily choose b = 0, then we obtain Landau's [1] potential with extrema at φ 0 = 0, ±a (see also Section 2.2.2 below).This choice is unjustifiable, and we are not going to adopt it.Instead, we shift V(φ) by a to the right, in order to place the first listed extremum at φ 0 = 0.The shift transforms the cusp-catastrophe Function (4) to where an additive constant has been dropped (eliminated by a vertical shift).In this function, the extrema have been shifted to φ 0 = 0, (a − b), (2a + b).We shall see that φ 0 = 0 corresponds to a local minimum of V(φ; a, b), and we are prepared to assume that the massless Higgs field occupies this minimum while waiting for a more stable state to open up and become accessible.The subject of accessibility of a new global minimum is very important in this regard; it is discussed further in Section 3 below.Next, we fix the third listed extremum to always be located at φ 0 = 1 by convention.Then, we set b = 1 − 2a , (6) and Equation ( 5) takes the form The extrema are now located at φ 0 = 0, (3a − 1), 1.When φ 0 = 1 is a global minimum, it represents the massive state of the Higgs field, and when this minimum becomes accessible at a critical point in the control parameter plane (a, b), the field will make the transition to a nonzero VEV (v = φ 0 = 1).The phase-transition path is described by Equation ( 6).Thus, the path is an oblique line that does not cross the apex a = b = 0 of the separatrix.
In a final step, we redefine the location of the second listed extremum This definition gives us a better handle on the location of this extremum.We want it to correspond to a local maximum (the location of an energy barrier that obstructs the phase transition for all k ∈ (0, 1].As such, k should be located between the other two extrema, viz. and then Equation ( 7) is rewritten in the final form Looking at this potential function, it is hard to imagine that it satisfies Thom's cuspcatastrophe theorem [4], but it does.Equation ( 10) is equivalent to the cusp-catastrophe potential (4) shifted by a to the right, where a = (k + 1)/3 and b = (1 − 2k)/3.The Higgs potential ( 10) is plotted across the transition path {k = 1 → 0} in Figure 2. The critical points (φ 0 , V(φ 0 )) of the potential V(φ; k) are Thus, the height of the energy barrier is ∆V = k 3 (2 − k)/3 and 0 ≤ ∆V ≤ 1/3.For k ≤ 1/2, once at the top of the barrier, a system will dissipate an amount of energy equal to ∆E = (1 − k) 3 (1 + k)/3 during its settling to the global minimum on the right side.This amount is maximized at the critical point k = 0 for which ∆E max = 1/3.10)) with 0 ≤ k ≤ 1 placing the energy barrier between the other two extrema.The phase transition path is described by {k = 1 → 0}.
A system initially located at (0, 0) may undergo a phase transition to φ 0 = 1 only when this state becomes the global minimum (for 0 < k < 1/2) and (quantum tunneling aside) only if a finite perturbation provides the free energy required for climbing over the top of the intervening energy barrier.The critical point of the phase transition occurs for k = 0, when the diminishing barrier finally disappears, and the system moves spontaneously to φ 0 = 1.For k < 0 (dashed curve), a new local minimum opens up at φ < 0, but the system remains at the global minimum φ 0 = 1 for all k min < k < 0, where k min represents the state at absolute zero-here, as usual, we think of k as proportional to the temperature difference (T − T c ).
For k = 1/2, where the two minima have equal depth, the barrier height is ∆V = 1/16, and an equal amount of energy, if gained from external perturbations, will be dissipated away (∆E = 1/16) during the transition from the top of the barrier to one of the two stable states.The k = 1/2 stage is important because it is the first instance along the evolutionary path {k = 1 → 0}, where another stable state (φ 0 = 1) becomes available to a system located at φ 0 = 0, although, barring a sufficiently strong nonlinear perturbation, the new state is not dynamically accessible because of the intervening barrier [28][29][30][31][32].In Section 3, we discuss the types of viable phase transitions along the latter path segment {k = 1/2 → 0}, where the energy barrier continues to diminish with decreasing k.

Method 2: Implementing a Shortcut
An alternative derivation of Equation ( 10) that dispenses with details and formalities is as follows.
We return to Equation (1) for the derivative V (φ), which we copy here for convenience: Following Landau's assumption (Cφ 3 ≡ 0; [1]), we eliminate the quadratic term from V , in which case the sum of the three zeros is set to zero and c = −a − b: We shift φ by a to the right to place one extremum always at φ 0 = 0: We constrain the control parameters by 2a + b = 1 (or by a − b = 1) to place another extremum always at a fixed location φ 0 = 1: We redefine −3a + 1 by Equation ( 8) to simplify the location of the remaining extremum: Integrating with respect to φ, we obtain the form (10) for V(φ; k).

Utilizing a Familiarity Heuristic
Perhaps surprisingly, the steps taken in the shortcut above can all be avoided by utilizing a familiarity heuristic [33].
The final result can be written down in just two steps, without proof or investigation of its validity, by simply recalling that we are interested in static potentials which we can use to demonstrate phase transitions.Such potentials must generally exhibit three extrema, two fixed minima (φ 0 = 0, 1) representing the initial and final stationary states, and a maximum representing an obstacle or barrier that separates the two states.Therefore, Equation ( 16) can be written down ab initio, and then it can be integrated to yield the potential V(φ; k) shown in Equation (10).
The problem with this extremely fast, albeit heuristic approach is, of course, that we cannot then formally justify the potential obtained by intuition and familiarity with nature's phase transitions [1][2][3].This problem is solved by the lengthy derivation given in Section 2.1 above.

Looking Back to Landau's Theory of Phase Transitions
By contrast, Landau's phase-transition theory can be formulated in the same context (Equations ( 12)-( 16)) as follows.
Control parameter c is replaced by −(a + b) in Equation ( 12) to eliminate the φ 2 term (no cubic term in the potential).Then, b is set to zero in Equation ( 13) (in disagreement with catastrophe theory, which was not known at that time), resulting in the overconstrained form Integrating with respect to φ, we obtain the final form which is depicted in Figure 1 for a = 0, 0.5, 1 (corresponding to m = 0, 1, 2; here, a = m/2).
As was discussed in Section 1 and summarized in the caption of Figure 1, these potential curves do not form an evolutionary path in the control parameter plane (a, b) with varying values of the remaining control parameter {|a| = 0 → 0.5 → 1 → • • • }.

The Higgs Phase Transition
For the purposes of our discussion, we rewrite the canonical potential Function (4) of the Higgs field in Thom's equivalent form of a cusp catastrophe, viz.
where the control parameters (A, B) are functions of two of the roots (a, b) of The third root c is not independent, i.e., c = −(a + b) (Equation ( 3)).If we set any one of these three roots to zero, then B = 0 and the perturbation of the germ φ 4 takes a specialized even form that cannot describe quantitatively any phase transition since the phenomenon occurs naturally for general perturbations of no particular symmetry.Moreover, B = 0 fixes the maximum of the bistable potential V(φ) to point (0, 0), where the system finds itself at the onset of the phase transition.Thus, the system is unstable and has no choice but to evolve.This setup is clearly problematic. 4 Physical reasoning [2-4,7-9,15-17, [26][27][28][29][30][31][32] formally requires that the system be located at a stable minimum of the potential at all times before the second-order critical point T c is reached, and that this minimum become degenerate for T = T c (i.e., an inflection point) and progressively a (no longer relevant) local maximum for T < T c , as in Figure 2 (the case T < T c corresponds to the k = −0.25 curve).This figure also shows the physical reason for the occurrence of the second-order phase transition for k = 0: the energy barrier that separates the two stable states diminishes as T → T + c (k → 0 + ) and disappears altogether for T = T c ; in fact, it is the merging of this maximum with the minimum at φ 0 = 0 that makes the critical point T = T c degenerate.This smooth process makes sense, as the minimum that initially hosts the system switches gradually, first to an inflection point, and then to a maximum.
The control parameters A and B do not vary independently along the evolutionary path.Therefore, setting B = 0 for all values of A in Equation ( 19) (as in Landau's theory) is prohibited.This can be proven as follows: Using Equations ( 6) and ( 8), we express the control parameters (A, B) of the canonical cusp catastrophe (19) as functions of k, viz. and we can see that B = 0 in k ∈ [0, 1] only for a single point, k = 1/2, for which A = −1/2. 5 Now, eliminating k between these two equations, we find that the control parameters (A, B) are related along the path {k = 1 → 0} by where A ∈ [−2/3, −1/2] and B ∈ [−8/27, 8/27].This curve effectively constrains the evolutionary path in the (A, B) plane; the constraint reveals the presence of an integral of motion (i.e., a conserved quantity) during the evolution, as was determined in astrophysical first-order and second-order phase transitions [29][30][31][32].
Figure 3 shows the (A, B) plane of the cusp potential (19) and the evolutionary path {k = 1 → 0} that lies entirely within the separatrix (the fold curve 8A 3 + 27B 2 = 0) and terminates at the critical point k = 0, where the phase transition occurs spontaneously.(Note the degenerate inflection point at φ 0 = 0 in the inset of the Higgs potential V(φ; 0).)For k = 0, the coordinates are (A, B) = (−2/3, −8/27).Because this point also lies on the separatrix, this is the first demonstration of the so-called "delay convention" [3][4][5][6]37] in a second-order phase transition.Except that the delay down to k = 0 does not occur by convention here, it is a calculated outcome in the evolution of V(φ; k) depicted in Figure 2.This brings the discussion to the other convention 6 commonly used in catastrophe theory, the so-called "Maxwell convention" [3,6,28,37] used in first-order phase transitions.

The Maxwell Convention and Chemical Reactions
The Maxwell convention singles out the point with k = 1/2 in the middle of the path shown in Figure 3 as a viable phase-transition point because the two minima seen in the V(φ; 0.5) inset have the same depth [37,38].This is utter speculation that came about because a system was thought to already be at the top of the energy barrier.From the top, both minima are accessible with equal probabilities of transition, and the two stable states coexist.This setup and assumptions are basically the same as in the Higgs field which can also transition to the two stable states (one with negative VEV) with equal probability, according to Landau's theory [1].But, as we explained above, placing a system at a local maximum at the transition point is unphysical; so, we proceed to describe and clarify the evolution of nonspontaneous first-order phase transitions and Maxwell's rule under the action of external perturbations in the control parameter plane of Figure 3.  4) and (19).The lower insets are borrowed from the potential curves V(φ; k) of Figure 2. The separatrix 8A 3 + 27B 2 = 0 is shown in blue color.
Maxwell's rule [38] (the basis for the Maxwell convention) identifies the point (k = 1/2) along the path {k = 1 → 0} in Figure 3, in which the two stable minima attain equal depths.The system initially occupies the left minimum and, as the evolution proceeds along the segment {k = 1/2 → 0}, it cannot generally access the other stable state because of the intervening energy barrier.Thus, Maxwell's rule simply captures the first instance that another stable state becomes available, but not necessarily accessible.Only external perturbations can induce such a nonspontaneous transition of system parts and sectors for k ≤ 1/2, if they are sufficiently strong, and then the two phases will coexist.Thus, chances are that such a discontinuous transition of sectors may occur at a value smaller than k = 1/2 because the barrier height decreases along the segment k < 1/2 (Figure 2).This is precisely what takes place in the chemical reactions that use catalysts [39][40][41][42][43]; catalysts lower the so-called activation energy barrier, thereby increasing the reaction rates (i.e., they induce a first-order phase transition in parts of the reactants), without actually being consumed.Lowering the energy barrier is a mechanism used in catalyzed chemical reactions.An alternative mechanism is to perturb the reactants by supplying excess heat.In this pathway, the barrier remains intact, but the reactants absorb the energy, and more constituents go over the top of the barrier to the other state that hosts the products of the reaction.

Overcoming the Energy Barrier
The above chemical reaction mechanisms fit rigorously into our framework of firstand second-order phase transitions (Figures 2 and 3).A spontaneous reaction occurs when there is no barrier (k = 0); and a catalyzed reaction or a heat-driven barrier jump occurs for k ≤ 1/2, but only under the action of perturbations supplying the necessary energy.An example is shown in Figure 4 for k = 0.4 in Equation (10).We consider a system oscillating initially about minimum I under the action of external perturbations.Since k < 1/2, the second minimum S that became available for k = 1/2 is now the global minimum and the energy barrier ∆V has decreased in height past the Maxwell point.
We distinguish three cases in the free-energy diagram sketched in Figure 4: (a) If the Gibbs free energy ∆E gained by the parts of the system is not sufficient to push any part up to at least point P (or B), then the perturbed system remains in the neighborhood of point I.(b) If, on the other hand, ∆E = ∆V in some parts, then these parts displaced to point P can overcome barrier B and roll over to the new global stable state S [2,3,[28][29][30][31][32].Then, the two phases, I and S, coexist [34][35][36].
(c) Furthermore, if ∆E > ∆V in some perturbed parts (displaced, e.g., up to point J), then these parts no longer recognize barrier B and collapse to the deep minimum S on a dynamical time [3,44,45].
In most physical systems undergoing phase transitions, the evolutionary paths just outlined cannot be obtained analytically as functions of time because the partial differential equations of motion cannot be solved analytically [46,47].For this reason, researchers in the past have approached the subject either by studying only the stationary minima of the potentials (e.g., [28][29][30][31]), thereafter assuming that quasistatic evolution takes place between the stable states; or by relating the stationary states to analytic non-dissipative solitary-wave solutions conserving localized finite energy (the topological method of "orbits"; [47][48][49]); or by numerical simulations (e.g., [2,16,17,46]).

Star-Forming Phase Transitions
Although the outcome of the above evolutionary scenarios is the same in Figure 4 (the settling of at least parts of a system into stable state S), the dynamics is quite different.The difference was recognized long ago in the context of star formation in giant molecular clouds, first by Whitworth [28] and subsequently by Tohline [2,3,[29][30][31], although the same ideas had been previously explored in various related contexts [50][51][52].Whitworth [28] described a perturbed diffuse molecular cloud region bound by external pressure, that reaches over time [3] (point P in Figure 4, as "preunstable," a condition that differentiates it from a region strongly compressed and displaced to point J, where it becomes Jeans unstable [44] and subject to dynamical collapse down to the compact (high-density) stellar state (point S in Figure 4).Tohline [30] recognized that the path (PIBS in Figure 4) highlights a slower phase transition (distinct from dynamical Jeans collapse) capable of producing stars of much lower masses (albeit over much longer timescales [3]), as compared to the famous Jeans critical mass [44], the hallmark of dynamical star formation since 1902 and for years to come [53].10) with k = 0.4) in which a system is displaced from equilibrium and oscillates about the local minimum I under the action of external perturbations.
If the system gets displaced to P, it gains enough energy (∆E = ∆V) to roll over the top of barrier B and down to the global minimum state S, as indicated by the arrows.Thus, the system at P is "preunstable" [28] and undergoes a first-order phase transition [2,3,[29][30][31]).If the system gets displaced to any point of higher energy (e.g., at J), then it becomes dynamically (Jeans) unstable [44] (it no longer recognizes the energy barrier at B), and collapses to the global minimum state S. If the system never gains enough energy to overcome the barrier, then it will remain near point I until k = 0, where the barrier disappears.Then, points I and B merge to an inflection point and the second-order phase transition to S is spontaneous [5,29,45].

Peculiar λ-Transitions
To complete the discussion of the various types of phase transitions encountered in nature, we should mention that some phase transitions do not fit into the modern classification scheme [35,36].Most puzzling among them are the so-called λ-transitions [8,32,[34][35][36]54,55] that may or may not [8,35] have infinite specific heat at the critical point (e.g., at the λ-point T c = 2.18 K of the superfluid liquid 4 He [34][35][36] or at the order-disorder critical point T c = 739 K of the β-brass Cu-Zn alloy [34,56]).Bose-Einstein condensation of an ideal Bose gas [34] and astrophysical binary fission and ring formation [32] are also types of λ-transitions, and the various types are all linked together only by spontaneous breaking of the topology [8,32] (the symmetry may break or not, and the specific heat may diverge or not [8,35,[54][55][56]).
From the viewpoint of the energetics of discontinuous λ-transitions, we know that a total of five extrema (not all of them isolated) are involved in the Gibbs free-energy function [32], which places these transitions along paths in the higher-order butterfly catastrophe [55], if the free energy is a continuous function of the order parameter [32].The main characteristic of the underlying potential function is an energy barrier that progressively becomes taller as T → T + c , and then, it suddenly disappears just past the critical point T = T c (see, e.g., Figures 3-5 in Ref. [32]).This astonishing behavior of the free-energy barrier in astrophysical systems exhibiting topology-breaking phase transitions [57][58][59][60][61][62][63] remains under investigation to this day (see Refs. [2,3,8,[29][30][31][32] for more details).

Conclusions
(a) Figure 1 shows the potential functions for Landau's phenomenological theory of second-order phase transitions [1], including that ascribed to the Higgs field.As the control parameter is increased (to simulate time evolution), the potential develops two features that render it unphysical: two symmetric stable global minima appear on either side of the local maximum at φ = 0, and they both continue to move away from φ = 0 in time.Thus, assuming an initial state at φ = 0, as usual [4-6], the phase transition does not produce a unique or a universal final state.But these properties are required for the Higgs field at present [9][10][11][12][13]20]; and uniqueness of the final state is required for many phase transitions in solids and fluids [8,[29][30][31][32][34][35][36]45,54].(b) Figure 2 shows how to get around the problems highlighted above.The potential functions of the cusp catastrophe [4] are all shifted so that one minimum is always at φ = 0 (representing the initial state) and another minimum is constrained to always be at φ = 1, making this (final) state universal [20].The linear term of the cusp catastrophe (Equation (4)) precludes the appearance of another minimum at φ < 0 before the second-order critical point at k = 0 is reached; and when such a minimum finally appears for k < 0, it is incapable of influencing the dynamics of the phase transition that took place spontaneously already for k = 0.
The maximum that appears for 0 < k < 1 represents a free-energy barrier between the two isolated stable states.External perturbations may drive a system from φ = 0 to φ = 1, if they supply the requisite energy to overcome the intervening barrier (a first-order phase transition [2,3,30,31]), otherwise the system remains oscillating about φ = 0.As the control parameter decreases toward k = 0 (Figure 2), the barrier becomes shorter (just as in catalyzed chemical reactions [39][40][41][42][43]), and a second-order phase transition appears at k = 0, where the barrier disappears [29,32,45].We believe that such a transition occurred in the massless Higgs field when it acquired its uniquely positive universal VEV [16,17,20] because we cannot imagine vacuum perturbations strong enough to overcome the energy barrier in the interval 0 < k ≤ 1/2.Before the phase transition occurs, the Higgs field is induced to executing small amplitude oscillations about the minimum at φ = 0 that represents the equilibrium VEV of the massless state.Such oscillations generate evanescent particles with both positive and negative masses that do not survive long into the future.The observed particles of our times were all assigned masses after the Higgs field had settled to its universal VEV of 246.22 GeV [20,64].(c) Figure 3 shows the evolutionary path in the control plane of the cusp catastrophe.The path remains within the fold lines of the separatrix at all times (even for 0 > k > 1), and exhibits a Maxwell critical point [37,38] for k = 1/2 and a second-order critical point for k = 0.As the green curve in Figure 2 shows, a first-order stable minimum at φ = 1 becomes available for k = 1/2, but it is not necessarily accessible to a system located at φ = 0 via a first-order phase transition due to the intervening energy barrier.For the Higgs field, such a continuous line of first-order phase transitions was observed in the past (before the discovery of the Higgs boson and the measurement of its mass; see Refs.[16,17,46] and references therein); although the lattice simulations were using Landau's even-symmetric potential and a Higgs mass of no more than 73.3 ± 6.4 GeV.For higher Higgs masses in non-perturbative simulations, the second-order critical point was replaced by a smooth crossover to the final massive state [16,17,46].These doubtful results must have their origin in the unphysical potential used, and new simulations are needed to revisit the true nature of the Higgs phase transition.(d) Figure 4 shows a schematic illustration of the various aspects of first-order phase transitions capable of overcoming the intervening energy barrier [3].Basically, there are two separate evolutionary modes depending on the amount of energy deposited by acting external perturbations over time: (i) a strongly perturbed system (point J in Figure 4) is not impeded by the barrier any longer and makes the transition to the final equilibrium state S on a dynamical time [30,31,44]; and (ii) a system oscillating about point I, and perturbed gradually upward to point P (or B), gains enough energy to jump over the top of the barrier B and down to the final equilibrium state S [28][29][30].
where k ≡ a − b and ≡ 3a + b − 1. Integrating Equation (A2) with respect to x, we find the potential function We see now that the choice of x 0 = 1 has limited the control space to only two independent parameters (k, ).This choice, which has been overlooked for generations, is necessary to create and define another stable state, so that we can apply this potential to actual physical systems.(The initial stable state created by the shift is also fixed at x 0 = 0.) We must say at this point that any arbitrary paths drawn in the deceiving general three-dimensional swallowtail control space are meaningless, in the sense that physical systems do not evolve unconstrained along such paths that keep moving the goalposts (see Section 1).
We also choose the extremum x 0 = k to be between 0 and 1, that is, to serve as an energy barrier between the two stable states.Thus, 0 ≤ k ≤ 1, allowing for N = 3 possible k locations in the interval k ∈ [0, 1] (0, 1, and in-between).Now, can be located anywhere on the x-axis, so there are N = 17 possible locations for the pair (k, ).Of those, the extrema x 0 = 0, 1 are local minima in only one case, in which < 0 (N = 3 cases, if we also count the degeneracies k = 0, 1).Therefore, only the case < 0 is of interest to phase transitions along the path {k = 1 → 0}.Now, the isolated extremum x 0 = < 0 is always a local maximum, and it can vary just as k varies within its own interval.But variation in does not change the qualitative properties of the transition, so we can assume here for demonstration purposes that is a negative constant along the considered transition path.In a physical system, however, the variation in will have to be determined from the physical parameters of the system itself.
Phase transitions.-Anillustration with constant = −1 (fixed) is shown in Figure A1.The transition proceeds on the right half of this diagram just as it does for the cusp potential in the main text.The second-order critical point appears for k 2 = 0 (the inflection point at x 0 = 0 on the magenta curve).The first-order critical point appears for k 1 = 8/15 (two equal-depth minima on the green curve), as determined from the equation for = −1 (note that k 1 → 1/2 as → ±∞, and the cusp catastrophe is fully recovered).Left energy barrier.-Point= −1 marks the location of another energy barrier on the left side of the diagram, and this barrier may be important in some applications concerned with systems falling apart: Assuming that a system (initially executing small-amplitude oscillations about x 0 = 0) can somehow be induced to climb over the top of this barrier (before it settles to the stable minimum x 0 = 1), then this system is doomed; it will certainly be destroyed since there is no other minimum of the potential available in the region x < .This path is however of no interest in customary applications of the swallowtail catastrophe, in which researchers are studying phase transitions terminating at stable states [4-7,65-68], such as x 0 = 1 in Figure A1.triple point along this phase-transition path.Then, it is easy to recognize that the control parameters (k, , m) must be related for a triple point to appear in the potential (A7).Their relationships are expressed by the conditions that V( ) = V(1) ≡ 0 . (A9) Here, and m are interchangeable parameters, so we chose the third minimum to be located at x 0 = .The V ≡ 0 equal-depth conditions (A9) then require that = −1 and m = −k = −1/ √ 3 , (A10) where now k and m are interchangeable (but we break the symmetry by choosing 0 ≤ k ≤ 1, as usual).It is easy to prove then that V(x) at the triple point is an even function of x, and this is why the third minimum (x 0 = or m) must be located at x 0 = −1.In fact, the potential V(x) at the triple point takes the simple form This reduced butterfly potential that exhibits a triple point (i.e., three minima of equal depth) is illustrated in Figure A4, where we chose k > 0 and = −1.The choice m = −1 is of course an alternative, and then the labels and m switch places in Figure A4.With the equal-depth minima set at x 0 = 0, ±1, then the extrema x 0 = k, m represent energy barriers of equal height.Finally, for evolutionary {k = 1 → 0} that exhibits a triple point, the following general relations hold along the path: = 3km and m = −k ; (A12) so, only one of the control parameters (k, , m) of the butterfly turns out to be independent in this model.The triple point occurs for = −1 along this path which, in terms of , is described by { = −3 → 0} and terminates at the second-order critical point = 0. Figure A5 provides an illustration of this phase transition.The control parameters of the potential Function (A7) have been reduced to functions of k by using the relations (A12), viz.
= −3k 2 and m = −k, and then k is the only independent parameter along the evolutionary path.

Figure 2 .
Figure 2. Cusp catastrophe in the potential V(φ; k) (Equation (10)) with 0 ≤ k ≤ 1 placing the energy barrier between the other two extrema.The phase transition path is described by {k = 1 → 0}.A system initially located at (0, 0) may undergo a phase transition to φ 0 = 1 only when this state becomes the global minimum (for 0 < k < 1/2) and (quantum tunneling aside) only if a finite perturbation provides the free energy required for climbing over the top of the intervening energy barrier.The critical point of the phase transition occurs for k = 0, when the diminishing barrier finally disappears, and the system moves spontaneously to φ 0 = 1.For k < 0 (dashed curve), a new local minimum opens up at φ < 0, but the system remains at the global minimum φ 0 = 1 for all k min < k < 0, where k min represents the state at absolute zero-here, as usual, we think of k as proportional to the temperature difference (T − T c ).

Figure 3 .
Figure 3.The evolutionary path {k = 1 → 0} on the control parameter plane (A, B) of the cusp catastrophe, where A and B are the coefficients seen in Equations (4) and(19).The lower insets are borrowed from the potential curves V(φ; k) of Figure2.The separatrix 8A 3 + 27B 2 = 0 is shown in blue color.

Figure 4 .
Figure 4. Bistable potential (Equation(10) with k = 0.4) in which a system is displaced from equilibrium and oscillates about the local minimum I under the action of external perturbations.If the system gets displaced to P, it gains enough energy (∆E = ∆V) to roll over the top of barrier B and down to the global minimum state S, as indicated by the arrows.Thus, the system at P is "preunstable"[28] and undergoes a first-order phase transition[2,3,[29][30][31]).If the system gets displaced to any point of higher energy (e.g., at J), then it becomes dynamically (Jeans) unstable[44] (it no longer recognizes the energy barrier at B), and collapses to the global minimum state S. If the system never gains enough energy to overcome the barrier, then it will remain near point I until k = 0, where the barrier disappears.Then, points I and B merge to an inflection point and the second-order phase transition to S is spontaneous[5,29,45].

Figure A1 .
Figure A1.Swallowtail potential functions for the phase-transition path {k = 1 → 0}.Parameter k controls the location of the right barrier which disappears for k 2 = 0. Parameter controls the left barrier, which we have fixed at x 0 = −1 for simplicity.The left barrier becomes shorter as → 0 − .A system that somehow is induced to overcome this barrier before it manages to settle to the x 0 = 1 minimum will fall apart.

Figure A2 .Figure A3 .
Figure A2.Butterfly potential functions for the phase-transition path {k = 1 → 0} with = m = −1.Parameter k controls the location of the barrier, which disappears for k 2 = 0. Parameters , m < 0 control extrema that develop in the x < 0 region.In the = m case shown here, the extrema degenerate to an inflection point at x 0 = −1.The height of this inflection point decreases as = m → 0 − .

3 Figure A4 .
Figure A4.Triple point of a butterfly potential with = m.The potential is now an even function (Equation (A11)) with no independent control parameter (Equation (A10)).With x 0 = k chosen to lie in (0, 1) to provide a barrier, then interchangeable parameters , m < 0 provide the locations of the two isolated extrema on the left side (i.e., x 0 = and x 0 = m).