Parameters of State in the Global Thermodynamics of Binary Ideal Gas Mixtures in a Stationary Heat Flow

In this paper, we formulate the first law of global thermodynamics for stationary states of the binary ideal gas mixture subjected to heat flow. We map the non-uniform system onto the uniform one and show that the internal energy U(S*,V,N1,N2,f1*,f2*) is the function of the following parameters of state: a non-equilibrium entropy S*, volume V, number of particles of the first component, N1, number of particles of the second component N2 and the renormalized degrees of freedom. The parameters f1*,f2*, N1,N2 satisfy the relation (N1/(N1+N2))f1*/f1+(N2/(N1+N2))f2*/f2=1 (f1 and f2 are the degrees of freedom for each component respectively). Thus, only 5 parameters of state describe the non-equilibrium state of the binary mixture in the heat flow. We calculate the non-equilibrium entropy S* and new thermodynamic parameters of state f1*,f2* explicitly. The latter are responsible for heat generation due to the concentration gradients. The theory reduces to equilibrium thermodynamics, when the heat flux goes to zero. As in equilibrium thermodynamics, the steady-state fundamental equation also leads to the thermodynamic Maxwell relations for measurable steady-state properties.


I. INTRODUCTION
In the classical thermodynamics the internal energy, U (S, V, N ), of a one-component ideal gas is a function of three parameters: entropy, S, volume V and the number of particles N .Each parameter of state represents one independent way of system's energy exchange with the external world.For fixed N there are two ways of energy change: heat and work.The infinitesimal change of the internal energy satisfies the equation: dU = T dS − pdV , where the first term is the heat and the last is the work term.Untill recently no such description was available for systems in non-equilibrium states, subjected to energy flow.
The classical theory of irreversible (non-equilibrium) thermodynamics [1] is based on three differential nonlinear equations representing conservation of mass, momentum (Navier-Stokes equation) and energy.These conservation laws are supplemented by the assumption of the local equilibrium, corresponding local equations of state, and constitutive relations between fluxes and thermodynamic forces.The solutions of these equations are given in terms of velocity, v(r, t), temperature, T (r, t) and number density of particles, n(r, t) profiles.In the stationary states the profiles and fluxes depend only on the position in space r, but not on time, t.The assumption of the local equilibrium is crucial in the formulation of the local equations of state.For the ideal gas it is justified, but contested for interacting systems [2].The assumption of local equilibrium for the ideal gas is valid for small temperature gradient i.e. such that, l f p |∇T | /T ≪ 1, for the mean free path of the molecules, l f p [3].At the pressure of 1 bar at room temperature, the mean free path is of the order of l f p ≈ 100nm.In such conditions the assumption of local equilibrium breaks down only for temperature gradients higher than 10 7 K/cm.Thus this assumption is even satisfied inside the Sun!This local description contains the first law of thermodynamics, but in a different form than that given in the equilibrium thermodynamics and based on a few global parameters of state.In our recent paper [4] we provided the latter, i.e., global thermodynamic description for the ideal gas in a heat flow.
We have previously [4] studied the one-component ideal gas in a heat flow between two parallel walls at distance L, kept at two different temperatures T 1 (z = 0) > T 2 (z = L).The local equilibrium gives the local pressure, p(z), and the internal energy per unit volume u(z): and with Boltzmann constant k B , particle number density n (z), and the temperature T (z) at position z.We have shown rigorously that when we integrate both equations over the volume of the system we can formulate the global thermodynamics with the internal energy as a function of a few parameters of state.After integration we get: and After definition of system's temperature, T * , In this way we made a mapping of the non-uniform system into the uniform one.It is worth mentioning that the local equilibrium assumption is not needed.If we had different local equations of state we could still integrate them over the volume and perform aforementioned mapping.We observe that obtained equations have the same form as in equilibrium for the ideal gas at temperature T * .We demanded such form because after the mapping we treat the system as the uniform one as in equilibrium.Moreover the obtained equations of state must reduce the equilibrium equations of state when the heat flux is zero.Now we define the internal energy as a function of three parameters of state U (S * , V, N ), with the thermodynamic relation: This mapping gives us the same formal structure as we know from equilibrium.The entropy S * is responsible for the net heat that enters or leaves the system [5][6][7] and changes the internal energy.In general the heat flows through the system all the time without changing the internal energy.Upon any process we would like to know how much of the heat transferred to the system changes the internal energy.This heat is called net heat and is given in the differential form by T * dS * .S * is only part of the total entropy of the system.The total entropy is given by: Only S * governs net heat in the system (heat absorbed/released in the system).∆S, controls the dissipative background and solely depends on the entropy production given by [1] where κ is the heat conductivity.The difference between the total entropy and S * vanishes, ∆S (U, A, L, T 2 /T 1 ) → 0, when the system approaches the equilibrium state, T 2 /T 1 → 1.Therefore, S * becomes in this limit the equilibrium entropy.Nonetheless the formal dependence of S * on U and V at non-equilibrium state in a heat flow is the same as at equilibrium.
In this paper we want to apply the same mapping procedure to the binary mixture of ideal gases.The big difference between previous work and the current one is the fact that apart from the number density profile and the temperature profile we have additional profiles of the number densities of each component.As we shall see, these additional profiles lead to new parameters of state in the non-equilibrium state, which have no direct counterpart at equilibrium.The purpose of this work is to formulate the first law of global thermodynamics for ideal gas binary mixture in the heat flow.The paper is organized as following: in section II we recall the equilibrium properties of the ideal gas mixture.In section III we discuss this mixture enclosed between two parallel walls at different temperatures and solve the equations of irreversible thermodynamics.In section IV we define all parameters of state for the mixture and perform the mapping of non-uniform system into the uniform one.We introduced the first law of non-equilibrium thermodynamics, which follows from these parameters, and discuss some of its consequences in section V. We determine the discuss the results in section VI.

II. PRELIMINARIES
We consider a binary mixture of ideal gases enclosed between two parallel walls separated from each other by a distance L in the z direction.The volume of a system is V = AL, where A is an area in the x− y plane.The components of the mixture have number densities n 1 and n 2 , such that n = N/V = n 1 + n 2 is the total density.The equation of state of ideal gas at pressure p and temperature T p = nk B T , where k B is the Boltzmann constant, can be written as a sum of partial pressures which is the Dalton's law [8].
The internal energy density (per volume) is the sum of the internal energy density u i (i = 1, 2) of the two components separately considered: From the classical equipartition theorem applied to ideal gas it follows that where f i , i = 1, 2 are translational and rotational degrees of freedom of the i component.f i takes the value 3 for monoatomic, 5 for diatomic, and 6 for polyatomic gas.
Often, the ideal gas law ( 9) is written using dimensionless specific heat capacity at constant volume c The total energy and entropy are given by (which reduces to U = V u and S = N s for homogeneous systems) where the equilibrium entropy density ns = ns (0) + ns mix of a binary mixture is Here ) is the entropy per particle of the two components separately considered.s mix is the mean-field expression for the mixing entropy.For ideal gas one has where Φ i is independent of the thermodynamic state of the gas and has dimension of n 2/fi /T .For a monoatomic ideal gas, a quantum mechanical theory of Sacur-Tetrode predicts that the constant Φ i depends only upon the mass of the gas particle [9,10].
From fundamental thermodynamic relation in terms of entropy, the differential dS has the following form Using expressions for the entropy (Eqs.( 11) and ( 12)), the chemical potential of each component of a mixture is obtained as µi where is the chemical potential of the component i individually considered; it depends only on temperature T and pressure p: We define the chemical potential difference as In the absence of external fields, the equilibrium ideal gas is homogeneous.

III. IDEAL GAS MIXTURE IN HEAT FLOW
Now, we introduce heat flow into the system by setting different temperatures on the walls, i.e., In the non-equilibrium state induced by this boundary condition, the system becomes inhomogeneous and one has to consider spatially varying temperature T (r), pressure p (r), density n (r) and difference in number densities ϕ (r).In the hydrodynamic limit, the time evolution of the binary mixture is given by the conservation laws for n, ϕ, momentum, and energy supplemented by the assumption of local equilibrium, relations between thermodynamic forces and fluxes, and thermodynamic equations of state [1].Note that ideal gas satisfies the local equilibrium exactly.This fact describes the strict absence of spatial correlations between particles, which is, of course, not the case for non-ideal systems.
We focus here on a stationary state with vanishing gas velocity field and constant pressure across the system.The latter, together with n (r) and ϕ (r), follows from conservation of mass and momentum.The assumptions v = 0 and p (r) = p are in agreement with our previous simulation studies [11][12][13] (for gas-liquid evaporating systems [11,12] and for Lennard-Jones fluid volumetrically heated [13]), which show that mechanical equilibrium is established very fast (in comparison to heat flow).Thus, the balance equations for ϕ and energy simplify to where is the diffusion current and J q is the heat current.We choose the following phenomeno- Schematic illustration of a binary ideal gas mixture between two parallel walls held at different temperatures T2 and T1.The resulting temperature gradient induces a gradient of fractions xi of the mixture components.Solid (green) lines indicate profiles of x1(z/L).Two cases are shown: (i) f1 > f2 for which the profile is given by Eq. ( 26) and (2) f1 = f2 for which the profile is given by Eq. (28).For the latter case, we took α12 = 0.3432, which corresponds to the Ne-He pair at T1 = 300 K.In both cases the reduced temperature gradient s = (T2 − T1)/T1 was chosen to be 0.5 and N1 = N2.
logical expressions for J d and J q [1]: In steady state, J d = 0 because the system has no contact with particles reservoir.First, we ignore the thermodiffusion effect, i.e., we assume that L 12 = 0. We also neglect the other cross-term by setting L 21 = 0.As a result, Eq. ( 19) reduces to and The latter yields a linear temperature profile (the system is translationally invariant in x, y directions) From Eq. ( 20) we determine the relationship between the local fraction x 1 (z) of the component 1 and the temperature.The local equilibrium allows us to use the local form of Eq. ( 16) for µ/T ; it depends on T (z), x 1 (z), x 2 (z), and the pressure p, which is constant throughout the system.Thus where h 1 is the enthalpy per particle and ∇ T is the gradient at constant temperature; it involves the derivative of the chemical potential with respect to composition, i.e., Assuming that the reference state has temperature T 1 and using the fact that h i = H i /N i = (f i /2 + 1)k B T for the ideal gas we find where Note that if f 1 = f 2 , temperature inhomogeneity does not induce spatial variation of the composition, unless we take into account the effect of thermodiffusion.In that case we have generally Many numerical results as well as measurements of thermodiffusion coefficients are reported in the literature (see e.g.[14,16,17] and references therein).They indicate that these coefficients are functions of composition and temperature.Within the framework of the generalized Stefan-Maxwell thermodiffusion equations for f 1 = f 2 , one obtains exactly the same equation as above, but with 2L 12 /(L 11 k B T ) replaced by α 12 = A 12 /D 12 .A ij is the Newman-Soret thermal diffusion coefficient and D ij is the Stefan-Maxwell diffusivity (see e.g., Eq. (6.8) from Ref. [17]).It is argued that coefficients A ij = D T i /ρ i − D T j /ρ j , where D T i is the thermal diffusivity coefficient of component i and ρ i is its mass density, are more approximately constant with respect to composition in binary systems [20].Assuming that the ratio α 12 does not depend on both temperature and composition, the solution of Eq. ( 27) for the profile x 1 (z) is This solution clearly shows that the difference f 1 − f 2 is a relevant parameter.An analogy can be drawn here with the order parameter: for non-zero f 1 − f 2 , the profile x 1 (z) is qualitatively different from that for the zero difference f 1 − f 2 .Specifically, for f 1 − f 2 > 0, the concentration of the first component of the mixture (with more degrees of freedom) increases towards the hotter wall, while for f 1 = f 2 the opposite is true (see Fig. 1).This also holds if, due to the concentration and temperature dependence of the thermodiffusion coefficients, the functional form of x 1 -dependence on T is different than that given by Eq. (28).For noble gas pairs at T = 300 K, the value of α 12 is equal to 0.3432 for Ne-He pair, 0.1741 for Ar-Ne pair, or 0.0262 for Xe-Kr pair [21,22].If we assume that Onsager coefficients L 12 and L 11 are independent on the concentration and temperature [18], the solution of Eq. ( 27) becomes The density profile follows from the local equation of state and the condition of constant pressure at steady state, n (z) = p/k B T (z).For a given number of particles, Finally, the total energy of a mixture is the sum of the energies of the both components of the mixture (see Eqs. (10), (8), and ( 9)) Thus, the set of independent parameters controlling the stationary state is T 1 , T 2 , A, L, N, N 1 , f 1 , f 2 .

IV. NON-EQUILIBRIUM PARAMETERS AND FUNCTIONS OF STATE
One can rewrite the total energy (per volume) as follows where This has the same form as in equilibrium but with temperature T replaced by T * and f 1 , f 2 replaced by parameters f * 1 , f * 2 , which we call the "effective" degrees of freedom.From Eq. ( 31) it follows that 34) where we have used n(z)T (z) = p/k B = const.The parameter T * can be determined from the requirement that the pressure (35) Then, Eq. ( 35) can be written as which is to say that partial pressures in steady state are related to the non-equilibrium energies in the same way as in equilibrium (see Eq. ( 9)) but with f i replaced by f * i .Solving Eq. ( 36) with the use of Eq. ( 30) we find T * can be interpreted as an "average" temperature in the non-equilibrium steady state.Eliminating p in expressions for f * 1 , f * 2 we find: Using x 2 (z) = 1 − x 1 (z) we can see that the new parameters of state f * 1 , f * 2 are not independent but obey the following relation: We note that even if both components of the mixture have the same number of degrees of freedom f 1 = f 2 , the effective parameters f * 1 , are not equal in non-equilibrium steady states.Using Eqs. ( 26) and ( 22) for the profiles of the fraction x 1 (z) of the component 1 and temperature T (z), we can express the steady state variable f * 1 in terms of the control parameters.The explicit formulas can be obtained for the case of f 1 = f 2 in terms of the special functions as presented in Appendix A. For small, reduced temperature gradients s = (T 2 − T 1 )/T 1 and for N 1 = N 2 we find If N 1 = N 2 , the expression for f * 1 in this limit is more complicated, but has a similar structure in terms of dependence on s.The coefficients are functions of N 1 /N 2 and f 1 , f 2 but not of T 1 as shown in Appendix A. In the case of f 1 = f 2 , the integral over the fraction profile x 1 (z) given by Eq. ( 28) has the same form as for f 1 = f 2 (given in Appendix A), but with f 1 −f 2 replaced by −α 12 .For a profile given by Eq. ( 29) the integral cannot be expressed analytically.
In the next step in the construction of global thermodynamics, we note that because the equations (32) and (36) have the same structure as the equilibrium equation of state, we may formally write: with Equation ( 41) has the functional form of the equilibrium fundamental relation for a binary ideal gas mixture with f i replaced by f * i (compare Eqs ( 11) and ( 12)).We will treat Eqs (41)-(43) as a definition of the nonequilibrium steady state entropy S * (U, V, N 1 , N 2 , f * 1 , f * 2 ) that also provides the fundamental relation for the nonequilibrium steady state.Note that in comparison with the equilibrium entropy, S * depends on additional two state parameters f * 1 and f * 2 .As a consequence, for the difference dS * we find: where and (46) Equations ( 40) and (37) provide the effective parameters of state, T * and f * 1 for given control parameters T 1 , T 2 and N 1 , N 2 at fixed N 1 + N 2 .

V. FIRST LAW AND ITS CONSEQUENCES
We now consider the change of the total internal energy dU in our system.From Eq. (44) we have As argued in Ref. [4], the change in internal energy dU in the very slow process of moving from one steady state to another by varying the control parameters (so that the pressure remains homogeneous) is given by where dU = dQ is the net heat entering the system during this transition, and for the fixed number of mixture components the mechanical work done is The equation ( 48) is therefore the first law of thermodynamics for non-equilibrium steady states where the net heat identified by Eqs (48) and (44) (at constant N 1 and N 2 ) is We note that the net heat flow during the transition between two steady states is a combination of two exact differentials: the effective entropy dS * and the effective degrees of freedom df * 1 (from Eq. (39) it follows that . This is contrary to the equilibrium thermodynamics, in which heat depends only on temperature and the change in entropy.
At equilibrium, the most experimentally available thermodynamic quantities are response functions such as heat capacities, compressibility, susceptibility or chemical response functions.They follow from the first law of thermodynamics and the equilibrium fundamental relationship.Using the symmetry of the second derivatives, i.e., Maxwell relations, one can express one response function in terms of others.This is very useful as, for example, the heat capacity can be determined by measurements of other quantities, such as isothermal compressibility [8].
Once we have established the fundamental relationship for steady states of binary mixtures of ideal gases, we can generalize the equilibrium response functions to steady-state response functions.First, we consider thermal response function, i.e., the heat capacity C, which is a measure of amount of heat needed to raise the temperature of a system by a given amount.Generally, it is defined as derivative, C = dQ/dT , but depending on which independent variables are fixed during the measurements, one has different heat capacities.For example, the heat capacity at constant volume and number of components {N j } is defined as C V = dQ/dT | V,{N j } , whereas C p = dQ/dT | p,{N j } is the heat capacity at constant pressure and {N j }.While in the equilibrium state at a constant volume and fixed {N j } and {f j } we can change only the temperature of the system, in the nonequilibrium steady state we have more possibilities: we can independently vary T 1 and T 2 , or equivalently T 2 and the reduced difference s = (T 2 −T 1 )/T 1 .If one changes T 2 at fixed s then the state parameter f * 1 remains constant (see Appendix A and Eq. ( 40)) we can define steadystate heat capacities for constant volume and pressure as follows and Concerning mechanical response functions, we can generalize isothermal compressibility and thermal expansivity to the steady-states as follows Because the fundamental equation of steady-state has the same form as in equilibrium, the Maxwell relations for the path with s = const and f * 1 = const provide connection between the thermal and mechanical response function, which is the same as in the equilibrium, i.e., In order to determine C * p one need to measure the excess heat due to the small change of T 2 , which is in principle possible due to the recent development of the appropriate experimental techniques [23].The coefficient κ * T,Nj should be measured by changing the pressure at fixed both temperatures T 1 and T 2 , while the coefficient α * p,Nj should be determined by varying T 2 at fixed reduced gradient s.

VI. SUMMARY
The internal energy of the binary ideal gas mixture in a heat flow is the function of 5 parameters of state U (S * , V, N 1 , N 2 , f * ), irrespective of the number of boundary conditions.Each parameter of state defines one independent way of changing the internal energy of the system.The parameter f * is responsible for additional net heat, not included in T * dS * .T * dS * is not the differential form describing total net heat in the system.The total net heat is given not only by the changes in S * but also in f * .The same observation was made for the van der Waals gas in the heat flow.In this case the total net heat was given not only by S * , but also by the renormalized in the mapping procedure of two parameters of state a * , b * describing attractive interactions and excluded volume in the van der Waals gas [24].These parameters are constant at equilibrium since they are material parameters that define interactions in a particular system, but in non-equilibrium van der Waals gas they are parameters of state, which change the energy due to the change of density profiles.In our case of binary mixture the new parameter of state emerged from concentration profiles.Due to existence of profiles of different physical quantities, we expect all material parameters to become state parameters in the non-equilibrium systems. where For N 1 = N 2 we obtain where .
For small reduced temperature gradients s = (T 2 −T 1 )/T 1 across the system we find for where B(s; p, q) is incomplete Euler beta function [25].
For N 1 = N 2 the above expression reduces to Eq. (40).