The super long-time decay of velocity fluctuations in a two-dimensional fluid

We have calculated the velocity autocorrelation function for a tracer particle in a model two-dimensional fluid. The fluid was represented by a lattice Boltzmann equation with imposed fluctuations. By choosing a low Boltzmann diffusion coefficient for the tracer, the diverging contribution to the diffusion coefficient can be made to exceed the Boltzmann value even at short times. We were thus able to find evidence for the renormalized, or ‘super long-time’, decay of the VACF in a two-dimensional fluid. We find quantitative evidence for the 1/t√ln(t) decay predicted by theory.     

Dissipative particle dynamics study of velocity autocorrelation function and self-diffusion coefficient in terms of interaction potential strength

This research focuses on numerically investigating the self-diffusion coefficient and velocity autocorrelation function (VACF) of a dissipative particle dynamics (DPD) fluid as a function of the conservative interaction strength. Analytic solutions to VACF and self-diffusion coefficients in DPD were obtained by many researchers in some restricted cases including ideal gases, without the account of conservative force. As departure from the ideal gas conditions are accentuated with increasing the relative proportion of conservative force, it is anticipated that the VACF should gradually deviate from its normally expected exponentially decay. This trend is confirmed through numerical simulations and an expression in terms of the conservative force parameter, density and temperature is proposed for the self-diffusion coefficient. As it concerned the VACF, the equivalent Langevin equation describing Brownian motion of particles with a harmonic potential is adapted to the problem and reveals an exponentially decaying oscillatory pattern influenced by the conservative force parameter, dissipative parameter and temperature. Although the proposed model for obtaining the self-diffusion coefficient with consideration of the conservative force could not be verified due to computational complexities, nonetheless the Arrhenius dependency of the self-diffusion coefficient to temperature and pressure permits to certify our model over a definite range of DPD parameters.     

Stochastic One-Dimensional Lorentz Gas on a Lattice

We study a one-dimensional stochastic Lorentz gas where a light particle moves in a fixed array of nonidentical random scatterers arranged in a lattice. Each scatterer is characterized by a random transmission/reflection coefficient. We consider the case when the transmission coefficients of the scatterers are independent identically distributed random variables. A symbolic program is presented which generates the exact velocity autocorrelation function (VACF) in terms of the moments of the transmission coefficients. The VACF is found for different types of disorder for times up to 20 collision times. We then consider a specific type of disorder: a two-state Lorentz gas in which two types of scatterers are arranged randomly in a lattice. Then a lattice point is occupied by a scatterer whose transmission coefficient is with probability p or + with probability 1–p. A perturbation expansion with respect to is derived. The ² term in this expansion shows that the VACF oscillates with time, the period of oscillation being twice the time of flight from one scatterer to its nearest neighbor. The coarse-grained VACF decays for long times like t ^–3/2, which is similar to the decay of the VACF of the random Lorentz gas with a single type of scatterer. The perturbation results and the exact ones (found up to 20 collision times) show good agreement.     

Measurement of angular dependence of M X-ray production cross-sections in Re,Bi and U at 5.96 keV

The M X-ray production differential cross sections in Re, Bi and U elements have been measured at the 5.96 keV incident photon energy in an angular range 135°–155°. The measurements were performed using a ⁵⁵Fe source and a Si(Li) detector. The present results contradict the predictions of Cooper and Zare [Atomic Collision Processes, Gordon and Breach, New York (1969)] and experimental results of Kumar et al. [J. Phys. B: At. Mol. Opt. 34, 613 (2001)]. that, after photoionization of inner shells, the vacancy state has equal population of magnetic substates and the subsequent X-ray emission is isotropic, but confirm the predictions of the calculations of Flügge et al. [Phys. Rev. Lett. 29, 7 (1972)] and experimental results of Sharma and Allawadhi [J. Phys. B: At. Mol. Opt. 32, 2343 (1999)] and Ertugrul [Nucl. Instrum. Meth. B 119, 345 (1996)]. Total M X-ray production cross sections from the decay at the 5.96 keV photon energies are found to be in good agreement with the calculated theoretical results using the theoretical values of M shell photoionization cross section.     

A hybrid formalism combining fluctuating hydrodynamics and generalized Langevin dynamics for the simulation of nanoparticle thermal motion in an incompressible fluid medium

A novel hybrid scheme based on Markovian fluctuating hydrodynamics of the fluid and a non-Markovian Langevin dynamics with the Ornstein-Uhlenbeck noise perturbing the translational and rotational equations of motion of the nanoparticle is employed to study the thermal motion of a nanoparticle in an incompressible Newtonian fluid medium. A direct numerical simulation adopting an arbitrary Lagrangian-Eulerian (ALE) based finite element method (FEM) is employed in simulating the thermal motion of a particle suspended in the fluid confined in a cylindrical vessel. The results for thermal equilibrium between the particle and the fluid are validated by comparing the numerically predicted temperature of the nanoparticle with that obtained from the equipartition theorem. The nature of the hydrodynamic interactions is verified by comparing the velocity autocorrelation function (VACF) and mean squared displacement (MSD) with well-known analytical results. For nanoparticle motion in an incompressible fluid, the fluctuating hydrodynamics approach resolves the hydrodynamics correctly but does not impose the correct equipartition of energy based on the nanoparticle mass because of the added mass of the displaced fluid. In contrast, the Langevin approach with an appropriate memory is able to show the correct equipartition of energy, but not the correct short- and long-time hydrodynamic correlations. Using our hybrid approach presented here, we show for the first time, that we can simultaneously satisfy the equipartition theorem and the (short- and long-time) hydrodynamic correlations. In effect, this results in a thermostat that also simultaneously preserves the true hydrodynamic correlations. The significance of this result is that our new algorithm provides a robust computational approach to explore nanoparticle motion in arbitrary geometries and flow fields, while simultaneously enabling us to study carrier adhesion mediated by biological reactions (receptor-ligand interactions) at the vessel wall at a specified finite temperature.     

Analysis of the quantum numbers J(PC) of the X(3872) particle

We present an analysis of angular distributions and correlations of the X(3872) particle in the exclusive decay mode X(3872)-->J/psipi+ pi- with J/psi-->mu+ mu-. We use 780 pb-1 of data from pp[over ] collisions at sqrt[s]=1.96 TeV collected with the CDF II detector at the Fermilab Tevatron. We derive constraints on spin, parity, and charge conjugation parity of the X(3872) particle by comparing measured angular distributions of the decay products with predictions for different J(PC) hypotheses. The assignments J(PC)=1++ and 2-+ are the only ones consistent with the data.     

Long-time tails of the velocity autocorrelation functions for the triangular periodic Lorentz gas

We present numrical results on the velocity autocorrelation function (VACF)C(t)=<ν(t)·ν(0)> for the periodic Lorentz gas on a two-dimensional triangular lattice as a function of the radiusR of the hard disk scatterers on the lattice. Our results for the unbounded horizon case confirm 1/t decay of the VACF for long times (out to 100 times the mean free time between collisions) and provide strong support for the conjecture by Friedman and Martin that the 1/t decay is due to long free paths along which a moving particle does not scatter up to timet. Even after new sets of long free paths become available forR<1/4, we continue to find good agreement between numerical results and an analytically estimated 1/t decay. For the bounded horizon case , our numerical VACFs decay exponentially, although it is difficult to discriminate among pure exponential decay, exponential decay with prefactor, and stretched exponential decay.     

Diffusion within a near-critical nanopore fluid

Results are presented for grand canonical Monte Carlo (GCMC) and both equilibrium and non-equilibrium molecular dynamics simulations (EMD and NEMD) conducted over a range of densities and temperatures that span the two-phase coexistence and supercritical regions for a pure fluid adsorbed within a model crystalline nanopore. The GCMC simulations provided the low temperature coexistence points for the open pore fluid and were used to locate the capillary critical temperature for the system. The equilibrium configurational states obtained from these simulations were then used as input data for the EMD simulations in which the self-diffusion coefficients were computed using the Einstein equation. NEMD colour diffusion simulations were also conducted to validate the use of a system averaged Einstein analysis for this inhomogeneous fluid. In all cases excellent agreement was observed between the equilibrium (linear response theory) predictions for the diffusivities and non-equilibrium colour diffusivities. The simulation results are also compared with a recently published quasi-hydrodynamic theory of Pozhar and Gubbins (Pozhar, L. A., and Gubbins, K. E., 1993, J. Chem. Phys., 99, 8970; 1997, Phys. Rev. E, 56, 5367.). The model fluid and the nature of the fluid wall interactions employed conform to the decomposition of the particle–particle interaction potential explicitly used by Pozhar and Gubbins. The local self-diffusivity was calculated from the local fluid–fluid and fluid wall hard core collision frequencies. While this theory provides reasonable results at moderate pore fluid densities, poor agreement is observed in the low density limit.     

Modeling of super-extreme events: An application to the hierarchical Weierstrass-Mandelbrot Continuous-time Random Walk

We analytically demonstrate and numerically simulate two utmost cases of dragon-kings’ impact on the (unnormalized) velocity autocorrelation function (VACF) of a complex time series generated by stochastic random walker. The first type of dragon-kings corresponds to a sustained drift whose duration time is much longer than that of any other event. The second type of dragon-kings takes the form of an abrupt shock whose amplitude velocity is much larger than those corresponding to any other event. The stochastic process in which the dragon-kings occur corresponds to an enhanced diffusion generated within the hierarchical Weierstrass-Mandelbrot Continuous-time Random Walk (WM-CTRW) formalism. Our analytical formulae enable a detailed study of the impact of the two super-extreme events on the VACF calculated for a given random walk realization on the form of upward deviations from the background power law decay present in the absence of dragon-kings. This allows us to provide a unambiguous distinction between the super-extreme dragon-kings and ‘normal’ extreme “black swans”. The results illustrate diagnostic that could be useful for the analysis of extreme and super-extreme events in real empirical time series.     

Observation of nearly perfect irrotational flow in normal and superfluid strongly interacting Fermi gases

We study the hydrodynamic expansion of a rotating strongly interacting Fermi gas by releasing a cigar-shaped cloud with a known angular momentum from an optical trap. As the aspect ratio of the expanding cloud approaches unity, the angular velocity increases, indicating quenching of the moment of inertia I to as low as 0.05 of the rigid body value I(rig). Remarkably, we observe this behavior in both the superfluid and collisional normal fluid regimes, which obey nearly identical zero-viscosity irrotational hydrodynamics. We attribute irrotational flow in the normal fluid to a decay of the rotational part of the stream velocity during expansion, which occurs when the shear viscosity is negligible. Using conservation of angular momentum, we directly observe a fundamental result of irrotational hydrodynamics, I/I(rig) = delta2, where delta is the deformation parameter of the cloud.     

Numerical entropy and phason elastic constants of plane random tilings with any 2D-fold symmetry

We perform Transition matrix Monte Carlo simulations to evaluate the entropy of rhombus tilings with fixed polygonal boundaries and 2D-fold rotational symmetry. We estimate the large-size limit of this entropy for D=4 to 10. We confirm analytic predictions of [N. Destainville et al., J. Stat. Phys. 120, 799 (2005) and M. Widom et al., J. Stat. Phys. 120, 837 (2005)], in particular that the large size and large D limits commute, and that entropy becomes insensible to size, phason strain and boundary conditions at large D. We are able to infer finite D and finite size scalings of entropy. We also show that phason elastic constants can be estimated for any D by measuring the relevant perpendicular space fluctuations.     

Linear and angular velocity autocorrelation functions for two-dimensional systems of repulsive and Lennard-Jones diatomics

The method of molecular dynamics has been employed to study the autocorrelation functions (ACF's) of the linear velocity and of the angular momentum for two two-dimensional systems of diatomic molecules. The anisotropic potential used is built from four atom-atom interactions. These are either purely repulsive (R), or repulso-attractive (LJ). Emphasis has been put on the long-time behaviour; for this we have used large systems (1600 molecules). The density lies in the intermediate range.

The linear velocity ACF exhibits, for system R, an initial decay which is exponential and whose characteristic time has been checked to be of the order of the mean ‘collision’ time. For system LJ (where collisions cannot be usefully defined) the decay does not follow a simple law, and lies above a single exponential. After an intermediate regime which lasts a long time-span for both systems, the ACF's adopt a t ^-1 decay in agreement with theoretical predictions based on hydrodynamics or other approaches.

The mean square displacement has also been extracted from the data. Its ratio to time, i.e. the apparent coefficient of diffusion, increases logarithmically with time, so that the two kinds of data are in accord. This asymptotic regime starts after about 11 collision times for the repulsive system, and about as soon for the other one.

The angular momentum ACF's, after an initial decay which is closer to an exponential for system R than for system LJ, adopt a t ^-3 behaviour for a very long time-span; this possibly masks a final t ^-2 decay predicted by several theories based on the coupling of the intrinsic momentum with the vorticity of the fluid.     

On the use of multiple interpolation series in scaled particle theory: improved predictions and limitations

We propose further modifications to recent formulations of the scaled particle theory (SPT) of hard particle fluids. In any formulation of SPT, a number of conditions are applied to the central function G. As the exact form of G is unknown, a form must be assumed. Representing G with more than one series has recently been shown (Siderius and Corti, J. Chem. Phys. 127, 144502, 2007) to lead to improved planar surface tension predictions with little change in the already accurate pressure, excess chemical potential and fluid structure predictions. Another condition on G that includes the curvature of the radial distribution at contact was not, however, invoked in this previous work. We incorporate here this additional condition into the multiple series approach, resulting in further improvements to the SPT predictions. We also consider a number of permutations of the multiple interpolation series to explore the limitations of this approach and again find that increasing the number of series will not by itself lead to better predictions of the properties of the hard sphere fluid.     

“Shallow-water” and “deep-water” approximations in the theory of the disruptive instability of thin current-carrying layers

Within the framework of the two-fluid hydrodynamics of plasmas it is shown that the problem instability of a thin current-carrying layer admits two limiting cases which allow analytic solutions and complement one another. These limits are analogous to the well-known shallow-water and deep-water approximations in the fluid mechanical “wave-breaking” instability. In this case, the long-wave limit coincides with the “quasi-Chaplygin” dynamic system of Bulanov and Sasorov, Fiz Plazmy 4, 746 (1978) [Sov. J. Plasma Phys. 4, 418 (1978)], while the short-wavelength limit corresponds to the phenomenological model of Trubnikov, Usp. Fiz. Nauk 160, 167 (1990) [Sov. Phys. Usp. 33, 87 (1990)], for the clumping of “elementary” currents. In the latter case, strong collapse is unavoidable with the appearance of current filaments that trap a finite current. Zh. éksp. Teor. Fiz. 113, 1313–1318 (April 1998)     

Structure factor of a hard-core fluid with short-range Yukawa attraction: analytical FMSA theory against Monte Carlo simulations

ABSTRACT

Analytical solution of the first-order Ornstein–Zernike equation known as the first-order mean spherical approximation (FMSA) theory due to Tang and Lu [J. Chem. Phys. 99, 9828 (1993)] is used to write down a closed equation for the static structure factor of the hard-sphere fluid with a short-range Yukawa attraction. Calculations are performed for a Yukawa decay exponent that corresponds to a range of attraction that does not exceed the first coordination shell of Lennard-Jones-like simple fluids. By comparison with Monte Carlo simulation data it is shown that the analytical FMSA equation for the static structure factor is of the same or even of superior accuracy as that within the seminumerical mean spherical approximation theory.     

On isotropic turbulence in the dark fluid universe

As a first part of this work, experimental information about the decay of isotropic turbulence in ordinary hydrodynamics, [`(u²(t))] μ t^-6/5\overline{\mathbf{u}^{2}(t)}\propto t^{-6/5}, is used as input in FRW equations in order to investigate how an initial fraction f of turbulent kinetic energy in the cosmic fluid influences the cosmological development in the late, quintessence/phantom, universe. First order perturbative theory to the first order in f is employed. It turns out that both in the Hubble factor and in the energy density, the influence from the turbulence fades away at late times. The divergences in these quantities near the Big Rip behave essentially as in a non-turbulent fluid. However, for the scale factor, the turbulence modification turns out to diverge logarithmically. As a second part of our work, we consider the full FRW equation in which the turbulent part of the dark energy is accounted for by a separate term. It is demonstrated that turbulence occurrence may change the future universe evolution due to dissipation of dark energy. For instance, the phantom-dominated universe becomes asymptotically a de Sitter one in the future, thus avoiding the Big Rip singularity.     

用细胞自动机方法模拟计算二维流体的平衡与非平衡相关函数

本文采用一个含有温度的简化点阵气模型,通过细胞自动机方法,模拟计算了服从于外部温度梯度的二维流体中的平衡与非平衡空间相关函数。数值结果表明,用本文方法所计算的结果与由涨落流体动力理论的预言在是性上是完全符合的。     

Economic mapping to the renormalization group scaling of stock markets

We make an attempt to map a simple economically motivated model for price evolution [J. Phys. A 33, 3637 (2000)] to the phenomenological renormalization group scaling of stock markets. This mapping gives insight into the critical exponents and the renormalization group predictions for the log-periodic oscillations preceding some stock market crashes from the perspective of non-linear changes in `the level of stock'. Received 7 August 2000