Bond fluctuation model

Last updated

The BFM (bond fluctuation model or bond fluctuation method) is a lattice model for simulating the conformation and dynamics of polymer systems. There are two versions of the BFM used: The earlier version was first introduced by I. Carmesin and Kurt Kremer in 1988, [1] and the later version by J. Scott Shaffer in 1994. [2] Conversion between models is possible. [3]

Contents

Model

Carmesin and Kremer version

In this model the monomers are represented by cubes on a regular cubic lattice with each cube occupying eight lattice positions. Each lattice position can only be occupied by one monomer in order to model excluded volume. The monomers are connected by a bond vector, which is taken from a set of typically 108 allowed vectors. There are different definitions for this vector set. One example for a bond vector set is made up from the six base vectors below using permutation and sign variation of the three vector components in each direction:

The resulting bond lengths are and .

The combination of bond vector set and monomer shape in this model ensures that polymer chains cannot cross each other, without explicit test of the local topology.

The basic movement of a monomer cube takes place along the lattice axes

so that each of the possible bond vectors can be realized. [4]

Shaffer's version

As in the case of the Carmesin-Kremer BFM, the Shaffer BFM is also constructed on a simple-cubic lattice. However, the lattice points, or vertices of each cube are the sites that can be occupied by a monomer. Each lattice point can be occupied by one monomer only. Successive monomers along a polymer backbone are connected by bond vectors. The allowed bond vectors must be one of: (a) A cube edge (b) A face diagonal or (c) A solid diagonal. The resulting bond lengths are . In addition to the bond length constraint, polymers should not be allowed to cross. This is done most efficiently by the use of a secondary lattice which is twice as fine as the original lattice. The secondary lattice tracks the midpoints of the bonds in the system, and forbids the overlap of bond midpoints. This effectively leads to disallowing polymers from crossing each other.

Monte Carlo step

In both versions of the BFM, a single attempt to move one monomer consists of the following steps which are standard for Monte Carlo methods:

  1. Select a monomer m and a direction randomly
  2. Check list of conditions (see below)
  3. If all conditions are fulfilled, perform move

The conditions to perform a move can be subdivided into mandatory and optional ones.

Mandatory conditions for CarmesinKremer BFM

  1. Four lattice sites next to monomer m in the direction d are empty.
  2. The move does not lead to bonds that are not contained in the bond vector set.

Mandatory conditions for Shaffer BFM

  1. The lattice site to which the chosen monomer is going to be moved is empty.
  2. The move does not lead to bonds that are not contained in the bond vector set.
  3. The move does not lead to overlapping of bond midpoints.

Optional conditions

If the move leads to an energetic difference for example due to an electric field or an adsorbing force to the walls. In this case a Metropolis algorithm is applied: The Metropolis rate which is defined as

is compared to a random number r from the interval [0, 1). If the Metropolis rate is smaller than r the move is rejected, otherwise it is accepted.

The number of Monte Carlo steps of the total system is defined as:

Notes

  1. Carmesin, I.; Kremer, Kurt (1988). "The bond fluctuation method: a new effective algorithm for the dynamics of polymers in all spatial dimensions". Macromolecules. 21 (9): 2819–2823. Bibcode:1988MaMol..21.2819C. doi:10.1021/ma00187a030. ISSN   0024-9297.
  2. Shaffer, J. Scott (1994). "Effects of chain topology on polymer dynamics: Bulk melts". The Journal of Chemical Physics. 101 (5): 4205–4213. Bibcode:1994JChPh.101.4205S. doi:10.1063/1.467470. ISSN   0021-9606.
  3. Subramanian, Gopinath; Shanbhag, Sachin (2008). "On the relationship between two popular lattice models for polymer melts". The Journal of Chemical Physics. 129 (14): 144904. Bibcode:2008JChPh.129n4904S. doi:10.1063/1.2992047. ISSN   0021-9606. PMID   19045165.
  4. Deutsch, H. P.; Binder, K. (1991). "Interdiffusion and self-diffusion in polymer mixtures: A Monte Carlo study". The Journal of Chemical Physics. 94 (3): 2294. Bibcode:1991JChPh..94.2294D. doi:10.1063/1.459901. ISSN   0021-9606.

Related Research Articles

Tetrahedron Polyhedron with 4 faces

In geometry, a tetrahedron, also known as a triangular pyramid, is a polyhedron composed of four triangular faces, six straight edges, and four vertex corners. The tetrahedron is the simplest of all the ordinary convex polyhedra and the only one that has fewer than 5 faces.

Moment of inertia Scalar measure of the rotational inertia with respect to a fixed axis of rotation

The moment of inertia, otherwise known as the mass moment of inertia, angular mass, or most accurately, rotational inertia, of a rigid body is a quantity that determines the torque needed for a desired angular acceleration about a rotational axis, akin to how mass determines the force needed for a desired acceleration. It depends on the body's mass distribution and the axis chosen, with larger moments requiring more torque to change the body's rate of rotation.

Polymer physics is the field of physics that studies polymers, their fluctuations, mechanical properties, as well as the kinetics of reactions involving degradation and polymerisation of polymers and monomers respectively.

Translation (geometry)

In Euclidean geometry, a translation is a geometric transformation that moves every point of a figure or a space by the same distance in a given direction. A translation can also be interpreted as the addition of a constant vector to every point, or as shifting the origin of the coordinate system. In a Euclidean space, any translation is an isometry.

Root system Geometric arrangements of points, foundational to Lie theory

In mathematics, a root system is a configuration of vectors in a Euclidean space satisfying certain geometrical properties. The concept is fundamental in the theory of Lie groups and Lie algebras, especially the classification and representation theory of semisimple Lie algebras. Since Lie groups and Lie algebras have become important in many parts of mathematics during the twentieth century, the apparently special nature of root systems belies the number of areas in which they are applied. Further, the classification scheme for root systems, by Dynkin diagrams, occurs in parts of mathematics with no overt connection to Lie theory. Finally, root systems are important for their own sake, as in spectral graph theory.

Reciprocal lattice Fourier transform of real-space lattices, important in solid-state physics

In physics, the reciprocal lattice represents the Fourier transform of another lattice. In normal usage, the initial lattice is usually a periodic spatial function in real-space and is also known as the direct lattice. While the direct lattice exists in real-space and is what one would commonly understand as a physical lattice, the reciprocal lattice exists in reciprocal space. The reciprocal lattice of a reciprocal lattice is equivalent to the original direct lattice, because the defining equations are symmetrical with respect to the vectors in real and reciprocal space. Mathematically, direct and reciprocal lattice vectors represent covariant and contravariant vectors, respectively.

The Lenstra–Lenstra–Lovász (LLL) lattice basis reduction algorithm is a polynomial time lattice reduction algorithm invented by Arjen Lenstra, Hendrik Lenstra and László Lovász in 1982. Given a basis with n-dimensional integer coordinates, for a lattice L with , the LLL algorithm calculates an LLL-reduced lattice basis in time

An ideal chain is the simplest model to describe polymers, such as nucleic acids and proteins. It only assumes a polymer as a random walk and neglects any kind of interactions among monomers. Although it is simple, its generality gives insight about the physics of polymers.

Flory–Huggins solution theory Lattice model of polymer solutions

Flory–Huggins solution theory is a lattice model of the thermodynamics of polymer solutions which takes account of the great dissimilarity in molecular sizes in adapting the usual expression for the entropy of mixing. The result is an equation for the Gibbs free energy change for mixing a polymer with a solvent. Although it makes simplifying assumptions, it generates useful results for interpreting experiments.

In solid-state physics, the tight-binding model is an approach to the calculation of electronic band structure using an approximate set of wave functions based upon superposition of wave functions for isolated atoms located at each atomic site. The method is closely related to the LCAO method used in chemistry. Tight-binding models are applied to a wide variety of solids. The model gives good qualitative results in many cases and can be combined with other models that give better results where the tight-binding model fails. Though the tight-binding model is a one-electron model, the model also provides a basis for more advanced calculations like the calculation of surface states and application to various kinds of many-body problem and quasiparticle calculations.

Cauchy stress tensor Representation of mechanical stress at every point within a deformed 3D object

In continuum mechanics, the Cauchy stress tensor, true stress tensor, or simply called the stress tensor is a second order tensor named after Augustin-Louis Cauchy. The tensor consists of nine components that completely define the state of stress at a point inside a material in the deformed state, placement, or configuration. The tensor relates a unit-length direction vector n to the traction vector T(n) across an imaginary surface perpendicular to n:

Linear dynamical systems are dynamical systems whose evaluation functions are linear. While dynamical systems, in general, do not have closed-form solutions, linear dynamical systems can be solved exactly, and they have a rich set of mathematical properties. Linear systems can also be used to understand the qualitative behavior of general dynamical systems, by calculating the equilibrium points of the system and approximating it as a linear system around each such point.

In condensed matter physics and crystallography, the static structure factor is a mathematical description of how a material scatters incident radiation. The structure factor is a critical tool in the interpretation of scattering patterns obtained in X-ray, electron and neutron diffraction experiments.

In statistical mechanics, the transfer-matrix method is a mathematical technique which is used to write the partition function into a simpler form. It was introduced in 1941 by Hans Kramers and Gregory Wannier. In many one dimensional lattice models, the partition function is first written as an n-fold summation over each possible microstate, and also contains an additional summation of each component's contribution to the energy of the system within each microstate.

In computational chemistry, a constraint algorithm is a method for satisfying the Newtonian motion of a rigid body which consists of mass points. A restraint algorithm is used to ensure that the distance between mass points is maintained. The general steps involved are: (i) choose novel unconstrained coordinates, (ii) introduce explicit constraint forces, (iii) minimize constraint forces implicitly by the technique of Lagrange multipliers or projection methods.

A polymer field theory is a statistical field theory describing the statistical behavior of a neutral or charged polymer system. It can be derived by transforming the partition function from its standard many-dimensional integral representation over the particle degrees of freedom in a functional integral representation over an auxiliary field function, using either the Hubbard–Stratonovich transformation or the delta-functional transformation. Computer simulations based on polymer field theories have been shown to deliver useful results, for example to calculate the structures and properties of polymer solutions, polymer melts and thermoplastics.

The transmission-line matrix (TLM) method is a space and time discretising method for computation of electromagnetic fields. It is based on the analogy between the electromagnetic field and a mesh of transmission lines. The TLM method allows the computation of complex three-dimensional electromagnetic structures and has proven to be one of the most powerful time-domain methods along with the finite difference time domain (FDTD) method.

In pure and applied mathematics, particularly quantum mechanics and computer graphics and their applications, a spherical basis is the basis used to express spherical tensors. The spherical basis closely relates to the description of angular momentum in quantum mechanics and spherical harmonic functions.

In pure and applied mathematics, quantum mechanics and computer graphics, a tensor operator generalizes the notion of operators which are scalars and vectors. A special class of these are spherical tensor operators which apply the notion of the spherical basis and spherical harmonics. The spherical basis closely relates to the description of angular momentum in quantum mechanics and spherical harmonic functions. The coordinate-free generalization of a tensor operator is known as a representation operator.

Path integrals in polymer science

A polymer is a macromolecule, composed of many similar or identical repeated subunits. Polymers are common in, but not limited to, organic media. They range from familiar synthetic plastics to natural biopolymers such as DNA and proteins. Their unique elongated molecular structure produces unique physical properties, including toughness, viscoelasticity, and a tendency to form glasses and semicrystalline structures. The modern concept of polymers as covalently bonded macromolecular structures was proposed in 1920 by Hermann Staudinger. One sub-field in the study of polymers is polymer physics. As a part of soft matter studies, Polymer physics concerns itself with the study of mechanical properties and focuses on the perspective of condensed matter physics.