Classical Langevin dynamics of a charged particle moving on a sphere and diamagnetism: A surprise

It is generally known that the orbital diamagnetism of a classical system of charged particles in thermal equilibrium is identically zero —the Bohr-van Leeuwen theorem. Physically, this null result derives from the exact cancellation of the orbital diamagnetic moment associated with the complete cyclotron orbits of the charged particles by the paramagnetic moment subtended by the incomplete orbits skipping the boundary in the opposite sense. Motivated by this crucial but subtle role of the boundary, we have simulated here the case of a finite but unbounded system, namely that of a charged particle moving on the surface of a sphere in the presence of an externally applied uniform magnetic field. Following a real space-time approach based on the classical Langevin equation, we have computed the orbital magnetic moment that now indeed turns out to be non-zero and has the diamagnetic sign. To the best of our knowledge, this is the first report of the possibility of finite classical diamagnetism in principle, and it is due to the avoided cancellation.

Brownian motion of spins; generalized spin Langevin equation

We derive the Langevin equations for a spin interacting with a heat bath, starting from a fully dynamical treatment. The obtained equations are non-Markovian with multiplicative fluctuations and concommitant dissipative terms obeying the fluctuation-dissipation theorem. In the Markovian limit our equations reduce to the phenomenological equations proposed by Kubo and Hashitsume. The perturbative treatment on our equations lead to Landau-Lifshitz equations and to other known results in the literature.

Fluctuation-dissipation relation for nonlinear Langevin equations

It is shown that the fluctuation-dissipation theorem is satisfied by the solutions of a general set of nonlinear Langevin equations with a quadratic free-energy functional (constant susceptibility) and field-dependent kinetic coefficients, provided the kinetic coefficients satisfy the Onsager reciprocal relations for the irreversible terms and the antisymmetry relations for the reversible terms. The analysis employs a perturbation expansion of the nonlinear terms, and a functional integral calculation of the correlation and response functions, and it is shown that the fluctuation-dissipation relation is satisfied at each order in the expansion.

Internal noise-driven generalized Langevin equation from a nonlocal continuum model

Starting with a micropolar formulation, known to account for nonlocal microstructural effects at the continuum level, a generalized Langevin equation (GLE) for a particle, describing the predominant motion of a localized region through a single displacement degree of freedom, is derived. The GLE features a memory-dependent multiplicative or internal noise, which appears upon recognizing that the microrotation variables possess randomness owing to an uncertainty principle. Unlike its classical version, the present GLE qualitatively reproduces the experimentally measured fluctuations in the steady-state mean square displacement of scattering centers in a polyvinyl alcohol slab. The origin of the fluctuations is traced to nonlocal spatial interactions within the continuum, a phenomenon that is ubiquitous across a broad class of response regimes in solids and fluids. This renders the proposed GLE a potentially useful model in such cases.

Fluctuation phenomena in electrochemistry part I. The formalism

A unified approach to the problem of electrochemical fluctuations is presented. On the basis of the Langevin procedure, primary noise sources are introduced in the basic phenomenological equations and a discussion of the secondary noise sources arising in the expressions for the power spectra of currents is given.

The fully implicit stochastic-alpha method for stiff stochastic differential equations

A fully implicit integration method for stochastic differential equations with significant multiplicative noise and stiffness in both the drift and diffusion coefficients has been constructed, analyzed and illustrated with numerical examples in this work. The method has strong order 1.0 consistency and has user-selectable parameters that allow the user to expand the stability region of the method to cover almost the entire drift-diffusion stability plane. The large stability region enables the method to take computationally efficient time steps. A system of chemical Langevin equations simulated with the method illustrates its computational efficiency.

A Model of Anomalous Chain Translocation Dynamics

A model of polymer translocation based on the stochastic dynamics of the number of monomers on one side of a pore-containing surface is formulated in terms of a one-dimensional generalized Langevin equation, in which the random force is assumed to be characterized by long-ranged temporal correlations. The model is introduced to rationalize anomalies in measured and simulated values of the average time of passage through the pore, which in general cannot be satisfactorily accounted for by simple Brownian diffusion mechanisms. Calculations are presented of the mean first passage time for barrier crossing and of the mean square displacement of a monomeric segment, in the limits of strong and weak diffusive bias. The calculations produce estimates of the exponents in various scaling relations that are in satisfactory agreement with available data.

Polymer melt dynamics: Microscopic roots of fractional viscoelasticity

The rheological properties of polymer melts and other complex macromolecular fluids are often successfully modeled by phenomenological constitutive equations containing fractional differential operators. We suggest a molecular basis for such fractional equations in terms of the generalized Langevin equation (GLE) that underlies the renormalized Rouse model developed by Schweizer [J. Chem. Phys. 91, 5802 (1989)]. The GLE describes the dynamics of the segments of a tagged chain under the action of random forces originating in the fast fluctuations of the surrounding polymer matrix. By representing these random forces as fractional Gaussian noise, and transforming the GLE into an equivalent diffusion equation for the density of the tagged chain segments, we obtain an analytical expression for the dynamic shear relaxation modulus G(t), which we then show decays as a power law in time. This power-law relaxation is the root of fractional viscoelastic behavior.

Observation of a power-law memory kernel for fluctuations within a single protein molecule

The fluctuation of the distance between a fluorescein-tyrosine pair within a single protein complex was directly monitored in real time by photoinduced electron transfer and found to be a stationary, time-reversible, and non-Markovian Gaussian process. Within the generalized Langevin equation formalism, we experimentally determine the memory kernel K(t), which is proportional to the autocorrelation function of the random fluctuating force. K(t) is a power-law decay, t(-0.51 +/- 0.07) in a broad range of time scales (10(-3)-10 s). Such a long-time memory effect could have implications for protein functions.

Hysteresis in model spin systems

The authors study the hysteretic response of model spin systems to periodic time-varying fields H(t) as a function of the amplitude H0 and the frequency Omega . At fixed H0, they find conventional, squarish hysteresis loops at low Omega , and rounded, roughly elliptical loops at high Omega , in agreement with experiment. For the O(N to infinity ), d=3, ( Phi 2)2 model with Langevin dynamics, they find a novel scaling behaviour for the area A of the hysteresis loop, of the form (valid for low fields) A approximately=H0066 Omega 0.33.

Temperature dependent free energy surface of polymer folding from equilibrium and quench studies

Langevin dynamics simulation studies have been employed to calculate the temperature dependent free energy surface and folding characteristics of a 500 monomer long linear alkane (polyethylene) chain with a realistic interaction potential. Both equilibrium and temperature quench simulation studies have been carried out. Using the shape anisotropy parameter (S) of the folded molecule as the order parameter, we find a weakly first order phase transition between the high-temperature molten globule and low-temperature rodlike crystalline states separated by a small barrier of the order of k(B)T. Near the melting temperature (580 K), we observe an intriguing intermittent fluctuation with pronounced ``1/f noise characteristics'' between these two states with large difference in shape and structure. We have also studied the possibilities of different pathways of folding to states much below the melting point. At 300 K starting from the all-trans linear configuration, the chain folds stepwise into a very regular fourfold crystallite with very high shape anisotropy. Whereas, when quenched from a high temperature (900 K) random coil regime, we identify a two step transition from the random coiled state to a molten globulelike state and, further, to a anisotropic rodlike state. The trajectory reveals an interesting coupling between the two order parameters, namely, radius of gyration (R-g) and the shape anisotropy parameter (S). The rodlike final state of the quench trajectory is characterized by lower shape anisotropy parameter and significantly larger number of gauche defects as compared to the final state obtained through equilibrium simulation starting from all-trans linear chain. The quench study shows indication of a nucleationlike pathway from the molten globule to the rodlike state involving an underlying rugged energy landscape. (C) 2010 American Institute of Physics. doi:10.1063/1.3509398]

Hysteresis in model spin system

The aim of this paper is to construct a nonequilibrium statistical‐mechanics theory to study hysteresis in ferromagnetic systems. We study the hysteretic response of model spin systems to periodic magnetic fields H(t) as a function of the amplitude H0 and frequency Ω. At fixed H0, we find conventional, squarelike hysteresis loops at low Ω, and rounded, roughly elliptical loops at high Ω, in agreement with experiments. For the O(N→∞), d=3, (Φ2)2 model with Langevin dynamics, we find a novel scaling behavior for the area A of the hysteresis loop, of the form A∝H0.660Ω0.33. We carry out a Monte Carlo simulation of the hysteretic response of the two‐dimensional, nearest‐neighbor, ferromagnetic Ising model. These results agree qualitatively with the results obtained for the O(N) model.

Nonlinear hydrodynamics of a hard-sphere fluid near the glass transition

We conduct a numerical study of the dynamic behavior of a dense hard-sphere fluid by deriving and integrating a set of Langevin equations. The statics of the system is described by a free-energy functional of the Ramakrishnan-Yussouff form. We find that the system exhibits glassy behavior as evidenced through a stretched exponential decay and a two-stage relaxation of the density correlation function. The characteristic times grow with increasing density according to the Vogel-Fulcher law. The wave-number dependence of the kinetics is extensively explored. The connection of our results with experiment, mode-coupling theory, and molecular-dynamics results is discussed.

A structural study of Li2S---B2S3 glasses by neutron diffraction

The radial distribution functions (RDFs) of five xLi2S.(1 - x)B2S3 glasses (x = 0.55, 0.60, 0.67, 0.71 and 0.75) have been determined from neutron diffraction experiments performed at the Institut Laue-Langevin, Grenoble. These glasses are prepared by casting a molten mixture of boron, sulphur and Li2S inside a controlled atmosphere glovebox. Addition of the Li2S Modifier is found gradually to suppress all peaks corresponding to interatomic distances > 3.5 angstrom, which implies that the structural entities present in these glasses become segmented, and therefore more ionic, as x increases. The assumption of the existence of four main structural entities based on four- and three-coordinated borons (the latter carrying bridging and/or non-bridging sulphurs) accounts for all the peaks present in the RDFs as a function of composition. Furthermore, in the most modified glass (x = 0.75), that which contains only 'isolated' BS33- triangles, there seems to be evidence for either octahedral or tetrahedral coordination of Li+ by S- ions