Estimating oceanic heat content change using isotherms

This paper presents a new analysis of ocean heat content changes over the last 50 yr using isotherms by calculating the mean temperature above the 148C isotherm and the depth of the 148C isotherm as separate variables. A new quantity called the ‘‘relative heat content’’ (‘‘RHC’’) is introduced, which represents the minimum local heat content change over time, relative to a fixed isotherm. It is shown how mean temperature and isotherm depth changes make separable and additive contributions to changes in RHC. Maps of RHC change between 1970 and 2000 show similar spatial patterns to a traditional fixed-depth ocean heat content change to 220 m. However, the separate contributions to RHC suggest a more spatially uniform contribution from warming above the isotherm, while isotherm depth changes show wind-driven signals, of which some are identifiable as being related to the North Atlantic Oscillation. The time series show that the warming contribution to RHC dominates the global trend, while the depth contribution only dominates on the basin scale in the North Atlantic. The RHC shows minima associated with the major volcanic eruptions (particularly in the Indian Ocean), and these are entirely contributed by mean temperature changes rather than isotherm depth changes. The depth change contributions to RHC are strongly affected by the recently reported XBT fall-rate bias, whereas the mean temperature contributions are not. Therefore, only the isotherm depth change contributions toRHCwill need to be reassessed as fall-rate-corrected data become available.

Analytical potentials for triatomic molecules from spectroscopic data V. Application to HOX (X=F, Cl, Br, I)

Analytical potential energy functions are reported for HOX (X=F, Cl, Br, I). The surface for HOF predicts two metastable minima as well as the equilibrium configuration. These correspond to HFO (bent) and OHF (linear). Ab initio calculations performed for the HOF surface confirm these predictions. Comparisons are drawn between the two sets of results, and a vibrational analysis is undertaken for the hydrogen bonded OHF species. For HOCl, one further minimum is predicted, corresponding to HClO (bent), the parameters for which compare favourably with those reported from ab initio studies. In contrast, only the equilibrium configurations are predicted to be stable for HOBr and HOI.

Analytical potentials for triatomic molecules from spectroscopic data VI. On choosing internal coordinates to describe triatomic potential energy surfaces with more than one minimum

Orthogonal internal coordinates are defined which have useful properties for constructing the potential energy functions of triatomic molecules with two or three minima on the surface. The coordinates are used to obtain ground state potentials of ClOO and HOF, both of which have three minima.

A benchmark method for global optimization problems in structure determination from powder diffraction data

Quasi-Newton-Raphson minimization and conjugate gradient minimization have been used to solve the crystal structures of famotidine form B and capsaicin from X-ray powder diffraction data and characterize the chi(2) agreement surfaces. One million quasi-Newton-Raphson minimizations found the famotidine global minimum with a frequency of ca 1 in 5000 and the capsaicin global minimum with a frequency of ca 1 in 10 000. These results, which are corroborated by conjugate gradient minimization, demonstrate the existence of numerous pathways from some of the highest points on these chi(2) agreement surfaces to the respective global minima, which are passable using only downhill moves. This important observation has significant ramifications for the development of improved structure determination algorithms.

The role of biotic and abiotic factors in evolution of ant dispersal in the milkwort family (Polygalaceae)

A phylogenetic approach was taken to investigate the evolutionary history of seed appendages in the plant family Polygalaceae (Fabales) and determine which factors might be associated with evolution of elaiosomes through comparisons to abiotic (climate) and biotic (ant species number and abundance) timelines. Molecular datasets from three plastid regions representing 160 species were used to reconstruct a phylogenetic tree of the order Fabales, focusing on Polygalaceae. Bayesian dating methods were used to estimate the age of the appearance of ant-dispersed elaiosomes in Polygalaceae, shown by likelihood optimizations to have a single origin in the family. Topology-based tests indicated a diversification rate shift associated with appearance of caruncular elaiosomes. We show that evolution of the caruncular elaiosome type currently associated with ant dispersal occurred 54.0-50.5 million year ago. This is long after an estimated increase in ant lineages in the Late Cretaceous based on molecular studies, but broadly concomitant with increasing global temperatures culminating in the Late Paleocene-Early Eocene thermal maxima. These results suggest that although most major ant clades were present when elaiosomes appeared, the environmental significance of elaiosomes may have been an important factor in success of elaiosome-bearing lineages. Ecological abundance of ants is perhaps more important than lineage numbers in determining significance of ant dispersal. Thus, our observation that elaiosomes predate increased ecological abundance of ants inferred from amber deposits could be indicative of an initial abiotic environmental function.

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.

An investigation of the germylene addition reaction, GeH2+C2H2: Time-resolved gas-phase kinetic studies and quantum chemical calculations of the reaction energy surface

Time resolved studies of germylene, GeH2, generated by laser flash photolysis of 3,4-dimethylgermacyclopentene-3, have been carried out to obtain rate constants for its bimolecular reaction with acetylene, C2H2. The reaction was studied in the gas-phase over the pressure range 1-100 Tort, with SF6 as bath gas, at 5 temperatures in the range 297-553 K. The reaction showed a very slight pressure dependence at higher temperatures. The high pressure rate constants (obtained by extrapolation at the three higher temperatures) gave the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1)) (-10.94 +/- 0.05) + (6.10 +/- 0.36 kJ mol(-1))/RTln10. These Arrhenius parameters are consistent with a fast reaction occurring at approximately 30% of the collision rate at 298 K. Quantum chemical calculations (both DFT and ab initio G2//B3LYP and G2//QCISD) of the GeC2H4 potential energy surface (PES), show that GeH2 + C2H2 react initially to form germirene which can isomerise to vinylgermylene with a relatively low barrier. RRKM modelling, based on a loose association transition state, but assuming vinylgermylene is the end product (used in combination with a weak collisional deactivation model) predicts a strong pressure dependence using the calculated energies, in conflict with the experimental evidence. The detailed GeC2H4 PES shows considerable complexity with ten other accessible stable minima (B3LYP level), the three most stable of which are all germylenes. Routes through this complex surface were examined in detail. The only product combination which appears capable of satisfying the (P-3) + C2H4.C2H4 was confirmed as a product by GC observed lack of a strong pressure dependence is Ge(P-3) + C2H4. C2H4 was confirmed as a product by GC analysis. Although the formation of these products are shown to be possible by singlet-triplet curve crossing during dissociation of 1-germiranylidene (1-germacyclopropylidene), it seems more likely (on thermochernical grounds) that the triplet biradical, (GeCH2CH2.)-Ge-., is the immediate product precursor. Comparisons are made with the reaction of SiH2 with C2H2.

Mapping fibre orientation in complex-shaped biological systems with micrometre resolution by scanning X-ray microdiffraction

A fully automated procedure to extract and to image local fibre orientation in biological tissues from scanning X-ray diffraction is presented. The preferred chitin fibre orientation in the flow sensing system of crickets is determined with high spatial resolution by applying synchrotron radiation based X-ray microbeam diffraction in conjunction with advanced sample sectioning using a UV micro-laser. The data analysis is based on an automated detection of azimuthal diffraction maxima after 2D convolution filtering (smoothing) of the 2D diffraction patterns. Under the assumption of crystallographic fibre symmetry around the morphological fibre axis, the evaluation method allows mapping the three-dimensional orientation of the fibre axes in space. The resulting two-dimensional maps of the local fibre orientations - together with the complex shape of the flow sensing system - may be useful for a better understanding of the mechanical optimization of such tissues.

Isolation, fractionation and characterisation of proteins from Mucuna bean

The chemical composition and fractional distribution of protein isolates prepared from species of Mucuna bean were studied. Using six different extraction media, the yield of protein based on the Kjeldahl procedure varied from 8% to 34%, and the protein content varied from 75% to 95%. When the yields were high, the colour of the isolates generally tended to be dark and unsatisfactory. Hence, the use of chemical treatments and high pressure processing were explored. The solubility maxima for the protein isolates in water were found to occur at pH values of 2.0 and 11.0, while the pH corresponding to minimum solubility (i.e. isoelectric region) occurred at pH values of 4.0 and 5.0. The total essential amino acid in the isolates ranged from 495 to 557 mg g(-1) protein, which compares favourably with the recommended level for pre-school and school children. Methionine and cysteine were the limiting amino acids. A key nutritional attribute of the protein isolates was its high lysine content. The isolate can therefore complement cereal-based foods which are deficient in lysine. The proteins mainly consisted of albumins, glutelins and globulins. Prolamins were only present in trace concentration (< 0.3%). Gel filtration chromatograms of the isolates indicated the presence of major protein fractions with molecular weights of 40 and 15 kDa, while gel electrophoresis (SDS-PAGE) indicated a major broad zone with molecular weights of 36 +/- 7 and 17.3 +/- 13 kDa. (c) 2006 Elsevier Ltd. All rights reserved.

Search space pruning and global optimization of multiple gravity assist trajectories with deep space manoeuvres

This paper deals with the design of optimal multiple gravity assist trajectories with deep space manoeuvres. A pruning method which considers the sequential nature of the problem is presented. The method locates feasible vectors using local optimization and applies a clustering algorithm to find reduced bounding boxes which can be used in a subsequent optimization step. Since multiple local minima remain within the pruned search space, the use of a global optimization method, such as Differential Evolution, is suggested for finding solutions which are likely to be close to the global optimum. Two case studies are presented.

Retrieval of global atmospheric electrical activity at a polluted urban site

The atmospheric electrical Potential Gradient (PG) arises from global thunderstorm activity, but surface measurements of the atmospheric Potential Gradient (PG) are influenced by global thunderstorms and local aerosol concentration changes. The local aerosol change can be monitored independently, and in some cases the concentration changes are closely related to PG changes. For these circumstances, a general theory to remove the local aerosol influence on PG measurements has been developed. Continuous measurements of PG and aerosol mass concentration were made during 24–31 Dec, 2005 within an urban environment at Reading, UK. The average diurnal variation of PG showed a double diurnal cycle, with maxima in the early morning and evening hours. The aerosol concentration has similar double maxima. Removing the aerosol using from the PG and aerosol correlation returns a single diurnal cycle, suggestive of the more global PG diurnal cycle.

Ferric and ferrous iron in nitroso-myoglobin: computer simulations of stable and metastable states and their infrared spectra

The binding of NO to iron is involved in the biological function of many heme proteins. Contrary to ligands like CO and O-2, which only bind to ferrous (Fe-II) iron, NO binds to both ferrous and ferric (Fe-II) iron. In a particular protein, the natural oxidation state can therefore be expected to be tailored to the required function. Herein, we present an ob initio potential-energy surface for ferric iron interacting with NO. This potential-energy surface exhibits three minima corresponding to eta'-NO coordination (the global minimum), eta(1)-ON coordination and eta(2) coordination. This contrasts with the potential-energy surface for Fe-II-NO, which ex- hibits only two minima (the eta(2) coordination mode for Fe-II is a transition state, not a minimum). In addition, the binding energies of NO are substantially larger for Fe-III than for Fe-II. We have performed molecular dynamics simulations for NO bound to ferric myoglobin (Mb(III)) and compare these with results obtained for Mb(II). Over the duration of our simulations (1.5 ns), all three binding modes are found to be stable at 200 K and transiently stable at 300 K, with eventual transformation to the eta(1)-NO global-minimum conformation. We discuss the implication of these results related to studies of rebinding processes in myoglobin.

Potential energy surface and moleculardynamics of MbNO: existence of an unsuspected FeON minimum

Ligands such as CO, O2, or NO are involved in the biological function of myoglobin. Here we investigate the energetics and dynamics of NO interacting with the Fe(II) heme group in native myoglobin using ab initio and molecular dynamics simulations. At the global minimum of the ab initio potential energy surface (PES), the binding energy of 23.4 kcal/mol and the Fe-NO structure compare well with the experimental results. Interestingly, the PES is found to exhibit two minima: There exists a metastable, linear Fe-O-N minimum in addition to the known, bent Fe-N-O global minimum conformation. Moreover, the T-shaped configuration is found to be a saddle point, in contrast to the corresponding minimum for NO interacting with Fe(III). To use the ab initio results for finite temperature molecular dynamics simulations, an analytical function was fitted to represent the Fe-NO interaction. The simulations show that the secondary minimum is dynamically stable up to 250 K and has a lifetime of several hundred picoseconds at 300 K. The difference in the topology of the heme-NO PES from that assumed previously (one deep, single Fe-NO minimum) suggests that it is important to use the full PES for a quantitative understanding of this system. Why the metastable state has not been observed in the many spectroscopic studies of myoglobin interacting with NO is discussed, and possible approaches to finding it are outlined.

A novel heteroditopic terpyridine-pincer ligand as building block for mono- and heterometallic Pd(II) and Ru(II) complexes

A palladium-catalyzed Stille coupling reaction was employed as a versatile method for the synthesis of a novel terpyridine-pincer (3, TPBr) bridging ligand, 4'-{4-BrC6H2(CH2NMe2)(2)-3,5}-2,2':6',2 ''-terpyridine. Mononuclear species [PdX(TP)] (X = Br, Cl), [Ru(TPBr)(tpy)](PF6)(2), and [Ru(TPBr)(2)](PF6)(2), synthesized by selective metalation of the NCNBr-pincer moiety or complexation of the terpyridine of the bifunctional ligand TPBr, were used as building blocks for the preparation of heterodi- and trimetallic complexes [Ru(TPPdCl)(tpy)](PF6)(2) (7) and [Ru(TPPdCl)(2)]-(PF6)(2) (8). The molecular structures in the solid state of [PdBr(TP)] (4a) and [Ru(TPBr)(2)](PF6)(2) (6) have been determined by single-crystal X-ray analysis. Electrochemical behavior and photophysical properties of the mono-and heterometallic complexes are described. All the above di- and trimetallic Ru complexes exhibit absorption bands attributable to (MLCT)-M-1 (Ru -> tpy) transitions. For the heteroleptic complexes, the transitions involving the unsubstituted tpy ligand are at a lower energy than the tpy moiety of the TPBr ligand. The absorption bands observed in the electronic spectra for TPBr and [PdCl(TP)] have been assigned with the aid of TD-DFT calculations. All complexes display weak emission both at room temperature and in a butyronitrile glass at 77 K. The considerable red shift of the emission maxima relative to the signal of the reference compound [Ru(tpy)(2)](2+) indicates stabilization of the luminescent (MLCT)-M-3 state. For the mono- and heterometallic complexes, electrochemical and spectroscopic studies (electronic absorption and emission spectra and luminescence lifetimes recorded at room temperature and 77 K in nitrile solvents), together with the information gained from IR spectroelectrochemical studies of the dimetallic complex [Ru(TPPdSCN)(tpy)](PF6)(2), are indicative of charge redistribution through the bridging ligand TPBr. The results are in line with a weak coupling between the {Ru(tpy)(2)} chromophoric unit and the (non)metalated NCN-pincer moiety.

Was UV spectral solar irradiance lower during the recent low sunspot minimum?

A detailed analysis is presented of solar UV spectral irradiance for the period between May 2003 and August 2005, when data are available from both the Solar Ultraviolet pectral Irradiance Monitor (SUSIM) instrument (on board the pper Atmosphere Research Satellite (UARS) spacecraft) and the Solar Stellar Irradiance Comparison Experiment (SOLSTICE) instrument (on board the Solar Radiation and Climate Experiment (SORCE) satellite). The ultimate aim is to develop a data composite that can be used to accurately determine any differences between the “exceptional” solar minimum at the end of solar cycle 23 and the previous minimum at the end of solar cycle 22 without having to rely on proxy data to set the long‐term change. SUSIM data are studied because they are the only data available in the “SOLSTICE gap” between the end of available UARS SOLSTICE data and the start of the SORCE data. At any one wavelength the two data sets are considered too dissimilar to be combined into a meaningful composite if any one of three correlations does not exceed a threshold of 0.8. This criterion removes all wavelengths except those in a small range between 156 nm and 208 nm, the longer wavelengths of which influence ozone production and heating in the lower stratosphere. Eight different methods are employed to intercalibrate the two data sequences. All methods give smaller changes between the minima than are seen when the data are not adjusted; however, correcting the SUSIM data to allow for an exponentially decaying offset drift gives a composite that is largely consistent with the unadjusted data from the SOLSTICE instruments on both UARS and SORCE and in which the recent minimum is consistently lower in the wave band studied.