Stability of hypothetical AgIICl2 polymorphs under high pressure, revisited: a computational study

A comparative computational study of stability of candidate structures for an as-yet unknown silver dichloride AgCl2 is presented. It is found that all considered candidates have a negative enthalpy of formation, but are unstable towards charge transfer and decomposition into silver(I) chloride and chlorine within the DFT and hybrid-DFT approaches in the entire studied pressure range. Within SCAN approach, several of the “true” AgIICl2 polymorphs (i.e. containing Ag(II) species) exhibit a region of stability below ca. 20 GPa. However, their stability with respect to aforementioned decomposition decreases with pressure by account of all three DFT methods, which suggests a limited possibility of high-pressure synthesis of AgCl2. Some common patterns in pressure-induced structural transitions observed in the studied systems also emerge, which further testify to an instability of hypothetical AgCl2 towards charge transfer and phase separation.

Additionally, the high-pressure, nanotubular polymorph of AgF 2 (referred to as HP-AgF 2 , space group Pbcn) was considered as a candidate 4 . Structures of selected polymorphs described above are presented in Fig. 1. Structures of all of the aforementioned candidate polymorphs are provided in Supplementary Information. Importantly, Ag(I)r and Ag(I)h polymorphs do not contain Ag(II) species and instead are composed of substructures of Ag I Cl and Cl 2 molecules. For the remaining five Ag II Cl 2 polymorphs containing paramagnetic d 9 silver cation, magnetic interactions were taken into account: • AgF 2 type: 2D antiferromagnetic (AFM) coupling within layers; • CuCl 2 type: AABB-type AFM coupling within chains, known to exist in CuCl 2 8 , and found to be the lowestenergy magnetic solution for AgCl 2 in this arrangement 6 ; • AgF 2 HP type: magnetic dimers coupled along ~ 180 degrees Ag-F-Ag bridges within nanotubes; • MnO 2 type: AABB-type AFM coupling-ferromagnetic (FM) between adjacent [AgCl 4 ] squares and AFM between pairs.
The AuCl 2 type derives from a structure of gold(I,III) chloride, a mixed-valence compound which consists of Au 4 Cl 8 molecules, each containing two Ag(I) and two Ag(III) species 9 . When this structure is taken as a starting point for geometry optimization, all three computational methods yield a nanotubular polymorph somewhat www.nature.com/scientificreports/ similar to the AgF 2 HP type. However, only the HSE06 approach is able to reproduce the mixed-valence nature of AuCl 2 and of the corresponding nanotubular polymorph of AgCl 2 derived from the former (as evidenced by two different coordination patterns of Ag sites in that solution) 6 . On the other hand, in the PBEsol + U and SCAN calculations, a magnetic model with dimers as in the AgF 2 HP type polymorphs was considered. Stability of the studied candidate structures was evaluated using two parameters: (a) enthalpy of formation (labelled henceforth as ΔH f ), according to a reaction: and (b) stability towards decomposition into AgCl and Cl 2 (ΔH r ), or more precisely, the enthalpy of reaction: Defined in this manner, both parameters indicate thermodynamic instability when positive. Figure 2 shows plots of stability of studied AgCl 2 candidate types in terms of ΔH r for all three methods used in this work. Enthalpy of Ag + Cl 2 mixture relative to AgCl + ½Cl 2 is also plotted for comparison.
The studied polymorphs exhibit negative (favorable) enthalpies of formation within the entire studied pressure range in all three computational approaches. This is indirectly visible in Fig. 2 as the fact that curves for those polymorphs lie below the curve for Ag + Cl 2 . On the other hand, they were found to be unstable in terms of ΔH r in the entire studied pressure range within PBEsol + U and HSE06. However, results of SCAN calculations indicate moderate stability of Ag II Cl 2 polymorphs below ca. 20 GPa. Ag(I)r polymorph is the most stable The initial set of calculations was performed within the PBEsol + U approach, as the least computationally demanding. It was found that the Ag(I)h type collapses upon compression to 10 GPa into the Ag(I)r structure and it was not considered any further in the analogous HSE06 and SCAN calculations. Similarly, MnO 2 type was also neglected past the PBEsol + U approach, as it was found to be the least stable in terms of ΔH r among the studied types. Therefore, these two polymorphs are not taken into account in Figs. 1, 2 and 3. Figure 3 compares relative volume (ΔV r ) of the studied polymorphs with respect to AgCl + ½Cl 2 mixture. Since the pV term becomes a considerable contribution to enthalpy at elevated pressures, the fact that most of the candidates for AgCl 2 considered here have a positive ΔV r can be seen as an important factor leading to their relative instability. However, it should be noted that only HSE06 accurately reproduces ambient-pressure volumes of AgCl and Cl 2 . PBEsol + U underestimates the volume of AgCl and Cl 2 , while SCAN-that of Cl 2 . Therefore, ΔV r values at 0 GPa in Fig. 3, and by extension-initial compressibilities-should be taken with a grain of salt. A noticeable dip at 10 GPa within the SCAN approach is likely a manifestation of this. The computationally demanding HSE06 results are likely to be the most correct.
Some of the features in Figs. 2 and 3, such as abrupt changes of relative energies or volumes, and apparent convergence of plots corresponding to different polymorphs are indicative of structural transitions. The nature and implication of those transition will be discussed in the next section.

Pressure-induced structural transitions.
As an introduction to analysis of structural transitions of AgCl 2 polymorphs, let us first discuss the Ag(I)r solution. As mentioned before, this structure is made up of subunits of rocksalt-type AgCl and of Cl 2 molecules. Within the studied pressure range, it undergoes structural www.nature.com/scientificreports/ rearrangements, which can be approximated as a sequence of deformations of the AgCl double layers, leading from a fundamentally NaCl-like coordination patterns to CsCl-like patterns, with increasing coordination number of Ag atoms. Importantly, it should be stressed again that this structure emerged as one of the lowest-energy solutions in an evolutionary algorithm structural search reported in the previous contribution 6 . It should not be treated as a viable candidate for the structure of AgCl 2 , but rather as a manifestation of the proclivity of the studied system towards charge transfer and phase separation. The case for instability of Ag II Cl 2 towards these processes is further strengthened by the fact that the Ag(I)r polymorph remains the most stable (with respect to ΔH r ) throughout the studied pressure range in both PBEsol + U and HSE06 approaches. A noticeable drop in ΔH r between 50 and 60 GPa for this polymorph in the SCAN approach was the reason for extending the studied pressure range to 100 GPa in this case. However, this drop is a result for NaCl-CsCl-like transition in the AgCl subphase, which is more abrupt than in analogous PBEsol + U and HSE06 calculations. No further phase transitions for Ag(I)r polymorph are observed up to 100 GPa. The CuCl 2 -type polymorph emerges as the most structurally robust in this study, as it does not undergo any collapse or drastic deformation in the studied pressure range, maintaining a relatively low ΔV r in all three methods. The only change to the structure of this polymorph occurs in terms of arrangement of chains relative to one another. Up to 10 GPa (in all three computational approaches), the chains are positioned as in Fig. 4a, resulting in octahedral coordination of Ag atoms (Fig. 4b). The octahedra are elongated by 29%, 29% and 25%, according to PBEsol, HSE06 and SCAN, respectively. By account of all three methods, the arrangement changes between 10 and 20 GPa, leading to a 4 + 4 coordination of Ag atoms, with the 4 inter-chain contacts longer by 23 to 33%, depending on the method (Fig. 4d). This is achieved in different ways: in PBEsol + U and HSE06 results, the chains move relative to one another in a direction perpendicular to direction of propagation (Fig. 4c). In SCAN approach, this is achieved through a change of one of the unit cell angles, which results in sliding the chains relative to each other. The resulting 4 + 4 coordination is the same in all cases (Fig. 4d), but the longer Cl contacts are aligned parallel (in SCAN) or perpendicular (in PBEsol + U and HSE06) to direction of propagation. Further compression to 30 GPa transforms the structure into that seen at 20 GPa in PBEsol + U and HSE06 pictures. All of these transitions can be seen as a means to achieve a more efficient packing of chains (as evidenced by increasing coordination number of Ag). Changes in local coordination of Ag atoms in this polymorph are plotted in the top panel of Fig. 6.
Another result of these rearrangements is a reduction of Cl…Cl distances between neighboring chains. These distances drop from ca. 3.8 to 2.7 Å between 0 and 60 GPa in HSE06, compared to 3.2 to 2.5 Å in PBEsol + U and 3.5 to 2.8 Å in SCAN. Recall that only HSE06 correctly reproduced the ambient pressure (low-temperature) volume of solid molecular Cl 2 at 0 GPa, so these results testify to the superiority of HSE06 in describing weak interactions compared to the other two methods utilized here (and free from the explicit van der Waals terms). Importantly, all of these distances are larger than the Cl-Cl bond in solid molecular Cl 2 , which remains at ca. 2.00 Å and contracts very little (less than 0.05 Å) with compression in results from all three methods. Additionally, this polymorph was further optimized at 100 GPa with PBEsol + U and SCAN methods, and does not undergo any structural modifications in that pressure range. AgF 2 -type polymorph undergoes substantial structural rearrangements with increasing pressure. At 0 GPa, it adopts a layered structure of Pbca space group (Fig. 5a), in which every Ag atom forms 4 in-layer bonds with Cl atoms, with two additional Cl atoms from adjacent layers together constituting a distorted octahedral coordination (Fig. 5b). The axial Cl contacts are noticeably longer (by 23%, 33% and 27% in PBEsol + U, HSE06 and SCAN, respectively) than equatorial ones. The same phenomenon is observed in AgF 2 and in general, elongated octahedral coordination is a well-documented phenomenon in silver(II) fluorides (which includes ternary compounds) 10 . This is usually attributed to Jahn-Teller effect, whereby a vibronic instability leads to elongation or contraction of bonds along one of the three axes of octahedron, which lowers the overall electronic energy. Upon compression, this elongation is reduced in AgCl 2 down to 1%-2% at 30 GPa in both PBEsol + U and SCAN  www.nature.com/scientificreports/ results. The structure at this point is reminiscent of high-pressure structure of PdF 2 11 , which-like AgF 2 -can be described as a distorted fluorite, but lacks Jahn-Teller elongation of [PdF 6 ] octahedra and is cubic (Pa-3). Further increase of pressure (40 GPa in PBEsol + U and 50 GPa in SCAN) leads to a structural transition into a polymorph consisting of 1D chains similar to those found in CuCl 2 -type polymorph (Fig. 5d). The symmetry of this collapsed chain polymorph corresponds to Pnma space group. In the HSE06 picture, the contraction is less pronounced-down to 10% at 30 GPa, and the transition to the chain polymorph occurs via a different layered structure observed at 40 GPa, where the layers become more corrugated and more separated from each other (Fig. 5c). Ag atoms retain an approximately octahedral coordination, but the axial Cl contacts are now within the same layer as [AgCl 4 ] squares. The elongation of octahedra due to Jahn-Teller effect is still noticeable (9%). Further compression to 50 GPa leads to a collapse to chains as in the other two methods. Changes in local coordination of Ag atoms in this polymorph are plotted in the bottom panel of Fig. 6.
The transitions described above can be seen as an abrupt drop in relative enthalpy in Fig. 2. In the previous contribution discussing relative stability of AgCl 2 candidate structures, AgF 2 -type emerged as the most likely candidate at higher pressures due to its comparatively low molar volume among the considered structures 6 .
However, it appears that compression of AgF 2 -type produces a lot of strain in the structure, as evidenced by a strong increase of ΔH r (Fig. 2), which is released through the aforementioned transition.
Arguably the most interesting pattern of pressure-induced transitions can be observed in the AuCl 2 -type nanotubular polymorph. The differences in outcomes of compression between the three computational approaches are the most pronounced for this system, although upon closer look we can identify their fundamental similarity. Recall that in the HSE06 approach, AuCl 2 -type is mixed-valent: Ag(I) species are connected to 3 Cl atoms in an approximately flat triangular pattern, while the Ag(III) sites appear as [AgCl 4 ] square units, which are analogous to those in AgF 3 12 . At 0 GPa, the triangular contacts average 2.53 Å, while the bonds in square units are 2.29 Å, which is even shorter than for Ag(II) in [AgCl 4 ] square in CuCl 2 -type and AgF 2 -type polymorphs at the same pressure. This supports the assignment of the sites as Ag(I) and as Ag(III), respectively. Compression to 15 GPa induces a change in local coordination of the Ag(I) species, which picks up 4 Cl atoms along the axis perpendicular to the plane of the former triangle, resulting in a 4 + 3 coordination, with an average bond length of 2.61 Å. Meanwhile, the Ag(III) subunit retains a square coordination with a shorter average bond length of 2.27 Å. Further compression to 30 GPa leads to a rearrangement of nanotubes, which are now formed by a different combination of Ag and Cl atoms; importantly, all Ag atoms are coordinated by 4 Cl atoms in an approximately square-planar manner, with an average for the formerly Ag(III) sites at 2.33 Å and the formerly Ag(I) ones-2.38 Å. This convergence of the two sites in terms of local coordination points to a comproportionation process, whereby all Ag sites are now nominally Ag(II) species. The transition described above is shown in Fig. 7.
A very different scenario is observed in the SCAN picture ( Fig. 8). At 0 GPa, the two Ag sites are already equivalent, but upon compression to 10 GPa, the local square-planar coordination of one of them is rotated by 90 degrees i.e., two of its four nearest Cl neighbors are substituted for another two, which leads to connections between nanotubes in the c direction. Further compression to 30 GPa results in an inward contraction of individual nanotubes. The local coordination of the Ag sites marked with red circle changes to more uniformly octahedral-a fifth Cl atom is picked up from a neighboring nanotube along the axis perpendicular to the [AgCl 4 ] plane. The average length of the five Ag-Cl bonds for this Ag site is 2.51 Å (2.45-2.55 Å). The sixth nearest neighbor, on the opposite side of the former [AgCl 4 ] is actually another Ag atom, located at a distance of 2.70 Å. (This connection is marked with a dashed line in Fig. 8 . This leads us to infer a formation of a weak Ag-Ag interaction in this polymorph, which is a very interesting finding and reminiscent of a recent work, where a silver subchloride Ag 8 Cl 6 was predicted, featuring [Ag 6 ] subunits within its structure 13 . Further compression above 50 GPa leads to a substantial structural rearrangement: AgCl-like chains with square crosssection are formed, interspersed with Cl 2 molecules.  www.nature.com/scientificreports/ A similar picture emerges from the PBEsol + U approach (Fig. 9). As in SCAN results, compression to 10 GPa leads to a rearrangement of local coordination of half of Ag sites. The subsequent contraction of nanotubes and formation of short Ag…Ag contacts is observed at a lower pressure of 20 GPa (compared to 30 GPa in SCAN). Upon further pressure increase to 30 GPa, a structure similar to Ag(I)r polymorph is formed, consisting of double layers of rocksalt-like AgCl interspersed with Cl atoms. The difference is that here, no discernible Cl-Cl molecules are formed-the distance between Cl atoms lying between AgCl layers is ca. 2.3 Å, compared to ca. 2.0 Å in solid molecular Cl 2 and in Ag(I)r polymorph. This is likely an artificial result, which will be further discussed in the "Electronic structure" section. Finally, compression to 40 GPa leads to rearrangement within AgCl layers, which increases the coordination number of Ag from 6 to 7 and the coordination environment resembles that in CsCl. Cl 2 molecules between the AgCl layers, characteristic of Ag(I)r polymorph, can also be discerned. Indeed, as can be seen in Fig. 2, this solution also converges with Ag(I)r polymorph in terms of ΔH r .
As in the case of Ag(I)r solutions, formation of separate domains of AgCl and Cl 2 during structural transitions of AuCl 2 -type in PBEsol + U and SCAN picture, can be interpreted as another manifestation of the system's tendency for AgCl + ½Cl 2 phase separation, rather than as viable structural candidates. However, it should be pointed out that sodium and potassium chlorides with exotic stoichiometries as e.g. NaCl or Na 3 Cl, have been predicted in the past 14,15 . Nanotubular AgF 2 -HP-type polymorph undergoes a structural collapse above 30 GPa in PBEsol + U results and above 60 GPa in SCAN results. The final structures do not resemble any of those discussed above; rather, they feature domains of connections between Ag and Cl atoms which do form any extended and discernible pattern, and are instead interspersed with Cl 2 molecules. Since those solutions are consistently very high in relative enthalpy (ΔH r ), they will not be further analyzed here. (However, an example of the resulting structures can be found in Supplementary Information).
Dynamical stability of candidate AgCl 2 structures. In addition to the analysis of relative enthalpy of hypothetical polymorphs of AgCl 2 , we also investigated stability of their crystal structures by calculating (at the PBEsol + U level) phonon frequencies in Γ point of the first Brillouin zone, in the same unit cells as presented above. We found that the Ag(I)r and AuCl 2 types exhibit no imaginary phonon frequencies at Γ point in the studied range from 0 to 60 GPa. CuCl 2 -type polymorph exhibits two imaginary frequencies only at 0 GPa, which, when followed along their normal coordinates and relaxed, lead to the same structure as the original solution at 0 GPa (within 0.01 meV/FU in terms of total energy), so these may be considered as artifacts. We have, however, found structural instabilities in AgF 2 -type and nanotubular HP-AgF 2 -type, which we discuss below. Figure 10 compares pressure dependence of ΔH r for the original (high symmetry) and phonon-relaxed (lower symmetry) structures in the two cases mentioned.
The layered AgF 2 -type polymorph is dynamically stable at 0 GPa, but at 10 GPa, an imaginary mode of B 1u symmetry emerges, which leads to a structure analogous to the previously reported non-centrosymmetric (Pca2 1 ), moderate-pressure polymorph of AgF 2 4 . This solution retains the layered arrangement, but the Ag atoms diverge from their position within the plane formed by surrounding Cl atoms (Fig. 11). The consequence of this transition is that the local coordination of Ag atoms changes from 4 + 2 (deformed octahedron) to 4 + 3 (Fig. 11), which likely minimizes repulsion between one of the lobes of the filled d(z 2 ) orbitals and ligands. We proceeded to optimize this new solution in the same pressure range as the original AgF 2 (0-60 GPa, every 10 GPa). At 0 GPa, the original Pbca remains more stable (by ca. 0.01 eV/FU in terms of ΔH r ). At 10 GPa, it becomes only slightly favored (by 0.004 eV/FU in terms of ΔH r ), and more favored at higher pressure up to 30 GPa. The Pca2 1 retains dynamic stability (no imaginary phonons in Γ) until 30 GPa. At higher pressure-like the original Pbca solution-it collapses into a 1D structure, which is higher in energy than the previously discussed Pnma structure resulting from the collapse of Pbca (Fig. 10). The new collapsed structure features chains of AgCl interspersed with chains of Cl atoms along the c direction (the structure can be found in ESI).
Similarly, nanotubular HP-AgF 2 type exhibits an imaginary mode of B 1g symmetry already at 0 GPa. Following a normal coordinate of that mode and relaxing the structure, we obtain a lower-symmetry nanotubular structure www.nature.com/scientificreports/ (P2 1 /c space group), which is more stable that the original in terms of ΔH r by 0.06 eV/FU already at 0 GPa. Likewise, we re-optimized this new phonon-relaxed solution in the 0-60 GPa range. The P2 1 /c solution remains more stable that the original Pbcn in terms of ΔH r up until its own collapse above 20 GPa (a drop between 20 and 30 GPa can be seen for the purple dashed line in Fig. 10). It also remains dynamically stable (lack of imaginary phonons in Γ) within that range. Figure 12 presents the change in local coordination experienced by Ag atoms upon aforementioned relaxation at 0 GPa. Importantly, for the Ag(I)r-type as well as CuCl 2 -type polymorphs at 0 GPa (i.e. the two most stable ones in terms of ΔH r ) we also calculated full phonon dispersion curves at the PBEsol + U level (in 2 × 2 × 2 and 1 × 2 × 2 supercells, respectively), which are shown in Supplementary Information. The results point to lack of imaginary phonons, which indicates metastability of these proposed structures. This leaves chance for the conceivable existence of these phases if prepared by some mildly exothermic reaction.
Electronic properties. As inferred from the previous paper 6 , antiferromagnetic superexchange in Ag II Cl 2 can, in principle, be expected to be strong, since Ag-Cl bonding in this hypothetical compound would likely be more covalent in nature than in its AgF 2 counterpart. Of course, AgCl-the only currently known binary combination of silver and chlorine-is an ionic solid, as is AgF. However, previous studies of AgF 2 demonstrated a covalent character of Ag-F bonding in that compound, evidenced by X-ray photoelectron spectroscopy 16 and by optical spectra 17 . In the former study, covalency increased in the sequence AgF → AgF 2 → AgF 3 . Therefore, it is reasonable to expect a similar trend in AgCl x compounds. On the other hand, Clanions are larger and more diffuse and are therefore softer Lewis bases than Fanions, which makes them more vulnerable to the strongly  www.nature.com/scientificreports/ oxidizing properties of Ag(II) cations. In addition, increasing pressure and the consequent reduction of interatomic distances increases orbital overlap, leading to broadening of electronic bands, which ultimately results in metallization of most known compounds (both ionic and covalent) 7 . Electronic properties of the studied AgCl 2 candidate polymorphs were scrutinized in terms of (a) magnetic moments on Ag atoms and (b) fundamental band gap at the Fermi level in electronic density of states (eDOS) graphs. As it turns out, changes in the two parameters are strongly correlated in that the pressure at which the band gap closes coincides with disappearance of magnetic moment on Ag atoms.
An example of eDOS plots-for AgF 2 -type and CuCl 2 -type-at 0 GPa and comparison between the three computational methods is presented in Fig. 10. One noticeable feature is the composition of conduction band. In principle, AgCl 2 , just like the known AgF 2 , is expected to be a charge-transfer insulator, where the band gap arises between filled nonmetal states and empty metal states (upper Hubbard band, UHB) 18 . While this is certainly the case in AgCl 2 , we can see a substantial admixing of Cl states to the conduction band, with an almost perfect overlap and approximately equal contributions from Ag and Cl states. This indicates a strong covalence of the Ag-Cl bonds that is comparable or indeed even stronger than in AgF 2 16 . It should also be pointed out that "insulator" in this case refers to a non-zero band gap resulting from electronic correlation as per the aforementioned Zaanen-Sawatzky-Allen model 18 . Clearly, with a band gap in the range 0.2-1.5 eV (depending on the method) (Fig. 13), the two AgCl 2 polymorphs in question can be more accurately described as semiconductors. Table 1 compares metallization pressure for selected polymorphs and between three methods. Metallization pressure is defined here as the lowest pressure point at which a solution exhibits null magnetic moments and null band gap. It can be seen that this value is the highest in HSE06 results. This coincides with the observation that the band gap is the largest for AgF 2 -type and CuCl 2 -type polymorphs in the HSE06 picture (Fig. 13). The differences between methods likely stem from the inclusion of exchange correlation in the hybrid-DFT-type functional, while in GGA-type PBEsol + U approach, electronic correlation, which is crucial for modelling open-subshell systems such as Ag(II) compounds, is only taken into account through U and J parameters (mentioned in "Methods" www.nature.com/scientificreports/ section). HSE06 results can likely be considered the most accurate, since-in general-hybrid functionals like HSE06 are very well suited for predictions of band gaps in solids (compared to GGA and even meta-GGA) 19 .
The phase-separated Ag(I)r polymorph remains insulating within the studied pressure range and by account of all three methods. Structural transitions into phase-separated solutions observed for AuCl 2 -type polymorph, which were discussed in the previous section, are also associated with reopening of the band gap. That is because constituent parts of these solutions-sublattices made up of ionic AgCl and of molecular Cl 2 -are insulators. The fact that the solution for AuCl 2 -type polymorph at 30 GPa in PBEsol + U picture retains metallic character after the afomentioned transition stems from the presence of dangling, unpaired Cl atoms on the outside of AgCl layers. Such arrangement should be unstable towards a Peierls distortion and formation of Cl 2 molecules, which is indeed what happens upon further compression to 40 GPa, and in SCAN picture. This testifies to the relatively poor suitability of PBEsol + U for describing electronic correlation.

Discussion
The Ag + Cl 2 phase diagram, which was studied here, consists of substances featuring a broad spectrum of chemical bonding types: metallic Ag, molecular Cl 2 , ionic AgCl, as well as a variety of more or less covalently bonded polymorphs of AgCl 2 , all of which have their own specific challenges when it comes to accurate theoretical description. Since AgCl 2 has not yet been observed experimentally, it is difficult to judge which of the three methods-GGA DFT (PBEsol + U), hybrid-DFT (HSE06) and meta-GGA DFT (SCAN)-is the most trustworthy for the assessment of stability of AgCl 2 at higher pressures-in particular, whether the region of stability seen in SCAN picture can be considered a reasonable finding. Our conclusions can to some extent be based on the results that these methods have provided when previously applied to other similar systems. To the best of our knowledge, a comparative computational study of an extended solid system comprised of a transition metal halide-like this one-has not been previously reported in literature. However, computations with hybrid functionals have been shown to provide excellent agreement with experimental data for e.g. ionic halides, where they are able to accurately reproduce lattice constants and bulk moduli of these solids 20 . In general, hybrid functionals appear to have the potential to be the most accurate among DFT methods, superior to both GGA and meta-GGA 21 . Additionally, in the case of our system, HSE06 functional is the only one of the three methods which correctly reproduces the volume and intermolecular distances in solid Cl 2 , which are determined by weak interactions between molecules.
Having said that, it is important to note that the overall picture which emerges from data presented here is remarkably consistent across the three computational methods, differing mostly in terms of pressure at which structural and electronic transitions are observed. Transitions of CuCl 2 and AgF 2 types reveal the tendency of the system to avoid repulsion between axial Cl atoms and filled d(z 2 ) orbital of Ag atoms-through relative displacement of chains in the former and through transition from a layered 2D structure into 1D chains in the latter. A similar tendency is seen in pressure-induced phase transitions of AgF 2 , where the high-pressure nanotubular structure can be viewed as a means both to maximize coordination number and to minimize the repulsion from the Ag d(z 2 ) lone pair 4 . The transitions of AuCl 2 -type and the relative stability of Ag(I)r polymorph firmly indicate that Ag II Cl 2 is unstable towards charge transfer and phase separation into AgCl and Cl 2 . Fundamentally, phase transitions observed in the studied polymorphs unfold upon pressure-induced decrease in distance between 1D or 2D structural constituents (chains, nanotubes, layers), which entails overcoming weak repulsive interactions between them. Among the three methods utilized here, the HSE06 functional, for reasons summarized above, is most likely best suited for description of those interactions.
Results reported in this work shed some light on properties and prospects of synthesis of the hypothetical AgCl 2 . Most importantly, our study shows that high pressure does not stabilize the candidate polymorphs derived from the previous study-in particular, we see that the prediction of stability of AgF 2 -type polymorph, which was based on a reasonable extrapolation from ambient-pressure data, nevertheless turned out not to be false 6 .
Of course, this work does not exhaust the list of possible candidates for the structure of AgCl 2 . Although the studied candidates generally retain a positive enthalpy with respect to decomposition into AgCl + ½Cl 2 , it is worth noting that the CuCl 2 -type chain structure does not collapse into phase-separated polymorph even at 100 GPa in PBEsol + U and SCAN results. Analogous HSE06 calculations at 100 GPa were not performed, but based on the observation made here that transition pressures are reliably the highest within HSE06 approach, it is reasonable to predict that such collapse would not be seen in HSE06 picture at 100 GPa, either. The apparent lack of such transition pathway could mean that, in principle, obtaining it as a metastable phase could be possible Table 1. Pressure of metallization for different polymorphs of AgCl 2 . a Although transformed into chains similar to CuCl 2 -type, the polymorph retains residual magnetic moment and non-zero band gap at the maximal studied pressure of 60 GPa. Magnetic coupling is FM intra-chain and AFM inter-chain. b Reopening of band gap occurs at higher pressures-see text. c This polymorph is mixed-valent in HSE06 picture and is non-magnetic at 0 GPa.  22,23 , all consistently suggest that Ag II Cl 2 , if ever obtained, would likely be metastable. Although the current study was focused on AgCl 2 stoichiometry only, the results of this study permit us to extrapolate observed trends towards the AgCl 3 . The tendency of AgCl 2 stoichiometry to undergo decomposition to (AgCl)(Cl 2 ) ½ may suggest that even at larger Cl contents, i.e. for AgCl 3 stoichiometry, one will observe phase separation to (AgCl)(Cl 2 ), or Ag + (Cl 3 -) 24 . A similar result is indicated by theoretical study for isolated AgCl 3 molecules in the gas phase 25 .

Methods
Calculations were carried out using VASP software [26][27][28][29][30] . Overall, three different computational methods were utilized (underlined are the names by which they are referred to throughout this work): I. PBEsol + U approach: GGA-type Perdew-Burke-Ernzerhof functional adapted for solids (PBEsol) 31 was used, additionally taking into account Coulombic interactions between d electrons through U and J parameters 32 explicitly set to 5 eV and 1 eV, respectively, and with correction for van der Waals interactions 33 . Plane-wave cutoff energy was set to 800 eV. k-space sampling of ca. 2π × 0.04 Å −1 was used for optimization and a denser k-spacing of ca. 2π × 0.03 Å −1 was used for electronic density of states (eDOS) calculations. On an example structure of CuCl 2 -type, this corresponds to 2 × 7 × 5 and 3 × 11 × 7 grid, respectively. Self-consistent-field convergence criterion was set to 10 -7 eV. Additionally at this level, dynamical stability of resulting structures was assessed by calculating phonon frequencies at Gamma point of the Brillouin zone. Structures were further optimized to minimize forces acting on atoms to ≤ 10 -5 eV/Å, and then Hessian matrix was constructed through finite differences approach as implemented in VASP. Phonon dispersion curves for selected systems were calculated and generated using PHONOPY interface 34,35 . PBEsol + U is the cheapest among the three methods, but has been successfully utilized for prediction of high-pressure structures of Ag(II) compounds in the past 2-4,36 . II. HSE06 approach: Hybrid-DFT HSE06 functional was utilized 37 . Due to higher computational load of this method, a coarser k-space sampling of ca. 2π × 0.05 Å −1 was used (e.g. 4 × 4 × 3 grid in AgF 2 -type) and plane-wave cut-off was set to 520 eV, with self-consistent-field convergence criterion of 10 -7 eV. Although the cutoff energy is lower than in the other two methods, our convergence tests have shown that increasing it to e.g. 600 eV does not lead to a substantial improvement (i.e. the resulting total energy change is ca. 0.0001 eV/FU, which is of no practical consequence in our data). HSE06 is the most computationally demanding of the three methods, but provides the most accurate description of electronic properties 19 . III. SCAN approach: Meta-GGA-type, strongly constrained and appropriately normed (SCAN) functional was used 38 , with correction for van der Waals interactions 39 . Plane-wave cut-off energy was set to 800 eV and k-space sampling of ca. 2π × 0.04 Å −1 was used, with self-consistent-field convergence criterion of 5 × 10 -7 eV. A denser k-spacing of ca. 2π × 0.03 Å −1 was used for eDOS calculations. SCAN approach is slightly more expensive than PBEsol + U, but still much cheaper than HSE06. This method has previously provided a value of antiferromagnetic coupling constant in AgF 2 that is in best agreement with experimental data 2 .
The studied pressure range was 0 to 60 gigapascals (GPa), additionally extended to 100 GPa for SCAN calculations. Integration grids were not substantially increased in metallic solutions at higher pressures, but our convergence tests showed that an increase of grid density does not lead to a meaningful change in total energy (within 0.005 eV/FU) or in lattice constants after re-optimization (within 0.01 Å) at denser mesh. Pressure step was 10 GPa up to 60 GPa (or in some cases 15 GPa in HSE06 approach) and 20 GPa above 60 GPa. Cl 2 was considered in its solid polymorph with Cmca space group, which is known to be stable in the entire pressure range considered here 40 . Similarly, metallic silver is also stable in its fcc structure within that range 41 . Known phase transitions of AgCl were taken into account when calculating relative enthalpy of AgCl 2 polymorphs 42 . A primitive unit cell of the KOH-type, high-pressure polymorph of AgCl was utilized in calculations, since it can be used to accurately describe the continuous nature of NaCl-KOH-TlI-CsCl sequence of phase transitions of AgCl 42,43 .