By performing constant-pressure variable-cell ab initio molecular dynamics simulations we find a quadrupolar orthorhombic structure of Pca2 1 symmetry for the broken symmetry phase (phase II) of solid H 2 at T 0 and P 110 150 GPa. We present results for the equation of state, lattice parameters, and vibronic frequencies, in very good agreement with experimental observations. Anharmonic quantum corrections to the vibrational frequencies are estimated using available data on H 2 and D 2 . We assign the observed modes to specific symmetry representations. [S0031-9007(97)02938-4] PACS numbers: 62.50.+p, 64.30.+t, The quest for the structure of the high-pressure phases of hydrogen is a long standing one. Early predictions of an insulator-metal transition [1] led to a large body of experimental and theoretical work during the past 60 years. Metallization was not found as promptly as initially expected [2] but a rich phase diagram emerged. The current picture of the low and room-temperature phase diagram up to ϳ230 GPa is essentially based on the optical studies performed in diamond anvil cell (DAC) devices during the past decade [3,4].There is a consensus that hydrogen exhibits at least three different solid molecular phases. (1) At low pressures (,110 GPa for para-H 2 ) the centers of the H 2 molecules crystallize into an hcp lattice, but zero-point motion overcomes rotational energy barriers, leading to a free-rotator phase (phase I).(2) Between 110 and 150 GPa intermolecular interactions freeze the molecular rotations into an ordered broken-symmetry phase (BSP, or phase II). (3) Above ϳ150 GPa, a third phase (H-A, or phase III) is attained, whose structure is unknown. Here we focus specifically on phase II. Optical measurements in phase II indicate the presence of two, possibly three, infrared (IR) active modes, in constrast to phase I, where only one mode is observed. Raman spectra show a single peak, at a frequency ϳ10% lower than that of the IR modes. A consistent picture of the structural and dynamical properties of phase II is still lacking [5].On the theoretical side, most of the existing work consists of static total energy calculations within the local density approximation (LDA). Zero-point energy (ZPE) of the protons has been in a few cases included a posteriori based on frozen phonon calculations [6,7]. Hexagonal close packed structures with two and four molecules per unit cell appear to be the strongest candidates for the ground structure of phase II, but the relative orientation of the molecules is uncertain [8][9][10]. Cubic structures were also suggested [9]. This uncertainty persists, due to the incomplete optimization of the lattice parameters and the a priori selection of a particular space group symmetry in static calculations. Ab initio molecular dynamics simulations, which do not rely on the choice of a specific space group, have also been reported [11,12]. However, these earlier attempts were hampered by a fixed simulation cell and a poor Brillouin zone (BZ) sampling.We report extensive ab init...