Thermal reactivities of polynuclear aromatic hydrocarbons and alkyl derivatives

The HMO free valence, the PMO Dewar reactivity number, and the structure count ratio are correlated with thermal reactivities of polynuclear aromatic hydrocarbons and derivatives. The compounds are classified as thermally reactive if carbonaceous residues are observed upon heating to 750° (Experimental, Ref. 13). Radical intermediates are assumed to mediate the pyrolytic reactivities. The structure count ratio which is the ratio SC (radical intermediate) ÷ SC (reactant) is recommended as a simple working method for correlating and predicting thermal reactivity.     

How to compute isomerization energies of organic molecules with quantum chemical methods

The reaction energies for 34 typical organic isomerizations including oxygen and nitrogen heteroatoms are investigated with modern quantum chemical methods that have the perspective of also being applicable to large systems. The experimental reaction enthalpies are corrected for vibrational and thermal effects, and the thus derived "experimental" reaction energies are compared to corresponding theoretical data. A series of standard AO basis sets in combination with second-order perturbation theory (MP2, SCS-MP2), conventional density functionals (e.g., PBE, TPSS, B3-LYP, MPW1K, BMK), and new perturbative functionals (B2-PLYP, mPW2-PLYP) are tested. In three cases, obvious errors of the experimental values could be detected, and accurate coupled-cluster [CCSD(T)] reference values have been used instead. It is found that only triple-zeta quality AO basis sets provide results close enough to the basis set limit and that sets like the popular 6-31G(d) should be avoided in accurate work. Augmentation of small basis sets with diffuse functions has a notable effect in B3-LYP calculations that is attributed to intramolecular basis set superposition error and covers basic deficiencies of the functional. The new methods based on perturbation theory (SCS-MP2, X2-PLYP) are found to be clearly superior to many other approaches; that is, they provide mean absolute deviations of less than 1.2 kcal mol-1 and only a few (<10%) outliers. The best performance in the group of conventional functionals is found for the highly parametrized BMK hybrid meta-GGA. Contrary to accepted opinion, hybrid density functionals offer no real advantage over simple GGAs. For reasonably large AO basis sets, results of poor quality are obtained with the popular B3-LYP functional that cannot be recommended for thermochemical applications in organic chemistry. The results of this study are complementary to often used benchmarks based on atomization energies and should guide chemists in their search for accurate and efficient computational thermochemistry methods.     

Comparison of different extraction methods for the determination of polycyclic aromatic hydrocarbons in soil

Two different extraction methods for the determination of polycyclic aromatic hydrocarbons in soil are compared: the extraction in combination with ultrasonic treatment, and the Soxhlet extraction method according to DIN-draft 38414 Part 21. Different types of real soil were extracted and analysed by HPLC with diode-array detector and fluorescence detection. The results show that the efficiency of the ultrasonic method is comparable to the Soxhlet method.     

Comparison of static first hyperpolarizabilities calculated with various quantum mechanical methods

The prediction of nonlinear electro-optic (EO) behavior of molecules with quantum methods is the first step in the development of organic-based electro-optic devices. Typical EO molecules may require calculations with several hundred electrons, which prevents all but the fastest methods (semiempirical and density functional theory (DFT)) from being used for EO estimation. To test the reliability of these methods, we compare dipole moments, polarizabilities, and first-order hyperpolarizabilities for a wide range of structures of experimental interest with Hartree-Fock (HF), intermediate neglect of differential overlap (INDO), and DFT methods. The relative merits of molecules are consistently predictable with every method.     

Polythiol binding to biologically relevant metal ions

The coordination properties of three peptides with CXXC motif: Ac-GCASCDNCRACKK-NH(2), Ac-GCASCDNCRAAKK-NH(2) and Ac-GCASCDNARAAKK-NH(2) as donors of four, three and two thiol ligands for Ni(2+),Cd(2+), Zn(2+) and Bi(3+) were studied by potentiometric titrations, UV-Vis and CD spectra measurements. Since the stability of the complexes is closely connected with the amount of the metal-bound cysteine sulfurs, competition plots of the complexes of peptides with 2, 3 and 4 cysteines further prove the involvement of all thiols in the metal ion binding. Furthermore, the sulfur-bound zinc complexes appear to be much more stable than the sulfur-bound nickel ones. The stabilities of the studied complexes decreases in the series Bi(3+) ? Cd(2+) > Zn(2+) > Ni(2+).     

Comparison of two simulation methods to compute solvation free energies and partition coefficients

The thermodynamic integration (TI) and expanded ensemble (EE) methods are used here to calculate the hydration free energy in water, the solvation free energy in 1‐octanol, and the octanol‐water partition coefficient for a six compounds of varying functionality using the optimized potentials for liquid simulations (OPLS) all‐atom (AA) force field parameters and atomic charges. Both methods use the molecular dynamics algorithm as a primary component of the simulation protocol, and both have found wide applications in fields such as the calculation of activity coefficients, phase behavior, and partition coefficients. Both methods result in solvation free energies and 1‐octanol/water partition coefficients with average absolute deviations (AAD) from experimental data to within 4 kJ/mol and 0.5 log units, respectively. Here, we find that in simulations the OPLS‐AA force field parameters (with fixed charges) can reproduce solvation free energies of solutes in 1‐octanol with AAD of about half that for the solute hydration free energies using a extended simple point charge (SPC/E) model of water. The computational efficiency of the two simulation methods are compared based on the time (in nanoseconds) required to obtain similar standard deviations in the solvation free energies and 1‐octanol/water partition coefficients. By this analysis, the EE method is found to be a factor of nine more efficient than the TI algorithm. For both methods, solvation free energy calculations in 1‐octanol consume roughly an order of magnitude more CPU hours than the hydration free energy calculations. © 2012 Wiley Periodicals, Inc.     

Chromatographic methods for carcinogenic/mutagenic nitropolycyclic aromatic hydrocarbons

Nitropolycyclic aromatic hydrocarbons (NPAHs), which are known to be carcinogenic and/or mutagenic, are considered to be one of the air pollutants that cause lung cancer. In the last two decades, a number of sensitive and selective methods have been developed for the determination of NPAHs and related compounds in environmental and biological samples. This paper describes the state of the art of the methods and applications.     

A quantum mechanical test of diophantine methods

A blend of Haselgrove's method and the biased selection method for evaluating multidimensional integrals was tested. The results were mixed. The error estimate varied from being proportional to 1/N when N was less than ca. 60,000 to being proportional to 1/√N when N was greater than 60,000. Also, for N greater than 60,000, the error estimate was one-half the error estimate given by biased selection alone. These numbers should be compared with the 10,000 points used to find an optimum set of Haselgrove's parameters. It is reasonable to expect that if 100,000 points were used in the optimization of Haselgrove's parameters that the above results would be found with 60,000 replaced by 600,000.     

A benchmark quantum chemical study of the stacking interaction between larger polycondensed aromatic hydrocarbons

Large-scale electronic structure calculations were performed for the interaction energy between coronene, C₂₄H₁₂ with circumcoronene, C₅₄H₁₈, and between two circumcoronene molecules, in order to get a picture of the interaction between larger graphene sheets. Most calculations were performed at the SCS-MP2 level but we have corrected them for higher-order correlation effects using a calculation on the coronene-circumcoronene system at the quadratic CI, QCISD(T) level. Our best estimate for the interaction energy between coronene and circumcoronene is 32.1?kcal/mol. We estimate the binding of coronene on a graphite surface to be 37.4 or 1.56?kcal/mol per carbon atom (67.5?meV/C atom). This is also our estimate for the exfoliation energy of graphite. It is higher than most previous theoretical estimates. The SCS-MP2 method which reproduces the CCSD(T) and QCISD(T) values very well for smaller aromatic hydrocarbons, e.g., for the benzene dimer, increasingly overestimates dispersion as the bandgap (the HOMO-LUMO separation) decreases. The barrier to the sliding motion of coronene on circumcoronene is 0.45?kcal/mol, and for two circumcoronene molecules 1.85?kcal/mol (0.018 and 0.034?kcal/mol per C atom, respectively). This means that larger graphenes cannot easily glide over each other.     

Synthesis of some biologically relevant beta-C-glycoconjugates

[reaction: see text] An esterification-RCM approach to a variety of biologically relevant beta-C-glycoconjugates is reported herein. A range of carboxylic acids were coupled with several different olefin alcohols 1 to provide esters 3. The esters were then converted to the final ring-closed product 6 in three steps in 49-60% overall yield. The formed compounds are biologically relevant and serve as stable carbohydrate mimics of the corresponding O-glycosides.     

Radical hydroxylation of aromatic hydrocarbons examined by MINDO/3 methods

MINDO/3 calculations have been made on the potential-energy surfaces for the attachment of OH^. radicals to benzene (1) and naphthalene (2) in the vapor state. The activation energies of these reactions are calculated as 88 and 58 kJ/mole. while the enthalpies at 298K are calculated as –211 and –199 kJ/mol. The transition states in (1) and (2) lie closer to the reagents than the products on the reaction coordinate, while (1) has an earlier transition state than does (2). The transition states in these reactions have high dipole moments: 3.1 and 3.6 D, respectively, which are due to charge transfer from the hydrocarbons to the OH^.. Quantum-chemical calculations and kinetic data on the reactions of aromatic hydrocarbons with OH^. in aqueous solution indicate that the mechanism is probably not one involving electron transfer and a rate-limiting stage in the attachment. These processes are of high performance because the radicals are of high stability, while polar effects determine the selectivity.Translated from Teoreticheskaya i Éksperimental'naya Khimiya, Vol. 22, No. 1, pp. 20–26, January–February, 1986.     

Characterization of C5 hydrocarbons relevant to catalysis

A recent in situ infrared study on the selective hydrogenation of C5 dienes and monoenes over a Pd/Al(2)O(3) catalyst only reported incomplete vibrational assignments for some of the reagents, intermediates, and products encountered in that study. This work uses a combination of infrared absorption spectroscopy, Raman, and inelastic neutron scattering to characterize the vibrational spectra of pentane, 1-pentene, cis- and trans-2-pentene, cis- and trans-1,3-pentadiene, 1,4-pentadiene, cyclopentane, and cyclopentene. Ab initio calculations of the potential energy surface, geometry, and vibrational transition energies were performed and simulations of the vibrational spectra compared to the experimental data. Complete vibrational assignments for the majority of the molecules are presented. The potential for using gas-phase infrared measurements for studying heterogeneously catalyzed gas-phase reactions is also briefly considered.     

A supramolecular hydrogel that responds to biologically relevant stimuli

First demonstration of heat and pH-responsive hydrogel of SDS and a zwitterionic amphiphile, sodium N-(n-dodecyl-2-aminoethanoyl)-L-valinate with very low minimum gelation concentration.     

Comparison of linear-scaling semiempirical methods and combined quantum mechanical/molecular mechanical methods for enzymic reactions. II. An energy decomposition analysis

QM/MM methods have been developed as a computationally feasible solution to QM simulation of chemical processes, such as enzyme-catalyzed reactions, within a more approximate MM representation of the condensed-phase environment. However, there has been no independent method for checking the quality of this representation, especially for highly nonisotropic protein environments such as those surrounding enzyme active sites. Hence, the validity of QM/MM methods is largely untested. Here we use the possibility of performing all-QM calculations at the semiempirical PM3 level with a linear-scaling method (MOZYME) to assess the performance of a QM/MM method (PM3/AMBER94 force field). Using two model pathways for the hydride-ion transfer reaction of the enzyme dihydrofolate reductase studied previously (Titmuss et al., Chem Phys Lett 2000, 320, 169-176), we have analyzed the reaction energy contributions (QM, QM/MM, and MM) from the QM/MM results and compared them with analogous-region components calculated via an energy partitioning scheme implemented into MOZYME. This analysis further divided the MOZYME components into Coulomb, resonance and exchange energy terms. For the model in which the MM coordinates are kept fixed during the reaction, we find that the MOZYME and QM/MM total energy profiles agree very well, but that there are significant differences in the energy components. Most significantly there is a large change (approximately 16 kcal/mol) in the MOZYME MM component due to polarization of the MM region surrounding the active site, and which arises mostly from MM atoms close to (<10 A) the active-site QM region, which is not modelled explicitly by our QM/MM method. However, for the model where the MM coordinates are allowed to vary during the reaction, we find large differences in the MOZYME and QM/MM total energy profiles, with a discrepancy of 52 kcal/mol between the relative reaction (product-reactant) energies. This is largely due to a difference in the MM energies of 58 kcal/mol, of which we can attribute approximately 40 kcal/mol to geometry effects in the MM region and the remainder, as before, to MM region polarization. Contrary to the fixed-geometry model, there is no correlation of the MM energy changes with distance from the QM region, nor are they contributed by only a few residues. Overall, the results suggest that merely extending the size of the QM region in the QM/MM calculation is not a universal solution to the MOZYME- and QM/MM-method differences. They also suggest that attaching physical significance to MOZYME Coulomb, resonance and exchange components is problematic. Although we conclude that it would be possible to reparameterize the QM/MM force field to reproduce MOZYME energies, a better way to account for both the effects of the protein environment and known deficiencies in semiempirical methods would be to parameterize the force field based on data from DFT or ab initio QM linear-scaling calculations. Such a force field could be used efficiently in MD simulations to calculate free energies.