Solvation states of HCl in mixed ether:acid crystals: a computational study

Acid solvation states are investigated in the recently discovered mixed ether:acid crystalline solids. The solids are simulated using on-the-fly molecular dynamics as implemented in the density functional code QUICKSTEP employing Gaussian basis sets. The solids are shown to display a remarkably broad range of acid solvation states, depending on the ether:acid ratio, including proton sharing in the 1:1 case, proton transfer to the ether in 1:2, and perturbed molecular acid in 1:6. The observed variation of the infrared spectra with the composition is accounted for qualitatively with the help of the calculations.     

A new molecular-dynamics based approach for molecular crystal structure search

A new molecular-dynamics based approach is proposed to search for candidate crystal structures of molecular solids. The procedure is based on the observation of spontaneous transitions between ordered and disordered states in molecular-dynamics simulations of an artificial periodic system with a small unit cell. In such a way only the most stable structures are automatically selected. The method can be applied to the solution of crystal structures from low-quality or very complex diffraction data. Tests are presented for H2O-ice polymorphs.     

Analytic solution of the mean spherical approximation for a multi-component plasma

The mean spherical approximation for a mixture of charged hard spheres in a uniform neutralizing background is solved analytically. The factor correlation functions and the excess thermodynamic properties are explicitly expressed through a single parameter, which can be obtained by solving an algebraic equation.     

A recipe for the computation of the free energy barrier and the lowest free energy path of concerted reactions

The recently introduced hills method (Proc. Natl. Acad. Sci. U.S.A. 2002, 99, 12562) is a powerful tool to compute the multidimensional free energy surface of intrinsically concerted reactions. We have extended this method by focusing our attention on localizing the lowest free energy path that connects the stable reactant and product states. This path represents the most probable reaction mechanism, similar to the zero temperature intrinsic reaction coordinate, but also includes finite temperature effects. The transformation of the multidimensional problem to a one-dimensional reaction coordinate allows for accurate convergence of the free energy profile along the lowest free energy path using standard free energy methods. Here we apply the hills method, our lowest free energy path search algorithm, and umbrella sampling to the prototype S(N)2 reaction. The hills method replaces the in many cases difficult problem of finding a good reaction coordinate with choosing relatively simple collective variables, such as the bond lengths of the broken and formed chemical bonds. The second part of the paper presents a guide to using the hills method, in which we test and fine-tune the method for optimal accuracy and efficiency using the umbrella sampling results as a reference.     

Investigating the polymorphism in PR179: a combined crystal structure prediction and metadynamics study

The polymorphism of an industrial important pigment (PR179) was studied with a combination of standard crystal structure prediction and metadynamics. The former provided a starting set of candidate polymorphs whose structural and thermal stability were then probed by metadynamics. Moreover, metadynamics allowed for exploring the free energy surface to look for other possible polymorphs that were not included in the original set. Our calculations indicate that two structures have a high structural stability and are therefore good candidates to be found in experiments. The lower energy phase of the two indeed corresponds to the known polymorph, and we suggest that the other might be a metastable polymorph not yet experimentally discovered.     

The thermal stability of lattice-energy minima of 5-fluorouracil: metadynamics as an aid to polymorph prediction

This paper reports a novel methodology for the free-energy minimization of crystal structures exhibiting strong, anisotropic interactions due to hydrogen bonding. The geometry of the thermally expanded cell was calculated by exploiting the dependence of the free-energy derivatives with respect to cell lengths and angles on the average pressure tensor computed in short molecular dynamics simulations. All dynamic simulations were performed with an elaborate anisotropic potential based on a distributed multipole analysis of the isolated molecule charge density. Changes in structure were monitored via simulated X-ray diffraction patterns. The methodology was used to minimize the free energy at ambient conditions of a set of experimental and hypothetical 5-fluorouracil crystal structures, generated in a search for lattice-energy minima with the same model potential. Our results demonstrate that the majority ( approximately 75%) of lattice-energy minima are thermally stable at ambient conditions, and hence, the free-energy (like the lattice-energy) surface is complex and highly undulating. Metadynamics trajectories (Laio, A.; Parrinello, M. Proc. Natl. Acad. Sci. U.S.A. 2002, 99, 12562) started from the free-energy minima only produced transitions that preserved the hydrogen-bonding motif, and thus, further developments are needed for this method to efficiently explore such free-energy surfaces. The existence of so many free-energy minima, with large barriers for the alteration of the hydrogen-bonding motif, is consistent with the range of motifs observed in crystal structures of 5-fluorouracil and other 5-substituted uracils.     

DFT research on the dehydroxylation reaction of pyrophyllite 1. First-principle molecular dynamics simulations

The dehydroxylation of pyrophyllite involves the reaction of OH groups and elimination of water molecules through two possible mechanisms, one involving the bridging hydroxyl groups of an octahedral Al (3+) pair and the other two hydroxyl groups reacting across the dioctahedral vacancy. First-principles molecular dynamics simulations at the density functional theory level are used together with the metadynamics algorithm to explore the free-energy surface (FES) of the initial step of the dehydroxylation. We observe that the two possible dehydroxylation mechanisms yield similar activation energies at 0 K, but at high temperatures, the cross mechanism has lower free energy than that of the on-site one. The dehydroxylation process produces different semidehydroxylated intermediates that should be taken into account. The role of the temperature in favoring a dehydroxylation nonconcerted chain mechanism over another is here elucidated, and a novel competitive mechanism, which is assisted by the structural apical oxygens in the high-temperature regime, is proposed.