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

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.

A model of a mixed helium-hydrogenic atmosphere of the isolated neutron star 1E1207.4-5209 is proposed. It is based on the assumption that the main components in the atmosphere are the exotic molecular ions $He_2^{3+}$ and $H_3^{2+}$ with the presence of $He^{+}, (HeH)^{2+}, H_2^{+}$ under a surface magnetic field $\approx 4.4 \times 10^{13}$ G. In addition to two absorption features observed by {\it Chandra} observatory (Sanwal et al, 2002) the model predicts one more narrow absorption feature at $\approx 400$ eV.

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

In the framework of a variational method with a single trial function an accurate study of the lowest gerade 1g and ungerade 1u electronic states of the molecular ion $H_2^+$ in a magnetic field is performed. Magnetic field ranges from 0 to 4.414x10^13 G and orientations of the molecular axis with respect to the magnetic line 0 \leq theta \leq 90 are considered. A one-parameter gauge dependent vector potential is used in the Hamiltonian, which is finally variationally optimized. A well pronounced minimum on the total energy surface of the (ppe) system in both 1g and 1u states is found for all magnetic fields and orientations studied. It is shown that for both states the parallel configuration (theta=0) at equilibrium always corresponds to the minimal total energy. It is found that for a given theta for both states the magnetic field growth is always accompanied by an increase in the total and binding energies as well as a shrinking of the equilibrium distance. We demonstrate that for B\gtrsim 1.8x10^11 G the molecular ion can dissociate, H2+ -> H + p, over a certain range of orientations (theta_cr \le theta \le 90), where the minimal thetaa_cr \simeq 25 occurs for the strongest magnetic field studied, B=4.414x10^13 G. For B<10^12 G the ion H2+ in 1g, 1u states is the most compact, being in the perpendicular configuration (theta=90), whereas for B\gtrsim 10^12 this occurs for an angle < 90. For the 1g state in any orientation, with the magnetic field growth at B \sim 10^11 G, a two-peak electronic distribution changes to single-peak one.

The quadrupole moment of a hydrogen atom in a magnetic field for field strengths from 0 to 4.414e13 G is calculated by two different methods. The first method is variational, and based on a single trial function. The second method deals with a solution of the Schroedinger equation in the form of a linear combination of Landau orbitals.

In framework of a variational method the molecular ion $H_2^+$ in a magnetic field is studied. An optimal form of the vector potential corresponding to a given magnetic field (gauge fixing) is chosen variationally. It is shown that for any magnetic field strength as well as for any orientation of the molecular axis the system $(ppe)$ possesses a minimum in the potential energy. The stable configuration always corresponds to elongation along the magnetic line. However, for magnetic fields $B \gtrsim 5\times 10^{11} G$ and some orientations the ion $H_2^+$ becomes unstable decaying to $H{-atom} + p$.

Using a variational method we study a sequence of the one-electron atomic and molecular-type systems H, H_2^+, H_3^(2+) and H_4^(3+) in the presence of a homogeneous magnetic field ranging B = 0 - 4.414x10^{13} G. These systems are taken as a linear configuration aligned with the magnetic lines. For H_3^(2+) the potential energy surface has a minimum for B\sim 10^{11} G which deepens with growth of the magnetic field strength (JETP Lett. 69, 844 (1999)); for B \gtrsim 10^{12} G the minimum of the potential energy surface becomes sufficiently deep to have longitudinal vibrational state. We demonstrate that for the (ppppe) system the potential energy surface at B \gtrsim 4.414x10^{13} G develops a minimum, indicating the possible existence of exotic molecular ion H_4^(3+). We find that for almost all accessible magnetic fields H_2^+ is the most bound one-electron linear system while for magnetic fields B \gtrsim 10^{13} G the molecular ion H_3^(2+) becomes the most bound.

Using the variational method it is shown that for magnetic fields $B\geq 10^{11}$ G there can exist a molecular ion $H_3^{++}$.

Using a single trial function we perform an accurate calculation of the ground state $1\sigma_g$ of the hydrogenic molecular ion $H^+_2$ in a constant uniform magnetic field ranging $0-10^{13}$ G. We show that this trial function also makes it possible to study the negative parity ground state $1\sigma_u$. We obtain that over the whole range of magnetic fields studied, the calculated binding energies are in most cases larger than binding energies obtained previously by other authors using different methods.