# Singlet fission in heterogeneous lycopene aggregates

Chloe Magne, Vasyl Veremeienko, Roxanne Bercy, Minh-Huong Ha-Thi, Ana Andreea Arteni, Andrew A Pascal, Mikas Vengris, Thomas Pino, Bruno Robert, Manuel J. Llansola-Portoles

► To cite this version:

Chloe Magne, Vasyl Veremeienko, Roxanne Bercy, Minh-Huong Ha-Thi, Ana Andreea Arteni, et al.: Singlet fission in heterogeneous lycopene aggregates. Scientific Reports, 2025, 15 (1), pp.5593. 10.1038/s41598-025-88220-z. http://dx.doi.org/10.1038/s41598-025-88220-z

HAL Id: hal-04972159

https://hal.science/hal-04972159v2

Submitted on 28 Feb 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN **Singlet fission in heterogeneous lycopene aggregates**

Chloe Magne<fnref n="1" />,⁴, Vasyl Veremeienko<fnref n="1" />,⁴, Roxanne Bercy<fnref n="1" />, Minh-Huong Ha-Thi², Ana A. Arteni<fnref n="1" />, Andrew A. Pascal<fnref n="1" />, Mikas Vengris³, Thomas Pino², Bruno Robert<fnref n="1" /> & Manuel J. Llansola-Portoles<fnref n="1" />

We have prepared lycopene aggregates with low scattering in an acetone-water suspension. The aggregates exhibit highly distorted absorption, extending from the UV up to 568 nm, as a result of strong excitonic interactions. We have investigated the structural organization of these aggregates by resonance Raman and TEM, revealing that the lycopene aggregates are not homogeneous, containing at least four different aggregate species. Transient absorption measurements upon excitation at 355, 515, and 570 nm, to sub-select these different species, reveal significant differences in dynamics between each of the aggregate types. The strong excitonic interactions produce extremely distorted transient electronic signatures, which do not allow an unequivocal identification of the excited states at times shorter than 60 ps. However, these experiments demonstrate that all the lycopene aggregated species form long-living triplets via singlet fission.

**Keywords** Self-assembled lycopene, Singlet exciton fission, Resonance Raman, Transient absorption

Carotenoids are a large family of tetraterpenoid derivatives with a remarkably complex excited-state manifold and dynamics. The lowest-lying excited state, S<sub>1</sub> (2<sup>1</sup>A<sup>-</sup>), is absorption-silent, displaying the same symmetry as the ground state S<sub>0</sub> (1<sup>1</sup>A<sup>-</sup>)<sup>1</sup>; the strong absorption of carotenoids arises from a transition from S<sub>0</sub> to the second excited state, S<sub>2</sub> (1<sup>1</sup>B<sub>u</sub><sup>+</sup>). The energy of the S<sub>0</sub>→S<sub>2</sub> transition depends on the length of the C=C conjugated chain, as well as on the polarisability of the environment<sup>2,3</sup>. In the simplest cases, this excited S<sub>2</sub> state decays by internal conversion (<200 fs) to S<sub>1</sub>, which itself decays to the ground-state S<sub>0</sub> by internal conversion in several picoseconds<sup>4-6</sup>. However, further “dark” states, such as S<sup>*</sup>, are proposed to account for the intricate network of relaxation pathways observed<sup>7-10</sup>. In monomeric carotenoids, the ultrafast deactivation of the excited states prevents the formation of triplets via intersystem crossing. Aggregated carotenoids, on the other hand, are able to generate triplets in an ultrafast manner through singlet fission. Carotenoids constitute the only known family of natural molecules capable of singlet fission in both artificial aggregates<sup>11-18</sup> and biological systems like photosynthetic antennae<sup>19-21</sup> or chromoplasts<sup>22,23</sup>.

The mechanisms governing singlet fission in synthetic compounds such as acenes and rylene is reasonably well understood<sup>24-27</sup>. However, there is a lack of consensus on the mechanisms and even the excited states involved for this process in carotenoids. For zeaxanthin aggregates created by adding water to organic solvent, for example, S<sub>2</sub> has been proposed as the parent state for triplet formation<sup>17</sup>, whereas a hot-S<sub>1</sub> state is proposed for zeaxanthin aggregates in ethanol/THF<sup>16</sup>. Zeaxanthin in lipid membranes forms triplets in a few ps, suggesting S<sub>1</sub>/S<sup>-</sup> as the parent state<sup>28</sup>. For β-carotene aggregates, the triplet parent state is attributed to S<sub>2</sub> for β-carotene micelles<sup>15</sup>, whereas the attribution was inconclusive for β-carotene aggregates in bovine serum<sup>14</sup>. Astaxanthin aggregates in acetone/water mixtures are proposed to have the S<sub>2</sub> state as triplet parent<sup>12</sup>, but no conclusions could be drawn for its aggregates in hydrated dimethyl sulfoxide, despite their absorption spectra being similar<sup>29</sup>. Fucoxanthin aggregates prepared in DMS/water or ethanol/water mixtures exhibit different time evolution with changes of excitation wavelength from 460 to 440 nm<sup>13</sup>. Finally, lycopene aggregates have been proposed to have S<sup>*</sup> as a parent state<sup>11,30</sup> and more recently a vibrational study on the excited states lead to the conclusion that a charge transfer (CT) state is the parent state for singlet fission<sup>31</sup>. The discrepancies are also extended to singlet fission in natural systems, where S<sup>*</sup> has been proposed for long carotenoids in photosynthetic proteins<sup>19-21</sup>, the vibrationally hot-S<sub>1</sub> for daffodil chromoplasts<sup>23</sup>, and an unknown S-like state for tomato chromoplasts<sup>22</sup>.

Treatment of the carotenoid aggregates as a single entity with just one set of spectroscopic properties may be one of the causes of these contradictory results. Changing the preparations conditions (e.g. temperature, solvent mixture, etc...) produces carotenoid aggregates with different properties. However, a few studies have shown

<footnote n="1">¹CEA, CNRS, Institute for Integrative Biology of the Cell (I2BC), Université Paris-Saclay, 91190 Gif-sur-Yvette, France. ²CNRS, Institut Des Sciences Moléculaires d'Orsay, Université Paris-Saclay, 91405 Orsay, France. ³Laser Research Center, Faculty of Physics, Vilnius University, Sauletekio Ave. 10, 10223 Vilnius, Lithuania. ⁴Chloe Magne and Vasyl Veremeienko have contributed equally to this work. ☒email: manuel.llansola@i2bc.paris-saclay.fr</footnote>
that each of the samples obtained using various protocols may contain more than one species of aggregate. These include characterization by resonance Raman of daffodil chromoplasts²³ and transient absorption spectroscopy in astaxanthin aggregates¹² and zeaxanthin aggregates¹⁷. Failure to consider that samples may contain a mixture of different carotenoid aggregates with different properties may lead to incorrect conclusions. Carotenoids are known to change their energetic relaxation pathways dramatically upon minor changes in conformation (e.g. photosynthetic proteins)³²,³³. The diversity of interpretations could also arise from an oversimplified classification of the aggregates. The traditional description for aggregates relies on Kasha theory³⁴, where J- and H-aggregates are assigned on the basis of red or blue absorption shifts, respectively. However, it is currently accepted that this absorption shift is not a strong marker for describing the type of aggregates³⁵. For example, helically-organized H-aggregates of lutein derivatives present a large blue shift for strong excitonic interactions (short molecular distance), and gradually redshift with changes in the molecular distance and twisting angle³⁶. Such red-shifted H-aggregates would previously have been identified (incorrectly) as J-aggregates³⁷,³⁸. Moreover, modelling approaches which take into account vibronic coupling and intermolecular charge transfer have shown that the characterization of aggregates and their excited state manifold is far more complex than the binary H/J determination, namely: Hj-, hJ-, or Herzberg–Teller J/H-aggregates³⁵. The lack of a description of the excited state manifold for carotenoids makes the creation of a model a daunting task.

In this work, we have produced lycopene aggregates in stable acetone-water suspension with low scattering, allowing us to perform spectroscopic measurements in the UV region of the spectrum, and explore the photo-dynamics upon photoexcitation in three different regions. Combined with resonance Raman spectroscopy we highlight the importance of aggregates heterogeneities in the interpretation of carotenoids aggregates photophysics.

## Results

The absorption spectrum of monomeric lycopene in *n*-hexane exhibits three sharp bands, located at 442.4, 470.4 and 501.7 nm, respectively, corresponding to the 0–2, 0–1 and 0–0 vibronic levels of the S<sub>0</sub> → S<sub>2</sub> electronic transition (Fig. 1, solid blue line). These absorption transitions redshift with increasing polarizability of the environment, the 0–0 transition reaching 543 nm in carbon disulfide (Fig. 1, dotted blue line)³. Lycopene aggregates (red curve in Fig. 1) display a very different absorption spectrum, featuring a pronounced band peaking at 355 nm with a shoulder around 375 nm, then two bands at ca 414 and 481 nm, and a very redshifted band at 568 nm. A long absorption tail extending up to 700 nm is also visible. Lycopene aggregates, prepared by different protocols or formed naturally in chromoplasts, lead invariably to similar spectral features, as previously observed in tomato chromoplasts<sup>22,39</sup>, lycopene films<sup>22</sup>, and lycopene crystalloids<sup>39,40</sup>. There is a consensus that these features are due to excitonic interactions<sup>41,42</sup>, but no clear description of the origin of these bands has been proposed. The peak at 355 nm is an unmistakable signature of lycopene forming strong excitonic interactions in H-type aggregates, whereas the bands with small blue or red shifts (481 and 514 nm) and the largely redshifted 568 nm band remain controversial.

Resonance Raman scattering spectroscopy is a highly sensitive tool to investigate multiple species in a complex matrix. In the *resonance* condition, when the energy of the incident light matches the electronic transition of a sub-population in the complex mixture, contributions from this subset of molecules may be selectively observed, and the narrowing of the electronic transitions at low temperatures increases this selectivity<sup>43,44</sup>. This technique

<figure><img src="image_3.png" /><figcaption>**Fig. 1.** Room temperature absorption spectra of monomeric lycopene in *n*-hexane (solid blue line), in CS<sub>2</sub> (dotted blue line), and lycopene aggregates in acetone: water 1:1 suspension (solid red line).</figcaption></figure>
thus allows for the observation of the vibrational properties of different lycopene species in a mixture, provided that they display different absorption properties. Resonance Raman spectra of carotenoids contain four main groups of bands, $v_1$ through $v_4$, each providing specific information on the conformation and/or configuration of the scattering molecule(s)$^2$. In the high-frequency region, the intense $v_1$ band around 1500–1530 cm$^{-1}$ arises from stretching vibrations of the conjugated C=C bonds in the polyene chain. In carotenoid monomers, its precise frequency depends on the effective conjugation length of the molecule$^3$, which is related to the energy of the ground-state. In aggregates, it can be used to probe intermolecular interactions, on the basis of their ground-state properties$^{41}$. The group of bands around 1160 cm$^{-1}$, called $v_2$, arises from stretching vibrations of C–C bonds coupled with C–H in-plane bending modes. These modes constitute a fingerprint region for determining carotenoid isomerization states (*trans/cis*)$^{45–47}$. The $v_4$ mode around 960 cm$^{-1}$ arises from C–H out-of-plane wagging motions coupled with C=C torsional modes$^{48}$, and this region is an indicator of out-of-plane distortions of the carotenoid conjugated chain$^{48,49}$.

The 77 K resonance Raman spectrum of lycopene in diethyl ether (polarizability 0.230) is shown in Fig. 2a. The $v_1$ band is observed at 1517 cm$^{-1}$, a frequency expected for a carotenoid containing eleven conjugated C=C bonds (it may downshift up to 2 cm$^{-1}$ for highly polar solvents). The $v_2$ region displays a main band at 1155 cm$^{-1}$ and is typical for lycopene in the all *trans* configuration. Finally, the $v_4$ region exhibits a small, structure-less envelope characteristic of carotenoid in a relaxed conformation in solution. Figure 2b shows the equivalent spectrum for lycopene aggregates measured using five different excitation wavelengths (363.8, 488.0, 501.7, 514.5, and 577.0 nm) (for room temperature resonance Raman, see supporting information, figure S1). The position of the $v_1$ band of lycopene aggregates exhibits a clear dependence on the excitation wavelength, peaking at 1511, 1529, 1526, 1522 and 1511 cm$^{-1}$ for 363.8, 488.0, 501.7, 514.5, and 577.0 nm excitation, respectively. At 488.0, 501.7 and 514.5 nm, it is accompanied by a small shoulder at 1511 cm$^{-1}$ whose intensity decreases as the excitation moves to the blue. Monomeric lycopene may absorb at 488.0, 501.7 and 514.5 nm, but the position of the corresponding $v_1$ bands (1529, 1526, and 1522 cm$^{-1}$) are shifted to higher frequencies relative to lycopene in ethyl acetate (1517 cm$^{-1}$), indicating that no lycopene monomers are present in the sample$^3$. The $v_2$ region retains the characteristic profile of all-*trans* lycopene, indicating that the aggregation process induces no isomerization of the molecule. The $v_4$ region exhibits sharp vibrational modes at 953, 958, and 966 cm$^{-1}$ for 577.0 nm and 363.8 nm (although 966 cm$^{-1}$ is somewhat obscured by the solvent band for 363.8 nm excitation). The vibrational modes in the $v_4$ region observed for 514.5 nm excitation are significantly lower in intensity, and appear at 950, 955, and 963 cm$^{-1}$. The absorption at 514 nm thus arises from a different lycopene species, for which the packing twists/bends their backbone to a lesser extent and induces a different configuration (different position of the vibrational modes). The $v_4$ regions observed for 488.0 nm and 501.7 nm excitation are similar to carotenoids in solution, indicating the presence of carotenoids with an unstrained conformation in these aggregates.

The resonance Raman data thus reveal the presence of several different lycopene aggregated species. Three different species are responsible for the absorption bands at 488.0, 501.7 and 514.5 nm, referred to below as species 1, 2 and 3, respectively. The species 1 and 2 are in a relaxed conformation, whereas species 3 is somewhat twisted/bent. Assignment of the absorption bands at 355 and 568 nm needs a more detailed analysis. The same $v_1$ frequency (1511 cm$^{-1}$) is observed for excitations in both transitions (363.8 and 577.0 nm, respectively). This could imply that these bands are two different electronic transitions of the same species. While the ratio between

<figure><img src="image_4.png" /><figcaption>**Fig. 2.** Resonance Raman spectra at 77 K of lycopene monomer in diethyl ether (a), and lycopene crystalloids excited at 363.8, 488.0, 501.7, 514.5 and 577.0 nm—asterisks mark acetone solvent bands (b).</figcaption></figure>
the $v_3$ band at 1005 cm⁻¹ and the band at 958 cm⁻¹ appears different between 577.0 and 363.8 nm excitation, this could originate from the different excited states used to induce resonance. These spectra can thus be interpreted in two ways. First, only one strongly-coupled H-type species is present in the mixture (species 4), corresponding to a highly twisted/bent lycopene displaying two electronic transitions at 577.0 and 363.8 nm. Alternatively, two strongly-coupled H-type species are present in the sample (species 5 and 6), absorbing at 355 and 568 nm, respectively, with the former packing the lycopene in a somewhat more relaxed conformation than the latter. The properties of species 1, 2 and 3 (as well as species 6 if there) are consistent with helical organizations of H-type aggregates, which have been shown to produce an absorption redshift depending on the twisting angle³⁶. Species 5, if present, would correspond to pure H-type aggregates with a strong absorption blue-shift. On the other hand, the possibility of species 4, presenting absorption at 577.0 nm and 363.8 nm, is highly puzzling. Although carotenoid H-aggregates have been reported with similar electronic transitions, it was never formally shown that both transitions arose from the same aggregate, as the homogeneity of the aggregate preparation was not verified. We propose that the characteristics of species 4 would correspond to the properties predicted by calculations obtained for medium-strength coupled H-aggregates in carotenes³⁶, where there is a main blue-shifted absorption, along with a weak red-shifted peak arising from Herzberg-Teller coupling (considered as a forbidden transition for H-aggregates)³⁵. Note that from these Raman spectra, we can only determine the *minimum* number of species present in the sample, and cannot exclude the existence of any number of minor species we do not detect.

Transmission electron microscope (TEM) images of lycopene aggregates appear to indicate at least two different aggregate morphologies (Fig. 3)—a crystal-like rod structure (length 2 µm and diameter 100 nm) and a more irregular geometry (diameter 50–200 nm). This is the first time that small, stable lycopene assemblies in solution have been characterized by TEM at high concentration (ca. 10⁴ M). Small rod-like structures have been found for three different H-aggregates of zeaxanthin⁴². It was proposed that the OH groups on the terminal rings of zeaxanthin and lutein facilitate the organization of these aggregates³⁶,⁴². The lack of any such groups in lycopene should thus increase the possibility of different types of association. The observed crystal-like morphologies are disordered, and do not allow a prediction to be made about the symmetry of the unit cell. TEM images do not contain enough information to determine whether each structure contains one single lycopene organization, or there are several sub-domains within the same aggregate.

Femtosecond-to-nanosecond transient absorption data were obtained to investigate the photophysics of the three populations of lycopene aggregate detected by resonance Raman (363.8 nm, 514.5 nm, and 577.0 nm). An overview of the dataset acquired at 355 nm excitation is presented in Fig. 4A, along with kinetic traces measured at selected wavelengths in Fig. 4B. The corresponding data for other excitation wavelengths are available in Supporting Information Fig. S3 and S4. Major differences are observed in the shapes of the kinetic traces and a multitude of different peaks and valleys is observed in the contour plot, both indicating that lycopene aggregates exhibit a highly complex dynamic behavior. Different bands not only appear, disappear and reappear at different timescales, they also exhibit spectral shifts, resulting in kinetic traces with as many as six different minima and maxima occurring at different delay times (see *e.g.* the induced absorption band at 327 nm).

A more systematic data presentation is given in Fig. 5A, where the same three datasets are plotted in the form of time-gated spectra and compared with evolution-associated difference spectra (EADS)—the latter resulting from a global fit of the data using a sequential kinetic scheme.

At first, observation of the datasets collected at three different excitation wavelengths suggests that the main features of the observed dynamics are similar for all three excitation wavelengths. However, significant differences are present in the details of decay kinetics and the intensity ratios of different bands. Let us first classify the different spectral regions by assigning them with names to aid the discussion of the dynamics in detail.

Starting from the high energy side (Fig. 5A,B,C), the difference spectra consist of an induced absorption (IA) band peaking around 330 nm (hereafter referred to as IA-330), and a ground state bleach (GSB) corresponding to the strong 355 nm absorption seen in the steady-state spectrum (denoted GSB-355). Further to the red, IA is again observed as a broad band stretching up to ca. 450 nm (IA-400), from where it is again replaced by a GSB signal exhibiting a clear vibronic-like structure—similar to the series of maxima observed in the red part of the absorption spectrum (GSB-550). Finally, the early signals at wavelengths beyond 580 nm are dominated by IA, and replaced by long-lived negative signals (IA-GSB-600).

A detailed summary of the fs-ns transient absorption dynamics is presented in the form of EADS in Fig. 5D,E,F. The presented spectra are a result of fitting the data to a sequential kinetic model ($A→B→C→...$). The error on the time constants in the ps range is of the order of 25%, and the last constant in the ns range cannot be estimated precisely due to the limited experimental time window. The spectrum of an initial component (lifetime of ca. 150 fs) was strongly contaminated by coherent-coupling artifacts, and was therefore omitted from the presentation. The quality of the fit can be assessed from the traces depicted in Fig. 4B, and the similar data for other excitation wavelengths in supporting Information Figs. S3 and S4. We first note that the evolution time constants produced by the fitting are quite similar for all three datasets, allowing us to directly compare the spectra of different components between the three-excitation wavelengths.

In the IA-300 region, the induced absorption band exhibits a 0.5-ps growth, and then decays almost to zero over the next 1.1 ps evolution step (see transition from black to red in Fig. 5D,E,F). This is in contrast with all the other spectral regions (GSB-355, IA-400, GSB-550 and IA-GSB-600), where the absorption difference decays twofold in 0.5 ps, but in the subsequent 1.1 ps remains approximately constant, albeit with small spectral shifts. These shifts are best noticeable in the GSB-355 region, where the GSB maximum first shifts to the red, and then back to the blue. The spectral evolution on ca. 10-ps timescale produces the most dramatic spectral change (green-to-blue in Fig. 5D,E,F). IA-330 reappears (slightly red-shifted compared to the 1.1-ps spectrum), while GSB-355 shifts noticeably to the red and develops a negative tail, replacing the initially-observed positive
<figure><img src="image_5.png" /><figcaption>Fig. 3. TEM of lycopene aggregates at two different magnifications: scale bar (a) 500 nm, and (b) 1 µm.</figcaption></figure>

IA-400 with negative GSB (the latter is clearly discerned in the blue EADS at wavelengths up to 400 nm). This is accompanied by the disappearance of IA in the IA-GSB-600 spectral region, and the emergence of a negative GSB signal. Within ca. 50–60 ps, the blue EADS is replaced by cyan EADS, where—once again—most of IA-330 is lost, the negative peak of GSB-355 is again shifted to the blue, and IA-400 reappears, reminiscent of the signal observed at the 10 ps evolutionary step. In the GSB-550 region, the vibronic bands shift to the red, and the red GSB tail in the IA-GSB-600 range becomes more pronounced. The last step in spectral evolution is mostly the loss of signal amplitude. However, it also shifts the GSB-355 band to the red (again), loses almost all of the IA signal in the IA-400 region, and further red-shifts the GSB-550 vibronic-like bands and the red GSB tail in IA-GSB-600. This plethora of growths and decays, and shifts in both directions, results in the complicated spectro-temporal picture shown in Fig. 4. Note that unlike lycopene in solution but similarly to lycopene aggregates in films, we are left with long-lived absorption features with bleaching above 580 nm and at 355 nm, with an induced absorption band peaking at ca. 370 nm.

The amplitudes of the different bands exhibit an interesting feature for 570 nm excitation: the signals in GSB-500 and IA-GSB-600 are significantly weaker in this dataset, compared to excitation at 355 and 515 nm. This is obvious in Fig. 5A,B,C, where the signals obtained at 570 nm (C) had to be multiplied by a factor of 5 to be readable on the same scale. This observation is counterintuitive in the context of both steady-state absorption and resonance Raman data, where both 570 and 355 nm bands seemed to be related with each other and different to the species absorbing at 515 nm.
<figure><img src="image_6.png" /><figcaption>**Fig. 4.** Femtosecond-to-nanosecond transient absorption data of lycopene aggregates upon excitation at 355 nm. (A) contour plot of the dataset. (B) kinetic traces measured at different wavelengths (indicated on the graphs) along with the results of the global fit to the data (solid lines). Horizontal dashed lines in panel A indicate the wavelengths at which the kinetic traces were taken. Note that the signals measured above 450 nm were multiplied by a factor of 3 to aid viewing.</figcaption></figure>

*Transient absorption spectra in the ns-to-µs range* were obtained to determine the fate of the excitation energy at longer timescales, with excitation wavelengths at 355, 514 and 570 nm. An overview of the datasets, along with the calculated EADS and selected kinetic traces, are shown in Fig. 6 for 355, 514 and 570 nm excitation. These EADS were the result of global analysis using an evolutionary model, and reflect a cascade of events with the following time constants: 20 ns → 2.3 µs → 128 µs. The evolution of intermediate species is the same for all the excitation wavelengths (355, 514 and 570 nm). The first EADS at 20 ns (black spectra) is convolved with the excitation laser, so it has been fixed to the value obtained from the fs-to-ns experiments. It shows two negative features peaking below 400 nm and in the 510–800 nm region, both of which can be assigned to GSB, whereas a featureless EADS appears in the 410–510 nm region. The observed spectrum is similar to the final EADS of the femtosecond dataset. The band composition of the second and third EADS are identical, showing a negative region around 400 nm that can be attributed GSB of the strong absorption peak at 355 nm. In the 440–580 nm region, there is ESA with valleys at 480 and 520 nm, which may correspond to the 0–2 and 0–1 vibronic peaks of the S<sub>0</sub> → S<sub>2</sub> electronic transitions, respectively. In monomeric lycopene, the T<sub>1</sub> → T<sub>n</sub> transition presents an ESA maximum circa 520 nm and decays over 12 µs<sup>50,51</sup>. The results for the lycopene aggregates suggest that there are at least two different triplet species formed. The species associated to the first EADS has been identified as a triplet in tomato chromoplasts<sup>22</sup>. The second triplet species (second, third EADS) resembles the triplets produced in monomeric lycopene by photosensitization<sup>51</sup>. The only significant difference observed here is the ground state bleaching band observed below 400 nm.

Discussion

The data from resonance Raman measurements shows that the v₁ frequency shifts with different excitation wavelengths, indicating that strong excitonic coupling between lycopene molecules affects the vibrational frequencies. From the analysis of the v₁ and v₂ regions, we observe that the sample contains at least four, and probably five species of lycopene aggregates. These observations are further supported by the complicated dynamics of femtosecond-to-nanosecond transient absorption. While it is, in principle, conceivable that a single type of aggregate or crystal could exhibit the number of evolutionary steps observed, it is hard to see how so many steps would result in the rise and fall of so many bands at more or less the same wavelengths. A more plausible explanation is that multiple types of aggregate are present in the sample with similarly (but not identically) positioned spectral bands, exhibiting overlapping dynamics.
<figure><img src="image_7.png" /><figcaption>**Fig. 5.** Dispersion-corrected time-gated pump-probe spectra of lycopene aggregates excited at 355 nm (A), 515 nm (B) and 570 nm (C). The times at which the spectra were taken are indicated on the graphs. Panels (D, E and F) show the EADS for the same datasets along with time constants, associated with different evolution steps, as *per* the legends.</figcaption></figure>

Although relaxation schemes involving singlet fission have been proposed to explain the formation of triplets in astaxanthin and lutein/violaxanthin aggregates¹²,²³, the current data seems too far from the established monomeric carotenoid dynamics to draw a straightforward connectivity scheme explaining all the experimental observations. The attempts at describing dynamics would first require a structural model that would reproduce the absorption spectrum in Fig. 1B. More specifically, the presence of aggregate absorption bands both red-shifted and blue-shifted from the main transition should be carefully addressed. Traditionally, such bands are attributed to J- and H-aggregates, resulting from two different types of transition dipole orientations. The quasi-one-dimensional structure of the π-conjugated chain of carotenoid molecules determines the direction of their S₀ → S₂ transition dipole along the molecular backbone. Hence, the two bands *must* result from (at least) two types of orientation (and therefore of coupling) of lycopene molecules in an aggregate. Both H- and J-type aggregates have been observed for astaxanthin ⁵², but they occur at different acetone-water ratios, allowing their separate investigation. Additionally, Fuciman et al. ²⁹ have reported that two spectrally-distinct H-aggregates of astaxanthin can be formed, and both produce long-lived triplet states, presumably via singlet fission. In line with these observations, the presence of a heterogeneous mixture of similar types of aggregates in the lycopene assemblies presented here would explain why the observed spectroscopic behavior is so complex.

Assuming our multi-aggregation hypothesis is correct, the observed spectroscopic and spectro-temporal features are a result of the following physical phenomena, taking place in multiple aggregate types:

1. Different exciton dynamics and intramolecular vibrational energy redistribution (IVR) within a single aggregate type (including energy equilibration within the excitonically-coupled S₂ manifold), and subsequent re-
<figure><img src="image_8.png" /><figcaption>**Fig. 6.** Nanosecond-to-microsecond transient absorption data of lycopene aggregates at 355, 515 and 570 nm excitation: (a, d, g) dataset overview; (b, e, h) EADS estimated by global fitting, using 3-component evolutionary model; (c, f, i) selected kinetics and fitting.</figcaption></figure>

laxation to S₁, S*, and other optically-dark states present in carotenoids. This would account for the redshifts of GSB bands, observed around 355 and 550–600 nm.

2. Vibrational cooling of the hot S₁ state, typically occurring in monomeric carotenoids on a sub-ps timescale ⁵³, which manifests itself as a narrowing and blue-shift of the IA band peaking at ca. 570 nm in lycopene monomers.

3. Excitation energy transfer between the different aggregate types. This would decay the GSB bands characteristic of one aggregate type and induce bleaching of the other bands, resulting in the delayed appearance of GSB signals.

4. Singlet-fission, occurring with different efficiencies and at different rates in different aggregates, eventually resulting in triplet states. Depending on the details of this mechanism, singlet fission can result in the growth of GSB signals, when two triplet excitations are produced from a singlet.

5. Equilibration, including triplet energy transfer, and inhomogeneous decay of produced triplet states via internal conversion, triplet-triplet annihilation, etc.

We have already hypothesized three coexisting moieties with at least five (S₂, hot S₁, S₁, S*, T₁) interconnected excited states. This is more than enough to over-parametrize the datasets that are satisfactorily described by the 8 time constants (five in the fs-ns range plus three in the ns-µs range). The electronic signatures for the excited states are highly distorted from the well-characterized signatures for monomeric or weakly-interacting carotenoids. It is reasonable that the excited states experiment a similar degree of aggregation-associated spectral distortion as that seen for the ground states in the absorption spectrum. However, based on the lifetimes obtained it is tempting to propose a tentative assignment as follows: < 500 fs EADS to S₂-like state, ca. 1p EADS to hot-S₁-like state, 10–11 ps EADS to S₁-like state, 46–60 ps EADS to S*-like state, which would then be the parent state for triplets, and the ca 550 ps EADS to interacting triplet states. Despite of not finding any evidence of CT-state we
cannot rule this possibility, since it can be masked by the plethora of processes evolving. To better understand the spectroscopic properties of lycopene aggregates, it is necessary to design an alternative experimental protocol where each aggregate type can be observed selectively. This would allow us to propose more sophisticated models than the current linear scheme. The selective purification of each species of lycopene aggregate would clearly be ideal, but this may well turn out to be extremely complex given the number of different aggregate types formed here and reported in the literature. An alternative would be to use time-resolved resonance Raman spectroscopy sub-selecting each of the populations, which should be sensitive enough to separate and distinguish the different excited states involved. Indeed, the selectivity of resonance Raman spectroscopy, and its sensitivity to the different arrangements of lycopene aggregates in their ground state, is illustrated very clearly in Fig. 2. Whatever the case, the presence of several aggregation-domains, within the large lycopene aggregates observed by TEM, remains the best hypothesis to explain the obtained spectroscopic data.

## Conclusions

We have prepared stable aggregates of lycopene in an acetone-water mixture and investigated their absorption properties, resonance spectra, and ultrafast transient absorption dynamics. Our results reveal highly complex dynamics, involving spectral bathochromic and hypsochromic shifts, and the decay and reappearance of multiple spectral bands. Eventually, the dynamics result in long-lived triplet states formed via singlet-fission. We postulate that the dynamics are the result non-homogeneous aggregates, as observed by TEM, which contain multiple, connected subdomains of H-aggregates, each with their own set of excited state dynamics. If this hypothesis is validated by further structure-sensitive experiments, multi-aggregation may prove to be an essential factor in understanding the properties of other carotenoid aggregates.

## Materials and methods

*Preparation of lycopene aggregates.* Tomato paste was used as raw material to obtain monomeric lycopene, extracted and purified as described elsewhere²². Briefly, tomato paste (25 g) was mixed with Milli-Q water (75 mL) and then dichloromethane (250 mL), then stirred at room temperature for 30 min. The mixture was filtered using a Büchner funnel, and the solid residue was rinsed three times with dichloromethane (50 mL each). After obtaining the filtrate, the aqueous phase was extracted with dichloromethane (3 × 75 mL). The combined organic phases were evaporated under reduced pressure. The resulting solid was purified by column chromatography on silica gel (hexane, followed by hexane/acetone 95/5 up to hexane/acetone 90/10 as gradient). The red solution was dried and then dissolved in a minimal amount of hexane and cooled at 7 °C for one hour, followed by storage at – 20 °C overnight to promote crystallization. Finally, the lycopene crystals were filtered using a Büchner funnel and rinsed. We improved the protocol to generate reproducible lycopene aggregates with low scattering, and high stability in suspension, as follows: distilled water (Millipore system 18.2 Ω cm) was added dropwise to a solution of 50 mg/L lycopene in absolute acetone (≥99%) under weak sonication at room temperature, to reach an acetone/water ratio of 80/20 (v/v). The aggregate solution (ca. 25 mg/L) was kept under sonication for 30 min, until the absorption remains stable.

*UV-Vis absorption spectra* were measured on a Varian Cary E5 Double-beam scanning spectrophotometer (Agilent), using a 4 mm path-length cuvette.

*Electron microscopy* The samples for electron microscopy were negatively stained with 2% uranyl acetate (w/v) on glow-discharged, carbon-coated copper grids. Data collection was performed using a Tecnai Spirit transmission electron microscope (TFS) equipped with a LaB6 filament, operating at 100 kV. Images were recorded on a K2 Base camera (Gatan/Ametek, 4kx4k) at 21,000 or 4400 magnification (pixel size at specimen level—1.9 and 0.83 nm, respectively).

*Resonance Raman spectra* were recorded at room temperature and 77 K, with laser excitations at 363.8, 488.0, 501.7, and 514.5 nm obtained from Coherent Ar+ (Sabre) and at 577.0 nm with a Genesis CX STM laser (Coherent). Output laser powers of 10–100 mW were attenuated to <5 mW at the sample and using unfocused laser beams. Scattered light was collected at 90° to the incident light, and focused into a Jobin-Yvon U1000 double-grating spectrometer (1800 grooves/mm) equipped with a red-sensitive, back-illuminated, LN2-cooled CCD camera. Sample stability and integrity were assessed based on the stability of the Raman signal.

*Nano-to-millisecond transient absorption* experiments were performed on a home-made nanosecond transient absorption setup, using a broadband OPO (optical parametric oscillator, ~4 ns temporal width) as pump laser, pumped by the 3rd harmonic of a Nd:YAG laser. A LEUKOS STM-2-UV super continuum laser (350–2400 nm, <1 ns temporal width) was used as probe laser. A system of lenses and mirrors shape the laser beams for maximum, homogeneous overlap on the sample. Analysis of the excited-state dynamics of self-assembled lycopene was performed using excitation at 570 nm, with an energy of 1.2 mJ per pump pulse. The frequency of the pump pulse was set at 10 Hz, while that of the probe pulse was reduced to 20 Hz using a rotating chopper. The sample was installed at ~45° relative to the pump and probe beams. After passing through the sample, the probe pulse was injected into a bundle of optical fibres (200 µm core), whose round entrance cross-section transitioned to a linear alignment at the exit in order to fit the entrance of the spectrometer. The latter was in turn coupled with an intensified charge coupled device (ICCD) camera. A reference spectrum of the probe was measured for each laser pulse, allowing correction for fluctuations. The time delay between pump and probe pulses was synchronized electronically through an in-house electronic pulse generator, to monitor temporally- and spectrally-resolved absorption spectra.

*Femto-to-nanosecond time-resolved transient absorption* spectra were recorded on a commercial transient absorption spectrometer (Harpia, Light Conversion). Pump and probe pulses for the spectrometer were derived from an amplified Ti:Sapphire laser (Libra, Coherent). Tunable pump pulses with center wavelength at 490 nm were obtained from an optical parametric amplifier (Topas-800, Light Conversion). A white light
continuum, used to probe the absorption changes in the 350–750 nm range, was generated in CaF₂ from the fundamental output of the laser. The time resolution of the instrument is ca. 120 fs. Excitation energies were set to approximately 600 nJ pulse.

Global analysis of time-resolved spectra was performed using commercial CarpetView data analysis software (Light Conversion).

## Data availability

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Received: 8 October 2024; Accepted: 24 January 2025

Published online: 15 February 2025

## References

1. Tavan, P. & Schulten, K. Electronic excitations in finite and infinite polyenes. *Phys. Rev. B* **36**, 4337-4358. https://doi.org/10.1103/PhysRevB.36.4337 (1987).

2. Llansola-Portoles, M. J., Pascal, A. A. & Robert, B. Electronic and vibrational properties of carotenoids: From in vitro to in vivo. *J. R. Soc. Interface* https://doi.org/10.1098/rsif.2017.0504 (2017).

3. Mendes-Pinto, M. M. et al. Electronic absorption and ground state structure of carotenoid molecules. *J. Phys. Chem. B* **117**, 11015-11021. https://doi.org/10.1021/jp309908r (2013).

4. Frank, H. A. & Christensen, R. L. in *Carotenoids: Volume 4: Natural Functions* (eds George Britton, Synnøve Liaaen-Jensen, & Hanspete Pfander) 167-188 (Birkhäuser Basel, 2008).

5. Frank, H. A., Young, A. J., Britton, G. & Cogdell, R. J. *Advances in Photosynthesis* (Kluwer Academic Publishing, 1999).

6. Hörvin Billsten, H., Zigmantas, D., Sundström, V. & Polivka, T. Dynamics of vibrational relaxation in the S1 state of carotenoids having 11 conjugated C=C bonds. *Chem. Phys. Lett.* **355**, 465-470. https://doi.org/10.1016/S0009-2614(02)00268-3 (2002).

7. Polivka, T. & Sundström, V. Dark excited states of carotenoids: Consensus and controversy. *Chem. Phys. Lett.* **477**, 1-11. https://doi.org/10.1016/j.cplett.2009.06.011 (2009).

8. Buckup, T. & Motzkus, M. Multidimensional time-resolved spectroscopy of vibrational coherence in biopolyenes. *Annu. Rev. Phys. Chem.* **65**, 39-57 (2014).

9. Balevičius, V., Abramavicius, D., Polivka, T., Galстian Pour, A. & Hauer, J. A unified picture of S\* in carotenoids. *J. Phys. Chem. Lett.* **7**, 3347-3352. https://doi.org/10.1021/acs.jpclett.6b01455 (2016).

10. Jailaubekov, A. E. et al. Deconstructing the excited-state dynamics of β-carotene in solution. *J. Phys. Chem. A* **115**, 3905-3916. https://doi.org/10.1021/jp1082906 (2011).

11. Kundu, A. & Dasgupta, J. Photogeneration of long-lived triplet states through singlet fission in lycopene H-aggregates. *J. Phys. Chem. Lett.* **12**, 1468-1474. https://doi.org/10.1021/acs.jpclett.0c03301 (2021).

12. Musser, A. J. et al. The nature of singlet exciton fission in carotenoid aggregates. *J. Am. Chem. Soc.* **137**, 5130-5139. https://doi.org/10.1021/jacs.5b01130 (2015).

13. Zuo, J. et al. Excited state properties of fucoxanthin aggregates. *Chem. Res. Chin. Univ.* **35**, 627-635. https://doi.org/10.1007/s40242-019-9097-2 (2019).

14. Chang, H.-T. et al. Singlet fission reaction of light-exposed β-carotene bound to bovine serum albumin. A novel mechanism in protection of light-exposed tissue by dietary carotenoids. *J. Agric. Food Chem.* **65**, 6058-6062. https://doi.org/10.1021/acs.jafc.7b01616 (2017).

15. Zhang, D. et al. Structure and excitation dynamics of β-carotene aggregates in Cetyltrimethylammonium bromide micelle. *Chem. Res. Chin. Univ.* **34**, 643-648. https://doi.org/10.1007/s40242-018-7379-8 (2018).

16. Wang, C. & Tauber, M. J. High-yield singlet fission in a zeaxanthin aggregate observed by picosecond resonance Raman spectroscopy. *J. Am. Chem. Soc.* **132**, 13988-13991. https://doi.org/10.1021/ja102851m (2010).

17. Wang, C., Angelella, M., Kuo, C.-H. & Tauber, M. J. Singlet fission in carotenoid aggregates: Insights from transient absorption spectroscopy. *Proceedings of SPIE* **8459**, 05-13. https://doi.org/10.1117/12.958612 (2012).

18. Billsten, H. H., Sundström, V. & Polivka, T. Self-assembled aggregates of the carotenoid zeaxanthin: Time-resolved study of excited states. *J. Phys. Chem. A* **109**, 1521-1529. https://doi.org/10.1021/jp044847j (2005).

19. Gradinaru, C. C. et al. An unusual pathway of excitation energy deactivation in carotenoids: Singlet-to-triplet conversion on an ultrafast timescale in a photosynthetic antenna. *Proc. Natl. Acad. Sci.* **98**, 2364-2369. https://doi.org/10.1073/pnas.051501298 (2001).

20. Papagiannakis, E. et al. Light harvesting by carotenoids incorporated into the b850 light-harvesting complex from *Rhodobacter sphaeroides* R-26.1: Excited-state relaxation, ultrafast triplet formation, and energy transfer to bacteriochlorophyll. *J. Phys. Chem. B* **107**, 5642-5649. https://doi.org/10.1021/jp027174i (2003).

21. Papagiannakis, E., Kennis, J. T. M., van Stokkum, I. H. M., Cogdell, R. J. & van Grondelle, R. An alternative carotenoid-to-bacteriochlorophyll energy transfer pathway in photosynthetic light harvesting. *Proc. Natl. Acad. Sci.* **99**, 6017-6022. https://doi.org/10.1073/pnas.092626599 (2002).

22. Llansola-Portoles, M. J. et al. Lycopene crystalloids exhibit singlet exciton fission in tomatoes. *Phys. Chem. Chem. Phys.* **20**, 8640-8646. https://doi.org/10.1039/C7CP08460A (2018).

23. Quaranta, A. et al. Singlet fission in naturally-organized carotenoid molecules. *Phys. Chem. Chem. Phys.* **23**, 4768-4776. https://doi.org/10.1039/D0CP04493H (2021).

24. Smith, M. B. & Michl, J. Singlet fission. *Chem. Rev.* (Washington) **110**, 6891-6936. https://doi.org/10.1021/cr1002613 (2010).

25. Pandya, R. et al. Optical projection and spatial separation of spin-entangled triplet pairs from the S1 (21 Ag-) State of Pi-conjugated systems. *Chem* **6**, 2826-2851. https://doi.org/10.1016/j.chempr.2020.09.011 (2020).

26. Casillas, R. et al. Molecular insights and concepts to engineer singlet fission energy conversion devices. *Energy Environ. Sci.* **13**, 2741-2804. https://doi.org/10.1039/D0EE00495B (2020).

27. Magne, C. et al. Perylene-derivative singlet exciton fission in water solution. *Chem. Sci.* **15**, 17831-17842. https://doi.org/10.1039/D4SC04732J (2024).

28. Wang, C., Schlamadinger, D. E., Desai, V. & Tauber, M. J. Triplet excitons of carotenoids formed by singlet fission in a membrane. *ChemPhysChem* **12**, 2891-2894. https://doi.org/10.1002/cphc.201100571 (2011).

29. Fuciman, M., Durchan, M., Šlouf, V., Kešan, G. & Polivka, T. Excited-state dynamics of astaxanthin aggregates. *Chem. Phys. Lett.* **568-569**, 21-25. https://doi.org/10.1016/j.cplett.2013.03.009 (2013).

30. Barford, W. Singlet fission in lycopene H-aggregates. *J. Phys. Chem. Lett.* https://doi.org/10.1021/acs.jpclett.3c02435 (2023).

31. Peng, B., Wang, Z., Jiang, J., Huang, Y. & Liu, W. Investigation of ultrafast intermediate states during singlet fission in lycopene H-aggregate using femtosecond stimulated Raman spectroscopy. *J. Chem. Phys.* https://doi.org/10.1063/50200802 (2024).

32. Ruban, A. V. et al. Identification of a mechanism of photoprotective energy dissipation in higher plants. *Nature* **450**, 575-578. https://doi.org/10.1038/nature06262 (2007).
33. Pascal, A. A. et al. Molecular basis of photoprotection and control of photosynthetic light-harvesting. *Nature* **436**, 134–137 (2005).

34. Kasha, M. Energy transfer mechanisms and the molecular exciton model for molecular aggregates. *Radiat. Res.* **20**, 55–70. https://doi.org/10.2307/3571331 (1963).

35. Hestand, N. J. & Spano, F. C. Expanded theory of H- and J-molecular aggregates: The effects of vibronic coupling and intermolecular charge transfer. *Chem. Rev.* (Washington, DC, U. S.) **118**, 7069–7163. https://doi.org/10.1021/acs.chemrev.7b00581 (2018).

36. Spano, F. C. Analysis of the UV/Vis and CD spectral line shapes of carotenoid assemblies: Spectral signatures of chiral H-aggregates. *J. Am. Chem. Soc.* **131**, 4267–4278. https://doi.org/10.1021/ja806853v (2009).

37. Lu, L., Liu, Y., Wei, L., Wu, F. & Xu, Z. Structures and exciton dynamics of aggregated lutein and zeaxanthin in aqueous media. *J. Lumin.* **222**, 117099. https://doi.org/10.1016/j.jlumin.2020.117099 (2020).

38. Hempel, J., Schädle, C. N., Leptihn, S., Carle, R. & Schweiggert, R. M. Structure related aggregation behavior of carotenoids and carotenoid esters. *J. Photochem. Photobiol. A Chem.* **317**, 161–174. https://doi.org/10.1016/j.jphotochem.2015.10.024 (2016).

39. Ishigaki, M. et al. Unvelling the aggregation of lycopene in vitro and in vivo: UV–Vis, resonance Raman, and Raman imaging studies. *J. Phys. Chem. B* **121**, 8046–8057. https://doi.org/10.1021/acs.jpcb.7b04814 (2017).

40. Wang, L., Du, Z., Li, R. & Wu, D. Supramolecular aggregates of lycopene. *Dyes Pigm.* **65**, 15–19. https://doi.org/10.1016/j.dyepig.2004.05.012 (2005).

41. Salares, V. R., Young, N. M., Carey, P. R. & Bernstein, H. J. Excited state (excitation) interactions in polyene aggregates. Resonance Raman and absorption spectroscopic evidence. *J. Raman Spectrosc.* **6**, 282–288. https://doi.org/10.1002/jrs.1250060605 (1977).

42. Wang, C., Berg, C. J., Hsu, C.-C., Merrill, B. A. & Tauber, M. J. Characterization of carotenoid aggregates by steady-state optical spectroscopy. *J. Phys. Chem. B* **116**, 10617–10630. https://doi.org/10.1021/jp3069514 (2012).

43. Streckaite, S. et al. Modeling dynamic conformations of organic molecules: Alkyne carotenoids in solution. *J. Phys. Chem. A* **124**, 2792–2801. https://doi.org/10.1021/acs.jSEA.9b11536 (2020).

44. Llansola-Portoles, M. J., Pascal, A. A. & Robert, B. in *Methods Enzymol*. Vol. 674 (ed Eleanore T. Wurtzel) 113–135 (Academic Press, Cambridge, 2022).

45. Koyama, Y., Takatsuka, I., Nakata, M. & Tasumi, M. Raman and infrared spectra of the all-trans, 7-cis, 9-cis, 13-cis and 15-cis isomers of β-carotene: Key bands distinguishing stretched or terminal-bent configurations form central-bent configurations. *J. Raman Spectrosc.* **19**, 37–49. https://doi.org/10.1002/jrs.1250190107 (1988).

46. Koyama, Y., Takii, T., Saiki, K. & Tsukida, K. Configuration of the carotenoid in the reaction centers of photosynthetic bacteria. 2. Comparison of the resonance Raman lines of the reaction centers with those of the 14 different cis-trans isomers of β-carotene. *Photobiochem. Photobiophys.* **5**, 139–150 (1983).

47. Koyama, Y. et al. Configuration of the carotenoid in the reaction centers of photosynthetic bacteria. Comparison of the resonance Raman spectrum of the reaction center of *Rhodopseudomonas sphaeroides* G1C with those of cis-trans isomers of β-carotene. *Biochim Biophys. Acta, Bioenerg.* **680**, 109–118. https://doi.org/10.1016/0005-2728(82)90001-9 (1982).

48. Rimai, L., Heyde, M. E. & Gill, D. Vibrational spectra of some carotenoids and related linear polyenes. Raman spectroscopic study. *J. Am. Chem. Soc.* **95**, 4493–4501. https://doi.org/10.1021/ja00795a005 (1973).

49. Lutz, M., Szponarski, W., Berger, G., Robert, B. & Neumann, J.-M. The stereoisomerization of bacterial, reaction-center-bound carotenoids revisited: An electronic absorption, resonance Raman and NMR study. *Biochem. Biophys. Acta* **894**, 423–433. https://doi.org/10.1016/0005-2728(87)90121-6 (1987).

50. Burke, M., Land, E. J., McGarvey, D. J. & Truscott, T. G. Carotenoid triplet state lifetimes. *J. Photochem. Photobiol. B* **59**, 132–138. https://doi.org/10.1016/S1011-1344(00)00150-0 (2000).

51. Land, E. J., Sykey, A. & George Truscott, T. The in vitro photochemistry of biological molecules II. The triplet states of β-carotene and lycopene excited by pulse radiolysis. *Photochem. Photobiol.* **13**, 311–320. https://doi.org/10.1111/j.1751-1097.1971.tb06119.x (1971).

52. Köpsel, C. et al. Structure investigations on assembled astaxanthin molecules. *J. Mol. Struct.* **750**, 109–115. https://doi.org/10.1016/j.molstruc.2005.02.038 (2005).

53. de Weerd, F. L., van Stokkum, I. H. M. & van Grondelle, R. Subpicosecond dynamics in the excited state absorption of all-trans-β-Carotene. *Chem. Phys. Lett.* **354**, 38–43. https://doi.org/10.1016/S0009-2614(02)00095-7 (2002).

## Acknowledgements

This work benefited from the Resonance Raman and CryoEM platforms of I2BC, supported by the French Infrastructure for Integrated Structural Biology (FRISBI) ANR-10-INSB-05. The work was supported by the Agence Nationale de la Recherche (ANR) (EXCIT Grant No: ANR-20-CE11-0022, SINGLEFISSION Grant No: ANR-23-CE29-0007 and FISCIENCY Grant No: ANR-23-CE50-009) and the CNRS as part of its 80PRIME interdisciplinary program. M.V. acknowledges “the Universities Excellence Initiative” programme by the Ministry of Education, Science and Sports of the Republic of Lithuania under the agreement with the Research Council of Lithuania (project No. S-A-UEI-23-6).

## Author contribution

Chloe Magne and Vasyl Veremeienko performed the extraction of pigments, synthesis of aggregates, Raman measurements, TEM measurements, and ns-to-microsecond measurements. Roxanne Bercy contributed to the extraction of pigments and synthesis of aggregates. Minh-Huong Ha-Thi conducted ns-to-microsecond measurements, while Ana A. Arteni carried out TEM measurements. Mikas Vengris performed femtosecond TA measurements and contributed to writing the paper. Andrew A. Pascal, Thomas Pino, Bruno Robert, and Manuel J. Llansola-Portoles conceived the idea, wrote the manuscript, and analyzed the data. All authors contributed to the writing and revision of the final version of the manuscript.

## Competing interests

The authors have no competing (include financial AND non-financial) interests as defined by Nature Research, or other interests that might be perceived to influence the results and/or discussion reported in this paper.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-025-88220-z.

Correspondence and requests for materials should be addressed to M.J.L.-P.

Reprints and permissions information is available at www.nature.com/reprints.
**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access** This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. 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://creativecommons.org/licenses/by-nc-nd/4.0/.

© The Author(s) 2025