Vibrational dynamics and structure of natural feldspars

The present study gives a contribution to the knowledge on the Na-feldspar and plagioclases, extending the database of the Raman spectra of plagioclases with different chemical compositions and structural orders. This information may be used for the future planetary explorations by “rovers”, for the investigation of ceramics nanocrystal materials and for the mineralogical phase identification in sediments. Na-feldspar and plagioclase solid solution have been investigated by Raman spectroscopy in order to determine the relationships between the vibrational changes and the plagioclase crystal chemistry and structure. We focused on the Raman micro-spectroscopy technique, being a non-destructive method, suited for contactless analysis with high spatial resolution. Chemical and structural analyses have been performed on natural samples to test the usefulness of Raman spectroscopy as a tool in the study of the pressure-induced structural deformations, the disordering processes due to change in the Al-Si distribution in the tetrahedral sites and, finally, in the determination of the anorthitic content (Anx) in plagioclase minerals. All the predicted 39 Ag Raman active modes have been identified and assigned to specific patterns of atomic vibrational motion. A detailed comparison between experimental and computed Raman spectra has been performed and previous assignments have been revised, solving some discrepancies reported in recent literature. The ab initio calculation at the hybrid HF/DFT level with the WC1LYP Hamiltonian has proven to give excellent agreement between calculated and experimentally measured Raman wavenumbers and intensities in triclinic minerals. A short digression on the 36 infrared active modes of Na-feldspar has been done too. The identification of all 39 computed Raman modes in the experimentally measured spectra of the fully ordered Na-feldspar, known as low albite, along with the detailed description of each vibrational mode, has been essential to extend the comparative analysis to the high pressure and high temperature structural forms of albite, which reflect the physical–chemical conditions of the hosting rocks. The understanding of feldspar structure response to pressure and temperature is crucial in order to constrain crustal behaviour. The compressional behaviour of the Na-feldspar has been investigated for the first time by Raman spectroscopy. The absence of phase transitions and the occurrence of two secondary compression mechanisms acting at different pressures have been confirmed. Moreover, Raman data suggest that the internal structural changes are confined to a small pressure interval, localized around 6 GPa, not spread out from 4 to 8 GPa as suggested by previous X-rays studies on elasticity. The dominant compression mechanisms act via tetrahedral tilting, while the T-O bond lengths remain nearly constant at moderate compressional regimes. At the spectroscopic level, this leads to the strong pressure dependencies of T-O-T bending modes, as found for the four modes at 478, 508, 578 and 815 cm-1. The Al-Si distribution in the tetrahedral sites affects also the Raman spectrum of Na-feldspar. In particular, peak broadening is more sensitive than peak position to changes in the degree of order. Raman spectroscopy is found to be a good probe for local ordering, in particular being sensitive to the first annealing steps, when the macroscopic order parameter is still high. Even though Raman data are scattered and there are outliers in the estimated values of the degree of order, the average peak linewidths of the Na-feldspar characteristic doublet band, labelled here as υa and υb, as a function of the order parameter Qod show interesting trends: both peak linewidths linearly increase until saturation. From Qod values lower than 0.6, peak broadening is no more affected by the Al-Si distribution. Moreover, the disordering process is found to be heterogeneous. SC-XRD and Raman data have suggested an inter-crystalline inhomogeneity of the samples, i.e., the presence of regions with different defect density on the micrometric scale. Finally, the influence of Ca-Na substitution in the plagioclase Raman spectra has been investigated. Raman spectra have been collected on a series of well characterized natural, low structural plagioclases. The variations of the Raman modes as a function of the chemical composition and the structural order have been determined. The number of the observed Raman bands at each composition gives information about the unit-cell symmetry: moving away from the C1 structures, the number of the Raman bands enhances, as the number of formula units in the unit cell increases. The modification from an “albite-like” Raman spectrum to a more “anorthite-like” spectrum occurs from sample An78 onwards, which coincides with the appearance of c reflections in the diffraction patterns of the samples. The evolution of the Raman bands υa and υb displays two changes in slope at ~An45 and ~An75: the first one occurs between e2 and e1 plagioclases, the latter separates e1 and I1 plagioclases with only b reflections in their diffraction patterns from I1 and P1 samples having b and c reflections too. The first variation represents exactly the e2→e1 phase transitions, whereas the second one corresponds in good approximation to the C1→I1 transition, which has been determined at ~An70 by previous works. The I1→P1 phase transition in the anorthite-rich side of the solid solution is not highlighted in the collected Raman spectra. Variations in peak broadening provide insights into the behaviour of the order parameter on a local scale, suggesting an increase in the structural disorder within the solid solution, as the structures have to incorporate more Al atoms to balance the change from monovalent to divalent cations. All the information acquired on these natural plagioclases has been used to produce a protocol able to give a preliminary estimation of the chemical composition of an unknown plagioclase from its Raman spectrum. Two calibration curves, one for albite-rich plagioclases and the other one for the anorthite-rich plagioclases, have been proposed by relating the peak linewidth of the most intense Raman band υa and the An content. It has been pointed out that the dependence of the composition from the linewidth can be obtained only for low structural plagioclases with a degree of order not far away from the references. The proposed tool has been tested on three mineralogical samples, two of meteoric origin and one of volcanic origin. Chemical compositions by Raman spectroscopy compare well, within an error of about 10%, with those obtained by elemental techniques. Further analyses on plagioclases with unknown composition will be necessary to validate the suggested method and introduce it as routine tool for the determination of the chemical composition from Raman data in planetary missions.

Computer simulation study of collective dynamics in the glass former Ca(NO3)(2)center dot 4H(2)O

Time correlation functions of current fluctuations were calculated by molecular dynamics (MD) simulations in order to investigate sound waves of high wavevectors in the glass-forming liquid Ca(NO3)(2)center dot 4H(2)O. Dispersion curves, omega(k), were obtained for longitudinal (LA) and transverse acoustic (TA) modes, and also for longitudinal optic (LO) modes. Spectra of LA modes calculated by MD simulations were modeled by a viscoelastic model within the memory function framework. The viscoelastic model is used to rationalize the change of slope taking place at k similar to 0.3 angstrom(-1) in the omega(k) curve of acoustic modes. For still larger wavevectors, mixing of acoustic and optic modes is observed. Partial time correlation functions of longitudinal mass currents were calculated separately for the ions and the water molecules. The wavevector dependence of excitation energies of the corresponding partial LA modes indicates the coexistence of a relatively stiff subsystem made of cations and anions, and a softer subsystem made of water molecules. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4751548]

Ultrafast detection and control of molecular dynamics

Many elementary chemical and physical processes such as the breaking of a chemical bond or the vibrational motion of atoms within a molecule take place on a femtosecond (fs = 10−15 s) or picosecond (ps = 10−12 s) time scale. It is now possible to monitor these events as a function of time with temporal resolution well below 100 fs. This capability is based on the pump-probe technique where one optical pulse triggers a reaction and a second delayed optical pulse probes the changes that ensue. To illustrate this capability, the dynamics of ligand motion within a protein are presented. Moving beyond casual observation of a reaction to active control of its outcome requires additional experimental and theoretical effort. To illustrate the concept of control, the effect of optical pulse duration on the vibrational dynamics of a tri-atomic molecule are discussed. The experimental and theoretical resources currently available are poised to make the dream of reaction control a reality for certain molecular systems.

Inferring a weighted elastic network from partial unfolding with coarse-grained simulations

A number of studies have demonstrated that simple elastic network models can reproduce experimental B-factors, providing insights into the structure-function properties of proteins. Here, we report a study on how to improve an elastic network model and explore its performance by predicting the experimental B-factors. Elastic network models are built on the experimental C coordinates, and they only take the pairs of C atoms within a given cutoff distance r(c) into account. These models describe the interactions by elastic springs with the same force constant. We have developed a method based on numerical simulations with a simple coarse-grained force field, to attribute weights to these spring constants. This method considers the time that two C atoms remain connected in the network during partial unfolding, establishing a means of measuring the strength of each link. We examined two different coarse-grained force fields and explored the computation of these weights by unfolding the native structures. Proteins 2014; 82:119-129. (c) 2013 Wiley Periodicals, Inc.

C–H local modes in cyclobutene. II. Laser photoacoustic studies 10 000–17 000 cm−1. Vibrational structure and C–H local mode dynamics

In part I of this study [Baggott, Clase, and Mills, Spectrochim. Acta Part A 42, 319 (1986)] we presented FTIR spectra of gas phase cyclobutene and modeled the v=1–3 stretching states of both olefinic and methylenic C–H bonds in terms of a local mode model. In this paper we present some improvements to our original model and make use of recently derived ‘‘x,K relations’’ to find the equivalent normal mode descriptions. The use of both the local mode and normal mode approaches to modeling the vibrational structure is described in some detail. We present evidence for Fermi resonance interactions between the methylenic C–H stretch overtones and ring C–C stretch vibrations, revealed in laser photoacoustic spectra in the v=4–6 region. An approximate model vibrational Hamiltonian is proposed to explain the observed structure and is used to calculate the dynamics of the C–H stretch local mode decay resulting from interaction with lower frequency ring modes. The implications of our experimental and theoretical studies for mode‐selective photochemistry are discussed briefly.

Molecular dynamics simulations of silicate and borate glasses and melts : structure, diffusion dynamics and vibrational properties

Molecular dynamics simulations of silicate and borate glasses and melts: Structure, diffusion dynamics and vibrational properties. In this work computer simulations of the model glass formers SiO2 and B2O3 are presented, using the techniques of classical molecular dynamics (MD) simulations and quantum mechanical calculations, based on density functional theory (DFT). The latter limits the system size to about 100−200 atoms. SiO2 and B2O3 are the two most important network formers for industrial applications of oxide glasses. Glass samples are generated by means of a quench from the melt with classical MD simulations and a subsequent structural relaxation with DFT forces. In addition, full ab initio quenches are carried out with a significantly faster cooling rate. In principle, the structural properties are in good agreement with experimental results from neutron and X-ray scattering, in all cases. A special focus is on the study of vibrational properties, as they give access to low-temperature thermodynamic properties. The vibrational spectra are calculated by the so-called ”frozen phonon” method. In all cases, the DFT curves show an acceptable agreement with experimental results of inelastic neutron scattering. In case of the model glass former B2O3, a new classical interaction potential is parametrized, based on the liquid trajectory of an ab initio MD simulation at 2300 K. In this course, a structural fitting routine is used. The inclusion of 3-body angular interactions leads to a significantly improved agreement of the liquid properties of the classical MD and ab initio MD simulations. However, the generated glass structures, in all cases, show a significantly lower fraction of 3-membered planar boroxol rings as predicted by experimental results (f=60%-80%). The largest boroxol ring fraction of f=15±5% is observed in the full ab initio quenches from 2300 K. In case of SiO2, the glass structures after the quantum mechanical relaxation are the basis for calculations of the linear thermal expansion coefficient αL(T), employing the quasi-harmonic approximation. The striking observation is a change change of sign of αL(T) going along with a temperature range of negative αL(T) at low temperatures, which is in good agreement with experimental results.

Parameterizing a Simple Potential Energy Function Appropriate for Vibrational Frequencies Using Molecular Dynamics

A Maple scheme for quickly parameterizing vibrational potential energy functions is presented. As an example, the potential energy function's parameters for the vibrational motions in H_2O_2 are obtained assuming the simplest potential energy function. This paper was originally written as a research paper, but rejected by the referees. It is therefore being edited into an ``educational'' paper for student usage.

Vibrational population relaxation of carbon monoxide in the heme pocket of photolyzed carbonmonoxy myoglobin: Comparison of time-resolved mid-IR absorbance experiments and molecular dynamics simulations

The vibrational energy relaxation of carbon monoxide in the heme pocket of sperm whale myoglobin was studied by using molecular dynamics simulation and normal mode analysis methods. Molecular dynamics trajectories of solvated myoglobin were run at 300 K for both the δ- and ɛ-tautomers of the distal His-64. Vibrational population relaxation times of 335 ± 115 ps for the δ-tautomer and 640 ± 185 ps for the ɛ-tautomer were estimated by using the Landau–Teller model. Normal mode analysis was used to identify those protein residues that act as the primary “doorway” modes in the vibrational relaxation of the oscillator. Although the CO relaxation rates in both the ɛ- and δ-tautomers are similar in magnitude, the simulations predict that the vibrational relaxation of the CO is faster in the δ-tautomer with the distal His playing an important role in the energy relaxation mechanism. Time-resolved mid-IR absorbance measurements were performed on photolyzed carbonmonoxy hemoglobin (Hb13CO). From these measurements, a T1 time of 600 ± 150 ps was determined. The simulation and experimental estimates are compared and discussed.

Real and imaginary parts of the vibrational relaxation of acetonitrile in its electrolyte solutions:new results for the dynamics of solvent molecules

Isotropic scattering Raman spectra of liquid acetonitrile (AN) solutions of LiBF4 and NaI at various temperatures and concentrations have been investigated. For the first time imaginary as well as real parts of the solvent vibrational correlation functions have been extracted from the spectra. Such imaginary parts are currently an important component of modern theories of vibrational relaxation in liquids. This investigation thus provides the first experimental data on imaginary parts of a correlation function in AN solutions. Using the fitting algorithm we recently developed, statistically confident models for the Raman spectra were deduced. The parameters of the band shapes, with an additional correction, of the ν2 AN vibration (CN stretching), together with their confidence intervals are also reported for the first time. It is shown that three distinct species, with lifetimes greater than ∼10−13 s, of the AN molecules can be detected in solutions containing Li+ and Na+. These species are attributed to AN molecules directly solvating cations; the single oriented and polarised molecules interleaving the cation and anion of a Solvent Shared Ion Pair (SShIP); and molecules solvating anions. These last are considered to be equivalent to the next layer of solvent molecules, because the CN end of the molecule is distant from the anion and thus less affected by the ionic charge compared with the anion situation. Calculations showed that at the concentrations employed, 1 and 0.3 M, there were essentially no other solvent molecules remaining that could be considered as bulk solvent. Calculations also showed that the internuclear distance in these solutions supported the proposal that the ionic entity dominating in solution was the SShIP, and other evidence was adduced that confirmed the absence of Contact Ion Pairs at these concentrations. The parameters of the shape of the vibrational correlation functions of all three species are reported. The parameters of intramolecular anharmonic coupling between the potential surfaces in AN and the dynamics of the intermolecular environment fluctuations and intermolecular energy transfer are presented. These results will assist investigations made at higher and lower concentrations, when additional species and interactions with AN molecules will be present.

Quantum signatures of chaos in the dynamics of a trapped ion

We show how a nonlinear chaotic system, the parametrically kicked nonlinear oscillator, may be realized in the dynamics of a trapped, laser-cooled ion, interacting with a sequence of standing-wave pulses. Unlike the original optical scheme [G. J. Milburn and C.A. Holmes, Phys. Rev. A 44, 4704 (1991)], the trapped ion enables strongly quantum dynamics with minimal dissipation. This should permit an experimental test of one of the quantum signatures of chaos: irregular collapse and revival dynamics of the average vibrational energy.

Dinámica y cinética de los procesos físicos y químicos del radical OH con compuestos orgánicos volátiles (COVs) de relevancia atmosférica. Dynamics and kinetics of the physical and chemical processes of OH radicals with atmospheric relevant volatile organic compounds (VOCs)

El objetivo del presente proyecto es estudiar los procesos físicos y químicos del radical OH con compuestos orgánicos volátiles (COVs), con los cuales sea factible la formación de agregados de van der Waals (vdW) responsables de la curvatura en los gráficos de Arrhenius, empleando técnicas modernas, complementarias entre si y novedosas en el país. El problema será abordado desde tres perspectivas complementarias: 1) estudios cinéticos, 2) estudios mecanísticos y de distribución de productos y 3) estudios de la dinámica de los procesos físicos y químicos. La finalidad es alcanzar una mejor comprensión de los mecanismos que intervienen en el comportamiento químico de especies presentes en la atmósfera y obtener datos cinéticos de alta calidad que puedan alimentar modelos computacionales capaces de describir la composición de la atmósfera, presente y futura. Los objetivos son estudiar: 1) mediante fotólisis láser pulsada con detección por fluorescencia inducida por láser (PLP-LIF), en reactores de flujo, la cinética de reacción del radical OH(v”=0) con COVs que presentan gráficos de Arrhenius curvos con energías de activación negativas, tales como alcoholes insaturados, alquenos halogenados, éteres halogenados, ésteres alifáticos; 2) en una cámara de simulación de condiciones atmosféricas de gran volumen (4500 L), la identidad y el rendimiento de productos de las reacciones mencionadas, a fines de evaluar su impacto atmosférico y dilucidar los mecanismos de reacción; 3) mediante haces moleculares y espectroscopía láser, la estructura y reactividad de complejos de vdW entre alcoholes insaturados o aromáticos (cresoles) y el radical OH, como modelo de los aductos propuestos como responsables de la desviación al comportamiento de Arrhenius de las reacciones mencionadas; 4) mediante PLP-LIF y expansiones supersónicas, las constantes específicas estado a estado (ksts) de relajación/reacción del radical OH(v”=1-4) vibracionalmente excitado con los COVs mencionados. Los resultados experimentales obtenidos serán contrastados con cálculos ab-initio de estructura electrónica, los cuales apoyarán las interpretaciones, permitirán proponer estructuras de estados de transición y aductos colisionales, como así también calcular las frecuencias de vibración de los complejos de vdW para su posterior asignación en los espectros LIF y REMPI. Asimismo, los mecanismos de reacción propuestos y los parámetros cinéticos medidos experimentalmente serán comparados con aquellos obtenidos por cálculos teóricos. The aim of this project is to study the physical and chemical processes of OH radicals with volatile organic compounds (VOCs) with which the formation of van der Waals (vdW) clusters, responsible for the observed curvature in the Arrhenius plots, might be feasible. The problem will be addressed as follow : 1) kinetic studies; 2) products distribution and mechanistic studies and 3) dynamical studies of the physical and chemical processes. The purpose is to obtain a better understanding of the mechanisms that govern the chemical behavior of species present in the atmosphere and to obtain high quality kinetic data to be used as input to computational models. We will study: 1) the reaction kinetics of OH (v”=0) radicals with VOCs such as unsaturated alcohols, halogenated alkenes, halogenated ethers, aliphatic esters, which show curved Arrhenius plots and negative activation energies, by PLP-LIF, in flow systems; 2) in a large volume (4500 L) atmospheric simulation chamber, reaction products yields in order to evaluate their atmospheric impact and reaction mechanisms; 3) using molecular beams and laser spectroscopy, the structure and reactivity of the vdW complexes formed between the unsaturated or aromatic alcohols and the OH radicals as a model of the adducts proposed as responsible for the non-Arrhenius behavior; 4) the specific state-to-state relaxation/reaction rate constants (ksts) of the vibrationally excited OH (v”=1-4) radical with the VOCs by PLP-LIF and supersonic expansions. Ab-initio calculations will be carried out to support the interpretation of the experimental results, to obtain the transition state and collisional adducts structures, as well as to calculate the vibrational frequencies of the vdW complexes to assign to the LIF and REMPI spectra. Also, the proposed reaction mechanisms and the experimentally measured kinetic parameters will be compared with those obtained from theoretical calculations.

Treatment of nonlinear optical properties due to large amplitude anharmonic vibrational motions: umbrella motion in NH3

A general reduced dimensionality finite field nuclear relaxation method for calculating vibrational nonlinear optical properties of molecules with large contributions due to anharmonic motions is introduced. In an initial application to the umbrella (inversion) motion of NH3 it is found that difficulties associated with a conventional single well treatment are overcome and that the particular definition of the inversion coordinate is not important. Future applications are described

Additional compact formulas for vibrational dynamic dipole polarizabilities and hyperpolarizabilities

Compact expressions, complete through second order in electrical and/or mechanical anharmonicity, are given for the dynamic dipole vibrational polarizability and dynamic first and second vibrational hyperpolarizabilities. Certain contributions not previously formulated are now included

Anharmonicity contributions to the vibrational second hyperpolarizability of conjugated oligomers

Restricted Hartree-Fock 6-31G calculations of electrical and mechanical anharmonicity contributions to the longitudinal vibrational second hyperpolarizability have been carried out for eight homologous series of conjugated oligomers - polyacetylene, polyyne, polydiacetylene, polybutatriene, polycumulene, polysilane, polymethineimine, and polypyrrole. To draw conclusions about the limiting infinite polymer behavior, chains containing up to 12 heavy atoms along the conjugated backbone were considered. In general, the vibrational hyperpolarizabilities are substantial in comparison with their static electronic counterparts for the dc-Kerr and degenerate four-wave mixing processes (as well as for static fields) but not for electric field-induced second harmonic generation or third harmonic generation. Anharmonicity terms due to nuclear relaxation are important for the dc-Kerr effect (and for the static hyperpolarizability) in the σ-conjugated polymer, polysilane, as well as the nonplanar π systems polymethineimine and polypyrrole. Restricting polypyrrole to be planar, as it is in the crystal phase, causes these anharmonic terms to become negligible. When the same restriction is applied to polymethineimine the effect is reduced but remains quantitatively significant due to the first-order contribution. We conclude that anharmonicity associated with nuclear relaxation can be ignored, for semiquantitative purposes, in planar π-conjugated polymers. The role of zero-point vibrational averaging remains to be evaluated

Linear response functions for a vibrational configuration interaction state

Linear response functions are implemented for a vibrational configuration interaction state allowing accurate analytical calculations of pure vibrational contributions to dynamical polarizabilities. Sample calculations are presented for the pure vibrational contributions to the polarizabilities of water and formaldehyde. We discuss the convergence of the results with respect to various details of the vibrational wave function description as well as the potential and property surfaces. We also analyze the frequency dependence of the linear response function and the effect of accounting phenomenologically for the finite lifetime of the excited vibrational states. Finally, we compare the analytical response approach to a sum-over-states approach