Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Selecting methods to solve multi-well master equations

There are several competing methods commonly used to solve energy grained master equations describing gas-phase reactive systems. When it comes to selecting an appropriate method for any particular problem, there is little guidance in the literature. In this paper we directly compare several variants of spectral and numerical integration methods from the point of view of computer time required to calculate the solution and the range of temperature and pressure conditions under which the methods are successful. The test case used in the comparison is an important reaction in combustion chemistry and incorporates reversible and irreversible bimolecular reaction steps as well as isomerizations between multiple unimolecular species. While the numerical integration of the ODE with a stiff ODE integrator is not the fastest method overall, it is the fastest method applicable to all conditions.

A Lanczos-powered implementation of the Faber polynomial quantum time propagator for reaction probabilities

Complementing our recent work on subspace wavepacket propagation [Chem. Phys. Lett. 336 (2001) 149], we introduce a Lanczos-based implementation of the Faber polynomial quantum long-time propagator. The original version [J. Chem. Phys. 101 (1994) 10493] implicitly handles non-Hermitian Hamiltonians, that is, those perturbed by imaginary absorbing potentials to handle unwanted reflection effects. However, like many wavepacket propagation schemes, it encounters a bottleneck associated with dense matrix-vector multiplications. Our implementation seeks to reduce the quantity of such costly operations without sacrificing numerical accuracy. For some benchmark scattering problems, our approach compares favourably with the original. (C) 2004 Elsevier B.V. All rights reserved.

Symmetry contaminations in reactive scattering through long-lived collision complexes

We extend our Lanczos subspace time-independent wave packet method [J. Chem. Phys. 116 (2002) 2354] to investigate the issue of symmetry contaminations for the challenging deep-well H + O-2 reaction. Our central objective is to address the issue of whether significant symmetry contamination can occur if a wavepacket initially possessing the correct O-O exchange symmetry is propagated over tens of thousands of recursive steps using a basis which does not explicitly enforce the correct symmetry, and if so how seriously this affects the results. We find that symmetry contamination does exist where the symmetry constraint is not explicitly enforced in the basis. While it affects individual resonances and the associated peak amplitudes, the overall shape of the more averaged quantities such as total reaction probabilities and vibrational branching ratios are not seriously affected. (C) 2004 Elsevier B.V. All rights reserved.

On the microscopic mechanism of carbon gasification: A theoretical study

In this paper we report the results of ab initio calculations on the energetics and kinetics of oxygen-driven carbon gasification reactions using a small model cluster, with full characterisation of the stationary points on the reaction paths. We show that previously unconsidered pathways present significantly reduced barriers to reaction and must be considered as alternative viable paths. At least two electronic spin states of the model cluster must be considered for a complete description. (C) 2004 Elsevier Ltd. All rights reserved.

Photodissociation of benzene under collision-free conditions: An ab initio/Rice-Ramsperger-Kassel-Marcus study

The ab initio/Rice-Ramsperger-Kassel-Marcus (RRKM) approach has been applied to investigate the photodissociation mechanism of benzene at various wavelengths upon absorption of one or two UV photons followed by internal conversion into the ground electronic state. Reaction pathways leading to various decomposition products have been mapped out at the G2M level and then the RRKM and microcanonical variational transition state theories have been applied to compute rate constants for individual reaction steps. Relative product yields (branching ratios) for C6H5+H, C6H4+H-2, C4H4+C2H2, C4H2+C2H4, C3H3+C3H3, C5H3+CH3, and C4H3+C2H3 have been calculated subsequently using both numerical integration of kinetic master equations and the steady-state approach. The results show that upon absorption of a 248 nm photon dissociation is too slow to be observable in molecular beam experiments. In photodissociation at 193 nm, the dominant dissociation channel is H atom elimination (99.6%) and the minor reaction channel is H-2 elimination, with the branching ratio of only 0.4%. The calculated lifetime of benzene at 193 nm is about 11 mus, in excellent agreement with the experimental value of 10 mus. At 157 nm, the H loss remains the dominant channel but its branching ratio decreases to 97.5%, while that for H-2 elimination increases to 2.1%. The other channels leading to C3H3+C3H3, C5H3+CH3, C4H4+C2H2, and C4H3+C2H3 play insignificant role but might be observed. For photodissociation upon absorption of two UV photons occurring through the neutral hot benzene mechanism excluding dissociative ionization, we predict that the C6H5+H channel should be less dominant, while the contribution of C6H4+H-2 and the C3H3+C3H3, CH3+C5H3, and C4H3+C2H3 radical channels should significantly increase. (C) 2004 American Institute of Physics.

Full S matrix calculation via a single real-symmetric Lanczos recursion: The Lanczos artificial boundary inhomogeneity method

We present an efficient and robust method for the calculation of all S matrix elements (elastic, inelastic, and reactive) over an arbitrary energy range from a single real-symmetric Lanczos recursion. Our new method transforms the fundamental equations associated with Light's artificial boundary inhomogeneity approach [J. Chem. Phys. 102, 3262 (1995)] from the primary representation (original grid or basis representation of the Hamiltonian or its function) into a single tridiagonal Lanczos representation, thereby affording an iterative version of the original algorithm with greatly superior scaling properties. The method has important advantages over existing iterative quantum dynamical scattering methods: (a) the numerically intensive matrix propagation proceeds with real symmetric algebra, which is inherently more stable than its complex symmetric counterpart; (b) no complex absorbing potential or real damping operator is required, saving much of the exterior grid space which is commonly needed to support these operators and also removing the associated parameter dependence. Test calculations are presented for the collinear H+H-2 reaction, revealing excellent performance characteristics. (C) 2004 American Institute of Physics.

Recent developments in the quantum dynamical characterization of unimolecular resonances

We give a selective review of quantum mechanical methods for calculating and characterizing resonances in small molecular systems, with an emphasis on recent progress in Chebyshev and Lanczos iterative methods. Two archetypal molecular systems are discussed: isolated resonances in HCO, which exhibit regular mode and state specificity, and overlapping resonances in strongly bound HO2, which exhibit irregular and chaotic behavior. Recent progresses for non-zero total angular momentum J calculations of resonances including parallel computing models are also included and future directions in this field are discussed.

Transition-state-theory calculations for reactions of O(3P) with halogenated olefins

The potential energy surfaces for the reactions of atomic oxygen in its ground electronic state, O(P-3), with the olefins: CF2=CCl2 and CF2=CF - CF3, have been characterized using ab initio molecular orbital calculations. Geometry optimization and vibrational frequency calculations were performed for reactants, transition states and products at the MP2 and QCISD levels of theory using the 6-31G(d) basis set. This database was then used to calculate the rate constants by means of Transition-State-Theory. To obtain a better reference and to test the reliability of the activation barriers we have also carried out computations using the CCSD(T)(fc)/6-311Gdagger, MP4(SDQ)(fc)/CBSB4 and MP2(fc)/CBSB3 single point energy calculations at both of the above levels of theory, as well as with the composite CBS-RAD procedure ( P. M. Mayer, C. J. Parkinson, D. M. Smith and L. Radom, J. Chem. Phys., 1998, 108, 604) and a modi. cation of this approach, called: CBS-RAD( MP2, MP2). It was found that the kinetic parameters obtained in this work particularly with the CBS-RAD ( MP2, MP2) procedure are in reasonable agreement with the experimental values. For both reactions it is found that the channels leading to the olefin double-bond addition predominates with respect to any other reaction pathway. However, on account of the different substituents in the alkenes we have located, at all levels of theory, two transition states for each reaction. Moreover, we have found that, for the reactions studied, a correlation exists between the activation energies and the electronic structure of the transition states which can explain the influence of the substituent effect on the reactivity of the halo-olefins.

Dynamic assessment of articulation during lingual fatigue in myasthenia gravis

The present study aimed to investigate how induced lingual fatigue affected lingual strength, articulatory kinematics, and perceptual speech features in CS, a 51-year-old female with active myasthenia gravis (MG), and three age and gender matched control participants, Lingual fatigue was elicited via a series of endurance tasks using a tongue pressure bulb. Following each endurance task, the participants performed a speech task containing the phonemes /k/, /t/, and /j/ that was recorded with an electromagnetic articulograph, followed by a lingual strength assessment using a tongue pressure bulb. Participants repeated this schedule over five phases and kinematic and strength changes during each phase were compared to baseline measurements. All of CSs significant kinematic changes occurred during the final fatigue phase compared to 27.3% of the control group's kinematic changes occurring during this phase, suggesting the kinematic changes associated with fatigue were not accelerated in CS. The endurance tasks also elicited different kinematic effects for CSs anterior and posterior tongue segments. While CS exhibited mostly similar kinematic and perceptual changes to the control group, some of CS's perceptual transcriptions for /k/ and kinematic changes were not replicated, indicating that some different perceptual and physiological consequences to CS's speech were elicited by the endurance tasks.

Microbial diversity within the water column of a larval rearing system for the ornate rock lobster (Panulirus ornatus)

The ornate tropical rock lobster, Panulirus ornatus has substantial potential as an aquaculture species though disease outbreaks during the animal's extended larval lifecycle are major constraints for success. In order to effectively address such disease-related issues, an improved understanding of the composition and dynamics of the microbial communities in the larval rearing tanks is required. This study used flow cytometry and molecular microbial techniques (clone libraries and denaturing gradient gel electrophoresis (DGGE)) to quantify and characterise the microbial community of the water column in the early stages (developmental stage I-II) of a P. ornatus larval rearing system. DGGE analysis of a 5000 L larval rearing trial demonstrated a dynamic microbial community with distinct changes in the community structure after initial stocking (day I to day 2) and from day 4 to day 5, after which the structure was relatively stable. Flow cytometry analysis of water samples taken over the duration of the trial demonstrated a major increase in bacterial load leading up to and peaking on the first day of the initial larval moult (day 7), before markedly decreasing prior to when > 50% of larvae moulted (day 9). A clone library of a day 10 water sample taken following a mass larval mortality event reflected high microbial diversity confirmed by statistical analysis indices. Sequences retrieved from both clone library and DGGE analyses were dominated by gamma- and alpha-Proteobacteria affiliated organisms with additional sequences affiliated with beta- and epsilon-Proteobacteria, Bacteroidetes, Cytophagales and Chlamydiales groups. Vibrio affiliated species were commonly retrieved in the clone library, though absent from DGGE analysis.