This is supplementary material for an erratum of an article in J. Chem. Phys. It replaces the entire original supplementary material.

Accurate modeling of high-temperature hypersonic flows in the atmosphere requires consideration of collision-induced dissociation of molecular species and energy transfer between the translational and internal modes of the gas molecules. Here, we describe a study of the N2 + N2⟶N2 + 2N and N2 + N2⟶4N nitrogen dissociation reactions using the quasiclassical trajectory (QCT) method. The simulations used a new potential energy surface for the N4 system; the surface is an improved version of one that was presented previously. In the QCT calculations, initial conditions were determined based on a two-temperature model that approximately separates the translational-rotational temperature from the vibrational temperature of the N2 diatoms. Five values from 8000 K to 30,000 K were considered for each of the two temperatures. Over 2.4 × 10(9) trajectories were calculated. We present results for ensemble-averaged dissociation rate constants as functions of the translational-rotational temperature T and the vibrational temperature T(v). The rate constant depends more strongly on T when T(v) is low, and it depends more strongly on T(v) when T is low. Quasibound reactant states contribute significantly to the rate constants, as do exchange processes at higher temperatures. We discuss two sets of runs in detail: an equilibrium test set in which T = T(v) and a nonequilibrium test set in which T(v) < T. In the equilibrium test set, high-v and moderately-low-j molecules contribute most significantly to the overall dissociation rate, and this state specificity becomes stronger as the temperature decreases. Dissociating trajectories tend to result in a major loss of vibrational energy and a minor loss of rotational energy. In the nonequilibrium test set, as T(v) decreases while T is fixed, higher-j molecules contribute more significantly to the dissociation rate, dissociating trajectories tend to result in a greater rotational energy loss, and the dissociation probability's dependence on v weakens. In this way, as T(v) decreases, rotational energy appears to compensate for the decline in average vibrational energy in promoting dissociation. In both the equilibrium and nonequilibrium test sets, in every case, the average total internal energy loss in the dissociating trajectories is between 10.2 and 11.0 eV, slightly larger than the equilibrium potential energy change of N2 dissociation.

Our ability to understand and simulate the reactions catalyzed by iron depends strongly on our ability to predict the relative energetics of spin states. In this work, we studied the electronic structures of Fe ion, gaseous FeO and 14 iron complexes using Kohn-Sham density functional theory with particular focus on determining the ground spin state of these species as well as the magnitudes of relevant spin-state energy splittings. The 14 iron complexes investigated in this work have hexacoordinate geometries of which seven are Fe(ii), five are Fe(iii) and two are Fe(iv) complexes. These are calculated using 20 exchange-correlation functionals. In particular, we use a local spin density approximation (LSDA) - GVWN5, four generalized gradient approximations (GGAs) - BLYP, PBE, OPBE and OLYP, two non-separable gradient approximations (NGAs) - GAM and N12, two meta-GGAs - M06-L and M11-L, a meta-NGA - MN15-L, five hybrid GGAs - B3LYP, B3LYP*, PBE0, B97-3 and SOGGA11-X, four hybrid meta-GGAs - M06, PW6B95, MPW1B95 and M08-SO and a hybrid meta-NGA - MN15. The density functional results are compared to reference data, which include experimental results as well as the results of diffusion Monte Carlo (DMC) calculations and ligand field theory estimates from the literature. For the Fe ion, all functionals except M11-L correctly predict the ground spin state to be quintet. However, quantitatively, most of the functionals are not close to the experimentally determined spin-state splitting energies. For FeO all functionals predict quintet to be the ground spin state. For the 14 iron complexes, the hybrid functionals B3LYP, MPW1B95 and MN15 correctly predict the ground spin state of 13 out of 14 complexes and PW6B95 gets all the 14 complexes right. The local functionals, OPBE, OLYP and M06-L, predict the correct ground spin state for 12 out of 14 complexes. Two of the tested functionals are not recommended to be used for this type of study, in particular M08-SO and M11-L, because M08-SO systematically overstabilizes the high spin state, and M11-L systematically overstabilizes the low spin state.

No abstract

Quantitative simulations of electronically nonadiabatic molecular processes require both accurate dynamics algorithms and accurate electronic structure information. Direct semiclassical nonadiabatic dynamics is expensive due to the high cost of electronic structure calculations, and hence it is limited to small systems, limited ensemble averaging, ultrafast processes, and/or electronic structure methods that are only semiquantitatively accurate. The cost of dynamics calculations can be made manageable if analytic fits are made to the electronic structure data, and such fits are most conveniently carried out in a diabatic representation because the surfaces are smooth and the couplings between states are smooth scalar functions. Diabatic representations, unlike the adiabatic ones produced by most electronic structure methods, are not unique, and finding suitable diabatic representations often involves time-consuming nonsystematic diabatization steps. The biggest drawback of using diabatic bases is that it can require large amounts of effort to perform a globally consistent diabatization, and one of our goals has been to develop methods to do this efficiently and automatically. In this Feature Article, we introduce the mathematical framework of diabatic representations, and we discuss diabatization methods, including adiabatic-to-diabatic transformations and recent progress toward the goal of automatization.

The site-specific basicities of imatinib (Gleevec, a new signal transduction inhibitor drug of chronic myeloid leukemia) and two of its fragment compounds were quantitated in terms of protonation macroconstants, microconstants, and group constants by NMR-pH and pH-potentiometric titrations. Sequential protonation of imatinib follows the N(34), N(11), N(31), N(13) order, in which N(11) and N(31) show commensurable basicity, but negligible intramolecular interaction. Fragment compounds include two "halves" of imatinib, and their moiety-specific basicities confirm the NMR-based protonation sequence of the parent compound. NMR-pH profiles, macro- and/or microscopic protonation schemes, and species-specific distribution diagrams are presented. On the basis of these data, imatinib is shown to be predominantly neutral, monocationic, and tricationic at intestinal, blood, and gastric pH, respectively. The molecular hypotheses on imatinib binding to the Bcr-Abl oncogene fusion protein are interpreted at the site-specific level in view of the moiety basicities of imatinib.

Multiconfiguration pair-density functional theory (MC-PDFT) has previously been applied successfully to carry out ground-state and excited-state calculations. However, because it includes no interaction between electronic states, MC-PDFT calculations in which...

We present a global ground-state triplet potential energy surface for the N2O2 system that is suitable for treating high-energy vibrational-rotational energy transfer and collision-induced dissociation. The surface is based on multi-state complete-active-space second-order perturbation theory/minimally augmented correlation-consistent polarized valence triple-zeta electronic structure calculations plus dynamically scaled external correlation. In the multireference calculations, the active space has 14 electrons in 12 orbitals. The calculations cover nine arrangements corresponding to dissociative diatom-diatom collisions of N2, O2, and nitric oxide (NO), the interaction of a triatomic molecule (N2O and NO2) with the fourth atom, and the interaction of a diatomic molecule with a single atom (i.e., the triatomic subsystems). The global ground-state potential energy surface was obtained by fitting the many-body interaction to 54 889 electronic structure data points with a fitting function that is a permutationally invariant polynomial in terms of bond-order functions of the six interatomic distances.