{'fraction_non_alphanumeric': 0.03320867397228504, 'fraction_numerical': 0.01865541833362324, 'mean_word_length': 4.454538266919671, 'pattern_counts': {'":': 0, '<': 7, '<?xml version=': 0, '>': 5, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The Role of Twist in Kinked Flux Rope Emergence and Delta-Spot Formation

It has been observationally well established that the magnetic configurations most favorable for producing energetic flaring events reside in delta-spots, a class of sunspots defined as having opposite polarity umbrae sharing a common penumbra. They are frequently characterized by extreme compactness, strong rotation and anti-Hale orientation. Numerous studies have shown that nearly all of the largest solar flares originate in delta-spots, making the understanding of these structures a fundamental step in predicting space weather. Despite their important influence on the space environment, surprisingly little is understood about the origin and behavior of delta-spots. In this paper, we perform a systematic study of the behavior of emerging flux ropes to test a theoretical model for the formation of delta-spots: the kink instability of emerging flux ropes. We simulated the emergence of highly twisted, kink-unstable flux ropes from the convection zone into the corona, and compared their photospheric properties to those of emerged weakly twisted, kink-stable flux ropes. We show that the photospheric manifestations of the emergence of highly twisted flux ropes closely match the observed properties of delta-spots, and we discuss the resulting implications for observations. Our results strongly support and extend previous theoretical work that suggested that the kink instability of emerging flux ropes is a promising candidate to explain delta-spot formation, as it reproduces their key characteristics very well.


INTRODUCTION
One of the most important features of the solar atmosphere, from the standpoint of space weather forecasting, is the subset of sunspots known as δ-spots. Classified as having opposite polarity umbrae that share a common penumbra (Künzel 1965), δ-spots exhibit, by definition, extremely compact photospheric flux distributions. They are often observed to rotate rapidly after emergence (Kurokawa 1987(Kurokawa , 1991Tanaka 1991;López Fuentes et al. 2000, 2003, contain highly twisted magnetic fields (López Fuentes et al. 2000, 2003Holder et al. 2004), and violate Hale's law (Smith & Howard 1968;Zirin & Liggett 1987). Hale's law states that for a given 11-year solar cycle, the majority of sunspot pairs observed in the Northern (Southern) solar hemisphere will be oriented with the negative (positive) polarity spot to the East (West). This ordering flips with each successive solar cycle. δ-spots are overwhelmingly associated with explosive flares, being responsible for > 90% of Xclass flares during the last two solar cycles (Tanaka 1991;Shi & Wang 1994;Sammis et al. 2000;Guo et al. 2014). Conversely, over 90% of the δ-spots studied by Tanaka (1991) were flare active. The details of why δ-spots manifest their unique properties, however, are poorly understood. Understanding the observed behavior of δ-spots is critical for predicting their explosive flares.
Numerical studies explaining δ-spot formation have typically come in four varieties: 1) The kink instability of emerging flux ropes (Linton et al. 1998Fan et al. 1999;Takasao et al. 2015;Toriumi & Takasao 2017), in which a highly twisted flux rope emerges with a deformed axis, producing a complex active region (AR) that manifests a strongly sheared polarity inversion line across the entire AR. Takasao et al. (2015) simulated one such region, NOAA 11429, and found that it was likely formed through the emergence of a kink-unstable flux rope. 2) The splitting off of a portion of a subsurface flux rope to create small parasitic polarities (Jaeggli 2016;Toriumi & Takasao 2017).
In this paper, we will focus on the first mechanism, i.e., the generation of δ-spots via the kink instability of emerging flux ropes. Much of the theoretical underpinning of the theory of the kink instability of twisted flux ropes was developed for applications in tokamak physics by Shafranov (1957) and Kruskal et al. (1958), and then in solar physics (Gold & Hoyle 1960;Anzer 1968;Hood & Priest 1980;Mikic et al. 1990). In highly twisted flux ropes, the kink instability converts twist (rotation of field lines about the flux rope's central axis) into writhe (a deformation of the axis itself; Moffatt & Ricca 1992). The kink instability sets in when the destabilizing magnetic pressure of the azimuthal field encircling the flux rope overwhelms the stabilizing magnetic tension in the central axial field. This occurs when flux ropes are highly twisted (e.g., Shafranov 1957;Kruskal et al. 1958), and, crucially for interpreting observations and predictive modeling, the helical sense of the kink will have the same sense as that of the flux rope twist (Linton et al. 1996). Linton et al. (1999) modeled the kinking of a twisted flux rope and argued that, if such a flux rope were to emerge through the photosphere, the kink would cause it to emerge with a significantly tilted configuration (suggestive of anti-Haleness), then rotate and reorient itself into a Hale configuration, in agreement with observed evolution of δ-spots (López Fuentes et al. 2003). Taken together, simulations of highly twisted flux ropes and the analytical work of Linton et al. (1996Linton et al. ( , 1998Linton et al. ( , 1999 suggest that the kink instability of emerging flux ropes is a promising mechanism to explain the formation and behavior of a significant proportion of δ-spots.
An open question is: what is the initial twist on these emerging flux ropes? One way to estimate flux rope twist is by measuring the twist parameter α = (∇ × B)/B. At the photospheric level, α can either be measured directly through its normal component, α n = J n /B n , with n denoting the component normal to the photosphere, J the current density and B the magnetic field (e.g., Leka et al. 1996) or indirectly by fitting a coronal constantα force free model to the observations (often denoted α best , e.g., Pevtsov et al. 1994Pevtsov et al. , 1995Longcope et al. 1998Longcope et al. , 1999Holder et al. 2004). Observations of α n averaged over each individual AR area find that α n ranges from ±1 × 10 −10 to±1 × 10 −7 m −1 (Leka et al. 1996;Otsuji et al. 2015), while linear force free fits of α best find a similar range of 0 to 5 × 10 −8 m −1 (Pevtsov et al. 1994(Pevtsov et al. , 1995. Analytical arguments suggest that the twist must be large enough for the flux ropes to maintain coherence inside the convection zone (Emonet & Moreno-Insertis 1998). Longcope et al. (1999) found that if this minimum twist were transferred from the convection zone directly to the photosphere, the corresponding peak value of α should be larger than the observed average ranges of α n or α best . Numerical models of emerging flux ropes have typically used initial twists larger than that required to survive the convection zone, and the twist needed for a kink is yet higher still. The question, therefore, is how to explain the relatively low observed values of photospheric α n .
Several theoretical arguments can help resolve this difficulty. First, both α and α n vary significantly across the initial subsurface flux rope, with α peaking at the center and vanishing at the edges, and α n even reversing sign beyond a certain radius. Thus, taking the average of either quantity will lead to values much smaller than the peak value. Second, the assumption that α (or α n ) is transferred directly from the convection zone to the photosphere is unlikely to hold up during the dynamic emergence of the flux rope, as significant reconnection, diffusion, and kinking are expected to reduce the twist on the emerging flux rope. It has also been demonstrated in several flux emergence simulations (e.g., Fan 2009;Leake et al. 2013), that a significant fraction of the twist, in the form of α, is trapped below the photosphere, so that the observed local α at the photosphere is about an order of magnitude lower than the local α in the convection zone portion of that same field line.
One consequence of this reasoning is that it implies that it may not be possible to infer the flux rope's initial twist directly from α n or α best . Nevertheless, it is clear that some ARs do have significantly larger values of α best or α n than others (Pevtsov et al. 1994(Pevtsov et al. , 1995Longcope et al. 1998;Holder et al. 2004). The ARs showing the largest values of either α n or α best have not been system-atically studied in observations, but an individual study found α best = 2.5 × 10 −8 m −1 in δ-spots as compared to 1.1 × 10 −8 m −1 in typical ARs (Holder et al. 2004). This suggests that the flux ropes that form δ-spots are more twisted than typical emerging flux ropes, and may explain these higher values of α best . In addition, due to the kink instability, highly twisted flux ropes are expected to manifest several other features at the photospheric level that suggest they may be the source of δ-spots on the photosphere.
The kink-instability model of δ-spot formation predicts four testable consequences for the photospheric and coronal properties of the resulting ARs.
First, since the kink instability occurs only for highly twisted flux ropes, the emerging AR should exhibit high twist. Thus, the kink-instability model naturally explains the high twists observed in δ-spots (Holder et al. 2004).
Second, the conversion of twist into writhe manifests as a strong rotation of the flux distribution on the photosphere. Crucially, the sense of this rotation will be determined by the sign of the flux rope twist , wherein a right (left) handed twist will give a right (left) handed kink, which, upon passing through the photosphere, will translate into a counter-clockwise (clockwise) rotation. This provides a testable prediction for δspot formation. Thus, the kink of the emerging flux rope is a promising mechanism for explaining the observed rotation of δ-spots on the photosphere (Kazachenko et al. 2010;Vemareddy et al. 2016;Wang et al. 2016).
Third, the kinking motion should break the cylindrical symmetry of the flux rope about its initial axis, moving the flux rope legs off of the initial axis. In its nonlinear phase, the kink instability should cause the rising loop to fold over on itself, causing the legs of the rising flux rope to be held very close to each other, or even, in extreme cases, to form a knot (Linton et al. 1998). The predicted result is that, on the photospheric level, opposite polarity regions of the emerging AR would not separate to any appreciable distance. As a result, the 'compactness' of the photospheric flux distribution in δ-spots may readily explained by the kink instability model of δ-spot formation (Howard 1996;López Fuentes et al. 2000, 2003. Finally, at the coronal level, the highly twisted field in the kinked structure should contain a tremendous amount of free magnetic energy, making it susceptible to intense energy release. This would then explain the association of δ-spots with X-class flares (Tanaka 1991;Shi & Wang 1994;Sammis et al. 2000;Guo et al. 2014).
The picture that emerges, therefore, is that at the photosphere, kinked flux ropes are expected to be compact (if they have kinked enough to bring the legs of the flux rope close together), to be twisted (since kinking requires high twist), to rotate (with the kinking motion manifesting itself as a rotation on the photosphere) and, at the coronal level, to flare (due to high levels of twist, which corresponds to free energy, and knots, which represents topological complexity). This is consistent with the observed properties of δ-spots, which are known to be compact, twisted, and strongly rotating and flaring. Thus, several predictions of the kinking flux rope model are eminently testable via numerical simulations. But although this is strong theoretical support for the kink instability model of δ-spot formation, very few numer-ical studies have tested this idea. With the exception of a few simulations (Matsumoto et al. 1998;Fan et al. 1999;Takasao et al. 2015;Toriumi & Takasao 2017), the overwhelming majority of flux emergence studies have focused on non-kinking, weakly twisted flux ropes. Furthermore, very few authors Murray et al. 2006;Sturrock & Hood 2016) have performed parameter studies to investigate how varying the initial flux rope parameters affects the photospheric flux distribution in emerging flux ropes. Of these, only one ) used a kink-unstable flux rope, and this study did not model the actual emergence. Thus, although significant theoretical work has argued that δ-spots are the photospheric manifestations of the emergence of kinked flux ropes rising from the convection zone into the corona, no systematic numerical study of the emergence of such kinked flux ropes has been undertaken.
In this paper, we perform a parameter study of the emergence of twisted flux ropes from the convection zone into the corona. We vary only a single parameter: the flux rope's initial dimensionless twist, which determines the stability of the flux rope. We study twists ranging from kink-stable to marginally stable to unstable and explore the consequences of the emergence of these flux ropes at the photospheric level.
In these equations, ρ is mass density, T is temperature, P is thermal pressure, γ is the ratio of specific heats, v is velocity, B is magnetic field, g is the gravitational acceleration, η is the magnetic diffusivity and t is time. We close the equations via the ideal gas equation, where R = 8.26×10 7 dyn cm K −1 g −1 is the gas constant. In all of the simulations presented below, we choose g = 2.7 × 10 4 cm s −1 and η = 10 10 cm 2 s −1 . We set up a stratified atmosphere, shown in Figure 1, with a temperature profile of the form in the convection zone, with x the vertical direction. The temperature gradient dT /dx has the constant value −T ph / cz , and the temperature T varies on the characteristic length scale cz . We use the exact solution for an adiabatically stratified Cartesian atmosphere: (2.7) This profile attaches to the bottom of a constant temperature region having temperature T ph starting at height x ph , which we take to be the photosphere. At the top of the constant (minimum) temperature region, at height x tr , the transition region attaches, having a profile of the form where T cor is the asymptotic temperature of the corona (x − x tr tr ) for a characteristic length scale tr . Background magnetic-field free hydrostatic density and pressure profiles are then calculated from Equations 2.2, 2.5, 2.6 and 2.8, with v = B = 0 initially and the photospheric temperature and pressure pinned at T ph = 5.1 × 10 3 K, and P ph = 1.14 × 10 5 dyn cm −2 (cf. Figure 1). These choices set the photospheric density, ρ ph = 2.7 × 10 −7 g cm −3 , pressure scale height H p = RT ph /g = 150 km, sound speed C s = γP/ρ = 8.4 km s −1 , and convection zone length scale cz = 0.383 Mm,. We choose x ph = 0.0 Mm, x tr = 1.75 Mm, tr = 0.3 Mm and T cor = 10 6 K. The simulation box has size L x ×L y ×L z = 25×20×20 Mm, with x representing the vertical direction and y the direction parallel to the initial flux rope axis. Our boundaries are periodic in the horizontal directions, closed at the bottom and open at the top. Our grid is specified such that the smallest grid cell is 0.039 Mm in each direction. We placed a flux rope with magnetic field profile B y (r) = B 0 e −r 2 /a 2 B φ (r) = ζ r a B y (r) (2.9) in pressure balance at a depth x = −d. Here r is the radial distance from the initial flux rope axis, φ is the azimuthal coordinate, d = 3 Mm, B 0 = 6.5 kG is the flux rope's initial magnetic field strength, and a = 0.3 Mm is its radius, so that the total axial flux in the flux rope is Φ 0 = 1.8 × 10 19 Mx. The parameter ζ is the dimensionless twist of the flux rope. Twisted flux ropes are kink unstable when their twist ζ exceeds the critical twist ζ c . For the profile given in Equation 2.9, Linton et al. (1996) showed that ζ c = 1.
To ensure hydrostatic equilibrium is maintained when this flux rope is added to the atmosphere, the field-free atmospheric pressure p 0 (x) is modified by an amount with c = 3 × 10 10 cm s −1 the speed of light. The equilibrium pressure is then To initialize the emergence, we impose a density perturbation of the form (Fan 2009): where λ = 1.2 Mm, so that ρ eq (x, y, z) = ρ 0 (x) + ρ 1 (x, y, z). (2.14) This generates an initial buoyant section of the flux rope, centered at y = 0 and with an extent of ∆y ≈ 2λ. We perform five simulations with these parameters, varying the value of the dimensionless twist ζ, while keeping the radius fixed. We set ζ = {0.5, 1, 1.5, 2, 4}, varying its value from kink-stable to marginally stable to kinkunstable.

RESULTS
The photospheric properties of each of the simulations are both qualitatively and quantitatively different. We first present qualitative results from each simulation individually before comparing the different cases quantitatively.
3.1. ζ = 0.5 3.1.1. Subsurface Evolution The first simulation we present is the initially kink stable flux rope, which starts with ζ = 0.5. The density perturbation creates a buoyant section in the middle of the flux rope, causing it to rise in an Ω-shaped loop, seen in the left panels of Figure 2. Striations, shown in closeup in Figure 3, form initially along the rising apex, but later extend down into the legs of the flux rope. These striations appear to form when field lines move apart to allow plasma to pass through, suggestive of the magnetic Rayleigh-Taylor instability (Duan et al. 2018). The topdown view shown in the right panels of Figure 2 demonstrates that there is a slight rotation of the apex of the rising loop out of its initial plane, suggestive of kinking behavior. This is consistent with the finding that even flux ropes with ζ < 1 initially may, during the course of their rise, become kink-unstable , though this is expected to occur due to an increase in the flux rope's radius, which is not observed here.

Photospheric Evolution
The flux rope breaks through the photosphere at t ≈ 2955 s, producing two opposite polarity regions whose centroids would be connected by a line oriented approximately along the z-direction, as seen in Figure 4a. The flux rope's axial field, oriented approximately in the ydirection throughout the emergence, is encircled by an azimuthal field which, at the apex of the loop, is oriented partially along the z-direction. It is this field that emerges first, producing the polarities seen in Figure 4a. The striations that formed along the flux rope result in the formation of the second, smaller, set of polarities seen in Fig. 4a when the flux rope hits the photosphere. The large and small polarities combine, rotate slightly and spread out as the legs of the flux rope break through the photosphere, forming two polarities with B x ∼ ±500 G at around t = 3295 s ( Figure 4b). The flux weighted centers of mass of each polarities are depicted with the red circles in panel b. They are defined by with the distance between them being Here B ± x is the positive/negative magnetic field normal to the photospheric surface S ph and s = yŷ + zẑ is the radial coordinate along the photosphere. L ± was evaluated by integrating over the entire active region, defined as locations at , but the results presented below were qualitatively unchanged for B thresh x = 20 G. In Figures 5a,b we plot several field lines traced from the photosphere, and view them from the side and from above. These field lines are colored by the component of their field along the original axial (y) direction. The field that emerges into the corona is weakly twisted, and a sheared arcade is evident between the two polarities later in the emergence ( Figure 5c). As the legs of the flux rope emerge, the two main polarities continue to separate, and weak flux towards the edges of the flux rope becomes scattered between the two main polarities ( Figure 4c). The vectors overplotted in Figure 4d show the horizontal field at x = 0. The horizontal field is twisted in the two main polarities, and is much stronger than in the region between the two polarities ( Figure 4d). As can be seen from Figure 5d, the resulting configuration is approximately an inverted U-loop in the corona, with plasma flowing down the legs of the emerged flux rope, and with the apex of each field line comprising the inverted U-loop rising higher into the corona.
3.2. ζ = 1 3.2.1. Subsurface Evolution Next, we present the emergence of an initially marginally kink-stable ζ = 1 flux rope. The subsurface evolution of this flux rope, shown in Figure 6, is qualitatively similar to that of the ζ = 0.5 case. An Ω-loop rises through the convection zone, with striations forming first along the rising loop and then developing along the flux rope legs. The Ω-loop starts to rotate gently around t ≈ 1670 s, and is slightly offset from the z = 0 axis when it hits the photosphere, suggesting that this flux rope may be undergoing the kink instability.

Photospheric Evolution
The flux rope breaks through the photosphere at t ≈ 3625 s forming two well defined polarities (Figure 7a), whose centroids are again connected by a line oriented approximately along the z-direction due to the twist field emerging first. Figures 8a,b show a set of field lines soon after the initial emergence from the side and from above, respectively, again colored by the axial field, B y . They show the two polarities being formed by the apex of the emerging flux rope. This flux rope emerges much more coherently than the previous flux rope because the larger value of ζ helps prevent the break up of the flux during its rise through the stratified atmosphere. These polarities gently rotate and separate, with magnetic field strengths approaching B x ∼ ±1 kG around t = 4190 s ( Figure  7b). Interestingly, small opposite polarity regions (yellow arrows) are visible near the primary polarities (blue arrows). At the coronal level, this rotation and separation results in the appearance of a flux rope comprised of sigmoidal field lines along with overlying quasi-potential loops ( Figure 8c). Dispersed flux develops near and between the two main polarities as they continue to separate. The horizontal field is primarily twisting around the two main polarities, and is extremely weak between them at t = 4990 s (Figures 7c,d). The end result is, like the ζ = 0.5 case, an inverted U-loop coronal arcade ( Figure 8d), with plasma draining along the loop while the loop itself expands upward, carrying plasma with it.
3.3. ζ = 1.5 3.3.1. Subsurface Evolution We now present a simulation showing an emerging flux rope whose initial twist exceeds the kink instability threshold. The rise of this flux rope through the convection zone is shown in Figure 9. The initial Ω-loop formed by the density perturbation quickly begins to coherently rotate out of its initial plane, converting twist into writhe. This rotation is much more well-defined than in the previous two cases. Around t ≈ 2000 s, the apex of the Ω-loop becomes oriented along the z-direction, meaning that it has rotated much further than the two lower twist cases. Striations again form along the rising Ω-loop and extend down the legs of the rising flux rope, although they are not as clearly visible in this simulation as in the previous two.

Photospheric Evolution
The flux rope begins to emerge through the photosphere at t ≈ 2265 s as a pair of elongated, crescentshaped polarities ( Figure 10a). Figures 11a,b show the field lines of the flux rope at this time from the side and from above, the former showing the shape formed as a result of the kink. The regions of red B y indicate locations where the axis of the flux rope has completely reversed from its initial direction as a result of its deformation during the kink. As the emergence proceeds, two pairs of opposite polarities form at the photosphere, and eventually develop into two primary polarities and two secondary, elongated polarities (denoted by the blue and yellow arrows, respectively, in Figure 10c). Low-lying sheared arcades connecting the adjacent primary and secondary polarities can be seen in Figure 11c, along with sigmoidal field lines connecting the two primary polarities. Overlying quasi-potential fields are also apparent. The offset of these two pairs of polarities from each other in the z-direction is a result of the two legs of the flux rope becoming offset from each other due to the kink. As a result, the sigmoidal field lines in Figure 11c are oriented primarily along the z-direction, rather than the y-direction, as was observed in the previous two cases. The field lines at this location, comprised of the outer part of the apex of the emerging loop, are concave-down, and plasma is able to drain down the loops, allowing this field to emerge. The secondary pair of polarities arise due to the oblique angle of incidence between the flux rope's legs and the photosphere. Highly twisted field lines in each leg thread the photosphere twice, creating two adjacent polarities (Figures 11c,d). This results in the formation of sheared field lines connecting nearby primary and secondary polarities. This primary and secondary polarity structure will be seen again in higher twist simulations and will be explained in more detail in §3.6. In the next stage, shown in Figure 11d, field lines originating closer to the inside of each leg begin to emerge. Due to the kink, these field lines in each leg are brought into contact and are able to reconnect, creating concave-up field lines at the chromospheric level. As a result of their concave-up shape, the field lines in this stage of emergence are unable to lift the heavy plasma in order to emerge from the chromosphere into the corona. This creates a dip in these field lines, seen in Figure 11d.
The photospheric signatures seen in this simulation can be compared to those reported by Takasao et al. (2015), who emerge a flux rope with the same value of ζ, but with a longer initial density perturbation, λ = 2.55 Mm, versus λ = 1.2 Mm used here. In particular, by comparing Figure 11d here with Figure 7 of Takasao et al. (2015), one can see that the secondary polarities that form in our simulation as a result of the oblique angle of incidence of the flux rope's legs onto the photosphere are also present in the simulation of Takasao et al. (2015). These features are visible both early in their emergence (t = 320τ ) and later (t = 340τ ) as a diffuse component of the normal magnetic field separating the two primary polarities. Although Takasao et al. (2015) do not discuss these features, we believe they are likely also formed as a result of the oblique angle of incidence of the legs of the flux rope onto the photosphere.
On the other hand, the simulation reported by Takasao et al. (2015) displays, at t = 330τ and t = 340τ , two long, thin polarities that are pressed up against each other between the two primary polarities. Our photosphere does not, however, show any evidence of such long, thin polarities that are pressed up against each other. Takasao et al. (2015) find that these thin polarities form as a result of plasma downflows submerging previously emerged flux. We looked for evidence of such downflows during the middle stages of our emergence but did not see any.
During the latter phases of the emergence (Figures 10c,d), the primary polarities elongate, reminiscent of the behavior of the so-called 'magnetic tongues' (e.g., Luoni et al. 2011;Poisson et al. 2015Poisson et al. , 2016, which originate from the azimuthal magnetic field component of the emerging flux rope. The horizontal magnetic field becomes strongly sheared along the polarity inversion lines separating the primary and secondary polarities ( Figure  10d,e), since this field is created by short field lines encircling each leg of the flux rope. The final configuration of the field after 25 min of emergence is a set of overlying loops that spread out laterally and up, as well as a lower lying flux rope onto which is collecting plasma in its concave up dip, as can be seen by the v x shading in Figure 11e. An interesting feature of this emergence simulation is its remarkable symmetry. Transforming the photospheric normal magnetic field by i.e, rotating it by 180 • and inverting its sign, leaves the photospheric distribution nearly unchanged for the duration of the simulation. This is not the case in the previ-ous two simulations, where the emerging flux ropes lost their symmetry as they were broken up during their rise through the convection zone. The ζ = 1.5 flux rope, on the other hand, stayed coherent throughout its rise, evidently due to the strong azimuthal field supporting the flux rope against perturbations. Since the kink instability is inherently symmetric, the flux rope and, as a result, the photospheric distribution, maintained its invariance for the most part under this transformation. Slight deviations from perfect symmetry may be due to the onset of reconnection inside the flux rope.
3.4. ζ = 2 3.4.1. Subsurface Evolution We now present the emergence of a flux rope with ζ = 2. In Figure 12, we show snapshots of its rise through the convection zone at several times. Although the initial density perturbation generates a rising Ω-loop, this behavior is extremely short lived. Already by t ≈ 1000 s, the Ω-loop has started to rotate out of its initial plane.
The Ω-loop becomes aligned to the z-axis by t ≈ 1605 s, and striations only become visible very late in the flux rope's rise through the convection zone.

Photospheric Evolution
The flux rope penetrates the photosphere at t ≈ 1745 s. Like the ζ = 1.5 emergence, the two polarities formed by the emerging flux rope develop into elongated, crescentshaped polarities (Figure 13a). Figures 14a,b show the field lines that emerge first, with the B y shading demonstrating the shape and almost complete revolution of the flux rope axis by the time it reaches the photosphere. These first two polarities elongate and make a circular shell around two more sets of polarities that subsequently emerge around t = 2225 s (blue/yellow arrows in Figure  13b), forming a four-leaf clover pattern that is qualitatively similar to the pattern observed in the ζ = 1.5 evolution ( Figure 10b). Once again, the primary polarities are connected by sigmoidal field lines oriented along the z-direction, and the secondary polarities are connected to the primary polarities by short, sheared field lines encircling each leg of the emerging flux rope (Figure 14c), with the atmospheric portion of the field line appearing sheared and the subsurface part appearing twisted. Both sets of polarities approach each other and elongate (Figure 13c), forming short, quasi-potential loops and longer sheared arcades in the corona (Figure 14d). Each polarity then rotates about itself, into the configuration shown in Figure 13d, where the primary and secondary polarities contain horizontal field that is highly sheared along the polarity inversion lines (Figure 13e). This stage of the photospheric evolution, with two primary polarities separated by two secondary, elongated polarities, is also qualitatively similar to the photospheric evolution seen in the ζ = 1.5 case (Figure 10c). The coronal field lines at the end of the simulation are much more complicated than in any of the previous simulations. They are extremely compact, and low-lying field lines appear to be interacting quite strongly with overlying field lines (Figure 14e,f). In this case, in spite of the complexity of the photospheric field, the transformation in Equation 3.3 leaves the photospheric distribution almost entirely unchanged for a large portion of the simulation.
3.5. ζ = 4 3.5.1. Subsurface Evolution Here we present the emergence of a flux rope with ζ = 4. To our knowledge, this is the most highly twisted flux rope that has been emerged in simulations to date. The rise of this flux rope through the convection zone is shown in Figure 15. The rising Ω-loop generated by the mass density deficit is very quickly affected by the onset of the kink instability. As early as t ≈ 505 s, the Ω-loop has rotated out of its initial plane. By t ≈ 705 s, the apex has formed a complicated knot-like structure, and certainly no longer resembles an Ω-loop. As the apex approaches the photosphere, the structure becomes completely aligned with the z-axis. Just before the flux rope emerges, at t ≈ 1205 s, the apex of the rising portion of the flux rope resembles a knot, and appears as an Sshaped structure when looking down on it. The peak of the rising flux rope is highly compact, and very different in structure than the kinks seen in the previous two simulations.

Photospheric Evolution
This flux rope reaches the photosphere at t ≈ 1315 s, where it emerges in the shape of a figure-8 (Figure 16a). Figures 17a,b show the field lines at and below the photosphere at this instant. There is clearly an extremely compact, highly tangled structure emerging into the corona. There is also some evidence that a double-overhand knot has formed from a succession of reconnection events during the rise of the flux rope, which may be a general result for an extreme kink (Linton et al. 1998). This is shown in Figure 17c, where a single red field line is shown looping around on itself in a way that cannot be untangled without cutting the ends. In Figure 17d, this field line is shown wrapped around a single yellow field line that wraps itself into an extremely compact bundle. The next stage of the evolution, seen at t = 1565 s in Figures 16b  and 17e shows many of the same features as the ζ = 1.5 and ζ = 2 simulations: primary polarities connected to each other via sigmoidal coronal field lines and connected to elongated secondary polarities by short, sheared field lines encircling the two offset legs of the emerging flux ropes. In this case, however, the field lines are much more complicated, since the compactness of the emerging flux rope results in long, wandering field lines in the corona. Nevertheless, even the third stage of the evolution, shown in Figures 16c and 17f, at t = 1895 s, shows two primary polarities and two elongated, secondary polarities (indicated by the blue and yellow arrows, respectively, in Figure 16c). This time, however, there is also another set of two elongated polarities encircling the central set. Numerous polarity inversion lines exist at the photosphere at this stage and persist until the end of the simulation, with strongly sheared horizontal fields between the primary and secondary polarities ( Figure  16d,e). The field lines in the corona are extremely tangled, and no coherent structures can be discerned (Figures 17g,h). While some regions on the photosphere remain wholly unchanged by the symmetry transformation Equation 3.3, some of the dispersed flux in the interior of the active region is not invariant under this transformation, likely due to copious reconnection reorganizing the field. Additionally, the polarities seen in Figure 16d near (y, z) = (8.5, 0) are due to the emergence of secondary buoyant sections along the rising flux rope. These sections are produced by the strong non-linearity of the kink mode in such a highly twisted flux rope, where it produces kinks along the entire axis of the flux rope (these can be seen in Figures 15g,h).

The Role of Twist in Flux Emergence
For flux ropes with ζ ≤ 1, the subsurface behavior is determined, to a large extent, by the initial density perturbation. An Ω-loop is formed due to the Gaussian shape of the density perturbation (Eq. 2.13), and this Ω-loop rises, mostly undisturbed, through the convection zone. As the flux ropes continue to rise, and the helical field lines emerge into the corona, the legs become slanted toward the photosphere, hitting it at an angle. This behavior turns out to have important implications for the signatures observed at the photosphere. When the flux rope reaches the photosphere, it forms two distinct polarities resembling simple, bipolar active regions which gently rotate and separate as the flux rope emerges, with dispersed, weak field appearing between the two polarities. This evolution is typical of many previous simulations (Manchester et al. 2004;Murray et al. 2006;Archontis et al. 2009;Leake et al. 2013).
The behavior of flux ropes with ζ > 1, however, is determined less by the initial density perturbation, and more by the properties of the kink instability. In the convection zone, these flux ropes undergo a deformation of their central axis, whereby the apex of the Ω-loop rotates to form an -shaped loop, and the legs of the rising structure become offset from the plane formed by the initial Ω-loop. Like the weakly twisted flux ropes, the legs of the highly twisted flux ropes become slanted, and hit the photosphere at an angle. A common feature of the photospheric evolution of the ζ ≥ 1 flux ropes was their apparent four-fold structure, with two strong, primary polarities separated by two elongated, secondary polarities of alternating sign. Sigmoidal field lines connected the two primary polarities, while the secondary polarities were connected to their adjacent primary polarities with short, highly sheared field lines. The sigmoidal field lines connecting the primary polarities were due to the connection between the two off-center legs of the flux rope, while the field lines connecting the secondary to the primary polarities were due to the field lines encircling each individual leg of the flux rope. An obvious question is: why was this behavior not observed in the lowest twist simulation with ζ = 0.5, despite the legs of that flux rope being incident on the photosphere at an oblique angle as well?
A simple explanation is that these signatures come from the photosphere cutting through a twisted flux rope at an oblique angle with respect to its axis. Figure 18 illustrates the concept. Two twisted flux ropes, with magnetic field structures defined by Equation 2.9 cut through the photospheric plane as shown in Figure 18, panels a and b. The axial magnetic field of the two flux ropes is specified to be oppositely directed, representing the two legs of an emerging Ω-loop flux rope, and the flux ropes are incident on the photosphere at an oblique angle, taken to be 45 • in this example. The direction of field lines through the photospheric plane depends on both their twist and angle of incidence, relative to the normal to the photosphere. For sufficiently strongly twisted field lines, or for large angles of incidence, field lines will go through the photospheric plane not once but twice (or even more), in opposite directions on either side of the flux rope axis, producing oppositely signed vertical magnetic field. This can be seen in the zoom-in in panel b, where the red arrows point to example field lines that go through the photospheric plane twice. Figures 18cg show the photospheric vertical magnetic field due to the two flux ropes shown in Figures 18a and b for different values of the twist. As the twist of each flux rope increases, a second polarity appears in the vertical magnetic field distribution of each flux rope.
Quantitatively, this reversal in the photospheric vertical field, B x , is due to the fact that both the axial and azimuthal components of the flux rope have a vertical component when projected onto the photosphere (for a non-zero angle of incidence of the flux rope axis). On one side of the axis of each flux rope, the projection of the azimuthal field into the vertical direction will be in the same direction as the projection of the axial field into the vertical direction and will add constructively to increase B x . On the other side of the axis of each flux rope, the azimuthal field projection will be in the opposite direction of the axial field projection and will add destructively, reducing B x . If the azimuthal field projection is strong enough or the angle of incidence is sufficiently oblique, the azimuthal field projection will dominate the axial field projection on one side of each flux rope's axis, and will produce a photospheric signature with different signs of B x on either side of each flux rope's axis (assuming, of course, that the flux rope radius extends far enough). As a result, the photospheric distribution of B x will consist of two oppositely signed polarities for each flux rope leg, for a total of two pairs of polarities. This is precisely what is observed in the high twist simulations. Each flux rope leg, in addition to having a strong azimuthal field, rises through the photosphere at an oblique angle, resulting in a dominant, primary polarity, and a smaller, secondary polarity. In the lowest twist flux rope, the twist and angle of incidence together were insufficient to create the four-fold structure observed for the high twist simulations.
A lower limit for the magnitude of the flux rope twist required for seeing the four-fold structure can be calculated mathematically by projecting a flux rope oriented in some arbitrary direction onto the photospheric plane. This is done as follows. A constant twist flux rope, such as the one given in Equation 2.9, oriented along some arbitrary axisξ has a magnetic field given by for the cutoff radius R with axial and azimuthal components where ρ is the radial coordinate, ξ is the direction along the flux rope and ψ is the polar coordinate around the flux rope perimeter. If the flux rope intersects the photosphere at an angle θ with respect to the vertical x-axis, i.e.,x ·ξ = cos θ, x ·ψ = ± sin θ, (3.6) the axial and azimuthal fields of the flux rope produce the normal photospheric field B x (y, z) given by (3.7) Here the ± sign represents the different direction at which the azimuthal field B ψ intersects the x-axis on either side of the flux rope axis. On one side of the flux rope axis,x ·ψ > 0, the projection of the azimuthal field B ψ adds to the projection of the axial field B ξ , so B x will be positive/negative if B ξ is positive/negative. On the other side of the flux rope axis,x ·ψ < 0, the azimuthal field B ψ projects in the opposite x-direction than the axial field B ξ , and in this case, the projection of the magnetic field will change sign, from positive/negative to negative/positive, if ρ > a ζ cot |θ| (3.8) The smallest twist for which B x will change sign at the photosphere at ρ = a is ζ = cot |θ|. (3.9) For the flux ropes shown in Figure 18, θ = ±45 • , which implies that the normal component of the photospheric field will change sign at ρ = a for an inclined flux rope with ζ = 1. Note, however, that the flux rope extends out beyond ρ = a, up to ρ = R, with R = 6a. Thus a reversal of B x will still occur outside ρ = a, as can be seen in Figure 18d. Even the lowest twist case shown in Figure  18c displays a weak reversal at ρ = 2a, though it is not visible with the color scale in the Figure. In these cases, the in-plane horizontal vector field is primarily parallel to the polarity inversion line separating the secondary polarity from the adjacent primary polarity, exactly as observed in the high twist simulations presented above. Interestingly, it is aligned parallel to the polarity inversion line, but it is mostly concentrated off the polarity inversion line, except for ζ = 2. As ζ increases, the azimuthal field increasingly overwhelms the axial field on one side of the flux rope axis, creating the reversals in B x evident in Figure 18d-f. Thus, the primary and secondary polarities are a direct consequence of flux ropes with large values of ζ being incident on the photosphere at an oblique angle. The mathematical arguments presented above explain why these features of the photospheric field do not appear for the lowest twist flux rope and why the photospheric signatures in the higher twist simulations display at least four, rather than two, strong polarities; two for each leg of the flux rope.
Of course, the above arguments do not take into account distortions of the flux ropes that occur during their rise through the convection zone. As the flux ropes rise through the convection zone and untwist themselves in the corona, their profile deviates from the simple form given in Equation 3.4, resulting in photospheric features that are more complicated than the simple picture presented in Figure 18, but the signatures observed at the photosphere are well represented by simply having two oppositely directed flux ropes oriented at angle to the photosphere. Figures 19-22 show the evolution of various photospheric parameters for each of the simulations as a function of time since emergence, defined as the time at which a single photospheric pixel first exceeds B thresh x = 30 G. For each of the quantities described below, we integrate over the entire photospheric plane, excluding pixels where |B x | < B thresh x . Hence, even regions far from the primary polarities, such as the secondary polarities in Figure 16d located at (y, z) = (8.5, 0), are included in the calculation.

Simulation Comparison
We write the unsigned flux at the photosphere as (3.10) and plot Φ ph scaled by the initial axial flux in the flux rope We calculate the average vertical current density, signed and unsigned, taken from the x-component of Equation 2.11,J x . We calculate the average vertical twist parameterᾱ is a commonly used proxy for the twist at the photospheric level (Pevtsov et al. 1994). For completeness, we also calculate the average twist parameter, (3.14) despite the fact that is not measurable in observations. We calculate the distance between the flux weighted centers of mass L com = |L com | (cf. Equation 3.2) and the angle between the flux weighted centers of mass and the initial flux rope axis Finally,we calculate the current neutralization ratio (Török et al. 2014;Vemareddy et al. 2015;Liu et al. 2017) where I ret and I dir are the return and direct currents, respectively, and are defined as (3.18) R has been suggested as an important parameter in determining flaring behavior (Dalmasse et al. 2015;Kontogiannis et al. 2017;Liu et al. 2017). We also investigate the energetics of each simulation by defining the magnetic energy in the entire simulation domain at the start of each simulation as (3.19) and the magnetic, potential and free energies in the coronal volume (x > 0), respectively, as follows: We calculate the scaled parameters (3.23) In the above expressions, B p is the potential magnetic field in the corona, which satisfies the equation The separation distance between the flux weighted centers of mass is one of the defining characteristics known to differentiate δ-spots from other sunspot group types (Künzel 1965). The evolution of the separation distance between the flux weighted centers of mass in our simulations, shown with the red dots in Figures 4b, 7b, 10c, 13c, and 16c, displays clearly different behavior between the initially weakly and highly twisted flux ropes, as seen in the top panel of Figure 19. The more weakly twisted flux ropes undergo significant separation, in one case more than double that of the highly twisted flux ropes. In the weakly twisted cases, there is clearly a general trend by the opposite polarities to move apart for the majority of the emergence. The kink-unstable cases, in contrast, have a peak separation distance of only 2 − 2.5 Mm, after which the separation does not increase, and may even decrease. Thus, the flux ropes with higher initial twist display, at later times in the emergence process, a much more compact photospheric flux distribution, in agreement with the observed behavior of δ-spots. Based on this plot, there appears to be a bimodal distribution of separation distances between highly and weakly twisted flux ropes. It will be interesting to check this with observations.
Rotation of the active region is another photospheric observable that is one of the defining characteristics of δ-spots (Kurokawa 1987). The angle θ com between the line connecting the flux weighted centers of mass and the initial flux rope axis changes quite a bit for each simulation, as shown in the bottom panel of Figure 19. The lower twist flux ropes emerge with θ com ∼ 90 • , corresponding to the fact that the azimuthal flux encircling the axial flux emerges first. The polarities then rotate back around to ∼ 20 • , corresponding to the emergence of almost purely axial flux in the legs of the flux rope. The kink-unstable flux ropes, meanwhile, emerge with θ com greater than around 180 • , before rotating to a value near 70 • (for ζ = 1.5 and ζ = 2) and ∼ 15 • for ζ = 4. The total rotation undergone by the emerging flux ropes increases with increasing initial twist. If θ com < 90 • is assumed to be the angle of emergence that obeys Hale's law, as it should be for a typical weakly twisted emerging flux rope, then the three highest twist flux ropes obviously violate Hale's law, in agreement with the observed behavior of δ-spots, and the two lowest twist flux ropes nearly violate Hale's law. In both cases, the active region rotates into a configuration that satisfies Hale's law. It should be noted, however, that this calculation is obviously affected by the presence of the secondary polarities created by the projection of the twisted flux onto the photosphere, possibly disguising an ever larger contrast between weakly and strongly twisted flux ropes.
One of the simplest quantities to measure observationally at the photospheric level, the unsigned magnetic flux, plotted as a function of time for each simulation in the top panel of Figure 20, is commonly measured to characterize an emerging active region (e.g., Sun & Norton 2017). In our simulations, although each flux rope has, initially, the same axial flux, the total flux at the photospheric level is determined by both the axial and azimuthal fluxes. As a result, Φ sc is largest for ζ = 4, and is smallest for ζ = 0.5 and ζ = 1. Interestingly, between about 4-10 minutes after emergence (emergence is defined here as the time when B x first exceeds 30 G on the photosphere), or about 25% of the total simulation time, the ζ = 1.5 case has less flux going through the photosphere than ζ = 1. In general, however, after equal emergence times, we find more flux at the photospheric level for higher initial flux rope twist.
In the bottom panel of Figure 20, we plot the flux emergence rate against the emerged flux, following the method of Norton et al. (2017). Here, the emerged flux is defined as with Φ 2 = 0.9Φ max and Φ 1 = 0.1Φ max , where Φ max is the peak emerged flux. Thus Φ em is the amount of flux that has emerged during the interval ∆t = t 2 − t 1 . The flux emergence rate is defined aṡ where ∆t = t 2 − t 1 is the interval of time during which the flux Φ em is emerged. Plotted this way, the flux emergence rate increases linearly with increasing flux (initial twist). The fluxes observed here (10 19−20 Mx) are within the observed range, as are the emergence rates (10 20−21 Mx hr −1 ). However, the rate of emergence for a given flux is an order of magnitude higher than observed by Norton et al. (2017), though in line with simulated results (Cheung et al. 2007(Cheung et al. , 2008. The magnitude of the vertical current density, often used as a proxy for energetic processes such as magnetic reconnection (Knizhnik et al. 2018), is a fundamental quantity that can be determined directly through observations. Spikes in current density have been shown to coincide with X-class flares in δ-spots (Vemareddy et al. 2015). In our simulations, the average unsigned vertical current densityJ x strongly oscillates during each simulation, as shown in Figure 21a. There does not seem to be any dependence ofJ x on the initial flux rope twist. At certain times during the emergence, the ζ = 0.5 or the ζ = 1.5 cases have higher values ofJ x than the other cases, making it difficult to say anything definitive about the emerging structure from these quantities.J x is different, in general, than the magnitude of the full vector current density, Equation 2.11, but the latter cannot be calculated from photospheric vector magnetograms. In our simulations, however, J is easily measurable, and we find qualitatively and quantitatively very similar behavior ofJ x andJ (not shown here). We also measure the average signed vertical current density |J x | for each simulation, shown in Figure 21b. There is some evidence that flux ropes with higher twist have larger average signed current densities than do flux ropes with smaller twist, but this breaks down around t ≈ 24 min into the emergence.
The current neutralization ratio, easily measurable at the photosphere, has been found to be further from unity for flaring active regions than for non-flaring active regions (Liu et al. 2017). R = 1 corresponds to a neutralized current at the photosphere, i.e., the direct current is balanced by the return current. Given the propensity of δ-spots for flaring behavior, it is natural to ask whether this deviation from unity is observed in our simulations. The bottom panel of Figure 21 shows that the currents in our simulation are quite non-neutralized, and R is not obviously dependent on the flux rope's initial twist. There seems to be a general trend of rapid decrease followed by gradual increase of the value of R for all of the flux ropes, but the specific ordering does not seem to depend on ζ. During some phases of emergence, ζ = 2 has higher values of R than other simulations, with ζ = 0.5 or ζ = 4 having the lowest values, and at other times ζ = 4 has the highest value and ζ = 0.5 or ζ = 1.5 has the lowest value. In all cases, however, currents do not seem to be neutralized at the photospheric level. α x has traditionally been the parameter used to study active region twist (Pevtsov et al. 1994(Pevtsov et al. , 1995Leka et al. 1996;Holder et al. 2004) since it can be calculated di-rectly from observations using photospheric vector magnetograms. However, in our simulations, the evolution ofᾱ x , shown in Figure 22a, does not display any obvious twist-dependent behavior. It is not generally the case that a largerᾱ x corresponds to a larger value of ζ. In fact, towards the end of the simulation, the ζ = 0.5 case has the largest value ofᾱ x . The ζ = 1 and ζ = 4 flux ropes have the smallest values ofᾱ x throughout a significant fraction of the emergence process. Although each flux rope continuously emerges, the value ofᾱ x does not continuously increase, instead showing peaks and troughs during the course of each emergence. This would seem to indicate that,ᾱ x does not, by itself, reveal anything about the structure of the emerging flux rope. Although α x is the quantity that can be measured from photospheric magnetogram, the twist parameter is actually α as defined in Equation 3.14. This quantity can be measured directly in our simulations, though not in observations. The time evolution ofᾱ is plotted in Figure 22b. Notably, the plots ofᾱ do not resemble the plots ofᾱ x , except that they share a lack of a clear dependence on twist. This indicates thatᾱ x may not be a good proxy forᾱ at the non-force free photosphere.
From the figures, it appears that little, if anything, can be ascertained from the temporal evolution of either α x orᾱ, since the photosphere is not force free, and J × B = 0 is expected to be a poor approximation. However, there are certain locations where the photosphere can be taken to be force free. These locations can be determined as follows. The current density and magnetic field are related via |J × B| = |J||B| sin θ, (3.27) (3.28) with θ the angle between J and B. Therefore and Thus, locations on the photosphere where Σ 1 (Σ 2 ) is large (small) correspond to locations where the magnetic field is close to force free. Therefore, we define the 'force free twist parameters' as (3.32) In panels c,d of Figure 22, we plot the average unsigned vertical and full force free twist parameters for each simulation, having chosen Σ thresh 1 = 0.8. Here, in contrast to panels a) and b), one can see that there is a discernable correlation between theᾱ x andᾱ. However, there is still no discernable correlation between ζ and the magnitude of eitherᾱ x orᾱ. Especially early in the emergence process, the ζ = 1.5 and ζ = 2 flux ropes have the largest values, while the ζ = 0.5 and ζ = 1 flux ropes have the smallest values. The ζ = 4 flux rope seems to fall somewhere in the middle, destroying any dependence of these quantities on initial twist. This results in the average photospheric α, for all values of ζ, clustering around a narrow range, particularly later in the emergence. For example,ᾱ x clusters within the range 10−30×10 −8 m −1 for these simulations. The fact that this is the observed result from flux ropes whose original peak alpha ranged from 330 − 2700 × 10 −8 m −1 illustrates how different the photospheric α can be from the initial, peak α of a convection zone flux rope.
In Figure 23, we plot α x , α, α x,f f , and α f f at a representative snapshot of the ζ = 1 simulation. Immediately obvious from these figures is the complexity of the signal in α x and α. No single sign of either dominates. Meanwhile, the sign of the photospheric distribution of α x,f f and α f f is overwhelmingly positive, reflecting the fact that α in the initial flux rope was positive. It is evident, from looking at these figures, that a simple average over the twist parameter, as is done in several papers (e.g. Pevtsov et al. 1994Pevtsov et al. , 1995Leka et al. 1996;Holder et al. 2004) is bound to be affected by the cancellation of positive/negative values of approximately equal magnitude. Thus, averaging over α x (or even α if it could be measured), will produce a value of minimal significance. The Figure also address the question of whether α x is a reasonable proxy for α. It is clear that α x is much noisier than α and dominates near polarity inversion lines outside the main polarities, whereas α is relatively smooth in the interior of the active region (near y = z = 0). Comparing α x and α reveals that they do not match well, especially in the interior of the active regions. On the other hand, comparison of α x,f f and α f f reveals significantly more similarity, indicating that α x,f f is a good proxy for α f f , supporting the expectation that α x should be a good proxy for α in force-free fields. This calls into question the usefulness of measurements of α x at the photosphere as a proxy for α, since identifying which regions are force-free in observations is, at present, a challenge. This result also justifies, a posteriori, our choices of ζ which are significantly larger than the values required to explain observations (Pevtsov et al. 1995;Longcope et al. 1999). Since there appears to be no relationship between the values of either α x or α observed at the photosphere and the initial flux rope twist, the twist inferred from α using theoretical arguments may not necessarily be representative of its true value. Each of the highly twisted flux ropes initially had a much larger α, but the conversion of twist into writhe as a result of the kink instability decreased α to values very similar to those measured for the lower twist flux ropes. Furthermore, the salt-andpepper nature of the distribution of photospheric α and α x means that any strong positive values may be canceled by strong negative values, decreasing the overall average. Most of the quantities described above are observationally measurable from photospheric data alone. However, some of the most important quantities, from the standpoint of energy release, are the magnetic, potential, and free energies in the corona. Although the magnetic and free energies rely on knowledge of the magnetic field in the entire corona, rather than just the photosphere, they are directly related to how much energy can be released in a flaring event. To understand how these quantities scale with initial flux rope twist, we plot the magnetic and potential energies in Figure 24. Neither the magnetic nor the potential energies show a tendency to be larger for large initial twist and smaller for small initial twist. While the ζ = 4 simulation has the largest magnetic energy for the duration of the simulation, the ζ = 0.5 flux rope has the second most, and the ζ = 1.5 flux rope has the least throughout the simulation. The ζ = 1.5 and ζ = 2 flux ropes have somewhat less free energy than the ζ = 0.5 flux rope. Both the ζ = 0.5 and ζ = 4 flux ropes have approximately constant free energy starting around 10 min after the initial emergence, while the ζ = 1.5 flux rope has approximately constant free energy starting around 15 min after the initial emergence. The ζ = 1 and ζ = 2 flux ropes continue to increase their free energy through their emergence. The free energies measured here are quite large: about 80% of the total magnetic energy is in the form of free energy. Meanwhile, the nonlinear force-free extrapolations of DeRosa et al. (2009) of a stable AR magnetogram found that approximately 20% of the total magnetic energy in the region was in the form of free magnetic energy. This large difference in free-energy content could arise from the assumptions inherent in the nonlinear force-free extrapolations, which use manifestly non force-free photospheric magnetograms as boundary conditions, or because the AR studied by DeRosa et al. (2009) did not contain a δ-spot and was not flare active during the observation, and even later did not produce anything larger than a C8 flare. It is possible that a similar analysis conducted on a δ-spot containing AR would find that a much larger fraction of the total magnetic energy was stored as free energy.

DISCUSSION
The lack of understanding of the formation of δ-spots is one of the most important impediments preventing proper prediction of solar energetic events. The emergence of highly twisted, kink unstable flux ropes has been hypothesized to be the source of a large percentage of δ-spots Linton et al. 1996;. In this work, we test this model using numerical simulations of the emergence of both initially kink-stable and -unstable flux ropes from the convection zone, through the photosphere, into the corona. We focused mainly on AR properties that can be measured using line-of-sight and vector magnetogram observations and related them to the subsurface properties of emerging flux ropes.
By calculating the flux weighted centers of mass of each polarity, we demonstrated that initially highly twisted flux ropes remain much more compact, and rotate much more than, their lower twist, kink-stable counterparts. We also showed that while low twist flux ropes obey Hale's law, kink unstable flux ropes appear to violate it. Our results demonstrate that at the photosphere, emerging kinked flux ropes behave in a manner that is both qualitatively and quantitatively similar to the most welldocumented behavior of δ-spots (Künzel 1965;Smith & Howard 1968;Kurokawa 1987;Zirin & Liggett 1987;Tanaka 1991;López Fuentes et al. 2000, 2003. However, quantitative statements of exactly what is meant by 'compactness' and 'rotation' are sparse, so it is diffi-cult to compare with observations. On the other hand, our work also indicates that quantitative measurements of many δ-spot properties are unlikely to produce results which are vastly different than those measured in simple sunspot groups. Parameters like α x , J x , R, all of which are measurable from photospheric observations, without recourse to models of the coronal field, do not reveal major consistent differences for highly twisted versus weakly twisted flux ropes. Photospheric measurements of these parameters may reveal higher values for δ-spots, but it is also possible that they may reveal lower values, depending on the instant at which the measurement was taken. Plotting the time evolution of these parameters may also reveal differences, but the results presented here suggest that, with the exceptions of the separation distance between opposite polarities and rotation angle of the opposite polarities, there is no simple dependence of many measurables on the flux rope's initial twist. Even the free energy in the coronal volume, thought to be a proxy for energetic flaring events, does not show an obvious dependence on the initial flux rope twist, despite the much more complicated coronal topology evident in the higher twist simulations. In contrast, Figure 8 in Toriumi et al. (2011) shows that, for initially weakly twisted flux ropes, there appears to be a dependence of magnetic energy on initial twist. However, the differences in the numerical setups between those simulations and the ones presented here make a direct comparison challenging.
The lack of dependence of our results on the initial twist arises due to several effects. First, the nature of the kink instability is to reduce the field line twist by converting twist to write (Berger 1984;Linton et al. 1996), so that once a flux rope has kinked, its effective twist, and therefore α, is smaller than its initial value, and may turn out to be comparable to that in a lower twist flux rope. Indeed, it is possible that the lower twist flux rope need not, itself, be kink stable for this result to hold. Since it started out with a smaller value of α, the lower twist flux rope may, as a result of a kink, experience a smaller decrease in α than the higher twist flux rope, resulting in the two flux ropes having comparable values of α. Regardless, it is clear that α is not simply transferred directly from the convection zone to the photosphere, but instead evolves dynamically during the emergence process, explaining its decrease from convection zone to photospheric values. Second, the rapid expansion of the magnetic field once it reaches the coronal level causes a significant decrease in the value of α at the photospheric level, since the expansion of the field into the corona occurs faster than Alfvén waves can spread the twist along the field line, creating a gradient in α along the expanding field line. This effect has also been observed in previous high twist simulations (Fan 2009;Leake et al. 2013). Finally, reconnection inside the flux rope itself could destroy a lot of the initial twist structure, perhaps through an inverse cascade of magnetic helicity (Antiochos 2013;Knizhnik et al. 2015Knizhnik et al. , 2017a, leaving a relatively weakly twisted internal structure, forming active regions that have lost much of the information about the structure of the initial flux rope. The simulations presented here show that the dynamics of emerging highly twisted flux ropes cause copious internal reconnection to occur, as evidenced by the for-mation of concave up loops and knots during the emergence process. Such complicated topological structures can only be formed, in the absence of an external magnetic field, by internal reconnection. This internal reconnection would likely manifest itself as flaring behavior on the Sun. We conclude that flaring internal to an emerging kinking active region can be strong, independent of the state (or even absence) of external fields around it, in agreement with observations showing that δ-spots can produce X-class flares within themselves (Zirin & Tanaka 1973;Zirin & Liggett 1987;Wang et al. 1991;Schmieder et al. 1994  The mass density (blue), pressure (orange) and temperature (red) are scaled by their photospheric (x = 0) values. Shown in black is the magnetic pressure profile for the ζ = 1 flux rope. The convection zone is located at x/Hp < 0, the temperature minimum region at 0 < x/Hp < 12, and the corona begins at x/Hp = 25, being separated from the temperature minimum portion of the atmosphere by a short transition region.      . Field lines, overplotted on photospheric magnetograms, at various stages of the ζ = 1.5 simulation. In a) and b) field lines, colored by By, are seen from the side and above, respectively. In e) field lines are shown at the final state of the simulation, colored by vx. Figure 17. Field lines, overplotted on photospheric magnetograms, at various stages of the ζ = 4 simulation. In a) and b) field lines, colored by By, are seen from the side and above, respectively. In c), a single knotted field line is shown, with several other field lines shown wrapped up in the knotted field line in d). In g) field lines are shown at the final state of the simulation, colored by vx. Figure 18. a) Two flux ropes, with a peak field strength of B 0 = 6.5 kG, oriented at ±45 • to the vertical. b) Zoom in of a), with red arrows denoting example field lines which cross the photosphere twice. Bx (greyscale) and B h (colored vectors) due to two flux ropes oriented at a 45 • angle to the vertical with c) ζ = 0.5, d) ζ = 1, e) ζ = 1.5, f) ζ = 2, g) ζ = 4. The color table applies to panels c-g.