author: J. Di Francesco, S. Sadavoy, F. Motte, N. Schneider, M. Hennemann, T. Csengeri, S. Bontemps, Z. Balog, A. Zavagno, Ph. André, P. Saraceno, M. Griffin, A. Men'shchikov, A. Abergel, J.-P. Baluteau, J.-Ph. Bernard, P. Cox, L. Deharveng, P. Didelon, A.-M. di Giorgio, P. Hargrave, M. Huang, J. Kirk, S. Leeks, J. Z. Li, A. Marston, P. Martin, V. Minier, S. Molinari, G. Olofsson, P. Persi, S. Pezzuto, D. Russeil, M. Sauvage, B. Sibthorpe, L. Spinoglio, L. Testi, D. Teyssier, R. Vavrek, D. Ward-Thompson, G. White, C. Wilson, and A. Woodcraft
date: Received 2010 March 31 ; accepted 2010 May 7
title: Small-scale structure in the Rosette molecular cloud
  revealed by *Herschel*[^1]

We present a preliminary analysis of the small-scale structure found in new 70-520 $\mu$m continuum maps of the Rosette molecular cloud (RMC), obtained with the SPIRE and PACS instruments of the *Herschel* Space Observatory. We find 473 clumps within the RMC using a new structure identification algorithm, with sizes up to $\sim$`<!-- -->`{=html}1.0 pc in diameter. A comparison with recent *Spitzer* maps reveals that 371 clumps are "starless" (without an associated young stellar object), while 102 are "protostellar." Using the respective values of dust temperature, we determine the clumps have masses ($M_{C}$) over the range -0.75 $\leq$ log ($M_{C}$/M$_{\odot}$) $\leq$ 2.50. Linear fits to the high-mass tails of the resulting clump mass spectra (CMS) have slopes that are consistent with those found for high-mass clumps identified in CO emission by other groups.

# Introduction

Stars form within molecular clouds, after some fraction of cloud material is first condensed into smaller-scale structures. How this process unfolds is not well understood, however. For example, what conditions within clouds drive the formation of small-scale structures that in turn produce stars of any given mass, e.g., high-mass stars. Detailed observations of the small-scale structure within molecular clouds should provide valuable insight into how stars of various masses form out of dense material.

The need for sensitive probes of small-scale structure in molecular clouds is the impetus behind the "*Herschel* OB Young Stellar objects" (HOBYS) Key Project (see Motte et al. 2010). The HOBYS team is currently using the ESA *Herschel* Space Observatory (Pilbratt et al. 2010) to obtain wide-field maps of $\sim$`<!-- -->`{=html}15 molecular clouds which are forming high-mass stars within 3 kpc of the Sun. These clouds are being observed with both the *Herschel* SPIRE (Griffin et al. 2010) and PACS (Poglitsch et al. 2010) instruments, ultimately to obtain wide maps of these clouds at 70-520 $\mu$m at diffraction-limited resolutions of $\sim$`<!-- -->`{=html}18$^{\prime\prime}$ $\times$ ($\lambda$/250 $\mu$m). Such data sample thermal emission from cold dust grains mixed with molecular gas within the cloud. Since continuum emission at submillimetre wavelengths is typically optically thin, detections of continuum emission from dust at high resolution can trace the organization of mass in clouds on small scales. By sensitively sampling emission across the peaks of the spectral energy distributions (SEDs) of this dust, SPIRE and PACS data can constrain both the column density and temperature of the dust and provide unparalleled censuses of the small-scale structure in molecular clouds.

In this paper, we describe the small-scale structure detected with SPIRE and PACS in the Rosette molecular cloud (RMC), the first of the HOBYS target list that was observed as part of the early *Herschel* "science demonstration phase" (SDP) campaign. The RMC is a rich location of star formation at a distance of 1.6 kpc (Johnson 1962; Pérez et al. 1987; Park & Sung 2002). The cloud is southeast of NGC 2244, the Rosette Nebula, and the expanding HII region associated with NGC 2244 has begun to interact with it. The structure of material within the RMC has been previously examined by Williams, Blitz & Stark (1995), Schneider et al. (1998) and Dent et al. (2009) using various CO observations. Alternatively, the populations of young stellar objects (YSOs) within the RMC have been studied by Phelps & Lada (1998), Li & Smith (2005), Poulton et al. (2008), and Román-Zúñiga et al. (2008) using infrared observations. Notably, these studies identified 10 embedded clusters within the RMC, i.e., PL01-07 and REFL08-10. A recent review of observations of the RMC has been made by Román-Zúñiga & Lada (2008; RZL).

The SPIRE and PACS observations described in this paper provide a preliminary look into how *Herschel* data trace small-scale structures in the RMC, and the relationships between these structures and YSOs. Analyses of other aspects within the SDP data of the RMC can be found in the accompanying papers by Motte et al. 2010, Henneman et al. 2010, and Schneider et al. 2010.)

# Observations

A complete description of the observations described here can be found in the accompanying paper by Motte et al. (2010). In short, an extinction map of the RMC was first constructed for planning purposes using 2MASS data and the AvMAP algorithm. AvMAP derives extinction maps with typical resolutions of $\sim$`<!-- -->`{=html}2$^{\prime}$ using colour excesses in a manner akin to the NICER algorithm of Lombardi & Alves (2001). AvMAP, however, also uses a stellar model to predict foreground source densities; see Schneider et al. (2010a). Using the AvMAP data as a guide, an area of the RMC $\approx$ 1 deg$^{2}$ in extent encompassing areas southeast of NGC 2244 with $A_{V}$ $>$ 5 was mapped in common with SPIRE and PACS in parallel scan mode on 2009 October 20. Maps were obtained simultaneously at 70 $\mu$m and 170 $\mu$m with PACS and at 250 $\mu$m, 350 $\mu$m, and 520 $\mu$m with SPIRE.

The SPIRE and PACS data were reduced separately using common routines of the *Herschel* Interactive Processing Environment[^2] (HIPE; see Ott 2010). To facilitate the identification and analysis of the cold, small-scale structure in the RMC, i.e., dense clumps, the PACS maps were Gaussian smoothed to the 18$^{\prime\prime}$ FWHM resolution of the 250 $\mu$m SPIRE map. A full analysis of the unsmoothed PACS maps will be provided in the near future by Sadavoy et al. (2010b, in preparation).

# Results

Figure 1 (available electronically) shows the 350 $\mu$m map of the RMC, revealing a rich amount of structure in continuum emission. (See Motte et al. for maps at all five wavelengths.) Several bright knots are seen embedded in large filaments, and each can be associated with known infrared clusters (see Figure 1, also RZL). Ubiquitous, faint diffuse emission is also seen, and in many cases this appears filamentary as well. Of special note in Figure 1 is the filamentary arc seen in the northwestern corner, indicating the "ridge" of the RMC (+ "Core D"), i.e., where the adjacent NGC 2244 HII region (the Rosette Nebula) is interacting with the molecular cloud.

From the (smoothed) PACS and SPIRE maps, 473 continuum objects were identified using the *getsources* algorithm; an accompanying paper in this volume by Men'shchikov et al. (2010) describes *getsources* in more detail. In brief, *getsources* first decomposes maps into images over a range of scales and objects are identified as discrete emission regions in each image. Objects are then defined using information on many scales and at many wavelengths. For example, positions and extents (i.e., "footprints") are determined from the first and second moments of intensities in images where all scales are recombined. Integrated flux densities are calculated by summing the background-subtracted emission under the footprint of an object. Sizes are well-approximated by 2D Gaussian fits.

For simplicity, we choose to identify all objects found by *getsources* as "clumps"[^3] Figure 1 shows the locations of all 473 clumps identified in the RMC overlaid onto the 350 $\mu$m map. Most are found within the brighter filaments of the RMC, with neighbouring clumps tracing the linear structures of underlying filaments. Several "isolated" clumps (i.e., those without close neighbours), however, are also seen associated with fainter emission. Relatively few clumps are seen inside the northwestern "ridge" and coincident with the NGC 2244 HII region.

Figure 1 also shows which clumps in the RMC contain YSOs ("protostellar clumps") and which do not ("starless clumps"). These determinations were made using *Spitzer* IRAC and MIPS data of the region (see Poulton et al. 2008; Balog et al. 2010, in preparation). The criteria used to determine which clumps contain YSOs are based on those developed recently by Sadavoy et al. (2010a) for small-scale structure detected at 850 $\mu$m in the less-distant Orion, Taurus, Perseus and Ophiuchus molecular clouds. For this comparison, potential YSOs had to be detected in all four IRAC bands or at 24 $\mu$m. Their colours and brightness were checked against colour criteria to distinguish them as either YSOs or other objects (e.g., background galaxies or stars). In addition, their positions relative to clumps were checked for coincidence within areas defined by the 2D Gaussian fits to the clump extents. We do not distinguish here between protostellar Classes of YSOs in clumps; see Henneman et al. for further discussion.

In total, we find 371 starless clumps and 102 protostellar clumps in the RMC. Most protostellar clumps contain one YSO, but $\sim$`<!-- -->`{=html}35% contain more than one. Figure 1 suggests that brighter clumps tend to be protostellar, while fainter clumps tend to be starless. (Exceptions to both cases are seen, however.) As expected, the protostellar clumps are associated with the embedded clusters known in the cloud, i.e., PL01-07 and REFL08-10. Many starless clumps appear to be associated with the RMC ridge and the area between it and the RMC "central core" (i.e., the concentration of bright protostellar clumps associated with the embedded clusters PL04-06 and REFL08 in the centre of Figure 1). Relatively few ($\sim$`<!-- -->`{=html}15) clumps are seen behind the ridge itself, and they all appear to be starless. Excluding areas inside and outside the ridge, the starless and protostellar clumps appear well-mixed, suggesting no "age gradient" across the RMC. (Schneider et al. 2010b suggest, however, that the most massive cores in the RMC decrease in age with increased distance from NGC 2244.)

Figure 2 shows distributions of the flux and size of clumps in the RMC, as a whole and when separated into starless or protostellar clumps. Figure 2a shows histograms of 350 $\mu$m integrated fluxes, $S_{350}$, of the clumps, demonstrating how the brightest are indeed protostellar (see Figure 1). Starless and protostellar clumps with $S_{350}$ $\approx$ 5 Jy are similar in number, but clumps fainter than 5 Jy are predominantly starless. For reference (see §4 below), a clump of 5 Jy flux at 350 $\mu$m in the RMC has a mass of $\sim$`<!-- -->`{=html}7 M$_{\odot}$, assuming $T_{dust}$ = 20 K. Figure 2b shows size histograms of each population where size is the geometric mean of the observed major and minor axes of the 2D Gaussian fit to each clump extent. (These sizes are not deconvolved from the finite beam.) Clumps of each type have sizes that extend from the resolution limit up to $\sim$`<!-- -->`{=html}0.65 pc. The few larger clumps (i.e., 0.7-1.0 pc) are all protostellar, however.

Figure 3 shows distributions of the local extinction and temperatures of clumps within the cloud, again as a whole and when separated into starless or protostellar clumps. Figure 3a shows the $A_{V}$ values associated with the clumps, as determined from the low-resolution (i.e., $\sim$`<!-- -->`{=html}2$^{\prime}$) extinction maps. For a given clump, its associated low-resolution $A_{V}$ can be used as a proxy for its depth within its parent cloud (see Johnstone, Di Francesco & Kirk 2004). Clumps of each type are found across the same range of extinction ($A_{V}$ $\approx$ 3-20). No clumps of either type are found at $A_{V} <$ 2, however. Starless clumps are found in greater numbers at lower extinctions ($A_{V} <$ 12), but both types are found in similar numbers at higher extinctions ($A_{V} \geq$ 12). The number of starless clumps peaks at $A_{V}$ $\approx$ 6, while that of protostellar clumps peaks at $A_{V}$ $\approx$ 13-14. Though starless and protostellar clumps are found over the same range of $A_{V}$, the difference in peak $A_{V}$ may be due in part to protostellar clumps having higher masses on average than starless clumps (see §4 below). Assuming similar sizes, the resulting relatively higher column densities could mean relatively higher local $A_{V}$ values, though the difference in the AvMAP and *Herschel* resolutions is large.

Figure 3b shows the distributions of the clump temperatures ($T_{dust}$) found across the RMC using all five *Herschel* wavelengths and a method described by Schneider et al. (2010b) Temperatures could be found for only 402 of the 473 identified clumps (i.e., 323 starless and 79 protostellar clumps), and we limit further discussion only to those clumps. Dust temperatures for these clumps range between 10 K and 30 K. Interestingly, starless clumps appear to be generally warmer than protostellar clumps, a surprise as the latter could be expected to be warmer given internal heating by YSOs. The difference here, however, stems from the large population of starless clumps located to the northwest of the RMC (see Figure 1). At this location, the dust is quite warm owing to its proximity to the O stars in NGC 2244.

# Clump mass spectra

The mass for each clump, $M_{C}$, was determined using its observed 350 $\mu$m flux and the relation $M_{C}$ $=$ $S_{\nu}d^{2}$/$\kappa_{\nu}B_{\nu}$($T_{dust}$) = 0.073($S_{\nu}/\rm{Jy}$)($d/\rm{kpc}$)$^{2}$($\kappa_{\nu}/0.07$ cm$^{2}$g$^{-1}$)$^{-1}$\[exp(41.1/$T_{dust}$)-1\] M$_{\odot}$. To account for variations in temperature across the RMC, we assume for each clump its respective $T_{dust}$ value when determining its mass, rather than averages of the distributions shown in Figure 3b. We also use $\kappa_{350}$ = 0.07 cm$^{2}$ g$^{-1}$for a dust opacity, consistent with the dust prescription of Ossenkopf & Henning (1994), and a distance of 1.6 kpc.

Figure 4 shows the distributions of mass for the total, starless, and protostellar clump populations, i.e., their clump mass spectra (CMS). Masses for all clumps range across $\sim$`<!-- -->`{=html}3.5 orders of magnitude, i.e., -0.75 $\leq$ log ($M_{C}$/M$_{\odot}$) $\leq$ 2.50. Evidence for differences in mass between the populations is seen, however, with starless clumps having on average lower masses than protostellar clumps. Moreover, the starless population peaks at log ($M_{C}$/M$_{\odot}$) = 0.25, while the protostellar population peaks at log ($M_{C}$/M$_{\odot}$) = 1.0. This difference is likely because the starless clumps in the RMC are both warmer and have lower submillimetre fluxes on average than the protostellar clumps (see Figures 2 and 3). Note also that there are fewer higher mass starless clumps than protostellar clumps, e.g., over the high-mass range of log ($M_{C}$/M$_{\odot}$) = 1.75-2.50, only two starless clumps are found but 11 protostellar clumps.

Table 1 lists the linear slopes ($\alpha$) fit to the high-mass tails of the respective CMS, using the least-squares method and assuming Poissonian errors in each mass bin. Table 1 also lists the mass ranges used to compute the slopes. In general, slopes were fit between the bin containing the peak of the respective CMS and the bin of the highest mass in the sample. The $\alpha$ = -0.65 $\pm$ 0.01 for all high-mass clumps is slightly more than the value of -0.8 found by Dent et al. (2009) from CO 3-2 observations over a 5 deg$^{2}$ map of the entire Rosette complex. The value found here, however, is quite similar to the -0.6 slope found by Schneider et al. (1998) for higher clump masses using CO (2-1) maps of the Rosette. Both these groups used different clump identification algorithms than used here, i.e., *clumpfind* by Dent et al. and *gaussclumps* by Schneider et al. (1998). Moreover, we have determined masses from clumps we identified here, not over the extents of the clumps found by either group using CO.

```latex
\centering{
\begin{tabular}{ccc}
\\
\hline
\hline
Clump Type & $\alpha$   &  Mass Range \\
\hline
All        & -0.65 $\pm$ 0.01 & 0.25 $\leq$ log~($M_{C}$/M$_{\odot}$) $\leq$ 2.50 \\
Starless   & -0.82 $\pm$ 0.02 & 0.25 $\leq$ log~($M_{C}$/M$_{\odot}$) $\leq$ 2.25 \\
Protostellar & -0.80 $\pm$ 0.02 & 1.00 $\leq$ log~($M_{C}$/M$_{\odot}$) $\geq$ 2.50 \\
\hline
\end{tabular}
}
```

After dividing the clumps by type, our starless and protostellar clumps have similar slopes to their high-mass distributions, i.e., $\alpha$ = -0.8, but over different mass ranges. Indeed, these slopes are more similar to that found by Dent et al. The high-mass tails of the all-clump CMS has a shallower slope than do just starless or protostellar clumps. This difference stems from the difference in the mass ranges used to measure slopes. Also, the numbers of protostellar clumps peak at a higher mass than do those of starless clumps, and they dominate the all-clump CMS at even higher masses.

With single fits, all three CMS are shallower at higher masses than seen in the high-mass tails of mass spectra of small-scale structures within closer clouds. Indeed, the slopes of such tails are seen to be consistent with the -1.35 slope of the Salpeter IMF (e.g., in Aquila; see André et al. 2010). For comparison with the RMC, a Salpeter slope is also shown in Figure 4. This difference may come simply from resolution; namely, in closer clouds, sites of individual star formation, ("cores") can be resolved. For example, typical core separations in Ophiuchus at 125 pc were found to be $\sim$`<!-- -->`{=html}0.03 pc by Motte, André & Neri (1998). Though our minimum linear resolution (at 250 $\mu$m) is $\sim$`<!-- -->`{=html}0.12 pc, or about the expected size of isolated cores, we may not be able to resolve individual cores in the RMC owing to crowding on scales $<$ 0.12 pc. Indeed, such crowding may cause the resulting mass spectra to be artificially shallow. Higher resolution data will be critical for testing this possibility.

Without appropriate spectral line data, we cannot estimate the states of gravitaitonal stability for the clumps in the RMC, and determine which may be bound or which may be transient structures. For example, we cannot presently guess which starless clumps will likely engage in future star formation. It is clear from Figure 1, however, that the *Herschel* data have revealed that star formation is ongoing in the RMC and that we are seeing only a snapshot of a very dynamical situation.

# Summary

We have conducted a preliminary analysis of small-scale emission within the Rosette molecular cloud, as observed by the *Herschel* Space Observatory. From SPIRE and smoothed PACS images, we identified 473 clumps using *getsources*, a new algorithm that identifies structure over multiple wavelengths and scales. Using *Spitzer* data, we classified these into 371 starless and 102 protostellar clumps. The clumps have dust temperatures of 10-30 K, with starless clumps surprisingly warmer on average than the protostellar clumps, owing to the proximity of many to NGC 2244. Masses were determined for each clump, revealing a range of -0.75 $\leq$ log ($M_{C}$/M$_{\odot}$) $\leq$ 2.50. Linear least-squares fits to the high-mass tails of the CMS reveal slopes that are consistent with slopes fit to mass spectra of clumps identified through CO in the RMC. Such slopes are shallower than those seen in the high-mass tails of the mass spectra of small-scale continuum structures ("cores") in closer clouds, possibly from crowding of such objects on scales smaller than probed here.

[^1]: *Herschel* is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA.

[^2]: HIPE is a joint development by the Herschel Science Ground Segment Consortium, consisting of ESA, the NASA Herschel Science Center, and the HIFI, PACS and SPIRE consortia.

[^3]: The term "clump" has been used previously to describe sites of grouped or clustered star formation, as opposed to sites of single- (or multiple-) star formation, i.e., "cores"; see Williams, Blitz & McKee (2000). We use "clump" here to describe all discrete condensations of material of size $<$ 1 pc within the larger cloud.


author: Bryan Ford [^1]; Rainer Böhme [^2]
bibliography: rationality.bib
title: Rationality is Self-Defeating in Permissionless Systems [^3]

# Introduction

Many blockchain and cryptocurrency fans seem to prefer building and analyzing decentralized systems in a rational or "greedy behavior" failure model, rather than a Byzantine or "arbitrary behavior" failure model. Many of the same blockchain and cryptocurrency fans also like open, permissionless systems like Bitcoin and Ethereum, which anyone can join and participate in using weak identities such as anonymous cryptography key pairs.

What most of these heavily-overlapping sets of fans do not seem to realize, however, is that rationality assumptions are self-defeating in open permissionless systems with weak identities. A fairly simple metacircular argument—a kind of "Gödel's incompleteness theorem for rationality"—shows that for any system $S$ that makes *any* behavioral assumption, including but not limited to a rationality assumption, a rational attacker both exists and *has an incentive* to defeat that behavioral assumption, thereby violating that assumption and exhibiting Byzantine behavior from the perspective of the system.

As a quick summary of the argument we will expand below, suppose a permissionless system like Bitcoin is secure against rational attacks, but has some weakness against irrational Byzantine attacks in which the attacker would lose money. Because the system is open, permissionless, and exists within a larger ecosystem, a rational attacker can find ways to "bet against" Bitcoin's security in *other* financially-connected systems (e. g., Ethereum), making a profit *outside of* Bitcoin on this attack against Bitcoin. An attack that appears irrational in the context of Bitcoin may be perfectly rational in the context of the larger ecosystem.

For this reason, an open permissionless system designed to be secure only against rational adversaries is actually just *insecure*, unless it remains secure even when the "rational" participants become fully Byzantine. Given this, one might as well have designed the permissionless system in a Byzantine model in the first place. The rationality assumption offers no actual benefit, but merely can make an insecure system appear secure under flawed analysis.

This manuscript is based partly on ideas in the second author's lecture  at the BDLT Summer School in Vienna . While formalizing the argument would require some effort, we thought it would be worth at least sketching the argument intuitively for the public record.

# Threat Modeling: Honest, Byzantine, and Rational Participants

In designing or analyzing the security of any decentralized system, we must define the system's *threat model*, and in particular our assumptions about the behaviors of the participants in the system. An *honest*, *correct*, or *altruistic* participant is one that we assume to follow the system's protocol rules as specified, hence representing a "well-behaved" participant exhibiting no adversarial behavior.

A *Byzantine* participant, named after the Byzantine Generals Problem , is one we make *no* assumptions about. A Byzantine participant can behave in *arbitrary* fashion, without restriction, and hence by definition represents the strongest possible adversary.

We would like to build systems that could withstand *all* participants being Byzantine, but this appears fundamentally impossible. We therefore in practice have to make threshold security assumptions, such as that over two-thirds of the participants in classic Byzantine consensus protocols are honest, or that the participants controlling over half the hash power in Bitcoin are well-behaved.

Even with threshold assumptions, however, building systems that resist Byzantine behavior is extremely difficult, and the resulting systems are often much more complex and inefficient than systems tolerating weaker adversaries. We may therefore be tempted to improve a design's simplicity or efficiency by making stronger assumptions about the behavior of adversarial participants, effectively weakening the assumed adversary.

One such popular assumption, especially in economic circles, is *rationality*. In essence, we assume that rational participants may deviate from the rules in arbitrary ways but *only when doing so is in their economic self-interest*, improving their expected rewards—usually but not always financial—in comparison with following the rules honestly (cf. Fig. ).

By assuming that adversarial participants are rational rather than Byzantine, we need not secure the system against *all* possible participant behaviors, such as against participants who pay money with no reward merely to sow chaos and destruction. Instead, we merely need to prove that the system is *incentive compatible*, for example by showing that its rules represent a Nash equilibrium, in which deviations from the equilibrium will not give participants a greater financial reward.

Besides simplicity and efficiency, another appeal of rationality assumptions is the promise of *strengthening* the system's security by lowering the threshold of participants we assume to be fully honest. To circumvent the classical Byzantine consensus requirement that fewer than one third of participants may be faulty, for example, we might hope to tolerate closer to 50 %, or even 100 %, of participants being "adversarial" if we assume they are rational and not Byzantine. Work on the Byzantine-Altruistic-Rational (BAR) model  and $(k,t)$-robustness  exemplifies this goal, which sometimes appears achievable in closed systems with strong identities. But a direct implication of our metacircular argument is that an *open* system cannot generally be secure if all participants are either Byzantine or rational.

# Assumptions Underlying the Argument

The metacircular argument makes three main assumptions.

#### Assumption 1

First, the system $S$ under consideration is open and permissionless, allowing anyone to join and participate in the system using only weak, anonymous identities such as bare cryptographic key pairs. Identities in $S$ need not even be costless  provided their price is modest: the argument still works even if $S$ imposes membership fees or requires new wallet keys to be "mined," for example. Proof-of-Work cryptocurrencies such as Bitcoin and Ethereum, Proof-of-Stake systems such as Algorand  and Ouroboros , and most other permissionless systems seem to satisfy this openness property. Because participation is open to anyone globally and can be anonymous, we cannot reasonably expect police or governments to protect $S$ from attack: even if they wanted to and considered it their job, they would not be able to find or discipline a smart rational attacker who might be attacking from anywhere around the globe, especially from a country with weak international agreements and extradition rules. Thus, $S$ must "stand on its own," by successfully either withstanding or disincentivizing attacks coming from anywhere. (And it will turn out that merely disincentivizing such attacks is impossible.)

#### Assumption 2

Second, the system $S$ does not control a majority of total economic power or value in the world: i. e., it is not totally economically dominant from a global perspective. Instead, there may be (and probably are) actors outside of $S$ who, if rationally incentivized to do so, can at least temporarily muster an amount of economic power outside of $S$ comparable to or greater than the economic value within or controlled by $S$. In other words, we assume that $S$ is not the "biggest fish in the ocean." Given that there can be at most one globally dominant economic system at a time, it seems neither useful nor advisable to design systems that are secure only when they are the biggest fish in the ocean, because almost always they are not.

#### Assumption 3

Third, the system $S$ actually *leverages* in some fashion the behavioral assumption(s) it makes on participants, such as a rationality assumption. That is, we assume there exist one or more (arbitrary) behavioral strategies that $S$ assumes some participants *will not* follow, such as economically-losing behaviors in the case of rationality. Further, we assume there exists such an assumption-violating strategy that will cause $S$ to malfunction or otherwise deviate observably from its correct operation. In fact, we need not assume that this deviant behavior will *always* succeed in breaking $S$, but only that it will non-negligibly *raise the probability* of $S$ failing. If this were not the case, and $S$ in fact operates correctly, securely, and indistinguishably from its ideal even if participants do violate their behavioral assumptions, then $S$ is actually Byzantine secure after all. In that case, $S$ is not actually benefiting from its assumptions about participant behavior, which are redundant and thus may be simply discarded.

# The Metacircular Argument: Rational Attacks on Rationality

Suppose permissionless system $S$ is launched, and operates smoothly for some time, with all participants conforming to $S$'s assumptions about them. Because $S$ is permissionless (assumption 1) and exists in a larger open world (assumption 2), new rational participants may arrive at any time, attracted by $S$'s success and presumably growing economic value provided there is an opportunity to profit from doing so.

Consider a particular newly-arriving participant $P$. $P$ could of course play by the rules $S$ assumes of $P$, in which case the greatest immediate economic benefit $P$ could derive from participating in $S$ is some fraction of the total economic value currently embodied in $S$ (e. g., its market cap). For most realistic permissionless systems embodying strong founders' or early-adopters' rewards, if $P$ is not one of the original founders of $S$ but arrives substantially after launch, then $P$'s near-term payoff prospectives from joining $S$ is likely bounded to a fairly *small* fraction of $S$'s total value. But what if there were another strategy $P$ could take, for perfectly *rational* and economically-motivated reasons, by which $P$ could in relatively short order acquire a *large* fraction of $S$'s total value?

Because $S$ is permissionless and operating in a larger open world, $P$ is not confined to operating exclusively within the boundaries of $S$. $P$ can also make use of facilities external to $S$. By assumption 2, $P$ may in particular have access to, or be able to borrow temporarily, financial resources comparable to or larger than the total value of $S$.

Suppose the facilities external to $S$ include another Ethereum-like cryptocurrency $S'$, which includes a smart contract facility with which decentralized exchanges, futures markets, and the like may be implemented. (This is not really a separate assumption because even if $S'$ did not already exist, $P$ could create and launch it, given sufficient economic resources under assumption 2.) Further, suppose that someone (perhaps $P$) has created on external system $S'$ a decentralized exchange, futures market, or any other mechanism by which tokens representing shares of the value of $S$ may be traded or speculated upon in the context of $S'$: e. g., a series of tradeable Ethereum tokens pegged to $S$'s cryptocurrency or stake units (see Fig. ).

Now suppose participant $P$ finds some behavioral strategy that system $S$ depends on participants *not* exhibiting, and that will observably break $S$—or even that just *might* break $S$ with significant non-negligible probability. Assumption 3 above guarantees the existence of such a behavioral strategy, unless $S$'s rationality assumptions were in fact redundant and worthless. $P$ must merely be clever enough to find and implement such a strategy. It is possible this strategy might first require $P$ to pretend to be one or more well-behaved participants of $S$ for a while, to build up the necessary reputation or otherwise get correctly positioned in $S$'s state space; a bit of patience and persistence on $P$'s part will satisfy this requirement. $P$ may also have to "buy into" $S$ enough to surmount any entry costs or stake thresholds $S$ might impose; the external funds $P$ can invoke or borrow by assumption 2 can satisfy this requirement, and are bounded by the total value of $S$. In general, $S$'s openness by assumption 1 and the existence of a correctness-violating strategy by assumption 3 ensures that there exists some course of action and supply of external resources by which $P$ can position itself to violate $S$'s behavioral assumption.

In addition to infiltrating and positioning itself within $S$, $P$ also invokes or borrows enough external funds and uses them to short-sell (bet against) shares of $S$'s value massively in the context of the external system $S'$, which (unlike $S$) $P$ trusts will remain operational and hold its value independently of $S$. Provided $P$ reaches this short-selling position gradually and carefully enough to avoid revealing its strategy early, the funds $P$ must invoke or borrow for this purpose must be bounded by some fraction of the total economic value of $S$. And provided there are at least some participants and/or observers of $S$ who believe that $S$ is secure and will remain operating correctly, and are willing to bet to that effect on $S'$, $P$ will eventually be able to build its short position.

Finally, once $P$ is positioned correctly within both $S$ and $S'$, $P$ then launches its assumption-violating behavior in $S$ that will observably cause $S$ to fail as per assumption 2. This might manifest as a denial-of-service attack, a correctness attack, or in any other fashion. The only requirement is that $P$'s behavior creates an *observable* failure, which a nontrivial number of the existing participants in $S$ believed would not happen because they believed in $S$ and its threat model. The fact that $S$ is now observed to be broken, and its basic design assumptions manifestly violated, causes the shares of $S$'s value to drop precipitously on external market $S'$, on which $P$ takes a handsome profit. Perhaps $S$ recovers and continues, or perhaps it fails entirely—but either way, $P$ has essentially transferred a significant fraction of system $S$'s economic value from system $S$ itself to $P$'s own short-sold position on external market $S'$. And to do so, $P$ needed only to find a way—any way—to *surprise* all those who believed $S$ was secure and that its threat model accurately modeled $S$'s real-world participants.

Even if $P$'s assumption-violating behavioral strategy does not break $S$ with perfect reliability, but only with some probability, $P$ can still create an *expectation* of positive profit from its attack by hedging its bets appropriately on $S'$. $P$ does not need a perfect attack, but merely needs to possess the *correct* knowledge that $S$'s failure probability is much higher than the other participants in $S$ believe it to be—because only $P$ knows that (and precisely when) it will violate $S$'s design assumptions to create that higher failure probability. Furthermore, even if $P$'s attack fails, and the vulnerability it exploits is quickly detected and patched, $P$ may still profit marginally from the market's adjustment to a realization that $S$'s failure probability was (even temporarily) higher than most of $S$'s participants thought it was.

Within the context of system $S$, $P$'s behavior manifests as Byzantine behavior, specifically violating the assumptions $S$'s designers thought participants would not exhibit and thus excluded from $S$'s threat model. Considered in the larger context of the external world in which $S$ is embedded, however, including the external trading system $S'$, $P$'s behavior is perfectly rational and economically-motivated. Thus, the very rationality of $P$ in the larger open world is precisely what motivates $P$ to break, and profit from, $S$'s ill-considered assumption that its participants would behave rationally.

# Implications for Practical Systems

This type of financial attack is by no means entirely theoretical or limited to fully-digital systems such as cryptocurrencies. In our scenario, $P$ is essentially playing a game closely-analogous to the investors in credit default swaps who both contributed to, and profited handsomely from, the 2007–2008 financial crisis. as covered more recently in the film "The Big Short."

In the cryptocurrency space, some real-world attacks we are seeing—such as increasingly-common 51 % attacks —might be viewed as special cases of this metacircular attack on rationality. It is often claimed that large proof-of-work miners (or proof-of-stake holders) will not attempt 51 % attacks because doing so would undermine the value of the cryptocurrency in which they by definition hold a large stake, and hence would be "irrational." But this argument falls apart if the attack allows the large stakeholder to reap rewards outside the attacked system, e. g., by defrauding exchanges or selling $S$ short in other systems.

Externally-motivated attacks on cryptocurrencies have been predicted before in the form of virtual protest or "Occupy Bitcoin" attacks , Goldfinger attacks  puzzle transaction attacks , merged mining attacks , hostile blockchain takeovers , and out-of-band variants of pay-to-win attacks . All these attacks are specific instances of our argument. They have been presented in the literature as open yet solvable challenges. We are not aware, however, of any prior attempt to summarize the lessons learned and formulate a general impossibility statement.

For most practical systems, we do not even know if they are incentive compatible in the absence of an external system $S'$—i. e., where assumption 2 is violated—and probably they are not. Almost all game-theoretic treatments of (parts of) the Bitcoin protocol deliver negative results. Many attacks against specific cryptocurrency system designs are known to be profitable in expectation, such as transaction withholding , empty block mining , selfish mining , block withholding , stubborn mining , fork after withholding , and whale attacks  It is likely thanks only to frictions such as risk aversion and other costs that we rarely observe such attacks in large deployed systems. Many specific attacks do not even depend on assumption 1, underlining the fact that rationality is not a silver bullet even where this metacircular argument does not apply. Where it does apply, it is more general and effectively *guarantees* the existence of attacks against *all* open systems that assume participants are rational.

Another related observation is that financial markets on derivatives of a system $S$ mature in the external world (e. g., $S'$) as $S$ grows and becomes more relevant. So in some sense, systems built on the rationality assumption are temporarily more secure only until they become fat enough targets to be eaten by their own success. We can see this effect, for example, in the growing and increasingly liquid market for hash power, which effectively thwarts Nakamoto's  (or Dwork's ) rule of thumb that the ratio of processors to individuals varies in a small band. Such dynamics happen in the real world, too. But there they have traditionally taken centuries or decades while in cryptocurrency space everything happens in time-lapse.

# Limitations of the Argument

This argument is of course currently only a rough and informal sketch. An enterprising student might wish to try formalizing it, or maybe someone has already done so but we are unaware of it.

The metacircular argument certainly does not apply to all cryptocurrencies or decentralized systems. In a permissioned system, for example, in which a closed group of participants are strongly-identified and subject to legal and contractual agreements with each other, one can hope that the threat of lawsuits for arbitrarily-large damages will keep rational participants incentivized to behave correctly. Similarly, in a national cryptocurrency, which might be relatively open but only to citizens of a given country, and which require verified identities with which the police can expect to track down and jail misbehaving participants, this metacircular argument does not necessarily apply.

Apart from police enforcement, rationality assumptions may be weakened in other ways to circumvent the metacircular argument. For example, an open system might be designed according to a "weak rationality" assumption that users need incentives to join the system in the first place (e. g., mining rewards in Bitcoin), but that after having become stakeholders, most will then behave honestly. In this case, rational incentives serve only as a tool for system growth, but become irrelevant and equivalent to a strong honesty assumption in terms of the internal security of the system itself.

# Conclusion: Irrationality Can Be Rational

What many in the cryptocurrency community seem to want is a system that is both permissionless and tolerant of strongly-rational behavior—either beyond the thresholds a similar a Byzantine system would tolerate (such as a rational majority), or by deriving some simplicity or efficiency benefit from assuming rationality. But in an open world in which the permissionless system is not the only game in town, a potential *perfectly rational* attacker can always exist, or appear at any time, whose entirely rational behavior is precisely to profit from bringing the system down by violating its assumptions on participant behavior.

So if you think you have designed a permissionless decentralized system that is cleverly secured based on rationality assumptions, you haven't. You have merely obfuscated the rational attacker's motive and opportunity to profit outside your system from breaking your rationality assumptions. The only practical way to eliminate this threat appears to be either to close the system and require strong identities and police protection, or else secure the system against arbitrary Byzantine behavior, thereby rendering rationality assumptions redundant and useless for security (cf. Fig. ).

### Acknowledgements

We wish to thank Jeff Allen, Ittay Eyal, Damir Filipovic, Patrik Keller, Alexander Lipton, Andrew Miller, and Haoqian Zhang for helpful feedback on earlier drafts of this work.

[^1]: École Polytechnique Fédérale de Lausanne, Switzerland, <a href="firstname.lastname@example.com" class="uri">firstname.lastname@example.com</a>

[^2]: Universität Innsbruck, Austria, <a href="email@example.com" class="uri">firstname.lastname@example.com</a>

[^3]: First posted on September 23, 2019 at <https://bford.info/2019/09/23/rational/>


abstract: We use CANDELS imaging, 3D-HST spectroscopy, and Chandra X-ray data to investigate if active galactic nuclei (AGNs) are preferentially fueled by violent disk instabilities funneling gas into galaxy centers at $1.3<z<2.4$. We select galaxies undergoing gravitational instabilities using the number of clumps and degree of patchiness as proxies. The CANDELS visual classification system is used to identify 44 clumpy disk galaxies, along with mass-matched comparison samples of smooth and intermediate morphology galaxies. We note that, despite being being mass-matched and having similar star formation rates, the smoother galaxies tend to be smaller disks with more prominent bulges compared to the clumpy galaxies. The lack of smooth extended disks is probably a general feature of the $z \sim
    2$ galaxy population, and means we cannot directly compare with the clumpy and smooth extended disks observed at lower redshift. We find that $z \sim 2$ clumpy galaxies have slightly enhanced AGN fractions selected by integrated line ratios (in the mass-excitation method), but the spatially resolved line ratios indicate this is likely due to extended phenomena rather than nuclear AGNs. Meanwhile the X-ray data show that clumpy, smooth, and intermediate galaxies have nearly indistinguishable AGN fractions derived from both individual detections and stacked non-detections. The data demonstrate that AGN fueling modes at $z \sim 1.85$ - whether violent disk instabilities or secular processes - are as efficient in smooth galaxies as they are in clumpy galaxies.
author: Jonathan R. Trump, Guillermo Barro, Stéphanie Juneau, Benjamin J. Weiner, Bin Luo, Gabriel B. Brammer, Eric F. Bell, W. Niel Brandt, Avishai Dekel, Yicheng Guo, Philip F. Hopkins, David C. Koo, Dale D. Kocevski, Daniel H. McIntosh, Ivelina Momcheva, S. M. Faber, Henry C. Ferguson, Norman A. Grogin, Jeyhan Kartaltepe, Anton M. Koekemoer, Jennifer Lotz, Michael Maseda, Mark Mozena, Kirpal Nandra, David J. Rosario, and Gregory R. Zeimann
title: No More Active Galactic Nuclei in Clumpy Disks Than in Smooth Galaxies at $z \sim 2$ in CANDELS / 3D-HST

# Introduction

The observed correlations between the mass of a supermassive black hole (SMBH) and the properties of its host galaxy bulge imply that SMBH and bulge growth may be linked. To maintain the SMBH-bulge relations over a galaxy's long evolutionary history, rapid SMBH accretion in active galactic nucleus (AGN) phases probably coincide with periods of high star formation rate in the host galaxy. Indeed, rapidly accreting AGNs are observed to be predominantly located in the rapidly star-forming galaxies, from the local universe to higher redshifts . But the detailed physical processes behind the coupled growth of SMBHs and galaxies remain mysterious.

A key difficulty for any AGN/galaxy coevolution model is in efficiently funneling gas down to the SMBH sphere of influence. Gas-rich major mergers of galaxies are an effective means to accomplish this, simultaneously (or near-simultaneously) igniting both a starburst and a luminous AGN . However many observations indicate that AGNs do not prefer merger remnant hosts , and major mergers are likely to fuel only nearby AGNs or the rare population of very luminous quasars .

Violent disk instabilities might provide another solution to the gas inflow problem . These models are supported by the observation that massive galaxies at $z \sim 2$ frequently have most of their star formation in massive ($10^8$-$10^9$ $M_\odot$, $100-500$ pc) clumps within irregular galaxies , unlike the much smoother morphologies of similar-mass star-forming galaxies in the local universe. The kinematics of these clumps indicate that they are not accreted as minor mergers, but instead form as in-situ gravitational instabilities . The high gas fractions observed in $z \sim 2$ galaxies , likely accreted as cold gas along filaments of the cosmic web , would naturally result in highly turbulent star-forming clumps .

Violent disk instabilities might lead to clumps rapidly falling into the center of a galaxy and efficiently fueling a luminous AGN . This remains a point of debate, however, as most simulations lack the resolution and detailed physics to follow the gas all the way into a galaxy's center. Observations indicate that star-forming clumps at $z \sim 2$ typically make up only a small fraction of their galaxy's stellar mass, suggesting that they may be too short-lived to migrate all the way into a galaxy center . Indeed, simulations show that clumps might be destroyed by feedback or exhausted by star formation before reaching the nuclear AGN. On the other hand, other simulations find that clumps are sufficiently long-lived to form central bulges and fuel AGNs even with feedback and mass loss , and the gas between clumps, rather than the clumps themselves, may be more important for large-scale inflows anyway . In support of AGN fueling from violent disk instabilities, found direct observational evidence for a higher AGN fraction in a small sample (14) of $z \sim 0.7$ clumpy galaxies compared to smooth (non-clumpy) extended disks matched in stellar mass and redshift.

Here we extend the search for a connection between clumps and AGNs to $1.3<z<2.4$, using data from the Cosmic Assembly Near-Infrared Legacy Survey . This redshift range is particularly interesting because it represents the peak of both cosmic star formation and the AGN luminosity function . Cold gas inflows and resultant violent disk instabilities are also most prominent in simulations at $z \sim 2$ . And in contrast to the local universe, AGNs at $z>1$ may be unique in growing before a host galaxy bulge develops, with evidence from both evolution in BH-bulge relations and from stacked detection of AGN signatures in low-mass disk galaxies . The gravitational torques provided by violent disk instabilities provide a potential way to fuel an AGN in a disk galaxy even while the bulge is still under construction.

Section 2 describes the combination of Hubble Space Telescope ( *HST*) Wide Field Camera 3 (WFC3) near-infrared (IR) imaging and spectroscopy used to construct and study our samples of $z \sim 2$ galaxies with clumpy, smooth, and intermediate morphologies. We identify AGNs using X-ray emission, emission-line diagnostics, and spatially resolved emission-line ratios: as Section 3 describes, this includes a broad range of AGN properties and avoids contamination by shock-dominated galaxies. We put these selection techniques together in Section 4, and demonstrate that clumpy galaxies are not more likely to host AGNs than smoother galaxies. Section 5 concludes with a discussion of what these results mean for AGN/galaxy coevolution at $z
\sim 2$. Throughout the paper we use a standard $\Lambda$CDM cosmology with $h_0=0.7$.

```latex
\begin{deluxetable*}{ccrrrrrr}
%\rotate
%\tabletypesize{\scriptsize}
  \tablecolumns{8}
  \tablecaption{Galaxy Properties\label{tbl:catalog}}
  \tablehead{
    \colhead{ID} &
    \colhead{Morphology} &
    \colhead{RA} &
    \colhead{Dec} &
    \colhead{$z$} &
    \colhead{$f(\hbox{{\rm H}$\beta$})$} &
    \colhead{$f(\hbox{[{\rm O}\kern 0.1em{\sc iii}]})$} &
    \colhead{$\log(M_*)$} \\
    \colhead{-} &
    \colhead{-} &
    \colhead{(deg)} &
    \colhead{(J2000)} &
    \colhead{-} &
    \colhead{$10^{-18}$~erg~s$^{-1}$~cm$^{-2}$} &
    \colhead{$10^{-18}$~erg~s$^{-1}$~cm$^{-2}$} &
    \colhead{$(\log(M_{\odot}))$} }
  \startdata
 7502 &  clumpy  &  53.07221603 & -27.84033012 &  1.61 & $99.95\pm28.56$ & $83.47\pm24.51$ &  9.63 \\
 7897 &  clumpy  &  53.14831924 & -27.83686638 &  1.90 & $23.49\pm12.22$ & $27.09\pm8.42$ &  9.82 \\
 7907 &  clumpy  &  53.15573502 & -27.83712769 &  1.39 & $18.95\pm13.73$ & $169.79\pm21.69$ &  9.67 \\
 7930 &  clumpy  &  53.15451050 & -27.83648491 &  2.03 & $14.20\pm10.43$ & $43.64\pm10.26$ &  9.41 \\
 8206\tablenotemark{a} &  clumpy  &  53.14358902 & -27.83471107 &  1.99 & $10.82\pm5.73$ & $39.97\pm10.54$ &  9.85 \\
 5023 &  smooth  &  53.10656738 & -27.86481476 &  1.90 & $<3.31$ & $18.62\pm5.92$ &  9.74 \\
 5728 &  smooth  &  53.08922577 & -27.85734558 &  2.03 & $16.80\pm7.95$ & $79.73\pm11.31$ & 10.37 \\
 6278\tablenotemark{b} &  smooth  &  53.06018448 & -27.85304642 &  1.54 & $<2.73$ & $240.61\pm21.63$ & 10.79 \\
 6678 &  smooth  &  53.07507324 & -27.84802055 &  1.73 & $37.84\pm7.38$ & $21.91\pm6.77$ & 10.14 \\
 7806 &  smooth  &  53.16970062 & -27.83803558 &  1.94 & $8.18\pm7.71$ & $38.97\pm11.10$ & 10.37 \\
 5552 &  intermed  &  53.06779480 & -27.85925293 &  1.55 & $15.67\pm14.78$ & $73.40\pm16.57$ &  9.81 \\
 6895 &  intermed  &  53.09495926 & -27.84582710 &  1.56 & $<4.50$ & $49.29\pm8.75$ &  9.44 \\
 7952 &  intermed  &  53.18961334 & -27.83626175 &  2.09 & $39.79\pm13.92$ & $52.87\pm9.57$ &  9.89 \\
 8124 &  intermed  &  53.10293198 & -27.83484268 &  1.38 & $55.60\pm21.62$ & $102.88\pm15.12$ &  9.18 \\
 8307 &  intermed  &  53.07443237 & -27.83325577 &  1.54 & $<6.95$ & $55.32\pm10.17$ &  9.26  \enddata
 \tablecomments{Weak \hbox{{\rm H}$\beta$}\ lines detected below a 1$\sigma$ threshold
   are treated as upper limits (using the 1$\sigma$ error as the
   limit).  The full catalog of 44 clumpy galaxies, 41 smooth
   galaxies, and 35 intermediate galaxies appears as a
   machine-readable table in the electronic version.}
 \tablenotetext{a}{X-ray detected, but in the soft band only and
   consistent with emission from a star-forming galaxy \citep{xue11}.}
 \tablenotetext{b}{X-ray detected and classified as an AGN by
   \citet{xue11}.}
\end{deluxetable*}
```

# Observational Data

We select a sample of 44 clumpy galaxies from the Great Observatories Origins Deep Survey South region of CANDELS. For comparison, we also construct mass-matched samples of 41 smooth (non-clumpy) and 35 intermediate galaxies. All galaxies have $H<24$ (to ensure reliable classification of clumpiness) and have $\hbox{[{\rm O}\kern 0.1em{\sc iii}]}$ detected at the 3$\sigma$ level (for reliable AGN line ratio diagnostics) in the redshift range $1.3<z<2.4$. The effective limits on star formation rate (SFR) and stellar mass ($M_*$) caused by these flux limits are discussed in Section 2.3. Each morphology category has a median redshift of $z=1.85$. Redshifts and emission line measurements come from *HST*/WFC3 grism spectroscopy taken by the 3D-HST survey .

The observational data are described below, with particular attention to the methods for clumpiness classification and AGN identification. The derived data for the clumpy, smooth, and intermediate galaxies are presented in Table .

## Visual Morphologies: Clumpy, Smooth, and Intermediate

Our samples of clumpy, smooth, and intermediate galaxies come from the 4-epoch CANDELS GOODS-S visual classification catalog . Each galaxy has high-resolution *HST* imaging in the rest-frame UV from ACS F606W and F850LP and in the rest-frame optical from WFC3 F125W and F160W . Galaxies were classified on the clumpiness/patchiness grid shown in Figure . Here "clumpiness" is a measure of the number of compact knots, while "patchiness" refers to the more diffuse irregularities within galaxies.

We use visual classification because automated selection of clumps and patchiness remains a difficult computational problem . Inspectors were instructed to focus on the bluest passband for clumpiness classification, as clumps are typically most evident in in the rest-UV . All galaxies in our sample (from the 4-epoch catalog of ) were inspected by at least five classifiers, and we combine the different classifications into a single averaged result for each galaxy.

We use both clumpiness and patchiness as indicators of violent disk instabilities within a galaxy. This is motivated by the presence of both compact/round clumps and diffuse/elongated structures in simulations of $z \sim 2$ galaxies undergoing violent disk instabilities . In addition, if star-forming clumps are enshrouded by dust, their emission would likely be reprocessed in a patchy morphology. Observations show that galaxies appearing clumpy in the UV are merely patchy morphologies in the near-infrared: adding dust to such a system would cause it to appear patchy in the UV rather than clumpy. We also note that find better agreement between automated clump-finding and visual classification when patchiness is included in the visual clumpiness. Therefore we combine both axes in Figure for the identification of clumpy galaxies.

The classifications on each axis are summed to give each galaxy a clumpiness+patchiness value ranging from 0 to 4. Here 0 represents the top left (no clumpiness, no patchiness) and 4 represents the lower right (3+ clumps and maximally patchy). The individual classification results are averaged together to assign "total clumpiness" parameter $C$, which we use to define the galaxies most and least likely to be dominated by violent disk instabilities. We define $C \geq 2$ galaxies as "clumpy," $C<1$ galaxies as "smooth," and $1 \leq C<2$ galaxies as "intermediate." These divisions are chosen such that each category has roughly the same number of galaxies, for ease of comparing AGN fractions across morphologies.

We compare the visual classification $C$ with the automated number of clumps from in Figure . There is a loose correlation between the two methods: the visually-identified smooth galaxies tend to have the fewest automated number of clumps, while the visually clumpy galaxies typically have the most automated number of clumps. There is some scatter due to the imperfection of both methods. include a more detailed comparison of automated clump-finding with the same visual classifications used here, finding $\sim$`<!-- -->`{=html}75% agreement between the two methods (see their Appendix).

Figure shows *HST* $iJH$ (rest UV-optical) color composite images for several representative clumpy, smooth, and intermediate galaxies in our sample. Images of all galaxies are additionally shown in the Appendix (Figures , , and ).

## HST/WFC3 G141 Slitless Grism

The GOODS-S region of CANDELS has near-complete spectroscopic coverage in the near-IR from publicly available 2-orbit *HST*/WFC3 G141 grism observations taken by the 3D-HST survey . For redshifts and line measurements of individual galaxies, the data were reduced using the `aXe` software to produce 2D and 1D wavelength- and flux-calibrated spectra at $1.1<\lambda<1.7\mu$m. We used the `specpro` IDL software[^1] to visually inspect and determine cross-correlation redshifts for the samples of clumpy, smooth, and intermediate galaxies. Spectra with significant ($>$`<!-- -->`{=html}10%) contamination (from neighboring objects) in the + region are rejected from each sample: this occurs for about $\sim$`<!-- -->`{=html}15% of galaxies.

Line fluxes and ratios of  and $\lambda$`<!-- -->`{=html}5007Å were also measured from the `aXe`-reduced 3D-HST spectra. The low resolution of the WFC3 slitless grism ($R \simeq 130$ and 46.5Å/pixel, worse for extended sources) means that the  and $\lambda$`<!-- -->`{=html}4959,5007Å lines are somewhat blended. We simultaneously fit all three lines with Gaussians, constraining each line to be fit within (rest-frame) 50Å ($\sim$`<!-- -->`{=html}2-3 pixels) of the line center. We do not constrain the width of each line because the spatial broadening of the slitless grism frequently causes unusual profiles which can differ for each line. The  line flux is measured directly from the Gaussian fit, while the $\lambda$`<!-- -->`{=html}5007 flux is measured as $3/4$ of the total flux from both blended  Gaussians . If we instead fix the two $\hbox{[{\rm O}\kern 0.1em{\sc iii}]}$ Gaussians to have total fluxes with a 1:3 ratio, the measured line fluxes do not significantly change. Examples of the line fits are shown in Figure .

Uncertainties in the line ratios are computed by re-fitting continua and Gaussians on 10000 realizations of the resampled data. Our sample includes only galaxies with $\lambda$`<!-- -->`{=html}5007Å fluxes measured at the 3$\sigma$ level. We do not set any requirement on the $\hbox{{\rm H}$\beta$}$ line flux measurements, and $\hbox{{\rm H}$\beta$}$ line fluxes less than the 1$\sigma$ error are treated as upper limits (using the 1$\sigma$ error as the limit). In general, the unique asymmetric shape of the blended $\hbox{[{\rm O}\kern 0.1em{\sc iii}]}$ lines enables secure redshifts even when $\hbox{{\rm H}$\beta$}$ is poorly detected.

To study spatially resolved line ratios, we used separate reductions of the WFC3 G141 data from the 3D-HST pipeline, as described in . The 3D-HST reduction method produces superior spatially resolved spectra because it interlaces rather than drizzles: this mitigates the correlated noise associated with drizzling in the final high-resolution ($0\farcs06$/pixel) combined 2D spectra.

## Stellar Masses and Star Formation Rates

We calculate stellar masses and star formation rates for the clumpy, smooth, and intermediate galaxies using the extensive UV / optical / IR photometry in GOODS-S, with 18 bands including 8 with high-resolution *HST* imaging . First, the spectroscopic redshift from the *HST*/WFC3 grism is used to shift the observed photometry to the rest-frame. We then fit the rest-frame spectral energy distribution (SED) with models that include a initial mass function, exponentially declining star-formation histories, and a extinction law. Stellar mass is given by the best-fit model. Example SED fits (for the same galaxies shown in Figure are shown in Figure .

There is some evidence that $z \sim 2$ galaxies are more likely to have constant star-formation histories rather than the exponentially declining models used here. In practice, our best-fit SEDs generally have $\tau$ values which are very long (parameterizing the star formation history as $SFR \sim e^{-t/\tau}$). We also tested the effects of forcing a constant star formation history, and found the masses to change by only $\lesssim$`<!-- -->`{=html}0.1 dex in all cases.

Meanwhile star formation rates (SFRs) are calculated following the method of when mid- and far-IR data are available, and SED-fitting following the method of otherwise. IR emission from an intrinsically luminous AGN might contaminate the IR with a hot dust bump , influencing both the stellar mass and SFR estimates. This occurs for one of our objects: the clumpy galaxy and X-ray AGN 16985, shown in Figure . The remainder of the sample have well-fit galaxy SEDs without evidence for significant AGN contribution. This is unsurprising, as only intrinsically luminous ($log(L_X) \gtrsim 44$) AGNs tend to have hot dust IR emission which dominates over the galaxy emission .

Figure (top panel) shows the mass distributions for each of the clumpy, smooth, and intermediate galaxy samples. The three samples are selected to be matched in stellar mass: starting with the clumpy galaxy sample, we constructed the smooth and intermediate galaxy samples to have a similar mass distribution. All three samples have median and mean $M_* = 10^{9.8}M_\odot$.

Our galaxies are plotted in $SFR-M_*$ space in Figure (bottom panel). For comparison we also show the SFR and $M_*$ of the larger population of GOODS-S galaxies with photometric redshifts calculated by in the same $1.3<z<2.4$ redshift range. Our $z \sim 1.85$ galaxies are limited to $M_* \gtrsim 10^9 M_\odot$ by the $H<24$ selection criterion. The requirement for  detection at the 3$\sigma$ level imposes a flux limit of $f(\hbox{[{\rm O}\kern 0.1em{\sc iii}]}) \gtrsim 5 \times 10^{-17}$ erg s$^{-1}$ cm$^{-2}$ in the 2-orbit WFC3 grism data. Although we estimate SFR from the SED rather than from emission lines, this $\hbox{[{\rm O}\kern 0.1em{\sc iii}]}$ line flux translates roughly to the observed SFR limit of SFR$\gtrsim 10$ M$_\odot$ yr$^{-1}$ (assuming $f(\hbox{{\rm H}$\alpha$}) \sim f(\hbox{[{\rm O}\kern 0.1em{\sc iii}]})$ and using the relation). For the lower-mass half of the sample ($M_*<10^{9.8}M_\odot$), this SFR limit restricts the galaxies to be starbursting and above the star-forming "main sequence" . On the other hand the higher-mass half of the sample ($M_*>10^{9.8}M_\odot$) generally lies on the main sequence for star-forming galaxies at $1.3<z<2.4$. The only notable exception is a few high-mass smooth galaxies below the main sequence which might be quenching.

## X-ray Data

The CANDELS GOODS-S field also contains 4 Ms of *Chandra* X-ray data . The 4 Ms depth corresponds to an on-axis flux limit of $3.2 \times 10^{-17}$ erg $s^{-1}$ cm$^{-2}$, which corresponds to a luminosity limit of $L_X>10^{41.9}$ erg s$^{-1}$ at our sample's median redshift of $z=1.85$. We use the classifications of to separate X-ray galaxies (with X-ray detections consistent with the $L_X-SFR$ relation, ) from the more luminous and hard-spectrum AGNs. The X-ray data for undetected sources in each morphology category are also stacked, as discussed in Section 3.3.

# AGNs in Clumpy Galaxies?

We compare the AGN fraction of clumpy galaxies with the AGN fraction among the mass-matched smooth and intermediate galaxies. AGN detection is accomplished in three different ways: with standard line ratio selection via the "mass-excitation" method , with spatially resolved line ratios, and with X-rays (both for individual sources and stacks of each morphological category).

## Mass Excitation Diagnostic for AGNs

It has long been known that the high-ionization emission of AGNs results in a different emission line signature than observed in typical  regions associated with star formation (SF): in particular, AGN narrow line regions tend to exhibit higher ratios of collisionally excited "forbidden" lines to hydrogen recombination lines . The classic and AGN/SF diagnostics use the ratios of $f(\hbox{[{\rm O}\kern 0.1em{\sc iii}]}\lambda5007)/f(\hbox{{\rm H}$\beta$})$ vs. $f(\hbox{[{\rm N}\kern 0.1em{\sc ii}]}\lambda6584)/f(\hbox{{\rm H}$\alpha$})$ or $f(\hbox{[{\rm S}\kern 0.1em{\sc ii}]}\lambda6718+6731)/f(\hbox{{\rm H}$\alpha$})$: the small wavelength separation of each line pair means that these ratios are effectively insensitive to reddening.

The low resolution of the WFC3 grism means that the  and  lines are not resolved from one another, and so the standard BPT or VO87 diagnostics cannot be used on our data. Instead we use the "mass-excitation" (MEx) method , which uses the $\hbox{[{\rm O}\kern 0.1em{\sc iii}]}/\hbox{{\rm H}$\beta$}$ ratio with the stellar mass $M_*$ to separate AGN and SF galaxies. Functionally, the MEx method uses the correlation between mass and metallicity to separate low-metallicity galaxies from AGN, both of which exhibit similarly high $\hbox{[{\rm O}\kern 0.1em{\sc iii}]}/\hbox{{\rm H}$\beta$}$ ratios. Compared to X-ray or infrared surveys, line ratio selection is often more sensitive to obscured and moderately accreting AGN, but it also less reliable . There is also some debate whether line ratio diagnostics remain applicable in $z>1$ galaxies at all, since higher redshift galaxies may have different gas properties , changing their line ratios in the absence of AGN . also demonstrate that the different selection effects in high- and low-redshift samples have important effects on their observed MEx distributions. For this reason we avoid comparing the $z \sim 1.5$ galaxy samples with local objects. Instead we simply compare clumpy, smooth, and intermediate galaxies all within the same redshift range.

Figure shows the MEx diagram for the clumpy, smooth, and intermediate galaxies. The solid line in each panel shows the empirical line from dividing AGN and SF galaxies, adjusted for the redshift ($z \sim 1.85$) and line flux limit ($f(\hbox{[{\rm O}\kern 0.1em{\sc iii}]}) \gtrsim 5 \times 10^{-17}$ erg s$^{-1}$ cm$^{-2}$): $$y = 0.375/(m-10.4)+1.14~{\rm if}~m \le 9.88$$ $$y = 290.2 - 76.34m + 6.69m^2 - 0.1955m^3~{\rm otherwise.}$$ Here $y=\log(\hbox{[{\rm O}\kern 0.1em{\sc iii}]}/\hbox{{\rm H}$\beta$})$ and $m = \log(M_*/M_\odot)-0.51$. Assuming that all galaxies above this line are AGNs, clumpy galaxies seem to have the highest AGN fraction (12/44, 27%), followed by smooth galaxies (11/41, 27%) and intermediate galaxies (7/35, 20%).

However the simple binary classification of counting objects above and below the AGN/SF line is not entirely appropriate, given that most galaxies have emission lines with composite contribution from both SMBH accretion and  regions. The MEx diagram is much more suited to a probabilistic classification approach, as introduced by . The MEx probabilities, updated for use at $z>1$ by and assuming an error of 0.2 dex in $M_*$, indicate that $39_{-6}^{+8}\%$ of clumpy, $36_{-6}^{+8}\%$ of smooth, and $35_{-7}^{+9}\%$ of intermediate galaxies are AGN-dominated. Given the uncertainty of the MEx diagram at $z>1$ , these probabilities do not necessarily represent absolute AGN fractions, but they remain useful for comparing relative AGN fractions among the three morphology classes. The three AGN fractions are all consistent with one another, and the clumpy galaxies have only marginally ($<$`<!-- -->`{=html}1$\sigma$) more AGNs than smoother galaxies. In the next two subsections we show that the higher line ratios in clumpy galaxies are likely an effect of high ionization in extended regions (due to shocks or dense star-forming clumps) rather than nuclear AGN.

## Spatially Resolved Line Ratios

The traditional BPT and MEx methods use line ratios integrated over the entire galaxy, and thus may be diluted by star formation and/or affected by non-nuclear ionization that has nothing to do with AGN activity. In particular, the violent disk instabilities of clumpy galaxies may lead to high-velocity shocks or dense $\hbox{{\rm H}\kern 0.1em{\sc ii}}$ regions which can mimic AGN-like line ratios . For this reason we use the spatial resolution of the WFC3 grism spectroscopy to go beyond integrated line ratios and investigate line ratio gradients, following . Spatially resolved line ratios can reveal if AGNs selected by integrated line ratios are the product of nuclear AGNs or extended phenomena (like shocks or high-ionization star formation).

The two-dimensional (2D) spectra of individual sources generally lack the signal-to-noise for well-measured spatially resolved line ratios, so we stack spectra of galaxies in each of the clumpy, smooth, and intermediate galaxy samples. Two stacks are constructed for each morphology class: one with all galaxies, and another using only galaxies classified as AGNs on the MEx diagram in Section 3.1. Centers of galaxies are defined by the SExtractor coordinates in the F160W detection image, translated to a spectral trace in the dispersed G141 2D spectrum using well-calibrated polynomials from `aXe` which are typically accurate to $\lesssim$`<!-- -->`{=html}0.2 pixels. The image centroid and resultant spectral trace may not coincide with the brightest (continuum or emission line) flux of a galaxy, particularly for clumpy galaxies: we investigate centering on the brightest flux rather than image centroids in Section 4.1. Each galaxy spectrum is also normalized by its total $\hbox{[{\rm O}\kern 0.1em{\sc iii}]}$ flux so that the stack is equal-weighted and not dominated by the most $\hbox{[{\rm O}\kern 0.1em{\sc iii}]}$-luminous AGNs. (This normalization means that the stacked line ratios are slightly lower than the mean of the individual ratios.) The stacked 2D spectra of all galaxies in each morphology class are shown in Figure . It is immediately evident that the clumpy galaxies are typically larger than the intermediate or smooth galaxies: for a galaxy to be classified as clumpy, it must be large enough for the clumps to be resolved in the *HST* data. Still, the smooth and intermediate galaxies are not so small as to be unresolved, and they have enough signal beyond the central three pixels to have well-measured extended (non-nuclear) spectra. We return to a discussion of the different sizes of the clumpy and smooth galaxies in Section 4.2.

We separately extract nuclear and extended one-dimensional (1D) spectra from the stacked data, shown in Figure . Here "nuclear" is defined as the central 3 pixels, corresponding to $0\farcs18$ across and a $\sim$`<!-- -->`{=html}1 kpc radius at $1.3<z<2.4$. The "extended" spectra are extracted from pixels 3-6 on either side of the trace, translating to an annulus extending radially from $\sim$`<!-- -->`{=html}2-4 kpc. It is possible that some emission from the AGN narrow line region extends into the extended aperture , and the $\sim$`<!-- -->`{=html}1 kpc nuclear region is also likely to include some galaxy starlight. Although our regions will not perfectly disentangle AGN and galaxy light, the nuclear region preferentially includes AGN light and the extended region includes more galaxy starlight. Thus the spatially resolved line ratios are likely to be significantly more sensitive than the integrated line ratios.

```latex
\begin{deluxetable}{lrr}
  \tablecolumns{3}
  \tablecaption{Stacked Spectra Nuclear and Extended $\log(\hbox{[{\rm O}\kern 0.1em{\sc iii}]}/\hbox{{\rm H}$\beta$})$
    \label{tbl:stackratios}}
  \tablehead{
    \colhead{Stack} &
    \colhead{Nuclear} &
    \colhead{Extended} }
  \startdata
    All Clumpy & $3.5\pm0.8$ & $4.5\pm1.2$ \\
    All Smooth & $3.2\pm0.4$ & $2.6\pm1.2$ \\
    All Intermediate & $3.9\pm1.1$ & $3.3\pm1.3$ \\
    MEx AGN Clumpy & $7.5\pm1.5$ & $7.1\pm2.0$ \\
    MEx AGN Smooth & $8.3\pm1.4$ & $5.8\pm1.3$ \\
    MEx AGN Intermediate & $11.8\pm2.1$ & $9.8\pm2.3$ \\
    \hbox{[{\rm O}\kern 0.1em{\sc iii}]}-centered Clumpy AGNs & $5.5\pm1.2$ & $10.5\pm1.8$
  \enddata
\end{deluxetable}
```

The  and  emission lines are measured from each nuclear and extended spectrum following the same method as in Section 2.2: three Gaussians were simultaneously fit to the continuum-subtracted spectra, with the $\lambda$`<!-- -->`{=html}5007 flux given as 3/4 of the total $\lambda$`<!-- -->`{=html}4959+5007 emission. Figure shows the nuclear and extended $\hbox{[{\rm O}\kern 0.1em{\sc iii}]}/\hbox{{\rm H}$\beta$}$ ratios for the stacks of all clumpy, smooth, and intermediate galaxies in the MEx diagram (using the median stellar mass of each sample with small offsets between the points). Figure shows the nuclear and extended line ratios using the stacked spectra of MEx AGN only. The nuclear and extended line ratio measurements for each stack are also presented in Table .

In the stacks of all galaxies in Figure , the nuclear and extended spectra are nearly indistinguishable. This indicates that there is not a dominant AGN population in any of the morphology classes. Nuclear AGNs begin to marginally emerge in Figure in the stacked spectra of MEx AGNs among smooth galaxies. This implies that AGN classification using the updated MEx method is effective for smooth galaxies at $z \sim
2$, with their stacked spectra indicating nuclear AGN. The clumpy galaxies classified as MEx AGN, however, have essentially no difference in nuclear and extended $\hbox{[{\rm O}\kern 0.1em{\sc iii}]}/\hbox{{\rm H}$\beta$}$ ratios. The high $\hbox{[{\rm O}\kern 0.1em{\sc iii}]}/\hbox{{\rm H}$\beta$}$ ratios observed in the integrated spectra of clumpy galaxies are probably caused by extended phenomena like shocks or high-ionization star formation, both of which are likely in dense star forming clumps . Spatial line ratios indicate that clumpy galaxies may actually be *less* likely to host nuclear AGNs than smooth or intermediate galaxies. We further investigate the possibility of off-nuclear AGNs in Section 4.1.

## X-ray AGN

X-rays tend to be the most efficient and least contaminated indicator of AGN activity, and so we also use the X-ray data to quantify the AGN fraction among clumpy, smooth, and intermediate morphology galaxies. All galaxies are matched to the 4 Ms *Chandra* catalog. In total, 5/44 clumpy, 4/41 smooth, and 2/35 intermediate galaxies have X-ray counterparts. X-ray detection is not sufficient for AGN classification, however, as luminous starburst (SFR$\sim$`<!-- -->`{=html}300 M$_\odot$/yr) galaxies can have enough X-ray binaries to exceed an integrated luminosity of $L_X \sim 10^{42}$ . classify the GOODS-S X-ray sources into galaxies consistent with pure star formation or systems likely to be AGN on the basis of their spectral shape, X-ray to optical flux ratio, and X-ray luminosity. These classifications reveal that 2/48 clumpy, 3/41 smooth, and 2/35 intermediate galaxies are likely to be X-ray AGNs. The detection fractions and full-band luminosities of the X-ray AGNs are broadly consistent across the morphology categories, ranging from $10^{42}-10^{44}$ erg s$^{-1}$. X-ray data for the detected galaxies, as well as for the stacked non-detections (discussed below), are presented in Table .

```latex
\begin{deluxetable*}{crrrrc}
  \tablecolumns{6}
  \tablecaption{X-ray Data: Detections and Stacked Non-detections\label{tbl:xray}}
  \tablehead{
    \colhead{Galaxy ID} &
    \colhead{Redshift} &
    \colhead{$f$(0.5-2~keV)} &
    \colhead{$f$(2-8~keV)} &
    \colhead{$f$(0.5-8~keV)} &
    \colhead{Classification} \\
    \colhead{} &
    \colhead{} &
    \colhead{(erg~s$^{-1}$~cm$^{-2}$)} &
    \colhead{(erg~s$^{-1}$~cm$^{-2}$)} &
    \colhead{(erg~s$^{-1}$~cm$^{-2}$)} &
    \colhead{} }
  \startdata
8206  & 1.99 & $2.24 \times 10^{-17}$ & $<1.02 \times 10^{-16}$ & $<7.13 \times 10^{-17}$ & Galaxy \\
8409  & 1.82 & $2.26 \times 10^{-17}$ & $<1.21 \times 10^{-16}$ & $ 5.18 \times 10^{-17}$ & AGN \\
9258  & 2.02 & $1.64 \times 10^{-17}$ & $<1.05 \times 10^{-16}$ & $<6.58 \times 10^{-17}$ & Galaxy \\
9474  & 2.12 & $2.59 \times 10^{-17}$ & $<1.39 \times 10^{-16}$ & $ 7.37 \times 10^{-17}$ & Galaxy \\
10493 & 1.83 & $8.00 \times 10^{-17}$ &  $2.11 \times 10^{-16}$ & $ 2.89 \times 10^{-16}$ & AGN \\
6278  & 1.54 & $4.70 \times 10^{-17}$ &  $6.58 \times 10^{-16}$ & $ 6.75 \times 10^{-16}$ & AGN \\
9493  & 1.47 & $2.72 \times 10^{-17}$ & $<1.14 \times 10^{-16}$ & $ 6.64 \times 10^{-17}$ & Galaxy \\
10650 & 1.31 & $1.28 \times 10^{-16}$ &  $1.33 \times 10^{-16}$ & $ 2.56 \times 10^{-16}$ & AGN \\
16985 & 1.73 & $2.99 \times 10^{-16}$ &  $3.95 \times 10^{-15}$ & $ 4.03 \times 10^{-15}$ & AGN \\
9245  & 2.07 & $1.59 \times 10^{-17}$ &  $9.04 \times 10^{-17}$ & $ 8.01 \times 10^{-17}$ & AGN \\
 18315 & 2.32 & $1.57 \times 10^{-16}$ &  $6.75 \times 10^{-16}$ & $ 8.18 \times 10^{-16}$ & AGN \\
\hline
Stacked Clumpy & -- & $4.3\pm0.9 \times 10^{-18}$ & $<1.1 \times 10^{-17}$ & $9.2\pm2.8 \times 10^{-18}$ \\
Stacked Smooth & -- & $4.0\pm0.9 \times 10^{-18}$ & $<1.6 \times 10^{-17}$ & $8.8\pm2.9 \times 10^{-18}$ \\
Stacked Intermediate & -- & $<4.5 \times 10^{-18}$ & $<9.3 \times
10^{-18}$ & $<9.1 \times 10^{-18}$ \\
\hline
Stacked Clumpy MEx AGN & -- & $5.2\pm1.7 \times 10^{-18}$ & $<2.4 \times 10^{-17}$ & $13.9\pm5.1 \times 10^{-18}$ \\
Stacked Smooth MEx AGN & -- & $<5.2 \times 10^{-18}$ & $<2.6 \times 10^{-17}$ & $<17.7 \times 10^{-18}$ \\
Stacked Intermediate MEx AGN & -- & $<4.6 \times 10^{-18}$ & $<3.3 \times 10^{-17}$ & $<15.7 \times 10^{-18}$
  \enddata
\tablecomments{X-ray fluxes and AGN/galaxy classifications from
  \citet{xue11}.  X-ray stacking following the method in \citet{luo11}.}
\end{deluxetable*}
```

Figure presents the $L_X-SFR$ relationship for X-ray detected sources, for both hard-band (2–8 keV) and soft-band (0.5–2 keV) luminosities. Also shown is the $L_X-SFR$ relation for star-forming galaxies without AGN, $\log(L_X)=[39.6\pm0.4]+\log(SFR)$ , with $L_X$ the full-band luminosity in erg/s and $SFR$ in M$_\odot$/yr. This relation is translated to the soft and hard bands by assuming a spectral index of $\Gamma=1.8$: the resultant relation is consistent with the hard-band $L_X-SFR/M_*$ relation of . A few X-ray sources which are undetected in the hard band lie near the soft-band $L_X-SFR$ relation and might not host AGN, but most have sufficient X-ray luminosities to be robustly considered as having X-ray emission powered by an AGN.

We also performed X-ray stacking of the three samples in the full band, soft band (0.5–2 keV), and hard band using the 4 Ms CDF-S data . X-ray detected galaxies and galaxies near any detected X-ray source (within twice the soft-band 90% encircled-energy aperture radius of the X-ray source) were excluded from the stacking. The final clumpy, smooth, and intermediate samples used in the stacking contain 40, 38, and 33 galaxies, respectively. We adopted the same stacking procedure as described in Section 3.1 of . Briefly, we extracted total (source plus background) counts for each galaxy within a 3-diameter circular aperture centered on its optical position. The corresponding background counts within this aperture were determined with a Monte Carlo approach which randomly (avoiding known X-ray sources) places 1000 apertures within a 1-radius circle of the optical position to measure the mean background. The total counts ($S$) and background counts ($B$) for the stacked sample were derived by summing the counts for individual sources. The net source counts are then given by $(S-B)$, and the S/N is calculated as $(S-B)/\sqrt{B}$. A stacked signal is considered undetected if it has a binomial no-source probability of $P_{\rm
  b}>0.01$ (equivalent to ${\rm S/N}\lesssim2.6\sigma$): for these sources we derive 90% confidence upper limits using the Bayesian method of . Aperture corrections were applied when converting the source count rates to fluxes and luminosities.

The stacked luminosities are shown with the median SFR of each subset in Figure . The small number of objects in each morphology subset means that none are well-detected in the hard band, and only the clumpy and smooth stacks are detected in the soft-band. Due to these poor hard band detections, the hardness ratios are unconstrained in all three stacks. The stacked X-ray luminosities (and upper limits) for all three samples lie only marginally above the $L_X-SFR$ line for SF galaxies, indicating a weak or sub-dominant AGN contribution in the X-ray undetected galaxies. This is consistent with the location of most galaxies in the star-forming region of the MEx diagram (Figure ), and with the lack of nuclear AGN in the resolved line ratios among all galaxies (Figure ). We also stack the MEx AGNs in each morphology class, but the numbers of objects are too small to make any meaningful conclusions.

Clumpy, smooth, and intermediate galaxies all have similar detected X-ray AGN fractions ($\sim$`<!-- -->`{=html}5%) and stacked X-ray luminosities for undetected sources ($\log(L_X) \sim 41.4$). In agreement with the spatially resolved line ratios, the X-ray data show essentially no differences for clumpy, smooth, and intermediate morphology galaxies at $z \sim 2$.

# Discussion

Both spatially resolved line ratios and X-ray data (detections and stacked non-detections) indicate that there is little or no difference in AGN fraction between clumpy galaxies and those with smoother morphologies. While integrated line ratios in the MEx method suggest a slightly higher AGN fraction among clumpy galaxies, the lack of a nuclear AGN signature suggests that this is due to extended phenomena like dense star-forming clumps rather than black hole growth. Galaxies at $1.3<z<2.4$ in each morphology category have similar AGN fractions of $\sim$`<!-- -->`{=html}5–30%, depending on whether X-ray data or spatially resolved line ratios are used for AGN selection.

We discuss several implications of these results below.

## Off-Nuclear AGN in Clumpy Galaxies?

In Section 3.2 we argued that the high $\hbox{[{\rm O}\kern 0.1em{\sc iii}]}/\hbox{{\rm H}$\beta$}$ ratio in the extended regions of clumpy galaxies indicates shocks or high-ionization  regions rather than AGN activity . However, it is instead possible that some clumpy galaxies host off-nuclear AGN . We test for the presence of off-nuclear AGN by creating a 2D stack of clumpy galaxies classified as MEx-AGNs, centering on the $\hbox{[{\rm O}\kern 0.1em{\sc iii}]}$ emission lines rather than the continuum trace. The brightest knot of $\hbox{[{\rm O}\kern 0.1em{\sc iii}]}$ emission would be most likely to correspond to any off-nuclear AGN[^2]. The "nuclear" region of the $\hbox{[{\rm O}\kern 0.1em{\sc iii}]}$-centered stack corresponds to the brightest $\hbox{[{\rm O}\kern 0.1em{\sc iii}]}$ knots rather than the galaxy center, and should therefore be where off-nuclear AGNs are most likely to reside.

Figure compares the nuclear and extended $\hbox{[{\rm O}\kern 0.1em{\sc iii}]}/\hbox{{\rm H}$\beta$}$ ratios for the continuum-centered stack with the $\hbox{[{\rm O}\kern 0.1em{\sc iii}]}$-centered stack of clumpy galaxies. The right panel suggests that the brightest $\hbox{[{\rm O}\kern 0.1em{\sc iii}]}$ knots are somewhat *l*ess likely to host AGNs than the galaxy centers. The lack of widespread off-nuclear AGNs is also supported by the X-ray data, which indicate that clumpy galaxies do not host a greater number of X-ray AGNs (whether nuclear or off-nuclear) than smoother mass-matched galaxies. We conclude that there is unlikely to be a large population of off-nuclear AGNs in clumpy galaxies misclassified by the spatially resolved line ratio method.

## Comparison to Lower Redshift

At first glance, our $z \sim 1.85$ result is in apparent contrast to the preference for AGN in clumpy galaxies at $z \sim 0.7$ observed by . Our clumpy galaxy sample has very similar images and morphologies to theirs, as seen in a comparison of Figure with Figure 4 of . And although we study a much larger sample (three times more clumpy galaxies), performed a comprehensive bootstrap analysis which indicates that sample size is very unlikely to lead to the difference. Our studies differ greatly, however, in the comparison samples of smooth galaxies.

required that both their clumpy and smooth galaxy samples be extended disks. Meanwhile we impose no requirements on morphology beyond the clumpiness classification, and Figure demonstrates that our $z \sim 2$ clumpy and smooth galaxies are quite different in size and index. These quantities were measured using GALFIT and presented by . Clumpy galaxies are much more likely to be extended disks, with a median semi-major axis of $0\farcs46$ and 41/44 galaxies having disk-dominated Sérsic indices ($n<2.5$). The smooth galaxies are still mostly disk-dominated, as 30/41 have $n<2.5$, but they are significantly smaller (with a median semi-major axis of $0\farcs20$) and have higher Sérsic indices.

It is interesting that so few smooth extended disks are observed in CANDELS and 3D-HST. The data are sensitive to the $z \sim 2$ star-forming main sequence at $M_*/M_\odot>10^{9.8}$ (Figure ), and so if smooth extended disks exist they must be low-mass and/or weakly star-forming. and argue that ordered disks are rare in $z>1$ galaxies of any mass, suggesting that the lack of smooth extended disks at $z \sim 2$ is a consequence of galaxy evolution rather than selection effects. Regardless, the lack of smooth extended disks at $z \sim 2$ means we cannot perform the same comparison between clumpy and smooth extended disks as at $z \sim 0.7$. Instead we conclude that clumpy extended galaxies have the same AGN fraction as more compact and higher-Sersic smooth galaxies matched in stellar mass.

## Fueling AGNs at $z \sim 2$

Our results demonstrate that clumpy galaxies at $z \sim 2$ are no more likely to host AGNs than are smooth galaxies of the same stellar mass and similar star formation rate. What does this indicate about AGN fueling modes at $z \sim 2$?

argue that compact star-forming galaxies (also called blue nuggets) might actually be constructed by a history of violent disk instabilities. Our smooth galaxies, which tend to be more compact and bulgy than the clumpy extended disks, have sizes and Sersic indices consistent with this category . The observed similarity in AGN fraction between $z
\sim 2$ clumpy and smooth galaxies would then indicate that nuclear inflow onto the AGN must persist until after the galaxy becomes stable and the star-forming clumps are no longer visible. In this scenario the violent disk instabilities are responsible for building a gas reservoir in the galaxy's center. Since the actual nuclear inflow onto the AGN persists into the smooth compact galaxy phase, it must not depend on violent instabilities in the larger disk. So while the gas reservoir is deposited by violent disk instabilities, the AGN fueling within $\sim$`<!-- -->`{=html}1 kpc must instead be driven by secular processes.

On the other hand, the disk-dominated nature of the smooth galaxies might imply that they are unlikely to have experienced violent mixing in the past from extreme disk instabilities (or major mergers). It is also not clear that the smooth galaxies are older than the clumpy galaxies, given their similar position in the SFR-M\* plane (excepting the most massive galaxies: see Figure ). This might imply that even the large-scale fueling occurs via secular processes, perhaps due to stochastic turbulence from the high gas fractions which are common in $z \sim 2$ galaxies . Indeed, predict that this turbulence-driven stochastic fueling is the dominant mode of SMBH growth for $z \sim 2$ galaxies with $M_{BH}<10^{7}M_\odot$ (or $M_*<10^{10}M_\odot$): the same mass range occupied by much of our sample. Assuming that the smooth and clumpy morphologies in our sample trace different fueling modes, then our data best support the scenario of , where AGN growth depends not on fueling mode, but on gas fraction.

# Summary

We use spatially resolved line ratios and X-ray data to demonstrate that $z \sim 1.85$ clumpy galaxies are no more likely to host AGNs than are smoother galaxies at the same redshift. While integrated line ratios indicate a higher AGN fraction in clumpy galaxies, the spatially resolved emission lines show that this is likely due to extended phenomena (like shocks or high-ionization $\hbox{{\rm H}\kern 0.1em{\sc ii}}$ regions) rather than nuclear (or off-nuclear) AGNs. While the smooth galaxies have the same masses and nearly the same star formation rates, they are somewhat smaller and not quite as disk-dominated as the clumpy galaxies. Whatever process drives AGN fueling in the smooth galaxies is just as efficient as are violent disk instabilities in clumpy galaxies at $z \sim 2$.

We thank Frederic Bournaud for valuable discussion which significantly improved this work. Support was provided by NASA through Hubble Fellowship grant \#51330 awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA under contract NAS 5-26555. JRT and the authors from UCSC also acknowledge support form NASA HST grants GO-12060.10-A and AR-12822.03, Chandra grant G08-9129A, and NSF grant AST-0808133. SJ acknowledges financial support from the E.C. through an ERC grant StG-257720. This work made use of the Rainbow Cosmological Surveys Database, which is operated by the Universidad Complutense de Madrid (UCM), partnered with the University of California Observatories at Santa Cruz (UCO/Lick, UCSC).

# Galaxy Images

Color-composite *HST* $iJH$ images of all clumpy, smooth, and intermediate galaxies are shown in Figures , , and .

[^1]: The specpro package is available at `http://specpro.caltech.edu`

[^2]: Our assumption that the brightest $\hbox{[{\rm O}\kern 0.1em{\sc iii}]}$ knot corresponds to an AGN is inaccurate if there is large-scale dust asymmetrically obscuring an AGN narrow-line region.


abstract: Non-randomized treatment effect models are widely used for the assessment of treatment effects in various fields and in particular social science disciplines like political science, psychometry, psychology. More specifically, these are situations where treatment is assigned to an individual based on some of their characteristics (e.g. scholarship is allocated based on merit or antihypertensive treatments are allocated based on blood pressure level) instead of being allocated randomly, as is the case, for example, in randomized clinical trials. Popular methods that have been largely employed till date for estimation of such treatment effects suffer from slow rates of convergence (i.e. slower than $\sqrt{n}$). In this paper, we present a new model coined SCENTS: Score Explained Non-Randomized Treatment Systems, and a corresponding method that allows estimation of the treatment effect at $\sqrt{n}$ rate in the presence of fairly general forms of confoundedness, when the 'score' variable on whose basis treatment is assigned can be explained via certain feature measurements of the individuals under study. We show that our estimator is asymptotically normal in general and semi-parametrically efficient under normal errors. We further extend our analysis to high dimensional covariates and propose a $\sqrt n$ consistent and asymptotically normal estimator based on a de-biasing procedure. Our analysis for the high dimensional incarnation can be readily extended to analyze partial linear models in the presence of noisy variables corresponding to the non-linear part of the model, where the noise can be correlated with the variables corresponding to the linear part. We analyze two real datasets via our method and compare our results with those obtained by using previous approaches. We conclude this paper with a discussion on some possible extensions of our approach.
author: Debarghya Mukherjee, Moulinath Banerjee and Ya'acov Ritov

Department of Statistics, University of Michigan
bibliography: mybib.bib
title: Estimation of a score-explained non-randomized treatment effect in fixed and high dimensions

# List of covariates

```latex
\fbox{%
\begin{tabular*}{0.9\textwidth}{*{2}{l}}
vshr-islam 1994 & Islamic vote share 1994 \\
 \hline
partycount & Number of parties receiving votes 1994 \\
 \hline
lpop1994 & Log Population in 1994 \\
 \hline
merkezi & District center \\
 \hline
merkezp & Province center \\
 \hline
subbuyuk & Sub-metro center \\
 \hline
buyuk & Metro center \\
 \hline
ageshr19 & Population share below 19 in 2000 \\
 \hline
ageshr60 & Population share below 60 in 2000 \\
 \hline
sexr & Gender ratio in 2000 \\
 \hline
shhs & Household size in 2000 \\
 \hline
i89 & Indicator to Islamic mayor in 1989 \\
 \hline
 partycount & Number of parties in the election 1994 \\
 \hline
 hischshr1520m & Share Men aged 15-20 with High School Education in 2000
\end{tabular*}}
```
```latex
\fbox{%
\begin{tabular*}{0.9\textwidth}{*{2}{l}}
Col. names & Explanation  \\
\hline
hsgrade\_pct & High school grade in percentage \\
\hline
totcredits\_year1 & Total credits taken in first year \\
\hline
loc\_campus1 & Indicator whether the student in from Campus 1 \\
\hline
loc\_campus2 & Indicator whether the student in from Campus 2 \\
\hline
male & Indicator of whether the student is male  \\
\hline
bpl\_north\_america & Whether the birth place in North America \\
\hline
age\_at\_entry & Age of the student when they entered the college \\
\hline
english & Indicator of whether the student is native english speaker \\
\end{tabular*}}
```