@unpublished{14597, abstract = {Phase-field models such as the Allen-Cahn equation may give rise to the formation and evolution of geometric shapes, a phenomenon that may be analyzed rigorously in suitable scaling regimes. In its sharp-interface limit, the vectorial Allen-Cahn equation with a potential with N≥3 distinct minima has been conjectured to describe the evolution of branched interfaces by multiphase mean curvature flow. In the present work, we give a rigorous proof for this statement in two and three ambient dimensions and for a suitable class of potentials: As long as a strong solution to multiphase mean curvature flow exists, solutions to the vectorial Allen-Cahn equation with well-prepared initial data converge towards multiphase mean curvature flow in the limit of vanishing interface width parameter ε↘0. We even establish the rate of convergence O(ε1/2). Our approach is based on the gradient flow structure of the Allen-Cahn equation and its limiting motion: Building on the recent concept of "gradient flow calibrations" for multiphase mean curvature flow, we introduce a notion of relative entropy for the vectorial Allen-Cahn equation with multi-well potential. This enables us to overcome the limitations of other approaches, e.g. avoiding the need for a stability analysis of the Allen-Cahn operator or additional convergence hypotheses for the energy at positive times.}, author = {Fischer, Julian L and Marveggio, Alice}, booktitle = {arXiv}, title = {{Quantitative convergence of the vectorial Allen-Cahn equation towards multiphase mean curvature flow}}, doi = {10.48550/ARXIV.2203.17143}, year = {2022}, } @article{10547, abstract = {We establish global-in-time existence results for thermodynamically consistent reaction-(cross-)diffusion systems coupled to an equation describing heat transfer. Our main interest is to model species-dependent diffusivities, while at the same time ensuring thermodynamic consistency. A key difficulty of the non-isothermal case lies in the intrinsic presence of cross-diffusion type phenomena like the Soret and the Dufour effect: due to the temperature/energy dependence of the thermodynamic equilibria, a nonvanishing temperature gradient may drive a concentration flux even in a situation with constant concentrations; likewise, a nonvanishing concentration gradient may drive a heat flux even in a case of spatially constant temperature. We use time discretisation and regularisation techniques and derive a priori estimates based on a suitable entropy and the associated entropy production. Renormalised solutions are used in cases where non-integrable diffusion fluxes or reaction terms appear.}, author = {Fischer, Julian L and Hopf, Katharina and Kniely, Michael and Mielke, Alexander}, issn = {0036-1410}, journal = {SIAM Journal on Mathematical Analysis}, keywords = {Energy-Reaction-Diffusion Systems, Cross Diffusion, Global-In-Time Existence of Weak/Renormalised Solutions, Entropy Method, Onsager System, Soret/Dufour Effect}, number = {1}, pages = {220--267}, publisher = {Society for Industrial and Applied Mathematics}, title = {{Global existence analysis of energy-reaction-diffusion systems}}, doi = {10.1137/20M1387237}, volume = {54}, year = {2022}, } @article{10548, abstract = {Consider a linear elliptic partial differential equation in divergence form with a random coefficient field. The solution operator displays fluctuations around its expectation. The recently developed pathwise theory of fluctuations in stochastic homogenization reduces the characterization of these fluctuations to those of the so-called standard homogenization commutator. In this contribution, we investigate the scaling limit of this key quantity: starting from a Gaussian-like coefficient field with possibly strong correlations, we establish the convergence of the rescaled commutator to a fractional Gaussian field, depending on the decay of correlations of the coefficient field, and we investigate the (non)degeneracy of the limit. This extends to general dimension $d\ge1$ previous results so far limited to dimension $d=1$, and to the continuum setting with strong correlations recent results in the discrete iid case.}, author = {Duerinckx, Mitia and Fischer, Julian L and Gloria, Antoine}, issn = {1050-5164}, journal = {Annals of applied probability}, number = {2}, pages = {1179--1209}, publisher = {Institute of Mathematical Statistics}, title = {{Scaling limit of the homogenization commutator for Gaussian coefficient fields}}, doi = {10.1214/21-AAP1705}, volume = {32}, year = {2022}, } @article{11701, abstract = {In this paper we develop a new approach to nonlinear stochastic partial differential equations with Gaussian noise. Our aim is to provide an abstract framework which is applicable to a large class of SPDEs and includes many important cases of nonlinear parabolic problems which are of quasi- or semilinear type. This first part is on local existence and well-posedness. A second part in preparation is on blow-up criteria and regularization. Our theory is formulated in an Lp-setting, and because of this we can deal with nonlinearities in a very efficient way. Applications to several concrete problems and their quasilinear variants are given. This includes Burgers' equation, the Allen–Cahn equation, the Cahn–Hilliard equation, reaction–diffusion equations, and the porous media equation. The interplay of the nonlinearities and the critical spaces of initial data leads to new results and insights for these SPDEs. The proofs are based on recent developments in maximal regularity theory for the linearized problem for deterministic and stochastic evolution equations. In particular, our theory can be seen as a stochastic version of the theory of critical spaces due to Prüss–Simonett–Wilke (2018). Sharp weighted time-regularity allow us to deal with rough initial values and obtain instantaneous regularization results. The abstract well-posedness results are obtained by a combination of several sophisticated splitting and truncation arguments.}, author = {Agresti, Antonio and Veraar, Mark}, issn = {1361-6544}, journal = {Nonlinearity}, number = {8}, pages = {4100--4210}, publisher = {IOP Publishing}, title = {{Nonlinear parabolic stochastic evolution equations in critical spaces Part I. Stochastic maximal regularity and local existence}}, doi = {10.1088/1361-6544/abd613}, volume = {35}, year = {2022}, } @article{11842, abstract = {We consider the flow of two viscous and incompressible fluids within a bounded domain modeled by means of a two-phase Navier–Stokes system. The two fluids are assumed to be immiscible, meaning that they are separated by an interface. With respect to the motion of the interface, we consider pure transport by the fluid flow. Along the boundary of the domain, a complete slip boundary condition for the fluid velocities and a constant ninety degree contact angle condition for the interface are assumed. In the present work, we devise for the resulting evolution problem a suitable weak solution concept based on the framework of varifolds and establish as the main result a weak-strong uniqueness principle in 2D. The proof is based on a relative entropy argument and requires a non-trivial further development of ideas from the recent work of Fischer and the first author (Arch. Ration. Mech. Anal. 236, 2020) to incorporate the contact angle condition. To focus on the effects of the necessarily singular geometry of the evolving fluid domains, we work for simplicity in the regime of same viscosities for the two fluids.}, author = {Hensel, Sebastian and Marveggio, Alice}, issn = {1422-6952}, journal = {Journal of Mathematical Fluid Mechanics}, number = {3}, publisher = {Springer Nature}, title = {{Weak-strong uniqueness for the Navier–Stokes equation for two fluids with ninety degree contact angle and same viscosities}}, doi = {10.1007/s00021-022-00722-2}, volume = {24}, year = {2022}, } @article{11858, abstract = {This paper is a continuation of Part I of this project, where we developed a new local well-posedness theory for nonlinear stochastic PDEs with Gaussian noise. In the current Part II we consider blow-up criteria and regularization phenomena. As in Part I we can allow nonlinearities with polynomial growth and rough initial values from critical spaces. In the first main result we obtain several new blow-up criteria for quasi- and semilinear stochastic evolution equations. In particular, for semilinear equations we obtain a Serrin type blow-up criterium, which extends a recent result of Prüss–Simonett–Wilke (J Differ Equ 264(3):2028–2074, 2018) to the stochastic setting. Blow-up criteria can be used to prove global well-posedness for SPDEs. As in Part I, maximal regularity techniques and weights in time play a central role in the proofs. Our second contribution is a new method to bootstrap Sobolev and Hölder regularity in time and space, which does not require smoothness of the initial data. The blow-up criteria are at the basis of these new methods. Moreover, in applications the bootstrap results can be combined with our blow-up criteria, to obtain efficient ways to prove global existence. This gives new results even in classical 𝐿2-settings, which we illustrate for a concrete SPDE. In future works in preparation we apply the results of the current paper to obtain global well-posedness results and regularity for several concrete SPDEs. These include stochastic Navier–Stokes equations, reaction– diffusion equations and the Allen–Cahn equation. Our setting allows to put these SPDEs into a more flexible framework, where less restrictions on the nonlinearities are needed, and we are able to treat rough initial values from critical spaces. Moreover, we will obtain higher-order regularity results.}, author = {Agresti, Antonio and Veraar, Mark}, issn = {1424-3202}, journal = {Journal of Evolution Equations}, keywords = {Mathematics (miscellaneous)}, number = {2}, publisher = {Springer Nature}, title = {{Nonlinear parabolic stochastic evolution equations in critical spaces part II}}, doi = {10.1007/s00028-022-00786-7}, volume = {22}, year = {2022}, } @article{12079, abstract = {We extend the recent rigorous convergence result of Abels and Moser (SIAM J Math Anal 54(1):114–172, 2022. https://doi.org/10.1137/21M1424925) concerning convergence rates for solutions of the Allen–Cahn equation with a nonlinear Robin boundary condition towards evolution by mean curvature flow with constant contact angle. More precisely, in the present work we manage to remove the perturbative assumption on the contact angle being close to 90∘. We establish under usual double-well type assumptions on the potential and for a certain class of boundary energy densities the sub-optimal convergence rate of order ε12 for general contact angles α∈(0,π). For a very specific form of the boundary energy density, we even obtain from our methods a sharp convergence rate of order ε; again for general contact angles α∈(0,π). Our proof deviates from the popular strategy based on rigorous asymptotic expansions and stability estimates for the linearized Allen–Cahn operator. Instead, we follow the recent approach by Fischer et al. (SIAM J Math Anal 52(6):6222–6233, 2020. https://doi.org/10.1137/20M1322182), thus relying on a relative entropy technique. We develop a careful adaptation of their approach in order to encode the constant contact angle condition. In fact, we perform this task at the level of the notion of gradient flow calibrations. This concept was recently introduced in the context of weak-strong uniqueness for multiphase mean curvature flow by Fischer et al. (arXiv:2003.05478v2).}, author = {Hensel, Sebastian and Moser, Maximilian}, issn = {1432-0835}, journal = {Calculus of Variations and Partial Differential Equations}, number = {6}, publisher = {Springer Nature}, title = {{Convergence rates for the Allen–Cahn equation with boundary contact energy: The non-perturbative regime}}, doi = {10.1007/s00526-022-02307-3}, volume = {61}, year = {2022}, } @article{8792, abstract = {This paper is concerned with a non-isothermal Cahn-Hilliard model based on a microforce balance. The model was derived by A. Miranville and G. Schimperna starting from the two fundamental laws of Thermodynamics, following M. Gurtin's two-scale approach. The main working assumptions are made on the behaviour of the heat flux as the absolute temperature tends to zero and to infinity. A suitable Ginzburg-Landau free energy is considered. Global-in-time existence for the initial-boundary value problem associated to the entropy formulation and, in a subcase, also to the weak formulation of the model is proved by deriving suitable a priori estimates and by showing weak sequential stability of families of approximating solutions. At last, some highlights are given regarding a possible approximation scheme compatible with the a-priori estimates available for the system.}, author = {Marveggio, Alice and Schimperna, Giulio}, issn = {10902732}, journal = {Journal of Differential Equations}, number = {2}, pages = {924--970}, publisher = {Elsevier}, title = {{On a non-isothermal Cahn-Hilliard model based on a microforce balance}}, doi = {10.1016/j.jde.2020.10.030}, volume = {274}, year = {2021}, } @article{9240, abstract = {A stochastic PDE, describing mesoscopic fluctuations in systems of weakly interacting inertial particles of finite volume, is proposed and analysed in any finite dimension . It is a regularised and inertial version of the Dean–Kawasaki model. A high-probability well-posedness theory for this model is developed. This theory improves significantly on the spatial scaling restrictions imposed in an earlier work of the same authors, which applied only to significantly larger particles in one dimension. The well-posedness theory now applies in d-dimensions when the particle-width ϵ is proportional to for and N is the number of particles. This scaling is optimal in a certain Sobolev norm. Key tools of the analysis are fractional Sobolev spaces, sharp bounds on Bessel functions, separability of the regularisation in the d-spatial dimensions, and use of the Faà di Bruno's formula.}, author = {Cornalba, Federico and Shardlow, Tony and Zimmer, Johannes}, issn = {1090-2732}, journal = {Journal of Differential Equations}, number = {5}, pages = {253--283}, publisher = {Elsevier}, title = {{Well-posedness for a regularised inertial Dean–Kawasaki model for slender particles in several space dimensions}}, doi = {10.1016/j.jde.2021.02.048}, volume = {284}, year = {2021}, } @article{9307, abstract = {We establish finite time extinction with probability one for weak solutions of the Cauchy–Dirichlet problem for the 1D stochastic porous medium equation with Stratonovich transport noise and compactly supported smooth initial datum. Heuristically, this is expected to hold because Brownian motion has average spread rate O(t12) whereas the support of solutions to the deterministic PME grows only with rate O(t1m+1). The rigorous proof relies on a contraction principle up to time-dependent shift for Wong–Zakai type approximations, the transformation to a deterministic PME with two copies of a Brownian path as the lateral boundary, and techniques from the theory of viscosity solutions.}, author = {Hensel, Sebastian}, issn = {2194-041X}, journal = {Stochastics and Partial Differential Equations: Analysis and Computations}, pages = {892–939}, publisher = {Springer Nature}, title = {{Finite time extinction for the 1D stochastic porous medium equation with transport noise}}, doi = {10.1007/s40072-021-00188-9}, volume = {9}, year = {2021}, } @article{9335, abstract = {Various degenerate diffusion equations exhibit a waiting time phenomenon: depending on the “flatness” of the compactly supported initial datum at the boundary of the support, the support of the solution may not expand for a certain amount of time. We show that this phenomenon is captured by particular Lagrangian discretizations of the porous medium and the thin film equations, and we obtain sufficient criteria for the occurrence of waiting times that are consistent with the known ones for the original PDEs. For the spatially discrete solution, the waiting time phenomenon refers to a deviation of the edge of support from its original position by a quantity comparable to the mesh width, over a mesh-independent time interval. Our proof is based on estimates on the fluid velocity in Lagrangian coordinates. Combining weighted entropy estimates with an iteration technique à la Stampacchia leads to upper bounds on free boundary propagation. Numerical simulations show that the phenomenon is already clearly visible for relatively coarse discretizations.}, author = {Fischer, Julian L and Matthes, Daniel}, issn = {0036-1429}, journal = {SIAM Journal on Numerical Analysis}, number = {1}, pages = {60--87}, publisher = {Society for Industrial and Applied Mathematics}, title = {{The waiting time phenomenon in spatially discretized porous medium and thin film equations}}, doi = {10.1137/19M1300017}, volume = {59}, year = {2021}, } @article{9352, abstract = {This paper provides an a priori error analysis of a localized orthogonal decomposition method for the numerical stochastic homogenization of a model random diffusion problem. If the uniformly elliptic and bounded random coefficient field of the model problem is stationary and satisfies a quantitative decorrelation assumption in the form of the spectral gap inequality, then the expected $L^2$ error of the method can be estimated, up to logarithmic factors, by $H+(\varepsilon/H)^{d/2}$, $\varepsilon$ being the small correlation length of the random coefficient and $H$ the width of the coarse finite element mesh that determines the spatial resolution. The proof bridges recent results of numerical homogenization and quantitative stochastic homogenization.}, author = {Fischer, Julian L and Gallistl, Dietmar and Peterseim, Dietmar}, issn = {0036-1429}, journal = {SIAM Journal on Numerical Analysis}, number = {2}, pages = {660--674}, publisher = {Society for Industrial and Applied Mathematics}, title = {{A priori error analysis of a numerical stochastic homogenization method}}, doi = {10.1137/19M1308992}, volume = {59}, year = {2021}, } @article{10005, abstract = {We study systems of nonlinear partial differential equations of parabolic type, in which the elliptic operator is replaced by the first-order divergence operator acting on a flux function, which is related to the spatial gradient of the unknown through an additional implicit equation. This setting, broad enough in terms of applications, significantly expands the paradigm of nonlinear parabolic problems. Formulating four conditions concerning the form of the implicit equation, we first show that these conditions describe a maximal monotone p-coercive graph. We then establish the global-in-time and large-data existence of a (weak) solution and its uniqueness. To this end, we adopt and significantly generalize Minty’s method of monotone mappings. A unified theory, containing several novel tools, is developed in a way to be tractable from the point of view of numerical approximations.}, author = {Bulíček, Miroslav and Maringová, Erika and Málek, Josef}, issn = {1793-6314}, journal = {Mathematical Models and Methods in Applied Sciences}, keywords = {Nonlinear parabolic systems, implicit constitutive theory, weak solutions, existence, uniqueness}, number = {09}, publisher = {World Scientific}, title = {{On nonlinear problems of parabolic type with implicit constitutive equations involving flux}}, doi = {10.1142/S0218202521500457}, volume = {31}, year = {2021}, } @phdthesis{10007, abstract = {The present thesis is concerned with the derivation of weak-strong uniqueness principles for curvature driven interface evolution problems not satisfying a comparison principle. The specific examples being treated are two-phase Navier-Stokes flow with surface tension, modeling the evolution of two incompressible, viscous and immiscible fluids separated by a sharp interface, and multiphase mean curvature flow, which serves as an idealized model for the motion of grain boundaries in an annealing polycrystalline material. Our main results - obtained in joint works with Julian Fischer, Tim Laux and Theresa M. Simon - state that prior to the formation of geometric singularities due to topology changes, the weak solution concept of Abels (Interfaces Free Bound. 9, 2007) to two-phase Navier-Stokes flow with surface tension and the weak solution concept of Laux and Otto (Calc. Var. Partial Differential Equations 55, 2016) to multiphase mean curvature flow (for networks in R^2 or double bubbles in R^3) represents the unique solution to these interface evolution problems within the class of classical solutions, respectively. To the best of the author's knowledge, for interface evolution problems not admitting a geometric comparison principle the derivation of a weak-strong uniqueness principle represented an open problem, so that the works contained in the present thesis constitute the first positive results in this direction. The key ingredient of our approach consists of the introduction of a novel concept of relative entropies for a class of curvature driven interface evolution problems, for which the associated energy contains an interfacial contribution being proportional to the surface area of the evolving (network of) interface(s). The interfacial part of the relative entropy gives sufficient control on the interface error between a weak and a classical solution, and its time evolution can be computed, at least in principle, for any energy dissipating weak solution concept. A resulting stability estimate for the relative entropy essentially entails the above mentioned weak-strong uniqueness principles. The present thesis contains a detailed introduction to our relative entropy approach, which in particular highlights potential applications to other problems in curvature driven interface evolution not treated in this thesis.}, author = {Hensel, Sebastian}, issn = {2663-337X}, pages = {300}, publisher = {Institute of Science and Technology Austria}, title = {{Curvature driven interface evolution: Uniqueness properties of weak solution concepts}}, doi = {10.15479/at:ista:10007}, year = {2021}, } @unpublished{10011, abstract = {We propose a new weak solution concept for (two-phase) mean curvature flow which enjoys both (unconditional) existence and (weak-strong) uniqueness properties. These solutions are evolving varifolds, just as in Brakke's formulation, but are coupled to the phase volumes by a simple transport equation. First, we show that, in the exact same setup as in Ilmanen's proof [J. Differential Geom. 38, 417-461, (1993)], any limit point of solutions to the Allen-Cahn equation is a varifold solution in our sense. Second, we prove that any calibrated flow in the sense of Fischer et al. [arXiv:2003.05478] - and hence any classical solution to mean curvature flow - is unique in the class of our new varifold solutions. This is in sharp contrast to the case of Brakke flows, which a priori may disappear at any given time and are therefore fatally non-unique. Finally, we propose an extension of the solution concept to the multi-phase case which is at least guaranteed to satisfy a weak-strong uniqueness principle.}, author = {Hensel, Sebastian and Laux, Tim}, booktitle = {arXiv}, keywords = {Mean curvature flow, gradient flows, varifolds, weak solutions, weak-strong uniqueness, calibrated geometry, gradient-flow calibrations}, title = {{A new varifold solution concept for mean curvature flow: Convergence of the Allen-Cahn equation and weak-strong uniqueness}}, doi = {10.48550/arXiv.2109.04233}, year = {2021}, } @unpublished{10013, abstract = {We derive a weak-strong uniqueness principle for BV solutions to multiphase mean curvature flow of triple line clusters in three dimensions. Our proof is based on the explicit construction of a gradient-flow calibration in the sense of the recent work of Fischer et al. [arXiv:2003.05478] for any such cluster. This extends the two-dimensional construction to the three-dimensional case of surfaces meeting along triple junctions.}, author = {Hensel, Sebastian and Laux, Tim}, booktitle = {arXiv}, title = {{Weak-strong uniqueness for the mean curvature flow of double bubbles}}, doi = {10.48550/arXiv.2108.01733}, year = {2021}, } @unpublished{10174, abstract = {Quantitative stochastic homogenization of linear elliptic operators is by now well-understood. In this contribution we move forward to the nonlinear setting of monotone operators with p-growth. This first work is dedicated to a quantitative two-scale expansion result. Fluctuations will be addressed in companion articles. By treating the range of exponents 2≤p<∞ in dimensions d≤3, we are able to consider genuinely nonlinear elliptic equations and systems such as −∇⋅A(x)(1+|∇u|p−2)∇u=f (with A random, non-necessarily symmetric) for the first time. When going from p=2 to p>2, the main difficulty is to analyze the associated linearized operator, whose coefficients are degenerate, unbounded, and depend on the random input A via the solution of a nonlinear equation. One of our main achievements is the control of this intricate nonlinear dependence, leading to annealed Meyers' estimates for the linearized operator, which are key to the quantitative two-scale expansion result.}, author = {Clozeau, Nicolas and Gloria, Antoine}, booktitle = {arXiv}, title = {{Quantitative nonlinear homogenization: control of oscillations}}, year = {2021}, } @article{10549, abstract = {We derive optimal-order homogenization rates for random nonlinear elliptic PDEs with monotone nonlinearity in the uniformly elliptic case. More precisely, for a random monotone operator on \mathbb {R}^d with stationary law (that is spatially homogeneous statistics) and fast decay of correlations on scales larger than the microscale \varepsilon >0, we establish homogenization error estimates of the order \varepsilon in case d\geqq 3, and of the order \varepsilon |\log \varepsilon |^{1/2} in case d=2. Previous results in nonlinear stochastic homogenization have been limited to a small algebraic rate of convergence \varepsilon ^\delta . We also establish error estimates for the approximation of the homogenized operator by the method of representative volumes of the order (L/\varepsilon )^{-d/2} for a representative volume of size L. Our results also hold in the case of systems for which a (small-scale) C^{1,\alpha } regularity theory is available.}, author = {Fischer, Julian L and Neukamm, Stefan}, issn = {1432-0673}, journal = {Archive for Rational Mechanics and Analysis}, keywords = {Mechanical Engineering, Mathematics (miscellaneous), Analysis}, number = {1}, pages = {343--452}, publisher = {Springer Nature}, title = {{Optimal homogenization rates in stochastic homogenization of nonlinear uniformly elliptic equations and systems}}, doi = {10.1007/s00205-021-01686-9}, volume = {242}, year = {2021}, } @article{10575, abstract = {The choice of the boundary conditions in mechanical problems has to reflect the interaction of the considered material with the surface. Still the assumption of the no-slip condition is preferred in order to avoid boundary terms in the analysis and slipping effects are usually overlooked. Besides the “static slip models”, there are phenomena that are not accurately described by them, e.g. at the moment when the slip changes rapidly, the wall shear stress and the slip can exhibit a sudden overshoot and subsequent relaxation. When these effects become significant, the so-called dynamic slip phenomenon occurs. We develop a mathematical analysis of Navier–Stokes-like problems with a dynamic slip boundary condition, which requires a proper generalization of the Gelfand triplet and the corresponding function space setting.}, author = {Abbatiello, Anna and Bulíček, Miroslav and Maringová, Erika}, issn = {1793-6314}, journal = {Mathematical Models and Methods in Applied Sciences}, number = {11}, pages = {2165--2212}, publisher = {World Scientific Publishing}, title = {{On the dynamic slip boundary condition for Navier-Stokes-like problems}}, doi = {10.1142/S0218202521500470}, volume = {31}, year = {2021}, } @article{8697, abstract = {In the computation of the material properties of random alloys, the method of 'special quasirandom structures' attempts to approximate the properties of the alloy on a finite volume with higher accuracy by replicating certain statistics of the random atomic lattice in the finite volume as accurately as possible. In the present work, we provide a rigorous justification for a variant of this method in the framework of the Thomas–Fermi–von Weizsäcker (TFW) model. Our approach is based on a recent analysis of a related variance reduction method in stochastic homogenization of linear elliptic PDEs and the locality properties of the TFW model. Concerning the latter, we extend an exponential locality result by Nazar and Ortner to include point charges, a result that may be of independent interest.}, author = {Fischer, Julian L and Kniely, Michael}, issn = {13616544}, journal = {Nonlinearity}, number = {11}, pages = {5733--5772}, publisher = {IOP Publishing}, title = {{Variance reduction for effective energies of random lattices in the Thomas-Fermi-von Weizsäcker model}}, doi = {10.1088/1361-6544/ab9728}, volume = {33}, year = {2020}, }