Ground-state properties of metallic solids from ab initio coupled-cluster theory

Verena A. Neufeld,Hong-Zhou Ye,Timothy C. Berkelbach

DOI: https://doi.org/10.48550/arXiv.2204.01563

2022-04-04

Abstract:Metallic solids are a challenging target for wavefunction-based electronic structure theories and have not been studied in great detail by such methods. Here, we use coupled-cluster theory with single and double excitations (CCSD) to study the structure of solid lithium and aluminum using optimized Gaussian basis sets. We calculate the equilibrium lattice constant, bulk modulus, and cohesive energy and compare them to experimental values, finding accuracy comparable to common density functionals. Because the quantum chemical "gold standard" CCSD(T) (CCSD with perturbative triple excitations) is inapplicable to metals in the thermodynamic limit, we test two approximate improvements to CCSD, which are found to improve the predicted cohesive energies.

Materials Science,Chemical Physics,Computational Physics

What problem does this paper attempt to address?

This paper aims to solve the problem of calculating the ground - state properties of metallic solids, especially to study the structural and energy characteristics of two simple metals, lithium and aluminum, through ab initio coupled - cluster theory. Specifically, the paper mainly solves the following problems:

1. **Challenges in the Application of Wave Function Methods in Metallic Solids**:
Metallic solids are a challenging research object for wave - function - based electronic structure theories, because many successful wave - function methods rely on finite - order perturbation theory, and the correlation energies of these methods usually diverge in the thermodynamic limit. Therefore, these materials have not been studied in detail by these methods yet.

2. **Studying Metallic Solids Using the CCSD Method**:
The paper uses the coupled - cluster theory with single and double excitations (CCSD) to study the structures of body - centered cubic (BCC) lithium and face - centered cubic (FCC) aluminum. The authors calculated the equilibrium lattice constants, bulk moduli, and binding energies, and compared the results with experimental values, finding that the accuracy is comparable to that of the commonly used density functionals.

3. **Improving the CCSD Method to Increase Prediction Accuracy**:
Since the "gold standard of quantum chemistry" CCSD(T) is not applicable in the thermodynamic limit of metals, the authors tested two approximate improvement methods for CCSD and found that these improvements can increase the predicted binding energies.

4. **Solutions to Technical Problems**:
- **Basis Set Errors**: Optimize Gaussian - type orbital (GTO) basis sets to reduce the total energy and the condition number of the overlap matrix.
- **Finite - Size Errors**: Use dense Brillouin - zone sampling and twisted boundary conditions and perform extrapolation.

5. **Benchmarking and Validation**:
The authors verified the effectiveness and accuracy of the proposed methods by comparing with experimental data and the results of other theoretical methods.

### Formula Summary

- **Cost Function for Optimizing Basis Sets**:
\[
\text{cost} = E_{\text{HF}} + E^{(2)}_c + \gamma \ln(\text{cond}(S))
\]
where \(E_{\text{HF}}\) is the Hartree - Fock energy, \(E^{(2)}_c\) is the second - order Møller - Plesset perturbation theory (MP2) correlation energy, \(S\) is the periodic overlap matrix, and \(\gamma = 10^{- 4}E_h\).

- **Binding Energy Equation**:
\[
E_{\text{coh}} = E_{\text{total}} - N E_{\text{atom}}
\]
where \(E_{\text{total}}\) is the total energy of the system, \(N\) is the number of atoms, and \(E_{\text{atom}}\) is the energy of an isolated atom.

- **Birch - Murnaghan Equation of State**:
\[
E(V) = E_0 + \frac{9}{8}B_0 V_0\left[\left(\frac{V_0}{V}\right)^{\frac{7}{3}}-\left(\frac{V_0}{V}\right)^{\frac{5}{3}}\right]+\frac{9}{16}(B_0' - 4)B_0 V_0\left[\left(\frac{V_0}{V}\right)^{\frac{7}{3}}-\left(\frac{V_0}{V}\right)^{\frac{5}{3}}\right]^2
\]
where \(E(V)\) is the total energy at volume \(V\), \(E_0\) is the energy at the equilibrium volume \(V_0\), \(B_0\) is the bulk modulus, and \(B_0'\) is the derivative of the bulk modulus.

Through these methods, the paper has successfully solved the key problems in calculating the ground - state properties of metallic solids and provided a valuable reference for future research.