We propose the Particle Swarm Optimization (PSO) as an alternative method for locating periodic orbits in a three--dimensional (3D) model of barred galaxies. We develop an appropriate scheme that transforms the problem of finding periodic orbits into the problem of detecting global minimizers of a function, which is defined on the Poincar\'{e} Surface of Section (PSS) of the Hamiltonian system. By combining the PSO method with deflection techniques, we succeeded in tracing systematically several periodic orbits of the system. The method succeeded in tracing the initial conditions of periodic orbits in cases where Newton iterative techniques had difficulties. In particular, we found families of 2D and 3D periodic orbits associated with the inner 8:1 to 12:1 resonances, between the radial 4:1 and corotation resonances of our 3D Ferrers bar model. The main advantages of the proposed algorithm is its simplicity, its ability to work using function values solely, as well as its ability to locate many periodic orbits per run at a given Jacobian constant.

We use the Smaller Alignment Index (SALI) to distinguish rapidly and with certainty between ordered and chaotic motion in Hamiltonian flows. This distinction is based on the different behavior of the SALI for the two cases: the index fluctuates around non--zero values for ordered orbits, while it tends rapidly to zero for chaotic orbits. We present a detailed study of SALI's behavior for chaotic orbits and show that in this case the SALI exponentially converges to zero, following a time rate depending on the difference of the two largest Lyapunov exponents $\sigma_1$, $\sigma_2$ i.e. $SALI \propto e^{-(\sigma_1-\sigma_2)t}$. Exploiting the advantages of the SALI method, we demonstrate how one can rapidly identify even tiny regions of order or chaos in the phase space of Hamiltonian systems of 2 and 3 degrees of freedom.