We investigate the Taylor-Couette system where the radius ratio is close to unity. Systematically increasing the Reynolds number, we observe a number of previously known transitions, such as one from the classical Taylor vortex flow (TVF) to wavy vortex flow (WVF) and the transition to fully developed turbulence. Prior to the onset of turbulence, we observe intermittent bursting patterns of localized turbulent patches, confirming the experimentally observed pattern of very short wavelength bursts (VSWBs). A striking finding is that, for a Reynolds number larger than that for the onset of VSWBs, a new type of intermittently bursting behavior emerges: patterns of azimuthally closed rings of various orders. We call them ring-bursting patterns, which surround the cylinder completely but remain localized and separated in the axial direction through nonturbulent wavy structures. We employ a number of quantitative measures including the cross-flow energy to characterize the ring-bursting patterns and to distinguish them from the background flow. These patterns are interesting because they do not occur in the wide-gap Taylor-Couette flow systems. The narrow-gap regime is less studied but certainly deserves further attention to gain deeper insights into complex flow dynamics in fluids.

The energy loss of light ions in dense plasmas is investigated with special focus on low to medium projectile energies, i.e., at velocities where the maximum of the stopping power occurs. In this region, exceptionally large theoretical uncertainties remain and no conclusive experimental data are available. We perform simulations of beam-plasma configurations well suited for an experimental test of ion energy loss in highly ionized, laser-generated carbon plasmas. The plasma parameters are extracted from two-dimensional hydrodynamic simulations, and a Monte Carlo calculation of the charge-state distribution of the projectile ion beam determines the dynamics of the ion charge state over the whole plasma profile. We show that the discrepancies in the energy loss predicted by different theoretical models are as high as 20–30%, making these theories well distinguishable in suitable experiments.

We present molecular dynamics (MD) calculations of shear viscosity for asymmetric mixed plasma for thermodynamic conditions relevant to astrophysical and inertial confinement fusion plasmas. Specifically, we consider mixtures of deuterium and argon at temperatures of 100–500 eV and a number density of ${10}^{25}$ ions/cc. The motion of 30 000–120 000 ions is simulated in which the ions interact via the Yukawa (screened Coulomb) potential. The electric field of the electrons is included in this effective interaction; the electrons are not simulated explicitly. Shear viscosity is calculated using the Green-Kubo approach with an integral of the shear stress autocorrelation function, a quantity calculated in the equilibrium MD simulations. We systematically study different mixtures through a series of simulations with increasing fraction of the minority high-$Z$ element (Ar) in the D-Ar plasma mixture. In the more weakly coupled plasmas, at 500 eV and low Ar fractions, results from MD compare very well with Chapman-Enskog kinetic results. In the more strongly coupled plasmas, the kinetic theory does not agree well with the MD results. We develop a simple model that interpolates between classical kinetic theories at weak coupling and the Murillo Yukawa viscosity model at higher coupling. This hybrid kinetics-MD viscosity model agrees well with the MD results over the conditions simulated, ranging from moderately weakly coupled to moderately strongly coupled asymmetric plasma mixtures.

We have studied Markov processes on denumerable state space and continuous time. We found that all these processes are connected via gauge transformations. We have used this result before as a method to resolve equations, included the case in a previous work in which the sample space is time-dependent [Phys. Rev. E 90, 022125 (2014)]. We found a general solution through dilation of the state space, although the prior probability distribution of the states defined in this new space takes smaller values with respect to that in the initial problem. The gauge (local) group of dilations modifies the distribution on the dilated space to restore the original process. In this work, we show how the Markov process in general could be linked via gauge (local) transformations, and we present some illustrative examples for this result.

We consider local heating of a part of a two-dimensional bilobal enclosure of a varying cross section confining a system of overdamped Brownian particles. Since varying cross section in higher dimension results in an entropic potential in lower dimension, local heating alters the relative stability of the entropic states. We show that this blow-torch effect modifies the entropic potential in a significant way so that the resultant effective entropic potential carries both the features of variation of width of the confinement and variation of temperature along the direction of transport. The reduced probability distribution along the direction of transport calculated by full numerical simulations in two dimensions agrees well with our analytical findings. The extent of population transfer in the steady state quantified in terms of the integrated probability of residence of the particles in either of the two lobes exhibits interesting variation with the mean position of the heated region. Our study reveals that heating around two particular zones of a given lobe maximizes population transfer to the other.

We show that the complex connectivity of percolating networks of nanoparticles provides a natural solid-state system in which bottom-up assembly provides a route to realization of neuromorphic behavior. Below the percolation threshold the networks comprise groups of particles separated by tunnel gaps; an applied voltage causes atomic scale wires to form in the gaps, and we show that the avalanche of switching events that occurs is similar to potentiation in biological neural systems. We characterize the level of potentiation in the percolating system as a function of the surface coverage of nanoparticles and other experimentally relevant variables, and compare our results with those from biological systems. The complex percolating structure and the electric field driven switching mechanism provide several potential advantages in comparison to previously reported solid-state neuromorphic systems.

A granular front emerges whenever the free-surface flow of a concentrated suspension spontaneously alters its internal structure, exhibiting a higher concentration of particles close to its front. This is a common and yet unexplained phenomenon, which is usually believed to be the result of fluid convection in combination with particle size segregation. However, suspensions composed of uniformly sized particles also develop a granular front. Within a large rotating drum, a stationary recirculating avalanche is generated. The flowing material is a mixture of a viscoplastic fluid obtained from a kaolin-water dispersion with spherical ceramic particles denser than the fluid. The goal is to mimic the composition of many common granular-fluid materials, such as fresh concrete or debris flow. In these materials, granular and fluid phases have the natural tendency to separate due to particle settling. However, through the shearing caused by the rotation of the drum, a reorganization of the phases is induced, leading to the formation of a granular front. By tuning the particle concentration and the drum velocity, it is possible to control this phenomenon. The setting is reproduced in a numerical environment, where the fluid is solved by a lattice-Boltzmann method, and the particles are explicitly represented using the discrete element method. The simulations confirm the findings of the experiments, and provide insight into the internal mechanisms. Comparing the time scale of particle settling with the one of particle recirculation, a nondimensional number is defined, and is found to be effective in predicting the formation of a granular front.

We consider a clustered network with small-world subnetworks of inhibitory fast spiking interneurons and investigate the effect of intermodular connection on the emergence of fast sparsely synchronized rhythms by varying both the intermodular coupling strength $J_{\mathrm{inter}}$ and the average number of intermodular links per interneuron $M_{\mathrm{syn}}^{\left(\mathrm{inter}\right)}$. In contrast to the case of nonclustered networks, two kinds of sparsely synchronized states such as modular and global synchronization are found. For the case of modular sparse synchronization, the population behavior reveals the modular structure, because the intramodular dynamics of subnetworks make some mismatching. On the other hand, in the case of global sparse synchronization, the population behavior is globally identical, independently of the cluster structure, because the intramodular dynamics of subnetworks make perfect matching. We introduce a realistic cross-correlation modularity measure, representing the matching degree between the instantaneous subpopulation spike rates of the subnetworks, and examine whether the sparse synchronization is global or modular. Depending on its magnitude, the intermodular coupling strength $J_{\mathrm{inter}}$ seems to play “dual” roles for the pacing between spikes in each subnetwork. For large $J_{\mathrm{inter}}$, due to strong inhibition it plays a destructive role to “spoil” the pacing between spikes, while for small $J_{\mathrm{inter}}$ it plays a constructive role to “favor” the pacing between spikes. Through competition between the constructive and the destructive roles of $J_{\mathrm{inter}}$, there exists an intermediate optimal $J_{\mathrm{inter}}$ at which the pacing degree between spikes becomes maximal. In contrast, the average number of intermodular links per interneuron $M_{\mathrm{syn}}^{\left(\mathrm{inter}\right)}$ seems to play a role just to favor the pacing between spikes. With increasing $M_{\mathrm{syn}}^{\left(\mathrm{inter}\right)}$, the pacing degree between spikes increases monotonically thanks to the increase in the degree of effectiveness of global communication between spikes. Furthermore, we employ the realistic sub- and whole-population order parameters, based on the instantaneous sub- and whole-population spike rates, to determine the threshold values for the synchronization-unsynchronization transition in the sub- and whole populations, and the degrees of global and modular sparse synchronization are also measured in terms of the realistic sub- and whole-population statistical-mechanical spiking measures defined by considering both the occupation and the pacing degrees of spikes. It is expected that our results could have implications for the role of the brain plasticity in some functional behaviors associated with population synchronization.

In this work we investigate a modified version of the Daisyworld model, originally introduced by Lovelock and Watson to describe in a simple way the interactions between an Earth-like planet, its biosphere, and the incoming solar radiation. Here a spatial dependency on latitude is included, and both a variable heat diffusivity along latitudes and a simple greenhouse effect description are introduced in the model. We show that the spatial interactions between the variables of the system can locally stabilize the coexistence of the two vegetation types. The feedback on albedo is able to generate equilibrium solutions which can efficiently self-regulate the planet climate, even for values of the solar luminosity relatively far from the current Earth conditions.

We investigate the survivor distributions of a spatially extended model of competitive dynamics in different geometries. The model consists of a deterministic dynamical system of individual agents at specified nodes, which might or might not survive the predatory dynamics: all stochasticity is brought in by the initial state. Every such initial state leads to a unique and extended pattern of survivors and nonsurvivors, which is known as an attractor of the dynamics. We show that the number of such attractors grows exponentially with system size, so that their exact characterization is limited to only very small systems. Given this, we construct an analytical approach based on inhomogeneous mean-field theory to calculate survival probabilities for arbitrary networks. This powerful (albeit approximate) approach shows how universality arises in survivor distributions via a key concept—the dynamical fugacity. Remarkably, in the large-mass limit, the survivor probability of a node becomes independent of network geometry and assumes a simple form which depends only on its mass and degree.

We examine Ising models with heat-bath dynamics on directed networks. Our simulations show that Ising models on directed triangular and simple cubic lattices undergo a phase transition that most likely belongs to the Ising universality class. On the directed square lattice the model remains paramagnetic at any positive temperature as already reported in some previous studies. We also examine random directed graphs and show that contrary to undirected ones, percolation of directed bonds does not guarantee ferromagnetic ordering. Only above a certain threshold can a random directed graph support finite-temperature ferromagnetic ordering. Such behavior is found also for out-homogeneous random graphs, but in this case the analysis of magnetic and percolative properties can be done exactly. Directed random graphs also differ from undirected ones with respect to zero-temperature freezing. Only at low connectivity do they remain trapped in a disordered configuration. Above a certain threshold, however, the zero-temperature dynamics quickly drives the model toward a broken symmetry (magnetized) state. Only above this threshold, which is almost twice as large as the percolation threshold, do we expect the Ising model to have a positive critical temperature. With a very good accuracy, the behavior on directed random graphs is reproduced within a certain approximate scheme.

In vertically vibrated dense suspensions, several localized structures have been discovered, such as heaps, stable holes, expanding holes, and replicating holes. Because an inclined free fluid surface is difficult to maintain because of gravitational pressure, the mechanism of those structures is not understood intuitively. In this paper, as a candidate for the driving mechanism, we focus on the boundary condition on a solid wall: the slip-nonslip switching boundary condition in synchronization with vertical vibration. By applying the lubrication approximation, we derived the time evolution equation of the fluid thickness from the Oldroyd-B fluid model. In our model we show that the initially flat fluid layer becomes unstable in a subcritical manner, and heaps and convectional flow appear. The obtained results are consistent with those observed experimentally. We also find that heaps climb a slope when the bottom is slightly inclined. We show that viscoelasticity enhances heap formation and climbing of a heap on the slope.

A high-order compact finite-difference lattice Boltzmann method (CFDLBM) is extended and applied to accurately simulate two-phase liquid-vapor flows with high density ratios. Herein, the He-Shan-Doolen-type lattice Boltzmann multiphase model is used and the spatial derivatives in the resulting equations are discretized by using the fourth-order compact finite-difference scheme and the temporal term is discretized with the fourth-order Runge-Kutta scheme to provide an accurate and efficient two-phase flow solver. A high-order spectral-type low-pass compact nonlinear filter is used to regularize the numerical solution and remove spurious waves generated by flow nonlinearities in smooth regions and at the same time to remove the numerical oscillations in the interfacial region between the two phases. Three discontinuity-detecting sensors for properly switching between a second-order and a higher-order filter are applied and assessed. It is shown that the filtering technique used can be conveniently adopted to reduce the spurious numerical effects and improve the numerical stability of the CFDLBM implemented. A sensitivity study is also conducted to evaluate the effects of grid size and the filtering procedure implemented on the accuracy and performance of the solution. The accuracy and efficiency of the proposed solution procedure based on the compact finite-difference LBM are examined by solving different two-phase systems. Five test cases considered herein for validating the results of the two-phase flows are an equilibrium state of a planar interface in a liquid-vapor system, a droplet suspended in the gaseous phase, a liquid droplet located between two parallel wettable surfaces, the coalescence of two droplets, and a phase separation in a liquid-vapor system at different conditions. Numerical results are also presented for the coexistence curve and the verification of the Laplace law. Results obtained are in good agreement with the analytical solutions and also the numerical results reported in the literature. The study shows that the present solution methodology is robust, efficient, and accurate for solving two-phase liquid-vapor flow problems even at high density ratios.

According to Rayleigh's criterion, rotating flows are linearly stable when their specific angular momentum increases radially outward. The celebrated magnetorotational instability opens a way to destabilize those flows, as long as the angular velocity is decreasing outward. Using a local approximation we demonstrate that even flows with very steep positive shear can be destabilized by azimuthal magnetic fields which are current free within the fluid. We illustrate the transition of this instability to a rotationally enhanced kink-type instability in the case of a homogeneous current in the fluid, and discuss the prospects for observing it in a magnetized Taylor-Couette flow.

The special macroscopic properties of liquid water stem from its structure as a complex network of molecules connected by hydrogen bonds. While the dynamics of single molecules within this network has been extensively investigated, only little attention has been paid to the closed loops (meshes) of hydrogen-bonded molecules which determine the network topology. Using molecular dynamics simulations we analyze the size, shape, geometrical arrangement, and dynamical stability of loops containing up to 10 hydrogen bonds. We find that six-membered loops in liquid water even at room temperature retain a striking similarity with the well-known structure of ice. Analyzing the network dynamics we find that rings of more than five hydrogen bonds are stabilized compared to a random collection containing the same number of single bonds. We finally show that in the vicinity of hydrophobic and hydrophilic interfaces loops arrange in a preferred orientation.

We study the structure of quasiperiodic Lorentz gases, i.e., particles bouncing elastically off fixed obstacles arranged in quasiperiodic lattices. By employing a construction to embed such structures into a higher-dimensional periodic hyperlattice, we give a simple and efficient algorithm for numerical simulation of the dynamics of these systems. This same construction shows that quasiperiodic Lorentz gases generically exhibit a regime with infinite horizon, that is, empty channels through which the particles move without colliding, when the obstacles are small enough; in this case, the distribution of free paths is asymptotically a power law with exponent $-3$, as expected from infinite-horizon periodic Lorentz gases. For the critical radius at which these channels disappear, however, a new regime with locally finite horizon arises, where this distribution has an unexpected exponent of $-5$, previously observed only in a Lorentz gas formed by superposing three incommensurable periodic lattices in the Boltzmann-Grad limit where the radius of the obstacles tends to zero.

We investigate the nonequilibrium behavior of dense, plastic-crystalline suspensions of mildly anisotropic colloidal hard dumbbells under the action of an oscillatory shear field by employing Brownian dynamics computer simulations. In particular, we extend previous investigations, where we uncovered nonequilibrium phase transitions, to other aspect ratios and to a larger nonequilibrium parameter space, that is, a wider range of strains and shear frequencies. We compare and discuss selected results in the context of scattering and rheological experiments. Both simulations and experiments demonstrate that the previously found transitions from the plastic crystal phase with increasing shear strain also occur at other aspect ratios. We explore the transition behavior in the strain-frequency phase and summarize it in a nonequilibrium phase diagram. Additionally, the experimental rheology results hint at a slowing down of the colloidal dynamics with higher aspect ratio.

We employ molecular dynamics simulations to study the wetting and evaporation of salt-water nanodroplets on platinum surfaces. Our results show that the contact angle of the droplets increases with the salt concentration. To verify this, a second simulation system of a thin salt-water film on a platinum surface is used to calculate the various surface tensions. We find that both the solid-liquid and liquid-vapor surface tensions increase with salt concentration and as a result these cause an increase in the contact angle. However, the evaporation rate of salt-water droplets decreases as the salt concentration increases, due to the hydration of salt ions. When the water molecules have all evaporated from the droplet, two forms of salt crystals are deposited, clump and ringlike, depending on the solid-liquid interaction strength and the evaporation rate. To form salt crystals in a ring, it is crucial that there is a pinned stage in the evaporation process, during which salt ions can move from the center to the rim of the droplets. With a stronger solid-liquid interaction strength, a slower evaporation rate, and a higher salt concentration, a complete salt crystal ring can be deposited on the surface.