This study investigates the convection flow of nano-encapsulated phase change material (NEPCM)-water mixture in an evacuated tube solar collector manifold, focusing on the effects of magnetohydrodynamic (MHD) double-diffusive convection and exothermic reaction. Computational fluid dynamics simulations using the Galerkin finite element method were employed to analyze temperature and concentration distributions, fluid flow, melting zone of PCM nanocapsules, Nusselt number, Sherwood number, entropy generation, and thermal performance. The numerical analysis considers a range of dimensionless parameters, including Rayleigh number (Ra=103−105), Lewis number (Le=0.1−10), buoyancy ratio (Nz=1−5), nanoparticle concentration (ϕ=0.01−0.035), fusion temperature (ϴf=0.1−0.9), Stefan number (Ste=0.1−0.9), Hartmann number (Ha=0−50), magnetic field inclination angle (γ=0°−90°), and Frank-Kamenetskii number (FK=0−2.5). Results show that increasing Ra enhances the average Nusselt number (Nuav) by up to 156.1 % and the average Sherwood number (Shav) by up to 104.5 %, while increasing the total entropy generation by up to 13,155 %. Increasing FK reduces Nuav by up to 42.2 % but increases Shav by up to 13.1 %. The Lewis number and buoyancy ratio significantly influence the hydrothermal performance, with Nuav exhibiting non-monotonic behavior and Shav increasing by up to 240.9 % as Le increases. The nanoparticle concentration enhances Nuav by up to 49.7 %. Interestingly, the magnetic field effects are less pronounced than in previous studies, with Ha having a minor impact on Nuav and Shav. These findings have significant implications for the design and optimization of evacuated tube solar collectors and other thermal energy storage systems, potentially leading to improved efficiency and performance in real-world applications.