Nuclear quantum effects on the high pressure melting of dense lithium

Data-driven prediction of complex crystal structures of dense lithium

Lithium (Li) is a prototypical simple metal at ambient conditions, but exhibits remarkable changes in structural and electronic properties under compression. There has been intense debate about the structure of dense Li, and recent experiments offered fresh evidence for yet undetermined crystalline phases near the enigmatic melting minimum region in the pressure-temperature phase diagram of Li. Here, we report on an extensive exploration of the energy landscape of Li using an advanced crystal structure search method combined with a machine-learning approach, which greatly expands the scale of structure search, leading to the prediction of four complex Li crystal structures containing up to 192 atoms in the unit cell that are energetically competitive with known Li structures. These findings provide a viable solution to the observed yet unidentified crystalline phases of Li, and showcase the predictive power of the global structure search method for discovering complex crystal structures in conjunction with accurate machine learning potentials.

Nuclear Quantum Effects on the Charge-Density Wave Transition in NbX2 (X = S, Se)

Understanding the origin of charge-density wave (CDW) instability is important for manipulating novel collective electronic states. Many layered transition metal dichalcogenides (TMDs) share similarity in the structural and electronic instability, giving rise to diverse CDW phases and superconductivity. It is still puzzling that even isostructural and isoelectronic TMDs show distinct CDW features. For instance, bulk NbSe₂ exhibits CDW order at low temperature, while bulk NbS₂ displays no CDW instability. The CDW transitions in single-layer NbS₂ and NbSe₂ are also different. In the classic limit, we investigate the electron correlation effects on the dimensionality dependence of the CDW ordering. By performing ab initio path integral molecular dynamics simulations and comparative analyses, we further revealed significant nuclear quantum effects in these systems. Specifically, the quantum motion of sulfur anions significantly reduces the CDW transition temperature in both bulk and single-layer NbS₂, resulting in distinct CDW features in the NbS₂ and NbSe₂ systems.

Prediction of Stable Ground-State Binary Sodium-Potassium Interalkalis under High Pressures

The complex structures and electronic properties of alkali metals and their alloys provide a natural laboratory for studying the interelectronic interactions of metals under compression. A recent theoretical study (J. Phys. Chem. Lett. 2019, 10, 3006) predicted an interesting pressure-induced decomposition-recombination behavior of the Na₂K compound over a pressure range of 10-500 GPa. However, a subsequent experiment (Phys. Rev. B 2020, 101, 224108) reported the formation of NaK rather than Na₂K at pressures above 5.9 GPa. To address this discordance, we study the chemical stability of different stoichiometries of Na_xK (x = 1/4, 1/3, 1/2, 2/3, 3/4, 4/3, 3/2, and 1-4) by an effective structure searching method combined with first-principles calculations. Na₂K is calculated to be unstable at 5-35 GPa due to the decomposition reaction Na₂K → NaK + Na, coinciding well with the experiment. NaK undergoes a combination-decomposition-recombination process accompanied by an opposite charge-transfer behavior between Na and K with pressure. Besides NaK, two hitherto unknown compounds NaK₃ and Na₃K₂ are uncovered. NaK₃ is a typical metallic alloy, while Na₃K₂ is an electride with strong interstitial electron localization.

Decomposition and Recombination of Binary Interalkali Na2K at High Pressures

Under compression, "'simple'" alkali metals and their alloys exhibit complex structural and electronic properties, leading to fundamental interest in their high-pressure behaviors. Here, the swarm-intelligence structure-searching method was employed to identify the high-pressure phases of binary interalkali Na₂K, which has long been known to possess a MgZn₂-Laves phase at ambient pressure, but the high-pressure behavior remains elusive. We uncovered four new structures over a pressure range of 10-500 GPa, although the compound was found to become unstable upon decomposition into Na and K from 37 to 273 GPa. In phases before decomposition, the electrons were gradually delocalized with an increase in pressure and there was charge transfer from K to Na, whereas in phases after recombination, the electrons were gradually localized into the interstitials of the crystals, showing the unexpected opposite trend of charge transfer from Na to K, remarkably, where K was found to exhibit an oxidation state beyond the -1 valence state. The results can improve our understanding of the interaction and evolution of s electrons under compression.

Proton Migration in Hybrid Lead Iodide Perovskites: From Classical Hopping to Deep Quantum Tunneling

The organic-inorganic halide perovskites (OIHPs) have shown enormous potential for solar cells, while problems like the current-voltage hysteresis and the long-term instability have seriously hindered their applications. Ion migrations are believed to be relevant. But the atomistic details still remain unclear. Here we study the migrations of ions in CH₃NH₃PbI₃ (MAPbI₃) at varying temperatures ( T's), using combined experimental and first-principle theoretical methods. Classical hopping of the iodide ions is the main migration mechanism at moderate T's. Below ∼270 K, the kinetic constant for ionic migration still shows an Arrenhius dependency, but the much lower activation energy is attributed to the migration of H⁺. A gradual classical-to-quantum transition takes place between ∼140 and ∼80 K. Below ∼80 K, the kinetic constant becomes T-independent, suggesting that deep quantum tunneling of H⁺ takes over. This study gives direct experimental evidence for the migrations of H⁺s in MAPbI₃ and confirms their quantum nature.

The collective and quantum nature of proton transfer in the cyclic water tetramer on NaCl(001)

Proton tunneling is an elementary process in the dynamics of hydrogen-bonded systems. Collective tunneling is known to exist for a long time. Atomistic investigations of this mechanism in realistic systems, however, are scarce. Using a combination of ab initio theoretical and high-resolution experimental methods, we investigate the role played by the protons on the chirality switching of a water tetramer on NaCl(001). Our scanning tunneling spectroscopies show that partial deuteration of the H₂O tetramer with only one D₂O leads to a significant suppression of the chirality switching rate at a cryogenic temperature (T), indicating that the chirality switches by tunneling in a concerted manner. Theoretical simulations, in the meantime, support this picture by presenting a much smaller free-energy barrier for the translational collective proton tunneling mode than other chirality switching modes at low T. During this analysis, the virial energy provides a reasonable estimator for the description of the nuclear quantum effects when a traditional thermodynamic integration method cannot be used, which could be employed in future studies of similar problems. Given the high-dimensional nature of realistic systems and the topology of the hydrogen-bonded network, collective proton tunneling may exist more ubiquitously than expected. Systems of this kind can serve as ideal platforms for studies of this mechanism, easily accessible to high-resolution experimental measurements.

A method of solid-solid phase equilibrium calculation by molecular dynamics

A method for evaluation of solid-solid phase equilibrium curves in molecular dynamics simulation for a given model of interatomic interaction is proposed. The method allows to calculate entropies of crystal phases and provides an accuracy comparable with that of the thermodynamic integration method by Frenkel and Ladd while it is much simpler in realization and less intense computationally. The accuracy of the proposed method was demonstrated in MD calculations of entropies for EAM potential for iron and for MEAM potential for beryllium. The bcc-hcp equilibrium curves for iron calculated for the EAM potential by the thermodynamic integration method and by the proposed one agree quite well.

Predicted reentrant melting of dense hydrogen at ultra-high pressures

The phase diagram of hydrogen is one of the most important challenges in high-pressure physics and astrophysics. Especially, the melting of dense hydrogen is complicated by dimer dissociation, metallization and nuclear quantum effect of protons, which together lead to a cold melting of dense hydrogen when above 500 GPa. Nonetheless, the variation of the melting curve at higher pressures is virtually uncharted. Here we report that using ab initio molecular dynamics and path integral simulations based on density functional theory, a new atomic phase is discovered, which gives an uplifting melting curve of dense hydrogen when beyond 2 TPa, and results in a reentrant solid-liquid transition before entering the Wigner crystalline phase of protons. The findings greatly extend the phase diagram of dense hydrogen, and put metallic hydrogen into the group of alkali metals, with its melting curve closely resembling those of lithium and sodium.

Nuclear quantum effects of hydrogen bonds probed by tip-enhanced inelastic electron tunneling

We report the quantitative assessment of nuclear quantum effects on the strength of a single hydrogen bond formed at a water-salt interface, using tip-enhanced inelastic electron tunneling spectroscopy based on a scanning tunneling microscope. The inelastic scattering cross section was resonantly enhanced by "gating" the frontier orbitals of water via a chlorine-terminated tip, so the hydrogen-bonding strength can be determined with high accuracy from the red shift in the oxygen-hydrogen stretching frequency of water. Isotopic substitution experiments combined with quantum simulations reveal that the anharmonic quantum fluctuations of hydrogen nuclei weaken the weak hydrogen bonds and strengthen the relatively strong ones. However, this trend can be completely reversed when a hydrogen bond is strongly coupled to the polar atomic sites of the surface.

The phase diagram of high-pressure superionic ice

Superionic ice is a special group of ice phases at high temperature and pressure, which may exist in ice-rich planets and exoplanets. In superionic ice liquid hydrogen coexists with a crystalline oxygen sublattice. At high pressures, the properties of superionic ice are largely unknown. Here we report evidence that from 280 GPa to 1.3 TPa, there are several competing phases within the close-packed oxygen sublattice. At even higher pressure, the close-packed structure of the oxygen sublattice becomes unstable to a new unusual superionic phase in which the oxygen sublattice takes the P2₁/c symmetry. We also discover that higher pressure phases have lower transition temperatures. The diffusive hydrogen in the P2₁/c superionic phase shows strong anisotropic behaviour and forms a quasi-two-dimensional liquid. The ionic conductivity changes abruptly in the solid to close-packed superionic phase transition, but continuously in the solid to P2₁/c superionic phase transition.

At high pressure, water forms superionic ice with an oxygen lattice and melted liquid hydrogens, which could exist on ice-rich planets. Here, Sun et al. predict a new phase of superionic ice where the hydrogens preferentially diffuse in two-dimensions within oxygen superlattice with the P2₁/c symmetry.

Chemical bonding in hydrogen and lithium under pressure

Though hydrogen and lithium have been assigned a common column of the periodic table, their crystalline states under common conditions are drastically different: the former at temperatures where it is crystalline is a molecular insulator, whereas the latter is a metal that takes on simple structures. On compression, however, the two come to share some structural and other similarities associated with the insulator-to-metal and metal-to-insulator transitions, respectively. To gain a deeper understanding of differences and parallels in the behaviors of compressed hydrogen and lithium, we performed an ab initio comparative study of these systems in selected identical structures. Both elements undergo a continuous pressure-induced s-p electronic transition, though this is at a much earlier stage of development for H. The valence charge density accumulates in interstitial regions in Li but not in H in structures examined over the same range of compression. Moreover, the valence charge density distributions or electron localization functions for the same arrangement of atoms mirror each other as one proceeds from one element to the other. Application of the virial theorem shows that the kinetic and potential energies jump across the first-order phase transitions in H and Li are opposite in sign because of non-local effects in the Li pseudopotential. Finally, the common tendency of compressed H and Li to adopt three-fold coordinated structures as found is explained by the fact that such structures are capable of yielding a profound pseudogap in the electronic densities of states at the Fermi level, thereby reducing the kinetic energy. These results have implications for the phase diagrams of these elements and also for the search for new structures with novel properties.