“…During the dynamics simulation, the exchange–correlation energy was approximated by generalized gradient approximation (GGA) and Perdew‐Burke‐Ernzerhof (PBE) functional, which calculated the electrostatic, exchange, and induced interaction energies well, [ 23,24 ] and GGA–PBE functional was used to explore the structure and properties of alloy clusters in many studies on clusters. [ 25–27 ] The system was set to canonical ensemble (NVT) system, the temperature was 300 K, the time step was set to 1 fs, and the total simulation time was 100 ps. In the process of geometry optimization and properties calculation, GGA–BLYP functional was used to approximate the exchange correlation energy.…”