The consequences of diffractively spreading beams in ultrafast THz spectroscopy

There are established methods for calculating optical constants from measurements using a broadband terahertz (THz) source. Applications to ultrafast THz spectroscopy have adopted the key assumption that the THz beam is treated as a normal incidence plane-wave. We show that this assumption results in a frequency-dependent systematic error, which is compounded by distortion of the beam on introduction of the sample.

Parallel performance and scalability experiments with the Danish Eulerian Model on the EPCC supercomputers

The Danish Eulerian Model (DEM) is a powerful air pollution model, designed to calculate the concentrations of various dangerous species over a large geographical region (e.g. Europe). It takes into account the main physical and chemical processes between these species, the actual meteorological conditions, emissions, etc.. This is a huge computational task and requires significant resources of storage and CPU time. Parallel computing is essential for the efficient practical use of the model. Some efficient parallel versions of the model were created over the past several years. A suitable parallel version of DEM by using the Message Passing Interface library (AIPI) was implemented on two powerful supercomputers of the EPCC - Edinburgh, available via the HPC-Europa programme for transnational access to research infrastructures in EC: a Sun Fire E15K and an IBM HPCx cluster. Although the implementation is in principal, the same for both supercomputers, few modifications had to be done for successful porting of the code on the IBM HPCx cluster. Performance analysis and parallel optimization was done next. Results from bench marking experiments will be presented in this paper. Another set of experiments was carried out in order to investigate the sensitivity of the model to variation of some chemical rate constants in the chemical submodel. Certain modifications of the code were necessary to be done in accordance with this task. The obtained results will be used for further sensitivity analysis Studies by using Monte Carlo simulation.

Global asymptotic stability in a class of Putnam-type equations

In this paper, we initiate the study of a class of Putnam-type equation of the form x(n-1) = A(1)x(n) + A(2)x(n-1) + A(3)x(n-2)x(n-3) + A(4)/B(1)x(n)x(n-1) + B(2)x(n-2) + B(3)x(n-3) + B-4 n = 0, 1, 2,..., where A(1), A(2), A(3), A(4), B-1, B-2, B-3, B-4 are positive constants with A(1) + A(2) + A(3) + A(4) = B-1 + B-2 + B-3 + B-4, x(-3), x(-2), x(-1), x(0) are positive numbers. A sufficient condition is given for the global asymptotic stability of the equilibrium point c = 1 of such equations. (c) 2005 Elsevier Ltd. All rights reserved.

Two-point boundary value problems for linear evolution equations

We study boundary value problems for a linear evolution equation with spatial derivatives of arbitrary order, on the domain 0 < x < L, 0 < t < T, with L and T positive nite constants. We present a general method for identifying well-posed problems, as well as for constructing an explicit representation of the solution of such problems. This representation has explicit x and t dependence, and it consists of an integral in the k-complex plane and of a discrete sum. As illustrative examples we solve some two-point boundary value problems for the equations iqt + qxx = 0 and qt + qxxx = 0.

Hydrogen bonded supramolecular elastomers : correlating hydrogen bonding strength with morphology and rheology

A series of six low molecular weight elastomers with hydrogen bonding end-groups have been designed, synthesised and studied. The poly(urethane) based elastomers all contained essentially the same hard block content (ca. 11%) and differ only in the nature of their end-groups. Solution state 1H NMR spectroscopic analysis of model compounds featuring the end-groups demonstrate that they all exhibit very low binding constants, in the range 1.4 to 45.0 M-1 in CDCl3, yet the corresponding elastomers each possess a markedly different nanoscale morphology and rheology in the bulk. We are able to correlate small variations of the binding constant of the end-groups with dramatic changes in the bulk properties of the elastomers. These results provide an important insight into the way in which weak non-covalent interactions can be utilized to afford a range of self-assembled polyurethane based materials that feature different morphologies.

Induced-fit binding of π-electron donor substrates to macrocyclic aromatic ether-imide-sulfones: a versatile approach to molecular assembly

Novel macrocyclic receptors which bind electron-donor aromatic substrates via π-stacking donor- acceptor interactions are obtained by cyclo-imidization of an amine-functionalized arylether-sulfone with pyromellitic- and 1,4,5,8-naphthalene-tetracarboxylic dianhydrides. These macrocycles complex with a wide variety of π-donor substrates including tetrathiafulvalene, naphthalene, anthracene, pyrene, perylene, and functional derivatives of these polycyclic hydrocarbons. The resulting supramolecular assemblies range from simple 1:1 complexes, to [2]- and [3]-pseudorotaxanes, and even (as a result of crystallographic disorder) an apparent polyrotaxane. Direct, five-component self-assembly of a metal-centred [3]pseudorotaxane is also observed, on complexation of a macrocyclic ether-imide with 8-hydroxyquinoline in the presence of palladium(II) ions. Binding studies in solution were carried out by 1H NMR and UV-visible spectroscopy, and the stoichiometries of binding were confirmed by Job plots based on charge-transfer absorption bands. The highest association constants are found for strong π-donor guests with large surface-areas, notably perylene and 1-hydroxypyrene, for which Ka values of 1.4 x 103 and 2.3 x 103 M-1 respectively are found. Single crystal X-ray analyses of the receptors and their derived complexes reveal large, induced-fit distortions of the macrocyclic frameworks as a result of complexation. These structures provide compelling evidence for the existence of strong, attractive forces between the electronically-complementary aromatic π-systems of host and guest.

Sterically controlled recognition of macromolecular sequence information by molecular tweezers

Sequence-specific binding is demonstrated between pyrene-based tweezer molecules and soluble, high molar mass copolyimides. The binding involves complementary pi - pi stacking interactions, polymer chain-folding, and hydrogen bonding and is extremely sensitive to the steric environment around the pyromellitimide binding-site. A detailed picture of the intermolecular interactions involved has been obtained through single-crystal X-ray studies of tweezer complexes with model diimides. Ring-current magnetic shielding of polyimide protons by the pyrene '' arms '' of the tweezer molecule induces large complexation shifts of the corresponding H-1 NMR resonances, enabling specific triplet sequences to be identified by their complexation shifts. Extended comonomer sequences (triplets of triplets in which the monomer residues differ only by the presence or absence of a methyl group) can be '' read '' by a mechanism which involves multiple binding of tweezer molecules to adjacent diimide residues within the copolymer chain. The adjacent-binding model for sequence recognition has been validated by two conceptually different sets of tweezer binding experiments. One approach compares sequence-recognition events for copolyimides having either restricted or unrestricted triple-triplet sequences, and the other makes use of copolymers containing both strongly binding and completely nonbinding diimide residues. In all cases the nature and relative proportions of triple-triplet sequences predicted by the adjacent-binding model are fully consistent with the observed H-1 NMR data.

The interplay between geometry, electronic structure, and reactivity of Cu-Ni bimetallic (111) surfaces

The mutual influence of surface geometry (e.g. lattice parameters, morphology) and electronic structure is discussed for Cu-Ni bimetallic (111) surfaces. It is found that on flat surfaces the electronic d-states of the adlayer experience very little influence from the substrate electronic structure which is due to their large separation in binding energies and the close match of Cu and Ni lattice constants. Using carbon monoxide and benzene as probe molecules, it is found that in most cases the reactivity of Cu or Ni adlayers is very similar to the corresponding (111) single crystal surfaces. Exceptions are the adsorption of CO on submonolayers of Cu on Ni(111) and the dissociation of benzene on Ni/Cu(111) which is very different from Ni(111). These differences are related to geometric factors influencing the adsorption on these surfaces.

Thermodynamic analysis of snowball Earth hysteresis experiment: Efficiency, entropy production and irreversibility

We present an extensive thermodynamic analysis of a hysteresis experiment performed on a simplified yet Earth-like climate model. We slowly vary the solar constant by 20% around the present value and detect that for a large range of values of the solar constant the realization of snowball or of regular climate conditions depends on the history of the system. Using recent results on the global climate thermodynamics, we show that the two regimes feature radically different properties. The efficiency of the climate machine monotonically increases with decreasing solar constant in present climate conditions, whereas the opposite takes place in snowball conditions. Instead, entropy production is monotonically increasing with the solar constant in both branches of climate conditions, and its value is about four times larger in the warm branch than in the corresponding cold state. Finally, the degree of irreversibility of the system, measured as the fraction of excess entropy production due to irreversible heat transport processes, is much higher in the warm climate conditions, with an explosive growth in the upper range of the considered values of solar constants. Whereas in the cold climate regime a dominating role is played by changes in the meridional albedo contrast, in the warm climate regime changes in the intensity of latent heat fluxes are crucial for determining the observed properties. This substantiates the importance of addressing correctly the variations of the hydrological cycle in a changing climate. An interpretation of the climate transitions at the tipping points based upon macro-scale thermodynamic properties is also proposed. Our results support the adoption of a new generation of diagnostic tools based on the second law of thermodynamics for auditing climate models and outline a set of parametrizations to be used in conceptual and intermediate-complexity models or for the reconstruction of the past climate conditions. Copyright © 2010 Royal Meteorological Society

Studying reactive processes with classical dynamics: rebinding dynamics in MbNO

A new surface-crossing algorithm suitable for describing bond-breaking and bond-forming processes in molecular dynamics simulations is presented. The method is formulated for two intersecting potential energy manifolds which dissociate to different adiabatic states. During simulations, crossings are detected by monitoring an energy criterion. If fulfilled, the two manifolds are mixed over a finite number of time steps, after which the system is propagated on the second adiabat and the crossing is carried out with probability one. The algorithm is extensively tested (almost 0.5 mu s of total simulation time) for the rebinding of NO to myoglobin. The unbound surface ((FeNO)-N-...) is represented using a standard force field, whereas the bound surface (Fe-NO) is described by an ab initio potential energy surface. The rebinding is found to be nonexponential in time, in agreement with experimental studies, and can be described using two time constants. Depending on the asymptotic energy separation between the manifolds, the short rebinding timescale is between 1 and 9 ps, whereas the longer timescale is about an order of magnitude larger. NO molecules which do not rebind within 1 ns are typically found in the Xenon-4 pocket, indicating the high affinity of NO to this region in the protein.

Kinetic studies of reactions of the nitrate radical (NO3) with peroxy radicals (RO2): an indirect source of OH at night?

A discharge-flow system, coupled to cavity-enhanced absorption spectroscopy (CEAS) detection systems for NO3 at lambda = 662 nm and NO2 at lambda = 404 nm, was used to investigate the kinetics of the reactions of NO3 with eight peroxy radicals at P similar to 5 Torr and T similar to 295 K. Values of the rate constants obtained were (k/10(-12) cm(3) molecule(-1) s(-1)): CH3O2 (1.1 +/- 0.5), C2H5O2 (2.3 +/- 0.7), CH2FO2 (1.4 +/- 0.9), CH2ClO2 (3.8(-2.6)(+1.4)), c-C5H9O2 (1.2(-0.5)(+1.1)), c-C6H11O2 (1.9 +/- 0.7), CF3O2 (0.62 +/- 0.17) and CF3CFO2CF3 (0.24 +/- 0.13). We explore possible relationships between k and the orbital energies of the reactants. We also provide a brief discussion of the potential impact of the reactions of NO3 with RO2 on the chemistry of the night-time atmosphere.

Diorganoruthenium complexes incorporating noninnocent [C6H2(CH2ER2)2-3,5]22-(E = N, P) Bis-Pincer bridging ligands: synthesis, spectroelectrochemistry, and DFT studies

The dinuclear complex [(tpy)Ru-II(PCP-PCP)Ru-II(tPY)]Cl-2 (bridging PCP-PCP = 3,3',5,5'-tetrakis(diphenylphosphinomethyl)biphenyl, [C6H2(CH2PPh2)(2)-3,5](2)(2-)) was prepared via a transcyclometalation reaction of the bis-pincer ligand [PC(H)P-PC(H)P] and the Ru(II) precursor [Ru(NCN)(tpy)]Cl (NCN = [C6H3(CH2NMe2)(2)-2,6](-)) followed by a reaction with 2,2':6',2 ''-terpyridine (tpy). Electrochemical and spectroscopic properties of [(tpy)Ru-II(PCP-PCP)Ru-II(tPY)]Cl-2 are compared with those of the closely related [(tpy)Ru-II(NCN-NCN)Ru-II(tpy)](PF6)(2) (NCN-NCN = [C6H2(CH2- NMe2)(2)-3,5](2)(2-)) obtained by two-electron reduction of [(tpy)Ru-III(NCN-NCN)Ru-III(tpy)](PF6)(4). The molecular structure of the latter complex has been determined by single-crystal X-ray structure determination. One-electron reduction of [(tpy)Ru-III(NCN-NCN)Ru-III(tpy)](PF6)(4) and one-electron oxidation of [(tpy)Ru-II(PCP-PCP)RUII(tpy)]Cl-2 yielded the mixed-valence species [(tpy)Ru-III(NCN-NCN)RUII(tpy)](3+) and [(tpy)Ru-III(PCP-PCP)RUII(tpy)](3+), respectively. The comproportionation equilibrium constants K-c (900 and 748 for [(tpy)Ru-III(NCN-NCN)Ru-III(tpy)](4+) and [(tpy)Ru-II(PCP-PCP)RUII(tpy)](2+), respectively) determined from cyclic voltammetric data reveal comparable stability of the [Ru-III-Ru-II] state of both complexes. Spectroelectrochemical measurements and near-infrared (NIR) spectroscopy were employed to further characterize the different redox states with special focus on the mixed-valence species and their NIR bands. Analysis of these bands in the framework of Hush theory indicates that the mixed-valence complexes [(tpy)Ru-III(PCP-PCP)RUII(tpy)](3+) and [(tpy)Ru-III(NCN-NCN)RUII(tpy)](3+) belong to strongly coupled borderline Class II/Class III and intrinsically coupled Class III systems, respectively. Preliminary DFT calculations suggest that extensive delocalization of the spin density over the metal centers and the bridging ligand exists. TD-DFT calculations then suggested a substantial MLCT character of the NIR electronic transitions. The results obtained in this study point to a decreased metal-metal electronic interaction accommodated by the double-cyclometalated bis-pincer bridge when strong sigma-donor NMe2 groups are replaced by weak sigma-donor, pi-acceptor PPh2 groups

Kinetic studies of reactions of organosilylenes: what have they taught us?

Rate constants for bimolecular reactions, obtained through time-resolved kinetic studies both in the gas and liquid phases are reviewed. Data for reactions of MeSiH, PhSiH, ClSiH, SiCl2, SiMe2, MeSiPh, SiPh2 and SiMes(2) are covered. Where possible, substituent effects relative to SiH2 have been obtained. These demonstrate widely varying effects between different types of reaction, which aids mechanistic understanding. Reactivities are high for all silylenes, but substituents can reduce them by both electronic and steric effects. The gas and liquid phase data (mainly for SiMe2) are compared and appear to be reasonably consistent. This review, although detailed, is not comprehensive.

Gas-phase kinetics of chlorosilylene reactions II. ClSiH + C2H4: absolute rate measurements and quantum chemical and RRKM calculations for the prototype pi addition reaction

Time-resolved studies of chlorosilylene, ClSiH, generated by the 193 nm laser flash photolysis of 1-chloro-1-silacyclopent-3-ene, are carried out to obtain rate constants for its bimolecular reaction with ethene, C2H4, in the gas-phase. The reaction is studied over the pressure range 0.13-13.3 kPa (with added SF6) at five temperatures in the range 296-562 K. The second order rate constants, obtained by extrapolation to the high pressure limits at each temperature, fitted the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1))=(-10.55 +/- 0.10) + (3.86 +/- 0.70) kJ mol(-1)/RT ln10. The Arrhenius parameters correspond to a loose transition state and the rate constant at room temperature is 43% of that for SiH2 + C2H4, showing that the deactivating effect of Cl-for-H substitution in the silylene is not large. Quantum chemical calculations of the potential energy surface for this reaction at the G3MP2//B3LYP level show that, as well as 1-chlorosilirane, ethylchlorosilylene is a viable product. The calculations reveal how the added effect of the Cl atom on the divalent state stabilisation of ClSiH influences the course of this reaction. RRKM calculations of the reaction pressure dependence suggest that ethylchlorosilylene should be the main product. The results are compared and contrasted with those of SiH2 and SiCl2 with C2H4.

Time-resolved gas-phase kinetic, quantum chemical and RRKM studies of reactions of silylene with cyclic ethers

Time-resolved kinetic studies of silylene, SiH2, generated by laser flash photolysis of phenylsilane, have been carried out to obtain rate constants for its bimolecular reactions with oxirane, oxetane, and tetrahydrofuran (THF). The reactions were studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas, at four or five temperatures in the range 294-605 K. All three reactions showed pressure dependences characteristic of third-body-assisted association reactions with, surprisingly, SiH2 + oxirane showing the least and SiH2 + THF showing the most pressure dependence. The second-order rate constants obtained by extrapolation to the high-pressure limits at each temperature fitted the Arrhenius equations where the error limits are single standard deviations: log(k(oxirane)(infinity)/cm(3) molecule(-1) s(-1)) = (-11.03 +/- 0.07) + (5.70 +/- 0.51) kJ mol(-1)/RT In 10 log(k(oxetane)(infinity)/cm(3) molecule(-1) s(-1)) = (-11.17 +/- 0.11) + (9.04 +/- 0.78) kJ mol(-1)/RT In 10 log(k(THF)(infinity)/cm(3) molecule(-1) s(-1)) = (-10.59 +/- 0.10) + (5.76 +/- 0.65) kJ mol(-1)/RT In 10 Binding-energy values of 77, 97, and 92 kJ mol(-1) have been obtained for the donor-acceptor complexes of SiH2 with oxirane, oxetane, and THF, respectively, by means of quantum chemical (ab initio) calculations carried Out at the G3 level. The use of these values to model the pressure dependences of these reactions, via RRKM theory, provided a good fit only in the case of SiH2 + THF. The lack of fit in the other two cases is attributed to further reaction pathways for the association complexes of SiH2 with oxirane and oxetane. The finding of ethene as a product of the SiH2 + oxirane reaction supports a pathway leading to H2Si=O + C2H4 predicted by the theoretical calculations of Apeloig and Sklenak.