A multifaceted first-principles approach utilizing density functional theory, evolutionary algorithms, and lattice dynamics was used to construct the phase diagram of silicon up to 4 TPa and 26000 K. These calculations predicted that (1) an anomalous sequence of face-centered cubic (fcc) to body-centered cubic (bcc) to simple cubic (sc) crystalline phase transitions occur at pressures of 2.87 TPa and 3.89 TPa, respectively, along the cold curve; (2) the orthorhombic phases of Imma and Cmce-16 appear on the phase diagram only when the anharmonic contribution to the Gibbs free energy is taken into account; and (3) a substantial change in the slope of the principal Hugoniot is observed if the anharmonic free energy of the cubic diamond phase is considered.