\DeclareSourcemap\maps

[datatype=bibtex] \map \step[fieldset=title, null]

Two-exponential decay of Acridine Orange

Francesco Giacosa corresponding author: [email protected] Jan Kochanowski University, ul. Żeromskiego 5, 25-369 Kielce, Poland Goethe University, Theodor-W.-Adorno-Platz 1, 60323, Frankfurt am Main, Germany Anna Kolbus Jan Kochanowski University, ul. Żeromskiego 5, 25-369 Kielce, Poland Krzysztof Kyziol Jan Kochanowski University, ul. Żeromskiego 5, 25-369 Kielce, Poland Magdalena Plodowska Jan Kochanowski University, ul. Żeromskiego 5, 25-369 Kielce, Poland Milena Piotrowska Jan Kochanowski University, ul. Żeromskiego 5, 25-369 Kielce, Poland Karol Szary Jan Kochanowski University, ul. Żeromskiego 5, 25-369 Kielce, Poland Arthur Vereijken Jan Kochanowski University, ul. Żeromskiego 5, 25-369 Kielce, Poland
(November 6, 2025)
Abstract

In this work, we experimentally study the fluorescence decay of Acridine Orange at late times, in order to test whether a late-time power-law behaviour emerges — a feature expected to be very small but consistent with quantum mechanical and quantum field theoretical predictions. Using two distinct photon detectors, we find that the data are well described by a sum of two exponential functions with lifetimes τ1=1.7331±0.001\tau_{1}=1.7331\pm 0.001 ns and τ2=5.948±0.012\tau_{2}=5.948\pm 0.012 ns, in agreement with values reported in the literature. While no deviation from the exponential decay law is observed, this study serves as a reliable test for the experimental setup and enables a precise determination of the sample lifetimes.

keywords:
quantum mechanics, fluorescence, time-correlated single photon counting spectroscopy, lifetime

Introduction: The standard radioactive decay law describes the change in number of elements in a sample as N(t)=N0et/τ,N(t)=N_{0}\,e^{-t/\tau}, an exponential function. However, this successful phenomenological formulation is not strictly reproduced within the framework of Quantum Mechanics (QM) and Quantum Field Theory (QFT) according to which N(t)=N0P(t)N(t)=N_{0}P(t), where P(t)P(t) is the survival probability as a function of time tt of a single unstable quantum state. The quantity P(t)P(t) emerges as the modulus squared of the Fourier transform of the energy distribution ρ(E)\rho(E) of the unstable quantum state/particle, e.g. [LFonda_1978]:

P(t)=|EthdEρ(E)eiEt|2 ,P(t)=\left|\int_{E_{th}}^{\infty}\mathrm{d}E\,\rho(E)e^{-i\frac{E}{\hbar}t}\right|^{2}\;\text{ ,} (1)

where EthE_{th} is the lowest admissible energy. Only for EthE_{th}\rightarrow-\infty and ρ(E)\rho(E) being a Breit-Wigner distribution, the decay is purely exponential. If, instead, the expectation value of the Hamiltonian operator HH is finite, then P(0)=0P^{\prime}(0)=0. If also the H2H^{2}-expectation is such, P(t)P(t) can be approximated at short times as P(t)1t2/tZ2P(t)\simeq 1-t^{2}/t_{Z}^{2} with τZ=/σE\tau_{Z}=\hbar/\sigma_{E}, where σE\sigma_{E} is the standard deviation of the energy distribution. This deviation from the exponential decay law at short times is linked to the quantum Zeno effect (QZE), i.e., the freezing in its initial state if probed sufficiently often [Misra:1976by, 1998PhRvA..57.1509S, Facchi:2008nrb, 2002OptCo.211..235B, Kofman:2000gle]. On the other hand, at large times the survival probability P(t)P(t) is well described by a power law, P(t)t(β1)P(t)\sim t^{-(\beta-1)}, leading to the decay rate intensity I(t)tβI(t)\sim t^{-\beta}. This is a direct consequence of the presence of the energy threshold EthE_{th} [khalfin]. Both mentioned effects are predicted in the framework of QM [LFonda_1978, dicus2002, Urbanowski1994] and also quantum field theory (QFT) [Giacosa:2011xa, Giacosa:2010br, Giacosa:2021hgl, Facchi:1999ik] and emerge from their very first principles. However, the deviations are typically very small for elementary systems/transitions, as the ‘parade’ example of the 2P2P-1S1S transition of the HH-atom shows: short-time deviations take place at time 108τ10^{-8}\tau and late-time deviations for 125τ\sim 125\tau (τ=1.595\tau=1.595 ns), thus very early and very late, rendering an observation extremely challenging [Facchi:1998abc, Giacosa:2024yhn].

At short times, the experimental verification of deviations from the exponential law was measured in the study of quantum tunneling of sodium atoms within an optical potential [1997Natur.387..575W, 2001PhRvL..87d0402F]. (For indirect evidence, see the photonic-waveguide array experiment reported in Ref. [crespi]; strongly decay systems also display deviations from the Breit-Wigner function [ALEPH:2005qgp, Giacosa:2010br], in turn implying a nonexponential behaviour which, however, cannot be directly measured due to the extremely short times (1023\sim 10^{-23} s) involved.)

For long times, the well-known work of Ref. [rothe] reports a power law intensity I(t)I(t) for various chemical compounds decaying via fluorescence [lakowicz]. Very recently, a power law was confirmed by our group for the case of erythrosine B upon using two distinct photodetectors [Giacosa:2025knr]. In both cases, dissolved fluorophores were employed; such systems, even if not as simple as ‘more elementary quantum states’ because of inhomogeneous broadening of the emission spectrum resulting from interaction of the compound with the solvent, are thought to still be described by the inherent quantum Eq. (1).

In this work, we report on the late-time study of a different substance, acridine orange, measured using two photon detectors acting in different ranges. As we shall show, the decay rate can be well described by the sum of two exponential functions, in agreement with a mixture of two quantum states with two distinct lifetimes. The rather precise determination of lifetime(s) is in agreement with values reported in the literature. Moreover, even if no ‘quantum’ late-time memory effect can be seen, this study offers a valuable test of our experimental approach aimed at investigating the late-time decay law.

Basic features of fluorescence. Excited states of fluorophores provide typical examples of unstable quantum states. These molecules may absorb a photon, that implies an excitation of an electron to higher energy levels denoted by S1,S2,S3S_{1},S_{2},S_{3}, … . Then, this electron undergoes a series of rapid transitions via a non-radiative process referred to as vibrational relaxation. After a short time (1012s\sim 10^{-12}\;\mathrm{s}), the electron is in the lowest vibrational level of the excited electronic state. Subsequently, it may decay into lower electronic states, usually via internal conversion (IC) – which is also a non-radiative process – or via the radiative fluorescence process, which we detect. This transition usually occurs from the lowest vibrational level of the first excited electronic state (S1S_{1}) to one of many vibrational levels of the electronic ground state (S0S_{0}). According to the Franck-Condon rule, the most probable transitions are those for which the initial and final electronic wave functions overlap the most. Typically, the dynamics of the fluorescence takes place in a timescale of the order of nanoseconds.

Description of the device. In the performed experiment, the Time-Correlated Single Photon Counting (TCSPC) setup available at the Institute of Biology, Faculty of Natural Sciences, Jan Kochanowski University in Kielce (Poland) was used. It consists of PicoQuant Laser Combining Unit (LCU), Nikon Eclipse Ti-E Inverted Confocal Microscope, detection system, and PicoQuant PicoHarp 300 TCSPC module. The laser combining unit is composed of two picosecond laser diodes – PicoQuant LDH-D-C-440 and PicoQuant LDH-D-C-485 emitting at, respectively, 438 nm and 485 nm with spectral bandwidth between 2 and 8 nm. However, only the latter one was actually used during the measurements. Laser diodes serve as pulsed excitation sources. They were operating at 10 MHz frequency (100 ns interval between pulses) and the width of the pulse at half maximum (FWHM) is less than 120 ps.

The detection setup is composed of two identical PicoQuant PMA Hybrid 40 detectors. They are separated with the dichroic mirror and every detector is coupled to an individual bandpass filter – FF01-520/35 and ET600/50. Thus, two independent detection channels can be distinguished that focus at different parts of the spectrum – Channel 1 detects photons of wavelength ranging from 485 to 555 nm and Channel 2 photons from 550 to 650 nm. Fluorescent samples are placed on the stage of the Nikon Eclipse Ti-E Inverted Confocal Microscope. Thus, the microscope plays an important role in the formation of the excitation and emission beams and their guidance. In our case, it is equipped with one optical fibre transferring the laser pulse from the LCU to the device and another one extracting the fluorescence signal towards the detection system. The last important piece is the PicoHarp 300 TCSPC module which is responsible for timing events, Analog-to-Digital conversion of the signal, and, finally, assigning particular events to appropriate time bins of the histogram.

Description of the experiment. The choice of acridine orange was motivated by the shape of the absorption spectrum that matches the available excitation wavelengths of the LCU, reasonable quantum yields, and its common use in fluorescence studies, offering a valuable test of our device.

As a first step, the fluorophore was dissolved in water reaching a concentration equal to 105mol/dm310^{-5}\;\mathrm{mol/dm^{3}}. After transferring the sample to the confocal microscope, one real-time test measurement was performed to locate the maximum signal-to-noise ratio. This initial procedure was followed by three 10-minute measurements, each of them made on a different part of the sample. The geometrical conditions remained similar throughout the process. Thus, the obtained fluorescence intensity decay curves could be combined into one final data set for an individual detection channel.

Results for Acridine Orange. In the process of analysis, the functions describing the decay rate were obtained by a χ2\chi^{2}-approach. The standard deviation was estimated with Poisson statistics. Thus, for a given rate I(t)I(t), the χ2\chi^{2} reads:

χ2=n=1Nt[InI(tn)]2I(tn),\chi^{2}=\sum_{n=1}^{N_{t}}\frac{\left[I_{n}-I(t_{n})\right]^{2}}{I(t_{n})}\;, (2)

where NtN_{t} is the total number of data points. A single exponential function was tested, but it cannot describe the data. Hence, we employed two models – a sum of two exponential functions and a nonexponential late-time power law. In the first model, it is assumed that there are two types of initial excited states which decay exponentially, characterised by different lifetimes (τ1\tau_{1} and τ2\tau_{2}). The second encodes a late-time QM-inspired power law.

Table I: Model functions employed in the analysis; bb is the background. The time t0=2.24t_{0}=2.24 ns corresponds (for both channels) to the maximum of the intensity (and is not a fit parameter).
Model Fluorescence intensity I(t)I(t) Fit parameters
Two-exponential I(t)=C1exp(tt0τ1)+C2exp(tt0τ2)+bI(t)=C_{1}\exp\left(-\frac{t-t_{0}}{\tau_{1}}\right)+C_{2}\exp\left(-\frac{t-t_{0}}{\tau_{2}}\right)+b χ2(C1,τ1,C2,τ2,b)\chi^{2}(C_{1},\tau_{1},C_{2},\tau_{2},b)
Nonexponential I(t)=Cexp(tt0τ)+Cp(tt0)β+bI(t)=C\exp\left(-\frac{t-t_{0}}{\tau}\right)+C_{p}\,(t-t_{0})^{-\beta}+b χ2(C0,τ,Cp,β,b)\chi^{2}(C_{0},\tau,C_{p},\beta,b)

For completeness, we recall that the fluorescence intensity I(t)I(t) is related to the survival probability P(t)P(t) of Eq. (1) as

I(t)=dNdt=N0P(t).I(t)=-\frac{\mathrm{d}N}{\mathrm{d}t}=-N_{0}\,P^{\prime}(t)\;. (3)

.

Table II: Fitting results for Acridine Orange measurements.
Fitting Range: 3.200 – 96.968 ns
Two-exponential model
Channel χν2\chi_{\nu}^{2} C1C_{1} τ1\tau_{1} [ns] C2C_{2} τ2\tau_{2} [ns] bb
1 1.0471 278 967 1.7333 6 371.3 5.9459 21.99
2 1.0247 150 898 1.7326 8 983.5 5.9493 42.07
Nonexponential model
Channel χν2\chi_{\nu}^{2} CC τ\tau [ns] CpC_{p} β\beta [ns] bb
1 11.0881 220 113 1.9360 39 458 1.7386 -1.404
2 14.7270 93 100 2.2649 40 747 2.421 2.421
Refer to caption
Refer to caption
Figure 1: Fluorescence intensity for both photon detectors (channel 1: left, channel 2: right) - comparison between data and two-exponential fitting function.

Discussions. The values of the fit parameters that minimise χ2\chi^{2} were found numerically111Since the excitation is not instantaneous, we start the fit when the early-time effects are not relevant. In the case of the nonexponential model, this is also required by the fact that the power function is not a valid approximation in the early-time domain., see the results in Tab. II. The following comments are in order:

  • For acridine orange the two-exponential model works well, see the direct comparison of the data and the fit in Fig. 1. The two-exponential behaviour of the decay was previously observed in aqueous solutions of Acridine Orange in sodium dodecyl sulfate (SDS) as a consequence of formation of the micellar structures due to interaction with the detergent – see Refs. [AO1983, AO1988]. The values of both lifetimes are also in the highlighted range (1-2 ns for the main and >3 ns for the other lifetime). This fact suggests that we observe similar aggregates. Interestingly, the deviations from single exponential decay law were also observed earlier in the case of proflavine, which is an acridine’s derivative, bound to DNA [Georghiou] and for acridine orange itself [KUBOTA1977279]. However, in the latter case this effect was caused mainly by the response of the instrument, and the authors observed a single-lifetime exponential decay after deconvolution of the signal.

  • Following a conservative estimate we quote the following lifetime determinations (inverse-variance weighted mean):

    τ1=1.7331±0.001ns , τ2=5.948±0.012ns .\tau_{1}=1.7331\pm 0.001~\text{ns , }\tau_{2}=5.948\pm 0.012~\text{ns .} (4)

    The lifetimes determined by the different detectors are consistent with each other: for Channel 1, the standard deviations are 0.0012ns0.0012~\text{ns} for τ1\tau_{1} and 0.0196ns0.0196~\text{ns} for τ2\tau_{2}; for Channel 2 they are 0.0021ns0.0021~\text{ns} and 0.0156ns0.0156~\text{ns}, respectively. The errors of the fit parameters were taken as the diagonal elements of the covariance matrix (the inverse of the Hesse matrix) [cowan]. As summarized in Table II, the absolute channel-to-channel differences are |Δτ1|=0.0007ns|\Delta\tau_{1}|=0.0007~\text{ns} and |Δτ2|=0.0034ns|\Delta\tau_{2}|=0.0034~\text{ns}, both comfortably below the quoted conservative error bars. A noticeable change is observed only in the C1/C2C_{1}/C_{2} ratio, attributable to the spectral shape and to the position of the spectral window associated with each detection channel.

  • The nonexponential model does not offer an acceptable data description, see Tab. II. Thus, no sign of a QM power-law is seen in the present data for acridine orange, even if we follow the sample up to more than 10τ110\;\tau_{1}. (Note, the eventual appearance of a power law after about 10 lifetimes [rothe, Giacosa:2025knr] implies the existence of an effective threshold EthE_{t}h quite close to the resonance peak [Giacosa:2025knr]. Yet, the exact nature and emergence of this threshold is not fully understood, so at present one cannot predict the turnover time but can only measure it.)

  • In Ref. [Giacosa:2025knr], using the same setup we could measure a power-law for erythrosine B. The comparison of that system with the one of this paper allows us to understand which features favour the detection of the late-time power-law: (i) erythrosine B main lifetime is shorter (τ1=0.45\tau_{1}=0.45 ns), thus we could measure the decay curve I(t)I(t) for a larger t/τ1t/\tau_{1} fraction. In general, a dye lifetime less than 1 ns is definitely preferable for our purposes. (ii) Erythrosine B is dominated by a single lifetime plus a late-time power law. Acridine Orange exhibits two lifetimes originating from a two-state admixture, the longer one dominating at times where a late-time power-law behaviour would possibly emerge.

  • For erythrosine B, a two-exponential fit was tested, but it is definitely worse than the power-law one. Moreover, the two detectors would measure incompatible values for τ2\tau_{2} [Giacosa:2025knr], which makes a two-state interpretation untenable. Indeed, also the power-law coefficients β\beta measured by both channels differ, but this is consistent with basic QM and QFT expectations [Giacosa:2021hgl, Giacosa:2011xa].

Conclusions. In this paper, we have presented decay–intensity measurements of acridine orange, whose fluorescence is well described by the sum of two exponential components, consistently confirmed by both photon detectors. Although no late-time power-law behaviour has been observed, this study provides a valuable validation of our experimental setup for late-time decay measurements, as discussed in Ref. [Giacosa:2025knr], and defines the experimental conditions necessary for future searches of possible quantum-mechanical (and quantum-field-theoretical) deviations from the exponential decay law.

Acknowledgments This work was supported by the Polish Minister of Science under the ‘Regional Excellence Initiative’ program (project RID/SP/0015/2024/01, sub-projects RID/2024/LIDER/08 and RID/2025/LIDER/02).