  In the previous paper [Yamada, Chaos, Solitons $\&$ Fractals, {\bf
109},99(2018)], we investigated localization properties of one-dimensional
disordered electronic system with long-range correlation generated by modified
Bernoulli (MB) map. In the present paper, we report localization properties of
phonon in disordered harmonic chains generated by the MB map. Here we show that
Lyapunov exponent becomes positive definite for almost all frequencies $\omega$
except $\omega=0$, and the $B-$dependence changes to exponential decrease for
$B > 2 $, where $B$ is a correlation parameter of the MB map. The distribution
of the Lyapunov exponent of the phonon amplitude has a slow convergence,
different from that of uncorrelated disordered systems obeying a normal
central-limit theorem. Moreover, we calculate the phonon dynamics in the MB
chains. We show that the time-dependence of spread in the phonon amplitude and
energy wave packet changes from that in the disordered chain to that in the
periodic one, as the correlation parameter $B$ increases.
