Parcel theory in three dimensions and the calculation of SCAPE

Slantwise convective available potential energy (SCAPE) is a measure of the degree to which the atmosphere is unstable to conditional symmetric instability (CSI). It has, until now, been defined by parcel theory in which the atmosphere is assumed to be nonevolving and balanced, that is, two-dimensional. When applying this two-dimensional theory to three-dimensional evolving flows, these assumptions can be interpreted as an implicit assumption that a timescale separation exists between a relatively rapid timescale for slantwise ascent and a slower timescale for the development of the system. An approximate extension of parcel theory to three dimensions is derived and it is shown that calculations of SCAPE based on the assumption of relatively rapid slantwise ascent can be qualitatively in error. For a case study example of a developing extratropical cyclone, SCAPE calculated along trajectories determined without assuming the existence of the timescale separation show large SCAPE values for parcels ascending from the warm sector and along the warm front. These parcels ascend into the cloud head within which there is some evidence consistent with the release of CSI from observational and model cross sections. This region of high SCAPE was not found for calculations along the relatively rapidly ascending trajectories determined by assuming the existence of the timescale separation.

The F test for comparing two normal variances: correct and incorrect calculation of the two-sided p-value?

This article illustrates that not all statistical software packages are correctly calculating a p-value for the classical F test comparison of two independent Normal variances. This is illustrated with a simple example, and the reasons why are discussed. Eight different software packages are considered.

The determination of molecular properties from MULTIMODE with an application to the calculation of Franck-Condon factors for photoionization of CF3 to CF3+

Extensions to the code MULTIMODE to obtain rovibrational wave functions and properties are described. An application of these new capabilities is made to a calculation of the Franck-Condon factors for photoionization of CF3 to CF3+. These calculations make use of a new, full-dimensional ab initio potential energy surface, which is also described here.

A variational method for the calculation of spin-rovibronic energy levels of any triatomic molecule in an electronic triplet state

The triatomic spin-rovibronic variational code RVIB3 has been extended to include the effect of two uncoupled electrons, for both (3)Sigma(-) and (3)Pi (Renner-Teller) electronic states. The spin-orbital-rotational kinetic energy is included in the usual way, via terms (J+L+S). The phenomenological terms AL.S and lambda 2/3(3S(z)(2)) are introduced to reproduce the 3 spin-orbit and spin-spin splittings, respectively. Calculations are performed to evaluate the spin-rovibronic energy levels of CCO (X) over tilde (3) Sigma(-) and CCO (A) over tilde (3) Pi for which the Born-Oppenheimer potentials are derived from high-accuracy ab initio calculations.

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.

A DFT study of free radicals formed from artemisinin and related compounds

Ab initio calculations using density functional theory have shown that the reactions that occur between artemisinin, 1, a cyclic trioxane active against malaria, and some metal ions and complexes lead to a series of radicals which are probably responsible for its therapeutic activity. In particular it has been shown that the interaction of Fe(H) with artemisinin causes the O-O bond to be broken as indeed does Fe(III) and Cu(I), while Zn(II) does not. Calculations were carried out with Fe(II) in several different forms including the bare ion, [Fe(H2O)(5)](2+) and [FeP(Im)] (P, porphyrin; Im, imadazole) and similar results were obtained. The resulting oxygen-based radicals are readily converted to more stable carbon-based radicals and/or. stable products. Similar radicals and products are also formed from two simple model trioxanes 2 and 3 that show little or no therapeutic action against malaria although some subtle differences were obtained. This suggests that the scaffold surrounding the pharmacophore may be involved in molecular recognition events allowing efficient uptake of this trioxane warhead into the parasite. (C) 2004 Elsevier B.V. All rights reserved.

Antimalarial drugs based on artemisinin: DFT calculations on the principal reactions

Theoretical calculations have been carried out on the interactions of several endoperoxides which are potential antimalarials, including the clinically useful artemisinin, with two possible sources of iron in the parasite, namely the hexa-aquo ferrous ion [Fe(H2O)(6)](2+) and haeme. DFT calculations show that the reactions of all endoperoxides considered, with both sources of iron, initially generate a Fe-O bond followed by cleavage of the O-O bond to oxygen radical species. Subsequently, they can be transformed into carbon-centred radicals of greater stability. However, with [Fe(H2O)(6)](2+) as the iron source, the oxygen-centred radical species are more likely to react further akin to Fenton's reagent, whereby iron salts encourage hydrogen peroxide to act as an oxidizing agent, and that solvent plays a major role. In contrast, when reacting with haeme, the oxygen-centred radicals interconvert to more stable carbon-centred radicals, which can then alkylate haeme. Subsequent cleavage of the Fe-O bond leads to stable and inactive antimalarial products. These results indicate that the reactivity of the endoperoxides as antimalarials is greater with iron hexahydrates for radical-mediated damage as opposed to haeme, which leads to unreactive species. Since only nanomolar quantities of hydrated metal ions could catalyse the reactions leading to damage to the parasites, this could be an alternative or competitive reaction responsible for the antimalarial activity. (c) 2005 Elsevier B.V. All rights reserved.

A DFT study of propylene polymerization using neutral salicyladiminato nickel(II) and palladium(II) as catalysts

Propylene polymerization using salicyladiminato metal catalalysts has been studied using density functional theory at the B3LYP/LANL2DZ level. In particular, the effects on the reaction mechanisms of changing the metal from Pd(II) to Ni(II) have been investigated. While the reaction mechanisms involving the salicyladiminato Ni(II) catalyst have been found to be similar to those established previously for the salicyladiminato Pd(II) catalyst, the nickel catalyst was found to differentiate the trans-O intermediate from the trans-.N intermediate with an energy difference of 46.63 U mol(-1) significantly more than the palladium catalyst for which the energy difference was calculated as 35.82 kJ mol(-1). The energy difference between the trans-O configuration and the trans-N configuration is decreased significantly when combining a molecule of propylene with the catalyst. For the Ni catalyst, the trans-O isomer is more stable than the trans-N isomer to a greater extent than for Pd, so that the insertion of propylene from 20 is relatively less favoured for Ni than for Pd. It is predicted that the mechanism of isomerization from 20 to 2N through a rotational transition state TS2O2N is more appropriate for the Ni catalyst system. The palladium system shows a larger preference for pi-coordination than its nickel counterpart, although the latter possesses a lower reaction barrier. It was found that the occupation of the trans-O position in the asymmetric salicyladiminato catalyst is also more favored by the alkene as it is by the alkyl so that insertion of the alkene may always start from a particular configuration so that specific products are obtained. (c) 2005 Elsevier B.V. All rights reserved.

Tuning of redox potential and visible absorption band of ruthenium(II) complexes of (benzimidazolyl) derivatives: synthesis, characterization, spectroscopic and redox properties, X-ray structures and DFT calculations

Six ruthenium(II) complexes have been prepared using the tridentate ligands 2,6-bis(benzimidazolyl) pyridine and bis(2-benzimidazolyl methyl) amine and having 2,2'-bipyridine, 2,2':6',2 ''-terpyridine, PPh3, MeCN and chloride as coligands. The crystal structures of three of the complexes trans-[Ru(bbpH(2))(PPh3)(2)(CH3CN)I(ClO4)(2) center dot 2H(2)O (2), [Ru(bbpH(2))(bpy)Cl]ClO4 (3) and [Ru(bbpH(2))(terpy)](ClO4)(2) (4) are also reported. The complexes show visible region absorption at 402-517 nm, indicating that it is possible to tune the visible region absorption by varying the ancillary ligand. Luminescence behavior of the complexes has been studied both at RT and at liquid nitrogen temperature (LNT). Luminescence of the complexes is found to be insensitive to the presence of dioxygen. Two of the complexes [Ru(bbpH(2))(bpy)Cl]ClO4 (3) and [Ru(bbpH(2))(terpy]ClO4)(2) (4) show RT emission in the NIR region, having lifetime, quantum yield and radiative constant values suitable for their application as NIR emitter in the solid state devices. The DFT calculations on these two complexes indicate that the metal t(2g) electrons are appreciably delocalized over the ligand backbone. (C) 2006 Elsevier B.V. All rights reserved.

Theoretical calculation of absorption intensities of C2H and C2D

The theory of dipole-allowed absorption intensities in triatomic molecules is presented for systems with three close-lying electronic states of doublet multiplicity. Its derivation is within the framework of a recently developed variational method [CARTER, S., HANDY, N. C., PUZZARINI, C., TARRONI, R., and PALMIERI, P., 2000, Molec. Phys., 98,1967]. The method has been applied to the calculation of the infrared absorption spectrum of the C2H radical and its deuterated isotopomer for energies up to 10000 cm(-1) above the ground state, using highly accurate ab initio diabatic potential energy and dipole moment surfaces. The calculated spectra agree very well with those recorded experimentally in a neon matrix [FORNEY, D., JACOX, M. E., and THOMPSON, W. E., 1995, J. molee. Spectrosc., 170, 178] and assignments in the high energy region of the IR spectra are proposed for the first time.

The Structure of the chiral Pt{531} surface: a combined LEED and DFT study

The structure of the chiral kinked Pt{531} surface has been determined by low-energy electron diffraction intensity-versus-energy (LEED-IV) analysis and density functional theory (DFT). Large contractions and expansions of the vertical interlayer distances with respect to the bulk-terminated surface geometry were found for the first six layers (LEED: d(12) = 0.44 angstrom, d(23) = 0.69 angstrom, d(34) = 0.49 angstrom, d(45) = 0.95 angstrom, d(56) = 0.56 angstrom; DFT: d(12) = 0.51 angstrom, d(23) = 0.55 angstrom, d(34) = 0.74 angstrom, d(45) = 0.78 angstrom, d(56) = 0.63 angstrom; d(bulk) = 0.66 angstrom). Energy-dependent cancellations of LEED spots over unusually large energy ranges, up to 100 eV, can be explained by surface roughness and reproduced by applying a model involving 0.25 ML of vacancies and adatoms in the scattering calculations. The agreement between the results from LEED and DFT is not as good as in other cases, which could be due to this roughness of the real surface.

Ligand dynamics in myoglobin: calculation of infrared spectra for photodissociated NO

We present molecular dynamics simulations of the photodissociated state of MbNO performed at 300 K using a fluctuating charge model for the nitric oxide (NO) ligand. After dissociation, NO is observed to remain mainly in the centre of the distal haem pocket, although some movement towards the primary docking site and the xenon-4 pocket can be seen. We calculate the NO infrared spectrum for the photodissociated ligand within the haem pocket and find a narrow peak in the range 1915-1922 cm(-1). The resulting blue shift of 1 to 8 cm(-1) compared to gas-phase NO is much smaller than the red shifts calculated and observed for carbon monoxide (CO) in Mb. A small splitting, due to NO in the xenon-4 pocket, is also observed. At lower temperatures, the spectra and conformational space explored by the ligand remain largely unchanged, but the electrostatic interactions with residue His64 become increasingly significant in determining the details of the ligand orientation within the distal haem pocket. The investigation of the effect of the L29F mutation reveals significant differences between the behaviour of NO and that of CO, and suggests a coupling between the ligand and the protein dynamics due to the different ligand dipole moments.

Efficient on the fly calculation of time correlation functions in computer simulations

Time correlation functions yield profound information about the dynamics of a physical system and hence are frequently calculated in computer simulations. For systems whose dynamics span a wide range of time, currently used methods require significant computer time and memory. In this paper, we discuss the multiple-tau correlator method for the efficient calculation of accurate time correlation functions on the fly during computer simulations. The multiple-tau correlator is efficacious in terms of computational requirements and can be tuned to the desired level of accuracy. Further, we derive estimates for the error arising from the use of the multiple-tau correlator and extend it for use in the calculation of mean-square particle displacements and dynamic structure factors. The method described here, in hardware implementation, is routinely used in light scattering experiments but has not yet found widespread use in computer simulations.

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

Molecular structures and electronic transitions of 3,6-Diphenyl-1,2-dithiin and its radical cation: a spectroelectrochemical and DFT study

One-electron oxidation of 3,6-diphenyl-1,2-dithiin yields the corresponding radical cation. The product is stable at low temperatures and can be distinguished by a triplet EPR signal. Cyclic voltammetric, UV-vis spectroelectrochemical, and DFT studies were performed to elucidate its molecular structure and electronic properties. Time-dependent DFT calculations reproduce appreciably well the UV-vis spectral changes observed during the oxidation. The results reveal a moderately twisted structure of the 1,2-dithiin heterocycle in the radical cation.