  We study the decay and oscillations of Majorana fermion wavefunctions and
ground state (GS) fermion parity in one-dimensional topological superconducting
lattice systems. Using a Majorana transfer matrix method, we find that Majorana
wavefunction properties are encoded in the associated Lyapunov exponent, which
in turn is the sum of two independent components: a `superconducting component'
which characterizes the gap induced decay, and the `normal component', which
determines the oscillations and response to chemical potential configurations.
The topological phase transition separating phases with and without Majorana
end modes is seen to be a cancellation of these two components. We show that
Majorana wavefunction oscillations are completely determined by an underlying
non-superconducting tight-binding model and are solely responsible for GS
fermion parity switches in finite-sized systems. These observations enable us
to analytically chart out wavefunction oscillations, the resultant GS parity
configuration as a function of parameter space in uniform wires, and special
parity switch points where degenerate zero energy Majorana modes are restored
in spite of finite size effects. For disordered wires, we find that band
oscillations are completely washed out leading to a second localization length
for the Majorana mode and the remnant oscillations are randomized as per
Anderson localization physics in normal systems. Our transfer matrix method
further allows us to i) reproduce known results on the scaling of mid-gap
Majorana states and demonstrate the origin of its log-normal distribution, ii)
identify contrasting behavior of disorder-dependent GS parity switches for the
cases of even versus odd number of lattice sites, and iii) chart out the GS
parity configuration and associated parity switch points as a function of
disorder strength.
