{'fraction_non_alphanumeric': 0.04133275507702751, 'fraction_numerical': 0.0073727678992325536, 'mean_word_length': 4.663551900897053, 'pattern_counts': {'":': 1, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 8, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Holism as the empirical significance of symmetries∗

Not all symmetries are on a par. For instance, within Newtonian mechanics, we seem to have a good grasp on the empirical significance of boosts, by applying it to subsystems. This is exemplified by the thought experiment known as Galileo’s ship: the inertial state of motion of a ship is immaterial to how events unfold in the cabin, but is registered in the values of relational quantities such as the distance and velocity of the ship relative to the shore. But the significance of gauge symmetries seems less clear. For example, can gauge transformations in Yang-Mills theory—taken as mere descriptive redundancy—exhibit a similar relational empirical significance as the boosts of Galileo’s ship? This question has been debated in the last fifteen years in philosophy of physics. I will argue that the answer is ‘yes’, but only for a finite subset of gauge transformations, and under special conditions. Under those conditions, we can mathematically identify empirical significance with a failure of supervenience: the state of the Universe is not uniquely determined by the intrinsic state of its isolated subsystems. Empirical significance is therefore encoded in those relations between subsystems that stand apart from their intrinsic states.


Overview of the debate and my position within it
In its broadest terms, a symmetry is a transformation of a system which preserves the values of a relevant (usually large) set of physical quantities. Of course, this broad idea is made precise in various different ways: for example as a map in the space of states, or on the set of quantities; and as a map that must respect the system's dynamics, e.g. by mapping solutions to solutions or even by preserving the value of the Lagrangian functional on the states.
The broad idea is also associated with various debates. 1 For example, should we say that a symmetry transformation applied to the whole universe cannot yield a different physical state of affairs? And relatedly: should we prefer a reduced i.e. quotiented formalism, so that if presented with a state space S partitioned into the orbits of a group of symmetries G acting on S, we prefer the reduced state space whose elements are the orbits, i.e. [s] ∈ S/ ∼ ? (where s ∼ s (s, s ∈ S) means that s and s are related by a symmetry transformation: s = ξ · s, for ξ ∈ G and · some action of G on S, and square brackets denote an entire equivalence class).
These "defining features" of symmetries are of central concern for one recent philosophical debate. More specifically, the debate is about whether gauge symmetries can have a direct empirical significance. Of course, all hands agree that symmetries have various important empirical implications. The obvious examples come from the Noether theorems: the restrictions on the equations of motion entailed by Noether's second theorem, and the (approximately) conserved charges given by Noether's first theorem. In other words, symmetries imply the (extraordinary) facts that charges are conserved.
Accordingly, in this debate, such familiar implications are often called the 'indirect' empirical significance of symmetries (IES); and through them, symmetries carry immense explanatory power.
But some familiar symmetries of the whole Universe, such as velocity boosts in classical or relativistic mechanics (Galilean or Lorentz transformations), have a direct empirical significance when applied solely to subsystems. Thus Galileo's famous thought-experiment about the ship-that a process involving some set of relevant physical quantities in the cabin below decks proceeds in exactly the same way whether or not the ship is moving uniformly relative to the shore-shows that sub-system boosts have a direct, albeit relational, empirical significance. For though the inertial state of motion of the ship is undetectable to experimenters confined to the cabin, yet the entire system, composed of ship and sea 2 registers the difference between two such motions, namely in the different relative velocities of the ship to the water. 3 Such examples rely on what are called 'external symmetries', i.e. symmetries which shift spacetime points around.
So the question arises: Can other symmetries-especially gauge symmetries-have a similar direct empirical significance when applied to subsystems?
For gauge symmetries are normally taken to encode descriptive redundancy: a view I will endorse. That is, they arise in a formalism that uses more variables than there are physical degrees of freedom in the dynamical system described. (They are also internal: unlike a boost or spatial translation, they do not shift spacetime points around).
This descriptive redundancy means that the natural answer to our question is 'No'. For surely, while a "freedom to redescribe" may have some indirect empirical implications, 4 it could not have the content needed for a direct empirical significance, like the one illustrated by Galileo's ship. This 'No' answer was developed in detail by Brading and Brown (Brading & Brown, 2004) (henceforth BB) in response to various discussions such as Kosso (Kosso, 2000). They take themselves-I think rightly, in this respect-to be articulating the traditional or orthodox answer.
The 'Yes' answer has been argued for by Greaves and Wallace (Greaves & Wallace, 2014) (henceforth GW), building on (Healey, 2009). I will agree with some aspects of both BB's and GW's analysis of symmetries. But, unlike either GW or BB, I will recast the topic to focus on gauge-invariant information about-i.e. states ofregions. My own conclusion will be that only a finite subset of gauge-transformations, usually called 'global' (but here called 'rigid', cf. section 2.2), can have DES.
Before we can summarize the shape of the debate, we need a couple of definitions, which I will now informally sketch (exact definitions will appear in later sections). First, the broad notion of 'direct empirical significance' is a matter of the existence of transformations of the universe possessing the following two properties (articulated in this way by BB): (i) the transformation should lead to an empirically different scenario, and (ii) the transformation should be a symmetry of the subsystem in question (e.g. Galileo's ship).
If such transformations exist, the symmetries of the theory in question are said to exhibit direct empirical significance (henceforth 'DES'). The empirical significance is to be witnessed by observers that lie outside the subsystem-it cannot be detected by those confined within it. Therefore, DES combines an inside and an outside perspective and, in this limited sense, acquires an epistemic dimension, or at least one that considers physical information as it is intrinsically accessible within a subsystem.
BB argue that gauge symmetries cannot exhibit DES according to (i) and (ii), while GW argue that they can. I myself will argue for a 'Yes' answer, but will approach the question in terms of gauge-invariant information.
By thus proceeding in terms of gauge-invariant information, I will identify DES as defined by (i) and (ii) above with a failure of Global Supervenience on Subsystems (GSS): 5 this failure is a form of holism-hence my title.
I will show that relational DES occurs when the gauge-invariant global state fails to supervene on the collection of intrinsic, gauge-invariant local states of the components of some arbitrary partition (of space or spacetime). Here I should make it clear that I am not claiming ontological priority for the subsystems composing the whole. Subsystems don't exist "before" the whole. The division of the Universe into subsystems is not mandatory, but it appears in item (ii): so we must consider what physical information is intrinsic to a subsystem when evaluating the direct empirical significance of symmetries. This is why, when environment and subsystems are on a par, DES can be rephrased a matter of global supervenience.
GSS is upheld when the intrinsic physical states of those subsystems composing the whole uniquely determine the physical state of the whole, without the requirement for additional relational information about those subsystems. When GSS fails, there can either be many physical states of the whole which are formed from the same physical states of the individual subsystems-in which case one is missing some relational information-or there can be no valid states of the whole, in which case the subsystems states are incompatible.

External sophistication
If the standard notion of DES in (i) and (ii) is to be identified with a failure of GSS, we first need to develop a physically meaningful notion of composition of those subsystems that possess descriptive redundancy.
Thus, we are led to revisit one other important debate in the philosophy of gauge, already mentioned at the start of the previous section. Namely, given a theory whose set of universe-descriptions-'states'-is partitioned by a group of symmetries, we can take one of two attitudes: (a) Reduction:try to write down a reduced theory whose states correspond to the cells of the partition; or (b) Sophistication about symmetries 6 :-resist quotienting the given theory, but take two symmetry-related states to be isomorphic.
I'll advocate a third position (c), which applies only in the presence of subsystems. Thus, for the entire universe, I unequivocally endorse reduction. In the same vein, I will also assume that the theory in question empirically discerns two different states s 1 and s 2 of the universe if and only if [s 1 ] = [s 2 ].
For subsystems, the question is more subtle, for there are two perspectives we can take: one from the inside, or intrinsic; and one from the outside, or extrinsic, in accord with the definition of DES. If we are interested in discriminating between (intrinsically) distinct physical possibilities, then surely those states which cannot be (intrinsically) discriminated are to be counted as one, and so 'reduction' still applies.
But, when we combine the subsystem with the rest of the world, we are required to exploit subsystem symmetries in a real physical sense: as emphasized first by Rovelli (Rovelli, 2014), reduced representations of subsystems cannot be straightforwardly coupled to each other. For coupling, we need to keep gauge-variant elements in the theory. In (Gomes, 2019a), it was similarly argued that reduction should only be endorsed for the entire universe: coupling regional states may require a re-expression of the states as particular gauge representatives of the physical states; 7 and thus, for regions, gauge-information should not be entirely eliminated.
Therefore, I have argued that we should delineate a third attitude (Gomes, 2019a): (c) External sophistication and internal reduction of subsystem symmetries:-Fix unique representations of the intrinsic physical states of the subsystems (i.e. from the internal perspective), but then allow these representations flexibility from an external perspective, as e.g. required for the smooth coupling of the states of subsystems.
Take the example used in (Rovelli, 2014): a non-relativistic classical system of N particles with translational invariance. From the intrinsic perspective of the subsystem, one could eliminate redundancy by taking the inter-particle distances as a new, autonomous set of coordinates, but this would leave no 'handle' for other subsystems 6 Sophistication has long been advocated for diffeomorphisms and spacetime metrics (see (Pooley, 2013) and references therein). The nomenclature was originally used for sophisticated substantivalism: points of spacetime may have identity, but this identity comes only through the complex web of inter-relations between different fields of the theory, and is in this way entirely dependent on the state. Some general features of this position have more recently also been suggested for gauge theories (Caulton, 2015;Dewar, 2017;Gomes, 2019a). Dewar describes it thus: "Whereas a reduced theory converts a class of symmetry-related models into a single model, sophistication converts a class of symmetry-related models into a class of isomorphic models." Some philosophers (e.g. (Healey, 2007, Sec. 4.2) and (Maudlin, 1998)) have resisted the analogy between the descriptive redundancy of the metric in gravitational theories and of the connection in gauge theories. Pace these philosophers, I see no reason for their resistance, but that is a topic for another day. 7 Rovelli focused on the coupling between different types of fields, or, in the finite-dimensional case, on the coupling of two different particle systems (Rovelli, 2014). In (Gomes, 2019a) I extended that requirement to the coupling of fields in regions. See also (Dougherty, 2017) for a 'stack-theoretic' argument emphasizing the problems of reduction for the coupling of regions: his notion of separability requires the preservation of gaugerelated representations, to be kept as isomorphic but not identified. That is, I construe Dougherty's view as a defense of position (b) motivated by the composition of subsystems (but using stack-theory). Also (implicitly) using a stack-theoretic approach, (Nguyen, Teh, & Wells, 2018) emphasize that gauge transformations are not just re-descriptions, but also parameterize the different ways in which regional field spaces can be composed into global field spaces.
to couple to. The joint system of autonomous coordinates for two such sets of particles (say N 1 'red' and N 2 'green') cannot express different ways of composing the subsystems-whether the center of mass of the 'reds' are five or ten feet away from the center of mass of the 'greens' along some direction does not register in these variables. From a degree of freedom count, we have clearly gone overboard: we have eliminated six degrees of freedom of the joint system-the position of the center of mass of 'reds' (three degrees of freedom) and the position of center of mass of 'greens' (three degrees of freedom)-when only three were eliminable: the position of center of mass of {reds} ∪ { greens}).
On the other hand, fixing the isolated subsystem's coordinates by reference to its center of mass, while leaving the center of mass embedded in Euclidean space, still affords us enough flexibility to characterize both the subsystems intrinsic degrees of freedom and a rigid subsystem translation with respect to another subsystem. This is a very simple example of a "covariant gauge-fixing", and in practice, it is how we implement option (c).
In my analysis of DES in the context of holism, the flexibility allowed by option (c), 'External sophistication' for short, is employed for melding the subsystems' physical content into the physical content of the joint state. 8 Option (c) allows us to have our cake and eat it too: we can both parametrize the intrinsic physical possibilities of the subsystems in a one-to-one manner, and yet keep track of those degrees of freedom that would be redundant from the intrinsic point of view but which must be retained, to be pressed into service for composing subsystems (Gomes, 2019a).

Non-technical summary
Note that the difference between the two scenarios above-between composing the reduced subsystems and composing the covariantly gauge-fixed subsystems-involves the space of possibilities of the subsystems. The mismatch is only apparent when we consider all the possible intrinsic physical states of the subsystems, say, n in number, and the physical states of the universe that are compatible with each n-tuple of physical subsystem states. That is, it is only apparent in the same context of GSS, or lack thereof.
Indeed, a given state of the entire universe will be taken to uniquely fix the physical subsystem states that constitute it, as well as any extra relational information between them. Given first the state of the universe, there is no question about composition: only of decomposition, which is straightforward. But GSS is not about decomposition, and neither is DES, even though this confusion permeates the debate, as we will briefly see in Section 4. To avoid this confusion, it is best to start with the information that is intrinsically available to the subsystems and assess which universal physical states can emerge from composition.
To be clear, in this paper, I will only countenance the case where a given physically allowed state of the universe decomposes into physically allowed states of its subsystems. I know of no examples violating this assumption. Therefore, we will focus on the more interesting of the two cases of failure of GS: namely, the one where, given just the intrinsic physical states of the subsystems, there are physically distinct possibilities to join these states into a state of the Universe. That is, the relation between states of the Universe and states of its subsystems are many-to-one, because there is relevant relational information that cannot be registered intrinsically within each subsystem. In these cases, we will say there is Global Non-Supervenience on Subsystems (GNSS).
Schematically: if the subsystems are 'sea' and 'ship', and there are equivalence relations, ∼, applicable to states of subsystems and of the whole, and given the physical (i.e. "gauge-invariant") states [s sea ], [s ship ] and [s sea and ship ], there is a many-to-one relation, encoded by the set I: and [s sea and ship ] (i) = [s sea and ship ] (i ) if and only if i = i . In this case, the set I that parameterizes the many-to-one relation is the (inhomogeneous) Galilean group (which is a semi-direct product ( ) of boosts and the group of translations and rotations).
One can thus see that the variety of states of the whole is not encoded in either subsystem: 9 it is encoded in the relations between the two subsystems, as denoted by ∪ (i) .
In certain situations, such as in Galileo's ship, there is remarkable order to this variety of physical states of the whole, an order also encoded in the structure of I. Namely, each element of this variety can be transformed into another by a subsystem symmetry which does not extend beyond the boundary of the subsystem. In other words, I carries the structure of finite-dimensional symmetry group of the subsystem. For Galileo's ship, these are the Galilean transformations; and for the gauge theory, we will see that they are (sub)groups of the Lie group characterizing the theory. In this manner, DES becomes a matter of GNSS.
In sum: technicalities apart, my main claim is that both Galilean boost symmetry for particle systems and gauge symmetry for certain field theories carry direct empirical significance through GNSS. This holism is empirically significant, since it registers physical-i.e. gauge-invariant-differences in the entire system and we take such differences to lead to empirically distinguishable universes. Moreover, the implied under-determination of the physical state of the whole universe by the physical state of its subsystems is encoded in a subsystem symmetry, but only as seen from the 'outside perspective'; again in accord with the above construal of DES.
2 Direct empirical significance I start, in section 2.1, by construing DES in terms of properties of transformations of the Universe. This description of DES immediately runs into some cumbersome notation when applied to gauge theories. Therefore, in section 2.2 I introduce a new terminology which better distinguishes the relevant categories of transformations. Then, having got the right nomenclature for addressing DES in the context of gauge theories, in section 2.3 I apply it to re-express the debate in these better terms. In section 2.4, I then proceed to offer an appetizer of my criticism of BB and GW's construals of DES, and supplant those construals with my own. Thus here I describe the relation between GNSS and DES. 9 Agreed, several approximations must be in place for this statement: for surely, with the right equipment (such as a window), the person within the cabin could discern movement of the ship from within, and different movements of the ship could create different sorts of eddies and turbulence in the sea. This sort of idealization is ubiquitous in physics, and generally unproblematic. This is echoed in (Greaves & Wallace, 2014, p. 52, footnote 7): " 'Dynamical isolation' is of course approximate in practice, as no proper subsystem can in practice be perfectly isolated. It is also relative to the observation capabilities of relevant observers (if Galileo had included a GPS tracker in his list of the cabin's accoutrements, things would have been rather different)." See also (Earman, 2019) for the utility of approximate symmetries in physics.
2.1 DES as a transformation of the universe BB frame the definition of DES in terms of two conditions. First, a transformation cannot be a symmetry of the entire universe, otherwise it would not have any direct empirical significance. But second, it needs to act as a symmetry for subsystems, otherwise the transformation in question could hardly be called a symmetry. Thus Brading and Brown define: Definition 1 (Direct Empirical Significance (DES) as a transformation) A symmetry has direct empirical significance if it is specified by a transformation that satisfies the following conditions: 10 1. Transformation Condition: the transformation must yield an empirically different scenario. In our words: the transformation in question is not a symmetry of the world as a whole.
2. Subsystem Symmetry Condition: The evolution of the untransformed and transformed subsystems must be empirically indistinguishable from the interior point of view. In our words: the transformation should count as a symmetry when restricted to the subsystems composing the entire system.This is the subsystem symmetry with DES. 11 Thus in the example of Galileo's ship, the entire system-both ship and sea-is in different states if the ship is heading through calm waters towards the North-West at 10km/h or towards the South at 20km/h. The entire system thus satisfies the first condition (Transformation). Nonetheless, inside the cabin, you would not be able to distinguish the two scenarios (cf. footnote 9): so the subsystem satisfies the second condition.
As is clear from the above, the empirical significance cannot be witnessed by observers within the subsystem and, in that way, DES combines an inside and an outside perspective: inside it considers intrinsic physical information that is accessible within the subsystem; outside it considers an overall change in the state of the universe.

Two distinctions
At this point in the discussion, standard terminology gets in the way of clarity. When used in conjunction with subsystem-Universe distinctions, the words 'local' and 'global' acquire other possible meanings, and may pull intuition in different directions (see p. 10 Both Teh (Teh, 2016) and GW add a condition of dynamical isolation between the two subsystems. Teh takes this to justify an asymptotic treatment for the subsystem in question. We won't need to make this isolation condition explicit: it emerges from the criteria.
11 Although BB distinguish a subsystem and its environment, and thus have only the singular 'subsystem' in their definition, it is customary to focus on relational DES, i.e. DES with the environment taken on a par with the subsystem in question. Thereby, the environment is taken as just one more subsystem, and transformations of the environment are to be considered just as much as transformations of the subsystem. In particular, this demotion of the environment to subsystem status means that one excludes the exclusively external relations of the environment from its state. A non-relational definition of DES would not require condition 2 to apply to the environment (cf. footnote 15), that is, it would include transformations that are not symmetries of the intrinsic state of the environment. The standard argument against non-relational DES is that only relational DES has a principled connection between a subsystem symmetry and physically distinct universes. To glean the difficulty with non-relational DES, consider the following example: in the Galilieo ship scenario, imagine a transformation that arbitrarily changes the configuration of the interior of the beach-taken as part of the environment-but otherwise keeps the subsystem-intrinsic and all other relational information untouched. If definition 1 did not require the physical state of the environment to be preserved, such a transformation would be included and therefore correspond to a symmetry with DES. But an arbitrarily changing beach has little to do with symmetries. See also footnote 25. 648 of (Brading & Brown, 2004)). Therefore, it is useful to introduce a nomenclature that distinguishes these meanings.
• Universal : A universal transformation is one that applies to the world as a whole.
The set of universal transformations may depend on an infinite or finite number of parameters.
• Regional : A regional transformation is one that applies only to a subsystem of the world. 12 The set of regional symmetry transformations may also depend on an infinite or finite number of paramaters.
• Malleable (aka 'local'): A malleable transformation depends on an infinite number of parameters: e.g. it is specified by an arbitrary smooth function over a given manifold or region of a manifold. Here, the usual label is 'local'. But using 'local' invites confusion with the above category, 'regional'. I will therefore prefer the term 'malleable'. A malleable symmetry transformation can be either regional or universal.
• Rigid (aka 'global'): A 'rigid' symmetry transformation depends only on a finite number of parameters. This is to be contrasted with malleable. The usual label is 'global'. But again, this term invites confusion, namely with the above category, 'universal'. So I will prefer the term 'rigid'. A rigid symmetry transformation can be either regional or universal.
Therefore a symmetry transformation may lie in any of the following four combinations of the above categories: regional and rigid, regional and malleable, universal and malleable, or universal and rigid. 13 Regional transformations are under-studied in the physics literature, but are known to hide many surprises: see e.g. (Regge & Teitelboim, 1974;Balachandran, Chandar, & Momen, 1996;Donnelly & Freidel, 2016;Gomes, Hopfmüller, & Riello, 2019;Gomes, 2019a). As to universal transformations, the rigid ones are familiar; they are associated with the standard treatments of Noether's first theorem, and thus correspond to conserved charges (Noether, 1917;Olver, 1986;Brading & Brown, 2000;Butterfield, 2007;Kosmann-Schwarzbach & Schwarzbach, 2011). The malleable universal transformations are associated to constraints, or relations between the equations of motion (such as the Hamiltonian constraints or the Bianchi identities of general relativity and the Gauss constraint in electromagnetism). 14 But since these constraints guarantee charge conservation as well, we note that it was argued in (Gomes, 2019b)-in the nomenclature introduced above-that the real power of the Noether theorems arises from weakening a rigid symmetry to a malleable one, which thereby enforces compatibility between charge conservation and the dynamics of the fields that the charges couple to.

The debate re-expressed
Using this nomenclature, we can re-express Definition 1 and better address the subtleties of applying it to gauge symmetry. Thus Definition 1 says that DES arises if there are transformations that are not universal symmetries and yet whose restrictions are regional symmetries. The question is which, if any, regional malleable symmetry can be obtained in this way, and thus be awarded DES.
Finite-dimensional theories, i.e. ones which do not involve fields, generally only have rigid symmetries (such as translations, etc.). In those cases, the strictly regional symmetries can give different values for appropriate physical quantities, viz. relational quantities relating the transformed subsystem to the rest of the universe. This is of course what Galileo's ship illustrates. In this case, the clear distinction between universal and regional rigid symmetries is illustrated in an uncontroversial case of DES.
But the situation for malleable symmetries seems different.
In certain examples, the generators of malleable symmetries are spacetime vector fields; in others, they are Lie-algebra-valued scalar fields, acting on an internal space over each spacetime point. In any case, it is easy to imagine a malleable symmetry acting on a region of spacetime and not on another, which can thereby serve as the environment, or reference system. In this case, the malleable transformation should smoothly tend to the identity at the boundary between the regions, lest it create discontinuities in the fields. But then, it seems we could suitably extend any such regional malleable symmetry to the rest of the Universe simply with the identity transformation. The conjunction of the two regional transformations-one that tends to the identity at the boundary and the other the identity on the rest of the universe-would be a universal malleable symmetry, and thus could not have DES.
According to BB (Brading & Brown, 2004), this is precisely the case: the environment can be assumed to be untransformed, and therefore, to avoid discontinuities and due to their malleability, gauge symmetries would violate Definition 1's first condition (Transformation), i.e. the requirement of an empirically distinguishable scenario, and so would have no empirical significance in the way that regional rigid symmetries do. As BB write: "Thus, a transformation applied to one subsystem will involve the other subsystem, even if only because the transformation of the gauge field goes smoothly to the identity. In conclusion, there can be no analogue of the Galilean ship experiment for local gauge transformations, and therefore local gauge symmetry has only indirect empirical significance (being a property of the equations of motion)." (p. 657) GW articulate DES for gauge theory in a manner that fosters DES for gauge symmetries. They focus on subsystems as given by regions, and thereby identify transformations possessing properties 1 and 2 of Definition 1 by first formulating the putative effects of such transformations on the gauge fields in these regions.
In particular, they focus their attention on relational DES. That means they consider the environment to be on a par with the subsystem in question (cf. footnote 11). Thus, in particular, the transformation in Definition 1 must obey property 2 -it must also be a symmetry of the environment of the subsystem. 15 In this case we can diagnose DES as originating in the relations between the subsystems (and, ultimately, as I will show, in a failure of supervenience of the global state on the intrinsic states of the subsystems).
More precisely, for a given subsystem state s, they claim the relational DES transformations are in 1-1 correspondence with the following quotient between two groups of transformations: where G(s |∂ ) are the gauge transformations of the region which preserve the state s at the boundary of the region, and G Id are the gauge transformations of the region which are the identity at the boundary. Here the equivalence class between two transformations ξ, ξ ∈ G(s |∂ ) is taken as ξ ∼ ξ iff ξ = ηξ for some η ∈ G Id . 16 The rough idea is that even if certain transformations would not preserve all possible states at the boundary, they will preserve some of those states. 17 So, G Id is a subgroup of G(s |∂ ), and one would like to 'factorize out' from those transformations those that would preserve all states at the boundary (i.e. the boundary-identity transformations): so one defines the quotient group G(s |∂ )/G Id . Several assumptions go into the results of GW and BB, and we will challenge some of these once we have described our own results (see Section 4).

My own position in the debate: rigid variety
Overall, I will argue for a position not considered by either the GW or the BB camp: an appropriate selection of rigid regional symmetries-but not all the malleable ones!can retain direct empirical significance (DES) in both the finite-dimensional case and in the field-theoretic case. In very specific circumstances, and according to a precise method, the rigid symmetries will be identified among the malleable ones-they are the ones that leave the gauge potential invariant but which shift the matter fields, and they will be the only ones lying in the kernel of a (configuration-dependent, in the non-Abelian case) elliptic differential operator-and they will be transformation with DES as per Definition 1.
Rigid variety and relational DES In the following, to make matters concrete, the fieldtheories I will focus on are general classical Yang-Mills theories in the presence of matter. The spacetime fields in question will be the standard, smooth gauge potentials, A, and charged scalars, ψ, valued in the appropriate vector spaces, which I discuss in more detail below, in Section 3.1. I denote the doublet of these two fields by: ϕ = (A, ψ), and the space of such doublets by Φ ϕ. The subsystems will consist of regions in the manifold, to which the fields get restricted. Thus the subsystems in questions are regional, and thus we refer to Global Supervenience on Regions (to avoid unnecessarily detailed acronyms, we will denote these also by GSS). This initial set-up lace, 2014, p.86,87). See footnotes 11 and 25 here for more on the 'principled connection' and the treatment of the environment as solely a reference and not a subsystem; and (Ramirez & Teh, 2019) for one possible interpretation of the term 'principled connection' for non-relational DES. 16 A group is just a set closed under an associative invertible binary operation. That is, if ξ, ξ ∈ G, then ξξ ∈ G, (ξξ )ξ = ξ(ξ ξ ) and for all ξ ∈ G there exists a ξ −1 such that ξ −1 ξ =Id, where ξId= ξ (here we will not need to distinguish left and right inverses). The quotient is well-defined for normal subgroups: namely, G/H is well-defined as a group if given η ∈ H, ξηξ −1 ∈ H for all ξ ∈ G. This holds in the example above.
17 One could think of it as follows: some state s might be 'periodic' in what it says about the boundary; so that a transformation with the same period at the boundary will fix the state s at the boundary even though the transformation is not the identity there (i.e does not fix all states s at the boundary. That is: when the transformation is restricted to boundary, it does not fix each state thus restricted to boundary).
is standard in the debate about the DES of gauge symmetries, and is applicable to all of the approaches considered here.
DES combines an inside and an outside perspective: inside it considers physical information that is intrinsically accessible to subsystems; outside it considers an overall change in the state of the universe. Definition 1 can be formulated as a failure of GSS because its requirement 2 is about information that is intrinsically available to a subsystem. Therefore, for relational DES, a transformation that does not change the intrinsic physical states of its subsystems and yet changes the physical state of the Universe must be changing the relations between the subsystems. Therefore, for DES to exist, there must be a physical variety of universes which are made up from the same (physically) subsystem states.
Accordingly, within this context of Yang-Mills fields, I define: 18 Definition 2 (GSS) Given a manifold Σ, that is decomposed as Σ = Σ + ∪ Σ − , along the boundary ∂Σ ± = Σ + ∩ Σ − =: S; given a universal field supported on ϕ and the regional fields, ϕ ± supported on Σ ± , GSS holds just in case the joint gauge-invariant contents of ϕ ± is compatible with a unique gauge-invariant content of ϕ. That is, the doublet of regional physical (i.e. gauge-invariant) states ([ϕ + ], [ϕ − ]) uniquely determines a valid physical state [ϕ] for the field over the entire manifold Σ.
More formally, we label each legitimate/physically possible composition of the two given regional states to form a physically possible universal state by i, with i belonging to some index set I, which can depend on the component states. Thus, formally, failure of GSS corresponds to the set I having either none or more than one element (if I = ∅, the two states are incompatible). As in the more general case of subsystems (as opposed to regions), I will not countenance the possibility that the relation between universe and regional physical states is one-to-many (which would amount to the universe state being a coarsegraining of the conjunction of regional states), nor that there are valid physical states of the universe whose restrictions to regions are not themselves physically valid (I can see no plausible scenario in which that occurs).
Here we will focus on the case in which I has many elements, and, as in the more general nomenclature for subsystems, we will call this Global Non-Supervenience on Regions (GNSS). Thus, denoting the gauge-equivalence class by square brackets, indicating the ith composition of states by ∪ (i) and the resulting universal physical state as [ϕ (i) ], 19 we write or, in terms of the standard Yang-Mills potential and matter fields: So here i is neither a spacetime index nor necessarily related to a Lie-algebra index: it is just an element of an index set I([ϕ ± ])-whose dependence on the given pair [ϕ ± ] will be henceforth omitted-and ∪ S i represents the i-th valid gluing, i.e. composition, of the two gauge-invariant data [ϕ ± ] along S. The global fields ϕ are in the same class of differentiability as the regional ones (albeit the latter will lie on manifolds with boundary).
If I is empty there is no possible gluing, i.e. the regional gauge-invariant states are incompatible and cannot conjoin into a universal physical state (regional incompatibility). If I has a single element, the gluing is unique, and then there is GSS. If otherwise, i.e. if I has more than one element, the universal physical state is undetermined just by the regional physical states: more information about relations between the subsystems is needed, and there is GNSS. In this case, we will also say the universal state possesses residual variety.
Note that as it stands Definition 2 is in line with both the Galileo's ship analogy and with the idea of gauge transformations as mere re-description. The analogy states that, in general, the physical states [ϕ sea ] and [ϕ ship ] can be glued in a variety of ways. Definition 2 is essentially relational: any variety will be solely a variety of relations between the subsystems. In the ship case, this variety is classified by Galilean transformations, i.e. I has a 1-1 correspondence with the Galilean group, as we saw in (1.1). This example illustrates how a transformation taking [ϕ] altering the physical (G-invariant) universal state, could recover DES as defined by Definition 1.
In the following, we will see that there are circumstances in which Yang-Mills subsystems indeed admit transformations with DES due to GNSS in this full sense. Namely, for certain regional gauge-invariant data which can be glued together i.e. composed to give a physically possible universal state, there remains a residual variety of universal gauge-invariant data obtained from this gluing. This variety is parametrized by regional rigid transformations, not regional malleable transformations, and is encoded by the external action of a finite-dimensional Lie group on a subsystem.
In other words, in some cases there is a 'regional rigid symmetries'-worth of universal states which are regionally gauge-equivalent to each other but have a relational physical distinctness, i.e. are physically distinct (not universally gauge-equivalent) due to relational differences. The particular structure of I will be that of a rigid subgroup of the set of gauge transformations, but it will depend on the regional field content in each case. It is the field content which identifies the rigid subgroups of the full infinite-dimensional malleable group that have DES.
In Yang-Mills theory, GNSS will occur only in conjunction with those conditions which are necessary for the existence of conserved global charges, as related to the rigid subgroup I. 20 This procedure thus establishes a link between what is known as an indirect consequence of gauge-the conservation of charges-and a 'direct' one (DES).

Finding GNSS
In this section I will explore the definition of DES as under-determination of universal gauge-invariant data from regional gauge-invariant data. That is, I will explore DES according to GNSS as described in Definition 2. I will illustrate this for Abelian gauge theories (i.e. electromagnetism with a scalar field). 21 For this simple case, I 20 Namely, it will occur only for reducible configurations-those which have stabilizing gauge transformation (analogous to non-trivial Killing fields for a spacetime metric)-in which case I is the group of reducibility parameters (analogous to the isometry group of a spacetime metric) (Barnich & Brandt, 2002). 21 I will only discuss the non-Abelian case in the Appendix. In that case, a concession must be made: due to the non-linear character of the theory, gluing takes place at the perturbative level, and so we must specify which underlying configuration is to suffer the perturbations. (So the index set I would have to be written as I([ϕ ± ]) as envisaged just after equation (2.3)). Nonetheless, the formalism transforms covariantly with respect to gauge transformations of the perturbed configuration and one is able to retain, for the non-Abelian context, all the interesting results obtained in this Section. This sounds like the BRST treatment of gauge theories (see will explicitly show that the whole state is underdetermined by the regional states, and that the ensuing variety of universal states is equivalent to a copy of (i.e. is parametrized by) the charge group (U(1)). I will thus prove my main claim that there is a regional (or subsystem) rigid group of symmetries with (relational) DES, emerging from residual variety, as per Definition 2.
I will organize this section as follows. Section 3.1 will introduce the necessary notation. Section 3.2 sets up the remaining tools for the procedure, giving an introduction to our use of gauge-fixing and gluing. In Section 3.3 I describe precisely how 'external sophistication' is deployed to facilitate the gluing of the regional physical states. In that section, I vindicate my main claims on the connection between supervenience and relational DES: viz, that there is a rigid group of regional symmetries parametrizing by the residual variety of universal physical states which are composed by identical regional physical states.

General Notation
We are given a manifold Σ, which for our illustrative purposes in this Subsection will represent a space(time) endowed with an Euclidean metric. 22 I will also assume Σ is closed, that is, compact and without boundary. Given a charge group G, i.e. the finitedimensional Lie-group characterizing the theory, the group of gauge transformations is G = C ∞ (Σ, G). The gauge field A and its gauge-transformed A g are given by Lie-algebra valued space(time) 1-forms.
In the main text of the paper, I will only consider the simpler case of Maxwell electrodynamics coupled to scalar Klein-Gordon theory. Thus, for the Abelian case, the structure, or charge, group is G = U (1) and we write A g = A + igrad(g) and ψ g = exp(ig)ψ. (3.1) Here A ∈ C ∞ (T * Σ, R) (a smooth Lie-algebra valued 1-form on Σ) and g ∈ C ∞ (Σ, R) = C ∞ (Σ, Lie(G)) =: Lie(G), an infinitesimal gauge transformation. This means I am taking a slight short-cut and working with the Lie-algebra u(1) R as opposed to the group U (1). But this distinction is unimportant for our context -the Abelian caseand purposes:-finding GNSS. 23 The non-Abelian case is treated in the Appendix. I will also assume that the manifold Σ is endowed with a Riemannian metric, g ij , and that it is decomposed into two regions (cf. Figure 1): Σ ± , with boundary ∂Σ ± =: S = Σ + ∩ Σ − . The Σ − piece is what is usually labeled 'the environment', but here will play the role of another subsystem (as for relational DES; cf. footnote 11). For now, each of Σ, Σ ± is assumed topologically trivial; and, if any of these manifolds is not compact, then all the fields on them will be restricted to have some suitably fast fall-off rate. 24 (Henneaux & Teitelboim, 1992) for a review), by which one retains global transformations even if one eliminates the degeneracy in the propagator of the theory through a perturbative gauge-fixing. And indeed, the tools used in this work in the non-Abelian context recover the properties of BRST ghosts; see (Gomes & Riello, 2017) and (Gomes et al., 2019, Sec 3.1) for more on this recovery. 22 The Euclidean metric provides a simpler interpretation of the results on gluing, in the next section. See  for more on this topic, and for how one can import the relevant results to the Lorentzian signature, "3+1" context. In any case, the philosophy of physics literature on DES ignores the contrast between Euclidean and Lorentzian signatures. 23 In the Abelian case we will be mostly concerned with, the relation between the Lie algebra element g ∈ u(1) and the group element, ξ ∈ G, is: g = −i ln ξ, or ξ = exp(ig). This translation can be applied at any point in the following computations. Using the Lie-algebra rather than group is useful in translating our results to the non-Abelian case, since there we cannot work directly with the group. 24 In Figure 1, one could think of Σ − as a collar around Σ + ; this would not block our treatment. The only Figure 1: The two spacetime regions Σ ± , separated by a manifold with codimension one, S.
I will denote the regional, unquotiented configuration spaces of each field sector (gauge field A, matter field ψ, and doublet (gauge and matter fields) ϕ) as A ± , Ψ ± and Φ ± (respectively). I will omit the subscript ± for the corresponding universal configuration spaces (i.e. A, Ψ, Φ, respectively). The restricted groups of gauge transformations (i.e. smooth maps from regions of the manifold into G) will be denoted in analogous fashion: G ± = C ∞ (Σ ± , G), and all abstract quotient spaces are denoted by square brackets, as in [Φ ± ] := Φ ± /G ± , and [Φ] := Φ/G.
If the fields compose smoothly, the left hand sides in the following equation are both smooth fields: Here Θ + and Θ − are the characteristic functions of the regions Σ + and Σ − i.e. they are distributions: unity in the region, and zero outside, with some conventional value at the boundary which is immaterial for our purposes. Smoothness requires equality of the following quantities at S: (∂ n A + ) |S = (∂ n A − ) |S and mutatis mutandis for ψ in place of A; where the superscript n denotes all derivatives: first order, second order, i.e. n = 1, 2, etc.; and for any combination of independent directions, i.e. x, y, z parametrizing Σ. We will understand all equalities at S, written for simplicity like A +|S = A −|S , to be such an equality.

DES as GNSS: set-up
At first sight, the question we face in this Section is much simpler than the previous literature's construals of DES: how do we make sense of the equation (cf. (2.2)): But there are two related obstacles. First, the elements of [A], being orbits of the gauge field space A under the action of G, are not directly parametrizable. Second, the only criterion for gluing quotients employs representatives, i.e. elements of A. That is, there is no composition of physical states that is not formulated in terms of the composition of representatives.
The first step in overcoming these issues is to "gauge-fix": that is, to select a class of elements of A which uniquely represent elements of [A] (cf. (Gomes, 2019a, Sec. 2)). A selective class fixes further conditions which the representatives must satisfy. The second step will consist in exploiting external sophistication, as advocated in Section 1.2.
complication would be to then consider further boundary conditions on Σ − , and so on. We therefore restrict our attention to the case where Σ − encompasses the "rest of the universe", i.e. the entire "environment" in the language of GW.
Given [A ± ] ∈ [A ± ], for any two representatives A ± ∈ A of [A ± ], the condition of composition can then be translated into the following gluing condition: there exist gauge transformations, g ± ∈ G ± , such that the gauge-transformed representatives smoothly join (cf. (3.1)): ( 3.4) If there are also matter fields that are non-zero on S, the corresponding relation to (3.4) is: ψ To ensure that we retain the full physical content of the regional states and the capacity to relate arbitrary configurations, there must be no prior restrictions on the gauge transformations at the boundary. 25 This is crucial, both conceptually-why should the "redundant descriptive fluff" at the boundary be any different than in the bulk?-as well as technically. Our focus on gauge-invariant quantities thus allows a pure gauge discontinuity at the boundary.
In other words, from the perspective of the subsystem-instrinsic informationas in condition 2 of Definition 1-the only criteria for the composition of regional representatives is whether under some appropriate regional gauge transformations the resulting doublets of regional representatives smoothly compose.

Reduction, sophistication, and gluing
This is the most important technical Section of the paper, in which I illustrate the conceptual tools developed so far. In Section 3.3.1, I will show how one reduces, or fixes a selective class of representatives for the gauge fields, through a gauge-fixing. 26 In Section 3.3.2 I will describe how external sophistication is important for gluing. And finally, in Section 3.3.3 I describe the residual physical variety.

The projection h introduced
Given the regions Σ ± , we will consider two states, ϕ ± = (A ± , ψ ± ). We need to uniquely characterize the physical content of such states.
Let us first focus on the gauge-fields, A ± . We will return to the matter fields in 3.3.3. Thus we must fix unique representatives of A ± , through a projection: 27 where h ± A will uniquely represent (i.e. is in the image of all elements of) the equivalence class [A ± ].
25 While such a truncation is standard in the literature concerning gauge theory in asymptotic regions (cf. (Regge & Teitelboim, 1974;Balachandran et al., 1996;Giulini, 1995;Strominger, 2018)), in that context there are no subsystems that should be glued back together, and the environment is not on a par with the subsystem. See (Gomes, 2020) for a treatment of finite regions with such an assumption, there called 'the externalist notion of boundary' and (Riello, 2019) for a recovery of the asymptotic results, including the so-called soft charges, using the present framework. 26 In the Abelian case, the covariance property under transformations of the perturbed configuration need not be flagged explicitly, cf. footnote 27 below. 27 Here h stands for horizontal. Although we will not need to introduce the entire field-space principal fiber bundle formalism here , this is where 'horizontal' comes from: horizontal directions are essentially a choice of non-gauge directions in field space transforming covariantly along the fiber. The word 'horizontal' is more appropriate in the non-Abelian case: in the Abelian, horizontal directions are integrable, and correspond rather to a (covariant) foliation of the field-space A, i.e. to a G-covariant family of gauge-fixings.
But we must maintain our ability to describe the possible gluing of the h ± A through an analogue of (3.4). Thus h ± A must be gauge-invariant intrinsically-reflecting internal reduction-while still allowing gauge transformations to act, "extrinsically", on them-reflecting external sophistication. That is, in the following, I endorse 'reduction' (cf. Section 1.2) for determining the regional and universal physical content. But, to describe gluing, I endorse 'external sophistication': i.e. allowing all the different representations of the same regional physical content to be counted as isomorphic but not identified from the outside view.
Technically, these demands imply we should look for a projection h : A → A, as opposed to a reduction, red : A → [A]. In (Gomes, 2019a(Gomes, , 2020) the construal of a gauge-fixing as a projection h, and not as a quotient, was argued to be fundamental for the gluing of regions: for both h and red are required to be gauge-invariant with respect to gauge-transformations on the common domain, A, i.e. red(A g ) = red(A) as well as h(A g ) = h(A), but only the projection h allows further transformations to be enacted on its range.
First, without loss of generality, we can find linear projections h ± : A ± → A ± such that their images satisfy: div where s is the normal to the boundary S, and · is induced by the inner product on Σ. It might seem surprising that we can restrict our attention to such constrained boundary conditions on A ± and yet still encompass the entire gamut of possible regional physical states. The reason for this is that any other state-including ones with boundary behavior different from (3.6)-differs from such a h ± A by a unique regional gauge transformation; again, a regional gauge transformation that is possibly non-trivial at the boundary.
There are many ways of gluing (as we will see in Section 3.3.2 below), which will give rise to many possible universal representatives of the state. But GNSS refers to the univesal physical state as well as to the regional physical states. Thus, after gluing, we need to resolve the mapping between regional and universal physical states, and so we need to determine physical universal states from glued regional physical states. That is, we need to eliminate the plethora of possible extrinsic gauge transformations as well as the intrinsic ones.
To discern or distinguish the possible universal physical states, we apply the corresponding gauge fixing for closed-compact without boundary-regions. In the absence of boundaries, i.e. for the universal state A, instead of (3.6) we have: and the significance of (3.7) reduces to the familiar statement that the condition div(h A ) = 0; (3.9) is a bona-fide (partial) gauge condition, called the (Euclidean) Lorentz gauge, or, in the non-Abelian setting, the (perturbative) covariant Landau gauge. It is 'partial' because a different choice, related to h by a constant shift, would still satisfy (3.9). It implies that the representative of the universal physical state is only determined up to a global phase shift. This will be important in what is to come. The condition (3.9) will fix the extrinsic gauge transformations, so that we find the sought-for correspondence between regional physical states and universal physical states.
Let us add some detail to this procedure.

Internal reduction and external sophistication: option (c) realized.
We know that each h and h ± yields a unique element in each orbit because of the way the projection h : A → A works by exploiting gauge transformations. For instance, in the global case, for U(1): And a similar σ ± , with analogous properties, exists for for h ± obeying the regional equations (3.7).
Therefore we can now distinguish two sorts of action of G: a subsystem-intrinsic action and a subsystem-extrinsic one. Subsystem-intrinsic transformations will map between the members of the same equivalence class, whereas the extrinsic ones act as transformations between the representatives of these equivalence classes. 28 In more detail, we have: Subsystem-intrinsic gauge transformations Given h ± : A ± → A ± , where domain and range are seen as distinct, but isomorphic, spaces, a subsystem-intrinsic gauge-transformation is just a gauge-transformation acting on the domain of h. It maps between members of the same equivalence class.
The label 'intrinsic' stands in opposition to 'extrinsic'. Intrinsic gauge transformations are all that is needed for a unique description of the entire Universe, since there is nothing external to the entire Universe.
But if we have more than one subsystem and we want to satisfy the gluing condition (3.4), we may need to change the representative of the equivalence class [A ± ]-from the outside, as it were.

Subsystem-extrinsic gauge transformations
We can define subsystem-extrinsic gauge transformations g ext , as those transformations which act on the range of h as Of course such a transformed field would no longer satisfy (3.7).
28 It is instructive to compare the two possibilities of action of G to the use of homotopy type theory (HoTT) in gauge theory, as advocated by (Ladyman, 2015). Ladyman says HoTT "both (a) distinguishes states conceived of differently even if they are subsequently identified, and (b) distinguishes the identity map from non-trivial transformations that nonetheless might be regarded as delivering an identical state". Here we have two sorts of transformations: the subsystem-intrinsic one, A → A g , which does not change h[A]-satisfying Ladyman's (b)-, and the subsystem-extrinsic one, that does the work of Ladyman's (a).
That is, subsystem-intrinsic gauge transformations are defined as those acting on the field configurations in the domain of the projection, whereas the subsystem extrinsic act on its range. Once we have eliminated redundancy and fixed a 1-1 correspondence with [A ± ], the image of h ± , i.e. h ± [A ± ] ⊂ A ± , is invariant with respect to gauge transformations acting on its domain, but we can still change representatives by acting on its range, A ± .
Gluing We are given the physical content of the regional configurations as (in terms of) their projected representatives h ± A , and while these representatives h ± A might not smoothly join, they may still jointly correspond to a physically possible universal state. The existence of subsystem-extrinsic gauge transformations smoothening out the transition between h + A and h − A is a necessary and sufficient condition for their compatibility.
That is, the h ± A determine whether they can be smoothly joined by subsystemextrinsic gauge transformations. Following (3.4), the condition is that subsystemextrinsic gauge transformations g ext ± exist such that (in spacetime index-free notation 29 ): (3.14) which is the appropriate rewriting of the gluing condition (3.4). However, as mentioned, this is not enough to establish a correspondece between universal and (the doublet of) regional physical states: although we have the physical representatives on the regional side, we do not yet have them on the universal side of the correspondence.

Establishing the main claims
Finally, we must know which of these gluings give rise to physically distinct universal configurations. Thus, we need to eliminate the redundancy of subsystem-extrinsic gauge transformations and get a unique representative of the universal state. We first write: where g ext ± obey (3.14), and then we apply the reduction of the universal state, through (3.9), i.e. demanding that div(h A ) = 0 (we assume the universe has no boundary). Imposition of this divergenceless condition on (3.15) almost uniquely fixes the solution: In the U(1) case, there is an under-determination by a constant extrinsic gauge transformation in each region; let us call this indeterminacy (c + , c − ).
is fixed up to the addition of constants, c ± ∈ R. We note that as expressed, this is an indeterminacy at the level of the Lie-algebra; for the Lie-group, we would get exp(ic ± ) ∈ [0, 2π]. 29 Using indices, the equation is: In both Abelian and non-Abelian case, g ext ± depends on h ± A only through their difference at the boundary: (h + A −h − A ) S (in the non-Abelian case, each g ext depends also on its respective regional gauge field, e.g. g ext . For illustration purposes, I display the solution here: where the subscript S denotes operators and quantities intrinsic (i.e. pulled-back) to the boundary surface S (and since normal components of h ± A match, is a harmonic function on (respectively) Σ ± with Neumann boundary condition ∂ n ζ u (±) = u, and R is the Dirichlet-to-Neumann operator. For the meaning of these operators, and also the analogous solution for the general non-Abelian Yang-Mills gauge theories, see (Gomes & Riello, 2019, Sec. 4).
To get to the main claims of the paper, we must now include charged matter fields in our description.
Charged matter fields First, note that the discussion so far focused on the gauge fields, A. But including matter fields is not difficult: since matter fields 'co-rotate' with A ± , they just get "taken for the ride" by the fixing of σ[A]. That is, since the σ ± are unique, a given doublet, ϕ ± = (A ± , ψ ± ) is also projected to a unique representative: where σ ± is the regional version of (3.11). 31 Since we find unique representatives for the full regional physical field content, we can identify (using the notation'≡') the content with the representative: and take these as the starting point for gluing. Inputting (3.17) into the universal state (3.15) and solving for (3.9), we obtain g ext ± [h + , h − ] + c ± . Thus, we have: where the term (+ ↔ −) is identical to the first, mutatis mutandis for + ↔ −, and I have omitted the dependence of g ext ± on h ± A . The important point to understand is that, although the degenaracy in gluing, c ± , has no effect on the gauge field, it will add a phase, or rotate, the matter part.
But now we must consider two cases: either the ψ ± vanish at S, or they don't. Assume first that ψ ± vanish at S. We then have no further constraints and thereby obtain a 2-parameter family of universal states, (3.18), parametrized by exp(ic ± ). In other words, given regional physical states [ϕ ± ], here a conjunction of e.g. a transverse projection of the photons and a Dirac dressing of the charges, we can glue them to form a two-parameter collection of universal states.
Naively, this would give us two copies of U(1), parametrizing the universal physical states that are compatible with the regional physical contents. But of course, if c + ≡ c − , we have a global constant phase shift (which precisely matches the expected leftover under-determination of the universal physical state [ϕ] by (3.9)) which does not change the universal physical state. Therefore, we are left with a residual physical variety parametrized by the difference, c := c + − c − , which is insensitive to any global phase shift. To sum up, if the matter fields vanish at the boundary, we obtain a residual variety that is isomorphic to a single copy of U(1) and in fact can be identified with relative, regional, rigid phase shifts.
Assume now that ψ ± do not vanish at S. Then not every boundary value is allowed: the two sides must differ by a phase at the boundary, and this phase difference needs to match the gauge transformation required to glue the gauge potentials. That is, by (3.5), exp i((g ext In the asymptotic flat case, the universal functional exp (iσ[A])ψ is known as the "Dirac dressed electron" (Dirac, 1955). It is an electron that is "dressed" by an appropriate Coulombic tail, rendering the electron also gauge-invariant (as can be easily checked from the gauge-covariant transformation properties of σ). The function h ± is known as the "radiative" projection of the photon. Here we have extended both notions to finite bounded regions (see (Gomes et  where (g ext + + σ + ) − (− ↔ +) is a fixed functional of the gauge fields, and c ± , which were entirely under-determined previously by the gauge fields A, will now completely fix the quantity (c + − c − ) by (3.19). Thus, if the matter fields do not vanish at the boundary (and are compatible with each other and with the gauge fields, satisfying (3.19)), there is no variety left, since the compatibility equation (3.19) completely fixes the difference c.
Finally, as a corollary of these constructions, we are able to state our main result of this Section for gluing physical states: Theorem 1 (Rigid variety for U(1)) For electromagnetism as coupled to a Klein-Gordon scalar field in a simply-connected universe: given the physical content of two regions, [ϕ ± ], for matter vanishing at the boundary but not in the bulk of the regions, the universal state is underdetermined, resulting in a residual variety parametrized by an element of U(1). In the notation of (2.2) and Definition 2: where the particular action of U(1) is that which leaves the gauge-fields invariant, but not the matter fields. 32 Thus we have found, in the case of electromagnetism coupled to a Klein-Gordon scalar field, a rigid variety through GNSS, as per Definition 2 in Section 2.4. Moreover, by fixing a reference subsystem (Σ − ), we can construe this U(1) as acting as a group of (rigid), regional symmetries over Σ + , that shift the phase of matter but do not affect the gauge potentials (as in the 't Hooft beam-splitter, 33 ('t Hooft, 1980)). That is, we have constant phase shifts.
In direct analogy with Galileo's ship-where we also recognize GNSS as encoded by the (external) action of a finite-dimensional group-we have fully vindicated our main claims: namely, that GNSS is a source of empirically significant subsystem symmetries.
Agreed, the non-Abelian Yang-Mills case is more complicated: non-linearities render the corresponding I of the equation corresponding to (3.20) dependent on the physical state. Nonetheless, our constructions are valid at a perturbative level, i.e. one needs to first fix a ground state and then perturb it. If the perturbed state is the 'vacuum', i.e. [A * ] for A * = 0, then we recover the full Lie algebra of the gauge group through the analogue of (3.20). In that case too, the particular rigid subgroup is identified as the one that leaves the perturbations invariant but not the matter sources.
Section 4.2, I briefly gloss GW's arguments for DES. In Section 4.4, I will show, contra GW, that the lack of DES for general malleable symmetries still allows for the context-dependent identification of some of their rigid subgroups-which do have DES.

Brading and Brown: regional gauge transformations must match
As the quotation at the end of Section 2.3 illustrates, BB take any two regional gauge transformations which fail to coincide at the boundary to be disallowed. In discussing the t'Hooft beam splitter (cf. footnote 33), they conclude: "The only remaining option is to consider a region where the wavefunction can be decomposed into two spatially separated components, and then to apply local gauge transformation to one region (i.e. to the component of wavefunction in that region, along with the electromagnetic potential in the region) and not to the other. But then either the transformation of the electromagnetic potential results in the potential being discontinuous at the boundary between the 'two subsystems', in which case the relative phase relations of the two components are undefined (it is meaningless to ask what the relative phase relations are), or the electromagnetic potential remains continuous, in which case what we have is a special case of a local gauge transformation on the entire system." (p. 656) BB are right in one respect: once you have a universal configuration one can not apply gauge transformations which are discontinuous (e.g. produce a delta function) at the interface. 34 On the other hand, it is also true that one could have regional configurations being acted upon by regional gauge transformations which don't match at the boundary. 35 In this second instance, one aims to take the regional subsystem intrinsically and have non-matching gauge transformations applied to them before gluing.
In other words, in their quotation, BB ignore the fact that the formulation of DES, through property 2 of Definition 1, requires only information intrinsic to a subsystemwhich is why we related it to supervenience of the universe on its subsystems. All hands agree that, from the viewpoint of the universe one cannot have a gauge transformation which is discontinuous at the boundary. From this viewpoint it is true that g + = g |Σ + and g − = g |Σ − ; i.e. that the regional gauge transformations are mere restrictions of a universal gauge transformation. And they are right: such an assumption would pre-empt any search for DES.
But these conclusions are unwarranted. For the topic of supervenience, one starts from the regional states and then composes them. From this perspective, it is the effect of the regional gauge transformations that matter.
GW spot this error, and assert that what should be fundamental is only the continuity of the glued gauge and matter fields, A, ψ, not of phase shifts. In their words (but my notation) (p. 83): The key to seeing why this argument fails is noting that what is given, when we are given the pre and post-transformed states of the universe, is not a function from spacetime to the gauge group, but merely the effect of whatever transformation is being performed on the particular pretransformation (universe) state (ψ, A). And if this particular ψ happens to vanish on the overlap region [the boundary S], then nothing about the corresponding gauge transformations g + , g − can be 'read off' from their effects on the wavefunction in that region [assuming they are constant near S]. It is therefore possible that the universe transformation being performed might correspond to the effect of (say) some constant gauge transformation g + in Σ + , and a different constant gauge transformation g − in Σ − , so that there is no way of patching g + , g − together to obtain a single smooth function from the whole of spacetime to the gauge group. [my italics] I fully agree with this verdict, as far as it goes. It means one considers the effect of the regional g ± on the subsystems from the intrinsic point of view. One does not take g ± as the regional projections of a discontinuous g. But I believe GW are not entirely consistent in applying this approach, as we will shortly see.
Nonetheless, GW are one step closer to an analysis based solely on physical (gaugeinvariat concepts than are BB: they take the representatives of the subsystem fields to be important, not the gauge transformations themselves. I will now develop their view, and, in particular, will parse the italicized text above .

Greaves and Wallace: non-matching gauge transformations encode DES
To recapitulate: GW claim the relational DES transformations are in 1-1 correspondence with the following quotient of two infinite-dimensional groups: where G(s |∂ ) are the gauge transformations of the region which preserve the state ϕ at the boundary of the region, and G Id are the gauge transformations of the region which are the identity at the boundary.
Here is the gist of the argument (reflected in the above quote) leading from Definition 1 to (4.1) (and reflected in the above quote): certain gauge transformations are not the identity at the boundary and yet they may keep particular boundary states invariant. One can use such regional transformations to obtain a different universal state, but the composition of the regional gauge transformations is not itself smooth and therefore does not count as a gauge transformation relating the initial and final universal configurations. GW's mistake is that they implement no criteria to establish whether the regional and universal states are indeed physically distinct. As I will show below, this omission allows us to find a simple counter-example, in which (4.1) is non-trivial and yet the transformations constructed above fail to yield physically distinct states, and thus do not satisfy condition 1 of Definition 1.
In more detail, let (A ± , ψ ± ) be two configurations, one in each region Σ ± , that join smoothly, and such that ψ ± = 0. Therefore we have the initial universal configuration As a condition for these representatives to smoothly compose, we must have: A +|S = A −|S (where equality here includes equality of derivatives at S). But A +|S has a stabilizer: at the Lie algebra level, this is g +|S = c + = 0 such that grad(g + ) |S = 0 and where c + is some constant on S. Therefore, the configuratioñ is still smooth, since g + doesn't change the value of A + at the boundary. But with g − = 0, 36 is not a smooth (infinitesimal) gauge transformation (because g + doesn't vanish at S). GW would conclude from this thatÃ is not a gauge-transformation of A, and that, therefore,Ã and A are physically distinct. Moreover, they observe, such a g + corresponds to an element of (4.1) and essentially, they claim, the same construction would apply for any other element of this quotient group.
But without dealing with gauge-invariant quantities, GW have no warrant to conclude thatÃ and A are physically distinct, and therefore no warrant to conclude that condition 1 of Definition 1 is satisfied. In fact, it is only true thatÃ as it is written, doesn't appear to be a gauge transformation of A. But we can explicitly construct a gauge transformation relatingÃ and A as follows: Letg + := g + − c + . Now Ag + + = A g + + (because the constant part has a trivial action on A + ), and so: But now g :=g + Θ + + g − Θ − =g + Θ + is a smooth gauge transformation (sinceg + vanishes at S). And thereforeÃ = A g . So [Ã] = [A]: that is, the physical, gauge-invariant, states are identical and the transformation cannot have empirical significance (condition 1 of Definition 1 fails). 37 (Although the above argument was formulated explicitly for electromagnetism, its extension to the more general cases is straightforward.)

GW and BB's error
In my view, neither GW nor BB could have obtained the right characterization of DES, for Definition 2's characterization of DES in terms of physical i.e. gauge-invariant, states was not articulated by either group. Indeed, they both explicitly endorse GSS. On this topic, GW write: GSS': "Firstly, in doing so we make the assumption that knowing the state of the subsystem and its environment suffices to specify the state of the total system. [...] For example, it is true for Yang-Mills gauge theories in the connection formalism but not in the holonomy formalism." (p. 67) This assumption, which I have labeled as an alternative statement of GSS, should not be confused with its converse, which, all hands agree, holds for all of these systems. Namely, the converse assumption-that the state of the whole, uniquely determines the states of the subsystems-is not under dispute since the state of the whole includes all relational information and more. But as we have stressed: given just the intrinsic physical states of the subsystems, it is not necessarily the case that there is just one way of putting them together; their conjunction may lack necessary relational information.
The exception GW make for the holonomy formalism is telling here: how can physical significance depend on the choice of variables? Indeed, for us it is immaterial: for electromagnetism, it is possible, using holonomies, to recover precisely the same 36 That is, GW take g − ≡ 0: they take the environment gauge transformations to be the identity, or to serve 'as a reference' for the gauge transformations of the subsystem, and the original g + to be the identity. This is a bit confusing, since stipulating the value of g does not usually fix a selective class: given A, one cannot use this condition to assess whether A belongs to the selective class. Here, since we are only interested in a counter-example to substantiate our criticism of their construal of DES, we will ignore this point (which is not an issue for our formulation through GNSS). 37 If there was matter in the bulk, our counter-example would fail, because ψg + + = ψ g+ + . But the source of such DES would be just what I described here through Theorem 1: namely, as required from the Theorem, there is a stabilizer of the gauge potential but not of the matter field. results as in Theorem 1. For GW, it fails precisely because the holonomy formalism deals with gauge-invariant observables. 38

How can a rigid subgroup have DES?
Now it is our turn to defend our construction from one of GW's arguments, who judge it impossible that a 'global' (i.e. rigid) gauge transformations, being a subset of the 'local' (i.e. malleable), may acquire DES, while none of the malleable do.
As expressed in equation (4.1) for relational DES, GW only claim a group isomorphism between G GW DES , whatever it may be, and a quotient group, "empirical symmetries correspond 1-1 [...] to elements of a quotient group" (p. 75). 39 But GW do not see the quotient nature of their result as problematic. In fact, they see it as exonerating their notion of relational DES from the charges they make against the orthodox view on DES: In any theory that has a malleable symmetry group, the rigid symmetries remain as a subgroup of that malleable symmetry group. (For example, in general relativity, the rigid translations and boosts form a subgroup of the group of all diffeomorphisms.) It is therefore logically impossible that all rigid symmetries, but no malleable symmetries, can have direct empirical significance; (p. 61) (again the quotation is adapted to my nomenclature). That is, their accusation against the orthodox view is that it attributes relational DES only to certain rigid subgroups of malleable groups. Recognizing the quotient group as being isomorphic (in 1-1 relation) to the symmetries exhibiting DES, GW argue, assuages this concern. 40 38 Here I comment on the relation between GNSS and Myrvold's 'global patchy-non-separability' (Myrvold, 2010, p.427), which they articulated for electromagnetism in holonomy variables (cf. footnote 5). In electromagnetism, given the space of loops (smooth embeddings γ : S 1 → Σ), one can form a basis of gauge-invariant quantities by the holonomies, exp (i γ A) (this can be accomplished more generally for non-Abelian theories using Wilson loops, (Barrett, 1991)). For simply connected regions Σ like ours, by composing regional loops γ ± ∈ Σ ± going in opposite directions at the boundary S it is true, as Myrvold argues, that we recover the gauge-invariant holonomy corresponding to a larger loop γ not contained in either region. Therefore any universal holonomy corresponds to a single doublet of regional holonomies. According to Myrvold, separability fails only for non-simply connected manifolds, where the holonomies of γ ± cannot recover the universal holonomy of γ.
I have two comments to make on the relation to the present work: (i) In the absence of matter and in the Abelian setting, non-trivial topology indeed is the only source of GNSS, and we recover Myrvold's conclusions. But with matter, we can close off curves which are not loops-which only change by gauge transformations at their ends-by placing charges to cap off the curves, thereby obtaining gauge-invariant holonomies (cf. (Gomes & Riello, 2019, Sec. 4.3.2)). These sourced holonomies correspond to the residual variety we found for electromagnetism in Theorem 1, which occurs even for simply-connected Σ when charged matter is present in the regions Σ ± .
(ii) Unfortunately, Myrvold's loop composition doesn't work in the same way for the non-Abelian theory: although the appropriate regional loops themselves will compose as curves in the manifold, Wilson loops-giving the gauge-invariant content of the holonomies-involve traces, and the traces ruin the composition properties: the corresponding regional gauge-invariant quantities do not compose. Our construction (see appendix A.2) gets around that.
The conclusions, in Myrvold's nomenclature, are then: global patchy separability fails for non-Abelian theory if and only if: the manifold is non-simply connected, or charged matter is present inside the regions and furthermore this charged matter obeys (perturbative) regional conservation laws (i.e. and the perturbed gauge-field has stabilizers). 39 The claim is based on intuitions for gauge theories in the asymptotic regime. But the situation in the asymptotic regime, or for truncated configuration spaces, is more complicated: there, one does not consider environment and subsystem to be on a par as we do for relational DES; the environment truncates the fields and the gauge transformations at the boundary of the subsystem. In that case, we can obtain a group of symmetries with DES that is isomorphic to the quotient and has no action on the subsystem states. These results, obtained in (Gomes, 2020, Sec 4.2.2), largely recover those presented by GW (erroneously) for relational DES in (4.1), and more.
40 " The reason lies in a structural difference between the original problematic claims and our replacements: But is it really true that one cannot endow significance solely to a subgroup of a group? What are the rigid translations and boosts of a generic spacetime metric? Poincaré transformations are not well-defined (i.e. geometrically defined) in a generic background metric. They are defined by Killing fields for a Minkowski metric. In that sense, that subgroup is physically distinguished, but only in particular circumstances: it is only meaningful for a Minkowski metric; 41 in specific backgrounds (e.g. Minkowski), there will be a physically well-defined rigid subgroup of the malleable transformations, but this subgroup will be effaced once one moves to generic backgrounds.
Admittedly: if one focuses just on the group itself, and not on its action on states, one indeed cannot "pluck out" a rigid subgroup in any meaningful way. But contra GW, it is entirely possible to associate DES only to certain physically meaningful rigid symmetries; the meaning is acquired through their action on the fields. The transformations are the only ones that leave the gauge-field, but not the matter fields, invariant. 42 That means that in any background, malleable transformations that don't belong to these subgroups would lack DES. That is, rigid symmetries may have DES in a given special background and yet lose that significance for a generic background. Such subgroups are "plucked out" from the surrounding malleable group by satisfying certain equations, e.g. the Killing equations, which are themselves physically significant.
Such subgroups of rigid symmetries are usually called stabilizers, and the configurations they stabilize are called reducible (cf. footnotes 17 and 20). And the same concepts apply to gauge theories. In the Abelian case all A ∈ A are reducible, and they all possess the same stabilizers, namely, the constant gauge transformation. The non-Abelian case is much more similar to the spacetime case (see A.3 for a brief description): generic A are not reducible, but some are. Being reducible is a physically significant fact: in the quotient space [A], the orbits of reducible configurations are qualitatively different than the generic orbits. 43 rather than holding that all elements of a (malleable) symmetry group have one property while elements of some subgroup thereof have a contradictory property, we hold that (for any given subsystem) there is a subgroup of subsystem symmetry transformations (the 'interior' ones) that cannot have empirical significance, and that it is elements of the quotient of the larger group by this subgroup that are candidates for correspondence to physical operations. There is thus no object of which we assert both that it does, and that it does not, have some given property." (p. 87, my italics). A rigid subgroup corresponds to physical operations through its action on the configuration variables (e.g. boosts on the ship, or global changes of phase). But how the quotient is to do the same is not specified. And it is easy to see that if one insists only on the 1-1 relation between the quotient group and the group with DES-i.e. the one which could act as physical operations on a subsystem-there is no unique way of realizing the said physical operations (as there is with Galileo's ship); even if the quotient is finite-dimensional, it will have a continuous infinity of equivalent representations on the region. 41 If one is thinking not in terms of active diffeomorphisms, but of coordinate transformations, then indeed, one can single out translations and boosts, but only with respect to that coordinate system. More broadly, there are generically no constant gauge-transformations: they usually require a global section to be defined. One should also note that although for principal fiber bundles one has a natural action of the charge group, G, this is not the case for associated bundles. One cannot define a "constant" action of the gauge group: it can only be constant with respect to a given section (see (Kobayashi & Nomizu, 1963)).
42 Indeed, such special configurations are part of the structure of the physical quotient of the configuration space of the theory by the gauge transformations, Φ/G: the quotient is not a manifold, but a stratified manifold, and the orbits whose fields are stabilized by subgroup of G are the 'strata '. 43 Note that the notion of reducibility covaries with the notion of stabilizers. Namely, if a configuration ϕ is reducible, with a given stabilizer f ∈ G, then for any given g ∈ G, ϕ g will also be stabilized by gf g −1 . This qualitative difference between the orbits renders the quotient space into a stratified manifold: i.e. a space formed by a concatenation of boundaries (see (Kondracki & Rogulski, 1983;Fischer, 1970) and, for a philosopher-friendly description of stabilizers, and their relation to conserved charges, (Gomes, 2019a, Secs. 3.3.5-7))

Summary
Broadly speaking, in this paper I have explored the role of GNSS in the context of gauge theories. I dissected the meaning and occurrence of "empirically significant subsystem symmetries", whose existence and characteristics are still matters of debate in the philosophy of physics. To make matters concrete, I have focused on Yang-Mills theory (and, in more detail, electromagnetism) and restricted subsystems to be demarcated by spacetime regions.
In order to better adjudicate the debate between the opposing sides-represented by Greaves and Wallace (GW) on one side, and Brading and Brown (BB) on the other-it was first necessary to clear the ground by introducing new nomenclature. The standard nomenclature of local and global gauge symmetries is perfectly adequate if the system under study is the entire Universe, as is usually the case. But if one wants to discuss subsystems, and needs to distinguish between symmetries acting also at these different levels, the standard nomenclature becomes awkward. The awkwardness is apparent when one refers, e.g.: to 'a global subsystem symmetry'. The new nomenclature uses four labels to express two logically independent distinctions: namely (i) whether symmetries act only regionally or universally, and (ii) whether they depend on an infinite or a finite set of parameters, i.e. whether they are rigid or malleable when acting throughout a spacetime region or throughout the universe. It also disentangles possible confusions with 'non-local and local' functions.
With these definitions in place, and in the specific setting of Yang-Mills theory in bounded regions, I found myself in agreement with GW in their criticism of BB. Namely, BB prematurely dismissed the possibility of DES tranformations by assuming that such a transformation would always create discontinuities in the boundary between the regions. The difference between GW and BB can be briefly summarized thus: one should be concerned with smooth composition of regional gauge and matter fields (GW), as opposed to a smooth composition of gauge transformations (BB), which are not physical.
But the gluing of the gauge and matter fields is also not physical; only the composition of their gauge equivalence classes is. This misunderstanding is reflected in GW's explicit assumption of supervenience on regions (cf. GSS' in Section 4.3), that is, that regional states should uniquely define the universal state. GW admit this condition fails for gauge theories in certain gauge-invariant bases (holonomies for the Abelian case), but it even clearly fails for Galileo's ship. Thus GSS seems inconsistent with an appropriate definition of DES for gauge theories.
My interpretation of GW and BB's oversight is that they forget aspects of the Definition 1 of DES: they do not check whether condition 1 is satisfied explicitly (by employing gauge-invariance information) and they do not consider condition 2 as involving physical information that is intrinsic to a subsystem or region.
Having noticed these issues, we could summarize the journey from the present results to GW and from there to BB as follows: one should be concerned with the smooth composition of the physical data, as given by gauge-equivalence classes, and not with smooth composition of gauge and matter fields through given representatives (GW), and much less with the smooth composition of gauge transformations (BB).
Finding the criteria for smooth composition of the information contained in the equivalence classes [ϕ ± ] still presented a challenge. The challenge was: can their physical content be composed into a universal physically valid state?
The way we got around the problem of gluing physical content was to use an 'exter-nally sophisticated view of symmetries'. 44 Namely, I employed 'reduction', or better, 'projection', as a means to identify the regional and universal physical contentseach non-locally determined within its corresponding domain,-and then I employed sophistication for gluing. We found that for electromagnetism as coupled to a scalar Klein-Gordon fieldwhen the Klein Gordon fields were taken to vanish at the boundary in between the regions, but not in the bulk within them-there can be multiple universal physical states formed by gluing the same regional physical states. That is, a failure of uniqueness in the gluing creates the gap which gives a residual variety of universal physical states. In other words, a variety of universal physical states can be built from the same regional physical states. These regional and universal states fulfill our Definition 2 for GNSS, and thereby also correspond to purely relational DES, as described in Definition 1.
Even though I have used 'external sophistication', Theorem 1, including its significance for DES, is completely compatible with viewing gauge degrees of freedom as "descriptive fluff". The fact that the theory admits such a particular sort of redundancy is related to the particular sort of non-locality of its gauge-independent degrees of freedom.
In other words, gauge-invariant quantities are to some extent non-local (L9 in Earman's classificatory scheme (Earman, 1987)), which means there is in principle "room to explore" between the whole and the sum of the parts. And it turns out that the known aspects of non-locality in gauge theory match the non-local aspects of GNSS.
One can see this is as follows. Gauss constraints are the defining characteristic of gauge theories: they stipulate that charges couple to fields in such a way that their conservation laws are dynamically respected. But Gauss's law can clearly fluster cluster decomposition: once one knows the electric flux around a closed surface, one knows precisely the amount of charge within it. Quantities measured on the totality of the boundary are therefore not independent of other quantities measured in the bulk. Of course, this type of "synchronic non-locality" is not causal; it only represents "non-locally possessed" properties (in the language of (Belot, 1998;Healey, 2007)). That is, there should be room for holism, which this paper has drawn on.
The work of  reported here precisely delineates the sort of non-locality involved. For generic regional states in Yang-Mills theory, one does not have GNSS: in most circumstances we can describe the physical state of the whole by describing the physical state of its composing regions. When the regional states do not uniquely determine the universal state, GNSS ensues; and with it, DES as I have defined it. That is, it is a particular non-locality that is responsible for the gap between the regional gauge-invariant information and the universal gauge-invariant information; it is this gap from which DES emerges.
Moreover, as expected, the group of symmetries with DES has intimate connections with the charge group of the theory: For both Abelian and non-Abelian theories the residual variety of universal physical states can be parametrized by sub-algebras of Lie(G). And these rigid symmetries with DES are in 1-1 correspondence with the possible (covariantly) conserved charges in the region (cf. (Gomes & Riello, 2019, Sec. 4.3.2)). In the Galileo's ship example, the variety is given by the action of the Galilean group, 'Boosts Euclidean,' on the subsystem, and the respective charges are the conserved momenta.
What these findings suggest is that (for simply connected regions, 45 and at the perturbative level in the case of the non-Abelian theory) the only failure of holism of gauge theories is precisely that related to the total covariantly conserved charges within any region. That is, the only failure is due to the Gauss law.
In the landscape of the debate, we thus locate ourselves somewhere between the two opposing views. On the side of BB and the orthodoxy, 46 we find that indeed only rigid, but no malleable symmetries may have direct empirical significance in gauge theory. In fact, surprisingly, the rigid ones correspond to 'the global' gauge transformations, G (for non-Abelian: only for very particular perturbed configurations, that can carry associated covariantly conserved charges). On the side of GW, we find that gauge theories may indeed harbor symmetries that have direct empirical significance. But to force my conclusion into either pigeonhole, or even a combination, would be to shave off some of the important subtleties of this situation.
In a few words: for the physical, i.e. gauge-invariant, content of the gauge fields, the whole can be more than the sum of the parts. Their difference manifests itself in a direct empirical significance of (subgroups of) the charge group, and are in 1-1 relation to the conserved charges of the theory. The same can be said for Galileo's ship scenario, where the difference manifests itself in the direct empirical significance of the Galilean group.

The role of the subsystem
The conclusion we have arrived at says that DES are not inextricable from other effects of symmetries. In particular, I have related direct empirical significance (DES) to a failure of supervenience on subsystems (GNSS), and then stated that GNSS is related to the indirect significance of symmetry (labeled IES), like conservation of charges.
Given the sort of holism we have seen in GNSS, it thus might be conceptually preferable to use global information about charges and other superselected quantities pertaining to IES to parametrize the possible gluings of regions in a gauge-free manner.
In that case, we may not want to invoke external sophistication.
Certainly in practice we will often lack the relevant global information, and so it makes sense to use a local gauge-fixed parametrization to do calculations, make predictions, etc. Similarly, insofar as we only have a perturbative handle on the gluing procedure in the non-Abelian case, and our current perturbation theory relies on compactly localized fields, there is reason to use gauge-fixed compactly-localized fields instead of gauge-invariant non-compactly localized fields, and external sophistication gives a conceptual underpinning to the ensuing gluing procedures. But these sorts of concerns are probably best seen as practical or technical barriers, not as in-principle reasons to adopt external sophistication.
Indeed, from a global perspective external sophistication loses its warrant. But then, so does the empirical relation between conserved charges and rigid symmetries, since DES necessarily involves an external perspective. 45 For both Abelian and non-Abelian Yang-Mills, for non-simply connected manifolds we get something beyond G: namely, we also get the first homology group, whose dimension is finite and given by the so-called first (equivariant) Betti number. Such a contribution from non-trivial topology was to be expected from the source of our results: holism. 46 In fact, in this point, we roughtly align with Kosso (Kosso, 2000)-but not always: only when covariantly conserved charges exist-and not with Brown and Brading, who want to claim that even rigid internal symmetries have no DES.

Acknowledgements
I would like to thank Aldo Riello, my collaborator in all things gauge-related, for many fruitful and insightful discussions and for co-developing the entire technical apparatus employed in this analysis. And similarly, I would like to thank Jeremy Butterfield, for many discussions, for many requests for clarification, for patiently guiding me in all things philosophy-related, and for giving immensely valuable and thorough feedback on this paper.

APPENDIX A Non-Abelian Yang-Mills
A.1 Notation for Yang-Mills theory coupled to matter The general case works with a finite-dimensional charge group, e.g. G = SU(N ), with Lie algebra g := Lie(G), e.g. g = su(N ). Given the charge group we define the group of gauge transformations G = C ∞ (Σ, G), with composition given by pointwise action of G, i.e. (gg )(x) := g(x)g (x), and the respective infinitesimal version, Lie(G) = C ∞ (Σ, g) (with pointwise Lie algebra commutator). An element of G is a map: g(·) : Σ x → g(x) ∈ G. The gauge fields and its gauge-transformed version A g are given by Lie-algebra valued space(time) 1-forms.
We define, for dx µ a basis of 1-forms over Σ and τ α a basis of the Lie algebra g Infinitesimally, i.e. for an infinitesimal gauge-transformation ξ ∈ Lie(G), the gauge field transforms as For full generality we can introduce charged fermions in a fundamental representation: so for some vector space W , the 4-component Dirac spinor field (i.e. in C 4 ) as: ψ ∈ C ∞ (Σ, C 4 ⊗ W ) and ψ g = gψ.
And I will write the joint configuration as ϕ = (A, ψ) I will denote the regional, unquotiented configuration spaces of each field sector (gauge field and matter, respectively) as And, for the joint configuration spaces, I write Φ ± = {ϕ ± := (A ± , ψ ± ) ∈ A ± × Ψ ± }, writing Φ, Ψ, A. I will omit the subscript ± for the corresponding universal configuration spaces. The restricted groups of gauge transformations will be denoted in analogous fashion: G ± = C ∞ (Σ ± , G), and all abstract quotient spaces are denoted by the square brackets, as in [Φ ± ] := Φ ± /G ± , and [Φ] := Φ/G.
We can solve these emerging conditions (see Section 4 in ), obtaining: where, in each region, ζ φ stands for a covariantly harmonic function, satisfying D 2 ζ = 0, with Neumann boundary conditions at S given by s i D i ζ |S = φ, and where the subscript S denotes "intrinsic to S", and the intrinsic coordinates to S are given by c, and where R is the so-called Dirichlet-to-Neumman operator. Briefly, R functions as follows: a given harmonic function with Dirichlet conditions-these conditions are the input of R-will possess a certain normal derivative at the boundary; i.e. will induce certain Neumann conditions there-these conditions are the output of R. That is, let ζ u be a harmonic function with Neumann boundary condition, then for ζ u a covariantly harmonic function with Dirichlet boundary condition, (ζ u ) |S = u, the Dirichlet to Neumann operator R is defined as ζ u = ζ R(u) (i.e. it finds the harmonic function with Neumann condition corresponding to one with a Dirichlet condition).

A.3 Non-Abelian GNSS
There are two sources of under-determination of solutions to (A.5): one topological and one from homogeneous fields. The first arises from the possibly non-trivial first (equivariant) homology group of Σ, we look at the first in Section A.3.1 and at the second in Section A.3.2.

A.3.1 Topological holism
We may have a universal field δA satifying (A.4) which, when restricted to each Σ ± is pure gauge, i.e. of the form Dξ ± . In other words, the room for discrepancy is equivalent to broken homology cycles. This is easier to see in the Abelian case, where δA would be a one form which is exact in the simply connected patches, but only closed in the entire manifold (see , Section 4.7 for an example). By the Poincaré lemma, this occurs if and only if the topology of the manifold is nontrivial. In other words, universal physical processes in this case may not come from regional physical processes; some physical processes are universal/global in nature. But of course for topologically non-trivial manifolds we should expect the whole to contain more information than the sum of the parts, foiling Antiholism! The arising GNSS would contribute with the suitable homology group to the variety set I in (2.2). Such a contribution has a finite number of generators (given by the first Betti number, i.e. the rank of the first homology group). 47

A.3.2 Stabilizer ambiguity
Barring non-trivial topology, is there another source of under-determination of ξ ± also in the non-Abelian case? To arrive at an answer, note that we are using the gauge-fields as a reference. It is through them that we fix the ξ and parametrize the reduced/quotient configuration space [Φ] = Φ/G. Therefore one source of GNSS would be a possible under-determination of ξ ± .
This under-determination occurs if and only if there exist χ ± ∈ Lie(G) such that Dχ ± = ∂χ ± + [A ± , χ ± ] = 0 (A.7) for A = A + Θ + + A − Θ − being the smooth field configuration around which we are considering perturbations. 48 This is, again, the "Killing" (or infinitesimal stabilizer) equation for gauge transformations. Namely, under-determination of ξ ± occurs if and only if there exist χ ± which stabilize the entire A ± , in which case we could use either ξ ± or ξ ± + χ ± for gluing. Such stabilizers are generalizations of the constant potential shifts in electromagnetism. Here the entire regional field must be stabilized by χ ± , and this singles out from G + the regional elements which can potentially exhibit DES.
As with Killing directions for spacetime geometries, such elements are hard to come by: they're generically trivial (for non-Abelian charge groups G), and, when non-trivial, they are generated by a finite-dimensional basis.
In sum: there may still be some small variety in ξ ± , which is innocuous as far as the gauge field is concerned. However, as we will shortly see, this variety can be relationally felt by regional matter fields.

A.3.3 Obtaining variety from the stabilizer ambiguity
When we include matter fields, the universal field is composed as: δψ = (δψ + + ξ + ψ + )Θ + + (δψ − + ξ − ψ) (A.8) Since we are not using the matter fields to parametrize the quotient/moduli/reduced configuration spaces, there are no analogues of equations (3.7) and (3.9) to be implemented for it. But smoothness of the field does impose certain conditions. In fact, different choices of stabilizer, adding χ + = χ − to ξ + , ξ − , respectively, will render the matter fields incompatible at the boundary. 49 Therefore, if the matter fields do not vanish at the boundary, we can assume there is again regional determinacy as per Definition 2. So we assume they vanish there; a condition postulated by GW as "dynamical isolation" of the two regions. And also, if the matter field vanishes not only at the boundary, but everywhere, no GNSS, i.e. GSS arises, yet again. Thus, I assume: the matter fields vanish at the boundary, but not in the bulk of the regions.
To illustrate the precise emergence of the Lie algebra g, we resort to a 'best-case' scenario. That is, apart from the conditions stipulated above for the matter fields, we will take the gauge field configuration around which we are perturbing to be the 'vacuum', A = 0. In that case, the gauge covariant differential just becomes the standard differential (as in the Abelian, electromagnetic case), i.e. D → ∂, and the stabilizer equation becomes (∂χ α (x))τ α = 0. Then not only is the space of stabilizers χ a finite-dimensional vector space, closed under commutation, but for this case it forms a Lie algebra homomorphic to g.
Finally, we have arrived at our destination: I have shown that the variety set I, at least infinitesimally and in the best-case scenario, can in some sense recover the charge group G, i.e. I have characterized the infinitesimal version of (3.3), with: But in the Abelian case, the tools utilized here do allow an integration to the finite setting. 50 Moving away from the best-case scenario of A = 0: to obtain GNSS the base configuration A must still be somewhat homogenous; it must admit internal directions which leave it unchanged. But as mentioned above, a general theorem about Killing fields (and stabilizers), shows that Killing directions are always generated by a finite dimensional basis; in the case of internal gauge transformations, this basis consists of at most dim(g) elements, closed under commutation, etc. Thus even away from the best-case, we obtain I's isomorphic to sub-algebras of g (we would also have to replace τ α in the definition of δψ (i) for a choice of basis {χ α } of the appropriate Killing, or stabilizing fields).
Note that, as had to be the case, according to Definition 2, there exists a regional variety implicit in equation (A.10), and this variety recovers also the definition of DES through transformations, if the transformations referred to in Definition 1 are taken to be e.g.: (δA, δψ) → (δA, δψ (i) ).