In this work, electrostatic forces between charged nanogels are explored through coarse-grained simulations. These simulations allow us to explicitly consider the complex topology of these nanoparticles and provide reliable force values to examine highly charged nanogels of a few tens of nanometers. The results obtained here clearly reveal that the electrostatic interactions between these nanoparticles are not governed by the net charge of the nanogel, which includes not only the charge of the polymer network but also the charge of ions inside. Thus two theoretical procedures for predicting effective charges are also proposed and investigated. Both provide predictions of the same order and capture the behavior found for the effective charge obtained from simulations.

Self-propelled particles phase-separate into coexisting dense and dilute regions above a critical density. The statistical nature of their stochastic motion lends itself to various theories that predict the onset of phase separation. However, these theories are ill-equipped to describe such behavior when noise becomes negligible. To overcome this limitation, we present a predictive model that relies on two density-dependent timescales: ${\tau }_F$, the mean time particles spend between collisions; and ${\tau }_C$, the mean lifetime of a collision. We show that only when ${\tau }_F<{\tau }_C$ do collisions last long enough to develop a growing cluster and initiate phase separation. Using both analytical calculations and active particle simulations, we measure these timescales and determine the critical density for phase separation in both two and three dimensions.

We obtain numerical solutions of the two-fluid fluctuating compressible Navier-Stokes (FCNS) equations, which consistently account for thermal fluctuations from meso- to macroscales, in order to study the effect of such fluctuations on the mixing behavior in the Richtmyer-Meshkov instability (RMI). The numerical method used was successfully verified in two stages: for the deterministic fluxes by comparison against air-${\mathrm{SF}}_6$ RMI experiment, and for the stochastic terms by comparison against the direct simulation Monte Carlo results for He-Ar RMI. We present results from fluctuating hydrodynamic RMI simulations for three He-Ar systems having length scales with decreasing order of magnitude that span from macroscopic to mesoscopic, with different levels of thermal fluctuations characterized by a nondimensional Boltzmann number (Bo). For a multidimensional FCNS system on a regular Cartesian grid, when using a discretization of a space-time stochastic flux $Z(x,t)$ of the form $Z(x,t)\to 1/\sqrt{h▵t}N(ih,n\mathrm{\Delta }t)$ for spatial interval $h$, time interval $\mathrm{\Delta }t$, $h$, and Gaussian noise $N$ should be greater than $h_0$, with $h_0$ corresponding to a cell volume that contains a sufficient number of molecules of the fluid such that the fluctuations are physically meaningful and produce the right equilibrium spectrum. For the mesoscale RMI systems simulated, it was desirable to use a cell size smaller than this limit in order to resolve the viscous shock. This was achieved by using a modified regularization of the noise term via $Z\left(x,t\right)\to 1/\sqrt{▵tmax\left(h^3,h_0^3\right)}N\left(ih,n\mathrm{\Delta }t\right)$, with $h_0=\xi h\forall h<h_0,\xi \in {\mathbb{Z}}^{+}$. Our simulations show that for systems with $\text{Bo}\ll 1$ deterministic mixing behavior emerges as the ensemble-averaged behavior of several fluctuating instances, whereas when $\text{Bo}\approx 1$, a deviation from deterministic behavior is observed. For all cases, the FCNS solution provides bounds on the growth rate of the amplitude of the mixing layer.

A computational fluid dynamic investigation has been carried out to study the dynamics of a moving compound droplet inside a tube. The motions associated with such a droplet is uncovered by solving the axisymmetric Navier-Stokes equations in which the spatiotemporal evolution of a pair of twin-deformable interfaces has been tracked employing the volume-of-fluid approach. The deformations at the interfaces and their subsequent dynamics are found to be stimulated by the subtle interplay between the capillary and viscous forces. The simulations uncover that when a compound drop composed of concentric inner and outer interfaces migrates inside a tube, initially in the unsteady domain of evolution, the inner drop shifts away from the concentric position to reach a morphology of constant eccentricity at the steady state. The coupled motions of the droplets in the unsteady regime causes a continuous deformation of the inner and outer interfaces to obtain a configuration with a (an) prolate (oblate) shaped outer (inner) interface. The magnitudes of capillary number and viscosity ratio are found to have significant influence on the temporal evolution of the interfacial deformations as well as the eccentricity of the droplets. Further, the simulations uncover that, following the asymmetric deformation of the interfaces, the migrating compound droplet can undergo an uncommon breakup stimulated by a rather irregular pinch-off of the outer shell. The breakup is found to initiate with the thinning of the outer shell followed by the pinch-off. Interestingly, the kinetics of the thinning of outer shell is found to follow two distinct power-law regimes—a swiftly thinning stage at the onset followed by a rate limiting stage before pinch-off, which eventually leads to the uncommon breakup of the migrating compound droplets.

It is widely assumed that thermionic current flow through a plasma is limited by a “space-charge-limited” (SCL) cathode sheath that consumes the hot cathode's negative bias and accelerates upstream ions into the cathode. Here, we formulate a fundamentally different current-limited mode. In the “inverse” mode, the potentials of both electrodes are above the plasma potential, so that the plasma ions are confined. The bias is consumed by the anode sheath. There is no potential gradient in the neutral plasma region from resistivity or presheath. The inverse cathode sheath pulls some thermoelectrons back to the cathode, thereby limiting the circuit current. Thermoelectrons entering the zero-field plasma region that undergo collisions may also be sent back to the cathode, further attenuating the circuit current. In planar geometry, the plasma density is shown to vary linearly across the electrode gap. A continuum kinetic planar plasma diode simulation model is set up to compare the properties of current modes with classical, conventional SCL, and inverse cathode sheaths. SCL modes can exist only if charge-exchange collisions are turned off in the potential well of the virtual cathode to prevent ion trapping. With the collisions, the current-limited equilibrium must be inverse. Inverse operating modes should therefore be present or possible in many plasma devices that rely on hot cathodes. Evidence from past experiments is discussed. The inverse mode may offer opportunities to minimize sputtering and power consumption that were not previously explored due to the common assumption of SCL sheaths.

We consider the Brownian motion of a collection of particles each with an additional degree of freedom. The degree of freedom of a particle (or, in general, a molecule) can assume distinct values corresponding to certain states or conformations. The time evolution of the additional degree of freedom of a particle is guided by those of its neighbors as well as the temperature of the system. We show that the local averaging over these degrees of freedom results in emergence of a collective order in the dynamics in the form of selection or dominance of one of the isomers leading to a symmetry-broken state. Our statistical model captures the basic features of homochirality, e.g., autocatalysis and chiral inhibition.

We design an irreversible worm algorithm for the zero-field ferromagnetic Ising model by using the lifting technique. We study the dynamic critical behavior of an energylike observable on both the complete graph and toroidal grids, and compare our findings with reversible algorithms such as the Prokof'ev-Svistunov worm algorithm. Our results show that the lifted worm algorithm improves the dynamic exponent of the energylike observable on the complete graph and leads to a significant constant improvement on toroidal grids.

Quantum interference is at the heart of what sets the quantum and classical worlds apart. We demonstrate that quantum interference effects involving a many-body working medium is responsible for genuinely nonclassical features in the performance of a quantum heat engine. The features with which quantum interference manifests itself in the work output of the engine depends strongly on the extensive nature of the working medium. While identifying the class of work substances that optimize the performance of the engine, our results shed light on the optimal size of such media of quantum workers to maximize the work output and efficiency of quantum energy machines.

The self-avoiding walk on the square site-diluted correlated percolation lattice is considered. The Ising model is employed to realize the spatial correlations of the metric space. As a well-accepted result, the (generalized) Flory's mean-field relation is tested to measure the effect of correlation. After exploring a perturbative Fokker-Planck-like equation, we apply an enriched Rosenbluth Monte Carlo method to study the problem. To be more precise, the winding angle analysis is also performed from which the diffusivity parameter of Schramm-Loewner evolution theory ($\kappa$) is extracted. We find that at the critical Ising (host) system, the exponents are in agreement with Flory's approximation. For the off-critical Ising system, we find also a behavior for the fractal dimension of the walker trace in terms of the correlation length of the Ising system $\xi \left(T\right)$, i.e., $D_F^{\text{SAW}}\left(T\right)-D_F^{\text{SAW}}\left(T_c\right)\sim \frac{1}{\sqrt{\xi \left(T\right)}}$.

We apply a recently developed semiclassical theory of short periodic orbits to the continuously open quantum tribaker map. In this paradigmatic system the trajectories are partially bounced back according to continuous reflectivity functions. This is relevant in many situations that include optical microresonators and more complicated boundary conditions. In a perturbative regime, the shortest periodic orbits belonging to the classical repeller of the open map—a cantor set given by a region of exactly zero reflectivity—prove to be extremely robust in supporting a set of long-lived resonances of the continuously open quantum maps. Moreover, for steplike functions a significant reduction in the number needed is obtained, similarly to the completely open situation. This happens despite a strong change in the spectral properties when compared to the discontinuous reflectivity case. In order to give a more realistic interpretation of these results we compare with a Fresnel-type reflectivity function.

The mechanisms driving the collective movement of cells remain poorly understood. To contribute toward resolving this mystery, a model was formulated to theoretically explore the possible functions of polarized cell-cell adhesion in collective cell migration. The model consists of an amoeba cell with polarized cell-cell adhesion, which is controlled by positive feedback with cell motion. This model cell has no persistent propulsion and therefore exhibits a simple random walk when in isolation. However, at high density, these cells acquire collective propulsion and form ordered movement. This result suggests that cell-cell adhesion has a potential function, which induces collective propulsion with persistence.

We present the effects of structure and stiffness of block copolymers on the interfacial properties of an immiscible homopolymer blend. Diblock and two-arm grafted copolymers with variation in stiffness are modeled using coarse-grained molecular dynamics to compare the compatibilization efficiency, i.e., reduction of interfacial tension. Overall, grafted copolymers are located more compactly at the interface and show better compatibilization efficiency than diblock copolymers. In addition, an increase in the stiffness for one of the blocks of the diblock copolymers causes unusual inhomogeneous interfacial coverage due to bundle formation. However, an increase in the stiffness for one of blocks of the grafted copolymers prevents the bundle formation due to the branched chain. As a result, homogeneous interfacial coverage of homopolymer blends is realized with significant reduction of interfacial tension which makes grafted copolymer a better candidate for the compatibilizer of immiscible homopolymer blend.

Rheology of a dilute cohesive granular gas is theoretically and numerically studied. The flow curve between the shear viscosity and the shear rate is derived from the inelastic Boltzmann equation for particles having square-well potentials in a simple shear flow. It is found that (i) the stable uniformly sheared state only exists above a critical shear rate and (ii) the viscosity in the uniformly sheared flow is almost identical to that for uniformly sheared flow of hard core granular particles. Below the critical shear rate, clusters grow with time, in which the viscosity can be approximated by that for the hard-core fluids if we replace the diameter of the particle by the mean diameter of clusters.

We carry out a direct comparison of experimental and numerical realizations of the exact same granular system as it undergoes shear jamming. We adjust the numerical methods used to optimally represent the experimental settings and outcomes up to microscopic contact force dynamics. Measures presented here range from microscopic through mesoscopic to systemwide characteristics of the system. Topological properties of the mesoscopic force networks provide a key link between microscales and macroscales. We report two main findings: (1) The number of particles in the packing that have at least two contacts is a good predictor for the mechanical state of the system, regardless of strain history and packing density. All measures explored in both experiments and numerics, including stress-tensor-derived measures and contact numbers depend in a universal manner on the fraction of nonrattler particles, $f_{\mathrm{NR}}$. (2) The force network topology also tends to show this universality, yet the shape of the master curve depends much more on the details of the numerical simulations. In particular we show that adding force noise to the numerical data set can significantly alter the topological features in the data. We conclude that both $f_{\mathrm{NR}}$ and topological metrics are useful measures to consider when quantifying the state of a granular system.

The full melting of a two-dimensional plasma crystal was induced in a principally stable monolayer by localized laser stimulation. Two distinct behaviors of the crystal after laser stimulation were observed depending on the amount of injected energy: (i) below a well-defined threshold, the laser melted area recrystallized; (ii) above the threshold, it expanded outwards in a similar fashion to mode-coupling instability-induced melting, rapidly destroying the crystalline order of the whole complex plasma monolayer. The reported experimental observations are due to the fluid mode-coupling instability, which can pump energy into the particle monolayer at a rate surpassing the heat transport and damping rates in the energetic localized melted spot, resulting in its further growth. This behavior exhibits remarkable similarities with impulsive spot heating in ordinary reactive matter.

We theoretically study artificial light harvesting by a Möbius ring. When the donors in the ring are dimerized, the energies of the donor ring are split into two subbands. Because of the nontrivial Möbius boundary condition, both the photon and acceptor are coupled to all collective-excitation modes in the donor ring. Therefore, the quantum dynamics in the light harvesting is subtly influenced by dimerization in the Möbius ring. It is discovered that energy transfer is more efficient in a dimerized ring than that in an equally spaced ring. This discovery is also confirmed by a calculation with the perturbation theory, which is equivalent to the Wigner-Weisskopf approximation. Our findings may be beneficial to the optimal design of artificial light harvesting.

Synchronization and stability under periodic oscillatory driving are well understood, but little is known about the effects of aperiodic driving, despite its abundance in nature. Here, we consider oscillators subject to driving with slowly varying frequency, and investigate both short-term and long-term stability properties. For a phase oscillator, we find that, counterintuitively, such variation is guaranteed to enlarge the Arnold tongue in parameter space. Using analytical and numerical methods that provide information on time-variable dynamical properties, we find that the growth of the Arnold tongue is specifically due to the growth of a region of intermittent synchronization where trajectories alternate between short-term stability and short-term neutral stability, giving rise to stability on average. We also present examples of higher-dimensional nonlinear oscillators where a similar stabilization phenomenon is numerically observed. Our findings help support the case that in general, deterministic nonautonomous perturbation is a very good candidate for stabilizing complex dynamics.

The present paper proposes a scheme for controlling wave segments in excitable media. This scheme consists of two phases: in the first phase, a simple mathematical model for wave segments is derived using only the time series data of input and output signals for the media; in the second phase, the model derived in the first phase is used in an advanced control technique. We demonstrate with numerical simulations of the Oregonator model that this scheme performs better than a conventional control scheme.

Analyzing the long-term behaviors (attractors) of dynamic models of biological networks can provide valuable insight. We propose a general method that can find the attractors of multilevel discrete dynamical systems by extending a method that finds the attractors of a Boolean network model. The previous method is based on finding stable motifs, subgraphs whose nodes’ states can stabilize on their own. We extend the framework from binary states to any finite discrete levels by creating a virtual node for each level of a multilevel node, and describing each virtual node with a quasi-Boolean function. We then create an expanded representation of the multilevel network, find multilevel stable motifs and oscillating motifs, and identify attractors by successive network reduction. In this way, we find both fixed point attractors and complex attractors. We implemented an algorithm, which we test and validate on representative synthetic networks and on published multilevel models of biological networks. Despite its primary motivation to analyze biological networks, our motif-based method is general and can be applied to any finite discrete dynamical system.

Net fluid flow of soap films induced by (ac or dc) electric fields in asymmetric frames is presented. Previous experiments of controllable soap film flow required the simultaneous use of an electrical current passing through the film and an external electric field or the use of nonuniform ac electric fields. Here a single voltage difference generates both the electrical current going through the film and the electric field that actuates on the charge induced on the film. The film is set into global motion due to the broken symmetry that appears by the use of asymmetric frames. If symmetric frames are used, the film flow is not steady but time dependent and irregular. Finally, we study numerically these film flows by employing the model of charge induction in ohmic liquids.