Magnetostriction-Driven Muon Localization in an Antiferromagnetic Oxide.

Magnetostriction results from the coupling between magnetic and elastic degrees of freedom. Though it is associated with a relatively small energy, we show that it plays an important role in determining the site of an implanted muon, so that the energetically favorable site can switch on crossing a magnetic phase transition. This surprising effect is demonstrated in the cubic rocksalt antiferromagnet MnO which undergoes a magnetostriction-driven rhombohedral distortion at the Néel temperature T_{N}=118 K. Above T_{N}, the muon becomes delocalized around a network of equivalent sites, but below T_{N} the distortion lifts the degeneracy between these equivalent sites. Our first-principles simulations based on Hubbard-corrected density-functional theory and molecular dynamics are consistent with the experimental data and help to resolve a long-standing puzzle regarding muon data on MnO, as well as having wider applicability to other magnetic oxides.

Observation of the most H2-dense filled ice under high pressure.

Hydrogen hydrates are among the basic constituents of our solar system's outer planets, some of their moons, as well Neptune-like exo-planets. The details of their high-pressure phases and their thermodynamic conditions of formation and stability are fundamental information for establishing the presence of hydrogen hydrates in the interior of those celestial bodies, for example, against the presence of the pure components (water ice and molecular hydrogen). Here, we report a synthesis path and experimental observation, by X-ray diffraction and Raman spectroscopy measurements, of the most H[Formula: see text]-dense phase of hydrogen hydrate so far reported, namely the compound 3 (or C[Formula: see text]). The detailed characterisation of this hydrogen-filled ice, based on the crystal structure of cubic ice I (ice I[Formula: see text]), is performed by comparing the experimental observations with first-principles calculations based on density functional theory and the stochastic self-consistent harmonic approximation. We observe that the extreme (up to 90 GPa and likely beyond) pressure stability of this hydrate phase is due to the close-packed geometry of the hydrogen molecules caged in the ice I[Formula: see text] skeleton.

First-Principles Thermodynamics of CsSnI3.

CsSnI3 is a promising ecofriendly solution for energy harvesting technologies. It exists at room temperature in either a black perovskite polymorph or a yellow 1D double-chain, which irreversibly deteriorates in the air. In this work, we unveil the relative thermodynamic stability between the two structures with a first-principles sampling of the CsSnI3 finite-temperature phase diagram, discovering how it is driven by anomalously large quantum and anharmonic ionic fluctuations. Thanks to a comprehensive treatment of anharmonicity, the simulations deliver a remarkable agreement with known experimental data for the transition temperatures of the orthorhombic, rhombohedral, and cubic perovskite structures and the thermal expansion coefficient. We disclose how the perovskite polymorphs are the ground state above 270 K and discover an abnormal decrease in heat capacity upon heating in the cubic black perovskite. Our results also significantly downplay the Cs+ rattling modes' contribution to mechanical instability. The remarkable agreement with experiments validates our methodology, which can be systematically applied to all metal halides.

Understanding anharmonic effects on hydrogen desorption characteristics of MgnH2n nanoclusters by ab initio trained deep neural network.

Magnesium hydride (MgH2) has been widely studied for effective hydrogen storage. However, its bulk desorption temperature (553 K) is deemed too high for practical applications. Besides doping, a strategy to decrease such reaction energy for releasing hydrogen is the use of MgH2-based nanoparticles (NPs). Here, we investigate first the thermodynamic properties of MgnH2n NPs (n < 10) from first-principles, in particular by assessing the anharmonic effects on the enthalpy, entropy and thermal expansion by means of the stochastic self consistent harmonic approximation (SSCHA). This method goes beyond previous approaches, typically based on molecular mechanics and the quasi-harmonic approximation, allowing the ab initio calculation of the fully-anharmonic free energy. We find an almost linear dependence on temperature of the interatomic bond lengths - with a relative variation of few percent over 300 K - alongside with a bond distance decrease of the Mg-H bonds. In order to increase the size of MgnH2n NPs toward experiments of hydrogen desorption we devise a computationally effective machine learning model trained to accurately determine the forces and total energies (i.e. the potential energy surfaces), integrating the latter with the SSCHA model to fully include the anharmonic effects. We find a significative decrease of the H-desorption temperature for sub-nanometric clusters MgnH2n with n ≤ 10, with a non-negligible, although little effect due to anharmonicities (up to 10%).

The microscopic origin of the anomalous isotopic properties of ice relies on the strong quantum anharmonic regime of atomic vibration.

Water ice is a unique material presenting intriguing physical properties, such as negative thermal expansion and anomalous volume isotope effect (VIE). They arise from the interplay between weak hydrogen bonds and nuclear quantum fluctuations, making theoretical calculations challenging. Here, we employ the stochastic self-consistent harmonic approximation to investigate how thermal and quantum fluctuations affect the physical properties of ice XI with ab initio accuracy. Regarding the anomalous VIE, our work reveals that quantum effects on hydrogen are so strong to be in a nonlinear regime: When progressively increasing the mass of hydrogen from protium to infinity (classical limit), the volume first expands and then contracts, with a maximum slightly above the mass of tritium. We observe an anharmonic renormalization of about 10% in the bending and stretching phonon frequencies probed in IR and Raman experiments. For the first time, we report an accurate comparison of the low-energy phonon dispersion with the experimental data, possible only thanks to high-level accuracy in the electronic correlation and nuclear quantum and thermal fluctuations, paving the way for the study of thermal transport in ice from first-principles and the simulation of ice under pressure.

Dominant Role of Quantum Anharmonicity in the Stability and Optical Properties of Infinite Linear Acetylenic Carbon Chains.

Carbyne, an infinite-length straight chain of carbon atoms, is supposed to undergo a second order phase transition from the metallic bond-symmetric cumulene (═C═C═)∞ toward the distorted insulating polyyne chain (-C≡C-)∞ displaying bond-length alternation. However, recent synthesis of ultra long carbon chains (∼6000 atoms, [Nat. Mater., 2016, 15, 634]) did not show any phase transition and detected only the polyyne phase, in agreement with previous experiments on capped finite carbon chains. Here, by performing first-principles calculations, we show that quantum-anharmonicity reduces the energy gain of the polyyne phase with respect to the cumulene one by 71%. The magnitude of the bond-length alternation increases by increasing temperature, in stark contrast with a second order phase transition, confining the cumulene-to-polyyne transition to extremely high and unphysical temperatures. Finally, we predict that a high temperature insulator-to-metal transition occurs in the polyyne phase confined in insulating nanotubes with sufficiently large dielectric constant due to a giant quantum-anharmonic bandgap renormalization.

TOLOMEO, a Novel Machine Learning Algorithm to Measure Information and Order in Correlated Networks and Predict Their State.

We present ToloMEo (TOpoLogical netwOrk Maximum Entropy Optimization), a program implemented in C and Python that exploits a maximum entropy algorithm to evaluate network topological information. ToloMEo can study any system defined on a connected network where nodes can assume N discrete values by approximating the system probability distribution with a Pottz Hamiltonian on a graph. The software computes entropy through a thermodynamic integration from the mean-field solution to the final distribution. The nature of the algorithm guarantees that the evaluated entropy is variational (i.e., it always provides an upper bound to the exact entropy). The program also performs machine learning, inferring the system's behavior providing the probability of unknown states of the network. These features make our method very general and applicable to a broad class of problems. Here, we focus on three different cases of study: (i) an agent-based model of a minimal ecosystem defined on a square lattice, where we show how topological entropy captures a crossover between hunting behaviors; (ii) an example of image processing, where starting from discretized pictures of cell populations we extract information about the ordering and interactions between cell types and reconstruct the most likely positions of cells when data are missing; and (iii) an application to recurrent neural networks, in which we measure the information stored in different realizations of the Hopfield model, extending our method to describe dynamical out-of-equilibrium processes.

The stochastic self-consistent harmonic approximation: calculating vibrational properties of materials with full quantum and anharmonic effects.

The efficient and accurate calculation of how ionic quantum and thermal fluctuations impact the free energy of a crystal, its atomic structure, and phonon spectrum is one of the main challenges of solid state physics, especially when strong anharmonicy invalidates any perturbative approach. To tackle this problem, we present the implementation on a modular Python code of the stochastic self-consistent harmonic approximation (SSCHA) method. This technique rigorously describes the full thermodynamics of crystals accounting for nuclear quantum and thermal anharmonic fluctuations. The approach requires the evaluation of the Born-Oppenheimer energy, as well as its derivatives with respect to ionic positions (forces) and cell parameters (stress tensor) in supercells, which can be provided, for instance, by first principles density-functional-theory codes. The method performs crystal geometry relaxation on the quantum free energy landscape, optimizing the free energy with respect to all degrees of freedom of the crystal structure. It can be used to determine the phase diagram of any crystal at finite temperature. It enables the calculation of phase boundaries for both first-order and second-order phase transitions from the Hessian of the free energy. Finally, the code can also compute the anharmonic phonon spectra, including the phonon linewidths, as well as phonon spectral functions. We review the theoretical framework of the SSCHA and its dynamical extension, making particular emphasis on the physical inter pretation of the variables present in the theory that can enlighten the comparison with any other anharmonic theory. A modular and flexible Python environment is used for the implementation, which allows for a clean interaction with other packages. We briefly present a toy-model calculation to illustrate the potential of the code. Several applications of the method in superconducting hydrides, charge-density-wave materials, and thermoelectric compounds are also reviewed.

Anomalous High-Temperature Superconductivity in YH6.

Pressure-stabilized hydrides are a new rapidly growing class of high-temperature superconductors, which is believed to be described within the conventional phonon-mediated mechanism of coupling. Here, the synthesis of one of the best-known high-TC superconductors-yttrium hexahydride I m 3 ¯ m -YH6 is reported, which displays a superconducting transition at ≈224 K at 166 GPa. The extrapolated upper critical magnetic field Bc2 (0) of YH6 is surprisingly high: 116-158 T, which is 2-2.5 times larger than the calculated value. A pronounced shift of TC in yttrium deuteride YD6 with the isotope coefficient 0.4 supports the phonon-assisted superconductivity. Current-voltage measurements show that the critical current IC and its density JC may exceed 1.75 A and 3500 A mm-2 at 4 K, respectively, which is higher than that of the commercial superconductors, such as NbTi and YBCO. The results of superconducting density functional theory (SCDFT) and anharmonic calculations, together with anomalously high critical magnetic field, suggest notable departures of the superconducting properties from the conventional Migdal-Eliashberg and Bardeen-Cooper-Schrieffer theories, and presence of an additional mechanism of superconductivity.

van der Waals driven anharmonic melting of the 3D charge density wave in VSe2.

Understanding of charge-density wave (CDW) phases is a main challenge in condensed matter due to their presence in high-Tc superconductors or transition metal dichalcogenides (TMDs). Among TMDs, the origin of the CDW in VSe2 remains highly debated. Here, by means of inelastic x-ray scattering and first-principles calculations, we show that the CDW transition is driven by the collapse at 110 K of an acoustic mode at qCDW = (2.25 0 0.7) r.l.u. The softening starts below 225 K and expands over a wide region of the Brillouin zone, identifying the electron-phonon interaction as the driving force of the CDW. This is supported by our calculations that determine a large momentum-dependence of the electron-phonon matrix-elements that peak at the CDW wave vector. Our first-principles anharmonic calculations reproduce the temperature dependence of the soft mode and the TCDW onset only when considering the out-of-plane van der Waals interactions, which reveal crucial for the melting of the CDW phase.

Weak Dimensionality Dependence and Dominant Role of Ionic Fluctuations in the Charge-Density-Wave Transition of NbSe_{2}.

Contradictory experiments have been reported about the dimensionality effect on the charge-density-wave transition in 2H NbSe_{2}. While scanning tunneling experiments on single layers grown by molecular beam epitaxy measure a charge-density-wave transition temperature in the monolayer similar to the bulk, around 33 K, Raman experiments on exfoliated samples observe a large enhancement of the transition temperature up to 145 K. By employing a nonperturbative approach to deal with anharmonicity, we calculate from first principles the temperature dependence of the phonon spectra both for bulk and monolayer. In both cases, the charge-density-wave transition temperature is estimated as the temperature at which the phonon energy of the mode driving the structural instability vanishes. The obtained transition temperature in the bulk is around 59 K, in rather good agreement with experiments, and it is just slightly increased in the single-layer limit to 73 K, showing the weak dependence of the transition on dimensionality. Environmental factors could motivate the disagreement between the transition temperatures reported by experiments. Our analysis also demonstrates the predominance of ionic fluctuations over electronic ones in the melting of the charge-density-wave order.

Anharmonicity and Doping Melt the Charge Density Wave in Single-Layer TiSe2.

Low-dimensional systems with a vanishing band gap and a large electron-hole interaction have been proposed to be unstable toward exciton formation. As the exciton binding energy increases in low dimension, conventional wisdom suggests that excitonic insulators should be more stable in 2D than in 3D. Here we study the effects of the electron-hole interaction and anharmonicity in single-layer TiSe2. We find that, contrary to the bulk case and to the generally accepted picture, in single-layer TiSe2, the electron-hole exchange interaction is much smaller in 2D than in 3D and it has weak effects on phonon spectra. By calculating anharmonic phonon spectra within the stochastic self-consistent harmonic approximation, we obtain TCDW ≈ 440 K for an isolated and undoped single layer and TCDW ≈ 364 K for an electron-doping n = 4.6 × 1013 cm-2, close to the experimental result of 200-280 K on supported samples. Our work demonstrates that anharmonicity and doping melt the charge density wave in single-layer TiSe2.

Quantum crystal structure in the 250-kelvin superconducting lanthanum hydride.

The discovery of superconductivity at 200 kelvin in the hydrogen sulfide system at high pressures1 demonstrated the potential of hydrogen-rich materials as high-temperature superconductors. Recent theoretical predictions of rare-earth hydrides with hydrogen cages2,3 and the subsequent synthesis of LaH10 with a superconducting critical temperature (Tc) of 250 kelvin4,5 have placed these materials on the verge of achieving the long-standing goal of room-temperature superconductivity. Electrical and X-ray diffraction measurements have revealed a weakly pressure-dependent Tc for LaH10 between 137 and 218 gigapascals in a structure that has a face-centred cubic arrangement of lanthanum atoms5. Here we show that quantum atomic fluctuations stabilize a highly symmetrical [Formula: see text] crystal structure over this pressure range. The structure is consistent with experimental findings and has a very large electron-phonon coupling constant of 3.5. Although ab initio classical calculations predict that this [Formula: see text] structure undergoes distortion at pressures below 230 gigapascals2,3, yielding a complex energy landscape, the inclusion of quantum effects suggests that it is the true ground-state structure. The agreement between the calculated and experimental Tc values further indicates that this phase is responsible for the superconductivity observed at 250 kelvin. The relevance of quantum fluctuations calls into question many of the crystal structure predictions that have been made for hydrides within a classical approach and that currently guide the experimental quest for room-temperature superconductivity6-8. Furthermore, we find that quantum effects are crucial for the stabilization of solids with high electron-phonon coupling constants that could otherwise be destabilized by the large electron-phonon interaction9, thus reducing the pressures required for their synthesis.

Quantum Enhancement of Charge Density Wave in NbS2 in the Two-Dimensional Limit.

At ambient pressure, bulk 2H-NbS2 displays no charge density wave instability, which is at odds with the isostructural and isoelectronic compounds 2H-NbSe2, 2H-TaS2, and 2H-TaSe2, and in disagreement with harmonic calculations. Contradictory experimental results have been reported in supported single layers, as 1H-NbS2 on Au(111) does not display a charge density wave, whereas 1H-NbS2 on 6H-SiC(0001) endures a 3 × 3 reconstruction. Here, by carrying out quantum anharmonic calculations from first-principles, we evaluate the temperature dependence of phonon spectra in NbS2 bulk and single layer as a function of pressure/strain. For bulk 2H-NbS2, we find excellent agreement with inelastic X-ray spectra and demonstrate the removal of charge ordering due to anharmonicity. In the two-dimensional limit, we find an enhanced tendency toward charge density wave order. Freestanding 1H-NbS2 undergoes a 3 × 3 reconstruction, in agreement with data on 6H-SiC(0001) supported samples. Moreover, as strains smaller than 0.5% in the lattice parameter are enough to completely remove the 3 × 3 superstructure, deposition of 1H-NbS2 on flexible substrates or a small charge transfer via field-effect could lead to devices with dynamical switching on/off of charge order.

Phonon Collapse and Second-Order Phase Transition in Thermoelectric SnSe.

Since 2014 the layered semiconductor SnSe in the high-temperature Cmcm phase is known to be the most efficient intrinsic thermoelectric material. Making use of first-principles calculations we show that its vibrational and thermal transport properties are determined by huge nonperturbative anharmonic effects. We show that the transition from the Cmcm phase to the low-symmetry Pnma is a second-order phase transition driven by the collapse of a zone border phonon, whose frequency vanishes at the transition temperature. Our calculations show that the spectral function of the in-plane vibrational modes are strongly anomalous with shoulders and double-peak structures. We calculate the lattice thermal conductivity obtaining good agreement with experiments only when nonperturbative anharmonic scattering is included. Our results suggest that the good thermoelectric efficiency of SnSe is strongly affected by the nonperturbative anharmonicity.

Manipulating Impulsive Stimulated Raman Spectroscopy with a Chirped Probe Pulse.

Photophysical and photochemical processes are often dominated by molecular vibrations in various electronic states. Dissecting the corresponding, often overlapping, spectroscopic signals from different electronic states is a challenge hampering their interpretation. Here we address impulsive stimulated Raman spectroscopy (ISRS), a powerful technique able to coherently stimulate and record Raman-active modes using broadband pulses. Using a quantum-mechanical treatment of the ISRS process, we show the mode-specific way the various spectral components of the broadband probe contribute to the signal generated at a given wavelength. We experimentally demonstrate how to manipulate the signal by varying the probe chirp and the phase-matching across the sample, thereby affecting the relative phase between the various contributions to the signal. These novel control knobs allow us to selectively enhance desired vibrational features and distinguish spectral components arising from different excited states.