Jeff Leiding

First Principles Dynamics and Minimum Energy Pathways for Mechanochemical Ring Opening of Cyclobutene

We use ab initio steered molecular dynamics to investigate the mechanically induced ring opening of cyclobutene. We show that the dynamical results can be considered in terms of a force-modified potential energy surface (FMPES). We show how the minimal energy paths for the two possible competing conrotatory and disrotatory ring-opening reactions are affected by external force. We also locate minimal energy pathways in the presence of applied external force and show that the reactant, product, and transition state geometries are altered by the application of external force. The largest effects are on the transition state geometries and barrier heights. Our results provide a framework for future investigations of the role of external force on chemical reactivity.

The First Row Anomaly and Recoupled Pair Bonding in the Halides of the Late p-Block Elements

The dramatic differences between the properties of molecules formed from the late p-block elements of the first row of the periodic table (N–F) and those of the corresponding elements in subsequent rows is well recognized as the first row anomaly. Certain properties of the atoms, such as the relative energies and spatial extents of the ns and np orbitals, can explain some of these differences, but not others. In this Account, we summarize the results of our recent computational studies of the halides of the late p-block elements. Our studies point to a single underlying cause for many of these differences: the ability of the late p-block elements in the second and subsequent rows of the periodic table to form recoupled pair bonds and recoupled pair bond dyads with very electronegative ligands. Recoupled pair bonds form when an electron in a singly occupied ligand orbital recouples the pair of electrons in a doubly occupied lone pair orbital on the central atom, leading to a central atom-ligand bond. Recoupled pair bond dyads occur when a second ligand forms a bond with the orbital left over from the initial recoupled pair bond. Recoupled pair bonds and recoupled pair bond dyads enable the late p-block elements to form remarkably stable hypervalent compounds such as PF5 and SF6 and lead to unexpected excited states in smaller halides of the late p-block elements such as SF and SF2. Recoupled pair bonding also causes the Fn–1X–F bond energies to oscillate dramatically once the normal valences of the central atoms have been satisfied. In addition, recoupled pair bonding provides a lower-energy pathway for inversion in heavily fluorinated compounds (PF3 and PF2H, but not PH2F and PH3) and leads to unusual intermediates and products in reactions involving halogens and late p-block element compounds, such as (CH3)2S + F2. Although this Account focuses on the halides of the second row, late p-block elements, recoupled pair bonds and recoupled pair bond dyads are important in the chemistry of p-block elements beyond the second row (As, Se, and Br) and for compounds of these elements with other very electronegative ligands, such as OH and O. Knowledge of recoupled pair bonding is thus critical to understanding the properties and reactivity of molecules containing the late p-block elements beyond the first row.

Bonding in SCln (n = 1-6): a quantum chemical study.

Following a previous study of bonding and isomerism in the SF(n) and singly chloro-substituted SF(n-1)Cl (n = 1-6) series, we describe bonding in the ground and low-lying excited states of the completely substituted series, SCl(n) (n = 1-6). All structures were characterized at least at the RCCSD(T)/aug-cc-pV(Q+d)Z level of theory. Both differences and similarities were observed between SCl(n) and our previous results on SF(n-1)Cl and SF(n). Several minimum structures that exist in SF(n) and SF(n-1)Cl are absent in SCl(n). For example, the optimized structure of SCl(2)((3)A(2)) is a transition state in C(s) symmetry, whereas the analogous states are minima in SF(n) and SF(n-1)Cl. Second, we found a continuation of a trend discovered in the SF(n-1)Cl series, where Cl substitution has a destabilizing effect that weakens bonds with respect to SF(n). This effect is much stronger in the SCl(n) series than it is in the SF(n-1)Cl series, which is why SCl(2) is the most stable observed species in the family and why SCl(4), SCl(5), and SCl(6) are unstable (SCl(n-2) + Cl(2) additions are endothermic for n = 4-6).

Theoretical studies of the excited doublet states of SF and SCl and singlet states of SF2, SFCl, and SCl2.

In previous work, we reported that the lowest-lying excited states of SF, SCl, SF(2), SFCl, and SCl(2) have recoupled pair bonds. In this study, we examine the analogous low-spin states--the (2)Σ(-) and (2)Δ states of SF and SCl and the excited singlet states of SF(2), SCl(2), and SFCl--which also possess recoupled pair bonds. In contrast to the excited states treated previously, the states studied in the present work have the same spin multiplicities as their respective ground states and are thus potentially observable via electronic excitation. Of particular interest are the minima on the (1)A″ potential energy surface of SFCl corresponding to bond-stretch isomers analogous to those found on the (3)A″ surface. In addition, we discovered that the first two excited states ((1)A″) accessible via vertical excitations from the ground state of SFCl have the electronic structure of the bond-stretch isomers. Thus, electronic excitation spectroscopic studies of SFCl could reveal a signature of the bond-stretch isomers. We will also present limited data on the lowest singlet Rydberg states of the triatomic species. Calculations were performed at the MRCI+Q/aug-cc-pV(Q+d,5+d)Z levels of theory.

Insights into the unusual barrierless reaction between two closed shell molecules, (CH3)2S + F2, and its H2S + F2 analogue: role of recoupled pair bonding.

Early flowtube studies showed that (CH(3))(2)S (DMS) reacted very rapidly with F(2); hydrogen sulfide (H(2)S), however, did not. Recent crossed molecular beam studies found no barrier to the reaction between DMS and F(2) to form CH(2)S(F)CH(3) + HF. At higher collision energies, a second product channel yielding (CH(3))(2)S-F + F was identified. Both reaction channels proceed through an intermediate with an unusual (CH(3))(2)S-F-F bond structure. Curiously, these experimental studies have found no evidence of direct F(2) addition to DMS, resulting in (CH(3))(2)SF(2), despite the fact that the isomer in which both fluorines occupy axial positions is the lowest energy product. We have characterized both reactions, H(2)S + F(2) and DMS + F(2), with high-level ab initio and generalized valence bond calculations. We found that recoupled pair bonding accounts for the structure and stability of the intermediates present in both reactions. Further, all sulfur products possess recoupled pair bonds with CH(2)S(F)CH(3) having an unusual recoupled pair bond dyad involving π bonding. In addition to explaining why DMS reacts readily with F(2) while H(2)S does not, we have studied the pathways for direct F(2) addition to both sulfide species and found that (for (CH(3))(2)S + F(2)) the CH(2)S(F)CH(3) + HF channel dominates the potential energy surface, effectively blocking access to F(2) addition. In the H(2)S + F(2) system, the energy of the transition state for formation of H(2)SF(2) lies very close to the H(2)SF + F asymptote, making the potential pathway a roaming atom mechanism.

Optimization of Semiempirical Quantum Chemistry Methods via Multiobjective Genetic Algorithms: Accurate Photodynamics for Larger Molecules and Longer Time Scales

Excited-state photodynamics is important in numerous varieties of important materials applications (e.g., liquid crystal display, light emitting diode), pharmaceuticals, and chemical manufacturing processing. We study the effectiveness of multiobjective genetic and evolutionary algorithms in multiscaling excited-state direct photodynamics via rapid reparameterization of semiempirical methods. Using a very limited set of ab initio and experimental data, semiempirical parameters are reoptimized to provide globally accurate potential energy surfaces, thereby eliminating the need for expensive ab initio dynamics simulations. Through reoptimization, excited-state energetics are predicted accurately via semiempirical methods, while retaining accurate ground-state predictions. In our initial study of small photo-excited molecules, our results show that the multiobjective evolutionary algorithm consistently yields solutions that are significantly better—up to 384% lower error in the energy and 87% lower error in the energy-gradient—than those reported previously. As verified with direct quantum dynamical calculations, multiple high-quality parameter sets obtained via genetic algorithms show near-ideal behavior on critical and untested excited-state geometries. The results demonstrate that the reparameterization via evolutionary algorithms is a promising way to extend direct dynamics simulations of photochemistry to multi-picosecond time scales and to larger molecules, with promise in more application beyond simple molecular chemistry.

Multiobjective genetic algorithms for multiscaling excited state direct dynamics in photochemistry

This paper studies the effectiveness of multiobjective genetic and evolutionary algorithms in multiscaling excited state direct dynamics in photochemistry via rapid reparameterization of semiempirical methods. Using a very limited set of ab initio and experimental data, semiempirical parameters are reoptimized to provide globally accurate potential energy surfaces, thereby eliminating the need for full-fledged ab initio dynamics simulations, which are very expensive. Through reoptimization of the semiempirical methods, excited-state energetics are predicted accurately, while retaining accurate ground-state predictions. The results show that the multiobjective evolutionary algorithm consistently yields solutions that are significantly better---up to 230% lower error in the energy and 86.5% lower error in the energy-gradient---than those reported in the literature. Multiple high-quality parameter sets are obtained that are verified with quantum dynamical calculations, which show near-ideal behavior on critical and untested excited state geometries. The results demonstrate that the reparameterization strategy via evolutionary algorithms is a promising way to extend direct dynamics simulations of photochemistry to multi-picosecond time scales.