Correlation Effects in Kinetics of One-Dimensional Atomic Systems

The paper is devoted to the analysis of the correlation effects andmanifestations of general properties of 1D systems (such as spatial heterogeneity that is associated with strong density fluctuations, the lack of phase transitions, the presence of frozen disorder, confinement, and blocked movement of nuclear particle by its neighbours) in nonequilibrium phenomena by considering the four examples. The anomalous transport in zeolite channels is considered. The mechanism of the transport may appear in carbon nanotubes and MOF structures, relaxation, mechanical properties, and stability of nonequilibrium states of free chains of metal atoms, non-Einstein atomic mobility in 1D atomic systems. Also we discuss atomic transport and separation of two-component mixture of atoms in a 1D system—a zeolitemembranewith subnanometer channels.We discuss the atomic transport and separation of two-component mixture of atoms in a 1D system—zeolite membrane with subnanometer channels. These phenomena are described by the response function method for nonequilibrium systems of arbitrary density that allows us to calculate the dynamic response function and the spectrum of relaxation of density fluctuations 1D atomic system.

According to modern concepts, the common feature of one-dimensional objects as statistical systems is spatially inhomogeneous state and density fluctuations.These fluctuations are caused by the crucial role of many-particle correlations in 1D systems.This feature is connected with the absence of phase transitions in 1D systems of particles of one kind with short-range potential of the interparticle interaction (see van Hove Theorem [28,29]).In this case, for certain values of the density the susceptibility of 1D system has a maximum but stays finite in contrast to 3D systems, where susceptibility becomes infinite in vapour-liquid and liquid-solid phase transitions [30].
The absence of phase transition means that stable nuclei of a new phase are not formed, and the one-dimensional clusters nuclei that are formed in a system have a finite lifetime and are the density fluctuations.1D systems are characterized by frozen disorder because the atomic particles cannot be exchanged with their neighbours and block the movement of each particle neighbors.For kinetic phenomena in atomic 1D confinement system, the role of confinement is not only to preserve the geometry of the system, but also the system relaxation by phonons.Relaxation of the system in confinement due to strong interparticle correlations depends on the ratio of the energy of particle interaction with confinement and interparticle interaction [31].
This paper is devoted to the analysis of correlation effects and manifestations of these general properties of nonequilibrium phenomena in the four examples.The anomalous transport in zeolite channels is considered.The mechanism of the transport may appear in carbon nanotubes and MOF structures, relaxation, mechanical properties, and stability of nonequilibrium states of free chains of metal atoms, non-Einstein atomic mobility in 1D atomic systems.Also we discuss atomic transport and separation of two-component mixture of atoms in a 1D system-a zeolite membrane with subnanometer channels.
These phenomena are described below using the response function method for nonequilibrium systems of arbitrary density [31].This method allows us to calculate the dynamic response function and the spectrum of relaxation of density fluctuations of 1D atomic system.The spectrum contains diffusion and hydrodynamic modes.In the latter case, the relaxation of the system and the nuclear transport occurs through the distribution of density waves and not by the diffusion of particles.This relaxation mechanism becomes the determining by increasing the density that allows one to explain the abnormally high rate of atomic transport in subnanometer channels.The change of relaxation mechanism by increase of the density is associated with the observed non-Einstein character of mobility of atomic particles (the law for the mean square displacement  2 = 2 is violated ( is the diffusion coefficient, and  is the time of observation) [32]).The found spectrum of relaxation allows to calculate the probability of decay of 1D free clusters of various lengths [30] and to explain the dependence of the stability of clusters of metal atoms on the chain atoms properties, the mechanical state of a cluster under the influence of an external force, and the number of particles in a stable cluster.These values were measured simultaneously in experiments [13].
The response function method was used to describe the practically important transport of two-component mixture in membranes with subnanometer channels.
It was shown that a phase transition in an inhomogeneous state may occur in 1D system of particles of two kinds with the increase of filling factor.Such a phase transition occurs due to an effective particle attraction and is characterised by the formation of clusters of a finite size [33].The amplitude of the fluctuations in the density corresponding to the formation of such a state is the order parameter.An equation for the order parameter for a mixture of particles in the channels is derived.In contrast to the well-known phenomenological Landau-Ginzburg equation [34], the coefficients of the equation are calculated for a known interparticle interaction potential depending on the density and temperature.The analysis of this equation allows us to describe the observed increase in the separation factor with increasing flux through 1D channels at their closing by a stable clusters with a certain concentration of the components and determine the dependence of the separation on the concentration, temperature, and density.

Molecular Transport in Subnanometer
Channels.Molecular transport in nanometer and subnanometer channels in porous bodies is currently drawing a great deal of interest from the standpoint of fundamental science [35][36][37][38][39] and because of the many applications of membrane and nanotechnologies in various fields ranging from nuclear power to ecology [40,41].Polycrystalline ceramic membranes consisting of complex oxides (zeolites), possessing subnanometer channels ranging in diameter from 0.3 to 1.4 nm, have been synthesized recently [40].On account of the high selectivity of molecular transport in zeolite membranes, as compared with the wellknown polymeric membranes [41][42][43], new technologies for separation, reprocessing, and utilization of materials are being intensively developed on the basis of zeolite membranes [41][42][43].
When the channel diameters in membranes decrease to the nanometer scale, molecular transport is determined by the Knudsen flow in the central zone of the channel, free of interaction of molecules with walls, and particle diffusion in the surface force field [39].In subnanometer channels, the interaction potentials between molecules and the opposite walls overlap, and molecular transport occurs under conditions of a constant interaction of molecules with a solid.Consequently, the diffusion coefficient in the limit of small filling numbers of a channel is determined by the relaxation of particles on phonons and surface defects [37].For molecules, with the exception of light particles (H, He), the channel walls are impenetrable, and consequently molecular transport is possible only along the channel axes.In this sense, it is different from diffusion in solids and may be assumed to be one-dimensional.In contrast to surface diffusion in channels with diameter  > 1nm, a fundamentally new property of molecular transport should appear in subnanometer channels.As the filling factor of a channel or the diameter () of molecules increases, so that  < 2, each molecule can block the motion of two molecules.Consequently, it can be expected that as the diameter of the molecules or the external gas pressure increases, the molecular flux in membranes with subnanometer channels should decrease.However, experiments have shown that the diffusion flux of a number of molecules (CH 4 , C 2 H 6 , and others) in a ZSM-5 zeolite membrane with channel diameters 0.54-0.57nm does not decrease but rather increases as the external gas pressure increases [43].It has also been established that the diffusion coefficient for these gases increases by more than an order of magnitude as the filling factor of the channels with molecules increases.It has been found that for a number of gases, the temperature dependence of the flux possesses a maximum and a minimum.
In this section, the transport of a single-component molecular gas (system of atomic particles) in subnanometerdiameter channels is investigated theoretically.

Sorption Isotherm in a 1D
Channel.Let us consider a surface in contact with an ideal single-component gas with temperature  and pressure .Let us assume that particles located on the surface do not interact with one another.We will also assume that the energy of the gas molecules on the surface is  0 .Let there be  adsorption centers on the surface.Then the average number  1 of particles adsorbed on the surface as a function of the gas pressure and temperature has the form (Langmuir isotherm) where  is the mass of a gas molecule, ℎ is Planck's constant, and  =  −1 is the reciprocal of the temperature.The method used to derive the Langmuir isotherms (1) admits wide extensions.Specifically, this method can be used to solve the problem of the filling factor of a cylindrical channel with diameter  taking into account the interparticle interaction in the channel.Let the channel diameter be comparable to the maximum diameter of a gas molecule.Let us consider the equilibrium of the gas with the surface on which  of the channels described previously emerge.Let  1 be the binding energy of a particle at the entrance to the channel.If the energy  1 is negative, then it is energetically favorable for the gas molecule to enter the channel.Let  be the total number of particles in a channel of length ,  the total number of particles in the channel and on the surface, and  0 the number of settling locations in the channel.Then the partition function for the grand canonical ensemble, taking into account the interaction of the gas particles in the channel, is ( Here  2 is the binding energy of the particles in the channel and  int () is the partition function, corresponding to taking account of the interaction of the gas particles in the channel.
Since  int depends on the number  of particles in the channel, it is impossible to calculate the partition function (2) in the grand canonical ensemble in the general case.However, in the problem of filling of a channel, the states for which the number of particles in the channel  ≫ 1 should make the main contribution to the partition function (2).Then the quantity  int () can be replaced by the partition function of interacting particles in the channel, calculated with the average number  of particles in the channel: The approximation (3) corresponds to the thermodynamic limit ( ≫ 1) for gas (system of atomic particles) molecules located in the channel.Using the relation (3), the partition function (2) can be easily calculated.Using the condition for equilibrium between the gas and the surface, we obtain for the average number  of particles in the channel Here  int (, ) is the free energy of interaction per atomic particle in the channel.Instead of the number  0 of settling locations in the channel, it is convenient to introduce the average distance between settling locations in the channel ( = / 0 ) and to replace  by the filling factor  ≡ / of particles in the channel.Then, we obtain from (4) an equation determining the pressure and temperature dependences of the filling factor  of the channel: The relations (4) show that the problem of obtaining the equation for the adsorption (, ), determining the filling factor at various pressures and temperatures of the gas above the surface, reduces to calculating the quantity  int (, ) ≡ − ln  int (, ), determining the "correction" to the pressure as a result of the interaction of the particles in the channel.Thus, to obtain the adsorption isotherms in the system under study, it is necessary to calculate the total free energy of the particles in a subnanometer channel taking into account their interaction with the channel wall and with one another.For this, we will consider  particles in a channel whose size is comparable to the average particle diameter.In the general case, the total potential energy of such a system can be written as Here (r  − r  ) is the potential energy of the pair interaction of the particles located at points with coordinates r  and r  ; (r  ) is the interaction energy between a particle located at the point r  and the channel walls.We will transform the expression for the pair interaction potential of the particles using the fact that the channel diameter is comparable to the particle diameter.Then the interaction potential depends only on the particle coordinates along the channel .We write the expression for (r  − r  ) in the form  (r  − r  ) ≡  (  −   ,   −   ,   −   ) ≃  (  −   , 0, 0) ≡ Φ (  −   ) .
This relation holds when the channel diameter is comparable to the particle diameter, and the pair interaction potential (r  − r  ) of the particles is not a long-range potential.
It is convenient to represent the interaction energy between a particle and the channel walls in the form  (r  ) =  (r   ,   ) .
Here r   is a dimensionless two-dimensional radius vector in a plane perpendicular to the channel.It is convenient to write the partition function of such a system in the form Here The integration over the particle coordinates in a plane perpendicular to the channel can be performed exactly: Here Ψ(  ) is the effective potential in which a particle moves in the channel as a result of the interaction of the particle with the walls: exp [−Ψ ()] ≡ ∫ exp [− (r  , )] r  .
Let us assume that the energy of a gas molecule interacting with the channel surface is constant and equal to  1 = Ψ( min ) = Ψ(0).Then The energy  1 physically corresponds to the binding energy of a particle in the surface potential Ψ().Thus, under the assumptions made previously, the gas in the channel can be assumed to be one-dimensional (the channel diameter is comparable to the maximum diameter of an atomic particle).
It is well known [44,45] that the partition function of a one-dimensional gas with an arbitrary interaction potential can be calculated exactly under certain assumptions, which are formulated in the following.Indeed, let us consider an equilibrium system of  particles in a channel which possesses only one degree of freedom per particle that are located in the segment [0, ] of the  axis.The total partition function of the system has the form where Here Φ() is the pair interaction potential of the particles in the channel.To calculate  int , the cofactor responsible for the partition function of an ideal gas must be eliminated from the partition function (14).In our problem, the cofactor corresponding to a gas of particles in a channel without a pair interaction potential between the particles is already included in the expression (2).Thus, since the partition function of a one-dimensional gas without an interaction has the form we obtain from ( 2) Thus, the problem has been reduced to calculating the configuration integral In studying the configuration space of a system of particles in a one-dimensional channel, it should be kept in mind that the quantities   are not independent.They are related as where  is the total length of the channel.Physically, this relation corresponds to the impossibility of particles penetrating one another and the condition of "blocking" of particles with diameter  in a channel whose diameter  < 2.An explicit expression for the configuration integral can be obtained if the explicit form of the pair interaction potential of the particles is known.We will consider a simple but nontrivial case: a system of hard spheres with diameter .The expression for the interatomic interaction energy in this case is We obtain for  int in the hard-sphere approximation Thus, in the hard-sphere approximation, the adsorption isotherm for a gas in a subnanometer channel has the form The relation (22) makes it possible to construct the dependence of the filling factor on the gas pressure for various molecules.It should be kept in mind that the sign of the energy  1 determines the possibility or impossibility of a particle entering a channel: for  1 > 0, it is energetically favorable for a gas particle to enter the channel for any external pressure greater than  0 .For  1 < 0, the particles must overcome a potential barrier to enter the channel.The configuration integral can be calculated exactly for an arbitrary interaction between the particles in the channels [45].This makes it possible to obtain an equation of state of a 1D gas for an arbitrary interaction.Thus, for an interparticle interaction potential of the form where  is the effective radius of attraction, the equation of state of a one-dimensional gas is In the limits  → 0 and  → 0, this equation becomes the equation of state of a system of hard spheres.The isotherms calculated numerically starting from (24) are presented in Figure 1.It is also evident in this figure that at high temperatures the gas behaves almost as an ideal gas ( 1D ∝ ).At low temperatures the isotherms seem to consist of two parts.For high density  ≃ const, which is typical for a condensed phase, whereas at pressure less than a characteristic value a gas-like phase for which  1D ∝  appears.
It can be shown that at a transition from one regime to another, the free energy of the system has no singularities.For this reason, there is no exact analogy with phase transitions.This assertion agrees completely with the Landau-van Hove Theorem [44], according to which any one-dimensional model of a gas with a finite interaction radius does not undergo phase transitions.On the other hand, it should be expected that for sufficiently low temperatures, as the filling factor  increases, the gas in the channel tends to form clusters whose size increases with the filling factor.The number and size of the clusters grow with the total number of particles in the system in a manner so as to ensure the existence of such clusters in the thermodynamic limit, where  → ∞,  → ∞, and / = const.This agrees completely with the formation of nuclei of a new phase in threedimensional systems with first-order phase transitions [34,46].On the other hand, in our one-dimensional system, in contrast to three-dimensional systems, a continuous behavior of thermodynamic quantities determined as derivatives of the free energy of the system with respect to the temperature, pressure, and filling factor should be expected in the entire range of variation of the thermodynamic parameters.For example, the specific heat calculated using the relation (24) will possess a maximum at the "critical" point.This fact has been observed experimentally [35], which attests to the applicability of the one-dimensional model employed previously.

Transport in a
Dense 1D System.The gas flux through a membrane is the main experimentally measurable parameter (Section 2.1.3).Particle transport is usually described using a relation between the outgoing particle flux and the parameters of the problem, such as, for example, the difference of the gas pressures on different sides of the membrane.This description is based on Fick's relation between the gas flux density J, the gradient of the gas concentration, and the diffusion coefficient D, using the conservation law for the number of particles: In the relations (25) and (26), J is the flux, D is the diffusion coefficient, and ∇ is the gradient of the particle concentration.The relation (25) is actually a definition of the diffusion coefficient D. Thus, the problem of determining the particle flux through a channel reduces to calculating the diffusion coefficient.A more general method of describing transport is to calculate the characteristic relaxation time of density fluctuations arising or specially created in the system under study [47].Formally, the problem of determining this time reduces to calculating the characteristic frequency or spectrum (k) for the system under study.Thus, when the particle density is low and there is no interparticle interaction, an explicit expression for the spectrum (k) can be easily obtained from the relation (25).Indeed, let  0 be the average particle concentration in the system.From (25), we obtain an equation determining the dynamics of the density fluctuations (r, ) =  −  0 : Switching to a Fourier representation of the density fluctuations we obtain an expression determining the reciprocal of the relaxation time of the density fluctuations  in the case at hand: It is obvious that a spectrum of this type is characteristic for systems where relaxation to equilibrium occurs by diffusion.
Let us consider a system of interacting particles with arbitrary density.Assume that we have been able to calculate the characteristic relaxation time  −1 of density fluctuations in the system.If the spectrum of relaxation times in such a system has the form then it is natural to consider the quantity (k,  0 ) to be the diffusion coefficient in the system, provided that in the limit of low densities  0 the quantity (k,  0 ) becomes the diffusion coefficient for a low-density system of particles: lim This definition of the diffusion coefficient means that the system possesses a diffusion characteristic mode.Proceeding from the definition (30), to calculate the diffusion coefficient (k,  0 ), it is necessary to know the characteristic modes of a weakly nonequilibrium system.If these characteristic modes have the form (30), then (k,  0 ) will be the diffusion coefficient of the system.It is obvious that the diffusion coefficient in the general case is nonlocal and is a functional of the density of the diffusing particles.The expression for the flux of diffusing particles will have a form generalizing equation (25): The equation, obtained using this relation, for the particle number density in the system is in general a nonlinear integrodifferential equation.When studying the diffusion of a single-component gas in subnanometer channels, it should be kept in mind that the pair interaction potential between the particles is of short range, and no phase transitions occur in the system of particles in the channel (see the preceding section).In this case, it can be assumed that the diffusion coefficient (r, r  , (r, )) is local: The corresponding expression for the flux of diffusing particles and the equation for their density become For low densities, the relation ( 34) becomes (25), which is the first term in the expansion of the flux in odd powers of the density gradient.Thus, for systems where the relaxation spectrum of the fluctuations is of the form (30), the diffusion coefficient can be determined starting from the explicit form of this spectrum.The expression obtained for the diffusion coefficient in this case is a generalization of the expressions obtained for the diffusion coefficient in various models, specifically, in the Maxwell-Stefan model [43].
It is convenient to calculate the relaxation spectrum using the response functions [31,48,49].Let us formulate the basis of this formalism.Consider the system of the atomic particles, which can interact with each other as well as the environment such as liquid, solid, or surface.Let us show, following [49], that many-particles distribution functions of such a system are the density functionals.Assume that each particle at the coordinate r interacts with time-dependent external field V(r, ).In order to describe such a system, one can use Liouville equation for many-particle distribution function ( 1 ,  2 , . . .,   , ), where   = (r  , p  ) are the coordinates and th momentum of the particles.Solving the Liouville equation with the fixed initial distribution function and for the different external potentials (r 1 , r 2 , . . ., ) = ∑  V(r  , ), one can obtain one-to-one conformance between ( 1 ,  2 , . . .,   , ) and (r 1 , r 2 , . . ., ).On the other hand, introducing the density of the particles where Γ is the phase space element, and n is the density operator, one can obtain one-to-one conformance between the external field (r 1 , r 2 , . . ., ) and the density (r, ).
From this it directly follows that many-particles distribution functions are density functionals: Following [49], let us write the action for the Liouville equation in the form where  is the Lagrangian of the system, which is equal to the difference between the kinetic and potential energy of the particles.It follows from ( 36) and (37) that the action  is the density functional: Functional Δ[] has an extrema for the equilibrium state of the system [48,49].Due to the linearity of the (38) by the potential energy, system of the particles interacting with the external field V(r, ) corresponds to the functional Δ[, V], which can be expressed as Functional (39) has an extrema for the stationary or quasistationary (r, ) =  0 (r, ) ≡  (0) state of the system [48,49] and represents the time average of the free energy functional [(r, ), ]: The functional Δ[] was for the first time introduced in such a form in the papers [31,48].Equation (40) determines the density of the atomic particles (r, ).Let us express the free energy using the response functions defined as follows.
Consider the system in the external field V ext (r, ) (external field amplitude).The response function (r, r  , ,   ) of such a system couples the density fluctuation (r, ) with the external field V ext (r, ), which generate this fluctuation [31,48] at  → 0: It is obvious that in the case of the linear response and  → 0, the response function (r, r  , ,   ) is independent on .Let us rewrite the definition (42) in the operator form Such a form (43) will be used further in the paper.In the presence of the external field, the functional Δ[, V] at  → 0 is determined by the relation (39): V ext (r, )  (r, ) r . ( Here Δ[] is the functional at  = 0, determined by the formula (38).In order to obtain linear response function, let us use the relation (39).One can rewrite the solution of the (39) at  → 0 with the sum (r, ) =  (0) (r, ) + (r, ).In this case, one can obtain, in the first order and with respect to (39), or At  → 0, the quantity is independent on .Comparison of the relations ( 43) and (46) gives the following form of the response function (r, r  , ,   ): It follows from the functional Δ[] definition that it can be expressed as Here Δ 0 [] corresponds to the functional of the particles without their interaction with each other, and Δ int [] corresponds to the functional (38), which represents particleparticle interactions.Two times variation of the (49), with respect to (48), gives the following: where is the response function of the noninteracting particles It is shown in the paper [31,48] that the function (r, r  , ,   ) is connected to the effective lagged interaction between the particles.Multiplying the relation (45) on the  in the left and on the  0 in the right, one can obtain the equation on the response function.In the case of the homogeneous quasi-stationary state of the atomic system, the functions  0 (r, r  , ,   ), (r, r  , ,   ) j (r, r  , ,   ) depend only on the differences r − r  j  −   , so (50) can be rewritten as where (k, ),  0 (k, ) j (k, ) are the Fourier transform of the response functions (r − r  ,  −   ),  0 (r − r  ,  −   ), and (r − r  ,  −   ), respectively; for example,  0 (r − r  ,  −   ) = ∫ k   −(−  )+k(r−r  )  0 (k, ) .(54) Equation (53) allows us (when the  function is known) to determine the response function of the system of interacting particles and relaxation spectrum of the system [31,48].
In order to determine the relaxation spectrum, let us note, that the expression  −1  = V ext follows directly from (43).Due to the response function definition (as the response of the system on the arbitrary small external field), one can write  −1  = 0.This equation has the nontrivial solutions  ̸ = 0 only if  −1 = 0. So, we can write that the relaxation spectrum of the system (k) is determined by the relation [31,48] We will use the response function formalism to calculate the response function and the diffusion coefficient of a dense system of atomic particles in a subnanometer channel.We will calculate first the response function  0 (, ) and the relaxation spectrum of such a gas when there is no interaction between the particles.We write the diffusion equation for the filling factor  of the channel with molecules in the presence of a weak perturbing external field V ext (, ) in the form Here  0 is the diffusion coefficient for noninteracting particles.Let  0 be the equilibrium filling factor of the channel.We will seek the solution of (56) in the form Then, we obtain from (56), assuming the external field  ext to be weak, an equation for (, ): In deriving (58), only first-order infinitesimals in  and V ext were retained in all terms.We will seek the solution of (58) in the form of a Fourier integral Substituting the expressions (59) into (58) and solving the resulting linear equation for (, ), we obtain  (, ) =  0 (, ) V ext (, ) . ( Here is the response function of a gas of noninteracting diffusing particles.Indeed, as  → 0, the function  0 (, ) reduces to the well-known response function of an ideal equilibrium gas with density  0 at temperature  [38,47]: The relaxation spectrum of such a gas is determined from the relation (55) and has the form In accordance with the definition (26), the quantity  0 is the diffusion coefficient.The relations obtained make it possible to calculate the relaxation spectrum () by solving (54).Thus, we find for the response function of a one-dimensional system Using the relations ( 61) and ( 64), we obtain an equation determining the relaxation spectrum of the system under study: Equation ( 65) can be solved in a general form in the quasistatic case  → 0. Indeed, in this case, we find from ( 65) It is convenient to rewrite the relation (67), introducing the pair distribution function Then we obtain from relations (53) and ( 61) Substituting the expression (69) into (67) gives a relation determining the relaxation spectrum of a dense gas in a onedimensional channel: We will now calculate the diffusion coefficient of a gas in a one-dimensional channel taking into account the interaction between hard-sphere particles.In this case, the correlation function  2 () =  2 0 [1 + ]()] has been calculated exactly [45] with an arbitrary filling factor of the channel: Using the relations ( 70) and ( 71), we find It follows from (72) that the general equation describing particle transport in a dense 1D system, in accordance with (34), has the form This equation can be obtained from (72) if the characteristic local equilibration times in a dense system of particles in a channel are short compared with the characteristic propagation times of disturbances along such a system.For the cases considered in the following, this assumption obviously holds.We will consider solutions of (73) in the case where complete equilibrium has been established in the channel and the filling factor  =  0 (, ), with the function  0 (, ) being determined by the relation (23).Transport in a channel in this case will be determined by the dynamics of the motion of disturbances of the filling factor.The solution of (73) in this case should be sought in the form  =  0 (, ) + (, ).It follows from (73) that the relaxation spectrum of disturbances in this case has the form (72), where  =  0 (, ).Thus, taking account of the interaction of hard-sphere particles in a channel does not change the character of the relaxation of its weakly nonequilibrium state: the relaxation spectrum (72) remains a diffusion spectrum with diffusion coefficient ( 0 ).However, it is important to underscore that transport in the channel in this case is a collective effect and proceeds via transport of disturbances of the equilibrium density  0 (, ).
The diffusion coefficient ( 0 ) in this case is the diffusion coefficient of disturbances of the equilibrium density.It is convenient to rewrite the relation (72) for the diffusion coefficient in a different form.Using the Arrhenius character of the diffusion coefficient for noninteracting particles ( D0 is proportional to the product of the squared lattice constant of the wall material of the channel and the relaxation frequency of particles in a channel on defects and lattice phonons and is calculated in [37], and  ≈ Ψ( max ) (13) is the activation energy of diffusing of noninteracting particles in the channel), we will rewrite (72) in the form This way of writing the diffusion coefficient for a gas of interacting particles in a channel makes it possible to give a physical interpretation for the change in the diffusion coefficient accompanying a change in the filling factor  0 .The relation (75) shows that the interaction of particles in a channel decreases the activation energy  of the motion of particles, even in the hard-sphere model.For an interaction between hard-sphere particles such that there is no direct attraction between the particles, the diffusion activation energy decreases with increasing filling factor  0 as a result of the effective interaction (see the following).Physically, this corresponds to a change because of the presence of another particle in a neighboring potential well, in the parameters of the potential in which a gas particle moves.Since the gas particles are assumed to be indistinguishable, diffusion with the filling factor  0 ∼ 1 (when the effective diffusion activation energy Ẽ( 0 ) becomes comparable to the temperature of the system) can be interpreted as a transfer of "excitation" of the density along a chain of close gas particles.
It is obvious that the motion of such an "excitation" will occur with substantial velocities.This leads to a large increase in the diffusion coefficient when  0 ∼ 1.For  = 0, the possibility of such a diffusion mechanism is played out in [36].We note that for all gases investigated, the barrier  is different from 0. This is indicated by the presence of a maximum in the temperature dependence of the flux in a zeolite membrane [43].It can be assumed that the increase in the diffusion coefficient in a channel observed with filling factors  0 ≥ 0.5 is related to the formation of clusters in the channel whose sizes increase with the filling factor.The transport of gas in the channel containing such clusters is determined by the motion of "excitation" in a finite-size cluster.The formation of clusters in a gas consisting of particles with a hard-sphere pair interaction potential is related to the well-known [46] appearance of an effective attraction between such particles.This effect is manifested in the appearance of a maximum in the correlation function  2 (), describing the probability of finding the "first" particle at a distance  from the "second" particle.
The relaxation spectra of the system for an arbitrary wave vector  can be found from (70) and (71).It is found that the spectrum () contains a real part, which corresponds to transport of an "excitation" in the system.The imaginary part of the spectrum obtained has a minimum at a definite value of the wave vector  =  min , which depends on the filling factor  0 .It is natural to interpret the value  = 2/ min ( 0 ) as the characteristic size of a cluster with a given filling factor and the value  = Im  −1 ( min ( 0 )) as its characteristic lifetime.As the filling factor increases, the value of the effective attraction between the particles increases [46]; this could lead to the appearance of clusters consisting of two or more particles.However, it is found that as a result of the one-dimensionality of the channel, the lifetime of the clusters that are formed is finite, which corresponds to the absence, as noted previously, of a phase transition in such systems.The dependences of the lifetime and size of the clusters formed on the filling factor  0 are presented in Figure 2.
It is evident in the figure that increasing the filling factor increases the lifetime and the size of the clusters formed.In finite-size channels, for a definite value of the filling factor   0 =   , a cluster equal in size to the length of the channel is formed.It is obvious that the diffusion coefficient in such a channel increases without bound as  0 →   , even though the lifetime of such a cluster is finite.If the range of relaxation frequencies () is known, then an expression that defines the relaxation of the Fourier components of the particle density in the channel (, ) in the case of arbitrary density (, ) is given by ṅ (, ) =  ()  (, ) . ( From (76), we can get that the equation for the amplitude of the fluctuations can be written as  ṅ (, ) =  ()  (, ) .
Equation ( 77) describes the relaxation of  component of the density fluctuation at arbitrary wave vector.In particular, at  ̸ = 0, this equation allows to describe the relaxation of density fluctuations and propagation of disturbances on the cluster of finite size in the case of the formation of the cluster.At  → 0, the equation describes the relaxation of density fluctuations on large spatial scales.This value is related to the macroscopic flux.Let us write the equation of continuity to calculate the flux : After applying the Fourier transform, (78) takes the form Using ( 76) and (79), one can get for ()  (, ) =  (, )  ()  .
Using the relation (70) we find that ]() and therefore the spectrum contains both real and imaginary parts, corresponding to the "hydrodynamic" (real part of ()) and diffusive (imaginary part of ()) transport of particles.Therefore, (80) for the flux (, ) contains terms relating to hydrodynamic transport of particles in a cluster and diffusive transport of particles in the space between clusters.We note that the response function method can be used to take into account the interaction of particles in a channel for a potential that is different from a hard-spheres interaction, similarly to the way this was done previously.
Analysis of the experimental data showed that for high filling factors  0 ∼ 1, the dependence of the diffusion coefficient ( 0 ), calculated from the relation (72), leads to a discrepancy between theory and experiment (Section 2.1.3).This discrepancy could be due, in our opinion, to the abovedescribed influence of a finite channel length, the asphericity of the gas molecules, and the possible breakdown of one dimensionality of the problem.These effects lead to the following dependence of the diffusion coefficient on the filling factor: Here  is a coefficient that takes into account the finite size of the channel, the possible deviation of the channel from one-dimensionality, and the asphericity of the molecules.The relations ( 72), (81), and ( 52) make it possible to describe transport in a subnanometer channel for different values of parameters such as the gas pressure and the pressure difference outside the membrane, the temperature of the membrane, and the type of gas molecules.As an example, we present the dependence, calculated according to ( 72), (81), and ( 23), of the diffusion coefficient in a subnanometer channel on the relative temperature / 0 for various relative pressures / 0 (Figure 3).It is evident in the figure that the character of the temperature dependence of the relative diffusion coefficient in subnanometer channels is substantially different at different pressures .For example, if the pressure is high (/ 0 > 1), the dependence of / 0 on / 0 has a pronounced maximum at temperatures  ∼ 0.5 0 (curve 1).As pressure decreases, the magnitude of this maximum decreases and the maximum itself shifts into the region of lower temperatures (curves 2, 3).A further decrease of pressure to values / 0 < 1 results in the appearance of a minimum in the dependence of / 0 on / 0 in the temperature range  ∼ 0.5 0 .The maximum occurs at  ∼ 0.3 0 .Since for / 0 < 1, the flux depends on  and the average distance between particles in the channel is much greater than the diameter of the particles, the increase in the diffusion coefficient for  ∼ 0.5 0 is due to the temperature dependence of the diffusion coefficient  0 on individual molecules in the channel.At high temperatures, all curves saturate at a value corresponding to / 0 = 1.In the next section, the theoretical laws obtained previously will be compared with the experimental data.

Analysis of the Experimental Data: Comparison of Theory with
Experiment.Quite extensive experimental investigations of the sorption and transport properties of a series of organic molecules (CH 4 , C 2 H 6 , -C 4 H 8 , -C 4 H 8 , and others) and the inert gases Ar, Ne, and Kr in ZSM-5 zeolite membranes (MFI, silicate-1) have now been performed.The results are presented in [42,50,51].According to [52], ZSM-5 zeolite membranes have a complicated chemical and crystalline structure.The chemical structure of zeolite membranes is given by the formula For small values of , zeolite membranes of this type are called silicate-1 membranes.The crystalline structure of ZSM-5 membranes has been studied quite well.It consists of a 3D structure, consisting of sinusoidal channels with a circular cross section (0.54 ± 0.02 nm) parallel to the  axis [53], which intersect straight channels with an elliptical cross section (0.57-0.58) × (0.51 − 09.52) nm 2 parallel to the  axis [010] [50,54].Cavities ∼0.9 nm in size form at the intersection of the channels.The sorption capacity of ZSM-5 zeolite is determined by the number of sorbed (entering the channels) molecules per cell of a crystal.According to [50,55], the cell parameters are  = 2.007 nm,  = 1.992 nm, and  = 1.342 nm.A single cell is a structure consisting of four segments of 0.46 nm linear channels, four segments of 0.66 nm sinusoidal channels, and four 0.54 nm intersections.
Depending on the structure of a molecule, sorption of one or two molecules per intersection is possible.In [50], the sorption capacity of ZSM-5 was measured for a number of molecules.It was concluded that the molecules CH For nitrogen-type molecules, two molecules can be arranged at an intersection, and therefore molecules can change places at an intersection (absence of blocking).Nonetheless, since there are significantly more nitrogen molecules in a channel (24 molecules per unit cell [50]), the one-dimensional model can also be used for this gas.Analysis of the experimental data presented [43,[55][56][57][58] shows that all gases investigated can be conditionally divided into two groups.The first group of gases, containing, specifically, Ar, Kr, Ne, and CH 4 , is characterized by a linear pressure dependence of the filling factor ; that is, Henry's law holds for them [43].For these gases, the pressure dependence of the flux is nearly linear, and the diffusion coefficient is essentially independent of the filling factor [55].
The temperature dependence of the flux is characterized by the presence of a minimum at temperatures  ∼ 400 K (e.g., for Ar and Kr), while a maximum is not observed in the diffusion coefficient [57,58].Here neon, whose diffusion coefficient increases with temperature [58], is an exception.
The second group of gases, containing, specifically, i-C 4 H 10 , C 2 H 4 , C 2 H 6 , and C 4 H 8 , is characterized by a dependence of the filling factor  on the pressure  in the form of a curve that saturates [43].For these gases, the pressure dependence of the flux is likewise characterized by a curve that emerges at a value that is pressure independent [55].For this group of gases, the diffusion coefficient depends strongly on the filling factor [56].The temperature dependence of the flux is characterized by the existence of a maximum and minimum at high ( > 500 K) temperature [57].The model, constructed in this paper, of the sorption and transport properties makes it possible to describe these dependences of transport in onedimensional channels on the basis of general assumptions.Figure 4 shows the dependence of the diffusion coefficient  on the filling factor  (solid line), calculated using (81) and the experimental data obtained in [57] for various gases.We note that the dependence (81) presented in Figure 4 takes into account the asphericity of the molecules, equal to  = 0.86.For inert gases and methane (first group of gases), the diffusion coefficient measured in the range of filling factors  0 < 0.5 depends weakly on  0 .A large increase in the diffusion coefficient at  0 ∼ 1 is observed for gases of the type C 2 H 6 and -C 4 H 8 (second group of gases).This behavior of the diffusion coefficient () can be satisfactorily described on the basis of the model proposed, which takes account of the pair interaction of particles with a very simple form of the interaction potential for hard spheres.It is obvious that this coefficient is different for different molecules.
The physical mechanism leading to an increase of the diffusion coefficient consists in a decrease of the activation barriers for diffusion, which is due to an interaction of the particles for high filling factors.Indeed, it follows from the relation (75) for the effective diffusion activation energy Ẽ( 0 ) that in this case, when the filling factor is low,  0 ≪ 1, gas diffusion in the channel can be studied in the singleparticle approximation neglecting the pair interaction of the gas particles with one another.As the filling factor increases, the diffusion activation energy decreases as a result of the effective attraction of the gas particles to one another.This can be interpreted as the formation of widely separated clusters.Diffusion in each cluster can be treated as a transfer of "excitation" along a chain of close gas particles.A further increase of the filling factor increases the cluster lifetime and decreases the average intercluster distance and, in consequence, increases the diffusion coefficient.For  0 ∼ 1, the diffusion process can be regarded as a motion of this "excitation" along the entire channel, when the arrival of a particle at the channel entrance results in the end-most particle leaving the channel.
Figure 5(a) shows the temperature dependences of the gas flux through a membrane for ethane and argon, calculated using the formulas ( 5) and ( 81) and a relation following from (81):  ∝ ( 0 (, )) 0 (, ).These dependences were compared with the experimental values obtained in [42].It is evident that the theoretical and experimental results agree satisfactorily with one another.The temperature dependence can also be understood on the basis of the model developed in the present work.Indeed, it follows from the relation ( 5) that at low temperatures ( → 0), the filling factor  0 ∼ 1.
The effective diffusion activation energy Ẽ( 0 ) is  0 at zero temperature (56).The diffusion coefficient of the gas particles in the channel (81) approaches zero as a result of the "freezing out" of the thermal motion of the particles in the potential field of the channel surface.Increasing the temperature decreases Ẽ( 0 ) without appreciably changing the filling factor.The value of Ẽ( 0 ) reaches a minimum at a certain temperature  0 .It is obvious that the diffusion coefficient reaches its maximum value precisely at this temperature.It follows from ( 54) and ( 75) that a further increase of temperature decreases the filling factor  0 , and this is accompanied by an increase in the effective diffusion activation energy right up to the value  0 which is attained for  0 (, ) ≈  0 ≪ 1 at some temperature  =  1 ().The temperature  1 corresponds to a minimum in the function ().A further increase of the diffusion coefficient with increasing temperature corresponds to the variation of the coefficient according to the Arrhenius law  ∝ exp(− 0 /).Figure 5(b) shows the temperature dependence of the flux for argon, calculated using (81) and the data from [57].We note that the theory developed in the present paper predicts that the presence of a minimum at  ∼ 400 K should be accompanied by the appearance of a maximum at  ∼ 200 K.In this connection, for further elaboration of the theory, it is of interest to make an experimental search for a maximum in the temperature dependence of the flux for argon near  ∼ 200 K.Thus, two types of behavior of particles in a channel, which have different diffusion mechanisms, can be distinguished in the pressure and temperature ranges investigated.For the first group, this is the diffusion of single particles in a channel, and for the second group it is diffusion as a result of a collective interaction of particles in a completely filled channel.Such behavior of single-component gases is a direct consequence of the real one-dimensionality of subnanometer channels, where molecules cannot change places with one another.The proposed model is based on two basic assumptions: (1) the pair interaction between gas particles plays the decisive role in the description of the state and transport phenomena in a one-component gas in subnanometer channels, and (2) as a result of the fact that the channel diameter is comparable to the diameter of the gas molecules, such channels can be assumed to be one dimensional.The one-dimensionality of the system studied is actually due to the fact that the gas molecules do not have classical transverse degrees of freedom.This difference from the conventional systems plays a fundamental role in the analysis of all phenomena described above.
that an experimental consequence of the existence of such solutions will be the appearance of substantial fluctuations of the gas flux through a membrane.The proposed model satisfactorily describes the equilibrium properties of a gas in a channel (specifically, the filling factor) and the weakly nonequilibrium properties of the gas, such as the diffusion coefficient and flux.It should be noted that a satisfactory description of the equilibrium properties, for example, the adsorption isotherm, has been obtained previously on the basis of a number of phenomenological models [42].A qualitative analysis of the maximum (but not minimum) in the temperature dependence of the relative diffusion coefficient is presented in [43].A model describing the dependence of the diffusion coefficient on the filling factor on the basis of a oneparticle approximation by introducing into the diffusion rate jumps of individual particles a phenomenological correction factor  = (1 −  0 ) −1 has been proposed in [36].However, even though each of the phenomena indicated previously can be described on the basis of individual models, the approximations used in so doing either are not physically substantiated or completely contradict one another.

Particle Mobility and Diffusion in 1D
System.Although, in many cases, the parameters of 1D systems can be calculated exactly, their complete theoretical description has not yet been accomplished.1D systems belong to the class of strongly fluctuating systems [59]; thus, in contrast to systems of higher dimensionalities, they cannot unambiguously be described only by mean values such as mean density and compressibility [59].The transport properties of 1D systems can also be rather specific (see [43,57,60] and Section 2.1).
Recent ab initio calculations [61] have shown that water flux through hydrophobic nanotubes is more than one order of magnitude higher than the estimated ordinary diffusion flux.The permeation experiments on zeolite membranes [57] have shown that for gases, whose particle diameter is larger than half the zeolite channel diameter, the increase of the diffusion coefficient with density is more than one order of magnitude larger than would be expected for one-particle diffusion.On the other hand, there are also experimental results, for example, obtained using nuclear magnetic resonance (NMR) and quasi-elastic neutron scattering (QENS) techniques, indicating that the particle mobility in 1D systems can decrease and the mean square displacement (MSD) can become proportional not to time, as it is in the case of a normal (Einstein) diffusion, but to the square root of time, ⟨ 2 ⟩ ∼  1/2 (so-called single-file diffusion mode (SFD)) [43,[62][63][64][65][66][67].Moreover, the experiments [64,66] show that both the Einstein and SFD diffusion modes can exist for the same particles but at different densities and observation times.Furthermore, for certain investigated gases, a density increase leads to a transition from the Einstein diffusion to the SFD mode [64].Thus some experiments suggest that the particle mobility in 1D systems must decrease with increasing density, whereas others show that diffusion accelerates as density becomes larger.An explanation for this contradiction is absent at present.Diffusion in 1D systems has been studied in numerous theoretical works [36,57,[67][68][69][70][71][72][73].The dependence ⟨ 2 ⟩ ∼  1/2 has been obtained for a 1D system within the screened-gas model, which partly takes into account correlations between the particles [68][69][70][71][72][73].It has been suggested [70][71][72][73] that the appearance of the two diffusion modes, characterized by the different time dependences of MSD, can be explained by a transition from the one-particle diffusion mode, occurring at short times, to SFD appearing at longer times when the particle mobility is limited by its neighbors.In this case at  → ∞ the dependence of MSD is always ⟨ 2 ⟩ ∼  1/2 .Although this approach can describe the experimentally observed dependences of MSD, corresponding to SFD, it fails to explain the increase of the diffusion coefficient and transport acceleration observed in dense 1D systems [43].The physical reason for this failing is that the proposed model [60] considered the system of individual tagged particles, whereas in experiments one does not follow the diffusion of a particular particle but rather deals with the system of many of identical particles.We propose a description for the particle diffusion [57] and mobility [43,[62][63][64][65] in 1D systems of identical particles with blocking effect (BE) when the motion of a particle is blocked by other particles or clusters.According to our model, different mobility and diffusion transport mechanisms can be realized in a 1D system depending on its density.At low densities, an ordinary diffusion of single particles takes place.As density increases, a formation of 1D clusters of a finite size and lifetime takes place.In this case, there are two competing mechanisms prevalent in the system: a barrier-free propagation of density perturbations ("collective diffusion" (CD)) [60] and a blocking effect.When BE dominates, the SFD mode is realized, whereas when the CD mechanism prevails, we predict an essential transport acceleration, as compared to the oneparticle diffusion mode.In this case, MSD is Einstein-like with the diffusion coefficient increasing with density.The proposed model consequently takes into account all the timespace correlations occurring in a 1D system composed of identical particles with an arbitrary interaction potential.This model is able to describe the experimentally observed MSD dependencies at different densities as well as the transport acceleration with increasing density of a 1D system.To describe the transport properties of a 1D system, we apply the response function method (see Section 2.1 and [31,32,48]).This approach allows one to consistently take into account many-body correlations and obtain the diffusion coefficient and MSD for a 1D system of arbitrary density.Within response function method, we calculate the response function (, ), which determines the Fourier image of the density fluctuation (, ) occurring in the 1D system in the external field of frequency  and wave number .The relaxation spectrum of a density fluctuation (()) can be calculated from the relation  −1 (, ) = 0. Note that in response function method, (, ) can be interpreted as a distribution function if  ≫  2 /, where  is the particle diameter and  is the one-particle diffusion coefficient [48].A density fluctuation (, ) in a 1D system is determined by the relation (87), and the MSD can be expressed as It is evident from ( 87) and ( 83) that MSD is determined by the relaxation spectrum ().The evaluation of () is reduced to the solution of Dyson equation for the response function (55) and for the 1D system takes the form (70), where ](, ) is the Fourier image of the pair correlation function of interacting particles.The function ](, ) can be obtained for a 1D system with an arbitrary interaction potential Φ() [45]: Here  is the temperature.The last expression in (84) helps us to obtain the dependence of pressure on the filling factor (()).It follows from ( 70) and ( 84) that at  → 0, the relaxation spectrum () describes the one-particle diffusion () = − 0  2 .For  ≫ 1, the Fourier image of the pair correlation function ](, ) approaches zero for any arbitrary short-range potential Φ().This corresponds to the absence of correlations between the particles at the length scale  <  0 , where  0 is the average particle path between two collisions.In this case, the relaxation spectrum () also describes the one-particle diffusion mode () = − 0  2 .The analysis of the relaxation spectrum performed in Section 2.1 for the hard-core potential has shown that due to collective effects and density fluctuations, the 1D system can transform into a nonhomogeneous state as density increases.This state is characterized by the formation of 1D clusters with a finite lifetime (Section 2.3) (Section 2.1).In this case at high density ( → 1), long distances ( → 0), and  → ∞ (see ( 70), ( 84)) the spectrum becomes diffusion-like with the diffusion coefficient increasing with density (72) that corresponds to the CD mode existing in a nonhomogeneous, clustered 1D system [62].In order to describe the relaxation of density fluctuations at times  ≫  2 / 0 (that corresponds to the longtime limit  → ∞) for an arbitrary filling factor , one should do series expansions of ( 70) and ( 84) with respect to the parameter  ≪ 1 and keep the expansion terms following the first one.It is convenient to rewrite this expansion with respect to parameter  0 .Having done this, the relaxation spectrum for  ≪ 1/ can be written as where the expansion coefficients are () =  0 (0, ),  1 () = ]  (0, ) 2 (0, ), and () =  2 ((1/2)]  (0, ) 2 (0, ) − ] 2 (0, ) 3 (0, )).The values with apostrophes represent the corresponding derivatives of ](, ) with respect to wave number , calculated at  = 0. Equations ( 83) and ( 85) can be used to evaluate the time dependence of ⟨ 2 ⟩ at any length scale, starting from a diffusion path of just few particle jumps up to those comparable to the total length of the channel.Note that the real part of the spectrum Re (), which at  → 1 corresponds to the sound mode, describes the collective effect of density perturbation propagation ((, ) ∼ ∫  −()   = ∫  1 (, ) 2 (, )  ,  1 ∼   Re () ), and its imaginary part Im () determines the diffusion mode for the relaxation of a density fluctuation taking into account the length-time correlations between the particles ( 2 ∼  − Im () ).Equations ( 87), (83), and (85) allow us to study the mobility and distribution function of the particles in a 1D system for different length/time scales, filling factors and interaction potentials Φ().Figure 6 presents the dependences of the distribution function on the dimensionless coordinate  = / 0 , calculated for different dimensionless times   = / ( ∼  2 0 / 0 ) using the squarewell interaction potential with the well depth  ∼  ( = / ∼ 1) and width .
It is important to mention that the obtained results do not depend on the specific form of the pair interaction potential Φ(); we have obtained similar results using rectangular potentials, the Lennard-Jones and Morse potentials.We also note that the shape of the distribution function is independent on density as both dimensionless coordinate and time are determined by the filling factor.The distribution functions at different time scales clearly differ from each other.The distribution is a Gauss-like function at short times ( 2 ≪   < 1,  ≪ || < 1) that corresponds to the one-particle diffusion (Figure 6).At   ∼ 1, || ∼ 1, the distribution function varies faster than the Gauss distribution (Figure 6) that corresponds to limited mobility of a particle due to its blocking by other particles.At longer times (  ≫ 1, || ≫ 1), the distribution function becomes Gauss-like again (Figure 6) as the transport in a dense 1D system is due to the CD mechanism [60].In this case, the relaxation spectrum is diffusion-like with a density-dependent diffusion coefficient MSDs calculated at different time and length scales (Figure 7) using the distribution functions (Figure 6) noticeably differ from each other confirming the existence of several diffusion modes in the system.At short length-time scales ( < , ⟨ 2 ⟩ ≪  2 0 ), one should observe the one-particle diffusion mode with ⟨ 2 ⟩ ∼  (Figure 7(a)).When ⟨ 2 ⟩ is comparable to  2 0 and time  ∼ , the correlations between the particles become essential that lead to the appearance of blocking effect.In this case, MSD has a non-Einstein character (Figure 7(b)).At longer times  ≫ , ⟨ 2 ⟩ ≫  2 0 , the particle mobility is mostly determined by the second term in (85).This leads to the reappearance of the Einstein-type 110 dependence of MSD, ⟨ 2 ⟩ = 2(), where () is an effective diffusion coefficient.Notice that the characteristic times at which different diffusion modes can be observed depend on the density of the system.Indeed, the average distance between the particles and, consequently, the time  ∼ 1/ 2 decrease with increasing density.
Further we compare our results to those of the QENS, NMR [43,[62][63][64][65], and permeation [57] experiments as well as to the results of the molecular dynamics (MD) calculations [71] summarized in Table 1.Notice that such a system can be considered as one dimensional with the presence of BE as the diameters of the experimentally investigated molecules are less than the doubled diameter of the membrane channel.Our estimates for zeolite ZSM-48 obtained using the data from [71] show that for the low filling factor ( = 0.11) and times ∼ 10 ns, the normal diffusion mode, ⟨ 2 ⟩ ∼ , should occur while at higher density ( = 0.48) the dependence ⟨ 2 ⟩ ∼  1/2 should be observed, which is in a perfect agreement with the experimental data (Table 1).In Table 1, experimental and MD results are presented;  and  are the particle and membrane channel diameters, respectively,  is a temperature, and  is a degree of filling of the channel.The fifth column represents the data on the motion mode of a particle in a channel ("SFD" means  2 ∼   ,  < 1, "normal" corresponds to  2 ∼ ).The observation time of corresponding modes is represented in the sixth column.
The experimental results obtained for ethane molecules in zeolite AlPO 4 -5 and for methane molecules in zeolite ZSM-48 (Table 1) are successfully described by our model.Indeed, using (85), one can write MSD for the SFD mode at () > 0 as Equation (86) with () = 1 is identical to the relation derived in [62] with the accuracy of the numeric coefficient √ 32.Note that () describes the correlations effects in a 1D system of interacting particles.Using (86) and the parameters for AlPO 4 -5 [62], we obtain ⟨ 2 ⟩ ∼ 3 m for times ∼100 ms that is in good agreement with the values reported in [62,64].At short times  ≤  2 0 / ∼ 10 −8 s, the main contribution comes from the first term in (85), and, in this case, the normal diffusion mode is realized.Our results obtained using (83) and ( 85) at  = 1 and the experimental data for ethane molecules in zeolite AlPO 4 -5 are compared in Figure 8.
The dependence of the diffusion coefficient on the filling factor was investigated for several gases (CH 4 , C 2 H 6 , and C 3 H 8 ) in silicalite-1 membrane at times  > 10 s [57].The diffusion coefficients for these gases were found to increase more than 20 times as the filling factor varied from 0 to 0.9.This observation can be explained by the collective effect of density perturbation propagation through a 1D clustered system.The calculated dependence of the relative diffusion coefficient / 0 = (Im ()/ 0  2 ) =0 on the     55), ( 70), (87), and (85); points represent the experimental data for the diffusion of ethane in zeolite AlPO 4 -5 [61].
The insert shows the dependence of the relative diffusion coefficient on the filling factor obtained for several gases (solid line) using ( 84) and ( 85) and / = 1.Points are the experimental data from [57].
filling factor  in comparison with experimental results is presented in Figure 8 (insert).The figure shows that the diffusion coefficient is only weakly dependent on the filling factor for  < 0.5 (Figure 8 insert) but it rapidly increases when the filling factor becomes larger.This behavior of () can be easily understood as at low filling factors, the particles inside the membrane channel do not interact with each other, no clusters are formed, and particle motion is purely diffusive.
The increase of the diffusion coefficient for  > 0.5 is due to the cluster formation inside the channel and activation of the CD mode (Section 2.1).

Formation of High-and Low-Density Clusters in a 1D
System.In recent experiments [74], the one-dimensional (1D) Au chains formed on the (110) NiAl surface have been investigated at low temperature ( = 12K).This surface exhibits 1D grooves made of the rows of Ni atoms, which suppress the transverse movements of Au adatoms allowing them to move only along the groove.The measured Au-Au distances in such chains were about 2.9 Å that is very close to those in bulk gold [74].The properties of 1D Au chains were also investigated in experiments [8,12,[75][76][77] where they were either pulled out of a gold surface by an STM tip or created in a process of thin film etching.Such chains were stable during the observation time and consisted of several gold atoms.The Au-Au distances in some of them appeared to be strikingly large, up to 3.5-4.0Å [75][76][77], whereas in others much smaller Au-Au separations (∼2.5-2.9Å) were observed [8,12,[78][79][80].These results suggest that a 1D gold chain can exist in two different states characterized by different interatomic distances.The ab initio calculations performed at zero temperature have shown that a 1D gold wire should break when the Au-Au separations exceed ∼2.9 Å [81][82][83][84].Being in agreement with the experiments reporting the Au-Au distances of 2.5-2.9Å, this result, however, does not explain the existence of large Au-Au separations.It has lately been demonstrated that the presence of hydrogen or other light atoms in the wire structure can stabilize gold chains and significantly alter their properties [3,4,[82][83][84][85][86][87][88].Assuming a somewhat ordered distribution of impurities within the wire, one can easily explain the experimentally observed large Au-Au separations [6,82,83,87].Although this finding offers an intriguing possibility to modify the wire properties by an intentional incorporation of impurities (Im), a spontaneous formation of locally ordered Au-Im-Au structures seems unlikely and remains to be proven experimentally.Another explanation has been proposed by Sun et al. [89] who suggested that the large Au-Au distances might result from chain melting.It is known, however, that phase transitions do not occur in 1D systems [34]; thus, the proposed model should be taken with precaution.
Monoatomic chains can disintegrate due to atomic movements in both longitudinal and transverse directions caused by thermal fluctuations.The probability of disintegration is determined by the characteristic frequencies of atomic vibrations in the chain, and the energy barrier atoms need to be overcome in order to leave the chain, which is comparable with the energy per atom in the chain.A change of interatomic distances due to transversal or longitudinal motion of one of the atoms leads to a transition to a nonequilibrium state, and the energy per atom increases.As the energy barriers in the longitudinal and transversal directions are comparable [90], the disintegration probabilities for these two mechanisms appear to be comparable as well.Therefore, further on we consider only one mechanism occurring due to longitudinal fluctuations.
The probability of disintegration of a chain due to longitudinal fluctuations is determined by the frequency spectrum of relaxation of density fluctuations (relaxation spectrum) of the quasi-1D system.To study the relaxation spectrum of a quasi-1D system, we apply the previously stated response functions for density functional method, which allows one to describe the equilibrium and dynamical properties of a system if its interaction potential is known.Within this method, one can find the response function (, ) that determines the relation between the density fluctuation (, ) and the external field as well as all the space-time correlations in the system.Additionally, one can find the free energy of the quasi-1D system taking into account density fluctuations ( [38] and Section 2.1).A density fluctuation appearing in the quasi-1D system can be determined as Here  is the modulus of the wave number, and () is the amplitude of the density fluctuation, which in the analysis of cluster relaxation can be arbitrary and adequately considered to be a constant for  in the interval { min = 2/;  max = 2/}.Clusters forming in the system are characterized by the parameter (, ), which is defined as a Fourier integral (87) in the interval { min ;  max }, where  is the cluster length, and  is the atomic diameter.The cluster lifetime () is determined by the maximal magnitude of the imaginary part of the spectrum () ( −1 = Im ()) in the interval  min <  <  max .The relaxation spectrum of the 1D system () is determined by (70).If the relaxation of such a system is determined by diffusion with the coefficient  0 , then  0 () = − 0  2 . =  0  is a dimensionless density. 0 is the average number of particles per length unit.](, ) is a Fourier image of the pair correlation function for interacting particles, which is related to the Fourier component of the pair distribution (()) as ]() = () − (), where () is the Dirac delta-function.Note that the divergence must be removed from () before any analysis.Otherwise reverse Fourier transformation and the transition from wave numbers to coordinates are impossible due to the integral divergence.For a quasi-1D system, the pair distribution can be calculated exactly for an arbitrary short-range interaction potential Φ() [44]: Here  1D and  are 1D pressure and temperature, respectively.Equation ( 88) allows one to evaluate the pair distribution and consequently the relaxation spectrum of the 1D system.We have performed an extensive analysis, which has shown that the properties of quasi-1D clusters are determined by the presence of the repulsive and attractive terms in the potential but not by the specific form of Φ().In particular, we have obtained similar results using rectangular potentials, and the Lennard-Jones as well as Morse potentials.We notice that this finding correlates well with the known results of the liquid state theory [45,46,91].For the sake of simplicity, we present here the calculations of 1D pressure  1D () ( = (/) is dimensionless density) and pair distribution done for the square-well potential with the well width  and depth  0 : The performed calculations give the  1D () dependence: Equations ( 88) and ( 90) allow one to analyze the possible states of the quasi-1D system of the interacting particles of an arbitrary density.They also let us to estimate the relaxation time of a density fluctuation (cluster) as well as the interatomic distances in the clusters.We choose to perform this analysis using the parameters of the pair potential (89),  0 = 1.1 eV,  = 2.5 Å, and  = 2.5 Å, corresponding to the gold atoms [80].For these parameters, we calculate the compressibility  = −(/ 1D ) (Figure 9), pair distribution (Figure 10), and dimensionless spectrum Ω ≡ ()/ 0 () (Figure 11).At low densities ( < 0.67), the compressibility increases with increasing density and it approaches its first maximum at  =   = 0.18 (Figure 9).The compressibility has a second maximum at  =  (1)    = 0.71.The derivative of the compressibility equals 0 at extrema.Therefore, the presence of extrema can indicate changes in the system state.As we show in the following, the presence of maxima does not depend on the specific form of the pair interaction potential Φ(), but it is determined by the  0 / ratio ( 0 is a characteristic depth of an arbitrary potential).The second maximum of the compressibility disappears if  0 / ≪ 1.When  0 / ≫ 1, the second maximum always exists if the condition  ≥  is satisfied regardless of the specific form of the potential Φ().Also we have to note that the ratio  0 / provides the existence of the second maximum whereas the position of this maximum does not depend on the magnitude of this ratio.This allows us to conclude that in the limiting case  0 / ≫ 1, the variation of the temperature of a system does not lead to any qualitative changes in a state of a system.
To demonstrate how the physical state of the system changes at the compressibility maxima, let us analyze the pair distribution () that determines the probability to find atoms at a certain distance from each other.It has been shown in Section 2.1 that the function () → 1 when  → ∞.The distances between the maxima of () are equal to 1, 2, 3, and so forth, atomic diameters.Hence, the separation between these maxima corresponds to the interatomic distance in the clusters [45].
In the vicinity of the first compressibility maximum (  = 0.18), the distance between the peaks of () and, correspondingly, the interatomic distances in the clusters are ∼3.6 Å (Figure 10(a)).At densities close to the second maximum, interatomic distance in the cluster ≤3 Å (Figure 10(b)).As it was mentioned above in the limiting case  0 / ≫ 1, the variation of the temperature of a system does not lead to any changes in the state of this system.Therefore, at  = 300 and 12 K, the magnitude of the gold cluster interatomic distance will be the same.
As it has been noticed previously the qualitative shape and the distance between the peaks of the pair distribution function are independent of the concrete form of the interaction potential.The pair distributions () calculated for  = 0.15 with square-well potential (solid line) and Morse potential for the Au atoms (dashed line) are compared in Figure 11.
Figure 11 demonstrates that the distances between the peaks are the same for both potentials, but the peaks for the Morse potential are higher compared to the ones obtained with the square-well potential.This difference is due to the different characteristic length scales of the Morse and squarewell potentials.The characteristic length scale of the Morse potential is larger than that of the square-well potential that leads to a larger atomic attraction at long distances.Note that in the case of the Morse potential  0 in ( 89) is a potential well depth and  +  is the characteristic length scale of the potential.
Further, we analyze the stability of the 1D clusters.The change of the free energy of the 1D system (ΔΦ) due to the cluster disintegration determines the probability of such a disintegration [34].ΔΦ can be divided into the following contributions: ΔΦ 0 that is the change of the free energy due to the pair interactions in the cluster; ΔΦ 1 that is the change of the free energy coming from changes in the cluster electronic states; and the fluctuation part ΔΦ  that is an average fluctuation of the free energy of the 1D system caused by cluster disintegration calculated for the stationary state of the 1D system.Therefore, the probability of cluster disintegration can be written as To estimate the cluster lifetime  ∼  −1 , we analyze the contributions to ΔΦ.Since there are no phase transitions in 1D systems, ΔΦ  is always positive and its magnitude is determined by the spectrum Ω() [46] (Figure 12): It follows from (92) that ΔΦ  , which is related to the probability of cluster disintegration due to the loss of one or several atoms, is determined by the maximal magnitude of Ω() in the interval ( min = 2/;  max = 2/).In the vicinity of the first maximum of the compressibility at  = 0.2 (Figure 9), ΔΦ  ∼ Ω(  ).Here   is the wave number corresponding to the maximum of Ω(),    ∼ 2 (Figure 12(a)), and  is the number of particles in the cluster.It follows from (91) that in this case the magnitude of   is large for all the clusters with the length exceeding several atomic radii.In the vicinity of the second maximum of the compressibility at  = 0.7 (Figure 9), ΔΦ  becomes larger as the cluster length increases (Figure 12(b)).In this case for all clusters with the length exceeding   , for which Ω(2/  ) ≫ 1, ΔΦ  ∼ Ω(2/  ) and the probability   is large.Thus at all densities, only clusters of a certain maximal length can exhibit a noticeable lifetime, and the formation of long clusters appears to be unfavorable.ΔΦ  is the sum of the changes of the free energy of the pair interaction between atoms in the cluster  0 and the free energy of the cluster electronic states  1 (ΔΦ  / =  0 +  1 ).In the limiting case  ≪ | 1 |,  1 , is negative [92] and it is determined by the overlap of the atomic orbitals.This overlap is determined by the penetrability coefficient Γ of the potential barrier separating the electronic potential wells of the neighboring atoms (Section 2.4).In the quasi-classical approximation, Γ depends exponentially on the width and height of the potential barrier (Section 2.4).It has been shown [92] that due to collective effects, the penetrability coefficient Γ depends not only on the interatomic distances, but also on the difference between the maximum and minimum values (Δ) of the pair distribution () (Figure 10).The interatomic distances in the clusters formed at the densities about  ∼ 0.2 are larger than those for the clusters existing at  ∼ 0.7 (Figure 10), but Δ( = 0.75) > Δ( = 0.2).Using the approach of [92], we find that in this case the magnitudes for the clusters of both kinds appear to be similar.Thus, for our analysis of cluster stability, we can use the value | 1 | ∼ 0.5 eV obtained from the measurements of the cluster conductivity in the scanning tunneling spectroscopy (STS) experiments for the Au chains on the NiAl surface [74].The energy  0 per atom can be written as  0 ∼ (/2) ∫ ∞ 0 Φ()() [46].Calculations show that  0 is negative and | 0 | ∼ 1 eV and | 0 | ∼ 0.7 eV for the clusters at  ∼ 0.7 and  ∼ 0.2, respectively.
Let us now compare the lifetime and interatomic distances calculated for the 1D clusters to the known properties of the 1D chains fabricated in experiments [74][75][76][77][78][79][80].The results of the STS experiments [74] as well as the ab initio study of the electronic structure of gold wires on (110) NiAl [93] suggest that the main role of the substrate is to be a natural template insuring the one-dimensionality of the Au chains.This allows one to suggest that the observed Au-Au distances (∼2.9 Å) in such chains to a large extent might be determined by the one-dimensionality of the system.Taking into account the minimal possible distance between Au atoms (∼2.5 Å [81]), we can estimate the density in these chains:  = (2.5/2.9)= 0.87.Thus, the fabrication of Au chains on (110) NiAl can be understood as an artificial formation of 1D clusters with the density  fulfilling the condition  >  (1)   ( ( = 0.74) to 2.7 Å ( = 0.9) or 2.5 Å ( = 1).The lifetime of such clusters of the length  can be estimated as the characteristic time ( ∼   exp(/)) during which an atom stays in the potential well of the depth  = | 0 | + | 1 | − Ω(2/) (see (91)).Here   ∼  −1 0 ( = 2/) ∼ 10 −12 s.The calculations of the spectrum Ω() at  = 0.87 show that, for example, for the cluster of the length  ∼ 15 Å, containing ∼5 Au atoms Ω(2/) ∼ 10.In this case using the estimates obtained previously (| 1 ∼ 0.5 eV, | 1 | ∼ 1 eV), one can estimate the lifetime of such chains to be  > 1000 s ( = 300 K).Notice that in the absence of collective effects and attractive interaction between atoms in a quasi-1D system the lifetime of the chain of  ∼ 15 Å would be very short,  ∼  −1 ( = 2/) ∼  2 / ∼ 10 −9 s (here  ∼ 10 −5 s m 2 /s is the diffusion coefficient for Au atoms at the surface of NiAl [94] at  = 300 K).
The formation of suspended Au chains by pulling them out of the Au surface or by etching an Au film with an electron beam [2,8,12,[75][76][77][78][79] corresponds to conditions where considerable tensile stress is applied to the chain.The presence of an external force corresponds to the appearance of the one-dimensional pressure  1D = −(/), where  is a free energy of the 1D chain.One can see from (90) that 1D chains of finite length are affected by two types of forces: forces corresponding to the effective repulsion of the atoms from each other as a result of their movement (the first term in (90)) and forces corresponding to their mutual attraction that is determined by interaction potential of atoms with each other.The attractive force | V 1D | depends on depth of potential well and radius of the potential  and equals  V 1D = −( 2 / 2 ) in the van der Waals approximation.Thus a cluster is in a compressed state.The value of compression is The presence of such a compression can lead to loss of stability of the cluster in the transverse direction and formation of a zigzag configuration of a chain.Therefore, to form a stable linear chain, it is necessary to attach an external force to the end of a chain to compensate the compression | V 1D |.The presence of an external force corresponds to the appearance of one-dimensional pressure | V 1D | = −(/), where  is the free energy of a onedimensional chain.Let us estimate the lifetime of these chains.
Within the framework of these assumptions, we can also estimate minimal lifetime of chains of other metals, such as Ag, Ni, and Pt.Due to the lack of experimental data on the characteristic values of  0 for these materials let us assume where  is the binding energy per atom for dimer [95].The calculated values of minimal lifetime of one-dimensional chains are presented in Table 2.As one can see, the material of a one-dimensional chain greatly affects its stability.
Let us estimate the minimal length of stable chains and the force that is necessary for their formation.Note that if we consider the chain as a rod of length  with a radius equal to the atomic radius , then such a rod loses its stability when it is compressed with a force greater than   =  3  4 / 2 [96], where  is the Youngs modulus.So, if the compression pressure of the chain is | V 1D | >   , then the chain is unstable.Therefore, the minimum compressive force that ensures stability of the chain is | V 1D | =   .Chains of metal atoms have a conductivity comparable to the metal one [13], so the range of the interatomic potential is comparable to the length of the chain  = , where  ∼ 1 and varies from metal to metal.Relations | V 1D | = ( 2 / 2 ) and   = ( 3  4 / 2 ) allow to estimate the minimum number of atoms in a cluster  min = ( 4 /) 1/3 , where  is an interatomic distance in the cluster.The value | V 1D | is the force that compresses a stable cluster without external influences on it.In this case, | V 1D | = ( min /) is a force  0 required to break the cluster with the minimum number of atoms in it.Table 2 shows the results of the estimations of the minimum number of particles in the chain  min and the corresponding force required to break the chains  0th for chains of different material.It was assumed that  Ag = 1 for the chains of Ag.For the rest of metal chains, this coefficient was estimated by using the relative conductivities   =  Ag (  / Ag ),  = Au, Ni, Pt.

Transport of a Two-Component Mixture in One-Dimensional Channels.
The dependence of transport selectivity on the pressure in a two-component system of atomic particles is nonmonotonic: the selectivity attains its maximal value and then decreases instead of increasing.Thus, the mechanism of diffusion enhancement in 1D channel, which was proposed in Section 2.1, is inapplicable for two-component mixtures.
The transport of a two-component mixture in subnanometer channels was considered earlier [51,97] on the basis of the generalized phenomenological Stephan-Maxwell equation.The authors of these publications used the dependence of the chemical potential on the fill factor, taking into account the finite size of particles, but disregarding their interaction.The dependences obtained in [51,97] also indicate a monotonous increase in the diffusion coefficients, fluxes, and selectivity upon an increase in  and, hence, fail to describe the experimental data for two-component mixtures.

Ground State of a System in 1D Channels.
Let us consider the surface of a porous body in contact with a twocomponent gas mixture at temperature  and pressure .Suppose that  adsorption centers are located on the surface.We assume that the particles on the outer surface do not interact with one another.We also assume that the energy of a gas molecule on the surface is equal to  ()  0 ,  = 1, 2, depending on the species of the molecule.We also suppose that  cylindrical channels ( ≫ 1) of diameter  and length  emerge at the surface.We assume that the diameter of a cylinder is comparable to the maximal diameter of the gas molecule.Let us assume that  ()  1 is the binding energy of the th particle at the mouth of a channel,  ()  1 is the number of particles of the th species above the membrane,  () is the total number of particles of the th species in the channel,  () is the total number of th particles in the channel and on the surface, and  0 is the number of "seats" in the channel.Then the partition function of the grand canonical ensemble taking into account the interaction of system of atomic particles in the channel has the form 1 ( (1) ) −  (1) ] ( (1) −  (1) )! [ − ( (1) −  (1) )]! ×  0 !exp ( () 2  () )  () !( 0 −  () )! exp ( () 1 )  int ( (1) ,  (2) ) , Here, the second sum corresponds to summation over configurations  =  (1) ,  () 2 is the binding energy of particles of the th component in the channel,  is the chemical potential, and  int ( (1) ,  (2) ) is the partition function corresponding to the inclusion of particle interaction in the channel.Considering that the main contribution to the partition function comes from states with a large number of particles ( () ≫ 1), we can replace partition function  int ( (1) ,  (2) ) by the partition function for particles interacting in the channel, calculated for the mean value  () of the number of particles of each species in the channel.From the physics standpoint, this corresponds to "averaging over channels, " when a single channel (whose state is calculated by averaging the parameters of particles in  channels) is considered instead of the large number of channels.It should be noted that this approximation is possible since the number of channels having a diameter of  < 1 nm and emerging on a square centimeter of the zeolite membrane surface is large:  ∼  −2  ∼ 10 14 cm −2 ≫ 1, where   is the characteristic distance between the channels on the membrane surface (Section 2.1).Using relation (93) and following the method proposed in Section 2.1, we obtain for the adsorption isotherm   av  =    +  (1)   0 where  = / 0 is the mean distance between the seats in a channel,  is the concentration of the first component in the gas phase above the membrane,  1 and  2 are the concentrations of the components in a channel,  1 and  2 are the diameters of molecules of the first and second species,  av ≡ ( 1 +  2 )/2 is the average diameter,  is the fill factor of the channel, and  int has the meaning of the free energy of interaction per particle in the system of atomic particles in the channel.The concentrations of the components in the channel are chosen so that the following relation is satisfied: In accordance with relations (94), the calculation of the adsorption isotherm is reduced to the calculation of the value of  int : where () is the potential energy of interaction of particles of species  separated by distance  from each other.It is convenient for subsequent computations to introduce, instead of coordinates   of individual particles in a channel, "pair" coordinates    , where  = 1,  = 1 corresponds to the location of two particles of the first species in the vicinity of point   ;  = 2,  = 2 corresponds to the same for particles of the second species, and so on.Obviously, the coordinates are connected through the relation Here   has the meaning of the total number of pairs of closely spaced particles of species  and  and  1 ( 2 ,  3 ) is the effective length "occupied" by all pairs  11 ( 12 ,  22 ) of particles.Considering that  1 +  2 +  3 =  and, as a consequence, the impossibility of mutual penetration of particles through one another, we can reduce integration with respect to coordinates  to integration with respect to coordinates   , which corresponds to summation over all possible configurations of pairs of particles of the first and second species in the channel.Consequently, we obtain the following expression for   from (97): Using the delta function representation in the form of a line integral [45] and carrying out integration with respect to    independently for different values of indices  and , we obtain Integrating this relation with respect to   for the simplest form of the intermolecular interaction potential of the hardsphere type we obtain Here we have used the definition   =     ( is the total number of particles in a channel) and the representation of a factorial in terms of the gamma function, ! = Γ( + 1).
Eliminating the factor corresponding to an ideal system of atomic particles from relation (102) and substituting the result into (96), we obtain the following expression for  int : As  1 → 0 or  2 → 0, expression ( 103) is transformed into the corresponding expression for the free energy of a one-component system of atomic particles in a 1D channel (Section 2.1).From relations ( 103) and ( 94) for the adsorption isotherm of a two-component system of atomic particles, we finally derive the expressions Using these relations and specifying the external conditions, we can obtain the fill factor and the component concentration in the channel.For  = 0 and  = 1, isotherm (104) transforms into the isotherm for a one-component system of atomic particles, which was obtained in Section 2.1.Let us now analyze the ground state of a mixture in 1D channels.It is known (Section 2.1) that short-lived clusters can be formed in 1D channels in the case of a one-component system of atomic particles.The formation of clusters can apparently be expected in the case of a two-component system of atomic particles also.In order to describe possible transition to an inhomogeneous clustered state, we must calculate and analyze the free energy of the system taking into account fluctuations and the equation determining the amplitude of these fluctuations in the system (equation for the order parameter), which follows from the free energy.Since molecular transport depends on the lifetime of clusters, the analysis of the equation for the order parameter aims at determining the characteristic lifetime of clusters as a function of the channel fill factor and, hence, as a fucntion of pressure, temperature, and composition of the mixture above the membrane.In order to describe the ground state of the 1D system and relaxation process to the ground state, one is to know the free energy  of the 1D system and two-component functional Δ[  ],  = 1, 2, determined by the relations (39), (44).The functional Δ[  ] can be written in the way of ( 49): Here Δ 0 [  ] is the functional of the non-interacting system and Δ int [  ] corresponds to the particles interaction: where  0 is the free energy per temperature unit of the ideal (non-interacting) system,  1 is the fluctuation part of the free energy due to interatomic interaction in the mixture.Quantity Δ can be expressed in terms of the response function  , (,   , ,   ) of the two-component system, which is defined as where   (, ) is the density fluctuation of particles of species  and  ext (  ,   ) is a weak external field depending on time.Here and in the following, we assume summation over recurrent indices.Response function   (, ) can be derived by varying directly the mean value of functional Δ over (, ) [38]: Using this relation, we can obtain the expansion of functional Δ into a series in deviations  of density from its mean value in the case of a multicomponent system.In the Fourier representation, the result has the form Here   (, ) is the response function of the multicomponent system in the Fourier representation, and the asterisk indicates complex conjugation.It follows from the fluctuation-dissipative theorem [98] that response function   (, ) is directly proportional to the pair distribution; this enables us to calculate   (, ) for a 1D system with an arbitrary density, for which the pair distribution is known exactly [45].To calculate Δ, we use the local approximation that makes it possible to replace functional derivatives by ordinary derivatives [38]: Using this relation, we obtain from (109) ) . ( Here   (0)(, ) and   (, ) are the partial response functions for noninteracting and interacting particles, respectively.The quantity   (, ) describes the effective dynamic interaction between particles and is defined as Premultiplying expression (112) by ( β(0) (, )) −1 and postmultiplying it by β−1 (, ), we obtain Over long time periods ( ≫  mom ), function   (, ) can be represented in the form In view of homogeneity and the absence of correlations in a system of noninteracting particles, we have β(0) (, 0) = β(0) (0, 0) =   , where with  being the entropy of the two-component system of atomic particles divided by temperature in the absence of the intermolecular interaction.In this case, we obtain the following expression for the effective interaction from (115): Substituting this expression into (114), we obtain the following relation for response function β(, ): In order to calculate the response function β(, 0) for a system of interacting particles, we can use the fluctuationdissipative theorem connecting response function β(, 0) with the pair correlation function for a system of interacting particles: Here ]  () is the pair correlation function divided by temperature [46].Relation ( 119) is a generalization of the relation (68) for one-component systems to two-component systems.
Thus, the calculation of the response function can be reduced to the calculation of the pair correlation function or function   called a pair distribution [46] and connected with the pair correlation function through the relation ]  () =   () −      () . (120) We will calculate the pair distribution using the method described in [45] and determine  11 ().Functions  12 () and  22 () are calculated similarly.Pair distribution  11 () can be written in the form [45] where the function Ψ where  ,, () is the configuration integral corresponding to the presence of  pairs of particles of the first species,  pairs of particles of the first and second species, and  pairs of particles of the second species in a channel.For  =  11 ,  =  12 , and  =  22 , the expression for configuration integral  ,, () coincides with (98) [45]: Here  has the meaning of a coordinate and   () is the Laplace mapping of function   () ≡ exp(−  ()/).In the case of an intermolecular interaction potential in the form of that in the hard-sphere model, we have (124) Substituting these relations into (123), we obtain the following expression for the configuration integral: Here  has the meaning of one-dimensional "pressure." Substituting relation (125) into (122), we obtain the following expression for pair distribution  11 (): In the long-wave approximation (0 ≤  ≤ 2/ max{ 1 ,  2 }), expression (126) contains no poles and we obtain Functions  12 () and  22 () can be calculated similarly.Substituting expressions (127) and ( 119) into (118), we obtain the final expression for response function β(, 0): ) . ( Here is the free energy of a two-component system of atomic particles with interaction, referred to temperature, which can be calculated by differentiating configuration integral   twice with respect to concentration.Function β(, ) can be obtained by generalizing the response function of the onecomponent system of atomic particles, derived in Section 2.1, to the case of a two-component mixture: Here  1 ,  2 ,  are the diffusion coefficients and the relaxation frequency spectra for pure components.Substituting relations (128) and ( 131) into (118), we obtain the following expression for response function (, ): It is convenient for subsequent analysis to pass to new variables.For a two-component system, the relation  1 + 2 = 1 holds and quantity  can be represented in the form where  can be interpreted as an order parameter.Substituting (133) into (111), we obtain the expression for functional Δ: where Series (134) can be summed and written in the form The formal proof of the correctness of representing quantity Δ in form (136) without analyzing the convergence of the corresponding series follows from the expansion of expression (136) into a power series of order parameter .This expansion exactly reproduces series (134) (at least, to within terms on the order of  5 ).Using relations (132), (135), and (136), we can obtain the dependence of functional Δ on the wave vector and order parameter of the system.The minimum of functional Δ determines the ground state of the system.Substituting relation (132) into (136), evaluating the derivative of Δ with respect to the order parameter, and applying the inverse Fourier transformation in frequency  to the resulting relation, we obtain where Equation ( 137) has the same form as the equation usually used for the order parameter.The role of the order parameter is played in this case by the quantity defined by relation (133).It should be noted that the order parameter here could be either positive or negative.The cases when  > 0 and  < 0 correspond to the density fluctuation of the first and the second components, respectively.The role of the diffusion coefficient is played by quantity D depending on the fill factor and the concentration of components in channels via coefficients   and   defined by relations (116) and (128).It should be noted that the concentration of the components in a channel (and its fill factor) could be determined unambiguously from isotherm (104) proceeding from the pressure, temperature, and composition of the system of atomic particles.Analysis of functional  1 as a function of order parameter  and wave vector  makes it possible to investigate the possibility for a transition of the system to an inhomogeneous state.For example, when the minimum of functional  1 is attained for  ̸ = 0,  = 0, a conventional phase transition to a homogeneous state takes place [99].The minimum at  = 0,  ≡   ̸ = 0 corresponds to the propagation of a density wave over a distance of  ∼  −1  .In the case when the minimum is attained at  ̸ = 0 and  ̸ = 0, a transition to an inhomogeneous state with clusters formed in the system is realized [99].Local minima here determine metastable states.Figure 13 shows the graphs illustrating the dependence of  1 on order parameter   and wave vector  for a mixture in which one of the gases (first) is a strong sorbate, while the other gas is a weak sorbate.The chosen parameters ( = 3.8 Å,  1 = 4.3 Å,  2 = 3.6 Å,  1 = 0.61 eV,  2 = 0.38 eV,  = 2000 kPa, and  = 300 K) correspond to a methane-butane molecular mixture [51] and make it possible to describe adsorption and fluxes of these one-component gases in zeolite membranes.It should be noted that the choice of parameters , , and  determines unambiguously the filling factor  and concentrations  1 and  2 of gases in a channel.For low concentrations of the first component in a gas mixture and for fill factors  ≤ 0.7, functional  1 has only one minimum at  = 0,  = 0 (Figures 13(a) and 13(b)).This corresponds to a homogeneous state of the system.An increase in the concentration of the first component leads to the emergence of local minima in functional  1 .

Journal of Nanomaterials
In this case, the global minimum is shifted to point  =   > 0,  =   .In this case, the ground state of the system becomes spatially inhomogeneous, which corresponds to the formation of clusters in the system.This situation will be analyzed in greater detail later.Equation (137) for the order parameter also allows us to calculate the characteristic relaxation times for the emerging density fluctuation depending on external conditions.For this purpose, we expand functional  1 (  , ) into a series in the vicinity of   = 0 to within the first nonvanishing term: Substituting (139) into (137), we obtain The quantity  is the lifetime of the th mode of the density fluctuation emerging in the vicinity of   = 0.The curves describing the lifetime of the emerging fluctuation as a function of the wave vector for different mixture compositions are shown in Figure 14.It can be seen from Figure 14 that, for a low concentration of the highly adsorbed component, the lifetime decreases monotonically with increasing wave vector, while the free energy has a minimum at  = 0.
From the physics standpoint, this means that the state of the system in a channel for a given composition of the mixture is homogeneous and the growth of clusters is energetically disadvantageous.It should be noted that the fill factor of the channel is small in this case, but short-lived clusters with the lifetime can be formed due to density fluctuations.Here  =  0 exp(−  /) is the diffusion coefficient of solitary particles,   is the diffusion activation energy, and  is the characteristic size of fluctuations.Characteristic size  is determined, as in Section 2.1, from the position of the extremum of the imaginary part of the spectrum divided by the square of the wave vector.It should be noted that the value of  increases with the fill factor of the channel Section 2.1.Expression (142) can be derived by expanding order parameter (137) into a series in the vicinity of  = 0.
For low concentrations, expression (141) coincides with the characteristic diffusion decay time for density fluctuations.An increase in the concentration of the first component of the mixture to  = 0.3 elevates the fill factor of the channel to  = 0.74 (see Figure 13(c)), while  1 acquires two more minima for  ∼ 1. Figure 15 shows the dependences of the free energy on the order parameter for different values of the wave vector at  = 0.3.It can be seen from the figure that the free energy minimum is attained at point  = 0 for values of the wave vector 0 <  < 0.4.For  ≈ 0.5, functional  has two minima (at  (0)  = 0 and  (1)  = 0.35), and the state with  (1)   is separated from the state with  (0)  by a potential barrier.Since the global minimum is attained at  =  (1)   , the state with  (0)  is metastable.The ground state of the system is attained at  =  (1)   and is clustered.It should be noted that, in accordance with formula (133), clusters of the first component are formed in the channel since  (1)   > 0. The lifetime of the clusters can be determined by passing from (137) for the order parameter to a stochastic differential equation (Langevin equation) and the corresponding Fokker-Plank equation via the introduction of additive noise.Steady-state solutions to the Fokker-Planck equation define the probability for the system being in the state described by order parameter  (1)   .For the lifetime, we have The decay of a cluster formed in the system occurs via overcoming an energy barrier of height .From the standpoint of physics, this can be interpreted as an effective increase in the diffusion activation energy for the particles in the cluster.However, for  = 0.3, the barrier height is small ( ∼ 0.2; see Figure 13(c)), and lifetime (142) is comparable to time (141).A further increase in the first component concentration in the mixture leads to negative values of the lifetime for a certain value of  =   < 1 (see Figure 14(c)).In this case, clusters with a size of  ∼  −1  are formed in the channel; for  <   , fluctuations decay over a finite time, while for  >   , the emerging density fluctuations evolve.Such a situation is typical of transitions to an inhomogeneous condensate state [99].
The curves describing the dependence of the free energy on the order parameter for different values of the wave vector are shown in Figure 16.Since functional  1 has a single minimum at  = 0 for  ∼ 0, the state of particles in the channel is homogeneous.However, as the wave vector increases, curve  1 () acquires two more minima separated from the state with  = 0 by potential barriers.A further increase in the wave vector leads to disappearance of the barriers between the state with  = 0 and the states corresponding to two other minima of  1 , and the state with  = 0 becomes unstable.
Thus, for a low concentration of the highly adsorbed component, the fill factor of the channel is small and the state of the system in the channel is homogeneous.As the concentration of the highly adsorbed component increases, the fill factor of the channel increases, and short-lived clusters of the highly adsorbed component can be formed in the channel.A further increase in the concentration (and the fill factor of the channel) leads to an increase in the lifetime of the clusters formed in the channel.At a certain concentration of the highly adsorbed component, the clusters in the channel become long lived and the ground state of the system becomes clustered.
The formation of clusters of one of the components in a channel may strongly affect the transport of the gas through the membrane (see Section 2.1).In the next chapter, partial gas fluxes through the membrane will be calculated, and the effect of clusters on the mechanisms of transport will be analyzed.

Transport in a Two-Component Gas in a High-Density
1D System.In order to calculate partial fluxes and to analyze the transport, we will use the approach proposed in Section 2.1.It will be shown in the following that the evaluation of the flux can be reduced to the calculation of the relaxation frequency spectra (k) for density fluctuations of the components.We can judge the mechanism of particle transport in a subnanometer channel from the type of the dependence of the spectrum on the wave vector of the system.We can write the expression determining the relaxation of the Fourier component (, ) of the number density of particles in a channel in the case of an arbitrary density (, ) under the conditions of slight deviation from equilibrium (see Section 2. where () is the relaxation frequency spectrum for the system under study; in the case of diffusion of noninteracting particles, this spectrum has the form [98]  =  0 () = − 2 .
(145) Equation ( 144) describes relaxation of the th component of density fluctuation for an arbitrary value of the wave vector.In particular, this equation for  ̸ = 0 makes it possible to describe the relaxation of density fluctuations and the propagation of a perturbation over a finite-size cluster in the case of its formation.For  → 0, (144) describes the relaxation of density fluctuation on a large spatial scale.This quantity is associated with macroscopic fluxes.In order to calculate the fluxes, we write the continuity equation For partial fluxes, we have where   is the density of the th component and   is the corresponding spectrum.The total gas flux is defined as the sum of partial fluxes of the components: It follows from (149) that partial fluxes are determined by relaxation frequency spectra   ().Thus, the problem is reduced to computing these spectra.The relaxation frequency spectra can be determined from the condition of the existence of nonzero density fluctuations for each component in an arbitrarily weak external field.Consequently, applying the Fourier transformation to (107), we find that the relaxation frequency spectrum can be determined by solving the following system of homogeneous equations: Here matrix  −1 (, ) of the response functions is defined by relation (132).In view of -function singularities for  → 0 that emerge in pair correlation functions (120) appearing in (132), we must calculate spectra   () and the values of   ( = 0), after which the renormalization procedure must be carried out.Using (132) and passing to the limit  → 0 for spectra   ( = 0), On the other hand, substituting (132) into (152), we obtain the following expressions for relaxation frequency spectra   () of the system under investigation: ] .We will carry out the renormalization procedure for spectra (154) taking into account relations (153), that is, imposing the requirement ω =  (0)   ( = 0) .
Taking into account the value of  (0)  , we find from relations (153) and ( 154 ] . (160) Using relations (104), (149), and (161) and passing to the limit for  → 0, we can obtain the dependence of partial fluxes of the components on the pressure, temperature, and composition of the mixture.For  ̸ = 0, expressions (144) and (160) enable us to analyze the relaxation mechanism of the emerging density fluctuation with a characteristic size of  ≈ 2/.
In order to analyze the experimentally observed partial fluxes, it is convenient to pass to the coordinate representation.Spectra (160) have the real and imaginary parts since pair correlation function ]  () is a complex quantity in accordance with relations (120) and (125).Separating the real and imaginary parts of the spectra and applying the inverse Fourier transformation to expression (149), we obtain Here   is the diffusion coefficient of a component,  is the channel length, and Ψ is the term emerging due to the mutual effect on the mixture components in a channel for large fill factors.It should be noted that, since  ≫   , where   is the characteristic size of clusters, passing to the limit for  → 2/ in relation (161) corresponds to averaging over characteristic scales of inhomogeneities in the case when the ground state of the system is clustered.Separating the terms linear in the wave vector from the real part of the spectrum, we obtain the following expression for Ψ: Here the first term describes the transport of molecules of the th component, which is induced by the effective intermolecular interaction, while the second term has the meaning of the "drag" effect familiar in the kinetics of mixtures.Substituting (162) into (161), we obtain the final expression for partial fluxes: Thus, it follows from this equation that the partial flux is the sum of three terms.The first term corresponds to diffusion transport, while the second term emerges due to field diffusion.These two terms in the  representation can be combined into one by introducing the effective diffusion coefficient It was shown in the previous section that, for large filling factors, clusters are formed in the system.The mechanism of density relaxation may be different for the case of transport over the characteristic scale , when 2 av / ≪ 1 and  → 0 (which corresponds to macroscopic transport of gas components through the channel), and for the case when  ̸ = 0, which corresponds to relaxation over distances comparable to the size of clusters formed in the system.To prove this, we must analyze dependences Δ]  ().In order to avoid cumbersome formulas, we consider this problem in greater detail in the limiting case of  1 = 1, which corresponds to a one-component mixture.
Passing to the limit for  1 → 1 in (160), we obtain the relaxation frequency spectrum for a one-component system; it is easy to obtain (70).It should be noted that, in accordance with the relation (70), the pair distribution and, hence, the pair correlation function associated with it depend not on the wave vector , but on the product : ] () ≡ ] () . ( Let us now write relation (164) in the form Separating the real and imaginary parts from relation (165), we obtain Here we have used the fact that This relation follows from the known relation [46] describing the increase in the compressibility with density.Using relation (165), we now expand Δ] to within first-order terms: The value of ]  (0) is finite for  ̸ = 0.For example, using direct expansion into the Taylor series, we can obtain Substituting relation ( 170) into (167), we get To calculate the fluxes for high densities ( ∼ 1), we must consider two limiting transitions in relation (172): for  → 1 and  → 0. From the standpoint of physics, the transition for  → 1 and  ̸ = 0 corresponds to the analysis of transport in a dense cluster.The transition for  → 0 and arbitrary  is equivalent to analysis of diffusion in a channel of length  ≫   , where   is the characteristic size of a cluster.For  → 1, we obtain, instead of relation (172), This relation shows that spectrum () corresponds to the hydrodynamic mode [46], for which quantity  0 /]  (0) is the effective velocity of sound.Passing in relation (172) to the limit for  → 0 and retaining the lowest order in wave number , we obtain This relation shows that the spectrum in this case is of the diffusion type with the diffusion coefficient which increases indefinitely as  → 1.
It should be noted that, in accordance with (172), both relaxation mechanisms (hydrodynamic and diffusion) operate in the system when  ̸ = 0 and  ̸ = 1.From the standpoint of physics, this corresponds to the diffusion transport between clusters (diffusion mode), over which density perturbations propagate (hydrodynamic mode).In spite of different interpretations of the sequence of limiting transitions, the calculated dependences of fluxes on external conditions are found to be equivalent from the physics standpoint.This is due to the fact that the hydrodynamic component is manifested in the second case in the increase in the effective diffusion coefficient  ≈  0 /(1 − ) 2 for high values of fill factor .
In the two-component case, the transition for  → 0 corresponds to the transition from relation (160) to (153).
It can be seen from relation (153) that the obtained spectra are of the diffusion type since  (0)  = −   2 ,  = 1,2, while quantities   and   are independent of .The limiting transition for  → 1 in the two-component case involves considerable computational difficulties.However, the above analysis shows that the relaxation frequency spectra for the components also acquire a hydrodynamic mode in view of the dependence of pair correlation function (120) on factor  for  → 1.
Depending on the sequence of limiting transition, a decisive role in relation ( 163) is played either by the term corresponding to diffusion or the nongradient (hydrodynamic) part of the flux.Since the gas flux for arbitrary fill factors is measured in experiments on the penetrability of the membrane over large distance ( ≫   ,  → 0), the limiting transition for  → 0 should be performed first.In this case, the effects associated with the formation of clusters are taken into account in the dependence of the diffusion coefficient on the concentration and fill factor of the channel.Numerical calculations show that the last term in relation ( 160) is insignificant in the entire range of concentrations and fill factors of the channel.Thus, the relaxation frequency spectrum of a two-component system can be treated as a diffusion spectrum in the entire range of concentrations and fill factors, where the role of diffusion coefficient is played by a quantity taking into account clustering of the components in the channel.
It should be noted that for calculating response function (132), relaxation frequency spectra (160), and, as a consequence, partial fluxes (163), we used the fluctuationdissipative theorem (119) presuming the homogeneous ground state of the system.In the case when the ground state of the system is clustered (see Figure 13 where (, ) is the probability of a particle being located at point .In accordance with the detailed balancing principle, in the case when states  and   differ from the equilibrium state insignificantly, we have where Δ is the energy difference between states  and  (178) Applying the Fourier transformation to this relation and using (151), we obtain the relaxation frequency spectrum in the form Thus, the application of relation (176) instead of (132) leads to the emergence of an additional exponential term in relation (163): where Δ is the depth of the potential well occupied by the system.In the case of small fill factors of the channel, the state of the system is homogeneous and Δ = 0 (see Figure 13(a)).
Then relation (180) transforms into (163).As the fill factor of the channel increases, the state of the system becomes clustered.For Δ ≪ , the exponential in relation (180) can be disregarded, and the flux virtually coincides with the flux calculated for a homogeneous state.This is due to the fact that expression (142) for the lifetime of clusters practically coincides with expression (141) for the lifetime of clusters in a homogeneous state.In this case, as noted previously, density excitation propagates via a cluster, leading to an increase in the effective diffusion coefficient.A further increase in the fill factor of the channel increases the value of Δ.The flux decreases thereby, which can be explained by an increase in the lifetime of the clusters formed in the channel.As a result, the mechanism of excitation transport via a cluster is not realized.Thus, for small fill factors of the channel, the transport in the system follows the diffusion mode.This can be demonstrated by passing to the limit of small fill factors  in relation (160).Then the second term in the denominator becomes insignificant and the spectra are reduced to diffusion spectra (153).An increase in the fill factor results in the formation of clusters.It was shown in the previous section that, depending on the fill factor of the channel, the clusters formed may be either short lived or stable in the case when the ground state of the system is clustered.For short-lived clusters (see Figure 13(b)), the transport between clusters occurs through diffusion, while the transport over clusters occurs through a rapid barrier-free transfer of density excitation, and the arrival of a particle from one side of a cluster leads to the emergence of a particle from the other side [100].In this case, the spectra contain both the diffusion and the hydrodynamic modes, and the effective diffusion coefficient increases in accordance with (175).An increase in the fill factor of the channel reduces the distance between clusters and increases their size.Since the transport over a cluster is faster than the diffusion transport, the effective diffusion coefficient and the partial flux increase.However, an increase in the fill factor also increases the lifetime of the clusters formed in the channel; as a result, the rate of excitation transfer over a cluster decreases.In the case when a cluster is stable (see Figure 13(d)), the transport of excitation over the cluster is ruled out, and the presence of such a cluster in the channel leads to blocking of transport in the system.
Using relations (160) and (180), we can calculate the partial gas fluxes proceeding from the data on the penetrability for pure components.Figures 17(a) and 17(b) show the dependence of partial fluxes on the mixture composition in the case when both gases are slightly adsorbed.The energies of interaction and the diameters of particles correspond to a methane-argon gas mixture ( = 3.8 Å,  1 = 3.6 Å,  2 = 3.0 Å,  1 = 0.38 eV,  2 = 0.25 eV [51],  = 100 kPa, and  = 300 K).
It can be seen from Figure 17(a) that an increase in the concentration of the first component for fixed pressure and temperature reduces the filling factor with the second component, while the filling factor with the first component increases thereby.A decrease in the total fill factor of the channel in the concentration range 0 <  < 0.4 of first component takes place due to depletion of the channel in the second component (Figure 17(b)).An increase in the total fill factor of the channel for  > 0.4 is associated with preferred enrichment of the channel in the first component (see Figure 17(b)).The partial flux of the first component in this case increases monotonically due to an increase in the filling factor with the first component (see Figure 17(a)), while the partial flux of the second component decreases due to a decrease in the filling factor with the second component.The transport is of the diffusion type, and clusters are not formed in the channel in view of small fill factors of the channel ( ≤ 0.25 for any composition of the mixture).
Figures 17(c) and 17(d) show the dependence of partial fluxes on the mixture composition in the case when one of the gases (the first) is a strong sorbate, while the other gas is a weak sorbate.The energy of interaction and the diameter of particles correspond to a butane-methane gas mixture ( = 3.8 Å,  1 = 4.3 Å,  2 = 3.6 Å,  1 = 0.61 eV,  2 = 0.38 eV [51],  = 2000 kPa, and  = 300 K).The behavior of the total and partial fill factors of the channel upon an increase in the concentration of the first (highly adsorbed) component in this case coincides qualitatively with the case of slightly adsorbed gases described previously: the partial filling factor with the first (highly adsorbed) component increases monotonically against the background of a monotonic decrease in the filling factor with the second (slightly adsorbed) component (see Figure 17(c)).
The behavior of partial fluxes in this case is less trivial.It can be seen from the figure that the partial flux of the first component increases in the concentration interval 0 <  < 0.2.This is due to the fact that, in accordance with relation (141), short-lived clusters (see Figure 13(a)) whose size increases with concentration are formed for such concentrations and fill factors of the channel.The transport over clusters occurs via a barrier-free transfer of density excitation, which increases the effective diffusion coefficient.This process is similar to the transport in one-component systems for high fill factors (see Section 2.1).It should be noted, however, that the formation of clusters and the mechanism of transfer of density excitation through a cluster in a twocomponent mixture can be realized for smaller fill factors ( = 0.71; see Figures 13(b), 17(c), and 17(d)) as compared to the one-component case ( ≥ 0.8).An increase in the concentration of the first component increases the lifetime ∼50 m and a random orientation of crystals was deposited on the substrate.Zeolites have a complex crystalline structure based on silicon oxide with admixtures of sodium and aluminum [51].The crystalline structure of zeolite ZSM-5 is formed by straight channels with an elliptical cross section (0.57 × 0.52 nm 2 ) intersecting sinusoidal channels with a circular cross section of diameter 0.54 nm (Section 2.1).The measured adsorption capacity and the calculated fraction of molecules in the channel intersections make it possible to treat the channels as one dimensional for the molecular mixtures under study (Section 2.1).The experiments [51,97] were carried out at temperatures from 300 to 700 K under pressures from 25 to 500 kPa on the external side of the membrane.The purity of the mixture components was higher than 99.95%.The mixture composition at the membrane exit was measured on a quadrupole mass spectrometer with a sensitivity of approximately 25 ppm when a Faraday cell was used for recording.The selectivity measured in the experiments was determined as the relative change in concentration  of -butane molecules and concentration 1 −  of methane molecules at the membrane entrance () and exit (): It should be noted that, for a constant composition of the mixture at the membrane entrance, the selectivity is proportional to the ratio of partial fluxes of the mixture components at the membrane exit: It was found experimentally [51] that the partial flux of the highly adsorbed gas -C 4 H 10 at  = 300 K under a total pressure of  = 100 kPa changes considerably (by a factor of several units) in the presence of the slightly adsorbed component CH 4 as compared to the flux of a pure gas, while the partial flux of the slightly adsorbed gas changes by two orders of magnitude in the presence of the strongly absorbed component.Depending on the mixture composition, selectivity  varies and attains its maximum value  max = 380 for a concentration ratio of 5 : 95 at the membrane entrance.It was shown [97] that the dependences of partial fluxes on the mixture composition and external pressure for the mixtures studied (C 2 H 6 -CH 4 and C 3 H 8 -CH 4 ) are monotonic in the range of experimental pressures.The partial flux of methane decreases in the presence of the second component (ethane or propane), while the partial flux of the second component increases with its concentration in the mixture.The selectivity of both mixtures increases monotonically with the concentration of the more highly adsorbed component in the mixture.The dependences of the partial fluxes on the total pressure of the mixture are also monotonic, but the pressure dependence of selectivity in ethane for the mixture C 2 H 6 -CH 4 has a peak at a pressure of  ≈ 300 kPa [97].Similar dependences of partial fluxes were observed in experiments with another mixture of slightly adsorbed gases (CO 2 -N 2 ) [53] (Figure 18).The theory developed in the previous sections enables us to describe peculiarities in the transport of binary gas mixtures in zeolite membranes.In the case of a mixture of slightly adsorbed gases, when the fill factor of a channel is small ( ≪ 1) for any mixture composition, the system relaxes in accordance with the diffusion mechanism corresponding to diffusion of solitary particles without the formation of clusters in a channel.In this case, in accordance with (104), (160), and (180), the partial flux of the first component decreases monotonically upon an increase in the concentration of the second component, while the partial flux of the second component increases with its concentration in the mixture.
Figure 18 shows that the theoretical curves calculated by formulas (104), (160), and (180) are in good agreement with the experimental data borrowed from [53].The mean distance  between the seats in a channel was estimated proceeding from the value corresponding to a close packing of a unit cell of a ZSM-5 membrane with nitrogen (Section 2.1).
The filling factor for mixtures with a strongly absorbed component depends on the mixture composition; it is significant and can be as high as  ∼ 1.In this case, in accordance with the arguments put forth in Section 2.4.1, the dependence of fluxes on the mixture composition is nonmonotonic due to the formation of clusters in a channel.
Figure 19 shows the dependence of the ratio of partial fluxes of gases for the mixture -C 4 H 10 -CH 4 on the butane concentration in the mixture under a pressure of  = 100 kPa at temperature  = 300 K.It can be seen from the figure that the ratio of partial fluxes increases with the concentration of the highly adsorbed component (butane) in the concentration range 0 <  < 0.35.Analysis of experimental data [51] and calculations based on formulas (104), (160), and (180) show that the methane flux for a butane concentration of  > 0.05 changes insignificantly upon an increase in the butane concentration in the mixture.In this case, the behavior of the flux ratio is completely determined by the behavior of the butane flux.The increase in the butane flux is due to the fact that, in accordance with formula (141), short-lived butane clusters are formed for concentrations 0 <  < 0.35 and for the fill factors corresponding to these concentrations (see Figure 13(a)); the size of these clusters increases with the concentration of the highly adsorbed component.The transport over clusters (see Figure 13(b)) occurs via barrier-free transfer of density excitation, leading to an increase in the effective diffusion coefficient for butane, which is observed for  < 0.35.The increase in the butane concentration in the mixture increases the lifetime of the clusters formed in the system (see Figures 14(b) and 15(c)).With increasing lifetime, the transport over the clusters slows down and the butane flux decreases for concentrations  > 0.35.A further increase in the concentration results in the formation of stable clusters in a channel (see Figures 14(c) and 16(c)), leading to a rapid (exponential) decrease in the butane flux upon an increase in its concentration in the mixture.
Thus, it follows from the theory that the butane flux becomes exponentially small for a butane concentration of  > 0.6 in the mixture, and the transport of particles must be blocked.It should be noted that calculations based on formulas (104), (160), and (180) for butane concentrations in the interval 0 <  < 0.05 show that the methane flux decreases approximately by an order of magnitude upon an increase in the butane concentration, which is also in qualitative agreement with the experimental data [51].The numerical difference between the predictions of the theory developed previously and the experimental data in this range of butane concentrations are due to the fact that relations (104), (160), and (180) were derived for the model potential of intermolecular interaction of the type of interaction (101) between hard spheres.It is well known that the potentials of hard spheres correctly describe the behavior of real systems for high densities [46,101].In order to describe the behavior of the system in the entire range of fill factors, potentials of the Lennard-Jones type, which take into account the attraction between particles even in the zeroth order in density, should be used.However, although the application of this type of potentials of intermolecular interaction does not change the results qualitatively, it leads to the formation of butane clusters at lower filling factors and, hence, to a sharper decrease in the methane flux upon an increase in the butane concentration in the concentration range under investigation.
Relations (104), ( 160), (180), and (181) also make it possible to derive the dependence of the selectivity on the pressure, temperature and composition of the mixture.Figure 20(a) shows the dependence of the ethane selectivity for the C 2 H 6 -CH 4 mixture with a concentration ratio of 50 : 50 at the entrance on the total pressure of the mixture [97].It can be seen that the dependence is nonmonotonic.This is due to the fact that the diffusion coefficient first increases due to an increase in the pressure and fill factor of the channel as a result of the formation of short-lived clusters, and the partial ethane flux increases (Figure 20(b)).The increase in the total fill factor of the channel is determined by the increase in the partial fill factors for both components.This process continues until the channel filling with ethane slows down.In this case, the fill factor of the channel increases mainly due to channel filling with methane.At such pressures, the increase in the ethane flux slows down (see Figure 20(b)), while the methane flux continues to increase (Figure 20(c)).The theoretical dependences plotted in the figures show that a further increase in pressure due to an increase in the lifetime of ethane clusters formed in the channel must lead to a decrease in the ethane flux (Figure 20(b)).Thus, the analysis of experimental data proves that the theory constructed here describes the experimental results satisfactorily.
One of the main results of the proposed theory is the prediction of nonmonotonic dependences of partial fluxes and selectivities on the mixture composition and pressure (see Figures 19 and 20), while conventional models [51,53,97] (in particular, the generalized Maxwell-Stefan equation [51,97]) lead to monotonic dependences of these quantities on the same parameters (see Figure 20(a)).The qualitative discrepancy between the dependences in question is due to the fact that conventional models [51,53,97] take into account only the finite size of particles, in which the interaction between particles in a channel is disregarded.It is well known [46,47,101], however, that the interaction between particles in dense systems plays a decisive role both in the construction of the equations of state of the system and in the description of transport.For example, the inclusion of interaction between particles of the hard-sphere type leads to the emergence of peaks in the pair correlation function at distances equal to one, two, three, and so forth particle diameters [44,46,101].This fact indicates the existence of "effective" attraction between particles and necessitates the inclusion of cluster formation in the description of the behavior of the system at high densities, when it passes to a spatially inhomogeneous state.It was demonstrated previously that this leads to the emergence of nonmonotonic dependences of fluxes and selectivities on external parameters.It should be noted that the stabilization of clusters in 1D systems by the second component at high fill factors, which has been investigated here, was obtained numerically in a recent work [83].It should be borne in mind, however, that, in accordance with our previous results (Section 2.1), clusters with a large but finite lifetime can be formed even in a onecomponent system with high filling factors.It is 1D gold clusters of this type that were apparently observed in recent experiments under the ultrahigh vacuum conditions.

Results and Discussion
The analysis of the experimental data showed that for most gases investigated, the sorption and transport can be treated as a process in a one-dimensional (1D) system.In statistical physics, the theoretical models of 1D systems are classical examples of exactly solvable models [44].It is well known [34,46] that in 1D systems, a phase transition (condensation) does not occur as density increases.This means that, specifically, no critical nuclei are present in the system, and stable nuclei of a new phase with macroscopically long lifetimes do not appear in the system.At the same time, the system of 1D channels in zeolite membranes consists of channels of finite length, where for a sufficiently high filling factors, clusters with sizes comparable to the channel length can form.Thus, the description of molecular transport in zeolite membranes reduces to describing transport in a 1D system, where there are strong density fluctuations with a finite lifetime of clusters.The density functional method [38] makes it possible to calculate the spectrum of density fluctuations and the diffusion coefficient for particles with arbitrary density and arbitrary laws of interaction with one another and with the channel walls.Since the flux measured in the experiments performed in [41] is determined by the diffusion coefficient and the filling factor of the channels in the membrane, the sorption isotherm of a one-component molecular gas is calculated separately.For the one-dimensional channels considered in the present paper, this problem is solved exactly for an arbitrary intermolecular interaction potential.A new feature here is an analysis of the dependence of the lifetime and size of the 1D clusters on the filling factor of the channel with particles.A description of molecular transport for an arbitrary filling factor () of a channel is performed.It was found that as the filling factor approaches 1, so that blocking of the relative motion of the particles becomes substantial, the diffusion coefficient increases.Following [31], a new diffusion mechanism is observed.It makes it possible to explain the transition from activated diffusion of single particles in a channel with low filling factors to collective fast barrier-free diffusion, which consists of the propagation of density disturbances for large values of  0 .In the hardsphere approximation, as  1 increases, the well-known [46] effective attraction of the molecules-hard spheres-results in vanishing of the energy barrier for diffusion along the channel axis.Another consequence of the effective attraction of molecules is the formation of molecular clusters in the channel, which have a finite lifetime because of the onedimensionality of the system.The size and lifetime of the clusters increase as  increases.The diffusion of particles in clusters is described as a barrier-free process of propagation of density disturbances.The pressure, temperature, and filling factor dependences of the flux of molecules make it possible, even in the high-sphere model, to describe the experimental data currently known to the authors.
We have shown that contradictory experimental data, reporting the presence of the SFD mode and transport acceleration with the increasing density, can be described within the proposed model consistently accounting for correlation effects in a strongly fluctuating 1D system of identical particles.The molecular transport in a 1D system can change with time going from the Einstein diffusion mode to SFD mode and then again to a normal diffusion at large observation times.Whereas the sequence of these modes is general, their actual realization depends on the concrete characteristics of the system.For example, at very low densities one expects to observe only a normal diffusion mode.If we consider a system of weakly interacting particles ( ≤ 1), we can see that at low densities  ≤ 0.2 and time scale  <  2 0 / 0 , the particles (3.5 ÷ 4.0 Å) or smaller Au-Au separations (2.5 ÷ 2.9 Å) can have an observable lifetime.Also we have estimated the minimal number of atoms in a 1D chain and a force required to break a chain with a minimal number of atoms in it for different metals (Au, Ag, Ni, and Pt).Thus we show that the methods of statistical mechanics, in some cases, can be used to investigate the properties of finite systems.In particular, these methods allow one to explain experimentally observed anomalously large interatomic distances in monoatomic Au chains and to predict its maximal length and to estimate its lifetime.The described approach was applied to analyze the interatomic distances and lifetimes of quasi-1D structures of different materials such as Ag, Ni, and Pt.The obtained results successfully describe the experimental data.

Conclusions
The paper is devoted to the analysis of the correlation effects and manifestations of general properties of 1D systems (spatial heterogeneity that is associated with strong density fluctuations, the lack of phase transitions, the presence of frozen disorder, confinement, and blocked movement of nuclear particle by its neighbours) in nonequilibrium phenomena by considering the four examples.The anomalous transport in zeolite channels is considered.The mechanism of the transport in carbon nanotubes and MOF structures, relaxation, mechanical properties, and stability of nonequilibrium states of free chains of metal atoms, non-Einstein atomic mobility in 1D atomic systems is considered.Description of molecular transport in zeolite membranes is reduced to the description of transport in 1D system in which strong density fluctuations with a finite lifetime of the clusters are presented.Response function method for density functional allows us to calculate the spectrum of density fluctuations and the diffusion coefficient of the particle with arbitrary density and arbitrary laws of interaction between particles and between a particle and a wall of the channel.Systematic account of correlation effects leads to the following physical picture of transport and mobility in strong fluctuating 1D confinement system.
Thus, if the interaction between particles in the 1D system is low ( ≤ 1), at low densities  ≤ 0.1 at times  <  2 0 / 0 particles move independently of each other, and, therefore, their mobility is of the Einstein type ⟨ 2 ⟩ ∼ .At the times  >  2 0 / 0 , the motion of the particles is limited due to blocking by other particles that results in the SFD regime ⟨ 2 ⟩ ∼  1/2 .1D system becomes spatially inhomogeneous due to collective effects in it by increasing of density of system.This state is characterized by the presence of clusters with a finite lifetime.In such a system, the diffusion of a particle is the result of its moving in a confined space between clusters with subsequent transfer of momentum to a cluster.This effect results in a disintegration of a cluster and a disconnection of the particles from the other side of the cluster.For identical particles, this process is a result of diffusion due to the "excitation" transfer along the chain of closely spaced particles and leads to an increase of a diffusion coefficient with increasing density of particles in 1D clustered system.
A further increase of the filling factor leads to the increase of a cluster lifetime, to the simultaneous decrease in the average distance between these clusters, and to the increase of diffusion coefficient.In this case, the dependence of the mean square displacement of particles at the large times ( → ∞) is of the Einstein type with an increased diffusion coefficient () ∼  0 /(1 − ) 2 ,  ∼ 1.If the pair interaction between particles is large, the cluster formation occurs at low densities, so that all the discussed above becomes true for small filling factors zeolite channels also.
Cluster of atomic particles that forms in 1D systems allows to describe and explain the phenomena observed in experiments, such as the formation of atomic chains with two different interatomic distances, mechanical properties, and stability of the nonequilibrium free chains, and allows to calculate the minimal number of atomic particles of different metals in free chains.In this case, it is not necessary to use impurity atoms, the presence of which requires experimental confirmation in the chain to explain the large interatomic distances in the 1D chains of gold.
The description of the practically important transport of two-component mixture in membranes with subnanometer channels by the response function method allowed to show that in 1D system of two types of particles at high filling factors, a phase transition in an inhomogeneous state with the formation of clusters of a finite size components may occur due to an effective particle attraction.
In the framework of the response function method, the equation of motion of the order parameter of the system is derived.The order parameter is the Fourier component of the density fluctuations.Analysis of this equation shows that two-component mixture in the channels passes into a spatially inhomogeneous state with the increase of  that results in formation of a short-lived clusters, the size and the lifetime of which increase with increasing of .
A new result here was the occurrence of a minimum of a dependence of fluctuation part of the free energy on the wave vector and the Fourier component of the order parameter.Thus, at high filling factors in the two-component mixture, there is a new feature: the clusters of a certain size in the channels stabilizes by a potential barrier due to the effective attraction of particles.
Such a situation is typical for the transition to an inhomogeneous state.The lifetime of the clusters increases exponentially according to the Arrhenius law, and the clusters can be locked for transport of the particles of the mixture at sufficiently low temperature channels.The observed increase of the separation factor with conductivity of 1D channels in case of blocking a channel by stable clusters with a certain concentration is described.The dependence of the separation coefficient on the concentration, temperature, and density is defined.
A comparison of the theory with experimental data and results is discussed.The resulting fluxes depending on the composition of the mixture (the filling factor) and selectivity depending on pressure describe the experimental data presented in the literature.

Figure 1 :
Figure 1: Reduced "pressure" ( 1D ) (24) of a one-dimensional gas in a channel versus the filling factor () of the channel at various temperatures  = / 0 .

Figure 9 :
Figure 9: The compressibility versus density for Au atoms at room temperature.

Figure 10 :Figure 11 :
Figure 10: Pair distribution of gold atoms at room temperature.

Figure 16 :
Figure 16: Dependence of free energy  1 on the order parameter for different values of the wave vector for  = 0.45:  = 0.1 (a), 0.3 (b), and 0.5 (c).

Table 1 :
Experimental and MD results:  and  are the particle and membrane channel diameters, respectively.

Table 2 :
Minimal chain lifetime , minimum number of atoms in a stable chain  min , and  0th ,  0 are the estimated and experimental values of the force required to break the cluster with a minimum number of atoms in it correspondingly.