Een reflectiecaustiek is de heldere kromme die een gekromde spiegel uit gereflecteerd licht tevoorschijn tovert. Je hebt er een gezien als de scherpe boog in een koffiekopje, en je kunt er een maken met elke gladde reflector die niet vlak is. Het overzicht van caustieken introduceert het idee naast breking, en deze pagina volgt de reflectiehelft van dat verhaal volledig. We beginnen bij de manier waarop één straal stuitert, bouwen het wiskundige object dat een hele familie van stuiters beschrijft, passen het toe op de twee spiegels uit de interactieve explorers en eindigen met het klassieke resultaat dat de ideale koffiekopkromme een nefroïde is.

De enige voorkennis die je nodig hebt, is de reflectieregel van gelijke hoeken en een beetje vectorrekenen. Als je vectoren kunt optellen en een inwendig product kunt uitrekenen, kun je elke stap hieronder volgen.

De Reflectiewet als Richtingsupdate

Begin met één straal. Die reist in een bepaalde richting naar een oppervlak, raakt het oppervlak in een punt en vertrekt in een nieuwe richting. De regel van gelijke hoeken zegt dat de invallende en uitgaande straal dezelfde hoek maken met de oppervlaktenormaal, de denkbeeldige pijl die loodrecht uit het oppervlak steekt. Reflectie is dus een regel die een richting vervangt: geef hem een invallende richting en een normaal, en hij geeft de uitgaande richting terug.

In vectorvorm is de regel kort. Voor een invallende eenheidsrichting $\mathbf{i}$ en een eenheidsoppervlaktenormaal $\mathbf{n}$ geldt

$$\mathbf{r} = \mathbf{i} - 2(\mathbf{i} \cdot \mathbf{n}),\mathbf{n}$$

waarbij $\mathbf{r}$ de uitgaande eenheidsrichting is. Lees de termen één voor één. Het inwendige product $\mathbf{i} \cdot \mathbf{n}$ meet hoeveel van de invallende richting langs de normaal wijst, en het loopt van $-1$ tot $1$. Dat getal met $2$ en met $\mathbf{n}$ vermenigvuldigen bouwt een vector die het deel van $\mathbf{i}$ dat naar de normaal toe wijst precies opheft en daarna omkeert, want twee keer een component aftrekken is hetzelfde als het teken ervan omdraaien. Het zijwaartse deel van de richting, de component loodrecht op de normaal, blijft onaangeroerd.

Daarom reproduceert de regel het gedrag van gelijke hoeken. Splits de invallende richting in een stuk langs de normaal en een stuk er dwars op. Reflectie klapt het stuk langs de normaal om en behoudt het stuk er dwars op, en één component omklappen is hetzelfde als de hele pijl spiegelen over een lijn die langs het oppervlak loopt. De formule doet die omklap in één stap, zonder hoeken of goniometrie.

Twee eigenschappen zijn belangrijk voor de rest van het artikel. Ten eerste komt alleen de normaal op het trefpunt de regel binnen, dus elk punt van een oppervlak werkt op zichzelf. Ten tweede doet het teken van de normaal er niet toe, want $\mathbf{n}$ vervangen door $-\mathbf{n}$ laat de formule onveranderd. Die vrijheid laat ons de richting kiezen die voor de geometrie het handigst is.

Een platte spiegel heeft overal dezelfde normaal. Elke straal ondergaat dezelfde omklap, dus gelijkmatig verdeelde stralen blijven gelijkmatig verdeeld en de familie concentreert zich nooit. Kromming is het ingrediënt dat een caustiek mogelijk maakt. Op een gekromd oppervlak verandert de normaal van punt tot punt, en naburige gereflecteerde stralen blijven niet langer evenwijdig. Sommige draaien naar elkaar toe, kruisen elkaar en hopen licht op tot een heldere lijn.

Waar Naburige Gereflecteerde Stralen Samenvouwen

Eén straal alleen zegt niets over helderheid. Het patroon zit in de familie: zend een dichte waaier van stralen uit en kijk waar ze terechtkomen. Op een gekromde spiegel kunnen we de familie indexeren met het punt waar elke straal het oppervlak raakt. Laat $\mathbf{P}(t)$ dat trefpunt zijn als functie van één parameter $t$, en laat $\mathbf{r}(t)$ de gereflecteerde richting daar zijn. Elke straal is dan de lijn

$$\mathbf{L}(t, s) = \mathbf{P}(t) + s,\mathbf{r}(t)$$

voor afstand $s$ gemeten vanaf het oppervlak. Verschillende waarden van $t$ geven verschillende stralen. Als twee naburige stralen, bij $t$ en $t + \Delta t$, elkaar ergens kruisen, dan zijn twee dicht bij elkaar gelegen delen van de lichtbundel door dezelfde plek geperst. Terwijl $\Delta t$ naar nul krimpt, bewegen de kruisingspunten naar één gladde kromme die alle stralen schampen.

Die kromme is de omhullende van de stralenfamilie, en dat is de caustiek. Om hem te vinden vragen we naar het punt op de straal bij $t$ waar de straal raakt aan de omhullende. Bewegen langs de familie verandert zowel het trefpunt als de richting, dus de voorwaarde is dat de verandering van $\mathbf{P}(t)$ naar $\mathbf{P}(t) + \Delta t,\mathbf{P}'(t)$ en de verandering in richting evenwijdig blijven aan de straal zelf. In twee dimensies is dat een determinantvoorwaarde, waarbij $\det(\mathbf{a}, \mathbf{b}) = a_x b_y - a_y b_x$:

$$\det!\left(\mathbf{P}'(t) + s,\mathbf{r}'(t),\ \mathbf{r}(t)\right) = 0$$

De accenten betekenen afgeleiden naar $t$. De determinant splitsen en oplossen naar de afstand langs de straal geeft direct het omhullingspunt:

$$s = -\frac{\det(\mathbf{P}'(t), \mathbf{r}(t))}{\det(\mathbf{r}'(t), \mathbf{r}(t))}, \qquad \mathbf{C}(t) = \mathbf{P}(t) + s(t),\mathbf{r}(t)$$

Het punt $\mathbf{C}(t)$ is de caustiek. Elke straal in de familie raakt deze kromme in precies één punt en loopt aan weerszijden ervan verder, dus de kromme is een vouw en geen muur. Waar de vouw scherp buigt, hopen de stralen zich nog dichter op, en daarom leest een caustiek als een lijn met heldere, zachte randen in plaats van een harde grens.

De explorers lossen de determinantvoorwaarde niet symbolisch op. Ze traceren een eindig aantal stralen en markeren de kruising van elk aangrenzend paar, wat hetzelfde idee is met een eindige $\Delta t$. Elke stip is een rechtelijnsbenadering van één punt op de vouw. Verhoog het aantal stralen en de stippen dringen samen langs dezelfde kromme, zodat de stippellijn naar de omhullende convergeert. Wanneer twee aangrenzende stralen toevallig evenwijdig zijn, wordt hun kruising naar oneindig geduwd en slaat de code dat paar gewoon over.

De Explorer met de Parabolische Spiegel

De eerste explorer gebruikt een spiegel waarvan de vorm gemakkelijk op te schrijven is en waarvan de normalen gemakkelijk te differentiëren zijn, zodat de aandacht bij de straalgeometrie blijft in plaats van bij de algebra van het oppervlak.

Plaats de top van de parabool in de oorsprong en laat de spiegel

$$\mathbf{P}(x) = (x, c x^2)$$

waarbij het enkele getal $c$ de krommingsparameter is. Een grotere $c$ buigt de spiegel scherper. De raakrichting is de afgeleide $\mathbf{P}'(x) = (1, 2cx)$, en een vector loodrecht op de raaklijn is $(2cx, -1)$, want hun inwendige product is nul. Na normaliseren naar eenheidslengte is de normaal

$$\mathbf{n}(x) = \frac{(2cx, -1)}{\sqrt{1 + (2cx)^2}}$$

Met de bron in een punt $S$ is de invallende richting naar het spiegelpunt $\mathbf{i}(x) = \bigl(\mathbf{P}(x) - S\bigr) / |\mathbf{P}(x) - S|$. $\mathbf{i}$ en $\mathbf{n}$ in de reflectieformule stoppen levert $\mathbf{r}(x)$ op voor elke gesamplede $x$.

In de explorer welft de parabool zich als een brug over de bron, met de top bovenaan het canvas en de bron eronder aan de holle kant. De blauwe lijnen zijn de invallende stralen van de oranje stip naar hun trefpunten, en de rode lijnen zijn de stralen na reflectie. Elke paarse stip markeert waar twee naburige rode stralen elkaar kruisen. In deze explorer worden $120$ stralen gesampled op gelijkmatig verdeelde waarden van $x$ over de spiegel, en de kruisingspunten worden bijgesneden tot het gebied voor de spiegel zodat de schatting zichtbaar blijft.

Een uitgewerkte straal maakt de getallen concreet. Neem $c = 0.0024$ en een spiegelpunt op $x = 100$, zodat de hoogte $c x^2 = 24$ is. De helling van de raaklijn is $2cx = 0.48$, en de eenheidsnormaal is $(0.48, -1)/1.1092 = (0.4327, -0.9016)$. Zet de bron op $(0, 240)$, op de as van de spiegel, $240$ eenheden van de top aan de kant van de bron. De invallende richting is dan evenredig met $(100, 24 - 240) = (100, -216)$, wat normaliseert naar $(0.4202, -0.9076)$. Het inwendige product met de normaal is $0.4202 \cdot 0.4327 + (-0.9076)(-0.9016) = 1.0001$, dus de straal komt bijna precies langs de normaal aan. Tweemaal de normaalcomponent aftrekken stuurt hem terug als $\mathbf{r} = (-0.445, 0.896)$, een richting die hem terug over de as draagt. De straal uit de naburige sample helt net iets anders, en de twee ontmoeten elkaar in een punt dat een paarse stip wordt.

De explorer maakt drie verbanden zichtbaar.

De krommingsschuif bepaalt hoe snel de normalen van de ene sample naar de volgende draaien. Bij de standaardwaarde $c = 0.0024$ vormen de kruisingspunten een smalle, puntige boog waarvan de punt tussen de spiegel en de bron ligt, met de twee armen die voorbij de bron uitwaaieren. Verhoog $c$ naar $0.0045$ en de punt klimt van ongeveer een derde van het canvas naar net onder de spiegel, terwijl de boog dramatisch breder wordt; de paarse stippen strekken zich nu uit over honderden pixels in plaats van een paar dozijn. Verlaag $c$ richting $0.0015$ en de spiegel wordt bijna vlak, de gereflecteerde stralen blijven dicht bij evenwijdig, en de kruisingen verhuizen voorbij de rand van het tafereel, waardoor het canvas bijna leeg blijft.

De positie van de bron bepaalt waar de vouw ligt. De bron omlaag slepen, verder van de spiegel, tilt de punt naar de spiegel toe en verbreedt de boog, dezelfde trend als het verhogen van de kromming, maar van de andere kant bereikt. Hem omhoog slepen, dichter bij de spiegel, laat de punt zakken en vernauwt de vouw tot de kruisingen van de onderkant van het canvas af drijven.

De bron zijwaarts slepen breekt de symmetrie. Verplaats de stip naar rechts en de hele vouw verschuift naar links, weg van de bron leunend en uitrekkend tot een langere, meer open kromme. Met een niet-gecentreerde bron komen de twee armen niet meer overeen, en de puntige tip ligt niet langer op de as van de spiegel. De spiegel vouwt de stralen nog steeds, maar de vouw is scheef.

De Koffiekopkromme en de Nefroïde

De tweede explorer vervangt de parabool door een cirkel, de vorm van de wand binnen een cilindrisch kopje.

Hier is de bron een punt binnen het kopje, en stralen verlaten het in elke richting in plaats van aan te komen vanaf een vaste verzameling oppervlaktepunten. De normaal van de cirkelvormige wand in een trefpunt wijst recht langs de straal, dus de reflectiegeometrie is ongewoon schoon. De explorer sampled $140$ richtingen rond de bron, vindt het eerste snijpunt van elke straal met de cirkel, reflecteert die en traceert ook een tweede stuiter, hoewel alleen de reflecties van de eerste stuiter de paarse kruisingsstippen voeden. Een kleine duw voorbij elk trefpunt voorkomt dat hetzelfde snijpunt op het volgende segment opnieuw wordt gedetecteerd.

Dicht bij het midden van het kopje vouwt de gereflecteerde familie zich tot een nette, ruwweg symmetrische kromme. Sleep de bron naar de wand en de kromme drijft naar de tegenovergestelde kant en scheeft, precies zoals bij de niet-gecentreerde parabool. Sleep hem omhoog of omlaag en de veranderende invalshoek hervormt de vouw. De heldere kromme blijft bestaan voor elke bronpositie binnen het kopje, omdat de normalen van de cirkel de gereflecteerde richtingen blijven sorteren. De vouw is niet het kopje dat een beeld van het licht toont; het is het kopje dat de stralenfamilie herschikt.

De Ideale Kromme Is een Nefroïde

De schoonste versie van het kopjesprobleem vervangt de puntbron door een bundel evenwijdige stralen, alsof het licht van de zon of een zeer verre lamp komt. Zet de cirkel in de oorsprong met straal $R$ en laat de invallende richting $\mathbf{i} = (1, 0)$ zijn. Een straal die reflecteert in het punt

$$\mathbf{P}(\phi) = R,(\cos\phi,\ \sin\phi)$$

ontmoet een radiale eenheidsnormaal $\mathbf{n}(\phi) = (\cos\phi,\ \sin\phi)$. Het inwendige product $\mathbf{i} \cdot \mathbf{n}$ is $\cos\phi$, dus de reflectieformule geeft

$$\mathbf{r}(\phi) = (1,0) - 2\cos\phi,(\cos\phi,\ \sin\phi) = (-\cos 2\phi,\ -\sin 2\phi)$$

De familie van gereflecteerde stralen is $\mathbf{L}(\phi, s) = \mathbf{P}(\phi) + s,\mathbf{r}(\phi)$. Pas nu de omhullingsvoorwaarde van eerder toe. De afgeleiden zijn $\mathbf{P}' = R(-\sin\phi, \cos\phi)$ en $\mathbf{r}' = (2\sin 2\phi, -2\cos 2\phi)$. De twee determinanten hebben eenvoudige waarden:

$$\det(\mathbf{P}', \mathbf{r}) = R\cos\phi, \qquad \det(\mathbf{r}', \mathbf{r}) = -2$$

Het oplossen van $s = -\det(\mathbf{P}', \mathbf{r}) / \det(\mathbf{r}', \mathbf{r})$ geeft $s = \tfrac{R}{2}\cos\phi$, dus het omhullingspunt is

$$\mathbf{C}(\phi) = \mathbf{P}(\phi) + \tfrac{R}{2}\cos\phi;\mathbf{r}(\phi)$$

De twee coördinaten uitwerken:

$$C_x = R\cos\phi - \tfrac{R}{2}\cos\phi\cos 2\phi, \qquad C_y = R\sin\phi - \tfrac{R}{2}\cos\phi\sin 2\phi = R\sin^3\phi$$

De $y$-coördinaat vereenvoudigt omdat $\sin 2\phi = 2\sin\phi\cos\phi$, zodat de twee termen samengaan tot $R\sin\phi,(1 - \cos^2\phi)$. De $x$-coördinaat vereenvoudigt met de dubbele-hoekidentiteit $\cos 2\phi = 2\cos^2\phi - 1$. Beide klappen samen tot een bekende vorm:

$$\mathbf{C}(\phi) = \frac{R}{4}\bigl(3\cos\phi - \cos 3\phi,\ 3\sin\phi - \sin 3\phi\bigr)$$

Dit is de nefroïde, de niervormige kromme die een evenwijdige bundel binnen een cirkel reflecteert. De twee keerpunten liggen op $(\pm R/2, 0)$, binnen de cirkel op de halve straal, en de kromme schampt de cirkel zelf in de twee punten $(0, \pm R)$. Een echte bundel verlicht slechts een deel van de wand, dus de heldere boog die je in een kopje ziet is het deel van de nefroïde dat door dat verlichte stuk wordt gevormd, niet de hele gesloten kromme. Echte kopjes zijn ook geen perfecte cirkels en echt licht is niet perfect evenwijdig, en daarom is de waargenomen kromme een licht vervormde nefroïde. Verplaats de bron van oneindig naar de rand van de cirkel en dezelfde constructie levert een cardioïde op, de naaste verwant van de nefroïde.

Een Uitgewerkt Punt op de Nefroïde

Kies een reflectiepunt op $\phi = 60^\circ$ op een cirkel met straal $R = 1$ en volg de vier stappen.

1. Het trefpunt is $\mathbf{P} = (\cos 60^\circ, \sin 60^\circ) = (0.5, 0.8660)$.
2. De invallende richting is $(1, 0)$, en de radiale normaal is $\mathbf{n} = (0.5, 0.8660)$, dus $\mathbf{i} \cdot \mathbf{n} = 0.5$.
3. De gereflecteerde richting is $\mathbf{r} = (1,0) - 2(0.5)(0.5, 0.8660) = (0.5, -0.8660)$. De omhullingsafstand is $s = \tfrac{R}{2}\cos 60^\circ = 0.25$.
4. Het caustiekpunt is $\mathbf{C} = (0.5, 0.8660) + 0.25,(0.5, -0.8660) = (0.625, 0.6495)$.

Ter controle tegen de nefroïdeformule: $\cos 180^\circ = -1$ en $\sin 180^\circ = 0$, dus $C_x = \tfrac{1}{4}(3 \cdot 0.5 - (-1)) = 0.625$ en $C_y = \tfrac{1}{4}(3 \cdot 0.8660 - 0) = 0.6495$. Beide routes komen overeen. De gereflecteerde straal bij $\phi = 60^\circ$ schampt de caustiek in $(0.625, 0.6495)$, precies zoals de omhullingsvoorwaarde beloofde.

Randgevallen en Limieten

Verschillende situaties laten de vouw zich gedragen op manieren die het plaatje op het eerste gezicht niet laat zien.

Wanneer de bron binnen het brandpunt van een spiegel komt, divergeren de gereflecteerde stralen in plaats van te convergeren en wordt de caustiek virtueel: de stralen lijken alleen maar uit een vouw achter het oppervlak te komen, dus er is geen voorwaarts kruisingspunt om te markeren. Voor een parabool $y = c x^2$ ligt het brandpunt op een afstand $1/(4c)$ van de top. Op de zachtste stand, $c = 0.0015$, ligt dat brandpunt ongeveer $167$ pixels van de top, en een bron die dichterbij staat levert helemaal geen voorwaartse kruisingen op.

Een leeg canvas kan ook iets minder dramatisch betekenen. Beide explorers snijden kruisingspunten bij tot het zichtbare gebied, en een echte vouw kan simpelweg buiten het kader liggen. Met een bijna vlakke spiegel en de bron ver genoeg weg liggen de kruisingspunten honderden pixels onder de onderrand. De kromming verhogen of de bron verder van de spiegel verplaatsen is meestal genoeg om de vouw weer in beeld te trekken.

Evenwijdige aangrenzende stralen zijn het volgende randgeval. Wanneer twee naburige gereflecteerde richtingen precies samenvallen, ligt hun kruising op oneindig en levert het paar geen stip op. Dit gebeurt aan de buitenste uiteinden van het gesamplede bereik van de parabool, waar de normalen het snelst draaien en de gereflecteerde stralen door een richting roteren die een andere straal ook aanneemt. Het kruisende paar valt gewoon uit de schatting weg, waardoor een kleine opening in de stippellijn ontstaat.

Niet-gecentreerde bronnen halen de symmetrie weg die de vouw gemakkelijk leesbaar maakt. Met de bron op de as heeft de caustiek één puntige tip en twee bijpassende armen. Verplaats de bron zijwaarts en de tip draait weg van de as, de ene arm wordt langer dan de andere, en de vouw kan ver naar één kant zwaaien. Hetzelfde gebeurt in een echt kopje wanneer je het licht opzij houdt: de heldere kromme schuift naar de tegenovergestelde wand en rekt uit.

Schampende inval is de laatste limiet om in gedachten te houden. Nabij de rand van het kopje of de uiteinden van de parabool kan een straal het oppervlak bijna rakend treffen. De gereflecteerde richting is dan extreem gevoelig voor de normaal, dus een piepkleine fout in de vorm van het oppervlak zwaait de uitgaande straal over een grote hoek. Dit is een fysiek effect, niet alleen een numeriek effect, en het is de reden dat de buitenste delen van een caustiek waziger kunnen lijken dan het midden.

Samplen en Normalen in de Praktijk

Als je deze plaatjes wilt reproduceren, doen twee keuzes het meeste werk.

De eerste is hoe je de stralenfamilie sampled. De parabolische explorer sampled gelijkmatig verdeelde punten langs de spiegel, de juiste keuze wanneer de normaal met de oppervlaktepositie verandert, omdat elk deel van het oppervlak evenveel gewicht krijgt. De kopjesexplorer sampled juist gelijkmatig verdeelde richtingen rond de bron, de juiste keuze voor een puntbron die in alle richtingen straalt. De sampling afstemmen op de bron van variatie houdt de stippen gelijkmatig verdeeld langs de caustiek in plaats van samen te klonteren waar de parameter langzaam verandert.

De tweede is het aantal stralen. Met te weinig stralen zijn de paarse stippen spaarzaam en leest de vouw als een gebroken lijn. Met meer stralen dringen de stippen samen op de omhullende en wordt de kromme gladder, tegen de prijs van meer berekening en zwaardere overdraw. De explorers gebruiken $120$ stralen voor de spiegel en $140$ voor het kopje, genoeg om de vorm leesbaar te maken zonder hem in stippen te begraven. Het aantal verdubbelen verdikt vooral de bestaande kromme in plaats van nieuwe structuur te onthullen.

De nauwkeurigheid van de normaal is bij schampende hoeken belangrijker dan het aantal stralen. Beide explorers berekenen de normaal analytisch, de ene uit de afgeleide van de parabool en de andere uit de straal van de cirkel, in plaats van hem te schatten uit naburige samples. Een normaal via eindige verschillen versterkt kleine fouten stilletjes precies daar waar de uitgaande richting het gevoeligst is, en dat is waar de caustiek het scherpst is. De exacte normaal gebruiken houdt die fout buiten het plaatje.

Houd ten slotte in gedachten wat het geometrische model weglaat. Het tellen van straalkruisingen laat zien waar licht zich concentreert, niet hoe helder het wordt. Een fysisch nauwkeurig beeld zou elke straal wegen met een Fresnelterm die afhangt van de hoek, van de afgelegde afstand en van hoe het oppervlak licht absorbeert. Die factoren veranderen de balans van helderheid over de caustiek, maar ze verplaatsen de vouw zelf niet, die alleen door geometrie wordt bepaald.

Samenvatting

Een reflectiecaustiek begint met één lokale regel. De reflectiewet van gelijke hoeken zet een invallende richting om in een uitgaande met alleen de oppervlaktenormaal op het trefpunt, compact geschreven als $\mathbf{r} = \mathbf{i} - 2(\mathbf{i} \cdot \mathbf{n})\mathbf{n}$. Een gekromd oppervlak geeft elk punt een net iets andere normaal, dus naburige gereflecteerde stralen draaien naar elkaar toe. Waar ze elkaar kruisen, hoopt licht zich op, en in de limiet van oneindig veel stralen tekenen die kruisingen de omhullende, de caustiek.

De twee interactieve spiegels tonen hetzelfde mechanisme op verschillende schalen. De parabolische spiegel maakt de normaal en de gereflecteerde straal gemakkelijk te berekenen, en hij laat zien hoe kromming, bronhoogte en bronverschuiving de vouw elk hervormen. Het cirkelvormige kopje voegt het schoonste exacte resultaat toe: voor evenwijdige invallende stralen is de caustiek de nefroïde, met keerpunten op de halve straal van het kopje, in een paar regels afgeleid uit de omhullingsvoorwaarde. Een puntbron binnen het kopje levert een vervormde versie van dezelfde vouw op, en dat is wat je in een echt koffiekopje ziet.

Het overzicht van caustieken behandelt de brekingskant van hetzelfde verhaal, waar de wet van Snell de reflectieregel vervangt en een golvend wateroppervlak zonlicht omzet in de verschuivende netwerken op een zwembadvloer.