I. P. Omelyan

Optimized Verlet-like algorithms for molecular dynamics simulations.

Explicit velocity- and position-Verlet-like algorithms of the second order are proposed to integrate the equations of motion in many-body systems. The algorithms are derived on the basis of an extended decomposition scheme at the presence of a free parameter. The nonzero value for this parameter is obtained by reducing the influence of truncated terms to a minimum. As a result, the proposed algorithms appear to be more efficient than the original Verlet versions that correspond to a particular case when the introduced parameter is equal to zero. Like the original versions, the extended counterparts are symplectic and time reversible, but lead to an improved accuracy in the generated solutions at the same overall computational costs. The advantages of the optimized algorithms are demonstrated in molecular dynamics simulations of a Lennard-Jones fluid.

Construction of high-order force-gradient algorithms for integration of motion in classical and quantum systems

A consequent approach is proposed to construct symplectic force-gradient algorithms of arbitrarily high orders in the time step for precise integration of motion in classical and quantum mechanics simulations. Within this approach the basic algorithms are first derived up to the eighth order by direct decompositions of exponential propagators and further collected using an advanced composition scheme to obtain the algorithms of higher orders. Contrary to the scheme proposed by Chin and Kidwell [Phys. Rev. E 62, 8746 (2000)], where high-order algorithms are introduced by standard iterations of a force-gradient integrator of order four, the present method allows one to reduce the total number of expensive force and its gradient evaluations to a minimum. At the same time, the precision of the integration increases significantly, especially with increasing the order of the generated schemes. The algorithms are tested in molecular dynamics and celestial mechanics simulations. It is shown, in particular, that the efficiency of the advanced fourth-order-based algorithms is better approximately in factors 5 to 1000 for orders 4 to 12, respectively. The results corresponding to sixth- and eighth-order-based composition schemes are also presented up to the sixteenth order. For orders 14 and 16, such highly precise schemes, at considerably smaller computational costs, allow to reduce unphysical deviations in the total energy up in 100 000 times with respect to those of the standard fourth-order-based iteration approach.

ALGORITHM FOR NUMERICAL INTEGRATION OF THE RIGID-BODY EQUATIONS OF MOTION

An algorithm for numerical integration of the rigid-body equations of motion is proposed. The algorithm uses the leapfrog scheme and the quantities involved are angular velocities and orientational variables that can be expressed in terms of either principal axes or quaternions. Due to specific features of the algorithm, orthonormality and unit norms of the orientational variables are integrals of motion, despite an approximate character of the produced trajectories. It is shown that the method presented appears to be the most efficient among all such algorithms known.

Optimized Forest-Ruth- and Suzuki-like algorithms for integration of motion in many-body systems

Abstract An approach is proposed to improve the efficiency of fourth-order algorithms for numerical integration of the equations of motion in molecular dynamics simulations. The approach is based on an extension of the decomposition scheme by introducing extra evolution subpropagators. The extended set of parameters of the integration is then determined by reducing the norm of truncation terms to a minimum. In such a way, we derive new explicit symplectic Forest–Ruth- and Suzuki-like integrators and present them in time-reversible velocity and position forms. It is proven that these optimized integrators lead to the best accuracy in the calculations at the same computational cost among all possible algorithms of the fourth order from a given decomposition class. It is shown also that the Forest–Ruth-like algorithms, which are based on direct decomposition of exponential propagators, provide better optimization than their Suzuki-like counterparts which represent compositions of second-order schemes. In particular, using our optimized Forest–Ruth-like algorithms allows us to increase the efficiency of the computations by more than ten times with respect to that of the original integrator by Forest and Ruth, and by approximately five times with respect to Suzukis approach. The theoretical predictions are confirmed in molecular dynamics simulations of a Lennard–Jones fluid. A special case of the optimization of the proposed Forest–Ruth-like algorithms to celestial mechanics simulations is considered as well.

Algorithm for molecular dynamics simulations of spin liquids

A new symplectic time-reversible algorithm for numerical integration of the equations of motion in magnetic liquids is proposed. It is tested and applied to molecular dynamics simulations of a Heisenberg spin fluid. We show that the algorithm exactly conserves spin lengths and can be used with much larger time steps than those inherent in standard predictor-corrector schemes. The results obtained for time correlation functions demonstrate the evident dynamic interplay between the liquid and magnetic subsystems.

Numerical integration of the equations of motion for rigid polyatomics: The matrix method

A new scheme for numerical integration of motion for classical systems composed of rigid polyatomic molecules is proposed. The scheme is based on a matrix representation of the rotational degrees of freedom. The equations of motion are integrated within the Verlet framework in velocity form. It is shown that, contrary to previous methods, in the approach introduced the rigidity of molecules can be conserved automatically without any additional transformations. A comparison of various techniques with respect to numerical stability is made.

Phase diagrams of classical spin fluids: The influence of an external magnetic field on the liquid-gas transition

The influence of an external magnetic field on the liquid-gas phase transition in Ising, XY, and Heisenberg spin fluid models is studied using a modified mean field theory and Gibbs ensemble Monte Carlo simulations. It is demonstrated that the theory is able to reproduce quantitatively all characteristic features of the field dependence of the critical temperature T(c)(H) for all the three models. These features include a monotonic decrease of T(c) with rising H in the case of the Ising fluid as well as a more complicated nonmonotonic behavior for the XY and Heisenberg models. The nonmonotonicity consists in a decrease of T(c) with increasing H at weak external fields, an increase of T(c) with rising H in the strong field regime, and the existence of a minimum in T(c)(H) at intermediate values of H. Analytical expressions for T(c)(H) in the large field limit are presented as well. The paramagnetic-ferromagnetic phase transition is also considered in simulations and described within the mean field theory.

Ferromagnetic phase transition in a Heisenberg fluid: Monte Carlo simulations and Fisher corrections to scaling.

The magnetic phase transition in a Heisenberg fluid is studied by means of the finite size scaling technique. We find that even for larger systems, considered in an ensemble with fixed density, the critical exponents show deviations from the expected lattice values similar to those obtained previously. This puzzle is clarified by proving the importance of the leading correction to the scaling that appears due to Fisher renormalization with the critical exponent equal to the absolute value of the specific heat exponent alpha. The appearance of such new corretions to scaling is a general feature of systems with constraints.

XY Spin Fluid in an External Magnetic Field

A method of integral equations is developed to study anisotropic fluids with planar spins in an external field. As a result, the calculations for these systems appear to be no more difficult than those for ordinary homogeneous liquids. The approach proposed is applied to the ferromagnetic XY spin fluid in a magnetic field using a soft mean spherical closure and the Born-Green-Yvon equation. This provides an accurate reproduction of the complicated phase diagram behavior obtained by cumbersome Gibbs ensemble simulation and multiple histogram reweighting techniques.

Wavevector- and frequency-dependent shear viscosity of water: the modified collective mode approach and molecular dynamics calculations

The transverse momentum time autocorrelation functions and wavevectorand frequency-dependent shear viscosity are calculated for an interaction site model of water using a modified collective mode approach and molecular dynamics simulations. The modified mode approach is based on a formulation which consistently takes into account non-Markovian effects into the kinetic memory kernels. As is demonstrated by comparing the theory results with the molecular dynamics data, the entire frequency dependence of the shear viscosity can be reproduced quantitatively over the whole wavelength range in terms of six generalized collective modes employing the kinetic memory kernel in the non-Markovian approximation of the third order. It is also shown that the results corresponding to the exact atomic and abbreviated molecular descriptions may differ considerably. In the infinite wavevector regime the dynamic correlations are completely determined by a single free motion of the molecules.