ePrints@IIScePrints@IISc Home | About | Browse | Latest Additions | Advanced Search | Contact | Help

Physically consistent simulation of mesoscale chemical kinetics: The non-negative FIS-alpha method

Dana, Saswati and Raha, Soumyendu (2011) Physically consistent simulation of mesoscale chemical kinetics: The non-negative FIS-alpha method. In: Journal of Computational Physics, 230 (24). pp. 8813-8834.

[img] PDF
Phusically.pdf - Published Version
Restricted to Registered users only

Download (3MB) | Request a copy
Official URL: http://dx.doi.org/10.1016/j.jcp.2011.07.032

Abstract

Biochemical pathways involving chemical kinetics in medium concentrations (i.e., at mesoscale) of the reacting molecules can be approximated as chemical Langevin equations (CLE) systems. We address the physically consistent non-negative simulation of the CLE sample paths as well as the issue of non-Lipschitz diffusion coefficients when a species approaches depletion and any stiffness due to faster reactions. The non-negative Fully Implicit Stochastic alpha (FIS alpha) method in which stopped reaction channels due to depleted reactants are deleted until a reactant concentration rises again, for non-negativity preservation and in which a positive definite Jacobian is maintained to deal with possible stiffness, is proposed and analysed. The method is illustrated with the computation of active Protein Kinase C response in the Protein Kinase C pathway. (C) 2011 Elsevier Inc. All rights reserved.

Item Type: Journal Article
Additional Information: Copyright of this article belongs to Elsevier Science.
Keywords: Chemical Langevin equations;Stiff stochastic differential equations;Meso-scale kinetics;Implicit method of numerical integration
Department/Centre: Division of Interdisciplinary Research > Supercomputer Education & Research Centre
Depositing User: Id for Latest eprints
Date Deposited: 21 Dec 2011 09:07
Last Modified: 21 Dec 2011 09:07
URI: http://eprints.iisc.ac.in/id/eprint/42631

Actions (login required)

View Item View Item