Estudio computacional cuanto-mecánico de la interacción del hidrógeno molecular [H 2 ] con el ZnO 4 (1,4-bencenodicarboxilato) 3 [MOF-5]

El potencial del ZnO4(1,4-bencenodicarboxilato)3 [MOF-5] para almacenar hidrógeno vía adsorción es evaluado a través de cálculos cuanto-mecánicos de la energía de interacción entre moléculas de H2 y diferentes sitios de la estructura cristalina de este material. En una primera fase del estudio, la energía de interacción es calculada en modelos periódi­cos usando el funcional B3LYP como nivel de teoría junto con un set de funciones base localizadas de tipo Gaussiano que es usado para expandir la función de onda del siste­ma. Posteriormente, con el objetivo de obtener resultados que describan mejor el proceso, se incluye a los valores calculados la contribución de las fuerzas de dispersión mediante el método P-ONIOMusando cálculos MP2 en modelos moleculares extraídos de la estructura cristalina

On the use of symmetry in configurational analysis for the simulation of disordered solids

The starting point for a quantum mechanical investigation of disordered systems usually implies calculations on a limited subset of configurations, generated by defining either the composition of interest or a set of compositions ranging from one end member to another, within an appropriate supercell of the primitive cell of the pure compound. The way in which symmetry can be used in the identification of symmetry independent configurations (SICs) is discussed here. First, Pólya's enumeration theory is adopted to determine the number of SICs, in the case of both varying and fixed composition, for colors numbering two or higher. Then, De Bruijn's generalization is presented, which allows analysis of the case where the colors are symmetry related, e.g. spin up and down in magnetic systems. In spite of their efficiency in counting SICs, neither Pólya's nor De Bruijn's theory helps in solving the difficult problem of identifying the complete list of SICs. Representative SICs are obtained by adopting an orderly generation approach, based on lexicographic ordering, which offers the advantage of avoiding the (computationally expensive) analysis and storage of all the possible configurations. When the number of colors increases, this strategy can be combined with the surjective resolution principle, which permits the efficient generation of SICs of a problem in |R| colors starting from the ones obtained for the (|R| − 1)-colors case. The whole scheme is documented by means of three examples: the abstract case of a square with C4v symmetry and the real cases of the garnet and olivine mineral families

Symmetry and random sampling of symmetry independent configurations for the simulation of disordered solids

A symmetry-adapted algorithm producing uniformly at random the set of symmetry independent configurations (SICs) in disordered crystalline systems or solid solutions is presented here. Starting from Pólya's formula, the role of the conjugacy classes of the symmetry group in uniform random sampling is shown. SICs can be obtained for all the possible compositions or for a chosen one, and symmetry constraints can be applied. The approach yields the multiplicity of the SICs and allows us to operate configurational statistics in the reduced space of the SICs. The present low-memory demanding implementation is briefly sketched. The probability of finding a given SIC or a subset of SICs is discussed as a function of the number of draws and their precise estimate is given. The method is illustrated by application to a binary series of carbonates and to the binary spinel solid solution Mg(Al,Fe)2O4

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