We investigate the vibrational modes of quasi-two-dimensional disordered colloidal packings of hard colloidal spheres with short-range attractions as a function of packing fraction. Certain properties of the vibrational density of states (vDOS) are shown to correlate with the density and structure of the samples (i.e., in sparsely versus densely packed samples). Specifically, a crossover from dense glassy to sparse gel-like states is suggested by an excess of phonon modes at low frequency and by a variation in the slope of the vDOS with frequency at low frequency. This change in phonon mode distribution is demonstrated to arise largely from localized vibrations that involve individual and/or small clusters of particles with few local bonds. Conventional order parameters and void statistics did not exhibit obvious gel-glass signatures as a function of volume fraction. These mode behaviors and accompanying structural insights offer a potentially new set of indicators for identification of glass-gel transitions and for assignment of gel-like versus glass-like character to a disordered solid material.

A simple model is constructed to study the phase diagram and thermodynamic properties of methanol, which is described as a dimer of an apolar sphere mimicking the methyl group and a sphere with core-softened potential as the hydroxyl group. Performing classical Monte Carlo simulations, we obtained the phase diagram, showing a second critical point between two different liquid phases. Evaluating systems with a different number of particles, we extrapolate to infinite size in accordance with Ising universality class to obtain bulk values for critical temperature, pressure, and density. Strong evidence that the structure of the liquid changes upon transition from high- to low-density phase was provided. From the experimentally determined hydrogen bond strength and length in methanol and water, we propose where the second critical point of methanol should be.

The temporal fluctuations in the near plume of a 100-A ${\mathrm{LaB}}_6$ hollow cathode are experimentally investigated. A probe array is employed to measure the amplitude and dispersion of axial modes in the plume, and these properties are examined parametrically as a function of cathode operating conditions. The onset of ion acoustic turbulence is observed at high current and is characterized by a power spectrum that exhibits a cutoff at low frequency and an inverse dependence on frequency at high values. The amplitude of the turbulence is found to decrease with flow rate but to depend nonmonotonically on discharge current. Estimates of the anomalous collision frequency based on experimental measurements indicate that the ion acoustic turbulence collision frequency can exceed the classical rate at high discharge current densities by nearly two orders of magnitude.

The local approach to computing electrostatic interactions proposed by Maggs and adapted by Rottler and Pasichnyk for molecular-dynamics simulations is extended to situations where the dielectric background medium is inhomogeneous. We furthermore correct a problem of the original algorithm related to the correct treatment of the global dipole moment, provide an error estimate for the accuracy of the algorithm, and suggest a different form of the treatment of the self-energy problem. Our implementation is highly scalable on many cores, and we have validated and compared its performance against theoretical predictions and simulation data obtained by other algorithmic approaches.

Distinctive orderings and phase diagram structures are found, from renormalization-group theory, for odd $q$-state clock spin-glass models in $d=3$ dimensions. These models exhibit asymmetric phase diagrams, as is also the case for quantum Heisenberg spin-glass models. No finite-temperature spin-glass phase occurs. For all odd $q\ge 5$, algebraically ordered antiferromagnetic phases occur. One such phase is dominant and occurs for all $q\ge 5$. Other such phases occupy small low-temperature portions of the phase diagrams and occur for $5\le q\le 15$. All algebraically ordered phases have the same structure, determined by an attractive finite-temperature sink fixed point where a dominant and a subdominant pair states have the only nonzero Boltzmann weights. The phase transition critical exponents quickly saturate to the high $q$ value.

A study of the multifractal characteristics of the structure formed by magnetic particles in a dilute magnetorheological fluid is presented. A quasi-two-dimensional magnetorheological fluid sample is simultaneously subjected to a static magnetic field and a sinusoidal magnetic field transverse to each other. We analyzed the singularity spectrum $f\left(\alpha \right)$ and the generalized dimension $D\left(q\right)$ of the whole structure to characterize the distribution of the aggregates under several conditions of particle concentration, magnetic field intensities, and liquid viscosity. We also obtained the fractal dimension $D_g$, calculated from the radius of gyration of the chains, to describe the internal distribution of the particles. We present a thermodynamic interpretation of the multifractal analysis, and based on this, we discussed the characteristics of the structure formed by the particles and its relation with previous studies of the average chain length. We have found that this method is useful to quantitatively describe the structure of magnetorheological fluids, especially in systems with high particle concentration where the aggregates are more complex than simple chains or columns.

We extend the kinematic matrix (“kinematrix”) formalism [Phys. Rev. E 89, 062304 (2014).], which via simple matrix algebra accesses ensemble properties of self-propellers influenced by uncorrelated noise, to treat Gaussian correlated noises. This extension brings into reach many real-world biological and biomimetic self-propellers for which inertia is significant. Applying the formalism, we analyze in detail ensemble behaviors of a 2D self-propeller with velocity fluctuations and orientation evolution driven by an Ornstein-Uhlenbeck process. On the basis of exact results, a variety of dynamical regimes determined by the inertial, speed-fluctuation, orientational diffusion, and emergent disorientation time scales are delineated and discussed.

Nematic and cholesteric liquid crystals are three-dimensional fluids that possess long-range orientational ordering and can support both topological defects and chiral superstructures. Implications of this ordering remain unexplored even for simple dynamic processes such as the ones found in so-called “fall experiments,” or motion of a spherical inclusion under the effects of gravity. Here we show that elastic and surface anchoring interactions prompt periodic dynamics of colloidal microparticles in confined cholesterics when gravity acts along the helical axis. We explore elastic interactions between colloidal microparticles and confining surfaces as well as with an aligned ground-state helical structure of cholesterics for different sizes of spheres relative to the cholesteric pitch, demonstrating unexpected departures from Stokes-like behavior at very low Reynolds numbers. We characterize metastable localization of microspheres under the effects of elastic and surface anchoring periodic potential landscapes seen by moving spheres, demonstrating the important roles played by anchoring memory, confinement, and topological defect transformation. These experimental findings are consistent with the results of numerical modeling performed through minimizing the total free energy due to colloidal inclusions at different locations along the helical axis and with respect to the confining substrates. A potential application emerging from this work is colloidal sorting based on particle shapes and sizes.

We use a molecular-motor model previously proposed for a nematic cell with an azo-dye monolayer to calculate the director orientation when light is normally impinged on the cell. We consider an initial planar configuration for which one of the surrounding plates, which we call the reference plate, is submitted to a hard-anchoring boundary condition. The other confining plate has a coating monolayer of azo-dye molecules such that the change of the orientation of azo-dye isomers, due to light, causes changes in the nematic director. The boundary conditions on both plates along with the optical field determine the director configuration in the bulk. The existence of periodic solutions for the density of isomers in trans and cis states, corresponding to weak optical fields, has been discussed in the literature. Using a similar approach, we find an approximate expression for the density of isomers, written in terms of the director angle, which allows us to close the equation for the director configuration on the boundary having a photosensitive plate. We decouple the director's angle and the isomer densities by assuming extremely different temporal time scales between them. We show for a given sample that switching times inversely depend on the trans-cis transition rate of photoexcitation whereas relaxation times do not depend on it. On the other hand, switching and relaxation times linearly depend on effective surface viscosity values. Our model allows us to estimate surface viscosity values.

The identification of early warning signals for regime shifts in ecosystems is of crucial importance given their impact in terms of economic and social effects. We present here the results of a theoretical study on the desertification transition in semiarid ecosystems under external stress. We performed numerical simulations based on a stochastic cellular automaton model, and we studied the dynamics of the vegetation clusters in terms of percolation theory, assumed as an effective tool for analyzing the geometrical properties of the clusters. Focusing on the role played by the strength of external stresses, measured by the mortality rate $m$, we followed the progressive degradation of the ecosystem for increasing $m$, identifying different stages: first, the fragmentation transition occurring at relatively low values of $m$, then the desertification transition at higher mortality rates, and finally the full desertification transition corresponding to the extinction of the vegetation and the almost complete degradation of the soil, attained at the maximum value of $m$. For each transition we calculated the spanning probabilities as functions of $m$ and the percolation thresholds according to different spanning criteria. The identification of the different thresholds is proposed as an useful tool for monitoring the increasing degradation of real-world finite-size systems. Moreover, we studied the time fluctuations of the sizes of the biggest clusters of vegetated and nonvegetated cells over the entire range of mortality values. The change of sign in the skewness of the size distributions, occurring at the fragmentation threshold for the biggest vegetation cluster and at the desertification threshold for the nonvegetated cluster, offers new early warning signals for desertification. Other new and robust indicators are given by the maxima of the root-mean-square deviation of the distributions, which are attained respectively inside the fragmentation interval, for the vegetated biggest cluster, and inside the desertification interval, for the nonvegetated cluster.

We suggest an approach to probing intermittency corrections to the Kolmogorov law in turbulent flows based on the autoregressive moving-average modeling of turbulent time series. We introduce an index ${\rm Y}$ that measures the distance from a Kolmogorov-Obukhov model in the autoregressive moving-average model space. Applying our analysis to particle image velocimetry and laser Doppler velocimetry measurements in a von Kármán swirling flow, we show that ${\rm Y}$ is proportional to traditional intermittency corrections computed from structure functions. Therefore, it provides the same information, using much shorter time series. We conclude that ${\rm Y}$ is a suitable index to reconstruct intermittency in experimental turbulent fields.

Quantum feedback networks have been introduced in quantum optics as a framework for constructing arbitrary networks of quantum mechanical systems connected by unidirectional quantum optical fields, and has allowed for a system theoretic approach to open quantum optics systems. Our aim here is to establish a network theory for quantum transport systems where typically the mediating fields between systems are bidirectional. Mathematically, this leads us to study quantum feedback networks where fields arrive at ports in input-output pairs, making it a special case of the unidirectional theory where inputs and outputs are paired. However, it is conceptually important to develop this theory in the context of quantum transport theory—the resulting theory extends traditional approaches which tend to view the components in quantum transport as scatterers for the various fields, in the process allowing us to consider emission and absorption of field quanta by these components. The quantum feedback network theory is applicable to both Bose and Fermi fields, moreover, it applies to nonlinear dynamics for the component systems. We advance the general theory, but study the case of linear passive quantum components in some detail.

We study the asymptotic and preasymptotic diffusive properties of Brownian particles in channels whose section varies periodically in space. The effective diffusion coefficient $D_{\mathrm{eff}}$ is numerically determined by the asymptotic behavior of the root mean square displacement in different geometries, considering even cases of steep variations of the channel boundaries. Moreover, we compared the numerical results to the predictions from the various corrections proposed in the literature to the well known Fick-Jacobs approximation. Building an effective one-dimensional equation for the longitudinal diffusion, we obtain an approximation for the effective diffusion coefficient. Such a result goes beyond a perturbation approach, and it is in good agreement with the actual values obtained by the numerical simulations. We discuss also the preasymptotic diffusion which is observed up to a crossover time whose value, in the presence of strong spatial variation of the channel cross section, can be very large. In addition, we show how the Einstein's relation between the mean drift induced by a small external field and the mean square displacement of the unperturbed system is valid in both asymptotic and preasymptotic regimes.

We investigate the totally asymmetric simple exclusion process on closed and directed random regular networks, which is a simple model of active transport in the one-dimensional segments coupled by junctions. By a pair mean-field theory and detailed numerical analyses, it is found that the correlations at junctions induce two notable deviations from the simple mean-field theory, which neglects these correlations: (1) the narrower range of particle density for phase coexistence and (2) the algebraic decay of density profile with exponent $1/2$ even outside the maximal-current phase. We show that these anomalies are attributable to the effective slow bonds formed by the network junctions.

Turbulent transport in magnetic fusion plasmas can be significantly suppressed by Reynolds-stress-induced zonal flows, allowing effective plasma confinement. We present experimental evidence of spatiotemporal correlation between small-scale turbulence-induced Reynolds stress and large-scale zonal flow production in the $\mathbf{E}\times \mathbf{B}$ driven hydrodynamic spectral condensation. We show that Reynolds stress is generated effectively by anisotropic vorticity structures possessing collective tilt angle. The maximum amplitude of the tilt, the Reynolds stress, and the mean zonal flow production coincide with the transition time of the velocity field, indicating a key role of turbulence-induced Reynolds stress in the condensation of the flow. The analysis of the energy transfer between turbulence and zonal flow shows coherent oscillations with $\pi /2$ phase delay, thus indicating a predator-prey-like interaction between zonal flow and turbulence.

We investigate the dynamics of localized solutions of the relativistic cold-fluid plasma model in the small but finite amplitude limit, for slightly overcritical plasma density. Adopting a multiple scale analysis, we derive a perturbed nonlinear Schrödinger equation that describes the evolution of the envelope of circularly polarized electromagnetic field. Retaining terms up to fifth order in the small perturbation parameter, we derive a self-consistent framework for the description of the plasma response in the presence of localized electromagnetic field. The formalism is applied to standing electromagnetic soliton interactions and the results are validated by simulations of the full cold-fluid model. To lowest order, a cubic nonlinear Schrödinger equation with a focusing nonlinearity is recovered. Classical quasiparticle theory is used to obtain analytical estimates for the collision time and minimum distance of approach between solitons. For larger soliton amplitudes the inclusion of the fifth-order terms is essential for a qualitatively correct description of soliton interactions. The defocusing quintic nonlinearity leads to inelastic soliton collisions, while bound states of solitons do not persist under perturbations in the initial phase or amplitude.

We study the Brownian motion of microbeads immersed in water and in a viscoelastic wormlike micelles solution by optical trapping interferometry and diffusing wave spectroscopy. Through the mean-square displacement obtained from both techniques, we deduce the mechanical properties of the fluids at high frequencies by explicitly accounting for inertia effects of the particle and the surrounding fluid at short time scales. For wormlike micelle solutions, we recover the 3/4 scaling exponent for the loss modulus over two decades in frequency as predicted by the theory for semiflexible polymers.

We study the pinning dynamics of magnetic flux (vortex) lines in a disordered type-II superconductor. Using numerical simulations of a directed elastic line model, we extract the pinning time distributions of vortex line segments. We compare different model implementations for the disorder in the surrounding medium: discrete, localized pinning potential wells that are either attractive and repulsive or purely attractive, and whose strengths are drawn from a Gaussian distribution; as well as continuous Gaussian random potential landscapes. We find that both schemes yield power-law distributions in the pinned phase as predicted by extreme-event statistics, yet they differ significantly in their effective scaling exponents and their short-time behavior.

We numerically simulate the diffusion of overdampd pointlike Janus particles along narrow two-dimensional periodically corrugated channels with reflecting walls. The self-propulsion velocity of the particle is assumed to rotate subject to an intrinsic bias modeled by a torque. Breaking the mirror symmetry of the channel with respect to its axis suffices to generate a directed particle flow with orientation and magnitude which depend on the channel geometry and the particle swimming properties. This means that chiral microswimmers drift autonomously along a narrow channel under more general asymmetry conditions than previously reported, a property of potential impact on their fabrication and technological applications.