How does basis set superposition error change the potential surfaces for hydrogen-bonded dimers?

We describe a simple method to automate the geometric optimization of molecular orbital calculations of supermolecules on potential surfaces that are corrected for basis set superposition error using the counterpoise (CP) method. This method is applied to the H-bonding complexes HF/HCN, HF/H2O, and HCCH/H2O using the 6-31G(d,p) and D95 + + (d,p) basis sets at both the Hartree-Fock and second-order Møller-Plesset levels. We report the interaction energies, geometries, and vibrational frequencies of these complexes on the CP-optimized surfaces; and compare them with similar values calculated using traditional methods, including the (more traditional) single point CP correction. Upon optimization on the CP-corrected surface, the interaction energies become more negative (before vibrational corrections) and the H-bonding stretching vibrations decrease in all cases. The extent of the effects vary from extremely small to quite large depending on the complex and the calculational method. The relative magnitudes of the vibrational corrections cannot be predicted from the H-bond stretching frequencies alone

A new perspective on Southern hemisphere storm-tracks.

A detailed view of Southern Hemisphere storm tracks is obtained based on the application of filtered variance and modern feature-tracking techniques to a wide range of 45-yr European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERA-40) data. It has been checked that the conclusions drawn in this study are valid even if data from only the satellite era are used. The emphasis of the paper is on the winter season, but results for the four seasons are also discussed. Both upper- and lower-tropospheric fields are used. The tracking analysis focuses on systems that last longer than 2 days and are mobile (move more than 1000 km). Many of the results support previous ideas about the storm tracks, but some new insights are also obtained. In the summer there is a rather circular, strong, deep high-latitude storm track. In winter the high-latitude storm track is more asymmetric with a spiral from the Atlantic and Indian Oceans in toward Antarctica and a subtropical jet–related lower-latitude storm track over the Pacific, again tending to spiral poleward. At all times of the year, maximum storm activity in the higher-latitude storm track is in the Atlantic and Indian Ocean regions. In the winter upper troposphere, the relative importance of, and interplay between, the subtropical and subpolar storm tracks is discussed. The genesis, lysis, and growth rate of lower-tropospheric winter cyclones together lead to a vivid picture of their behavior that is summarized as a set of overlapping plates, each composed of cyclone life cycles. Systems in each plate appear to feed the genesis in the next plate through downstream development in the upper-troposphere spiral storm track. In the lee of the Andes in South America, there is cyclogenesis associated with the subtropical jet and also, poleward of this, cyclogenesis largely associated with system decay on the upslope and regeneration on the downslope. The genesis and lysis of cyclones and anticyclones have a definite spatial relationship with each other and with the Andes. At 500 hPa, their relative longitudinal positions are consistent with vortex-stretching ideas for simple flow over a large-scale mountain. Cyclonic systems near Antarctica have generally spiraled in from lower latitudes. However, cyclogenesis associated with mobile cyclones occurs around the Antarctic coast with an interesting genesis maximum over the sea ice near 150°E. The South Pacific storm track emerges clearly from the tracking as a coherent deep feature spiraling from Australia to southern South America. A feature of the summer season is the genesis of eastward-moving cyclonic systems near the tropic of Capricorn off Brazil, in the central Pacific and, to a lesser extent, off Madagascar, followed by movement along the southwest flanks of the subtropical anticyclones and contribution to the “convergence zone” cloud bands seen in these regions.

Microscopic Mechanisms of Strain Hardening in Glassy Polymers

The mechanisms underlying the increase in stress for large mechanical strains of a polymer glass, quantified by the strain-hardening modulus, are still poorly understood. In the present paper we aim to elucidate this matter and present new mechanisms. Molecular-dynamics simulations of two polymers with very different strain-hardening moduli (polycarbonate and polystyrene) have been carried out. Nonaffine displacements occur because of steric hindrances and connectivity constraints. We argue that it is not necessary to introduce the concept of entanglements to understand strain hardening, but that hardening is rather coupled with the increase in the rate of nonaffine particle displacements. This rate increases faster for polycarbonate, which has the higher strain-hardening modulus. Also more nonaffine chain stretching is present for polycarbonate. It is shown that the inner distances of such a nonaffinely deformed chain can be well described by the inner distances of the worm-like chain, but with an effective stiffness length (equal to the Kuhn length for an infinite worm-like chain) that increases during deformation. It originates from the finite extensibility of the chain. In this way the increase in nonaffine particle displacement can be understood as resulting from an increase in the effective stiffness length of the perturbed chain during deformation, so that at larger strains a higher rate of plastic events in terms of nonaffine displacement is necessary, causing in turn the observed strain hardening in polymer glasses.

Modelling European winter wind storm losses in current and future climate

Severe wind storms are one of the major natural hazards in the extratropics and inflict substantial economic damages and even casualties. Insured storm-related losses depend on (i) the frequency, nature and dynamics of storms, (ii) the vulnerability of the values at risk, (iii) the geographical distribution of these values, and (iv) the particular conditions of the risk transfer. It is thus of great importance to assess the impact of climate change on future storm losses. To this end, the current study employs—to our knowledge for the first time—a coupled approach, using output from high-resolution regional climate model scenarios for the European sector to drive an operational insurance loss model. An ensemble of coupled climate-damage scenarios is used to provide an estimate of the inherent uncertainties. Output of two state-of-the-art global climate models (HadAM3, ECHAM5) is used for present (1961–1990) and future climates (2071–2100, SRES A2 scenario). These serve as boundary data for two nested regional climate models with a sophisticated gust parametrizations (CLM, CHRM). For validation and calibration purposes, an additional simulation is undertaken with the CHRM driven by the ERA40 reanalysis. The operational insurance model (Swiss Re) uses a European-wide damage function, an average vulnerability curve for all risk types, and contains the actual value distribution of a complete European market portfolio. The coupling between climate and damage models is based on daily maxima of 10 m gust winds, and the strategy adopted consists of three main steps: (i) development and application of a pragmatic selection criterion to retrieve significant storm events, (ii) generation of a probabilistic event set using a Monte-Carlo approach in the hazard module of the insurance model, and (iii) calibration of the simulated annual expected losses with a historic loss data base. The climate models considered agree regarding an increase in the intensity of extreme storms in a band across central Europe (stretching from southern UK and northern France to Denmark, northern Germany into eastern Europe). This effect increases with event strength, and rare storms show the largest climate change sensitivity, but are also beset with the largest uncertainties. Wind gusts decrease over northern Scandinavia and Southern Europe. Highest intra-ensemble variability is simulated for Ireland, the UK, the Mediterranean, and parts of Eastern Europe. The resulting changes on European-wide losses over the 110-year period are positive for all layers and all model runs considered and amount to 44% (annual expected loss), 23% (10 years loss), 50% (30 years loss), and 104% (100 years loss). There is a disproportionate increase in losses for rare high-impact events. The changes result from increases in both severity and frequency of wind gusts. Considerable geographical variability of the expected losses exists, with Denmark and Germany experiencing the largest loss increases (116% and 114%, respectively). All countries considered except for Ireland (−22%) experience some loss increases. Some ramifications of these results for the socio-economic sector are discussed, and future avenues for research are highlighted. The technique introduced in this study and its application to realistic market portfolios offer exciting prospects for future research on the impact of climate change that is relevant for policy makers, scientists and economists.

Chronology for the Aegean Late Bronze Age 1700-1400 BC

Radiocarbon (carbon-14) data from the Aegean Bronze Age 1700-1400 B.C. show that the Santorini (Thera) eruption must have occurred in the late 17th century B.C. By using carbon-14 dates from the surrounding region, cultural phases, and Bayesian statistical analysis, we established a chronology for the initial Aegean Late Bronze Age cultural phases (Late Minoan IA, IB, and II). This chronology contrasts with conventional archaeological dates and cultural synthesis: stretching out the Late Minoan IA, IB, and II phases by similar to 100 years and requiring reassessment of standard interpretations of associations between the Egyptian and Near Eastern historical dates and phases and those in the Aegean and Cyprus in the mid-second millennium B.C.

Monte Carlo phase diagram for diblock copolymer melts

A partial phase diagram is constructed for diblock copolymer melts using lattice-based Monte Carlo simulations. This is done by locating the order-disorder transition (ODT) with the aid of a recently proposed order parameter and identifying the ordered phase over a wide range of copolymer compositions (0.2 <= f <= 0.8). Consistent with experiments, the disordered phase is found to exhibit direct first-order transitions to each of the ordered morphologies. This includes the spontaneous formation of a perforated-lamellar phase, which presumably forms in place of the gyroid morphology due to finite-size and/or nonequilibrium effects. Also included in our study is a detailed examination of disordered cylinder-forming (f=0.3) diblock copolymers, revealing a substantial degree of pretransitional chain stretching and short-range order that set in well before the ODT, as observed previously in analogous studies on lamellar-forming (f=0.5) molecules. (c) 2006 American Institute of Physics.

Moist singular vectors and the predictability of some high impact European cyclones

The ECMWF full-physics and dry singular vector (SV) packages, using a dry energy norm and a 1-day optimization time, are applied to four high impact European cyclones of recent years that were almost universally badly forecast in the short range. It is shown that these full-physics SVs are much more relevant to severe cyclonic development than those based on dry dynamics plus boundary layer alone. The crucial extra ingredient is the representation of large-scale latent heat release. The severe winter storms all have a long, nearly straight region of high baroclinicity stretching across the Atlantic towards Europe, with a tongue of very high moisture content on its equatorward flank. In each case some of the final-time top SV structures pick out the region of the actual storm. The initial structures were generally located in the mid- to low troposphere. Forecasts based on initial conditions perturbed by moist SVs with opposite signs and various amplitudes show the range of possible 1-day outcomes for reasonable magnitudes of forecast error. In each case one of the perturbation structures gave a forecast very much closer to the actual storm than the control forecast. Deductions are made about the predictability of high-impact extratropical cyclone events. Implications are drawn for the short-range forecast problem and suggestions made for one practicable way to approach short-range ensemble forecasting. Copyright © 2005 Royal Meteorological Society.

On the relationship of normal modes to local modes in molecular vibrations

A simple model for the effective vibrational hamiltonian of the XH stretching vibrations in H2O, NH3 and CH4 is considered, based on a morse potential function for the bond stretches plus potential and kinetic energy coupling between pairs of bond oscillators. It is shown that this model can be set up as a matrix in local mode basis functions, or as a matrix in normal mode basis functions, leading to identical results. The energy levels obtained exhibit normal mode patterns at low vibrational excitation, and local mode patterns at high excitation. When the hamiltonian is set up in the normal mode basis it is shown that Darling-Dennison resonances must be included, and simple relations are found to exist between the xrs, gtt, and Krrss anharmonic constants (where the Darling-Dennison coefficients are denoted K) due to their contributions from morse anharmonicity in the bond stretches. The importance of the Darling-Dennison resonances is stressed. The relationship of the two alternative representations of this local mode/normal mode model are investigated, and the potential uses and limitations of the model are discussed.

A model for stretch-bend interactions in molecular force fields

Interaction force constants between bond-stretching and angle-bending co-ordinates in polyatomic molecules have been attributed, by some authors, to changes of hybridization due to orbital-following of the bending co-ordinate, and consequent changes of bond length due to the change of hybridization. A method is described for using this model quantitatively to reduce the number of independent force constants in the potential function of a polyatomic molecule, by relating stretch-bend interaction constants to the corresponding diagonal stretching constants. It is proposed to call this model the Hybrid Orbital Force Field. The model is applied to the tetrahedral four co-ordinated carbon atom (as in methane) and to the trigonal planar three coordinated carbon atom (as in formaldehyde).

Vibration–rotation variational calculations: precise results on HCN up to 25 000 cm−1

Variation calculations of the vibration–rotation energy levels of many isotopomers of HCN are reported, for J=0, 1, and 2, extending up to approximately 8 quanta of each of the stretching vibrations and 14 quanta of the bending mode. The force field, which is represented as a polynomial expansion in Morse coordinates for the bond stretches and even powers of the angle bend, has been refined by least squares to fit simultaneously all observed data on the Σ and Π state vibrational energies, and the Σ state rotational constants, for both HCN and DCN. The observed vibrational energies are fitted to roughly ±0.5 cm−1, and the rotational constants to roughly ±0.0001 cm−1. The force field has been used to predict the vibration rotation spectra of many isotopomers of HCN up to 25 000 cm−1. The results are consistent with the axis‐switching assignments of some weak overtone bands reported recently by Jonas, Yang, and Wodtke, and they also fit and provide the assignment for recent observations by Romanini and Lehmann of very weak absorption bands above 20 000 cm−1.

Local modes and x,K relations for ethene and propadiene

The results recently obtained by Mills and Robiette on local-mode effects in H2O, NH3 and CH4 type molecules are extended to ethene (C2H4) and propadiene (C3H4) type molecules. General relations among the anharmonic xrs constants and the Darling-Dennison Krrss constants for the stretching vibrations are derived, called “x,K relations”, which allow local-mode effects to be generated by adding the appropriate anharmonic and Darling-Dennison constants to the familiar normal-mode model of molecular vibrations. The general utility of x,K relations is discussed, and the results are reviewed for the molecular types for which they have so far been derived.

Born-Oppenheimer failure in the separation of low frequency vibrations

Changes in the effective potential function of a low-frequency large-amplitude molecular vibration, resulting from excitation of a high-frequency vibration, are discussed. It is shown that in some situations a significant contribution to such changes may arise from failure of the Born-Oppenheimer separation of the low-frequency mode. In the particular example of the HF dimer, recent evidence that the tunneling barrier increases on exciting either of the H-stretching vibrations is probably due to this effect.

Excited vibrational state microwave spectra, structure, and force-field of trifluorosilane

The J = 2−1 microwave spectrum of six isotopic species of HSiF3 has been observed and assigned in excited states of five of the six fundamental vibrations. The assignment is based on relative intensities, double resonance experiments, and trial anharmonic force constant calculations. Analysis of the spectra leads to experimental values for five of the constants, all three l-doubling constants qt, one Fermi resonance constant φ233, and one zeta constant. The harmonic force field has been refined to all the available data on vibration wavenumbers, centrifugal distortion constants, and zeta constants. The cubic anharmonic force field has been refined to the data on and qt constants, using two models: a valence force model with two cubic force constants for SiH and SiF stretching, and a more sophisticated model. With the help of these calculations, the following equilibrium structure has been determined: re(SiH) = 1.4468(±5) Å, re(SiF) = 1.5624(±1) Å, HSiF = 110.64(±3)°,