For some time now, Solvency II requires that insurance companies calculate minimum capital requirements to face the risk of insolvency, either in accordance with the Standard Formula or using a full or partial Internal Model. An Internal Model must be based on a market-consistent valuation of assets and liabilities at a 1-year time span, where a real-world probabilistic structure is used for the first year of projection. In this paper, we describe the major risks of a non-life insurance company, i.e. the non-life underwriting risk and market risk, and their interactions, focusing on the non-life premium risk, equity risk, and interest rate risk. This analysis is made using some well-known stochastic models in the financial-actuarial literature and practical insurance business, i.e. the Collective Risk Model for non-life premium risk, the Geometric Brownian Motion for equity risk, and a real-world version of the G2++ Model for interest rate risk, where parameters are calibrated on current and real market data. Finally, we illustrate a case study on a single-line and a multi-line insurance company in order to see how the risk drivers behave in both a stand-alone and an aggregate framework.

We derive closed-form solutions to some discounted optimal stopping problems related to the perpetual American cancellable dividend-paying put and call option pricing problems in an extension of the Black–Merton–Scholes model. The cancellation times are assumed to occur when the underlying risky asset price process hits some unobservable random thresholds. The optimal stopping times are shown to be the first times at which the asset price reaches stochastic boundaries depending on the current values of its running maximum and minimum processes. The proof is based on the reduction of the original optimal stopping problems to the associated free-boundary problems and the solution of the latter problems by means of the smooth-fit and modified normal-reflection conditions. We show that the optimal stopping boundaries are characterised as the maximal and minimal solutions of certain first-order nonlinear ordinary differential equations.

We present closed-form solutions to some discounted optimal stopping problems for the running maximum of a geometric Brownian motion with payoffs switching according to the dynamics of a continuous-time Markov chain with two states. The proof is based on the reduction of the original problems to the equivalent free-boundary problems and the solution of the latter problems by means of the smooth-fit and normal-reflection conditions. We show that the optimal stopping boundaries are determined as the maximal solutions of the associated two-dimensional systems of first-order nonlinear ordinary differential equations. The obtained results are related to the valuation of real switching lookback options with fixed and floating sunk costs in the Black–Merton–Scholes model.

In two recent papers (Maccone 2013, 2014) as well as in the book (Maccone 2012), this author described the Evolution of life on Earth over the last 3.5 billion years as a lognormal stochastic process in the increasing number of living Species. In (Maccone 2012, 2013), the process used was ‘Geometric Brownian Motion’ (GBM), largely used in Financial Mathematics (Black-Sholes models). The GBM mean value, also called ‘the trend’, always is an exponential in time and this fact corresponds to the so-called ‘Malthusian growth’ typical of population genetics. In (Maccone 2014), the author made an important generalization of his theory by extending it to lognormal stochastic processes having an arbitrary trend mL(t), rather than just a simple exponential trend as the GBM have.

The author named ‘Evo-SETI’ (Evolution and SETI) his theory inasmuch as it may be used not only to describe the full evolution of life on Earth from RNA to modern human societies, but also the possible evolution of life on exoplanets, thus leading to SETI, the current Search for ExtraTerrestrial Intelligence. In the Evo-SETI Theory, the life of a living being (let it be a cell or an animal or a human or a Civilization of humans or even an ET Civilization) is represented by a b-lognormal, i.e. a lognormal probability density function starting at a precise instant b (‘birth’) then increasing up to a peak-time p, then decreasing to a senility-time s (the descending inflexion point) and then continuing as a straight line down to the death-time d (‘finite b-lognormal’).

1. (1) Having so said, the present paper describes the further mathematical advances made by this author in 2014–2015, and is divided in two halves: Part One, devoted to new mathematical results about the History of Civilizations as b-lognormals, and

2. (2) Part Two, about the applications of the Evo-SETI Theory to the Molecular Clock, well known to evolutionary geneticists since 50 years: the idea is that our EvoEntropy grows linearly in time just as the molecular clock.

   1. (a) Summarizing the new results contained in this paper: In Part One, we start from the History Formulae already given in (Maccone 2012, 2013) and improve them by showing that it is possible to determine the b-lognormal not only by assigning its birth, senility and death, but rather by assigning birth, peak and death (BPD Theorem: no assigned senility). This is precisely what usually happens in History, when the life of a VIP is summarized by giving birth time, death time, and the date of the peak of activity in between them, from which the senility may then be calculated (approximately only, not exactly). One might even conceive a b-scalene (triangle) probability density just centred on these three points (b, p, d) and we derive the relevant equations. As for the uniform distribution between birth and death only, that is clearly the minimal description of someone's life, we compare it with both the b-lognormal and the b-scalene by comparing the Shannon Entropy of each, which is the measure of how much information each of them conveys. Finally we prove that the Central Limit Theorem (CLT) of Statistics becomes a new ‘E-Pluribus-Unum’ Theorem of the Evo-SETI Theory, giving formulae by which it is possible to find the b-lognormal of the History of a Civilization C if the lives of its Citizens Ci are known, even if only in the form of birth and death for the vast majority of the Citizens.

   2. (b) In Part Two, we firstly prove the crucial Peak-Locus Theorem for any given trend mL(t) and not just for the GBM exponential. Then we show that the resulting Evo-Entropy grows exactly linearly in time if the trend is the exponential GMB trend.

   3. (c) In addition, three Appendixes (online) with all the relevant mathematical proofs are attached to this paper. They are written in the Maxima language, and Maxima is a symbolic manipulator that may be downloaded for free from the web.

In conclusion, this paper further increases the huge mathematical spectrum of applications of the Evo-SETI Theory to prepare Humans for the first Contact with an Extra-Terrestrial Civilization.

Financial contracts with options that allow the holder to extend the contract maturity by paying an additional fixed amount have found many applications in finance. Closed-form solutions for the price of these options have appeared in the literature for the case when the contract for the underlying asset follows a geometric Brownian motion with constant interest rate, volatility and nonnegative dividend yield. In this paper, option price is derived for the case of the underlying asset that follows a geometric Brownian motion with time-dependent drift and volatility, which is more important for real life applications. The option price formulae are derived for the case of a drift that includes nonnegative or negative dividend. The latter yields a solution type that is new to the literature. A negative dividend corresponds to a negative foreign interest rate for foreign exchange options, or storage costs for commodity options. It may also appear in pricing options with transaction costs or real options, where the drift is larger than the interest rate.

In a series of recent papers and in a book, this author put forward a mathematical model capable of embracing the search for extra-terrestrial intelligence (SETI), Darwinian Evolution and Human History into a single, unified statistical picture, concisely called Evo-SETI. The relevant mathematical tools are:

1. (1) Geometric Brownian motion (GBM), the stochastic process representing evolution as the stochastic increase of the number of species living on Earth over the last 3.5 billion years. This GBM is well known in the mathematics of finances (Black–Sholes models). Its main features are that its probability density function (pdf) is a lognormal pdf, and its mean value is either an increasing or, more rarely, decreasing exponential function of the time.

2. (2) The probability distributions known as b-lognormals, i.e. lognormals starting at a certain positive instant b>0 rather than at the origin. These b-lognormals were then forced by us to have their peak value located on the exponential mean-value curve of the GBM (Peak-Locus theorem). In the framework of Darwinian Evolution, the resulting mathematical construction was shown to be what evolutionary biologists call Cladistics.

3. (3) The (Shannon) entropy of such b-lognormals is then seen to represent the ‘degree of progress’ reached by each living organism or by each big set of living organisms, like historic human civilizations. Having understood this fact, human history may then be cast into the language of b-lognormals that are more and more organized in time (i.e. having smaller and smaller entropy, or smaller and smaller ‘chaos’), and have their peaks on the increasing GBM exponential. This exponential is thus the ‘trend of progress’ in human history.

4. (4) All these results also match with SETI in that the statistical Drake equation (generalization of the ordinary Drake equation to encompass statistics) leads just to the lognormal distribution as the probability distribution for the number of extra-terrestrial civilizations existing in the Galaxy (as a consequence of the central limit theorem of statistics).

5. (5) But the most striking new result is that the well-known ‘Molecular Clock of Evolution’, namely the ‘constant rate of Evolution at the molecular level’ as shown by Kimura's Neutral Theory of Molecular Evolution, identifies with growth rate of the entropy of our Evo-SETI model, because they both grew linearly in time since the origin of life.

6. (6) Furthermore, we apply our Evo-SETI model to lognormal stochastic processes other than GBMs. For instance, we provide two models for the mass extinctions that occurred in the past: (a) one based on GBMs and (b) the other based on a parabolic mean value capable of covering both the extinction and the subsequent recovery of life forms.

7. (7) Finally, we show that the Markov & Korotayev (2007, 2008) model for Darwinian Evolution identifies with an Evo-SETI model for which the mean value of the underlying lognormal stochastic process is a cubic function of the time.

We extend the classical empirical interpolation method [M. Barrault, Y. Maday, N.C.Nguyen and A.T. Patera, An empirical interpolation method: application to efficientreduced-basis discretization of partial differential equations. Compt. Rend. Math.Anal. Num. 339 (2004) 667–672] to a weighted empiricalinterpolation method in order to approximate nonlinear parametric functions with weightedparameters, e.g. random variables obeying various probabilitydistributions. A priori convergence analysis is provided for the proposedmethod and the error bound by Kolmogorov N-width is improved from the recent work [Y.Maday, N.C. Nguyen, A.T. Patera and G.S.H. Pau, A general, multipurpose interpolationprocedure: the magic points. Commun. Pure Appl. Anal. 8(2009) 383–404]. We apply our method to geometric Brownian motion, exponentialKarhunen–Loève expansion and reduced basis approximation of non-affine stochastic ellipticequations. We demonstrate its improved accuracy and efficiency over the empiricalinterpolation method, as well as sparse grid stochastic collocation method.

This note contains two main results. (i) (Discrete time) Suppose that S is a martingale whose marginal laws agree with a geometric simple random walk. (In financial terms, let S be a risk-neutral asset price and suppose that the initial option prices agree with the Cox-Ross-Rubinstein binomial tree model.) Then S is a geometric simple random walk. (ii) (Continuous time) Suppose that S=S0eσ X-σ2〈 X〉/2 is a continuous martingale whose marginal laws agree with a geometric Brownian motion. (In financial terms, let S be a risk-neutral asset price and suppose that the initial option prices agree with the Black-Scholes model with volatility σ>0.) Then there exists a Brownian motion W such that Xt=Wt+o(t1/4+ ε) as t↑∞ for any ε> 0.

In this paper we propose a new mathematical model capable of merging Darwinian Evolution, Human History and SETI into a single mathematical scheme:

(1) Darwinian Evolution over the last 3.5 billion years is defined as one particular realization of a certain stochastic process called Geometric Brownian Motion (GBM). This GBM yields the fluctuations in time of the number of species living on Earth. Its mean value curve is an increasing exponential curve, i.e. the exponential growth of Evolution.

(2) In 2008 this author provided the statistical generalization of the Drake equation yielding the number N of communicating ET civilizations in the Galaxy. N was shown to follow the lognormal probability distribution.

(3) We call “b-lognormals” those lognormals starting at any positive time b (“birth”) larger than zero. Then the exponential growth curve becomes the geometric locus of the peaks of a one-parameter family of b-lognormals: this is our way to re-define Cladistics.

(4) b-lognormals may be also be interpreted as the lifespan of any living being (a cell, or an animal, a plant, a human, or even the historic lifetime of any civilization). Applying this new mathematical apparatus to Human History, leads to the discovery of the exponential progress between Ancient Greece and the current USA as the envelope of all b-lognormals of Western Civilizations over a period of 2500 years.

(5) We then invoke Shannon's Information Theory. The b-lognormals' entropy turns out to be the index of “development level” reached by each historic civilization. We thus get a numerical estimate of the entropy difference between any two civilizations, like the Aztec-Spaniard difference in 1519.

(6) In conclusion, we have derived a mathematical scheme capable of estimating how much more advanced than Humans an Alien Civilization will be when the SETI scientists will detect the first hints about ETs.

Burn-in is a widely used engineering procedure useful for eliminating ‘weak’ items and consequently improving the quality of remaining items. The quality of items can be measured via various performance characteristics. In the present paper we develop new performance criteria for the burn-in method. Our criteria not only take into account the reliability of an item, they also incorporate covariates.

Suppose that {Xs, 0 ≤ s ≤ T} is an m-dimensional geometric Brownian motion with drift, μ is a bounded positive Borel measure on [0,T], and ϕ : ℝm → [0,∞) is a (weighted) lq(ℝm)-norm, 1 ≤ q ≤ ∞. The purpose of this paper is to study the distribution and the moments of the random variable Y given by the Lp(μ)-norm, 1 ≤ p ≤ ∞, of the function s ↦ ϕ(Xs), 0 ≤ s ≤ T. By using various geometric inequalities in Wiener space, this paper gives upper and lower bounds for the distribution function of Y and proves that the distribution function is log-concave and absolutely continuous on every open subset of the distribution's support. Moreover, the paper derives tail probabilities, presents sharp moment inequalities, and shows that Y is indetermined by its moments. The paper will also discuss the so-called moment-matching method for the pricing of Asian-styled basket options.

The maximality principle [6] is shown to be valid in some examples of discounted optimal stopping problems for the maximum process. In each of these examples explicit formulas for the value functions are derived and the optimal stopping times are displayed. In particular, in the framework of the Black-Scholes model, the fair prices of two lookback options with infinite horizon are calculated. The main aim of the paper is to show that in each considered example the optimal stopping boundary satisfies the maximality principle and that the value function can be determined explicitly.

We study strong asymptotic properties of two types of integral functionals of geometric stochastic processes. These integral functionals are of interest in financial modelling, yielding various option pricings, annuities, etc., by appropriate selection of the processes in their respective integrands. We show that under fairly general conditions on the latter processes the logs of the integral functionals themselves asymptotically behave like appropriate sup functionals of the processes in the exponents of their respective integrands. We illustrate the possible use and applications of these strong invariance theorems by listing and elaborating on several examples.

Explicit formulas are found for the payoff and the optimal stopping strategy of the optimal stopping problem supτE (max0≤t≤τXt − c τ), where X = (Xt)t≥0 is geometric Brownian motion with drift μ and volatility σ > 0, and the supremum is taken over all stopping times for X. The payoff is shown to be finite, if and only if μ < 0. The optimal stopping time is given by τ* = inf {t > 0 | Xt = g* (max0≤t≤sXs)} where s ↦ g*(s) is the maximal solution of the (nonlinear) differential equation under the condition 0 < g(s) < s, where Δ = 1 − 2μ / σ2 and K = Δ σ2 / 2c. The estimate is established g*(s) ∼ ((Δ − 1) / K Δ)1 / Δs1−1/Δ as s → ∞. Applying these results we prove the following maximal inequality: where τ may be any stopping time for X. This extends the well-known identity E (supt>0Xt) = 1 − (σ 2 / 2 μ) and is shown to be sharp. The method of proof relies upon a smooth pasting guess (for the Stephan problem with moving boundary) and the Itô–Tanaka formula (being applied two-dimensionally). The key point and main novelty in our approach is the maximality principle for the moving boundary (the optimal stopping boundary is the maximal solution of the differential equation obtained by a smooth pasting guess). We think that this principle is by itself of theoretical and practical interest.