We present an analytical treatment of the dissipative-stochastic dynamics of a charged classical particle confined biharmonically in a plane with a uniform static magnetic field directed perpendicular to the plane. The stochastic dynamics gives a steady state in the long-time limit. We have examined the orbital magnetic effect of introducing a parametrized deviation $(\eta -1)$ from the second fluctuation-dissipation relation that connects the driving noise and the frictional memory kernel in the standard Langevin dynamics. The main result obtained here is that the moving charged particle generates a finite orbital magnetic moment in the steady state, and that the moment shows a crossover from para- to diamagnetic sign as the parameter $\eta$ is varied. It is zero for $\eta =1$ that makes the steady state correspond to equilibrium, as it should. The magnitude of the orbital magnetic moment turns out to be a nonmonotonic function of the applied magnetic field, tending to zero in the limit of an infinitely large as well as an infinitesimally small magnetic field. These results are discussed in the context of the classic Bohr–van Leeuwen theorem on the absence of classical orbital diamagnetism. Possible realization is also briefly discussed.

A renormalization approach is used to derive an analytic formula for the total current crossing the reactive surface of a Cayley tree of cylindrical tubes under a Helmholtz-type approximation to the full diffusion-reaction problem. We provide analytic conditions for the emergence of a plateau in the current—a region of maximum insensitivity of the current to variations in either the reaction rate (membrane permeability) or the diffusivity. The occurrence of such a plateau is associated with a partial screening regime wherein most of the active surface is screened to incoming diffusing particles. Large trees trade efficiency for fault tolerance, a valuable feature which may provide robustness to mammalian respiratory systems and tolerance to catalytic poisoning in chemical reactors.

The Boltzmann equation for inelastic Maxwell models is considered in order to investigate the dynamics of an impurity (or intruder) immersed in a granular gas driven by a uniform shear flow. The analysis is based on an exact solution of the Boltzmann equation for a granular binary mixture. It applies for conditions arbitrarily far from equilibrium (arbitrary values of the shear rate $a$) and for arbitrary values of the parameters of the mixture (particle masses $m_i$, mole fractions $x_i$, and coefficients of restitution ${\alpha }_{ij}$). In the tracer limit where the mole fraction of the intruder species vanishes, a nonequilibrium phase transition takes place. We thereby identify ordered phases where the intruder bears a finite contribution to the properties of the mixture, in a region of parameter space that is worked out in detail. These findings extend previous results obtained for ordinary Maxwell gases, and further show that dissipation leads to new ordered phases.

We explore the topology of configuration spaces of hard disks experimentally and show that several changes in the topology can already be observed with a small number of particles. The results illustrate a theorem of Baryshnikov, Bubenik, and Kahle [2] that critical points correspond to configurations of disks with balanced mechanical stresses and suggest conjectures about the asymptotic topology as the number of disks tends to infinity.

Off-equilibrium dynamics of a three-dimensional lattice model with nearest- and next-nearest-neighbors exclusions is studied. At equilibrium, the model undergoes a first-order fluid-solid transition. Nonequilibrium filling, through random sequential adsorption with diffusion, creates amorphous structures and terminates at a disordered state with random closest packing density that lies in the equilibrium solid regime. The approach toward random closest packing is characterized by two distinct power-law regimes, reflecting the formation of small densely packed grains in the long-time regime of the filling process. We then study the fixed-density relaxation of these amorphous structures toward the solid phase. The route to crystallization for the high-density, deeply supercooled regime is shown to deviate from the simple grain growth proposed by classical nucleation theory. Our measurements suggest that relaxation in this regime is driven mainly by coalescence of neighboring crystallized grains which exist in the initial amorphous state.

A microscale shear cell is used to study the formation of parabolic focal conic defects in the thermotropic smectic-$A$ liquid crystal 8CB (4-octyl-4′-cyanobiphenyl). Defects are produced by four distinct methods: by the application of dilatational strain alone, by shear flow alone, by dilatational strain and subsequent shear flow, and by the simultaneous application of dilatational strain and shear flow. We confirm that defects originate within the bulk, consistent with the previously suggested undulation instability mechanism. In the presence of a shear flow, we observe that defect formation requires micrometer-level dilatations, whose magnitude depends on the sample thickness. The size and ordering of both disordered and ordered defect arrays is quantified using a pair distribution function. Deviations from the predictions of linear stability theory are observed that have not been reported previously. For example, defects form a square array with greater ordering in the principal flow direction. Ordering due to shear flow does not change the average defect size. It has been shown previously that the principal defect sizes of ordered defects scale differently with sample thickness than the wavelength of the small amplitude undulations. We find that disordered defects show a similar deviation from this predicted wavelength.

During the last decade, the development of nanomaterials to penetrate inside living cells has been the focus of a large number of studies, with applications for the biomedical field. However, the further dynamics of these nanomaterials inside the cells is dictated by the intracellular environment and in particular its mechanical properties. The mechanical characteristics of the cell interior can be probed with either active or passive microrheological approaches. However, active intracellular microrheology is still in its infancy, owing to the difficulty of inserting probes that can be manipulated by external forces. Here we review recent active microrheology studies using magnetic nanoprobes inserted into endosomes or phagosomes as useful approaches for measuring frequency-dependent viscoelasticity, for mapping the viscoelastic landscape, as well as for identifying the contribution of individual cytoskeleton components and the influence of cell motility. The results of such direct measurements challenge the validity of more typical passive approaches in which the spontaneous displacement of embedded nanoprobes is measured. Here we discuss that one must distinguish probes suitable for use in conditions of thermal equilibrium, whose movements reflect the mechanical environment from probes that interact actively with the cytoplasm and cytoskeleton, in a state of nonequilibrium for which fluctuation-dissipation theorem no longer holds. However, when data on these probes’ viscoelastic microenvironment is available, such passive probe movements can yield useful information on the forces responsible for intracellular activity.

For two-phase flow in porous media, the natural medium heterogeneity necessarily gives rise to capillary nonequilibrium effects. The relaxation to the equilibrium is a slow process which should be introduced in macroscopic flow models. Many nonequilibrium models are based on a phenomenological approach. At the same time there exists a rigorous mathematical way to develop the nonequilibrium equations. Its formalism, developed by Bourgeat and Panfilov [Computational Geosciences 2, 191 (1998)], is based on the homogenization of the microscale flow equations over medium heterogeneities. In contrast with the mentioned paper, in which the case of a sufficiently fast relaxation was analyzed, we consider the case of long relaxation, which leads to the appearance of long-term memory on the macroscale. Due to coupling between the nonlinearity and nonlocality in time, the macroscopic model remains, however, incompletely homogenized, in the general case. At the same time, frequently only the relationship for the nonequilibrium capillary pressure is of interest for applications. In the present paper, we obtain such an exact relationship in two different independent forms for the case of long-term memory. This relationship is more general than that obtained by Bourgeat and Panfilov. In addition, we prove the comparison theorem which determines the upper and lower bounds for the macroscopic model. These bounds represent linear flow models, which are completely homogenized. The results obtained are illustrated by numerical simulations.

For cavitation clouds generated in a standing sound wave driven by an ultrasonic horn, the nonlinear acoustic wave equation governing cavitation dynamics is numerically solved together with the bubble motion equation under an approximation. This conceptual calculation can qualitatively reproduce the observed characteristics of cavitation.

A three-dimensional instability analysis of an electrified non-Newtonian liquid jet is performed in this paper so as to understand the competition between viscoelastic stresses, electric force, and surface tension. The analysis employs the leaky-dielectric model to account for the charge transport process and the Oldroyd-B model to take the viscoelastic effect into account. Results show that the viscoelastic stresses play a stabilizing role while electrification destabilizes the disturbances with shorter wavelength and higher azimuthal wave numbers. However, the viscoelastic stabilizing effect is not so significant on the sinuous mode (with azimuthal wave number $m$ $=$ 1) compared with the other modes, which thus promotes sinuous disturbances. This result is consistent with the presence of bending motion in most experimental observations.

Cylindrical nonlinear optics is a burgeoning research area which describes cylindrical electromagnetic wave propagation in nonlinear media. Finding new exact solutions for different types of nonlinearity and inhomogeneity to describe cylindrical electromagnetic wave propagation is of great interest and meaningful for theory and application. This paper gives exact solutions for the cylindrical nonlinear Maxwell equations and presents an interesting connection between the exact solutions for different cylindrical nonlinear Maxwell equations. We also provide some examples and discussion to show the application of the results we obtained. Our results provide the basis for solving complex systems of nonlinearity and inhomogeneity with simple systems.

The time-dependent Schrödinger equation is solved by using an explicitly time-dependent basis. This approach allows efficient reflection-free time propagation of the wave function. The applicability of the method is illustrated by solving various time-dependent problems including the calculation of the above threshold ionization of a model atom and the optical absorption spectrum of a sodium dimer.

We use the method of dimensional continuation to isolate singularities in integrals containing products of Green's functions or their derivatives. Rules for the extraction of the finite part of so-called hypersingular integrals are developed, which should be useful in methods based on boundary integral techniques in science and engineering. In applications to potential theory, electromagnetic scattering, and crack dynamics in continuum mechanics, boundary integrals now can be readily evaluated using computational techniques without recourse to complex analysis or contour distortions since the hypersingularities occurring in intermediate steps of the computations can be isolated and ignored while taking the finite parts of the integrals into account in a consistent manner. We have also identified new forms of the Dirac $\delta$ function in $D$ dimensions, which are useful and convenient in the calculations. A summary of the integrable singular integrals is given in tabular form. We extend the considerations to a wider class of Green's functions and present a theorem, with additional results arising from it, that shows that hypersingular integrals associated with three-dimensional potential problems can be reduced to one-dimensional finite integrals rather than two-dimensional integrals, again leading to direct evaluations in such cases. These calculations are compared with existing results to show the efficacy of the approach.

Recently there has been interest in developing efficient ways to model heterogeneous surface reactions with hybrid computational models that couple a kinetic Monte Carlo (KMC) model for a surface to a finite-difference model for bulk diffusion in a continuous domain. We consider two representative problems that validate a hybrid method and show that this method captures the combined effects of nonlinearity and stochasticity. We first validate a simple deposition-dissolution model with a linear rate showing that the KMC-continuum hybrid agrees with both a fully deterministic model and its analytical solution. We then study a deposition-dissolution model including competitive adsorption, which leads to a nonlinear rate, and show that in this case the KMC-continuum hybrid and fully deterministic simulations do not agree. However, we are able to identify the difference as a natural result of the stochasticity coming from the KMC surface process. Because KMC captures inherent fluctuations, we consider it to be more realistic than a purely deterministic model. Therefore, we consider the KMC-continuum hybrid to be more representative of a real system.

We study the site-bond percolation on a hierarchical scale-free network, namely, the decorated (2,2)-flower, by using the renormalization group technique. The phase diagram essentially depends on the fraction of occupied sites. Surprisingly, when each site is unoccupied even with a small probability, the system permits neither the percolating phase nor the nonpercolating phase, but rather only critical phases. Although the order parameter always remains zero, a transition still exists between the critical phases that is characterized by the value of the fractal exponent, which measures the degree of criticality; the system changes from one critical state to another with the jump of the fractal exponent at the transition point. The phase boundary depends on the fraction of occupied sites. When the fraction of unoccupied sites exceeds a certain value, the transition line between the critical phases disappears, and a unique critical phase remains.

In this paper we present a first-principles analysis of the nonequilibrium work distribution and the free energy difference of a quantum system interacting with a general environment (with arbitrary spectral density and for all temperatures) based on a well-understood microphysics (quantum Brownian motion) model under the conditions stipulated by the Jarzynski equality [Jarzynski, Phys. Rev. Lett. 78, 2690 (1997)] and Crooks’ fluctuation theorem [Crooks, Phys. Rev. E 60, 2721 (1999)] (in short, fluctuation theorems, FTs). We use the decoherent histories conceptual framework to explain how the notion of trajectories in a quantum system can be made viable and use the environment-induced decoherence scheme to assess the strength of noise that could provide sufficient decoherence to warrant the use of trajectories to define work in open quantum systems. From the solutions to the Langevin equation governing the stochastic dynamics of such systems we were able to produce formal expressions for these quantities entering in the FTs and from them prove explicitly the validity of the FTs at the high temperature limit. At low temperatures our general results would enable one to identify the range of parameters where FTs may not hold or need be expressed differently. We explain the relation between classical and quantum FTs and the advantage of this microphysics open-system approach over the phenomenological modeling and energy-level calculations for substitute closed quantum systems.

We employ a mean-field approximation to study the Ising model with aperiodic modulation of its interactions in one spatial direction. Two different values for the exchange constant, $J_A$ and $J_B$, are present, according to the Fibonacci sequence. We calculate the pseudocritical temperatures for finite systems and extrapolate them to the thermodynamic limit. We explicitly obtain the exponents $\beta$, $\delta$, and $\gamma$ and, from the usual scaling relations for anisotropic models at the upper critical dimension (assumed to be 4 for the model we treat), we calculate $\alpha$, $\nu$, ${\nu }_{\parallel }$, $\eta$, and ${\eta }_{\parallel }$. Within the framework of a renormalization-group approach, the Fibonacci sequence is a marginal one and we obtain exponents that depend on the ratio $r\equiv J_B/J_A$, as expected; however, the scaling relation $\gamma =\beta (\delta -1)$ is obeyed for all values of $r$ we studied. We characterize some thermodynamic functions as log-periodic functions of their arguments, as expected for aperiodic-modulated models, and obtain precise values for the exponents from this characterization.

We have simulated the relative shear motion of both neutral and polyelectrolyte end-grafted polymer brushes using molecular dynamics. The flexible neutral polymer brush is treated as a bead-spring model, and the polyelectrolyte brush is treated the same way except that each bead is charged and there are counterions present to neutralize the charge. We investigated the friction coefficient, monomer density, and brush penetration for both polyelectrolyte and neutral brushes with both equal grafting density and equal normal force under good solvent conditions. We found that polyelectrolyte brushes had a smaller friction coefficient and monomer penetration than neutral polymer brushes with the identical grafting density and chain length, and the polyelectrolyte brushes supported a much higher normal load than the neutral brushes for the same degree of compression. Charged and neutral brushes with their grafting densities chosen so that they support the same load exhibited approximately the same degree of interpenetration, but the polyelectrolyte brush exhibited a significantly lower friction coefficient. We present evidence that the reason for this is that the extra normal force contribution provided by the counterion osmotic pressure that exists for polyelectrolyte brushes permits them to support the same load as an identical neutral polymer brush of higher grafting density. Because of the resulting lower monomer density for the charged brushes, fewer monomer collisions take place per unit time, resulting in a lower friction coefficient.

Many eukaryotic and prokaryotic cells can exhibit remarkable sensing ability under small gradients of chemical compounds. In this study, we approach this phenomenon by considering the contribution of multiple ligands to the chemical kinetics within the Michaelis-Menten model. This work was inspired by the recent theoretical findings of Hu et al. [Phys. Rev. Lett. 105, 048104 (2010)]. Our treatment with practical binding energies and chemical potentials provides results that are consistent with experimental observations.