A non-relativistic classification of charged molecular hydrogenic, helium and mixed helium-hydrogenic chains with one or two electrons which can exist in a strong magnetic field $B \lesssim 10^{16} $G is given. It is shown that for both $1e-2e$ cases at the strongest studied magnetic fields the longest hydrogenic chain contains at most five protons indicating to the existence of the $\rm{H}_5^{4+}$ and $\rm{H}_5^{3+}$ ions, respectively. In the case of the helium chains the longest chains can exist at the strongest studied magnetic fields with three and four $\al-$particles for $1e-2e$ cases, respectively. For mixed helium-hydrogenic chains the number of heavy centers can reach five for highest magnetic fields studied. In general, for a fixed magnetic field two-electron chains are more bound than one-electron ones.

A first detailed study of the ground state of the H$_3^+$ molecular ion in linear configuration, parallel to a magnetic field direction, and its low-lying $\Si,\Pi,\De$ states is carried out for magnetic fields $B=0-4.414 \times 10^{13} $G in the Born-Oppenheimer approximation. The variational method is employed with a single trial function which includes electronic correlation in the form $\exp{(\ga r_{12})}$, where $\ga$ is a variational parameter. It is shown that the quantum numbers of the state of the lowest total energy (ground state) depend on the magnetic field strength. The ground state evolves from the spin-singlet ${}^1\Si_g$ state for weak magnetic fields $B \lesssim 5 \times 10^{8} $G to a weakly-bound spin-triplet ${}^3\Si_u$ state for intermediate fields and, eventually, to a spin-triplet $^3\Pi_u$ state for $5 \times 10^{10} \lesssim B \lesssim 4.414 \times 10^{13} $G. Local stability of the linear parallel configuration with respect to possible small deviations is checked.

A detailed study of the ground state of the Coulomb system $(\al \al e e)$ which corresponds to the $He_2^{2+}$ molecular ion in a magnetic field $B=0-4.414 \times 10^{13} $ G in parallel configuration (infinitely massive $\al-$particles are situated along a magnetic field line) is presented. The variational method is employed using a trial function which includes electronic correlation in the form $\exp{(\gamma r_{12})}$ where $\gamma$ is a variational parameter. It is shown that the quantum numbers of the lowest total energy state depend on the magnetic field strength. It evolves from the spin-singlet ${}^1\Si_g$ metastable state at $0 \leq B \lesssim 0.85$ a.u. to a repulsive spin-triplet ${}^3\Si_u$ state for $0.85 {a.u.} \lesssim B \lesssim 1100$ a.u. and, finally, to a strongly bound spin-triplet ${}^3\Pi_u$ state at stronger fields $B \gtrsim 1100$ a.u. The lowest vibrational energy for the latter case is calculated.

An accurate study of the lowest $1\sigma_g$ and the low-lying excited $1\sigma_u$, $1\pi_{u,g}$, $1\delta_{g,u}$ electronic states of the exotic molecular ion $H_3^{2+}$ in linear configuration parallel to a magnetic field is carried out. The magnetic field ranges from $10^{10}$ G up to $4.414 \times 10^{13}$ G where non-relativistic considerations are justified. The variational method is exploited and the {\it same} trial function is used for different magnetic fields. It is shown that the states of positive $z$-parity $1\sigma_g, 1\pi_u, 1\delta_{g}$ are bound states of the $H_3^{2+}$ exotic ion for all magnetic fields studied. We also demonstrate that for magnetic fields $B\gtrsim 2.35\times 10^{12}$ G the potential energy surface well corresponding to the $1\sigma_g$ state contains at least one longitudinal vibrational state. It is also shown that the negative $z$-parity states $1\sigma_u, 1\pi_g, 1\delta_{u}$, are purely repulsive in the whole range of magnetic fields studied, $B=10^{10}- 4.414 \times 10^{13}$ G.