Aditi Mallik,Keith Runge,James W. Dufty,Hai-Ping Cheng
Abstract:The approximate representation of a quantum solid as an equivalent composite semi-classical solid is considered for insulating materials. The composite is comprised of point ions moving on a potential energy surface. In the classical bulk domain this potential energy is represented by pair potentials constructed to give the same structure and elastic properties as the underlying quantum solid. In a small local quantum domain the potential is determined from a detailed quantum calculation of the electronic structure. The primary new ingredients are 1) a determination of the pair potential from quantum data for equilibrium and strained structures, 2) development of pseudo-atoms for a realistic treatment of charge densities where bonds have been broken to define the quantum domain, and 3) inclusion of polarization effects on the quantum domain due to its environment. This formal structure is illustrated in detail for an silica nanorod. For each configuration considered, the charge density of the entire solid is calculated quantum mechanically to provide the reference by which to judge the accuracy of the <a class="link-external link-http" href="http://modeling.It" rel="external noopener nofollow">this http URL</a> is then shown that the quantum rod, the rod constructed from the classical pair potentials, and the composite classical/quantum rod all have the same equilibrium structure and response to elastic strain. The accuracy of the modeling is shown to apply for two quite different quantum chemical methods for the underlying quantum mechanics: transfer Hamiltonian and density functional methods.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: How to accurately describe complex phenomena in materials through the combination of quantum mechanics (QM) and classical mechanics (CM) in multi - scale modeling. Specifically, the paper focuses on insulating materials (such as silica) and proposes a new multi - scale modeling scheme, aiming to solve the following problems:
1. **Smooth transition between the quantum and classical domains**:
- In multi - scale modeling, different parts of the system need to be described with methods of different precisions. For example, in certain regions of the material (such as crack tips), the breaking of bonds between atoms causes significant deformation of the electron cloud, and detailed quantum - mechanical methods are required; while in areas far from these regions, classical - mechanical methods can be used. How to establish a smooth and accurate connection between these two different descriptions is a challenge.
2. **Treatment of dangling bonds and charge density**:
- When dividing the quantum and classical domains, covalent bonds may be cut, resulting in dangling bonds in the quantum domain, which in turn affects the accuracy of the charge density. The paper proposes to use "pseudo - atoms" and polarization effects to deal with these problems, ensuring that the charge density and force field in the quantum domain can be accurately reproduced.
3. **Construction of classical potential functions**:
- In order to make the description of the classical domain as close as possible to the real physical behavior, the paper proposes a new method to construct classical potential functions. These potential functions are trained based on data obtained from quantum - mechanical calculations to ensure that they can accurately reproduce the structure and elastic properties of the material.
4. **Model verification**:
- The paper verifies the effectiveness of the proposed multi - scale modeling scheme through a series of tests. For example, for silica nanorods, the paper shows that the structural and elastic properties of the quantum - mechanical, pure classical - mechanical, and composite quantum / classical - mechanical models are almost the same in the near - equilibrium state, proving the reliability of the model.
### Formula summary
- **Hamiltonians of ions and electrons**:
\[
H = H_i + H_e + U_{ie}
\]
where \(H_i\) and \(H_e\) are the Hamiltonians of ions and electrons in an isolated system respectively, and \(U_{ie}\) is the interaction term between ions and electrons:
\[
U_{ie}=\int dr dr' \frac{\rho_i(r)\rho_e(r')}{|r - r'|}
\]
- **Classical ionic potential energy function**:
\[
V_{ic}=\frac{1}{2}\int_C dr \rho_i(r)\left(\int_C dr' \frac{1}{|r - r'|}(\rho_i(r')-\delta(r - r')+\rho_{ec}(r',t))+\int_Q dr' \frac{1}{|r - r'|}(\rho_i(r')+\rho_{eq}(r',t))\right)
\]
- **Ionic potential energy function in the quantum domain**:
\[
V_{iq}=\frac{1}{2}\int_Q dr \rho_i(r)\left(\int_Q dr' \frac{1}{|r - r'|}(\rho_i(r')-\delta(r - r')+\rho_{eq}(r,t))+\int_C dr' \frac{1}{|r - r'|}(\rho_i(r')+\rho_{ec}(r',t))\right)
\]
Through these formulas and methods, the paper successfully solves some key problems in multi - scale modeling and provides a valuable reference for future research.