Simulation methods with extended stability for stiff biochemical kinetics

Background: With increasing computer power, simulating the dynamics of complex systems in chemistry and biology is becoming increasingly routine. The modelling of individual reactions in (bio)chemical systems involves a large number of random events that can be simulated by the stochastic simulation algorithm (SSA). The key quantity is the step size, or waiting time, τ, whose value inversely depends on the size of the propensities of the different channel reactions and which needs to be re-evaluated after every firing event. Such a discrete event simulation may be extremely expensive, in particular for stiff systems where τ can be very short due to the fast kinetics of some of the channel reactions. Several alternative methods have been put forward to increase the integration step size. The so-called τ-leap approach takes a larger step size by allowing all the reactions to fire, from a Poisson or Binomial distribution, within that step. Although the expected value for the different species in the reactive system is maintained with respect to more precise methods, the variance at steady state can suffer from large errors as τ grows. Results: In this paper we extend Poisson τ-leap methods to a general class of Runge-Kutta (RK) τ-leap methods. We show that with the proper selection of the coefficients, the variance of the extended τ-leap can be well-behaved, leading to significantly larger step sizes.Conclusions: The benefit of adapting the extended method to the use of RK frameworks is clear in terms of speed of calculation, as the number of evaluations of the Poisson distribution is still one set per time step, as in the original τ-leap method. The approach paves the way to explore new multiscale methods to simulate (bio)chemical systems.

Chebyshev-type quadrature formulas for new weight classes

We give Chebyshev-type quadrature formulas for certain new weight classes. These formulas are of highest possible degree when the number of nodes is a power of 2. We also describe the nodes in a constructive way, which is important for applications. One of our motivations to consider these type of problems is the Faraday cage phenomenon for discrete charges as discussed by J. Korevaar and his colleagues.

Additional compact formulas for vibrational dynamic dipole polarizabilities and hyperpolarizabilities

Compact expressions, complete through second order in electrical and/or mechanical anharmonicity, are given for the dynamic dipole vibrational polarizability and dynamic first and second vibrational hyperpolarizabilities. Certain contributions not previously formulated are now included

Macaulay style formulas for sparse resultants

We present formulas for computing the resultant of sparse polyno- mials as a quotient of two determinants, the denominator being a minor of the numerator. These formulas extend the original formulation given by Macaulay for homogeneous polynomials.

Restricted walks in regular trees

Let T be the Cayley graph of a finitely generated free group F. Given two vertices in T consider all the walks of a given length between these vertices that at a certain time must follow a number of predetermined steps. We give formulas for the number of such walks by expressing the problem in terms of equations in F and solving the corresponding equations.

Elements of algebraic geometry and the positive theory of partially commutative groups

The first main result of the paper is a criterion for a partially commutative group G to be a domain. It allows us to reduce the study of algebraic sets over G to the study of irreducible algebraic sets, and reduce the elementary theory of G (of a coordinate group over G) to the elementary theories of the direct factors of G (to the elementary theory of coordinate groups of irreducible algebraic sets). Then we establish normal forms for quantifier-free formulas over a non-abelian directly indecomposable partially commutative group H. Analogously to the case of free groups, we introduce the notion of a generalised equation and prove that the positive theory of H has quantifier elimination and that arbitrary first-order formulas lift from H to H * F, where F is a free group of finite rank. As a consequence, the positive theory of an arbitrary partially commutative group is decidable.

Weighted limits in simplicial homotopy theory

We extend the theory of Quillen adjunctions by combining ideas of homotopical algebra and of enriched category theory. Our results describe how the formulas for homotopy colimits of Bousfield and Kan arise from general formulas describing the derived functor of the weighted colimit functor.

Horospheres in hyperbolic geometry

In this paper we investigate the role of horospheres in Integral Geometry and Differential Geometry. In particular we study envelopes of families of horocycles by means of “support maps”. We define invariant “linear combinations” of support maps or curves. Finally we obtain Gauss-Bonnet type formulas and Chern-Lashof type inequalities.

Fum i miralls: trobant l'ordre en un món caòtic

El calidoscopi és una joguina que converteix fragments de vidre acolorit en patrons infinitament variats mitjançant l'ús de miralls. Hi ha quelcom satisfactori en els patrons –els necessitem per crear significat, i trobem bellesa en la seva simetria. No obstant, els sistemes complexos desobeeixen les fórmules simples. Per exemple, és impossible descriure una cosa tan prosaica com el comportament del fum mentre ascendeix fent servir les matemàtiques clàssiques. Agafant en préstec idees de diferents disciplines, podem considerar l'alfabetització informacional com una invitació per als nostres estudiants per jugar amb la imprevisible varietat de l'experiència humana i descobrir patrons que els resultin plens de significat.

The Peace Processes in the Basque Country and Northern Ireland (1994-2006): a Comparative Approach

The peace process in Northern Ireland demonstrates that new sovereignty formulas need to be explored in order to meet the demands of the populations and territories in conflict. The profound transformation of the classic symbolic elements of the nation-state within the context of the European Union has greatly contributed to the prospects for a resolution of this old conflict. Today’s discussions are focused on the search for instruments of shared sovereignty that are adapted to a complex and plural social reality. This new approach for finding a solution to the Irish conflict is particularly relevant to the Basque debate about formulating creative and modern solutions to similar conflicts over identity and sovereignty. The notion of shared sovereignty implemented in Northern Ireland –a formula for complex interdependent relations– is of significant relevance to the broader international community and is likely to become an increasingly potent and transcendent model for conflict resolution and peace building.

Anàlisi de la qualitat de l'aigua: estudi dels impactes humans sobre la qualitat de l'aigua en pous i basses del Parc del Garraf i rodalies

Amb la finalitat de conèixer l’estat de qualitat de les aigües de les basses i pous del Parc del Garraf, s’analitzen una sèrie de paràmetres fisico-químics en 17 estacions de mostreig prèviament seleccionades, distribuïdes en zones amb diferents tipologies d’ús del sòl. La base de l’anàlisi ha estat la integració d’informació provinent de diferents fonts. Mitjançant l’elaboració de taules i gràfics, la generació de cartografia i el tractament estadístic de les dades, s’ha procedit a la tria de punts de mostreig i s’ha obtingut un inventari que ha permès la interpretació global dels resultats, facilitant la diagnosi. El procés de tractament de dades inclou la confecció d’un índex de qualitat de les aigües (ICA) propi, no vinculant, a partir de fórmules genèriques de normalització i ponderació de valors. Durant la realització de la diagnosi s’han detectat pertorbacions puntuals en determinats paràmetres corresponents a contaminacions locals, en diferents estacions de mostreig. Aquestes pertorbacions s’han relacionat amb la situació dels pous i les basses al Garraf i les tipologies d’ús del sòl de cada zona. El diagnòstic de pertorbacions ha orientat les propostes de millora aplicables que s’han dividit en tres classes segons el nivell d’actuació. Aquestes incideixen principalment, en la millora de la informació disponible, l’aplicació de l’agricultura ecològica, l’explotació sostenible dels aqüífers i la realització d’estudis globals i/o locals, més complets i exhaustius.

Sharpened lower bounds for cut elimination

We present sharpened lower bounds on the size of cut free proofs for first-order logic. Prior lower bounds for eliminating cuts from a proof established superexponential lower bounds as a stack of exponentials, with the height of the stack proportional to the maximum depth d of the formulas in the original proof. Our new lower bounds remove the constant of proportionality, giving an exponential stack of height equal to d − O(1). The proof method is based on more efficiently expressing the Gentzen-Solovay cut formulas as low depth formulas.

Haar wavelets-based approach for quantifying credit portfolio losses

This paper proposes a new methodology to compute Value at Risk (VaR) for quantifying losses in credit portfolios. We approximate the cumulative distribution of the loss function by a finite combination of Haar wavelet basis functions and calculate the coefficients of the approximation by inverting its Laplace transform. The Wavelet Approximation (WA) method is specially suitable for non-smooth distributions, often arising in small or concentrated portfolios, when the hypothesis of the Basel II formulas are violated. To test the methodology we consider the Vasicek one-factor portfolio credit loss model as our model framework. WA is an accurate, robust and fast method, allowing to estimate VaR much more quickly than with a Monte Carlo (MC) method at the same level of accuracy and reliability.

The β-Meixner model

We propose to approximate the Meixner model by a member of the B–family introduced in [Kuz10a]. The advantage of such approximations are the semi–explicit formulas for the running extrema under the B–family processes which enables us to produce more efficient algorithms for certain path dependent options.

The β-Meixner model

We propose to approximate the Meixner model by a member of the B-family introduced in [Kuz10a]. The advantage of such approximations are the semi-explicit formulas for the running extrema under the B-family processes which enables us to produce more efficient algorithms for certain path dependent options.