We have used the galaxy modeling algorithm \texttt{GravPot16}, to explore the more probable orbital elements of a sample of 64 selected N-rich stars across the Milky Way. Using the newly measured proper motions from \texttt{Gaia} DR2 with existing line-of-sight velocities from APOGEE-2 survey and spectrophotometric distance estimations from the \texttt{StarHorse}. We adopted a set of high-resolution particle simulations evolved in the same steady-state Galactic potential model with a bar, in order to identify the groups of N-rich stars that have a high probability of belonging to the bulge/bar, disk, and stellar halo component. We find that the vast majority of the N-rich stars show typically maximum height from the Galactic plane below 3 kpc, and develop rather eccentric orbits (\textit{e}$>$0.5), which means these stars appear to have bulge/bar-like and/or halo-like orbits. We also show that $\sim66$\% of the selected N-rich stars currently lives in the inner Galaxy inside the corotation radius (C.R.), whilst $\sim14$\% of the N-rich star resides in halo-like orbits. Among the N-rich in the inner Galaxy, $\sim27\%$ of them share orbital properties in the boundary between bulge/bar and disk, depending on the bar pattern speeds. Our dynamical analysis also indicates that some of the N-rich are likely halo interlopers and therefore suggest that halo contamination is not insignificant within the bulge area.

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 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.

We performed a series of 3D N-body simulations where the initial conditions were chosen to get two sets of models; unbarred and barred ones. In this work, we analyze the growth of spirals and bar structures using 1D, and 2D Fourier Transforms FT methods. Spectrograms and diagrams of the amplitude of the Fourier coefficients as a function of time, radius and pitch angle show that the general morphology, of our modeled galaxies, is due to the superposition of structures which have different values of pitch angle and number of arms. Also, we made in barred models a geometric classification of orbits from the bar reference frame showing that the barred potential and the Lagrangian points $L_4$ and $L_5$ catch approximately one-third of the total disk mass.