School of Chemistry, University of Nottingham, University Park, Nottingham, NG7 2RD, United Kingdom
|
|
- Frederick Cunningham
- 5 years ago
- Views:
Transcription
1 pubs.acs.org/jpcb Calculation of Partition Functions and Free Energies of a Binary Mixture Using the Energy Partitioning Method: Application to Carbon Dioxide and Methane Haina Do,* Jonathan D. Hirst, and Richard J. Wheatley School of Cheistry, University of Nottingha, University Park, Nottingha, NG7 2RD, United Kingdo ABSTRACT: It is challenging to copute the partition function (Q) for systes with enorous configurational spaces, such as fluids. Recently, we developed a Monte Carlo technique (an energy partitioning ethod) for coputing Q [J. Che. Phys. 2011, 135, ]. In this paper, we use this approach to copute the partition function of a binary fluid ixture (carbon dioxide + ethane); this allows us to obtain the Helholtz free energy (F) via F = k B T ln Q and the Gibbs free energy (G) via G = F + pv. We then utilize G to obtain the coexisting ole fraction curves. The cheical potential of each species is also obtained. At the vapor liquid equilibriu condition, the cheical potential of ethane significantly increases, while that of carbon dioxide slightly decreases, as the pressure increases along an isother. Since Q is obtained fro the density of states, which is independent of the teperature, equilibriu therodynaic properties at any condition can be obtained by varying the total coposition and volue of the syste. Our ethodology can be adapted to explore the free energies of other binary ixtures in general and of those containing CO 2 in particular. Since the ethod gives access to the free energy and cheical potentials, it will be useful in any other applications. 1. INTRODUCTION Knowing the free energy allows any physical and cheical phenoena to be predicted. The calculation of free energy fro theory has been an active field for any decades, yet reains a challenging proble. Methods such as free energy perturbation theory 1 3 and therodynaic integration 3,4 are frequently used to deterine the free energy difference between two states (states I and II for exaple), and rely on calculations along a path fro state I to state II. These states can be characterized by either different potential energy functions or different structures. Systes where states I and II are far apart, with large and coplex changes in structure, can require a coplicated path between states and can be prohibitively expensive in ters of coputing effort. Thus, it is iportant to develop ethods that can provide the free energy for each state I and II independently; in this case, the free energy difference can be calculated even for significantly different states because the integration path is avoided. Methods such as inia ining, 5 7 hypothetical scanning, 8 11 haronic reference state, 12,13 and nondynaic growth have ade proising progress toward this goal. However, the efficiency of these techniques deteriorates rapidly with the size and the coplexity of the syste. Elegant approaches for calculating free energies via the density of states have been developed in the past few decades, including ubrella sapling, 17 the weighted histogra analysis ethod (WHAM), 18,19 transition atrix, 20 uticanonical, 21 and Wang Landau 22 sapling. The idea behind these techniques is either to perfor ultiple siulations, whose saplings overlap in configurational space and connect the together (reweight) or to perfor the siulations in biased ensebles to achieve broader sapling of particular states that are rarely visited. If the Wang Landau technique is eployed for systes with an infinite energy range, such as fluids, one often has to choose a finite range of energy (cutting off the high-energy range) either via trial and error or by calculation WHAM also uses a finite range of energies. Restricting the density of states to a finite energy range eans that the partition function can only be coputed to within a ultiplicative constant. To circuvent these liitations, we have recently developed an energy partitioning Monte Carlo (MC) technique, 27 which can calculate the absolute density of configurational energy states with no high-energy cutoff, fro which the excess partition function, and hence the excess free energy of a olecular fluid, can be obtained. The absolute free energy can be calculated by adding the ideal gas free energy to the excess free energy. Phase equilibriu properties of odel systes and real fluids are of treendous iportance to scientists and engineers. Since the easureent of phase equilibriu properties is tie consuing and expensive, coputer siulation based on olecular odeling is a proising alternative. Free energy can be used to predict the vapor liquid equilibriu properties and Received: Deceber 16, 2011 Revised: March 15, 2012 XXXX Aerican Cheical Society A
2 Figure 1. Sketch (not to scale) explaining the partitioning of the density of states (a) at the start of the siulation, where there is only one energy subdivision; (b) after placing the first energy boundary; (c) after placing the second energy boundary; and (d) after placing the nth energy boundary. w(e) is the weighting function and int Ω(E) is the integrated noralized density of states. any other therodynaic properties, including critical phenoena and pvt data for single phases. Of particular interest to us are ixtures that contain carbon dioxide (CO 2 ) Systes of CO 2 in ulticoponent ixtures with light paraffins are iportant to the natural gas industry, since it has becoe coon practice to process gas streas with oderate to high levels of CO 2. Supercritical CO 2 is a novel solvent that has attracted uch attention. 34 A property of interest is the unusually high solubility of fluorinated hydrocarbons in supercritical CO 2. Mixtures of hydrofluorocarbons and CO 2 are possible replaceents for ozonedepleting refrigerants. Thus, knowledge of the phase equilibriu properties of these ixtures can help optiize their perforance in industrial processes. 35 In this paper, we calculate the partition function and free energy of the binary ixture CO 2 +CH 4. Mixtures of CO 2 with other sall olecules (and other binary ixtures) can be studied analogously. Using the free energy, the vapor liquid equilibriu properties and the critical points of the ixture are calculated. The cheical potentials of each species are extracted fro the free energy. The critical point is observed fro a plot of the Helholtz free energy versus volue. The paper is organized as follows. In section 2, we describe the theoretical background of the energy partitioning approach. In section 3, we outline the siulation procedure. In section 4, we present results for the partition function, free energy, and vapor liquid equilibriu properties of the CO 2 +CH 4 ixture. We obtain the vapor liquid equilibriu properties fro the free energies, and copare our results with experiental data. We also copare our results with a conventional MC (Gibbs enseble) siulation technique and deonstrate the advantages of our ethod copared to this approach. Finally, in section 5, we give concluding rearks and outline future work. 2. METHOD 2.1. Free Energy and Partition Function. The free energy of a syste in the canonical enseble is given by F = k B T ln Q(N,V,T), where Q(N,V,T) is the partition function of the syste, which is the integral of the Boltzann factor exp( βe) over particle positions (r N ) and oenta (p N ). As B
3 the integration over particle oenta can be solved exactly, the canonical partition function can be expressed as 1 Q( N, V, T) =!Λ exp[ βe( r N )] dr N 3 N N (1) where N is the nuber of particles in the syste. The factor 1/ N! coes fro the fact that particles are indistinguishable, Λ is the theral de Broglie wavelength, E is configurational energy, β is 1/k B T, k B is the Boltzann constant, and T is the teperature of the syste. Unlike other equilibriu therodynaic properties, such as the internal energy, that can be estiated as a tie or enseble average, the Helholtz free energy F is related directly to the volue of the configurational space (as shown in eq 1), and calculating F still reains a challenge in olecular siulation. Equation 1 can be rewritten in the for of reduced coordinates (s = r/l) as N V 1 Q N V T =!Λ β N N (,, ) exp[ E( s )]ds 3N N 0 where L is the length of the cubic siulation box, which contains N olecules, and V is the volue of the box. The ter V N /(N!Λ 3N ) in eq 2 is the translational partition function of an ideal gas and can be calculated analytically. Thus, we seek to copute the excess part of the partition function (Q ex ) 1 N N Q ex = exp[ βe( s )] ds 0 The probability of finding an energy E is proportional to the density of states, Ω(E), where Ω(E) = 0 1 δ(e E(s N )) ds N. Thus, Q ex can be expressed ore conveniently as Q ex = exp( βe) Ω( E) de Ω( E)dE where the integral is over all possible energies of the syste. Equation 4 shows that if the density of states is known, the excess partition function can be calculated fro it. Our goal is to calculate the density of states and hence Q ex Calculation of the Density of States. The ai of the energy partitioning ethod is to divide the energy range recursively into subdivisions (indexed ), such that the E integrated noralized density of states 1 E Ω(E) de/ Ω(E) de is 1/2 1 for the first energy subdivision ( =1, E 1 E E 0, E 0 = ) (Figure 1d), 1/2 2 for the second subdivision ( =2,E 2 E E 1 ), and so on down to 1/2 n for the two lowest-energy subdivisions ( = n, E n E E n 1 and = n +1, E E n ). At the start of the siulation, there is only one energy subdivision (Figure 1a). All MC oves are accepted at this stage (rando sapling), and the sapled energies are saved. After a predeterined nuber of MC oves (usually several thousand), the energy range E is divided: the first energy boundary E 1 is set equal to the edian configurational energy. The choice of the nuber of MC steps used in each division of the energy is explained in section 3. After the first division of the energy (Figure 1b), there are two energy subdivisions with the sae integrated density of states (1/2 and 1/2). We then throw away all the sapled energies and ove on to the second division of the energy. (2) (3) (4) During the second division of the energy (Figure 1b), MC oves are accepted based on a biased weighting function w(e) = 4, where is the energy subdivision into which the configurational energy E falls so the weighting functions for the energy subdivisions 1 (E > E 1 ) and 2 (E E 1 ) are 4 1 and 4 2, respectively. These biased weighting functions are necessary to push the syste toward the low-energy region; the energy subdivision 2 is visited four ties as often on average as the energy subdivision 1. After the sae predeterined nuber of MC oves as above, the second energy boundary E 2 is set equal to the edian configurational energy found in the lowestenergy subdivision E E 1. This produces three energy subdivisions with the integrated density of states equal to 1/2, 1/4, and 1/4 for the highest energy, iddle, and lowest energy, respectively (Figure 1c). Then the sapled energies are discarded once again and the third division of the energy starts. During the third division of the energy (Figure 1c), the weighting for each energy subdivision is 4 1,4 2, and 4 3 for the highest-energy, iddle, and lowest-energy subdivision respectively, which eans that the lowest-energy subdivision (subdivision 3) is visited four and eight ties as often on average as the subdivisions 2 and 1 respectively. At the end of the MC procedure, the third energy boundary E 3 is set equal to the edian configurational energy found in the lowest-energy subdivision E E 2. There are now four energy subdivisions with the integrated density of states equal to 1/2, 1/4, 1/8, and 1/8 for the highest-energy to the lowest-energy subdivisions respectively. The procedure continues iteratively. In general, for the nth division of the density of states (to produce energy boundary E n, n 2), MC sapling is perfored with a weighting function w(e) =4 for the previously calculated subdivisions, 1 n (subdivision 1 being the highest energy). The weighting function is essential to speed up the siulations; it ensures that about 2/3 or ore of the configurations of the syste fall into the current lowest-energy subdivision. 27 The probability of accepting a ove fro an old state with configurational energy E old that falls into subdivision old,toa new state with configurational energy E new that falls into subdivision new, is given by we ( new) P(old new) = in 1, = in 1, we ( old) 4 new 4 old (5) At the end of the predeterined nuber of MC oves, the energy boundary E n is set equal to the edian configurational energy found in the lowest-energy subdivision E E n 1 and all sapled energies are discarded. The siulation is repeated for the next division of the energy until it is terinated. The density of states Ω(E) strongly increases with energy (except at very high energy), exp( βe) strongly decreases with energy, and the function Ω(E) exp( βe), which is integrated in eq 4 to get Q, has a axiu at an average energy. Since this average energy is not known, the nuber of energy boundaries is not fixed in advance. It is found during the siulation using the criterion that when the integrated function [ [Ω(E) exp( βe) de] ] for the current lowest-energy (th) subdivision is uch saller than its axiu value over all energy subdivisions, [ [Ω(E) exp( βe) de] ax ], the siulation can be terinated, as nothing would be gained by further subdividing the density of states. A stopping criterion of C
4 [ [Ω(E) exp( βe) de] ]=10 9 [ [Ω(E) exp( βe) de] ax ] is used. This is teperature dependent, and we use the lowest teperature of interest, in order to cover all the iportant parts of the energy range. Once the siulation is finished, the excess partition function is obtained fro the noralized integrated density of states as (cf. eq 4) n+ E = 1 1 ( 1 Ω( E)d E)exp( β E ) E Q ex = n+ E = 1 1 ( 1 Ω( E)d E) E = n+ = exp( β E ) + = 1 n 1 2 (6) where E =(E + E 1 )/2. For the first energy subdivision ( = 1), E is set equal to E 1, and for the last energy subdivision, E is set equal to E n. Also, in the last energy subdivision, 2 is replaced by 2 n, as the noralized integrated density of states of the two lowest-energy subdivisions are the sae and equal to 2 n. Errors produced in this application by introducing E are negligible. However, ore care would be needed if very few energy subdivisions were found Interolecular Potentials. The quality of the results of a coputer siulation depends on the potential odels describing the interactions between the olecules of the studied substances. Much effort has been devoted to the developent of an accurate potential for CO 2. In our siulations, the rigid fixed-point charge eleentary physics odel (EPM) is eployed, due to its widespread use. 36 In general, one would not expect force fields to predict properties accurately at high pressure. However, the EPM odel perfors reasonably well under oderately high pressure (<200 bar). 37,38 To enable cobination with the CO 2 odel, potentials for CH 4 with the sae functional for are required. The transferable potentials for phase equilibriu with explicit hydrogen atos (TraPPE-EH) odel 39 is used for CH 4 : the CH 4 olecule has five Lennard-Jones interaction sites, which are located at the carbon ato and the centers of the CH bonds. For Lennard-Jones interactions between unlike atos, the Lorentz Berthelot cobining rules are used. 3. SIMULATION DETAILS A series of siulations are perfored at different total ole fractions and densities to calculate the excess partition functions of the ixture under these conditions. The results are a set of Q ex at different copositions ranging fro the vapor phase to the liquid phase (passing through the two-phase region). The Gibbs free energy is obtained directly fro Q ex (G = k B T ln Q ex + pv) and utilized to obtain the coexisting ole fraction curves at constant pressure. More details about this are given in section 4. At the start of each siulation, 300 olecules with the desired total ole fraction are inserted randoly into a cubic periodic box of fixed volue MC steps are eployed for each partition of the energy. Each MC step involves either a translational ove or a rotational ove of a single olecule with the sae probability. At the beginning of the siulation, the axiu displaceent of a translational ove is set to L/4 and that of a rotational ove is set to 180. The acceptance rates of both oves are 100% during the first division of the energy (as iplied by the weighting function). They then decrease as the nuber of energy subdivisions,, increases. Once becoes large enough for the acceptance rates to drop to 30%, the axiu displaceents are adjusted autoatically after each energy division to keep acceptance rates of about 30% throughout the rest of the siulation. A spherical cutoff of half of the length of the siulation box is used to truncate the potential energy. Thus, a long-range correction for the dispersion r 6 ter (tail correction) 3 is used. The r 12 ter decays rapidly with distance, so a correction for this ter is unnecessary. The CO 2 and CH 4 olecules are nonpolar; the quadrupole quadrupole interactions are sufficiently short range that a long-range correction for the electrostatic energy is not needed. This has been exained and confired in our previous study on this syste. 31 The efficiency and accuracy of the siulations depend on the choice of the nuber of MC steps used in each partition of the energy. Fewer MC steps give faster calculations but lower accuracy. Figure 2 shows the effect on the accuracy of the Figure 2. Effect on the partitioning ethod of the nuber of MC steps used in each partitioning of the energy (50% CO 2 and 50% CH 4 at nuber density = Å 3 ). (a) Errors in selecting the energy boundary in the high-energy region. (b) Errors in selecting the energy boundary in the low-energy region. (c) Errors in the excess partition function (at K). The total nuber of energy subdivisions under these conditions is about The standard deviations are deterined by perforing a set of 15 different siulations. D
5 ethod of the nuber of MC steps used to deterine each energy boundary. In the early state of the partitioning process (the first few hundred partitions), or fewer MC steps could be used (Figure 2a), but as the nuber of energy subdivisions grows, at least MC steps are needed, as the siulations take longer to equilibrate (Figure 2b). If fewer than MC steps are eployed in the low-energy region, a systeatic error occurs. The excess partition function is ost accurately calculated with at least MC steps (Figure 2c). In the rest of this work, MC steps are used. To speed up the calculations, fewer MC steps could be used in the highenergy region than in the low-energy region, but we have not investigated this possibility further. In the first iteration of the siulation, all MC steps are in the lowest-energy subdivision, which spans the entire energy range (Figure 1a). In the second iteration, the weighting between the energy subdivision 1 ( = 1) and the energy subdivision 2 ( = 2) is 1:4 and, therefore, we expect 1/5 of the total MC steps ( 9000) to be in the high-energy subdivision and 4/5 of the total MC steps ( ) in the low-energy subdivision. The energy boundary E 2 is drawn fro the MC steps in the low-energy subdivision. In the third iteration, the ratio between the weighting of the energy subdivisions 1, 2, and 3 (lowest-energy) is 1:4:16 and the integrated density of states is 2:1:1. Thus, we expect to collect 2/22 ( 4091), 4/22 ( 8182), and 16/22 ( ) of the total MC steps in the energy subdivisions 1, 2, and 3, respectively. The energy boundary E 3 is the edian configurational energy found in the MC steps in the lowestenergy subdivision. As the nuber of energy subdivisions increases, about 2/3 ( ) MC steps are found in the current lowest-energy subdivision, as dictated by the weighting function, and the energy boundary E n is set equal to the edian configurational energy found in the lowest-energy subdivision. 4. RESULTS AND DISCUSSION The calculated density of states is used to obtain the partition function, free energy, vapor liquid equilibriu properties, and cheical potential of the binary ixture CO 2 +CH 4. After the partitioning process, the entire energy range (for a fixed coposition and volue) has been discretized into n + 1 energy subdivisions (see Figure 1). The excess partition function is obtained using eq 6, and the excess Helholtz free energy (F ex ) is calculated using F ex = k B T ln Q ex. The density dependence of the ideal gas free energy (F id ) is calculated using F id = k B T(x CO2 ln ρ CO2 + x CH4 ln ρ CH4 ), where x is the ole fraction and ρ is the nuber density. The densitydependent part of the absolute free energy is the su of the excess part and the ideal gas part (F = F ex + F id ). The Gibbs free energy is G = F + pv, where p is the pressure and V is the volue of the syste. Figure 3 shows the Gibbs free energy per particle versus volue per particle for CO 2 +CH 4 ixtures at K, bar and K, bar. This shows the preference for each phase of the syste under different conditions. For exaple, at K and bar (Figure 3a), the liquid phase is favored (a lower iniu in the free energy) when there is no CH 4 present. When the ole fraction of CH 4 is 0.2, two inia are observed, but the liquid phase is still favored. As the total ole fraction of CH 4 increases, the vapor phase becoes favored. The liquid phase is less favored and disappears when the ole fraction of CH 4 reaches about 0.6; i.e., the critical coposition is reached for this teperature. Figure 3. Gibbs free energy per particle versus volue per particle of CO 2 +CH 4 at different copositions: (a) K and bar and (b) K and bar. The standard errors are deterined by perforing a set of 15 different siulations for each volue. In ost cases, the errors in G/N are between 0.4% and 0.7%. For the sake of clarity, the plots at the copositions of CH 4 equal to 0.0, 0.2, 0.4, 0.6, and 0.8 are offset by 5, 3.5, 2.5, 1.5, and 0.5 kj/ol, respectively. Siilar observations can also be ade at other conditions, for exaple, at K and bar (Figure 3b). Using the inia in Figure 3, a plot of the Gibbs free energy versus the coposition for both phases is constructed, fro which the vapor liquid equilibriu copositions are deterined (Figure 4). The coexisting copositions are deterined by constructing a double coon tangent connecting both phases. The cheical potential for each species at any coposition can be calculated by taking G and dg/dx fro Figure 4 and solving the siultaneous equations μ + μ = x1 1 x2 2 G =μ μ d G /dx1 1 2 (7) for μ 1 and μ 2, where μ is the cheical potential and can be either vapor or liquid phase. The phase diagra (pressure versus coposition) of the binary ixture CO 2 +CH 4 is extracted fro Figure 4 and plotted in Figure 5, with results fro experient 40 and fro the Gibbs enseble technique. 31 The Gibbs enseble technique is teperature-dependent. Its efficiency depends on the nuber of state points that are required to be calculated. E
6 Figure 4. Gibbs free energy per particle versus CH 4 coposition for CO 2 +CH 4 : (a) K and bar and (b) K and bar. The standard deviations are deterined and shown as error bars by perforing a set of 15 different siulations for each volue. The energy partitioning ethod, on the other hand, is teperature-independent and does not suffer fro this shortcoing. The vapor liquid equilibriu properties coputed using both siulation techniques agree well to within the statistical uncertainties of the siulations apart fro a couple of points at low pressure. This discrepancy is, perhaps, due to the errors in fitting the curves and placing the double tangent lines. Reasonable agreeent between experiental data and siulations is achieved for the liquid phase at both teperatures. However, the siulations slightly overestiate the solubility of CH 4 in the vapor phase, which is due to the liitation of the force fields. This has been discussed in our previous work. 31 Due to the finite size effect, the Gibbs enseble technique cannot be used near the critical point. Using the energy partitioning ethod, we observe the occurrence of the critical point by constructing a plot of the Helholtz free energy versus the volue of the syste. Figure 6 shows such plots for ixtures of CO 2 +CH 4 at and K. The critical point occurs when the second and third derivatives d n F/dV n are both zero, and the critical pressure is df/dv at the sae point. The occurrence of the critical point is observed visually. The critical coposition at K is x CH and the critical pressure is approxiately 70 bar. The critical coposition at K is x CH4 0.55, and the critical pressure is approxiately 81 bar. These results were obtained using ole fractions of CH 4 separated by 0.02 (not shown in Figure 6 for clarity), and the uncertainties in the critical coposition and pressure are estiated to be 0.02 and 5 bar, respectively. The critical pressure can also be extrapolated fro Figure 5. Phase diagra of the CO 2 +CH 4 syste at (a) K and (b) K. Energy partitioning ethod (squares with error bars) versus experiental data (pluses) and Gibbs enseble technique (crosses with error bars). Stars indicate the critical points calculated fro the energy partitioning ethod. the siulated data by plotting p versus x vap x liq. At the critical point, x vap x liq is equal to zero. These results agree with those given above. The apparent critical points of a finite syste depend on the syste size, V. These paraeters obey a scaling law behavior with V To investigate the sensitivity of the calculated critical points with the syste size used, we have calculated the critical points (critical pressures and copositions) for systes with double and half of the size of the studied syste. We found that the apparent critical pressures and CH 4 ole fractions of sall systes are slightly higher than those of bigger systes. However, these differences are still within the uncertainties of the estiation of the critical points. Figure 7 shows the cheical potentials of each species in the ixture CO 2 +CH 4 at the vapor liquid equilibriu condition, calculated using eq 7. At both teperatures, the cheical potential of CH 4 increases steeply while that of CO 2 decreases gradually as the pressure increases. The gradient of the cheical potential of CH 4 decreases as the pressure increases. This indicates that at the vapor liquid equilibriu condition the cheical potential of CH 4 is ore sensitive to pressure than that of CO 2 and that the solubility of CH 4 in CO 2 increases with the pressure. The rate of change of the cheical potentials with respect to the pressure of a given species (species 1, for exaple) at twophase equilibriu can be expressed in ters of the ole F
7 Figure 6. Helholtz free energy per particle versus volue per particle of the ixture CO 2 +CH 4 at different copositions at (a) K and (b) K. The standard errors are deterined by perforing a set of 15 different siulations for each volue. In ost cases, the errors are between 0.4% and 0.7%. For the sake of clarity, the plots in (a) at the copositions of CH 4 equal to 0.5, 0.54, 0.58, 0.62, and 0.64 are offset by 1.0, 0.8, 0.6, 0.4, and 0.2 kj/ol, respectively. The plots in (b) at the copositions of CH 4 equal to 0.3, 0.4, 0.5, and 0.55 are offset by 0.8, 0.6, 0.4, and 0.2 kj/ol, respectively. fraction of the other species in both phases and the olar volue (v) of each phase as β dμ β 1 x = 2v x2 v dp β x2 x2 (8) Thus, since the ole fractions and the volues of each phase are known (Figures 3 and 4), dμ 1 /dp can be calculated. For exaple, dμ CH4 /dp at K and bar is about 870 c 3 / ol and that calculated using eq 8 gives 890 c 3 /ol. The difference between these two nubers coes fro the uncertainties in placing the tangent line. Equation 8 can also be used to interpret the behavior of the cheical potentials shown in Figure 7. The rate of change of the cheical potential of CH 4 with respect to the pressure is (x liq CO2 v vap x vap CO2 v liq )/ liq (x CO2 x vap CO2 ). Since the olar volue of vapor is larger than that of liquid and the ole fraction of CO 2 in the liquid is always greater than in the vapor, this quantity is always positive. At low pressure, the olar volue of the vapor phase is uch larger than that of the liquid phase, and this results in a steep slope. As the pressure increases, the volue of the vapor phase decreases and also the ole fractions of CO 2 in both phases get closer to each other, which decreases the agnitude of dμ CH4 / dp. This is observed fro Figure 7. A siilar analysis can be done for dμ CH4 /dp, which equals (x liq CH4 v vap x vap CH4 v liq liq )/(x CH4 vap ). Although the volue of the vapor phase is uch greater x CH4 Figure 7. Cheical potentials of each species in the ixture CO 2 + CH 4 at the vapor liquid equilibriu condition versus pressure at (a) K and (b) K. than that of the liquid phase, the ole fraction of CH 4 in the liquid phase is saller than that in the vapor phase. Therefore, x liq CH4 v vap x vap CH4 v liq liq is positive, while x CH4 x vap CH4 is negative. Thus, dμ CH2 /dp is negative, which is observed in Figure 7b. 5. CONCLUSION In this paper, we have calculated the density of states, partition function, free energy, and cheical potential of the binary ixture CO 2 + CH 4 using an energy partitioning MC technique. The Gibbs free energy of the binary ixture is coputed as a function of the volue and the coposition. Knowing the Gibbs free energy allows us to predict the vapor liquid equilibriu properties, which agree quite well with experiental data and siulations using the Gibbs enseble technique. One of the advantages of the energy partitioning ethod is that it allows us to observe the occurrence of the critical coposition by constructing plots of the Helholtz free energy versus the volue. We also obtain the cheical potentials of each species in the ixture directly. These quantities are useful, but difficult to copute using other ethods. At the vapor liquid equilibriu condition, we find that the rate of change of the cheical potentials with respect to pressure of CH 4 increases significantly, whereas that of CO 2 slightly decreases. The cheical potential of CH 4 varies significantly with pressure, while that of CO 2 does not vary uch. This is understandable in ters of the coposition and the volue of both phases. The energy partitioning ethod can be applied to other binary ixtures and can be easily extended to study ternary ixtures and the free energy of transfer between two solvent phases. Standard ethods for easuring free energy differences, G
8 such as therodynaic integration and free energy perturbation, are not directly applicable to calculations of the free energy of transfer. 46 A cobination of the Gibbs enseble ethod, Wido s test particle 47 and the configurational biased MC technique 48,49 have been used to tackle this proble for alkanes. 46,50 Using the energy partitioning ethod, the free energy of transfer could be obtained fro the cheical potential. Since the ethod gives access to the free energy, it will be useful in any other applications. Also, the ethod can be extended to study systes with discrete energy levels, for which there is also a large nuber of applications. AUTHOR INFORMATION Corresponding Author *E-ail: haina.do@nottingha.ac.uk. Notes The authors declare no copeting financial interest. ACKNOWLEDGMENTS We thank the University of Nottingha High Perforance Coputing facility for providing coputer recourses and the Engineering and Physical Sciences Research Council (EPSRC) for funding (Grant No. EP/E06082X). We thank Prof. Martyn Poliakoff and Dr. Jie Ke for useful discussions. H.D. is grateful to the EPSRC for a Ph.D. Plus Fellowship. REFERENCES (1) Henderson, D.; Barker, J. A. Phys. Rev. A 1970, 1, (2) Zwanzig, R. W. J. Che. Phys. 1954, 22, (3) Frenkel, D.; Sit, B. Understanding Molecular Siulation: fro Algoriths to Applications, 2nd ed.; Acadeic Press: San Diego, CA, (4) Hansen, J.; Verlet, L. Phys. Rev. 1969, 184, (5) Head, M. S.; Given, J. A.; Gilson, M. K. J. Phys. Che. A 1997, 101, (6) Krisvov, S. V.; Karplus, M. J. J. Che. Phys. 2002, 117, (7) Evans, D. A.; Wales, D. J. J. Che. Phys. 2003, 119, (8) White, R. P.; Meirovitch, H. J. Che. Phys. 2003, 119, (9) White, R. P.; Meirovitch, H. Proc. Natl. Acad. Sci. U.S.A 2004, 101, (10) Meirovitch, H. Phys. Rev. A 1985, 32, (11) Meirovitch, H. J. Che. Phys. 1999, 111, (12) Tyka, M. D.; Clarke, A. R.; Sessions, R. B. J. Phys. Che. B 2006, 110, (13) Tyka, M. D.; Sessions, R. B.; Clarke, A. R. J. Phys. Che. B 2007, 111, (14) Ytreberg, F. M.; Zuckeran, D. M. J. Che. Phys. 2006, 124, (15) Zhang, X.; Maonov, A. B.; Zuckeran, D. M. J. Coput. Che. 2009, 30, (16) Bhatt, D.; Zuckeran, D. M. J. Phys. Che. 2009, 131, (17) Torrie, G. M.; Valleau, J. P. J. Coput. Phys. 1977, 23, (18) Ferrenberg, A. M.; Swendsen, R. H. Phys. Rev. Lett. 1988, 61, (19) Ferrenberg, A. M.; Swendsen, R. H. Phys. Rev. Lett. 1989, 63, (20) Wang, J. S.; Tay, T. K.; Swendsen, R. H. Phys. Rev. Lett. 1999, 82, (21) Berg, B.; Neuhaus, T. Phys. Lett. B 1991, 267, (22) Wang, F. G.; Landau, D. P. Phys. Rev. Lett. 2001, 86, (23) Ganzenuller, G.; Cap, P. J. J. Che. Phys. 2007, 127, H (24) Shell, M. S.; Debenedetti, P. G.; Panagiotopoulos, A. Z. Phys. Rev. E 2002, 66, (25) Shell, M. S.; Debenedetti, P. G.; Panagiotopoulos, A. Z. J. Che. Phys. 2003, 119, (26) Yan, Q.; Faller, R.; de Pablo, J. J. J. Che. Phys. 2002, 116, (27) Do, H.; Hirst, J. D.; Wheatley, R. J. J. Che. Phys. 2011, 135, (28) Oakley, M. T.; Wheatley, R. J. J. Che. Phys. 2009, 130, (29) Oakley, M. T.; Do, H.; Wheatley, R. J. Fluid Phase Equilib. 2009, 290, (30) Oakley, M. T.; Do, H.; Hirst, J. D.; Wheatley, R. J. J. Che. Phys. 2011, 134, (31) Do, H.; Wheatley, R. J.; Hirst, J. D. J. Phys. Che. B 2010, 114, (32) Do, H.; Wheatley, R. J.; Hirst, J. D. Phys. Che. Che. Phys. 2010, 12, (33) Do, H.; Wheatley, R. J.; Hirst, J. D. Phys. Che. Che. Phys. 2011, 13, (34) Poliakoff, M.; King, P. Nature 2001, 412, (35) Skaroutsos, I.; Hunt, P. A. J. Phys. Che. B 2010, 114, (36) Harris, J. G.; Yung, K. H. J. Phys. Che. 1995, 99, (37) Vorholz, J.; Hariseadis, V. I.; Rupf, B.; Panagiotopoulos, A. Z.; Maurer, G. Fluid Phase Equilib. 2000, 170, (38) Liu, Y.; Panagiotopoulos, A. Z.; Debenedetti, P. G. J. Phys. Che. B 2011, 115, (39) Chen, B.; Siepann, J. I. J. Phys. Che. B 1999, 103, (40) Wei, M. S. W.; Brown, T. S.; Kidnay, A. J.; Sloan, E. D. J. Che. Eng. Data 1995, 40, (41) Chen, J. H.; Fisher, M. E.; Nickel, B. G. Phys. Rev. Lett. 1982, 48, (42) Ferrenberg, A. M.; Landau, D. P. Phys. Rev. B 1991, 44, (43) Wilding, N. B. Phys. Rev. E 1995, 52, (44) Wilding, N. B. Phys. Rev. E 1997, 55, (45) Potoff, J. J.; Panagiotopoulos, A. Z. J. Che. Phys. 1998, 109, (46) Martin, M. G.; Siepann, J. I. J. A. Che. Soc. 1997, 119, (47) Wido, B. J. Che. Phys. 1963, 39, (48) Siepann, J. I.; Frenkel, D. Mol. Phys. 1992, 68, (49) Frenkel, D.; Mooji, G. C. A. M.; Sit, B. J. Phys.: Condens. Matter 1992, 4, (50) Martin, M. G.; Siepann, J. I. Theor. Che. Acc. 1998, 99,
PHY 171. Lecture 14. (February 16, 2012)
PHY 171 Lecture 14 (February 16, 212) In the last lecture, we looked at a quantitative connection between acroscopic and icroscopic quantities by deriving an expression for pressure based on the assuptions
More informationDeveloped Correlations for Prediction of The Enthalpies of Saturated Vapor Liquid Coexisting Phases
Nahrain University, College of Engineering Journal (NUCEJ) Vol.13 No.2, 2010 pp.116-128 Developed Correlations for Prediction of he Enthalpies of Saturated Vapor Liquid Coexisting Phases Mahoud Oar bdullah
More informationDensity and structure of undercooled liquid titanium
Article Condensed Matter Physics March 2012 Vol.57 No.7: 719 723 doi: 10.1007/s11434-011-4945-6 Density and structure of undercooled liquid titaniu WANG HaiPeng, YANG ShangJing & WEI BingBo * Departent
More informationModel Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon
Model Fitting CURM Background Material, Fall 014 Dr. Doreen De Leon 1 Introduction Given a set of data points, we often want to fit a selected odel or type to the data (e.g., we suspect an exponential
More informationClassical systems in equilibrium
35 Classical systes in equilibriu Ideal gas Distinguishable particles Here we assue that every particle can be labeled by an index i... and distinguished fro any other particle by its label if not by any
More informationExtension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels
Extension of CSRSM for the Paraetric Study of the Face Stability of Pressurized Tunnels Guilhe Mollon 1, Daniel Dias 2, and Abdul-Haid Soubra 3, M.ASCE 1 LGCIE, INSA Lyon, Université de Lyon, Doaine scientifique
More informationGAUTENG DEPARTMENT OF EDUCATION SENIOR SECONDARY INTERVENTION PROGRAMME. PHYSICAL SCIENCE Grade 11 SESSION 11 (LEARNER NOTES)
PYSICAL SCIENCE Grade 11 SESSION 11 (LEARNER NOTES) MOLE CONCEPT, STOICIOMETRIC CALCULATIONS Learner Note: The ole concept is carried forward to calculations in the acid and base section, as well as in
More informationarxiv: v1 [cond-mat.stat-mech] 22 Dec 2017
Notes on the Hybrid Monte Carlo Method arxiv:1712.08278v1 [cond-at.stat-ech] 22 Dec 2017 Jerey C. Paler, 1, Air Haji-Akbari, 2, Rakesh S. Singh, 2 Fausto Martelli, 3 Roberto Car, 3 Athanassios Z. Panagiotopoulos,
More informationA method to determine relative stroke detection efficiencies from multiplicity distributions
A ethod to deterine relative stroke detection eiciencies ro ultiplicity distributions Schulz W. and Cuins K. 2. Austrian Lightning Detection and Inoration Syste (ALDIS), Kahlenberger Str.2A, 90 Vienna,
More informationOcean 420 Physical Processes in the Ocean Project 1: Hydrostatic Balance, Advection and Diffusion Answers
Ocean 40 Physical Processes in the Ocean Project 1: Hydrostatic Balance, Advection and Diffusion Answers 1. Hydrostatic Balance a) Set all of the levels on one of the coluns to the lowest possible density.
More informationChapter 4: Hypothesis of Diffusion-Limited Growth
Suary This section derives a useful equation to predict quantu dot size evolution under typical organoetallic synthesis conditions that are used to achieve narrow size distributions. Assuing diffusion-controlled
More informationANALYTICAL INVESTIGATION AND PARAMETRIC STUDY OF LATERAL IMPACT BEHAVIOR OF PRESSURIZED PIPELINES AND INFLUENCE OF INTERNAL PRESSURE
DRAFT Proceedings of the ASME 014 International Mechanical Engineering Congress & Exposition IMECE014 Noveber 14-0, 014, Montreal, Quebec, Canada IMECE014-36371 ANALYTICAL INVESTIGATION AND PARAMETRIC
More informationTHERMAL ENDURANCE OF UNREINFORCED UNSATURATED POLYESTERS AND VINYL ESTER RESINS = (1) ln = COMPOSITES & POLYCON 2009
Aerican Coposites Manufacturers Association January 15-17, 29 Tapa, FL USA Abstract THERMAL ENDURANCE OF UNREINFORCED UNSATURATED POLYESTERS AND VINYL ESTER RESINS by Thore M. Klaveness, Reichhold AS In
More informationThe Thermal Conductivity Theory of Non-uniform Granular Flow and the Mechanism Analysis
Coun. Theor. Phys. Beijing, China) 40 00) pp. 49 498 c International Acadeic Publishers Vol. 40, No. 4, October 5, 00 The Theral Conductivity Theory of Non-unifor Granular Flow and the Mechanis Analysis
More informationExperimental Design For Model Discrimination And Precise Parameter Estimation In WDS Analysis
City University of New York (CUNY) CUNY Acadeic Works International Conference on Hydroinforatics 8-1-2014 Experiental Design For Model Discriination And Precise Paraeter Estiation In WDS Analysis Giovanna
More informationRecommended Reading. Entropy/Second law Thermodynamics
Lecture 7. Entropy and the second law of therodynaics. Recoended Reading Entropy/econd law herodynaics http://en wikipedia http://en.wikipedia.org/wiki/entropy http://2ndlaw.oxy.edu/index.htl. his site
More informationChapter 12. Quantum gases Microcanonical ensemble
Chapter 2 Quantu gases In classical statistical echanics, we evaluated therodynaic relations often for an ideal gas, which approxiates a real gas in the highly diluted liit. An iportant difference between
More informationKinetic Molecular Theory of Ideal Gases
Lecture -3. Kinetic Molecular Theory of Ideal Gases Last Lecture. IGL is a purely epirical law - solely the consequence of experiental obserations Explains the behaior of gases oer a liited range of conditions.
More informationAP Physics Thermodynamics Wrap-up
AP Physics herodynaics Wrap-up Here are your basic equations for therodynaics. here s a bunch of the. 3 his equation converts teperature fro Fahrenheit to Celsius. his is the rate of heat transfer for
More informationSpine Fin Efficiency A Three Sided Pyramidal Fin of Equilateral Triangular Cross-Sectional Area
Proceedings of the 006 WSEAS/IASME International Conference on Heat and Mass Transfer, Miai, Florida, USA, January 18-0, 006 (pp13-18) Spine Fin Efficiency A Three Sided Pyraidal Fin of Equilateral Triangular
More informationln P 1 saturation = T ln P 2 saturation = T
More Tutorial at www.littledubdoctor.co Physical Cheistry Answer each question in the space provided; use back of page if extra space is needed. Answer questions so the grader can READILY understand your
More informationKinetic Theory of Gases: Elementary Ideas
Kinetic Theory of Gases: Eleentary Ideas 17th February 2010 1 Kinetic Theory: A Discussion Based on a Siplified iew of the Motion of Gases 1.1 Pressure: Consul Engel and Reid Ch. 33.1) for a discussion
More informationA Self-Organizing Model for Logical Regression Jerry Farlow 1 University of Maine. (1900 words)
1 A Self-Organizing Model for Logical Regression Jerry Farlow 1 University of Maine (1900 words) Contact: Jerry Farlow Dept of Matheatics Univeristy of Maine Orono, ME 04469 Tel (07) 866-3540 Eail: farlow@ath.uaine.edu
More informationKINETIC THEORY. Contents
KINETIC THEORY This brief paper on inetic theory deals with three topics: the hypotheses on which the theory is founded, the calculation of pressure and absolute teperature of an ideal gas and the principal
More informationStatistical associating fluid theory-dimer of the solid phase of the pearl-necklace model
Indian Journal of Pure & Applied Physics Vol., epteber 00, pp 688-69 tatistical associating fluid theory-dier of the solid phase of the pearl-necklace odel Vivek Kuar ingh & K N Khanna Departent of Physics,
More informationOBJECTIVES INTRODUCTION
M7 Chapter 3 Section 1 OBJECTIVES Suarize data using easures of central tendency, such as the ean, edian, ode, and idrange. Describe data using the easures of variation, such as the range, variance, and
More informationNonmonotonic Networks. a. IRST, I Povo (Trento) Italy, b. Univ. of Trento, Physics Dept., I Povo (Trento) Italy
Storage Capacity and Dynaics of Nononotonic Networks Bruno Crespi a and Ignazio Lazzizzera b a. IRST, I-38050 Povo (Trento) Italy, b. Univ. of Trento, Physics Dept., I-38050 Povo (Trento) Italy INFN Gruppo
More informationCONTINUOUS THERMODYNAMICS FINITE DIFFUSION MODEL FOR MULTICOMPONENT FUEL SPRAY EVAPORATION
CONTINUOUS THERMODYNAMICS FINITE DIFFUSION MODEL FOR MULTICOMPONENT FUEL SPRAY EVAPORATION Dongyao Wang 1 and Chia-fon F. Lee Departent of Theoretical and Applied Mechanics 1 Departent of Mechanical and
More information1 The properties of gases The perfect gas
1 The properties of gases 1A The perfect gas Answers to discussion questions 1A. The partial pressure of a gas in a ixture of gases is the pressure the gas would exert if it occupied alone the sae container
More informationFeature Extraction Techniques
Feature Extraction Techniques Unsupervised Learning II Feature Extraction Unsupervised ethods can also be used to find features which can be useful for categorization. There are unsupervised ethods that
More informationI affirm that I have never given nor received aid on this examination. I understand that cheating in the exam will result in a grade F for the class.
Che340 hysical Cheistry for Biocheists Exa 3 Apr 5, 0 Your Nae _ I affir that I have never given nor received aid on this exaination. I understand that cheating in the exa will result in a grade F for
More informationMeasuring Temperature with a Silicon Diode
Measuring Teperature with a Silicon Diode Due to the high sensitivity, nearly linear response, and easy availability, we will use a 1N4148 diode for the teperature transducer in our easureents 10 Analysis
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 11 Jan 2007
Transport and Helfand oents in the Lennard-Jones fluid. II. Theral conductivity arxiv:cond-at/7125v1 [cond-at.stat-ech] 11 Jan 27 S. Viscardy, J. Servantie, and P. Gaspard Center for Nonlinear Phenoena
More information1.1 Heat and Mass transfer in daily life and process/mechanical engineering Heat transfer in daily life: Heating Cooling Cooking
1. Introduction 1.1 Heat and Mass transfer in daily life and process/echanical engineering Heat transfer in daily life: Heating Cooling Cooking ransfer of heat along a teperature difference fro one syste
More informationKinetic Theory of Gases: Elementary Ideas
Kinetic Theory of Gases: Eleentary Ideas 9th February 011 1 Kinetic Theory: A Discussion Based on a Siplified iew of the Motion of Gases 1.1 Pressure: Consul Engel and Reid Ch. 33.1) for a discussion of
More informationMolecular dynamic simulation of the melting and solidification processes of argon
Journal of Mechanical Science and Technology 23 (200) 1563~1570 Journal of Mechanical Science and Technology www.springerlink.co/content/1738-44x DOI 10.1007/s12206-00-0418-0 Molecular dynaic siulation
More informationThe Thermal Dependence and Urea Concentration Dependence of Rnase A Denaturant Transition
The Theral Dependence and Urea Concentration Dependence of Rnase A Denaturant Transition Bin LI Departent of Physics & Astronoy, University of Pittsburgh, Pittsburgh, PA 15260, U.S.A Feb.20 th, 2001 Abstract:
More information1 (40) Gravitational Systems Two heavy spherical (radius 0.05R) objects are located at fixed positions along
(40) Gravitational Systes Two heavy spherical (radius 0.05) objects are located at fixed positions along 2M 2M 0 an axis in space. The first ass is centered at r = 0 and has a ass of 2M. The second ass
More informationAll Excuses must be taken to 233 Loomis before 4:15, Monday, April 30.
Miscellaneous Notes he end is near don t get behind. All Excuses ust be taken to 233 Loois before 4:15, Monday, April 30. he PHYS 213 final exa ties are * 8-10 AM, Monday, May 7 * 8-10 AM, uesday, May
More informationAnnealing contour Monte Carlo algorithm for structure optimization in an off-lattice protein model
JOURNAL OF CHEMICAL PHYSICS VOLUME 120, NUMBER 14 8 APRIL 2004 Annealing contour Monte Carlo algorith for structure optiization in an off-lattice protein odel Faing Liang a) Departent of Statistics, Texas
More informationIN modern society that various systems have become more
Developent of Reliability Function in -Coponent Standby Redundant Syste with Priority Based on Maxiu Entropy Principle Ryosuke Hirata, Ikuo Arizono, Ryosuke Toohiro, Satoshi Oigawa, and Yasuhiko Takeoto
More informationKeywords: Estimator, Bias, Mean-squared error, normality, generalized Pareto distribution
Testing approxiate norality of an estiator using the estiated MSE and bias with an application to the shape paraeter of the generalized Pareto distribution J. Martin van Zyl Abstract In this work the norality
More informationMolecular dynamics algorithm for multiple time scales: Systems with long range forces
Molecular dynaics algorith for ultiple tie scales: Systes with long range forces Mark E. Tuckeran@ Bruce J. Berne Departent of Cheistry, Colubia University, New York, New York 10027 Glenn J. Martyna Departent
More informationDiffusion time-scale invariance, randomization processes, and memory effects in Lennard-Jones liquids
PHYSICAL REVIEW E 68, 52 23 Diffusion tie-scale invariance, randoization processes, and eory effects in Lennard-Jones liquids Renat M. Yuletyev* and Anatolii V. Mokshin Departent of Physics, Kazan State
More informationChapter 6 1-D Continuous Groups
Chapter 6 1-D Continuous Groups Continuous groups consist of group eleents labelled by one or ore continuous variables, say a 1, a 2,, a r, where each variable has a well- defined range. This chapter explores:
More informationMcMillan Mayer theory for solvent effects in inhomogeneous systems: Calculation of interaction pressure in aqueous electrical double layers
McMillan Mayer theory for solvent effects in inhoogeneous systes: Calculation of interaction pressure in aqueous electrical double layers Roland Kjellander, Alexander P. Lyubartsev, and Stjepan Marčelja
More informationarxiv: v1 [cond-mat.stat-mech] 7 Jul 2009
arxiv:97.114v1 [cond-at.stat-ech] 7 Jul 29 Finding the optiu activation energy in DNA breathing dynaics: A Siulated Annealing approach Pinaki Chaudhury 1, Ralf Metzler 2 and Suan K Banik 3 1 Departent
More informationNumerical Solution of the MRLW Equation Using Finite Difference Method. 1 Introduction
ISSN 1749-3889 print, 1749-3897 online International Journal of Nonlinear Science Vol.1401 No.3,pp.355-361 Nuerical Solution of the MRLW Equation Using Finite Difference Method Pınar Keskin, Dursun Irk
More informationThe Transactional Nature of Quantum Information
The Transactional Nature of Quantu Inforation Subhash Kak Departent of Coputer Science Oklahoa State University Stillwater, OK 7478 ABSTRACT Inforation, in its counications sense, is a transactional property.
More informationOStudy of Real Gas Behavior: Ideality of CO 2 Gas
OStudy of Real Gas Behavior: Ideality of CO Gas Subitted: March, 014 CHEM 457, Section Departent of Cheistry, The Pennsylvania State University, University Park, PA 1680 Jessica Slavejkov Bashayer Aldakkan,
More information13 Harmonic oscillator revisited: Dirac s approach and introduction to Second Quantization
3 Haronic oscillator revisited: Dirac s approach and introduction to Second Quantization. Dirac cae up with a ore elegant way to solve the haronic oscillator proble. We will now study this approach. The
More informationEXPLORING PHASE SPACES OF BIOMOLECULES WITH MONTE CARLO METHODS
Journal of Optoelectronics and Advanced Materials Vol. 7, o. 3, June 2005, p. 1563-1571 EXPLORIG PHASE SPACES OF BIOMOLECULES WITH MOTE CARLO METHODS A. Bu u * ational Institute for Research and Developent
More informationAn analytical relation between relaxation time spectrum and molecular weight distribution
An analytical relation between relaxation tie spectru and olecular weight distribution Wolfgang Thi, Christian Friedrich, a) Michael Marth, and Josef Honerkap b) Freiburger Materialforschungszentru, Stefan-Meier-Straße
More informationKinetic Molecular Theory of. IGL is a purely empirical law - solely the
Lecture -3. Kinetic Molecular Theory of Ideal Gases Last Lecture. IGL is a purely epirical law - solely the consequence of experiental obserations Explains the behaior of gases oer a liited range of conditions.
More informationA Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair
Proceedings of the 6th SEAS International Conference on Siulation, Modelling and Optiization, Lisbon, Portugal, Septeber -4, 006 0 A Siplified Analytical Approach for Efficiency Evaluation of the eaving
More informationNational 5 Summary Notes
North Berwick High School Departent of Physics National 5 Suary Notes Unit 3 Energy National 5 Physics: Electricity and Energy 1 Throughout the Course, appropriate attention should be given to units, prefixes
More informationSolidification of Porous Material under Natural Convection by Three Phases Modeling
Solidification of Porous Material under Natural Convection by Three Phases Modeling Hassan Basirat Tabrizi, Meber, IAENG and F. Sadeghpour Abstract The perforance of natural convective flow over a rectangular
More informationElectronic Supplementary Information Proton Conductivity in Mixed-Conducting BSFZ Perovskite from Thermogravimetric Relaxation
Electronic Suppleentary Material (ESI) for Physical Cheistry Cheical Physics. This journal is the wner Societies 014 Z ig / k Z / k phase / Electronic Suppleentary Inforation Proton Conductivity in Mixed-Conducting
More informationGeneral Physical Chemistry I
General Physical Cheistry I Lecture 12 Aleksey Kocherzhenko Aril 2, 2015" Last tie " Gibbs free energy" In order to analyze the sontaneity of cheical reactions, we need to calculate the entroy changes
More informationOptical Properties of Plasmas of High-Z Elements
Forschungszentru Karlsruhe Techni und Uwelt Wissenschaftlishe Berichte FZK Optical Properties of Plasas of High-Z Eleents V.Tolach 1, G.Miloshevsy 1, H.Würz Project Kernfusion 1 Heat and Mass Transfer
More informationLONG-TERM PREDICTIVE VALUE INTERVAL WITH THE FUZZY TIME SERIES
Journal of Marine Science and Technology, Vol 19, No 5, pp 509-513 (2011) 509 LONG-TERM PREDICTIVE VALUE INTERVAL WITH THE FUZZY TIME SERIES Ming-Tao Chou* Key words: fuzzy tie series, fuzzy forecasting,
More informationCOS 424: Interacting with Data. Written Exercises
COS 424: Interacting with Data Hoework #4 Spring 2007 Regression Due: Wednesday, April 18 Written Exercises See the course website for iportant inforation about collaboration and late policies, as well
More informationMolecular Speeds. Real Gasses. Ideal Gas Law. Reasonable. Why the breakdown? P-V Diagram. Using moles. Using molecules
Kinetic Theory of Gases Connect icroscopic properties (kinetic energy and oentu) of olecules to acroscopic state properties of a gas (teperature and pressure). P v v 3 3 3 But K v and P kt K v kt Teperature
More informationNon-Parametric Non-Line-of-Sight Identification 1
Non-Paraetric Non-Line-of-Sight Identification Sinan Gezici, Hisashi Kobayashi and H. Vincent Poor Departent of Electrical Engineering School of Engineering and Applied Science Princeton University, Princeton,
More informationTopic 5a Introduction to Curve Fitting & Linear Regression
/7/08 Course Instructor Dr. Rayond C. Rup Oice: A 337 Phone: (95) 747 6958 E ail: rcrup@utep.edu opic 5a Introduction to Curve Fitting & Linear Regression EE 4386/530 Coputational ethods in EE Outline
More informationCh 12: Variations on Backpropagation
Ch 2: Variations on Backpropagation The basic backpropagation algorith is too slow for ost practical applications. It ay take days or weeks of coputer tie. We deonstrate why the backpropagation algorith
More informationSF Chemical Kinetics.
SF Cheical Kinetics. Lecture 5. Microscopic theory of cheical reaction inetics. Microscopic theories of cheical reaction inetics. basic ai is to calculate the rate constant for a cheical reaction fro first
More informationSPECTRUM sensing is a core concept of cognitive radio
World Acadey of Science, Engineering and Technology International Journal of Electronics and Counication Engineering Vol:6, o:2, 202 Efficient Detection Using Sequential Probability Ratio Test in Mobile
More informationProbability Distributions
Probability Distributions In Chapter, we ephasized the central role played by probability theory in the solution of pattern recognition probles. We turn now to an exploration of soe particular exaples
More informationI. Concepts and Definitions. I. Concepts and Definitions
F. Properties of a syste (we use the to calculate changes in energy) 1. A property is a characteristic of a syste that can be given a nuerical value without considering the history of the syste. Exaples
More informationOptimizing energy potentials for success in protein tertiary structure prediction Ting-Lan Chiu 1 and Richard A Goldstein 1,2
Research Paper 223 Optiizing energy potentials for success in protein tertiary structure prediction Ting-Lan Chiu 1 and Richard A Goldstein 1,2 Background: Success in solving the protein structure prediction
More informationIs the seismic moment frequency relation universal?
Geophys. J. Int. (2) 142, 193 198 Is the seisic oent frequency relation universal? C. Godano1 and F. Pingue2 1 Dipartiento di Scienze Abientali, Seconda Università di Napoli, via V ivaldi 43, 811 Caserta,
More informationANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER
IEPC 003-0034 ANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER A. Bober, M. Guelan Asher Space Research Institute, Technion-Israel Institute of Technology, 3000 Haifa, Israel
More informationPhysical Chemistry I for Biochemists Chem340. Lecture 32 (4/4/11)
Physical Cheistry I for Biocheists Che340 Lecture 32 (4/4/11) Yoshitaka Ishii Ch8.8-8.12 If you have a note 33, skip printing p2-3. 8.5 he Gibbs-Duhe Equation In Ch 6, we learned dg = -Sd + VdP + i dn
More informationDepartment of Electronic and Optical Engineering, Ordnance Engineering College, Shijiazhuang, , China
6th International Conference on Machinery, Materials, Environent, Biotechnology and Coputer (MMEBC 06) Solving Multi-Sensor Multi-Target Assignent Proble Based on Copositive Cobat Efficiency and QPSO Algorith
More informationCrystallization of Supercooled Liquid Elements Induced by Superclusters Containing Magic Atom Numbers Abstract: Keywords: 1.
Crystallization of Supercooled Liquid Eleents Induced by Superclusters Containing Magic Ato Nubers Robert F. Tournier, CRETA /CNRS, Université Joseph Fourier, B.P. 166, 804 Grenoble cedex 09, France. E-ail:
More informationPhysical Chemistry I for Biochemists Chem340. Lecture 26 (3/14/11)
Physical Cheistry I or Biocheists Che340 Lecture 26 (3/14/11) Yoshitaka Ishii Ch 7.2, 7.4-5, & 7.10 Announceent Exa 2 this Friday. Please be well prepared! HW average 80-85. You will probably have one
More informationSharp Time Data Tradeoffs for Linear Inverse Problems
Sharp Tie Data Tradeoffs for Linear Inverse Probles Saet Oyak Benjain Recht Mahdi Soltanolkotabi January 016 Abstract In this paper we characterize sharp tie-data tradeoffs for optiization probles used
More informationANALYSIS OF THE EFFECT OF THE CHEMICAL SPECIES CONCENTRATIONS ON THE RADIATION HEAT TRANSFER IN PARTICIPATING GASES USING A MONTE CARLO METHODOLOGY
Proceedings of the 11 th Brazilian Congress of Theral Sciences and Engineering -- ECIT 2006 Braz. Soc. of Mechanical Sciences and Engineering -- ABCM, Curitiba, Brazil,- Dec. 5-8, 2006 Paper CIT06-0366
More informationRegular article Maximum radius of convergence perturbation theory: test calculations on Be, Ne, H 2 and HF
Theor Che Acc (2003) 110: 185 189 DOI 10.1007/s00214-003-0473-z Regular article Maxiu radius of convergence perturbation theory: test calculations on Be, Ne, H 2 and HF Kotaro Yokoyaa 1, Haruyuki Nakano
More informationN-Point. DFTs of Two Length-N Real Sequences
Coputation of the DFT of In ost practical applications, sequences of interest are real In such cases, the syetry properties of the DFT given in Table 5. can be exploited to ake the DFT coputations ore
More informationReading from Young & Freedman: For this topic, read the introduction to chapter 25 and sections 25.1 to 25.3 & 25.6.
PHY10 Electricity Topic 6 (Lectures 9 & 10) Electric Current and Resistance n this topic, we will cover: 1) Current in a conductor ) Resistivity 3) Resistance 4) Oh s Law 5) The Drude Model of conduction
More informationSupport Vector Machine Classification of Uncertain and Imbalanced data using Robust Optimization
Recent Researches in Coputer Science Support Vector Machine Classification of Uncertain and Ibalanced data using Robust Optiization RAGHAV PAT, THEODORE B. TRAFALIS, KASH BARKER School of Industrial Engineering
More informationInspection; structural health monitoring; reliability; Bayesian analysis; updating; decision analysis; value of information
Cite as: Straub D. (2014). Value of inforation analysis with structural reliability ethods. Structural Safety, 49: 75-86. Value of Inforation Analysis with Structural Reliability Methods Daniel Straub
More information5.60 Thermodynamics & Kinetics Spring 2008
MIT OpenCourseWare http://ocw.it.edu 5.60 Therodynaics & Kinetics Spring 2008 For inforation about citing these aterials or our Ters of Use, visit: http://ocw.it.edu/ters. 1 Enzye Catalysis Readings: SAB,
More informationReconstruction of the electron density of molecules with single-axis alignment
SLAC-PUB-143 Reconstruction of the electron density of olecules with single-axis alignent Ditri Starodub* a, John C. H. Spence b, Dilano K. Saldin c a Stanford PULSE Institute, SLAC National Accelerator
More informationAnalyzing Simulation Results
Analyzing Siulation Results Dr. John Mellor-Cruey Departent of Coputer Science Rice University johnc@cs.rice.edu COMP 528 Lecture 20 31 March 2005 Topics for Today Model verification Model validation Transient
More informationSupporting information for Self-assembly of multicomponent structures in and out of equilibrium
Supporting inforation for Self-assebly of ulticoponent structures in and out of equilibriu Stephen Whitela 1, Rebecca Schulan 2, Lester Hedges 1 1 Molecular Foundry, Lawrence Berkeley National Laboratory,
More informationModulation of Harmonic Emission Spectra from Intense Laser-Plasma Interactions
Modulation of Haronic Eission Spectra fro Intense Laser-Plasa Interactions T.J.M. Boyd and R. Ondarza-Rovira 2 Centre for Physics, University of Essex, Wivenhoe Park, Colchester CO4 3SQ, U.K. 2 ININ, A.P.
More informationŞtefan ŞTEFĂNESCU * is the minimum global value for the function h (x)
7Applying Nelder Mead s Optiization Algorith APPLYING NELDER MEAD S OPTIMIZATION ALGORITHM FOR MULTIPLE GLOBAL MINIMA Abstract Ştefan ŞTEFĂNESCU * The iterative deterinistic optiization ethod could not
More informationChapter 1 Introduction and Kinetics of Particles
Chapter 1 Introduction and Kinetics of Particles 1.1 Introduction There are two ain approaches in siulating the transport equations (heat, ass, and oentu), continuu and discrete. In continuu approach,
More informationIII.H Zeroth Order Hydrodynamics
III.H Zeroth Order Hydrodynaics As a first approxiation, we shall assue that in local equilibriu, the density f 1 at each point in space can be represented as in eq.iii.56, i.e. f 0 1 p, q, t = n q, t
More informationSynchronization in large directed networks of coupled phase oscillators
CHAOS 16, 015107 2005 Synchronization in large directed networks of coupled phase oscillators Juan G. Restrepo a Institute for Research in Electronics and Applied Physics, University of Maryland, College
More informationInteractive Markov Models of Evolutionary Algorithms
Cleveland State University EngagedScholarship@CSU Electrical Engineering & Coputer Science Faculty Publications Electrical Engineering & Coputer Science Departent 2015 Interactive Markov Models of Evolutionary
More information1 Brownian motion and the Langevin equation
Figure 1: The robust appearance of Robert Brown (1773 1858) 1 Brownian otion and the Langevin equation In 1827, while exaining pollen grains and the spores of osses suspended in water under a icroscope,
More informationSolving initial value problems by residual power series method
Theoretical Matheatics & Applications, vol.3, no.1, 13, 199-1 ISSN: 179-9687 (print), 179-979 (online) Scienpress Ltd, 13 Solving initial value probles by residual power series ethod Mohaed H. Al-Sadi
More informationIntelligent Systems: Reasoning and Recognition. Perceptrons and Support Vector Machines
Intelligent Systes: Reasoning and Recognition Jaes L. Crowley osig 1 Winter Seester 2018 Lesson 6 27 February 2018 Outline Perceptrons and Support Vector achines Notation...2 Linear odels...3 Lines, Planes
More informationPattern Recognition and Machine Learning. Learning and Evaluation for Pattern Recognition
Pattern Recognition and Machine Learning Jaes L. Crowley ENSIMAG 3 - MMIS Fall Seester 2017 Lesson 1 4 October 2017 Outline Learning and Evaluation for Pattern Recognition Notation...2 1. The Pattern Recognition
More informationEMPIRICAL COMPLEXITY ANALYSIS OF A MILP-APPROACH FOR OPTIMIZATION OF HYBRID SYSTEMS
EMPIRICAL COMPLEXITY ANALYSIS OF A MILP-APPROACH FOR OPTIMIZATION OF HYBRID SYSTEMS Jochen Till, Sebastian Engell, Sebastian Panek, and Olaf Stursberg Process Control Lab (CT-AST), University of Dortund,
More informationQ5 We know that a mass at the end of a spring when displaced will perform simple m harmonic oscillations with a period given by T = 2!
Chapter 4.1 Q1 n oscillation is any otion in which the displaceent of a particle fro a fixed point keeps changing direction and there is a periodicity in the otion i.e. the otion repeats in soe way. In
More information