Bulk and surface energetics of lithium hydride crystal: benchmarks from quantum Monte Carlo and quantum chemistry

We show how accurate benchmark values of the surface formation energy of crystalline lithium hydride can be computed by the complementary techniques of quantum Monte Carlo (QMC) and wavefunction-based molecular quantum chemistry. To demonstrate the high accuracy of the QMC techniques, we present a detailed study of the energetics of the bulk LiH crystal, using both pseudopotential and all-electron approaches. We show that the equilibrium lattice parameter agrees with experiment to within 0.03 %, which is around the experimental uncertainty, and the cohesive energy agrees to within around 10 meV per formula unit. QMC in periodic slab geometry is used to compute the formation energy of the LiH (001) surface, and we show that the value can be accurately converged with respect to slab thickness and other technical parameters. The quantum chemistry calculations build on the recently developed hierarchical scheme for computing the correlation energy of a crystal to high precision. We show that the hierarchical scheme allows the accurate calculation of the surface formation energy, and we present results that are well converged with respect to basis set and with respect to the level of correlation treatment. The QMC and hierarchical results for the surface formation energy agree to within about 1 %.Comment: 16 pages, 4 figure

Elastic and magnetic effects on the infrared phonon spectra of MnF2

We measured the temperature dependent infrared reflectivity spectra of MnF2 between 4 K and room temperature. We show that the phonon spectrum undergoes a strong renormalization at TN. The ab-initio calculation we performed on this compound accurately predict the magnitude and the direction of the phonon parameters changes across the antiferromagnetic transition, showing that they are mainly induced by the magnetic order. In this material, we found that the dielectric constant is mostly from phonon origin. The large change in the lattice parameters with temperature seen by X-ray diffraction as well as the A2u phonon softening below TN indicate that magnetic order induced distortions in MnF2 are compatible with the ferroelectric instabilities observed in TiO2, FeF2 and other rutile-type fluorides. This study also shows the anomalous temperature evolution of the lower energy Eu mode in the paramagnetic phase, which can be compared to that of the B1g one seen by Raman spectroscopy in many isostructural materials. This was interpreted as being a precursor of a phase transition from rutile to CaCl2 structure which was observed under pressure in ZnF2.Comment: 8 pages, 8 figures, updated version accepted in PR

Electronic structure of GdN, and the influence of exact exchange

GdN bulk is studied with the local density approximation, on the Hartree-Fock level, and on the level of the hybrid functional B3LYP. A local basis set formalism is used, as implemented in the present CRYSTAL06 release. It is demonstrated that the code is technically capable of treating this system with its 4f electrons explicitly, i.e. out of the core. The band structure at the level of the local density approximation is in good agreement with earlier calculations and is found to be half-metallic. The Hartree-Fock band structure is insulating with a large gap. Interestingly, three solutions were found at the B3LYP level. The lowest of them is insulating for majority spin, and the Fermi surface for minority spin consists only of points, resulting in a very low density of states around the Fermi level.Comment: to appear in J. Phys.: Condensed Matte

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

The infrared vibrational spectrum of andradite-grossular solid solutions: A quantum mechanical simulation

Infrared spectroscopy is a powerful technique for the characterization of minerals, permitting insights into their structural and thermodynamic properties. The intrinsic complexity of mineral solid solutions makes the interpretation of their spectroscopic data a challenging task. In this work, the IR vibrational spectra of andradite-grossular (Ca3Fe2Si3O12–Ca3Al2Si3O12) solid solutions were simulated at the ab initio level with the CRYSTAL09 code by using a large all-electron Gaussian-type basis set and the B3LYP hybrid functional. All the 23 symmetry-independent configurations resulting from the substitution of 1 to 8 Fe atoms with Al atoms in the 16a octahedral site of the andradite primitive cell were considered. The IR active transverse optical frequencies and their intensities were computed. Graphical representation of the spectra, animation of the modes and isotopic substitution of the cations were used as additional interpretation tools. The dominant high-frequency modes, corresponding to Si-O stretching motions, show a simple linear behavior of both frequencies and intensities with respect to the binary composition; this trend is related to the linear behavior of the mean lattice parameter. Also the frequencies of the low-energy bands show, roughly speaking, a linear dependence on composition; however, the behavior of the dominant intensities is more complicated and strongly connected to the Al and Fe atomic fraction. When considering different possible structures at fixed composition, some spectral features display a dependence upon short-range Y cation ordering. Overall, we show how ab initio calculations permit to analyze complex systems such as solid solutions, establishing relations among structure and properties and providing critical and robust interpretations to the experimental findings

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

Vibrational analysis of paraelectric–ferroelectric transition of linbo3: An ab-initio quantum mechanical treatment

The phase transitions between paraelectric (PE) and ferroelectric (FE) isomorph phases of LiNbO3 have been investigated quantum mechanically by using a Gaussian-type basis set, the B3LYP hybrid functional and the CRYSTAL17 code. The structural, electronic and vibrational properties of the two phases are analyzed. The vibrational frequencies evaluated at the Γ point indicate that the paraelectric phase is unstable, with a complex saddle point with four negative eigenvalues. The energy scan of the A2u mode at −215 cm−1 (i215) shows a dumbbell potential with two symmetric minima. The isotopic substitution, performed on the Li and Nb atoms, allows interpretation of the nontrivial mechanism of the phase transition. The ferroelectric phase is more stable than the paraelectric one by 0.32 eV