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

Pricing stock loans under the Lˊevy-α-stable process with jumps

  • Received: 16 September 2022 Revised: 20 November 2022 Accepted: 21 November 2022 Published: 30 November 2022
  • In this paper, the pricing of stock loans when the underlying follows a Lˊevy-α-stable process with jumps is considered. Under this complicated model, the stock loan value satisfies a fractional-partial-integro-differential equation (FPIDE) with a free boundary. The difficulties in solving the resulting FPIDE system are caused by the non-localness of the fractional-integro differential operator, together with the nonlinearity resulting from the early exercise opportunity of stock loans. Despite so many difficulties, we have managed to propose a preconditioned conjugate gradient normal residual (PCGNR) method to price efficiently the stock loan under such a complicated model. In the proposed approach, the moving pricing domain is successfully dealt with by introducing a penalty term, however, the semi-globalness of the fractional-integro operator is elegantly handled by the PCGNR method together with the fast Fourier transform (FFT) technique. Remarkably, we show both theoretically and numerically that the solution determined from the fixed domain problem by the current method is always above the intrinsic value of the corresponding option. Numerical experiments suggest the accuracy and advantage of the current approach over other methods that can be compared. Based on the numerical results, a quantitative discussion on the impacts of key parameters is also provided.

    Citation: Congyin Fan, Wenting Chen, Bing Feng. Pricing stock loans under the Lˊevy-α-stable process with jumps[J]. Networks and Heterogeneous Media, 2023, 18(1): 191-211. doi: 10.3934/nhm.2023007

    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
  • In this paper, the pricing of stock loans when the underlying follows a Lˊevy-α-stable process with jumps is considered. Under this complicated model, the stock loan value satisfies a fractional-partial-integro-differential equation (FPIDE) with a free boundary. The difficulties in solving the resulting FPIDE system are caused by the non-localness of the fractional-integro differential operator, together with the nonlinearity resulting from the early exercise opportunity of stock loans. Despite so many difficulties, we have managed to propose a preconditioned conjugate gradient normal residual (PCGNR) method to price efficiently the stock loan under such a complicated model. In the proposed approach, the moving pricing domain is successfully dealt with by introducing a penalty term, however, the semi-globalness of the fractional-integro operator is elegantly handled by the PCGNR method together with the fast Fourier transform (FFT) technique. Remarkably, we show both theoretically and numerically that the solution determined from the fixed domain problem by the current method is always above the intrinsic value of the corresponding option. Numerical experiments suggest the accuracy and advantage of the current approach over other methods that can be compared. Based on the numerical results, a quantitative discussion on the impacts of key parameters is also provided.



    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 is a Young measure on , if for all , is a Borel probability measure on (equipped with the weak topology) and

    where denotes the Borel--algebra on (equipped with the weak topology).

    Theorem 3.12. Let denote a bounded sequence in .Then there exists a subsequence (still denoted by ) and a Young measure on such that for all ,

    and

    Moreover, we have

    Finally, if there exists measurable and for -a.a. , then up to extraction of a further subsequence (still denoted by ) we have

    (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 of sequences in and suppose that:

    (i)There exists s.t. for -a.a. we have .

    (ii)There exists a sequence in s.t. for a.a. .

    (iii)There exists a bounded sequence in such that for a.a. .

    Then 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 be such that for all , is -measurable and for a.a. , is convex. Let denote a bounded sequence in that generates a Young measure on in the sense of Theorem 3.12.Suppose that , is uniformly integrable. Then

    (11)

    (For the proof see Section 3.3.1).

    Remark 3. In [18,Lemma 5.1] it is shown that satisfying the assumptions of Lemma 3.14 satisfies the following property: For -a.a. we have: For any sequence in it holds

    (12)

    We first recall some notions and results of Balder's theory for Young measures [4]. Throughout this section is assumed to be a separable, complete metric space with metric .

    Definition 3.15. ● We say a function is measurable, if it is -measurable where denotes the Borel--algebra in .

    ● A function is called a normal integrand, if is -measurable, and for all the function is lower semicontinuous.

    ● A sequence of measurable functions is called tight, if there exists a normal integrand such that for every the function has compact sublevels in and .

    ● A Young measure in is a family of Borel probability measures on such that for all the map is -measurable.

    Theorem 3.16.([4,Theorem I]). Let denote a tight sequence of measurable functions. Then there exists a subsequence, still indexed by , and a Young measure such that for every normal integrand we have

    (13)

    provided that the negative part is uniformly integrable.Moreover, for -a.a. the measure is supported in the set of all cluster points of , i.e. in (where the closure is taken in ).

    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 the set of all triples with , , . endowed with the metric

    is a complete, separable metric space.

    (ii)For we denote by the set of all triples such that

    (14)

    Then is a closed subspace of .

    (iii)Let , and . Then the function in the representation (14) of is unique, and

    (15)

    (iv)For the function ,

    is lower semicontinuous on .

    (v)For all with and we have and . Likewise, for all with and we have and .

    (vi)For all the set is compact in .

    (vii)Let and let denote a bounded sequence in . Then the triple defines a sequence in . Moreover, we have in as if and only if for some , , and .

    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 denote a sequence in that converges in to some . We need to show that . By passing to a subsequence, it suffices to study the following three cases: for all , for all , and while for all .

    Case 1: for all .

    W.l.o.g. we may assume that . Hence, there exist with . Since , we conclude that is bounded in . We thus may pass to a subsequence (not relabeled) such that weakly in , and

    (16)

    Moreover, in the metric space implies pointwise convergence on , and thus for all we have . We thus conclude that . Since dense, we deduce that weakly in for the entire sequence. On the other hand the properties of the shift imply that for any we have in . Hence, for any and we have

    and thus (by linearity) .

    Case 2: for all .

    In this case there exist a bounded sequence in with for . By passing to a subsequence we may assume that weakly in for some with

    (17)

    This implies that in . Hence, and we conclude that .

    Case 3: for all and .

    There exists a bounded sequence in . Thanks to Lemma 3.5, by passing to a subsequence we may assume that for some with

    (18)

    Furthermore, Lemma 3.6 implies that in , and thus . We conclude that .

    (iii)We first argue that the representation (14) of by is unique. In the case suppose that satisfy . Then for all we have , and since dense, we conclude that . In the case the statement follows by a similar argument from the fact that is dense .

    To see (15) let and be related by (14). Since (resp. ) is dense in (resp. ), we have

    (iv)Let denote a sequence in that converges in to a limit . By (ii) we have . For let in or denote the representation of 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 weakly converges to (in or ), or . In any of these cases the claimed lower semicontinuity of follows from , , 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 denote a sequence in . Let in or denote the (unique) representative of in the sense of (14). Suppose that . Then and are bounded sequences in , and (where stands short for either or ). Thus we may pass to a subsequence such that , , and one of the following three cases applies:

    ● Case 1: . In that case we conclude (after passing to a further subsequence) that weakly in , and thus in .

    ● Case 2: for all . In that case we conclude (after passing to a further subsequence) that weakly in , and thus in .

    ● Case 3: for all and . In that case we conclude (after passing to a further subsequence) that , and thus in .

    In all of these cases we deduce that , and in .

    (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 , , etc. be defined as in Lemma 3.17.

    Step 1. (Identification of with a tight -valued sequence). Since , by Fubini's theorem, we have for -a.a. . By modifying on a null-set in (which does not alter two-scale limits in the mean), we may assume w.l.o.g. that for all . Consider the measurable function defined as

    We claim that is tight. To that end consider the integrand defined by

    From Lemma 3.17 (iv) and (vi) we deduce that is a normal integrand and has compact sublevels for all . Moreover, for all we have and thus . Hence,

    We conclude that is tight.

    Step 2. (Compactness and definition of ). By appealing to Theorem 3.16 there exists a subsequence (still denoted by ) and a Young measure that is generated by . Let denote the first component of , i.e. the Young measure on characterized for by

    for all continuous and bounded, where denotes the projection to the first component. By Balder's theorem, is concentrated on the limit points of . By Lemma 3.17 we deduce that for all any limit point of has the form where and is a -two-scale limit of a subsequence of . Thus, is supported on which in particular is a subset of . Since is an isometric isomorphism (by the Riesz-Frechét theorem), we conclude that , (for all Borel sets where is equipped with the weak topology) defines a Young measure on , and for all , is supported on .

    Step 3. (Lower semicontinuity estimate). Note that ,

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

    In view of Lemma 3.17 we have for , and thus the left-hand side turns into . Thanks to the definition of the right-hand side turns into .

    Step 4. (Identification of the two-scale limit in the mean). Let . Then ,

    defines a normal integrand. Since for -a.a. , we deduce that is uniformly integrable. Thus, (13) applied to and the definition of imply that

    Set . Then Fubini's theorem yields

    Since dense, we conclude that .

    Step 5. Recovery of quenched two-scale convergence. Suppose that is a delta distribution on , say for some measurable . Note that ,

    is a normal integrand and is uniformly integrable. Thus, (13) yields

    Thus, there exists a subsequence (not relabeled) such that for a.a. . In view of Lemma 3.6 this implies that for a.a. .

    Proof of Lemma 3.14. Step 1. Representation of the functional by a lower semicontinuous integrand on .

    For all and we write for the unique representation in (resp. ) of in the sense of (14). We thus may define for and the integrand

    We extend to by , and define for . We claim that is lower semicontinuous for all . It suffices to consider and a convergent sequence in . For brevity we only consider the (interesting) case when . Set . By construction we have

    and

    Since and , Lemma 3.17 (vi) implies that , and since satisfies 12 from Remark 3, we conclude that , and thus is a normal integrand.

    Step 2. Conclusion. As in Step 1 of the proof of Theorem 3.12 we may associate with the sequence a sequence of measurable functions that (after passing to a subsequence that we do not relabel) generates a Young measure on . Since by assumption generates the Young measure on , we deduce that the first component satisfies for any Borel set . Applying (13) to the integrand of Step 1, yields

    Proof of Lemma 3.13. By (ii) and (iii) the sequence is bounded in and thus we can pass to a subsequence such that generates a Young measure . Set and note that Theorem 3.12 implies that weakly two-scale in the mean. On the other hand the theorem implies that concentrates on the quenched two-scale cluster points of (for a.a. ). Hence, in view of (i) we conclude that for a.a. the measure is a Dirac measure concentrated on , and thus a.e. in .

    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 . A realization of a stationary and ergodic random checkerboard is given by

    where is the integer part of and is the center of the checkerboard chosen uniformly from . For , we may consider the map given by for and is -periodically extended. The push forward of the map 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 and be open and bounded. We consider and assume:

    (A1) is -measurable for all .

    (A2) is convex for a.a. .

    (A3)There exists such that

    for a.a. and all .

    We consider the functional

    (19)

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

    (20)

    Theorem 4.1 (Two-scale homogenization). Let and be open and bounded. Assume (A1)-(A3).

    (i)(Compactness and liminf inequality.) Let be such that.There exist and a subsequence (not relabeled) such that

    (21)
    (22)

    (ii)(Limsup inequality.) Let . There exists a sequence such that

    (23)
    (24)

    Proof of Theorem 4.1. (i) The Poincaré inequality and (A3) imply that is bounded in . By Proposition 1 (ii) there exist and such that (21) holds. Also, note that weakly in and , which implies that also has boundary values, i.e., . Finally, we note that, see [19,Proposition 3.5 (i)],

    (25)

    and thus using the convexity of we conclude

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

    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 be a minimizer of the functional

    where and with and . By a standard argument from the theory of -convergence Theorem 4.1 (cf. [34,Corollary 7.2]) implies that there exist a subsequence (not relabeled), , and such that , , and

    where is given by . This, in particular, rigorously justifies the formal two-scale expansion .

    Remark 5 (Uniqueness). If 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 that satisfies (A1), (A2) and the following assumption (as a replacement of (A3)):

    (A3')There exists and a random variable such that

    (26)

    and

    for a.a. and all .

    Moreover, we assume that is ergodic. For we consider the following functional

    for and otherwise. Here denotes the closure of w.r.t. the weighted norm

    Recently, in [29,20,21] it shown that Mosco-converges to the functional

    for and otherwise, where is given by the homogenization formula,

    (27)

    for and . Moreover, it is shown that is a normal convex integrand that satisfies a standard -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 . It is thus easy to see that the homogenized integrand satisfies -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 in (27) is non-degenerate, in particular, in [29,Theorem 3.1] it is shown that there exists a constant such that for all and we have

    (28)

    One of the difficulties in the proof of the homogenization result for is due to the fact that the domain of the functionals are -dependent. Moreover, assumption (A3') only yields equicoercivity in , while the limit is properly defined on . Therefore, in practice it is convenient to regularize the problem: For we consider the regularized homogenization formula

    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

    for and otherwise. Furthermore, thanks to (A3'), it is relatively easy to see that this regularization is consistent:

    Lemma 4.2. Let and be open and bounded. Assume (A1), (A2) and (A3'). Then, for all and , we have

    (29)

    Moreover, Mosco converges to as , i.e., the following statements hold:

    (i)If weakly in , then

    (ii)For any there exists a sequence such that

    Proof. Let and . Since , we have . On the other hand, we consider a minimizing sequence in (27), e.g.,

    We have

    Letting first and then , we conclude (29).

    We further consider a sequence such that weakly in as . We assume without loss of generality that. This, in particular, with the help of (28) and the Poincaré inequality implies that . Thus, up to a subsequence, we have weakly in . Using this, we obtain

    The first inequality follows by (29) and the second is a consequence of the fact that is convex and of Fatou's Lemma. We conclude that (i) holds.

    If , we simply choose . On the other hand, for , (29) and the dominated convergence theorem yield

    This means that (ii) holds.

    In the following we introduce a variance regularization of the original functional that removes the degeneracy of the problem and thus can be analyzed by the standard strategy of Section 4.1. For , we consider ,

    (30)

    for and 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 , the functional is equicoercive on :

    Lemma 4.3. Let and be open and bounded. Assume (A1) and (A3'). Then there exists such that, for all , it holds

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

    where we use the notation . Furthermore, using (A3'), we conclude that

    In the end, using the variance-regularization we obtain

    This concludes the proof.

    The regularization on the -level is also consistent. In particular, we show that in the limit , we recover . We discuss the mean functionals and , since the former does not admit a well-defined pointwise evaluation in 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 can be conducted on the level of typical realizations, that was in fact studied in [29,20,21].

    Lemma 4.4. Let and be open and bounded. Assume (A1), (A2) and (A3'). Then, Mosco converges to as i.e., the following statements hold:

    (i)If weakly in , then

    (ii)For any there exists a sequence such that

    Proof. (i) Let be a sequence such that weakly in . Without loss of generality we assume that . This and the proof of Lemma 4.3 imply that the sequence is bounded in with the notation . This means that, up to a subsequence, we have weakly in for some . Thus, for an arbitrary , we have

    This means that converges weakly in and since weakly in we may conclude that weakly in . This yields

    (ii) For an arbitrary , we find a sequence such that, for ,

    Using this and the dominated convergence theorem, we conclude that

    This in turn yields

    We extract a diagonal sequence as such that satisfies strongly in and . This concludes the proof.

    The homogenization of the regularized functional boils down to a very similar simple argumentation as in Section 4.1.

    Theorem 4.5. Let and be open and bounded. Assume (A1), (A2) and (A3'). For all , as , Mosco converges to in the following sense:

    (i)Let be such that. Then there exist and a subsequence (not relabeled) such that

    (ii)If , and weakly in , then

    (iii)For any , there exists a sequence such that

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

    (ii) Let weakly in . We may assume without loss of generality that . In this case, Lemma 4.3 implies that is bounded in . We may proceed analogously to Theorem 4.1 and Remark 5 to obtain

    (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:

    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 is ergodic. For we define ,

    with satisfying (A1)-(A3). The goal of this section is to relate two-scale limits of "mean"-minimizers, i.e. functions that minimize , with limits of "quenched"-minimizers, i.e. families of minimizers to in . We also remark that if is strictly convex and may be identified since minimizers of both functionals and 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

    where the homogenized integrand is given for and by

    (31)

    In particular, we may obtain an analogous statement to Theorem 4.1 where we replace with . The proof of this follows analogously with the only difference that in the construction of the recovery sequence we first need to find such that . This is done by a usual measurable selection argument, cf. [34,Theorem 7.6].

    Theorem 4.6. Let , be open and bounded, and be ergodic. Assume (A1)-(A3). Let be a minimizer of . Then there exists a subsequence such that generates a Young measure in in the sense of Theorem 3.12, and for -a.a. , concentrates on the set of minimizers of the limit functional. Moreover, if is strictly convex for all and -a.a. , then the minimizer of and the minimizer of are unique, and for -a.a. we have (for a not relabeled subsequence)

    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 two-scale converges in the mean to a limit of the form with , and for a.a. the set of quenched -two-scale cluster points is contained in . In the strictly convex case we further obtain that where is the unique minimizer to . Note, however, that our argument (that extracts quenched two-scale limits from the sequence of "mean" minimizers) involves an exceptional -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 such that for all the minimizer to weakly converges in to the deterministic minimizer of the reduced functional for any sequence .

    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 ).

    Since is a sequence of minimizers, by Corollary 4 there exists a subsequence (not relabeled) and minimizers of such that that , , and

    (32)

    In particular, the sequence is bounded in . By Theorem 3.12 we may pass to a further subsequence (not relabeled) such that generates a Young measure on . Since is supported on the set of quenched -two-scale cluster points of , we deduce from Lemma 3.10 that the support of is contained in which is a closed subspace of . Moreover, thanks to the relation of the generated Young measure and stochastic two-scale convergence in the mean, we have . Furthermore, Lemma 3.14 implies that

    In view of (32) and the fact that is supported in , we conclude that

    Since , we have , and thus we conclude that for -a.a. , concentrates on .

    Step 2. (The strictly convex case).

    The uniqueness of and is clear. From Step 1 we thus conclude that where . Theorem 3.12 implies that (for -a.a. ). By Lemma 3.14 we have for -a.a. ,

    On the other hand, since minimizes , we deduce by a standard argument that for -a.a. ,

    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] W. Chen, L. Xu, S. P. Zhu, Stock loan valuation under a stochstic interest rate model, Comput. Math. Appl., 70 (2015), 1757–1771. https://doi.org/10.1016/j.camwa.2015.07.019 doi: 10.1016/j.camwa.2015.07.019
    [2] J. Xia, X. Y. Zhou, Stock loans, Math. Financ., 17 (2007), 307–317. https://doi.org/10.1111/j.1467-9965.2006.00305.x doi: 10.1111/j.1467-9965.2006.00305.x
    [3] Z. Liang, W. Wu, S. Jiang, Stock loan with automatic termination clause, cap and margin, Comput. Math. Appl., 60 (2010), 3160–3176. https://doi.org/10.1016/j.camwa.2010.10.021 doi: 10.1016/j.camwa.2010.10.021
    [4] N. Cai, W. Zhang, Regime classification and stock loan valuation, Oper. Res., 68 (2020), 965–983. https://doi.org/10.1287/opre.2019.1934 doi: 10.1287/opre.2019.1934
    [5] M. Dai, Z. Q. Xu, Optimal redeeming strategy of stock loans with finite maturity, Math. Financ., 21 (2011), 775–793. https://doi.org/10.1111/j.1467-9965.2010.00449.x doi: 10.1111/j.1467-9965.2010.00449.x
    [6] X. Lu, E. R. M. Putri, Semi-analytic valuation of stock loans with finite maturity, Commun. Nonlinear. Sci., 27 (2015), 206–215. https://doi.org/10.1016/j.cnsns.2015.03.007 doi: 10.1016/j.cnsns.2015.03.007
    [7] D. Prager, Q. Zhang, Stock loan valuation under a regime-switching model with mean-reverting and finite maturity, J. Syst. Sci. Complex., 23 (2010), 572–583. https://doi.org/10.1007/s11424-010-0146-7 doi: 10.1007/s11424-010-0146-7
    [8] T. W. Wong, H. Y. Wong, Stochastic volatility asymptotics of stock loans: valuation and optimal stopping, J. Math. Anal. Appl., 394 (2012), 337–346. https://doi.org/10.1016/j.jmaa.2012.04.067 doi: 10.1016/j.jmaa.2012.04.067
    [9] Z. Zhou, J. Ma, X. Gao, Convergence of iterative Laplace transform methods for a system of fractional PDEs and PIDEs arising in option pricing, E. Asian. J. Appl. Math., 8 (2018), 782–808. https://doi.org/10.4208/eajam.130218.290618 doi: 10.4208/eajam.130218.290618
    [10] W. Chen, X. Xu, S. P. Zhu, Analytically pricing double barrier options based on a time-fractional Black-Scholes equation, Comput. Math. Appl., 69 (2015), 1407–1419. https://doi.org/10.1016/j.camwa.2015.03.025 doi: 10.1016/j.camwa.2015.03.025
    [11] P. Carr, H. Geman, D. B. Madan, Y. Marc, The fine structure of asset returns: an empirical investigation, J. Bus., 75 (2002), 305–332. https://doi.org/10.1086/338705 doi: 10.1086/338705
    [12] A. Cartea, D. del-Castillo-Negrete, Fractional diffusion models of option prices in markets with jumps, Physica A, 374 (2007), 749–763. https://doi.org/10.1016/j.physa.2006.08.071 doi: 10.1016/j.physa.2006.08.071
    [13] W. Chen, X. Xu, S. P. Zhu, Analytically pricing European-style options under the modified Black-Scholes equation with a spatial-fractional derivative, Q. Appl. Math., 72 (2014), 597–611. https://doi.org/10.1090/S0033-569X-2014-01373-2 doi: 10.1090/S0033-569X-2014-01373-2
    [14] W. Chen, M. Du, X. Xu, An explicit closed-form analytical solution for European options under the CGMY model, Commun. Nonlinear. Sci., 42 (2017), 285–297. https://doi.org/10.1016/j.cnsns.2016.05.026 doi: 10.1016/j.cnsns.2016.05.026
    [15] W. Chen, S. Lin, Option pricing under the KoBol model, Anziam. J., 60 (2018), 175–190. https://doi.org/10.1017/S1446181118000196 doi: 10.1017/S1446181118000196
    [16] W. Chen, X. Xu, S. P. Zhu, A predictor-corrector approach for pricing American options under the finite moment log-stable model, Appl. Numer. Math., 97 (2015), 15–29. https://doi.org/10.1016/j.apnum.2015.06.004 doi: 10.1016/j.apnum.2015.06.004
    [17] X. M. Gu, T. Z. Huang, C. C. Ji, B. Carpentieri, A. A. Alikhanov, Fast iterative method with a second-order implicit difference scheme for time-space fractional convection-diffusion equation. J. Sci. Comput., 72 (2017): 957–985. https://doi.org/10.1007/s10915-017-0388-9 doi: 10.1007/s10915-017-0388-9
    [18] X. M. Gu, Y. L. Zhao, X. L. Zhao, B. Carpentieri, Y. Y. Huang, A note on parallel preconditioning for the all-at-once solution of Riesz fractional diffusion equations, Numer. Math-Theory. Me., 14 (2021): 893–919. https://doi.org/10.48550/arXiv.2003.07020 doi: 10.48550/arXiv.2003.07020
    [19] X. Chen, D. Ding, S. L. Lei, W. Wang, A fast preconditioned iterative method for two-dimensional options pricing under fractional differential models. Comput. Math. Appl., 79 (2020): 440–456. https://doi.org/10.1016/j.camwa.2019.07.010 doi: 10.1016/j.camwa.2019.07.010
    [20] S. L. Lei, W. Wang, X. Chen, D. Ding, A fast preconditioned penalty method for American options pricing under regime-switching tempered fractional diffusion models, J. Sci. Comput., 75 (2018): 1633–1655. https://doi.org/10.1007/s10915-017-0602-9 doi: 10.1007/s10915-017-0602-9
    [21] W. Wang, X. Chen, D. Ding, S. L. Lei, Circulant preconditioning technique for barrier options pricing under fractional diffusion models, Int. J. Comput. Math., 92 (2015): 2596–2614. http://dx.doi.org/10.1080/00207160.2015.1077948 doi: 10.1080/00207160.2015.1077948
    [22] X. Chen, W. Wang, D. Ding, S. L. Lei, A fast preconditioned policy iteration method for solving the tempered fractional HJB equation governing American options valuation, Comput. Math. Appl., 73 (2017): 1932–1944. http://dx.doi.org/10.1016/j.camwa.2017.02.040 doi: 10.1016/j.camwa.2017.02.040
    [23] S. G. Kou, A jump-diffusion model for option pricing, Manage. Sci., 48 (2002), 1086–1101. https://doi.org/10.1287/mnsc.48.8.1086.166 doi: 10.1287/mnsc.48.8.1086.166
    [24] B. F. Nielsen, O. Skavhaug, A. Tveito, Penalty and front-fixing methods for the numerical solution of American option problems, J. Comput. Financ., 4 (2002), 69–97. https://doi.org/10.21314/jcf.2002.084 doi: 10.21314/jcf.2002.084
    [25] S. Chen, F. Liu, X. Jiang, I. Turner, V. Anh, A fast semi-implicit difference method for a nonlinear two-sided space-fractional diffusion equation with variable diffusivity coefficients, Appl. Math. Comput., 257 (2015), 591–601. https://doi.org/10.1016/j.amc.2014.08.031 doi: 10.1016/j.amc.2014.08.031
    [26] H. Hejazi, T. Moroney, F. Liu, Stability and convergence of a finite volume method for the space fractional advection-dispersion equation, J. Comput. Appl. Math., 255 (2014), 684–697. https://doi.org/10.1016/j.cam.2013.06.039 doi: 10.1016/j.cam.2013.06.039
    [27] O. Marom, E. Momoniat, A comparison of numerical solutions of fractional diffusion models in finance, Nonlinear. Anal-Real., 10 (2009), 3435–3442. https://doi.org/10.1016/j.nonrwa.2008.10.066 doi: 10.1016/j.nonrwa.2008.10.066
    [28] N. Cai, S. G. Kou, Option pricing under a mixed-exponential jump diffusion model, Manage. Sci., 57 (2011), 2067–2081. https://doi.org/10.1287/mnsc.1110.1393 doi: 10.1287/mnsc.1110.1393
    [29] Y. L. Zhao, P. Y. Zhu, X. M. Gu, X. L. Zhao, H. Y. Jian, A preconditioning technique for all-at-once system from the nonlinear tempered fractional diffusion equation, J. Sci. Comput., 83 (2020): 10. https://doi.org/10.1007/s10915-020-01193-1 doi: 10.1007/s10915-020-01193-1
    [30] S. L. Lei, H. W. Sun, A circulant preconditioner for fractional diffusion equations, J. Comput. Phys., 242 (2013), 715–725. https://doi.org/10.1016/j.jcp.2013.02.025 doi: 10.1016/j.jcp.2013.02.025
  • 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
  • © 2023 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(1871) PDF downloads(60) Cited by(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog