A first-principles approach called the self-consistent quasiharmonic approximation (SC-QHA) method is formulated to calculate the thermal expansion, thermomechanics, and thermodynamic functions of solids at finite temperatures with both high efficiency and accuracy. The SC-QHA method requires fewer phonon calculations than the conventional QHA method, and also facilitates the convenient analysis of the microscopic origins of macroscopic thermal phenomena. The superior performance of the SC-QHA method is systematically examined by comparing it with the conventional QHA method and experimental measurements on silicon, diamond, and alumina. It is then used to study the effects of pressure on the anharmonic lattice properties of diamond and alumina. The thermal expansion and thermomechanics of Ca3Ti2O7, which is a recently discovered important ferroelectric ceramic with a complex crystal structure that is computationally challenging for the conventional QHA method, are also calculated using the formulated SC-QHA method. The SC-QHA method can significantly reduce the computational expense for various quasiharmonic thermal properties especially when there are a large number of structures to consider or when the solid is structurally complex. It is anticipated that the algorithm will be useful for a variety of fields, including oxidation, corrosion, high-pressure physics, ferroelectrics, and high-throughput structure screening when temperature effects are required to accurately describe realistic properties.

Alloys are widely used in catalysts and structural materials. The nature of chemical bonding and the origin of alloy formation energies, defect energies, and interfacial properties have not been well understood to date but are critical to material performance. In this contribution, we explain the polar nature of chemical bonding and an implementation in classical and reactive atomistic simulations to understand such properties more quantitatively. Electronegativity differences between metal atoms lead to polar bonding, and exothermic alloy formation energies are related to charge transfer between the different elements. These differences can be quantified by atomic charges using pairwise charge increments, determined by matching the computed alloy formation energy to experimentally measured alloy formation energies using pair potentials for the pure metals. The polar character of alloys is comparable to organic molecules and partially ionic minerals, for example, AlNi and AlNi3 alloys assume significant atomic charges of ±0.40e and +0.60e/–0.20e, respectively. The subsequent analysis of defect sites and defect energies using force-field-based calculations shows excellent agreement with calculations using density functional theory and embedded atom models (EAM). The formation of vacancy and antisite defects is characterized by a redistribution of charge in the first shell of neighbor atoms in the classical models whereby electroneutrality is maintained and charge increments correlate with differences in electronegativity. The proposed atomic charges represent internal dipole and multipole moments, consistent with existing definitions for organic and inorganic compounds and with the extended Born model (Heinz, H.; Suter, U. W. J. Phys. Chem. B 2004, 108 (47), 18341–18352). The method can be applied to any alloy and has a reproducibility of ±10%. In contrast, quantum mechanical charge schemes remain associated with deviations exceeding ±100%. The atomic charges for alloys provide a simple initial measure for the internal electronic structure, surface adsorption of molecules, and reactivity in catalysis and corrosion. The models are compatible with the Interface force field (IFF), CHARMM, AMBER, OPLS-AA, PCFF, CVFF, and GROMOS for reliable atomistic simulations of alloys and their interfaces with minerals and electrolytes from the nanometer scale to the micrometer scale.