Learning to absorb light: light organizes, matter remembers

FR | EN

Research theme 07

Learning to absorb light: light organizes, matter remembers

Matter does not merely receive light: it responds to it. Under light’s action, matter transforms, organizes itself, and gradually reshapes its own way of absorbing energy. From mineral formations to plants, this capacity to let an order emerge that fits the flows running through it seems to belong to matter well before life itself appeared. What we call learning might then be understood differently: as a system’s capacity to inscribe, in its own structure, the effects of what it has experienced.

Complexity, self-organization and machine learning for ultrafast laser processing. With Fayad Ali Banna, Eduardo Brandao, Anthony Nakhoul, Rémi Emonet and Marc Sebban, in collaboration with the Inria MALICE team.

Cycle linking far-from-equilibrium dynamics, dissipative structures, frozen pattern and increased absorptivity
The loop that makes it a learning process. Each double pulse drives the surface far from equilibrium, dissipative instabilities select the morphologies that absorb best, resolidification freezes them, and the stored complexity raises the absorptivity seen by the next pulse. Nothing here is optimization in the algorithmic sense: the feedback is entirely physical, and the memory is the relief itself.

01A surface that keeps a record of the light it has seen

SEM images of nickel surfaces irradiated with cross-polarized double pulses
Fifty pulses, four routes. Nickel irradiated with cross-polarized double pulses, imaged as the number of double pulses grows from 0 to 50, for four combinations of fluence and interpulse delay. Each series runs through the same sequence, diffuse roughness, oriented ripples, a two-dimensional network, then disorder, but reaches it at its own pace.

The experiment is deliberately simple. Nickel is irradiated with cross-polarized double pulses, two orthogonal linear polarizations fired in succession with a controlled delay, and the surface is imaged after an increasing number of pulses, up to fifty, for several combinations of fluence and interpulse delay. What appears is a sequence: first a diffuse roughness, then oriented ripples, then a well-formed two-dimensional network, then a breakdown into disorder.

The important point is that the surface at pulse N is the initial condition for pulse N+1, and not in a trivial way. The frozen relief determines where the next pulse deposits its energy, which determines where the melt flows, which determines the next relief. This is the same electromagnetic feedback described in the theme on nanoscale electromagnetism, but followed here over tens of pulses and read as an information process rather than a pattern-formation one.

02Two numbers: how much is absorbed, and how organized it is

The first quantity is the absorbed energy, computed from the measured height map by solving the electromagnetic problem on the real surface and integrating the deposited energy. It is averaged over six polarization angles so that it measures the propensity of the morphology to absorb, not the answer to one particular orientation.

The second quantity has to capture organization without rewarding either perfect order or pure noise. A perfect grating and white noise are both uninteresting; what matters is the regime in between.

equation

Energy absorbed by the surface at pulse N, obtained by integrating over the surface the depth-integrated absorption computed from the measured topography. Normalized by the value for a flat surface, so that different fluences can be compared.

equation

Statistical complexity in the López-Ruiz, Mancini and Calbet sense: Shannon entropy times disequilibrium. Here in its Taylor form, adapted to surface height maps.

Statistical complexity does exactly that. It is the product of the Shannon entropy, which grows with randomness, by the disequilibrium, which measures the distance from a uniform distribution and therefore grows with structure. The product is maximal for a system that is neither too ordered nor too random, which is precisely the signature of a dissipative structure.

Plotting the two quantities against the number of pulses for four independent experimental series gives the central result of this work: they rise together, plateau together, and then part company.

03Four phases, and the criterion that separates them

Absorbed energy and complexity versus number of double pulses, with mutual information
Two observables, four phases. Absorbed energy (solid) and complexity (dashed) against the number of double pulses for the four series, with the mutual information below. Phases I to III form the learning regime, response, iterative learning and memory stabilization, while phase IV is the chaotic regime where the relief breaks down and the two quantities lose their correlation.

The trajectory splits into four phases, and the boundaries are not drawn by eye. They follow from the mutual information between absorbed energy and complexity, computed cumulatively as the number of pulses grows, which measures how strongly the two are coupled including through nonlinear dependencies.

In the response phase the surface merely detects the stimulus: roughness appears, hydrodynamic instabilities are transient, and both quantities stay flat. The iterative learning phase begins when two conditions are met at once, a roughness amplitude large enough to create significant near-field inhomogeneities, and thermal conditions allowing Rayleigh-Bénard-Marangoni convection to develop within the lifetime of the molten layer. Absorption and complexity then climb together.

The memory stabilization phase is a kinetic arrest: the time needed for the convective pattern to reorganize exceeds the melt lifetime, so the morphology freezes in an optimized configuration and complexity reaches a plateau. Finally the destruction phase sets in when the deposited energy exceeds what structure-preserving processes can absorb. Enhanced absorption at the optimized nanostructures produces thermal runaway at the hotspots, material is removed, and the two quantities decouple: absorption keeps rising while complexity falls. The drop in mutual information marks that transition, and it is the boundary between a learning regime and a chaotic one.

The system overlearns. Prolonged exposure produces a state so finely optimized for absorption that it destroys itself.

04Why the trajectory has a direction

A history-dependent process is not automatically a learning one. What makes this trajectory directional rather than a random walk in configuration space is a physical feedback loop with a definite sign.

Each pulse deposits energy with a spatial distribution set by the existing topography. Where absorption is locally higher, the melt is hotter, the thermal gradients are stronger and the material transport is faster, so those modulations grow at a higher rate. When resolidification freezes the relief, the geometry that enhanced absorption is the one that survives, and it enhances absorption again at the next pulse. Morphologies that couple light in are therefore selected, not by an optimizer but by the dissipation itself.

The loop runs: far-from-equilibrium dynamics, dissipative instability, freezing by resolidification, increase in complexity, increase in absorptivity, and back to far-from-equilibrium dynamics. It has a close analogue outside the laboratory in snow penitentes, where pointed structures emerge under solar irradiation because differential absorption and sublimation couple geometry to energy dissipation in exactly the same way.

05The surface answers best to the polarization that made it

Effect of the laser polarization angle on the absorbed energy of the optimal patterns
The surface remembers an orientation. For each series, the height map of the pattern that maximizes the absorption of the next pulse, then the energy deposited by that next pulse for polarization angles from 0 to 150 degrees, and finally the integrated absorbed energy as a function of pulse number and angle. The optimum is not isotropic: it favours the polarization that built it.
Polarization direction rotating over a laser-textured surface, with the absorbed energy turning with it
Sweeping the polarization over a surface that already carries a pattern. The absorbed energy turns with the field, and the responses at 0 and 90 degrees bound the whole family of intermediate angles.

If the surface really stores information about the light it has seen, then it should respond preferentially to that light. Testing this means taking a pattern that maximizes absorption and recomputing the absorption for six polarization angles.

Structures parallel to the field absorb more, and the peaks of highest curvature dominate while the valleys stay dark. As the number of pulses grows, curves that were superimposed separate, and the two experimental polarization angles, the ones actually used to write the structure, become the upper and lower bounds of the whole family. The surface has become anisotropic in its sensitivity, and the anisotropy is aligned with its own history.

This is the behaviour that non-associative learning calls sensitization: repeated exposure to a stimulus increases the response to that stimulus. In an organism it comes from neuroplasticity. Here it comes from a relief whose orientation was written by the field it now responds to.

06How fragile is the optimum?

Relative change in absorbed energy when parts of the optimal pattern are replaced
Damaging the memory on purpose. Relative change in absorbed energy against the Frobenius distance between the experimental pattern and an artificial version, obtained either by transplanting a patch from another surface of the same series or by inserting a patch generated with a Swift-Hohenberg solver. Even large geometrical edits cost only a few percent of absorption, so what is stored is a statistical organization rather than one particular arrangement of ridges.

An optimum reached by a physical process should degrade when perturbed, and the way it degrades says what kind of optimum it is. Two perturbations were applied to the pattern that maximizes absorption in each series. The first replaces a patch by another patch taken from a different pattern of the same series: physically plausible material, but locally violating the self-organization law that produced the rest.

The result is unambiguous. Whatever the amplitude of the change, the absorbed energy falls below the reference. The configuration is finely tuned, and a local deviation from the law that generated it costs performance.

equation

Relative change in absorbed energy when a pattern is perturbed, plotted against the Frobenius norm of the difference between the original and the modified height maps.

equation

Swift-Hohenberg equation, the amplitude equation used to generate synthetic patches. It reproduces the pattern-selection dynamics of the thermoconvective instability.

The second perturbation replaces the patch by one generated with a Swift-Hohenberg solver, an amplitude equation that models the hydrodynamic fluctuations and reproduces the same class of patterns. Here the absorbed energy scatters around the reference rather than below it, and stays close to it whatever the amplitude of the change.

The two tests together give the sharpest statement of the paper. The surface is not merely rough in a way that happens to absorb well: it is organized according to a specific dynamical law, and only modifications that respect that law preserve its performance.

07Learning the absorption operator

There is a second, more literal sense in which learning enters this problem. Everything above depends on knowing where the light is absorbed on a given surface, and that map has to be recomputed every time the surface changes. Solving Maxwell’s equations by FDTD on a 320×320 height map takes minutes; following a pulse train, over many realizations, needs tens of thousands of such solves. The absorption step, not the hydrodynamics, is the bottleneck.

Rather than accelerating the solver, we learn the operator it implements (the map that takes a topography to the absorbed energy density it produces), by minimizing the discrepancy with the solver over a dataset of surfaces.

Equation
Building the dataset. The solver produces pairs of topography and absorbed energy; each pair is duplicated at 0° and 90° so that the network sees...
Building the dataset. The solver produces pairs of topography and absorbed energy; each pair is duplicated at 0° and 90° so that the network sees the polarization direction as an explicit symmetry of the data rather than having to infer it.
Topography, FDTD absorbed energy, learned operator prediction and relative error for four surfaces
Four surfaces, one operator. Left to right: the topography, the absorbed energy computed by FDTD, the prediction of the learned operator, and the relative error, essentially featureless on the scale of the field itself. Each row is a different morphology drawn from the test set.

Row after row, the learned operator reproduces the solver to within a few percent on every morphology in the dataset, and evaluates a map more than a hundred thousand times faster. The error panel is the one that matters: it stays flat and structureless, so the surrogate is not simply smoothing the fine detail away. That alone is what makes a stochastic, pulse-by-pulse study of self-organization possible at all.

The interesting part, however, is which architecture wins. Four were put in competition, chosen because they encode two opposite hypotheses about the range of the light–matter coupling. A U-Net is a purely local operator: stacked 3×3 convolutions with a finite receptive field, a few tens of nanometres on the surface. A Fourier Neural Operator is a global one: it multiplies in reciprocal space and can in principle carry a coupling between any two points of the map. Between them sit an FNO with added convolutions, and a U-shaped Neural Operator combining both ideas.

The plain local network wins, in accuracy and in cost. That is not a disappointment for the elegant architectures: it is a physical result. At scales far below the wavelength, the absorbed energy at a point is fixed by the topography in its immediate neighbourhood. The light here is a near field, and the operator that describes it is short-ranged.

The winning architecture: an encoder–decoder of 3×3 convolutions with skip connections, mapping a 320×320 topography directly onto the absorbed...
The winning architecture: an encoder–decoder of 3×3 convolutions with skip connections, mapping a 320×320 topography directly onto the absorbed energy map. Its receptive field is finite by construction.
The architecture built for nonlocality: a U-shaped Neural Operator whose Fourier layers apply a global integral kernel in reciprocal space. Here...
The architecture built for nonlocality: a U-shaped Neural Operator whose Fourier layers apply a global integral kernel in reciprocal space. Here it has nothing to carry.

By comparing architectures with different receptive fields, the network therefore measures the range of the electromagnetic coupling on a rough surface, and confirms from data what the near-field analysis argues from the fields themselves. This is the clearest example on this site of a network used not as a predictor but as an instrument.

08What learning means here, and what it does not

Learning is used in the broad sense: a history-dependent evolution that improves performance with experience. Nothing in this system computes, and no analogy with a specific machine-learning architecture is intended. The optimization strategy is inherent to the physics, to Rayleigh-Bénard-Marangoni instabilities coupled to electromagnetic feedback.

What the work does provide is a way of measuring that process rather than describing it qualitatively, and a way of predicting it. Learning the partial differential equation behind an observed sequence of patterns, and using complexity as the quantity to be guided, turns self-organization into something one can steer towards a target morphology instead of merely observing. That is the practical end of this theme, and it connects directly to the electromagnetic and atomistic models developed in the other three.

Key publications

  1. Photonic learning in ultrafast laser-induced complexityF. Ali Banna, E. Brandao, A. Nakhoul, R. Emonet, M. Sebban, J.-P. Colombier, Ultrafast Science, under revision (2026).
  2. Learning complexity to guide light-induced self-organized nanopatternsE. Brandao et al., Physical Review Letters 130, 226201 (2023). doi.org/10.1103/PhysRevLett.130.226201
  3. Learning PDE to model self-organization of matterE. Brandao et al., Entropy 24, 1096 (2022). doi.org/10.3390/e24081096
  4. Beyond the microscale: advances in surface nanopatterning by laser-driven self-organizationA. Nakhoul, J.-P. Colombier, Laser & Photonics Reviews (2024). doi.org/10.1002/lpor.202301019
← Back to research themes | All publications →