It is shown that the addition of small amounts of microscopic rods in a viscous fluid at low Reynolds number causes a significant increase of the flow resistance. Numerical simulations of the dynamics of the solution reveal that this phenomenon is associated to a transition from laminar to chaotic flow. Polymer stresses give rise to flow instabilities which, in turn, perturb the alignment of the rods. This coupled dynamics results in the activation of a wide range of scales, which enhances the mixing efficiency of viscous flows.

A discrete Boltzmann model (DBM) is proposed to probe the Rayleigh-Taylor instability (RTI) in two-component compressible flows. Each species has a flexible specific-heat ratio and is described by one discrete Boltzmann equation (DBE). Independent discrete velocities are adopted for the two DBEs. The collision and force terms in the DBE account for the molecular collision and external force, respectively. Two types of force terms are exploited. In addition to recovering the modified Navier-Stokes equations in the hydrodynamic limit, the DBM has the capability of capturing detailed nonequilibrium effects. Furthermore, we use the DBM to investigate the dynamic process of the RTI. The invariants of tensors for nonequilibrium effects are presented and studied. For low Reynolds numbers, both global nonequilibrium manifestations and the growth rate of the entropy of mixing show three stages (i.e., the reducing, increasing, and then decreasing trends) in the evolution of the RTI. On the other hand, the early reducing tendency is suppressed and even eliminated for high Reynolds numbers. Relevant physical mechanisms are analyzed and discussed.

It is shown that by fitting a Markovian quantum master equation to the numerical solution of the time-dependent Schrödinger equation of a system of two spin-1/2 particles interacting with a bath of up to 34 spin-1/2 particles, the former can describe the dynamics of the two-spin system rather well. The fitting procedure that yields this Markovian quantum master equation accounts for all non-Markovian effects in as much the general structure of this equation allows and yields a description that is incompatible with the Lindblad equation.

In this paper, we give an alternative derivation for the cascaded lattice Boltzmann method (CLBM) within a general multiple-relaxation-time (MRT) framework by introducing a shift matrix. When the shift matrix is a unit matrix, the CLBM degrades into an MRT LBM. Based on this, a consistent forcing scheme is developed for the CLBM. The consistency of the nonslip rule, the second-order convergence rate in space, and the property of isotropy for the consistent forcing scheme is demonstrated through numerical simulations of several canonical problems. Several existing forcing schemes previously used in the CLBM are also examined. The study clarifies the relation between MRT LBM and CLBM under a general framework.

The fluctuation-dissipation relation, a central result in nonequilibrium statistical physics, relates equilibrium fluctuations in a system to its linear response to external forces. Here we provide a direct experimental verification of this relation for viscously coupled oscillators, as realized by a pair of optically trapped colloidal particles. A theoretical analysis, in which interactions mediated by slow viscous flow are represented by nonlocal friction tensors, matches experimental results and reveals a frequency maximum in the amplitude of the mutual response which is a sensitive function of the trap stiffnesses and the friction tensors. This allows for its location and width to be tuned and suggests the utility of the trap setup for accurate two-point microrheology.

We study stochastic thermodynamics for a quantum system of interest whose dynamics is described by a completely positive trace-preserving (CPTP) map as a result of its interaction with a thermal bath. We define CPTP maps with equilibrium as CPTP maps with an invariant state such that the entropy production due to the action of the map on the invariant state vanishes. Thermal maps are a subgroup of CPTP maps with equilibrium. In general, for CPTP maps, the thermodynamic quantities, such as the entropy production or work performed on the system, depend on the combined state of the system plus its environment. We show that these quantities can be written in terms of system properties for maps with equilibrium. The relations that we obtain are valid for arbitrary coupling strengths between the system and the thermal bath. The fluctuations of thermodynamic quantities are considered in the framework of a two-point measurement scheme. We derive the entropy production fluctuation theorem for general maps and a fluctuation relation for the stochastic work on a system that starts in the Gibbs state. Some simplifications for the probability distributions in the case of maps with equilibrium are presented. We illustrate our results by considering spin $\frac{1}{2}$ systems under thermal maps, nonthermal maps with equilibrium, maps with nonequilibrium steady states, and concatenations of them. Finally, and as an important application, we consider a particular limit in which the concatenation of maps generates a continuous time evolution in Lindblad form for the system of interest, and we show that the concept of maps with and without equilibrium translates into Lindblad equations with and without quantum detailed balance, respectively. The consequences for the thermodynamic quantities in this limit are discussed.

We propose a method that makes use of the nonlinear properties of some hypothetical microscopic solid material as the working substance for a microscopic machine. The protocols used are simple (step and elliptic) and allow us to obtain the work and heat exchanged between machine and reservoirs. We calculate the work for a nonlinear single-particle machine that can be treated perturbingly. We obtain the instantaneous work and heat for the machine undergoing cycles that mimic the Carnot and multireservoir protocols. The work calculations are then extended to high values of the nonlinear parameter yielding the quasistatic limit, which is verified numerically. The model we propose is fluctuation driven and we can study in detail its thermostatistics, namely, the work distribution both per cycle and instantaneous and the corresponding fluctuation relations.

Electric pulses in biological cells (action potentials) have been reported to be accompanied by a propagating cell-surface deformation with a nanoscale amplitude. Typically, this cell surface is covered by external layers of polymer material (extracellular matrix, cell wall material, etc.). It was recently demonstrated in excitable plant cells (Chara braunii) that the rigid external layer (cell wall) hinders the underlying deformation. When the cell membrane was separated from the cell wall by osmosis, a mechanical deformation, in the micrometer range, was observed upon excitation of the cell. The underlying mechanism of this mechanical pulse has, to date, remained elusive. Herein we report that Chara cells can undergo a pearling instability, and when the pearled fragments were excited even larger and more regular cell shape changes were observed ($\sim 10–100\mu \mathrm{m}$ in amplitude). These transient cellular deformations were captured by a curvature model that is based on three parameters: surface tension, bending rigidity, and pressure difference across the surface. In this paper these parameters are extracted by curve-fitting to the experimental cellular shapes at rest and during excitation. This is a necessary step to identify the mechanical parameters that change during an action potential.

We explore the phenomenon of the recently discovered inverse transfer of energy from small to large scales in decaying magnetohydrodynamical turbulence by Brandenburg et al. [Phys. Rev. Lett. 114, 075001 (2015)], even for nonhelical magnetic fields. For this investigation we mainly employ the Pencil Code performing a parameter study, where we vary the Prandtl number, the kinematic viscosity, and the initial spectrum. We find that to get a decay that exhibits this inverse transfer, large Reynolds numbers ($O\sim {10}^3$) are needed and low Prandtl numbers of the order unity $\text{Pr}=1$ are preferred. Compared to helical MHD turbulence, though, the inverse transfer is much less efficient in transferring magnetic energy to larger scales than the well-known effect of the inverse cascade. Hence, applying the inverse transfer to the magnetic field evolution in the Early Universe, we question whether the nonhelical inverse transfer is effective enough to explain the observed void magnetic fields if a magnetogenesis scenario during the electroweak phase transition is assumed.

We report shock-wave-like structures that are strikingly different from previously observed fingering instabilities, which occur in a two-layer system of miscible fluids reacting by a second-order reaction $A+B\to S$ in a vertical Hele-Shaw cell. While the traditional analysis expects the occurrence of a diffusion-controlled convection, we show both experimentally and theoretically that the exothermic neutralization reaction can also trigger a wave with a perfectly planar front and nearly discontinuous change in density across the front. This wave propagates fast compared with the characteristic diffusion times and separates the motionless fluid and the area with anomalously intense convective mixing. We explain its mechanism and introduce a new dimensionless parameter, which allows to predict the appearance of such a pattern in other systems. Moreover, we show that our governing equations, taken in the inviscid limit, are formally analogous to well-known shallow-water equations and adiabatic gas flow equations. Based on this analogy, we define the critical velocity for the onset of the shock wave which is found to be in the perfect agreement with the experiments.

We apply recent advances in machine learning and computer vision to a central problem in materials informatics: the statistical representation of microstructural images. We use activations in a pretrained convolutional neural network to provide a high-dimensional characterization of a set of synthetic microstructural images. Next, we use manifold learning to obtain a low-dimensional embedding of this statistical characterization. We show that the low-dimensional embedding extracts the parameters used to generate the images. According to a variety of metrics, the convolutional neural network method yields dramatically better embeddings than the analogous method derived from two-point correlations alone.

We write a master equation for the distributions related to hidden Markov processes (HMPs) and solve it using a functional equation. Thus the solution of HMPs is mapped exactly to the solution of the functional equation. For a general case the latter can be solved only numerically. We derive an exact expression for the entropy of HMPs. Our expression for the entropy is an alternative to the ones given before by the solution of integral equations. The exact solution is possible because actually the model can be considered as a generalized random walk on a one-dimensional strip. While we give the solution for the two second-order matrices, our solution can be easily generalized for the $L$ values of the Markov process and $M$ values of observables: We should be able to solve a system of $L$ functional equations in the space of dimension $M-1$.

We report experimental evidence for a crossover between a liquidlike state and a gaslike state in fluid methane ($\mathrm{C}{\mathrm{H}}_4$). This crossover is observed in all of our experiments, up to a temperature of 397 K, 2.1 times the critical temperature of methane. The crossover has been characterized with both Raman spectroscopy and x-ray diffraction in a number of separate experiments, and confirmed to be reversible. We associate this crossover with the Frenkel line—a recently hypothesized crossover in dynamic properties of fluids extending to arbitrarily high pressure and temperature, dividing the phase diagram into separate regions where the fluid possesses liquidlike and gaslike properties. On the liquidlike side the Raman-active vibration increases in frequency linearly as pressure is increased, as expected due to the repulsive interaction between adjacent molecules. On the gaslike side this competes with the attractive van der Waals potential leading the vibration frequency to decrease as pressure is increased.

Fire ant aggregations are active materials composed of individual constituents that are able to transform internal energy into work. We find using rheology and direct visualization that the aggregation undergoes activity cycles that affect the mechanical properties of the system. When the activity is high, the aggregation approximately equally stores and dissipates energy, it is more homogeneous, and exerts a high outward force. When the activity is low, the aggregation is predominantly elastic, it is more heterogeneous, and it exerts a small outward force. We rationalize our results using a simple kinetic model where the number of active ants within the aggregation is the essential quantity.

The theoretical description of the thermodynamics of water is challenged by the structural transition towards tetrahedral symmetry at ambient conditions. As perturbation theories typically assume a spherically symmetric reference fluid, they are incapable of accurately describing the liquid properties of water at ambient conditions. In this paper we address this problem by introducing the concept of an associated reference perturbation theory (APT). In APT we treat the reference fluid as an associating hard sphere fluid which transitions to tetrahedral symmetry in the fully hydrogen bonded limit. We calculate this transition in a theoretically self-consistent manner without appealing to molecular simulations. This associated reference provides the reference fluid for a second order Barker-Henderson perturbative treatment of the long-range attractions. We demonstrate that this approach gives a significantly improved description of water as compared to standard perturbation theories.

I propose a continuum theory for active nematic gels, defined as fluids or suspensions of orientable rodlike objects endowed with active dynamics, that is based on symmetry arguments and compatibility with thermodynamics. The starting point is our recent theory that models (passive) nematic liquid crystals as relaxing nematic elastomers. The interplay between viscoelastic response and active dynamics of the microscopic constituents is naturally taken into account. By contrast with standard theories, activity is not introduced as an additional term of the stress tensor, but it is added as an external remodeling force that competes with the passive relaxation dynamics and drags the system out of equilibrium. In a simple one-dimensional channel geometry, we show that the interaction between nonuniform nematic order and activity results in either a spontaneous flow of particles or a self-organization into subchannels flowing in opposite directions.

Linear interaction analysis (LIA) is employed to investigate the interaction of reactive and nonreactive shock waves with isotropic vortical turbulence. The analysis is carried out, through Laplace-transform technique, accounting for long-time effects of vortical disturbances on the burnt-gas flow in the fast-reaction limit, where the reaction-region thickness is significantly small in comparison with the most representative turbulent length scales. Results provided by the opposite slow-reaction limit are also recollected. The reactive case is here restricted to situations where the overdriven detonation front does not exhibit self-induced oscillations nor inherent instabilities. The interaction of the planar detonation with a monochromatic pattern of perturbations is addressed first, and then a Fourier superposition for three-dimensional isotropic turbulent fields is employed to provide integral formulas for the amplification of the kinetic energy, enstrophy, and anisotropy downstream. Transitory evolution is also provided for single-frequency disturbances. In addition, further effects associated to the reaction rate, which have not been included in LIA, are studied through direct numerical simulations. The numerical computations, based on WENO-BO4-type scheme, provide spatial profiles of the turbulent structures downstream for four different conditions that include nonreacting shock waves, unstable reacting shock (sufficiently high activation energy), and stable reacting shocks for different detonation thicknesses. Effects of the propagation Mach number, chemical heat release, and burn rate are analyzed.