The role of fixed degrees of freedom in soft or granular matter systems has broad applicability and theoretical interest. Here we address questions of the geometrical role that a scaffolding of fixed particles plays in tuning the threshold volume fraction and force network in the vicinity of jamming. Our two-dimensional simulated system consists of soft particles and fixed “pins,” both of which harmonically repel overlaps. On the one hand, we find that many of the critical scalings associated with jamming in the absence of pins continue to hold in the presence of even dense pin latices. On the other hand, the presence of pins lowers the jamming threshold in a universal way at low pin densities and a geometry-dependent manner at high pin densities, producing packings with lower densities and fewer contacts between particles. The onset of strong lattice dependence coincides with the development of bond-orientational order. Furthermore, the presence of pins dramatically modifies the network of forces, with both unusually weak and unusually strong forces becoming more abundant. The spatial organization of this force network depends on pin geometry and is described in detail. Using persistent homology, we demonstrate that pins modify the topology of the network. Finally, we observe clear signatures of this developing bond-orientational order and broad force distribution in the elastic moduli which characterize the linear response of these packings to strain.

In this paper, a high-order solver combining the flux reconstruction (FR) method and the thermal lattice Boltzmann flux solver (FRTLBFS) is developed for accurately and efficiently simulating incompressible thermal flows. The conservative differential equations recovered from Chapman-Enskog analysis of the thermal lattice Boltzmann equation are solved by the high-order FR method. The thermal lattice Boltzmann method is only applied to reconstruct the local solution used for evaluating fluxes at the solution and flux points. Unlike the traditional Navier-Stokes–Boussinesq (NSB) solvers where the inviscid and viscous terms are treated separately, the inviscid and viscous fluxes in the current FRTLBFS are coupled and computed uniformly. In comparison with the recently developed high-order flux reconstruction thermal lattice Boltzmann method, the FRTLBFS holds advantages such as high-order accuracy, good stability, and compactness but is more efficient and low storage, since only macroscopic flow variables including density, velocity, and temperature are stored and evolved. In addition, the physical boundary conditions in FRTLBFS can be directly implemented by using the same method as in conventional NSB solvers. Numerical validations of the proposed method are implemented by simulating (a) the porous plate problem, (b) natural convection in a square cavity, (c) unsteady natural convection in a tall cavity, and (d) thermal lid-driven cavity flow. Numerical results demonstrate that the present solver is an attractive tool to simulate incompressible thermal flows due to its high-order accuracy, stability, and low memory cost.

Coupled rotors can spontaneously synchronize, giving rise to a plethora of intriguing dynamics. We present here a pair of spiral waves as two synchronizing rotors, coupled by diffusion. The spirals are pinned to unexcitable obstacles, which enables us to modify their frequencies and restrain their drift. In experiments with the Belousov-Zhabotinsky reaction, we show that two counterrotating spiral rotors, pinned to circular heterogeneities, can synchronize in frequency and phase. The nature of the phase synchronization varies depending on the difference in their characteristic frequencies. We observe in-phase and out-of-phase synchronization, lag synchronization, and phase resetting across the experiments. The time required for the two spirals to synchronize is found to depend upon the relative size of their pinning obstacles and the distance separating them. This distance can also modify the phase lag of the two rotors upon synchronization. Our experimental observations are reproduced and explained further on the basis of numerical simulations of an excitable reaction-diffusion model.

We study influence maximization on temporal networks. This is a special setting where the influence function is not submodular, and there is no optimality guarantee for solutions achieved via greedy optimization. We perform an exhaustive analysis on both real and synthetic networks. We show that the influence function of randomly sampled sets of seeds often violates the necessary conditions for submodularity. However, when sets of seeds are selected according to the greedy optimization strategy, the influence function behaves effectively as a submodular function. Specifically, violations of the necessary conditions for submodularity are never observed in real networks, and only rarely in synthetic ones. The direct comparison with exact solutions obtained via brute-force search indicates that the greedy strategy provides approximate solutions that are well within the optimality gap guaranteed for strictly submodular functions. Greedy optimization appears, therefore, to be an effective strategy for the maximization of influence on temporal networks.

Due to significant computational expense, discrete element method simulations of jammed packings of size-dispersed spheres with size ratios greater than 1:10 have remained elusive, limiting the correspondence between simulations and real-world granular materials with large size dispersity. Invoking a recently developed neighbor binning algorithm, we generate mechanically stable jammed packings of frictionless spheres with power-law size distributions containing up to nearly 4 000 000 particles with size ratios up to 1:100. By systematically varying the width and exponent of the underlying power laws, we analyze the role of particle size distributions on the structure of jammed packings. The densest packings are obtained for size distributions that balance the relative abundance of large-large and small-small particle contacts. Although the proportion of rattler particles and mean coordination number strongly depend on the size distribution, the mean coordination of nonrattler particles attains the frictionless isostatic value of six in all cases. The size distribution of nonrattler particles that participate in the load-bearing network exhibits no dependence on the width of the total particle size distribution beyond a critical particle size for low-magnitude exponent power laws. This signifies that only particles with sizes greater than the critical particle size contribute to the mechanical stability. However, for high-magnitude exponent power laws, all particle sizes participate in the mechanical stability of the packing.

The Maxwell equations-based 3D-analytical solution for the terahertz (THz) half-cycle electromagnetic wave transition radiation pulse has been found. This solution describes generation and propagation of transition radiation into free space from laser-produced relativistic electron bunch which crosses a target-vacuum interface as a result of ultrashort laser pulse interaction with a thin high-conductivity target. The analytical solution found complements the theory of laser initiated transition radiation. It describes the THz wave half-cycle pulse at the arbitrary distance from a target surface including near-field zone rather than its standard far-field characterization. The analytical research has also been supplemented with the 3D simulations using the finite-difference time-domain method, which makes it possible for description of much wider spatial domain as compared to that from the particle-in-cell approach. The presented result sheds light fundamentally on the interference of the electron bunch field and the generated THz field of broadband transition radiation from laser-plasma interaction. The latter is studied for a long time in the experiments with solid density plasma and the theory developed may inspire to targeted measurements and investigations of unique super intense half-cycle THz radiation waves near the laser target.

The friction of an incident particle interacting with the background molecules is a cornerstone in the nonequilibrium dynamics and statistics. It is reported that the Stokes force may fail while the Brown particle's size is small enough. In this work, the mean collision force of a small classical particle moving through the rarefied gases is analyzed by the direct calculation of the mean decrease of particle's velocity by elastic collisions. As an example, a whole velocity space applicable mean collision force in Maxwell gas is obtained. A self-consistent solution is further provided based on the numerical simulations. Within the low speed limit, comparison of the friction and the Stokes force has been demonstrated. Although the two forces are both proportional to the speed of the particle, their coefficients are different. Unlike the linear speed dependence of Stokes force, the linear behavior in rarefied gases is broken with increasing the speed of incident particle, and a quadratic speed dependence is resulted in high speed. This work clarifies the nonequilibrium dynamics of microscopic particles moving in rarefied gases, and can improve our microscopic understanding of the collision force.

We simulate the transport of colloidal particles driven by a static and homogeneous drift force, and subject to the optical potential created by two lattices of optical tweezers. The lattices of optical tweezers are parallel to each other, shifted, and rotated by a twist angle. Due to a negative interference between the potential of the two lattices, flat channels appear in the total optical potential. At specific twist angles, known as magic angles, the flat channels percolate the entire system and the colloidal particles can then be transported using a weak external drift force. We characterize the transport in both square and hexagonal lattices of twisted optical tweezers.

We show that the Schnakenberg's entropy production rate in a master equation is lower bounded by a function of the weight of the Markov graph, here defined as the sum of the absolute values of probability currents over the edges. The result is valid for time-dependent nonequilibrium entropy production rates. Moreover, in a general framework, we prove a theorem showing that the Kullback-Leibler divergence between distributions $P\left(s\right)$ and $P^{\prime }\left(s\right):=P(m\left(s\right))$, where $m$ is an involution, $m(m\left(s\right))=s$, is lower bounded by a function of the total variation of $P$ and $P^{\prime }$, for any $m$. The bound is tight and it improves on Pinsker's inequality for this setup. This result illustrates a connection between nonequilibrium thermodynamics and graph theory with interesting applications.

We analyze the effect of small-amplitude noise and heterogeneity in a network of coupled excitable oscillators with strong timescale separation. Using mean-field analysis, we uncover the mechanism of a nontrivial effect—diversity-induced decoherence (DIDC)—in which heterogeneity modulates the mechanism of self-induced stochastic resonance to inhibit the coherence of oscillations. We argue that DIDC may offer one possible mechanism via which, in excitable neural systems, generic heterogeneity and background noise can synergistically prevent unwanted resonances that may be related to hyperkinetic movement disorders.

Classically, the quantity of contact area $A_R$ between two bodies is considered a proxy for the force of friction. However, bond density across the interface—quality of contact—is also relevant, and contemporary debate often centers around the relative importance of these two factors. In this work, we demonstrate that a third factor, often overlooked, plays a significant role in static frictional strength: The spatial distribution of contact. We perform static friction measurements, $\mu$, on three pairs of solid blocks while imaging the contact plane. By using linear regression on hundreds of image-$\mu$ pairs, we are able to predict future friction measurements with three to seven times better accuracy than existing benchmarks, including total quantity of contact area. Our model has no access to quality of contact, and we therefore conclude that a large portion of the interfacial state is encoded in the spatial distribution of contact, rather than its quality or quantity.

We investigate quantum Otto engine and refrigeration cycles of a time-dependent harmonic oscillator operating under the conditions of maximum $\mathrm{\Omega }$ function, a trade-off objective function which represents a compromise between energy benefits and losses for a specific job, for both adiabatic and nonadiabatic (sudden) frequency modulations. We derive analytical expressions for the efficiency and coefficient of performance of the Otto cycle. For the case of adiabatic driving, we point out that in the low-temperature regime, the harmonic Otto engine (refrigerator) can be mapped to Feynman's ratchet and pawl model which is a steady-state classical heat engine. For the sudden switch of frequencies, we obtain loop-like behavior of the efficiency-work curve, which is characteristic of irreversible heat engines. Finally, we discuss the behavior of cooling power at maximum $\mathrm{\Omega }$ function.

Many natural and man-made systems require suitable feedback to function properly. In this study, we aim to investigate the impact of additional complex conjugate feedback on globally coupled Stuart-Landau oscillators. We find that this additional feedback results in the onset of symmetry breaking clusters and out-of-phase clusters. Interestingly, we also find the existence of explosive amplitude death along with disparate multistable states. We characterize the first-order transition to explosive death through the amplitude order parameter and show that the transition from oscillatory to death state indeed shows a hysteresis nature. Further, we map the global dynamical transitions in the parametric spaces. In addition, to understand the existence of multistabilities and their transitions, we analyze the bifurcation scenarios of the reduced model and also explore their basin stability. Our study will shed light on the emergent dynamics in the presence of additional feedback.

We study the effect of Kerr anharmonicity on the symmetry-breaking phenomena of coupled quantum oscillators. Two types of symmetry-breaking processes are studied, namely, the inhomogeneous steady state (or quantum oscillation death state) and quantum chimera state. Remarkably, it is found that Kerr nonlinearity hinders the process of symmetry breaking in both the cases. We establish our results using direct simulation of the quantum master equation and analysis of the stochastic semiclassical model. Interestingly, in the case of quantum oscillation death, an increase in the strength of Kerr nonlinearity tends to favor the symmetry and at the same time decreases the degree of quantum mechanical entanglement. This paper presents a useful mean to control and engineer symmetry-breaking states for quantum technology.

We address various regimes of charge transfer in DNA within the framework of the Peyrard-Bishop-Holstein model and analyze them from the standpoint of the characteristic size and timescales of the electronic and vibrational subsystems. It is demonstrated that a polaron is an unstable configuration within a broad range of temperatures and therefore polaronic contribution to the charge transport is irrelevant. We put forward an alternative fluctuation-governed charge transfer mechanism and show that the charge transfer can be quasiballistic at low temperatures, diffusive or mixed at intermediate temperatures, and subdiffusive close to the DNA denaturation transition point. Dynamic fluctuations in the vibrational subsystem is the key ingredient of our proposed mechanism which allows for explanation of all charge transfer regimes at finite temperatures. In particular, we demonstrate that in the most relevant regime of high temperatures (above the aqueous environment freezing point), the electron dynamics is completely governed by relatively slow fluctuations of the mechanical subsystem. We argue also that our proposed analysis methods and mechanisms can be relevant for the charge transfer in other organic systems, such as conjugated polymers, molecular aggregates, $\alpha$-helices, etc.

It is shown that the normalized rms fluctuation of the shear modulus on the medium-range order scale in glasses correlates with fragility: the higher fragility, the smaller the fluctuation amplitude. The latter is calculated within the heterogeneous elasticity theory using the data on the boson peak in glasses. On a smaller scale corresponding to cooperative structural relaxation, the normalized rms fluctuation of the infinite-frequency shear modulus was estimated using the data on the decoupling of viscosity and diffusion in supercooled liquids. These fluctuations are much smaller in amplitude, and, in contrast, they increase with increasing fragility. Extrapolation predicts intersection of both rms fluctuations and disappearing of the boson peak at the upper limit to fragility ≈180.

In this work, we consider a general consistent and conservative phase-field model for the incompressible two-phase flows. In this model, not only the Cahn-Hilliard or Allen-Cahn equation can be adopted, but also the mass and the momentum fluxes in the Navier-Stokes equations are reformulated such that the consistency of reduction, consistency of mass and momentum transport, and the consistency of mass conservation are satisfied. We further develop a lattice Boltzmann (LB) method, and show that through the direct Taylor expansion, the present LB method can correctly recover the consistent and conservative phase-field model. Additionally, if the divergence of the extra momentum flux is seen as a force term, the extra force in the present LB method would include another term which has not been considered in the previous LB methods. To quantitatively evaluate the incompressibility and the consistency of the mass conservation, two statistical variables are introduced in the study of the deformation of a square droplet, and the results show that the present LB method is more accurate. The layered Poiseuille flow and a droplet spreading on an ideal wall are further investigated, and the numerical results are in good agreement with the analytical solutions. Finally, the problems of the Rayleigh-Taylor instability, a single rising bubble, and the dam break with the high Reynolds numbers and/or large density ratios are studied, and it is found that the present consistent and conservative LB method is robust for such complex two-phase flows.

In a recent series of papers, an exact combinatorial solution was claimed for a variant of the so-called Marcus-Lushnikov model of aggregation. In this model, a finite number of aggregates are initially assumed to be present in the form of monomers. At each time step, two aggregates are chosen according to certain size-dependent probabilities and irreversibly joined to form an aggregate of higher mass. The claimed result given an expression for the full probability distribution over all possible size distributions in terms of the so-called Bell polynomials. In this paper, we develop the asymptotics of this solution in order to check whether the exact solution yields correct expressions for the average cluster size distribution as obtained from the Smoluchowski equations. The answer is surprisingly involved: For the generic case of an arbitrary reaction rate, it is negative, but for the so-called classical rate kernels, constant, additive, and multiplicative, the solutions obtained are indeed exact. On the other hand, for the multiplicative kernel, a discrepancy is found in the full solution between the combinatorial solution and the exact solution. The reasons for this puzzling pattern of agreement and disagreement are unclear. A better understanding of the combinatorial solution's derivation is needed, the better to understand its range of validity.

We revisit quantum Maxwell demon in thermodynamic feedback cycle in the steady-state regime. We derive a generalized version of the Clausius inequality for a finite-time steady feedback cycle with a single heat bath. It is shown to be tighter than previously known ones, and allows us to clarify that feedback control is necessary to violate the standard Clausius inequality.