Recently, many orbital studies in barred galaxy potentials have revealed the existence of orbits which are not trapped around x1 tree orbits, but could be potentially appropriate building blocks for bars. These findings question the uniqueness of the x1 family as the standard paradigm of orbital motion in galactic bars. The main goal of this paper is to investigate the role that such orbits could play in shaping the morphology of bars. We trace the morphological patterns appearing in the face-on and edge-on views of the non-periodic orbits presented in these studies and we show that they are introduced in the system by second type bifurcations of x1. For this purpose, we use a typical 3D Ferrers bar model and follow the radial and vertical bifurcations of the x1 family considered as being multi periodic, with multi = 2, 3, 5. The variation of the stability indices of x1 in the multi = 2, 3 cases give us also the four- and six- periodic orbits, respectively. We tabulate these orbits including all information necessary to assess their role as appropriate building blocks. We discuss their stability and their extent, as well as their size and morphological evolution, as a function of energy. We conclude that even the most important of the multi-periodic orbits presented in Tables 2 to 5 are less appropriate building blocks for bars than the families of the x1 tree at the same energy.

A method which we have developed for determining corotation radii, has allowed us to map in detail the radial resonant structures of barred spiral galaxies. Here we have combined this information with new determinations of the bar strength and the pitch angle of the innermost segment of the spiral arms to find relationships between these parameters of relevance to the dynamical evolution of the galaxies. We show how (1) the bar mass fraction, (2) the scaled bar angular momentum, (3) the pitch angle, and (4) the shear parameter vary along the Hubble sequence, and we also plot along the Hubble sequence (5) the scaled bar length, (6) the ratio of bar corotation radius to bar length, (7) the scaled bar pattern speed, and (8) the bar strength. It is of interest to note that the parameters (2), (5), (6), (7), and (8) all show breaks in their behaviour at type Scd. We find that bars with high shear have only small pitch angles, while bars with large pitch angles must have low shear; we also find a generally inverse trend of pitch angle with bar strength. An inference which at first seems counter-intuitive is that the most massive bars rotate most slowly but have the largest angular momenta. Among a further set of detailed results we pick out here the 2:1 ratio between the number of spiral arms and the number of corotations ouside that of the bar. These results give a guideline to theories of disc-bar evolution.

The PERLAS potential has been successfully used in many studies related with the dynamics of the spiral arms \textit{on} the equatorial plane of normal (non-barred) spiral galaxies. In the present work we extend these studies by investigating the three-dimensional dynamics of the spiral arms in the same type of potential. We consider a typical open, logarithmic, spiral pattern of pitch angle 25$^{\circ}$ and we examine the stellar orbits that can support it as the ratio of the masses of the spiral over the disk component ($M_{s}/M_{d}$) varies. We indicate the families of `three-dimensional' periodic orbits that act as the backbone of the spiral structure and we discuss their stability in the models we present. We study further the quasi-periodic and non-periodic orbits in general that follow spiral-supporting orbits as the $M_{s}/M_{d}$ ratio increases. We find that a bisymmetric spiral with 25$^{\circ}$ pitch angle is better supported by orbits in models with $0.03\lessapprox M_{s}/M_{d} \lessapprox 0.07$. In these cases a strong spiral pattern is supported between the radial 2:1 and 4:1 resonances, while local enhancements of the imposed spirals are encountered in some models between 4:1 and corotation. A characteristic bar-like structure is observed in all models at radii smaller than the radius of the 2:1 resonance.

We present and discuss the orbital content of a rather unusual rotating barred galaxy model, in which the three-dimensional (3D) family, bifurcating from x1 at the 2:1 vertical resonance with the known "frown-smile" side-on morphology, is unstable. Our goal is to study the differences that occur in the phase space structure at the vertical 2:1 resonance region in this case, with respect to the known, well studied, standard case, in which the families with the frown-smile profiles are stable and support an X-shaped morphology. The potential used in the study originates in a frozen snapshot of an $N$-body simulation in which a fast bar has evolved. We follow the evolution of the vertical stability of the central family of periodic orbits as a function of the energy (Jacobi constant) and we investigate the phase space content by means of spaces of section. The two bifurcating families at the vertical 2:1 resonance region of the new model change their stability with respect to that of most studied analytic potentials. The structure in the side-on view that is directly supported by the trapping of quasi-periodic orbits around 3D stable periodic orbits has now an infinity symbol (i.e. $\infty$-type) profile. However, the available sticky orbits can reinforce other types of side-on morphologies as well. In the new model, the dynamical mechanism of trapping quasi-periodic orbits around the 3D stable periodic orbits that build the peanut, supports the $\infty$-type profile. The same mechanism in the standard case supports the X shape with the frown-smile orbits. Nevertheless, in both cases (i.e. in the new and in the standard model) a combination of 3D quasi-periodic orbits around the stable x1 family with sticky orbits can support a profile reminiscent of the shape of the orbits of the 3D unstable family existing in each model.

We investigate regular and chaotic two-dimensional (2D) and three-dimensional (3D) orbits of stars in models of a galactic potential consisting in a disk, a halo and a bar, to find the origin of boxy components, which are part of the bar or (almost) the bar itself. Our models originate in snapshots of an N-body simulation, which develops a strong bar. We consider three snapshots of the simulation and for the orbital study we treat each snapshot independently, as an autonomous Hamiltonian system. The calculated corotation-to-bar-length ratios indicate that in all three cases the bar rotates slowly, while the orientation of the orbits of the main family of periodic orbits changes along its characteristic. We characterize the orbits as regular, sticky, or chaotic after integrating them for a 10 Gyr period by using the GALI$_2$ index. Boxiness in the equatorial plane is associated either with quasi-periodic orbits in the outer parts of stability islands, or with sticky orbits around them, which can be found in a large range of energies. We indicate the location of such orbits in diagrams, which include the characteristic of the main family. They are always found about the transition region from order to chaos. By perturbing such orbits in the vertical direction we find a class of 3D non-periodic orbits, which have boxy projections both in their face-on and side-on views.

Using N-body simulations we study the structures induced on a galactic disc by repeated flybys of a companion in decaying eccentric orbit around the disc. Our system is composed by a stellar disc, bulge and live dark matter halo, and we study the system's dynamical response to a sequence of a companion's flybys, when we vary i) the disc's temperature (parameterized by Toomre's Q-parameter) and ii) the companion's mass and initial orbit. We use a new 3D Cartesian grid code: MAIN (Mesh-adaptive Approximate Inverse N-body solver). The main features of MAIN are reviewed, with emphasis on the use of a new Symmetric Factored Approximate Sparse Inverse (SFASI) matrix in conjunction with the multigrid method that allows the efficient solution of Poisson's equation in three space variables. We find that: i) companions need to be assigned initial masses in a rather narrow window of values in order to produce significant and more long-standing non-axisymmetric structures (bars and spirals) in the main galaxy's disc by the repeated flyby mechanism. ii) a crucial phenomenon is the antagonism between companion-excited and self-excited modes on the disc. Values of $Q >1.5$ are needed in order to allow for the growth of the companion-excited modes to prevail over the the growth of the disc's self-excited modes. iii) We give evidence that the companion-induced spiral structure is best represented by a density wave with pattern speed nearly constant in a region extending from the ILR to a radius close to, but inside, corotation.

We investigate the orbital dynamics of a \textit{barred-spiral} model when the system is rotating slowly and corotation is located beyond the end of the spiral arms. In the characteristic of the central family of periodic orbits we find a "bistable region". In the response model we observe a ring surrounding the bar and spiral arms starting tangential to the ring. This is a morphology resembling barred-spiral systems with inner rings. However, the dynamics associated with this structure in the case we study is different from that of a typical bar ending close to corotation. The ring of our model is round, or rather elongated perpendicular to the bar. It is associated with a folding (an "S" shaped feature) of the characteristic of the central family, which is typical in bistable bifurcations. Along the "S" part of the characteristic we have a change in the orientation of the periodic orbits from a x1-type to a x2-type morphology. The orbits populated in the response model change rather abruptly their orientation when reaching the lowest energy of the "S". The spirals of the model follow a standard "precessing ellipses flow" and the orbits building them have energies beyond the "S" region. The bar is structured mainly by sticky orbits from regions around the stability islands of the central family. This leads to the appearance of X-features in the bars \textit{on} the galactic plane. Such a bar morphology appears in the unsharp-masked images of some moderately inclined galaxies.

The investigation of the phase-space properties of structures encountered in a dynamical system is essential for understanding their formation and enhancement. In the present paper we explore the phase space in energy intervals where we have orbits that act as building blocks for boxy-peanut (b/p) and "{\sf X}-shaped" structures in rotating potentials of galactic type. We underline the significance of the rotational tori around the 3D families x1v1 and x1v1$^{\prime}$ that have been bifurcated from the planar x1 family. These tori play a multiple role: (i) They belong to quasi-periodic orbits that reinforce the local density. (ii) They act as obstacles for the diffusion of chaotic orbits and (iii) they attract a large number of chaotic orbits that become sticky to them. There are also bifurcations of unstable families (x1v2, x1v2$^{\prime}$). Their unstable asymptotic curves wind around the x1v1 and x1v1$^{\prime}$ tori generating orbits with hybrid morphologies between that of x1v1 and x1v2. In addition, a new family of multiplicity 2, called x1mul2, is found to be important for the peanut construction. Our work shows also that there are peanut-supporting orbits before the vertical ILR. Non-periodic orbits associated with the x1 family secure this contribution as well as the support of b/p structures at several other energy intervals. Non-linear phenomena associated with complex instability of single and double multiplicity families of periodic orbits show that these structures are not interrupted in regions where such orbits prevail. Depending on the main mechanism behind their formation, boxy bulges exhibit different morphological features. Finally our analysis indicates that "X" features shaped by orbits in the neighbourhood of x1v1 and x1v1$^{\prime}$ periodic orbits are pronounced only in side-on or nearly end-on views of the bar.

We study the dynamical mechanisms that reinforce the formation of boxy structures in the \textit{inner} regions, roughly in the middle, of bars observed nearly \textit{face-on}. Outer boxiness, at the ends of the bars, is usually associated with orbits at the inner, radial 4:1 resonance region and can be studied with 2D dynamics. However, in the middle of the bar dominate 3D orbits that give boxy/peanut bulges in the edge-on views of the models. In the present paper we show that 3D quasi-periodic, as well as 3D chaotic orbits sticky to the x1v1 and x1v1$^{\prime}$ tori, especially from the Inner Lindblad Resonance (ILR) region, have boxy projections on the equatorial plane of the bar. The majority of vertically perturbed 2D orbits, initially on the equatorial plane in the ILR resonance region, enhance boxy features in face-on bars. Orbits that build a bar by supporting sharp "{\sf X}" features in their side-on views at energies \textit{beyond} the ILR, may also have a double boxy character. If populated, the extent of the inner boxiness along the major axis is about the same with that of the peanut supporting orbits in the side-on views. At any rate these orbits do not obscure the observation of the boxy orbits of the ILR region in the face-on views, as they contribute more to the surface density at the sides of the bar than to their central parts.

We study the dynamics in the neighborhood of fixed points in a 4D symplectic map by means of the color and rotation method. We compare the results with the corresponding cases encountered in galactic type potentials and we find that they are in good agreement. The fact that the 4D phase space close to fixed points is similar to the 4D representations of the surfaces of section close to periodic orbits, indicates an archetypical 4D pattern for each kind of (in)stability, not only in 3D autonomous Hamiltonian systems with galactic type potentials but for a larger class of dynamical systems. This pattern is successfully visualized with the method we use in the paper.

We study in detail the structure of phase space in the neighborhood of stable periodic orbits in a rotating 3D potential of galactic type. We have used the color and rotation method to investigate the properties of the invariant tori in the 4D spaces of section. We compare our results with those of previous works and we describe the morphology of the rotational, as well as of the tube tori in the 4D space. We find sticky chaotic orbits in the immediate neighborhood of sets of invariant tori surrounding 3D stable periodic orbits. Particularly useful for galactic dynamics is the behavior of chaotic orbits trapped for long time between 4D invariant tori. We find that they support during this time the same structure as the quasi-periodic orbits around the stable periodic orbits, contributing however to a local increase of the dispersion of velocities. Finally we find that the tube tori do not appear in the 3D projections of the spaces of section in the axisymmetric Hamiltonian we examined.

We study the dynamics in the neighborhood of simple and double unstable periodic orbits in a rotating 3D autonomous Hamiltonian system of galactic type. In order to visualize the four dimensional spaces of section we use the method of color and rotation. We investigate the structure of the invariant manifolds that we found in the neighborhood of simple and double unstable periodic orbits in the 4D spaces of section. We consider orbits in the neighborhood of the families x1v2, belonging to the x1 tree, and the z-axis (the rotational axis of our system). Close to the transition points from stability to simple instability, in the neighborhood of the bifurcated simple unstable x1v2 periodic orbits we encounter the phenomenon of stickiness as the asymptotic curves of the unstable manifold surround regions of the phase space occupied by rotational tori existing in the region. For larger energies, away from the bifurcating point, the consequents of the chaotic orbits form clouds of points with mixing of color in their 4D representations. In the case of double instability, close to x1v2 orbits, we find clouds of points in the four dimensional spaces of section. However, in some cases of double unstable periodic orbits belonging to the z-axis family we can visualize the associated unstable eigensurface. Chaotic orbits close to the periodic orbit remain sticky to this surface for long times (of the order of a Hubble time or more). Among the orbits we studied we found those close to the double unstable orbits of the x1v2 family having the largest diffusion speed.

This paper discusses phase space structures encountered in the neighborhood of periodic orbits with high order multiplicity in a 3D autonomous Hamiltonian system with a potential of galactic type. We consider 4D spaces of section and we use the method of color and rotation [Patsis and Zachilas 1994] in order to visualize them. As examples we use the case of two orbits, one 2-periodic and one 7-periodic. We investigate the structure of multiple tori around them in the 4D surface of section and in addition we study the orbital behavior in the neighborhood of the corresponding simple unstable periodic orbits. By considering initially a few consequents in the neighborhood of the orbits in both cases we find a structure in the space of section, which is in direct correspondence with what is observed in a resonance zone of a 2D autonomous Hamiltonian system. However, in our 3D case we have instead of stability islands rotational tori, while the chaotic zone connecting the points of the unstable periodic orbit is replaced by filaments extending in 4D following a smooth color variation. For more intersections, the consequents of the orbit which started in the neighborhood of the unstable periodic orbit, diffuse in phase space and form a cloud that occupies a large volume surrounding the region containing the rotational tori. In this cloud the colors of the points are mixed. The same structures have been observed in the neighborhood of all m-periodic orbits we have examined in the system. This indicates a generic behavior.

We study the orbital behavior at the neighborhood of complex unstable periodic orbits in a 3D autonomous Hamiltonian system of galactic type. At a transition of a family of periodic orbits from stability to complex instability (also known as Hamiltonian Hopf Bifurcation) the four eigenvalues of the stable periodic orbits move out of the unit circle. Then the periodic orbits become complex unstable. In this paper we first integrate initial conditions close to the ones of a complex unstable periodic orbit, which is close to the transition point. Then, we plot the consequents of the corresponding orbit in a 4D surface of section. To visualize this surface of section we use the method of color and rotation [Patsis and Zachilas 1994]. We find that the consequents are contained in 2D "confined tori". Then, we investigate the structure of the phase space in the neighborhood of complex unstable periodic orbits, which are further away from the transition point. In these cases we observe clouds of points in the 4D surfaces of section. The transition between the two types of orbital behavior is abrupt.

In a series of papers we study the stellar dynamics of the grand design barred-spiral galaxy NGC~1300. In the first paper of this series we estimate the gravitational potential and we give it in a form suitable to be used in dynamical studies. The estimation is done directly from near-infrared observations. Since the 3D distribution of the luminous matter is unknown, we construct three different general models for the potential corresponding to three different assumptions for the geometry of the system, representing limiting cases. A pure 2D disc, a cylindrical geometry (thick disc) and a third case, where a spherical geometry is assumed to apply for the major part of the bar. For the potential of the disc component on the galactic plane a Fourier decomposition method is used, that allows us to express it as a sum of trigonometric terms. Both even and odd components are considered, so that the estimated potential accounts also for the observed asymmetries in the morphology. For the amplitudes of the trigonometric terms a smoothed cubic interpolation scheme is used. The total potential in each model may include two additional terms (Plummer spheres) representing a central mass concentration and a dark halo component, respectively. In all examined models, the relative force perturbation points to a strongly nonlinear gravitational field, which ranges from 0.45 to 0.8 of the axisymmetric background with the pure 2D being the most nonlinear one. We present the topological distributions of the stable and unstable Lagrangian points as a function of the pattern speed $(\Omega_p)$. In all three models there is a range of $\Omega_p$ values, where we find multiple stationary points whose stability affects the overall dynamics of the system.

We present the orbital analysis of four response models, that succeed in reproducing morphological features of NGC 1300. Two of them assume a planar (2D) geometry with $\Omega_p$=22 and 16 \ksk respectively. The two others assume a cylindrical (thick) disc and rotate with the same pattern speeds as the 2D models. These response models reproduce most successfully main morphological features of NGC 1300 among a large number of models, as became evident in a previous study. Our main result is the discovery of three new dynamical mechanisms that can support structures in a barred-spiral grand design system. These mechanisms are presented in characteristic cases, where these dynamical phenomena take place. They refer firstly to the support of a strong bar, of ansae type, almost solely by chaotic orbits, then to the support of spirals by chaotic orbits that for a certain number of pat tern revolutions follow an n:1 (n=7,8) morphology, and finally to the support of spiral arms by a combination of orbits trapped around L$_{4,5}$ and sticky chaotic orbits with the same Jacobi constant. We have encountered these dynamical phenomena in a large fraction of the cases we studied as we varied the parameters of our general models, without forcing in some way their appearance. This suggests that they could be responsible for the observed morphologies of many barred-spiral galaxies. Comparing our response models among themselves we find that the NGC 130 0 morphology is best described by a thick disc model for the bar region and a 2D disc model for the spirals, with both components rotating with the same pattern speed $\Omega_p$=16 \ksk !. In such a case, the whole structure is included inside the corotation of the system. The bar is supported mainly by regular orbits, while the spirals are supported by chaotic orbits.

We study the stellar response in a spectrum of potentials describing the barred spiral galaxy NGC 1300. These potentials have been presented in a previous paper and correspond to three different assumptions as regards the geometry of the galaxy. For each potential we consider a wide range of $\Omega_p$ pattern speed values. Our goal is to discover the geometries and the $\Omega_p$ supporting specific morphological features of NGC 1300. For this purpose we use the method of response models. In order to compare the images of NGC 1300 with the density maps of our models, we define a new index which is a generalization of the Hausdorff distance. This index helps us to find out quantitatively which cases reproduce specific features of NGC 1300 in an objective way. Furthermore, we construct alternative models following a Schwarzschild type technique. By this method we vary the weights of the various energy levels, and thus the orbital contribution of each energy, in order to minimize the differences between the response density and that deduced from the surface density of the galaxy, under certain assumptions. We find that the models corresponding to $\Omega_p\approx16$\ksk and $\Omega_p\approx22$\ksk are able to reproduce efficiently certain morphological features of NGC 1300, with each one having its advantages and drawbacks.

The pattern speeds of spiral galaxies are closely related to the flow of material in their disks. Flows that follow the `precessing ellipses' paradigm (see e.g., Kalnajs 1973) are likely associated with slowly rotating spirals, which have corotation beyond their end. Such a flow can be secured by material trapped around stable, elliptical, x_1 periodic orbits precessing as their Jacobi constant varies. Contrarily, if part of the spiral arms is located at a corotation region then the spiral structure has to `survive' in chaotic regions. Barred-spiral systems with a single pattern speed and a bar ending before, but close to, corotation are candidates for having spirals supported by stars in chaotic motion. In this work we review the flows we have found in response models for various types of spiral potentials and indicate the cases, where order or chaos shapes the observed morphologies.

The vertical profiles of disc galaxies are built by the material trapped around stable periodic orbits, which form their "skeletons". According to this, the knowledge of the stability of the main families of periodic orbits in appropriate 3D models, can predict possible morphologies for edge-on disc galaxies. In a pilot survey we compare the orbital structures which lead to the appearance of "peanuts" and "X"-like features with the edge-on profiles of three disc galaxies (IC 2531, NGC 4013 and UGC 2048). The subtraction from the images of a model representing the axisymmetric component of the galaxies reveals the contribution of the non-axisymmetric terms. We find a direct correspondence between the orbital profiles of 3D bars in models and the observed main morphological features of the residuals. We also apply a simple unsharp masking technique in order to study the sharpest features of the images. Our basic conclusion is that the morphology of the boxy "bulges" of these galaxies can be explained by considering disc material trapped around stable 3D periodic orbits. In most models these building-blocks periodic orbits are bifurcated from the planar central family of a non-axisymmetric component, usually a bar, at low order vertical resonances. In such a case the boxy "bulges" are parts of bars seen edge-on. For the three galaxies we study the families associated with the "peanut" or "X"-shape morphology are most probably bifurcations at the vertical 2/1 or 4/1 resonance.

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 argue that inner rings in barred spiral galaxies are associated with specific 2D and 3D families of periodic orbits located just beyond the end of the bar. These are families located between the inner radial ultraharmonic 4:1 resonance and corotation. They are found in the upper part of a type-2 gap of the x1 characteristic, and can account for the observed ring morphologies without any help from families of the x1-tree. Due to the evolution of the stability of all these families, the ring shapes that are favored are mainly ovals, as well as polygons with `corners' on the minor axis, on the sides of the bar. On the other hand pentagonal rings, or rings of the NGC 7020 type hexagon, should be less probable. The orbits that make the rings belong in their vast majority to 3D families of periodic orbits and orbits trapped around them.

We study the conditions that favour boxiness of isodensities in the face-on views of orbital 3D models for barred galaxies. Using orbital weighted profiles we show that boxiness is in general a composite effect that appears when one considers stable orbits belonging to several families of periodic orbits. 3D orbits that are introduced due to vertical instabilities, play a crucial role in the face-on profiles and enhance their rectangularity. This happens because at the 4:1 radial resonance region we have several orbits with boxy face-on projections, instead of few rectangular-like x1 orbits, which, in a fair fraction of the models studied so far, are unstable at this region. Massive bars are characterized by rectangular-like orbits. However, we find that it is the pattern speed that affects most the elongation of the boxy feature, in the sense that fast bars are more elongated than slow ones. Boxiness in intermediate distances between the center of the model and the end of the bar can be attributed to x1v1 orbits, or to a combination of families related to the radial 3:1 resonance.

Boxy edge-on profiles can be accounted for not only in models of barred galaxies, but also in models of normal (non-barred) galaxies. Thus, the presence of a bar is not a sine qua non condition for the appearance of this feature, as often assumed. We show that a `boxy' or a `peanut' structure in the central parts of a model is due to the presence of vertical resonances at which stable families of periodic orbits bifurcate from the planar x1 family. The orbits of these families reach in their projections on the equatorial plane a maximum distance from the center, beyond which they increase their mean radii by increasing only their deviations from the equatorial plane. The resulting orbital profiles are `stair-type' and constitute the backbone for the observed boxy structures in edge-on views of $N$-body models and, we believe, in edge-on views of disc galaxies. Since the existence of vertical resonances is independent of barred or spiral perturbations in the disc, `boxy' profiles may appear also in almost axisymmetric cases.

We investigate the orbital structure in a class of 3D models of barred galaxies. We consider different values of the pattern speed, of the strength of the bar and of the parameters of the central bulge of the galactic model. The morphology of the stable orbits in the bar region is associated with the degree of folding of the x1-characteristic. This folding is larger for lower values of the pattern speed. The elongation of rectangular-like orbits belonging to x1 and to x1-originated families depends mainly on the pattern speed. The detailed investigation of the trees of bifurcating families in the various models shows that major building blocks of 3D bars can be supplied by families initially introduced as unstable in the system, but becoming stable at another energy interval. In some models without radial and vertical 2:1 resonances we find, except for the x1 and x1-originated families, also families related to the z-axis orbits, which support the bar. Bifurcations of the x2 family can build a secondary 3D bar along the minor axis of the main bar. This is favoured in the slow rotating bar case.

In this series of papers we investigate the orbital structure of 3D models representing barred galaxies. In the present introductory paper we use a fiducial case to describe all families of periodic orbits that may play a role in the morphology of three-dimensional bars. We show that, in a 3D bar, the backbone of the orbital structure is not just the x1 family, as in 2D models, but a tree of 2D and 3D families bifurcating from x1. Besides the main tree we have also found another group of families of lesser importance around the radial 3:1 resonance. The families of this group bifurcate from x1 and influence the dynamics of the system only locally. We also find that 3D orbits elongated along the bar minor axis can be formed by bifurcations of the planar x2 family. They can support 3D bar-like structures along the minor axis of the main bar. Banana-like orbits around the stable Lagrangian points build a forest of 2D and 3D families as well. The importance of the 3D x1-tree families at the outer parts of the bar depends critically on whether they are introduced in the system as bifurcations in $z$ or in $\dot{z}$.