Mon, 11 Oct 2004
15:45
DH 3rd floor SR

Joe Doob (1910-2004)

Professor N H Bingham
(University of Sheffield)
Abstract

Joe Doob, who died recently aged 94, was the last survivor of the

founding fathers of probability. Doob was best known for his work on

martingales, and for his classic book, Stochastic Processes (1953).

The talk will combine an appreciation of Doob's work and legacy with

reminiscences of Doob the man. (I was fortunate to be a colleague of

Doob from 1975-6, and to get to know him well during that year.)

Following Doob's passing, the mantle of greatest living probabilist

descends on the shoulders of Kiyosi Ito (b. 1915), alas now a sick

man.

Mon, 11 Oct 2004
14:15
DH 3rd floor SR

Stochastic individual processes and approximations in the Darwinian evolution

Professor Sylvie Meleard
(Universite Paris 10)
Abstract

We are interested in a microscopic stochastic description of a

population of discrete individuals characterized by one adaptive

trait. The population is modeled as a stochastic point process whose

generator captures the probabilistic dynamics over continuous time of

birth, mutation and death, as influenced by each individual's trait

values, and interactions between individuals. An offspring usually

inherits the trait values of her progenitor, except when a mutation

causes the offspring to take an instantaneous mutation step at birth

to new trait values. Once this point process is in place, the quest

for tractable approximations can follow different mathematical paths,

which differ in the normalization they assume (taking limit on

population size , rescaling time) and in the nature of the

corresponding approximation models: integro or integro-differential

equations, superprocesses. In particular cases, we consider the long

time behaviour for the stochastic or deterministic models.

Tue, 21 Sep 2004
11:30
DH 3rd floor SR

Theory and simulation of the shielding of emitting dust particles

Gian Luca Delzanno
Abstract

The role of electron emission (either thermionic, secondary or

photoelectric) in charging an object immersed in a plasma is

investigated, both theoretically and numerically.

In fact, recent work [1] has shown how electron emission can

fundamentally affect the shielding potential around the object. In

particular, depending on the physical parameters of the system (that

were chosen such to correspond to common experimental conditions), the

shielding potential can develop an attractive potential well.

The conditions for the formation of the well will be reviewed, based

on a theoretical model of electron emission from the

grain. Furthermore, simulations will be presented regarding specific

laboratory, space and astrophysical applications.

[1] G.L. Delzanno, G. Lapenta, M. Rosenberg, Phys. Rev.

Lett., 92, 035002 (2004).

Thu, 16 Sep 2004
12:00
DH 3rd floor SR

Quantifying Damage: Comparing Models to Data

Prof Donald L Turcotte
(University of California)
Abstract
Time delays are associated with rock fracture and earthquakes. The delay associated with the initiation of a single fracture can be attributed to stress corrosion and a critical stress intensity factor [1]. Usually, however, the fracture of a brittle material, such as rock, results from the coalescence and growth of micro cracks. Another example of time delays in rock is the systematic delay before the occurrence of earthquake aftershocks. There is also a systematic time delay associated with rate-and-state friction. One important question is whether these time delays are related. Another important question is whether the time delays are thermally activated. In many cases systematic scaling laws apply to the time delays. An example is Omori92s law for the temporal decay of after shock activity. Experiments on the fracture of fiber board panels, subjected instantaneously to a load show a systematic power-law decrease in the delay time to failure as a function of the difference between the applied stress and a yield stress [2,3]. These experiments also show a power-law increase in energy associated with acoustic emissions prior to rupture. The frequency-strength statistics of the acoustic emissions also satisfy the power-law Gutenberg-Richter scaling. Damage mechanics and dynamic fibre-bundle models provide an empirical basis for understanding the systematic time delays in rock fracture and seismicity [4-7]. We show that these approachesgive identical results when applied to fracture, and explain the scaling obtained in the fibre board experiments. These approaches also give Omori92s type law. The question of precursory activation prior to rock bursts and earthquakes is also discussed. [1] Freund, L. B. 1990. Dynamic Fracture Mechanics, Cambridge University Press, Cambridge.20
[2] Guarino, A., Garcimartin, A., and Ciliberto, S. 1998. An experimental test of the critical behaviour of fracture precursors. Eur. Phys. J.; B6:13-24.20
[3] Guarino, A., Ciliberto, S., and Garcimartin, A. 1999. Failure time and micro crack nucleation. Europhys. Lett.; 47: 456.20
[4] Kachanov, L. M. 1986. Introduction to Continuum Damage Mechanics, Martinus Nijhoff, Dordrecht, Netherlands.20
[5] Krajcinovic, D. 1996. Damage Mechanics, Elsevier, Amsterdam.20
[6] Turcotte, D. L., Newman, W. I., and Shcherbakov, R. 2002. Micro- and macroscopic models of rock fracture, Geophys. J. Int.; 152: 718-728.
[7] Shcherbakov, R. and Turcotte, D. L. 2003. Damage and self-similarity in fracture. Theor. and Appl. Fracture Mech.; 39: 245-258.
Fri, 18 Jun 2004
14:15
DH 3rd floor SR

Analytic Approximation to Loss Distributions of Heterogeneous Portfolios

Harry Zheng
(Imperial College, London)
Abstract

In this talk we discuss the analytic approximation to the loss

distribution of large conditionally independent heterogeneous portfolios. The

loss distribution is approximated by the expectation of some normal

distributions, which provides good overall approximation as well as tail

approximation. The computation is simple and fast as only numerical

integration is needed. The analytic approximation provides an excellent

alternative to some well-known approximation methods. We illustrate these

points with examples, including a bond portfolio with correlated default risk

and interest rate risk. We give an analytic expression for the expected

shortfall and show that VaR and CVaR can be easily computed by solving a

linear programming problem where VaR is the optimal solution and CVaR is the

optimal value.

Thu, 17 Jun 2004

14:00 - 15:00
Comlab

Generating good meshes and inverting good matrices

Prof Gilbert Strang
(MIT)
Abstract

An essential first step in many problems of numerical analysis and

computer graphics is to cover a region with a reasonably regular mesh.

We describe a short MATLAB code that begins with a "distance function"

to describe the region: $d(x)$ is the distance to the boundary

(with d

Tue, 15 Jun 2004

14:00 - 15:00
Rutherford Appleton Laboratory, nr Didcot

Fast and high quality display of large relational information with an introduction to recent advances in mathematica

Dr Yifan Hu
(Wolfram Research)
Abstract

The talk will start with an introduction to recent development in Mathematica, with emphasis on numerical computing. This will be followed by a discussion of graph drawing algorithms for the display of relational information, in particular force directed algorithms. The talk will show that by employing multilevel approach and octree data structure, it is possible to achieve fast display of very large relational information, without compromising the quality.

Mon, 14 Jun 2004
14:15
DH 3rd floor SR

Completing Stochastic Volatility Models with Variance Swaps

Chris Potter
(Oxford)
Abstract

Complete stochastic volatility models provide prices and

hedges. There are a number of complete models which jointly model an

underlying and one or more vanilla options written on it (for example

see Lyons, Schonbucher, Babbar and Davis). However, any consistent

model describing the volatility of options requires a complex

dependence of the volatility of the option on its strike. To date we

do not have a clear approach to selecting a model for the volatility

of these options