Random attractors via pathwise mild solutions for stochastic parabolic evolution equations

We investigate the longtime behavior of stochastic partial differential equations (SPDEs) with differential operators that depend on time and the underlying probability space. In particular, we consider stochastic parabolic evolution problems in Banach spaces with additive noise and prove the existence of random exponential attractors. These are compact random sets of finite fractal dimension that contain the global random attractor and are attracting at an exponential rate. In order to apply the framework of random dynamical systems, we use the concept of pathwise mild solutions.


Introduction
Our aim is to study the longtime dynamics of stochastic evolution equations using an approach that is different from the classical one. Namely, instead of transforming the SPDE into a random PDE, we work with solutions that are defined pathwise, see [5,26]. We consider parabolic problems with random differential operators and use a pathwise representation formula to show that the solution operator generates a random dynamical system and to prove that it possesses random attractors of finite fractal dimension.
In particular, let X be a separable Banach space and let ( , F, P) be a complete probability space with filtration (F t ) t∈R . We consider stochastic parabolic evolution equations of the form du(t) = (A(t, ω)u(t) + F(u(t)))dt + σ dW (t), (1.1) random coefficients. The resulting random PDEs can be studied by deterministic techniques. This method has been applied to a large variety of PDEs, mainly for equations perturbed by additive noise or by a particular linear multiplicative noise, e.g., see [4,16,17,28,31] and the references therein. However, for more general situations such a change of variables is not always known or cannot be performed. In [15], using the theory of strongly monotone operators, a strictly stationary solution of the equation was constructed. This allows to transform SPDEs of the form (1.1) into a family of random PDEs. Using this ansatz, the existence of random attractors was shown in [15] for a class of SPDEs including equations such as (1.1). Here, we follow a different approach. We aim to use the notion of pathwise mild solutions introduced by Pronk and Veraar in [26] to establish the existence of global and exponential random attractors for (1.1). So far, only few results concerning the existence of exponential attractors for SPDE have been obtained. In particular, let (U (t, s, ω)) t≥s,ω∈ be the stochastic evolution system generated by the family (A(t, ω)) t∈R,ω∈ , then the pathwise mild solution of (1.1) is defined as where, for simplicity, we omit to write the dependency of A and U on ω. This formula is motivated by formally applying integration by parts for the stochastic integral, and it, indeed, yields a pathwise representation for the solution, see [26]. Instead, if one directly used the classical mild formulation of SPDEs to define a solution, the resulting stochastic integral would not be well defined (see Sect. 2

.3).
Our aim is to show that problem (1.1) generates a random dynamical system using the concept of pathwise mild solutions and to prove the existence of random attractors. We will not only consider (global) random attractors, but also show that random exponential attractors exist. In particular, the existence of random exponential attractors immediately implies the existence and finite fractal dimension of the (global) random attractor. To this end, we employ a general existence result for random exponential attractors in [7] which turns out to be easily applicable in our setting.
Stochasticity plays an important role in many real world applications. Complex systems in physics, engineering or biology can be described by PDEs with coefficients that depend on stochastic processes. These random terms quantify the lack of knowledge of certain parameters in the equation or reflect external fluctuations. Problem (1.1) is a semilinear parabolic problem where the coefficients of the differential operators (A(t, ω)) t∈R,ω∈ depend on a stochastic process with suitable properties, and the equation is perturbed by additive noise. A related but simpler setting is random parabolic equations of the form du(t) = (A(t, ω)u(t) + F(t, ω, u(t)))dt. (1. 2) The longtime behavior of such random evolution equations has been investigated using the random dynamical system approach in [6,20,21,27]. To this end, the following structure of the random generators was assumed, where ( , F, P, (θ t ) t∈R ) is an ergodic metric dynamical system. In this context, results concerning invariant manifolds [6,20], principal Lyapunov exponents [21] and the stability of equilibria [27] have been obtained. Random evolution equations of the form (1.2) arise in several applications. For instance, setting A(θ t ω)u := u+a(θ t ω)u and F(t, ω, u) := −a(θ t ω)u 2 , with a suitable measurable function a : → (0, ∞), we recover a random version of the Fisher-KPP equation, which was analyzed in [27].
In the present work, we consider equations of the form (1.1), i.e., we perturb a semilinear nonautonomous random parabolic equation by an infinite-dimensional Brownian motion and investigate the existence of random attractors.
The outline of our paper is as follows. In Sect. 2, we collect basic notions and results from the theory of random dynamical systems and nonautonomous stochastic evolution equations and recall an existence result for random exponential attractors. In Sect. 3, we formulate and prove our main results. First, we show that under suitable assumptions on A, F and W , the solution operator corresponding to (1.1) generates a random dynamical system. Then, we establish the existence of an absorbing set and verify the so-called smoothing property of the random dynamical system. These properties allow us to conclude the existence of random exponential attractors in Theorem 3.8 and to derive upper bounds for their fractal dimension. As a consequence, the (global) random attractor exists and its fractal dimension is finite. In Sect. 4, we provide explicit examples of nonautonomous random differential operators satisfying our hypotheses and point out potential applications of our results.
Our paper provides a first, simple example that illustrates how the concept of pathwise mild solutions can be used to show the existence of global and exponential random attractors for SPDEs with random differential operators. Numerous extensions are imaginable. In particular, in future works we plan to relax the assumptions on the nonlinear term F and to consider Problem (1.1) with multiplicative noise. Another interesting aspect would be to investigate the regularity of random attractors.

Preliminaries
We first collect some basic notions and results from the theory of random dynamical systems, which are mainly taken from [3,10,28]. Then, in Sect. 2.2, we state a general existence theorem for random exponential attractors which was proven in [7]. In Sect. 2.3, we recall the notion of pathwise mild solutions for stochastic evolution equations introduced in [26].

Random dynamical systems and random attractors
In order to quantify uncertainty, we describe an appropriate model of the noise. If not further specified, ( , F, P) denotes a probability space. Moreover, X is a separable and reflexive Banach space and · X denotes the norm in X .
be a family of P-preserving transformations (meaning that θ t P = P for all t ∈ R) with the following properties: where B(R) denotes the Borel σ -algebra. Then, the quadruple ( , F, P, (θ t ) t∈R ) is called a metric dynamical system.
(b) We always assume that P is ergodic with respect to (θ t ) t∈R , i.e., any (θ t ) t∈Rinvariant subset has either zero or full measure.
Our aim is to introduce a metric dynamical system associated with a two-sided X -valued Wiener process.
For the sake of completeness, we recall the construction of such a process if X is not a Hilbert space. First, we introduce an auxiliary separable Hilbert space H and denote by (W H (t)) t≥0 an H -cylindrical Brownian motion, i.e., ( where (γ n ) n∈N is a sequence of independent standard Gaussian random variables and (e n ) n∈N is an orthonormal basis in H . If X is isomorphic to H , then the previous condition means that G is a Hilbert-Schmidt operator (notation: G ∈ L 2 (H )). In this framework, letting (ẽ n ) n∈N be an orthonormal basis of (kerG) ⊥ , we know according to [30,Prop. 8.8] that the series converges almost surely and defines an X -valued Brownian motion. Its covariance operator is given by tGG * , where G * denotes the adjoint. Moreover, every X -valued Brownian motion can be obtained in this way. Again, if X is isomorphic to H and G ∈ L 2 (H ) with G L 2 (H ) = Tr(GG * ), then the previous definition entails a traceclass Wiener process. Finally, we extend this to a two-sided process in the standard way.
To obtain a metric dynamical system associated with such a process, we let C 0 (R; X ) denote the set of continuous X -valued functions which are zero at zero equipped with the compact open topology. We take P as the Wiener measure on B(C 0 (R; X )) having a covariance operator Q on X . Then, Kolmogorov's theorem about the existence of a continuous version yields the canonical probability space (C 0 (R; X ), B(C 0 (R; X )), P). Moreover, to obtain an ergodic metric dynamical system we introduce the Wiener shift, which is defined as (2.1) Throughout this manuscript, θ t ω(·) will always denote the Wiener shift. We now recall the definition of a random dynamical system.
The second property is referred to as the cocycle property and generalizes the semigroup property. In fact, if ϕ is independent of ω, (ii) reduces exactly to the semigroup property, i.e., ϕ(t + τ, x) = ϕ(t, ϕ(τ, x)). For random dynamical systems, the evolution of the noise (θ t ω) has additionally to be taken into account.
Under suitable assumptions, the solution operator of a random differential equation generates a random dynamical system. Stochastic (partial) differential equations are more involved since stochastic integrals are defined almost surely, though the cocycle property must hold for all ω.
Referring to the monograph by Arnold [3], it is well known that stochastic (ordinary) differential equations generate random dynamical systems under suitable assumptions on the coefficients. This is due to the flow property, see [19], which can be deduced from Kolmogorov's theorem about the existence of a (Hölder-) continuous random field with a finite-dimensional parameter range. Here, the parameters of this random field are the time and the nonrandom initial data.
Whether an SPDE generates a random dynamical system has been a long-standing open problem, since Kolmogorov's theorem breaks down for random fields parameterized by infinite-dimensional Hilbert spaces, see [23]. As a consequence, the question of how a random dynamical system can be obtained from an SPDE is not trivial, since solutions are only defined almost surely which is insufficient for the cocycle property. In particular, there exist exceptional sets which depend on the initial condition, and if more than countably many exceptional sets occur, it is unclear how the random dynamical system can be defined. This problem was fully solved only under restrictive assumptions on the structure of the noise. More precisely, for SPDEs with additive or linear multiplicative noise, there are standard transformations which reduce these SPDEs in PDEs with random coefficients. Since random PDEs can be solved pathwise, the generation of the random dynamical system is straightforward.
Before we recall the notions of global and exponential random attractors, we need to introduce the class of tempered random sets. From now on, in this Sect. 2.2, when stating properties involving a random parameter, we assume, unless otherwise specified, that they hold on a (θ t ) t∈R -invariant subset of of full measure, i.e., there exists a (θ t ) t∈R -invariant subset 0 ⊂ of full measure such that the property holds for all ω ∈ 0 . To simplify notations, we denote 0 again by .
The random set B is bounded (or compact) if the sets B(ω) ⊂ X are bounded (or compact) for all ω ∈ .
Here and in the sequel, we denote by D the collection of tempered random sets in X . Definition 2.6. Let ϕ be a random dynamical system on X . A random set {A(ω)} ω∈ ∈ D is called a D-random (pullback) attractor for ϕ if the following properties are satisfied: The following theorem provides a criterion for the existence of random attractors, see Theorem 4 in [11]. The uniqueness follows from Corollary 1 in [11].
Theorem 2.7. There exists a D-random (pullback) attractor for ϕ if and only if there exists a compact random set that pullback attracts all random sets D ∈ D. Moreover, the random (pullback) attractor is unique.
One way of proving the existence of the random attractor, that in addition implies its finite fractal dimension, is to show that a random exponential attractor exists. Exponential attractors are compact subsets of finite fractal dimension that contain the global attractor and are attracting at an exponential rate. This notion was first introduced for semigroups in the autonomous deterministic setting [12] and has later been extended for nonautonomous and random dynamical systems, see [7,9] and the references therein.
Here, we consider so-called nonautonomous random exponential attractors, see [7]. While random exponential attractors in the strict sense are positively ϕ-invariant, nonautonomous random exponential attractors are positively ϕ-invariant in the weaker, nonautonomous sense. To construct exponential attractors for time-continuous random dynamical systems that are positively ϕ-invariant typically requires the Hölder continuity in time of the cocycle which is a restrictive assumption. However, if we relax the invariance property and consider nonautonomous random exponential attractors instead, only the Lipschitz continuity of the cocycle in space is needed. In fact, the construction can be essentially simplified, we obtain better bounds for the fractal dimension and the assumption of Hölder continuity in time can be omitted, see [7]. Even though we could prove the Hölder continuity in time for the cocycle for our particular problem, we omit it since it has no added value for our main results and would lead to weaker bounds for the fractal dimension. Definition 2.8. A nonautonomous tempered random set {M(t, ω)} t∈R,ω∈ is called a nonautonomous D-random (pullback) exponential attractor for ϕ if there existst > 0 such that M(t +t, ω) = M(t, ω) for all t ∈ R, ω ∈ , and the following properties are satisfied: (c) {M(t, ω)} t∈R,ω∈ is pullback D-attracting at an exponential rate, i.e., there exists α > 0 such that We recall that the fractal dimension of a precompact subset M ⊂ X is defined as where N ε (M) denotes the minimal number of ε-balls in X with centers in M needed to cover the set M. By Theorem 2.7, the existence of a nonautonomous random exponential attractor immediately implies that the (global) random attractor exists. Moreover, the global random attractor is contained in the random exponential attractor, and hence, its fractal dimension is finite.
Existence proofs for global and exponential random attractors are typically based on the existence of a pullback D-absorbing set for ϕ.
The following condition is convenient to show the existence of an absorbing set. Namely, if for every x ∈ D(θ −t ω), D ∈ D and ω ∈ , it holds that where ρ(ω) > 0 for every ω ∈ , then the ball B(ω) := B(0, ρ(ω) + δ) centered in 0 with radius ρ(ω) + δ for some constant δ > 0 is a random absorbing set. For further details and applications, see [5,28]. Instead of considering random exponential attractors which is typically more involved and requires to verify additional properties of the cocycle, the existence of random attractors is frequently shown using the following result, see Theorem 2.1 in [28]. Theorem 2.10. Let ϕ be a continuous random dynamical system on X over ( , F, P, (θ t ) t∈R ). Suppose that {B(ω)} ω∈ is a compact random absorbing set for ϕ in D. Then, ϕ has a unique D-random attractor {A(ω)} ω∈ which is given by We could apply Theorem 2.10 to prove the existence of a random attractor for our particular problem. However, showing that a nonautonomous random exponential attractor exists does not only imply the existence of the random attractor, but also its finite fractal dimension. Moreover, it turns out to be even simpler in our case than applying Theorem 2.10. To this end, we use an existence result for random exponential attractors obtained in [7] that we recall in the next subsection.

An existence result for random exponential attractors
The existence result for random pullback exponential attractors is based on an auxiliary normed space that is compactly embedded into the phase space and the entropy properties of this embedding. We recall some notions and results that we will need in the sequel, see also [7][8][9].
The (Kolmogorov) ε-entropy of a precompact subset M of a Banach space X is defined as where N X ε (M) denotes the minimal number of ε-balls in X with centers in M needed to cover the set M. It was first introduced by Kolmogorov and Tihomirov [14]. The order of growth of H X ε (M) as ε tends to zero is a measure for the massiveness of the set M in X , even if its fractal dimension is infinite.
If X and Y are Banach spaces such that the embedding Y → X is compact, we use the notation where B Y (0, 1) denotes the closed unit ball in Y .
Remark 2.11. The ε-entropy is related to the entropy numbersê k for the embedding Y → X, which are defined bŷ If the embedding is compact, thenê k is finite for all k ∈ N. For certain function spaces, the entropy numbers can explicitly be estimated (see [13]). For instance, if D ⊂ R n is a smooth bounded domain, then the embedding of the Sobolev spaces is compact if l 1 > l 2 and l 1 n − 1 p 1 > l 2 n − 1 p 2 . Moreover, the entropy numbers grow polynomially, namelyê k k − l 1 −l 2 n (see Theorem 2, Section 3.3.3 in [13]), and consequently, for some constant c > 0. Here, we write f g, if there exist positive constants c 1 and c 2 such that The following existence result for nonautonomous random pullback exponential attractors is a special case of the main result in [7]. In fact, we formulate a simplified version that suffices for the parabolic stochastic evolution problem we consider. In particular, we assume that the cocycle is uniformly Lipschitz continuous and satisfies the smoothing property with a constant that is independent of ω. More generally, one can allow that the constants depend on the random parameter ω and that the cocycle is asymptotically compact, i.e., it is the sum of a mapping satisfying the smoothing property and a contraction. Theorem 2.12. Let ϕ be a random dynamical system in a separable Banach space X , and let D denote the universe of tempered random sets. Moreover, we assume that the following properties hold for all ω ∈ : (H 1 ) Compact embedding: There exists another separable Banach space Y that is compactly and densely embedded into X .

Pathwise mild solutions for parabolic SPDEs
Let := {(s, t) ∈ R 2 : s ≤ t}, X be a separable, reflexive, type 2 Banach space and ( , F , P) be a probability space. Similarly to [26], we consider nonautonomous SPDEs of the form where A = {A(t, ω)} t∈R,ω∈ is a family of time-dependent random differential operators. Intuitively, this means that the differential operator depends on a stochastic processes, in a meaningful way which will be specified later. We aim to investigate the longtime behavior of (2.3) using a random dynamical systems approach. First, we recall sufficient conditions that ensure that the family A generates a parabolic stochastic evolution system, see [26]. In particular, we make the following assumptions concerning measurability, sectoriality and Hölder continuity of the operators.  A(t, ω)) for all t ∈ R, ω ∈ . (A1) The mapping A : R × → L(D A , X ) is strongly measurable and adapted. (A2) There exists ϑ ∈ (π, π 2 ) and M > 0 such that ϑ := {μ ∈: |arg μ| < ϑ} ⊂ ρ(A(t, ω)) and (A3) There exists ν ∈ (0, 1] and a mapping C : → X such that where we assume that C(ω) is uniformly bounded with respect to ω, see [26].
Assumptions (A2) and (A3) are referred to in the literature as the Kato-Tanabe assumptions, compare [2], p. 55, or [24], p. 150, and are common in the context of nonautonomous evolution equations. Since the constants in (A2) and (A3) are uniformly bounded w.r.t. ω, all constants arising in the estimates below do not dependent on ω.

Assumption 2. (AC) We assume that the operators
(2.5) (Drift) The nonlinearity F : X → X is globally Lipschitz continuous, i.e., there exists a constant C F > 0 such that This implies a linear growth condition on F. Namely, there exist a positive constant C F such that Definition 2.13. Let X 1 and X 2 be two separable Banach spaces. A random operator L : × X 1 → X 2 is called strongly measurable if the mapping ω → L(ω)x,ω ∈ , is a random variable on X 2 for every x ∈ X 1 . Theorem 2.14. There exists a unique parabolic evolution system U : × → L(X ) with the following properties: for all 0 ≤ r ≤ s ≤ t, ω ∈ . (5) The evolution system U : × → L(X ) is strongly measurable in the uniform operator topology. Moreover, for every t ≥ s, the mapping ω → U (t, s, ω) ∈ L(X ) is strongly F t -measurable in the uniform operator topology.
To prove the existence of random attractors, we need additional smoothing properties of the evolution system. The following properties and estimates were shown in Lemmas 2.6 and 2.7 in [26]. The exponential decay is a consequence of our assumption (U ).

Lemma 2.15. We assume that the family of adjoint operators A * (t, ω) satisfies (A3)
with exponent ν * > 0. Then, for every t > 0, the mapping s → U (t, s, ω) belongs to C 1 ([0, t); L(X )), and for all x ∈ D A one has Moreover, for α ∈ [0, 1] and η ∈ (0, 1) there exist positive constants C α , C α,η such that the following estimates hold for t > s andω ∈ : To shorten notations, in the sequel we omit the ω-dependence of A and U if there is no danger of confusion. The classical mild formulation of the SPDE (2.3) is The formula is motivated by formally applying integration by parts for the stochastic integral, and, as shown in [26], it indeed yields a pathwise representation for the solution.
Our aim is to show the existence of random attractors for SPDEs using this concept of pathwise mild solutions. It allows us to study random attractors without transforming the SPDE into a random PDE, as it is typically done.
Remark 2. 16. We emphasize that the concept of pathwise mild solutions also applies if σ is not constant, see [26,Sec.5]. In this case, the solution of (2.3) is given by (2.11) However, it is not possible to obtain a random dynamical system in this case, due to the presence of the stochastic integrals in (2.10) and (2.11) which are not defined in a pathwise sense. Consequently, this representation formula does not hold for every ω ∈ . We aim to investigate this issue in a future work.
Recalling that W is an X β -valued Wiener process, we introduce the canonical probability space and identify W (t, ω) =: ω(t), for ω ∈ . Moreover, together with the Wiener shift, we obtain, analogously as in Sect. 2.1, the ergodic metric dynamical system ( , F, P, (θ t ) t∈R ).
In the following, ( , F, P) always denotes the probability space (2.12).

Random dynamical system and absorbing set
Since we consider SPDEs with time-dependent random differential operators, we need to impose additional structural assumptions in order to use the framework of random dynamical systems, see [6,20].

Assumption 3.
(RDS) We assume that the generators depend on t and ω in the following way: A(t, ω) = A(θ t ω) for all t ∈ R and ω ∈ . (3.1) This assumption is needed to obtain the cocycle property. In this case, the group property of the metric dynamical system implies hat Moreover, one can easily show that A(θ t ω) generates a random dynamical system, i.e., the solution operator corresponding to the linear evolution equation forms a random dynamical system.

Theorem 3.1. Let U : × → L(X ) be the evolution operator generated by
is a random dynamical system.
Intuitively, this means that starting at time s on the ω-fiber of the noise and letting time t > 0 pass lead to the same state as starting at time zero on the shifted θ s -fiber of the noise and letting the system evolve for time t. At the level of random generators, U (t + s, s, ω) is obtained from A(θ t ω) which is the same as A(θ s θ t−s ω) due to the properties of the metric dynamical system. Therefore, the cocycle property is satisfied. The measurability ofŨ follows from Theorem 2.14, Property 5.

Remark 3.2.
We obtain the measurability of U directly from the results in [26]. Alternatively, one can show the measurability ofŨ as in [20,Lem. 14] using Yosida approximations of A(ω). Here, one argues that the evolution operators corresponding to these approximations are strongly measurable and then pass to the limit. These arguments exploit the structural assumption (3.1). The proof of the measurability in [26] is more involved and holds under more general assumptions.
We give a standard example of a random nonautonomous generator and its corresponding evolution operator. Example 3.3. A simple example for an operator that satisfies our assumptions is a random perturbation of a uniformly elliptic operator A (in a smooth bounded domain with homogeneous Dirichlet boundary conditions) by a real-valued Ornstein-Uhlenbeck process, which is the stationary solution of the Langevin equation Here, μ > 0 and W is a two-sided real-valued Brownian motion. We denote by ( , F , P) its associated canonical probability space and make the identification W (t, ω) := ω(t) for ω ∈ . Then, we have that In this case, the parabolic evolution operator generated by A + z(θ t ω) is where (T (t)) t≥0 is the analytic C 0 -semigroup generated by A. We have This simple example illustrates that the formalism we introduced above is meaningful. Further examples of random time-dependent generators are provided in Sect. 4. For additional applications, we refer to [26], and to [6,20] in the context of random dynamical systems.
We now prove the existence of random attractors for SPDEs of the form using pathwise mild solutions as defined in (2.9).

Remark 3.4.
• We emphasize that the SPDE (3.4) cannot be transformed into a PDE with random coefficients using the stationary Ornstein-Uhlenbeck process, since the convolution is not defined and one has to make sense of it using the integration by parts formula (3.5) A different transformation based on a strictly stationary solution has been introduced in [15].
• Another approach would be to subtract the noise, i.e., to introduce the change of variables v := u − σ ω. This would formally lead to the random PDE The mild solution of (3.6) would be given by However, we would need to justify that this mild solution is well defined, and the noise also interacts with the nonlinear term. In order to avoid these difficulties, we work with pathwise mild solutions.
Using (3.5), the representation formula of a solution for (3.4) reads as Here, we used in the last line that t, s, ω)A(s, ω).
Remark 3.5. We emphasize that the pathwise mild solution concept is applicable also under weaker assumptions on the noise, for instance, if ω takes values in a suitable extrapolation space [25, Section 3.1]. Moreover, the formal computations made in (3.8) can be justified even if ω ∈ D A . In fact, according to [26,Theorem 4.9] we know that the pathwise mild solution is equivalent to the weak solution of (3.4). For simplicity, we test the linear part (i.e., for F ≡ 0) of (3.4) with x * ∈ D A * := D((A * (t))). This yields where ·, · denotes the dual pairing. Plugging the identity which holds for ω(·) ∈ X (see [26,Section 4.4]), into (3.9) entails Lemma 3.6. The solution operator corresponding to (3.4) generates a random dynamical system ϕ : Proof. We only verify the cocycle property. The continuity is straightforward, and the measurability of ϕ follows from the measurability ofŨ . Let s, t ≥ 0. Using (3.8), we have Using that one immediately gets θ s ω, ϕ(s, ω, u 0 )).
This proves the statement.
Next, we show the existence of an absorbing set. We recall that λ > cC F as assumed in (Drift). Here, λ, c and C F are the constants in (2.5) and (2.6).
From now on, the properties and statements hold for all ω ∈ 0 where 0 ⊂ is the set of all ω that have subexponential growth. The set 0 is (θ t ) t∈R -invariant and has full measure, see, e.g., [5]. To simplify notations, in the sequel, we will denote 0 again by . Lemma 3.7. The random dynamical system ϕ has a pullback absorbing set.
Proof. We verify (2.2). To this end, we use the estimates in Lemma 2.15 and Gronwall's lemma. We observe that ω ∈ X β implies that for some constantĉ > 0.
It is convenient to work with the representation formula (3.7). Assuming that the fiber is given for θ −t ω, we obtain where C 1−β denotes the constant in Lemma 2.15. Hence, combining the estimates we obtain we can rewrite the previous inequality as (3.10) Applying Gronwall's lemma to e λt u(t) X , we infer that e λt u(t) X ≤ e λt γ (t) + cC F t 0 e λs γ (s)e cC F (t−s) ds, and multiplying with e −λt , we obtain This estimate allows us to determine the pullback absorbing set. First, note that all terms in γ are well defined for the limit t → ∞, due to the subexponential growth of ω(t) X β , and consequently, (3.12) We now focus on the second term in (3.11), The first and the third terms are bounded by Vol. 21 (2021) Random attractors via pathwise mild solutions 2653 The second one can be estimated by Finally, for the last term, we observe that In conclusion, using all the previous estimates in (3.11) we have is a pullback absorbing set for our random dynamical system. This expression is natural, since we can immediately see the influence of the linear part, nonlinear term and noise intensity. The previous integrals are well defined due to the subexponential growth of ω. More precisely, the set of all ω ∈ that have subexponential growth is invariant and has full measure. The temperedness of the absorbing set can be verified as in [5,Lem. 3.7].

Existence and finite fractal dimension of random attractors
We now apply Theorem 2.12 to deduce the existence of nonautonomous random exponential attractors for the random dynamical system ϕ.
We remark that κ is determined by the constant C η in Lemma 2.15, the Lipschitz constant C F of F and the constants c and λ in (2.5).
Proof. We verify the hypotheses in Theorem 2.12.
(H 1 ) By Assumptions 2 (AC), this property holds for the spaces X and Y = X η , for arbitrary η ∈ (0, 1). (H 2 ) This was shown in Lemma 3.7. In fact, B(ω) = B(0, ρ(ω) + δ), for some δ > 0, is pullback D-absorbing and B ∈ D. Moreover, the absorbing time fulfills the condition in (H 1 ). In fact, let D ∈ D, then Hence, Gronwall's lemma implies that (H 3 ) Finally, we use the Lipschitz continuity in (H 4 ) to verify the smoothing property for the spaces X and Y = X η . Lett > 0. We estimate the difference of the solutions in the X η -norm, where we used the Lipschitz continuity of ϕ in X in the last step. Hence, the smoothing property holds with the smoothing constant The smoothing property holds for anyt > 0, and consequently, (H 3 ) is satisfied.
An immediate consequence is the existence and finite fractal dimension of the global random attractor. Corollary 3.9. There exists a unique global random attractor for ϕ, and its fractal dimension is bounded by for all η ∈ (0, 1), where κ is the constant given in Theorem 3.8.
The existence of a nonautonomous random exponential attractor implies that the global random attractor exists and that its fractal dimension is finite. We point out that in our particular case, it is, in fact, easier to consider random exponential attractors than to deduce the existence of the global random attractor from Theorem 2.10. In fact, we have shown that a tempered absorbing set exists, but the theorem requires the existence of a compact absorbing set. To show how this can be established and to indicate that the proof is indeed more involved than verifying the hypotheses of Theorem 2.12, we provide the following lemma, even if we do not use it to prove our main results. Proof. The proof is based on compact embeddings of fractional power spaces. Let η > 0 be such that 0 < η < β ≤ 1. It suffices to derive uniform estimates of the solutions w.r.t. the X η -norm since X η is compactly embedded into X . Let u 0 ∈ B(θ −T B ω). We observe that To estimate these terms, we use that u 0 ∈ B(θ −T B ω) and that B is a pullback absorbing set. The first and second terms yield the following expressions, For the generalized convolution, we obtain Finally, we estimate the drift term, where we used that u 0 ∈ B(θ −T B ω) and the absorbing property of B, i.e., ϕ(r + T B , θ −T B ω, u 0 ) ⊂ B(θ r ω).

Examples
We provide examples of differential operators A that satisfy our assumptions in the previous sections. The canonical examples are uniformly elliptic operators with random, time-dependent coefficients. Such operators have been investigated in the context of SPDEs, see [25,26] and the references therein. In the framework of random dynamical systems, several properties of the evolution system generated by such operators have been analyzed, including results on the spectral theory and principal Lyapunov exponents [21,22,29], stable and unstable manifolds and multiplicative ergodic theorems [6,20].
We consider the random partial differential operators in the Banach space X := L p (G) for 2 ≤ p < ∞, where G ⊂ R n is a bounded open domain with smooth boundary ∂G (see also [24,Sec. 7.6]). We recall that in our case, the differential operator A(t, ω) depends on time t ∈ R and the random parameter ω ∈ in the following way Example 4.1. Let m ∈ N and A be the random partial differential operator with homogeneous Dirichlet boundary conditions, D k u = 0 on ∂G for |k| < m.
We make the following assumptions on the coefficients.
1. The operator A is uniformly strongly elliptic in G, i.e., there exists a constant c > 0 such that a k (θ t ω, x)ξ k ≥ c|ξ | 2m for all t ∈ R, x ∈ G, ξ ∈ R n .
We define the L p -realization of A(·, ·, D) by  We verify now Assumptions (A0)-(A3) and (AC). It is well known that A p (ω) generates a compact analytic semigroup in L p (G) for every ω ∈ . Moreover, the mapping ω → A p (ω)v is measurable for every smooth function v ∈ C ∞ (G). This entails the measurability of the mapping ω → A p (ω)v for every v ∈ L p (G). Consequently, the Assumptions (A0)-(A2) and (AC) are satisfied. We only need to show the Hölder continuity of the mapping t → A p (θ t ω) to verify (A3). therein. Several PDEs where the linear part has this structure have been investigated, see, e.g., [27] where a random nonautonomous version of the Fisher-KPP equation is considered. The asymptotic dynamics of the solutions as t tends to infinity has been characterized depending on the behavior of a.