Premium
Solving the unimolecular master equation with a weighted subspace projection method
Author(s) -
Frankcombe Terry J.,
Smith Sean C.
Publication year - 2000
Publication title -
journal of computational chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.907
H-Index - 188
eISSN - 1096-987X
pISSN - 0192-8651
DOI - 10.1002/(sici)1096-987x(200006)21:8<592::aid-jcc2>3.0.co;2-2
Subject(s) - eigenvalues and eigenvectors , subspace topology , weighting , transpose , mathematics , projection (relational algebra) , matrix (chemical analysis) , equivalence (formal languages) , master equation , dissociation (chemistry) , mathematical analysis , algorithm , chemistry , pure mathematics , quantum mechanics , physics , quantum , chromatography , acoustics
A weighted subspace projection method for solving the unimolecular master equation over a wide range of temperatures and pressures is developed. Sample calculations modeling the dissociation of ethane at 300 K and pressures as low as 0.65 Torr demonstrates the utility of the method in regimes where standard projection methods fail. For the sample calculations the weighted Arnoldi method was able to reliably calculate the smallest eigenvalue of the rate matrix in excellent agreement with calculations using the Nesbet algorithm. Extremely small eigenvalues of the order of −10 −48 could be calculated without difficulty. The formal equivalence between various weighting schemes and common matrix transformations is shown. The point that merely taking the transpose of the rate matrix can be extremely beneficial is made, commenting on the relationship between the left and right eigenvectors of the rate matrix. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 592–606, 2000