Backward Euler and other methods for simulating molecular fluids

Wu, Jian; Watts, Robert O.
September 1995
Journal of Chemical Physics;9/1/1995, Vol. 103 Issue 9, p3718
Academic Journal
Two implementations of the backward Euler method for simulating molecular fluids are compared with Brownian dynamics and molecular dynamics simulations of a single diatomic molecule, liquid argon, a single butane molecule, and liquid butane. By comparison with standard molecular dynamics results, backward Euler simulations give different thermodynamic properties for liquids; predict liquid structures which are too solidlike; and incorrectly represent dynamical relaxation processes. The backward Euler methods allows longer time steps to be used in simulations at the cost of an energy minimization at every time step. Even when time steps more than 20 times larger than that required for energy and momentum conservation are used, neither implementation of the backward Euler algorithm is more accurate than standard molecular dynamics calculations with the same time step. We conclude that the new methods offer no computational advantages over more usual methods for simulating molecular fluids and that they often predict incorrect results. © 1995 American Institute of Physics.


Related Articles

  • Brownian dynamics simulation of rotational correlation functions for a three-body macromolecular model. Ying, Ruoxian; Peters, Michael H. // Journal of Chemical Physics;7/15/1991, Vol. 95 Issue 2, p1234 

    A former investigation [R. Ying and M. H. Peters, J. Chem. Phys. 91, 1287 (1989)] on the rotational dynamics of rigid and partially flexible macromolecules in solution by means of Brownian dynamics simulation methods is extended to a trimer or three-body macromolecular model. We present...

  • Rotational Brownian motion of a pair of linear molecules or dipoles with anisotropic interaction. Felderhof, B. U. // Journal of Chemical Physics;8/22/2002, Vol. 117 Issue 8, p3583 

    The rotational Smoluchowski equation for the orientational distribution function of two identical linear molecules or dipoles with anisotropic interaction is solved by expansion in terms of a limited set of eigenfunctions of the operator for free rotational diffusion. The equilibrium...

  • Revisiting angular jump models of molecular reorientations in viscous liquids. Diezemann, Gregor; Sillescu, Hans // Journal of Chemical Physics;7/15/1999, Vol. 111 Issue 3, p1126 

    Studies the model of isotropic reorientational Brownian motion via finite angular steps to the case in which molecular reorientations by different angles take place around arbitrary axes. Analytical expressions for experimentally relevant correlation functions; Comparison of results with...

  • Torque constraints for modeling the behavior of rigid and semirigid macromolecules. Ying, Ruoxian; Peters, Michael H. // Journal of Chemical Physics;7/15/1989, Vol. 91 Issue 2, p1287 

    This work investigates the rotational dynamics of rigid and partially flexible macromolecules in solution by means of Brownian dynamics simulation techniques. In a previous study by Diaz et al. [J. Chem. Phys. 87, 6021 (1987)] on the rotational correlation functions of rigid models, the effect...

  • Krylov subspace methods for computing hydrodynamic interactions in Brownian dynamics simulations. Ando, Tadashi; Chow, Edmond; Saad, Yousef; Skolnick, Jeffrey // Journal of Chemical Physics;8/14/2012, Vol. 137 Issue 6, p064106 

    Hydrodynamic interactions play an important role in the dynamics of macromolecules. The most common way to take into account hydrodynamic effects in molecular simulations is in the context of a Brownian dynamics simulation. However, the calculation of correlated Brownian noise vectors in these...

  • Computer investigations on the mechanism of the influence of quencher concentration on the rate of simple bimolecular reaction. Litniewski, Marek // Journal of Chemical Physics;11/7/2006, Vol. 125 Issue 17, p174501 

    Molecular dynamics investigations on the influence of the concentration of B (quencher) on the rate coefficient, k(t), for the reaction A+B→C+B are continued [M. Litniewski, J. Chem. Phys. 123, 124506 (2005); 124, 114501 (2006)]. The problem is investigated by analyzing the excess in the...

  • Role of Brownian motion hydrodynamics on nanofluid thermal conductivity. Evans, William; Fish, Jacob; Keblinski, Pawel // Applied Physics Letters;2/27/2006, Vol. 88 Issue 9, p093116 

    We use a kinetic theory based analysis of heat flow in fluid suspensions of solid nanoparticles (nanofluids) to demonstrate that the hydrodynamics effects associated with Brownian motion have only a minor effect on the thermal conductivity of the nanofluid. This analysis is supported by the...

  • Comment on constrained molecular dynamics of macromolecules. Toxvaerd, So\ren // Journal of Chemical Physics;11/15/1987, Vol. 87 Issue 10, p6140 

    Molecular dynamics calculations for a free decane with and without bond and bond angle constraints show that the effect of freezing the bond angles is bigger—or of the same order of magnitude as the solvent effect from other decane molecules in a corresponding liquid decane. Especially...

  • Brownian dynamics simulation of uniform comb polymers in two dimensions. Bishop, Marvin; Saltiel, Craig J. // Journal of Chemical Physics;7/15/1992, Vol. 97 Issue 2, p1471 

    Brownian dynamics simulations are made for continuum, excluded, and nonexcluded volume, uniform H-comb polymers in two dimensions. It is found that the end-to-end distance and radius of gyration of the inner branches are expanded more than those of the outer branches for excluded volume combs....


Read the Article


Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics