The existence and stability of the linear hydrogenic chain H$_3$ and H${}_2^-$ in a strong magnetic field is established. Variational calculations for H$_3$ and H${}_2^-$ are carried out in magnetic fields in the range $10^{11}\leq B \leq 10^{13}\,$G with 17-parametric (13-parametric for H${}_2^-$), physically adequate trial function. Protons are assumed infinitely massive, fixed along the magnetic line. States with total spin projection $S_z=-3/2$ and magnetic quantum numbers $M=-3,-4,-5$ are studied. It is shown that for both H$_3$ and H${}_2^-$ the lowest energy state corresponds to $M=-3$ in the whole range of magnetic fields studied. As for a magnetic field $B \gtrsim 10^{11}\,$G both H$_3$ and H${}_2^-$ exist as metastable states, becoming stable for $B \geq 1.9 \times 10^{11}\,$G and for $B \geq 2.7 \times 10^{11}\,$G, respectively. The excited states $^4(-4)^+$, $^4(-5)^+$ of ${\rm H}_3$ and H${}_2^-$ appear at magnetic fields $B > 7 \times 10^{11}$ and $10^{12}$\,G, respectively.

The existence of bound states of the system $(\al,e,e,e)$ in a magnetic field $B$ is studied using the variational method. It is shown that for $B \gtrsim 0.13\,{\rm a.u.}$ this system gets bound with total energy below the one of the $(\al,e,e)$ system. It manifests the existence of the stable He$^-$ atomic ion. Its ground state is a spin-doublet $^2(-1)^{+}$ at $0.74\, {\rm a.u.} \gtrsim B \gtrsim 0.13\, {\rm a.u.}$ and it becomes a spin-quartet $^4(-3)^{+}$ for larger magnetic fields. For $0.8\, {\rm a.u.} \gtrsim B \gtrsim 0.7\, {\rm a.u.}$ the He$^-$ ion has two (stable) bound states $^2(-1)^{+}$ and $^4(-3)^{+}$.

Three simple $7-, (7+3)-, 10-$parametric trial functions for the ${\rm H}_3^+$ molecular ion are presented. Each of them provides subsequently the most accurate approximation for the Born-Oppenheimer ground state energy among several-parametric trial functions. These trial functions are chosen following a criterion of physical adequacy and includes the electronic correlation in the exponential form $\sim\exp{(\gamma r_{12})}$, where $\gamma$ is a variational parameter. The Born-Oppenheimer energy is found to be $E=-1.340 34, -1.340 73, -1.341 59$\,a.u., respectively, for optimal equilateral triangular configuration of protons with the equilibrium interproton distance $R=1.65$\,a.u. The variational energy agrees in three significant digits (s.d.) with most accurate results available at present as well as for major expectation values.

The one-electron Li-containing Coulomb systems of atomic type $(li, e)$ and molecular type $(li, li, e)$, $(li, \alpha, e)$ and $(li, p, e)$ are studied in the presence of a strong magnetic field $B \leq 10^{7}$ a.u. in the non-relativistic framework. They are considered at the Born-Oppenheimer approximation of zero order (infinitely massive centers) within the parallel configuration (molecular axis parallel to the magnetic field). The variational and Lagrange-mesh methods are employed in complement to each other. It is demonstrated that the molecular systems ${\rm LiH}^{3+}$, ${\rm LiHe}^{4+}$ and ${\rm Li}_{2}^{5+}$ can exist for sufficiently strong magnetic fields $B \gtrsim 10^{4}$ a.u. and that ${\rm Li}_{2}^{5+}$ can even be stable at magnetic fields typical of magnetars.

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 detailed study of the lowest states $1s_0, 2p_{-1}, 2p_0$ of the hydrogen atom placed in a magnetic field $B\in(0-4.414\times 10^{13} {\rm G})$ and their electromagnetic transitions ($1s_{0} \leftrightarrow 2p_{-1}$ and $ 1s_{0} \leftrightarrow 2p_{0}$) is carried out in the Born Oppenheimer approximation. The variational method is used with a physically motivated recipe to design simple trial functions applicable to the whole domain of magnetic fields. We show that the proposed functions yield very accurate results for the ionization (binding) energies. Dipole and oscillator strengths are in good agreement with results by Ruder {\em et al.} \cite{Ruderbook} although we observe deviations up to $\sim 30%$ for the oscillator strength of the (linearly polarized) electromagnetic transition $1s_{0} \leftrightarrow 2p_{0}$ at strong magnetic fields $B\gtrsim 1000$ a.u.

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.

The Coulombic systems $(\al p e)$ and $(\al\al e)$, $(\al p p e)$, $(\al \al p e)$ and $(Li^{3+} Li^{3+} e)$ placed in a magnetic field $B \gtrsim 10^{11} {G}$ are studied. It is demonstrated a theoretical existence of the exotic ion $(He H)^{2+}$ for $B\gtrsim 5\times 10^{12} {G}$ in parallel configuration (the magnetic field is directed along internuclear axis) as optimal as well as its excited states $1\pi, 1\delta$. As for the exotic ion ${He}_{2}^{3+}$ it is shown that in spite of strong electrostatic repulsion of $\al-$particles this ion can also exist for $B \gtrsim 100 {a.u.} (= 2.35\times 10^{11} {G})$ in parallel configuration as optimal in the states $1\si_g$ (ground state), $1\pi_u, 1\delta_g$. Upon appearance both ions are unstable towards dissociation with $He^+$ in the final state but with very large lifetime. However, at $B\gtrsim 10000$ a.u. the ion $(He H)^{2+}$ becomes stable, while at $B\gtrsim 1000$ a.u. the ion ${He}_{2}^{3+}$ becomes stable. With a magnetic field growth, both exotic ions become more and more tightly bound and compact, their lowest rotational and vibrational energies grow. At the edge of applicability of non-relativistic approximation, $B \sim 4.414 \times 10^{13}$ G, there are indications that three more exotic linear ions $(H-He-H)^{3+}$, $(He-H-He)^{4+}$ and even $Li_2^{5+}$ in parallel configuration may also occur.

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.

A model of a hydrogenic content of atmosphere of the isolated neutron star 1E1207.4-5209 is proposed. It is based on the assumption that the main component in the atmosphere is the exotic molecular ion $H_3^{2+}$ and that there exists a magnetic field in the range of $(4 \pm 2) \times 10^{14}$ G. Photoionization $H_3^{2+} \rar e + 3p$ and photodissociation $H_3^{2+} \rar H + 2p$ correspond to two absorption features at 0.7 KeV and 1.4 KeV, respectively, discovered by {\it Chandra} observatory (Sanwal et al, 2002). The model predicts one more absorption feature at 80-150 eV corresponding to photodissociation $H_3^{2+} \rar H_2^+ + p$.

As a continuation of our previous work ({\it Phys. Rev. A68, 012504 (2003)}) an accurate study of the lowest $1\si_g$ and the low-lying excited $1\si_u$, $2\si_g$, $1\pi_{u,g}$, $1\de_{g,u}$ electronic states of the molecular ion $H_2^+$ is made. Since the parallel configuration where the molecular axis coincides with the magnetic field direction is optimal, this is the only configuration which is considered. The variational method is applied and the {\it same} trial function is used for different magnetic fields. The magnetic field ranges from $10^9 G$ to $4.414 \times 10^{13} G$ where non-relativistic considerations are justified. Particular attention is paid to the $1\si_u$ state which was studied for an arbitrary inclination. For this state a one-parameter vector potential is used which is then variationally optimized.

A detailed quantitative analysis of the system $(ppe)$ placed in magnetic field ranging from $0 - 4.414 \times 10^{13} G$ is presented. The present study is focused on the question of the existence of the molecular ion $H_2^{+}$ in a magnetic field. As a tool, a variational method with an optimization of the form of the vector potential (optimal gauge fixing) is used. It is shown that in the domain of applicability of the non-relativistic approximation the system $(ppe)$ in the Born-Oppenheimer approximation has a well-pronounced minimum in the total energy at a finite interproton distance for $B \lesssim 10^{11} G$, thus manifesting the existence of $H_2^{+}$. For $B \gtrsim 10^{11} G$ and large inclinations (of the molecular axis with respect to the magnetic line) the minimum disappears and hence the molecular ion $H_2^{+}$ does not exist. It is shown that the most stable configuration of $H_2^{+}$ always corresponds to protons situated along the magnetic line. With magnetic field growth the ion $H_2^{+}$ becomes more and more tightly bound and compact, and the electronic distribution evolves from a two-peak to a one-peak pattern. The domain of inclinations where the $H_2^{+}$ ion exists reduces with magnetic field increase and finally becomes $0^o - 25^o$ at $B = 4.414 \times 10^{13} G$. Phase transition type behavior of variational parameters for some interproton distances related to the beginning of the chemical reaction $H_2^{+} \leftrightarrow H + p$ is found.

The existence of the molecular ion $H_3^{++}$ in a magnetic field in a triangular configuration is revised. Variational method with an optimization of the form of the vector potential (gauge fixing) is used. It is shown that in the range of magnetic fields $10^8 < B < 10^{11} G$ the system $(pppe)$, with protons forming an equilateral triangle perpendicular to magnetic line, has a well-pronounced minimum in the total energy. This configuration is unstable towards decays $H{-atom} + p + p$ and $H_2^+ + p$. Triangular configuration of $H_3^{++}$ complements $H_3^{++}$ in the linear configuration which exists for $B \gtrsim 10^{10} G$.