Zjawisko to tłumaczy się wielofotonową absorpcją. Może ono zachodzić, gdy moc monochromatycznego liniowo spolaryzowanego światła osiąga bardzo duże wartości. Oznacza to, że w wyniku superpozycji pól elektromagnetycznych poszczególnych fotonów powstaje wyjątkowo silne pole elektromagnetyczne. Efekt taki można uzyskać dzięki użyciu światła laserowego. Energię tego procesu określa wzór:
gdzie oznacza liczbę fotonów, których energia została przejęta przez pojedynczy elektron. Ze względu na oddziaływanie z jonem macierzystym i zderzenia, widmo elektronów nie ma postaci ostrych pików, lecz jest dość szerokie i stałe (plateau). Energie fotoelektronów powstających w procesach wielofotonowych mogą osiągać wartości porównywalne z wartościami energii β-elektronów, chociaż ich widmo jest zupełnie inne.
Jonizacje powyżej progu można wyjaśnić rozwiązując równanie Schrödingera w sposób przybliżony. Równanie Schrödingera dla elektronu swobodnego w polu fali elektromagnetycznej w jednym wymiarze w cechowaniu promieniowania jest dane przez
Równanie Schrödingera dla elektronu w polu fali i w polu potencjału atomowego będzie dane przez
gdzie jest hamiltonianem elektronu swobodnego. Dodając i odejmując energię stanu podstawowego, z którego będzie jonizowany elektron otrzymujemy równanie
Ponieważ w stanie podstawowym energia kinetyczna elektronu jest równa energii całkowitej z przeciwnym znakiem (twierdzenie o wiriale) i tylko ona zostanie po szybkim usunięciu elektronu, pomijamy w tym równaniu sumę dla wszystkich i otrzymujemy równanie przybliżone
gdzie jedyna pozostałość po potencjale atomowym jest stała.
Równanie to można rozwiązać wykorzystując rozwiązania dla elektronu swobodnego i rozkładając stan podstawowy na składowe Fouriera:
z
Równanie to ma więc rozwiązanie
Widmo jonizacji otrzymujemy ze wzoru
mówiącego ile składowej fali płaskiej elektronu swobodnego o danej energii kinetycznej jest pod koniec procesu jonizacji, gdzie jest funkcją uśredniającą detektora pomiarowego, np.
Rozkładając czynnik
z uogólnionymi funkcjami Bessela zdefiniowanymi przez transformatę odwrotną otrzymujemy
( jest transformatą Fouriera funkcji detektora), czyli sumą ostrych lub rozmytych maksimów zlokalizowanych wokół warunku energii emitowanych elektronów w zależności od szybkości, tzn. od parametru uśredniania detektora