{"id":"5119c730-c0af-435e-a22f-9dc9bca6e832","arxiv_id":"2502.03318","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In Ni/2D material/Ni vertical junctions, the layers next to the metal become metallic through hybridization, while the middle layer of a trilayer retains its intrinsic electronic properties, so conductivity falls sharply with thickness.","lead":"This paper simulates vertical junctions made of nickel electrodes around graphene, boron nitride, and two phases of molybdenum diselenide, using density functional theory and Green's function transport calculations. It finds that the 2D layers touching the metal become metallic, while the central layer of a three-layer stack keeps its intrinsic electronic character and controls the current.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Spin polarization of the Ni electrodes is not specified in the Methods; if the DFT/NEGF runs are spin-unpolarized, the interfacial hybridization that metallizes the 2D layers and sets the barrier heights is not a trustworthy basis for the central claim.","rationale":"The reader's weakest assumption focused on no-vdW/PBE distances, which is a valid concern but not the most load-bearing one for this paper's central claim. The claim is explicitly about Ni-induced orbital hybridization metallizing the 2D layers and setting the tunnel barriers. For ferromagnetic Ni, spin polarization changes the d-band structure and its energy alignment with the 2D states, so an undocumented spin treatment is a concrete, testable threat to the mechanism itself. I do not call the authors dishonest; the omission may be an oversight in the manuscript, but it must be settled before the central claim is accepted. I keep the reader's CONDITIONAL verdict rather than escalating to REJECT because the internal logic is coherent and a single spin-polarized rerun can decide the issue. The reader's concern about vdW/PBE is complementary: both are DFT approximations at the interface, hence 'partial' agreement. The paper's own admission that computed conductivities exceed experiment by orders of magnitude reinforces that the interface description is the fragile part, not the qualitative thickness trend.","tokens_in":10256,"tokens_out":11153,"duration_ms":116183,"concrete_test":"Check the SIESTA input files for spin-polarized flags (e.g., collinear spin with initial Ni magnetic moments); if absent, rerun the monolayer and trilayer hBN/Ni and 2H-MoSe2/Ni cases with spin-polarized PBE using the same pseudopotentials and basis, then compare (i) the Ni-2D distances in Table 1, (ii) the projected density of states at the Fermi level in Figs. 2-3, and (iii) the conductivities in Tables 2-3. If the gap states and conductivity rankings survive, the concern is settled; if the interfacial metallization disappears or the conductivities shift by more than an order of magnitude, the central claim depends on the nonmagnetic approximation and the verdict should be downgraded.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central mechanism is orbital hybridization between Ni d states and the 2D material states, which is claimed to metallize the interfacial layers and determine the tunneling barriers (Figs. 2-4, Tables 1-2). The Methods (Section 2.1 and Tables A1-A2) list PBE, pseudopotentials, basis sets, and smearing, and explicitly discuss excluding vdW and SOC, but never state whether spin polarization is included. Ni(111) is strongly ferromagnetic: its d bands are exchange-split, and the spin-resolved density of states at the Fermi level controls hybridization with C/B/N/Se states. If the calculations are spin-unpolarized, the gap states in Figs. 2-3 and the barriers in Fig. 4 arise from a nonmagnetic Ni electronic structure, and the conductivity rankings in Tables 2-3 could change or even reverse. This is a more direct threat than vdW distances because the proposed mechanism is explicitly electronic: hybridization with Ni d states. The authors report no test of this approximation, and their own comparison with experiment shows computed conductivities orders of magnitude above measured values, so the missing spin treatment cannot be dismissed as a small quantitative effect.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript reports first-principles DFT/NEGF simulations of vertical Ni(111)/2D-material/Ni(111) junctions for graphene, hBN, 1T- and 2H-MoSe2, for monolayer and multilayer (2-4 layers) thicknesses, considering two stacking modes. The central claim is that the 2D layers in direct contact with Ni electrodes become metallic through orbital hybridization with Ni d states, so an insulating monolayer of hBN or semiconducting monolayer of 2H-MoSe2 conducts, while in trilayers the middle layer retains its intrinsic electronic structure and acts as the tunneling barrier, causing a sharp decrease of conductivity with layer number. The authors also provide geometric parameters (distances, interaction energies), electronic structure (DoS, Hartree potentials), and computed 2D conductivities.","tokens_in":10416,"tokens_out":8384,"duration_ms":72447,"significance":"If the central claim is robust, the work has practical implications for understanding the high-resistance state of atomristors and for designing vertical 2D-material devices, since it implies that interface hybridization, not the isolated material's band gap, controls out-of-plane transport. The paper is transparent about its approximations: it explicitly lists the pseudopotentials, basis sets, K-grids, and the exclusion of vdW and SOC, and it acknowledges that the computed three-layer conductivities exceed experimental values by orders of magnitude. It also provides a systematic comparison across four materials and several thicknesses, which is valuable as a qualitative benchmark.","major_comments":[{"comment":"The manuscript never states whether spin polarization is included in the DFT and NEGF calculations. Ni(111) is strongly ferromagnetic, with a magnetic moment of approximately 0.6 μB per atom, and the exchange-split Ni d bands are precisely the states that hybridize with the C, B, N, and Se orbitals to create the interfacial gap states (Figs. 2–5) and set the barrier heights (Fig. 4 and Table 1). If the simulations are spin-unpolarized, the central claim that monolayer hBN and 2H-MoSe2 are metallized by hybridization, and the resulting conductivity rankings in Tables 2 and 3, would rest on a nonmagnetic Ni electronic structure. Please state the spin treatment explicitly and provide at least one spin-polarized test (e.g., for the Ni|hBN|Ni junction) showing that the gap states and barrier heights are qualitatively unchanged.","section":"Section 2.1 (Methods)"},{"comment":"The exclusion of van der Waals corrections is justified by the assertion that strong coupling dominates, but no supporting calculation is provided. The electrode–2D material distances in Table 1 are very short (e.g., 2.04–2.08 Å for the strongly coupled stacking of graphene and hBN), and the tunneling barriers in Fig. 4 and the conductivities in Table 2 depend sensitively on these distances. A vdW-corrected calculation for at least one system (e.g., Ni|Gr|Ni or Ni|hBN|Ni) is needed to establish that the hybridization-driven metallization and the barrier heights are robust against the choice of exchange-correlation functional; otherwise the agreement with the qualitative picture may be coincidental.","section":"Section 2.1 and Table 1"},{"comment":"The paper reports convergence tolerances for the self-consistent DFT loop (10 mÅ, 5×10−3 e/Å3, 100 meV) but gives no convergence tests for the transport quantities themselves, such as the K-grid for TRANSIESTA/TBtrans, the number of bias points used for the spline interpolation, or the basis-set size. The monolayer conductivities in Table 2 span only a factor of 2.5, so without error estimates it is unclear whether the ordering Gr > 1T-MoSe2 > 2H-MoSe2 > hBN is numerically significant. Please provide at least one convergence check (e.g., doubling the TRANSIESTA K-grid and/or basis size for one device) and report the resulting change in conductivity.","section":"Section 2.1 and Tables 2–3"},{"comment":"The text in Section 3.1 states that the DoS of the Ni atom plane closest to the 2D material is also provided, with coinciding peaks proving hybridization, but Fig. 2 shows only the total DoS of the 2D material (and the free-standing reference). The Ni-projected DoS does not appear anywhere in the paper. Either add the Ni-plane DoS to Fig. 2 (or Fig. 3) or remove the sentence, because the claim of hybridization-driven metallization is currently supported only by the appearance of gap states in the 2D material, not by direct evidence of Ni–2D orbital overlap.","section":"Section 3.1 and Fig. 2"}],"minor_comments":[{"comment":"The rows for 1T-MoSe2 are missing for 2L and 4L; please state in the text whether these systems were not calculated, and why.","section":"Table 3"},{"comment":"The statement that the middle layer 'preserves its intrinsic electronic structure' would be better supported by overlaying the free-standing DoS in the same panel; currently the comparison is only implicit.","section":"Fig. 5"},{"comment":"The SISL Python library is used for post-processing but is not cited; please add the appropriate reference (e.g., N. Papior, sisl, Zenodo, 2020).","section":"References"},{"comment":"The comparison with experiment is only qualitative ('still much higher than experimentally observed [9]'). Please report the experimental conductivity values from ref. [9] so the reader can gauge the discrepancy quantitatively.","section":"Section 3.2"},{"comment":"The phrase 'two phases of molybdenum diselenide' could be made clearer by naming 1T and 2H explicitly in the abstract itself.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":"The most pressing issue is the missing statement on spin polarization; if the Ni electrodes were treated as nonmagnetic, the central hybridization mechanism could be qualitatively wrong. I encourage the editor to require the authors to clarify this before publication. The paper is within scope for a mesoscopic/electronic-structure journal, and the authors are honest about the limitations."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a solid, useful comparative DFT/NEGF study of Ni/2D/Ni vertical junctions, and the central qualitative message—interfacial layers metallize, inner layers control tunneling—is convincing. But the manuscript never says whether the Ni calculations include spin polarization, and that is not a detail: Ni(111) is ferromagnetic, and the d-band exchange splitting is precisely what sets the hybridization that creates the gap states and barrier heights. If the runs are spin-unpolarized, the quantitative rankings in Tables 2-3 could shift or even reverse. The authors should test this or justify it.\n\nWhat is genuinely new: the systematic dataset across graphene, hBN, 2H- and 1T-MoSe2, with two stackings each and up to four layers, all in the same simulation setup. The PAO basis is optimized against plane-wave band structures in a ±3 eV window, which is a real effort toward reproducibility. There are no fitted parameters; the conductivities come directly from DFT/NEGF. The thickness trend—monolayer conducts, trilayer does not—is consistent across all four materials and matches the DoS and potential-barrier plots. That is a clean, useful result for atomristor and vertical-interconnect modeling.\n\nSoft spots, in order of severity. First, the missing spin treatment. The Methods list PBE, pseudopotentials, basis, smearing, but never spin. For a Ni(111) contact this is a load-bearing omission. Second, vdW corrections are excluded with the argument that strong coupling dominates; that may be true, but the interface distances in Table 1 are the parameters that set the hybridization, and no test is reported. Third, no convergence checks or error bars are given for any quantity; the conductivities are single numbers. Fourth, the authors themselves note the 3-layer conductivities are orders of magnitude above experiment; they attribute this to absent experimental factors, which is plausible but unquantified. Fifth, only the energetically favorable stacking is used for transport, which could bias cross-material comparisons. No input files or data are released.\n\nNone of this kills the qualitative conclusion; the layer-resolved DoS and potential profiles are internally consistent. But the paper's own quantitative claims and design rules are conditional on the missing checks. A referee should ask for a spin-polarized test, a vdW test, convergence data, and at least the key input files. This is a paper worth engaging, not desk rejecting.","headline":"Useful comparative DFT/NEGF dataset on Ni/2D/Ni junctions, but the missing spin treatment of Ni is a real gap; qualitative message holds, quantitative rankings do not.","tokens_in":11007,"tokens_out":3313,"would_cite":true,"duration_ms":30370,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In Ni-contacted vertical junctions, a monolayer of an insulating or semiconducting 2D material becomes metallic and conducts.","keywords":["vertical transport","metal-2D interface","orbital hybridization","density functional theory","non-equilibrium Green's function","hBN","MoSe2","atomristor"],"falsifier":"Measure the zero-bias conductance of a clean Ni|hBN|Ni junction with one, two, and three hBN layers and compare with the simulated values: the monolayer should be conductive and the drop from monolayer to trilayer should be roughly three orders of magnitude; alternatively, repeat the DFT with vdW-corrected functionals and check whether the interfacial distance and barrier heights change enough to reverse the conductivity ranking.","tokens_in":10018,"feed_emoji":"⚡","tokens_out":4611,"duration_ms":41369,"temperature":0.7,"pith_summary":"This paper uses density functional theory and non-equilibrium Green's function transport simulations to establish that in nickel-contacted vertical junctions, the electronic structure of the 2D layer is set by the contact, not by the isolated material. A single monolayer of insulating hBN or semiconducting 2H-MoSe2 becomes metallic through orbital hybridization with the nickel electrodes, with new states appearing inside the gap, so the monolayer conducts. In trilayers, the outer layers are metallized and act as passivating contacts while the middle layer keeps its intrinsic band structure and becomes the tunneling barrier, which explains why conductivity drops sharply with thickness. The result matters for atomristors and vertical beyond-CMOS devices because it says the high-resistance state is controlled by the inner-layer barrier and by interface hybridization rather than by the intrinsic gap of the 2D material alone.","feed_headline":"A single hBN monolayer turns conductive between nickel contacts","feed_subtitle":"Contacts metallize the first layer; in trilayers the intrinsic middle layer blocks current.","key_machinery":"The mechanism is orbital hybridization at the metal/2D interface: nickel d states hybridize with the p states of carbon, nitrogen, boron, or selenium, producing gap states and metallizing the first 2D layer, which then acts as a passivating contact. The counteracting element is the intrinsic middle layer in multilayers, which keeps its isolated band gap and imposes a tunneling barrier that controls the high-resistance state. The authors quantify both through layer-resolved densities of states, planar-averaged Hartree potentials, and conductance normalized by device area.","core_discovery":"The central claim is that out-of-plane transport through metal/2D material/metal junctions is governed by the metal-induced hybridization at the interfaces and by the intrinsic electronic structure of the inner layers, not by the band structure of the isolated monolayer. For monolayers of hBN and 2H-MoSe2, the paper finds that contact with Ni(111) creates hybridized states inside the gap, making the layer metallic and conductive; for graphene the interface coupling further lowers the tunneling barrier and gives the highest conductivity, while hBN has the highest barrier and the lowest conductivity. For trilayers, the layers in direct contact with the electrodes behave like the monolayer case while the middle layer remains essentially intrinsic, so the conductivity falls by orders of magnitude as thickness increases from one to four layers. This is why the paper concludes that a single hBN monolayer is conductive inside non-passivated metallic contacts despite being an insulator when isolated.","pith_inferences":["If the hybridization picture is right, the contact metal choice is as important as the 2D material: metals with lower reactivity or weaker d-p coupling should preserve the intrinsic gap better and restore a true tunnel barrier at monolayer thickness.","The strong stacking sensitivity of graphene and hBN suggests that rotational alignment between the two electrodes could be a knob to engineer barrier height, testable by twisting the contacts relative to each other.","The discrepancy with experiments points to a testable hypothesis: intentional contamination or oxidation at the interface should increase the effective barrier and recover the experimentally low conductance, which could be checked by controlled exposure studies."],"forward_implications":["A monolayer of an insulating or semiconducting 2D material between non-passivated metal contacts is not a tunneling barrier; it conducts.","The high-resistance state of vertical devices is set by the number of 2D layers, through the intrinsic middle layers, so thickness can tune ON/OFF behavior.","Stacking orientation matters most when the supercell is small (graphene, hBN): atomic alignment can change interaction energy and interface distance by up to about 0.8 Å.","For semiconducting materials, conductivity should go to zero with increasing thickness, while for metallic 2D materials it saturates to a finite value set by interlayer coupling.","The simulated trilayer conductivities exceed experimental values by orders of magnitude, implying that real devices have additional resistances such as multilayer regions, residues, or oxide."],"supporting_citations":[{"why":"Provides the experimental atomristor devices whose high-resistance state the simulations are compared with.","marker":"[9]"},{"why":"SIESTA is the DFT code used for all structural relaxations and electronic structure calculations.","marker":"[22]"},{"why":"TRANSIESTA and TBtrans are the non-equilibrium Green's function implementations used to compute conductance.","marker":"[23]"},{"why":"The PBE exchange-correlation functional is the DFT approximation underlying all computed properties.","marker":"[24]"},{"why":"Supplies the basis-optimization protocol that matches plane-wave band structures within a set error tolerance.","marker":"[27]"},{"why":"The PseudoDojo norm-conserving pseudopotentials are used with SIESTA for the DFT calculations.","marker":"[25, 26]"}],"fun_headline_variants":["Insulating hBN monolayer turns metallic under nickel contacts","Metal contacts turn single hBN layer conductive, trilayers block","One layer conducts, three layers block: metal/2D/metal junctions","Metal-induced states make hBN monolayer conductive in junctions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The computed interfacial distances, hybridization strengths, and barrier heights come from DFT with the PBE functional and no van der Waals corrections, so if those distances or gaps are off, the conductivity ordering could change.","fun_headline_variants_meta":{"raw":{"variants":["Insulating hBN monolayer turns metallic under nickel contacts","Metal contacts turn single hBN layer conductive, trilayers block","One layer conducts, three layers block: metal/2D/metal junctions","Metal-induced states make hBN monolayer conductive in junctions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000295,"raw_usage":{"total_tokens":1666,"prompt_tokens":846,"completion_tokens":820,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":462,"completion_tokens_details":{"reasoning_tokens":750}},"tokens_in":462,"tokens_out":820,"duration_ms":7065,"temperature":1.0,"reasoning_tokens":750,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T05:08:19.642638+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the zero-bias conductance of a clean Ni|hBN|Ni junction with one, two, and three hBN layers and compare with the simulated values: the monolayer should be conductive and the drop from monolayer to trilayer should be roughly three orders of magnitude; alternatively, repeat the DFT with vdW-corrected functionals and check whether the interfacial distance and barrier heights change enough to reverse the conductivity ranking.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the experimental atomristor devices whose high-resistance state the simulations are compared with."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"SIESTA is the DFT code used for all structural relaxations and electronic structure calculations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"TRANSIESTA and TBtrans are the non-equilibrium Green's function implementations used to compute conductance."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The PBE exchange-correlation functional is the DFT approximation underlying all computed properties."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the basis-optimization protocol that matches plane-wave band structures within a set error tolerance."}],"review_version":1}