Youcef Kehal, Khireddine Nouicer

We study neutron star configurations in a teleparallel gravity model featuring a scalar field coupled to both matter and torsion. In the Einstein frame, the theory includes a derivative coupling between the scalar field and the torsion vector, together with a conformal matter coupling \(A(φ)=\exp(βφ^{2}/2)\). Static and slowly rotating neutron-star solutions are constructed for realistic equations of state, focusing on the APR and MS1 equations of state. Scalarized solutions appear only within a finite range of central densities and correspond to localized deviations from the general-relativistic mass--radius and mass--central-density relations. The onset and extent of scalarization depend on the equation of state and on the strength of the derivative torsional interaction, which can either enhance or suppress scalarization relative to the general-relativistic scalarized branch. At high central densities, scalarization is quenched and the solutions approach the general-relativistic limit, remaining bounded even for strong torsional couplings. No scalarized solutions are found in the absence of matter coupling (\(β=0\)). The normalized scalar charge follows trends consistent with the global mass relations, indicating an intermediate scalarized regime suppressed at high compactness. For slowly rotating stars, the moment of inertia depends systematically on the torsional coupling and the equation of state, with stiffer equations yielding larger values. These results highlight the potential of neutron-star radius and rotational measurements to test teleparallel scalarization scenarios.