Comparing the subdivion modifier with the subdivisions of Geometry Nodes I found that the smoothing was better when using the nodes.
Checking the Advanced values of the modifier I have seen that Use Limit Surface affects the normals when activated although the geometry remains the same. I understand the limit improves performance, but in this case it seems more of a smoothing issue.
Should I report it as a bug?