We experimentally study quasi-two-dimensional dilute granular flow around intruders whose shape, size, and relative impact speed are systematically varied. Direct measurement of the flow field reveals that three in-principle independent measurements of the nonuniformity of the flow field are in fact all linearly related: (1) granular temperature, (2) flow-field divergence, and (3) shear-strain rate. The shock front is defined as the local maxima in each of these measurements. The shape of the shock front is well described by an inverted catenary and is driven by the formation of a dynamic arch during steady flow. We find universality in the functional form of the shock front within the range of experimental values probed. Changing the intruder size, concavity, and impact speed only results in a scaling and shifting of the shock front. We independently measure the horizontal lift force on the intruder and find that it can be understood as a result of the interplay between the shock profile and the intruder shape.

The triumph of heat engines is their ability to convert the disordered energy of thermal sources into useful mechanical motion. In recent years, much effort has been devoted to generalizing thermodynamic notions to the quantum regime, partly motivated by the promise of surpassing classical heat engines. Here, we instead adopt a bottom-up approach: we propose a realistic autonomous heat engine that can serve as a test bed for quantum effects in the context of thermodynamics. Our model draws inspiration from actual piston engines and is built from closed-system Hamiltonians and weak bath coupling terms. We analytically derive the performance of the engine in the classical regime via a set of nonlinear Langevin equations. In the quantum case, we perform numerical simulations of the master equation. Finally, we perform a dynamic and thermodynamic analysis of the engine's behavior for several parameter regimes in both the classical and quantum case and find that the latter exhibits a consistently lower efficiency due to additional noise.

This paper shows that, in a pair of one-dimensional complex Ginzburg-Landau (CGL) systems, diffusive connections can induce amplitude death. Stability analysis of a spatially uniform steady state in coupled CGL systems reveals that amplitude death never occurs in a pair of identical CGL systems coupled by no-delay connection, but can occur in the case of delay connection. Moreover, amplitude death never occurs in coupled identical CGL systems with zero nominal frequency. Based on these analytical results, we propose a procedure for designing the connection delay time and the coupling strength to induce spatial-robust stabilization, that is, a stabilization of the steady state for any system size and any boundary condition. Numerical simulations are performed to confirm the analytical results.

We illustrate a counterintuitive effect of an additive stochastic force, which acts independently on each element of an ensemble of globally coupled oscillators. We show that a very small white noise not only broadens the clusters, wherever they are induced by the deterministic forces, but can also stabilize a linearly unstable collective periodic regime: self-consistent partial synchrony. With the help of microscopic simulations we are able to identify two noise-induced bifurcations. A macroscopic analysis, based on a perturbative solution of the associated nonlinear Fokker-Planck equation, confirms the numerical studies and allows determining the eigenvalues of the stability problem. We finally argue about the generality of the phenomenon.

We investigate how leaders emerge as a consequence of the competitive dynamics between coupled papers in a model citation network. Every paper is allocated an initial fitness depending on its intrinsic quality. Its fitness then involves dynamically as a consequence of the competition between itself and all the other papers in the field. It picks up citations as a result of this adaptive dynamics, becoming a leader if it has the highest citation count at a given time. Extensive analytical and numerical investigations of this model suggest the existence of a universal phase diagram, divided into regions of weak and strong coupling. In the former, we find an “extended” and rather structureless distribution of citation counts among many fit papers; leaders are not necessarily those with the maximal fitness at any given time. By contrast, the strong-coupling region is characterized by a strongly hierarchical distribution of citation counts, that are “localized” among only a few extremely fit papers, and exhibit strong history-to-history fluctuations, as a result of the complex dynamics among papers in the tail of the fitness distribution.

Proper neural connectivity inference has become essential for understanding cognitive processes associated with human brain function. Its efficacy is often hampered by the curse of dimensionality. In the electroencephalogram case, which is a noninvasive electrophysiological monitoring technique to record electrical activity of the brain, a possible way around this is to replace multichannel electrode information with dipole reconstructed data. We use a method based on maximum entropy and the renormalization group to infer the position of the sources, whose success hinges on transmitting information from low- to high-resolution representations of the cortex. The performance of this method compares favorably to other available source inference algorithms, which are ranked here in terms of their performance with respect to directed connectivity inference by using artificially generated dynamic data. We examine some representative scenarios comprising different numbers of dynamically connected dipoles over distinct cortical surface positions and under different sensor noise impairment levels. The overall conclusion is that inverse problem solutions do not affect the correct inference of the direction of the flow of information as long as the equivalent dipole sources are correctly found.

Many cells exploit the bending or rotation of flagellar filaments in order to self-propel in viscous fluids. While appropriate theoretical modeling is available to capture flagella locomotion in simple, Newtonian fluids, formidable computations are required to address theoretically their locomotion in complex, nonlinear fluids, e.g., mucus. Based on experimental measurements for the motion of rigid rods in non-Newtonian fluids and on the classical Carreau fluid model, we propose empirical extensions of the classical Newtonian resistive-force theory to model the waving of slender filaments in non-Newtonian fluids. By assuming the flow near the flagellum to be locally Newtonian, we propose a self-consistent way to estimate the typical shear rate in the fluid, which we then use to construct correction factors to the Newtonian local drag coefficients. The resulting non-Newtonian resistive-force theory, while empirical, is consistent with the Newtonian limit, and with the experiments. We then use our models to address waving locomotion in non-Newtonian fluids and show that the resulting swimming speeds are systematically lowered, a result which we are able to capture asymptotically and to interpret physically. An application of the models to recent experimental results on the locomotion of Caenorhabditis elegans in polymeric solutions shows reasonable agreement and thus captures the main physics of swimming in shear-thinning fluids.

We investigate experimentally the influence of a fixed obstacle on gas rising in a dense suspension. Air is injected at a constant flow rate by a single nozzle at the bottom center of a Hele-Shaw cell. Without obstacles, previous works have shown that a fluidized zone is formed with a parabolic shape, with a central air channel and two granular convection rolls on its sides. Here, we quantify the influence of the obstacle's shape, size, and height on the location and dynamics of the central air channel. Different regimes are reported: the air channel can simply deviate (stable), or it can switch sides over time (unstable), leading to two signatures not only above the obstacle, but sometimes also below it. This feedback also influences the channel deviation when bypassing the obstacle. A wake of less or no motion is reported above the largest obstacles as well as the maximum probability of gas location, which can be interesting for practical applications. The existence of a critical height $h_c\simeq 7$ cm is discussed and compared with the existence of an air finger that develops from the injection nozzle and is stable in time. A dimensionless number describing the transition between air fingering and fracturing makes it possible to predict the channel's stability.

This article describes a phase-field model for an isothermal multicomponent, multiphase system which avoids implicit interfacial energy contributions by starting from a grand potential formulation. A method is developed for incorporating arbitrary forms of the equilibrium thermodynamic potentials in all phases to determine an explicit relationship between chemical potentials and species concentrations. The model incorporates variable densities between adjacent phases, defect migration, and dependence of internal pressure on object dimensions ranging from the macro- to nanoscale. A demonstrative simulation of an overpressurized nanoscopic intragranular bubble in nuclear fuel migrating to a grain boundary under kinetically limited vacancy diffusion is shown.

Applications of the transient Boltzmann transport equation (BTE) have undergone much investigation, such as radiative heat transfer and neutron transport. This paper provides a lattice Boltzmann model to efficiently resolve the multidimensional transient BTE. For a higher angular resolution, enough transport directions are considered while the transient BTE in each direction is treated as a conservation law equation and solved independently. Both macroscopic equations recovered from a Chapman-Enskog expansion and simulated results of typical benchmark problems show not only the second-order accuracy but also the flexibility and applicability of the proposed lattice Boltzmann model. This approach may contribute a powerful technique for the parallel simulation of large-scale engineering and some alternative perspectives for solving the nonlinear transport problem further.

We consider the large-time dynamics of one-dimensional processes involving adsorption and desorption of extended hard-core particles (dimers, trimers, ..., $k$-mers), while interacting through their constituent monomers. Desorption can occur whether or not these latter adsorbed together, which leads to reconstitution of $k$-mers and the appearance of sectors of motion with nonlocal conservation laws for $k\ge 3$. Dynamic exponents of the sector including the empty chain are evaluated by finite-size scaling analyses of the relaxation times embodied in the spectral gaps of evolution operators. For attractive interactions it is found that in the low-temperature limit such time scales converge to those of the Glauber dynamics, thus suggesting a diffusive universality class for $k\ge 2$. This is also tested by simulated quenches down to $T=0$, where a common scaling function emerges. By contrast, under repulsive interactions the low-temperature dynamics is characterized by metastable states which decay subdiffusively to a highly degenerate and partially jammed phase.

We study two identical FitzHugh-Nagumo oscillators which are coupled with one or two different time delays. If only a single-delay coupling is used, the length of the delay determines whether the synchronization manifold is transversally stable or unstable, exhibiting mixed-mode or chaotic oscillations in which the small amplitude oscillations are always in phase but the large amplitude oscillations are in phase or out of phase, respectively. For two delays we find an intricate dynamics which comprises an irregular alteration of small amplitude oscillations, in-phase and out-of-phase large amplitude oscillations, also called extreme events. This transient chaotic dynamics is sandwiched between a bubbling transition and a blowout bifurcation.

Many neurons in the auditory system of the brain must encode periodic signals. These neurons under periodic stimulation display rich dynamical states including mode locking and chaotic responses. Periodic stimuli such as sinusoidal waves and amplitude modulated sounds can lead to various forms of $n:m$ mode-locked states, in which a neuron fires $n$ action potentials per $m$ cycles of the stimulus. Here, we study mode-locking in the Izhikevich neurons, a reduced model of the Hodgkin–Huxley neurons. The Izhikevich model is much simpler in terms of the dimension of the coupled nonlinear differential equations compared with other existing models, but excellent for generating the complex spiking patterns observed in real neurons. We obtained the regions of existence of the various mode-locked states on the frequency-amplitude plane, called Arnold tongues, for the Izhikevich neurons. Arnold tongue analysis provides useful insight into the organization of mode-locking behavior of neurons under periodic forcing. We find these tongues for both class-1 and class-2 excitable neurons in both deterministic and noisy regimes.

The present interdisciplinary study combines electromagnetics and combustion to unveil an original and basic experiment displaying a spontaneous flame instability that is mitigated as the non-premixed sooting flame experiences a magnetic perturbation. This magnetic instability mitigation is reproduced by direct numerical simulations to be further elucidated by a flow stability analysis. A key role in the stabilization process is attributed to the momentum and thermochemistry coupling that the magnetic force, acting mainly on paramagnetic oxygen, contributes to sustain. The spatial local stability analysis based on the numerical simulations shows that the magnetic field tends to reduce the growth rates of small flame perturbations.

A data-science approach to solving the ill-conditioned inverse problem for analytical continuation is proposed. The root of the problem lies in the fact that even tiny noise of imaginary-time input data has a serious impact on the inferred real-frequency spectra. By means of a modern regularization technique, we eliminate redundant degrees of freedom that essentially carry the noise, leaving only relevant information unaffected by the noise. The resultant spectrum is represented with minimal bases and thus a stable analytical continuation is achieved. This framework further provides a tool for analyzing to what extent the Monte Carlo data need to be accurate to resolve details of an expected spectral function.

We derive hydrodynamics of paradigmatic conserved-mass transport processes on a ring. The systems, governed by chipping, diffusion, and coalescence of masses, eventually reach a nonequilibrium steady state, having nontrivial correlations, with steady-state measures in most cases not known. In these processes, we analytically calculate two transport coefficients, bulk-diffusion coefficient and conductivity. Remarkably, the two transport coefficients obey an equilibrium-like Einstein relation even when the microscopic dynamics violates detailed balance and systems are far from equilibrium. Moreover, we show, using a macroscopic fluctuation theory, that the probability of large deviation in density, obtained from the above hydrodynamics, is in complete agreement with the same derived earlier by Das et al. [Phys. Rev. E 93, 062135 (2016)] using an additivity property.

In any general cycle of measurement, feedback, and erasure, the measurement will reduce the entropy of the system when information about the state is obtained, while erasure, according to Landauer's principle, is accompanied by a corresponding increase in entropy due to the compression of logical and physical phase space. The total process can in principle be fully reversible. A measurement error reduces the information obtained and the entropy decrease in the system. The erasure still gives the same increase in entropy, and the total process is irreversible. Another consequence of measurement error is that a bad feedback is applied, which further increases the entropy production if the proper protocol adapted to the expected error rate is not applied. We consider the effect of measurement error on a realistic single-electron box Szilard engine, and we find the optimal protocol for the cycle as a function of the desired power $P$ and error $\varepsilon$.

Infrared and dielectric spectroscopic techniques are used to investigate the characteristics of two chiral smectics, namely, 1,1,3,3,5,5,5-heptamethyltrisiloxane 1-[$4^{\prime }$-(undecyl-1-oxy)-4-biphenyl(S,S)-2-chloro-3-methylpentanoate] ($\mathrm{MS}{\mathrm{i}}_3\mathrm{M}{\mathrm{R}}_{11}$) and tricarbosilane-hexyloxy-benzoic acid (S)-4′-(1-methyl-hexyloxy)-3′-nitro-biphenyl-4-yl ester (W599). The orientational features and the field dependencies of the apparent tilt angle and the dichroic ratio for homogeneous planar-aligned samples were calculated from the absorbance profiles obtained at different temperatures especially in the smectic-A* phase of these liquid crystals. The dichroic ratios of the C-C phenyl ring stretching vibrations were considered for the determination of the tilt angle at different temperatures and different voltages. The low values of the order parameter obtained with and without an electric field applied across the cell in the $\mathrm{Sm}-A^{*}$ phase for both smectics are consistent with the de Vries concept. The generalized Langevin-Debye model introduced in the literature for explaining the electro-optical response has been applied to the results from infrared spectroscopy. The results show that the dipole moment of the tilt-correlated domain diverges as the transition temperature from $\mathrm{Sm}-A^{*}$ to $\mathrm{Sm}-C^{*}$ is approached. The Debye-Langevin model is found to be extremely effective in confirming some of the conclusions of the de Vries chiral smectics and gives additional results on the order parameter and the dichroic ratio as a function of the field across the cell. Dielectric spectroscopy finds large dipolar fluctuations in the $\mathrm{Sm}-A^{*}$ phase for both compounds and again these confirm their de Vries behavior.