Photodissociation dynamics of the NO dimer. I. Theoretical overview of the ultraviolet singlet excited states

NO (A) Rotational State Distributions from Photodissociation of the N2-NO Complex

We have recorded the resonance-enhanced multiphoton ionization spectrum for NO (A) products from photodissociation of the N₂-NO complex. We made measurements at excitation energies ranging from 28 to 758 cm^-1 above the threshold to produce NO (A) + N₂ (X) products, and the resulting spectra reveal the NO (A) rotational states formed during dissociation, allowing us to determine the rotational state distribution. At the lowest available energies, 28 and 50 cm^-1 above threshold, we observed contributions from NO (A) rotational states that exceed the available energy and must originate from excitation due to hotbands of the complex. At all higher energies, we did not observe any energetically disallowed NO (A) rotational states, and for all available energies above 259 cm^-1 the observed rotational transitions do not extend to the maximum allowed by energy conservation. Furthermore, the observed distributions were typically biased toward low rotational states, in contrast with expectations from vibrational predissociation. From the rotational state distributions, we determined the average fraction of energy partitioned into NO (A) rotation, f_{NO rot, ave}, to be 0.088 at the highest available energy, and this fraction increased as the available energy decreased. By combining the average NO (A) rotational energy along with the average center-of-mass translational energy from our previous work, we determined the average rotational energy for the undetected N₂ (X) photoproduct. The results showed that the N₂ fragment has a higher average rotational energy relative to the NO fragment. Finally, we found that the NO (A) rotational state distribution was colder than expected for a statistical dissociation.

Two Cycling Centers in One Molecule: Communication by Through-Bond Interactions and Entanglement of the Unpaired Electrons

Many applications in quantum information science (QIS) rely on the ability to laser-cool molecules. The scope of applications can be expanded if laser-coolable molecules possess two or more cycling centers, i.e., moieties capable of scattering photons via multiple absorption-emission events. Here we employ the equation-of-motion coupled-cluster method for double electron attachment (EOM-DEA-CCSD) to study the electronic structure of hypermetallic molecules with two alkaline-earth metals connected by an acetylene linker. The electronic structure of the molecules is similar to that of two separated alkali metals; however, the interaction between the two electrons is weak and largely dominated by through-bond interactions. The communication between the two cycling centers is quantified by the extent of the entanglement of the two unpaired electrons associated with the two cycling centers. This contribution highlights the rich electronic structure of hypermetallic molecules that may advance various applications in QIS and beyond.

The gas-phase structure of the asymmetric, trans-dinitrogen tetroxide (N2O4), formed by dimerization of nitrogen dioxide (NO2), from rotational spectroscopy and ab initio quantum chemistry

We report the first experimental gas-phase observation of an asymmetric, trans-N₂O₄ formed by the dimerization of NO₂. In additional to the dominant ¹⁴N₂¹⁶O₄ species, rotational transitions have been observed for all species with single ¹⁵N and ¹⁸O substitutions as well as several multiply substituted isotopologues. These transitions were used to determine a complete substitution structure as well as an r₀ structure from the fitted zero-point averaged rotational constants. The determined structure is found to be that of an ON-O-NO₂ linkage with the shared oxygen atom closer to the NO₂ than the NO (1.42 vs 1.61 Å). The structure is found to be nearly planar with a trans O-N-O-N linkage. From the spectra of the ¹⁴N¹⁵NO₄ species, we were able to determine the nuclear quadrupole coupling constants for each specific nitrogen atom. The equilibrium structure determined by ab initio quantum chemistry calculations is in excellent agreement with the experimentally determined structure. No spectral evidence of the predicted asymmetric, cis-N₂O₄ was found in the spectra.

Ultrafast Dissociation of Metastable CO^{2+} in a Dimer

We triply ionize the van der Waals bound carbon monoxide dimer with intense ultrashort pulses and study the breakup channel (CO)_{2}^{3+}→C^{+}+O^{+}+CO^{+}. The fragments are recorded in a cold target recoil ion momentum spectrometer. We observe a fast CO^{2+} dissociation channel in the dimer, which does not exist for the monomer. We found that a nearby charge breaks the symmetry of a X^{3}Π state of CO^{2+} and induces an avoided crossing that allows a fast dissociation. Calculation on the full dimer complex shows the coupling of different charge states, as predicted from excimer theory, gives rise to electronic state components not present in the monomer, thereby enabling fast dissociation with higher kinetic energy release. These results demonstrate that the electronic structure of molecular cluster complexes can give rise to dynamics that is qualitatively different from that observed in the component monomers.

Ab initio and semi-empirical Molecular Dynamics simulations of chemical reactions in isolated molecules and in clusters

Recent progress in “on-the-fly” trajectory simulations of molecular reactions, using different electronic structure methods is discussed, with analysis of the insights that such calculations can provide and of the strengths and limitations of the algorithms available.

Comparison of the performance of exact-exchange-based density functional methods

How to describe nondynamic electron correlation is still a major challenge to density functional theory (DFT). Recent models designed particularly for this problem, such as Becke'05 (B05) and Perdew-Staroverov-Tao-Scuseria (PSTS) functionals employ the exact-exchange density, the efficient calculation of which is technically quite challenging. We have recently implemented self-consistently the B05 functional based on an efficient resolution-identity (RI) technique. In this study, we report a self-consistent RI implementation of the PSTS functional. In contrast to its original implementation, our version brings no limitation on the choice of the basis set. We have also implemented the Mori-Sanchez-Cohen-Yang-2 (MCY2) functional, another recent DFT method that includes full exact exchange. The performance of PSTS, B05, and MCY2 is validated on thermochemistry, reaction barriers, and dissociation energy curves, with an emphasis on nondynamic correlation effects in the discussion. All three methods perform rather well in general, B05 and MCY2 being on average somewhat better than PSTS. We include also results with other functionals that represent various aspects of the development in this field in recent years, including B3LYP, M06-HF, M06-2X, ωB97X, and TPSSh. The performance of the heavy-parameterized functionals M06-2X and ωB97X is on average better than that of B05, MCY2, and PSTS for standard thermodynamic properties and reactions, while the latter functionals do better in hydrogen abstraction reactions and dissociation processes. In particular, B05 is found to be the only functional that yields qualitatively correct dissociation curves for two-center symmetric radicals like He(2)(+). Finally, we compare the performance of all these functionals on a strongly correlated exemplary case system, the NO dimer. Only PSTS, B05, and MCY2 describe the system qualitatively correctly. Overall, this new type of functionals show good promise of overcoming some of the difficulties DFT encounters for systems with strong nondynamic correlation.

Improved self-consistent and resolution-of-identity approximated Becke'05 density functional model of nondynamic electron correlation

In a recent letter [E. Proynov, Y. Shao, and J. Kong, Chem. Phys. Lett. 493, 381 (2010)], Becke's B05 model of nondynamic electron correlation in density functional theory was implemented self-consistently with computational efficiency (the "SCF-RI-B05" scheme). Important modifications of the algorithm were done in order to make the self-consistency feasible. In the present work, we give a complete account of the SCF-RI-B05 algorithm, including all the formulae for the analytical representation of the B05 functional and for its self-consistent field (SCF) potential. The average performance of the SCF-RI-B05 method reported in the above letter was somewhat less accurate, compared to the original B05 implementation, mainly because the parameters of the original B05 model were optimized with post-local-spin-density calculations. In this work, we report improved atomization energies with SCF-RI-B05, based on a SCF re-optimization of its four linear parameters. The re-optimized SCF-RI-B05 scheme is validated also on reaction barriers, and on the subtle energetics of NO dimer, an exemplary system of strong nondynamic correlation. It yields both the binding energy and the singlet-triplet splitting of the NO dimer correctly, and close to the benchmarks reported in the literature.

Experimental investigation of the absolute enthalpies of formation of 2,3-, 2,4-, and 3,4-pyridynes

The absolute enthalpies of formation of 3,4-, 2,3-, and/or 2,4-didehydropyridines (3,4-, 2,3- and 2,4-pyridynes) have been determined by using energy-resolved collision-induced dissociation of deprotonated 2- and 3-chloropyridines. Bracketing experiments find the gas-phase acidities of 2- and 3-chloropyridines to be 383 ± 2 and 378 ± 2 kcal/mol, respectively. Whereas deprotonation of 3-chloropyridine leads to formation of a single ion isomer, deprotonation of the 2-chloro isomer results in a nearly 60:40 mixture of regioisomers. The enthalpy of formation of 3,4-pyridyne is measured to be 121 ± 3 kcal/mol by using the chloride dissociation energy for deprotonated 3-chloropyridine. The structure of the product formed upon dissociation of the ion from 2-chloropyridine cannot be unequivocally assigned because of the isomeric mixture of reactant ions and the fact that the potential neutral products (2,3-pyridyne and 2,4-pyridyne) are predicted by high level spin-flip coupled-cluster calculations to be nearly the same in energy. Consequently, the enthalpies of formation for both neutral products are assigned to be 130 ± 3 kcal/mol. Comparison of the enthalpies of dehydrogenation of benzene and pyridine indicates that the nitrogen in the pyridine ring does not have any effect on the stability of the aryne triple bond in 3,4-pyridyne, destabilizes the aryne triple bond in 2,3-pyridyne, and stabilizes the 1,3-interaction in 2,4-pyridyne compared to that in m-benzyne. Natural bond order calculations show that the effects on the 2,3- and 2,4-pyridynes result from polarization of the electrons caused by interaction with the lone pair. The polarization in 2,4-pyridyne is stabilizing because it creates a 1,2-interaction between the nitrogen and dehydrocarbons that is stronger than the 1,3-interaction between the dehydrocarbons.

Photoelectron imaging of NCCCN(-): The triplet ground state and the singlet-triplet splitting of dicyanocarbene

The photoelectron spectra of NCCCN(-) have been measured at 355 and 266 nm by means of photoelectron imaging. The spectra show two distinct features, corresponding to the ground and first excited states of dycianocarbene. With support from theoretical calculations using the spin-flip coupled-cluster methods, the ground electronic state of HCCCN is assigned as a triplet state, while the first excited state is a closed-shell singlet. The photoelectron band corresponding to the triplet is broad and congested, indicating a large geometry change between the anion and neutral. A single sharp feature of the singlet band suggests that the geometry of the excited neutral is similar to that of the anion. In agreement with these observations, theoretical calculations show that the neutral triplet state is either linear or quasilinear (X (3)B(1) or (3)Sigma(g) (-)), while the closed-shell singlet (a (1)A(1)) geometry is strongly bent, similar to the anion structure. The adiabatic electron binding energy of the closed-shell singlet is measured to be 3.72+/-0.02 eV. The best estimate of the origin of the triplet band gives an experimental upper bound of the adiabatic electron affinity of NCCCN, EA</=3.25+/-0.05 eV, while the Franck-Condon modeling yields an estimate of EA(NCCCN)=3.20+/-0.05 eV. From these results, the singlet-triplet splitting is estimated to be DeltaE(ST)(X (3)B(1)/(3)Sigma(g) (-)-a (1)A(1))=0.52+/-0.05 eV (12.0+/-1.2 kcal/mol).

Efficient self-consistent DFT calculation of nondynamic correlation based on the B05 method

Becke's B05 method of describing nondynamic electron correlation in Density Functional Theory is implemented self-consistently with computational efficiency. Important modifications of the method are proposed in order to make the self-consistency feasible. Resolution-of-identity technique is used to reduce dramatically the cost of the required exact-exchange energy density. The method is briefly validated on a variety of properties. It describes accurately for the first time the subtle energetics of the NO dimer, an exemplary system of strong nondynamic correlation. The efficient algorithm for the exact-exchange energy density can be applied to other functionals that use this quantity.

Reactivity of the NO dimer: on the role of the triplet electronic states

We performed ab initio calculations to investigate the lowest triplet and singlet electronic states of the NO dimer and their mutual spin-orbit couplings. The electronic structure calculations are done using multiconfigurational approaches and a large diffuse basis set. A high density of electronic states is found favoring their mutual interactions by vibronic and spin-orbit couplings. We used our potential curves and spin-orbit couplings to discuss the mechanisms for the IR dissociation and the UV photodissociation of the NO dimer and the electronic de-excitation of NO (A(2)Sigma(+)) after collision with NO (X(2)Pi). For these reactions, multistep pathways, which involve the long-range and the molecular regions of the potential energy surfaces of the triplet and singlet electronic states of N(2)O(2), are suggested. A qualitative agreement between our findings and previous experimental assumptions is found.

Photofragment Spectroscopy and Predissociation Dynamics of Weakly Bound Molecules

Photofragment spectroscopy is combined with imaging techniques and time-resolved measurements of photoions and photoelectrons to explore the predissociation dynamics of weakly bound molecules. Recent experimental advances include measurements of pair-correlated distributions, in which energy disposal in one cofragment is correlated with a state-selected level of the other fragment, and femtosecond pump-probe experiments, in some cases with coincidence detection. An application in which coincident measurements are carried out in the molecular frame is also described. To illustrate these state-selective and time-resolved techniques, we review two recent applications: (a) the photoinitiated dissociation of the covalently bound NO dimer on the ground and excited electronic states and the role of state couplings and (b) the state-selected vibrational predissociation of hydrogen-bonded acetylene dimers with HCl (acid) and ammonia (base) and the importance of angular momentum constraints. We highlight the crucial role of theoretical models in interpreting results.

Equation-of-motion coupled-cluster methods for open-shell and electronically excited species: the Hitchhiker's guide to Fock space

The equation-of-motion coupled-cluster (EOM-CC) approach is a versatile electronic-structure tool that allows one to describe a variety of multiconfigurational wave functions within single-reference formalism. This review provides a guide to established EOM methods illustrated by examples that demonstrate the types of target states currently accessible by EOM. It focuses on applications of EOM-CC to electronically excited and open-shell species. The examples emphasize EOM's advantages for selected situations often perceived as multireference cases [e.g., interacting states of different nature, Jahn-Teller (JT) and pseudo-JT states, dense manifolds of ionized states, diradicals, and triradicals]. I also discuss limitations and caveats and offer practical solutions to some problematic situations. The review also touches on some formal aspects of the theory and important current developments.

Photochemistry studied with ab initio orbital-correlation and state-correlation plots: Classic cyclobutene ring opening, and the reaction of N2 with photoexcited O2

Pericyclic reaction theory arose from ideas presented in 1965, based on orbital-energy correlation diagrams (Woodward and Hoffmann) and state-energy correlation diagrams (Longuet-Higgins and Abrahamson). Here we have used ab initio complete-active-space self-consistent field (CASSCF) calculations to generate such diagrams. First we present diagrams for the classic case of cyclobutene ring opening, to demonstrate agreement between the CASSCF results and the classic diagrams of both Woodward/Hoffmann and Longuet-Higgins/Abrahamson. Then we present diagrams for the more difficult cases of N(2) + photoexcited O(2), to produce either 2 NO or NNO + O. These N(2) + O(2) cases feature significant electron reorganization, for which elementary pencil-and-paper diagrams are less accurate. Finally, the benefits and limitations of such diagrams for predicting photochemistry are briefly discussed.