Replica Field Theory for a Generalized Franz–Parisi Potential of Inhomogeneous Glassy Systems: New Closure and the Associated Self-Consistent Equation

On approaching the dynamical transition temperature, supercooled liquids show heterogeneity over space and time. Static replica theory investigates the dynamical crossover in terms of the free energy landscape (FEL). Two kinds of static approaches have provided a self-consistent equation for determining this crossover, similar to the mode coupling theory for glassy dynamics. One uses the Morita–Hiroike formalism of the liquid state theory, whereas the other relies on the density functional theory (DFT). Each of the two approaches has advantages in terms of perturbative field theory. Here, we develop a replica field theory that has the benefits from both formulations. We introduce the generalized Franz–Parisi potential to formulate a correlation functional. Considering fluctuations around an inhomogeneous density determined by the Ramakrishnan–Yussouf DFT, we find a new closure as the stability condition of the correlation functional. The closure leads to the self-consistent equation involving the triplet direct correlation function. The present field theory further helps us study the FEL beyond the mean-field approximation.


I. INTRODUCTION
Glass, an amorphous solid with elasticity, has a microscopic structure in which localized particles oscillate around their mean positions of a random lattice [1][2][3][4][5][6].The spatial randomness is self-generated by the particle localization that breaks translational symmetry.A remarkable feature of the random structure is that glass microscopically lacks the long-range order and is similar to liquid in terms of density-density correlations [1][2][3][4][5][6].As a precursor to the random structure of glass, supercooled liquids show heterogeneity over space and time [1][2][3][4][5][6][7][8][9][10][11].The dynamical heterogeneity and facilitation [4,[7][8][9][10][11][12][13][14] emerges on approaching the dynamical transition temperature (T d ), accompanied by the crossover from relaxational to activated dynamics.For T > T d , transport is not collective on a large scale.For T → T d , the system gets stuck in a glassy metastable state, and the dynamical behaviors of supercooled liquids exhibit features such as a two-step decay with a first relaxation (β-relaxation) to a plateau followed by a stretched exponential relaxation (α-relaxation) of density fluctuations.Along with this dynamical crossover, the dynamics become progressively heterogeneous and correlated in space.For T < T d , the relaxation times of the two-step decay increase rapidly despite slight changes in the disordered microstructure.
Various theories have tried to explain the dynamical heterogeneity and facilitation, as well as the dynamical crossover at T d .These include either elasticity theory or kinetically constrained models focusing on the dynamical facilitation [4,[12][13][14], and the mode coupling theory (MCT) [15,16], a dynamical theory relevant when approaching T d from the liquid phase.The MCT describes the onset of the two-step relaxation above T d and predicts the divergence of βrelaxation time at T d .Extension to the inhomogeneous MCT further allows us to describe a growing dynamical heterogeneity using a time-dependent three-or four-point correlation function [2,[7][8][9][10][11]].Yet the dynamical transition temperature T d is higher than the glass transition temperature observed in simulation and experimental studies.An interpretation of this discrepancy is that a mean-field description of the MCT is beyond the scope of the barrierdominated dynamics between metastable states, though applicable to the relaxation dynamics within a metastable state.The divergent behavior of the two-step decay predicted by the MCT at T d becomes incomplete because of the activated events remaining in actual liquids for T ≤ T d [15,16].
The activation dynamics dominant below T d are due to transitions between metastable states [15][16][17][18][19][20].Therefore, the dynamical crossover implies the emergence of many metastable states at T d or the appearance of a free energy landscape (FEL) characterized by an exponentially large number of metastable states below T d .From the thermodynamic point of view, we can describe the characteristic of the FEL using the configurational entropy obtained from the logarithm of the number of metastable states [1,2,[17][18][19][20].The Adam-Gibbs relation provides results in quantitative agreement with simulation and experimental results, relating the drastic changes in the relaxation time and viscosity to the decrease of the configurational entropy on approaching the glass transition temperature [17][18][19][20][21][22].For example, simulation studies on mixtures interacting via the Lennard-Jones potential and its repulsive counterpart, the WCA one, demonstrate that these systems exhibit quite different dynamics despite having nearly identical structures [23][24][25][26][27][28].Such a large difference in the dynamics is ascribable to a considerable gap between the configurational entropies while making a slight difference between the two-point correlation functions.Previous investigations confirmed that the configurational entropies associated with correlation functions differ greatly between the Lennard-Jones and WCA mixtures despite the structural similarity, therefore predicting the distinct dynamical behaviors from the Adam-Gibbs relation [23][24][25][26][27][28].
The Franz-Parisi (FP) potential obtained in the RHNC approximation serves as a starting point for this paper.The FP potential [65][66][67][68][69][70][71][72][73][74][75] is a function of overlap Q, a weighted average over the system of the two-point correlation function, and plays the same role as the Landau free energy of a global parameter Q that indicates a distance between the two copies in configuration space.Theoretical and simulation studies have demonstrated that the FP potential reproduces the temperature evolution of FELs, just like the Landau free energy [65][66][67][68][69][70][71][72][73][74][75].With decreasing temperature, the FP potential develops a secondary minimum for Q > 0 representing a metastable state.Considering Q = 0 in the liquid state, we can see that the potential difference, V (Q) − V (0), corresponds to the entropic cost of localizing the system in a single metastable state (i.e., the configurational entropy).
In this paper, we generalize the FP potential by fixing an inter-replica correlation function instead of the overlap Q.We formulate the generalized FP potential by developing a new framework that is beneficial to investigate the FEL while considering inhomogeneous supercooled liquids with the help of field theoretical method.A field theory combining the DFT [53][54][55][56][57][58][59][60][61] and replica theory [17-20, 29-38, 53-57] forms the basis of our framework.There are two requirements to be satisfied by the field theory and the associated functional.The first requirement is that the developed framework can consider inhomogeneous systems.The second requirement is that the generalized FP functional applies to non-equilibrium states away from metastable states.To meet the requirements, this paper presents the correlation functional theory that provides the generalized FP potential functional without going through the Morita-Hiroike functional [33][34][35][36][37][38][62][63][64].The generalized FP potential has three features as a functional of density and correlation function.First, this potential is a functional of metastable density that becomes equal to that of the DFT in the limit of m → 1.Second, the field-theoretical perturbation method allows us to have a new correlation functional different from the Morita-Hiroike one while maintaining consistency with the liquid theory in that the approximate form reduces to the RHNC functional.Last, the potential functional of a given inter-replica correlation function has a minimum where a new closure reducible to the RHNC approximation [33-38, 62, 63] holds.A remarkable result is that an approximation of the new closure yields the self-consistent equation for a non-ergodicity parameter that includes the triplet direct correlation function (DCF) [62,63,[76][77][78], similar to that formulated by either the MCT [15,16] or the replica theory [37,56,57], respectively.
The paper is organized as follows.In Section II, we define the generalized FP potential.Comparison between the generalized and original FP potentials clarifies what we modify through the generalization.Section III summarizes the theoretical results consisting of four parts as follows: relation for obtaining the generalized FP potential from the grand potential of m-replica system with inter-replica correlation function fixed (Result 1); functional form of the constrained grand potential (Result 2); new closure for two-point correlation function (Result 3); the associated self-consistent equation for a non-ergodicity parameter (Result 4).We obtain the generalized FP potential from Result 2 with the help of the relation in Result 1.The extremum condition of this potential yields a new closure in Result 3. It also turns out that a self-consistent equation obtained from an approximate form of the closure involves the triplet DCF as presented in Result 4. In Section IV, we calculate the perturbative terms using a strong-coupling perturbation theory developed for obtaining Result 2 (see Appendix B).In the saddle-point approximation, the strong-coupling perturbation theory provides the correlation functional form of the constrained grand potential given in Result 2. In Section V, we make some concluding remarks.

II. GENERALIZED FRANZ-PARISI (FP) POTENTIAL
We generalize the FP potential in comparison with its original definition.
The FP potential V (Q) is obtained in two steps [65][66][67][68][69][70][71][72][73][74][75].First, we fix a reference configuration ρ 1 of replica 1, which plays the role of quenched variable in the effective potential V (Q, ρ 1 ) as seen from the following definition: where Ca denotes (1/N !) • • • dr a,1 • • • dr a,N in the canonical ensemble of replica a for a ≥ 2, β the inverse thermal energy (k B T ) −1 , and U a ( ρ 1 , ρ a ) the interaction energy of replica a that is the sum of intra-replica interaction energy U intra ( ρ a ) and inter-replica one U inter ( ρ 1 , ρ a ): where using the intra-replica interaction potential v(r) and the inter-replica one v(r).It is noted that the effective potential V (Q, ρ 1 ) is defined in the absence of inter-replica interactions as represented by Equation (3).Next, we perform the canonical average of V (Q, ρ 1 ) over all possible choices for the reference configuration with the statistical weight p eq ( ρ 1 ) as follows: p eq ( ρ 1 ) = e −βUintra( ρ1) C1 e −βUintra( ρ1) . ( The replica trick allows us to calculate Equation (7), thus obtaining the FP potential V (Q) of the Landau type.

B. Generalization
Here, we introduce a generalized FP potential W ( G) as a functional of prescribed correlation function G(r, r ′ ), instead of the overlap Q.In terms of the Landau theory, we consider a local order parameter, instead of the global one.We use the grand canonical ensemble represented by the following operator: where the chemical potential βµ in units of k B T determines the most probable number N * , thereby providing the uniform density ρ = N * /V common to each replica with volume V .Given a reference configuration C 1 of replica 1, we have the interaction energy U a ( ρ 1 , ρ a ) of replica a, providing the grand potential ω a ( ρ 1 ) of replica a as follows: where the functional integral representation in Equation ( 10) is obtained from multiplying the right-hand side (rhs) of the first line in Equation ( 10) by the following identity: Equation ( 12) implies that The relation ( 13) at b = 1 represents that only the density distribution ρ Equations ( 10)-( 14) reveal that the field-theoretical formulation of the effective potential W ( G, ρ 1 ) can be developed as follows [79][80][81][82][83][84][85][86][87][88]: where In the last line of Equation ( 16), we have the ideal gas entropy defined by The generalized FP potential W ( G) is obtained from the grand canonical average of W ( G, ρ 1 ) for the reference configuration as follows: similar to Equations ( 7) and (8).Equation ( 18) clarifies that a given configuration ρ 1 plays a role of quenched disorder to another replica a [65][66][67][68][69][70][71][72][73][74][75].Since we consider all possible configurations of ρ 1 , the statistical weight P eq ( ρ 1 ) is of the Boltzmann form as well as p eq ( ρ 1 ) in Equation ( 8).Several remarks on Equations ( 13)-( 19) are in order: • Equation ( 14) tells us that a prescribed correlation field G(r, r ′ ) represents a product ρ 1 (r)ρ a (r ′ ) of two instantaneous density distributions in different replicas, or a statistical realization of density-density correlation function [62,63].
• To perform the configurational integral Tr a in the second line on the rhs of Equation ( 16), it is indispensable to introduce the Fourier transform representation of the delta functional using the functional integral over the one-body potential field, which is dual to the density field ρ a (r) [84][85][86][87][88].The ideal gas entropy given by Equation ( 17) appears in the last line of Equation ( 16) due to the saddle-point approximation of the one-body potential field [84][85][86][87][88].
• When different replica particles form complexes because of the attractive inter-replica interactions between them (i.e., v(r) < 0), we have G(r)/ρ 2 ≫ 1 in an overlapped region (e.g., |r| ≤ a), thereby providing a significant value of overlap Q that is greater than the random overlap obtained from G(r) = ρ 2 .The glassy state preserves an overlapped state due to frozen configurations of particles even after the attractive inter-replica interactions are switched off (i.e., v(r) → 0).The generalized FP potential W ( G) is available to explore such an overlapped state that is locally stable.

III. MAIN RESULTS
We present four sets of main results based on the strong-coupling perturbation theory (see Appendix B for details).Figure 1 summarizes the results schematically.
It is readily seen from Equations ( 16) and ( 20) that the constrained grand potential Ω The replica trick allows us to have the relation between the constrained grand potential Ω (m) v and the generalized FP potential W ( G): which is the first result (Result 1; see Appendix A for the detailed derivation).It is noted that the conventional replica trick proves the necessity of m → 1 to consider the quenched type of the FP formalism, though it has been physically motivated to take the limit of m → 1 based on the Monasson formalism [18,22,75].In Result 2, we provide the correlation functional form of the constrained grand potential Ω (m) v ( G). Section IV will sketch how the perturbative field theory developed in Appendix B yields the correlation functional given in Result 2.
Let us consider the inhomogeneous system characterized by the mean-field density ρ * a (r) satisfying where c ab (r) denotes the two-point DCF (simply called DCF) between replica a and replica b.Here we suppose that a given function G(r, r ′ ) imposed on the inter-replica correlation between replica 1 and replica a (a ≥ 2) is expressed as using a statistical realization of inter-replica radial distribution function g(r) or inter-replica total correlation function (TCF) h(r) ≡ g(r) − 1 [62,63].Namely, Equations ( 14) and ( 25) imply the constraint, on ρ 1 (r)ρ a (r ′ ) which is a statistical realization of density-density correlation [62,63] as mentioned above.Equation (26) includes the trivial inter-replica constraints as follows: one constraint, ρ 2 = ρ 1 (r)ρ a (r ′ ) (i.e., h(r − r ′ ) = 0), forces the two-replica system to maintain uniformity without inter-replica correlations, whereas another constraint, 0 = ρ 1 (r)ρ a (r ′ ) (i.e., h(r − r ′ ) = −1), imposes a region where two particles of different replicas exclude each other.
In Section III C, we will see that the metastable TCF h * (r − r ′ ) corresponds to the TCF obtained from averaging over statistical realizations of instantaneous density-density correlation ρ 1 (r)ρ a (r ′ ) consistently with Equation ( 26) as well as the liquid-state theory [62,63].Let h(r) and c(r) be the correlation matrices of TCFs and DCFs, respectively.The intra-and inter-replica matrix elements vary, depending on whether replica 1 is included or not: when setting χ(r) = h(r) or c(r) with the subscripts of their matrix elements denoting a pair of replicas, χ 11 (r) = χ 1 (r) and χ aa (r) = χ(r) for a ≥ 2, whereas χ 1a (r) = χ a1 (r) = χ(r) for a ≥ 2 and χ ab (r) = χ ′ (r) for a = b and a, b ≥ 2. As a consequence, we see from Equation ( 24) that where ρ * a (r) = ρ * (r) for a ≥ 2. It is noted that the metastable density distribution ρ * 1 (r) reduces to that from the Ramakrishnan-Yussouf density functional [58][59][60][61]: in the limit of m → 1.
The second result (Result 2) can be obtained using the perturbative field theory at strong coupling (see Appendix B).It will be shown in Section IV that the constrained grand potential is of the following functional form: where the matrix elements of v ′ has a non-zero potential v 1a (r) = v a1 (r) = v(r) between replica 1 and replica a that enforces Equation ( 26) without the constraint ∆ a ( G, ρ), and the matrix element of S is given by S ab (r) = δ ab δ(r) + ρ * a (r 0 )h ab (r).It is noted that the last line of Equation ( 32) is reduced to the RHNC functional of h(r) in the approximation of e h(r) − g(r) ≈ h 2 (r)/2 [33-38, 62-64, 85, 92].The stationary condition of W ( G) given by Equation ( 23) can be written as It is found from Equation (32) that where the third line of Equation ( 34) is obtained from the derivative of the logarithmic term in the fifth line of Equation (32) with respect to m using the Laplace expansion of |S| along the first row as follows: where use has been made of the cofactor expansion in calculating ∂ |S(r)| /∂m.It follows from Equation ( 34) that the stationary condition (33) becomes where the subscript 1 has been dropped because of the indistinguishability of all replicas in the limits of m → 1 and v(r) → 0, δ h/δρ * and its inverse are ignored, and the first term on the rhs is an approximate form obtained from the third line of Equation (34) (see Appendix F for the detailed derivation).We can easily verify the equivalence between Equation ( 36) and the following closure: which corresponds to the third result (Result 3), a new closure in the context of the liquid-state theory [62,63].Two remarks on Equations ( 33), (37) and (38) are in order: • Equation ( 33) is valid when a metastable state at h * (r) = h(r) is stable in the vanishing limit of the interreplica interaction potential (i.e., v(r) → 0); otherwise, transitions between basins occur in the FEL and the inter-replica correlations disappear, thereby amounting to g * (r) = 1 + h * (r) = 1, the trivial solution to Equation (33).In other words, the new closure (37) applies to the metastable state defined by Equation (33).
• The bridge function B(r) given by Equation ( 38) is approximated by B(r) = − h 2 * (r)/2, which coincides with the main term of either the soft mean spherical approximation (MSA) or various approximations used for hard-sphere systems [63,85].

D. Result 4: Self-Consistent Equation for the Non-Ergodicity Parameter
In the fourth result (Result 4), we restrict ourselves to uniform systems in Fourier space.We introduce the non-ergodicity parameter f (k) by relating the inter-replica TCF h * (k) to the intra-replica structure factor S(k) = 1 + ρh * (k) [15,16,[33][34][35][36][37][38]: We need to find an approximation of the closure (37) that is available to obtain the self-consistent equation including terms up to quadratic order in the non-ergodicity parameter f (k).It is appropriate for this purpose to expand the rhs of the closure (37), providing Equation ( 40) reads in Fourier space when making the approximation of B(r) ≈ −h 2 * (r)/2 as remarked after Equation (38).Meanwhile, the neglect of inhomogeneity (i.e., ρ * (r) = ρ) allows us to express the Fourier transform of the Ornstein-Zernike Equation ( 31) at m = 1 as using the non-ergodicity parameter f (k) defined by Equation (39).

IV. DERIVATION PROCESS OF RESULT 2
This section presents a scheme to obtain Result 2 based on the strong-coupling perturbation theory (see Appendix B).To this end, we focus on how to perform the functional integrals over one-body and two-body potential fields appearing in Equations (B10), (B18), (B27) and (B30)-(B32).

A. One-Body Potential Field (1): Evaluating Equation (B30) in the Saddle-Point Approximation
We see from Equation (B29) that the saddle-point equation δH mf (φ)/δφ| φa=iψ * a = 0 in Equation (B30) gives Substituting Equations (B13) and (B24) into Equation ( 47), we have We can verify that Equation (48) transforms to Equation ( 24) by setting ρ * a (r) = e βµ−ψ * a (r) .Let F mf (−k B T c, ρ * ) be the mean-field free energy defined by where U (v, ρ) has been defined in Equation ( 21) and S id (ρ * ) denotes the sum of ideal gas entropy S id a (ρ * a ) given by Equation ( 17): Plugging Equation (24) into Equations ( 49) and ( 50), we find (see also Appendix C for details of the last equality).The quadratic terms due to fluctuations around the saddle-point path iψ * are written as In the last equality of Equation ( 52), use has been made of the following relation: which is equivalent to the inhomogeneous Ornstein-Zernike Equations ( 29) as confirmed in Appendix D. It is found from Equations ( 51) and ( 52) that the saddle-point approximation of Equation (B30) yields Equations ( B12)-(B14) further imply that Equation ( 54) is transformed into We will use the last line on the rhs of Equation ( 55) as a reference form in evaluating βF (ν) − βF (ν = 0) given by Equation (B31).It follows from Equations (B30) and ( 55) that βΦ (m) where N h /N c is related to the determinant of the matrix, S = c −1 h, as and the matrix element of S is given by ensuring that |S| = |c −1 h| ≥ 0. Replacing r and r ′ by r 0 and r 0 − r, respectively, in Equation ( 58), we have in agreement with the logarithmic term in Equation (32).(1) where the subscript ϕ denotes the following average: according to Equation ( 55) (see Appendix E for the detailed derivation of Equation ( 60)).
C. Two-Body Potential Field: Derivation of Result 2: Rearrangements in the Mean-Field Approximation of Equation ( 63) There are two remaining steps toward obtaining Equation (32): the first step is to evaluate the ν-functional integral given by Equation ( 63), and the second step is to rearrange the interaction energy when adding the last two terms on the rhs of Equation (B3) to U (−k B T c, ρ * ).
First, let us evaluate the ν-field integral given by Equation ( 63) in the mean-field approximation.Equations ( 63) and ( 64) provide the saddle-point equation as follows: giving similar to a closure in the liquid-state theory [62,63] though given correlation functions do not necessarily satisfy any closure, other than the Ornstein-Zernike equation.Substituting Equation (66) into Equation ( 64), we obtain or due to Equations ( 63) and (64).

V. CONCLUDING REMARKS
The generalized FP potential W ( G) as a functional of given TCF h(r) is similar to the original FP potential [65][66][67][68][69][70][71][72][73][74][75] in that both have constraints on inter-replica correlations.The difference is that the generalized FP potential adopts a local order parameter instead of a global order one, the overlap Q (see Equation ( 1)), used in the original FP potential V (Q).Upon reviewing the formulation of W ( G) presented so far, we find two essentials for the fieldtheoretical achievements.The former lies in the variational method described in Appendix B 1, whereas Equation (B8) represents the latter.The details follow: • Unconstrained grand potential mimicking inter-replica correlations: At first, we consider a coupled m-replica system that reproduces a given distribution of the inter-replica TCF h(r) without constraints.We tune the inter-replica interaction potential v(r) to mimic the inter-replica correlations.From evaluating the free-energy functional without constraints in the Gaussian approximation, we obtain the same functional form as the random phase approximation (RPA) in terms of the liquid-state theory [62,63]; however, the density distribution is different.The variational method presented in Appendix B 1 justifies the input of the density distribution given by Equation (27), which converges to that of the Ramakrishnan-Yussouf density functional theory [61] in the limit of m → 1 as demonstrated in Equation (28).
• Evaluating the difference between the constrained and unconstrained grand potentials: Next, we take the freeenergy functional of the unconstrained system as a reference energy.Equation (B8) indicates that the fieldtheoretical formulation focuses on the free energy difference between the constrained and unconstrained freeenergy functionals.The strong-coupling expansion method developed in Appendix B 3 allows us to evaluate this difference in Sections IV B and IV C. Thus, we obtain Equation ( 68), the constraint-associated free energy difference as a functional of inter-replica TCF h(r) and density distribution ρ * (r) determined by the Ramakrishnan-Yussouf theory [61].
Combination of Equations ( 23) and (34) gives the difference between the generalized FP potentials at zero and a finite value of the inter-replica TCFs as follows: The potential difference in Equation ( 72) arises from the entropic cost of localizing the system in an arbitrary state.It is noted, however, that the closure given by Equations ( 37) and ( 38) applies only to Equation ( 72) in a metastable state characterized by h * (r), which is in contrast to the Morita-Hiroike functional covering only the inter-replica TCF that necessarily satisfies the conventional closure [62,63] of the liquid-state theory due to the Legendre-transform-based formalism.That is, the generalized FP potential expressed as Equation ( 72) has a characteristic inherited from the original FP theory, a Landau-type theory relevant to investigate the FEL.Furthermore, Equation (72) represents that our study provides the basis of Ginzburg-Landau-type theory [79] as an extension of Landau-type one: the generalized FP potential W ( G) as a functional of local order parameter h(r) is a natural extension of the FP potential V (Q) as a function of the global order parameter Q.
The stationary Equation (33) reveals that the new closure (37) corresponds to the mean-field equation of W ( G) given by Equation (72).The closure (37) gives the self-consistent Equation (43), similar to the previous one that predicts a dynamical transition [37,56,57]; we need to quantitatively assess the validity of Equation (43) in terms of the dynamical transitions in simulation models.Equation (10) further suggests that we can go beyond the mean-field approximation as is the case with the Ginzburg-Landau-type theory: the greatest advantage of our replica field theory is to systematically improve the self-consistent equation by considering fluctuations of inter-replica correlation field h(r).It remains to be addressed whether the modified self-consistent equation explains the blurring of dynamical transition into a crossover from relaxational to activated dynamics.
There is a caveat, turning our attention to the stability condition on h(r): translational and rotational symmetries are broken in frozen phases.The violation becomes evident by expanding h(r) around that at the uniform density as follows [58][59][60]: We also have a non-perturbative approach to avoid the difficulty using a globally weighted density ρ WD in the interreplica TCF: h(r−r ′ ; ρ * (r)) = h(r−r ′ ; ρ WD ), according to the modified weighted density functional approximation [43,48].Therefore, the functional derivative in Equation (36), or the new closure (37), holds approximately when either neglecting the second and higher-order terms in Equation (73) or finding ρ WD .
The new closure (37) in a metastable state provides the self-consistent Equation (43) for the non-ergodicity parameter f (k).The present field theory has demonstrated the necessity to consider higher-order contributions in the perturbative treatment for obtaining the self-consistent Equation (43) with a kernel containing the triplet DCF [76][77][78]: we obtain Equation ( 43) by adopting the approximate bridge function B(r) = − h 2 (r)/2 beyond the RHNC approximation of B(r) = 0.For comparison, we would like to mention two previous replica approaches to provide the triplet DCF in the self-consistent equation [37,56,57].The first approach considers the perturbative contribution to the replicated HNC functional along the liquid-state theory [37].The Legendre-transform-based method allows us to calculate the third order in h(r) concerning the Morita-Hiroike functional.Meanwhile, the second method considers the third-order term in density difference ρ * (r) − ρ by taking the Ramakrishnan-Yussouf functional of the DFT as a reference form [56,57]. Consequently, both perturbation methods amount to having the triplet DCF in the kernel of the self-consistent equation.This agreement implies the equivalence between the replicated HNC and Ramakrishnan-Yussouf approximations, consistent with the conventional results of the liquid-state theory [62].
Our scheme bears similarity to the Legendre-transform-based theory [33][34][35][36][37][38] rather than the DFT [53][54][55][56][57].However, more elaborate input from the DFT [58][59][60] is also to be investigated, which is particularly necessary to investigate the glass transition in thin polymer films [93,94]; for example, we can improve the Ramakrishnan-Yussouf approximation by performing a variational evaluation beyond the Gaussian approximation (see Appendix B 1).Furthermore, our replica theory has two additional features arising from the field-theoretical treatment of the inter-replica TCF h(r) and the associated two-body interaction potential iν a (r).First, we can systematically consider fluctuations around the mean-field potential field, ν * a (r) = iu(r), given by Equation ( 65), which is the same relation as that of the Legendre-transform-based method [33-38, 64, 86].Second, we can develop the replica field theory to include TCF fluctuations around the metastable field h * (r) as described above.Thus, the present field-theoretical formalism opens up promising avenues to advance studies on the dynamical heterogeneity in terms of the correlation function of TCF fluctuations (i.e., the so-called four-point correlation function [7][8][9][10]) as well as the FEL that includes fluctuations around a metastable state. where Dφ a ≡ Dφ, the Gaussian integration over the ρ-field yields the normalization factor N c written as and Without particles (i.e., ρ = 0), we have consistent with the trivial result U (−k B T c, ρ = 0) = 0 in Equation (B11) (see also Appendix B 5).Meanwhile, the configurational integral represented by Tr provides the perturbative contribution, H 1 (ν, φ), given by e −βH1(ν,φ) = Tr e dr i φ(r)• ρ(r)− m a=2 Defining the functional, Equations (B15)-(B18) indicate that the remaining task is to properly evaluate the perturbative contribution H 1 (ν, φ) in Equation (B16) at strong coupling.

FIG. 1 :
FIG.1:A schematic summary of the main results colored orange.In addition, functional variables are colored green, and underlying potentials blue or gray.

mb=1
Dρ b , Tr ≡ m b=1 Tr b , the matrix elements of v are v ab (r) = 0 (a = b) and v aa

B. Result 2 :
The Constrained Grand Potential Functional of m Replicas in an Inhomogeneous State

C. Result 3 :
New Closure Obtained from the Generalized FP Potential