The deviations of the glass transition temperature ($T_g$) in thin films of an amorphous conjugated polymer poly(9,9-dioctylfluorene-co-$N$-(4-butylphenyl)diphenylamine) (TFB) are reported. Monotonic and nonmonotonic $T_g$ deviations are observed in TFB thin films supported on Si-SiOx and poly(3,4-ethylenedioxythiophene) poly(styrenesulfonate) (PEDOT:PSS), respectively. A three-layer model is developed to fit both monotonic and nonmonotonic $T_g$ deviations in these films. A 5-nm PEDOT:PSS capping layer was not found to be effective to remove the free-surface effect in Si-SiOx supported TFB films.

The effect of overpasses in the Biham-Middleton-Levine traffic flow model with random and parallel update rules has been studied. An overpass is a site that can be occupied simultaneously by an eastbound car and a northbound one. Under periodic boundary conditions, both self-organized and random patterns are observed in the free-flowing phase of the parallel update model, while only the random pattern is observed in the random update model. We have developed mean-field analysis for the moving phase of the random update model, which agrees with the simulation results well. An intermediate phase is observed in which some cars could pass through the jamming cluster due to the existence of free paths in the random update model. Two intermediate states are observed in the parallel update model, which have been ignored in previous studies. The intermediate phases in which the jamming skeleton is only oriented along the diagonal line in both models have been analyzed, with the analyses agreeing well with the simulation results. With the increase of overpass ratio, the jamming phase and the intermediate phases disappear in succession for both models. Under open boundary conditions, the system exhibits only two phases when the ratio of overpasses is below a threshold in the random update model. When the ratio of the overpass is close to 1, three phases could be observed, similar to the totally asymmetric simple exclusion process model. The dependence of the average velocity, the density, and the flow rate on the injection probability in the moving phase has also been obtained through mean-field analysis. The results of the parallel model under open boundary conditions are similar to that of the random update model.

Nonequilibrium self-organized patterns formed by particles interacting through competing range interaction are driven over a substrate by an external force. We show that, with increasing driving force, the preexisted static patterns evolve into dynamic patterns either via disordered phase or depinned patterns or via the formation of nonequilibrium stripes. Strikingly, the stripes are formed either in the direction of the driving force or in the transverse direction, depending on the pinning strength. The revealed dynamical patterns are summarized in a dynamical phase diagram.

We consider radial sine-Gordon kinks in two, three, and higher dimensions. A full two-dimensional simulation showing that azimuthal perturbations remain small allows us to reduce the problem to the one-dimensional radial sine-Gordon equation. We solve this equation on an interval $[r_0,r_1]$ and absorb all outgoing radiation. As the kink shrinks toward $r_0$, before the collision, its motion is well described by a simple law derived from the conservation of energy. In two dimensions for $r_0\le 2$, the collision disintegrates the kink into a fast breather, while for $r_0\ge 4$ we obtain a kink-breather metastable state where breathers are shed at each kink “return.” In three and higher dimensions $d$, an additional kink-oscillon state appears for small $r_0$. On the application side, the kink disintegration opens the way for new types of terahertz microwave generators.

We present a semiclassical approximation to the scattering wave function $\Psi (\mathbf{r},k)$ for an open quantum billiard, which is based on the reconstruction of the Feynman path integral. We demonstrate its remarkable numerical accuracy for the open rectangular billiard and show that the convergence of the semiclassical wave function to the full quantum state is controlled by the mean path length or equivalently the dwell time for a given scattering state. In the numerical implementation a cutoff length in the maximum path length or, equivalently, a maximum dwell time ${\tau }_{\max }$ included implies a finite energy resolution $\Delta E\sim {\tau }_{\max }^{-1}$. Possible applications include leaky billiards and systems with decoherence present.

The evolution of the interface separating a conduit of light, viscous fluid rising buoyantly through a heavy, more viscous, exterior fluid at small Reynolds numbers is governed by the interplay between nonlinearity and dispersion. Previous authors have proposed an approximate model equation based on physical arguments, but a precise theoretical treatment for this two-fluid system with a free boundary is lacking. Here, a derivation of the interfacial equation via a multiple scales, perturbation technique is presented. Perturbations about a state of vertically uniform, laminar conduit flow are considered in the context of the Navier-Stokes equations with appropriate boundary conditions. The ratio of interior to exterior viscosities is the small parameter used in the asymptotic analysis, which leads systematically to a maximal balance between buoyancy driven, nonlinear self-steepening and viscous, interfacial stress induced, nonlinear dispersion. This results in a scalar, nonlinear partial differential equation describing large amplitude dynamics of the cross-sectional area of the intrusive fluid conduit, in agreement with previous derivations. The leading order behavior of the two-fluid system is completely characterized in terms of the interfacial dynamics. The regime of model validity is characterized and shown to agree with previous experimental studies. Viscous fluid conduits provide a robust setting for the study of nonlinear, dispersive wave phenomena.

We investigate the spontaneous oscillations of drops levitated above an air cushion, eventually inducing a breaking of axisymmetry and the appearance of “star drops”. This is strongly reminiscent of the Leidenfrost stars that are observed for drops floating above a hot substrate. The key advantage of this work is that we inject the airflow at a constant rate below the drop, thus eliminating thermal effects and allowing for a better control of the flow rate. We perform experiments with drops of different viscosities and observe stable states, oscillations, and chimney instabilities. We find that for a given drop size the instability appears above a critical flow rate, where the latter is largest for small drops. All these observations are reproduced by numerical simulations, where we treat the drop using potential flow and the gas as a viscous lubrication layer. Qualitatively, the onset of instability agrees with the experimental results, although the typical flow rates are too large by a factor 10. Our results demonstrate that thermal effects are not important for the formation of star drops and strongly suggest a purely hydrodynamic mechanism for the formation of Leidenfrost stars.

Models of porous media are often applied to relatively small systems, which leads not only to system-size-dependent results, but also to phenomena that would be absent in larger systems. Here we investigate one such finite-size effect: anisotropy of the permeability tensor. We show that a nonzero angle between the external body force and macroscopic flux vector exists in three-dimensional periodic models of sizes commonly used in computer simulations and propose a criterion, based on the ratio of the system size to the grain size, for this phenomenon to be relevant or negligible. The finite-size anisotropy of the porous matrix induces a pressure gradient perpendicular to the axis of a porous duct and we analyze how this effect scales with the system and grain sizes. We also analyze how the size of the representative elementary volume (REV) for anisotropy compares with the REV for permeability.

Space platforms such as the Space Shuttle and International Space Station have been considered an ideal environment for production of protein and semiconductor crystals of superior quality due to the negligible gravity-induced convection. Although it was believed that under microgravity environment diffusive mass transport would dominate the growth of the crystals, some related experiments have not shown satisfactory results possibly due to the movement of the growing crystals in fluid cells caused by small vibrations present in the space platforms called $g$-jitter. In ground-based experiments, there have been clear observations of attraction and repulsion of a solid particle with respect to a nearby wall of the fluid cell due to small vibrations. The present work is a numerical investigation on the physical mechanisms responsible for the repulsion force, which has been predicted to increase with the cell vibration frequency and amplitude, as well as the fluid viscosity. Moreover, the simulations have revealed that the repulsion force occurs mostly due to the increased pressure in the narrow gap between the particle and the nearest wall.

We study the self-modulation of a circularly polarized Alfvén wave in a strongly magnetized relativistic electron-positron plasma with finite temperature. This nonlinear wave corresponds to an exact solution of the equations, with a dispersion relation that has two branches. For a large magnetic field, the Alfvén branch has two different zones, which we call the normal dispersion zone (where $d\omega /dk>0$) and the anomalous dispersion zone (where $d\omega /dk<0$). A nonlinear Schrödinger equation is derived in the normal dispersion zone of the Alfvén wave, where the wave envelope can evolve as a periodic wave train or as a solitary wave, depending on the initial condition. The maximum growth rate of the modulational instability decreases as the temperature is increased. We also study the Alfvén wave propagation in the anomalous dispersion zone, where a nonlinear wave equation is obtained. However, in this zone the wave envelope can evolve only as a periodic wave train.

We discuss Lehmann effects and rotato-electricity for liquid crystalline phases made of achiral molecules. We point out that for static and dynamic Lehmann effects to exist, it is not necessary to have chiral molecules provided the overall structure has macroscopic chirality. This question is of direct relevance for liquid crystalline phases formed by bent-core molecules provided they have a sufficiently low symmetry. This includes systems which break parity symmetry and have overall $C_2$ or $C_1$ symmetry. We point out that for liquid crystalline gels and elastomers one should be able to observe rotato-electricity for systems with macroscopic chirality. Rotatoelectricity is associated with the relative rotations between two subsystems, namely, between the network and the director, in an external electric field. Candidates include gels and even monolayers prepared from bent-core molecules with sufficiently low symmetry.

Recently, a critical test of the Navier-Stokes-Fourier equations for compressible fluid continua was proposed [H. Brenner, Phys. Rev. E 87, 013014 (2013)]. It was shown that the equations of bivelocity hydrodynamics imply that a compressible fluid in an isolated rotating circular cylinder attains a nonequilibrium steady state with a nonuniform temperature increasing radially with distance from the axis. We demonstrate that statistical mechanical arguments, involving Hamiltonian dynamics and ergodicity due to irregularity of the wall, lead instead to a thermal equilibrium state with uniform temperature. This is the situation to be expected in experiment.

We discuss a method of organizing incoherent motion of a colloidal suspension to produce synchronized, coherent motion, extending the discussion of our recent Letter [Moths and Witten, Phys. Rev. Lett. 110, 028301 (2013)]. The method does not require interaction between the objects. Instead, the effect is controlled by the “twist matrix” which gives the angular velocity of an asymmetric object in a fluid resulting from a weak external force. We analyze the two types of forcing considered in the Letter: a force alternating between two directions and a continuously rotating force. For the alternating force, we justify the claim of the Letter that under appropriate forcing conditions, the orientational entropy of the objects decreases indefinitely with time, on average. We provide a bound on that rate in terms of the twist matrix. For the case of rotating force, we derive conditions for phased-locked motion of the objects to the force and prove that there is only one stable phase-locked orientation under these conditions. We find numerically that the fastest alignment typically occurs for tilt angles of order unity. We discuss how the alignment effect scales with the object size for external forcing caused by gravity or an electric field. Under practical forcing conditions we estimate that the alignment should persist despite rotational diffusion for objects larger than about 10 microns. Potential misalignment owing to hydrodynamic interaction of the objects is estimated to be negligible at volume fractions smaller than about ${10}^{-4.5}$ (${10}^{-3}$) when the forcing is gravitational (electrophoretic).

In the context of the random first order transition theory we use an extended mode coupling theory of the glass transition that includes activated events to account for spatiotemporal structures in rejuvenating glasses. We numerically solve fluctuating dynamical equations for mobility and fictive temperature fields which capture both mobility generation through activated events and facilitation effects. Upon rejuvenating, a source of high mobility at a glass surface initiates a growth front of mobility which propagates into the unstable low mobility region. The speed of the front quantitatively agrees with experiments on the rejuvenation of ultrastable glasses, which “melt” from their surface.

Dynamics of the sub-Ohmic spin-boson model is examined using three numerical approaches, namely the Dirac-Frenkel time-dependent variation with the Davydov D${}_1$ ansatz, the adaptive time-dependent density matrix renormalization group method within the representation of orthogonal polynomials, and a perturbative approach based on a unitary transformation. In order to probe the validity regimes of the three approaches, we study the dynamics of a qubit coupled to a bosonic bath with and without a local field. Comparison of the up-state population evolution shows that the three approaches are in agreement in the weak-coupling regime but exhibit marked differences when the coupling strength is large. The Davydov D${}_1$ ansatz and the time-dependent density matrix renormalization group can both be reliably employed in the weak-coupling regime, while the former is also valid in the strong-coupling regime as judged by how faithfully the trial state follows the Schrödinger equation. We further explore the bipartite entanglement dynamics between two qubits coupled with individual bosonic baths which reveals entanglement sudden death and revival.

There exists a plethora of mathematical models for cluster growth and/or aggregation on regular lattices. Almost all suffer from inherent anisotropy caused by the regular lattice upon which they are grown. We analyze the little-known model for stochastic cluster growth on a regular lattice first introduced by Ferreira Jr. and Alves [J. Stat. Mech. Theo. & Exp. (2006) P11007], which produces circular clusters with no discernible anisotropy. We demonstrate that even in the noise-reduced limit the clusters remain circular. We adapt the model by introducing a specific rearrangement algorithm so that, rather than adding elements to the cluster from the outside (corresponding to apical growth), our model uses mitosis-like cell splitting events to increase the cluster size. We analyze the surface scaling properties of our model and compare it to the behavior of more traditional models. In “$1+1$” dimensions we discover and explore a new, nonmonotonic surface thickness scaling relationship which differs significantly from the Family-Vicsek scaling relationship. This suggests that, for models whose clusters do not grow through particle additions which are solely dependent on surface considerations, the traditional classification into “universality classes” may not be appropriate.

We consider grain growth and stagnation in polycrystalline microstructures. From the phase-field crystal modeling of the coarsening dynamics, we identify a transition from a grain-growth stagnation upon deep quenching below the melting temperature $T_m$ to a continuous coarsening at shallower quenching near $T_m$. The grain evolution is mediated by local grain rotations. In the deep quenching regime, the grain assembly typically reaches a metastable state where the kinetic barrier for recrystallization across boundaries is too large and grain rotation with subsequent coalescence or boundary motion is infeasible. For quenching near $T_m$, we find that the grain growth depends on the average rate of grain rotation, and follows a power-law behavior with time, with a scaling exponent that depends on the quenching depth.

We consider an ensemble of fully connected networks of $N$ oscillators coupled harmonically with random springs and show, using random-matrix-theory considerations, that both the average phonon heat current and its variance are scale invariant and take universal values in the large $N$ limit. These anomalous mesoscopic fluctuations is the hallmark of strong correlations between normal modes.

The mean-field dynamics of an electric dipole moment in a deterministic and a fluctuating electric field is solved to obtain the average over fluctuations of the dipole moment and the angular momentum as a function of time for a Gaussian white-noise stochastic electric field. The components of the average electric dipole moment and the average angular momentum along the deterministic electric-field direction do not decay to zero, despite fluctuations in all three components of the electric field. This is in contrast to the decay of the average over fluctuations of a magnetic moment in a stochastic magnetic field with Gaussian white noise in all three components. The components of the average electric dipole moment and the average angular momentum perpendicular to the deterministic electric-field direction oscillate with time but decay to zero, and their variance grows with time.

We present a theory of effective electrostatic interactions in polydisperse suspensions of charged macroions, generalizing to mixtures a theory previously developed for monodisperse suspensions. Combining linear response theory with a random phase approximation for microion correlations, we coarse grain the microion degrees of freedom to derive general expressions for effective macroion-macroion pair potentials and a one-body volume energy. For model mixtures of charged hard-sphere colloids, we give explicit analytical expressions. The resulting effective pair potentials have the same general form as predicted by linearized Poisson-Boltzmann theory, but consistently incorporate dependence on macroion density and excluded volume via the Debye screening constant. The volume energy, which depends on the average macroion density, contributes to the free energy and so can influence thermodynamic properties of deionized suspensions. To validate the theory, we compute radial distribution functions of binary mixtures of oppositely charged colloidal macroions from molecular dynamics simulations of the coarse-grained model (with implicit microions), taking effective pair potentials as input. Our results agree closely with corresponding results from more computationally intensive Monte Carlo simulations of the primitive model (with explicit microions). Simulations of a mixture with large size and charge asymmetries indicate that charged nanoparticles can enhance electrostatic screening of charged colloids. The theory presented here lays a foundation for future large-scale modeling of complex mixtures of charged colloids, nanoparticles, and polyelectrolytes.