-
A Critical Drift-Diffusion Equation: Connections to the Diffusion on $\textbf{SL}(2)$
Authors:
Peter Morfe,
Felix Otto,
Christian Wagner
Abstract:
In this note, we connect two seemingly unrelated objects: On the one hand is a two-dimensional drift-diffusion process $X$ with divergence-free and time-independent drift $b$. The drift is given by a stationary Gaussian ensemble, and we focus on the critical case where a small-scale cut-off is necessary for well-posedness and the large-scale cancellations lead to a borderline super-diffusive behav…
▽ More
In this note, we connect two seemingly unrelated objects: On the one hand is a two-dimensional drift-diffusion process $X$ with divergence-free and time-independent drift $b$. The drift is given by a stationary Gaussian ensemble, and we focus on the critical case where a small-scale cut-off is necessary for well-posedness and the large-scale cancellations lead to a borderline super-diffusive behavior. On the other hand is the natural diffusion $F$ on the Lie group $\textbf{SL}(2)$ of matrices of determinant one. As a consequence of this connection, the strongly non-Gaussian character of $F$ transmits to how $X$ depends on its starting point.
△ Less
Submitted 21 October, 2024;
originally announced October 2024.
-
Repetitive Finite Automata With Translucent Letters
Authors:
František Mráz,
Friedrich Otto
Abstract:
Here we propose an extension of the (deterministic and the nondeterministic) finite automaton with translucent letters (DFAwtl and NFAwtl), which lies between these automata and their non-returning variants (that is, the nr-DFAwtl and the nr-NFAwtl). This new model works like a DFAwtl or an NFAwtl, but on seeing the end-of-tape marker, it may change its internal state and continue with its computa…
▽ More
Here we propose an extension of the (deterministic and the nondeterministic) finite automaton with translucent letters (DFAwtl and NFAwtl), which lies between these automata and their non-returning variants (that is, the nr-DFAwtl and the nr-NFAwtl). This new model works like a DFAwtl or an NFAwtl, but on seeing the end-of-tape marker, it may change its internal state and continue with its computation instead of just ending it, accepting or rejecting. This new type of automaton is called a repetitive deterministic or nondeterministic finite automaton with translucent letters (RDFAwtl or RNFAwtl). In the deterministic case, the new model is strictly more expressive than the DFAwtl, but less expressive than the nr-DFAwtl, while in the nondeterministic case, the new model is equivalent to the NFAwtl.
△ Less
Submitted 10 September, 2024;
originally announced September 2024.
-
A Partition Function Estimator
Authors:
Ying-Chih Chiang,
Frank Otto,
Jonathan W. Essex
Abstract:
We propose a simple estimator that allows to calculate the absolute value of a system's partition function from a finite sampling of its canonical ensemble. The estimator utilizes a volume correction term to compensate the effect that the finite sampling cannot cover the whole configuration space. As a proof of concept, the estimator is applied to calculate the partition function for several model…
▽ More
We propose a simple estimator that allows to calculate the absolute value of a system's partition function from a finite sampling of its canonical ensemble. The estimator utilizes a volume correction term to compensate the effect that the finite sampling cannot cover the whole configuration space. As a proof of concept, the estimator is applied to calculate the partition function for several model systems, and the results are compared with the numerically exact solutions. Excellent agreement is found, demonstrating that a solution for an efficient calculation of partition functions is possible.
△ Less
Submitted 4 September, 2024;
originally announced September 2024.
-
A critical drift-diffusion equation: intermittent behavior
Authors:
Felix Otto,
Christian Wagner
Abstract:
We consider a drift-diffusion process with a time-independent and divergence-free random drift that is of white-noise character. We are interested in the critical case of two space dimensions, where one has to impose a small-scale cut-off for well-posedness, and is interested in the marginally super-diffusive behavior on large scales.
In the presence of an (artificial) large-scale cut-off at sca…
▽ More
We consider a drift-diffusion process with a time-independent and divergence-free random drift that is of white-noise character. We are interested in the critical case of two space dimensions, where one has to impose a small-scale cut-off for well-posedness, and is interested in the marginally super-diffusive behavior on large scales.
In the presence of an (artificial) large-scale cut-off at scale L, as a consequence of standard stochastic homogenization theory, there exist harmonic coordinates with a stationary gradient $F_L$; the merit of these coordinates being that under their lens, the drift-diffusion process turns into a martingale.
It has recently been established that the second moments diverge as $\mathbb{E}|F_L|^2\sim\sqrt{\ln L}$ for $L\uparrow\infty$. We quantitatively show that in this limit, and in the regime of small Péclet number, $|F_L|^2/\mathbb{E}|F_L|^2$ is not equi-integrable, and that $\mathbb{E}|{\rm det}F_L|/\mathbb{E}|F_L|^2 $ is small. Hence the Jacobian matrix of the harmonic coordinates is very peaked and non-conformal.
We establish this asymptotic behavior by characterizing a proxy $\tilde F_L$ introduced in previous work as the solution of an Itô SDE w. r. t. the variable $\ln L$, and which implements the concept of a scale-by-scale homogenization based on a variance decomposition and admits an efficient calculus. For this proxy, we establish $\mathbb{E}|\tilde F_L|^4\gg(\mathbb{E}|\tilde F_L|^2)^2$ and $\mathbb{E}({\rm det}\tilde F_L-1)^2\ll 1$. In view of the former property, we assimilate this phenomenon to intermittency. In fact, $\tilde F_L$ behaves like a tensorial stochastic exponential, and as a field can be assimilated to multiplicative Gaussian chaos.
△ Less
Submitted 6 July, 2024; v1 submitted 21 April, 2024;
originally announced April 2024.
-
The Gaussian free-field as a stream function: continuum version of the scale-by-scale homogenization result
Authors:
Peter Morfe,
Felix Otto,
Christian Wagner
Abstract:
This note is about a drift-diffusion process $X$ with a time-independent, divergence-free drift $b$, where $b$ is a smooth Gaussian field that decorrelates over large scales. In two space dimensions, this just fails to fall into the standard theory of stochastic homogenization, and leads to a borderline super-diffusive behavior. In a previous paper by Chatzigeorgiou, Morfe, Otto, and Wang (2022),…
▽ More
This note is about a drift-diffusion process $X$ with a time-independent, divergence-free drift $b$, where $b$ is a smooth Gaussian field that decorrelates over large scales. In two space dimensions, this just fails to fall into the standard theory of stochastic homogenization, and leads to a borderline super-diffusive behavior. In a previous paper by Chatzigeorgiou, Morfe, Otto, and Wang (2022), precise asymptotics of the annealed second moments of $X$ were derived by characterizing the asymptotics of the effective diffusivity $λ_L$ in terms of an artificially introduced large-scale cut-off $L$. The latter was carried out by a scale-by-scale homogenization, and implemented by monitoring the corrector $φ_L$ for geometrically increasing cut-off scales $L^+=ML$. In fact, proxies $(\tildeφ_L,\tildeσ_L)$ for the corrector and flux corrector were introduced incrementally and the residuum $f_L$ estimated.
In this short supplementary note, we reproduce the arguments of the above paper in the continuum setting of $M\downarrow 1$. This has the advantage that the definition of the proxies $(\tildeφ_L,\tildeσ_L)$ becomes more transparent -- it is given by a simple Itô SDE with $\ln L$ acting as a time variable. It also has the advantage that the residuum $f_L$, which is a martingale, can be efficiently and precisely estimated by Itô calculus. This relies on the characterization of the quadratic variation of the (infinite-dimensional) Gaussian driver.
△ Less
Submitted 31 March, 2024;
originally announced April 2024.
-
Vlearn: Off-Policy Learning with Efficient State-Value Function Estimation
Authors:
Fabian Otto,
Philipp Becker,
Ngo Anh Vien,
Gerhard Neumann
Abstract:
Existing off-policy reinforcement learning algorithms often rely on an explicit state-action-value function representation, which can be problematic in high-dimensional action spaces due to the curse of dimensionality. This reliance results in data inefficiency as maintaining a state-action-value function in such spaces is challenging. We present an efficient approach that utilizes only a state-va…
▽ More
Existing off-policy reinforcement learning algorithms often rely on an explicit state-action-value function representation, which can be problematic in high-dimensional action spaces due to the curse of dimensionality. This reliance results in data inefficiency as maintaining a state-action-value function in such spaces is challenging. We present an efficient approach that utilizes only a state-value function as the critic for off-policy deep reinforcement learning. This approach, which we refer to as Vlearn, effectively circumvents the limitations of existing methods by eliminating the necessity for an explicit state-action-value function. To this end, we introduce a novel importance sampling loss for learning deep value functions from off-policy data. While this is common for linear methods, it has not been combined with deep value function networks. This transfer to deep methods is not straightforward and requires novel design choices such as robust policy updates, twin value function networks to avoid an optimization bias, and importance weight clipping. We also present a novel analysis of the variance of our estimate compared to commonly used importance sampling estimators such as V-trace. Our approach improves sample complexity as well as final performance and ensures consistent and robust performance across various benchmark tasks. Eliminating the state-action-value function in Vlearn facilitates a streamlined learning process, enabling more effective exploration and exploitation in complex environments.
△ Less
Submitted 20 June, 2024; v1 submitted 7 March, 2024;
originally announced March 2024.
-
Open the Black Box: Step-based Policy Updates for Temporally-Correlated Episodic Reinforcement Learning
Authors:
Ge Li,
Hongyi Zhou,
Dominik Roth,
Serge Thilges,
Fabian Otto,
Rudolf Lioutikov,
Gerhard Neumann
Abstract:
Current advancements in reinforcement learning (RL) have predominantly focused on learning step-based policies that generate actions for each perceived state. While these methods efficiently leverage step information from environmental interaction, they often ignore the temporal correlation between actions, resulting in inefficient exploration and unsmooth trajectories that are challenging to impl…
▽ More
Current advancements in reinforcement learning (RL) have predominantly focused on learning step-based policies that generate actions for each perceived state. While these methods efficiently leverage step information from environmental interaction, they often ignore the temporal correlation between actions, resulting in inefficient exploration and unsmooth trajectories that are challenging to implement on real hardware. Episodic RL (ERL) seeks to overcome these challenges by exploring in parameters space that capture the correlation of actions. However, these approaches typically compromise data efficiency, as they treat trajectories as opaque \emph{black boxes}. In this work, we introduce a novel ERL algorithm, Temporally-Correlated Episodic RL (TCE), which effectively utilizes step information in episodic policy updates, opening the 'black box' in existing ERL methods while retaining the smooth and consistent exploration in parameter space. TCE synergistically combines the advantages of step-based and episodic RL, achieving comparable performance to recent ERL methods while maintaining data efficiency akin to state-of-the-art (SoTA) step-based RL.
△ Less
Submitted 21 January, 2024;
originally announced January 2024.
-
Rotational covariance restricts available quantum states
Authors:
Fynn Otto,
Konrad Szymański
Abstract:
Quantum states of angular momentum and spin generally are not invariant under rotations of the reference frame. Therefore, they can be used as a resource of relative orientation, which is encoded in the asymmetry of the state under consideration. In this paper we introduce the analytical characterization of the rotational information by parameterizing the group characteristic function by polynomia…
▽ More
Quantum states of angular momentum and spin generally are not invariant under rotations of the reference frame. Therefore, they can be used as a resource of relative orientation, which is encoded in the asymmetry of the state under consideration. In this paper we introduce the analytical characterization of the rotational information by parameterizing the group characteristic function by polynomial functions. By doing so, we show that the set of states achievable through transformations lacking a reference frame (rotationally covariant ones) admits an analytical characterization and can be studied through the use of semidefinite optimization techniques. We demonstrate the developed methods via examples, and provide a physical scenario in which a reference-independent operation performs a metrologically useful operation: the preparation of a state of light improving interferometer sensitivity, which equivalently can be realized as a postprocessing step.
△ Less
Submitted 14 August, 2024; v1 submitted 11 January, 2024;
originally announced January 2024.
-
Lecture notes on Malliavin calculus in regularity structures
Authors:
Lucas Broux,
Felix Otto,
Markus Tempelmayr
Abstract:
Malliavin calculus provides a characterization of the centered model in regularity structures that is stable under removing the small-scale cut-off. In conjunction with a spectral gap inequality, it yields the stochastic estimates of the model.
This becomes transparent on the level of a notion of model that parameterizes the solution manifold, and thus is indexed by multi-indices rather than tre…
▽ More
Malliavin calculus provides a characterization of the centered model in regularity structures that is stable under removing the small-scale cut-off. In conjunction with a spectral gap inequality, it yields the stochastic estimates of the model.
This becomes transparent on the level of a notion of model that parameterizes the solution manifold, and thus is indexed by multi-indices rather than trees, and which allows for a more geometric than combinatorial perspective. In these lecture notes, this is carried out for a PDE with heat operator, a cubic nonlinearity, and driven by additive noise, reminiscent of the stochastic quantization of the Euclidean $φ^4$ model.
More precisely, we informally motivate our notion of the model $(Π,Γ)$ as charts and transition maps, respectively, of the nonlinear solution manifold. These geometric objects are algebrized in terms of formal power series, and their algebra automorphisms. We will assimilate the directional Malliavin derivative to a tangent vector of the solution manifold. This means that it can be treated as a modelled distribution, thereby connecting stochastic model estimates to pathwise solution theory, with its analytic tools of reconstruction and integration. We unroll an inductive calculus that in an automated way applies to the full subcritical regime.
△ Less
Submitted 4 October, 2024; v1 submitted 11 January, 2024;
originally announced January 2024.
-
Lifshitz Transition and Band Structure Evolution in Alkali Metal Intercalated 1Tprime-MoTe2
Authors:
Joohyung Park,
Ayan N. Batyrkhanov,
Jonas Brandhoff,
Marco Gruenewald,
Felix Otto,
Maximilian Schaal,
Saban Hus,
Torsten Fritz,
Florian Göltl,
An-Ping Li,
Oliver L. A. Monti
Abstract:
MoTe2 is a paradigmatic van der Waals layered semimetal with two energetically close electronic phases, the topologically trivial 1Tprime and the low-temperature Td type-II Weyl semimetal phase. The ability to manipulate this phase transition, perhaps towards occurring near room temperature, would open new avenues for harnessing the full potential of Weyl semimetals for high-efficiency electronic…
▽ More
MoTe2 is a paradigmatic van der Waals layered semimetal with two energetically close electronic phases, the topologically trivial 1Tprime and the low-temperature Td type-II Weyl semimetal phase. The ability to manipulate this phase transition, perhaps towards occurring near room temperature, would open new avenues for harnessing the full potential of Weyl semimetals for high-efficiency electronic and spintronic applications. Here, we show that potassium dosing on 1Tprime-MoTe2 induces a Lifshitz transition by a combination of angle-resolved photoemission spectroscopy, scanning tunneling microscopy, x-ray spectroscopy and density functional theory. While the electronic structure shifts rigidly for small concentrations of K, MoTe2 undergoes significant band structure renormalization for larger concentrations. Our results demonstrate that the origin of this electronic structure change stems from alkali metal intercalation. We show that these profound changes are caused by effectively decoupling the 2D sheets, bringing K-intercalated 1Tprime-MoTe2 to the quasi-2D limit, but do not cause a topological phase transition.
△ Less
Submitted 31 January, 2024; v1 submitted 23 December, 2023;
originally announced December 2023.
-
Almost sharp rates of convergence for the average cost and displacement in the optimal matching problem
Authors:
Michael Goldman,
Martin Huesmann,
Felix Otto
Abstract:
In this note we prove estimates for the average cost in the quadratic optimal transport problem on the two-dimensional flat torus which are optimal up to a double logarithm. We also prove sharp estimates on the displacement. This is based on the combination of a post-processing of our quantitative linearization result together with a quasi-orthogonality property.
In this note we prove estimates for the average cost in the quadratic optimal transport problem on the two-dimensional flat torus which are optimal up to a double logarithm. We also prove sharp estimates on the displacement. This is based on the combination of a post-processing of our quantitative linearization result together with a quasi-orthogonality property.
△ Less
Submitted 13 December, 2023;
originally announced December 2023.
-
There is no stationary $p$-cyclically monotone Poisson matching in 2D
Authors:
Martin Huesmann,
Francesco Mattesini,
Felix Otto
Abstract:
We show that for $p>1$ there is no $p$-cyclically monotone stationary matching of two independent Poisson processes in dimension $d=2$. The proof combines the $p$-harmonic approximation result from \cite[Theorem 1.1]{koch23} with local asymptotics for the two-dimensional matching problem. Moreover, we prove a.s. local upper bounds of the correct order in the case $p>1$, which, to the best of our k…
▽ More
We show that for $p>1$ there is no $p$-cyclically monotone stationary matching of two independent Poisson processes in dimension $d=2$. The proof combines the $p$-harmonic approximation result from \cite[Theorem 1.1]{koch23} with local asymptotics for the two-dimensional matching problem. Moreover, we prove a.s. local upper bounds of the correct order in the case $p>1$, which, to the best of our knowledge, are not readily available in the current literature.
△ Less
Submitted 29 November, 2023;
originally announced November 2023.
-
Convergence to the planar interface for a nonlocal free-boundary evolution
Authors:
Felix Otto,
Richard Schubert,
Maria G. Westdickenberg
Abstract:
We capture optimal decay for the Mullins-Sekerka evolution, a nonlocal, parabolic free boundary problem from materials science. Our main result establishes convergence of BV solutions to the planar profile in the physically relevant case of ambient space dimension three. Far from assuming small or well-prepared initial data, we allow for initial interfaces that do not have graph structure and are…
▽ More
We capture optimal decay for the Mullins-Sekerka evolution, a nonlocal, parabolic free boundary problem from materials science. Our main result establishes convergence of BV solutions to the planar profile in the physically relevant case of ambient space dimension three. Far from assuming small or well-prepared initial data, we allow for initial interfaces that do not have graph structure and are not connected, hence explicitly including the regime of Ostwald ripening. In terms only of initially finite (not small) excess mass and excess surface energy, we establish that the surface becomes a Lipschitz graph within a fixed timescale (quantitatively estimated) and remains trapped within this setting. To obtain the graph structure, we leverage regularity results from geometric measure theory. At the same time, we extend a duality method previously employed for one-dimensional PDE problems to higher dimensional, nonlocal geometric evolutions. Optimal algebraic decay rates of excess energy, dissipation, and graph height are obtained.
△ Less
Submitted 2 June, 2024; v1 submitted 25 September, 2023;
originally announced September 2023.
-
MP3: Movement Primitive-Based (Re-)Planning Policy
Authors:
Fabian Otto,
Hongyi Zhou,
Onur Celik,
Ge Li,
Rudolf Lioutikov,
Gerhard Neumann
Abstract:
We introduce a novel deep reinforcement learning (RL) approach called Movement Primitive-based Planning Policy (MP3). By integrating movement primitives (MPs) into the deep RL framework, MP3 enables the generation of smooth trajectories throughout the whole learning process while effectively learning from sparse and non-Markovian rewards. Additionally, MP3 maintains the capability to adapt to chan…
▽ More
We introduce a novel deep reinforcement learning (RL) approach called Movement Primitive-based Planning Policy (MP3). By integrating movement primitives (MPs) into the deep RL framework, MP3 enables the generation of smooth trajectories throughout the whole learning process while effectively learning from sparse and non-Markovian rewards. Additionally, MP3 maintains the capability to adapt to changes in the environment during execution. Although many early successes in robot RL have been achieved by combining RL with MPs, these approaches are often limited to learning single stroke-based motions, lacking the ability to adapt to task variations or adjust motions during execution. Building upon our previous work, which introduced an episode-based RL method for the non-linear adaptation of MP parameters to different task variations, this paper extends the approach to incorporating replanning strategies. This allows adaptation of the MP parameters throughout motion execution, addressing the lack of online motion adaptation in stochastic domains requiring feedback. We compared our approach against state-of-the-art deep RL and RL with MPs methods. The results demonstrated improved performance in sophisticated, sparse reward settings and in domains requiring replanning.
△ Less
Submitted 2 July, 2023; v1 submitted 22 June, 2023;
originally announced June 2023.
-
Lecture notes on the harmonic approximation to quadratic optimal transport
Authors:
Lukas Koch,
Felix Otto
Abstract:
These lecture notes present the quantitative harmonic approximation result for quadratic optimal transport and general measures obtained by Goldman and Otto. The aim is to give a clear presentation of the proof of the main theorem with more motivations, less PDE machinery, and a number of simplifications.
These lecture notes present the quantitative harmonic approximation result for quadratic optimal transport and general measures obtained by Goldman and Otto. The aim is to give a clear presentation of the proof of the main theorem with more motivations, less PDE machinery, and a number of simplifications.
△ Less
Submitted 25 March, 2023;
originally announced March 2023.
-
Combining Reconstruction and Contrastive Methods for Multimodal Representations in RL
Authors:
Philipp Becker,
Sebastian Mossburger,
Fabian Otto,
Gerhard Neumann
Abstract:
Learning self-supervised representations using reconstruction or contrastive losses improves performance and sample complexity of image-based and multimodal reinforcement learning (RL). Here, different self-supervised loss functions have distinct advantages and limitations depending on the information density of the underlying sensor modality. Reconstruction provides strong learning signals but is…
▽ More
Learning self-supervised representations using reconstruction or contrastive losses improves performance and sample complexity of image-based and multimodal reinforcement learning (RL). Here, different self-supervised loss functions have distinct advantages and limitations depending on the information density of the underlying sensor modality. Reconstruction provides strong learning signals but is susceptible to distractions and spurious information. While contrastive approaches can ignore those, they may fail to capture all relevant details and can lead to representation collapse. For multimodal RL, this suggests that different modalities should be treated differently based on the amount of distractions in the signal. We propose Contrastive Reconstructive Aggregated representation Learning (CoRAL), a unified framework enabling us to choose the most appropriate self-supervised loss for each sensor modality and allowing the representation to better focus on relevant aspects. We evaluate CoRAL's benefits on a wide range of tasks with images containing distractions or occlusions, a new locomotion suite, and a challenging manipulation suite with visually realistic distractions. Our results show that learning a multimodal representation by combining contrastive and reconstruction-based losses can significantly improve performance and solve tasks that are out of reach for more naive representation learning approaches and other recent baselines.
△ Less
Submitted 26 June, 2024; v1 submitted 10 February, 2023;
originally announced February 2023.
-
Lecture notes on tree-free regularity structures
Authors:
Felix Otto,
Kihoon Seong,
Markus Tempelmayr
Abstract:
These lecture notes are intended as reader's digest of recent work on a diagram-free approach to the renormalized centered model in Hairer's regularity structures. More precisely, it is about the stochastic estimates of the centered model, based on Malliavin calculus and a spectral gap assumption. We focus on a specific parabolic partial differential equation in quasi-linear form driven by (white)…
▽ More
These lecture notes are intended as reader's digest of recent work on a diagram-free approach to the renormalized centered model in Hairer's regularity structures. More precisely, it is about the stochastic estimates of the centered model, based on Malliavin calculus and a spectral gap assumption. We focus on a specific parabolic partial differential equation in quasi-linear form driven by (white) noise.
We follow a natural renormalization strategy based on preserving symmetries, and carefully introduce Hairer's notion of a centered model, which provides the coefficients in a formal series expansion of a general solution. We explain how the Malliavin derivative in conjunction with Hairer's re-expansion map allows to reformulate this definition in a way that is stable under removing the small-scale regularization.
A few exemplary proofs are provided, both of analytic and of algebraic character. The working horse of the analytic arguments is an ``annealed'' Schauder estimate and related Liouville principle, which is provided. The algebra of formal power series, in variables that play the role of coordinates of the solution manifold, and its algebra morphisms are the key algebraic objects.
△ Less
Submitted 21 September, 2023; v1 submitted 2 January, 2023;
originally announced January 2023.
-
The Gaussian free-field as a stream function: asymptotics of effective diffusivity in infra-red cut-off
Authors:
Georgiana Chatzigeorgiou,
Peter Morfe,
Felix Otto,
Lihan Wang
Abstract:
We analyze the large-time asymptotics of a passive tracer with drift equal to the curl of the Gaussian free field in two dimensions with ultra-violet cut-off at scale unity. We prove that the mean-squared displacement scales like $t \sqrt{\ln t}$, as predicted in the physics literature and recently almost proved by the work of Cannizzaro, Haunschmidt-Sibitz, and Toninelli (2022), which uses mathem…
▽ More
We analyze the large-time asymptotics of a passive tracer with drift equal to the curl of the Gaussian free field in two dimensions with ultra-violet cut-off at scale unity. We prove that the mean-squared displacement scales like $t \sqrt{\ln t}$, as predicted in the physics literature and recently almost proved by the work of Cannizzaro, Haunschmidt-Sibitz, and Toninelli (2022), which uses mathematical-physics type analysis in Fock space. Our approach involves studying the effective diffusivity $λ_{L}$ of the process with an infra-red cut-off at scale $L$, and is based on techniques from stochastic homogenization.
△ Less
Submitted 24 May, 2023; v1 submitted 29 December, 2022;
originally announced December 2022.
-
Atomically Thin Metal-Dielectric Heterostructures by Atomic Layer Deposition
Authors:
Pallabi Paul,
Paul Schmitt,
Vilborg Vala Sigurjonsdottir,
Kevin Hanemann,
Nadja Felde,
Sven Schroeder,
Felix Otto,
Marco Gruenewald,
Torsten Fritz,
Vladimir Roddatis,
Andreas Tuennermann,
Adriana Szeghalmi
Abstract:
Heterostructures increasingly attracted attention over the past several years to enable various optoelectronic and photonic applications. In this work, atomically thin interfaces of Ir/Al2O3 heterostructures compatible with micro-optoelectronic technologies are reported. Their structural and optical properties were determined by spectroscopic and microscopic techniques (XRR, XPS, HRTEM, spectrosco…
▽ More
Heterostructures increasingly attracted attention over the past several years to enable various optoelectronic and photonic applications. In this work, atomically thin interfaces of Ir/Al2O3 heterostructures compatible with micro-optoelectronic technologies are reported. Their structural and optical properties were determined by spectroscopic and microscopic techniques (XRR, XPS, HRTEM, spectroscopic ellipsometry, and UV/VIS/NIR spectrophotometry). The XRR and HRTEM analyses reveal a layer-by-layer growth mechanism of Ir in atomic scale heterostructures, which is different from the typical island-type growth of metals on dielectrics. Alongside, XPS investigations imply the formation of Ir-O-Al bonding at the interfaces for lower Ir concentrations, in contrast to the nanoparticle core-shell structure formation. Precisely tuning the ratio of the constituents ensures the control of the dispersion profile along with a transition from effective dielectric to metallic heterostructures. The Ir coating thickness was varied ranging from a few Å to films of about 7 nm in the heterostructures. The transition has been observed in the structures containing individual Ir coating thicknesses of about 2-4 nm. Following this, show epsilon-near-zero metamaterials with tunable dielectric constants by precisely varying the composition of such heterostructures. Overall, a comprehensive study on structural and optical properties of the metal-dielectric interface of Ir/Al2O3 heterostructures was addressed indicating an extension of the material portfolio available for optical system design.
△ Less
Submitted 16 December, 2022;
originally announced December 2022.
-
Heterostructure films of SiO$_2$ and HfO$_2$ for high power laser optics prepared by plasma-enhanced atomic layer deposition
Authors:
Shawon Alam,
Pallabi Paul,
Vivek Beladiya,
Olaf Stenzel,
Marcus Trost,
Steffen Wilbrandt,
Sven Schroeder,
Gabor Matthaeus,
Stefan Nolte,
Sebastian Riese,
Felix Otto,
Torsten Fritz,
Alexander Gottwald,
Adriana Szeghalmi
Abstract:
Absorption losses and laser-induced damage threshold (LIDT) are considered as the major constraint for the development of optical coatings for high-power laser optics. Such coatings require paramount properties like low losses due to optical absorption, high mechanical stability, and enhanced damage resistance to withstand high-intensity laser pulses. In this work, heterostructure films were devel…
▽ More
Absorption losses and laser-induced damage threshold (LIDT) are considered as the major constraint for the development of optical coatings for high-power laser optics. Such coatings require paramount properties like low losses due to optical absorption, high mechanical stability, and enhanced damage resistance to withstand high-intensity laser pulses. In this work, heterostructure films were developed by the intermixing of SiO$_2$ and HfO$_2$ using plasma-enhanced atomic layer deposition (PEALD) technique. Thin film characterization techniques such as spectroscopic ellipsometry, spectrophotometry, substrate curvature measurements, x-ray reflectivity, and Fourier transform infrared spectroscopy were employed for extracting optical constants, spectral inter-pretation, residual stress, layer formation, and functional groups present in the heterostructures, respectively. These heterostructures demonstrate tunable refractive index, bandgap, and improved optical losses and LIDT properties. The films were incorporated into antireflection coatings (multilayer stacks and graded index coatings) and the LIDT was determined at 355 nm wavelength by the R-on-1 method. Optical absorptions at the reported wavelengths were characterized using photothermal common-path interferometry and laser-induced deflection techniques.
△ Less
Submitted 16 December, 2022;
originally announced December 2022.
-
Deep Black-Box Reinforcement Learning with Movement Primitives
Authors:
Fabian Otto,
Onur Celik,
Hongyi Zhou,
Hanna Ziesche,
Ngo Anh Vien,
Gerhard Neumann
Abstract:
\Episode-based reinforcement learning (ERL) algorithms treat reinforcement learning (RL) as a black-box optimization problem where we learn to select a parameter vector of a controller, often represented as a movement primitive, for a given task descriptor called a context. ERL offers several distinct benefits in comparison to step-based RL. It generates smooth control trajectories, can handle non…
▽ More
\Episode-based reinforcement learning (ERL) algorithms treat reinforcement learning (RL) as a black-box optimization problem where we learn to select a parameter vector of a controller, often represented as a movement primitive, for a given task descriptor called a context. ERL offers several distinct benefits in comparison to step-based RL. It generates smooth control trajectories, can handle non-Markovian reward definitions, and the resulting exploration in parameter space is well suited for solving sparse reward settings. Yet, the high dimensionality of the movement primitive parameters has so far hampered the effective use of deep RL methods. In this paper, we present a new algorithm for deep ERL. It is based on differentiable trust region layers, a successful on-policy deep RL algorithm. These layers allow us to specify trust regions for the policy update that are solved exactly for each state using convex optimization, which enables policies learning with the high precision required for the ERL. We compare our ERL algorithm to state-of-the-art step-based algorithms in many complex simulated robotic control tasks. In doing so, we investigate different reward formulations - dense, sparse, and non-Markovian. While step-based algorithms perform well only on dense rewards, ERL performs favorably on sparse and non-Markovian rewards. Moreover, our results show that the sparse and the non-Markovian rewards are also often better suited to define the desired behavior, allowing us to obtain considerably higher quality policies compared to step-based RL.
△ Less
Submitted 18 October, 2022;
originally announced October 2022.
-
ProDMPs: A Unified Perspective on Dynamic and Probabilistic Movement Primitives
Authors:
Ge Li,
Zeqi Jin,
Michael Volpp,
Fabian Otto,
Rudolf Lioutikov,
Gerhard Neumann
Abstract:
Movement Primitives (MPs) are a well-known concept to represent and generate modular trajectories. MPs can be broadly categorized into two types: (a) dynamics-based approaches that generate smooth trajectories from any initial state, e. g., Dynamic Movement Primitives (DMPs), and (b) probabilistic approaches that capture higher-order statistics of the motion, e. g., Probabilistic Movement Primitiv…
▽ More
Movement Primitives (MPs) are a well-known concept to represent and generate modular trajectories. MPs can be broadly categorized into two types: (a) dynamics-based approaches that generate smooth trajectories from any initial state, e. g., Dynamic Movement Primitives (DMPs), and (b) probabilistic approaches that capture higher-order statistics of the motion, e. g., Probabilistic Movement Primitives (ProMPs). To date, however, there is no method that unifies both, i. e. that can generate smooth trajectories from an arbitrary initial state while capturing higher-order statistics. In this paper, we introduce a unified perspective of both approaches by solving the ODE underlying the DMPs. We convert expensive online numerical integration of DMPs into basis functions that can be computed offline. These basis functions can be used to represent trajectories or trajectory distributions similar to ProMPs while maintaining all the properties of dynamical systems. Since we inherit the properties of both methodologies, we call our proposed model Probabilistic Dynamic Movement Primitives (ProDMPs). Additionally, we embed ProDMPs in deep neural network architecture and propose a new cost function for efficient end-to-end learning of higher-order trajectory statistics. To this end, we leverage Bayesian Aggregation for non-linear iterative conditioning on sensory inputs. Our proposed model achieves smooth trajectory generation, goal-attractor convergence, correlation analysis, non-linear conditioning, and online re-planing in one framework.
△ Less
Submitted 4 October, 2022;
originally announced October 2022.
-
Non-Returning Finite Automata With Translucent Letters
Authors:
František Mráz,
Friedrich Otto
Abstract:
Here we propose a variant of the nondeterministic finite automaton with translucent letters (NFAwtl) which, after reading and deleting a letter, does not return to the left end of its tape, but rather continues from the position of the letter just deleted. When the end-of-tape marker is reached, our automaton can decide whether to accept, to reject, or to continue, which means that it again reads…
▽ More
Here we propose a variant of the nondeterministic finite automaton with translucent letters (NFAwtl) which, after reading and deleting a letter, does not return to the left end of its tape, but rather continues from the position of the letter just deleted. When the end-of-tape marker is reached, our automaton can decide whether to accept, to reject, or to continue, which means that it again reads the remaining tape contents from the beginning. This type of automaton, called a non-returning finite automaton with translucent letters or an nrNFAwtl, is strictly more expressive than the NFAwtl. We study the expressive capacity of this type of automaton and that of its deterministic variant. Also we are interested in closure properties of the resulting classes of languages and in decision problems.
△ Less
Submitted 31 August, 2022;
originally announced August 2022.
-
A tree-free approach to regularity structures: The regular case for quasi-linear equations
Authors:
Pablo Linares,
Felix Otto
Abstract:
We give a motivation and gentle introduction into the regularity structure and model introduced by Otto, Sauer, Smith and Weber, which fall into the framework of Hairer, but have a greedier index set than the one given by trees. We do this here for a simple quasi-linear parabolic equation and assume that the driving noise is so regular that no renormalization is needed. We introduce the abstract m…
▽ More
We give a motivation and gentle introduction into the regularity structure and model introduced by Otto, Sauer, Smith and Weber, which fall into the framework of Hairer, but have a greedier index set than the one given by trees. We do this here for a simple quasi-linear parabolic equation and assume that the driving noise is so regular that no renormalization is needed. We introduce the abstract model space $\mathsf{T}$ and its grading, the pre-model $\mathbfΠ$, the centered model $Π_x$, the structure group $\mathsf{G}$, and the re-centering transformations $Γ_{xy}$. Using integration and reconstruction, we establish the desired estimates on $Π_x$ and $Γ_{xy}$, which here are deterministic since we deal with the regular case.
△ Less
Submitted 21 July, 2022;
originally announced July 2022.
-
Bias in the representative volume element method: periodize the ensemble instead of its realizations
Authors:
Nicolas Clozeau,
Marc Josien,
Felix Otto,
Qiang Xu
Abstract:
We study the Representative Volume Element (RVE) method, which is a method to approximately infer the effective behavior $a_{\text{hom}}$ of a stationary random medium. The latter is described by a coefficient field $a(x)$ generated from a given ensemble $\langle\cdot\rangle$ and the corresponding linear elliptic operator $-\nabla\cdot a\nabla$. In line with the theory of homogenization, the metho…
▽ More
We study the Representative Volume Element (RVE) method, which is a method to approximately infer the effective behavior $a_{\text{hom}}$ of a stationary random medium. The latter is described by a coefficient field $a(x)$ generated from a given ensemble $\langle\cdot\rangle$ and the corresponding linear elliptic operator $-\nabla\cdot a\nabla$. In line with the theory of homogenization, the method proceeds by computing $d = 3$ correctors (d denoting the space dimension).To be numerically tractable, this computation has to be done on a finite domain: the so-called "representative" volume element, i. e. a large box with, say, periodic boundary conditions. The main message of this article is: Periodize the ensemble instead of its realizations. By this we mean that it is better to sample from a suitably periodized ensemble than to periodically extend the restriction of a realization $a(x)$ from the whole-space ensemble $\langle\cdot\rangle$. We make this point by investigating the bias (or systematic error), i. e. the difference between $a_{\text{hom}}$ and the expected value of the RVE method, in terms of its scaling w. r. t. the lateral size $L$ of the box. In case of periodizing $a(x)$, we heuristically argue that this error is generically $O(L^{-1})$. In case of a suitable periodization of $\langle\cdot\rangle$, we rigorously show that it is $O(L^{-d})$. In fact, we give a characterization of the leading-order error term for both strategies, and argue that even in the isotropic case it is generically non-degenerate. We carry out the rigorous analysis in the convenient setting of ensembles $\langle\cdot\rangle$ of Gaussian type with integrable covariance, which allow for a straightforward periodization and which make the Price theorem and the Malliavin calculus available for optimal stochastic estimates of correctors.
△ Less
Submitted 27 May, 2022;
originally announced May 2022.
-
A diagram-free approach to the stochastic estimates in regularity structures
Authors:
Pablo Linares,
Felix Otto,
Markus Tempelmayr,
Pavlos Tsatsoulis
Abstract:
In this paper, we explore the version of Hairer's regularity structures based on a greedier index set than trees, as introduced by Otto, Sauer, Smith and Weber. More precisely, we construct and stochastically estimate the renormalized model avoiding the use of Feynman diagrams but still in a fully automated, i. e. inductive way. This is carried out for a class of quasi-linear parabolic PDEs driven…
▽ More
In this paper, we explore the version of Hairer's regularity structures based on a greedier index set than trees, as introduced by Otto, Sauer, Smith and Weber. More precisely, we construct and stochastically estimate the renormalized model avoiding the use of Feynman diagrams but still in a fully automated, i. e. inductive way. This is carried out for a class of quasi-linear parabolic PDEs driven by noise in the full singular but renormalizable range. We assume a spectral gap inequality on the (not necessarily Gaussian) noise ensemble. The resulting control on the variance of the model naturally complements its vanishing expectation arising from the BPHZ-choice of renormalization. We capture the gain in regularity on the level of the Malliavin derivative of the model by describing it as a modelled distribution. Symmetry is an important guiding principle and built-in on the level of the renormalization Ansatz. Our approach is analytic and top-down rather than combinatorial and bottom-up.
△ Less
Submitted 30 May, 2024; v1 submitted 20 December, 2021;
originally announced December 2021.
-
There is no stationary cyclically monotone Poisson matching in 2d
Authors:
Martin Huesmann,
Francesco Mattesini,
Felix Otto
Abstract:
We show that there is no cyclically monotone stationary matching of two independent Poisson processes in dimension $d=2$. The proof combines the harmonic approximation result from \cite{GHO} with local asymptotics for the two-dimensional matching problem for which we give a new self-contained proof using martingale arguments.
We show that there is no cyclically monotone stationary matching of two independent Poisson processes in dimension $d=2$. The proof combines the harmonic approximation result from \cite{GHO} with local asymptotics for the two-dimensional matching problem for which we give a new self-contained proof using martingale arguments.
△ Less
Submitted 28 September, 2021;
originally announced September 2021.
-
Thermodynamically consistent and positivity-preserving discretization of the thin-film equation with thermal noise
Authors:
Benjamin Gess,
Rishabh S. Gvalani,
Florian Kunick,
Felix Otto
Abstract:
In micro-fluidics not only does capillarity dominate but also thermal fluctuations become important. On the level of the lubrication approximation, this leads to a quasi-linear fourth-order parabolic equation for the film height $h$ driven by space-time white noise. The gradient flow structure of its deterministic counterpart, the thin-film equation, which encodes the balance between driving capil…
▽ More
In micro-fluidics not only does capillarity dominate but also thermal fluctuations become important. On the level of the lubrication approximation, this leads to a quasi-linear fourth-order parabolic equation for the film height $h$ driven by space-time white noise. The gradient flow structure of its deterministic counterpart, the thin-film equation, which encodes the balance between driving capillary and limiting viscous forces, provides the guidance for the thermodynamically consistent introduction of fluctuations. We follow this route on the level of a spatial discretization of the gradient flow structure. Starting from an energetically conformal finite-element (FE) discretization, we point out that the numerical mobility function introduced by Grün and Rumpf can be interpreted as a discretization of the metric tensor in the sense of a mixed FE method with lumping. While this discretization was devised in order to preserve the so-called entropy estimate, we use this to show that the resulting high-dimensional stochastic differential equation (SDE) preserves pathwise and pointwise strict positivity, at least in case of the physically relevant mobility function arising from the no-slip boundary condition. As a consequence, this discretization gives rise to a consistent invariant measure, namely a discretization of the Brownian excursion (up to the volume constraint), and thus features an entropic repulsion. The price to pay over more naive discretizations is that when writing the SDE in Itô's form, which is the basis for the Euler-Mayurama time discretization, a correction term appears. To conclude, we perform various numerical experiments to compare the behavior of our discretization to that of the more naive finite difference discretization of the equation.
△ Less
Submitted 12 February, 2022; v1 submitted 13 September, 2021;
originally announced September 2021.
-
Optimal artificial boundary conditions based on second-order correctors for three dimensional random elliptic media
Authors:
Jianfeng Lu,
Felix Otto,
Lihan Wang
Abstract:
We are interested in numerical algorithms for computing the electrical field generated by a charge distribution localized on scale $\ell$ in an infinite heterogeneous medium, in a situation where the medium is only known in a box of diameter $L\gg\ell$ around the support of the charge. We propose a boundary condition that with overwhelming probability is (near) optimal with respect to scaling in t…
▽ More
We are interested in numerical algorithms for computing the electrical field generated by a charge distribution localized on scale $\ell$ in an infinite heterogeneous medium, in a situation where the medium is only known in a box of diameter $L\gg\ell$ around the support of the charge. We propose a boundary condition that with overwhelming probability is (near) optimal with respect to scaling in terms of $\ell$ and $L$, in the setting where the medium is a sample from a stationary ensemble with a finite range of dependence (set to be unity and with the assumption that $\ell \gg 1$). The boundary condition is motivated by quantitative stochastic homogenization that allows for a multipole expansion [BGO20].
This work extends [LO21], the algorithm in which is optimal in two dimension, and thus we need to take quadrupoles, next to dipoles, into account. This in turn relies on stochastic estimates of second-order, next to first-order, correctors. These estimates are provided for finite range ensembles under consideration, based on an extension of the semi-group approach of [GO15].
△ Less
Submitted 11 January, 2024; v1 submitted 3 September, 2021;
originally announced September 2021.
-
A priori bounds for quasi-linear SPDEs in the full sub-critical regime
Authors:
Felix Otto,
Jonas Sauer,
Scott Smith,
Hendrik Weber
Abstract:
This paper is concerned with quasi-linear parabolic equations driven by an additive forcing $ξ\in C^{α-2}$, in the full sub-critical regime $α\in (0,1)$. We are inspired by Hairer's regularity structures, however we work with a more parsimonious model indexed by multi-indices rather than trees. This allows us to capture additional symmetries which play a crucial role in our analysis. Assuming boun…
▽ More
This paper is concerned with quasi-linear parabolic equations driven by an additive forcing $ξ\in C^{α-2}$, in the full sub-critical regime $α\in (0,1)$. We are inspired by Hairer's regularity structures, however we work with a more parsimonious model indexed by multi-indices rather than trees. This allows us to capture additional symmetries which play a crucial role in our analysis. Assuming bounds on this model, which is modified in agreement with the concept of algebraic renormalization, we prove local a priori estimates on solutions to the quasi-linear equations modified by the corresponding counter terms.
△ Less
Submitted 26 March, 2024; v1 submitted 19 March, 2021;
originally announced March 2021.
-
Representative volume elements for matrix-inclusion composites -- a computational study on periodizing the ensemble
Authors:
Matti Schneider,
Marc Josien,
Felix Otto
Abstract:
We investigate volume-element sampling strategies for the stochastic homogenization of particle-reinforced composites and show, via computational experiments, that an improper treatment of particles intersecting the boundary of the computational cell may affect the accuracy of the computed effective properties. Motivated by recent results on a superior convergence rate of the systematic error for…
▽ More
We investigate volume-element sampling strategies for the stochastic homogenization of particle-reinforced composites and show, via computational experiments, that an improper treatment of particles intersecting the boundary of the computational cell may affect the accuracy of the computed effective properties. Motivated by recent results on a superior convergence rate of the systematic error for periodized ensembles compared to taking snapshots of ensembles, we conduct computational experiments for microstructures with circular, spherical and cylindrical inclusions and monitor the systematic errors in the effective thermal conductivity for snapshots of ensembles compared to working with microstructures sampled from periodized ensembles. We observe that the standard deviation of the apparent properties computed on microstructures sampled from the periodized ensembles decays at the scaling expected from the central limit theorem. In contrast, the standard deviation for the snapshot ensembles shows an inferior decay rate at high filler content. The latter effect is caused by additional long-range correlations that necessarily appear in particle-reinforced composites at high, industrially relevant, volume fractions. Periodized ensembles, however, appear to be less affected by these correlations. Our findings provide guidelines for working with digital volume images of material microstructures and the design of representative volume elements for computational homogenization.
△ Less
Submitted 13 March, 2021;
originally announced March 2021.
-
The structure group for quasi-linear equations via universal enveloping algebras
Authors:
Pablo Linares,
Felix Otto,
Markus Tempelmayr
Abstract:
We consider the approach of replacing trees by multi-indices as an index set of the abstract model space $\mathsf{T}$ introduced by Otto, Sauer, Smith and Weber to tackle quasi-linear singular SPDEs. We show that this approach is consistent with the postulates of regularity structures when it comes to the structure group $\mathsf{G}$. In particular, $\mathsf{G}\subset{\rm Aut}(\mathsf{T})$ arises…
▽ More
We consider the approach of replacing trees by multi-indices as an index set of the abstract model space $\mathsf{T}$ introduced by Otto, Sauer, Smith and Weber to tackle quasi-linear singular SPDEs. We show that this approach is consistent with the postulates of regularity structures when it comes to the structure group $\mathsf{G}$. In particular, $\mathsf{G}\subset{\rm Aut}(\mathsf{T})$ arises from a Hopf algebra $\mathsf{T}^+$ and a comodule $Δ\colon\mathsf{T}\rightarrow \mathsf{T}^+\otimes\mathsf{T}$. In fact, this approach, where the dual $\mathsf{T}^*$ of the abstract model space $\mathsf{T}$ naturally embeds into a formal power series algebra, allows to interpret $\mathsf{G}^*\subset{\rm Aut}(\mathsf{T}^*)$ as a Lie group arising from a Lie algebra $\mathsf{L} \subset{\rm End}(\mathsf{T}^*)$ consisting of derivations on this power series algebra. These derivations in turn are the infinitesimal generators of two actions on the space of pairs (nonlinearities, functions of space-time mod constants). These actions are shift of space-time and tilt by space-time polynomials. The Hopf algebra $\mathsf{T}^+$ arises from a coordinate representation of the universal enveloping algebra ${\rm U}(\mathsf{L})$ of the Lie algebra $\mathsf{L}$. The coordinates are determined by an underlying pre-Lie algebra structure of the derived algebra of $\mathsf{L}$. Strong finiteness properties, which are enforced by gradedness and the restrictive definition of $\mathsf{T}$, allow for this purely algebraic construction of $\mathsf{G}$. We also argue that there exist pre-Lie algebra and Hopf algebra morphisms between our structure and the tree-based one in the cases of branched rough paths (Grossman-Larson, Connes-Kreimer) and of the generalized parabolic Anderson model.
△ Less
Submitted 1 February, 2023; v1 submitted 6 March, 2021;
originally announced March 2021.
-
Differentiable Trust Region Layers for Deep Reinforcement Learning
Authors:
Fabian Otto,
Philipp Becker,
Ngo Anh Vien,
Hanna Carolin Ziesche,
Gerhard Neumann
Abstract:
Trust region methods are a popular tool in reinforcement learning as they yield robust policy updates in continuous and discrete action spaces. However, enforcing such trust regions in deep reinforcement learning is difficult. Hence, many approaches, such as Trust Region Policy Optimization (TRPO) and Proximal Policy Optimization (PPO), are based on approximations. Due to those approximations, the…
▽ More
Trust region methods are a popular tool in reinforcement learning as they yield robust policy updates in continuous and discrete action spaces. However, enforcing such trust regions in deep reinforcement learning is difficult. Hence, many approaches, such as Trust Region Policy Optimization (TRPO) and Proximal Policy Optimization (PPO), are based on approximations. Due to those approximations, they violate the constraints or fail to find the optimal solution within the trust region. Moreover, they are difficult to implement, often lack sufficient exploration, and have been shown to depend on seemingly unrelated implementation choices. In this work, we propose differentiable neural network layers to enforce trust regions for deep Gaussian policies via closed-form projections. Unlike existing methods, those layers formalize trust regions for each state individually and can complement existing reinforcement learning algorithms. We derive trust region projections based on the Kullback-Leibler divergence, the Wasserstein L2 distance, and the Frobenius norm for Gaussian distributions. We empirically demonstrate that those projection layers achieve similar or better results than existing methods while being almost agnostic to specific implementation choices. The code is available at https://git.io/Jthb0.
△ Less
Submitted 9 March, 2021; v1 submitted 22 January, 2021;
originally announced January 2021.
-
Variational methods for a singular SPDE yielding the universality of the magnetization ripple
Authors:
Radu Ignat,
Felix Otto,
Tobias Ried,
Pavlos Tsatsoulis
Abstract:
The magnetization ripple is a microstructure formed in thin ferromagnetic films. It can be described by minimizers of a nonconvex energy functional leading to a nonlocal and nonlinear elliptic SPDE in two dimensions driven by white noise, which is singular. We address the universal character of the magnetization ripple using variational methods based on $Γ$-convergence. Due to the infinite energy…
▽ More
The magnetization ripple is a microstructure formed in thin ferromagnetic films. It can be described by minimizers of a nonconvex energy functional leading to a nonlocal and nonlinear elliptic SPDE in two dimensions driven by white noise, which is singular. We address the universal character of the magnetization ripple using variational methods based on $Γ$-convergence. Due to the infinite energy of the system, the (random) energy functional has to be renormalized. Using the topology of $Γ$-convergence, we give a sense to the law of the renormalized functional that is independent of the way white noise is approximated. More precisely, this universality holds in the class of (not necessarily Gaussian) approximations to white noise satisfying the spectral gap inequality, which allows us to obtain sharp stochastic estimates. As a corollary, we obtain the existence of minimizers with optimal regularity.
△ Less
Submitted 4 October, 2023; v1 submitted 25 October, 2020;
originally announced October 2020.
-
Variational approach to regularity of optimal transport maps: general cost functions
Authors:
Felix Otto,
Maxime Prod'homme,
Tobias Ried
Abstract:
We extend the variational approach to regularity for optimal transport maps initiated by Goldman and the first author to the case of general cost functions. Our main result is an $ε$-regularity result for optimal transport maps between Hölder continuous densities slightly more quantitative than the result by De Philippis-Figalli. One of the new contributions is the use of almost-minimality: if the…
▽ More
We extend the variational approach to regularity for optimal transport maps initiated by Goldman and the first author to the case of general cost functions. Our main result is an $ε$-regularity result for optimal transport maps between Hölder continuous densities slightly more quantitative than the result by De Philippis-Figalli. One of the new contributions is the use of almost-minimality: if the cost is quantitatively close to the Euclidean cost function, a minimizer for the optimal transport problem with general cost is an almost-minimizer for the one with quadratic cost. This further highlights the connection between our variational approach and De Giorgi's strategy for $ε$-regularity of minimal surfaces.
△ Less
Submitted 9 July, 2021; v1 submitted 22 July, 2020;
originally announced July 2020.
-
The annealed Calderon-Zygmund estimate as convenient tool in quantitative stochastic homogenization
Authors:
Marc Josien,
Felix Otto
Abstract:
This article is about the quantitative homogenization theory of linear elliptic equations in divergence form with random coefficients. We derive gradient estimates on the homogenization error, i.e. on the difference between the actual solution and the two-scale expansion of the homogenized solution, both in terms of strong norms (oscillation) and weak norms (fluctuation). These estimates are optim…
▽ More
This article is about the quantitative homogenization theory of linear elliptic equations in divergence form with random coefficients. We derive gradient estimates on the homogenization error, i.e. on the difference between the actual solution and the two-scale expansion of the homogenized solution, both in terms of strong norms (oscillation) and weak norms (fluctuation). These estimates are optimal in terms of scaling in the ratio between the microscopic and the macroscopic scale. The purpose of this article is to highlight the usage of the recently introduced annealed Calderon-Zygmund (CZ) estimates in obtaining the above, previously known, error estimates. Moreover, the article provides a novel proof of these annealed CZ estimate that completely avoids quenched regularity theory, but rather relies on functional analysis. It is based on the observation that even on the level of operator norms, the Helmholtz projection is close to the one for the homogenized coefficient (for which annealed CZ estimates are easily obtained). In this article, we strive for simple proofs, and thus restrict ourselves to ensembles of coefficient fields that are local transformations of Gaussian random fields with integrable correlations and Hölder continuous realizations. As in earlier work, we use the natural objects from the general theory of homogenization, like the (potential and flux) correctors and the homogenization commutator. Both oscillation and fluctuation estimates rely on a sensitivity calculus, i.e. on estimating how sensitively the quantity of interest does depend on an infinitesimal change in the coefficient field, which is fed into the Spectral Gap inequality. In this article, the annealed CZ estimate is the only form in which elliptic regularity theory enters.
△ Less
Submitted 18 May, 2020;
originally announced May 2020.
-
Uniform energy distribution in a pattern-forming system of surface charges
Authors:
Katarina Bellova,
Antoine Julia,
Felix Otto
Abstract:
We consider a variational model for a charge density $u\in\{-1,1\}$ on a (hyper)plane, with a short-range attraction coming from the interfacial energy and a long-range repulsion coming from the electrostatic energy. This competition leads to pattern formation. We prove that the interfacial energy density is (asymptotically) equidistributed at scales large compared to the scale of the pattern. We…
▽ More
We consider a variational model for a charge density $u\in\{-1,1\}$ on a (hyper)plane, with a short-range attraction coming from the interfacial energy and a long-range repulsion coming from the electrostatic energy. This competition leads to pattern formation. We prove that the interfacial energy density is (asymptotically) equidistributed at scales large compared to the scale of the pattern. We follow the strategy laid out in [G. Alberti, R. Choksi, F. Otto, Uniform energy distribution for an isoperimetric problem with long-range interactions, J. A.M.S.]. The challenge comes from the reduced screening capabilities of surface charges compared to the volume charges considered in that paper.
△ Less
Submitted 2 April, 2020;
originally announced April 2020.
-
Sharp boundary $\varepsilon$-regularity of optimal transport maps
Authors:
Tatsuya Miura,
Felix Otto
Abstract:
In this paper we develop a boundary $\varepsilon$-regularity theory for optimal transport maps between bounded open sets with $C^{1,α}$-boundary. Our main result asserts sharp $C^{1,α}$-regularity of transport maps at the boundary in form of a linear estimate under certain assumptions: The main quantitative assumptions are that the local nondimensionalized transport cost is small and that the boun…
▽ More
In this paper we develop a boundary $\varepsilon$-regularity theory for optimal transport maps between bounded open sets with $C^{1,α}$-boundary. Our main result asserts sharp $C^{1,α}$-regularity of transport maps at the boundary in form of a linear estimate under certain assumptions: The main quantitative assumptions are that the local nondimensionalized transport cost is small and that the boundaries are locally almost flat in $C^{1,α}$. Our method is completely variational and builds on the recently developed interior regularity theory.
△ Less
Submitted 13 February, 2021; v1 submitted 20 February, 2020;
originally announced February 2020.
-
On the existence of the Green function for elliptic systems in divergence form
Authors:
Arianna Giunti,
Felix Otto
Abstract:
We study the existence of the Green function for an elliptic system in divergence form $-\nabla\cdot a\nabla$ in $\mathbb{R}^d$, with $d>2$. The tensor field $a=a(x)$ is only assumed to be bounded and $λ$-coercive. For almost every point $y \in \mathbb{R}^d$, the existence of a Green's function $G(a; \cdot, y)$ centered in $y$ has been proven in [J. Conlon, A. Giunti and F.Otto, "Green's function…
▽ More
We study the existence of the Green function for an elliptic system in divergence form $-\nabla\cdot a\nabla$ in $\mathbb{R}^d$, with $d>2$. The tensor field $a=a(x)$ is only assumed to be bounded and $λ$-coercive. For almost every point $y \in \mathbb{R}^d$, the existence of a Green's function $G(a; \cdot, y)$ centered in $y$ has been proven in [J. Conlon, A. Giunti and F.Otto, "Green's function for elliptic systems: Delmotte-Deuschel bounds", 2017]. In this paper, we show that the set of points $y \in \mathbb{R}^d$ for which $G(a; \cdot, y)$ does not exist has zero $p$-capacity, for an exponent $p >2$ depending only on the dimension $d$ and the ellipticity ratio of $a$.
△ Less
Submitted 6 June, 2020; v1 submitted 5 November, 2019;
originally announced November 2019.
-
The thresholding scheme for mean curvature flow and de Giorgi's ideas for minimizing movements
Authors:
Tim Laux,
Felix Otto
Abstract:
We consider the thresholding scheme and explore its connection to De Giorgi's ideas on gradient flows in metric spaces; here applied to mean curvature flow as the steepest descent of the interfacial area.
The basis of our analysis is the observation by Esedoglu and the second author that thresholding can be interpreted as a minimizing movements scheme for an energy that approximates the interfac…
▽ More
We consider the thresholding scheme and explore its connection to De Giorgi's ideas on gradient flows in metric spaces; here applied to mean curvature flow as the steepest descent of the interfacial area.
The basis of our analysis is the observation by Esedoglu and the second author that thresholding can be interpreted as a minimizing movements scheme for an energy that approximates the interfacial area.
De Giorgi's framework provides an optimal energy dissipation relation for the scheme in which we pass to the limit to derive a dissipation-based weak formulation of mean curvature flow.
Although applicable in the general setting of arbitrary networks, here we restrict ourselves to the case of a single interface, which allows for a compact, self-contained presentation.
△ Less
Submitted 24 October, 2019;
originally announced October 2019.
-
Quantitative estimates in stochastic homogenization for correlated coefficient fields
Authors:
Antoine Gloria,
Stefan Neukamm,
Felix Otto
Abstract:
This paper is about the homogenization of linear elliptic operators in divergence form with stationary random coefficients that have only slowly decaying correlations. It deduces optimal estimates of the homogenization error from optimal growth estimates of the (extended) corrector. In line with the heuristics, there are transitions at dimension $d=2$, and for a correlation-decay exponent $β=2$; w…
▽ More
This paper is about the homogenization of linear elliptic operators in divergence form with stationary random coefficients that have only slowly decaying correlations. It deduces optimal estimates of the homogenization error from optimal growth estimates of the (extended) corrector. In line with the heuristics, there are transitions at dimension $d=2$, and for a correlation-decay exponent $β=2$; we capture the correct power of logarithms coming from these two sources of criticality.
The decay of correlations is sharply encoded in terms of a multiscale logarithmic Sobolev inequality (LSI) for the ensemble under consideration --- the results would fail if correlation decay were encoded in terms of an $α$-mixing condition. Among other ensembles popular in modelling of random media, this class includes coefficient fields that are local transformations of stationary Gaussian fields.
The optimal growth of the corrector $φ$ is derived from bounding the size of spatial averages $F=\int g\cdot\nablaφ$ of its gradient. This in turn is done by a (deterministic) sensitivity estimate of $F$, that is, by estimating the functional derivative $\frac{\partial F}{\partial a}$ of $F$ w.~r.~t.~the coefficient field $a$. Appealing to the LSI in form of concentration of measure yields a stochastic estimate on $F$. The sensitivity argument relies on a large-scale Schauder theory for the heterogeneous elliptic operator $-\nabla\cdot a\nabla$. The treatment allows for non-symmetric $a$ and for systems like linear elasticity.
△ Less
Submitted 12 October, 2019;
originally announced October 2019.
-
Quantitative linearization results for the Monge-Ampère equation
Authors:
Michael Goldman,
Martin Huesmann,
Felix Otto
Abstract:
This paper is about quantitative linearization results for the Monge-Ampère equation with rough data. We develop a large-scale regularity theory and prove that if a measure $μ$ is close to the Lebesgue measure in Wasserstein distance at all scales, then the displacement of the macroscopic optimal coupling is quantitatively close at all scales to the gradient of the solution of the corresponding Po…
▽ More
This paper is about quantitative linearization results for the Monge-Ampère equation with rough data. We develop a large-scale regularity theory and prove that if a measure $μ$ is close to the Lebesgue measure in Wasserstein distance at all scales, then the displacement of the macroscopic optimal coupling is quantitatively close at all scales to the gradient of the solution of the corresponding Poisson equation. The main ingredient we use is a harmonic approximation result for the optimal transport plan between arbitrary measures. This is used in a Campanato iteration which transfers the information through the scales.
△ Less
Submitted 30 April, 2021; v1 submitted 23 May, 2019;
originally announced May 2019.
-
Numerical study in stochastic homogenization for elliptic PDEs: convergence rate in the size of representative volume elements
Authors:
Venera Khoromskaia,
Boris N. Khoromskij,
Felix Otto
Abstract:
We describe the numerical scheme for the discretization and solution of 2D elliptic equations with strongly varying piecewise constant coefficients arising in the stochastic homogenization of multiscale composite materials. An efficient stiffness matrix generation scheme based on assembling the local Kronecker product matrices is introduced. The resulting large linear systems of equations are solv…
▽ More
We describe the numerical scheme for the discretization and solution of 2D elliptic equations with strongly varying piecewise constant coefficients arising in the stochastic homogenization of multiscale composite materials. An efficient stiffness matrix generation scheme based on assembling the local Kronecker product matrices is introduced. The resulting large linear systems of equations are solved by the preconditioned CG iteration with a convergence rate that is independent of the grid size and the variation in jumping coefficients (contrast). Using this solver we numerically investigate the convergence of the Representative Volume Element (RVE) method in stochastic homogenization that extracts the effective behavior of the random coefficient field. Our numerical experiments confirm the asymptotic convergence rate of systematic error and standard deviation in the size of RVE rigorously established in [6]. The asymptotic behavior of covariances of the homogenized matrix in the form of a quartic tensor is also studied numerically. Our approach allows laptop computation of sufficiently large number of stochastic realizations even for large sizes of the RVE.
△ Less
Submitted 28 March, 2019;
originally announced March 2019.
-
Higher-order pathwise theory of fluctuations in stochastic homogenization
Authors:
Mitia Duerinckx,
Felix Otto
Abstract:
We consider linear elliptic equations in divergence form with stationary random coefficients of integrable correlations. We characterize the fluctuations of a macroscopic observable of a solution to relative order $\frac{d}{2}$, where $d$ is the spatial dimension; the fluctuations turn out to be Gaussian. As for previous work on the leading order, this higher-order characterization relies on a pat…
▽ More
We consider linear elliptic equations in divergence form with stationary random coefficients of integrable correlations. We characterize the fluctuations of a macroscopic observable of a solution to relative order $\frac{d}{2}$, where $d$ is the spatial dimension; the fluctuations turn out to be Gaussian. As for previous work on the leading order, this higher-order characterization relies on a pathwise proximity of the macroscopic fluctuations of a general solution to those of the (higher-order) correctors, via a (higher-order) two-scale expansion injected into the homogenization commutator, thus confirming the scope of this notion. This higher-order generalization sheds a clearer light on the algebraic structure of the higher-order versions of correctors, flux correctors, two-scale expansions, and homogenization commutators. It reveals that in the same way as this algebra provides a higher-order theory for microscopic spatial oscillations, it also provides a higher-order theory for macroscopic random fluctuations, although both phenomena are not directly related. We focus on the model framework of an underlying Gaussian ensemble, which allows for an efficient use of (second-order) Malliavin calculus for stochastic estimates. On the technical side, we introduce annealed Calderón-Zygmund estimates for the elliptic operator with random coefficients, which conveniently upgrade the known quenched large-scale estimates.
△ Less
Submitted 24 October, 2019; v1 submitted 6 March, 2019;
originally announced March 2019.
-
Relaxation Augmented Free Energy Perturbation
Authors:
Ying-Chih Chiang,
Livia B. Pártay,
Guanglian Li,
Christopher Cave-Ayland,
Marley L. Samways,
Frank Otto,
Jonathan W. Essex
Abstract:
Inspired by the recent development on calculating the free energy change via a relaxation process [Nat. Phys. 14, 842 (2018)], we investigate the role of heat released in an irreversible relaxation following a large perturbation. Utilizing a derivation without microscopic reversibility, we arrive at a new free energy estimator that employs a volume term to account for missing important rare events…
▽ More
Inspired by the recent development on calculating the free energy change via a relaxation process [Nat. Phys. 14, 842 (2018)], we investigate the role of heat released in an irreversible relaxation following a large perturbation. Utilizing a derivation without microscopic reversibility, we arrive at a new free energy estimator that employs a volume term to account for missing important rare events. Applications to harmonic oscillators and particle insertion in Lennard-Jones fluid agree well with the (numerical) exact solutions. Our study hence suggests an alternative interpretation to the insufficient sampling problem in free energy calculations.
△ Less
Submitted 29 September, 2019; v1 submitted 3 November, 2018;
originally announced November 2018.
-
A large-scale regularity theory for the Monge-Ampere equation with rough data and application to the optimal matching problem
Authors:
Michael Goldman,
Martin Huesmann,
Felix Otto
Abstract:
The aim of this paper is to obtain quantitative bounds for solutions to the optimal matching problem in dimension two. These bounds show that up to a logarithmically divergent shift, the optimal transport maps are close to be the identity at every scale. These bounds allow us to pass to the limit as the system size goes to infinity and construct a locally optimal coupling between the Lebesgue meas…
▽ More
The aim of this paper is to obtain quantitative bounds for solutions to the optimal matching problem in dimension two. These bounds show that up to a logarithmically divergent shift, the optimal transport maps are close to be the identity at every scale. These bounds allow us to pass to the limit as the system size goes to infinity and construct a locally optimal coupling between the Lebesgue measure and the Poisson point process which retains the stationarity properties of the Poisson point process only at the level of second-order differences. Our quantitative bounds are obtained through a Campanato iteration scheme based on a deterministic and a stochastic ingredient. The deterministic part, which can be seen as our main contribution, is a regularity result for Monge-Ampère equations with rough right-hand side. Since we believe that it could be useful in other contexts, we prove it for general space dimensions. The stochastic part is a concentration result for the optimal matching problem which builds on previous work by Ambrosio, Stra and Trevisan.
△ Less
Submitted 28 August, 2018;
originally announced August 2018.
-
Robustness of the pathwise structure of fluctuations in stochastic homogenization
Authors:
Mitia Duerinckx,
Antoine Gloria,
Felix Otto
Abstract:
We consider a linear elliptic system in divergence form with random coefficients and study the random fluctuations of large-scale averages of the field and the flux of the solution operator. In the context of the random conductance model, we developed in a previous work a theory of fluctuations based on the notion of homogenization commutator: we proved that the two-scale expansion of this special…
▽ More
We consider a linear elliptic system in divergence form with random coefficients and study the random fluctuations of large-scale averages of the field and the flux of the solution operator. In the context of the random conductance model, we developed in a previous work a theory of fluctuations based on the notion of homogenization commutator: we proved that the two-scale expansion of this special quantity is accurate at leading order in the fluctuation scaling when averaged on large scales (as opposed to the two-scale expansion of the solution operator taken separately) and that the large-scale fluctuations of the field and the flux of the solution operator can be recovered from those of the commutator. This implies that the large-scale fluctuations of the commutator of the corrector drive all other large-scale fluctuations to leading order, which we refer to as the pathwise structure of fluctuations in stochastic homogenization. In the present contribution we extend this result in two directions: we treat continuum elliptic (possibly non-symmetric) systems and allow for strongly correlated coefficient fields (Gaussian-like with a covariance function that can display an arbitrarily slow algebraic decay at infinity). Our main result shows in this general setting that the two-scale expansion of the homogenization commutator is still accurate to leading order when averaged on large scales, which illustrates the robustness of the pathwise structure of fluctuations.
△ Less
Submitted 23 October, 2019; v1 submitted 31 July, 2018;
originally announced July 2018.
-
Toward a quantitative theory of the hydrodynamic limit
Authors:
Deniz Dizdar,
Georg Menz,
Felix Otto,
Tianqi Wu
Abstract:
This article provides non-trivial technical ingredients for the article "The quantitative hydrodynamic limit of the Kawasaki dynamics" by the same authors. In that work a quantitative version of the hydrodynamic limit is deduced using a refinement of the two-scale approach. In this work, we deduce the strict convexity of the coarse-grained Hamiltonian, a uniform logarithmic Sobolev inequality for…
▽ More
This article provides non-trivial technical ingredients for the article "The quantitative hydrodynamic limit of the Kawasaki dynamics" by the same authors. In that work a quantitative version of the hydrodynamic limit is deduced using a refinement of the two-scale approach. In this work, we deduce the strict convexity of the coarse-grained Hamiltonian, a uniform logarithmic Sobolev inequality for the canonical ensemble and the convergence of the coarse-grained Hamiltonian to the macroscopic free energy of the system. We deduce those results following an approach developed by Grunewald, Otto, Villani and Westdickenberg. Because in our setting the associated coarse-graining operator is non-local, the arguments are much more subtle and need additional ingredients like the Brascamp-Lieb inequality and a multivariate local central-limit theorem.
△ Less
Submitted 27 July, 2018; v1 submitted 25 July, 2018;
originally announced July 2018.
-
The quantitative hydrodynamic limit of the Kawasaki dynamics
Authors:
Deniz Dizdar,
Georg Menz,
Felix Otto,
Tianqi Wu
Abstract:
We derive for the first time in the literature a rate of convergence in the hydrodynamic limit of the Kawasaki dynamics for a one-dimensional lattice system. We use an adaptation of the two-scale approach. The main difference to the original two-scale approach is that the observables on the mesoscopic level are described by a projection onto splines of second order, and not by a projection onto pi…
▽ More
We derive for the first time in the literature a rate of convergence in the hydrodynamic limit of the Kawasaki dynamics for a one-dimensional lattice system. We use an adaptation of the two-scale approach. The main difference to the original two-scale approach is that the observables on the mesoscopic level are described by a projection onto splines of second order, and not by a projection onto piecewise constant functions. This allows us to use a more natural definition of the mesoscopic dynamics, which yields a better rate of convergence than the original two-scale approach.
△ Less
Submitted 27 July, 2018; v1 submitted 25 July, 2018;
originally announced July 2018.
-
Optimal $L^1$-type relaxation rates for the Cahn-Hilliard equation on the line
Authors:
Felix Otto,
Sebastian Scholtes,
Maria G. Westdickenberg
Abstract:
In this paper we derive optimal algebraic-in-time relaxation rates to the kink for the Cahn-Hilliard equation on the line. We assume that the initial data have a finite distance---in terms of either a first moment or the excess mass---to a kink profile and capture the decay rate of the energy and the perturbation. Our tools include Nash-type inequalities, duality arguments, and Schauder estimates.
In this paper we derive optimal algebraic-in-time relaxation rates to the kink for the Cahn-Hilliard equation on the line. We assume that the initial data have a finite distance---in terms of either a first moment or the excess mass---to a kink profile and capture the decay rate of the energy and the perturbation. Our tools include Nash-type inequalities, duality arguments, and Schauder estimates.
△ Less
Submitted 10 July, 2019; v1 submitted 7 June, 2018;
originally announced June 2018.