In a recent breakthrough experiment [Nature (London) 614, 445 (2023)10.1038/s41586-022-05585-1], signatures of Majorana zero modes have been observed in tunnel spectroscopy for a minimal Kitaev chain constructed from coupled quantum dots. However, as Ising anyons, Majoranas' most fundamental property of non-Abelian statistics is yet to be detected. Moreover, the minimal Kitaev chain is qualitatively different from topological superconductors in that it supports Majoranas only at a sweet spot. Therefore, it is not obvious whether non-Abelian characteristics such as braiding and fusion can be demonstrated in this platform with a reasonable level of robustness. In this work, we theoretically propose a protocol for detecting the Majorana fusion rules in an artificial Kitaev chain consisting of four quantum dots. In contrast with the previous proposals for semiconductor-superconductor hybrid nanowire platforms, here we do not rely on mesoscopic superconducting islands, which are difficult to implement in quantum dot chains. To show the robustness of the fusion protocol, we discuss the effects of three types of realistic imperfections on the fusion outcomes, e.g., diabatic errors, dephasing errors, and calibration errors. We also propose a fermion parity readout scheme using quantum capacitance. Our work will shed light on future experiments on detecting the non-Abelian properties of Majorana modes in a quantum dot chain.

Majorana zero modes can appear at the wire ends of a one-dimensional topological superconductor and manifest themselves as a quantized zero-bias conductance peak in the tunneling spectroscopy of normal-superconductor junctions. However, in superconductor-semiconductor hybrid nanowires, zero-bias conductance peaks may arise owing to topologically trivial mechanisms as well, mimicking the Majorana-induced topological peak in many aspects. In this work, we systematically investigate the characteristics of zero-bias conductance peaks for topological Majorana bound states, trivial quasi-Majorana bound states and low-energy Andreev bound states arising from smooth potential variations and disorder-induced subgap bound states. Our focus is on the conductance peak value (i.e., equal to, greater than, or less than 2e2/h), as well as the robustness (plateau- or spike-like) against the tuning parameters (e.g., the magnetic field and tunneling gate voltage) for zero-bias peaks arising from the different mechanisms. We find that for Majoranas and quasi-Majoranas, the zero-bias peak values are no more than 2e2/h, and a quantized conductance plateau forms generically as a function of parameters. By contrast, for conductance peaks due to low-energy Andreev bound states or disorder-induced bound states, the peak values may exceed 2e2/h, and a conductance plateau is rarely observed unless through careful postselection and fine-tuning. Our findings should shed light on the interpretation of experimental measurements on the tunneling spectroscopy of normal-superconductor junctions of hybrid Majorana nanowires.