Internal standardization in graphite furnace atomic absorption spectrometry: Comparative use of As and Ge to minimize matrix effects on Se determination in milk

Arsenic and germanium have been evaluated as internal standards to minimize matrix effects on the direct determination of selenium in milk by graphite furnace atomic absorption spectrometry (GFAAS) using tubes with integrated platform, pre-treated with W together with I'd as chemical modifier. The efficiency of As and Ge as internal standards for 25 mu g L-1 Se plus 500 mu g (L)-1 As or Ge in diluted (1 + 9 v/v) milk plus 1.0% (v/v) HNO3 was evaluated by means of correlation graphs plotted from the normalized absorbance signals (n = 20) of internal standard (axis gamma) versus analyte (axis x). The equations that describe the linear regression were: A(As)= - 0.004 +/- 0.019 +/- 1.02 + 0.019 A(Se) (r=0.9967 +/- 0.005); A(Ge)= - 0.0 17 +/- 0.015 + 1.01 +/- 0.015 A(Se) (r = 0.9978 +/- 0.004). Samples and reference solutions were automatically spiked with 500 mu g L-1 Ge or As and 1.0% (v/v) HNO3 by the autosampler. For 20 mu L of aqueous standard solutions, analytical curves in the 5.00-40.0 mu g L-1 Se range were established using the ratio of Se absorbance to internal standard absorbance (A(Se)A(IS)) versus analyte concentration, and good linear correlations were obtained. The characteristic mass was 40 pg Se. Limits of detection were 0.55 and 0.40 mu g L-1 with As and Ge as the internal standard, respectively. Relative standard deviations (RSD) for a sample containing 25 mu g L-1 Se were 1.2% and 1.0% (n = 12) using As and Ge, respectively. The RSD without internal standardization was about 6%. The accuracy of the proposed method was evaluated by an addition-recovery experiment and all recovered values were in the 99-105% range with IS and in the 70-80% range without IS. Using Ge as the internal standard, results of analysis of standard reference materials were in agreement with certified values at a 95% confidence level. The selenium concentration for 10 analyzed milk samples varied from 5.0 to 20 mu g L-1. (c) 2005 Elsevier B.V. All rights reserved.

Two-matrix string model as constrained (2+1)-dimensional integrable system

We show that the 2-matrix string model corresponds to a coupled system of 2 + 1-dimensional KP and modified KP ((m)KP2+1) integrable equations subject to a specific symmetry constraint. The latter together with the Miura-Konopelchenko map for (m)KP2+1 are the continuum incarnation of the matrix string equation. The (m)KP2+1 Miura and Backhand transformations are natural consequences of the underlying lattice structure. The constrained (m)KP2+1 system is equivalent to a 1 + 1-dimensional generalized KP-KdV hierarchy related to graded SL(3,1). We provide an explicit representation of this hierarchy, including the associated W(2,1)-algebra of the second Hamiltonian structure, in terms of free currents.

Tau-functions and dressing transformations for zero-curvature affine integrable equations

The solutions of a large class of hierarchies of zero-curvature equations that includes Toda- and KdV-type hierarchies are investigated. All these hierarchies are constructed from affine (twisted or untwisted) Kac-Moody algebras g. Their common feature is that they have some special vacuum solutions corresponding to Lax operators lying in some Abelian (up to the central term) subalgebra of g; in some interesting cases such subalgebras are of the Heisenberg type. Using the dressing transformation method, the solutions in the orbit of those vacuum solutions are constructed in a uniform way. Then, the generalized tau-functions for those hierarchies are defined as an alternative set of variables corresponding to certain matrix elements evaluated in the integrable highest-weight representations of g. Such definition of tau-functions applies for any level of the representation, and it is independent of its realization (vertex operator or not). The particular important cases of generalized mKdV and KdV hierarchies as well as the Abelian and non-Abelian affine Toda theories are discussed in detail. © 1997 American Institute of Physics.

Equivalence of the Duffin-Kemmer-Petiau and Klein-Gordon-Fock equations

A strict proof of the equivalence of the Duffin-Kemmer-Petiau and Klein-Gordon Fock theories is presented for physical S-matrix elements in the case of charged scalar particles minimally interacting with an external or quantized electromagnetic field. The Hamiltonian canonical approach to the Duffin - Kemmer Petiau theory is first developed in both the component and the matrix form. The theory is then quantized through the construction of the generating functional for the Green's functions, and the physical matrix elements of the S-matrix are proved to be relativistic invariants. The equivalence of the two theories is then proved for the matrix elements of the scattered scalar particles using the reduction formulas of Lehmann, Symanzik, and Zimmermann and for the many-photon Green's functions.

Stability analysis of Crank-Nicolson and Euler schemes for time-dependent diffusion equations

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Preconditioning strategies for the numerical solution of convection-diffusion partial differential equations

Il trattamento numerico dell'equazione di convezione-diffusione con le relative condizioni al bordo, comporta la risoluzione di sistemi lineari algebrici di grandi dimensioni in cui la matrice dei coefficienti è non simmetrica. Risolutori iterativi basati sul sottospazio di Krylov sono ampiamente utilizzati per questi sistemi lineari la cui risoluzione risulta particolarmente impegnativa nel caso di convezione dominante. In questa tesi vengono analizzate alcune strategie di precondizionamento, atte ad accelerare la convergenza di questi metodi iterativi. Vengono confrontati sperimentalmente precondizionatori molto noti come ILU e iterazioni di tipo inner-outer flessibile. Nel caso in cui i coefficienti del termine di convezione siano a variabili separabili, proponiamo una nuova strategia di precondizionamento basata sull'approssimazione, mediante equazione matriciale, dell'operatore differenziale di convezione-diffusione. L'azione di questo nuovo precondizionatore sfrutta in modo opportuno recenti risolutori efficienti per equazioni matriciali lineari. Vengono riportati numerosi esperimenti numerici per studiare la dipendenza della performance dei diversi risolutori dalla scelta del termine di convezione, e dai parametri di discretizzazione.

Interface discontinuity factors in the modal eigenspace of the multigroup diffusion matrix

Interface discontinuity factors based on the Generalized Equivalence Theory are commonly used in nodal homogenized diffusion calculations so that diffusion average values approximate heterogeneous higher order solutions. In this paper, an additional form of interface correction factors is presented in the frame of the Analytic Coarse Mesh Finite Difference Method (ACMFD), based on a correction of the modal fluxes instead of the physical fluxes. In the ACMFD formulation, implemented in COBAYA3 code, the coupled multigroup diffusion equations inside a homogenized region are reduced to a set of uncoupled modal equations through diagonalization of the multigroup diffusion matrix. Then, physical fluxes are transformed into modal fluxes in the eigenspace of the diffusion matrix. It is possible to introduce interface flux discontinuity jumps as the difference of heterogeneous and homogeneous modal fluxes instead of introducing interface discontinuity factors as the ratio of heterogeneous and homogeneous physical fluxes. The formulation in the modal space has been implemented in COBAYA3 code and assessed by comparison with solutions using classical interface discontinuity factors in the physical space

A dynamic programming approach to the formulation and solution of finite element equations

A method for formulating and algorithmically solving the equations of finite element problems is presented. The method starts with a parametric partition of the domain in juxtaposed strips that permits sweeping the whole region by a sequential addition (or removal) of adjacent strips. The solution of the difference equations constructed over that grid proceeds along with the addition removal of strips in a manner resembling the transfer matrix approach, except that different rules of composition that lead to numerically stable algorithms are used for the stiffness matrices of the strips. Dynamic programming and invariant imbedding ideas underlie the construction of such rules of composition. Among other features of interest, the present methodology provides to some extent the analyst's control over the type and quantity of data to be computed. In particular, the one-sweep method presented in Section 9, with no apparent counterpart in standard methods, appears to be very efficient insofar as time and storage is concerned. The paper ends with the presentation of a numerical example

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.

Optical Bloch equations and enhanced decay of Rabi oscillations in strong driving fields

Optical Bloch equations are widely used for describing dynamics in a system consisting molecules, electromagnetic waves, and a thermal bath. We analyze applicability of these equations to a single molecule imbedded in a solid matrix. Classical Bloch equations and the limits of their applicability are derived from more general master equations. Simple and intuitively appealing picture based on stochastic Bloch equations shows that at low temperatures, contrary to common believes, a strong driving field can not only suppress but can also increase decay rates of Rabi oscillations. A physical system where predicted effects can be observed experimentally is suggested. (c) 2005 Elsevier B.V. All rights reserved.

Boundary-Value Problems for almost Nonlinear Singularly Perturbed Systems of Ordinary Differential Equations

A boundary-value problems for almost nonlinear singularly perturbed systems of ordinary differential equations are considered. An asymptotic solution is constructed under some assumption and using boundary functions and generalized inverse matrix and projectors.

Some Properties of Mittag-Leffler Functions and Matrix-Variate Analogues: A Statistical Perspective

Mathematical Subject Classification 2010:26A33, 33E99, 15A52, 62E15.

Numerical Solution of Fractional Diffusion-Wave Equation with two Space Variables by Matrix Method

Mathematics Subject Classi¯cation 2010: 26A33, 65D25, 65M06, 65Z05.

Matrix-Variate Statistical Distributions and Fractional Calculus

MSC 2010: 15A15, 15A52, 33C60, 33E12, 44A20, 62E15 Dedicated to Professor R. Gorenflo on the occasion of his 80th birthday