Additional Forces

So far, we have considered how gradients drive the flows of particles, and how chemical reactions can affect those flows by serving as sources or sinks. Diffusing particles can also be subjected to directional external forces. Unlike sources and sinks, forces and fluxes are vectorial. 

For example, gravity draws particles toward the bottom of a beaker of water. Forces accelerate particles according to Newton s laws. But for particles in liquids, acceleration is short-lived. The velocity v times a proportionality constant called the friction coefficient equals the applied force /, 

For proteins or colloids, this frictional regime is reached in nanoseconds or less (see Example 18.2), so you are justified in using Equation (18.41) except for the very fastest processes. Substituting Equation (18.41) into Equation (18.1) gives the flux of particles Jap that results from an applied force /, 

If particles are subject to both a concentration gradient and an applied force, then the two fluxes will sum and Fick s First Law Equation (18.2) generalizes to 

The corresponding generaiization of the diffusion equation, called the Smolu-chow ski equation, is given by combining Equations (18.8) and (18.43)  

---

The interest in vesicles as models for cell biomembranes has led to much work on the interactions within and between lipid layers. The primary contributions to vesicle stability and curvature include those familiar to us already, the electrostatic interactions between charged head groups (Chapter V) and the van der Waals interaction between layers (Chapter VI). An additional force due to thermal fluctuations in membranes produces a steric repulsion between membranes known as the Helfrich or undulation interaction. This force has been quantified by Sackmann and co-workers using reflection interference contrast microscopy to monitor vesicles weakly adhering to a solid substrate [78]. Membrane fluctuation forces may influence the interactions between proteins embedded in them [79]. Finally, in balance with these forces, bending elasticity helps determine shape transitions [80], interactions between inclusions [81], aggregation of membrane junctions [82], and unbinding of pinched membranes [83]. Specific interactions between membrane embedded receptors add an additional complication to biomembrane behavior. These have been stud-... 

Unfortunately, discretization methods with large step sizes applied to such problems tend to miss this additional force term [3]. Furthermore, even if the implicit midpoint method is applied to a formulation in local coordinates, similar problems occur [3]. Since the midpoint scheme and its variants (6) and (7) are basically identical in local coordinates, the same problem can be expected for the energy conserving method (6). To demonstrate this, let us consider the following modified model problem [13] ... 

Following the procedure defined by (23) the fourth order SISM for MD simulations written explicitly can be found In ref. [22]. In the fourth order SISM additional steps in the algorithm occur due to additional force evaluations. 

In general, Laiigeviii dynamics sim illation s run much the same as nioleciilar dynamics simulations. There are differences due Lo the presence of additional forces. Most of the earlier discussions (see pages 69-yO an d p. 3 10-327 of this man ual) on simulation parameters and strategies for molecular dyn amics also apply to Lan gevin dynamics exceptions and additional con sideraiion s are noted below. 

Solenoid Valves The electric solenoid valve has tw o output states. Wlien sufficient electric current is supplied to the coil, an internal armature moves against a spring to an extreme position. This motion causes an attached pneumatic or hvdraiilic valve to operate. Wlien current is removed, the spring returns the armature and the attached solenoid valve to the deenergized position. An intermediate pilot stage is sometimes used when additional force is required to operate the main solenoid valve. Generallv, solenoid valves are used to pressurize or vent the actuator casing for on/off control-valve application and safetv shutdown applications. 

As discussed later, the enclosure of an IPB may carry induced currents up to 95% of the current through the main conductors. Accordingly, the enclosure is designed to carry longitudinal parasitic currents up to 90-95% of the rated current of the main busbars. The cross-sectional area of the enclosure is therefore maintained almost equal to and even more than the main conductors to account for the dissipation of heat of the main conductors through the enclosure only, unless an additional forced cooling system is also adopted. The outdoors part of the enclosure exposed to atmospheric conditions is also subjected to solar radiation. Provision must be made to dissipate this additional heat, from the enclosure. 

When the isobars are curved, an additional force, a centrifugal force outward from the center of curvature, enters into the balance of forces. In the case of curvature around low pressure, a balance of forces occurs when the pressure gradient force equals the sum of the coriolis and centrifugal forces (Fig. 17-12) and the wind continues parallel to the isobars. In the... 

The following additional force fields are currently available in serious professional modelling packages. 

Force. The action of one body on another. This action will cause a change in the motion of the first body unless counteracted by an additional force or forces. A force may be produced either by actual contact or remotely (gravitation, electrostatics, magnetism, etc.). Force is a vector quantity. 

Water is flowing at 5 m/s in a 50 mm diameter pipe which incorporates a 90° bend, as shown in Figure 2.8. What is the additional force to which a retaining bracket will be subjected, as a result of the momentum changes in the liquid, if it is arranged symmetrically in the pipe bend ... 

Equilibrium Theory of Fluid Structure. In all the theoretical work reported herein, we assume that the particles Interact with pair additive forces whose pair potentials can be approximated by... 

While nonbonded atom pairs will typically not come within 1A of each other, it is possible for covalently bound pairs, either directly bounds, as in 1-2 pairs, or at the vertices of an angle, as in 1-3 pairs. Accordingly it may be considered desirable to omit the 1-2 and 1-3 dipole-dipole interactions as is commonly performed on additive force fields for the Coulombic and van der Waals terms. However, it has been shown that inclusion of the 1-2 and 1-3 dipole-dipole interactions is required to achieve anistropic molecular polarizabilites when using isotropic atomic polariz-abilites [50], For example, in a Drude model of benzene in which isotropic polarization was included on the carbons only inclusion of the 1-2 and 1-3 dipole-dipole interactions along with the appropriate damping of those interactions allowed for reproduction of the anisotropic molecular polarizability of the molecule [64], Thus, it may be considered desirable to include these short range interactions in a polarizable force field. 

At every joint, a force of 50 N is applied down and 50 N left, simulating gravity and a mild horizontal force. Additional forces are defined by the fitness function in question. 

When the network chains contain ionic groups, there will be additional forces that affect their swelling properties. Translational entropy of counterions, Coulomb interactions, and ion pair multiplets are forces that lead to interesting phenomena in ion-containing gels. These phenomena were studied in detail by Khokhlov and collaborators [74-77]. The free energy of the networks used by this group is... 

As noted previously, rate accelerations imposed by the cycloamyloses may be competitively inhibited by the addition of inert reagents to the reaction medium. The inhibitor, by competing with the substrate for the cycloamylose cavity, effectively removes a fraction of the catalyst from the reaction coordinate. This observation lends additional force to the mechanism illustrated in scheme I. 

In the Langevin dynamics simulations [33], there is an additional force Sr, acting on the zth bead, where r, is its position and... 

The additional term A( — (x)) is needed to enforce the constraint (x) = constant. This corresponds to an additional force equal to —V(A( — (x))). A is chosen such that (x) = . Therefore the force can be expressed more simply as AV . Then, the interpretation is quite natural. In order to enforce the constraint we are applying a force parallel to V which opposes the mechanical force acting on . 

Based on the ionic radii, nine of the alkali halides should not have the sodium chloride structure. However, only three, CsCl, CsBr, and Csl, do not have the sodium chloride structure. This means that the hard sphere approach to ionic arrangement is inadequate. It should be mentioned that it does predict the correct arrangement of ions in the majority of cases. It is a guide, not an infallible rule. One of the factors that is not included is related to the fact that the electron clouds of ions have some ability to be deformed. This electronic polarizability leads to additional forces of the types that were discussed in the previous chapter. Distorting the electron cloud of an anion leads to part of its electron density being drawn toward the cations surrounding it. In essence, there is some sharing of electron density as a result. Thus the bond has become partially covalent. 

Further cost reductions are possible if the pressure sensor can be fastened to the vibrating suds container instead of requiring an additional component Pressure sensors normally measure the force exerting pressure on a defined surface. Depending on the principle of the particular pressure sensor, it will require an accelerating mass which, when accelerated, produces an additional force and the resulting acceleration signal. Three possible examples are described below. 

---