BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Daniel Matthes

Daniel Matthes appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.APmath.NAmath.DSmath.PRmath-phmath.MPmath.OCcond-mat.quant-gasmath.DGmath.FANumerical Analysisquant-ph

Catalog footprint

What is connected

20works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2021arXiv

Barycenters for the Hellinger--Kantorovich distance over $\mathbb{R}^d$

We study the barycenter of the Hellinger--Kantorovich metric over non-negative measures on compact, convex subsets of $\mathbb{R}^d$. The article establishes existence, uniqueness (under suitable assumptions) and equivalence between a coupled-two-marginal and a multi-marginal formulation. We analyze the HK barycenter between Dirac measures in detail, and find that it differs substantially from the Wasserstein barycenter by exhibiting a local `clustering' behaviour, depending on the length scale of the input measures. In applications it makes sense to simultaneously consider all choices of this scale, leading to a 1-parameter family of barycenters. We demonstrate the usefulness of this family by analyzing point clouds sampled from a mixture of Gaussians and inferring the number and location of the underlying Gaussians.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Construction of a two-phase flow with singular energy by gradient flow methods

We prove the existence of weak solutions to a system of two diffusion equations that are coupled by a pointwise volume constraint. The time evolution is given by gradient dynamics for a free energy functional. Our primary example is a model for the demixing of polymers, the corresponding energy is the one of Flory, Huggins and deGennes. Due to the non-locality in the equations, the dynamics considered here is qualitatively different from the one found in the formally related Cahn-Hilliard equations. Our angle of attack is from the theory of optimal mass transport, that is, we consider the evolution equations for the two components as two gradient flows in the Wasserstein distance with one joint energy functional that has the volume constraint built in. The main difference to our previous work arXiv:1712.06446 is the nonlinearity of the energy density in the gradient part, which becomes singular at the interface between pure and mixed phases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Lagrangian schemes for Wasserstein gradient flows

This paper reviews different numerical methods for specific examples of Wasserstein gradient flows: we focus on nonlinear Fokker-Planck equations,but also discuss discretizations of the parabolic-elliptic Keller-Segel model and of the fourth order thin film equation. The methods under review are of Lagrangian nature, that is, the numerical approximations trace the characteristics of the underlying transport equation rather than solving the evolution equation for the mass density directly. The two main approaches are based on integrating the equation for the Lagrangian maps on the one hand, and on solution of coupled ODEs for individual mass particles on the other hand.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A fully discrete variational scheme for solving nonlinear Fokker-Planck equations in higher space dimensions

We introduce a novel spatio-temporal discretization for nonlinear Fokker-Planck equations on the multi-dimensional unit cube. This discretization is based on two structural properties of these equations: the first is the representation as a gradient flow of an entropy functional in the $L^2$-Wasserstein metric, the second is the Lagrangian nature, meaning that solutions can be written as the push forward transformation of the initial density under suitable flow maps. The resulting numerical scheme is entropy diminishing and mass conserving. Further, the scheme is weakly stable, which allows us to prove convergence under certain regularity assumptions. Finally, we present results from numerical experiments in space dimension $d=2$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Computing Coherent Sets using the Fokker-Planck Equation

We perform a numerical approximation of coherent sets in finite-dimensional smooth dynamical systems by computing singular vectors of the transfer operator for a stochastically perturbed flow. This operator is obtained by solution of a discretized Fokker-Planck equation. For numerical implementation, we employ spectral collocation methods and an exponential time differentiation scheme. We experimentally compare our approach to the more classical method by Ulam that is based on integration of the exact transfer operator.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Discrete approximation of the minimizing movement scheme for evolution equations of Wasserstein gradient flow type with nonlinear mobility

We propose a fully discrete variational scheme for nonlinear evolution equations with gradient flow structure on the space of finite Radon measures on an interval with respect to a generalized version of the Wasserstein distance with nonlinear mobility. Our scheme relies on a spatially discrete approximation of the semi-discrete (in time) minimizing movement scheme for gradient flows. Performing a finite-volume discretization of the continuity equation appearing in the definition of the distance, we obtain a finite-dimensional convex minimization problem usable as an iterative scheme. We prove that solutions to the spatially discrete minimization problem converge to solutions of the spatially continuous original minimizing movement scheme using the theory of $Γ$-convergence, and hence obtain convergence to a weak solution of the evolution equation in the continuous-time limit if the minimizing movement scheme converges. We illustrate our result with numerical simulations for several second- and fourth-order equations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Existence of solutions for a class of fourth order cross-diffusion systems of gradient flow type

This article is concerned with the existence and the long time behavior of weak solutions to certain coupled systems of fourth-order degenerate parabolic equations of gradient flow type. The underlying metric is a Wasserstein-like transportation distance for vector-valued functions, with nonlinear mobilities in each component. Under the hypothesis of (flat) convexity of the driving free energy functional, weak solutions are constructed by means of the variational minimizing movement scheme for metric gradient flows. The essential regularity estimates are derived by variational methods.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Constructing solutions to the Björling problem for isothermic surfaces by structure preserving discretization

In this article, we study an analog of the Björling problem for isothermic surfaces (that are more general than minimal surfaces): given a real analytic curve $γ$ in ${\mathbb R}^3$, and two analytic non-vanishing orthogonal vector fields $v$ and $w$ along $γ$, find an isothermic surface that is tangent to $γ$ and that has $v$ and $w$ as principal directions of curvature. We prove that solutions to that problem can be obtained by constructing a family of discrete isothermic surfaces (in the sense of Bobenko and Pinkall) from data that is sampled along $γ$, and passing to the limit of vanishing mesh size. The proof relies on a rephrasing of the Gauss-Codazzi-system as analytic Cauchy problem and an in-depth-analysis of its discretization which is induced from the geometry of discrete isothermic surfaces. The discrete-to-continuous limit is carried out for the Christoffel and the Darboux transformations as well.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Convergence of a fully discrete variational scheme for a thin-film equation

This paper is concerned with a rigorous convergence analysis of a fully discrete Lagrangian scheme for the Hele-Shaw flow, which is the fourth order thin-film equation with linear mobility in one space dimension. The discretization is based on the equation's gradient flow structure in the $L^2$-Wasserstein metric. Apart from its Lagrangian character --- which guarantees positivity and mass conservation --- the main feature of our discretization is that it dissipates both the Dirichlet energy and the logarithmic entropy. The interplay between these two dissipations paves the way to proving convergence of the discrete approximations to a weak solution in the discrete-to-continuous limit. Thanks to the time-implicit character of the scheme, no CFL-type condition is needed. Numerical experiments illustrate the practicability of the scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Existence of weak solutions to a class of fourth order partial differential equations with Wasserstein gradient structure

We prove the global-in-time existence of nonnegative weak solutions to a class of fourth order partial differential equations on a convex bounded domain in arbitrary spatial dimensions. Our proof relies on the formal gradient flow structure of the equation with respect to the $L^2$-Wasserstein distance on the space of probability measures. We construct a weak solution by approximation via the time-discrete minimizing movement scheme; necessary compactness estimates are derived by entropy-dissipation methods. Our theory essentially comprises the thin film and Derrida-Lebowitz-Speer-Spohn equations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Long-time behavior of a finite volume discretization for a fourth order diffusion equation

We consider a non-standard finite-volume discretization of a strongly non-linear fourth order diffusion equation on the $d$-dimensional cube, for arbitrary $d \geq 1$. The scheme preserves two important structural properties of the equation: the first is the interpretation as a gradient flow in a mass transportation metric, and the second is an intimate relation to a linear Fokker-Planck equation. Thanks to these structural properties, the scheme possesses two discrete Lyapunov functionals. These functionals approximate the entropy and the Fisher information, respectively, and their dissipation rates converge to the optimal ones in the discrete-to-continuous limit. Using the dissipation, we derive estimates on the long-time asymptotics of the discrete solutions. Finally, we present results from numerical experiments which indicate that our discretization is able to capture significant features of the complex original dynamics, even with a rather coarse spatial resolution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Optimal control of Bose-Einstein condensates in three dimensions

Ultracold gases promise many applications in quantum metrology, simulation and computation. In this context, optimal control theory (OCT) provides a versatile framework for the efficient preparation of complex quantum states. However, due to the high computational cost, OCT of ultracold gases has so far mostly been applied to one-dimensional (1D) problems. Here, we realize computationally efficient OCT of the Gross-Pitaevskii equation (GPE) to manipulate Bose-Einstein condensates in all three spatial dimensions. We study various realistic experimental applications where 1D simulations can only be applied approximately or not at all. Moreover, we provide a stringent mathematical footing for our scheme and carefully study the creation of elementary excitations and their minimization using multiple control parameters. The results are directly applicable to recent experiments and might thus be of immediate use in the ongoing effort to employ the properties of the quantum world for technological applications.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A convergent Lagrangian discretization for a nonlinear fourth order equation

A fully discrete Lagrangian scheme for numerical solution of the nonlinear fourth order DLSS equation in one space dimension is analyzed. The discretization is based on the equation's gradient flow structure in the $L^2$-Wasserstein metric. We prove that the discrete solutions are strictly positive and mass conserving. Further, they dissipate both the Fisher information and the logarithmic entropy. Numerical experiments illustrate the practicability of the scheme. Our main result is a proof of convergence of fully discrete to weak solutions in the limit of vanishing mesh size. Convergence is obtained for arbitrary non-negative initial data with finite entropy, without any CFL type condition. The key ingredient in the proof is a discretized version of the classical entropy dissipation estimate.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Asymptotic Behavior of Gradient Flows Driven by Nonlocal Power Repulsion and Attraction Potentials in One Dimension

We study the long time behavior of the Wasserstein gradient flow for an energy functional consisting of two components: particles are attracted to a fixed profile $ω$ by means of an interaction kernel $ψ_a(z)=|z|^{q_a}$,and they repel each other by means of another kernel $ψ_r(z)=|z|^{q_r}$. We focus on the case of one space dimension and assume that $1\le q_r\le q_a\le 2$. Our main result is that the flow converges to an equilibrium if either $q_r<q_a$ or $1\le q_r=q_a\le4/3$,and if the solution has the same (conserved) mass as the reference state $ω$. In the cases $q_r=1$ and $q_r=2$, we are able to discuss the behavior for different masses as well, and we explicitly identify the equilibrium state, which is independent of the initial condition. Our proofs heavily use the inverse distribution function of the solution.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Exponential convergence to equilibrium in a coupled gradient flow system modelling chemotaxis

We study a system of two coupled nonlinear parabolic equations. It constitutes a variant of the Keller-Segel model for chemotaxis, i.e. it models the behaviour of a population of bacteria that interact by means of a signalling substance. We assume an external confinement for the bacteria and a nonlinear dependency of the chemotactic drift on the signalling substance concentration. We perform an analysis of existence and long-time behaviour of solutions based on the underlying gradient flow structure of the system. The result is that, for a wide class of initial conditions, weak solutions exist globally in time and converge exponentially fast to the unique stationary state under suitable assumptions on the convexity of the confinement and the strength of the coupling.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Transport distances and geodesic convexity for systems of degenerate diffusion equations

We introduce Wasserstein-like dynamical transport distances between vector-valued densities on the real line. The mobility function from the scalar theory is replaced by a mobility matrix, that is subject to positivity and concavity conditions. Our primary motivation is to cast certain systems of nonlinear parabolic evolution equations into the variational framework of gradient flows. In the first part of the paper, we investigate the structural properties of the new class of distances like geodesic completeness. The second part is devoted to the identification of $λ$-geodesically convex functionals and their $λ$-contractive gradient flows. One of our results is a generalized McCann condition for geodesic convexity of the internal energy. In the third part, the existence of weak solutions to a certain class of degenerate diffusion systems with drift is shown. Even if the underlying energy function is not geodesically convex w.r.t. our new distance, the construction of a weak solution is still possible using de Giorgi's minimizing movement scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Infinite energy solutions to inelastic homogeneous Boltzmann equation

This paper is concerned with the existence, shape and dynamical stability of infinite-energy equilibria for a general class of spatially homogeneous kinetic equations in space dimensions $d \geq 3$. Our results cover in particular Bobylëv's model for inelastic Maxwell molecules. First, we show under certain conditions on the collision kernel, that there exists an index $α\in(0,2)$ such that the equation possesses a nontrivial stationary solution, which is a scale mixture of radially symmetric $α$-stable laws. We also characterize the mixing distribution as the fixed point of a smoothing transformation. Second, we prove that any transient solution that emerges from the NDA of some (not necessarily radial symmetric) $α$-stable distribution converges to an equilibrium. The key element of the convergence proof is an application of the central limit theorem to a representation of the transient solution as a weighted sum of i.i.d. random vectors.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cahn-Hilliard and Thin Film equations with nonlinear mobility as gradient flows in weighted-Wasserstein metrics

In this paper, we establish a novel approach to proving existence of non-negative weak solutions for degenerate parabolic equations of fourth order, like the Cahn-Hilliard and certain thin film equations. The considered evolution equations are in the form of a gradient flow for a perturbed Dirichlet energy with respect to a Wasserstein-like transport metric, and weak solutions are obtained as curves of maximal slope. Our main assumption is that the mobility of the particles is a concave function of their spatial density. A qualitative difference of our approach to previous ones is that essential properties of the solution - non-negativity, conservation of the total mass and dissipation of the energy - are automatically guaranteed by the construction from minimizing movements in the energy landscape.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Curves of steepest descent are entropy solutions for a class of degenerate convection-diffusion equations

We consider a nonlinear degenerate convection-diffusion equation with inhomogeneous convection and prove that its entropy solutions in the sense of Kružkov are obtained as the - a posteriori unique - limit points of the JKO variational approximation scheme for an associated gradient flow in the $L^2$-Wasserstein space. The equation lacks the necessary convexity properties which would allow to deduce well-posedness of the initial value problem by the abstract theory of metric gradient flows. Instead, we prove the entropy inequality directly by variational methods and conclude uniqueness by doubling of the variables.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Homogeneous kinetic equations for probabilistic linear collisions in multiple space dimensions

We analyze the convergence to equilibrium in a family of Kac-like kinetic equations in multiple space dimensions. These equations describe the change of the velocity distribution in a spatially homogeneous gas due to binary collisions between the particles. We consider a general linear mechanism for the exchange of the particles' momenta, with interaction coefficients that are random matrices with a distribution that is {independent} of the velocities of the colliding particles. Applying a synthesis of probabilistic methods and Fourier analysis, we are able to identify sufficient conditions for the existence and uniqueness of a stationary state, we characterize this stationary state as a mixture of Gaussian distributions, and we prove equilibration of transient solutions under minimal hypotheses on the initial conditions. In particular, we are able to classify the high-energy tails of the stationary distribution, which might be of Pareto type. We also discuss several examples to which our theory applies, among them models with a non-symmetric stationary state.

0 citations0 reviewsNo code linkedOpen access