Mesoscale simulations of organized convection: importance of convective equilibrium

The validity of convective parametrization breaks down at the resolution of mesoscale models, and the success of parametrized versus explicit treatments of convection is likely to depend on the large-scale environment. In this paper we examine the hypothesis that a key feature determining the sensitivity to the environment is whether the forcing of convection is sufficiently homogeneous and slowly varying that the convection can be considered to be in equilibrium. Two case studies of mesoscale convective systems over the UK, one where equilibrium conditions are expected and one where equilibrium is unlikely, are simulated using a mesoscale forecasting model. The time evolution of area-average convective available potential energy and the time evolution and magnitude of the timescale of convective adjustment are consistent with the hypothesis of equilibrium for case 1 and non-equilibrium for case 2. For each case, three experiments are performed with different partitionings between parametrized and explicit convection: fully parametrized convection, fully explicit convection and a simulation with significant amounts of both. In the equilibrium case, bulk properties of the convection such as area-integrated rain rates are insensitive to the treatment of convection. However, the detailed structure of the precipitation field changes; the simulation with parametrized convection behaves well and produces a smooth field that follows the forcing region, and the simulation with explicit convection has a small number of localized intense regions of precipitation that track with the mid-levelflow. For the non-equilibrium case, bulk properties of the convection such as area-integrated rain rates are sensitive to the treatment of convection. The simulation with explicit convection behaves similarly to the equilibrium case with a few localized precipitation regions. In contrast, the cumulus parametrization fails dramatically and develops intense propagating bows of precipitation that were not observed. The simulations with both parametrized and explicit convection follow the pattern seen in the other experiments, with a transition over the duration of the run from parametrized to explicit precipitation. The impact of convection on the large-scaleflow, as measured by upper-level wind and potential-vorticity perturbations, is very sensitive to the partitioning of convection for both cases. © Royal Meteorological Society, 2006. Contributions by P. A. Clark and M. E. B. Gray are Crown Copyright.

Oxygen precipitation in Czochralski grown silicon heat treated at 525 degrees C

The diffusion of interstitial oxygen In silicon at 525 degrees C is studied using time-of-flight small-angle neutron scattering (SANS) to separate the elastic scattering from oxygen-containing aggregates from the inelastic scattering from neutron-phonon interactions. The growth of oxygen-containing aggregates as a function of time gives a diffusion coefficient, D, calculated from Ham's theory, that is I factor of similar to 3.8 +/- 1.4 times higher than that expected by extrapolation of higher and lower temperature data (D = 0.13 exp(-2.53 eV kT(-1)) cm(2) s(-1)). This result confirms previous observations of enhanced diffusion at intermediate temperatures (400 degrees C-650 degrees C) although the magnitude of the enhancement we find is Much smaller than that reported by some others.

The equilibrium structure of HCN

The equilibrium structure of HCN has been determined from the previously published ground state rotational constants of eight isotopomers by using (B0‐Be) values obtained from a variational calculation of the vibration–rotation spectrum. The results are re(CH)=1.065 01(8) Å, and re(CN)=1.153 24(2) Å.

A potential energy surface for the ground state of formaldehyde, H2CO(1A1)

A model potential energy function for the ground state of H2CO has been derived which covers the whole space of the six internal coordinates. This potential reproduces the experimental energy, geometry and quadratic force field of formaldehyde, and dissociates correctly to all possible atom, diatom and triatom fragments. Thus there are good reasons for believing it to be close to the true potential energy surface except in regions where both hydrogen atoms are close to the oxygen. It leads to the prediction that there should be a metastable singlet hydroxycarbene HCOH which has a planar trans structure and an energy of 2•31 eV above that of equilibrium formaldehyde. The reaction path for dissociation into H2 + CO is predicted to pass through a low symmetry transition state with an activation energy of 4•8 eV. Both of these predictions are in good agreement with recently published ab initio calculations.

The infrared spectrum, equilibrium structure and harmonic and anharmonic force field of thioborine, HBS

Infrared spectra of the two stretching fundamentals of both HBS and DBS have been observed, using a continuous flow system through a multiple reflection long path cell at a pressure around 1 Torr and a Nicolet Fourier Transform spectrometer with a resolution of about 0•1 cm-1. The v3 BS stretching fundamental of DBS, near 1140 cm-1, is observed in strong Fermi resonance with the overtone of the bend 2v2. The bending fundamental v2 has not been observed and must be a very weak band. The analysis of the results in conjunction with earlier work gives the equilibrium structure (re(BH) = 1•1698(12) , re(BS) = 1•5978(3) ) and the harmonic and anharmonic force field.

The equilibrium structure of monofluoroacetylene

The equilibrium rotational constants Be of HCCF and DCCF have been determined from the ground state rotational constants B0, by determining the αr constants for all five fundamentals from the high-resolution vibrational—rotation spectrum making appropriate corrections for the effects of Fermi resonance. By combination with results from the 13C isotopomers and the recent ab initio calculations by Botschwina (Chem. Phys. Lett., 209 (1993) 117), the equilibrium structure is deduced to be: re(CH) = 1.0555(15) Å, re(CC) = 1.1955(8) Å and re(CF) = 1.2781(8) Å.

Inclusion of nonlinear demand-supply relationships within large-scale partial equilibrium linear programming models

This paper presents a new method for the inclusion of nonlinear demand and supply relationships within a linear programming model. An existing method for this purpose is described first and its shortcomings are pointed out before showing how the new approach overcomes those difficulties and how it provides a more accurate and 'smooth' (rather than a kinked) approximation of the nonlinear functions as well as dealing with equilibrium under perfect competition instead of handling just the monopolistic situation. The workings of the proposed method are illustrated by extending a previously available sectoral model for the UK agriculture.

Study of the lactation curve in dairy cattle on farms in Central Mexico

Accurate knowledge of lactation curves has an important relevance to management and research of dairy production systems. A number of equations have been proposed to describe the lactation curve, the most widely applied being the gamma equation. The objective of this work was to compare and evaluate candidate functions for their predictive ability in describing lactation curves from central Mexican dairy cows reared under 2 contrasting management systems. Five equations were considered: Gaines ( exponential decay), Wood ( gamma equation), Rook ( Michaelis-Menten x exponential), and 2 more mechanistic ones (Dijkstra and Pollott). A database consisting of 701 and 1283 records of cows in small-scale and intensive systems, respectively, was used in the analysis. Before analysis, the database was divided into 6 groups representing first, second, and third and higher parity cows in both systems. In all cases except second and above parity cows in small-scale systems, all models improved on the Gaines equation. The Wood equation explained much of the variation, but its parameters do not have direct biological interpretation. Although the Rook equation fitted the data well, some of the parameter estimates were not significant. The Dijkstra equation consistently gave better predictions, and its parameters were usually statistically significant and lend themselves to physiological interpretation. As such, the differences between systems and parity could be explained due to variations in theoretical initial milk production at parturition, specific rates of secretory cell proliferation and death, and rate of decay, all of which are parameters in the model. The Pollott equation, although containing the most biology, was found to be over-parameterized and resulted in nonsignificant parameter estimates. For central Mexican dairy cows, the Dijkstra equation was the best option to use in describing the lactation curve.

Structural and equilibrium studies on mixed ligand complexes of Co(II), Ni(II), Cu(II) and Zn(II) with N-(2-benzimidazolyl)methyliminodiacetic acid and typical N,N donor ligands

Mixed ligand complexes: [Co(L)(bipy)] (.) 3H(2)O (1), [Ni(L)(phen)] (.) H2O (2), [Cu(L)(phen)] (.) 3H(2)O (3) and [Zn(L)(bipy)] (.) 3H(2)O (4), where L2- = two -COOH deprotonated dianion of N-(2-benzimidazolyl)methyliminodiacetic acid (H(2)bzimida, hereafter, H,L), bipy = 2,2' bipyridine and phen = 1,10-phenanthroline have been isolated and characterized by elemental analysis, spectral and magnetic measurements and thermal studies. Single crystal X-ray diffraction studies show octahedral geometry for 1, 2 and 4 and square pyramidal geometry for 3. Equilibrium studies in aqueous solution (ionic strength I = 10(-1) mol dm(-3) (NaNO3), at 25 +/- 1 degrees C) using different molar proportions of M(II):H2L:B, where M = Co, Ni, Cu and Zn and B = phen, bipy and en (ethylene diamine), however, provides evidence of formation of mononuclear and binuclear binary and mixed ligand complexes: M(L), M(H-1L)(-), M(B)(2+), M(L)(B), M(H-1L)(B)(-), M-2(H-1L)(OH), (B)M(H-1L)M(B)(+), where H-1L3- represents two -COOH and the benzimidazole NI-H deprotonated quadridentate (O-, N, O-, N), or, quinquedentate (O-, N, O-, N, N-) function of the coordinated ligand H,L. Binuclear mixed ligand complex formation equilibria: M(L)(B) + M(B)(2+) = (B)M(H-1L)M(B)(+) + H+ is favoured with higher pi-acidity of the B ligands. For Co(II), Ni(II) and Cu(II), these equilibria are accompanied by blue shift of the electronic absorption maxima of M(II) ions, as a negatively charged bridging benzimidazolate moiety provides stronger ligand field than a neutral one. Solution stability of the mixed ligand complexes are in the expected order: Co(II) < Ni(II) < Cu(II) > Zn(II). The Delta logK(M) values are less negetive than their statistical values, indicating favoured formation of the mixed ligand complexes over the binary ones. (c) 2005 Elsevier B.V. All rights reserved.

Equilibrium studies on mixed ligand complex formation of Co(II), Ni(II), Cu(II) and Zn(II) with N-(2-hydroxybenzyl)-L-histidine (H(2)hb-L-his) and typical N,N donor ligands: Crystal structure of [Ni(hb-L-his)(bipyridine)]center dot H2O complex

Equilibrium study on complex formation of Co(II), Ni(II), Cu(II) and Zn(II), hereafter M(II), with the quadridentate (O-, N, O-, N) donor ligand, N-(2-hydroxybenzyl)-L-histidine (H(2)hb-L-his, hereafter H2L), in the absence and in the presence of typical (N, N) donor bidentate ligands, 1,10 phenanthroline(phen), 2, 2'-bipyridine(bipy), ethylenediamine(en), hereafter B, in aqueous solution at 25 +/- 1 degrees C was done at a fixed ionic strength, I = 0.1 mol dm(-3) (NaNO3) by combined pH-metric, UV-Vis and EPR measurements provide evidence for the formation of mononuclear and dinuclear binary and mixed ligand complexes of the types: M(L), M(L)(2)(2-), M-2(L)(2+), M-2(H-1L)(+), M(L)(B), (B)M(H-1L)M(B)(+). The imidazole moiety of the ligand is found to act as a bridging bidentate ligand in the dinuclear M-2(L)(2+), M-2(H-1L)(+) and (B)M(H-1L)M(B)(+) complexes, using its N-3 atom and N1-H deprotonated moiety. Stability constants of the complexes provide evidence of discrimination of Cu(II) from the other M(II) ions by this ligand. Solid complexes: [Ni(L)(H2O)(2)] (1), [Cu(L)(H2O)] (2), and [Ni(L)(bipy)] (.) H2O (3) have been isolated and characterized by various physicochemical studies. Single crystal X-ray diffraction of the ternary complex, 3, shows an octahedral [(O-,N,N,O-)(N,N)] geometry with extensive pi-pi stacking of the aromatic rings and H-bonding with imidazole (N1-H), secondary amino N-atom, the lattice H2O molecule, and the carboxylate and phenolate O-atoms. (c) 2006 Elsevier B.V. All rights reserved.

Time-resolved gas-phase kinetic and quantum chemical studies of the reaction of silylene with oxygen

Time-resolved kinetic studies of the reaction of silylene, SiH2, generated by laser flash photolysis of phenylsilane, have been carried out to obtain rate constants for its bimolecular reaction with O-2. The reaction was studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas, at five temperatures in the range 297-600 K. The second order rate constants at 10 Torr were fitted to the Arrhenius equation: log(k/cm(3) molecule(-1) s(-1)) = (-11.08 +/- 0.04) + (1.57 +/- 0.32 kJ mol(-1))/RT ln10 The decrease in rate constant values with increasing temperature, although systematic is very small. The rate constants showed slight increases in value with pressure at each temperature, but this was scarcely beyond experimental uncertainty. From estimates of Lennard-Jones collision rates, this reaction is occurring at ca. 1 in 20 collisions, almost independent of pressure and temperature. Ab initio calculations at the G3 level backed further by multi-configurational (MC) SCF calculations, augmented by second order perturbation theory (MRMP2), support a mechanism in which the initial adduct, H2SiOO, formed in the triplet state (T), undergoes intersystem crossing to the more stable singlet state (S) prior to further low energy isomerisation processes leading, via a sequence of steps, ultimately to dissociation products of which the lowest energy pair are H2O + SiO. The decomposition of the intermediate cyclo-siladioxirane, via O-O bond fission, plays an important role in the overall process. The bottleneck for the overall process appears to be the T -> S process in H2SiOO. This process has a small spin orbit coupling matrix element, consistent with an estimate of its rate constant of 1 x 10(9) s(-1) obtained with the aid of RRKM theory. This interpretation preserves the idea that, as in its reactions in general, SiH2 initially reacts at the encounter rate with O-2. The low values for the secondary reaction barriers on the potential energy surface account for the lack of an observed pressure dependence. Some comparisons are drawn with the reactions of CH2 + O-2 and SiCl2 + O-2.

Observation of carbonyl oxygen-pi interaction in a metal complex

X-ray crystallography shows that [Ag2L2(H2O)(2)](BF4)(2) where L is a 1:1 condensate of 1,2-diphenylethane-1,2-dione and 2-(2-aminoethyl pyridine), contains an Ag(I)-Ag(I) bond of length 2.979(2) angstrom and an angular, intraligand interaction of the keto O with the pi cloud of the pyridine moiety (O-pyridine centroid = 3.12 angstrom). Model MP2/6-311++G(d,p) calculations indicate that the observed lone pair-pi type interaction is stabilising and not merely a tolerated short contact. (C) 2008 Elsevier B.V. All rights reserved.

Synthesis of chiral 3,4,5,6-tetrahydro-1,4-thiazin-2-ones (thiamorpholin-2-ones)-novel heterocycles possessing enhanced carbonyl electrophilicity over their oxygen counterparts

We report herein the first synthesis of chiral derivatives possessing the 1,4-thiazinone core. As predicted, the thiolactone is more susceptible to nucleophilic attack than the equivalent lactone system.