We demonstrate that the current-voltage (I-V) characteristics of resistively and capacitively shunted Josephson junctions (RCSJ) hosting localized subgap Majorana states provides a phase sensitive method for their detection. The I-V characteristics of such RCSJs, in contrast to their resistively shunted counterparts, exhibit subharmonic odd Shapiro steps. These steps, owing to their subharmonic nature, exhibit qualitatively different properties compared to harmonic odd steps of conventional junctions. In addition, the RCSJs hosting Majorana bound states also display an additional sequence of steps in the devil staircase structure seen in their I-V characteristics; such sequence of steps make their I-V characteristics qualitatively distinct from that of their conventional counterparts. A similar study for RCSJs with graphene superconducting junctions hosting Diraclike quasiparticles reveals that the Shapiro step width in their I-V curves bears a signature of the transmission resonance phenomenon of their underlying Dirac quasiparticles; consequently, these step widths exhibit a π periodic oscillatory behavior with variation of the junction barrier potential. We discuss experiments which can test our theory.

This paper proposes a method of modeling the dynamic properties of multi-valley semiconductors. The model is applied to the relevant materials GaN, AlN, and InN, which are now known by the general name of III-nitrides. The method is distinguished by economical use of computational resources without significant loss of accuracy and the possibility of application for both dynamic time-dependent tasks and the fields variable in space. The proposed approach is based on solving a system of differential equations, which are known as relaxation ones, and derived from the Boltzmann kinetic equation in the approximation of relaxation time by the function of distribution over k-space. Unlike the conventional system of equations for the concentration of carriers, their pulse and energy, we have used, instead of the energy relaxation equation, an equation of electronic temperature as a measure of the energy of the chaotic motion only. Relaxation times are defined not as integral values from the static characteristics of the material but the averaging of quantum-mechanic speeds for certain types of scattering is used. Averaging was carried out according to the Maxwellian distribution function in the approximation of electronic temperature, as a result of which various mechanisms of dispersion of carriers are taken into consideration through specific relaxation times. The system of equations includes equations in partial derivatives from time and coordinates, which makes it possible to investigate the pulse properties of the examined materials. In particular, the dynamic effect of the "overshoot" in drift velocity and a spatial "ballistic transport" of carriers. The use of Fourier transforms of pulse dependence of the drift carrier velocity to calculate maximum conductivity frequencies is considered. It has been shown that the limit frequencies are hundreds of gigahertz and, for aluminum nitride, exceed a thousand gigahertz