Background The benchmark dose (BMD) method calculates the dose associated with a specific change in response based on a specific dose-response relationship. Compared with the traditional no observed adverse effect level (NOAEL) method, the BMD method has many advantages, and the 95% lower confidence limit of benchmark dose lower limit (BMDL) is recommended to replace NOAEL in deriving biological exposure limits. No authority has yet published any health-based guideline for rare earth elements. Objective To evaluate genotoxicity threshold induced by acute exposure to neodymium nitrate in mice using BMD modeling through micronucleus test and comet assay. Methods SPF grade mice (n=90) were randomly divided into nine groups, including seven neodymium nitrate exposure groups, one control group (distilled water), and one positive control group (200 mg·kg−1 ethyl methanesulfonate), 10 mice in each group, half male and half female. The seven dose groups were fed by gavage with different concentrations of neodymium nitrate solution (male: 14, 27, 39, 55, 77, 109, and 219 mg·kg−1; female: 24, 49, 69, 97, 138, 195, and 389 mg·kg−1) twice at an interval of 21 h. Three hours after the last exposure, the animals were neutralized by cervical dislocation. The bone marrow of mice femur was taken to calculate the micronucleus rate of bone marrow cells, and the liver and stomach were taken for comet test. Results The best fitting models for the increase of polychromatophil micronucleus rate in bone marrow of female and male mice induced by neodymium nitrate were the exponential 4 model and the hill model, respectively. The BMD and the BMDL of female mice were calculated to be 31.37 mg·kg−1 and 21.90 mg·kg−1, and those of male mice were calculated to be 58.62 mg·kg−1 and 54.31 mg·kg−1, respectively. The best fitting models for DNA damage induced by neodymium nitrate in female and male mouse hepatocytes were the exponential 5 model and the exponential 4 model, respectively, and the calculated BMD and BMDL were 27.15 mg·kg−1 and 11.99 mg·kg−1 for female mice, and 16.28 mg·kg−1 and 10.47 mg·kg−1 for male mice, respectively. The hill model was the best fitting model for DNA damage of gastric adenocytes in both female and male mice, and the calculated BMD and BMDL were 36.73 mg·kg−1 and 19.92 mg·kg−1 for female mice, and 24.74 mg·kg−1 and 14.08 mg·kg−1 for male mice, respectively. Conclusion Taken the micronucleus rate of bone marrow cells, DNA damage of liver cells and gastric gland cells as the end points of genotoxicity, the BMDL of neodymium nitrate is 10.47 mg·kg−1, which can be used as the threshold of genotoxic effects induced by acute exposure to neodymium nitrate in mice.