Phase Behavior of Elemental Aluminum Using Monte Carlo Simulations
|
|
- Steven Hoover
- 6 years ago
- Views:
Transcription
1 J. Phys. Chem. B 2006, 110, Phase Behavior of Elemental Aluminum Using Monte Carlo Simulations Divesh Bhatt, Nathan E. Schultz, Ahren W. Jasper, J. Ilja Siepmann,* and Donald G. Truhlar* Department of Chemistry and Supercomputing Institute, UniVersity of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota ReceiVed: July 3, 2006; In Final Form: September 15, 2006 Monte Carlo simulations are presented for two models of aluminum: an embedded-atom model and an explicit many-body model. Vapor/liquid coexistence curves are determined using Gibbs ensemble Monte Carlo simulations. The normal boiling points predicted by both models are somewhat higher (by about 10%) than the experimental value. Isothermal constant-stress simulations are used to simulate solid Al from 300 K to the triple point. The solid structures are at least metastable in the face-centered cubic configuration, and the specific heat is determined to be lower than the experimental value. The melting point for the embeddedatom model determined via thermodynamic integration along a pseudo-supercritical path is approximately 20% higher than the experimental value. Introduction Aluminum is of technological importance not only as a lightweight, rust-resistant structural material, but also as an ingredient for high-energy fuels and, potentially, as a hydrogen storage device. 1 In many applications (e.g., the controlled growth of Al nanoparticles), precise knowledge of the thermodynamic properties of Al, such as the saturated vapor pressure over a large range of temperatures up to the critical point, is pivotal. The high-temperature thermodynamic properties of metals in general, and of Al in particular, are not known, and the present article reports computer simulations performed on analytical potential models to determine such properties. In this context, it becomes imperative that these analytical potential energy functions be validated before applying them for prediction of experimentally difficult-to-determine properties. Density-functional theory (DFT) 2 is used quite extensively to predict potential energy functions, and the PBE0 functional 3 (also called PBEh) has been shown to provide accurate potential energy functions for Al clusters. 4 Recently, we and co-workers 5 presented analytical potential energy functions that were validated against DFT results for Al clusters and nanoparticles. At 0 K, an accurate potential function and a calculation of the zero-point energy suffice to give a reasonably complete and accurate thermodynamic description of a cluster, nanoparticle, or solid, but at finite temperatures, the entropy is important in determining thermodynamic properties, and statistical mechanical methods must be employed. 6 However, DFT often is prohibitively costly for use in finite-temperature statistical mechanical simulations, 7 and analytic potential models become useful. Accordingly, phase equilibrium calculations on bulk systems with known thermodynamic properties provide an additional and important way to further validate the analytic potential energy functions, which then can be used to calculate bulk properties under conditions in which they are not wellknown experimentally, 8 as well to calculate finite-temperature properties of large clusters and nanoparticles that have been recalcitrant to experimental size-selected measurement. * To whom correspondence should be addressed. (J.I.S.) siepmann@chem.umn.edu; (D.G.T) truhlar@umn.edu. Thus, the main aim of this manuscript is to test the applicability of two analytical potential energy functions for Al, previously validated against DFT results at 0 K, against known experimental results such as the normal boiling point and the melting point. Previously, we calculated the vapor-liquid coexistence properties of two embedded-atom potential energy functions and highlighted the sensitivity of the phase diagram to the force field parametrization. 8 While an embedded-atom potential fitted to solid-state data 9 yields an unsatisfactory description of the vapor-liquid coexistence curve (VLCC), another potential function fitted to clusters and nanoparticles of various sizes 10 gives an accurate description of the VLCC at lower temperatures (where experimental data are available) and allows for the prediction of the critical point of Al. 8 This manuscript is organized as follows. First, the potential models that are used for aluminum are described. Then, the details of the simulation methodologies used to simulate various thermodynamic properties of Al are given. This is followed by results and detailed discussions of the vapor-liquid equilibria, simulations of solid Al structures, the solid-vapor equilibria, and the melting point. The final section summarizes the key conclusions. Potential Models for Aluminum We explore two different potential models for Al. The first is an example of an embedded-atom (EAM) model. 11 The total energy of a system of N atoms interacting via an EAM potential is given by 5,9-11 N U ) i F(F i ) + i>j φ(r ij ) (1) and where r ij is the distance between two atoms i and j, and the first term in eq 1 (the embedding energy) is model dependent. For the EAM model employed by Mei and Davenport, 9 this function is given by /jp064169w CCC: $ American Chemical Society Published on Web 11/29/2006
2 26136 J. Phys. Chem. B, Vol. 110, No. 51, 2006 Bhatt et al. F(F) ) -E c[ 1 - R β ( ln F F e)]( F + F e)r/β φ 0 s m exp[-( m - 1)γ] m)1 [ 1 + ( m - 1)δ - m δ β ( ln F F e)]( F F e) mγ/β (2) with and where x ) (r - r n )/(r c - r n ). The second term in eq 1 is a pairwise interaction given by The parameters used in eqs 1-4 are given in Table 1. Manybody interactions appear in this potential in the densitydependent term, F(F), although all terms are functions of the set of pairwise additive distances. Thus, the potential is a function of interpair distances but is not pairwise additive. This form of the EAM potential was used by Mei and Davenport 9 with different parameters than those given in Table 1 to calculate the melting point of Al as 800 ( 9 K. Recently, Jasper et. al. 5 have reparametrized the EAM potential for a large data set of Al cluster and nanoparticle energies. Table 1 represents the result of this reparametrization, which is called NP-B in ref 5, but is simply called EAM here. The second Al potential that is considered in this work is an explicit many-body potential (EMB), also parametrized in ref 5, where it is called NP-A. The form of this potential is given by 5 where the explicit pair potentials, u 2 and u 2, are given by and where Y ij ) r ij - r e and Y ij ) r ij - r e. Explicit multi-body terms appear through the screening function, f ij S, and the coordination number function,f ij CN. The screening function is where F i ) (r ij )f(r ij ) (3) j(*i) 5 c l f(r ij ) )F e 12( r l 0 r ij) (4) l)0 {1 r e r n (r) ) (1 - x) 3 (1 + 3x + 6x 2 ) r n < r < r c (5) 0 r g r c φ(r) )-φ 0 [ 1 + δ ( r r 0-1)] exp [ -γ ( r r 0-1)] (6) U ) u 2 (r ij ) - u 2 (r ij )[f S ij + f CN ij ] (7) i>j i>j u 2 (r ij ) )-D e (1 + a 1 Y ij + a 2 Y ij 2 + a 3 Y ij 3 ) exp(-a 1 Y ij ) (8) u 2 (r ij ) )-D e (1 + a 1 Y ij + a 2 Y 2 ij + a 3 Y ij3 ) exp(-a 1 Y ij ) (9) f ij S ) tan h(χ ij S ) (10) The coordination number function is where χ S ij ) κ 1 exp[-κ 2 (r ik + r jk ) κ 3 κ /r 3 ij ] (11) k*i,j f ij CN ) d(1 - G ij ) (12) 1 1 G ij ) 1 + ( g i - f g (r ij ) γ g 0 ) 1 + ( g j - f g (r ij ) g 0 ) g i is the effective coordination number of atom i defined by and the weighting function, f g is {exp( f g (r ij ) ) γ 1 + γ 1 γ 2 r ij - γ 2) if r ij < γ 2 (15) 0 if r ij g γ 2 The parameters 5 of eqs 5-12 are given in Table 2. We note that the parameters for u 2 were adjusted to fit data for Al 2, and a large data set including Al N clusters and nanoparticles with N ) was then used to optimize the parameters for u 2, f S, and f CN. 5 The many-body terms were designed to vanish at large separations, giving the EMB potential the correct two-body limit. Thus, the EMB potential reproduces the dimer binding energy very well with an error of only 0.01 ev. 5 In contrast, the EAM potential, which was fully optimized over the entire data set including the dimer data, overbinds Al 2 by 0.38 ev. 5 Nevertheless, the fitting errors measured over the entire data set are similar for the EAM and EMB models (0.05 ev/atom and 0.03 ev/atom, respectively). Simulation Methods Vapor-Liquid Coexistence. The vapor-liquid coexistence curve is determined using Gibbs-ensemble Monte Carlo (GEMC) simulations, 12,13 employing translations, aggregation-volumebias Monte Carlo 14 moves to sample clustering in the vapor phase, volume exchanges, and configuration-bias Monte Carlo swaps. 15,16 For the EAM potential, 350 Al atoms with periodic boundary conditions are employed. The size of the vapor box TABLE 1: Parameters for the EAM Potential constant value unit E c ev φ ev r Å R β γ δ c(1) c(2) c(3) c(4) c(5) cc(6) s(1) s(2) s(3) r n 1.75r 0 Å r c 1.95r 0 Å γ (13) g i ) f g (r ik ) (14) k*i
3 Behavior of Al Using Monte Carlo Simulations J. Phys. Chem. B, Vol. 110, No. 51, TABLE 2: Parameters for the EMB Potential constant value unit constant value unit D e ev r e Å a Å -1 a Å -2 a Å -3 a u b u2 6.5 Å κ κ κ a u γ d g g γ Å D e ev r e Å a Å -1 a Å -2 a Å -3 a u2b b u2b 6.5 Å is adjusted such that about 50 atoms are in the vapor phase. Because a liquid structure is used as the initial configuration, these GEMC can be extended to temperatures below the triple point (i.e., into a region in which the liquid phase is only metastable compared to the solid phase). 17 Averages are collected for cycles after allowing equilibration cycles. Each cycle consists of N moves (where N is the total number of atoms in the simulation boxes of the Gibbs ensemble). Evaluation of the EAM potential scales as N 2, whereas that for the EMB potential scales as N 3 ; thus a smaller system size was used for the EMB calculations. For the EMB model, the simulation includes 150 atoms of which about 20 atoms are in the vapor phase. For the EMB potential, equilibration cycles are followed by production cycles during which the averages are collected. For the EAM potential, additional GEMC simulations were performed at the two highest temperatures for a system of 800 Al atoms. As noted in our previous work, 8 the effect of increasing the box size on the vapor-liquid coexistence properties calculated with the EAM potential was negligible; because of its more expensive nature, such a study was not conducted for the EMB potential. Periodic Solid Structures. Isothermal constant-stress Monte Carlo simulations 18 are used to explore the face-centered cubic (FCC) solid for both potentials. In these simulations, Monte Carlo moves allow for sampling of the cell parameters (i.e., the three lengths and three angles that describe the periodic simulation cell) in accordance to a constant external stress, as first suggested by Parrinello and Rahman for molecular dynamics simulations. 19,20 The initial starting structure is FCC, which has the lowest lattice energy, and the noncubic primitive cell was replicated to yield simulation boxes of a suitable size. For the EAM potential, 392 atoms are in the simulation cell and cycles are used, whereas 252 atoms are simulated for cycles for the EMB potential. Solid-Vapor Coexistence. To determine the solid/vapor coexistence for Al, we chose the starting structures to be the equilibrated constant-stress structures generated above. The solid-slab GEMC method of Chen et al. 17 was employed. Thus, a vapor space (of the thickness of the solid slab) was added to the exposed 111 surface of a solid slab on either side (thus, tripling the length of the box containing the solid slab). Because of the addition of the vapor space, the exchange of atoms between the solid structure and the second vapor box of the Gibbs-ensemble is greatly facilitated by allowing the exchange to occur from the surface of the solid slab. 17 The size of the vapor box is chosen such that about 10-15% of the atoms are in the vapor box. Melting Point. To determine the melting point, thermodynamic integration 21 along a pseudo-supercritical path is performed. In this method, 22,23 the average volumes of the two bulk phases (V liq and V sol, in the case of melting) are determined using isothermal-isobaric simulations at the melting pressure of interest (usually 1 atm) and a temperature at which both phases are (meta)stable. Subsequently, the Helmholtz free energy difference between the two phases at their average volumes is calculated via thermodynamic integration along a pseudosupercritical path. This path is divided into three distinct stages. Only a brief description is given here, and further details are available elsewhere. 22,23 In stage A, the liquid is transformed into a weakly attracting fluid at a system volume V liq. Accordingly, a series of canonical ensemble simulations are performed in which the strength of intermolecular interactions is gradually scaled down. In stage B, the intermolecular interactions are kept at their scaled-down values, and an external potential that acts at the lattice sites of an ordered solid is gradually turned on. The form of this external potential is Gaussian (ae -br2,asin the original references, with a ) 3.45 ev and b ) 0.76 Å 2 ; the values of a and b are somewhat arbitrary, but result in gradual ordering of the atoms into a lattice as the external potential is turned on). Additionally, the volume is also changed from V liq to V sol in this stage. In the final stage C, the external potential is gradually turned off, and the full intermolecular potential is gradually restored. Thus, at the end of stage C, a regular bulk solid phase is obtained. The Gibbs free energy difference between the liquid and the solid phase at the particular state point (given by T and p) is obtained as G(T, p) G liq - G sol ) F ex A + F ex B + F ex C + F id + p(v liq - V sol ) (16) where F i ex is the difference in excess Helmholtz free energies at the beginning and the end of stage i calculated from thermodynamic integration, and F id is the ideal part that is calculated analytically (and equals RT ln(f s /F l ), where F s is the density of the solid at the given T and p, and F l is the liquid density under the same conditions). A system size of 256 atoms with periodic boundary conditions was used. A cubic simulation box was used for the simulations of the liquid and solid phases with the latter having an FCC structure. For each stage, thermodynamic integration is performed using points (with a higher point density in the region where the integrand is rapidly varying). Once the Gibbs free energy difference is calculated at a given T and p, itis evaluated at other T using multiple-histogram reweighting 24 for each bulk phase. Results and Discussion Vapor Liquid Equilibria. The vapor-liquid coexistence curve (VLCC) for the EAM potential is obtained for a temperature range from 1100 to 5250 K. At lower temperatures, the particle exchange moves are very inefficient. For the EMB potential, the VLCC is investigated between 1200 and 5500 K.
4 26138 J. Phys. Chem. B, Vol. 110, No. 51, 2006 Bhatt et al. Figure 1. Vapor-liquid coexistence curves for the EAM (diamonds) and EMB (circles) potentials. The corresponding closed symbols represent the respective critical points. Figure 2. Clausius-Clapeyron plot for the EAM and the EMB potentials. The dashed line through each set is a least-squares fit. Figure 1 shows the temperature-density phase diagram for liquid-vapor coexistence with diamonds for the EAM potential and circles for the EMB potential. The open symbols represent the GEMC results, whereas the filled ones are the critical points obtained by the method discussed next. As noted in our previous work, 8 two different methods are used to obtain the critical properties. In the first method, the critical temperature is obtained from the scaling law 21 F l -F v (T c - T) β (17) where F l and F v are the liquid and vapor coexistence densities, respectively, T c is the critical temperature, and the scaling coefficient β is the critical exponent 25 with a value of Additionally, the critical density F c is obtained from the law of rectilinear diameters 21 F l +F v )F 2 c + A(T - T c ) (18) where A is a slope obtained by fitting. The second method uses additional terms in the scaling and rectilinear laws to account for deviations from the principle of corresponding states. 30 As in our previous work, 8 simulated temperatures above the normal boiling point were used in the fits to obtain the critical properties. The average value of the critical temperature for the EAM potential from the two methods is 6299 ( 48 K, and the critical density is 707 ( 60 kg/m 3. Critical properties for the EMB potential are obtained in the manner described above as T c ) 7075 ( 45 K and F c ) 538 ( 8 kg/m 3. As compared to the VLCC of the EAM potential, the VLCC of the EMB potential shows a lower saturated liquid density at low and intermediate temperatures, but the curves cross at about 5500 K. Moreover, the vapor coexistence densities are higher for the EAM potential. This results in a higher T c for the EMB potential. The normal (i.e., 1 atm) boiling temperature, T b, is usually obtained from the simulations using the Clausius-Clapeyron equation 31,32 dlnp v d (1/T) )- H lv R (19) where p v is the vapor pressure () F v RT), and H lv is the heat of vaporization. However, the form of the EAM potential does not permit a straightforward calculation of the pressure via the virial route, 33 because although the force on a given atom can be computed, this force cannot be decomposed simply into pairwise additive terms. To overcome this problem, an alternate method 26 that utilizes the thermodynamic definition of pressure Figure 3. Instantaneous orientational order parameter as a function of MC cycles for bulk structures at three different temperatures obtained from constant-stress simulations at 1 atm. All simulations were started from the FCC structure. is used, p ) p id - U (20) V where p id is the ideal-gas contribution to the pressure, U is the potential energy of the system, and the averaging is done in a system with fixed volume (such as a canonical ensemble). Accordingly, additional NVT simulations are performed at the average GEMC vapor densities, and the pressure is calculated using eq 20. For a given configuration in the canonical ensemble simulations, U/ V is calculated by performing small test volume changes. 26 The resulting plot of p v as a function of 1/T for the temperatures closest to the normal boiling point is shown in Figure 2 for both the potentials. The experimental data is depicted by the solid line. The temperature corresponding to p v ) 1 atm gives T b ) 2993 ( 8 K. Compared to the experimental normal boiling point of 2792 K, 27 the EAM potential overestimates T b for Al by about 7%. For the EMB potential, a similar procedure yields T b to be 3097 ( 9 K. Thus, compared to the more expensive EMB potential, the EAM potential gives a slightly better estimate of the normal boiling point. However, a reliable value of the experimental T c is unavailable 8 to make similar comparison of the two Al potentials for a wider range of the vapor-liquid coexistence. Constant-Stress Simulations. Constant-stress simulations for the FCC solid are performed at 1 atm and several different temperatures for the EAM potential. Although the experimental melting point of Al is K, 27 constant-stress simulations were carried out for temperatures ranging from 900 to 1300 K. Figure 3 depicts the orientational order parameter, Q 6, 28,29 as a
5 Behavior of Al Using Monte Carlo Simulations J. Phys. Chem. B, Vol. 110, No. 51, Figure 5. Stacking of the atoms in three adjacent 111 planes at 900 K obtained from the constant-stress simulations for bulk periodic systems. The open circles represent the first layer, the gray diamonds are in the second layer, and the black triangles are in the third layer. Figure 4. Bulk density profiles along the direction perpendicular to the 111 plane obtained from constant-stress simulations for the EAM Al potential at a pressure of 1 atm and three different temperatures. function of the number of Monte Carlo (MC) cycles for 900, 1250, and 1300 K. In calculating Q 6, the neighbors are defined as pairs of atoms that are closer than the first minimum in the radial distribution functions. For a crystalline solid, Q 6 is approximately 0.5, and it approaches a value of 0 in the thermodynamic limit for a homogeneous liquid. 29 From the values of Q 6 in Figure 3, it is clear that the structures remain crystalline at 900 and 1250 K. On the other hand, as the simulation proceeds at 1300 K, the starting structure loses its crystallinity. For comparison, instantaneous density profiles (after cycles) for the EAM potential are shown for 1300, 1250, and 900 K in Figure 4. These density profiles are collected in bins of width 0.25 Å. In each figure, z is the direction perpendicular to the 111 plane. Thus, the spikes in the density profiles at 900 K (Figure 4, panel c) and 1250 K (Figure 4, panel b) indicate a solidlike structure with the atoms on sets of parallel 111 planes. The spacing between the spikes is approximately 2.8 Å, corresponding to the lattice spacing between adjacent Al atoms. Upon increasing the temperature to 1250 K, the density spikes broaden, as expected. However, at 1300 K, the density profile indicates the onset of melting. For 900 K, the layered structure is shown in more detail in Figure 5; the ABC stacking representative of the an FCC lattice can be clearly seen. This FCC stacking is maintained for all the solid structures. It should be pointed out that the preservation of the FCC structure does not imply that the FCC structure is the most stable one. It does indicate, however, that the FCC structure is probably at least metastable. Table 3 shows the densities and enthalpies of the solid structures as functions of temperature. At 300 K, the density of EAM Al is almost identical to the experimental density of Al (at 293 K, the experimental density is g/cm 3 ). Also, because all of the simulations are performed at 1 atm, the specific heat, C p ) ( H/ T) p, can be evaluated. The filled circles of Figure 6 depict the enthalpies as a function of temperature. Given the wide range of T simulated, it is unlikely that C p can be approximated as a constant; this is confirmed by the Figure 6. Enthalpy of the bulk FCC solid as a function of temperature for the EAM potential at 1 atm. A quadratic fit is shown by the dashed line with H s ) T + ( )T 2. TABLE 3: Solid Densities and Enthalpies for Different Temperatures Using the EAM Potential T (k) density (g/cm 3 ) enthalpy (kj/mol) nonlinearity of the enthalpy-temperature plot in Figure 6. However, the enthalpy can be well-approximated by a quadratic fit, and as expected, the specific heat is an increasing function of temperature. From a quadratic fit to the data, the specific heat is estimated to be 15.2 ( 0.4 J/mol K at 300 K and 19.6 ( 0.6 J/mol K at 900 K. Compared to the experimental values 27 of 24.2 and 32.6 J/mol K at 300 and 900 K, respectively, the heat capacity for the EAM potential is smaller by about 40% at both temperatures. For the EMB potential, one constant-stress simulation at 1000 K was performed. As observed for the EAM potential, the FCC structure for the EMB potential is at least metastable at 1000 K and has an enthalpy of kj/mol. This is larger in magnitude by about 1 kj/mol than the enthalpy value for the
6 26140 J. Phys. Chem. B, Vol. 110, No. 51, 2006 Bhatt et al. Figure 7. Instantaneous orientational order parameter as a function of MC cycles for structures at two different temperatures obtained from solid-slab simulations. Figure 8. Stacking of the atoms on on three adjacent 111 planes at 1075 K for the solid-slab simulation. Symbols as in Figure 5. EAM potential at 1000 K. In contrast, the enthalpy of vaporization at 1200 K for the EAM potential is larger in magnitude than that for the EMB potential. Solid-Vapor Equilibria. As outlined earlier, the solid structures generated along the 111 surface are surrounded by vapor, and solid-slab GEMC simulations are performed to calculate the solid/vapor coexistence. At lower T (below 950 K), there are a few to no particle exchange moves between the solid slab and the vapor boxes. The instantaneous value of Q 6 is shown as a function of the number of MC cycles for two different temperatures in Figure 7. At 1075 K, the evolution of Q 6 shows that the system is a crystalline solid throughout the entire length of the simulation. In contrast, the structure at 1100 K melts, resulting in a lower limiting value of Q 6. Because the current system is in solid/vapor equilibrium, the pressure on the solid structure is the vapor pressure of EAM Al at that temperature, in contrast to the above constant-stress simulations at 1 atm. However, for this range of pressure, the solid structure is still FCC at 1075 K, as can be seen in Figure 8 in which the atoms in three adjacent layers are shown by different symbols (similar to Figure 5). The temperature at which the solid slab melts (i.e., 1100 K) can be contrasted with the temperature of 1250 K mentioned above for which the periodic solid structures (without solid/ vapor interface in the simulation box) are at least metastable. The presence of a solid/vapor interface allows for surface Figure 9. Clausius-Clapeyron plot for the solid-vapor coexistence for the EAM potential represented by the closed circles. The dashed line is the best fit to the sublimation pressure, and the open circles are the liquid-vapor data. melting to occur, and the free energy barrier from the solid to the liquid phase is greatly lowered. Figure 9 shows the Clausius-Clapeyron plot for sublimation and vaporization vapor pressures. Because of the very low coexistence vapor density, the vapor pressure for the solidvapor coexistence at a given T is obtained by treating the vapor phase as an ideal gas instead of performing additional NVT simulation at average vapor density (as in the case of vapor/ liquid coexistence). In principle, the triple point can be determined from the intersection of the sublimation and the boiling vapor pressure lines on such a plot. 17 However, the slopes of the Clausius-Clapeyron plot for the solid-vapor coexistence are almost identical to the slope of the liquid-vapor coexistence. This prevents an accurate determination of the triple point for the EAM potential using this method. An additional issue that arises in such solid-slab simulations is whether the slab is thick enough to allow for an interior region with bulk properties because surface relaxation and melting can be observed for slabs. 17 Effects of the slab thickness were not explored in this work because of the very small heat of fusion. Thus, as described in the following section, the melting point is obtained using a different route. Melting Point. Isobaric-isothermal ensemble simulations at 1100 K and 1 atm with the EAM potential yield a box length of Å for the cubic liquid box containing 256 atoms and a length of Å for the cubic FCC solid box with 256 atoms. Using these two volumes for the bulk condensed phases, explicit thermodynamic integration is performed at 1100 K yielding F A ex ) ( 0.02 kj/mol, F B ex ) ( 0.05 kj/mol, and F C ex ) ( 0.04 kj/mol. The error in these values are due to uncertainties in the values of the integrand at each integration point and do not represent the errors due to the discrete nature of the integration itself. Using these values in eq 13, we find that G(1100 K, 1 atm) ) 0.19 ( 0.11 kj/mol. Accordingly, the solid phase is more stable than the liquid phase at 1100 K and 1 atm, and the normal melting point of EAM Al is higher than 1100 K. Multiple histogram reweighting simulations in the isobaricisothermal ensemble for each bulk phase were performed at four different temperatures (1080, 1100, 1120, and 1140 K) and 1 atm. Figure 10 shows the Gibbs free energies (relative to the respective values at 1100 K) for both the solid and the liquid phases. It must be noted here that the reference values for free energies for the two phases are different, and the purpose of Figure 10 is to highlight that the free energy of the solid phase decreases more rapidly than that of the liquid phase. Combining Figure 10 with the value of G(1100 K, 1atm) given above, yields Figure 11, which depicts the difference in Gibbs free
7 Behavior of Al Using Monte Carlo Simulations J. Phys. Chem. B, Vol. 110, No. 51, Figure 10. Gibbs free energy (relative to the Gibbs free energy at 1100 K and 1 atm) for both the FCC solid and the liquid phases of the EAM Al as functions of temperature. Figure 11. Difference in the Gibbs free energies of the liquid and the solid phases as a function of temperature at 1 atm. energies of liquid and the solid phase as a function of T. Figure 11 shows that G(T, 1atm) changes sign at 1122 K. Thus, the normal melting point of EAM Al is 1122 ( 13 K in which the error bar is based on the uncertainty of G(1100 K, 1atm). This is higher than the experimental value of 933 K by approximately 20%. For the more expensive EMB potential, explicit thermodynamic integration is not performed to determine the melting point. However, the similarities of the enthalpies for the solid and liquid phases at 1000 and 1200 K, respectively, obtained for the EMB and EAM potentials (see above) give an indication that the EMB potential would lead to a similar overestimation of the melting point. The effect of pressure on the melting point was studied by a similar procedure for two different pressures. In one case, explicit thermodynamic integration, performed at 1000 K and atm (vapor pressure of the liquid at 1100 K, as obtained from GEMC simulations), yields G to be 1.17 ( 0.12 kj/mol. In combination with multiple histogram reweighting (at 1000, 1040, 1070, 1100, and 1140 K and at the given pressure) the melting point at the above stated pressure is obtained as 1132 ( 15 K. This value is statistically the same as that at 1 atm. In the second case, explicit thermodynamic integration is performed at 1100 K and at atm; G then is calculated to be 0.99 ( 0.10 kj/mol. Multiple histograms at 1050, 1100, 1150, and 1200 K result in a melting point of 1215 ( 13 K. Accordingly, an increase in pressure to values encountered inside the Earth s crust leads to a significant increase in the melting point. Vocadlo and Alfe 33 and Alfe et al. 34 calculated the melting curve of Al for pressures (up to 150 GPa) using ab initio molecular dynamics simulations. They calculated the liquidphase free energy by constructing a thermodynamically reversibleally path to a Lennard-Jones fluid (and subsequently using the Lennard-Jones equation of state of Johnson et al. 35 ). The free energy of the solid was computed by referencing to a harmonic crystal. As in the current work, a significant increase in the melting point was observed with an increase in pressure. 33,34 Further, the slope of the melting curve was calculated to be approximately 80 K/GPa, 33 a value that agrees well with the current work (an increase in the melting point of approximately 100 K with an increase in pressure by 1 GPa). The form of the EAM potential for Al is the same as that of the original Mei-Davenport embedded-atom (MDEAM) potential 9 for Al that is parametrized to bulk solid-state data. Thus, a comparison of the results obtained with EAM potential to those obtained with their original parameters (called the MDEAM potential) is helpful in determining the effect of parameterizing potentials not only to match bulk solid-state data but also to reproduce accurate energies of clusters and nanoparticles of various sizes (from dimer to bulk). Mei and Davenport calculated a value of 800 K for the melting point of the MDEAM potential; this underestimates the experimental value by 14%, whereas the present value of 1122 K is a 21% overestimate. In a previous study, we reported that the MDEAM potential dramatically underestimates both the boiling point and the heat of vaporization. 8 Conclusions The thermodynamic properties of Al are calculated for two nonpairwise-additive potentials. Vapor-liquid equilibria are determined using Gibbs-ensemble Monte Carlo. The more expensive EMB potential shows a lower saturated vapor density than the EAM potential for a wide range of temperatures. Accordingly, the critical temperature of the EMB potential is higher than that of the EAM potential. From the vapor pressures, the normal boiling point is determined for both the potentials. In accord with the higher value of the critical temperature, the normal boiling point of the EMB potential is higher than that for the EAM potential. In turn, both the potentials overestimate the experimental normal boiling point. Constant-stress simulations in the solid phase show that the FCC structure remains stable over a range of temperatures and that the density of Al at ambient conditions using the EAM potential is very close to the experimental value. The specific heat of the solid phase increases with the temperature and is lower than the experimental value. Solid-slab Gibbs-ensemble simulations are performed to determine the solid-vapor coexistence of the EAM potential. For this purpose, the 111 facet of the FCC structure is exposed to the vapor phase. Compared to a periodic solid (without any vapor phase surrounding it) of the constant-stress simulations, surface melting results in a significantly lower melting temperature. However, an accurate determination of the triple point using the vapor pressure vs temperature plot is precluded by almost identical enthalpies of vaporization and of sublimation. Thermodynamic integration along a pseudo-supercritical path gives a value, 1122 K, of the melting point of EAM Al that is approximately 20% higher than the experimental value. This study highlights the efficacy of analytical potentials parametrized to reproduce accurate DFT energies for clusters of various sizes 5 in making useful predictions of bulk thermodynamic phase behavior. While an earlier embedded-atom model by Mei and Davenport 9 yields very inaccurate results for dimers and smaller clusters that results in a significant error for the vapor-liquid coexistence curve due to formation of dimers in the vapor phase, 8 the present EAM potential yields significantly improved binding energies for smaller clusters, a small fraction
8 26142 J. Phys. Chem. B, Vol. 110, No. 51, 2006 Bhatt et al. of aggregates in the saturated vapor phase, and hence a fairly accurate heat of vaporization and boiling point. The more expensive EMB potential yields similar results for the vaporliquid coexistence curve, but its functional form is more physical, and it performs better than EAM for small clusters. Acknowledgment. This work was supported in part by the Defense-University Research Institute in Nanotechnology (DUR- INT) of the U.S. Army Research Office under agreement number DAAD and by the National Science Foundation under grants CTS , CHE , and ITR (VLab). A Research Scholar Award for D.B. from the Minnesota Supercomputing Institute is gratefully acknowledged. References and Notes (1) Zuttel, A.; Wenger, P.; Sudan, P.; Mauron, P.; Orimo, S. Mater. Sci. Eng. 2004, B108, 9. (2) Kohn, W. ReV. Mod. Phys. 1999, 71, (3) Adamo, C.; Barone, V. J. Chem. Phys. 1999, 110, (4) Schultz, N. E.; Staszewska, G.; Staszewski, P.; Truhlar, D. G. J. Phys. Chem. B 2004, 108, (5) Jasper, A. W.; Schultz, N. E.; Truhlar, D. G. J. Phys. Chem. B 2005, 109, (6) Hill, J.-R.; Subramanian, L.; Maiti, A. Molecular Modeling Techniques in Material Science; CRC Press: Boca Raton, Florida, (7) McGrath, M. J.; Siepmann, J. I.; Kuo, I.-F. W.; Mundy, C. J.; VandeVondele, J.; Hutter, J.; Mohamed, F.; Krack, M. J. Phys. Chem. A 2006, 110, 646. (8) Bhatt, D.; Jasper, A. W.; Schultz, N. E.; Siepmann, J. I.; Truhlar, D. G. J. Am. Chem. Soc. 2006, 128, (9) Mei, J.; Davenport, J. W. Phys. ReV. B1992, 46, 21. (10) Jasper, A. W.; Staszewski, P.; Staszewska, G,; Schultz, N. E.; Truhlar, D. G. J. Phys. Chem. B 2004, 108, (11) Daw, M. S.; Baskes, M. I. Phys. ReV. Lett. 1983, 85, (12) Panagiotopoulos, A. Z. Mol. Phys. 1987, 61, 813. (13) Panagiotopoulos, A. Z.; Quirke, N.; Stapleton, M.; Tildesley, D. J. Mol. Phys. 1988, 63, 527. (14) Chen, B.; Siepmann, J. I. J. Phys. Chem. B 2000, 104, (15) Siepmann, J. I.; Frenkel, D. Mol. Phys. 1992, 75, 59. (16) Mooij, G. C. A. M.; Frenkel, D.; Smit, B. J. Phys.: Condens. Matter 1992, 4, L255. (17) Chen, B.; Siepmann, J. I.; Klein, M. L. J. Phys. Chem. B 2001, 105, (18) Najafabadi, R.; Yip, S. Scr. Metall. 1983, 17, (19) Parrinello, M.; Rahman, A. Phys. ReV. Lett. 1980, 45, (20) Parrinello, M.; Rahman, A. J. Appl. Phys. 1981, 52, (21) Frenkel, D.; Smit, B. Understanding Molecular Simulations; Academic Press: New York, (22) Grochola, G. J. Chem. Phys. 2004, 120, (23) Eike, D. M.; Brennecke, J. F.; Maginn, E. J. J. Chem. Phys. 2005, 122, (24) Ferrenberg, A. M.; Swendsen, R. H. Phys. ReV. Lett. 1988, 61, (25) Rowlinson, J. S.; Widom, B. Molecular Theory of Capillarity; Clarendon Press: Oxford, England, (26) Hummer, G.; Jensen, N. G.; Neumann, M. J. Chem. Phys. 1998, (27) CRC Handbook of Chemistry and Physics, 85th ed.; CRC Press: Boca Raton, Florida, 2004; Section 4. (28) Steinhardt, P. J.; Nelson, D. R.; Ronchetti, M. Phys. ReV. B1983, 28, 784. (29) van Duijneveldt, J. S.; Frenkel, D. J. Chem. Phys. 1992, 96, (30) Sengers, J. V.; Level-Sengers, J. M. H. In Progress in Liquid Physics; Croxton, C. A., Ed.; Wiley: Chichester, England, 1978; pp (31) Clausius, R. Poggendorff s Annalen der Physik und Chemie 1850, 79, 368. (32) Smith, J. M.; Van Ness, H. C. Introduction to Chemical Engineering Thermodynamics, international ed.; McGraw-Hill: New York, 1987; Chapter 6. (33) Vocadlo, L.; Alfe, D. Phys. ReV. B2002, 65, (34) Alfe, D.; Vocadlo, L.; Price, G. D.; Gillan, M. J. J. Phys.: Condens. Matter 2004, 16, S973. (35) Johnson, K.; Zollweg, J. A.; Gubbins, K. E. Mol. Phys. 1993, 78, 591.
Melting line of the Lennard-Jones system, infinite size, and full potential
THE JOURNAL OF CHEMICAL PHYSICS 127, 104504 2007 Melting line of the Lennard-Jones system, infinite size, and full potential Ethan A. Mastny a and Juan J. de Pablo b Chemical and Biological Engineering
More informationSupplemental Material for Temperature-sensitive colloidal phase behavior induced by critical Casimir forces
Supplemental Material for Temperature-sensitive colloidal phase behavior induced by critical Casimir forces Minh Triet Dang, 1 Ana Vila Verde, 2 Van Duc Nguyen, 1 Peter G. Bolhuis, 3 and Peter Schall 1
More informationUB association bias algorithm applied to the simulation of hydrogen fluoride
Fluid Phase Equilibria 194 197 (2002) 249 256 UB association bias algorithm applied to the simulation of hydrogen fluoride Scott Wierzchowski, David A. Kofke Department of Chemical Engineering, University
More informationThermodynamics of Three-phase Equilibrium in Lennard Jones System with a Simplified Equation of State
23 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 28 (2014) Thermodynamics of Three-phase Equilibrium in Lennard Jones System with a Simplified Equation of State Yosuke
More informationChapter 12. Insert picture from First page of chapter. Intermolecular Forces and the Physical Properties of Liquids and Solids
Chapter 12 Insert picture from First page of chapter Intermolecular Forces and the Physical Properties of Liquids and Solids Copyright McGraw-Hill 2009 1 12.1 Intermolecular Forces Intermolecular forces
More informationGibbs ensemble simulation of phase equilibrium in the hard core two-yukawa fluid model for the Lennard-Jones fluid
MOLECULAR PHYSICS, 1989, VOL. 68, No. 3, 629-635 Gibbs ensemble simulation of phase equilibrium in the hard core two-yukawa fluid model for the Lennard-Jones fluid by E. N. RUDISILL and P. T. CUMMINGS
More informationPhase transitions of quadrupolar fluids
Phase transitions of quadrupolar fluids Seamus F. O Shea Department of Chemistry, University of Lethbridge, Lethbridge, Alberta, Canada, T1K 3M4 Girija S. Dubey Brookhaven National Laboratory, Upton, New
More informationChapter 7 PHASE EQUILIBRIUM IN A ONE-COMPONENT SYSTEM
Chapter 7 PHASE EQUILIBRIUM IN A ONE-COMPONENT SYSTEM 7.1 INTRODUCTION The intensive thermodynamic properties of a system are temperature, pressure, and the chemical potentials of the various species occurring
More informationCE 530 Molecular Simulation
CE 530 Molecular Simulation Lecture 20 Phase Equilibria David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Thermodynamic Phase Equilibria Certain thermodynamic states
More informationMolecular dynamics simulation of the liquid vapor interface: The Lennard-Jones fluid
Molecular dynamics simulation of the liquid vapor interface: The Lennard-Jones fluid Matthias Mecke and Jochen Winkelmann a) Institut für Physikalische Chemie, Universität Halle-Wittenberg, Geusaer Str.,
More informationPhase Equilibria of binary mixtures by Molecular Simulation and PR-EOS: Methane + Xenon and Xenon + Ethane
International Journal of ChemTech Research CODEN( USA): IJCRGG ISSN : 0974-4290 Vol.5, No.6, pp 2975-2979, Oct-Dec 2013 Phase Equilibria of binary mixtures by Molecular Simulation and PR-EOS: Methane +
More informationChemical Potential of Benzene Fluid from Monte Carlo Simulation with Anisotropic United Atom Model
Chemical Potential of Benzene Fluid from Monte Carlo Simulation with Anisotropic United Atom Model Mahfuzh Huda, 1 Siti Mariyah Ulfa, 1 Lukman Hakim 1 * 1 Department of Chemistry, Faculty of Mathematic
More informationThermodynamic behaviour of mixtures containing CO 2. A molecular simulation study
Thermodynamic behaviour of mixtures containing. A molecular simulation study V. Lachet, C. Nieto-Draghi, B. Creton (IFPEN) Å. Ervik, G. Skaugen, Ø. Wilhelmsen, M. Hammer (SINTEF) Introduction quality issues
More informationOrder-parameter-based Monte Carlo simulation of crystallization
THE JOURNAL OF CHEMICAL PHYSICS 124, 134102 2006 Order-parameter-based Monte Carlo simulation of crystallization Manan Chopra Department of Chemical Engineering, University of Wisconsin, Madison, Wisconsin
More informationCHEMISTRY The Molecular Nature of Matter and Change
CHEMISTRY The Molecular Nature of Matter and Change Third Edition Chapter 12 Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display. Chapter 11 INTERMOLECULAR FORCES
More informationElastic constants and the effect of strain on monovacancy concentration in fcc hard-sphere crystals
PHYSICAL REVIEW B 70, 214113 (2004) Elastic constants and the effect of strain on monovacancy concentration in fcc hard-sphere crystals Sang Kyu Kwak and David A. Kofke Department of Chemical and Biological
More informationVapor liquid phase coexistence curves for Morse fluids
Fluid Phase Equilibria 248 (2006) 1 6 Vapor liquid phase coexistence curves for Morse fluids Jayant K. Singh a,, Jhumpa Adhikari b, Sang Kyu Kwak c a Department of Chemical Engineering, Indian Institute
More informationPressure-enthalpy driven molecular dynamics for thermodynamic property calculation II: applications
Fluid Phase Equilibria 200 (2002) 93 110 Pressure-enthalpy driven molecular dynamics for thermodynamic property calculation II: applications Loukas I. Kioupis, Gaurav Arya, Edward J. Maginn Department
More informationMatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing
MatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing Due Thursday Feb. 21 at 5pm in Durand 110. Evan Reed In this homework,
More informationChem 112 Dr. Kevin Moore
Chem 112 Dr. Kevin Moore Gas Liquid Solid Polar Covalent Bond Partial Separation of Charge Electronegativity: H 2.1 Cl 3.0 H Cl δ + δ - Dipole Moment measure of the net polarity in a molecule Q Q magnitude
More informationMOLECULAR DYNAMICS SIMULATION OF HETEROGENEOUS NUCLEATION OF LIQUID DROPLET ON SOLID SURFACE
MOLECULAR DYNAMICS SIMULATION OF HETEROGENEOUS NUCLEATION OF LIQUID DROPLET ON SOLID SURFACE Tatsuto Kimura* and Shigeo Maruyama** *Department of Mechanical Engineering, The University of Tokyo 7-- Hongo,
More informationThe Second Law of Thermodynamics (Chapter 4)
The Second Law of Thermodynamics (Chapter 4) First Law: Energy of universe is constant: ΔE system = - ΔE surroundings Second Law: New variable, S, entropy. Changes in S, ΔS, tell us which processes made
More informationPhase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany
Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals Phase Equilibria Phase diagrams and classical thermodynamics
More informationPROPERTIES OF PURE SUBSTANCES. Chapter 3. Mehmet Kanoglu. Thermodynamics: An Engineering Approach, 6 th Edition. Yunus A. Cengel, Michael A.
Thermodynamics: An Engineering Approach, 6 th Edition Yunus A. Cengel, Michael A. Boles McGraw-Hill, 2008 Chapter 3 PROPERTIES OF PURE SUBSTANCES Mehmet Kanoglu Copyright The McGraw-Hill Companies, Inc.
More informationChapter 3 PROPERTIES OF PURE SUBSTANCES
Thermodynamics: An Engineering Approach Seventh Edition in SI Units Yunus A. Cengel, Michael A. Boles McGraw-Hill, 2011 Chapter 3 PROPERTIES OF PURE SUBSTANCES Copyright The McGraw-Hill Companies, Inc.
More informationCHAPTER ELEVEN KINETIC MOLECULAR THEORY OF LIQUIDS AND SOLIDS KINETIC MOLECULAR THEORY OF LIQUIDS AND SOLIDS
CHAPTER ELEVEN AND LIQUIDS AND SOLIDS KINETIC MOLECULAR THEORY OF LIQUIDS AND SOLIDS Differences between condensed states and gases? KINETIC MOLECULAR THEORY OF LIQUIDS AND SOLIDS Phase Homogeneous part
More informationTransferable Potentials for Phase Equilibria. 1. United-Atom Description of n-alkanes
J. Phys. Chem. B 1998, 102, 2569-2577 2569 Transferable Potentials for Phase Equilibria. 1. United-Atom Description of n-alkanes Marcus G. Martin and J. Ilja Siepmann* Department of Chemistry, UniVersity
More information2. As gas P increases and/or T is lowered, intermolecular forces become significant, and deviations from ideal gas laws occur (van der Waal equation).
A. Introduction. (Section 11.1) CHAPTER 11: STATES OF MATTER, LIQUIDS AND SOLIDS 1. Gases are easily treated mathematically because molecules behave independently. 2. As gas P increases and/or T is lowered,
More informationThermodynamics I. Properties of Pure Substances
Thermodynamics I Properties of Pure Substances Dr.-Eng. Zayed Al-Hamamre 1 Content Pure substance Phases of a pure substance Phase-change processes of pure substances o Compressed liquid, Saturated liquid,
More informationStates of Matter; Liquids and Solids. Condensation - change of a gas to either the solid or liquid state
States of Matter; Liquids and Solids Phase transitions - a change in substance from one state to another Melting - change from a solid to a liquid state Freezing - change of a liquid to the solid state
More informationLiquids and Solids. Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display.
Liquids and Solids Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display. 1 Gases, Liquids and Solids Gases are compressible fluids. They have no proper volume and proper
More informationAn Extended van der Waals Equation of State Based on Molecular Dynamics Simulation
J. Comput. Chem. Jpn., Vol. 8, o. 3, pp. 97 14 (9) c 9 Society of Computer Chemistry, Japan An Extended van der Waals Equation of State Based on Molecular Dynamics Simulation Yosuke KATAOKA* and Yuri YAMADA
More informationChapter 3 PROPERTIES OF PURE SUBSTANCES
Thermodynamics: An Engineering Approach Seventh Edition Yunus A. Cengel, Michael A. Boles McGraw-Hill, 2011 Chapter 3 PROPERTIES OF PURE SUBSTANCES Copyright The McGraw-Hill Companies, Inc. Permission
More informationOn the Calculation of the Chemical Potential. Using the Particle Deletion Scheme
On the Calculation of the Chemical Potential Using the Particle Deletion Scheme Georgios C. Boulougouris,2, Ioannis G. Economou and Doros. Theodorou,3,* Molecular Modelling of Materials Laboratory, Institute
More informationTitle Super- and subcritical hydration of Thermodynamics of hydration Author(s) Matubayasi, N; Nakahara, M Citation JOURNAL OF CHEMICAL PHYSICS (2000), 8109 Issue Date 2000-05-08 URL http://hdl.handle.net/2433/50350
More informationAn Introduction to Two Phase Molecular Dynamics Simulation
An Introduction to Two Phase Molecular Dynamics Simulation David Keffer Department of Materials Science & Engineering University of Tennessee, Knoxville date begun: April 19, 2016 date last updated: April
More informationChapter 8 Phase Diagram, Relative Stability of Solid, Liquid, and Gas
Chapter 8 Phase Diagram, Relative Stability of Solid, Liquid, and Gas Three states of matter: solid, liquid, gas (plasma) At low T: Solid is most stable. At high T: liquid or gas is most stable. Ex: Most
More informationDIRECT EVALUATION OF VAPOUR-LIQUID EQUILIBRIA BY MOLECULAR DYNAMICS USING GIBBS-DUHEM INTEGRATION
Molecular Simulation, 1996, Vol. 17, pp. 21-39 Reprints available directly from the publisher Photocopying permitted by license only 0 1996 OPA (Overseas Publishers Association) Amsterdam B.V. Published
More informationMonte Carlo simulations of alloy segregation in PtAg octahedral nanoparticles
Monte Carlo simulations of alloy segregation in PtAg octahedral nanoparticles Louis C. Jones 6/8/12 Abstract Simulations were carried out to investigate phase segregation of insoluble alloy nanoparticles
More information= = 10.1 mol. Molar Enthalpies of Vaporization (at Boiling Point) Molar Enthalpy of Vaporization (kj/mol)
Ch 11 (Sections 11.1 11.5) Liquid Phase Volume and Density - Liquid and solid are condensed phases and their volumes are not simple to calculate. - This is different from gases, which have volumes that
More informationCHEM-UA 652: Thermodynamics and Kinetics
1 CHEM-UA 652: Thermodynamics and Kinetics Notes for Lecture 13 I. PHASE DIAGRAMS The different phases of substances are characterized by different ranges of thermodynamic variables in which these phasesarethestablephases.
More informationThermodynamics of solids 5. Unary systems. Kwangheon Park Kyung Hee University Department of Nuclear Engineering
Thermodynamics of solids 5. Unary systems Kwangheon ark Kyung Hee University Department of Nuclear Engineering 5.1. Unary heterogeneous system definition Unary system: one component system. Unary heterogeneous
More informationATOMISTIC MODELING OF DIFFUSION IN ALUMINUM
ATOMISTIC MODELING OF DIFFUSION IN ALUMINUM S. GRABOWSKI, K. KADAU and P. ENTEL Theoretische Physik, Gerhard-Mercator-Universität Duisburg, 47048 Duisburg, Germany (Received...) Abstract We present molecular-dynamics
More informationPhase Diagrams. NC State University
Chemistry 433 Lecture 18 Phase Diagrams NC State University Definition of a phase diagram A phase diagram is a representation of the states of matter, solid, liquid, or gas as a function of temperature
More informationChapter 2 Experimental sources of intermolecular potentials
Chapter 2 Experimental sources of intermolecular potentials 2.1 Overview thermodynamical properties: heat of vaporization (Trouton s rule) crystal structures ionic crystals rare gas solids physico-chemical
More informationCHAPTER 4 Physical Transformations of Pure Substances.
I. Generalities. CHAPTER 4 Physical Transformations of Pure Substances. A. Definitions: 1. A phase of a substance is a form of matter that is uniform throughout in chemical composition and physical state.
More informationA New Intermolecular Potential Model for the n- Alkane Homologous Series
A New Intermolecular Potential Model for the n- Alkane Homologous Series Jeffrey R. Errington and Athanassios Z. Panagiotopoulos School of Chemical Engineering Cornell University Ithaca, NY 14853 and Institute
More informationChapter 12 INTERMOLECULAR FORCES. Covalent Radius and van der Waals Radius. Intraand. Intermolecular Forces. ½ the distance of non-bonded
Chapter 2 INTERMOLECULAR FORCES Intraand Intermolecular Forces Covalent Radius and van der Waals Radius ½ the distance of bonded ½ the distance of non-bonded Dipole Dipole Interactions Covalent and van
More informationIntermolecular Forces and Liquids and Solids. Chapter 11. Copyright The McGraw Hill Companies, Inc. Permission required for
Intermolecular Forces and Liquids and Solids Chapter 11 Copyright The McGraw Hill Companies, Inc. Permission required for 1 A phase is a homogeneous part of the system in contact with other parts of the
More informationarxiv: v1 [cond-mat.soft] 20 Jun 2008
Accurate determination of crystal structures based on averaged local bond order parameters Wolfgang Lechner and Christoph Dellago Faculty of hysics, University of Vienna, Boltzmanngasse, 19 Vienna, Austria
More informationMonte Carlo Study of the Crystalline and Amorphous NaK Alloy
Available online at www.sciencedirect.com ScienceDirect Procedia Computer Science 08C (207) 25 22 International Conference on Computational Science, ICCS 207, 2-4 June 207, Zurich, Switzerland Monte Carlo
More informationPhase Change (State Change): A change in physical form but not the chemical identity of a substance.
CHM 123 Chapter 11 11.1-11.2 Phase change, evaporation, vapor pressure, and boiling point Phase Change (State Change): A change in physical form but not the chemical identity of a substance. Heat (Enthalpy)
More informationNucleation rate (m -3 s -1 ) Radius of water nano droplet (Å) 1e+00 1e-64 1e-128 1e-192 1e-256
Supplementary Figures Nucleation rate (m -3 s -1 ) 1e+00 1e-64 1e-128 1e-192 1e-256 Calculated R in bulk water Calculated R in droplet Modified CNT 20 30 40 50 60 70 Radius of water nano droplet (Å) Supplementary
More informationRustam Z. Khaliullin University of Zürich
Rustam Z. Khaliullin University of Zürich Molecular dynamics (MD) MD is a computational method for simulating time evolution of a collection of interacting atoms by numerically integrating Newton s equation
More informationChapter 11. Intermolecular Forces and Liquids & Solids
Chapter 11 Intermolecular Forces and Liquids & Solids The Kinetic Molecular Theory of Liquids & Solids Gases vs. Liquids & Solids difference is distance between molecules Liquids Molecules close together;
More informationThermodynamic expansion of nucleation free-energy barrier and size of critical nucleus near the vapor-liquid coexistence
JOURNAL OF CHEMICAL PHYSICS VOLUME 110, NUMBER 7 15 FEBRUARY 1999 hermodynamic expansion of nucleation free-energy barrier and size of critical nucleus near the vapor-liquid coexistence Kenichiro Koga
More informationSupporting Online Material (1)
Supporting Online Material The density functional theory (DFT) calculations were carried out using the dacapo code (http://www.fysik.dtu.dk/campos), and the RPBE (1) generalized gradient correction (GGA)
More informationFrom Atoms to Materials: Predictive Theory and Simulations
From Atoms to Materials: Predictive Theory and Simulations Week 3 Lecture 4 Potentials for metals and semiconductors Ale Strachan strachan@purdue.edu School of Materials Engineering & Birck anotechnology
More informationCARBON 2004 Providence, Rhode Island. Adsorption of Flexible n-butane and n-hexane on Graphitized Thermal Carbon Black and in Slit Pores
CARBON Providence, Rhode Island Adsorption of Flexible n-butane and n-hexane on Graphitized Thermal Carbon Black and in Slit Pores D. D. Do* and H. D. Do, University of Queensland, St. Lucia, Qld 7, Australia
More information- intermolecular forces forces that exist between molecules
Chapter 11: Intermolecular Forces, Liquids, and Solids - intermolecular forces forces that exist between molecules 11.1 A Molecular Comparison of Liquids and Solids - gases - average kinetic energy of
More informationChapter 11 Intermolecular Forces, Liquids, and Solids
Chapter 11 Intermolecular Forces, Liquids, and Solids Dissolution of an ionic compound States of Matter The fundamental difference between states of matter is the distance between particles. States of
More informationChapter 11 Intermolecular Forces, Liquids, and Solids
Chemistry, The Central Science, 11th edition Theodore L. Brown, H. Eugene LeMay, Jr., and Bruce E. Bursten Chapter 11 Intermolecular Forces, Liquids, and Solids John D. Bookstaver St. Charles Community
More information2. As gas P increases and/or T is lowered, intermolecular forces become significant, and deviations from ideal gas laws occur (van der Waal equation).
A. Introduction. (Section 11.1) CHAPTER 11: STATES OF MATTER, LIQUIDS AND SOLIDS 1. Gases are easily treated mathematically because molecules behave independently. 2. As gas P increases and/or T is lowered,
More informationIntermolecular Forces and Liquids and Solids
Intermolecular Forces and Liquids and Solids Chapter 11 Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display. A phase is a homogeneous part of the system in contact
More informationSupplementary Information
Supplementary Information Supplementary Figure 1: Electronic Kohn-Sham potential profile of a charged monolayer MoTe 2 calculated using PBE-DFT. Plotted is the averaged electronic Kohn- Sham potential
More informationFor an incompressible β and k = 0, Equations (6.28) and (6.29) become:
Internal Energy and Entropy as Functions of T and V These are general equations relating the internal energy and entropy of homogeneous fluids of constant composition to temperature and volume. Equation
More informationGas-liquid phase separation in oppositely charged colloids: stability and interfacial tension
7 Gas-liquid phase separation in oppositely charged colloids: stability and interfacial tension We study the phase behaviour and the interfacial tension of the screened Coulomb (Yukawa) restricted primitive
More informationProperty of liquid and Phase Diagram for EN 2017
Property of Liquid and Phase Diagram by Assist.Prof.Dr.Choosak Poonsawat choosak@kku.ac.th FB: เคม อ ช ศ กด You should know - Intermolecular forces - Properties of Liquids - Phase Diagram Properties of
More informationMolecular Simulation Study of the Vapor Liquid Interfacial Behavior of a Dimer-forming Associating Fluid
Molecular Simulation, Vol. 30 (6), 15 May 2004, pp. 343 351 Molecular Simulation Study of the Vapor Liquid Interfacial Behavior of a Dimer-forming Associating Fluid JAYANT K. SINGH and DAVID A. KOFKE*
More informationPhysical transformations of pure substances Boiling, freezing, and the conversion of graphite to diamond examples of phase transitions changes of
Physical transformations of pure substances Boiling, freezing, and the conversion of graphite to diamond examples of phase transitions changes of phase without change of chemical composition. In this chapter
More informationChapter 12 Intermolecular Forces of Attraction
Chapter 12 Intermolecular Forces of Attraction Intermolecular Forces Attractive or Repulsive Forces between molecules. Molecule - - - - - - Molecule Intramolecular Forces bonding forces within the molecule.
More informationStates of matter Part 2
Physical Pharmacy Lecture 2 States of matter Part 2 Assistant Lecturer in Pharmaceutics Overview The Liquid State General properties Liquefaction of gases Vapor pressure of liquids Boiling point The Solid
More informationBasic Thermodynamics Module 1
Basic Thermodynamics Module 1 Lecture 9: Thermodynamic Properties of Fluids Thermodynamic Properties of fluids Most useful properties: Properties like pressure, volume and temperature which can be measured
More informationLecture Notes 2: Physical Equilibria Phase Diagrams
Lecture Notes 2: Physical Equilibria Phase Diagrams There are number of graphical means to help to understand the relationships between the different phases of a particular substance. The first thing we
More informationComparison of different mixing rules for prediction of density and residual internal energy of binary and ternary Lennard Jones mixtures
Fluid Phase Equilibria 178 (2001) 87 95 Comparison of different mixing rules for prediction of density and residual internal energy of binary and ternary Lennard Jones mixtures Jian Chen a,, Jian-Guo Mi
More informationChemical Potential. Combining the First and Second Laws for a closed system, Considering (extensive properties)
Chemical Potential Combining the First and Second Laws for a closed system, Considering (extensive properties) du = TdS pdv Hence For an open system, that is, one that can gain or lose mass, U will also
More informationAn EAM potential for the dynamical simulation of Ni-Al alloys
J. At. Mol. Sci. doi: 10.4208/jams.022310.031210a Vol. 1, No. 3, pp. 253-261 August 2010 An EAM potential for the dynamical simulation of Ni-Al alloys Jian-Hua Zhang, Shun-Qing Wu, Yu-Hua Wen, and Zi-Zhong
More informationModeling aluminum using molecular dynamics simulation
Journal of Materials and Environmental Sciences ISSN : 2028-2508 CODEN : JMESCN Copyright 2017, University of Mohammed Premier Oujda Morocco J. Mater. Environ. Sci., 2018, Volume 9, Issue 1, Page 93-99
More informationStudy of Phase Transition in Pure Zirconium using Monte Carlo Simulation
Study of Phase Transition in Pure Zirconium using Monte Carlo Simulation Wathid Assawasunthonnet, Abhinav Jain Department of Physics University of Illinois assawas1@illinois.edu Urbana, IL Abstract NPT
More informationCritical Properties of Isobaric Processes of Lennard-Jones Gases
Critical Properties of Isobaric Processes of Lennard-Jones Gases Akira Matsumoto Department of Material Sciences, College of Integrated Arts Sciences, Osaka Prefecture University, Sakai, Osaka, 599-8531,
More informationChapter 11. Intermolecular Forces and Liquids and Solids. Chemistry, Raymond Chang 10th edition, 2010 McGraw-Hill
Chemistry, Raymond Chang 10th edition, 2010 McGraw-Hill Chapter 11 Intermolecular Forces and Liquids and Solids Ahmad Aqel Ifseisi Assistant Professor of Analytical Chemistry College of Science, Department
More informationMonte Carlo methods for phase equilibria of uids 2 1. Introduction The phase behavior of uids and their mixtures is of central importance to many tech
TOPICAL REVIEW Monte Carlo methods for phase equilibria of uids Athanassios Z. Panagiotopoulos z Institute for Physical Science and Technology and Department of Chemical Engineering, University of Maryland,
More informationSelf-referential Monte Carlo method for calculating the free energy of crystalline solids
Self-referential Monte Carlo method for calculating the free energy of crystalline solids M. B. Sweatman* Department of Chemical and Process Engineering, University of Strathclyde, Glasgow, G1 1XJ, United
More informationCHEMISTRY. Chapter 11 Intermolecular Forces Liquids and Solids
CHEMISTRY The Central Science 8 th Edition Chapter 11 Liquids and Solids Kozet YAPSAKLI States of Matter difference between states of matter is the distance between particles. In the solid and liquid states
More informationIntermolecular Forces and Liquids and Solids
Intermolecular Forces and Liquids and Solids Chapter 11 Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display. 1 A phase is a homogeneous part of the system in contact
More informationMichael W. Mahoney Department of Physics, Yale University, New Haven, Connecticut 06520
JOURNAL OF CHEMICAL PHYSICS VOLUME 115, NUMBER 23 15 DECEMBER 2001 Quantum, intramolecular flexibility, and polarizability effects on the reproduction of the density anomaly of liquid water by simple potential
More informationChapter 3 PROPERTIES OF PURE SUBSTANCES SUMMARY
Chapter 3 PROPERTIES OF PURE SUBSTANCES SUMMARY PURE SUBSTANCE Pure substance: A substance that has a fixed chemical composition throughout. Compressed liquid (sub-cooled liquid): A substance that it is
More informationIdeal Gas Behavior. NC State University
Chemistry 331 Lecture 6 Ideal Gas Behavior NC State University Macroscopic variables P, T Pressure is a force per unit area (P= F/A) The force arises from the change in momentum as particles hit an object
More informationComputer simulation methods (2) Dr. Vania Calandrini
Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing
More informationChapter 5. On-line resource
Chapter 5 The water-air heterogeneous system On-line resource on-line analytical system that portrays the thermodynamic properties of water vapor and many other gases http://webbook.nist.gov/chemistry/fluid/
More informationMonte Carlo strategies for first-principles simulations of elemental systems
Monte Carlo strategies for first-principles simulations of elemental systems Lev Gelb Department of Materials Science and Engineering, University of Texas at Dallas XSEDE12 Lev Gelb (UT Dallas) Monte Carlo
More informationStructural and Morphological Transitions in Gold Nanorods: A Computer Simulation Study
9214 J. Phys. Chem. B 2003, 107, 9214-9219 Structural and Morphological Transitions in Gold Nanorods: A Computer Simulation Study Yanting Wang, and Christoph Dellago*, Department of Physics and Astronomy
More informationLecture 14. Phases of Pure Substances (Ch.5)
ecture 4. hases of ure Substances (Ch.5) Up to now we have dealt almost exclusively with systems consisting of a single phase. In this lecture, we will learn how more complicated, multiphase systems can
More informationSupporting Information: Surface Polarons Reducing Overpotentials in. the Oxygen Evolution Reaction
Supporting Information: Surface Polarons Reducing Overpotentials in the Oxygen Evolution Reaction Patrick Gono Julia Wiktor Francesco Ambrosio and Alfredo Pasquarello Chaire de Simulation à l Echelle Atomique
More informationTitle Theory of solutions in the energy r of the molecular flexibility Author(s) Matubayasi, N; Nakahara, M Citation JOURNAL OF CHEMICAL PHYSICS (2003), 9702 Issue Date 2003-11-08 URL http://hdl.handle.net/2433/50354
More informationAt this point, we've developed the tools and basic concepts necessary to apply
18 Lecture 0 At this point, we've developed the tools and basic concepts necessary to apply thermodynamics to a number of different systems, with the ultimate goal of describing chemically reacting systems.
More informationSupporting Online Material for
www.sciencemag.org/cgi/content/full/310/5753/1480/dc1 Supporting Online Material for Electrowetting in Carbon Nanotubes J. Y. Chen, A. Kutana, C. P. Collier,* K. P. Giapis* *To whom correspondence should
More informationIntermolecular Forces and Liquids and Solids
PowerPoint Lecture Presentation by J. David Robertson University of Missouri Intermolecular Forces and Liquids and Solids Chapter 11 Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationProblem Set 10 Solutions
Chemistry 360 Dr Jean M Standard Problem Set 10 Solutions 1 Sketch (roughly to scale) a phase diagram for molecular oxygen given the following information: the triple point occurs at 543 K and 114 torr;
More informationTemperature C. Heat Added (Joules)
Now let s apply the heat stuff to real-world stuff like phase changes and the energy or cost it takes to carry it out. A heating curve...a plot of temperature of a substance vs heat added to a substance.
More information