I am trying to implement an algorithm I created to get the sun rise, transit and set accurately, but in order to do that I need to get the Earth-Moon barycenter (EMB) perihelion, and not the Earth perihelion. I tried to get such time using the Horizons System, but it resulted in weird outputs (it showed that the perihelion is around December 1st, 2021, but it should be closer to the Earth perihelion at January 4th, 2022). Are there alternatives to get the EMB perihelion time?
From a comment by the OP,
How do I set it to the center of the Sun?
Select Vector table as the ephemeris type. Choose the Earth-Moon barycenter as the target body and @sun as the coordinate center. The only value you're interested in in the output is the range rate. Look for when it changes from negative to positive if you're looking for perihelion passage time, or from positive to negative if you're looking for aphelion passage time.