Time evolution in the unimolecular master equation at low temperatures: full spectral solution with scalable iterative methods and high precision

A new method is presented to determine an accurate eigendecomposition of difficult low temperature unimolecular master equation problems. Based on a generalisation of the Nesbet method, the new method is capable of achieving complete spectral resolution of the master equation matrix with relative accuracy in the eigenvectors. The method is applied to a test case of the decomposition of ethane at 300 K from a microcanonical initial population with energy transfer modelled by both Ergodic Collision Theory and the exponential-down model. The fact that quadruple precision (16-byte) arithmetic is required irrespective of the eigensolution method used is demonstrated. (C) 2001 Elsevier Science B.V. All rights reserved.

Calculation of product state distributions from resonance decay via Lanczos subspace filter diagonalization: Application to HO2

Resonance phenomena associated with the unimolecular dissociation of HO2 have been investigated quantum-mechanically by the Lanczos homogeneous filter diagonalization (LHFD) method. The calculated resonance energies, rates (widths), and product state distributions are compared to results from an autocorrelation function-based filter diagonalization (ACFFD) method. For calculating resonance wave functions via ACFFD, an analytical expression for the expansion coefficients of the modified Chebyshev polynomials is introduced. Both dissociation rates and product state distributions of O-2 show strong fluctuations, indicating the dissociation of HO2 is essentially irregular. (C) 2001 American Institute of Physics.

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).

Limits to squeezing in the degenerate optical parametric oscillator

We develop a systematic theory of quantum fluctuations in the driven optical parametric oscillator, including the region near threshold. This allows us to treat the limits imposed by nonlinearities to quantum squeezing and noise reduction in this nonequilibrium quantum phase transition. In particular, we compute the squeezing spectrum near threshold and calculate the optimum value. We find that the optimal noise reduction occurs at different driving fields, depending on the ratio of damping rates. The largest spectral noise reductions are predicted to occur with a very high-Q second-harmonic cavity. Our analytic results agree well with stochastic numerical simulations. We also compare the results obtained in the positive-P representation, as a fully quantum-mechanical calculation, with the truncated Wigner phase-space equation, also known as the semiclassical theory.

Wavelet-based image segment representation

An efficient representation method for arbitrarily shaped image segments is proposed. This method includes a smart way to select wavelet basis to approximate the given image segment, with improved image quality and reduced computational load.

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.

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.

Hydrofluoric acid pre-treatment for improving 13C CPMAS NMR spectral quality of forest soils in south-east Queensland, Australia

Hydrofluoric acid (HF) was used to pre-treat forest soils of south-east Queensland for assessing the effectiveness of iron (Fe) removal, carbon (C) composition using C-13 cross-polarisation (CP) with magic-angle-spinning (MAS) nuclear magnetic resonance (NMR) before and after the HF pre-treatment, and the improvement of C-13 CPMAS NMR spectra. Soil samples were collected from 4 experimental sites of different soil types, harvest residue management or prescribed burning, and tree species. More than 86% of Fe was in all soil types removed by the HF treatment. The C-13 NMR spectral quality was improved with increased resolution, especially in the alkyl C and O-alkyl C regions, and reduced NMR run-time (1-5 h per sample compared with >20 h per sample without the pre-treatment). The C composition appeared to alter slightly after the pre-treatment, but this might be largely due to improved spectrometer conditions and increased resolution leading to more accurate NMR spectral integration. Organic C recovery after HF pre-treatment varied with soil types and forest management, and soluble soil organic matter (SOM) could be lost during the pre-treatment. The Fourier Transform-Infrared (FT-IR) spectra of HF extracts indicated the preferential removal of carboxylic C groups during the pre-treatment, but this could also be due to adsorbed water on the mineral matter. The NMR spectra revealed some changes in C composition and quality due to residue management and decomposition. Overall, the HF treatment was a useful pre-treatment for obtaining semi-quantitative C-13 CPMAS NMR spectra of subtropical Australian forest soils.

Color Models and Illumination Changes

Color model representation allows characterizing in a quantitative manner, any defined color spectrum of visible light, i.e. with a wavelength between 400nm and 700nm. To accomplish that, each model, or color space, is associated with a function that allows mapping the spectral power distribution of the visible electromagnetic radiation, in a space defined by a set of discrete values that quantify the color components composing the model. Some color spaces are sensitive to changes in lighting conditions. Others assure the preservation of certain chromatic features, remaining immune to these changes. Therefore, it becomes necessary to identify the strengths and weaknesses of each model in order to justify the adoption of color spaces in image processing and analysis techniques. This chapter will address the topic of digital imaging, main standards and formats. Next we will set the mathematical model of the image acquisition sensor response, which enables assessment of the various color spaces, with the aim of determining their invariance to illumination changes.

Eigenvalue computations in the context of data-sparse approximations of integral operators

In this work, we consider the numerical solution of a large eigenvalue problem resulting from a finite rank discretization of an integral operator. We are interested in computing a few eigenpairs, with an iterative method, so a matrix representation that allows for fast matrix-vector products is required. Hierarchical matrices are appropriate for this setting, and also provide cheap LU decompositions required in the spectral transformation technique. We illustrate the use of freely available software tools to address the problem, in particular SLEPc for the eigensolvers and HLib for the construction of H-matrices. The numerical tests are performed using an astrophysics application. Results show the benefits of the data-sparse representation compared to standard storage schemes, in terms of computational cost as well as memory requirements.

Error bounds for low-rank approximations of the first exponential integral kernel

A hierarchical matrix is an efficient data-sparse representation of a matrix, especially useful for large dimensional problems. It consists of low-rank subblocks leading to low memory requirements as well as inexpensive computational costs. In this work, we discuss the use of the hierarchical matrix technique in the numerical solution of a large scale eigenvalue problem arising from a finite rank discretization of an integral operator. The operator is of convolution type, it is defined through the first exponential-integral function and, hence, it is weakly singular. We develop analytical expressions for the approximate degenerate kernels and deduce error upper bounds for these approximations. Some computational results illustrating the efficiency and robustness of the approach are presented.

Grid structure impact in sparse point representation of derivatives

In the Sparse Point Representation (SPR) method the principle is to retain the function data indicated by significant interpolatory wavelet coefficients, which are defined as interpolation errors by means of an interpolating subdivision scheme. Typically, a SPR grid is coarse in smooth regions, and refined close to irregularities. Furthermore, the computation of partial derivatives of a function from the information of its SPR content is performed in two steps. The first one is a refinement procedure to extend the SPR by the inclusion of new interpolated point values in a security zone. Then, for points in the refined grid, such derivatives are approximated by uniform finite differences, using a step size proportional to each point local scale. If required neighboring stencils are not present in the grid, the corresponding missing point values are approximated from coarser scales using the interpolating subdivision scheme. Using the cubic interpolation subdivision scheme, we demonstrate that such adaptive finite differences can be formulated in terms of a collocation scheme based on the wavelet expansion associated to the SPR. For this purpose, we prove some results concerning the local behavior of such wavelet reconstruction operators, which stand for SPR grids having appropriate structures. This statement implies that the adaptive finite difference scheme and the one using the step size of the finest level produce the same result at SPR grid points. Consequently, in addition to the refinement strategy, our analysis indicates that some care must be taken concerning the grid structure, in order to keep the truncation error under a certain accuracy limit. Illustrating results are presented for 2D Maxwell's equation numerical solutions.