Open access peer-reviewed chapter

Quantum Perturbation Theory in Fluid Mixtures

By S. M. Motevalli and M. Azimi

Submitted: June 12th 2012Reviewed: October 25th 2012Published: April 3rd 2013

DOI: 10.5772/54056

Downloaded: 1642

1. Introduction

Experimental assessment of macroscopic thermo-dynamical parameters under extreme conditions is almost impossible and very expensive. Therefore, theoretical EOS for further experiments or evaluation is inevitable. In spite of other efficient methods of calculation such as integral equations and computer simulations, we have used perturbation theory because of its extensive qualities. Moreover, other methods are more time consuming than perturbation theories. When one wants to deal with realistic intermolecular interactions, the problem of deriving the thermodynamic and structural properties of the system becomes rather formidable. Thus, perturbation theories of liquid have been devised since the mid-20th century. Thermodynamic perturbation theory offers a molecular, as opposed to continuum approach to the prediction of fluid thermodynamic properties. Although, perturbation predictions are not expected to rival those of advanced integral-equations or large scale computer simulations methods, they are far more numerically efficient than the latter approaches and often produced comparably accurate results.

Dealing with light species such as Heand H2at low temperature and high densities makes it necessary taking into account quantum mechanical effects. Quantum rules and shapes related with the electronic orbital change completely the macroscopic properties.

Furthermore, for this fluid mixture, the quantum effect has been exerted in terms of first order quantum mechanical correction term in the Wigner-Kirkwood expansion. This term by generalizing the Wigner-Kirkwood correction for one component fluid to binary mixture produce acceptable results in comparison with simulation and other experimental data. Since utilizing Wigner-Kirkwood expansion in temperatures below 50 K bears diverges, we preferred to restrict our investigations in ranges above those temperatures from 50 to 4000 degrees. In these regions our calculations provide more acceptable results in comparison with other studies.

This term make a negligible contribution under high temperatures conditions. Taking into account various contributions, we have utilized an improved version of the equation of state to study the Helmholtz free energy F, to investigate the effects of P and T on thermodynamic properties of helium and hydrogen isotopes mixtures over a wide range of densities. We also have studied effects of concentrations of each component on macroscopic parameters. In addition, comparisons among various perturbation and ideal parts have been presented in logarithmic diagrams for different densities and concentrations for evaluation of perturbation terms validity in respect to variables ranges.

The first section is dedicated to a brief description of Wigner expansion which leads to derivation of first quantum correction term in free energy. With the intention of describing effects of quantum correction term we have explained theoretical method of our calculations in the frame work of statistical perturbation theory of free energy in section two. In section three we have depicted diagrams resulted from our theoretical evaluations and gave a brief explanation for them. In section four we have focused on the description of our calculations and its usages in different areas. Finally, some applications of this study have been introduced in the last section.

2. Quantum correction term

Considering quantum system of Nidentical particles of mass mconfined to the region of Λwith the interacting potential ofU. This structure is considered in υ-dimension space (Rυ). In the absence of external fields the Hamiltonian of particles is given as

H=12m(i)2+U(r)E1

Where, is the Plank constant. The equilibrium statistical mechanics of the particle system is studied in the canonical ensemble at the temperature T(or, alternatively, the inverse temperature β=1/kBTwith kBbeing Boltzmann’s constant). Quantum effects will be considered via de Broglie wavelength λ=β/m. For a typical microscopic length of particles l, for sufficiently small dimensionless parameter λ/lsemi-classical regime is dominant. In such system Boltzmann density in configuration space rcan be expanded in powers of λ2within the well-known Wigner-Kirkwood expansion [1, 2]. In the case of an inverse-power-law repulsive potential V(r)=V0(a/r)nfrom the range 1<n<, the Wigner-Kirkwood expansion turns out to be analytic in λ2[3]. In the hard-core limit n, this expansion is not further correct and one has the non-analyticity of type (λ2)1/2, as was shown in numerous analytic studies [4-7]. In contrast to the bulk case, the resulting Boltzmann density involves also position dependent terms which are non-analytic in λ. Under some condition about the classical density profile, the analyticity in λis restored by integrating the Boltzmann density over configuration space.

2.1. Wigner-Kirkwood expansion

To have an analytical equation for quantum effects in fluid we must derive partition function of it. In approximating partition function we need to evaluate Boltzmann density. Consequently having an expansion of quantum correction terms it is necessary to expand Boltzmann density. Considering system of Nparticles in the infinite space in standard Wigner-Kirkwood expansion [1, 2] fermions or boson exchange effects between quantum particles have been neglected. In the “bulk” regime, equilibrium quantities of this system in the nearly classical regime can be expanded in powers of h2. In this section, we review briefly the derivation of this expansion for utilizing it in statistical perturbation framework. The Boltzmann density Bβin configuration space rcan be formally written in the basis of plane waves as a υN-dimensional integral defined in an infinite domain Rυ:

Bβ=r|eβH|r=dP(2π)υNe(i/)p.reβHe(i/)p.rE2

Where p=(p1,p2,p3,...)is the υN-dimensional momentum vector. Instead of considering we take the Laplace transform of this operator with respect to the inverse temperature β,

0βdβeβHeβz=1H+ZE3

Via integrating equation 2 in respect to βwe have

0βdβr|eβH|reβz=dP(2π)υNe(i/)p.r1H+Ze(i/)p.rE4

Let us introduce following definition

H+z=D+QE5

That Qand Drespectively represent

Q=12m(i)212mP2D=12mP2+U(r)+zE6

One can expand

1H+Z=1D1DQ1D+1DQ1DQ1D...E7

Q, operates in the following manner

Q[f(r)e(i/)p.r]=e(i/)p.r[imp.+22m2]f(r)E8

And then we can find that

e(i/)p.r1H+ze(i/)p.r=1Dn=0{[imp.+22m2]1D}ndP(2π)υNe(i/)p.r1H+ze(i/)p.r=dP(2π)υN1Dn=0{[imp.+22m2]1D}n=0βdβr|eβH|reβz=n=0dP(2π)υN1D{[imp.+22m2]1D}nE9

So we have expanded series in 2nwhich enable us power series of n. It remains to define 1/Dj

1Dj=0dβ1(j1)!βj1eβD=0dβeβz1(j1)!βj1eβ[p2/2m+U(r)]E10

and finally integrating on the momentum variables p, the Boltzmann density in configuration space is obtained as the series

r|eβH|r=n=0Bβ(n)(r),E11

where

Bβ(0)(r)=1(2πλ)υNeβUBβ(1)(r)=1(2πλ)υNeβUλ2[β42U+β26(U)2],Bβ(2)(r)=1(2πλ)υNeβU{λ2[β62Uβ28(U)2]+Oλ4}We conclude that the quantum Boltzmann density in configuration space is given, to orderλ2, byr|eβH|r=1(2πλ)υNeβU{1+λ2[β122U+β224(U)2]+Oλ4}=1(2πλ)υN{eβU[1λ2β242U]+λ2242eβU+Oλ4}E12

Integrating Boltzmann density ignoring exchange effects over configuration space will result in partition function of fluids mixture.

Zqu=1N!Vdrr|eβH|rE13

Substituting the λ-expansion of the Boltzmann density (12A) into formula (13), the quantum partition function takes the expansion form

Zqu=1N!Λdr1(2πλ)υN{eβU[1λ2β242U]1+λ2242eβU+Oλ4}E14

For expressing macroscopic physical quantities, one defines the quantum average of a function f(r)as follows

fqu=1ZquN!Λdrr|eβH|rf(r)E15

At the one-particle level, one introduces the particle density

nqu(r)=j=1Nδ(rrj)quE16

At the two-particle level, the two-body density is given by

nqu(2)(r,r)=j,k=1jkNδ(rrj)δ(rrj)quE17

And the pair distribution function

gqu(r,r)=nqu(2)(r,r)nqu(r)nqu(r)E18

The classical partition function and the classical average of a function f(r)are defined as follows

Z=1N!Λdr(2πλ)υNeβU(r)E19
f=1ZN!Λdr(2πλ)υNeβU(r)f(r)E20

Consequently with the definition of equation 19 one can derive below equation for Zqu

Zqu=Z{1λ2β242U+Oλ4}E21
βFqu=ln(Zqu)E22
ln(Zqu)=ln(Z)+ln(1λ2β242U+Oλ4+...)E23

Since we have ln(1x)=xx2/2...we can expand the second term in the right side. By means of equation 18 in deriving 2Uwe can have explicit formula for the second term of which indicates the first term of Wigner-Kirkwood correction part that is consist of the second derivative of potential function that leads to below equation for quantum correction term with the number density of nwe have

Fqu(1)=h2NAnβ96π2mσ02U(r)g(r)drE24

g(r)represents radial distribution function, which is a measure of the spatial structure of the particles in reference system, is the expected number of particles at a distance r. NAis Avogadro constant and σ0is the distance in which potential function effectively tend to zero.

2.2. Free energy

Generalizing to multi-component system we have [8]

Fqu(1)=h2NAnβ96π2i,jcicjmijσij02uij(r)gij(r)4πr2Vij¯drE25

m11=m1,m22=m2,m12=c1m1+c2m2

miis the ith particle’s mass. Vij¯is the average molecular volume. Distribution function defines probability of finding particle at particular point r. In many literatures that have studied distribution function found it more versatile to use Laplace transform of this function G(s).

Gij(s)=0rgij(r)esrdrE26

In this chapter the two formula which use RDF, we will encounter below integral equation that need expansion.

I=σ0rϕ(r)rg(r)dr=σ031xϕ(x)xg(x)dx=σ03{0x2ϕ(x)g(x)dx01x2ϕ(x)g(x)dx}E27

On the right side of above equation from the right in the first equation we approximate distribution function with its values at contact points. This choice has been resulted from the behavior of molecules of which their repulsive interactions dominate their attractive potential. However, for the second term (I) we will use change in integrals to employ Laplace transform of RDF instead of RDF directly.

xg(x)=12πiγiγiG(s)esxdsE28

Substituting above equation in Iwe have

I=0φ(s)G(s)dsE29

Where φ(s)represents

φ(s)=12πiγiγixϕ(x)esxdxE30

That indicates inverse Laplace of xϕ(x). So it suffices to just define inverse Laplace of potential function multiplied by x.

Therefore, Using Laplace transform of RDF G(s)[9] quantum correction term for DY potential turn out to be

FQ=h2NAnβ24πi,jcicjεijAijVij¯mijσij0(λij2eλijG(λijσij0)υij2eυijG(υijσij0))E31

ciis the i particle’s concentration and nrepresents number density. Aij, λijand υijare controlling parameters of double Yukawa(DY). εijis the attractive well depth of mutual interacting potential.

3. Framework

The derivation of the thermodynamic and structural properties of a fluid system becomes a rather difficult problem when one wants to deal with realistic intermolecular interactions. For that reason, since the mid-20th century, simplifying attempts to (approximately) solve this problem have been devised, among which the perturbation theories of liquids have played a prominent role [10]. In this instance, the key idea is to express the actual potential in terms of a reference potential (that in terms of Ross perturbation theory Helmholtz free energy is expressed as of the “unperturbed” system) plus a correction term. This in turn implies that the thermodynamic and structural properties of the real system may be expressed in terms of those of the reference system which, of course, should be known. In the case of two component fluids, a natural choice for the reference system is the hard-sphere fluid, even for this simple system the thermodynamic and structural properties are known only approximately. Let us now consider a system defined by a pair interaction potential u(r). The usual perturbation expansion for the Helmholtz free energy, F, to first order in β=1/kβT, with Tbeing the absolute temperature and kBbeing the Boltzmann constant, leads to F. Common starting point of many thermodynamic perturbation theories is an expansion of the Helmholtz free energy, the resulting first-order prediction for a fluid composed of particles helium and hydrogen is given via the following equation.

F=Ft+FQ+FHB+FidE32

The terms respectively are perturbation, Quantum, hard convex body and ideal terms. Perturbation term due to long range attraction of potential is given by [10]

Ft=2πni,jcicjσij0uij(r)gijHS(r,ρ,σij)4πr2Vij¯drE33

Via Laplace transform of RDF (rgijHS) in calculation of first order perturbation contribution due to long-ranged attraction for DY potential we can employ below equation:

Ft=kTi,jcicjεijσij0AijV¯ij(eλijG(λijσij0)eυijG(υijσij0))δFtE34

Vij¯the average molecular volume defined as:

Vij¯=1+(ni1)σii3[32(σii2+σij2)li12(li)332[(σii+σjj)2li2]12sin1[liσii+σjj](σij2)]E35

Where nidefine the number of element in a molecule, liis distance of centre to centre for each molecule. δFtcorresponds to the interval of [σij,σij0]which long range attractive range is not further applicable. Consequently, we prefer to use the contact value of hard sphere RDF at r=σij. By this approach we can express this term as:

δFtn2kTi,jcicjσijσij0uij(r)gijHS(σij)4πr2Vij¯drE36

gijHS(σij)is the contact value of radial distribution function. σijstands for separation distance at contact between the centers of two interacting fluid particles, with species iand j. Although via minimization of Helmholtz free energy we can achieve value for hard sphere diameter, we preferred to use its analytical form due to its practical approach [17]. Hard sphere diameter will be calculated by means of Barker-Henderson equation as a function of interacting potential and temperature. Using Gauss-Legendre qudrature integration method we are able to evaluate its values numerically.

σij=0σij0(1exp(βuijDY(r)))drE37

FHB,Helmholtz free energy of hard convex body is given by following equation:

FHB=amix(FHS+Fnd)E38

Non-sphericity parameter amixfor the scaling theory [11] is defined as

amix=13πi,jcicjVijef(Vijeff)(Vijeff)i,jcicjVijeff,Vijeff=π6σij3Vij¯E39

(Vijeff)and (Vijeff)are the first and second partial derivatives of Vijeffwith respect to σiiand σjj. From Boublik, Mansoori, Carnahan, Starling, Leland (BMCSL) [12, 8] with correction term of Barrio [13] on EOS, the Helmholtz free energy, FHSfor hard sphere term becomes:

FHSKT=η3[ξ1+(2η3)ξ2]1η3+η3ξ3(1η3)2+(ξ3+2ξ21)ln(1η3),

ξ1=3η1η2η0η3,ξ2=η1η2η32(η4z1+η0z2),ξ3=η23η0η33

z1=2c1c2σ11σ22(σ11σ22σ11+σ22)z2=c1c2σ11σ223(σ112σ222)E40

The correction term due to nonadditivity of the hard sphere diameter is the first order perturbation correction [14]

Fnd=kTπnc1c2(σ11+σ22)(σ11+σ222σ12)g12HS(σ12)E41

In Eq. (41), g12HS(σ12)refer to as hard sphere radial distribution function at r=σ12contact point by inclusion of Barrio and Solana correction on equation of state of BMCSL. Undoubtedly, the availability of the analytical HS RDF obtained from the solution to the corresponding Percus–Yevick (PY) equation represented a major step toward the successful application of the perturbation theory of liquids to more realistic inter-particle potentials. However, the lack of thermodynamic consistency between the virial and compressibility routes to the equation of state present in the PY approximation (as well as in other integral equation theories) is a drawback that may question the results derived from its use within a perturbation treatment. Fortunately, for our purposes, another analytical approximation for the RDF of the HS fluid, which avoids the thermodynamic consistency problem, has been more recently derived [15, 16]. We used improved RDF that yields exact asymptotic expression for the thermodynamic properties. However, we have used improved version of RDF that yields exact asymptotic expression for the thermodynamic properties. This have been derived by inclusion of Barrio and Solana correction on EOS of BMCSL at r=σ12[9]

gijHS(σij)=gijBMCSL(σij)+gijBS(σij)g12BMCSL(σ12)=11η3+3η2(1η3)2σiiσjjσii+σjj+2η22(1η3)3(σiiσjjσii+σjj)2,g12BS(σ12)=1δijci2η1η2(1η3)2(σiiσjjσii+σjj)2(σ11σ22)(δij+(1δij)σ22σ11)E42

ηi=π6njcjσjji

δijis the Kronecker delta function. For additive mixtures σijis arithmetic mean of hard-core diameters of each species. Otherwise, the system is said to be non-additive.

The ideal free energy with Nparticle for the atomic and molecular components of fluid mixture are given by,

Fid(n,T)=32ln(h22πkTm1c1m2c2)+lnn+icilnci1E43

Compressibility factor of ideal term is one and ZHBwould be estimated with the following derivation of related Helmholtz free energy

ZHB=nnFHBkTE44

For the perturbation term due to long rage attraction of potential tail employing (44) we will have

Zt=2πnkTi,jcicjεijσij0AijV¯ij(eλij(G(λijσij0)nnG(λijσij0))eυij(G(υijσij0)nnG(υijσij0)))δZt

δZtn2kTi,jcicj(gijHS(σij)+nngijHS(σij))σijσij0uijDY(r)4πr2Vij¯dr

Numerical integration has been used for calculation of δZtin the range of [σij,σij0]. Expressions for first order perturbation and quantum correction term of compressibility factor are achievable via applying (44) for the free energy part of the quantum correction term.

ZQ=h2NAnβ224πi,jcicjεijAijVij¯mijσij0(λij2eλij(G(λijσij0)nnG(λijσij0))υij2eυij(G(υijσij0)nnG(υijσij0)))E45

Summation over compressibility factors gives the total pressure of mixture

P=nkT(1+ZHB+Zt+ZQ)E46

Defining Gibbs free energy provides information at critical points of phase stability diagram. Concavity and convexity of Gibbs diagram indicates if mixture is in one phase or not,

G=F+NnPE47

Furthermore, Gibbs excess free energy is an appropriate measure in the definition of phase stability. Negative values for this energy describe stable state. This is expressed as

Gxs=GiciGi0NkTicilnciE48

That Gi0represents the Gibbs free energy of pure fluid of species i. Concentration-concentration structure factor is defined as

Scc(0)=NkT(2c2G)1T,P,NE49

Compairing this equation with Sccidenable us to define degree of hetero-coordination. In a given composition if Scc(0)<<Sccidthen unlike atoms tend to pair as nearest neighbors (hetero-coordination) and when Scc(0)>>Sccidthen like atoms are preferred as a neighbor.

3.1. Potentials

It is convenient to consider interacting potential with short-range sharply repulsive and longer-range attractive tail and treat them within a combined potential. The most practical method for the repulsive term of potential is the hard-sphere model with the benefit of preventing particles overlap. Furthermore, attractive or repulsive tails may be included using a perturbation theory. It is incontrovertible to generalize this potential to multi-component mixtures. This behavior is conveyed in double Yukawa (DY) potential which provides accurate thermodynamic properties of fluid in low temperatures and high density [18, 19]. At first we define DY potential as its effects on pressure of HeH2mixture has been studied in this work, written as:

uijDY(r)=εijAijσij0r[eλij(1r/σij0)eυij(1r/σij0)]E50

Aij, λij, υijare controlling parameters. These parameters for Heand H2are listed in table 1 with their reference [20].

He−HeHe-H2H2-H2
σ02.6342.9702.978
A2.5482.8013.179
ε/kB10.5715.5036.40
υ3.3363.3863.211
λ12.20410.9549.083

Table 1.

Potential parameters for He, H2interactions for DY potential [20].

For the atomic and molecular fluids studies in this mixture, these particles interact via a exponential six (exp-6) or Double Yukawa (DY) potential energy function [20]. The fluids considered in this work are binary mixtures that their constituents are spherical particles of two species, iand j, interacting via pair potential uij(r).

uijexp6(r)={r<σεi,jαijαij6(66αijexp(αij(1rσmin,ij))(σmin,ijr)6)r>σE51

So we consider two-component fluid interacting via Buckingham potential uij(r)between molecules of types iand j. This potential is more realistic than square-well or Yukawa type potential for hydrogen isotope’s mixture [21] at high temperatures. Because of same atomic structure of hydrogen and its isotopes, the three constant of potential are same for hydrogen isotopes. These constants obtained experimentally from molecular scattering [22].

σminindicate the range of interaction and the parameter αregulates the stiffness of repulsion. For hydrogen and helium type atoms these parameters have been organized in Table 2. It is well known that the long range attractive part of exp-6 potential is similar to Lenard–Jones potential.

In view of the energy equation (32), one can readily obtain equation for total pressure and different contributions to pressure from standard derivation of respective Helmholtz free energy. By the exp-6 potential, we have computed the Helmholtz free energy. The ten-point Gausses quadrature has been used to calculate integrals in quantum correction and perturbation contribution. The calculated pressure for D2+T2fluid mixture with equal mole fraction and at temperature of T=100oKis showed logarithmically in figure 1. As it is clear from this figure, the effect of hard sphere term of pressure in given rang of temperature is significant and the range of pressure variation is wider than ideal part. As it is mentioned earlier the difference between isotopes is simply related to the neutron number in each nucleus and correction due to difference in mass which involves in non-additive correction term doesn’t affect sensibly. In addition, this figure shows the predicted equimolar surface of the deuterium and tritium mixture for quantum correction term. This part is the most significant contribution at low temperature and varying smoothly in higher temperature. At very high densities, perturbation term contribution increases sharply with reducing density. Also, terms, Pt, PQ, PHSand Pid, tend to infinity as ρ.

He−HeHe-H2H2-H2
α13.1012.711.1
ε/kB10.8015.5036.40
σmin0.296730.3370.343

Table 2.

Potential parameters for He, H2 interactions for exp-6 potential [20]

Figure 1.

Different contributions of pressure as a function of reduced density for T=100oK for fluid mixture of deuterium and tritium

4. Results

For helium-hydrogen mixtures different parts of pressure due to correction terms and ideal parts have been showed in figure 2 at T=100oK. Ideal pressure at reduced density of approximately zero, to about 0.25 rises drastically. However, afterward it soars gently up to 100M (pa). Pressure due to hard sphere is the most significant contribution except that it is less than perturbation part at value of 1.5 for reduced density. Effects of perturbation and Quantum correction are important in high densities. In low densities, these contributions are insignificant and may possibly be ignored. Non-additive part has been caused by dissimilarity of particles which surges steadily from the beginning.

Figure 2.

Different contribution of correction terms on pressure of helium-hydrogen mixture at T=100, che=0.5 vs. reduce density

Gibbs excess free energy which is a measure for indicating phase stability of matters has been depicted in figure 3. Stability is limited to the areas that Gibbs excess free energy tends to negative values. This figure explains that stability rages for helium-hydrogen mixture at room temperature is confide in the boundaries in which helium concentration is less than 0.1.

Figure 3.

Gibbs excess free energy for helium-hydrogen mixture

Table 3 presents a comparison between results of pressure from this work using DY potential in place of exp-6, Monte–Carlo simulations and additionally study of reference [23] Obviously, there are appreciable adaption among our investigation results and MC which proves validity of our calculations. As Table 3, exhibits in low temperatures DY potential have more consistent results in comparison with exp-6. However, values of pressure extracted using DY potential cannot adjust with simulation resembling exp-6. Moreover, at higher temperatures after T=1000oK, DY potential is not good choice for evaluating EOS of hydrogen and helium mixture. We clarify our deduction presenting comparison between effects of these two potentials over wide ranges of temperatures.

T(K)cHeρ*ηPMCP[19]P[23]PDY[19]
3000.251.1010.4332.30902.70391.96642.8678
3000.51.1010.4001.85601.70011.57291.8402
3000.751.1010.3671.42401.28161.31601.3887
10000.51.2230.3354.51004.42054.10944.9406
10000.751.2230.3073.71503.51903.59043.9328
40000.51.3760.24712.430012.083212.101414.154
40000.51.5720.28216.330016.448516.472019.859

Table 3.

Comparison of values of pressure results from our study [19], Monte-Carlo simulation [24] and Isam Ali’s study [23].

Figure 4.

Comparision of efect of DY and EXP-6 potential on pressure of mixture in che=0.5, T=300 vs. Reduced density

Providing evidence of gradual divergence of DY and exp-6 potentials, a comparative figure has been made in figure 4 for helium-hydrogen mixtures. This figure shows more steepening effects of DY on total pressure of this mixture. Both potentials engender increase in pressure, except that, Buckingham affects moderately on pressure increase. The exp-6’s more steady behavior makes it adjustable with previous studies and MC simulation.

In figures 5, 6, 7, 8 we tried to give information about effects of quantum correction term on total pressure of helium-hydrogen and deuterium-tritium mixtures at the high reduced density of 1.3. This correction term has been plotted in 3-dimensional diagram in figure 5. This term is approximately zero for temperatures higher than 200 (K). Figure 6 represents that for hydrogen rich mixture at low temperature due to quantum effects pressure rise is significant. For effectual discussion on the effects of this term we have described PQ/Pin figure 7 for helium-hydrogen and in figure 8 for deuterium tritium mixture. For the third picture increase in pressure is similar to what have been elaborated for figures 5 and 6. For figure 8 this manner remains analogous to helium-hydrogen mixture and temperatures next to 100 (K). However, for temperature lower than this it would behave inversely. For this range any increase in tritium concentration bears decrease in pressure.

Figure 5.

Pressure of quantum correction term at ρs=1.3 for helium-hydrogen mixture.

Figure 6.

Total pressure from 50 K at ρS=1.3 for helium-hydrogen mixture.

Figure 7.

Fraction of quantum perturbation term to total pressure for helium-hydrogen mixture.

Figure 8.

Fraction of quantum perturbation term to total pressure for Deuterium-Tritium mixture.

5. Conclusion

An Equation of state of hydrogen–helium mixture has been studied up to 90G (pa) pressure and temperature equal to 4000◦K. We have used perturbation theory as an adequate theory for describing EOS of fluid mixtures. As well, by using this theory we can add extra distributive terms as perturb part which makes it more applicable than other theories. Considering this advantage, we can spread it out with additional terms for investigation on other states of matter like plasma in the direction of compares with experimental data. Otherwise, using simulation methods, for evaluating our theoretical results. Such as ab initio simulations with the code VASP,[25] which combines classical molecular dynamics simulation for the ions with electrons, behave in quantum mechanical system by means of finite temperature density functional theory [26]. In this chapter, two potentials have been presented, which we have used them for hydrogen isotopes and helium, and their mixtures. By means of comparison with Monte Carlo simulation and results of refrence [14] in Table 3 we could prove that exp-6 potential is more beneficial than DY in wider ranges of variables, since its application in this theory shows more convergent results in comparison with MC simulation [24]. Also exp-6 potential is a good choice of potential since it allows us to elevate temperature and density [28]. But as hydrogen molecules dissociation occurs [28] for pressures more than 100G (pa), this effect must be accounted. Therefore, we have restricted ourselves to pressures below 100G (pa).

Furthermore, we have used Wertheim RDF which enables us to use this EOS for extended values of temperature. As well, we have compared different contributions of pressure to represent which one is more effective in different density and temperature regimes. By finding the most effective parts of pressure contributions in each ranges of independent variables (Temperature, reduced density, mole fraction), we can omit the less significant parts which are considered ignorable in value, to decrease unnecessary efforts. Likewise, we can speculate from Fig. 1 that in low temperature and high densities, long range perturbation term has the most significant effect in comparison with other parts. On the other hand, hard sphere part can be assumed as the most noticeable part in high temperature ranges. Moreover, comparison of DY and exp-6 potentials effects, on pressure of this mixture has been studied to express benefits of using exp-6 potential for higher temperatures and densities. Additionally, as it is obvious in high temperature and density difference between effects of two potentials are considerable for this equimolar mixture. This discriminating property makes exp-6 potential preferable.

Furthermore, this approach has been used to evaluate EOS of D2+T2mixture. Also, we have used this method to compare different contribution parts of pressure. These comparisons indicate that in low temperature quantum effects are more important, however in high temperatures, hard sphere part is the most effective. The last two three dimensional diagrams reveals the importance of quantum term in comparison with total pressure. However, for temperatures below 100 (k) for deuterium-tritium mixture negative pressure express that in low tritium concentration, deuterium rich fluid tend to consolidate.

6. Applications

One of the topics which can count on a great deal of interest from both theoretical and experimental physics is research in fluid mixture properties. These interests, not only comprise in the wide abundance of mixtures in our everyday life and in our universe but also the surprising new phenomena which were detected in the laboratories responsible for this increased attention. Mixtures, in general, have a much richer phase diagram than their pure constituents and various effects can be observed only in multi-component systems.

These kinds of studies have allowed a more complete modeling of mixture and consequently a better prediction and a more accurate calculation of thermodynamic quantities of mixture, such as activity coefficient, partial molar volume, phase behavior, local composition in general and have promoted a deeper understanding of the microscopic structure of mixtures.

Furthermore, for astronomical applications it is known that most of giant gas planets are like Jupiter is consisting primarily of hydrogen and helium. Modeling the interior of such planets requires an accurate equation of state for hydrogen-helium mixtures at high pressure and temperature conditions similar to those in planetary interiors [29]. Thus, the characterization of such system by statistical perturbation calculations will help us to answer questions concerning the inner structure of planets, their origin and evolution [29, 30].

In addition, in perturbation consideration of plasma via chemical picture, perturbation corrections will be included by means of additional free energy correction terms. Therefore, in considering transition behavior of molecular fluid to fully ionized plasma these terms are suitable in studying the neutral interaction parts. Consequently this will help us in studying inertial confinement fusion [31] and considering plasma as a fluid mixture in tokomak [32].

© 2013 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

S. M. Motevalli and M. Azimi (April 3rd 2013). Quantum Perturbation Theory in Fluid Mixtures, Advances in Quantum Mechanics, Paul Bracken, IntechOpen, DOI: 10.5772/54056. Available from:

chapter statistics

1642total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Quantal Cumulant Mechanics as Extended Ehrenfest Theorem

By Yasuteru Shigeta

Related Book

First chapter

Mutiple Hopf Bifurcation on Center Manifold

By Qinlong Wang, Bo Sang and Wentao Huang

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us