The relationship between macroscopic shear yield criteria and local stress distributions in deformed polymer glasses is investigated via molecular dynamics simulations on different scales of coarse-graining. Macroscopic shear stresses at the yield point obey a pressure-modified von Mises (pmvM) criterion for many different loading conditions and strain rates. Average local stresses in small volume elements obey the same yield criterion for volumes containing approx. 100 atoms or more. Qualitatively different behavior is observed on smaller scales: the average octahedral atomic shear stress has a simple linear relationship to hydrostatic pressure regardless of macroscopic stress state and failure mode. Local plastic events are identified through a threshold in the mean-squared nonaffine displacement and compared to the local stress state. We find that the pmvM criterion only predicts local yield events when stress and displacements are averaged over at least 100 atoms. By contrast, macroscopic shear yield criteria appear to lose their ability to predict plastic activity on the atomic scale.

We study the system of cylindrical colloids adhering to an originally flat fluid membrane, on the basis of a full treatment of the Helfrich model. Our approach allows for numerical calculation of the free energy in both shallow- and deep-wrapping conformations. We show that the free energy of two cylinders adhering to the same side of a membrane has two branches corresponding to shallow and deep wrapping and that the system of two cylinders adhering to opposite sides of a membrane can undergo a first-order phase transition between two membrane-mediated attractive states.

We determine the particulate transport properties of fluid membranes with nontrivial geometries that are surrounded by viscous Newtonian solvents. Previously, this problem in membrane hydrodynamics was discussed for the case of flat membranes by Saffman and Delbrück [P. G. Saffman and M. Delbrück, Proc. Natl. Acad. Sci. U.S.A. 72, 3111 (1975)]. We review and develop the formalism necessary to consider the hydrodynamics of membranes with arbitrary curvature and show that the effect of local geometry is twofold. First, local Gaussian curvature introduces in-plane viscous stresses even for situations in which the velocity field is coordinate-independent. Secondly, even in the absence of Gaussian curvature, the geometry of the membrane modifies the momentum transport between the bulk fluids and the membrane. We illustrate these effects by examining in detail the mobilities of particles bound to spherical and cylindrical membranes. These two examples provide experimentally testable predictions for particulate mobilities and membrane velocity fields on giant unilamellar vesicles and membrane tethers. Finally, we use the examples of spherical and cylindrical membranes to demonstrate how the global geometry and topology of the membrane influences the membrane velocities and the particle mobilities.

In a recent experiment [Müller et al., Phys. Rev. A 79, 031804(R) (2009)] reported a splitting of the stochastic resonance peak, which they attributed to the asymmetry of an effective double-well restoring potential in their optomechanical read-out device. We show here that such an effect, though smaller than reported, is indeed consistent with a characterization of stochastic resonance as a synchronization phenomenon, while it proves elusive in terms of spectral quantifiers.

We study the effect of quantum fluctuations on the critical behavior of the Ising ferromagnetic phase transitions that do not belong to the mean-field universality class. A model system is considered, in which Ising spins are placed on the nodes of a scale-free network. Our Monte Carlo analysis shows that the critical exponents differ from those of mean-field phase transitions when degree exponent $\gamma$ is in the range $3<\gamma <5$. This confirms earlier analytic calculations based on ansatzes and approximation methods. As we apply quantum fluctuations by means of a magnetic field perpendicular to the Ising spin direction, the transition temperature $T_c$ decreases with increasing magnetic field strength. We find, however, that the quantum fluctuations do not alter the critical exponents and the universality class remains unchanged.

On the basis of the tokamap, characteristic features of magnetic field lines and zeroth-order guiding-center particle motion in the whole body of a magnetically confined plasma, e.g., a tokamak plasma, are investigated. It is shown that the tokamap exhibits a poloidal transport that can be considered as a Hamiltonian ratchet. In a situation with partially chaotic magnetic field lines the locking of the averaged poloidal velocity occurs to a value that does not depend precisely on the initial conditions. The so-called sum rule predicts the mean velocity in agreement with the observed magnitude. Possible consequences for the onset of poloidal plasma rotation in ergodized plasmas are elucidated.

Harmonic moments are integrals of integer powers of $z=x+iy$ over a domain. Here, the domain is an exterior of a bubble of air growing in an oil layer between two horizontal closely spaced plates. Harmonic moments are a natural basis for such Laplacian growth phenomena because, unlike other representations, these moments linearize the zero surface tension problem [S. Richardson, J. Fluid Mech. 56, 609 (1972)], so that all moments except the lowest one (the area of the bubble) are conserved in time. In our experiments, we directly determine the harmonic moments and show that for nonzero surface tension, all moments (except the lowest one) decay in time rather than exhibiting the divergences of other representations. Further, we derive an expression that relates the derivative of the $k^{th}$ harmonic moment $M_k$ to measurable quantities (surface tension, viscosity, the distance between the plates, and a line integral over the contour encompassing the growing bubble). The laboratory observations are in good accord with the expression we derive for $d{M}_k/dt$, which is proportional to the surface tension; thus in the zero surface tension limit, the moments (above $k=0$) are all conserved, in accord with Richardson’s theory. In addition, from the measurements of the time evolution of the harmonic moments we obtain a value for the surface tension that is within 20% of the accepted value. In conclusion, our analysis and laboratory observations demonstrate that an interface dynamics description in terms of harmonic moments is physically realizable and robust.

The concept of preferential concentration in the transport of inertial particles by random velocity fields is extended to account for the possibility of zero correlation time and compressibility of the velocity field. It is shown that, in the case of an uncorrelated in time random velocity field, preferential concentration takes the form of a condition on the field history leading to the current particle positions. This generalized form of preferential concentration appears to be a necessary condition for clustering in the uncorrelated in time case. The standard interpretation of preferential concentration is recovered considering local time averages of the velocity field. In the compressible case, preferential concentration occurs in regions of negative divergence of the field. In the incompressible case, it occurs in regions of simultaneously high strain and negative field skewness.

We investigate theoretically the onset of capillary-gravity waves created by a small object moving at the water-air interface. It is well established that, for straight uniform motion, no steady waves appear at velocities below the minimum phase velocity $c_{\text{min}}=23\text{cm}/\text{s}$. At higher velocities, the emission of capillary-gravity waves creates an additional drag force. The behavior of this force near the critical velocity is still poorly understood. A linear-response theory where the object is replaced by an effective pressure source predicts a singular behavior for the wave drag. However, experimental data tend to indicate a more continuous transition. In this paper, we show that a proper treatment of the flow equations around the obstacle can regularize wave emission, even in the linear wave approximation, thereby ensuring a continuous behavior of the drag force.

Taking advantage of a closed-form generalized Maxwell distribution function [P. Asinari and I. V. Karlin, Phys. Rev. E 79, 036703 (2009)] and splitting the relaxation to the equilibrium in two steps, an entropic quasiequilibrium (EQE) kinetic model is proposed for the simulation of low Mach number flows, which enjoys both the $H$ theorem and a free-tunable parameter for controlling the bulk viscosity in such a way as to enhance numerical stability in the incompressible flow limit. Moreover, the proposed model admits a simplification based on a proper expansion in the low Mach number limit (LQE model). The lattice Boltzmann implementation of both the EQE and LQE is as simple as that of the standard lattice Bhatnagar-Gross-Krook (LBGK) method, and practical details are reported. Extensive numerical testing with the lid driven cavity flow in two dimensions is presented in order to verify the enhancement of the stability region. The proposed models achieve the same accuracy as the LBGK method with much rougher meshes, leading to an effective computational speed-up of almost three times for EQE and of more than four times for the LQE. Three-dimensional extension of EQE and LQE is also discussed.

A free vortex in excitable media can be displaced and removed by a wave train. However, simple physical arguments suggest that vortices anchored to large inexcitable obstacles cannot be removed similarly. We show that unpinning of vortices attached to obstacles smaller than the core radius of the free vortex is possible through pacing. The wave-train frequency necessary for unpinning increases with the obstacle size and we present a geometric explanation of this dependence. Our model-independent results suggest that decreasing excitability of the medium can facilitate pacing-induced removal of vortices in cardiac tissue.

Simulation of the deformation of polymers below their glass transition through molecular dynamics provides an useful route to correlate their molecular architecture to deformation behavior. However, present computational capabilities severely restrict the time and length scales that can be simulated when detailed models of these macromolecules are used. Coarse-graining techniques for macromolecular structures intend to make bigger and longer simulations possible by grouping atoms into superatoms and devising ways of determining reasonable force fields for the superatoms in a manner that retains essential macromolecular features relevant to the process under study but jettisons unnecessary details. In this work we systematically develop a coarse-graining scheme aimed at simulating uniaxial stress-strain behavior of polymers below their glass transition. The scheme involves a two step process of obtaining the coarse grained force field parameters above glass transition. This seems to be enough to obtain “faithful” stress-strain responses after quenching to below the glass transition temperature. We apply the scheme developed to a commercially important polymer polystyrene, derive its complete force field parameters and thus demonstrate the effectiveness of the technique.

Exclusion processes on a regular lattice are used to model many biological and physical systems at a discrete level. The average properties of an exclusion process may be described by a continuum model given by a partial differential equation. We combine a general class of contact interactions with an exclusion process. We determine that many different types of contact interactions at the agent-level always give rise to a nonlinear diffusion equation, with a vast variety of diffusion functions $D\left(C\right)$. We find that these functions may be dependent on the chosen lattice and the defined neighborhood of the contact interactions. Mild to moderate contact interaction strength generally results in good agreement between discrete and continuum models, while strong interactions often show discrepancies between the two, particularly when $D\left(C\right)$ takes on negative values. We present a measure to predict the goodness of fit between the discrete and continuous model, and thus the validity of the continuum description of a motile, contact-interacting population of agents. This work has implications for modeling cell motility and interpreting cell motility assays, giving the ability to incorporate biologically realistic cell-cell interactions and develop global measures of discrete microscopic data.

Using the lateral Casimir force—a manifestation of the quantum fluctuations of the electromagnetic field between objects with corrugated surfaces—as the main force transduction mechanism, a nanomechanical device with rich dynamical behaviors is proposed. The device is made of two parallel racks that are moving in the same direction and a pinion in the middle that couples with both racks via the noncontact lateral Casimir force. The built-in frustration in the device causes it to be very sensitive and react dramatically to minute changes in the geometrical parameters and initial conditions of the system. The noncontact nature of the proposed device could help with the ubiquitous wear problem in nanoscale mechanical systems.

We investigate the role of equilibrium methods and stress transfer range in describing the process of damage. We find that equilibrium approaches are not applicable to the description of damage and the catastrophic failure mechanism if the stress transfer is short ranged. In the long-range limit, equilibrium methods apply only if the healing mechanism associated with ruptured elements is instantaneous. Furthermore we find that the nature of the catastrophic failure depends strongly on the stress transfer range. Long-range transfer systems have a failure mechanism that resembles nucleation. In short-range stress transfer systems, the catastrophic failure is a continuous process that, in some respects, resembles a critical point.

A numerical study has been performed by using a combined Euler and Lagrangian method on the convective heat transfer of Cu and ${\text{Al}}_2{\text{O}}_3$ nanofluids under the turbulent flow conditions. The effects of volume fraction of nanoparticles, nanoparticle sizes, and nanoparticle material are investigated. The mechanism of convective heat transfer enhancement in nanofluids has also been investigated, by studying the influence of particle dispersion and two-way interaction between fluid and particle temperature. The results show significant enhancement of heat transfer of nanofluids. The numerical data are compared with the correlation data of the experiments and reasonably good agreement is achieved.

The evolution of traveling and solitary waves in Bose-Einstein condensates (BECs) with a time-dependent scattering length in an attractive/repulsive parabolic potential is studied. The homogeneous balance principle and the $F$-expansion technique are used to solve the one-dimensional Gross-Pitaevskii equation with time-varying coefficients. We obtained three classes of new exact traveling wave and localized solutions. Our results demonstrate that the BEC solitary wave solutions can be manipulated and controlled by the time-dependent scattering length.

We present a simple method to quantitatively image the dielectric permittivity of soft materials at nanoscale using electrostatic force microscopy (EFM) by means of the double pass method. The EFM experiments are based on the measurement of the frequency shifts of the oscillating tip biased at two different voltages. A numerical treatment based on the equivalent charge method allows extracting the values of the dielectric permittivity at each image point. This method can be applied with no restrictions of film thickness and tip radius. This method has been applied to image the morphology and the nanodielectric properties of a model polymer blend of polystyrene and poly(vinyl acetate).

We consider the dynamics of a disordered $p$-spin model with $p=2$, analyzing the dynamics within Ruelle’s thermodynamic formalism, We use an indicator of the dynamical activity to construct the relevant dynamical Gibbs ensemble. We prove that the dynamics in the low-temperature (spin-glass) phase of the model take place at a second-order phase transition between dynamically active and inactive trajectories. We also show that the same behavior is found in a related model of a three-dimensional ferromagnet.

We study large deviations of the dynamical activity in the random orthogonal model. This is a fully connected spin-glass model with one-step replica symmetry-breaking behavior, consistent with the random first-order transition scenario for structural glasses. We show that this model displays dynamical (space-time) phase transitions between active and inactive phases, as demonstrated by singularities in large deviation functions. We argue that such transitions are generic in systems with long-lived metastable states.