Recent reactor neutrino experiments have shown anomalous results in both the reactor flux and spectrum measurements. Compared the measurements, Reactor neutrino flux shows a 6% deficit while reactor netrino spectrum illustrates a bump-like structure at around 5 MeV region. In this talk we will employ the methods of both theoretical model prediction and global neutrino data analysis to expore all the possible solutions to these reactor anomalies. We will discuss drawbacks of the model predictions within the standard model, and the possible new physics solution beyond the standard model.
 Ab initio calculations of reactor antineutrino fluxes with exact lepton wave functions, Dong-Liang Fang, Yu-Feng Li, Di Zhang,arXiv:2001.01689.
 New Realization of the Conversion Calculation for Reactor Antineutrino Fluxes
Yu-Feng Li, Di Zhang, arXiv:1904.07791.
 KATRIN bound on 3+1 active-sterile neutrino mixing and the reactor antineutrino anomaly, C. Giunti, Y.F. Li, Y.Y. Zhang, arXiv:1912.12956.
 Diagnosing the Reactor Antineutrino Anomaly with Global Antineutrino Flux Data, C. Giunti, Y. F. Li, B. R. Littlejohn, P. T. Surukuchi, arXiv:1901.01807.