I'm trying to simulate the Kepler-444 exoplanetary system using data from this website. The problem is that there is no argument of periapsis parameter for the Kepler-444 system, so while I think I'm getting the correct shape of the orbits, their periapsises and apoapsis are all aligned in the same direction, as you can see in the screen dump below.
So my question is how, or if, I can get the periapsis to be correctly aligned with the data that is available on this website for these planets (unfortunately I can't link to a query, but you can search Kepler-444 and the planets should pop up in the table)? The way I'm doing it now is that I take the semimajor axis and the eccentricity data from the website I linked to, and then I calculate the periapsis for each planet and then I set the x axis to be that value while I let y and z be zero. Once there, I use the vis-viva equation to calculate the velocity vectors for the planet, and then I iterate the simulation forwards in time and get the shapes of the orbits. With the solar system planets, I take the argument of periapsis and the inclination of the orbit and rotate the position and velocity vectors that I derived accordingly, and this way I end up with orbits that have the right shape and their periapsises are correctly aligned, too. Naturally I don't get the right position of a planet in its orbit at a given time using this approach, but for now that is not important; I just want the orbits to have the correct shape and to have their periapsises correctly aligned.
I know very little about orbital mechanics, so any help would be much appreciated so that I can learn something and improve the quality of my simulation.