Global Potential Energy Surfaces by Compressed-State Multistate Pair-Density Functional Theory: The Lowest Doublet States Responsible for the N(4Su)+C2(a 3 Pi u) -> CN(X 2 Sigma+) + C(3Pg) Reaction

Global potential energy surfaces (PESs) for the 1 2A ' and 1 2A '' states of the C2N system responsible for the N(4Su) + C2(a 3 pi u) -> CN(X 2 sigma+) + C(3Pg) reaction are mapped using compressed-state multistate pair-density functional theory (CMSPDFT), which is a multi-state version of multiconfiguration pair density functional theory (MC-PDFT). Calculations are also performed at selected geometries by explicitly correlated multi reference configuration interaction with quadruple corrections, MRCI-F12+Q, and the comparison of the two sets of calculations shows that CMS-PDFT describes the globally reactive PESs well, including the bond-breaking asymptotes. We conclude that CMS-PDFT is an efficient method for constructing PESs for strongly correlated reactive systems. The PESs for producing CN + C are found to be barrierless and proceed through intermediate complexes. The CMS-PDFT PESs were fitted with a neural network method, and quasiclassical trajectories were computed on the resulting analytic PESs. These trajectories predict that the reaction produces vibrationally excited CN.