On numerical stationary distribution of overdamped Langevin equation in harmonic system

Efficient numerical algorithm for stochastic differential equation has been an important object in the research of statistical physics and mathematics for a long time.In this work we study the highly accurate numerical algorithm for the overdamped Langevin equation.In particular,our interest is in the behaviour of the numerical schemes for solving the overdamped Langevin equation in the harmonic system.Based on the large friction limit of the underdamped Langevin dynamic scheme,three algorithms for overdamped Langevin equation are obtained.We derive the explicit expression of the stationary distribution of each algorithm by analysing the discrete time trajectory for both one-dimensional case and multi-dimensional case.The accuracy of the stationary distribution of each algorithm is illustrated by comparing with the exact Boltzmann distribution.Our results demonstrate that the“BAOA-limit”algorithm generates an accurate distribution of the harmonic system in a canonical ensemble,within a stable range of time interval.The other algorithms do not produce the exact distribution of the harmonic system.

Nonequilibrium Statistical Physics Subject to the Anomalous Langevin Equation in Liouville Space

Considering that thermodynarmic irreversibility and hydrodynamic equations can not be derived rigorously and unifiedly from the Liouville equations, the anomalous Langevin equation in the Liouville space is proposed as a fundamental equation of statistical physics. This equation reflects that the law of motion of particles obeying reversible, deterministic laws in dynamics becomes irreversible and stochastic in thermodynamics. From this the fundamental equations of nonequilibrium thermodynamics, the principle of entropy increase and the theorem of minimum entropy production have been derived. The hydrodynamic equations, such as the generalized Navier-Stokes equation and the mass drift-diffusion equation etc. have been derived rigorously from the kinetic kinetic equation which is reduced from the anomalous Langevin equation in Liouville space. All these are unified and self consistent. But it is difficult to prove that entropy production density σ can never be negative everywhere for all the isolated inhomogeneous systems far from equilibrium.

From Chemical Langevin Equations to Fokker-Planck Equation: Application of Hodge Decomposition and Klein-Kramers Equation

The stochastic systems without detailed balance are common in various chemical reaction systems, such as metabolic network systems. In studies of these systems, the concept of potential landscape is useful However, what are the su^cient and necessary conditions of the existence of the potential function is still an open problem. Use Hodge decomposition theorem in differential form theory, we focus on the general chemical Langevin equations, which reitect complex chemical reaction systems. We analysis the conditions for the existence of potential landscape of the systems. By mapping the stochastic differential equations to a Hamiltonian mechanical system, we obtain the Fokker-Planck equation of the chemical reaction systems. The obtained Fokker-Planck equation can be used in further studies of other steady properties of complex chemical reaction systems, such as their steady state entropies.

Analytical solution of a multidimensional Langevin equation at high friction limits and probability passing over a two-dimensional saddle

The analytical solution of a multidimensional Langevin equation at the overdamping limit is obtained and the probability of particles passing over a two-dimensional saddle point is discussed. These results may break a path for studying further the fusion in superheavy elements synthesis.

Temperature dependent fission fragment distribution in the Langevin equation

The temperature dependent width of the fission fragment distributions was simulated in the Langevin equation by taking two-parameter exponential form of the fission fragment mass variance at scission point for each fission event. The result can reproduce experimental data well, and it permits to make reliable estimate for unmeasured product yields near symmetry fission.

Electron Transfer Between Molecules: Study with Langevin Equations

In this paper, we investigate the electron transfer(ET) in donor-acceptor model. The Langevin equation with random forces is used. The oscillations of the primary states observed in experimental data have been shown with this approach. And other features on the dependence of the rate of ET on temperature, free energy, and reorganization energy have also been clearly shown.

Dynamic of the Dielectric Nano-Particle in Optical Tweezer Using Counter-Propagating Pulsed Laser Beams

In this article, the dynamical process of the dielectric particle in the optical tweezer using the counter-propagating Gaussian pulses is investigated by the Langevin equation concerning the Brownian motion. The temporal stabilities of particle is simulated. The influence of the duration, repetition period and delay time between pulses on stability is discussed.

The Third Law of Quantum Thermodynamics in the Presence of Anomalous Couplings

The quantum thermodynamic functions of a harmonic oscillator coupled to a heat bath through velocity-dependent coupling are obtained analytically. It is shown that both the free energy and the entropy decay fast with the temperature in relation to that of the usual coupling from. This implies that the velocity-dependent coupling helps to ensure the third law of thermodynamics.

Derivation of persistent time for anisotropic migration of cells

Cell migration plays an essential role in a wide variety of physiological and pathological processes. In this paper we numerically discuss the properties of an anisotropic persistent random walk （APRW） model, in which two different and independent persistent times are assumed for cell migrations in the x-and y-axis directions. An intrinsic orthogonal coordinates with the primary and non-primary directions can be defined for each migration trajectory based on the singular vector decomposition method. Our simulation results show that the decay time of single exponential distribution of velocity auto-correlation function （VACF） in the primary direction is actually the large persistent time of the APRW model, and the small decay time of double exponential VACF in the non-primary direction equals the small persistent time of the APRW model. Thus, we propose that the two persistent times of anisotropic migration of cells can be properly estimated by discussing the VACFs of trajectory projected to the primary and non-primary directions.

Roles of Isospin in Evaporation Residue Cross Section as a Probe of Nuclear Dissipation

The influence of isospin on the excess of evaporation residue cross section over its standard statistical-model value for nuclei ^194pb, ^200Pb, and ^206pb is studied via a Langevin equation coupled with a statistical decay model. The magnitude of this excess for a low-isospin fissioning nucleus is shown to be larger and its dependence on the nuclear viscosity coefficient to be stronger than those of a high-isospin fissioning nucleus. These results suggest that to obtain a more accurate information of viscosity coefficient inside the saddle point by measuring evaporation residue cross sections, we had better choose those compound systems with small isospin.

Intrinsic noise analysis and stochastic simulation on transforming growth factor beta signal pathway

A typical biological cell lives in a small volmne at room temperature; the noise effect on the cell signal transduction pathway may play an important role in its dynamics. Here, using the transforming growth factor-β signal transduction pathway as an example, we report our stochastic simulations of the dynamics of the pathway and introduce a linear noise approximation method to calculate the transient intrinsic noise of pathway components. We compare the numerical solutions of the linear noise approximation with the statistic results of chemical Langevin equations, and find that they are quantitatively in agreement with the other. When transforming growth factor-β dose decreases to a low level, the time evolution of noise fluctuation of nuclear Smad2-Smad4 complex indicates the abnormal enhancement in the transient signal activation process.

Effects of N/Z on Spin Distribution of Evaporation Residue Cross Section as a Probe of Nuclear Dissipation

The spin distribution of the evaporation residue cross section of nuclei ^194pb, ^200Pb, ^206Pb, and ^200 Os are calculated via a Langevin equation coupled with a statistical decay model. It is shown that with increasing the neutronto-proton ratio （N/Z） of the system, the sensitivity of the significantly. Moreover, for ^200Os this spin distribution is no spin distribution to the nuclear dissipation is decreased longer sensitive to the nuclear dissipation. These results suggest that to obtain a more accurate pre-saddle viscosity coefficient through the measurement of the evaporation residue spin distribution, it is best to yield those compound systems with low N/Z.

A Detailed Study of Pre-scission γ Emission as a Probe of Nuclear Dissipation

Using a Langevin equation coupled with a statistical model, we calculate pre-scission giant dipole resonance （GDR） γ-ray multiplicity of nuclei 194 pb, 200Pb, 206Pb, and 200 Os. It is demonstrated that with increasing the isospin asymmetry of these fissioning nuclei the sensitivity of the emitted γ multiplicity to the nuclear viscosity coefficient is decreased significantly. For 200Os nuc/eus, this γ-ray emission is no longer sensitive to the magnitude of the viscosity coefficient. In addition, the effect of the isospin asymmetry on the γ rays as a probe of nuclear dissipation is reduced with increasing angular momentum. These results suggest that to obtain a more accurate information of the viscosity coemfficient by the measurement of pre-scission GDR γ-ray multiplicity it is better to choose those compound systems with small isospin asymmetry and low spin.

On the nonclassical dynamics of cavity-assisted four-channel nonlinear coupler

The non-classical properties of light propagating in four-channel Kerr waveguides, confined in an optical cavity, are studied. The solution to the Hamiltonian of field operators is obtained semi-analytically by using symmetrically ordered phase-space representation. Full quantum analysis of the input coherent fields displays a strong transition of photon property between the super-Poissonian and sub-Poissonian statistics. It is found that the cavity-assisted multichannel system exhibits enhanced squeezing both in single-and compound-mode. This multichannel system may be utilized as an efficient quantumlight generator.

Diffusion induced by bounded noise in a two-dimensional coupled memory system

The diffusion behavior driven by bounded noise under the influence of a coupled harmonic potential is investigated in a two-dimensional coupled-damped model. With the help of the Laplace analysis we obtain exact descriptions for a particle’s two-time dynamics which is subjected to a coupled harmonic potential and a coupled damping. The time lag is used to describe the velocity autocorrelation function and mean square displacement of the diffusing particle. The diffusion behavior for the time lag is also discussed with respect to the coupled items and the amplitude of bounded noise.

Dynamic of Polystyrene Microsphere Linking to DNA Molecule under Optical Tweezer

Up to now there are many works investigated the stretched state of deoxyribonucleic axid （DNA） by optical trap. The main purpose of mentioned works is to find out the elastic force-extension characterisitic approximatelly based on the experimental parameters. In experiments, the optical force is seem as an apply force, but not an own elastic force of DNA molecules. To use sufficiently optical tweezer to manipulate and hold the DNA molecules keeping in stretched state, it is necessary to use the elastic as one force contributing to process trapping it, and consequently, the term describing elastic force must be present in the general Langevin equation （GLE） describing the dynamic of driven bead linking to DNA molecule. In this paper, the finite difference equation of GLE is derived for simulation the dynamic of bead linking to DNA molecule embedded in the fluid. The process manipulation of polystyrene bead to tweezer center is simulated and discussed. From results, the role of Brownian force, elastic force and optical force on process trapping driven bead is evaluated. Moreover, the trapping time, velocity of bead are discussed, too.

A New Insight into Energy Distribution of Electrons in Fuel-Rod Gap in VVER-1000 Nuclear Reactor

In order to calculate the electron energy distribution in the fuel rod gap of a VVER- 1000 nuclear reactor, the Fokker-Planck equation （FPE） governing the non-equilibrium behavior of electrons passing through the fuel-rod gap as an absorber has been solved in this paper. Besides, the Monte Carlo Geant4 code was employed to simulate the electron migration in the fuel-rod gap and the energy distribution of electrons was found. As for the results, the accuracy of the FPE was compared to the Geant4 code outcomes and a satisfactory agreement was found. Also, different percentage of the volatile and noble gas fission fragments produced in fission reactions in fuel rod, i.e. Krypton, Xenon, Iodine, Bromine, Rubidium and Cesium were employed so as to investigate their effects on the electrons＇ energy distribution. The present results show that most of the electrons in the fuel rod＇s gap were within the thermal energy limitation and the tail of the electron energy distribution was far from a Maxwellian distribution. The interesting outcome was that the electron energy distribution is slightly increased due to the accumulation of fission fragments in the gap. It should be noted that solving the FPE for the energy straggling electrons that are penetrating into the fuel-rod gap in the VVER-1000 nuclear reactor has been carried out for the first time using an analytical approach.

NONLINEAR LANGEVIN MODEL WITH PRODUCT STOCHASTICITY FOR BIOLOGICAL NETWORKS： THE CASE OF THE SCHNAKENBERG MODEL

Langevin equation is widely used to study the stochastic effects in molecular networks, as it often approximates well the underlying chemical master equation. However, frequently it is not clear when such an approximation is applicable and when it breaks down. This paper studies the simple Schnakenberg model consisting of three reversible reactions and two molecular species whose concentrations vary. To reduce the residual errors from the conventional formulation of the Langevin equation, the authors propose to explicitly model the effective coupling between macroscopic concentrations of different molecular species. The results show that this formulation is effective in correcting residual errors from the original uncoupled Langevin equation and can approximate the underlying chemical master equation very accurately.

Langevin study of neutron emission in the reactions ^(16)O+^(181)Ta and ^(19)F+^(178)Hf

The pre-scission neutrons measured in the reactions ^16O＋^181Ta and ^19F＋^178Hf are studied via a Langevin equation coupled with a statistical decay model. We find that because of the mass asymmetry of different entrance channels, the spin distributions of compound nuclei would be different, consequently, the measured neutrons in these two reactions would also different. This means that the entrance channel will affect the particle emission in the fission process of hot nuclei.

Investigation of fission fragments average spin based on four dimensional Langevin dynamical model

We applied the four dimensional Langevin dynamical model to investigate the average spin of fission fragments. Elongation, neck thickness, asymmetry parameter, and the orientation degree of freedom（K coordinate）are the four dimensions of the dynamical model. We assume that the collective modes depend on the emission angle of the fragments, then different parameters related to the average spin of fission fragments are calculated dynamically.The angle dependence of average spin of fission fragments is investigated by calculating the spin at angles 90？and165？. Also, the obtained results based on the transition state model at scission point are presented. One can obtain better agreement between the results of the dynamical model and experimental data in comparison with the results of the transition state model.