Linear-response calculations within the linearized augmented plane-wave method

Calculation of phonons in real-space density functional theory

We present an accurate and efficient formulation for the calculation of phonons in real-space Kohn-Sham density functional theory. Specifically, employing a local exchange-correlation functional, norm-conserving pseudopotential in the Kleinman-Bylander representation, and local form for the electrostatics, we derive expressions for the dynamical matrix and associated Sternheimer equation that are particularly amenable to the real-space finite-difference method, within the framework of density functional perturbation theory. In particular, the formulation is applicable to insulating as well as metallic systems of any dimensionality, enabling the efficient and accurate treatment of semi-infinite and bulk systems alike, for both orthogonal and nonorthogonal cells. We also develop an implementation of the proposed formulation within the high-order finite-difference method. Through representative examples, we verify the accuracy of the computed phonon dispersion curves and density of states, demonstrating excellent agreement with established plane-wave results.

OpenCL-accelerated first-principles calculations of all-electron quantum perturbations on HPC resources

We have proposed, for the first time, an OpenCL implementation for the all-electron density-functional perturbation theory (DFPT) calculations in FHI-aims, which can effectively compute all its time-consuming simulation stages, i.e., the real-space integration of the response density, the Poisson solver for the calculation of the electrostatic potential, and the response Hamiltonian matrix, by utilizing various heterogeneous accelerators. Furthermore, to fully exploit the massively parallel computing capabilities, we have performed a series of general-purpose graphics processing unit (GPGPU)-targeted optimizations that significantly improved the execution efficiency by reducing register requirements, branch divergence, and memory transactions. Evaluations on the Sugon supercomputer have shown that notable speedups can be achieved across various materials.

Identification of local phase of nanoscale BaTiO₃ powders by high-resolution electron energy loss spectroscopy

The electron energy loss spectroscopy (EELS) technique was applied to investigate the local variation in the phase of barium titanate (BaTiO₃) ceramics. It was found that the fine structure of the titanium L₂,₃ edge and their satellite peaks were sensitively varied with the tetragonal-cubic phase transition. The peak splitting of Ti-L₃ edge of tetragonal-phased BaTiO₃ ceramics was widened because of the increased crystal field effect compared with that of cubic-phased BaTiO₃. In case of nanoscale BaTiO₃ powders, the L₃ edge splitting of the core region was found to be smaller than that of the shell region. The energy gap between peaks t₂g and eg varied from 2.36 to 1.94 eV with changing the probe position from 1 to 20 nm from the surface. These results suggest that the EELS technique can be used to identify the local phase of sintered BaTiO₃ ceramics.

An introduction to “Computational Crystallography”

The currently available methods for the computation of structures and their properties are reviewed. After a brief introduction into some common technical aspects, the capabilities and limitations of the most commonly used approaches are discussed. Examples are given to show the state of the art in Computational “Crystallography”, and possible future developments are outlined

Linear response calculations of lattice dynamics in strongly correlated systems

We introduce a new linear response method to study the lattice dynamics of materials with strong correlations. It is based on a combination of dynamical mean field theory of strongly correlated electrons and the local density functional theory of electronic structure. We apply the method to study the phonon dispersions of Mott insulators NiO and MnO in their paramagnetic insulating state not accessible by local density functionals. Our results are in good agreement with experiment.