Quantum and classical physics can be used for mathematical computations that are hard to tackle by conventional electronics. Very recently, optical Ising machines have been demonstrated for computing the minima of spin Hamiltonians, paving the way to new ultra-fast hardware for machine learning. However, the proposed systems are either tricky to scale or involve a limited number of spins. We design and experimentally demonstrate a large-scale optical Ising machine based on a simple setup with a spatial light modulator. By encoding the spin variables in a binary phase modulation of the field, we show that light propagation can be tailored to minimize an Ising Hamiltonian with spin couplings set by input amplitude modulation and a feedback scheme. We realize configurations with thousands of spins that settle in the ground state in a low-temperature ferromagnetic-like phase with all-to-all and tunable pairwise interactions. Our results open the route to classical and quantum photonic Ising machines that exploit light spatial degrees of freedom for parallel processing of a vast number of spins with programmable couplings.A large number of internal states characterizes complex systems from biology to social science. The fact that the number of these states grows exponentially with the system size hampers large-scale computational possibilities. Complex optimization problems involving these models are in many cases classified as NP-hard and cannot be tackled efficiently by standard computing architectures. A broad class of such computationally intractable problems maps to the search of the ground state of a classical system of interacting spins: the minimization of an Ising Hamiltonian with specific spin couplings [1-3].Growing research interest is emerging towards physical and artificial systems that evolve according to an Ising Hamiltonian and enable to find the optimal combinatorial solution by the ground state observed in the experiment. Quantum and classical Ising systems have been realized by trapped atoms [4,5], single photons [6], superconducting circuits [7], electromechanical modes [8], nanomagnets [9] and polariton condensates [10]. In optics, spin-glass dynamics have been observed in random lasers [11, 12], multimodal cavities [13, 14] , coupled laser lattices [15], beam filamentation [16] and nonlinear wave propagation in disordered media [17]. These photonic systems host thousands of optical spins, but the spin variables are not easy to access and controlling their interaction is challenging.Novel photonic platforms with numerous and easily accessible spins are particularly relevant for computation. Optical computing machines offer high-speed and parallelization. Various authors reported coherent Ising machines based on time-multiplexed optical parametric oscillators finding approximate solutions to optimization problems with several nodes [18][19][20][21][22][23][24]. Others proposed nanophotonic circuits to implement any small-scale spin systems directly on a programmable chip [25][26][27]. Matrix opera...

We study artificial neural networks with nonlinear waves as a computing reservoir. We discuss universality and the conditions to learn a dataset in terms of output channels and nonlinearity. A feed-forward three-layer model, with an encoding input layer, a wave layer, and a decoding readout, behaves as a conventional neural network in approximating mathematical functions, realworld datasets, and universal Boolean gates. The rank of the transmission matrix has a fundamental role in assessing the learning abilities of the wave. For a given set of training points, a threshold nonlinearity for universal interpolation exists. When considering the nonlinear Schrödinger equation, the use of highly nonlinear regimes implies that solitons, rogue, and shock waves do have a leading role in training and computing. Our results may enable the realization of novel machine learning devices by using diverse physical systems, as nonlinear optics, hydrodynamics, polaritonics, and Bose-Einstein condensates. The application of these concepts to photonics opens the way to a large class of accelerators and new computational paradigms. In complex wave systems, as multimodal fibers, integrated optical circuits, random, topological devices, and metasurfaces, nonlinear waves can be employed to perform computation and solve complex combinatorial optimization.

One of the most controversial phenomena in nonlinear dynamics is the reappearance of initial conditions. Celebrated as the Fermi-Pasta-Ulam-Tsingou problem, the attempt to understand how these recurrences form during the complex evolution that leads to equilibrium has deeply influenced the entire development of nonlinear science. The enigma is rendered even more intriguing by the fact that integrable models predict recurrence as exact solutions, but the difficulties involved in upholding integrability for a sufficiently long dynamic has not allowed a quantitative experimental validation. In natural processes, coupling with the environment rapidly leads to thermalization, and finding nonlinear multimodal systems presenting multiple returns is a long-standing open challenge. Here, we report the observation of more than three Fermi-Pasta-Ulam-Tsingou recurrences for nonlinear optical spatial waves and demonstrate the control of the recurrent behavior through the phase and amplitude of the initial field. The recurrence period and phase shift are found to be in remarkable agreement with the exact recurrent solution of the nonlinear Schrödinger equation, while the recurrent behavior disappears as integrability is lost. These results identify the origin of the recurrence in the integrability of the underlying dynamics and allow us to achieve one of the basic aspirations of nonlinear dynamics: the reconstruction, after several return cycles, of the exact initial condition of the system, ultimately proving that the complex evolution can be accurately predicted in experimental conditions.

Rogue waves are observed as light propagates in the extreme nonlinear regime that occurs when a photorefractive ferroelectric crystal is undergoing a structural phase transition. The transmitted spatial light distribution contains bright localized spots of anomalously large intensity that follow a signature long-tail statistics that disappears as the nonlinearity is weakened. The isolated wave events form as out-of-equilibrium response and disorder enhance the Kerr-saturated nonlinearity at the critical point. Self-similarity associable to the individual observed filaments and numerical simulations of the generalized nonlinear Schrödinger equation suggests that dynamics of soliton fusions and scale invariance can microscopically play an important role in the observed rogue intensities and statistics.

Ising machines are novel computing devices for the energy minimization of Ising models. These combinatorial optimization problems are of paramount importance for science and technology, but remain difficult to tackle on large scale by conventional electronics. Recently, various photonics-based Ising machines demonstrated fast computing of a Ising ground state by data processing through multiple temporal or spatial optical channels. Experimental noise acts as a detrimental effect in many of these devices. On the contrary, here we demonstrate that an optimal noise level enhances the performance of spatial-photonic Ising machines on frustrated spin problems. By controlling the error rate at the detection, we introduce a noisy-feedback mechanism in an Ising machine based on spatial light modulation. We investigate the device performance on systems with hundreds of individually-addressable spins with all-to-all couplings and we found an increased success probability at a specific noise level. The optimal noise amplitude depends on graph properties and size, thus indicating an additional tunable parameter helpful in exploring complex energy landscapes and in avoiding getting stuck in local minima. Our experimental results identify noise as a potentially valuable resource for optical computing. This concept, which also holds in different nanophotonic neural networks, may be crucial in developing novel hardware with optics-enabled parallel architecture for large-scale optimizations.

From optics to hydrodynamics, shock and rogue waves are widespread. Although they appear as distinct phenomena, new theories state that transitions between extreme waves are allowed. However, these have never been experimentally observed because of the lack of control strategies.We introduce a new concept of nonlinear wave topological control, based on the one-to-one correspondence between the number of wave packet oscillating phases and the genus of toroidal surfaces associated with the nonlinear Schrödinger equation solutions by the Riemann theta function. We prove it experimentally by reporting the first observation of supervised transitions between extreme waves with different genera, like the continuous transition from dispersive shock to rogue waves. Specifically, we use a parametric time-dependent nonlinearity to shape the asymptotic wave genus. We consider the box problem in a focusing Kerr-like photorefractive medium and tailor time-dependent propagation coefficients, as nonlinearity and dispersion, to explore each region in the state-diagram and include all the dynamic phases in the nonlinear wave propagation.Our result is the first example of the topological control of integrable nonlinear waves. This new technique casts light on dispersive shock waves and rogue wave generation, and can be extended to other nonlinear phenomena, from classical to quantum ones. The outcome is not only important for fundamental studies and control of extreme nonlinear waves, but can be also applied to spatial beam shaping for microscopy, medicine and spectroscopy, and to the broadband coherent light generation.

Combinatorial optimization problems are crucial for widespread applications but remain difficult to solve on a large scale with conventional hardware. Novel optical platforms, known as coherent or photonic Ising machines, are attracting considerable attention as accelerators on optimization tasks formulable as Ising models. Annealing is a well-known technique based on adiabatic evolution for finding optimal solutions in classical and quantum systems made by atoms, electrons, or photons. Although various Ising machines employ annealing in some form, adiabatic computing on optical settings has been only partially investigated. Here, we realize the adiabatic evolution of frustrated Ising models with 100 spins programmed by spatial light modulation. We use holographic and optical control to change the spin couplings adiabatically, and exploit experimental noise to explore the energy landscape. Annealing enhances the convergence to the Ising ground state and allows to find the problem solution with probability close to unity. Our results demonstrate a photonic scheme for combinatorial optimization in analogy with adiabatic quantum algorithms and classical annealing methods but enforced by optical vector-matrix multiplications and scalable photonic technology.

No abstract