{'fraction_non_alphanumeric': 0.06470529017969497, 'fraction_numerical': 0.015377258770825993, 'mean_word_length': 3.6295001747640687, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 27, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 17, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Singular Stochastic Allen-Cahn equations with dynamic boundary conditions

We prove a well-posedness result for stochastic Allen-Cahn type equations in a bounded domain coupled with generic boundary conditions. The (nonlinear) flux at the boundary aims at describing the interactions with the hard walls and is motivated by some recent literature in physics. The singular character of the drift part allows for a large class of maximal monotone operators, generalizing the usual double-well potentials. One of the main novelties of the paper is the absence of any growth condition on the drift term of the evolution, neither on the domain nor on the boundary. A well-posedness result for variational solutions of the system is presented using a priori estimates as well as monotonicity and compactness techniques. A vanishing viscosity argument for the dynamic on the boundary is also presented.


Introduction
Allen-Cahn type equations were introduced within the Van der Waals theory of phase transitions as a basic model to describe the evolution of a two-phase fluid. Generally, the unknown process is called (non-conserved) order parameter and represents the normalized density of one of the two involved phases. The starting point of the theory is the definition of a free energy functional associated to the order parameter x(·), which is given by (1.1) jumps between the two pure phases, penalizing the variation of x. Then, the Allen-Cahn equation can be viewed as the L 2 -gradient flow of the Ginzburg-Landau free energy (1.1), i.e. the semilinear parabolic PDE of the form In order to model the thermal fluctuation of the system it is quite natural to perturb the equation with a random force and, from a phenomenological point of view, the choice of space time white noise seems vary natural. Unfortunately, nonlinear equations such as (1.2) become ill-posed in space dimensions N ≥ 2 as soon as a white noise term is added. A classical way to bypass the problem is to smooth out the noise via a suitable covariance operator. Given D ⊆ R N (N ≥ 2) a smooth bounded domain with smooth boundary Γ, what we end up with is the so called stochastic Allen-Cahn equation: with a given initial datum where T > 0 is a fixed final time, W is a cylindrical Wiener process on a separable Hilbert space U and B is a Hilbert-Schmidt operator from U to L 2 (D) depending on x as well.
In the classical literature on the deterministic and stochastic Allen-Cahn equation, the order parameter is usually assumed to satisfy homogeneous Neumann conditions on the boundary Γ, which may be interpreted as a null interaction of the phase-transition with the hard walls. In this setting, well-posedness results for the stochastic Allen-Cahn equation can be found e.g. in the monograph [13], within the framework of dissipative SPDEs.
More recently, the class of energy functionals has been extended in several models to take into account also a possible interaction of the phase transition phenomenon with the hard walls. The main idea is to require that the free energy functional could (possibly) penalize the variation of x and the pure states also on the boundary Γ: to this aim, if we denote with ∇ Γ the surface gradient on the boundary, the form of the functional becomes where ε ≥ 0 is fixed and F Γ : R → [0, +∞) is another smooth double-well potential acting on the boundary. Arguing as before, the L 2 -gradient flow of this generalized free energy describes the following system where the symbol ∂ n denotes the outward normal derivative on Γ and ∆ Γ is the usual Laplace-Beltrami operator. Note that the presence of a free energy term on the boundary leads to non-standard dynamic conditions on the boundary (i.e. that involve the time derivative of x on Γ), in contrast with the classical homogeneous Neumann conditions for x. Let us point out that in the last few years there has been a lot of interest in describing phase separation phenomena in confined systems under general boundary conditions, see e.g. [2], [15].
Finally, since the physical system may be subject also to a thermal fluctuation in the boundary (due both to the diffusion from the interior of D and to a boundary noise) the natural idea is to perturb the equation satisfied by x on Γ the same way that we have done for the one in the interior of D. These considerations lead to consider stochastic boundary conditions of the type where here W Γ is a cylindrical Wiener process on another separable Hilbert space U Γ , independent of W , and B Γ is a random time-dependent Hilbert-Schmidt operator from U Γ to L 2 (Γ) of multiplicative type.
In the present paper we are interested in studying the stochastic system arising from the equations (1.3) and (1.5) from a more general mathematical perspective. The main extension that we carry out concerns the form of the double-well potentials: more precisely, instead of assuming that F and F Γ are smooth functions on R, we simply require that F = j + G and F Γ = j Γ + G Γ , where j, j Γ : R → [0, +∞) are given convex functions with subdifferentials β = ∂j and β Γ = ∂j Γ everywhere defined, respectively, and G, G Γ are smooth functions with Lipschitz differentials π, π Γ , respectively. This means essentially that we are looking at the double-well potentials F and F Γ as sufficiently smooth concave perturbations of convex potentials. Bearing in mind these considerations, in the present paper we are concerned with the following system: dx t − ∆x t dt + β(x t ) dt + π(x t ) dt ∋ B(t, x t ) dW t in (0, T ) × D (1. 6) x = y in (0, T ) × Γ (1.7) dy t + ∂ n x t dt − ε∆ Γ y t dt + β Γ (y t ) dt + π Γ (y t ) dt ∋ B Γ (t, y t ) dW Γ t in (0, T ) × Γ (1.8) x(0) = x 0 in D (1.9) (1.10) More specifically, we aim at proving well-posedness for problem (1.6)-(1.10) both in the case ε > 0 and ε = 0, as well as a suitable continuity of the solutions with respect to the parameter ε ≥ 0, under no restrictive growth assumptions on the potentials.
Mathematical results concerning the Allen-Cahn (and similarly Cahn-Hillard) equation have been obtained recently in the framework of generalized (deterministic) boundary dynamics. Let us mention e.g. [9,11,12,16,22] and the references therein. On the contrary, not much is known on the stochastic counterpart, where one is interested in the dynamical impact of a noise term on the boundary. In this direction we have to mention [5,7], where the authors study (nonlinear) diffusion problems with stochastic boundary conditions in a variational framework, and [33] where long-time properties of the Cahn-Hillard equation are investigated.
Concerning general well-posedness results for stochastic PDEs, let us mention [18] and the references therein, where unique existence of analytical strong solutions for a large class of SPDEs of gradient type is exhibited. In that case, a crucial hypothesis used by the author is the sub-homogeneous character of the potential, which unfortunately forbids e.g. exponential growth. In order to avoid any growth condition, in the present paper we only ask that D(β) = R. Still, this hypothesis does not seem to be optimal: it is not needed for the well-posedness of deterministic systems (see e.g. [9]), whereas in the stochastic formulation seems to be essential, at least in the approach we develop. Actually, in our case of interest β = ∂j with j being a convex potential, this restriction on the domain is the most general assumption in literature: it was considered for the first time in [3] in a problem related to existence of semilinear Laplace-driven stochastic equations and then in [23] when studying well-posedness for a class of abstract semilinear SPDEs with singular drift, avoiding any conditions on the growth of β. Again, let us remark that we are not able to consider a graph of the formβ = ∂I [−1,1] , but only an approximation of it defined everywhere in R. With respect to the result obtained in [23], tailored for a large class of singular dissipative SPDEs, here we focus on a precise choice of diffusion operator, the Laplacian. This is motivated by the physical description of the model, but greater generality can be achieved without any substantial change in the proof. Within this framework, the key idea is to get good estimates in expectation and produce the pathwise counterpart in a set Ω ′ of probability 1 using a suitable regularization on the noise.
The strategy of the proof is as follows. We start by rewriting the system with additive noise as an equation for the pair (x, y) in the product space L 2 (D) × L 2 (Γ). Here we develop a variational approach à la Krylov, Rozovskiȋ and Pardoux. In particular we define a suitable Gelfand triple and we smooth out the equation via Yosida approximations of the singular part. In this way, the approximated system satisfies the usual assumptions and a version of the Itô formula can be applied. We derive estimates in expectation of the solution as well as of the monotone maps. By compactness we pass to the limit pathwise to get a candidate limit equation and we identify the drift part as an element of the maximal monotone graphs β and β Γ . Then we recover uniqueness of the solution which is essential to infer measurability in ω of the limit. At last we generalize the well-posedness result also to noises of multiplicative type using a standard fixed-point argument.
A crucial step to get uniqueness of the solution is the application of the Itô formula, for which a suitable smoothing of the equation is required. A classical way of proceeding is to apply the resolvent operator of the diffusion to the equation itself. In our case of interest, rewriting the equation in the product space L 2 (D)×L 2 (Γ), the diffusion operator is not "standard" and has the form where y = τ x is the trace part. Hence, given a couple of functions (f, g), one needs to study the smoothing effect of the resolvent (I + δC ε ) −1 through an ad hoc regularity analysis of the associated elliptic system Precisely, what we show is the ultracontractivity of (I + δC ε ) −1 from To this aim, we generalize a classical regularity result by Stampacchia for elliptic equations with homogeneous boundary conditions contained in [31] and subsequently prove a version of the maximum principle with data (f, g) ∈ L 1 (D) × L 1 (Γ). Let us note that the study of this operator forces us to impose some additional constraints on the relative growth of β and β Γ , which are indeed quite natural from the point of view of the physical applications. All the results mentioned above are collected in the Appendix.
Throughout the paper we use the parameter ε to indicate the presence of the diffusion operator ε∆ Γ on the boundary. The presence/absence of this term creates a gap between the effective domain of C ε , when ε > 0 and ε = 0. Form the physical point of view, a vanishing viscosity argument on the boundary becomes interesting as it is related to the formation of sharp interfaces between the two phases. What we show is the well-posedness of the problem in the singular case ε = 0 as well as the continuous dependence of the solutions to the system (1.6), (1.8) with respect to the variation of ε ≥ 0.
The paper is organised as follows: in the first section we introduce the notations, assumptions and we present the main results. Section 3 is devoted to proving the wellposedness of the system with additive noise: here, we study the approximated equation and we pass to the limit using compactness arguments. In Sections 4 and 5 we extend the previous result to the Allen-Cahn equation with multiplicative noise and we study the asymptotic behaviour of the system as ε → 0. Finally, in the Appendix, we derive the smoothing properties of the diffusion operator C ε .

Notation, setting and main results
In the section we state the notation that we use and the precise assumptions of the work; moreover, the concept of solution and the main results are presented.

Notation
Throughout the paper, (Ω, F , F, P) is a filtered probability space, with the filtration F = (F t ) t∈[0,T ] satisfying the so-called "usual conditions" (i.e. it is saturated and right continuous). As we have anticipated, D ⊆ R N is a smooth bounded domain with smooth boundary Γ and T > 0 is a fixed time.
If U is a Banach space, for any t ∈ [0, T ] we use the classical notations L p (Ω, F t , P; U) and L p (0, T ; U) for the classes of U-valued p-Bochner-integrable functions on Ω and (0, T ), respectively (without specifying the σ-algebra and the probability measure if t = T ). The symbol C 0 w ([0, T ]; U) denotes the space of continuous functions from [0, T ] to the space U endowed with the weak topology. Moreover, if U 1 and U 2 are separable Hilbert Spaces, we may write L (U 1 , U 2 ) and L 2 (U 1 , U 2 ) to indicate the spaces of the linear continuous operators and Hilbert-Schmidt operators from U 1 to U 2 , respectively.
For every a, b ≥ 0, we use the classical notation a b to mean that there exists a positive constant C such that a ≤ Cb.

Assumptions
We precise here the assumptions that are in order throughout the work.
Since β and β Γ are everywhere defined, j and j Γ are actually continuous. Moreover, the convex conjugates of j and j Γ , i.e. We need to make some further hypotheses on j and j Γ . Firstly, we require a symmetry property for the growth of the two potentials at infinity, namely which is very common in literature (see [3,6,[23][24][25]28]). Secondly, a natural assumption to make is that which means essentially that j and j Γ control each other at +∞. If we keep in mind the physical interpretation of the problem, (H1) is very reasonable and can be reinterpreted as the requirement that the nonlinear flux on the boundary is of the same type as the one in the interior of the domain.
However, note that condition (H1) is much stronger than the corresponding one in the deterministic case, in which it is sufficient to assume just one of the two inequalities (see [9]). Consequently, for sake of completeness, it is worth introducing two other possible hypotheses, in which the potentials are allowed to have different growth at infinity, provided that they are bounded by specific polynomial functions: Let us comment on these conditions, focusing in particular on the cases N = 2, 3, which are the most interesting in terms of applications. In (H2) we are requiring that j is controlled by j Γ and that j Γ is bounded by a polynomial of degree six if N = 3, or by any generic polynomial if N = 2. The second hypothesis requires instead that j Γ is controlled by j, and it depends on wether we are working with ε > 0 or ε = 0. In the case ε > 0, we can assume that j has any polynomial growth if N = 3, or any arbitrary growth if N = 2. In the case ε = 0, j has to be controlled by a polynomial of degree four if N = 3, or by any generic polynomial if N = 2. In particular, note that the classical double-well potentials of degree four are included in the interesting cases N = 2, 3.
Polynomial growths of this type for the potentials have been widely used in the deterministic setting, also in the framework of Cahn-Hilliard and quasilinear equations. Among the great literature, we can mention the works [10,14,16,17,30] and the references therein.
We introduce also the multivalued operator and the proper, convex and lower semicontinuous function it is not difficult to check that γ is maximal monotone on R 2 and 0 ∈ γ(0) , D(γ) = R 2 , γ = ∂k , k(0) = 0 .
Finally, we define which is Lipschitz continuous with Lipschitz constant C P := max{C π , C π Γ }.
Assumptions on the variational setting. For any ε ≥ 0, we define the spaces endowed with their natural norms · H , · H Γ , · V and · V ε Γ , respectively. We also use the notations (·, ·) H , (·, ·) H Γ , ·, · V and ·, · V ε Γ for the standard scalar products of H and H Γ and the duality pairings between V * and V , (V ε Γ ) * and V ε Γ , respectively. In this way, for every ε ≥ 0, (V, H, V * ) and (V ε Γ , H Γ , (V ε Γ ) * ) are Hilbert triplets with compact inclusions V c ֒→ H and V ε Γ c ֒→ H Γ . The variational setting is obtained considering the product spaces Let us recall that, using the continuity of τ : H 1 (D) → H 1/2 (Γ), an equivalent norm on V ε is given by Notice that the above equivalence is not uniform in ε. It is clear that H is a Hilbert space with respect to the scalar product and that, for every ε ≥ 0, V ε is included in H continuously and densely, so that (V ε , H, V * ε ) is a Hilbert triplet. Now, it is natural to introduce the operator Then it is immediate to check that for every (ϕ, ψ) ∈ V ε we have Analysing separately the cases ε > 0 and ε = 0 and using the fact that τ : is continuous, it is not difficult to check that the previous expression defines a linear continuous functional on V ε for every ε ≥ 0. Hence, the operator A ε can be extended to In the sequel, we will denote by C ε : so that we have the representation A ε = C ε + P, where we have used the same symbol for P and its corresponding Lipschitz operator induced on H.
Assumptions on the noises. Let W and W Γ be two independent cylindrical Wiener processes on two separable Hilbert spaces U and U Γ , respectively. We introduce Then, setting we have that W is a cylindrical Wiener process on U and B is progressively measurable. Moreover, we assume that B and B Γ are Lipschitz-continuous and at most with linear growth in their third arguments, uniformly on Ω × [0, T ], i.e. that there exists a positive constant C such that Then, it is clear that the same hypotheses hold also for B in its corresponding spaces.
Assumption on the initial datum. We assume that the initial datum satisfies

Formulation of the problem and main results
In this setting, we can write the SPDE of the joint process (x t , y t ) as follows: Definition 2.2. We say that problem (2.1) is well-posed for a given ε ≥ 0 if for any initial datum (x 0 , y 0 ) ∈ L 2 (Ω, F 0 , P; H) there exists a unique strong solution to (2.1) in the sense of Definition 2.1 and the following solution map is Lipschitz-continuous: and (x, y, ξ, ξ Γ ) be the unique strong solutions to (2.1) with additive noise given by Theorem 2.3 in the cases ε > 0 and ε = 0, respectively. Then for every sequence {ε k } k∈N , with ε k ց 0 as k → ∞, we have, as k → ∞, (2.14)

Well-posedness
First of all, we prove well-posedness for the problem (2.1) with additive and more regular noise. Namely, let Z be a separable Hilbert space such that Z ֒→ V ε ∩ L ∞ (D) × L ∞ (Γ) (such a Z exists thanks to the Sobolev embeddings theorems) and consider the problem The hypothesis on B will be removed at the end of the section. Existence of a solution is proved using a suitable approximation on the equation and then passing to the limit using compactness results and monotonicity arguments. Continuous dependence on the data is obtained using an appropriate version of Itô's formula.
Throughout the section, ε ≥ 0 is a fixed constant, so that the argument fits both to the case ε = 0 and ε > 0 at the same time; when two different approaches are needed, we will specify it explicitly.

The approximated problem
For any λ ∈ (0, 1), let β λ , β Γ,λ , j λ and j Γ,λ denote the Yosida approximations of the graphs β and β Γ and the Moreau regularizations of the functions j and j Γ , respectively. With this notation, it is a standard matter to check that the Yosida approximation of γ and the Moreau regularization of k are given by We consider the following approximated problem: For sake of simplicity, let us use the notation so that we can write the approximated problem as We recall some properties of the operator A ε,λ in the following lemma.
Lemma 3.1. The operator A ε,λ : V ε → V * ε is hemicontinuous, weakly monotone, weakly coercive and bounded. More specifically, there exist two positive constants C and C ε,λ , with the first being independent of ε and λ, such that the following conditions hold for any Proof. Firstly, let (r 1 , s 1 ), (r 2 , s 2 ), (r 3 , s 3 ) ∈ V ε : for any t ∈ R we have and by the Lipschitz continuity of π, π Γ , β λ and β Γ,λ , the right-hand side is a continuous function of t. Hence, A ε,λ is hemicontinuous. Secondly, using the Lipschitz continuity of π and π Γ and the monotonicity of β λ and β Γ,λ , we have H , from which the weak monotonicity. Moreover, using the Lipschitz continuity of π and π Γ , a similar computation leads to Finally, for any (ϕ, ψ) ∈ V ε , by the Lipschitz continuity of π, π Γ , β λ and β Γ,λ , using the Hölder inequality and renominating the positive constant C at each passage we have from which the boundedness follows.
The previous lemma ensures that the approximated problem is well-posed according to the classical variational approach by Pardoux, Krylov and Rozovskiȋ (see [21,26,27]) in the Gelfand triple (V ε , H, V * ε ). Hence, for any λ ∈ (0, 1) there is a unique strong solution to the approximated problem (3.3)-(3.4). Moreover we can exhibit some a priori estimates, as it is shown in the following lemmata.
Lemma 3.2. There exists a constant K such that the following inequality holds Proof. The proof relies on the application of the version of Itô formula introduced in [21]. Precisely, we have for every t ∈ [0, T ] and P-almost surely that Using the Lipschitzianity of π and π Γ and the weak coercivity of A ε,λ , we deduce that for a positive constant C independent of λ. Taking the supremum in time and expectation, thanks to the Gronwall lemma we get A direct consequence of the Burkholder-Davis-Gundy and Young inequalities (the reader can refer to [23, Lem. 3.1]) is that for every δ > 0 we have From the arbitrariness of δ, we conclude choosing δ small enough.
Lemma 3.3. There exists Ω ′ ∈ F with P(Ω ′ ) = 1 such that, for every ω ∈ Ω ′ , there exists a positive constant K = K(ω) satisfying Proof. To shorten the notation, we will use the notation B · W to mean the stochastic integral of B with respect to W. The approximated equation can be written as Moreover, thanks to (3.2) and the choice of Z we have Let then Ω ′ ∈ F with P(Ω ′ ) = 1 such that the two previous relations hold and fix ω ∈ Ω ′ . Testing (deterministically) the first equation by (x λ , y λ ) − B · W, we get Rearranging the terms, using the regularity of B · W, the Young inequality and the Lipschitzianity of π and π Γ we infer that , by continuity of j and j Γ we have that j(2(B · W )) ∈ L 1 ((0, T ) × D) and j Γ (2(B Γ · W Γ )) ∈ L 1 ((0, T ) × Γ). Hence, recalling that on the left-hand side rearranging the terms and using the Gronwall lemma we can conclude.
Proposition 3.4. For any ω ∈ Ω ′ we can extract a subsequence λ ′ = λ ′ (ω) of λ for which the following convergences hold as λ ′ → 0: Proof. Let us fix ω ∈ Ω ′ . The first two convergences follow from Lemma 3.3. Regarding the third one, recall that for any u ∈ R we have β λ (u) ∈ ∂j ((I + λβ) −1 u) = β ((I + λβ) −1 u) and β λ (u)(I + λβ) −1 u ≥ 0, so that thanks to the contraction property of the resolvent operator. Thanks to Lemma 3.3 and since j * is superlinear at infinity, for any ω ∈ Ω ′ the sequence (β λ (x λ (ω))) λ turns out to be weakly relatively compact in L 1 ((0, T ) × D) thanks to the de la Vallée Poussin criterion along with Dunford-Pettis theorem. Hence, we can extract a subsequence λ ′ (ω) which satisfies the required convergence. The same reasoning can be applied to β Γ,λ to get the fourth convergence statement. Finally, it remains to show the strong convergence of (x λ ′ , y λ ′ ) in L 2 (0, T ; H). To this end, going back to from Lemma 3.3, the boundedness of A ε and the fact that V * ε ֒→ Z * , we have Hence, by Lemma 3.2, d dt ((x λ , y λ ) − B · W) L 1 (0,T ;Z * ) is uniformly bounded in λ and we can apply Simon's theorem (see [29,Cor. 4,p. 85]) to get that (x λ , y λ ) is relatively compact in L 2 (0, T ; H). Then the weak convergence of ( which is the required convergence.

The limit problem
Now we are ready for the proof of the well-posedness of the equation (3.1), where the noise enter the system in an additive fashion. We divide the proof in several steps: Identification of the limit. Fix ω ∈ Ω ′ : in the sequel we do not emphasize the ω-dependence as no confusion can arise. By Proposition 3.4, (x λ ′ , y λ ′ ) → (x, y) strongly in L 2 (0, T ; H), and (x λ ′ (t), y λ ′ (t)) → (x(t), y(t)) for a.e. t ∈ [0, T ] up to passing to a further subsequence. As for the A ε -part we write A ε = C ε + P. Firstly, we employ weak . This is straightforward setting φ 0 ∈ V ε and choosing as a test function φ := s → 1 [0,t] (s)φ 0 ∈ L 2 (0, T ; V ε ). Secondly, we simply use the strong convergence (x λ ′ , y λ ′ ) → (x, y) in L 2 (0, T ; H) and the Lipschitz continuity of P to get t 0 P(x λ ′ (s), y λ ′ (s)) ds → for all t ∈ [0, T ]. Summing up all the previous convergences, we get the limit equation It remains to show that (ξ, ξ Γ ) belongs to (β(x), β Γ (y)) a.e. in (0, T )×D. Let use the strong convergence of (x λ ′ , y λ ′ ) → (x, y) in L 2 (0, T ; H) to extract a subsequence (still denoted with λ ′ ) so that both (I + λ ′ β) −1 x λ ′ and (I + λ ′ β Γ ) −1 y λ ′ would converge a.e. in (0, T ) × D to x and y, respectively. At this point remember that because of Lemma 3.2 and the fact that ω ∈ Ω ′ . Then we are in position to apply Brezis' lemma [8,Thm. 18,p. 126] and we get the identification ξ ∈ β(x). The same reasoning holds for the monotone operator β Γ on the boundary and the claim is proved.
Finally, using the lower semicontinuity of the convex integrands for the weak convergence and Lemma 3.3, we have that the limit solution satisfies where C 1 , C 2 are constants which depend only on ω.
Uniqueness. Here we show a conditional uniqueness result for equation (3.5) with ω fixed. By contradiction, suppose that (x i , y i , ξ i , ξ Γ,i ) satisfy equation (3.5) and are such that and it is enough to show that (x,ȳ) = 0 and (ξ,ξ Γ ) = 0 for every t ∈ [0, T ]. The idea is to smooth out the equation, recover some structure condition and then use strong convergence to pass to the limit in the approximation. To this end, we multiply the above equation by a power (high enough) of the resolvent of C ε . Remember that the strong and that we are able to show regularizing properties of the resolvent operator associated to C ε , which are presented in detail in the Appendix. In particular, we can prove (see Corollary A.6) that there is m ∈ N so that (I + δC ε ) −m is ultracontractive, i.e. that it is well-defined, linear and continuous from L 1 (D) × L 1 (Γ) to L ∞ (D) × L ∞ (Γ). If we use the notation (u, v) δ = (u δ , v δ ) := (I + δC ε ) −m (u, v) for any (u, v) for which it makes sense and for any δ > 0, then we have for every t ∈ [0, T ]. Applying energy estimates (with ω still fixed) of the form e −rt · 2 H , with r > 0, we get By the monotonicity of C ε , the Lipschitz continuity of P and the Gronwall lemma we have Thanks to Lemma A.7 we know that (x,ȳ) δ (t) → (x,ȳ)(t) in H for all t ∈ [0, T ]. This yields the convergence of (x,ȳ) δ (t) 2 H towards (x,ȳ)(t) 2 H . Concerning the second term, let us carefully rewrite it using the projection of the resolvent map defined in the Appendix.
Regularity in ω. So far we have worked with ω ∈ Ω ′ fixed. Now we are interested in showing the regularity of the solution with respect to ω, starting from the issue of measurability. We follow the same ideas developed in [23], which we briefly resume here for the convenience of the reader, with obvious adaptations to our setting. The key point is the uniqueness result obtained for the solution, which guarantees that the subsequence λ ′ actually does not depend on ω ∈ Ω ′ . The proof relies on a standard argument: from any subsequence of λ we can extract a further subsequence λ ′ , depending on ω, such that all the convergences obtained in Proposition 3.4 take place. Since the limit is unique, then the same results holds for the entire sequence λ which does not depend on ω anymore.
Let us show some details for F λ ; the same technique can be adopted for F Γ,λ . Firstly, for any l ∈ L ∞ (Ω), setting j 0 (·) := j * (·/M) and M = 1/ ( g L ∞ ((0,T )×D) ∨ 1)( l L ∞ (Ω) ∨ 1) , we use Jensen's inequality to get The last term is bounded due to Lemma 3.2: then by the de la Vallée Poussin criterion the sequence F λ l is uniformly integrable. This yields the strong convergence F λ l → F l in L 1 (Ω) thanks to Vitali's theorem. From the arbitrariness of the function l we get weak convergence in L 1 (Ω) of (F λ ) λ , thus the convergence of ξ λ to ξ in L 1 (Ω × (0, T ) × D). Applying Mazur's lemma we extract a convex combination (ξ λ ) λ of (ξ λ ) λ which converge strongly to ξ in L 1 (Ω × (0, T ) × D). To conclude note thatξ λ are predictable (as finite convex combination of predictable processes), hence the limit ξ is a predictable L 1 (D)valued process.
Once we have the measurability of the limit solution (x, y, ξ, ξ Γ ), we are interested in showing some estimates in expectation. To this end, we just use lower semicontinuity of the norms with respect to the weak convergence and Fatou's lemma: . Thanks to Lemma 3.2 and Proposition 3.4 the right hand side of all the above inequalities is bounded uniformly in λ, hence and we have the required regularity. Finally, as we did at the beginning of Subsection 3.2, the weak lower semicontinuity of the convex integrands, Lemma 3.2 and the fact that the subsequence λ does not depend on ω ensure also that . This completes the proof of existence of solutions for the problem with additive and more regular noise.

Continuous dependence on the data
Here we prove a continuous dependence result for equation (3.1) for every ε ≥ 0.

Extension to general additive noise
Here we conclude the proof of existence of a solution to (3.1), with general additive noise: namely, we remove the hypothesis (3.2) and we assume just that To do it, we follow the same ideas as in [23,Prop. 5.1]. Using the notation of the Appendix, for any δ ∈ (0, 1) and m given by Corollary A.6 we have Hence, the problem with additive noise B δ admits a unique solution (x δ , y δ , ξ δ , ξ δ Γ ) for what we have already proved. Since by contraction of the resolvent B δ L 2 (Ω;L 2 (0,T ;H)) ≤ B L 2 (Ω;L 2 (0,T ;H)) for every δ ∈ (0, 1), by Lemma 3.2 and the weak lower semicontinuity of the norms we have that for a positive constant K, for every δ ∈ (0, 1). Moreover, by the continuous dependence result proved in the previous section, for every δ, η ∈ (0, 1) we have , since j * and j * Γ are superlinear, the de la Vallée-Poussin and Dunford-Pettis theorems ensure that {ξ δ } δ and {ξ δ Γ } δ are relatively weakly compact in L 1 (Ω × (0, T ) × D) and L 1 (Ω × (0, T ) × Γ), respectively. Taking this information into account, we deduce that (along a subsequence) in L 2 (Ω; L ∞ (0, T ; H)) ∩ L 2 (Ω; L 2 (0, T ; V ε )) , Using these convergences, the weak lower semicontinuity of the convex integrands and the strong-weak closure of the maximal monotone graphs as in the passage to the limit in λ, it is not difficult to check that (x, y, ξ, ξ Γ ) is a strong solution for the problem with additive noise.

Well-posedness with multiplicative noise
We collect here the conclusion of the proof of Theorem 2.3, showing that the original problem (2.1) is well-posed also with multiplicative noise.
It is clear that any fixed point (x, y) for S ε (together with its corresponding (ξ, ξ Γ )) is a strong solution to (2.1). If we introduce for any p ∈ [1, +∞] the norms on L p (0, T ) given by f L p r (0,T ) := f r L p (0,T ) , for f ∈ L p (0, T ), using (3.8) and the Lipschitz continuity of B, it is immediate to check that where the implicit constant is independent of r. Consequently, there exists r large enough such that S ε is a strict contraction on L 2 (Ω; L 2 r (0, T ; H)). By Banach fixed point theorem, there exists a unique fixed point (x, y) for S ε : this implies that (2.1) has a unique strong solution.

Singular Stochastic Allen-Cahn equations
Consequently, it is sufficient to prove Theorem 2.4 for B δ , with δ ∈ (0, 1) fixed. For this reason, it is not restrictive to assume (3.2) here.
The last thing that we have to show is that the convergences (2.10)-(2.14) hold. To this aim, note that (2.10) and (2.14) coincide with (4.2) and (4.6), respectively. Moreover, (2.11)-(2.13) follow from (4.3)-(4.5) using the fact that ε ′ is independent of ω and Vitali's convergence theorem, via a similar argument to the one performed in Section 3.2 when dealing with the regularity in ω.

A An auxiliary result
Let ε ≥ 0 be fixed. We introduce the operator where It is clear C ε is a well-defined linear operator on H if ε > 0. If ε = 0, this is still true since the conditions −∆u ∈ H and u ∈ H 3/2 (D) imply that ∂ n u ∈ H Γ by [20,Thm. 2.27]. Note that C ε is the strong formulation of the linear component of the operator A ε on H.
Proof. It is immediate to see that C ε is monotone. Hence, we only have to check that R(I + C ε ) = H. Let then (f, g) ∈ H: since C ε is coercive on V ε , by the Lax-Milgram lemma there exists a couple (u, v) ∈ V ε such that Taking ϕ ∈ C ∞ c (D) and ψ = 0 in the previous variational formulation, one can see that u satisfies in D the following PDE (in the sense of distributions on D) We recover immediately that also −∆u ∈ H by difference, and the previous equation holds in H. Moreover, for any ψ ∈ C ∞ (Γ), since the trace operator τ has a right-inverse which is linear and continuous from H k−1/2 (Γ) to H k (D) for every k ≥ 1 (see [20,Thm. 2.24]), thanks to the usual Sobolev embeddings results there is ϕ ∈ C ∞ (D) such that ψ = τ ϕ: choosing (ϕ, ψ) in (A.9), integrating by parts on D and taking (A.10) into account, we see that v satisfies (in the sense of distributions on the boundary Γ) Assume now ε > 0. By definition of V ε , we know that v = τ u ∈ H 1 (Γ): hence, thanks to [20,Thm. 3.2] we also deduce that u ∈ H 3/2 (D). Using the result contained in [20,Thm. 2.27] and the facts that u ∈ H 3/2 (D) and −∆u ∈ H, we deduce that ∂ n u ∈ H Γ as well. By comparison in (A.11) we have that v − ε∆ Γ v ∈ H Γ , so that v ∈ H 2 (Γ) by elliptic regularity on the boundary and (A.11) holds in H Γ . Finally, since ∆u ∈ H and v ∈ H 2 (Γ) ֒→ H 3/2 (Γ), thanks again to [20,Thm. 3.2] we deduce that u ∈ H 2 (D). Hence, (u, v) ∈ D(C ε ) and this completes the proof in the case ε > 0; moreover, note that the two results [20, Thm. 2.27 and 3.2] that we have used ensure also that Assume ε = 0. By definition of V ε , now we have u ∈ H 1 (D) and v ∈ H 1/2 (Γ). Moreover, by comparison in (A.11) we have ∂ n u ∈ H Γ : thanks to [20,Thm. 3.2] and the facts that ∆u ∈ H and ∂ n u ∈ H Γ , we deduce that u ∈ H 3/2 (D), and consequently also v ∈ H 1 (Γ). Hence, (u, v) ∈ D(C ε ) and the proof is complete. As before, owing to the results [20, Thm. 2.27 and 3.2] we also have that Remark A.2. Note that the first part of the proof of the previous lemma ensures that C ε can be extended to a linear operator C ε : V ε → V * ε , given by the left-hand side of (A.9), which is still maximal monotone in the sense of Minty-Browder theory (see [4,Ch. 2]). Lemma A.3 (Regularity). Let δ > 0 and k ∈ N. Then, Proof. It is not restrictive to assume that δ = 1. In Lemma A.1 we have already proved the case k = 0: let us only show the case k = 1, since for a general k it follows by induction.
Proof. Since δ > 0 is fixed throughout the proof, we do not use notation for the dependence on δ of the quantities that we introduce. Given (f, g) ∈ L 1 (D) × L 1 (Γ), let us consider {(f n , g n )} n∈N ⊆ H such that (f n , g n ) → (f, g) in L 1 (D) × L 1 (Γ) as n → ∞, and (u n , v n ) := (I + δC ε ) −1 (f n , g n ). Moreover, let {ρ k } k∈N be a sequence of smooth Lipschitzcontinuous increasing functions on R approximating pointwise the maximal monotone graph For example, one can take ρ k (x) = tanh(kx), x ∈ R. Now, we consider equation (A.9) with respect to n, m ∈ N, take the difference and test by ( Using monotonicity and the fact that |ρ k | ≤ 1, letting k → ∞, by the dominated convergence theorem it is immediate to see that This proves that (I + δC ε ) −1 extends uniquely to a linear contraction on L 1 (D) × L 1 (Γ) and the first part of the lemma is proved. Let us focus on the second part. First of all, we need to prove an auxiliary result, which is a generalization of the classical elliptic regularity theorems by Stampacchia (see [31]). Namely, for every p > N and h 0 , . . . , h N ∈ L p (D), by the Lax-Milgram lemma there is a weak solution (z, w) ∈ V ε such that for every (ϕ, ψ) ∈ V ε . Let us prove that (z, w) ∈ L ∞ (D) × L ∞ (Γ) and that there exists C > 0 such that For every k ∈ N, we introduce the Lipschitz function Using the Young inequality, the fact that |G k (z)| ≤ |z| and the monotonicity of G k , we deduce that which can be rewritten as Let us consider first the case N ≥ 3. If we set 2 * := 2N N −2 and 2 * := 2 * 2 * −1 = 2N N +2 , the Sobolev embedding H 1 (D) ֒→ L 2 * (D) on the left-hand side and the Hölder inequality on the right-hand side yield for a positive constant C, from which, thanks to the Young inequality we have Using the fact that h i ∈ L p (D) for i = 0, . . . , N, that p > 2 * (since p > N) we deduce Renominating the constant C, since it is not restrictive to assume that |A k | < 1, it follows Now, using the fact that p > N, it is a standard matter to see that If N = 2, then we know that H 1 (D) ֒→ L r (D) for all r ∈ [1, +∞): using this fact, we repeat the same argument replacing 2 * and 2 * by an arbitrary q ∈ (1, 2) and its conjugate exponent q ′ = q q−1 , respectively. With such a choice, the same computations yield It is easily seen that α > 1 if and only if q < 2p p+2 : since the fact that p > 2 implies that 2p p+2 ∈ (1, 2), we can choose q ∈ (1, 2p p+2 ), getting α > 1 also in the case N = 2, as desired. By [31,Lem. 4.1], we can conclude that z ∈ L ∞ (D) and z L ∞ (D) ≤ C N i=0 h i L p (D) , suitably renominating the positive constant C. Moreover, since w = τ z, we also have that w ∈ L ∞ (Γ) and w L ∞ (Γ) ≤ z L ∞ (D) . We are now ready to complete the proof of the lemma. Testing (A.13) by (u n , v n ), recalling the definition of (u n , v n ) we have taking into account that h 0 , . . . , h N ∈ L p (D) are arbitrary, we deduce that where q := p p−1 is the conjugate exponent of p. Since (f n , g n ) → (f, g) in L 1 (D) × L 1 (Γ), recalling that the operator C ε is linear and p > N, we have that u n → u in W 1,q (D) for every q ∈ [1, N N −1 ), and consequently v n → v in W 1−1/q,q (Γ). This ensures that v = τ u; moreover, letting n → ∞ we have (u, v) W 1,q (D)×W 1−1/q,q (Γ) ≤ C (f, g) L 1 (D)×L 1 (Γ) , from which the thesis follows.
It follows (for a subsequence, which we still denote by δ) that where by a standard density argument (u, v) = (f, g). Moreover, we also have that which implies that (u δ , v δ ) → (f, g) in H for the original sequence. If (f, g) ∈ L 1 (D) × L 1 (Γ), we introduce {(f n , g n )} n∈N ⊆ H such that (f n , g n ) → (f, g) in L 1 (D) × L 1 (Γ) as n → ∞: let (u n,δ , v n,δ ) := (I + δC ε ) −1 (f n , g n ). Using the fact that (I + δC ε ) −1 is a contraction on L 1 (D) × L 1 (Γ) (see Lemma A.4), we have ≤ 2 (f, g) − (f n , g n ) L 1 (D)×L 1 (Γ) + C (u n,δ , v n,δ ) − (f n , g n ) H for a positive constant C independent of n and δ. Now, for any η > 0, there is n ∈ N such that the first term on the right-hand side of the previous expression is controlled by η: for such an n, thanks to what we have already proved, there is δ such that the second term is less or equal than η. Hence, the right-hand side can be made smaller than 2η and the claim is proved.
The first thesis follows summing the two inequalities, while the second can be easily proved with the other (obvious) choice of L 2 ε,δ and G 2 ε,δ .