The proximity of a superconductor   to the surface of a topological insulator (TI) gives rise to a topological superconductor (TS) characterized by the Majorana zero modes. Recently, vortices and Majorana fermions in a magnetic field have been reported  in heterostructures of Bi2Te3/NbSe2. Additionally, Majorana fermions have been studied in Abrikosov lattices   . The Majorana zero modes are neutral excitations which, in an Abrikosov vortex lattice, become a gapless dispersive band    . In a variety of materials, a scanning tunneling microscope (STM) is used to detect charge tunneling.
A number of methods based on sound waves have been used to investigate superconductors     . Ultrasound attenuation studies   , and investigations of the p-wave superconductor Sr2RuO2 have been carried out in ref.  . In the late fifties ultrasound attenuation techniques were used to measure the temperature dependence of the superconducting gap      and recently the techniques have been applied to liquid 3He  and to 4He by  where in the presence of a dislocation we have a direction where low energy excirtations are free to move in analogy with the Majoranas modes in topological superconductors.
The purpose of this paper is to demonstrate that piezoelectric transducers can be used to detect Majorana fermions in an Abrikosov vortex lattice and the dissipative flow trough 4He can be explained.The tunneling amplitude of the Majorana fermions gives rise to a dispersive band  which is detectable. We compute the stress viscosity as a response to an applied strain velocity field   . The explicit dependence of the strain field on the system is obtained from a coordinate transformation  . The linear stress response theory  used for a TS provides information about the Majorana fermions. In order to demonstrate the detection of Majorana fermions we consider a TS       and consider the dissipative flow trough 4He.
In this paper, we have derived the following specific results: (a) We have obtained the vortex lattice solution for a p-wave superconductor. (b) We have derived the stress-strain Hamiltonian and have computed the stress viscosity using the linear response theory. (c) We have identified the sound analog of the Andreev crossed reflection and have obtained the viscosity equivalent to the crossed reflection conductance. (d)We consider the analogues of the Majorana mode in solid 4He.
The structure of the paper is as follows: In Section 2, we show that for an attractive interaction on the surface of a TI a TS is obtained. In the presence of an Abrikosov vortex lattice dispersive Majorana fermions are formed. In Section 3, we present the formation of the Abrikosov vortex lattice. We find dispersive Majorana fermions and quasi-particles in the vortex lattice. A new solution for the p-wave Abrikosov vortex lattice is obtained and discussed in detail. Section 4 is devoted to the derivation of the viscosity tensor for a TS. In Section 5, we compare our results to the one obtained by the ultrasound attenuation technique. In Section 6, we compute the transverse impedance for the TS in a magnetic field. The transverse impedance provides distinct information about the Abrikosov vortex lattice and the Majorana fermions. Our numerical estimates indicate that the signal should be experimentally detectable. Section 7 is devoted to studies of sound wave of the Andreev crossed reflection induced in one dimension by a piezoelectric transducer. In Section 8 we consider the disipative flow of a superfluid of 4He through a solid 4He  . Section 9 contains our main con- clusions.
2. Formation of Majorana Fermions on the Surface of a TI with Attractive Interactions
In this section, we review the formation of a p-wave superconductor on the surface of a TI with attractive interactions. For two space dimensions, the quantum Hall system and the p-wave superconductor    are cha- racterized by the first Chern integer number (which means that the integral of the Berry curvature over a closed manifold is quantized in units of)  . In the presence of an attractive interaction (due to the electron-phonon interaction or proximity to another superconductor), on the surface of a TI a two-dimensional TS emerges. The proximity of a superconductor  to the three dimensional TI gives rise to Majorana zero modes on the surface of the TI. Recently it was reported that the application of a magnetic field on the hete- rostructure Bi2Te3/NbSe2 induces an Abrikosov vortex lattice  .
We propose that a realization of the model introduced in  emerges from the TI surface Hamiltonian in the presence of an attractive interaction. Here, denote Pauli matrices and, are wave vector components. We express the pairing interaction in terms of the field where are the TI surface spinors for the
conduction band, (T stands for
transpose), and is a momentum scale. This representation generates the linear derivatives of the pairing field. For a positive chemical potential in the presence of a magnetic field, it gives rise to a superconductor with vortices. The attractive interaction expressed in terms of the TI spinors gives rise to the p-wave Hamiltonian  which in our case also includes vortices.
Here is the vector potential and is a coupling constant. The pairing field depends on the phase which includes a multivalued part. We perform the gauge transformation:,
and the Hamiltonian is replaced by. The pairing field
has points, , where vanishes, , , where is the multivalued phase.
As a result of the gauge transformation the fermion operators, are replaced by, (and the Hamiltonian is replaced by).
The Hamiltonian without the condensation energy is expressed
in terms of the particle-hole Pauli matrices, and. We introduce the two-component spinor and find:
The spinor  contains two parts, the non-zero mode and the zero mode (Majorana fermions), .
3. TS Abrikosov Vortex Lattice
Next we discuss in detail the non-zero and zero modes of the Abrikosov lattice in a TS.
a) Non-zero modes
In this section we consider the non-zero modes for an Abrikosov vortex lattice in the presence of a magnetic field. The experimental work on the TS Bi2Te3/ NbSe2  shows that an Abrikosov vortex lattice is formed. A vortex lattice is stabilized for superconductors when the penetration depth of the magnetic field is larger than the coherence length   . We find for a single vortex a string-like solution for the effective magnetic field. For the mag- netic field vanishes (d is the vortex lattice constant and is the magnetic penetration depth). Since we are interested in the long distance behavior we can approximate the magnetic field for by a constant field, which is the spatial average around the vortex core with a radius of d.
Following refs.   we solve the p-wave Hamiltonian in a periodic mag- netic field b. The periodicity being. The periodic spinor solution is given by:
, where and
is a two component spinor which is given by the eigenfunction of the p-wave Hamiltonian. We consider a square lattice and the solution is periodic in the y direction with the periodicity. The periodicity in the x direction is achieved by demanding the in- variance of the Hamiltonian under the transformation,. The system has a finite extention L in the x direction. Therefore the value of the momentum in the y direction must be restricted to (to ensure that the states lie in the box).
Due to the spinor structure of the solution it is convenient to solve the problem in momentum space and at the end to impose the periodicity of the
wave function. We represent the periodic solution as:. We use the momentum representation, and find that for the p-wave Hamiltonian in the magnetic field b:
where is the y component of the vector potential in the Landau gauge for the periodic magnetic field. The eigenvectors and eigenvalues are computed next. The ground state is given in terms of the ground state energy of the Harmonic
We note that is dimensionless. At this stage we impose the periodic boundary conditions. This results in replacing, with the
condition,. As a result the eigenspinors and eigenvectors are given in terms of the integers:
The spinor determines the non-zero mode fermion fields. We introduce the annihilation and creation operators, with respect to the exact ground state, which allows us to write:
We have replaced the discrete sites by the coordinate x in order to consider even and odd rows that we need to introduce in the matrix and double the dimension of the spinor. (This needs to be done in order to have the same dimension for the zero and nonzero modes). Taking in consideration also the particle-hole symmetry we use the representation, which is four- dimensional; acts on the particle-hole space and acts on the even and odd row in real space.
b) Zero modes
In this section we consider the zero modes for an Abrikosov vortex lattice. In the absence of the vortex lattice the zero-mode solutions for the Hamiltonian in Equation (2) are given by:
. The function obeys the normalization condition. Due to
the charge conjugation property of the Hamiltonian, the zero modes are Majorana modes. The solution obtained here is similar to the solution given in refs.   for the p-wave superconductors. The explicit form of the kinetic
operator determines the exact form of the amplitude for the zero modes  . We consider the case where the gauge transformed Hamiltonian has Majorana zero modes. The Majorana operators obey, for,
In the second stage we want to discuss the effect of the vortex lattice on the localized Majorana fermions. The effect of the vortex lattice is to delocalize the Majorana fermions and form dispersive Majorana bands. This is a result due to ref.  which showed that Majorana fermions enclose a flux which depends on the number of vortices on a closed polygon. The flux on a polygon of n vortices
is (n is the number of vortices in the polygon)  . For a square vortex lattice with four vortices () per plaquette the flux will be. For the Majorana case we restrict ourselves to plaquettes with four vortices (per
plaquette)  . We consider the effect of the overlap between Majorana fermions given by matrix element   .
where introduces the phase on the bond and determines the flux of per plaquette and is the overlap between the Majorana fermions which is determined from Eqution (7). (The minimum energy for the Abrikosov vortices is obtained for a triangular lattice  . For a square lattice the energy is less favorable, but it is simpler to analyze.)
We choose a gauge for which the hopping constant along columns has positive sign and alternating signs between adjacent rows. This means that, the phase is zero along columns, has positive sign () and alternating signs between adjacent rows ()  . As a result we obtain two flat bands and a third band that is dispersive and gapless    comprising Majorana fermions with the eigenvalues:
For the zero modes we find that the representation of the zero modes is given in terms of the zero mode Majorana operators in the continuum representation (is replaced with the continuum coordinate x, even and odd rows are introduced with the help of the matrix) and spinor eigenfunctions are given in momentum space:
The four components of the spinor in Equation (10) are given by:
We replace the discrete sites by the coordinate x, therefore we introduce the matrix to double the dimension of the spinor. Taking in consideration also the particle-hole symmetry we use the representation which is four-dimensional; acts on the particle-hole space and acts on the even and odd row. We note that a triangular lattice for the Majorana modes has been considered in refs.   .
4. Viscosity Tensor for the Topological Superconductor
In this section we introduce the theory for the dissipative viscosity. In Sec. V and Sec. VI we use this theory to investigate the Abrikosov lattice. The physics of solids  and fluids  provides us the relation between the stress tensor, strain field and the velocity strain field  . We now use this description for quantum fluids in a solid. The combination of the stress tensor resulting from a strain field and the dissipative part of the stress determines the equation:
The strain field is given in terms of the lattice deformation,
,. The viscosity tensor is given by,
and is separated into two parts, , is the symmetric part and is the antisymmetric part  . From the Onsager relations  we know that when the time reversal symmetry is violated, like in the quantum Hall system and the p-wave superconductor case     , we have. For with we have a situation where the strain field generates stress in the perpendicular direction. The stress tensor for quantum fluids is obtained from the invariance of the Lagrangian under an arbitrary local coordinate transformation     . The explicit dependence of the strain field on the lattice deformation determines the coordinate transformation  from which we obtain the stress tensor. In the presence of an elastic deformation the coordinates transform in the fol- lowing way:  (see the Appendix).
This coordinate transformation allows the identification of the stress tensor. Using the invariance of the spinless fields under the coordinate deformation gives,
For the p-wave Hamiltonian given in Equation (2) we find the momentum density and stress tensor induced by the strain fields ,. To linear order in the deformation strain field we obtain the response stress field. Using the invariance, Equation (13), with respect to the coordinate transformation determines the strain-stress Hamiltonian:
In the absence of disclinations in Equation (14).
Using the explicit dependence of the spinor in Equation (2) allows us to represent the stress fields:
We compute the viscosity stress tensor and the dissipative viscosity tensor in quantum fluids. We vary the p-wave Hamiltonian
by a linear strain field, and find from Equation (6) that the variation with respect to satisfies the continuity equation:
Using the linear response theory  with respect to the strain-stress Ha- miltonian given in Equation (14) we obtain:
where is the stress in the Heisenberg representation computed with respect the ground state of the Hamiltonian in Equation (2). Following ref.  we obtain the relation:
The dissipative part of the viscosity tensor is obtained after the analytic continuation. Equation (18) is computed using Wick’s theorem  for the stress which is expressed in terms of the spinor in Equation (15).
5. Application of the Viscosity Tensor to Ultrasound Attenuation
In this section we compare our calculation to the existing ultrasound attenuation method given in the literature    . For the topological superconductors we need to work with the spinor. Due to the dispersive nature of the zero modes, we will have new contributions to the absorption from the mixed pairs and.
In a superconductor the electron-phonon interaction couples to longitudinal as well as transverse phonons. The coupling of the electrons to the transverse phonons in the superconducting phase is less understood. We show that by applying a transverse strain we can obtain non-diagonal response for stress; in this way we study the transverse effect of phonons.
The ultrasound attenuation method measures the superconducting gap. The single particle contribution to the absorption in a superconductor is given by
(is the absorption for the superconductor and is
the absorption in the normal phase). In our case we have contributions from the electrons and the Majorana modes. Due to the magnetic vortex lattice the absorption is given by discrete summations instead of the integration of quasi- particle density of states. The absorption of transverse phonons is obtained from the response of non-diagonal strain tensor. Theoretically we express the strain tensor in terms of the normal modes of the harmonic crystal  [depends only on the crystal phonons] without requiring knowledge of the explicit electron-phonon interaction. The strain tensor is represented in terms of the normal phonon operators and (are the two phonon polarizations for a phonon in the i direction given by the vector in the orthogonal direction to the vector):
In order to compute the dissipation for the vortex lattice Hamiltonian with the eigenvalues and wave function given in Equation ((5), (6)) stress-strain Hamiltonian (the representation of the Hamiltonian in Equation (14)). We need to do it for the non-zero mode part and zero mode. For the non zero mode we have from Equation (5) the representation,
Where is the real space representation of given in Equation (4). Next we use the invariance of the spinors under the coordinates transformation following Equation (13) and the metric transformation given by the Jacobian transformation given in Appendix we obtain the the stress-strain coupling for the non-zero modes. We perform a similar derivation for the zero mode part.
As a result we find to lowest order the strain stress Hamiltonian:
The contribution to the stress-strain Hamiltonian for the no-zero modes is given by the change of the metric given by the Jacobian of the coordinate transformation e (see the Appendix). For the zero mode we have two con- tribution one from the metric change e and the first order derivative originates from the Fourier tranform of the eigenvalue. This will be our new stress-strain Hamitonian which replaces Equation (15).
Using first order perturbation theory we compute the ultrasound attenuation in agreement with refs.    .
Following ref.  , the transverse sound absorption and longitudinal absorption can be represented in terms of the viscosity tensor and defined in Equation (18). In two dimensions we
have for the transverse absorption, where . The longitudinal absorption is given by where. The information
about the crystal enters through the sound velocity (transverse), (longitudinal), crystal density and quantum fluid viscosity. For example the viscosity terms are computed according to Equation (18) using the mode expansion of the spinor which allows us to compute the longitudinal absorption: represents the non-zero mode part [the index means that the two fields which contribute to absorption are only non-zero modes, the symbol means that we have contributions only from and]. Similarly, represents the mixed contribution, a zero mode and a non-zero mode (one field contains the zero mode and the second contains the non-zero mode) whereas represents the contribution when both fields are zero modes.
Using the imaginary time order operator in the imaginary time re- presentation  we find: for the tensor is given as , and .
The particle-hole contribution is given by (the particle- particle contributions are neglected and the symbol means particle-hole). The mixed terms particle-Majorana and hole-Majorana are given by . Using Wick's theorem with the spinor representation given in Equation ((5), (10)) we compute and.
a) Non-zero modes absorption
We consider first the absorption for the particle-hole in the absence of Majorana modes. We find for the dissipative viscosity:
where is the quasi-particle dispersion given in Equation (11). Here denotes the scattering life-time for the quasi-particles, d is the lattice separation between the vortices, is the frequency of the transducer strain field and T is the temperature.
For an s-wave superconductor the absorption agrees with the results given in the literature,. For the present case with the dispersion
, the absorption is controlled by the magnetic field with the ground state energy [see Equation (11)].
b) Absorption due to Majorana modes
The Majorana modes give rise to the particle-hole (Majorana) absorption and the particle-particle (Majorana) absorption. This notation means that the absorption is controlled by a zero and a non-zero mode, “” and the non-zero mode is either a particle-hole “” or a particle-particle “” channel. Both absorptions are controlled by the tunneling amplitude of the dispersive Majorana mode. From the particle-hole (Majorana) absorption we have the combination where the nonzero mode is a particle and the Majorana operator in the fermionic re- presentation [see Equation (5)] is a hole or vice versa. We introduce the scattering life time () and find for the representation:
where is the dispersion of the Majorana fermions given in Equation (4).
The particle-particle (Majorana) absorption (this is the case that a non-zero mode particle and a zero mode Majorana are created) is given by:
where is the dispersion of the Majorana fermions given in Equation (4). From the theory of electromagnetic paramagnetic response in supercon- ductors [see ref.  , Equations (8.47)-(8.50)] we can see the similarity with our results, Equations ((21), (22)). In our case the electric field is replaced by the strain velocity.
6. Transverse Impedance for Topological Superconductors
The transverse impedance is given by:.
It describes the response of the stress field in the direction to an applied strain field in the direction similar to the Hall effect. In the frequency space we have:
where is the real dissipative part and is the ima- ginary part. We now compute the transverse impedance for the TS Abrikosov vortex lattice.
Similar studies have been performed for the superfluid Helium 3He phase. The authors in ref.  have measured the superfluid acoustic impedance of 3He-B coated with a wall of several layers of 4He. The measurement has been performed using the resonance frequency of an ac-cut transducer which oscillates in a shear or longitudinal mode. The coating was used to enhance the specularity of quasi-particle scattering by the wall. In our case we do not have a rough wall and do not approximate the scattering by quasi-classical theory with a random S-matrix. Instead, we use an oscillating wall and compute the dissipative vis- cosity using the linear response theory given by Equations (14)-(18). We make use of the full spectrum of the zero and non-zero modes given in Equations (5)-(11).
a) Impedance for non-zero modes
The dissipative quasi-particle contribution in the absence of Majorana modes is:
In this subsection the superscript () is implicit in the expression for. Here is the TS pairing field, is the momentum of the vortex
lattice with the vortex separation distance, is the frequency of the applied strain field and is the scattering life-time. We find .
The gap parameter controls the dissipative stress. Using Equation (24) with the gap parameter, momentum, and vortex lattice we have and find: , with varying from 1 to 6.
Figure 1 shows that the sound dissipative impedance is controlled by the absorption edge condition [here is the temperature and is the ground state energy determined by the magnetic field b]. Using
the explicit formula, given in Equation (11) we can determine from the impe- dance the magnetic field b and the vortex lattice constant d.
The absorption is given in units of. We plot for the case and as a
function of temperature (corresponds to 1 Kelvin) and the vortex separation is. The function is shown for three different cases: The thin line gives the absorption for the ground state energy
, the thickest line represents the absorption for the ground state energy. The line in between describes
Figure 1. The dissipative part for the particle-hole contribution, , with only shown here.
the situation for. We observe that the absorption edge temperature scales with the energy as a function of the magnetic field and Fermi energy.
b) Impedance for Majorana modes
The Majorana contribution given by the particle-hole [see Equation (21)] depends on the tunneling amplitude. We consider the ground state energies (thin line), (thick line) and (intermediate ground state energy).
The Majorana contribution for the particle-particle part [see Equation (22)] for the same values of and tunneling amplitude as in Figure 2 shows an absorption for low temperatures. Comparing Figure 1 with Figure 2, Figure 3 we observe that the Majorana fermion gives rise to absorption at low temperatures, in a region where the particle-hole absorption (Figure 1) is absent.
Since impurities are always present, it is important to know how to dif- ferentiate between the impurity absorption and the Majorana fermions. Impu- rities will give rise to absorption for frequencies and for temperatures; on the other hand, the Majorana absorption persists at (see Figure 3). The total impedance is given by the sum of contributions in the three Figures 1-3. We therefore conclude that the information about the Majorana modes, the magnetic field and tunneling amplitude can be ob- tained from the viscosity stress measurement.
7. The Sound Wave Analog of Andreev Crossed Reflection
Next we comment on the sound wave analog of the Andreev reflection. Two
Figure 2. The dissipative impedance for Majorana particle-hole contribution , with only shown here. The range of temperature is (temperature corresponds to 0.1 meV and the structure for is an artifact of the numerics).
Figure 3. The dissipative impedance for Majorana particle contribution as a function of (thin line), , (thick line). Only is shown here. The range of temperature is (temperature corresponds to 0.1 meV and the structure for is an artifact of the numerics).
Majorana modes located at the two ends of a p-wave superconductor wire are detectable by piezoelectric transducers representing the sound equivalent of the two-leads experiments which measure the Andreev crossed reflection   . We demonstrate that the same equations which were obtained for the Andreev crossed reflection   induced by a voltage between the two tips are obtained for a sound wave which creates a time-dependent lattice deformation
. Here is the sound deformation in the vicinity of each tip. The lattice deformation acts as a bias field. The voltage field in the two-tip experiment is replaced by a bias field for the sound wave case.
We follow the derivation given in refs.   . For a p-wave (or equivalently a one-dimensional wire with spin-orbit interaction in the proximity of an s-wave superconductor and a magnetic field) with length we have two Majorana
modes localized at and. The fermions for the two tips are represented by and . (Here, , and are the right and left chiral fermions and is the Fermi momentum of the electrons).
Following ref.  we integrate the Majorana fermions and obtain the coupling between the two tips: where is the overlap energy between the two Majoranas. Ignoring the oscillating
terms allows us to simplify the form of. In order to study the response to sound waves we replace, where is the sound deformation induced by the transducer. The deformation field is a function of the transducer frequency. The velocity strain field is given by. The derivative of the effective Hamiltonian with respect to the strain velocity determines the momentum density,
Here is the correlation between the tips; when a voltage is applied between the tips this correlation represents the current operator  .
From the momentum density we compute the dissipative viscosity and the time ordered correlation function,. Following Equation (18) we obtain the viscosity. The viscosity is equivalent to the Andreev crossed reflection conductance obtained when voltage is applied between the tips   . [To compare the two correlation functions we need to replace with and with.]
8. The Dissipative Superfluid Flow through Solid 4He
The series of experiments performed by  show a flow of a superfluid through solid 4He. Our starting point is the model for solid Helium in the broken symmetry phase including phonons and localized quasiparticles pinned by the solid:
We propose to explain the experiment  by applying two principles introduced in this paper : a-The coordinate transformation, b-The identification of the analog of the Majorana modes. We consider a model with two leads connected to a superfluid reservoirs which are described in the Bosonic language of a Luttinger liquid. The Majorana tunneling term is identify with the low energy bosons confined around the dislocation.(The low energy excitation is the analog of the Majorana). Around the dislocation the bosons are free to move in a direction perpendicular to the two-dimensional plane of the dislocations  . The bosons are disordered in the vicinity of the dislocation line. Due to the bending of the dislocation line the bosons binding energy depends on the local curvature. As a result we have a situation where the superfluid reservoirs are coupled to low energy excitations; due to the varying curvature they can be treated as disordered bosons (the density waves will not contribute, due to the large energy needed to excite regular bosons).
First we need to understand the effect of the coordinate transformation: In an ideal crystal the low momentum will give rise to the quasi particles excitations, the momentum corresponding to the inverse lattice separation will give rise to the density wave crystal Hamiltonian with the periodicity (is the crystal periodicity). A supersolid is not formed since the low energy quasi-particles are pinned by the mass of the density excitations. The coupling between the density excitation and the low energy quasi-particles is given by:
Here is the reciprocal lattice vector and are the phonons. is the bosonic density which provides the large pinning mass and prevents the exitations of the quasiparticles. In the vicinity of an edge dislocation the term vanishes. To see this consider an edge dislocation at:
, , , where is the Burgers vector. Per- forming a spatial integration with respect y in thevecinity of the dislocation we find:
As result we have:
The dislocation causes the replacement of the kinetic energy where is the metric tensor intro- duced by the edge dislocation. In the z direction the metric tensor is. As a result the crystal has low energy excitations in the z direction. For a finite density of dislocations we will have a direction perpendicular to the surface of the dislocations where quasi particles can propagate. For each dis- location we will have one dimensional low energy solution. Since the path is not a straight line the curvature of the dislocation line gives rise to bound states (realistic modeling of the free path boson is waveguide tube with a narrow width around the dislocation). The distribution of the curvature for different line dislocations can be considered as a source of disorder. As a result the model of the solid is restricted to the low energy paths. The paths which are not dislocations have high resistivity. Since the paths are in parallel we can neglect high resistivity paths. The conducting pats of the solid is given by, (to treat exactly the disorder we need to use the replicas trick). Restricting ourselves to the free bosons confined to the dislocations line we have:
where is a randomly distributed variable, due to the random curvature. The model for the superfluid reservoirs and is:
. For the right reservoir we have with replaced by and chemical potential.
The coupling between the low energy excitations and the two reservoirs is:
The disipative flow is given by a similar formula to the formula for the Andreev tunneling. The wires are replaced by the reservoirs and Majorana tunneling Hamiltonian by the tunneling of through the path caused by dislocations. As a result we propose the viscosity tensor given by the correlation of the two reservoirs superfluids which tunnels through the solid:
Here is the ground state in the absence of the tunneling Hamiltonian.
The quantity is to comparable with with the Hallock experimental results discussed by  .
In the first part of this paper, we derived the spinor solution for an Abrikosov vortex lattice in a topological superconductor. We then obtained the zero and non-zero mode wave functions. These results have been used to compute the dissipative viscosity, which is obtained as a stress response to an applied velocity strain field. Experimentally one uses two transducers, one for measuring the stress response and the second transducer to generate the strain field. We find in addition to the particle-hole contribution, a viscosity term which reflects the presence of Majorana fermions. In the second part, we analyse the dissipative flow through Probing the p-wave wire with a sound wave one thus finds an effect similar to the Andreev crossed reflection. In last part, we consider the dissipative superfluid flow through solid 4He. We show that the dislocation forms an effective one dimensional zero mode.
Appendix: Invariance of the Hamiltonian under the Coordinate Transformation: Derivation of the Stress-Strain Hamiltonian
The viscosity tensor is obtained from the linear response of electrons in an external field  . In order to accomplish this task we need to identify the elastic analog of the external electromagnetic field. This is accomplished using the invariance of the action under the coordinate transformation.
The unperturbed crystal is described by the coordinates and the deformed crystal by the coordinates. The distortion of the crystal is given by which is caused either by the phonons of the crystal or by an external force. We use the system to describe the orthonormal coordinates for the deformed crystal with the basis vector frame. The unperturbed crystal is described by the Cartesian coordinates with the basis frame vectors. The two coordinate systems in the two frames are related  :
where i represents the Cartesian coordinates and a represents the deformed crystal. When the deformation of the crystal vanishes we have the relation. This allows us to introduce the non-relativistic transformation of the derivatives:
The metric integration for the deformed space is given in terms of the crystal deformation vector field. That is, where e is the Jacobian of the coordinate transformation which for the two-dimensional case is given by:. We replace (the exact relation is given by the matrix equation). We
introduce the notation, and.
When the excitations are caused by the phonons, we use the phonon spectrum of the crystal (in the harmonic representation). Due to the compatibility conditions  we have and for a crystal with disclinations we have.
 Xu, J.P., Liu, C., Wang, M.X., Ge, J., Liu, Z.L., Yang, X., Chen, Y., Liu, Y., Xu, Z.A., Gao, C.L., Qian, D., Zhang, F.C. and Jia, J.F. (2014) Physical Review Letters, 112, Article ID: 217001.
 Magnusson, E.B., Williams, B.H., Mannenti, R., Nerisisyan, M.-S., Peterrer, M.J., Ardavan, A. and Leek, P.J. (2015) Applied Physics Letters, 106, Article ID: 063509.
 Murakawa, S., Tamura, Y., Wada, Y., Wasai, M., Saitoh, M., Aoki, Y., Nomura, R., Okuda, Y., Nagato, Y., Yamamoto, M., Higashitani, S. and Nagai, K. (2009) Physical Review Letters, 103, Article ID: 155301.