We use the Navier-Stokes-Fourier constitutive equations to study plane shock waves in dilute gases. It is shown that the experimental information on the normalized density profiles can be fit by using the so-called soft sphere model, in which the viscosity and thermal conductivity are proportional to a power of the temperature.

The ion-ion dynamical structure factor and the equation of state of warm dense aluminum in a two-temperature quasiequilibrium state, with the electron temperature higher than the ion temperature, are investigated using molecular-dynamics simulations based on ion-ion pair potentials constructed from a neutral pseudoatom model. Such pair potentials based on density functional theory are parameter-free and depend directly on the electron temperature and indirectly on the ion temperature, enabling efficient computation of two-temperature properties. Comparison with ab initio simulations and with other average-atom calculations for equilibrium aluminum shows good agreement, justifying a study of quasiequilibrium situations. Analyzing the van Hove function, we find that ion-ion correlations vanish in a time significantly smaller than the electron-ion relaxation time so that dynamical properties have a physical meaning for the quasiequilibrium state. A significant increase in the speed of sound is predicted from the modification of the dispersion relation of the ion acoustic mode as the electron temperature is increased. The two-temperature equation of state including the free energy, internal energy, and pressure is also presented.

In this paper, a general propagation lattice Boltzmann model is proposed for nonlinear advection-diffusion equations (NADEs), and the Chapman-Enskog analysis shows that the NADEs with variable coefficients can be recovered correctly from the present model. We also perform some simulations of the linear advection-diffusion equation, nonlinear heat conduction equation, NADEs with anisotropic diffusion, and variable coefficients to test the present model, and find that the numerical results agree well with the corresponding analytical solutions. Moreover, it is also shown that by properly adjusting the two free parameters introduced into the propagation step, the present model could be more stable and more accurate than the standard lattice Bhatnagar-Gross-Krook model.

Random sequential adsorption (RSA) is a time-dependent packing process, in which particles of certain shapes are randomly and sequentially placed into an empty space without overlap. In the infinite-time limit, the density approaches a “saturation” limit. Although this limit has attracted particular research interest, the majority of past studies could only probe this limit by extrapolation. We have previously found an algorithm to reach this limit using finite computational time for spherical particles and could thus determine the saturation density of spheres with high accuracy. In this paper, we generalize this algorithm to generate saturated RSA packings of two-dimensional polygons. We also calculate the saturation density for regular polygons of three to ten sides and obtain results that are consistent with previous, extrapolation-based studies.

We investigate dilution effects on the classical spin-ice materials such as ${\mathrm{Ho}}_2{\mathrm{Ti}}_2{\mathrm{O}}_7$ and ${\mathrm{Dy}}_2{\mathrm{Ti}}_2{\mathrm{O}}_7$. In particular, we derive a formula of the thermodynamic quantities as functions of the temperature and a nonmagnetic ion concentration based on a Husimi-cactus approximation. We find that the formula predicts a dilution-induced crossover from the cooperative to the conventional paramagnets in a ground state, and that it also reproduces the “generalized Pauling's entropy” given by Ke et al. To verify the formula from a numerical viewpoint, we compare these results with Monte Carlo simulation calculation data, and then find good agreement for all parameter values.

We study avalanches in the Kuramoto model, defined as excursions of the order parameter due to ephemeral episodes of synchronization. We present scaling collapses of the avalanche sizes, durations, heights, and temporal profiles, extracting scaling exponents, exponent relations, and scaling functions that are shown to be consistent with the scaling behavior of the power spectrum, a quantity independent of our particular definition of an avalanche. A comprehensive scaling picture of the noise in the subcritical finite-$N$ Kuramoto model is developed, linking this undriven system to a larger class of driven avalanching systems.

We investigate the dynamics of highly polydisperse finite granular chains. From the spatiospectral properties of small vibrations, we identify which particular single-particle displacements lead to energy localization. Then, we address a fundamental question: Do granular nonlinearities and the resulting chaotic dynamics destroy this energy localization? Our numerical simulations show that for moderate nonlinearities, the overall system behaves chaotically, and spreading of energy occurs. However, long-lasting chaotic energy localization is observed for particular single-particle excitations in the presence of the nonsmooth nonlinearities. On the other hand, for sufficiently strong nonlinearities, the granular chain reaches energy equipartition. In this case, an equilibrium chaotic state is reached independent of the initial position excitation.

We generalize a binary majority-vote model on adaptive networks to its plurality-vote counterpart and analyze the timescale to consensus when voters are given more than two options. When opinions are uniformly distributed in the population of voters in the initial state, we find that the timescale to consensus is shorter than the binary vote model from both numerical simulations and mathematical analysis using the master equation for the three-state plurality-vote model. When intervention such as opinion conversion is allowed, as in the case of sudden change of mind of voter for any reason, the effort needed to push the fragmented three-opinion population in the thermodynamic limit to the consensus state, measured in minimal intervention cost, is less than that needed to push a polarized two-opinion population to the consensus state, when the degree ($p$) of homophily is less than 0.8. For a finite system, the fragmented three-opinion population will spontaneously reach the consensus state, with faster time to consensus, compared to polarized two-opinion population, for a broad range of $p$.

The biological processes of cellular decision making and differentiation involve a plethora of signaling pathways and gene regulatory circuits. These networks in turn exhibit a multitude of motifs playing crucial parts in regulating network activity. Here we compare the topological placement of motifs in gene regulatory and signaling networks and observe that it suggests different evolutionary strategies in motif distribution for distinct cellular subnetworks.

A molecular field theory and coarse-grained computer simulations with dissipative particle dynamics have been used to study the spontaneous orientational ordering of anisotropic nanoparticles in the lamellar and hexagonal phases of diblock copolymers and the effect of nanoparticles on the phase behavior of these systems. Both the molecular theory and computer simulations indicate that strongly anisotropic nanoparticles are ordered orientationally mainly in the boundary region between the domains and the nematic order parameter possesses opposite signs in adjacent domains. The orientational order is induced by the boundary and by the interaction between nanoparticles and the monomer units in different domains. In simulations, sufficiently long and strongly selective nanoparticles are ordered also inside the domains. The nematic order parameter and local concentration profiles of nanoparticles have been calculated numerically using the model of a nanoparticle with two interaction centers and also determined using the results of computer simulations. A number of phase diagrams have been obtained which illustrate the effect of nanoparticle selectivity and molar fraction of the stability ranges of various phases. Different morphologies have been identified by analyzing the static structure factor and a phase diagram has been constructed in coordinates' nanoparticle concentration–copolymer composition. Orientational ordering of even a small fraction of nanoparticles may result in a significant increase of the dielectric anisotropy of a polymer nanocomposite, which is important for various applications.

We consider a lattice gas in spaces of dimensionality $\mathcal{D}=1,2,3$. The particles are subject to a hardcore exclusion interaction and an attractive pair interaction that satisfies Gauss' law as do Newtonian gravity in $\mathcal{D}=3$, a logarithmic potential in $\mathcal{D}=2$, and a distance-independent force in $\mathcal{D}=1$. Under mild additional assumptions regarding symmetry and fluctuations we investigate equilibrium states of self-gravitating material clusters, in particular radial density profiles for closed and open systems. We present exact analytic results in several instances and high-precision numerical data in others. The density profile of a cluster with finite mass is found to exhibit exponential decay in $\mathcal{D}=1$ and power-law decay in $\mathcal{D}=2$ with temperature-dependent exponents in both cases. In $\mathcal{D}=2$ the gas evaporates in a continuous transition at a nonzero critical temperature. We describe clusters of infinite mass in $\mathcal{D}=3$ with a density profile consisting of three layers (core, shell, halo) and an algebraic large-distance asymptotic decay. In $\mathcal{D}=3$ a cluster of finite mass can be stabilized at $T>0$ via confinement to a sphere of finite radius. In some parameter regime, the gas thus enclosed undergoes a discontinuous transition between distinct density profiles. For the free energy needed to identify the equilibrium state we introduce a construction of gravitational self-energy that works in all $\mathcal{D}$ for the lattice gas. The decay rate of the density profile of an open cluster is shown to transform via a stretched exponential for $1<\mathcal{D}<2$, whereas it crosses over from one power-law at intermediate distances to a different power-law at larger distances for $2<\mathcal{D}<3$.

External forces acting on a microswimmer can feed back on its self-propulsion mechanism. We discuss this load response for a generic microswimmer that swims by cyclic shape changes. We show that the change in cycle frequency is proportional to the Lighthill efficiency of self-propulsion. As a specific example, we consider Najafi's three-sphere swimmer. The force-velocity relation of a microswimmer implies a correction for a formal superposition principle for active and passive motion.

We investigate the effects of flow on the dynamics of ferromagnetic nematic liquid crystals. As a model, we study the coupled dynamics of the magnetization, $\mathbf{M}$, the director field, $\mathbf{n}$, associated with the liquid crystalline orientational order, and the velocity field, $\mathbf{v}$. We evaluate how simple shear flow in a ferromagnetic nematic is modified in the presence of small external magnetic fields, and we make experimentally testable predictions for the resulting effective shear viscosity: an increase by a factor of 2 in a magnetic field of about 20 mT. Flow alignment, a characteristic feature of classical uniaxial nematic liquid crystals, is analyzed for ferromagnetic nematics for the two cases of magnetization in or perpendicular to the shear plane. In the former case, we find that small in-plane magnetic fields are sufficient to suppress tumbling and thus that the boundary between flow alignment and tumbling can be controlled easily. In the latter case, we furthermore find a possibility of flow alignment in a regime for which one obtains tumbling for the pure nematic component. We derive the analogs of the three Miesowicz viscosities well-known from usual nematic liquid crystals, corresponding to nine different configurations. Combinations of these can be used to determine several dynamic coefficients experimentally.

We provide a criterion based on graph duality to predict whether the three-state Potts antiferromagnet on a plane quadrangulation has a zero- or finite-temperature critical point, and its universality class. The former case occurs for quadrangulations of self-dual type, and the zero-temperature critical point has central charge $c=1$. The latter case occurs for quadrangulations of non-self-dual type, and the critical point belongs to the universality class of the three-state Potts ferromagnet. We have tested this criterion against high-precision computations on four lattices of each type, with very good agreement. We have also found that the Wang-Swendsen-Kotecký algorithm has no critical slowing-down in the former case, and critical slowing-down in the latter.

Knotted fields enrich a variety of physical phenomena, ranging from fluid flows, electromagnetic fields, to textures of ordered media. Maxwell's electrostatic equations, whose vacuum solution is mathematically known as a harmonic field, provide an ideal setting to explore the role of domain topology in determining physical fields in confined space. In this work, we show the uniqueness of a harmonic field in knotted tubes, and reduce the construction of a harmonic field to a Neumann boundary value problem. By analyzing the harmonic field in typical knotted tubes, we identify the torsion driven transition from bipolar to vortex patterns. We also analogously extend our discussion to the organization of liquid crystal textures in knotted tubes. These results further our understanding about the general role of topology in shaping a physical field in confined space, and may find applications in the control of physical fields by manipulation of surface topology.

We derive open quantum maps from periodically kicked scattering systems and discuss the computation of their resonance spectra in terms of theoretically grounded methods, such as complex scaling and sufficiently weak absorbing potentials. In contrast, we also show that current implementations of open quantum maps, based on strong absorptive or even projective openings, fail to produce the resonance spectra of kicked scattering systems. This comparison pinpoints flaws in current implementations of open quantum maps, namely, the inability to separate resonance eigenvalues from the continuum as well as the presence of diffraction effects due to strong absorption. The reported deviations from the true resonance spectra appear, even if the openings do not affect the classical trapped set, and become appreciable for shorter-lived resonances, e.g., those associated with chaotic orbits. This makes the open quantum maps, which we derive in this paper, a valuable alternative for future explorations of quantum-chaotic scattering systems, for example, in the context of the fractal Weyl law. The results are illustrated for a quantum map model whose classical dynamics exhibits key features of ionization and a trapped set which is organized by a topological horseshoe.

We study cluster synchronization in networks and show that the stability of all possible cluster synchronization patterns depends on a small set of Lyapunov exponents. Our approach can be applied to clusters corresponding to both orbital partitions of the network nodes (symmetry-cluster synchronization) and equitable partitions of the network nodes (input-cluster synchronization). Our results are verified experimentally in networks of coupled optoelectronic oscillators.