Motility-induced coexistence of a hot liquid and a cold gas

If two phases exist at the same time, such as a gas and a liquid, they have the same temperature. This fundamental law of equilibrium physics is known to apply even to many non-equilibrium systems. However, recently, there has been much attention in the finding that inertial self-propelled particles like Janus colloids in a plasma or microflyers could self-organize into a hot gas-like phase that coexists with a colder liquid-like phase. Here, we show that a kinetic temperature difference across coexisting phases can occur even in equilibrium systems when adding generic (overdamped) self-propelled particles. In particular, we consider mixtures of overdamped active and inertial passive Brownian particles and show that when they phase separate into a dense and a dilute phase, both phases have different kinetic temperatures. Surprisingly, we find that the dense phase (liquid) cannot only be colder but also hotter than the dilute phase (gas). This effect hinges on correlated motions where active particles collectively push and heat up passive ones primarily within the dense phase. Our results answer the fundamental question if a non-equilibrium gas can be colder than a coexisting liquid and create a route to equip matter with self-organized domains of different kinetic temperatures.


INTRODUCTION
We are all used to the experience that a gas is often hotter than a liquid of the same material.For example, to evaporate water from a pot in our kitchens, we need to increase its temperature.Then, at some point, vapor molecules rapidly escape from the liquid and distribute in the surrounding air.This experience that vapor emerges when increasing the temperature of a liquid has played a key role throughout human history: It was an essential ingredient, e.g., for the development of the steam engine [1], and it is key to technological applications like distillation techniques [2,3] or physical vapor deposition [4,5] as well as to natural spectacles such as geysers [6].The central exception from the experience that gases are hotter than liquids of the same material occurs when two phases, e.g., a gas and a liquid, coexist at the same time.Then they share the same temperature.This is guaranteed by the fundamental laws of statistical mechanics and thermodynamics for all equilibrium systems and it also applies to some non-equilibrium systems [7][8][9][10][11].Intuitively, this is plausible since any type of temperature gradient seems to evoke an energy flow evening out an initial temperature gradient.Despite this, very recently, it was found that at phase coexistence in certain active systems consisting of active particles which consume energy from their environment to propel themselves [12][13][14], the dilute (gas-like) phase is hotter by up to one or two orders of magnitude compared to the dense (liquid-like) phase in terms of the kinetic temperature [11,15,16].(Here, the kinetic temperature is defined as the mean kinetic energy per particle, which is equivalent to other temperature definitions in equilibrium.)While this complies with our intuition that gases are often hotter than liquids, it is in stark contrast to the situation in equilibrium systems and the expectation that any temperature difference should evoke an energy flux that balances it out.By now, such a temperature difference across coexisting phases has been shown to occur for a variety of temperature definitions that all coincide in equilibrium.It has been observed, e.g., for the kinetic temperature [11], the effective temperature [15] as well as for tracer-based temperature definitions [16,17] in systems undergoing motility-induced phase separation (MIPS) [7][8][9][18][19][20][21][22][23][24][25][26][27][28][29][30], i.e., in systems of particles that self-organize into a dilute (gas) and a coexisting dense (liquid) phase.The mechanism underlying the emergence of a temperature difference across coexisting phases hinges on the consumption of energy at the level of the active particles when undergoing frequent collisions within the dense phase.This mechanism crucially requires inertia [11,15,16], whereas overdamped active particles show the same kinetic temperature in coexisting phases [11].The requirement of inertia restricts the observation of different coexisting temperatures to a special class of active systems and precludes its experimental observation in generic microswimmer experiments.
In the present work, we explore the possibility to achieve a kinetic temperature difference across coexisting phases in mixtures of two components that on their own would not lead to a temperature difference: an ordinary equilibrium system made of inertial passive Brownian tracer particles such as granular particles or colloidal plasmas and overdamped active Brownian particles like bacteria or synthetic microswimmers.Our exploration leads to the following central insights: First, we show that when the mixture undergoes MIPS, the passive particles in the dense and the dilute phase indeed can have different kinetic temperatures (and different Maxwell-Boltzmann temperatures, which are defined based on the width of the velocity distribution).This demonstrates that kinetic temperature differences in coexisting phases can occur in a broader class of systems than what was anticipated so far.Second, we find that not only the gas can be hotter -but, counterintuitivelyalso the dense phase can be hotter than the dilute phase.This appears particularly surprising since the current understanding of the mechanism leading to different temperatures across coexisting phases in pure active systems hinges on the idea that frequent directional changes due to collisions lead to a local loss of kinetic energy (similarly as inelastic collisions do in granular systems [31][32][33][34][35][36][37][38][39][40][41][42][43]).
Such collisions are more frequent in dense regions suggesting that the dense phase is always colder than the dilute one, which coincides with all previous observations [11,15,16,[31][32][33][34][35][36][37][38][39][40][41][42]44].We find that this mechanism also applies to the passive tracers of the active-passive mixtures studied here in a certain parameter regime in which the tracers are trapped within dense clusters by surrounding active particles leading to low tracer temperatures in dense regions analogously as in the singlespecies case of inertial active particles.However, surprisingly, we find that this effect can also be reverted in mixtures of active and passive particles.This is because for strong self-propulsion, the active particles persistently push passive ones forward even within the dense phase, which can overcompensate the slowing of the latter ones due to (isotropic) collisions with other particles.
Hence, the passive particles can achieve a higher (kinetic) temperature in the dense phase than in the surrounding gas, where correlated active-passive motions occur less frequently and last shorter.

Model
We study a mixture of active and passive particles in two spatial dimensions, in which the active (passive) particles are represented by the active (passive) Brownian particle [ABP (PBP)] model [11,26,30,[55][56][57].While the motion of the active particles is overdamped, the passive species is significantly heavier (inertial).For simplicity, we consider active and passive particles with the same size and drag coefficients [58] but different material density.However, note that the key effects which we discuss in the following are similar for particles with significantly different sizes and drag coefficients, as we shall see.The particles are represented by (slightly soft) spheres and the dynamics of the active particles is made overdamped by choosing a very small mass m a and a small moment of inertia I = m a σ 2 a /10 (corresponding to a rigid sphere).The active particles feature an effective self-propulsion force F SP,i = γ t v 0 p i (t), where v 0 , p i denote the (terminal) self-propulsion speed and the orientation p i (t) = (cos ϕ i (t), sin ϕ i (t)) of the i-th active particle (i = 1, 2, ..., N a ), respectively.The position r i and the orientation angle ϕ i of the i-th active particle evolve according to dr i /dt = v i and dϕ i /dt = ω i , respectively, where the velocity v i and the angular velocity ω i evolve as Here, T b represents the bath temperature, γ t and γ r are the translational and rotational drag coefficients, respectively, and k B is the Boltzmann constant.To have access to a well-defined instantaneous particle velocity, we explicitly account for inertia for the active species but choose a very small mass to stay in the overdamped regime.Notice that using overdamped Langevin equations instead (m a = 0) essentially yields the same results (Fig. S13, Supplementary Information).
numerically by using LAMMPS [60] (see Methods for details).Finally, we define the Péclet number, which measures the relative importance of self-propulsion and diffusion, by Pe = v 0 / √ 2D r D t , where denote the translational and the rotational diffusion coefficients of the active particles, respectively.

Coexistence of a hot gas and a cold liquid
Let us first consider an initially uniform distribution of an overdamped mixture of active and passive particles [58,61,62].In our simulations with Pe = 100, an area fraction of φ tot = 0.5, and a fraction of x a = 0.6 active particles, we observe that the active and passive particles aggregate and form persistent clusters despite the fact that they interact purely repulsively.These clusters are motility induced [7-9, 18-26, 29, 30] and continuously grow (coarsen), ultimately leading to a phase separated state comprising a dense liquid-like region that coexists with a dilute gas phase (Fig. 1a-d and Movie S1, Supplementary Information), which is in agreement with previous studies [58,61].As for systems of active overdamped particles alone [11], we find that the active and passive particles in both phases have the same kinetic temperature (shown in Fig. 1d for the passive particles).Here, following Refs.[11,16,63,64], we define the temperature of the particles based on their kinetic energy as k B T a/p kin = m a/p (v − ⟨v⟩) 2 /2, which is well-defined also in nonequilibrium systems [65].(Note that the phenomena which we report occur similarly if using other temperature definitions such as the Maxwell-Boltzmann temperature, as further discussed below.)Let us now explore if the situation found for the overdamped mixture changes when replacing the overdamped tracers with (heavier) underdamped ones (Fig. 1e-h).Then, at the level of the structures that emerge, not much changes in our simulations: We still observe the formation of small clusters, which is followed by coarsening, ultimately leading to complete phase separation.However, when exploring the kinetic temperature of the passive particles within the steady state, we find that, remarkably, the passive particles in the dense phase are colder than in the dilute phase.The temperature ratio of the two phases is highly significant and amounts to approximately 2.5 (Fig. 1h and Movie S2, Supplementary Information).While this temperature difference is similar to what has previously been seen in underdamped active particles [11,15,16] and driven granular particles [31][32][33][34][35][36][37][38][39][40][41][42]44], its emergence in the present setup is surprising since it is well known that neither the overdamped active particles [11] nor the underdamped tracers alone [66,67] would result in a kinetic temperature difference across coexisting phases.Accordingly, the temperature difference must arise from the interactions of the two species.To understand this in detail, it first might be tempting to start from the common understanding of kinetic temperature differences in granular systems or purely active systems made of inertial ABPs such as Janus colloids in a plasma [68], microflyers [69], or beetles at interfaces [70] which relates the emergence of a temperature difference to an enhanced energy dissipation in the dense phase at the level of the particles.The latter could occur due to inelastic collisions as for granular particles [39] or due to multiple collisions between which drag forces transfer energy from the particles to the surrounding liquid as for active particles [11,16].However, in the present case of a mixture, collisions between active and passive particles have a different effect in the dense and in the dilute phase.In the dense phase, the motion of the passive particles is constricted by the surrounding clustered ABPs (see, e.g., Fig. 1e), which accumulate mostly at the border of the clusters, similarly as in completely overdamped mixtures [58,61], and which cause an effective attraction between the passive tracers by pushing them together [49,71,72] (see also supplementary text and Fig. S12, Supplementary Information).Therefore, the passive particles cannot move much in the dense phase and have a lower kinetic energy there compared to the dilute phase.This is also visible in the velocity distribution of the passive particles, which narrows for increasing x a in the dense phase (Fig. S1, Supplementary Information).In contrast, in the dilute phase, when active particles collide with passive particles, they can persistently push passive particles forward and accelerate them such that their kinetic energy increases above the energy they would have due to the surrounding heat bath.Such a correlated active-passive dynamics (Fig. 4b) heats up the passive particles in the dilute phase and leads to a broader velocity distribution (Fig. 3c) at intermediate Pe such as Pe = 100.

Hot liquid-like droplets in a cold gas
Since the observed temperature differences are activity induced, one might expect that the temperature gradient further increases when enhancing the self-propulsion speed of the active particles, i.e., when increasing Pe.Surprisingly, however, in many cases, the opposite is true.For example, for fractions x a = 0.3, 0.6, or 0.9 of ABPs, we find that the kinetic temperature difference is largest for some intermediate Pe and then decreases essentially monotonously with increasing Pe (Fig. 2) before it even reverts and we obtain dense liquid-like droplets that are hotter than the surrounding gas.As time evolves, these droplets grow (coarsening) leading to larger and larger clusters, ultimately resulting in a single hot and dense cluster that persists over time.Exemplarily, we show typical snapshots for the case Pe = 400, x a = 0.9 in Fig. 1i-l (see also Movie S3 and Fig. S2, Supplementary Information).In panel l, one can clearly see that the liquid in the center of the figure is hotter than the surrounding gas.Such a coexistence of a hot liquid-like droplet and a cold gas -in terms of the kinetic temperature -is in stark contrast to what has been found for underdamped active particles [11,15,16] and for driven granular particles [31,32,[39][40][41][42].Note that for very large liquid-like droplets containing significantly more than about 10 4 particles, it may happen that not the entire droplets are hot but only a certain layer at their boundaries.The emergence of a hot dense droplet also contrasts with the intuitive picture given above that hinges on the idea that active particles can efficiently push forward and accelerate passive particles only in low density regions.Therefore, the key question that guides our explorations in the following is: What is the mechanism allowing for a coexistence of hot liquid-like droplets and a colder gas?
Mechanism: Correlated active-passive dynamics heats tracers in the dense phase.
We now explore the mechanism underlying our previous observations that in mixtures of overdamped ABPs and inertial passive particles dense liquid-like droplets are persistently hotter than the surrounding gas at large Pe.To this end, we first analyze the velocity distribution of the passive tracers in the uniform regime at x a = 0.2 and in the phase-separated regime at x a = 0.8, which broadens as Pe increases (Fig. 3a-c).Strikingly, if and only if the active particles are sufficiently fast (Pe ≳ 200), the velocity distribution broadens more in the dense phase than in the dilute phase (Fig. 3b-d).This means that increasing the speed of the active particles (i.e., increasing Pe) has a much stronger effect on the speed of the passive particles in the dense regime (where collisions are more frequent) than in the dilute regime, which ultimately leads to hot liquid-like droplets.What remains open at this stage is why the velocity distribution broadens faster for passive particles in the dense regime than in the dilute regime (only) if the Péclet number is large.
To answer this question, we now explore the power balance of the passive particles in the dense and the dilute phase.As we will see, this power balance will point us to correlations between active and passive particles which lead to hot liquid-like droplets at large Pe.To obtain a power balance equation for the passive particles, we first multiply Eq. ( 3) by v j and take the ensemble average.
in the steady state, where F int,j = − N n=1 n̸ =j ∇ r j u (r nj ) is the total interaction force on particle j and r nj = |r n − r j |.If we now compare the power balance for particles in the dense and in the gas phase, we can express the kinetic temperature difference as This central equation leads to two important conclusions: First, the kinetic temperature difference between the dense and the gas phase is proportional to m p /γ t , which vanishes if the passive particles are overdamped in accordance with our simulations (Figs.1d and S3, Supplementary Information).Interestingly, the same proportionality has also been observed for a single-component system consisting of inertial ABPs, where it has been observed that the dense phase is always colder than the gas phase [11].Second, the temperature difference depends on the interaction between the particles given by the term ⟨v • F int ⟩, which measures how strongly interactions push passive particles forward in their direction of motion.From the probability distribution of the individual values v • F int that contribute to the mean (Fig. 4a,e), we obtain significant differences between the dense and the gas phase at large values which determine the sign of the temperature difference: At intermediate Pe, e.g., Pe = 80 (Fig. 4a), large values of v • F int are more frequent in the gas phase than in the dense phase (see also Fig. 4d).That is, events in which the interaction force and the velocity of the passive particles are aligned and large (e.g., if an ABP is pushing a passive particle forward [73]) are more frequent in the gas phase than in the dense phase, in which the particles have significantly less space to move and accelerate.In contrast, at large Pe, such events are more frequent in the dense phase finally leading to the coexistence of hot liquid-like droplets with a colder gas (Fig. 4e,h).Intuitively, this is because at very large Pe, ABPs can (collectively) push passive particles forward over relatively long periods of time even in the dense phase without being stopped by collisions with other particles due to the strong effective self-propulsion force (cf.Movie S4, Supplementary Information).These correlated particle dynamics are exemplarily shown in Fig. 4b,f and schematically visualized in Fig. 5.The correlated dynamics of active and passive particles also lead to a long ballistic regime in the mean-square displacement of the passive particles at intermediate times (similar as for a completely overdamped mixture [61]) before the dynamics of the passive particles becomes diffusive again (Fig. S4f, Supplementary Information).
Finally, we can ask why the temperature difference between the hot liquid-like droplets and the cold gas is larger at large x a .This is because at large x a , the active particles accumulate in the dense phase and induce stronger collective motions in that place when they are many (see also Fig. S5, Supplementary Information).Conversely, the fraction of active particles in the surrounding gas does not depend much on x a , and hence, the collision rate in the gas does not increase with x a .

Non-equilibrium state diagram
Having seen that the coexistence of hot liquid-like droplets and a cold gas requires sufficiently fast self propulsion of the active particles, i.e., large Pe, we now examine the parameter dependence more systematically.Therefore, we explore the non-equilibrium state diagram by varying Pe ∈ [0, 400] and x a ∈ [0.0, 1.0] at a constant area fraction φ tot = 0.5.The transition line between the uniform and the MIPS regime is obtained by analyzing the local area fraction , where H is the Heaviside step function and σ j the diameter of particle j, calculated from averages over circles of radius R = 5σ.Its distribution is unimodal in the uniform regime and bimodal in the coexistence regime allowing to distinguish between the uniform and the MIPS regime [25,30,74,75] (Fig. S6, Supplementary Information).We distinguish between the passive particles in the dense and the dilute phase in the steady state and calculate their mean kinetic temperature (see Methods and Fig. S7, Supplementary Information, for details).
The system phase separates for large enough fraction of active particles x a and large enough Pe (Fig. 6).At small Pe, the transition line approximately follows the transition line of a purely overdamped mixture as obtained in Ref. [58], which reads x (critical) a ∝ 1/(φ tot Pe).However, at large Pe, the partially underdamped system requires a larger fraction of active particles to undergo MIPS than the purely overdamped system, which can be understood as a consequence of inertial effects: At large Pe, passive particles are typically fast when they collide with an ABP.Due to their inertia, the passive particles slow down only gradually and sometimes even push aggregated ABPs apart, which can destroy small aggregations.This effect is particularly pronounced for large Pe and opposes the onset of MIPS.Hence, compared to a completely overdamped system, a larger fraction of active particles is required to initiate MIPS at large Pe.
The different kinetic temperatures in the dense and the dilute phase are indicated by the colors in Fig. 6.It can be seen that the temperature difference between the dense and the dilute phase strongly depends on both x a and Pe: In accordance to the mechanism which we have discussed in the previous section, we find that for intermediate Pe, the dense phase shows a lower kinetic temperature than the dilute phase with a maximum temperature difference around (Pe, x a ) ≈ (80, 0.7) (red circle in Fig. 6).For large Pe and large x a , the kinetic temperature difference changes its sign, indicated by the squares in Fig. 6a, i.e., hot liquid-like droplets coexist with a cold gas.
The latter occurs at lower Pe for increasing x a because the overall energy transfer from the active to the passive particles is larger for large x a only in the dense phase, where the active particles increasingly accumulate as x a increases (see also Figs.S4 and S5, Supplementary Information).This can also be seen from the parameter dependence of the kinetic temperature of the passive particles (Fig. S8, Supplementary Information): The kinetic temperature increases with increasing x a (and increasing Pe) in the dense phase but shows a maximum at intermediate x a in the gas phase, where the fraction of active particles hardly increases when increasing x a beyond a certain point.

Role of inertia
Inertia of the passive particles is a key ingredient to observe coexisting temperatures.This can be seen in Fig. S3 (Supplementary Information) and from Eq. ( 5): The temperature difference is proportional to the ratio m p /γ t .Thus, in the overdamped limit m p /γ t → 0, the temperature difference vanishes (Fig. 1d) because the passive particles react instantaneously to acting forces.Thus, their motion, and hence, also their kinetic temperature, is dominated by diffusion [11].In contrast, sufficiently heavy (inertial) tracer particles can store the energy gained during collisions with active particles as kinetic energy such that their kinetic temperature is not determined by diffusion alone, which is fully consistent with our simulation data and previous literature [11,15,16,31,32,[39][40][41][42].Increasing inertia does also lead to a significant violation of the equipartition theorem both in the dense and the gas phase (Fig. S9, Supplementary Information), which indicates that the system is increasingly far away from equilibrium when increasing inertia of the passive particles.

Role of the particle size
For simplicity, we have considered active and passive particles with the same size and the same drag coefficients so far but with significantly different material density.Now, we show that persistent temperature differences also occur when the passive particles are significantly larger than the active ones.We have varied the size ratio s = σ p /σ a ∈ {1, 2, 3, 4, 6, 8, 10} keeping σ a as well as m p , m a , and φ tot fixed.The fraction x a is chosen such that the area fraction of active particles is approximately 0.5 for small size ratios.For large size ratios, we kept x a = 0.99 fixed to ensure that enough passive particles are inside the system.For the drag coefficient of the passive particles, we choose γt = σ p γ t /σ a .Our results are exemplarily shown in Fig. 7 for Pe = 100.Here, we observe a persistent kinetic temperature difference between the passive particles in the dense and the gas phase even for significantly different particle sizes (Fig. 7h).This temperature difference is also visible in the velocity distributions, which are broader in the gas phase compared to the dense phase (Fig. 7e,f).Hence, the observation of a cold dense phase that coexists with a hotter surrounding gas persists even for large size ratios.The opposite case, i.e., hot liquid-like droplets coexisting with a colder gas, is also robust and leads to a temperature difference of T for σ p /σ a = 10 at Pe = 400, x a = 0.99, and φ tot = 0.70 for example.Notice however, that very large passive particles tend to accumulate in the dilute phase especially at large Pe making it challenging to calculate a precise value of the temperature difference.Interestingly, this is in contrast to purely overdamped mixtures, where large size ratios support the formation of large passive-particle clusters [49,72].
How representative is the kinetic temperature?So far, following Refs.[11,16,63,64], we have used the kinetic energy of the particles to define a kinetic temperature as a measure for the temperature.The kinetic temperature has frequently been used for granular systems [31,32,[39][40][41][42][76][77][78] and is also well-defined in non-equilibrium systems [65].In equilibrium, the kinetic temperature is equal to the thermodynamic temperature [66].In the binary mixtures of active and passive particles studied in the present work, the kinetic temperature of the passive tracer particles, which measures the velocity fluctuations, has two contributions: one from the thermal Brownian motion and one originating from collisions with surrounding active and passive particles.From the previously discussed results, we know that the latter cause the kinetic temperature difference between passive particles in the dense and the gas and N a + N p = 5 × 10 4 for σ p /σ a = 6, 8, 10 (other parameters as in Fig. 2).
phase.Additionally, we analyzed the velocity distribution of the passive particles in the dense and the gas phase.The variance of this distribution (Fig. 3d) exhibits the same behavior as the kinetic temperature.Remarkably, the velocity distributions are approximately Gaussian for sufficiently large Pe (Fig. 3a-c).We exploit this to define a Maxwell-Boltzmann temperature T MB by fitting a Maxwell-Boltzmann distribution to the velocity distribution with one free fit parameter k B T MB .
For the data shown in Fig. 3b,c at Pe = 400 we obtain T MB /T bath = 3.6 × 10 2 (liquid-like droplets) and T MB /T bath = 1.9 × 10 2 (gas).This shows that mixtures of inertial tracers and overdamped ABPs can lead to self-organized hot liquid-like droplets that coexist with a colder gas also in terms of the Maxwell-Boltzmann temperature.
Since both the kinetic temperature and the Maxwell-Boltzmann temperature are sensitive to local collective motion patterns of the particles and a sensible measure for the temperature of the particles should measure their independent motion, we now explore if spatial velocity correlations of the passive tracer particles are crucial for the emergence of a temperature difference.For that we calculate the spatial velocity correlation function [79] As shown in Fig. 8b (and Movie S4, Supplementary Information), velocity correlations are indeed present between the passive particles in the dense phase over a significant spatial range indicating that collective motion might strongly influence the kinetic temperatures.In fact, we find that the mean distance between the passive particles in the dense phase calculated from a Voronoi tessellation is given by approximately 4.6σ for the case shown in Fig. 8a-c and therefore, much smaller than the length scale of the velocity correlations (Fig. 8b).To exclude that such collective motions are required to achieve a coexistence of a hot liquid and a cold gas, we have performed simulations with a very low fraction of passive particles such that their typical distances to each other are significantly longer than the velocity correlations within the dense phase.Concretely, we did a simulation with 10 5 particles and x a = 0.996 at Pe = 400, which again shows MIPS and a significant temperature difference between the passive particles in the dense droplets and the surrounding gas (Fig. 8d-f and Tab.S1, Supplementary Information).In contrast to the previous scenario shown in Fig. 8a-c, the correlations between passive particles are now significantly reduced, and the mean distance between passive particles in the dense phase is approximately 41σ, i.e., larger than the length scale of the velocity correlations (Fig. 8e).Hence, in this parameter regime, the temperature calculation is not much influenced by local collective motion of the passive particles, but remarkably, the passive particles in the dense phase still have a higher temperature than the passive particles in the gas phase.This is shown in Fig. 8f and in Tab.S1 (Supplementary Information).These results show that the coexistence of hot liquid-like droplets with a colder gas is really induced by the interactions between the active and the passive particles in the dense phase and should also occur for other (fluctuation-based) temperature definitions that are not sensitive to collective motions of the passive particles.
To explicitly see this, we additionally calculated a relative kinetic temperature by using the relative velocity of each particle to the mean velocity of particles in the vicinity where ⟨v⟩ R denotes the mean velocity of all particles in a circle of radius R = 5σ around the tagged particle.As shown exemplarily in Tab.S1 (Supplementary Information), the temperature difference is also visible for the relative kinetic temperature, and thus, it is not only a consequence of the observed collective motion but rather a pure effect of the particle interactions.As a result, the key phenomenon of the present work -the coexistence of hot liquid-like droplets and a cold gas -is robust with respect to the choice of the definition of the particle temperature.For a discussion regarding the role of the solvent we refer the reader to Ref. [16].

DISCUSSION
Mixing overdamped active Brownian particles and underdamped passive Brownian particles leads to a persistent kinetic temperature difference between the dense and the dilute phase when the system undergoes motility-induced phase separation.This temperature difference emerges despite the fact that each of the two components on their own would show a uniform temperature profile.Counterintuitively, the dilute gas-like phase is not always hotter than the dense liquidlike phase but at large Péclet number and fraction of active particles, hot liquid-like droplets can coexist with a cold gas.This temperature reversal results from the competition of two effects: The trapping of passive particles in the dense cluster provokes a cold liquid whereas the emergence of persistent correlated active-passive particle trajectories in the dense phase primarily heats up the liquid.While the latter effect has not been known in the literature so far, we have shown that it can even overcome the previously discussed trapping effect and lead to the coexistence of a cold gas and hot liquid-like droplets.This phenomenon is robust with respect to the choice of definition of particle temperature and particle-size effects at least up to a size ratio of 10.For even larger size ratios, it can happen that all inertial passive particles remain in the gas phase, and hence, no temperature difference can be observed.Besides their conceptual relevance, our results open a route to create a persistent temperature profile in systems like dusty plasmas or passive granulates by inserting overdamped active particles like bacteria, algae, or synthetic colloidal microswimmers.

Simulation details
The interaction between the particles is modeled by the Weeks-Chandler-Anderson (WCA) potential [59] u (1)-( 3)] are solved numerically in a quadratic box with periodic boundary conditions and with a time step ∆t = 10 −6 τ p using LAMMPS [60] first for a time of 100τ p to reach a steady state and afterwards for a time of 900τ p for computing time averages of observables in the steady state.

Kinetic temperature calculation for dense and dilute phases
To calculate the kinetic temperature of the dense and the dilute phase separately, we distinguish between passive particles in the dense and the gas phase by identifying the largest cluster in the system using the criterion that two particles belong to the same cluster if their distance to each other is smaller than the cutoff distance r c = 2 1/6 σ of the WCA potential.Then, all particles in the largest cluster are considered as the dense phase and all other particles as the gas phase (Fig. S7, Supplementary Information).Finally, the kinetic temperature of the passive particles in the dense and the gas phase is obtained by averaging over all passive particles in the dense phase and all passive particles in the gas phase, respectively.correspond to the two simulations shown in Fig. 8a-c and d-f in the main text, respectively.For the latter, the velocity distribution is not Gaussian, and therefore, T MB cannot be determined.
T kin /T bath T kin,rel /T bath T MB /T bath A (Fig. 8a-c

FIG. 1 .
FIG. 1. Kinetic temperature and area fraction at coexistence.From left to right we show a snapshot of the particle positions in the steady state (a, e, and i), the local area fraction of active (b, f, and j) and passive (c, g, and k) particles, and the coarse-grained kinetic temperature field of the passive tracer particles (d, h, and l) averaged over time in the steady state, respectively.Panels i-l are slightly zoomed in and the black and white solid lines are guides to the eye denoting the border of the dense phase.Parameters:

FIG. 3 .
FIG. 3. Velocity distribution of passive tracers.Distribution of the x component v x of the passive particles' velocities a in the uniform state at x a = 0.2 and b,c in the phase-separated state at x a = 0.8 for passive particles in the dense and the gas phase, respectively, averaged over time in the steady state.Different Pe values are given in the key (parameters as in Fig. 2).The black dashed lines are Gaussian fits.d Variance of v x as function of Pe for the three cases shown in panels a-c averaged over 5 independent ensembles in the steady state.

FIG. 4 .
FIG. 4. Correlated particle dynamics.a,e Distribution of v • F int for the passive tracers (PBPs) averaged over time in the steady state.The dashed vertical lines mark the corresponding mean values.b Exemplary trajectories of an active and a passive particle in the dilute phase (gray dashed and blue solid line, respectively), where the active particle pushes the passive particle forward such that v and F int are aligned and v • F int is large.f Exemplary trajectory of a passive particle in a liquid-like droplet pushed forward as a result of correlated dynamics of the active particles.c,g Snapshots of the corresponding simulations in the steady state.d,h Corresponding coarse-grained values of v • F int .The black and white solid lines in panels c,d,g, and h indicate the border of the dense phase.Parameters: Pe = 80 and x a = 0.7 (a,c, and d), Pe = 100 and x a = 0.6 (b), Pe = 400 and x a = 0.9 (e-h); other parameters as in Fig. 2.

FIG. 5 .
FIG. 5. Schematic illustration of the mechanism.a At low Pe, particles are trapped in the dense phase and have no space to speed up.Long correlated trajectories where active particles push passive ones forward occur only in the dilute phase.Thus, passive particles are faster in the dilute phase.b At large Pe (fast self-propulsion), long correlated trajectories where active particles push passive ones forward occur even in the dense phase (and are supported by collective motion of the active particles).The more frequent collisions in combination with the collective motion of the ABPs in the dense phase lead to faster passive particles in liquid-like droplets compared to the surrounding gas.

FIG. 6 .
FIG. 6. Non-equilibrium state diagram for 2 × 10 4 particles.a (I) uniform state, (II a ) hot gas and cold liquid coexistence, and (II b ) cold gas and hot liquid-like droplets coexistence.The colors denote the (normalized) kinetic temperature difference of the passive particles in the dense and the gas phase.The black dashed line denotes the transition line as obtained in Ref. [58] for a purely overdamped mixture.Panel b shows a zoomed version of panel a and the red circle denotes the maximum temperature difference in regime II a .All data is averaged over a time interval ∆t = 900τ p in the steady state (parameters as in Fig.

FIG. 7 .
FIG. 7. Effect of the particle size.a-d Simulation snapshots with four different sizes of the passive particles as given in the panel titles.The passive particles are colored with respect to their kinetic temperature.e,f Velocity distribution of the passive particles in the dense and the gas phase, respectively, for different size ratios as indicated in the key.g Mean-square displacement of the active (dashed lines) and passive (solid lines) particles.h Corresponding kinetic temperatures of passive particles in the dense and the gas phase.Parameters: Pe = 100, γt = σ p γ t /σ a , φ tot = 0.7, N a + N p = 2 × 10 4 for σ p /σ a = 1, 2, 3, 4

FIG. 8 .
FIG. 8. Spatial velocity correlations.a-c Snapshot in the steady state, spatial velocity correlation function of passive particles [Eq.(6)], and velocity distribution of passive particles in the dense and the gas phase, respectively, for a simulation of N = 20 000 particles at x a = 0.90.The dashed and dotted lines in panel b denote the mean distance d between passive particles obtained from a Voronoi tessellation in the dense phase ( d ≈ 4.6σ) and in the gas phase ( d ≈ 5.8σ), respectively.d-f Same as a-c but for a simulation of N = 100 000 particles at x a = 0.996.The mean distance between the passive particles denoted by the dashed and dotted line in panel e is d ≈ 41σ (liquid) and d ≈ 21σ (gas), respectively.All data has been averaged over time in the steady state.Simulation parameters: Pe = 400 (other parameters as in Fig. 2).

6 +
ϵ, r nl /σ ≤ 2 1/6 0, else with particle diameter σ, strength ϵ, and r nl = |r n − r l |.For the simulations with active and passive particles of different diameters, the effective diameter for the interaction between the active and passive particles is chosen as σ ap = (σ a + σ p )/2, where σ a and σ p denote the diameters of the active and passive particles, respectively.In all simulations, we fix m a /(γ t τ p ) = 5 × 10 −5 , I/(γ r τ p ) = 5 × 10 −6 to recover overdamped dynamics for the active particles[11].For the passive particles, we fix m p /(γ t τ p ) = 5 × 10 −2 with the persistence time τ p = 1/D r .Furthermore, we set ϵ = 10k B T b , γt = γ t , and σ a = σ p = σ = D t /D r (unless otherwise indicated), and we use systems with N = N a + N p particles.We choose γ t = γ r /σ 2 and vary Pe and the fraction x a = N a /(N a + N p ) of the active particles.The total area fraction φ tot = (N a + N p )πσ 2 /(4A) is set to φ tot = 0.5, where A denotes the area of the simulation box.The Langevin equations [Eqs.

FIG. S1 .
FIG. S1.Velocity distribution.Distribution of the x component of the velocities of the passive particles a in the uniform state at Pe = 10 and in the MIPS state at Pe = 100 b in the dense phase and c in the gas phase for different values of x a as given in the key (other parameters as in Fig. 2 in the main text).The black dashed lines are Gaussian fits showing that the distributions are clearly non-Gaussian in the MIPS state.