Effects of phosphorus intake on phosphorus flow in growing pigs: Application and comparison of two models

A comparison of the models of Vitti et al. (2000, J. Anim. Sci. 78, 2706-2712) and Fernandez (1995c, Livest. Prod. Sci. 41, 255-261) was carried out using two data sets on growing pigs as input. The two models compared were based on similar basic principles, although their aims and calculations differed. The Vitti model employs the rate:state formalism and describes phosphorus (P) flow between four pools representing P content in gut, blood, bone and soft tissue in growing goats. The Fernandez model describes flow and fractional recirculation between P pools in gut, blood and bone in growing pigs. The results from both models showed similar trends for P absorption from gut to blood and net retention in bone with increasing P intake, with the exception of the 65 kg results from Date Set 2 calculated using the FernAndez model. Endogenous loss from blood back to gut increased faster with increasing P intake in the FernAndez than in the Vitti model for Data Set 1. However, for Data Set 2, endogenous loss increased with increasing P intake using the Vitti model, but decreased when calculated using the FernAndez model. Incorporation of P into bone was not influenced by intake in the FernAndez model, while in the Vitti model there was an increasing trend. The FernAndez model produced a pattern of decreasing resorption in bone with increasing P intake, with one of the data sets, which was not observed when using the Vitti model. The pigs maintained their P homeostasis in blood by regulation of P excretion in urine. (c) 2005 Elsevier Ltd. All rights reserved.

From a discrete to a continuum model of cell dynamics in one dimension

Multiscale modeling is emerging as one of the key challenges in mathematical biology. However, the recent rapid increase in the number of modeling methodologies being used to describe cell populations has raised a number of interesting questions. For example, at the cellular scale, how can the appropriate discrete cell-level model be identified in a given context? Additionally, how can the many phenomenological assumptions used in the derivation of models at the continuum scale be related to individual cell behavior? In order to begin to address such questions, we consider a discrete one-dimensional cell-based model in which cells are assumed to interact via linear springs. From the discrete equations of motion, the continuous Rouse [P. E. Rouse, J. Chem. Phys. 21, 1272 (1953)] model is obtained. This formalism readily allows the definition of a cell number density for which a nonlinear "fast" diffusion equation is derived. Excellent agreement is demonstrated between the continuum and discrete models. Subsequently, via the incorporation of cell division, we demonstrate that the derived nonlinear diffusion model is robust to the inclusion of more realistic biological detail. In the limit of stiff springs, where cells can be considered to be incompressible, we show that cell velocity can be directly related to cell production. This assumption is frequently made in the literature but our derivation places limits on its validity. Finally, the model is compared with a model of a similar form recently derived for a different discrete cell-based model and it is shown how the different diffusion coefficients can be understood in terms of the underlying assumptions about cell behavior in the respective discrete models.

Calculation of converged rovibrational energies and partition function for methane using vibrational-rotational configuration interaction

The rovibration partition function of CH4 was calculated in the temperature range of 100-1000 K using well-converged energy levels that were calculated by vibrational-rotational configuration interaction using the Watson Hamiltonian for total angular momenta J=0-50 and the MULTIMODE computer program. The configuration state functions are products of ground-state occupied and virtual modals obtained using the vibrational self-consistent field method. The Gilbert and Jordan potential energy surface was used for the calculations. The resulting partition function was used to test the harmonic oscillator approximation and the separable-rotation approximation. The harmonic oscillator, rigid-rotator approximation is in error by a factor of 2.3 at 300 K, but we also propose a separable-rotation approximation that is accurate within 2% from 100 to 1000 K. (C) 2004 American Institute of Physics.

Protonation of phenylboronic acid: comparison of G3B3 and G2MP2 methods

G3B3 and G2MP2 calculations using Gaussian 03 have been carried out to investigate the protonation preferences for phenylboronic acid. All nine heavy atoms have been protonated in turn. With both methodologies, the two lowest protonation energies are obtained with the proton located either at the ipso carbon atom or at a hydroxyl oxygen atom. Within the G3B3 formalism, the lowest-energy configuration by 4.3 kcal . mol(-1) is found when the proton is located at the ipso carbon, rather than at the electronegative oxygen atom. In the resulting structure, the phenyl ring has lost a significant amount of aromaticity. By contrast, calculations with G2MP2 show that protonation at the hydroxyl oxygen atom is favored by 7.7 kcal . mol(-1). Calculations using the polarizable continuum model (PCM) solvent method also give preference to protonation at the oxygen atom when water is used as the solvent. The preference for protonation at the ipso carbon found by the more accurate G3B3 method is unexpected and its implications in Suzuki coupling are discussed. (C) 2006 Wiley Periodicals, Inc.

Theoretical rovibrational line intensities in the electronic ground state of ozone

First-principles calculations of absolute line intensities and rovibrational energies of ozone (O-16(3)) are reported using potential energy and electric dipole moment functions calculated by the internally contracted MRCI approach. The rovibrational energies and eigenfunctions (up to about 8500 cm(-1) and J = 64) were obtained variationally with an exact Hamiltonian in internal valence coordinates. More than 4.8 x 10(6) electric dipole transition matrix elements were calculated for the absolute rovibrational line intensities. They are compared with the values of the HITRAN database. The purely rotational absolute line intensities in the (000) state and the rovibrational intensities for the (001)-(000) band agree to within about 0.3 to 1% for the (0 10)-(000) band to within about 3 to 4%. Excellent agreement with experiment is also achieved for low-lying overtone and combination bands. Inconsistencies are found for the (100)-(000) band overlapping with the antisymmetric stretching fundamental and also for the (002)-(000) antisymmetric stretching overtone. The generated dipole moment function can be used for predicting the absorption intensities in any of the heavier isotopomers, hot bands or the rates of spontaneous emission.

Structure, and thermal, magnetic and chemical properties of copper(II)iodoacetate monohydrate

Reaction of iodoacetic acid with cupric carbonate in water in dimmed light yields green Cu(ICH2COO)(2 center dot)H2O (1). From X-ray crystallography, it is found to be a tetra-acetato bridged copper(II) dimer with the water molecules occupying the apical positions. In thermogravimetry, the coordinated water molecules are lost in the temperature range 50-100 degrees C. From magnetic susceptibility measurements in the temperature range 300-1.8 K, the exchange coupling constant J is found to be -142(1) cm(-1) and g = 2.18(2) with the spin Hamiltonian H = -2J{S-Cu1 center dot S-Cu2}. It reacts with 2,2'-bipyridine (bpy) to yield [Cu(bpy)(2)I]I. It oxidises thiophenol to Ph-S-S-Ph under dry N-2 atmosphere.

Large vibrational variational calculations using 'multimode' and an iterative diagonalization technique

We have favoured the variational (secular equation) method for the determination of the (ro-) vibrational energy levels of polyatomic molecules. We use predominantly the Watson Hamiltonian in normal coordinates and an associated given potential in the variational code 'Multimode'. The dominant cost is the construction and diagonalization of matrices of ever-increasing size. Here we address this problem, using pertubation theory to select dominant expansion terms within the Davidson-Liu iterative diagonalization method. Our chosen example is the twelve-mode molecule methanol, for which we have an ab initio representation of the potential which includes the internal rotational motion of the OH group relative to CH3. Our new algorithm allows us to obtain converged energy levels for matrices of dimensions in excess of 100 000.

Two new mu-(1,3-azido)-bridged polymers: alternating single and double bridges in a 1D nickel(II) complex and uniform bridge in a 2D copper(II) complex: Syntheses, single-crystal structures and magnetic studies

Two polymeric azido bridged complexes [Ni2L2(N-3)(3)](n)(ClO4). (1) and [Cu(bpdS)(2)(N-3)],(ClO4),(H2O)(2.5n) (2) [L = Schiff base, obtained from the condensation of pyridine-2-aldehyde with N,N,2,2-tetramethyl-1,3-propanediamine; bpds = 4,4'-bipyridyl disulfide] have been synthesized and their crystal structures have been determined. Complex 1, C26H42ClN15Ni2O4, crystallizes in a triclinic system, space group P1 with a 8.089(13), b = 9.392(14), c = 12.267(18) angstrom, a = 107.28(l), b 95.95(1), gamma = 96.92(1)degrees and Z = 2; complex 2, C20H21ClCuN7O6.5S4, crystallizes in an orthorhombic system, space group Pnna with a = 10.839(14), b = 13.208(17), c = 19.75(2) angstrom and Z = 4. The crystal structure of I consists of 1D polymers of nickel(L) units, alternatively connected by single and double bridging mu-(1,3-N-3) ligand with isolated perchlorate anions. Variable temperature magnetic susceptibility data of the complex have been measured and the fitting,of magnetic data was carried out applying the Borris-Almenar formula for such types of alternating one-dimensional S = 1 systems, based on the Hamiltonian H = -J Sigma(S2iS2i-1 + aS(2i)S(2i+1)). The best-fit parameters obtained are J = -106.7 +/- 2 cm(-1); a = 0.82 +/- 0.02; g = 2.21 +/- 0.02. Complex 2 is a 2D network of 4,4 topology with the nodes occupied by the Cu-II ions, and the edges formed by single azide and double bpds connectors. The perchlorate anions are located between pairs of bpds. The magnetic data have been fitted considering the complex as a pseudo-one-dimensional system, with all copper((II)) atoms linked by [mu(1,3-azido) bridging ligands at axial positions (long Cu...N-3 distances) since the coupling through long bpds is almost nil. The best-fit parameters obtained with this model are J = -1.21 +/- 0.2 cm(-1), g 2.14 +/- 0.02. (c) Wiley-VCH Verlag GmbH & Co. KGaA, 69451 Weinheim, Germany, 2005).

Influence of counter anions on the structures and magnetic properties of trinuclear Cu(II) complexes containing a mu(3)-OH core and Schiff base ligands

Four trinuclear Cu(II) complexes, [(CuL1)(3)(mu(3)-OH)](NO3)(2) (1), [(CuL2)(3)(mu(3)-OH)](I)(2)center dot H2O (2), [(CuL3)(3)(mu(3)-OH)](I)(2) (3) and [(CuL1)(3)(mu(3)-OH)][(CuI3)-I-1] (4), where HL1 (8-amino-4-methyl-5-azaoct-3-en-2-one), HL2 [7-amino-4-methyl-5-azaoct-3-en-2-one] and HL3 [7-amino-4-methyl-5-azahept-3-en-2- one] are the three tridentate Schiff bases, have been synthesized and structurally characterized by X-ray crystallography. All four complexes contain a partial cubane core, [(CuL)(3)(mu(3)-OH)](2+) in which the three [CuL] subunits are interconnected through two types of oxygen bridges afforded by the oxygen atoms of the ligands and the central OH- group. The copper(II) ions are in a distorted square-pyramidal environment. The equatorial plane consists of the bridging oxygen of the central OH- group together with three atoms (N, N, O) from the Schiff base. The oxygen atom of the Schiff base also coordinates to the axial position of Cu(II) of another subunit to form the cyclic trimer. Magnetic susceptibilities have been determined for these complexes over the temperature range of 2-300 K. The isotropic Hamiltonian, H = -J(12)S(1)S(2) - J(13)S(1)S(3) - J(23)S(2)S(3) has been used to interpret the magnetic data. The best fit parameters obtained are: J = - 54.98 cm(-1) g = 2.24 for 1; J = - 56.66 cm(-1), g = 2.19 for 2; J = -44.39 cm(-1), g = 2.16 for 3; J = - 89.92 cm(-1), g = 2.25 for 4. The EPR data at low temperature indicate that the phenomenon of spin frustration occurs for complexes 1-3. (c) 2007 Elsevier B.V. All rights reserved.

Incorporation of a sodium ion guest in the host of copper(II)-Schiff-base complexes: Structural characterization and magnetic study

Three copper(II) complexes, [CuL1], [CuL2] and [CuL3] where L-1, L-2 and L-3 are the tetradentate di-Schiff-base ligands prepared by the condensation of acetylacetone and appropriate diamines (e.g. 1,2-diaminoethane, 1,2-diaminopropane and 1,3-diaminopropane, respectively) in 2:1 ratios, have been prepared. These complexes act as host molecules and include a guest sodium ion by coordinating through the oxygen atoms to result in corresponding new trinuclear complexes, [(CuL1)(2)Na(ClO4)(H2O)][CuL1], [(CuL2)(2)Na(ClO4)(H2O)] (2) and [(CuL3)(2)Na(ClO4)(H2O)] (3) when crystallized from methanol solution containing sodium perchlorate. All three complexes have been characterized by single crystal X-ray crystallography. In all the complexes, the sodium cation has a six-coordinate distorted octahedral environment being bonded to four oxygen atoms from two Schiff-base complexes of Cu(II) in addition to a perchlorate anion and a water molecule. The copper atoms are four coordinate in a square planar environment being bonded to two oxygen atoms and two nitrogen atoms of the Schiff-base ligand. The variable temperature susceptibilities for complexes 1-3 were measured over the range 2-300 K. The isotropic Hamiltonian, H = g(1)beta HS1 + g(2)beta HS2 + J(12)S(1)S(2) + g(3)beta HS3 for complex 1 and H = g(1)beta HS1 + g2 beta HS2 +J(12)S(1)S(2) for complexes 2 and 3 has been used to interpret the magnetic data. The best fit parameters obtained are: g(1) = g(2) = 2.07(0), J = - 1.09(1) cm(-1) for complex 1, g(1) = g(2) = 2.06(0), J = -0.55(1) cm(-1) for complex 2 and g1 = g2 = 2.07(0).J = -0.80(1) cm(-1) for 3. Electrochemical studies displayed an irreversible Cu(II)/Cu(I) one-electron reduction process. (C) 2008 Elsevier Ltd. All rights reserved.

Cu-II acetate complexes involving N,N,O donor Schiff base ligands: Mono-atomic oxygen bridged dimers and alternating chains of the dimers and Cu-2(OAc)(4)

Two tridentate N,N,O donor Schiff bases, HL1 (4-(2-ethylamino-ethylimino)-pentan-2-one) and HL2 (3-(2-amino-propylimino)-1-phenyl-butan-1-one) on reaction with Cu-II acetate in presence of triethyl amine yielded two basal-apical, mono-atomic acetate oxygen-bridging dimeric copper(II) complexes, [Cu2L21(OAc)(2)] (1), [Cu2L22(OAc)(2)] (2). Whereas two other similar tridentate ligands HL3 (4-(2-amino-propylimino)-pentane-2-one) and HL3 (3-(2-amino-ethylimino)-1-phenyl-butan-1-one) under the same conditions produced a mixture of the corresponding dinners and a one-dimensional alternating chain of the dimer and copper acetate moiety, [Cu4L23(OAc)(6)](n) (3) and [Cu4L24(OAc)(6)](n) (4), formed by a very rare mu(3) bridging mode of the acetate ion. All four complexes (1-4) have been characterized by X-ray crystallography. The isotropic Hamiltonian, H = -JS(1)S(2) has been used to interpret the magnetic data. Magnetic measurements of 1 and 2 in the temperature range 2-300 K reveal a very weak antiferromagnetic coupling for both complexes U = -0.56 and -1.19 cm(-1) for 1 and 2, respectively). (C) 2008 Elsevier Ltd. All rights reserved.

The vibrations and tunnelling of malonaldehyde on a Moller-Plesset surface

The vibrations and tunnelling motion of malonaldehyde have been studied in their full dimensionality using an internal coordinate path Hamiltonian. In this representation there is one large amplitude internal coordinate s and 3N - 7 (=20) normal coordinates Q which are orthogonal to the large amplitude motion at all points. It is crucial that a high accuracy potential energy surface is used in order to obtain a good representation for the tunneling motion; we use a Moller-Plesset (MP2) surface. Our methodology is variational, that is we diagonalize a sufficiently large matrix in order to obtain the required vibrational levels, so an exact representation for the kinetic energy operator is used. In a harmonic valley representation (s, Q) complete convergence of the normal coordinate motions and the internal coordinate motions has been obtained; for the anharmonic valley in which we use two- and three-body terms in the surface (s, Q(1), Q(2)), we also obtain complete convergence. Our final computed stretching fundamentals are deficient because our potential energy surface is truncated at quartic terms in the normal coordinates, but our lower fundamentals are good.

Optimal kinematic control of an autonomous underwater vehicle

This note investigates the motion control of an autonomous underwater vehicle (AUV). The AUV is modeled as a nonholonomic system as any lateral motion of a conventional, slender AUV is quickly damped out. The problem is formulated as an optimal kinematic control problem on the Euclidean Group of Motions SE(3), where the cost function to be minimized is equal to the integral of a quadratic function of the velocity components. An application of the Maximum Principle to this optimal control problem yields the appropriate Hamiltonian and the corresponding vector fields give the necessary conditions for optimality. For a special case of the cost function, the necessary conditions for optimality can be characterized more easily and we proceed to investigate its solutions. Finally, it is shown that a particular set of optimal motions trace helical paths. Throughout this note we highlight a particular case where the quadratic cost function is weighted in such a way that it equates to the Lagrangian (kinetic energy) of the AUV. For this case, the regular extremal curves are constrained to equate to the AUV's components of momentum and the resulting vector fields are the d'Alembert-Lagrange equations in Hamiltonian form.

Generalized matching networks and their properties

In this paper, we introduce two kinds of graphs: the generalized matching networks (GMNs) and the recursive generalized matching networks (RGMNs). The former generalize the hypercube-like networks (HLNs), while the latter include the generalized cubes and the star graphs. We prove that a GMN on a family of k-connected building graphs is -connected. We then prove that a GMN on a family of Hamiltonian-connected building graphs having at least three vertices each is Hamiltonian-connected. Our conclusions generalize some previously known results.

Locally twisted cubes are 4-pancyclic

The locally twisted cube is a newly introduced interconnection network for parallel computing. Ring embedding is an important issue for evaluating the performance of an interconnection network. In this paper, we investigate the problem of embedding rings into a locally twisted cube. Our main contribution is to find that, for each integer l is an element of (4,5,...,2(n)}, a ring of length I can be embedded into an n-dimensional locally twisted cube so that both the dilation and the load factor are one. As a result, a locally twisted cube is Hamiltonian. We conclude that a locally twisted cube is superior to a hypercube in terms of ring embedding capability. (C) 2004 Elsevier Ltd. All rights reserved.