Computational searches for structures of solid oxygen under pressures in the multi TPa range have been carried out using density-functional-theory methods. We find that molecular oxygen persists to about 1.9 TPa at which it transforms into a semiconducting square spiral-like polymeric structure (I41/acd) with a band gap of ∼3.0 eV. Solid oxygen forms a metallic zig-zag chain-like structure (Cmcm) at about 3.0 TPa, but the chains in each layer gradually merge as the pressure is increased and a structure of F mmm symmetry forms at about 9.5 TPa in which each atom has four nearest neighbors. The superconducting properties of molecular oxygen do not vary much with compression, although the structure becomes more symmetric. The electronic properties of oxygen have a complex evolution with pressure, swapping between insulating, semiconducting and metallic. Among the simple molecules studied at high pressures, oxygen has attracted particular attention due to its fundamental importance and intriguing properties. [1][2][3][4][5] For example, oxygen is the only known elemental molecule that exhibits magnetism, which adds substantial complexity to its phase diagram. When liquid oxygen is cooled at ambient pressure it undergoes a sequence of transitions to the γ, β and α solid phases at 54.39, 43.76 and 23.88 K, respectively [6]. Upon increase of pressure the monoclinic α-O 2 (C2/m) phase transforms into an orthorhombic δ-O 2 (F mmm) phase at about 3 GPa and to -O 2 at about 10 GPa. The structure of -O 2 has only recently been solved by x-ray diffraction (XRD) studies of single crystal [3] and powder [4] samples, although its vibrational spectrum was reported more than 20 years earlier [7]. The unit cell of -O 2 has C2/m symmetry and contains four O 2 molecules forming O 8 units. Both α and δ-O 2 are antiferromagnetic and the magnetic collapse at the δ-transition was predicted using molecular dynamics simulations [8]. This breakdown of the long-range antiferromagnetic order at about 8 GPa was recently observed in a neutron scattering experiment [2]. Density functional theory (DFT) studies have found another chain-like structure to be energetically slightly more favorable than the O 8 structure [9,10], although it has not been observed in experiments.Recent experiments [5] have shown that insulating -O 2 remains stable up to about 96 GPa before undergoing a continuous displacive and isosymmetric transition to the ζ phase, in agreement with earlier predictions [10]. The metallic ζ phase [11,12] has C2/m symmetry [2] and superconducts at temperatures below 0.6 K [1]. Manybody perturbation theory GW calculations [13,14] have suggested that the metal-insulator transition occurs at lower pressures than the measured -ζ transition pressure of 96 GPa. The above information, however, covers only a small part of the phase diagram and rather little is known about pure oxygen at pressures above 100 GPa.A previous DFT study reported that oxygen molecules persist to at least 250 GPa [10]. It is interesting to speculate about the highest p...