How viscous bubbles collapse: topological and symmetry-breaking instabilities in curvature-driven hydrodynamics

The duality between the mechanical equilibrium of elastic bodies and non-inertial flow in viscous liquids has been a guiding principle in decades of research. However, this duality is broken whenever the Gaussian curvature of viscous films evolves rapidly in time. In such a case the film evolves through a non-inertial yet geometrically nonlinear surface dynamics, which has remained largely unexplored. We reveal the driver of such dynamics as the flow of currents of Gaussian curvature of the evolving surface, rather than the existence of an energetically favored target metric, as in the elastic analogue. Focusing on the prototypical example of a bubble collapsing after rapid depressurization, we show that the spherically-shaped viscous film evolves via a topological instability, whose elastic analogue is forbidden by the existence of an energy barrier. The topological transition brings about compression within the film, triggering another, symmetry breaking instability and radial wrinkles that grow in amplitude and invade the flattening film. Building on an analogy between surface currents of Gaussian curvature and polarization in electrostatic conductors, we obtain quantitative predictions for the evolution of the flattening film. We propose that this classical dynamics has ramifications in the emerging field of viscous hydrodynamics of strongly correlated electrons in two-dimensional materials.Comment: 17 pages (10 pages of appendices). A demonstration video is included as an ancillary fil

Dynamical susceptibility near a long-wavelength critical point with a nonconserved order parameter

We study the dynamic response of a two-dimensional system of itinerant fermions in the vicinity of a uniform ($\mathbf{Q}=0$) Ising nematic quantum critical point of $d-$wave symmetry. The nematic order parameter is not a conserved quantity, and this permits a nonzero value of the fermionic polarization in the $d-$wave channel even for vanishing momentum and finite frequency: $\Pi(\mathbf{q} = 0,\Omega_m) \neq 0$. For weak coupling between the fermions and the nematic order parameter (i.e. the coupling is small compared to the Fermi energy), we perturbatively compute $\Pi (\mathbf{q} = 0,\Omega_m) \neq 0$ over a parametrically broad range of frequencies where the fermionic self-energy $\Sigma (\omega)$ is irrelevant, and use Eliashberg theory to compute $\Pi (\mathbf{q} = 0,\Omega_m)$ in the non-Fermi liquid regime at smaller frequencies, where $\Sigma (\omega) > \omega$. We find that $\Pi(\mathbf{q}=0,\Omega)$ is a constant, plus a frequency dependent correction that goes as $|\Omega|$ at high frequencies, crossing over to $|\Omega|^{1/3}$ at lower frequencies. The $|\Omega|^{1/3}$ scaling holds also in a non-Fermi liquid regime. The non-vanishing of $\Pi (\mathbf{q}=0, \Omega)$ gives rise to additional structure in the imaginary part of the nematic susceptibility $\chi^{''} (\mathbf{q}, \Omega)$ at $\Omega > v_F q$, in marked contrast to the behavior of the susceptibility for a conserved order parameter. This additional structure may be detected in Raman scattering experiments in the $d-$wave geometry

Dynamical susceptibility of a near-critical non-conserved order parameter and B2g Raman response in Fe-based superconductors

We analyze the dynamical response of a two-dimensional system of itinerant fermions coupled to a scalar boson $\phi$, which undergoes a continuous transition towards nematic order with $d-$wave form-factor. We consider two cases: (a) when $\phi$ is a soft collective mode of fermions near a Pomeranchuk instability, and (b) when it is an independent critical degree of freedom, such as a composite spin order parameter near an Ising-nematic transition. In both cases, the order-parameter is not a conserved quantity and the $d-$wave fermionic polarization $\Pi (q, \Omega)$ remains finite even at $q=0$. The polarization $\Pi (0, \Omega)$ has similar behavior in the two cases, but the relations between $\Pi (0, \Omega)$ and the bosonic susceptibility $\chi (0, \Omega)$ are different, leading to different forms of $\chi^{\prime \prime} (0, \Omega)$, as measured by Raman scattering. We compare our results with polarization-resolved Raman data for the Fe-based superconductors FeSe$_{1-x}$S$_x$, NaFe$_{1-x}$Co$_x$As and BaFe$_2$As$_2$. We argue that the data for FeSe$_{1-x}$S$_x$ are well described within Pomeranchuk scenario, while the data for NaFe$_{1-x}$Co$_x$As and BaFe$_2$As$_2$ are better described within the "independent" scenario involving a composite spin order