Due to the growing complexity of problems encountered in the financial industry and energy markets, the research on single and multiple stopping problems, and more general control problems remains an important challenge. In particular, there is an increasing demand for solutions to realworld optimal control problems which are both accurate and effective.
Highlights
New Monte Carlo procedures are developed that allow for efficient evaluation of upper bounds for callable financial products. These procedures are based on regression methods for (approximate) identification of the Doob martingale part of the corresponding Snell envelope in the context of optimal stopping and, more generally, optimal control. Unlike the industrial standard AndersenBroadie type algorithms these new approaches do not require nested simulation. In particular, new dual martingale representations were developed for general multiple stopping and control problems that straightforwardly allow for a simulation based algorithmic solution. On the other hand, it is shown that standard nested simulation algorithms, such as simulation based policy iteration and AndersenBroadie type dual algorithms, can also be significantly improved by introducing Multilevel Monte Carlo on the level of inner simulations.Publications
Monographs

D. Belomestny, J. Schoenmakers, Advanced SimulationBased Methods for Optimal Stopping and Control: With Applications in Finance, Macmillan Publishers Ltd., London, 2018, 364 pages, (Monograph Published), DOI 10.1057/9781137033512 .

J.G.M. Schoenmakers, Robust Libor Modelling and Pricing of Derivative Products, Chapman & Hall CRC Press, 2005, 202 pages, (Monograph Published).
Articles in Refereed Journals

M. Redmann, P. Kürschner, An output error bound for timelimited balanced truncation, Systems & Control Letters, 121 (2018), pp. 16, DOI 10.1016/j.sysconle.2018.08.004 .
Abstract
When solving partial differential equations numerically, usually a high order spatial discretization is needed. Model order reduction (MOR) techniques are often used to reduce the order of spatiallydiscretized systems and hence reduce computational complexity. A particular MOR technique to obtain a reduced order model (ROM) is balanced truncation (BT). However, if one aims at finding a good ROM on a certain finite time interval only, timelimited BT (TLBT) can be a more accurate alternative. So far, no error bound on TLBT has been proved. In this paper, we close this gap in the theory by providing an output error bound for TLBT with two different representations. The performance of the error bound is then shown in several numerical experiments. 
M. Redmann, Type II balanced truncation for deterministic bilinear control systems, SIAM Journal on Control and Optimization, 56 (2018), pp. 25932612, DOI 10.1137/17M1147962 .
Abstract
When solving partial differential equations numerically, usually a high order spatial discretisation is needed. Model order reduction (MOR) techniques are often used to reduce the order of spatiallydiscretised systems and hence reduce computational complexity. A particular MOR technique to obtain a reduced order model (ROM) is balanced truncation (BT), a method which has been extensively studied for deterministic linear systems. As socalled type I BT it has already been extended to bilinear equations, an important subclass of nonlinear systems. We provide an alternative generalisation of the linear setting to bilinear systems which is called type II BT. The Gramians that we propose in this context contain information about the control. It turns out that the new approach delivers energy bounds which are not just valid in a small neighbourhood of zero. Furthermore, we provide an ℋ∞error bound which so far is not known when applying type I BT to bilinear systems. 
M. Redmann, Type II singular perturbation approximation for linear systems with Lévy noise, SIAM Journal on Control and Optimization, 56 (2018), pp. 21202158, DOI 10.1137/17M113160X .
Abstract
When solving linear stochastic partial differential equations numerically, usually a high order spatial discretisation is needed. Model order reduction (MOR) techniques are often used to reduce the order of spatiallydiscretised systems and hence reduce computational complexity. A particular MOR technique to obtain a reduced order model (ROM) is singular perturbation approximation (SPA), a method which has been extensively studied for deterministic systems. As socalled type I SPA it has already been extended to stochastic equations. We provide an alternative generalisation of the deterministic setting to linear systems with Lévy noise which is called type II SPA. It turns out that the ROM from applying type II SPA has better properties than the one of using type I SPA. In this paper, we provide new energy interpretations for stochastic reachability Gramians, show the preservation of mean square stability in the ROM by type II SPA and prove two different error bounds for type II SPA when applied to Lévy driven systems 
M. Redmann, Type II singular perturbation approximation for linear systems with Lévy noise, SIAM Journal on Control and Optimization, 56 (2018), pp. 21202158, DOI 10.1137/17M113160X .
Abstract
When solving linear stochastic partial differential equations numerically, usually a high order spatial discretisation is needed. Model order reduction (MOR) techniques are often used to reduce the order of spatiallydiscretised systems and hence reduce computational complexity. A particular MOR technique to obtain a reduced order model (ROM) is singular perturbation approximation (SPA), a method which has been extensively studied for deterministic systems. As socalled type I SPA it has already been extended to stochastic equations. We provide an alternative generalisation of the deterministic setting to linear systems with Lévy noise which is called type II SPA. It turns out that the ROM from applying type II SPA has better properties than the one of using type I SPA. In this paper, we provide new energy interpretations for stochastic reachability Gramians, show the preservation of mean square stability in the ROM by type II SPA and prove two different error bounds for type II SPA when applied to Lévy driven systems 
V. Krätschmer, M. Ladkau, R.J.A. Laeven, J.G.M. Schoenmakers, M. Stadje, Optimal stopping under uncertainty in drift and jump intensity, Mathematics of Operations Research, (2018), published online on 09.08.2018, urlhttps://doi.org/101287/moor.2017.0899., DOI 10.1287/moor.2017.0899 .
Abstract
This paper studies the optimal stopping problem in the presence of model uncertainty (ambiguity). We develop a method to practically solve this problem in a general setting, allowing for general timeconsistent ambiguity averse preferences and general payoff processes driven by jumpdiffusions. Our method consists of three steps. First, we construct a suitable Doob martingale associated with the solution to the optimal stopping problem %represented by the Snell envelope using backward stochastic calculus. Second, we employ this martingale to construct an approximated upper bound to the solution using duality. Third, we introduce backwardforward simulation to obtain a genuine upper bound to the solution, which converges to the true solution asymptotically. We analyze the asymptotic behavior and convergence properties of our method. We illustrate the generality and applicability of our method and the potentially significant impact of ambiguity to optimal stopping in a few examples. 
M. Hintermüller, M. Hinze, Ch. Kahle, T. Keil, A goaloriented dualweighted adaptive finite element approach for the optimal control of a nonsmooth CahnHilliardNavierStokes system, Optimization and Engineering. International Multidisciplinary Journal to Promote Optimization Theory & Applications in Engineering Sciences, (2018), published online on 23.06.2018, DOI 10.1007/s1108101893936 .
Abstract
This paper is concerned with the development and implementation of an adaptive solution algorithm for the optimal control of a timediscrete CahnHilliardNavierStokes system with variable densities. The free energy density associated to the CahnHilliard system incorporates the doubleobstacle potential which yields an optimal control problem for a family of coupled systems in each time instant of a variational inequality of fourth order and the NavierStokes equation. A dualweighed residual approach for goaloriented adaptive finite elements is presented which is based on the concept of Cstationarity. The overall error representation depends on primal residual weighted by approximate dual quantities and vice versa as well as various complementary mismatch errors. Details on the numerical realization of the adaptive concept and a report on numerical tests are given. 
D. Belomestny, R. Hildebrand, J.G.M. Schoenmakers, Optimal stopping via pathwise dual empirical maximisation, Applied Mathematics and Optimization. An International Journal with Applications to Stochastics, pp. published online on 08.11.2017, urlhttps://doi.org/10.1007/s0024501794549, DOI 10.1007/s0024501794549 .
Abstract
The optimal stopping problem arising in the pricing of American options can be tackled by the so called dual martingale approach. In this approach, a dual problem is formulated over the space of martingales. A feasible solution of the dual problem yields an upper bound for the solution of the original primal problem. In practice, the optimization is performed over a finitedimensional subspace of martingales. A sample of paths of the underlying stochastic process is produced by a MonteCarlo simulation, and the expectation is replaced by the empirical mean. As a rule the resulting optimization problem, which can be written as a linear program, yields a martingale such that the variance of the obtained estimator can be large. In order to decrease this variance, a penalizing term can be added to the objective function of the pathwise optimization problem. In this paper, we provide a rigorous analysis of the optimization problems obtained by adding different penalty functions. In particular, a convergence analysis implies that it is better to minimize the empirical maximum instead of the empirical mean. Numerical simulations confirm the variance reduction effect of the new approach. 
F. Dickmann, N. Schweizer, Faster comparison of stopping times by nested conditional Monte Carlo, Journal of Computational Finance, 20 (2016), (24).

D. Belomestny, F. Dickmann, T. Nagapetyan, Pricing Bermudan options via multilevel approximation methods, SIAM Journal on Financial Mathematics, ISSN 1945497X, 6 (2015), pp. 448466.
Abstract
In this article we propose a novel approach to reducing the computational complexity of various approximation methods for pricing discrete time American or Bermudan options. Given a sequence of continuation values estimates corresponding to different levels of spatial approximation, we propose a multilevel low biased estimate for the price of the option. It turns out that the resulting complexity gain can be of order ? ?1 with ? denoting the desired precision. The performance of the proposed multilevel algorithms is illustrated by a numerical example. 
D. Belomestny, M. Ladkau, J.G.M. Schoenmakers, Simulation based policy iteration for American style derivatives  A multilevel approach, SIAM ASA J. Uncertainty Quantification, 3 (2015), pp. 460483.
Abstract
This paper presents a novel approach to reduce the complexity of simulation based policy iteration methods for pricing American options. Typically, Monte Carlo construction of an improved policy gives rise to a nested simulation algorithm for the price of the American product. In this respect our new approach uses the multilevel idea in the context of the inner simulations required, where each level corresponds to a specific number of inner simulations. A thorough analysis of the crucial convergence rates in the respective multilevel policy improvement algorithm is presented. A detailed complexity analysis shows that a significant reduction in computational effort can be achieved in comparison to standard Monte Carlo based policy iteration. 
CH. Bender, J.G.M. Schoenmakers, J. Zhang, Dual representations for general multiple stopping problems, Mathematical Finance. An International Journal of Mathematics, Statistics and Financial Economics, 25 (2015), pp. 339370.
Abstract
In this paper, we study the dual representation for generalized multiple stopping problems, hence the pricing problem of general multiple exercise options. We derive a dual representation which allows for cashflows which are subject to volume constraints modeled by integer valued adapted processes and refraction periods modeled by stopping times. As such, this extends the works by Schoenmakers [2010], Bender [2011a], Bender [2011b], Aleksandrov and Hambly [2010] and Meinshausen and Hambly [2004] on multiple exercise options, which either take into consideration a refraction period or volume constraints, but not both simultaneously. We also allow more flexible cashflow structures than the additive structure in the above references. For example some exponential utility problems are covered by our setting. We supplement the theoretical results with an explicit Monte Carlo algorithm for constructing confidence intervals for the price of multiple exercise options and exemplify it by a numerical study on the pricing of a swing option in an electricity market. 
G. Milshteyn, V. Spokoiny, Construction of meanselffinancing strategies for European options under regimeswitching, SIAM Journal on Financial Mathematics, ISSN 1945497X, 5 (2014), pp. 532556.
Abstract
The paper focuses on the problem of pricing and hedging a European contingent claim for an incomplete market model, in which evolution of price processes for a saving account and stocks depends on an observable Markov chain. The pricing function is evaluated using the martingale approach. The equivalent martingale measure is introduced in a way that the Markov chain remains the historical one. Due to the Markovian structure of the considered model, the pricing function satisfies the Cauchy problem for a system of linear parabolic partial differential equations. It is shown that any European contingent claim is attainable using a generalized replicating strategy which is selffinancing in mean. For such a strategy, apart from the initial endowment and trading, some additional funds are required both stepwise at the jump moments of the Markov chain and continuously between the jump moments. The connection of the considered pricing and hedging problems with partial differential equations is very useful for computations. 
S. Balder, A. Mahayni, J.G.M. Schoenmakers, Primaldual linear Monte Carlo algorithm for multiple stopping  An application to flexible caps, Quantitative Finance, 13 (2013), pp. 10031013.
Abstract
In this paper we consider the valuation of Bermudan callable derivatives with multiple exercise rights. We present in this context a new primaldual linear Monte Carlo algorithm that allows for efficient simulation of lower and upper price bounds without using nested simulations (hence the terminology). The algorithm is essentially an extension of a primaldual Monte Carlo algorithm for standard Bermudan options proposed in Schoenmakers et al (2011), to the case of multiple exercise rights. In particular, the algorithm constructs upwardly a system of dual martingales to be plugged into the dual representation of Schoenmakers (2010). At each level the respective martingale is constructed via a backward regression procedure starting at the last exercise date. The thus constructed martingales are finally used to compute an upper price bound. At the same time, the algorithm also provides approximate continuation functions which may be used to construct a price lower bound. The algorithm is applied to the pricing of flexible caps in a Hull White (1990) model setup. The simple model choice allows for comparison of the computed price bounds with the exact price which is obtained by means of a trinomial tree implementation. As a result, we obtain tight price bounds for the considered application. Moreover, the algorithm is generically designed for multidimensional problems and is tractable to implement. 
D. Belomestny, J.G.M. Schoenmakers, F. Dickmann, Multilevel dual approach for pricing American style derivatives, Finance and Stochastics, 17 (2013), pp. 717742.
Abstract
In this article we propose a novel approach to reduce the computational complexity of the dual method for pricing American options. We consider a sequence of martingales that converges to a given target martingale and decompose the original dual representation into a sum of representations that correspond to different levels of approximation to the target martingale. By next replacing in each representation true conditional expectations with their Monte Carlo estimates, we arrive at what one may call a multilevel dual Monte Carlo algorithm. The analysis of this algorithm reveals that the computational complexity of getting the corresponding target upper bound, due to the target martingale, can be significantly reduced. In particular, it turns out that using our new approach, we may construct a multilevel version of the wellknown nested Monte Carlo algorithm of Andersen and Broadie (2004) that is, regarding complexity, virtually equivalent to a nonnested algorithm. The performance of this multilevel algorithm is illustrated by a numerical example. 
J.G.M. Schoenmakers, J. Zhang, J. Huang, Optimal dual martingales, their analysis and application to new algorithms for Bermudan products, SIAM Journal on Financial Mathematics, ISSN 1945497X, 4 (2013), pp. 86116.
Abstract
In this paper we introduce and study the concept of optimal and surely optimal dual martingales in the context of dual valuation of Bermudan options. We provide a theorem which give conditions for a martingale to be surely optimal, and a stability theorem concerning martingales which are near to be surely optimal in a sense. Guided by these theorems we develop a regression based backward construction of such a martingale in a Wiener environment. In turn this martingale may be utilized for computing upper bounds by nonnested Monte Carlo. As a byproduct, the algorithm also provides approximations to continuation values of the product, which in turn determine a stopping policy. Hence, we obtain lower bounds at the same time. The proposed algorithm is pure dual in the sense that it doesn't require an (input) approximation to the Snell envelope, is quite easy to implement, and in a numerical study we show that, regarding the computed upper bounds, it is comparable with the method of Belomestny, et. al. (2009). 
A. Mahayni, J.G.M. Schoenmakers, Minimum return guarantees with funds switching rights  An optimal stopping problem, Journal of Economic Dynamics & Control, 35 (2012), pp. 18801897.
Abstract
Recently, there is a growing trend to offer guarantee products where the investor is allowed to shift her account/investment value between multiple funds. The switching right is granted a finite number per year, i.e. it is American style with multiple exercise possibilities. In consequence, the pricing and the risk management is based on the switching strategy which maximizes the value of the guarantee put option. We analyze the optimal stopping problem in the case of one switching right within different model classes and compare the exact price with the lower price bound implied by the optimal deterministic switching time. We show that, within the class of logprice processes with independent increments, the stopping problem is solved by a deterministic stopping time if (and only if) the price process is in addition continuous. Thus, in a sense, the Black & Scholes model is the only (meaningful) pricing model where the lower price bound gives the exact price. It turns out that even moderate deviations from the Black & Scholes model assumptions give a lower price bound which is really below the exact price. This is illustrated by means of a stylized stochastic volatility model setup. 
J.G.M. Schoenmakers, A pure martingale dual for multiple stopping, Finance and Stochastics, 16 (2012), pp. 319334.
Abstract
In this paper we present a dual representation for the multiple stopping problem, hence multiple exercise options. As such it is a natural generalization of the method in Rogers (2002) and Haugh and Kogan (2004) for the standard stopping problem for American options. We consider this representation as the real dual as it is solely expressed in terms of an infimum over martingales rather than an infimum over martingales and stopping times as in Meinshausen and Hambly (2004). For the multiple dual representation we present three Monte Carlo simulation algorithms which require only one degree of nesting. 
V. Krätschmer, J.G.M. Schoenmakers, Representations for optimal stopping under dynamic monetary utility functionals, SIAM Journal on Financial Mathematics, ISSN 1945497X, 1 (2010), pp. 811832.
Abstract
In this paper we consider the optimal stopping problem for general dynamic monetary utility functionals. Sufficient conditions for the Bellman principle and the existence of optimal stopping times are provided. Particular attention is payed to representations which allow for a numerical treatment in real situations. To this aim, generalizations of standard evaluation methods like policy iteration, dual and consumption based approaches are developed in the context of general dynamic monetary utility functionals. As a result, it turns out that the possibility of a particular generalization depends on specific properties of the utility functional under consideration. 
D. Belomestny, L. Rüschendorf, M. Urusov, Optimal stopping of integral functionals and a ``noloss'' free boundary formulation, SIAM Journal of Theory of Probability and Its Applications, 54 (2010), pp. 1428.

D. Belomestny, A. Kolodko, J.G.M. Schoenmakers, Regression methods for stochastic control problems and their convergence analysis, SIAM Journal on Control and Optimization, 48 (2010), pp. 35623588.
Abstract
In this paper we develop several regression algorithms for solving general stochastic optimal control problems via Monte Carlo. This type of algorithms is particulary useful for problems with a highdimensional state space and complex dependence structure of the underlying Markov process with respect to some control. The main idea behind the algorithms is to simulate a set of trajectories under some reference measure and to use the Bellman principle combined with fast methods for approximating conditional expectations and functional optimization. Theoretical properties of the presented algorithms are investigated and the convergence to the optimal solution is proved under mild assumptions. Finally, we present numerical results for the problem of pricing a highdimensional Bermudan basket option under transaction costs in a financial market with a large investor. 
D. Belomestny, G.N. Milstein, J.G.M. Schoenmakers, Sensitivities for Bermudan options by regression methods, Decisions in Economics and Finance. A Journal of Applied Mathematics, 33 (2010), pp. 117138.
Abstract
In this article we propose several pathwise and finite difference based methods for calculating sensitivities of Bermudan options using regression methods and Monte Carlo simulation. These methods rely on conditional probabilistic representations which allows, in combination with a regression approach, an efficient simultaneous computation of sensitivities at all initial positions. Assuming that the price of a Bermudan option can be evaluated sufficiently accurate, we develop a method for constructing deltas based on least squares. We finally propose a testing procedure for assessing the performance of the developed methods. 
R.L. Loeffen, J.F. Renaud, De Finetti's optimal dividends problem with an affine penalty function at ruin, Insurance: Mathematics & Economics, 46 (2010), pp. 98108.

D. Belomestny, Ch. Bender, J.G.M. Schoenmakers, True upper bounds for Bermudan products via nonnested Monte Carlo, Mathematical Finance. An International Journal of Mathematics, Statistics and Financial Economics, 19 (2009), pp. 5371.
Abstract
We present a generic nonnested Monte Carlo procedure for computing true upper bounds for Bermudan products, given an approximation of the Snell envelope. The pleonastic “true” stresses that, by construction, the estimator is biased above the Snell envelope. The key idea is a regression estimator for the Doob martingale part of the approximative Snell envelope, which preserves the martingale property. The so constructed martingale may be employed for computing dual upper bounds without nested simulation. In general, this martingale can also be used as a control variate for simulation of conditional expectations. In this context, we develop a variance reduced version of the nested primaldual estimator (Anderson & Broadie (2004)) and nested consumption based (Belomestny & Milstein (2006)) methods . Numerical experiments indicate the efficiency of the nonnested Monte Carlo algorithm and the variance reduced nested one. 
D. Belomestny, G.N. Milstein, V. Spokoiny, Regression methods in pricing American and Bermudan options using consumption processes, Quantitative Finance, 9 (2009), pp. 315327.
Abstract
Here we develop methods for efficient pricing multidimensional discretetime American and Bermudan options by using regression based algorithms together with a new approach towards constructing upper bounds for the price of the option. Applying sample space with payoffs at the optimal stopping times, we propose sequential estimates for continuation values, values of the consumption process, and stopping times on the sample paths. The approach admits constructing both low and upper bounds for the price by Monte Carlo simulations. The methods are illustrated by pricing Bermudan swaptions and snowballs in the Libor market model. 
R.L. Loeffen, An optimal dividends problem with a terminal value for spectrally negative Lévy processes with a completely monotone jump density, Journal of Applied Probability, 46 (2009), pp. 8598.

R.L. Loeffen, An optimal dividends problem with transaction costs for spectrally negative Lévy processes, Insurance: Mathematics & Economics, 45 (2009), pp. 4148.

J. Kampen, A. Kolodko, J.G.M. Schoenmakers, Monte Carlo Greeks for financial products via approximative transition densities, SIAM Journal on Scientific Computing, 31 (2008), pp. 122.

CH. Bender, A. Kolodko, J.G.M. Schoenmakers, Enhanced policy iteration for American options via scenario selection, Quantitative Finance, 8 (2008), pp. 135146.
Abstract
In Kolodko & Schoenmakers (2004) and Bender & Schoenmakers (2004) a policy iteration was introduced which allows to achieve tight lower approximations of the price for early exercise options via a nested MonteCarlo simulation in a Markovian setting. In this paper we enhance the algorithm by a scenario selection method. It is demonstrated by numerical examples that the scenario selection can significantly reduce the number of actually performed inner simulations, and thus can heavily speed up the method (up to factor 10 in some examples). Moreover, it is shown that the modified algorithm retains the desirable properties of the original one such as the monotone improvement property, termination after a finite number of iteration steps, and numerical stability. 
A. Kolodko, J.G.M. Schoenmakers, Iterative construction of the optimal Bermudan stopping time, Finance and Stochastics, 10 (2006), pp. 2749.
Abstract
We present an iterative procedure for computing the optimal Bermudan stopping time, hence the Bermudan Snell envelope. The method produces an increasing sequence of approximations of the Snell envelope from below, which coincide with the Snell envelope after finitely many steps. Then, by duality, the method induces a convergent sequence of upper bounds as well. In a Markovian setting the presented procedure allows to calculate approximative solutions with only a few nestings of conditional expectations and is therefore tailormade for a plain Monte Carlo implementation. The method may be considered generic for all discrete optimal stopping problems. The power of the procedure is demonstrated for Bermudan swaptions in a full factor LIBOR market model. 
CH. Bender, A. Kolodko, J.G.M. Schoenmakers, Iterating cancelable snowballs and related exotics, Risk Magazine, 9 (2006), pp. 126130.
Abstract
Effective valuation procedures for callable exotics are a thorny problem. Standard methods reveal limitations in pricing manydimensional and pathdependent products, such as cancellable snowballs. Christian Bender, Anastasia Kolodko and John Schoenmakers ally these methods with their recent iterative methodology to fill the final gap. 
CH. Bender, A. Kolodko, J.G.M. Schoenmakers, Policy iteration for American options: Overview, Monte Carlo Methods and Applications, 12 (2006), pp. 347362.
Abstract
This paper is an overview of recent results by Kolodko and Schoenmakers (2006), Bender and Schoenmakers (2006) on the evaluation of options with early exercise opportunities via policy improvement. Stability is discussed and simulation results based on plain Monte Carlo estimators for conditional expectations are presented. 
CH. Bender, J.G.M. Schoenmakers, An iterative method for multiple stopping: Convergence and stability, Advances in Applied Probability, 38 (2006), pp. 729749, DOI 10.1239/aap/1158684999 .
Abstract
We present a new iterative procedure for solving the multiple stopping problem in discrete time and discuss the stability of the algorithm. The algorithm produces monotonically increasing approximations of the Snell envelope, which coincide with the Snell envelope after finitely many steps. Contrary to backward dynamic programming, the algorithm allows to calculate approximative solutions with only a few nestings of conditional expectations and is, therefore, tailormade for a plain MonteCarlo implementation. 
A. Kolodko, J.G.M. Schoenmakers, Upper bounds for Bermudan style derivatives, Monte Carlo Methods and Applications, 10 (2004), pp. 331343.
Abstract
Based on a duality approach for Monte Carlo construction of upper bounds for American/Bermudan derivatives (Rogers, Haugh & Kogan), we present a new algorithm for computing dual upper bounds in a more e?cient way. The method is applied to Bermudan swaptions in the context of a LIBOR market model, where the dual upper bound is constructed from the maximum of still alive swaptions. We give a numerical comparison with Andersen's lower bound method. 
G.N. Milstein, O. Reiss, J.G.M. Schoenmakers, A new Monte Carlo method for American options, International Journal of Theoretical and Applied Finance, 7 (2004), pp. 591614, DOI 10.1142/S0219024904002554 .
Abstract
We introduce a new Monte Carlo method for constructing the exercise boundary of an American option in a generalized BlackScholes framework. Based on a known exercise boundary, it is shown how to price and hedge the American option by Monte Carlo simulation of suitable probabilistic representations in connection with the respective parabolic boundary value problem. The method presented is supported by numerical experiments. 
R. Liptser, A.Y. Veretennikov, V. Spokoiny, FreidlinWentzell type moderate deviations for smooth processes, Markov Processes and Related Fields, 8 (2002), pp. 611636.
Contributions to Collected Editions

D. Becherer, J.G.M. Schoenmakers, E3  Stochastic simulation methods for optimal stopping and control  Towards multilevel approaches, in: MATHEON  Mathematics for Key Technologies, M. Grötschel, D. Hömberg, J. Sprekels, V. Mehrmann ET AL., eds., 1 of EMS Series in Industrial and Applied Mathematics, European Mathematical Society Publishing House, Zurich, 2014, pp. 317331.

P. Čížek, V. Spokoiny, Varying coefficient GARCH models, in: Handbook of Financial Time Series, T.G. Andersen, R.A. Davis, J.P. Kreiss, T. Mikosch, eds., Springer, Berlin, 2009, pp. 169185.

CH. Bender, A. Kolodko, J.G.M. Schoenmakers, Policy iteration method for American options, in: Proceedings of 4th Actuarial and Financial Mathematics Day, M. Vanmaele, A. De Schepper, J. Dhaene, H. Reynaerts, W. Schoutens, P. Van Goethem, eds., Royal Flemish Academy of Belgium for Sciences and Arts, Brussels, 2006, pp. 3145.
Preprints, Reports, Technical Reports

CH. Bayer, M. Redmann, J.G.M. Schoenmakers, Dynamic programming for optimal stopping via pseudoregression, Preprint no. 2532, WIAS, Berlin, 2018, DOI 10.20347/WIAS.PREPRINT.2532 .
Abstract, PDF (337 kByte)
We introduce new variants of classical regressionbased algorithms for optimal stopping problems based on computation of regression coefficients by Monte Carlo approximation of the corresponding L^{2} inner products instead of the leastsquares error functional. Coupled with new proposals for simulation of the underlying samples, we call the approach "pseudo regression". We show that the approach leads to asymptotically smaller errors, as well as less computational cost. The analysis is justified by numerical examples. 
D. Belomestny, J.G.M. Schoenmakers, V. Spokoiny, Y. Tavyrikov, Optimal stopping via deeply boosted backward regression, Preprint no. 2530, WIAS, Berlin, 2018, DOI 10.20347/WIAS.PREPRINT.2530 .
Abstract, PDF (209 kByte)
In this note we propose a new approach towards solving numerically optimal stopping problems via boosted regression based Monte Carlo algorithms. The main idea of the method is to boost standard linear regression algorithms in each backward induction step by adding new basis functions based on previously estimated continuation values. The proposed methodology is illustrated by several numerical examples from finance. 
P. Goyal, M. Redmann, Towards timelimited $H_2$optimal model order reduction, Preprint no. 2441, WIAS, Berlin, 2017, DOI 10.20347/WIAS.PREPRINT.2441 .
Abstract, PDF (363 kByte)
In order to solve partial differential equations numerically and accurately, a high order spatial discretization is usually needed. Model order reduction (MOR) techniques are often used to reduce the order of spatiallydiscretized systems and hence reduce computational complexity. A particular class of MOR techniques are H_2optimal methods such as the iterative rational Krylov subspace algorithm (IRKA) and related schemes. However, these methods are used to obtain good approximations on a infinite timehorizon. Thus, in this work, our main goal is to discuss MOR schemes for timelimited linear systems. For this, we propose an alternative timelimited H_2norm and show its connection with the timelimited Gramians. We then provide firstorder optimality conditions for an optimal reduced order model (ROM) with respect to the timelimited H_2norm. Based on these optimality conditions, we propose an iterative scheme which upon convergences aims at satisfying these conditions. Then, we analyze how far away the obtained ROM is from satisfying the optimality conditions. We test the efficiency of the proposed iterative scheme using various numerical examples and illustrate that the newly proposed iterative method can lead to a better reducedorder compared to unrestricted IRKA in the time interval of interest. 
R. Hildebrand, J.G.M. Schoenmakers, J. Zhang, F. Dickmann, Regression based duality approach to optimal control with application to hydro electricity storage, Preprint no. 2330, WIAS, Berlin, 2016, DOI 10.5072/WIAS.PREPRINT.2330 .
Abstract, PDF (341 kByte)
In this paper we consider the problem of optimal control of stochastic processes. We employ the dual martingale method brought forward in [Brown, Smith, and Sun, 2010]. The martingale constituting the solution of the dual problem is determined by linear regression within a MonteCarlo approach. We apply the solution algorithm to a model of a hydro electricity storage and production system coupled with a model of the electricity wholesale market. 
CH. Heinemann, K. Sturm, Shape optimisation for a class of semilinear variational inequalities with applications to damage models, Preprint no. 2209, WIAS, Berlin, 2016.
Abstract, PDF (590 kByte)
The present contribution investigates shape optimisation problems for a class of semilinear elliptic variational inequalities with Neumann boundary conditions. Sensitivity estimates and material derivatives are firstly derived in an abstract operator setting where the operators are defined on polyhedral subsets of reflexive Banach spaces. The results are then refined for variational inequalities arising from minimisation problems for certain convex energy functionals considered over upper obstacle sets in $H^1$. One particularity is that we allow for dynamic obstacle functions which may arise from another optimisation problems. We prove a strong convergence property for the material derivative and establish stateshape derivatives under regularity assumptions. Finally, as a concrete application from continuum mechanics, we show how the dynamic obstacle case can be used to treat shape optimisation problems for timediscretised brittle damage models for elastic solids. We derive a necessary optimality system for optimal shapes whose state variables approximate desired damage patterns and/or displacement fields.
Talks, Poster

M. Redmann, Type II singular perturbation approximation for linear systems with Levy noise, London Mathematical Society  EPSRC Durham Symposium: Model Order Reduction, Durham University, Department of Mathematical Sciences, UK, August 14, 2017.

J.G.M. Schoenmakers, Optimal stopping and control via approximative dynamic programming, Workshop on Mathematics of Deep Learning 2017, September 13  15, 2017, WIAS, Berlin, September 13, 2017.

R. Hildebrand, Periodic discrete dynamical systems and copositive matrices with circulant zero patterns, Optimization Without Borders, February 7  12, 2016, Les Houches, France, February 11, 2016.

J.G.M. Schoenmakers, Multilevel dual evaluation and multilevel policy iteration for optimal stopping/American options, Advances in Financial Mathematics, January 7  10, 2014, l'Institut Louis Bachelier, Paris, France, January 9, 2014.

R. Hildebrand, Generic fractal structure of the optimal synthesis in problems with affine multidimensional control, European Control Conference 2013, July 17  19, 2013, Eidgenoessische Technische Hochschule Zürich, Switzerland, July 19, 2013.

M. Ladkau, Multilevel policy iteration for pricing American options, 26th European Conference on Operational Research, June 30  July 4, 2013, Università La Sapienza, Rome, Italy, July 2, 2013.

M. Ladkau, Multilevel policy iteration for pricing american options, Workshop on Stochastic Models and Control, March 18  22, 2013, HumboldtUniversität zu Berlin, March 21, 2013.

M. Ladkau, Multilevel policy iteration for pricing American options, PreMoLab: MoscowBerlin Stochastic and Predictive Modeling, May 31  June 1, 2012, Russian Academy of Sciences, Institute for Information Transmission Problems (Kharkevich Institute), Moscow, June 1, 2012.

J. Zhang, Dual representations for general multiple stopping problems, Young Researchers Meeting on BSDEs, July 2  4, 2012, Oxford University, Man Institute, UK, July 2, 2012.

J. Zhang, Dual representations for general multiple stopping problems, PreMoLab: MoscowBerlin Stochastic and Predictive Modeling, May 31  June 1, 2012, Russian Academy of Sciences, Institute for Information Transmission Problems (Kharkevich Institute), Moscow, June 1, 2012.

J. Zhang, Dual representations for general multiple stopping problems, Oberseminar Finanz und Versicherungsmathematik, LudwigMaximiliansUniversität München, Mathematisches Institut, July 17, 2012.

J. Zhang, Dual representations for generalized multiple stopping problems, 10th German Probalility and Statistics Days 2012, March 6  9, 2012, Johannes Gutenberg Universität Mainz, March 7, 2012.

J.G.M. Schoenmakers, Multilevel dual approach for pricing American options, Miniworkshop CWIEUR Backward Stochastic Differential Equations (BSDE's), January 17  18, 2012, Eindhoven University of Technology, European Institute for Statistics, Probability, Stochastic Operations Research and its Applications (EURANDOM), January 18, 2012.

J.G.M. Schoenmakers, Optimal dual martingales, their analysis and application to new algorithms for Bermudan products, 10th German Probalility and Statistic Days 2012, March 6  8, 2012, Johannes Gutenberg Universität Mainz, March 7, 2012.

J.G.M. Schoenmakers , Multilevel dual approach for pricing American options, PreMoLab: MoscowBerlin Stochastic and Predictive Modeling, May 31  June 1, 2012, Russian Academy of Sciences, Institute for Information Transmission Problems (Kharkevich Institute), Moscow, June 1, 2012.

J.G.M. Schoenmakers , Multilevel primal and dual approaches for pricing American options, 21st International Symposium on Mathematical Programming (ISMP), August 20  24, 2012, Technische Universität Berlin, August 24, 2012.

M. Ladkau, A new multifactor stochastic volatility model with displacement, International Workshop on Numerical Algorithms in Computational Finance, July 20  22, 2011, Goethe Universität Frankfurt, Goethe Center for Scientific Computing (GCSC), July 21, 2011.

J. Zhang, Dual representations for general multiple stopping problems, Forschungsseminar Professor Dr. Michael Kupper, Humboldt Universität zu Berlin, December 14, 2011.

R.L. Loeffen, Some discretetime stochastic optimal control problems in finance, May 12  June 3, 2011, University of Wroclaw, Faculty of Mathematics and Computer Science, Poland.

J.G.M. Schoenmakers, New dual methods for single and multiple exercise option, Universität Ulm, Institut für Numerische Mathematik, May 27, 2011.

J.G.M. Schoenmakers, New dual methods for single and multiple exercise options, Workshop ``Quantitative Methods in Financial and Insurance Mathematics'', April 18  21, 2011, Lorentz Center, Leiden, Netherlands, April 21, 2011.

J.G.M. Schoenmakers , New dual methods for single and multiple exercise options, International Workshop on Numerical Algorithms in Computational Finance, July 20  22, 2011, Goethe Universität Frankfurt, Goethe Center for Scientific Computing, July 22, 2011.

V. Krätschmer, Representations for optimal stopping under dynamic monetary utility functionals, Leipziger Stochastik Tage, March 1  5, 2010, Universität Leipzig, Fakultät für Mathematik und Informatik, March 3, 2010.

D. Belomestny, On the rates of convergence of simulationbased optimization algorithms for optimal stopping problems, Leipziger StochastikTage, March 1  5, 2010, Universität Leipzig, Fakultät für Mathematik und Informatik, March 3, 2010.

R.L. Loeffen, Absolute ruin in the insurance risk model of OrnsteinUhlenbeck type, 24th European Conference on Operational Research (EURO XXIV LISBON), July 11  14, 2010, Universidade de Lisboa, Faculdade de Ciéncias, Portugal, July 14, 2010.

J.G.M. Schoenmakers, Advanced Libor modeling, Postbank Bonn, February 25, 2010.

J.G.M. Schoenmakers, On three innovations in financial modeling, Colloquium, University of Twente, Faculty of Electrical Engineering, Mathematics and Computer Science, Netherlands, August 24, 2010.

J.G.M. Schoenmakers, The real multiple dual, Leipziger Stochastik Tage, March 1  5, 2010, Universität Leipzig, Fakultät für Mathematik und Informatik, March 2, 2010.

B. Metzger, The discrete GrossPitaevskii model and condensation in the single particle ground state, BerlinLeipzig Seminar on Analysis and Probability Theory, November 27, 2009, Technische Universität Berlin, Institut für Mathematik, November 27, 2009.

D. Belomestny, Regression methods for stochastic control problems and their convergence analysis, Fourth General Conference on Advanced Mathematical Methods in Finance, May 4  10, 2009, University of Oslo, Norway, May 9, 2009.

R. Loeffen, De Finetti's optimal dividends problem with an affine penalty function at ruin, Université Libre de Bruxelles, Département de Mathématiques, Belgium, November 19, 2009.

CH. Meyer, Finite element error analysis for stateconstrained optimal control of the Stokes problem, 80th Annual Meeting of the International Association of Applied Mathematics and Mechanics (GAMM 2009), Section ``Optimization of Differential Equations'', February 9  13, 2009, Gdansk University of Technology, Poland, February 10, 2009.

A. Mielke, BV solutions and parametrized solutions for rateindependent systems, 80th Annual Meeting of the International Association of Applied Mathematics and Mechanics (GAMM 2009), Session ``Applied Analysis'', February 9  13, 2009, Gdansk University of Technology, Poland, February 10, 2009.

J.G.M. Schoenmakers, Monte Carlo methods for pricing of complex structured callable derivatives, RheinMain Arbeitskreis em Mathematics of Computation, Johann Wolfgang GoetheUniversität Frankfurt, Villa mboxGiersch, January 16, 2009.

J.G.M. Schoenmakers, Regression methods for stochastic control problems and their convergence analysis, Workshop ``Computational Finance'', August 10  12, 2009, Kyoto University, Faculty of Sciences, Japan, August 10, 2009.

J.G.M. Schoenmakers, Statistical and numerical methods for evaluation for financial derivates and risk, Center Days 2009 (DFG Research Center scshape Matheon), March 30  April 1, 2009, Technische Universität Berlin, March 31, 2009.

J.G.M. Schoenmakers, Monte Carlo methods for pricing of complex structured callable derivatives, 2nd International Conference on Numerical Methods for Finance, June 4  5, 2008, Institute for Numerical Computation and Analysis, Dublin, Ireland, June 4, 2008.

J.G.M. Schoenmakers, New Monte Carlo methods for pricing highdimensional callable derivatives, Conference on Numerical Methods in Finance, June 26  27, 2008, Università di Udine, Italy, June 27, 2008.

J.G.M. Schoenmakers, Pricing and hedging exotic interest rate derivatives with Monte Carlo simulation, finance master class$^rm TM$ workshop, March 19  20, 2008, Concentric Italy, Milan.

J.G.M. Schoenmakers, Regression methods for highdimensional Bermudan derivatives and stochastic control problems, Conference on Numerical Methods for American and Bermudan Options, October 17  18, 2008, Wolfgang Pauli Institute (WPI), Fakultät für Mathematik, Vienna, Austria, October 17, 2008.

J.G.M. Schoenmakers, Enhanced policy iteration via scenario selection, International Multidisciplinary Workshop on Stochastic Modeling, June 25  29, 2007, Sevilla, Spain, June 27, 2007.

J.G.M. Schoenmakers, Iterative procedures for the Bermudan stopping problem, International Multidisciplinary Workshop on Stochastic Modeling, June 25  29, 2007, Sevilla, Spain, June 26, 2007.

J.G.M. Schoenmakers, Policy iteration for American/Bermudan style derivatives, 6th Winter School on Mathematical Finance, January 22  24, 2007, CongresHotel De Werelt, Lunteren, Netherlands, January 23, 2007.

J.G.M. Schoenmakers, Robust Libor modelling and calibration, International Multidisciplinary Workshop on Stochastic Modeling, June 25  29, 2007, Sevilla, Spain, June 29, 2007.

J.G.M. Schoenmakers, True upper bounds for Bermudan style derivatives, International Multidisciplinary Workshop on Stochastic Modeling, June 25  29, 2007, Sevilla, Spain, June 28, 2007.

J.G.M. Schoenmakers, Interest rate modelling: Practical calibration and implementation techniques, June 15  16, 2006, Risk, London, UK.

J.G.M. Schoenmakers, Iterative Methoden zur Bewertung komplex strukturierter Finanzderivate mit vorzeitigen Ausübungsrechten, WIASDay, WIAS, Berlin, February 24, 2006.

J.G.M. Schoenmakers, Iterative methods for complex structured callable products, 7th GOR Workshop on Financial Optimization and Optimal Pricing Strategies, May 22  23, 2006, BASF AG, Ludwigshafen, May 23, 2006.

J.G.M. Schoenmakers, Iterative procedures for the Bermudan stopping problem, 42nd Dutch Mathematical Congress, March 27  28, 2006, Delft University of Technology, Faculty of Electrical Engineering, Mathematics and Computer Science, Netherlands, March 27, 2006.

J.G.M. Schoenmakers, Iterative construction of the optimal Bermudan stopping time, Frankfurt MathFinance Workshop ``Derivatives and risk management in theory and practice'', April 14  15, 2005, HfB  Business School of Finance and Management, Frankfurt am Main, April 14, 2005.

J.G.M. Schoenmakers, Robust Libor modelling and pricing of derivative products, Delft University of Technology, Netherlands, June 9, 2005.

J.G.M. Schoenmakers, Interactive construction of the optimal Bermudan stopping time, Eidgenössische Technische Hochschule Zürich, Institut für Mathematik, Switzerland, May 6, 2004.

J.G.M. Schoenmakers, Interest rate modelling  Practical calibration and implementation techniques, Incisive Media Events, Hong Kong, China, December 8, 2004.

J.G.M. Schoenmakers, Iterative construction of the optimal Bermudan stopping time, 2nd IASTED International Conference on Financial Engineering and Applications (FEA 2004), November 8  10, 2004, Cambridge, USA, November 8, 2004.

J.G.M. Schoenmakers, New Monte Carlo Methods for American and Bermudan style derivatives, Delft University of Technology, Faculty of Information Technology and Systems Numerical Analysis Group, Netherlands, January 30, 2004.

J.G.M. Schoenmakers, Monte Carlo methods for pricing and hedging American options, IV IMACS Seminar on Monte Carlo Methods; (Workshop Financial Models and Simulation), September 15  19, 2003, Berlin, September 19, 2003.

J.G.M. Schoenmakers, Monte Carlo simulation of Bermudan derivatives by dual upper bounds, Miniworkshop ``Risikomaße und ihre Anwendungen'', HumboldtUniversität zu Berlin, December 1, 2003.

J.G.M. Schoenmakers, Monte Carlo simulation of Bermudan derivatives by dual upper bounds, Graduiertenkolleg Angewandte Algorithmische Mathematik, Workshop on the Interface of Numerical Analysis, Optimisation and Applications, November 13  14, 2003, Technische Universität München, November 14, 2003.

J.G.M. Schoenmakers, Robust calibration of LIBOR market models, Petit Dejeuner de la Finance, November 4  5, 2003, Paris, November 5, 2003.

J. Polzehl, Structural adaptive smoothing and its applications in imaging and time series, Uppsala University, Department of Mathematics, Sweden, May 2, 2002.
External Preprints

CH. Bayer, R. Tempone, S. Wolfers, Pricing american options by exercise rate optimization, Preprint no. arXiv:1809.07300, Cornell University Library, arXiv.org, 2018.
Abstract
We present a novel method for the numerical pricing of American options based on Monte Carlo simulation and optimization of exercise strategies. Previous solutions to this problem either explicitly or implicitly determine socalled optimal emphexercise regions, which consist of points in time and space at which the option is exercised. In contrast, our method determines emphexercise rates of randomized exercise strategies. We show that the supremum of the corresponding stochastic optimization problem provides the correct option price. By integrating analytically over the random exercise decision, we obtain an objective function that is differentiable with respect to perturbations of the exercise rate even for finitely many sample paths. Starting in a neutral strategy with constant exercise rate then allows us to globally optimize this function in a gradual manner. Numerical experiments on vanilla put options in the multivariate BlackScholes model and preliminary theoretical analysis underline the efficiency of our method both with respect to the number of timediscretization steps and the required number of degrees of freedom in the parametrization of exercise rates. Finally, the flexibility of our method is demonstrated by numerical experiments on max call options in the BlackScholes model and vanilla put options in Heston model and the nonMarkovian rough Bergomi model. 
J.G.M. Schoenmakers, J. Huang, Optimal dual martingales, their analysis and application to new algorithms for Bermudan products, Preprint no. 1825944, Social Science Research Network (SSRN) Working Paper Series, 2011.

A. Kyprianou, R.L. Loeffen, J.L. Perez, Optimal control with absolutely continuous strategies for spectrally negative Levy processes, Preprint no. arXiv:1008.2363, Cornell University Library, arXiv.org, 2010.
Abstract
In the last few years there has been renewed interest in the classical control problem of de Finetti for the case that underlying source of randomness is a spectrally negative Levy process. In particular a significant step forward is made in an article of Loeffen where it is shown that a natural and very general condition on the underlying Levy process which allows one to proceed with the analysis of the associated HamiltonJacobiBellman equation is that its Levy measure is absolutely continuous, having completely monotone density. In this paper we consider de Finetti's control problem but now with the restriction that control strategies are absolutely continuous with respect to Lebesgue measure. This problem has been considered by Asmussen and Taksar, Jeanblanc and Shiryaev and Boguslavskaya in the diffusive case and Gerber and Shiu for the case of a CramerLundberg process with exponentially distributed jumps. We show the robustness of the condition that the underlying Levy measure has a completely monotone density and establish an explicit optimal strategy for this case that envelopes the aforementioned existing results. The explicit optimal strategy in question is the socalled refraction strategy. 
P. Friz, M. Caruana, H. Oberhauser, A (rough) pathwise approach to fully nonlinear stochastic partial differential equations, Report no. RICAM Report no. 200903, Johann Radon Institute for Computational and Applied Mathematics, 2009.