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.

A combination of interpolation methods and local saddle-point search algorithms is probably the most efficient way of finding transition states in chemical reactions. Interpolation methods such as the growing-string method and the nudged-elastic band are able to find an approximation to the minimum-energy pathway and thereby provide a good initial guess for a transition state and imaginary mode connecting both reactant and product states. Since interpolation methods employ usually just a small number of configurations and converge slowly close to the minimum-energy pathway, local methods such as partitioned rational function optimization methods using either exact or approximate Hessians or minimum-mode-following methods such as the dimer or the Lanczos method have to be used to converge to the transition state. A modification to the original dimer method proposed by [Henkelman and Jonnson J. Chem. Phys. 111, 7010 (1999)] is presented, reducing the number of gradient calculations per cycle from six to four gradients or three gradients and one energy, and significantly improves the overall performance of the algorithm on quantum-chemical potential-energy surfaces, where forces are subject to numerical noise. A comparison is made between the dimer methods and the well-established partitioned rational function optimization methods for finding transition states after the use of interpolation methods. Results for 24 different small- to medium-sized chemical reactions covering a wide range of structural types demonstrate that the improved dimer method is an efficient alternative saddle-point search algorithm on medium-sized to large systems and is often even able to find transition states when partitioned rational function optimization methods fail to converge.

Our civilization relies on synthetic polymers for all aspects of modern life; yet, inefficient recycling and extremely slow environmental degradation of plastics are causing increasing concern about their widespread use. After a single use, many of these materials are currently treated as waste, underutilizing their inherent chemical and energy value. In this study, energy-rich polyethylene (PE) macromolecules are catalytically transformed into value-added products by hydrogenolysis using well-dispersed Pt nanoparticles (NPs) supported on SrTiO3 perovskite nanocuboids by atomic layer deposition. Pt/SrTiO3 completely converts PE (Mn = 8000–158,000 Da) or a single-use plastic bag (Mn = 31,000 Da) into high-quality liquid products, such as lubricants and waxes, characterized by a narrow distribution of oligomeric chains, at 170 psi H2 and 300 °C under solvent-free conditions for reaction durations up to 96 h. The binding of PE onto the catalyst surface contributes to the number averaged molecular weight (Mn) and the narrow polydispersity (Đ) of the final liquid product. Solid-state nuclear magnetic resonance of 13C-enriched PE adsorption studies and density functional theory computations suggest that PE adsorption is more favorable on Pt sites than that on the SrTiO3 support. Smaller Pt NPs with higher concentrations of undercoordinated Pt sites over-hydrogenolyzed PE to undesired light hydrocarbons.

Overconsumption of single-use plastics is creating a global waste catastrophe, with widespread environmental, economic, and health-related consequences. Inspired by the benefits of processive enzyme-catalyzed conversions of biomacromolecules and guided by spectroscopic interrogations of conformation and dynamics of polymer-surface interactions, we have developed the selective hydrogenolysis of high density polyethylene into a narrow distribution of diesel and lubricant-range alkanes catalyzed by an ordered, mesoporous shell/active site/core catalyst architecture. Solid-state nuclear magnetic resonance investigations of polymer chains adsorbed onto solid materials reveal that an appropriately ordered, porous support orients polymer chains into an all-anti conformation, while measurements of polymer dynamics reveal that long hydrocarbon macromolecules readily move within the pores, with a subsequent escape being inhibited by polymer-surface interactions. These interactions and dynamic behavior resemble the binding and translocation of macromolecules in the catalytic cleft of processive enzymes. Thus, transfer of these features to a mesoporous silica material incorporating catalytic platinum sites for carbon-carbon bond hydrogenolysis of polyethylene provides a reliable stream of alkane products through a processive process.

Interpolation methods such as the nudged elastic band and string methods are widely used for calculating minimum energy pathways and transition states for chemical reactions. Both methods require an initial guess for the reaction pathway. A poorly chosen initial guess can cause slow convergence, convergence to an incorrect pathway, or even failed electronic structure force calculations along the guessed pathway. This paper presents a growing string method that can find minimum energy pathways and transition states without the requirement of an initial guess for the pathway. The growing string begins as two string fragments, one associated with the reactants and the other with the products. Each string fragment is grown separately until the fragments converge. Once the two fragments join, the full string moves toward the minimum energy pathway according to the algorithm for the string method. This paper compares the growing string method to the string method and to the nudged elastic band method using the alanine dipeptide rearrangement as an example. In this example, for which the linearly interpolated guess is far from the minimum energy pathway, the growing string method finds the saddle point with significantly fewer electronic structure force calculations than the string method or the nudged elastic band method.

Reaction pathways and kinetics governing the Rucatalyzed hydrogenation of levulinic acid (LA) in the aqueous phase to form γ-valerolactone (GVL) were considered in a packed bed reactor. GVL can be produced by two distinct hydrogenation pathways; however, over Ru/C at temperatures below 423 K, it forms exclusively via intramolecular esterification of 4-hydroxypentanoic acid (HPA). Over Ru/C, reasonable hydrogenation rates of LA to HPA were observed at near-ambient temperatures (e.g., 0.08 s −1 at 323 K), but GVL selectivities are poor (<5%) under these conditions. Apparent barriers for LA hydrogenation and HPA esterification are 48 and 70 kJ mol −1 , respectively, and GVL selectivity improves at higher temperatures alongside increasing mass transfer limitations in 45−90 μm catalyst particles. Reactivity and selectivity trends in LA hydrogenation below 343 K are well-described by an empirical kinetic model capturing sequential hydrogenation and esterification. Coupling stacked beds of Ru/ C and Amberlyst-15 delivers high GVL yields (∼80%) at near ambient temperatures (323 K) and practical residence times.

In many applications of multilevel/multiscale methods, an active zone must be modeled by a high-level electronic structure method, while a larger environmental zone can be safely modeled by a lower-level electronic structure method, molecular mechanics, or an analytic potential energy function. In some cases though, the active zone must be redefined as a function of simulation time. Examples include a reactive moiety diffusing through a liquid or solid, a dislocation propagating through a material, or solvent molecules in a second coordination sphere (which is environmental) exchanging with solvent molecules in an active first coordination shell. In this article, we present a procedure for combining the levels smoothly and efficiently in such systems in which atoms or groups of atoms move between high-level and low-level zones. The method dynamically partitions the system into the high-level and low-level zones and, unlike previous algorithms, removes all discontinuities in the potential energy and force whenever atoms or groups of atoms cross boundaries and change zones. The new adaptive partitioning (AP) method is compared to Rode's "hot spot" method and Morokuma's "ONIOM-XS" method that were designed for multilevel molecular dynamics (MD) simulations. MD simulations in the microcanonical ensemble show that the AP method conserves both total energy and momentum, while the ONIOM-XS method fails to conserve total energy and the hot spot method fails to conserve both total energy and momentum. Two versions of the AP method are presented, one scaling as O(2N) and one with linear scaling in N, where N is the number of groups in a buffer zone separating the active high-level zone from the environmental low-level zone. The AP method is also extended to systems with multiple high-level zones to allow, for example, the study of ions and counterions in solution using the multilevel approach.

The reaction mechanism for nitrous oxide decomposition has been studied on hydrated and dehydrated mononuclear iron sites in Fe-ZSM-5 using density functional theory. In total, 46 different surface species with different spin states (spin multiplicity M(S) = 4 or 6) and 63 elementary reactions were considered. Heats of adsorption, activation barriers, reaction rates, and minimum energy pathways were determined. The approximate minimum energy pathways and transition states were calculated using the "growing string method" and a modified "dimer method". Spin surface crossing (e.g., O(2) desorption) was considered. The minimum potential energy structure on the seam of two potential energy surfaces was determined with a multiplier penalty function algorithm by Powell and approximate rates of spin surface crossings were calculated. It was found that nitrous oxide decomposition is first order with respect to nitrous oxide concentration and zero order with respect to oxygen concentration. Water impurities in the gas stream have a strong inhibiting effect. In the concentration range of 1-100 ppb, the presence of water vapor influences the surface composition and the apparent rate coefficient. This is especially relevant in the temperature range of 600-700 K where most experimental kinetic studies are performed. Apparent activation barriers determined over this temperature range vary from 28.4 (1 ppb H(2)O) to 54.8 kcal/mol (100 ppb H(2)O). These results give an explanation why different research groups and different catalyst pretreatments often result in very different activation barriers and preexponential factors. Altogether perfect agreement with experimental results could be achieved.