Hilfe berechnung reaktionsraten kernfusionsprozessen innerhalb sternen

Gleichung für die Reaktionsrate: R= n1 * n2 * vr * σ n1 und n2 sind die Dichten der Stoßpartner vr die Relativgeschwindigkeit σ der Wirkungsquerschnitt, in Abhängigkeit von vr Wie muss ich vorgehen, wenn ich daraus die Reaktionsrate des p-p-Zyklus z.B. innerhalb der Sonne berechnen möchte. Wo finde ich geeignete Zahlenwerte, um diese Rechnung durchzuführen? Letztendlich möchte ich also berechnen können wie viele Fusionsprozesse in einem bestimmten Stern pro Sekunde stattfinden. Konnte selbst leider kein geeignetes Material finden.