Adams-Type Methods for the Numerical Solution of Stochastic Ordinary Differential Equations

Stochastic differential equations (SDEs) arise fi om physical systems where the parameters describing the system can only be estimated or are subject to noise. There has been much work done recently on developing numerical methods for solving SDEs. This paper will focus on stability issues and variable stepsize implementation techniques for numerically solving SDEs effectively. (C) 2000 Elsevier Science B.V. All rights reserved.

High strong order methods for non-commutative stochastic ordinary differential equation systems and the Magnus formula

In recent years considerable attention has been paid to the numerical solution of stochastic ordinary differential equations (SODEs), as SODEs are often more appropriate than their deterministic counterparts in many modelling situations. However, unlike the deterministic case numerical methods for SODEs are considerably less sophisticated due to the difficulty in representing the (possibly large number of) random variable approximations to the stochastic integrals. Although Burrage and Burrage [High strong order explicit Runge-Kutta methods for stochastic ordinary differential equations, Applied Numerical Mathematics 22 (1996) 81-101] were able to construct strong local order 1.5 stochastic Runge-Kutta methods for certain cases, it is known that all extant stochastic Runge-Kutta methods suffer an order reduction down to strong order 0.5 if there is non-commutativity between the functions associated with the multiple Wiener processes. This order reduction down to that of the Euler-Maruyama method imposes severe difficulties in obtaining meaningful solutions in a reasonable time frame and this paper attempts to circumvent these difficulties by some new techniques. An additional difficulty in solving SODEs arises even in the Linear case since it is not possible to write the solution analytically in terms of matrix exponentials unless there is a commutativity property between the functions associated with the multiple Wiener processes. Thus in this present paper first the work of Magnus [On the exponential solution of differential equations for a linear operator, Communications on Pure and Applied Mathematics 7 (1954) 649-673] (applied to deterministic non-commutative Linear problems) will be applied to non-commutative linear SODEs and methods of strong order 1.5 for arbitrary, linear, non-commutative SODE systems will be constructed - hence giving an accurate approximation to the general linear problem. Secondly, for general nonlinear non-commutative systems with an arbitrary number (d) of Wiener processes it is shown that strong local order I Runge-Kutta methods with d + 1 stages can be constructed by evaluated a set of Lie brackets as well as the standard function evaluations. A method is then constructed which can be efficiently implemented in a parallel environment for this arbitrary number of Wiener processes. Finally some numerical results are presented which illustrate the efficacy of these approaches. (C) 1999 Elsevier Science B.V. All rights reserved.

General order conditions for stochastic Runge-Kutta methods for both commuting and non-commuting stochastic ordinary differential equation systems

In many modeling situations in which parameter values can only be estimated or are subject to noise, the appropriate mathematical representation is a stochastic ordinary differential equation (SODE). However, unlike the deterministic case in which there are suites of sophisticated numerical methods, numerical methods for SODEs are much less sophisticated. Until a recent paper by K. Burrage and P.M. Burrage (1996), the highest strong order of a stochastic Runge-Kutta method was one. But K. Burrage and P.M. Burrage (1996) showed that by including additional random variable terms representing approximations to the higher order Stratonovich (or Ito) integrals, higher order methods could be constructed. However, this analysis applied only to the one Wiener process case. In this paper, it will be shown that in the multiple Wiener process case all known stochastic Runge-Kutta methods can suffer a severe order reduction if there is non-commutativity between the functions associated with the Wiener processes. Importantly, however, it is also suggested how this order can be repaired if certain commutator operators are included in the Runge-Kutta formulation. (C) 1998 Elsevier Science B.V. and IMACS. All rights reserved.

A bound on the maximum strong order of stochastic Runge-Kutta methods for stochastic ordinary differential equations

In Burrage and Burrage [1] it was shown that by introducing a very general formulation for stochastic Runge-Kutta methods, the previous strong order barrier of order one could be broken without having to use higher derivative terms. In particular, methods of strong order 1.5 were developed in which a Stratonovich integral of order one and one of order two were present in the formulation. In this present paper, general order results are proven about the maximum attainable strong order of these stochastic Runge-Kutta methods (SRKs) in terms of the order of the Stratonovich integrals appearing in the Runge-Kutta formulation. In particular, it will be shown that if an s-stage SRK contains Stratonovich integrals up to order p then the strong order of the SRK cannot exceed min{(p + 1)/2, (s - 1)/2), p greater than or equal to 2, s greater than or equal to 3 or 1 if p = 1.

Regulating journalists? The Finkelstein Review, the Convergence Review, and News Media Regulation in Australia

This paper identifies two major forces driving change in media policy worldwide: media convergence, and renewed concerns about media ethics, with the latter seen in the U.K. Leveson Inquiry. It focuses on two major public inquiries in Australia during 2011-2012 – the Independent Media Inquiry (Finkelstein Review) and the Convergence Review – and the issues raised about future regulation of journalism and news standards. Drawing upon perspectives from media theory, it observes the strong influence of social responsibility theories of the media in the Finkelstein Review, and the adverse reaction these received from those arguing from Fourth Estate/free press perspectives, which were also consistent with the longstanding opposition of Australian newspaper proprietors to government regulation. It also discusses the approaches taken in the Convergence Review to regulating for news standards, in light of the complexities arising from media convergence. The paper concludes with consideration of the fast-changing environment in which such proposals to transform media regulation are being considered, including the crisis of news media organisation business models, as seen in Australia with major layoffs of journalists from the leading print media publications.

Improved estimates of percolation and anisotropic permeability from 3-D X-ray microtomography using stochastic analyses and visualization

X-ray microtomography (micro-CT) with micron resolution enables new ways of characterizing microstructures and opens pathways for forward calculations of multiscale rock properties. A quantitative characterization of the microstructure is the first step in this challenge. We developed a new approach to extract scale-dependent characteristics of porosity, percolation, and anisotropic permeability from 3-D microstructural models of rocks. The Hoshen-Kopelman algorithm of percolation theory is employed for a standard percolation analysis. The anisotropy of permeability is calculated by means of the star volume distribution approach. The local porosity distribution and local percolation probability are obtained by using the local porosity theory. Additionally, the local anisotropy distribution is defined and analyzed through two empirical probability density functions, the isotropy index and the elongation index. For such a high-resolution data set, the typical data sizes of the CT images are on the order of gigabytes to tens of gigabytes; thus an extremely large number of calculations are required. To resolve this large memory problem parallelization in OpenMP was used to optimally harness the shared memory infrastructure on cache coherent Non-Uniform Memory Access architecture machines such as the iVEC SGI Altix 3700Bx2 Supercomputer. We see adequate visualization of the results as an important element in this first pioneering study.

Convergence of corporate social responsibility and corporate governance in weak economies : the case of Bangladesh

The convergence of corporate social responsibility (CSR) and corporate governance (CG) has changed the corporate accountability mechanism. This has developed a socially responsible ‘corporate self-regulation’, a synthesis of governance and responsibility in the companies of strong economies. However, unlike in the strong economies, this convergence has not been visible in the companies of weak economies, where the civil society groups are unorganised, regulatory agencies are either ineffective or corrupt and the media and non-governmental organisations do not mirror the corporate conscience. Using the case of Bangladesh, this article investigates the convergence between CSR and CG in the self-regulation of companies in a less vigilant environment.

The Convergence Review and the principle of Australians' access to and opportunities for participation

The recent Australian Convergence Review’s second principle states: “Australians should have access to and opportunities for participation in a diverse mix of services, voices, views and information”. However, in failing to define its own use and understanding of the terms ‘access’ and ‘participation’ the Convergence Review exposes itself to criticism. These terms would no doubt be made unambiguously clear should the Review’s recommendations move towards policy, and this paper contributes to this discussion by framing access and participation, from the perspective of the ‘produser’ (Bruns, 2008), around three separate but related issues: the failure to frame the discussion that will be undertaken by the Australian Law Reform Commission’s 2012 2013 Copyright Inquiry; the prioritising of the market over and above media accountability and the health of the public sphere; and the missed opportunity to develop a national framework for digital literacy and advanced digital citizenry.

Convergence in the gender wage gap in Australia over the 1980s : identifying the role of counteracting forces via the Juhn, Murphy and Pierce decomposition

The paper utilises the Juhn Murphy and Pierce (1991) decomposition to shed light on the pattern of slow male-female wage convergance in Australia over the 1980s. The analysis allows one to distinguish between the role of wage structure and genderspecific effects. The central question addressed is whether rising wage inequality counteracted the forces of increased female investment in labour market skills, i.e. education and experience. The conclusion is that in contrast to the US and the UK, Australian women do not appear to have been swimming against a tide of adverse wage structure changes.

Convergence and divergence during the adaptation to similar environments by an Australian Groundsel

Adaptation to replicate environments is often achieved through similar phenotypic solutions. Whether selection also produces convergent genomic changes in these situations remains largely unknown. The variable groundsel, Senecio lautus, is an excellent system to investigate the genetic underpinnings of convergent evolution, because morphologically similar forms of these plants have adapted to the same environments along the coast of Australia. We compared range-wide patterns of genomic divergence in natural populations of this plant and searched for regions putatively affected by natural selection. Our results indicate that environmental adaptation followed complex genetic trajectories, affecting multiple loci, implying both the parallel recruitment of the same alleles and the divergence of completely different genomic regions across geography. An analysis of the biological functions of candidate genes suggests that adaptation to coastal environments may have occurred through the recruitment of different genes participating in similar processes. The relatively low genetic convergence that characterizes the parallel evolution of S. lautus forms suggests that evolution is more constrained at higher levels of biological organization.

The gender wage gap in Australia - The path of future convergence

The paper attempts to project the future trend of the gender wage gap in Australia up to 2031. The empirical analysis utilises the Income Distribution Survey (1996) together with Australian Bureau of Statistics (ABS) demographic projections. The methodology combines the ABS projections with assumptions relating to the evolution of educational attainment in order to project the future distribution of human capital skills and consequently the future size of the gender wage gap. The analysis suggests that female relative pay will continue to rise up to 2031. However, gender wage convergence will be relatively slow, with a substantial gap remaining in 2031.

Distinguishing between mean-field, moment dynamics and stochastic descriptions of birth-death-movement processes

Mathematical descriptions of birth–death–movement processes are often calibrated to measurements from cell biology experiments to quantify tissue growth rates. Here we describe and analyze a discrete model of a birth–death-movement process applied to a typical two–dimensional cell biology experiment. We present three different descriptions of the system: (i) a standard mean–field description which neglects correlation effects and clustering; (ii) a moment dynamics description which approximately incorporates correlation and clustering effects, and; (iii) averaged data from repeated discrete simulations which directly incorporates correlation and clustering effects. Comparing these three descriptions indicates that the mean–field and moment dynamics approaches are valid only for certain parameter regimes, and that both these descriptions fail to make accurate predictions of the system for sufficiently fast birth and death rates where the effects of spatial correlations and clustering are sufficiently strong. Without any method to distinguish between the parameter regimes where these three descriptions are valid, it is possible that either the mean–field or moment dynamics model could be calibrated to experimental data under inappropriate conditions, leading to errors in parameter estimation. In this work we demonstrate that a simple measurement of agent clustering and correlation, based on coordination number data, provides an indirect measure of agent correlation and clustering effects, and can therefore be used to make a distinction between the validity of the different descriptions of the birth–death–movement process.

Heat-producing crust regulation of subsurface temperatures : a stochastic model re-evaluation of the geothermal potential in western Queensland, Australia

A large subsurface, elevated temperature anomaly is well documented in Central Australia. High Heat Producing Granites (HHPGs) intersected by drilling at Innamincka are often assumed to be the dominant cause of the elevated subsurface temperatures, although their presence in other parts of the temperature anomaly has not been confirmed. Geological controls on the temperature anomaly remain poorly understood. Additionally, methods previously used to predict temperature at 5 km depth in this area are simplistic and possibly do not give an accurate representation of the true distribution and magnitude of the temperature anomaly. Here we re-evaluate the geological controls on geothermal potential in the Queensland part of the temperature anomaly using a stochastic thermal model. The results illustrate that the temperature distribution is most sensitive to the thermal conductivity structure of the top 5 km. Furthermore, the results indicate the presence of silicic crust enriched in heat producing elements between and 40 km.

Structure-preserving Runge-Kutta methods for stochastic Hamiltonian equations with additive noise

There has been considerable recent work on the development of energy conserving one-step methods that are not symplectic. Here we extend these ideas to stochastic Hamiltonian problems with additive noise and show that there are classes of Runge-Kutta methods that are very effective in preserving the expectation of the Hamiltonian, but care has to be taken in how the Wiener increments are sampled at each timestep. Some numerical simulations illustrate the performance of these methods.

A stochastic thermostat algorithm for coarse-grained thermomechanical modeling of large-scale soft matters : theory and application to microfilaments

As all-atom molecular dynamics method is limited by its enormous computational cost, various coarse-grained strategies have been developed to extend the length scale of soft matters in the modeling of mechanical behaviors. However, the classical thermostat algorithm in highly coarse-grained molecular dynamics method would underestimate the thermodynamic behaviors of soft matters (e.g. microfilaments in cells), which can weaken the ability of materials to overcome local energy traps in granular modeling. Based on all-atom molecular dynamics modeling of microfilament fragments (G-actin clusters), a new stochastic thermostat algorithm is developed to retain the representation of thermodynamic properties of microfilaments at extra coarse-grained level. The accuracy of this stochastic thermostat algorithm is validated by all-atom MD simulation. This new stochastic thermostat algorithm provides an efficient way to investigate the thermomechanical properties of large-scale soft matters.