How to unravel and solve soil fertility problems

This book provides a way for farmers in developing countries to benefit from scientific knowledge on plant nutrition and soil fertility. Specifically, it will help farmers recognise and deal with shortages or excesses of chemical elements.

Matrix elements of U(2n) generators in a multishell spin-orbit basis. I. The del-operator MEs in a two-shell composite Gelfand-Paldus basis

This is the first in a series of three articles which aimed to derive the matrix elements of the U(2n) generators in a multishell spin-orbit basis. This is a basis appropriate to many-electron systems which have a natural partitioning of the orbital space and where also spin-dependent terms are included in the Hamiltonian. The method is based on a new spin-dependent unitary group approach to the many-electron correlation problem due to Gould and Paldus [M. D. Gould and J. Paldus, J. Chem. Phys. 92, 7394, (1990)]. In this approach, the matrix elements of the U(2n) generators in the U(n) x U(2)-adapted electronic Gelfand basis are determined by the matrix elements of a single Ll(n) adjoint tensor operator called the del-operator, denoted by Delta(j)(i) (1 less than or equal to i, j less than or equal to n). Delta or del is a polynomial of degree two in the U(n) matrix E = [E-j(i)]. The approach of Gould and Paldus is based on the transformation properties of the U(2n) generators as an adjoint tensor operator of U(n) x U(2) and application of the Wigner-Eckart theorem. Hence, to generalize this approach, we need to obtain formulas for the complete set of adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis. The nonzero shift coefficients are uniquely determined and may he evaluated by the methods of Gould et al. [see the above reference]. In this article, we define zero-shift adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis which are appropriate to the many-electron problem. By definition, these are proportional to the corresponding two-shell del-operator matrix elements, and it is shown that the Racah factorization lemma applies. Formulas for these coefficients are then obtained by application of the Racah factorization lemma. The zero-shift adjoint reduced Wigner coefficients required for this procedure are evaluated first. All these coefficients are needed later for the multishell case, which leads directly to the two-shell del-operator matrix elements. Finally, we discuss an application to charge and spin densities in a two-shell molecular system. (C) 1998 John Wiley & Sons.

Matrix elements of U(2n) generators in a multishell spin-orbit basis. II. The two-shell nonzero shift ACCs and U(2n) generator MEs

This is the second in a series of articles whose ultimate goal is the evaluation of the matrix elements (MEs) of the U(2n) generators in a multishell spin-orbit basis. This extends the existing unitary group approach to spin-dependent configuration interaction (CI) and many-body perturbation theory calculations on molecules to systems where there is a natural partitioning of the electronic orbital space. As a necessary preliminary to obtaining the U(2n) generator MEs in a multishell spin-orbit basis, we must obtain a complete set of adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis. The zero-shift coefficients were obtained in the first article of the series. in this article, we evaluate the nonzero shift adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis. We then demonstrate that the one-shell versions of these coefficients may be obtained by taking the Gelfand-Tsetlin limit of the two-shell formulas. These coefficients,together with the zero-shift types, then enable us to write down formulas for the U(2n) generator matrix elements in a two-shell spin-orbit basis. Ultimately, the results of the series may be used to determine the many-electron density matrices for a partitioned system. (C) 1998 John Wiley & Sons, Inc.

Matrix elements of U(2n) generators in a multishell spin-orbit basis. III. General formulas

This is the third and final article in a series directed toward the evaluation of the U(2n) generator matrix elements (MEs) in a multishell spin/orbit basis. Such a basis is required for many-electron systems possessing a partitioned orbital space and where spin-dependence is important. The approach taken is based on the transformation properties of the U(2n) generators as an adjoint tensor operator of U(n) x U(2) and application of the Wigner-Eckart theorem. A complete set of adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis (which is appropriate to the many-electron problem) were obtained in the first and second articles of this series. Ln the first article we defined zero-shift coupling coefficients. These are proportional to the corresponding two-shell del-operator matrix elements. See P. J. Burton and and M. D. Gould, J. Chem. Phys., 104, 5112 (1996), for a discussion of the del-operator and its properties. Ln the second article of the series, the nonzero shift coupling coefficients were derived. Having obtained all the necessary coefficients, we now apply the formalism developed above to obtain the U(2n) generator MEs in a multishell spin-orbit basis. The methods used are based on the work of Gould et al. (see the above reference). (C) 1998 John Wiley & Sons, Inc.

Modeling ex vivo hematopoiesis using chemical engineering metaphors

Ex vivo hematopoiesis is increasingly used for clinical applications. Models of ex vivo hematopoiesis are required to better understand the complex dynamics and to optimize hematopoietic culture processes. A general mathematical modeling framework is developed which uses traditional chemical engineering metaphors to describe the complex hematopoietic dynamics. Tanks and tubular reactors are used to describe the (pseudo-) stochastic and deterministic elements of hematopoiesis, respectively. Cells at any point in the differentiation process can belong to either an immobilized, inert phase (quiescent cells) or a mobile, active phase (cycling cells). The model describes five processes: (1) flow (differentiation), (2) autocatalytic formation (growth),(3) degradation (death), (4) phase transition from immobilized to mobile phase (quiescent to cycling transition), and (5) phase transition from mobile to immobilized phase (cycling to quiescent transition). The modeling framework is illustrated with an example concerning the effect of TGF-beta 1 on erythropoiesis. (C) 1998 Published by Elsevier Science Ltd. All rights reserved.

Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

Numerical modelling of chemical effects of magma solidification problems in porous rocks

The solidification of intruded magma in porous rocks can result in the following two consequences: (1) the heat release due to the solidification of the interface between the rock and intruded magma and (2) the mass release of the volatile fluids in the region where the intruded magma is solidified into the rock. Traditionally, the intruded magma solidification problem is treated as a moving interface (i.e. the solidification interface between the rock and intruded magma) problem to consider these consequences in conventional numerical methods. This paper presents an alternative new approach to simulate thermal and chemical consequences/effects of magma intrusion in geological systems, which are composed of porous rocks. In the proposed new approach and algorithm, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with the proposed mass source and physically equivalent heat source. The major advantage in using the proposed equivalent algorithm is that a fixed mesh of finite elements with a variable integration time-step can be employed to simulate the consequences and effects of the intruded magma solidification using the conventional finite element method. The correctness and usefulness of the proposed equivalent algorithm have been demonstrated by a benchmark magma solidification problem. Copyright (c) 2005 John Wiley & Sons, Ltd.

Lanczos subspace filter diagonalization: Homogeneous recursive filtering and a low-storage method for the calculation of matrix elements

We develop a new iterative filter diagonalization (FD) scheme based on Lanczos subspaces and demonstrate its application to the calculation of bound-state and resonance eigenvalues. The new scheme combines the Lanczos three-term vector recursion for the generation of a tridiagonal representation of the Hamiltonian with a three-term scalar recursion to generate filtered states within the Lanczos representation. Eigenstates in the energy windows of interest can then be obtained by solving a small generalized eigenvalue problem in the subspace spanned by the filtered states. The scalar filtering recursion is based on the homogeneous eigenvalue equation of the tridiagonal representation of the Hamiltonian, and is simpler and more efficient than our previous quasi-minimum-residual filter diagonalization (QMRFD) scheme (H. G. Yu and S. C. Smith, Chem. Phys. Lett., 1998, 283, 69), which was based on solving for the action of the Green operator via an inhomogeneous equation. A low-storage method for the construction of Hamiltonian and overlap matrix elements in the filtered-basis representation is devised, in which contributions to the matrix elements are computed simultaneously as the recursion proceeds, allowing coefficients of the filtered states to be discarded once their contribution has been evaluated. Application to the HO2 system shows that the new scheme is highly efficient and can generate eigenvalues with the same numerical accuracy as the basic Lanczos algorithm.

Combined chemical separation of Lu, Hf, Sm, Nd, and REEs from a single rock digest: Precise and accurate isotope Determinations of Lu-Hf and Sm-Nd using multicollector-ICPMS

A combined procedure for separating Lu, Hf, Sm, Nd, and rare earth elements (REEs) from a single sample digest is presented. The procedure consists of the following five steps: (1) sample dissolution via sodium peroxide sintering; (2) separation of the high field strength elements from the REEs and other matrix elements by a HF-free anion-exchange column procedure; (3) purification of Hf on a cation-exchange resin; (4) separation of REEs from other matrix elements by cation exchange; (5) Lu, Sm, and Nd separation from the other REEs by reversed-phase ion chromatography. Analytical reproducibilities of Sm-Nd and Lu-Hf isotope systematics are demonstrated for standard solutions and international rock reference materials. Results show overall good reproducibilities for Sm-Nd systematics independent of the rock type analyzed. For the Lu-Hf systematics, the reproducibility of the parent/daughter ratio is much better for JB-1 (basalt) than for two analyzed felsic crustal rocks (DR-N and an Archaean granitoid). It is demonstrated that this poorer reproducibility of the Lu/Hf ratio is truly caused by sample heterogeneity; thus, results are geologically reasonable.

Characterisation of early Archaean chemical sediments by trace element signatures

Rare earth element (REE) plus yttrium (Y) patterns of modem seawater have characteristic features that can be used as chemical fingerprints. Reliable proxies for marine REE + Y chemistry have been demonstrated from a large geological time span, including Archaean banded iron formation (BIF), stromatolitic limestone, Phanerozoic reef carbonate and Holocene microbialite. Here we present new REE + Y data for two distinct suites of early Archaean (ca. 3.7-3.8 Ga) metamorphosed rocks from southern West Greenland, whose interrelationships, if any, have been much debated in recent literature. The first suite comprises mangetite-quartz BIF, magnetite-carbonate BIF and banded magnetite-rich quartz rock, mostly from the Isua Greenstone Belt (IGB). The REE + Y patterns, particularly diagnostic anomalies (Ce/Ce*, Pr/Pr*), are closely related to those of published seawater proxies. The second suite includes banded quartz-pyroxene-amphibole +/- garnet rocks with minor magnetite from the so-called Akilia Association enclaves (in early Archaean granitoid gneisses) of the coastal region, some 150 km southwest of the IGB. Rocks of this type from one much publicised and highly debated locality (the island of Akilia) have been identified by some workers [Nature 384 (1996) 55; Geochim. Cosmochim. Acta 61 (1997) 2475] as BIF-facies, and their C-13-depleted signature in trace graphite interpreted as a proxy for earliest life on Earth. However, REE + Y patterns of the Akilia Association suite (except for one probably genuine magnetite-rich BIF from Ugpik) are inconsistent with a seawater origin. We agree with published geological and geochemical (including REE) work [Science 296 (2002) 1448] that most of the analysed Akilia rocks are not chemical sediments, and that C-isotopes in such rocks therefore cannot be used as biological proxies. Application of the REE + Y discriminant for the above two rock suites has been facilitated in this study by the use of MC-ICP technique which yields a more complete and precise REE + Y spectrum than was available in many previous studies. (C) 2004 Elsevier B.V. All rights reserved.

The Equilibrium Flux Method for the calculation of flows with non-equilbrium chemical reactions

The Equilibrium Flux Method [1] is a kinetic theory based finite volume method for calculating the flow of a compressible ideal gas. It is shown here that, in effect, the method solves the Euler equations with added pseudo-dissipative terms and that it is a natural upwinding scheme. The method can be easily modified so that the flow of a chemically reacting gas mixture can be calculated. Results from the method for a one-dimensional non-equilibrium reacting flow are shown to agree well with a conventional continuum solution. Results are also presented for the calculation of a plane two-dimensional flow, at hypersonic speed, of a dissociating gas around a blunt-nosed body.

Triple objective team mentoring: achieving learning objectives with chemical engineering students

A sophisticated style of mentoring has been found to be essential to support engineering student teams undertaking technically demanding, real-world problems as part of a Project-Centred Curriculum (PCC) at The University of Queensland. The term ‘triple-objective’ mentoring was coined to define mentoring that addresses not only the student’s technical goal achievement but also their time and team management. This is achieved through a number of formal mentor meetings that are informed by a confidential instrument which requires students to individually reflect on team processes prior to the meeting, and a checklist of technical requirements against which the interim student team progress and achievements are assessed. Triple-objective mentoring requires significant time input and coordination by the academic but has been shown to ensure effective student team work and learning undiminished by team dysfunction. Student feedback shows they value the process and agree that the tools developed to support the process are effective in developing and assessing team work and skills with average scores mostly above 3 on a four point scale.

Chemical treatment of Escherichia coli. II. Direct extraction of recombinant protein from cytoplasmic inclusion bodies in intact cells

A method is presented for the direct extraction of the recombinant protein Long-R-3-IGF-I from inclusion bodies located in the cytoplasm of intact Escherichia coli cells. Chemical treatment with 6M urea, 3 mM EDTA, and 20 mM dithiothreitol (DTT) at pH 9.0 proved an effective combination for extracting recombinant protein from intact cells. Comparable levels of Long-R-3-IGF-I were recovered by direct extraction as achieved by in vitro dissolution following mechanical disruption. However, the purity of directly extracted recombinant protein was lower due to contamination by bacterial cell components. The kinetics of direct extraction are described using a first-order equation with the time constant of 3 min. Urea appears important for permeabilization of the cell and dissolution of the inclusion body. Conversely, EDTA is involved in permeabilization of the cell wall and DTT enhances protein release. pH proved to be important with lower levels of protein release achieved at low pH values (

Prediction of the chemical composition of lamb carcasses from multi-frequency impedance data

Multi-frequency bioimpedance analysis (MFBIA) was used to determine the impedance, reactance and resistance of 103 lamb carcasses (17.1-34.2 kg) immediately after slaughter and evisceration. Carcasses were halved, frozen and one half subsequently homogenized and analysed for water, crude protein and fat content. Three measures of carcass length were obtained. Diagonal length between the electrodes (right side biceps femoris to left side of neck) explained a greater proportion of the variance in water mass than did estimates of spinal length and was selected for use in the index L-2/Z to predict the mass of chemical components in the carcass. Use of impedance (Z) measured at the characteristic frequency (Z(c)) instead of 50 kHz (Z(50)) did not improve the power of the model to predict the mass of water, protein or fat in the carcass. While L-2/Z(50) explained a significant proportion of variation in the masses of body water (r(2) 0.64), protein (r(2) 0.34) and fat (r(2) 0.35), its inclusion in multi-variate indices offered small or no increases in predictive capacity when hot carcass weight (HCW) and a measure of rib fat-depth (GR) were present in the model. Optimized equations were able to account for 65-90 % of the variance observed in the weight of chemical components in the carcass. It is concluded that single frequency impedance data do not provide better prediction of carcass composition than can be obtained from measures of HCW and GR. Indices of intracellular water mass derived from impedance at zero frequency and the characteristic frequency explained a similar proportion of the variance in carcass protein mass as did the index L-2/Z(50).