This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Work and entropy of mixing in isolated quantum systems

Budhaditya Bhattacharjee budhadityab@ibs.re.kr Center for Theoretical Physics of Complex Systems, Institute for Basic Science (IBS), Daejeon - 34126, Korea    Rohit Kishan Ray rkray@vt.edu Center for Theoretical Physics of Complex Systems, Institute for Basic Science (IBS), Daejeon - 34126, Korea Department of Material Science and Engineering, Virginia Tech, Blacksburg, VA 24061, USA    Dominik Šafránek dsafranekibs@gmail.com Center for Theoretical Physics of Complex Systems, Institute for Basic Science (IBS), Daejeon - 34126, Korea
(September 29, 2025)
Abstract

The mixing of two different gases is one of the most common natural phenomena, with applications ranging from CO2\text{CO}_{2} capture to water purification. Traditionally, mixing is analyzed in the context of local thermal equilibrium, where systems exchange energy with a heat bath. Here, we study mixing in an isolated system with potentially non-equilibrium initial states, characterized solely by macroscopic observables. We identify the entropy of mixing as a special case of observational entropy within an observer-dependent framework, where both entropy and extractable work depend on the resolution of measurement. This approach naturally resolves the Gibbs mixing paradox in quantum systems: while an observer experiences a discontinuous increase in entropy upon learning of the existence of two particle types, this knowledge does not provide an advantage in work extraction if the types of particles remain operationally indistinguishable in their measurements. Finally, we derive a Landauer-like bound on the difference in energy extracted by two observers, where an ‘observational temperature’ emerges, determined by the accessible information. These results provide a foundation for rigorously determining the energy required to unmix in non-equilibrium settings and extend beyond quantum systems, offering insights into the thermodynamics of isolated classical gases.

I Introduction

Mixing two types of particles is one of the most common occurrences in nature. Whether it is the down-to-earth act of putting milk into our coffee or the release of CO2\text{CO}_{2} into the atmosphere, the underlying physical laws remain the same. Such processes also have practical implications—most notably in determining how much energy is needed to separate two compounds once mixed. For gases, this minimal energy has a fundamental bound, revealing that capturing CO2\text{CO}_{2} before it is released is far more efficient than attempting to remove it afterward [1]. A similar principle applies to other modern environmental challenges, such as desalination [2] or filtering microplastics [3] and hormones [4] from our water supply, where proactive measures can significantly reduce both costs and ecological impacts.

In classical physics, where gases are assumed to be in thermal equilibrium at a well-defined temperature, the energy cost of mixing is well understood. However, in systems that do not equilibrate or are not connected to a thermal bath the usual assumptions break down. Consequently, it becomes unclear how much energy is required to separate the gases. This challenge is especially pronounced in isolated quantum systems, which are now routinely realized and studied in experimental settings [5, 6, 7].

These questions are also fundamentally observer-dependent: to study how much energy it costs to separate the two gases, one must first define what it means “to separate.” Typically, this notion carries a connotation of spatial separation, such as one type of gas being confined within a box while the other remains outside. Even in this case, however, the energy cost depends on the size of the box, which is determined by the observer or the agent attempting to separate the gases [8]. This definition is inherently tied to what the agent can achieve, shaped by their specific capabilities and limitations in carrying out the separation [9].

This interplay between the observer’s experimental capabilities and knowledge lies at the heart of the Gibbs mixing paradox (GMP). Consider two observers, one who can distinguish the two gases and the other who cannot. Recognizing the difference, the first observer assigns a higher entropy to the observed state, based on the larger number of microstates in compliance with any macroscopic observation. The second observer perceives no distinction and assigns a lower entropy, asserting that all particles are of the same type. Crucially, the first observer’s knowledge provides an advantage—if the gases are spatially separated, they can extract energy by allowing the gases to mix (classically, through semi-permeable membranes). The second observer, unaware of the distinction, cannot exploit this opportunity. As E.T. Jaynes noted [9]: “We would observe, however, that the number of fish that you can catch is an ‘objective experimental fact’; yet it depends on how much ‘subjective’ information you have about the behavior of fish.”

While the role of different observers is widely acknowledged in these discussions, their treatment remains somewhat ambiguous in the absence of a formal framework defining what constitutes an “observer.” Instead, the situation is typically described qualitatively, with explanations stating that one observer can distinguish the two gases while the other cannot. From this premise, the appropriate definition of equilibrium entropy is either postulated or argued for, but it is not systematically deduced from the original statement. This lack of rigorous derivation leaves the connection between the observer’s capabilities and the entropy definition somewhat imprecise. Does “being able to distinguish” imply the ability to probe each particle individually and determine its type, or does it refer to probing a specific region of space to ascertain the proportion of particles of each type within it? These two interpretations are fundamentally different and lead to distinct outcomes, both in terms of the observer’s knowledge and the resulting thermodynamic description of the system.

This also highlights another aspect of the paradox: the abrupt, discontinuous change in entropy that occurs when particles are distinguishable, even if differences between them are minimal. The mere potential for distinction leads to a significantly larger associated entropy. This apparent paradox arises from treating entropy as a fundamental property of the physical system. However, if entropy is instead understood as an informational measure—reflecting the observer’s perceived ignorance about the system and defining the limits of their ability to extract energy—the paradox dissolves. In this interpretation, entropy is relational, shaped by the observer’s knowledge and capabilities rather than an intrinsic characteristic of the system.

In this paper, we formalize the concept of an observer and their capabilities in precise mathematical terms. This systematic approach removes the need for overly intuitive reasoning. Instead, it allows us to directly approach the question with higher clarity and precision. Thereby it allows us to quantify the additional energy that an observer with a higher resolution can extract from a system compared to an observer with a lower resolution. As an application, we analyze the difference in extracted energy between two observers in the Gibbs mixing scenario, in the generic 1D lattice model with fermions of two colors. Additionally, we focus on isolated quantum systems and compute the difference in energy that an observer can extract before and after the two gases mix, thereby establishing a lower bound on the energy required to separate them again.

Interestingly, this result bears a strong mathematical resemblance to Landauer’s bound, with a key difference: instead of the bath temperature, an “observational temperature” emerges, defined by the observer’s experimental capabilities and limitations. Furthermore, since equivalent reasoning applies to classical systems, the insights gained extend beyond the quantum case. The framework allows us to discuss the energy required to separate isolated classical gases out of equilibrium, broadening the applicability of the results.

First, we review the background including mixing and observational entropy (Sec. II). Then, we derive the difference between extracted energy by two observers (Sec. III), with a specific result resembling Landauer’s bound for observers that have similar capabilities (Sec. III.1), followed by the difference in extracted energy by the same observer at different times (Sec. III.2). In Sec. IV, we analyze entropy in the Gibbs mixing scenario in quantum systems for an observer that can distinguish the particles (Sec. IV.1) and three versions of an observer who cannot (Sec. IV.2IV.3IV.4), and compare them (Sec. IV.5). We also analyze the change in entropy for these observers when an isolated quantum system evolves IV.6, showing that for all observers with various degrees of ignorance, the observed change in entropy is the same. Then we combine the previous results to show the difference between extracted energy in the Gibbs mixing scenario (Sec. V). In particular, we outline the assumptions under which we prove that an observer unaware of the two particle types extracts the same amount of energy as one who is aware, thereby resolving the Gibbs mixing paradox. We also include two simulations: first, the difference between two observers, and second, the difference for the same observer but at different times. Finally, we conclude (Sec. VI).

Refer to caption
Figure 1: A schematic representation of the Gibbs mixing paradox. Rick can distinguish the particles in the two partitions based on their colors, whereas Morty cannot. As a result, when the partition is removed, Rick observes a net increase in entropy due to the mixing of two distinct gases, while Morty, unable to perceive the distinction, does not observe the same change in entropy.

II Background

Consider two ideal gases—initially separated by a partition, both with identical volume, pressure, and temperature—that are allowed to mix. The Gibbs Mixing Paradox (GMP) arises because of the differing observer abilities. An Observer Rick, who can distinguish the gases notices an entropy increase from mixing. Morty, on the other hand, cannot distinguish them and thus attributes no entropy change. This discontinuity depends on the gases being distinguishable, challenging classical thermodynamic notions of entropy (see Fig. 1). Note here, that we consider distinguishability between the gases rather than their constituent individual particles, which may be fundamentally indistinguishable in quantum physics. In second quantization, this corresponds to being able to distinguish different spatial modes or different values of spin.

Gibbs Mixing Paradox. Jaynes [9] framed entropy as a subjective measure tied to an observer’s coarse-graining of microstates. For Rick, reverting the system to its original state requires work proportional to the entropy change; for Morty, simply reinserting the partition suffices. This perspective anchors entropy in information. Building on Jaynes’ perspective, many authors have explored the GMP from both philosophical and physical standpoints [10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26]. These approaches can be further categorized. Philosophical approaches focus on particle identity and the foundation of entropy. For example, Maslov [14] emphasizes contextual particle distinguishability, while Ihnatovych [24] critiques classical entropy additivity. From an information-theoretic perspective, Lairez [26] aligns with Jaynes. Operational approaches resolve the discontinuity pragmatically. Allahverdyan and Nieuwenhuizen [11] replace coarse-graining with ergotropy (extractable work), showing it decreases continuously as gases become indistinguishable. Peters [19] attributes entropy changes to uncertainty in particle locations. Also, in agreement with Jaynes, Yadin et al. [27] reconcile observer-dependent work extraction with respect to a single heat bath in a first-quantization description. With a setup similar to the Hong-Ou-Mandel experiment, they argue that Morty can also extract work, contrary to the predictions of the classical GMP.

Mixing Beyond GMP. Mixing has also been studied independently. Yoshida and Nakagawa [28] quantify free energy changes in small quantum systems, while Takakura [29] explore entropy in generalized probabilistic theories. Sasa et al. [30] analyze quasistatic work in quantum thermal processes. These works highlight broader thermodynamic implications of mixing beyond classical paradoxes.

Observational entropy. Von Neumann argued that von Neumann entropy does not accurately represent thermodynamic entropy, as it assigns zero entropy to all pure states and fails to capture the information loss that arises when only macroscopic degrees of freedom are observed. As an alternative, he introduced macroscopic entropy [31, 32], which uses state-counting similar to Boltzmann entropy. This concept has recently gained renewed interest and been generalized to observational entropy [33, 34, 35, 36, 37, 38, 39], leading to a range of applications [40, 41, 42, 43, 44, 45, 46, 47, 48].

For a projective measurement defining a complete set of orthogonal projectors, known as coarse-graining 𝒞={Pi}{\mathcal{C}}=\{{P}_{i}\}, observational entropy is defined as

S𝒞=ipilnpi+ipilnVi.S_{\mathcal{C}}=-\sum_{i}p_{i}\ln p_{i}+\sum_{i}p_{i}\ln V_{i}. (1)

Here, pi=tr[Piρ]p_{i}=\tr[{P}_{i}{\rho}] represents the probability of the system being in macrostate ii, while Vi=tr[Pi]V_{i}=\tr[{P}_{i}] denotes the volume of the macrostate, i.e., the number of microstates (pure states) contained within it. Each projector also corresponds to a subspace, defined by its support. Thus, we can equivalently write that a coarse-graining corresponds to a partition of the Hilbert space into subspaces, =ii{\mathcal{H}}=\bigoplus_{i}{\mathcal{H}}_{i}, where Vi=dimiV_{i}=\dim{\mathcal{H}}_{i}. The first term in observational entropy corresponds to the Shannon entropy, capturing the uncertainty in identifying the macrostate to which the system belongs. The second term represents the mean Boltzmann entropy, quantifying the residual uncertainty in determining the specific microstate due to the coarse-grained nature of the measurement.

quantum

  Refer to caption

classical

  Refer to caption
Figure 2: Sketch of the extraction protocol in quantum (top panels) and classical (bottom panels) isolated systems. First, a coarse measurement is performed, which provides probability of finding the system in a given macrostate. In the quantum case, we sketch the Hilbert space with the measurement basis given by a set of projectors. In the example, we have chosen the case where the observer cannot distinguish the two excited state from each other but can distinguish them from the ground state, corresponding to coarse-graining 𝒞={|gg|,|e1e1|+|e2e2|}{\mathcal{C}}=\{\outerproduct{g}{g},\outerproduct{e_{1}}{e_{1}}+\outerproduct{e_{2}}{e_{2}}\}. The first projector projects onto the x axis and the second to the y-z plane. We start from an initial state |ψ=13|e1+23|e2\ket{\psi}=\frac{1}{\sqrt{3}}\ket{e_{1}}+\sqrt{\frac{2}{3}}\ket{e_{2}} which is unknown to the observer, represented by a vector lying in y-z plane. First, a random unitary U~\tilde{U} is applied, which effectively spreads this state over a macrostate (i.e., on average, as described by operation 𝒰~=U~U~𝑑μ(U~)\mathcal{\tilde{U}}=\int\tilde{U}\bullet\tilde{U}^{\dagger}d\mu(\tilde{U})). The resulting state ρ1=12(|e1e1|+|e2e2|){\rho}_{1}=\frac{1}{2}(\outerproduct{e_{1}}{e_{1}}+\outerproduct{e_{2}}{e_{2}}) is known to the observer, represented by a circle in the y-z plane. Finally, an extraction unitary UU is applied, which lowers the energy of the state, extracting it in the process, yielding ρ2=12(|gg|+|e1e1|){\rho}_{2}=\frac{1}{2}(\outerproduct{g}{g}+\outerproduct{e_{1}}{e_{1}}) (represented by a circle lying in the x-y plane). In the classical case, we sketch the phase space with the coarse-graining given by dividing it into cells (dashed grid). Starting from some unknown initial state, given by a point in phase space (blue dot), a random canonical transformation is applied, spreading the state over the macrostate (blue box). Then, a canonical extraction transformation is applied which reduces its energy (blue elliptic disk). This represents the example of a classical harmonic oscillator for which the energy shells form ellipses, with the smallest energy concentrated in the middle.

III Difference in the work extracted between two observers

In this section, we derive the difference in work extracted by two observers with different observational capabilities.

Consider two observers, Rick and Morty, with two measuring devices of different precision. These measuring devices correspond to coarse-grainings 𝒞R{\mathcal{C}_{\mathrm{R}}} and 𝒞M{\mathcal{C}_{\mathrm{M}}}, We define the difference in observational entropies as

ΔSS𝒞MS𝒞R,\Delta S\equiv S_{\mathcal{C}_{\mathrm{M}}}-S_{{\mathcal{C}_{\mathrm{R}}}}, (2)

We assume that the measurement of Rick has a higher precision as measured by lower observational entropy, meaning that ΔS0\Delta S\geq 0. For what follows, we require this assumption only for the particular quantum state from which the observers will extract energy, while for other states this inequality may be reversed. This will make all other quantities of type Δ\Delta introduced below positive.

The work extracted by performing a joint extraction on a large number of copies of the system, which is characterized by the measurement probabilities in a single coarse-grained basis 𝒞{\mathcal{C}}, is given by observational ergotropy [49],

W𝒞=tr[H(ρρβ)].W_{\mathcal{C}}^{\infty}=\tr[{H}({\rho}-{\rho}_{\beta})]. (3)

There, ρβ=exp(βH)/Z{\rho}_{\beta}=\exp(-\beta{H})/Z is the thermal state with inverse temperature β=1/kT\beta=1/kT implicitly defined by the von Neumann entropy, as a solution to SvN(ρβ)=S𝒞S_{{\mathrm{vN}}}({\rho}_{\beta})=S_{\mathcal{C}}, and kk is the Boltzmann constant. We call temperature TT obtained by solving this equation an observational temperature. Expression tr[Hρ]\tr[{H}{\rho}] represents the initial energy of the state from which the energy is extracted, while tr[Hρβ]\tr[{H}{\rho}_{\beta}] is the average energy of the final state after the extraction protocol is employed.

Type of system Heat capacity Work difference in units of kTΔSkT\Delta S
Ideal gas 32Nk\frac{3}{2}Nk 113ΔSN+227(ΔSN)21-\frac{1}{3}\frac{\Delta S}{N}+\frac{2}{27}\left(\frac{\Delta S}{N}\right)^{2}
Debye model at low temperatures 324Nk(TTD)3324Nk\left(\frac{T}{T_{D}}\right)^{3} 112ΔS324NTD3T313(ΔS324NTD3T3)21-\frac{1}{2}\frac{\Delta S}{324N}\frac{T_{D}^{3}}{T^{3}}-\frac{1}{3}\left(\frac{\Delta S}{324N}\frac{T_{D}^{3}}{T^{3}}\right)^{2}
Quantum Critical Metals π6mNT(1+ANT1/3)\frac{\pi}{6}mNT\left(1+\frac{A}{NT^{1/3}}\right) 13kπmΔSNT(1+ANT1/3)+2k2π2m2(1+NT1/3A)(ΔSNT(1+ANT1/3))21-\frac{3k}{\pi m}\frac{\Delta S}{NT\left(1+\frac{A}{NT^{1/3}}\right)}+\frac{2k^{2}}{\pi^{2}m^{2}\left(1+\frac{NT^{1/3}}{A}\right)}\left(\frac{\Delta S}{NT\left(1+\frac{A}{NT^{1/3}}\right)}\right)^{2}
Liquid Helium above the lambda point Nk(A+Btα)Nk(A+Bt^{-\alpha}) 112ΔSN(A+Btα)+16(1+(1+t)α(1+ABtα)t)(ΔSN(A+Btα))21-\frac{1}{2}\frac{\Delta S}{N(A+Bt^{-\alpha})}+\frac{1}{6}\left(1+\frac{(1+t)\alpha}{(1+\frac{A}{B}t^{\alpha})t}\right)\left(\frac{\Delta S}{N(A+Bt^{-\alpha})}\right)^{2}
s-wave superconductor at low temperatures ANkTeΔ/kTANk\sqrt{T}e^{-\Delta/kT} 112ΔSeΔ/kTANT+112(12ΔkT)(ΔSeΔ/kTANT)21-\frac{1}{2}\frac{\Delta Se^{\Delta/kT}}{AN\sqrt{T}}+\frac{1}{12}\left(1-\frac{2\Delta}{kT}\right)\left(\frac{\Delta Se^{\Delta/kT}}{AN\sqrt{T}}\right)^{2}
Table 1: Work difference as per Eq. 7 in units of kTΔSkT\Delta S for a variety of classical and quantum gases. AA, BB, Δ\Delta, mm and TDT_{D} are constants that depend on the specific model and TT denotes observational temperature. For the liquid helium we defined t=T/Tc1t=T/T_{c}-1 with T>TcT>T_{c}. For liquid Helium, the second factor diverges as the observational temperature approaches the critical temperature. For others, the expansion diverges as temperature converges to absolute zero. As a result, the expansion holds only for very small entropy differences around these critical points. Additionally, this indicates that there is a critical amount of information, given by S𝒞=SvN(ρβc)S_{\mathcal{C}}=S_{{\mathrm{vN}}}({\rho}_{\beta_{c}}), which the observer is required to possess in order to extract significantly more energy.

The work extraction protocol is partially random, meaning that one first employs a randomizing unitary which washes away the information about the degrees of freedom that cannot be obtained from the coarse measurement 𝒞{\mathcal{C}}, and then the global work extraction unitary that is optimized for the average state. In the language of quantum reference frames, the first step is known as “twirling” and encodes the lack of a quantum reference frame [50]. Observational ergotropy, Eq. (3), then represents the average extracted work upon repetition of this protocol many times. See Fig. 2 for the sketch of this extraction protocol. For more details see Refs. [49, 51, 44].

The difference in extracted energies by the two observers is equal to,

ΔWW𝒞RW𝒞M=tr[H(ρβMρβR)],\Delta W\equiv W_{{\mathcal{C}_{\mathrm{R}}}}^{\infty}-W_{{\mathcal{C}_{\mathrm{M}}}}^{\infty}=\tr[{H}({\rho}_{{\beta_{\mathrm{M}}}}-{\rho}_{{\beta_{\mathrm{R}}}})], (4)

where SvN(ρβR)=S𝒞RS_{{\mathrm{vN}}}({\rho}_{{\beta_{\mathrm{R}}}})=S_{\mathcal{C}_{\mathrm{R}}} and SvN(ρβM)=S𝒞MS_{{\mathrm{vN}}}({\rho}_{{\beta_{\mathrm{M}}}})=S_{{\mathcal{C}_{\mathrm{M}}}}.

This formula holds completely generally, for any two coarse-grainings. The positive difference in entropies ΔS\Delta S implies that work extracted ΔW\Delta W is also positive.

The formula for the observational ergotropy fundamentally depends on the initial energy, which is unknowable to the observers given that they obtain only limited information by measuring in their basis, 𝒞R{\mathcal{C}_{\mathrm{R}}} and 𝒞M{\mathcal{C}_{\mathrm{M}}}, respectively. It can, however, be estimated from these measurements [52]. On the other hand, the formula for the difference no longer depends on this initial energy, only on the respective observational entropies, which allows us to study it in more detail.

Note that in some scenarios, especially when the initial state energy is small, execution of the protocol may cost more energy than is being extracted for both observers. The positive difference in extracted work ΔW\Delta W then means that Rick has lost less energy than Morty by executing the protocol.

III.1 Different observers with similar capabilities

We show that when the difference in observational entropies is small, the extracted work difference is given by an expression that strikingly resembles Landauer’s principle.

When ΔS\Delta S is small, then the difference between corresponding temperatures are small, too. We define inverse observational temperatures of Rick and Morty implicitly via SvN(ρβR)=S𝒞RS_{{\mathrm{vN}}}({\rho}_{{\beta_{\mathrm{R}}}})=S_{\mathcal{C}_{\mathrm{R}}} and SvN(ρβM)=S𝒞MS_{{\mathrm{vN}}}({\rho}_{{\beta_{\mathrm{M}}}})=S_{{\mathcal{C}_{\mathrm{M}}}}, respectively. The temperature difference

ΔββRβM,\Delta\beta\equiv{\beta_{\mathrm{R}}}-{\beta_{\mathrm{M}}}, (5)

is positive for positive ΔS\Delta S, implying TRTMTT_{R}\leq T_{M}\equiv T (meaning that TT is the observational temperature related to the higher of the two observational entropies). Further, the temperature difference vanishes as ΔS\Delta S vanishes.

We derive the work difference as an expansion in small ΔS\Delta S 111In the example of the mixing paradox below, in the initial situation of the macroscopic non-equilibirum state with equal-sized boxes and non-dilute scenario, ΔS\Delta S is about a third of SS, so this condition is not satisfied. Including higher order terms will yield better agreement with the exact formula, Eq. (4), and we employ empirical approaches to do that. For the unequal sized boxes, or sufficiently dilute gas, this condition will generally be satisfied.. To do that, we first insert the identities for observational entropy below Eq. (4) into Eq. (2) and expand it in the powers of Δβ\Delta\beta up to the third order coefficient. Then, we solve this equation for Δβ\Delta\beta while ignoring the higher order terms, yielding an expansion Δβ(ΔS)\Delta\beta(\Delta S). Second, we expand the work difference, Eq. (4), in terms of Δβ\Delta\beta.

Inserting the first expansion into the second, we obtain an expression for the ergotropy difference,

ΔW=kTΔS(1(kT)22vΔS+(3kT+X)(kT)36v2ΔS2+).\Delta W\!=\!kT\Delta S\!\left(1\!-\!\frac{(kT)^{2}}{2v}\Delta S\!+\!\frac{(3kT\!+\!X)(kT)^{3}}{6v^{2}}\Delta S^{2}\!+{\!\cdot\!\cdot\!\cdot\!}\right)\!. (6)

Here, v=dEdβ=E2E2v=-\frac{d\langle E\rangle}{d\beta}=\langle E^{2}\rangle-\langle E\rangle^{2} is the variance in energy of the thermal state corresponding to temperature β\beta, and X=dlnvdβ=(E3+3EE22E3)/vX=\frac{d\ln v}{d\beta}=(-\langle E^{3}\rangle+3\langle E\rangle\langle E^{2}\rangle-2\langle E\rangle^{3})/v, a term which encodes the third moment in energy.

Realizing that the variance is related to the heat capacity CE=E/T{C_{\!E}}=\partial E/\partial T through v=dEdβ=CEkT2v=-\frac{d\langle E\rangle}{d\beta}={C_{\!E}}kT^{2}, we obtain a more physically-relevant expression for the ergotropy difference,

ΔW=kTΔS(112ΔSC~E+16(1lnC~ETT)(ΔSC~E)2+).\Delta W\!\!=\!kT\Delta S\!\left(\!1\!-\!\frac{1}{2}\frac{\Delta S}{{\tilde{C}_{\!E}}}\!+\!\frac{1}{6}\!\left(1\!-\!\frac{\partial\ln{\tilde{C}_{\!E}}}{\partial T}T\right)\!\left(\frac{\Delta S}{{\tilde{C}_{\!E}}}\right)^{\!2}\!\!\!+{\!\cdot\!\cdot\!\cdot\!}\right)\!. (7)

This represents our main result. C~E=CE/k{\tilde{C}_{\!E}}={C_{\!E}}/k is the dimensionless heat capacity, which measures the total number of degrees of freedom that contribute to the heat capacity.

The zeroth order term strikingly resembles Landauer’s principle about the minimal required energy to erase a bit of data. There are a few notable differences:

First, this formula represents the maximal extractable work given the information at hand, which is quantified by observational entropy. The work is extracted by applying the extraction unitary, which is dictated by the measurement outcome. This is similar to the Szilard engine. Also there, the position of the particle is determined by the measurement. Based on the outcome, the weight is attached to the piston at the correct side and is lifted. The maximal extracted work for the Szilard engine has been extensively studied 222See reviews [73, 74, 75] and papers for quantum [76, 77, 78, 79, 80, 81, 82] (Maxwell’s Demon, Szilard Engine and Landauer Principle), classical [83, 84, 85], and active matter [86, 87, 88] systems..

The difference to the Szilard engine is that this result applies to isolated systems, while the Szilard engine operates at a fixed temperature of the surrounding thermal bath TT. Here, in contrast, the temperature TT is not a physical temperature of a thermal bath, but rather the “informational” temperature corresponding to the observational entropy. To be more precise, we identified the temperature TT as the higher of the two observational temperatures, see Eq. (5). Given that we assume the entropy difference ΔS\Delta S to be small, any difference between the observational temperatures remains small as well.

Physically, the observational temperature is the temperature of the resulting final thermal state when the randomized extraction protocol is employed. As already mentioned above, if the mean energy of such a thermal state is higher than the initial mean energy, employing the protocol results in a net energy investment. If it is lower, energy is extracted. A positive ergotropy difference then means that Rick has gained ΔW\Delta W amount more energy (or lost less energy) than Morty when employing the protocol.

The first correction term is always negative, meaning that kTΔSkT\Delta S bounds the work difference ΔW\Delta W from above. This term also implies that the work difference is the largest when (kT)2v(kT)^{2}\ll v, meaning that energy fluctuations are much larger than the thermal energy scale. Equivalently, these are the situations when the heat capacity grows very large or diverges. This happens, for example, for systems near a critical point or in a low-temperature regime with large degeneracy [55, 56, 57, 58, 59].

The second-order correction can be both positive or negative, and it will generally depend on how the heat capacity scales with temperature. Considering power law C~ETα{\tilde{C}_{\!E}}\propto T^{\alpha}, the correction is negative for a superlinear scaling α>1\alpha>1, meaning that the number of degrees of freedom contributing to the internal energy grows faster than linearly with temperature. It is positive otherwise. A few examples of classical and quantum gases with varying heat capacities are illustrated in Table 1.

Finally, note that there was not much particularly quantum about this calculation, apart from the unitary extraction when employing the protocol. Observational entropy can be consistently defined also for classical systems [60]. Furthermore, classical equivalents of unitary operations are symplectic (canonical) transformations, which preserve all phase-space volume, the symplectic product between the two pase-space vectors, and the relationship between conjugate variables (the Poisson brackets), sketched in Fig. 2. One then infers that this result also applies identically to classical isolated systems of indistinguishable particles.

III.2 Work extracted at different times

Equivalent calculations can also be performed for the same observer who extracts work at different times. Consider an initial state ρ0{\rho}_{0}, describing, for example, particles not being mixed yet. After some time, the system evolves into state ρ1{\rho}_{1}, e.g., a state where particles have mixed already. The ergotropy difference can, in that case, also describe how much the same observer can extract at different times. Defining two observational entropies and the corresponding thermal state SvN(ρβi)=S𝒞(ρi)S_{{\mathrm{vN}}}({\rho}_{\beta_{i}})=S_{\mathcal{C}}({\rho}_{i}), i=0,1i=0,1, the ergotropy difference is given by ΔW=tr[H(ρβ1ρβ0)]\Delta W=\tr[{H}({\rho}_{\beta_{1}}-{\rho}_{\beta_{0}})]. The same approximate expression, Eq. (7), also holds.

This indicates that an observer at an earlier time, when the particles are separated, can extract ΔW\Delta W more energy than at a later time when all particles are mixed. Then, one could argue, it must cost at least as much energy to separate the two gases, meaning that the ergotropy difference also gives the lower bound on the amount of work the observer needs to perform to macroscopically separate these gases. This would hold generally, not only for mixing gases: to undo the macroscopic changes observed by an observer with an ability to do a macroscopic measurement 𝒞{\mathcal{C}}, they must spend at least ΔW\Delta W amount of energy to do so.

While this is intuitively reasonable, there is a caveat: the difference in extracted work assumes the same protocol is applied at both times. However, a different protocol might reverse the macroscopic changes at a lower cost. In principle, such a protocol must exist—if the full quantum state were known, the reversal could be performed without extra energy. Morally, we would argue that with limited information, no protocol should achieve this reversal at an average cost lower than ΔW\Delta W. A rigorous proof is left for future work.

IV Observational entropy of quantum mixing

We derive the entropy difference ΔS\Delta S, which we can then input into the formula for the ergotropy difference, for a quantum version of the Gibbs mixing paradox.

We consider observers illustrated in Fig. 3: Rick, who is able to determine the total number of particles of each color in both the left and the right side of the box; and three clones of Morty, who can identify the total number of particles in each part of the box but not their color, but have varying degrees of knowledge about colors of particles in the full system. Morty 1: knows the total number of blue particles in the full system; Morty 2: does not know the total number of particles in the full system, but is aware that these two types exist; and Morty 3: is not even aware of the existence of the two types of colors and as a result, will associate incorrect volumes to the observed macrostates.

The ability to determine or identify is described by a specific coarse-graining representing a macroscopic measurement that goes into the definition of observational entropy that either of them associates to the observed system.

In the following calculations, we consider both the blue and the red particles to be either fermions or hard-core bosons, which will affect the counting, however, a very similar derivation can also be done with bosons. The results for these will coincide when the gas is sufficiently dilute.

We consider the total Hilbert space to be that of any possible configuration of any number of blue and red particles on the lattice of LL sites – i.e., spatial modes in second quantization. Here, color could be a spin degree of freedom (immutable within the observers capability; evolved with a spin-preserving Hamiltonian). The actual size of the total Hilbert space will not actually matter, as long as every state that we consider can be contained within it. For instance, when considering Hamiltonians that preserve the total number of particles of each type, we can restrict to a Hilbert space that respects this symmetry. However, both of these approaches will give the same results since these come from counting the dimensions of the macrostates and not the dimension of the full Hilbert space.

For simplicity, we assume that both sides are of the same lattice size, ww. The non-symmetric case, of LAL_{A} and LBL_{B} being the size of the left and the right side, respectively, can be obtained by substituting lnwalnLA+(1a)lnLB\ln w\rightarrow a\ln L_{A}+(1-a)\ln L_{B} into the results below (aa is defined below).

Refer to caption
Figure 3: The schematic for counting the states for Rick. In this illustration, the number of sites on each side is ω=5\omega=5. Following the notation in the text, we can see that the number of blue particles is N+=3N_{+}=3, the number of red particles is N=3N_{-}=3, the number of particles on the left side is NA=3N_{A}=3, and the number of particles on the right NB=3N_{B}=3. Colors mix as the system evolves and Rick observes an increase in entropy.

IV.1 Rick

Rick can observe the total number of particles of each color on both the left and the right side of the box; see Fig. 3 (top panel). The non-symmetric case is derived in Appendix B.5.

Rick’s macrostate is given by four numbers, the total number of blue and red particles, on the left and on the right side. This macrostate is mathematically described by a subspace,

R=+AA+BB.{\mathcal{H}}^{\mathrm{R}}=\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.16796875,0.3125,0.640625}\definecolor[named]{pgfstrokecolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@color@rgb@stroke{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@moveto{-8.39357pt}{-7.39755pt}\pgfsys@curveto{-8.39357pt}{-6.14043pt}{-9.41266pt}{-5.12134pt}{-10.66978pt}{-5.12134pt}\pgfsys@curveto{-11.92691pt}{-5.12134pt}{-12.946pt}{-6.14043pt}{-12.946pt}{-7.39755pt}\pgfsys@curveto{-12.946pt}{-8.65468pt}{-11.92691pt}{-9.67377pt}{-10.66978pt}{-9.67377pt}\pgfsys@curveto{-9.41266pt}{-9.67377pt}{-8.39357pt}{-8.65468pt}{-8.39357pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{-14.22638pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{{\mathcal{H}}_{+}^{A}}\otimes\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.9375,0.23046875,0.15234375}\definecolor[named]{pgfstrokecolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@color@rgb@stroke{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@moveto{-8.39357pt}{-7.39755pt}\pgfsys@curveto{-8.39357pt}{-6.14043pt}{-9.41266pt}{-5.12134pt}{-10.66978pt}{-5.12134pt}\pgfsys@curveto{-11.92691pt}{-5.12134pt}{-12.946pt}{-6.14043pt}{-12.946pt}{-7.39755pt}\pgfsys@curveto{-12.946pt}{-8.65468pt}{-11.92691pt}{-9.67377pt}{-10.66978pt}{-9.67377pt}\pgfsys@curveto{-9.41266pt}{-9.67377pt}{-8.39357pt}{-8.65468pt}{-8.39357pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{-14.22638pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{{\mathcal{H}}_{-}^{A}}\otimes\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.16796875,0.3125,0.640625}\definecolor[named]{pgfstrokecolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@color@rgb@stroke{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@moveto{-1.28038pt}{-7.39755pt}\pgfsys@curveto{-1.28038pt}{-6.14043pt}{-2.29947pt}{-5.12134pt}{-3.5566pt}{-5.12134pt}\pgfsys@curveto{-4.81372pt}{-5.12134pt}{-5.83281pt}{-6.14043pt}{-5.83281pt}{-7.39755pt}\pgfsys@curveto{-5.83281pt}{-8.65468pt}{-4.81372pt}{-9.67377pt}{-3.5566pt}{-9.67377pt}\pgfsys@curveto{-2.29947pt}{-9.67377pt}{-1.28038pt}{-8.65468pt}{-1.28038pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{0.0pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{{\mathcal{H}}_{+}^{B}}\otimes\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.9375,0.23046875,0.15234375}\definecolor[named]{pgfstrokecolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@color@rgb@stroke{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@moveto{-1.28038pt}{-7.39755pt}\pgfsys@curveto{-1.28038pt}{-6.14043pt}{-2.29947pt}{-5.12134pt}{-3.5566pt}{-5.12134pt}\pgfsys@curveto{-4.81372pt}{-5.12134pt}{-5.83281pt}{-6.14043pt}{-5.83281pt}{-7.39755pt}\pgfsys@curveto{-5.83281pt}{-8.65468pt}{-4.81372pt}{-9.67377pt}{-3.5566pt}{-9.67377pt}\pgfsys@curveto{-2.29947pt}{-9.67377pt}{-1.28038pt}{-8.65468pt}{-1.28038pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{0.0pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{{\mathcal{H}}_{-}^{B}}\ . (8)

The Hilbert spaces on the right correspond to blue (+) particles on the left side, red (-) on the left, blue on the right, and red on the right, from left to right.

To simplify calculations, we introduce an equivalent set of different numbers: ii, the number of blue particles on the left side; N+N_{+}, the total number of blue particles; NAN_{A}, the total number of particles on the left; and NN, the total number of particles. From these, the total number of red particles, NN_{-}, and the total number of particles on the right, NBN_{B}, is determined by N++N=NA+NB=NN_{+}+N_{-}=N_{A}+N_{B}=N.

The volume of this macrostate, parametrized by the tuple (i,N+,NA,N)(i,N_{+},N_{A},N), is given by

VR=(wi)(wNAi)(wN+i)(wNN+(NAi)).V^{\mathrm{R}}=\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.16796875,0.3125,0.640625}\definecolor[named]{pgfstrokecolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@color@rgb@stroke{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@moveto{-8.39357pt}{-7.39755pt}\pgfsys@curveto{-8.39357pt}{-6.14043pt}{-9.41266pt}{-5.12134pt}{-10.66978pt}{-5.12134pt}\pgfsys@curveto{-11.92691pt}{-5.12134pt}{-12.946pt}{-6.14043pt}{-12.946pt}{-7.39755pt}\pgfsys@curveto{-12.946pt}{-8.65468pt}{-11.92691pt}{-9.67377pt}{-10.66978pt}{-9.67377pt}\pgfsys@curveto{-9.41266pt}{-9.67377pt}{-8.39357pt}{-8.65468pt}{-8.39357pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{-14.22638pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(w\\ i)}\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.9375,0.23046875,0.15234375}\definecolor[named]{pgfstrokecolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@color@rgb@stroke{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@moveto{-8.39357pt}{-7.39755pt}\pgfsys@curveto{-8.39357pt}{-6.14043pt}{-9.41266pt}{-5.12134pt}{-10.66978pt}{-5.12134pt}\pgfsys@curveto{-11.92691pt}{-5.12134pt}{-12.946pt}{-6.14043pt}{-12.946pt}{-7.39755pt}\pgfsys@curveto{-12.946pt}{-8.65468pt}{-11.92691pt}{-9.67377pt}{-10.66978pt}{-9.67377pt}\pgfsys@curveto{-9.41266pt}{-9.67377pt}{-8.39357pt}{-8.65468pt}{-8.39357pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{-14.22638pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(w\\ N_{A}-i)}\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.16796875,0.3125,0.640625}\definecolor[named]{pgfstrokecolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@color@rgb@stroke{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@moveto{-1.28038pt}{-7.39755pt}\pgfsys@curveto{-1.28038pt}{-6.14043pt}{-2.29947pt}{-5.12134pt}{-3.5566pt}{-5.12134pt}\pgfsys@curveto{-4.81372pt}{-5.12134pt}{-5.83281pt}{-6.14043pt}{-5.83281pt}{-7.39755pt}\pgfsys@curveto{-5.83281pt}{-8.65468pt}{-4.81372pt}{-9.67377pt}{-3.5566pt}{-9.67377pt}\pgfsys@curveto{-2.29947pt}{-9.67377pt}{-1.28038pt}{-8.65468pt}{-1.28038pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{0.0pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(w\\ N_{+}-i)}\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.9375,0.23046875,0.15234375}\definecolor[named]{pgfstrokecolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@color@rgb@stroke{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@moveto{-1.28038pt}{-7.39755pt}\pgfsys@curveto{-1.28038pt}{-6.14043pt}{-2.29947pt}{-5.12134pt}{-3.5566pt}{-5.12134pt}\pgfsys@curveto{-4.81372pt}{-5.12134pt}{-5.83281pt}{-6.14043pt}{-5.83281pt}{-7.39755pt}\pgfsys@curveto{-5.83281pt}{-8.65468pt}{-4.81372pt}{-9.67377pt}{-3.5566pt}{-9.67377pt}\pgfsys@curveto{-2.29947pt}{-9.67377pt}{-1.28038pt}{-8.65468pt}{-1.28038pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{0.0pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(w\\ N-N_{+}-(N_{A}-i))}. (9)

The binomial coefficients are the dimensions of the respective Hilbert spaces in Eq. (8), from left to right. We also denote this volume by a subscript “init” since we assume that when modelling time evolution, the system starts fully in this macrostate.

Assuming a diluted case, in which ww is significantly larger than other variables, we can apply Stirling’s approximation, as detailed in Appendix B.1. We derive

lnVR=N(lnwlnN+S(q)+1)12ln((2πN)4j=14qjmissing).\ln V^{\mathrm{R}}=N(\ln w-\ln N+S(q)+1)-\frac{1}{2}\ln\Big((2\pi N)^{4}\prod_{j=1}^{4}q_{j}\Big{missing}). (10)

Here, q1=piq_{1}=p_{i}, q2=apiq_{2}=a-p_{i}, q3=rpiq_{3}=r-p_{i}, and q4=1ra+piq_{4}=1-r-a+p_{i}, S(q)=iqilnqiS(q)=-\sum_{i}q_{i}\ln q_{i} denotes the Shannon entropy. We defined pi=i/Np_{i}=i/N as the proportion of blue particles on the left (out of all particles). We also define proportion of all particles on the left, a=NA/Na=N_{A}/N, and proportion of the blue particles out of all r=N+/Nr=N_{+}/N.

The last term of Eq. 10 represents logarithmic corrections. In those, the possible divergences when some qj=0q_{j}=0 result from the use of Stirling’s approximation. Upon closer inspection of these boundary cases one finds that the general formula can still be used by subtracting 12ln(2πNqj)-\frac{1}{2}\ln(2\pi Nq_{j}) for each qj=0q_{j}=0, see the end of Appendix B.1.

IV.2 Morty 1

Morty 1 can observe three macroscopic degrees of freedom: the total number of particles on the left NAN_{A}, the total number of particles on the right NNAN-N_{A}, and the total number of blue particles, N+N_{+}.

We again have N++N=NA+NB=NN_{+}+N_{-}=N_{A}+N_{B}=N. Without loss of generality, we define NAN_{A} and N+N_{+} to be the smaller of the two, i.e., NANBN_{A}\leq N_{B} and N+NN_{+}\leq N_{-}. This guarantees that NNA+N+N\geq N_{A}+N_{+} for the expression to be well defined.

M1=i=0min{N+,NA}R.{\mathcal{H}}^{\mathrm{M}_{1}}=\bigoplus_{i=0}^{\min\{N_{+},N_{A}\}}{\mathcal{H}}^{\mathrm{R}}. (11)

The volume corresponding this macrostate, parametrized by a tuple (N+,NA,N)(N_{+},N_{A},N), is given by the sum over all smaller macrostates (of Rick) that lead to the same macroscopic observation for Morty,

VM1=i=0min{N+,NA}VRV^{\mathrm{M}_{1}}=\sum_{i=0}^{\min\{N_{+},N_{A}\}}V^{\mathrm{R}} (12)

Assuming the dilute limit, when taking the logarithm, the dominant term will be only the largest term of the sum, while all others will contribute only as a logarithmic correction. To compute this correction, we use the Laplace method which involves approximating the binomials with a Gaussian around the maximum and approximating the sum with an integral. Finally, we obtain

lnVM1=N(lnwlnN+S(a)+S(r)+1)32ln(2πN).\begin{split}\ln V^{\mathrm{M}_{1}}&=N\left(\ln w-\ln N+S(a)+S(r)+1\right)\\ &-\frac{3}{2}\ln(2\pi N).\end{split} (13)

Recalling definitions a=NA/Na=N_{A}/N and r=N+/Nr=N_{+}/N, S(a)=alna(1a)ln(1a)S(a)=-a\ln a-(1-a)\ln(1-a) is the binary Shannon entropy.

See Appendix B.2 for the derivation and Appendix B.6 for the non-symmetric case.

IV.3 Morty 2

Morty 2 can observe two macroscopic degrees of freedom: the total number of particles on the left, NAN_{A} and the total number of particles on the right, NNAN-N_{A}. However, he is aware that two types of particles exist and could be in the system.

His macrostates are parametrized by a couple (NA,N)(N_{A},N). They are given by subspaces with corresponding volumes defined as their dimension,

M2=N+=0NM1,VM2=N+=0NVM1=2N+=0N/2i=0min{N+,NA}VR.\begin{split}{\mathcal{H}}^{\mathrm{M}_{2}}&=\bigoplus_{N_{+}=0}^{N}{\mathcal{H}}^{\mathrm{M}_{1}},\\ V^{\mathrm{M}_{2}}&=\sum_{N_{+}=0}^{N}V^{\mathrm{M}_{1}}=2\sum_{N_{+}=0}^{N/2}\sum_{i=0}^{\min\{N_{+},N_{A}\}}V^{\mathrm{R}}.\end{split} (14)

The right hand side of the volumes was split into two parts, due to the assumption of N+NN_{+}\leq N_{-} in Eq. (11).

Using the Laplace method on two variables, ii and N+N_{+}, we obtain

lnVM2=N(lnwlnN+S(a)+ln2+1)ln(4πN).\begin{split}\ln V^{\mathrm{M}_{2}}&=N\left(\ln w-\ln N+S(a)+\ln 2+1\right)-\ln(4\pi N).\end{split} (15)

(See Appendix B.3.) One can observe that up to the logarithmic correction, the leading order is the same as for Morty 1, with r=1/2r=1/2.

IV.4 Morty 3 - the ignorant

Finally, we consider the case of “ignorant” Morty, who not only cannot distinguish the two colors of particles, but is also not aware that these two colors of particles even exist. As a result, he also does not know how many particles of each type are in the system that he is observing. The only information he can access is the number of particles on the left and right sides, respectively.

This means that when calculating the associated observational entropy, he will assign different volumes to the macrostates than those corresponding to the actual number of microstates they contain, as given by Eq. (14).

From his perspective, the macrostate corresponding to having NAN_{A} particles on the left and NNAN-N_{A} particles on the right will have a volume given by

VM3=(wNA)(wNNA).V^{\mathrm{M}_{3}}=\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0,0,0}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@gray@stroke{0}\pgfsys@invoke{ }\pgfsys@color@gray@fill{0}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0,0,0}\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@moveto{-8.39357pt}{-7.39755pt}\pgfsys@curveto{-8.39357pt}{-6.14043pt}{-9.41266pt}{-5.12134pt}{-10.66978pt}{-5.12134pt}\pgfsys@curveto{-11.92691pt}{-5.12134pt}{-12.946pt}{-6.14043pt}{-12.946pt}{-7.39755pt}\pgfsys@curveto{-12.946pt}{-8.65468pt}{-11.92691pt}{-9.67377pt}{-10.66978pt}{-9.67377pt}\pgfsys@curveto{-9.41266pt}{-9.67377pt}{-8.39357pt}{-8.65468pt}{-8.39357pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{-14.22638pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(w\\ N_{A})}\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0,0,0}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@gray@stroke{0}\pgfsys@invoke{ }\pgfsys@color@gray@fill{0}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0,0,0}\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@moveto{-1.28038pt}{-7.39755pt}\pgfsys@curveto{-1.28038pt}{-6.14043pt}{-2.29947pt}{-5.12134pt}{-3.5566pt}{-5.12134pt}\pgfsys@curveto{-4.81372pt}{-5.12134pt}{-5.83281pt}{-6.14043pt}{-5.83281pt}{-7.39755pt}\pgfsys@curveto{-5.83281pt}{-8.65468pt}{-4.81372pt}{-9.67377pt}{-3.5566pt}{-9.67377pt}\pgfsys@curveto{-2.29947pt}{-9.67377pt}{-1.28038pt}{-8.65468pt}{-1.28038pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{0.0pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(w\\ N-N_{A})}. (16)

We obtain

lnVM3=N(lnwlnN+S(a)+1)12ln((2πN)2a(1a)).\begin{split}\ln V^{\mathrm{M}_{3}}&=N(\ln w-\ln N+S(a)+1)\\ &-\frac{1}{2}\ln((2\pi N)^{2}a(1-a)).\end{split} (17)

(See Appendix B.4.) The macrostate volume associated to the same observation is smaller than that of Morty 1 and Morty 2, by a factor of NS(r)NS(r) and Nln2N\ln 2, respectively. This is because Morty 3 has an incorrect judgement on the underlying physical Hilbert space. It is not because his uncertainty is, in reality, smaller, but because he thinks it is smaller. The factors are precisely those related to the existence of two types of particles, which cannot be distinguished locally. In the next section, we show, however, that under some reasonable assumptions, he will extract the same amount of energy as Morty 2 does when he executes the extraction protocol with his limited information. This amount is further consistent with his computed value of entropy. For ease of comparison, see Table 2 listing the abilities of each of the observers.

Observer Can access Cannot access Hilbert space
Rick ii, N+N_{+}, NAN_{A}, NBN_{B} R=+AA+BB{\mathcal{H}}^{\mathrm{R}}={\mathcal{H}}_{+}^{A}\otimes{\mathcal{H}}_{-}^{A}\otimes{\mathcal{H}}_{+}^{B}\otimes{\mathcal{H}}_{-}^{B}
Morty 1 N+N_{+}, NAN_{A}, NBN_{B} ii M1=i=0min{N+,NA}R{\mathcal{H}}^{\mathrm{M_{1}}}=\bigoplus_{i=0}^{\min\{N_{+},N_{A}\}}{\mathcal{H}}^{\mathrm{R}}
Morty 2 NAN_{A}, NBN_{B}, ii, N+N_{+} M2=N+=0NM1{\mathcal{H}}^{\mathrm{M}_{2}}=\bigoplus_{N_{+}=0}^{N}{\mathcal{H}}^{\mathrm{M}_{1}}
Morty 3 NAN_{A}, NBN_{B}
ii, N+N_{+}
(not aware)
M3={AB,(perceived)M2,(actual){\mathcal{H}}^{\mathrm{M}_{3}}=\begin{cases}{\mathcal{H}}^{A}\otimes{\mathcal{H}}^{B},&\text{(perceived)}\\ {\mathcal{H}}^{\mathrm{M}_{2}},&\text{(actual)}\end{cases}
Table 2: Different type of observer for the purpose of this work. We have tabulated the various macroscopic observables accessible (or inaccessible) to each of the observers, and their effective Hilbert spaces. Use of besides the observables of Morty 2 is to remind the reader that Morty 2 has knowledge of the existence of two colors, but he is unable to measure them (see the text in IV.1-IV.4). The labels denote: ii the number of blue particles on the left; N+N_{+} the total number of blue particles; NAN_{A} and NBN_{B}, the total number of particles on the left and the right, respectively.

IV.5 Entropy difference between observers

Consider a state which is contained within one of the macrostates of Rick, i.e., it has a definite macro-observable values (i,N+,NA,N)(i,N_{+},N_{A},N). For such a state, only one of the probabilities coming into the definition of observational entropy is equal to one while all others are zero. The same holds for Morty. Consequently, ignoring logarithmic terms, for Morty 1 we have

ΔS1=S𝒞M1S𝒞R=lnVM1lnVR=N(S(a)+S(r)S(q)).\begin{split}\Delta S_{1}&=S_{{\mathcal{C}_{\mathrm{M}}}_{1}}-S_{{\mathcal{C}_{\mathrm{R}}}}=\ln V^{\mathrm{M}_{1}}-\ln V^{\mathrm{R}}\\ &=N(S(a)+S(r)-S(q)).\end{split} (18)

For Morty 2, we have

ΔS2=S𝒞M2S𝒞R=lnVM2lnVR=N(S(a)+ln2S(q)).\begin{split}\Delta S_{2}&=S_{{\mathcal{C}_{\mathrm{M}}}_{2}}-S_{{\mathcal{C}_{\mathrm{R}}}}=\ln V^{\mathrm{M}_{2}}-\ln V^{\mathrm{R}}\\ &=N(S(a)+\ln 2-S(q)).\end{split} (19)

When there is the same number of blue and red particles, r=1/2r=1/2, while all blue are on the left and all red are on the right side, a=pi=1/2a=p_{i}=1/2, we obtain the long-expected result

ΔS1=ΔS2=Nln2,\Delta S_{1}=\Delta S_{2}=N\ln 2, (20)

which holds for both Morty 1 and 2. Later, we will insert these entropy differences into the formula for ergotropy difference, Eq. (7).

It would not be fair to compare Morty 3’s entropy with Rick’s, who clearly has a deeper understanding of the system. Rick may assign a higher entropy to the same state, but only because he accounts for internal degrees of freedom that Morty 3, due to his incorrect judgment, does not recognize. As noted above and elaborated below, this discrepancy does not prevent Morty 3 from correctly predicting the amount of extractable work—though it will generally be lower than what Rick could extract, given his fuller knowledge.

IV.6 Difference in entropy over time

In the following, we consider the initial state to be in one of the macrostates of Rick having a definite macro-observable values (i,N+,NA,N)(i,N_{+},N_{A},N). Further, we assume that the Hamiltonian preserves the total number of particles of each color and that the evolution will cause the particles to uniformly spread.

Rick. Assuming that the evolution will cause the particles to uniformly spread, in the long-time limit, the particles will uniformly fill the entire accessible Hilbert space, with probabilities proportional to the size of each macrostate. Mathematically, this means plR=VlR/Vp_{l}^{\mathrm{R}}=V_{l}^{\mathrm{R}}/V, where l=(i,N+,NA,N)l=(i,N_{+},N_{A},N) denotes the label of each macrostate, VlRV_{l}^{\mathrm{R}} is defined in Eq. (9), and the dimension of the entire accessible Hilbert space is,

V=(2wN+)(2wNN+).V=\matrixquantity(2w\\ N_{+})\matrixquantity(2w\\ N-N_{+}). (21)

This means that the final observational entropy equals

S𝒞Rfin=l(VlR/V)lnVlR/VVlR=lnV=N(ln2wlnN+S(r)+1)12ln(2πN)2r(1r).\begin{split}S_{{\mathcal{C}}_{R}}^{\mathrm{fin}}&=-\sum_{l}(V_{l}^{\mathrm{R}}/V)\ln\frac{V_{l}^{\mathrm{R}}/V}{V_{l}^{\mathrm{R}}}=\ln V\\ &=N(\ln 2w-\ln N+S(r)+1)-\frac{1}{2}\ln(2\pi N)^{2}r(1\!-\!r).\end{split} (22)

Ignoring the logarithmic terms, this gives us the difference between the final and the initial entropy,

ΔSR=S𝒞RfinS𝒞R=lnVVR=N(ln2+S(r)S(q)).\Delta S^{\mathrm{R}}=S_{{\mathcal{C}}_{R}}^{\mathrm{fin}}-S_{{\mathcal{C}}_{R}}=\ln\frac{V}{V^{\mathrm{R}}}=N(\ln 2+S(r)-S(q)). (23)

For the standard scenario of Gibbs mixing paradox depicted in Fig. 3, in which all blue particles start on the left, while the same number of the red particles starts on the right, corresponding to the initial configuration pi=1p_{i}=1 and r=a=1/2r=a=1/2, we have S(q)=ln2S(q)=\ln 2 and S(r)=ln2S(r)=\ln 2. This means that Rick will observe growth in entropy

ΔSR=Nln2,\Delta S^{\mathrm{R}}=N\ln 2, (24)

corresponding to the standard entropy of mixing, as the system thermalizes.

Morty. Assuming that the evolution will cause the particles to uniformly spread, for all clones of Morty 1 and 2, in the long-time limit the final entropy will be the same as Rick’s, Eq. (22), since they share the same Hilbert space. For Morty 3, it will be given by the reduced Hilbert space with dimension Vred.=(2wN)V_{\mathrm{red.}}=\binom{2w}{N},

S𝒞M3fin=lnVred.=N(ln2wlnN+1)12ln2πN.S_{{\mathcal{C}_{\mathrm{M}}}_{3}}^{\mathrm{fin}}=\ln V_{\mathrm{red.}}=N(\ln 2w-\ln N+1)-\frac{1}{2}\ln 2\pi N. (25)

Ignoring logarithmic terms, this gives entropy growth,

ΔSM1=S𝒞M1finS𝒞M1=lnVVM1=N(ln2S(a)),ΔSM2=lnVVM2=N(ln2S(a)),ΔSM3=lnVred.VM3=N(ln2S(a)).\begin{split}\Delta S^{\mathrm{M}_{1}}&=S_{{\mathcal{C}_{\mathrm{M}}}_{1}}^{\mathrm{fin}}-S_{{\mathcal{C}_{\mathrm{M}}}_{1}}=\ln\frac{V}{V^{\mathrm{M}_{1}}}=N(\ln 2-S(a)),\\ \Delta S^{\mathrm{M}_{2}}&=\ln\frac{V}{V^{\mathrm{M}_{2}}}=N(\ln 2-S(a)),\\ \Delta S^{\mathrm{M}_{3}}&=\ln\frac{V_{\mathrm{red.}}}{V^{\mathrm{M}_{3}}}=N(\ln 2-S(a)).\end{split} (26)

In the standard Gibbs mixing scenario, the initial state starts with configuration a=1/2a=1/2. This means that none of the clones will not observe a change in entropy,

ΔSM1=ΔSM2=ΔSM3=0.\Delta S^{\mathrm{M}_{1}}=\Delta S^{\mathrm{M}_{2}}=\Delta S^{\mathrm{M}_{3}}=0. (27)

Thus, irrespective of the correct or incorrect microstate counting, the difference in entropy of the final and the initial state is the same for all. For most applications in thermodynamics it does not matter where the initial–baseline entropy is; all that matters is the entropy difference. Thus, there is fundamentally no Gibbs paradox, considering that Morty 1 and 2, who have a proper understanding of the underlying Hilbert space, and Morty 3, who does not, will observe the same increase in entropy.

Finally, we take a look at the difference in entropy assigned by different versions of Morty,

S𝒞M1S𝒞M3=NS(r),S𝒞M2S𝒞M3=Nln2.\begin{split}S_{{\mathcal{C}_{\mathrm{M}}}_{1}}-S_{{\mathcal{C}_{\mathrm{M}}}_{3}}&=NS(r),\\ S_{{\mathcal{C}_{\mathrm{M}}}_{2}}-S_{{\mathcal{C}_{\mathrm{M}}}_{3}}&=N\ln 2.\end{split} (28)

For a system that contains only one type of particle, we have r=0r=0, Morty 3’s judgment becomes correct, and entropies associated by Morty 1 and 3 will coincide, recovering the correct limit. Since the system is conserving the number of particles of each type, the proportion of types rr is time-independent. In other words, for the entire time evolution, entropies associated to the system by the clones of Morty differ only by a time-independent value.

Refer to caption
Figure 4: Ergotropy difference between Rick, who can distinguish the particles, and Morty 1, who cannot, for variable system sizes. We fix k=1k=1. We compare the zeroth (ΔW0=TΔS\Delta W_{0}=T\Delta S), the first (ΔW1=TΔS(1+ΔS/2CE)\Delta W_{1}=T\Delta S(1+\Delta S/2{C_{\!E}})), and the second (ΔW2\Delta W_{2}, Eq. (7)) order expansion of ergotropy difference, with the exact value ΔW\Delta W, Eq. (4). ΔWav=(ΔW1+ΔW2)/2\Delta W_{\mathrm{av}}=(\Delta W_{1}+\Delta W_{2})/2 is the average value, which approximates the exact value the best. This is reasonable considering that the alternating signs of the expansion imply that the real value must fall somewhere in between ΔW1\Delta W_{1} and ΔW2\Delta W_{2}. As another heuristic approach, we employ a “smart” average by treating the expansion as a geometric series in the higher-order terms. This leads to ΔWav,2ΔW2w2q1+q\Delta W_{\mathrm{av},2}\equiv\Delta W_{2}-w_{2}\frac{q}{1+q}. The coefficients are defined as w1=ΔW1ΔW0w_{1}=\Delta W_{1}-\Delta W_{0}, w2=ΔW2ΔW1w_{2}=\Delta W_{2}-\Delta W_{1}, q=w2/w1q=w2/w1. We show the case of symmetric box sizes, where NA=NB=2N_{A}=N_{B}=2 and LA=LB=L/2L_{A}=L_{B}=L/2, and the non-symmetric case (inset), where NA=3N_{A}=3, NB=2N_{B}=2, LA=6,,17L_{A}=6,\dots,17, LB=NBNALAL_{B}=\lfloor\frac{N_{B}}{N_{A}}L_{A}\rfloor, and L=LA+LBL=L_{A}+L_{B}. All blue particles are on the left side, while all red ones are on the right. The average value closely matches the exact value. Additionally, the non-symmetric case shows clear convergence, confirming the expansion of the ergotropy difference.
Refer to caption
Figure 5: Observational ergotropy and entropy evolution through the mixing of a quantum gas. Evolved with Hamiltonian, Eq. (30), and initial configuration LA=8L_{A}=8, LB=6L_{B}=6, NA,B=LA,B/2N_{A,B}=L_{A,B}/2, with all blue (red) particles starting on the left (right). We plot the exact ergotropy difference and two approximate values: the average defined in Fig. (4) and the “smart” average, defined in caption 4. The smart average approximates the the exact value well for small values of ΔS\Delta S. Notice a small peak (dip) in ergotropy (entropy) around value t=4.2t=4.2, illustrating that the ergotropy is inversely related to the entropy difference (this follows from Eq. (3)). The initial and final entropy values are given by Rick’s initial configuration, Sinit=ln(84)(63)S^{\mathrm{init}}=\ln\binom{8}{4}\binom{6}{3}, and by the dimension of the accessible Hilbert space, Sfin=ln(144)(143)S^{\mathrm{fin}}=\ln\binom{14}{4}\binom{14}{3}, respectively.

V Work difference in quantum mixing

Inserting the entropy difference, Eq. (20), into Eq. (7), we obtain the difference between the amount of work Rick and Morty can extract in both quantum and classical Gibbs mixing scenario,

ΔW=NkTln2(1ln22c~E+16(1lnc~ETT)(ln2c~E)2+).\Delta W\!\!=\!NkT\ln\!2\!\left(\!1\!-\!\frac{\ln\!2}{2\tilde{c}_{E}}\!+\!\frac{1}{6}\!\left(1\!-\!\frac{\partial\ln\tilde{c}_{E}}{\partial T}T\right)\!\left(\frac{\ln\!2}{\tilde{c}_{E}}\right)^{\!2}\!\!\!+{\!\cdot\!\cdot\!\cdot\!}\right)\!. (29)

where c~E=CE/(Nk)\tilde{c}_{E}={C_{\!E}}/(Nk) is the dimensionless specific heat. We reiterate that temperature TT is not a temperature of a thermal bath. Instead, it is the observational temperature given by S(ρT)=S𝒞MS({\rho}_{T})=S_{{\mathcal{C}}_{M}}. For an ideal gas described in Table 1, we obtain ΔW=NkTln2(113ln2+227(ln2)2+)0.805×NkTln2\Delta W=NkT\ln 2(1-\frac{1}{3}\ln 2+\frac{2}{27}(\ln 2)^{2}+\cdots)\approx 0.805\times NkT\ln 2.

We observe that the formula for Rick’s entropy difference at different times, Eq. (24), is the same as the aforementioned entropy difference. Following the arguments of Section III.2, this shows that Rick will extract ΔW\Delta W less energy if he chooses to perform the extraction at a later time, as opposed to immediately.

While different types of Morty associate different entropy, this difference remains constant throughout the evolution, Eq. (28). Consequently, the increase in entropy will be the same for all of them, Eq. (26). The extracted work difference depends only on this entropy increase. There is no advantage for either Morty 1 or 2; they will extract the same amount of energy.

Morty 3 requires much more care. This is because he lacks a correct understanding of the underlying Hilbert space, resulting from his unawareness of the existence of the two particles, due to his inability to distinguish between different colors (which is his defining feature). This we explained in Section IV.4. As a result, he further misjudges the actual unitary he applies to the system when attempting to extract work, misjudges the measurement he performs, and would infer a density matrix from quantum tomography that deviates from the true one. However, in Appendix C we prove that under some reasonable assumptions, he will extract the same amount of energy as Morty 2. This includes three consistency assumptions: one on the density matrix, stating that the actual density is macroscopically consistent with the one he would infer by performing any measurement that is unable to distinguish the colors of the particles, which we call an admissible measurement. Second, on the measurements, saying that the probability of distribution outcomes he predicts is the same as the actual one when he performs any admissible measurement. Third, on the unitary operations, saying that the actual unitary measurement is consistent with what is expected and detectable by admissible measurements. Furthermore, we have the super-selection rule assumptions, which state that it is not possible to have superpositions between different types of particles. This means that states such as |0++|0\ket{0+}+\ket{-0} are not allowed, nor can they be created by applying either measurement or unitary operations, or by letting the system evolve naturally under the Hamiltonian. So far, these assumptions have been relatively weak. Finally, we have one non-trivial assumption that will not be true for many systems, namely that the Hamiltonian does not depend on the color of the particles but only on the configuration—describing whether the particle is or is not at a given lattice site, irrespective of its color. This means that every two states with the same configuration, such as |0+\ket{0+} and |0\ket{0-} corresponding to |01\ket{01}, must have the same energy. This ensures that Morty 3 correctly understands the system’s energy, given the constraints of what he is able to measure.

With these assumptions, Morty 3 extracts exactly the same amount of energy as Morty 2. Thus, the difference in extracted work between the two times will be the same.

To conclude, there is no advantage for either clone of Morty. In other words, it does not matter whether Morty knows or does not know whether there are two types of particles in the system rather than just one, or how many blue particles there are in total. In either case, he will extract the same amount of energy. This resolves the Gibbs mixing paradox: while absolute values of associated entropy might be different, there is no difference in terms of what different clones of Morty can practically do.

We provide two simulations of the ergotropy difference: Static, applied only to the initial state which is localized in one of Rick’s macrostates but for varying system sizes, see Fig. 4; and dynamic, where the system starts in an initial macrostate of Rick but is then evolved with Hamiltonian that preserves the total number of particles of each kind and distributes them uniformly in the long-time limit, see Fig. 5.

The Hamiltonian is given by

H=it1cici+1+v1nini+1+it2cici+2+v2nini+2+h.c.\begin{split}H&=\sum_{i}t_{1}c^{\dagger}_{i}c_{i+1}+v_{1}n_{i}n_{i+1}\\ &+\sum_{i}t_{2}c^{\dagger}_{i}c_{i+2}+v_{2}n_{i}n_{i+2}+\text{h.c.}\end{split} (30)

with periodic boundary conditions. cic_{i}^{\dagger} and cic_{i} represent creation and annihiliation operators at site ii, respectively, and nin_{i} represents the number operator. The numerical computations are done in the chaotic phase, with t1=v1=1t_{1}=v_{1}=1 and t2=v2=0.96t_{2}=v_{2}=0.96.

To minimize the error, we compute the entropy exactly from the binomial coefficients, Eqs. (9) and (12). Further, we compute the ergotropy difference by inserting this entropy difference into the exact formula, Eq. (4), and then from varying levels of approximation using the perturbative expansion, Eq. (7).

VI Conclusion

In this work, we investigated how two observers—each characterized by a distinct ability to distinguish macrostates—assign entropy to an isolated quantum or classical system, and how this observer-dependent entropy influences the energy they can extract. Focusing on the difference in extractable energy between the two observers, we found that its perturbative expansion closely resembles Landauer’s bound, with higher-order corrections determined by the ratio of entropy to heat capacity. Notably, rather than the usual temperature of a thermal bath, an “observational temperature” emerges, governed by the observers’ knowledge of the system.

In the case of critical systems, we identified a critical amount of knowledge—reached when the observational temperature equals the critical temperature—beyond which the amount of extractable work jumps discontinuously. We further demonstrated that the same formula for the extracted energy difference applies when considering a single observer extracting energy at different times.

We then applied these ideas to the Gibbs mixing paradox, in which one observer can distinguish two types of particles while another cannot. A closer examination revealed three distinct subtypes of the latter observer, each associating a different amount of entropy to the system.

However, all of them observe the same entropy increase. For the two observers who are aware of the existence of two types of particles, this means that they extract the same amount of work. For the third, who is not aware of their existence, we assume consistency of their description, superselection rules, and independence of the Hamiltonian on different particle types. With these assumptions, we show that they extract exactly the same amount of energy as the second observer, who is aware of the existence of two particle types, but unable to distinghuish them. This resolves the Gibbs mixing paradox as follows. Even though an observer who discovers the presence of “two colors” of particles experiences a discontinuous jump in the entropy they assign to the system (counterintuitively, knowing more but associating higher entropy!), this new knowledge does not provide a practical advantage for work extraction if the particles remain indistinguishable in their measurements.

We further argued that the loss of extractable energy caused by mixing must also place a lower bound on the energy required to separate the particles, though there are subtleties in reversing the protocol that merit future investigation. Our analysis was restricted to projective measurements; generalizing to POVMs [38, 61], including non-commuting measurements [62, 63] would be an interesting avenue for follow-up work. While free energy is the relevant quantity in standard mixing scenarios, we showed that “observational ergotropy” plays an analogous role in the isolated systems considered here. It would be worthwhile to explore whether these concepts overlap, potentially through studying local equilibrium states where a combination of the local particle number (as considered here) and the local energy coarse-grainings define a “dynamical thermodynamic entropy [64].” Further, as dynamical thermodynamic entropy corresponds with thermodynamic entropy for systems in equilibrium, it would be interesting to compare the work difference derived from this work and that from the adiabatic expansion of ideal gas. This would serve as a check for the consistency of the derived theory. Another promising direction is to examine Fisher information associated with temperature measurements, given its relation to the energy variance that appears in the perturbative expansion of the extracted work.

Finally, we propose two experimental platforms which allow the manipulation of two conserved, distinguishable particle states to test these ideas in a quantum setting. First, we can use Bose-Einstein condensates (BECs), e.g. of Rb87{}^{87}\text{Rb}, where hyperfine states (as “colors”) are spatially separated by applying a magnetic field to induce Zeeman splitting [65, 66, 67]. The evolution conserves the density of each type, and mixing can be created by modulating magnetic field [68]. The second promising platform is a 2D optical lattice with Lithium atoms in hyperfine states to denote |\ket{\uparrow} and |\ket{\downarrow}. Each lattice site holds either no atom, one atom of either spin, or two atoms of opposite spin. Using a spin-removal beam and site-resolved imaging [69, 70], we can determine the spin number, type, and mixing at each site, enabling the computation of observational entropy and extractable work. To realize the work-extraction protocol in practice would require further modifications—such as developing controls for the random and extraction unitaries—or adapting the protocol to fit experimental constraints.

Acknowledgments. B.B., R.K.R., and D.Š. acknowledge the support from the Institute for Basic Science in Korea (IBS-R024-D1). Many thanks to Felix C. Binder for his insightful comments and questions, which have shaped the direction of this research and improved its overall presentation.

Appendix A Derivation of the work difference between observers with similar capabilities

First, we derive expansion of Δβ\Delta\beta in terms of ΔS\Delta S. We do this by expanding Eq. (2) as,

16S′′′Δβ3+12S′′Δβ2+SΔβ+ΔS=0.\frac{1}{6}S^{\prime\prime\prime}\Delta\beta^{3}+\frac{1}{2}S^{\prime\prime}\Delta\beta^{2}+S^{\prime}\Delta\beta+\Delta S=0. (31)

There, the coefficients represent the derivative of the von Neumann entropies, S(k)dkSvN(ρβ)dβkS^{(k)}\equiv\frac{d^{k}S_{{\mathrm{vN}}}({\rho}_{\beta})}{d\beta^{k}}. Explicitly, we derive

S\displaystyle S^{\prime} =βv\displaystyle=-\beta v (32)
S′′\displaystyle S^{\prime\prime} =vβx\displaystyle=-v-\beta x (33)
S′′′\displaystyle S^{\prime\prime\prime} =2xβy,\displaystyle=-2x-\beta y, (34)

where

v\displaystyle v =dEdβ=E2E2,\displaystyle=-\frac{d\langle E\rangle}{d\beta}=\langle E^{2}\rangle-\langle E\rangle^{2}, (35)
x\displaystyle x =dvdβ=E3+3EE22E3,\displaystyle=\frac{dv}{d\beta}=-\langle E^{3}\rangle+3\langle E\rangle\langle E^{2}\rangle-2\langle E\rangle^{3}, (36)
y\displaystyle y =dxdβ=E44EE33E22\displaystyle=\frac{dx}{d\beta}=\langle E^{4}\rangle-4\langle E\rangle\langle E^{3}\rangle-3\langle E^{2}\rangle^{2} (37)
+12E2E26E4,\displaystyle\ \ \ \ \ \ \ \ \;\!+12\langle E\rangle^{2}\langle E^{2}\rangle-6\langle E\rangle^{4}, (38)

are the variance of the thermal state, the first derivative of the variance, and the second derivative, respectively. To derive this, we have used SvN(ρβ)=lnZ+βES_{{\mathrm{vN}}}({\rho}_{\beta})=\ln Z+\beta\langle E\rangle, dZdβ=EZ\frac{dZ}{d\beta}=-\langle E\rangle Z, and dEkdβ=Ek+1+EEk\frac{d\langle E^{k}\rangle}{d\beta}=-\langle E^{k+1}\rangle+\langle E\rangle\langle E^{k}\rangle.

We rewrite Eq. (31) as

β(2X+βY)6b3+1+βX2b2+bδ=0.\frac{\beta(2X+\beta Y)}{6}b^{3}+\frac{1+\beta X}{2}b^{2}+b-\delta=0. (39)

We have defined b=Δβ/βb=\Delta\beta/\beta, X=x/vX=x/v, Y=y/vY=y/v, and δ=ΔS/(vβ2)\delta=\Delta S/(v\beta^{2}). We look for the solution of in the small expansion of δ\delta, and derive

Δβ=βδ(11+βX2δ+3+4βX+3β2X2β2Y6δ2+).\begin{split}\Delta\beta\!=\!\beta{\delta}\left(1\!-\!\frac{1\!+\!\beta X}{2}{\delta}\!+\!\frac{3\!+\!4\beta X\!+\!3\beta^{2}X^{2}\!-\!\beta^{2}Y}{6}{\delta}^{2}\!+{\!\cdot\!\cdot\!\cdot\!}\right).\\ \end{split} (40)

Similarly, we expand the work difference in the expansion of Δβ\Delta\beta as,

ΔW=vΔβ(1+12XΔβ+16YΔβ2+).\Delta W=v\Delta\beta\left(1+\frac{1}{2}X\Delta\beta+\frac{1}{6}Y\Delta\beta^{2}+\cdots\right). (41)

There, we have used ΔW=EβEβR\Delta W=\langle E\rangle_{\beta}-\langle E\rangle_{{\beta_{\mathrm{R}}}} (note that ββM\beta\equiv{\beta_{\mathrm{M}}}), where EEβ=tr[Hρβ]\langle E\rangle\equiv\langle E\rangle_{\beta}=\tr[{H}{\rho}_{\beta}], and EβR=EvΔβ12xΔβ216yΔβ3+\langle E\rangle_{{\beta_{\mathrm{R}}}}=\langle E\rangle-v\Delta\beta-\frac{1}{2}x\Delta\beta^{2}-\frac{1}{6}y\Delta\beta^{3}+{\!\cdot\!\cdot\!\cdot\!}. Combining the above two equations, we obtain our final result in the expansion of δ\delta,

ΔW=vβδ(1δ2+3+Xβ6δ2+).\Delta W=v\beta\delta\left(1-\frac{\delta}{2}+\frac{3+X\beta}{6}\delta^{2}+\cdots\right). (42)

Inserting the substitutions back together with β=1/(kT)\beta=1/(kT), we obtain

ΔW=kTΔS(1(kT)22vΔS+(3kT+X)(kT)36v2ΔS2+).\Delta W\!=\!kT\Delta S\!\left(1\!-\!\frac{(kT)^{2}}{2v}\Delta S\!+\!\frac{(3kT\!+\!X)(kT)^{3}}{6v^{2}}\Delta S^{2}\!+{\!\cdot\!\cdot\!\cdot\!}\right)\!. (43)

which is Eq. (6). Note that the final result does not depend on the fourth moment of energy, YY, as its contributions cancel out.

The variance of the thermal state is connected to the heat capacity,

v=dEdβ=ETTβ=CE1kββ=CEkβ2=CEkT2,v=-\frac{d\langle E\rangle}{d\beta}=-\frac{\partial\langle E\rangle}{\partial T}\frac{\partial T}{\partial\beta}=-{C_{\!E}}\frac{\partial\tfrac{1}{k\beta}}{\partial\beta}=\frac{{C_{\!E}}}{k\beta^{2}}={C_{\!E}}kT^{2}, (44)

where TT is the temperature, and CE=E/T{C_{\!E}}=\partial E/\partial T is the heat capacity. Similarly, we derive an expression for xx as

x=dvdβ=(CEkT2)TkT2=CET(kT2)22CEk2T3.x=\frac{dv}{d\beta}=-\frac{\partial({C_{\!E}}kT^{2})}{\partial T}kT^{2}=-\frac{\partial{C_{\!E}}}{\partial T}(kT^{2})^{2}-2{C_{\!E}}k^{2}T^{3}. (45)

This gives

X=xv=1CECETkT22kT=ln(CE/k)TkT22kT.X=\frac{x}{v}=-\frac{1}{{C_{\!E}}}\frac{\partial{C_{\!E}}}{\partial T}kT^{2}-2kT=-\frac{\partial\ln({C_{\!E}}/k)}{\partial T}kT^{2}-2kT. (46)

Finally, we obtain Eq. (7),

ΔW=kTΔS(1kΔS2CE+1ln(CE/k)TT6(kΔSCE)2+).\Delta W=\!kT\!\Delta S\left(1\!-\!\frac{k\Delta S}{2{C_{\!E}}}\!+\!\frac{1\!-\!\frac{\partial\ln({C_{\!E}}/k)}{\partial T}T}{6}\left(\frac{k\Delta S}{{C_{\!E}}}\right)^{2}\!\!\!\!+{\!\cdot\!\cdot\!\cdot\!}\right)\!. (47)

Appendix B Volumes of Rick and Morty macrostates

B.1 Rick

Rick can observe both the number and type of particle in each part. His macrostate is given by four numbers: ii, denoting the number of blue particles on the left; N+N_{+}, the total number of blue particles; NAN_{A}, the total number of particles on the left; and NN, the total number of particles. From these, the total number of red particles, NN_{-}, and the total number of particles on the right, NBN_{B}, can be derived using N++N=NA+NB=NN_{+}+N_{-}=N_{A}+N_{B}=N.

Volume this macrostate is given by

VR=(wi)(wNAi)(wN+i)(wNN+(NAi)).V^{\mathrm{R}}=\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.16796875,0.3125,0.640625}\definecolor[named]{pgfstrokecolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@color@rgb@stroke{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@moveto{-8.39357pt}{-7.39755pt}\pgfsys@curveto{-8.39357pt}{-6.14043pt}{-9.41266pt}{-5.12134pt}{-10.66978pt}{-5.12134pt}\pgfsys@curveto{-11.92691pt}{-5.12134pt}{-12.946pt}{-6.14043pt}{-12.946pt}{-7.39755pt}\pgfsys@curveto{-12.946pt}{-8.65468pt}{-11.92691pt}{-9.67377pt}{-10.66978pt}{-9.67377pt}\pgfsys@curveto{-9.41266pt}{-9.67377pt}{-8.39357pt}{-8.65468pt}{-8.39357pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{-14.22638pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(w\\ i)}\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.9375,0.23046875,0.15234375}\definecolor[named]{pgfstrokecolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@color@rgb@stroke{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@moveto{-8.39357pt}{-7.39755pt}\pgfsys@curveto{-8.39357pt}{-6.14043pt}{-9.41266pt}{-5.12134pt}{-10.66978pt}{-5.12134pt}\pgfsys@curveto{-11.92691pt}{-5.12134pt}{-12.946pt}{-6.14043pt}{-12.946pt}{-7.39755pt}\pgfsys@curveto{-12.946pt}{-8.65468pt}{-11.92691pt}{-9.67377pt}{-10.66978pt}{-9.67377pt}\pgfsys@curveto{-9.41266pt}{-9.67377pt}{-8.39357pt}{-8.65468pt}{-8.39357pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{-14.22638pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(w\\ N_{A}-i)}\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.16796875,0.3125,0.640625}\definecolor[named]{pgfstrokecolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@color@rgb@stroke{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@moveto{-1.28038pt}{-7.39755pt}\pgfsys@curveto{-1.28038pt}{-6.14043pt}{-2.29947pt}{-5.12134pt}{-3.5566pt}{-5.12134pt}\pgfsys@curveto{-4.81372pt}{-5.12134pt}{-5.83281pt}{-6.14043pt}{-5.83281pt}{-7.39755pt}\pgfsys@curveto{-5.83281pt}{-8.65468pt}{-4.81372pt}{-9.67377pt}{-3.5566pt}{-9.67377pt}\pgfsys@curveto{-2.29947pt}{-9.67377pt}{-1.28038pt}{-8.65468pt}{-1.28038pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{0.0pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(w\\ N_{+}-i)}\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.9375,0.23046875,0.15234375}\definecolor[named]{pgfstrokecolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@color@rgb@stroke{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@moveto{-1.28038pt}{-7.39755pt}\pgfsys@curveto{-1.28038pt}{-6.14043pt}{-2.29947pt}{-5.12134pt}{-3.5566pt}{-5.12134pt}\pgfsys@curveto{-4.81372pt}{-5.12134pt}{-5.83281pt}{-6.14043pt}{-5.83281pt}{-7.39755pt}\pgfsys@curveto{-5.83281pt}{-8.65468pt}{-4.81372pt}{-9.67377pt}{-3.5566pt}{-9.67377pt}\pgfsys@curveto{-2.29947pt}{-9.67377pt}{-1.28038pt}{-8.65468pt}{-1.28038pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{0.0pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(w\\ N-N_{+}-(N_{A}-i))}. (48)

The Binomial coefficients correspond to counting the number of microstates for blue (+) particles on the left side, red (-) on the left, blue on the right, and red on the right, from left to right.

We compute this volume by using Stirling’s approximation,

ln(wk)=wlnwklnk(wk)ln(wk)+12lnw2πk(wk),\begin{split}\ln\matrixquantity(w\\ k)&=w\ln w-k\ln k-(w-k)\ln(w-k)\\ &+\frac{1}{2}\ln\frac{w}{2\pi k(w-k)},\end{split} (49)

Assuming that ww is significantly larger than other parameters, we approximate terms such as

ln(wk)lnwkw.\ln(w-k)\approx\ln w-\frac{k}{w}. (50)

This gives

ln(wk)=wlnwklnk(wk)(lnwkw)12ln2πk+12kw,=klnwklnk+kk2w12ln2πk+12kw,klnwklnk+k12ln2πk,\begin{split}\ln\matrixquantity(w\\ k)&=w\ln w-k\ln k-(w-k)\left(\ln w-\frac{k}{w}\right)\\ &-\frac{1}{2}\ln 2\pi k+\frac{1}{2}\frac{k}{w},\\ &=k\ln w-k\ln k+k-\frac{k^{2}}{w}-\frac{1}{2}\ln 2\pi k+\frac{1}{2}\frac{k}{w},\\ &\approx k\ln w-k\ln k+k-\frac{1}{2}\ln 2\pi k,\end{split} (51)

where in the last term we neglected terms divided by ww.

For any 0<p10<p\leq 1, we can conveniently rewrite Eq. (51) as

ln(wpN)=pN(lnwlnNlnp+1)12ln2πpN.\ln\matrixquantity(w\\ p\;\!N)=pN(\ln w-\ln N-\ln p+1)-\frac{1}{2}\ln 2\pi pN. (52)

We define pi=i/Np_{i}=i/N as the proportional amount of blue particles on the left out of all particles. We also define the proportion of the total particles on the left, a=NA/Na=N_{A}/N, and the proportion of the blue particles r=N+/Nr=N_{+}/N.

We obtain

lnVR=N(lnwlnN+S(q)+1)12ln((2πN)4j=14qjmissing)\ln V^{\mathrm{R}}=N(\ln w-\ln N+S(q)+1)-\frac{1}{2}\ln\Big((2\pi N)^{4}\prod_{j=1}^{4}q_{j}\Big{missing}) (53)

where q1=piq_{1}=p_{i}, q2=apiq_{2}=a-p_{i}, q3=rpiq_{3}=r-p_{i}, and q4=1ra+piq_{4}=1-r-a+p_{i}.

Note that for aa or rr equal to zero, the constant term is not defined. However, neither is Eq. (51) for k=0k=0 due to the limits of the validity of Stirling’s approximation. However, for any tiny but non-zero aa and rr, the equation is valid and very precise. For example, when a=1/Na=1/N corresponding to NA=1N_{A}=1, the problematic term becomes lna(1a)N=ln(11/N)1/N\ln a(1-a)N=\ln(1-1/N)\approx 1/N. Consequently, terms will be small compared to the dominant NlnwN\ln w term.

To investigate this situation in detail, consider the situation in which pi=0p_{i}=0. As a result, the volume consists of three terms,

VR=(wNA)(wN+)(wNN+NA)V^{\mathrm{R}}=\matrixquantity(w\\ N_{A})\matrixquantity(w\\ N_{+})\matrixquantity(w\\ N-N_{+}-N_{A}) (54)

We analogously derive

lnVinitR(i=0)=N(lnwlnN+S(q)+1)12ln((2πN)3j=24qjmissing),\begin{split}\ln V^{\mathrm{R}}_{\mathrm{init}}(i=0)&=N(\ln w-\ln N+S(q)+1)\\ &-\frac{1}{2}\ln\Big((2\pi N)^{3}\prod_{j=2}^{4}q_{j}\Big{missing}),\end{split} (55)

where q2=aq_{2}=a, q3=rq_{3}=r, and q4=1raq_{4}=1-r-a. Thus, the general form is still valid and equal to Eq. (53), considering that one takes out 12ln(2πNqj)-\frac{1}{2}\ln(2\pi Nq_{j}) for each qj=0q_{j}=0.

B.2 Morty 1

Morty can observe only two macroscopic degrees of freedom, that is, the total number of particles on the left NAN_{A}, and the total number of particles on the right, which we denote as the total number of particles minus those on the left, NNAN-N_{A}.

We again have N++N=NA+NB=NN_{+}+N_{-}=N_{A}+N_{B}=N. Without loss of generatity we define NAN_{A} and N+N_{+} to be the smaller of the two, i.e., NANBN_{A}\leq N_{B} and N+NN_{+}\leq N_{-}. This guarantees that NNA+N+N\geq N_{A}+N_{+} for the expression to be well defined. The volume corresponding the this macrostate is the sum over all smaller macrostates (of Rick) that lead to the same macroscopic observation for Morty,

VM1=i=0min{N+,NA}(wi)(wNAi)(wN+i)(wNN+(NAi))=:i=0min{N+,NA}f(i).\begin{split}&V^{\mathrm{M}_{1}}=\\ &\sum_{i=0}^{\min\{N_{+},N_{A}\}}\!\!\!\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.16796875,0.3125,0.640625}\definecolor[named]{pgfstrokecolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@color@rgb@stroke{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@moveto{-8.39357pt}{-7.39755pt}\pgfsys@curveto{-8.39357pt}{-6.14043pt}{-9.41266pt}{-5.12134pt}{-10.66978pt}{-5.12134pt}\pgfsys@curveto{-11.92691pt}{-5.12134pt}{-12.946pt}{-6.14043pt}{-12.946pt}{-7.39755pt}\pgfsys@curveto{-12.946pt}{-8.65468pt}{-11.92691pt}{-9.67377pt}{-10.66978pt}{-9.67377pt}\pgfsys@curveto{-9.41266pt}{-9.67377pt}{-8.39357pt}{-8.65468pt}{-8.39357pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{-14.22638pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(w\\ i)}\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.9375,0.23046875,0.15234375}\definecolor[named]{pgfstrokecolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@color@rgb@stroke{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@moveto{-8.39357pt}{-7.39755pt}\pgfsys@curveto{-8.39357pt}{-6.14043pt}{-9.41266pt}{-5.12134pt}{-10.66978pt}{-5.12134pt}\pgfsys@curveto{-11.92691pt}{-5.12134pt}{-12.946pt}{-6.14043pt}{-12.946pt}{-7.39755pt}\pgfsys@curveto{-12.946pt}{-8.65468pt}{-11.92691pt}{-9.67377pt}{-10.66978pt}{-9.67377pt}\pgfsys@curveto{-9.41266pt}{-9.67377pt}{-8.39357pt}{-8.65468pt}{-8.39357pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{-14.22638pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(w\\ N_{A}\!-\!i)}\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.16796875,0.3125,0.640625}\definecolor[named]{pgfstrokecolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@color@rgb@stroke{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@moveto{-1.28038pt}{-7.39755pt}\pgfsys@curveto{-1.28038pt}{-6.14043pt}{-2.29947pt}{-5.12134pt}{-3.5566pt}{-5.12134pt}\pgfsys@curveto{-4.81372pt}{-5.12134pt}{-5.83281pt}{-6.14043pt}{-5.83281pt}{-7.39755pt}\pgfsys@curveto{-5.83281pt}{-8.65468pt}{-4.81372pt}{-9.67377pt}{-3.5566pt}{-9.67377pt}\pgfsys@curveto{-2.29947pt}{-9.67377pt}{-1.28038pt}{-8.65468pt}{-1.28038pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{0.0pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(w\\ N_{+}\!-\!i)}\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.9375,0.23046875,0.15234375}\definecolor[named]{pgfstrokecolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@color@rgb@stroke{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@moveto{-1.28038pt}{-7.39755pt}\pgfsys@curveto{-1.28038pt}{-6.14043pt}{-2.29947pt}{-5.12134pt}{-3.5566pt}{-5.12134pt}\pgfsys@curveto{-4.81372pt}{-5.12134pt}{-5.83281pt}{-6.14043pt}{-5.83281pt}{-7.39755pt}\pgfsys@curveto{-5.83281pt}{-8.65468pt}{-4.81372pt}{-9.67377pt}{-3.5566pt}{-9.67377pt}\pgfsys@curveto{-2.29947pt}{-9.67377pt}{-1.28038pt}{-8.65468pt}{-1.28038pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{0.0pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(w\\ N\!-\!N_{+}\!-\!(N_{A}\!-\!i))}\\ &=:\sum_{i=0}^{\min\{N_{+},N_{A}\}}f(i).\end{split} (56)

We find the largest terms of the sum assuming that ww is much larger than the rest. Then, using Stirling’s approximation, Eq. (51), we express lnf(i)\ln f(i). Setting the derivative of this function to zero, we find imax=N+NANi_{\max}=\frac{N_{+}N_{A}}{N}. Clearly, 0imaxN+,NA0\leq i_{\max}\leq N_{+},N_{A}, so this value gives a valid maximum. Since logarithm is a strictly growing function, imaxi_{\max} also gives the maximum of f(i)f(i).

Recalling definitions a=NA/Na=N_{A}/N and r=N+/Nr=N_{+}/N, we can elegantly write imax=arNi_{\max}=arN. Using Eq. (52), we obtain the maximal term,

lnfmax=N(lnwlnN+S(a)+S(r)+1)12ln((2πN)4a(1a)r(1r)).\begin{split}\ln f_{\max}&=N\left(\ln w-\ln N+S(a)+S(r)+1\right)\\ &-\frac{1}{2}\ln((2\pi N)^{4}a(1-a)r(1-r)).\end{split} (57)

Next, we include the effect of the sum, which will contribute at most logarithmically in NN. We use the Laplace method to approximate the binomial coefficients around the maximum and then estimate the sum by computing the integral. Mathematically, we have

f(i)=fmaxexp((iimax)22σ2)f(i)=f_{\max}\exp\left(-\frac{(i-i_{\max})^{2}}{2\sigma^{2}}\right) (58)

We compute the second derivative to compute the variance,

σ2=(d2lnfdi2|i=imax)1=Na(1a)r(1r).\sigma^{2}=-\left(\left.\frac{d^{2}\ln f}{di^{2}}\right|_{i=i_{\max}}\right)^{-1}=Na(1-a)r(1-r). (59)

Finally, we obtain the sum,

VM1=+f(i)=fmax2πNa(1a)r(1r).\begin{split}V^{\mathrm{M}_{1}}&=\int_{-\infty}^{+\infty}f(i)=f_{\max}\sqrt{2\pi Na(1-a)r(1-r)}.\end{split} (60)

Collecting all results, we obtain the macrostate volume,

lnVM1=N(lnwlnN+S(a)+S(r)+1)32ln(2πN).\begin{split}\ln V^{\mathrm{M}_{1}}&=N\left(\ln w-\ln N+S(a)+S(r)+1\right)\\ &-\frac{3}{2}\ln(2\pi N).\end{split} (61)

B.3 Morty 2

The macrostate volume of Morty 2, parametrized by a couple (NA,N)(N_{A},N), is given by

VM2=2N+=0N/2i=0min{N+,NA}f(i,N+).V^{\mathrm{M}_{2}}=2\sum_{N_{+}=0}^{N/2}\sum_{i=0}^{\min\{N_{+},N_{A}\}}f(i,N_{+}). (62)

where f(i,N+)f(i,N_{+}) is a function identical to that of f(i)f(i) in Eq. (56) and equal to the Rick’s macrostate volume, Eq. (48).

The maximum of f(i,N+)f(i,N_{+}) is reached at imax=aN/2i_{\max}=aN/2 and R+max=N/2R_{+\max}=N/2, and given by

lnfmax=N(lnwlnN+S(a)+ln2+1)12ln((2πN)44a(1a)).\begin{split}\ln f_{\max}&=N\left(\ln w-\ln N+S(a)+\ln 2+1\right)\\ &-\frac{1}{2}\ln(\frac{(2\pi N)^{4}}{4}a(1-a)).\end{split} (63)

(Obtained by inserting r=1/2r=1/2 into Eq. (63).)

We got a Gaussian approximation

f(i,N+)fmaxexp(12(ΔiΔN+)H(ΔiΔN+)),f(i,N_{+})\approx f_{\max}\exp\left(\frac{1}{2}\begin{pmatrix}\Delta i&\Delta N_{+}\end{pmatrix}H\begin{pmatrix}\Delta i\\ \Delta N_{+}\end{pmatrix}\right), (64)

where Δi=iimax\Delta i=i-i_{\max}, ΔN+=N+R+max\Delta N_{+}=N_{+}-R_{+\max}, and the Hessian is defined as,

H=(2lnfi22lnfiN+2lnfN+i2lnfN+2)|(imax,R+max)=4N(1a(1a)11a11a1a(1a)).H=\left.\begin{pmatrix}\displaystyle\frac{\partial^{2}\ln f}{\partial i^{2}}&\displaystyle\frac{\partial^{2}\ln f}{\partial i\partial N_{+}}\\[12.0pt] \displaystyle\frac{\partial^{2}\ln f}{\partial N_{+}\partial i}&\displaystyle\frac{\partial^{2}\ln f}{\partial N_{+}^{2}}\end{pmatrix}\right|_{(i_{\max},R_{+\max})}\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!=\frac{4}{N}\begin{pmatrix}\frac{-1}{a(1-a)}&\frac{1}{1-a}\\ \frac{1}{1-a}&\frac{-1}{a(1-a)}\\ \end{pmatrix}. (65)

Considering that the function is symmetric around the maximum R+maxR_{+\max}, we can approximate

VM2=20N/20min{N+,NA}f(i,N+)𝑑i𝑑N+.++f(i,N+)𝑑i𝑑N+=fmax2πdet(H)=fmaxπN2a(1a).\begin{split}V^{\mathrm{M}_{2}}&=2\int_{0}^{N/2}\int_{0}^{\min\{N_{+},N_{A}\}}f(i,N_{+})\ di\,dN_{+}.\\ &\approx\int_{-\infty}^{+\infty}\int_{-\infty}^{+\infty}f(i,N_{+})\ di\,dN_{+}\\ &=f_{\max}\frac{2\pi}{\sqrt{\mathrm{det}(-H)}}=f_{\max}\frac{\pi N}{2}\sqrt{a(1-a)}.\end{split} (66)

We obtain

lnVM2=N(lnwlnN+S(a)+ln2+1)ln(4πN).\begin{split}\ln V^{\mathrm{M}_{2}}&=N\left(\ln w-\ln N+S(a)+\ln 2+1\right)-\ln(4\pi N).\end{split} (67)

B.4 Morty 3

Finally, we consider the case of “ignorant” Morty, who not only cannot distiguish the two colors of particles, but also is not aware that these two colors of particles even exist. As a result, he also does not know how many particles of each type are in the system that he is observing. The only information he can obtain is how many particles are on the left and on the right side.

This means that when calculating the associated observational entropy, he will associate different volumes to these macrostates than their values given by the actual number of microstates that are contained within, given by Eq. (62).

From his point of view, the macrostate associated with having ii particles on the left and NiN-i particles on the right is given by

VM3=(wNA)(wNNA).V^{\mathrm{M}_{3}}=\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0,0,0}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@gray@stroke{0}\pgfsys@invoke{ }\pgfsys@color@gray@fill{0}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0,0,0}\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@moveto{-8.39357pt}{-7.39755pt}\pgfsys@curveto{-8.39357pt}{-6.14043pt}{-9.41266pt}{-5.12134pt}{-10.66978pt}{-5.12134pt}\pgfsys@curveto{-11.92691pt}{-5.12134pt}{-12.946pt}{-6.14043pt}{-12.946pt}{-7.39755pt}\pgfsys@curveto{-12.946pt}{-8.65468pt}{-11.92691pt}{-9.67377pt}{-10.66978pt}{-9.67377pt}\pgfsys@curveto{-9.41266pt}{-9.67377pt}{-8.39357pt}{-8.65468pt}{-8.39357pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{-14.22638pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(w\\ N_{A})}\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0,0,0}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@gray@stroke{0}\pgfsys@invoke{ }\pgfsys@color@gray@fill{0}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0,0,0}\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@moveto{-1.28038pt}{-7.39755pt}\pgfsys@curveto{-1.28038pt}{-6.14043pt}{-2.29947pt}{-5.12134pt}{-3.5566pt}{-5.12134pt}\pgfsys@curveto{-4.81372pt}{-5.12134pt}{-5.83281pt}{-6.14043pt}{-5.83281pt}{-7.39755pt}\pgfsys@curveto{-5.83281pt}{-8.65468pt}{-4.81372pt}{-9.67377pt}{-3.5566pt}{-9.67377pt}\pgfsys@curveto{-2.29947pt}{-9.67377pt}{-1.28038pt}{-8.65468pt}{-1.28038pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{0.0pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(w\\ N-N_{A})}. (68)

Using Eq. (52), we obtain

lnVM3=N(lnwlnN+S(a)+1)12ln((2πN)2a(1a))\begin{split}\ln V^{\mathrm{M}_{3}}&=N(\ln w-\ln N+S(a)+1)\\ &-\frac{1}{2}\ln((2\pi N)^{2}a(1-a))\end{split} (69)

B.5 Rick - non-symmetric case

Finally, let us consider a non-symmetric case of Rick, where the macrostate volume is given by

VR=(LAi)(LANAi)(LBN+i)(LBNN+(NAi))=(LApiN)(LA(api)N)(LB(rpi)N)(LB(1ra+pi)N)=(LAq1N)(LAq2N)(LBq3N)(LBq4N).\begin{split}&V^{\mathrm{R}}=\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.16796875,0.3125,0.640625}\definecolor[named]{pgfstrokecolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@color@rgb@stroke{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@moveto{-8.39357pt}{-7.39755pt}\pgfsys@curveto{-8.39357pt}{-6.14043pt}{-9.41266pt}{-5.12134pt}{-10.66978pt}{-5.12134pt}\pgfsys@curveto{-11.92691pt}{-5.12134pt}{-12.946pt}{-6.14043pt}{-12.946pt}{-7.39755pt}\pgfsys@curveto{-12.946pt}{-8.65468pt}{-11.92691pt}{-9.67377pt}{-10.66978pt}{-9.67377pt}\pgfsys@curveto{-9.41266pt}{-9.67377pt}{-8.39357pt}{-8.65468pt}{-8.39357pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{-14.22638pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(L_{A}\\ i)}\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.9375,0.23046875,0.15234375}\definecolor[named]{pgfstrokecolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@color@rgb@stroke{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@moveto{-8.39357pt}{-7.39755pt}\pgfsys@curveto{-8.39357pt}{-6.14043pt}{-9.41266pt}{-5.12134pt}{-10.66978pt}{-5.12134pt}\pgfsys@curveto{-11.92691pt}{-5.12134pt}{-12.946pt}{-6.14043pt}{-12.946pt}{-7.39755pt}\pgfsys@curveto{-12.946pt}{-8.65468pt}{-11.92691pt}{-9.67377pt}{-10.66978pt}{-9.67377pt}\pgfsys@curveto{-9.41266pt}{-9.67377pt}{-8.39357pt}{-8.65468pt}{-8.39357pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{-14.22638pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(L_{A}\\ N_{A}\!-\!i)}\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.16796875,0.3125,0.640625}\definecolor[named]{pgfstrokecolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@color@rgb@stroke{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@moveto{-1.28038pt}{-7.39755pt}\pgfsys@curveto{-1.28038pt}{-6.14043pt}{-2.29947pt}{-5.12134pt}{-3.5566pt}{-5.12134pt}\pgfsys@curveto{-4.81372pt}{-5.12134pt}{-5.83281pt}{-6.14043pt}{-5.83281pt}{-7.39755pt}\pgfsys@curveto{-5.83281pt}{-8.65468pt}{-4.81372pt}{-9.67377pt}{-3.5566pt}{-9.67377pt}\pgfsys@curveto{-2.29947pt}{-9.67377pt}{-1.28038pt}{-8.65468pt}{-1.28038pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{0.0pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(L_{B}\\ N_{+}\!-\!i)}\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.9375,0.23046875,0.15234375}\definecolor[named]{pgfstrokecolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@color@rgb@stroke{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@moveto{-1.28038pt}{-7.39755pt}\pgfsys@curveto{-1.28038pt}{-6.14043pt}{-2.29947pt}{-5.12134pt}{-3.5566pt}{-5.12134pt}\pgfsys@curveto{-4.81372pt}{-5.12134pt}{-5.83281pt}{-6.14043pt}{-5.83281pt}{-7.39755pt}\pgfsys@curveto{-5.83281pt}{-8.65468pt}{-4.81372pt}{-9.67377pt}{-3.5566pt}{-9.67377pt}\pgfsys@curveto{-2.29947pt}{-9.67377pt}{-1.28038pt}{-8.65468pt}{-1.28038pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{0.0pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(L_{B}\\ N\!-\!N_{+}\!-\!(N_{A}\!-\!i))}\\ &={\matrixquantity(L_{A}\\ p_{i}N)}{\matrixquantity(L_{A}\\ (a\!-\!p_{i})N)}{\matrixquantity(L_{B}\\ (r\!-\!p_{i})N)}{\matrixquantity(L_{B}\\ (1\!-\!r\!-\!a\!+\!p_{i})N)}\\ &={\matrixquantity(L_{A}\\ q_{1}N)}{\matrixquantity(L_{A}\\ q_{2}N)}{\matrixquantity(L_{B}\\ q_{3}N)}{\matrixquantity(L_{B}\\ q_{4}N)}.\end{split} (70)

Using Eq. (52), we find

lnVR=q1N(lnLAlnNlnq1+1)12ln2πNq1+q2N(lnLAlnNlnq2+1)12ln2πNq2+q3N(lnLBlnNlnq3+1)12ln2πNq3+q4N(lnLBlnNlnq4+1)12ln2πNq4=N(alnLA+(1a)lnLBlnN+S(q)+1)12ln((2πN)4j=14qjmissing).\begin{split}\ln V^{\mathrm{R}}&=q_{1}N(\ln L_{A}-\ln N-\ln q_{1}+1)-\frac{1}{2}\ln 2\pi Nq_{1}\\ &+q_{2}N(\ln L_{A}-\ln N-\ln q_{2}+1)-\frac{1}{2}\ln 2\pi Nq_{2}\\ &+q_{3}N(\ln L_{B}-\ln N-\ln q_{3}+1)-\frac{1}{2}\ln 2\pi Nq_{3}\\ &+q_{4}N(\ln L_{B}-\ln N-\ln q_{4}+1)-\frac{1}{2}\ln 2\pi Nq_{4}\\ &=N(a\ln L_{A}+(1\!-\!a)\ln L_{B}-\ln N+S(q)+1)\\ &-\frac{1}{2}\ln\Big((2\pi N)^{4}\prod_{j=1}^{4}q_{j}\Big{missing}).\end{split} (71)

If the gas is uniformly distributed with all the blue particles are on the left and all the red particles on the right, we have LA=aLL_{A}=aL, LB=(1a)LL_{B}=(1-a)L, q1=aq_{1}=a, q2=q3=0q_{2}=q_{3}=0 (which also means r=ar=a), and q4=1aq_{4}=1-a. This gives

lnVinitR=N(lnLlnN+1)12ln((2πN)2a(1a)missing).\begin{split}\ln V_{\mathrm{init}}^{\mathrm{R}}=N(\ln L-\ln N+1)-\frac{1}{2}\ln\Big((2\pi N)^{2}a(1-a)\Big{missing}).\end{split} (72)

When the particles are allow to mix they fill the space uniformly, which corresponds to

S𝒞Rfin=lnV=ln((LaN))(L(1a)N)=N(lnLlnN+S(a)+1)12ln((2πN)2a(1a)missing).\begin{split}S_{{\mathcal{C}}_{R}}^{\mathrm{fin}}&=\ln V=\ln{\matrixquantity(L\\ aN)}{\matrixquantity(L\\ (1-a)N)}\\ &=N(\ln L-\ln N+S(a)+1)\\ &-\frac{1}{2}\ln\Big((2\pi N)^{2}a(1-a)\Big{missing}).\end{split} (73)

The growth in entropy is the standard entropy of mixing,

ΔSR=S𝒞RfinS𝒞Rinit=lnVVinitR=NS(a).\begin{split}\Delta S^{\mathrm{R}}&=S_{{\mathcal{C}}_{R}}^{\mathrm{fin}}-S_{{\mathcal{C}}_{R}}^{\mathrm{init}}=\ln\frac{V}{V_{\mathrm{init}}^{R}}=NS(a).\end{split} (74)

B.6 Morty 1 - non-symmetric case

VM1=i=0min{N+,NA}(LAi)(LANAi)(LBN+i)(LBNN+(NAi))=:i=0min{N+,NA}f(i)fmax2πσ2\begin{split}&V^{\mathrm{M}_{1}}=\\ &\sum_{i=0}^{\min\{N_{+},N_{A}\}}\!\!\!\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.16796875,0.3125,0.640625}\definecolor[named]{pgfstrokecolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@color@rgb@stroke{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@moveto{-8.39357pt}{-7.39755pt}\pgfsys@curveto{-8.39357pt}{-6.14043pt}{-9.41266pt}{-5.12134pt}{-10.66978pt}{-5.12134pt}\pgfsys@curveto{-11.92691pt}{-5.12134pt}{-12.946pt}{-6.14043pt}{-12.946pt}{-7.39755pt}\pgfsys@curveto{-12.946pt}{-8.65468pt}{-11.92691pt}{-9.67377pt}{-10.66978pt}{-9.67377pt}\pgfsys@curveto{-9.41266pt}{-9.67377pt}{-8.39357pt}{-8.65468pt}{-8.39357pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{-14.22638pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(L_{A}\\ i)}\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.9375,0.23046875,0.15234375}\definecolor[named]{pgfstrokecolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@color@rgb@stroke{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@moveto{-8.39357pt}{-7.39755pt}\pgfsys@curveto{-8.39357pt}{-6.14043pt}{-9.41266pt}{-5.12134pt}{-10.66978pt}{-5.12134pt}\pgfsys@curveto{-11.92691pt}{-5.12134pt}{-12.946pt}{-6.14043pt}{-12.946pt}{-7.39755pt}\pgfsys@curveto{-12.946pt}{-8.65468pt}{-11.92691pt}{-9.67377pt}{-10.66978pt}{-9.67377pt}\pgfsys@curveto{-9.41266pt}{-9.67377pt}{-8.39357pt}{-8.65468pt}{-8.39357pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-10.66978pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{-14.22638pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(L_{A}\\ N_{A}\!-\!i)}\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.16796875,0.3125,0.640625}\definecolor[named]{pgfstrokecolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@color@rgb@stroke{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.16796875}{0.3125}{0.640625}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.16796875,0.3125,0.640625}\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@moveto{-1.28038pt}{-7.39755pt}\pgfsys@curveto{-1.28038pt}{-6.14043pt}{-2.29947pt}{-5.12134pt}{-3.5566pt}{-5.12134pt}\pgfsys@curveto{-4.81372pt}{-5.12134pt}{-5.83281pt}{-6.14043pt}{-5.83281pt}{-7.39755pt}\pgfsys@curveto{-5.83281pt}{-8.65468pt}{-4.81372pt}{-9.67377pt}{-3.5566pt}{-9.67377pt}\pgfsys@curveto{-2.29947pt}{-9.67377pt}{-1.28038pt}{-8.65468pt}{-1.28038pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{0.0pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(L_{B}\\ N_{+}\!-\!i)}\underset{\leavevmode\hbox to14.63pt{\vbox to7.51pt{\pgfpicture\makeatletter\hbox{\qquad\lower-11.58092pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@invoke{ }\pgfsys@endscope\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{}} {}{{}}{}{{{}}{}{}{}{}{}{}{}{}}\pgfsys@beginscope\pgfsys@invoke{ }\color[rgb]{0.9375,0.23046875,0.15234375}\definecolor[named]{pgfstrokecolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@color@rgb@stroke{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0.9375}{0.23046875}{0.15234375}\pgfsys@invoke{ }\definecolor[named]{pgffillcolor}{rgb}{0.9375,0.23046875,0.15234375}\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@moveto{-1.28038pt}{-7.39755pt}\pgfsys@curveto{-1.28038pt}{-6.14043pt}{-2.29947pt}{-5.12134pt}{-3.5566pt}{-5.12134pt}\pgfsys@curveto{-4.81372pt}{-5.12134pt}{-5.83281pt}{-6.14043pt}{-5.83281pt}{-7.39755pt}\pgfsys@curveto{-5.83281pt}{-8.65468pt}{-4.81372pt}{-9.67377pt}{-3.5566pt}{-9.67377pt}\pgfsys@curveto{-2.29947pt}{-9.67377pt}{-1.28038pt}{-8.65468pt}{-1.28038pt}{-7.39755pt}\pgfsys@closepath\pgfsys@moveto{-3.5566pt}{-7.39755pt}\pgfsys@fill\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{{}}{} {}{}{}\pgfsys@moveto{-14.22638pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-11.38092pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{-7.11319pt}{-11.38092pt}\pgfsys@lineto{-7.11319pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } {}{{}}{} {}{}{}\pgfsys@moveto{0.0pt}{-11.38092pt}\pgfsys@lineto{0.0pt}{-4.26773pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{ }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{ }\pgfsys@endscope\hss}}\endpgfpicture}}}{\matrixquantity(L_{B}\\ N\!-\!N_{+}\!-\!(N_{A}\!-\!i))}\\ &=:\sum_{i=0}^{\min\{N_{+},N_{A}\}}f(i)\approx f_{\max}\sqrt{2\pi\sigma^{2}}\end{split} (75)

where fmaxf_{\max} is given again by i=N+NAN=arNi=\frac{N_{+}N_{A}}{N}=arN, and the variance is given again by

σ2=(d2lnfdi2|i=imax)1=Na(1a)r(1r).\sigma^{2}=-\left(\left.\frac{d^{2}\ln f}{di^{2}}\right|_{i=i_{\max}}\right)^{-1}=Na(1-a)r(1-r). (76)

We obtain fmaxf_{\max} by inserting pi=arp_{i}=ar into Eq. (71) and derive

lnVM1=N(alnLA+(1a)lnLBlnN+S(a)+S(r)+1)32ln(2πN).\begin{split}\ln V^{\mathrm{M}_{1}}&=N\big{(}a\ln L_{A}+(1\!-\!a)\ln L_{B}-\ln N\\ &+S(a)+S(r)+1\big{)}-\frac{3}{2}\ln(2\pi N).\end{split} (77)

For the same initial and final state as considered in Rick’s non-symmetric case, we have

lnVinitM1=N(lnLlnN+S(a)+1)32ln(2πN),\begin{split}\ln V_{\mathrm{init}}^{\mathrm{M}_{1}}&=N\big{(}\ln L-\ln N+S(a)+1\big{)}-\frac{3}{2}\ln(2\pi N),\end{split} (78)

while the final entropy is the same as Rick’s, Eq. (73). Thus, the entropy stays almost constant during the time evolution, with only a logarithmic increase coming from the increased uncertainty about the exact number of particles on each side,

ΔSM1=lnVVinitM1=32ln(2πNmissing)12ln((2πN)2a(1a)missing)=12ln2πNa(1a).\begin{split}\Delta S^{\mathrm{M}_{1}}&=\ln\frac{V}{V_{\mathrm{init}}^{\mathrm{M}_{1}}}=\frac{3}{2}\ln\Big(2\pi N\Big{missing})-\frac{1}{2}\ln\Big((2\pi N)^{2}a(1-a)\Big{missing})\\ &=\frac{1}{2}\ln\frac{2\pi N}{a(1-a)}.\end{split} (79)

Analysis for the non-symmetric case for Morty 2 and 3 proceeds analogously, with the same results as the symmetric case with a substitution lnwalnLA+(1a)lnLB\ln w\rightarrow a\ln L_{A}+(1\!-\!a)\ln L_{B}.

Appendix C Morty 3 work extraction

In this section, we investigate what it means when Morty 3 attempts to extract work using his incorrect judgment of the situation.

For that, we will have to assume some structure about what he thinks he is doing versus what he is actually doing, as his understanding of the situation is limited.

Let us summarize what we know about Morty 3. He believes there is only one type of particle in the system—he sees only black particles and is unaware that the system may in fact consist of a mixture of red and blue particles. A similar situation would be that of a physicist working with spin-1/2 particles before the discovery of spin.

As a result, he not only underestimates the system’s entropy—judging it by Eq. (17)—but also misidentifies the Hilbert space the particles occupy. He assigns a state corresponding to what he believes to be black particles, whereas in reality, the system is in a different state composed of red and blue particles. Consequently, when he performs unitary operations to extract work, he believes they act only on black particles, while in fact they affect a mixture of red and blue particles.

Given this, we construct a model of the actual system that aligns with Morty 3’s incomplete description of reality, based on specific assumptions. Additionally, we must identify the underlying transformation that is truly implemented when he believes he is performing an operation. Our approach must be consistent—the actual transformation should reproduce the outcome he intends to achieve.

C.1 Consistency, assumptions, and the resulting structure

Mathematically, we define consistency as follows: for all measurements accessible to him, the outcome probabilities before and after his unitary operation must coincide with those generated by the actual state under our model transformation (see Eqs. (83) and  (84) below).

Consider that ρ{\rho} and UU represent the actual state and transformation, while ρM{\rho}_{M} and UMU_{M}, respectively, denote his corresponding perception of the same. Also, consider 𝒞={Pi}{\mathcal{C}}=\{{P}_{i}\} the actual measurement, while 𝒞M={PiM}{\mathcal{C}}_{M}=\{{P}_{i}^{M}\} is the measurement he thinks he performs.

We begin by assuming the consistency of measurements (a1)(a1). This assumption consists of two parts: First, we assume that Morty 3 cannot perform a measurement that distinguishes the colors of the particles in any way. In other words, the set of admissible measurements is restricted to those that cannot distinguish between the blue and the red particles. Second, we assume a physically motivated correspondence between the perceived and actual measurements:

(a1)\displaystyle(a1) 𝒞M𝒞.\displaystyle{\mathcal{C}}_{M}\longrightarrow{\mathcal{C}}. (80)

For example, consider a two-site chain containing a single red particle, which appears black to Morty 3. From his perspective, the measurement of the particle’s position—whether it is on the left or right site—is given by

𝒞M={|1010|,|0101|}.{\mathcal{C}}_{M}=\{\outerproduct{10}{10},\outerproduct{01}{01}\}. (81)

In reality, the particle could be either red or blue, corresponding to an extended Hilbert space. Thus, the measurement Morty performs effectively projects onto the subspaces where some particle—regardless of its internal type—is located on the left or right site.

𝒞={|+0+0|+|00|,|0+0+|+|00|},\begin{split}{\mathcal{C}}=\{&\outerproduct{+0}{+0}+\outerproduct{-0}{-0},\\ &\outerproduct{0+}{0+}+\outerproduct{0-}{0-}\},\end{split} (82)

This measurement can be expressed using the formalism of Morty 2, who has the correct understanding of the Hilbert space. It illustrates that both internal states, |+0\ket{+0} and |0\ket{-0}, which are indistinguishable to Morty 3, will be registered as the same outcome by him. This can be seen as follows: denoting the projector onto the ”left” position as PL=|+0+0|+|00|{P}_{L}=\outerproduct{+0}{+0}+\outerproduct{-0}{-0} (seen on the left), we find that pL=+0|PL|+0=0|PL|0=1p_{L}=\bra{+0}{P}_{L}\ket{+0}=\bra{-0}{P}_{L}\ket{-0}=1. We will formalize this assumption in more detail later and extend it to arbitrary projective measurements.

Second, we assume consistency of density matrices (assumption (a2)(a2).). For any admissible perceived measurement 𝒞M{\mathcal{C}}_{M}, with corresponding actual measurement 𝒞{\mathcal{C}}, the actual density matrix ρ{\rho} must reproduce the same outcome distribution as the perceived density matrix ρM{\rho}_{M},

(a2)\displaystyle(a2) ρ𝒞\displaystyle{\rho}\overset{{\mathcal{C}}}{\longrightarrow}\ pi𝒞MρM,\displaystyle p_{i}\overset{{\mathcal{C}}_{M}}{\longleftarrow}{\rho}_{M}, (83)

where pi=tr[ρPi]=tr[ρMPiM]p_{i}=\tr[{\rho}{P}_{i}]=\tr[{\rho}_{M}{P}_{i}^{M}].

Third, we assume consistency of unitary operations (assumption (a3)(a3)). If Morty 3 believes he is applying a unitary UMU_{M}, then the actual unitary UU applied must produce a state that is indistinguishable—under any measurement that cannot resolve particle colors—from the state he believes he has prepared,

(a3)\displaystyle(a3) UρU𝒞\displaystyle U{\rho}U^{\dagger}\overset{{\mathcal{C}}}{\longrightarrow}\ qi𝒞MUMρMUM,\displaystyle q_{i}\overset{{\mathcal{C}}_{M}}{\longleftarrow}U_{M}{\rho}_{M}U_{M}^{\dagger}, (84)

where qi=tr[UρUPi]=tr[UMρMUMPiM]q_{i}=\tr[U{\rho}U^{\dagger}{P}_{i}]=\tr[U_{M}{\rho}_{M}U_{M}^{\dagger}{P}_{i}^{M}].

The next set of assumptions concerns the admissible states, allowed evolutions, and applicable unitary operators, in order to uphold superselection rules between different particle types. We restrict ourselves to the situation where there is a fixed and known total number of particles NN, as assumed in the main text, implying the actual Hilbert space (which is the Hilbert space of Morty 2, who cannot distinguish the colors, but knows that two colors exist),

(a4)\displaystyle(a4) M2\displaystyle{\mathcal{H}}\equiv{\mathcal{H}}^{\mathrm{M}_{2}} =N+=0N+N+NN+.\displaystyle=\bigoplus_{N_{+}=0}^{N}{\mathcal{H}}_{+}^{N_{+}}\otimes{\mathcal{H}}_{-}^{N-N_{+}}. (85)

(Recall definition (14).)

Further, we assume that the actual density matrix does not admit superpositions of states that involve differing number of particles of each type, such as (|0+|+0)/2(\ket{0-}+\ket{+0})/\sqrt{2} 333We denote |0\ket{0-} to express that in the first position there are no particles, and in the second position there is a red particle and so on., which we can write as

(a5)\displaystyle(a5) ρ=N+=0NλN+ρN+.\displaystyle{\rho}=\bigoplus_{N_{+}=0}^{N}\lambda_{N_{+}}{\rho}_{N_{+}}. (86)

where N+=0NλN+=1\sum_{N_{+}=0}^{N}\lambda_{N_{+}}=1, where ρN+{\rho}_{N_{+}} is a density matrix—linear operator—acting on +N+NN+{\mathcal{H}}_{+}^{N_{+}}\otimes{\mathcal{H}}_{-}^{N-N_{+}}. In other words, it can only admit classical mixtures between these subspaces, thereby forming a super-selection rule.

We further assume that the unitary that he performs cannot change the color of the particles, meaning that it acts independently on both the red and the blue particle Hilbert spaces, which is mathematically described as factorizing

(a6)\displaystyle(a6) U=N+=0NUN,N+.\displaystyle U=\bigoplus_{N_{+}=0}^{N}U^{N,N_{+}}. (87)

where every unitary UN,N+U^{N,N_{+}} acts only at the subspace +N+NN+{\mathcal{H}}_{+}^{N_{+}}\otimes{\mathcal{H}}_{-}^{N-N_{+}}. This also means that UU is block-diagonal in the full Hilbert space {\mathcal{H}}. This ensures that the super-selection cannot be broken by Morty 3 applying a unitary. Intuitively, this ensures that Morty 3, by virtue of his ignorance, does not wield more operational ability than Rick. We constructed Rick’s Hilbert space to be constrained by the fixed number of red and blue particles, which implicitly restricts the set of extraction unitaries he can perform. However, because Morty 3’s ignorance grants access to a larger Hilbert space, an unconstrained unitary could, in principle, swap the particle colors. To preserve the intended hierarchy between the agents, we ensure that Rick remains the more capable one.

Furthermore, let HM{H}_{M} denote the Hamiltonian used by Morty 3 to describe the system’s evolution and to assign energy to its states. This Hamiltonian does not account for red and blue particles, as Morty 3 is unaware of their existence. In contrast, let H{H} denote the actual Hamiltonian, which includes a full description of particles of both colors. We assume that the total (actual) Hamiltonian conserves the number of particles of each color,

(a7)\displaystyle(a7) H=N+=0NHN,N+.\displaystyle{H}=\bigoplus_{N_{+}=0}^{N}H^{N,N_{+}}. (88)

This means that the blocks defined by NN and N+N_{+} act and evolve independently of each other. In practice, (a6)(a6) and (a7)(a7) impose conservation of the probability of finding the state in each subspace +N+NN+{\mathcal{H}}_{+}^{N_{+}}\otimes{\mathcal{H}}_{-}^{N-N_{+}}.

Finally, we assume that (a8)(a8) the actual Hamiltonian is independent of particle color and consistent with Morty 3’s description. For example, if |011+|110\ket{011}+\ket{110} is an energy eigenvector of HM{H}_{M}, then |0++|+0\ket{0+-}+\ket{-+0} must also be an eigenvector of H{H} with the same energy. In other words, the Hamiltonian depends only on the configuration—that is, the arrangement of particles on the lattice—not on their color. This ensures that Morty can correctly determine the energy of a state by inspecting its configuration alone. It is further assumed that Morty 3 possesses a Hamiltonian which, although incorrect in structure, assigns correct energies to all relevant states.

For this, it will be useful to define a linear mapping 𝒦{\mathcal{K}} between the Morty 2 and Morty 3 Hilbert spaces, {\mathcal{H}} and MM3{\mathcal{H}}^{M}\equiv{\mathcal{H}}^{\mathrm{M}_{3}}, defined by its action on the basis vectors which erases the information about the type of particle, i.e., for example

𝒦(|00++0)=|001101{\mathcal{K}}(\ket{00++0-})=\ket{001101} (89)

We call the vector on the right-hand side configuration and on the left-hand side inner configuration. The linearity means that for two inner configuration states c~1,c~1\tilde{c}_{1},\tilde{c}_{1}\in{\mathcal{H}}, and a complex number α\alpha we have

𝒦(αc~1+c~2)=α𝒦(c~1)+𝒦(c~1).{\mathcal{K}}(\alpha\tilde{c}_{1}+\tilde{c}_{2})=\alpha{\mathcal{K}}(\tilde{c}_{1})+{\mathcal{K}}(\tilde{c}_{1}). (90)

Given a vector cMc\in{\mathcal{H}}^{M} , we define a one-dimensional Hilbert space cM=span{c}{\mathcal{H}}_{c}^{M}=\mathrm{span}\{c\}. We define the pre-image of this Hilbert space as,

c𝒦1(cM)={c~|𝒦(c~)cM}.{\mathcal{H}}_{c}\equiv{\mathcal{K}}^{-1}({\mathcal{H}}_{c}^{M})=\{\tilde{c}\in{\mathcal{H}}|{\mathcal{K}}(\tilde{c})\in{\mathcal{H}}_{c}^{M}\}. (91)

For example, for c1=(|01+|10)/2c_{1}=(\ket{01}+\ket{10})/\sqrt{2}, we have

c1=span{(|0++|+0)/2,(|0++|0)/2,(|0+|+0)/2,(|0+|0)/2}.\begin{split}{\mathcal{H}}_{c_{1}}=\mathrm{span}\{&(\ket{0+}+\ket{+0})/\sqrt{2},(\ket{0+}+\ket{-0})/\sqrt{2},\\ &(\ket{0-}+\ket{+0})/\sqrt{2},(\ket{0-}+\ket{-0})/\sqrt{2}\}.\end{split} (92)

For c2=(|01|10)/2c_{2}=(\ket{01}-\ket{10})/\sqrt{2}, we have

c2=span{(|0+|+0)/2,(|0+|0)/2,(|0|+0)/2,(|0|0)/2}.\begin{split}{\mathcal{H}}_{c_{2}}=\mathrm{span}\{&(\ket{0+}-\ket{+0})/\sqrt{2},(\ket{0+}-\ket{-0})/\sqrt{2},\\ &(\ket{0-}-\ket{+0})/\sqrt{2},(\ket{0-}-\ket{-0})/\sqrt{2}\}.\end{split} (93)

c{\mathcal{H}}_{c} is a Hilbert space. To prove that it is a vector space, assume a linear combination |ψ=iαi|ϕi\ket{\psi}=\sum_{i}\alpha_{i}\ket{\phi_{i}}, where |ϕic\ket{\phi_{i}}\in{\mathcal{H}}_{c}. Then from the linearity we have,

𝒦|ψ=𝒦(iαi|ϕi)=iαi𝒦(|ϕi)=ciαiβicM.\begin{split}{\mathcal{K}}\ket{\psi}&={\mathcal{K}}\big{(}\sum_{i}\alpha_{i}\ket{\phi_{i}}\big{)}=\sum_{i}\alpha_{i}{\mathcal{K}}(\ket{\phi_{i}})\\ &=c\sum_{i}\alpha_{i}\beta_{i}\in{\mathcal{H}}_{c}^{M}.\end{split} (94)

Additionally, it is equipped with an inner product induced from {\mathcal{H}}, and it is finite. Thus, it is complete.

Note, however, that surprisingly also (|+0+|0)c1(\ket{+0}+\ket{-0})\in{\mathcal{H}}_{c_{1}} and (|0+|0)c2(\ket{0+}-\ket{0-})\in{\mathcal{H}}_{c_{2}} are in the kernel, corresponding to zero vector in c1M{\mathcal{H}}_{c_{1}}^{M} and c2M{\mathcal{H}}_{c_{2}}^{M}, respectively. There is only a finite number of such non-trivial vectors in the kernel, which is a set of measure zero, since it requires a very specific combination of parameters so that iαiβi=0\sum_{i}\alpha_{i}\beta_{i}=0. This is because the vectors are always characterized by a combination that maps onto the same vector, such as |+|+\ket{+-}-\ket{-+} or similar, which maps onto |11|11=0\ket{11}-\ket{11}=0. First, physical states always contain some noise, so this pathological situation should not occur in practice. Second, the situation where the state is 0 corresponds to not having a system; thus, we will ignore it. The second observation is that although c1c_{1} and c2c_{2} are orthogonal, not every two vectors in |ψ1c1\ket{\psi_{1}}\in{\mathcal{H}}_{c_{1}} and |ψ1c2\ket{\psi_{1}}\in{\mathcal{H}}_{c_{2}} are. E.g., (|0++|0)/2(\ket{0+}+\ket{-0})/\sqrt{2} and (|0|0)/2(\ket{0-}-\ket{-0})/\sqrt{2} have squared overlap 1/21/2. Even if we restrict ourselves to a subspace +N+NN+{\mathcal{H}}_{+}^{N_{+}}\otimes{\mathcal{H}}_{-}^{N-N_{+}}, still, one can find two vectors created as preimages from orthogonal vectors which have non-zero overlap. If the preimages were, however, created in different subspaces +N+NN+{\mathcal{H}}_{+}^{N_{+}}\otimes{\mathcal{H}}_{-}^{N-N_{+}}, they are orthogonal by default.

Consider an orthonormal basis {c}\{c\}, so that M=span{c}{\mathcal{H}}^{M}=\mathrm{span}\{c\}, following the argument after Eq. . Then the subspaces created from the pre-images combine into the entire Hilbert space,

=cc.{\mathcal{H}}=\bigoplus_{c}{\mathcal{H}}_{c}. (95)

As we mentioned above, states from different c{\mathcal{H}}_{c} are not orthogonal. However, they are still linearly independent. To prove this, one can assume |ψc\ket{\psi}\in{\mathcal{H}}_{c^{\prime}}, can be written as a linear combination of vectors from c{\mathcal{H}}_{c}, ccc\neq c^{\prime}. However, if that is possible, then 𝒦(|ψ)=αc{\mathcal{K}}(\ket{\psi})=\alpha c for some number α\alpha, meaning that |ψc\ket{\psi}\in{\mathcal{H}}_{c}, which is in contradiction with the assumption that ccc\neq c^{\prime}. This means that the direct sum is well-defined. Additionally, one can always find an orthonormal basis of each c{\mathcal{H}}_{c} so that each vectors from there are orthogonal to each other. This, however, doesn’t necessarily imply that they will be orthogonal to vectors from another subspace, c{\mathcal{H}}_{c^{\prime}}, unless there is some additional structure involved (such as when different subspaces are eigensubspaces of a Hermitian operator).

Assumption (a8)(a8) says that the Hamiltonian HN,N+H^{N,N_{+}} does not depend on the inner configuration. Consider that Morty 3 has a Hamiltonian

HM=EE|EE|,H_{M}=\sum_{E}E\outerproduct{E}{E}, (96)

which judges the energy correctly. Energy eigenstates |E\ket{E} are a superposition of different configurations. This means that all eigenstates of the actual Hamiltonian that map onto the same eigenstate of Morty 3’s Hamiltonian must correspond to the same energy. In other words, combining this with assumption (a7)(a7), Eq. (88), the actual Hamiltonian is given by

HN,N+=EEPEN,N+,\displaystyle H^{N,N_{+}}=\sum_{E}E{P}_{E}^{N,N_{+}}, (97)

where PEN,N+{P}_{E}^{N,N_{+}} is a projector onto a subspace that is the overlap of two vector spaces

|E+N+NN+.{\mathcal{H}}_{\ket{E}}\cap{\mathcal{H}}_{+}^{N_{+}}\otimes{\mathcal{H}}_{-}^{N-N_{+}}. (98)

The overlap means that superselection rules still apply during evolution, i.e., superpositions of states that involve differing numbers of particles of each type, such as (|0+|+0)/2(\ket{0-}+\ket{+0})/\sqrt{2}, are not created while the system evolves.

Similarly, assuming that Morty thinks that he applies a unitary with spectral decomposition

UM=uu|uu|,U_{M}=\sum_{u}u\outerproduct{u}{u}, (99)

the unitary he actually applies will be given by

UN,N+=uuUuN,N+,U^{N,N_{+}}=\sum_{u}u\,U_{u}^{N,N_{+}}, (100)

inserted into (87). UuN,N+U_{u}^{N,N_{+}} is an operator acting on

+N+NN+,{\mathcal{H}}_{+}^{N_{+}}\otimes{\mathcal{H}}_{-}^{N-N_{+}}, (101)

which acts as a unitary on its subspace |u{\mathcal{H}}_{\ket{u}}, while zero on its complement,

UuN,N+UuN,N+0|u.U_{u}^{N,N_{+}}\equiv U_{u}^{N,N_{+}}\oplus 0_{{\mathcal{H}}_{\ket{u}}^{\perp}}. (102)

The reasoning behind this is that only those internal states which are indistinguishable to Morty 2 can be freely mapped onto each other via unitary transformations—such changes will go unnoticed by him. From his perspective, the evolution appears exactly as if the unitary UMU_{M} had been applied. If this were not the case, he could, in principle, perform a measurement revealing that the unitary did not act as expected, violating assumption (a3)(a3) and contradicting Eq. (84). The structure of the unitary operator (100) reflects the structure of (99): the same way the rank-1 projector |uu|\outerproduct{u}{u} acts as an identity unitary operator on the one-dimensional Hilbert space spanned by itself, it acts as a zero operator on all other vectors |u|u\ket{u^{\prime}}\perp\ket{u}. The inner structure (102) reflects that.

This definition leads to,

UuN,N+UuN,N+={I|u0|u,u=u,0,uu.U_{u}^{N,N_{+}}U_{u^{\prime}}^{N,N_{+}{\dagger}}=\begin{cases}{I}_{{\mathcal{H}}_{\ket{u}}}\oplus 0_{{\mathcal{H}}_{\ket{u}}^{\perp}},&u=u^{\prime},\\ 0,&u\neq u^{\prime}.\end{cases} (103)

Using this, we can prove unitarity,

UN,N+UN,N+=uuUuN,N+uuUuN,N+=u,uδuuuuI=u|u|2I|u0|u=uI|u0|u=I,\begin{split}U^{N,N_{+}}U^{N,N_{+}{\dagger}}&=\sum_{u}u\,U_{u}^{N,N_{+}}\sum_{u^{\prime}}u^{\prime}\,U_{u^{\prime}}^{N,N_{+}}\\ &=\sum_{u,u^{\prime}}\delta_{uu^{\prime}}uu^{\prime}{I}=\sum_{u}\absolutevalue{u}^{2}{I}_{{\mathcal{H}}_{\ket{u}}}\oplus 0_{{\mathcal{H}}_{\ket{u}}^{\perp}}\\ &=\sum_{u}{I}_{{\mathcal{H}}_{\ket{u}}}\oplus 0_{{\mathcal{H}}_{\ket{u}}^{\perp}}={I},\end{split} (104)

because |u|2=1\absolutevalue{u}^{2}=1. The meaning of UuN,N+U_{u}^{N,N_{+}} is that it can mix all vectors within the subspace in which Morty 3 is unable to observe any changes, while outwardly appearing to act exactly as UMU_{M} does.

Finally, if the density matrix describing the system in the Morty 3 framework is

ρM=rr|rr|,{\rho}_{M}=\sum_{r}r\outerproduct{r}{r}, (105)

then the actual density matrix is given by inserting

ρN+=rrρrN,N+{\rho}_{N_{+}}=\sum_{r}r{\rho}_{r}^{N,N_{+}} (106)

into Eq. (86). ρrN,N+{\rho}_{r}^{N,N_{+}} is a linear operator acting on

+N+NN+.{\mathcal{H}}_{+}^{N_{+}}\otimes{\mathcal{H}}_{-}^{N-N_{+}}. (107)

which acts as density matrix on its subspace |r{\mathcal{H}}_{\ket{r}} while zero on its complement,

ρrN,N+ρrN,N+0|r.{\rho}_{r}^{N,N_{+}}\equiv{\rho}_{r}^{N,N_{+}}\oplus 0_{{\mathcal{H}}_{\ket{r}}^{\perp}}. (108)

This construction ensures that outwardly the ρ{\rho} appears exactly as ρM{\rho}_{M}, since all the integral degrees of freedom (described by ρr{\rho}_{r}) map onto the same external (described by |rr|\outerproduct{r}{r}) visible to Morty 3, while the external probabilities, rr, which are observable by Morty 3, stay the same. If it had some other structure, it would be distinguishable by some measurement, and thus in contradition with (a2)(a2), Eq. (83).

Finally, using the map 𝒦{\mathcal{K}} we can rephrase assumption (a1)(a1) and derive the prescription for the measurement Morty 3 actually performs. Given Morty 3’s 𝒞M={PiM}{\mathcal{C}}_{M}=\{{P}_{i}^{M}\} description of the measurement he thinks he performs, the actual measurement he performs is

(a1)\displaystyle(a1) 𝒞={Pi}.\displaystyle{\mathcal{C}}=\{{P}_{i}\}. (109)

There, Pi{P}_{i} is the projector onto

PiM:=k|ψikM,{\mathcal{H}}_{{P}_{i}^{M}}:=\bigoplus_{k}{\mathcal{H}}_{\ket{\psi_{ik}^{M}}}, (110)

where we use the spectral decomposition into orthonormal projectors, PiM=k|ψikMψikM|{P}_{i}^{M}=\sum_{k}\outerproduct{\psi_{ik}^{M}}{\psi_{ik}^{M}}. Note that |ψik|ψikM\ket{\psi_{ik}}\in{\mathcal{H}}_{\ket{\psi_{ik}^{M}}} and |ψik|ψikM\ket{\psi_{ik^{\prime}}}\in{\mathcal{H}}_{\ket{\psi_{ik^{\prime}}^{M}}} are not necessarily orthogonal for kkk\neq k^{\prime}, however, one can always employ an orthogonalization procedure from which Pi{P}_{i} is explicitly formed. Further, subspaces PiM{\mathcal{H}}_{{P}_{i}^{M}} and PiM{\mathcal{H}}_{{P}_{i^{\prime}}^{M}} are not necessarily orthogonal for iii\neq i^{\prime}, but they will be if they correspond to different eigenvalues of some normal operator, such as a Hamiltonian or a unitary operator. (We use this several times below.) We could additionally assume that Morty 3 is unable to project onto superpositions of inner configurations—analogous to our assumption on the Hamiltonian in Eq. (88)—to prevent the creation of forbidden superpositions, such as |++|\ket{+}+\ket{-}, through measurement. However, this assumption is unnecessary to prove the statement we are trying to reach.

One can verify that this definition indeed corresponds to our earlier, less formal prescription given in Eq. (82). As an additional example, consider a rank-2 projector corresponding to a measurement that detects a total particle number of one.

P1M=|1010|+|0101|.{P}_{1}^{M}=\outerproduct{10}{10}+\outerproduct{01}{01}. (111)

The actual projector corresponding to this is

P1=|+0+0|+|0+0+|+|00|+|00|.\begin{split}{P}_{1}&=\outerproduct{+0}{+0}+\outerproduct{0+}{0+}+\outerproduct{-0}{-0}+\outerproduct{0-}{0-}.\end{split} (112)

Finally, consider a projector which projects onto a superposition of number states, e.g.,

PiM=|ψMψM||ψM=a1|011+a2|110.P_{i}^{M}=\outerproduct{\psi^{M}}{\psi^{M}}\quad\ket{\psi^{M}}=a_{1}\ket{011}+a_{2}\ket{110}. (113)

When restricting ourselves to the assumptions (a4a7)(a4-a7) which forbids superpositions such as a1|0+a2|++0a_{1}\ket{0--}+a_{2}\ket{++0}, the corresponding Hilbert space is,

|ψM=123,{\mathcal{H}}_{\ket{\psi^{M}}}={\mathcal{H}}_{1}\oplus{\mathcal{H}}_{2}\oplus{\mathcal{H}}_{3}, (114)

where

1=span{a1|0+++a2|++0}2=span{a1|0+a2|0}3=span{a1|0++a2|+0,a1|0++a2|+0,a1|0++a2|+0,a1|0++a2|+0}.\begin{split}{\mathcal{H}}_{1}=\mathrm{span}\{&a_{1}\ket{0++}+a_{2}\ket{++0}\}\\ {\mathcal{H}}_{2}=\mathrm{span}\{&a_{1}\ket{0--}+a_{2}\ket{--0}\}\\ {\mathcal{H}}_{3}=\mathrm{span}\{&a_{1}\ket{0+-}+a_{2}\ket{+-0},\\ &a_{1}\ket{0+-}+a_{2}\ket{-+0},\\ &a_{1}\ket{0-+}+a_{2}\ket{+-0},\\ &a_{1}\ket{0-+}+a_{2}\ket{-+0}\}.\\ \end{split} (115)

and allowed states follow the superselection rule ρ=iPiρPi{\rho}=\sum_{i}{P}_{i}{\rho}{P}_{i}, where Pi{P}_{i} projects onto i{\mathcal{H}}_{i}. As mentioned above, this is an example of where not all of these basis vectors are orthogonal, but can be made orthogonal by orthogonalization procedure.

As a side note, this can be further generalized to any observable: consider Morty 3 has an observable, OM=ooPoMO^{M}=\sum_{o}o{P}_{o}^{M}. The corresponding actual observable is O=ooPoO=\sum_{o}o{P}_{o}, where Po{P}_{o} is the projector onto PoM{\mathcal{H}}_{{P}_{o}^{M}}. The Hamiltonian, Eq. (97), already follows this pattern, with an additional assumption of Eq. (88).

C.2 Extracted work by Morty 3

Now we have most of the ingredients for deriving Morty 3’s actual extracted work.

We are interested in whether Morty 3’s estimate of the extracted work matches the actual extracted work, given the set of consistency assumptions (a1(a1a3)a3), the superselection rule stating that no superposition between different types of particles is allowed, and the assumption that such superselection sectors are preserved under both the dynamics and applied unitaries (a4(a4a7)a7). Finally, we assume that Morty 3 can correctly judge the system’s energy (a8)(a8), meaning that the Hamiltonian depends only on the configuration and not on the colors of the particles.

Morty 3’s estimate of work, for a single copy of the state of the system (i.e., without simultaneous estimation from multiple copies) is [72],

W𝒞M=U~Mtr[HM(ρUMextU~MρU~MUMext)]μ(U~M)=tr[HM(ρMπMcg)],\begin{split}W_{{\mathcal{C}}_{M}}&=\int_{\tilde{U}_{M}}\tr[{H}_{M}(\rho-U_{M}^{\mathrm{ext}}\tilde{U}_{M}{\rho}\tilde{U}_{M}^{\dagger}U_{M}^{\mathrm{ext}{\dagger}})]\mu(\tilde{U}_{M})\\ &=\tr[{H}_{M}({\rho}_{M}-\pi_{M}^{\mathrm{cg}})],\end{split} (116)

where

πMcg=UMextρMcgUMext,\pi_{M}^{\mathrm{cg}}=U_{M}^{\mathrm{ext}}{\rho}_{M}^{\mathrm{cg}}U_{M}^{\mathrm{ext}{\dagger}}, (117)

is a passive state to ρMcg=ipiViMPiM=U~MρMU~Mμ(U~M){\rho}_{M}^{{\mathrm{cg}}}=\sum_{i}\frac{p_{i}}{V_{i}^{M}}{P}_{i}^{M}=\int\tilde{U}_{M}{\rho}_{M}\tilde{U}_{M}^{\dagger}\mu(\tilde{U}_{M}), where ΠiM\Pi_{i}^{M}, where 𝒞M={PiM}{\mathcal{C}}_{M}=\{{P}_{i}^{M}\}, ViM=tr[PiM]V_{i}^{M}=\tr[{P}_{i}^{M}]. U~\tilde{U} is a random unitary randomizing vectors in each subspace i{\mathcal{H}}_{i}, onto which PiM{P}_{i}^{M} projects, picked with the Haar measure, and UMext.U_{M}^{\mathrm{ext.}} is the extraction unitary which reorders the eigenvalues from the largest to the smallest and associates them to the energy eigenvectors from the lowest to the highest energy. UM=UMextU~MU_{M}=U_{M}^{\mathrm{ext}}\tilde{U}_{M} is the full unitary that Morty 3 thinks he applies on the system, consisting both of the random and the extraction unitary.

In reality, however, he performs a measurement in the basis 𝒞={Pi}{\mathcal{C}}=\{{P}_{i}\}, where projectors Pi{P}_{i} include all the different colors that lead to the same configurations as vectors is Pi{P}_{i}. This is mathematically described in Eq. (109), denoted 𝒞M𝒞{\mathcal{C}}_{M}\rightarrow{\mathcal{C}}. His actual state of the system is ρMρ{\rho}_{M}\rightarrow{\rho}, Eq. (106). The random unitary he actually applies is UMUU_{M}\rightarrow U, Eq. (100), which consists both of the random and the extraction unitary, U=UextU~U=U^{\mathrm{ext}}\tilde{U}, U~MU~\tilde{U}_{M}\rightarrow\tilde{U}, UMextUextU_{M}^{\mathrm{ext}}\rightarrow U^{\mathrm{ext}}. The actual Hamiltonian is HMH{H}_{M}\rightarrow{H}, Eq. (97).

The actual energy he will extract on average is then

W𝒞=U~tr[H(ρUρU)]ν(U~),\begin{split}&W_{{\mathcal{C}}}=\int_{\tilde{U}}\tr[{H}({\rho}-U{\rho}U^{\dagger})]\nu(\tilde{U}),\end{split} (118)

where ν\nu is a measure over unitaries U~\tilde{U}. We examine each term individually. First, however, we prove the following lemma, which will be useful later.

Lemma 1.

Let |ψM\ket{\psi_{M}} an arbitrary but normalized vector in Morty 3’s description, a normalized vector |ψ|ψM+N+NN+\ket{\psi}\in{\mathcal{H}}_{\ket{\psi_{M}}}\cap{\mathcal{H}}_{+}^{N_{+}}\otimes{\mathcal{H}}_{-}^{N-N_{+}}, and PE{P}_{E} is a projector onto |E+N+NN+{\mathcal{H}}_{\ket{E}}\cap{\mathcal{H}}_{+}^{N_{+}}\otimes{\mathcal{H}}_{-}^{N-N_{+}}. Then

tr[PE|ψψ|]=|E|ψM|2.\tr[{P}_{E}\outerproduct{\psi}{\psi}]=\absolutevalue{\innerproduct{E}{\psi_{M}}}^{2}. (119)
Proof.

We expand |ψM\ket{\psi_{M}} in the basis of |E\ket{E},

|ψM=EαE|E,αE=E|ψM.\ket{\psi_{M}}=\sum_{E}\alpha_{E}\ket{E},\quad\alpha_{E}=\innerproduct{E}{\psi_{M}}. (120)

This means that whatever |ψ\ket{\psi} is, it must have a form of

|ψ=EαE|ψE.\ket{\psi}=\sum_{E}\alpha_{E}\ket{\psi_{E}}. (121)

where |ψEE\ket{\psi_{E}}\in{\mathcal{H}}_{E}. This can be proven by contradiction by applying the operator 𝒦{\mathcal{K}} and using orthogonality of |E\ket{E}. Since |ψE\ket{\psi_{E}} are eigenvectors of a Hermitian operator (Hamiltonian H{H}) corresponding to different eigenvalues, they must be orthogonal to each other. More specifically, we have PE|ψE=0{P}_{E}\ket{\psi_{E^{\prime}}}=0 for any EEE^{\prime}\neq E. This implies

PE|ψ=αE|ψE,{P}_{E}\ket{\psi}=\alpha_{E}\ket{\psi_{E}}, (122)

from which follows,

tr[PE|ψψ|]=ψ|PE2|ψ=|αE|2=|E|ψM|2.\tr[{P}_{E}\outerproduct{\psi}{\psi}]=\bra{\psi}{P}_{E}^{2}\ket{\psi}=\absolutevalue{\alpha_{E}}^{2}=\absolutevalue{\innerproduct{E}{\psi_{M}}}^{2}. (123)

We also make one useful observation.

Corollary 1.

The Lemma (1) also holds for any other normal operator acting on the actual Hilbert space, rather than Hamiltonian H=EEPE{H}=\sum_{E}E{P}_{E}, such as a unitary operator U=uuPuU=\sum_{u}u{P}_{u}.

Proof.

Any normal operator has orthogonal eigenspaces, which the only central assumption in Eq. (122). ∎

For the first term in Eq. (118), we have

U~tr[Hρ]ν(U~)=tr[Hρ]=tr[N+=0NHN,N+N+=0NλN+ρN+]=N+=0NλN+tr[HN,N+ρN+]=N+=0NλN+tr[EEPEN,N+rrρrN,N+]=N+=0NErErλN+tr[PEN,N+ρrN,N+]=N+=0NErErλN+|E|r|2=ErEr|E|r|2=tr[HMρM].\begin{split}\int_{\tilde{U}}\tr[{H}{\rho}]\nu(\tilde{U})&=\tr[{H}{\rho}]\\ &=\tr[\bigoplus_{N_{+}=0}^{N}H^{N,N_{+}}\bigoplus_{N_{+}=0}^{N}\lambda_{N_{+}}{\rho}_{N_{+}}]\\ &=\sum_{N_{+}=0}^{N}\lambda_{N_{+}}\tr[H^{N,N_{+}}{\rho}_{N_{+}}]\\ &=\sum_{N_{+}=0}^{N}\lambda_{N_{+}}\tr[\sum_{E}E{P}_{E}^{N,N_{+}}\sum_{r}r{\rho}_{r}^{N,N_{+}}]\\ &=\sum_{N_{+}=0}^{N}\sum_{E}\sum_{r}Er\lambda_{N_{+}}\tr[{P}_{E}^{N,N_{+}}{\rho}_{r}^{N,N_{+}}]\\ &=\sum_{N_{+}=0}^{N}\sum_{E}\sum_{r}Er\lambda_{N_{+}}\absolutevalue{\innerproduct{E}{r}}^{2}\\ &=\sum_{E}\sum_{r}Er\absolutevalue{\innerproduct{E}{r}}^{2}\\ &=\tr[{H}_{M}{\rho}_{M}].\end{split} (124)

We have used

tr[PEN,N+ρrN,N+]=|E|r|2.\tr[{P}_{E}^{N,N_{+}}{\rho}_{r}^{N,N_{+}}]=\absolutevalue{\innerproduct{E}{r}}^{2}. (125)

This follows from Eq. (108), which says that ρrN,N+{\rho}_{r}^{N,N_{+}} has support only on |r{\mathcal{H}}_{\ket{r}}, meaning that it has a spectral decomposition

ρrN,N+=jβjr|ψjr,N,N+ψjr,N,N+|,{\rho}_{r}^{N,N_{+}}=\sum_{j}\beta_{j}^{r}\outerproduct{\psi_{j}^{r,N,N_{+}}}{\psi_{j}^{r,N,N_{+}}}, (126)

|ψjr,N,N+|r\ket{\psi_{j}^{r,N,N_{+}}}\in{\mathcal{H}}_{\ket{r}}. Applying Lemma 1 yields

tr[PEN,N+ρrN,N+]=jβjrtr[PEN,N+|ψjr,N,N+ψjr,N,N+|]=jβjrtr[PEN,N+|E|r|2]=|E|r|2.\begin{split}\tr[{P}_{E}^{N,N_{+}}{\rho}_{r}^{N,N_{+}}]&\!=\!\sum_{j}\beta_{j}^{r}\tr[{P}_{E}^{N,N_{+}}\!\outerproduct{\psi_{j}^{r,N,N_{+}}}{\psi_{j}^{r,N,N_{+}}\!}]\\ &\!=\!\sum_{j}\beta_{j}^{r}\tr[{P}_{E}^{N,N_{+}}\absolutevalue{\innerproduct{E}{r}}^{2}]=\absolutevalue{\innerproduct{E}{r}}^{2}.\end{split} (127)

Before we can address the second term, we need to discuss the measure μ(U~)\mu(\tilde{U}). From the definition of the protocol, Morty 3 chooses a Haar measure over the unitary group μ(U~)\mu(\tilde{U}), where unitaries act on the Hilbert space 3{\mathcal{H}}_{3}, describing a chain of NN colorless particles. In reality, he applies some other random unitary on the Hilbert space, with some arbitrary measure over the degrees of freedom he cannot observe. Thus, to define the corresponding measure, we look for any measure ν\nu of the unitary group U~()=:X\tilde{U}({\mathcal{H}})=:X such that it reduces to the Haar measure μ\mu for U~M(3)=:Y\tilde{U}_{M}({\mathcal{H}}_{3})=:Y.

Let us define map f:XYf:X\rightarrow Y, using the prescription of Eq. (100),

U~=N+=0Nu~u~U~u~N,N+𝑓U~M=u~u~|u~u~|.\tilde{U}=\bigoplus_{N_{+}=0}^{N}\sum_{\tilde{u}}{\tilde{u}}\,\tilde{U}_{\tilde{u}}^{N,N_{+}}\overset{f}{\longrightarrow}\tilde{U}_{M}=\sum_{\tilde{u}}{\tilde{u}}\outerproduct{\tilde{u}}{\tilde{u}}. (128)

For a measurable subset AXA\subset X, we define the measure

ν(A)=Yκ(U~M,Af1(U~M))𝑑μ(U~M),\nu(A)=\int_{Y}\kappa(\tilde{U}_{M},A\cap f^{-1}(\tilde{U}_{M}))\,d\mu(\tilde{U}_{M}), (129)

where f1(U~M)f^{-1}(\tilde{U}_{M}) is a fiber and κ(U~M,)\kappa(\tilde{U}_{M},\bullet) is a conditional probability measure supported on that fiber. We require that

κ(U~M,f1(U~M))=1.\kappa(\tilde{U}_{M},f^{-1}(\tilde{U}_{M}))=1. (130)

This conditional probability encodes an arbitrary choice over the unobservable degrees of freedom, while the normalization condition ensures that integrating over them reproduces the Haar measure. This can be directly confirmed by integrating ν\nu over the preimage of a measurable subset BYB\subset Y,

ν(f1(B))\displaystyle\nu(f^{-1}(B)) =Yκ(U~M,f1(B)f1(U~M))𝑑μ(U~M)\displaystyle=\int_{Y}\kappa(\tilde{U}_{M},f^{-1}(B)\cap f^{-1}(\tilde{U}_{M}))\,d\mu(\tilde{U}_{M})
=Y𝟏B(U~M)κ(U~M,f1(U~M))𝑑μ(U~M)\displaystyle=\int_{Y}\mathbf{1}_{B}(\tilde{U}_{M})\cdot\kappa(\tilde{U}_{M},f^{-1}(\tilde{U}_{M}))\,d\mu(\tilde{U}_{M})
=Bκ(U~M,f1(U~M))=1𝑑μ(U~M)\displaystyle=\int_{B}\underbrace{\kappa(\tilde{U}_{M},f^{-1}(\tilde{U}_{M}))}_{=1}\,d\mu(\tilde{U}_{M})
=B𝑑μ(U~M)=μ(B).\displaystyle=\int_{B}d\mu(\tilde{U}_{M})=\mu(B).

Here, we used the indicator function 𝟏B:Y{0,1}\mathbf{1}_{B}:Y\to\{0,1\}, which selects whether a given fiber f1(U~M)f^{-1}(\tilde{U}_{M}) lies in the pre-image f1(B)f^{-1}(B),

𝟏B(y)={1,if yB,0,if yB.\mathbf{1}_{B}(y)=\begin{cases}1,&\text{if }y\in B,\\ 0,&\text{if }y\notin B.\end{cases} (131)

In summary, this construction ensures that the pushforward of ν\nu under ff recovers the original Haar measure,

fν=μ.f_{*}\nu=\mu. (132)

We will also need the following Lemma.

Lemma 2.

Let H=EEPE{H}=\sum_{E}E{P}_{E} be an actual Hamiltonian corresponding to Hamiltonian HM=EE|EE|{H}_{M}=\sum_{E}E\outerproduct{E}{E} in the Morty 3’s frame, i.e., PE{P}_{E} projects onto |E{\mathcal{H}}_{\ket{E}}. Let P~E:=UPEU\tilde{P}_{E}:=U^{\dagger}{P}_{E}U be a projector of a transformed Hamiltonian H~=UHU=EEP~E\tilde{H}=U^{\dagger}{H}U=\sum_{E}E\tilde{P}_{E}. There,

U=uuUu0|uU=\sum_{u}u\,U_{u}\oplus 0_{{\mathcal{H}}_{\ket{u}}^{\perp}} (133)

is an actual unitary operator corresponding to unitary operator UM=u|uu|U_{M}=\sum u\outerproduct{u}{u} in the Morty 3’s frame, given by construction in Eq. (102), where UuU_{u} is a unitary operator acting on subspace |u{\mathcal{H}}_{\ket{u}}. Then projectors P~E\tilde{P}_{E} project onto subspace UM|E{\mathcal{H}}_{U_{M}^{\dagger}\ket{E}}.

Proof.

We explicitly express the definition,

UM|E={|ψ|𝒦(|ψ)=αUM|E,α}.{\mathcal{H}}_{U_{M}^{\dagger}\ket{E}}=\big{\{}\ket{\psi}\,\big{|}\,{\mathcal{K}}(\ket{\psi})=\alpha U_{M}^{\dagger}\ket{E},\ \alpha\in\mathbb{C}\big{\}}. (134)

We have

UM|E=uu|uu|E,U_{M}^{\dagger}\ket{E}=\sum_{u}u^{*}\ket{u}\innerproduct{u}{E}, (135)

where both |u\ket{u} and |E\ket{E} are normalized and form orthogonal bases. Due to linearity of 𝒦{\mathcal{K}}, any vector |ψUM|E\ket{\psi}\in{\mathcal{H}}_{U_{M}^{\dagger}\ket{E}} must have a form

|ψ=αuuu|E|ψu,\ket{\psi}=\alpha\sum_{u}u^{*}\innerproduct{u}{E}\ket{\psi_{u}}, (136)

where |ψu|u\ket{\psi_{u}}\in{\mathcal{H}}_{\ket{u}} is normalized so that 𝒦(|ψu)=|u{\mathcal{K}}(\ket{\psi_{u}})=\ket{u}. To show that, let us define projector onto |u{\mathcal{H}}_{\ket{u}} as Pu{P}_{u}. Then, according to Lemma 1 combined with Corollary 1,

tr[Pu|ψuψu|]=|u|u|2=1.\tr[{P}_{u}\outerproduct{\psi_{u}}{\psi_{u}}]=\absolutevalue{\innerproduct{u}{u}}^{2}=1. (137)

Considering Pu|ψu=|ψu{P}_{u}\ket{\psi_{u}}=\ket{\psi_{u}}, we can rewrite this as,

|ψu|ψu|2=1.\absolutevalue{\innerproduct{\psi_{u}}{\psi_{u}}}^{2}=1. (138)

Thus, for 𝒦(|ψu)=|u{\mathcal{K}}(\ket{\psi_{u}})=\ket{u} to hold, |ψu\ket{\psi_{u}} must be also normalized.

Consider a spectral decomposition into orthonormal eigenvectors,

PE=k|ψEkψEk|{P}_{E}=\sum_{k}\outerproduct{\psi_{E}^{k}}{\psi_{E}^{k}} (139)

By definition, we have

P~E=kU|ψEkψEk|U\tilde{P}_{E}=\sum_{k}U^{\dagger}\outerproduct{\psi_{E}^{k}}{\psi_{E}^{k}}U (140)

where all eigenvectors U|ψEkU^{\dagger}\ket{\psi_{E}^{k}} are orthogonal and normalized. Consider a single eigenvector of this operator,

U|ψEk=(uuUu0|u)|ψEk=uuUuPu|ψEk=(1)uuu|EeiϕuEkUuPu|ψEkPu|ψEk=uuu|E|ψuEk,\begin{split}U^{\dagger}\ket{\psi_{E}^{k}}&=\Big{(}\sum_{u}u^{*}\,U_{u}\oplus 0_{{\mathcal{H}}_{\ket{u}}^{\perp}}\Big{)}\ket{\psi_{E}^{k}}\\ &=\sum_{u}u^{*}\,U_{u}{P}_{u}\ket{\psi_{E}^{k}}\\ &\overset{(1)}{=}\sum_{u}u^{*}\,\innerproduct{u}{E}\frac{e^{i\phi_{uEk}}U_{u}{P}_{u}\ket{\psi_{E}^{k}}}{\norm{{P}_{u}\ket{\psi_{E}^{k}}}}\\ &=\sum_{u}u^{*}\,\innerproduct{u}{E}\ket{\psi_{uE}^{k}},\\ \end{split} (141)

where |ψuEku\ket{\psi_{uE}^{k}}\in{\mathcal{H}}_{u} is a normalized vector. The step (1)(1) follows from the normalization: using Lemma 1 combined with Corollary 1, we have

UuPu|ψEk2=tr[Pu|ψEkψEk|]=|u|E|2,\norm{U_{u}{P}_{u}\ket{\psi_{E}^{k}}}^{2}=\tr[{P}_{u}\outerproduct{\psi_{E}^{k}}{\psi_{E}^{k}}]=\absolutevalue{\innerproduct{u}{E}}^{2}, (142)

and the complex phase eiϕuEke^{i\phi_{uEk}} ensures that u|EeiϕuEk=|u|E|\innerproduct{u}{E}e^{i\phi_{uEk}}=\absolutevalue{\innerproduct{u}{E}}. Since every eigenvector U|ψEkU^{\dagger}\ket{\psi_{E}^{k}} of P~E\tilde{P}_{E} has a form of Eq. (136), is normalized and orthogonal, and |ψuEk\ket{\psi_{uE}^{k}} are normalized, they form an orthonormal basis of UM|E{\mathcal{H}}_{U_{M}^{\dagger}\ket{E}}. In other words, P~E\tilde{P}_{E} projects onto this subspace. ∎

Finally, we address the second term,

U~tr[HUρU]ν(U~)=U~tr[HUextU~ρU~Uext]ν(U~)\displaystyle\int_{\tilde{U}}\tr[{H}U{\rho}U^{\dagger}]\nu(\tilde{U})=\int_{\tilde{U}}\tr[{H}U^{\mathrm{ext}}\tilde{U}{\rho}\tilde{U}^{\dagger}U^{\mathrm{ext}{\dagger}}]\nu(\tilde{U})
=U~tr[N+=0NHN,N+N+=0NUN,N+,extN+=0NU~N,N+N+=0NλN+ρN+N+=0NU~N,N+N+=0NUN,N+,ext]ν(U~)\displaystyle=\int_{\tilde{U}}\tr[\bigoplus_{N_{+}=0}^{N}H^{N,N_{+}}\bigoplus_{N_{+}=0}^{N}U^{N,N_{+},\mathrm{ext}}\bigoplus_{N_{+}=0}^{N}\tilde{U}^{N,N_{+}}\bigoplus_{N_{+}=0}^{N}\lambda_{N_{+}}{\rho}_{N_{+}}\bigoplus_{N_{+}=0}^{N}\tilde{U}^{N,N_{+}{\dagger}}\bigoplus_{N_{+}=0}^{N}U^{N,N_{+},\mathrm{ext}{\dagger}}\Bigg{]}\nu(\tilde{U})
=U~N+=0NλN+tr[HN,N+UN,N+,extU~N,N+ρN+U~N,N+UN,N+,ext]ν(U~)\displaystyle=\int_{\tilde{U}}\sum_{N_{+}=0}^{N}\lambda_{N_{+}}\tr[H^{N,N_{+}}U^{N,N_{+},\mathrm{ext}}\tilde{U}^{N,N_{+}}{\rho}_{N_{+}}\tilde{U}^{N,N_{+}{\dagger}}U^{N,N_{+},\mathrm{ext}{\dagger}}\Big{]}\nu(\tilde{U})
=U~N+=0NλN+tr[EEPEN,N+UN,N+,extU~N,N+rrρrN,N+U~N,N+UN,N+,ext]ν(U~)\displaystyle=\int_{\tilde{U}}\sum_{N_{+}=0}^{N}\lambda_{N_{+}}\tr[\sum_{E}E{P}_{E}^{N,N_{+}}U^{N,N_{+},\mathrm{ext}}\tilde{U}^{N,N_{+}}\sum_{r}r{\rho}_{r}^{N,N_{+}}\tilde{U}^{N,N_{+}{\dagger}}U^{N,N_{+},\mathrm{ext}{\dagger}}\Big{]}\nu(\tilde{U})
=U~N+=0NλN+ErErtr[U~N,N+UN,N+,extPEN,N+UN,N+,extU~N,N+ρrN,N+]ν(U~)\displaystyle=\int_{\tilde{U}}\sum_{N_{+}=0}^{N}\lambda_{N_{+}}\sum_{E}\sum_{r}Er\tr[\tilde{U}^{N,N_{+}{\dagger}}U^{N,N_{+},\mathrm{ext}{\dagger}}{P}_{E}^{N,N_{+}}U^{N,N_{+},\mathrm{ext}}\tilde{U}^{N,N_{+}}{\rho}_{r}^{N,N_{+}}\Big{]}\nu(\tilde{U})
=U~N+=0NλN+ErErtr[P~EN,N+ρrN,N+]ν(U~)\displaystyle=\int_{\tilde{U}}\sum_{N_{+}=0}^{N}\lambda_{N_{+}}\sum_{E}\sum_{r}Er\tr[\tilde{P}_{E}^{N,N_{+}}{\rho}_{r}^{N,N_{+}}\Big{]}\nu(\tilde{U}) (143)
=(1)U~N+=0NλN+ErEr|r|U~MUMext|E|2ν(U~)\displaystyle\overset{(1)}{=}\int_{\tilde{U}}\sum_{N_{+}=0}^{N}\lambda_{N_{+}}\sum_{E}\sum_{r}Er\absolutevalue{\bra{r}\tilde{U}_{M}^{\dagger}U_{M}^{\mathrm{ext}{\dagger}}\ket{E}}^{2}\,\nu(\tilde{U})
=(2)U~ErEr|r|U~MUMext|E|2ν(U~)\displaystyle\overset{(2)}{=}\int_{\tilde{U}}\sum_{E}\sum_{r}Er\absolutevalue{\bra{r}\tilde{U}_{M}^{\dagger}U_{M}^{\mathrm{ext}{\dagger}}\ket{E}}^{2}\,\nu(\tilde{U})
=(3)U~M(f1(U~M)ErEr|r|U~MUMext|E|2κ(U~M,dU~))μ(U~M)\displaystyle\overset{(3)}{=}\int_{\tilde{U}_{M}}\left(\int_{f^{-1}(\tilde{U}_{M})}\sum_{E}\sum_{r}Er\absolutevalue{\bra{r}\tilde{U}_{M}^{\dagger}U_{M}^{\mathrm{ext}{\dagger}}\ket{E}}^{2}\,\kappa(\tilde{U}_{M},d\tilde{U})\right)\mu(\tilde{U}_{M})
=(4)U~MErEr|r|U~MUMext|E|2μ(U~M)\displaystyle\overset{(4)}{=}\int_{\tilde{U}_{M}}\sum_{E}\sum_{r}Er\absolutevalue{\bra{r}\tilde{U}_{M}^{\dagger}U_{M}^{\mathrm{ext}{\dagger}}\ket{E}}^{2}\,\mu(\tilde{U}_{M})
=U~Mtr[HMUMextU~MρMU~MUMext]μ(U~M).\displaystyle=\int_{\tilde{U}_{M}}\tr[{H}_{M}U_{M}^{\mathrm{ext}}\tilde{U}_{M}{\rho}_{M}\tilde{U}_{M}^{\dagger}U_{M}^{\mathrm{ext}{\dagger}}]\mu(\tilde{U}_{M}).

There we have defined coarse-graining

P~EN,N+=U~N,N+UN,N+,extPEN,N+UN,N+,extU~N,N+,\tilde{P}_{E}^{N,N_{+}}=\tilde{U}^{N,N_{+}{\dagger}}U^{N,N_{+},\mathrm{ext}{\dagger}}{P}_{E}^{N,N_{+}}U^{N,N_{+},\mathrm{ext}}\tilde{U}^{N,N_{+}}, (144)

corresponding to Hamiltonian H~N,N+=EP~EN,N+\tilde{H}^{N,N_{+}}=\sum E\tilde{P}_{E}^{N,N_{+}}, where, according to Lemma 2, each P~EN,N+\tilde{P}_{E}^{N,N_{+}} projects onto subspace

UMUMext|E+N+NN+.{\mathcal{H}}_{U_{M}^{\dagger}U_{M}^{\mathrm{ext}{\dagger}}\ket{E}}\cap{\mathcal{H}}_{+}^{N_{+}}\otimes{\mathcal{H}}_{-}^{N-N_{+}}. (145)

This is because both UU and UextU^{\mathrm{ext}} follow the same structure as required by the Lemma, so the statement follows by applying the Lemma twice. In step (1)(1), we have used this property together with Eq. (125). In step (2)(2), we have summed the eigenvalues, N+=0NλN+=1\sum_{N_{+}=0}^{N}\lambda_{N_{+}}=1, because the summand no longer depends on N+N_{+}. In step (3)(3), we have explicitly substituted the measure ν\nu using Eq. (129). In step (4)(4), we trivially integrated over the inner degrees of freedom, since the integrand does not depend on the fiber variables, and the conditional measure κ\kappa is normalized on each fiber.

Inspecting Eqs. (116) and (118), we have proven that

W𝒞=W𝒞M.W_{\mathcal{C}}=W_{{\mathcal{C}}_{M}}. (146)

Thus, given the assumptions (a1(a1a8)a8), Morty 3 extracts exactly the amount of energy he predicts, without the knowledge of the existence of two types of particles.

The average work extracted per copy when extracting from multiple copies simultaneously, which is defined by [72],

W𝒞=limL1LU~tr[HLρL]tr[HLUextU~LρLU~LUext]μ(U~),\begin{split}W_{{\mathcal{C}}}=\lim_{L\rightarrow\infty}\frac{1}{L}\int_{\tilde{U}}&\tr[{H}^{L}\rho^{\otimes L}\big{]}\\ -&\tr[{H}^{L}U^{\mathrm{ext}}\tilde{U}^{\otimes L}{\rho}^{\otimes L}\tilde{U}^{\otimes L{\dagger}}U^{\mathrm{ext}{\dagger}}\big{]}\mu(\tilde{U}),\end{split} (147)

where HL=HI+IH{H}^{L}={H}\otimes{I}\otimes\cdots+{I}\otimes{H}\otimes\cdots, proceeds analogously.

References