Polymerization of epoxidized triglycerides with fluorosulfonic acid

The use of triglycerides as agri-based renewable raw materials for the development of new products is highly desirable in view of uncertain future petroleum prices. A new method of polymerizing epoxidized soybean oil has been devised with the use of fluorosulfonic acid. Depending on the reaction conditions, one can get a viscous oil or a solid as the reaction product. The reaction mechanism has been studied through a systematic examination of the fluorosulfonic acid-initiated reaction of epoxidized methyl oleate and epoxidized methyl linoleate. Through ¹³C NMR spectroscopy, the major species derived from fluorosulfonic acid-initiated polymerization of epoxidized soybean oil have been determined. The effects of temperature, initiator dosage, and reaction time have been studied and shown to affect the nature of polymer products obtained and the distribution of different chemical species present.     

Protonation and deprotonation of hydroxamic acids. An MO ab initio study

Ab initio molecular orbital calculations with 4-31G//4-31G, 6-31G^*//4-31G and 6-31+G//4-31G basis sets have been used to examine the structure, relative energy, protonation and deprotonation of a series of seven hydroxamic acids in the gas phase. The results show that hydroxamic acids are predominantly in the E-TS form and that the most probable protonation site is the carbonyl oxygen atom, while deprotonation proceeds by loss of NH hydrogen.     

Density functional theory for efficient ab initio molecular dynamics simulations in solution

We present a density functional for first-principles molecular dynamics simulations that includes the electrostatic effects of a continuous dielectric medium. It allows for numerical simulations of molecules in solution in a model polar solvent. We propose a smooth dielectric model function to model solvation into water and demonstrate its good numerical properties for total energy calculations and constant energy molecular dynamics.     

An ab initio study of pyruvic acid

The geometries and vibrational frequencies of two conformers of pyruvic acid have been obtained at the ab initio second order Möller-Plesset level of theory using the 6-311++G^** basis set. While the calculated geometries have been compared to the experimental microwave data, the vibrational frequencies have been assigned, using the experimental gas phase IR spectra of 13 isotopes of pyruvic acid by a recently developed scaling procedure (IRPROG). An attempt has been made to explain the stability of the eclipsed conformation over the staggered conformation of pyruvic acid by taking account of the molecular orbitals.     

Variable ζ calculation for self-consistent screening parameters in ab initio molecular theory

A resolution of Roothaan's HF–SCF–LCAO–MO equations is proposed in which atomic orbital exponents (ζ) are made dependent on the molecular charge distribution and included in the self-consistent scheme. Screening parameters so obtained are self-consistent with the molecular orbital coefficients and compare closely to optimum orbital exponents found by other methods. The technique is applied to the ground, lowest positive, and lowest negative ion states of the hydride series LiH, BH, and HF. Calculated potential curves are used to define purely theoretical values for the vertical and adiabatic ionization energies and electron affinities. Predictions are compared to experimental values where available.     

Protonation acidity constants for benzotoluidides in sulfuric acid solutions

The protonation of o-, m-and p-benzotoluidide in sulfuric acid solutions is studied by UV spectroscopy in the 190–350 nm region. Principal component analysis is applied to estimate the contributions of the effect of protonation and the medium effect. For the substances studied in this work, the first principal component (PC) captures about 98 % of the variance and the second PC ∼100 % of the cumulative percentage variance in the 210–350 nm region. The same spectral region is used for calculation of the ionization ratio from the coefficients of the first PC and mole fractions of the base and its conjugate acid. Using these data and Hammett's equation (pK _{BH +} = H _X + log I), the pK _{BH +} values for the protonation reaction are obtained. The dissociation constants as well as the solvent parameters m* (∼0.43) and ϕ (∼0.60) are calculated using the Excess Acidity Method (-pK _{BH +}= 2.28–2.30) and the Bunnett-Olsen Method (-pK _{BH +}= 2.24–2.28). The probable sites of protonation are discussed.     

Insights into amine-based CO<Subscript>2</Subscript> capture: an ab initio self-consistent reaction field investigation

Ab initio many-body perturbation theory (MP2/6-311++G(,dp)), density functional theory (B3LYP/6-31++G(d,p)) and self-consistent reaction field (IEF-PCM UA HF/6-31G(d)) calculations have been used to study the CO₂ capture reagents NH₃, 2-hydroxyethylamine (MEA), diaminoethane (EN), 2-amino-1-propanol (2A1P), diethanolamine (DEA), N-methyl-2-hydroxyethylamine (N-methylMEA), 2-amino-2-methyl-1-propanol (AMP), trishydroxymethylaminomethane (tris), piperazine (PZ) and piperidine (PD). This study involved full conformational searches of the capture amines in their native and protonated forms, and their carbamic acid and carbamate derivatives. Using this data, we were able to compute Boltzmann-averaged thermodynamic values for the amines, carbamates and carbamic acid derivatives, as well as equilibrium constants for a series of ‘universal’ aqueous capture reactions. Important findings include (i) relative pK _a values for the carbamic acid derivatives are a useful measure of carbamate stability, due to a particular chemical resonance which is also manifest in short computed N–CO₂H bonds at both levels of theory, (ii) the computational results for sterically hindered amines such as AMP and tris are consistent with these species forming carbamates which readily hydrolyse and (iii) the amine-catalysed reaction between OH⁻ and CO₂ to generate bicarbonate correlates with amine pK _a. Thermodynamic data from the ab initio computations predicts that the heterocycles PD and PZ and the acyclic sorbent EN are good choices for a capture solvent. AMP and tris perform poorly in comparison.     

Protonation of 18-crown-6 by dichloropicric acid in anhydrous and in water saturated 1,2-dichloroethane. Homoconjugation of dichloropicric acid

From conductometric and UV-VIS spectrophotometric studies of the reaction between 18-crown-6 (L) and dichloropicric acid (HA) in dry and water saturated 1,2-dichloroethane, it has been concluded that formation of a 1:1 homoconjugate HA ₂ ^– accompanies the simple protonation of L, viz, L+HALH⁺A^– and L+2HALH⁺HA ₂ ^– . The electrolytes LH⁺A^– and LH⁺HA ₂ ^– are extensively, or practically completely dissociated in both solvents under the experimental conditions. The specie LH⁺A^– appears to be a contact ion pair in DCE. The stability constant of HA ₂ ^– in the dry solvent, 5.7×10³ mol^–1-cm³, is some 10^2.4 times that in propylene carbonate reflecting the difference in H-bond accepting capacity of the two solvents. Hydration of HA, A^– and HA ₂ ^– in wet dichloroethane is negligible or slight. As expected, LH⁺ is rather strongly hydrated, the ratio of the hydration constants of LH⁺ and L being about 1×10¹.     

The Diels-Alder reaction of ethene and 1,3-butadiene: an extended multireference ab initio investigation.

The concerted and stepwise mechanisms of the Diels-Alder reaction between 1,3-butadiene and ethene have been investigated using highly correlated multireference methods (MRAQCC) and extended basis sets. Full MRAQCC geometry optimizations have been performed in all cases. The best estimate for the energy barrier of the Diels-Alder reaction is 22 kcalmol(-1). Anti- and gauche-out minima for the biradical structures and corresponding fragmentation saddle points have been determined. The biradical anti fragmentation saddle point is located 6.5 kcalmol(-1) above the concerted saddle point. The gauche-in structure does not correspond to a local minimum, but leads on geometry optimization directly to cyclohexene.     

Protonation of thymine,cytosine, adenine,and guanine DNA nucleic acid bases: Theoretical investigation into the framework of density functional theory

Gradient-corrected density functional computations with triple-zeta-type basis sets were performed to determine the preferred protonation site and the absolute gas-phase proton affinities of the most stable tautomer of the DNA bases thymine (T), cytosine (C), adenine (A), and guanine (G). Charge distribution, bond orders, and molecular electrostatic potentials were considered to rationalize the obtained results. The vibrational frequencies and the contribution of the zero-point energies were also computed. Significant geometrical changes in bond lengths and angles near the protonation sites were found. At 298 K, proton affinities values of 208.8 (T), 229.1 (C), 225.8 (A), and 230.3 (G) kcal/mol were obtained in agreement with experimental results. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 989–1000, 1998     

Ionic dissociations of chlorosulfonic acid in microsolvated clusters: A density functional theory and ab initio MO study

Ionic dissociation of chlorosulfonic acid (HSO3Cl) in the molecular clusters HSO3Cl-(H2O)n (n = 1-4) and HSO3Cl-NH3-(H2O)n (n = 0-3) was investigated by density functional theory and ab initio molecular orbital theory. The equilibrium structures, binding energies, and thermodynamic properties, such as relative enthalpy and relative Gibbs free energy, and were calculated using the hybrid density func- tional (B3LYP) method and the second order M?ller-Plesset approximation (MP2) method with the 6-311 G** basis set. Chlorosulfonic acid was found to require a minimum of three water molecules for ionization to occur and at least one water molecule to protonate ammonia. The corresponding clusters with fewer water molecules were found to be strongly hydrogen-bonded. The related properties and acid strength of chlorosulfonic acid were discussed and compared to the acid strengths of perchloric acid and sulfuric acid in the context of clusters with ammonia and water. The relative stabilities of these clusters were also investigated.     

DFT and ab initio quantum chemical studies on p-cyanobenzoic acid

The Fourier transform infrared (FTIR) and FT-Raman spectra of p-cyanobenzoic acid (CBA) have been recorded in the range 4000-400 and 4000-100 cm(-1), respectively. The complete vibrational assignment and analysis of the fundamental modes of the compound were carried out using the observed FTIR and FT-Raman data. The vibrational frequencies determined experimentally were compared with theoretical wavenumbers obtained from ab initio HF and DFT-B3LYP gradient calculations employing 6-31G**, 6-311++G** and cc-pVTZ basis sets for the optimised geometry of the compound. The geometry and normal modes of vibration obtained from the HF and DFT methods are in good agreement with the experimental data. The normal coordinate analysis was also carried out with ab initio force fields utilising Wilson's FG matrix method. The interactions of cyano and carboxylic acid groups with the skeletal vibrational modes were investigated.     

Development of Abraham model expressions for predicting the enthalpies of solvation of solutes dissolved in acetic acid

Experimental data have been compiled from the published chemical and engineering literature on the enthalpies of solvation for 92 different inorganic gases and organic vapours in acetic acid. The compiled data are used to derive mathematical correlations based on the Abraham solvation parameter model. The derived expressions describe the experimental solvation enthalpies in acetic acid to within a standard deviation (SD) of 2.2 kJ mol^?1. Principal Component Analysis (PCA) on the five equation coefficients from a derived Abraham model correlation shows that acetic acid does not resemble hydroxylic solvents in terms of enthalpic interactions, but is more akin to moderately polar solvents such as ethyl acetate or acetone.     

Applications of self-consistent field theory in polymer systems

1Introduction Owing to the specificity of the long chain,polymers present complexity and versatility.These molecules in the system can be various in their topological struc-tures,such as linear,star,comb or circle structures;meanwhile they can be polymeri…     

Aggregate morphologies of amphiphilic ABC triblock copolymer in dilute solution using self-consistent field theory

The complex microstructures of amphiphilic ABC linear triblock copolymers in which one of the end blocks is relatively short and hydrophilic, and the other two blocks B and C are hydrophobic in a dilute solution, have been investigated by the real-space implementation of self-consistent field theory (SCFT) in two dimensions (2D). In contrast to diblock copolymers in solution, the aggregation of triblock copolymers are more complicated due to the presence of the second hydrophobic blocks and, hence, big ranges of parameter space controlling the morphology. By tailoring the hydrophobic degree and its difference between the blocks B and C, the various shapes of vesicles, circlelike and linelike micelles possibly corresponding to spherelike, and rodlike micelles in 3D, and especially, peanutlike micelles not found in diblock copolymers are observed. The transition from vesicles to circlelike micelles occurs with increasing the hydrophobicity of the blocks B and C, while the transition from circlelike micelles to linelike micelles or from the mixture of micelles and vesicles to the long linelike micelles takes place when the repulsive interaction of the end hydrophobic block C is stronger than that of the middle hydrophobic block B. Furthermore, it is favorable for dispersion of the block copolymer in the solvent into aggregates when the repulsion of the solvent to the end hydrophobic block is larger than that of the solvent to the middle hydrophobic block. Especially when the bulk block copolymers are in a weak segregation regime, the competition between the microphase separation and macrophase separation exists and the large compound micelle-like aggregates are found due to the macrophase separation with increasing the hydrophobic degree of blocks B and C, which is absent in diblock copolymer solution. The simulation results successfully reproduce the existing experimental ones.     

Update on phosphate and charged post‐translationally modified amino acid parameters in the GROMOS force field

In this study, we propose newly derived parameters for phosphate ions in the context of the GROMOS force field parameter sets. The non‐bonded parameters used up to now lead to a hydration free energy, which renders the dihydrogen phosphate ion too hydrophobic when compared to experimentally derived values, making a reparametrization of the phosphate moiety necessary. Phosphate species are of great importance in biomolecular simulations not only because of their crucial role in the backbone of nucleic acids but also as they represent one of the most important types of post‐translational modifications to protein side‐chains and are an integral part in many lipids. Our re‐parametrization of the free dihydrogen phosphate (H PO ) and three derivatives (methyl phosphate, dimethyl phosphate, and phenyl phosphate) leads, in conjunction with the previously updated charged side‐chains in the GROMOS parameter set 54A8, to new nucleic acid backbone parameters and a 54A8 version of the widely used GROMOS protein post‐translational modification parameter set. © 2017 Wiley Periodicals, Inc.     

Study of singlet-triplet transitions in the ozone molecule using the multiconfigurational self-consistent field theory

The intensities and dipole moments of the lower singlet-triplet transitions ³ A ₂ ← X ¹ A ₁ and ³ B ₁ ← X ¹ A ₁ in the ozone molecule were calculated by the multiconfiguration self-consistent field theory with the quadratic response function. The results of calculations of the intensities of singlet-triplet transitions using different basis sets and complete active spaces were compared. The assignment of the ³ A ₂ ← X ¹ A ₁ transition in the ozone spectrum to the Wulf band is discussed.     

Implicit solvation in the self-consistent mean field theory method: sidechain modelling and prediction of folding free energies of protein mutants

The Atomic Solvation Parameter (ASP) model is one of the simplest models of solvation, in which the solvation free energy of a molecule is proportional to the solvent accessible surface area (SAS) of its atoms. However, until now this model had not been incorporated into the Self-Consistent Mean Field Theory (SCMFT) method for modelling sidechain conformations in proteins. The reason for this is that SAS is a many-body quantity and, thus, it is not obvious how to define it within the Mean Field (MF) framework, where multiple copies of each sidechain exist simultaneously. Here, we present a method for incorporating an SAS-based potential, such as the ASP model, into SCMFT. The theory on which the method is based is exact within the MF framework, that is, it does not depend on a pairwise or any other approximation of SAS. Therefore, SAS can be calculated to arbitrary accuracy. The method is computationally very efficient: only 7.6% slower on average than the method without solvation. We applied the method to the prediction of sidechain conformation, using as a test set high-quality solution structures of 11 proteins. Solvation was found to substantially improve the prediction accuracy of well-defined surface sidechains. We also investigated whether the methodology can be applied to prediction of folding free energies of protein mutants, using a set of barnase mutants. For apolar mutants, the modest correlation observed between calculated and observed folding free energies without solvation improved substantially when solvation was included, allowing the prediction of trends in the folding free energies of this type of mutants. For polar mutants, correlation was not significant even with solvation. Several other factors also responsible for the correlation were identified and analysed. From this analysis, future directions for applying and improving the present methodology are discussed.     

One- and three-parameter density functional theory and high level ab initio theory study of the CH+CH disproportionation reaction

Complete basis set (CBS) ab initio computational studies were performed with the target being to explore the CH+CH potential energy surface. Several closed and open shell intermediates were located on the potential energy surface. Computed enthalpies for the branching reactions, as well as heats of formation are in excellent agreement. Although CBS computed energies are of high quality, this computational study is not capable of predicting the branching product ratio due to fact that neither the MP2 nor the 6-311G(2d,2p) basis set are sufficient to locate the reactant complexes and the transition state structures for the hydrogen and carbon transfer reactions in the reaction complexes. To properly explore the CH+CH potential energy surface a much higher ab initio theory level is required.