Zoledronic acid (ZL) has been used widely for treating skeletal diseases because of its high potency in inhibiting bone resorption. A detailed understanding of its physicochemical characteristics may be of great significance in both medicinal chemistry and structural biology for the design of novel bisphosphonates with higher activity. In the present work, the monoclinic (IM) and triclinic (IT) polymorphs of ZL in the gas phase and the aqueous phase were studied by density functional theory (DFT) method at the B3LYP/6-311++G** level. The polarizable continuum model (PCM) was employed to study the solvent effect on structures and properties. The optimized IM and IT conformations in both phases are in reasonable agreement with the experimental structures with the overall mean absolute percent deviation (MAPD%) less than 3.1 %. The presence of intramolecular hydrogen bond within both conformations was identified in the solvent. The IR spectra were simulated and assigned in detail, which agreed well with the experimental data. The intramolecular hydrogen bonding interactions resulted in the shift of vibrational frequencies of hydroxyl to the low band by 12-22 cm(-1) and 24-26 cm(-1) for IM and IT conformations, respectively. Their thermodynamic properties were also calculated based on the harmonic vibrational analysis, including standard heat capacity (C(°)p,m), entropy (S(°)m), and enthalpy (H(°)m). The molecular stability, hydrogen bonding interaction and other electronic properties have been further analyzed by the natural bond orbital (NBO), atoms in molecules (AIM), molecular electrostatic potential (MEP) and frontier molecular orbital (FMO) analysis.