mini|rechts|Exemplarische Lichtkurve des Asteroiden (250) Bettina vom 10. September 1980 mini|rechts|Durch konvexe Inversion berechnetes 3D-Modell von (250) Bettina Konvexe Inversion, auch Lichtkurveninversion genannt, ist ein Verfahren, um aus photometrischen Daten, wie Lichtkurven, physikalische Modelle von atmosphärelosen Körpern im Sonnensystem, wie Monden oder Asteroiden, zu berechnen, d. h. ihren Rotationszustand und dreidimensionale Gestalt.
Vorgeschichte Bereits kurz nach der Entdeckung der ersten Asteroiden zu Beginn des 19. Jahrhunderts wurde an ihnen eine veränderliche Helligkeit festgestellt. Bei diesen Beobachtungen wurden aber zunächst Effekte beobachtet, die auf einem veränderten Phasenwinkel beruhten. Dabei verändert sich im Laufe einiger Tage (vom Asteroiden aus gesehen) der Winkel zwischen Sonne und Erde und damit auch der Rückstrahlwinkel des zur Erde reflektierten Sonnenlichts und (von der Erde aus gesehen) die Phasengestalt des Asteroiden (vergleichbar mit den Mondphasen). Erst gegen Ende des 19. Jahrhunderts wurden systematisch photometrische Beobachtungen an Asteroiden durchgeführt, bei denen z. B. festgestellt wurde, dass der Phasenkoeffizient (die Helligkeitsänderung pro 1° Veränderung des Phasenwinkels) bei verschiedenen Asteroiden unterschiedlich groß ist.
Zusätzlich wurden aber auch kurzfristige Helligkeitsveränderungen beobachtet. Im Jahr 1899 wurde erstmals beim Asteroiden (345) Tercidina eine Periodizität dieser Veränderungen im Bereich von 3,8 Stunden abgeleitet, im Jahr 1901 konnte Egon von Oppolzer bei (433) Eros sogar eine enorme Änderung der Scheinbaren Helligkeit von fast 2 mag innerhalb von 2,6 Stunden feststellen. Genauere Untersuchungen zeigten geringe Unterschiede in den aufeinanderfolgenden Wellen, so dass schließlich auf eine vollständige Periode von 5,2704 h geschlossen wurde. Bei den Beobachtungen in den folgenden Jahren konnte aber auch bemerkt werden, dass die Amplitude dieser Helligkeitsschwankungen in manchen Jahren wesentlich kleiner und in manchen Jahren nahezu völlig verschwunden war.
Als wahrscheinlichste Ursache für diese regelmäßigen Helligkeitsveränderungen wurde zunächst eine Rotation des Asteroiden und eine gefleckte Oberfläche mit Bereichen stark unterschiedlicher Albedo vermutet. Als Grund für die unterschiedlichen Amplituden der Schwankungen wurden inhomogene Verteilungen dieser Flecken angenommen, die mal mehr und mal weniger von der Erde aus zu sehen wären. Es erwies sich aber, dass diese Theorie ebenso wie die Vermutung einer Bedeckungsveränderlichkeit durch einen umlaufenden Satelliten (vergleichbar mit bedeckungsveränderlichen Sternen) nicht hinreichend ist, um so starke Helligkeitsänderungen, wie sie bei (433) Eros beobachtet wurden, alleine zu erklären.M. Harwood: Variations in the Light of Asteroids. In: Harvard College Observatory Circular. Nr. 269, 1924, S. 1–15, (PDF; 490 kB).
Obwohl auch solche Effekte in besonderen Fällen eine Rolle spielen können, erwies es sich jedoch, dass die Helligkeitsschwankungen in den meisten Fällen auf eine rotierende, unregelmäßige Gestalt der Asteroiden in Verbindung mit einer zur Ebene der Ekliptik geneigten Rotationsachse zurückzuführen sind. Dadurch ändert sich der Querschnitt des sichtbaren und beleuchteten Teils ihrer Oberfläche und somit auch ihre Helligkeit mit der Zeit. Diese Helligkeitsänderung kann als unregelmäßige, aber periodische Lichtkurve registriert werden, aus der sich dann oft direkt die Rotationsperiode des Asteroiden bestimmen lässt. Der Rückschluss auf die verursachende Gestalt des Asteroiden ist aber ungleich schwieriger: Ein nahezu kugelförmiger Asteroid wäre beispielsweise konstant hell, während ein länglicher Asteroid, der um die kurze Achse rotiert, bei Betrachtung senkrecht zur Rotationsachse große Helligkeitsschwankungen und bei Betrachtung aus Polrichtung nur kleine Schwankungen aufweisen würde. Versuche, aus der Form und Veränderlichkeit der Lichtkurve eines Asteroiden Rückschlüsse auf seine Gestalt zu ziehen, führten für (433) Eros z. B. bereits 1938 zur Ableitung der Achsenverhältnisse eines dreiachsig-ellipsoidischen GestaltmodellsF. E. Roach, L. G. Stoddard: A Photoelectric Light-Curve of Eros. In: The Astrophysical Journal. Band 88, 1938, S. 305–312, doi:10.1086/143984 (PDF; 331 kB). und 1952 zur Bestimmung der Lage der Rotationsachse.M. Beyer: Der Lichtwechsel und die Lage der Rotationsachse des Planeten 433 Eros während der Opposition 1951–52. In: Astronomische Nachrichten. Band 281, Nr. 7, 1952, S. 121–130, doi:10.1002/asna.19522810705 (PDF; 482 kB).
Methode der konvexen Inversion Das Problem der Lichtkurveninversion wurde erstmals im Jahr 1906 adressiert durch den US-amerikanischen Astronomen Henry Norris Russell (1877–1957) an der Princeton University in New Jersey. Er formulierte eine mathematische Theorie für das Problem und kam zum Schluss, dass durch die Beobachtung eines veränderlichen Asteroiden während der Opposition in allen Abschnitten seiner Umlaufbahn festgestellt werden kann, ob seine Lichtkurve allein durch seine Rotation erklärt werden kann und falls ja, ob der Asteroid eine absorbierende Atmosphäre besitzt, keine konvexe Form aufweist, eine gefleckte Oberfläche hat oder ob diese Annahmen unzutreffend sind. Außerdem sei es immer möglich, die Lage der Rotationsachse (aber nicht des Drehsinns) zu ermitteln. Aber es wäre wegen der unbekannten Helligkeitsverteilung auf der Oberfläche unmöglich, die Form des Asteroiden zu bestimmen.H. N. Russell: On the light variations of asteroids and satellites. In: The Astrophysical Journal. Band 24, Nr. 1, 1906, S. 1–18, doi:10.1086/141361 (PDF; 888 kB). Russells Schlussfolgerungen waren somit ziemlich pessimistisch, was jedoch an der sehr begrenzten Auswahl an Beobachtungsgeometrien und an den etwas restriktiven mathematischen Methoden lag, die er verwendete.
In der Folgezeit gab es nur wenige Arbeiten zu dem Thema, erst in den 1980er Jahren wurden wieder Untersuchungen über die Auswirkungen verschiedener Form- und Oberflächeneigenschaften auf Lichtkurven durchgeführt. Eine Untersuchung von 1992 präsentierte dann eine Theorie und Verfahren zur Lichtkurveninversion, die zur Bestimmung der dreidimensionalen Form und/oder der Albedoverteilung auf der Oberfläche eines Körpers aus photometrischen Daten verwendet werden konnten – unter der Annahme einer streng konvexen Form des Körpers.M. Kaasalainen, L. Lamberg, K. Lumme, E. Bowell: Interpretation of lightcurves of atmosphereless bodies. I. General theory and new inversion schemes. In: Astronomy & Astrophysics. Band 259, Nr. 1, 1992, S. 318–332, (PDF; 1,43 MB).M. Kaasalainen, L. Lamberg, K. Lumme: Interpretation of lightcurves of atmosphereless bodies. II. Practical aspects of inversion. In: Astronomy & Astrophysics. Band 259, Nr. 1, 1992, S. 333–340, (PDF; 740 kB). Die vollständige und praktikable Methode der Lichtkurveninversion (LI) wurde dann bis 2001 von Mikko Kaasalainen (1965–2020), Johanna Torppa und Karri Muinonen (* 1961) am Observatorium und Astrophysikalischen Labor der Universität Helsinki entwickelt.M. Kaasalainen, J. Torppa: Optimization Methods for Asteroid Lightcurve Inversion. I. Shape Determination. In: Icarus. Band 153, Nr. 1, 2001, S. 24–36, doi:10.1006/icar.2001.6673.M. Kaasalainen, J. Torppa, K. Muinonen: Optimization Methods for Asteroid Lightcurve Inversion. II. The Complete Inverse Problem. In: Icarus. Band 153, Nr. 1, 2001, S. 37–51, doi:10.1006/icar.2001.6674.
Mathematischer Hintergrund Das Problem der konvexen Inversion kann durch folgende Gleichung beschrieben werden:
,
wobei der Vektor der beobachteten Helligkeiten ist, der über die Matrix mit dem Vektor verknüpft ist, der die zu lösenden Parameter enthält. Dieser Vektor kann die Krümmungsfunktion des Objekts beschreiben, welche die Form eindeutig bestimmt, er kann aber auch die Albedoverteilung oder das Produkt aus beiden darstellen.
Die Inversion umfasst drei Schritte. Zunächst wird eine Funktion (oder mehrere Funktionen) bestimmt, die Informationen sowohl über die Form als auch über die Albedo-Variation eines Asteroiden enthält. Dieser Schritt ist durchführbar, sofern die Oberfläche streng konvex ist und eine ausreichende Anzahl von Lichtkurven bei verschiedenen Beobachtungsgeometrien und Phasenwinkeln ungleich null vorliegt. Außerdem muss die funktionale Form des Oberflächenstreugesetzes bekannt und von geeignetem Typ sein (z. B. Lambertsches Gesetz oder Lommel-Seeliger-Gesetz). Dieses Inversionsproblem ist mathematisch schlecht gestellt, d. h. kleine Fehler in den Daten können große Auswirkungen auf die Ergebnisse haben. Auch die Anzahl und der Bereich der Beobachtungsgeometrien haben einen signifikanten Einfluss auf die Inversion. Die Verwendung positiv definiter Größen kann jedoch die scheinbare Schlechtgestelltheit des Problems effektiv beseitigen.
Im zweiten Schritt werden separate Ausdrücke für den Kehrwert der Gaußschen Krümmung und der Albedo-Verteilung aus den im ersten Schritt gewonnenen Informationen abgeleitet. Dies ist möglich, wenn die funktionale Form des Streugesetzes als Funktion der Albedo bekannt ist und eine geeignete Form aufweist.
Im dritten Schritt wird das nicht-triviale Problem der Bestimmung des Radiusvektors der Oberfläche aus der Gaußschen Krümmung mithilfe von iterativen Optimierungsverfahren gelöst. Dabei ergeben sich auch Informationen über die Rotationsachse und die Rotationsperiode.
Praktische Durchführung Eingangsgrößen für die konvexe Inversion sind photometrische Daten zu dem jeweiligen Asteroiden, die möglichst nicht nur als relative, sondern als genaue scheinbare Helligkeiten vorliegen müssen. Zu jedem Datensatz muss auch genau die Position des Asteroiden auf seiner Umlaufbahn sowie der Phasenwinkel bekannt sein. Quelle für die photometrischen Daten können durch Amateurastronomen oder an professionellen Observatorien durchgeführte Messungen sowie die Auswertung von Archiven oder (auch ältere) wissenschaftliche Veröffentlichungen sein. In neuerer Zeit kommen dafür auch die Daten von Weltraumteleskopen infrage, wie Hubble, Spitzer, Gaia, Transiting Exoplanet Survey Satellite (TESS), Hipparcos und anderen, sowie die veröffentlichten Ergebnisse großer Himmelsdurchmusterungen, wie die Catalina Sky Survey, die Siding Spring Survey, die Mount Lemmon Survey, die Palomar Transient Factory (PTF), die Zwicky Transient Facility (ZTF), das Asteroid Terrestrial-impact Last Alert System (ATLAS), die All-Sky Automated Survey for Supernovae (ASAS-SN), LONEOS und weitere.
Diese Daten können sowohl in „dichter“ (dense) Form, wie z. B. als detaillierte Lichtkurven, aber auch in Form von „verstreuten“ (sparse) Werten mit großem zeitlichen Abstand vorliegen. Ursprünglich wurde nur dichte Photometrie verwendet. Etwa 20 solche dichten Lichtkurven von mindestens vier oder fünf Erscheinungen sind für eine eindeutige Formbestimmung erforderlich. Mit diesem Ansatz waren zunächst etwa 100 Asteroidenmodelle abgeleitet worden. Um die Anzahl der Modelle deutlich zu erhöhen, wurden dann zusätzlich auch verstreute Daten, z. B. aus archivierten Messungen des United States Naval Observatory (USNO) in Arizona verwendet.