Neutron-deuteron scattering revisited with the EKM chiral nuclear force and the WPCD method

We revisit the neutron-deuteron scattering using the Wave-Packet Continuum Discretization (WPCD) method with the EKM chiral nuclear force at various chiral orders. We rederive the permutation operator and solve the Faddeev-AGS equations directly, without rewriting the initial Faddeev kernel $tG_0$ and introducing pseudo-states, thereby rendering the approach easily extendable to a relativistic framework. We find that up to the next-to-next-to-next-to-leading order (N$^3$LO), although one can well describe the differential cross sections, one cannot resolve the long-standing $A_y$ puzzle, consistent with previous studies. The fact that the N$^3$LO chiral forces can well describe the $NN$ phase shifts and the results obtained with the EKM and Idaho N$^3$LO chiral forces agree with each other underscores the need for further investigations to resolve the $A_y$ puzzle, e.g., considering three-body forces or relativistic effects.