A new method for deriving atomic charges and dipoles for n, -alkanes: investigation of transferability and geometry dependence

Atomic charge models are known to be unsatisfactory for representing the ab initio electrostatic potential (ESP) of n, -alkanes. A new method for deriving atomic charges and dipoles is proposed and applied to n-alkanes ranging from C₄ to C₁₀. Electrostatic parameters found by this method reproduce accurately the ab initio ESP. The issues of transferability and conformational dependence are also addressed by introducing charges and dipoles taken from a truncated distributed multipole analysis, in the same spirit as the restrained electrostatic potential method. A transferable model is proposed for larger alkanes (>C₁₀). We also estimate the error made when using a set of Boltzmann-weighted electrostatic parameters for all conformers. The reduced number of electrostatic sites considered in our model makes it suitable for computer simulation of liquid n-alkanes.     

Experimental and theoretical study of the mechanism for the kinetic of elimination of methyl carbazate in the gas phase

The elimination kinetic of methyl carbazate in the gas phase was determined in a static system over the temperature range of 340–390 °C and pressure range of 47–118 Torr. The reaction is homogeneous, unimolecular, and obeys a first order rate law. The decomposition products are methyl amine, nitrous acid, and CO gas. The variation of the rate coefficients with temperatures is given by the Arrhenius expression: log k₁ (s^?1) = (11.56 ± 0.34) ? (180.7 ± 4.1) kJ mol^?1(2.303 RT)^?1. The estimated kinetics and thermodynamics parameters are in good agreement to the experimental values using B3LYP/6‐31G (d,p), and MP2/6‐31G (d,p) levels of theory. These calculations imply a molecular mechanism involving a concerted non‐synchronous quasi three‐membered ring cyclic transition state to give an unstable intermediate, 1,2‐oxaziridin‐3‐one. Bond order analysis and natural charges implies that polarization of O (alkyl)? C (alkyl) bond of the ester is rate determining in this reaction. Copyright © 2008 John Wiley & Sons, Ltd.     

Raman excitation of (+)‐(R)‐methyloxirane and its origin of optical activity via bond polarizabilities

We studied the bond polarizabilities of chiral (+)‐(R)‐methyloxirane from its Raman intensities. The bond polarizabilities provide much information concerning the electronic structure of its nonresonant Raman‐excited virtual state. At the initial moment of Raman excitation by the 514.5 nm laser, the tendency of the excited charges (mapped out by the bond polarizabilities) is to spread to the methine bond near the stereogenic center and its triangular oxirane skeleton. Thereby, the coupling of the electric dipole induced by the excited charges in the methine bond and the magnetic moment vibrationally induced by the electric current in the triangular oxirane skeleton as the molecule vibrates is shown to be the key factor leading to its significant Raman chirality. When the final stage of Raman relaxation is approached, the relative magnitudes of the bond polarizabilities are congruent to the bond electronic densities of the ground state, which are otherwise by the theoretical quantities via the quantum chemical calculation. During Raman relaxation, we found that the polarizabilities of the peripheral C H bonds relax faster than the rest, as indicated by their relaxation characteristic times. Copyright © 2010 John Wiley & Sons, Ltd.     

Aluminum alkoxide complexes from highly substituted 3,4‐epoxy alcohols: An NMR and computational study

The organoaluminum mediated epoxide ring opening of epoxy alcohols is a key step in the oxirane‐based approach for polypropionate synthesis. However, this reaction has shown unanticipated regioselectivities when applied to 2‐methyl‐3,4‐epoxy alcohols. In order to gain mechanistic insight into the factors controlling the epoxide ring opening process, diastereomeric 2‐methyl‐3,4‐epoxy alcohols were reacted with triethylaluminum in order to identify the aluminum complexes formed by these systems. Different epoxide–aluminum complexes were calculated using ab initio HF/[13s7p/11s5p] and B3PW91/6‐31G** gauge‐including atomic orbital (GIAO) methods and compared to the experimental NMR data. The calculated and experimental data correlates with the aluminum dimer complex (TIPSOCH₂CH(OAlEt₂)CH(CH)₃CHCH(O)(AlEt₃))₂ (VIII) for the systems favoring the nucleophilic attack at the external C4 epoxide carbon, while an unusual trialuminum species TIPSOCH₂CH(OAlEt₂)CH(CH)₃CHCH(O)(AlEt₃)₂ (X) is consistent with the systems favoring the internal C3 attack. The ²⁷Al NMR data established the tetracoordinated nature of the aluminum metal in all alkoxy aluminum intermediates, while the ¹³C NMR data provided insight into the aluminum‐oxygen coordination. The formation of the complexes was dictated by the stereochemical disposition of the substituents. These complexes are different from the generally accepted bidentate intermediates proposed for 2,3‐epoxy alcohols and simpler 3,4‐epoxy alcohols. Copyright © 2012 John Wiley & Sons, Ltd.     

Substituent effect investigation of 3‐(2, 4‐dichlorophenyl)‐1‐(4′‐X‐phenyl)‐2‐propen‐1‐one. Part 1. Correlation analysis of 13C NMR chemical shifts

A series of substituted chlorinated chalcones namely, 3‐(2,4‐dichlorophenyl)‐1‐(4′‐X‐phenyl)‐2‐propen‐1‐one, have been synthesized, X being H, NH₂, OMe, Me, F, Cl, CO₂Et, CN, and NO₂. Dual substituent parameter (DSP) models of ¹³C NMR chemical shift (CS) have revealed that π‐polarization concept could be utilized to explain the reverse field effect at CO, the enhanced substituent field effect at CO, C‐2, and C‐5, and the decreased sensitivity of substituent field effect at C‐6. Chlorine atoms dipole direction at the benzylidene ring either enhances or reduces substituent effect depending on how they couple with the substituent dipole at the probe site. The correlation of ¹³C NMR CS of C‐2, C‐5, and C‐6 with σ and σ indicates that chlorine atoms in the benzylidine ring deplete the ring from charges. Both MSP of Hammett and DSP of Taft ¹³C NMR CS models give similar trends of substituent effects at C‐2, C‐5, and C‐6. However, the former fail to give a significant correlation for CO and C‐6 ¹³C NMR CS. MSP of σ_q and DSP of Taft and Reynolds models significantly correlated ¹³C NMR CS of C_β. MSP of σ_q fails to correlate C‐1′ ¹³C NMR CS. Investigation of ¹³C NMR CS of non‐chlorinated chalcones series: 3‐phenyl‐1‐(4′‐X‐phenyl)‐2‐propen‐1‐one has revealed similar trends of substituent effects as in the chlorinated chalcones series for C‐1′, CO, C_α, and C_β. In contrast, the substituent effect of the non‐chlorinated chalcone series at C‐2, C‐5, and C‐6 did not correlate with any substituent constant. Copyright © 2010 John Wiley & Sons, Ltd.     

Unusual conformational preferences of 1,3‐dimethyl‐3‐isopropoxy‐3‐silapiperidine

The conformational analysis of the first representative of the Si‐alkoxy substituted six‐membered Si,N‐heterocycles, 1,3‐dimethyl‐3‐isopropoxy‐3‐silapiperidine, was performed by low‐temperature ¹H and ¹³C NMR spectroscopy and DFT theoretical calculations. In contrast to the expectations from the conformational energies of methyl and alkoxy substituents, the Me_axi‐PrO_eq conformer was found to predominate in the conformational equilibrium in the ratio Me_axi‐PrO_eq : Me_eqi‐PrO_ax of ca. 2 : 1 as from the ¹H and ¹³C NMR study. The thermodynamic parameters obtained by the complete line shape analysis showed that the main contribution to the barrier to ring inversion originates from the entropy term of the free energy of activation. Copyright © 2012 John Wiley & Sons, Ltd.     

Spectroscopic (FT‐IR,FT‐Raman,UV, 1H,and 13C NMR) and theoretical studies of m‐anisic acid and lithium,sodium, potassium,rubidium, and caesium m‐anisates

In order to understand the nature of the interactions of biologically important ligands, it is necessary to carry out the physico‐chemical studies of these compounds with their biological targets (e.g., receptors in the cell or important cell components). Results of this study make it possible to predict some properties of a molecule, such as its reactivity, durability of complex compounds, and kinship to enzymes. In this paper the effect of alkali metal cations (Li, Na, K, Rb, and Cs) on the electronic structure of m‐methoxybenzoic acid (m‐anisic acid) was studied. The experimental IR (in solid state and solution), Raman, UV (in solid state and solution), ¹H, and ¹³C NMR spectra of m‐methoxybenzoic acid, and its salts were registered, assigned, and analyzed. Some of the obtained results were compared with published data for o‐anisic acid and o‐anisates. The structures of anisic acid and Li, Na, and K m‐anisates were optimized at the B3LYP/6‐311++G** level. The IR, ¹H, and ¹³C NMR spectra and NPA, ChelpG, and MK atomic charges were calculated. The change of metal along with the series: Li → Na → K → Rb → Cs caused: (1) the change in the electronic charge distribution in anisate anion that is seen via the occurrence of the systematic shifts of several bands in the experimental and theoretical IR and Raman spectra of anisates; (2) systematic ¹H and ¹³C NMR chemical shifts; (3) hypsochromic shifts in UV spectra of salts as compared to ligands. Copyright © 2009 John Wiley & Sons, Ltd.     

Nucleophile‐electrophile interactions in the reaction of oxiranes with carboxylic acids in the presence of tertiary amines

The nucleophile‐electrophile interactions in the reaction system “N,N‐dimethylaniline – acetic acid – epichlorohydrin” have been investigated using kinetic methods and computer modeling. The observed orders of reactions have been determined for the overall reaction as well as for individual stages. The kinetic equations have been proposed; the activation parameters of the reactions have been evaluated. The behavior of the initial (amine) and intermediate (carboxylate) nucleophiles has been detailed in reaction pathway. Amine reacts with oxirane activated by acidic reagent while carboxylate‐anion—with both activated and nonactivated epichlorohydrin. The mechanism of oxirane ring opening by acid reagent in the presence of tertiary amine has been proposed, which comprise parallel‐consecutive compensation stages of reaction product formation. It has been demonstrated that the observed reaction order with respect to acid reactant depends upon the nature of electrophilic reagent (activated/nonactivated oxirane) and the ratio of the rates of compensation stages.     

Possibility of non‐adiabatic level crossing by DFT study of tautomerism and potential energy surfaces in of 3‐hydroxy‐5‐(pyrimidin‐2‐yl)‐2H‐pyrrol‐2‐one and its tautomer

3‐Hydroxy‐5‐(pyrimidin‐2‐yl)‐2H‐pyrrol‐2‐one (HYPO, T1) and 2‐hydroxy‐5‐(pyrimidine‐2‐yl)‐3H‐pyrrole‐3‐one (HYPO, T2) have designed in this research to study potential energy curves for their dynamic motions and possibility of crossing between levels. Study of tautomerism shows that T1 tautomer is more stable than T2 (about 5.83 kJ/mol). Dynamic study of possible motions show rate constants (highest possible) equal to 8.82 M/s for tautomerism, 1.70 × 10⁹ M/s for relative rotation of ring (rr) and 3.67 × 10⁶ M/s for rotation of OH bond (br). Moreover, variations of orbital populations, NBO charges, hybridations, and acceptor–donor interactions in IRC steps have been investigated to study the possibility of non‐adiabatic crossing between tautomerism and ring rotation potential energy curves. The data showed that in spite of the fact that these two potentials share three common points, these two potential curves cannot have non‐adiabatic crossing because of different symmetries and a large difference between their barrier energies. Copyright © 2010 John Wiley & Sons, Ltd.     

In situ Grafting onto Silica Surface with Epoxidized Natural Rubber via Solid State Method

An in situ solid state grafting reaction between epoxidized natural rubber (ENR) and silica was performed in a Haake internal mixer. Resulting ENR‐grafted silica was characterized by Fourier transform infrared spectroscopy (FTIR) and thermogravimetric analysis (TGA) measurements. Based on these results, it was concluded the silanol groups (Si‐OH) of silica caused the ring opening of ENR oxirane rings so that ENR was grafted onto the silica surface. Transmission electron microscopy (TEM) photographs showed ENR‐grafted silica had better dispersibility and smaller aggregates compared with the original silica. Dynamical mechanical analysis (DMA) of vulcanized rubber compounds contained ENR‐grafted silica showed the glass transition temperature (T _g) of grafted ENR molecules shifted to higher temperature, from ?3°C to 20°C, indicating the mobility of ENR was greatly restricted. As a result, the compounds containing ENR‐grafted silica have higher hysteresis, and can be applied in a much wider field, such as damping materials, tires of racing cars, and so on.     

Basicity of trifluoromethylsulfonylformamidines. DFT and FTIR study and NBO analysis

The effect of the electron–acceptor substituent CF₃SO₂ at the imine nitrogen atom on the basicity and the electron distribution in N,N‐alkylformamidines ( 1 , 2 , 3 , 4 , 5 ) was studied experimentally by the FTIR spectroscopy and theoretically at the DFT (B3LYP/6‐311+G(d,p)) level of theory, including the natural bond orbital (NBO) analysis. The calculated proton affinities of the imine nitrogen atom and the sulfonyl oxygen (PA_N′ and PA_O) depend on the atomic charges, the C?N′ and N′―S bond polarity and on the energy of interaction of the amine nitrogen and the oxygen lone pairs with antibonding π* and σ*‐orbitals. The basicity of the imine nitrogen atom is increased with the increase of the electron‐donating power of the substituent at the amine nitrogen atom due to stronger interaction n_N → π*_C?N′, but is decreased for the electron‐withdrawing groups MeSO₂ and CF₃SO₂ at the imine nitrogen atom in spite of the increase of this conjugation. Protonation of ( 1 , 2 , 3 , 4 , 5 ) in CH₂Cl₂ solution in the presence of CF₃SO₃H occurs at the imine nitrogen atom, while the formation of hydrogen bonds with 4‐fluorophenol takes place at the sulfonyl oxygen atom, whose basicity is lower than that of N,N′‐dimethylmethanesulfonamide but higher than of N,N′‐dimethyltrifluoromethanesulfonamide. Copyright © 2015 John Wiley & Sons, Ltd.     

Theoretical calculations for neighboring group participation in gas‐phase elimination kinetics of 2‐hydroxyphenethyl chloride and 2‐methoxyphenethyl chloride

Theoretical calculation of the kinetics and mechanisms of gas‐phase elimination of 2‐hydroxyphenethyl chloride and 2‐methoxyphenethyl chloride has been carried out at the MP2/6‐31G(d,p), B3LYP/6‐31G(d,p), B3LYP/6‐31 + G(d,p), B3PW91/6‐31G(d,p) and CCSD(T) levels of the theory. The two substrates undergo parallel elimination reactions. The first process of elimination appears to proceed through a three‐membered cyclic transition state by the anchimeric assistance of the aromatic ring to produce the corresponding styrene product and HCl. The second process of elimination occurs through a five‐membered cyclic transition state by participation of the oxygen of o‐OH or the o‐OCH₃ to yield in both cases benzohydrofuran. The B3PW91/6‐31G(d,p) method was found to be in good agreement with the experimental kinetic and thermodynamic parameters for both substrates in the two reaction channels. However, some differences in the performance of the different methods are observed. NBO analysis of the pyrolysis of both phenethyl chlorides implies a C? Cl bond polarization, in the sense of C^δ+…Cl^δ?, which is a rate‐determining step for both parallel reactions. Synchronicity parameters imply polar transition states of these elimination reactions. Copyright © 2009 John Wiley & Sons, Ltd.     

Inner‐shell chemical shift of DNA/RNA bases and inheritance from their parent purine and pyrimidine

Inner‐shell electronic structures, properties and ionization spectra of DNA/RNA bases are studied with respect to their parent pyrimidine and purine species. Density functional theory B3LYP/aug‐cc‐pVTZ has been employed to produce the geometries of the bases, whereas LB94/et‐pVQZ//B3LYP/aug‐cc‐pVTZ is used to calculate site‐related Hirshfeld charges and core (vertical) ionization energies, as well as inner‐shell spectra of C1s, N1s and O1s for DNA/RNA bases and their parent pyrimidine and purine species. The site‐dependent variations of properties indicate the changes and inheritance of chemical environment when pyrimidine and purine become substituted. In general, although the changes are site‐dependent, they are also ring‐dependent. Pyrimidine bases change less significantly with respect to their parent pyrimidine than the purine bases with respect to their parent purine. Pyrimidine bases such as uracil, thymine and cytosine inherit certain properties from their parent pyrimidine, such as the Hirshfeld charge distributions and the order of core ionization energy level etc. No particular sites in the pyrimidine derivatives are engaged with a dramatic chemical shift nor with energy crossings to other sites. For the core shell spectra, the purine bases inherit very little from their parent purine, and guanine exhibits the least similarities to the parent among all the DNA/RNA bases.     

Chemically triggered C–ON bond homolysis in alkoxyamines. Part 7. Remote polar effect

In a recent work (Org. Lett. 2012, 14, 358–361), we showed that the activation by benzylation of alkoxyamine 1 (diethyl (1‐(tert‐butyl(1‐(pyridin‐4‐yl)ethoxy)amino)‐2,2‐dimethylpropyl)phosphonate) afforded a surprisingly large C–ON bond homolysis rate constant k_d. Taking advantage of the easy preparation of para‐X‐benzyl‐activated alkoxyamines 2 and of the presence of a shielding methylene group between the two aromatic moieties, we investigated the long range (10 bonds between the X group and the C–ON bond) polar effect for X = H, F, OMe, CN, NO₂, NMe₂, ⁺NHMe₂,Br^?. It was observed that the effect was weak (4‐fold) and mainly due to the zwiterionic mesomeric forms generated by the presence of group X on the para position, i.e. k_d increased for CN and NO₂ and decreased for OMe, NMe₂ and ⁺NMe₂H,Br^?. DFT calculations at the B3LYP/6‐31G(d,p) level were performed to determine orbital interactions (natural bond orbital (NBO) analysis), Mulliken and NBO charges which support the reactivity described. Copyright © 2014 John Wiley & Sons, Ltd.     

An optimal point-charge model for molecular electrostatic potentials

Motivated by the Legendre expansion of the electrostatic potential (ESP), we propose a method for obtaining atomic point-charges for a molecule based on reproducing the low-order multipole moments of the system. The resulting multipole-derived charges (MDCs) are well defined, do not require sampling of the ESP at grid points around the molecule and provide excellent reproduction of the electrostatic potential. No constraints are placed on the magnitude of the atomic charges.     

Physical parameters for atomic inner‐shell photoionization processes and analytical applications: a status report

The atomic inner‐shell vacancy decay processes comprising of radiative and non‐radiative transitions are characterized by the physical parameters, namely, the photoionization cross‐sections; X‐ray, Auger and Coster–Kronig (CK) transition rates; fluorescence and CK yields; and the vacancy transfer probabilities. These parameters are required to calculate the K‐shell and L_i (i = 1–3)/M_i (i = 1–5) sub‐shell X‐ray production cross‐sections and relative intensities which, in turn, are needed for different analytical applications. This report intended to provide a detailed account of the currently available data sets of different physical parameters for use in various analytical applications. Copyright © 2011 John Wiley & Sons, Ltd.     

Pseudo‐homogeneous kinetic model for esterification of acetic acid with propanol isomers over dodecatungstophosphoric acid supported on montmorillonite k10

Esterification of acetic acid with propanol isomers such as isopropanol and n‐propanol was carried out over dodecatungstophosphoric acid (DTPA), dodecamolybdophosphoric acid ammonium salt hydrate, and sodium tungstate hydrated purified supported on montmorillonite K10, which were characterized by powder X‐ray diffraction, Brunauer–Emmett–Teller, and temperature programmed ammonia desorption. A pseudo‐homogeneous (P‐H) kinetic model was established for esterification of acetic acid with propanol isomers over DTPA supported on montmorillonite K10. Effects of various parameters such as reaction time, speed of agitation, particle size, temperature, percent catalyst loading, molar ratio and mixture of propanol isomer were investigated in detail. The 20% (w/w) DTPA/K10 was found to be an optimum solid catalyst with 82% n‐propanol and 53% isopropanol conversion with 100% selectivity toward propyl acetate. The 20% (w/w) DTPA/K10 catalyst was found to be reusable for three cycles. The reaction follows second‐order kinetics with activation energies of 25.53 kJ mol^?1 and 28.15 kJ mol^?1 for isopropanol and n‐propanol, respectively. Pseudo‐homogeneous kinetic model fitted with R² value of trend line 0.999. This implies that esterification reaction is kinetically controlled owing to high activation energy. Copyright © 2013 John Wiley & Sons, Ltd.