In the present work, three-dimensional molecular dynamics simulation is carried out to elucidate the nanoindentation behaviors of CuZr Bulk metallic glasses (BMGs). The substrate indenter system is modeled using hybrid interatomic potentials including both many-body Finnis Sinclair (FS) and two-body Morse potentials. A spherical rigid indenter (diameter = 60(1 = 10-10 m)) is employed to simulate the indentation process. Three samples of BMGs including Cu25Zr75 , Cu50Zr50 , and Cu75Zr25 are designed and the metallic glasses are formed by rapid cooling from the melt state at about 2000 K. The radial distribution functions are analyzed to reveal the dynamical evolution of the structure of the atoms with different compositions and different cooling rates. The mechanical behavior can be well understood in terms of load-depth curves and Hardness-depth curves during the nanoindentation process. Our results indicate a positive linear relationship between the hardness and the Cu concentration of the BMG sample. To reveal the importance of cooling rate provided during the processing of BMGs, we investigate the indentation behaviors of Cu50Zr50 at three different quenching rates. Nanoindentation results and radial distribution function (RDF) curves at room temperature indicate that a sample can be made harder and more stable by slowing down the quenching rate. In the present work, three-dimensional molecular dynamics simulation is carried out to elucidate the nanoindentation behaviors of CuZr Bulk metallic glasses (BMGs). The substrate indenter system is modeled using hybrid interatomic potentials including both many-body Finnis Sinclair (FS) and two The three samples of BMGs including Cu25Zr75, Cu50Zr50, and Cu75Zr25 are designed and the metallic glasses are formed by rapid cooling from the melt state at about 2000 K. The radial distribution functions are analyzed to reveal the dynamical evolution of the structure of the atoms with different compositions and different cooling rates. The mechanical behavior can be well understood in terms of load- depth curves and Hardness-depth curves during the nanoindentation process. Our results indicate a positive linear relationship between the hardness and the Cu concentra tion of the BMG sample; To reveal the importance of cooling rate provided during the processing of BMGs, we investigate the indentation behaviors of Cu50Zr50 at three different quenching rates. Nanoindentation results and radial distribution function (RDF) curves at room temperature indicate that a sample can be made harder and more stable by slowing down the quenching rate.