Results and plans for nuclear spectroscopy of superheavy nuclei: the Lund perspective

A survey of nuclear structure results obtained along α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document}-decay chains stemming from isotopes of element 114, Fl, and element 115, Mc, is presented. Measures taken to ensure a reliable hand-over of nuclear spectroscopy output to nuclear structure theory are outlined. Possible links between decay chains starting from element 117, Ts, and element 115, Mc, are commented on. The experiments were conducted with the TASISpec set-up placed in the focal plane of the TASCA separator at the GSI Helmholtzzentrum für Schwerionenforschung GmbH in Darmstadt, Germany. Plans to improve the decay-spectroscopy set-up towards future superheavy endeavours are described.


Introduction
In the wake of the discovery of new superheavy elements, nuclear spectroscopy experiments aim at providing benchmarks for nuclear structure theory at the uppermost end of the nuclear chart. A central and long-standing question concerns nuclear binding or nuclear stability of superheavy nuclei. This relates to the underlying nuclear shell structure, i.e., the evolution of the nuclear shell structure as a function of, for example, proton and neutron number. In other words: What is the next (spherical) doubly-magic nucleus beyond 208 Pb?
By 1969 elements up to Z = 104 had been discovered experimentally. Calculations by, for instance, nuclear theorist Sven Gösta Nilsson and collaborators in Lund described a superheavy 'Island of Stability' around nuclei with proton number Z = 114 and neutron number N = 184 [1]. Calculations by Meyers and Swiatecki [2] as well as Mosel and Greiner [3] had provided similar predictions. By the year 2000 elements up to Z = 116 and N = 176 had been discovered in experiments and a new class of nuclear models based on nuclear density-functional theory started to make predica e-mail: Dirk.Rudolph@nuclear.lu.se (corresponding author) tions for how heavy elements that can exist [4,5]. This new class of models has indicated the importance of density redistribution effects that can move the energy gaps in the nuclear spectra. One such effect is that the electrostatic repulsion between the protons can form a region of reduced density in the centre of the nucleus that can lead to decisive changes in the spectra. Owing to this interplay of subtle effects, a more nuanced picture has emerged where models that make similar predictions for lighter nuclei differ in the predictions of heavy or superheavy nuclei (see, e.g., [6,7] and references therein). Further, the inclusion of beyond mean-field methods can probe deformation softness or shape coexistence in the full β-γ plane [8][9][10].
As a consequence, to avoid extrapolations over some 30 or so mass units, the experimental task is to provide theory with spectroscopic data to probe their models or model parametrisations at or near element Z = 114, flerovium. Interesting quantities are precise Q α values and the corresponding decay times. Another crucial piece of data concerns γ and X-ray radiation that can be emitted following an α-decay. Since the X-ray energies are characteristic for the different elements [11] they can provide confirmation of the proton number of the nucleus being studied and, hence, the chemical element [12]. The γ rays on the other hand give access to spectral properties of the nuclei and inform us about their internal nuclear structure. For nuclei with odd particle numbers it is essential to evaluate the α decays to different excited states of the daughter nucleus in order to determine if the daughter will be left in an excited state. It is for such cases that one has a chance of observing the subsequent electromagnetic decay providing us with in-depth information directly on the heaviest elements. Figuratively speaking, the idea is to convert the discovery-related α-decay squares on the chart of nuclei (see, e.g., [13][14][15] and references therein) into decay schemes which carry crucial information to improve the description of the strong force in the heaviest nuclear systems.
The main experimental challenge is posed by the tiny production and thus observation rates. Therefore, it is of utmost importance to have a spectroscopy set-up capable of comprehensively measuring α decays, γ rays, x rays, as well as electrons from internal conversion in coincidence, with high efficiency, and with low energy thresholds. In Sect. 2, the TASISpec set-up developed in Lund is described [16], together with a reminder of spectroscopic tools used and applied to ensure high-quality spectroscopic data handed over to nuclear theory. Sections 3 and 4 review the main results of the spectroscopy experiments conducted on element 114, flerovium, and element 115, moscovium. Based on the latter, suggestions for linking decay chains associated with isotopes of moscovium and tennessine are presented in Sect. 5. This is followed in Sect. 6 by an outlook toward a novel decay station being prepared in Lund, named Lundium . A brief summary concludes this contribution.

Experimental and spectroscopic tools
The activities of the Lund group have been largely embedded in the superheavy element research programme at the GSI Helmholtzzentrum für Schwerionenforschung in Darmstadt, Germany, and more recently also in the frame of FAIR Phase-0. Both commissioning, preparatory, and the two main experiments aimed at spectroscopy of decay chains stemming from isotopes of flerovium and moscovium, respectively, were conducted in the focal plane of the gas-filled TransActinide Separator and Chemistry Apparatus (TASCA) [17,18]. The idea of a compact cube of silicon strip detectors for detection of implantation and charged-particle decays of (super)heavy nuclei, surrounded by large composite germanium detectors for highly efficient coincidence measurements with γ rays or x rays, had its origin in the so-called small-image mode of TASCA -hence the acronym TASISpec, TASCA Small-Image mode Spectroscopy. Figure 1 shows a photograph of the 'silicon cube' of TASISpec. For more information on TASISpec in general see [16]. Further TASISpec details for either of the two main experiments can be found in the PhD theses of Forsberg [20] and Såmark-Roth [21], respectively.
During TASCA and TASISpec commissioning runs, it had turned out that background in the focal plane reduced significantly more for high-transmission than small-image mode settings of TASCA. Though this implied a minor loss in measured [22] and simulated [23] TASCA transmission into TASISpec, detector rate and spectral quality were judged more important. Specific TASISpec focussing schemes, which involve asymmetric settings of the two TASCA quadrupole magnets [17], were then tested and established [20,22].
The experimental concepts of TASISpec and applied data-analysis techniques are similar to those outlined for Mounted directly at the focal plane of TASCA, superheavy ions can reach the cube through a (hidden) tube of 70 mm diameter. Implantation and charged-particle decays are registered in the DSSD seen on top of the picture. Particles emitted from the implantation spot into beam direction can be caught by four box DSSDs. Miniature 50-coaxial cables guide the detector signals to 32-channel preamplifiers [19] which are directly connected to the TASISpec chamber. For the 'element 114 experiment' in 2019/2020, a sixth DSSD was positioned right behind the implantation DSSD to veto β-decay events and other types of radioactive background. All DSSD are 60 × 60 mm 2 in size and have 32 strips on each side, i.e., a nominal pixelation of 1024. Two adjacent strips are combined for box detectors, reducing this number to 256 for each box DSSD. During experiments, the DSSD cube is surrounded by a 0.5-1.0 mm thin aluminium housing, and each DSSD is backed by a composite germanium detector of 4crystal clover and/or 7-crystal cluster type [16] other decay stations focussing on multi-coincidence highresolution decay spectroscopy in recent or this Topical Issue on heavy and superheavy nuclei (see, e.g., [24][25][26][27] and references therein). TASISpec was one of the first decay stations to be upgraded into a fully pixelated system of doublesided silicon strip detectors (DSSDs) (cf. Fig. 1). This enables tracking α-decay particles (or fission fragments) leaving the implantation DSSD and reaching one of the box DSSDs. Once dead-layer thicknesses of the two engaged pixels are known, decay energies can be reconstructed with almost as good resolution as if the α particle had moved deeper into the implantation DSSD. Relevant dead-layer scans were conducted for all TASISpec DSSD with standard α and electron radioactive sources [20,28]. TASISpec was also one of the first decay stations to employ digital electronics, i.e., to directly register the output signals of the preamplifiers in typically 70-80 μs long traces at 50-60 MHz and 12-14-bit depth. Developed at the GSI-EE Department, 16-channel FEBEX3 modules were used. Early on, this by now common way of digitization had been done solely for the 32 p-side channels of the implanta-tion DSSD, later-on for all DSSD strips, up to 256 in total, and more recently also for the germanium detector crystals, the latter with more elaborate FEBEX4 hardware and firmware. Inspecting these traces has a multitude of advantages compared with standard analog electronics. One of specific relevance for decay chains originating from flerovium and moscovium, with α-decay energies in the interval 9.0-10.5 MeV, is that false candidate events due to fast α emitters from transfer-reaction channels can be readily traced and disregarded immediately: A full-energy α detection in the implantation DSSD of some 8.5-9.5 MeV, preceded or followed by an escape-α pile-up event of some 0.5-1.0 MeV, also deposited in the implantation DSSD and within less than a few microseconds, leads to summing in conventional electronics and could thus be misinterpreted as being part of a superheavy element decay chain (cf. chap. 6 in [20]). In turn, dedicated analysis of such events can reveal new structure information for these transfer-reaction channels [29].
From a spectroscopy perspective, the availability of the preamplifier pulse allows to probe and develop codes for various digital filters in order to optimize the channel-by-channel energy resolution of the DSSDs [20,21,29]. Prior to experiments, energy calibration of the DSSD cube is done with the help of external α and electron standard sources, of course using and cross-checking pixel-by-pixel dead-layer thicknesses. For α decays with energy E α , the energy deposited by the implanted recoil has to be accounted for as well: A fraction, r , of the recoil energy contributes to the detected energy, E det -unless there is further summing with electrons from coincident internal conversion decay, of course. For the 'element 115 experiment', a constant r = 0.5 was used [20]. For the 'element 114 experiment', a more sophisticated study of a number of cases of implanted nuclei produced in preparatory experiments yielded r = r (Z 2 ) ≈ 0.40 − 0.60 for Z ≈ 85 − 115 [21,30]. Either way, the methods were tested on ample known α emitters in the range E α ≈ 7 − 12 MeV and found to agree with tabulated values within < 10 keV [20,21,29]. Finally, over the period of week-long experiments, electronic gains can drift. Therefore, data was cut into reasonably sized portions and such drifts accounted for if deemed necessary. For this purpose, the modest but notable amount of α decays of transfer-reaction background detected in the focal plane of TASCA proved useful, including events at E α = 11.66 MeV originating from 212m Po (cf. Fig. 6.9 in [20] and Fig. 3 in [30]). Simulations -or virtual experiments -using the Geant4 toolkit have entered the nuclear physics domain some time ago. Therefore, TASISpec has been coded in detail in Geant4 [31] while allowing for adoptions in geometry and detector types for a given TASISpec set-up, for instance, for the 'ele-ment 115 experiment' [32,33]. It is also possible to conduct virtual Geant4 experiments in 'real time' and compare virtual Geant4 data with data taken in real experiments. This feature has proven very helpful for experiments in which TASISpec had been placed behind ion traps, such that the decay of a certain isotope with a given mass, i.e., the decay of a selected single quantum state, could be followed over time [34] -regardless whether that be α decay, β decay, or proton emission [35][36][37]. The most important asset of the Geant4 simulations for superheavy element spectroscopy studies, however, is the capability of self-consistency control of the decay schemes derived from the measured α, electron, and photon spectra, as well as coincidences among these (see, e.g., [32,33] and below). Using the decay schemes as physics input for Geant4, virtual experiments provide expected spectra, which within given statistical limits must match the observed ones in order to validate the very same suggested decay scheme. Figures 3 in [32] or [38] provide good examples of both a calibration aspect (α-γ coincidence in 211 Bi) or discriminating E1 or M1 character of γ rays seen in 276 Mt (Z = 109) in coincidence with the α decay of 280 Rg (Z = 111). Besides this core feature of Geant4 simulations, they can also prove helpful in specific data assessments. For example, energies and depths of ions entering the implantation DSSD can be studied, in correlation with the partial energy loss of reconstructed αparticle events. This can provide additional evidence to disregard a candidate flerovium decay chain and allows to deduce a depth profile of ion implantations (cf. Fig. 5.6 in [21]).
An important facet of the nuclear spectroscopy experiments at GSI, and in particular for those in the focal plane of TASCA, is the 25% duty cycle of the UNILAC, with 5 ms beam on and 15 ms beam off. With TASISpec necessarily placed inside the TASCA experiment cave, nuclear spectroscopy with germanium detectors is hardly possible during beam-on periods due to too high radiation levels. In turn, for most of the time there is no primary beam, which allows for proper decay-spectroscopy measurements. Furthermore, detecting recoil-α correlations as they are recorded online allows to trigger an electrostatic chopper to bend the primary beam off axis in less than 20 μs. The beam-off periods can then be prolonged as long as needed to detect the remaining part of a superheavy element decay chain in a rather quiet radiation environment. It is noteworthy to point out that the 15-ms beam-off periods can help tremendously in identifying proper initial recoil-α candidates, which was most relevant for the 'element 114 experiment' due to the rather long halflives of 288,289 Fl, which are on the order of a second.
Finally, hardly any of the above mentioned spectroscopic tools or measures are particularly new or unique for TASISpec. Nevertheless, they should in one way or another be applied to any experiment claiming spectroscopic relevance for superheavy element research.

Decay chains associated with Fl isotopes
An experiment proposal targeting decay chains stemming from 289 Fl was approved at GSI in 2017. Originally scheduled for summer 2018, it could first be conducted in two parts in 2019 and 2020, respectively, due to accelerator problems. Four out of five Compex modules were delivered and used, while the set-up was otherwise a fine-tuned version of TASISpec [30]. More details of the set-up, the experiment itself, and the specific analysis steps and tools are described comprehensively in [21]. In brief, the reaction 48 Ca+ 244 Pu was focused upon, while one out of four target segments in 2019 held enriched 242 Pu instead. A total beam integral of 6.0(4) · 10 18 was collected, leading to the observation of 32 flerovium candidate chains. Of these, two were associated with 286 Fl and twelve with 288 Fl [39], as well as fifteen with N = 175 289 Fl [21,40]. Three additional decay chains of possible 289 Fl origin were disregarded because of too much similarity with random background [21,41]. Thus, 29 flerovium decay chains were observed in less than three weeks of beam time, i.e., a little more than ten per week at average beam currents of 4 · 10 12 /s at the time. Accounting for target thicknesses of 0.80(1) mg/cm 2 244 Pu, recoil transmission into TASISpec, as well as detection and identification probabilities [21], this corresponds to production cross sections of ≈ 5 pb for 288 Fl and 289 Fl, respectively. The exact numbers [21,41] agree well with previous measurements summarized in Fig. 11 of [42], and they help to further delineate the shape of the excitation functions for 3n and 4n evaporation channels.
Turning to admittedly unexpected spectroscopy highlights of the even-even isotopes 286,288 Fl [39], Fig. 2a shows an αdecay energy spectrum for the decay step 288 Fl → 284 Cn for the 32 previous world-data events, originating from several discovery and confirmation experiments conducted at different facilities [42][43][44][45][46][47][48][49]. To take into account the various uncertainties of individual events, for a single entry, a Gaussian with integral one and a width compliant with its measured uncertainty was added into the spectrum. One may think that the rather broad doublet structure, which in overall shape is at variance with a common ground-state to ground-state α decay line of an even-even nucleus, relates to somewhat different calibrations of the different experiments. However, the puzzling ≈ 100-keV split [50] was reported independently in all previous experiments, let it be via direct or indirect population of 288 Fl. Figure 2b shows the same spectrum for the twelve chains observed with TASISpec and associated with 288 Fl. The difference in spectroscopic quality is striking, and the above mentioned puzzle is resolved: A single α-decay line at 9.92(1) MeV with a full-width at half maximum of 35 keV emerges from eleven available events. This establishes the decay characteristics of 288 Fl to Q α = 10.06(1) MeV and T 1/2 = T 1/2,α = 0.65 ( 12 8 ) s, in line with earlier values. The event plotted in yellow in Fig. 2b marks the decay chain which comprises a second α decay with E α = 9.33(1) MeV, prior to a fast fission event after only Δt = 518 μs. This chain is highlighted in [39] as it marks the first firm full-energy observation of a weak (≈2%) Q α = 9.46(1)-MeV decay branch from 284 Cn into 280 Ds. Adding the known Q α = 10.79(4) MeV for the 292 Lv → 288 Fl α decay, a smooth Q α sequence across Z = 114 could be experimentally established for the first time. This is illustrated in Fig. 3a and compared with Po-Pb-Hg sequences in Fig. 3b. Note that these sequences cross Z = 82 some 20 neutrons away from doubly-magic 208 Pb. Obviously, there is hardly any kink seen for Z = 114, while it is characteristic for any pronounced shell gap, i.e., magic numbers, such as the well-known one at Z = 82. In essence, that observation rules out all nuclear structure models predicting a distinct magic shell gap at Z = 114.
The second spectroscopic highlight concerns the observation of an α-electron coincidence event in one of the two decay chains starting from 286 Fl [39]. Since the E e = 0.36(1)-MeV conversion electron [21,41] was detected in one of the box DSSD, the detected E = 9.60(1) MeV in the implantation DSSD likely corresponds to the sum of the α decay into an excited state in 282 Cn plus a few tens of keV attributed to atomic electron relaxation processes. Electron binding energies for the K and L shells in Z = 112 copernicium are ≈ 190 keV and ≈ 40 keV, respectively [51]. Based on the known average of E α = 10.19(3) MeV extracted from in total fifteen of 29 hitherto observed decay chains [39], and thus Q α = 10.33(3) MeV for the 286 Fl ground-state to ground-state decay, two reasonable options for a decay scenario emerge for the α-electron coincident event [39,41]. They are displayed in Fig. 4a,b, respectively. Using Q α ≈ 9.70 − 9.75 MeV, b α = 1/29, and T 1/2 = 121 ms, a hindrance factor close to unity can be derived [39,52]. This clearly favours a 0 + spin-parity assignment to the excited state. In turn, this excited 0 + 2 state can either directly decay into the ground-state via a K -converted E0 transition, or via an E2 cascade through the first excited 2 + 1 state expected at ≈ 0.2 MeV excitation energy. In that case, the observed electron stems from L conversion, while the (converted) electromagnetic 2 + 1 → 0 + 1 decay escaped detection. For the first scenario, E α = 9.60 MeV (Q α = 9.74 MeV), E e = 0.37 MeV, and E BE,K = 0.19 MeV sum to E tot = 10.16 MeV (Q tot = 10.30 MeV), which is within uncertainties in line with previously measured or averaged values. One can note that E α = 10.16(1) MeV was measured with TASISpec for the second, full-energy αdecay event of 286 Fl [39]. For the scenario shown in Fig. 4b, MeV, plus some 30-40 keV atomic electron relaxation energy sum to Q tot = 10.33 MeV, which is also within uncertainties in line with measured or average values.
Both decay scenarios were probed against Geant4 simulations, and both decay scenarios were assessed compatible with data, though finally with a slight preference for the scenario shown in Fig. 4a [21,41]. Nevertheless, and irre-  [21,40] spective of the decay scenario, there is an excited 0 + state at E x ≈ 0.6 MeV in 282 Cn, which requires a theoretical explanation. Based on contemporary beyond mean-field calculations exploring the full β-γ deformation plane [9,10,53], shape coexistence is suggested and a pronounced magic shell gap at Z = 114 judged unlikely once again [39]. The global data set on 289 Fl nuclei was roughly doubled, but for the first time the decay chains were interrogated with an α-photon and α-electron coincidence set-up [21,40]. Historically, after some early isotopic re-assignments, five decay chains starting from 289 Fl were reported in experiments conducted at the Dubna Gas-Filled Recoil Separator (DGFRS) [43,45,54]. Four congruent decay chains were observed in TASCA experiments a few years later [42], including one recoil-α-α-α chain reaching 277 Hs [47]. Indeed, this type of chain was confirmed twice in the TASISpec data [21,40]. Four additional 289 Fl decay chains are part of the singleatom chemistry studies conducted behind TASCA [55,56], and nine indirect observations of 289 Fl and its decay were documented following the α decay of 293 Lv [46,48,49,54]. Based on rather broad α-decay energy distributions of the existing measurements, which are displayed in Fig. 5a, b, and predicted single-particle states (cf. Fig. 6) one can anticipate both α-decay fine structure and the potential for X-ray fingerprinting [12] along the decay chain of 289 Fl.
The experimental energy spectra of the TASISpec measurement for the α-decay steps 289 Fl→ 285 Cn and 285 Cn→ 281 Ds are displayed in Fig. 5c, d, respectively. Based on the rather large number of events for a single experiment, superior energy resolution, the observation of one α-photon-electron triple coincidence plus one α-photon and α-electron coincidence each, and following the self-consistency check with Geant4 simulations, the decay sequence sketched in Fig. 6a can be derived [21,40,41]: The main decay branch of 289 Fl with Q α = 9.94(2) MeV, which is indicated by blue arrows in Fig. 6a, proceeds (most likely) into a very low lying state at ≈ 30 keV in 285 Cn. Summing of the α-decay energy and internal conversion electron energies explains the main peak in Fig. 5c, the FWHM of which is somewhat larger than expected for a single plain α decay like the one shown in Fig. 2b for even-even 288 Fl. The triplet peak structure at some 9.1 MeV in Fig. 5d can only be explained by an α decay into an excited state at some 120-160 keV, which cascades to the ground state via a second, low lying excited state. All possible electromagnetic decays need to be dominated by internal conversion. A tentative third state in 281 Ds at about 300 keV excitation energy relates to the 8.92(1)-0.25(1)-MeV α-electron coincidence event, marked with a in Fig. 5d. Depending on whether one puts more weight on the two photon energies [40] or the two conversion electron energies (manuscript VI in [21]) one ends up with numerically slightly different sequences of excitation energies (35-124-266 keV or 60-160-340 keV). However, and including the ground state in the count, the decisive message is the disclosure of a triplet (or quadruplet) of closely lying states of presumably the same parity and mutually connected by highly converted M1 and E2 transitions in 281 Ds. The three decay chains reaching 277 Hs have similar but not the same α energies measured [40,47]. This suggests a Q α = 8.75(2)-MeV decay into an excited state of 277 Hs at ≈ 90 keV, followed by a highly converted M1 transition into the ground state. A second, presumably independent Q α = 9.70(3)-9.17(6)-MeV sequence, indicated by the red arrows in Fig. 6a, can be distinguished from the main branch. Note, however, that implantation-decay time correlations are consistent with both branches starting from the same (ground) state of 289 Fl: Combining all 25 known data points, T 1/2 = 2.3( 6 4 ) s is obtained. Those 18 data points, which can be associated with the main branch, yield T 1/2 = 2.5( 8 5 ) s, while the three data points for the low α-energy branch imply T 1/2 = 2.9( 39 10 ) s [21,40]. There are many examples available from many different types of nuclear structure models and parametrisations for ground-state to ground-state Q α values, i.e., masses of superheavy nuclei. At variance, published sequences of excited (single-particle) states and connecting α-decay probabilities are surprisingly scarce in pertinent literature. Figure 6b-d show three available predictions for the 289 Fl decay chain. Typical for odd-N (odd-Z ) nuclei in the superheavy regime is the competition between high-Ω Nilsson intruder orbitals from the spherical N = 7, j 15/2 (N = 6, i 13/2 ) shell and low-spin natural parity N = 6 (N = 5) states of opposite parity. An interesting side note is that occupancy or nonoccupancy of these high-Ω intruders plays a central role in  [21,40]. The relative energy between the ground and first excited state in 289 Fl is unknown (and can be zero). Dashed levels and α-decay connections are tentative. See text and [21,40] for additional information. b-d are predicted single-particle energies along 289 Fl decay chains and likely α-decay transitions connecting high-Ω states (red) and low-Ω states (blue), respectively. Dashed lines are indicative for alternative branches. The decay schemes use numbers or are adopted from information provided by [57][58][59]. For the calculations shown in a and b, asymptotic Nilsson labels and calculated (ground-state to ground-state) Q α values are available "the interplay of single-particle shell structure, ground-state deformation and dynamic quadrupole collectivity" which can lead to deformation softness or shape coexistence [8]. For the α-decay sequences of odd-A superheavy nuclei, the competition between high-Ω and low-Ω levels implies a high probability for two (or more) α-decay branches, because larger angular momentum difference between initial and final states must be compensated by a significant difference in decay energy to be competitive. In this respect, the superfluid tunneling model [60] is an interesting tool to investigate this kind of competitiveness, with excellent results on even-even [61] and odd-A and odd-odd nuclei [62]. In view of the decay of 289 Fl shown in Fig. 6b Fig. 6a]. In the model prediction shown in Fig. 6c, the [611]1/2 + state is the predicted ground state, while the N = 7 high-Ω levels are presumably beyond the applied 400-keV cut in excitation energy in the relevant Fig. 2 in [58]. These states do appear in the spectrum of the lighter nuclei though. The other difference is the appearance of low-lying negative-parity states, and the α-decay pattern of that prediction is likely to follow a mix of low-spin positive and low-spin negative parity states. In principle, this can give rise to electromagnetic transitions, similar to experiment. Nevertheless, a [611]1/2 + → [761]3/2 − transition requires a significant change in the underlying configurations. Hence, though Q-values are close to experiment and though there is energetically room for electromagnetic transitions, in particular the decay step 289 Fl→ 285 Cn is difficult to bring in accordance with experiment. There is no option for the separate low-Q α sequence either. Finally, the predictions of the symmetry conserving configuration mixing, i.e., beyond mean-field calculations in the full β − γ plane, provide an almost too good match with experiment following the anticipated decay path along positive-parity low-spin states [40,59]. Here, only the relatively large energy shift predicted for the 15/2 − intruder state is seemingly problematic. Interestingly, the very same model is also capable of explaining the smooth Q α sequence across 288 Fl and the low-lying state in 282 Ds (see above and [39]).
Clearly, the spectroscopy results obtained with TASISpec on the decay chain of 289 Fl carry significant direct informa-tion to select between nuclear structure calculations of the heaviest elements.

Decay chains associated with Mc isotopes
The TASISpec 'element 115 experiment' was conducted in 2012. In total, 30 decay chains were associated with the isotopes 287 Mc, 288 Mc, and 289 Mc [63,64]. Comprehensive information on the preparations and the experiment itself can be found in [20,22]. At about the same time, summarizing information on 37 decay chains observed in a series of experiments at Joint Institute for Nuclear Research (JINR) in Dubna and associated with isotopes of Mc was presented in [65], including the first four decay chains seen already in 2003 and published shortly after [66,67]. Undoubtedly, element 115 was discovered at Dubna. In 2015, another 46 α-decay chains associated with 288,289 Mc were reported as a result of a single experiment conducted at Lawrence Berkeley National Laboratory [68]. The latter experiment was, alike TASISpec, sensitive to α-photon coincidences. All three data sets are based on the same fusion-evaporation reaction, namely 48 Ca+ 243 Am, at by and large comparable beam energies at mid target.
The vast majority of the reported decay chains, 99 out of 113, comprise five α-decay steps prior to concluding fission events. They are usually denoted 'long chains'. The remaining 'short chains' have either one (four chains) or two (ten chains) α-decay steps detected prior to fission events. Three out of these 99 long chains are commonly assigned to start from 287 Mc. This is based on various energy-energy, energy-time, and time-time correlations (cf. Fig. 2 in [69]), as well as the standard single-decay step 'Schmidt test' [70] (cf. Fig. 1 and Table 1 in [71]), and the more elaborate congruence check along n-step long decay chains [64,72]. The fact that the two 287 Mc decay chains observed in Dubna experiments were produced at comparatively high beam energies supported their assignment to the 4n-evaporation channel.
The average values of the 287 Mc chain are shown in Fig. 7a. Similarly, basic and advanced statistical methods probing the bulk of the long decay chains [69,71,72], while also applied to beam-energy dependent subsets, did not show any significant sign of not associating these 96 decay chains with the decay of one single radioactivity [64], i.e., starting from one (ground) state of odd-odd 288 Mc. This isotope assignment was more recently confirmed by the first direct mass-number measurement of superheavy nuclei [73]. The average values of these 96 288 Mc chains are shown in Fig. 7c. Nevertheless, a note of caution might be appropriate here: Despite the majority of the data having been taken at beam energies near the anticipated peak cross section of the 3n evaporation channel, 288   From a nuclear structure point of view, the observation [63] and confirmation [68] of a 237-keV E1 γ ray in 276 Mt provides an intriguing benchmark for nuclear structure theory. It relates to the relative positioning of high-Ω intruder orbitals vs. low-Ω natural orbitals already discussed in connection with the observed fine structure along 289  Accordingly, at least one of these combinations of Nilsson orbitals must be close to the Fermi surface in 276 Mt. This has been nicely illustrated with calculations using contemporary Skyrme energy density functionals [74]. The parameter set UNEDF1 SO , which was optimized for heavy Z ≈ 100, A ≈ 250 nuclei, fails to place any of the two combinations at the Fermi surface for 276 Mt. In turn, the standard UNEDF1 parameter set does suggest the E1 in the proton sector. This is discussed in detail in [74] and highlighted in connection with Fig. 2 in [38].
K X-ray candidates at proper photon energies are present in spectra coincident with the 280 Rg [68,71] and 276 Mt α decays [63,68,71]. However, unequivocal fingerprinting of the proton number like, for instance, for rutherfordium [12] is not possible with the at present still limited number of αphoton coincidence events. The reasons are either the seemingly complex γ decay scheme of 272 Bh [63,68,76], or the relatively small conversion coefficients of the E1 transitions in 276 Mt [63,68]. In either case, they can lead to background at the K X-ray energies of interest due to Compton scattered events [71].
For consistency, the fourteen short recoil-α(-α)-fission chains are labelled according to [64,72,77], i.e., D1-D4 for four such chains reported from Dubna [65], T1-T7 for the seven chains detected with TASISpec [64], and B1-B3 for the three chains observed in the Berkeley experiment [68]. Based on standard statistical assessments of correlation times per single decay step [69,70,77] and more elaborate tests of correlation times along decay chains [20,64,72,77], it has been shown that the hypothesis of this set of fourteen decay chains forming a single congruent ensemble can be falsified with an error probability of < 1%. The straight-forward solution is to single out decay chain D3 and associate it with the 2n evaporation channel 289 Mc, as it is done in Fig. 7f. This action is further supported by the reported α-decay energies of D3; for the first decay step, it is lower than the vast majority of decay chains associated with 287,288,289 Mc, and it lies clearly outside the common range of energies for the second decay step (see, for instance, Fig. 1 in [64], Fig. 1 in [77], or Fig. 2 in [69]).
The correlation times of the remaining thirteen short chains were found to be compatible with representing one ensemble [64,72,77], and two extreme but in terms of correlation-time analyses sound cases have been studied and presented in depth in [64]: All thirteen chains either form a separate, second 289 Mc decay chain, or all thirteen chains could be readily combined with the long decay chain of 288 Mc. This leads to either small fission or β + /EC branches of 284 Nh and 280 Rg, in the latter case followed by fast fission of the even-even daughters 284  Third, there could be even more decay paths and connections [62].
The nuclear structure aspects of having two similar but different decay paths of 289 Mc have been discussed previously (see Fig. 3 in [64]). In accordance with the discussions of finestructure pathways of 289 Fl (cf. Fig. 6) and 293 Ts (cf. Fig. 11), high-Ω N = 6 [606]13/2 + or [615]11/2 + positive-parity Nilsson orbitals are predicted to form one decay path, and low-Ω N = 5 negative-parity orbitals are suggested to form the other one. Other models predict [78], or at least indicate [79], the same pattern of at least two possible decay paths of 289 Mc, while details like Q α values, half lives, hindrance factors, or overall complexity are model dependent and eventually require detailed assessments of α-decay probabilities [62].
As expected, predictions are comparable for 287 Mc, but with less evidence for two parallel decay paths for some of them. Figure 8 shows the few available examples, representing in panel Fig. 8a a macroscopic-microscopic approach based on the Woods-Saxon potential [80], in panel Fig. 8b a modified two-center shell-model approach adopted from Fig. 1 of [78], and in panel Fig. 8c the Skyrme energy-density functional UNEDF1 using numbers from Tables IV-VI of The macroscopic-microscopic model predicts at least two different α decays starting from either the [606]13/2 + level at E x = 0.29 MeV or the [541]1/2 − ground state of 287 Mc. The latter populates low-spin negative-parity states in the daughter 283 Nh, and its [512]3/2 − ground state is predicted to continue to α decay into ground or excited states of the same or similar Nilsson orbitals. On average, the predicted Q α values [81] are close to experiment, including the drop in the last decay step. Note that Nilsson singleparticle states in 267 Db are lacking in [80]. There is a risk Table 1 Overview of correlation-time analyses of single decay steps according to [70] of the first four decay steps of the long tennessine chains previously assigned to start from 294 Ts (L294). The ensemble S293A is taken from the second column of Table 2 for completeness. The combination with the moscovium chain D3, which is also singled out in Fig. 7f, is presented in the rightmost column. The tag 'H' indi-cates a failed 'Schmidt test'. The labels of the ensembles are defined in the text and Fig. 9 [74], and Q α values are not given there either.

How to connect moscovium and tennessine?
Hitherto, in total 24 decay chains originating from isotopes of element 117, tennessine, have been reported [82][83][84][85][86][87]. The first decay chains were observed at Dubna via the reaction 48 Ca+ 249 Bk [82,83]. The enriched 249 Bk target material was provided by Oak Ridge National Laboratory. Given the rather short half life of T 1/2 = 330(4) y, the production, extraction,  Fig. 7f. Panel e shows average values for the remaining eleven short chains, S293B, which is combined with the four moscovium chains shown in Fig. 7e. See text for more discussion and isotopic enrichment of the required amount of tens of milligrams is a major effort [88,89], awarded with the name of the element a few years later [90,91].
In connection with a second round of provision of 249 Bk material from Oak Ridge, experiments at Dubna and GSI confirmed the early findings [84][85][86][87]. Similar to the decay chains starting from moscovium isotopes (cf. Fig. 7), 'long' and 'short' decay chains from tennessine have been observed. The six long chains, four from Dubna experiments [85], two from the GSI experiment [86], comprise up to seven α decays. Their correlation times are statistically sound [72] and the αdecay sequence is commonly assigned to the odd-odd isotope 294 Ts. This ensemble of six chains is labelled 'L294' in the following. It is summarized in Table 1, and it is illustrated in Fig. 9c.
The situation is more complicated for the by now in total eighteen reported short decay chains. Sixteen were observed in Dubna [85] but only two at GSI [87]. All of them have been assigned to start from the isotope 293 Ts. They comprise either three α-decay steps, i.e., conclude with a spontaneous fission of 281 Rg, or four α-decay steps, i.e., conclude with a spontaneous fission of 277 Mt. The former refer to fourteen Dubna chains, while the latter relate to two chains from Dubna and both chains from GSI.
Decay chains populated via 249 Bk( 48 Ca,4n) 293 Ts feed those populated via 243 Am( 48 Ca,2n) 289 Mc. This kind of cross bombardment provides one out of several means [92] for isotope assignments and element discovery stories [93,94]. In fact, the official approval case of discovery of element 115 and element 117, published early 2016, relies primarily if not exclusively on the above mentioned cross bombardment. Unfortunately, the case is built upon the four in themselves non-congruent [72,77] short moscovium decay chains reported from Dubna [65], in combination with a subset of ten out of the sixteen short tennessine decay chains reported from Dubna [85]. Notably, neither these ten nor all sixteen or eighteen short tennessine decay chains are compatible with a single radioactivity either [72,77]. For instance, already the first decay step does not pass the 'Schmidt test', as shown in Table for the corresponding ensemble labelled 'S293'. In essence, the two elements were approved based on a hypothesis 1 which can be falsified on the 2-σ level [72,77]. Moreover, the approval case also ignored the short moscovium chains detected at GSI and Berkeley, which were (made) available prior to approval [68,95], and which clearly show that the Dubna chain D3 is an outlier [77] [cf. Sect. 3 and Fig. 7f]. In fact, the task is to investigate ensemble S293 in more detail to identify a spot for the peculiar moscovium decay chain D3, and as a result suggest an alternative but statistically sound connection of decay chains of the two elements based on current world data.
There are many correlations between decay energies and correlation times between the first few decay steps of tennessine chains that can be studied. A selection of four correlations is displayed in Fig. 10. Based on such correlations and subsequently probing any single of the short tennessine chains against a reference chain formed by the congruent long tennessine chains, in the same fashion as it had been done 1 It is important to state that the discussion in this contribution does not intend to challenge the discovery, naming rights, or the element names of moscovium and tennessine. There is clear evidence and there are many good arguments that isotopes of these two elements were discovered in joint Russian-US experiments conducted at JINR Dubna. In fact, a statistically sound cross-bombardment link between the two elements is suggested here. Further, congruent long decay chains reported first from Dubna have been confirmed in several independent experiments at GSI and Berkeley. The point is that many aspects of the report forming the basis of official approval are scientifically questionable. That report is a good case for a university course in science ethics -go ahead and identify mathematical and methodological flaws and debatable use of previously published material.  Table 2). Correlation times reported for recoil-α events of isotopes of tennessine are plotted on the y axis, while the x axes refer to their decay energies (a) as well as the correlation times between Ts-Mc events (b), Ts-Nh events (c), and Ts-Rg events (d). In case of missing α events along the chain, subsequent correlation times were corrected for the half life of the missing decay step. This has not been done for the numerical analysis in the two tables for moscovium chains [64], it is suggestive to split ensemble S293 into S293A (seven chains) and S293B (eleven chains). Details are provided in Table 2, while S293A is also listed in Table 1 for completeness. In the notation of [77], ensemble S293A consists of S01, S04, S05, S08, S09, S10, and S15. Figure 10a shows that on average somewhat higher decay energies correlate with somewhat smaller correlation times, and that the entries of ensemble S293A (blue) are similar to those of ensemble L294 (black). The difference between S293A and S293B (green) becomes more pronounced when looking at the time-time correlations in panels b-d: Data points in green tend to be located in the lower left triangle, all the others towards the upper right section in each of the three plots. In the advanced statistical test, chains S01, S04, S05, S09, and S10 show values similar to the long chains themselves, while the overlap is somewhat weaker for S08. S15 (missing first α decay) is a rather arbitrary choice based on its reported decay energies. If one compares column S293A in Table 1 with the reference column L294, they obviously fit well together both in terms of derived half lives and decayenergy ranges, i.e., they can very well be considered forming the same decay chain with an ≈ 50% fission branch at roentgenium. This option is illustrated in Fig. 9b-d. Moreover, this mixed ensemble M29X provides an excellent explanation for decay chain D3 associated with 289 Mc. This can also be seen in the respective columns of Table 1. Here, all five D3 observables, i.e., three correlation times and two decay energies, are found compatible with ensemble M29X. The final suggestion for these 6 + 7 + 1 = 14 chains is displayed in Fig. 9a. For chain D3, either representing the 1n-channel of the reaction 48 Ca+ 243 Am or being due to a reaction on 244 Am target nuclei is unrealistic. Therefore, the combined ensemble L29X forms a decay sequence of 293 Ts rather than 294 Ts. A possible counter argument arises from the reported mid-target beam energies and thus initial 297 Ts compound nucleus excitation energies. However, and folding in GSI data on element 117, the combined cross-section plot shown as Fig. 8 in [87] turns out to be ambiguous. One may recall that setting up the planned beam energy and then maintaining it over long periods and presumably different beam tunes can be a challenging task at cyclotron facilities.
The remaining eleven short tennessine decay chains, S293B in Table 2 and Fig. 9e, form a consistent ensemble in themselves, i.e., there is no need for an additional split. In Table 2 and Fig. 9f they are combined with four chains associated with 289 Mc, the latter based on E α (Mc)< 10.4 MeV, as mentioned earlier. These four chains end in 281 Rg, like seven out of the eleven chains of the S293B ensemble. It may therefore be tempting to separate those four S293B chains ending in 277 Mt. However, statistically it is not a major problem that four out of eleven tennessine decay chains yield an α-decay branch of 281 Rg, but none of the four moscovium chains. With the alternative association of chain (a) (c) (b) Fig. 11 Predicted single-particle energies for 293 Ts and 289 Mc and likely α-decay transitions connecting high-Ω positive-parity states (red) and negative-parity states (blue), respectively. The decay schemes use numbers or are adopted from information provided by [74,78,80]. Asymptotic Nilsson labels are used, and calculated (ground-state to ground-state) Q α values are included provided availability [78,98] Tables 1 and 2 , respectively, both starting with 293 Ts, are one possible solution of the task outlined earlier, but neither an exclusive one, nor is the final word claimed to be spoken.
Under the assumption that there are two parallel decay chains starting from 293 Ts, they can be confronted with the few available calculations of single-particle Nilsson states of 293 Ts and its α-decay daughter 289 Mc. The predictions of three models are illustrated in Fig. 11. The microscopicmacroscopic approach of Parkhomenko and Sobiczewski [80] Fig. 8b]. This implies that excited high-Ω positive-parity states can decay by γ -ray emission into the [505]9/2 − ground state of 293 Ts, which readily α decays into its analog state in 289 Mc. However, different to the 287 Mc→ 283 Nh decay discussed in connection with Fig. 8b, once electromagnetic decays within 293 Ts reach its isomeric [510]1/2 − level, this gives rise to second α-decay branch. The predicted Q α values are close to experiment. Calculations using the Skyrme energy density functional UNEDF1 suggest a practically degenerate 'ground-state pair' of [510]1/2 − and [606]13/2 + states, which will both α decay. The low-Ω negative-parity state has its low-lying analog Nilsson state to decay into, while the [606]13/2 + level encounters some angular momentum and energy hindrance, such that an additional branching may occur. In conclusion, all three models rather support than reject the suggested presence of two parallel α-decay sequences starting from 293 Ts.

From TASISpec to LUNDIUM
There are several components to further advance the nuclear spectroscopy opportunities for the heaviest elements. With the Lund group being engaged in the superheavy element research programme at GSI/FAIR, the new superconducting continuous-wave heavy-ion linear accelerator HELIAC [99] is going to deliver an increase in ion-beam rate and thus superheavy element production and observation rate. Provided that actinide targets can withstand the ion-beam related heating, also an increase in their thickness beyond 1 mg/cm 2 might become an option. However, given that these isotopically enriched materials are, besides being radioactive, extremely rare and precious, this is a major technological challenge (see, for example, [100] or Chap. 4 in [101]).
With these efforts and with the real cost of beam-time hours in mind, it is mandatory to also consider upgrades of the decay station. This can involve improvement of passive shielding and inclusion of active shielding to prevent unwanted background events in the silicon and germanium detectors. A higher granularity in the implantation DSSD pixelation allows for higher rates or prolonged correlation periods, respectively. Last but not least, a more compact germanium detector arrangement could increase α-photon coincidence efficiency. The method must rely on the extraction of as much experimental spectroscopic information as possible for every single one of these extremely rarely produced superheavy nuclei.
These ideas are realized in the Lundium decay station, a project financed by the Knut and Alice Wallenberg foundation. The major investment and technological development concerns a novel type of composited germanium detectors named Compex [102]. As sketched in Fig. 12a, TASISpec uses existing cluster (7 Ge crystals) and clover (4) Fig. 12 Evolution of set-ups towards the new Lundium decay station. All drawings are generated from Geant4 simulations performed by Sarmiento [31]. a Positioning of germanium-detector crystals (coloured shapes) and the pixelated silicon-detector box (semi-transparent white) of the existing TASISpec arrangement [16]. b Set-up using four custommade four-crystal Compex detectors [102] behind the box DSSD and one former Euroball cluster detector backing up the implantation DSSD. c Envisaged germanium detector configuration of the Lundium chamber (cf. Fig. 14). The seven-crystal cluster is replaced by another Compex module, but with an electrical fourfold segmentation of each crystal. See text for details nium detectors, which were available at GSI or made available through European networks. During experiments, the pixelated silicon cube detecting implanted superheavy elements and their α decay and fission is enclosed in a thin aluminium housing. The germanium detectors are thus outside the dilute helium atmosphere of the gas-filled separator TASCA. This also means that low-energy photon detection is impaired by unwanted material causing significant loss of photon detection efficiency -in particular for superheavyelement L X-rays at 30-35 keV (cf. Fig. 13). Despite the close proximity to the silicon cube, there is still a distance of a few centimetres between silicon detectors and germanium crystals. Nevertheless, TASISpec reaches ≈ 40% efficiency for K X-rays [16,31].
Several Geant4 simulations investigating different germanium set-ups were conducted and led to configurations drawn in Fig. 12b, c. Common for them is the replacement of the four old clover detectors with novel, custom-made Compex detectors. These comprise four 5×5×5 cm 3 cubes of encapsulated high-purity germanium crystals next to each other, and they use a modern type of electromechanical cryocooler [102]. Additionally, they can be placed inside the vacuum of a new Lundium detection chamber, as illustrated in Fig. 14. No window at all besides the encapsulation is required, which maximizes low-energy photon detection. The difference in panels (b) and (c) is the type of germanium detector placed behind the implantation silicon detector: Panel (b), Compex-4, shows an existing cluster. In fact, with all germanium detectors still surrounding the TASISpec DSSD cube, this was the set-up used for the 'element 114 experiment' [21]. Positioning a fifth new Compex-type detector behind the implantation detector yields yet another decisive improve- (e) (d) Fig. 13 a Gain in photon detection efficiency for (low-energy) photons comparing the set-ups illustrated in Fig. 12b (Compex-4 plus Cluster, blue dots) and Fig. 12c (Compex-5S, cyan stars, Lundium configuration) with the TASISpec version shown in Fig. 12a (dashed black line). Gamma-ray decays of excited superheavy elements are expected at energies of about 50-400 keV. Note the logarithmic scale on the x-axis with a focus on low-energy photon detection efficiencies. The X-ray gain factor, i.e., the figure-of-merit, is defined as the product of the gain in detection efficiency for a coincidence between L and K α X rays from a superheavy element. Panels b-e provide Geant4 simulations of photon spectra confronting the existing TASISpec set-up (black) with the Compex-5S system (cyan). The physics input is the simplified test case shown in the inset, namely an electromagnetic E1 (b) or M1 (c) 237-keV decay in meitnerium, Z = 109. Panels d and e show L X rays in coincidence with K α2 and K α1 transitions, respectively, i.e., the decisive measure for a potential Z -fingerprint of a new superheavy element ment in photon detection efficiency. In addition, electrical four-fold segmentation of each of the four crystals of this special Compex detector provides a significant improvement of spectral quality -this Compex-5S configuration became therefore the choice for the Lundium set-up.
As mentioned earlier, the crux of the new concept, however, is to place the Compex detectors in vacuum directly behind the silicon detectors. The idea is to remove the end-cap of the cryostat housing, which is possible due to anticipated operation in vacuum. This is not possible with any of the existing detectors, and this is the essential improvement that will lead to the highly desirable increase in relative detection efficiency for low photon energies. In Fig. 13a the relative gain in efficiency as a function of photon energy is shown with respect to the existing TASISpec (dashed black horizontal line). The data points in Fig. 13a are derived from comprehensive Geant4 simulations [31] of 133 Ba and 152 Eu source ions implanted into the silicon detector with a depth and position profile representative of superheavy element experiment conditions. Overall, an increase in relative efficiency of K ,rel = 1.2-1.3 at superheavy element K X-ray energies and L ,rel = 1.7-1.8 at L X-ray energies is achieved for Compex-4, as shown by the blue circles in Fig. 13a. The product of those two numbers, F = K ,rel · L ,rel ∼2.1 for Compex-4 with respect to F 0 =1.0 for TASISpec, may serve as the spectroscopic figure-of-merit of a given detector configuration. It represents to first order the increase in performance for detecting a coincidence between subsequent atomic K α and a L x ray following an electromagnetic decay of an excited nuclear state by inner conversion. With fluorescence yields [103,104] of more than 90% and 50%, respectively, the coincident observation of two X-ray fingerprinting energies is considered to provide the best possible evidence for the element number, Z , on (very) few-event statistics, based on the spectral analysis on the right hand side of Fig. 13.
Closing in with a fifth Compex detector right behind the implantation silicon detector implies that each of its four germanium crystals covers a rather large solid angle. This gives rise to a problem: possible energy summing due to the simultaneous detection of more than one photon in the same crystal. For instance, instead of measuring a 35-keV L X-ray in one crystal and a coincident 140-keV K α X-ray in a second crystal, a sum of 175 keV would be registered in one of the closed-in crystals. In more general terms, summing leads to unwanted, spurious peaks at summed energies of intense photon transitions of the respective source and has to be kept to a reasonable minimum. Therefore, a fourfold electrical segmentation of each of the four crystals of Compex module 5 is foreseen. The figure-of-merit increases to F = 3.8 for Compex-5S.
On the right hand side of Fig. 13, Geant4 simulations of photon spectra show the order-of-magnitude increase of (Xray) coincidence spectroscopy with Compex-5S (cyan) compared to TASISpec (black). The physics input to the simulations is the simplified test case shown in the inset, namely an electromagnetic E1 (b) or M1 (c) 237-keV decay in meitner- ium. In case of M1 decay, the transitions are highly converted, i.e., instead of the γ -ray line at 237 [panel (b)], K -x rays dominate the spectrum in panel (c). In line with increased relative efficiencies illustrated in Fig. 13(a) statistics are roughly doubled for Compex-5S vs. TASISpec. Most importantly, however, Compex-5S allows for X-ray coincidence spectroscopy: Panels (d) and (e) reveal a significant number of L X rays in coincidence with K α2 and K α1 transitions, respectively. Such simultaneous X-ray energy measurements are unique for an element and thus the most decisive fingerprint for singleor few-event element determination. Chances are almost an order of magnitude larger for Compex-5S compared with TASISpec and reach several per cent for heavier and heavier superheavy elements. Figure 14 provides a view into the new Lundium chamber. The downstream segmented Compex module is hidden inside an anti-Compton shield assembled from 28 BGO crystals sized 19.1 × 19.1 × 100.0 mm 3 [105]. Each of the four Compex detectors behind a box DSSD is protected against radiation from the beam direction with an anti-Compton assembly of twelve BGO crystals, so-called 'half-pockets' [105]. They are not shown in Fig. 14 for clarity. Compared with TASISpec, the strip pitch of the implantation DSSD is reduced from 1.8 to 1.0 mm, i.e., the number of pixelation is increased by more than a factor of 3. These detectors and implantation-veto correlations were subject to recent BSc thesis work [106,107]. There was no spectroscopic reason identified to increase the pixelation of the box DSSDs.
For the real Lundium chamber, we expect the delivery of the segmented Compex module in 2022, likewise a few vacuum pieces and BGO read-out electronics. It will then be subject to intense testing and characterization with standard α, electron, and γ -ray sources.

Summary
To summarize, a brief overview of experimental and data analysis tools used in connection with TASISpec experiments has been presented. For more details, see in particular the PhD theses of Forsberg and Såmark-Roth [20,21]. Concerning envisaged improvements of detection, the Lundium decay station employs novel concepts and technologies for compact high-resolution germanium detectors with unprecedentedly high detection efficiencies and precision placed directly behind the silicon detectors registering superheavy ion implantation and their charged-particle decays. The hope is to provide a spectroscopy set-up which ensures improved direct access to the nuclear structure of the heaviest manmade atomic nuclei, as well as X-ray fingerprinting their atomic number.
Concerning element discovery, alternative but statistically sound possible links between decay chains of 293 Ts and 289 Mc have been suggested. Clearly, more systematic data of spectroscopic quality is required. Coincidences between α particles and electrons or photons are likely to add important additional information to finally unravel the tangles of these chains.
Spectroscopic results on superheavy nuclei have been obtained which are capable of benchmarking nuclear structure models. These include several series of precise Q α values across Z = 114, flerovium. The experimental findings imply that the next pronounced magic number beyond Z = 82, lead, is not at Z = 114, which is in line with many contemporary mean-field and beyond mean-field model calculations. Triaxial shapes, shape coexistence, and shape driving intruder orbitals are seemingly at play at Z = 114. High-Ω intruder orbitals are also involved in issuing a 237-keV E1 transition in 276 Mt, the presence of which near the Fermi surface of 276 Mt should be also be demonstrated by any nuclear structure model of relevance. The measured αdecay fine-structure pattern in even-Z -odd-N nuclei along the 289 Fl decay chain is likely to engage a multitude of N = 6 low-Ω positive-parity orbitals rather than N = 7 high-Ω intruder orbitals. These are suggested to form an independent decay path. Starting from the somewhat provisional discussion of Figs. 6, 8, and 11, time has come to assess model and model parametrisations in greater detail along predicted single-particle states, collective states on top of these if applicable, and more elaborate α-decay probabilities or hindrance factors. As a result, more robust predictions should emerge for experimentally unexplored territory.
The more recent spectroscopic achievements along flerovium decay chains were commented on in the referee report of [39]: "The wealth of physics discussion that a single event of this kind can generate is both stimulating and satisfying ... and demonstrates the robust promise for superheavy spectroscopy in the near future". Technical developments on beam intensity, target integrity, and detection efficiency should allow to "wring out tantalizing physics from compound nucleus production data where cross-sections are in the picobarn range".
Acknowledgements Discussions with and enduring superheavy engagement from Lund colleagues present and past are gratefully acknowledged; U. Forsberg and A. Såmark-Roth, D.M. Cox and L.G. Sarmiento, C. Fahlander and P. Golubev, to name but a few, likewise colleagues from the Division of Mathematical Physics. The experiments would not have been possible without decisive support from the international TASCA collaboration. In Lund, we are particularly indebted to invaluable efforts of the technical staff at GSI Darmstadt and the radiochemistry laboratory at JGU Mainz, as well as the workshops at the University of Cologne. Viable chats with M. Krug during our sabbaticals should not go unmentioned.

Note added in proofs:
Recently, the observation of additional 55 long and six short decay chains associated with isotopes of moscovium was reported [108]. The decay characteristics of the long chains are considered consistent with those previously assigned to start from 288 Mc (cf. Fig. 7b and c and [69,71]), but an ∼50% a-decay branch of 268 Db is newly proposed [108]. The decay characteristics of the six new short chains are detailed in their Supplemental Material [108]. In view of earlier assessments [64,69,71,72] (cf. Sec. 4), the new decay chains 1, 2 and 5, 6 are compliant with the nine decay chains shown in Fig. 7d, i.e., they can start from 288 Mc and exit in 280 Rg and 284 Nh, respectively. Decay chains 3, 4 are candidates for the ensembles displayed in Fig. 7e or Fig. 7d. Overall, the scenarios discussed and suggested in Secs. 4 and 5 remain unchanged.
Funding Information Open access funding provided by Lund University. TASISpec and the Lundium upgrade were made possible by a series of contributions from the Royal Physiographic Society of Lund, the Crafoord foundation, and major support from the Knut and Alice Wallenberg foundation (KAW 2015.0021).

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: This review article relies on data already published elsewhere. See appropriate references throughout the article.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.