Nicholas J. Mosey

Rotationally invariant ab initio evaluation of Coulomb and exchange parameters for DFT+U calculations

Conventional density functional theory (DFT) fails for strongly correlated electron systems due to large intra-atomic self-interaction errors. The DFT+U method provides a means of overcoming these errors through the use of a parametrized potential that employs an exact treatment of quantum mechanical exchange interactions. The parameters that enter into this potential correspond to the spherically averaged intra-atomic Coulomb (U) and exchange (J) interactions. Recently, we developed an ab initio approach for evaluating these parameters on the basis of unrestricted Hartree-Fock (UHF) theory, which has the advantage of being free of self-interaction errors and does not require experimental input [Mosey and Carter, Phys. Rev. B 76, 155123 (2007)]. In this work, we build on that method to develop a more robust and convenient ab initio approach for evaluating U and J. The new technique employs a relationship between U and J and the Coulomb and exchange integrals evaluated using the entire set of UHF molecular orbitals (MOs) for the system. Employing the entire set of UHF MOs renders the method rotationally invariant and eliminates the difficulty in selecting unambiguously the MOs that correspond to localized states. These aspects overcome two significant deficiencies of our earlier method. The new technique is used to evaluate U and J for Cr(2)O(3), FeO, and Fe(2)O(3). The resulting values of U-J are close to empirical estimates of this quantity for each of these materials and are also similar to results of constrained DFT calculations. DFT+U calculations using the ab initio parameters yield results that are in good agreement with experiment. As such, this method offers a means of performing accurate and fully predictive DFT+U calculations of strongly correlated electron materials.

A generalization of the charge equilibration method for nonmetallic materials.

Assigning effective atomic charges that properly reproduce the electrostatic fields of molecules is a crucial step in the construction of accurate interatomic potentials. We propose a new approach to calculate these charges, which as previous approaches are, is based on the idea of charge equilibration. However, we only allow charge to flow between covalently bonded neighbors by using the concept of so-called split charges. The semiempirical fit parameters in our approach do not only reflect atomic properties (electronegativity and atomic hardness) but also bond-dependent properties. The new method contains two popular but hitherto disjunct approaches as limiting cases. We apply our methodology to a set of molecules containing the elements silicon, carbon, oxygen, and hydrogen. Effective charges derived from electrostatic potential surfaces can be predicted more than twice as accurately as with previous works, at the expense of one additional fit parameter per bond type controlling the polarizability between two bonded atoms. Additional bond-type parameters can be introduced, but barely improve the results. An increase in accuracy of only 30% over existing techniques is achieved when predicting Mulliken charges. However, this could be improved with additional bond-type parameters.

Prediction of reaction barriers and force-induced instabilities under mechanochemical conditions with an approximate model: A case study of the ring opening of 1,3-cyclohexadiene

Mechanochemistry, the use of mechanical stresses to activate chemical reactions, has emerged as a topic of significant interest. The present study examines the use of an approximate model for the prediction of reaction barriers under mechanochemical conditions using the ring opening of 1,3-cyclohexadiene along conrotatory and disrotatory directions as a specific test case. To do this, reaction barriers are evaluated using quantum chemical methods with an external force applied between various pairs of atoms. The results show that the consequent effects on the barrier exhibit a significant dependence on the locations of the atoms used to apply the external force, and in some cases, force-induced instabilities occur that alter the fundamental nature of the reaction pathway. The ability of an approximate model based on a second-order expansion of the force-modified potential energy with respect to nuclear coordinates to reproduce this behavior is then assessed. Good agreement between the results obtained through the quantum chemical calculations and approximate model is attained when force-induced instabilities do not occur. In addition, a strategy for predicting when such instabilities occur is presented and found to yield results that are in qualitative agreement with the quantum chemical calculations. Finally, the response of the system to the external force is interpreted in terms of the parameters entering the model, which correspond to interatomic distances and stiffnesses, and possibly sheds lights on ways to design molecules that exhibit a desired chemical response to mechanochemical conditions.

Simple direct formation of self-assembled N-heterocyclic carbene monolayers on gold and their application in biosensing

The formation of organic films on gold employing N-heterocyclic carbenes (NHCs) has been previously shown to be a useful strategy for generating stable organic films. However, NHCs or NHC precursors typically require inert atmosphere and harsh conditions for their generation and use. Herein we describe the use of benzimidazolium hydrogen carbonates as bench stable solid precursors for the preparation of NHC films in solution or by vapour-phase deposition from the solid state. The ability to prepare these films by vapour-phase deposition permitted the analysis of the films by a variety of surface science techniques, resulting in the first measurement of NHC desorption energy (158±10 kJ mol−1) and confirmation that the NHC sits upright on the surface. The use of these films in surface plasmon resonance-type biosensing is described, where they provide specific advantages versus traditional thiol-based films.

Spirooxazine to merooxazine interconversion in the presence and absence of zinc: approach to a bistable photochemical switch.

A spironaphthoxazine (SO) photoswitch was synthesized, and its photochromic behaviors were investigated. SO underwent reversible ring-opening/closure isomerization between a spirocyclic isomer (closed form) and a merocyanine (MO isomer, open form) upon ultraviolet light irradiation. For the model SO in this work, the thermal equilibrium is substantially shifted toward the spirocyclic isomer even at -30.0 °C. However, addition of zinc, as Zn(ClO(4))(2), exerted an important effect on the thermal reversion process from the open (MO) to the closed form (SO). Kinetic analysis showed that thermal reversion with zinc is retarded more than 13-fold, significantly improving bistability. Moreover, introduction of zinc to the spirooxazine-merooxazine (SO-MO) system resulted in a new absorption band readily distinguishable from the bands arising from spirooxazine and merooxazine. For the first time, to the best of our knowledge, the microscopic rate constants for: MO photogeneration from SO (k(1)), thermal reversion of MO to SO (k(2)), complexation of MO with zinc (k(3)) and for dissociation of the complex, MO-Zn (k(4)), as well as for the ionization equilibria of Zn(ClO(4))(2) have been evaluated. The preferred transoid structures of MO and those of MO-Zn derived from the preferred MO structures are considered. Although the kinetic study does not permit elucidation of the nature of zinc binding to MO to give MO-Zn, nor the precursor isomers of MO, a DFT calculational study in progress should shed light on the structure and relative stability of these essential intermediates.

Mechanistic and computational study of a palladacycle-catalyzed decomposition of a series of neutral phosphorothioate triesters in methanol.

The methanolytic cleavage of a series of O,O-dimethyl O-aryl phosphorothioates (1a−g) catalyzed by a C,N-palladacycle, (2-[N,N-dimethylamino(methyl)phenyl]-C1,N)(pyridine) palladium(II) triflate (3), at 25 °C and sspH 11.7 in methanol is reported, along with data for the methanolytic cleavage of 1a−g. The methoxide reaction gives a linear log k2−OMe vs sspKa (phenol leaving group) Brønsted plot having a gradient of βlg = −0.47 ± 0.03, suggesting about 34% cleavage of the P−OAr bond in the transition state. On the other hand, the 3-catalyzed cleavage of 1 gives a Brønsted plot with a downward break at sspKa (phenol) 13, signifying a change in the rate-limiting step in the catalyzed reaction, with the two wings having βlg values of 0.0 ± 0.03 and −1.93 ± 0.06. The rate-limiting step for good substrates with low leaving group sspKa values is proposed to be substrate/pyridine exchange on the palladacycle, while for substrates with poor leaving groups, the rate-limiting step is a chemical one with extensive cleavage of the P−OAr bond. DFT calculations support this process and also identify two intermediates, namely, one where substrate/pyridine interchange has occurred to give the palladacycle coordinated to substrate through the S═P linkage and to methoxide (6) and another where intramolecular methoxide attack has occurred on the P═S unit to give a five-coordinate phosphorane (7) doubly coordinated to Pd via the S− and through a bridging methoxide linked to P and Pd. Attempts to identify the existence of the phosphorane by 31P NMR in a d4-methanol solution containing 10 mM each of 3, trimethyl phosphorothioate (a very slow cleaving substrate), and methoxide proved unsuccessful, instead showing that the phosphorothioate was slowly converted to trimethyl phosphate, with the palladacycle decomposing to Pd0 and free pyridine. These results provide the first reported example where a palladacycle-promoted solvolysis reaction exhibits a break in the Brønsted plot signifying at least one intermediate, while the DFT calculations provide further insight into a more complex mechanism involving two intermediates.

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.

Methanolysis of thioamide promoted by a simple palladacycle is accelerated by 10(8) over the methoxide-catalyzed reaction.

Palladacycle 1 catalyzes the methanolytic cleavage of N-methyl-N-(4-nitrophenyl)thiobenzamide (4) via a mechanism involving formation of a Pd-bound tetrahedral intermediate (TI). The rate constant for decomposition of the complex formed between 1, methoxide, and 4 is 9.3 s(-1) at 25 °C; this reaction produces methyl thiobenzoate and N-methyl-4-nitroaniline. The ratio of the second-order rate constant for the catalyzed reaction, given as k(cat)/K(d), relative to that of the methoxide-promoted reaction is 3 × 10(8), representing a very large catalysis of thioamide bond cleavage by a synthetic metal complex.

On the pressure-induced loss of crystallinity in orthophosphates of zinc and calcium

A recently suggested mechanism for the stress memory of various metal phosphates is investigated experimentally. Based on first-principles simulations [N. J. Mosey et al., Science 307, 1612 (2005)], it had been argued that atoms with flexible coordination, such as zinc or heavy-metal cations, act as network-forming agents, undergoing irreversible pressure-induced changes in bonding that lead to increased connectivity between phosphate anions. In the present study, orthophosphates of zinc and calcium were exposed to high pressures on surfaces and in diamond anvil cells. An additional set of first-principles simulations was accomplished on alpha-orthophosphate of zinc, which suggested that this material was already cross-linked before compression but that it nevertheless underwent a reversible coordination change under pressure in agreement with the experimental results presented here. Raman spectra indicate an irreversible, pressure-induced loss of long-range crystallinity. The pressures required to induce these changes are around 7 GPa for the zinc phosphates, while they are close to 21 GPa for the calcium phosphates. Hydrogenation of the metal phosphate lowers the threshold pressure by approximately 2-3 GPa in both cases. Moreover, alpha-orthophosphate of zinc could be partially amorphisized under nonisotropic pressure on copper foils.

N-heterocyclic carbene self-assembled monolayers on copper and gold : dramatic effect of wingtip groups on binding, orientation and assembly

Abstract Self‐assembled monolayers of N‐heterocyclic carbenes (NHCs) on copper are reported. The monolayer structure is highly dependent on the N,N‐substituents on the NHC. On both Cu(111) and Au(111), bulky isopropyl substituents force the NHC to bind perpendicular to the metal surface while methyl‐ or ethyl‐substituted NHCs lie flat. Temperature‐programmed desorption studies show that the NHC binds to Cu(111) with a desorption energy of E des=152±10 kJ mol−1. NHCs that bind upright desorb cleanly, while flat‐lying NHCs decompose leaving adsorbed organic residues. Scanning tunneling microscopy of methylated NHCs reveals arrays of covalently linked dimers which transform into adsorbed (NHC)2Cu species by extraction of a copper atom from the surface after annealing.