The behavior of metals at high pressure is of great importance to the fields of shock physics, geophysics, astrophysics, and nuclear materials. We study here bulk crystalline aluminum from firstprinciples at pressures up to 2500 GPa -soon within reach of laser-based experimental facilities. Our simulations use density-functional theory and density-functional perturbation theory in the localdensity and generalized-gradient approximations. Notably, the two different exchange-correlation functionals predict very similar results for the f cc → hcp, f cc → bcc, and hcp → bcc transition pressures, around 175 GPa, 275 GPa, and 380 GPa respectively. In addition, our results indicate that core overlaps become noticeable only beyond pressures of 1200 GPa. From the phonon dispersions of the fcc phase at increasing pressure, we predict a softening of the lowest transverse acoustic vibrational mode along the [110] direction, which corresponds to a Born instability of the fcc phase around 725 GPa.First-principles calculations have proved useful to the fields of geophysics, 1 astrophysics, 2 and nuclear materials. 3 Aluminum, being cubic close-packed and having no d-shell electrons, is a prototype for theoretical predictions and understanding the high-pressure behavior of simple metals. 4 Currently the National Ignition Facility 5 at LLNL is expected to achieve shockless compression 6 of metals up to 2000 GPa. This new facility may provide rapid advancements to high-pressure physics and could partner very successfully with theoretical studies.The equation of state (EOS) and phase stability of aluminum were first studied from first-principles in the early 1980s. 7,8,9 In all cases the predicted phase sequence was f cc → hcp → bcc, but predictions differed greatly in the transition pressures. Several other calculations within the local-density approximation (LDA) 10 or the generalized-gradient approximation (GGA) 11,12 have since then been performed, with a predicted static (i.e. without the phonon contribution) f cc → hcp transition pressure of 205 ± 20 GPa 10 in LDA and 170 GPa 11 and 192 GPa 12 in GGA. These discrepancies are more notable for the hcp → bcc transition pressure: 565 ± 60 GPa 10 in LDA versus 360 GPa 11 in GGA, leaving significant uncertainties open. Theoretical work on the vibrational properties of aluminum also suggests for the f cc → hcp transition a transition pressure higher than the static one. 11,12 Elastic properties 13,14 and the absolute strength under tension 15 have also been calculated; the latter results are of particular interest as they demonstrate the important role vibrational modes play in determining mechanical stability and suggest that shear failure modes are inherent in aluminum.Experimentally, the equation of state of aluminum at high pressures was studied by shock-compression 16 at pressures above the predicted maximum for the f cc → hcp phase boundary, 10 but a transition was not observed. However recent diamond anvil cell experiments observed a f cc → hcp transition at 217 ± 10 GPa 17 ...