Rod-Like Virus-Based Multiarm Colloidal Molecules

Abstract: We report on the construction of multiarm colloidal molecules by tip-linking filamentous bacteriophages, functionalized either by biological engineering or chemical conjugation. The affinity for streptavidin of a genetically modified vector phage displaying Strep-tags fused to one end of the viral particle is measured by determining the dissociation constant, Kd. In order to improve both the colloidal stability and the efficiency of the self-assembly process, a biotinylation protocol having a chemical yield higher than 90% is presented to regioselectively functionalize the cystein residues located at one end of the bacteriophages. For both viral systems, a theoretical comparison is performed by developing a quantitative model of the self-assembly and interaction of the modified viruses with streptavidin compounds, which accurately accounts for our experimental results. Multiarm colloidal structures of different valencies are then produced by conjugation of these tip-functionalized viruses with streptavidin activated nanoparticles. We succeed to form stable virus based colloidal molecules, whose number of arms, called valency, is solely controlled by tuning the molar excess. Thanks to a fluorescent labeling of the viral arms, the dynamics of such systems is also presented in real time by fluorescence microscopy.

\section{Introduction}

The design of structured materials at the mesoscopic scale is a challenge that can be addressed with two main strategies usually referred as top-down and bottom-up approaches.\cite{Rev1,Rev2} While the top-down methods consist in the elaboration of patterns and in their progressive size reduction, the bottom-up approach requires the use of nano- or micro-scaled building blocks for which the self-assembly can be driven by entropic or enthalpic contributions allowing for producing hierarchical architectures.\cite{Frenkel} Specifically, most of reported studies have focused on sphere-based assemblies with specific sizes and/or functionalizations opening the field of so-called ``colloidal molecules''.\cite{vanB,Man,Man2} To go beyond the spherical symmetry and developing more complex architectures, some investigations have been performed using nanorods as elementary building blocks.\cite{Salant,Manna,Huang,Fraden,Solomon}
Mainly two distinct paths have been considered for creating multiarm colloidal molecules: the core-first and the arm-first approaches. In the former case, the core is a multivalent initiator from which the arms can be grown. In this way, multivalent-branched nanocrystals \cite{Manna2} or polymer stars can be obtained by a specific crystal growth or by the polymerization of monomers around the core.\cite{Pol1,Pol2,ER3,ER4,ER5,ER6,ER7} This method allows for a control over the length and/or the width of each constituent but limits the construction of objects to a specific valency and a tailored size.
The second approach consists in the self-assembly of linear constituents into multivalent star- or flower-like structures. This technique requires an end-to-end assembly of the particles which can be achieved in different ways such as DNA templating,\cite{Dujardin,Pal} amphiphilic interactions,\cite{Nie,ER1,ER2} biorecognition \cite{Salant,Huang,Sweeney,Caswell,Linder} or tip-linking reaction of linear polymer chains.\cite{

Pol3,Pol4}
In this paper, we focus on the end-to-end self-assembly of filamentous viruses and aim to obtain multiarm structures of tunable valency (\textit{i.e.} variable number of arms). Such viral rod-like particles can be seen both as elementary building blocks monodisperse in size and shape, and as versatile scaffolds for biological engineering and chemical functionalization of their coat proteins.\cite{Molek,Ng,Cao} Prior studies have reported the use of bacteriophages, mainly modified by phage display, that self-assemble into chains of various length,\cite{Sweeney,Hess} rings \cite{Nam} or radial structures.\cite{Huang,Fraden} These works are based on the genetic engineering of viruses with mutated tip proteins (minor proteins p3 and p9, corresponding respectively to the two extremities of the phage) onto which specific polypeptidic sequences are fused. The formation of chains or rings is then driven either by the complementarity of the involved polypeptidic sequences \cite{Sweeney} or by the introduction of an external crosslinker.\cite{Nam} Star-like structures were also achieved by using the affinity of biotin for streptavidin \cite{Fraden} or
by biorecognition of a specific motif Histidine-Proline-Glutamine (HPQ),\cite{Devlin} called Strep-tag, with streptavidin coated nanoparticles.\cite{Huang} Specifically, a mutant called M13-AntiStreptavidin (M13AS), initially isolated through the screening of a phage display library, has been already used in templating hierarchically organized hybrid materials.\cite{SWLee,Huang,SWLee2016} Here, we study this M13AS mutant displaying Strep-tags at one end by first quantifying its affinity for streptavidin with the determination of its dissociation constant in solution, $K_{d}$.
Because of the relatively high dissociation constant $K_{d}$ found, limiting therefore the self-assembly efficiency for the formation of colloidally stable multiarm structures, we have developed a chemical alternative approach based on the biotinylation of cystein residues present at one tip (p3 proteins) of the so-called M13C7C bacteriophages. The protocol relies on the thiol based bioconjugation, and it has been optimized to be both highly specific (no other coat protein modified) and efficient (yield of about 92\%). We then discuss a quantitative model accounting for the interaction  with streptavidin of both the biologically and the chemically functionalized phages for various initial conditions.
When exposed to a dispersion of streptavidin coated nanoparticles, biotinylated phages self-assemble into colloidally stable multiarm structures, as the affinity of streptavidin for biotin is several orders of magnitude higher than for Strep-tags. Moreover, these multiarm molecules display a tunable valency which can be continuously monitored by varying the relative molar excess of biotinylated phages and streptavidin activated nanoparticles. This approach makes also possible a second chemical functionalization of the phages by their body labeling with fluorescent dyes, allowing for the \textit{in situ} observation by optical microscopy of their self-assembly into multiarm colloidal molecules.

\section*{Results and discussion}

\subsection*{M13AS Binding Affinity and Biotinylation Yield of M13C7C-B}

In our mixtures of M13 viruses (either AntiStreptavidin phages, M13AS or biotinylated phages, M13C7C-B) and streptavidin coated nanoparticles (NP) exhibiting $q$ streptavidin molecules (Strep) per NP, the principle of mass conservation applied on virus and streptavidin species provides the two following relations:

\begin{equation}
[M13]_{0} =
[M13]_{free} + [M13]_{bound}
\label{eq:M13ASmass}
\end{equation}

with [M13]$_{0}$ the initial virus concentration, [M13]$_{free}$ and [M13]$_{bound}$ the respective concentrations at equilibrium of free and bound (or reacted) viruses with streptavidin coated nanoparticles;

\begin{equation}
q\times [NP]_{0} =
\left
(q \times [NP]_{free} + (q-n)\times [NP]_{bound}
\right) + n\times [NP]_{bound}
\label{eq:NPmass}
\end{equation}

where $n$ is the number of viruses bound per NP ($1\leq n \leq q$), [NP]$_{0}$ and [NP]$_{free}$ the respective initial and unreacted concentrations of nanoparticles, and [NP]$_{bound}$ the concentration of nanoparticles conjugated with at least one M13 virus (Fig. \ref{Fig1}). The binding affinity of a ligand for a receptor is given by the dissociation constant, $K_{d}$, which results from the following equilibrium in solution:\cite{Sanders}

\begin{equation}
M13_{free} + Strep_{free}  \rightleftharpoons M13_{bound}
\nonumber
\end{equation}

and is defined as:

\begin{equation}
K_{d} \equiv
\frac{[M13]_{free} \times [Strep]_{free}}{[M13]_{bound}}
\label{eq:Kd}
\end{equation}

 where [Strep]$_{free}=q \times [NP]_{free} + (q-n)\times [NP]_{bound}$. Qualitatively, the lower the $K_{d}$ value is, the more stable the complex is and thus the stronger the affinity.
 In the case of the interaction between biotin and streptavidin protein,
 $K_{d}$ reaches a value of 10$^{-15}$~M, considered the strongest non-covalent bond found in nature.\cite{Thermo}
 The Histidine-Proline-Glutamine (HPQ) motif is a well-known Strep-tag with however a lower affinity for streptavidin than biotin, typically found in the $\mu$M range.\cite{Devlin} However, the affinity of M13AS phage where the HPQ motif is displayed on each of its five p3 proteins, has not yet been determined.
 Using Eqs. \ref{eq:M13ASmass} and \ref{eq:NPmass}, let's rewrite $K_{d}$ as a function of an unique variable, the fraction of bound nanoparticles $f(NP)\equiv [NP]_{bound}/[NP]_{0}$:

 \begin{equation}
 K_{d} =
 \left(
 \frac{[M13]_{0}}{n \times f(NP)} - [NP]_0
 \right)
 \times
 \left(
 q-n \times f(NP)
 \right)
 \label{eq:KdNP}
 \end{equation}

 Thanks to Eq. \ref{eq:KdNP}, the M13AS-streptavidin dissociation constant, $K_{d}$, can be determined by counting by TEM both $f(NP)$ and $n$, the average number of viruses per nanoparticle having at least one virus bound (defined also as the valency of the self-assembled structure). Note that using this model, a precise estimation of $K_{d}$ can only be performed if no steric hindrance limits the binding of viruses to nanoparticles. This means practically that $n \ll q$,
 as experimentally confirmed in Fig. \ref{Fig-counts}(a) and (b) for which $n\simeq~1.5$ (see Table \ref{tab:NP}).

\begin{figure}
  \includegraphics[width=0.7\columnwidth]{Figure1new.eps}
  \caption{\label{Fig1} Transmission electron microscopy (TEM) images of M13AS ((a) and (c)) and M13C7C-B ((b) and (d)) viruses complexed with streptavidin-coated nanoparticles (NP). (a),(b) Single viruses and (c),(d) pairs of viruses bound to a nanoparticle. The scale bars represent 100~nm. (e) Schematic representation of two viruses ($n=2$) bound to a NP exhibiting $q=$8 streptavidin proteins at its surface.}
\end{figure}

\begin{figure}
	\includegraphics[width=1\columnwidth]{Fig-counts.eps}
	\caption{\label{Fig-counts} Distribution of the number of bound viruses per
	nanoparticle in suspensions of (a) M13AS at initial concentration of 10 mg/mL (5.4  $\times 10^{-7}$~M) with a molar excess of virus per nanoparticle $\alpha=$~41 (Sample \textbf{1}), (b) [M13AS]$_0$=1 mg/mL (5.4 $\times 10^{-8}$~M) with $\alpha=$~10 (Sample \textbf{2}), (c) biotinylated viruses at [M13C7C-B]$_0$=0.1 mg/mL (5.4 $\times 10^{-9}$~M) with $\alpha$=1.28 (Sample \textbf{3}). $N$ indicates the total number of NP counted by TEM for establishing the statistics of each sample. }
\end{figure}

As M13AS viruses are in large excess compared to NP (see Table \ref{tab:NP}), only the concentrations at equilibrium of free and bound NP as well as the average number of viruses grafted per NP, $n$, can been counted by TEM on a given sample ($[M13AS]_{free}$ being too high to be determined at the same time).
In order to determine $K_{d}$,

two specific mixtures of M13AS and streptavidin-coated nanoparticles (Samples \textbf{1} and \textbf{2}, described in Table \ref{tab:NP}) were realized with different initial conditions, for which the molar excess is introduced, $\alpha$:

\begin{equation}
\alpha \equiv
\frac{[M13]_{0}}{[NP]_{0}}
\label{eq:alpha}
\end{equation}

\begin{table}
\caption{\label{tab:NP} Summary of the experimental conditions and measurements obtained by TEM for the three mixtures of viruses and streptavidin coated nanoparticles. Samples \textbf{1} and \textbf{2} correspond to biologically engineered M13AS phages and sample \textbf{3} to tip-biotinylated M13C7C-B viruses. The fraction of bound nanoparticles, $f(NP)$, and the average number of viruses per bead, $n$, result from TEM counting shown in the distributions plotted in Fig. \ref{Fig-counts}. The dissociation constant, $K_{d}$, is calculated for samples \textbf{1} and \textbf{2} according to Eq. \ref{eq:KdNP}, and is taken from literature \cite{Thermo} for sample \textbf{3}.}
{\setlength{\extrarowheight}{10pt}
\begin{tabular}
	{c c c c}

  \multicolumn{1}{c}
  {Sample} & \textbf{1} & \textbf{2} & \textbf{3} \\
  \hline \hline
  Virus & M13AS & M13AS & M13C7C-B \\

  $[M13]_{0}$~(M)& 5.4~$\times$~10$^{-7}$ & 5.4~$\times$~10$^{-8}$ & 5.4~$\times$~10$^{-9}$ \\

  $[NP]_{0}$~(M)& 1.3~$\times$~10$^{-8}$ & 5.4~$\times$~10$^{-9}$ &  4.2~$\times$~10$^{-9}$ \\

  $\alpha=[M13]_{0}/[NP]_{0}$& 41 & 10 & 1.28\\

  $f(NP)$ & 0.27 & 0.11 & 0.61\\

  $n$ & 1.62 & 1.31 & 1.94 \\

  $K_{d}$~(M)& 9~$\times$~10$^{-6}$ & 3~$\times$~10$^{-6}$ & 10$^{-15}$ \\
 \hline
\end{tabular}}
\end{table}

Both mixtures were observed in TEM (Fig. \ref{Fig1}) and the distributions of the number of viruses grafted per bead is provided in Fig. \ref{Fig-counts}(a) and (b) for samples \textbf{1} and \textbf{2}, respectively. The mean values of $n$ and $f(NP)$ are deduced from these distributions, and consequently
$K_{d}$ is obtained thanks to Eq. \ref{eq:KdNP}. A similar dissociation constant (within the experimental margin of error) $K_{d}$=6~($\pm$~3)~$\times$~10$^{-6}$~M is found for both samples (Table \ref{tab:NP}). This average value is also in good agreement with the one reported in literature for the Histidine-Proline-Glutamine Strep-tag of about 1~$\mu$M.\cite{Devlin} The slightly lower binding affinity found in our case can be explained by the following reason: The determination of binding affinity is usually performed with the target compound (here: streptavidin) deposited on a solid surface, whereas our experiments were performed in bulk with dispersions of streptavidin coated NPs. These NPs not only diffuse slower than molecular streptavidin, but also induce steric hindrance as soon as some viruses are bound making more difficult further binding of other viruses.
Nevertheless, this relative high value of the dissociation constant, for which $K_d\gg[M13AS]_0\geq[NP]_0$ (Table \ref{tab:NP}), does not allow for a tight and non-reversible binding. A first alternative would be to design
a Strep-tag with an affinity in the nM range, providing then a suitable efficiency for high enough initial virus ($>$~1mg/mL) and bead concentrations. The second alternative is to work with biotin derivatives ($K_d\simeq 10^{-15}$~M), which have been chemically grafted at one of the virus tips (see \textit{Materials and methods}). In order to determine the biotinylation yield of the M13C7C-B phages, sample \textbf{3} has been prepared (Table \ref{tab:NP}), where phages and streptavidin activated nanoparticles have been introduced in similar amount ($\alpha\simeq 1$,) allowing for their simultaneous counting by TEM (Fig. \ref{Fig1}). Specifically, among a total number of 229 counted viruses, 210 were bound by their tip to nanoparticles and 19 were found free, giving a fraction of bound M13C7C-B of $f(M13)=0.92$ and therefore a yield of biotinylation of $92\%$. It is worth mentioning that the functionalization of the virus tips by biotin is not only efficient, but also highly selective as shown by the absence of any virus based by-products. Moreover, the specificity of the reaction between beads and functionalized phages has been checked by a control experiment with raw M13 viruses, and, as expected no interaction is observed with nanoparticles (see Figure S3 in Supporting Information).

\subsection*{Quantitative Model for the Formation of Virus/NP Structures}

\begin{figure}[!h]
	\includegraphics[width=0.5\columnwidth]{Model.eps}\\
	\caption{\label{Fig:Model} Fraction of bound nanoparticles, $f(NP)$, as modelized by Eq. \ref{eq:model}, as a function of the molar excess, $\alpha$, obtained for different average numbers of viruses per NP, $n$. Results calculated for: (a) M13AS virus strain with K$_{d}\simeq6\times 10^{-6}$~M for two different initial virus concentrations, [M13AS]$_{0}$;
		(b) biotinylated viruses, M13C7C-B, with K$_{d}= 10^{-15}$~M. In this case, the results are independent of the initial virus concentration, as [M13AS]$_{0}\gg $K$_{d}$, and only depend on the molar excess $\alpha$. The cyan, blue and green open square symbols represent the samples \textbf{1}, \textbf{2} and \textbf{3}, respectively, as described in Table \ref{tab:NP}. The green, orange, blue and pink full square symbols represent the four experimental samples of respective average valency $n=2.1, 4.2, 6.1$ and 8.3, as reported in the next Section. }
\end{figure}
Rewriting the expression of $K_{d}$ given in Eq. \ref{eq:KdNP}, the fraction of bound nanoparticles $f(NP)$ can be accounted by the following quadratic equation:

\begin{equation}
{f(NP)}^2 - \frac{1}{n}
\left(
\frac{\alpha K_d }{[M13]_0}  + q + \alpha
\right)f(NP) + \frac{q \alpha}{n^2}
=0
\label{eq:model}
\end{equation}

Eq. \ref{eq:model} can be easily solved analytically to provide $f(NP)$ for different initial conditions.
The results are plotted in Fig. \ref{Fig:Model}. When the binding affinity for streptavidin is low, as for M13AS viruses, the self-assembly yield to form multiarm structures ($n\geq2$) remains weak ($f(NP)<25\%$), whatever the initial virus and NP concentrations (Fig. \ref{Fig:Model}(a)). At a fixed value of the dissociation constant $K_{d}$, only the initial conditions affect the equilibrium, as illustrated by the two experimental data points corresponding to samples \textbf{1} and \textbf{2}, which are in very good agreement with the predictions of the model.
Fig. \ref{Fig:Model}(b) shows that strongly increasing the binding affinity favors the formation of multiarm structures, for which a total reaction can be obtained. More importantly, the valency of multiarm self-assemblies can be controlled and tuned \textit{only} by varying the molar excess $\alpha$. Finally, an outstanding agreement between our model and sample \textbf{3} is found  (Fig. \ref{Fig:Model}(b)), especially considering the absence of any free parameter. In conclusion, our model quantitatively accounts for the self-assembly behavior of our two experimental systems, and it confirms that biotinylated phages are the most suitable one to form multiarm colloidal molecules, as described in the next paragraph.

\subsection*{Construction of Multiarm Colloidal Molecules}

\begin{figure}
	\includegraphics[width=1\columnwidth]{FigureTEM.eps}\\
	\caption{\label{Fig:M13C7CBS} Virus based multiarm self-assemblies observed by TEM obtained from the reaction between tip-biotinylated phages (M13C7C-B) and streptavidin coated nanoparticles.
		The number of arms can be continuously tuned, as illustrated with structures of valency $n=$2 ((a), dimer), 4 ((b), tetramer), 6 ((c), hexamer) and 8 ((d), octamer). (e) Large field of view with three multiarm molecules of average valency $n=$8. Inset: zoom-in of the multiarm structure core. Except in the inset, the scale bars represent 500~nm.}
\end{figure}

\begin{figure}
	\includegraphics[width=0.9\columnwidth]{M13C7C-B-NP.eps}\\
	\caption{\label{Fig:M13C7CBSStat} (a) Distribution of the number of viruses per nanoparticle for four different samples, at fixed biotinylated virus concentration, [M13C7C-B]$_0$, and increasing nanoparticle concentration, [NP]$_0$. The arithmetic mean, or average valency $n$, is indicated and the solid lines correspond to Gaussian fits. $N$ is the total number of NP counted by TEM for establishing the statistics of each sample. (b) The average valency $n$ increases linearly with the molar excess, and saturates when $n$ reaches the maximum number of available streptavidin proteins per NP, $q \simeq 8$. The error bars are the standard deviation obtained from the Gaussian fits.
	}
\end{figure}

\begin{figure}
	\includegraphics[width=0.6\columnwidth]{FigureFluo.eps}\\
	\caption{\label{Fig:M13C7CFluo} Optical fluorescence microscopy images showing multiarm colloidal molecules formed by fluorescently labeled M13C7C-B phages. Different valencies, $n$, can be observed \textit{in situ} in aqueous dispersions.
		The scale bar, identical for the six images, represents 1~$\mu m$.}
\end{figure}

Multiarm colloidal molecules were achieved by preparing four samples of M13C7C-B phages mixed with aqueous dispersions of streptavidin coated nanoparticles (see \textit{Materials and methods} section). As shown in Fig. \ref{Fig:M13C7CBS} and in Supplementary Fig. 4, the rod-based colloidal molecules are radial structures possessing a single central NP surrounded by up to $n=8$ arms formed by tip-linked M13C7C-B functionalized phages.
Figure \ref{Fig:M13C7CBSStat} displays the increase of the average valency $n$ of multiarm assemblies when the initial concentration of NP
is decreased at a fixed virus concentration.
This behavior is \textit{quantitatively} consistent with our theoretical model, as plotted in Fig. \ref{Fig:Model}(b) where the full symbols represent the four samples of colloidal molecules reported in Fig. \ref{Fig:M13C7CBSStat}. We experimentally confirm that the molar excess is the key-parameter to tune the valency of the self-assemblies. After a linear increase, the average valency $n$ saturates in Fig. \ref{Fig:M13C7CBSStat}(b) when it reaches the maximum number of available streptavidin binding sites per bead, $q$. At this saturated value, $n$ becomes independent of the initial molar excess $\alpha$ (for $\alpha \geq q$), as predicted by the model (Fig. \ref{Fig:Model}(b)).
In the results shown in Fig. \ref{Fig:M13C7CBSStat}(a), all samples present a similar distribution of virus populations (except for the lowest average valency $n=2.1$ which is asymmetric), whose width, defined as the standard deviation of the corresponding Gaussian fits, is $\delta n \pm 2$. This width of the average valency is narrow enough
to obtain slightly polydisperse (in terms of number of arms) batches of star-like colloidal molecules for the four different prepared samples (corresponding to the four different initial NP concentrations).
Thanks to an efficient and regio-selective biochemical functionalization of the M13C7C phages followed by a conjugation with streptavidin activated beads, the large scale production of multiarm structures of well-defined and tunable valency has been achieved. Furthermore, biotinylated phages are an interesting system to work with because of their ability to be twice chemically modified without losing their functionality. We have succeeded in labeling with fluorescent dyes the virus body, targeting mainly the p8 major coat proteins. Labeled M13C7C-B phages can therefore be visualized at the single particle level by optical fluorescence microscopy (Fig. \ref{Fig:M13C7CFluo}), allowing for \textit{in situ}, \textit{i.e.} at equilibrium, sample observation as shown in the movie presented in Supporting Information.
If the structures are identical in shape to those seen by TEM (Fig.
\ref{Fig:M13C7CBS}), information on the dynamics of the different multiarm structures are here available. Remarkably, some flexibility can be noticed between the arms of the self-assemblies, contrary for instance to multipod structures obtained by crystal growth.\cite{Solomon} This flexibility in the link between bound viruses is distinguishable from the Brownian motion of the multiarm structures themselves. Such \textit{in situ} observations
provide the experimental equilibrium state of virus based multiarm structures.

\section*{Conclusions}

In this paper, we have studied the self-assembly of functionalized filamentous viruses into multiarm colloidal structures. First, we used an engineered phage (M13AS), modified at one end to display a polypeptidic sequence (HPQ motif) exhibiting affinity for streptavidin. We determined the dissociation constant of the system in solution to be  $K_{d} = 6 \times$ 10$^{-6}$~M consistent with literature value for this protein tag.
In order to strengthen the affinity for streptavidin, a selective and efficient biotinylation protocol of one of the phage tips has been successfully performed. When conjugated with streptavidin activated nanoparticles, these chemically modified phages self-assemble into colloidally stable multiarm structures, whose valency has been continuously tuned from 1 to 8 by varying the molar excess. An additional chemical functionalization of the phage body by fluorescent dyes has been done to optically visualize the multiarm colloidal molecules at the single particle scale and to get \textit{in situ} information on their dynamics.
The results presented here show how regio-functionalized viruses can be used as building blocks for mesostructured materials and for the development of a next generation of rod based colloidal particles.

\section*{Materials and methods}

M13-AntiStreptavidin phage presents Strep-tags on each of its five p3 proteins localized at one tip. M13AS was kindly provided by S.-W. Lee from the University of Berkeley (USA). M13AS was isolated through the screening of a phage display library and several rounds of selection for streptavidin, with the following inserted sequence on the p3 proteins, identified by DNA sequencing (see Figure S1 in Supporting Information for the complete sequencing): AC\textbf{\underline{HPQ}}GPPCGGGS.\cite{Yoo} The specific motif Histidine-Proline-Glutamine (HPQ),\cite{Devlin} called Strep-tag, is known for its affinity with streptavidin. Large amplification of these phages of molecular weight $M_w=1.85\times10^7$~g/mol was performed using standard biological procedures.\cite{Sambrook} Before further use, M13AS phages were extensively dialyzed against PBS buffer (pH 7.5, I=20~mM).

A second filamentous phage, called M13C7C and displaying cystein residues on the exposed part of each p3 protein (cystein amino acids otherwise absent from the other coat proteins of the phage), has been used in this study and bioconjugated with maleimide activated biotin. Infecting solution was purchased from New England Biolabs (NEB, USA), and corresponds to a mixture of several M13C7C clones. The amino acid sequences inserted on each of the five p3 proteins of the viruses display 7 random amino acids flanked by two oxidized cysteins which form a disulfide bridge. After phage titration, we isolated a plaque consisting of a colony of bacteria infected by a single virus clone and thus a unique amino acid sequence (see Figure S2 in Supporting Information for the complete DNA sequencing). M13C7C phages were also grown and purified using standard biological protocols.\cite{Sambrook}

In order to have free thiols available for conjugation with maleimide activated compounds, a reduction step of the disulfide bridges is required. This was done by introducing a reducing agent, as Tris(2-carboxyethyl)phosphine hydrochloride (TCEP, Thermo Scientific) used at a concentration around 2~mM.\cite{Heinis, Heinis2} To prevent any re-oxidization of the free thiols and the reformation of the disulfide bridges, all the reaction process occurred in degassed media and in the presence of 1 mM of ethylenediaminetetraacetic acid (EDTA, Sigma Aldrich).\cite{Heinis} The viruses were first depleted by the addition of 20~\% (v/v) of a PEG-8000/NaCl mixture (200 g/L of PEG-8000, 2.5 M of NaCl, degassed), then redispersed in 50 mM Tris (+ 1 mM EDTA, degassed) to have a final concentration of about 1 mg/mL.
The dispersion of phages in presence of TCEP was then stirred (400~rpm) at 4$^{{\rm o}}$C during 48 hours.\cite{Ng,Derda2}
After the reduction reaction, viruses were purified by a centrifugation step (13.5~kg during 20 minutes) followed by 6 rounds of 10~min dialysis against PBS buffer (pH 6.5, I=350~mM + 1~mM EDTA, degassed). Maleimide-PEG2-biotin (Thermo Scientific) was dissolved in distilled water (Maximum solubility: 25~mg/mL). A 1000-fold molar excess was taken per virus and added drop-wise to the suspension. The mixture was kept under inert atmosphere thanks to the initial bubbling of Argon, protected from light and stirred during 2 hours at room temperature.
The resulting phages after biotinylation, called M13C7C-B, were then dialyzed against Bis-Tris (I=20 mM pH 7.4, 0.02\% NaN3) and purified by rinsing with the same buffer in a 100k MWCO filter tube (Amicon Ultra, Merck Millipore).\\

The fluorescent labeling of the M13C7C-B was achieved
in PBS buffer (pH 7.8, I = 350~mM) at a concentration of 0.2~mg/mL. An excess of 3 dyes (either Alexa488-NHS ester activated, Molecular Probes, or Fluorescein-NHS ester) per p8 main coat protein (about 3000 proteins per virus) first dissolved in N,N-dimethyleformamide (DMF) was added to the suspension. Particular care was taken that the proportion of DMF does not exceed 20\% of the total volume in order to avoid any denaturation of the phages by the organic solvent. The mixture was protected from light and kept under stirring at room temperature during 2~hours. The resulting viruses
were then extensively dialyzed against Bis-Tris buffer (pH 7.4, I=20~mM).

To construct multiarm structures with M13C7C-B as well as to estimate the binding affinity of M13AS with streptavidin, the viruses were mixed with 30~nm streptavidin coated iron oxide nanoparticles (Ocean NanoTech) for 14 hours at room temperature. The number of available streptavidin proteins per bead for conjugation with phages, $q$, is estimated to be of about 8 (maximum mean valency observed). Rather large nanoparticles have been chosen in order to make easier their counting by transmission electron microscopy (TEM) experiments. In particular, non-specific interactions between viruses can then be easily discriminated, especially when compared with the use of free molecular streptavidin whose relative small size makes tricky its observation by TEM. Note that some streptavidin release occurred with time due to nanoparticle aging (life time of 3 months).

Samples were then observed by TEM on a Hitachi H-600 microscope operating at 75~kV and images were recorded with an AMT CCD camera. The diluted virus suspension (typically 10$^{-2}$ - 10$^{-3}$~mg/mL) was settled onto freshly treated (O$_2$ plasma, K1050X, Quorum Technologies)
200-mesh formvar/carbon-coated grids (Agar) and stained with 2\% (w/w) uranyl acetate.

The fluorescent labeled virus particles were visualized using an inverted microscope (IX-71, Olympus) equipped with high numerical aperture (NA) oil objective (100x PlanAPO NA 1.40, Olympus) and a Neo sCMOS camera (Andor).

\section*{Associated content}

Supporting Information: DNA sequencing of the two bacteriophages, large field of view images by TEM of the control experiment with raw viruses and of the multiarm structures, and movie of multiarm colloidal molecules observed by fluorescence microscopy.

\section*{Author Contributions}
E.G. designed the project and supervised the research. C.W. and A.C. prepared the biotinylated phages and studied their self-assembly, M.T. prepared the samples and performed the experiments on M13AS viruses, C.W. and A.R. performed the fluorescence microscopy experiments. E.G. and A.C. developed the theoretical model, and wrote the manuscript.

\begin{acknowledgments}

This research was supported by the French National Research Agency (ANR) through the project AURORE. We thank
Pr S.-W Lee
for the generous gift of M13AS and P. van der Schoot for useful discussion. A.R. acknowledges support from the European Union's Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie Grant Agreement No 641839.

\end{acknowledgments}


Study of Silicon Photomultiplier Performance at Different Temperatures

Abstract: Decreasing the operation temperature of a Silicon Photo-Multiplier (SiPM) leads to a drop in its dark noise. Some experiments consider cold temperatures as an option for low noise applications of SiPM. One of those is the TAO detector, which requires operation at $T\approx -50~^\circ$C. A significant dependence of the Photon Detection Efficiency (PDE) of a SiPM on different temperatures was reported with a drastic drop around this temperature. In this paper, we present studies of performance for two samples of SiPMs from Hamamatsu and AdvanSID(FBK) companies in a broad temperature range. No significant difference for the PDE was observed.

\section{Introduction}
\label{sec:introduction}
A Silicon photomultiplier (SiPM) is an array of small avalanche photodiodes, called pixels, that operate above the breakdown voltage (overvoltage (OV)) in a Geiger mode. The breakdown voltage is the minimum reverse bias voltage that leads to self-sustaining avalanche multiplication in every SiPM's pixel. Such a structure allows for detecting and working with very low-intensity light in a proportional mode with a range limited by the total number of pixels. SiPMs are being introduced to many different applications. Nowadays they are widely used in experimental physics and medical diagnostic techniques. The main advantages with respect to conventional photomultipliers are higher detection efficiency, less sensitivity to magnetic fields, compact dimensions, and gain of up to 10$^6$. On the other hand, there are drawbacks of SiPMs: high dark rate, optical crosstalk and temperature dependence of some parameters.

Currently, a satellite experiment for the JUNO neutrino experiment called TAO is under consideration\cite{tao1}. It is proposed that TAO will precisely measure the primary neutrino spectrum from the nuclear reactors. The TAO detector is a spherical vessel filled with a liquid scintillator covered by photosensors. To achieve the required energy resolution of less than  $2~\%$ at 1 MeV for the TAO experiment one needs to maximize the light detection efficiency and minimize the noise level. SiPM is a good candidate as a photosensor for TAO since it has a higher photon detection efficiency(PDE) compared to conventional PMTs. To reduce the high noise level of the SiPM it is proposed to operate the TAO detector at low temperature $T\approx -50~^\circ$C~\cite{tao2}. But there is a relevant study \cite{italian} where a significant drop of the PDE for 1$\times$1 mm$^2$ FBK-IRST SiPM type T6-V1-PD  at $T\approx-50~^\circ$C$~(\approx220~$K) was observed (see \cref{fig:fig1}). Such a drop may totally wash out the advantage of SiPM usage. Assuming the effect might exist for other SiPM models we have studied two SiPM samples: Hamamatsu MPPC-S13360-6025CS~\cite{ham} (single element), AdvanSID ASD-NUV4S-P-4x4TD~\cite{FBK} (SiPM array).

The paper is organized as follows. In \cref{sec:setup} we describe experimental setup and methods that we used in our study. Next \cref{sec:results} presents results for both SiPM samples at different temperatures. And finally, in  \cref{sec:summary} we draw our conclusions.

\begin{figure}[!h]
   \centering
  \includegraphics[width=230px]{italian01.jpg} \\
   \caption{Normalized photon detection efficiency as a function of temperature for fixed overvoltage of 2V and three wavelengths (300 nm, 500 nm, and 800 nm) for 1$\times$1 mm$^2$ FBK-IRST SiPM type T6-V1-PD~\cite{italian}.}
  \label{fig:fig1}
\end{figure}

\section{Experimental setup}
\label{sec:setup}

In our measurements, a Dewar's vessel is used as a cryostat to keep Liquid Nitrogen (LN). Nitrogen vapor produces a gradient of temperature at different heights inside the Dewar. By moving the SiPM one can provide different temperatures in a broad range from room conditions down to LN environment. To prevent any light leakage into the Dewar's vessel it is placed inside a dark room that is completely insulated from exterior light.

\subsection{Incident light stability}

To provide a valid study one needs to guarantee that light intensity is stable in time and with environmental changes. In our measurements, we use a stabilized LED light source with $\lambda=$ 425~nm wavelength produced by the HVSys company~\cite{HV_sys}. The light from the LED is transported to the SiPM by means of a plastic optical fiber. The light traveling through the optical fiber may vary at level of about 10~\% with temperature changes. It might be driven by changing the optical properties between the fiber's core and the cladding. To reduce such a behavior, temperature stabilization is implemented along the fiber's length.

\begin{figure}[htbp]
  \begin{minipage}[ht]{0.29\linewidth}
  \label{fig:fig_setup_a}
  \centering
   \subcaptionbox{}
          {\includegraphics[width=\linewidth]{dewar.png}}
  \end{minipage}
  \hfill
  \begin{minipage}[ht]{0.69\linewidth}
    \label{fig:fig_setup_b}
  \centering
   \subcaptionbox{}
          {\includegraphics[width=\linewidth]{assembly.pdf}}

  \end{minipage}
  \hfill
  \caption{Schematic view of the experimental setup: a - Dipped probe into the Dewar's vessel, b - Probe schematic with SiPM's bed and light guides  are mounted in a copper tube.}
  \label{fig:fig_setup}
\end{figure}

\subsection{Setup description}
The cryostat vessel of $\sim$30 liters is filled to about 1/3 with Liquid Nitrogen (see \cref{fig:fig_setup}a). An assembly of a light delivery system  (see \cref{fig:fig_setup}b) with the SiPM is moving along the cryostat depth. The light delivery system is an optical fiber bundle of 6 fibers which is placed inside a copper pipe that screens it from temperature changes. The light is sent to the SiPM or reflective plate (calibration) through the central fiber. Then the light reflects back from the SiPM and travels to the other fibers which are used to monitor light stability with temperature changes. The pipe is wound by a heating cable with feedback provided by a thermal sensor that is placed inside the pipe. The heating cable keeps the fiber bundle at a stable temperature. The assembly is additionally thermally insulated on the outside.

A reflective plate that replaces the silicon dice (SiPM) for the light stability calibration is placed in front of the fiber bundle in order to reflect the light back through the other 5 fibers. In this measurement, we use quite high light intensity which is enough to monitor the reflected light by means of a PIN-photodiode ThorLabs FDS100-CAL~\cite{thorlabs} and a Keithley 6487 Picoammeter/Voltage source. The measurement gives an estimation of the light stability for different environmental temperatures. Overall light stability in our system is reached at the level of 1\%.

The reflection plate is then replaced with a single SiPM which is stuck to the aluminum substrate (bed) with a thermal grease with an embedded thermosensor on its backside. For the AdvanSID SiPM array we use a copper tub. The SiPM and thermosensor are placed inside the tub which is wound with heating cable(see~\cref{fig:tub}). In both cases, the temperature stability is guaranteed with precision better than $1^\circ$C.

\begin{figure}[!h]
   \centering
  \includegraphics[width=300px]{copper_tub.png} \\
   \caption{Copper tub is wound with heating cable and PCB with ASD-NUV4S-P-4x4TD}
  \label{fig:tub}
\end{figure}

\subsection{Readout chain and data acquisition}
A readout chain in our measurements consists of a SiPM power supply source, an amplifier and a high-speed digitizer(see~\cref{fig:readout_chain}). Apart from the readout chain, two power supply sources are utilized to supply the heating cables. All of this equipment is located outside the darkroom, except the warm amplifier that sits close to the Dewar.

To provide a bias voltage to the SiPM, the Keithley 6487 is used. The source automatically measures the current-voltage (I-V) curves of the SiPM that are used for a breakdown voltage determination.

A signal from the SiPM is digitized by means of DRS4 \cite{drs} evaluation board that can sample a signal within a range of (0.7-5)~GS/s with about 11.5 effective bits precision. It operates in external-trigger mode and is synchronized with the LED trigger signal. The saved data are analyzed by software using the ROOT package \cite{root}. We put a high-gain amplifier into the readout chain  to increase the sensitivity of the ADC by 225 times in order to distinguish single photoelectrons.

\begin{figure}[!h]
   \centering
  \includegraphics[width=270px]{readout_chain.png} \\
   \caption{Readout chain circuit.}
  \label{fig:readout_chain}
\end{figure}

\subsection{Data analysis}

By integrating waveforms within a fixed window (gate) we can extract charge\footnote{A SiPM's charge spectrum usually provides better single electron resolution compared to its amplitude spectrum (see~\cref{fig:fig_spectra}b). Usage of the charge spectra allows better extraction of pedestal events.} spectrum of  a SiPM. We acquire $N$-events for signal - (LED is ON), see~\cref{fig:fig_spectra}a, and $D$-events for dark (LED is OFF) spectra. Using the number of events in the pedestal for signal $N_0$ (blue area under the pedestal peak) and for dark $D_0$ spectra (gray area given in the same range) the estimation of the average number of photoelectrons $\hat{\mu}$ can be found \cite{anfimov}
\begin{equation}
    \label{eq:mu_00}
   \hat{\mu} = - \ln{\left(\frac{N_0}{N}\cdot\frac{D}{D_0}\right)} = -\ln{\left(\frac{\hat{p}_{\xi_0}}{\hat{p}_{\lambda_0}}\right)} = \hat{\xi} - \hat{\lambda},
\end{equation}

where $\hat{p}_{\xi_0}$, $\hat{p}_{\lambda_0}$ - estimators of the pedestal probability in signal and dark  spectra,  and  $\hat{\xi} = -\ln(\hat{p}_{\xi_0})$, $\hat{\lambda} = -\ln(\hat{p}_{\lambda_0})$ - estimators of the average number of photoelectrons in signal and dark spectra, respectively.
Statistical error in this case could be evaluated \cite{anfimov}

\begin{equation}
    \label{eq:dispersion_mu_00_noise}
      \frac{\hat{\sigma}_{\hat{\mu}}}{\hat{\mu}} \approx \frac{1}{\sqrt{N}}\sqrt{\frac{e^{\hat{\xi}}+ e^{\hat{\lambda}}-2}{(\hat{\xi}-\hat{\lambda})^2}}.
\end{equation}

\begin{figure}[htbp]
  \begin{minipage}[ht]{\linewidth}
  \label{fig:fig_pde_1}
  \centering
   \subcaptionbox{}
          {\includegraphics[width=\linewidth]{charge_spectrum.jpeg}}
  \end{minipage}
  \hfill
  \begin{minipage}[ht]{\linewidth}
    \label{fig:fig_pde_2}
  \centering
   \subcaptionbox{}
          {\includegraphics[width=\linewidth]{amplitude_spectrum.jpeg}}

  \end{minipage}

  \hfill
  \caption{Typical charge and amplitude spectra of SiPMs: a - charge spectrum illustrate the pedestal method for evaluation of average number of photoelectrons. Blue area and $N_0$ - pedestal events on signal spectrum, gray area and $D_0$ - pedestal events on dark spectrum (incorporated picture); b - amplitude spectrum of SiPM and an illustration of the threshold extraction method.}
  \label{fig:fig_spectra}
\end{figure}

 To compare parameters at different temperatures one should operate SiPM with similar overvoltage. A way to find the breakdown voltage is to analyze the IV-curve of a SiPM. There are several approaches of breakdown voltage extraction. In our study, we use the Inverse Logarithmic Derivative (ILD) test described in \cite{chmill}. The ILD for IV-curve reads

\begin{equation}
\label{eq:ILD}
   f(V) =  \left[\frac{1}{I}\cdot \frac{dI}{dV}\right]^{-1}
\end{equation}
and the minimum of the $f(V)$ yields the inflection point which defines the breakdown voltage. We also test the quadratic fit for AndvanSid SiPM which was used in the paper~\cite{italian}. Both approaches demonstrate good consistency of the breakdown voltage extraction (see~\cref{table:bd}).  Measured voltages for both SiPMs using the ILD test are shown in~\cref{fig:Vbd}.

\begin{figure}[!h]
   \centering
  \includegraphics[width=350px]{Vbr_vs_temp_last.pdf} \\
   \caption{Breakdown voltages for Hamamatsu and AdvanSid SiPMs vs temperature using ILD test (see~\cref{eq:ILD}).}
  \label{fig:Vbd}
\end{figure}

\begin{table}[h]
\caption{\textbf{Breakdown voltages for the AdvanSid SiPM at different temperature points extracted by ILD test and quadratic fit}}
 \label{table:bd}
\centering
 \scalebox{0.8}{
 \begin{tabular}{c | l l l}
 \hline
 \hline
\bf{SiPM type}  &  \multicolumn{3}{c}{\bf{Breakdown voltage $V_{bd}$ at different temperatures}} \\
 \hline
\multirow{2}{*}{ASD-NUV4S-P-4x4TD (ILD test)}  & 26.57 V[ $23~^o$C] & 26.02 V[ $2~^o$C] & 25.39 V[$-22~^o$C] \\
\cline{2-4}
  & 24.58 V[$-52~^o$C] & 23.83 V[$-84~^o$C]\\
\hline
\hline
\multirow{2}{*}{ASD-NUV4S-P-4x4TD (quadratic fit)} & 26.55 V[ $23~^o$C] & 26.00 V[ $2~^o$C] & 25.39 V[$-22~^o$C] \\
\cline{2-4}
  & 24.57 V[$-52~^o$C] & 23.84 V[$-84~^o$C]\\
\hline
\hline
\end{tabular}}
\end{table}

 We evaluated the Dark Count Rate (DCR) of the SiPM using the following procedure. By processing raw waveforms we count a number of dark pulses (LED is OFF) that exceed some fixed threshold. We processed $w = 10000$ waveforms within $W = 1.024~\mu$s duration window\footnote{In our case, SiPM pulses $\sim100$~ns (FWHM) which can introduce a small fraction of pulses miscounting. We took it into account shrinking the window by the pulse width.}  (1024 samples at 1~GS/s sampling rate).  The overall counting window is $w \times W$. The DCR rate $R$ could be found from the number of dark pulses $N$ that have been counted in the overall window as
 \begin{equation}
    R = \frac{N}{w\times W}
    \label{eq:dcr}
\end{equation}
 The number of the random pulses $N$ follows Poisson distribution\footnote{Standard deviation $\sigma_{N} = \sqrt{N}$} and a standard deviation of the DCR using \cref{eq:dcr} is

  \begin{equation}
     \sigma_R = \frac{\sqrt{N}}{w\times W}.
    \label{eq:err}
\end{equation}

 The threshold is obtained as an empirical compromise between efficiency and noise pick up. By measuring SiPM's amplitude spectrum (see.~\cref{fig:fig_spectra}b) we defined the optimal threshold\footnote{We assume signal and noise pulses are of similar amplitude.} $T$ as
\begin{equation}
    \text{T} = A_1 - 2\sigma_{A_1},
    \label{eq:thr}
\end{equation}
where $A_1$ and $\sigma_{A_1}$ - the amplitude and its standard deviation of the single photoelectron signal, respectively.  Decreasing the threshold leads to significant noise pick up at a small gain (overvoltage). In our case (see \cref{eq:thr}) only $\sim$2,5~\% of SiPM's signals are miscounted in the dark rate evaluation procedure.\footnote{Assuming single pixel's amplitudes follow Normal distribution}

 \begin{figure}[ht]
  \begin{minipage}[ht]{0.5\linewidth}
  \label{fig:fig_pde_1}
  \centering
   \subcaptionbox{}
          {\includegraphics[width=\linewidth]{pde_ham_datasheet.pdf}}
  \end{minipage}
  \hfill
  \begin{minipage}[ht]{0.5\linewidth}
    \label{fig:fig_pde_2}
  \centering
   \subcaptionbox{}
          {\includegraphics[width=\linewidth]{pde_fbk_datasheet.pdf}}

  \end{minipage}
  \hfill
  \caption{PDE as a function of wavelength: a -- for Hamamatsu SiPM at overvoltage OV = 5 V\cite{ham}, b -- for AdvanSid SiPM at OV = 6 V\cite{FBK}}
  \label{fig:fig_pde}
\end{figure}

Measurements were performed for both SiPM samples at different temperatures and in the range of operating overvoltages. To present the PDE more illustratively we normalized measured relative PDE ($\hat{\mu}$) to 25~\% (see \cref{fig:fig_pde}) for Hamamatsu SiPM  ($V_{OV}=5.0~V$ at $T=23.0~^\circ$C) and to 42~\% \footnote{On \cref{fig:fig_pde}b the PDE is a bit higher than 40~\% at $V_{OV}$=6~V. For our sample, we were not able to reach $V_{OV}$ above 5.68~V because of high noise level. And we assumed the PDE=42~\% at this point.} for ASD-NUV4S-P-4x4TD ($V_{OV}=5.68~V$ at $T=23.0~^\circ$C).

\section{Results}
\label{sec:results}

The PDE and the DCR dependencies on overvoltage for two SiPMs were studied at the wide range of temperatures from $T=23.0~^\circ$C (room temperature) and even below the $T=-52.0~^\circ$C which is close to the TAO operational temperature.

To compare the results of the PDE behavior for all temperatures we plot the PDE as a function of overvoltage (see.~\cref{fig:fig_ov}). One can see that there is no significant temperature dependence of the PDE for both samples.

The DCR for SiPMs(see.~\cref{fig:fig_dcr}) decreases rapidly and reaches a drop by a factor of several orders of magnitude for low temperatures compared to room conditions. This fact is one of the major considerations for the SiPM usage at $T=-50~^\circ$C for the TAO detector. As can be seen from~\cref{fig:fig_dcr} for low overvoltages at negative temperatures we cannot obtain DCR values\footnote{Small number of dark pulses comparing to a noise pick up at a very low threshold}.

\begin{figure}[htbp]
  \begin{minipage}[ht]{0.5\linewidth}
  \label{fig:fig_pde_1}
  \centering
   \subcaptionbox{}
          {\includegraphics[width=\linewidth]{pde_vs_ov_ham_update_last.pdf}}
  \end{minipage}
  \hfill
  \begin{minipage}[ht]{0.5\linewidth}
    \label{fig:fig_pde_2}
  \centering
   \subcaptionbox{}
          {\includegraphics[width=\linewidth]{pde_vs_ov_fbk_last.pdf}}

  \end{minipage}
  \hfill
  \caption{PDE as a function of overvoltage: a -- for Hamamatsu SiPM, b -- for AdvanSid SiPM}
  \label{fig:fig_ov}
\end{figure}

\begin{figure}[htbp]
  \begin{minipage}[ht]{1.0\linewidth}
  \label{fig:fig_dcr1}
  \centering
   \subcaptionbox{}
          {\includegraphics[width=\linewidth]{dcr_vs_ov_ham_v2.pdf}}
  \end{minipage}
  \hfill
  \begin{minipage}[ht]{1.0\linewidth}
    \label{fig:fig_dcr2}
  \centering
   \subcaptionbox{}
          {\includegraphics[width=\linewidth]{dcr_vs_ov_fbk_v5.pdf}}

  \end{minipage}
  \hfill
  \caption{Evaluated DCR as a function of overvoltage: a -- for Hamamatsu SiPM, b -- for AdvanSid SiPM}
  \label{fig:fig_dcr}
\end{figure}

\section{Summary}
\label{sec:summary}
In general, the setup and technique demonstrated in this paper allows for measurements with different types of SiPMs in a wide range of temperatures from Liquid Nitrogen (77~K) up to the Room temperature.
We studied the Photon Detection Efficiency of SiPMs in a broad temperature range including $T = -52.0~^\circ$C ($\sim 220$~K) which is in the considered region of the TAO detector. As a result, two types of SiPMs - Hamamatsu S13360-6025CS and AdvanSID ASD-NUV4S-P-4x4TD have no significant difference in the PDE. These measurements do not exclude the behavior of the PDE in a broader temperature range and demonstrate an absence of difference for particular temperature points.\footnote{And incident light wavelength $\lambda=$425~nm.}

The result presented in the study \cite{italian} that demonstrates the PDE drop in the temperature region around $T\sim$ 200~K to 220~K have not been found for our SiPM samples. Therefore, the SiPM usage in the TAO experiment is relevant and might have advantages comparing to the conventional PMTs. The choice of the proper type of SiPM for TAO requires a complex study of the SiPM's parameters at the TAO operational temperature.

\section{Acknowledgments}
The authors are indebted to Dr. Dmitriy Naumov, prof. Alexander~Olshevskiy and Christopher Kullenberg from the Joint Institute for Nuclear Research for their great help in the preparation of this paper.
\label{sec:acknowledgments}

Instrumentation and Detectors (physics.ins-det); High Energy Physics - Experiment (hep-ex)


YSO disk structure and planetary signatures

Abstract: Detailed radial/vertical structure modeling of observations of disks around Young Stellar Objects (YSOs) can provide information on the physical conditions and on the characteristics of the gas and dust in their interiors. We describe recent results of self-consistent modeling of spectral energy distributions, optical and infrared images, and millimeter fluxes of YSOs. We discuss observations and interpretations of the initial stages of planet formation, including indications of dust growth and settling in disks around young stars. We also discuss how the unprecedented resolution and sensitivity of ALMA may help us study the interior of the innermost disks, a region unaccessible with present day instrumentation, and witness the very first stages of planet formation.

\section{Introduction}

Since the late 1970s, it was known that
the low mass pre-main sequence  T Tauri stars were surrounded
by disks. For instance,
the lack of correlation between the (near infrared) excess and the
extinction had led to the suggestion that the cool material producing the excess
was located in a geometrically flat configuration
(Cohen \& Kuhi 1979).
However, it was not until the first {\it IRAS} observations of TTS appeared
that the large mid and far infrared excesses were attributed to
accretion disks (Rucinski 1985), in agreement with the original
suggestion of Lynden-Bell \& Pringle (1974). Masses and radii for these
disks were estimated from submillimeter and millimeter observations
around the late 1980's. The indirect inference for disks from all
these observations was dramatically confirmed in this
decade, when disks around young stars have been directly
imaged by instruments on board of {\it HTS}. Disks are seen
against the bright background of the Orion nebula
as silhouettes (McCaughrean \& O'Dell 1996; McCaughrean et al. 1998;
Stauffer et al. 1999), and in the favorable edge-on situations,
where the direct light from the star is blocked by the
disks, so the faint glow of stellar light scattered by dust
grains on the disk surface can be seen.
A number of cases have been presented so far;
the best known and analyzed include HH30 (Burrows et al. 1996)
and HK Tau/c
(Stapelfeldt et al 1998).
These disks are the predecessor of planetary systems.

T Tauri disks are composed of gas and dust. The mass in the
disk is of the order of a few 0.01 $M_{\sun}$, and matter is
still
accreting onto the star. On the other extreme, in disks with planetary
systems,  most of the mass is on the planets, a few 0.001 $M_{\sun}$,
and only small amount of dust ($<< 10^{-6} M_{\sun}$)
remains.
The transition between these
two stages is of upmost interest. What are the
processes that lead to planet formation and in which timescales?

In this contribution, we argue that the
first signs of the transition to planetary disks
appear in the inner disk, within a few AU from the
star, and then when this
transition starts, it happens very fast. We will
review evidences provided by two groups of
different ages, 1-2 Myr and 10 Myr. We then will
argue that present and future observations in
the optical and infrared cannot tell us what
is happening near the midplane of the inner regions of the gas and
dust disks characterizing the initial stages before
the transition. Only ALMA, in the highest resolution
configuration, will allow us to scrutinize the  innermost
disk and to obtain information on the processes
occurring to the dust in these regions leading
to the formation of planets.

\begin{figure}[ht]
\plotfiddle{calvetn_fig1.ps}{2.5in}{0}{50.0}{40.0}{-150.0}{-70.0}
\end{figure}

\noindent
Fig. 1. Fluxes of classical T Tauri stars in Taurus-Auriga.
From D'Alessio et al. (1999b.) Fluxes are normalized at H.
The solid line is the mean SED constructed from these data,
including upper limits at 100$\mu$m. The dashed line does
not include these upper limits.

\section{Indications of dust evolution}

Figure 1 shows fluxes for all classical T Tauri stars (CTTS)
in Taurus-Auriga, scaled at H (D'Alessio et al. 1999b).
The large infrared excesses due to disk emission
are apparent.
Disk emission depends on the temperature of the disk, which
in turn, is determined by the amount of heating. Two agents
can heat the disk: viscous dissipation related to accretion, and
stellar irradiation. For the typical values of the mass accretion
rates in the disks around CTTS (Gullbring et al. 1997; Hartmann et al.
1998), irradiation from the central star plays a
dominant role in the heating of the disk (Kenyon
\& Hartmann 1987, KH87;
D'Alessio et al. 1998, 1999a,b,c).
Heating by irradiation
will depend on the shape of the disk surface, since
disks with flared surfaces can capture more
stellar radiation
than flat disks (KH87). In turn,
the scale height of the disk, and therefore
the height of the surface, depends on the disk temperature.
A self-consistent stability analysis of irradiated
disks for typical CTTS parameters indicates that
a unique, stable solution for the surface exists,
and corresponds to a flared surface
(D'Alessio et al. 1999a).
This means that
disk emission will be determined essentially by
how much energy comes into the disk, set by
the star, and how much of this energy is actually
absorbed and re-emitted by the disk, which depends
on the dust in the disk.

Most stars in Taurus
shown in Figure 1 have spectral types between K7 and M3,
so they have similar amounts of stellar irradiation.
The large dispersion observed in the SEDs must
essentially be due then to differences in the dust properties
in the disk, indicating that even though the stars
have similar ages (1 - 2 Myr, Gomez et al. 1993),
there is a
large diversity in the conditions of the dust in these disks.

Figure 2 shows the SEDs for two ``transition
cases'' in Taurus: DI Tau (Meyer et al. 1997) and
V819 Tau (M. Meyer, personal communication). These
stars are weak TTS, that is, they are not accreting
mass onto the star. The indications of disk accretion,
i.e., broad emission line profiles, UV excesses,
which characterize CTTS are
not present in these stars, nor are near-infrared
excesses. However, mid-infrared and even millimeter
flux excesses persists, indicating that remnants of the disk
are still there. Some other transition cases may exist
and remain to be
discovered, but not in large numbers; the total
number of transition cases may be less than 10, which when compared
to the total number of TTS in Taurus, $\sim$ 100,
indicates that the loss of the disk happens fast.
The actual processes causing disk dissipation
are not well known, and
binary companions may play a crucial role in
clearing the disks; but in any event, the processes
leading to disk dissipation occur in very short
timescales, compared to the lifetime of a typical T Tauri star.

\begin{figure}[ht]
\vskip -1in
\plotfiddle{calvetn_fig2.ps}{2.5in}{0}{50.0}{40.0}{-150.0}{-70.0}
\end{figure}

\noindent
Fig. 2. SEDs of two ``transition'' T Tauri stars in Taurus:
DI Tau (right) and V819 Tau (left) (data from Kenyon
\& Hartmann 1995). These are
weak T Tauri stars, where accretion has already stopped,
but still with mid and far infrared excesses.
The median SED in Taurus (solid line) and
photospheric fluxes (dashed line) are shown for comparison.

\vskip 0.3in

Figure 3 shows SEDs for stars in the TW Hya
association, which has an estimated age of
10 Myr (Webb et al. 1999) from Jayawardhana et al. (1999).
Although TW Hya itself has been known for many
years (Herbig 1978; Rucinski \& Krautter 1983;
de la Reza et al. 1989) as an ``isolated T Tauri star'', it
is not until recently that this extremely interesting
association has begun to be extensively studied.
The main problem has been finding its members, since
molecular gas is no longer associated with the
association.
Most of the members found are turning out to be
weak TTS, but a few of them still show indications
of disk emission, similar to the mean SED
in Taurus in the mid and
far-IR and in the millimeter (Figure 3). Again, in some cases,
a deficit in the near-IR compared to longer wavelengths
seem to be present.

The evidence provided by young
stars in the two extremes of the temporal range
where disks are expected to dissipate and planets
to form ($\sim$ 10 Myr, Podosek \& Cassen 1994)
indicates that disks can last for a very long
time, but once clearing starts, it occurs in very short
time-scales. The first signs of disk clearing
seem to appear in the near-IR, indicating that
the higher temperature regions in the inner disk
are the first to clear up. This is is consistent
with the evidence provided by debris disks; for
instance, in the youngest debris disk known,
around the A0 star
HR 4796A, $\sim$ 8 Myr (Jayawardhana et al. 1998;
Koerner et al. 1998), a possible member of the TW Hya association,
disk emission is only apparent beyond 10 microns.
The inner disk has been cleared out maybe
by a planet or planets, consistent with
the images in the mid-IR and in the near-IR from
{\it HST} (Telesco et al. 1999; Schneider et al. 1999).

\vskip 0.3in

\begin{figure}[ht]
\plotfiddle{calvetn_fig3.ps}{2.5in}{0}{50.0}{40.0}{-150.0}{-70.0}
\end{figure}

\noindent
Fig. 3. SEDs of stars in the TW Hya association.
From Jayawardhana et al. (1999).
The solid curve is the median SED in Taurus,
and the dashed lines indicate the quartiles.
The dotted line indicates the photosphere.

\vskip 0.3in

The process of inner disk clearing must start
with dust growth and settling towards the midplane.
Eventually, a planet (or planets) forms
in the disk, and
a gap is opened at the position of the
planet; the inner disk is accreted onto the
star, while a considerable fraction of the
remaining mass of the disk goes into planet
formation. Disks around CTTS
must be in the very first stages of this
sequence; once planets form
and gaps appear, the evolution happens very
fast, and it may prove difficult to find
disks with gaps, still retaining the
inner disk.
It may bear more results to study the early phases,
previous to planet formation,
when the star still retains
its original gas plus dust disk,
searching for indication of dust coagulation
and settling.

\vskip 0.3in

\begin{figure}[ht]
\plotfiddle{calvetn_fig4a.ps}{1.5in}{0}{25.0}{25.0}{-80.0}{-30.0}
\plotfiddle{calvetn_fig4b.ps}{1.5in}{0}{25.0}{25.0}{-150.0}{-60.0}
\plotfiddle{calvetn_fig4c.ps}{1.5in}{0}{25.0}{25.0}{-0.0}{60.0}
\end{figure}

\vskip -0.9 true in

\noindent
Fig. 4. Upper panel: (a) Disk cumulative flux as a
function of radius at three wavelengths:
10 $\mu$m (solid), 0.5 mm (dotted) and 1 mm (dashed).
The disk model has typical
parameters for CTTS: mass 0.5 $M_{\sun}$, radius 2 $R_{\sun}$,
mass accretion rate $10^{-8} M_{\sun} \, yr^{-1}$, and
effective temperature
4000K.
Lower left: (b) Characteristic optical depths in the disk:
$\tau(1\mu$m), where stellar radiation is absorbed (dotted);
$\tau_{Ross}$, the optical depth at the local radiation (solid).
The disk is optically thick to the stellar radiation
even in the outer regions where it becomes optically
thin (to its own radiation).
Lower right: (c) Characteristic heights in the disk:
$z_s$, surface where $\tau(1\mu$m) $\sim$ 1 (solid);
the disk scale height (dashed);
$z_{phot}$, surface where $\tau_{Ross} \sim$ 1
(defined where the disk is optically thick)(long dashes).

\section{The Tauri disks}

Self-consistent radial/vertical calculations of disk
structure for typical CTTS parameters, including viscous
and irradiation heating, have been carried out by D'Alessio
et al. (1998, 1999b,c). These calculations show that
the flux at $ \lambda \le$ 10 $\mu$m, where the first
indications of disk clearing appear, comes from the region
inside a few AU of the disk. Figure 4a shows the
cumulative flux at three wavelengths, normalized
at the maximum, indicating that only in the millimeter
range the whole disk is contributing to the flux.

\begin{figure}[ht]
\plotfiddle{calvetn_fig5a.ps}{3.0in}{0}{35.0}{43.0}{-180.0}{-110.0}
\plotfiddle{calvetn_fig5b.ps}{1.5in}{0}{30.0}{28.0}{20.0}{111.0}
\end{figure}

\vskip -1.8in
\noindent
Fig. 5. Left: {\it HST} scattered light model images
at 0.8$\mu$m for edge-on
disks, corresponding to dust well mixed with gas and
different size distribution and maximum size (see text).
Right: SEDs for models with images in the left size.
$a_{max} = 0.3 \mu$m (solid), 10 $\mu$m (dashed), and 100 $\mu$m
(dotted).  From D'Alessio et al. (1999c).

\vskip 0.3in

Figure 4b shows characteristic optical depths in the
disk. The optical depth at $\lambda \sim 1 \mu$m,
at the peak of the stellar radiation, is always
much larger than 1,
since dust opacity is high at optical and
near-IR wavelengths. In contrast, the Rosseland
mean optical depth, which corresponds
to the optical depth
to the local radiation, falls below 1 beyond $\sim$
10 AU. This means
that the disk is thick
to stellar radiation even when it is optically
thin to the local radiation (for instance, in the outer disk.)
This last property implies that the disk
can efficiently capture stellar radiation
up to a few hundred AU, of the order of
the estimated radii of CTTS (Beckwith et al. 1990;
Dutrey et al. 1996). It also implies that
scattered light images of edge-on disks
observed with {\it HST} should show the flared
shape characteristic of the surface where
stellar irradiation is being captured, as
it is indeed the case. Figure 4c shows this
surface, $z_s$, as well as the scale height
and $z_{phot}$, the
height of the optically thick disk (obviously defined
only in the region where the disk is thick).

\vskip -1.5in

\begin{figure}[ht]
\plotfiddle{ calvetn_fig6a.ps}{4.0in}{0}{45.0}{35.0}{-200.0}{-28.0}
\plotfiddle{calvetn_fig6b.ps}{1.5in}{0}{30.0}{28.0}{20.0}{111.0}
\end{figure}

\vskip -1.9in
\noindent
Fig. 6. Left: Model scattered light {\it HST} image at 0.9 {\AA}
for typical CTTS parameters.
Right: isocontours of image on left (solid), superimposed
to disk emission observed with ALMA
at 0.9 mm with
0.01" resolution (dotted).

\vskip 0.3in

{\it HST} scattered light images and
SEDs of CTTS may
already be given indication of dust growth
in the disk. Figure 5 shows calculations
of images and SEDs for
a mixture of astronomical silicates
and graphite (Draine \& Lee 1987, DL87), with
a size distribution $n(a) \propto a^{-p}$,
with $p =$ 3.5 (MRN, typical of interstellar medium dust,
DL 87) and $p =$ 2.5 (indicating some degree of
coagulation, Miyake \& Nakagawa 1995). The
minimum size of the dust is $a_{min} = 0.005 \mu$m (DL87)
and the maximum size is $a_{max}$ = 0.3 $\mu$m (DL87, ISM
dust),
10 $\mu$m,  and
100 $\mu$m, for fixed dust to gas ratio. As the
maximum grain size increases, the opacity decreases in the
optical and near-IR, as more of the dust is locked
in bigger grains; at the same time, the opacity
increases in the millimeter range. The effect
of grain growth on scattered light images is
shown in Figure 5, left. ISM dust produces images
where the widths of the
lanes separating the bright nebulae above and below the
disk are of the order of $\sim$ 100 AU (at Taurus),
which is larger than the width of the
observed edge-on disks (D'Alessio et al. 1999b).
The width of the lane decreases as $a_{max}$
increases, as a result of the decreased opacity
in the optical and the consequent lowering of the
surface where $\tau (1 \mu$m) $\sim 1$, producing
better matches to the observations.
At the same time,
mixtures with larger grains produce much better
agreement between the typical CTTS disk model
fluxes and
the median Taurus SED (Figure 5, right).

\section{Scrutiny of the deep inner disk: only with ALMA}

{\it HST} observations, even at $\sim$ 0.05" resolution,
can only probe the upper regions of the disk, because
the CTTS disk is {\it always optically thick to optical and
near-IR radiation}. Only submillimeter and
millimeter observations can reach the mid-plane of the
disk. But present day interferometers can
only give us a resolution of at most $\sim$ 40 AU at the
closest star formation regions. This implies
that we have no direct
observations to constraint the detailed structure of the
regions near the midplane in the inner disk
regions. ALMA at the highest
resolution, $\sim$ 0.01", will constitute the
perfect complement to {\it HST} observations, allowing us
to reach the true scale height of the disk.
To illustrate this,
Figure 6 compares typical CTTS model calculations for
a scattered light {\it HST} edge-on disk image at 0.9 {\AA} and
the corresponding ALMA image at 0.9 mm with
0.01" resolution.

Moreover, with the highest ALMA resolution, 0.01" $\sim$ 1.4 AU
at Taurus, we will be able to probe different
heights at the innermost disk regions observing
at different wavelengths. Figure 7, left, shows
the heights in the disk where $\tau \sim 1$
at different wavelengths. In the submillimeter
range and beyond, the disk opacity is low enough
to allow inspection of regions closer and closer
to the midplane. These observations will allow
us for the first time to {\it effectively map the inner disk
structure}.

\vskip -0.8in

\begin{figure}[ht]
\plotfiddle{calvetn_fig7a.ps}{2.5in}{0}{30.0}{30.0}{-180.0}{-80.0}
\plotfiddle{calvetn_fig7b.ps}{2.5in}{0}{30.0}{30.0}{20.0}{111.0}
\end{figure}

\vskip -1.9in
\noindent
Fig. 7. Left: Height in AU at which the optical depth
becomes 1 at given $\lambda$, i.e., height of formation
of radiation at $\lambda$, for the case
of gas and dust well mixed and $a_{max} = 100 \mu$m.
Light solid: 0.5$\mu$m,
dotted: 0.5 mm, dashed: 1 mm.
Right: Brightness temperature at different $\lambda$,
corresponding to heights in left. The heavy solid
line is the midplane temperature.

\vskip 0.3in

The temperature in the inner disk, $T(z,R)$
depends on local conditions. Using the diffusion
approximation, $dT/dz \propto \kappa F$. The
opacity $\kappa$ will depend on the dust
properties, type, size, and vertical distribution.
The flux generated at each height, $F$, will
depend on the local viscous generation
of energy. Figure 7, right,
shows the brightness temperature
at several wavelengths,
for the case
of dust with $a_{max} = 100 \mu$m, well mixed
with gas, and viscous heating at each height.
Because of the large opacity,
the temperature in the mid-plane is much
higher than the surface temperature. Brightness temperatures
in the millimeter are consequently much higher than
those obtained if the disk were isothermal.
However, the actual brightness temperatures
may be very different
if there is a vertically-dependent dust opacity,
due to dust coagulation and settling.
Moreover, they may be very different if
a ``dead zone'' exists near the mid-plane,
where viscous generation of energy is not
occurring because there is not enough
ionization to maintain the magneto rotational
instability (MRI, Gammie 1996; Glassdgold et al. 1997).
The importance of this dead zone for planet
formation is that inhibition of the MRI
produces an increase of the
density near the midplane (since matter cannot
be transported inwards towards the star),
favoring gravitational instabilities and
potentially leading to planet formation.

\section{Conclusions}

The present evidence from disks around 1 to 10 Myr
old stars indicates that dust growth and
settling leading to planet formation and disk clearing happens fast
and starts near the midplane in the few inner AU of the disk.
Since the inner disk regions are optically thick
to radiation shortwards of $\sim$ 0.5 mm, only ALMA
in the highest resolution configuration, corresponding to
0.01" or 1.4 AU at the closest star-forming regions,
will allow us to reach the midplane of the planet-forming
regions in the disks and study the
processes that precede the formation of these bodies.

\acknowledgements
This work was supported in part by NASA grants NAGW 2306 and NAG5-4282.


Mottness versus unit-cell doubling as the driver of the insulating state in 1T-TaS2

Abstract: If a material with an odd number of electrons per unit cell turns out to be insulating, Mott localisation may be invoked as an explanation. This is widely accepted for the layered compound 1T-TaS2, which has a low-temperature insulating phase comprising charge order clusters with 13 unpaired orbitals each. But if the stacking of layers doubles up the unit cell to include an even number of orbitals, the nature of the insulating state is ambiguous. Here, scanning tunnelling microscopy (STM) reveals two distinct terminations of the charge order in 1T-TaS2, the sign of such a double-layer stacking pattern. However, spectroscopy at both terminations allows us to disentangle unit-cell doubling effects and determine that Mott localisation alone is enough to drive gap formation. We also observe the collapse of Mottness at an extrinsically restacked termination, demonstrating that the microscopic mechanism of insulator-metal transitions lies in degrees of freedom of interlayer stacking.

\subsection*{Methods}
Crystals of 1\textit{T}-TaS$_{2}$ were synthesized using a chemical vapor transport method described previously \cite{Tani1981}, with 2 \% excess S. Samples were cleaved in ultra-high vacuum ($\sim$10$^{-10}$ Torr) at 77 K and quickly inserted into a Unisoku 1300 low temperature STM system, of which the STM head has been replaced with a homemade one \cite{Hanaguri2006}. All measurements were performed at a temperature of 1.5 K. STM tips were formed using electro-chemical etching of tungsten wire, and after insertion into UHV, were cleaned and characterized using field ion microscopy followed by careful conditioning on a clean Cu(111) surface. STM topography images were collected in constant-current mode. For conductance curves and spectroscopic mapping on 1\textit{T}-TaS$_{2}$ surfaces, the lock-in technique with a bias modulation of amplitude $V_{mod}$ = 10 mV and frequency $f_{mod}$ = 617.3 Hz was used.

The $dI/dV$ spectroscopy data in Figs. 1 \& 2, including the spatially resolved data shown in Fig. 2b, were normalized according to the current signal at $V$ = 500 mV. This somewhat compensates for the large difference in raw signal intensity caused by differing tip-sample distances while scanning at $V$ = 250 mV over areas with large or small spectral gaps.

\subsection*{Acknowledgements}
We are grateful to Y. Kohsaka, T. Machida and P. A. Lee for helpful discussions. C.J.B. gratefully acknowledges support from RIKEN's SPDR fellowship. This work was supported in part by JSPS KAKENHI grant numbers JP18K13511, JP19H00653 and JP19H01855.

\subsection*{Author contributions}
T.H. and Y.I. conceived the project, and M.Y. and Y.I. synthesized the 1\textit{T}-TaS$_{2}$ crystals. C.J.B. performed the STM measurements with assistance from T.H., and prepared the manuscript with input from all authors.

\section{Supplementary discussion}

\subsection{I. Additional side-by-side observations of distinct terminations.}

In this work, 24 surfaces were investigated and the density-of-states (DOS) spectrum with the larger gap ($\sim$150 meV), and with the smaller gap ($\sim$50 meV) were observed 18 and 6 times, respectively, considering only the spectrum found at the landing zone of the tip in each case. At 3 surfaces the type of spectrum was seen to change upon crossing some form of boundary after moving away from the landing zone. One example, a pair of single-layer steps, is described in the main text (Figs. 2 \& 3). The others, a CDW domain wall (DW), and another single-layer step, are described here.

Domain walls -- in-plane phase discontinuities in the CDW -- were often found at pristine cleaved surfaces. In principle, such DWs can offer opportunities to investigate the effects of inter-layer stacking: If one CDW layer with a DW is overlaid atop another layer without one, the phase jump \textit{must} result in different inter-layer stacking on either side, and be accompanied by a change in the DOS if the stacking vectors on each side are symmetrically inequivalent (e.g. $\textbf{T}_{A}$ versus $\textbf{T}_{C}$). In practice however, the same type of spectrum (either large or small gap) was almost always observed on both sides of DWs. This may be because the DWs extend downwards below the surface such that the surface layer and buried layers all have matching phase jumps, or because the change in stacking from one side to another is merely from one among the group of symmetrically equivalent $\textbf{T}_{C}$ vectors to another. However, on one occasion, shown in Fig. S1 below, the two types of spectrum were observable in the same field of view, in different domains separated by an unusual DW.

\begin{figure}[!hb]
\centering
\includegraphics[scale=1]{FIGS1_v3.eps}
\caption{\label{fig:1} \textbf{Change in electronic structure with change of stacking across a buried DW.} (\textbf{a}) A topographic image showing multiple DWs found at a pristine cleaved surface ($V$ = 500 mV, $I_{set}$ = 10 pA, scale bar 20 nm). (\textbf{b}) The displacement map $\textit{u}_{r,\phi}(\textbf{x,y})$ (with displacement taken relative to a point in the bottom-right corner) highlighting the in-plane phase jumps of the CDW pattern of the surface layer at the bright DWs. At the position of the buried DW, marked by the black dotted line, there is no displacement. (\textbf{c}) Spatially resolved tunnelling spectroscopy taken along the long axis of the red tinted rectangle in (a) \& (b), crossing both top layer DW and the buried one. A change in DOS is observed only at the buried DW. (\textbf{d}) Representative tunnelling spectra acquired in Regions I--III. (\textbf{e}) A schematic cartoon of the possible stacking configurations in each of the three regions.}
\end{figure}

In Figure S1(a), ordinary DWs appear as bright linear features, along with another feature appearing as a dark line or depression (indicated by the white arrow). Here we use a partial implementation of the Lawler-Fujita algorithm \cite{Lawler2010} in order to visualize the in-plane phase discontinuities in the surface CDW pattern. Essentially the lock-in phase detection technique is applied to the image, yielding a map of the spatially varying phase discrepancy (displacement) between the observed CDW modulation and an ideal reference modulation. Figure S1(b) shows the resulting displacement field $\textit{u}_{r,\phi}(\textbf{x,y})$. (The displacement is taken with reference to the bottom-right corner of the image, and expressed in polar coordinates.) The sharp changes in hue signify in-plane phase jumps. Importantly, although phase jumps are clear for the bright DWs, none is seen for the dark meandering feature (marked by a black dotted line). We therefore interpret this dark feature as a \textit{buried} DW, with a continuous-phase layer of the CDW overlaid on top. Such a feature must be accompanied by a change of the inter-layer stacking.

\begin{figure}[!ht]
\centering
\includegraphics[scale=1]{FIGS2_v2.eps}
\caption{\label{fig:2} \textbf{Partial replication of step-terrace measurements in Figs 2 \& 3 accompanying the main text.} (\textbf{a}) A topographic image showing two terraces separated by a single-layer step ($V$ = -500 mV, $I_{set}$ = 500 pA, scale bar 5 nm). (\textbf{b}) Map of the displacement field $\textit{u}_{r,\phi}(\textbf{x,y})$ (with displacement taken with respect to a reference point in the top-left corner), showing an absence of in-plane phase jump (displacement) from the upper to the lower terrace. (\textbf{c}) Tunneling spectroscopy along the path marked by the red-tinted rectangle in (a), averaged over the rectangle's short axis, and (\textbf{d}) typical spectra acquired on each terrace. (\textbf{e}) Schematic showing the correspondence between the above observations and the \textit{ACAC} stacking pattern.}
\end{figure}

Figure S1(c) shows the spatially resolved tunnelling spectroscopy collected along the red-tinted rectangle shown in (a) \& (b). The relatively small gap is seen on the left-hand-side (Region I), and does not change when crossing the ordinary DW into Region II. The spectrum does change when crossing the \textit{buried} DW in to Region III however. This provides strong evidence explicitly linking the change of DOS to a change of the inter-layer stacking vector \textbf{T}, although the exact stacking vectors cannot be inferred only from this measurement. Representative spectra from each of the Regions I--III are shown in Fig. S1(d).

A possible scenario which could realise the above observations in shown in Figure S1(e). The top and second layer CDW patterns are depicted with blue and red star motifs, respectively. The DW in the top layer (blue dashes) changes the inter-layer stacking from one of the group of symmetrically equivalent $\textbf{T}_{C}$ vectors to another (from Region I to II). The DW in the lower layer (red dashes) changes the stacking from $\textbf{T}_{C}$ to $\textbf{T}_{A}$ (from Region II to III), which would induce the observed change in DOS from the small gap (unpaired SD layer) to the large gap (intact BL).

Figure S2 shows a partial replication of the results in Figs 2 \& 3 accompanying the main text, acquired at the cleaved surface of a different sample. This provides supporting evidence for the distinct DOS spectra of the two terminations of the BL stacking pattern. An array of white dots (shown in Fig. S2(a)), along with the maps of displacement field $\textit{u}_{r,\phi}(\textbf{x,y})$ (Fig. S2(b)), show an absence of a displacement in the in-plane projected CDW pattern at the step. Within the \textit{ACAC} stacking scheme, we interpret this as a step from the top of an intact BL (left-hand-side), down to the top of an unpaired layer (right-hand-side), formerly the bottom half of the same BL. The concurrent change in DOS from the larger to the smaller gap, shown in Fig S2(c) \& (d), is consistent with the findings described in the main text for Figs 2 \& 3. Figure S2(e) gives a schematic showing a possible surface configuration of SD layers for this case.

\subsection{II. Determination of SD cluster orientation shown in Figs. 3(c) \& (d).}

At the step-terrace formation discussed in the main text, the approximate in-plane displacements of Regions 1 \& 2 with respect to the underlying bilayer (BL) seen in Regions 3 \& 4, are observed. But in order to interpret these displacements in terms of stacking vectors \textbf{T}, the relative angular orientations of the atomic and CDW lattices must be determined.

The $\sqrt{13}\times\sqrt{13}$ \textit{R}13.9$^\circ$ superstructure of Star-of-David (SD) clusters is chiral, and can manifest in two configurations, shown in Fig. S3(a) \& (b), related to each other by a mirror symmetry operation. They can be described with angles $\pm13.9^\circ$ between the atomic and CDW lattice vectors.

\begin{figure*}
\centering
\includegraphics[scale=1]{FIGS3_v1.eps}
\caption{\label{fig:3} \textbf{Determination of relative orientations of atomic and CDW lattices.} (\textbf{a}) and (\textbf{b}) Depictions of two possible orientations of the underlying atomic lattice with respect to the CDW lattice. (\textbf{c}) Atomically resolved STM topography acquired in Region 1 of Figs. 2 \& 3 discussed in the main text. (\textbf{d}) The 2D Fourier transform of the topography. First order reciprocal lattice peaks for the CDW are circled in blue. The two possible sets of locations for the atomic Bragg peaks are circled in black and red. (\textbf{e}) Line cuts through the Fourier transform at $\pm13.9^\circ$ with respect to the axis of the first order CDW reciprocal lattice peaks. The inset depicts the likely S surface lattice positions as red dots, for comparison with panel (a).}
\end{figure*}

Atomically resolved topography acquired in the same domain as Region 1 of Figs. 2 \& 3 accompanying the main text, shown in Fig. S3(c), became obtainable after the STM tip was contaminated by unstable debris. In the corresponding 2D Fourier transform image (Fig. S3(d)), the first order peaks of the CDW reciprocal lattice are circled in blue. Higher order peaks of the extended CDW reciprocal lattice appear over the entire \textbf{q}-space field of view. The atomic Bragg peaks are expected to lie on top of one set of these higher order CDW peaks, at the positions marked by sets of black or red circles in Fig. S3(d). Taking linecuts at $\pm13.9^\circ$ rotation away from the direction of a first order CDW reciprocal lattice peak (Fig. S3(e)), we find that the intensity at the expected position of the atomic Bragg peak is higher at $+13.9^\circ$ than at $-13.9^\circ$. The corresponding real-space atomic lattice is depicted in the inset, with red dots. In the 1\textit{T}-TaS$_{2}$ structure the same relative orientation between the CDW and atomic lattices can be assumed for each of the three observed terraces. At this point the in-plane in-plane displacements discussed in the main text can be interpreted in terms of the atomic lattice sites within the SD clusters depicted in panel (a), and the stacking vectors \textbf{T} can be determined.

Strongly Correlated Electrons (cond-mat.str-el); Materials Science (cond-mat.mtrl-sci)