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.

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.

The effect of impaired thyroid function during development on the mechanical properties of avian bone

Thyroid hormones show fluctuating levels during the post-hatching development of birds. In this paper we report the results of the first mechanical tests to quantify the effect of hypothyroidism, during post-natal development, on the skeletal properties of a precocial bird, the barnacle goose, as determined by microhardness testing. The effect of hypothyroidism is tissue-specific; bone from the femora of birds is not significantly affected by induced hypothyroidism, however, there is a strong positive relationship between the levels of circulating thyroid hormones and the mechanical properties of bone from humeri. In the barnacle goose the development of the wing skeleton and musculature depends on an increase in circulating thyroid hormones and our analysis shows that, in its absence, the mechanical competence of the bone mineral itself is reduced in addition to the decreased bone length and muscle development previously reported in the literature. (C) 2004 Wiley-Liss, Inc.

Mechanical effects of plant cell wall enzymes on xyloglucan/cellulose composites

Xyloglucan-acting enzymes are believed to have effects on type I primary plant cell wall mechanical properties. In order to get a better understanding of these effects, a range of enzymes with different in vitro modes of action were tested against cell wall analogues (bio-composite materials based on Acetobacter xylinus cellulose and xyloglucan). Tomato pericarp xyloglucan endo transglycosylase (tXET) and nasturtium seed xyloglucanase (nXGase) were produced heterologously in Pichia pastoris. Their action against the cell wall analogues was compared with that of a commercial preparation of Trichoderma endo-glucanase (EndoGase). Both 'hydrolytic' enzymes (nXGase and EndoGase) were able to depolymerise not only the cross-link xyloglucan fraction but also the surface-bound fraction. Consequent major changes in cellulose fibril architecture were observed. In mechanical terms, removal of xyloglucan cross-links from composites resulted in increased stiffness (at high strain) and decreased visco-elasticity with similar extensibility. On the other hand, true transglycosylase activity (tXET) did not affect the cellulose/xyloglucan ratio. No change in composite stiffness or extensibility resulted, but a significant increase in creep behaviour was observed in the presence of active tXET. These results provide direct in vitro evidence for the involvement of cell wall xyloglucan-specific enzymes in mechanical changes underlying plant cell wall re-modelling and growth processes. Mechanical consequences of tXET action are shown to be complimentary to those of cucumber expansin.

Comparison of mechanical properties of tension and opposite wood in Populus

In this paper we focused on the differences of mechanical properties of tension and normal wood of 1-year-old poplar trees, artificially tilted. Elastic and fracture properties have been measured and linked to the anatomy. Tension wood is well known because it prevents good surface finishing and leads to difficulties with sawing. We studied three main mechanical properties: young modulus, energy of cutting and longitudinal residual strain of maturation (with strain gauges) because of their importance in wood technology. Moreover, this work takes place in a larger project of study, the phenomena of axes re-orientation in trees (allowing by the production of reaction wood), where these data are required for biomechanical modelling. The results show that tension wood has a higher young modulus, needs a higher energy to be cut and exhibited a higher level of longitudinal residual strain of maturation than those of normal wood. The results suggest that these differences require deeper analysis of the wood than anatomy: measurement of microfibril orientation in the S2 layer and also the lignin composition in monomeric units.

Interaction of heat-moisture conditions and physical properties in oat processing: I. Mechanical properties of steamed oat groats

Research interest in oats has focussed on their nutritional value, but there have been few studies of their food processing. Heat treatment is characteristic of oat processing, as it is needed to inactivate lipase and to facilitate flaking. A Texture Analyser was used to characterise the mechanical properties of unkilned and kilned oat groats after steaming and tempering in an oven for 30, 60 and 90 min at 80, 95 and 110 degrees C. Maximum force, number of peaks before maximum and final force after 5s hold were used to characterise the behaviour of the groats during compression. Kilned groats were larger and softer before steaming. After steaming and tempering, the moisture content of the kilned groats was higher than for unkilned groats. Hot, steamed oats were softer than cold, unsteamed groats, indicated by a decrease in maximum force from 59 to 55 N, and there was no significant difference between kilned and unkilned groats. However, higher temperatures during tempering increased maximum force. These results suggest that mild steam treatment yields softer oat groats, whereas cold or over-treated groats tend to be harder. (c) 2007 Elsevier Ltd. All rights reserved.

Influence of some processing and storage conditions on the mechanical properties of oat flakes

Flake breakage and texture are important quality, criteria in oat flakes. These properties are determined by the mechanical properties of the flakes, which may be influenced by process variables such as kilning and flake thickness. A pin deformation method was used to measure the rupture force of individual oat flakes at different water activities. The monolayer value of ground oat flakes ranged from 5.83 to 684 g/100 g dry matter Thick flakes were strongest, requiring 3.4 N to rupture the flake compared to 2.2 N for the thin flakes. Water softened the flakes, causing a decrease in rupture force from 3.6 N to 2.4 N as water activity; increased from 0.115 to 0.848. Kilning had a significant effect on flake thickness but not on the mechanical properties. This study suggests that oat flakes should be stored at water activity 0.4 or less as there is a sharp loss of flake strength above this point.

Mechanical properties of oats and oat products

The aim of this review is to illustrate how physical properties are important to food processing and quality. Three food products, flakes, porridge and bread, in addition to oat groats are used to show the influence of water and heat-treatments on the mechanical properties. The hydrothermal history of ingredients is shown to affect product quality. Water acts as a plasticiser and solvent in these foods, whilst heat modifies the conformation and interactions of macromolecular components. Structure as well as chemical composition is shown to govern texture.

Self-assembled healable materials through a combination of pi-pi stacking and hydrogen bonding

The discovery of polymers with stimuli responsive physical properties is a rapidly expanding area of research. At the forefront of the field are self-healing polymers, which, when fractured can regain the mechanical properties of the material either autonomically, or in response to a stimulus. It has long been known that it is possible to promote healing in conventional thermoplastics by heating the fracture zone above the Tg of the polymer under pressure. This process requires reptation and subsequent re-entanglement of macromolecules across the fracture void, which serves to bridge, and ‘heal’ the crack. The timescale for this mechanism is highly dependent on the molecular weight of the polymer being studied. This process is in contrast to that required to affect healing in supramolecular polymers such as the plasticised, hydrogen bonded elastomer reported by Leibler et al. The disparity in bond energies between the non-covalent and covalent bonds within supramolecular polymers results in fractures propagating through scission of the comparatively weak supramolecular interactions, rather than through breaking the stronger, covalent bonds. Thus, during the healing process the macromolecules surrounding the fracture site only need sufficient energy to re-engage their supramolecular interactions in order to regenerate the strength of the pristine material. Herein we describe the design, synthesis and optimization of a new class of supramolecular polymer blends that harness the reversible nature of pi-pi stacking and hydrogen bonding interactions to produce self-supporting films with facile healable characteristics.

A statistical mechanical approach for the computation of the climatic response to general forcings

The climate belongs to the class of non-equilibrium forced and dissipative systems, for which most results of quasi-equilibrium statistical mechanics, including the fluctuation-dissipation theorem, do not apply. In this paper we show for the first time how the Ruelle linear response theory, developed for studying rigorously the impact of perturbations on general observables of non-equilibrium statistical mechanical systems, can be applied with great success to analyze the climatic response to general forcings. The crucial value of the Ruelle theory lies in the fact that it allows to compute the response of the system in terms of expectation values of explicit and computable functions of the phase space averaged over the invariant measure of the unperturbed state. We choose as test bed a classical version of the Lorenz 96 model, which, in spite of its simplicity, has a well-recognized prototypical value as it is a spatially extended one-dimensional model and presents the basic ingredients, such as dissipation, advection and the presence of an external forcing, of the actual atmosphere. We recapitulate the main aspects of the general response theory and propose some new general results. We then analyze the frequency dependence of the response of both local and global observables to perturbations having localized as well as global spatial patterns. We derive analytically several properties of the corresponding susceptibilities, such as asymptotic behavior, validity of Kramers-Kronig relations, and sum rules, whose main ingredient is the causality principle. We show that all the coefficients of the leading asymptotic expansions as well as the integral constraints can be written as linear function of parameters that describe the unperturbed properties of the system, such as its average energy. Some newly obtained empirical closure equations for such parameters allow to define such properties as an explicit function of the unperturbed forcing parameter alone for a general class of chaotic Lorenz 96 models. We then verify the theoretical predictions from the outputs of the simulations up to a high degree of precision. The theory is used to explain differences in the response of local and global observables, to define the intensive properties of the system, which do not depend on the spatial resolution of the Lorenz 96 model, and to generalize the concept of climate sensitivity to all time scales. We also show how to reconstruct the linear Green function, which maps perturbations of general time patterns into changes in the expectation value of the considered observable for finite as well as infinite time. Finally, we propose a simple yet general methodology to study general Climate Change problems on virtually any time scale by resorting to only well selected simulations, and by taking full advantage of ensemble methods. The specific case of globally averaged surface temperature response to a general pattern of change of the CO2 concentration is discussed. We believe that the proposed approach may constitute a mathematically rigorous and practically very effective way to approach the problem of climate sensitivity, climate prediction, and climate change from a radically new perspective.

Molecular dynamics simulations of proteins: can the explicit water model be varied?

In molecular mechanics simulations of biological systems, the solvation water is typically represented by a default water model which is an integral part of the force field. Indeed, protein nonbonding parameters are chosen in order to obtain a balance between water-water and protein-water interactions and hence a reliable description of protein solvation. However, less attention has been paid to the question of whether the water model provides a reliable description of the water properties under the chosen simulation conditions, for which more accurate water models often exist. Here we consider the case of the CHARMM protein force field, which was parametrized for use with a modified TIP3P model. Using quantum mechanical and molecular mechanical calculations, we investigate whether the CHARMM force field can be used with other water models: TIP4P and TIP5P. Solvation properties of N-methylacetamide (NMA), other small solute molecules, and a small protein are examined. The results indicate differences in binding energies and minimum energy geometries, especially for TIP5P, but the overall description of solvation is found to be similar for all models tested. The results provide an indication that molecular mechanics simulations with the CHARMM force field can be performed with water models other than TIP3P, thus enabling an improved description of the solvent water properties.

Structure–activity relations (SARs) for gas-phase reactions of NO3, OH and O3 with alkenes: an update

We present updated structure-activity relations (SARs) for the prediction of rate coefficients for gas-phase reactions with alkenes of the major atmospheric oxidants NO3, OH and O-3. Such SARs provide one way of incorporating essential information about reactivity into atmospheric models. Rate coefficients obtained from correlations relating the logarithms of the rate coefficients to the energies of the highest occupied molecular orbitals (HOMOs) of the alkenes were used to refine the SARs. SARs have an advantage for the user over the direct application of the correlations in that knowledge of the structure of the alkene of interest is sufficient to estimate rate coefficients, and no quantum-mechanical calculations need to be performed. A comparison of the values predicted by the SARs with experimental data where they exist allowed us to assess the reliability of our method.