We explore the unknown structure of phosphorus in phase IV (P-IV phase) based on first-principles calculations using the metadynamics simulation method. Starting from the simple cubic structure, we find a new modulated structure of the monoclinic lattice. The modulation is crucial to the stability of the structure. Through refining the structure further by changing the modulation period, we find the structure whose x-ray powder diffraction pattern is in best agreement with the experimental pattern. We expect that the modulation period of the structure in the P-IV phase is very close to that found in this study and probably incommensurate. DOI: 10.1103/PhysRevLett.96.095502 PACS numbers: 61.50.Ah, 62.50.+p, 64.70.Kb Recent progress in high-pressure physics has enhanced our recognition of a wide variety of crystal structures. Development of high-pressure techniques has also enabled the identification of structures that are stabilized only in a narrow pressure range. Interesting structures were found unexpectedly through high-pressure experiments. For instance, modulated structures are often found in the highpressure phases of elements. Lattice modulations have been found in group Vb elements, including As, Sb, and Bi Scarcity of experimental ultrahigh-pressure data restricts high-pressure studies. Thus, researchers often encounter difficulties in the identification of a crystal structure on the basis of experimental data alone. A theoretical approach provides additional information on the same problem. First-principles theory for determining crystal structures is believed to be sufficiently accurate. However, the limitations of computational resources sometimes impede full structure searches.We will focus on the case of phase IV of the phosphorus (P-IV) phase. Observation of the P-IV phase was first reported by Akahama et al. [8] in 1999. In the sequence of pressure-induced transformations, the simple cubic (sc) phase (P-III) appears at 10 GPa at low temperature. Akahama et al. [8] reported the appearance of a simple hexagonal (sh) phase, i.e., the P-V phase, which stabilizes above 137 GPa, and an intermediate phase, i.e., the P-IV phase, between sc and sh on the basis of x-ray diffraction data. At even higher pressures, the bcc structure (P-VI) has been theoretically predicted [9] and later identified in an experiment at 262 GPa [10]. The structure of phase IV, however, has not been identified experimentally. Ordinary Rietveld analysis based on a knowledge of the monoclinic symmetry alone has not been successful, presumably owing to the complexity of the lattice. Thus, we must guess the crystal structure or a pseudocrystal.Several structures have been tested as candidate structures for P-IV. Ahuja considered a structure of space group Imma [11]. Ehlers and Christensen studied relative stability of the Ba-IV structure against sc and sh in the pressure range from 100 to 200 GPa [12]. The Ba-IV structure is a kind of modulated structure. Despite these extensive studies, the structure of P-IV remains unidentifi...