RPP design strategy and characterization

MPA characterizes a perturbative nonlinear polarization response to excitation intensity, where a n-photon absorption process is governed by the imaginary part of the (2n-1)-th order nonlinear susceptibility \({\chi }^{\left(2n-1\right)}(n\,=2,\,3,\,4\ldots )\). Accordingly, the third-order (\({\chi }^{(3)}\)), fifth-order (\({\chi }^{(5)}\)), and seventh-order (\({\chi }^{(7)}\)) nonlinear susceptibilities are responsible for 2PA, 3PA, and 4PA processes, respectively, and so forth34. Thus, to realize high-efficiency MPA nonlinearities under moderate laser intensities to avoid phototoxicity and material damage, a crucial approach is to enhance the imaginary part of high-order \({\chi }^{(n)}\). In the context of MPA processes, an approximation of \({\chi }^{(n)}\) can be simplified as \({\chi }^{\left(n\right)} \sim {\mu }^{n+1}/{{(\hslash \omega }_{0})}^{n} \sim {\left[{P}^{(1)}\right]}^{n+1}\), where \(\mu\) is the dipole moment, \({\omega }_{0}\) is the applied field frequency, and \({P}^{(1)}\) is the linear polarization34. Therefore, a larger dipole moment inherently leads to increased \({P}^{(1)}\), thereby enhancing multiphoton nonlinearities.

Further considering the contribution of excitonic dipoles to optical nonlinearity17,18,19,20,21,22,23,24,25,26,27,28,29,30, designing RPPs capable of high-order multiphoton upconversion should prioritize the following aspects: (1) Utilizing strongly electronegative halide anions to induce a large bandgap in VIS region and facilitate high-order MPA. (2) Employing thin quantum well (i.e., few-layered inorganic sheets) with strong quantum confinement to enhance exciton binding energy and prevent exciton dissociation under intense excitation. (3) Maximizing the contrast of dielectric constants between the inorganic sheets and organic cations to reduce charge screening, thereby enhancing Coulomb attraction and enlarging exciton binding energy35. (4) Establishing strong interlayer hydrogen bonds and well-matched lattices at the abrupt organic-inorganic interface for ensuring high environment stability and robust resistance to damage under intense excitation.

Based on all the concerns mentioned above, we have designed three layered RPP structures: (CA)2(MA)Pb2I7 (CMPI, orthorhombic), (i-BA)2(MA)2Pb3Br10 (i-BMPB, orthorhombic), and (F-BEA)2PbBr4 (F-BPB, monoclinic), as illustrated in Fig. 1a (please see Supplementary Note 1 for detailed crystal structure characterizations). These three crystal structures feature inorganic frameworks composed of corner-sharing PbX6 (X = Br, I) octahedra that link two layers of hydrophobic organic spacer cations (CA, i-BA or F-BEA) along the a-axis through strong N-H ∙ ∙ ∙ X hydrogen bonds. All the structures exhibit inverse symmetry in bc plane, with about 45o between the Pb-X bonding and in-plane b (or c)-axis directions (the left panels of Fig. S2a–c). Such an interlayer arrangement constitutes wells and barriers, forming a two-dimensional (2D) superlattice motif that enables a platform for tailoring individual bandgaps within the RGB gamut range (the middle panels of Fig. S2a-c). Further, theoretical analyses on the partial density of states and band alignments certify the crucial role of inorganic layers in establishing the natural Type-I QWs of these three RPPs (the right panels of Fig. S2a-c). Moreover, owing to the high dielectric contrast between the well and barrier layers, the confinement potential energies of these three structures are significantly higher than those of many reported RPPs with bandgaps in RGB gamut range (Figs. S2 and S3; detailed data are listed in Table S2), which is advantageous for stabilizing excitonic behaviors.

Fig. 1: Structural properties and the principle of MPA upconversion for layered CMPI, i-BMPB, and F-BPB.

Fig. 1: Structural properties and the principle of MPA upconversion for layered CMPI, i-BMPB, and F-BPB.

a Diagram of the MQW nanostructures for top view in ab plane with alternative arrangements of inorganic perovskite sheets (well) and organic cation layers (barrier). b Linear optical absorption spectra with strong excitonic absorption peaks. c Comparison of the exciton binding energies (309 meV for CMPI, 246 meV for i-BMPB, and 332 meV for F-BPB) with typical perovskites with RGB gamut bandgaps. References: MAPbI336, MAPbBr337, MAPbI3-xClx38, FAPbBr339, CsPbI340, CsPbBr341, CsPbI2Br40, CsPbCl342, (BA)2(MA)Pb3I1043, (BIA)2PbI444, (BZA)2PbI444, (R-MBA)2PbI445, (p-FPEA)2PbBr446, (BA)2SnI447, (TEA)2SnI448, and (R-MBA)2Pb0.8Sn0.2I445. d Top panel: microscopic mechanism for the enhancement of MPA nonlinearities enabled by dynamic BOED in the case of 6PA process. Bottom panel: comparison of the physical models between atomic/molecular dipoles and excitonic dipoles

After synthesizing bulk single crystals of the three RPPs (see Materials and methods for details), their linear optical absorption spectra were measured as shown in Fig. 1b. Figure 1b manifests distinct absorption peaks at 556 (CMPI), 462 (i-BMPB) and 397 nm (F-BPB), corresponding to the lowest energy of excitonic transition35. Their linear photoluminescence (PL) spectra are centered at 640, 530, and 442 nm, respectively (Fig. S4), well located at the RGB gamut. From these measurements, the optical bandgap energies are estimated to be approximately 1.8 (CMPI), 2.2 (i-BMPB), and 2.8 eV (F-BPB). As shown in Fig. 1c, the measured exciton binding energies Eb are 309 meV for CMPI, 246 meV for i-BMPB, and 332 meV for F-BPB (see Supplementary Note 3 for detailed measurement). These values are substantially larger than those reported for prominent 2D/3D perovskites with RGB gamut bandgaps36,37,38,39,40,41,42,43,44,45,46,47,48 (detailed data are listed in Table S3), indicating their great potential in excitonic multiphoton nonlinearity.

Origin of giant multiphoton nonlinearities

Before delving into the experimental investigation on the MPA properties of these RPPs, it is essential to discuss the physical origins of the possible strong multiphoton nonlinearities in a robust excitonic system. According to classic theories regarding resonant optical nonlinearities of excitonic dipoles in systems with strong dielectric or quantum confinement17,18,19, we anticipate that the microscopic MPA response of the as-grown RPPs is linked to the coherent dynamics of excitonic dipoles when subjected to femtosecond (fs) laser driving.

To elucidate the physical origin of giant multiphoton nonlinearities, the high-order MPA (6PA) induced interband transitions and upconversion PL process are taken as an example. As schematically shown in Fig. 1d, in the initial state of laser pulse excitation, some valence electrons in the well are excited from the valence to conduction band through virtual intermediate levels via atomic/molecular dipole polarization in crystals. Simultaneously, exciton states are formed with tight Coulomb attraction between the excited electrons and created valence holes. This band-to-band correlation can induce significant polarization of the excitonic dipoles. Laser pulses then drive these excitons moving within valence and conduction bands.

In relatively weak excitonic systems, such as 3D perovskites with low exciton binding energies, strong excitation can drive excitons to higher energy levels, causing them to dissociate into free electrons and holes. In the absence of significant scattering, these electrons and holes can be accelerated to move periodically within their respective bands and eventually diffract at the Brillouin zone boundaries. This delocalized dynamic is periodical at Bloch frequency and known as intraband Bloch oscillations of electrons (holes)49,50,51, producing nonlinear intraband currents that are widely recognized as the dominate source of high-harmonic generation in solids49,50,51,52,53.

However, it is important to note that the designed RPPs here exhibit strong excitonic characteristics with large exciton binding energies, making excitons dissociation difficult even under strong laser excitation. In this regime, the electron and hole remain bound and move together as a unit following the laser-field driving. Consequently, the Bloch oscillations become localized in the vicinity of conduction band minima (CBM) and valence band maxima (VBM), governed by the combined effect of the laser-field driving and excitonic Coulomb binding. Here, we highlight this local motion of excitons as BOED, in which the coherent superposition of excitonic dipoles results in a new source of collective linear polarization moment \({P}^{(1)}\). The component of the electric dipole moment parallel to the electric field vector (i.e., the polarization direction) of pump laser can contribute an enhancement of each-order nonlinear susceptibility for MPA, amplifying electron transition, producing a high density of excitonic dipoles, and eventually leading to giant multiphoton nonlinearities. Exciton recombination occurs through decoherence from the Bloch states near the CBM and VBM, followed by relaxation of electrons from the conduction band to the excitonic energy level and subsequently back to the valence band. This process is accompanied by the emission of a cascade of upconversion photons.

It is noteworthy that in conventional nonlinear optics for a crystal system, as depicted in the bottom figure of Fig. 1d, the generated atomic/molecular dipoles exhibit similar polarization strength due to the consistent optical and crystal field regulations. However, in the BOED scenario, excitonic dipoles experience different polarization strengths because they occupy various energy states during the oscillations along the bands. This distinction suggests that nonlinearities governed by excitonic dipole moments should be more pronounced at higher-order MPA, overcoming the limitations of conventional nonlinearities. This scenario is further supported by the stronger enhancement effect of upconverted PL at higher-order MPA, as discussed in the subsequent section.

Measurement of giant multiphoton nonlinearities in RPP bulks

The experimental setup for characterizing multiphoton nonlinearities of the three as-grown bulk RPPs is presented in Fig. 2a. Stable MPA upconversion PL was observed with broadband tuning of pump wavelength (λp) from NIR-I to NIR-IV windows (800–2300 nm). The corresponding PL spectra and photographs are shown in Fig. 2b(I)-d(I) (λp = 2300 nm) and Fig. S6 (other pump wavelengths). The peaks of the PL spectra for CMPI, i-BMPB, and F-BPB are centered at 625, 525, and 440 nm, respectively, well within the RGB gamut. Further, the dependence of PL intensity on pump intensity (Ip) was studied under NIR-II/III/IV excitation (Fig. 2b(II)–d(II)). The power-law fitting, denoted by exponential function \({I}_{{\rm{p}}}^{\eta }\) (\(\eta \approx n\) for \(n\)-photon absorption), proves that the PL-pump dependence of CMPI, i-BMPB, and F-BPB fulfills the perturbative nonlinear-optics characteristics of 2PA-4PA, 3PA-5PA, and 3PA-6PA processes, respectively54,55,56,57. Notably, the minimum pump intensities (\({I}_{\min }\)) for 3PA-6PA are significantly lower than those of reported materials55,56,57,58,59,60,61,62,63,64 and well below the bio-damage threshold65,66,67, as demonstrated in Fig. 2e (detailed data are listed in Table S4), highlighting their potential for biophotonics applications.

Fig. 2: MPA-induced RGB upconversion emission from the bulks of CMPI, i-BMPB and F-BPB.

Fig. 2: MPA-induced RGB upconversion emission from the bulks of CMPI, i-BMPB and F-BPB.

a Schematic diagram of the experimental setup. b–d (I) PL spectra under different pump intensities (λp = 2300 nm). Inset: photographs of upconversion PL from individual RPP bulks under NIR-IV pumping (λp = 2300 nm). Orange arrows indicate the wave vector (k) and the crystallographic b- and c-axes. (II) Logarithmic plot of the MPA PL peak strength dependence on pump intensity at different wavelengths within NIR-II/III/IV windows. The solid lines represent fitting with the power-law equation \({I}_{{\rm{p}}}^{\eta }\), where the extracted values of \(\eta\) follow an \(\eta \approx n\) dependence for \(n\)-photon absorption. This verifies a perturbative scaling of nonlinear optics for MPA process. e Comparison of the minimum pump intensity in 3PA-6PA processes with representative results reported from other materials. The pump intensity \({I}_{\min }\) refers to the required minimum pump intensity to excite the MPA-induced upconversion PL. References: MAPbBr355, CsPbBr355, MAPbBr3/(OA)2PbBr455, Zn2(SDC)2(An2Py)56, Zn2(SDC)2(An2Py)•anthracene56, CsPbCl357, Cs3Cu2I558, CsCu2I358, ZnO59, ZJU-28\(\supset\)MAPbBr360, CsPbBr3@ZIF-861, FAPbBr3/CsPbBr362, SpL-363, and MAPbCl364. The red horizontal bar indicates the bio-damage threshold referring to ref. 65,66,67. f PL intensity variation with time under long-time fs laser irradiation (λp = 1200 nm and Ip = 7.2 GW/cm2). The duration is 120 min with an interval of 3.5 min

To accurately investigate the high-order nonlinearities related to MPA, open aperture (OA) Z-scan experiments were performed to determine the absorption coefficients and absorption cross-sections for MPA processes in the three RPP bulks (see details in Supplementary Note 5). As anticipated, the results for 3PA-6PA (Table 1) surpass those of almost all the reported MPA materials to date (detailed comparison is provided in Table S5). The absorption cross-sections for 5PA (9.1 × 10−134 cm10s4photon−4 for F-BPB) and 6PA (4.2 × 10−163 cm12s5photon−5 for F-BPB) are approximately three and six orders of magnitude larger than those for reported MPA materials, respectively55,59,63.

Table 1 Multiphoton absorption cross-section σn in NIR-II/III/IV windows

Furthermore, we monitored the MPA upconversion PL intensities of the three RPP bulks under long-time (120 min) intense pumping (λp = 1200 nm, Ip = 7.2 GW/cm2) to evaluate the MPA stability (see details in Fig. S7). As depicted in Fig. 2f, the PL intensities remain ~75% (CMPI), ~71% (i-BMPB), and ~72% (F-BPB) of the initial values, indicating superior laser-damage resistance. Additionally, these RPPs exhibit excellent long-term environmental stability, retaining over 90% of their initial PL intensity after storage in a vacuum box for 15 months (Fig. S8). These giant multiphoton nonlinearities are attributed to the BOED dynamics, supported by the experimental observations of anomalous MPA evolutions in anisotropy strength as discussed in the following.

Anomalous anisotropy of MPA upconversion PL

The exceptional single-crystalline nature of the three RPP bulks allows us to further investigate the PL polarization dependence. The anisotropic evolutions of MPA upconversion PL under different crystallographic orientation in bc plane were measured by employing a half-wave plate to rotate the angle θ between the c-axis direction of the samples and the linear-polarized direction of the pump laser (Fig. 2a). The results reveal significant difference between one-photon absorption (1PA) and MPA processes. The 1PA PL remains isotropic regardless of pump strength (see details in Fig. S14), but the MPA upconversion PL exhibits exotic anisotropic behavior that depends on the pump strength.

To illustrate this anisotropic evolution in MPA upconversion PL, we used i-BMPB as an example and employed both a Ti:sapphire oscillator and an amplifier to realize a wide tuning range of pump intensity at λp = 800 nm. The 2PA PL strength exhibits an anomalous evolution from isotropy (at Ip = 35.8 kW/cm2) to anisotropy with increasing pump strength. This anisotropy displays distinct in-plane four-fold symmetry with maxima at θ = 45°, 135°, 225°, and 315° (see the corresponding polar plots in Fig. 3a). Further increasing the pump strength leads to a gradual enhancement of anisotropy strength, with the extinction ratio ε (defined as the ratio of PL intensity between θ = 45o and 0o) increasing from 1 to ~1.5 (Fig. 3a, c).

Fig. 3: Anomalous evolutions of MPA anisotropy in the bulk RPPs.

Fig. 3: Anomalous evolutions of MPA anisotropy in the bulk RPPs.

a Representative polar plots demonstrate the orientation evolutions of 2PA PL (λp = 800 nm) in i-BMPB for large-scale variation of pump intensity. The minima at 0o and 180o along c-axis, 90o and 270o along b-axis; while the maxima at 45o, 135o, 225o, and 315o corresponds to the PbX6 octahedra corner-connecting directions. b Orientation evolutions of 5PA PL (λp = 2300 nm) in i-BMPB for available scale of pump intensity. c Retrieved extinction ratios as a function of pump intensity in different MPA processes for i-BMPB. d Orientation evolutions of MPA PL in i-BMPB at different pump wavelengths under a constant pump intensity of 5.1 GW/cm2. The green hollow dots in (a), (b), and (d) are experimental data, while the red lines depict the simulation results based on BOED. e Extinction ratios as a function of pump wavelength under a constant pump intensity of 5.1 GW/cm2 for the three RPPs. The solid lines are linear fittings for the experimental data in each order of MPA processes

For higher-order MPA, similar behavior of growing anisotropy was recorded. Polar plots of 5PA (Fig. 3b) and 4PA (Fig. S16c) with corresponding ε values presented in Fig. 3c demonstrate this behavior. Furthermore, ε increases monotonically with the pump wavelength, even within the same MPA order (i.e., the same order of optical nonlinearity), as shown in Fig. 3d and the middle panel in Fig. 3e. The polarization-dependent absorbance exhibits the same anomalous anisotropic behavior as the PL (the corresponding polar plots are given in Fig. S16a–c). This indicates that the anisotropic PL evolutions result directly from the anisotropy of MPA process, without significant post-absorption energy-loss pathways such as defect-state-induced luminescence. Likewise, very similar anisotropy evolutions were experimentally observed in CMPI and F-BPB (the top and bottom panels of Fig. 3e and Fig. S19, see Figs. S17 and S18 for detailed polar plots).

These exotic anisotropy evolutions of MPA upconversion PL cannot be simply interpreted by conventional nonlinear optics based on atomic/molecular dipole moments34, as discussed below. (1) The symmetry of atomic/molecular dipole polarization depends on the inversion symmetry of crystal structure, which means that the integrated nonlinear susceptibility in each crystallographic direction remains constant under laser excitation. In this regard, the PL symmetry after MPA should not change with pump intensity at a fixed pump wavelength, thus should not induce transition from isotropy to anisotropy. (2) Conventional laws in nonlinear optics also suggest that longer pump wavelengths (i.e., lower photon energies) result in weaker polarization of atomic/molecular dipoles, and thus a lower extinction ratio. This is contrary to our experimental observations, where an increased extinction ratio is observed with longer wavelengths, indicating a need for a new understanding of high-order nonlinear effects in our experiments.

To address the anisotropy of MPA upconversion PL within these RPPs, we explore the role of BOED in determining the symmetry of MPA. As previously discussed, BOED is localized within the energy bands, making the anisotropy of excitonic nonlinearities dependent upon the band structure. First-principles calculations (Fig. S20) reveal that the energy bands of the RPPs exhibit a clear four-fold symmetry from CBM and VBM (Γ point) to the Brillouin zone boundary with anisotropy dispersion. Also, we simulated the interaction between the laser field and excitonic dynamics in RPPs by solving two-band semiconductor Bloch equations (SBEs) integrated with excitonic dipoles and anisotropy dispersion of energy bands (see Materials and methods for details). Most of the excitons generated through MPA processes initially localize around Γ point because the energy difference between CBM and VBM is minimum and the corresponding interband transition probability is maximum. The excitons can acquire momentum from the laser field and then oscillate in k space, with the additional momentum determined by vector potential \(A\left(t\right)=-{\int }_{0}^{t}E({t}^{{\prime} }){{dt}}^{{{\prime} }}\) of laser field.

Deduced from Eqs. 1 and 2 in Materials and methods, the maximum vector potential \({A}_{\max }\propto {E}_{0}{\lambda }_{{\rm{p}}}\), where \({E}_{0}\) is the amplitude of the laser electric field. For short pump wavelength (i.e., low-order MPA regime) and low pump intensity, the \(A(t)\) value is small, leading to motion and localization of the excitonic dipoles very close to Γ point. Given that the energy bands of the three RPPs are approximately isotropic within the region very close to Γ point (Fig. S20), the MPA dynamics behave as isotropic or weak anisotropic. With increasing the laser intensity or pump wavelength (i.e., higher-order MPA regime), \(A(t)\) becomes larger, driving the excitons into anisotropic k space to facilitate BOED tracing along the bands at θ = 45°, 135°, 225°, and 315°, thereby generating PL dynamics from MPA upconversion with strong anisotropy. Theoretical simulations of BOED along these anisotropic energy bands align with the experimentally observed anomalous anisotropic evolutions of MPA upconversion PL (results for i-BMPB are presented in Fig. 3a, b, d and Fig. S16c, results for CMPI and F-BPB are shown in Figs. S17 and S18).

Full-color MPA upconversion in nanoscale RPPs

To take full advantage of the giant multiphoton nonlinearities enabling RGB upconversion, we explore full-color upconversion using nanostructured RPPs to meet the demands of microimaging applications. High-yield preparation of low-dimensional nanocrystals (NCs) was conveniently realized through ultrasonic exfoliation of the three RPP bulks (see Materials and methods for details). Transmission electron microscope (TEM) images reveal that the mean sizes of the as-prepared NCs of CMPI, i-BMPB, and F-BPB are 45, 13, and 70 nm, respectively (upper figures of Fig. 4a, with corresponding selected area electron diffraction (SAED) patterns presented in the bottom figures of Fig. 4a). It is important to emphasize that obtaining clear high-resolution TEM images for characterizing crystal structure in more micro scale is challenging due to the high resistance of the three RPPs induced by robust excitonic effects (approximately 109 Ω·cm, measured by the space charge limited current method). Despite this, the intrinsic single-crystal structures of the RPPs are preserved in the NCs since they are merely physically exfoliated from the single-crystal bulks without any chemical re-treatment.

Fig. 4: Full-color MPA upconversion PL emission from NCs of CMPI, i-BMPB, F-BPB and their composites.

Fig. 4: Full-color MPA upconversion PL emission from NCs of CMPI, i-BMPB, F-BPB and their composites.

a Top panels: TEM images of CMPI, i-BMPB, and F-BPB NCs. Bottom panels: corresponding size distributions and SAED patterns. b Full-color MPA PL spectra and photographs including eight main VIS colors and white color pumped by NIR-III/IV lasers. The marked ratios denote the mixing molar ratios of the three NCs in each composite. The inset numbers in brackets refer to the CIE1931 coordinates. c Schematics for four types of binary/ternary NC composites. d, e SEM image and EDS elemental mapping of C/B composite. f CIE coordinates (circle) as marked in (b) plotted on a chromaticity diagram. The full-color MPA upconversion encompasses a broader range of color gamut (solid line) compared to that of sRGB (dotted line)

As anticipated, the three NC samples also demonstrated efficient high-order MPA upconversion under NIR-III/IV laser excitation (see Fig. 4b for the MPA PL spectra and photographs). The PL peaks exhibit a slightly blue shift within 5 nm with respect to their bulk counterparts due to the enhanced quantum confinement of low-dimensional nanostructures, yet they still lie within the standard RGB (sRGB) gamut. The stronger quantum confinement of the NCs also results in larger MPA cross-sections compared to the RPP bulks (Table 1). The 5PA and 6PA absorption cross-sections of F-BPB NCs are five and eight orders of magnitude larger, respectively, than those reported for other MPA materials55,59,63 (see detailed comparison in Fig. S13 and Table S5).

By mixing the RGB primary colors with appropriate ratios of the three NC samples, a wide range of additive colors can be achieved. As presented in Fig. 4b, we prepared several mixing suspensions based on the RGB intensity ratios in the CIE color model. These suspensions include CMPI/i-BMPB (C/B), CMPI/F-BPB (C/F), i-BMPB/F-BPB (B/F), and CMPI/F-BPB/i-BMPB (C/F/B), enabling MPA upconversion PL in five main colors (orange, yellow, cyan, magenta, and violet) and in white color under NIR-III/IV excitation. Further, as schematically illustrated in Fig. 4c, a series of binary and ternary composites were fabricated by depositing different mixing suspensions on a substrate, followed by rapid drying to form dense NC films.

Here we choose the C/B sample as a representative to investigate the properties of the composites. As depicted in Fig. 4d, e, the C/B composites exhibit well-formed structures as observed in scanning electron microscope (SEM) image and energy-dispersive X-ray spectroscopy (EDS) elemental mapping. It should be noted that, similar to the TEM measurement, the SEM resolution is also hindered by the high resistance of these RPPs. Thus, the SEM images in Fig. 4d, e were recorded by employing larger-size NC samples (prepared under shorter ultrasonic exfoliation time) than those in Fig. 4a. But it can be inferred that successful formation of nanocomposites with smaller NCs has been achieved considering the identical fabrication processes. The SEM and EDS images of C/F and B/F also reveal their composite morphology, indicating successful implementation of composite fabrication (Fig. S21). Furthermore, time-resolved photoluminescence (TRPL) measurements (Figs. S22 and S23) provide clear evidence of strong interactions within the composites. The large-bandgap NCs exhibit a shortened PL lifetime, while the small-bandgap NCs show a prolonged lifetime, indicating efficient resonance energy transfer. These results not only confirm the successful formation of nanoscale composites but also help to explain slight spectral blue-shift observed in Fig. 4b.

The MPA upconversion PL colors of these NC composites are in line with those of mixing suspensions, according to the calculations of CIE coordinates marked in Fig. 4b. Figure 4f depicts that the MPA upconversion covers an RGB triangle’s area with 7.7% more perceptible colors than sRGB region (the calculation method is given in Supplementary Note 8). Moreover, owing to the high structural stability of the three RPPs, the composites exhibit good PL stability and long-term damage resistance even under intense excitation (Ip = 9.2 GW/cm2) for 90 min (we take C/F as an example to demonstrate the stability measurement in Fig. S24). These exceptional performances suggest that these NC composites are promising candidates for full-color upconversion displays of NIR-III/IV at the nanoscale. Overall, these results provide strong evidence for the presence of BOED and determine its crucial role in enhancing MPA nonlinearities. We speculate that BOED could be a universal dynamic in robust excitonic systems, which may pave the way for developing material physics towards robust high-order optical nonlinearities.

Applications of full-color MPA upconversion nonlinearities

Inspired by the robust full-color MPA upconversion, we explored its potential applications in upconversion displays. Initially, we deposited single NCs as well as binary and ternary NC composites onto a substrate to demonstrate a proof of concept for full-color upconversion displays with point-by-point scanning the substrate under NIR-III/IV laser pumping (pulse width ~120 fs, repetition rate 5 kHz). To the best of our knowledge, this is the first reported instance of full-color MPA upconversion displays in NIR-III/IV windows (Fig. 5a, see detailed implementation methods in Supplementary Note 9). This innovation may hold promise for advancing full-color NIR-III/IV bioimaging in non-invasive in vivo microscopy applications by using these NCs as fluorophores. Although the presence of Pb element in the RPPs poses a challenge for bioimaging applications, this issue can be addressed by encapsulating the NC with silica or biocompatible core-shell structures as certified by a lot of successful works55,68,69,70,71,72,73. In particular, the solvent-free encapsulation strategies are demonstrated to be suitable for the soluble RPPs, such as in-situ photopolymerization of silicones68 and vapor-phase deposition of dense silica69 or fluorinated silane layers70.

Fig. 5: Paradigm of high-order MPA upconversion in full-color upconversion displays and multicolor information encryption in NIR-III/IV windows.

Fig. 5: Paradigm of high-order MPA upconversion in full-color upconversion displays and multicolor information encryption in NIR-III/IV windows.

a Implementation procedure of full-color upconversion displays prepared using single and composite NCs. b Concept of using NIR-III/IV fs laser to decode information. c Illustration of high-confidentiality multi-layer decryption by using multi-color MPA, wherein decoding information requires both NIR-III/IV pumping and spectrum filtering

The giant high-order multiphoton nonlinearities also offer an intriguing opportunity for applications in information decoding based on the upconversion display strategy. As shown in Fig. 5b, utilizing NIR-III/IV fs lasers scanning to excite 4PA-6PA effects, a mobile phone equipped with a standard silicon CCD can readily identify a distinct quick-response (QR) code fabricated by i-BMPB NCs. In contrast, excitation with common VIS continuous-wave (CW) lasers results in weak PL, and background noise making code retrieval impossible. Even NIR nanosecond (ns) laser fails to excite high-order MPA PL, resulting in an all-black picture. Besides, VIS and NIR-I/II fs lasers often cause CCD saturation because they can easily induce efficient 1PA and 2PA effects of CCD.

The potential of colorful MPA for multi-layer information encryption based on various MPA upconversion is further demonstrated in Fig. 5c. Compared to common method with NIR-I/II light, this scheme based on NIR-III/IV MPA allows for clearer decoding and decryption with higher spatial resolution, profited from the better scattering suppression of longer wavelengths. This increased resolution enables higher information storage density and superior potential for super-resolution image encryption74. Moreover, this specificity offers enhanced information security compared to conventional methods that use common CW lasers.