The capabilities of the Crystal14program are presented, and the improvements made with respect to the previous Crystal09version discussed. Crystal14is anab initiocode that uses a Gaussian‐type basis set: both pseudopotential and all‐electron strategies are permitted; the latter is not much more expensive than the former up to the first‐second transition metal rows of the periodic table. A variety of density functionals is available, including as an extreme case Hartree–Fock; hybrids of various nature (global, range‐separated, double) can be used. In particular, a very efficient implementation of global hybrids, such as popular B3LYP and PBE0 prescriptions, allows for such calculations to be performed at relatively low computational cost. The program can treat on the same grounds zero‐dimensional (molecules), one‐dimensional (polymers), two‐dimensional (slabs), as well as three‐dimensional (3D; crystals) systems. No spurious 3D periodicity is required for low‐dimensional systems as happens when plane‐waves are used as a basis set. Symmetry is fully exploited at all steps of the calculation; this permits, for example, to investigate nanotubes of increasing radius at a nearly constant cost (better than linear scaling!) or to perform self‐consistent‐field (SCF) calculations on fullerenes as large as (10,10), with 6000 atoms, 84,000 atomic orbitals, and 20 SCF cycles, on a single core in one day. Three versions of the code exist, serial, parallel, and massive‐parallel. In the second one, the most relevant matrices are duplicated, whereas in the third one the matrices in reciprocal space are distributed for diagonalization. All the relevant vectors are now dynamically allocated and deallocated after use, making Crystal14much more agile than the previous version, in which they were statically allocated. The program now fits more easily in low‐memory machines (as many supercomputers nowadays are). Crystal14can be used on parallel machines up to a high number of cores (benchmarks up to 10,240 cores are documented) with good scalability, the main limitation remaining the diagonalization step. Many tensorial properties can be evaluated in a fully automated way by using a single input keyword: elastic, piezoelectric, photoelastic, dielectric, as well as first and second hyperpolarizabilies, electric field gradients, Born tensors and so forth. Many tools permit a complete analysis of the vibrational properties of crystalline compounds. The infrared and Raman intensities are now computed analytically and related spectra can be generated. Isotopic shifts are easily evaluated, frequencies of only a fragment of a large system computed and nuclear contribution to the dielectric tensor determined. New algorithms have been devised for the investigation of solid solutions and disordered systems. The topological analysis of the electron charge density, according to the Quantum Theory of Atoms in Molecules, is now incorporated in the code via the integrated merge of theTopondpackage. Electron correlation can be evaluated at the Möller–Plesset second‐order level (namely MP2) and a set of double‐hybrids are presently available via the integrated merge with the Cryscorprogram.

The Crystal program for quantum-mechanical simulations of materials has been bridging the realm of molecular quantum chemistry to the realm of solid state physics for many years, since its first public version released back in 1988. This peculiarity stems from the use of atom-centered basis functions within a linear combination of atomic orbitals (LCAO) approach and from the corresponding efficiency in the evaluation of the exact Fock exchange series. In particular, this has led to the implementation of a rich variety of hybrid density functional approximations since 1998. Nowadays, it is acknowledged by a broad community of solid state chemists and physicists that the inclusion of a fraction of Fock exchange in the exchange-correlation potential of the density functional theory is key to a better description of many properties of materials (electronic, magnetic, mechanical, spintronic, lattice-dynamical, etc.). Here, the main developments made to the program in the last five years (i.e., since the previous release, Crystal17) are presented and some of their most noteworthy applications reviewed.

CRYSCOR is a periodic post-Hartree–Fock program based on local functions in direct space, i.e., Wannier functions and projected atomic orbitals. It uses atom centered Gaussians as basis functions. The Hartree–Fock reference, as well as symmetry information, is provided by the CRYSTAL program. CRYSCOR presently features an efficient and parallel implementation of periodic local second order Møller–Plesset perturbation theory (MP2), which allows us to study 1D-, 2D- and 3D-periodic systems beyond 1000 basis functions per unit cell. Apart from the correlation energy also the MP2 density matrix, and from that the Compton profile, are available. Very recently, a new module for calculating excitonic band gaps at the uncorrelated Configuration-Interaction-Singles (CIS) level has been added. Other advancements include new extrapolation techniques for calculating surface adsorption on semi-infinite solids. In this paper the diverse features and recent advances of the present CRYSCOR version are illustrated by exemplary applications to various systems: the adsorption of an argon monolayer on the MgO (100) surface, the rolling energy of a boron nitride nanoscroll, the relative stability of different aluminosilicates, the inclusion energy of methane in methane–ice-clathrates, and the effect of electron correlation on charge and momentum density of α-quartz. Furthermore, we present some first tentative CIS results for excitonic band gaps of simple 3D-crystals, and their dependence on the diffuseness of the basis set.