What (science for) adaptation to climate change in Colombian agriculture? A commentary on “A way forward on adaptation to climate change in Colombian agriculture: perspectives towards 2050” by J. Ramirez-Villegas, M. Salazar, A. Jarvis, C. E. Navarro-Valcines.

Climate change is putting Colombian agriculture under significant stress and, if no adaptation is made, the latter will be severely impacted during the next decades. Ramirez-Villegas et al. (2012) set out a government-led, top-down, techno-scientific proposal for a way forward by which Colombian agriculture could adapt to climate change. However, this proposal largely overlooks the root causes of vulnerability of Colombian agriculture, and of smallholders in particular. I discuss some of the hidden assumptions underpinning this proposal and of the arguments employed by Ramirez-Villegas et al., based on existing literature on Colombian agriculture and the wider scientific debate on adaptation to climate change. While technical measures may play an important role in the adaptation of Colombian agriculture to climate change, I question whether the actions listed in the proposal alone and specifically for smallholders, truly represent priority issues. I suggest that by i) looking at vulnerability before adaptation, ii) contextualising climate change as one of multiple exposures, and iii) truly putting smallholders at the centre of adaptation, i.e. to learn about and with them, different and perhaps more urgent priorities for action can be identified. Ultimately, I argue that what is at stake is not only a list of adaptation measures but, more importantly, the scientific approach from which priorities for action are identified. In this respect, I propose that transformative rather than technical fix adaptation represents a better approach for Colombian agriculture and smallholders in particular, in the face of climate change.

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

Linear viscoelasticity from molecular dynamics simulation of entangled polymers

The linear viscoelastic (LVE) spectrum is one of the primary fingerprints of polymer solutions and melts, carrying information about most relaxation processes in the system. Many single chain theories and models start with predicting the LVE spectrum to validate their assumptions. However, until now, no reliable linear stress relaxation data were available from simulations of multichain systems. In this work, we propose a new efficient way to calculate a wide variety of correlation functions and mean-square displacements during simulations without significant additional CPU cost. Using this method, we calculate stress−stress autocorrelation functions for a simple bead−spring model of polymer melt for a wide range of chain lengths, densities, temperatures, and chain stiffnesses. The obtained stress−stress autocorrelation functions were compared with the single chain slip−spring model in order to obtain entanglement related parameters, such as the plateau modulus or the molecular weight between entanglements. Then, the dependence of the plateau modulus on the packing length is discussed. We have also identified three different contributions to the stress relaxation:  bond length relaxation, colloidal and polymeric. Their dependence on the density and the temperature is demonstrated for short unentangled systems without inertia.

Hierarchical modeling of entangled polymers

We demonstrate that it is possible to link multi-chain molecular dynamics simulations with the tube model using a single chain slip-links model as a bridge. This hierarchical approach allows significant speed up of simulations, permitting us to span the time scales relevant for a comparison with the tube theory. Fitting the mean-square displacement of individual monomers in molecular dynamics simulations with the slip-spring model, we show that it is possible to predict the stress relaxation. Then, we analyze the stress relaxation from slip-spring simulations in the framework of the tube theory. In the absence of constraint release, we establish that the relaxation modulus can be decomposed as the sum of contributions from fast and longitudinal Rouse modes, and tube survival. Finally, we discuss some open questions regarding possible future directions that could be profitable in rendering the tube model quantitative, even for mildly entangled polymers

LDA or GGA? A combined experimental inelastic neutron scattering and ab initio lattice dynamics study of alkali metal hydrides

In a previous work, we carried out inelastic neutron scattering (INS) spectroscopy experiments and preliminary first principles calculations on alkali metal hydrides. The complete series of alkali metal hydrides, LiH, NaH, KH, RbH and CsH was measured in the high-resolution TOSCA INS spectrometer at ISIS. Here, we present the results of ab initio electronic structure calculations of the properties of the alkali metal hydrides using both the local density approximation (LDA) and the generalized gradient approximation (GGA), using the Perdew–Burke–Ernzerhof (PBE) parameterization. Properties calculated were lattice parameters, bulk moduli, dielectric constants, effective charges, electronic densities and inelastic neutron scattering (INS) spectra. We took advantage of the currently available computer power to use full lattice dynamics theory to calculate thermodynamic properties for these materials. For the alkali metal hydrides (LiH, NaH, KH, RbH and CsH) using lattice dynamics, we found that the INS spectra calculated using LDA agreed better with the experimental data than the spectra calculated using GGA. Both zero-point effects and thermal contributions to free energies had an important effect on INS and several thermodynamic properties.

The vibrational spectrum and ultimate modulus of polyethylene

We have performed the first completely ab initio lattice dynamics calculation of the full orthorhombic cell of polyethylene using periodic density functional theory in the local density approximation (LDA) and the generalized gradient approximation (GGA). Contrary to current perceptions, we show that LDA accurately describes the structure whereas GGA fails. We emphasize that there is no parametrization of the results. We then rigorously tested our calculation by computing the phonon dispersion curves across the entire Brillouin zone and comparing them to the vibrational spectra, in particular the inelastic neutron scattering (INS) spectra, of polyethylene (both polycrystalline and aligned) and perdeuteriopolyethylene. The F-point frequencies (where the infrared and Raman active modes occur) are in good agreement with the latest low temperature data. The near-perfect reproduction of the INS spectra, gives confidence in the results and allows Lis to deduce a number of physical properties including the elastic moduli, parallel and perpendicular to the chain. We find that the Young's modulus for an infinitely long, perfectly crystalline polyethylene is 360.2 GPa at 0 K. The highest experimental value is 324 GPa, indicating that current high modulus fibers are similar to 90% of their maximum possible strength.

Molecular aluminum hydrides identified by inelastic neutron scattering during H-2 regeneration of catalyst-doped NaAlH4

Catalyst-doped sodium aluminum hydrides have been intensively studied as solid hydrogen carriers for onboard proton-exchange membrane (PEM) fuel cells. Although the importance of catalyst choice in enhancing kinetics for both hydrogen uptake and release of this hydride material has long been recognized, the nature of the active species and the mechanism of catalytic action are unclear. We have shown by inelastic neutron scattering (INS) spectroscopy that a volatile molecular aluminum hydride is formed during the early stage of H-2 re-eneration of a depleted, catalyst-doped sodium aluminum hydride. Computational modeling of the INS spectra suggested the formation of AlH3 and oligomers (AlH3)(n) (Al2H6, Al3H9, and Al4H12 clusters), which are pertinent to the mechanism of hydrogen storage. This paper demonstrates, for the first time, the existence of these volatile species.

Inelastic neutron scattering in spectroscopic studies of hydrogen on carbon-supported catalysts-experimental spectra and computed spectra of model systems

Inelastic neutron scattering spectroscopy has been used to observe and characterise hydrogen on the carbon component of a Pt/C catalyst. INS provides the complete vibration spectrum of coronene, regarded as a molecular model of a graphite layer. The vibrational modes are assigned with the aid of ab initio density functional theory calculations and the INS spectra by the a-CLIMAX program. A spectrum for which the H modes of coronene have been computationally suppressed, a carbon-only coronene spectrum, is a better representation of the spectrum of a graphite layer than is coronene itself. Dihydrogen dosing of a Pt/C catalyst caused amplification of the surface modes of carbon, an effect described as H riding on carbon. From the enhancement of the low energy carbon modes (100-600 cm(-1)) it is concluded that spillover hydrogen becomes attached to dangling bonds at the edges of graphitic regions of the carbon support. (C) 2003 Elsevier Science B.V. All rights reserved.

Hydrogen spillover on carbon-supported metal catalysts studied by inelastic neutron scattering. Surface vibrational states and hydrogen riding modes

Hydrogen spillover on carbon-supported precious metal catalysts has been investigated with inelastic neutron scattering (INS) spectroscopy. The aim, which was fully realized, was to identify spillover hydrogen on the carbon support. The inelastic neutron scattering spectra of Pt/C, Ru/C, and PtRu/C fuel cell catalysts dosed with hydrogen were determined in two sets of experiments: with the catalyst in the neutron beam and, using an annular cell, with carbon in the beam and catalyst pellets at the edge of the cell excluded from the beam. The vibrational modes observed in the INS spectra were assigned with reference to the INS of a polycyclic aromatic hydrocarbon, coronene, taken as a molecular model of a graphite layer, and with the aid of computational modeling. Two forms of spillover hydrogen were identified: H at edge sites of a graphite layer (formed after ambient dissociative chemisorption of H-2), and a weakly bound layer of mobile H atoms (formed by surface diffusion of H atoms after dissociative chernisorption of H-2 at 500 K). The INS spectra exhibited characteristic riding modes of H on carbon and on Pt or Ru. In these riding modes H atoms move in phase with vibrations of the carbon and metal lattices. The lattice modes are amplified by neutron scattering from the H atoms attached to lattice atoms. Uptake of hydrogen, and spillover, was greater for the Ru containing catalysts than for the Pt/C catalyst. The INS experiments have thus directly demonstrated H spillover to the carbon support of these metal catalysts.

Hydrogen adsorption in a copper ZSM5 zeolite: an inelastic neutron scattering study

Currently microporous oxidic materials including zeolites are attracting interest as potential hydrogen storage materials. Understanding how molecular hydrogen interacts with these materials is important in the rational development of hydrogen storage materials and is also challenging theoretically. In this paper, we present an incoherent inelastic neutron scattering (INS) study of the adsorption of molecular hydrogen and hydrogen deuteride (HD) in a copper substituted ZSM5 zeolite varying the hydrogen dosage and temperature. We have demonstrated how inelastic neutron scattering can help us understand the interaction of H-2 molecules with a binding site in a particular microporous material, Cu ZSM5, and by implication of other similar materials. The H-2 molecule is bound as a single species lying parallel with the surface. As H-2 dosing increases, lateral interactions between the adsorbed H-2 molecules become apparent. With rising temperature of measurement up to 70 K (the limit of our experiments), H-2 molecules remain bound to the surface equivalent to a liquid or solid H-2 phase. The implication is that hydrogen is bound rather strongly in Cu ZSM5. Using the simple model for the anisotropic interaction to calculate the energy levels splitting, we found that the measured rotational constant of the hydrogen molecule is reduced as a consequence of adsorption by the Cu ZSM5. From the decrease in total signal intensity with increasing temperature, we were able to observe the conversion of para-hydrogen into ortho-hydrogen at paramagnetic centres and so determine the fraction of paramagnetic sites occupied by hydrogen molecules, ca. 60%. (c) 2006 Elsevier B.V. All rights reserved.

Dihydrogen in cation-substituted zeolites X - an inelastic neutron scattering study

An inelastic neutron scattering (INS) study of the rotational - vibrational spectrum of dihydrogen sorbed by zeolite X having substituted sodium, calcium and zinc cations is reported. The rotational - vibrational spectrum of H-2 was observed at low energy transfer ( below ca. 25 meV, 202 cm(-1)); the vibration was that of the H-2 molecule against the binding site (H-2 - X, not H - H). The vibration frequency was proportional to the polarising power of the cation (Na+ < Ca2+ < Zn2+). Polarisation of the H-2 molecule dominated the interaction of H-2 with this binding site. The total scattering intensity was proportional to the dihydrogen dose. However the vibrational intensities became constant at ca. 0.3 wt% showing that the H-2 binding sites had saturated. Additional dihydrogen appeared as unbound or weakly bound dihydrogen exhibiting recoil.

Dihydrogen in zeolite CaX: an inelastic neutron scattering study

We report an inelastic neutron scattering (INS) study of the rotational–vibrational spectrum of dihydrogen sorbed by zeolite CaX. In the low energy (<200 cm−1) INS spectrum of adsorbed H2 we observe the rotational–vibrational spectrum of H2, where the vibration is that of the H2 molecule against the binding site (i.e. H2–X, not H–H). We have observed for the first time the vibrational overtones of the hydrogen molecule against the adsorption surface up to sixth order. These vibrations are usually forbidden in INS spectroscopy because of the selection rules imposed by the spin flip event required. In our case we are able to observe such a vibration because the rotational transition J(1 ← 0) convolutes the vibrational spectrum. This paper reports the effect for the first time.