In the Leidenfrost effect a small drop of fluid is levitated, above a sufficiently hot surface, on a persistent vapor layer generated by evaporation from the drop. The vapor layer thermally insulates the drop from the surface leading to extraordinarily long drop lifetimes. The top-view shape of the levitated drops can exhibit persistent starlike vibrations. I extend recent work [Burton et al., Phys. Rev. Lett. 109, 074301 (2012)] to study the bottom surface of the drop using interference imaging. In this work I use a high-speed camera and automated image analysis to image, locate, and classify the interference fringes. From the interference fringes I reconstruct the shape and height profile of the rim where the drop is closest to the surface. I measure the drop-size dependence of the planar vibrational mode frequencies, which agree well with previous work. I observe a distinct breathing mode in the average radius of the drop, the frequency of which scales differently with drop size than the other modes. This breathing mode can be tightly coupled to a vertical motion of the drop. I further observe a qualitative difference in the structure and dynamics of the vertical profile of the rim between large and small drops.

We analyze both theoretically and experimentally the dynamical behavior of an isothermal axisymmetric liquid bridge close to the minimum-volume stability limit. First, the nature of this stability limit is investigated experimentally by determining the liquid bridge response to a mass force pulse for volumes just above that limit. In our experiments, the liquid bridge breakup takes place only when the critical volume is surpassed and is never triggered by the mass force pulse. Second, the growth of the small-amplitude perturbation mode initiating the liquid bridge breakage is measured experimentally and calculated from the linearized Navier-Stokes equations. The results of the linear stability analysis allow one to explain why liquid bridges with volumes just above the stability limit are so robust. Finally, the nonlinear process leading to the liquid bridge breakup is described from both experimental data and the solution of the full Navier-Stokes equations. Special attention is paid to the free-surface pinchoff. The results show that the flow becomes universal (independent of both the initial and boundary conditions) sufficiently close to that singularity and suggest that the transition from the inviscid to the viscous regime is about to take place in the final stage of both the experiments and numerical simulations.

We propose a modified Poisson-Nernst-Planck (PNP) model to investigate charge transport in electrolytes of inhomogeneous dielectric environment. The model includes the ionic polarization due to the dielectric inhomogeneity and the ion-ion correlation. This is achieved by the self energy of test ions through solving a generalized Debye-Hückel (DH) equation. We develop numerical methods for the system composed of the PNP and DH equations. Particularly, toward the numerical challenge of solving the high-dimensional DH equation, we developed an analytical WKB approximation and a numerical approach based on the selective inversion of sparse matrices. The model and numerical methods are validated by simulating the charge diffusion in electrolytes between two electrodes, for which effects of dielectrics and correlation are investigated by comparing the results with the prediction by the classical PNP theory. We find that, at the length scale of the interface separation comparable to the Bjerrum length, the results of the modified equations are significantly different from the classical PNP predictions mostly due to the dielectric effect. It is also shown that when the ion self energy is in weak or mediate strength, the WKB approximation presents a high accuracy, compared to precise finite-difference results.

We theoretically study energy pumping processes in an electrical circuit with avalanche diodes, where non-Gaussian athermal noise plays a crucial role. We show that a positive amount of energy (work) can be extracted by an external manipulation of the circuit in a cyclic way, even when the system is spatially symmetric. We discuss the properties of the energy pumping process for both quasistatic and finite-time cases, and analytically obtain formulas for the amounts of the work and the power. Our results demonstrate the significance of the non-Gaussianity in energetics of electrical circuits.

We derive the Fokker-Planck equation for multivariable Langevin equations with cross-correlated Gaussian white noises for an arbitrary interpretation of the stochastic differential equation. We formulate the conditions when the solution of the Fokker-Planck equation does not depend on which stochastic calculus is adopted. Further, we derive an equivalent multivariable Ito stochastic differential equation for each possible interpretation of the multivariable Langevin equation. To demonstrate the usefulness and significance of these general results, we consider the motion of Brownian particles. We study in detail the stability conditions for harmonic oscillators with two white noises, one of which is additive, random forcing, and the other, which accounts for fluctuations of either the damping or the spring coefficient, is multiplicative. We analyze the role of cross correlation in terms of the different noise interpretations and confirm the theoretical predictions via numerical simulations. We stress the interest of our results for numerical simulations of stochastic differential equations with an arbitrary interpretation of the stochastic integrals.

We investigate the onset of the discontinuous percolation transition in small-world hyperbolic networks by studying the systems-size scaling of the typical largest cluster approaching the transition, $p\nearrow p_c$. To this end, we determine the average size of the largest cluster $\langle s_{\max }\rangle \sim N^{\Psi \left(p\right)}$ in the thermodynamic limit using real-space renormalization of cluster-generating functions for bond and site percolation in several models of hyperbolic networks that provide exact results. We determine that all our models conform to the recently predicted behavior regarding the growth of the largest cluster, which found diverging, albeit subextensive, clusters spanning the system with finite probability well below $p_c$ and at most quadratic corrections to unity in $\Psi \left(p\right)$ for $p\nearrow p_c$. Our study suggests a large universality in the cluster formation on small-world hyperbolic networks and the potential for an alternative mechanism in the cluster formation dynamics at the onset of discontinuous percolation transitions.

Pair potentials that are bounded at the origin provide an accurate description of the effective interaction for many systems of dissolved soft macromolecules (e.g., flexible dendrimers). Using numerical free-energy calculations, we reconstruct the equilibrium phase diagram of a system of particles interacting through a potential that brings together a Gaussian repulsion with a much weaker Gaussian attraction, close to the thermodynamic stability threshold. Compared to the purely repulsive model, only the reentrant branch of the melting line survives, since for lower densities solidification is overridden by liquid-vapor separation. As a result, the phase diagram of the system recalls that of water up to moderate (i.e., a few tens of MPa) pressures. Upon superimposing a suitable hard core on the double-Gaussian potential, a further transition to a more compact solid phase is induced at high pressure, which might be regarded as the analog of the ice I-to-ice III transition in water.

Agglomeration processes occur in many different realms of science, such as colloid and aerosol formation or formation of bacterial colonies. We study the influence of primary-particle density in agglomerate structures using diffusion-controlled Monte Carlo simulations with realistic space scales through different regimes (diffusion-limited aggregation and diffusion-limited colloid aggregation). The equivalence of Monte Carlo time steps to real time scales is given by Hirsch's hydrodynamical theory of Brownian motion. Agglomerate behavior at different time stages of the simulations suggests that three indices (the fractal exponent, the coordination number, and the eccentricity index) characterize agglomerate geometry. Using these indices, we have found that the initial density of primary particles greatly influences the final structure of the agglomerate, as observed in recent experimental works.

The existence, stability properties, and bifurcation diagram of the nematic umbilical defects is studied. Close to the Fréedericksz transition of nematic liquid crystals with negative anisotropic dielectric constant and homeotropic anchoring, an anisotropic Ginzburg-Landau equation for the amplitude of the tilt of the director away from the vertical axis is derived by taking the three-dimensional (3D) to 2D limit of the Frank-Oseen model. The anisotropic Ginzburg-Landau equation allows us to reveal the mechanism of symmetry breaking of nematic umbilical defects. The positive defect is fully characterized as a function of the anisotropy, while the negative defect is characterized perturbatively. Numerical simulations show quite good agreement with the analytical results.

The interfacial pattern formation problem in an injection-driven radial Hele-Shaw flow is studied for the situation in which a Newtonian fluid of negligible viscosity displaces a viscous non-Newtonian power-law fluid. By utilizing a Darcy-law-like formulation, we tackle the fluid-fluid interface evolution problem perturbatively, and we derive second-order mode-coupling equations that describe the time evolution of the perturbation amplitudes. This allows us to investigate analytically how the non-Newtonian nature of the dislocated fluid determines the morphology of the emerging interfacial patterns. If the pushed fluid is shear-thinning, our results indicate the development of side-branching structures. On the other hand, if the displaced fluid is shear-thickening, one detects the formation of petal-like shapes, markedly characterized by strong tip-splitting events. Finally, a time-dependent injection protocol is presented that is able to restrain finger proliferation via side-branching and tip-splitting. This permits the emergence of symmetric $n$-fold interfacial shapes for which the number of fingers remains fixed as time progresses. This procedure generalizes existing controlling strategies for purely Newtonian flow circumstances to the case of a non-Newtonian, displaced power-law fluid.

We develop a semirelativistic quantum fluid theory based on the expansion of the Dirac Hamiltonian to second order in $1/c$. By making use of the Madelung representation of the wave function, we derive a set of hydrodynamic equations that comprises a continuity equation, an Euler equation for the mean velocity, and an evolution equation for the electron spin density. This hydrodynamic model is then applied to study the dynamics of a dense and weakly relativistic electron plasma. In particular, we investigate the impact of the quantum-relativistic spin effects on the Faraday rotation in a one-dimensional plasma slab irradiated by an x-ray laser source.

We investigate subpicosecond dynamics of warm dense hydrogen at the XUV free-electron laser facility (FLASH) at DESY (Hamburg). Ultrafast impulsive electron heating is initiated by a $\le 300$-fs short x-ray burst of 92-eV photon energy. A second pulse probes the sample via x-ray scattering at jitter-free variable time delay. We show that the initial molecular structure dissociates within $(0.9\pm 0.2)$ ps, allowing us to infer the energy transfer rate between electrons and ions. We evaluate Saha and Thomas-Fermi ionization models in radiation hydrodynamics simulations, predicting plasma parameters that are subsequently used to calculate the static structure factor. A conductivity model for partially ionized plasma is validated by two-temperature density-functional theory coupled to molecular dynamic simulations and agrees with the experimental data. Our results provide important insights and the needed experimental data on transport properties of dense plasmas.

Generalized Bloch wave functions of bulk structures, which are composed of not only propagating waves but also decaying and growing evanescent waves, are known to be essential for defining the open boundary conditions in the calculations of the electronic surface states and scattering wave functions of surface and junction structures. Electronic complex band structures being derived from the generalized Bloch wave functions are also essential for studying bound states of the surface and junction structures, which do not appear in conventional band structures. We present a novel calculation method to obtain the generalized Bloch wave functions of periodic bulk structures by solving a generalized eigenvalue problem, whose dimension is drastically reduced in comparison with the conventional generalized eigenvalue problem derived by Fujimoto and Hirose [Phys. Rev. B 67, 195315 (2003)]. The generalized eigenvalue problem derived in this work is even mathematically equivalent to the conventional one, and, thus, we reduce computational cost for solving the eigenvalue problem considerably without any approximation and losing the strictness of the formulations. To exhibit the performance of the present method, we demonstrate practical calculations of electronic complex band structures and electron transport properties of Al and Cu nanoscale systems. Moreover, employing atom-structured electrodes and jellium-approximated ones for both of the Al and Si monatomic chains, we investigate how much the electron transport properties are unphysically affected by the jellium parts.

Nonreciprocal shape deformations can drive inertialess cellular swimming, as first explored by Taylor and Lighthill in the 1950s, for the small-amplitude squirming of a planar and a spherical surface, respectively. Lighthill's squirmer, in particular, has been extensively studied for large wave numbers in the context of ciliated microbes. The maximal power efficiency for small-amplitude planar squirming motility is well characterized and degenerate, with nonunique optimal swimming strokes. We explicitly show that this degeneracy is retained at high wave numbers for the small-amplitude spherical squirmer such as a ciliated microbe but is broken for low wave numbers. Hence further complexity emerges in parameter regimes outside that of ciliate swimming even at small amplitudes. Large-amplitude squirming also characterizes more recent observations of large-amplitude/low-wave-number membrane deformations driving the motility of Euglena, neutrophils, and Dictyostelium discoideum. Thus boundary element numerical methods are used to explore swimming with increased deformation amplitudes, especially in the context of power efficiency and swimming performance. As radial squirming amplitudes are increased, small-amplitude linearized theories can be unreliable even for nominally low deformation amplitudes. Furthermore, even for a simple single-mode metachronal wave, a highly motile and efficient large-deformation/small-wave-number swimming modality arises, which can surpass theoretical limitations of purely tangential squirming given a constrained surface deformation velocity.

We present analytic and numeric results for percolation in a network formed of interdependent spatially embedded networks. We show results for a treelike and a random regular network of networks each with (i) unconstrained dependency links and (ii) dependency links restricted to a maximum Euclidean length $r$. Analytic results are given for each network of networks with spatially unconstrained dependency links and compared to simulations. For the case of two fully interdependent spatially embedded networks it was found [Li et al., Phys. Rev. Lett. 108, 228702 (2012)] that the system undergoes a first-order phase transition only for $r>r_c\approx 8$. We find here that for treelike networks of networks (composed of $n$ networks) $r_c$ significantly decreases as $n$ increases and rapidly $\left(n\ge 11\right)$ reaches its limiting value of 1. For cases where the dependencies form loops, such as in random regular networks, we show analytically and confirm through simulations that there is a certain fraction of dependent nodes, $q_{\max }$, above which the entire network structure collapses even if a single node is removed. The value of $q_{\max }$ decreases quickly with $m$, the degree of the random regular network of networks. Our results show the extreme sensitivity of coupled spatial networks and emphasize the susceptibility of these networks to sudden collapse. The theory proposed here requires only numerical knowledge about the percolation behavior of a single network and therefore can be used to find the robustness of any network of networks where the profile of percolation of a singe network is known numerically.

Spatiotemporal patterns, indicating the spatiotemporal variability of individual abundance, are a pronounced scenario in ecological interactions. Most of the existing models for spatiotemporal patterns treat species as homogeneous groups of individuals with average characteristics by ignoring intraspecific physiological variations at the individual level. Here we explore the impacts of size variation within species resulting from individual ontogeny, on the emergence of spatiotemporal patterns in a fully size-structured population model. We found that size dependency of animal's diffusivity greatly promotes the formation of spatiotemporal patterns, by creating regular spatiotemporal patterns out of temporal chaos. We also found that size-dependent diffusion can substitute large-amplitude base harmonics with spatiotemporal patterns with lower amplitude oscillations but with enriched harmonics. Finally, we found that the single-generation cycle is more likely to drive spatiotemporal patterns compared to predator-prey cycles, meaning that the mechanism of Hopf bifurcation might be more common than hitherto appreciated since the former cycle is more widespread than the latter in case of interacting populations. Due to the ubiquity of individual ontogeny in natural ecosystems we conclude that diffusion variability within populations is a significant driving force for the emergence of spatiotemporal patterns. Our results offer a perspective on self-organized phenomena, and pave a way to understand such phenomena in systems organized as complex ecological networks.

One of the most efficient actuation methods of robotic microswimmers for biomedical applications is by applying time-varying external magnetic fields. In order to improve the design of the swimmer and optimize its performance, one needs to develop simple theoretical models that enable explicit analysis of the swimmer's dynamics. This paper studies the dynamics of a simple microswimmer model with two magnetized links connected by an elastic joint, which undergoes planar undulations induced by an oscillating magnetic field. The nonlinear dynamics of the microswimmer is formulated by assuming Stokes flow and using resistive force theory to calculate the viscous drag forces. Key effects that enable the swimmer to overcome the scallop theorem and generate net propulsion are identified, including violation of front-back symmetry. Assuming small oscillation amplitude, approximate solution is derived by using perturbation expansion, and leading-order expressions for the swimmer's displacement per cycle $X$ and average speed $V$ are obtained. Optimal actuation frequencies that maximize $X$ or $V$ are found for given swimmer's parameters. An ultimate optimal choice of swimmer's parameters and actuation frequency is found, for which the average swimming speed $V$ attains a global maximum. Finally, the theoretical predictions of optimal performance values are validated by comparison to reported experimental results of magnetic microswimmers.