---
abstract: 'Vertical integration of two-dimensional (2D) van der Waals (vdW) materials with different quantum ground states is predicted to lead to novel electronic properties that are not found in the constituent layers. Here, we present the direct synthesis of superconductor-magnet hybrid heterostructures by combining superconducting niobium diselenide (NbSe$_2$) with the monolayer (ML) vanadium diselenide (VSe$_2$). More significantly, the in-situ growth in ultra-high vacuum (UHV) allows to produce a clean and an atomically sharp interfaces. Combining different characterization techniques and density-functional theory (DFT) calculations, we investigate the electronic and magnetic properties of VSe$_2$ on NbSe$_2$. Low temperature scanning tunneling microscopy (STM) measurements demonstrate a reduction of the superconducting gap on VSe$_2$ layer. This together with the lack of charge density wave signatures indicates magnetization of the sheet, but not of a conventional itinerant ferromagnet.'
author:
- Shawulienu Kezilebieke
- Md Nurul Huda
- Paul Dreher
- Ilkka Manninen
- Yifan Zhou
- Jani Sainio
- Rhodri Mansell
- 'Miguel M. Ugeda'
- Sebastiaan van Dijken
- 'Hannu-Pekka Komsa'
- Peter Liljeroth
bibliography:
- 'vse2.bib'
title: 'Electronic and Magnetic Characterization of Epitaxial VSe$_2$ Monolayers on Superconducting NbSe$_2$'
---

There has been a surge of interest in designer materials that would realize electronic responses not found in naturally occurring materials. There are many routes towards this goal and they are all being explored vigorously: e.g., artificial atomic lattices[@Gomes2012DesignerGraphene; @Slot2017ExperimentalLattice; @Drost2017; @Khajetoorians2019designer; @Yan2019engineered], atomically precise graphene nanoribbons[@Talirz2016On-SurfaceNanoribbons; @Groning2018EngineeringNanoribbons; @Rizzo2018TopologicalNanoribbons], and controlled heterostructures of two-dimensional materials [@Geim2013; @Novoselov2016_review; @Huang2017_CrI3; @Gong2017_Cr2Ge2Te6; @Bonilla2018; @Cao2018UnconventionalSuperlattices; @Gibertini2019]. The designer concept is well illustrated in systems combining magnetism with superconductivity to realize topological superconductivity [@Zhang2011_RMP; @Mourik2012science; @Nadj-Perge2014; @Sato2017_review; @Zhou2019]. Individual magnetic impurity atoms give rise to so-called Yu-Shiba-Rusinov states [@Heinrich2018_review], which can be coupled in extended structures to give rise to bands (inside the superconducting gap). Eventually, the system can be driven into a topological phase in the presence of certain spin textures or spin-orbit coupling [@Nadj-Perge2014; @Ojanen2014_prb; @Kim2018; @Roentynen2015; @Li2016; @Rachel2017]. Topological superconductors are a distinct form of matter that is predicted to host boundary Majorana fermions. Experimental realization of Majorana fermions is exciting in its own right, but this is compounded by the proposal that systems with non-abelian statistics can serve as the basis for topological quantum computation[@Nayak2008; @Alicea2012; @Wilczek2009]. Experimentally, these systems have been realized in one-dimensional chains of magnetic adatoms on the surface of s-wave superconductors[@Nadj-Perge2014; @Kim2018; @Ruby2015] and this has been extended to two-dimensional systems by using magnetic metal islands [@Menard2017_NatComm; @palacio]. However, these types of systems can be sensitive to disorder and interface engineering through, *e.g.*, the use of an atomically thin separation layer, might be required [@palacio].

Issues with interface inhomogeneities can potentially be avoided in van der Waals (vdW) heterostructures, where the different layers interact only through vdW forces [@Geim2013; @Novoselov2016_review]. Layered materials that remain magnetic down to the monolayer (ML) limit have been recently demonstrated [@Huang2017_CrI3; @Gong2017_Cr2Ge2Te6]. CrI$_3$ and Cr$_2$Ge$_2$Te$_6$ have very desirable magnetic properties; however, they are not readily amenable to in-situ growth in ultra-high vacuum (UHV) leading to clean edges and interfaces (through, *e.g.*, molecular-beam epitaxy (MBE)). Very recently, monolayer magnetism was suggested in the transition metal dichalcogenide (TMD) vanadium diselenide (VSe$_2$), which can be readily grown using MBE on various layered materials [@Bonilla2018]. Later reports have questioned the existence of magnetism in VSe$_2$ as no magnetic signal was detected in X-ray magnetic circular dichroism (XMCD) experiments [@Feng2018_NanoLett; @Wong2019_frustrated]. Angle-resolved photoemission spectroscopy revealed an enhanced charge-density wave (CDW) transition at a higher temperature than in the bulk, and it was suggested that the presence of CDW driven by Fermi-surface nesting removes the usual mechanism for achieving a magnetic ground state [@Chen2018_PRL; @Duvjir2018_NanoLett; @Feng2018_NanoLett; @Coelho_2019jpcc]. While there still is no consensus on the nature of the possible magnetic ground state of VSe$_2$, it is clear that there is a delicate balance between different competing interacting states and phases in monolayer transition-metal dichalcogenides, which may also depend on the nature of the substrate.[@Feng2018_NanoLett; @Coelho_2019jpcc; @Wong2019_frustrated; @Sugawara2019_PRB; @Kim2019_arXiv; @Fumega2019; @Yu2019chemically]

Combining 2D magnetic and superconducting TMDs would create a promising material platform for investigating the coexistence of superconductivity, magnetism and the resulting emergent quantum phases of matter. The inherent lack of surface bonding sites due to the layered nature of these materials prevents chemical bonding between the layers and results in a better control of the interfaces. We report growth of single layer vanadium diselenide (VSe$_2$) by molecular beam epitaxy on superconducting niobium diselenide (NbSe$_2$) substrate and study the magnetic and electrical properties of these heterostructures. MBE growth under UHV conditions facilitates the formation of clean edges and interfaces and we characterize the electronic structure of the resulting hybrid materials down to the atomic scale using low-temperature scanning tunneling microscopy (STM) and spectroscopy (STS). Our results give further experimental information on the magnetic properties of VSe$_2$ and demonstrate a clean and controllable platform for creating superconducting-magnetic hybrid TMD materials with great potential of integrating TMDs into future electronics devices.

Results and discussion
======================

![Growth of VSe$_2$ on NbSe$_2$. (a) Large scale STM image of submonolayer VSe$_2$ on NbSe$_2$. (b) Line-profile along the blue line shown in panel a. (c,d) Atomically resolved images on VSe$_2$ (c) and NbSe$_2$ (d). (e) Line-profiles along the lines in panels c and d (VSe$_2$ (black), NbSe$_2$ (red)) showing the atomic periodicities and the charge-density wave modulation on the NbSe$_2$ substrate. (f) Computed structure of VSe$_2$ on NbSe$_2$. []{data-label="fig1"}](Fig1.jpg){width="95.00000%"}

VSe$_2$ was grown on NbSe$_2$ by MBE and results are illustrated in Fig. \[fig1\] (see Methods for details). Briefly, vanadium was evaporated under excess flux of selenium onto a NbSe$_2$ crystal cleaved in-situ in UHV and held at $T = 520-540$ K during the growth. The samples were characterized *in-situ* by STM and X-ray photoelectron spectroscopy (XPS). In addition, after capping the films with a thick Se layer, the samples were characterized *ex-situ* by temperature dependent magnetization measurements. Fig. \[fig1\]a shows STM characterization of sub-monolayer VSe$_2$ films on NbSe$_2$ substrate. VSe$_2$ grows atomically smooth, large uniform ML islands. Higher coverages result in the formation of a second layer. The profile along the blue line in Fig. \[fig1\]a shows that the apparent height of the VSe$_2$ film is 6.5 Å(Fig. \[fig1\]b) consistent with the unit cell height [@Bayard1976]. Atomically resolved STM images of the VSe$_2$ monolayer and NbSe$_2$ crystal surface are shown in Figs. \[fig1\]c and \[fig1\]d, respectively. While NbSe$_2$ shows the well-known $3\times3$ charge-density wave modulation in the atomic contrast, we do not detect a charge density wave on VSe$_2$ (even at temperature of $T=4.2$ K). This is in contrast to reports on HOPG and bilayer graphene substrates[@Bonilla2018; @Feng2018_NanoLett; @Chen2018_PRL]. The lattice constants can be measured from the atomically resolved images, as depicted in Fig. \[fig1\]e. This yields values of 3.5 Å and 3.4 Å  for VSe$_2$ and NbSe$_2$, respectively. These values match well with previous experimental results[@Bayard1976; @Duvjir2018_NanoLett; @Liu2018].

XPS was used to study the chemical composition of the as-grown VSe$_2$ films on NbSe$_2$. Characteristic peaks of V, Se and Nb are found in the XPS spectra (see Supporting Information (SI) Fig. S1). The binding energies of the V 2p$_{3/2}$ peak at 513.7 eV and the Se 3d$_{5/2}$ peak at 53.4 eV are similar to those previously observed for VSe$_2$ on HOPG[@Liu2018; @Bonilla2018]. The Nb 3d$_{5/2}$ peak is found at 203.5 eV which is typical for NbSe$_2$. Both VSe$_2$ and NbSe$_2$ have similar selenium binding energies and thus they cannot be resolved from the Se 3d spectrum [@Wang2017]. The V:Se:Nb stoichiometry estimated from the peak areas was roughly 1:5:2 (for a 0.6 ML VSe$_2$ film). No other elements, such as possible magnetic impurities, were detected above the detection limit of $\sim$1 atomic percent.

We have complemented the experiments by density functional theory (DFT) calculations (see Methods and SI for details). Fig. \[fig1\]f shows the fully relaxed geometry of VSe$_2$/NbSe$_2$ heterostructure from the side and top views. The energetically most favourable stacking has the lower layer Se atoms of VSe$_2$ on top of hollow site of the NbSe$_2$ (2.94 Å from Se in NbSe$_2$ to Se in VSe$_2$) and V on top of Nb, with a distance of 6.16 Å from Nb to V.

![(a) Typical long-range experimental d$I$/d$V$ spectra on a ML VSe$_2$ and on the NbSe$_2$ substrate. (b) Simulated STS obtained by integrating LDOS at a constant height on top of the NbSe$_2$ surface. The solid and dashed lines show different values of the energy broadening. (c) Calculated density of states of bulk NbSe$_2$. Density of state plots show both the total DOS and the DOS projected to the metal and Se atoms. []{data-label="fig2"}](Fig2.png){width="50.00000%"}

We have probed the electronic structure of single-layer VSe$_2$/NbSe$_2$ heterostructure by both scanning tunnelling spectroscopy (STS) and DFT. Fig. \[fig2\]a shows typical d$I$/d$V$ spectra taken on the ML VSe$_2$ and on the bulk NbSe$_2$ substrate over a large bias range. We will first focus on the NbSe$_2$ response. At positive bias (empty states) region, the most pronounced features on bulk NbSe$_2$ are the broad resonances at $\sim0.3$ V and $\sim1.8$ V, while at negative bias the d$I$/d$V$ signal is broad and rising. The measured d$I$/d$V$ spectrum on NbSe$_2$ is in agreement with earlier STS studies on bulk NbSe$_2$[@Arguello2014]. The features in our d$I$/d$V$ spectroscopy also match the simulated spectra on a 3-layer slab of NbSe$_2$ (Fig. \[fig2\]b) and can be compared with the bulk density of states (Fig. \[fig2\]c). The first resonance at positive bias arises from the Nb-derived band, while the broad feature at negative bias overlaps with the mostly selenium derived bands below $E_\mathrm{F}$. On VSe$_2$, at positive bias, there are pronounced features close to the Fermi level and also at 0.9 V and 1.45 V. At negative bias, we observe a peak at -0.5 V, and a shallow feature at $\sim -1$ V. These VSe$_2$ features are discussed in more detail below.

![(a) Calculated band structure of VSe$_2$ (non-magnetic state). (b,c) Calculated PDOS (b) and simulated d$I$/d$V$ spectra (with two different broadenings, panel c) for the non-magnetic and charge-density wave states. (d) Typical long-range experimental d$I$/d$V$ spectra on 1 and 2 ML VSe$_2$ on NbSe$_2$. (e,f) Simulated d$I$/d$V$ spectra (with two different broadenings, panel e) and calculated PDOS (f) for the ferromagnetic ground state. []{data-label="fig3"}](Fig3.jpg){width="95.00000%"}

The electronic structure and the simulated and experimental STS are shown in Fig. \[fig3\]. First focusing on the experimental spectra for monolayer and bilayer VSe$_2$ in Fig. \[fig3\]d, we observe peaks close to the Fermi level at both positive and negative bias and they are more pronounced for the bilayer compared to the monolayer. Their energy spacing (gap) is $\sim 0.2$ eV with an abrupt edge at positive energy (peak p1) and smoother edge at negative energy that develops to peaks labelled with n1/n2. At larger positive or negative bias, several peaks can be distinguished with the peak positions shifting between the mono- and bilayer spectra.

Starting from the non-magnetic (NM) phase, the band structure shows the partially filled d-band (Fig. \[fig3\]a). There is a flat region between $\Gamma$ and K-points, which also happens to coincide with the Fermi-level. This leads to a strong peak at the Fermi-level in the DOS and also in the simulated STS, which obviously contrasts the experimental spectra. NM phase is also unstable in calculations and may either develop CDW or ferromagnetism, depending on conditions (strain, doping, defects) and on computational parameters (see SI). In the CDW phase, DOS exhibits a clear pseudogap formation, although not at the Fermi-level, but just above it. While the simulated STS would be relatively consistent with the experimental one, we do not detect CDW in the STM topography and also do not observe a hard gap at the Fermi-level as reported in Ref. .

In the ferromagnetic (FM) phase, the spin-up and -down band structures are only shifted in energy as in itinerant FM. Consequently also the DOS and simulated STS show splitting of the peaks w.r.t. NM phase. Due to the splitting of the DOS peak at Fermi-level, the lower branch could explain the peaks n1/n2 with the higher branch falling slightly above Fermi-level (feature p1 in the experimental spectra). The FM phase also seems to yield better agreement between the simulated and experimental STS for the peaks p2 and p3 in the monolayer sample. Overall, in the monolayer sample, several features in the STS spectra seem to have better match with the FM phase than the CDW phase.

![Magnetic measurements on ML VSe$_2$ on NbSe$_2$. (a) Background corrected magnetization curves loops taken at $T=10-300$ K. (b) The temperature dependencies of $M_\mathrm{s}$ and $H_\mathrm{c}$.[]{data-label="fig4"}](Fig4.png){width="50.00000%"}

After the electronic characterization of the samples, we will next focus on their magnetic properties. To explore the magnetic properties of ML VSe$_2$ on NbSe$_2$, we carried out magnetization measurements at various sample temperatures (See Methods for details). All the VSe$_2$ samples measured showed an in-plane magnetic response similar to that shown in Fig. \[fig4\]a, where the loops have a vanishing coercivity and remanence, but show saturation around $200-300$ mT [@Duvjir2018_NanoLett]. A paramagnetic background has been subtracted from this data (see SI Fig. S14; Fig. S15 shows a direct comparison between the substrate and VSe$_2$ responses). As shown in Fig. \[fig4\]b, the coercive field is vey small and it and saturation magnetization are practically independent of temperature in the range of $T=10-300$ K. This behaviour is inconsistent with standard ferromagnetism. While the mechanism is different here, there are other systems that exhibit phenomenologically similar magnetic responses [@Coey_2016; @Coey2019].

It has been argued that there is no conventional ferromagnetism in VSe$_2$ on graphite and graphene substrates[@Chen2018_PRL; @Duvjir2018_NanoLett; @Feng2018_NanoLett; @Wong2019_frustrated; @Fumega2019]. Here, one potentially important difference between our samples and those on graphite and graphene is the absence of a charge-density wave: we find no evidence of CDW in either the STM images or the d$I$/d$V$ spectroscopy (see below). The stability of the CDW according to our DFT calculations is extensively discussed in the SI. Furthermore, in previous studies [@Bonilla2018], monolayer VSe$_2$ sample exhibits a maximum in $M_\mathrm{s}$ and $H_\mathrm{c}$ in a range around $T=100$ K and this non-monotonic behaviour is ascribed to the CDW transition. This is in contrast to our data on $M_\mathrm{s}$ and $H_\mathrm{c}$ as a function of temperature. The absence of CDW transation is also independently determined from the magnetization measurement under zero-field-cooled (ZFC) and field-cooled (FC) regimes, where we do not observe noticeable changes in the slope of the ZFC and FC curves due to the CDW (see SI Fig. S16).

The magnetic response of the sample below the superconducting transition temperature ($T_c = 7.2$ K) is dominated by the superconducting NbSe$_2$ substrate. Magnetization $M(H)$ curves measured at 2 K for bulk NbSe$_2$ and VSe$_2$/NbSe$_2$ (see SI Fig. S16) show superconducting diamagnetic behavior as reported previously[@Valle2017; @Galvan2003]. Moreover, there is no obvious difference between these magnetization curves, which indicates that the signal is mostly dominated by the bulk NbSe$_2$. This conclusion is further supported by the temperature dependent magnetization curves for the zero-field cooling (ZFC) of a bulk NbSe$_2$ and VSe$_2$/NbSe$_2$ heterostructure (see SI). They show a rapid decrease at the onset of the diamagnetic signal below the critical temperature $T_\mathrm{c}$ for both samples and, again, there is no obvious change in $T_\mathrm{c}$, which further indicates that the dominance of the bulk signals. It is worth to mention that $T_\mathrm{c}$ is very sensitive to the magnetic doping, with studies suggesting that $T_\mathrm{c}$ drops rapidly upon metal atom doping[@Naik2013; @Pervin2017]. This suggests that we do not lose any vanadium due to intercalation at the normal growth temperatures. Increasing the growth temperature to $T>300$ $^\circ$C results in intercalation of vanadium, which is clearly seen in atomically resolved images of the NbSe$_2$ surface and results in the loss of the long range order of the CDW [@Chatterjee2015] (see SI Fig. S17 for details).

There is a particular interest in the interaction between the superconducting substrate and the magnetic layer both in terms of the proximity effect induced in the single layer VSe$_2$ and conversely, the effect of the magnetic layer on the underlying superconductor. The superconducting proximity effect can be used to spontaneously drive a non-superconducting material (normal metal) into superconductivity, however, this picture is altered when the superconductor makes a contact with a magnetic layer. In the case of a ferromagnet, the superconducting order parameter is expected to decay exponentially with a very short coherence length $ \zeta _F$ (typically some nm) at the superconductor-ferromagnet (SF) interface. Moreover, the pairing potential $\Delta _p$ inside the ferromagnet shows a strong oscillatory and damped behaviour due to the internal exchange field of the ferromagnet [@Linder2015; @Buzdin2005]. This type of phenomena could also occur in our hybrid VSe$_2$/NbSe$_2$ layers and it will allow us to shed some light on the nature of the magnetism in VSe$_2$ on the atomic scale.

![d$I$/d$V$ spectroscopy on VSe$_@$ layers on NbSe$_2$. (a) Spectra measured on clean NbSe$_2$ and on 1 and 2 ML of VSe$_2$ (spectra offset vertically for clarity). (b) The measured energy gap as a function of the VSe$_2$ thickness. (c) Spectra measured at different locations, including over the edge of a VSe$_2$ island. (d) Topographic STM image and a constant-height d$I$/d$V$ map at zero bias over an edge of a VSe$_2$ island.[]{data-label="fig5"}](Fig5.jpg){width="95.00000%"}

Fig. \[fig5\]a shows the d$I$/d$V$ spectra measured on the NbSe$_2$ substrate as well as VSe$_2$ layers with different thicknesses. On NbSe$_2$, we observe a typical superconducting gap: a pronounced dip in the DOS at the Fermi level and coherence peaks on both sides of the gap[@hudson1999thesis; @Noat2015_nbse2]. The spectra measured on the VSe$_2$ films with thicknesses of 1 ML and 2 ML also shows a superconducting gap, but the gap width is significantly reduced compared to the bare NbSe$_ 2$ substrate. We do not observe oscillatory behaviour of the pairing potential as bilayer VSe$_2$ is not sufficiently thick for this. To further quantify the reduction of the SC gap width, we use a simple analytical relation between a superconducting gap $\Delta$ and decay length $\lambda$ [@Reeg2016; @Belzig1996; @Floetotto2018]: $\Delta \approx \Delta _\mathrm{NbSe_2} e^{-d /\lambda}$, where the $d$ is the VSe$_2$ film thickness and $\Delta_\mathrm{NbSe_2}$ is the gap of bare NbSe$_2$, respectively. We extract the apparent gap widths from the spectra shown in Fig. \[fig5\]a by fitting them (see SI Fig. S18 for details) and plot the results in Fig. \[fig5\]b. The decrease of the energy gap follows an exponential dependence with an decay length of $\lambda = 1$ nm (1 ML thickness is roughly 0.65 nm). This is a measure of the quasi-particle coherence length associated with Andreev reflections. However, this decay length is much shorter than recent experimental results on Bi$_2$Se$_3$ on NbSe$_2$[@Wang2012]: they observe a decay length of $\lambda=2$ nm, which is shown as a dashed blue line in Fig. \[fig5\]b. In the case of Bi$_2$Se$_3$ on Nb, an even much higher value of $\lambda=8.4$ nm was reported[@Floetotto2018]. The faster decay we observe in VSe$_2$ is most easily explained by magnetism of the VSe$_2$ layer. In this case, one would expect a shorter coherence length $ \zeta_F$ governed by the magnetization in the VSe$_2$ layer and not by the diffusion (which is a case in Bi$_2$Se$_3$/NbSe$_2$). The reduced gap is very uniform within the VSe$_2$ islands and also between different islands, which suggests that it is related to bulk properties of ML VSe$_2$ and not simply an effect arising from impurities, vacancies or structural imperfections (e.g. layer edges).

We have also probed the spatially dependent spectroscopic response over the edges of the VSe$_2$ islands and typical spectra are shown in Fig. \[fig5\]c (more results in the SI Fig. S19). The spectra evolve from the typical gapped structure over the NbSe$_2$ into a sharp peak at zero bias at the edge of the VSe$_2$ island. This feature is very localized at the edge of the VSe$_2$ layer. Furthermore, it is inhomogeneously distributed along the edges of VSe$_2$ islands and there are strong intensity variations as illustrated in Fig. \[fig5\]d. In addition to the spatial distribution of the zero bias peak, its width is also strongly position dependent. We observe both zero bias peaks that are confined within the superconducting gap (e.g. Fig. \[fig5\]c), but on some other locations (see SI Fig. S19), its width can be a couple of times larger than the superconducting gap width. Features inside the superconducting gap could arise from Yu-Shiba-Rusinov bands or topological edge modes, but the broader peaks suggest the presence of free unpaired spins giving rise to the Kondo effect[@Franke2011_science; @Heinrich2018_review]. In any case, the reason is likely related to the changes of the gap structure of the underlying superconductor due to local magnetic fields arising from the edges of the VSe$_2$ layer. These results suggest importance of edge effects in the magnetism of the VSe$_2$ layer.

Conclusions {#conclusions .unnumbered}
===========

In conclusion, we have demonstrated high-quality epitaxial growth of VSe$_2$-NbSe$_2$ hybrid structures using MBE. We have observed significant and spatially uniform reduction of the superconducting gap of the NbSe$_2$ substrate on the VSe$_2$ islands with the reduction being thickness dependent and stronger on bilayer VSe$_2$. This would be most naturally explained to result from magnetization of the VSe$_2$ layer. The other electronic and magnetic characterization results are also more consistent with magnetization than with charge density waves. Finally, we observe strongly position-dependent, enhanced d$I$/d$V$ intensity at the Fermi level around the edges of the VSe$_2$ layer suggesting that the atomic-scale structural details of the edge of monolayer VSe$_2$ may contribute to its unusual magnetic response. Finally, our work suggests that it will be possible to combine 2D TMDs with different quantum ground states to stimulate new work in the field of 2D-TMDs hybrids.

Methods
=======

**Sample preparation.** Monolayer VSe$_2$ was grown on NbSe$_2$ substrates by e-beam evaporation of V (99.8 %, Goodfellow Cambridge Ltd.) and simultaneous deposition of atomic Se (99.99 %, Sigma-Aldrich) from a Knudsen cell under ultra-high vacuum conditions (UHV, base pressure $\sim 10^{-10}$). The NbSe$_2$ substrate (HQ Graphene) was cleaved in vacuum and subsequently annealed in ultra-high vacuum at $T=600$ K for 1 h before film growth. VSe$_2$ was grown at a substrate temperature of $T=520-540$ K. The growth was carried out under Se-rich conditions. The excess selenium desorbs from the substrate since the substrate temperature was higher than the evaporation temperature of selenium atoms ($T=393$ K). The samples were either characterized in-situ by XPS and STM or capped with Se by deposition of Se at room temperature before transferring them out of the UHV system.

**XPS measurements.** XPS measurements (Surface Science Instruments SSX-100 spectrometer) were performed using monochromated Al K-alpha radiation with X-ray power of 200 W, a pass energy of 100 eV and a measurement spot size of 1 mm.

**STM measurements.** The STM experiments (Unisoku USM-1300) of the samples were performed at $T=4$ K. STM images were taken in the constant current mode. d$I$/d$V$ spectra were recorded by standard lock-in detection while sweeping the sample bias in an open feedback loop configuration, with a peak-to-peak bias modulation of 5 mV (long-range spectra) or 0.1 mV (short-range spectra of the superconducting gap) at a frequency of 709 Hz.

**Magnetic characterization.** Magnetic characterization was carried out using a Quantum Design Dynacool PPMS operating as a vibrating sample magnetometer (VSM). For magnetization measurements the sample is attached to a quartz rod.

**DFT calculations.** All density-functional theory calculations are carried out in the plane-wave basis in the projector augmented wave framework as implemented in VASP [@kres1; @kres2; @PAW]. In all calculations, we use 500 eV cutoff and $k$-point sampling corresponding to $24\times24$ mesh in the primitive cell. High $k$-point mesh is required to correctly describe e.g. the CDW phases. Further computational details are given in the SI.

V tends to exhibit strong Coulomb correlations, which usually necessitates using either hybrid functionals or +U. It was shown in Ref. [@Fumega2019], that depending on the U-parameter, the monolayer is either in nonmagnetic CDW state or in ferromagnetic no-CDW state. This balance is also affected by the number of electrons [@Fumega2019] and likely also by strain. Further details on the computational results on the interplay between magnetism and CDW stability is given in the SI.

This research made use of the Aalto Nanomicroscopy Center (Aalto NMC) facilities and was supported by the European Research Council (ERC-2017-AdG no. 788185 “Artificial Designer Materials”) and Academy of Finland (Academy professor funding no. 318995 and 320555, Academy postdoctoral researcher no. 309975, and Academy research fellow no. 311058). Our DFT calculations were performed using computer resources within the Aalto University School of Science “Science-IT” project and the Finnish CSC-IT Center for Science.

More experimental and computational details.
