In search of simple structures in climate: Simplifying EOFs

Empirical orthogonal functions (EOFs) are widely used in climate research to identify dominant patterns of variability and to reduce the dimensionality of climate data. EOFs, however, can be difficult to interpret. Rotated empirical orthogonal functions (REOFs) have been proposed as more physical entities with simpler patterns than EOFs. This study presents a new approach for finding climate patterns with simple structures that overcomes the problems encountered with rotation. The method achieves simplicity of the patterns by using the main properties of EOFs and REOFs simultaneously. Orthogonal patterns that maximise variance subject to a constraint that induces a form of simplicity are found. The simplified empirical orthogonal function (SEOF) patterns, being more 'local'. are constrained to have zero loadings outside the main centre of action. The method is applied to winter Northern Hemisphere (NH) monthly mean sea level pressure (SLP) reanalyses over the period 1948-2000. The 'simplified' leading patterns of variability are identified and compared to the leading patterns obtained from EOFs and REOFs. Copyright (C) 2005 Royal Meteorological Society.

Vibration-rotation bands of diazomethane

Some absorption bands of diazomethane vapour between 1950-3500 cm-1 have been measured with very high resolving power. The rotational structure of two parallel bands and of one perpendicular band has been resolved, and approximate values have been determined for the rotational constants. The results are consistent with the geometrical structure usually accepted for this molecule. A peculiarity in the results for the band near 2100 cm-1, together with other facts, leads to the suggestion that a tautomeric form of this molecule exists, HCN=NH, being an isoelectronic analogue of hydrazoic acid.

Multi-model assessment of stratospheric ozone return dates and ozone recovery in CCMVal-2 models

Projections of stratospheric ozone from a suite of chemistry-climate models (CCMs) have been analyzed. In addition to a reference simulation where anthropogenic halogenated ozone depleting substances (ODSs) and greenhouse gases (GHGs) vary with time, sensitivity simulations with either ODS or GHG concentrations fixed at 1960 levels were performed to disaggregate the drivers of projected ozone changes. These simulations were also used to assess the two distinct milestones of ozone returning to historical values (ozone return dates) and ozone no longer being influenced by ODSs (full ozone recovery). The date of ozone returning to historical values does not indicate complete recovery from ODSs in most cases, because GHG-induced changes accelerate or decelerate ozone changes in many regions. In the upper stratosphere where CO2-induced stratospheric cooling increases ozone, full ozone recovery is projected to not likely have occurred by 2100 even though ozone returns to its 1980 or even 1960 levels well before (~2025 and 2040, respectively). In contrast, in the tropical lower stratosphere ozone decreases continuously from 1960 to 2100 due to projected increases in tropical upwelling, while by around 2040 it is already very likely that full recovery from the effects of ODSs has occurred, although ODS concentrations are still elevated by this date. In the midlatitude lower stratosphere the evolution differs from that in the tropics, and rather than a steady decrease in ozone, first a decrease in ozone is simulated from 1960 to 2000, which is then followed by a steady increase through the 21st century. Ozone in the midlatitude lower stratosphere returns to 1980 levels by ~2045 in the Northern Hemisphere (NH) and by ~2055 in the Southern Hemisphere (SH), and full ozone recovery is likely reached by 2100 in both hemispheres. Overall, in all regions except the tropical lower stratosphere, full ozone recovery from ODSs occurs significantly later than the return of total column ozone to its 1980 level. The latest return of total column ozone is projected to occur over Antarctica (~2045–2060) whereas it is not likely that full ozone recovery is reached by the end of the 21st century in this region. Arctic total column ozone is projected to return to 1980 levels well before polar stratospheric halogen loading does so (~2025–2030 for total column ozone, cf. 2050–2070 for Cly+60×Bry) and it is likely that full recovery of total column ozone from the effects of ODSs has occurred by ~2035. In contrast to the Antarctic, by 2100 Arctic total column ozone is projected to be above 1960 levels, but not in the fixed GHG simulation, indicating that climate change plays a significant role.

Metal complexes of a tetraazacyclophane: solution and molecular modelling studies

An alternative synthetic approach to yield the compound 2,3,5,6,8,9,11,14-octahydrobenzo[1][ 1,4,7,10]tetraazacyclotetradecine (bz[14]N-4) is presented. The protonation constants of bz[14]N-4 and the stability constants of its complexes with Ni2+, Cu2+, Zn2+, Cd2+, Pb2+ were determined in H2O at 25degreesC with ionic strength 0.10 mol dm(-3) in KNO3 and they were compared with structurally related macrocycles cyclam (1,4,8,11-tetraazacyclotetradecane) and cyclen (1,4,7,10-tetraazacyclododecane). These studies indicate that only 1 : 1 ( M : L) species are formed in solution, and the ligand exhibits a high affinity for larger ions such as Cd2+ and Pb2+. The X-ray study of [bz[14]N4H3](3+) shows that an inclusion compound with a chloride counter-anion is formed through NH...Cl hydrogen bonds. Spectroscopic data in solution ( electronic and NMR spectra) showed that the macrocycle adopts a planar arrangement upon metal complexation. Molecular mechanics calculations reveal that in spite of the presence of the benzene ring in the macrocyclic framework this ligand can encapsulate metal ions with different stereo-electronic sizes in square planar arrangements. Our results indicate that the presence of the benzene ring in the backbone of the bz[14]N-4 confers a coordination behaviour intermediate between that of cyclam and cyclen.

Reactions of silylene with unreactive molecules. 2. nitrogen: gas-phase kinetic and theoretical studies

Time-resolved studies of the reaction of silylene, SiH2, with N-2 have been attempted at 296, 417, and 484 K, using laser flash photolysis to generate and monitor SiH2. No conclusive evidence for reaction could be found even with pressures of N-2 of 500 Torr. This enables us to set upper limits of ca. 3 x 10(-15) cm(3) molecule(-1) s(-1) for the second-order rate constants. A lower limit for the activation energy, E-a, of ca. 47 kJ mol(-1) is also derived. Ab initio calculations at the G3 level indicate that the only SiH2N2 species of lower energy than the separated reactants is the H2Si...N-2 donor-acceptor (ylid) species with a relative enthalpy of -26 kJ mol(-1), insufficient for observation of reaction under the experimental conditions. Ten bound species on the SiH2N2 surface were found and their energies calculated as well as those of the potential dissociation products: HSiN + NH((3)Sigma(-)) and HNSi + NH((3)Sigma(-)). Additionally two of the transition states involving cyclic-SiH2N2 (siladiazirine) were explored. It appears that siladiazirine is neither thermodynamically nor kinetically stable. The findings indicate that Si-N-d bonds (where N-d is double-bonded nitrogen) are not particularly strong. An unexpected cyclic intermediate was found in the isomerization of silaisocyanamide to silacyanamide.

Bis[bis(methoxycarbimido)aminato]copper(II) 1-methylpyrrolidin-2-one disolvate

The title compound, [Cu(C4H8N3O2)(2)]center dot 2C(5)H(9)NO, consists of a neutral copper complex, in which the Cu II centre coordinates to two bis(methoxycarbimido) aminate ligands, solvated by two molecules of 1-methylpyrrolidin-2-one. The complex is planar and centrosymmetric, with the Cu II centre occupying a crystallographic inversion centre and adopting approximately square-planar geometry. N-H center dot center dot center dot O hydrogen-bonding interactions exist between the amine NH groups of the ligands and the O atoms of the 1-methylpyrrolidin-2-one molecules. The associated units pack to form sheets.

Solid state annular tautomerism in a molecule containing two imidazole moieties

The X-ray crystal structure of the 1:2 condensate (1) of hydrazine hydrate and 4-methyl-imidazole-5-carboxaldehyde has been determined. The molecule is centrosymmetric crystallising in the space group Fddd with cell dimensions: a = 10.557(14), b = 17.062(22), c = 24.759(27) angstrom. Fourier map shows that the NH hydrogen atom of each imidazole moiety has equal possibility of occupying any of its two ring N atoms. This poses the possibility of finding three tautomers in 1 in the solid state. Consideration of the H-bonding pattern observed in 1 and related B3LYP/6-311+G(2d, p) calculations show that only two tautomers are present in the solid state. The situation is compared with that in the structure of 4(5)-nitro-5(4)-methoxy-imidazole reported previously by Kubicki.

Acid dissociation of 2,2 '-dipyridylamine in non-aqueous medium when chelated to some Ru(II)N-4 cores

[Ru(2,2'-bipyridine)(2)(Hdpa)](BF4)(2) center dot 2H(2)O (1), [Ru(1,10-phenanthroline)(2)(Hdpa)] (PF6)(2) center dot CH2Cl2 (2) and [Ru(4,4,4',4'-tetramethyl-2,2'- bisoxazoline)(2)(Hdpa)] (PF6)(2) (3) are synthesized where Hdpa is 2,2'-dipyridylamine. The X-ray crystal structures of 1 and 2 have been determined. Hdpa in 1 and 2 is found to bind the metal via the two pyridyl N ends. Comparing the NMR spectra in DMSO-d(6), it is concluded that 3 has a similar structure. The pK(a) values (for the dissociation of the NH proton in Hdpa) of free Hdpa and its complexes are determined in acetonitrile by exploiting molar conductance. These correlate linearly with the chemical shift of the NH proton in the respective entities. (C) 2007 Elsevier B.V. All rights reserved.

Design of a new foldamer turn-linker-turn in acyclic hexapeptides and formation of channels through self-assembly

Single crystal X-ray diffraction studies and solvent dependent NMR titration reveal that the designed pepticles I and 11, Boc-Xx(1)-Aib(2)-Yy(3)-NH(CH2)(2)NH-Yy(3)-Aib(2)-Xx(1)-Boc, where Xx and Yy are lie and Leu in peptide I and Leu and Val in peptide 11, respectively, fold into a turn-linker-turn (T-L-T) conformation both in the solid state and in solution. In the crystalline state the T-L-T foldamers; of peptide I and II self-assemble to form a three-dimensional framework of channels. The insides of the channels are hydrophilic and found to contain solvent CHCl3 hydrogen bonded to exposed C=O of Aib located at the turn regions. (c) 2008 Elsevier B.V. All rights reserved.

Infrared study of intercomponent interactions in a switchable hydrogen-bonded rotaxane

The macrocycle in rotaxane 1 is preferentially hydrogen bonded to the succinamide station in the neutral form, but can be moved to the naphthalimide station by one-electron reduction of the latter. The hydrogen bonding between the amide NH groups of the macrocycle and the C=O groups in the binding stations in the thread was studied with IR spectroscopy in different solvents in both states. In addition, the solvent effect on the vibrational frequencies was analyzed; a correlation with the solvent acceptor number (AN) was observed. The conformational switching upon reduction could be detected by monitoring the hydrogen-bond-induced shifts of the v(CO) frequencies of the C=O groups of the succinamide and the reduced naphthalimide stations. The macrocycle was found to shield the encapsulated station from the solvent: wavenumbers of v(CO) bands of the C=O groups residing inside the macrocycle cavity remain unaffected by the solvent polarity.

Overlapping double turn conformations adopted by tetrapeptides containing non-coded alpha-Amino Isobutyric Acid (AIB) and formation of tape-like structures through supramolecular helix mediated self-assembly

Single crystal X-ray diffraction studies and solvent dependent H-1 NMR titrations reveal that a set of four tetrapeptides with general formula Boc-Xx(1)-Aib(2)-Yy(3)-Zz(4)-OMe, where Xx, Yy and Zz are coded L- amino acids, adopt equivalent conformations that can be described as overlapping double turn conformations stabilized by two 4 -> 1 intramolecular hydrogen bonds between Yy(3)-NH and Boc C=O and Zz(4)-NH and Xx(1)C=O. In the crystalline state, the double turn structures are packed in head-to-tail fashion through intermolecular hydrogen bonds to create supramolecular helical structures. Field emission scanning electron microscopic (FE-SEM) images of the tetrapeptides in the solid state reveal that they can form flat tape-like structures. The results establish that synthetic Aib containing supramolecular helices can form highly ordered self-aggregated amyloid plaque like human amylin.

Design and theoretical study of nickel catalysts for syndiotactic polyolefins

A nickel catalyst was modeled with ligand L-2, [ NH = CH-CH = CH-O](-), which should have potential use as a syndiotactic polyolefin catalyst, and the reaction mechanism was studied by theoretical calculations using the density functional method at the B3LYP/ LANL2MB level. The mechanism involves the formation of the intermediate [(NiLMe)-Me-2](+), in which the metal occuples a T-shaped geometry. - This intermediate has two possible structures with the methyl group trans either to the oxygen or to the nitrogen atom of L-2. The results show that both structures can lead to the desired product via similar reaction paths, A and B. Thus, the polymerization could be considered as taking place either with the alkyl group occupying the position trans to the Ni-O or trans to the Ni-N bond in the catalyst. The polymerization process thus favors the catalysis of syndiotactic polyolefins. The syndiotactic synthesis effects could also be enhanced by variations in the ligand substituents. From energy considerations, we can conclude that it is more favorable for the methyl group to occupy the trans-O position to form a complex than to occupy the trans-N position. From bond length considerations, it is also more favoured for ethene to occupy the trans-O position than to occupy the trans-N position.

New polynuclear Mo-Fe complexes with ferrocenylamidobenzimidazole ligands

The reaction between [Mo(eta(3)-C3H5)(CO)(2)(NCMe)(2)Br] (1) and the ferrocenylamidobenzimidazole ligands FcCO(NH(2)benzim) (L1) and (FcCO)(2)(NHbenzim) (L2) led to a binuclear (2) and a trinuclear (3) Mo-Fe complex, respectively. The single-crystal X-ray structure of [Mo(eta(3)-C3H5)(CO)(2)(L2)Br] [L2 = {[(eta(5)-C5H5)Fe(eta(5)-C5H4CO)](2)(2-NH-benzimidazol-yl)}] shows that L2 is coordinated to the endo Mo(eta(3)-C3H5)(CO)(2) group in a kappa(2)-N,O-bidentate chelating fashion whereas the Mo-II centre displays a pseudooctahedral environment with Br occupying an equatorial position. Complex 2 was formulated as [MO(eta(3)-C3H5)(CO)(2)(L1)Br] on the basis of a combination of spectroscopic data, elemental analysis, conductivity and DFT calculations. L1 acts as a kappa(2)-N,N-bidentate ligand. In both L1 and L2, the HOMOs are mainly localised on iron while the C=O bond(s) contribute to the LUMO(s) and the next highest energy orbitals are Fe-allyl antibonding orbitals. When the ligands bind to Mo(eta(3)-C3H5)(CO)(2)Br, the greatest difference is that Mo becomes the strongest contributor to the HOMO. Electrochemical studies show that, in complex 2, no electronic interaction exists between the two ferrocenyl ligands and that the first electron has been removed from the Mo-II-centred HOMO. (c) Wiley-VCH Verlag GmbH & Co. KGaA.

Characterising the variability and extremes of the stratospheric polar vortices using 2D moment analysis

The mean state, variability and extreme variability of the stratospheric polar vortices, with an emphasis on the Northern Hemisphere vortex, are examined using 2-dimensional moment analysis and Extreme Value Theory (EVT). The use of moments as an analysis to ol gives rise to information about the vortex area, centroid latitude, aspect ratio and kurtosis. The application of EVT to these moment derived quantaties allows the extreme variability of the vortex to be assessed. The data used for this study is ECMWF ERA-40 potential vorticity fields on interpolated isentropic surfaces that range from 450K-1450K. Analyses show that the most extreme vortex variability occurs most commonly in late January and early February, consistent with when most planetary wave driving from the troposphere is observed. Composites around sudden stratospheric warming (SSW) events reveal that the moment diagnostics evolve in statistically different ways between vortex splitting events and vortex displacement events, in contrast to the traditional diagnostics. Histograms of the vortex diagnostics on the 850K (∼10hPa) surface over the 1958-2001 period are fitted with parametric distributions, and show that SSW events comprise the majority of data in the tails of the distributions. The distribution of each diagnostic is computed on various surfaces throughout the depth of the stratosphere, and shows that in general the vortex becomes more circular with higher filamentation at the upper levels. The Northern Hemisphere (NH) and Southern Hemisphere (SH) vortices are also compared through the analysis of their respective vortex diagnostics, and confirm that the SH vortex is less variable and lacks extreme events compared to the NH vortex. Finally extreme value theory is used to statistically mo del the vortex diagnostics and make inferences about the underlying dynamics of the polar vortices.

Algorithmic statistical analysis of electrophysiological data for the investigation of structure-activity relationship in single neurons

We are developing computational tools supporting the detailed analysis of the dependence of neural electrophysiological response on dendritic morphology. We approach this problem by combining simulations of faithful models of neurons (experimental real life morphological data with known models of channel kinetics) with algorithmic extraction of morphological and physiological parameters and statistical analysis. In this paper, we present the novel method for an automatic recognition of spike trains in voltage traces, which eliminates the need for human intervention. This enables classification of waveforms with consistent criteria across all the analyzed traces and so it amounts to reduction of the noise in the data. This method allows for an automatic extraction of relevant physiological parameters necessary for further statistical analysis. In order to illustrate the usefulness of this procedure to analyze voltage traces, we characterized the influence of the somatic current injection level on several electrophysiological parameters in a set of modeled neurons. This application suggests that such an algorithmic processing of physiological data extracts parameters in a suitable form for further investigation of structure-activity relationship in single neurons.