Dynamics of driven dissipative Frenkel-Kontorova model is examined by using largest Lyapunov exponent computational technique. Obtained results show that besides the usual way where behavior of the system in the presence of external forces is studied by analyzing its dynamical response function, the largest Lyapunov exponent analysis can represent a very convenient tool to examine system dynamics. In the dc driven systems, the critical depinning force for particular structure could be estimated by computing the largest Lyapunov exponent. In the dc+ac driven systems, if the substrate potential is the standard sinusoidal one, calculation of the largest Lyapunov exponent offers a more sensitive way to detect the presence of Shapiro steps. When the amplitude of the ac force is varied the behavior of the largest Lyapunov exponent in the pinned regime completely reflects the behavior of Shapiro steps and the critical depinning force, in particular, it represents the mirror image of the amplitude dependence of critical depinning force. This points out an advantage of this technique since by calculating the largest Lyapunov exponent in the pinned regime we can get an insight into the dynamics of the system when driving forces are applied.

The devil's staircase structure arising from the complete mode locking of an entirely nonchaotic system, the overdamped dc+ac driven Frenkel-Kontorova model with deformable substrate potential, was observed. Even though no chaos was found, a hierarchical ordering of the Shapiro steps was made possible through the use of a previously introduced continued fraction formula. The absence of chaos, deduced here from Lyapunov exponent analyses, can be attributed to the overdamped character and the Middleton no-passing rule. A comparative analysis of a one-dimensional stack of Josephson junctions confirmed the disappearance of chaos with increasing dissipation. Other common dynamic features were also identified through this comparison. A detailed analysis of the amplitude dependence of the Shapiro steps revealed that only for the case of a purely sinusoidal substrate potential did the relative sizes of the steps follow a Farey sequence. For nonsinusoidal (deformed) potentials, the symmetry of the Stern-Brocot tree, depicting all members of particular Farey sequence, was seen to be increasingly broken, with certain steps being more prominent and their relative sizes not following the Farey rule.

The largest Lyapunov exponent has been examined in the dynamical-mode locking phenomena of the ac+dc driven dissipative Frenkel-Kontorova model with deformable substrate potential. Due to deformation, large fractional and higher order subharmonic steps appear in the response function of the system. Computation of the largest Lyapunov exponent as a way to verify their presence led to the observation of the Farey sequence. In the standard regime, the appearance of half-integer and other subharmonic steps between the large harmonic steps, and their relative sizes follow the Farey construction. In the nonstandard regime, however, the half-integer steps are larger than harmonic ones, and Farey construction is only present in the appearance of higher order subharmonic steps. The examination of Lyapunov exponents has also shown that regardless of the substrate potential or deformation, there was no chaos in the system.