We study large networks of parametric oscillators as heuristic solvers of random Ising models. In these networks, known as coherent Ising machines, the model to be solved is encoded in the coupling between the oscillators, and a solution is offered by the steady state of the network. This approach relies on the assumption that mode competition steers the network to the ground-state solution of the Ising model. By considering a broad family of frustrated Ising models, we show that the most efficient mode does not correspond generically to the ground state of the Ising model. We infer that networks of parametric oscillators close to threshold are intrinsically not Ising solvers. Nevertheless, the network can find the correct solution if the oscillators are driven sufficiently above threshold, in a regime where nonlinearities play a predominant role. We find that for all probed instances of the model, the network converges to the ground state of the Ising model with a finite probability.

Marangoni instabilities can emerge when a liquid interface is subjected to a concentration or temperature gradient. It is generally believed that for these instabilities bulk effects like buoyancy are negligible compared to interfacial forces, especially on small scales. Consequently, the effect of a stable stratification on the Marangoni instability has hitherto been ignored. Here, however, we show that they can matter. We report, for an immiscible drop immersed in a stably stratified ethanol-water mixture, a new type of oscillatory solutal Marangoni instability that is triggered once the stratification has reached a critical value. We experimentally explore the parameter space spanned by the stratification strength and the drop size and theoretically explain the observed crossover from levitating to bouncing by balancing the advection and diffusion around the drop. Finally, the effect of the stable stratification on the Marangoni instability is surprisingly strongly amplified in confined geometries, leading to an earlier onset.

Symmetry breaking (SB) of fluid-structure interaction problems plays an important role in our understanding of animals’ locomotive and sensing behaviors. In this Letter, we study the SB of flexible filaments clamped at one end and placed in a spanwise periodic array in Stokes flow. The equilibrium state of the filament along the streamwise direction loses stability and experiences two-dimensional and then three-dimensional SBs as the spanwise distance increases, or as the filament rigidity reduces. For slightly deformed filaments, the viscous and pressure forces are commensurate, while for extremely deformed filaments the viscous force becomes dominant.

Synthetic dimension provides a new platform for realizing topological photonic devices. Here, we propose a method to realize a rainbow concentrator of topological photonic states based on the synthetic dimension concept. The synthetic dimension is constructed using a translational degree of freedom of the nanostructures inside the unit cell of a two-dimensional photonic crystal. The translational deformation induces a nontrivial topology in the synthetic dimension, which gives rise to robust interface states at different frequencies. The topological rainbow can trap states with different frequencies, controlled by tuning the spatial modulation of interface state group velocities. The operation frequency as well as the bandwidth of the topological rainbow can be easily tuned by controlling the band gap of the photonic crystal. The topological principle can be applied to photonic crystals of any symmetry and arbitrary material composition, as long as a complete band gap exists. This Letter provides a new and general scheme for the realization of a topological rainbow concentrator and will be useful for the development of topological photonic devices.

We report on a theoretical analysis of the Floquet topological crystalline phases in driven one-dimensional photonic crystals mediated by second-order optical nonlinearity. We define the photonic Berry connection and photonic polarization in such systems using different methods and prove their equivalence. We present two examples of topological phase transitions in which two Floquet bands cross and open new gaps under the driving field. Finally, we analyze the physical consequences of each topological phase transition by examining edge states and filling anomalies. Our study presents routes toward the realization of robust reconfigurable photonic cavities with topologically protected light confinement.

Wavelength selection in reaction-diffusion systems can be understood as a coarsening process that is interrupted by counteracting processes at certain wavelengths. We first show that coarsening in mass-conserving systems is driven by self-amplifying mass transport between neighboring high-density domains. We derive a general coarsening criterion and show that coarsening is generically uninterrupted in two-component systems that conserve mass. The theory is then generalized to study interrupted coarsening and anticoarsening due to weakly broken mass conservation, providing a general path to analyze wavelength selection in pattern formation far from equilibrium.

The transitions between two states of a bistable system are investigated experimentally and analyzed in the framework of rare-event statistics. Considering a disk pendulum swept by a flow in a wind tunnel, bistability between two aerodynamic branches is observed, with spontaneous transitions from one branch to the other. The waiting times before spontaneous transition are distributed following a double exponential as a function of the control parameter, spanning 4 orders of magnitude in time, for both transitions. Inspired by a model originally applied to the transition to turbulence, we show that, for the disk pendulum, both transitions are controlled by rare events of the aerodynamic forces acting on the disk which we propose to link in particular to the vortex shedding-induced fluctuations. Beyond the aerodynamic aspects, this work has interesting fundamental outcomes regarding the broad field of rare events in out-of-equilibrium systems.

Conventionally neutral atmospheric boundary layers (CNBLs), which are characterized with zero surface potential temperature flux and capped by an inversion of potential temperature, are frequently encountered in nature. Therefore, predicting the wind speed profiles of CNBLs is relevant for weather forecasting, climate modeling, and wind energy applications. However, previous attempts to predict the velocity profiles in CNBLs have had limited success due to the complicated interplay between buoyancy, shear, and Coriolis effects. Here, we utilize ideas from the classical Monin-Obukhov similarity theory in combination with a local scaling hypothesis to derive an analytic expression for the stability correction function $\psi =-c_{\psi }(z/L{)}^{1/2}$, where $c_{\psi }=4.2$ is an empirical constant, $z$ is the height above ground, and $L$ is the local Obukhov length based on potential temperature flux at that height, for CNBLs. An analytic expression for this flux is also derived using dimensional analysis and a perturbation method approach. We find that the derived profile agrees excellently with the velocity profile in the entire boundary layer obtained from high-fidelity large eddy simulations of typical CNBLs.

Chimera states have attracted significant attention as symmetry-broken states exhibiting the unexpected coexistence of coherence and incoherence. Despite the valuable insights gained from analyzing specific systems, an understanding of the general physical mechanism underlying the emergence of chimeras is still lacking. Here, we show that many stable chimeras arise because coherence in part of the system is sustained by incoherence in the rest of the system. This mechanism may be regarded as a deterministic analog of noise-induced synchronization and is shown to underlie the emergence of strong chimeras. These are chimera states whose coherent domain is formed by identically synchronized oscillators. Recognizing this mechanism offers a new meaning to the interpretation that chimeras are a natural link between coherence and incoherence.

The scalings of the Rayleigh-Taylor instability are studied numerically for porous media flows when the denser fluid lying on top of the less dense one is also much more viscous. We show that, above a critical value of the viscosity ratio $M$, a symmetry breaking of the buoyancy-driven fingers is observed as they extend much further downward than upward. The asymmetry ratio scales as $M^{1/2}$ while the asymptotic flux across the initial contact line, quantifying the mixing between the two fluids, scales as $M^{-1/2}$. A new fingering mechanism induced by large viscosity contrasts is identified and shows good agreement with experimentally observed dynamics.

The geometric Pancharatnam-Berry (PB) phase not only is of physical interest but also has wide applications ranging from condensed-matter physics to photonics. Space-varying PB phases based on inhomogeneously anisotropic media have previously been used effectively for spin photon manipulation. Here we demonstrate a novel wave-vector-varying PB phase that arises naturally in the transmission and reflection processes in homogeneous media for paraxial beams with small incident angles. The eigenpolarization states of the transmission and reflection processes are determined by the local wave vectors of the incident beam. The small incident angle breaks the rotational symmetry and induces a PB phase that varies linearly with the transverse wave vector, resulting in the photonic spin Hall effect (PSHE). This new PSHE can address the contradiction between spin separation and energy efficiency in the conventional PSHE associated with the Rytov-Vladimirskii-Berry phase, allowing spin photons to be separated completely with a spin separation up to 2.2 times beam waist and a highest energy efficiency of 86%. The spin separation dynamics is visualized by wave coupling equations in a uniaxial crystal, where the centroid positions of the spin photons can be doubled due to the conservation of the angular momentum. Our findings can greatly deepen the understanding in the geometric phase and spin-orbit coupling, paving the way for practical applications of the PSHE.

We present a metastructure architecture with a bistable microstructure that enables extreme broadband frequency conversion. We use numerical and experimental tools to unveil the relationship between input excitations at the unit cell level and output responses at the macrostructural level. We identify soliton-lattice mode resonances resulting in input-independent energy transfer into desired metabeam vibration modes as long as transition waves are triggered within the metastructure. We observe both low-to-high and high-to-low incommensurate frequency interactions in the metabeams, thus enabling energy exchange between bands 2 orders of magnitude apart. This behavior generalizes fluxon-cavity mode resonance in superconducting electronics, providing a general method to extreme frequency conversion in mechanics. Importantly, the introduced architecture allows for expanding the metamaterials design paradigm by fundamentally breaking the dependence of macroscopic dynamics on the unit cell properties. The resulting input-independent nature implies potential applications in broadband frequency regulation and energy transduction.

We solve the advection-diffusion equation for a stochastically stationary passive scalar $\theta$, in conjunction with forced 3D Navier-Stokes equations, using direct numerical simulations in periodic domains of various sizes, the largest being ${8192}^3$. The Taylor-scale Reynolds number varies in the range 140–650 and the Schmidt number $\mathrm{Sc}\equiv \nu /D$ in the range 1–512, where $\nu$ is the kinematic viscosity of the fluid and $D$ is the molecular diffusivity of $\theta$. Our results show that turbulence becomes an ineffective mixer when Sc is large. First, the mean scalar dissipation rate $⟨\chi ⟩=2D⟨|\nabla \theta {|}^2⟩$, when suitably nondimensionalized, decreases as $1/\log \mathrm{Sc}$. Second, 1D cuts through the scalar field indicate increasing density of sharp fronts on larger scales, oscillating with large excursions leading to reduced mixing, and additionally suggesting weakening of scalar variance flux across the scales. The scaling exponents of the scalar structure functions in the inertial-convective range appear to saturate with respect to the moment order and the saturation exponent approaches unity as Sc increases, qualitatively consistent with 1D cuts of the scalar.

We uncover a novel and robust phenomenon that causes the gradual self-replication of spatiotemporal Kerr cavity patterns in cylindrical microresonators. These patterns are inherently synchronized multifrequency combs. Under proper conditions, the axially localized nature of the patterns leads to a fundamental drift instability that induces transitions among patterns with a different number of rows. Self-replications, thus, result in the stepwise addition or removal of individual combs along the cylinder’s axis. Transitions occur in a fully reversible and, consequently, deterministic way. The phenomenon puts forward a novel paradigm for Kerr frequency comb formation and reveals important insights into the physics of multidimensional nonlinear patterns.