Effect of headgroup size, charge, and solvent structure on polymer−micelle interactions, studied by molecular dynamics simulations

We performed atomistic molecular dynamics simulations of anionic and cationic micelles in the presence of poly(ethylene oxide) (PEO) to understand why nonionic water-soluble polymers such as PEO interact strongly with anionic micelles but only weakly with cationic micelles. Our micelles include sodium n-dodecyl sulfate (SDS), n-dodecyl trimethylammonium chloride (DTAC), n-dodecyl ammonium chloride (DAC), and micelles in which we artificially reverse the sign of partial charges in SDS and DTAC. We observe that the polymer interacts hydrophobically with anionic SDS but only weakly with cationic DTAC and DAC, in agreement with experiment. However, the polymer also interacts with the artificial anionic DTAC but fails to interact hydrophobically with the artificial cationic SDS, illustrating that large headgroup size does not explain the weak polymer interaction with cationic micelles. In addition, we observe through simulation that this preference for interaction with anionic micelles still exists in a dipolar "dumbbell" solvent, indicating that water structure and hydrogen bonding alone cannot explain this preferential interaction. Our simulations suggest that direct electrostatic interactions between the micelle and polymer explain the preference for interaction with anionic micelles, even though the polymer overall carries no net charge. This is possible given the asymmetric distribution of negative charges on smaller atoms and positive charges oil larger units in the polymer chain.

Exploring the climatic impact of the continental vegetation on the Mezosoic atmospheric CO2 and climate history

In this contribution, we continue our exploration of the factors defining the Mesozoic climatic history. We improve the Earth system model GEOCLIM designed for long term climate and geochemical reconstructions by adding the explicit calculation of the biome dynamics using the LPJ model. The coupled GEOCLIM-LPJ model thus allows the simultaneous calculation of the climate with a 2-D spatial resolution, the coeval atmospheric CO2, and the continental biome distribution. We found that accounting for the climatic role of the continental vegetation dynamics (albedo change, water cycle and surface roughness modulations) strongly affects the reconstructed geological climate. Indeed the calculated partial pressure of atmospheric CO2 over the Mesozoic is twice the value calculated when assuming a uniform constant vegetation. This increase in CO2 is triggered by a global cooling of the continents, itself triggered by a general increase in continental albedo owing to the development of desertic surfaces. This cooling reduces the CO2 consumption through silicate weathering, and hence results in a compensating increase in the atmospheric CO2 pressure. This study demonstrates that the impact of land plants on climate and hence on atmospheric CO2 is as important as their geochemical effect through the enhancement of chemical weathering of the continental surface. Our GEOCLIM-LPJ simulations also define a climatic baseline for the Mesozoic, around which exceptionally cool and warm events can be identified.

Magnetic coupling in trinuclear partial cubane copper(II) complexes with a hydroxo bridging core and peripheral phenoxo bridges from NNO donor Schiff base ligands

Three new trinuclear copper(II) complexes, [(CuL1)(3)(mu(3)-OH)](ClO4)(2)center dot 3.75H(2)O (1), [(CuL2)(3)(mu(3)-OH)](ClO4)(2) (2) and [(CuL3)(3)(mu(3)-OH)](BF4)(2)center dot 0.5CH(3)CN (3) have been synthesized from three tridentate Schiff bases HL1, HL2, and HL3 (HL1 = 2-[(2-amino-ethylimino)-methyl]-phenol, HL2 = 2-[(2-methylamino-ethylimino)-methyl]-phenol and HL3 = 2-[1-(2-dimethylamino-ethylimino)-ethyl]-phenol). The complexes are characterized by single-crystal X-ray diffraction analyses, IR, UV-vis and EPR spectroscopy, and variable-temperature magnetic measurements. All the compounds contain a partial cubane [Cu3O4] core consisting of the trinuclear unit [(CuL)(3)(mu(3)-OH)](2+) together with perchlorate or fluoroborate anions. In each of the complexes, the three copper atoms are five-coordinated with a distorted square-pyramidal geometry except in complex 1, in which one of the Cu-II ions of the trinuclear unit is six-coordinate being in addition weakly coordinated to one of the perchlorate anions. Variable-temperature magnetic measurements and EPR spectra indicate an antiferromagnetic exchange coupling between the CuII ions of complexes 1 and 2, while this turned out to be ferromagnetic for complex 3. Experimental values have been fitted according to an isotropic exchange Hamiltonian. Calculations based on Density Functional Theory have also been performed in order to estimate the exchange coupling constants in these three complexes. Both sets of values indicate similar trends and specially calculated J values establish a magneto-structural correlation between them and the Cu-O-Cu bond angle, in that the coupling is more ferromagnetic for smaller bond angle values.

Can climate trends be calculated from reanalysis data?

Several global quantities are computed from the ERA40 reanalysis for the period 1958-2001 and explored for trends. These are discussed in the context of changes to the global observing system. Temperature, integrated water vapor (IWV), and kinetic energy are considered. The ERA40 global mean temperature in the lower troposphere has a trend of +0.11 K per decade over the period of 1979-2001, which is slightly higher than the MSU measurements, but within the estimated error limit. For the period 1958 2001 the warming trend is 0.14 K per decade but this is likely to be an artifact of changes in the observing system. When this is corrected for, the warming trend is reduced to 0.10 K per decade. The global trend in IWV for the period 1979-2001 is +0.36 mm per decade. This is about twice as high as the trend determined from the Clausius-Clapeyron relation assuming conservation of relative humidity. It is also larger than results from free climate model integrations driven by the same observed sea surface temperature as used in ERA40. It is suggested that the large trend in IWV does not represent a genuine climate trend but an artifact caused by changes in the global observing system such as the use of SSM/I and more satellite soundings in later years. Recent results are in good agreement with GPS measurements. The IWV trend for the period 1958-2001 is still higher but reduced to +0.16 mm per decade when corrected for changes in the observing systems. Total kinetic energy shows an increasing global trend. Results from data assimilation experiments strongly suggest that this trend is also incorrect and mainly caused by the huge changes in the global observing system in 1979. When this is corrected for, no significant change in global kinetic energy from 1958 onward can be found.

Significance of cross correlations in the stress relaxation of polymer melts

According to linear response theory, all relaxation functions in the linear regime can be obtained using time correlation functions calculated under equilibrium. In this paper, we demonstrate that the cross correlations make a significant contribution to the partial stress relaxation functions in polymer melts. We present two illustrations in the context of polymer rheology using (1) Brownian dynamics simulations of a single chain model for entangled polymers, the slip-spring model, and (2) molecular dynamics simulations of a multichain model. Using the single chain model, we analyze the contribution of the confining potential to the stress relaxation and the plateau modulus. Although the idea is illustrated with a particular model, it applies to any single chain model that uses a potential to confine the motion of the chains. This leads us to question some of the assumptions behind the tube theory, especially the meaning of the entanglement molecular weight obtained from the plateau modulus. To shed some light on this issue, we study the contribution of the nonbonded excluded-volume interactions to the stress relaxation using the multichain model. The proportionality of the bonded/nonbonded contributions to the total stress relaxation (after a density dependent "colloidal" relaxation time) provides some insight into the success of the tube theory in spite of using questionable assumptions. The proportionality indicates that the shape of the relaxation spectrum can indeed be reproduced using the tube theory and the problem is reduced to that of finding the correct prefactor. (c) 2007 American Institute of Physics

Is the OECD acute worm toxicity test environmentally relevant? The effect of mineral form on calculated lead toxicity

In a series of experiments the toxicity of lead to worms in soil was determined following the draft OECD earthworm reproduction toxicity protocol except that lead was added as solid lead nitrate, carbonate and sulphide rather than as lead nitrate solution as would normally be the case. The compounds were added to the test soil to give lead concentrations of 625-12500 pg Pb g-1 of soil. Calculated toxicities of the lead decreased in the order nitrate > carbonate > sulphide, the same order as the decrease in the solubility of the metal compounds used. The 7-day LC50 (lethal concentration when 50% of the population is killed) for the nitrate was 5321 +/- 275 mug Pb g(-1) of soil and this did not change with time. The LC50 values for carbonate and sulphide could not be determined at the concentration ranges used. The only parameter sensitive enough to distinguish the toxicities of the three compounds was cocoon (egg) production. The EC50S for cocoon production (the concentration to produce a 50% reduction in cocoon production) were 993, 8604 and 10 246 mug Pb g(-1) of soil for lead nitrate, carbonate and sulphide, respectively. Standard toxicity tests need to take into account the form in which the contaminant is present in the soil to be of environmental relevance. (C) 2002 Elsevier Science Ltd. All rights reserved.

Hydraulic modeling of runoff over a rough surface under partial inundation

This study presents a numerical method to derive the Darcy- Weisbach friction coefficient for overland flow under partial inundation of surface roughness. To better account for the variable influence of roughness with varying levels of emergence, we model the flow over a network which evolves as the free surface rises. This network is constructed using a height numerical map, based on surface roughness data, and a discrete geometry skeletonization algorithm. By applying a hydraulic model to the flows through this network, local heads, velocities, and Froude and Reynolds numbers over the surface can be estimated. These quantities enable us to analyze the flow and ultimately to derive a bulk friction factor for flow over the entire surface which takes into account local variations in flow quantities. Results demonstrate that although the flow is laminar, head losses are chiefly inertial because of local flow disturbances. The results also emphasize that for conditions of partial inundation, flow resistance varies nonmonotonically but does generally increase with progressive roughness inundation.

Dehydration does not affect the radial pressures produced by the earthworm Aporrectodea caliginosa

As a soil dries, the earthworms in that soil dehydrate and become less active. Moisture stress may weaken an earthworm, lowering the radial pressure that the animal can produce. This possibility was investigated for the earthworm Aporrectodea caliginosa (Savigny). Pressures were compared for saturated earthworms (worms taken from saturated soil) and stressed earthworms (worms that had been partially dehydrated by leaving them in dry soil). A load cell was used to record the forces that earthworms produced as they moved through artificial burrows (holes that had been drilled through blocks of aluminium or Perspex). The radial pressure was calculated using the forces exerted and the dimensions of the artificial burrows. There was a negative correlation between burrow diameter and radial pressure, although radial pressure was independent of the length of the block through which the earthworms had burrowed. The highest radial pressures were produced by the anterior segments of the animal. Partial dehydration caused the earthworms to become quiescent, but did not decrease the radial pressure that the earthworms produced. It is suggested that coelomic fluid is retained in the anterior segments while the rest of the animal dehydrates. Dehydrated earthworms became lethargic, and we suggest that lethargy is due to the loss of coelomic fluid from the posterior segments. Coelomic fluid is known to be lost through dorsal pores. In burrowing species of earthworm such as Aporrectodea caliginosa, these pores are only present on the posterior segments.

Scale-invariant moving finite elements for nonlinear partial differential equations in two dimensions

A scale-invariant moving finite element method is proposed for the adaptive solution of nonlinear partial differential equations. The mesh movement is based on a finite element discretisation of a scale-invariant conservation principle incorporating a monitor function, while the time discretisation of the resulting system of ordinary differential equations is carried out using a scale-invariant time-stepping which yields uniform local accuracy in time. The accuracy and reliability of the algorithm are successfully tested against exact self-similar solutions where available, and otherwise against a state-of-the-art h-refinement scheme for solutions of a two-dimensional porous medium equation problem with a moving boundary. The monitor functions used are the dependent variable and a monitor related to the surface area of the solution manifold. (c) 2005 IMACS. Published by Elsevier B.V. All rights reserved.

Calculated rotational and rovibrational spectra of D2S and HDS

Rovibrational energy levels, transition frequencies, and linestrengths are computed variationally for the sulfur hydrides D2S and HDS, using ab initio potential energy and dipole surfaces. Wave-numbers for the pure rotational transitions agree to within 0.2 cm−1 of the experimental lines. For the fundamental vibrational transitions, the band origins for D2S are 860.4, 1900.6, and 1912.0 cm−1 for ν2, ν1, and ν3, respectively, compared with the corresponding experimental values of 855.4, 1896.4, and 1910.2 cm−1. For HDS, we compute ν2 to be 1039.4 cm−1, compared with the experimental value of 1032.7 cm−1. The relative merits of local and normal mode descriptions for the overtone stretching band origins are discussed. Our results confirm the local mode nature of the H2S, D2S, and HDS system.

Inclusion of nonlinear demand-supply relationships within large-scale partial equilibrium linear programming models

This paper presents a new method for the inclusion of nonlinear demand and supply relationships within a linear programming model. An existing method for this purpose is described first and its shortcomings are pointed out before showing how the new approach overcomes those difficulties and how it provides a more accurate and 'smooth' (rather than a kinked) approximation of the nonlinear functions as well as dealing with equilibrium under perfect competition instead of handling just the monopolistic situation. The workings of the proposed method are illustrated by extending a previously available sectoral model for the UK agriculture.

Forage production and nutrient availability in small-scale dairy systems in central Mexico using linear programming and partial budgeting

Small-scale dairy systems play an important role in the Mexican dairy sector and farm planning activities related to resource allocation have a significant impact on the profitability of such enterprises. Linear programming is a technique widely used for planning and ration formulation, and partial budgeting is a technique for assessing the impact of changes on the profitability of an enterprise. This study used both methods to optimise land use for forage production and nutrient availability, and to evaluate the economic impact of such changes in small-scale Mexican dairy systems. The model showed satisfactory performance when optimal solutions were compared with the traditional strategy. The strategy using fresh ryegrass, maize silage and oat hay, and the strategy using a combination of alfalfa hay, maize silage, fresh ryegrass and oat hay appeared attractive options for providing a better nutrient supply and maintaining a higher stocking rate throughout the year than the traditional strategy.