Anharmonic Vibrational States of Solids from DFT Calculations. Part II: Implementation of the VSCF and VCI Methods.

Two methods are implemented in the Crystal program for the calculation of anharmonic vibrational states of solids: the vibrational self-consistent field (VSCF) and the vibrational configuration-interaction (VCI). While the former is a mean-field approach, where each vibrational mode interacts with the average potential of the others, the latter allows for an explicit and complete account of mode-mode correlation. Both schemes are based on the representation of the adiabatic potential energy surface (PES) discussed in Part I, where the PES is expanded in a Taylor's series so as to include up to cubic and quartic terms.

Anharmonic Vibrational States of Solids from DFT Calculations. Part I: Description of the Potential Energy Surface.

A computational approach is presented to compute anharmonic vibrational states of solids from quantum-mechanical DFT calculations by taking into explicit account phonon-phonon couplings via the vibrational configuration interaction (VCI) method. The Born-Oppenheimer potential energy surface (PES) is expanded in a Taylor's series in terms of harmonic normal coordinates, centered at the equilibrium nuclear configuration, is truncated to quartic order, and contains one-mode, two-mode, and three-mode interatomic force constants. The description of the anharmonic terms of the PES involves the numerical evaluation of high-order energy derivatives (cubic and quartic in our case) with respect to nuclear displacements and constitutes the most computationally demanding step in the characterization of anharmonic vibrational states of materials.

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 electronic structure of MgO nanotubes. An ab initio quantum mechanical investigation.

The structural, vibrational and response properties of the (n,0) and (m,m) MgO nanotubes are computed by using a Gaussian type basis set, a hybrid functional (B3LYP) and the CRYSTAL09 code. Tubes in the range 6 ≤ n ≤ 140 and 3 ≤ m ≤ 70 were considered, being n = 2 × m the number of MgO units in the unit cell (so, the maximum number of atoms is 280). Tubes are built by rolling up the fully relaxed 2-D conventional cell (2 MgO units, with oxygen atoms protruding from the Mg plane alternately up and down by 0.

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.

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.

Nuclear motion effects on the density matrix of crystals: an ab initio Monte Carlo harmonic approach.

In the frame of the Born-Oppenheimer approximation, nuclear motions in crystals can be simulated rather accurately using a harmonic model. In turn, the electronic first-order density matrix (DM) can be expressed as the statistically weighted average over all its determinations each resulting from an instantaneous nuclear configuration. This model has been implemented in a computational scheme which adopts an ab initio one-electron (Hartree-Fock or Kohn-Sham) Hamiltonian in the CRYSTAL program.

A new massively parallel version of CRYSTAL for large systems on high performance computing architectures.

Fully ab initio treatment of complex solid systems needs computational software which is able to efficiently take advantage of the growing power of high performance computing (HPC) architectures. Recent improvements in CRYSTAL, a periodic ab initio code that uses a Gaussian basis set, allows treatment of very large unit cells for crystalline systems on HPC architectures with high parallel efficiency in terms of running time and memory requirements. The latter is a crucial point, due to the trend toward architectures relying on a very high number of cores with associated relatively low memory availability.

Static and dynamic coupled perturbed Hartree-Fock vibrational (hyper)polarizabilities of polyacetylene calculated by the finite field nuclear relaxation method.

The vibrational contribution to static and dynamic (hyper)polarizability tensors of polyacetylene are theoretically investigated. Calculations were carried out by the finite field nuclear relaxation (FF-NR) method for periodic systems, newly implemented in the CRYSTAL code, using the coupled perturbed Hartree-Fock scheme for the required electronic properties. The effect of the basis set is also explored, being particularly important for the non-periodic direction perpendicular to the polymer plane.

The first and second static electronic hyperpolarizabilities of zigzag boron nitride nanotubes. An ab initio approach through the coupled perturbed Kohn-Sham scheme.

The coupled perturbed Kohn-Sham (CPKS) computational scheme for the evaluation of electric susceptibility tensors in periodic systems, recently implemented in the CRYSTAL code, has been extended to third-order. It is, then, used to obtain static electronic hyperpolarizabilities of zigzag BN nanotubes for the first time. This procedure, which is fully analytic in all key steps, requires a double self-consistent treatment for taking into account the first- and second-order response of the system to the applied field.

Ab initio modeling of TiO2 nanotubes.

TiO(2) nanotubes constructed from a lepidocrocite-like TiO(2) layer were investigated with ab initio methods employing the periodic CRYSTAL code. The dependence of strain energies, structural and electronic properties on the tube diameter was investigated in the 18-57 A range. Nanotubes constructed by a (0,n) rollup proved to be the most stable at all diameters.