Abstract 1 Introduction 2 Preliminary Definitions and Models 3 Similarities between the aTAM and the syncTAM 4 Separating the aTAM and the syncTAM by Shape Production 5 Separating the aTAM and the syncTAM by 𝝉=𝟏 Computation References Appendix A Technical Appendix

Synchronous Versus Asynchronous Tile-Based Self-Assembly

Florent Becker ORCID Laboratoire d’Informatique Fondamentale d’Orléans (UR4022), Université d’Orléans, France Phillip Drake ORCID University of Arkansas, Fayetteville, AR, USA Matthew J. Patitz ORCID University of Arkansas, Fayetteville, AR, USA Trent A. Rogers ORCID University of Arkansas, Fayetteville, AR, USA
Abstract

In this paper we study the relationship between mathematical models of tile-based self-assembly which differ in terms of the synchronicity of tile additions. In the standard abstract Tile Assembly Model (aTAM), each step of assembly consists of a single tile being added to an assembly. At any given time, each location on the perimeter of an assembly to which a tile can legally bind is called a frontier location, and for each step of assembly one frontier location is randomly selected and a tile is added. In the Synchronous Tile Assembly Model (syncTAM), at each step of assembly every frontier location simultaneously receives a tile.

Our results show that while directed, non-cooperative syncTAM systems are capable of universal computation (while directed, non-cooperative aTAM systems are known not to be), and they are capable of building shapes that can’t be built within the aTAM, the non-cooperative aTAM is also capable of building shapes that can’t be built within the syncTAM even cooperatively. We show a variety of results that demonstrate the similarities and differences between these two models.

Keywords and phrases:
self-assembly, noncooperative self-assembly, models of computation, tile assembly systems
Funding:
Florent Becker: Research made at a public, state-funded university, thanks to the guarantee of full academic freedom and perennial research funding.
Phillip Drake: This author’s work was supported in part by NSF grant 2329908.
Matthew J. Patitz: This author’s work was supported in part by NSF grant 2329908.
Copyright and License:
[Uncaptioned image] © Florent Becker, Phillip Drake, Matthew J. Patitz, and Trent A. Rogers; licensed under Creative Commons License CC-BY 4.0
2012 ACM Subject Classification:
Theory of computation Models of computation
Editors:
Josie Schaeffer and Fei Zhang

1 Introduction

At the macroscale, humans have mastered the construction of structures with desired spatial and temporal properties. This mastery has provided us with not only entertainment, but also the tools necessary for our survival. As civilization shifts towards a more virtual world and the threats to our survival evolve, a mastery of matter at the macroscale is no longer sufficient. Unfortunately, controlled interaction with the universe at the nanoscale presents significant hurdles. We face major difficulties in observing the universe at the nanoscale, let alone manipulating matter at the this scale. Just as axiomatic models (such as the natural or real numbers) have greatly aided in our understanding (and manipulating) of the universe at the macroscale, axiomatic models of self-assembly have proven successful in allowing us to understand (and manipulate) the universe at the nanoscale. Indeed, such axiomatic models have already enabled the design and implementation of a variety of futuristic nanotechnologies which have a wide range of uses: molecular computers [28], drug delivery [26], and lithographically-patterned electronic/photonic devices [13].

Here, we focus on tile-based, axiomatic models of self-assembly. In these models, the atomic components are square “tiles” which have a “glue” on each of their edges. Growth begins from an initial “seed” assembly and proceeds via single tile attachments. On one end of the spectrum, these single tile attachments occur asynchronously and nondeterministically in a model called the abstract Tile Assembly Model (aTAM). And, on the other end of the spectrum, these single tile attachments occur completely synchronously so that every “frontier” location receives a tile at every step. The fully synchronous model has yet to be studied, so in this paper we introduce a formalization of this model called the synchronous tile assembly model (syncTAM). The syncTAM provides a formal framework for studying the power and limitations of fully-synchronous systems and gives us insight into how we can leverage these synchronization mechanisms to imbue our experimental systems with more power than their asynchronous counterparts.

Perhaps surprisingly, the aTAM gives rise to a variety of complex behaviors. For example, the aTAM is known to be computationally universal [27], intrinsically universal [6], and capable of assembling complex, fractal structures [2]. But, it also has its limitations. For example, the aTAM is limited in the types of interesting, infinite shapes it can assemble [2] and the class of directed, temperature-1 aTAM systems is incapable of computing [19]. Here we investigate whether the syncTAM can overcome some of the limitations of the aTAM and whether a fully-synchronous attachment regime introduces its own set of limitations.

Current ubiquitous technology is capable of physically realizing approximations of asynchronous models of self-assembly, such as the aTAM, but it is unclear if a fully synchronous model can be physically realized in an efficient manner. Inefficient implementations of more synchronous models currently exist. For example, in [25], experimentalists demonstrated the ability to implement a “staged” system where a target structure is assembled in stages by assembling different parts of the structure in separate test tubes and then combining the results to obtain a new substructure which feeds into the next stage. This staging can be thought of as a synchronization mechanism since (if we wait long enough) every assembly will be terminal before it is combined with others. To implement a fully synchronous model using this approach a stage is required for every step making this approach so tedious and laborious it is infeasible.

This work seeks to answer the following question: Is exploring ways to efficiently implement a more synchronous model (than the existing models which are able to be physically realized) a worthy endeavor? We seek to answer this question by comparing the relative powers of an asynchronous model of self-assembly (the aTAM) to a synchronous model of self-assembly (the syncTAM) in terms of shape-building power and computational power.

While the asynchronous mode of tile attachment is the default for self-assembly literature, the synchronous mode of update is the default update model for cellular automata. Fatès, Régnault, Schabanel and others [9, 10, 23, 5] have shown that introducing asynchronism in even the simplest cellular automata vastly alters their behavior. Likewise, Paulevé and Sené [22] show the importance of the choice of update mode in boolean networks. Due to its synchronous nature and the fact that a tile attachment depends on the surrounding neighbors, the class of syncTAM systems can be thought of as a restricted class of freezing cellular automata (CA). Indeed, the class of freezing CA only allows states to increase (according to some ordering), so restricting this class to only allow a single state change and a radius-1 Von Neumann neighborhood yields a class which roughly corresponds to syncTAM systems – it is a slightly more permissive superclass of the syncTAM. Consequently, the results from these models are transitive which provides an interesting mathematical motivation for the study of SyncTAM systems.

Just as comparing asynchronous and synchronous models has been a central theme in cellular automata, our results seek to compare the powers and limitations of the syncTAM with respect to the aTAM. To this end, we begin with formally defining the models used in this paper in Section 2, as well as the formal notions we use to compare the models. Section 3 states the elementary connections between the behavior of the same system in the aTAM versus the syncTAM. In Section 4, we show that there exists shapes that can be assembled in the temperature-1 syncTAM which cannot be assembled at any scale factor in the aTAM (at any temperature). Likewise, that section exposes a set of shapes which can be assembled nondeterministically by an aTAM system but cannot be assembled at any scale factor by a nondeterministic syncTAM system. In Section 5, we show that like the directed aTAM, the directed syncTAM is computationally universal, but unlike the directed, temperature-1 aTAM [19], the syncTAM is computationally universal at temperature 1. To show this, we show that the any temperature-2 compact zig-zag system in the aTAM, can be simulated by a temperature-1 syncTAM system. And, since the class of aTAM temperature-2 compact zig-zag systems is computationally universal [1], the class of temperature-1 syncTAM systems is computationally universal. While previous results about computing at temperature 1 in the aTAM or adjacent models led to very poorly connected assemblies [16, 11, 12, 21, 4], we also show in Section 5.2 that the syncTAM can compute using more well-connected structures.

To make it easier to understand and visualize our results, we have created a web-based simulator for the syncTAM [15] and have implemented all of the systems used in our results [8]. These can all be accessed from http://self-assembly.net/wiki/index.php/Synchronous_Tile_Assembly_Model_(syncTAM).

2 Preliminary Definitions and Models

In this section we define the models and terminology used throughout the paper.

2.1 The abstract Tile-Assembly Model

We work within the abstract Tile-Assembly Model[27] in 2 and 3 dimensions. The abbreviation aTAM refers to the 2D model. These definitions are borrowed from [14] and we note that [24] and [17] are good introductions to the model for unfamiliar readers.

Let be the set of nonnegative integers, for l,u, let l,u={klk<u}. For n, let [n]=0,n={0,1,,n2,n1}.

Fix Σ to be some alphabet with Σ its finite strings. A glue gΣ× consists of a finite string label and non-negative integer strength. There is a single glue of strength 0, referred to as the null glue. A tile type is a tuple t(Σ×)4, thought of as a unit square with a glue on each side. A tile set is a finite set of tile types. We always assume a finite set of tile types, but allow an infinite number of copies of each tile type to occupy locations in the 2 lattice, each called a tile. A glue set GT is the set of all glues associated with a tile set T. Given a direction d{N,E,S,W}=𝒟 and a tile set T, GT,d is the set of all glues on the edges of tiles in direction D and gt,d is the glue g on tile t in direction D. We write str(d)t to denote the strength of gt,d.

Given a tile set T, a configuration is an arrangement (possibly empty) of tiles in the lattice 2, i.e. a partial function α:2T. Two adjacent tiles in a configuration interact, or are bound or attached, if the glues on their abutting sides are equal (in both label and strength) and have positive strength. Each configuration α induces a binding graph Bα whose vertices are those points occupied by tiles, with an edge of weight s between two vertices if the corresponding tiles interact with strength s. An assembly is a configuration whose domain (as a graph) is connected and non-empty. The shape of X2 is {X+(x,y)|x,y}, i.e. the set of all its translations. The shape of an assembly α is the shape Sα of its domain: two assemblies have the same shape if they have the same domain modulo translation. An assembly α contains a shape S if there is sS such that sdomα. For some τ+, an assembly α is τ-stable if every cut of Bα has weight at least τ, i.e. a τ-stable assembly cannot be split into two pieces without separating bound tiles whose shared glues have cumulative strength τ. Given two assemblies α,β, we say α is a subassembly of β (denoted αβ) if domαdomβ and for all pdomα, α(p)=β(p) (i.e., they have tiles of the same types in all locations of domα).

A tile-assembly system (TAS) is a triple 𝒯=(T,σ,τ), where T is a tile set, σ is a finite τ-stable assembly called the seed assembly, and τ+ is called the binding threshold (a.k.a. temperature). If the seed σ consists of a single tile, i.e. |σ|=1, we say 𝒯 is singly-seeded. Given a TAS 𝒯=(T,σ,τ) and two τ-stable assemblies α and β, we say that α 𝒯-produces β in one step (written α1𝒯β) if αβ and |SβSα|=1. That is, α1𝒯β if β differs from α by the addition of a single tile. The 𝒯-frontier is the set 𝒯α=α1𝒯βSβSα of locations in which a tile could τ-stably attach to α. When 𝒯 is clear from context we simply refer to these as the frontier locations.

We use 𝒜T to denote the set of all assemblies of tiles in tile set T. Given a TAS 𝒯=(T,σ,τ), a sequence of k+{} assemblies α0,α1, over 𝒜T is called a 𝒯-assembly sequence if, for all 1i<k, αi11𝒯αi. The result limα of an assembly sequence α is the unique limiting assembly of the sequence. For finite assembly sequences, this is the final assembly; whereas for infinite assembly sequences, this is the assembly consisting of all tiles from any assembly in the sequence. We say that α 𝒯-produces β (denoted α𝒯β) if there is a 𝒯-assembly sequence starting with α whose result is β. We say α is 𝒯-producible if σ𝒯α and write 𝒜[𝒯] to denote the set of 𝒯-producible assemblies. We say α is 𝒯-terminal if α is τ-stable and there exists no assembly that is 𝒯-producible from α. We denote the set of 𝒯-producible and 𝒯-terminal assemblies by 𝒜[𝒯]. If |𝒜[𝒯]|=1, i.e., there is exactly one terminal assembly, we say that 𝒯 is directed. When 𝒯 is clear from context, we may omit 𝒯 from notation.

Given TAS 𝒯, we define 𝒮[𝒯]={Sαα𝒜[𝒯]} as the set of all shapes of terminal assemblies in 𝒯. Clearly, if 𝒯 is directed |𝒮[𝒯]|=1. However, even if 𝒯 is not directed it may be the case that |𝒮[𝒯]|=1 (see [3] for an example). In such a case we call the system shape directed.

Given S, a connected set of points in 2, we define a version of S scaled by factor c, and denote it by cS, as cS={(x,y)2(xc,yc)X}. Intuitively, cS is a version of S expanded by replacing each point of S by a c×c square of points. Likewise, for a set A of shapes, cA is defined as {cSSA}.

2.1.1 Zig-zag Tile Assembly Systems

We now provide the definition for a subset of aTAM systems which are known as compact zig-zag systems (originally defined in [4]). These systems have very simple dynamics (the tiles have uniform input and output sides) yet are capable of universal computation – for every Turing machine, there exists a compact zig-zag system which simulates it [4]. Due to their simple dynamics and Turing Universality, these systems are an ideal target for systems with simple mechanisms to simulate to demonstrate they are Turing Universality. This approach has become a paradigm in axiomatic, tile-based self-assembly [4, 16, 11, 12, 21], and we follow this same approach in Section 5.

The following definitions of zig-zag, and compact zig-zag systems are almost identical to those described in [21]. 𝒯=(T,σ,τ) is a zig-zag tile assembly system provided that (1)𝒯 is directed, (2) there is a single 𝒯-assembly sequence α𝒯, with 𝒜[𝒯]={α}, and (3) no tile initially binds to an assembly with its north glue. More intuitively, a zig-zag tile assembly system is a system which grows to the left or right, grows up some amount, and then continues growth again to the left or right. Again, as defined in [20], a tile assembly system 𝒯=(T,σ,τ) is a compact zig-zag tile assembly system iff 𝒜[α]={α} and for every x dom α, str(N)α(x) + str(S)α(x)<2τ and str(W)α(x) + str(E)α(x)<2τ. Informally, this can be thought of as a zig-zag tile assembly system where two τ-strength tile attachments in a row cannot occur in the same direction. This means a row can never extend by more than one tile over the previous row, and the system is only able to travel upwards one tile at a time before being require to zig-zag again.

2.2 The Synchronous Tile Assembly Model

The Synchronous Tile Assembly Model (syncTAM) is the variant of the aTAM where at every step of assembly, rather that a single location of the frontier being randomly selected to receive a tile, all frontier locations simultaneously receive a tile each. As the frontier of a system can grow arbitrarily large, the number of tiles added per assembly step may also grow arbitrarily large. As in the aTAM, if any single frontier location allows for τ-strength binding of tiles of multiple types, one of those types is non-deterministically chosen for each such location.

We formally define the syncTAM in the same way as the aTAM with the exception of the definition of a single step in the assembly sequence. Give a syncTAS 𝒮=(T,σ,τ) and two τ-stable assemblies α and β, we say that α 𝒮-produces β in one step, written α1𝒮β, if αβ and dom(β)dom(α)=𝒮(α). Note that like in the aTAM, we define 𝒮(α) to be the set of locations where a tile can τ-stably attach to α. The syncTAM inherits all other definitions and notation from the aTAM.

To distinguish between a tile assembly system (a.k.a. TAS) in the aTAM versus one in the syncTAM, we will use the term syncTAS to refer to a TAS in the syncTAM.

3 Similarities between the aTAM and the syncTAM

In this section we show that directed and shape-directed aTAM systems behave equivalently in the syncTAM. That is, given a tile set T, seed σ, and temperature τ, if (T,σ,τ) is used as an aTAM system and that system is directed (resp. shape directed), then when (T,σ,τ) is used as a syncTAM system the same unique terminal assembly (resp. assembly shape) is produced. Although these are relatively straightforward results, they demonstrate the types of systems, namely those with reduced nondeterminism, which are invariant to the differing dynamics of the models.

Lemma 1.

Let 𝒮=(T,σ,τ) be a syncTAM system and 𝒯=(T,σ,τ) be an aTAM system composed of the same components. If α and α are assemblies of 𝒮 such that α𝒮1α, then α𝒯α.

Proof.

Let 𝒮=(T,σ,τ) be a syncTAM system and 𝒯=(T,σ,τ) be an aTAM system composed of the same components. Furthermore, assume that α and α are assemblies of 𝒮 such that α𝒮1α.

Let F be the set of pairs consisting of a frontier location of α and a tile that can τ-stably attach to α at that location such that for every pair (p,t)F, α(p)=t. Intuitively, F is the set of all tiles (a tile type and position) such that adding the tiles of F to α yields α. Our construction of F ensures that all the positions contained in the pairs are distinct (that is, no two pairs contain the same position), and every tile can attach τ stably to α. Consequently, we can construct a valid 𝒯-assembly sequence from α to α through single tile attachments from the tiles in F. Indeed, let (fi)i=1n be an enumeration of the set F. Define α0=α and αi=αi1fi. Then α0α1αn is a valid assembly sequence since all tiles can τ-stably attach and none prevent the others from binding, and αn=α.

Section 3 shows that each single step of syncTAM assembly can be matched by a set of steps of an equivalent aTAM system. Thus, beginning from the same seed assembly, successive applications of that lemma show that any assembly that is producible in a syncTAM system is also producible in the equivalent aTAM system.

Lemma 2.

Let 𝒯=(T,σ,τ) be an arbitrary directed (resp. shape directed) aTAM system. Then the syncTAM system 𝒮=(T,σ,τ), composed of the same components, is also directed (resp. shape directed) and terminally produces exactly the same singular terminal assembly (resp. singular shape for all terminal assemblies) as 𝒯.

Proof.

Let 𝒯=(T,σ,τ) be an arbitrary directed aTAM system, and let 𝒮=(T,σ,τ) be composed of the same components. By definition, 𝒯 has exactly one terminal assembly which we denote by α. We claim that this is also the only terminal assembly of 𝒮. Suppose, for the sake of contradiction, this was not the case. That is, suppose α𝒜[𝒮] and αα. Let αi be a 𝒯-assembly sequence such that α0=σ and limαi=α and let αi be a 𝒮-assembly sequence such that α0=σ and limαi=α. Then it follows from Lemma 3 that αi is a valid 𝒯 assembly sequence which contradicts our assumption 𝒯 is directed.

The argument for the shape directed version of the lemma is similar, but instead of assuming (for the sake of contradiction) that αα as above, we assume for the sake of contradiction that SαSα. The rest of the argument remains the same.

Thus, the two lemmas presented here show that systems with limited nondeterminism produce the same assemblies (and/or shapes) in both the aTAM and the syncTAM. This completely changes in the general case: the additional power, and constraints, imposed by the global synchronization of the syncTAM also provide for separations between the models, which we now present in Sections 4 and 5. This culminates in Theorem 8 which demonstrates a shape S such that no aTAM system can assemble S, but there exists a syncTAM system which assembles S.

4 Separating the aTAM and the syncTAM by Shape Production

The difference between the aTAM and the syncTAM is manifest in the shapes an sets of shape each model is able to assemble. This section exhibits a set of shapes and a shape which are exclusive to one of the models, respectively the aTAM and the syncTAM. Much like how we ignore constants when studying the runtime of algorithms, often in tile-based self-assembly the size of a “pixel” in the shape is discounted. In our results, it will not just be the shape which is unobtainable in the “wrong” model, but all of its scaled versions. Likewise, while in each case, the “right” model can assemble those shapes without collaboration (i.e. at temperature 1), no amount of collaboration can help the wrong model in obtaining them.

4.1 Sets of Shapes that need the Asynchronism of the aTAM

The synchronism of the syncTAM can act as a constraint, making it impossible for a single syncTAS to match the dynamics of the more asynchronous aTAM. Namely, we exhibit a simple aTAM system that can create an infinite set of different shapes, while any syncTAS must be constrained to only forming a subset of them.

Impossibility proofs for the aTAM often rely on the Window Movie Lemma [20] (an overview of which can be found in Section A.2). Alas, it does not hold in the syncTAM, as the delay between successive attachments along the window can be used by either side to avoid being “fooled” by the change in the other. On the other hand, in the syncTAM, positions cannot remain attachable for an arbitrarily long time, which prevents the delaying tactics the aTAM can use. Hence, the fooling lemmas for the syncTAM differ from those of the aTAM. The following two lemmas capture this and will prove sufficient for the purpose of this paper.

The first limit of syncTAM systems is that they cannot leave positions attachable for an infinitely long time (and this holds for any infinite assembly sequence, which differs from the aTAM).

Lemma 3 (Completion).

Let 𝒮 be a syncTAM system. For any infinite assembly sequence α of 𝒮, limα is a terminal production: limα𝒜[𝒮].

Proof.

Assume z is a frontier location in limα. There is a subset N of the neighbors of z such that a tile could attach τ-stably to the tiles of limα within N. Let t be the last time a tile was attached in N. Then z belongs to the frontier of αt, has been filled at time t+1 and cannot be empty in limα.

Thus, the frontier of limα is empty and limα is a terminal production.

The second limitation of syncTAM systems is that if their productions can avoid some patterns for an arbitrarily long time, then they have a terminal production avoiding those patterns.

Lemma 4 (Compacity).

Let F be a set of finite shapes of 2. Let 𝒮 be a syncTAM system with, for each k, a production Pk𝒜[𝒮] producible in k assembly steps that does not contain F.

Then there is a terminal production T𝒜[𝒮] such that dom(T) does not contain any element of F.

Proof.

The attachment sequence (βi)i0 of 𝒮 (defined recursively as follows) is such that from each βi, infinitely many Pi can be 𝒮-produced:

  • β0 is the seed assembly, from which all the Pi can be produced,

  • for each i, since βi is a finite production, it has only finitely many descendants (productions d such that βi1𝒮d). Since βi 𝒮-produces infinitely productions Pj, one of its descendants must also produce infinitely many Pj; pick that descendant to be βi+1.

By Section 4.1, T=limβ is a terminal production, yet it cannot contain any element of F. If it did contain an XF, since X is finite, there would be a finite f such that βf contains X and that βf would not be compatible with any of the Pk.

Theorem 5.

There is an aTAM system 𝒜 such that for any syncTAM system 𝒮 and scale factor c1, 𝒮[𝒮]c𝒮[𝒜].

Proof.

Let 𝒜 be the aTAM system of Figure 1. As depicted in Figure 1, let A(d,r,u) for d0,r0,u4 be the subset of 2 defined by:

A(d,r,u)=(0,r×0)(0×d,0)(0,min(u+4,4)×{d})({3}×d,ud)
Figure 1: Left, the tile types of the temperature 1 aTAM System 𝒜. There is a line between two tile edges whenever their glue match, edges without incident lines have a 0-strength glue. The seed is the red s tile. Each of the tiles r,d,u can grow arbitrarily long tile segments, respectively to the right, upwards and downwards. To the right, a typical production A(u,d,r) of 𝒜, then the three sets of terminal assemblies of 𝒜: AΓ, Ad(d) and Aσ(d). The yellow parts at the top of the productions are the copies of Γ, the gray parts at the bottom, the copies of U4, and the encircled part at the top is the copy of B5.

The shape of any assembly of 𝒜 is of the form A(d,r,u) for some values of d,r and u. The quantity u starts from negative values so as to correspond to the number of u tiles of the assembly, as shown on Figure 1. The “3rd, 2nd, 1st and 0th u tiles” are the 4 tiles at the bottom of the assembly between the d and u arms. The values of d,r and u are not independent:

  • when d=0, u=4: the bottom part can only grow after the lower arm has grown

  • only one of u>d and r>3 can hold at once, because of the competition between the two arms for position (3,0).

The shapes of the terminal productions of 𝒜 are of one of the following forms:

AΓ = A(+,+,3)
Ad(d) = A(d,3,+)
Aσ(d) = A(d,+,d)
 Remark 6.

Having glanced at Figure 1, it may seem straightforward that 𝒜 cannot be simulated by a syncTAS 𝒮, assuming 𝒮 proceeds like A, with two concurrent arms starting from the top-left: the horizontal arm cannot wait for a downwards arm which may never be coming back to the meeting point; yet, it can never turn to the north without knowing that the other arm will come back. The remainder of the proof essentially makes that argument robust against 𝒮 placing its seed elsewhere, using its scale factor to synchronize the two arms, or any other shenanigans.  

Back to the proof, consider the following shapes:

  • Γ={(0,0),(0,1),(1,0)},

  • U4={(0,1),(3,1)}({0}×0,3), a 4-tile wide U,

  • B5=(0,4×{0}), a 5-tile horizontal bar.

Any terminal production of 𝒜 must contain Γ, and either B5 or a translated copy of U4.

Assume now there is a syncTAM system 𝒮 and an integer c such that 𝒮[𝒮]=c𝒮[𝒜]. Then for any terminal production T𝒮 of 𝒮, there is a terminal production T𝒜 of 𝒜 and a vector τ such that domT𝒮=cdomT𝒜+τ. Note that while c is the same for all productions, τ may be different for each of them.

All terminal productions of 𝒮 must have a translated copy of cΓ, and a translated copy of either cU4 or cB5. By compacity (Section 4.1), there is an integer t({cΓ}) such that any production at time tt({cΓ}) has (cΓ)+τ as a subset of its shape for some vector τ2. Likewise, there is an integer t({cU4,cB5}) such that any production at time tt({U4,B5}) has either (cU4)+τ or (cB5)+τ as a subset of its shape for some vector τ2.

Let T=max(t({cΓ}),t({cU4,cB5})), the shape of any production P of 𝒮 at time T contains: a copy of cΓ, and a copy of either cU4 or cB5. Because P was produced in time T, the distance between its copy of cΓ and its copy of cB5 or cU4 is less than 2T+1.

Let d=2T+1c+1. The terminal production Ad(d) of 𝒜 contains neither a copy of B5, nor a copy of U4 within distance 2T+1c of its copy of Γ: its unique Γ is at (0,0), while the leftmost tile of its U4 is at (0,d). Hence, there is no production of 𝒮 at time larger than T whose shape is contained in cAd(d). Therefore, there is no final (infinite) production of 𝒮 with shape cAd(d), which is a contradiction.

4.2 A Shape Demanding the Synchronization of the syncTAM to Self-assemble

This section exhibits a simple (infinite) shape which cannot be uniquely self-assembled in the asynchronous aTAM, while it can be assembled at temperature 1 by a relatively simple syncTAS. Essentially, the synchronous nature of the syncTAM allows for a (shape) directed syncTAS that creates a single shape using an infinite series of rows of tiles growing from opposite directions that are always guaranteed to meet near a central location, while the asynchronous aTAM cannot guarantee that both sides will always grow at the same speed and therefore the meeting locations may differ and the aTAM system must also be able to produce other, non-target shapes.

Consider the syncTAM system 𝒮 of Figure 2(a), and let V be the following subset of 2, as depicted on that figure:

V={(0,0)}LRk>0(BkTkBk+Tk+)

where:

the left branch

is the set L={(x,x),(x1,x)|x},

the right branch

is the set R={(x,x),(x+1,x)|x},

the left bar at height k

is the set Bk={(x,4k)|x4k+1,1}, and

the right bar at height k

is the set Bk+={(x,4k+1)|x0,4k}

the left thumbs

form the set Tk={(2x+1,4k+1)|x2k,0}

the right thumbs

form the set Tk+={(2x,4k)|x0,2k}

(a) The tiles of the temperature 1 syncTAM system 𝒮 (the seed is the yellow s tile), and a portion of the infinite shape V assembled by 𝒮. Growth begins from the seed tile at the bottom and proceeds with the diagonal “arms” growing upward to the left and right, and with red rows growing right from the left arm, and blue rows growing left from the right arm, at every 4th step. The green part is L, the orange and brown part is R, the blue parts are the Bk+Tk+, while the red parts are the BkTk.
(b) Using pumping in the aTAM on the two windows depicted at the top to produce an illegal production at the bottom: the right part has grown past the middle.
Figure 2: A syncTAM system assembling a shape that cannot be obtained in the aTAM.
 Remark 7.

The system 𝒮 is directed and V is the shape of its unique terminal assembly.  

Theorem 8.

There is no TAS system which uniquely assembles V at any scale factor: for all TAS systems 𝒜 and c1, 𝒮[𝒜]c𝒮[𝒮].

Proof.

Consider an integer c and a TAS system 𝒜 with g different glues. Given a producible assembly A in 𝒜 and a position z2, we say that A has entered p if c{z}domA.

As in [18], define a glue movie on a set of edges W of the grid to be the order of placement, position and glue type for each glue that appears along the window W in an assembly sequence. There are at most p=2(4c)!(4c)g possible glue movies of 𝒜 on a window of size 4c. Indeed, there are (4c)g potential scenarios for the tuple of glue types that lie along one side of the window, and there are (4c)! permutations of this tuple, and we get the 2 in the beginning of the expression from the fact that we must consider both sides of the windows. Assume without loss of generality that p>4 and thus p is even.

Consider an assembly sequence (α)i=0 of 𝒜 reaching a terminal production; the shape of limα is cV translated by a vector τ2. From now on, all coordinates are translated by τ, so that dom(limα)=cV.

Let t be the minimal t such that there are p+3 different positions of B4p+2 where αt has entered. Likewise, t+ is the time where p+3 different positions of B4p+2+ have been entered by αt+. Assume for now that t+<t, and let A=αt+.

Pose A0=A, and x0max=c(4p+2)2. Let Wx be a rectangle of height 4c, width p+1 with its lower-left corner at (x,(16p+7)c. In A0, on any window Wx with xx0maxp1,x0max there are no tiles touching the left, upper and lower edges of Wx, while there must be at least a tile inside Wx. Thus, there are only p different window movies possible for all windows Wx, while there are p+1 such windows. There must be a,b such that the Window Movie Lemma of [18] applies between Wa and Wb. This begets a production A1 in which at most p+3 positions have been entered in B4p+2, while p+4 positions have been entered in B4p+2+. This process can be iterated by posing ximax=x0maxi, until i=3cp. In A3cp, the right bar from the right extends past the midpoint, and the resulting shape is not a subset of V. We have obtained a terminal shape that is not included in V and a contradiction.

If t<t+, the reasoning is the same, with x varying at iteration i between ximin and ximin+p+1, with ximin=ximax. Then A3p extends straight past the midpoint.

5 Separating the aTAM and the syncTAM by 𝝉=𝟏 Computation

In this section we show that directed, temperature-1 syncTAM systems are computationally universal (while this is not the case for directed, temperature-1 aTAM systems [19]). To prove this, we show that any temperature-2, compact zig-zag aTAM system can be simulated by a temperature-1 syncTAM system. The result then follows from a previous result: For every Turing machine M, there exists a temperature-2 compact zig-zag aTAM system which simulates it [4]. This approach (often used to show a class of temperature-1 systems are computationally universal [16, 11, 12, 21, 4]) benefits from the fact that rigorous definitions are already in place defining simulation between the necessary models. We end the section by showcasing a class of computationally universal temperature-1 syncTAM systems wherein each row is connected to the next by two north-south attachments, a level of connectivity not seen in previous simulations of TMs by non-cooperative tile assembly systems.

5.1 Computation at Temperature 1 in the syncTAM

We first give a high-level overview of how a temperature-2, compact zig-zag aTAM system 𝒯 simulates a Turing machine M. At a high level, 𝒯 simulates M by exposing glues on the north of assembly which encode the symbols on the tape of M. One time step of a TM is simulated by growing a row which “zigs” (a row which grows to the right) and then a row which “zags” (a row which grows to the left). If the head transitions to the right during the step, this occurs during the growth of the “zig” row. A left transition, is performed during the growth of the “zag” row. In either case, the cooperative placement of tiles (that is, the requirement that a tile binds with two strength-1 glues) allows for information about the head of the state to be propagated through the east and west glues while information about the symbols stored in the tape is propagated through the south and north glues. See part (a) of Figure 5 for an example which highlights the movement of the head as an aTAM system simulates a Turing machine. For a more in depth explanation of how compact zig-zag systems simulate Turing machines, see [4].

The generally used technique (employed in [16, 11, 12, 21, 4]) is to get one of the two necessary inputs of a simulated cooperative binding via a glue attachment, and the other via carefully regulated growth around some sort of hindrance provided by geometry or glue repulsion. The subsets of tiles that perform the simulation of a single cooperative binding are referred to as bit-readers and bit-writers, i.e. tiles that grow a portion of an assembly capable of logically reading or writing bit values.

Theorem 9.

Let 𝒯=(T,σ,2) be an arbitrary compact zig-zag aTAM system. There exists a syncTAS system 𝒮=(S,σ,1) such that 𝒮 simulates 𝒯 at scale factor O(log|GT,N|)

Refer to caption
Figure 3: An example of writing and reading a bit in the syncTAM. The red tiles represent the bit-writers. (a) A value of 0 was first written by the red tiles. Later, growth from the west causes placement of the green tiles. The rightmost green tile initiates the growth of the yellow and blue paths. Since the yellow path is shorter it always arrives at the second yellow location first, blocking placement of a blue tile later. The yellow path proceeds, encoding that a value of 0 was read. (b) A value of 1 is encoded in the geometry of the red tiles. The northernmost red tile blocks the growth of the yellow path allowing the blue path to complete instead, encoding that a value of 1 was read.

The system 𝒮 simulates 𝒯 by assembling n×m blocks of tiles, which we call macrotiles. These macrotiles map to tiles in T under some macrotile representation function, R. We design 𝒮 so that it behaves exactly like 𝒯 under the macrotile representation function.

It is straightforward to design 𝒮 to simulate the strength-2 attachments in 𝒯, but simulating cooperative binding attachments (attachments where a tile binds with two strength-1 glues) in 𝒯 using 𝒮 is much less obvious. Suppose that we want to simulate the attachment of a tile which binds cooperatively via its west and south glues. At a high-level, our construction uses gadgets to encode the northern strength-1 glues on a preceding row as geometry. Information about the strength-1 east glues being simulated by the macrotile to the west is presented as a glue on the macrotile boundary. From this glue, gadgets assemble which “read” the existing geometry below. Details about how these gadgets read and write the geometry is shown in Figure 3. The end result of this is the exposure of a glue which is unique to the simulated glue from below and the simulated glue from the preceding tile in the row. Consequently, this glue allows a tile to attach that causes the macrotile to map to the correct tile t under the representation function, and it allows for a chain of bit-writers which encode information about the north glue of t to grow into the neighboring macrotile region to the north. After the bit-writers are assembled, a path of tiles grow to expose a glue corresponding to the west glue of t on the boundary of the next macrotile region.

We now provide a sketch of the proof.

Proof.

Let 𝒯=(T,σ,2) be a compact zig-zag TAS. We now describe how to construct a syncTAS 𝒮=(S,σs,1) which simulates 𝒯 (under the definition of simulation formally defined in Section A.1). Note that under the definition of simulation, the domain of a macrotile is square. Our construction is designed to work with a notion of simulation that allows for rectangular macrotile domains, but we could design our tile set to pad the lengths of the paths that it grows to account for square macrotile domains.

We construct 𝒮 so that the macrotiles over 𝒮 consist of translations of subassemblies, also referred to as gadgets, called bit-writers, bit-readers and translators. We define 0-bit-writers and 1-bit-writers to be any translation of the red subassemblies shown in Figure 3. We refer to these collectively as just bit-writers. A bit-writer encodes (in binary) the label of a strength 1, north glue of a tile in 𝒯 so that it can be “read” by a bit-reader contained in the macrotile region above it which will enable 𝒮 to determine which tile the macrotile should map to in 𝒯. Gadgets called 0-bit-readers and 1-bit-readers “read” this information (accomplished through growing exactly one of two paths depending on the geometry). We define these to be translations of the non-red subassemblies shown in Figure 3. Note that each of these gadgets grow from the same tile, but the existing geometry of the bit-writers (which were guaranteed to have previously completed growth) give rise to exactly one of the two subassemblies growing. A translator gadget is a translation of a subassembly that maps to a tile, denote this tile t, which binds with a strength 2 glue. This subassembly is simply a path of tiles that grows from the “beginning” of a macrotile to the “beginning” of the next macrotile. The exact beginning and end position of this path depends on the location of the strength 2 glue on t and the locations of its output glue. In the case that the tile to which a translator gadget belongs contains a strength-1 glue to its north, the translator may contain a chain of bit-writers in order to express this glue.

We now describe the idea behind bit-writers and bit-readers. At a high level, bit-readers attempt to grow two paths of tiles, say a 0-path and a 1-path. The bit-writers are designed to block exactly one of these paths, leaving the “surviving” path to constitute a read of its relevant bit. To accomplish this behavior, we design the 1-path to be exactly one tile longer than the 0-path, so if a previously assembled 0-bit-writer does not block the 0-path, then it will be one tile “ahead” of the 1-path, and thus block the 1-path, hence, a 0 is read. However, if the 1-bit-writer blocks the 0-path, then the 1-path will not be blocked, and a 1 a is read. Figure 3 shows the details of this mechanism. The glues in these surviving paths then “know” the bit that was encoded in the geometry by the bit-writer. It can use this knowledge for subsequent tasks such as assembling its own bit-writers.

We now describe the desired behavior of 𝒮, and then we describe how we can construct the tile set S to obtain this desired behavior. We describe the behavior of 𝒮 in two scenarios. The first scenario we consider is the scenario where 𝒮 needs to assemble a macrotile that maps to a tile (under the macrotile representation function) that attaches cooperatively. The second scenario we consider is the scenario where 𝒮 is growing a macrotile that “simulates” the attachment of a tile in 𝒯 via a strength-2 attachment.

In the following discussion, let t=((x,y),t), for tT, be a tile (i.e. a placed instance of a tile type) in a producible assembly of 𝒯, and let REG(x,y) be the macrotile region which maps to the tile at location (x,y) (in 𝒯) under the macrotile representation function. That is REG(x,y) is the set {(i,j)|(mxim(x+1)) and myjm(y+1)}. Then, intuitively, REG(x,y) is the region which maps to t, REG(x1,y) is the region which maps to the tile to the west of t, call this tile tw, and REG(x,y1) is the region which maps to the tile to the south of t, call this tile ts.

In the first case, assume t is a tile in a producible assembly of 𝒯 which binds into the assembly cooperatively using two strength-1 glues. Without loss of generality, we assume the two “input” glues t uses to bind to the assembly are on its west and south. The case where they are on its east and south is symmetric. We now describe how 𝒮 leverages bit-writers and bit-readers to assemble a macrotile in the macrotile region R that maps to t under the macrotile representation function. Before a bit-reader begins assembling in a macrotile region, our construction of 𝒮 ensures a chain of bit-writers has assembled in the southern region of the macrotile. The bit-writers are chained together such that the last tile in a bit-writer exposes a glue on its east so that another bit-writer can begin growing from this glue. This chain of bit-writers is assembled during the growth of the macrotile to the south, and it encodes the label of the north glue of ts in binary. For each strength 1 glue on the west of a tile in T, say gw, there is a unique chain of bit-readers that grows from a glue that corresponds to gw. So, a chain of bit-readers specific to the east glue of tw starts assembling from the southwest most corner of REG(x,y). Each bit-reader in this chain is unique to the binary sequence read up to that point, so that by the time the bit-reading chain is assembled, there is a tile placed at the end of the bit-reading chain that is specific to the east glue of tw and the north glue of ts. It is at this point that the macrotile representation function begins mapping the subassembly contained in REG(x,y) to t. A new path of tiles begins assembling from this glue which grows bit-writers in the macrotile region to its north, that is the region R(x,y+1), encoding the north output glue of t. Then a path grows from the end of this chain of bit-writers to place a glue on the border of the macrotile region directly to the west of R, that is the region R(x+1,y), which is specific to the east output glue of t. The process then repeats.

Note that log|GT,N| bits are necessary to express the label of a strength-1 north glue in 𝒯. As such, each macrotile that simulates a tile that has a strength-1 glue on its south must include log|GT,N| bit-readers, and each macrotile that expresses a strength-1 north glue must grow log|GT,N| bit writers. Therefore, with a bit-reader/writer length of 4, our macrotiles are of scale 4×O(log|GT,N|).

The second case we consider is where t is a tile in a producible assembly of 𝒯 that binds into the assembly via a strength-2 glue. In this case, the simulator simply grows a path of tiles (which are unique to that glue and, hence, uniquely determine the tile in T the macrotile needs to map to under the representation function) to place a glue at the beginning of the next macrotile region so that the growth process can repeat. This path may include a chain of bit-writers, in the case that the tile exhibits a strength-1 glue to its north. The exact path that is grown depends on the exact locations of the input and output glues on t.

Compact zig-zag systems are singly-seeded (consequently their seed tile exposes a strength-2 glue if growth proceeds) by definition, so to form the seed for our simulator, we simply place a tile so that it exposes a glue that begins the simulation process described above.

We now describe the tiles in S that lead to the desired behavior described above. For each east/west glue in T, we add tiles to S with unique glues that assemble a chain of bit-readers. Note that the glues for each bit-reader in the chain are distinct and we add tiles for each position in the chain which assemble a 0-bit-reader and a 1-bit-reader. In addition, for each tile in T, we add unique tiles to S that grow the bit-writer corresponding to the label of the strength-1 north glue (if it exists) and grows a path to the macrotile boundary exposing a glue corresponding to the east/west glue. An example macrotile, including both a chain of bit-readers and bit-writers is displayed in Figure 4(c). Finally, for each translator we add tiles to S that assemble the path of tiles composing the translator and have unique glues that do not interact with any other gadgets in 𝒮. The tiles required to construct the translators of an example compact zig-zag system are displayed in Figure 4(b).

Refer to caption
(a)
Refer to caption
(b)
Refer to caption
(c)
Figure 4: (a): A simple assembly producible by a compact zig-zag system, which contains all potential tiles that bind into the assembly via a strength-2 glue in a compact zig-zag system, subject to reflection. (b) A temperature-1 syncTAM system consisting of macrotiles that simulate the behavior of the assembly in (a). Each macrotile is enclosed with a black rectangle, and assembly starts from the bottom-left red tile, concluding with the placement of the white tile (which attaches to the orange macrotile). Note that in the case that the yellow tile were to contain a strength-2 glue to its west, then the red macrotile would be on the east of its macrotile location. (c) A macrotile that may bind to the orange, and yellow macrotiles, shown in (b). The chain of bit readers are colored blue and yellow, with 0-paths denoted by gold arrows, and 1-paths denoted by light blue arrows. Arrows indicate connected glues: the side of a tile at the tail of the arrow attaches to the abutting side of the next tile that the arrow points to, and tiles connect sequentially along the arrow. X’s denote glues to which a tile may be placed, but is blocked as a result of a tile that had previously been placed. For the sake of clarity, each bit which is read is associated with a label to its south, displaying the current bit string which either be encoded in the subsequent bit-reader, or used to resolve to a tile t𝒯.

The above construction correctly simulates 𝒯 provided that a bit-reader never attempts to grow into a macrotile region without a bit-writer chain. Note that provided a strength-1 glue is not exposed at the end of a row in 𝒯, we have constructed 𝒮 to essentially consist of one long path of tiles so that the bit-writer chain must be assembled before a bit-reading gadget reaches the macrotile region. To ensure no strength-1 glue is exposed at the end of a row in 𝒯, we modify 𝒯 so that the north glues of the tiles at the end of each row are marked with a special symbol. We leverage this special glue to further modify 𝒯 so that no strength-1 glue is ever exposed at the end of the row. Consequently, no bit-reader will ever attempt to “read” an empty macrotile region and the construction is correct.

An example compact zig-zag aTAM system, along with a syncTAM system that simulates it is displayed in Figure 5. Along with this, a script which may be used to convert an arbitrary compact zig-zag TAS to a syncTAS may be found at [8].

Corollary 10.

For every standard Turing machine M and input w, there exists a temperature-1 syncTAM system that simulates M on w.

This follows directly from Lemma 7 of [4] (i.e. that, given an arbitrary Turing machine M and input w, there exists a compact zig-zag that simulates M on w), and Theorem 9.

Refer to caption
(a)
Refer to caption
(b)
Figure 5: (a): An example compact zig-zag one-way infinite Turing machine simulation by an aTAM system. (b) The compact zig-zag system shown in (a) simulated by a temperature-1 syncTAS.

5.2 Greater Assembly Connectedness at Temperature 1

Refer to caption
Figure 6: A snapshot of the syncTAS system in Figure 5 modified to exhibit greater assembly connectedness. The assembly is shown alongside it’s reflected counterpart, growing from a single-tile seed, with a tile ’stitching’ both sides together (shown in green).

As previously mentioned, a variety of tile-based, temperature-1 models of self-assembly are computationally universal. Discrimination between types of tiles representing the correct input values from two different locations can be easily done using multiple glue bonds in a temperature-2 system. However, in a temperature-1 system this discrimination must be performed using some sort of hindrance that prevents incorrect tile types from binding while still always allowing the correct tile type to bind. This hindrance (often provided by geometry, as with the “blocker” tile in our construction for Section 5.1) allows for algorithmically correct growth but necessarily results in assemblies with extremely sparse connectivity. The resulting have a binding graph that is a single-tile-wide path that zig-zags and folds back over itself, row by row, throughout the entire assembly. Each row, of increasingly great length, is attached by a single glue bond to the row below it. This results in an assembly that would be extremely “floppy” if physically realized. In all others models where such temperature-1 computation is achieved, this floppiness appears unavoidable since any exposed glue that may be intended to bind to a later-growing portion of the assembly and provide greater stability could, at any time, initiate its own growth (which may or may not agree with correct later growth). However, the precise synchronization enabled by the syncTAM can ensure that two separately-growing portions of an assembly of the same size can be guaranteed to complete at the exact same time, and we can leverage this to provide greater connectivity.

In order to provide this greater connectivity, we must ensure that each row which is placed in the zig-zag TM grows in from both directions, this means that the binding graph will have two paths for each zig / zag of the system, which converge, and diverge in repeating intervals. To accomplish this, the compact zig-zag TM construction displayed in Figure 5 may be reflected along the y axis, and connected to the east, and west of a single-tile seed such that the two identical simulations assemble in parallel.

Due to the synchronous nature of the syncTAM, both identical simulations will arrive at a single-tile wide unoccupied column during the exact same step. All macrotiles which contain a side abutting this unoccupied column only belong to the tiles which grow on the far-east edge of the original aTAM system, and thus only appear abutting the unoccupied column. As such, those tiles may have the side abutting the unoccupied column modified such that they contain a unique “connecting” glue, to which a single “connector” tile type may connect, “stitching” both sides of the assembly together.

As a result of the connector tile stitching both sub-assemblies together along each row, each row will therefore be connected to the row previous by at least two north-south attachments, one on the original TM simulation, and one on its reflection. Consequentially, every row is connected to the previous row by two or more glues. A snapshot of this connected assembly is shown in Figure 6. As with the section previous, a script which takes as input a TM definition, and outputs a more connected syncTAS which simulates it may be found at [8].

References

  • [1] Gagan Aggarwal, Michael H. Goldwasser, Ming-Yang Kao, and Robert T. Schweller. Complexities for generalized models of self-assembly. In Proceedings of ACM-SIAM Symposium on Discrete Algorithms, 2004.
  • [2] Florent Becker, Daniel Hader, and Matthew J. Patitz. Strict Self-Assembly of Discrete Self-Similar Fractals in the abstract Tile Assembly Model, pages 2387–2466. SIAM, 2025. doi:10.1137/1.9781611978322.80.
  • [3] Nathaniel Bryans, Ehsan Chiniforooshan, David Doty, Lila Kari, and Shinnosuke Seki. The power of nondeterminism in self-assembly. Theory of Computing, 9:1–29, 2013. doi:10.4086/toc.2013.v009a001.
  • [4] Matthew Cook, Yunhui Fu, and Robert T. Schweller. Temperature 1 self-assembly: Deterministic assembly in 3D and probabilistic assembly in 2D. In SODA 2011: Proceedings of the 22nd Annual ACM-SIAM Symposium on Discrete Algorithms. SIAM, 2011. doi:10.1137/1.9781611973082.45.
  • [5] Isabel Donoso Leiva, Eric Goles, Martín Ríos-Wilson, and Sylvain Sené. Asymptotic (a)synchronism sensitivity and complexity of elementary cellular automata. In Proceedings of LATIN’24, volume 14579 of LNCS, pages 272–286. Springer, March 2024. doi:10.1007/978-3-031-55601-2_18.
  • [6] David Doty, Jack H. Lutz, Matthew J. Patitz, Robert T. Schweller, Scott M. Summers, and Damien Woods. The tile assembly model is intrinsically universal. In Proceedings of the 53rd Annual IEEE Symposium on Foundations of Computer Science, FOCS 2012, pages 302–310, 2012. doi:10.1109/FOCS.2012.76.
  • [7] Phillip Drake, Daniel Hader, and Matthew J Patitz. Simulation of the abstract tile assembly model using crisscross slats. In 30th International Conference on DNA Computing and Molecular Programming (DNA 30)(2024), pages 3–1. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2024.
  • [8] Phillip Drake and Matthew J. Patitz. Synchronous TAM wiki page and software downloads, 2025. URL: http://self-assembly.net/wiki/index.php/Synchronous_Tile_Assembly_Model_(syncTAM).
  • [9] Nazim Fatès. A guided tour of asynchronous cellular automata. In Jarkko Kari, Martin Kutrib, and Andreas Malcher, editors, Cellular Automata and Discrete Complex Systems, pages 15–30, Berlin, Heidelberg, 2013. Springer Berlin Heidelberg. doi:10.1007/978-3-642-40867-0_2.
  • [10] Nazim Fatès, Damien Regnault, Nicolas Schabanel, and Éric Thierry. Asynchronous behavior of double-quiescent elementary cellular automata. In José R. Correa, Alejandro Hevia, and Marcos Kiwi, editors, LATIN 2006: Theoretical Informatics, pages 455–466, Berlin, Heidelberg, 2006. Springer Berlin Heidelberg.
  • [11] Sándor P. Fekete, Jacob Hendricks, Matthew J. Patitz, Trent A. Rogers, and Robert T. Schweller. Universal computation with arbitrary polyomino tiles in non-cooperative self-assembly. In Proceedings of the Twenty-Sixth Annual ACM-SIAM Symposium on Discrete Algorithms (SODA 2015), San Diego, CA, USA, January 4-6, 2015, pages 148–167, 2015. doi:10.1137/1.9781611973730.12.
  • [12] Oscar Gilbert, Jacob Hendricks, Matthew J. Patitz, and Trent A. Rogers. Computing in continuous space with self-assembling polygonal tiles. In Proceedings of the Twenty-Seventh Annual ACM-SIAM Symposium on Discrete Algorithms (SODA 2016), Arlington, VA, USA January 10-12, 2016, pages 937–956, 2016.
  • [13] Ashwin Gopinath, Evan Miyazono, Andrei Faraon, and Paul WK Rothemund. Engineering and mapping nanocavity emission via precision placement of DNA origami. Nature, 2016.
  • [14] Daniel Hader, Aaron Koch, Matthew J. Patitz, and Michael Sharp. The impacts of dimensionality, diffusion, and directedness on intrinsic universality in the abstract tile assembly model. In Shuchi Chawla, editor, Proceedings of the 2020 ACM-SIAM Symposium on Discrete Algorithms, SODA 2020, January 5-8, 2020, pages 2607–2624, Salt Lake City, UT, USA, 2020. SIAM. doi:10.1137/1.9781611975994.159.
  • [15] Daniel Hader and Matthew J. Patitz. WebTAS (beta): A browser-based simulator for tile-assembly models, 2025. URL: http://self-assembly.net/software/WebTAS_beta/WebTAS_beta-latest/.
  • [16] Jacob Hendricks, Matthew J. Patitz, Trent A. Rogers, and Scott M. Summers. The power of duples (in self-assembly): It’s not so hip to be square. Theoretical Computer Science, 743:148–166, 2018. doi:10.1016/J.TCS.2015.12.008.
  • [17] James I. Lathrop, Jack H. Lutz, and Scott M. Summers. Strict self-assembly of discrete Sierpinski triangles. Theoretical Computer Science, 410:384–405, 2009. doi:10.1016/J.TCS.2008.09.062.
  • [18] Pierre-Étienne Meunier, Matthew J. Patitz, Scott M. Summers, Guillaume Theyssier, Andrew Winslow, and Damien Woods. Intrinsic universality in tile self-assembly requires cooperation. In Proceedings of the ACM-SIAM Symposium on Discrete Algorithms (SODA 2014), (Portland, OR, USA, January 5-7, 2014), pages 752–771, 2014. doi:10.1137/1.9781611973402.56.
  • [19] Pierre-Étienne Meunier and Damien Regnault. Directed Non-Cooperative Tile Assembly Is Decidable. In Matthew R. Lakin and Petr Šulc, editors, 27th International Conference on DNA Computing and Molecular Programming (DNA 27), volume 205 of Leibniz International Proceedings in Informatics (LIPIcs), pages 6:1–6:21, Dagstuhl, Germany, 2021. Schloss Dagstuhl – Leibniz-Zentrum für Informatik. doi:10.4230/LIPIcs.DNA.27.6.
  • [20] Pierre-Étienne Meunier and Damien Woods. The non-cooperative tile assembly model is not intrinsically universal or capable of bounded turing machine simulation. In Proceedings of the 49th Annual ACM SIGACT Symposium on Theory of Computing, STOC 2017, pages 328–341, New York, NY, USA, 2017. ACM. doi:10.1145/3055399.3055446.
  • [21] Matthew J. Patitz, Robert T. Schweller, and Scott M. Summers. Exact shapes and Turing universality at temperature 1 with a single negative glue. In Luca Cardelli and William M. Shih, editors, DNA Computing and Molecular Programming - 17th International Conference, DNA 17, Pasadena, CA, USA, September 19-23, 2011. Proceedings, volume 6937 of Lecture Notes in Computer Science, pages 175–189. Springer, 2011. doi:10.1007/978-3-642-23638-9_15.
  • [22] Loïc Paulevé and Sylvain Sené. Boolean networks and their dynamics: the impact of updates. (Chapter) Systems biology modelling and analysis: formal bioinformatics methods and tools, pages 173–250, November 2022. doi:10.1002/9781119716600.ch6.
  • [23] Damien Regnault, Nicolas Schabanel, and Éric Thierry. Progresses in the analysis of stochastic 2d cellular automata: A study of asynchronous 2d minority. Theoretical Computer Science, 410(47):4844–4855, 2009. doi:10.1016/j.tcs.2009.06.024.
  • [24] Paul W. K. Rothemund and Erik Winfree. The program-size complexity of self-assembled squares (extended abstract). In STOC ’00: Proceedings of the thirty-second annual ACM Symposium on Theory of Computing, pages 459–468, Portland, Oregon, United States, 2000. ACM. doi:10.1145/335305.335358.
  • [25] Grigory Tikhomirov, Philip Petersen, and Lulu Qian. Fractal assembly of micrometre-scale DNA origami arrays with arbitrary patterns. Nature, 552(7683):67–71, 2017.
  • [26] Ning Wang, Chang Yu, Tingting Xu, Dan Yao, Lingye Zhu, Zhifa Shen, and Xiaoying Huang. Self-assembly of dna nanostructure containing cell-specific aptamer as a precise drug delivery system for cancer therapy in non-small cell lung cancer. Journal of Nanobiotechnology, 20(1):1–19, 2022.
  • [27] Erik Winfree. Algorithmic Self-Assembly of DNA. PhD thesis, California Institute of Technology, June 1998.
  • [28] Damien Woods, David Doty, Cameron Myhrvold, Joy Hui, Felix Zhou, Peng Yin, and Erik Winfree. Diverse and robust molecular algorithms using reprogrammable dna self-assembly. Nature, 567(7748):366–372, 2019. doi:10.1038/S41586-019-1014-9.

Appendix A Technical Appendix

This Technical Appendix contains additional technical definitions and details related to the proofs in the main body of the paper.

A.1 Simulation Definition Details

In this section, we present the formal definitions of simulation between aTAM and syncTAM systems.

The definition of the syncTAM leads itself to a simulation definition that is equivalent to the definition of intrinsic simulation for a TAS simulating another TAS. In large part, our definitions come from [7].

From this point on, let 𝒯=(T𝒯,σ𝒯,τ𝒯) be a TAS, and let 𝒮=(T𝒮,σ𝒮,τ𝒮) be some syncTAS that simulates 𝒯. For each such simulation we fix a positive integer m called the scale factor with the intention that m×m blocks of locations from 𝒮 map to individual locations in 𝒯. In other words, simulation happens at scale so that by “zooming out” by a factor of m, the assemblies in 𝒮 resemble corresponding assemblies in 𝒯.

In the context of 𝒮, we call each m×m block of locations in 2 a macrotile location under the convention that the origin (0,0)2 occupies the south-westernmost location in one of the m×m blocks. That is, macrotiles locations are squares with southwest corners of the form (cx,cy) and northeast corners of the form (m(x+1)1,m(y+1)1) where x,y. Given an assembly α𝒜[𝒮], we use the notation x,yc[α] to refer to the set of tiles in α which occupy locations in the macrotile location whose southwest corner is (cx,cy). This set of tiles is referred to as the macrotile located at (x,y). Note that in this context, we keep track of the tiles that occupy a macrotile with coordinates relative to (x,y) so that two macrotiles which differ only by a translation are identical. We denote the set of all possible macrotiles of scale factor m made from tiles in T𝒮 as m[S].

A partial function R:m[T𝒮]T𝒯 is called a macrotile representation function from T𝒮 to T𝒯 if, for any macrotiles α,βm[T𝒮] where αβ and αdomR, then R(α)=R(β). In other words, a macrotile representation function may not map a macrotile to an individual tile type in T𝒯, but if it does, then any additional tile attachments do not change how the macrotile is mapped under R. In the context of some 𝒮-assembly sequence, a macrotile is said to resolve to a tile type tT𝒯 when an assembly maps under R to t, but the prior assembly is not in the domain of R.

From a macrotile representation function R, a function R:𝒜T𝒮𝒜T𝒯, called the assembly representation function, is induced which maps entire assemblies in 𝒮 to assemblies in 𝒯. This function is defined by applying the function R to each macrotile location containing slats. We also use the notation R1(α) to refer to the producible pre-image of an assembly α𝒜T. That is, R1(α)={α𝒜[𝒮]|R(α)=α} so that R1(α) includes every 𝒮-producible assembly mapping to α under R. We say that an assembly α𝒜S maps cleanly to an assembly α𝒜T under R if R(α)=α and no slats in α occupy any macrotile whose location is not adjacent (not including diagonally) to a resolved macrotile. Formally, if α maps cleanly to α, then for each non-empty macrotile x,ym[α], there exists some vector (u,v)2 such that (u,v)1 so that x+u,y+vm[α]domR.111Note that this definition requires that the seed assembly σ𝒮 resolves to at least one tile in 𝒯. This will suffice for the constructions of this paper, but if needed this definition can be relaxed to specifically allow for the seed assembly to grow before resolving. The definition of maps cleanly requires that any non-empty but unresolved macrotile has a resolved macrotile to its N, E, S. or W, thus ensuring that the growth of tiles is only occurring in macrotile locations that map to locations in 𝒯 adjacent to tiles, and thus with some potential to receive a tile.

Definition 11 (Equivalent Productions).

We say that 𝒮 and 𝒯 have equivalent productions (under R) and we write 𝒮𝒯 if the following conditions hold:

  1. 1.

    {R(α)|α𝒜[𝒮]}=𝒜[𝒯].

  2. 2.

    {R(α)|α𝒜[𝒮]}=𝒜[𝒯].

  3. 3.

    for all α𝒜[𝒮],α maps cleanly to R(α).

Definition 12 (Follows).

We say that 𝒯 follows 𝒮 (under R), and we write 𝒯𝒮β, for some α,β𝒜[𝒮], implies that R(α)𝒯R(β).

Definition 13 (Models).

We say that 𝒮 models 𝒯 (under R), written 𝒮R𝒯 if for every α𝒜[𝒯], there exists a non-empty subset ΠαR1(α), such that for all β𝒜[𝒯] where α𝒯β, the following conditions are satisfied:

  1. 1.

    for every αΠα, there exists βR1(β) such that α𝒮β

  2. 2.

    for every α′′R1(α) and β′′R1(β) where α′′𝒮β′′, there exists αΠα such that α𝒮α′′.

In Definition A.1 above, the set Πα is defined to be a set of assemblies representing α from which it is still possible to produce assemblies representing all possible β producible from α. Informally, the first condition specifies that all assemblies in Πα can produce some assembly representing any β producible from α, while the second condition specifies that any assembly α′′ representing α that may produce an assembly representing β is producible from an assembly in Πα. In this way, Πα represents a set of the earliest possible representations of α where no commitment has yet been made regarding the next simulated assembly. Requiring the existence of such a set Πα for every producible α ensures that non-determinism is faithfully simulated. That is, the simulation cannot simply “decide in advance” which tile attachments will occur.

Definition 14 (Simulates).

Given an aTAM system 𝒯=(T,σ,2), a syncTAM system 𝒮=(S,σ𝒮,c), and a macrotile representation function R from 𝒮 to 𝒯, we say that 𝒮 intrinsically simulates 𝒯 (under R) if 𝒮R𝒯 (they have equivalent productions), 𝒯R𝒮 and 𝒮R𝒯 (they have equivalent dynamics).

A.2 Window Movie Lemma

The Window Movie Lemma is a powerful tool often used in impossibility results for tile-based self-assembly models. Originally introduced in [18], it allow us to perform surgery on an assembly sequence to obtain a new producible assembly (assuming that some conditions are met). In its simplest form, the lemma states that if there are two enclosed areas of the plane that have the same “shape” and the same glues appear in the same order at the same position relative to these two enclosed shapes (that is, the same window movie), then the subassemblies produced in the two enclosed areas are interchangeable. See Figure 7 for a graphical representation of this simple case of the Window Movie Lemma.

Figure 7: An illustration of the window movie lemma. The squares represent tiles and the color of the tile is an indication of the glue type the tile exposes along the window w (represented by a red dotted box). The numbers on the tile represent the relative order in which the glues contained on the tile appear on the window as the assembly sequence is played forward. The Window Movie Lemma tells us that since αLαR is producible (shown top left) and βLβR is producible (show bottom left) and because the assembly sequence produces the same window movie sequence along w, the assemblies shown on the right are also producible.

The above surgical procedure can be extended in order to pump assemblies. Indeed, if there is a fixed width subassembly of sufficient length, then there is some window movie which repeats. We can then pump the subassembly between these two window movies indefinitely. See Figure 8 for a graphical representation showing how assemblies can be pumped via the Window Movie Lemma.

The following definitions are taken from [18], and replicated here for completeness. A window is an edge cut which partitions the lattice graph (2 in 2D or 3 in 3D) into two regions. Given some window w and some assembly sequence α in a TAS 𝒯, a window movie M is defined to be the ordered sequence of glues presented along w by tiles in 𝒯 during the assembly sequence α. Informally, the window w is a cut dividing two regions of tile locations and M is constructed by recording the glues which appear on the cut and their relative order as the assembly sequence α is played forward. More formally, a window movie is the sequence Mwα={(vi,gi)} of pairs of grid graph vertices vi and glues gi, given by order of appearance of the glues along window w during α. Furthermore, if k glues appear along w during the same assembly step in α, then these glues appear contiguously and are listed in lexicographical order of the unit vectors describing their orientation in Mwα.

Window Movie Lemma

Let α={αi} and β={βi} be assembly sequences in TAS 𝒯 and let α,β be the result assemblies of each respectively. Let w be a window that partitions α into two configurations αL and αR and let w=w+c be a translation of w that partitions β into two configurations βL and βR (with αL and βL being the configurations containing their respective seed tiles). Furthermore define Mwα and Mwβ to be the window movies for α,w and β,w respectively. Then if Mwα=Mwβ, the assemblies αLβR and βLαR (where βL=βLc and βR=βRc) are also producible.

Figure 8: An illustration of how the Window Movie Lemma can be leveraged to pump assemblies. The squares represent tiles and the color of the tile is an indication of the glue type the tile exposes along the window w (represented by a red dotted box). The numbers on the tile represent the relative order in which the glues contained on the tile appear on the window as the assembly sequence is played forward. The top assembly shows a fixed width ribbon grown sufficiently long. By the pigeon hole principle, there must exist a window that has the same window movie when translated to two distinct positions along the ribbon (as indicated by the same colored tiles with the same numbers on them occurring twice in the ribbon). Then, by the Window Movie Lemma, the assembly between these two translations of the window can be repeated indefinitely. This is illustrated in the bottom portion of the figure which shows that the portion of the ribbon between repeating window movies can be also be assembled beginning where the last repeating window movie is located.