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.

Structural evolution in a hydrothermal reaction between Nb2O5 and NaOH solution: From Nb2O5 grains to microporous Na2Nb2O6 center dot(2)/3H2O fibers and NaNbO3 cubes

Niobium pentoxide reacts actively with concentrate NaOH solution under hydrothermal conditions at as low as 120 degrees C. The reaction ruptures the corner-sharing of NbO7 decahedra and NbO6 octahedra in the reactant Nb2O5, yielding various niobates, and the structure and composition of the niobates depend on the reaction temperature and time. The morphological evolution of the solid products in the reaction at 180 degrees C is monitored via SEM: the fine Nb2O5 powder aggregates first to irregular bars, and then niobate fibers with an aspect ratio of hundreds form. The fibers are microporous molecular sieve with a monoclinic lattice, Na2Nb2O6 center dot(2)/3H2O. The fibers are a metastable intermediate of this reaction, and they completely convert to the final product NaNbO3 Cubes in the prolonged reaction of 1 h. This study demonstrates that by carefully optimizing the reaction condition, we can selectively fabricate niobate structures of high purity, including the delicate microporous fibers, through a direct reaction between concentrated NaOH solution and Nb2O5. This synthesis route is simple and suitable for the large-scale production of the fibers. The reaction first yields poorly crystallized niobates consisting of edge-sharing NbO6 octahedra, and then the microporous fibers crystallize and grow by assembling NbO6 octahedra or clusters of NbO6 octahedra and NaO6 units. Thus, the selection of the fibril or cubic product is achieved by control of reaction kinetics. Finally, niobates with different structures exhibit remarkable differences in light absorption and photoluminescence properties. Therefore, this study is of importance for developing new functional materials by the wet-chemistry process.

Electronic structure methods applied to gas-carbon reactions

A review is given on the fundamental studies of gas-carbon reactions using electronic structure methods in the last several decades. The three types of electronic structure methods including semi-empirical, ab initio and density functional theory, methods are briefly introduced first, followed by the studies on carbon reactions with hydrogen and oxygen-containing gases (non-catalysed and catalysed). The problems yet to solve and possible promising directions are discussed. (C) 2002 Elsevier Science Ltd. 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.

Converged quantum calculations of HO2 bound states and resonances for J=6 and 10

Bound and resonance states of HO2 are calculated quantum mechanically using both the Lanczos homogeneous filter diagonalization method and the real Chebyshev filter diagonalization method for nonzero total angular momentum J=6 and 10, using a parallel computing strategy. For bound states, agreement between the two methods is quite satisfactory; for resonances, while the energies are in good agreement, the widths are in general agreement. The quantum nonzero-J specific unimolecular dissociation rates for HO2 are also calculated. (C) 2004 American Institute of Physics.

Iterative quantum computations of HO2 bound states and resonances for J=4 and 5

Bound and resonance states of HO2 have been calculated by both the complex Lanczos homogeneous filter diagonalisation (LHFD) method(1,2) and the real Chebyshev filter diagonalisation method(3,4) for non-zero total angular momentum J = 4 and 5. For bound states, the agreement between the two methods is quite satisfactory; for resonances while the energies are in good agreement, the widths are only in general agreement. The relative performances of the two iterative FD methods have also been discussed in terms of efficiency as well as convergence behaviour for such a computationally challenging problem. A helicity quantum number Ohm assignment (within the helicity conserving approximation) is performed and the results indicate that Coriolis coupling becomes more important as J increases and the helicity conserving approximation is not a good one for the HO2 resonance states.

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.

Unimolecular Rovibrational Bound and Resonance States for Large Angular Momentum: J=20 Calculations for HO2

We explore the calculation of unimolecular bound states and resonances for deep-well species at large angular momentum using a Chebychev filter diagonalization scheme incorporating doubling of the autocorrelation function as presented recently by Neumaier and Mandelshtam [Phys. Rev. Lett. 86, 5031 (2001)]. The method has been employed to compute the challenging J=20 bound and resonance states for the HO2 system. The methodology has firstly been tested for J=2 in comparison with previous calculations, and then extended to J=20 using a parallel computing strategy. The quantum J-specific unimolecular dissociation rates for HO2-> H+O-2 in the energy range from 2.114 to 2.596 eV have been reported for the first time, and comparisons with the results of Troe and co-workers [J. Chem. Phys. 113, 11019 (2000) Phys. Chem. Chem. Phys. 2, 631 (2000)] from statistical adiabatic channel method/classical trajectory calculations have been made. For most of the energies, the reported statistical adiabatic channel method/classical trajectory rate constants agree well with the average of the fluctuating quantum-mechanical rates. Near the dissociation threshold, quantum rates fluctuate more severely, but their average is still in agreement with the statistical adiabatic channel method/classical trajectory results.

Observations of the reduction of FeO from slag by graphite, coke and coal char

The reduction of FeO from iron-saturated FeO-CaO-Al2O3-SiO2 slags by graphite, coke and coal char at 1 673 K has been investigated using a sessile drop technique. Metallographic analysis of samples quenched from the reaction temperature, and in situ observations of the reaction interface, reveal significant differences in the slag/carbon contact, and in the morphologies of the product iron and its composition; these differences were found to depend on the carbon type used in the reduction. In particular it has been shown that, in the case of graphite and coke, liquid Fe-C droplets were rapidly formed at the slag/C interface. Reactions of the slag with coal chars, in contrast, result predominantly in the formation of solid iron. These observations indicate that the reaction pathways, and hence reaction kinetics, are dependent on carbon type.

Discoloration and mineralization of Reactive Red HE-3B by heterogeneous photo-Fenton reaction

Discoloration and mineralization of Reactive Red HE-3B were studied by using a laponite clay-based Fe nanocomposite (Fe-Lap-RD) as a heterogeneous catalyst in the presence of H2O2 and UV light. Our experimental results clearly indicate that Fe-Lap-RD mainly consists of Fe2O3 (meghemite) and Fe2Si4O10(OH)2 (iron silicate hydroxide) which have tetragonal and monoclinic structures, respectively, and has a high specific surface area (472m(2) / g) as well as a high total pore volume (0.547 cm(3)/g). It was observed that discoloration of HE-3B undergoes a much faster kinetics than mineralization of HE-3B. It was also found that initial HE-3B concentration, H2O2 concentration, UV light wavelength and power, and Fe-Lap-RD catalyst loading are the four main factors that can significantly influence the mineralization of HE-3B. At optimal conditions, complete discoloration of 100 mg/L HE-3B can be achieved in 30 min and the total organic carbon removal ratio can attain 76% in 120 min, illustrating that Fe-Lap-RD has a high photo-catalytic activity in the photo-assisted discoloration and mineralization of HE-3B in the presence of UV light (254nm) and H2O2. (C) 2003 Elsevier Science Ltd. All rights reserved.

Degradation of azo-dye Orange II by a photoassisted Fenton reaction using a novel composite of iron oxide and silicate nanoparticles as a catalyst

A novel nanocomposite of iron oxide and silicate, prepared through a reaction between a solution of iron salt and a dispersion of Laponite clay, was used as a catalyst for the photoassisted Fenton degradation of azo-dye Orange II. This catalyst is much cheaper than the Nafion-based catalysts, and our results illustrate that it can significantly accelerate the degradation of Orange II under the irradiation of UV light (lambda = 254 nm). An advantage of the catalyst is its long-term stability that was confirmed through using the catalyst for multiple runs in the degradation of Orange II. The effects of the H2O2 molar concentration, solution pH, wavelength and power of the LTV light, catalyst loading, and initial Orange II concentration on the degradation of Orange 11 were studied in detail. In addition, it was also found that discoloration of Orange 11 undergoes a faster kinetics than mineralization of Orange II and 75% total organic carbons of 0.1 mM Orange II can be eliminated after 90 min in the presence of 1.0 g of Fe-nanocomposite/L, 4.8 mM H2O2, and 1 x 8W UVC.