Accurate calculations for the ground state of the molecular ions He$_2^{3+}$ and HeH$^{2+}$ placed in a strong magnetic field $B\gtrsim 10^{2}$ a.u. ($\approx 2.35 \times 10^{11}$G) using the Lagrange-mesh method are presented. The Born-Oppenheimer approximation of zero order (infinitely massive centers) and the parallel configuration (molecular axis parallel to the magnetic field) are considered. Total energies are found with 9-10 s.d. The obtained results show that the molecular ions He$_2^{3+}$ and HeH$^{2+}$ exist at $B > 100$\,a.u. and $B > 1000$\,a.u., respectively, as predicted in \cite{Tu:2007} while a saddle point in the potential curve appears for the first time at $B \sim 80$ a.u. and $B \sim 740$ a.u., respectively.

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