Helmholtz energy simulation

The Helmholtz and Gibbs energies on the other hand involve constant temperature and volume and constant temperature and pressure, respectively. Most experiments are done at constant Tandp, and most simulations at constant Tand V. Thus, we have now defined two functions of great practical use. In a spontaneous process at constant p and T or constant p and V, the Gibbs or Helmholtz energies, respectively, of the system decrease. These are, however, only other measures of the second law and imply that the total entropy of the system and the surroundings increases. 

The exponential in this equation involves the difference of two energies, rather than an energy itself, and as long as this is sufficiently small compared with k T, a typical simulation run is able to provide a good estimate of the difference in Helmholtz energy of A and B using eq. (11.28). 

In this review, we introduce another approach to study the multiscale structures of polymer materials based on a lattice model. We first show the development of a Helmholtz energy model of mixing for polymers based on close-packed lattice model by combining molecular simulation with statistical mechanics. Then, holes are introduced to account for the effect of pressure. Combined with WDA, this model of Helmholtz energy is further applied to develop a new lattice DFT to calculate the adsorption of polymers at solid-liquid interface. Finally, we develop a framework based on the strong segregation limit (SSL) theory to predict the morphologies of micro-phase separation of diblock copolymers confined in curved surfaces. 

For a binary Ising lattice, we introduced a nonrandom factor that was observed from simulation to have a linear relation with composition. The characteristic parameter of the linear relation was found by combining a series expansion and the infinite dilution properties. On this basis, an accurate expression for the Helmholtz energy of mixing... 

In a nanopore confinement, the Helmholtz energy of both symmetrical and asymmetrical concentric cylinder barrel phases can be obtained by taking d > Rcx in Equation (57) and Equation (60)—(64). In order to compare with the above MC simulation quantitatively, the same parameters were selected. In addition, the statistical bond length in SSL (a = 1.29) was calculated by an ensemble-average bond length over all the collected configurations in MC simulation. 

To establish the molecular thermodynamic model for uniform systems based on concepts from statistical mechanics, an effective method by combining statistical mechanics and molecular simulation has been recommended (Hu and Liu, 2006). Here, the role of molecular simulation is not limited to be a standard to test the reliability of models. More directly, a few simulation results are used to determine the analytical form and the corresponding coefficients of the models. It retains the rigor of statistical mechanics, while mathematical difficulties are avoided by using simulation results. The method is characterized by two steps (1) based on a statistical-mechanical derivation, an analytical expression is obtained first. The expression may contain unknown functions or coefficients because of mathematical difficulty or sometimes because of the introduced simplifications. (2) The form of the unknown functions or unknown coefficients is then determined by simulation results. For the adsorption of polymers at interfaces, simulation was used to test the validity of the weighting function of the WDA in DFT. For the meso-structure of a diblock copolymer melt confined in curved surfaces, we found from MC simulation that some more complex structures exist. From the information provided by simulation, these complex structures were approximated as a combination of simple structures. Then, the Helmholtz energy of these complex structures can be calculated by summing those of the different simple structures. 

To follow to actually carry out a TSA simulation a three-dimensional grid, with grid interval of about 0.2 A (5 -106 equispaced points in (132)) is built and the Helmholtz energies at all grid points are computed. Before this can be done in practice, a value for (A2) must be found. Then, local minima and the crest surfaces must be found, using the procedures given in (130,132,165). To study the dynamics of the penetrant molecules on the network of sites a Monte-Carlo procedure is employed, which is presented is some detail in (97). 

Figure 5 shows the Helmholtz energy as a function of distance between LKal4 (Ca center of mass (COM)) and the surface (frozen CIO atom) for each simulation listed in Table 1. Figure 5c shows the minimum peptide/surface distance for the control simulation is approximately 1 nm therefore, any minima in Fig. 5a, b below 1 nm represent binding to defective areas of the SAM surfaces. 

Our discussion so far has considered the calculation of Helmholtz free energies, which a obtained by performing simulations at constant NVT. For proper comparison with expe inental values we usually require the Gibbs free energy, G. Gibbs free energies are obtaini trorn a simulation at constant NPT. 

As noted above, it is very difficult to calculate entropic quantities with any reasonable accmacy within a finite simulation time. It is, however, possible to calculate differences in such quantities. Of special importance is the Gibbs free energy, as it is the natoal thermodynamical quantity under normal experimental conditions (constant temperature and pressme. Table 16.1), but we will illustrate the principle with the Helmholtz free energy instead. As indicated in eq. (16.1) the fundamental problem is the same. There are two commonly used methods for calculating differences in free energy Thermodynamic Perturbation and Thermodynamic Integration. 

---