Thermodynamics of Solution of Some Alkyl Acetates in Water

Summary.  The temperature dependence of the solubility of methyl acetate (MeOAc), ethyl acetate (EtOAc), 1-propyl acetate (1-PrOAc), 1-butyl acetate (1-BuOAc), 2-methyl-1-propyl acetate (iso-BuOAc), 2-butyl acetate (sec-BuOAc), and 2-methyl-2-propyl acetate (tert-BuOAc) in water in the temperature range from 298.2 to 318.2 K was determined. The experimental solubility data together with some literature values are presented as a function of temperature and given in analytical form. The solubilities of the investigated compounds at 298.2 K were correlated with the number of carbon atoms in the solute molecule. The standard thermodynamic functions for solubility, i.e. Gibbs free energy, enthalpy, heat capacity, and entropy, were calculated. The Gibbs free energy is positive and increases linearly with the number of carbon atoms, whereas the enthalpy and entropy are negative. The standard thermodynamic functions were converted to thermodynamic quantities for hydration. The Gibbs free energies of hydration of alkyl acetates, which are negative, are compared with the Gibbs free energy of hydration of some n-alkanes, 1-alcohols, and 1-alkylamines. The standard thermodynamic functions of hydration were analyzed using a modified version of the theoretical approach developed by Lee and Graziano [1]. Received October 16, 2000. Accepted May 14, 2001     

Monte Carlo simulations of the hydration of substituted benzenes with OPLS potential functions

Intermolecular potential functions have been developed for use in computer simulations of substituted benzenes. Previously reported optimized potentials for liquid simulations (OPLS) for benzene and organic functional groups were merged and tested by computing free energies of hydration for toluene, p-xylene, phenol, anisole, benzonitrile, p-cresol, hydroquinone, and p-dicyanobenzene. The calculations featured Monte Carlo simulations at 25°C and 1 atm with statistical perturbation theory. The average difference between the computed results and experimental data for the absolute free energies of hydration is 0.5 kcal/mol. The AM1-SM2 method is also found to perform well in predicting the free energies of hydration for the substituted benzenes. In addition, the Monte Carlo simulations provided details on the hydration of the substituted benzenes, in particular for the solute–water hydrogen bonding. © 1993 John Wiley & Sons, Inc.     

Characterization of the interaction energies for polystyrene blends with various methacrylate polymers

The binary interaction energies between styrene and various methacrylates were determined from newly examined phase boundaries with lattice–fluid theory. Because the blends of polystyrene (PS) and poly(cyclohexylmethacrylate) (PCHMA) were only miscible at high molecular weights when the blends were prepared by solution casting from tetrahydrofuran, we examined the miscibility of other blends by changing the molecular weights of PS or methacrylate polymers. On the basis of the phase‐separation temperature caused by the lower critical solution temperature, the miscibility of PS with the various methacrylates appeared to be in the order PCHMA > poly(n‐propyl‐methacrylate) (PnPMA) > poly(ethyl methacrylate) (PEMA) > poly(n‐butyl‐methacrylate) (PnBMA) > poly(iso‐butyl‐methacrylate) > poly(methyl methacrylate) (PMMA) > poly(tert‐butyl methacrylate), and the branching of butylmethacrylate appeared to decrease the miscibility with PS. The interaction energies between PS with various methacrylates obtained from phase boundaries with lattice–fluid theory reached minimum value corresponding to the styrene/n‐propylmethacrylate interaction. They were in the order PnPMA < PEMA < PCHMA < PnBMA < PMMA. The difference in the order of miscibility and interaction energies might be attributed to the terms related to the compressibility. The phase‐separation temperatures calculated with the interaction energies obtained here indicated that the PS/PEMA and PS/PnPMA blends at high molecular weights were miscible, whereas the PS/PnBMA blends were immiscible at high molecular weights. © 2000 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 38: 2666–2677, 2000     

Density functional theory studies on the dissociation energies of metallic salts: relationship between lattice and dissociation energies

The formation and physicochemical properties of polymer electrolytes strongly depend on the lattice energy of metal salts. An indirect but efficient way to estimate the lattice energy through the relationship between the heterolytic bond dissociation and lattice energies is proposed in this work. The heterolytic bond dissociation energies for alkali metal compounds were calculated theoretically using the Density Functional Theory (DFT) of B3LYP level with 6‐311+G(d,p) and 6‐311+G(2df,p) basis sets. For transition metal compounds, the same method was employed except for using the effective core potential (ECP) of LANL2DZ and SDD on transition metals for 6‐311+G(d,p) and 6‐311+G(2df,p) calculations, respectively. The dissociation energies calculated by 6‐311+G(2df,p) basis set combined with SDD basis set were better correlated with the experimental values with average error of ca. ±1.0% than those by 6‐311+G* combined with the LANL2DZ basis set. The relationship between dissociation and lattice energies was found to be fairly linear (r>0.98). Thus, this method can be used to estimate the lattice energy of an unknown ionic compound with reasonably high accuracy. We also found that the dissociation energies of transition metal salts were relatively larger than those of alkaline metal salts for comparable ionic radii. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 827–834, 2001     

Rapid approximate calculation of water binding free energies in the whole hydration domain of (bio)macromolecules

The evaluation of water binding free energies around solute molecules is important for the thermodynamic characterization of hydration or association processes. Here, a rapid approximate method to estimate water binding free energies around (bio)macromolecules from a single molecular dynamics simulation is presented. The basic idea is that endpoint free‐energy calculation methods are applied and the endpoint quantities are monitored on a three‐dimensional grid around the solute. Thus, a gridded map of water binding free energies around the solute is obtained, that is, from a single short simulation, a map of favorable and unfavorable water binding sites can be constructed. Among the employed free‐energy calculation methods, approaches involving endpoint information pertaining to actual thermodynamic integration calculations or endpoint information as exploited in the linear interaction energy method were examined. The accuracy of the approximate approaches was evaluated on the hydration of a cage‐like molecule representing either a nonpolar, polar, or charged water binding site and on α‐ and β‐cyclodextrin molecules. Among the tested approaches, the linear interaction energy method is considered the most viable approach. Applying the linear interaction energy method on the grid around the solute, a semi‐quantitative thermodynamic characterization of hydration around the whole solute is obtained. Disadvantages are the approximate nature of the method and a limited flexibility of the solute. © 2016 Wiley Periodicals, Inc.     

Experimentation with different thermodynamic cycles used for pKa calculations on carboxylic acids using complete basis set and Gaussian‐n models combined with CPCM continuum solvation methods

Complete basis set and Gaussian‐n methods were combined with Barone and Cossi's implementation of the polarizable conductor model (CPCM) continuum solvation methods to calculate pK_a values for six carboxylic acids. Four different thermodynamic cycles were considered in this work. An experimental value of ?264.61 kcal/mol for the free energy of solvation of H⁺, ΔG_s(H⁺), was combined with a value for G_gas(H⁺) of ?6.28 kcal/mol, to calculate pK_a values with cycle 1. The complete basis set gas‐phase methods used to calculate gas‐phase free energies are very accurate, with mean unsigned errors of 0.3 kcal/mol and standard deviations of 0.4 kcal/mol. The CPCM solvation calculations used to calculate condensed‐phase free energies are slightly less accurate than the gas‐phase models, and the best method has a mean unsigned error and standard deviation of 0.4 and 0.5 kcal/mol, respectively. Thermodynamic cycles that include an explicit water in the cycle are not accurate when the free energy of solvation of a water molecule is used, but appear to become accurate when the experimental free energy of vaporization of water is used. This apparent improvement is an artifact of the standard state used in the calculation. Geometry relaxation in solution does not improve the results when using these later cycles. The use of cycle 1 and the complete basis set models combined with the CPCM solvation methods yielded pK_a values accurate to less than half a pK_a unit. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Substituent effect on the cationic ring‐opening polymerization of acyloxymethyl five‐member cyclic dithiocarbonates

Five‐member cyclic dithiocarbonates were synthesized by the reactions of carbon disulfide with benzoic, p‐anisic, p‐chlorobenzoic, 1‐naphthalenecarboxylic, p‐nitrobenzoic, and p‐(tert‐butyl)benzoic glycidyl esters, and their cationic ring‐opening polymerizations were carried out with methyl trifluoromethane sulfonate and trifluoromethane sulfonic acid as initiators at room temperature to 80 °C. Polymers with number‐average molecular weights of 3400–24,900 were obtained in high yields, and their structures were estimated by NMR and IR spectroscopy. The monomers showed a clear difference in the polymerization rate according to the substituents. The rate of polymerization decreased in the order of p‐chlorobenzoic ≥ benzoic > 1‐naphthalenecarboxylic > p‐nitro‐benzoic > p‐tert‐butylbenzoic > p‐anisic. The data of the reaction kinetics, NMR studies, and molecular orbital calculations proved a plausible mechanism involving the participation of p‐substituted benzoyloxymethyl groups to stabilize the cationic propagating end. The polymers showed decomposition temperatures with 5% weight loss ranging from 200 to 260 °C. No glass‐transition temperatures for the polymers were observed below 200 °C by differential scanning calorimetry. © 2001 John Wiley & Sons, Inc. J Polym Sci Part A: Polym Chem 39: 3967–3980, 2001     

DFT plus PCM calculation for pairing specificity of Watson–Crick‐type bases in aqueous solutions

This study is aimed at explaining the preference for AT and CG pairings and the possible insertion of other tautomeric DNA base pairs such as G_enolT, that respect energetic and steric requirements including at least two hydrogen bonds and 11 ± 0.5Å distance between the 9‐CH₃ of purine and 5‐CH₃ of pyrimidine. The calculated free energy of formation ΔΔG at the DFT B3LYP/6‐31G*‐PCM/BEM level pointed out the CG and AT pairs as the most favored, followed closely by G_enolT, in good agreement with Michaelis–Menten first order kinetics (CG ≈ AT > G_enolT). Unusual DNA base pairs complexes such as AG (BEM) and CT (PCM) resulted to be stable, but it is very difficult to assume that they are likely to be included in the double strand DNA. The calculated enthalpy and dipole moments of isolated DNA bases agree well with experiment. The free energy of hydration, ΔG_hyd, was found to depend on the electrostatic term, while cavitation‐dispersion components are almost constant. The stability of DNA complexes in water resulted from PCM calculations is markedly influenced by the free energy of hydration. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Examination of several density functionals in numerical Kohn–Sham calculations for atoms

We studied several exchange‐only and exchange–correlation energy density functionals in numerical, i.e., basis‐set‐free, nonrelativistic Kohn–Sham calculations for closed‐shell ¹S states of atoms and atomic ions with N electrons, where 2≤N≤120. Accurate total energies are presented to serve as reference data for algebraic approaches, as do the numerical Hartree–Fock results, which are also provided. Gradient‐corrected exchange‐only functionals considerably improve the total energies obtained from the usual local density approximation, when compared to the Hartree–Fock results. Such an improvement due to gradient corrections is not seen in general for highest orbital energies, neither for exchange‐only results (to be compared with Hartree–Fock results), nor for exchange–correlation results (to be compared with experimental ionization energies). © 2001 John Wiley & Sons, Inc. Int J Quant Chem 82: 227–241, 2001     

Thermodynamics of Solution of Some Alkyl Acetates in Water

 The temperature dependence of the solubility of methyl acetate (MeOAc), ethyl acetate (EtOAc), 1-propyl acetate (1-PrOAc), 1-butyl acetate (1-BuOAc), 2-methyl-1-propyl acetate (iso-BuOAc), 2-butyl acetate (sec-BuOAc), and 2-methyl-2-propyl acetate (tert-BuOAc) in water in the temperature range from 298.2 to 318.2 K was determined. The experimental solubility data together with some literature values are presented as a function of temperature and given in analytical form. The solubilities of the investigated compounds at 298.2 K were correlated with the number of carbon atoms in the solute molecule. The standard thermodynamic functions for solubility, i.e. Gibbs free energy, enthalpy, heat capacity, and entropy, were calculated. The Gibbs free energy is positive and increases linearly with the number of carbon atoms, whereas the enthalpy and entropy are negative. The standard thermodynamic functions were converted to thermodynamic quantities for hydration. The Gibbs free energies of hydration of alkyl acetates, which are negative, are compared with the Gibbs free energy of hydration of some n-alkanes, 1-alcohols, and 1-alkylamines. The standard thermodynamic functions of hydration were analyzed using a modified version of the theoretical approach developed by Lee and Graziano [1].     

Calculation of sequence‐dependent free energies of hydration of dipeptides formed by alanine and glycine

The relative free energies of hydration of the dipeptides glycylalanine and alanyl‐glycine in their naturally occurring form have been calculated both for the zwitterionic and protonated species. Emphasis was laid on comparisons between the conventional cutoff method and the Particle Mesh Ewald method to account for possible differences in electrostatic contributions to the free energy. Furthermore, the convergence behavior of the total free energy and its individual contributions were examined. The results, obtained by means of the thermodynamic integration technique as implemented in the free energy module of the AMBER program suite, suggest that in aqueous solution glycylalanine is more stable than alanylglycine by 2.7 kcal/mol in the zwitterionic form and by 3.5 kcal/mol in the protonated form. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 846–860, 2001     

Computation of hydration free energies using a parameterized continuum model: Study of equilibrium geometries and reactive processes in water solution

A parameterized self-consistent reaction field model allowing computation of the total free energy of hydration of organic molecules at the ab initio level is presented. The approach uses electrostatic plus polarization energies calculated with the help of a continuum model. The remaining solvation free energy terms are obtained by a simple formula based on atomic parameters and atomic accessible surface areas (ASAs), which are determined with the ASA analytical algorithm. Analytical derivatives of the atomic surfaces areas have been implemented. The atomic parameters have been obtained by a linear regression fit of the calculated and experimental free energies of solution in water for a set of 35 molecules, leading to a standard deviation of 0.75 kcal/mol. Effects of nonelectrostatic terms on solute geometries, association energies, and activation barriers are illustrated. © 1996 by John Wiley & Sons, Inc.     

The investigation on the relative stability of different clusters of thiourea dioxide in water using gas phase quantum chemical calculations

The relative stability of different clusters of thiourea dioxide (TDO) in water is examined using gas phase quantum chemical calculations at the MP2 and B3LYP level with 6‐311++G(d,p) basis set. The possible equilibrium structures and other energetic and geometrical data of the thiourea dioxide clusters, TDO‐(H₂O)_n (n is the number of water molecules), are obtained. The calculation results show that a strong interaction exists between thiourea dioxide and water molecules, as indicated by the binding energies of the TDO clusters progressively increased by adding water molecules. PCM model is used to investigate solvent effect of TDO. We obtained a negative hydration energy of ?20.6 kcal mol^?1 and free‐energy change of ?21.0 kcal mol^?1 in hydration process. On the basis of increasing binding energies with adding water molecules and a negative hydration energy by PCM calculation, we conclude thiourea dioxide can dissolve in water molecules. Furthermore, the increases of the C? S bond distance by the addition of water molecules show that the strength of the C? S bonds is attenuated. We find that when the number of water molecules was up to 5, the C? S bonds of the clusters, TDO‐(H₂O)₅ and TDO‐(H₂O)₆ were ruptured. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Evaluation of the 1-octanol/water partition coefficient of nucleoside analogs via free energy estimated in quantum chemical calculations

The partition coefficients (logP) of nucleoside analogs determined by the difference in the free energies of hydration and solvation in water-saturated octanol using the thermodynamic integration method are reported. The logP values calculated in this approach are closer to the experimental values compared to other ab initio methods. Solvation free energy in water and octanol, free energy of cavity formation in water and Henry’s constants, and some other parameters are estimated at the density functional theory (DFT) and Hartree-Fock level with 6–31G*, 6–31G, and 6–31+G basis sets. Surface area, mass, refractivity, volume, polarizability, and dipole moment are calculated for some drugs with HF and DFT methods. The results show that log P decreases with the decrease in polarizability and the increase in dipole moment.     

Quantum-chemical study of the molecular forms of vanillin

Equilibrium structural parameters, dipole moments, relative energies, vibration spectra, and ionization potentials of various conformers of benzoid (B), p- and o-quinoid (Q), and deprotonated anionic forms of vanillin in a vacuum and in a polarizable medium with dielectric properties of water were calculated by the quantum-chemical method (U)B3LYP/cc-pVTZ. The intramolecular hydrogen bond in free vanillin benzoid molecule B1 is characterized by the energy of 4 kcal mol⁻¹ and by the barrier of cleavage 8 kcal mol⁻¹ that correspond to the conformational transition B1 → B2. The energy of the optimal p-quinoid conformer (Q1) is 22 kcal mol⁻¹ larger than the energy of B1 in a vacuum and 15 kcal mol⁻¹ in water. The p-quinoid forms greatly exceeds the benzoid forms by dipole moment, electrophilicity, and nucleophilicity. The structure calculated for the anion is intermediate between benzoid and p-quinoid forms, and the calculated energy of polarization stabilization in water is 57 kcal mol⁻¹.     

An application of the reaction field theory to hydrated metal cations in the framework of the MNDO,AM1, and PM3 methods

A reaction field theory, combined with the MNDO, AM1, and PM3 molecular orbital methods, was applied to hydration phenomena of metal cationic species. The first hydration shell was treated explicitly by using a supermolecular model, [M(H₂O)_n]^m+, and its surrounding medium was described with a continuum dielectric. Hydration free energies were evaluated as a sum of the contributions from the electrostatic interaction with the bulk medium, the hydrated cluster formation, the cavity formation, and the vaporization of water molecules forming the cluster. As a whole, calculated hydration energies were in good agreement with the corresponding experimental data over various kinds of metal cationic species. © 1995 by John Wiley & Sons, Inc.     

Quantum wave packet calculation of reaction probabilities,cross sections,and rate constants for the C(1D) + HD reaction

The time‐dependent real wave packet method has been used to study the C(¹D) + HD reaction. The state‐to‐state and state‐to‐all reactive scattering probabilities for a broad range of energies are calculated at zero total angular momentum. The probabilities for J > 0 are estimated from accurately computed J = 0 probabilities by using the J‐shifting approximation. The integral cross sections for a large energy range, and thermal rate constants are calculated. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Theoretical Studies of Inorganic Compounds. 19 1) Quantum Chemical Investigations of the Phosphane Complexes X3B‐PY3 and X3Al‐PY3 (X = H,F, Cl; Y = F,Cl, Me,CN)

We report about quantum chemical ab initio calculations at the MP2/6‐311+G(2d)//MP2/6‐31G(d) level and DFT calculations at BP86/TZP of the geometries and bond dissociation energies of the borane‐phosphane complexes X₃B‐PY₃ and the alane‐phosphane complexes X₃Al‐PY₃ (X = H, F, Cl; Y = F, Cl, Me, CN). The nature of the B‐P and Al‐P bonds is analyzed with a bond energy partitioning method. The calculated bond dissociation energies D_e of the borane adducts X₃B‐PY₃ show for the phosphane ligands the trend PMe₃ > PCl₃ ∼ PF₃ > P(CN)₃. A similar trend PMe₃ > PCl₃ > PF₃ > P(CN)₃ is predicted for the alane complexes X₃Al‐PY₃. The order of the Lewis acid strength of the boranes depends on the phosphane Lewis base. The boranes show with PMe₃ and PCl₃ the trend BH₃ > BCl₃ > BF₃ but with PF₃ and P(CN)₃ the order is BH₃ > BF₃ > BCl₃. The bond energies of the alane complexes show always the trend AlCl₃ ≥ AlF₃ > AlH₃. The bonding analysis shows that it is generally not possible to correlate the trend of the bond energies with one single factor which determines the bond strength. The preparation energy which is necessary to deform the Lewis acid and Lewis base from the equilibrium form to the geometry in the complex may have a strong influence on the bond energies. The intrinsic interaction energies may have a different order than the bond dissociation energies. The trend of the interaction energies are sometimes determined by a single factor (Pauli repulsion, electrostatic attraction or covalent bonding) but sometimes all components are important. The higher Lewis acid strength of BCl₃ compared with BF₃ in strongly bonded complexes is not caused by the deformation energy of the fragments but it is rather caused by the intrinsic interaction energy. P(CN)₃ is a weaker Lewis base than PF₃, PCl₃ and PMe₃ mainly because of its weaker electrostatic attraction. The complex H₃B‐P(CN)₃ is predicted to have a bond dissociation energy D_o = 14.8 kcal/mol which should be sufficient to synthesize the compound as the first adduct with the ligand P(CN)₃. The calculated bond energies at the BP86 level are in most cases very similar to the MP2 results. In a few cases significantly different absolute values have been found which are caused by the method and not by the quality of the basis set.     

Gas‐phase and liquid‐state properties of esters,nitriles, and nitro compounds with the OPLS‐AA force field

Nonbonded and torsional parameters for carboxylate esters, nitriles, and nitro compounds have been developed for the OPLS‐AA force field. In addition, torsional parameters for alkanes have been updated. These parameters were fit to reproduce ab initio gas‐phase structures and conformational energetics, experimental condensed‐phase structural and thermodynamic properties, and experimental free energies of hydration. The computed densities, heats of vaporization, and heat capacities for fifteen liquids are in excellent agreement with experimental values. The new parameters permit accurate molecular modeling of compounds containing a wider variety of functional groups, which are common in organic molecules and drugs. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1340–1352, 2001     

Investigation of the use of density functionals in second‐ and third‐row transition metal dimer calculations

We explore the use of density functionals in calculating the equilibrium distances, dissociation energies, and harmonic vibrational frequencies of the homonuclear diatomics of the second‐row transition metals, platinum, and gold. The outermost s–d interconfigurational energies (ICEs) and the outermost s and d ionization potentials (IPs) were also calculated for the second‐ and third‐row transition metal atoms. Compared with the first‐row transition metal dimer calculations (J Chem Phys 2000, 112, 545–553), the binding energies calculated using the combination of the Becke 1988 exchange and the one‐parameter progressive correlation (BOP) functional and Becke's three‐parameter hybrid (B3LYP) functional are in better agreement with the experiment. However, the pure BOP functional still gives the deep and narrow dissociation potential wells for the electron configurations containing high‐angular‐momentum open‐shell orbitals. Analysis of the s–d ICEs and the s and d IPs suggests that the overestimation may be due to the insufficient long‐range interaction between the outermost s and d orbitals in the exchange functional. The hybrid B3LYP functional seems to partly solve this problem for many systems by the incorporation of the Hartree–Fock exchange integral. However, this still leads to an erroneous energy gap between the configurations of fairly different spin multiplicity, probably because of the unbalance of exchange and correlation contributions. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1995–2009, 2001