Large vibrational variational calculations using 'multimode' and an iterative diagonalization technique

We have favoured the variational (secular equation) method for the determination of the (ro-) vibrational energy levels of polyatomic molecules. We use predominantly the Watson Hamiltonian in normal coordinates and an associated given potential in the variational code 'Multimode'. The dominant cost is the construction and diagonalization of matrices of ever-increasing size. Here we address this problem, using pertubation theory to select dominant expansion terms within the Davidson-Liu iterative diagonalization method. Our chosen example is the twelve-mode molecule methanol, for which we have an ab initio representation of the potential which includes the internal rotational motion of the OH group relative to CH3. Our new algorithm allows us to obtain converged energy levels for matrices of dimensions in excess of 100 000.

Full-dimensional vibrational calculations for H5O2+ using an ab initio potential energy surface

We report quantum diffusion Monte Carlo (DMC) and variational calculations in full dimensionality for selected vibrational states of H5O2+ using a new ab initio potential energy surface [X. Huang, B. Braams, and J. M. Bowman, J. Chem. Phys. 122, 044308 (2005)]. The energy and properties of the zero-point state are focused on in the rigorous DMC calculations. OH-stretch fundamentals are also calculated using "fixed-node" DMC calculations and variationally using two versions of the code MULTIMODE. These results are compared with infrared multiphoton dissociation measurements of Yeh [L. I. Yeh, M. Okumura, J. D. Myers, J. M. Price, and Y. T. Lee, J. Chem. Phys. 91, 7319 (1989)]. Some preliminary results for the energies of several modes of the shared hydrogen are also reported.

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.

Trinuclear cu(ll) complexes containing peripheral ketonic oxygen bridges and a mu(3)-OH core: Syntheses, crystal structures, spectroscopic and magnetic properties

Four new trinuclear copper(II) complexes, [(CuL1)(3)(mu(3)-OH)](ClO4)(2)center dot H2O (1), [(CuL2)(3)(mu(3)-OH)](CIO4)(2) (2), [(CuL3)(3)-(mu(3)-OH)](ClO4)(4)center dot H2O (3), and [(CuL4)(3)(mu(3)-OH)](ClO4)(2)center dot H2O (4), where HL1 = 8-amino-4,7,7-trimethyl-5-azaoct-3-en-2-one, HL2 = 7-amino-4-methyl-5-azaoct-3-en-2-one, HL3 = 7(ethylamino)-4-methyl-5-azahept-3-en-2-one, and HL4 = 4-methyl-7-(methylamino)-5-azahept-3-en-2-one, have been derived from the four tridentate Schiff bases (HL1, HL2, HL3, and HL4) and structurally characterized by X-ray crystallography. For all compounds, the cationic part is trinuclear with a CU3OH core held by three carbonyl oxygen bridges between each pair of copper(II) atoms. The copper atoms are five-coordinate with a distorted square-pyramidal geometry; the equatorial plane consists of the bridging oxygen atom of the central OH group together with three atoms (N, N, O) from one ligand whereas an oxygen atom of a second ligand occupies the axial position. Magnetic measurements have been performed in the 2-300 K temperature range. The experimental data could be satisfactorily reproduced by using an isotropic exchange model, H = -J(S1S2+S2S3+S1S3) yielding as best-fit parameters: J = -66.7 and g = 2.19 for 1, J = -36.6 and g = 2.20 for 2, J = -24.5 and g = 2.20 for 3, and J = -14.9 and g = 2.05 for 4. EPR spectra at low temperature show the existence of spin frustration in complexes 3 and 4, but it has not been possible to carry out calculations of the antisymmetric exchange parameter, G, from magnetic data. In frozen methanolic solution, at 4 K, hyperfine splitting in all complexes and spin frustration in complex 4 seem to be confirmed. ((c) Wiley-VCH Verlag GmbH & Co. KGaA, 69451 Weinheim, Germany, 2005)

Evaluation of data for atmospheric models: master equation/RRKM calculations on the combination reaction, BrO+NO2 -> BrONO2, a conundrum

Experimental data for the title reaction were modeled using master equation (ME)/RRKM methods based on the Multiwell suite of programs. The starting point for the exercise was the empirical fitting provided by the NASA (Sander, S. P.; Finlayson-Pitts, B. J.; Friedl, R. R.; Golden, D. M.; Huie, R. E.; Kolb, C. E.; Kurylo, M. J.; Molina, M. J.; Moortgat, G. K.; Orkin, V. L.; Ravishankara, A. R. Chemical Kinetics and Photochemical Data for Use in Atmospheric Studies, Evaluation Number 15; Jet Propulsion Laboratory: Pasadena, California, 2006)(1) and IUPAC (Atkinson, R.; Baulch, D. L.; Cox, R. A.: R. F. Hampson, J.; Kerr, J. A.; Rossi, M. J.; Troe, J. J. Phys. Chem. Ref. Data. 2000, 29, 167) 2 data evaluation panels, which represents the data in the experimental pressure ranges rather well. Despite the availability of quite reliable parameters for these calculations (molecular vibrational frequencies (Parthiban, S.; Lee, T. J. J. Chem. Phys. 2000, 113, 145)3 and a. value (Orlando, J. J.; Tyndall, G. S. J. Phys. Chem. 1996, 100,. 19398)4 of the bond dissociation energy, D-298(BrO-NO2) = 118 kJ mol(-1), corresponding to Delta H-0(circle) = 114.3 kJ mol(-1) at 0 K) and the use of RRKM/ME methods, fitting calculations to the reported data or the empirical equations was anything but straightforward. Using these molecular parameters resulted in a discrepancy between the calculations and the database of rate constants of a factor of ca. 4 at, or close to, the low-pressure limit. Agreement between calculation and experiment could be achieved in two ways, either by increasing Delta H-0(circle) to an unrealistically high value (149.3 kJ mol(-1)) or by increasing , the average energy transferred in a downward collision, to an unusually large value (> 5000 cm(-1)). The discrepancy could also be reduced by making all overall rotations fully active. The system was relatively insensitive to changing the moments of inertia in the transition state to increase the centrifugal effect. The possibility of involvement of BrOONO was tested and cannot account for the difficulties of fitting the data.

Full-dimensional quantum calculations of ground-state tunneling splitting of malonaldehyde using an accurate ab initio potential energy surface

Quantum calculations of the ground vibrational state tunneling splitting of H-atom and D-atom transfer in malonaldehyde are performed on a full-dimensional ab initio potential energy surface (PES). The PES is a fit to 11 147 near basis-set-limit frozen-core CCSD(T) electronic energies. This surface properly describes the invariance of the potential with respect to all permutations of identical atoms. The saddle-point barrier for the H-atom transfer on the PES is 4.1 kcal/mol, in excellent agreement with the reported ab initio value. Model one-dimensional and "exact" full-dimensional calculations of the splitting for H- and D-atom transfer are done using this PES. The tunneling splittings in full dimensionality are calculated using the unbiased "fixed-node" diffusion Monte Carlo (DMC) method in Cartesian and saddle-point normal coordinates. The ground-state tunneling splitting is found to be 21.6 cm(-1) in Cartesian coordinates and 22.6 cm(-1) in normal coordinates, with an uncertainty of 2-3 cm(-1). This splitting is also calculated based on a model which makes use of the exact single-well zero-point energy (ZPE) obtained with the MULTIMODE code and DMC ZPE and this calculation gives a tunneling splitting of 21-22 cm(-1). The corresponding computed splittings for the D-atom transfer are 3.0, 3.1, and 2-3 cm(-1). These calculated tunneling splittings agree with each other to within less than the standard uncertainties obtained with the DMC method used, which are between 2 and 3 cm(-1), and agree well with the experimental values of 21.6 and 2.9 cm(-1) for the H and D transfer, respectively. (C) 2008 American Institute of Physics.

MULTIMODE quantum calculations of intramolecular vibrational energies of the water dimer and trimer using ab initio-based potential energy surfaces

We report vibrational configuration interaction calculations of the monomer fundamentals of (H2O)(2), (D2O)(2), (H2O)(3), and (D2O)(3) using the code MULTIMODE and full dimensional ab initio-based global potential energies surfaces (PESs). For the dimer the HBB PES [Huang , J. Chem. Phys 128, 034312 (2008)] is used and for the trimer a new PES, reported here, is used. The salient properties of the new trimer PES are presented and compared to previous single-point calculations and the vibrational energies are compared with experiments. (C) 2008 American Institute of Physics.

Tests of MULTIMODE calculations of rovibrational energies of CH4

We report variational calculations of rovibrational energies of CH4 using the code MULTIMODE and an ab initio force field of Schwenke and Partridge. The systematic convergence of the energies with respect to the level of mode coupling is presented. Converged vibrational energies calculated using the five-mode representation of the potential for zero total angular momentum are compared with previous, benchmark calculations based on Radau coordinates using this force field for zero total angular momentum and for J = 1. Very good agreement with the previous benchmark calculations is found. (c) 2006 Elsevier B.V. All rights reserved.

Verification of design calculations of a wind catcher/tower natural ventilation system with performance testing in real building

A wind catcher/tower natural ventilation system was installed in a seminar room in the building of the School of Construction Management and Engineering, the University of Reading in the UK . Performance was analysed by means of ventilation tracer gas measurements, indoor climate measurements (temperature, humidity, CO2) and occupant surveys. In addition, the potential of simple design tools was evaluated by comparing observed ventilation results with those predicted by an explicit ventilation model and the AIDA implicit ventilation model. To support this analysis, external climate parameters (wind speed and direction, solar radiation, external temperature and humidity) were also monitored. The results showed the chosen ventilation design provided a substantially greater ventilation rate than an equivalent area of openable window. Also air quality parameters stayed within accepted norms while occupants expressed general satisfaction with the system and with comfort conditions. Night cooling was maximised by using the system in combination with openable windows. Comparisons of calculations with ventilation rate measurements showed that while AIDA gave reasonably correlated results with the monitored performance results, the widely used industry explicit model was found to over estimate the monitored ventilation rate.

Fast and accurate SCFT calculations for periodic block-copolymer morphologies using the spectral method with Anderson mixing

We study the numerical efficiency of solving the self-consistent field theory (SCFT) for periodic block-copolymer morphologies by combining the spectral method with Anderson mixing. Using AB diblock-copolymer melts as an example, we demonstrate that this approach can be orders of magnitude faster than competing methods, permitting precise calculations with relatively little computational cost. Moreover, our results raise significant doubts that the gyroid (G) phase extends to infinite $\chi N$. With the increased precision, we are also able to resolve subtle free-energy differences, allowing us to investigate the layer stacking in the perforated-lamellar (PL) phase and the lattice arrangement of the close-packed spherical (S$_{cp}$) phase. Furthermore, our study sheds light on the existence of the newly discovered Fddd (O$^{70}$) morphology, showing that conformational asymmetry has a significant effect on its stability.

Geodetic mass balance of Azarova glacier, Kodar Mountains, eastern Siberia, and its links to observed and projected climatic change

The Kodar Mountains in eastern Siberia accommodate 30 small, cold-based glaciers with a combined surface area of about 19 km2. Very little is known about these glaciers, with the first survey conducted in the late 1950s. In this paper, we use terrestrial photogrammetry to calculate changes in surface area, elevation, volume and geodetic mass balance of the Azarova Glacier between 1979 and 2007 and relate these to meteorological data from nearby Chara weather station (1938-2007). The glacier surface area declined by 20±6.9% and surface lowered on average by 20±1.8 m (mean thinning: 0.71 m a-1) resulting in a strongly negative cumulative and average mass balance of -18±1.6 m w.e. and -640±60 mm w.e.a-1 respectively. The July-August air temperature increased at a rate of 0.036oC a-1 between 1979 and 2007 and the 1980-2007 period was, on average, around 1oC warmer than 1938-1979. The regional climate projections for A2 and B2 CO2 emission scenarios developed using PRECIS regional climate model indicate that summer temperatures will increase in 2071–2100 by 2.6-4.7°C and 4.9-6.2°C respectively in comparison with 1961–1990. The annual total of solid precipitation will increase by 20% under B2 scenario but decline by 3% under A2 scenario. The length of the ablation season will extend from July–August to June-September. The Azarova Glacier exhibits high sensitivity to climatic warming due to its low elevation, exposure to comparatively high summer temperatures, and the absence of a compensating impact of cold season precipitation. Further summer warming and decline of solid precipitation projected under the A2 scenario will force Azarova to retreat further while impacts of an increase in solid precipitation projected under the B2 scenario require further investigation.

A model study of factors influencing projected changes in regional sea level over the 21st century

In addition to projected increases in global mean sea level over the 21st century, model simulations suggest there will also be changes in the regional distribution of sea level relative to the global mean. There is a considerable spread in the projected patterns of these changes by current models, as shown by the recent Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment (AR4). This spread has not reduced from that given by the Third Assessment models. Comparison with projections by ensembles of models based on a single structure supports an earlier suggestion that models of similar formulation give more similar patterns of sea level change. Analysing an AR4 ensemble of model projections under a business-as-usual scenario shows that steric changes (associated with subsurface ocean density changes) largely dominate the sea level pattern changes. The relative importance of subsurface temperature or salinity changes in contributing to this differs from region to region and, to an extent, from model-to-model. In general, thermosteric changes give the spatial variations in the Southern Ocean, halosteric changes dominate in the Arctic and strong compensation between thermosteric and halosteric changes characterises the Atlantic. The magnitude of sea level and component changes in the Atlantic appear to be linked to the amount of Atlantic meridional overturning circulation (MOC) weakening. When the MOC weakening is substantial, the Atlantic thermosteric patterns of change arise from a dominant role of ocean advective heat flux changes.

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

Syntheses, x-ray structures, photochemistry, redox properties, and DFT calculations of interconvertible fac- and mer-[Mn(SPS)(CO)3] isomers containing a flexible SPS-based pincer ligand

The lithium salt of the anionic SPS pincer ligand composed of a central hypervalent lambda(4)-phosphinine ring bearing two ortho-positioned diphenylphosphine sulfide side arms reacts with [Mn(CO)(5)Br] to give fac-[Mn(SPS)(CO)(3)], This isomer can be converted photochemicaily to mer-[Mn(SPS)(CO)(3)], with a very high quantum yield (0.80 +/- 0.05). The thermal backreaction is slow (taking ca. 8 h at room temperature), in contrast to rapid electrodecatalyzed mer-to-fac isomerization triggered by electrochemical reduction of mer-[Mn(SPS)(CO)(3)]. Both geometric isomers of [Mn(SPS)(CO)(3)] have been characterized by X-ray crystallography. Both isomers show luminescence from a low-lying (IL)-I-3 (SPS-based) excited state. The light emission of fac-[Mn(SPS)(CO)(3)] is largely quenched by the efficient photoisomerization occurring probably from a low-lying Mn-CO dissociative excited state. Density functional theory (DFT) and time-dependent DFT calculations describe the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) of fac- and mer-[Mn(CO)(3)(SPS)] as ligand-centered orbitals, largely localized on the phosphinine ring of the SPS pincer ligand. In line with the ligand nature of its frontier orbitals, fac-[Mn(SPS)(CO)(3)] is electrochemically reversibly oxidized and reduced to the corresponding radical cation and anion, respectively. The spectroscopic (electron paramagnetic resonance, IR, and UV-vis) characterization of the radical species provides other evidence for the localization of the redox steps on the SIPS ligand. The smaller HOMO-LUMO energy difference in the case of mer-[Mn(CO)(3)(SPS)], reflected in the electronic absorption and emission spectra, corresponds with its lower oxidation potential compared to that of the fac isomer. The thermodynamic instability of mer-[Mn(CO)(3)(SPS)], confirmed by the DFT calculations, increases upon one-electron reduction and oxidation of the complex.