Controlling 229Th isomeric state population in a VUV transparent crystal

The radioisotope thorium-229 (229Th) is renowned for its extraordinarily low-energy, long-lived nuclear first-excited state. This isomeric state can be excited by vacuum ultraviolet (VUV) lasers and 229Th has been proposed as a reference transition for ultra-precise nuclear clocks. To assess the feasibility and performance of the nuclear clock concept, time-controlled excitation and depopulation of the 229Th isomer are imperative. Here we report the population of the 229Th isomeric state through resonant X-ray pumping and detection of the radiative decay in a VUV transparent 229Th-doped CaF2 crystal. The decay half-life is measured to 447(25) s, with a transition wavelength of 148.18(42) nm and a radiative decay fraction consistent with unity. Furthermore, we report a new “X-ray quenching” effect which allows to de-populate the isomer on demand and effectively reduce the half-life. Such controlled quenching can be used to significantly speed up the interrogation cycle in future nuclear clock schemes.


Introduction
The exceptionally low-energy nuclear excited state of the isotope 229 Th has attracted considerable attention because of its potential use in ultra-precise optical clocks 1,2,11,12 .With only a few electron volts of excitation energy, this laser-accessible first excited state is predicted to have a long half-life exceeding 10 3 seconds (called an "isomer"), making it a candidate for a clock state.
Given that direct excitation experiments from the ground to the isomeric state 229m Th are ongoing worldwide 13,14 , a more precise determination of the isomer half-life, and a measure of radiative decay fraction in material environments is highly desirable.The spectroscopic studies on the isomer transition and the developments of the nuclear clock are currently focusing on two different physical implementations: Th 3+ in an ion trap 15 and thorium-doped ionic single crystals with a large band gap 6,7 .A valuable feature of the latter "solid-state nuclear clock" using thorium-doped crystals is that a much larger number of thorium atoms (≈ 10 18 cm −3 ) can be used than in an ion-trap-based clock.This allows the solid-state nuclear clock to be used as the fast-averaging reference clock, wherein the temporal stability converges with a short averaging time.
The development of a solid-state nuclear clock requires a superior-grade 229 Th-doped crystal, characterized by an adequate concentration, facilitating isomer excitation and detection of the isomeric decay, while suppressing non-radiative decay processes in the crystal.Theoretical and experimental studies on thorium-doped crystals such as Th:CaF 2 suggest that the thorium substitutes a Ca 2+ ion in charge state Th 4+ in the crystal 16,17 , realizing [Rn] electronic configuration.The absence of energetically near-by electronic states should suppress the internal conversion of 229m Th and make direct radiative VUV photon emission the dominating decay process.
One of the crucial matters to verify is whether the excited isomeric state indeed shows the predicted long radiative lifetimes also in the crystal environment.In 2023, a successful optical measurement of the radiative decay of 229m Th was reported for the first time 8 .The authors implanted the radioisotopes 229 Fr and 229 Ra into MgF 2 and CaF 2 crystals at the CERN-ISOLDE facility, where the isotopes undergo successive beta-decays to 229 Th with an uncertain branching into 229m Th.The wavelength of the detected photon from the isomer state was determined to be 148.71(41) nm by a VUV-spectrometer corresponding to an isomer energy of 8.338 (24) eV.The half-life of the isomer in the MgF 2 crystal was determined to be 670(102) s through the analysis of a time spectrum encompassing multiple beta decays and isomeric decay.If we restrict our discussion to cases where the charge state of thorium is such that it is 4+, i.e.where internal conversion decay of 229m Th is forbidden, the half-life of 229m Th is expected to be reduced by a factor n 3 (n being the refractive index of the crystal at the relevant wavelength) as compared to the vacuum half-life 18,19 .Influences of the shallow implan-tation depth (17 nm) compared with the wavelength of the VUV photon from 229m Th, realized in ref. 8 , require further investigation.
In this work, we synthesized the 229 Th-doped CaF 2 crystal through single crystal growth with a concentration of the order of 10 18 cm −3 (ref. 17).Subsequently, we excited the isomer state from the ground state of the doped 229 Th nucleus with a resonant X-ray beam, and observed the radiative decay to the ground state accompanied by the emission of a VUV photon.Using 229 Th doped crystals realizes a controlled chemical environment 17 in which the excitation to 229m Th does not introduce disturbances due to radioactive decay or implantation of the parent isotopes of 229 Th which could promote non-radiative decay paths.
The closed-loop character of our experiment, with time-controlled isomer excitation followed by VUV-signaldetection, is a precursor for solid-state nuclear clocks using laser excitation, and provides useful additional information for the excitation and decay dynamics of 229m Th.One such important finding is the accelerated decay of the isomer state during X-ray beam irradiation, presenting a potential key process for quenching the remaining isomer population to the ground state: "re-initializing" the interrogation cycle for solidstate nuclear clock operation.
The paper is structured as follows: Following this section, the experimental method is described in Sec.2.The analysis and results are presented in Sec.3.Finally, the discussions and conclusions are given in Sec.4.Additionally, six appendices provide detailed descriptions of the Th-doped crystal, the data acquisition scheme, the systematic errors in our analysis of the isomer half-life, detection efficiency and the signal wavelength.

Experiment
The experiment was performed at the BL19LXU beamline of SPring-8 20,21 .The data presented in this paper was taken in May and July 2023.The energy of the X-ray beam is tuned to the second excited state of the 229 Th nucleus at 29189.9 eV (see Fig. 1(a)).The X-ray beam is monochromatized to an FWHM bandwidth of 30 meV by a Si monochromator system consisting of Si(111), Si(660) and Si(880).The beam has a spot size of approximately 1.0 × 0.5 mm 2 at the target position with an intensity of 2 × 10 11 photons per second at maximum.The beam energy and intensity were regularly monitored by an absolute X-ray energy monitor 22 and free-air ionization chambers during the experiment.Before searching for the VUV-photon from 229m Th decay using the VUV-detection system, the absolute center energy of the X-ray beam was confirmed at the NRS-detection system (see Fig. 1(b)) by detecting the nuclear resonant scattering (NRS) signal from the second excited state of 229 Th nuclei with a separately fabricated thorium nitrate target 23 .This resonance curve of the X-ray signal is obtained (Fig. 2 measurement can be found in refs. 20,24,25 T229 Th-doped CaF 2 crystal targets were developed in the laboratory of TU-Wien 17 (see Appendices).Two cuboidal pieces of 229 Th:CaF 2 of appropriate size ∼ 1 mm 3 (named "Basso" and "Alto") were cut from the same crystal ingot and facets were polished to optical quality.The measured 229 Th activities were 20.0 kBq (17.2 kBq), yielding a 229 Th doping concentration of 4.0 × 10 18 (4.7 × 10 18 ) cm −3 for the "Basso" ("Alto") crystal, respectively.At such high concentrations, radiolysis during growth leads to a fluoride-deficiency of the obtained crystal, reducing VUV transparency to <1 % in the entire VUV region.Annealing the crystal under a fluoriderich atmosphere recovered a transmission of >40 % at 150 nm wavelength 26 .
The experimental setup for VUV-photon detection consists of two connected vacuum chambers ("irradiation" and "measurement") with a vacuum level of 10 −3 Pa (see Fig. 1 (b)).The target crystal is mounted on a moving holder (Fig. 1 (e)), which translates between the two chambers using a pneumatic stage.The crystal is first irradiated by the X-ray beam for a certain time to populate the 229m Th state in the irradiation chamber and is then moved to the measurement chamber within about one second.This transfer scheme is used to sup-press X-ray beam-induced optical backgrounds.
The 229 Th:CaF 2 crystal is positioned at the focal point of a parabolic mirror with a protected aluminum coating so that the emitted light is reflected and passes through a spectral filtering device as shown in Fig. 1 (c).The spectral filtering device consists of right-angle prisms with special dielectric coatings, so that only light in a specific wavelength region around 150 nm is reflected.The use of four consecutive rightangle prisms (Fig. 1(d)) with installation directions orthogonal to each other guides the signal beam to the detector while reducing the number of background events down to O(10 −6 ) in the calculation.A set of selectable band pass filters (BPFs) installed on a rotating wheel is introduced for spectroscopic measurements of the radiation emitted by the crystal.The filtered light is focused by a MgF 2 lens, and finally detected by a solar-blind photomultiplier tube (PMT, Hamamatsu R10454).The dark count rate of the PMT (called VUV-PMT) was reduced to less than 0.1 Hz by cooling it to −30 • C.
Another PMT (Hamamatsu R11265-203, called VETO-PMT), which is sensitive to ultraviolet (UV) and visible (VIS) light, is installed behind the first right-angle prism to implement temporal filtering of the optical signal (see Fig. occasional alpha-and beta-decay of 229 Th and its daughter elements is known to produce bursts of order O(10 4 ) photons per decay event, mainly in the UV region 27 , which produces a background also on the VUV-PMT.Events detected by the VUV-PMT that coincide with such bursts detected by the VETO-PMT are rejected using a timing coincidence of the two PMTs.
The waveform signals from the PMTs are amplified and then recorded by an oscilloscope (National Instruments, PXIe-5162).See Appendices for details of this signal processing.During the measurement, reference pulse signals are generated by a delay generator.They are combined with the VUV-PMT signal and used as the trigger channel.The reference pulse data are used to monitor the dead time of the data acquisition, and the remaining signal after background rejection.

Analysis and Results
The clear evidence for VUV photon emission from the 229 Th isomer is the enhancement of the photon counts in the VUV-PMT when tuning the X-ray energy to be resonant with the 29-keV level.Figure 2 (b) shows a resonance curve of the measured VUV-photon rate, using the Th:CaF 2 crystal and the VUV detection system.The crystal was irradiated for approximately 600 seconds for each X-ray energy point and VUV photons were counted for 1,800 seconds after each X-ray irradiation.The resonance energy is obtained by fitting a Gaus-sian function plus a constant offset value.Taking into account the systematic uncertainties of the X-ray energy calibration 22 , the peak energy is determined to be E = 29189.89± 0.07 eV.This energy coincides with that of the X-ray signal spectrum within the error, validating that the observed VUV signal is the result of isomeric decay.
The small difference in the resonance width (σ = 18.2(25) and 13.4(11) meV in Fig. 1(f) and (g), respectively) between these two resonance curves could be caused by the effects of the solid-state environment in the different targets: the resonance width is determined not only by the monochromator bandwidth but also by the target-dependent inelastic scattering processes corresponding to the phonon contributions in nuclear resonance scattering 28,29 .
The wavelength of the VUV radiation emitted by the 229 Th isomer was determined by measuring the transmission of this radiation through six different bandpass filters (BPFs).The BPFs were switched every 10 seconds with a switching time of less than two seconds during the 1,800 s measurement period following an X-ray pumping isomer excitation.Measuring all six BPF transmissions during single excitation-decay sequence suppresses possible detrimental effects of beam intensity fluctuations and crystal damage on the wavelength measurement.
Figure 3 shows the determination process of the radiation wavelength from the thorium isomer.Each measured transmission T meas (i) (i = 1 − 6) was obtained (Fig. 3(a)) by taking the signal ratio with and without the filters.The filter number #0 corresponding to a blank one (without the filter) were used once in one rotating cycle for this T meas (i)-calculation.The transmission spectra T BPF (i, λ) for the BPFs (Fig. 3(b)) were measured before and after the whole measurement campaign using a custom-built measurement system with a VUVmonochromator.
The wavelength of the emitted radiation is obtained by fitting the parameter λ to minimize the weighted quadratic sum of the deviations between T meas (i) and T BPF (i, λ), describing as From this procedure, the wavelength is determined to be 148.18± 0.38 (stat.)nm as shown in Fig. 3(c) and (d).The systematic uncertainties for T meas (i) and T BPF (i, λ) were estimated by off-line systematic studies of BPF transmission, wavelength calibration, and their reproducibility (see Appendices).Considering these systematic effects, the wavelength of the 229 Th isomer transition is determined as This value corresponds to an isomer energy of 8.367±0.024eV, and is consistent with recently reported energies within the error [8][9][10]30 . The xcitation dynamics can be observed by measuring the total isomer signal yield immediately after X-ray beam-off for various excitation periods (see Fig. 4 (a)).The subsequent radiative decay dynamics can be directly extracted from PMT time traces recorded after termination of the excitation beam (see Fig. 4 (b)).To eliminate effects of X-ray-induced crystal luminescence, the isomer signal is obtained by subtracting the temporal profile of the VUV-signal after the excitation period at two different beam energies: one on the resonance peak of Fig. 1(g) and one far from (about 0.10 eV) the resonance peak, called "on-resonance" and "off-resonance", respectively.Such a subtracted decay signal after the excitation period of the Basso crystal exhibits a clear exponential decay as expected from a nuclear process (Fig. 4(b)).The least squares fitting yields a half-life of T 1/2 = 431 ± 23 s.Taking into account several types of systematic uncertainty (see Appendices), the result is T 1/2 = 431 ± 23(stat.) ± 26(syst.)s.
The same experiment and analysis were performed using the Alto crystal to check consistency between samples with similar properties.The obtained half-life is T 1/2 = 459 ± 21(stat.) ± 23(syst.)s, which is consistent with the Basso-crystal within the error.The half-life value from the combined data of both crystals is presented as the final result: T 1/2 = 447 ± 16(stat.) ± 20(syst.)s. (3) We observe no change in the measured isomer half-life using crystals with different doping concentrations, indicating that the effects of re-absorption of the emitted VUV isomer photons ("light-trapping") can be neglected (see Appendices).We also observe no systematic change of the isomer decay half-life T 1/2 with parameters of the X-ray excitation beam; only the total amount of isomer state population is affected by the X-ray energy, photon flux, and spectral purity.This is in strong contrast to the excitation dynamics, where we clearly observe a scaling of the isomer population timescale with X-ray flux. Figure 4 shows an excitation-decay dynamics for an X-ray beam flux of 2 × 10 11 photons • s −1 .
The isomeric state population after the excitation time T ir is fitted by a build-up function T a factor of about ten times shorter than the isomer decay halflife T 1/2 .This indicates the presence of an additional isomer decay channel during the exposure of the Th-doped crystal to the X-ray source, which we refer to as "X-ray quenching".We have investigated the scaling of this effect with X-ray intensity and found an approximately linear scaling of the quenching factor T 1/2 /T (ir) 1/2 with X-ray beam flux as shown in Fig. 5.These observations raise the question on the probability of emitting a VUV photon in a single decay event from the 229 Th isomer back to the ground state (in the absence of the Xray beam).The 8 eV nuclear-excited state in the crystal may interact with the surrounding chemical environment.As a result, the isomer state has the possibility to decay to the ground state through non-radiative processes.
The radiative transition probability B (is) rad is estimated by comparing the experimentally measured number of VUVphotons with the one calculated, based on the production yield of 229m Th and the experimental detection efficiency (see Appendices).
For the parameters of Fig. 4, an excitation time of 600 s, an effective half-life T (ir) 1/2 = 39.2 s, we expect to excite in total (3.0-3.5)×10 5 nuclei to the isomeric state which will decay back to the ground state with the half-life T 1/2 = 447 s.Integrating the detected VUV count rate over 0-1800 seconds and dividing with the detection efficiency and trigger efficiency (ϵ det = 6.5 × 10 −3 and ϵ trig = 0.845) (see Appendices) yields a number of 3 × 10 5 emitters.
We conclude that the probability of radiative decay B (is) rad is consistent with one, the lower limit (90 % confidence level) for this probability is estimated to be 0.45 (0.37) for the Basso (Alto)-crystal, respectively.

Discussion and Conclusion
Our findings indicate that, in the absence of the X-ray beam, radiative decay under the emission of a VUV photon is the dominant de-excitation process for the 229 Th isomer in the crystal.However, more systematic studies are needed to clarify whether weak non-radiative decay channels exist and under which experimental conditions.The long half-life of T 1/2 = 447 s measured after the excitation period is considered to be close to the half-life of the radiative decay, indicating the predicted narrow transition linewidth.Since the isomeric transition is expected to be mainly magnetic dipole (M1) 31,32 , the measured half-life in the crystal is corrected by applying a factor of n 3 (ref. 18) to the half-life in vacuum if we assume the condition where internal conversion decay is forbidden.By using n = 1.588 at 148.2 nm for a CaF 2 crystal, the obtained half-life value for an isolated nucleus in vacuum is 1, 790 ± 64(stat.) ± 80(syst.)s, which is consistent with the recently reported isomer half-life of 1, 400 +600 −300 s of the trapped 229m Th 3+ ion 33 .
This value is as long as expected by several theoretical works.In fact, there are several reports for the isomeric half-life of O(10 3 ) s predicted using the microscopic nuclear model 31 , and the simple rotational model of the deformed nucleus using the relevant experimental data 34,35 .More specific studies comparing predicted and measured values are expected in the future.
The half-life obtained from measurements recently performed at CERN-ISOLDE was 607(102) s 8 .Possible reasons for this discrepancy are the effect of different substrates (MgF 2 and CaF 2 ), different environmental conditions (heavy-ion implantation vs. crystal doping), and the effect of different depths of the Th-isomer end-position with respect to the crystal surface.In the present experiment, the pumped 229m Th nuclei are distributed over the entire (illuminated) crystal volume and are expected to interact with a common crystal environment (refractive index).
The significantly reduced (flux-dependent) isomer half-life observed during X-ray beam irradiation suggests the existence of additional decay paths.Possible candidate mechanisms are couplings to electronic near-bandgap or defect states, or electronic bridge processes in the electron cascades following X-ray-induced core-electron ionization.This "isomerquenching" facilitates an artificial acceleration of nuclear isomeric decay, which has been discussed in other research fields, such as triggering the energy release of long-lived isomer states (sometimes called "isomer depletion") [36][37][38] .In atomic clock operation, the atomic state is initialized within milliseconds after measuring the excitation efficiency of the clock transition.However, there is currently no available method to initialize the isomer state to the ground within a short time scale for a solid-state nuclear clock.Identifying and controlling the present decay dynamics opens up the possibility to actively quench the remaining isomer populations, which can be used to speed up the interrogation sequences in future solid-state nuclear clocks based on 229 Th 6 .
After submission of this manuscript, resonant laser excitation of the Thorium-229 isomer was reported in ref. 39 .The excitation energy and radiative decay lifetime are compatible with values reported in this work.We find it interesting to note that the crystal X2 used in ref. 39 was cut from the same ingot as the "Alto" and "Basso" crystals used here.
In conclusion, we have excited the 229 Th nuclear clock isomer from the ground state using resonant X-ray pumping, and observed the subsequent VUV-photon emission of the radiative isomeric decay.We identify direct VUV photon emission as the predominant de-excitation process for 229 Th doped into bulk CaF 2 single crystals.The determined isomer excitation wavelength is consistent with recently reported values [8][9][10]30 . We eplored the temporal profiles of the excitation and decay dynamics and found substantially different half-lives, indicating qualitatively different coupling mechanisms between the isomer and the crystal environment.The X-ray-induced iso-mer quenching effect, as revealed in this study, is anticipated to become a pivotal element for the forthcoming operation of nuclear clocks.

Appendices A Th-doped CaF 2 target
To obtain high doping concentrations exceeding 10 18 cm −3 with the available restricted amounts of 229 Th source material, a miniaturized version to the commonly used vertical gradient freeze (VGM) method was developed by TU-Wien with support by the Fraunhofer institute for integrated systems and device technology (IISB).We grow 3.2 mm diameter, 11 mm long crystals by applying a steep temperature gradient (20 K/cm) across the starting growth material (dopant powder and single crystal seed) and dynamically control the liquid-solid interface.A typical growing speed is (<0.5 mm/h).The crystal growing device is kept under vacuum during the growth process, approximately 10 −4 mbar at the start of the growth.The vacuum prevents oxidation of the used graphite insulation and CaF 2 powder.
Crystal growing starts from a mechanically processed seed single crystals, which is milled down from commercially available bulk crystal (with known crystal orientation) to 3.2 mm cylinders of 11 mm height.From the top, a 2 mm diameter hole is drilled (5 mm) into the crystal to form a pocket for the doping material.
The doping material is 15 mg of 229 ThF 4 :PbF 2 :CaF 2 powder.PbF 2 acts as a scavenger for oxygen removal and as a carrier that facilitates the handling of the minuscule (micrograms) amounts of 229 ThF 4 during the wet chemistry preparation. 229Th is obtained as dried nitrate from Oak Ridge National Laboratory, an activity of 1.6 MBq was used to grow the crystal used in the reported experiments.It is dissolved in 0.1 M HNO 3 and mixed with PbF 2 and CaF 2 .Subsequently, 229 ThF 4 :PbF 2 :CaF 2 was precipitated by addition of hydrofluoric acid, washed and dried (80 • C, 4 days).
The doping material is transferred into the pocket of the seed crystal using a funnel and inserted into the VGF growing device.Crystal growth itself consists of 5 sections: 1) 18 hours of heating up the system to 800 • C, outgassing, and restoring pressure 2) 6 hours of scavenging oxygen through reaction with PbF 2 at constant temperature 3) 22 hours of melting the top half of the crystal and also slowly freezing it, reaching melting temperature around 1418 • C. The temperature (gradient) has to be controlled to 1 • C (0.1 • C/cm) level, respectively, to adjust the melting depth into the crystal to 1 mm precision 4) 18 hours of annealing the crystal at 1200 • C 5) 14 hours of cooling down.A vacuum of at least 10 −4 mbar is obtained before growth.During growth (especially during section 1) the pressure can go up to 10 −2 mbar.The complete growth process typically takes 3 days.
After growth, the crystal is cut to the required dimensions using a diamond-wire saw (80 microns) and polished using standard polishing disks and slurry.Details of the crystal fabrication and characterization of the optical properties can be found in ref. 17 To improve the VUV transmission, the crystal is heated again to ≈ 1150 • C in a separate device under CF 4 atmosphere at ambient pressure to counteract a fluoride deficiency that occurs due to radiolysis and losses during growth.

B DAQ scheme and background rejection
The data used in this work are derived from two voltage-time traces, simultaneously registered from the VUV-PMT and the VETO-PMT, recorded by an oscilloscope, that is triggered by an event on the VUV-PMT.During measurement, most of the triggering events are radioluminescence-induced background: 98.0 -99.9 % depending on the elapsed time from the stop of beam irradiation.Most of the radioluminescence background can be rejected using waveform information of the VETO-PMT.Basically, pulse height and timing information of radioluminescence are used for background rejection.
An example of a waveform of radioluminescence background events is shown in Fig. 6.At the beginning of the measure- ment, strong photoluminescence signals were measured that caused a shift from the ground level.This shift size depends on the crystal used, the elapsed time from the stop of beam irradiation, and the reduction of the crystal transmittance due to beam irradiation.Therefore, the radioluminescence rejection condition depends on each measurement and the elapsed time in the measurement.
The threshold for the VUV-PMT signal was set to -5 mV, the minimum value at which it can operate effectively.For a quantitative estimate of the obtained VUV isomer photon yield, it is important to evaluate the data loss resulting from the failure to capture the distribution of signals with pulse height less than -5 mV.
To determine the trigger efficiency, an analysis was carried out on the decrease of the pulse-height distribution of the trigger signal from near-threshold to zero.By performing this analysis with different threshold settings, we estimated the number of low pulse-height signals (-10 -0 mV) that were missed throughout the DAQ procedure at the threshold of -5 mV.The efficiency was then estimated to be ϵ trig = 84.5(14) %.The uncertainty in this estimation comes from the minimum edge value and the fitting error of the distribution function in the measured pulse height distribution.
The data acquisition efficiency ϵ DAQ (or dead time) was also estimated because the oscilloscope may not record all triggered data due to accidental coincidence with another event or oscilloscope deadtime.The clock pulse (frequency of 10 Hz or 100 Hz depending on the run condition) of the delay generator was used as the additional trigger for the data recording in order to estimate the DAQ efficiency.The measured timeaveraged efficiencies ranged from 98.31 to 99.95 % in various runs.
Another important source of possible data loss is due to the rejection of the true VUV signal by the VETO-cut procedure.This VETO-cut efficiency ϵ veto corresponding to the survival probability through the VETO-cut selection is estimated by the reduction of the clock trigger events by the VETO-cut selection.The efficiencies ranging from 27.0 to 32.7 % were measured in different runs.

C Systematic uncertainties in lifetime measurements of 229m Th
As described in the main text, the half-life of the isomeric state after the X-ray beam termination was determined by fitting the temporal spectrum (Fig. 4(b)) with the single exponential decay function.Three types of systematic uncertainty were considered in the analysis.One is the uncertainty in the constant offset parameter of the exponential decay function in the fitting procedure.This type of half-life uncertainty was evaluated to be 16 s. in the analysis of the spectrum Fig. 4(b).
Another one is the uncertainty arising from the choice of cut parameters for pulse-height selection in the VETO-PMT signals.This uncertainty was estimated to be 12 s. in the same spectrum.The last systematic uncertainty was estimated from the choice of the fitting time window.The possible maximum deviation was 16 s.if we use the fitting time window 200-1,780 s. instead of 0-1,800 s.We also included the 16 seconds for the systematic error.By using the quadratic sum of the above three systematic uncertainties, the result of the half-life value for the Basso crystal is T 1/2 = 431 ± 23(stat.) ± 26(syst.)s.
The same evaluation of the systematic error was applied to other crystal data.

D Lifetime measurement using low 229 Th density crystal -reabsorption effect -
If the absorption length in the crystal of the on-resonance VUV-photon from the isomer state is smaller than the crystal length, so-called radiation trapping occurs due to reabsorption of the photon by the ground state of other 229 Th nuclei 6 .Such a situation tends to prolong the observed halflife for the isomeric radiative decay.The efficiency of reabsorption (or reabsorption length) depends on the inhomogeneous broadening for the isomeric transition in the crystal, which is currently unknown.However, by comparing the measured radiative half-life of the isomer state for different 229 Th doping concentrations, the effect of such reabsorption can be investigated.We additionally measured the isomer half-life with a low-density Th:CaF 2 crystal of 7 × 10 17 cm −3 , which is one order of magnitude lower concentration compared the Basso/Alto-crystals used in the main result.The measured half-life is 471 ± 67(stat.) ± 59(syst.)s.This value is consistent with the measurement using the high 229 Th density crystal (Basso/Alto).This means we did not observe any reabsorption effect in the radiative decay of 229m Th by other 229 Th nuclei.

E Estimation of detection rate for VUVphoton -detection efficiencies -
In the present scheme of the radiative decay measurement from 229m Th, the 229 Th-doped CaF 2 target crystal is irradiated by the X-ray beam for a time period of T ir .The number of 229 Th isomer nuclei excited at the end of the excitation period T ir is calculated as where τ ir is the lifetime of the isomeric state inside the crystal target during beam irradiation (τ ir = T (ir) 1/2 / ln 2).R is is the production rate of 229m Th by the X-ray beam irradiation and is obtained by where R 2nd is the excitation rate to the second excited state, Br in tot is the branching ratio to the isomeric state ("in"-band transition in the rotational band structure of the 229 Th nucleus) including γ-decay and internal conversion.The parameter ϵ shift indicates the reduction of excitation cross section due to a shift from the on-resonant energy.The parameters of n Th , Φ, l x , and L are the number column density of 229 Th inside the crystal, the X-ray beam flux irradiating the crystals, the attenuation length of the X-ray beam in the crystal, and the crystal thickness, respectively.The n Th was estimated by measuring 193 keV γ-ray intensity following the α-decay of 229 Th by using a germanium detector.The effective cross section from the ground state to the second excited state of 229 Th, σ eff , is written as where λ 2nd , Γ cr γ are the wavelength and the cross-band transition width between the ground state and the second excited state of 229 Th, respectively.The σ x is the RMS energy width of the X-ray beam.The excitation width Γ cr γ and its uncertainty are described in ref. 20 .The branching ratio Br in tot is also described in this article, but the updated value in the ref 40 is used in the present analysis.The uncertainty of the branching ratio is estimated with the uncertainties of the relevant measured quantities.These parameter values and their uncertainties are summarized in Tab. 1.
where each parameter in this equation is explained in the following.The time origin t = 0 is defined as the timing when the target holder was moved in front of the focusing mirror.The parameter τ is the lifetime of the isomeric state inside the crystal in the absence of X-ray irradiation.The parameter B (is) rad is the probability of photon emission from 229m Th in the absence of X-ray irradiation, the efficiencies ϵ trig , ϵ DAQ and ϵ veto in the data processing are explained in section D.2.The parameter ϵ det is the total detection efficiency of the experimental apparatus, which is written as a product of several component efficiencies as where each term corresponds to the VUV-photon transmission of the crystal, geometrical efficiency, reflectance of the parabolic mirror, the reflectance of four dichroic mirror prisms, the transmission of the MgF 2 lens, and the quantum efficiency of the VUV-PMT, respectively.The ϵ crys , ϵ mir , ϵ prism , and ϵ lens were measured by using a deuterium light source and a VUV spectrometer.The geometrical efficiency, which was estimated using optical simulations, depends on ϵ crys and ϵ lens .Therefore, it is appropriate to evaluate these three detection efficiencies as a single term, which was estimated to be ϵ crys ϵ geom ϵ lens = 0.071±0.019.The remaining three efficiencies were estimated to be ϵ mir = 0.89 ± 0, 04, ϵ prism = 0.45 ± 0.02, and ϵ PMT = 0.22 ± 0.01.The parameter C is the correction factor describing the reduction in crystal transmission by beam-induced damage and target holder movement time.The reduction factor due to crystal damage was estimated by observing the decrease in radioluminescence-induced VUV-event rates over several measurement cycles and scaling this decrease to an irradiation time of 600 s, resulting in a value of 0.997 per one excitation period of 600 s.The above efficiency parameters are estimated at the wavelength of 148 nm.The estimated wavelength-dependence of the total detection efficiency ϵ det is plotted in Fig. 7.

F Wavelength determination
The transmission spectra of each BPF and the dichroicmirror set were measured with the custom-made setup consisting of a deuterium lamp (Heraeus, D200VUV), a VUVmonochromator (Vacuum & Optical Instruments, VMK-200-II), two collimators, and a photomultiplier (Hamamatsu, R6836).The systematic error in the wavelength determination arises from uncertainties in wavelength calibration for the monochromator.The absolute wavelength calibration of our assessment setup was performed using the discharge light source with the nitrogen and oxygen gas mixture.We used the several atomic transitions of O(I) and N(I) in the vicinity of 150 nm as listed in Tab. 2 for the calibration procedure.The uncertainty includes the error of the calibration itself and the effect of off-axis mounting of the light sources (deuterium lamp and N 2 /O 2 gas discharge).The measured systematic uncertainty in the wavelength calibration was conservatively estimated to be 0.1 nm.The wavelength of the radiation from the 229m Th was coarsely selected by the dichroic-mirror set placed upstream of the BPFs.This wavelength selection was investigated using different sets of dichroic mirrors without the BPFs. Figure 8 a shows the reflectance of three different dichroic mirror sets.For each dichroic mirror set, the signal yield was measured by changing the X-ray beam energy.The wavelength λ is is determined by minimizing the following χ 2 : where C and λ is are the two independent variables for the minimization.The quantities N data (i)/Φ(i) is the number of events measured in the dichroic mirror set i, which is normalized by beam flux density Φ(i).∆(N data (i)/Φ(i)) is its uncertainty.From this measurement, the wavelength is determined as λ is = 148.14± 0.40(stat.)± 0.53(syst.)nm.(12)   This wavelength value is consistent with the λ is measurement using band-pass filters.The wavelength of the radiation from the 229m Th is finally determined by comparing the transmission of this radiation through the six different BPFs and their transmission spectra as described in the main text.Therefore, the systematic uncertainties in wavelength determination come from both measurements of transmitted isomeric radiation and transmission spectra of BPFs.
One of the systematic uncertainties was determined by calculating the difference in wavelength between two spectra measured before and after the VUV search experiment.This uncertainty was found to be 0.05 nm.The uncertainty in the measurement of BPF-transmission spectra, which includes the position dependence inside each BPF and the measurement reproducibility, produces another systematic error in wavelength determination.This type of systematic uncertainty was measured to be 0.07 nm.Systematic uncertainty also comes from the analyzing procedure in the T meas (i)-determination, such as the selection of a time window in the radiative decay signal and how to integrate the multiple measurement runs.These uncertainties were estimated to be 0.09 nm and 0.15 nm, respectively.
By taking the quadratic sum of the above several systematic errors, we concluded that the final result of the measured wavelength with the BPF system is λ is = 148.18± 0.38(stat.)± 0.19(syst.)nm, (13)   as described in the main text.

Fig. 1 :
Fig. 1: Experimental setup.a, The nuclear levels of 229 Th and the scheme of isomer excitation and the observable signals (X-ray and VUV signals).b, Overview of the system for isomer excitation and the detection of radiative decay.c, Setup around the wavelength filtering devices inside the detection chamber; the moving target holder, the parabolic mirror, the prism set, the MgF 2 lens, the band-pass filters, the VUV-PMT and the VETO-PMT.d, The four consecutive prisms with dichroic mirror coating inside the wavelength filtering device.e, Enlarged view of the thorium-doped crystal and the crystal holder.The crystal was fixed by using stainless wires.

Fig. 2 :
Fig. 2: Measured resonance spectra.a, The resonance spectrum of X-ray fluorescence (NRS signals) obtained with the thorium nitrate target.b, The resonance spectrum of VUV-photon signals obtained with the 229 Th:CaF 2 crystal.The horizontal axes of a and b indicate the absolute X-ray beam energy (offset by 29,189 eV).The vertical axis is the number of PMT counts in 1,800 seconds after a beam irradiation time of approximately 600 seconds.The background components (the constant offset in each Gaussian) show the beam-induced luminescence and radioluminescence from each target.

Fig. 3 :
Fig.3: Determination of the isomer transition wavelength.a, The measured transmissions of the radiation from the thorium isomer through the band-pass filters (BPFs), expressed as T meas (i) (i =1-6).b, The transmission spectra T BPF (i, λ) for each BPF.c, The six overlaid plots between T BPF (i, λ) and T meas (i) for each BPF, where the T meas (i) are located at the best-fit wavelength (see the text).d, The calculated χ 2 distribution defined as eq.(1).

Fig. 4 :
Fig. 4: Time evolution of 229m Th production and its decay.a, The observed time evolution of isomer production.The data points show the observed VUV counts at the time of zero in the decay signal at each excitation period (X-ray beam irradiation time).The solid black curve indicates the result of fitting with a build-up function (see the text).The dashed blue line represents the expected time evolution if T (ir) 1/2 equals to T 1/2 .Here, we scaled the amplitude of this curve to match the fitted function at infinite irradiation time.b, The observed temporal profile of the isomeric decay signal after beam irradiation time of 600 s.The data points show the observed VUV counts at each elapsed time after the excitation beam was switched off.The blue solid line represents the result of fitting with an exponential decay function.The bottom row in (a) and (b) shows the corresponding residual plots for each fitting procedure.The error bars in these plots show only statistical errors.

Fig. 5 :
Fig. 5: Scaling of quenching factor T 1/2 /T (ir) 1/2 with Xray beam flux.The horizontal axis indicates the weighted average beam intensity in photons per second.The inset shows an enlarged view near the zero-flux region.The green diamond-shaped point indicates the quenching factor at zero flux, which is equal to one.The point with the flux of 18.6×10 10 /s corresponds to the excitation dynamics depicted in Fig. 4(b).The error bars represent the statistical uncertainty of one standard deviation.The solid line indicates a fit with a linear function where the intercept is fixed at one.

Fig. 6 :
Fig.6: An example of the waveform of the oscilloscope.Time origin t = 0 indicates the timing of the triggered event from the VUV-PMT.Some strong peaks in the waveform of the VETO-PMT are caused by radioluminescence.a, Waveform taken at the beginning of one measurement.b, Waveform taken at the end of one measurement.

Fig. 7 :
Fig.7: The wavelength dependence of estimated total detection efficiency of the isomeric radiation.

Fig. 8 :
Fig.8: Determination of wavelength using different dichroic mirrors.a, Reflectance of different dichroic mirror sets used for determination of the wavelength of radiative decay photon from 229m Th. b, The obtained energy spectra of VUV-signal with the three different dichroic mirror sets.Each point corresponds to a measurement result with a beam irradiation period of 600 seconds.Solid lines show the simultaneous result of a Gaussian fit with a constant background of the three spectra, where center and width are the common parameters.Color of each line in a corresponds to that in b.

Table 1 :
The parameter values and their uncertainties for each component.The values with a range mean the run-dependent values.Some uncertainty parameters also depend on experimental conditions.These are used for the estimation of detectable VUV-photon number from the 229m Th.

Table 2 :
41st of transition wavelengths used for calibration of the transmittance measurement system.The wavelength values are taken from NIST Atomic Spectra Database.41