Quantifying the potential and flux landscapes for nonequilibrium multiverse, a new scenario for time arrow

We propose a new scenario of nonequilibirum multiverse. We quantified the potential landscape and the flux landscape for the Bousso-Polchinski type of multiverse. The potential landscape quantifies the weight of each universe. When the terminal vacuum with zero (flat) or negative cosmological constant (AdS) have a chance to tunnel back to the normal universes with positive cosmological constant (dS) through the bounce suggested by the recent studies, the detailed balance of the populations of the multiverse can be broken. We found that the degree of the detailed balance breaking can be quantified by the underlying average flux and associated flux landscape, which gives arise to the dynamical origin of irreversibility and the arrow of time. We also showed that the steady state of the multiverse is maintained by the thermodynamic cost quantified by the entropy production rate which is associated to the flux. This gives arise to thermodynamic origin of time irreversibility. We show that terminal vacuum universes can have dominant weights or lowest potentials giving arise to a funnel shaped potential landscape, while terminal vacuum universes together with other normal universes including ours can form dominant cycles giving arise to a funnel shaped cycle flux landscape. This indicates that even our universe may not be distinct from others based on the probability measure, it may lie in the dominant cycle(s), leading to higher chance of being found. This may provide an additional way beyond the anthropological principle for identifying our universe.


I. INTRODUCTION
The conventional Big Bang cosmology theory faces difficulty for interpreting the homogeneity problem, the flatness problem, and the magnetic monopole problem and so on. Inflation theory was suggested to resolve these issues [1]. The essential picture of the inflation is that the universe went through a transient fast accelerating expansion phase soon after the Big Bang. The inflation is often thought as being driven by the vacuum energy/dark energy. If there are different vacuum states, different vacuum can tunnel to each other [2]. Following the inflation of the universe, the process of vacuum decay is reminiscent to the bubble formation in vaporization of liquid to gas phase transitions. However, there are issues associated to this old inflation model.
Although the false vacuum can decay to the true vacuum through the bubble nucleation, They cannot have the chance to collide and reheat the universe since they cannot catch up with the expansion of the rest of the inflationary universe. For this, a new inflation theory was proposed to have a slow and continuous transition from the false vacuum to the true vacuum within a single bubble (universe) [3]. The inflation theory predicted the amplitudes and the fluctuation spectrum of the cos- * Email: jin.wang.1@stonybrook.edu mic background radiation as the seeds for large structure formation in our universe consistent with the observations.
Furthermore, it was pointed out that the inflation can continue without ending in most part of the universe, giving arise to new bubble universes one after another. Thus, the inflation theory gives arise to a very different picture of the evolution of the universe, leading to a multiverse with eternal inflation [3,4].
On the other hand, in string theory, there are huge number of different kinds of universes with different vacuum and different coupling constants. In fact, there are estimated of 10 500 of such universes or vacuum. Some studies suggest that there is no preferred universe and all universes should be treated on an equal footing [4,5].
There are immediate questions about how our universe can be identified from the enormous amount of the possible other universes [4]. There is another closely related issue. According to the quantum field theory, the vacuum energy is huge. But according to the cosmological observation, the cosmological constant is tiny. The associated issue is why the cosmological constant of our universe is so small out of so many possibilities of the different universes [4]. Anthropological principle has been used to explain the observational existence of our universe [4]. For example, if we have somewhat different universe with different cosmological constant and coupling constants, the galaxies will not be formed properly, our human be-ing will not come to existence to be able to observe our current universe. Therefore, according to the anthropological principle, we happen to live in a universe which can produce ourselves and the observational universe we currently see. This has been used to explain why the cosmological constant of our universe is so small [4].
Despite of progresses being made, there are still challenging issues related to the multiverse picture. One is related to the observational evidences. Suggestions have been made to study the emergence of the black holes as the trace of the possible consequence of the collisions of different universes [6]. The clear evidences remain to be seen.
Another important issue is related to the time arrow. How can we determine the direction of time in our universe or the multiverse? Suggestions have been made to explain the time arrow by including the terminal universes with negative or zero cosmological constants [7]. These terminal universes cannot convert to the universes through the vacuum tunneling. But the other universes can tunnel to these terminal universes [7]. Therefore, the terminal universes act like sinks in the multiverse. This produces the irreversibility and the arrow of time. However, the time arrow can only last at finite time during the multiverse evolution. In the long time limit, the time arrow ceases to exist.
Recent studies suggested the possibilities of bounces avoiding singularities from a contracting universe (with negative cosmological constant) back to the expanding universe (with positive cosmological constant) [8]. Furthermore, there were also discussions on the possibility of expanding universe (with positive cosmological constant) born from the flat universe (with zero cosmological constant) [9][10][11][12][13][14]. If one takes this into consideration, the terminal vacuum (with zero or negative cosmological constant) will have small chances of coming back to the normal vacuum (with positive cosmological constant). The detailed balance can be broken in this case [8,14]. This suggests a need to study the multiverse problem from a nonequilibrium perspective.
In this work, we will study the multiverse evolution from the perspective of nonequilibrium physics. Without the assumption of the detailed balance among the vacuum tunneling switching of the universes, we show that the multiverse evolution is driven by two forces. One force is the underlying probability landscape of the multiverse quantified by the steady state probability distribution in multiverse state space. The other is the steady state probability flux which quantifies the degree of the detailed balance breaking. While the landscape attracts the multiverse down to their steady state basins of attractions, the flux provides a driving force for the cycle flow in multiverse state space. The steady state probability forms a potential or weight landscape while the steady state flux forms a flux landscape of cycles. We found that while the terminal vacuum has largest weight and therefore lowest potential, the largest flux cycles can include not only terminal vacuum but also normal vac-uum with positive cosmological constant including our universe. Therefore, although our universe might not be the highest probability one, it may lie in the highest probable flux cycle, which still gives a higher chance to be observed. On the other hand, the explicit detailed balance breaking characterized by the magnitude of the flux through the possibility of terminal vacuum tunneling back to normal universe provides a source of irreversibility and therefore the arrow of time at all times (even at long time limit). Moreover, the dynamical driving force in terms of flux gives arise to the entropy production rate. In fact, the entropy production rate provides the thermodynamic cost for maintaining the steady state of the universe, which is also an indication of the time arrow.

A. Multiverse description
To explore the multiverse evolution, let us first consider the evolution of the universe. In general, we can use Einstein's general relativity to describe the evolution of universe. Let us assume that the universe is homogeneous and isotropic on the large scales and treat matter as the perfect fluid [15]. From the Einstein's equations, we can obtain the Friedmann Equations (flat universe) for the evolution of the universe [15]: where a is the expansion factor and H is the Hubble constant representing the rate of expansion, ρ is the proper energy density and P is the pressure. If we assume that the evolution of the universe has an inflation driven by a constant vacuum energy (Λ), then H = Λ/3, a(t)∝exp(H t). The universe rapidly expands. Suppose there are two adjacent places, Inflation will drive them away from each other. Beyond H −1 , there is no causal relation [3]. Within the universe (H −1 ), there are chances of the birth of another universe through quantum tunneling such as the ones illustrated from brown to yellow in Fig 1 [6]. The yellow universe also grows rapidly, but it never catches up with the brown universe. This process continues with the birth of more other universes represented by different colors shown in Fig 1 [6]. As a result, this eternal inflation leads to the multiverse landscape, a foam of expanding bubble universes within bubble universes.
It was suggested that the cosmological constant Λ has a discrete spectrum based on the string theory [16], Here Λ bare is the bare cosmological constant. q m is the charge of the string compactified vacuum flux. The range of n m is an integer assumed to be within N is a positive integer. The state characterized by the cosmological constant Λ i has a degree of degeneracy D(i ) [17]: We often view the cosmological constant as the vacuum energy. Different cosmological constants Λ i correspond to different vacuum states. These different vacua can be classified into three categories. The vacuum with positive cosmological constant Λ i >0 is called the de-Sitter (dS) vacuum. The vacuum with zero cosmological constant Λ i =0 is called the Minkowski vacuum. The vacuum with negative cosmological constant Λ i <0 is called the antide-Sitter (AdS) vacuum. In this study, we use the Greek letters to label the dS vacuum and Latin letters to label any vacuum state. There are different definitions to describe the fractions of different vacua [19]. For simplicity, we consider the fraction of comoving volume (P i ). For the description by the fraction of proper volume, see [19,20]. The probabilistic evolution of the multiverse containing various vacua for dS vacua is determined by the master equation [14,17,21] where κ αβ is the transition rate or the transition probability per unit time for an observer who is currently in vacuum β to find himself in vacuum α, it is determined by [8,14,17,21] and Γ αβ is the tunneling transition rate per unit physical space time volume between two different dS vacuum, it is determined by [2,8,17,18] Γ αβ = e −S(β→α)+S β .
Here S β = −8π 2 /H 2 β is the Euclidean action for the vacuum state Λ β and S(β → α) is the Euclidean action for the tunneling trajectory β → α which satisfied S(β → α) = S(α → β) [2,8,17,18]. The physical meaning of equation (7) is clear. The change in the chance or the probability of the universe α being observed will be determined by the input from the chances of other universes decaying to the current universe (α) subtracting the output from decaying of the current universe (α) to the other universes. Considering the state with cosmological constant Λ α under the degree of degeneracy of D(α), the equation (8) should be modified as One can easily show that is the steady state solution of the master equation (7) and From equation (11), we can see that the factor S β→α in Γ αβ acts as a scale factor constant and therefore is not important. In addition, we know that S β = −8π 2 /H 2 β and H β = Λ β /3 [17], So we can simply set If all the universes are dS vacuum, the detailed balance is preserved. The resulting steady state is an equilibrium state. Since this is an equilibrium state, there is no emergence of time arrow. To resolve this issue, one can explore the effects of the emergence of the Minkowski vacuum or AdS vacuum.
The general master equation for the probability evolution of the multiverse involving all possible universes (dS, AdS and Minkowski) is given by In previous studies [7,17,21], the tunneling transition rates from Mikowski and AdS universes to dS universes are assumed to be zero. Recent studies show that there is a possibility of the tunneling transition from Minkowski or AdS vacuum to dS universes through the bounce [8][9][10][11][12][13][14]. For this study, we simply assume that: In equation (15), state "1" represents the AdS vacuum state. State "2" represents the Minkowski vacuum state. The Latin letters "i" represent any vacuum state. S 1,2,3,4 are certain small constants. In addition, we assume that the formula (13) can be generalized to

III. RESULTS
A. Potential landscape and flux as the driving forces for the evolution of multiverse 1. Potential landscape and flux decomposition for the driving force dictating the evolution of multiverse.
We can write the master equation in (14) for determining the probabilistic evolution of the multiverse involving various vacua in the form of [22] d P dt where Here M is the transition rate matrix representing the transition rates from one state (universe) to another while P (i) represents the fraction of comoving volume of vacuum Λ i . For the steady state solution (P ss i ), M P ss = 0, P ss represents the steady state probability distribution of the vacuum states of the multiverse. We define the potential landscape as [22,23] The potential landscape U can attract the multiverse to the steady state.
As seen clearly, the transition rate matrix M determines the evolution of the probability dynamics of the multiverse. We can decompose the driving force into the following form through the symmetrization and antisymmetrization decomposition [24]: and therefore where Here F ss ji = M ij P ss j − M ji P ss i is the local steady state probability flux between i and j. If all the local steady state flux is zero between any i and j states, then there is no net input or output to or from the system. Thus the detailed balance is preserved and the system is in equilibrium state. On the other hand, if any local steady state flux is not zero, then there is a net input or output to or from the system. Thus the detailed balance is broken and the system is in nonequilibrium state.
From this decomposition, we can see that the driving force for the probability dynamics is determined by two parts. ∆ is time-reversible and the detailed balance preserved part of the driving force. It is determined by the potential landscape difference or gradient. This is analogous to the usual equilibrium dynamics where the driving force is dictated by the gradient of the potential. Θ is time-irreversible. It is determined by the steady state probability flux. The steady state probability flux can be used to measure the degree of the detailed balance breaking and thus the degree of the time-irreversibility. Thus the flux quantifies the nonequilibrium (detailed balance breaking) part of the driving force of the probability evolution dynamics. If the system is in equilibrium state, Θ must be zero and equilibrium dynamics is determined by the landscape gradient alone. However, if flux is not zero, then the nonequilibrium dynamics is determined by both the landscape gradient and flux. This decomposition is in the discretized representation of the state space.
In the continuous state space, the probability evolution can be described by the Fokker-Planck equation, we can decompose the driving force into a landscape gradient and curl steady state probability flux [25][26][27]. Compared to the continuous representation, ∆ is related to the landscape gradient and Θ is related to the curl flux. These two parts determine the evolution of the multiverse. While the potential landscape attracts the multiverse into certain states and provides the stability of these attractors, the flux provides a driving force to form stable flow in the multiverse state space.
We can also decompose the driving force of the probability dynamics for the transition matrix M into two components in another mathematical rigorous form as [22]: where One can easily prove that This demonstrates that we can decompose the driving force (M ) for the probability evolution into two parts. One part preserves the detailed balance (D), this force is time-reversible, another part breaks down the detailed balance (C ) and this part is time-irreversibility [22]. Let us define [22] F ji ≡ P ss j M ji − P ss i M ij ≡ P ss j C ji − P ss i C ij (30) as the steady state flux between state i and j.
In the next, we will show that the steady state probability flux can be further decomposed to certain cycles to form the cycle landscape of the steady state probability flux.

Cycle fluxes forming flux landscapes.
Let us look at the flux component of the driving force for the evolution of the multiverse. We can define the flux directly from the definition of the master equation while the probability flux F ji is defined as The master equation can be interpreted as the local conservation equation for the probability. The evolution of the probability is equal to the net flux in or out.
On the other hand, at steady state, dP i /dt = 0. If F ji = 0, the net flux is zero. This corresponds to the detailed balance and equilibrium situation. If the F ji is not equal zero at the steady state, then the presence of the local net flux to the system indicates that the detailed balance is broken. The steady state probability flux breaking the detailed balance becomes F ss ji = P ss j κ ij − P ss i κ ji .
Since both F ss ij and F ss ji refer to the same net local flux, we can delete this kind of redundancy. For simplicity, one can delete the one which is smaller than zero to reach the following definition [28]: One can easily prove that The definition (34) included all net local fluxes without redundancy. One can easily prove that Except for the null matrix, for any matrix, if it has the property (36) and (37), then one can always decompose the flux into the flux cycles (or flux loops) [22,28]. The procedures are outlined as follows.
Suppose J i1i2 > 0, then from (37), one can find One can easily prove that J ij still has the property (36) and (37). Therefore, one can repeat the above process and get α 2 , r (2) , J (2) , α 3 , r (3) , J (3) ,... until J (M +1) = 0, where M is a finite positive integer. We now have Therefore, from the global perspective, the steady state flux can be decomposed to many cycles or loops of circulation fluxes. This forms the nonequilibrium flux landscape.
B. The potential landscape in the comoving coordinate of multiverse According to our definition of the potential landscape, the potential landscape spectrum is shown in Fig 2, Fig  2 is plotted under the choice of the parameters for the string vacuum and transition among AdS, Minkowski and dS vacuum as J = 7, q 1 = 2.5, q 2 = 5.185, q 3 = 5.155,  The rational of the parameter choice is as follows: First, we need to specify AdS and Minkowski vacuum states, respectively. Second, the transition rates from AdS or Minkowski vacuum to the dS vacuum is assumed to be small compared with the transition rates in the opposite direction. Lastly, we assume that the transition rate from Minkowski vacuum to AdS vacuum is larger than the rate in the opposite direction.
If we arrange all these states of the multiverse on a 2 dimensional plane, the potential landscape can be shown in We notice that there is a characteristic of the energy band (Fig 2). The reason for the presence of this band is that the string theory inspired parameters q 2,...,7 are close to each other. When these string theory implied parameters become closer, the resulting characteristics of the band become clearer. Under these string parameters, the cosmological constant spectrum has the characteristic of the band. This can lead to the band structure in the potential landscape.
If there is a significant gap between the minimum and average of the potential landscape spectrum compared with the dispersion or variation, the potential landscape topography can be biased towards the dominant vacuum at the bottom. One can use a characteristic ratio for landscape topography to measure the degree of this bias [22], defined as Here < U > is the average value of the potential landscape, < U 2 > is the average value of the square of the potential landscape, and U m is the lowest potential energy value. While the numerator represents the gap, the denominator represents the fluctuation characterized by the standard deviation.
RR thus represents the ratio of the gap between the lowest and the average potential of the universes against the fluctuations characterized by the standard deviation through the variances. While the gap can be viewed as the slope or the bias towards the dominant vacuum, the fluctuations through the variances can be viewed as the measure of the roughness or traps of the potential landscape. A large ratio of RR indicates a landscape with large bias towards the dominant universe against the roughness or traps. That is the landscape of the multiverse has a funneled shape towards the bottom terminal vacua (AdS and Mikowski). It also shows that the dominant state sitting at the bottom of the landscape is distinct and discriminant from the rest of the universes and therefore stable. In this sense AdS and the Minkowski vacua (universes) are usually more stable than the dS vacua (universes).
We found that when N =2, RR(U ) = 3.8166 and when N =3, RR(U ) = 12.6087. These two values are significantly larger than 1. This indicates that the AdS vacuum is dominant and stable against others. In this figure, the blue curve is for N =2 and the red curve is for N =3. We can see that when S 2 becomes smaller, the RR(U) at first do not change significantly and then becomes larger. This indicates that when S 2 become smaller, the dominant vacuum becomes more distinct and discriminant from others and therefore stable. It is worthwhile to point out that our numerical simulation shows that when S 2 changes from 10 −6 to 10 −14 , the dominant vacuum is not always the AdS vacuum. In addition, we noticed that these two curves with N =2 and N =3 have similar trends against the vacuum transition rates from AdS to the other universes.
Here, V f represents the value of the flux in a loop. The red line represents the cycle with the dominant flux.
There are altogether 5073 flux loops or cycles in this multiverse connecting different universes together. The flux landscape is illustrated in the right panel of the Fig  5. For the purpose of clear view, we do not present all the flux loops but only a few dominant ones. Here, different nodes represent different states or universes. The thickness of the arrows represents the magnitude of the circulation flux in the loop. The size of node represents the weight or the steady state probability of the node or the specific universe. The dominant black node represents the AdS vacuum. The second dominant black node represents the Minkowski vacuum. Other small gray nodes represent the dS vacuum. We noticed that although there are 128 states (N =2), there are only 127 nodes. The third state does not appear in any loop. This is not strange since under these parameters,J i3 = 0, for any "i". Therefore, the third state cannot appear in any loop. The cycle involving the red arrows represents the dominant loop.
Since the dominant flux loop stands out from the rest of the others, it represents a limit cycle oscillation in multiverse state space. The orange arrows represent the second dominant loop and the green arrows represent the third dominant loop, while the blue and purple arrows represent the fourth and fifth dominant loops, respectively. The second, third, fourth and fifth limit cycle oscillations can also emerge but with much less chances than the dominant one.
This gives us a new angle of looking at the organization of the universes. The universes in the multiverse forms a network. Through the tunneling connections, the network of the universes is organized in a hierarchical fashion with the bottom layer or ground flux state as the dominant flux cycle connecting certain universes to- gether. Then the hierarchical structure towers are built up by the subsequent less dominant flux cycles layer by layer or excited flux states. This new structure of the network universe indicates that the universes may emerge or function in a clustered fashion in the form of cycles.

Robustness of the flux landscape
Similar to the definition (40), we can use the definition [22] to measure the shape of the flux landscape. < U F > represents the average value of the flux landscape while < U F 2 > represents the fluctuations as the average value of the square of the flux landscape. The U F m represents the smallest value of U F . The RR is defined as the ratio of the gap between the minimum (dominant) flux potential and the average of the flux potential against the standard deviation characterizing the fluctuations in the fluxes. A high value of RR indicates a high discrimination of the dominant flux cycle/loop against the other ones, leading to the limit cycle oscillations among the universes in the dominant flux cycle/loop. In Fig 5, RR(U F ) = 11.1334, giving arise to a distinct dominant oscillation flux cycle/loop. Although the numerical results are different with different input parameters, the flux landscape can always emerge, as long as the detailed balance is broken and S 2 = 0. If S 2 = 0, this corresponds to terminal vacuum without the transitions to the other universes. Then the irreversibility comes from the time evolution and ceases to exist at long time steady state.
The flux landscape provides a new scenario and perspective for the evolution of our universe. First, due to the detailed balance breaking, the emergence of the flux gives arise to another driving force for the evolution dynamics of the multiverse. Second, the steady state flux can be decomposed to directional loops or cycles of fluxes. This generates the irreversibility and therefore the direction of time due to the intrinsic nonequilbrium nature of the detailed balance breaking. Note that this mechanism of time arrow is quite different from the one suggested before on the time evolution under terminal vacuum assumption where at the long time limit the evolution reaches the equilibrium preserving the detailed balance. Third, another distinct feature is that the dominant steady state probability flux cycles or loops with higher chances being seen not only involve the universes with negative or zero cosmological constant but also with the normal universes with positive cosmological constant. This indicates although the weight of our living universe may not stand out from the rest of the universes in the multiverse, it can be involved in a dominant flux cycle which increases its chance of being observed. This may help to pick up our living universe with less dependence on the anthropological principle for resolving the cosmological constant problem. Fourth, the global nature of the multiverse evolution dynamics is determined by both the underlying potential landscape and flux landscape. Fifth, the flux landscape gives arise a hierarchical structure of the organization of the universes in terms of flux cycles (or loops). Sixth, this gives a new picture of the universe evolution. From the viewpoint of the dominant cycle, the birth of our universe may be due to the transition from another universe and the death of our universe may be due to the transition to another universe. The birth and the death of the universe occur in a periodic fashion. In this sense our universe never dies. It repeats itself over and over again at a fixed amount of times by emerging from or transforming to the other universes on the dominant cycle, much like the case of a biological cell cycle [22]. From this perspective, our universe is eternal.

D. The irreversibility, the thermodynamic dissipation and arrow of time
To address the irreversibility or time arrow, let us study the irreversibility of a particular trajectory. If there exists the irreversibility of any one of the trajectories, then there is the irreversibility of the whole system. Given a particular trajectory specified as We consider the following functional [29] R = R({σ}, P (σ 0 ), P (σ n )) = ln[ κ σn,σn−1 ...κ σ2,σ1 κ σ1,σ0 P (σ 0 ) κ σ0,σ1 ...κ σ1,σ2 κ σn−1,σn P (σ n ) ] to measure the irreversibility of the trajectory (43) since this gives the ratio of forward transition probability against the backward transition probability. An equal forward and backward transition probability will result in an equilibrium with a zero value in R. This indicates the time reversibility. The magnitude of R gives a measure of how forward rates are different from the backward rates in time. Therefore, R is a measure of irreversibility or time arrow. Fig 6 shows the value of R or irreversibility versus the variation of the parameter S 2 for some trajectories where EPR is the entropy production rate [31]. the precise relation between EPR and the irreversibility revealed by the famous fluctuation theorem [32].We noted that This is the connection between the EPR and the steady state probability flux. The steady state flux gives the dynamical origin of the nonequilibriumness through the detailed balance breaking, while the steady state EPR provides the thermodynamic origin of the nonequilibirumness for maintaining the steady state. We can define that For the potential landscape in Fig 2, when N =2, < R > 1 = 2.3747 × 10 −9 ; when N =3, < R > 1 = 1.4419 × 10 −8 . Again, the entropy production rate is not zero indicating the nonequilibrium thermodynamic cost is needed to maintain the irreversibility and time arrow. Fig 7 shows the first (average) and second order (fluctuation) of < R > versus the variation of the parameter S 2 . Larger transition rate from AdS to other universes will give arise to more degree of detailed balance breaking The green curve is < R >2 versus the variation of -lg(S2) for N =2 multiverse. The red curve is < R >1 versus the variation of -lg(S2) for N =3 multiverse. The black curve is < R >2 versus the variation of -lg(S2) for N =3 multiverse. and therefore the arrow of time. The blue dotted curve and the green dash-dotted curve represent the < R > 1 and < R > 2 respectively for N =2. Similar two curves representing the < R > 1 and < R > 2 respectively are plotted for N =3.
Fig7 shows that the following approximation is a good one: Therefore, we can use both < R > and EPR to describe the irreversibility of the system. In other words, < R > and EPR can provide the thermodynamic measure of the arrow of time in comoving coordinate. Here, the origin of the time arrow is the detail-balance breaking. This is different with the time arrow in [33] which comes from the dynamic evolution of our universe. We should point out that when S 2 is strictly equal to zero, the AdS vacuum and the Minkowski vacuum become the terminal vacuum. They cannot decay to another vacuum. In this case, the formula (51) or (52) cannot be directly used to describe the arrow of time. Because in this case, the fraction of the AdS vacuum is 1 and others are zero. The formula (51) and (52) are ill-defined.

IV. CONCLUSIONS
In this work, we studied the nonequilibrium evolution of the multiverse in the comoving coordinate in detail. We uncovered that the driving force of the multiverse evolution in general is determined by the global landscape quantified by the steady state probability distribu-tion and the steady state probability flux. While the landscape attracts the multiverse to the steady state basins, the flux gives arise to cycles and loops associating certain universes together. The emergence of the flux loops provides a signature and a quantitative measure of the degree of detailed balance breakingnonequilibriumness. In other words, the landscape of the multiverse is uncovered and quantified in this work. However, the multiverse evolution dynamics is driven by not only the landscape but also the flux. This is very different from the conventional picture of the dynamics determined by the landscape gradient for equilibrium systems under detailed balance.
The vacuum of the universes with positive cosmological constant can be transformed from each other through tunneling. Recent studies show the possibilities of the bounce of contracting universe with negative cosmological constant back to the expanding universe with positive cosmological constant [9][10][11][12][13]. This can lead to the breaking of the detailed balance [8,14]. The flux associated with the detailed balance breaking provides the dynamical origin of the irreversibility and time arrow. On the other hand, the entropy production is associated with the flux. Thus, it provides the thermodynamic origin or cost for maintaining the irreversibility and time arrow. In contrast to the time evolution argument of the underlying detailed balance system with terminal vacuum for generating the irreversibility, the current approach emphasizes the possibility of the time arrow generated by the detailed balance breaking even at long time steady state. This provides a new mechanism for the time arrow for the multiverse.
Furthermore, the potential landscape shows a funneled shape dominated by the contracting universes with negative cosmological constant and flat universes with zero cosmological constant. On the other hand, the flux loops form a flux landscape. The flux landscape also shows a funneled shape dominated by certain flux cycles or loops. The dominant flux loops or cycles can involve not only the contracting and flat universes, but also expanding universes with positive cosmological constant. The dominant loop gives arise to the associations of certain different kinds of the universes together and this leads to the oscillation cycles among these universes in the multiverse evolution. In other words, a universe on the dominant cycle can appear due to the transition from the other universes on the same cycle and can also disappear due to the transition to the other universes on the same cycle. This occurs periodically. If our universe is on this dominant cycle, then it can be born due to the transition from another universe and die due to the transition to another universe in a periodic manner. In this sense, our universe never dies. The birth and death of our universe go through cycles. Therefore, the universes on the dominant cycle oscillate from one universe to another and repeatedly appear and disappear under a fixed period of times.
Moreover, although our living universe with small pos-itive cosmological constant does not necessarily have significantly higher probability compared to others in the multiverse, it can lie in the dominant flux cycles. If that is the case, then the chance of being observed can be significantly enhanced. Thus, this may provide a boost in addition to the anthropological principle for selecting our living universe.
The presence of the steady state probability flux breaks the detailed balance. This leads to the multiverse as an intrinsic nonequilibirum system. Even in steady state, the equilibrium is not reached since the detailed balance is not preserved and there is a net flux. In other words, the time arrow originated from this intrinsic nonequilibriumness exists at all times (even at the long time limit). This is in contrast to the case where the time arrow is generated during the evolution of the multiverse and ceases to exist at the long time limit.