Search

Ising machines are an emerging class of hardware that promises ultrafast and energy-efficient solutions to NP-hard combinatorial optimization problems. Spatial photonic Ising machines (SPIMs) exploit optical computing in free space to accelerate the computation, showcasing parallelism, scalability, and low power consumption. However, current SPIMs can implement only a restricted class of problems. This partial programmability is a critical limitation that hampers their benchmark. Achieving full programmability of the device while preserving its scalability is an open challenge. Here, we report a fully programmable SPIM achieved through a novel operation method based on the division of the focal plane. In our scheme, a general Ising problem is decomposed into a set of Mattis Hamiltonians, whose energies are simultaneously computed optically by measuring the intensity on different regions of the camera sensor. Exploiting this concept, we experimentally demonstrate the computation with high success probability of ground-state solutions of up to 32-spin Ising models on unweighted maximum cut graphs with and without ferromagnetic bias. Simulations of the hardware prove a favorable scaling of the accuracy with the number of spins. Our fully programmable SPIM enables the implementation of many quadratic unconstrained binary optimization problems, further establishing SPIMs as a leading paradigm in non von Neumann hardware.

We report the observation of a two-dimensional dam break flow of a photon fluid in a nonlinear optical crystal. By precisely shaping the amplitude and phase of the input wave, we investigate the transition from one-dimensional (1D) to two-dimensional (2D) nonlinear dynamics. We observe wave breaking in both transverse spatial dimensions with characteristic timescales determined by the aspect ratio of the input box-shaped field. The interaction of dispersive shock waves propagating in orthogonal directions gives rise to a 2D ensemble of solitons. Depending on the box size, we report the evidence of a dynamic phase characterized by a constant number of solitons, resembling a 1D solitons gas in integrable systems. We measure the statistical features of this gaseous-like phase. Our findings pave the way to the investigation of collective solitonic phenomena in two dimensions, demonstrating that the loss of integrability does not disrupt the dominant phenomenology.

We explore the potential of spatial-photonic Ising machines (SPIMs) to address computationally intensive Ising problems that employ low-rank and circulant coupling matrices. Our results indicate that the performance of SPIMs is critically affected by the rank and precision of the coupling matrices. By developing and assessing advanced decomposition techniques, we expand the range of problems SPIMs can solve, overcoming the limitations of traditional Mattis-type matrices. Our approach accommodates a diverse array of coupling matrices, including those with inherently low ranks, applicable to complex NP-complete problems. We explore the practical benefits of low-rank approximation in optimization tasks, particularly in financial optimization, to demonstrate the real-world applications of SPIMs. Finally, we evaluate the computational limitations imposed by SPIM hardware precision and suggest strategies to optimize the performance of these systems within these constraints., Comment: 15 pages, 7 figures

Non-deterministic polynomial-time (NP) problems are ubiquitous in almost every field of study. Recently, all-optical approaches have been explored for solving classic NP problems based on the spin-glass Ising Hamiltonian. However, obtaining programmable spin-couplings in large-scale optical Ising simulators, on the other hand, remains challenging. Here, we demonstrate control of the interaction length between user-defined parts of a fully-connected Ising system. This is achieved by exploiting the knowledge of the transmission matrix of a random medium and by using diffusers of various thickness. Finally, we exploit our spin-coupling control to observe replica-to-replica fluctuations and its analogy to standard replica symmetry breaking., Comment: arXiv admin note: text overlap with arXiv:2111.07893

States of light encoding multiple polarizations - vector beams - offer unique capabilities in metrology and communication. However, their practical application is limited by the lack of methods for measuring many polarizations in a scalable and compact way. Here we demonstrate polarimetry of vector beams in a single shot without any polarization optics. We map the beam polarization content into a spatial intensity distribution through multiple light scattering and exploit supervised learning for single-shot measurements of multiple polarizations. The method also allows us to classify beams with an unknown number of polarization modes, a functionality missing in conventional techniques. Our findings enable a fast and compact polarimeter for polarization-structured light, a universal tool that may radically impact optical devices for sensing, imaging, and computing.

Modern machine learning applications require huge artificial networks demanding in computational power and memory. Light-based platforms promise ultra-fast and energy-efficient hardware, which may help in realizing next-generation data processing devices. However, current photonic networks are limited by the number of input-output nodes that can be processed in a single shot. This restricted network capacity prevents their application to relevant large-scale problems such as natural language processing. Here, we realize a photonic processor with a capacity exceeding $1.5 \times 10^{10}$ optical nodes, more than one order of magnitude larger than any previous implementation, which enables photonic large-scale text encoding and classification. By exploiting the full three-dimensional structure of the optical field propagating in free space, we overcome the interpolation threshold and reach the over-parametrized region of machine learning, a condition that allows high-performance natural language processing with a minimal fraction of training points. Our results provide a novel solution to scale-up light-driven computing and open the route to photonic language processing., Comment: 12 pages, 4 figures, 1 table

The race to heuristically solve non-deterministic polynomial-time (NP) problems through efficient methods is ongoing. Recently, optics was demonstrated as a promising tool to find the ground state of a spin-glass Ising Hamiltonian, which represents an archetypal NP problem. However, achieving completely programmable spin couplings in these large-scale optical Ising simulators remains an open challenge. Here, by exploiting the knowledge of the transmission matrix of a random medium, we experimentally demonstrate the possibility of controlling the couplings of a fully connected Ising spin system. By further tailoring the input wavefront we showcase the possibility of modifying the Ising Hamiltonian both by accounting for an external magnetic field and by controlling the number of degenerate ground states and their properties and probabilities. Our results represent a relevant step toward the realisation of fully-programmable Ising machines on thin optical platforms, capable of solving complex spin-glass Hamiltonians on a large scale., Comment: 13 pages, 9 figures

Networks of optical oscillators simulating coupled Ising spins have been recently proposed as a heuristic platform to solve hard optimization problems. These networks, called coherent Ising machines (CIMs), exploit the fact that the collective nonlinear dynamics of coupled oscillators can drive the system close to the global minimum of the classical Ising Hamiltonian, encoded in the coupling matrix of the network. To date, realizations of large-scale CIMs have been demonstrated using hybrid optical-electronic setups, where optical oscillators simulating different spins are subject to electronic feedback mechanisms emulating their mutual interaction. While the optical evolution ensures an ultrafast computation, the electronic coupling represents a bottleneck that causes the computational time to severely depend on the system size. Here, we propose an all-optical scalable CIM with fully-programmable coupling. Our setup consists of an optical parametric amplifier with a spatial light modulator (SLM) within the parametric cavity. The spin variables are encoded in the binary phases of the optical wavefront of the signal beam at different spatial points, defined by the pixels of the SLM. We first discuss how different coupling topologies can be achieved by different configurations of the SLM, and then benchmark our setup with a numerical simulation that mimics the dynamics of the proposed machine. In our proposal, both the spin dynamics and the coupling are fully performed in parallel, paving the way towards the realization of size-independent ultrafast optical hardware for large-scale computation purposes., Comment: 7 pages, 3 figures

Photonic brain-inspired platforms are emerging as novel analog computing devices, enabling fast and energy-efficient operations for machine learning. These artificial neural networks generally require tailored optical elements, such as integrated photonic circuits, engineered diffractive layers, nanophotonic materials, or time-delay schemes, which are challenging to train or stabilize. Here we present a neuromorphic photonic scheme - photonic extreme learning machines - that can be implemented simply by using an optical encoder and coherent wave propagation in free space. We realize the concept through spatial light modulation of a laser beam, with the far field that acts as feature mapping space. We experimentally demonstrated learning from data on various classification and regression tasks, achieving accuracies comparable to digital extreme learning machines. Our findings point out an optical machine learning device that is easy-to-train, energetically efficient, scalable and fabrication-constraint free. The scheme can be generalized to a plethora of photonic systems, opening the route to real-time neuromorphic processing of optical data., Comment: 12 pages, 4 figures

We perform real-time stereoscopic wide-area imaging of a ferroelectric phase-transition in KTN:Li. Spontaneous polarization is observed to form a thermally hysteretic 3D lattice of mutually interlinked closed-flux patterns that spans the entire sample. Results are compatible with a supercrystal of $N=4$ topological texture defects arising as the three-fold spatial and one-fold time-inversion symmetries are simultaneously broken. Each lattice site of the texture supercrystal emerges as the projection in actual space of an S$^3$ hypersphere, an extended volume Hopf-link fabric able to screen both volume charge and ferroelectric strain.

Many developments in science and engineering depend on tackling complex optimizations on large scales. The challenge motivates intense search for specific computing hardware that takes advantage from quantum features, nonlinear dynamics, or photonics. A paradigmatic optimization problem is finding low-energy states in classical spin systems with fully-random interactions. To date no alternative computing platform can address such spin-glass problems on a large scale. Here we propose and realize an optical scalable spin-glass simulator based on spatial light modulation and multiple light scattering. By tailoring optical transmission through a disordered medium, we optically accelerate the computation of the ground state of large spin networks with all-to-all random couplings. Scaling of the operation time with the problem size demonstrates optical advantage over conventional computing. Our results point out optical vector-matrix multiplication as a tool for spin-glass problems and provide a general route towards large-scale computing that exploits speed, parallelism and coherence of light., Comment: 9 pages, 4 figures

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 ultra-fast computing of 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, we here 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 trapping into local minima. The result points out noise as a 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., Comment: 7 pages, 4 figures

Nonlinear response in a material increases with its index of refraction as $n^4$. Commonly, $n \sim$ 1 so that diffraction, dispersion, and chromatic walk-off limit nonlinear scattering. Ferroelectric crystals with a periodic 3D polarization structure overcome some of these constraints through versatile Cherenkov and quasi-phase-matching mechanisms. Three-dimensional self-structuring can also lead to a giant optical refraction. Here, we perform second-harmonic-generation experiments in KTN:Li in conditions of giant broadband refraction. Enhanced response causes wavelength conversion to occur in the form of bulk Cherenkov radiation without diffraction and chromatic walk-off, even in the presence of strong wave-vector mismatch and highly focused beams. The process occurs with a wide spectral acceptance of more than 100 nm in the near infrared spectrum, an ultra-wide angular acceptance of up to $\pm 40^{\circ}$, with no polarization selectivity, and can be tuned to allow bulk supercontinuum generation. Results pave the way to highly efficient and adaptable nonlinear optical devices with the promise of single-photon-to-single-photon nonlinear optics., Comment: 10 pages, 5 figures

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, real-world 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 Schroedinger 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., Comment: 5 pages, 4 figures

Dispersive shock waves in thermal optical media belong to the third-order nonlinear phenomena, whose intrinsic irreversibility is described by time asymmetric quantum mechanics. Recent studies demonstrated that nonlocal wave breaking evolves in an exponentially decaying dynamics ruled by the reversed harmonic oscillator, namely, the simplest irreversible quantum system in the rigged Hilbert spaces. The generalization of this theory to more complex scenarios is still an open question. In this work, we use a thermal third-order medium with an unprecedented giant Kerr coefficient, the M-Cresol/Nylon mixed solution, to access an extremely-nonlinear highly-nonlocal regime and realize anisotropic shock waves. We prove that a superposition of the Gamow vectors in an ad hoc rigged Hilbert space describes the nonlinear beam propagation beyond the shock point. Specifically, the resulting rigged Hilbert space is a tensorial product between the reversed and the standard harmonic oscillators spaces. The anisotropy turns out from the interaction of trapping and antitrapping potentials in perpendicular directions. Our work opens the way to a complete description of novel intriguing shock phenomena, and those mediated by extreme nonlinearities., Comment: 6 pages, 5 figures

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\"odinger 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., Comment: 11 pages, 5 figures

Dispersive shock waves are fascinating phenomena occurring when nonlinearity overwhelms linear effects, such as dispersion and diffraction. Many features of shock waves are still under investigation, as the interplay with noninstantaneity in temporal pulses transmission and nonlocality in spatial beams propagation. Despite the rich and vast literature on nonlinear waves in optical Kerr media, spatial dispersive shock waves in nonlocal materials deserve further attention for their unconventional properties. Indeed, they have been investigated in colloidal matter, chemical physics and biophotonics, for sensing and control of extreme phenomena. Here we review the last developed theoretical models and recent optical experiments on spatial dispersive shock waves in nonlocal media. Moreover, we discuss observations in novel versatile materials relevant for soft matter and biology., Comment: 27 pages, 16 figures, review paper

Novel machine learning computational tools open new perspectives for quantum information systems. Here we adopt the open-source programming library TensorFlow to design multi-level quantum gates including a computing reservoir represented by a random unitary matrix. In optics, the reservoir is a disordered medium or a multi-modal fiber. We show that trainable operators at the input and the readout enable one to realize multi-level gates. We study various qudit gates, including the scaling properties of the algorithms with the size of the reservoir. Despite an initial low slop learning stage, TensorFlow turns out to be an extremely versatile resource for designing gates with complex media, including different models that use spatial light modulators with quantized modulation levels., Comment: Added a new section and a new figure about implementation of the gates by a single spatial light modulator. 9 pages and 4 figures