The Role of Anharmonicity in Hydrogen-Bonded Systems: The Case of Water Clusters

The nature of vibrational anharmonicity has been examined for the case of small water clusters using second-order vibrational perturbation theory (VPT2) applied on second-order Møller–Plesset perturbation theory (MP2) potential energy surfaces. Using a training set of 16 water clusters (H2O)n=2–6,8,9 with a total of 723 vibrational modes, we determined scaling factors that map the harmonic frequencies onto anharmonic ones. The intermolecular modes were found to be substantially more anharmonic than intramolecular bending and stretching modes. Due to the varying levels of anharmonicity of the intermolecular and intramolecular modes, different frequency scaling factors for each region were necessary to achieve the highest accuracy. Furthermore, new scaling factors for zero-point vibrational energies (ZPVE) and vibrational corrections to the enthalpy (ΔHvib) and the entropy (Svib) have been determined. All the scaling factors reported in this study are different from previous works in that they are intended for hydrogen-bonded systems, while others were built using experimental frequencies of covalently bonded systems. An application of our scaling factors to the vibrational frequencies of water dimer and thermodynamic functions of 11 larger water clusters highlights the importance of anharmonic effects in hydrogen-bonded systems.

Structures, energetics and fragmentation pathways of CnH22+ carbodications

Potential energy curves have been calculated for CnH22+ (n = 2−9) ions and results have been compared with data on unimolecular charge-separation reactions obtained by Rabrenović and Beynon. Geometry-optimized, minimum energy, linear CnH22+ structures have been computed for ground and low-lying excited states. These carbodications exist in stable configurations with well depths greater than 3 eV. Decomposition pathways into singly charged fragment ions lead to products with computed kinetic energies in excess of 1 eV. A high degree of correlation exists between experimental information and results computed for linear CnH22+ structures having hydrogen atoms on each end. The exception involves C4H22+reactions where a low-lying doubly charged isomer must be invoked to rationalize the experimental data.

Modeling of temperature and turbidity in a natural lake and a reservoir connected by pumped-storage operations

Pumped-storage (PS) systems are used to store electric energy as potential energy for release during peak demand. We investigate the impacts of a planned 1000 MW PS scheme connecting Lago Bianco with Lago di Poschiavo (Switzerland) on temperature and particle mass concentration in both basins. The upper (turbid) basin is a reservoir receiving large amounts of fine particles from the partially glaciated watershed, while the lower basin is a much clearer natural lake. Stratification, temperature and particle concentrations in the two basins were simulated with and without PS for four different hydrological conditions and 27 years of meteorological forcing using the software CE-QUAL-W2. The simulations showed that the PS operations lead to an increase in temperature in both basins during most of the year. The increase is most pronounced (up to 4°C) in the upper hypolimnion of the natural lake toward the end of summer stratification and is partially due to frictional losses in the penstocks, pumps and turbines. The remainder of the warming is from intense coupling to the atmosphere while water resides in the shallower upper reservoir. These impacts are most pronounced during warm and dry years, when the upper reservoir is strongly heated and the effects are least concealed by floods. The exchange of water between the two basins relocates particles from the upper reservoir to the lower lake, where they accumulate during summer in the upper hypolimnion (10 to 20 mg L−1) but also to some extent decrease light availability in the trophic surface layer.

Photodissociation dynamics of tert-butylnitrite following excitation to the S(1) and S(2) states. A study by velocity-map ion-imaging and 3D-REMPI spectroscopy

Excitation of tert-butylnitrite into the first and second UV absorption bands leads to efficient dissociation into the fragment radicals NO and tert-butoxy in their electronic ground states (2)Π and (2)E, respectively. Velocity distributions and angular anisotropies for the NO fragment in several hundred rotational and vibrational quantum states were obtained by velocity-map imaging and the recently developed 3D-REMPI method. Excitation into the well resolved vibronic progression bands (k = 0, 1, 2) of the NO stretch mode in the S(1) ← S(0) transition produces NO fragments mostly in the vibrational state with v = k, with smaller fractions in v = k - 1 and v = k - 2. It is concluded that dissociation occurs on the purely repulsive PES of S(1) without barrier. All velocity distributions from photolysis via the S(1)(nπ*) state are monomodal and show high negative anisotropy (β ≈ -1). The rotational distributions peak near j = 30.5 irrespective of the vibronic state S(1)(k) excited and the vibrational state v of the NO fragment. On average 46% of the excess energy is converted to kinetic energy, 23% and 31% remain as internal energy in the NO fragment and the t-BuO radical, respectively. Photolysis via excitation into the S(2) ← S(0) transition at 227 nm yields NO fragments with about equal populations in v = 0 and v = 1. The rotational distributions have a single maximum near j = 59.5. The velocity distributions are monomodal with positive anisotropy β ≈ 0.8. The average fractions of the excess energy distributed into translation, internal energy of NO, and internal energy of t-BuO are 39%, 23%, and 38%, respectively. In all cases ∼8500 cm(-1) of energy remain in the internal degrees of freedom of the t-BuO fragment. This is mostly assigned to rotational energy. An ab initio calculation of the dynamic reaction path shows that not only the NO fragment but also the t-BuO fragment gain large angular momentum during dissociation on the purely repulsive potential energy surface of S(2).

Experimental and Calculated Spectra of π-Stacked Mild Charge-Transfer Complexes: Jet-Cooled Perylene·(Tetrachloroethene) n , n = 1,2

The S0 ↔ S1 spectra of the mild charge-transfer (CT) complexes perylene·tetrachloroethene (P·4ClE) and perylene·(tetrachloroethene)2 (P·(4ClE)2) are investigated by two-color resonant two-photon ionization (2C-R2PI) and dispersed fluorescence spectroscopy in supersonic jets. The S0 → S1 vibrationless transitions of P·4ClE and P·(4ClE)2 are shifted by δν = −451 and −858 cm–1 relative to perylene, translating to excited-state dissociation energy increases of 5.4 and 10.3 kJ/mol, respectively. The red shift is ∼30% larger than that of perylene·trans-1,2-dichloroethene; therefore, the increase in chlorination increases the excited-state stabilization and CT character of the interaction, but the electronic excitation remains largely confined to the perylene moiety. The 2C-R2PI and fluorescence spectra of P·4ClE exhibit strong progressions in the perylene intramolecular twist (1au) vibration (42 cm–1 in S0 and 55 cm–1 in S1), signaling that perylene deforms along its twist coordinate upon electronic excitation. The intermolecular stretching (Tz) and internal rotation (Rc) vibrations are weak; therefore, the P·4ClE intermolecular potential energy surface (IPES) changes little during the S0 ↔ S1 transition. The minimum-energy structures and inter- and intramolecular vibrational frequencies of P·4ClE and P·(4ClE)2 are calculated with the dispersion-corrected density functional theory (DFT) methods B97-D3, ωB97X-D, M06, and M06-2X and the spin-consistent-scaled (SCS) variant of the approximate second-order coupled-cluster method, SCS-CC2. All methods predict the global minima to be π-stacked centered coplanar structures with the long axis of tetrachloroethene rotated by τ ≈ 60° relative to the perylene long axis. The calculated binding energies are in the range of −D0 = 28–35 kJ/mol. A second minimum is predicted with τ ≈ 25°, with ∼1 kJ/mol smaller binding energy. Although both monomers are achiral, both the P·4ClE and P·(4ClE)2 complexes are chiral. The best agreement for adiabatic excitation energies and vibrational frequencies is observed for the ωB97X-D and M06-2X DFT methods.

The elusive S2 state, the S1/S2 splitting, and the excimer states of the benzene dimer

We observe the weak S 0 → S 2 transitions of the T-shaped benzene dimers (Bz)2 and (Bz-d 6)2 about 250 cm−1 and 220 cm−1 above their respective S 0 → S 1 electronic origins using two-color resonant two-photon ionization spectroscopy. Spin-component scaled (SCS) second-order approximate coupled-cluster (CC2) calculations predict that for the tipped T-shaped geometry, the S 0 → S 2 electronic oscillator strength f el (S 2) is ∼10 times smaller than f el (S 1) and the S 2 state lies ∼240 cm−1 above S 1, in excellent agreement with experiment. The S 0 → S 1 (ππ ∗) transition is mainly localized on the “stem” benzene, with a minor stem → cap charge-transfer contribution; the S 0 → S 2 transition is mainly localized on the “cap” benzene. The orbitals, electronic oscillator strengths f el (S 1) and f el (S 2), and transition frequencies depend strongly on the tipping angle ω between the two Bz moieties. The SCS-CC2 calculated S 1 and S 2 excitation energies at different T-shaped, stacked-parallel and parallel-displaced stationary points of the (Bz)2 ground-state surface allow to construct approximate S 1 and S 2 potential energy surfaces and reveal their relation to the “excimer” states at the stacked-parallel geometry. The f el (S 1) and f el (S 2) transition dipole moments at the C 2v -symmetric T-shape, parallel-displaced and stacked-parallel geometries are either zero or ∼10 times smaller than at the tipped T-shaped geometry. This unusual property of the S 0 → S 1 and S 0 → S 2 transition-dipole moment surfaces of (Bz)2 restricts its observation by electronic spectroscopy to the tipped and tilted T-shaped geometries; the other ground-state geometries are impossible or extremely difficult to observe. The S 0 → S 1/S 2 spectra of (Bz)2 are compared to those of imidazole ⋅ (Bz)2, which has a rigid triangular structure with a tilted (Bz)2 subunit. The S 0 → S 1/ S 2 transitions of imidazole-(benzene)2 lie at similar energies as those of (Bz)2, confirming our assignment of the (Bz)2 S 0 → S 2 transition.

Territorial Intelligence of Vulnerability Systems 2- Sustainable Modelling of Globalization Challenge

Globalization as progress of economic development has increased population socioeconomical vulnerability when unequal wealth distribution within economic development process constitutes the main rule, with widening the gap between rich and poors by environmental pricing. Econological vulnerability is therefore increasing too, as dangerous substance and techniques should produce polluted effluents and industrial or climatic risk increasing (Woloszyn, Quenault, Faburel, 2012). To illustrate and model this process, we propose to introduce an analogical induction-model to describe both vulnerability situations and associated resilience procedures. At this aim, we first develop a well-known late 80?s model of socio-economic crack-up, known as 'Silent Weapons for Quiet Wars', which presents economics as a social extension of natural energy systems. This last, also named 'E-model', is constituted by three passive components, potential energy, kinetic energy, and energy dissipation, thus allowing economical data to be treated as a thermodynamical system. To extend this model to social and ecological sustainability pillars, we propose to built an extended E(Economic)-S(Social)-O(Organic) model, based on the three previous components, as an open model considering feedbacks as evolution sources. An applicative illustration of this model will then be described, through this summer's american severe drought event analysis

Territorial Intelligence of Vulnerability Systems 2- Sustainable Modelling of Globalization Challenge

Globalization as progress of economic development has increased population socioeconomical vulnerability when unequal wealth distribution within economic development process constitutes the main rule, with widening the gap between rich and poors by environmental pricing. Econological vulnerability is therefore increasing too, as dangerous substance and techniques should produce polluted effluents and industrial or climatic risk increasing (Woloszyn, Quenault, Faburel, 2012). To illustrate and model this process, we propose to introduce an analogical induction-model to describe both vulnerability situations and associated resilience procedures. At this aim, we first develop a well-known late 80?s model of socio-economic crack-up, known as 'Silent Weapons for Quiet Wars', which presents economics as a social extension of natural energy systems. This last, also named 'E-model', is constituted by three passive components, potential energy, kinetic energy, and energy dissipation, thus allowing economical data to be treated as a thermodynamical system. To extend this model to social and ecological sustainability pillars, we propose to built an extended E(Economic)-S(Social)-O(Organic) model, based on the three previous components, as an open model considering feedbacks as evolution sources. An applicative illustration of this model will then be described, through this summer's american severe drought event analysis

Territorial Intelligence of Vulnerability Systems 2- Sustainable Modelling of Globalization Challenge

Globalization as progress of economic development has increased population socioeconomical vulnerability when unequal wealth distribution within economic development process constitutes the main rule, with widening the gap between rich and poors by environmental pricing. Econological vulnerability is therefore increasing too, as dangerous substance and techniques should produce polluted effluents and industrial or climatic risk increasing (Woloszyn, Quenault, Faburel, 2012). To illustrate and model this process, we propose to introduce an analogical induction-model to describe both vulnerability situations and associated resilience procedures. At this aim, we first develop a well-known late 80?s model of socio-economic crack-up, known as 'Silent Weapons for Quiet Wars', which presents economics as a social extension of natural energy systems. This last, also named 'E-model', is constituted by three passive components, potential energy, kinetic energy, and energy dissipation, thus allowing economical data to be treated as a thermodynamical system. To extend this model to social and ecological sustainability pillars, we propose to built an extended E(Economic)-S(Social)-O(Organic) model, based on the three previous components, as an open model considering feedbacks as evolution sources. An applicative illustration of this model will then be described, through this summer's american severe drought event analysis

Porewater concentrations, pH, single cell number, anaerobic oxidation of methane (AOM) rate, sulfate reduction (SR) rate of samples during METEOR cruise M76/3b at the REGAB cold seep site from the West African margin in 2008

The giant pockmark REGAB (West African margin, 3160 m water depth) is an active methane-emitting cold seep ecosystem, where the energy derived from microbially mediated oxidation of methane supports high biomass and diversity of chemosynthetic communities. Bare sediments interspersed with heterogeneous chemosynthetic assemblages of mytilid mussels, vesicomyid clams and siboglinid tubeworms form a complex seep ecosystem. To better understand if benthic bacterial communities reflect the patchy distribution of chemosynthetic fauna, all major chemosynthetic habitats at REGAB were investigated using an interdisciplinary approach combining porewater geochemistry, in situ quantification of fluxes and consumption of methane, as well bacterial community fingerprinting. This study revealed that sediments populated by different fauna assemblages show distinct biogeochemical activities and are associated with distinct sediment bacterial communities. The methane consumption and methane effluxes ranged over one to two orders of magnitude across habitats, and reached highest values at the mussel habitat, which hosted a different bacterial community compared to the other habitats. Clam assemblages had a profound impact on the sediment geochemistry, but less so on the bacterial community structure. Moreover, all clam assemblages at REGAB were restricted to sediments characterized by complete methane consumption in the seafloor, and intermediate biogeochemical activity. Overall, variations in the sediment geochemistry were reflected in the distribution of both fauna and microbial communities; and were mostly determined by methane flux.

CCSD(T) Study of CD3-O-CD3 and CH3-O-CD3 Far-Infrared Spectra

From a vibrationally corrected 3D potential energy surface determined with highly correlated ab initio calculations (CCSD(T)), the lowest vibrational energies of two dimethyl-ether isotopologues, 12CH3–16O–12CD3 (DME-d3) and 12CD3–16O–12CD3 (DME-d6), are computed variationally. The levels that can be populated at very low temperatures correspond to the COC-bending and the two methyl torsional modes. Molecular symmetry groups are used for the classification of levels and torsional splittings. DME-d6 belongs to the G36 group, as the most abundant isotopologue 12CH3–16O–12CH3 (DME-h6), while DME-d3 is a G18 species. Previous assignments of experimental Raman and far-infrared spectra are discussed from an effective Hamiltonian obtained after refining the ab initio parameters. Because a good agreement between calculated and experimental transition frequencies is reached, new assignments are proposed for various combination bands corresponding to the two deuterated isotopologues and for the 020 → 030 transition of DME-d6. Vibrationally corrected potential energy barriers, structural parameters, and anharmonic spectroscopic parameters are provided. For the 3N – 9 neglected vibrational modes, harmonic and anharmonic fundamental frequencies are obtained using second-order perturbation theory by means of CCSD and MP2 force fields. Fermi resonances between the COC-bending and the torsional modes modify DME-d3 intensities and the band positions of the torsional overtones.

Aplicación de la física estadística a los medios granulares

Esta tesis se centra en el estudio de medios granulares blandos y atascados mediante la aplicación de la física estadística. Esta aproximación se sitúa entre los tradicionales enfoques macro y micromecánicos: trata de establecer cuáles son las propiedades macroscópicas esperables de un sistema granular en base a un análisis de las propiedades de las partículas y las interacciones que se producen entre ellas y a una consideración de las restricciones macroscópicas del sistema. Para ello se utiliza la teoría estadística junto con algunos principios, conceptos y definiciones de la teoría de los medios continuos (campo de tensiones y deformaciones, energía potencial elástica, etc) y algunas técnicas de homogeneización. La interacción entre las partículas es analizada mediante las aportaciones de la teoría del contacto y de las fuerzas capilares (producidas por eventuales meniscos de líquido cuando el medio está húmedo). La idea básica de la mecánica estadística es que entre todas soluciones de un problema físico (como puede ser el ensamblaje en equilibrio estático de partículas de un medio granular) existe un conjunto que es compatible con el conocimiento macroscópico que tenemos del sistema (por ejemplo, su volumen, la tensión a la que está sometido, la energía potencial elástica que almacena, etc.). Este conjunto todavía contiene un número enorme de soluciones. Pues bien, si no hay ninguna información adicional es razonable pensar que no existe ningún motivo para que alguna de estas soluciones sea más probable que las demás. Entonces parece natural asignarles a todas ellas el mismo peso estadístico y construir una función matemática compatible. Actuando de este modo se obtiene cuál es la función de distribución más probable de algunas cantidades asociadas a las soluciones, para lo cual es muy importante asegurarse de que todas ellas son igualmente accesibles por el procedimiento de ensamblaje o protocolo. Este enfoque se desarrolló en sus orígenes para el estudio de los gases ideales pero se puede extender para sistemas no térmicos como los analizados en esta tesis. En este sentido el primer intento se produjo hace poco más de veinte años y es la colectividad de volumen. Desde entonces esta ha sido empleada y mejorada por muchos investigadores en todo el mundo, mientras que han surgido otras, como la de la energía o la del fuerza-momento (tensión multiplicada por volumen). Cada colectividad describe, en definitiva, conjuntos de soluciones caracterizados por diferentes restricciones macroscópicas, pero de todos ellos resultan distribuciones estadísticas de tipo Maxwell-Boltzmann y controladas por dichas restricciones. En base a estos trabajos previos, en esta tesis se ha adaptado el enfoque clásico de la física estadística para el caso de medios granulares blandos. Se ha propuesto un marco general para estudiar estas colectividades que se basa en la comparación de todas las posibles soluciones en un espacio matemático definido por las componentes del fuerza-momento y en unas funciones de densidad de estados. Este desarrollo teórico se complementa con resultados obtenidos mediante simulación de la compresión cíclica de sistemas granulares bidimensionales. Se utilizó para ello un método de dinámica molecular, MD (o DEM). Las simulaciones consideran una interacción mecánica elástica, lineal y amortiguada a la que se ha añadido, en algunos casos, la fuerza cohesiva producida por meniscos de agua. Se realizaron cálculos en serie y en paralelo. Los resultados no solo prueban que las funciones de distribución de las componentes de fuerza-momento del sistema sometido a un protocolo específico parecen ser universales, sino que también revelan que existen muchos aspectos computacionales que pueden determinar cuáles son las soluciones accesibles. This thesis focuses on the application of statistical mechanics for the study of static and jammed packings of soft granular media. Such approach lies between micro and macromechanics: it tries to establish what the expected macroscopic properties of a granular system are, by starting from a micromechanical analysis of the features of the particles, and the interactions between them, and by considering the macroscopic constraints of the system. To do that, statistics together with some principles, concepts and definitions of continuum mechanics (e.g. stress and strain fields, elastic potential energy, etc.) as well as some homogenization techniques are used. The interaction between the particles of a granular system is examined too and theories on contact and capillary forces (when the media are wet) are revisited. The basic idea of statistical mechanics is that among the solutions of a physical problem (e.g. the static arrangement of particles in mechanical equilibrium) there is a class that is compatible with our macroscopic knowledge of the system (volume, stress, elastic potential energy,...). This class still contains an enormous number of solutions. In the absence of further information there is not any a priori reason for favoring one of these more than any other. Hence we shall naturally construct the equilibrium function by assigning equal statistical weights to all the functions compatible with our requirements. This procedure leads to the most probable statistical distribution of some quantities, but it is necessary to guarantee that all the solutions are likely accessed. This approach was originally set up for the study of ideal gases, but it can be extended to non-thermal systems too. In this connection, the first attempt for granular systems was the volume ensemble, developed about 20 years ago. Since then, this model has been followed and improved upon by many researchers around the world, while other two approaches have also been set up: energy and force-moment (i.e. stress multiplied by volume) ensembles. Each ensemble is described by different macroscopic constraints but all of them result on a Maxwell-Boltzmann statistical distribution, which is precisely controlled by the respective constraints. According to this previous work, in this thesis the classical statistical mechanics approach is introduced and adapted to the case of soft granular media. A general framework, which includes these three ensembles and uses a force-moment phase space and a density of states function, is proposed. This theoretical development is complemented by molecular dynamics (or DEM) simulations of the cyclic compression of 2D granular systems. Simulations were carried out by considering spring-dashpot mechanical interactions and attractive capillary forces in some cases. They were run on single and parallel processors. Results not only prove that the statistical distributions of the force-moment components obtained with a specific protocol seem to be universal, but also that there are many computational issues that can determine what the attained packings or solutions are.

Highly correlated ab initio study of the far infrared spectra of methyl acetate

Highly correlated ab initio calculations (CCSD(T)) are used to compute gas phase spectroscopic parameters of three isotopologues of the methyl acetate (CH3COOCH3, CD3COOCH3, and CH3COOCD3), searching to help experimental assignments and astrophysical detections. The molecule shows two conformers cis and trans separated by a barrier of 4457 cm−1. The potential energy surface presents 18 minima that intertransform through three internal rotation motions. To analyze the far infrared spectrum at low temperatures, a three-dimensional Hamiltonian is solved variationally. The two methyl torsion barriers are calculated to be 99.2 cm−1 (C–CH3) and 413.1 cm−1 (O–CH3), for the cis-conformer. The three fundamental torsional band centers of CH3COOCH3 are predicted to lie at 63.7 cm−1 (C–CH3), 136.1 cm−1 (O–CH3), and 175.8 cm−1 (C–O torsion) providing torsional state separations. For the 27 vibrational modes, anharmonic fundamentals and rovibrational parameters are provided. Computed parameters are compared with those fitted using experimental data.