M. Waroquier

We demonstrate how a limited number of ab initio calculations in combination with a simple Debye model can predict a concentration- and temperature-dependent mixing enthalpy for a binary system. Fe-Mo is taken as a test case, and our predictions are compared with phase diagram information and a recently measured heat of solution for Mo in Fe. Crystallographic and magnetic information is calculated for the λ and μ intermetallic phases in the Fe-Mo phase diagram as well. The present methodology can be useful for making a quick survey of mixing enthalpies in a large set of binary systems, in particular in the dilute concentration ranges where tabulated data are often lacking and where calphad-style modeling is less reliable.

The concept of the atomic charge is extensively used to model the electrostatic properties of proteins. Atomic charges are not only the basis for the electrostatic energy term in biomolecular force fields but are also derived from quantum mechanical computations on protein fragments to get more insight into their electronic structure. Unfortunately there are many atomic charge schemes which lead to significantly different results, and it is not trivial to determine which scheme is most suitable for biomolecular studies. Therefore, we present an extensive methodological benchmark using a selection of atomic charge schemes [Mulliken, natural, restrained electrostatic potential, Hirshfeld-I, electronegativity equalization method (EEM), and split-charge equilibration (SQE)] applied to two sets of penta-alanine conformers. Our analysis clearly shows that Hirshfeld-I charges offer the best compromise between transferability (robustness with respect to conformational changes) and the ability to reproduce electrostatic properties of the penta-alanine. The benchmark also considers two charge equilibration models (EEM and SQE), which both clearly fail to describe the locally charged moieties in the zwitterionic form of penta-alanine. This issue is analyzed in detail because charge equilibration models are computationally much more attractive than the Hirshfeld-I scheme. Based on the latter analysis, a straightforward extension of the SQE model is proposed, SQE+Q⁰, that is suitable to describe biological systems bearing many locally charged functional groups.

The difference in reactivity between the activated 2-bromomethyl-1-tosylaziridine and the non-activated 1-benzyl-2-(bromomethyl)aziridine with respect to sodium methoxide was analyzed by means of DFT calculations within the supermolecule approach, taking into account explicit solvent molecules. In addition, the reactivity of epibromohydrin with regard to sodium methoxide was assessed as well. The barriers for direct displacement of bromide by methoxide in methanol are comparable for all three heterocyclic species under study. However, ring opening was found to be only feasible for the epoxide and the activated aziridine, and not for the non-activated aziridine. According to these computational analyses, the synthesis of chiral 2-substituted 1-tosylaziridines can take place with inversion (through ring opening/ring closure) or retention (through direct bromide displacement) of configuration upon treatment of the corresponding 2-(bromomethyl)aziridines with one equivalent of a nucleophile, whereas chiral 2-substituted 1-benzylaziridines are selectively obtained with retention of configuration (via direct bromide displacement). Furthermore, the computational results showed that explicit accounting for solvent molecules is required to describe the free energy profile correctly. To verify the computational findings experimentally, chiral 1-benzyl-2-(bromomethyl)aziridines and 2-bromomethyl-1-tosylaziridines were treated with sodium methoxide in methanol. The presented work concerning the reactivity of 2-bromomethyl-1-tosylaziridine stands in contrast to the behaviour of the corresponding 1-tosyl-2-(tosyloxymethyl)aziridine with respect to nucleophiles, which undergoes a clean ring-opening/ring-closure process with inversion of configuration at the asymmetric aziridine carbon atom.