Why does a microwave oven work? How does biological tissue absorb electromagnetic radiation? Astonishingly, we do not have a definite answer to these simple questions because the microscopic processes governing the absorption of electromagnetic waves by water are largely unclarified. This absorption can be quantified by dielectric loss spectra, which reveal a huge peak at a frequency of the exciting electric field of about 20 GHz and a gradual tailing off toward higher frequencies. The microscopic interpretation of such spectra is highly controversial and various superpositions of relaxation and resonance processes ascribed to single-molecule or molecule-cluster motions have been proposed for their analysis. By combining dielectric, microwave, THz, and far-infrared spectroscopy, here we provide nearly continuous temperature-dependent broadband spectra of water. Moreover, we find that corresponding spectra for aqueous solutions reveal the same features as pure water. However, in contrast to the latter, crystallization in these solutions can be avoided by supercooling. As different spectral contributions tend to disentangle at low temperatures, this enables us to deconvolute them when approaching the glass transition under cooling. We find that the overall spectral development, including the 20 GHz feature (employed for microwave heating), closely resembles the behavior known for common supercooled liquids. Thus water's absorption of electromagnetic waves at room temperature is not unusual but very similar to that of glass-forming liquids at elevated temperatures, deep in the low-viscosity liquid regime, and should be interpreted along similar lines.

We analyzed the optical band gaps for axially propagating electromagnetic waves throughout a metallic doped cholesteric elastomer. The composed medium is made of metallic nanospheres (silver) randomly dispersed in a cholesteric elastomer liquid crystal whose dielectric properties can be represented by a resonant effective uniaxial tensor. We found that the band gap properties of the periodic system greatly depend on the volume fraction of nanoparticles in the cholesteric elastomer. In particular, we observed a displacement of the reflection band for quite small fraction volumes whereas for larger values of this fraction there appears a secondary band in the higher frequency region. We also have calculated the transmittance and reflectance spectra for our system. These calculations verify the mentioned band structure and provide additional information about the polarization features of the radiation.

The study of elastic membranes carrying topological defects has a longstanding history, going back at least to the 1950s. When allowed to buckle in three-dimensional space, membranes with defects can totally relieve their in-plane strain, remaining with a bending energy, whose rigidity modulus is small compared to the stretching modulus. In this paper we study membranes with a single edge dislocation. We prove that the minimum bending energy associated with strain-free configurations diverges logarithmically with the size of the system.

We present results from extensive Monte Carlo simulations of polymer models where each lattice site can be visited by up to $K$ monomers and no restriction is imposed on the number of bonds on each lattice edge. These multiple monomer per site (MMS) models are investigated on the square and cubic lattices, for $K=2$ and 3, by associating Boltzmann weights ${\omega }_0=1$, ${\omega }_1=e^{{\beta }_1}$, and ${\omega }_2=e^{{\beta }_2}$ to sites visited by 1, 2, and 3 monomers, respectively. Two versions of the MMS models are considered for which immediate reversals of the walks are allowed (RA) or forbidden (RF). In contrast to previous simulations of these models, we find the same thermodynamic behavior for both RA and RF versions. In three dimensions, the phase diagrams, in space ${\beta }_2\times {\beta }_1$, are featured by coil and globule phases separated by a line of $\mathrm{\Theta }$ points, as thoroughly demonstrated by the metric ${\nu }_t$, crossover ${\varphi }_t$, and entropic ${\gamma }_t$ exponents. The existence of the $\mathrm{\Theta }$ lines is also confirmed by the second virial coefficient. This shows that no discontinuous collapse transition exists in these models, in contrast to previous claims based on a weak bimodality observed in some distributions, which indeed exists in a narrow region very close to the $\mathrm{\Theta }$ line when ${\beta }_1<0$. Interestingly, in two dimensions, only a crossover is found between the coil and globule phases.

In this work we revisit the nonmonotonic behavior (NMB) of synchronization time with velocity reported for systems of mobile pulse-coupled oscillators (PCOs). We devise a control parameter that allows us to predict in which range of velocities NMB may occur, also uncovering the conditions allowing us to establish the emergence of NMB based on specific features of the connectivity rule. Specifically, our results show that if the connectivity rule is such that the interaction patterns are sparse and, more importantly, include a large fraction of nonreciprocal interactions, then the system will display NMB. We furthermore provide a microscopic explanation relating the presence of such features of the connectivity patterns to the existence of local clusters unable to synchronize, termed frustrated clusters, for which we also give a precise definition in terms of simple graph concepts. We conclude that, if the probability of finding a frustrated cluster in a system of moving PCOs is high enough, NMB occurs in a predictable range of velocities.

We study the influence of receptor-ligand interactions on the force response of single polymer chains theoretically. The extension of the chain is modeled in terms of freely jointed chain or elastic freely jointed chain (EFJC) models. The situation involving noninteracting bonds is solved exactly, while effects of interactions are treated within a mean-field approximation. The form with shorter bonds governs the low force situation, while the form with longer bonds is relevant in the high force regime. We further discuss the accuracy of approximate relations, which were used to describe the response of the EFJC model.

We obtain several exact results for universal distributions involving the maximum of the ${\mathrm{Airy}}_2$ process minus a parabola and plus a Brownian motion, with applications to the one-dimensional Kardar-Parisi-Zhang (KPZ) stochastic growth universality class. This allows one to obtain (i) the universal limit, for large time separation, of the two-time height correlation for droplet initial conditions, e.g., $C_{\infty }={lim}_{t_2/t_1\to +\infty }{\overline{h\left(t_1\right)h\left(t_2\right)}}^c/{\overline{h{\left(t_1\right)}^2}}^c$, with $C_{\infty }\approx 0.623$, as well as conditional moments, which quantify ergodicity breaking in the time evolution; (ii) in the same limit, the distribution of the midpoint position $x\left(t_1\right)$ of a directed polymer of length $t_2$; and (iii) the height distribution in stationary KPZ with a step. These results are derived from the replica Bethe ansatz for the KPZ continuum equation, with a “decoupling assumption” in the large time limit. They agree and confirm, whenever they can be compared, with (i) our recent tail results for two-time KPZ with the work by de Nardis and Le Doussal [J. Stat. Mech. (2017) 053212], checked in experiments with the work by Takeuchi and co-workers [De Nardis et al., Phys. Rev. Lett. 118, 125701 (2017)] and (ii) a recent result of Maes and Thiery [J. Stat. Phys. 168, 937 (2017)] on midpoint position.

Micron-sized particles (microbeads) dispersed in a suspension of magnetic nanoparticles, i.e., ferrofluids, can be assembled into different types of structures upon application of an external magnetic field. This paper is devoted to theoretical modeling of a relative motion of a pair of microbeads (either soft ferromagnetic or diamagnetic) in the ferrofluid under the action of applied uniform magnetic field which induces magnetic moments in the microbeads making them attracting to each other. The model is based on a point-dipole approximation for the magnetic interactions between microbeads mediated by the ferrofluid; however, the ferrofluid is considered to possess an anisotropic magnetic permeability thanks to field-induced structuring of its nanoparticles. The model is tested against experimental results and shows generally better agreement with experiments than the model considering isotropic magnetic permeability of ferrofluids. The results could be useful for understanding kinetics of aggregation of microbeads suspended in a ferrofluid. From a broader perspective, the present study is believed to contribute to a general understanding of particle behaviors in anisotropic media.

Using experiments and simulations, we study the flow of soft particles through quasi-two-dimensional hoppers. The first experiment uses oil-in-water emulsion droplets in a thin sample chamber. Due to surfactants coating the droplets, they easily slide past each other, approximating soft frictionless disks. For these droplets, clogging at the hopper exit requires a narrow hopper opening only slightly larger than the droplet diameter. The second experiment uses soft hydrogel particles in a thin sample chamber, where we vary gravity by changing the tilt angle of the chamber. For reduced gravity, clogging becomes easier and can occur for larger hopper openings. Our simulations mimic the emulsion experiments and demonstrate that softness is a key factor controlling clogging: with stiffer particles or a weaker gravitational force, clogging is easier. The fractional amount a single particle is deformed under its own weight is a useful parameter measuring particle softness. Data from the simulation and hydrogel experiments collapse when compared using this parameter. Our results suggest that prior studies using hard particles were in a limit where the role of softness is negligible, which causes clogging to occur with significantly larger openings.

Theory and practical implementations for wake-free propulsion systems are proposed and proven with computational fluid dynamic modeling. Introduced earlier, the concept of active hydrodynamic metamaterials is advanced by introducing magnetohydrodynamic metamaterials, structures with custom-designed volumetric distribution of Lorentz forces acting on a conducting fluid. Distributions of volume forces leading to wake-free, laminar flows are designed using multivariate optimization. Theoretical indications are presented that such flows can be sustained at arbitrarily high Reynolds numbers. Moreover, it is shown that in the limit $\mathrm{Re}\gg {10}^2$, a fixed volume force distribution may lead to a forced laminar flow across a wide range of Re numbers, without the need to reconfigure the force-generating metamaterial. Power requirements for such a device are studied as a function of the fluid conductivity. Implications to the design of distributed propulsion systems underwater and in space are discussed.

Flow of particles of two different species through a narrow channel with solely two discrete spatial positions is analyzed with respect to the species' capability to cooperate or compete for transport. The origin of the latter arises from particle-channel and interparticle interactions within the channel, i.e., blocking the position of a particle, and its interaction with its neighbors in the channel. The variety of occupation options within the channel defines the state space. The transition dynamics within is considered as a continuous Markov process, which, in contrast to mean-field approaches, conserves explicitly spatial correlations. A strong repulsive interaction between particles of the same kind and a very attractive empty channel imply a strong entanglement of transport of both species. In the limiting case of perfect coupling, flows in state space are restricted to a cyclic subspace, where they become all equivalent in the steady state. In particular, this implies equal particle flows of the two species. Entanglement of transport implies that the species mutually exert entropic forces on each other. For parallel directed concentration gradients this implies that the species' ability to cooperate increases with the degree of entanglement. Thus, the gradient of one species reciprocally induces a higher flow of the other species when compared to that in its absence. The opposite holds for antiparallel gradients where species mutually hamper their transport. For a sufficient strong coupling, the species under the influence of the stronger concentration gradient drives the other against its gradient, i.e., the positive mixing entropy production of the driving species becomes the motor for the negative mixing entropy production of the driven one. The degree of effectiveness by which negative entropy production emerges at the cost of positive entropy production increases with the coupling strength. This becomes evident from location and connectivity of the sources of entropy production in state space.

We study the entropy production of systems out of equilibrium using random networks. We focus on systems with a finite number of states described by a master equation close to equilibrium. The dynamics are mapped into a network of states (nodes) connected by transition rates (links). Using this framework, we analyze the entropy production of ensembles of randomly generated networks owing to specific constraints (e.g., size or symmetries) and identify the most important parameters that determine its value. This analysis gives a null-model estimation for the entropy production that can be used for comparison with specific systems.

Motivated by identifying the origin of the bias dependence of tension propagation, we investigate methods for measuring tension propagation quantitatively in computer simulations of driven polymer translocation. Here, the motion of flexible polymer chains through a narrow pore is simulated using Langevin dynamics. We measure tension forces, bead velocities, bead distances, and bond angles along the polymer at all stages of translocation with unprecedented precision. Measurements are done at a standard temperature used in simulations and at zero temperature to pin down the effect of fluctuations. The measured quantities were found to give qualitatively similar characteristics, but the bias dependence could be determined only using tension force. We find that in the scaling relation $\tau \sim N^{\beta }{f}_d^{\alpha }$ for translocation time $\tau$, the polymer length $N$, and the bias force $f_d$, the increase of the exponent $\beta$ with bias is caused by center-of-mass diffusion of the polymer toward the pore on the cis side. We find that this diffusion also causes the exponent $\alpha$ to deviate from the ideal value $-1$. The bias dependence of $\beta$ was found to result from combination of diffusion and pore friction and so be relevant for polymers that are too short to be considered asymptotically long. The effect is relevant in experiments all of which are made using polymers whose lengths are far below the asymptotic limit. Thereby, our results also corroborate the theoretical prediction by Sakaue's theory [Polymers 8, 424 (2016)] that there should not be bias dependence of $\beta$ for asymptotically long polymers. By excluding fluctuations we also show that monomer crowding at the pore exit cannot have a measurable effect on translocation dynamics under realistic conditions.

A proposal to link the equation of state of a monocomponent hard-disk fluid to the equation of state of a polydisperse hard-disk mixture is presented. Event-driven molecular dynamics simulations are performed to obtain data for the compressibility factor of the monocomponent fluid and of 26 polydisperse mixtures with different size distributions. Those data are used to assess the proposal and to infer the values of the compressibility factor of the monocomponent hard-disk fluid in the metastable region from those of mixtures in the high-density region. The collapse of the curves for the different mixtures is excellent in the stable region. In the metastable regime, except for two mixtures in which crystallization is present, the outcome of the approach exhibits a rather good performance. The simulation results indicate that a (reduced) variance of the size distribution larger than about 0.01 is sufficient to avoid crystallization and explore the metastable fluid branch.

The macroscopic description of buoyancy-driven thermal convection in porous media is governed by advection-diffusion processes, which in the presence of thermophysical heterogeneities fail to predict the onset of thermal convection and the average rate of heat transfer. This work extends the classical model of heat transfer in porous media by including a fractional-order advective-dispersive term to account for the role of thermophysical heterogeneities in shifting the thermal instability point. The proposed fractional-order model overcomes limitations of the common closure approaches for the thermal dispersion term by replacing the diffusive assumption with a fractional-order model. Through a linear stability analysis and Galerkin procedure, we derive an analytical formula for the critical Rayleigh number as a function of the fractional model parameters. The resulting critical Rayleigh number reduces to the classical value in the absence of thermophysical heterogeneities when solid and fluid phases have similar thermal conductivities. Numerical simulations of the coupled flow equation with the fractional-order energy model near the primary bifurcation point confirm our analytical results. Moreover, data from pore-scale simulations are used to examine the potential of the proposed fractional-order model in predicting the amount of heat transfer across the porous enclosure. The linear stability and numerical results show that, unlike the classical thermal advection-dispersion models, the fractional-order model captures the advance and delay in the onset of convection in porous media and provides correct scalings for the average heat transfer in a thermophysically heterogeneous medium.

The role of an external field on capillary waves at the liquid-vapor interface of a dipolar fluid is investigated using molecular dynamics simulations. For fields parallel to the interface, the interfacial width squared increases linearly with respect to the logarithm of the size of the interface across all field strengths tested. The value of the slope decreases with increasing field strength, indicating that the field dampens the capillary waves. With the inclusion of the parallel field, the surface stiffness increases with increasing field strength faster than the surface tension. For fields perpendicular to the interface, the interfacial width squared is linear with respect to the logarithm of the size of the interface for small field strengths, and the surface stiffness is less than the surface tension. Above a critical field strength that decreases as the size of the interface increases, the interface becomes unstable due to the increased amplitude of the capillary waves.

Quantum phase transitions of a $q$-state Potts model in fractal lattices are studied using a continuous-time quantum Monte Carlo simulation technique. For small values of $q$, the transition is found to be second order and critical exponents of the quantum critical point are calculated. The dynamic critical exponent $z$ is found to be greater than one for all fractals studied, which is in contrast to integer-dimensional regular lattices. When $q$ is greater than a certain value $q_c$, the phase transition becomes first order, where $q_c$ depends on the lattice. Further analysis shows that the characteristics of phase transitions are more sensitive to the average number of nearest neighbors than the Hausdorff dimension or the order of ramification.