Theoretical Approaches for Understanding the Interplay Between Stress and Chemical Reactivity.

The use of mechanical stresses to induce chemical reactions has attracted significant interest in recent years. Computational modeling can play a significant role in developing a comprehensive understanding of the interplay between stresses and chemical reactivity. In this review, we discuss techniques for simulating chemical reactions occurring under mechanochemical conditions. The methods described are broadly divided into techniques that are appropriate for studying molecular mechanochemistry and those suited to modeling bulk mechanochemistry. In both cases, several different approaches are described and compared. Methods for examining molecular mechanochemistry are based on exploring the force-modified potential energy surface on which a molecule subjected to an external force moves. Meanwhile, it is suggested that condensed phase simulation methods typically used to study tribochemical reactions, i.e., those occurring in sliding contacts, can be adapted to study bulk mechanochemistry.

Dramatic substituent effects on the mechanisms of nucleophilic attack on Se-S bridges.

The reactions of XSeSX, XSeSY, and YSeSX (X, Y = CH3, NH2, OH, F) with F(-) and CN(-) nucleophiles have been investigated by means of B3PW91/6-311+G(2df,p) and G4 calculations. In systems where the two substituents are not identical (XSeSY), the more stable of the two possible isomers corresponds to those in which the most electronegative substituent is attached to Se. Nucleophilic attack takes place at Se, independent of the nature of the nucleophile, with the only exception being XSeSF (X = CH3 , NH2 , OH), in which case the attack occurs at S. In agreement with recent results for disulfide and diselenide linkages, the mechanisms leading to Se-S bond cleavage are not always the more favorable ones because for highly electronegative substituents the most favorable process is fission of the chalcogen-substituent bond. These dissimilarities in the observed reactivity pattern as a function of the electronegativity of the substituents are due to the fact that the σ-type Se-S antibonding orbital, which for low-electronegative substituents is the lowest unnoccupied molecular orbital (LUMO), becomes strongly destabilized when the electronegativity of the substituent increases, and is replaced by an antibonding π-type Se-X (or S-X) orbital. In contrast, however, with what has been found for disulfide and diselenide derivatives, the observed reactivity does not change with the nature of the nucleophile. The activation strain model provides interesting insight into these processes, showing that in most cases the activation barriers are the consequence of subtle differences in the strain or in the interaction energies.

Revealing unexpected mechanisms for nucleophilic attack on S-S and Se-Se bridges.

The reactivity of disulfide and diselenide derivatives towards F(-) and CN(-) nucleophiles has been investigated by means of B3PW91/6-311+G(2df,p) calculations. This theoretical survey shows that these processes, in contrast with the generally accepted view of disulfide and diselenide linkages, do not always lead to SS or SeSe bond cleavage. In fact, SS or SeSe bond fission is the most favorable process only when the substituents attached to the S or the Se atoms are not very electronegative. Highly electronegative substituents (X) strongly favor SX bond fission. This significant difference in the observed reactivity patterns is directly related to the change in the nature of the LUMO orbital of the disulfide or diselenide derivative as the electronegativity of the substituents increases. For weakly electronegative substituents, the LUMO is a σ-type SS (or SeSe) antibonding orbital, but as the electronegativity of the substituents increases the π-type SX antibonding orbital stabilizes and becomes the LUMO. The observed reactivity also changes with the nature of the nucleophile and with the S or Se atom that undergoes the nucleophilic attack in asymmetric disulfides and diselenides. The activation strain model provides interesting insights into these processes. There are significant similarities between the reactivity of disulfides and diselenides, although some dissimilarities are also observed, usually related to the different interaction energies between the fragments produced in the fragmentation process.

Mechanism of the Reduction of an Oxidized Glutathione Peroxidase Mimic with Thiols.

N,N-dimethylbenzylamine-2-selenol is a well-known, efficient glutathione peroxidase mimic. This compound reduces peroxides through a three-step catalytic mechanism, of which the first step has been well-characterized computationally. The mechanism for the reaction of N,N-dimethylbenzylamine-2-selenenic acid with a thiol, the second step in the catalytic cycle, is studied using reliable electronic structure techniques. Two different mechanisms are identified, using either a thiol or a deprotonated thiolate as the nucleophile. It is found that the lowest energy barrier mechanism incorporates two explicit solvent molecules to shuttle the thiol hydrogen to the leaving hydroxyl group, while the alternative mechanism using the thiolate has a barrier four times higher.

Substituent-controlled reactivity in the Nazarov cyclisation of allenyl vinyl ketones.

Alkyl substitution α to the ketone of an allenyl vinyl ketone enhances Nazarov reactivity by inhibiting alternative pathways involving the allene moiety and through electron donation and/or steric hindrance. This substitution pattern also accelerates Nazarov cyclisation by increasing the population of the reactive conformer and by stabilising the oxyallyl cation intermediate. Furthermore, α substitution by an alkyl group does not alter the regioselectivity of interrupted Nazarov reactions when the oxyallyl cation intermediate is intercepted by addition of an oxygen nucleophile, or by [4+3] cyclisation with acyclic dienes. The regioselectivity of the Nazarov process for allenyl vinyl ketones was determined to be a result of an electronic bias in the oxyallyl cation intermediate. Computational data are consistent with this observation.

Systematic study of the performance of density functional theory methods for prediction of energies and geometries of organoselenium compounds.

A variety of density functional theory (DFT) methods are paired with Pople basis sets of varying sizes and evaluated for use with organoselenium compounds. The ability of each method to predict reliable geometries and energies is determined through comparison with quadratic configuration interaction with single and double excitations (QCISD) results. The recommended procedure for accurate prediction of energies and geometries is to use the B3PW91 functional with the 6-311G(2df,p) basis set. The B3PW91/6-31G(d,p) level of theory gives almost identical geometries as larger basis sets, so geometries can be predicted at this level for computational efficiency.

Theoretical investigations on the reaction of monosubstituted tertiary-benzylamine selenols with hydrogen peroxide.

The effects of introducing electron-donating (NH(2), OCH(3), CH(3)) and electron-withdrawing (NO(2), CF(3), CN, F) groups to N,N-dimethylbenzylamine-2-selenol are studied to determine the effect of the selenium electron density on the efficiency of the reduction of hydrogen peroxide. Introducing substituents in the meta and para positions decreases or increases the energy barrier of the reaction in the expected way, due to changes in the electronic environment of the reacting selenium center. Ortho substituents are found to have a greater effect on the electronic environment of the selenium center, which is mitigated by changing the steric environment. Insight into the origins of the substituent effects is obtained through quantum theory of atoms in molecules (QTAIM) and electrostatic potential analysis.

Reduction of hydrogen peroxide by glutathione peroxidase mimics: reaction mechanism and energetics.

The reaction mechanism for the reduction of hydrogen peroxide by N,N-dimethylbenzylamine diselenide, its selenol analogue, and the charged analogues of the diselenide and selenol are elucidated using reliable electronic structure techniques. It is found that the reaction using the diselenide has a large Gibbs energy barrier of 173.5 kJ/mol. The cationic diselenide, with both amines protonated, shows a lower barrier of 103.5 kJ/mol. Both diselenide species show significant Se-Se bond lengthening upon oxidation. An unusual two-step mechanism is found for the selenol with barriers of 136.3 and 141.9 kJ/mol, respectively, showing that it is unlikely that the selenol is the active form. The zwitterion, selenolate, and protonated amine analogues of the selenol show one-step reactions with energy barriers of 82.7, 92.7, and 102.3 kJ/mol, respectively. The zwitterion of the selenol shows the most favorable reaction energies, which is in good agreement with proposed mechanisms for this reaction.