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 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 correlated Cr₂ dimer, exploring zeolite-catalysed 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.

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.

Nitrogen (N)-doped carbon materials were shown in recent studies to have promising catalytic activity for oxygen reduction reaction (ORR) as a metal-free alternative to platinum, but the underlying molecular mechanism or even the active sites for high catalytic efficiency are still missing or controversial both experimentally and theoretically. We report here the results of periodic density functional theory (DFT) calculations about the ORR at the edge of a graphene nanoribbon (GNR). The edge structure and doped-N near the edge are shown to enhance the oxygen adsorption, the first electron transfer, and also the selectivity toward the four-electron, rather than the two-electron, reduction pathway. We find that the outermost graphitic nitrogen site in particular gives the most desirable characteristics for improved ORR activity, and hence the active site. However, the latter graphitic nitrogen becomes pyridinic-like in the next electron and proton transfer reaction via the ring-opening of a cyclic C–N bond. This inter-conversion between the graphitic and pyridinic sites within a catalytic cycle may reconcile the controversy whether the pyridinic, graphitic, or both nitrogens are active sites.