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.

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 resonances and product state distributions for the unimolecular dissociation of H2S

Resonance phenomena associated with the unimolecular dissociation of H2S --> SH + H have been investigated quantum mechanically by the Lanczos homogeneous filter diagonalization method using a newly developed potential energy surface (J. Chem. Phys. 2001, 114, 320). Resonance energies, widths (rates), and product state distributions have been obtained. Both dissociation rates and product state distributions of SH show, strong fluctuations, indicating that the dissociation of H2S is essentially irregular. Statistical analysis of neighboring level spacing and width distributions also confirms this behavior. The dissociation rates and product state distributions are compared to the predictions of quantum phase space theory.

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.

Absorption and dispersion by a multiple driven two-level atom

We investigate the absorption and dispersion properties of a two-level atom driven by a polychromatic field. The driving field is composed of a strong resonant (carrier) frequency component and a large number of symmetrically detuned sideband fields (modulators). A rapid increase in the absorption at the central frequency and the collapse of the response of the system from multiple frequencies to a single frequency are predicted to occur when the Rabi frequency of the modulating fields is equal to the Rabi frequency of the carrier field. These are manifestations of the undressing or a disentanglement of the atomic and driving field states, that leads to a collapse of the atom to its ground state. Our calculation permits consideration of the question of the undressing of the driven atom by a multiple-modulated field and the predicted spectra offer a method of observing undressing. Moreover, we find that the absorption and dispersion spectra split into multiplets whose structures depend on the Rabi frequency of the modulating fields. The spectral features can jump between different resonance frequencies by changing the Rabi frequency of the modulating fields or their initial phases, which can have potential applications as a quantum frequency filter.

Mineral chemistry, whole-rock compositions, and petrogenesis of leg 176 gabbros: Data and discussion

We report mineral chemistry, whole-rock major element compositions, and trace element analyses on Hole 735B samples drilled and selected during Leg 176. We discuss these data, together with Leg 176 shipboard data and Leg 118 sample data from the literature, in terms of primary igneous petrogenesis. Despite mineral compositional variation in a given sample, major constituent minerals in Hole 735B gabbroic rocks display good chemical equilibrium as shown by significant correlations among Mg# (= Mg/[Mg+Fe2+]) of olivine, clinopyroxene, and orthopyroxene and An (=Ca/[Ca+Na]) of plagioclase. This indicates that the mineral assemblages olivine + plagioclase in troctolite, plagioclase + clinopyroxene in gabbro, plagioclases + clinopyroxene + olivine in olivine gabbro, and plagioclase + clinopyroxene + olivine + orthopyroxene in gabbronorite, and so on, have all coprecipitated from their respective parental melts. Fe-Ti oxides (ilmenite and titanomagnetite), which are ubiquitous in most of these rocks, are not in chemical equilibrium with olivine, clinopyroxene, and plagioclase, but precipitated later at lower temperatures. Disseminated oxides in some samples may have precipitated from trapped Fe-Ti–rich melts. Oxides that concentrate along shear bands/zones may mark zones of melt coalescence/transport expelled from the cumulate sequence as a result of compaction or filter pressing. Bulk Hole 735B is of cumulate composition. The most primitive olivine, with Fo = 0.842, in Hole 735B suggests that the most primitive melt parental to Hole 735B lithologies must have Mg# ≤ 0.637, which is significantly less than Mg# = 0.714 of bulk Hole 735B.

Simulation of the micro-physics of rocks using LSMearth

The particle-based Lattice Solid Model (LSM) was developed to provide a basis to study the physics of rocks and the nonlinear dynamics of earthquakes (MORA and PLACE, 1994; PLACE and MORA, 1999). A new modular and flexible LSM approach has been developed that allows different microphysics to be easily included in or removed from the model. The approach provides a virtual laboratory where numerical experiments can easily be set up and all measurable quantities visualised. The proposed approach provides a means to simulate complex phenomena such as fracturing or localisation processes, and enables the effect of different micro-physics on macroscopic behaviour to be studied. The initial 2-D model is extended to allow three-dimensional simulations to be performed and particles of different sizes to be specified. Numerical bi-axial compression experiments under different confining pressure are used to calibrate the model. By tuning the different microscopic parameters (such as coefficient of friction, microscopic strength and distribution of grain sizes), the macroscopic strength of the material and can be adjusted to be in agreement with laboratory experiments, and the orientation of fractures is consistent with the theoretical value predicted based on Mohr-Coulomb diagram. Simulations indicate that 3-D numerical models have different macroscopic properties than in 2-D and, hence, the model must be recalibrated for 3-D simulations. These numerical experiments illustrate that the new approach is capable of simulating typical rock fracture behaviour. The new model provides a basis to investigate nucleation, rupture and slip pulse propagation in complex fault zones without the previous model limitations of a regular low-level surface geometry and being restricted to two-dimensions.

Diverse synaptic mechanisms generate direction selectivity in the rabbit retina

The synaptic conductance of the On-Off direction-selective ganglion cells was measured during visual stimulation to determine whether the direction selectivity is a property of the circuitry presynaptic to the ganglion cells or is generated by postsynaptic interaction of excitatory and inhibitory inputs. Three synaptic asymmetries were identified that contribute to the generation of direction-selective responses: (1) a presynaptic mechanism producing stronger excitation in the preferred direction, (2) a presynaptic mechanism producing stronger inhibition in the opposite direction, and (3) postsynaptic interaction of excitation with spatially offset inhibition. Although the on- and off-responses showed the same directional tuning, the off-response was generated by all three mechanisms, whereas the on- response was generated primarily by the two presynaptic mechanisms. The results indicate that, within a single neuron, different strategies are used within distinct dendritic arbors to accomplish the same neural computation.

Selective ingestion of pelagic versus benthic algae by the cockle Cerastoderma edule (Linne, 1758)

The pre-ingestive selection of microphytobenthic algae by the cockle Cerastoderma edule was studied in comparison with diets containing the pelagic diatom Phaeodactylum tricomutum. Treatments with the different diets covered a range of seston concentrations and organic content similar to field conditions. Rejection rates of C. edule exposed to the different treatments were significantly correlated with the concentration of total particulate matter. No significant differences in total rejection rates were found between pelagic and benthic diets. Organic rejection rate was significantly correlated with particulate organic matter of the treatments and no significant differences were found between both diets. Selection efficiency was significantly correlated with particulate organic matter concentration in both diets and no significant differences were found between the diets. Analysis of the pseudofeces composition by flow cytometry from cockles exposed to a mixed diet of microphytobenthic algae and P. tricornutum, showed a preferential ingestion of the pelagic diatom. Benthic species, such as small pennates and Navicula sp., were preferentially ingested in comparison to larger microphytobenthic species. The largest microphytobenthic species, Cylindrotheca sp., was significantly rejected. In general, C. edule is an opportunistic filter feeder that takes advantage of both pelagic and benthic algal cells.

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.

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.

Discrete cyanide-bridged mixed-valence Co/Fe complexes: Outer-sphere redox behaviour

The outer-sphere redox behaviour of a series of [LnCoIII-NCFeII(CN)(5)](-) (L-n = n-membered pentadentate aza-macrocycle) complexes have been studied as a function of pH and oxidising agent. All the dinuclear complexes show a double protonation process at pH approximate to 2 that produces a shift in their UV/Vis spectra. Oxidation of the different non-protonated and diprotonated complexes has been carried out with peroxodisulfate, and of the non-protonated complexes also with trisoxalatocobaltate(III). The results are in agreement with predictions from the Marcus theory. The oxidation of [Fe(phen)(3)](3+) and [IrCl6](2-) is too fast to be measured, although for the latter the transient observation of the process has been achieved at pH = 0. The study of the kinetics of the outer-sphere redox process, with the S2O82- and [Co(ox)(3)](3-) oxidants, has been carried out as a function of pH, temperature, and pressure. As a whole, the values found for the activation volumes, entropies, and enthalpies are in the following margins, for the diprotonated and non-protonated dinuclear complexes, respectively: DeltaV(not equal) from 11 to 13 and 15 to 20 cm(3) mol(-1); DeltaS(not equal) from 110 to 30 and -60 to -90 J K-1 mol(-1); DeltaH(not equal) from 115 to 80 and 50 to 65 kJ.mol(-1). The thermal activation parameters are clearly dominated by the electrostriction occurring on outer-sphere precursor formation, while the trends found for the values of the volume of activation indicate an important degree of tuning due to the charge distribution during the electron transfer process. The special arrangement on the amine ligands in the isomer trans[(L14CoNCFeII)-N-III(CN)(5)](-) accounts for important differences in solvent-assisted hydrogen bonding occurring within the outer-sphere redox process, as has been established in redox reactions of similar compounds. ((C) Wiley-VCH Verlag GmbH & Co. KGaA, 69451 Weinheim, Germany, 2003).

Measuring the accumulated hazards of smoking: global and regional estimates for 2000

Objective: Current prevalence of smoking, even where data are available, is a poor proxy for cumulative hazards of smoking, which depend on several factors including the age at which smoking began, duration of smoking, number of cigarettes smoked per day, degree of inhalation, and cigarette characteristics such as tar and nicotine content or filter type. Methods: We extended the Peto-Lopez smoking impact ratio method to estimate accumulated hazards of smoking for different regions of the world. Lung cancer mortality data were obtained from the Global Burden of Disease mortality database. The American Cancer Society Cancer Prevention Study, phase 11 (CPS-II) with follow up for the years 1982 to 1988 was the reference population. For the global application of the method, never-smoker lung cancer mortality rates were chosen based on the estimated use of coal for household energy in each region. Results: Men in industrialised countries of Europe, North America, and the Western Pacific had the largest accumulated hazards of smoking. Young and middle age males in many regions of the developing world also had large smoking risks. The accumulated hazards of smoking for women were highest in North America followed by Europe. Conclusions: In the absence of detailed data on smoking prevalence and history, lung cancer mortality provides a robust indicator of the accumulated hazards of smoking. These hazards in developing countries are currently more concentrated among young and middle aged males.