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.

Calculation of bound and resonance states of HO2 for nonzero total angular momentum

Bound and resonance states of HO2 have been calculated quantum mechanically by the Lanczos homogeneous filter diagonalization method [Zhang and Smith, Phys. Chem. Chem. Phys. 3, 2282 (2001); J. Chem. Phys. 115, 5751 (2001)] for nonzero total angular momentum J = 1,2,3. For lower bound states, agreement between the results in this paper and previous work is quite satisfactory; while for high lying bound states and resonances these are the first reported results. A helicity quantum number V assignment (within the helicity conserving approximation) is performed and the results indicate that for lower bound states it is possible to assign the V quantum numbers unambiguously, but for resonances it is impossible to assign the V helicity quantum numbers due to strong mixing. In fact, for the high-lying bound states, the mixing has already appeared. These results indicate that the helicity conserving approximation is not good for the resonance state calculations and exact quantum calculations are needed to accurately describe the reaction dynamics for HO2 system. Analysis of the resonance widths shows that most of the resonances are overlapping and the interferences between them lead to large fluctuations from one resonance to another. In accord with the conclusions from earlier J = 0 calculations, this indicates that the dissociation of HO2 is essentially irregular. (C) 2003 American Institute of Physics.

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.

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.

Vascular remodelling in the internal mammary artery graft and association with elevated endothelin-1 expression

Introduction: The vasoconstricting peptide Endothelin-1 (ET-1) has been associated with atherosclerotic cardiovascular disease, AAA, hypertension and hypercholesterolemia. It is known to stimulate quiescent vascular smooth muscle cells (VSMC) into the growth cycle and has been linked to intimal thickening following endothelial injury and is associated with vessel wall remodelling in salt-sensitive hypertension models. Enhanced ET-1 expression has been reported in the internal mammary artery (IMA) and was markedly higher in patients undergoing cardiac bypass surgery who were diabetic and /or hypercholesterolemic. Aims: To firstly review the histopathology of the IMA and secondly, determine the relationship between ET-1 expression in this vessel and mitogenic activity in the medial VSMC. Methods: Vessel tissue collected at the time of CABG surgery was formalin-fixed and paraffin-embedded for histological investigation. Cross sections of the left distal IMAwere stained with Alcian Blue/Verhoeff’s van Gieson to assess medial degeneration and identify the elastic lamellae and picrosirius red to determine the collagen content (specifically type I and type III). Immunohistochemistry staining was used to assess VSMC growth (PCNA label), tissue ET-1 expression, VSMC (SMCa-actin) area and macrophage/monocyte (anti-CD68) infiltration. Quantitative analysis was performed to measure the VSMC area in relation to ET-1 staining. Results: Fifty-five IMA specimens from the CABG patients (10F; 45M; mean age 65 years) were collected for this study. Fourteen donor IMAspecimens were used as controls (7F; 7M; mean age 45 years). Significant medial hypertrophy, VSMC disorganisation and elastic lamellae destruction was detected in the CABG IMA. The amount of Alcian blue staining in the CABG IMA was almost double that of the control (31.85+/14.52% Vs 17.10+/9.96%, P= .0006). Total collagen and type I collagen content was significantly increased compared with controls (65.8+/18.3% Vs 33.7 + / 13.7%, P= .07), (14.2 + /10.0% Vs 4.8 + /2.8%, P= .01), respectively. Tissue ET-1 and PCNA labelling were also significantly elevated the CABG IMA specimens relative to the controls (69.99 + /18.74%Vs 23.33 + /20.53%, P= .0001, and 37.29 + /12.88% Vs 11.06 + /8.18, P= .0001), respectively. There was mild presence of macrophages and monocytes in both CABG and control tissue. Conclusions: The IMA from CABG patients has elevated levels of type I collagen in the extracellular matrix indicative of fibrosis and was coupled with deleterious structural remodelling. Abnormally high levels of ET-1 were measured in the medial SMC layer and was associated with VSMC growth but not related to any chronic inflammatory response within the vessel wall.

A director theory for visco-elastic folding instabilities in multilayered rock

A model for finely layered visco-elastic rock proposed by us in previous papers is revisited and generalized to include couple stresses. We begin with an outline of the governing equations for the standard continuum case and apply a computational simulation scheme suitable for problems involving very large deformations. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered beam under compression. We analyse folding up to 40% shortening. The standard continuum solution becomes unstable for extreme values of the shear/normal viscosity ratio. The instability is a consequence of the neglect of the bending stiffness/viscosity in the standard continuum model. We suggest considering these effects within the framework of a couple stress theory. Couple stress theories involve second order spatial derivatives of the velocities/displacements in the virtual work principle. To avoid C-1 continuity in the finite element formulation we introduce the spin of the cross sections of the individual layers as an independent variable and enforce equality to the spin of the unit normal vector to the layers (-the director of the layer system-) by means of a penalty method. We illustrate the convergence of the penalty method by means of numerical solutions of simple shears of an infinite layer for increasing values of the penalty parameter. For the shear problem we present solutions assuming that the internal layering is oriented orthogonal to the surfaces of the shear layer initially. For high values of the ratio of the normal-to the shear viscosity the deformation concentrates in thin bands around to the layer surfaces. The effect of couple stresses on the evolution of folds in layered structures is also investigated. (C) 2002 Elsevier Science Ltd. All rights reserved.

A methodology for current density calculations in high-frequency RF resonators

As nuclear magnetic resonance imaging and spectroscopy move inexorably toward higher field-strength magnets in search of improved signal-to-noise ratio, spectral resolution, and spatial resolution, the way in which radiofrequency (RF) probes are designed changes. At higher frequencies, resonant cavities become the favored RF ''coil'' type and may be built using streamline elements to reduce the inductance of the system. In modeling such systems, the quasi-static approach of assuming that current flows evenly in all conductor cross sections and that adjacent conductors do not affect each other becomes less reasonable. The proximity of RF conductors in resonators typically causes RF eddy currents to flow, whereby the current density in each rung is altered by the RF fields generated by nearby conductors. The proper understanding and prediction of how resonators will perform require a model of the current densities flowing in conducting sections, including all RF eddy current effects. Very few models of this type have been presented in the literature. This article presents an overview of one such model and of how it may be applied to a variety of resonators, both shielded and unshielded, circular, and elliptical, in cross section. Results are presented from a shielded head coil operating at 2 tesla. (C) 1997 John Wiley & Sons, Inc.

Proton transport model in the ionosphere .1. Multistream approach of the transport equations

The suprathermal particles, electrons and protons, coming from the magnetosphere and precipitating into the high-latitude atmosphere are an energy source of the Earth's ionosphere. They interact with ambient thermal gas through inelastic and elastic collisions. The physical quantities perturbed by these precipitations, such as the heating rate, the electron production rate, or the emission intensities, can be provided in solving the kinetic stationary Boltzmann equation. This equation yields particle fluxes as a function of altitude, energy, and pitch angle. While this equation has been solved through different ways for the electron transport and fully tested, the proton transport is more complicated. Because of charge-changing reactions, the latter is a set of two-coupled transport equations that must be solved: one for protons and the other for H atoms. We present here a new approach that solves the multistream proton/hydrogen transport equations encompassing the collision angular redistributions and the magnetic mirroring effect. In order to validate our model we discuss the energy conservation and we compare with another model under the same inputs and with rocket observations. The influence of the angular redistributions is discussed in a forthcoming paper.

Induction of smooth muscle cell nitric oxide synthase by human leukaemia inhibitory factor: Effects in vitro and in vivo

We have previously shown that human leukaemia inhibitory factor (hLIF) inhibits perivascular cuff-induced neointimal formation in the rabbit carotid artery. Since nitric oxide (NO) is a known inhibitor of smooth muscle growth, NO synthase (NOS) activity in the presence of hLIF was examined in vivo and in vitro. In rabbit aortic smooth muscle cell (SMC) culture, significant NOS activity was observed at 50 pg/ml hLIF, with maximal activity at 5 ng/ml. In the presence of the NOS inhibitor L-NAME, hLIF-induced activation of NOS was greatly decreased, however it was still 63-fold higher than in control (p < 0.05). SMC-DNA synthesis was significantly reduced (-47%) following incubation with hLIF plus L-arginine, the substrate required for NO production (p < 0.05), with no effect observed in the absence of L-arginine. Silastic cuff placement over the right carotid artery of rabbits resulted in a neointima 19.3 +/- 5.4% of total wall cross-sectional area, which was increased in the presence of L-NAME (27.0 +/- 2.0%; p < 0.05) and reduced in the presence of L-arginine (11.3 +/- 2.0%; p < 0.05). The effect of L-arginine was ameliorated by co-administration of L-NAME (16.4 +/- 1.5%). However, administration of L-NAME with hLIF had no effect on the potent inhibition of neointimal formation by hLIF (3.2 +/- 2.5 vs. 2.1 +/- 5.4%, respectively). Similarly, with hLIF administration, NOS activity in the cuffed carotid increased to 269.0 +/- 14.0% of saline-treated controls and remained significantly higher with coadministration of L-NAME (188.5 +/- 14.7%). These results indicate that hLIF causes superinduction of NO by SMC, and that it is, either partially or wholly, through this mechanism that hLIF is a potent inhibitor of neointimal formation in vivo and of smooth muscle proliferation in vitro.

Comments on alternative calculations of the broadening of spectral lines of neutral sodium by H-atom collisions

With the exception of the sodium D-lines, recent calculations of line broadening cross sections for several multiplets of sodium by Leininger et al (Leininger T, Gadea F X and Dickinson A 2000 J. Phys. B: At. Mol. Opt. Phys. 33 1805) are in substantial disagreement with cross sections interpolated from the tables of Anstee and O'Mara (Anstee and O'Mara 1995 Mon. Not. R. Astron. Soc. 276 859) and Barklem and O'Mara (Barklem P S and O'Mara B J 1997 Mon. Not. R. Astron. Soc. 290 102). The discrepancy is as large as a factor of 3 for the 3p-4d multiplet. The two theories are tested by using the results of each to synthesize lines in the solar spectrum. It is found that generally the data from the theory of Anstee, Barklem and O'Mara produce the best match to the observed solar spectrum. It is found, using a simple model for reflection of the optical electron by the potential barrier between the two atoms, that the reflection coefficient is too large for avoided crossings with the upper states of subordinate lines to contribute to line broadening, supporting the neglect of avoided ionic crossings by Anstee, Barklem and O'Mara for these lines. The large discrepancies between the two sets of calculations is a result of an approximate treatment of avoided ionic crossings for these lines by Leininger et al (Leininger T, Gadea F X and Dickinson A 2000 J. Phys. B: At. Mol. Opt. Phys. 33 1805).

Catheter effects in organ perfusion experiments

In a typical isolated organ perfusion experiment, a substance is injected upstream of an organ and then collected at some distance downstream. To reach the organ from the injection site, and then from the organ to the collector, a solute passes through catheters, usually tubes with circular cross-sections. Catheters cause distortion to the concentration-time profile of the perfusion. In this paper, we analyse catheter distribution kinetics from a mathematical point of view, develop the function most suitable for modeling this distribution and successfully apply this function to experimental data. (C) 2002 Academic Press.

Binomial leap methods for simulating stochastic chemical kinetics

This paper discusses efficient simulation methods for stochastic chemical kinetics. Based on the tau-leap and midpoint tau-leap methods of Gillespie [D. T. Gillespie, J. Chem. Phys. 115, 1716 (2001)], binomial random variables are used in these leap methods rather than Poisson random variables. The motivation for this approach is to improve the efficiency of the Poisson leap methods by using larger stepsizes. Unlike Poisson random variables whose range of sample values is from zero to infinity, binomial random variables have a finite range of sample values. This probabilistic property has been used to restrict possible reaction numbers and to avoid negative molecular numbers in stochastic simulations when larger stepsize is used. In this approach a binomial random variable is defined for a single reaction channel in order to keep the reaction number of this channel below the numbers of molecules that undergo this reaction channel. A sampling technique is also designed for the total reaction number of a reactant species that undergoes two or more reaction channels. Samples for the total reaction number are not greater than the molecular number of this species. In addition, probability properties of the binomial random variables provide stepsize conditions for restricting reaction numbers in a chosen time interval. These stepsize conditions are important properties of robust leap control strategies. Numerical results indicate that the proposed binomial leap methods can be applied to a wide range of chemical reaction systems with very good accuracy and significant improvement on efficiency over existing approaches. (C) 2004 American Institute of Physics.

Application of X-ray microtomography to numerical simulations of agglomerate breakage by distinct element method

Computer-aided tomography has been used for many years to provide significant information about the internal properties of an object, particularly in the medical fraternity. By reconstructing one-dimensional (ID) X-ray images, 2D cross-sections and 3D renders can provide a wealth of information about an object's internal structure. An extension of the methodology is reported here to enable the characterization of a model agglomerate structure. It is demonstrated that methods based on X-ray microtomography offer considerable potential in the validation and utilization of distinct element method simulations also examined.

Migration of crystals during the filling of semi-solid castings

In cold-chamber high-pressure die castings (HPDC), the microstructure consists of coarse externally solidified crystals (ESCs) that are commonly observed in the central region of cross sections. In the present work, controlled laboratory scale casting experiments have been conducted with particular emphasis on the flow and solidification conditions. An A356 aluminum alloy was used to produce castings by pouring semi-solid metal through a steel die. Microstructures similar to those encountered in HPDC have been produced and the resulting microstructure is found to depend on the melt and die temperature: (1) the fraction of ESCs determines the extent of migration to the central region; (2) a maximum packing determines the area fraction of ESCs in the center; and (3) the die temperature affects the position of the ESCs-a higher die temperature can induce a displaced ESC distribution. It is found that the n-figration of crystals to the central region requires a flow, which is constrained at all melt/die interfaces. Furthermore, potential lift mechanisms are discussed. An assessment of the Saffman lift force on individual particles shows it has no significant effect on the migration of ESCs.

Microhardness mapping and the hardness-yield strength relationship in high-pressure diecast magnesium alloy AZ91

Microhardness maps of cross-sections of high-pressure diecast test bars of AZ91 have been determined. Specimens with rectangular cross-sections, 1, 2 and 3 mm thick, or with a circular cross-section 6.4 mm in diameter, have been studied. The hardness is generally higher near the edges in all specimens, and more so near the corners of the rectangular specimens. The hardness at the center of the castings is generally lower, due to a coarser solidification microstructure and the concentration of porosity. The evidence confirms that the surface of the castings is harder than the core, but it does not support the concept of a skin with a sharp. and definable boundary. This harder layer is irregular in hardness and depth and is not equally hard on opposite sides of the casting. The mean hardness obtained by integrating the microhardness maps over the entire cross-section increased with decreasing thickness of the bars, and was found to be in good correlation with each bar's yield strength. (c) 2005 Elsevier B.V. All rights reserved.