We report the first measurement of the hindrance function for sedimentation or creaming of disk-shaped colloids via the analytical centrifugation. Disks align with the external flow right above a volume fraction of a few percent, and this effect is extremely sensitive to the aspect ratio of disks. Due to this alignment effect, disk sedimentation or creaming demonstrates distinct trends in dilute and semidilute regions.

Spherical Couette flow involves fluid sheared between concentric coaxially rotating spheres. Its scientific relevance lies not only in the simplicity of the system but also in its applicability to astrophysical objects such as atmospheres, oceans, and planetary cores. One common behavior in all rotating flows, including spherical Couette flow, is the presence of inertial modes, which are linear wave modes restored by the Coriolis force. Building on a previous identification of inertial modes in a laboratory spherical Couette cell, here we propose selection mechanisms to explain the presence of the particular modes we have observed. Mode selection depends on both amplification and damping. Our experimental observations are consistent with amplification and selection by over-reflection at a shear layer, and we would expect other spherical Couette devices to behave similarly. Damping effects, due in part to the presence of an inner sphere, add further constraints which are likely to play a role in mode selection in planetary atmospheres and cores, including the core of earth.

The velocity and stress autocorrelation decay in a dissipative particle dynamics ideal fluid model is analyzed in this paper. The autocorrelation functions are calculated at three different friction parameters and three different time steps using the well-known Groot/Warren algorithm and newer algorithms including self-consistent leap-frog, self-consistent velocity Verlet and Shardlow first and second order integrators. At low friction values, the velocity autocorrelation function decays exponentially at short times, shows slower-than exponential decay at intermediate times, and approaches zero at long times for all five integrators. As friction value increases, the deviation from exponential behavior occurs earlier and is more pronounced. At small time steps, all the integrators give identical decay profiles. As time step increases, there are qualitative and quantitative differences between the integrators. The stress correlation behavior is markedly different for the algorithms. The self-consistent velocity Verlet and the Shardlow algorithms show very similar stress autocorrelation decay with change in friction parameter, whereas the Groot/Warren and leap-frog schemes show variations at higher friction factors. Diffusion coefficients and shear viscosities are calculated using Green-Kubo integration of the velocity and stress autocorrelation functions. The diffusion coefficients match well-known theoretical results at low friction limits. Although the stress autocorrelation function is different for each integrator, fluctuates rapidly, and gives poor statistics for most of the cases, the calculated shear viscosities still fall within range of theoretical predictions and nonequilibrium studies.

We present a method for calculation of the potential and related physical quantities experienced by a particle traversing an aligned periodic complex atomic structure. Classical physics equations and the expansion of periodic functions as a Fourier series have been used for the calculation. Based on this method, we have developed the ECHARM program, which calculates one- and two-dimensional averaged physical quantities of interest along the main axes of any orthorhombic and tetragonal structure. For the case of cubic symmetry, the calculation holds for any orientation. Complex structures such as zeolites have been worked out to show the capability of the program.

Spectral properties of the weighted Laplace operator in the presence of fractal boundaries are numerically investigated for both Neumann and Dirichlet boundary conditions. This corresponds to the characterization of heat and mass transport in microchannels with irregular and rough surfaces induced by the microfabrication process. The axial velocity field with no-slip boundary conditions, representing the weighting function of the Laplace operator, influences the localization properties of the eigenfunctions and the scaling of the integrated density of state (IDOS) $N\left(\lambda \right)$. The results indicate that $N\left(\lambda \right)$ deviates from the form given by the modified Weyl-Berry-Lapidus conjecture as it shows a correction of $\Delta N\left(\lambda \right)\sim {\lambda }^{D_f/4}$ to the leading-order Weil term. Numerical results are presented for Koch and Koch snowflake fractal boundaries. The role of slip or no-slip boundary conditions of the velocity field on the IDOS is also investigated.

Multiarm star polymers are model systems with tunable intermediate colloid to polymerlike character, exhibiting rich phase behavior, internal relaxations, and flow properties. An important puzzle for several years has been the lack of clear experimental proof of crystalline states despite strong theoretical predictions. We present unambiguous evidence via multispeckle dynamic light scattering (MSDLS) and small-angle neutron scattering (SANS) for such crystallization in a solvent of intermediate quality. An unexpected speed up of the short-time star diffusion observed in MSDLS was attributed by SANS to crystallization, via aging, of the multiam star glass. This delayed glass to crystal transition establishes a pathway for star crystallization that might be generic in colloidal glasses.

In this paper we extend the concept of persistence, well defined for classical stochastic dynamics, to the context of quantum dynamics. We demonstrate the idea via quantum random walk and a successive measurement scheme, where persistence is defined as the time during which a given site remains unvisited by the walker. We also investigated the behavior of related quantities, e.g., the first-passage time and the succession probability (newly defined), etc. The study reveals power-law scaling behavior of these quantities with new exponents. Comparable features of the classical and the quantum walks are discussed.

The volume fluctuations in the steady state reached by a vibrated granular gas of hard particles confined by a movable piston on the top are investigated by means of event-driven simulations. Also, a compressibility factor, measuring the response in volume of the system to a change in the mass of the piston, is introduced and measured. From the second moment of the volume fluctuations and the compressibility factor, an effective temperature is defined by using the same relation as obeyed by equilibrium molecular systems. The interpretation of this effective temperature and its relationship with the granular temperature of the gas, and also with the velocity fluctuations of the movable piston, is discussed. It is found that the ratio of the temperature based on the volume fluctuations to the temperature based on the piston kinetic energy obeys simple dependencies on the inelasticity and on the piston-particle mass ratio.

We apply a formalism for deriving stochastic continuum equations associated with lattice models to obtain equations governing the transient regimes of epitaxial growth for various experimental scenarios and growth conditions. The first step of our methodology is the systematic transformation of the lattice model into a regularized stochastic equation of motion that provides initial conditions for differential renormalization-group (RG) equations for the coefficients in the regularized equation. The solutions of the RG equations then yield trajectories that describe the original model from the transient regimes, which are of primary experimental interest, to the eventual crossover to the asymptotically stable fixed point. We first consider regimes defined by the relative magnitude of deposition noise and diffusion noise. If the diffusion noise dominates, then the early stages of growth are described by the Mullins-Herring (MH) equation with conservative noise. This is the classic regime of molecular-beam epitaxy. If the diffusion and deposition noise are of comparable magnitude, the transient equation is the MH equation with nonconservative noise. This behavior has been observed in a recent report on the growth of aluminum on silicone oil surfaces [Z.-N. Fang et al., Thin Solid Films 517, 3408 (2009)]. Finally, the regime where deposition noise dominates over diffusion noise has been observed in computer simulations, but does not appear to have any direct experimental relevance. For initial conditions that consist of a flat surface, the Villain-Lai-Das Sarma (VLDS) equation with nonconservative noise is not appropriate for any transient regime. If, however, the initial surface is corrugated, the relative magnitudes of terms can be altered to the point where the VLDS equation with conservative noise does indeed describe transient growth. This is consistent with the experimental analysis of growth on patterned surfaces [H.-C. Kan et al., Phys. Rev. Lett. 92, 146101 (2004); T. Tadayyon-Eslami et al., Phys. Rev. Lett. 97, 126101 (2006)].

We investigated the insertion of transmembrane structures in a lipid bilayer and their interactions using self-consistent field theory. The lipids are coarse-grained on a united-atom level and consist of a phosphatidylcholinelike headgroup and two hydrophobic tails. The inclusions, acting as simple models for proteins that span biological membranes, are rigid rods (radius $R$) with a hydrophobic surface and hydrophilic end caps. The insertion free energy $\Omega$ of an individual rod is strongly regulated by the affinity between its hydrophobic surface and the lipid tails. This affinity also controls the best match of the hydrophobic length of the rod with that of the bilayer. The line tension $\tau \left(=\Omega /2\pi R\right)$ is practically independent of $R$. The perturbations in the bilayer as a function of distance from the inclusion, have the shape of a damped oscillation. The wavelength and decay length are related to the elastic properties of the bilayer and do not depend on $R$. These results are used to analyze how the lipid matrix affects the interaction between transmembrane objects, for computational reasons considering the limit of $R\to \infty$. Contributions on different length scales can be distinguished: (i) a long-range elastic interaction, which is an exponentially decaying oscillation; (ii) an exponentially decaying repulsion on an intermediate length scale, resulting from the loss of conformational entropy of the lipid tails; and (iii) a short-range interaction due to the finite compressibility of the lipid tails, which manifests either as a depletion attraction if there is no affinity between the tails and the inclusions’ surface or, otherwise, as an oscillatory structural force.

We present an experimental procedure to perform broadband microrheological measurements with optical tweezers. A generalized Langevin equation is adopted to relate the time-dependent trajectory of a particle in an imposed flow to the frequency-dependent moduli of the complex fluid. This procedure allows us to measure the material linear viscoelastic properties across the widest frequency range achievable with optical tweezers.

The effect of gas compression on spontaneous flame acceleration leading to deflagration-to-detonation transition is studied theoretically for small Reynolds number flame propagation from the closed end of a tube. The theory assumes weak compressibility through expansion in small Mach number. Results show that the flame front accelerates exponentially during the initial stage of propagation when the Mach number is negligible. With continuous increase in the flame velocity with respect to the tube wall, the flame-generated compression waves subsequently moderate the acceleration process by affecting the flame shape and velocity, as well as the flow driven by the flame.

The dynamic structure factor, which determines the Thomson scattering spectrum, is calculated via an extended Mermin approach. It incorporates the dynamical collision frequency as well as the local-field correction factor. This allows to study systematically the impact of electron-ion collisions as well as electron-electron correlations due to degeneracy and short-range interaction on the characteristics of the Thomson scattering signal. As such, the plasmon dispersion and damping width is calculated for a two-component plasma, where the electron subsystem is completely degenerate. Strong deviations of the plasmon resonance position due to the electron-electron correlations are observed at increasing Brueckner parameters $r_s$. These results are of paramount importance for the interpretation of collective Thomson scattering spectra, as the determination of the free electron density from the plasmon resonance position requires a precise theory of the plasmon dispersion. Implications due to different approximations for the electron-electron correlation, i.e., different forms of the one-component local-field correction, are discussed.

Here we extend the plasmonic cloaking technique to irregularly shaped objects with anisotropic scattering response. The scattering-cancellation approach to cloaking [A. Alù and N. Engheta, Phys. Rev. E 72, 016623 (2005)] has been extensively applied in the past to symmetrical geometries and canonical shapes. However, recent papers have raised some doubts concerning the fact that its use may not be as effective when dealing with strongly anisotropic and noncanonical geometries. Our goal here is to extend the plasmonic cloaking technique to irregular obstacles and to show that proper cloak design may provide a significant and uniform scattering reduction, independent of angle of incidence, position, and polarization of the illumination. We investigate how the volumetric effect of scattering cancellation provided by plasmonic media may drastically suppress the scattering for these irregular geometries independent of the illumination angle, and we shed some light on the physical mechanisms and the design rules at the basis of this cloaking technique when applied to objects whose scattering properties are dependent upon polarization and angle of incidence.

Stress-strain measurements and ultrasound propagation experiments in glass bead packs have been simultaneously conducted to characterize the stress-induced anisotropy under uniaxial loading. These measurements realized, respectively, with finite and incremental deformations of the granular assembly, are analyzed within the framework of the effective-medium theory based on the Hertz-Mindlin contact theory. Our work shows that both compressional and shear wave velocities and consequently the incremental elastic moduli agree fairly well with an effective-medium model developed by Johnson et al. [J. Appl. Mech. 65, 380 (1998)] for the oedometric test, but the anisotropic stress ratio resulting from finite deformation does not at all. As indicated by numerical simulations, the discrepancy may arise from the fact that the model does not properly allow the grains to relax from the affine motion approximation. Here we find that the interaction nature at the grain contact could also play a crucial role for the relevant prediction by the model; indeed, such discrepancy can be significantly reduced if the frictional resistance between grains is removed. Another main experimental finding is the influence of the inherent anisotropy of granular packs, realized by different protocols of the sample preparation. Our results reveal that compressional waves are more sensitive to the stress-induced anisotropy, whereas the shear waves are more sensitive to the fabric anisotropy.

A hierarchical divisive algorithm is proposed for identifying communities in complex networks. To that effect, the definition of community in the weak sense of Radicchi et al. [Proc. Natl. Acad. Sci. U.S.A. 101, 2658 (2004)] is extended into a criterion for a bipartition to be optimal: one seeks to maximize the minimum for both classes of the bipartition of the ratio of inner edges to cut edges. A mathematical program is used within a dichotomous search to do this in an optimal way for each bipartition. This includes an exact solution of the problem of detecting indivisible communities. The resulting hierarchical divisive algorithm is compared with exact modularity maximization on both artificial and real world data sets. For two problems of the former kind optimal solutions are found; for five problems of the latter kind the edge ratio algorithm always appears to be competitive. Moreover, it provides additional information in several cases, notably through the use of the dendrogram summarizing the resolution. Finally, both algorithms are compared on reduced versions of the data sets of Girvan and Newman [Proc. Natl. Acad. Sci. U.S.A. 99, 7821 (2002)] and of Lancichinetti et al. [Phys. Rev. E 78, 046110 (2008)]. Results for these instances appear to be comparable.

We consider a particular type of parametrically driven discrete Klein-Gordon system describing microdevices and nanodevices, with integrated electrical and mechanical functionality. Using a multiscale expansion method we reduce the system to a discrete nonlinear Schrödinger equation. Analytical and numerical calculations are performed to determine the existence and stability of fundamental bright and dark discrete solitons admitted by the Klein-Gordon system through the discrete Schrödinger equation. We show that a parametric driving can not only destabilize onsite bright solitons, but also stabilize intersite bright discrete solitons and onsite and intersite dark solitons. Most importantly, we show that there is a range of values of the driving coefficient for which dark solitons are stable, for any value of the coupling constant, i.e., oscillatory instabilities are totally suppressed. Stability windows of all the fundamental solitons are presented and approximations to the onset of instability are derived using perturbation theory, with accompanying numerical results. Numerical integrations of the Klein-Gordon equation are performed, confirming the relevance of our analysis.

A method is presented for generating compact fractal disordered media by generalizing the random midpoint displacement algorithm. The obtained structures are invasive stochastic fractals, with the Hurst exponent varying as a continuous parameter, as opposed to lacunar deterministic fractals, such as the Menger sponge. By employing the detrending moving average algorithm [A. Carbone, Phys. Rev. E 76, 056703 (2007)], the Hurst exponent of the generated structure can be subsequently checked. The fractality of such a structure is referred to a property defined over a three-dimensional topology rather than to the topology itself. Consequently, in this framework, the Hurst exponent should be intended as an estimator of compactness rather than of roughness. Applications can be envisaged for simulating and quantifying complex systems characterized by self-similar heterogeneity across space. For example, exploitation areas range from the design and control of multifunctional self-assembled artificial nanostructures and microstructures to the analysis and modeling of complex pattern formation in biology, environmental sciences, geomorphological sciences, etc.

On the basis of Monte Carlo simulations, the present work systematically investigates how conductivity exponents depend on the ratio of stick-stick junction resistance to stick resistance for two-dimensional stick percolation. Simulation results suggest that the critical conductivity exponent extracted from size-dependent conductivities of systems exactly at the percolation threshold is independent of the resistance ratio and has a constant value of $1.280\pm 0.014$. In contrast, the apparent conductivity exponent extracted from density-dependent conductivities of systems well above the percolation threshold monotonically varies with the resistance ratio, following an error function, and lies in the vicinity of the critical exponent.