Planetary cores consist of liquid metals (low Prandtl number Pr) that convect as the core cools. Here, we study nonlinear convection in a rotating (low Ekman number Ek) planetary core using a fully 3D direct numerical simulation. Near the critical thermal forcing (Rayleigh number Ra), convection onsets as thermal Rossby waves, but as Ra increases, this state is superseded by one dominated by advection. At moderate rotation, these states (here called the weak branch and strong branch, respectively) are smoothly connected. As the planetary core rotates faster, the smooth transition is replaced by hysteresis cycles and subcriticality until the weak branch disappears entirely and the strong branch onsets in a turbulent state at $\mathrm{Ek}<{10}^{-6}$. Here, the strong branch persists even as the thermal forcing drops well below the linear onset of convection ($\mathrm{Ra}=0.7{\mathrm{Ra}}_{\mathrm{crit}}$ in this study). We highlight the importance of the Reynolds stress, which is required for convection to subsist below the linear onset. In addition, the Péclet number is consistently above 10 in the strong branch. We further note the presence of a strong zonal flow that is nonetheless unimportant to the convective state. Our study suggests that, in the asymptotic regime of rapid rotation relevant for planetary interiors, thermal convection of liquid metals in a sphere onsets through a subcritical bifurcation.

We report that metallic electrodes are physically pitted during charge transfer events with water droplets or other conductive objects moving in strong electric fields ($>1\mathrm{kV}/\mathrm{cm}$). Post situ microscopic inspection of the electrode shows that an individual charge transfer event yields a crater approximately 1–$3\mu \mathrm{m}$ wide, often with features similar to a splash corona. We interpret the crater formation in terms of localized melting of the electrode via resistive heating concurrent with dielectric breakdown through the surrounding insulating fluid. A scaling analysis indicates that the crater diameter scales as the inverse cube root of the melting point temperature $T_m$ of the metal, in accord with measurements on several metals ($660°\mathrm{C}\le T_m\le 3414°\mathrm{C}$). The process of crater formation provides a possible explanation for the longstanding difficulty in quantitatively corroborating Maxwell’s prediction for the amount of charge acquired by spheres contacting a planar electrode.

Levitating droplets of liquid condensate are known to organize themselves into ordered arrays over hot liquid-gas interfaces. We report experimental observation of similar behavior over a dry heated solid surface. Even though the lifetime of the array is shorter in this case, its geometric characteristics are remarkably similar to the case of droplets levitating over liquid-gas interfaces. A simple model is developed that predicts the mechanisms of both droplet levitation and interdroplet interaction leading to pattern formation over a dry surface; the model is shown to be in good agreement with the experimental data. Using the insights from the new experiments, we are able to resolve some long-standing controversies pertaining to the mechanism of levitation of droplets over liquid-gas interfaces.

We document experimentally four different interactions of a laser-induced bubble and a free-settling particle, with different combinations of the geometric and physical parameters of the system. Our force balance model shows that four nondimensional factors involving the particle radius $a$, the maximum bubble radius $R_{\max }$, the initial separation distance $l_0$ between the particle center and the bubble center, the fluid viscosity ${\mu }_f$, and the particle and fluid densities ${\rho }_p$ and ${\rho }_f$, respectively, in detail $l_0/R_{\max }$, $a/R_{\max }$, ${\rho }_p/{\rho }_f$, and ${\mu }^{*}={\mu }_f{T}_c/{\rho }_f{R}_{\max }^2$, where $T_c=0.915R_{\max }\sqrt{{\rho }_f/(p_{\infty }-p_v)}$, influence the particle-bubble dynamics, and reasonably predict the maximum particle velocity and the limiting condition when the particle starts to “bounce off” the bubble during bubble growth. In particular, we also discover the high-speed ejection of the particle, and a cavity behind the particle, in cases when initially the particle is in very close proximity to the bubble. These observations offer new insights into the causal mechanism for the enhanced cavitation erosion in silt-laden water.

Cluster synchronization is a phenomenon in which a network self-organizes into a pattern of synchronized sets. It has been shown that diverse patterns of stable cluster synchronization can be captured by symmetries of the network. Here, we establish a theoretical basis to divide an arbitrary pattern of symmetry clusters into independently synchronizable cluster sets, in which the synchronization stability of the individual clusters in each set is decoupled from that in all the other sets. Using this framework, we suggest a new approach to find permanently stable chimera states by capturing two or more symmetry clusters—at least one stable and one unstable—that compose the entire fully symmetric network.

We examine the momentum, spin, and orbital angular momentum of structured monochromatic optical fields in dispersive inhomogeneous isotropic media. There are two bifurcations in this general problem: the Abraham-Minkowski dilemma and the kinetic (Poynting-like) versus canonical (spin-orbital) pictures. We show that the kinetic Abraham momentum describes the energy flux and group velocity of the wave in the medium. At the same time, we introduce novel canonical Minkowski-type momentum, spin, and orbital angular momentum densities of the field. These quantities exhibit fairly natural forms, analogous to the Brillouin energy density, as well as multiple advantages as compared with previously considered formalisms. As an example, we apply this general theory to inhomogeneous surface plasmon-polariton (SPP) waves at a metal-vacuum interface and show that SPPs carry a “supermomentum,” proportional to the wave vector $k_p>\omega /c$, and a transverse spin, which can change its sign depending on the frequency $\omega$.

Nonlinear optical interactions of light with materials originate in the microscopic response of the molecular constituents to excitation by an optical field, and are expressed by the first ($\beta$) and second ($\gamma$) hyperpolarizabilities. Upper bounds to these quantities were derived seventeen years ago using approximate, truncated state models that violated completeness and unitarity, and far exceed those achieved by potential optimization of analytical systems. This Letter determines the fundamental limits of the first and second hyperpolarizability tensors using Monte Carlo sampling of energy spectra and transition moments constrained by the diagonal Thomas-Reiche-Kuhn (TRK) sum rules and filtered by the off-diagonal TRK sum rules. The upper bounds of $\beta$ and $\gamma$ are determined from these quantities by applying error-refined extrapolation to perfect compliance with the sum rules. The method yields the largest diagonal component of the hyperpolarizabilities for an arbitrary number of interacting electrons in any number of dimensions. The new method provides design insight to the synthetic chemist and nanophysicist for approaching the limits. This analysis also reveals that the special cases which lead to divergent nonlinearities in the many-state catastrophe are not physically realizable.

To understand the behavior of composite fluid particles such as nucleated cells and double emulsions in flow, we study a finite-size particle encapsulated in a deforming droplet under shear flow as a model system. In addition to its concentric particle-droplet configuration, we numerically explore other eccentric and time-periodic equilibrium solutions, which emerge spontaneously via supercritical pitchfork and Hopf bifurcations. We present the loci of these solutions around the codimension-two point. We adopt a dynamic system approach to model and characterize the coupled behavior of the two bifurcations. By exploring the flow fields and hydrodynamic forces in detail, we identify the role of hydrodynamic particle-droplet interaction which gives rise to these bifurcations.

Many natural and engineering systems are simultaneously subjected to a driving force and a stabilizing force. The interplay between the two forces, especially for highly nonlinear systems such as fluid flow, often results in surprising features. Here we reveal such features in three different types of Rayleigh-Bénard (RB) convection, i.e., buoyancy-driven flow with the fluid density being affected by a scalar field. In the three cases different stabilizing forces are considered, namely (i) horizontal confinement, (ii) rotation around a vertical axis, and (iii) a second stabilizing scalar field. Despite the very different nature of the stabilizing forces and the corresponding equations of motion, at moderate strength we counterintuitively but consistently observe an enhancement in the flux, even though the flow motion is weaker than the original RB flow. The flux enhancement occurs in an intermediate regime in which the stabilizing force is strong enough to alter the flow structures in the bulk to a more organized morphology, yet not too strong to severely suppress the flow motions. Near the optimal transport enhancements all three systems exhibit a transition from a state in which the thermal boundary layer (BL) is nested inside the momentum BL to the one with the thermal BL being thicker than the momentum BL. The observed optimal transport enhancement is explained through an optimal coupling between the suction of hot or fresh fluid and the corresponding scalar fluctuations.

Subradiant excitations, originally predicted by Dicke, have posed a long-standing challenge in physics owing to their weak radiative coupling to environment. Here we engineer massive coherently driven classical subradiance in planar metamaterial arrays as a spatially extended eigenmode comprising over 1000 metamolecules. By comparing the near- and far-field response in large-scale numerical simulations with those in experimental observations we identify strong evidence for classically correlated multimetamolecule subradiant states that dominate the total excitation energy. We show that similar spatially extended many-body subradiance can also exist in plasmonic metamaterial arrays at optical frequencies.

We study the dynamics of a Bose-Einstein condensate in a Sinai-oscillator trap under a monochromatic driving force. Such a trap is formed by a harmonic potential and a repulsive disk located in the center vicinity corresponding to the first experiments of condensate formation by Ketterle and co-workers in 1995. We allow that the external driving allows us to model the regime of weak wave turbulence with the Kolmogorov energy flow from low to high energies. We show that in a certain regime of weak driving and weak nonlinearity such a turbulent energy flow is defeated by the Anderson localization that leads to localization of energy on low energy modes. This is in a drastic contrast to the random phase approximation leading to energy flow to high modes. A critical threshold is determined above which the turbulent flow to high energies becomes possible. We argue that this phenomenon can be studied with ultracold atoms in magneto-optical traps.

It is generally believed that when a single metallic nanowire is sufficiently small, it scatters like a point electric dipole. We show theoretically when a metallic nanowire is placed inside specially designed beams, the magnetic dipole contribution along with the electric dipole resonance can lead to unidirectional scattering in the far field, fulfilling Kerker’s condition. Remarkably, this far-field unidirectional scattering encodes information that is highly dependent on the nanowire’s deflection at a scale much smaller than the wavelength. The special roles of small but essential magnetic response along with the plasmonic resonance are highlighted for this extreme sensitivity as compared with the dielectric counterpart. In addition, the same essential role of the magnetic dipole contribution is also presented for a very small metallic nanosphere.

We study the chaoticity and the predictability of a turbulent flow on the basis of high-resolution direct numerical simulations at different Reynolds numbers. We find that the Lyapunov exponent of turbulence, which measures the exponential separation of two initially close solutions of the Navier-Stokes equations, grows with the Reynolds number of the flow, with an anomalous scaling exponent, larger than the one obtained on dimensional grounds. For large perturbations, the error is transferred to larger, slower scales, where it grows algebraically generating an “inverse cascade” of perturbations in the inertial range. In this regime, our simulations confirm the classical predictions based on closure models of turbulence. We show how to link chaoticity and predictability of a turbulent flow in terms of a finite size extension of the Lyapunov exponent.

In excitable media, chaotic dynamics governed by spiral or scroll waves is often not persistent but transient. Using extensive simulations employing different mathematical models we identify a specific type-II supertransient by an exponential increase of transient lifetimes with the system size in 2D and an investigation of the dynamics (number and lifetime of spiral waves, Kaplan-Yorke dimension). In 3D, simulations exhibit an increase of transient lifetimes and filament lengths only above a critical thickness. Finally, potential implications for understanding cardiac arrhythmias are discussed.

Evaporating droplets are known to show complex motion that has conventionally been explained by the Marangoni effect (flow induced by the gradient of surface tension). Here, we show that the droplet motion can be induced even in the absence of the Marangoni effect due to the gradient of the evaporation rate. We derive an equation for the velocity of a droplet subject to the nonuniform evaporation rate and nonuniform surface tension placed on an inert substrate, where the wettability is uniform and unchanged. The equation explains the previously observed attraction-repulsion-chasing behaviors of evaporating droplets.

We theoretically show that a dipole externally driven by a pulse with a lower-bounded temporal width, and placed inside a cylindrical hollow waveguide, can generate a train of arbitrarily short and focused electromagnetic pulses. The waveguide encloses vacuum with perfect electric conducting walls. A dipole driven by a single short pulse, which is properly engineered to exploit the linear spectral filtering of the cylindrical hollow waveguide, excites longitudinal waveguide modes that are coherently refocused at some particular instances of time, thereby producing arbitrarily short and focused electromagnetic pulses. We numerically show that such ultrafocused pulses persist outside the cylindrical waveguide at distances comparable to its radius.

Band-limited functions can oscillate locally at an arbitrarily fast rate through an interference phenomenon known as superoscillations. Using an optical pulse with a superoscillatory envelope we experimentally break the temporal Fourier-transform focusing limit with a temporal feature that is approximately three times shorter than the duration of a transform-limited Gaussian pulse having a comparable bandwidth while maintaining 30% visibility. We experimentally demonstrate the ability of such signals to achieve temporal superresolution and show numerically in which cases such pulses can outperform transform-limited pulses.