A comparative analysis of the astrophysical S factor and the reaction rate for the direct $ \alpha(d,\gamma)^{6}{\rm Li}$ capture reaction, and the primordial abundance of the $^6$Li element, resulting from two-body, three-body and combined cluster models is presented. It is shown that the two-body model, based on the exact-mass prescription, can not correctly describe the dependence of the isospin-forbidden E1 S factor on energy and does not reproduce the temperature dependence of the reaction rate from the direct LUNA data. It is demonstrated that the isospin-forbidden E1 astrophysical S factor is very sensitive to the orthogonalization procedure of Pauli-forbidden states within the three-body model. On the other hand, the E2 S factor does not depend on the orthogonalization method. This insures that the orthogonolizing pseudopotentials method yields a very good description of the LUNA collaboration's low-energy direct data. At the same time, the SUSY transformation significantly underestimates the data from the LUNA collaboration. On the other hand, the energy dependence of the E1 S factor are the same in both methods. The best description of the LUNA data for the astrophysical S factor and the reaction rates is obtained within the combined E1(three-body OPP)+E2(two-body) model. It yields a value of $(0.72 \pm 0.01) \times 10^{-14}$ for the $^6$Li/H primordial abundance ratio, consistent with the estimation $(0.80 \pm 0.18) \times 10^{-14}$ of the LUNA collaboration. For the $^6{\rm Li}/^7{\rm Li}$ abundance ratio an estimation $(1.40\pm 0.12)\times 10^{-5}$ is obtained in good agreement with the Standard Model prediction.

At the long-wavelength approximation, electric dipole transitions are forbidden between isospin-zero states. In an $\alpha+n+p$ model with $T = 1$ contributions, the $\alpha(d,\gamma)^6$Li astrophysical $S$-factor is in agreement with the experimental data of the LUNA collaboration, without adjustable parameter. The exact-masses prescription used to avoid the disappearance of $E1$ transitions in potential models is not founded at the microscopic level.

The reaction rates of the direct astrophysical capture processes $^{3}{\rm He}(\alpha, \gamma)^{7}{\rm Be}$ and $^{3}{\rm H}(\alpha,\gamma)^{7}{\rm Li}$, as well as the abundance of the $^{7}{\rm Li}$ element are estimated in the framework of a two-body potential model. The estimated $^{7}{\rm Li/H}$ abundance ratio of $^{7}{\rm Li/H}=(5.07\pm 0.14 )\times 10^{-10}$ is in a very good agreement with the recent measurement $^{7}{\rm Li/H}=(5.0\pm 0.3) \times 10^{-10}$ of the LUNA collaboration.

The astrophysical S-factor and reaction rate of the direct capture process $\alpha+d$ $\rightarrow$ $^6$Li + $\gamma$, as well as the abundance of the $^6$Li element are estimated in a three-body model. The initial state is factorized into the deuteron bound state and the $\alpha+d$ scattering state. The final nucleus $^6$Li(1+) is described as a three-body bound state $\alpha+n+p$ in the hyperspherical Lagrange-mesh method. Corrections to the asymptotics of the overlap integral in the S- and D-waves have been done for the E2 S-factor. The isospin forbidden E1 S-factor is calculated from the initial isosinglet states to the small isotriplet components of the final $^6$Li(1+) bound state. It is shown that the three-body model is able to reproduce the newest experimental data of the LUNA collaboration for the astrophysical S-factor and the reaction rates within the experimental error bars. The estimated $^6$Li/H abundance ratio of $(0.67 \pm 0.01)\times 10^{-14}$ is in a very good agreement with the recent measurement $(0.80 \pm 0.18)\times 10^{-14}$ of the LUNA collaboration.

The astrophysical $^{3}{\rm He}(\alpha, \gamma)^{7}{\rm Be}$ and $^{3}{\rm H}(\alpha, \gamma)^{7}{\rm Li}$ direct capture processes are studied in the framework of the two-body model with the potentials of a simple Gaussian form, which describe correctly the phase-shifts in the s-, p-, d-, and f-waves, as well as the binding energy and the asymptotic normalization constant of the ground $p_{3/2}$ and the first excited $p_{1/2}$ bound states. It is shown that the E1-transition from the initial s-wave to the final p-waves is strongly dominant in both capture reactions. On this basis the s-wave potential parameters are adjusted to reproduce the new data of the LUNA collaboration around 100 keV and the newest data at the Gamov peak estimated with the help of the observed neutrino fluxes from the Sun, $S_{34}$(23$^{+6}_{-5}$ keV)=0.548$\pm$0.054 keV b for the astrophysical S-factor of the capture process $^{3}{\rm He}(\alpha, \gamma)^{7}{\rm Be}$. The resulting model describes well the astrophysical S-factor in low-energy Big Bang nucleosynthesis region of 180-400 keV, however has a tendency to underestimate the data above 0.5 MeV. Two-body potentials, adjusted on the properties of the $^7$Be nucleus, $^3{\rm He}+\alpha$ elastic scattering data and the astrophysical S-factor of the $^{3}{\rm He}(\alpha, \gamma)^{7}{\rm Be}$ direct capture reaction, are able to reproduce the properties of the $^7$Li nucleus, the binding energies of the ground 3/2$^-$ and first excited 1/2$^-$ states, and phase shifts of the $^3 {\rm H}+\alpha$ elastic scattering in partial waves. Most importantly, these potential models can successfully describe both absolute value and energy dependence of the existing experimental data for the mirror astrophysical $^{3}{\rm H}(\alpha, \gamma)^{7}{\rm Li}$ capture reaction without any additional adjustment of the parameters.

The astrophysical capture process $\alpha+d$ $\rightarrow$ $^6$Li + $\gamma$ is studied in a three-body model. The initial state is factorized into the deuteron bound state and the $\alpha+d$ scattering state. The final nucleus $^6$Li(1+) is described as a three-body bound state $\alpha+n+p$ in the hyperspherical Lagrange-mesh method. The contribution of the E1 transition operator from the initial isosinglet states to the isotriplet components of the final state is estimated to be negligible. An estimation of the forbidden E1 transition to the isosinglet components of the final state is comparable with the corresponding results of the two-body model. However, the contribution of the E2 transition operator is found to be much smaller than the corresponding estimations of the two-body model. The three-body model perfectly matches the new experimental data of the LUNA collaboration with the spectroscopic factor 2.586 estimated from the bound-state wave functions of $^6$Li and deuteron.