Revisiting the birth of NaCl crystals using molecular dynamics simulation.

Molecular dynamics simulations were used to investigate the initial stage of phase separation mechanisms for an oversaturated electrolytic solution. We developed a low computational cost methodology to determine the simulation frames where the first ionic clusters are formed. By discretizing the simulation box, we obtain a density profile in the moments preceding and succeeding the nuclei's formation. The growth of the clusters identified with our methodology was analyzed until the end of the simulation. Calculation of the Steinhardt parameter showed symmetry of the solid, giving indications that the classical nucleation theory explains the mechanism of the solid formation. The methodology developed was useful for identifying phase separation mechanisms in the nucleation process. At lower concentrations, there was no formation of stable clusters. At intermediate concentrations, the analyses indicate a transition of phases in one stage, from a oversaturate electrolytic solution to a crystalline solid. At high concentration, a transition of phases in two stages, initially, is the formation of a dense liquid, and only after that, crystalline solid formed inside the dense liquid. The change in phase separation mechanism due to increasing oversaturation underscores the importance of precise determination of the driving force for phase separation and concentration limits for each mechanism.

The isobutylene-isobutane alkylation process in liquid HF revisited.

Details on the mechanism of HF catalyzed isobutylene-isobutane alkylation were investigated. On the basis of available experimental data and high-level quantum chemical calculations, a detailed reaction mechanism is proposed taking into account solvation effects of the medium. On the basis of our computational results, we explain why the density of the liquid media and stirring rates are the most important parameters to achieve maximum yield of alkylate, in agreement with experimental findings. The ab initio Car-Parrinello molecular dynamics calculations show that isobutylene is irreversibly protonated in the liquid HF medium at higher densities, leading to the ion pair formation, which is shown to be a minimum on the potential energy surface after optimization using periodic boundary conditions. The HF medium solvates preferentially the fluoride anion, which is found as solvated [FHF](-) or solvated F(-.)(HF)(3). On the other hand, the tert-butyl cation is weakly solvated, where the closest HF molecules appear at a distance of about 2.9 Angstrom with the fluorine termination of an HF chain.