A summary of the technical advances that are incorporated in the fourth major release of the Q-Chem quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and openshell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Møller-Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Some specific chemical examples include calculations on the strongly correlated Cr 2 dimer, exploring zeolitecatalysed ethane dehydrogenation, energy decomposition analysis of a charged ter-molecular complex arising from glycerol photoionisation, and natural transition orbitals for a Frenkel exciton state in a nine-unit model of a self-assembling nanotube.Keywords quantum chemistry, software, electronic structure theory, density functional theory, electron correlation, computational modelling, Q-Chem Disciplines Chemistry CommentsThis article is from Molecular Physics: An International Journal at the Interface Between Chemistry and Physics 113 (2015): 184, doi:10.1080/00268976.2014. RightsWorks produced by employees of the U.S. Government as part of their official duties are not copyrighted within the U.S. The content of this document is not copyrighted. Authors 185A summary of the technical advances that are incorporated in the fourth major release of the Q-CHEM quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and open-shell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Møller-Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Some specific chemical examples include calculations on the strongly corre...

Advances in theory and algorithms for electronic structure calculations must be incorporated into program packages to enable them to become routinely used by the broader chemical community. This work reviews advances made over the past five years or so that constitute the major improvements contained in a new release of the Q-Chem quantum chemistry package, together with illustrative timings and applications. Specific developments discussed include fast methods for density functional theory calculations, linear scaling evaluation of energies, NMR chemical shifts and electric properties, fast auxiliary basis function methods for correlated energies and gradients, equation-of-motion coupled cluster methods for ground and excited states, geminal wavefunctions, embedding methods and techniques for exploring potential energy surfaces.

High-level ab initio molecular orbital calculations are used to study the thermodynamics and electrochemistry relevant to the mechanism of atom transfer radical polymerization (ATRP). Homolytic bond dissociation energies (BDEs) and standard reduction potentials (SRPs) are reported for a series of alkyl halides (R-X; R = CH 2CN, CH(CH 3)CN, C(CH 3) 2CN, CH 2COOC 2H 5, CH(CH 3)COOCH 3, C(CH 3) 2COOCH 3, C(CH 3) 2COOC 2H 5, CH 2Ph, CH(CH 3)Ph, CH(CH 3)Cl, CH(CH 3)OCOCH 3, CH(Ph)COOCH 3, SO 2Ph, Ph; X = Cl, Br, I) both in the gas phase and in two common organic solvents, acetonitrile and dimethylformamide. The SRPs of the corresponding alkyl radicals, R (*), are also examined. The computational results are in a very good agreement with the experimental data. For all alkyl halides examined, it is found that, in the solution phase, one-electron reduction results in the fragmentation of the R-X bond to the corresponding alkyl radical and halide anion; hence it may be concluded that a hypothetical outer-sphere electron transfer (OSET) in ATRP should occur via concerted dissociative electron transfer rather than a two-step process with radical anion intermediates. Both the homolytic and heterolytic reactions are favored by electron-withdrawing substituents and/or those that stabilize the product alkyl radical, which explains why monomers such as acrylonitrile and styrene require less active ATRP catalysts than vinyl chloride and vinyl acetate. The rate constant of the hypothetical OSET reaction between bromoacetonitrile and Cu (I)/TPMA complex was estimated using Marcus theory for the electron-transfer processes. The estimated rate constant k OSET = approximately 10 (-11) M (-1) s (-1) is significantly smaller than the experimentally measured activation rate constant ( k ISET = approximately 82 M (-1) s (-1) at 25 degrees C in acetonitrile) for the concerted atom transfer mechanism (inner-sphere electron transfer, ISET), implying that the ISET mechanism is preferred. For monomers bearing electron-withdrawing groups, the one-electron reduction of the propagating alkyl radical to the carbanion is thermodynamically and kinetically favored over the one-electron reduction of the corresponding alkyl halide unless the monomer bears strong radical-stabilizing groups. Thus, for monomers such as acrylates, catalysts favoring ISET over OSET are required in order to avoid chain-breaking side reactions.

Standard reduction potentials, SRPs, of the halogen atoms have been calculated in water on the basis of an appropriate thermochemical cycle. Using the best up-to-date thermodynamic data available in the literature, we have calculated E(o)(X•/X-) values of 3.66, 2.59, 2.04, and 1.37 V vs SHE for F•, Cl•, Br•, and I•, respectively. Additionally, we have computed the SRPs of Cl•, Br•, and I• in acetonitrile (CH3CN) and dimethylformamide (DMF) by correcting the values obtained in water for the free energies of transfer of X• and X- from water to the nonaqueous solvent S and the intersolvent potential between water and S. From the values of E(o)(X •/X-) in CH(3)CN and DMF, the SRPs of a series of alkyl halides of relevance to atom transfer radical polymerization and other important processes such as pollution abatement have been calculated in these two solvents. This has been done with the aid of a thermochemical cycle involving the gas-phase homolytic dissociation of the C-X bond, solvation of RX, R•, and X•, and reduction of X• to X- in solution.

High-level ab initio molecular orbital theory is used to obtain benchmark values for the ferricenium/ferrocene (Fc(+)/Fc) couple, the IUPAC recommended reference electrode for nonaqueous solution. The gas-phase ionization energy of ferrocene is calculated using the high-level composite method, G3(MP2)-RAD, and two higher-level variants of this method. These latter methods incorporate corrections for core correlation and, in the case of the highest level considered, use (RO)CCSD(T)/6-311+G(d,p) in place of (RO)CCSD(T)/6-31G(d) as the base level of theory. All methods provide good agreement with one another and the corresponding experimental values. Solvation energies have been calculated using PCM, CPCM, SMD, and COSMO-RS. Using G3(MP2)-RAD-Full-TZ gas-phase energies and COSMO-RS solvation energies, the absolute redox potentials of the Fc(+)/Fc couple have been calculated as 4.988, 4.927, and 5.043 V in acetonitrile, 1,2-dichloroethane, and dimethylsulfoxide solutions, respectively.

The performance of a variety of DFT functionals (BLYP, PBE, B3LYP, B3P86, KMLYP, B1B95, MPWPW91, MPW1B95, BB1K, MPW1K, MPWB1K, and BMK), together with the ab initio methods RHF, RMP2, and G3(MP2)-RAD, and with ONIOM methods based on combinations of these procedures, is examined for calculating the enthalpies of a range of radical reactions. The systems studied include the bond dissociation energies (BDEs) of R-X (R = CH3, CH2F, CH2OH, CH2CN, CH2Ph, CH(CH3)Ph, C(CH3)2Ph; X = H, CH3, OCH3, OH, F), RCH(Ph)-X (R = CH3, CH3CH2, CH(CH3)2, C(CH3)3, CH2F, CH2OH, CH2CN; X = H, F), R-TEMPO (R = CH3, CH2CH3, CH(CH3)2, C(CH3)3, CH2CH2CH3, CH2F, CH2OH, CH2CN, CH(CN)CH3, CH(Cl)CH3; TEMPO = 2,2,6,6,-tetramethylpiperidin-1-yloxyl) and HM1M2-X (M1, M2 = CH2CH(CH3), CH2CH(COOCH3), CH2C(CH3)(COOCH3); X = Cl, Br), the beta-scission energies of RXCH2* and RCH2CHPh* (R = CH3, CH2CH3, CH(CH3)2, C(CH3)3; X = O, S, CH2), and the enthalpies of several radical addition, ring-opening, and hydrogen- and chlorine-transfer reactions. All of the DFT methods examined failed to provide an accurate description of the energetics of the radical reactions when compared with benchmark G3(MP2)-RAD values, with all methods tested showing unpredictable deviations of up to 40 kJ mol-1 or more in some cases. RMP2 also shows large deviations from G3(MP2)-RAD in the absolute values of the enthalpies of some types of reaction and, although it fares somewhat better than the DFT methods in modeling the relative values, it fails for substituents capable of strongly interacting with the unpaired electron. However, it is possible to obtain cost-effective accurate calculations for radical reactions using ONIOM-based procedures in which a high-level method, such as G3(MP2)-RAD, is only used to model the core reaction (which should contain all substituents alpha to the reaction center), and the full system is modeled using a lower-cost procedure such as RMP2.

Alpha-, beta-, gamma-, and delta-tocopherols can be oxidized in dry CH2Cl2 or CH3CN by one electron to form cation radicals that deprotonate to form the neutral phenoxyl radicals, which are then immediately further oxidized by one electron to the phenoxonium cations (an ECE electrochemical mechanism, where E signifies an electron transfer and C represents a chemical step, with the electrochemical mechanism having been determined by in situ spectroscopic analysis). The principal difference in the electrochemical behavior of the tocopherols relates to the stability of their associated phenoxonium cations. The phenoxonium cation of alpha-tocopherol is stable in solution for at least several hours, the phenoxonium cation of beta-tocopherol is stable for several minutes, and the phenoxonium cations of gamma- and delta-tocopherol are stable for <1 s. In dry CH2Cl2 containing >0.75 M acid (CF3COOH), the deprotonation reaction of the cation radicals can be completely inhibited resulting in the cyclic voltammetric behavior of the tocopherols appearing as chemically reversible one-electron oxidation processes (an E mechanism). In dry acid conditions, the cation radicals can be further oxidized by one electron to form the dications, which are unstable and immediately deprotonate. The high stability of the phenoxonium cation of alpha-tocopherol compared to the other tocopherols (and most other phenols) is a chemically important feature that may shed new light on understanding alpha-tocopherol's unique biological properties.

The mechanism of reductive cleavage of model alkyl halides (methyl 2-bromoisobutyrate, methyl 2-bromopropionate, and 1-bromo-1-chloroethane), used as initiators in living radical polymerization (LRP), has been investigated in acetonitrile using both experimental and computational methods. Both theoretical and experimental investigations have revealed that dissociative electron transfer to these alkyl halides proceeds exclusively via a concerted rather than stepwise manner. The reductive cleavage of all three alkyl halides requires a substantial activation barrier stemming mainly from the breaking C-X bond. The activation step during single electron transfer LRP (SET-LRP) was originally proposed to proceed via formation and decomposition of RX(•-) through an outer sphere electron transfer (OSET) process (Guliashvili, T.; Percec, V. J. Polym. Sci., Part A: Polym. Chem. 2007, 45, 1607). These radical anion intermediates were proposed to decompose via heterolytic rather than homolytic C-X bond dissociation. Here it is presented that injection of one electron into RX produces only a weakly associated charge-induced donor-acceptor type radical anion complex without any significant covalent σ type bond character between carbon-centered radical and associated anion leaving group. Therefore, neither homolytic nor heterolytic bond dissociation applies to the reductive cleavage of C-X in these alkyl halides inasmuch as a true radical anion does not form in the process. In addition, the whole mechanism of SET-LRP has to be revisited since it is based on presumed OSET involving intermediate RX(•-), which is shown here to be nonexistent.