A classic molecular dynamic melting and cooling approach was used to generate bulk amorphous ZrO2 using a reactive force field. The density and the structure characteristics of the generated bulk ZrO2 were analyzed and compared with literature simulations and experimental data. The density dependence on the cooling rate in the simulations was established. The generated bulk ZrO2 can further be used for theoretical studies of amorphous ZrO2 properties. (C) 2015 Elsevier B.V. All rights reserved.