Variational calculations of rovibrational states: a precise high-energy potential surface for HCN [and discussion]

We report the results of variational calculations of the rovibrational energy levels of HCN for J = 0, 1 and 2, where we reproduce all the ca. 100 observed vibrational states for all observed isotopic species, with energies up to 18000 cm$^{-1}$, to about $\pm $1 cm$^{-1}$, and the corresponding rotational constants to about $\pm $0.001 cm$^{-1}$. We use a hamiltonian expressed in internal coordinates r$_{1}$, r$_{2}$ and $\theta $, using the exact expression for the kinetic energy operator T obtained by direct transformation from the cartesian representation. The potential energy V is expressed as a polynomial expansion in the Morse coordinates y$_{i}$ for the bond stretches and the interbond angle $\theta $. The basis functions are built as products of appropriately scaled Morse functions in the bond-stretches and Legendre or associated Legendre polynomials of cos $\theta $ in the angle bend, and we evaluate matrix elements by Gauss quadrature. The hamiltonian matripx is factorized using the full rovibrational symmetry, and the basis is contracted to an optimized form; the dimensions of the final hamiltonian matrix vary from 240 $\times $ 240 to 1000 $\times $ 1000.We believe that our calculation is converged to better than 1 cm$^{-1}$ at 18 000 cm$^{-1}$. Our potential surface is expressed in terms of 31 parameters, about half of which have been refined by least squares to optimize the fit to the experimental data. The advantages and disadvantages and the future potential of calculations of this type are discussed.

Local mode and normal mode interpretations of the CH and CD stretching vibrational manifolds of 1,1-dideuteroethene

The =CH2 AND =CD2 stretching vibrational overtones of H2C=CD2 have been studied up to V= 6 and V= 3, respectively. We report their interpretation in terms of a transition from normal to local modes, involving Fermi resonance with the C=C stretching and CH2 scissoring vibrations. We discuss the alternative representation of the vibrational Hamiltonian matrix in local mode and normal mode basis functions, and conclude that the normal mode basis offers greater flexibility in representing small anharmonic couplings with other modes.

Vibration-rotation spectra of deuterated nitrous acid, DONO

High-resolution Fourier transform infrared spectra have been recorded and analyzed for the ν3, ν4, ν5, and ν6 fundamental bands of trans-DONO, and for the ν4 fundamental of cis-DONO. The spectral resolution was better than 0.01 cm−1, and the bands have been fitted using an asymmetric top Hamiltonian with a standard deviation of around 0.0006 cm−1.

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.