Alessandro Erba

Elastic properties of six silicate garnet end members from accurate ab initio simulations

AbstractThe elastic properties of six silicate garnet end members, among the most important rock-forming minerals, are investigated here for the first time via accurate ab initio theoretical simulations. The Crystal program is used, which works within periodic boundary conditions and allows for all-electron basis sets to be adopted. From the computed elastic tensor, Christoffel’s equation is solved along a set of crystallographic directions in order to fully characterize the seismic wave velocity anisotropy in such materials. Polycrystalline isotropic aggregate elastic properties are derived from the computed single-crystal data via the Voigt-Reuss-Hill averaging procedure. Transferability of the elastic properties from end members to their solid solutions with different chemical compositions is also addressed.

The vibrational spectrum of CaCO3 aragonite: A combined experimental and quantum-mechanical investigation

The vibrational properties of CaCO(3) aragonite have been investigated both theoretically, by using a quantum mechanical approach (all electron Gaussian type basis set and B3LYP HF-DFT hybrid functional, as implemented in the CRYSTAL code) and experimentally, by collecting polarized infrared (IR) reflectance and Raman spectra. The combined use of theory and experiment permits on the one hand to analyze the many subtle features of the measured spectra, on the other hand to evidentiate limits and deficiencies of both approaches. The full set of TO and LO IR active modes, their intensities, the dielectric tensor (in its static and high frequency components), and the optical indices have been determined, as well as the Raman frequencies. Tools such as isotopic substitution and graphical animation of the modes are available, that complement the analysis of the spectrum.

Quantum mechanical predictions to elucidate the anisotropic elastic properties of zeolitic imidazolate frameworks: ZIF-4 vs. ZIF-zni

We use ab initio density functional theory (DFT) to elucidate the mechanical properties of two topologically distinct zeolitic imidazolate framework (ZIF) materials: ZIF-4 and ZIF-zni, both of which have the same chemical composition of Zn(Im)2 [Im = C3H3N2−] and are constructed from an identical Zn—Im—Zn basic building block. The CRYSTAL code was used to compute the single-crystal elastic constants Cij of the (orthorhombic) ZIF-4 and (tetragonal) ZIF-zni structures at the PBE level of theory. Through tensorial analysis of the Cij, we reveal the three-dimensional representation surfaces of the Youngs modulus, shear modulus, Poissons ratio and linear compressibility, which enable us to describe the detailed elasticity behaviour and to pinpoint basic crystal structure–property correlations. Notably, we discover that ZIF-4 can potentially exhibit a negative Poissons ratio, thereby representing the first example of an ‘auxetic-ZIF’ to be identified to date. Furthermore, we show that our DFT predictions are consistent with recently reported experimental measurements of the Youngs and bulk moduli of such complex ZIF structures.

Anisotropic displacement parameters for molecular crystals from periodic Hartree–Fock and density functional theory calculations

Fully periodic Hartree–Fock and density functional theory calculations have been used to compute the anisotropic displacement parameters (ADPs) of molecular crystals at different temperatures by using the CRYSTAL code. Crystalline urea was adopted as a benchmark system to investigate the dependence on basis set and Hamiltonian. The results were compared with ADPs derived from neutron diffraction experiments. The approach can estimate the internal ADPs, corresponding to the contributions of high-frequency intramolecular vibrations, and for these internal contributions the results are almost independent of the basis set and Hamiltonian. Much larger variations and discrepancies from neutron diffraction experiments are seen for the external, low-frequency modes, which become dominant at higher temperatures. The approach was then tested on benzene and urotropine. Finally, ADPs of l-alanine were predicted at the B3LYP/6-31G(d,p) level of theory. The total ADPs, including low-frequency external modes, are underestimated, but can be brought into good agreement with the experimental ADPs by introducing a Gruneisen parameter, which partly accounts for anharmonicity of the potential energy surface, but likely also contains contributions from other deficiencies of the calculations.

Accurate dynamical structure factors from ab initio lattice dynamics: the case of crystalline silicon.

A fully ab initio technique is discussed for the determination of dynamical X‐ray structure factors (XSFs) of crystalline materials, which is based on a standard Debye–Waller (DW) harmonic lattice dynamical approach with all‐electron atom‐centered basis sets, periodic boundary conditions, and one‐electron Hamiltonians. This technique requires an accurate description of the lattice dynamics and the electron charge distribution of the system. The main theoretical parameters involved and final accuracy of the technique are discussed with respect to the experimental determinations of the XSFs at 298 K of crystalline silicon. An overall agreement factor of 0.47% between the ab initio predicted values and the experimental determinations is found. The best theoretical determination of the anisotropic displacement parameter, of silicon is here 60.55 × 10−4 Å2, corresponding to a DW factor B = 0.4781 Å2.

High pressure elastic properties of minerals from ab initio simulations: The case of pyrope, grossular and andradite silicate garnets

A computational strategy is devised for the accurate ab initio simulation of elastic properties of crystalline materials under pressure. The proposed scheme, based on the evaluation of the analytical stress tensor and on the automated computation of pressure-dependent elastic stiffness constants, is implemented in the CRYSTAL solid state quantum-chemical program. Elastic constants and related properties (bulk, shear and Young moduli, directional seismic wave velocities, elastic anisotropy index, Poissons ratio, etc.) can be computed for crystals of any space group of symmetry. We apply such a technique to the study of high-pressure elastic properties of three silicate garnet end-members (namely, pyrope, grossular, and andradite) which are of great geophysical interest, being among the most important rock-forming minerals. The reliability of this theoretical approach is proved by comparing with available experimental measurements. The description of high-pressure properties provided by several equations of state is also critically discussed.

The vibration properties of the (n,0) boron nitride nanotubes from ab initio quantum chemical simulations

The vibration spectrum of single-walled zigzag boron nitride (BN) nanotubes is simulated with an ab initio periodic quantum chemical method. The trend towards the hexagonal monolayer (h-BN) in the limit of large tube radius R is explored for a variety of properties related to the vibrational spectrum: vibration frequencies, infrared intensities, oscillator strengths, and vibration contributions to the polarizability tensor. The (n,0) family is investigated in the range from n = 6 (24 atoms in the unit cell and tube radius R = 2.5 Å) to n = 60 (240 atoms in the cell and R = 24.0 Å). Simulations are performed using the CRYSTAL program which fully exploits the rich symmetry of this class of one-dimensional periodic systems: 4n symmetry operators for the general (n,0) tube. Three sets of infrared active phonon bands are found in the spectrum. The first one lies in the 0-600 cm(-1) range and goes regularly to zero when R increases; the connection between these normal modes and the elastic and piezoelectric constants of h-BN is discussed. The second (600-800 cm(-1)) and third (1300-1600 cm(-1)) sets tend regularly, but with quite different speed, to the optical modes of the h-BN layer. The vibrational contribution of these modes to the two components (parallel and perpendicular) of the polarizability tensor is also discussed.

Periodic density functional theory and local-MP2 study of the librational modes of Ice XI.

Two periodic codes, CRYSTAL and CRYSCOR, are here used to simulate and characterize the librational modes of the nu(R) band of Ice XI: this band has been found experimentally to be the region of the vibrational spectrum of ordinary ice most affected by the transition from the proton-disordered (Ice Ih) to the proton-ordered (Ice XI) phase. With CRYSTAL, the problem is solved using Hartree-Fock (HF), pure Kohn-Sham (PW91) or hybrid (B3LYP) one-electron Hamiltonians: the harmonic approximation is employed to obtain the vibrational spectrum after optimizing the geometry. The B3LYP results are those in best agreement with the experiment. For a given crystalline geometry, CRYSCOR computes the energy per cell in an ab initio HF+MP2 approximation using a local-correlation approach; this technique is employed for recalculating the frequencies of the different modes identified by the B3LYP approach, by fully accounting for long range dispersive interactions. The effect of anharmonicity is evaluated separately for each mode both in the B3LYP and HF+MP2 case. The two approaches accurately reproduce the four-peak structure of the librational band. The harmonic B3LYP nu(R) bandwidth of 70 meV is lowered to 60 meV by anharmonic corrections, and becomes 57 meV in the HF+MP2 anharmonic calculation, in excellent agreement with the experimental IINS data (56-59 meV). The assignment of the librational modes is discussed.

DFT and Local-MP2 Periodic Study of the Structure and Stability of Two Proton-Ordered Polymorphs of Ice

The equilibrium geometry and the formation energy of two periodic polymorphs of Ice have been theoretically studied: the former (Ice XI, crystal group Cmc2(1)) is experimentally observed as the most stable structure at low temperature and pressure; the latter (crystal group Pna2(1)) is the simplest proton-ordered model of ordinary ice. With the Crystal code, the problem is solved using Hartree-Fock (HF), pure Kohn-Sham (PW91), or hybrid (B3LYP) one-electron Hamiltonians. The B3LYP results are those in best agreement with the experiment. Using the B3LYP-optimized geometry and starting from the corresponding HF Crystal solution, the energetics of the two polymorphs have been investigated at an ab initio MP2 level using the Cryscor code, based on a local-correlation approach: these calculations have allowed us not only to confirm the excellent B3LYP results as concerns the formation energy and the relative stability of the two structures but also to analyze the role in this respect of the intra- and intermolecular contributions to the correlation energy. Since both Crystal and Cryscor adopt a basis set of localized Gaussian-type functions and since very small energy differences are involved, utter attention has been paid to correcting for the basis set superposition error in the calculation of formation energies.

Quantum‐mechanical condensed matter simulations with CRYSTAL

The latest release of the Crystal program for solid‐state quantum‐mechanical ab initio simulations is presented. The program adopts atom‐centered Gaussian‐type functions as a basis set, which makes it possible to perform all‐electron as well as pseudopotential calculations. Systems of any periodicity can be treated at the same level of accuracy (from 0D molecules, clusters and nanocrystals, to 1D polymers, helices, nanorods, and nanotubes, to 2D monolayers and slab models for surfaces, to actual 3D bulk crystals), without any artificial repetition along nonperiodic directions for 0–2D systems. Density functional theory calculations can be performed with a variety of functionals belonging to several classes: local‐density (LDA), generalized‐gradient (GGA), meta‐GGA, global hybrid, range‐separated hybrid, and self‐consistent system‐specific hybrid. In particular, hybrid functionals can be used at a modest computational cost, comparable to that of pure LDA and GGA formulations, because of the efficient implementation of exact nonlocal Fock exchange. Both translational and point‐symmetry features are fully exploited at all steps of the calculation, thus drastically reducing the corresponding computational cost. The various properties computed encompass electronic structure (including magnetic spin‐polarized open‐shell systems, electron density analysis), geometry (including full or constrained optimization, transition‐state search), vibrational properties (frequencies, infrared and Raman intensities, phonon density of states), thermal properties (quasi‐harmonic approximation), linear and nonlinear optical properties (static and dynamic [hyper]polarizabilities), strain properties (elasticity, piezoelectricity, photoelasticity), electron transport properties (Boltzmann, transport across nanojunctions), as well as X‐ray and inelastic neutron spectra. The program is distributed in serial, parallel, and massively parallel versions. In this paper, the original developments that have been devised and implemented in the last 4 years (since the distribution of the previous public version, Crystal14, occurred in December 2013) are described.