Understanding solvent effects on structure and reactivity of organic intermediates: a Raman study

The effect of solvent on chemical reactivity has generally been explained on the basis of the dielectric constant and viscosity. However a number of spectroscopic studies, including UV-VIS, IR and Raman, has led to numerous empirical parameters to define solvent effect based on either solvating ability or polarity scale. These parameters include solvent polarizability, dipolarity, Lewis acidity and Lewis basicity, E-T(30), pi*, alpha, beta etc. However, from a structural point of view, we can separate solvation as static and dynamic processes. The static solvation basically relates to stabilization of the molecular structure by the solvent to attain the equilibrium structure, both in the intermediate and ground state. Dynamic solvation relates to solvent reorganization-induced dynamics prior to the structural reorganization to reach the equilibrium state. In this paper, we present (a) structural distortions induced by the solvent due to preferential solvation of the triplet excited state, and (b) the importance of dynamic solvation induced by vibronic coupling (pseudo-Jahn-Teller coupling). The examples include the effect of solvent on structure and reactivity of excited states of 2,2,2-trifluoroacetophenone (TFA). Based on the comparison of time resolved resonance Raman (TR3) data of TFA and other substituted acetophenone systems, it was found that change in solvent polarity indeed results in electronic state switching and structural changes in the excited state, which explains the trend in reactivity. Further, a TR3 study of fluoranil (FA) in the triplet excited state in solvents of varying polarities indicates that the structure of FA in the triplet excited state is determined by vibronic coupling effects and thus distorted structure. These experimental results have been well supported by density functional theoretical computational studies.

Free Energy Gap Dependence of the Electron-Transfer Rate from the Inverted to the Normal Region

numerical study of the free energy gap (FEG) dependence of the electron-transfer rate in polar solvents is presented. This study is based on the generalized multidimensional hybrid model, which not only includes the solvent polarization and the molecular vibration modes, but also the biphasic polar response of the solvent. The free energy gap dependence is found to be sensitive to several factors, including the solvent relaxation rate, the electronic coupling between the surfaces, the frequency of the high-frequency quantum vibrational mode, and the magnitude of the solvent reorganization energy. It is shown that in some cases solvent relaxation can play an important role even in the Marcus normal regime. The minimal hybrid model involves a large number of parameters, giving rise to a diverse non-Marcus FEG behavior which is often determined collectively by these parameters. The model gives the linear free energy gap dependence of the logarithmic rate over a substantial range of FEG, spanning from the normal to the inverted regime. However, even for favorable values of the relevant parameters, a linear free energy gap dependence of the rate could be obtained only over a range of 5000-6000 cm(-1) (compared to the experimentally observed range of 10000 cm(-1) reported by Benniston et al.). The present work suggests several extensions/generalizations of the hybrid model which might be necessary to fully understand the observed free energy gap dependence.

On the relationship of thermodynamic parameters with the buried surface area in protein-ligand complex formation

Prediction of thermodynamic parameters of protein-protein and antigen-antibody complex formation from high resolution structural parameters has recently received much attention, since an understanding of the contributions of different fundamental processes like hydrophobic interactions, hydrogen bonding, salt bridge formation, solvent reorganization etc. to the overall thermodynamic parameters and their relations with the structural parameters would lead to rational drug design. Using the results of the dissolution of hydrocarbons and other model compounds the changes in heat capacity (DeltaCp), enthalpy (DeltaH) and entropy (DeltaS) have been empirically correlated with the polar and apolar surface areas buried during the process of protein folding/unfolding and protein-ligand complex formation. In this regard, the polar and apolar surfaces removed from the solvent in a protein-ligand complex have been calculated from the experimentally observed values of changes in heat capacity (DeltaCp) and enthalpy (DeltaH) for protein-ligand complexes for which accurate thermodynamic and high resolution structural data are available, and the results have been compared with the x-ray crystallographic observations. Analyses of the available results show poor correlation between the thermodynamic and structural parameters. Probable reasons for this discrepancy are mostly related with the reorganization of water accompanying the reaction which is indeed proven by the analyses of the energetics of the binding of the wheat germ agglutinin to oligosaccharides.

Development and Application of Embedding Methods for the Simulation of Large Chemical Systems

The high computational cost of correlated wavefunction theory (WFT) calculations has motivated the development of numerous methods to partition the description of large chemical systems into smaller subsystem calculations. For example, WFT-in-DFT embedding methods facilitate the partitioning of a system into two subsystems: a subsystem A that is treated using an accurate WFT method, and a subsystem B that is treated using a more efficient Kohn-Sham density functional theory (KS-DFT) method. Representation of the interactions between subsystems is non-trivial, and often requires the use of approximate kinetic energy functionals or computationally challenging optimized effective potential calculations; however, it has recently been shown that these challenges can be eliminated through the use of a projection operator. This dissertation describes the development and application of embedding methods that enable accurate and efficient calculation of the properties of large chemical systems.

Chapter 1 introduces a method for efficiently performing projection-based WFT-in-DFT embedding calculations on large systems. This is accomplished by using a truncated basis set representation of the subsystem A wavefunction. We show that naive truncation of the basis set associated with subsystem A can lead to large numerical artifacts, and present an approach for systematically controlling these artifacts.

Chapter 2 describes the application of the projection-based embedding method to investigate the oxidative stability of lithium-ion batteries. We study the oxidation potentials of mixtures of ethylene carbonate (EC) and dimethyl carbonate (DMC) by using the projection-based embedding method to calculate the vertical ionization energy (IE) of individual molecules at the CCSD(T) level of theory, while explicitly accounting for the solvent using DFT. Interestingly, we reveal that large contributions to the solvation properties of DMC originate from quadrupolar interactions, resulting in a much larger solvent reorganization energy than that predicted using simple dielectric continuum models. Demonstration that the solvation properties of EC and DMC are governed by fundamentally different intermolecular interactions provides insight into key aspects of lithium-ion batteries, with relevance to electrolyte decomposition processes, solid-electrolyte interphase formation, and the local solvation environment of lithium cations.

Atomistic simulations of cosolvent/water mixtures

In this thesis, atomistic simulations are performed to investigate hydrophobic solvation and hydrophobic interactions in cosolvent/water binary mixtures. Many cosolvent/water binary mixtures exhibit non-ideal behavior caused by aggregation at the molecular scale level although they are stable and homogenous at the macroscopic scale. Force-field based atomistic simulations provide routes to relate atomistic-scale structure and interactions to thermodynamic solution properties. The predicted solution properties are however sensitive to the parameters used to describe the molecular interactions. In this thesis, a force field for tertiary butanol (TBA) and water mixtures is parameterized by making use of the Kirkwood-Buff theory of solution. The new force field is capable of describing the alcohol-alcohol, water-water and alcohol-water clustering in the solution as well as the solution components’ chemical potential derivatives in agreement with experimental data. With the new force field, the preferential solvation and the solvation thermodynamics of a hydrophobic solute in TBA/water mixtures have been studied. First, methane solvation at various TBA/water concentrations is discussed in terms of solvation free energy-, enthalpy- and entropy- changes, which have been compared to experimental data. We observed that the methane solvation free energy varies smoothly with the alcohol/water composition while the solvation enthalpies and entropies vary nonmonotonically. The latter occurs due to structural solvent reorganization contributions which are not present in the free energy change due to exact enthalpy-entropy compensation. It is therefore concluded that the enthalpy and entropy of solvation provide more detailed information on the reorganization of solvent molecules around the inserted solute. Hydrophobic interactions in binary urea/water mixtures are next discussed. This system is particularly relevant in biology (protein folding/unfolding), however, changes in the hydrophobic interaction induced by urea molecules are not well understood. In this thesis, this interaction has been studied by calculating the free energy (potential of mean force), enthalpy and entropy changes as a function of the solute-solute distance in water and in aqueous urea (6.9 M) solution. In chapter 5, the potential of mean force in both solution systems is analyzed in terms of its enthalpic and entropic contributions. In particular, contributions of solvent reorganization in the enthalpy and entropy changes are studied separately to better understand what are the changes in interactions in the system that contribute to the free energy of association of the nonpolar solutes. We observe that in aqueous urea the association between nonpolar solutes remains thermodynamically favorable (i.e., as it is the case in pure water). This observation contrasts a long-standing belief that clusters of nonpolar molecules dissolve completely in the presence of urea molecules. The consequences of our observations for the stability of proteins in concentrated urea solutions are discussed in the chapter 6 of the thesis.

Sequence-specific DNA-binding dominated by dehydration.

Fluorescence spectroscopy and isothermal titration calorimetry were used to study the thermodynamics of binding of the glucocorticoid receptor DNA-binding domain to four different, but similar, DNA-binding sites. The binding sites are two naturally occurring sites that differ in the composition of one base pair, i.e., an A-T to G-C mutation, and two sites containing chemical intermediates of these base pairs. The calorimetrically determined heat capacity change (Delta C(p)o(obs)) for glucocorticoid receptor DNA-binding domain binding agrees with that calculated for dehydration of solvent-accessible surface areas. A dominating effect of dehydration or solvent reorganization on the thermodynamics is also consistent with an observed linear relationship between observed enthalpy change (Delta Ho(obs)) and observed entropy change (Delta So(obs)) with a slope close to the experimental temperature. Comparisons with structural data allow us to rationalize individual differences between Delta Ho(obs) (and Delta So(obs)) for the four complexes. For instance, we find that the removal of a methyl group at the DNA-protein interface is enthalpically favorable but entropically unfavorable, which is consistent with a replacement by an ordered water molecule.

Pressure and temperature effects on metal-to-metal charge transfer in cyano-bridged Co-III-Fe-II complexes

The effects of pressure and temperature on the energy (E-op) of the metal-to-metal charge transfer (MMCT, Fe-II --> Co-III) transition of the cyano-bridged complexes trans - [(LCoNCFe)-Co-14(CN)(5)](-) and cis-[(LCoNCFe)-Co-14(CN)(5)](-) (where L-14 = 6-methyl-1,4,8,11-tetraazacyclotetradecan-6-amine) were examined. The changes in the redox potentials of the cobalt and iron metal centres with pressure and temperature were also examined and the results interpreted with Marcus Hush theory. The observed redox reaction volumes can mainly be accounted for in terms of localised electrostriction effects. The shifts in E-op due to both pressure and temperature were found to be less than the shifts in the energy difference (E degrees) between the Co-III-Fe-II and Co-II-Fe-III redox isomers. The pressure and temperature dependence of the reorganisational energy, as well as contributions arising from the different spin states of Co-II, are discussed in order to account for this trend. To study the effect of pressure on Co-III electronic absorption bands, a new cyano-bridged complex, trans - [(LCoNCCo)-Co-14(CN)(5)], was prepared and characterised spectroscopically and structurally. X-Ray crystallography revealed this complex to be isostructural with trans -[(LCoNCFe)-Co-14(CN)(5)] center dot 5H(2)O.

Microphase Separation of Mixed Polymer Brushes: Dependence of the Morphology on Grafting Density, Composition, Chain-Length Asymmetry, Solvent Quality, and Selectivity

Microphase separation of binary mixed A/B polymer brushes exposed to different solvents is studied using Single-Chain-in-Mean-Field simulations. Effects of solvent quality and selectivity, grafting density, composition, and chain-length asymmetry are systematically investigated, and diagrams of morphologies in various solvents are constructed as a function of grafting density and composition or chain-length asymmetry. The structure of the microphase segregated morphologies lacks long-range periodic order, and it is analyzed quantitatively Using Minkowski measures.

A nonmonotonic dependence of standard rate constant on reorganization energy for heterogeneous electron transfer processes on electrode surface

In the present work a nonmonotonic dependence of standard rate constant (k(0)) on reorganization energy (lambda) was discovered qualitatively from electron transfer (Marcus-Hush-Levich) theory for heterogeneous electron transfer processes on electrode surface. It was found that the nonmonotonic dependence of k(0) on lambda is another result, besides the disappearance of the famous Marcus inverted region, coming from the continuum of electronic states in electrode: with the increase of lambda, the states for both Process I and Process II ET processes all vary from nonadiabatic to adiabatic state continuously, and the lambda dependence of k(0) for Process I is monotonic thoroughly, while for Process II on electrode surface the lambda dependence of k(0) could show a nonmonotonicity.

Solvent-dependent modulation of metal–metal electronic interactions in a dinuclear cyanoruthenate complex: a detailed electrochemical, spectroscopic and computational study

The dinuclear complex [{Ru(CN)4}2(μ-bppz)]4− shows a strongly solvent-dependent metal–metal electronic interaction which allows the mixed-valence state to be switched from class 2 to class 3 by changing solvent from water to CH2Cl2. In CH2Cl2 the separation between the successive Ru(II)/Ru(III) redox couples is 350 mVand the IVCT band (from the UV/Vis/NIR spectroelectrochemistry) is characteristic of a borderline class II/III or class III mixed valence state. In water, the redox separation is only 110 mVand the much broader IVCT transition is characteristic of a class II mixed-valence state. This is consistent with the observation that raising and lowering the energy of the d(π) orbitals in CH2Cl2 or water, respectively, will decrease or increase the energy gap to the LUMO of the bppz bridging ligand, which provides the delocalisation pathway via electron-transfer. IR spectroelectrochemistry could only be carried out successfully in CH2Cl2 and revealed class III mixed-valence behaviour on the fast IR timescale. In contrast to this, time-resolved IR spectroscopy showed that the MLCTexcited state, which is formulated as RuIII(bppz˙−)RuII and can therefore be considered as a mixed-valence Ru(II)/Ru(III) complex with an intermediate bridging radical anion ligand, is localised on the IR timescale with spectroscopically distinct Ru(II) and Ru(III) termini. This is because the necessary electron-transfer via the bppz ligand is more difficult because of the additional electron on bppz˙− which raises the orbital through which electron exchange occurs in energy. DFT calculations reproduce the electronic spectra of the complex in all three Ru(II)/Ru(II), Ru(II)/Ru(III) and Ru(III)/Ru(III) calculations in both water and CH2Cl2 well as long as an explicit allowance is made for the presence of water molecules hydrogen-bonded to the cyanides in the model used. They also reproduce the excited-state IR spectra of both [Ru(CN)4(μ-bppz)]2– and [{Ru(CN)4}2(μ-bppz)]4− very well in both solvents. The reorganization of the water solvent shell indicates a possible dynamical reason for the longer life time of the triplet state in water compared to CH2Cl2.

Polyurethane (PU) scaffolds prepared by solvent casting/particulate leaching (SCPL) combined with centrifugation

This article reports an enhanced solvent casting/particulate (salt) leaching (SCPL) method developed for preparing three-dimensional porous polyurethane (PU) scaffolds for cardiac tissue engineering. The solvent for the preparation of the PU scaffolds was a mixture of dimethylformamide (DFM) and tetrahydrofuran (THF). The enhanced method involved the combination of a conventional SCPL method and a step of centrifugation, with the centrifugation being employed to improve the pore uniformity and the pore interconnectivity of scaffolds. Highly porous three-dimensional scaffolds with a well interconnected porous structure could be achieved at the polymer solution concentration of up to 20% by air or vacuum drying to remove the solvent. When the salt particle sizes of 212-295, 295-425, or 425-531 µm and a 15% w/v polymer solution concentration were used, the porosity of the scaffolds was between 83-92% and the compression moduli of the scaffolds were between 13 kPa and 28 kPa. Type I collagen acidic solution was introduced into the pores of a PU scaffold to coat the collagen onto the pore walls throughout the whole PU scaffold. The human aortic endothelial cells (HAECs) cultured in the collagen-coated PU scaffold for 2 weeks were observed by scanning electron microscopy (SEM). It was shown that the enhanced SCPL method and the collagen coating resulted in a spatially uniform distribution of cells throughout the collagen-coated PU scaffold.

Use of micro-ATR/FTIR imaging to study heterogeneous polymer oxidation by direct solvent casting onto the ATR IRE

A technique is described whereby micro-ATR/FTIR imaging can be used to follow polymer degradation reactions in situ in real time. The internal reflection element (IRE) assembly is removed from the ATR objective and polymer is solvent cast directly onto the IRE surface. The polymer is then subjected to degradation conditions and molecular structural changes monitored by periodically replacing the IRE assembly back in the ATR objective and collecting spectra which can be used to construct images. This approach has the benefit that the same part of the sample is always studied, and that contact by pressure which might damage the polymer surface is not required. The technique is demonstrated using the polymer Topas which was degraded by exposure to UVC light in air.

Asymptotic and numerical results for a model of solvent-dependent drug diffusion through polymeric spheres

A model for drug diffusion from a spherical polymeric drug delivery device is considered. The model contains two key features. The first is that solvent diffuses into the polymer, which then transitions from a glassy to a rubbery state. The interface between the two states of polymer is modelled as a moving boundary, whose speed is governed by a kinetic law; the same moving boundary problem arises in the one-phase limit of a Stefan problem with kinetic undercooling. The second feature is that drug diffuses only through the rubbery region, with a nonlinear diffusion coefficient that depends on the concentration of solvent. We analyse the model using both formal asymptotics and numerical computation, the latter by applying a front-fixing scheme with a finite volume method. Previous results are extended and comparisons are made with linear models that work well under certain parameter regimes. Finally, a model for a multi-layered drug delivery device is suggested, which allows for more flexible control of drug release.