Researcher profile

Martin Hutzenthaler

Martin Hutzenthaler contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
28works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

28 published item(s)

preprint2022arXiv

A path-dependent stochastic Gronwall inequality and strong convergence rate for stochastic functional differential equations

We derive a stochastic Gronwall lemma with suprema over the paths in the upper bound of the assumed affine-linear growth assumption. This allows applications to Itô processes with coefficients which depend on earlier time points such as stochastic delay equations or Euler-type approximations of stochastic differential equations. We apply our stochastic Gronwall lemma with path-suprema to stochastic functional differential equations and prove a strong convergence rate for coefficient functions which depend on path-suprema.

preprint2022arXiv

Costly defense traits in structured populations

We propose a model for the dynamics of frequencies of a costly defense trait. More precisely, we consider Lotka-Volterra-type models involving a prey (or host) population consisting of two types and a predator (or parasite) population, where one type of prey individuals - modeling carriers of a defense trait - is more effective in defending against the predators but has a weak reproductive disadvantage. Under certain assumptions we prove that the relative frequency of these defenders in the total prey population converges to spatially structured Wright-Fisher diffusions with frequency-dependent migration rates. For the many-demes limit (mean-field approximation) hereof, we show that the defense trait goes to fixation/extinction if and only if the selective disadvantage is smaller/larger than an explicit function of the ecological model parameters.

preprint2022arXiv

Deep neural networks overcome the curse of dimensionality in the numerical approximation of semilinear partial differential equations

We prove that deep neural networks are capable of approximating solutions of semilinear Kolmogorov PDE in the case of gradient-independent, Lipschitz-continuous nonlinearities, while the required number of parameters in the networks grow at most polynomially in both dimension $d \in \mathbb{N}$ and prescribed reciprocal accuracy $\varepsilon$. Previously, this has only been proven in the case of semilinear heat equations.

preprint2022arXiv

Multilevel Picard approximations for high-dimensional decoupled forward-backward stochastic differential equations

Backward stochastic differential equations (BSDEs) appear in numeruous applications. Classical approximation methods suffer from the curse of dimensionality and deep learning-based approximation methods are not known to converge to the BSDE solution. Recently, Hutzenthaler et al. (arXiv:2108.10602) introduced a new approximation method for BSDEs whose forward diffusion is Brownian motion and proved that this method converges with essentially optimal rate without suffering from the curse of dimensionality. The central object of this article is to extend this result to general forward diffusions. The main challenge is that we need to establish convergence in temporal-spatial Hölder norms since the forward diffusion cannot be sampled exactly in general.

preprint2022arXiv

Multilevel Picard approximations of high-dimensional semilinear partial differential equations with locally monotone coefficient functions

The full history recursive multilevel Picard approximation method for semilinear parabolic partial differential equations (PDEs) is the only method which provably overcomes the curse of dimensionality for general time horizons if the coefficient functions and the nonlinearity are globally Lipschitz continuous and the nonlinearity is gradient-independent. In this article we extend this result to locally monotone coefficient functions. Our results cover a range of semilinear PDEs with polynomial coefficient functions.

preprint2022arXiv

On moments and strong local Hölder regularity of solutions of stochastic differential equations and of their spatial derivative processes

Spatial differentiability of solutions of stochastic differential equations (SDEs) is a classical question in stochastic analysis. The case of coefficients with globally Lipschitz continuous derivatives is well understood in the literature. Counterexamples with smooth and bounded coefficients demonstrate that the non-globally Lipschitz case is more subtle. In this article we establish conditions, including a suitable local monotonicity property, which provide existence of continuously differentiable solutions of SDEs, moment estimates and strong local Hölder regularity.

preprint2022arXiv

On the speed of convergence of Picard iterations of backward stochastic differential equations

It is a well-established fact in the scientific literature that Picard iterations of backward stochastic differential equations with globally Lipschitz continuous nonlinearity converge at least exponentially fast to the solution. In this paper we prove that this convergence is in fact at least square-root factorially fast. We show for one example that no higher convergence speed is possible in general. Moreover, if the nonlinearity is $z$-independent, then the convergence is even factorially fast. Thus we reveal a phase transition in the speed of convergence of Picard iterations of backward stochastic differential equations.

preprint2022arXiv

Stopped Brownian-increment tamed Euler method

In this article we propose a new explicit Euler-type approximation method for stochastic differential equations (SDEs). In this method, Brownian increments in the recursion of the Euler method are replaced by suitable bounded functions of the Brownian increments. We prove strong convergence rate one-half for a large class of SDEs with polynomial coefficient functions whose local monotonicity constant grows at most like the logarithm of a Lyapunov-type function.

preprint2022arXiv

Strong convergence rate of Euler-Maruyama approximations in temporal-spatial Hölder-norms

Classical approximation results for stochastic differential equations analyze the $L^p$-distance between the exact solution and its Euler-Maruyama approximations. In this article we measure the error with temporal-spatial Hölder-norms. Our motivation for this are multigrid approximations of the exact solution viewed as a function of the starting point. We establish the classical strong convergence rate $0.5$ with respect to temporal-spatial Hölder-norms if the coefficient functions have bounded derivatives of first and second order.

preprint2022arXiv

The Kolmogorov backward equation for stochastic Burgers equations and for stochastic 2D-Navier-Stokes equations

In this book we establish under suitable assumptions the uniqueness and existence of viscosity solutions of Kolmogorov backward equations for stochastic partial differential equations (SPDEs). In addition, we show that this solution is the semigroup of the corresponding SPDE. This generalizes the Feynman-Kac formula to SPDEs and establishes a link between solutions of Kolmogorov equations and solutions of the corresponding SPDEs. In contrast to the literature we only assume that the nonlinear part of the drift is Lipschitz continuous on bounded sets (and not globally Lipschitz continuous) and we allow the diffusion coefficient to be degenerate and non-constant. In the last part of this book we apply our results to stochastic Burgers equations and to stochastic 2-D Navier-Stokes equations.

preprint2021arXiv

Full history recursive multilevel Picard approximations for ordinary differential equations with expectations

We consider ordinary differential equations (ODEs) which involve expectations of a random variable. These ODEs are special cases of McKean-Vlasov stochastic differential equations (SDEs). A plain vanilla Monte Carlo approximation method for such ODEs requires a computational cost of order $\varepsilon^{-3}$ to achieve a root-mean-square error of size $\varepsilon$. In this work we adapt recently introduced full history recursive multilevel Picard (MLP) algorithms to reduce this computational complexity. Our main result shows for every $δ>0$ that the proposed MLP approximation algorithm requires only a computational effort of order $\varepsilon^{-(2+δ)}$ to achieve a root-mean-square error of size $\varepsilon$.

preprint2021arXiv

Multilevel Picard approximations for McKean-Vlasov stochastic differential equations

In the literatur there exist approximation methods for McKean-Vlasov stochastic differential equations which have a computational effort of order $3$. In this article we introduce full-history recursive multilevel Picard (MLP) approximations for McKean-Vlasov stochastic differential equations. We prove that these MLP approximations have computational effort of order $2+$ which is essentially optimal in high dimensions.

preprint2020arXiv

Differentiability of semigroups of stochastic differential equations with Hölder-continuous diffusion coefficients

Differentiability of semigroups is useful for many applications. Here we focus on stochastic differential equations whose diffusion coefficient is the square root of a differentiable function but not differentiable itself. For every $m\in\{0,1,2\}$ we establish an upper bound for a $C^m$-norm of the semigroup of such a diffusion in terms of the $C^m$-norms of the drift coefficient and of the squared diffusion coefficient. The constants in our upper bound are often dimension-independent. Our estimates are thus suitable for analyzing certain high-dimensional and infinite-dimensional degenerate stochastic differential equations.

preprint2020arXiv

Overcoming the curse of dimensionality in the numerical approximation of semilinear parabolic partial differential equations

For a long time it is well-known that high-dimensional linear parabolic partial differential equations (PDEs) can be approximated by Monte Carlo methods with a computational effort which grows polynomially both in the dimension and in the reciprocal of the prescribed accuracy. In other words, linear PDEs do not suffer from the curse of dimensionality. For general semilinear PDEs with Lipschitz coefficients, however, it remained an open question whether these suffer from the curse of dimensionality. In this paper we partially solve this open problem. More precisely, we prove in the case of semilinear heat equations with gradient-independent and globally Lipschitz continuous nonlinearities that the computational effort of a variant of the recently introduced multilevel Picard approximations grows polynomially both in the dimension and in the reciprocal of the required accuracy.

preprint2020arXiv

Strong convergence rates on the whole probability space for space-time discrete numerical approximation schemes for stochastic Burgers equations

The main result of this article establishes strong convergence rates on the whole probability space for explicit space-time discrete numerical approximations for a class of stochastic evolution equations with possibly non-globally monotone coefficients such as stochastic Burgers equations with additive trace-class noise. The key idea in the proof of our main result is (i) to bring the classical Alekseev-Gröbner formula from deterministic analysis into play and (ii) to employ uniform exponential moment estimates for the numerical approximations.

preprint2019arXiv

A stochastic Gronwall inequality and applications to moments, strong completeness, strong local Lipschitz continuity, and perturbations

There are numerous applications of the classical (deterministic) Gronwall inequality. Recently, Michael Scheutzow discovered a stochastic Gronwall inequality which provides upper bounds for $p$-th moments, $p\in(0,1)$, of the supremum of nonnegative scalar continuous processes which satisfy a linear integral inequality. In this article we complement this with upper bounds for $p$-th moments, $p\in[2,\infty)$, of the supremum of general Itô processes which satisfy a suitable one-sided affine-linear growth condition. As example applications, we improve known results on strong local Lipschitz continuity in the starting point of solutions of stochastic differential equations (SDEs), on (exponential) moment estimates for SDEs, on strong completeness of SDEs, and on perturbation estimates for SDEs.

preprint2019arXiv

Convergence in Hölder norms with applications to Monte Carlo methods in infinite dimensions

We show that if a sequence of piecewise affine linear processes converges in the strong sense with a positive rate to a stochastic process which is strongly Hölder continuous in time, then this sequence converges in the strong sense even with respect to much stronger Hölder norms and the convergence rate is essentially reduced by the Hölder exponent. Our first application hereof establishes pathwise convergence rates for spectral Galerkin approximations of stochastic partial differential equations. Our second application derives strong convergence rates of multilevel Monte Carlo approximations of expectations of Banach space valued stochastic processes.

preprint2019arXiv

Overcoming the curse of dimensionality in the numerical approximation of Allen-Cahn partial differential equations via truncated full-history recursive multilevel Picard approximations

One of the most challenging problems in applied mathematics is the approximate solution of nonlinear partial differential equations (PDEs) in high dimensions. Standard deterministic approximation methods like finite differences or finite elements suffer from the curse of dimensionality in the sense that the computational effort grows exponentially in the dimension. In this work we overcome this difficulty in the case of reaction-diffusion type PDEs with a locally Lipschitz continuous coervice nonlinearity (such as Allen-Cahn PDEs) by introducing and analyzing truncated variants of the recently introduced full-history recursive multilevel Picard approximation schemes.

preprint2019arXiv

Propagation of chaos and the many-demes limit for weakly interacting diffusions in the sparse regime

Propagation of chaos is a well-studied phenomenon and shows that weakly interacting diffusions may become independent as the system size converges to infinity. Most of the literature focuses on the case of exchangeable systems where all involved diffusions have the same distribution and are "of the same size". In this paper, we analyze the case where only a few diffusions start outside of an accessible trap. Our main result shows that in this "sparse regime" the system of weakly interacting diffusions converges in distribution to a forest of excursions from the trap. In particular, initial independence propagates in the limit and results in a forest of independent trees.

preprint2017arXiv

Multi-level Picard approximations of high-dimensional semilinear parabolic differential equations with gradient-dependent nonlinearities

Parabolic partial differential equations (PDEs) and backward stochastic differential equations (BSDEs) have a wide range of applications. In particular, high-dimensional PDEs with gradient-dependent nonlinearities appear often in the state-of-the-art pricing and hedging of financial derivatives. In this article we prove that semilinear heat equations with gradient-dependent nonlinearities can be approximated under suitable assumptions with computational complexity that grows polynomially both in the dimension and the reciprocal of the accuracy.

preprint2017arXiv

Strong convergence of full-discrete nonlinearity-truncated accelerated exponential Euler-type approximations for stochastic Kuramoto-Sivashinsky equations

This article introduces and analyzes a new explicit, easily implementable, and full discrete accelerated exponential Euler-type approximation scheme for additive space-time white noise driven stochastic partial differential equations (SPDEs) with possibly non-globally monotone nonlinearities such as stochastic Kuramoto-Sivashinsky equations. The main result of this article proves that the proposed approximation scheme converges strongly and numerically weakly to the solution process of such an SPDE. Key ingredients in the proof of our convergence result are a suitable generalized coercivity-type condition, the specific design of the accelerated exponential Euler-type approximation scheme, and an application of Fernique's theorem.

preprint2016arXiv

Exponential integrability properties of numerical approximation processes for nonlinear stochastic differential equations

Exponential integrability properties of numerical approximations are a key tool for establishing positive rates of strong and numerically weak convergence for a large class of nonlinear stochastic differential equations. It turns out that well-known numerical approximation processes such as Euler-Maruyama approximations, linear-implicit Euler approximations, and some tamed Euler approximations from the literature rarely preserve exponential integrability properties of the exact solution. The main contribution of this article is to identify a class of stopped increment-tamed Euler approximations which preserve exponential integrability properties of the exact solution under minor additional assumptions on the involved functions.

preprint2013arXiv

Branching diffusions in random environment

We consider the diffusion approximation of branching processes in random environment (BPREs). This diffusion approximation is similar to and mathematically more tractable than BPREs. We obtain the exact asymptotic behavior of the survival probability. As in the case of BPREs, there is a phase transition in the subcritical regime due to different survival opportunities. In addition, we characterize the process conditioned to never go extinct and establish a backbone construction. In the strongly subcritical regime, mean offspring numbers are increased but still subcritical in the process conditioned to never go extinct. Here survival is solely due to an immortal individual, whose offspring are the ancestors of additional families. In the weakly subcritical regime, the mean offspring number is supercritical in the process conditioned to never go extinct. Thus this process survives with positive probability even if there was no immortal individual.

preprint2013arXiv

Ecological and genetic effects of introduced species on their native competitors

Species introductions to new habitats can cause a decline in the population size of competing native species and consequently also in their genetic diversity. We are interested in why these adverse effects are weak in some cases whereas in others the native species declines to the point of extinction. While the introduction rate and the growth rate of the introduced species in the new environment clearly have a positive relationship with invasion success and impact, the influence of competition is poorly understood. Here, we investigate how the intensity of interspecific competition influences the persistence time of a native species in the face of repeated and ongoing introductions of the nonnative species. We analyze two stochastic models: a model for the population dynamics of both species and a model that additionally includes the population genetics of the native species at a locus involved in its adaptation to a changing environment. Counterintuitively, both models predict that the persistence time of the native species is lowest for an intermediate intensity of competition. This phenomenon results from the opposing effects of competition at different stages of the invasion process: With increasing competition intensity more introduction events are needed until a new species can establish, but increasing competition also speeds up the exclusion of the native species by an established nonnative competitor. By comparing the ecological and the eco-genetic model, we detect and quantify a synergistic feedback between ecological and genetic effects.

preprint2013arXiv

Interacting diffusions and trees of excursions: convergence and comparison

We consider systems of interacting diffusions with local population regulation. Our main result shows that the total mass process of such a system is bounded above by the total mass process of a tree of excursions with appropriate drift and diffusion coefficients. As a corollary, this entails a sufficient, explicit condition for extinction of the total mass as time tends to infinity. On the way to our comparison result, we establish that systems of interacting diffusions with uniform migration between finitely many islands converge to a tree of excursions as the number of islands tends to infinity. In the special case of logistic branching, this leads to a duality between the tree of excursions and the solution of a McKean-Vlasov equation.

preprint2013arXiv

Supercritical branching diffusions in random environment

Supercritical branching processes in constant environment conditioned on eventual extinction are known to be subcritical branching processes. The case of random environment is more subtle. A supercritical branching diffusion in random environment (BDRE) conditioned on eventual extinction of the population is not a BDRE. However the quenched law of the population size of a supercritical BDRE conditioned on eventual extinction is equal to the quenched law of the population size of a subcritical BDRE. As a consequence, supercritical BDREs have a phase transition which is similar to a well-known phase transition of subcritical branching processes in random environment.

preprint2012arXiv

Strong convergence of an explicit numerical method for SDEs with nonglobally Lipschitz continuous coefficients

On the one hand, the explicit Euler scheme fails to converge strongly to the exact solution of a stochastic differential equation (SDE) with a superlinearly growing and globally one-sided Lipschitz continuous drift coefficient. On the other hand, the implicit Euler scheme is known to converge strongly to the exact solution of such an SDE. Implementations of the implicit Euler scheme, however, require additional computational effort. In this article we therefore propose an explicit and easily implementable numerical method for such an SDE and show that this method converges strongly with the standard order one-half to the exact solution of the SDE. Simulations reveal that this explicit strongly convergent numerical scheme is considerably faster than the implicit Euler scheme.

preprint2011arXiv

Convergence of the stochastic Euler scheme for locally Lipschitz coefficients

Stochastic differential equations are often simulated with the Monte Carlo Euler method. Convergence of this method is well understood in the case of globally Lipschitz continuous coefficients of the stochastic differential equation. The important case of superlinearly growing coefficients, however, has remained an open question. The main difficulty is that numerically weak convergence fails to hold in many cases of superlinearly growing coefficients. In this paper we overcome this difficulty and establish convergence of the Monte Carlo Euler method for a large class of one-dimensional stochastic differential equations whose drift functions have at most polynomial growth.