Abstract
During a band-gap-tuned semimetal-to-semiconductor transition, Coulomb attraction between electrons and holes can cause spontaneously formed excitons near the zero-band-gap point, or the Lifshitz transition point. This has become an important route to realize bulk excitonic insulators – an insulating ground state distinct from single-particle band insulators. How this route manifests from weak to strong coupling is not clear. In this work, using angle-resolved photoemission spectroscopy (ARPES) and high-resolution synchrotron x-ray diffraction (XRD), we investigate the broken symmetry state across the semimetal-to-semiconductor transition in a leading bulk excitonic insulator candidate system Ta2Ni(Se,S)5. A broken symmetry phase is found to be continuously suppressed from the semimetal side to the semiconductor side, contradicting the anticipated maximal excitonic instability around the Lifshitz transition. Bolstered by first-principles and model calculations, we find strong interband electron-phonon coupling to play a crucial role in the enhanced symmetry breaking on the semimetal side of the phase diagram. Our results not only provide insight into the longstanding debate of the nature of intertwined orders in Ta2NiSe5, but also establish a basis for exploring band-gap-tuned structural and electronic instabilities in strongly coupled systems.
Similar content being viewed by others
Introduction
When the indirect band gap of a material is continuously reduced from positive to negative, band theory predicts a semiconductor-to-semimetal Lifshitz transition. However, such a transition may never be reached in the presence of electron-hole Coulomb attraction. It is postulated that poor screening in low carrier density semimetals and strong electron-hole binding in narrow gap semiconductors can both result in the spontaneous formation of excitons1,2. The condensation of excitons results in an insulating ground state, the so-called excitonic insulator (EI), which is separated from the normal state by a dome-shaped phase boundary peaked around the Lifshitz point (see Fig. 1a)3. This simple setup became a primary platform to search for condensed excitons and other exotic correlated phases, as well as potential solid-state BCS-BEC crossover phenomena4,5,6,7. For example, adding electronic anisotropy can further narrow the EI dome due to reduced Coulomb attraction8,9, and sufficiently heavy holes may cause Wigner crystallization10. Moreover, band-gap-tuned BCS-BEC crossover has been recently suggested in 2D moiré systems11. However, interband electron-phonon interaction can also lead to a gapped charge-density-wave state that occurs concurrently with the exciton condensation4. It remains unclear how the excitonic instability and the density-wave instability interact with each other, and whether it is possible at all to distinguish them. Since no super-transport (charge or heat) is associated with the excitonic insulator phase, and its analogy to BCS superconductivity is only formal10, an effective experimental method to distinguish the Coulomb and lattice channels is to identify a tuning parameter that gives differentiating predictions of the phase diagram.
Band-gap-controlled materials that undergo semimetal-to-semiconductor transitions provide an ideal platform to investigate this issue. However, the experimental realization of such a transition in bulk systems has been challenging. Divalent metals under hydrostatic pressure and uniaxial strain were first proposed as promising candidates12, but experimental evidence was indirect and scarce without clear ways to separate Coulomb and electron-phonon interaction effects. Recent experiments suggest a potential exciton condensation in 1T-TiSe213, where the band gap tunability is limited. Quasi-one-dimensional (quasi-1D) ternary chalcogenide Ta2NiSe5 has recently emerged as another leading EI candidate, exhibiting an orthorhombic-to-monoclinic structural transition at Ts = 329 K, concomitant with a putative exciton condensation14,15. The system consists of parallel quasi-1D Ta and Ni chains (see Fig. 1c), where long-lived electron-hole pairs are supposed to dwell on Ta 5\({d}_{{x}^{2}-{y}^{2}}\) and Ni 3dxy orbitals. On the one hand, the sizable ground-state single-particle gap16,17,18,19, minute structural distortion14,20, and a “dome-like” temperature-pressure phase diagram10,15 are all ostensibly consistent with the predictions in a band-gap-tuned EI. On the other hand, strong electron-phonon coupling effects are also found to coexist with correlation effects20,21,22,23,24,25,26,27,28,29,30,31,32. With the presence of both interactions, the band structure alone no longer sufficiently dictates the excitonic insulator phase region, and it is unclear how the low-temperature broken symmetry state would evolve across the Lifshitz transition point. In Ta2Ni(Se,S)5, this question can be addressed since the value of the direct band gap can be effectively tuned by isovalent sulfur substitution of selenium15,33. However, the phase diagram of Ta2Ni(Se,S)5 – especially where the broken symmetry phase ends and where the normal state crosses from a gapless to a gapped electronic structure with S-doping – has seen vastly contradicting results14,15,27, mainly due to the lack of direct probes of the single-particle band gap and the order parameter of the broken symmetry.
In this work, using angle-resolved photoemission spectroscopy (ARPES) and high-resolution synchrotron x-ray diffraction (XRD), we directly probe the electronic and lattice structure of Ta2Ni(Se,S)5 system and map out its phase diagram across the semimetal-to-semiconductor. A broken symmetry phase is found to be continuously suppressed from the semimetal side to the semiconductor side, contradicting the anticipated maximal excitonic instability around the Lifshitz transition. Bolstered by first-principles and model calculations, we reveal the crucial role of strong interband electron-phonon coupling behind the phase diagram.
Results
The general impact of S-substitution on the low-energy electronic structure can be revealed through the normal-state (defined here as the high-temperature symmetric structure phase) resistivity measured along the chain direction, which slowly increases until an abrupt upturn around 70% doping (see Fig. 1b). Such an upturn indicates a potential Lifshitz transition on the normal-state band structures. Hence, S-doped Ta2NiSe5 constitutes an ideal platform to study band-gap-tuned semimetal-to-semiconductor transition in an electron-phonon coupled correlated system. We first track the evolution of the broken symmetry phase boundary in Ta2Ni(Se,S)5with direct measurements of the structural order parameter. In bulk excitonic insulators, exciton condensation is always tied to a simultaneous lattice distortion10. Employing high-resolution single-crystal synchrotron XRD, Fig. 2a characterizes the temperature- and S-substitution-dependent evolution of the lattice order parameter β (also see Supplementary Note 1), which is defined as the angle between the lattice a and c axes shown in Fig. 1c. We find the structural phase transition remains second-order throughout the entire S-doping range. The ground-state order parameter β(T = 0) monotonically decreases with increasing S-doping (Fig. 2b), and the structural transition temperature only approaches zero near full S-substitution (Ta2NiS5) (bottom panel in Fig. 2b). This is notably different from the electronic phase diagram inferred from resistivity (see Fig. 2b and Supplementary Note 1), which instead suggests an abrupt collapse of the broken symmetry phase at an intermediate S-doping level15. Therefore, it is imperative to obtain a direct electronic structure view to determine the normal state semimetal-to-semiconductor transition point pL along the S-doping axis.
Previous studies have reported conflicting values of pL from resistivity, Raman measurements15,27. Here, we conduct high-resolution ARPES measurement, which is a direct probe of electronic structure, to reveal the evolution of band overlap/gap of Ta2Ni(Se,S)5 across the semimetal-semiconductor phase transition. The high statistics of the measurements enable electronic state restoration up to 4.5 kBT above EF (see Methods), and the orbital selectivity through polarized X-ray beam, allows the direct comparison with density functional theory (DFT) calculations, both of which help us precisely locate pL. Figure 3 shows the photoemission spectra in both the normal and the broken symmetry state (also see Supplementary Note 1). As the S-doping level increases, the Ta 5\({d}_{{x}^{2}-{y}^{2}}\) orbital component (blue) gradually withdraws from the valence band top, which becomes dominated by the Ni 3dxy orbital (red). Pseudogap states, the strong depletion of spectra intensity around Fermi level19,20, are also observed, which account for the high-temperature insulating behavior in these semimetallic samples (see ref. 20 for detailed description). Beyond pL ( ~ 70% S-doping), the normal-state conduction and valence bands are completely separated in energy (see Fig. 3a) and, simultaneously in the broken-symmetry-state, the Ta orbital component largely disappears from the valence band top (see Fig. 3b). Fitting the normal-state band dispersions (see Supplementary Note 2 for details) quantifies the band gap (or band overlap in the case of negative values), which changes sign across pL (Fig. 3c). In addition, revealed through the “M-shape” dispersion of the valence band top in the broken symmetry state, the band back-bending momentum kF also approaches Γ near pL (Fig. 3d). Both results are consistent with a concurrent normal state semimetal-to-semiconductor transition. Our DFT results further confirm such a Lifshitz transition near pL via the band-gap-sign flip and the kF reduction at the single-particle level, without involving ad hoc electron-hole Coulomb attraction (see Fig. 3c, d and Supplementary Note 3).
The ARPES measurement of the semimetal-to-semiconductor transition doping is consistent with that deduced from resistivity measurement (Figs. 1b and 2b). However, the structural phase transition temperature Ts and the structural order parameter (approximated as \({\beta }_{\max }-9{0}^{\circ }\)) are found to decrease monotonically with increasing S-doping through our XRD measurements (Fig. 2b). This poses a marked departure from the archetypal band-gap-tuned excitonic insulator phase diagram, where the excitonic instability maximizes around pL (Fig. 1a). To highlight this contrast, Fig. 4a shows the evolution of both the measured structural transition temperature Ts and the experimentally determined Ta atom displacement ∣X∣ in Ta2Ni(Se,S)5. Both ∣X∣ and Ts appear insensitive to the system crossing pL, which leads to considerations beyond weak electron-hole Coulomb interaction. One possibility is a heavy-hole Wigner crystal phase10; but the electron and hole effective masses are both near the free electron mass in Ta2Ni(Se,S)533. Alternatively, strong Coulomb interaction may push the phase boundary towards the semiconductor side34, but the Hartree component will undo this effect20 (see also Supplementary Note 4). Considering the fact that strong electron-phonon coupling effect is revealed in both equilibrium and ultrafast pump-probe experiments in Ta2Ni(Se,S)520,21,22,23,24,25,26,27,28, we hereby investigate an experimentally informed minimal model with strong electron-phonon coupling, shown in Fig. 1c and Eq. (1) of Methods. This model contains both electron-hole Coulomb interactions and an interband electron-phonon coupling, which can well describe both the spectral properties of Ta2NiSe5 below and above Ts20. To examine the impact of strong electron-phonon coupling on the band-gap-tuned transition as a matter-of-principle question, we first set the inter-band Coulomb interaction to zero. The role of this Coulomb interaction primarily lies in a Hartree shift between two bands and can be mapped to a correction of Eg (see detailed discussions in Supplementary Note 4). The average phonon displacement ∣xph∣ at T = 0 is employed to approximate the order parameters ∣X∣ measured in the experiment (Fig. 4a). With the quantized phonon modes, a zero-point fluctuation leads to the zero-coupling baseline at ∣xph∣ = 1.
Solved with exact diagonalization, ∣xph∣ for various Eg and coupling strengths are shown in Fig. 4b. In the weak-coupling case (e.g. g = 14 meV), the lattice order parameter displays a dome-like structure similar to the EI phase diagram derived based on Fermi-surface instability in Fig. 1a. With increasing coupling strengths, the dome gradually evolves into a monotonic crossover. This evolution exists not only for the many-body measure of fluctuations, but also for order parameters when symmetry breaking is allowed (see Supplementary Note 4). The monotonic crossover driven by strong, non-perturbative interactions reflects the irrelevance of Fermi-surface instability (top panel in Fig. 4c). The coupling strength in Ta2NiSe5, estimated from both ARPES spectral fitting and first-principles deformation potentials20, is 50 ~ 60 meV, placing it in the strong-coupling regime. While more complicated models with multiple bands may be necessary to quantitatively reproduce the order parameter evolution, strong electron-phonon coupling qualitatively describes the evolution of the broken-symmetry phase across pL in experiments. Moreover, strong coupling also expands the momentum range over which the lattice becomes unstable35, opening up the opportunity for a q = 0 transition when 2kF is small.
We demonstrate that strong electron-phonon interaction can greatly enhance the structural symmetry breaking in a semimetal, and the evolution of the broken-symmetry phase boundary can be indifferent to the semimetal-to-semiconductor transition. On the specific debate over the nature of the symmetry breaking in Ta2Ni(Se,S)5 systems, our results suggest a substantial deviation from the electron-hole Coulomb attraction-based EI phase (bottom panel in Fig. 4c). Instead, a strong electron-phonon-coupling induced structural transition can explain the observed symmetry breaking. This conclusion is consistent with recent experimental evidence of the crucial role of phonons in the equilibrium and nonequilibrium properties of Ta2NiSe522,23. Meanwhile, it has been suggested recently that additional Coulomb and acoustic-phonon effects can also cooperatively contribute to the transition27. On the high-temperature fluctuations, previous investigations have identified a normal-state spectral “pseudogap”, characterized by an anomalous spectral weight accumulation near the conduction band bottom. This phenomenon is interpreted as the signature of preformed excitons due to strong electronic correlation or electron coupling to lattice fluctuations20,36. We find that such a spectral anomaly ubiquitously exists in the normal state on the semimetal side of Ta2Ni(Se,S)5, but rapidly disappears on the semiconductor side (Fig. 3a). Intriguingly, the boundary of the pseudogap, or preformed excitons, is almost temperature-independent, similar to the abruptly terminated pseudogap in hole-doped superconducting cuprates37,38. Last but not least, the heavily S-substituted Ta2Ni(Se,S)5 system provides a promising platform to study both solid state BCS-BEC crossover phenomenon39 and critical fluctuations of low-energy phonons. Intriguingly, a potential competing order is observed in Fig. 2a for S-dopings above pL, similar to the situation in iron-based high Tc superconductors, where multiple orders are found competing with each other at low temperatures, including nematicity, spin-density wave, and superconductivity40,41.
Methods
Sample synthesis and resistivity measurements
Single crystals of Ta2Ni(Se,S)5 were grown by the chemical vapor transfer method with iodine (I2) as the transport agent. Starting materials composed of Ta powder (99.99%), Ni powder (99.99%), Se powder (99.99%), and S powder (99.99%) with a nominal molar ratio 2:1:5(1-x):5x (x is the nominal doping level of S) were fully ground and mixed together inside the glovebox. About 50 mg of iodine (I2) was added to the mixture of the starting powder. The mixture was then vacuumed, back-filled with 1/3 Argon, and sealed inside a quartz tube with an inner diameter of 8 mm, an outer diameter of 12 mm, and a length of about 120 mm. The sealed quartz tube was placed horizontally inside a muffle furnace during the growth. The hot end reaction temperature was set to 950 °C, and the cold end was left in the air with the temperature stabilized at 850 °C. Long and thin single crystals were harvested by quenching the furnace in the air after one week of reaction. Excess iodine was removed from the surfaces of the crystals with ethanol. Electrical resistivity measurements were carried out on a commercial PPMS (Quantum Design) by the four-probe method with the current applied along the a-axis of the Ta2Ni(Se,S)5 single crystals.
Angle-resolved photoemission spectroscopy (ARPES)
Synchrotron-based ARPES measurements were performed at the beamline BL5-2 of Stanford Synchrotron Radiation Laboratory (SSRL), SLAC, USA, and the BL 7.0.2 of the Advanced Light Source (ALS), USA. The samples were cleaved in situ and measured under the ultra-high vacuum below 3 × 10−11 Torr. Data was collected by the DA30L analyzer. The total energy and angle resolutions were 10 meV and 0.2°, respectively. High statistics measurements were done with a signal-to-noise ratio of 100 ~ 150 at the Fermi level on the energy distribution curve taken at the Fermi momentum. This enables electronic structure restoration up to 4.5 ~ 5 kBT after Fermi function division38,42. Orbital content selection is based on the dipole transition matrix element effect, following the experimental geometry described in ref. 20.
Single crystal x-ray scattering
Hard X-ray single-crystal diffraction is carried out at the energy of 44 keV at the beamline QM2 of the Cornell High Energy Synchrotron Source (CHESS). The needle-like sample is chosen with a typical lateral dimension of 100 microns, then mounted with GE Varnish on a rotating pin before being placed in the beam. A Pilatus 6M 2D area detector is used to collect the diffraction pattern with the sample rotated 365° around three different axes at a 0.1° step and 0.1s/frame data rate at each temperature. The full 3D intensity cube is stacked and indexed with the beamline software.
First-principles DFT calculation
Ab initio calculations are performed using the Quantum ESPRESSO package43,44. The crystal structure relaxations are performed using the r2SCAN functional45 with a semiempirical Grimme’s DFT-D2 van der Waals correction46. S-doped compounds are simulated by the virtual crystal approximation (VCA), where the pseudopotentials of Se and S atoms are linearly interpolated (i.e. mixed together) according to the chemical composition. A 30 × 30 × 15 k-mesh and a 100 Ry wavefunction energy cut-off were used. The electronic structure calculations are performed using the same r2SCAN functional with the relaxed structures.
Two-band model for many-body simulations
We consider the two-band model with inter-band e-ph couplings to describe the physics in the Ta2Ni(Se,S)5 systems, whose Hamiltonian reads as47:
where \({c}_{k\sigma }^{{{{\dagger}}} }\) (ckσ) creates (annihilates) an electron at the conduction band (primarily Ta 5d) for momentum k and spin σ, with dispersion given by \({\varepsilon }_{k}^{c}\); and the \({f}_{k\sigma }^{{{{\dagger}}} }\) (fkσ) creates (annihilates) an electron at the valence band (primarily Ni 3d), with dispersion given by \({\varepsilon }_{k}^{v}\) (Fig. 1c). The \({n}_{i\sigma }^{c}\) and \({n}_{i\sigma }^{f}\) are the density operators for the conduction and valence band, respectively. We employ the band structures fitted from the Ta2NiSe5 experiments, as reported in ref. 20, which determine the valence and conduction band dispersion for a negative band gap Eg = − 0.3 eV. Due to the difficulty of fitting the full dispersions in the gapped Ta2Ni(Se,S)5 materials, we simplify the model using a rigid band separation controlled by Eg. That being said, the conduction and valence band structures for Ta2Ni(Se,S)5 read as
Here, Eg = 0 indicates the Lifshitz transition where \({\varepsilon }_{k=0}^{c}={\varepsilon }_{k=0}^{v}\), while Eg = 0.3 eV reflects the situation of Ta2NiS5.
To approximate the lattice distortion in a finite-size simulation, we employ the average displacement of phonons
where ρ is the density matrix of the e-ph system and \({\hat{x}}_{{{{{{{{\rm{ph}}}}}}}}}=({a}_{0}^{{{{\dagger}}} }+{a}_{0})/N\) is the uniform phonon displacement operator (In a finite-size system, the exact ground state does not spontaneously break symmetries, so that \(\langle {\hat{x}}_{{{{{{{{\rm{ph}}}}}}}}}\rangle \equiv 0\)).
Data availability
All data generated in this study have been deposited in the Figshare database with open access at the https://doi.org/10.6084/m9.figshare.24435481.
References
Mott, N. F. The transition to the metallic state. Philos. Mag. 6, 287 (1961).
Knox, R. Introduction to exciton physics, Collective Excitations in Solids, 183 (1983).
Keldysh, L. & Kopaev, Y. V. Possible instability of semimetallic state toward Coulomb interaction. Sov. Phys. Solid State 6, 2219 (1965).
Halperin, B. & Rice, T. The excitonic state at the semiconductor-semimetal transition. Solid State Phys. 21, 115 (1968).
Fedders, P. A. & Martin, P. C. Itinerant antiferromagnetism. Phys. Rev. 143, 245 (1966).
Arrott, A., Werner, S. & Kendrick, H. First-order magnetic phase change in chromium at 38.5 °C. Phys. Rev. Lett. 14, 1022 (1965).
Bronold, F. X. & Fehske, H. Possibility of an excitonic insulator at the semiconductor-semimetal transition. Phys. Rev. B 74, 165107 (2006).
Zittartz, J. Anisotropy effects in the excitonic insulator. Phys. Rev. 162, 752 (1967).
Kozlov, A. & Maksimov, L. The metal-dielectric divalent crystal phase transition. Sov. Phys. JETP 21, 790 (1965).
Halperin, B. & Rice, T. Possible anomalies at a semimetal-semiconductor transition. Rev. Mod. Phys. 40, 755 (1968).
Kim, H. et al. Spectroscopic signatures of strong correlations and unconventional superconductivity in twisted trilayer graphene, arXiv preprint arXiv:2109.12127 (2021).
Jérome, D., Rice, T. & Kohn, W. Excitonic insulator. Phys. Rev. 158, 462 (1967).
Kogar, A. et al. Signatures of exciton condensation in a transition metal dichalcogenide. Science 358, 1314 (2017).
Di Salvo, F. et al. Physical and structural properties of the new layered compounds Ta2NiS5 and Ta2NiSe5. J. Less-common Met. 116, 51 (1986).
Lu, Y. et al. Zero-gap semiconductor to excitonic insulator transition in Ta2NiSe5. Nat. Commun. 8, 14408 (2017).
Lee, J. et al. Strong interband interaction in the excitonic insulator phase of Ta2NiSe5. Phys. Rev. B 99, 075408 (2019).
Wakisaka, Y. et al. Excitonic insulator state in Ta2NiSe5 probed by photoemission spectroscopy. Phys. Rev. Lett. 103, 026402 (2009).
Wakisaka, Y. et al. Photoemission spectroscopy of Ta2NiSe5. J. Supercond. Nov. Magn. 25, 1231 (2012).
Seki, K. et al. Excitonic bose-einstein condensation in Ta2NiSe5 above room temperature. Phys. Rev. B 90, 155116 (2014).
Chen, C. et al. Role of electron-phonon coupling in excitonic insulator candidate Ta2NiSe5, Phys. Rev. Res. 5, 043089 (2023).
Mazza, G. et al. Nature of symmetry breaking at the excitonic insulator transition: Ta2NiSe5. Phys. Rev. Lett. 124, 197601 (2020).
Baldini, E. et al. The spontaneous symmetry breaking in Ta2NiSe5 is structural in nature. Proc. Natl. Acad. Sci. 120, e2221688120 (2023).
Watson, M. D. et al. Band hybridization at the semimetal-semiconductor transition of Ta2NiSe5 enabled by mirror-symmetry breaking. Phys. Rev. Res. 2, 013236 (2020).
Bretscher, H. M. et al. Ultrafast melting and recovery of collective order in the excitonic insulator Ta2NiSe5. Nat. Commun. 12, 1699 (2021).
Yan, J. et al. Strong electron-phonon coupling in the excitonic insulator Ta2NiSe5. Inorg. Chem. 58, 9036 (2019).
Volkov, P. A. et al. Critical charge fluctuations and emergent coherence in a strongly correlated excitonic insulator. npj Quantum Mater. 6, 52 (2021).
Volkov, P. A. et al. Failed excitonic quantum phase transition in Ta2Ni(Se1−xSx)5. Phys. Rev. B 104, L241103 (2021).
Kim, K. et al. Direct observation of excitonic instability in Ta2NiSe5. Nat. Commun. 12, 1969 (2021).
Larkin, T. et al. Giant exciton Fano resonance in quasi-one-dimensional Ta2NiSe5. Phys. Rev. B 95, 195144 (2017).
Larkin, T. et al. Infrared phonon spectra of quasi-one-dimensional Ta2NiSe5 and Ta2NiS5. Phys. Rev. B 98, 125113 (2018).
Werdehausen, D. et al. Coherent order parameter oscillations in the ground state of the excitonic insulator Ta2NiSe5. Sci. Adv. 4, eaap8652 (2018).
Bretscher, H. M. et al. Imaging the coherent propagation of collective modes in the excitonic insulator Ta2NiSe5 at room temperature. Sci. Adv. 7, eabd6147 (2021).
Saha, T. et al. Electronic band structure in pristine and sulfur-doped Ta2NiSe5. Phys. Rev. B 106, 075148 (2022).
Sugimoto, K., Nishimoto, S., Kaneko, T. & Ohta, Y. Strong coupling nature of the excitonic insulator state in Ta2NiSe5. Phys. Rev. Lett. 120, 247602 (2018).
Zhu, X., Cao, Y., Zhang, J., Plummer, E. & Guo, J. Classification of charge density waves based on their nature. Proc. Natl. Acad. Sci. 112, 2367 (2015).
Fukutani, K. et al. Detecting photoelectrons from spontaneously formed excitons. Nat. Phys. 17, 1024 (2021).
He, Y. et al. Rapid change of superconductivity and electron-phonon coupling through critical doping in Bi-2212. Science 362, 62 (2018).
Chen, S.-D. et al. Incoherent strange metal sharply bounded by a critical doping in Bi2212. Science 366, 1099 (2019).
Phan, V.-N., Becker, K. W. & Fehske, H. Spectral signatures of the BCS-BEC crossover in the excitonic insulator phase of the extended Falicov-Kimball model. Phys. Rev. B 81, 205117 (2010).
Yi, M. et al. Dynamic competition between spin-density wave order and superconductivity in underdoped Ba(1−x)KxFe2As2. Nat. Commun. 5, 3711 (2014).
Fernandes, R., Chubukov, A. & Schmalian, J. What drives nematic order in iron-based superconductors? Nat. Phys. 10, 97 (2014).
He, Y. et al. Superconducting fluctuations in overdoped Bi2Sr2CaCu2O8+δ. Phys. Rev. X 11, 031068 (2021).
Giannozzi, P. et al. QUANTUM ESPRESSO: A Modular and open-source software project for quantum simulations of materials. J. Condens. Matter Phys. 21, 395502 (2009).
Giannozzi, P. et al. Advanced capabilities for materials modelling with Quantum ESPRESSO. J. Condens. Matter Phys. 29, 465901 (2017).
Furness, J. W., Kaplan, A. D., Ning, J., Perdew, J. P. & Sun, J. Accurate and numerically efficient r2SCAN meta-generalized gradient approximation. J. Phys. Chem. Lett. 11, 8208 (2020).
Grimme, S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J. Comput. Chem. 27, 1787 (2006).
Kaneko, T., Toriyama, T., Konishi, T. & Ohta, Y. Orthorhombic-to-monoclinic phase transition of Ta2NiSe5 induced by the Bose-Einstein condensation of excitons. Phys. Rev. B 87, 035121 (2013).
Acknowledgements
The authors wish to thank A. F. Kemper, E. Demler, D. Y. Qiu, P. J. Guo, P. A. Volkov, L. I. Glazman for helpful discussions. Use of the Stanford Synchrotron Radiation Light Source, SLAC National Accelerator Laboratory, is supported by the US Department of Energy, Office of Science, Office of Basic Energy Sciences under Contract No. DE-AC02-76SF00515. M. H. and D. L. acknowledge the support of the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Material Sciences and Engineering, under Contract No. DE-AC02-76SF00515. This research used resources of the Advanced Light Source, a US DOE Office of Science User Facility under Contract No. DE-AC02-05CH11231. Research conducted at the Center for High-Energy X-ray Science (CHEXS) is supported by the National Science Foundation (BIO, ENG and MPS Directorates) under award DMR-1829070. Work at Lawrence Berkeley National Laboratory was funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Materials Sciences and Engineering Division under Contract No. DE-AC02-05-CH11231 within the Quantum Materials Program (KC2202) which provided the numerical simulations and within the Theory of Materials Program (KC2301) which provided the DFT calculations. The quantum many-body simulations (S.D. and Y.W.) are supported initially by U.S. Department of Energy, Office of Science, Basic Energy Sciences, under Early Career Award No. DE-SC0022874 and are completed with the support under award No. DE-SC0024524. The many-body and DFT simulations were performed on the Frontera and Stampede2 computing system, respectively, at the Texas Advanced Computing Center. The work at Yale University is partially supported by the National Science Foundation (NSF) under DMR-2239171.
Author information
Authors and Affiliations
Contributions
Y.H. and Y.W. conceived the project. C.C. carried out ARPES measurements and analysis with the assistance of Z.K., S.W., Y.H., Y.C.. Single crystals were synthesized and characterized by X.C. with guidance from R.B.. XRD measurements and analysis were performed by X.C. and Y.H. with the help of Z.K. and J.R.. DFT calculations were performed by W.T, Z.L. with guidance from S.L.. Model calculations were performed by S.D and Y.W.. K.S. and E.N carried out preliminary STM characterization. J.R., M.H., D.L., C.J., A.B., and E.R. maintain the XRD and ARPES beamlines. C.C., W.Y, W.T. and Y.H. wrote the manuscript with the input from all authors. All authors contributed to the scientific planning and discussion.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks Alberto de la Torre, and the other, anonymous, reviewers for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Chen, C., Tang, W., Chen, X. et al. Anomalous excitonic phase diagram in band-gap-tuned Ta2Ni(Se,S)5. Nat Commun 14, 7512 (2023). https://doi.org/10.1038/s41467-023-43365-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-023-43365-1
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.