Exploitation of symmetry in periodic Self-Consistent-Field ab initio calculations: application to large three-dimensional compounds

Symmetry can dramatically reduce the computational cost (running time and memory allocation) of Self-Consistent-Field ab initio calculations for crystalline systems. Crucial for running time is use of symmetry in the evaluation of one- and two-electron integrals, diagonalization of the Fock matrix at selected points in reciprocal space, reconstruction of the density matrix. As regards memory allocation, full square matrices (overlap, Fock and density) in the Atomic Orbital (AO) basis are avoided and a direct transformation from the packed AO to the SACO (Symmetry Adapted Crystalline Orbital) basis is performed, so that the largest matrix to be handled has the size of the largest sub-block in the latter basis. We here illustrate the effectiveness of this scheme, following recent advancements in the CRYSTAL code, concerning memory allocation and direct basis set transformation. Quantitative examples are given for large unit cell systems, such as zeolites (all-silica faujasite and silicalite MFI) and garnets (pyrope). It is shown that the full SCF of 3D systems containing up to 576 atoms and 11136 Atomic Orbitals in the cell can be run with a hybrid functional on a single core PC with 500 MB RAM in about 8 h. © 2014 Science China Press and Springer-Verlag Berlin Heidelberg

Exceptionally low shear modulus in a prototypical imidazole-based metal-organic framework.

Using Brillouin scattering, we measured the single-crystal elastic constants (C(ij)'s) of a prototypical metal-organic framework (MOF): zeolitic imidazolate framework (ZIF)-8 [Zn(2-methylimidazolate)(2)], which adopts a zeolitic sodalite topology and exhibits large porosity. Its C(ij)'s under ambient conditions are (in GPa) C(11)=9.522(7), C(12)=6.865(14), and C(44)=0.967(4). Tensorial analysis of the C(ij)'s reveals the complete picture of the anisotropic elasticity in cubic ZIF-8. We show that ZIF-8 has a remarkably low shear modulus G(min) < or approximately 1 GPa, which is the lowest yet reported for a single-crystalline extended solid. Using ab initio calculations, we demonstrate that ZIF-8's C(ij)'s can be reliably predicted, and its elastic deformation mechanism is linked to the pliant ZnN(4) tetrahedra. Our results shed new light on the role of elastic constants in establishing the structural stability of MOF materials and thus their suitability for practical applications

Periodic ab initio study of the electronic structure of alpha-Al2O3 and AIN(w) surfaces based on localized Wannier functions

The behavior of surface undercoordinated Al3+ cations as electron acceptors is investigated by means of periodic ab initio calculations. The (0001) surfaces of alpha-Al2O3 and AlN(w) were chosen as representative of the different Al coordination types. The electronic structure was analyzed through the valence localized Wannier functions that allow a straightforward description of the electron pairs in terms of chemical concepts. The relevance of the electron back-donation toward surface Al3+ cations in the reconstruction of the surfaces is discussed, and evidence was found that some pi-bond character is displayed by the surface Al-X bonds after relaxation. This explains the tendency to form surface Al-X planar layers in this kind of system

Synthesis of Aluminosilicate Natrolites and Control of Their Tetrahedral Atom Ordering

X1187sciescopu

In situ disorder-order transformation in synthetic gallosilicate zeolites with the NAT topology

Here, we report that synthetic gallosilicate molecular sieves with the NAT topology and Si/Ga ratios close to but slightly higher than 1.50 undergo an in situ transformation under their crystallization conditions. The materials have been studied ex situ by using powder X-ray diffraction, elemental and thermal analyses, and multinuclear MAS NMR. The transformation is characterized by a change in the distribution of Si and Ga of the NAT framework, from a quite (but not completely) disordered phase to a very highly (but not completely) ordered one, accompanied by a change from tetragonal to orthorhombic symmetry. During most of the solution-mediated transformation, no noticeable signs of fresh precipitation, phase segregation, or changes in the chemical composition are detected. Intermediate materials show variations in the degree of Si-Ga ordering and orthorhombic distortion and are not physical mixtures of the disordered and ordered phases. Ab initio calculations strongly suggest a preferential siting of Si in the tetrahedral sites involved in a smaller number of 4-rings in the NAT topology (i.e., the low multiplicity site). The cost of violations of Loewenstein's rule has also been calculated. For this topology and chemical composition the preferential siting and Loewenstein's rule drive together the system to the ordered configuration. A Monte Carlo sampling procedure affords a reasonable model for the initial, mainly disordered state, which fits well within the experimental disorder-order series

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(-)center dot(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