Anharmonicity Reveals the Tunability of the Charge Density Wave Orders in Monolayer VSe2

Anharmonicity Reveals the Tunability of the Charge Density Wave Orders in Monolayer VSe2 Adolfo Otero Fumega,*Josu Diego, Victor Pardo, Santiago Blanco-Canosa, and Ion Errea* Cite This: Nano Lett. 2023, 23, 1794−1800 Read Online ACCESS Metrics & More Article Recommendations * sı Supporting Information ABSTRACT: VSe2is a layered compound that has attracted great attention due to its proximity to a ferromagnetic state that is quenched by its charge density wave (CDW) phase. In the monolayer limit, unrelated experiments have reported different CDW orders with different transition temperatures, making this monolayer very controversial. Here we perform first-principles nonperturbative anharmonic phonon calculations in monolayer VSe2in order to estimate the CDW order and the corresponding transition temperature. They reveal that monolayer VSe2develops two independent charge density wave orders that compete as a function of strain. Variations of only 1.5% in the lattice parameter are enough to stabilize one order or the other. Moreover, we analyze the impact of external Lennard-Jones interactions, showing that these can act together with anharmonicity to suppress the CDW orders. Our results solve previous experimental contradictions, highlighting the high tunability and substrate dependency of the CDW orders of monolayer VSe2. KEYWORDS: anharmonic effects, charge density wave, competing orders, 2D materials, van der Waals interactions, strain Two-dimensional (2D) materials are an ideal platform to artificially engineer heterostructures with new functionalities due to the weak van der Waals bonding between layers. 1 Monolayers hosting symmetry-broken phases, such as superconductivity, 2,3 magnetism, 4−8 ferroelectricity, 9,10 charge density waves (CDWs), 11,12 or multiferroicity, 13,14 represent the most interesting building blocks to design novel phases of matter. One of the main challenges in the task of engineering novel functional materials with broken-symmetry monolayers is to overcome the restrictions imposed by the reduced dimensionality, 15,16 which may prevent the formation of these phases, and the competition between ordered phases due to the subtle interplay of different interactions. 17,18 For instance, CDW phases have been reported to destroy 19,20 or promote 21 2D ferromagnetism. VSe2is a paradigmatic example of this as, despite some early claims, 22 it is now clear both experimentally and theoretically that the CDW order quenches the emergence of itinerant ferromagnetism. 19,20,23−28 In its bulk form, VSe2 develops a commensurate 4 ×4×3 CDW phase below 110 K. 29 The CDW phase opens pseudogaps at the Fermi level impeding the emergence of ferromagnetism. 19 Inelastic X-ray scattering experiments and nonperturbative anharmonic phonon calculations have proven that the CDW transition is driven by the collapse of a low-energy acoustic mode and that the electron−phonon coupling is the origin of the instability, 30 as suggested as well by other quantitative models. 31 These anharmonic calculations have shown that van der Waals interactions are essential to melt the CDW and obtain a charge density wave temperature (TCDW) in agreement with experiments. This suggests that the CDW in the monolayer may also be characterized by similar phonon softening effects but with limited influence of van der Waals interactions. The main problem in the monolayer of VSe2is that the CDW is not fully understood yet, as unrelated experiments have reported distinct CDW orders with different transition temperatures. A nonmonotonic evolution of TCDW as a function of the number of layers has been reported in refs 32−34 but retaining an in-plane 4 ×4 modulation. A metastable phase with modulation has also been identified for the few-layer case. 35 In the purely 2D limit different CDW orders with nonequivalent modulations have been found. A 4 ×4 order was observed in VSe2films grown on bilayer graphene on top of SiC and on highly oriented pyrolytic graphite (HOPG) with a TCDW of K and a lattice parameter of a= 3.31 ±0.05 Å. 27 On the contrary, a modulation has been observed in VSe2samples grown on several substrates by molecular beam epitaxy by different groups, with a consistent TCDW = 220 K. 20,23 Some Received: November 21, 2022 Revised: February 22, 2023 Published: February 24, 2023 Letterpubs.acs.org/NanoLett © 2023 The Authors. Published by American Chemical Society 1794 https://doi.org/10.1021/acs.nanolett.2c04584 Nano Lett. 2023, 23, 1794−1800 Downloaded via UNIV DEL PAIS VASCO on April 26, 2023 at 15:55:58 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles. other orders have also been reported: a combination of and with a TCDW ∼135 K 25,26 and a 4 ×1 modulation with TCDW ∼350 K. 26,36 These experimental contradictions point to the presence of different competing CDW orders, which can lead to different low-temperature phases depending on the substrate. 26,36 By calculating the harmonic phonons of the VSe2monolayer within density functional theory (DFT), theoretical studies have also described the competition of different CDW orders and how strain can influence the ground state. 37 Harmonic phonon calculations, however, cannot explain that above TCDW the 1T phase is the ground state. In the presence of competing orders, only calculations considering anharmonicity can disentangle what is the CDW order and the transition temperature, as it has already been shown in different transition metal dichalcogenides (TMDs). 30,38−41 Therefore, in order to unveil the intrinsic CDW orders of monolayer VSe2 and how they are affected by external fields, a DFT study including anharmonicity is required. In this work, we present a theoretical analysis of the CDW orders arising in monolayer VSe2using nonperturbative anharmonic phonon calculations based on the stochastic selfconsistent harmonic approximation (SSCHA). 42−45 This formalism has been crucial to understand and characterize the CDWs in several TMDs 30,38−41 as it overcomes the limitations of the harmonic analysis, allowing to determine the dependence of the CDW order as a function of temperature. We demonstrate the emergence and competition of two intrinsic CDW orders in monolayer VSe2as for a lattice parameter of a= 3.35 Å a order dominates with TCDW = 217 K, while for slightly smaller a= 3.30 Å the 4 ×4 order prevails with TCDW = 223 K. Moreover, the nonperturbative anharmonic procedure allows us to demonstrate that the CDW can be suppressed by the inclusion of LennardJones energy terms, which might appear naturally or may be artificially induced by the interplay between the monolayer and a particular substrate. We start our analysis performing harmonic phonon calculations on monolayer VSe2. The normal state (NS) unit cell is shown in Figure 1a. The value of the experimentally reported lattice parameter a= 3.31 ±0.05 Å is in rather good agreement with the theoretical one of 3.35 Å obtained at the Perdew−Burke−Ernzerhof 46 level without considering the zero-point motion. 23,27 Therefore, in this study we perform calculations for two lattice parameters a= 3.35 Å and a= 3.30 Å, which provide a good representation of the experimental range. Density functional perturbation theory (DFPT) 47 is used to compute the harmonic phonon band structure for both lattice parameters (see the Supporting Information for a detailed description of the calculations). Both harmonic phonon bands (see Figure 1b) show two dominant instabilities at and . These are associated with the two intrinsic CDW orders of monolayer VSe2, with modulations shown in Figure 1c,d, that lower its Born− Oppenheimer energy. The instability at q1is associated with a supercell, while the one at q2leads to a 4 ×4 modulation. Both softened phonon modes have an out-ofplane component in the displacement vectors, as it is also the case of the CDW instability in the bulk form of this compound. In spite of providing the two intrinsic CDW orders, harmonic calculations do not suffice to predict which of these CDW orders is the dominant one or the associated transition temperature for each lattice parameter. In fact, the small change in the lattice parameter does not impact the weight of the instabilities. Our nonperturbative anharmonic calculations based on a free energy formalism within the SSCHA can give the answer to these questions (see the Supporting Information for a detailed description of the SSCHA method and the technical aspects of these calculations). Figure 2 shows the temperature evolution of the phonon band structure obtained with the SSCHA method for the two lattice parameters a= 3.35 Å (in red in the top panels of Figure 2) and a= 3.30 Å (in blue in the lower panels of Figure 2). At high enough temperature (T= 250 K) the 1T NS phase (Figure 1a) is dynamically stable for both lattice parameters as shown in Figures 2a,d, showing that anharmonicity melts the CDW phase as it happens in other TMDs. 30,38−41 By decreasing the temperature, the phonon modes associated Figure 1. (a) Normal-state structure of monolayer VSe2with lattice parameter a. V (Se) atoms are depicted in gray (green). (b) Harmonic phonon band structures of monolayer VSe2as a function of the lattice parameter, left (right) panel for a= 3.35 Å (a= 3.30 Å). Two dominant instabilities at and can be identified. (c,d) Intrinsic CDW orders with and 4 ×4 modulations associated with the instabilities at q1and q2can be identified in the harmonic phonon band structures. The displacement vectors associated with each CDW order are plotted as brown arrows. Planes perpendicular to the z-direction for V (Se) were plotted in gray (green) for a better characterization of the displacement vectors. Nano Letters pubs.acs.org/NanoLett Letter https://doi.org/10.1021/acs.nanolett.2c04584 Nano Lett. 2023, 23, 1794−1800 1795 with the CDW instabilities at q1and q2soften. In particular, for a= 3.35 Å at 200 K (Figure 2b) we can observe that the mode at becomes unstable, even if the one at q2remains stable. This result indicates that for this lattice parameter the CDW order dominates. However, for a= 3.30 Å at 200 K (Figure 2e) the phonon mode at is unstable, but not at q1. Therefore, for the smaller lattice parameter, the 4 ×4 CDW order is the dominant one. At low enough temperatures (Figures 2c,f) both q-vectors show unstable modes. However, note that this situation is not indicating that at low temperatures both CDW orders coexist, although it is a clear signature that the anharmonic free energy landscape becomes more complex. Once one of the CDW orders becomes stable when the temperature is decreased, the system collapses to it, and the analysis in terms of the anharmonic phonons of the NS phase is no longer useful to describe the evolution of each of the CDW phases at low temperatures. Nevertheless, the anharmonic phonons at low temperature shown in Figure 2c,f confirm that both q1and q2 are the intrinsic CDW orders of VSe2that can be accessed through a transition from the NS phase. To analyze in more detail the competition between the two CDW orders as a function of the lattice parameter, Figure 3a shows the temperature evolution of the frequency of the phonon mode that softens at q1and q2. For the larger lattice parameter, a= 3.35 Å, the frequency at q1becomes negative (imaginary) at higher temperature than at q2, and hence the CDW order is the dominant (left panel in Figure 3a in red). The opposite behavior is observed for the small lattice parameter a= 3.30 Å. The frequency at q2becomes negative at higher temperature than that at q1, and hence the 4 ×4 CDW order is dominant (right panel in Figure 3a in blue). From Figure 3a we can obtain the transition temperature for each lattice parameter: for a= 3.35 Å the order emerges at TCDW = 217 K, while for a= 3.30 Å the 4 ×4 order arises at TCDW = 223 K. Importantly, our anharmonic calculations including the zero point energy 44 predict an associated in-plane pressure of 0.7 GPa for a= 3.35 Å and 1.3 GPa for a= 3.30 Å. Considering that its in-plane pressure is lower, these results point out that the intrinsic CDW order in monolayer VSe2is with a TCDW = 217 K, which is in perfect agreement with the experiments on refs 20 and 23 and that the 4 ×4 order, which is the in-plane projection of the bulk 4 ×4×3 CDW order, appears only under strain. Our results provide an explanation for the different CDW orders observed for small variations (∼1.5%) of the lattice parameter. 23,27 Note that, eventually, other modulations could appear in monolayer VSe2 as experimentally reported. 25,26,36 However, our results show that the and 4 ×4 modulations are the intrinsic CDW orders and point out that those different modulations are a consequence of the interplay between the highly dynamically unstable NS of VSe2monolayer and the particular substrate. Having established the competition between the two intrinsic CDW orders of monolayer VSe2as a function of the lattice parameter, we study now the origin of these CDW orders. In order to do so, we use DFPT to compute both the nesting function η(q) (Figure 3b), which is given by (1) and the phonon line width associated with the electron− phonon interaction (see Figure 3c). (2) In eqs 1 and 2ϵnkis the energy of band nwith wavenumber k,ϵFthe Fermi energy, and Nis the number of kpoints in the sum over the first Brillouin zone (1BZ). The nesting function Figure 2. Nonperturbative anharmonic calculations of the NS of monolayer VSe2. (a−c) For lattice parameter a= 3.35 Å and temperatures of 250, 200, 100 K respectively. (d−f) For lattice parameter a= 3.30 Å and temperatures 250, 200, 150 K. The phonon melting occurs first at the q1(q2) point for a= 3.35 Å (a= 3.30 Å) as shown in panels (b) and (e). Nano Letters pubs.acs.org/NanoLett Letter https://doi.org/10.1021/acs.nanolett.2c04584 Nano Lett. 2023, 23, 1794−1800 1796 peaks for qindicates that nested regions of the Fermi surface connect and, thus, reveal if the instability emerges from a purely electronic instability. The equation for γμ(q) is very similar to the nesting function, but the value is weighted by the mode μand momentum qdependent electron−phonon matrix elements and, thus, reveals if the instability emerges from electron−phonon interactions. The electron−phonon line width is independent of the phonon frequency wμ(q) as the electron−phonon matrix elements scale as wμ(q)−1/2. It is worth noting that both quantities need to be computed in order to establish the origin of the CDW orders. For 2D systems like monolayer VSe2, despite the existence of nesting conditions at the qpoints associated with the CDW orders, a purely electronic picture does not suffice to produce a CDW order and electron−phonon interactions play a key role to drive the transition. 48 Therefore, the direct comparison between the nesting function and the electron−phonon line width allows us to establish which is the main driving force of the CDW orders in monolayer VSe2. We can see in Figure 3b that, for both lattice parameters the nesting function does not show any strong peak at the CDW vectors despite the existence of small shoulders near q1and q2. In a different way, Figure 3c shows that the phonon line width coming from the electron−phonon interaction abruptly peaks at both q1and q2for both lattice parameters, meaning that in all cases the enhancement comes from the mode and momentum dependence of the electron−phonon matrix elements. Therefore, the two intrinsic CDW orders developed by monoloyer VSe2are driven by the electron−phonon coupling, in agreement with the theoretical predictions for 2D systems. 48 Besides, the electron−phonon interaction also plays a key role in the CDW transition in 3D systems, as in bulk 1T-VSe2, in which case the presence of nesting is symbolic. 30 The analysis about the stability of the different CDW orders as a function of strain was performed with a nonlocal van der Waals density exchange-correlation functional. 49 The reason for this is that this functional allows us to properly describe both bulk and monolayer limits of VSe2, oppositely to the widely used GGA-PBE functional (see the Supporting Information for a more detailed description of the election of the exchange correlation functional to study the CDW orders of VSe2). In the latter case, a huge overestimation of TCDW occurs in bulk due to the lack of van der Waals interactions that cause a melting of the CDW phase. 30 Motivated by these results in the bulk, we explore here the effect that external van der Waals interactions may cause in the CDW orders of monolayer VSe2. These van der Waals interactions might naturally appear by proximity effect between the analyzed monolayer and other layers, such as the substrate or in van der Waals heterostructures. We will introduce these external interactions in a phenomenological way. This approach allows us to derive general conclusions for the pure effect of van der Waals interactions on CDW orders, i.e., factoring out effects that could appear in particular substrates, such as charge transfer. 40 First, to illustrate the effect of external van der Waals forces on the CDW, we make use of a simple one-dimensional double-well fourth order potential (3) where Aand Bare the coefficients of the different powers, ris the position of the atoms, and r0is the equilibrium atomic position in the high temperature phase. The CDW occurs when the free energy calculated with this potential is lower at r −r0≠0 than at r−r0= 0. Obviously the lower and wider the minimum of the well the more probable it is to find a brokensymmetry CDW order. We can now add to the potential V(r) a ELJ(r) Lennard-Jones energetic contribution to mimic the role played by external van der Waals interactions (4) where Cis the coefficient that controls the strength of the Lennard-Jones interactions, rcis a cutoff radius that prevents a divergence at r=r0, and D=C/rcis simply a constant that fixes the potential to be equal to 0 at r0. The effect of the LennardJones interactions on V(r) can be seen in Figure 4a for different values of C. We can observe that an increase in the strength of the Lennard-Jones interactions makes the potential shallower. Therefore, Lennard-Jones interactions tend to quench the low-temperature CDW phase and promote the high-temperature symmetric phase. Figure 3. (a) Temperature evolution of the frequencies of the softened modes at q1and q2as obtained in the SSCHA calculations. (b) Nesting function η(q). It does not clearly peak at q1and q2. (c) Phonon line width γμ(q) given by the electron−phonon interaction. Sizable peaks appear at the CDW q-vectors. All of the results for a= 3.35 Å (a= 3.30 Å) are shown in red (blue) in the left (right) panels. Nano Letters pubs.acs.org/NanoLett Letter https://doi.org/10.1021/acs.nanolett.2c04584 Nano Lett. 2023, 23, 1794−1800 1797 We can confirm this simple picture in the particular case of monolayer VSe2by including an energy term like the one shown in eq 4 through the Grimme’s semiempirical approach in our SSCHA calculations on top of the PBE functional. 50 Figure 4c shows the harmonic phonon band structure including energy contributions from Lennard-Jones interactions. (These calculations are for a= 3.35 Å, but qualitatively the results hold for any other lattice parameter. The strength of the Lennard-Jones interactions was set to the one considered by default by the Grimme’s semiempirical correction implemented in the QUANTUM ESPRESSO package 51,52 ) The instability at q1and q2is reduced compared to the harmonic bands without the Lennard-Jones contribution (see the Supporting Information for a better visualization of this effect at the harmonic level.). In this situation, anharmonic effects are able to suppress both CDW orders by stabilizing the softened phonons even at 0 K, as shown in Figure 4(c). This effect can be also analyzed in Figure 4(b), where the evolution of the frequencies of the softened modes at q1and q2as a function of temperature is shown. The frequencies remain stable at any temperature. Therefore, this demonstrates that the combination of Lennard-Jones interactions and anharmonicity can destroy the CDW orders, stabilizing the NS phase at low temperatures. Note that here we have considered a strength of the Lennard-Jones interactions that quenches both CDW orders. However, this strength could be modulated by the parameter Cin eq 4, providing simply a decrease of TCDW, but not a suppression of the CDW order, as reported for bulk. 30 This effect may also be related with the enhancement of the CDW order in the 2D limit, 34 in which this kind of interactions decreases. Therefore, not only strain but also Lennard-Jones interactions can explain the huge variability of transition temperatures and CDW orders reported in experiments where VSe2is grown on different substrates. 23,25−27 Finally, note that monolayer VSe2has attracted great attention due to its proximity to an itinerant ferromagnetic state, which is suppressed by the presence of CDW orders. 19,20 Our analysis suggests that a ferromagnetic state in monolayer VSe2may be possible if its intrinsic CDW orders are suppressed by external Lennard-Jones interactions. This tuning could be implemented by substrate engineering or by artificial design of van der Waals heterostructures. In particular, combining compounds that display CDW orders with ferroelectric materials, which provide strong dipolar interactions, might allow an electric control of CDW phases and the emergence of other competing orders. This mechanism to tune or destroy CDW orders is generic and could be extended to other similar systems, offering a novel platform to engineering new functional materials. In conclusion, to solve previous experimental contradictions found for the CDW orders of monolayer VSe2, in this work we analyze the CDW orders of monolayer VSe2using nonperturbative anharmonic phonon calculations that allow us to determine the CDW orders of this system and their corresponding transition temperatures. We analyze the effect of strain and external van der Waals interactions, since these two parameters are intrinsic to any experiment. We demonstrate the competition between two intrinsic CDW orders as a function of the lattice parameter. Variations of 1.5% in the lattice parameter are enough to drive the system from the to the 4 ×4 order. Transition temperatures on the order of 220 K have been found for both CDW orders, in very good agreement with experiments. We show that external Lennard-Jones interactions tend to weaken or even suppress the CDW orders. These results together help to understand the great variability of CDW orders and associated TCDW’s found in the experiments of monolayer VSe2. Moreover, they pave the way to tune CDW orders occurring in van der Waals materials, thus promoting competing orders that might arise in these systems. ■ASSOCIATED CONTENT * sı Supporting Information Detailed explanation of the computational methods, the election of the exchange correlation functional and the effect of external Lennard-Jones interactions at the harmonic level. The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.nanolett.2c04584. Detailed explanation of the computational methods, the election of the exchange correlation functional and the effect of external Lennard-Jones interactions at the harmonic level (PDF) ■AUTHOR INFORMATION Corresponding Authors Adolfo Otero Fumega −Department of Applied Physics, Aalto University, 02150 Espoo, Finland; orcid.org/00000002-3385-6409; Email: [email protected] Ion Errea −Fisika Aplikatua Saila, Gipuzkoako Ingeniaritza Eskola, University of the Basque Country (UPV/EHU), 20018 San Sebastián, Spain; Centro de Física de Materiales (CSIC-UPV/EHU), 20018 San Sebastián, Spain; Donostia Figure 4. (a) V(r) potential as a function of the atomic position r described by eqs 3 and 4in arbitrary units (a.u.) for different Cvalues. Particular values of A=−8, B=−A/2, and rc= 0.1 are considered. (b) Temperature evolution of the soften modes’ frequencies at q1and q2when anharmonicity and Lennard-Jones interactions are included. (c) Harmonic and anharmonic phonon band structures at T= 0 K for monolayer VSe2including Lennard-Jones contributions. Nano Letters pubs.acs.org/NanoLett Letter https://doi.org/10.1021/acs.nanolett.2c04584 Nano Lett. 2023, 23, 1794−1800 1798 International Physics Center (DIPC), 20018 San Sebastián, Spain; Email: [email protected] Authors Josu Diego −Fisika Aplikatua Saila, Gipuzkoako Ingeniaritza Eskola, University of the Basque Country (UPV/EHU), 20018 San Sebastián, Spain; Centro de Física de Materiales (CSIC-UPV/EHU), 20018 San Sebastián, Spain; orcid.org/0000-0002-8659-2144 Victor Pardo −Departamento de Física Aplicada, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain; Instituto de Materiais iMATUS, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain; orcid.org/0000-0002-4713-3519 Santiago Blanco-Canosa −Donostia International Physics Center (DIPC), 20018 San Sebastián, Spain; IKERBASQUE, Basque Foundation for Science, 48013 Bilbao, Spain Complete contact information is available at: https://pubs.acs.org/10.1021/acs.nanolett.2c04584 Notes The authors declare no competing financial interest. ■ACKNOWLEDGMENTS We acknowledge the computational resources provided by the CESGA and the Aalto Science-IT project. A.O.F. acknowledges the financial support received through the Academy of Finland Project No. 349696. J.D. thanks the Department of Education of the Basque Government for a predoctoral fellowship (Grant No. PRE-2020-1-0220). We thank the Ministry of Science and Education of Spain for financial support through the projects PGC2018-101334-A-C22, GC2018-101334-B-C21, PID2021-122609NB-C22. I.E. acknowledges funding from the Department of Education, Universities and Research of the Eusko Jaurlaritza, and the University of the Basque Country UPV/EHU (Grant No. IT1527-22). ■REFERENCES (1) Geim, A. K.; Grigorieva, I. V. Van der Waals heterostructures. Nature 2013,499, 419−425. (2) Ugeda, M. M.; Bradley, A. J.; Zhang, Y.; Onishi, S.; Chen, Y.; Ruan, W.; Ojeda-Aristizabal, C.; Ryu, H.; Edmonds, M. T.; Tsai, H.- Z.; et al. Characterization of collective ground states in single-layer NbSe2. Nat. Phys. 2016,12, 92−97. (3) Vano, V.; Ganguli, S. C.; Amini, M.; Yan, L.; Khosravian, M.; Chen, G.; Kezilebieke, S.; Lado, J. L.; Liljeroth, P. Evidence of nodal fwave superconductivity in monolayer 1H-TaS2with hidden order fluctuations. arXiv e-prints 2021,DOI: 10.48550/arXiv.2112.07316. (4) Lee, J.-U.; Lee, S.; Ryoo, J. H.; Kang, S.; Kim, T. Y.; Kim, P.; Park, C.-H.; Park, J.-G.; Cheong, H. Ising-Type Magnetic Ordering in Atomically Thin FePS3. Nano Lett. 2016,16, 7433−7438. (5) Huang, B.; Clark, G.; Navarro-Moratalla, E.; Klein, D. R.; Cheng, R.; Seyler, K. L.; Zhong, D.; Schmidgall, E.; McGuire, M. A.; Cobden, D. H.; Yao, W.; Xiao, D.; Jarillo-Herrero, P.; Xu, X. Layer-dependent ferromagnetism in a van der Waals crystal down to the monolayer limit. Nature 2017,546, 270−273. (6) Gong, C.; Li, L.; Li, Z.; Ji, H.; Stern, A.; Xia, Y.; Cao, T.; Bao, W.; Wang, C.; Wang, Y.; Qiu, Z. Q.; Cava, R. J.; Louie, S. G.; Xia, J.; Zhang, X. Discovery of intrinsic ferromagnetism in two-dimensional van der Waals crystals. Nature 2017,546, 265−269. (7) Fei, Z.; Huang, B.; Malinowski, P.; Wang, W.; Song, T.; Sanchez, J.; Yao, W.; Xiao, D.; Zhu, X.; May, A. F.; Wu, W.; Cobden, D. H.; Chu, J.-H.; Xu, X. Two-dimensional itinerant ferromagnetism in atomically thin Fe3GeTe2. Nat. Mater. 2018,17, 778−782. (8) Zhang, Z.; Shang, J.; Jiang, C.; Rasmita, A.; Gao, W.; Yu, T. Direct Photoluminescence Probing of Ferromagnetism in Monolayer Two-Dimensional CrBr3. Nano Lett. 2019,19, 3138−3142. (9) Cui, C.; et al. Intercorrelated In-Plane and Out-of-Plane Ferroelectricity in Ultrathin Two-Dimensional Layered Semiconductor In2Se3. Nano Lett. 2018,18, 1253−1258. (10) Yuan, S.; Luo, X.; Chan, H. L.; Xiao, C.; Dai, Y.; Xie, M.; Hao, J. Room-temperature ferroelectricity in MoTe2 down to the atomic monolayer limit. Nat. Commun. 2019,10, 1775. (11) Wang, Y.; Ren, J.; Li, J.; Wang, Y.; Peng, H.; Yu, P.; Duan, W.; Zhou, S. Evidence of charge density wave with anisotropic gap in a monolayer VTe2film. Phys. Rev. B 2019,100, 241404. (12) Chen, P.; Pai, W. W.; Chan, Y.-H.; Takayama, A.; Xu, C.-Z.; Karn, A.; Hasegawa, S.; Chou, M. Y.; Mo, S.-K.; Fedorov, A.-V.; Chiang, T.-C. Emergence of charge density waves and a pseudogap in single-layer TiTe2. Nat. Commun. 2017,8, 516. (13) Song, Q.; Occhialini, C. A.; Ergecen, E.; Ilyas, B.; Amoroso, D.; Barone, P.; Kapeghian, J.; Watanabe, K.; Taniguchi, T.; Botana, A. S.; Picozzi, S.; Gedik, N.; Comin, R. Evidence for a single-layer van der Waals multiferroic. Nature 2022,602, 601−605. (14) Fumega, A. O.; Lado, J. L. Microscopic origin of multiferroic order in monolayer NiI2.2D Materials 2022,9, 025010. (15) Hohenberg, P. C. Existence of Long-Range Order in One and Two Dimensions. Phys. Rev. 1967,158, 383−386. (16) Mermin, N. D.; Wagner, H. Absence of Ferromagnetism or Antiferromagnetism in Oneor Two-Dimensional Isotropic Heisenberg Models. Phys. Rev. Lett. 1966,17, 1133−1136. (17) Du, L.; Hasan, T.; Castellanos-Gomez, A.; Liu, G.-B.; Yao, Y.; Lau, C. N.; Sun, Z. Engineering symmetry breaking in 2D layered materials. Nature Reviews Physics 2021,3, 193−206. (18) Li, W.; Qian, X.; Li, J. Phase transitions in 2D materials. Nature Reviews Materials 2021,6, 829−846. (19) Fumega, A. O.; Gobbi, M.; Dreher, P.; Wan, W.; GonzálezOrellana, C.; Pena-Díaz, M.; Rogero, C.; Herrero-Martín, J.; Gargiani, P.; Ilyn, M.; Ugeda, M. M.; Pardo, V.; Blanco-Canosa, S. Absence of Ferromagnetism in VSe2 Caused by Its Charge Density Wave Phase. J. Phys. Chem. C 2019,123, 27802−27810. (20) Coelho, P. M.; Nguyen Cong, K.; Bonilla, M.; Kolekar, S.; Phan, M.-H.; Avila, J.; Asensio, M. C.; Oleynik, I. I.; Batzill, M. Charge Density Wave State Suppresses Ferromagnetic Ordering in VSe2Monolayers. J. Phys. Chem. C 2019,123, 14089−14096. (21) Otero Fumega, A.; Phillips, J.; Pardo, V. Controlled TwoDimensional Ferromagnetism in 1T−CrTe2: The Role of Charge Density Wave and Strain. J. Phys. Chem. C 2020,124, 21047−21053. (22) Bonilla, M.; Kolekar, S.; Ma, Y.; Diaz, H. C.; Kalappattil, V.; Das, R.; Eggers, T.; Gutierrez, H. R.; Phan, M.-H.; Batzill, M. Strong room-temperature ferromagnetism in VSe2 monolayers on van der Waals substrates. Nat. Nanotechnol. 2018,13, 289−293. (23) Chen, P.; Pai, W. W.; Chan, Y.-H.; Madhavan, V.; Chou, M. Y.; Mo, S.-K.; Fedorov, A.-V.; Chiang, T.-C. Unique Gap Structure and Symmetry of the Charge Density Wave in Single-Layer VSe2.Phys. Rev. Lett. 2018,121, 196402. (24) Lei, B.-H.; Singh, D. J. Identification of a low-energy metastable 1T-type phase for monolayer VSe2.Phys. Rev. B 2021,104, 125430. (25) Biswas, D.; et al. Ultrafast Triggering of Insulator−Metal Transition in Two-Dimensional VSe2. Nano Lett. 2021,21, 1968− 1975. (26) Duvjir, G.; et al. Emergence of a Metal−Insulator Transition and High-Temperature Charge-Density Waves in VSe2 at the Monolayer Limit. Nano Lett. 2018,18, 5432−5438. (27) Feng, J.; et al. Electronic Structure and Enhanced ChargeDensity Wave Order of Monolayer VSe2. Nano Lett. 2018,18, 4493− 4499. (28) Vinai, G.; Bigi, C.; Rajan, A.; Watson, M. D.; Lee, T.-L.; Mazzola, F.; Modesti, S.; Barua, S.; Ciomaga Hatnean, M.; Balakrishnan, G.; King, P. D. C.; Torelli, P.; Rossi, G.; Panaccione, Nano Letters pubs.acs.org/NanoLett Letter https://doi.org/10.1021/acs.nanolett.2c04584 Nano Lett. 2023, 23, 1794−1800 1799 G. Proximity-induced ferromagnetism and chemical reactivity in fewlayer VSe2heterostructures. Phys. Rev. B 2020,101, 035404. (29) Eaglesham, D. J.; Withers, R. L.; Bird, D. M. Charge-densitywave transitions in 1T-VSe2. Journal of Physics C: Solid State Physics 1986,19, 359−367. (30) Diego, J.; Said, A. H.; Mahatha, S. K.; Bianco, R.; Monacelli, L.; Calandra, M.; Mauri, F.; Rossnagel, K.; Errea, I.; Blanco-Canosa, S. van der Waals driven anharmonic melting of the 3D charge density wave in VSe2. Nat. Commun. 2021,12, 598. (31) Henke, J.; Flicker, F.; Laverock, J.; van Wezel, J. Charge order from structured coupling in VSe2.SciPost Phys. 2020,9, 56. (32) Xu, K.; Chen, P.; Li, X.; Wu, C.; Guo, Y.; Zhao, J.; Wu, X.; Xie, Y. Ultrathin Nanosheets of Vanadium Diselenide: A Metallic TwoDimensional Material with Ferromagnetic Charge-Density-Wave Behavior. Angew. Chem., Int. Ed. 2013,52, 10477−10481. (33) Yang, J.; Wang, W.; Liu, Y.; Du, H.; Ning, W.; Zheng, G.; Jin, C.; Han, Y.; Wang, N.; Yang, Z.; Tian, M.; Zhang, Y. Thickness dependence of the charge-density-wave transition temperature in VSe2. Appl. Phys. Lett. 2014,105, 063109. (34) Pásztor, A.; Scarfato, A.; Barreteau, C.; Giannini, E.; Renner, C. Dimensional crossover of the charge density wave transition in thin exfoliated VSe 2. 2D Materials 2017,4, 041005. (35) Zhang, D.; Ha, J.; Baek, H.; Chan, Y.-H.; Natterer, F. D.; Myers, A. F.; Schumacher, J. D.; Cullen, W. G.; Davydov, A. V.; Kuk, Y.; Chou, M. Y.; Zhitenev, N. B.; Stroscio, J. A. Strain engineering a 4a× √3a charge-density-wave phase in transition-metal dichalcogenide 1T −VSe2.Phys. Rev. Materials 2017,1, 024005. (36) Duvjir, G.; Choi, B. K.; Ly, T. T.; Lam, N. H.; Jang, K.; Dung, D. D.; Chang, Y. J.; Kim, J. Multiple charge density wave phases of monolayer VSe2 manifested by graphene substrates. Nanotechnology 2021,32, 364002. (37) Si, J. G.; Lu, W. J.; Wu, H. Y.; Lv, H. Y.; Liang, X.; Li, Q. J.; Sun, Y. P. Origin of the multiple charge density wave order in 1T− VSe2.Phys. Rev. B 2020,101, 235405. (38) Bianco, R.; Errea, I.; Monacelli, L.; Calandra, M.; Mauri, F. Quantum Enhancement of Charge Density Wave in NbS2 in the Two-Dimensional Limit. Nano Lett. 2019,19, 3098−3103. (39) Bianco, R.; Monacelli, L.; Calandra, M.; Mauri, F.; Errea, I. Weak Dimensionality Dependence and Dominant Role of Ionic Fluctuations in the Charge-Density-Wave Transition of NbSe2.Phys. Rev. Lett. 2020,125, 106101. (40) Zhou, J. S.; Monacelli, L.; Bianco, R.; Errea, I.; Mauri, F.; Calandra, M. Anharmonicity and Doping Melt the Charge Density Wave in Single-Layer TiSe2. Nano Lett. 2020,20, 4809−4815. (41) Sky Zhou, J.; Bianco, R.; Monacelli, L.; Errea, I.; Mauri, F.; Calandra, M. Theory of the thickness dependence of the charge density wave transition in 1 T-TiTe2. 2D Materials 2020,7, 045032. (42) Errea, I.; Calandra, M.; Mauri, F. Anharmonic free energies and phonon dispersions from the stochastic self-consistent harmonic approximation: Application to platinum and palladium hydrides. Phys. Rev. B 2014,89, 064302. (43) Bianco, R.; Errea, I.; Paulatto, L.; Calandra, M.; Mauri, F. Second-order structural phase transitions, free energy curvature, and temperature-dependent anharmonic phonons in the self-consistent harmonic approximation: Theory and stochastic implementation. Phys. Rev. B 2017,96, 014111. (44) Monacelli, L.; Errea, I.; Calandra, M.; Mauri, F. Pressure and stress tensor of complex anharmonic crystals within the stochastic self-consistent harmonic approximation. Phys. Rev. B 2018,98, 024106. (45) Monacelli, L.; Bianco, R.; Cherubini, M.; Calandra, M.; Errea, I.; Mauri, F. The stochastic self-consistent harmonic approximation: calculating vibrational properties of materials with full quantum and anharmonic effects. J. Phys.: Condens. Matter 2021,33, 363001. (46) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996,77, 3865−3868. (47) Baroni, S.; de Gironcoli, S.; Dal Corso, A.; Giannozzi, P. Phonons and related crystal properties from density-functional perturbation theory. Rev. Mod. Phys. 2001,73, 515−562. (48) Johannes, M. D.; Mazin, I. I. Fermi surface nesting and the origin of charge density waves in metals. Phys. Rev. B 2008,77, 165135. (49) Thonhauser, T.; Cooper, V. R.; Li, S.; Puzder, A.; Hyldgaard, P.; Langreth, D. C. Van der Waals density functional: Self-consistent potential and the nature of the van der Waals bond. Phys. Rev. B 2007, 76, 125112. (50) Grimme, S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J. Comput. Chem. 2006,27, 1787−1799. (51) Giannozzi, P.; et al. QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials. J. Phys.: Condens. Matter 2009,21, 395502. (52) Giannozzi, P.; et al. Advanced capabilities for materials modelling with Q uantum ESPRESSO. J. Phys.: Condens. Matter 2017,29, 465901. Nano Letters pubs.acs.org/NanoLett Letter https://doi.org/10.1021/acs.nanolett.2c04584 Nano Lett. 2023, 23, 1794−1800 1800 Recommended by ACS Anisotropic Phonon Scattering and Thermal Transport Property Induced by the Liquid-like Behavior of AgCrSe2 Chen Wang and Yue Chen APRIL 17, 2023 NANO LETTERS READ Evidences of Topological Surface States in the Nodal-Line Semimetal SnTaS2 Nanoflakes Wenshuai Gao, Mingliang Tian, et al. FEBRUARY 20, 2023 ACS NANO READ Theoretical Study of Intercalation Effects: Graphene and hBN Layers in Metal and Monolayer Black Phosphorus Contacts Na Shen. APRIL 18, 2023 ACS APPLIED ELECTRONIC MATERIALS READ Pixel Super-Resolved Lensless on-Chip Sensor with Scattering Multiplexing Xuyang Chang, Liheng Bian, et al. JANUARY 31, 2023 ACS PHOTONICS READ Get More Suggestions >