Efficient time-independent wave packet scattering calculations within a Lanczos subspace: H+O-2 (J=0) state-to-state reaction probabilities

An efficient Lanczos subspace method has been devised for calculating state-to-state reaction probabilities. The method recasts the time-independent wave packet Lippmann-Schwinger equation [Kouri , Chem. Phys. Lett. 203, 166 (1993)] inside a tridiagonal (Lanczos) representation in which action of the causal Green's operator is affected easily with a QR algorithm. The method is designed to yield all state-to-state reaction probabilities from a given reactant-channel wave packet using a single Lanczos subspace; the spectral properties of the tridiagonal Hamiltonian allow calculations to be undertaken at arbitrary energies within the spectral range of the initial wave packet. The method is applied to a H+O-2 system (J=0), and the results indicate the approach is accurate and stable. (C) 2002 American Institute of Physics.

State-to-state reaction probabilities within the quantum transition state framework

Rigorous quantum dynamics calculations of reaction rates and initial state-selected reaction probabilities of polyatomic reactions can be efficiently performed within the quantum transition state concept employing flux correlation functions and wave packet propagation utilizing the multi-configurational time-dependent Hartree approach. Here, analytical formulas and a numerical scheme extending this approach to the calculation of state-to-state reaction probabilities are presented. The formulas derived facilitate the use of three different dividing surfaces: two dividing surfaces located in the product and reactant asymptotic region facilitate full state resolution while a third dividing surface placed in the transition state region can be used to define an additional flux operator. The eigenstates of the corresponding thermal flux operator then correspond to vibrational states of the activated complex. Transforming these states to reactant and product coordinates and propagating them into the respective asymptotic region, the full scattering matrix can be obtained. To illustrate the new approach, test calculations study the D + H2(ν, j) → HD(ν′, j′) + H reaction for J = 0.

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.

A real symmetric Lanczos subspace implementation of quantum scattering using boundary inhomogeneities

We present an efficient and robust method for calculating state-to-state reaction probabilities utilising the Lanczos algorithm for a real symmetric Hamiltonian. The method recasts the time-independent Artificial Boundary Inhomogeneity technique recently introduced by Jang and Light (J. Chem. Phys. 102 (1995) 3262) into a tridiagonal (Lanczos) representation. The calculation proceeds at the cost of a single Lanczos propagation for each boundary inhomogeneity function and yields all state-to-state probabilities (elastic, inelastic and reactive) over an arbitrary energy range. The method is applied to the collinear H + H-2 reaction and the results demonstrate it is accurate and efficient in comparison with previous calculations. (C) 2002 Elsevier Science B.V. All rights reserved.

Lanczos subspace time-independent wave packet calculations of S(D-1)+H-2 reactive scattering

In this paper. we present the results of quantum dynamical simulations of the S (D-1) + H-2 insertion reaction on a newly developed potential energy surface (J. Chem. Phys. 2001, 114, 320). State-to-state reaction probabilities. product state distributions, and initial-state resolved cumulative reaction probabilities from a given incoming reactant channel are obtained from a time-independent wave packet analysis, performed within a single Lanczos subspace. Integral reaction cross sections are then estimated by J-shifting method and compared with the results from molecular beam experiment and QCT calculations.

Chebyshev real wave packet propagation: H+O-2 (J=0) state-to-state reactive scattering calculations

In this paper we explore the relative performance of two recently developed wave packet methodologies for reactive scattering, namely the real wave packet Chebyshev domain propagation of Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)] and the Lanczos subspace wave packet approach of Smith [J. Chem. Phys. 116, 2354 (2002); Chem. Phys. Lett. 336, 149 (2001)]. In the former method, a modified Schrodinger equation is employed to propagate the real part of the wave packet via the well-known Chebyshev iteration. While the time-dependent wave packet from the modified Schrodinger equation is different from that obtained using the standard Schrodinger equation, time-to-energy Fourier transformation yields wave functions which differ only trivially by normalization. In the Lanczos subspace approach the linear system of equations defining the action of the Green operator may be solved via either time-dependent or time-independent methods, both of which are extremely efficient due to the simple tridiagonal structure of the Hamiltonian in the Lanczos representation. The two different wave packet methods are applied to three dimensional reactive scattering of H+O-2 (total J=0). State-to-state reaction probabilities, product state distributions, as well as initial-state-resolved cumulative reaction probabilities are examined. (C) 2002 American Institute of Physics.

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. Future directions in this field are also discussed.

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.

Reaction-diffusion lattice gas: Theory and computer results

We report on the study of nonequilibrium ordering in the reaction-diffusion lattice gas. It is a kinetic model that relaxes towards steady states under the simultaneous competition of a thermally activated creation-annihilation $(reaction$) process at temperature T, and a diffusion process driven by a heat bath at temperature T?T. The phase diagram as one varies T and T, the system dimension d, the relative priori probabilities for the two processes, and their dynamical rates is investigated. We compare mean-field theory, new Monte Carlo data, and known exact results for some limiting cases. In particular, no evidence of Landau critical behavior is found numerically when d=2 for Metropolis rates but Onsager critical points and a variety of first-order phase transitions.

Detecting migrants in populations of Rhizoctonia solani anastomosis group 3 from potato in North Carolina using multilocus genotype probabilities

The relative contribution of migration of Rhizoctonia solani anastomosis group 3 (AG-3) on infested potato seed tubers originating from production areas in Canada, Maine, and Wisconsin (source population) to the genetic diversity and structure of populations of R. solani AG-3 in North Carolina (NC) soil (recipient population) was examined. The frequency of alleles detected by multilocus polymerase chain reaction-restriction fragment length polymorphisms, heterozygosity at individual loci, and gametic phase disequilibrium between all pairs of loci were determined for subpopulations of R. solani AG-3 from eight sources of potato seed tubers and from five soils in NC. Analysis of molecular variation revealed little variation between seed source and NC recipient soil populations or between subpopulations within each region. Analysis of population data with a Bayesian-based statistical method previously developed for detecting migration in human populations suggested that six multilocus genotypes from the NC soil population had a statistically significant probability of being migrants from the northern source population. The one-way (unidirectional) migration of genotypes of R. solani AG-3 into NC on infested potato seed tubers from Canada, Maine, and Wisconsin provides a plausible explanation for the lack of genetic subdivision (differentiation) between populations of the pathogen in NC soils or between the northern source and the NC recipient soil populations.

Influence Of Substrate Steps On The Catalytic Properties Of Pt Layers: The Ethanol Electrooxidation Reaction.

The ethanol oxidation reaction (EOR) is investigated on Pt/Au(hkl) electrodes. The Au(hkl) single crystals used belong to the [n(111)x(110)] family of planes. Pt is deposited following the galvanic exchange of a previously deposited Cu monolayer using a Pt(2+) solution. Deposition is not epitaxial and the defects on the underlying Au(hkl) substrates are partially transferred to the Pt films. Moreover, an additional (100)-step-like defect is formed, probably as a result of the strain resulting from the Pt and Au lattice mismatch. Regarding the EOR, both vicinal Pt/Au(hkl) surfaces exhibit a behavior that differs from that expected for stepped Pt; for instance, the smaller the step density on the underlying Au substrate, the greater the ability to break the CC bond in the ethanol molecule, as determined by in situ Fourier transform infrared spectroscopy measurements. Also, we found that the acetic acid production is favored as the terrace width decreases, thus reflecting the inefficiency of the surface array to cleave the ethanol molecule.

Compact Ag@fe3o4 Core-shell Nanoparticles By Means Of Single-step Thermal Decomposition Reaction.

A temperature pause introduced in a simple single-step thermal decomposition of iron, with the presence of silver seeds formed in the same reaction mixture, gives rise to novel compact heterostructures: brick-like Ag@Fe3O4 core-shell nanoparticles. This novel method is relatively easy to implement, and could contribute to overcome the challenge of obtaining a multifunctional heteroparticle in which a noble metal is surrounded by magnetite. Structural analyses of the samples show 4 nm silver nanoparticles wrapped within compact cubic external structures of Fe oxide, with curious rectangular shape. The magnetic properties indicate a near superparamagnetic like behavior with a weak hysteresis at room temperature. The value of the anisotropy involved makes these particles candidates to potential applications in nanomedicine.

Polyphenol-enriched Cocoa Protects The Diabetic Retina From Glial Reaction Through The Sirtuin Pathway.

Cocoa is rich in flavonoids, which are potent antioxidants with established benefits for cardiovascular health but unproven effects on neurodegeneration. Sirtuins (SIRTs), which make up a family of deacetylases, are thought to be sensitive to oxidation. In this study, the possible protective effects of cocoa in the diabetic retina were assessed. Rat Müller cells (rMCs) exposed to normal or high glucose (HG) or H2O2 were submitted to cocoa treatment in the presence or absence of SIRT-1 inhibitor and small interfering RNA The experimental animal study was conducted in streptozotocin-induced diabetic rats randomized to receive low-, intermediate-, or high-polyphenol cocoa treatments via daily gavage for 16 weeks (i.e., 0.12, 2.9 or 22.9 mg/kg/day of polyphenols). The rMCs exposed to HG or H2O2 exhibited increased glial fibrillary acidic protein (GFAP) and acetyl-RelA/p65 and decreased SIRT1 activity/expression. These effects were cancelled out by cocoa, which decreased reactive oxygen species production and PARP-1 activity, augmented the intracellular pool of NAD(+), and improved SIRT1 activity. The rat diabetic retinas displayed the early markers of retinopathy accompanied by markedly impaired electroretinogram. The presence of diabetes activated PARP-1 and lowered NAD(+) levels, resulting in SIRT1 impairment. This augmented acetyl RelA/p65 had the effect of up-regulated GFAP. Oral administration of polyphenol cocoa restored the above alterations in a dose-dependent manner. This study reveals that cocoa enriched with polyphenol improves the retinal SIRT-1 pathway, thereby protecting the retina from diabetic milieu insult.

FiltekTM Silorane and FiltekTM Supreme XT resins: tissue reaction after subcutaneous implantation in isogenic mice

The aim of this study was to evaluate the tissue compatibility of a silorane-based resin system (FiltekTM Silorane) and a methacrylate-based nanoparticle resin (FiltekTM Supreme XT) after implantation in the subcutaneous connective tissue of isogenic mice. One hundred and thirty five male isogenic BALB/c mice were randomly assigned to 12 experimental and 3 control groups, according to the implanted material and the experimental period of 7, 21 and 63 days. At the end of each period, the animals were killed and the tubes with the surrounding tissues were removed and processed for microscopic analysis. Samples were subjected to a descriptive and a semi-quantitative analyses using a 4-point scoring system (0-3) to evaluate the collagen fiber formation and inflammatory infiltrate. Data were statistically analyzed using the Kruskal Wallis test (?=0.05). The results showed that there was no significant difference between the experimental and control groups considering the three evaluation periods (p>0.05). The silorane-based and the methacrylate-based nanoparticle resins presented similar tissue response to that of the empty tube (control group) after subcutaneous implantation in isogenic mice.

Effect of phosphate-bonded investments on titanium reaction layer and crown fit

This study analyzed the reaction layer and measured the marginal crown fit of cast titanium applied to different phosphate-bonded investments, prepared under the following conditions (liquid concentration/casting temperature): Rema Exakt (RE) - 100%/237°C, 75%/287°C, Castorit Super C (CS)-100%/70°C, 75%/141°C and Rematitan Plus (RP)- 100%/430°C (special to titanium cast, as the control group). The reaction layer was studied using the Vickers hardness test, and analyzed by two way ANOVA and Tukey's HSD tests (α = 0.05). Digital photographs were taken of the crowns seated on the die, the misfit was measured using an image analysis system and One-way ANOVA, and Tukey's test was applied (α = 0.05). The hardness decreased from the surface (601.17 VHN) to 150 μm (204.03 VHN). The group CS 75%/141°C presented higher hardness than the other groups, revealing higher surface contamination, but there were no differences among the groups at measurements deeper than 150 μm. The castings made with CS - 100%/70°C presented the lowest levels of marginal misfit, followed by RE -100%/237°C. The conventional investments CS (100%) and RE (100%) showed better marginal fit than RP, but the CS (75%) had higher surface contamination.