Loading [MathJax]/extensions/TeX/boldsymbol.js
Research article Special Issues

Maximum degree and minimum degree spectral radii of some graph operations


  • Received: 30 April 2022 Revised: 25 June 2022 Accepted: 04 July 2022 Published: 18 July 2022
  • New results relating to the maximum and minimum degree spectral radii of generalized splitting and shadow graphs have been constructed on the basis of any regular graph, referred as base graph. In particular, we establish the relations of extreme degree spectral radii of generalized splitting and shadow graphs of any regular graph.

    Citation: Xiujun Zhang, Ahmad Bilal, M. Mobeen Munir, Hafiz Mutte ur Rehman. Maximum degree and minimum degree spectral radii of some graph operations[J]. Mathematical Biosciences and Engineering, 2022, 19(10): 10108-10121. doi: 10.3934/mbe.2022473

    Related Papers:

    [1] Mashael Maashi, Mohammed Abdullah Al-Hagery, Mohammed Rizwanullah, Azza Elneil Osman . Deep convolutional neural network-based Leveraging Lion Swarm Optimizer for gesture recognition and classification. AIMS Mathematics, 2024, 9(4): 9380-9393. doi: 10.3934/math.2024457
    [2] Simone Fiori . Coordinate-free Lie-group-based modeling and simulation of a submersible vehicle. AIMS Mathematics, 2024, 9(4): 10157-10184. doi: 10.3934/math.2024497
    [3] Youseef Alotaibi, Veera Ankalu. Vuyyuru . Electroencephalogram based face emotion recognition using multimodal fusion and 1-D convolution neural network (ID-CNN) classifier. AIMS Mathematics, 2023, 8(10): 22984-23002. doi: 10.3934/math.20231169
    [4] Nouf Almutiben, Ryad Ghanam, G. Thompson, Edward L. Boone . Symmetry analysis of the canonical connection on Lie groups: six-dimensional case with abelian nilradical and one-dimensional center. AIMS Mathematics, 2024, 9(6): 14504-14524. doi: 10.3934/math.2024705
    [5] Baiying He, Siyu Gao . The nonisospectral integrable hierarchies of three generalized Lie algebras. AIMS Mathematics, 2024, 9(10): 27361-27387. doi: 10.3934/math.20241329
    [6] David Delphenich . The role of pseudo-hypersurfaces in non-holonomic motion. AIMS Mathematics, 2020, 5(5): 4793-4829. doi: 10.3934/math.2020307
    [7] Junyuan Huang, Xueqing Chen, Zhiqi Chen, Ming Ding . On a conjecture on transposed Poisson n-Lie algebras. AIMS Mathematics, 2024, 9(3): 6709-6733. doi: 10.3934/math.2024327
    [8] Tamilvizhi Thanarajan, Youseef Alotaibi, Surendran Rajendran, Krishnaraj Nagappan . Improved wolf swarm optimization with deep-learning-based movement analysis and self-regulated human activity recognition. AIMS Mathematics, 2023, 8(5): 12520-12539. doi: 10.3934/math.2023629
    [9] He Yuan, Zhuo Liu . Lie n-centralizers of generalized matrix algebras. AIMS Mathematics, 2023, 8(6): 14609-14622. doi: 10.3934/math.2023747
    [10] Guangyu An, Xueli Zhang, Jun He, Wenhua Qian . Characterizations of local Lie derivations on von Neumann algebras. AIMS Mathematics, 2022, 7(5): 7519-7527. doi: 10.3934/math.2022422
  • New results relating to the maximum and minimum degree spectral radii of generalized splitting and shadow graphs have been constructed on the basis of any regular graph, referred as base graph. In particular, we establish the relations of extreme degree spectral radii of generalized splitting and shadow graphs of any regular graph.



    In this paper we compare quenched stochastic two-scale convergence [38] with the notion of stochastic unfolding [30,19], which is equivalent to stochastic two-scale convergence in the mean [6]. In particular, we introduce the concept of stochastic two-scale Young measures to relate quenched stochastic two-scale limits with the mean limit and discuss examples of convex homogenization problems that can be treated with two-scale convergence in the mean, but not conveniently in the quenched setting of two-scale convergence.

    Two-scale convergence has been introduced in [32,1,25] for homogenization problems (partial differential equations or variational problems) with periodic coefficients. The essence of two-scale convergence is that the two-scale limit of an oscillatory sequence captures oscillations that emerge along the sequence and that are to leading order periodic on a definite microscale, typically denoted by ε>0. It is especially well-suited for problems where oscillations of solutions solely stem from prescribed oscillations of the coefficients or the data. For instance, this is the case for equations with a monotone structure or convex variational problems. In contrast, problems that feature pattern formation to leading order (e.g., nonconvex variational problems or singular partial differential equations with non-convex domain) typically cannot be conveniently treated with two-scale convergence. Another well established method for periodic homogenization is periodic unfolding, see [9,35,27,10] as well as [36,3] for the periodic modulation method, which is related. These methods build on an isometric operator---the periodic unfolding (or dilation) operator. It allows us to embed oscillatory sequences into a larger two-scale space and to transform an oscillatory problem into an "unfolded" problem on the two-scale space. The latter often features a better separation of macro- and microscopic properties, which often is convenient for the analysis. We refer to [14,7,28,8,15,24,26] for various interesting applications of this method. Both notions are closely linked, since weak convergence of "unfolded" sequence in the two-scale space is equivalent to weak two-scale convergence, see [5].

    In this paper we are interested in stochastic homogenization, i.e. problems with random coefficients with a stationary distribution. The first stochastic homogenization result has been obtained by Papanicolaou and Varadhan in [33] (and independently by Kozlov [23]) for linear, elliptic equations with stationary and ergodic random coefficients on Rd. In their seminal paper, Papanicolaou and Varadhan introduce a functional analytic framework, which, by now, is the standard way to model random coefficients. We briefly recall it in the special case of convex integral functionals with quadratic growth: Let (Ω,F,P) denote a probability space of parameter fields ωΩ and let τx:ΩΩ, xRd, denote a measure preserving and ergodic group action, see Assumption 2.1 for details. A standard model for a convex, integral functional with a stationary, ergodic, random microstructure on scale ε>0 is then given by the functional Eωε:H1(Q)R{},

    Eωε(u)=QV(τxεω,u(x))f(x)u(x)dx

    where QRd denotes an open and bounded domain, fL2(Q), and V(ω,F) is an integrand that is measurable in ωΩ, convex in FRd, and satisfies a quadratic growth condition. A classical result [11] shows that in the homogenization limit ε0, the functionals Γ-converge to the homogenized functional Ehom:H1(Q)R{}, given by

    Ehom(u)=QVhom(u(x))f(x)u(x)dx,

    where Vhom is a deterministic, convex integrand and characterized by a homogenization formula, see (31) below. There are different natural choices for the topology when passing to this limit:

    ● In the mean setting, minimizers uωε of Eωε, ωΩ, are viewed as random fields (ω,x)uωε(x) in L2(Ω;H1(Q)) and one considers Γ-convergence of the averaged functional L2(Ω;H1(Q))uΩEε(u)dP w.r.t. strong convergence in L2(Ω×Q). In fact, the first result in stochastic homogenization [33] establishes convergence of solutions in this mean sense.

    ● In the quenched setting, one studies the limiting behavior of a minimizer uεH1(Q) of Eωε for fixed ωΩ. One then considers Γ-convergence of Eωε w.r.t. strong convergence in L2(Q) for P-a.a. ωΩ.

    Similarly, two variants of stochastic two-scale convergence have been introduced as generalizations of periodic two-scale convergence (for the sake of brevity, we restrict the following review to the Hilbert-space case p=2, and note that the following extends to Lp(Ω×Q) with p(1,)):

    ● In [6,2] the mean variant has been introduced as follows: We say that a sequence of random fields (uε)L2(Ω×Q) stochastically two-scale converges in the mean to uL2(Ω×Q), if

    limε0Ω×Quε(ω,x)φ(τxεω,x)dP(ω)dx=Ω×Qu(ω,x)φ(ω,x)dP(ω)dx, (1)

    for all admissible test functions φL2(Ω×Q), see Remark 1 for details.

    ● More recently, Zhikov and Pyatnitskii introduced in [38] a quenched variant: We say that a sequence (uε)L2(Q) quenched stochastically two-scale converges to uL2(Ω×Q) w.r.t. to a fixed parameter field ω0Ω, if

    limε0Quε(x)φ(τxεω0,x)dx=Ω×Qu(ω,x)φ(ω,x)dP(ω)dx,

    for all admissible test functions φL2(Ω×Q). Note that the two-scale limit u a priori depends on ω0. In fact, in [37] (see also [16]) quenched two-scale convergence has been introduced in a very general setting that includes the case of integration against random, rapidly oscillating measures, which naturally emerge when describing coefficients defined relative to random geometries. In this work, we restrict our considerations to the simplest case where the random measure is the Lebesgue measure.

    Similarly to the periodic case, stochastic two-scale convergence in the mean can be rephrased with help of a transformation operator, see [34,19,30], where the stochastic unfolding operator Tε:L2(Ω×Q)L2(Ω×Q),

    Tεu(ω,x)=u(τxεω,x), (2)

    has been introduced. As in the periodic case, it is a linear isometry and it turns out that for a bounded sequence (uε)L2(Ω×Q), stochastic two-scale convergence in the mean is equivalent to weak convergence of the unfolded sequence Tεuε. As we demonstrate below in Section 4.1, the stochastic unfolding method leads to a very economic and streamlined analysis of convex homogenization problems. Moreover, it allows us to derive two-scale functionals of the form E(u,χ)=ΩQV(ω,u(x)+χ(ω,x))dxdP as a Γ-limit of Eε, see Theorem 4.1 for details. In contrast to the periodic case, where the unfolding operator is an isometry from L2(Rd) to L2(Y×Rd) (with Y denoting the unit torus), in the random case it is not possible to interpret (2) as a continuous operator from L2(Q) to L2(Ω×Q). Therefore, quenched two-scale convergence cannot be characterized via stochastic unfolding directly.

    In the present paper we compare the different notions of stochastic two-scale convergence. Although the mean and quenched notion of two-scale convergence look quite similar, it is non-trivial to relate both. As a main result, we introduce stochastic two-scale Young measures as a tool to compare quenched and mean limits, see Theorem 3.12. The construction invokes a metric characterization of quenched stochastic two-scale convergence, which is a tool of independent interest, see Lemma 3.6. As an application we demonstrate how to lift a mean two-scale homogenization result to a quenched statement, see Section 4.3. Moreover, we present two examples that can only be conveniently treated with the mean notion of two-scale convergence. In the first example, see Section 4.1, the assumption of ergodicity is dropped (as it is natural in the context of periodic representative volume approximation schemes). In the second example we consider a model that invokes a mean field interaction in form of a variance-type regularization of a convex integral functional with degenerate growth, see Section 4.2.

    Structure of the paper. In the following section we present the standard setting for stochastic homogenization. In Section 3 we provide the main properties of the stochastic unfolding method, present the most important facts about quenched two-scale convergence and present our results about Young measures. In Section 4 we present examples of stochastic homogenization and applications of the methods developed in this paper.

    In the following we briefly recall the standard setting for stochastic homogenization. Throughout the entire paper we assume the following:

    Assumption 2.1. Let (Ω,F,P) be a complete and separable probability space. Let τ={τx}xRd denote a group of invertible measurable mappings τx:ΩΩ such that:

    (i)(Group property). τ0=Id and τx+y=τxτy for all x,yRd.

    (ii)(Measure preservation). P(τxE)=P(E) for all EF and xRd.

    (iii)(Measurability). (ω,x)τxω is (FL(Rd),F)-measurable, where L(Rd) denotes the Lebesgue σ-algebra.

    We write to denote the expectation ΩdP. By the separability assumption on the measure space it follows that Lp(Ω) is separable for p1. The proof of the following lemma is a direct consequence of Assumption 2.1, thus we omit it.

    Lemma 2.2 (Stationary extension). Let φ:ΩR be F-measurable. Let QRd be open and denote by L(Q) the corresponding Lebesgue σ-algebra. Then Sφ:Ω×QR, Sφ(ω,x):=φ(τxω) defines an FL(Q)-measurable function – called the stationary extension of φ. Moreover, if Q is bounded, for all 1p< the map S:Lp(Ω)Lp(Ω×Q) is a linear injection satisfying

    SφLp(Ω×Q)=|Q|1pφLp(Ω).

    We say (Ω,F,P,τ) is ergodic ( is ergodic), if

     every shift invariant AF (i.e. τxA=A for all xRd) satisfies P(A){0,1}.

    In this case the celebrated Birkhoff's ergodic theorem applies, which we recall in the following form:

    Theorem 2.3(Birkhoff's ergodic Theorem [12,Theorem 10.2.II]). Let be ergodic and φ:ΩR be integrable. Then for P-a.a. ωΩ it holds: Sφ(ω,) is locally integrable and for all open, bounded sets QRd we have

    limε0QSφ(ω,xε)dx=|Q|φ. (3)

    Furthermore, if φLp(Ω) with 1p, then for P-a.a. ωΩ it holds: Sφ(ω,)Lploc(Rd), and provided p< it holds Sφ(ω,ε)φ weakly in Lploc(Rd) as ε0.

    Stochastic gradient. For p(1,) consider the group of isometric operators {Ux:xRd} on Lp(Ω) defined by Uxφ(ω)=φ(τxω). This group is strongly continuous (see [22,Section 7.1]). For i=1,...,d, we consider the 1-parameter group of operators {Uhei:hR} and its infinitesimal generator Di:DiLp(Ω)Lp(Ω)

    Diφ=limh0Uheiφφh,

    which we refer to as stochastic derivative. Di is a linear and closed operator and its domain Di is dense in Lp(Ω). We set W1,p(Ω)=di=1Di and define for φW1,p(Ω) the stochastic gradient as Dφ=(D1φ,...,Ddφ). In this way, we obtain a linear, closed and densely defined operator D:W1,p(Ω)Lp(Ω)d, and we denote by

    Lppot(Ω):=¯R(D)Lp(Ω)d (4)

    the closure of the range of D in Lp(Ω)d. We denote the adjoint of D by D:DLq(Ω)dLq(Ω) where here and below q:=pp1 denotes the dual exponent. It is a linear, closed and densely defined operator (D is the domain of D). We define the subspace of shift invariant functions in Lp(Ω) by

    Lpinv(Ω)={φLp(Ω):Uxφ=φfor all xRd},

    and denote by Pinv:Lp(Ω)Lpinv(Ω) the conditional expectation with respect to the σ-algebra of shift invariant sets {AF:τxA=A for all xRd}. Pinv a contractive projection and for p=2 it coincides with the orthogonal projection onto L2inv(Ω). The following well-known equivalence holds:

     is ergodic  Lpinv(Ω)R  Pinvf=f.

    Random fields. We introduce function spaces for functions defined on Ω×Q as follows: For closed subspaces XLp(Ω) and YLp(Q), we denote by XY the closure of

    XaY:={ni=1φiηi:φiX,ηiY,nN}

    in Lp(Ω×Q). Note that in the case X=Lp(Ω) and Y=Lp(Q), we have XY=Lp(Ω×Q). Up to isometric isomorphisms, we may identify Lp(Ω×Q) with the Bochner spaces Lp(Ω;Lp(Q)) and Lp(Q;Lp(Ω)). Slightly abusing the notation, for closed subspaces XLp(Ω) and YW1,p(Q), we denote by XY the closure of

    XaY:={ni=1φiηi:φiX,ηiY,nN}

    in Lp(Ω;W1,p(Q)). In this regard, we may identify uLp(Ω)W1,p(Q) with the pair (u,u)Lp(Ω×Q)1+d. We mostly focus on the space Lp(Ω×Q) and the above notation is convenient for keeping track of its various subspaces.

    In the following we first discuss two notions of stochastic two-scale convergence and their connection through Young measures. In particular, Section 3.1 is devoted to the introduction of the stochastic unfolding operator and its most important properties. In Section 3.2 we discuss quenched two-scale convergence and its properties. Section 3.3 presents the results about Young measures.

    In the following we briefly introduce the stochastic unfolding operator and provide its main properties, for the proofs and detailed studies we refer to [19,30,31,34].

    Lemma 3.1([19,Lemma 3.1]). Let ε>0, 1<p<, q=pp1, and QRd be open. There exists a unique linear isometric isomorphism

    Tε:Lp(Ω×Q)Lp(Ω×Q)

    such that

    uLp(Ω)aLp(Q):(Tεu)(ω,x)=u(τxεω,x)a.e.inΩ×Q.

    Moreover, its adjoint is the unique linear isometric isomorphism Tε:Lq(Ω×Q)Lq(Ω×Q) that satisfies (Tεu)(ω,x)=u(τxεω,x) a.e. in Ω×Q for all uLq(Ω)aLq(Q), q:=pp1.

    Definition 3.2 (Unfolding and two-scale convergence in the mean). The operator Tε:Lp(Ω×Q)Lp(Ω×Q) in Lemma 3.1 is called the stochastic unfolding operator. We say that a sequence (uε)Lp(Ω×Q) weakly (strongly) two-scale converges in the mean in Lp(Ω×Q) to uLp(Ω×Q) if (as ε0)

    Tεuεu weakly (strongly) in Lp(Ω×Q).

    In this case we write uε2u (uε2u) in Lp(Ω×Q).

    Remark 1 (Equivalence to stochastic two-scale convergence in the mean). Stochastic two-scale convergence in the mean was introduced in [6]. In particular, it is said that a sequence of random fields uεLp(Ω×Q) stochastically two-scale converges in the mean if

    limε0Quε(ω,x)φ(τxεω,x)dx=Qu(ω,x)φ(ω,x)dx, (5)

    for any φLq(Ω×Q), q=pp1, that is admissible, i.e., in the sense that the transformation (ω,x)φ(τxεω,x) is well-defined. For a bounded sequence uεLp(Ω×Q), (5) is equivalent to Tεuεu weakly in Lp(Ω×Q), i.e., to weak stochastic two-scale convergence in the mean. Indeed, with help of Tε (and its adjoint) we might rephrase the integral on the left-hand side in (5) as

    Quε(Tεφ)dx=Q(Tεuε)φdx, (6)

    which proves the equivalence.

    We summarize some of the main properties:

    Proposition 1 (Main properties). Let p(1,), q=pp1 and QRd be open.

    (i)(Compactness, [19,Lemma 3.4].) If lim supε0uεLp(Ω×Q)<, then there exists a subsequence ε and uLp(Ω×Q) such that uε2u in Lp(Ω×Q).

    (ii)(Limits of gradients, [19,Proposition 3.7]) Let (uε) be a bounded sequence in Lp(Ω)W1,p(Q). Then, there exist uLpinv(Ω)W1,p(Q) and χLppot(Ω)Lp(Q) such that (up to a subsequence)

    uε2uinLp(Ω×Q),uε2u+χinLp(Ω×Q)d. (7)

    If, additionally, is ergodic, then u=Pinvu=uW1,p(Q) and uεu weakly in W1,p(Q).

    (iii)(Recovery sequences, [19,Lemma 4.3]) Let uLpinv(Ω)W1,p(Q) and χLppot(Ω)Lp(Q). There exists uεLp(Ω)W1,p(Q) such that

    uε2u,uε2u+χinLp(Ω×Q).

    If additionally uLpinv(Ω)W1,p0(Q), we have uεLp(Ω)W1,p0(Q).

    In this section, we recall the concept of quenched stochastic two-scale convergence (see [38,16]). The notion of quenched stochastic two-scale convergence is based on the individual ergodic theorem, see Theorem 2.3. We thus assume throughout this section that

    isergodic.

    Moreover, throughout this section we fix exponents p(1,), q:=pp1, and an open and bounded domain QRd. We denote by (Bp,Bp) the Banach space Lp(Ω×Q) and the associated norm, and we write (Bp) for the dual space. For the definition of quenched two-scale convergence we need to specify a suitable space of test-functions in Bq that is countably generated. To that end we fix sets DΩ and DQ such that

    DΩ is a countable set of bounded, measurable functions on (Ω,F) that contains the identity 1Ω1 and is dense in L1(Ω) (and thus in Lr(Ω) for any 1r<).

    DQC(¯Q) is a countable set that contains the identity 1Q1 and is dense in L1(Q) (and thus in Lr(Q) for any 1r<).

    We denote by

    A:={φ(ω,x)=φΩ(ω)φQ(x):φΩDΩ,φQDQ}

    the set of simple tensor products (a countable set), and by D0 the Q-linear span of A, i.e.

    D0:={mj=1λjφj:mN,λ1,,λmQ,φ1,,φmA}.

    We finally set

    D:=spanA=spanD0and¯D:=span(DQ)

    (the span of DQ seen as a subspace of D), and note that D and D0 are dense subsets of Bq, while the closure of ¯D in Bq is isometrically isomorphic to Lq(Q). Let us anticipate that D serves as our space of test-functions for stochastic two-scale convergence. As opposed to two-scale convergence in the mean, "quenched" stochastic two-scale convergence is defined relative to a fixed "admissible" realization ω0Ω. Throughout this section we denote by

    Ω0the set of admissible realizations;

    it is a set of full measure determined by the following lemma:

    Lemma 3.3. There exists a measurable set Ω0Ω with P(Ω0)=1 s.t. for all φ,φA, all ω0Ω0, and r{p,q} we have with (Tεφ)(ω,x):=φ(τxεω,x),

    lim supε0(Tεφ)(ω0,)Lr(Q)φBrandlimε0QTε(φφ)(ω0,x)dx=Q(φφ)(ω0,x)dx.

    Proof. This is a simple consequence of Theorem 2.3 and the fact that A is countable.

    For the rest of the section Ω0 is fixed according to Lemma 3.3.

    The idea of quenched stochastic two-scale convergence is similar to periodic two-scale convergence: We associate with a bounded sequence (uε)Lp(Q) and ω0Ω0, a sequence of linear functionals (uε) defined on D. We can pass (up to a subsequence) to a pointwise limit U, which is again a linear functional on D and which (thanks to Lemma 3.3) can be uniquely extended to a bounded linear functional on Bq. We then define the weak quenched ω0-two-scale limit of (uε) as the Riesz-representation uBp of U(Bq).

    Definition 3.4 (quenched two-scale limit, cf. [38,17]). Let (uε) be a sequence in Lp(Q), and let ω0Ω0 be fixed. We say that uε converges (weakly, quenched) ω0-two-scale to uBp, and write uε2ω0u, if the sequence uε is bounded in Lp(Q), and for all φD we have

    limε0Quε(x)(Tεφ)(ω0,x)dx=ΩQu(x,ω)φ(ω,x)dxdP(ω). (8)

    Lemma 3.5 (Compactness). Let (uε) be a bounded sequence in Lp(Q) and ω0Ω0. Then there exists a subsequence (still denoted by ε) and uBp such that uε2ω0u and

    uBplim infε0uεLp(Q), (9)

    and uεu weakly in Lp(Q).

    (For the proof see Section 3.2.1).

    For our purpose it is convenient to have a metric characterization of two-scale convergence.

    Lemma 3.6 (Metric characterization). (i)Let {φj}jN denote an enumeration of the countable set {φφBq:φD0}. The vector space Lin(D):={U:DRlinear} endowed with the metric

    d(U,V;Lin(D)):=jN2j|U(φj)V(φj)||U(φj)V(φj)|+1

    is complete and separable.

    (ii)Let ω0Ω0. Consider the maps

    Jω0ε:Lp(Q)Lin(D),(Jω0εu)(φ):=Qu(x)(Tεφ)(ω0,x)dx,J0:BpLin(D),(J0u)(φ):=Quφ.

    Then for any bounded sequence uε in Lp(Q) and any uBp we have uε2ω0u if and only if Jω0εuεJ0u in Lin(D).

    (For the proof see Section 3.2.1).

    Remark 2. Convergence in the metric space (Lin(D),d(,,Lin(D)) is equivalent to pointwise convergence. (Bq) is naturally embedded into the metric space by means of the restriction J:(Bq)Lin(D), JU=U|D. In particular, we deduce that for a bounded sequences (Uk) in (Bq) we have UkU if and only if JUkJU in the metric space. Likewise, Bp (resp. Lp(Q)) can be embedded into the metric space Lin(D) via J0 (resp. Jω0ε with ε>0 and ω0Ω0 arbitrary but fixed), and for a bounded sequence (uk) in Bp (resp. Lp(Q)) weak convergence in Bp (resp. Lp(Q)) is equivalent to convergence of (J0uk) (resp. (Jω0εuk)) in the metric space.

    Lemma 3.7 (Strong convergence implies quenched two-scale convergence). Let (uε) be a strongly convergent sequence in Lp(Q) with limit uLp(Q). Then for all ω0Ω0 we have uε2ω0u.

    (For the proof see Section 3.2.1).

    Definition 3.8 (set of quenched two-scale cluster points). For a bounded sequence (uε) in Lp(Q) and ω0Ω0 we denote by CP(ω0,(uε)) the set of all ω0-two-scale cluster points, i.e. the set of uBp with J0uk=1¯{Jω0εuε:ε<1k} where the closure is taken in the metric space (Lin(D),d(,;Lin(D))).

    We conclude this section with two elementary results on quenched stochastic two-scale convergence:

    Lemma 3.9 (Approximation of two-scale limits). Let uBp.Then for all ω0Ω0, there exists a sequenceuεLp(Q) such that uε2ω0u as ε0.

    (For the proof see Section 3.2.1).

    Similar to the slightly different setting in [17] one can prove the following result:

    Lemma 3.10 (Two-scale limits of gradients). Let (uε) be a sequence in W1,p(Q) and ω0Ω0. Then there exist a subsequence (not relabeled) and functions uW1,p(Q) and χLppot(Ω)Lp(Q) such that uεu weakly in W1,p(Q) and

    uε2ω0uanduε2ω0u+χasε0.

    Proof of Lemma 3.5. Set C0:=lim supε0uεLp(Q) and note that C0<. By passing to a subsequence (not relabeled) we may assume that C0=lim infε0uεLp(Q). Fix ω0Ω0. Define linear functionals UεLin(D) via

    Uε(φ):=Quε(x)(Tεφ)(ω0,x)dx.

    Note that for all φA, (uε(φ)) is a bounded sequence in R. Indeed, by Hölder's inequality and Lemma 3.3,

    lim supε0|uε(φ)|lim supε0uεLp(Q)Tεφ(ω0,)Lq(Q)C0φBq. (10)

    Since A is countable we can pass to a subsequence (not relabeled) such that uε(φ) converges for all φA. By linearity and since D=span(A), we conclude that uε(φ) converges for all φD, and U(φ):=limε0uε(φ) defines a linear functional on D. In view of (10) we have |U(φ)|C0φBq, and thus U admits a unique extension to a linear functional in (Bq). Let uBp denote its Riesz-representation. Then uε2ω0u, and

    uBp=U(Bq)C0=lim infε0uεLp(Q).

    Since 1ΩDΩ we conclude that for all φQDQ we have

    Quε(x)φQ(x)dx=uε(1ΩφQ)U(1ΩφQ)=Qu(ω,x)φQ(x)dx=Qu(x)φQ(x)dx.

    Since (uε) is bounded in Lp(Q), and DQLp(Q) is dense, we conclude that uεu weakly in Lp(Q).

    Proof of Lemma 3.6. We use the following notation in this proof A1:={φφBq:φD0}.

    (i) Argument for completeness: If (Uj) is a Cauchy sequence in Lin(D), then for all φA1, (Uj(φ)) is a Cauchy sequence in R. By linearity of the Uj's this implies that (Uj(φ)) is Cauchy in R for all φD. Hence, UjU pointwise in D and it is easy to check that U is linear. Furthermore, UjU pointwise in A1 implies UjU in the metric space.

    Argument for separability: Consider the (injective) map J:(Bq)Lin(D) where J(U) denotes the restriction of U to D. The map J is continuous, since for all U,V(Bq) and φA1 we have |(JU)(φ)(JV)(φ)|UV(Bq)φBq=UV(Bq) (recall that the test functions in A1 are normalized). Since (Bq) is separable (as a consequence of the assumption that F is countably generated), it suffices to show that the range R(J) of J is dense in Lin(D). To that end let ULin(D). For kN we denote by Uk(Bq) the unique linear functional that is equal to U on the the finite dimensional (and thus closed) subspace span{φ1,,φk}Bq (where {φj} denotes the enumeration of A1), and zero on the orthogonal complement in Bq. Then a direct calculation shows that d(U,J(Uk);Lin(D))j>k2j=2k. Since kN is arbitrary, we conclude that R(J)Lin(D) is dense.

    (ii) Let uε denote a bounded sequence in Lp(Q) and uBp. Then by definition, uε2ω0u is equivalent to Jω0εuεJ0u pointwise in D, and the latter is equivalent to convergence in the metric space Lin(D).

    Proof of Lemma 3.7. This follows from Hölder's inequality and Lemma 3.3, which imply for all φA the estimate

    limsupε0Q|(uε(x)u(x))Tεφ(ω0,x)|dxlimsupε0(uεuLp(Q)(Q|Tεφ(ω0,x)|qdx)1q)=0.

    Proof of Lemma 3.9. Since D (defined as in Lemma 3.6) is dense in Bp, for any δ>0 there exists vδD0 with uvδBpδ. Define vδ,ε(x):=Tεvδ(ω0,x). Let φD. Since vδ and φ (resp. vδφ) are by definition linear combinations of functions (resp. products of functions) in A, we deduce from Lemma 3.3 that (vδ,ε)ε is bounded in Lp(Q), and that

    Qvδ,εTεφ(ω0,x)=QTε(vδφ)(ω0,x)Qvδφ.

    By appealing to the metric characterization, we can rephrase the last convergence statement as d(Jω0εvδ,ε,J0vδ;Lin(D))0. By the triangle inequality we have

    d(Jω0εvδ,ε,J0u;Lin(D))d(Jω0εvδ,ε,J0vδ;Lin(D))+d(J0vδ,J0u;Lin(D)).

    The second term is bounded by vδuBpδ, while the first term vanishes for ε0. Hence, there exists a diagonal sequence uε:=vδ(ε),ε (bounded in Lp(Q)) such that there holds d(Jω0εuε,J0u;Lin(D))0. The latter implies uε2ω0u by Lemma 3.6.

    In this section we establish a relation between quenched two-scale convergence and two-scale convergence in the mean (in the sense of Definition 3.2). The relation is established by Young measures: We show that any bounded sequence uε in Bp – viewed as a functional acting on test-functions of the form Tεφ – generates (up to a subsequence) a Young measure on Bp that (a) concentrates on the quenched two-scale cluster points of uε, and (b) allows to represent the two-scale limit (in the mean) of uε. In entire Section 3.3 we assume that

     is ergodic.

    Also, throughout this section we fix exponents p(1,), q:=pp1, and an open and bounded domain QRd. Furthermore, we frequently use the objects and notations introduced in Section 3.2.

    Definition 3.11. We say \boldsymbol{\nu}: = \left\{ \nu_{\omega}\right\} _{\omega\in\Omega} is a Young measure on {\mathscr{B}}^p , if for all \omega\in\Omega , \nu_\omega is a Borel probability measure on {\mathscr{B}}^p (equipped with the weak topology) and

    \omega\mapsto\nu_{\omega}(B)\quad\mbox{is }\mbox{measurable for all }B\in {\mathcal{B}}( {\mathscr{B}}^p),

    where {\mathcal{B}}( {\mathscr{B}}^p) denotes the Borel- \sigma -algebra on {\mathscr{B}}^p (equipped with the weak topology).

    Theorem 3.12. Let u_{ \varepsilon} denote a bounded sequence in {\mathscr{B}}^p .Then there exists a subsequence (still denoted by \varepsilon ) and a Young measure \boldsymbol{\nu} on {\mathscr{B}}^p such that for all \omega_0\in\Omega_0 ,

    \nu_{\omega_0}\;\mathit{\mbox{is concentrated on}}\;{\mathscr{C\!P}}\left(\omega_0,\big(u_{ \varepsilon}(\omega_0,\cdot)\big)\right),

    and

    \liminf\limits_{ \varepsilon\to0}\Vert u_{ \varepsilon}\Vert_{ {\mathscr{B}}^p}^{p}\geq \int_{\Omega}\left(\int_{ {\mathscr{B}}^p}\left\Vert v\right\Vert _{ {\mathscr{B}}^p}^{p}\,d\nu_{\omega}(v)\right)\,dP(\omega).

    Moreover, we have

    u_{ \varepsilon} {\stackrel{2}{\rightharpoonup}} u\qquad{{where}}\;u: = \int_{\Omega}\int_{ {\mathscr{B}}^p}v\, d\nu_{\omega}(v)dP(\omega).

    Finally, if there exists \hat u:\Omega\to {\mathscr{B}}^p measurable and \nu_{\omega} = \delta_{\hat u(\omega)} for P -a.a. \omega\in\Omega , then up to extraction of a further subsequence (still denoted by \varepsilon ) we have

    u_{ \varepsilon}(\omega){\stackrel{2}{\rightharpoonup}_{{\omega}}}\hat u(\omega)\qquad{{for\; P -a.a.\; \omega\in\Omega }}.

    (For the proof see Section 3.3.1).

    In the opposite direction we observe that quenched two-scale convergence implies two-scale convergence in the mean in the following sense:

    Lemma 3.13. Consider a family \{(u_{ \varepsilon}^\omega)\}_{\omega\in\Omega} of sequences (u^\omega_{ \varepsilon}) in L^p(Q) and suppose that:

    (i)There exists u\in {\mathscr{B}}^p s.t. for P -a.a. \omega\in\Omega we have u_{ \varepsilon}^\omega{\stackrel{2}{\rightharpoonup}_{{\omega}}}u .

    (ii)There exists a sequence (\tilde u_{ \varepsilon}) in {\mathscr{B}}^p s.t. u_{ \varepsilon}^\omega(x) = \tilde u_{ \varepsilon}(\omega,x) for a.a. (\omega,x)\in\Omega\times Q .

    (iii)There exists a bounded sequence (\chi_{ \varepsilon}) in L^p(\Omega) such that \|u^\omega_{ \varepsilon}\|_{L^p(Q)}\leq\chi_{ \varepsilon}(\omega) for a.a. \omega\in\Omega .

    Then \tilde u_{ \varepsilon} \overset{2}{\rightharpoonup} u weakly two-scale (in the mean).

    (For the proof see Section 3.3.1).

    To compare homogenization of convex integral functionals w.r.t. stochastic two-scale convergence in the mean and in the quenched sense, we appeal to the following result:

    Lemma 3.14. Let h:\,\Omega\times Q\times \mathbb{R}^{d}\to \mathbb{R} be such that for all \xi\in \mathbb{R}^d , h(\cdot,\cdot,\xi) is \mathcal F\otimes {\mathcal{B}}( \mathbb{R}^{d}) -measurable and for a.a. (\omega,x)\in\Omega\times Q , h(\omega,x,\cdot) is convex. Let (u_{ \varepsilon}) denote a bounded sequence in {\mathscr{B}}^p that generates a Young measure \boldsymbol{\nu} on {\mathscr{B}}^p in the sense of Theorem 3.12.Suppose that h_{ \varepsilon}:\Omega\to \mathbb{R} , h_{ \varepsilon}(\omega): = -\int_Q\min\big\{0,h(\tau_{\frac{x}{ \varepsilon}}\omega,x,u_{ \varepsilon}(\omega,x))\big\}\,dx is uniformly integrable. Then

    \begin{array}{*{20}{c}} {\mathop {\lim \inf }\limits_{\varepsilon \to 0} \int_\Omega {\int_Q h } ({\tau _{\frac{x}{\varepsilon }}}\omega ,x,{u_\varepsilon }(\omega ,x)){\kern 1pt} dx{\kern 1pt} dP(\omega )}\\ { \ge \int_\Omega {\int_{{{\mathscr B}^p}} {\left( {\int_\Omega {\int_Q h } (\tilde \omega ,x,v(\tilde \omega ,x)){\kern 1pt} dx{\kern 1pt} dP(\tilde \omega )} \right)} } {\kern 1pt} d{\nu _\omega }(v){\kern 1pt} dP(\omega ).} \end{array} (11)

    (For the proof see Section 3.3.1).

    Remark 3. In [18,Lemma 5.1] it is shown that h satisfying the assumptions of Lemma 3.14 satisfies the following property: For P -a.a. \omega_0\in\Omega_0 we have: For any sequence (u_{ \varepsilon}) in L^p(Q) it holds

    \begin{equation} u_{ \varepsilon}{\stackrel{2}{\rightharpoonup}_{{\omega_0}}}u\quad\Rightarrow \quad \liminf\limits_{ \varepsilon\to0}\int_{Q}h(\tau_{\frac{x}{ \varepsilon}}\omega_0,x,u_{ \varepsilon}(x))dx\geq\int_{\Omega}\int_{Q}h(\omega,x,u(\omega,x))\,dx\,dP(\omega). \end{equation} (12)

    We first recall some notions and results of Balder's theory for Young measures [4]. Throughout this section {\mathscr M} is assumed to be a separable, complete metric space with metric d(\cdot,\cdot; {\mathscr M}) .

    Definition 3.15. ● We say a function s:\Omega\to {\mathscr M} is measurable, if it is \mathcal F-\mathcal B( {\mathscr M}) -measurable where \mathcal B( {\mathscr M}) denotes the Borel- \sigma -algebra in {\mathscr M} .

    ● A function h:\Omega\times {\mathscr M}\to(-\infty,+\infty] is called a normal integrand, if h is \mathcal F\otimes\mathcal B( {\mathscr M}) -measurable, and for all \omega\in\Omega the function h(\omega,\cdot): {\mathscr M}\to(-\infty,+\infty] is lower semicontinuous.

    ● A sequence s_{ \varepsilon} of measurable functions s_{ \varepsilon}:\Omega\to {\mathscr M} is called tight, if there exists a normal integrand h such that for every \omega\in\Omega the function h(\omega,\cdot) has compact sublevels in {\mathscr M} and \limsup_{ \varepsilon\to 0}\int_\Omega h(\omega,s_{ \varepsilon}(\omega))\,dP(\omega)<\infty .

    ● A Young measure in {\mathscr M} is a family \boldsymbol{\mu}: = \left\{ \mu_{\omega}\right\} _{\omega\in\Omega} of Borel probability measures on {\mathscr M} such that for all B\in\mathcal B( {\mathscr M}) the map \Omega\ni \omega\mapsto \mu_\omega(B)\in \mathbb{R} is \mathcal F -measurable.

    Theorem 3.16.([4,Theorem I]). Let s_{ \varepsilon}:\,\Omega\to {\mathscr M} denote a tight sequence of measurable functions. Then there exists a subsequence, still indexed by \varepsilon , and a Young measure {\boldsymbol\mu}:\Omega\to {\mathscr M} such that for every normal integrand h:\,\Omega\times {\mathscr M}\rightarrow (-\infty,+\infty] we have

    \begin{equation} \liminf\limits_{ \varepsilon\to0}\int_{\Omega}h(\omega,s_{ \varepsilon}(\omega))\,dP(\omega)\geq\int_{\Omega}\int_{ {\mathscr M}}h(\omega,\xi) d\mu_{\omega}(\xi)dP(\omega)\,, \end{equation} (13)

    provided that the negative part h^-_ \varepsilon(\cdot) = |\min\{0,h(\cdot,s_{ \varepsilon}(\cdot))\}| is uniformly integrable.Moreover, for P -a.a. \omega\in\Omega_0 the measure \mu_\omega is supported in the set of all cluster points of s_{ \varepsilon}(\omega) , i.e. in \bigcup_{k = 1}^\infty\overline{\{s_{ \varepsilon}(\omega)\,:\, \varepsilon<\frac{1}{k}\}} (where the closure is taken in {\mathscr M} ).

    In order to apply the above theorem we require an appropriate metric space in which two-scale convergent sequences and their limits embed:

    Lemma 3.17. (i)>We denote by {\mathscr M} the set of all triples (U, \varepsilon,r) with U\in{\mbox{Lin}}( {\mathscr D}) , \varepsilon\geq 0 , r\geq 0 . {\mathscr M} endowed with the metric

    \begin{equation*} d((U_1, \varepsilon_1,r_1),(U_2, \varepsilon_2,r_2); {\mathscr M}): = d(U_1,U_2;{\mbox{Lin}}( {\mathscr D}))+| \varepsilon_1- \varepsilon_2|+|r_1-r_2| \end{equation*}

    is a complete, separable metric space.

    (ii)For \omega_0\in\Omega_0 we denote by {\mathscr M}^{\omega_0} the set of all triples (U, \varepsilon,r)\in {\mathscr M} such that

    \begin{equation} U = \begin{cases} J^{\omega_0}_ \varepsilon u&\;\mathit{\text{for some}}\;u\in L^p(Q)\;\mathit{\text{with}}\;\|u\|_{L^p(Q)}\leq r\;\mathit{\text{in the case}}\; \varepsilon > 0,\\ J_0 u&\;\mathit{\text{for some}}\;u\in {\mathscr{B}}^p\;\mathit{\text{with}}\;\|u\|_{ {\mathscr{B}}^p}\leq r\;\mathit{\text{in the case}}\; \varepsilon = 0. \end{cases} \end{equation} (14)

    Then {\mathscr M}^{\omega_0} is a closed subspace of {\mathscr M} .

    (iii)Let \omega_0\in\Omega_0 , and (U, \varepsilon,r)\in {\mathscr M}^{\omega_0} . Then the function u in the representation (14) of U is unique, and

    \begin{equation} \begin{cases} \|u\|_{L^p(Q)} = \sup\limits_{\varphi\in\overline{ {\mathscr D}},\ \|\varphi\|_{ {\mathscr{B}}^q}\leq 1}|U(\varphi)|&\;\mathit{\text{if}}\; \varepsilon > 0,\\ \|u\|_{ {\mathscr{B}}^p} = \sup\limits_{\varphi\in {\mathscr D},\ \|\varphi\|_{ {\mathscr{B}}^q}\leq 1}|U(\varphi)|&\;\mathit{\text{if}}\; \varepsilon = 0. \end{cases} \end{equation} (15)

    (iv)For \omega_0\in\Omega_0 the function \|\cdot\|_{\omega_0}: {\mathscr M}^{\omega_0}\to[0,\infty) ,

    \begin{equation*} \|(U, \varepsilon,r)\|_{\omega_0}: = \begin{cases} \big(\sup\limits_{\varphi\in\overline{ {\mathscr D}},\ \|\varphi\|_{ {\mathscr{B}}^q}\leq 1}|U(\varphi)|^p+ \varepsilon+r^p\big)^{\frac{1}{p}}&\;\mathit{\text{if}}\;(U, \varepsilon,r)\in {\mathscr M}^{\omega_0},\, \varepsilon > 0,\\ \big(\sup\limits_{\varphi\in {\mathscr D},\ \|\varphi\|_{ {\mathscr{B}}^q}\leq 1}|U(\varphi)|^p+r^p\big)^\frac1p&\;\mathit{\text{if}}\;(U, \varepsilon,r)\in {\mathscr M}^{\omega_0},\, \varepsilon = 0,\\ \end{cases} \end{equation*}

    is lower semicontinuous on {\mathscr M}^{\omega_0} .

    (v)For all (u, \varepsilon) with u\in L^p(Q) and \varepsilon>0 we have s: = (J^{\omega_0}_ \varepsilon u, \varepsilon,\|u\|_{L^p(Q)})\in {\mathscr M}^{\omega_0} and \|s\|_{\omega_0} = \big(2\|u\|_{L^p(Q)}^p+ \varepsilon\big)^\frac1p . Likewise, for all (u, \varepsilon) with u\in {\mathscr{B}}^p and \varepsilon = 0 we have s = (J_0u, \varepsilon,\|u\|_{ {\mathscr{B}}^p}) and \|s\|_{\omega_0} = 2^\frac1p\|u\|_{ {\mathscr{B}}^p} .

    (vi)For all R<\infty the set \{(U, \varepsilon,r)\in {\mathscr M}^{\omega_0}\,:\,\|(U, \varepsilon,r)\|_{\omega_0}\leq R\} is compact in {\mathscr M} .

    (vii)Let \omega_0\in\Omega_0 and let u_{ \varepsilon} denote a bounded sequence in L^p(Q) . Then the triple s_{ \varepsilon}: = (J^{\omega_0}_{ \varepsilon} u_{ \varepsilon}, \varepsilon,\|u_{\varepsilon}\|_{L^p(Q)}) defines a sequence in {\mathscr M}^{\omega_0} . Moreover, we have s_{ \varepsilon}\to s_0 in {\mathscr M} as \varepsilon\to0 if and only if s_0 = (J_0u,0,r) for some u\in {\mathscr{B}}^p , r\geq\|u\|_{ {\mathscr{B}}^p} , and u_{ \varepsilon}{\stackrel{2}{\rightharpoonup}_{{\omega_0}}}u .

    Proof.(i)This is a direct consequence of Lemma 3.6 (i) and the fact that the product of complete, separable metric spaces remains complete and separable.

    (ii)Let s_k: = (U_k, \varepsilon_k,r_k) denote a sequence in {\mathscr M}^{\omega_0} that converges in {\mathscr M} to some s_0 = (U_0, \varepsilon_0,r_0) . We need to show that s_0\in {\mathscr M}^{\omega_0} . By passing to a subsequence, it suffices to study the following three cases: \varepsilon_k>0 for all k\in \mathbb{N}_0 , \varepsilon_k = 0 for all k\in \mathbb{N}_0 , and \varepsilon_0 = 0 while \varepsilon_k>0 for all k\in \mathbb{N} .

    Case 1: \varepsilon_k>0 for all k\in \mathbb{N}_0 .

    W.l.o.g. we may assume that \inf_k \varepsilon_k>0 . Hence, there exist u_k\in L^p(Q) with U_k = J^{\omega_0}_{ \varepsilon_k}u_k . Since r_k\to r , we conclude that (u_k) is bounded in L^p(\Omega) . We thus may pass to a subsequence (not relabeled) such that u_k \rightharpoonup u_0 weakly in L^p(Q) , and

    \begin{equation} \|u_0\|_{L^p(Q)}\leq \liminf\limits_{k}\|u_k\|_{L^p(Q)}\leq \lim\limits_k r_k = r_0. \end{equation} (16)

    Moreover, U_k\to U in the metric space \mbox{Lin}( {\mathscr D}) implies pointwise convergence on {\mathscr D} , and thus for all \varphi_Q\in {\mathscr D}_Q we have U_k(\mathbf 1_{\Omega}\varphi_Q) = \int_Qu_k\varphi_Q\to \int_Qu_0\varphi_Q . We thus conclude that U_0(\mathbf 1_{\Omega}\varphi_Q) = \int_Q u_0\varphi_Q . Since {\mathscr D}_Q\subset L^q(Q) dense, we deduce that u_k \rightharpoonup u_0 weakly in L^p(Q) for the entire sequence. On the other hand the properties of the shift \tau imply that for any \varphi_\Omega\in {\mathscr D}_\Omega we have \varphi_\Omega(\tau_{\frac{\cdot}{ \varepsilon_k}}\omega_0)\to\varphi_\Omega(\tau_{\frac{\cdot}{ \varepsilon_0}}\omega_0) in L^q(Q) . Hence, for any \varphi_\Omega\in {\mathscr D}_\Omega and \varphi_Q\in {\mathscr D}_Q we have

    \begin{align*} U_k(\varphi_\Omega\varphi_Q)& = \int_Q u_k(x)\varphi_Q(x)\varphi_\Omega(\tau_{\frac{x}{ \varepsilon_k}}\omega_0)\,dx\\ &\to \int_Q u_0(x)\varphi_Q(x)\varphi_\Omega(\tau_{\frac{x}{ \varepsilon_0}}\omega_0)\,dx = J^{\omega_0}_{ \varepsilon_0}(\varphi_\Omega\varphi_Q) \end{align*}

    and thus (by linearity) U_0 = J^{\omega_0}_{ \varepsilon_0}u_0 .

    Case 2: \varepsilon_k = 0 for all k\in \mathbb{N}_0 .

    In this case there exist a bounded sequence u_k in {\mathscr{B}}^p with U_k = J_0u_k for k\in \mathbb{N} . By passing to a subsequence we may assume that u_k \rightharpoonup u_0 weakly in {\mathscr{B}}^p for some u_0\in {\mathscr{B}}^p with

    \begin{equation} \|u_0\|_{ {\mathscr{B}}^p}\leq \liminf\limits_{k}\|u_{ \varepsilon_k}\|_{ {\mathscr{B}}^p}\leq r_0. \end{equation} (17)

    This implies that U_k = J_0u_k\to J_0u_0 in \mbox{Lin}( {\mathscr D}) . Hence, U_0 = J_0u_0 and we conclude that s_0\in {\mathscr M}^{\omega_0} .

    Case 3: \varepsilon_k>0 for all k\in \mathbb{N} and \varepsilon_0 = 0 .

    There exists a bounded sequence u_k in L^p(Q) . Thanks to Lemma 3.5, by passing to a subsequence we may assume that u_k{\stackrel{2}{\rightharpoonup}_{{\omega_0}}} u_0 for some u\in {\mathscr{B}}^p with

    \begin{equation} \|u_0\|_{ {\mathscr{B}}^p}\leq \liminf\limits_{k}\|u_k\|_{L^p(Q)}\leq r_0. \end{equation} (18)

    Furthermore, Lemma 3.6 implies that J^{\omega_0}_{ \varepsilon_k}u_k\to J_0u_0 in \mbox{Lin}( {\mathscr D}) , and thus U_0 = J_0u_0 . We conclude that s_0\in {\mathscr M}^{\omega_0} .

    (iii)We first argue that the representation (14) of U by u is unique. In the case \varepsilon>0 suppose that u,v\in L^p(Q) satisfy U = J^{\omega_0}_{ \varepsilon}u = J^{\omega_0}_{ \varepsilon}v . Then for all \varphi_Q\in {\mathscr D}_Q we have \int_Q(u-v)\varphi_Q = J^{\omega_0}_ \varepsilon u(\mathbf 1_\Omega\varphi_Q)-J^{\omega_0}_ \varepsilon v(\mathbf 1_\Omega\varphi_Q) = U(\mathbf 1_\Omega\varphi_Q)-U(\mathbf 1_\Omega\varphi_Q) = 0 , and since {\mathscr D}_Q\subset L^q(Q) dense, we conclude that u = v . In the case \varepsilon = 0 the statement follows by a similar argument from the fact that {\mathscr D} is dense {\mathscr{B}}^q .

    To see (15) let u and U be related by (14). Since \overline {\mathscr D} (resp. {\mathscr D} ) is dense in L^q(Q) (resp. {\mathscr{B}}^q ), we have

    \begin{equation*} \begin{cases} \|u\|_{L^p(Q)} = \sup\limits_{\varphi\in\overline{ {\mathscr D}},\ \|\varphi\|_{ {\mathscr{B}}^q}\leq 1}|\int_Qu\varphi\,dx\,dP| = \sup\limits_{\varphi\in\overline{ {\mathscr D}},\ \|\varphi\|_{ {\mathscr{B}}^q}\leq 1}|U(\varphi)|&\text{if } \varepsilon > 0,\\ \|u\|_{ {\mathscr{B}}^p} = \sup\limits_{\varphi\in {\mathscr D},\ \|\varphi\|_{ {\mathscr{B}}^q}\leq 1}|\int_{\Omega}\int_{Q}u\varphi\,dx\,dP| = \sup\limits_{\varphi\in {\mathscr D},\ \|\varphi\|_{ {\mathscr{B}}^q}\leq 1}|U(\varphi)|&\text{if } \varepsilon = 0. \end{cases} \end{equation*}

    (iv)Let s_k = (U_k, \varepsilon_k,r_k) denote a sequence in {\mathscr M}^{\omega_0} that converges in {\mathscr M} to a limit s_0 = (U_0, \varepsilon_0,r_0) . By (ii) we have s_0\in {\mathscr M}^{\omega_0} . For k\in \mathbb{N}_0 let u_k in L^p(Q) or {\mathscr{B}}^p denote the representation of U_k in the sense of (14). We may pass to a subsequence such that one of the three cases in (ii) applies and (as in (ii)) either u_k weakly converges to u_0 (in L^p(Q) or {\mathscr{B}}^p ), or u_k{\stackrel{2}{\rightharpoonup}_{{\omega_0}}}u_0 . In any of these cases the claimed lower semicontinuity of \|\cdot\|_{\omega_0} follows from \varepsilon_k\to \varepsilon_0 , r_k\to r_0 , and (15) in connection with one of the lower semicontinuity estimates (16) – (18).

    (v)This follows from the definition and duality argument (15).

    (vi)Let s_k denote a sequence in {\mathscr M}^{\omega_0} . Let u_k in L^p(Q) or {\mathscr{B}}^p denote the (unique) representative of U_k in the sense of (14). Suppose that \|s_k\|_{\omega_0}\leq R . Then (r_k) and ( \varepsilon_k) are bounded sequences in \mathbb{R}_{\geq 0} , and \sup_{k}\|u_k\|\leq \sup_kr_k<\infty (where \|\cdot\| stands short for either \|\cdot\|_{L^p(Q)} or \|\cdot\|_{ {\mathscr{B}}^p} ). Thus we may pass to a subsequence such that r_k\to r_0 , \varepsilon_k\to \varepsilon_0 , and one of the following three cases applies:

    ● Case 1: \inf_{k\in \mathbb{N}_0} \varepsilon_k>0 . In that case we conclude (after passing to a further subsequence) that u_k \rightharpoonup u_0 weakly in L^p(Q) , and thus U_k\to U_0 = J^{\omega_0}_{ \varepsilon_0}u_0 in \mbox{Lin}( {\mathscr D}) .

    ● Case 2: \varepsilon_k = 0 for all k\in \mathbb{N}_0 . In that case we conclude (after passing to a further subsequence) that u_k \rightharpoonup u_0 weakly in {\mathscr{B}}^p(Q) , and thus U_k\to U_0 = J_0u_0 in \mbox{Lin}( {\mathscr D}) .

    ● Case 3: \varepsilon_k>0 for all k\in \mathbb{N} and \varepsilon_0 = 0 . In that case we conclude (after passing to a further subsequence) that u_k{\stackrel{2}{\rightharpoonup}_{{\omega_0}}}u_0 , and thus U_k\to U_0 = J_0u_0 in \mbox{Lin}( {\mathscr D}) .

    In all of these cases we deduce that s_0 = (U_0, \varepsilon_0,r_0)\in {\mathscr M}^{\omega_0} , and s_k\to s_0 in {\mathscr M} .

    (vii)This is a direct consequence of (ii) – (vi), and Lemma 3.6.

    Now we are in position to prove Theorem 3.12

    Proof of Theorem 3.12. Let {\mathscr M} , {\mathscr M}^{\omega_0} , J^{\omega_0}_ \varepsilon etc. be defined as in Lemma 3.17.

    Step 1. (Identification of (u_{ \varepsilon}) with a tight {\mathscr M} -valued sequence). Since u_{ \varepsilon}\in {\mathscr{B}}^p , by Fubini's theorem, we have u_{ \varepsilon}(\omega,\cdot)\in L^p(Q) for P -a.a. \omega\in\Omega . By modifying u_{ \varepsilon} on a null-set in \Omega\times Q (which does not alter two-scale limits in the mean), we may assume w.l.o.g. that u_{ \varepsilon}(\omega,\cdot)\in\ L^p(Q) for all \omega\in\Omega . Consider the measurable function s_{ \varepsilon}:\Omega\to {\mathscr M} defined as

    \begin{equation*} s_{ \varepsilon}(\omega): = \begin{cases} \big(J^{\omega}_{ \varepsilon} u_{ \varepsilon}(\omega,\cdot), \varepsilon,\|u_{ \varepsilon}(\omega,\cdot)\|_{L^p(Q)}\big)&\text{if }\omega\in\Omega_0\\ (0,0,0)&\text{else.} \end{cases} \end{equation*}

    We claim that (s_{ \varepsilon}) is tight. To that end consider the integrand h:\Omega\times {\mathscr M}\to(-\infty,+\infty] defined by

    \begin{equation*} h(\omega,(U, \varepsilon,r)): = \begin{cases} \|(U, \varepsilon,r)\|_{\omega}^p&\text{if }\omega\in\Omega_0\text{ and }(U, \varepsilon,r)\in {\mathscr M}^{\omega},\\ +\infty&\text{else.} \end{cases} \end{equation*}

    From Lemma 3.17 (iv) and (vi) we deduce that h is a normal integrand and h(\omega,\cdot) has compact sublevels for all \omega\in\Omega . Moreover, for all \omega_0\in\Omega_0 we have s_{ \varepsilon}(\omega_0)\in {\mathscr M}^{\omega_0} and thus h(\omega_0,s_{ \varepsilon}(\omega_0)) = 2\|u_{ \varepsilon}(\omega_0,\cdot)\|^p_{L^p(Q)}+ \varepsilon . Hence,

    \begin{equation*} \int_\Omega h(\omega,s_{ \varepsilon}(\omega))\,dP(\omega) = 2\|u_{ \varepsilon}\|^p_{ {\mathscr{B}}^p}+ \varepsilon. \end{equation*}

    We conclude that (s_{ \varepsilon}) is tight.

    Step 2. (Compactness and definition of \boldsymbol{\nu} ). By appealing to Theorem 3.16 there exists a subsequence (still denoted by \varepsilon ) and a Young measure \boldsymbol{\mu} that is generated by (s_{ \varepsilon}) . Let \boldsymbol{\mu_1} denote the first component of \boldsymbol{\mu} , i.e. the Young measure on \mbox{Lin}( {\mathscr D}) characterized for \omega\in\Omega by

    \begin{equation*} \int_{\mbox{Lin}( {\mathscr D})}f(\xi)\,d\mu_{1,\omega}(\xi) = \int_{ {\mathscr M}}f(\xi_1)\,d\mu_\omega(\xi), \end{equation*}

    for all f:\mbox{Lin}( {\mathscr D})\to \mathbb{R} continuous and bounded, where {\mathscr M}\ni\xi = (\xi_1,\xi_2,\xi_3)\to\xi_1\in\mbox{Lin}( {\mathscr D}) denotes the projection to the first component. By Balder's theorem, \mu_\omega is concentrated on the limit points of (s_{ \varepsilon}(\omega)) . By Lemma 3.17 we deduce that for all \omega\in\Omega_0 any limit point s_0(\omega) of s_{ \varepsilon}(\omega) has the form s_0(\omega) = (J_0u,0,r) where 0\leq r<\infty and u\in {\mathscr{B}}^p is a \omega -two-scale limit of a subsequence of u_{ \varepsilon}(\omega,\cdot) . Thus, \mu_{1,\omega} is supported on \{J_0u\,:\,u\in {\mathscr{C\!P}}(\omega,(u_{ \varepsilon}(\omega,\cdot))\} which in particular is a subset of ( {\mathscr{B}}^q)^* . Since J_0: {\mathscr{B}}^p\to ( {\mathscr{B}}^q)^* is an isometric isomorphism (by the Riesz-Frechét theorem), we conclude that \boldsymbol{\nu} = \{\nu_\omega\}_{\omega\in\Omega} , \nu_\omega(B): = \mu_{1,\omega}(J_0B) (for all Borel sets B\subset {\mathscr{B}}^p where {\mathscr{B}}^p is equipped with the weak topology) defines a Young measure on {\mathscr{B}}^p , and for all \omega\in\Omega_0 , \nu_\omega is supported on {\mathscr{C\!P}}(\omega,(u_{ \varepsilon}(\omega,\cdot)) .

    Step 3. (Lower semicontinuity estimate). Note that h:\Omega\times {\mathscr M}\to[0,+\infty] ,

    \begin{equation*} h(\omega,(U, \varepsilon,r)): = \begin{cases} \sup_{\varphi\in\overline {\mathscr D},\,\|\varphi\|_{ {\mathscr{B}}^q}\leq 1}|U(\varphi)|^p&\text{if }\omega\in\Omega_0\text{ and }(U, \varepsilon,r)\in {\mathscr M}^{\omega}, \varepsilon > 0,\\ \sup_{\varphi\in {\mathscr D},\,\|\varphi\|_{ {\mathscr{B}}^q}\leq 1}|U(\varphi)|^p&\text{if }\omega\in\Omega_0\text{ and }(U, \varepsilon,r)\in {\mathscr M}^{\omega}, \varepsilon = 0,\\ +\infty&\text{else.} \end{cases} \end{equation*}

    defines a normal integrand, as can be seen as in the proof of Lemma 3.17. Thus Theorem 3.16 implies that

    \begin{equation*} \liminf\limits_{ \varepsilon\to 0}\int_\Omega h(\omega,s_{ \varepsilon}(\omega))\,dP(\omega)\geq \int_\Omega\int_{ {\mathscr M}}h(\omega,\xi)\,d\mu_\omega(\xi)dP(\omega). \end{equation*}

    In view of Lemma 3.17 we have \sup_{\varphi\in\overline {\mathscr D},\,\|\varphi\|_{ {\mathscr{B}}^q}\leq 1}|(J^\omega_ \varepsilon u_ \varepsilon)(\omega,\cdot))(\varphi)| = \|u_{ \varepsilon}(\omega,\cdot)\|_{L^p(Q)} for \omega\in\Omega_0 , and thus the left-hand side turns into \liminf_{ \varepsilon\to 0}\|u_{ \varepsilon}\|^p_{ {\mathscr{B}}^p} . Thanks to the definition of \boldsymbol{\nu} the right-hand side turns into \int_\Omega \int_{ {\mathscr{B}}^p}\|v\|_{ {\mathscr{B}}^p}^p\,d\nu_\omega(v)dP(\omega) .

    Step 4. (Identification of the two-scale limit in the mean). Let \varphi\in {\mathscr D}_0 . Then h:\Omega\times {\mathscr M}\to[0,+\infty] ,

    \begin{equation*} h(\omega,(U, \varepsilon,r)): = \begin{cases} U(\varphi)&\text{if }\omega\in\Omega_0,\,(U, \varepsilon,r)\in {\mathscr M}^\omega,\\ +\infty&\text{else.} \end{cases} \end{equation*}

    defines a normal integrand. Since h(\omega,s_{ \varepsilon}(\omega)) = \int_Qu_{ \varepsilon}(\omega,x) \mathcal{T}_{\varepsilon}^*\varphi(\omega,x)\,dx for P -a.a. \omega\in\Omega , we deduce that |h(\cdot, s_{ \varepsilon}(\cdot))| is uniformly integrable. Thus, (13) applied to \pm h and the definition of \boldsymbol{\nu} imply that

    \begin{eqnarray*} \lim\limits_{ \varepsilon\to 0}\int_\Omega\int_Qu_{ \varepsilon}(\omega,x)( \mathcal{T}_{\varepsilon}^*\varphi)(\omega,x)\,dx\,dP(\omega)& = & \lim\limits_{ \varepsilon\to 0}\int_\Omega h(\omega,s_{ \varepsilon}(\omega))\,dP(\omega)\\ & = &\int_\Omega\int_{ {\mathscr{B}}^p}h(\omega,v)\,d\nu_\omega(v)\,dP(\omega)\\ & = &\int_\Omega\int_{ {\mathscr{B}}^p}\langle {{\int_Qv\varphi}} \rangle\,d\nu_\omega(v)\,dP(\omega). \end{eqnarray*}

    Set u: = \int_\Omega\int_{ {\mathscr{B}}^p}v\,d\nu_\omega(v)dP(\omega)\in {\mathscr{B}}^p . Then Fubini's theorem yields

    \begin{eqnarray*} \lim\limits_{ \varepsilon\to 0}\int_\Omega\int_Qu_{ \varepsilon}(\omega,x)( \mathcal{T}_{\varepsilon}^*\varphi)(\omega,x)\,dx\,dP(\omega)& = & \langle {{\int_Qu\varphi}} \rangle. \end{eqnarray*}

    Since \mbox{span}( {\mathscr D}_0)\subset {\mathscr{B}}^q dense, we conclude that u_{ \varepsilon} {\stackrel{2}{\rightharpoonup}} u .

    Step 5. Recovery of quenched two-scale convergence. Suppose that \nu_\omega is a delta distribution on {\mathscr{B}}^p , say \nu_\omega = \delta_{\tilde u(\omega)} for some measurable \tilde u:\Omega\to {\mathscr{B}}^p . Note that h:\Omega\times {\mathscr M}\to[0,+\infty] ,

    \begin{equation*} h(\omega,(U, \varepsilon,r)): = -d(U,J_0\tilde u(\omega);\mbox{Lin}( {\mathscr D})) \end{equation*}

    is a normal integrand and |h(\cdot, s_{ \varepsilon}(\cdot))| is uniformly integrable. Thus, (13) yields

    \begin{eqnarray*} &&\limsup\limits_{ \varepsilon\to 0}\int_{\Omega} d(J^\omega_{ \varepsilon} u_{ \varepsilon}(\omega,\cdot),J_0\tilde u(\omega);\mbox{Lin}( {\mathscr D}))\,dP(\omega)\\ & = &-\liminf\limits_{ \varepsilon\to 0}\int_\Omega h(\omega,s_{ \varepsilon}(\omega))\,dP(\omega)\\ &\leq&-\int_\Omega\int_{ {\mathscr{B}}^p}h(\omega,J_0v)\,d\nu_\omega(v)\,dP(\omega) = -\int_\Omega h(\omega,J_0\tilde u(\omega))\,dP(\omega) = 0. \end{eqnarray*}

    Thus, there exists a subsequence (not relabeled) such that d(J^\omega_{ \varepsilon} u_{ \varepsilon}(\omega,\cdot),J_0\tilde u(\omega);\mbox{Lin}( {\mathscr D}))\to 0 for a.a. \omega\in\Omega_0 . In view of Lemma 3.6 this implies that u_{ \varepsilon}{\stackrel{2}{\rightharpoonup}_{{\omega}}}\tilde u(\omega) for a.a. \omega\in\Omega_0 .

    Proof of Lemma 3.14. Step 1. Representation of the functional by a lower semicontinuous integrand on {\mathscr M} .

    For all \omega_0\in\Omega_0 and s = (U, \varepsilon,r)\in {\mathscr M}^{\omega_0} we write \pi^{\omega_0}(s) for the unique representation u in {\mathscr{B}}^p (resp. L^p(Q) ) of U in the sense of (14). We thus may define for \omega\in\Omega_0 and s\in {\mathscr M}^{\omega_0} the integrand

    \begin{equation*} \overline h(\omega_0,s): = \begin{cases} \int_Qh(\tau_{\frac{x}{ \varepsilon}}\omega,x,(\pi^{\omega_0}s)(x))\,dx&\text{if }s = (U, \varepsilon,s)\text{ with } \varepsilon > 0,\\ \int_\Omega\int_Qh(\omega,x,(\pi^{\omega_0}s)(x))\,dx\,dP(\omega)&\text{if }s = (U, \varepsilon,s)\text{ with } \varepsilon = 0. \end{cases} \end{equation*}

    We extend \overline h(\omega_0,\cdot) to {\mathscr M} by +\infty , and define \overline h(\omega,\cdot)\equiv 0 for \omega\in\Omega\setminus\Omega_0 . We claim that \overline h(\omega,\cdot): {\mathscr M}\to(-\infty,+\infty] is lower semicontinuous for all \omega\in\Omega . It suffices to consider \omega_0\in\Omega_0 and a convergent sequence s_k = (U_k, \varepsilon_k,r_k) in {\mathscr M}^{\omega_0} . For brevity we only consider the (interesting) case when \varepsilon_k\downarrow \varepsilon_0 = 0 . Set u_k: = \pi^{\omega_0}(s_k) . By construction we have

    \begin{equation*} \overline h(\omega_0,s_k) = \int_Q h(\tau_{\frac{x}{ \varepsilon_k}}\omega_0,u_k(\omega_0,x))\,dx, \end{equation*}

    and

    \begin{equation*} \overline h(\omega_0,s_0) = \int_{\Omega}\int_Q h(\omega,x,u_0(\omega,x))\,dx\,dP(\omega). \end{equation*}

    Since s_k\to s_0 and \varepsilon_k\to 0 , Lemma 3.17 (vi) implies that u_k{\stackrel{2}{\rightharpoonup}_{{\omega_0}}}u_0 , and since h satisfies 12 from Remark 3, we conclude that \liminf\limits_{k} \overline h(\omega_0,s_k)\geq \overline h(\omega_0,s_0) , and thus \overline h is a normal integrand.

    Step 2. Conclusion. As in Step 1 of the proof of Theorem 3.12 we may associate with the sequence (u_{ \varepsilon}) a sequence of measurable functions s_{ \varepsilon}:\Omega\to {\mathscr M} that (after passing to a subsequence that we do not relabel) generates a Young measure \boldsymbol{\mu} on {\mathscr M} . Since by assumption u_{ \varepsilon} generates the Young measure {\boldsymbol \nu} on {\mathscr{B}}^p , we deduce that the first component \boldsymbol{\mu_1} satisfies \nu_\omega(B) = \mu_{\omega}(J_0B) for any Borel set B . Applying (13) to the integrand \overline h of Step 1, yields

    \begin{align*} \liminf\limits_{ \varepsilon\to 0}\int_\Omega\int_Q &h(\tau_{\frac{x}{ \varepsilon}}\omega_0,u_{ \varepsilon}(\omega_0,x))\,dx\,dP(\omega)\\ & = \liminf\limits_{ \varepsilon\to 0}\int_\Omega\overline h(\omega,s_{ \varepsilon}(\omega))\,dP(\omega)\\ &\geq\int_\Omega\int_{ {\mathscr M}}\overline h(\omega,\xi)\,d\mu_\omega(\xi)\,dP(\omega)\\ & = \int_\Omega\int_{ {\mathscr{B}}^p}\Big(\int_\Omega\int_Qh( \tilde{\omega},x,v( \tilde{\omega},x))\,dx\,dP( \tilde{\omega})\Big)\,d\nu_\omega(v)\,dP(\omega). \end{align*}

    Proof of Lemma 3.13. By (ii) and (iii) the sequence (\tilde u_{ \varepsilon}) is bounded in {\mathscr{B}}^p and thus we can pass to a subsequence such that (\tilde u_{ \varepsilon}) generates a Young measure \boldsymbol \nu . Set \tilde u: = \int_\Omega\int_{ {\mathscr{B}}^p}v\,d\nu_\omega(v)\,dP(\omega) and note that Theorem 3.12 implies that \tilde u_{ \varepsilon} \overset{2}{\rightharpoonup} \tilde u weakly two-scale in the mean. On the other hand the theorem implies that \nu_\omega concentrates on the quenched two-scale cluster points of (u^\omega_{ \varepsilon}) (for a.a. \omega\in\Omega ). Hence, in view of (i) we conclude that for a.a. \omega\in\Omega the measure \nu_\omega is a Dirac measure concentrated on u , and thus \tilde u = u a.e. in \Omega\times Q .

    In this section we revisit a standard model example of stochastic homogenization of integral functionals from the viewpoint of stochastic two-scale convergence and unfolding. In particular, we discuss two examples of convex homogenization problems that can be treated with stochastic two-scale convergence in the mean, but not with the quenched variant. In the first example in Section 4.1 the randomness is nonergodic and thus quenched two-scale convergence does not apply. In the second example, in Section 4.2, we consider a variance-regularization to treat a convex minimization problem with degenerate growth conditions. In these two examples we also demonstrate the simplicity of using the stochastic unfolding operator. Furthermore, in Section 4.3 we use the results of Section 3.3 to further reveal the structure of the previously obtained limits in the classical ergodic case with non-degenerate growth with help of Young measures. In particular, we show how to lift mean homogenization results to quenched statements.

    In this section we consider a nonergodic stationary medium. Such random ensembles arise naturally, e.g., in the context of periodic representative volume element (RVE) approximations, see [13]. For example, we may consider a family of i.i.d. random variables \left\lbrace {{\overline{\omega}(z)}} \right\rbrace_{z\in \mathbb{Z}^d} . A realization of a stationary and ergodic random checkerboard is given by

    \begin{equation*} \omega: \mathbb{R}^d \to \mathbb{R}, \quad \omega(x) = \sum\limits_{i\in \mathbb{Z}^d} \mathbf{1}_{i+y+\Box}(x)\overline{\omega}(\lfloor x \rfloor), \end{equation*}

    where \lfloor x\rfloor\in \mathbb{Z}^d is the integer part of x and y\in \Box is the center of the checkerboard chosen uniformly from \Box = [0,1)^d . For L\in \mathbb{N} , we may consider the map \pi_{L}: \omega \mapsto \omega_{L} given by \pi_{L}\omega(x) = \omega(x) for x\in [0,L)^d and \pi_{L}\omega is L -periodically extended. The push forward of the map \pi_{L} defines a stationary and nonergodic probability measure, that is a starting point in the periodic RVE method. Another standard example of a nonergodic structure may be obtained by considering a medium with a noncompatible quasiperiodic microstructure, see [38,Example 1.2].

    In this section we consider the following situation. Let p\in (1,\infty) and Q\subset \mathbb{R}^d be open and bounded. We consider V:\Omega\times Q\times \mathbb{R}^{{d}}\rightarrow \mathbb{R}^{{}} and assume:

    (A1) V(\cdot,\cdot, F) is \mathcal{F}\otimes \mathcal{L}(Q) -measurable for all F\in \mathbb{R}^{d} .

    (A2) V(\omega, x, \cdot) is convex for a.a. (\omega,x)\in \Omega\times Q .

    (A3)There exists C>0 such that

    \begin{equation*} \frac{1}{C}|F|^p-C\leq V(\omega, x, F) \leq C(|F|^p+1) \end{equation*}

    for a.a. (\omega,x) \in \Omega\times Q and all F\in \mathbb{R}^{{d}} .

    We consider the functional

    \begin{equation} \mathcal{E}_{\varepsilon}: L^p(\Omega)\otimes W^{1,p}_0(Q)\rightarrow \mathbb{R}^{{}},\quad \mathcal{E}_{\varepsilon}(u) = \langle {{\int_Q V(\tau_{\frac{x}{\varepsilon}}\omega, x,\nabla u(\omega,x))dx}} \rangle. \end{equation} (19)

    Under assumptions (A1)-(A3), in the limit \varepsilon\rightarrow 0 we obtain the two-scale functional

    \begin{align} \begin{split} & \mathcal{E}_0:\left({{L^p_{{ {\mathrm{inv}}}}(\Omega)\otimes W^{1,p}_0(Q)}}\right) \times \left({{L^p_{ {\mathrm{pot}}}(\Omega)\otimes L^p(Q)}}\right) \to \mathbb{R},\\ & \mathcal{E}_0(u,\chi) = \langle {{\int_Q V(\omega, x, \nabla u(\omega,x)+ \chi(\omega,x)) dx}} \rangle. \end{split} \end{align} (20)

    Theorem 4.1 (Two-scale homogenization). Let p\in (1,\infty) and Q\subset \mathbb{R}^d be open and bounded. Assume (A1)-(A3).

    (i)(Compactness and liminf inequality.) Let u_{ \varepsilon} \in L^p(\Omega)\otimes W^{1,p}_0(Q) be such that \limsup_{\varepsilon\rightarrow 0}\mathcal{E}_{\varepsilon}(u_{ \varepsilon})<\infty .There exist (u,\chi) \in \left({{L^p_{{ {\mathrm{inv}}}}(\Omega)\otimes W^{1,p}_0(Q)}}\right) \times \left({{L^p_{ {\mathrm{pot}}}(\Omega)\otimes L^p(Q)}}\right) and a subsequence (not relabeled) such that

    \begin{align} & u_{ \varepsilon} \overset{2}{\rightharpoonup} u \;{{in}}\;L^p(\Omega \times Q), \quad \nabla u_{ \varepsilon} \overset{2}{\rightharpoonup} \nabla u+\chi \;{{in}}\; L^p(\Omega \times Q), \end{align} (21)
    \begin{align} & \liminf\limits_{\varepsilon\to 0}\mathcal{E}_{ \varepsilon}(u_{ \varepsilon})\geq \mathcal{E}_0(u,\chi). \end{align} (22)

    (ii)(Limsup inequality.) Let (u,\chi)\in \left({{L^p_{{ {\mathrm{inv}}}}(\Omega)\otimes W^{1,p}_0(Q)}}\right) \times \left({{L^p_{ {\mathrm{pot}}}(\Omega)\otimes L^p(Q)}}\right) . There exists a sequence u_{ \varepsilon} \in L^p(\Omega)\otimes W^{1,p}_0(Q) such that

    \begin{align} & u_{ \varepsilon} \overset{2}{\rightarrow } u \;{{in}}\; L^p(\Omega \times Q), \quad \nabla u_{ \varepsilon} \overset{2}{\rightarrow } \nabla u+\chi \;{{in}}\; L^p(\Omega \times Q), \end{align} (23)
    \begin{align} & \limsup\limits_{\varepsilon\rightarrow 0}\mathcal{E}_{\varepsilon}(u_{ \varepsilon}) \leq \mathcal{E}_0(u,\chi). \end{align} (24)

    Proof of Theorem 4.1. (i) The Poincaré inequality and (A3) imply that u_{ \varepsilon} is bounded in L^p(\Omega)\otimes W^{1,p}(Q) . By Proposition 1 (ii) there exist u \in L^p_{{ {\mathrm{inv}}}}(\Omega)\otimes W^{1,p}(Q) and \chi \in L^p_{ {\mathrm{pot}}}(\Omega)\otimes L^p(Q) such that (21) holds. Also, note that P_{\mathrm{inv}}u_{ \varepsilon} \rightharpoonup u weakly in L^p(\Omega)\otimes W^{1,p}(Q) and P_{\mathrm{inv}}u_{ \varepsilon} \in L^p_{\mathrm{inv}}(\Omega)\otimes W^{1,p}_0(Q) , which implies that u also has 0 boundary values, i.e., u \in L^p_{\mathrm{inv}}(\Omega)\otimes W^{1,p}_0(Q) . Finally, we note that, see [19,Proposition 3.5 (i)],

    \begin{equation} \langle {{\int_{Q} V(\tau_{\frac{x}{\varepsilon}}\omega,x,v(\omega,x))}} \rangle = \langle {{\int_{Q} V(\omega,x, \mathcal{T}_{\varepsilon} v(\omega,x))}} \rangle \quad \text{for any }v \in L^p(\Omega\times Q), \end{equation} (25)

    and thus using the convexity of V we conclude

    \begin{equation*} \liminf\limits_{\varepsilon\to 0} \mathcal{E}_{ \varepsilon}(u_{ \varepsilon}) = \liminf\limits_{\varepsilon\to 0} \langle {{\int_{Q}V(\omega,x, \mathcal{T}_{\varepsilon} \nabla u_{ \varepsilon})}} \rangle\geq \mathcal{E}_0(u,\chi). \end{equation*}

    (ii) The existence of a sequence u_{ \varepsilon} with (23) follows from Proposition 1 (iii). Furthermore, (25) and the growth assumption (A3) yield

    \begin{equation*} \lim\limits_{\varepsilon\to 0}\mathcal{E}_{\varepsilon}(u_{ \varepsilon}) = \lim\limits_{\varepsilon\to 0} \langle {{\int_{Q}V(\omega,x, \mathcal{T}_{\varepsilon} \nabla u_{ \varepsilon})}} \rangle = \mathcal{E}_{0}(u,\chi). \end{equation*}

    This concludes the claim, in particular, we even show a stronger result stating convergence of the energy.

    Remark 4 (Convergence of minimizers). We consider the setting of Theorem 4.1. Let u_{ \varepsilon}\in L^p(\Omega)\otimes W^{1,p}_0(Q) be a minimizer of the functional

    \begin{equation*} \mathcal{I}_{\varepsilon}:L^p(\Omega)\otimes W^{1,p}_0(Q)\to \mathbb{R}, \quad \mathcal{I}_{\varepsilon}(u) = {\mathcal{E}}_{ \varepsilon}(u) - \langle {{\int_{Q}u_{ \varepsilon} f_{ \varepsilon} dx}} \rangle, \end{equation*}

    where f_{ \varepsilon} \in L^q(\Omega\times Q) and f_{ \varepsilon} \overset{2}{\to}f with f\in L^q(Q) and \frac{1}{p}+\frac{1}{q} = 1 . By a standard argument from the theory of \Gamma -convergence Theorem 4.1 (cf. [34,Corollary 7.2]) implies that there exist a subsequence (not relabeled), u\in L^p_{ {\mathrm{inv}}}(\Omega)\times W^{1,p}_0(Q) , and \chi\in L^p_{ {\mathrm{pot}}}(\Omega)\otimes L^p(Q) such that u_{ \varepsilon} \overset{2}{\rightharpoonup} u \text{ in } L^p(\Omega \times Q) , \nabla u_{ \varepsilon} \overset{2}{\rightharpoonup} \nabla u+\chi \text{ in } L^p(\Omega \times Q) , and

    \begin{equation*} \lim\limits_{ \varepsilon\to 0}\min\mathcal{I}_{ \varepsilon} = \lim\limits_{ \varepsilon\to 0}\mathcal{I}_{ \varepsilon}(u_{ \varepsilon}) = \mathcal{I}_0(u,\chi) = \min\mathcal{I}_0, \end{equation*}

    where \mathcal{I}_0: L^p_{\mathrm{inv}}(\Omega)\otimes W^{1,p}_0(Q)\to \mathbb{R} is given by \mathcal{I}_{0}(u) = \mathcal{E}_0(u)-\int_{Q}f u dx . This, in particular, rigorously justifies the formal two-scale expansion \nabla u_{ \varepsilon}(x) \approx \nabla u_{0}(\omega, x) + \chi( \tau_{\frac{x}{\varepsilon}}\omega, x) .

    Remark 5 (Uniqueness). If V(\omega, x,\cdot) is strictly convex the minimizers are unique and the convergence in the above remark holds for the entire sequence.

    In this section we consider homogenization of convex functionals with degenerate growth. More precisely, we consider an integrand V that satisfies (A1), (A2) and the following assumption (as a replacement of (A3)):

    (A3')There exists C>0 and a random variable \lambda\in L^1(\Omega) such that

    \begin{equation} \langle {{\lambda^{-\frac{1}{p-1}}}} \rangle^{p-1} < C \end{equation} (26)

    and

    \begin{equation*} \lambda(\omega)|F|^p-C\leq V(\omega, x, F) \leq C(\lambda(\omega)|F|^p+1) \end{equation*}

    for a.a. (\omega,x) \in \Omega\times Q and all F\in \mathbb{R}^{{d}} .

    Moreover, we assume that \langle {{\cdot}} \rangle is ergodic. For \varepsilon>0 we consider the following functional

    \begin{equation*} \mathcal{E}_{\varepsilon}:L^1(\Omega \times Q) \to \mathbb{R}\cup\left\lbrace {{\infty}} \right\rbrace, \qquad \mathcal{E}_{\varepsilon}(u) = \langle {{\int_{Q} V(\tau_{\frac{x}{\varepsilon}}\omega, x, \nabla u) dx}} \rangle, \end{equation*}

    for u \in X_{\varepsilon} and \mathcal{E}_{\varepsilon}(u) = \infty otherwise. Here X_{\varepsilon} denotes the closure of \left\lbrace {{u\in L^p(\Omega)\otimes W^{1,p}_0(Q)}} \right\rbrace w.r.t. the weighted norm

    \begin{equation*} \|u\|_{\lambda_{\varepsilon}} : = \langle {{\int_{Q}\lambda(\tau_{\frac{x}{\varepsilon}}\omega) |\nabla u|^p dx}} \rangle^{\frac{1}{p}}. \end{equation*}

    Recently, in [29,20,21] it shown that \mathcal{E}_\varepsilon Mosco-converges to the functional

    \begin{equation*} \mathcal E_{\hom}: L^1(Q)\to \mathbb{R}\cup \left\lbrace {{\infty}} \right\rbrace,\qquad \mathcal E_{\hom}(u): = \int_Q V_{\hom}(x,\nabla u(x))\,dx, \end{equation*}

    for u \in W^{1,p}_0(Q) and \mathcal E_{\hom}(u) = \infty otherwise, where V_{\hom}:Q\times \mathbb{R}^{d}\to \mathbb{R} is given by the homogenization formula,

    \begin{align} V_{\hom}(x,F) = \inf\limits_{\chi\in L^p_{ {\mathrm{pot}}}(\Omega)}\langle {{V(\omega, x,F+\chi(\omega))}} \rangle, \end{align} (27)

    for x\in Q and F\in \mathbb{R}^{d} . Moreover, it is shown that V_{\mathrm{hom}} is a normal convex integrand that satisfies a standard p -growth condition. Note that the assumption (A3') in comparison to (A3) makes a genuine difference in regard to the homogenization formula (27). In particular, in the setting of assumption (A3) minimizers are attained due to the coercivity of the underlying functional in L^p_{ {\mathrm{pot}}}(\Omega) . It is thus easy to see that the homogenized integrand satisfies p -growth condition as well, see Section 4.3 below. On the other hand, in the setting of this section assuming (A3'), (27) is a degenerate minimization problem and a priori minimizers will only have finite first moments. An additional argument is required to infer that V_{\hom} in (27) is non-degenerate, in particular, in [29,Theorem 3.1] it is shown that there exists a constant C'>0 such that for all x\in Q and F\in \mathbb{R}^d we have

    \begin{equation} \frac{1}{C'} |F|^p - C' \leq V_{\mathrm{hom}}(x,F) \leq C'\left({{|F|^p+1}}\right). \end{equation} (28)

    One of the difficulties in the proof of the homogenization result for \mathcal E_{\varepsilon} is due to the fact that the domain of the functionals are \varepsilon -dependent. Moreover, assumption (A3') only yields equicoercivity in L^1(\Omega)\otimes W^{1,1}_0(Q) , while the limit \mathcal E_{\hom} is properly defined on W^{1,p}(Q) . Therefore, in practice it is convenient to regularize the problem: For \delta>0 we consider the regularized homogenization formula

    \begin{equation*} V_{\hom,\delta}(x,F) = \inf\limits_{\chi\in L^p_{ {\mathrm{pot}}}(\Omega)}\langle {{V(\omega, x,F+\chi(\omega))+\delta|\chi(\omega)|^p}} \rangle. \end{equation*}

    It is simple to show that the infimum on the right-hand side is attained by a unique minimizer. We also consider the corresponding regularized homogenized integral functional

    \begin{equation*} \mathcal E_{\hom,\delta}:L^1(Q)\to \mathbb{R}\cup \left\lbrace {{\infty}} \right\rbrace,\qquad \mathcal E_{\hom,\delta}(u): = \int_QV_{\hom,\delta}(\nabla u)\,dx, \end{equation*}

    for u\in W^{1,p}_0(Q) and \mathcal{E}_{\hom, \delta}(u) = \infty otherwise. Furthermore, thanks to (A3'), it is relatively easy to see that this regularization is consistent:

    Lemma 4.2. Let p\in (1,\infty) and Q\subset \mathbb{R}^d be open and bounded. Assume (A1), (A2) and (A3'). Then, for all x\in Q and F\in \mathbb{R}^{d} , we have

    \begin{equation} \lim\limits_{\delta\to 0} V_{\hom,\delta}(x,F) = V_{\hom}(x,F). \end{equation} (29)

    Moreover, \mathcal{E}_{\hom, \delta} Mosco converges to \mathcal{E}_{\hom} as \delta \to 0 , i.e., the following statements hold:

    (i)If u_{\delta}\rightharpoonup u weakly in L^1(Q) , then

    \begin{equation*} \liminf\limits_{\delta \to 0}\mathcal{E}_{\mathrm{hom},\delta}(u_{\delta}) \geq \mathcal{E}_{\mathrm{hom}}(u). \end{equation*}

    (ii)For any u \in L^1(Q) there exists a sequence u_{\delta}\in L^1(Q) such that

    \begin{equation*} u_{\delta} \to u \quad \;\mathit{\text{strongly in}}\;L^1(Q), \quad \mathcal{E}_{\mathrm{hom},\delta}(u_{\delta})\to \mathcal{E}_{\mathrm{hom}}(u). \end{equation*}

    Proof. Let F\in \mathbb{R}^{d} and x\in Q . Since \delta>0 , we have V_{\hom,\delta}(x,F)\geq V_{\hom}(x,F) . On the other hand, we consider a minimizing sequence \chi_{\eta}\in L^p_{\mathrm{pot}}(\Omega) in (27), e.g.,

    \begin{equation*} \langle {{V(\omega, x, F+ \chi_{\eta})}} \rangle \leq V_{\hom}(x,F) + \eta. \end{equation*}

    We have

    \begin{equation*} V_{\hom,\delta}(x,F)\leq \langle {{V(\omega, x, F+ \chi_{\eta})+ \delta |\chi_{\eta}|^p}} \rangle \leq V_{\hom}(x,F)+ \eta + \delta \langle {{|\chi_{\eta}|^p}} \rangle. \end{equation*}

    Letting first \delta\to 0 and then \eta\to 0 , we conclude (29).

    We further consider a sequence u_{\delta} such that u_{\delta}\rightharpoonup u weakly in L^1(Q) as \delta \to 0 . We assume without loss of generality that \limsup_{\delta \to 0}\mathcal{E}_{\hom,\delta}(u_{\delta})< \infty . This, in particular, with the help of (28) and the Poincaré inequality implies that \limsup_{\delta \to 0} \|u_{\delta}\|_{W^{1,p}_0(Q)}< \infty . Thus, up to a subsequence, we have u_{\delta} \rightharpoonup u weakly in W^{1,p}_0(Q) . Using this, we obtain

    \begin{equation*} \liminf\limits_{\delta \to 0}\mathcal{E}_{\hom, \delta}(u_{\delta})\geq \liminf\limits_{\delta \to 0}\mathcal{E}_{\hom}(u_{\delta}) \geq \mathcal{E}_{\hom}(u). \end{equation*}

    The first inequality follows by (29) and the second is a consequence of the fact that V_{\hom}(x,\cdot) is convex and of Fatou's Lemma. We conclude that (i) holds.

    If u \notin \mathrm{dom}(\mathcal{E}_{\hom}) , we simply choose u_{\delta} = u . On the other hand, for u\in \mathrm{dom}(\mathcal{E}_{\hom}) = W^{1,p}_0(Q) , (29) and the dominated convergence theorem yield

    \begin{equation*} \lim\limits_{\delta\to 0}\mathcal{E}_{\hom, \delta}(u) = \mathcal{E}_{\hom}(u). \end{equation*}

    This means that (ii) holds.

    In the following we introduce a variance regularization of the original functional \mathcal E_{ \varepsilon} that removes the degeneracy of the problem and thus can be analyzed by the standard strategy of Section 4.1. For \delta>0 , we consider \mathcal{E}_{\varepsilon,\delta}:L^1(\Omega \times Q)\rightarrow \mathbb{R}^{{}} ,

    \begin{equation} \mathcal{E}_{\varepsilon,\delta}(u) = \langle {{\int_Q V(\tau_{\frac{x}{\varepsilon}}\omega, x,\nabla u(x))+\delta|\nabla u(x)-\langle {{\nabla u( x)}} \rangle|^pdx}} \rangle, \end{equation} (30)

    for u\in L^p(\Omega)\otimes W^{1,p}_0(Q) and \mathcal{E}_{\varepsilon, \delta} = \infty otherwise. Due to the structure of the additional term, we call it a variance-regularization and we note that it only becomes active for non-deterministic functions. For fixed \delta>0 , the functional \mathcal E_{\varepsilon,\delta} is equicoercive on L^p(\Omega)\otimes W^{1,p}_0(Q) :

    Lemma 4.3. Let p\in (1,\infty) and Q\subset \mathbb{R}^d be open and bounded. Assume (A1) and (A3'). Then there exists C = C(Q,p)>0 such that, for all u\in L^p(\Omega)\otimes W^{1,p}_0(Q) , it holds

    \begin{equation*} \langle {{\int_Q|\nabla u|}} \rangle^p+\delta\langle {{\int_Q|\nabla u|^p}} \rangle\leq C\big(\mathcal E_{\varepsilon,\delta}(u)+1\big). \end{equation*}

    Proof. By Jensen's and Hölder's inequalities we have

    \begin{equation*} \langle {{\int_Q|\nabla u|dx}} \rangle^p \leq |Q|^{p-1}\int_Q\langle {{|\nabla u|}} \rangle^p\leq|Q|^{p-1}\langle {{{\lambda_{\varepsilon}^{-\frac{1}{p-1}}}}} \rangle^{p-1}\,\langle {{\int_Q\lambda_{\varepsilon}|\nabla u|^p}} \rangle, \end{equation*}

    where we use the notation \lambda_{\varepsilon}(x,\omega) = \lambda(\tau_{\frac{x}{\varepsilon}}\omega) . Furthermore, using (A3'), we conclude that

    \begin{equation*} \langle {{\int_Q|\nabla u|dx}} \rangle^p \leq C(Q,p)\left({{\mathcal{E}_{\varepsilon,\delta}(u) + 1}}\right). \end{equation*}

    In the end, using the variance-regularization we obtain

    \begin{eqnarray*} 2^{-p}\langle {{\int_Q|\nabla u|^p}} \rangle&\leq&\langle {{\int_Q|\nabla u-\langle {{\nabla u}} \rangle|^p}} \rangle+\int_Q\langle {{|\nabla u|}} \rangle^p\\ &\leq& \frac{C}{\delta}\left({{\mathcal E_{\varepsilon,\delta}(u)+1}}\right)+C\big(\mathcal E_{\varepsilon,\delta}(u)+1\big). \end{eqnarray*}

    This concludes the proof.

    The regularization on the \varepsilon -level is also consistent. In particular, we show that in the limit \delta \to 0 , we recover \mathcal{E}_{\varepsilon} . We discuss the mean functionals \mathcal{E}_{\varepsilon,\delta} and \mathcal{E}_{\varepsilon} , since the former does not admit a well-defined pointwise evaluation in \omega for the reason of the nonlocal variance term. Also, for the same reason the quenched version of stochastic two-scale convergence is not suitable for this setting and we apply the unfolding procedure. On the other hand, the homogenization of \mathcal{E}_{\varepsilon} can be conducted on the level of typical realizations, that was in fact studied in [29,20,21].

    Lemma 4.4. Let p\in (1,\infty) and Q\subset \mathbb{R}^d be open and bounded. Assume (A1), (A2) and (A3'). Then, \mathcal{E}_{\varepsilon, \delta} Mosco converges to \mathcal{E}_{\varepsilon} as \delta \to 0 i.e., the following statements hold:

    (i)If u_{\delta}\rightharpoonup u weakly in L^1(\Omega \times Q) , then

    \begin{equation*} \liminf\limits_{\delta \to 0}\mathcal{E}_{\varepsilon, \delta}(u_{\delta}) \geq \mathcal{E}_{\varepsilon}(u). \end{equation*}

    (ii)For any u \in L^1(\Omega \times Q) there exists a sequence u_{\delta}\in L^1(\Omega \times Q) such that

    \begin{equation*} u_{\delta} \to u \quad {{strongly\;in}}\;L^1(\Omega \times Q), \quad \mathcal{E}_{\varepsilon,\delta}(u_{\delta})\to \mathcal{E}_{\varepsilon}(u). \end{equation*}

    Proof. (i) Let u_{\delta} be a sequence such that u_{\delta}\rightharpoonup u weakly in L^1(\Omega \times Q) . Without loss of generality we assume that \limsup_{\delta \to 0}\mathcal{E}_{\varepsilon, \delta}(u_{\delta}) <\infty . This and the proof of Lemma 4.3 imply that the sequence \lambda_{\varepsilon}^{\frac{1}{p}} \nabla u_{\delta} is bounded in L^p(\Omega \times Q) with the notation \lambda_{\varepsilon}(x,\omega) = \lambda(\tau_{\frac{x}{\varepsilon}}\omega) . This means that, up to a subsequence, we have \lambda_{\varepsilon}^{\frac{1}{p}} \nabla u_{\delta} \rightharpoonup \psi weakly in L^p(\Omega \times Q) for some \psi \in L^p(\Omega \times Q) . Thus, for an arbitrary \eta \in L^{\infty}(\Omega\times Q) , we have

    \begin{equation*} \langle {{\int_{Q}\nabla u_{\delta} \eta dx}} \rangle = \langle {{\int_{Q}\lambda_{\varepsilon}^{\frac{1}{p}}\nabla u_{\delta} \lambda_{\varepsilon}^{-\frac{1}{p}}\eta dx}} \rangle\to \langle {{\int_{Q}\psi \lambda_{\varepsilon}^{-\frac{1}{p}}\eta dx}} \rangle \quad \text{as }\varepsilon \to 0. \end{equation*}

    This means that \nabla u_{\delta} converges weakly in L^1(\Omega \times Q) and since u_{\delta}\rightharpoonup u weakly in L^1(\Omega\times Q) we may conclude that \nabla u_{\delta}\rightharpoonup \nabla u weakly in L^1(\Omega\times Q) . This yields

    \begin{equation*} \liminf\limits_{\delta \to 0}\mathcal{E}_{\varepsilon, \delta}(u_{\delta}) \geq \liminf\limits_{\delta \to 0} \mathcal{E}_{\varepsilon}(u_{\delta}) \geq \mathcal{E}_{\varepsilon}(u). \end{equation*}

    (ii) For an arbitrary u\in \mathrm{dom}(\mathcal{E}_{\varepsilon})\subset X_{\varepsilon} , we find a sequence u_{\eta}\in L^p(\Omega)\otimes W^{1,p}_0(Q) such that, for \eta \to 0 ,

    \begin{equation*} u_{\eta} \to u \quad \text{strongly in }L^1(\Omega)\otimes W^{1,1}_0(Q), \quad \langle {{\int_{Q}\lambda_{\varepsilon} |\nabla u_{\eta}-\nabla u|^p dx}} \rangle \to 0. \end{equation*}

    Using this and the dominated convergence theorem, we conclude that

    \begin{equation*} \lim\limits_{\eta \to 0}\mathcal{E}_{\varepsilon}(u_{\eta}) = \mathcal{E}_{\varepsilon}(u). \end{equation*}

    This in turn yields

    \begin{equation*} \limsup\limits_{\eta\to 0} \limsup\limits_{\delta\to 0} |\mathcal{E}_{\varepsilon, \delta}(u_{\eta})- \mathcal{E}_{\varepsilon}(u)| = 0. \end{equation*}

    We extract a diagonal sequence \eta(\delta)\to 0 as \delta \to 0 such that u_{\delta}: = u_{\eta(\delta)} satisfies u_{\delta}\to u strongly in L^1(\Omega \times Q) and \mathcal{E}_{\varepsilon,\delta}(u_{\delta})\to \mathcal{E}_{\varepsilon}(u) . This concludes the proof.

    The homogenization of the regularized functional \mathcal{E}_{\varepsilon, \delta} boils down to a very similar simple argumentation as in Section 4.1.

    Theorem 4.5. Let p\in (1,\infty) and Q\subset \mathbb{R}^d be open and bounded. Assume (A1), (A2) and (A3'). For all \delta>0 , as \varepsilon \to 0 , \mathcal{E}_{\varepsilon, \delta} Mosco converges to \mathcal{E}_{\mathrm{hom}, \delta} in the following sense:

    (i)Let u_{ \varepsilon} \in L^p(\Omega)\otimes W^{1,p}_0(Q) be such that \limsup_{\varepsilon\rightarrow 0}\mathcal{E}_{\varepsilon, \delta}(u_{ \varepsilon})<\infty . Then there exist (u,\chi) \in W^{1,p}_0(Q) \times \left({{L^p_{ {\mathrm{pot}}}(\Omega)\otimes L^p(Q)}}\right) and a subsequence (not relabeled) such that

    \begin{equation*} u_{ \varepsilon} \overset{2}{\rightharpoonup} u \;\mathit{\text{in}}\; L^p(\Omega \times Q), \quad \nabla u_{ \varepsilon} \overset{2}{\rightharpoonup} \nabla u+\chi \;\mathit{\text{in}}\; L^p(\Omega \times Q). \end{equation*}

    (ii)If u_{ \varepsilon} \in L^1(\Omega \times Q) , u \in L^1(Q) and \mathcal{T}_{\varepsilon} u_{\varepsilon} \rightharpoonup u weakly in L^1(\Omega\times Q) , then

    \begin{equation*} \liminf\limits_{\varepsilon\to 0}\mathcal{E}_{\varepsilon, \delta}(u_{ \varepsilon}) \geq \mathcal{E}_{\hom,\delta}(u). \end{equation*}

    (iii)For any u \in L^1(Q) , there exists a sequence u_{ \varepsilon} \in L^1(\Omega \times Q) such that

    \begin{equation*} \mathcal{T}_{\varepsilon} u_{ \varepsilon} \to u \quad \;\mathit{\text{strongly in}}\;L^1(\Omega\times Q), \quad \mathcal{E}_{\varepsilon,\delta}(u_{ \varepsilon})\to \mathcal{E}_{\hom, \delta}(u). \end{equation*}

    Proof. (i) The statement follows analogously to the proof of Theorem 4.1 (i).

    (ii) Let \mathcal{T}_{\varepsilon} u_{\varepsilon}\to u weakly in L^1(\Omega\times Q) . We may assume without loss of generality that \limsup_{\varepsilon \to 0}\mathcal{E}_{\varepsilon,\delta}(u_{\varepsilon})<\infty . In this case, Lemma 4.3 implies that u_{ \varepsilon} is bounded in L^p(\Omega)\otimes W^{1,p}_0(Q) . We may proceed analogously to Theorem 4.1 and Remark 5 to obtain

    \begin{equation*} \liminf\limits_{\varepsilon\to 0}\mathcal{E}_{\varepsilon,\delta}(u_{\varepsilon}) \geq \mathcal{E}_{\hom, \delta}(u). \end{equation*}

    (ii) This part is analogous to Theorem 4.1 and Remark 5.

    The results of Lemmas (4.2) and (4.4), Theorem (4.5) and [29, 20,21] can be summarized in the following commutative diagram:

    \begin{equation*} \begin{array}{ccc} \qquad\quad \mathcal E_{\varepsilon,\delta} & \stackrel{(\delta \to 0)}{\to}& \mathcal E_{\varepsilon}\qquad \quad \\ {(\varepsilon \to 0)}{8pt} \downarrow & &\downarrow {(\varepsilon \to 0)}{8pt}\\ \qquad\quad \mathcal E_{\hom,\delta}&\stackrel{(\delta \to 0)}{\to}&\mathcal E_{\hom}\qquad \quad \end{array} \end{equation*}

    The arrows denote Mosco convergence in the corresponding convergence regimes.

    In this section we demonstrate how to lift homogenization results w.r.t. two-scale convergence in the mean to quenched statements at the example of Section 4.1. Throughout this section we assume that \langle {{\cdot}} \rangle is ergodic. For \omega\in\Omega we define {\mathcal{E}}^{\omega}_{ \varepsilon}: W^{1,p}_0(Q)\to \mathbb{R}^{{}} ,

    \begin{equation*} {\mathcal{E}}^{\omega}_{ \varepsilon}(u): = \int_{Q}V\left(\tau_{\frac{x}{ \varepsilon}}\omega, x,\nabla u(x)\right)\,dx, \end{equation*}

    with V satisfying (A1)-(A3). The goal of this section is to relate two-scale limits of "mean"-minimizers, i.e. functions u_{ \varepsilon}\in L^p(\Omega)\otimes W^{1,p}_0(Q) that minimize {\mathcal{E}}_{ \varepsilon} , with limits of "quenched"-minimizers, i.e. families \{u_{ \varepsilon}(\omega)\}_{\omega\in\Omega} of minimizers to {\mathcal{E}}^{\omega}_{ \varepsilon} in W^{1,p}_0(Q) . We also remark that if V(\omega,x,\cdot) is strictly convex u_{ \varepsilon} and \left\lbrace {{u_{ \varepsilon}(\omega)}} \right\rbrace_{\omega\in\Omega} may be identified since minimizers of both functionals \mathcal{E}_{ \varepsilon} and \mathcal{E}_{ \varepsilon}^{\omega} are unique.

    Before presenting the main result of this section, we remark that in the ergodic case, the limit functional (20) reduces to a single-scale energy

    \begin{equation*} \mathcal{E}_{\hom}: W^{1,p}_0(Q) \rightarrow \mathbb{R}^{{}}, \quad \mathcal{E}_{\hom}(u) = \int_Q V_{\hom}(x,\nabla u(x))dx, \end{equation*}

    where the homogenized integrand V_{\hom} is given for x\in \mathbb{R}^d and F\in \mathbb{R}^{d} by

    \begin{align} V_{\hom}(x,F) = \inf\limits_{\chi\in L^p_{ {\mathrm{pot}}}(\Omega)}\langle {{V(\omega, x,F+\chi(\omega))}} \rangle. \end{align} (31)

    In particular, we may obtain an analogous statement to Theorem 4.1 where we replace \mathcal{E}_0 with \mathcal{E}_{\mathrm{hom}} . The proof of this follows analogously with the only difference that in the construction of the recovery sequence we first need to find \chi such that \mathcal{E}_0(u,\chi) = \mathcal{E}_{\mathrm{hom}}(u) . This is done by a usual measurable selection argument, cf. [34,Theorem 7.6].

    Theorem 4.6. Let p\in (1,\infty) , Q\subset \mathbb{R}^d be open and bounded, and \langle {{\cdot}} \rangle be ergodic. Assume (A1)-(A3). Let u_{ \varepsilon}\in L^{p}(\Omega)\otimes W_{0}^{1,p}(Q) be a minimizer of {\mathcal{E}}_{ \varepsilon} . Then there exists a subsequence such that (u_{ \varepsilon},\nabla u_{ \varepsilon}) generates a Young measure \boldsymbol{\nu} in {\mathscr{B}}: = ( {\mathscr{B}}^p)^{1+d} in the sense of Theorem 3.12, and for P -a.a. \omega\in\Omega , \nu_{\omega} concentrates on the set \big\{\,(u,\nabla u+\chi)\,:\, {\mathcal{E}}_0(u,\chi) = \min {\mathcal{E}}_0\,\big\} of minimizers of the limit functional. Moreover, if V(\omega,x,\cdot) is strictly convex for all x\in Q and P -a.a. \omega\in\Omega , then the minimizer u_{ \varepsilon} of {\mathcal{E}}_{ \varepsilon} and the minimizer (u,\chi) of {\mathcal{E}}_0 are unique, and for P -a.a. \omega\in\Omega we have (for a not relabeled subsequence)

    \begin{gather*} u_{ \varepsilon}(\omega,\cdot) \rightharpoonup u\;{{weakly\;in}}\;W^{1,p}(Q),\qquad u_{ \varepsilon}(\omega,\cdot){\stackrel{2}{\rightharpoonup}_{{\omega}}}u,\qquad\nabla u_{ \varepsilon}(\omega,\cdot){\stackrel{2}{\rightharpoonup}_{{\omega}}}\nabla u+\chi,\\ \;{{and}}\;\min {\mathcal{E}}^\omega_ \varepsilon = {\mathcal{E}}^\omega_ \varepsilon(u_{ \varepsilon}(\omega,\cdot))\to {\mathcal{E}}_0(u,\chi) = \min {\mathcal{E}}_0. \end{gather*}

    Remark 6 (Identification of quenched two-scale cluster points). If we combine Theorem 4.6 with the identification of the support of the Young measure in Theorem 3.12 we conclude the following: There exists a subsequence such that (u_{ \varepsilon},\nabla u_{ \varepsilon}) two-scale converges in the mean to a limit of the form (u_0,\nabla u_0+\chi_0) with {\mathcal{E}}_0(u_0,\chi_0) = \min {\mathcal{E}}_0 , and for a.a. \omega\in\Omega the set of quenched \omega -two-scale cluster points {\mathscr{C\!P}}(\omega, (u_{ \varepsilon}(\omega,\cdot),\nabla u_{ \varepsilon}(\omega,\cdot))) is contained in \big\{\,(u,\nabla u+\chi)\,:\, {\mathcal{E}}_0(u,\chi) = \min {\mathcal{E}}_0\,\big\} . In the strictly convex case we further obtain that {\mathscr{C\!P}}(\omega, (u_{ \varepsilon}(\omega,\cdot),\nabla u_{ \varepsilon}(\omega,\cdot))) = \{(u,\nabla u+\chi)\} where (u,\chi) is the unique minimizer to {\mathcal{E}}_0 . Note, however, that our argument (that extracts quenched two-scale limits from the sequence of "mean" minimizers) involves an exceptional P -null-set that a priori depends on the selected subsequence. This is in contrast to the classical result in [11] which is based on a subadditive ergodic theorem and states that there exists a set of full measure \Omega' such that for all \omega\in\Omega' the minimizer u_{ \varepsilon}^\omega to {\mathcal{E}}^{\omega}_ \varepsilon weakly converges in W^{1,p}(Q) to the deterministic minimizer u of the reduced functional {\mathcal{E}}_{\hom} for any sequence \varepsilon\to 0 .

    In the proof of Theorem 4.6 we combine homogenization in the mean in form of Theorem 4.1, the connection to quenched two-scale limits via Young measures in form of Theorem 3.12, and a recent result described in Remark 3 by Nesenenko and the first author.

    Proof of Theorem 4.6. Step 1. (Identification of the support of \boldsymbol{\nu} ).

    Since u_{ \varepsilon} is a sequence of minimizers, by Corollary 4 there exists a subsequence (not relabeled) and minimizers (u,\chi)\in W^{1,p}_0(Q)\times (L^p_{ {\mathrm{pot}}}(\Omega)\otimes L^p(Q)) of {\mathcal{E}}_0 such that that u_{ \varepsilon} \overset{2}{\rightharpoonup} u \text{ in } L^p(\Omega \times Q) , \nabla u_{ \varepsilon} \overset{2}{\rightharpoonup} \nabla u+\chi \text{ in } L^p(\Omega \times Q)^{d} , and

    \begin{equation} \lim\limits_{ \varepsilon\to 0}\min {\mathcal{E}}_{ \varepsilon} = \lim\limits_{ \varepsilon\to 0} {\mathcal{E}}_{ \varepsilon}(u_{ \varepsilon}) = {\mathcal{E}}_0(u,\chi) = \min {\mathcal{E}}_0. \end{equation} (32)

    In particular, the sequence (u_{ \varepsilon},\nabla u_{ \varepsilon}) is bounded in {\mathscr{B}} . By Theorem 3.12 we may pass to a further subsequence (not relabeled) such that (u_{ \varepsilon},\nabla u_{ \varepsilon}) generates a Young measure \boldsymbol{\nu} on {\mathscr{B}} . Since \nu_\omega is supported on the set of quenched \omega -two-scale cluster points of (u_{ \varepsilon}(\omega,\cdot),\nabla u_{ \varepsilon}(\omega,\cdot)) , we deduce from Lemma 3.10 that the support of \nu_\omega is contained in {\mathscr{B}}_0: = \{\xi = (\xi_1,\xi_2) = (u',\nabla u'+\chi')\,:\,u'\in W^{1,p}_0(Q),\,\chi\in L^p_{ {\mathrm{pot}}}(\Omega)\otimes L^p(Q)\} which is a closed subspace of {\mathscr{B}} . Moreover, thanks to the relation of the generated Young measure and stochastic two-scale convergence in the mean, we have (u,\chi) = \int_\Omega \int_{ {\mathscr{B}}_0}(\xi_1,\xi_2-\nabla\xi_1)\,\nu_\omega(d\xi)\,dP(\omega) . Furthermore, Lemma 3.14 implies that

    \begin{equation*} \lim\limits_{ \varepsilon\to 0} {\mathcal{E}}_ \varepsilon(u_{ \varepsilon})\geq \int_\Omega\int_{ {\mathscr{B}}}\Big(\int_\Omega\int_Q V( \tilde{\omega},x,\xi_2)\,dx\,dP( \tilde{\omega})\Big)\,\nu_\omega(d\xi)\,dP(\omega). \end{equation*}

    In view of (32) and the fact that \nu_\omega is supported in {\mathscr{B}}_0 , we conclude that

    \begin{equation*} \min {\mathcal{E}}_0\geq \int_\Omega\int_{ {\mathscr{B}}_0} {\mathcal{E}}_0(\xi_1,\xi_2-\nabla\xi_1)\,\nu_\omega(d\xi)\,dP(\omega)\geq \min {\mathcal{E}}_0\int_\Omega\int_{ {\mathscr{B}}_0}\nu_\omega(d\xi)dP(\omega). \end{equation*}

    Since \int_\Omega\int_{ {\mathscr{B}}_0}\nu_\omega(d\xi)dP(\omega) = 1 , we have \int_\Omega\int_{ {\mathscr{B}}_0}| {\mathcal{E}}_0(\xi_1,\xi_2-\nabla\xi_1)-\min {\mathcal{E}}_0|\,\nu_\omega(d\xi)\,dP(\omega) = 0 , and thus we conclude that for P -a.a. \omega\in\Omega_0 , \nu_\omega concentrates on \{(u,\nabla u+\chi)\,:\, {\mathcal{E}}_0(u,\chi) = \min {\mathcal{E}}_0\} .

    Step 2. (The strictly convex case).

    The uniqueness of u_{ \varepsilon} and (u,\chi) is clear. From Step 1 we thus conclude that \nu_\omega = \delta_{\xi} where \xi = (u,\nabla u+\chi) . Theorem 3.12 implies that (u_{ \varepsilon}(\omega,\cdot),\nabla u_{ \varepsilon}(\omega,\cdot)){\stackrel{2}{\rightharpoonup}_{{\omega}}}(u,\nabla u+\chi) (for P -a.a. \omega\in\Omega ). By Lemma 3.14 we have for P -a.a. \omega\in\Omega ,

    \begin{equation*} \liminf\limits_{ \varepsilon\to 0} {\mathcal{E}}^\omega_ \varepsilon(u_{ \varepsilon}(\omega,\cdot))\geq {\mathcal{E}}_0(u,\chi) = \min {\mathcal{E}}_0. \end{equation*}

    On the other hand, since u_{ \varepsilon}(\omega,\cdot) minimizes {\mathcal{E}}^\omega_ \varepsilon , we deduce by a standard argument that for P -a.a. \omega\in\Omega ,

    \begin{equation*} \lim\limits_{ \varepsilon\to 0}\min {\mathcal{E}}^\omega_ \varepsilon = \lim\limits_{ \varepsilon\to 0} {\mathcal{E}}^\omega_ \varepsilon(u_{ \varepsilon}(\omega,\cdot)) = {\mathcal{E}}_0(u,\chi) = \min {\mathcal{E}}_0. \end{equation*}

    The authors thank Alexander Mielke for fruitful discussions and valuable comments. MH has been funded by Deutsche Forschungsgemeinschaft (DFG) through grant CRC 1114 "Scaling Cascades in Complex Systems", Project C05 "Effective models for materials and interfaces with multiple scales". SN and MV acknowledge funding by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – project number 405009441.



    [1] I. Gutman, The energy of a graph, Ber. Math. Statist. Sekt. Forschungsz. Graz, 103 (1978), 1-22. https://doi.org/10.1002/9783527627981
    [2] D. Cvetkovic, M. Doob, H. Sachs, Spectra of Graphs Theory and Applications, Academic Press, 1982.
    [3] D. Cvetkovic, P. Rowlinson, The largest eigenvalue of a graph, Linear Multilinear Algebra, 28 (1990), 3-33. https://doi.org/10.1080/03081089008818026 doi: 10.1080/03081089008818026
    [4] Y. Hong, Upper bounds of the spectral radius of graphs in terms of genus, J. Combin. Theory Ser., 74 (1998), 153-159. https://doi.org/10.1006/jctb.1998.1837 doi: 10.1006/jctb.1998.1837
    [5] R. A. Brualdi, E. S. Solheid, On the spectral radius of complementary acyclic matrices of zeros and ones, SIAM J. Algebraic Discrete Methods, 7 (1986), 265-272. https://doi.org/10.1137/0607030 doi: 10.1137/0607030
    [6] M. Fiedler, V. Nikiforov, Spectral radius and Hamiltonicity of graphs, Linear Algebra Appl., 432 (2010), 2170-2173. https://doi.org/10.1016/j.laa.2009.01.005 doi: 10.1016/j.laa.2009.01.005
    [7] M. Lu, H. Liu, F. Tian, Spectral radius and Hamiltonian graphs, Linear Algebra Appl., 437 (2012), 1670-1674. https://doi.org/10.1016/j.laa.2012.05.021 doi: 10.1016/j.laa.2012.05.021
    [8] V. Nikiforov, More spectral bounds on the clique and independence numbers, J. Combin. Theory Ser., 99 (2009), 819-826. https://doi.org/10.1016/j.jctb.2009.01.003 doi: 10.1016/j.jctb.2009.01.003
    [9] V. Nikiforov, The spectral radius of graphs without paths and cycles of specified length, Linear Algebra Appl., 432 (2010), 2243-2256. https://doi.org/10.1016/j.laa.2009.05.023 doi: 10.1016/j.laa.2009.05.023
    [10] W. Yuan, B. Wang, M. Zhai, On the spectral radii of graphs without given cycles, Electron. J. Linear Algebra, 23 (2012), 599-606. https://doi.org/10.13001/1081-3810.1544 doi: 10.13001/1081-3810.1544
    [11] M. Zhai, B. Wang, Proof of a conjecture on the spectral radius of C4-free graphs, Linear Algebra Appl., 437 (2012), 1641-1647. https://doi.org/10.1016/j.laa.2012.05.006 doi: 10.1016/j.laa.2012.05.006
    [12] S. Butler, F. Chung, Small spectral gap in the combinatorial Laplacian implies Hamiltonian, Ann. Comb., 13 (2010), 403-412. https://doi.org/10.1007/s00026-009-0039-4 doi: 10.1007/s00026-009-0039-4
    [13] J. V. Heuvel, Hamilton cycles and eigenvalues of graphs, Linear Algebra Appl., 228 (1995), 723-730. https://doi.org/10.1016/0024-3795(95)00254-O doi: 10.1016/0024-3795(95)00254-O
    [14] M. Krivelevich, B. Sudakov, Sparse pseudo-random graphs are Hamiltonian, J. Graph Theory, 42 (2003), 17-33. https://doi.org/10.1002/jgt.10065 doi: 10.1002/jgt.10065
    [15] B. Mohar, A domain monotonicity theorem for graphs and hamiltonicity, Discrete Appl. Math., 36 (1992), 169-177. https://doi.org/10.1016/0166-218X(92)90230-8 doi: 10.1016/0166-218X(92)90230-8
    [16] B. Zhou, Signless Laplacian spectral radius and Hamiltonicity, Linear Algebra Appl., 432 (2010), 566-570. https://doi.org/10.1016/j.laa.2009.09.004 doi: 10.1016/j.laa.2009.09.004
    [17] R. A. Horn, C. R. Johnson, Topics in Matrix Analysis, Cambridge Univ. Press, 1991. https://doi.org/10.1017/CBO9780511840371
    [18] F. R. Gantmacher, The Theory of Matrices, Chelsea, 1959.
    [19] R. Balakrishnan, The energy of a graph, Lin. Algebra Appl., 387 (2004), 287-295. https://doi.org/10.1016/j.laa.2004.02.038
    [20] R. B. Bapat, S. Pati, Energy of a graph is never an odd integer, Bull. Kerala Math. Assoc., 1 (2004), 129-132. https://doi.org/10.2298/AADM0801118P doi: 10.2298/AADM0801118P
    [21] S. Pirzada, S. Gutman, Energy of a graph is never the square root of an odd integer, Appl. Anal. Discr. Math., 2 (2008), 118-121. https://doi.org/10.2298/AADM0801118P doi: 10.2298/AADM0801118P
    [22] O. Jones, Spectra of simple graphs, Whitman College, 13 (2013), 1-20.
    [23] S. Meenaksh, S. Lavanya, A Survey on Energy of Graphs, Ann. Pure Appl. Math., 8 (2014), 183-191.
    [24] V. Nikiforov, The energy of graphs and matrices, J. Math. Anal. Appl., 326 (2007), 1472-1475. https://doi.org/10.1016/j.jmaa.2006.03.072 doi: 10.1016/j.jmaa.2006.03.072
    [25] K. Samir, Vaidya, M. P. Kalpesh, Some new results on energy of graphs, MATCH Common. Math Comput. Chem, 77 (2017), 589-594.
    [26] K. Samir, Vaidya, M. P. Kalpesh, On energy m-splitting and m-shadow graphs, Far East J. Of Math. Sci., 102 (2017), 1571-1578. https://doi.org/10.17654/MS102081571 doi: 10.17654/MS102081571
    [27] J. B. Liu, M. Munir, A. Yousaf, A. Naseem, K. Ayub, Distance and adjacency energies of multi-level wheel networks, Mathematics, 7 (2019), 43. https://doi.org/10.3390/math7010043 doi: 10.3390/math7010043
    [28] Z. Q. Chu, M. Munir, A. Yousaf, M. I. Qureshi, J, B. Liu, Laplacian and signless laplacian spectra and energies of multi-step wheels, Math. Biosci. Eng., 17 (2020), 3649-3659. https://doi.org/10.3934/mbe.2020206 doi: 10.3934/mbe.2020206
    [29] J. B. Liu, X. F. Pan, F. T. Hu, F. HuF, On Asymptotic Laplacian-energy-like invariant of lattices, Appl. Math. Comput., 253 (2015), 205-214. https://doi.org/10.1016/j.amc.2014.12.035 doi: 10.1016/j.amc.2014.12.035
    [30] S. M. Hosamani, H. S. Ramane, On degree sum energy of a graph, Eur. J. Pure Appl. Math., 9 (2016), 340-345.
    [31] B. Basavanagoud, E. Chitra, On degree square sum energy of graphs, Int. J. Math. Its Appl., 6 (2018), 193-205.
    [32] N. J. Rad, A. Jahanbani, I. Gutman, Zagreb energy and Zagreb Estrada index of graphs, MATCH Commun. Math. Comput. Chem., 79 (2018), 371-386.
    [33] I. Gutman, B. Furtula, On graph energies and their application, Bulletin (Acadmie Serbe Des Sciences Et Des Arts. Classe Des Sciences Mathmatiques Et Naturelles. SciencesMath matiques), 44 (2019), 29-45.
    [34] X. Li, Y. Shi, I. Gutman, Graph Energy, Springer, New York, 2012. https://doi.org/10.1007/978-1-4614-4220-2
    [35] R. Balakrishnan, K. Ranganathan, A Textbook of Graph Theory, Springer, 2000. https://doi.org/10.1007/978-1-4419-8505-7
    [36] R. J. Wilson, Introduction to Graph Theory, Oliver and Boyd, 1972.
    [37] J. J. Sylvester, Chemistry and algebra, Nature, 17 (1877), 284. https://doi.org/10.1038/017284a0
    [38] A. T. Balaban, Applications of graph theory in chemistry, J. Chem. Inf. Comput. Sci., 25 (1985), 334-343. https://doi.org/10.1021/ci00047a033 doi: 10.1021/ci00047a033
    [39] M. RandiC, On characterization of molecular branching, J. Amer. Chem. Soc., 97 (1975), 6609-6615. https://doi.org/10.1021/ja00856a001 doi: 10.1021/ja00856a001
    [40] K. Yuge, Extended configurational polyhedra based on graph representation for crystalline solids, Trans. Material Res. Soc., 43 (2018), 233-236. https://doi.org/10.14723/tmrsj.43.233 doi: 10.14723/tmrsj.43.233
    [41] A. Dhanalakshmi, K. S. Rao, K. Sivakumar, Characterization of a-cyclodextrin using adjacency and distance matrix, Indian J. Sci., 12 (2015), 78-83.
    [42] J. Praznikar, M. Tomic, D. Turk, Validation and quality assessment of macromolecular structures using complex network analysis, Sci. Rep., 9 (2019), 1678. https://doi.org/10.1038/s41598-019-38658-9 doi: 10.1038/s41598-019-38658-9
    [43] H. Wu, Y. Zhang, W. Chen, Z. Mu, Comparative analysis of protein primary sequences with graph energy, Phys. A, 437 (2015), 249-262. https://doi.org/10.1016/j.physa.2015.04.017 doi: 10.1016/j.physa.2015.04.017
    [44] L. D. Paola, G. Mei, A. D. Venere, A. Giuliani, Exploring the stability of dimers through protein structure topology, Curr. Protein Peptide Sci., 17 (2016), 30-36. https://doi.org/10.2174/1389203716666150923104054 doi: 10.2174/1389203716666150923104054
    [45] D. Sun, C. Xu, Y. Zhang, A novel method of 2D graphical representation for proteins and its application, MATCH Commun. Math. Comput. Chem., 75 (2016), 431-446.
    [46] L. Yu, Y. Z. Hang, I. Gutman, Y. Shi, M. Dehmer, Protein sequence comparison based on physicochemical properties and positionfeature energy matrix, Sci. Rep., 7 (2017), 46237. https://doi.org/10.1038/srep46237 doi: 10.1038/srep46237
    [47] A. Giuliani, S. Filippi, M. Bertolaso, Why network approach can promote a new way of thinking in biology, Front. Genet., 5 (2014), 83. https://doi.org/10.3389/fgene.2014.00083 doi: 10.3389/fgene.2014.00083
    [48] J. Jiang, R. Zhang, L. Guo, W. Li, X. Cai, Network aggregation process in multilayer air transportation networks, Chin. Phys. Lett., 33 (2016), 108901. https://doi.org/10.1088/0256-307X/33/10/108901 doi: 10.1088/0256-307X/33/10/108901
    [49] T. A. Shatto, E. K. C. Etinkaya, Variations in graph energy: A measure for network resilience, in 2017 9th International Workshop on Resilient Networks Design and Modeling (RNDM), (2017), 1-7. https://doi.org/10.1109/RNDM.2017.8093019
    [50] M. Akram, S. Naz, Energy of Pythagorean fuzzy graphs with applications, Mathematics, 6 (2018), 136. https://doi.org/10.3390/math6080136 doi: 10.3390/math6080136
    [51] M. Gao, E. P. Lim, D. Lo, Network Data Mining and Analysis, World Scientific, 4 (2018), 1-4. https://doi.org/10.1142/9789813274969-0001
    [52] A. Pugliese, R. Nilchiani, Complexity analysis of fractionated spacecraft architectures, in AIAA SPACE and Astronautics Forum and Exposition, (2017), 5118-5126. https://doi.org/10.2514/6.2017-5118
    [53] C. Adiga, M. Smitha, On maximum degree energy of a graph, Int. J. Contemp. Math. Sci., 4 (2009), 385-396.
    [54] B. Basavanagoud, P. Jakkannavar, Minimum degree energy pf graphs, Electron. J. Math. Anal. Appl., 7 (2019), 230-243.
    [55] S. Alikhani, N. Ghanbari, Randic energy of specific graphs, Appl. Math and Comp, 269 (2015), 722-730. https://doi.org/10.1016/j.amc.2015.07.112 doi: 10.1016/j.amc.2015.07.112
    [56] M. R. Oboudi, Energy and seidel energy of graphs, MATCH Commun. Math. Comput. Chem., 75 (2016), 291-303.
  • This article has been cited by:

    1. Thuyen Dang, Yuliya Gorb, Silvia Jiménez Bolaños, Homogenization of a NonLinear Strongly Coupled Model of Magnetorheological Fluids, 2023, 55, 0036-1410, 102, 10.1137/22M1476939
    2. S. Yu. Dobrokhotov, V. E. Nazaikinskii, Averaging method for problems on quasiclassical asymptotics, 2024, 70, 2949-0618, 53, 10.22363/2413-3639-2024-70-1-53-76
    3. Patrick Dondl, Yongming Luo, Stefan Neukamm, Steve Wolff-Vorbeck, Efficient Uncertainty Quantification for Mechanical Properties of Randomly Perturbed Elastic Rods, 2024, 22, 1540-3459, 1267, 10.1137/23M1574233
    4. Sabine Haberland, Patrick Jaap, Stefan Neukamm, Oliver Sander, Mario Varga, Representative Volume Element Approximations in Elastoplastic Spring Networks, 2024, 22, 1540-3459, 588, 10.1137/23M156656X
    5. S. Yu. Dobrokhotov, V. E. Nazaikinskii, Homogenization Method for Problems on Quasiclassical Asymptotics, 2024, 1072-3374, 10.1007/s10958-024-07490-6
  • Reader Comments
  • © 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2960) PDF downloads(138) Cited by(9)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog